范寶倉 陳文貴 王海昆 李慶收

(中海油田服務股份有限公司物探事業(yè)部 天津 300451)

海底電纜Z分量橫波噪音特征分析及衰減方法*

范寶倉 陳文貴 王海昆 李慶收

(中海油田服務股份有限公司物探事業(yè)部 天津 300451)

在海底電纜(OBC)多分量地震資料的處理中，陸檢Z分量中除面波干擾外，還存在一種難消除的橫波側漏干擾噪音，這種橫波噪音具有矢量性、振幅強、頻帶窄、相速度低的特點，并且在十字交叉排列的共檢波點道集上表現(xiàn)為三維錐形特征。針對這種橫波噪音特征，結合地震采集施工特點，提出在十字交叉域共檢波點道集上的三維小波變換錐形濾波去噪方法。該方法通過小波變換，采用局部模擬噪聲，僅針對干繞波存在區(qū)域內進行，回避了傳統(tǒng)方法如F-K濾波等所帶來的有效信息丟失問題，而且進行自適應的減法運算消除干擾，增強了去噪的保真度，實現(xiàn)了真正意義上的針對性去噪處理。該方法在渤海某區(qū)塊實際資料的處理應用中取得了很好的噪音衰減效果，極大地提高了資料的信噪比。

海底電纜；Z分量；橫波噪音；特征分析；十字交叉排列；小波域；錐形濾波

目前海底電纜(OBC)多分量地震資料(水檢P分量，陸檢X、Y和Z分量)的處理流程已趨于成熟，其中用于壓制鬼波的水、陸檢合并(即P和Z分量合并)是必不可少的步驟[1]，但在渤海地區(qū)實際資料處理中，常常會發(fā)現(xiàn)陸檢資料(以下均指Z分量)的信噪比相對水檢而言普遍偏低。該地區(qū)陸檢資料干擾噪音中地滾波能量強、特征明顯，易于壓制，但陸檢單炮上還存在一種不規(guī)則的噪音(圖1)，在處理中很容易被當成隨機噪音而忽視它的影響。由于這種噪音在近道處與有效信號的能量和頻譜相近而很難區(qū)分，所以采用常規(guī)的壓制方法容易損傷有效信號；但如果為了保護有效信號而采取保守的去噪方式，往往又制約陸檢資料的信噪比提升，從而影響水、陸檢合并效果，并且很難再從后續(xù)的處理中解決所帶來的問題。針對上述問題，筆者通過分析該噪音的來源、產生機理、特征以及傳統(tǒng)去噪方法存在的問題，結合地震采集施工特點，提出在十字交叉域共檢波點道集上的三維小波變換錐形濾波方法去除噪音，并在渤海某區(qū)塊實際資料處理中取得了很好的噪音衰減效果，極大地提高了資料的信噪比。

圖1 渤海某測線地滾波壓制后的水、陸檢單炮Fig .1 SP gather after ground roll noise attenuation for a seismic line in Bohai sea

1 噪音特征分析

1.1 噪音來源

將圖1中的資料轉換為共檢波點道集(圖2)，陸檢所示的噪音則表現(xiàn)為規(guī)則的干擾，其近道處能量和頻譜與有效信號相近而很難區(qū)分，視速度明顯低于有效信號的視速度，而且對比發(fā)現(xiàn)該噪音形態(tài)特征與陸檢X、Y分量的資料極其相似。另外，從圖1、2中可以看出該噪音干擾只存在于陸檢，因此推測該噪音有可能是陸檢的X、Y分量的信息側漏帶來的，即陸檢橫波噪音。通過對該噪音的頻率成分進行分析(圖3)，發(fā)現(xiàn)噪音頻率主要集中于60 Hz以下，且表現(xiàn)為低速的線性或近雙曲的特征，其振幅和頻譜與有效波很相近，同樣表明該噪音為陸檢橫波噪音。

圖2 渤海某測線地滾波壓制后的水、陸檢共檢波點道集Fig .2 Common receiver gather after ground roll noise attenuation for a seismic line in Bohai sea

圖3 渤海某測線陸檢共檢波點道集頻率掃描Fig .3 Frequency scanning on common receiver gather for a seismic line in Bohai sea

1.2 噪音產生機理

海底電纜為四分量采集(圖4)，水檢P分量的檢波器為壓力檢波器，沒有方向性；陸檢為速度檢波器，具有方向性，包括Z、X、Y分量，其中X、Y分量平行于地面接收轉換波，而Z分量垂直于地面接收縱波資料[2]，且每個四分量檢波器單元都封裝在一個具有一定質量的扁平金屬外殼中[3]。

