曾華會(huì) 王小衛(wèi) 蘇 勤 雍運(yùn)動(dòng) 劉 桓 孟會(huì)杰

(中國(guó)石油勘探開發(fā)研究院西北分院，甘肅蘭州730030)

1 研究背景

中國(guó)油氣勘探與開發(fā)的持續(xù)進(jìn)行，使得在復(fù)雜地震地質(zhì)條件下尋找油氣藏已成為目前及今后陸上勘探的必然趨勢(shì)。隨著地震勘探難度日益增大，地表?xiàng)l件越來(lái)越復(fù)雜，既有可用炸藥震源的山地、沙漠區(qū)，又有村莊、城鎮(zhèn)、工業(yè)園、江河、水庫(kù)等“禁炮區(qū)”，因此野外采集常采用多種類型震源(如炸藥震源、可控震源、氣槍震源等)混合施工?？煽卣鹪淳哂邪踩h(huán)保、使用高效、方便操控等優(yōu)勢(shì)，它的應(yīng)用和發(fā)展不僅符合國(guó)家安全需要，更是環(huán)保要求；既能降低采集施工成本，還能在保證資料品質(zhì)基礎(chǔ)上大幅度提高現(xiàn)場(chǎng)施工效率，縮短項(xiàng)目周期。

在復(fù)雜障礙區(qū)應(yīng)用混合震源采集面臨諸多挑戰(zhàn)。受激發(fā)方式、激發(fā)井深、接收儀器、地表地質(zhì)條件等客觀因素影響，采集所得地震資料在子波、振幅能量、相位、頻率、信噪比、時(shí)差等方面均存在較大差異[1]。例如，炸藥震源子波是最小相位，能量較強(qiáng)，頻率較高； 可控震源子波為零相位，能量、頻率可以人為控制； 重錘震源子波為最小相位，能量較低，地表因素影響較嚴(yán)重[2]。這些差異使不同震源采集的資料在處理時(shí)無(wú)法同相疊加，嚴(yán)重降低了資料的信噪比和分辨率。如果采用常規(guī)處理流程，地震資料在兩種震源拼接處會(huì)存在較大差異，地震剖面存在波組特征較差、信噪比和分辨率較低等問(wèn)題[3-7]； 混合震源地震剖面存在較大時(shí)差(圖1)，兩者不閉合給后續(xù)地震資料精細(xì)處理、解釋帶來(lái)很大影響，如分層不清晰、不同(混合)震源之間時(shí)差被誤認(rèn)為是假斷裂、波阻抗反演的精度難以保證、地震屬性參數(shù)提取困難等。

目前，實(shí)際數(shù)據(jù)處理中解決混合震源不一致性問(wèn)題的主要方法有時(shí)移掃描法、反褶積參數(shù)調(diào)整法和匹配濾波法等。匹配濾波方法通常用于陸上不同震源數(shù)據(jù)的匹配，或淺海與陸地過(guò)渡區(qū)域地震數(shù)據(jù)匹配。唐恩德等[8]認(rèn)為在處理混合震源地震資料時(shí)，要優(yōu)化模型道及對(duì)單炮記錄子波進(jìn)行最小相位化處理以消除子波差異，再通過(guò)譜均衡等地表一致性技術(shù)處理振幅差異； 莊東海等[9]強(qiáng)調(diào)，針對(duì)混合震源資料，要做好相對(duì)振幅保持、提高信噪比和資料一致性等處理，而其中的關(guān)鍵是互均衡處理； 婁兵等[10]將子波整形及匹配濾波兩種方法應(yīng)用于不同震源的拼接處理，提高了地震資料信噪比和分辨率； 王西文等[11]將基于小波變換的地震子波處理技術(shù)應(yīng)用于均衡連片工區(qū)數(shù)據(jù)在頻率、能量、信噪比等方面差異，為地震解釋和反演提供高信噪比地震資料； 宋玉龍[12]在灘淺海地區(qū)利用匹配濾波技術(shù)解決由觀測(cè)系統(tǒng)、激發(fā)、接收等因素導(dǎo)致的海陸資料之間差異，利用疊前多點(diǎn)綜合串聯(lián)匹配濾波處理技術(shù)解決炸藥與可控震源波組特征、頻率、相位及能量不一致情況，較好地解決了資料銜接時(shí)不匹配問(wèn)題； 高少武等[13-14]采用匹配濾波技術(shù)處理炸藥震源與可控震源數(shù)據(jù)的銜接，使這兩種震源資料的振幅和相位趨于一致。

