一種基于數(shù)據(jù)重構(gòu)預測巖性體的方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及地球物理勘探中的地震資料解釋處理技術(shù)領(lǐng)域，具體的說是一種利用 數(shù)據(jù)體重構(gòu)來預測不同巖性體的方法。
【背景技術(shù)】
[0002] 目前，在油氣勘探中巖性預測是一個比較重要的環(huán)節(jié)，不同的巖性體對油氣的儲 集能力基本上不一樣，如高孔隙度的砂巖及溶孔發(fā)育的礁灰?guī)r對油氣具有良好的儲集能 力，所以在石油勘探中是相當重要的勘探目標。如在四川盆地探區(qū)中的海相及湖相地層中 尋找這類目標對該區(qū)的油氣勘探具有重要意義，鉆遇此類巖性的鉆井大多數(shù)經(jīng)測試均獲工 業(yè)氣流。
[0003] 常規(guī)的巖性體預測技術(shù)方法有很多種，基本上可分為疊前、疊后反演兩大類，經(jīng)過 反演或?qū)傩杂嬎憧傻玫结槍r性體預測的屬性數(shù)據(jù)體。但對于一個研究區(qū)內(nèi)有多種類型的 巖性體時，往往此時的單一屬性體不能很好地對目的層段中的多種巖性體進行描述，預測 結(jié)果往往具有很強的片面性;也有的方法是通過井上的巖性體、疊后波阻抗反演數(shù)據(jù)體進 行巖性體解釋，得到相關(guān)巖性體空間上的分布情況。一些專利如發(fā)明《基于數(shù)學形態(tài)學的巖 性體自動追蹤方法》（專利【申請?zhí)枴?01410154418.8)的專利公開提供一種基于數(shù)學形態(tài)學 的巖性體自動追蹤方法，該方法包括選取用于巖性體自動追蹤的地震屬性體，并進行數(shù)據(jù) 篩選;數(shù)據(jù)篩選結(jié)束后，檢測保留下來的數(shù)據(jù)點的空間位置信息，對這些數(shù)據(jù)點進行砂體分 類和編號;利用開運算對處理后得到的砂體邊界做處理，斷開兩個粘連砂體間的窄通道，實 現(xiàn)砂體的分離；以及根據(jù)最小砂體大小的設(shè)置，對應(yīng)用開運算處理之后得到的砂體進行重 新分類編號?？梢娊陙?，對巖性體預測的技術(shù)方法很多，研究相對比較成熟。但是，針對巖 性體的預測技術(shù)還不夠完善，主要體現(xiàn)在以下幾個方面：
[0004] (1)采用單一數(shù)據(jù)體進行不同巖性體的預測結(jié)果往往針對性不強，準確度不高。 [0005] (2)利用數(shù)據(jù)交會分析實施不同巖性體預測往往帶有人為主觀因素，并且預測的 準確率不高。
[0006] 為了針對上述傳統(tǒng)巖性體描述技術(shù)方法的不足，解決對同一工區(qū)內(nèi)的多種巖性體 進行相應(yīng)的區(qū)分及描述，本發(fā)明提供了一個相對簡單的技術(shù)流程來實現(xiàn)對工區(qū)內(nèi)的多種巖 性體進行預測，更方便、準確地了解不同巖性體的分布情況，從而更好的為油氣勘探服務(wù)。

【發(fā)明內(nèi)容】

