面向海洋剩余多次波衰減方法的研究

張如偉， 張寶金， 黃捍東， 文鵬飛

(1. 廣州海洋地質調查局 國土資源部海底礦產資源重點實驗室，廣州 510075；2.中國石油大學 油氣資源與探測國家重點實驗室，北京 102249)

?

面向海洋剩余多次波衰減方法的研究

張如偉1， 張寶金1， 黃捍東2， 文鵬飛1

(1. 廣州海洋地質調查局 國土資源部海底礦產資源重點實驗室，廣州 510075；2.中國石油大學 油氣資源與探測國家重點實驗室，北京 102249)

隨著近幾年長排列大容量震源地震勘探技術的廣泛采用，勘探目標逐漸趨向于中深部地層，使其成為目前研究南海深部構造的主要研究手段之一。但由于一些剩余多次波的影響，一直成為中深層成像與解釋的瓶頸。針對南海某海域的實際資料，在采用SRME技術與高精度Radon變換等常規衰減技術的基礎上，仍會存在一些殘余多次波，主要為崎嶇海底或者復雜構造引起的繞射多次波。針對剩余多次波的特點，在不同的頻帶范圍內，通過徑向中值濾波評估得到不同的“含信因子”用于異常噪音的衰減，同時為了提高數據保真度，引入逐步衰減的處理機制，層層檢驗衰減效果，最后結果表明,數據中剩余多次波能量被較好壓制，提高了中深層的信噪比，基底以下的深部構造得到了較好地展示。

多次波； 徑向中值濾波； 繞射信號； 保幅處理

0 引言

多次波一直是海洋地震處理的重點與難點，其存在著分布范圍廣、能量強、種類多等特點，是影響中深部地質構造成像的主要瓶頸[1]。從反射機理上來說，多次波有表層多次波、層間多次波和微屈多次波；從時差上來說，有長周期多次波與短周期多次波，海上地震資料還存在鳴震現象[2]。實際資料中一般是多種類型的多次波與一次波交叉在一起，難以識別和消除，如果在成像過程中處理不當，還會給后續地震資料的解釋帶來一系列問題[3]。

針對多次波的壓制方法，國內外眾多專家學者均進行了研究，Verschuur[4]與Weglein[5]基本將去除多次波的方法分為兩類：①為基于一次波與多次波在空間上的特性差異的去噪方法；②為基于多次波的周期性與可預測性的去噪方法，只是前者將預測反褶積劃為第二類，而后者將其劃為第一類；葉月明等[6]依據數據驅動對層間多次波的壓制進行研究，將多次波分解為三個波場分量，此三個波場通過互相關與褶積即可預測出層間多次波；李振春等[7]將波動理論與濾波理論向結合，基于角度域共成像道集實現多次波壓制；Dedem等[8]提出采用稀疏反演方法來預測表面相關多次波，Ketil等[9]在之后提出拋物線稀疏反演技術；王學軍等[10]針對阿曼Amal區塊地震數據中復雜多次波進行分析與壓制；郭夢秋等[11]分析了海上地震資料處理中的組合壓制多次波技術；謝宋雷等[12]對深水陡坡帶繞射多次波壓制的方法進行了研究；石穎等[13]提出基于波動方程預測和雙曲Radon變換聯合壓制表面多次波的方法；Mariusz等[14]針對海洋拖纜數據，采用上/下行波反褶積方法來壓制只有表面多次波，Joost等[15]在之后發展為多維度反褶積方法；劉華鋒等[16]利用改進的擴展多道匹配相減法壓制多次波；謝俊法等[17]將迭代拋物Radon變換法用于分離一次波與多次波。

作者在常規衰減多次波(SRME與高精度Radon變換)之后，分析了殘余多次波的特征，在地震數據分頻的基礎上，采用疊前道集的徑向中值濾波計算得到“含信因子”，并以此為基礎來衰減殘余異常噪聲。同時為了提高數據的保真度，引入逐步衰減的處理機制，每一步被衰減的噪音均須重復之前的處理步驟，最后全部回加到有效信號中。

1 剩余多次波特點