但在同一工區(qū)因復(fù)雜障礙區(qū)施工困難而采用混合震源時(shí)，兩種混合震源資料的重復(fù)段疊加剖面由于炮檢距范圍、頻率、疊加速度、靜校正等多種因素的影響，信噪比存在較大差異，難以準(zhǔn)確求取匹配算子； 另外，目前常規(guī)匹配濾波處理方法只是計(jì)算一個(gè)濾波因子，而不考慮地震資料信噪比。

圖1 可控震源與炸藥震源混合區(qū)疊加剖面

因此，本文在充分考慮信號(hào)和噪聲的基礎(chǔ)上，開展了混合震源高精度匹配處理技術(shù)研究，通過(guò)綜合應(yīng)用疊前平均振幅譜匹配、混合震源相位匹配、剩余時(shí)差校正等處理技術(shù)，消除了混合震源間存在的振幅、相位、頻率差異及時(shí)差不一致等問(wèn)題，同時(shí)兼顧地震子波的變化。將該方法應(yīng)用于四川盆地大川中北斜坡射洪—鹽亭三維復(fù)雜障礙區(qū)的地震資料處理，提高了不同震源資料匹配處理的成像品質(zhì)和地震資料信噪比，消除了假斷裂和假構(gòu)造等現(xiàn)象，可為精細(xì)地震解釋和反演提供高保真、高信噪比的地震資料。

2 研究區(qū)概況

研究區(qū)處于四川盆地川中古隆起平緩構(gòu)造區(qū)儀隴構(gòu)造群，以山地—丘陵地貌為主。區(qū)內(nèi)河流較發(fā)育，包括涪江、梓江及其支流，涪江江面寬為300～1000m，梓江江面寬為80～500m； 涉及中小型水庫(kù)12座，這些水庫(kù)多為當(dāng)?shù)匾患?jí)飲水水源，禁止布設(shè)水下檢波器； 區(qū)內(nèi)主要障礙包括：中華涪江濕地走廊自然保護(hù)區(qū)、射洪和鹽亭縣城及城鎮(zhèn)、水系、天然氣管線和水渠等地面地下設(shè)施、風(fēng)景區(qū)等。工區(qū)內(nèi)復(fù)合障礙面積大、范圍廣，嚴(yán)重影響觀測(cè)系統(tǒng)設(shè)計(jì)和物理點(diǎn)的均勻布設(shè)，因此采用可控震源與炸藥震源混合激發(fā)施工。

圖2 研究區(qū)障礙分布示意圖

3 混合震源高精度匹配處理方法

3.1 疊前平均振幅譜匹配處理技術(shù)

常規(guī)振幅匹配處理方法只計(jì)算一個(gè)濾波因子，而很少顧及地震資料信噪比，而地震數(shù)據(jù)一般包含信號(hào)和噪聲兩部分，在設(shè)計(jì)兩種混合震源地震數(shù)據(jù)匹配算子時(shí)，理想的條件是進(jìn)行信號(hào)之間的匹配，而不是噪聲，噪聲會(huì)影響所求取的匹配算子的穩(wěn)定性[15-16]。假設(shè)有x(t)和y(t)兩道地震數(shù)據(jù)，分別由信號(hào)sx(t)、sy(t)和噪ny(t)、ny(t)聲組成

x(t)=sx(t)+ny(t)

(1)

y(t)=sy(t)+ny(t)

(2)

其傅里葉變換為

X(f)=Sx(f)+Nx(f)

(3)

Y(f)=Sy(f)+Ny(f)

(4)

x(t)的自相關(guān)譜可寫為

(5)

式中帶上劃線為共軛復(fù)數(shù)。

同理可得x(t)與y(t)的互相關(guān)譜表達(dá)式

(6)

若是兩組地震道數(shù)據(jù)，那么一組x(t)的平均自相關(guān)譜用數(shù)學(xué)期望可表示為

