專利名稱：地震勘探烴類檢測方法
地震勘探烴類檢測方法，是一種利用數(shù)字地震記錄，通過適當(dāng)?shù)臄?shù)字處理來直接找油找氣的方法。
關(guān)于這種方法，大慶石油管理局地球物理勘探公司，在1987年11月1日申請一個名稱為“地震勘探烴類檢測方法”的專利，其請?zhí)枮?7107616。
在1989年8月10日我公司又申請名為“地震勘探烴類檢測方法”的專利，申請?zhí)枮?9105573.8。這兩項(xiàng)專利申請的基本內(nèi)容和公式為將野外數(shù)字磁帶記錄，在地震處理機(jī)上，做基本處理，使地震記錄道變成做過動靜校正的道、迭加道、迭偏道。在同一條測線上，對上述三種道任選一種或分選后靠近炮點(diǎn)道，按申請?zhí)枮?7107616公式處理在第i個采樣點(diǎn)，設(shè)
其中，k＝0，1，2，……，n。n是小于道的采樣點(diǎn)個數(shù)的正整數(shù)。
按申請?zhí)枮?9105573.8專利申請公式處理在第i個采樣點(diǎn)，設(shè)
其中，k＝0，1，2，……，n。n是小于道的采樣點(diǎn)個數(shù)的正整數(shù)。
在所選道的時窗內(nèi)，對每個采樣點(diǎn)都計算出λ′(t)。對所有的λ′(t)≥1的值分成等級，每一級用一個機(jī)器符號代表，形成符號道，把符號道輸出出來就是烴類檢測圖。
按申請?zhí)枮?9105573.8的專利申請的另一個公式處理在第i個采樣做
d＇高
與
d＇低
的符號比較，不妨設(shè)
其中，K＝0，1，2，……，n。n是小于道的采樣點(diǎn)個數(shù)的正整數(shù)。
取λ″(ti)的符號。若λ″(ti)＜0，用地震機(jī)上任一符號表示;λ″(ti)≥0，用空白表示。對每個采樣點(diǎn)都這做，形成符號道，把符號道輸出出來，也形成烴類檢測圖。
本發(fā)明仍根據(jù)專利申請?zhí)枮?7107616，89105573.8的“地震勘探烴類檢測方法”的理論依據(jù)，提出另一“地震勘探烴類檢測方法”的申請。
根據(jù)方法的理論依據(jù)，可畫出如附