實際采集中，由于海底并非完全水平等原因，檢波器通常會與海底耦合不理想，即檢波器與海底往往存在一定的角度，這樣陸檢Z分量就有可能接收到轉換波資料，也可以直接理解為側漏的橫波信息，所以在陸檢Z分量上就會存在橫波干擾噪音。另外，單炮資料是由不同檢波器接收的，由于檢波器之間存在較大耦合差異，導致噪音干擾表現(xiàn)不明顯，而共檢波點道集是來自于同一個檢波器，耦合情況基本相同，所以在共檢波點上表現(xiàn)為比較規(guī)律的特征。

1.3 噪音在三維空間上的特征

在地震資料處理中，對于同一位置的檢波點，由所有與之相關的炮點記錄組成并按照炮點位置依次正交排列的道集稱為一個共檢波點十字交叉正交排列道集(圖5a)，因此可以通過選取一組十字交叉排列數(shù)據(jù)來進一步了解該噪音在三維空間上的特征。從十字交叉排列數(shù)據(jù)的時間切片(圖5b)上可以看到，橫波干擾在近偏移距處比較強，在每個切片上呈圓形，由淺到深圓不斷變大，形成三維矢量錐形體特征。而且結合圖3的分析可知，該噪音還具有類似于面波的低相速度特征，且主要集中在某一頻段內。

圖4 海底電纜四分量采集示意圖Fig .4 Diagram of OBC 4C receiver gather

圖5 十字交叉排列道集示意圖及時間切片F(xiàn)ig .5 Diagram of regular cross-spreading system and time slices

2 噪音衰減方法

2.1 常規(guī)去噪方法存在的不足

針對上述橫波噪音，顯然僅僅在共炮點道集上進行隨機噪音衰減是遠遠不夠的。而常規(guī)去除該噪音的方法是在共檢波點道集上利用二維拉冬域去噪或F-K濾波去噪。若去噪?yún)?shù)應用較保守，去噪后在共檢波點道集往往會有較多的殘余噪音，而且殘余噪音主要由于去噪方法的限制而集中在近偏移距附近(圖6a)；若去噪?yún)?shù)應用較激進，雖然道集上表現(xiàn)良好(圖6b),但疊加剖面上會發(fā)現(xiàn)明顯損失了部分有效波(圖6c)；因此，這2種常規(guī)方法的去噪效果都不理想，均會影響后續(xù)的疊加及水、陸檢合并的效果。

圖6 渤海某測線常規(guī)方法的噪音衰減效果Fig .6 Noise attenuation with conventional method for a seismic line in Bohai sea

2.2 十字交叉域小波變換三維錐形濾波法及去噪效果

在海底電纜多分量地震采集作業(yè)中，目前絕大部分采用較規(guī)則的密集施工方式，即炮線和檢波點線是正交或平行的，炮點和檢波點較密集。因此對于抽取共檢波點十字交叉排列道集相對較容易，且具有分布相對較均勻的特點，從而避免了因地震記錄的稀疏問題帶來的假頻影響[4-5]。分析可知，該橫波噪音在十字交叉排列三維空間上表現(xiàn)出明顯的規(guī)律特征，并能夠與有效波區(qū)別，因此可以采用在正交十字排列的共檢波點道集上進行小波域三維錐形濾波去噪處理，這主要是由于小波域處理較傳統(tǒng)的傅里葉變換具有一定優(yōu)勢，它是采用局部模擬噪聲，僅針對干繞波存在區(qū)域內進行，回避了F-K濾波、最小二乘濾波等傳統(tǒng)方法所帶來的有效信息丟失問題[6-9]。此外，小波變換可以通過伸縮平移運算對信號函數(shù)逐步進行多尺度細化，最終達到高頻處時間細分、低頻處頻率細分，即自動適應時頻信號分析的要求，從而可聚焦到信號的任意細節(jié)[5],更為高效保真地壓制噪音。

本文提出的十字交叉域小波變換三維錐形濾波法去噪處理過程包括5個步驟：①抽取十字交叉排列的共檢波點道集(圖5)；②根據(jù)頻率速度分析提取橫波干擾噪音的特征(圖7)，并設計相應的小波濾波器；③通過高精度小波變換將原始數(shù)據(jù)分解為多個小波子帶[6](圖8a)；④根據(jù)噪音的主要存在區(qū)域在相應的子帶內進行特征提取，建立噪音模型[7](圖8b、c)，并設計濾波器減掉噪音干擾[4、8-11]；⑤通過小波域數(shù)據(jù)反變換回時間域，完成去噪過程[4]。