(7)

在平均意義下，信號(hào)與噪聲是不相關(guān)的，則有

(8)

同理可得

(9)

(10)

進(jìn)一步假設(shè)兩種數(shù)據(jù)的噪聲也是不相關(guān)的，那么平均互相關(guān)譜可變?yōu)?/p>

(11)

從該式可看出，平均互相關(guān)譜的相位等于兩組地震道數(shù)據(jù)的信號(hào)譜的相位差。

定義x(t)的噪聲與信號(hào)比值為

(12)

(13)

進(jìn)而可得

(14)

假設(shè)兩種數(shù)據(jù)的噪聲與信號(hào)比值相同，則有

(15)

定義兩組地震道的譜相關(guān)為

(16)

可得

(17)

若γ已知，那么可得到x(t)的信號(hào)譜為

(18)

顯然有

(19)

同樣可得到噪聲譜

(20)

為促使y(t)向x(t)匹配，定義匹配算子為

(21)

匹配算子包含兩項(xiàng)： 前一項(xiàng)是整形濾波器； 后一項(xiàng)是為了壓制噪聲，在噪聲較大的地震資料上也能得到穩(wěn)定的匹配算子。顯然，該式是匹配算子的振幅譜計(jì)算公式，取得的實(shí)際效果如圖3所示。

3.2 混合震源相位匹配處理技術(shù)

當(dāng)炸藥震源和可控震源兩種混合震源聯(lián)合處理時(shí)，由于兩種激發(fā)方式獲得的地震資料相位不統(tǒng)一，影響疊加剖面成像效果，降低了地震資料信噪比和分辨率[17-18]。最小平方脈沖反褶積、預(yù)測(cè)反褶積方法假設(shè)子波為最小相位，而可控震源子波更接近零相位，因此在地震處理前，需將可控震源子波最小相位化，具體步驟如下：

(1)通過(guò)可控震源掃描頻率、掃描長(zhǎng)度、坡度及掃描類型等參數(shù)，得到可控震源掃描信號(hào)；

(2)對(duì)可控震源掃描信號(hào)進(jìn)行自相關(guān)計(jì)算，得到可控震源的零相位子波；

(3)利用可控震源零相位子波，提取可控震源的最小相位化轉(zhuǎn)換算子；

(4)將轉(zhuǎn)換算子在疊前應(yīng)用于可控震源單炮，即完成相位匹配； 再將零相位可控震源子波地震數(shù)據(jù)轉(zhuǎn)化為最小相位震源子波地震數(shù)據(jù)。

圖4為可控震源子波最小相位化前、后單炮對(duì)比?？梢?jiàn)通過(guò)可控震源最小相位化處理，一方面較好地解決了與炸藥震源子波不一致問(wèn)題，即相位匹配后兩種震源都為最小相位子波，形成同相疊加，疊加剖面連續(xù)性變好，能量加強(qiáng)； 另一方面，最小相位化處理后可控震源單炮初至更光滑，更利于初至拾取，為后續(xù)近地表速度模型反演和靜校正計(jì)算奠定基礎(chǔ)。

圖3 同一位置炸藥震源與可控震源對(duì)比 (a)炸藥震源單炮； (b)可控震源振幅匹配前單炮； (c)可控震源振幅匹配后單炮

圖4 可控震源最小相位化前(a)、后(b)的單炮對(duì)比

3.3 疊前剩余時(shí)差校正方法

由于可控震源與井炮震源在激發(fā)深度、激發(fā)巖性、激發(fā)機(jī)制、子波形態(tài)、工作參數(shù)上存在較大差異，使兩者激發(fā)的子波在振幅、相位、頻率、時(shí)差上存在較大差異，通過(guò)平均振幅譜匹配和混合震源相位匹配兩項(xiàng)處理技術(shù)，較好地消除了混合震源在振幅、頻率、相位上的差異。但由于可控震源是在地面激發(fā)，而炸藥震源一般選擇在地下高速層頂界以下激發(fā)，可控震源掃描信號(hào)下傳過(guò)程中，受近地表低降速帶和吸收衰減的影響[19]，與炸藥震源之間仍存在剩余時(shí)差的影響，若剩余時(shí)差過(guò)大，會(huì)影響速度分析和同相疊加的精度，直接影響地震資料同相疊加的效果[2]。

