高密度寬方位地震資料處理技術研究進展

王兆磊， 公 亭， 李隆梅， 李國生， 趙志強, 顧小弟， 孟曉夢

(東方地球物理公司 研究院處理中心，涿州 072751)

?

高密度寬方位地震資料處理技術研究進展

王兆磊， 公 亭， 李隆梅， 李國生， 趙志強, 顧小弟， 孟曉夢

(東方地球物理公司 研究院處理中心，涿州 072751)

隨著可控震源滑動掃描、獨立同步掃描、高保真采集、點激發接收等高效采集技術的推出，高密度、寬方位地震勘探變得越來越經濟可行，但隨之也帶來了高效采集噪聲干擾、高保真采集數據分離等問題。這里著重論述了相關預測減去法諧波壓制、多域壓制鄰炮干擾、線性反演拆分混疊數據等方法，同時展示了擬全三維和寬方位數據的處理成果。通過運用以上技術取得了較好的信噪分離效果，體現了寬方位數據在解決成像和方位各向異性問題上的優勢。

高密度； 寬方位； 方位各向異性； 方位角道集

0 引言

近年來隨著可控震源滑動掃描(Slip-Sweep)、獨立同步掃描ISS(Independent Simultaneous Sweep)、高保真采集HFVS(High Fidelity Vibrator Seismic)、點激發接收等高效采集技術的推出，加之萬道以上地震儀、海量數據存儲、GPU等硬件設備的出現，地震勘探已經從窄方位常規采集，發展到寬方位或高密度采集，并且正在向高密度加全方位采集方向發展[1-3]。

新的高效采集技術在更經濟地獲得高密度寬方位地震資料的同時，也給地震數據處理帶來了新的問題，如可控震源的諧波壓制問題、ISS的鄰炮干擾問題、HFVS的數據分離問題等。而學者們也對地震數據處理提出了更高的要求，這些需求包括：更高的信噪比(如淺層薄層的連續性)、更高的分辨率(識別3 m～5 m斷層、5 m～10 m儲層)、更精確的成像(如碳酸鹽巖溶洞)、更豐富的描述方位各向異性信息(如裂縫)等。諸多的問題和需求，都促進了地震資料處理技術的發展。

1 高效采集資料相應處理技術

1.1 可控震源諧波壓制

由于可控震源的機械裝置、震動裝置和液壓伺服系統的非線性，以及震板與大地的耦合效應，使可控震源在向地下傳輸掃描信號的同時，不可避免地產生諧波干擾。而當使用滑動掃描升頻方式施工時，后一炮的諧波干擾會較強地污染到前一炮數據。針對這一問題國內、外學者進行了大量研究，并從各自角度提出多種壓制方法，主要包括：①相移濾波法，該方法通過設計一個純相位濾波器來壓制諧波，通常由于相關前數據，壓制諧波效果較好，但計算量較大；②相關預測減去法，該方法基于震源信號設計諧波預測算子，求取相關后記錄中的諧波干擾，然后將求得的諧波干擾從被干擾區域中剔除，具有諧波效果較好、計算量相對較小、實現簡便等特點；③分數傅里葉變換法，該方法借助分數傅里葉變換，使可控震源掃描信號在某個分數階上實現能量聚焦，在此基礎上設置濾波器壓制噪音，此方法實現較為復雜且容易引起數據畸變；④特征信號反褶積法，該方法利用震源力信號與相關前數據做反褶積，代替傳統的用掃描信號做相關，達到及壓縮震源子波又壓制諧波的目的，此方法實現原理簡單容易實現但對諧波壓制效果一般；⑤可控震源高保真采集數據分離技術，該方法采用反演代替相關，不但具有較高的保真度，同時可以利用諧波中的高頻信息，但需要采用特殊的施工方式[4- 5]。

這里研究了GeoEast軟件壓制諧波干擾所采用的相關預測減去法原理、步驟及其壓制效果[6]：

1)利用震源力信號，把基波H1和各次諧波Hi(i=2,3,4…)分離，并可以得到在頻率域的諧波預測算子：

(1)