圖1、2、3、4的示意圖。設(shè)地下有一儲油氣水平單層，它的反射模型如圖1所示。(射線[1]，[2]，[3]都是垂直層面并重合在一起的，為說明方便，而把它們分開畫出。)入射波方向[1]，波形函數(shù)b(t);層頂反射波方向[2]，波形函數(shù)b1(t);層底反射波方向[3]，波形函數(shù)b2(t)。b(t)、b1(t)、b2(t)的頻譜為B(f);B1(f);B2(f)。
圖2是水平單層頂、底的反射波形。[5]是b2(t)的形狀，[7]是b1(t)的形狀，[6]是采樣點(diǎn)ti的位置。
設(shè)此層反射記錄道主要是b1(t)、b2(t)組合的復(fù)合振動。反射記錄為它們的代數(shù)和。則d(t)＝b1(t)+b2(t)。
假設(shè)，在做地震處理時，在帶通濾波B(f1，f2，f3，f4)中，選f1＝5，f2＝7，f3＝12，f4＝15時，地震道被這樣低頻帶通濾波器濾波后d(t)變?yōu)閐′低(t)，b1(t)被濾波后變?yōu)閎′1(t)，b2(t)被濾波后變?yōu)閎′2(t)。圖3中[8]，[10]分別是b′2(t)，b′1(t)的波形。[9]是采樣點(diǎn)ti的位置。在帶通濾波B(f1，f2，f3，f4)中，選f1＝5，f2＝7，f3＝55，f4＝60，地震道被這樣寬頻帶濾波器濾波后d(t)變?yōu)閐′全(t)。
假設(shè)，在地震資料處理時，在帶通濾波B(f1，f2，f3，f4)中，選f1＝12，f2＝15，f3＝55，f4＝60時，地震道被這樣躲開10HZ附近的低頻的寬頻帶濾波器濾波后d(t)變?yōu)閐′高(t)。同時，b1(t)被濾波后變?yōu)閎″1(t)，b2(t)被濾波后變?yōu)閎″2(t)。
圖4中，[11]，[13]分別為b″2(t)，b″1(t)的波形。[12]是采樣點(diǎn)ti的位置。
新方法的理論分析。從示意圖3和圖4可看出在ti附近，d′低(t)與d′高(t)值的符號相反。在ti附近，設(shè)λ-(t)＝｜d′低(t)-d′高(t)｜＝｜d′低(t)｜+｜d′高(t)｜則有λ-(t)＞｜d′低(t)｜，λ-(t)＞｜d′高(t)｜，λ-(t)＞｜d′全(t)｜。
在ti附近，設(shè)λ+(t)＝｜d′低(t)+d′高(t)｜則有λ+(t)小于｜d′低(t)｜，｜d′高(t)｜中的至少一個值。如果，λ+(t)＞｜d′低(t)｜又λ+(t)＞｜d′高(t)｜，則在ti附近的層肯定不是油氣層。
在示意圖ti附近，用λ-(t)，λ+(t)，｜d′低(t)｜，｜d′高(t)｜，｜d′全(t)｜做比較(也可把所有絕對值符號都去掉做值比較)，能把儲油氣層做為異常帶找出來。
d′高(ti-j)是選的所有道被躲開10HZ附近低頻的寬帶濾波器濾波后的第i-j個采樣點(diǎn)的值。
例如，對迭加道，在帶通濾波為B(f1，f2，f3，f4)中，選f1＝12，f2＝15，f3＝55，f4＝60，對所選的每道做濾波。
d′全(ti-j)是選的所有道，在帶通濾波器B(f1，f2，f3，f4)中，選如f1＝5，f2＝7，f3＝55，f4＝60，地震道被這樣寬頻帶濾波器濾波后的第i-j個采樣點(diǎn)的值。
d′低(ti-j)是選的所有道被低頻帶通濾波器濾波后的第i-j個采樣點(diǎn)的值例如，對迭加道，在帶通濾波B(f1，f2，f3，f4)中，選f1＝6，f2＝8，f3＝12，f4＝14，對每個道做濾波。
根據(jù)以上分析，該方法可這樣描述。本方法的內(nèi)容是將野外數(shù)字磁帶記錄，在地震處理機(jī)上，做基本處理，使地震記錄道變成做過動靜校正的道、迭加道、迭偏道。在同一條測線上，對上述三種道任選一種或分選后靠近炮點(diǎn)道，做如下處理在第i個采樣點(diǎn)，設(shè)
其中，K＝0，1，2，……，n。n是小于道的采樣點(diǎn)個數(shù)的正整數(shù)。
(1)若λ-(ti)＞a1λ1(ti)，λ-(ti)＞a2λ2(ti)或λ1(ti)＞a3λ3(ti)，λ2(ti)＞a4λ3(ti)，或做λ4與λ5的符號比較符號相異時，則在ti采樣點(diǎn)保留原振幅值(用其他非零值也可)。否則，把在ti采樣點(diǎn)的值用零值代替。(也可反過來，當(dāng)前者成立時，在采樣點(diǎn)處用零值代替;后者成立用非零值代替或保留原振幅值。)在選定的時窗內(nèi)，對每道采樣點(diǎn)都這樣做。然后，用波形或波形加面積辦法顯示出來，形成烴類檢測圖。
或在ti采樣點(diǎn)，若λ-(ti)＞a1λ1(ti)，λ-(ti)＞a2λ2(ti)或當(dāng)λ1(ti)＞a3λ3(ti)，λ2(ti)＞a4λ3(ti)時，在ti處采樣點(diǎn)之值用任何機(jī)器符號代替;否則，用與上一個不同的符號代替ti處的采樣點(diǎn)之值。在選定的時窗內(nèi)，對每道每個采樣點(diǎn)都這樣做，形成符號道，把符號道輸出出來，也形成烴類檢測圖。(其中a1，a2，a3，a4＞0)a1，a2，a3，a4根據(jù)已知井，看烴類顯示符合已知油氣層情況來確定。
具體做圖時，對上述符號顯示，是用機(jī)器符號“0”顯示的，另一個不同符號是用空白顯示出來的，效果好。
(2)若λ1(ti)≥b1λ+(ti)或λ2(ti)≥b2λ+(ti)(或同時成立)，則在采樣點(diǎn)ti處保留原振幅值(用其他非零值也可);否則，把在ti采樣點(diǎn)用零值代替。(也可反過來，當(dāng)前者成立時，在采樣點(diǎn)處用零值代替;后者成立用非零值代替或保留原振幅值。)在選定的時窗內(nèi)，對每道采樣點(diǎn)都這樣做。然后，用波形或波形加面積辦法顯示出來，形成烴類檢測圖(其中b1＞0，b2＞0)。
或在ti采樣點(diǎn)，若λ1(ti)≥b1λ+(ti)或λ2(ti)≥b2λ+(ti)(或同時成立)，則在ti處采樣點(diǎn)之值用任何機(jī)器符號代替;否則，用與上一個不同的符號代替ti處的采樣點(diǎn)之值。在選定的時窗內(nèi)，對每道每個采樣點(diǎn)都這樣做，形成符號道，把符號道輸出出來，也形成烴類檢測圖。
當(dāng)λ1(ti)≥b1λ+(ti)或λ2(ti)≥b2λ+(ti)成立時(或同時成立)，從b1或b2開始，把λ′1(ti)，λ′2(ti)分成若干等級，然后用機(jī)器符號代表每一級(其中，λ′1(ti)＝λ1(ti)/λ+(ti)，λ′2(ti)＝λ2(ti)/λ+(ti))，形成符號道，輸出出來，成烴類檢測圖(可取λ′1(ti)，λ′2(ti)的倒數(shù))。
最簡單辦法，把λ′1(ti)，λ′2(ti)分成一級，當(dāng)條件成立時，在第i個采樣點(diǎn)用機(jī)器符號“0”顯示;否則用空白代替在第i個采樣點(diǎn)之值。
b1，b2的選擇，要根據(jù)已知井顯示符合油氣層情況來確定。
(3)對(1)中所述，當(dāng)n＝0時效果最好，此時λ1(ti)＝｜d′低(ti)｜，λ2(ti)＝｜d′高(ti)｜，λ-(ti)＝｜d′低(ti)-d′高(ti)｜，λ4(ti)＝d′低(ti)，λ5(ti)＝d′高(ti)，其他處理和顯示辦法與(1)中所述相同。
(4)對(2)中所述，當(dāng)n＝0時效果最好，此時λ+(ti)＝｜d′低(ti)+d′高(ti)｜。其他處理和顯示辦法與(2)相同。
(5)在第i個采樣點(diǎn)，設(shè)λ″1(ti)＝λ-(ti)/λ3(ti)，λ″2(ti)＝λ-(ti)/λ+(ti)(可取λ″1(ti)，λ″2(ti)的倒數(shù))當(dāng)條件λ″1(ti)≥a′2或λ″1(ti)≥a′1成立時，在第i個采樣點(diǎn)之值，用非零值代替或保存原振幅值;否則，用零值代替第i個采樣點(diǎn)之值。然后，用波形或波形加面積顯示，也成烴類檢測圖。
還可以當(dāng)條件成立時，可把λ″1(ti)從下界a′1開始，λ″2(ti)從下界a′2開始，分成若干個等級，用機(jī)器符號代表每一級，形成符號道，把符號道輸出出來，形成烴類檢測圖(其中，a′1，a′2＞0)。
最簡單的辦法，用迭加道n＝0時，當(dāng)條件成立時，用符號“0”代表第i個采樣點(diǎn)之值;否則用空白代替第i個采樣點(diǎn)之值。
如何予測地下儲油氣帶。用申請?zhí)枮?7107616，名為“地震勘探烴類檢測方法”的專利申請，所給頻率掃描的辦法，確定予測地下儲油氣帶的濾波器B(f1，f2，f3，f4)，B′(f′1，f′2，f′3，f′4)。然后確定d′高(t)相應(yīng)的濾波B(f″1，f″2，f″3，f″4)(也可按相似辦法，先確定B(f″1，f″2，f″3，f″4)，B′(f′1，f′2，f′3，f′4)，后確定B(f1，f2，f3，f4)，要保證f″1≥f3，是寬頻帶濾波器。不妨設(shè)，f″1＝f3，f″2＝f4，f″3＝f′3，f″4＝f′4。再用三個濾波器，對所有道做濾波，計算λ-(ti)，λ+(ti)，λ1(ti)，λ2(ti)，λ3(ti)，λ4(ti)，λ5(ti)。然后，按(1)，(2)，(3)，(4)，(5)所述辦法，把予測的油氣帶確定出來。
如何看檢測圖。圖的上方所示是X方向CDP點(diǎn)數(shù);縱向是時間t。根據(jù)已知井，已知儲油氣帶，烴類檢測圖所示儲油氣剖面，符合的最佳，圖的顯示強(qiáng)度最強(qiáng)或符號顯示最密集的為最佳。
本發(fā)明的顯示圖，特別用波形或波形加面積的顯示，是比較直觀的。對確定相應(yīng)迭加剖面上的儲油氣帶，比以前報的二項(xiàng)專利申請方便，但分辨能力差些。
圖5是宋站139.0線的部分水平迭加剖面，A處離測線的垂直方向的東2井931米，CDP點(diǎn)數(shù)為765。
圖6是宋站139.0線用波形加面積的烴類檢測圖。
從圖5看，東2井在箭頭S，T所指的二層儲油氣帶都得到了工業(yè)氣流。
在圖6上S和T箭頭所指異常帶就是儲油氣帶在烴類檢測圖上的顯示。這二個異常帶的位置、形態(tài)、范圍都符合地下兩個儲油氣帶。
權(quán)利要求
1.一種利用數(shù)字地震記錄資料，直接找油氣的地震勘探烴類檢測方法，先將野外數(shù)字磁帶記錄，在地震機(jī)上做基本處理，使地震記錄道變成做過動靜校正后的道，迭加道，迭偏道。其特征在于在同條測線上，任選上述三種道的一種或分選后靠近炮點(diǎn)道，按如下公式處理在第i個采樣點(diǎn)，設(shè)
其中，k=0，1，2……，n。n是小于每道采樣點(diǎn)個數(shù)的正整數(shù)。d′高(ti-j)是選的所有道被躲開10HZ附近低頻的寬帶通濾波器濾波后的第i-j個采樣點(diǎn)的值；例如，對迭加道，在帶通濾波B(f1，f2，f3，f4)中，選f1=12，f2=15，f3=55，f4=60，對每個道做濾波。d′全(ti-j)是選的所有道被寬頻帶通濾波器濾波后的第i-j個采樣點(diǎn)的值例如，對迭加道，在帶通濾波B(f1，f2，f3，f4)中，選f1=6，f2=8，f3=55，f4=60，對每個道做濾波。d′低(ti-j)是選的所有道被低頻帶通濾波器濾波后的第i-j個采樣點(diǎn)的值例如，對迭加道，在帶通濾波B(f1，f2，f3，f4)中，選f1=6，f2=8，f3=12，f=14，對每個道做濾波。對λ-(ti)，λ1(ti)，λ2(ti)，λ3(ti)做值的比較(也可把公式中的絕對值號都去掉，做值比較)加判斷條件，把相對迭加剖面的儲油氣帶剖面找出來。若λ-(ti)＞a1λ1(ti)，λ-(ti)＞a2λ2(ti)或λ1(ti)＞a3λ3(ti)，λ2(ti)＞a4λ3(ti)，或做λ4與λ5的符號比較，符號相異時，則在ti采樣點(diǎn)保留原振幅值(用其他非零值也可)。否則，把在ti采樣點(diǎn)的值用零值代替。(也可反過來，當(dāng)前者成立時，在采樣點(diǎn)處用零值代替；后者成立用非零值代替或保留原振幅值。)在選定的時窗內(nèi)，對每道采樣點(diǎn)都這樣做。然后，用波形或波形加面積辦法顯示出來，形成烴類檢測圖?；蛟趖i采樣點(diǎn)，若λ-(ti)＞a1λ1(ti)，λ-(ti)＞a2λ2(ti)或當(dāng)λ1(ti)＞a3λ3(ti)，λ2(ti)＞a4λ3(ti)時，在ti處采樣點(diǎn)之值用任何機(jī)器符號代替；否則，用與上一個不同的符號代替ti處的采樣點(diǎn)之值。在選定的時窗內(nèi)，對每道每個采樣點(diǎn)都這樣做，形成符號道，把符號道輸出出來，也形成烴類檢測圖。(其中a1，a2，a3，a4＞0)a1，a2，a3，a4根據(jù)已知井，看烴類顯示符號已知油氣層情況來確定。具體做圖時，對上述符號顯示，是用機(jī)器符號“O”顯示的，另一個不同符號是用空白顯示出來的，效果好。
2.據(jù)權(quán)利要求1所述的地震勘探烴類檢測方法，其特征在于在同條測線上，任選權(quán)利要求1所述的四種道的一種，按如下辦法處理設(shè)λ′1(ti)＝λ1(ti)/λ+(ti)，λ′2(ti)＝λ2(ti)/λ+(ti)(可取λ′1(ti)，λ′2(ti)的倒數(shù))若λ1(ti)≥b1λ+(ti)或λ2(ti)≥b2λ+(ti)(或同時成立)，則在采樣點(diǎn)ti處保留原振幅值(用其他非零值也可);否則，把在ti采樣點(diǎn)用零值代替。(也可反過來，當(dāng)前者成立時，在采樣點(diǎn)處用零值代替;后者成立用非零值代替或保留原振幅值。)在選定的時窗內(nèi)，對每道采樣點(diǎn)都這樣做。然后，用波形或波形加面積辦法顯示出來，形成烴類檢測圖(其中b1＞0，b2＞0)?；蛟趖i采樣點(diǎn)，若λ1(ti)≥b1λ+(ti)，λ2(ti)≥b2λ+(ti)(或同時成立)，則在ti處采樣點(diǎn)之值用任何機(jī)器符號代替;否則，用與上一個不同的符號代替ti處的采樣點(diǎn)之值。在選定的時窗內(nèi)，對每道每個采樣點(diǎn)都這樣做，形成符號道，把符號道輸出出來，也形成烴類檢測圖。當(dāng)λ1(ti)≥b1λ+(ti)或λ2(ti)≥b2λ+(ti)成立時(或同時成立)，從b1或b2開始，把λ′1(ti)，λ′2(ti)分成若干等級，然后用機(jī)器符號代表每一級，形成符號道，輸出出來，成烴類檢測圖。最簡單辦法，把λ′1(ti)，λ′2(ti)分成一級，當(dāng)條件成立時，在第i個采樣點(diǎn)用機(jī)器符號“0”顯示;否則用空白代替在第i個采樣點(diǎn)之值。b1，b2的選擇，要根據(jù)已知井顯示符合油氣層情況來確定。
3.據(jù)權(quán)利要求1所述的地震勘探烴類檢測方法，其特征在于當(dāng)n＝0時效果最好，此時λ1(ti)＝｜d′低(ti)｜，λ2(ti)＝｜d′高(ti)｜，λ-(ti)＝｜d′低(ti)-d′高(ti)｜，λ4(ti)＝d′低(ti)，λ5(ti)＝d′高(ti)，其他處理和顯示辦法與權(quán)利要求1.所述相同。
4.據(jù)權(quán)利要求2的地震勘探烴類檢測方法，其特征在于當(dāng)n＝0時效果最好，此時λ+(ti)＝｜d′低(ti)+d′高(ti)｜。其他處理和顯示辦法與權(quán)利要求2.相同。
5.據(jù)權(quán)利要求1所述的地震勘探烴類檢測方法，其特征在于在同條測線上，任選權(quán)利要求1所述的四種道的一種，按如下辦法處理在第i個采樣點(diǎn)，設(shè)λ″1(ti)＝λ-(ti)/λ3(ti)，λ″2(ti)＝λ-(ti)/λ+(ti)(可取λ″1(ti)，λ″2(ti)的倒數(shù))當(dāng)條件λ″2(ti)≥a′2或λ″1(ti)≥a′1成立時，在第i個采樣點(diǎn)之值，用非零值代替或保存原振幅值;否則，用零值代替第i個采樣點(diǎn)之值。然后，用波形或波形加面積顯示，也成烴類檢測圖。還可以當(dāng)條件成立時，可把λ″1(ti)從下界a′1開始，λ″2(ti)從下界a′2開始，分成若干個等級，用機(jī)器符號代表每一級，形成符號道，把符號道輸出出來，形成烴類檢測圖(其中，a″1，a″2＞0)。最簡單的辦法，用迭加道n＝0時，當(dāng)條件成立時，用符號“0”代表第i個采樣點(diǎn)之值;否則用空白代替第i個采樣點(diǎn)之值。如何予測地下儲油氣帶。用申請?zhí)枮?7107616，名為“地震勘探烴類檢測方法”的專利申請，所給頻率掃描的辦法，確定予測地下儲油氣帶的濾波器B(f1，f2，f3，f4)，B′(f′1，f′2，f′3，f′4)。然后確定d′高(t)相應(yīng)的濾波B(f″1，f″2，f″3，f″4)(也可按相似辦法，先確定B(f″1，f″2，f″3，f″4)，B′(f′1，f′2，f′3，f′4)，后確定B(f1，f2，f3，f4)，要保證f″1≥f3，是寬頻帶濾波器。不妨設(shè)，f″1＝f3，f″2＝f4，f″3＝f′3，f″4＝f′4。再用三個濾波器，對所有道做濾波，計算λ-(ti)，λ+(ti)，λ1(ti)，λ2(ti)，λ3(ti)，λ4(ti)，λ5(ti)。然后，按權(quán)利要求1，2，3，4，5所述辦法，把予測的油氣帶確定出來。如何看檢測圖。圖的上方所示是X方向CDP點(diǎn)數(shù);縱向是時間t。根據(jù)已知井，已知儲油氣帶，烴類檢測圖所示儲油氣剖面，符合的最佳，圖的顯示強(qiáng)度最強(qiáng)或符號顯示最密集的為最佳。
全文摘要
地震勘探烴類檢測方法，是一種用數(shù)字地震記錄來直接找油氣的方法。利用不同帶通濾波器，對做過基本處理的各種道濾波。做不同帶通濾波后的相同道上相同采樣點(diǎn)上的值做比較和判別，在選定的時窗內(nèi)對所有道都這樣做，把地下儲油氣帶相應(yīng)迭加剖面的位置找出并顯示出來，形成烴類檢測圖。此圖能把相對迭加剖面上儲油氣帶整個形態(tài)顯示出來。
文檔編號G01V1/28GK1056171SQ91101518
公開日1991年11月13日 申請日期1991年3月12日 優(yōu)先權(quán)日1991年3月12日
發(fā)明者黃克有 申請人:大慶石油管理局地球物理勘探公司