基于CRP道集的頻寬一致性子波反褶積技術

明　君,劉力輝,丁　燕,李　聰

（1.中海石油（中國）有限公司天津分公司，天津300452；2.北京諾克斯達石油科技有限公司，北京100192）

技術方法

基于CRP道集的頻寬一致性子波反褶積技術

明君1,劉力輝2,丁燕2,李聰2

（1.中海石油（中國）有限公司天津分公司，天津300452；2.北京諾克斯達石油科技有限公司，北京100192）

將疊后復賽譜域子波提取技術用于疊前CRP道集子波提取和疊前提高分辨率處理時，由于道集內遠近偏移距道頻寬不一致，影響了子波的空間一致性，也影響了疊前子波反褶積提高分辨率的效果。為此，對復賽譜域子波提取技術進行了改進，使其提取的子波具有空間一致性，并將改進后的復賽譜域子波提取技術與道集譜平衡技術相結合，形成了一種道集內頻寬一致性子波反褶積技術。該技術增加了子波反褶積希望輸出的人為定義，使其能適應低信噪比資料的疊前提高分辨率處理。

復賽譜域子波提取；譜平衡；頻寬一致的子波反褶積；疊前提高分辨率

0　引言

CRP道集優化是疊前反演和疊前屬性提取的基礎［1］，因此，基于CRP道集的優化處理越來越受到人們的重視。CRP道集優化主要包括提高信噪比、提高分辨率及增強保真度等方面的處理［2］，其中，提高分辨率處理是關鍵。AVO屬性分析的基礎是佐布里茲方程，而該方程是基于單一界面的，沒有考慮層間干涉，所以在薄層地區CRP道集優化中的提高分辨率處理對AVO屬性分析至關重要［3］。此外，CRP道集優化中的提高分辨率處理若能拓寬有效頻帶，對疊后進一步提高分辨率也具有重要意義，但疊前道集一般信噪比不高，是疊前提高分辨率處理的一個難點。

確定性子波反褶積、時變譜白化和反Q濾波是常規疊后提高分辨率處理的三大技術。對于疊前CRP道集上的提高分辨率處理，雖然時變譜白化和反Q濾波技術都具有一定的效果，但時變譜白化技術要求資料有一定的信噪比，而反Q濾波技術也受到多種因素的影響，如穩定性、噪聲和Q值誤差等。因此，針對低信噪比資料，在實際應用中尋找一種簡便且穩定的疊前提高分辨率處理技術，具有很重要的現實意義。確定性子波反褶積中的同態反褶積［4］是一種非線性變換，對子波的相位沒有要求，它通過將地震頻譜轉換成復賽譜，可有效地分離出地震子波和反射系數。該方法抗噪性相對較強，應用較廣，但將它用于疊前提高分辨率處理時往往效果不甚理想。這主要是由于其假設子波不隨時間和空間發生變化［5-7］，而實際疊前CRP道集資料中遠近偏移距道的子波頻帶寬度、振幅能量及相位均會隨著偏移距的變化而出現較大的改變［5，8］，顯示出明顯的空間不一致性，導致使用傳統的同態反褶積方法達不到理想的提高分辨率效果。為此，筆者對復賽譜域混合相位子波提取的同態反褶積方法進行改進，并將它與CRP道集的譜平衡技術相結合，形成一種新的頻寬一致性子波反褶積技術，用于疊前道集的提高分辨率處理。

1　方法原理

1.1復賽譜域子波提取技術

忽略噪聲項，地震記錄褶積模型在頻率域表示為

式中：x（ω），w（ω）和ξ（ω）分別為地震記錄、子波和反射系數的傅氏變換。

通過對式（1）取對數，可以將其轉化為線性系統，即

對式（2）作反傅氏變換，則有

式中：x（t），w（t）和ξ（ω）分別為地震記錄、子波和反射系數的時間序列的復賽譜序列。