其中：wi為權系數；ρi為振幅校正項；φi為相移項。

2) 經震源相關后，對初至后的有效信號S應用諧波預測算子P，可以求得初至前的諧波N。

N=S·P

(2)

3)將求得的諧波干擾從上一炮相應位置中減去，從而達到諧波干擾壓制的目的。

由于該方法可以針對可控震源相關后數據進行處理，因此具有計算效率高的優點。

圖1展示了該方法的壓制效果。其資料為東方地球物理公司2011年在吐哈盆地采集的國內第一塊三維滑動掃描地震數據，從圖1(a)所示的諧波干擾壓制前炮集記錄上，可以明顯地看到諧波干擾(綠色箭頭所示位置)，取受諧波影響較重的4101道做時頻譜，從中可以更清楚地看到2s以下數據受到較強的2階、3階、4階諧波干擾，在4s以下區域尤為嚴重。從圖1(b)諧波干擾壓制后炮集記錄及其時頻譜圖1(d)可以看出，采用諧波壓制技術后，各階諧波干擾得以明顯壓制。

從圖2諧波壓制前、后疊加剖面來看，諧波干擾壓制前疊加剖面2s以下區域受到較強諧波干擾，在頻譜上反映為40Hz以上高頻干擾。經諧波壓制處理后，從疊加剖面來看，諧波干擾基本消除，在相應的頻譜上也看不到諧波產生的高頻噪聲。

1.2 ISS鄰炮干擾壓制

對于采用ISS技術采集的地震資料，由于多組震源可以同時施工(不必等到上一炮記錄完成后，下一炮已經開始激發)，因此如圖3(a)所示，同一檢波點在同一時刻可以接收到來自不同炮點的信息，這種干擾稱為鄰炮干擾[7]。這種相鄰，不僅體現在兩炮的激發起始時間臨近(不大于記錄時間)，也體現在距離遠近，距離越近、時間間隔越短干擾越重。

圖1 諧波干擾壓制前后炮集記錄及其時頻譜分析

圖2 諧波干擾壓制前后疊加剖面及其頻譜分析

圖3 ISS鄰炮干擾壓制前后炮集記錄及減出噪聲

鄰炮干擾在野外采集時可以通過增加放炮時間間隔的方式減小影響，但這同時降低施工效率；也可以通過增大野外一次擺放排列的長度和寬度，使時間上相鄰的兩炮空間上距離增加，減小干擾影響范圍，但這需要有足夠多的設備投入，因此需在質量、效率、效益之間取得平衡。

對于室內處理，鄰炮干擾在炮域內成片出現具有很好的相干性，可以作為相干噪聲進行壓制，但對于不同炮來說，由于相鄰兩炮的放炮時間間隔、空間相對位置關系各不相同，很難通過統一參數設置達到較好壓制效果。在空間上每一炮都是一個獨立源，可以采用類似固定源的方式壓制。但對于高效采集數據來說，這種干擾的數目巨大，操作起來不現實。因此在實際處理時，通常采用將炮域數據變換到共炮檢距域、CMP域、OVT域或其他域中，使鄰炮干擾呈現出更好的隨機性，然后用中值濾波、異常振幅壓制等統計學方法進行壓制。圖3(b)即為壓制后再轉換回炮域的結果，從中可以看出，ISS鄰炮干擾得到較好壓制。

根據作者經驗，在壓制時應本著“寧輕勿重”的原則，避免傷及有效信號，因為對于高密度資料，即使在道集資料上殘留微弱噪聲，在后續的疊加或偏移過程中也會得到較好壓制，不會影響最終處理結果。

1.3 HFVS數據分離

在可控震源高保真采集(HFVS)施工過程中，同時使用多臺震源車多次激發，每臺震源在一指定的震源點位置或其周圍，并以不同于其他震源車的特有相位執行掃描，這樣就得到一個來自不同地面位置的不同震源車的多震源多次掃描數據集，再通過數據處理技術將其分離成不同震源的炮集記錄。

其分離原理是通過線性反演實現的。眾所周知，地震記錄可以表示為震源信號與大地脈沖響應的褶積，因此可控震源高保真采集所得多震源數據集可以表示為[8]：

