用于確定斷層展布格局的地震資料解釋方法和設(shè)備的制造方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及地震資料解釋領(lǐng)域，具體地，涉及一種用于確定斷層展布格局的地震資料解釋方法和設(shè)備。
【背景技術(shù)】
[0002]我國大陸經(jīng)歷多次板塊構(gòu)造旋回的疊加，地質(zhì)構(gòu)造復(fù)雜，尤其是斷層構(gòu)造發(fā)育。在能源、工程勘探中，首先需要對斷層構(gòu)造展布進行確定，以便科學(xué)確定能源開發(fā)、工程布置方案。地震資料解釋是用于確定構(gòu)造展布的一種重要方法。目前，常規(guī)的斷層展布地震解釋方法主要是:首先對目的層反射波進行層位追蹤，然后通過手工解釋地震剖面，畫出斷面形狀，層位與斷面的交點進行平面投影，通過多個交點圈出斷層，獲得斷層展布格局。這種方法需要較大的人工解釋工作，而且屬于一種二維的地震解釋方法，只能對較大的構(gòu)造有較好效果，工作效率較低。為了提高斷層構(gòu)造展布格局的解釋效率，一般解釋方法是通過地震屬性提取水平切片或者沿層切片。該方法需要兩個大步驟，一是地震屬性的計算，二是目的層地震反射波追蹤。地震屬性體(例如，方差體、地震振幅等參數(shù))的計算往往涉及較大的計算量，對于普通的計算機，則需要較長的時間，在高性能的機器才能獲得較短的計算時間。同時，該屬性體的計算結(jié)果主要是幫助人工識別斷點，在進行目的層地震反射波追蹤中，主要是人工解釋為主，因此在實際應(yīng)用中存在一定的局限性。

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