在海洋地震數據中，崎嶇海底與復雜構造會引起大量的繞射波，從而導致繞射多次也更加嚴重。由于繞射多次波一般不同于常規多次波的幾何特征，其不符合雙曲線頂點在近偏移距的規律，而是隨機分布，如圖1中的繞射多次波的反射頂點均位于遠偏移距，對于這樣的繞射多次能量，在采用常規的SRME技術與高精度Radon變換等常規處理流程[18]之后，剖面上還會殘留著較多的殘余多次波能量，如果不壓制而直接進行偏移處理，會出現明顯畫弧現象，極大程度上影響了中深層的成像效果。圖2為采用常規處理之后，疊加剖面上還會存在明顯的殘余多次波能量(CMP#7500-8500)，嚴重降低了數據整體的信噪比。

圖1 炮道集上的繞射多次波Fig.1 Diffraction multiples in the shot gather

雖然繞射多次波不同于常規多次波的顯著特征，但也比較容易設別，如圖1的6.3 s左右(黃框)與6.8 s左右(藍框)，其幾何特征類似于線性干擾，只是傾角不一致而已。其實大部分繞射多次波在局部基本以似線性的特征存在，而傾角卻又不同，作者只希望衰減某一些傾角的信號，將此作為噪音，采用徑向中值濾波即可較好地實現信噪分離，為其后的“含信因子”評估打下基礎。

圖2 常規流程壓制之后的剩余多次能量Fig.2 Residual multiples after convention processing flow(a)原始數據;(b)SRME處理后;(C)高精度Radon變換處理后

2 剩余繞射多次波衰減思路

圖3 面向剩余繞射多次波的二次衰減流程圖Fig.3 The flow of twice attenuation for residual diffraction multiples

基于逐步衰減的處理思路，在原始信號分頻的基礎上，通過徑向中值濾波處理，計算得到含信因子，并將含信因子作為權值來衡量衰減強度，在保持有效信號的前提下，最大可能地去除影響有效信號的干擾。為了達到保真的效果，可以進行多次衰減，最后將所有的數據相加，即為最終多次波衰減結果(圖3)。

假設疊前地震數據為X，通過分頻處理之后，不同頻段的疊前數據分別為X1、X2、X3，各自進行徑向中值濾波之后，則有信號道S1、S2、S3與噪音道

N1、N2、N3。該方法是一種非線性濾波技術，具有信噪分離效果好，保幅性較高，不降低分辨率等優勢[19]，主要利用了繞射多次波局部道線性化的特點，處理過程中主要包含兩個步驟：①局部徑向道求取；②局部徑向道噪音提取[20]。

含信因子W為有效信號在原始數據中的比重，計算公式如式(1)所示：

(1)

其中：S為信號道；N為噪音道；i為縱向采樣點；j為道數；m與k分別是計算窗口的橫向與縱向長度。

經過衰減之后的數據Y為：

Y= (X1-N1*W1)+(X2-N2*W2)+

(X3-N3*W3)

(2)

式中：W1、W2、W3分別為不同頻率段的含信因子，計算方式采用公式(1)。

同時可以對噪音道繼續進行分離，再計算含信因子，然后在對疊前數據進行二次衰減，在實際運算過程中，需要依據資料的情況進行分析衰減次數，但一般衰減2次～3次即可以到達保幅衰減效果。

3 衰減效果分析

這里選取了南海某海域的實際數據進行分析處理，該區域的地質目標主要為中深層地質特征研究， 而中深層往往是多種類型多次波發育的高發區域，圖4為研究區域發育的一些復雜的繞射多次波，基本很難發現有效信號的蹤跡，對繞射多次波的保幅衰減帶來一定的難度。

圖5為采用作者方法衰減殘余多次波之后的效果剖面對比。由圖5可以看出，在采用SRME與高精度Radon變換衰減(圖5(a))之后，剖面上還存在著一些較強的剩余多次波能量，使得中深層的信噪比比較低，如果不能夠衰減，將大大影響成像的質量，采用本文方法可以看出一些剩余多次被較好的衰減，剖面的信噪比得以明顯提高，8.5 s左右的有效信號更加清晰地展示出來。

圖6為應用作者方法衰減殘余繞射多次波前后偏移剖面效果對比圖。由圖6可以看出，如果不對這些殘余繞射多次波進行壓制，勢必會在偏移剖面上產生大量的畫弧現象，影響了有效信號的成像。而采用作者方法壓制剩余繞射多次波之后，偏移剖面上畫弧明顯減少，一些有效的中深地層反射特征展示出來(8.8 s～9.2 s)，利于地質人員的追蹤與解釋。

4 結論與認識

1)多次波發育嚴重是影響南海深部構造成像的主要障礙，采用常規流程往往會剩余一些多次能量，如果處理不當，將會大大降低剖面成像質量，并會帶來解釋的誤區。