疊前剩余時(shí)差校正主要是為了解決疊前平均振幅譜匹配處理和混合震源相位匹配處理后，不同震源在地震剖面拼接處仍存在較大時(shí)差的一種方法，具體做法是先對(duì)工區(qū)炸藥震源井深、近地表速度、厚度等信息進(jìn)行調(diào)查[20]，如圖5所示，在此基礎(chǔ)上對(duì)剩余時(shí)差進(jìn)行量化分析。

可控震源和炸藥震源不考慮炮檢距的影響，以垂直入射為例進(jìn)行時(shí)差分析，可得到剩余時(shí)差為

(22)

式中：h0、v0分別為低速帶厚度和速度；h1、v1對(duì)應(yīng)降速帶厚度和速度；h2為炸藥震源距高速層頂界距離；v2為高速層速度。

在研究區(qū)內(nèi)選取混合震源附近16口微測(cè)井進(jìn)行分析，根據(jù)微測(cè)井調(diào)查和炸藥震源井深信息，按式(22)可求得剩余時(shí)差約為14ms(圖6)。

除了受近地表因素影響，剩余時(shí)差還受炮檢距影響。因此，為提高剩余時(shí)差分析精度，在以上剩余時(shí)差分析的基礎(chǔ)上，再對(duì)炸藥震源和可控震源數(shù)據(jù)分別進(jìn)行疊加，在剖面上選取兩種混合震源的拼接處，分別求取可控震源和炸藥震源各自的疊加道數(shù)據(jù)xi(t)和yi(t)(道號(hào)i=1，2，…，N，其中N為疊加段的道數(shù))。

因?yàn)橹丿B部分是地下同一區(qū)域地層的相同反射信息，在地震剖面上具有相同的構(gòu)造形態(tài)，兩種地震記錄應(yīng)該具有良好的相關(guān)性。而互相關(guān)函數(shù)是在時(shí)間域比較兩個(gè)信號(hào)相似程度的一種方法，互相關(guān)結(jié)果可定量描述兩個(gè)信號(hào)的時(shí)間延遲，通過(guò)互相關(guān)分析，能定量求取可控震源和炸藥震源的剩余時(shí)差(圖7)約為18ms。通過(guò)疊后剩余時(shí)差統(tǒng)計(jì)，將該時(shí)差校正到可控震源的所有疊前數(shù)據(jù)上，較好地解決了兩種震源系統(tǒng)時(shí)差。

從剩余時(shí)差校正前、后的水平疊加對(duì)比剖面(圖8)可看出，經(jīng)過(guò)剩余時(shí)差校正處理，混合震源重疊區(qū)域疊加剖面信噪比和分辨率都得到顯著提升，侏羅系大安寨標(biāo)志層位不連續(xù)、混合震源時(shí)差誤認(rèn)為的假斷裂現(xiàn)象分別得到改善和厘清，有利于后續(xù)精細(xì)構(gòu)造解釋和圈閉準(zhǔn)確落實(shí)。

圖5 近地表的地質(zhì)結(jié)構(gòu)示意圖

圖6 可控震源和炸藥震源時(shí)差統(tǒng)計(jì)直方圖

圖7 炸藥震源和可控震源互相關(guān)曲線

圖8 剩余時(shí)差校正前(a)、后(b)的水平疊加剖面對(duì)比

4 實(shí)際資料應(yīng)用效果

川中古隆起北斜坡是海、陸相主要有利勘探領(lǐng)域和目標(biāo)區(qū)帶的疊合區(qū)，具有多層系立體勘探優(yōu)勢(shì)。本次實(shí)際數(shù)據(jù)取自大川中射洪—鹽亭三維工區(qū)，滿覆蓋面積為2052.1km2，跨越四川省遂寧、綿陽(yáng)、德陽(yáng)及南充四市，共涉及105個(gè)鄉(xiāng)鎮(zhèn)(圖2)。該區(qū)屬低丘陵山地地貌，整體地形較平緩，相對(duì)高程較小，但局部高程變化較大，海拔高度一般在350～720m，地勢(shì)低洼處多為涪江及梓江河沿岸，北部地勢(shì)相對(duì)較高。