[0003]本發(fā)明的目的是提供一種能夠快速確定斷層展布格局的地震資料解釋方法和設(shè)備。
[0004]為了實現(xiàn)上述目的，本發(fā)明提供一種用于確定斷層展布格局的地震資料解釋方法。該方法包括:步驟S1:對目的層反射波進行全三維層位自動追蹤，獲得原始層位信息；步驟S2:對所述原始層位信息進行平滑，獲得平滑層位信息；步驟S3:基于所述原始層位信息與所述平滑層位信息之間的差，得出反射波時間差信息；以及步驟S4:確定所述反射波時間差信息的各個波峰極值點和各個波谷極值點，并對所述各個波峰極值點進行擬合以及對所述各個波谷極值點進行擬合，得出斷層交線。
[0005]優(yōu)選地，采用互相關(guān)法對所述目的層反射波進行所述全三維層位自動追蹤，獲得所述原始層位信息。
[0006]優(yōu)選地,所述步驟SI包括:步驟SlOl:在一地震剖面上,確定所述目的層反射波在基準地震道上的位置，該位置作為種子點；步驟S102:設(shè)定一時窗；步驟S103:在所述地震剖面上的目標(biāo)地震道上，按照設(shè)定的采樣間隔移動所述時窗，并確定對應(yīng)的相關(guān)系數(shù)，其中，該相關(guān)系數(shù)表示在所述時窗未移動時所述基準地震道上被包括在該時窗內(nèi)的波形與在所述時窗移動之后所述目標(biāo)地震道上被包括在該時窗內(nèi)的波形之間的相關(guān)性；步驟S104:判斷所述相關(guān)系數(shù)是否滿足設(shè)定的閾值條件，在所述相關(guān)系數(shù)不滿足所述閾值條件的情況下，返回步驟S103繼續(xù)進行，直到所確定出的相關(guān)系數(shù)滿足所述閾值條件為止；在所述相關(guān)系數(shù)滿足所述閾值條件的情況下，進行步驟S105:基于所述采樣間隔和與該相關(guān)系數(shù)對應(yīng)的移動次數(shù)，確定所述種子點在所述目標(biāo)地震道上的移動量，以確定所述目的層反射波在所述目標(biāo)地震道上的位置；步驟S106:判斷是否已確定出所述目的層反射波在所述地震剖面上的所有地震道上的位置，若是，進行步驟S108:基于所述目的層反射波在所述地震剖面上的各個地震道上的位置，得出所述原始層位信息；若不是，進行步驟S107:更換目標(biāo)地震道，并復(fù)位所述時窗，之后返回步驟S103繼續(xù)進行，直到確定出所述目的層反射波在所述地震剖面上的所有地震道上的位置為止。
[0007]優(yōu)選地，所述步驟SI包括:步驟SlOl:在一地震剖面上，確定目的層反射波在基準地震道上的位置，該位置作為種子點；步驟S102:設(shè)定一時窗；步驟S103:在所述地震剖面上的目標(biāo)地震道上，按照設(shè)定的采樣間隔移動所述時窗，并確定對應(yīng)的相關(guān)系數(shù)，其中，該相關(guān)系數(shù)表示在所述時窗未移動時所述基準地震道上被包括在該時窗內(nèi)的波形與在所述時窗移動之后所述目標(biāo)地震道上被包括在該時窗內(nèi)的波形之間的相關(guān)性；步驟S104:判斷所述時窗在所述目標(biāo)地震道上的移動次數(shù)是否已達到最大移動次數(shù)，其中，該最大移動次數(shù)是基于設(shè)定的最大傾角時差和所述采樣間隔來確定的，在所述移動次數(shù)未達到所述最大移動次數(shù)的情況下，返回步驟S103繼續(xù)進行，直到所述移動次數(shù)達到所述最大移動次數(shù)為止；在所述移動次數(shù)達到所述最大移動次數(shù)的情況下，進行步驟S105:基于所述采樣間隔和與最大相關(guān)系數(shù)對應(yīng)的移動次數(shù)，確定所述種子點在所述目標(biāo)地震道上的移動量，以確定所述目的層反射波在所述目標(biāo)地震道上的位置；步驟S106:判斷是否已確定出所述目的層反射波在所述地震剖面上的所有地震道上的位置，若是，進行步驟S108:基于所述目的層反射波在所述地震剖面上的各個地震道上的位置，得出所述原始層位信息；若不是，進行步驟S107:更換目標(biāo)地震道，并復(fù)位所述時窗，之后返回步驟S103繼續(xù)進行，直到確定出所述目的層反射波在所述地震剖面上的所有地震道上的位置為止。
[0008]優(yōu)選地，該方法還包括:在所述移動次數(shù)達到所述最大移動次數(shù)的情況下，先進行步驟S109:判斷所述最大相關(guān)系數(shù)是否滿足設(shè)定的閾值條件，在所述最大相關(guān)系數(shù)滿足所述閾值條件的情況下，再進行步驟S105。
[0009]優(yōu)選地，采用最小二乘平滑法對所述原始層位信息進行平滑，獲得所述平滑層位信息。
[0010]優(yōu)選地，采用以下方法之一對所述各個波峰極值點進行擬合以及對所述各個波谷極值點進行擬合:最小二乘擬合法或螞蟻算法。
[0011]本發(fā)明還提供一種用于確定斷層展布格局的地震資料解釋設(shè)備。該設(shè)備包括:用于對目的層反射波進行全三維層位自動追蹤，獲得原始層位信息的裝置；用于對所述原始層位信息進行平滑，獲得平滑層位信息的裝置；用于基于所述原始層位信息與所述平滑層位信息之間的差，得出反射波時間差信息的裝置；以及用于確定所述反射波時間差信息的各個波峰極值點和各個波谷極值點，并對所述各個波峰極值點進行擬合以及對所述各個波谷極值點進行擬合，得出斷層交線的裝置。
[0012]優(yōu)選地，對所述目的層反射波進行所述全三維層位自動追蹤，獲得所述原始層位信息是采用互相關(guān)法進行的。
[0013]優(yōu)選地，所述用于對目的層反射波進行全三維層位自動追蹤，獲得原始層位信息的裝置包括:用于在一地震剖面上，確定所述目的層反射波在基準地震道上的位置的第一裝置，該位置作為種子點；用于設(shè)定一時窗的第二裝置；用于在所述地震剖面上的目標(biāo)地震道上，按照設(shè)定的采樣間隔移動所述時窗，并確定對應(yīng)的相關(guān)系數(shù)的第三裝置，其中，該相關(guān)系數(shù)表示在所述時窗未移動時所述基準地震道上被包括在該時窗內(nèi)的波形與在所述時窗移動之后所述目標(biāo)地震道上被包括在該時窗內(nèi)的波形之間的相關(guān)性；用于判斷所述相關(guān)系數(shù)是否滿足設(shè)定的閾值條件的第四裝置；用于在所述相關(guān)系數(shù)不滿足所述閾值條件的情況下，重復(fù)運行所述第三裝置，直到所確定出的相關(guān)系數(shù)滿足所述閾值條件為止的第五裝置；用于在所述相關(guān)系數(shù)滿足所述閾值條件的情況下，基于所述采樣間隔和與該相關(guān)系數(shù)對應(yīng)的移動次數(shù)，確定所述種子點在所述目標(biāo)地震道上的移動量，以確定所述目的層反射波在所述目標(biāo)地震道上的位置的第六裝置；用于判斷是否已確定出所述目的層反射波在所述地震剖面上的所有地震道上的位置的第七裝置；用于在尚未確定出所述目的層反射波在所述地震剖面上的所有地震道上的位置的情況下，更換目標(biāo)地震道，并復(fù)位所述時窗，之后重復(fù)運行所述第三裝置，直到確定出所述目的層反射波在所述地震剖面上的所有地震道上的位置為止的第八裝置；以及用于在已確定出所述目的層反射波在所述地震剖面上的所有地震道上的位置的情況下，基于所述目的層反射波在所述地震剖面上的各個地震道上的位置，得出所述原始層位信息的第九裝置。
[0014]優(yōu)選地，所述用