2)雖然剩余的繞射多次波的幾何特征不同于常規多次波，但在到道集上局部呈現似線性化特征，并在傾角上與有效信號不同，可以采用徑向中值濾波進行信噪分離。

3)“含信因子”的計算使噪音衰減更加合理化，而不是簡單的相減，更加利于后期的保幅處理。

4)為了有效地保證有效信號的保真性，引入逐步衰減的處理機制，每一步被衰減之后的噪音，再重復之前處理流程，挖掘其中遺留的有效能量，為中深層的保幅處理提供了基礎。

[1] 孫福利, 王真理, 郝天珧, 等. 南海南部深部結構的復雜構造地震成像[J]. 地球物理學報, 2011,54(12): 3210-3216.

SUN F L, WANG Z L, HAO T Y, et al. Seismic imaging of complicated deep structures in southern South China Sea[J]. Chinese Journal Geophysics, 2011, 54(12): 3210-3216.(In Chinese)

[2] 張軍華，王要森，鄭旭剛，等.海上地震資料多次波特征分析[J].石油地球物理勘探, 2009,44(5):574-577.

ZHANG J H, WANG Y S, ZENG X G, et al. Offshore seismic data multiple characteristics analysis[J]. Oil Geophysics Prospecting, 2009, 44(5):574-577.(In Chinese)

[3] YUAN S Y, WANG S X, LUO C M,et al. Simultaneous multitrace impedance inversion with transform-domain sparsity promotion[J]. Geophysics, 2015, 80(2): R71-R80.

[4] VERCHUUR.D.J, A.J.BERKOUT, C.P.A.WAPENAAR. Adaptive surface-related multiple elimination[J].Geophysics,1992,57:1166-1177.

[5] WEGLEIN.A.B, Multiple attenuation: an overview of recent advances and the road ahead[J].The Leading Edge, 1999:18(1):40-44.

[6] 葉月明,趙昌壘,姚根順,等.數據驅動型層間多次波預測方法研究[J].石油地球物理勘探, 2014,49(2):244-251.

YE Y M, ZHAO C L, YAO G S,et al. Study of data-driven interbed multiple prediction[J]. Oil Geophysics Prospecting, 2014, 49(2):244-251. (In Chinese)

[7] 李振春,李志娜,郭書娟,等.基于成像域的多次波壓制[J].石油地球物理勘探, 2013,48(6):903-909.

LI Z C, LI Z N, GUO S J, et al. Multiple attenuation based on image space[J]. Oil Geophysics Prospecting, 2013, 48(6):903-909.(In Chinese)

[8] DEDEM E. J, VERSCHUUR D. J. 3D surface-related multiple prediction: A sparse inversion approach[J].Geophysics, 2005, 70(3):31-43.

[9] KETIL H, ROGER S. 3D surface-related multiple elimination using parabolic sparse inversion[J].Geophysics, 2006, 71(6):145-152.

[10]王學軍,李虹,郭惠英,等.阿曼Amal區塊地震數據中復雜多次波的分析與壓制[J].石油地球物理勘探, 2012,47(6):894-900.

WANG X J, LI H, GUO H Y, et al. Analysis and suppression of complex multiples for land seismic data in Block Amal, Oman[J] . Oil Geophysics Prospecting, 2012, 47(6):894-900.(In Chinese)

[11]郭夢秋,趙彥良,左勝杰,等.海上地震資料處理中的組合壓制多次波技術[J].石油地球物理勘探, 2012,47(4):537-544.

GUO M Q, ZHAO Y L, ZUO S J, et al. Combined multiple attenuation in marine seismic data processing[J] . Oil Geophysics Prospecting, 2012, 47(4):537-544.(In Chinese)

[12]謝宋雷,劉伊克,王一博,等. 深水陡坡帶繞射多次波壓制方法研究[J].地球物理學報, 2013,56(9): 3118-3123.

XIE S L, LIU Y K, WANG Y B,et al. Study of diffracted multiple elimination on steep slope zone in deep water[J]. Chinese Journal Geophysics, 2013,56(9): 3118-3123.(In Chinese)