由于子波和反射系數序列“平滑度”的差別，子波的復賽譜一般在原點附近，而反射系數序列的復賽譜遠離原點［9-11］。基于這種認識，在復賽譜域設計低通濾波器就可以實現子波與反射系數的分離，達到子波提取的目的。相同振幅譜而相位譜不同的地震子波對應不同的最大或最小相位分量。在復賽譜域，遵照一定的規則，可將子波振幅譜分割為2部分，分別對應于子波的最大和最小相位分量，分割比例的不同決定了子波的不同相位，由此可以得到最小相位子波、零相位子波和一系列混合相位子波。根據地震記錄信噪比譜，設定地震記錄的希望輸出，對地震記錄進行確定性子波反褶積，再依據反褶積后地震記錄方差模的大小和分辨率情況，最終從一系列相同振幅譜而相位譜不同的地震子波中確定理想的混合相位子波。

1.2頻寬一致性子波反褶積技術

頻寬一致性子波反褶積技術主要是將改進后的復賽譜域子波提取技術與道集譜平衡技術進行了結合，其方法原理為：首先進行道集內子波提取及道集處理，在CRP道集內取標準道，采用Gabor變換的時頻譜對遠偏移距道進行譜平衡處理，并提取子波；其次，在道集間選取多個提取的子波形成子波集，采用“優化平均”法得到一個空間一致性子波；最后，根據井-震匹配和預測目的，人為定義希望輸出子波的振幅譜，并求取反子波，從而獲得分辨率得到了提高的疊前數據。

1.2.1疊前道集譜平衡技術

在實際應用中，采用同態反褶積方法對有些資料進行提高分辨率處理的效果好，而有些效果差。仔細分析后發現，提高分辨率效果差的資料大多都存在道集內遠近偏移距道頻寬不一致的現象，即遠偏移距道波形“發胖”［12］，這會造成統計子波振幅譜時空間的不一致性。對有些資料作道集內遠近偏移距道的頻譜分析時會發現，遠偏移距道的有效頻帶比近偏移距道的有效頻帶窄，即遠近偏移距道的頻寬不一致，而且有些差別還比較大。為此可采用譜平衡技術來提高復賽譜域子波提取的精度，從而提高反褶積效果。

所謂譜平衡技術［2］，就是將信噪比高的近、中偏移距道作為標準道，以標準道的時頻特征為標準，將其他偏移距道的時頻特征［13］向其靠攏，使得遠近偏移距道的頻寬保持一致。具體步驟為：①將欲補償的道集進行時頻展開，并參考標準道的時頻譜，在時頻譜上逐點補償遠偏移距道所丟失的高頻信息，從而得到遠近偏移距道頻寬一致的道集；②進行濾波及道均衡，使各偏移距道的振幅能量保持一致。

1.2.2復賽譜統計“優化平均”

將復賽譜域子波提取技術用于疊前道集時，其重點在于如何求取空間一致的地震子波的振幅譜，使得其提取的子波是一種具有代表性的綜合子波（圖1）。對于一個三維資料，在特定的目的層段，選用一定的線距、道距和道集內不同的偏移距道作為分析樣板，即可求得一個子波集。由于該子波集中不同地震道的信噪比與頻譜特征均不相同，因此得出的子波空間一致性往往不太理想。為得到一個能代表整體數據特征的“典型”子波，筆者采用一種“剔除平均”的方法，根據井-震匹配和整體反褶積的效果，人為優選出一種子波，這樣提取的子波不僅具有全區的代表性，而且用其所做的反褶積抗噪性較強。

對比圖1（a）與圖1（b）可以發現，多道統計平均后得到的子波振幅譜主頻有所提高，頻帶得到了拓寬。從圖1（c）可以看出，多道“優選”提取的子波旁瓣得到了壓制，采用這種優選的子波進行反褶積處理能夠達到提高分辨率的目的。

圖1　多道子波振幅譜與單道子波振幅譜對比Fig.1　Comparison between multi-trace wavelet amplitude spectrum and single-trace wavelet amplitude spectrum

1.2.3人為定義希望輸出