[0007] 本發(fā)明目的在于克服現(xiàn)有常規(guī)巖性體預測技術(shù)中的上述和其他缺點，可以準確地 預測出不同巖性體的平面及垂向上的分布狀態(tài)。為此，本發(fā)明提供一種利用數(shù)據(jù)體重構(gòu)來 預測不同巖性體分布情況的方法。
[0008] 本發(fā)明的技術(shù)方案是：
[0009] -種基于數(shù)據(jù)重構(gòu)預測巖性體的方法，它包括以下步驟：
[0010] S1、對待測區(qū)域的地震資料進行提取和反演，獲取該區(qū)域多種巖性體的預測數(shù)據(jù) 體，并獲得各巖性體所對應(yīng)的預測數(shù)據(jù)體的門檻值；
[0011] S2、對步驟1所確定門檻值后的各個預測數(shù)據(jù)體進行數(shù)據(jù)重構(gòu)后加權(quán)處理，得到一 個用于巖性預測的加權(quán)重構(gòu)數(shù)據(jù)體，根據(jù)前述加權(quán)重構(gòu)數(shù)據(jù)體獲取待測區(qū)域的巖性。
[0012] 本發(fā)明的步驟S1中，預測數(shù)據(jù)體為利用疊前或疊后地震資料反演或?qū)傩蕴崛『笏?得到的數(shù)據(jù)體，或者為兩個或以上的數(shù)據(jù)體之間的數(shù)學運算所得到的數(shù)據(jù)體，在用于巖性 體預測的數(shù)據(jù)體中，如果一個數(shù)據(jù)體能夠?qū)崿F(xiàn)巖性體預測時，則使用該數(shù)據(jù)體進行數(shù)據(jù)重 構(gòu)、加權(quán)計算；否則應(yīng)該使用兩個或兩個以上的數(shù)據(jù)體分別實施數(shù)據(jù)重構(gòu)、加權(quán)處理后，得 到一個加權(quán)重構(gòu)數(shù)據(jù)體用于巖性體預測。
[0013] 本發(fā)明的巖性體包括砂巖、泥巖、灰?guī)r、頁巖以及生物礁，前述巖性體的物理特征 各不相同，具有差異性。
[0014] 本發(fā)明的步驟S1中，不同巖性體的預測數(shù)據(jù)體門檻值設(shè)定包括以下步驟：
[0015] Sl-a、根據(jù)待測區(qū)域中各井的測井資料（由測井設(shè)備記錄）、巖性資料(錄井地質(zhì)資 料)獲取不同巖性體的深度范圍數(shù)據(jù)，并根據(jù)測井資料從預測數(shù)據(jù)體中提取相關(guān)巖性的預 測數(shù)據(jù)體；（根據(jù)各井的測井資料及巖性資料統(tǒng)計及確定不同巖性體的深度范圍數(shù)據(jù)，并根 據(jù)不同的巖性體的物理差異特征采取合適的地震資料、反演或?qū)傩杂嬎惴椒ㄌ崛∠嚓P(guān)的數(shù) 據(jù)體，要求相關(guān)預測數(shù)據(jù)體在數(shù)據(jù)值方面能將預測巖性體與其它非預測巖性體區(qū)分開來）
[0016] si-b、利用測井資料及巖性資料對各相關(guān)巖性的預測數(shù)據(jù)體進行井一震標定，確 定各井的時深關(guān)系；
[0017] S1-C、根據(jù)井一震標定所得到的時深關(guān)系將步驟si-a中獲取的不同巖性體的深度 范圍數(shù)據(jù)轉(zhuǎn)換到相關(guān)數(shù)據(jù)體的時間域中，得到各巖性體所對應(yīng)的預測數(shù)據(jù)體的數(shù)據(jù)值，并 綜合各井相關(guān)的同一巖性體的數(shù)據(jù)構(gòu)成數(shù)據(jù)集，從各巖性體的數(shù)據(jù)集中確定該巖性體預測 數(shù)據(jù)的門檻值范圍，依此類推，從而實現(xiàn)不同巖性體的門檻值設(shè)置。
[0018] 本發(fā)明的步驟Sl-c中：對于各巖性體，從數(shù)據(jù)集中確定出最大及最小門檻值，其中 門檻值求取方式為綜合各井相關(guān)的同一巖性體的數(shù)據(jù)構(gòu)成數(shù)據(jù)集{Jd，門檻值求取的計算 公式如下：
[0019] Jmax=max{Ji} (1)
[0020] Jmin=min{Ji} (2)
[0021] 式(1) 一 (2)中心為井中該巖性段所對應(yīng)數(shù)據(jù)體的第i個采樣點的數(shù)據(jù)值，Jmax為該 巖性體的最大門檻值，J min為該巖性體的最小門檻值。
[0022]本發(fā)明的步驟S2具體為：
[0023] S2_a、采用步驟Sl-c獲取的各巖性體的門檻值對步驟Sl-a獲取的相關(guān)巖性的預測 數(shù)據(jù)體進行重構(gòu)賦值得到重構(gòu)數(shù)據(jù)體:對于不同的巖性體將各自門檻值內(nèi)的數(shù)據(jù)值賦值為 某一常數(shù)值，不同的巖性體對應(yīng)不同的常數(shù)值，各個門檻值外的數(shù)據(jù)均賦值為零值，以示對 巖性體的數(shù)據(jù)進行區(qū)分。依此類推，得到一系列各個巖性體對應(yīng)的重構(gòu)數(shù)據(jù)體。
[0024] S2_b、對各個巖性體的重構(gòu)數(shù)據(jù)體進行加權(quán)處理，得到一個加權(quán)重構(gòu)數(shù)據(jù)體。（首 先，針對不同的巖性體預測使用分配好的數(shù)據(jù)體進行，并設(shè)置相關(guān)巖性體的門檻值。要求一 個數(shù)據(jù)體可以用于預測一種或以上的巖性體，也可將不同的巖性體分開使用不同的數(shù)據(jù)體 進行相關(guān)的預測計算，要求不同的數(shù)據(jù)體不能預測同一巖性體;然后，不同的數(shù)據(jù)體根據(jù)分 配好預測的巖性體的門檻值進行重構(gòu)賦值，門檻值外的數(shù)據(jù)設(shè)置為零值，門檻值內(nèi)的數(shù)據(jù) 設(shè)置為某一個常數(shù)值，不同巖性體的常數(shù)值不能相同，并且兩種或以上的代表巖性的常數(shù) 值進行加權(quán)計算時所得的結(jié)果不能與任一巖性體的常數(shù)值相等;最后，對各個數(shù)據(jù)體進行 重構(gòu)賦值后針對同一道同一個采樣點進行加權(quán)計算得到一個加權(quán)重構(gòu)數(shù)據(jù)體，利用該數(shù)據(jù) 體進行巖性體預測。若當加權(quán)后所得到的數(shù)據(jù)值不與任何一個巖性體的常數(shù)值一致時，應(yīng) 根據(jù)實際情況、地質(zhì)資料及井上巖性資料進行再次重新賦值，賦值為某一個預測巖性體的 常數(shù)值。一般情況下，應(yīng)根據(jù)該采樣點的上、下、左及右采樣點的常數(shù)值進行賦值為某一個 巖性體的常數(shù)值）
[0025] 本發(fā)明的步驟S2_a具體為:根據(jù)門檻值范圍實施數(shù)據(jù)體的數(shù)據(jù)重構(gòu)：
[0026] 若某一巖性的預測數(shù)據(jù)體中采樣點的數(shù)據(jù)值在該巖性體的門檻值范圍內(nèi)，則將該 數(shù)據(jù)值設(shè)置為某一個常數(shù)K值;若采樣點的數(shù)據(jù)值在門檻值范圍外，則將該采樣點的數(shù)據(jù)值 賦值為零即：
[0027] a、若某一個預測數(shù)據(jù)體的第i個采樣點的數(shù)據(jù)值為土，假設(shè)當土<1_時，J_為該 數(shù)據(jù)體針對某一預測巖性體的最小門檻值，則& = 〇;
[0028] b、假設(shè)當Jmi"JdJmax時，Jmax為該數(shù)據(jù)體針對該預測巖性體的最大門檻值，則心 =K，要求K值一般情況下為一常數(shù)值，不同巖性體的常數(shù)K值不能相同，并且K值之間的組合 后加權(quán)結(jié)果不能與任何一個巖性體設(shè)定的Κ值一致；
[0029] c、假設(shè)當 Ji > Jmax時，則 Ji = 0。
[0030] 8、根據(jù)權(quán)利要求6所述的基于數(shù)據(jù)重構(gòu)預測巖性體的方法，其特征是步驟S2_b具 體為:對各巖性體的重構(gòu)數(shù)據(jù)體進行加權(quán)處理，加權(quán)處理的計算公式如下：
[0031] Μ:?=Σ