(3)

其中：di(t)表示記錄到的第i次掃描的多組震源記錄；sij(t)表示第j臺震源第i次的震源掃描信號；ej(t)表示第j臺震源對應的大地脈沖響應。

因此在頻率域，利用N臺震源進行M次掃描的HFVS數據可以表示為以下面矩陣形式：

(4)

當M≥N即同一震源重復掃描次數大于等于震源臺數時，式(4)為一超定線性方程組，可以方便的求出其最小二乘解。

圖4展示了HFVS記錄分離效果，其中圖4(a)為四臺震源四次掃描的原始記錄，圖4(b)為分離后的結果。從圖4(b)中可以看出，分離后得到的數據，不但將四臺震源記錄分開，而且完成了以往通過互相關實現的可控震源信號壓縮工作。

2 不同地質需求的針對性處理技術

隨著勘探開發的逐步深入，學者們對地震數據處理提出了更高的要求。針對不同的地質需求，如何發揮高密度寬方位地震資料的優勢，是地震數據處理面臨的一大挑戰，通過兩個實例研究這方面的處理技術。

2.1 擬全三維數據融合處理技術

前人已經做過大量工作，很好地總結了寬方位地震數據的優點。諸如：①在地下介質復雜或未知的情況下，寬方位角采集比窄方位角更容易跨越地表障礙物和地下陰影帶，因此可以增加采集照明度，獲得較完整的地震波場；②炮檢對的三維疊前成像軌跡是橢球，因此寬方位角具有更高的陡傾角成像能力和更豐富的振幅成像信息，且寬方位角比窄方位角的成像分辨率更高；③由于寬方位角、窄方位角在炮點和檢波點的空間采樣特性不同，寬方位角成像的空間連續性優于窄方位角等。其實寬方位數據優于窄方位數據是不言而喻的，畢竟我們可以通過數據抽取將寬方位數據變成窄方位數據，但是現有的處理技術還不能將窄方位數據處理成可信的寬方位數據。

但是對于已有窄方位數據的地區，如何充分利用原有窄方位數據，在控制成本的前提下獲取高密度寬方位的地震資料？解決辦法之一就是采用擬全三維采集加融合處理技術。

如圖5所示某探區在原有12線30炮南北向采集數據基礎上，重新部署16線10炮東西向采集三維，兩者一并組成擬全三維數據。由于二次采集在設計時就應充分考慮到后續的融合處理，因此在激發接收方式上應盡量與上次保持一致，在觀察系統方面要保證面元一致、覆蓋次數基本相當。

圖4 HFVS記錄分離前后炮集數據

圖5 擬全三維觀測系統及玫瑰圖

經過多年的改革與發展，我國農村金融體制逐步完善，農村金融業取得了很大的發展成就。經歷了多次改革，我國逐步形成了由農業銀行、農村信用社、農業發展銀行、郵政儲蓄銀行、村鎮銀行和小額貸款公司等組成的農村金融體系。

以往的連片處理強調的是在橫向上“連成一片”，與此相區別擬全三維資料處理強調的將兩次采集數據在縱向上“融為一體”。因此，處理過程中每一個關鍵環節都要監控兩個融合數據體的一致性。

具體步驟為：①靜校正，首先應調查兩次采集是否存在系統時差，在消除系統時差的基礎上統一計算基準面靜校正以及剩余靜校正；②提高信噪比處理，詳細分析兩次采集各自的噪聲發育情況，若相同則采用統一的去噪流程，若不同則針對各自特點分別針對性壓制，確保兩次采集資料具有大致相當的信噪比；③振幅處理，采用地表一致性振幅補償消除因激發接收因素造成的兩次采集資料能量差異，在此基礎上分別提取兩次采集資料的沿目的層振幅屬性，如仍存在明顯差異，可采用倒數加權的方式加以消除；④提高分辨率處理，采用地表一致性反褶積消除因激發接收因素造成的兩次采集資料的子波差異，然后分別提取兩次采集資料的地震子波調查他們的一致性，必要時繼續做整形處理。

