裂縫確定方法和裝置與流程
2023-09-17 19:09:00 1
本申請涉及油氣勘探技術領域，特別涉及一種裂縫確定方法和裝置。
背景技術：
在油氣勘探領域，地層中發育的裂縫因為不僅可以作為儲集空間，還可以對地層中分散、孤立的孔隙和溶洞起連通作用，增加其有效孔隙度，使之成為優質儲層，所以準確地預測地下裂縫空間展布是油氣勘探過程的重要研究內容。
為了確定地層中的裂縫，現有的裂縫確定方法一般可以分為相干法、曲率法、螞蟻追蹤法、體積膨脹係數分析法和疊前方位各向異性分析法等。其中，相干法是通過地震道組合分析空間不連續性來預測複雜巖性地層裂縫，可以定性刻畫裂縫空間展布趨勢。曲率法是通過獲取曲率屬性，對巖層的彎曲、褶皺和裂縫、斷層等進行分析，進而確定地層裂縫空間展布。螞蟻追蹤法是基於螞蟻算法刻畫地下裂縫空間展布，可以定性描述中小尺度裂縫的空間展布形態，但空間展布趨勢不清晰。體積膨脹係數分析法是通過對比地層原始沉積形態體積和現今形態體積，以預測地下裂縫發育情況。疊前方位各向異性分析法是通過分析不同方位角地震屬性的差異，來預測地下裂縫發育情況。
但是，具體實施時，現有裂縫確定方法中的相干法、螞蟻追蹤法往往只能進行定性預測，不能開展定量化描述。而曲率法、體積膨脹係數分析法和疊前方位各向異性分析法則可以定量化預測複雜巖性地層裂縫發育，但對裂縫空間展布的預測精度比較低。此外，以上裂縫確定方法在應用過程中大多未考慮地層中巖性的變化。但實際上在許多地層中不同部位的巖性會有較大差別，而不同巖性又對應不同脆性。其中，所述脆性一般是指當外力達到一定限度時，材料發生無先兆的突然破壞且破壞時無明顯塑性變形的性質。所以脆性的不同，意味著不同巖性產生裂縫的難易程度不同，形成裂縫的強度也不相同。而現有單項裂縫預測技術均未考慮地層中巖性變化情況，無法準確地定量化預測複雜巖性地層的裂縫空間展布信息。因此，如何利用地震數據，準確地定量化預測地下複雜巖性地層的裂縫空間展布信息成為地震儲層描述綜合研究中亟待解決的問題。綜上可知，現有的裂縫確定方法往往存在不能準確地定量化預測複雜巖性地層的裂縫的技術問題。
針對上述問題，目前尚未提出有效的解決方案。
技術實現要素：
本申請實施例提供了一種裂縫確定方法和裝置，以解決現有裂縫確定方法中存在的不能準確地定量化預測複雜巖性地層裂縫的技術問題。
本申請實施例提供了一種裂縫確定方法，包括：
獲取待測目的層的地震資料和測井資料；
根據所述地震資料和所述測井資料，確定所述目的層的縱波阻抗反演數據體；
根據所述測井資料，確定所述目的層中的多個縱波阻抗門檻值；
根據所述多個縱波阻抗門檻值、所述測井資料和所述縱波阻抗反演數據體，確定所述目的層的脆性數據體；
根據所述地震資料，確定目的層的曲率數據體；
根據所述目的層的脆性數據體和所述目的層的曲率數據體，確定所述目的層的裂縫屬性數據體；
根據所述裂縫屬性數據體，確定所述目的層的裂縫。
在一個實施方式中，根據所述地震資料和所述測井資料，確定所述目的層的縱波阻抗反演數據體，包括：
根據所述地震資料，確定所述目的層的構造解釋層位；
根據所述測井資料和所述構造解釋層位，建立地震反演模型；
根據所述地震反演模型，得到所述縱波阻抗反演數據體。
在一個實施方式中，在根據地震資料和測井資料，確定目的層的縱波阻抗反演數據體之前，所述方法還包括：
對所述地震資料進行保幅處理；
對所述地震資料和所述測井資料進行井震標定。
在一個實施方式中，根據所述測井資料，確定所述目的層中的多個縱波阻抗門檻值，包括：
根據所述測井資料，獲得測井解釋結果；
根據所述測井解釋結果，確定多個縱波阻抗門檻值。
在一個實施方式中，根據所述多個縱波阻抗門檻值和所述測井資料，確定所述目的層的脆性數據體，包括：
根據所述測井資料，分別確定所述多個巖性中各個巖性的脆性因子；
根據所述各個巖性的脆性因子，確定各個巖性的脆性相對比值；
根據所述各個巖性的脆性相對比值、所述多個縱波阻抗門檻值和所述縱波阻抗反演數據體，確定所述脆性數據體。
在一個實施方式中，在多個巖性包括：第一巖性和第二巖性的情況下，按照以下公式，確定所述脆性數據體：
其中，Ii為所述脆性數據體中標號為i的測點的相對脆性值，a為預設數值，根據所述第一巖性和所述第二巖性的脆性相對比值確定，AIi為所述縱波阻抗反演數據體中標號為i的測點的縱波阻抗反演數據，AIm1為所述第一巖性和所述第二巖性之間的縱波阻抗門檻值。
在一個實施方式中，根據所述地震資料，確定目的層的曲率數據體，包括：
根據所述地震資料，確定所述目的層的各個測點的傾角；
根據所述地震資料和所述目的層的各個測點的傾角，確定所述各個測點的曲率；
根據所述各個測點的曲率，確定所述目的層的曲率數據體。
在一個實施方式中，根據所述地震資料，確定所述目的層的各個測點的傾角，包括：
根據所述地震資料，獲取所述各個測點的瞬時頻率和瞬時波數；
根據所述各個測點的瞬時頻率和瞬時波數，確定所述各個測點的傾角。
在一個實施方式中，確定所述各個測點的曲率，包括：分別計算所述各個測點的多個曲率，將所述多個曲率中符合預設要求的曲率作為所述測點的曲率。
在一個實施方式中，根據所述目的層的脆性數據體和所述目的層的曲率數據體，確定目的層的裂縫屬性數據體，包括：
將所述目的層的脆性數據體和所述目的層的曲率數據體進行叉乘處理，得到所述目的層的裂縫屬性數據體。
基於相同的發明構思，本申請實施例還提供了一種裂縫確定裝置，包括：
資料獲取模塊，用於獲取待測目的層的地震資料和測井資料；
第一確定模塊，用於根據地震資料和測井資料，確定目的層的縱波阻抗反演數據體；
第二確定模塊，用於根據所述測井資料，確定所述目的層中多個縱波阻抗門檻值；
第三確定模塊，用於根據所述多個縱波阻抗門檻值、所述測井資料和所述縱波阻抗反演數據體，確定所述目的層的脆性數據體；
第四確定模塊，用於根據所述地震資料，確定目的層的曲率數據體；
第五確定模塊，用於根據所述目的層的脆性數據體和所述目的層的曲率數據體，確定目的層的裂縫屬性數據體；
裂縫確定模塊，用於根據所述裂縫屬性數據體，確定目的層的裂縫。
在一個實施方式中，所述第三確定模塊包括：
脆性因子確定單元，用於根據所述測井資料，分別確定所述多個巖性中各個巖性的脆性因子；
脆性相對比值確定單元，用於根據所述各個巖性的脆性因子，確定各個巖性的脆性相對比值；
脆性數據體確定單元，用於根據所述各個巖性的脆性相對比值、所述多個縱波阻抗門檻值和所述縱波阻抗反演數據體，確定所述脆性數據體。
在一個實施方式中，所述第四確定模塊包括：
傾角確定單元，用於根據所述地震資料，確定所述目的層的各個測點的傾角；
曲率確定單元，用於根據所述地震資料和所述目的層的各個測點的傾角，確定所述各個測點的曲率；
曲率數據體確定單元，用於根據所述各個測點的曲率，確定所述目的層的曲率數據體。
在本申請實施例中，通過考慮地層中不同巖性的變化特點，將不同巖性的脆性相對關係進行定量化描述，構造出地層的脆性數據體，再結合地層的曲率數據體來定量化地確定裂縫，從而解決了現有的裂縫確定方法存在的不能準確地定量化預測複雜巖性地層裂縫的技術問題，達到了提高裂縫預測精度的技術效果。
附圖說明
為了更清楚地說明本申請實施例或現有技術中的技術方案，下面將對實施例或現有技術描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本申請中記載的一些實施例，對於本領域普通技術人員來講，在不付出創造性勞動性的前提下，還可以根據這些附圖獲得其他的附圖。
圖1是根據本申請實施例的裂縫確定方法的處理流程圖；
圖2是根據本申請實施例的裂縫確定裝置的組成結構圖；
圖3是應用本申請實施例提供的裂縫確定方法/裝置的流程示意圖；
圖4是應用本申請實施例提供的裂縫確定方法/裝置得到的Z段地層裂縫預測結果示意圖；
圖5是應用本申請實施例提供的裂縫確定方法/裝置通過曲率算法得到的Z段地層裂縫預測結果示意圖；
圖6是應用本申請實施例提供的裂縫確定方法/裝置通過使用彈性參數得到的Z段地層巖石相對脆性預測結果。
具體實施方式
為了使本技術領域的人員更好地理解本申請中的技術方案，下面將結合本申請實施例中的附圖，對本申請實施例中的技術方案進行清楚、完整地描述，顯然，所描述的實施例僅僅是本申請一部分實施例，而不是全部的實施例。基於本申請中的實施例，本領域普通技術人員在沒有做出創造性勞動前提下所獲得的所有其他實施例，都應當屬於本申請保護的範圍。
考慮到現有的裂縫確定方法，因為往往沒有充分考慮到不同巖性的不同脆性對於裂縫發育的具體影響，導致存在具體實施時，存在不能準確地定量化預測複雜巖性地層裂縫的技術問題。針對產生上述技術問題的根本原因，又考慮到實際油氣勘探開發過程中，研究人員常常使用地層巖石中脆性礦物含量或彈性參數來構建脆性評價參數，表徵地層的相對脆性程度。其中：方法一，利用脆性礦物含量評價巖石脆性，但需要大量測試巖石中各類礦物的含量，導致實際操作成本很高。方法二，由於彈性參數表徵巖石脆性可以通過常規地震資料和測井資料獲得，可以根據材料的應力-應變關係，利用表徵徑向形變量的彈性模量和表徵橫向形變量的泊松比表徵脆性。例如，高彈性模量和低泊松比表示脆性程度高。該方法相對易於操作，實際應用廣泛。但是，單純地通過地層脆性評價參數雖然可以對不同巖性地層的脆性進行評價，但還是無法預測地層裂縫的空間展布信息。綜上，本申請考慮可以引入不同巖性變化對地層中裂縫的影響，即可以通過定量化地描述地層中不同巖性間的脆性相對關係，建立地層的脆性數據體，再將地層的曲率數據體和脆性數據體結合來確定裂縫。從而解決了現有的裂縫確定方法存在的不能準確地定量化預測複雜巖性地層裂縫的技術問題，達到了提高確定裂縫精度的技術效果。
基於上述思考思路，本申請提供了一種裂縫確定方法。請參閱圖1。本申請提供的裂縫確定方法，可以包括以下步驟。
步驟101：獲取待測目的層的地震資料和測井資料。
步驟102：根據所述地震資料和所述測井資料，確定所述目的層的縱波阻抗反演數據體。
在一個實施方式中，為了確定目的層縱波反抗反演數據體，需要根據採集得到地震資料和測井資料進行處理分析。具體可以包括以下步驟：
S1：根據所述地震資料，確定目的層的構造解釋層位；
S2：根據所述測井資料和所述目的層的構造解釋層位，建立地震反演模型；
S3：根據所述地震反演模型，得到所述目的層的縱波阻抗反演數據體。
在一個實施方式中，為了使用地震資料和測井資料獲得目的層的縱波阻抗反演數據體，還需要先對採集得到的地震資料和測井資料進行預處理，具體可以包括：對所述地震資料進行保幅處理；對所述地震資料和所述測井資料進行井震標定。其中，所述保幅處理具體可以是保幅處理是指在地震資料處理過程中保持振幅的相對關係。當然，還可以根據具體情況和施工要求，選擇其他相關的處理方式對地震資料和測井資料進行預處理。例如，刪除地震資料和測井資料中的無效數據或者對地震資料和測井資料進行校正等等。對此，本申請不作限定。
步驟103：根據所述測井資料，確定所述目的層中的多個縱波阻抗門檻值。
在一個實施方式中，考慮到目標地層中可能會包含有多個不同巖性的區域，而不同巖性對於地層中裂縫的發育有著重要的影響。為了將地層中不同巖性的影響引入對裂縫的確定工作中，首先可以根據測井資料，設置多個巖性間的縱波阻抗門檻值。具體實施時，可以包括：
S1：根據所述測井資料，獲得測井解釋結果；
S2：根據所述測井解釋結果，確定多個縱波阻抗門檻值。
在本實施方式中，需要說明的是當目標地層有兩個不同巖性時，如：第一巖性和第二巖性。對應地，可以根據測井解釋結果，這裡可以根據目的層段第一巖性地層和第二巖性地層的測井縱波阻抗值分布範圍，劃定第一巖性和第二巖性間的縱波阻抗門檻值。當目標地層有三個不同巖性時，如：第一巖性、第二巖性和第三巖性。其中，第一巖性地層的測井縱波阻抗值分布範圍相對與第二巖性地層的測井縱波阻抗值分布範圍接近，第二巖性地層的測井縱波阻抗值分布範圍相對與第三巖性地層的測井縱波阻抗值分布範圍接近。對應地，可以根據測井解釋結果，即根據目的層段第一巖性地層、第二巖性地層和第三巖性地層的測井縱波阻抗值分布範圍，劃定第一巖性和第二巖性間的縱波阻抗門檻值，以及第二巖性和第三巖性間的縱波阻抗門檻值。類似的，對於地層中有四個不同巖性或者更多不同巖性的情況可以參照上述實施方法。對此，本申請不再贅述。
步驟104：根據所述多個縱波阻抗門檻值、所述測井資料和所述縱波阻抗反演數據體，確定所述目的層的脆性數據體。
在一個實施方式中，考慮到不同巖性的差異對裂縫發育的影響，具體實施時可以按照以下方法確定目的層中的不同巖性的脆性數據體：
S1：據所述測井資料，分別確定所述多個巖性中各個巖性脆性因子。
在本實施方式中，需要說明的是所述脆性是指當外力達到一定限度時，材料發生無先兆的突然破壞且破壞時無明顯塑性變形的性質。所以脆性不同意味著不同巖性產生裂縫的難易程度不同，形成裂縫的強度也不相同。
在一個實施方式中，考慮到利用脆性礦物含量評價巖石脆性需要大量測試巖石中各類礦物的含量，實際操作成本很高。本申請考慮可以根據地震資料和測井資料確定表徵徑向形變量的彈性模量和表徵橫向形變量的泊松比來評價脆性，即脆性因子。具體實施時，可以按照以下的方法處理。
據所述測井資料，按照以下公式，計算地層不同巖性的巖石脆性評價參數，即脆性因子IRBI：
μ＝ρVs2
IRBI＝E/υ
其中：E為楊氏模量，υ為泊松比，μ為剪切模量，λ為彈性模量，需要說明的是上述這4個參數為地層巖石的彈性參數；ρ、Vp、Vs分別為對應巖性地層的密度、縱波速度和橫波速度，相應數值，具體可以從對應目的層的測井資料中獲取。
S2：根據所述各個巖性的脆性因子，分別確定各個巖性的脆性相對比值。
在一個實施方式中，為了考慮不同巖性的影響，可以通過確定不同巖性間的脆性相對比值來將不同巖性的影響引入到具體的裂縫確定的過程中。具體地，可以將不同巖性的脆性因子相除，以得到的商值作為脆性相對比值。例如，目的層包括兩種巖性，第一巖性和第二巖性。其中第一巖性的反演數據體的脆性因子為1，第二巖性的反演數據體的脆性因子為2，則第一巖性和第二巖性的脆性相對比值就是用第一巖性的反演數據體的脆性因子1除以第二巖性的反演數據體的脆性因子2，得到的商值0.5。對於，地層包括其他數量的巖性情況可以參照兩種巖性情況處理。對此，本申請不再贅述。
S3：根據所述各個巖性的脆性相對比值、所述多個縱波阻抗門檻值和所述縱波阻抗反演數據體，確定所述脆性數據體。
在一個實施方式中，在多個巖性為兩個巖性，即：第一巖性和第二巖性的情況下，具體實施可以按照以下公式，確定所述脆性數據體：
其中，Ii為所述脆性數據體中標號為i的測點的相對脆性值，a為預設數值，根據所述第一巖性和所述第二巖性的脆性相對比值確定，AIi為所述縱波阻抗反演數據體中標號為i的測點的縱波阻抗反演數據，AIm1為所述第一巖性和所述第二巖性之間的縱波阻抗門檻值。
需要說明的是，多個巖性還可以是三個、四個甚至更多個不同巖性的情況，具體實施時都可以參照上述兩個巖性的情況進行處理。對此，本申請不作限定。
步驟105：根據所述地震資料，確定目的層的曲率數據體。
在一個實施方式中，考慮可以採用曲率法作為依據設計步驟先求取目的層的曲率數據體。具體實施，可以按照以下的步驟執。
S1：根據所述地震資料，確定所述目的層的各個測點的傾角。
在一個實施方式中，為了準確確定各個測點的傾角，具體執行可以包括：
S1-1：根據所述地震資料，獲取所述各個測點的瞬時頻率和瞬時波數；
S1-2：根據所述各個測點的瞬時頻率和瞬時波數，確定所述各個測點的傾角。
S2：根據所述地震資料和所述目的層的各個測點的傾角，確定所述各個測點的曲率。
在一個實施方式中，所述各個測點的曲率，包括：分別計算所述各個測點的多個曲率，將所述多個曲率中符合預設要求的曲率作為所述測點的曲率。例如，可以根據具體情況，從所述多個曲率中選擇各個測點的最大正曲率Kpos或最大負曲率為Kneg作為該測點的曲率。
S3：根據所述各個測點的曲率，確定所述目的層的曲率數據體。
在一個實施例中，已知地震資料中的在空間任意點的數據是一個時間標量，即為：u(t,x,y)，則瞬時頻率ω為：
其中：uH(t,x,y)為關於時間t的希爾伯特變換。
沿x方向和y方向的瞬時波數kx、ky分別是：
則傾角dip(px,qy)可由瞬時波數kx、ky和瞬時頻率ω計算得到：
其中：px、qy分別為x方向和y方向的傾角分量。
三維空間中某點的曲率值通過與其相鄰道和樣點的傾角值擬合出的空間曲面方程進行計算，根據最小二乘逼近原理為：
z(x,y)＝ax2+by2+cxy+dx+ey+f
對方程兩邊求微分，帶入傾角分量px、qy，從而得出方程係數為：
則對描述裂縫構造較為有效的最大正曲率Kpos和最大負曲率為Kneg可表示為：
將相對脆性數據體I與曲率屬性體K進行融合得到新的定量化裂縫預測屬性體Z，即：
Zi＝Ii×Ki i＝1,2,3…,N
其中，Zi為新的定量化裂縫預測屬性體的任意採樣點，Ii為相對脆性數據體對應位置的採樣點，Ki為曲率屬性體對應位置的採樣點,N為採樣點的總個數。
步驟106：根據所述目的層的脆性數據體和所述目的層的曲率數據體，確定目的層的裂縫屬性數據體。
在一個實施方式中，為了將不同巖性對裂縫發育的影響作用加入確定裂縫方法中，可以綜合利用得到目的層的脆性數據體和目的層的曲率數據體。具體實施時，可以包括將所述目的層的脆性數據體和所述目的層的曲率數據體進行叉乘處理，得到所述目的層的裂縫屬性數據體。從而，可以實現基於不同巖性的裂縫分析預測。
步驟107：根據所述裂縫屬性數據體，確定目的層的裂縫。
在一個實施方式中，為了確定目的層的裂縫。具體實施時，可以根據裂縫屬性數據體中各個測點的裂縫屬性數據，生成目的層各個測點的裂縫分布圖。其中，該分布圖可以清楚、準確地展示目的層內的裂縫展布信息。根據裂縫展布信息，可以進行後續的油氣勘探。需要說明的是除了上述列舉的根據裂縫屬性數據體生成裂縫分布圖外還可以通過其他的方式利用裂縫屬性數據體，確定目的層的裂縫。對此，本申請不作限定。
在本申請實施例中，通過考慮不同巖性變化對裂縫發育的影響，定量地描述不同巖性間的脆性相對關係，再將地層中的脆性相對關係與曲率法結合，從而解決了現有的裂縫確定方法存在的不能準確地定量化預測複雜巖性地層裂縫的技術問題，達到了提高確定裂縫準確度的技術效果。
基於同一發明構思，本發明實施例中還提供了一種裂縫確定裝置，如下面的實施例所述。由於裝置解決問題的原理與裂縫確定方法相似，因此裂縫確定裝置的實施可以參見裂縫確定方法的實施，重複之處不再贅述。以下所使用的，術語「單元」或者「模塊」可以實現預定功能的軟體和/或硬體的組合。儘管以下實施例所描述的裝置較佳地以軟體來實現，但是硬體，或者軟體和硬體的組合的實現也是可能並被構想的。請參閱圖2，是本發明實施例的裂縫確定裝置的一種組成結構圖，該裝置可以包括：資料獲取模塊201、第一確定模202、第二確定模塊203、第三確定模塊204、第四確定模塊205、第五確定模塊206和裂縫確定模塊207，下面對該結構進行具體說明。
資料獲取模塊201，用於獲取地震資料和測井資料。
第一確定模塊202，用於根據所述地震資料和所述測井資料，確定目的層的縱波阻抗反演數據體。
第二確定模塊203，用於根據所述測井資料，確定所述目的層中多個縱波阻抗門檻值。
第三確定模塊204，用於根據所述多個縱波阻抗門檻值、所述測井資料和所述目的層的縱波阻抗反演數據體，確定所述目的層的脆性數據體。
第四確定模塊205，用於根據所述地震資料，確定目的層的曲率數據體。
第五確定模塊206，用於根據所述目的層的脆性數據體和所述目的層的曲率數據體，確定目的層的裂縫屬性數據體。
裂縫確定模塊207，用於根據所述裂縫屬性數據體，確定目的層的裂縫。
在一個實施方式中，為了確定目的層的縱波阻抗反演數據體，第一確定模塊202具體可以包括：
解釋層位確定單元，用於根據所述地震資料，確定目的層的構造解釋層位；
地震反演模型建立單元，用於根據所述測井資料和所述目的層的構造解釋層位，建立地震反演模型；
反演數據體獲取單元，用於根據所述地震反演模型，得到所述目的層的縱波阻抗反演數據體。
在一個實施方式中，為了在使用前，對地震資料和測井資料進行預處理，所述第一確定模塊202還包括可以包括預處理單元，用於執行以下預處理：對所述地震資料進行保幅處理；和對所述地震資料和所述測井資料進行井震標定。
在一個實施方式中，為了確定多個縱波阻抗門檻值，第二確定模塊203具體可以包括：
測井解釋結果獲取單元，用於根據所述測井資料，獲得測井解釋結果；
門檻值確定單元，用於根據所述測井解釋結果，確定多個縱波阻抗門檻值。
在一個實施方式中，為了確定目的層的脆性數據體，所述第三確定模塊204具體可以包括：
脆性因子確定單元，用於根據所述測井資料，分別確定所述多個巖性中各個巖性的脆性因子；
脆性相對比值確定單元，用於根據所述各個巖性的脆性因子，確定各個巖性的脆性相對比值；
脆性數據體確定單元，用於根據所述各個巖性的脆性相對比值、所述多個縱波阻抗門檻值和所述縱波阻抗反演數據體，確定所述脆性數據體。
在一個實施方式中，在多個巖性為兩個巖性：第一巖性和第二巖性的情況下，脆性數據體確定單元具體可以按照以下公式，確定所述脆性數據體：
其中，Ii為所述脆性數據體中標號為i的測點的相對脆性值，a為預設數值，根據所述第一巖性和所述第二巖性的脆性相對比值確定，AIi為所述縱波阻抗反演數據體中標號為i的測點的縱波阻抗反演數據，AIm1為所述第一巖性和所述第二巖性之間的縱波阻抗門檻值。需要說明的是，參數a的取值範圍一般可以為0.01-100，具體可以根據脆性相對比值確定。例如，第一巖性和第二巖性的脆性相對比值為0.5，可以根據具體情況將參數a的預設數值設為0.5，這時，只要滿足AIi＜AIm1條件的脆性數據體中任意點測點的相對脆性值可以取為0.5。從而，可以引入不同巖性對裂縫發育產生的影響。
在一個實施方式中，為了確定目的層的曲率數據體，所述第四確定模塊205具體可以包括：
傾角確定單元，用於根據所述地震資料，確定所述目的層的各個測點的傾角；
曲率確定單元，用於根據所述地震資料和所述目的層的各個測點的傾角，確定所述各個測點的曲率；
曲率數據體確定單元，用於根據所述各個測點的曲率，確定所述目的層的曲率數據體。
在一個實施方式中，為了確定所述目的層的各個測點的傾角，所述傾角確定單元具體可以包括：
參數獲取子單元，用於根據所述地震資料，獲取所述各個測點的瞬時頻率和瞬時波數；
傾角確定子單元，用於根據所述各個測點的瞬時頻率和瞬時波數，確定所述各個測點的傾角。
在一個實施方式中，所述曲率確定單元可以分別計算所述各個測點的多個曲率，將所述多個曲率中符合預設要求的曲率作為所述測點的曲率。例如，可以根據具體情況，從所述多個曲率中選擇測點的最大正曲率Kpos或最大負曲率為Kneg作為該測點的曲率。
在一個實施方式中，為了將得到的目的層的曲率數據體和目的層的脆性數據體進行融合，所述第五確定模塊206具體可以將所述目的層的脆性數據體和所述目的層的曲率數據體進行叉乘處理，得到所述目的層的裂縫屬性數據體。
在一個實施方式中，為了確定目的層的裂縫，裂縫確定模塊207具體可以根據裂縫屬性數據體中各個測點的裂縫屬性數據，生成目的層各個測點的裂縫分布圖。其中，該分布圖可以清楚、準確地展示目的層內的裂縫展布信息。根據裂縫展布信息，可以進行後續的油氣勘探。需要說明的是除了上述列舉的根據裂縫屬性數據體生成裂縫分布圖外還可以通過其他的方式利用裂縫屬性數據體，確定目的層的裂縫。對此，本申請不作限定。
本說明書中的各個實施例均採用遞進的方式描述，各個實施例之間相同相似的部分互相參見即可，每個實施例重點說明的都是與其他實施例的不同之處。尤其，對於系統實施例而言，由於其基本相似於方法實施例，所以描述的比較簡單，相關之處參見方法實施例的部分說明即可。
需要說明的是，上述實施方式闡明的系統、裝置、模塊或單元，具體可以由計算機晶片或實體實現，或者由具有某種功能的產品來實現。為了描述的方便，在本說明書中，描述以上裝置時以功能分為各種單元分別描述。當然，在實施本申請時可以把各單元的功能在同一個或多個軟體和/或硬體中實現。
此外，在本說明書中，諸如第一和第二這樣的形容詞僅可以用於將一個元素或動作與另一元素或動作進行區分，而不必要求或暗示任何實際的這種關係或順序。在環境允許的情況下，參照元素或部件或步驟(等)不應解釋為局限於僅元素、部件、或步驟中的一個，而可以是元素、部件、或步驟中的一個或多個等。
從以上的描述中，可以看出，本申請實施例提供的裂縫確定方法和裝置。通過定量確定地層中不同巖性間的脆性相對關係將不同巖性變化對地層中裂縫發育的影響融入了通過曲率確定裂縫的方法，即根據地層的脆性數據體和地層的曲率數據體確定裂縫，從而解決了現有裂縫確定方法存在的不能準確地定量化預測複雜巖性地層裂縫的技術問題，達到了在定量化確定地層裂縫的同時，提高確定裂縫準確度的技術效果。
在一個具體實施場景，應用本申請提供裂縫確定方法/裝置對四川盆地Z段地層進行裂縫預測。具體實施流程，可以參閱圖3。主要包括以下七個步驟。
(1)收集保幅處理後的偏移疊加地震數據，地震工區內鑽穿目的層的測井數據，井-震標定後開展地震數據目的層構造解釋，建立基於解釋層位和標定後測井數據的地震反演模型，提取地震數據體子波，帶入地震反演模型進行反演獲得目的層段縱波阻抗反演數據體AI。
(2)根據測井解釋結果設定目的層段地層不同巖性間的縱波阻抗門檻值，目的層段地層包含灰巖和泥質灰巖兩種巖性，灰巖縱波阻抗值大等於13000，泥質灰巖縱波阻抗值小於13000，縱波阻抗單位為g/cm3·m/s，則目的層段地層巖性間縱波阻抗門檻值AImax為13000。
(3)計算地層不同巖性的脆性相對關係；首先計算地層不同巖性的巖石脆性評價參數-脆性因子IRBI，然後通過比較地層不同巖性的脆性因子計算其脆性相對關係，這裡脆性因子IRBI主要通過彈性參數表徵：
IRBI＝E/υ
其中：
μ＝ρVs2
式中：E為楊氏模量，υ為泊松比，μ為剪切模量，λ為彈性模量，這4個參數為地層巖石的彈性參數，ρ、Vp、Vs分別為對應巖性地層的密度、縱波速度和橫波速度；本工區目的層段灰巖的密度、縱波速度和橫波速度分別為2.7g/cm3、6300m/s和3300m/s，而泥質灰巖的密度、縱波速度和橫波速度分別為2.8g/cm3、5350m/s和2750m/s，則灰巖和泥質灰巖二者的脆性相對關係為1：0.7；
(4)根據巖性差異將目的層段縱波阻抗反演數據體AI轉化為相對脆性數據體I，這裡設定目的層段地層主要巖性灰巖脆性值為1，則泥質灰巖根據其與地層主要巖性灰巖的脆性相對關係設定為0.7；則
式中：Ii為相對脆性數據體I的任意採樣點，0.7為泥質灰巖與地層主要巖性灰巖的脆性相對數值，13000為目的層段地層灰巖和泥質灰巖兩種巖性間的縱波阻抗門檻值。
(5)求取地震數據體的曲率屬性。三維地震數據體先轉化為傾角數據體，再根據傾角數據體計算數據體內任意反射點的曲率。三維地震數據體的傾角可以由瞬時頻率和瞬時波數計算。已知三維地震數據體在空間任意點是時間標量u(t,x,y)，則瞬時頻率ω為：
式中：uH(t,x,y)為關於時間t的希爾伯特變換。
沿x方向和y方向的瞬時波數kx、ky分別是：
則傾角dip(px,qy)可由瞬時波數kx、ky和瞬時頻率ω計算得到：
式中：px、qy分別為x方向和y方向的傾角分量。
三維空間中某點的曲率值通過與其相鄰道和樣點的傾角值擬合出的空間曲面方程進行計算，根據最小二乘逼近原理為：
z(x,y)＝ax2+by2+cxy+dx+ey+f
對方程兩邊求微分，帶入傾角分量px、qy，從而得出方程係數為：
則對描述裂縫構造較為有效的最大正曲率Kpos可表示為：
Kpos＝(a+b)+[(a-b)2+c2]1/2
(6)將相對脆性數據體I與曲率屬性體K進行融合得到新的定量化裂縫預測屬性體Z
Zi＝Ii×Ki i＝1,2,3…,N
式中：Zi為新的定量化裂縫預測屬性體的任意採樣點，Ii為相對脆性數據體對應位置的採樣點，Ki為曲率屬性體對應位置的採樣點。
(7)利用最終得到的定量化裂縫預測屬性體預測地下裂縫的空間展布信息。
在本實施例中，最終Z段地層基於針對複雜巖性地層裂縫定量化預測方法得到的結果如圖4所示。圖5為利用原始地震資料經過曲率算法得到的裂縫預測結果。圖6為利用原始地震、測井資料得到的巖石相對脆性預測結果。本地區Z地層巖性主要為灰巖，在中東部和北部局部地區分布有泥質灰巖，整體上脆性大易於產生裂縫。對比圖4和圖5可知，本發明實施例提供的一種基於針對複雜巖性地層的裂縫定量化預測方法進步在於在地層裂縫預測過程中考慮到地層巖性變化情況，根據不同巖性的脆性評價參數設定相對脆性比值，並將其應用於地層的裂縫定量化描述過程中，使其更加符合實際情況，最終得到準確的裂縫預測結果。
綜上可知，本申請實施例提供的裂縫確定方法/裝置使用時的有益成果是：提供了一種針對複雜巖性地層的裂縫定量化預測方法，利用地震縱波阻抗反演數據門檻值區分不同巖性，將目的層段地層不同巖性空間展布範圍進行精細刻畫，然後用彈性參數將目的層段地層不同巖性的脆性相對關係進行定量化描述，構建地層相對脆性數據體，將其與曲率屬性體融合，生成的新屬性體可以準確地定量化描述地層裂縫空間展布信息，為儲層精細評價和油氣高效開發提供技術支撐。驗證了通過本申請提供的確定裂縫方法和裝置確實可以解決現有裂縫確定方法具體使用時存在的不能準確地定量化預測複雜巖性地層的裂縫的技術問題。同時，也驗證了本申請具體實施時具有如下幾個優點：
1、在裂縫預測過程中考慮到地層不同巖性的脆性特徵對裂縫發育強度的影響，使其最終的定量化預測結果更加符合地層裂縫發育的真實情況。
2、在地震縱波阻抗反演結果中設定門檻值，有效刻畫地層不同巖性的空間展布範圍，再通過彈性參數計算地層不同巖性的脆性相對關係，從而為準確地定量化描述複雜巖性地層的裂縫空間展布信息奠定基礎。
儘管本申請內容中提到不同的裂縫確定方法或裝置，但是，本申請並不局限於必須是行業標準或實施例所描述的情況等，某些行業標準或者使用自定義方式或實施例描述的實施基礎上略加修改後的實施方案也可以實現上述實施例相同、等同或相近、或變形後可預料的實施效果。應用這些修改或變形後的數據獲取、處理、輸出、判斷方式等的實施例，仍然可以屬於本申請的可選實施方案範圍之內。
雖然本申請提供了如實施例或流程圖所述的方法操作步驟，但基於常規或者無創造性的手段可以包括更多或者更少的操作步驟。實施例中列舉的步驟順序僅僅為眾多步驟執行順序中的一種方式，不代表唯一的執行順序。在實際中的裝置或客戶端產品執行時，可以按照實施例或者附圖所示的方法順序執行或者並行執行(例如並行處理器或者多線程處理的環境，甚至為分布式數據處理環境)。術語「包括」、「包含」或者其任何其他變體意在涵蓋非排他性的包含，從而使得包括一系列要素的過程、方法、產品或者設備不僅包括那些要素，而且還包括沒有明確列出的其他要素，或者是還包括為這種過程、方法、產品或者設備所固有的要素。在沒有更多限制的情況下，並不排除在包括所述要素的過程、方法、產品或者設備中還存在另外的相同或等同要素。
上述實施例闡明的裝置或模塊等，具體可以由計算機晶片或實體實現，或者由具有某種功能的產品來實現。為了描述的方便，描述以上裝置時以功能分為各種模塊分別描述。當然，在實施本申請時可以把各模塊的功能在同一個或多個軟體和/或硬體中實現，也可以將實現同一功能的模塊由多個子模塊的組合實現等。以上所描述的裝置實施例僅僅是示意性的，例如，所述模塊的劃分，僅僅為一種邏輯功能劃分，實際實現時可以有另外的劃分方式，例如多個模塊或組件可以結合或者可以集成到另一個系統，或一些特徵可以忽略，或不執行。
本領域技術人員也知道，除了以純計算機可讀程序代碼方式實現控制器以外，完全可以通過將方法步驟進行邏輯編程來使得控制器以邏輯門、開關、專用集成電路、可編程邏輯控制器和嵌入微控制器等的形式來實現相同功能。因此這種控制器可以被認為是一種硬體部件，而對其內部包括的用於實現各種功能的裝置也可以視為硬體部件內的結構。或者甚至，可以將用於實現各種功能的裝置視為既可以是實現方法的軟體模塊又可以是硬體部件內的結構。
本申請可以在由計算機執行的計算機可執行指令的一般上下文中描述，例如程序模塊。一般地，程序模塊包括執行特定任務或實現特定抽象數據類型的例程、程序、對象、組件、數據結構、類等等。也可以在分布式計算環境中實踐本申請，在這些分布式計算環境中，由通過通信網絡而被連接的遠程處理設備來執行任務。在分布式計算環境中，程序模塊可以位於包括存儲設備在內的本地和遠程計算機存儲介質中。
通過以上的實施方式的描述可知，本領域的技術人員可以清楚地了解到本申請可藉助軟體加必需的通用硬體平臺的方式來實現。基於這樣的理解，本申請的技術方案本質上或者說對現有技術做出貢獻的部分可以以軟體產品的形式體現出來，該計算機軟體產品可以存儲在存儲介質中，如ROM/RAM、磁碟、光碟等，包括若干指令用以使得一臺計算機設備(可以是個人計算機，移動終端，伺服器，或者網絡設備等)執行本申請各個實施例或者實施例的某些部分所述的方法。
本說明書中的各個實施例採用遞進的方式描述，各個實施例之間相同或相似的部分互相參見即可，每個實施例重點說明的都是與其他實施例的不同之處。本申請可用於眾多通用或專用的計算機系統環境或配置中。例如：個人計算機、伺服器計算機、手持設備或可攜式設備、平板型設備、多處理器系統、基於微處理器的系統、置頂盒、可編程的電子設備、網絡PC、小型計算機、大型計算機、包括以上任何系統或設備的分布式計算環境等等。
雖然通過實施例描繪了本申請，本領域普通技術人員知道，本申請有許多變形和變化而不脫離本申請的精神，希望所附的權利要求包括這些變形和變化而不脫離本申請。