圖7 原始共檢波點道集與對應的頻率-視速度Fig .7 Common receiver gather and frequency-apparent velocity panel

圖8 小波變換的子帶劃分及噪音模型Fig .8 Sub-band division and noise models in wavelet domain

渤海某區(qū)塊為海底電纜多分量地震采集區(qū)塊，其陸檢Z分量上存在明顯的橫波側漏干擾噪音(圖1、2)，圖9～11為該區(qū)塊某測線十字交叉排列的共檢波點小波變換三維錐形濾波方法應用后的道集質控圖，圖12為該測線十字交叉排列的共檢波點三維錐形濾波方法應用后的疊加質控圖。對比發(fā)現(xiàn)，無論是共炮點道集質控(圖9)，還是共檢波點道集質控(圖10)和三維質控(圖11)，橫波干擾噪音均得到了有效壓制；去噪前后的疊加剖面對比表明有效信號基本沒有損失，去噪效果較理想(圖12)。

圖13為該區(qū)塊某測線陸檢橫波噪音采用常規(guī)F-K濾波法和十字交叉排列三維錐形濾波法去噪后對水、陸檢合并的疊加效果對比及相應的自相關質控對比。很顯然，十字交叉排列三維錐形濾波法去噪效果更好，疊加剖面的信噪比更高，而自相關質控對比也說明了水陸檢合并的鬼波壓制效果更好，亦更有利于后續(xù)的水、陸檢合并處理。

圖9 本文新方法對渤海某測線橫波噪音衰減前后的共炮點道集Fig .9 Common shot gather before and after noise attenuation for a seismic line in Bohai sea with the suggested method

圖11 本文新方法對渤海某測線橫波噪音衰減前后三維立體顯示Fig .11 3D display before and after noise attenuation for a seismic line in Bohai sea with the suggested method

圖12 本文新方法對渤海某測線橫波噪音衰減前后疊加剖面Fig .12 Stack section before and after noise attenuation for a seismic line in Bohai sea with the suggested method

圖13 渤海某測線不同去噪方法的水陸檢合并效果對比Fig .13 Dual-sensor summation with different noise attenuation methods for a seismic line in Bohai sea

3 結論

1) 海底電纜陸檢Z分量中的橫波噪音主要是由于檢波器與海底的耦合不理想造成的橫波側漏，具有矢量性、振幅強、頻帶窄、相速度低的特點，在共炮點道集上表現(xiàn)為隨機噪音，而在十字交叉排列的共檢波點道集上表現(xiàn)為三維錐體特征。

2) 針對這種橫波噪音特征，結合地震采集施工特點，提出在十字交叉域共檢波點道集上的三維小波變換錐形濾波去噪方法。該方法通過小波變換，采用局部模擬噪聲，僅針對干繞波存在區(qū)域內進行，回避了傳統(tǒng)方法如F-K濾波、最小二乘濾波等所帶來的有效信息丟失問題，而且進行自適應的減法運算消除干擾，增強了去噪的保真度，實現(xiàn)了真正意義的針對性去噪處理。該方法在渤海地區(qū)的實際應用中取得了很好的噪音衰減效果，疊加剖面的信噪比更高，也更利于后續(xù)的水、陸檢合并處理。

[1] 周濱,龔旭東,高夢晗,等.海底電纜交叉鬼波化雙檢合并技術改進與應用[J].中國海上油氣,2015,27(1):49-52.DOI:10.11935/j.issn.1673-1506.2015.01.007. ZHOU Bin,GONG Xudong,GAO Menghan,et al.An improvement of the dual-sensor summation technique for cross-ghosting from OBC and its application[J].China Offshore Oil and Gas,2015,27(1):49-52.DOI:10.11935/j.issn.1673-1506.2015.01.007.

[2] 張文波,朱光明.海底電纜數(shù)據(jù)中壓力分量與垂直分量的分析與應用[J].地球科學與環(huán)境學報,2005,27(1):72-75. ZHANG Wenbo,ZHU Guangming.Analysis and application of pressure and vertical components of the ocean bottom cable data[J].Journal of Earth Sciences and Environment，2005,27(1):72-75.

[3] 王偉巍.海底地震勘探最新方法與技術發(fā)展[J].中國科技成果,2010,11(6):7-11. WANG Weiwei.The new method and technique progress about seismic exploration on seabed[J].China Science and Technology Achievements,2010,11(6):7-11.

[4] 李立紅.三維錐形濾波在稀疏十字交叉道集中的應用[J].天然氣技術與經濟，2011,5(5):25-29. LI Lihong.Application 3D cone filter on sparse cross gathers[J].Natural Gas Technology and Economy,2011,5(5):25-29.