通過以上步驟最大程度地消除了兩次采集差異，將其兩次采集數據融合為一個整體。

從圖6的偏移結果看，擬全三維數據通過新采集補充了原有數據東西向遠偏移距信息，解決了原來困擾處理人員聯絡測線方向(東西向)偏移不足的問題。

2.2 高密度寬方位資料分方位處理技術

我們都知道可以通過增加道密度的方式提高淺層信噪比與有效分辨率，并且高密度采集可以在保證覆蓋次數的前提下采用更小的面元，這樣可以提高空間分辨率。但是高密度數據的另一個潛在的優勢是可以在保證有效覆蓋次數與偏移距分布均勻的前提下，將寬方位數據劃分成更多個方位角道集，而這對于方向各向異性介質來說，可更好地研究振幅隨炮檢距和方位角的變化(AVOA)，增強識別斷層、裂隙和地層巖性變化的能力。

現有的高密度寬方位資料分方位處理技術主要包括[9-10]：①各向異性偏移技術，用于消除各個方向上由各向異性引起的遠偏移距動校誤差；②分方位速度分析加多方位網格層析技術，常規的網格層析沒有考慮到剩余時差與方位角的關系，多方位網格層析，在分方位拾取剩余時差的基礎上，統一分解計算，得到更為精細和考慮了方位各向異性的速度場；③方位時差校正技術，當不利用方位時差信息時，為了提高全方位資料的同相疊加程度，可以采用該方法消除方位各向異性引起的方位時差；④OVT域偏移技術，由于OVT域偏移所得道集可以保存方位角信息，避免了先分方位再偏移的麻煩，有利于多次分方位處理。

圖6 老數據過井十字剖面與擬全三維過井十字剖面

圖7 高密度全方位地震資料分方位處理時間切片

圖8高密度全方位地震資料疊后方位角道集

圖7為某探區高密度全方位資料，從0°到165°每隔15°分12個方位的疊前時間偏移3 s時間切片(每個方位角跨度為45°)。從圖7中可以明顯看出，方位角為30°的時間切片最能描述工區內主測線方向的裂縫，而與之垂直的方位角為120°的時間切片上，幾乎看不到裂縫現象，其他各方向時間切片描述該裂縫的能力在兩者之間漸變。

將以上12個方位的疊前時間偏移疊加數據，在同一面元內按方位角排序，組成疊后方位角道集。圖8為抽取的斷裂位置處道集，每一小段為一個CMP，其中每一道即為一個方位角(一個CMP內有12道，即為12個方位角，從零度開始間隔15°一個)。從圖8中同樣可以看到，在3 s左右的位置，30°的疊后地震道較其相鄰道有明顯突變，通過以上手段可以更好地提取高密度、寬方位地震資料的方位各向異性及裂縫信息。

3 認識

1)可控震源高效采集地震資料針對性處理技術，可以較好地解決由于高效采集給高密度、寬方位地震資料帶來的可控震源諧波干擾、ISS鄰炮干擾、HFVS的數據分離等問題。

2)通過針對性的處理技術，可以使高密度、寬方位地震資料具有更高的信噪比、分辨率、成像精度以及更豐富的描述方位各向異性信息。

3)對于尚無條件采用高效采集的地區，可以在原有的窄方位資料基礎上，通過擬全三維采集技術達到提高道密度、填補方位角的目的，并且通過融合處理技術可以實現更好的成像效果。

4)對于寬方位處理應盡可能消除方位各向異性，達到同相疊加；對于分方位處理應凸顯其方位各向異性信息。

[1] 凌云,高軍,孫德勝.可控震源在地震勘探中的應用前景與問題分析[J].石油物探,2008,47(5):425-438. LING Y, GAO J, SUN D SH. Analysis of vibroseis in seismic exploration and its application. Geophysical Prospecting for Petroleum,2008,47(5):425-438.(In Chinese)

[2] 劉振武,撒利明,董世泰,等.中國石油高密度地震技術的實踐與未來[J].石油勘探與開發,2009,36(2):129-135. LIU ZH W,SA L M,DONG SH T. Practices and expectation of high-density seismic exploration technology in CNPC. Petroleum Exploration and Development,2009,36(2):129-135.(In Chinese)