常規確定性子波反褶積和時變譜白化技術都具有抗噪能力低的局限，但子波反褶積有一個明顯的優勢是能定義希望輸出子波。最常規的希望輸出子波為Ricker子波，其通過控制Ricker子波主頻來控制希望輸出頻譜的特征，但Ricker子波的頻帶較窄，不能適應靈活、復雜的解釋性處理要求。為此，筆者采用了一種人為定義希望輸出頻譜的方法，即手工定義希望輸出的子波振幅譜（圖2），使其主頻偏向低頻，頻帶得到拓寬，具有俞氏子波的特征，這樣就可以壓制高頻噪聲的抬升，將反褶積后道集的信噪比控制在合理的范圍內，從而提高低信噪比資料的分辨率和AVO屬性的穩定性。

對比圖2（a）與圖2（b）可看出，通過人為定義希望輸出子波振幅譜進行頻寬一致性子波反褶積處理后的振幅譜主頻有所提高，頻帶得到了拓寬，高頻噪聲未抬升，表明頻寬一致性子波反褶積技術具有一定的抗噪能力。

圖2　人為定義希望輸出子波振幅譜Fig.2　Artificial definition of the desired output wavelet amplitude spectrum

2　模型試算

為說明遠近偏移距道頻寬不一致對提高分辨率效果的影響，設計了如下實驗模型：地層厚度依次取800m，32.5m，30m，20m，20m，20m，20m，60m；對應的縱波速度依次為3100m/s，2500m/s，2800m/s，3000m/s，2700m/s，3300m/s，3000m/s，3400m/s；橫波速度依次為1280m/s，1265m/s，1255.023m/s，1224.745m/s，1653.406m/s，1697.836m/s，1224.745m/s，1749.286m/s；密度依次為2.3 g/cm3，2.12 g/cm3，2.01 g/cm3，2.12 g/cm3，2.07 g/cm3，2.21 g/cm3，2.10 g/cm3，2.20 g/cm3；道間距為40 m，雷克子波的主頻為50 Hz。

動校正后的道集如圖3（a）所示，可看出遠偏移距道波形“發胖”。圖3（b）是對圖3（a）模型進行復賽譜域子波反褶積處理后的效果，可看出遠偏移距道的分辨率未得到有效提高。圖3（c）是對圖3（a）模型進行頻寬一致性反褶積處理后的效果，可以看出遠、近偏移距道的分辨率均得到了有效提高。

圖3　模型數據試算Fig.3　Numerical experiments

3　實際資料應用

渤中X油田位于渤海灣盆地渤中凹陷與黃河口凹陷的分界處，該區沙三段目的層段為湖底扇沉積，發育的微相類型為辮狀水道，主要巖性為含礫砂巖、砂礫巖和泥巖互層。工區內有20多口井，為一開發井區，斷塊復雜，且橫向沉積相變快，儲層厚度為3～5 m，而地震主頻僅為29 Hz。因此，地震資料分辨率低是該區儲層預測的主要難題。

圖4為渤中X油田的實際道集資料。從疊前道集上看，近偏移距道有多次波，地震資料信噪比較低，主頻僅約20 Hz，原始地震資料的分辨率遠達不到薄層反演的要求。此外，由于動校正拉伸使得道集資料中遠近偏移距道頻寬不一致，遠偏移距道波形“發胖”，能量降低，給傳統反褶積方法的應用造成嚴重困難。

圖4　渤中X油田原始CRP道集及相應頻譜Fig.4　The seismic profile of CRP gather and its frequency spectrum of Bozhong X oilfield

對渤中X油田的實際道集資料進行疊前提高分辨率處理，首先采用均值濾波對原始道集（圖4）作去噪優化處理，去噪后的道集如圖5（a）所示；然后對去噪后的道集作復賽譜域子波反褶積處理［圖5（b）］和頻寬一致性子波反褶積處理［圖5（c）］。由于近道信噪比低，所以將譜平衡的模板道選擇為15～20道。對比圖5（b）和圖5（c）可以看出，復賽譜域子波反褶積處理后遠偏移距道的分辨率未得到顯著提高，而頻寬一致性子波反褶積處理后遠、近偏移距道的分辨率均得到了顯著提高。

圖5　CRP道集提高分辨率效果的對比Fig.5　Comparison of increased resolution of CRP gathers