[13]石穎, 王維紅. 基于波動方程預測和雙曲Radon變換聯合壓制表面多次波[J]. 地球物理學報, 2012,55(9): 3115-3125.

SHI Y, WANG W H. Surface-related multiple suppression approach by combining wave equation prediction and hyperbolic radon transform[J]. Chinese Journal Geophysics, 2012, 55(09): 3115-3125.(In Chinese)

[14]MARIUSZ M, CLEMENT K, KRAGH E,et al. Attenuation of free-surface multiples by up/down deconvolution for marine towed-streamer data[J]. Geophysics, 2011, 76(6):129-138.

[15]JOOST NEUT, MARTIJN F, ROALD B. Data matching for free-surface multiple attenuation by multidimensional deconvolution[J]. Geophysical Journal International, 2012,191(2):743-750.

[16]劉華鋒,李慶春,王立明.利用改進的擴展多道匹配相減法壓制多次波[J].石油地球物理勘探, 2009, 44(3):270-275.

LIU H F, LI Q C, WANG L M. Suppression of multiples by improved expanding multi-channel match minus method[J]. Oil Geophysics Prospecting, 2009, 44(3):270-275.(In Chinese)

[17]謝俊法, 孫成禹, 韓文功.迭代拋物Radon變換法分離一次波與多次波[J].石油地球物理勘探, 2014,49(1):76-81.

XIE J F, SUN C Y, HAN W G. Iterative parabolic radon transform for primary and multiple separations [J]. Oil Geophysics Prospecting, 2014, 49(1):76-81. (In Chinese)

[18]張如偉，張寶金，文鵬飛，等.南海某工區復雜多次波衰減與中深層成像[C]. 中國地球物理學會第二十八屆年會,2012:748-749

ZHANG R W, ZHANG B J, WEN P F,et al. Complex multiple attenuation and middle-deep imaging in South China Sea[C]. The 28thannual meeting of Chinese geophysics, 2012:748-749.(In Chinese)

[19]ZHU WEIHONG, PANOS G. Linear noise attenuation using local radial trace median filtering[J]. The Leading Edge, 2004, 23(8): 728-730.

[20]黃雪繼，劉來祥，王永勝. 分頻徑向道中值濾波在地震資料處理中的應用[J]. 物探與化探，2012,36(2)：317-320.

HUANG X J, LIU L X, WANG Y S. The application of frequency-decomposed radial trace median filtering to seismic data processing[J]. Geophysical and Geochemical Exploration, 2012, 36(2): 317-320.(In Chinese)

The method of marine residual multiple attenuation

ZHANG Ru-wei1, ZHANG Bao-jin1, HUANG Han-dong2, WEN Pen-fei1

(1.Key Laboratory of Marine Mineral Resources, Ministry of Land and Resources , Guangzhou marine geological survey,Guangzhou 510760,China; 2.State Key Laboratory of Petroleum Resource and Prospecting,China University of Petroleum,Beijing 102249,China)

Because of using expanding spread and big source seismic acquisition, exploration targets gradually tend towards mid-depth strata, it will be an important research method for deep seated structure in South China Sea. But due to much interferences of residual multiple, that always disturb imagination and interpretation of mid-depth strata. In this paper, we apply actual data in south china sea, and the data always contain more residual multiples after using SRME technique and high accuracy Radon transform, which are diffractions of rough seafloor and complex construct. Based on the features of residual multiple, in the different extent of frequency band, this method will apply radial median filtering in calculating “signal-bearing factor” to attenuate the anomaly noise. Moreover, in order to improve the quality of data, we import a multistep attenuation processing mechanism to examine the attenuation effect. The final results indicate that residual multiples have been suppressed, and have enhanced signal-to-noise ratio of mid-depth strata, the deep seated structures beyond the basement have been better displayed.

multiple; radial median filtering; diffraction; amplitude-preserved processing

2015-06-30 改回日期：2015-07-22

國家自然科學基金(41176056)

張如偉(1984-)，男，碩士，主要從事于海洋地震資料處理與天然氣水合物預測研究工作，E-mail：cgszrw@163.com。

1001-1749(2016)05-0666-06

P 631.4

A

10.3969/j.issn.1001-1749.2016.05.15