[3] 凌云研究小組. 寬方位角地震勘探應用研究[J].石油地球物理勘探.2003,38(4): 350-357. LING Y study group. Application and study on wide-azimuth seismic exploration . Oil Geophysical Prospecting,2003,38(4): 350-357.(In Chinese)

[4] 曹務祥,張慕剛.滑動掃描諧波分析[J].石油地球物理勘探，2005,40(5):499-510. CAO W X,ZHANG M G. Slide-sweeping harmonic analysis. Oil Geophysical Prospecting,2005,40(5):499-510.(In Chinese)

[5] 魏國偉,張慕剛,魏鐵.可控震源滑動掃描采集方法及應用[J].石油地球物理勘探,2008,43(2):67-70. WEI G W, ZHANG M G, WEI T. Slide-sweeping acquisition method of vibroseis and application. Oil Geophysical Prospecting,2008,43(2):67-70.(In Chinese)

[6] 王寶彬,李合群,趙波.滑動掃描地震數據的諧波干擾壓制[J]. 石油地球物理勘探,2013,48(2):37-41. WANG B B,LI H Q,ZHAO B. Harmonic noise removal on vibroseis slip-sweep data . Oil Geophysical Prospecting,2013,48(2):37-41.(In Chinese)

[7] 魏鐵,于世東,于敏杰.可控震源噪聲分析[J].石油地球物理勘探,2008,43(2):38-44. WEI T, YU SH D, YU M J. Analysis of noise in vibroseis. Oil Geophysical Prospecting,2008,43(2):38-44.(In Chinese)

[8] 張留爭,趙志寬,楊曉昭,等. 采用428XL儀器實現高保真可控震源三維地震采集[J]. 物探裝備,2007,17(4):245-250. ZHANG L ZH,ZHAO ZH K,YANG X ZH, et al. Using 428XL seismograph to implement 3-D hi-fi vibroseis seismic acquisition. Equipment for Geophysical Prospecting,2007,17(4):245-250.(In Chinese)

[9] 張保慶,周輝,左黃金,等.寬方位地震資料處理技術及應用效果[J].石油地球物理勘探,2011,46(3):396-400. ZHANG B Q,ZHOU H,ZUO H J,et al. Wide azimuth data processing techniques and their applications . Oil Geophysical Prospecting,2011,46(3):396-400.(In Chiness)

[10]張軍華,朱煥,鄭旭剛,等.寬方位角地震勘探技術評述[J].石油地球物理勘探,2007, 42( 5):603-609. ZHANG J H, ZHU H, ZHENG X G,et al. Summary of wide-azimuth seismic exploration technique . Oil Geophysical Prospecting,2011,46(3):396-400.(In Chinese)

The new progress of high-density and width azimuth seismic data processing

WANG ZHao-lei, GONG Ting, LI Long-mei, LI Guo-sheng, ZHAO Zhi-qiang,GU Xiao-di, MENG Xiao-meng

(Geophysical Research Institute, BGP INC., China National Petroleum Corporation,Zhuozou 072751,China)

With the development of high-productivity vibrate, techniques such as slip-sweep, independent simultaneous sweep, high fidelity vibrator seismic acquisition method, seismic exploration with width azimuth and high-density data has become viable economically. This paper introduces key processing techniques of these vibrate is acquisition methods, which are harmonic attenuation for slip-sweep, interfered shots de-blend for ISS, and separation of high fidelity acquisition data. This paper also discusses targeted processing techniques according to different geologic demands and gives results.

high-density; width-azimuth; azimuth-anisotropy; azimuth-angle-gathers

2014-05-30 改回日期：2014-10-08

中國石油天然氣股份有限公司重大科技專項(2012E-34-11)

王兆磊(1978-)，男，博士，高級工程師，從事地震數據處理及處理方法研究工作，E-mail:wzrmxm@126.com。

1001-1749(2015)04-0465-07

P 631.4

A

10.3969/j.issn.1001-1749.2015.04.09