與去噪處理后的疊加剖面［圖6（a）］相比，復賽譜域子波反褶積處理后的疊加剖面［圖6（b）］變化不大，而頻寬一致性子波反褶積處理后的疊加剖面［圖6（c）］目的層主頻由原來的20 Hz增大至40 Hz，頻帶拓寬為53 Hz，可達到的時間分辨率為12.2 ms，為后續將沙三段振幅屬性用于沉積相等研究提供了依據［14-16］。

圖6　疊加剖面提高分辨率效果對比Fig.6　Comparison of increased resolution of poststack profiles

圖7　振幅切片提高分辨率效果對比Fig.7　Comparison of increased resolution of poststack slices

將提取的沙三段振幅切片與沉積相（圖7）進行對比后可看出：原始振幅切片［圖7（a）］信噪比低，沉積模式不清晰；去噪后振幅切片［圖7（b）］信噪比有所提高，連續性得到增強，但沉積模式仍然模糊不清；復賽譜域子波反褶積振幅切片［圖7（c）］能夠顯現一定的沉積模式；頻寬一致性子波反褶積振幅切片［圖7（d）］與沉積相的模式［圖7（e）］對應最好，扇三角洲前緣泥巖發育區、扇根和扇中主體輪廓均較清晰，更好地反映了沉積相的平面展布。因此，頻寬一致性子波反褶積技術對提高低信噪比資料的分辨率具有更好的效果。

4　結論

（1）實驗模型證明，遠近偏移距道頻寬不一致，采用常規復賽譜域子波提取技術不能有效提高CRP道集的分辨率。為此，對復賽譜域子波提取技術進行了改進，并將它與道集譜平衡技術相結合，形成了頻寬一致性子波反褶積技術。

（2）實際資料應用證明，利用頻寬一致性子波反褶積技術提取的子波具有空間一致性，能夠有效壓制高頻噪聲的抬升，對提高CRP道集低信噪比資料的分辨率具有較好的效果。

（References）：

［1］鮑熙杰，趙海波.CRP道集優化處理對疊前反演的影響［J］.石油地質與工程，2013，27（1）：39-41. Bao Xijie，Zhao Haibo.CRP gathers optimization effect on prestack inversion［J］.Petroleum Geology and Engineering，2013，27（1）：39-41.

［2］劉力輝，楊曉，丁燕，等.基于巖性預測的CRP道集優化處理［J］.石油物探，2013，52（5）：34-39. Liu Lihui，Yang Xiao，Ding Yan，et al.Advanced process of CRP gather based on the prediction of lithology［J］.Geophysical Prospecting for Petroleum，2013，52（5）：34-39.

［3］趙偉，陳小宏，李景葉.薄互層調諧效應對AVO的影響［J］.石油物探，2006，45（6）：570-573. Zhao Wei，Chen Xiaohong，Li Jingye.The influence of thin interbed tuning effect of AVO［J］.Geophysical Prospecting for Petroleum，2006，45（6）：570-573.

［4］王萬里，李國發，桂金詠.混合相位子波有色反褶積［J］.巖性油氣藏，2013，25（3）：82-86. Wang Wanli，Li Guofa，Gui Jinyong.Colored deconvolution of mixed phase wavelet［J］.Lithologic Reservoirs，2013，25（3）：82-86.

［5］Barnes A E.Another look at NMO stretch［J］.Geophysics，1992，57（5）：749-751.

［6］李秉富，魏長江.地震資料動校正拉伸畸變分析［J］.青島海洋大學學報，1994（3）：39-44. Li Bingfu，Wei Changjiang.An analysis of NMO correction lengthening［J］.Journal of Ocean University of Qingdao，1994（3）：39-44.

［7］趙波，丈政軍.除動校正拉伸影響的方法［J］.石油地球物理勘探，1995，30（3）：417-421. Zhao Bo，Zhang Zhengjun.A method for removing NMO correction lengthening［J］.Oil Geophysical Prospecting，1995，30（3）：417-421.

［8］黃緒德.反褶積與地震道反演［M］.北京：石油工業出版社，1992：256-278. Huang Xude.De-convolution and seismic inversion［M］.Beijing：Petroleum Industry Press，1992：256-278.