在現(xiàn)場(chǎng)觀測(cè)系統(tǒng)方面，由于受到射洪縣城及工業(yè)園區(qū)、鹽亭縣城、涪江及中華涪江濕地走廊自然保護(hù)區(qū)等的影響，另外梓江河穿越工區(qū)西南部，升鐘水渠、武引水渠，天然氣管線等貫穿工區(qū)，采用炸藥震源和可控震源混合震源進(jìn)行施工，而混合震源采集在復(fù)雜障礙區(qū)接收到的地震資料在子波、振幅能量、相位、頻率、信噪比、時(shí)差等方面均存在較大差異，若不進(jìn)行精細(xì)匹配處理，必將影響地震層位精細(xì)解釋和斷裂刻畫。

在充分考慮信號(hào)和噪聲的基礎(chǔ)上，開展混合震源高精度匹配處理，具體流程如圖9所示。在該復(fù)雜障礙區(qū)地震資料處理中，綜合應(yīng)用疊前平均振幅譜匹配、混合震源相位匹配、疊前剩余時(shí)差校正等處理技術(shù)，消除了混合震源間存在的振幅、相位、頻率、時(shí)差等不一致現(xiàn)象，提高了混合震源資料的成像品質(zhì)和剖面信噪比，為精細(xì)層位解釋、斷裂解釋、疊前反演和儲(chǔ)層預(yù)測(cè)提供了高保真地震資料。

圖10為可控震源單炮最小相位化處理前、后初至拾取對(duì)比圖，可見(jiàn)最小相位化處理后更利于準(zhǔn)確拾取初至、近地表速度模型反演及層析靜校正計(jì)算。

圖9 混合震源高精度匹配處理流程

圖11為混合震源高精度匹配處理前、后水平疊加對(duì)比剖面，可見(jiàn)匹配前可控震源振幅能量明顯強(qiáng)于炸藥震源資料，且兩種混合震源存在系統(tǒng)時(shí)差，在淺層侏羅系大安寨(J1dn)標(biāo)志層存在明顯錯(cuò)斷； 通過(guò)高精度匹配處理后，振幅能量趨于一致，標(biāo)志層錯(cuò)斷得到明顯消除，信噪比顯著提升。

圖12為混合震源高精度匹配處理前、后沿侏羅系大安寨地層的均方根振幅屬性平面圖?？梢?jiàn)在進(jìn)行高精度匹配處理后，可控震源與炸藥震源振幅能量關(guān)系一致性更好，為后續(xù)疊前保真處理奠定了資料基礎(chǔ)。圖13為混合震源高精度匹配處理前、后等時(shí)切片(1.4s)屬性平面圖，經(jīng)過(guò)高精度匹配處理，信噪比得到明顯提升，假的構(gòu)造形態(tài)得到了甄別，為精細(xì)構(gòu)造解釋和儲(chǔ)層預(yù)測(cè)提供了高保真數(shù)據(jù)。

圖13 混合震源高精度匹配處理前(a)、后(b)等時(shí)切片(1.4s)屬性平面圖

5 結(jié)束語(yǔ)

本文討論的混合震源高精度匹配處理方法，通過(guò)綜合應(yīng)用疊前平均振幅譜匹配、混合震源相位匹配及疊前剩余時(shí)差校正等處理技術(shù)，方便、快捷地實(shí)現(xiàn)了混合震源數(shù)據(jù)的一致性處理，較好地解決了不同震源激發(fā)(接收)造成的振幅、頻率、相位、時(shí)差等不一致問(wèn)題。該套技術(shù)具有較好的應(yīng)用前景，不但適用于陸上和海上混合震源采集資料，還可用于不同時(shí)間、不同采集因素施工獲得的多塊地震資料的連片處理，能顯著提高資料信噪比，淺、中、深層地震反射都能較好地拼接，實(shí)現(xiàn)了不同資料的同相疊加，并便于連續(xù)追蹤，為后續(xù)精細(xì)構(gòu)造解釋及儲(chǔ)層預(yù)測(cè)提供了高保真成果數(shù)據(jù)。