[5] MORLET J,ARENS G,FOURGEAU E,et al.Wave propagation and sampling theory-Part II:sampling theory and complex waves[J].Geophysics,1982,47(2):222-236.

[6] 王劍平,張捷.小波變換在數(shù)字圖像處理中的應用[J].現(xiàn)代電子技術.2011,34(1):91-94. WANG Jianping,ZHANG Jie.Application of wavelet transform in digital image processing[J].Modern Electronics Technique.2011,34(1):91-94.

[7] 黃雪繼,劉來祥,王永勝.分頻徑向道中值濾波在地震資料處理中的應用[J].物探與化探,2012,36(2):317-320. HUANG Xueji,LIU Laixiang,WANG Yongsheng.The application of frequency-decomposed radial trace median filtering to seismic data processing[J].Geophysical & Geochemical Exploration,2012,36(2):317-320.

[8] 胡雙玉，徐春梅，蔣波,等.基于正交子集的疊前噪音壓制技術[J].科學技術與工程，2010,10(16):3832-3836. HU Shuangyu,XU Chunmei,JIANG Bo，et al.Pre-stack noise attenuation based on orthogonal subset[J].Science Technology and Engineering,2010,10(16):3832-3836.

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

[10] 吳長江.三維錐形濾波器壓制面波的應用效果分析[J].工程地球物理學報，2010,7(1):15-18. WU Changjiang.Application of 3-D cone-shaped filter to surface wave suppression and its analysis[J].Chinese Journal of Engineering Geophysis,2010,7(1):15-18.

[11] 龔旭東，張建峰，王志亮，等.8L4S176T觀測系統(tǒng)在渤海CF區(qū)塊海底電纜地震采集中的應用[J].中國海上油氣，2012,24(2):27-31. GONG Xudong,ZHANG Jianfeng,WANG Zhiliang,et al.An application of 8L4S176T observation system to seismic acquisition of ocean-bottom cable in CF block,Bohai sea[J].China Offshore Oil and Gas,2012,24(2):27-31. 收稿日期：2016-10-19 改回日期：2017-01-12

(編輯：張喜林)

Characteristic analysis and attenuation of OBC geophoneZcomponent shear wave noise

FAN Baocang CHEN Wengui WANG Haikun LI Qingshou

(COSLGeophysicalDivision,Tianjin300451,China)

Besides surface wave interference, the geophoneZcomponent noise from the shear wave leakage is difficult to be attenuated in OBC multi-components seismic data processing. This noise has the characteristics of vector, strong amplitude, narrow frequency, low phase velocity, and 3D cone-shaped in cross-spread receiver gather. Considering the characteristics of this noise and seismic acquisition, a new attenuation method of 3D cone-shaped filter using wavelet transform in cross-spread receiver gather is introduced. This method is based on wavelet transform and local analog noise, only used in the area which includesSwave interference, and can avoid the effective information loss caused by conventional methods such as F-K filtering. Furthermore the adaptive subtraction attenuation can enhance the fidelity of de-noising, realizing the true sense of the targeted de-noising. The data processing and application in a block of the Bohai sea show the effect of noise attenuation, which greatly improves the SNR of the data.

OBC;Zcomponent; shear wave interference; characteristic analysis; cross-spread; wavelet domain; cone-shaped filter

*國家高技術研究發(fā)展計劃(863計劃)“深水高精度地震勘探系統(tǒng)成套化研制(編號：2012AA09A211)”、中國海洋石油總公司“十二五”科技重大項目“海上開發(fā)地震技術集成及應用研究(編號：CNOOC-KJ125 ZDXM 06 LTD-10-KFSC-14)”部分研究成果。

范寶倉，男，高級工程師，1998年畢業(yè)于西南石油學院應用地球物理專業(yè)，獲學士學位，主要從事地震資料處理工作。地址：天津市濱海新區(qū)中心北路1889號514信箱(郵編：300451)。E-mail：fanbc@cosl.com.cn。

1673-1506(2017)03-0031-09

10.11935/j.issn.1673-1506.2017.03.005

P631.4

A

范寶倉,陳文貴,王海昆，等.海底電纜Z分量橫波噪音特征分析及衰減方法[J].中國海上油氣，2017,29(3)：31-39.

FAN Baocang,CHEN Wengui,WANG Haikun,et al.Characteristic analysis and attenuation of OBC geophoneZcomponent shear wave noise[J].China Offshore Oil and Gas,2017,29(3):31-39.