［9］李國發，彭蘇萍，高日勝，等.復賽譜域提取混合相位子波的方法［J］.天然氣工業，2005，25（1）：85-87. Li Gufa，Peng Suping，Gao Risheng，et al.Method of mixed phase wavelet extraction in cepstrum domain［J］.Gas Industries，2005，25（1）：85-87.

［10］Ulrych T J.Application of homomorphic deconvolution to seismology［J］.Geophysics，1971，36（4）：650-660.

［11］Stoffa P L，Buhl P，Bryan G M.The application of homomorphic deconvolution to shallow-water marine seismology-PartⅠ：Models［J］.Geophysics，1974，39（4）：401-416.

［12］Xu Y，Chopra S.Improving AVO fidelity by NMO stretching and offset-dependent tuning corrections［J］.The Leading Edge，2007，26（12）：1548-1551.

［13］刁瑞，尚新民，芮擁軍，等.時頻域譜模擬反褶積方法研究［J］.巖性油氣藏，2013，25（1）：116-121. Diao Rui，Shang Xinmin，Rui Yongjun，et al.Study on spectral modeling deconvolution in time-frequency domain［J］.Lithologic Reservoirs，2013，25（1）：116-121.

［14］陳可洋，吳清嶺，李來林，等.松遼盆地三維地震資料連片處理關鍵技術及其應用效果分析［J］.巖性油氣藏，2012，24（2）：87-91. Chen Keyang，Wu Qingling，Li Lailin，et al.Key techniques and case study of Multi-survey joint processing of 3-D seismic data in Songliao Basin［J］.Lithologic Reservoirs，2012，24（2）：87-91.

［15］張猛剛，洪忠，竇玉壇，等.時頻分析在蘇里格地區含氣性檢測中的應用［J］.巖性油氣藏，2013，25（5）：76-85. Zhang Menggang，Hong Zhong，Dou Yutan，et al.Application of time-frequency analysis technology to the gas detection in Sulige area［J］.Lithologic Reservoirs，2013，25（5）：76-85.

［16］陶云光，袁剛，陳新安，等.二維地震連片相對保幅處理及儲層預測技術在塔西南地區的應用［J］.巖性油氣藏，2007，19（1）：96-100. Tao Yunguang，Yuan Gang，Chen Xin'an et al.Application of relative preserved amplitude of 2D Merging seismic data processing and reservoir prediction in Southwest Tarim［J］.Lithologic Reservoirs，2007，19（1）：96-100.

（本文編輯：王會玲）

Wavelet deconvolution with spatial consistency based on CRP gathers

Ming Jun1，Liu Lihui2，Ding Yan2，Li Cong2
（1.Tianjin Branch Company，CNOOC，Tianjin 300452，China；2.Beijing Rockstar Oil&Techology Company Ltd.，Beijing 100192，China）

When poststack wavelet extraction method of cepstrum domain is used for prestack wave extraction of CRP gathers and prestack raising resolution processing，the far and near bandwidth of the gather is not consistent，which influences the spatial consistency of wavelet and the effect of the resolution of prestack wavelet deconvolution.Therefore，this paper improved the wavelet extraction technology by using multi-cepstrum average，made the extraction wavelets have spatial consistency，and combined this improved wavelet extraction technology of cepstrum domain with spectral balance technology，to form a wavelet extraction technology which can keep the bandwidth consistent between the traces. This technology added artificial definition of the output of wavelet deconvolution，which makes it can adapt to prestack extraction frequency processing of the data with low signal to noise ratio.

wavelet extraction of cepstrum domain；spectral balance technology；wavelet deconvolution keeping bandwidth consistant；prestack raising resolution processing

P631.4

A

1673-8926（2015）06-0097-07

2015-07-20；

2015-09-17

中國海洋石油總公司“十二五”科技重大專項“渤海典型低孔低滲油藏勘探開發關鍵技術研究與實踐”（編號：CNOOC-KJ-125-ZDXM-07）資助

明君（1970-），男，博士，高級工程師，主要從事油田地球物理綜合解釋和方法研究。地址：（300452）天津市濱海新區塘沽閘北路1號609信箱渤海石油研究院。E-mail：mingjun@cnooc.com.cn。