多種頻率成分工業(yè)電干擾的壓制算法及其應用

陳可洋

（中國石油大慶油田有限責任公司勘探開發(fā)研究院，黑龍江大慶163712）

多種頻率成分工業(yè)電干擾的壓制算法及其應用

陳可洋

（中國石油大慶油田有限責任公司勘探開發(fā)研究院，黑龍江大慶163712）

為了解決現(xiàn)行工業(yè)化軟件中工業(yè)電干擾壓制處理的計算效率和處理精度較低的問題，提出并推導了一種壓制多種頻率成分工業(yè)電干擾的新算法。該算法采用余弦函數(shù)組合近似逼近工業(yè)電干擾，根據(jù)最小二乘法最優(yōu)逼近和歸一化互相關系數(shù)實現(xiàn)多種頻率成分工業(yè)電干擾的自動識別與自適應壓制。在理論合成數(shù)據(jù)和松遼盆地CD地區(qū)的地震資料工業(yè)電干擾壓制中進行應用，其結果表明，新算法的保幅性較好，信噪比較高，在勘探成熟區(qū)的高保真疊前地震資料預處理方面具有較好的推廣應用價值。

保幅預處理；工業(yè)電干擾；多種頻率成分；自動識別與自適應壓制；理論模型分析

0　引言

在松遼盆地勘探成熟區(qū)，地震資料品質(zhì)受到密集電網(wǎng)電磁場的影響，形成的工業(yè)電干擾將疊加在野外采集的地震資料中，這種工業(yè)電干擾的頻率約為50Hz，且與有效地震反射信號的優(yōu)勢頻帶重疊，進而降低了原始資料的信噪比。同時，工業(yè)電干擾主要表現(xiàn)為單頻特征，從淺層到深層，其相位、振幅和頻率均基本保持不變［1］。因此，如何高保真、高效地識別和壓制工業(yè)電干擾是當前地震資料保幅預處理中一個重要的技術問題［2］。早期采用的快速簡便的方法是帶限濾波法，其缺陷是丟失了工業(yè)電干擾頻率附近的有效地震信號，不利于保幅處理。為此，針對如何更好地壓制工業(yè)電干擾以滿足保幅預處理的要求，國內(nèi)外許多學者提出了多種新的處理思路：正余弦型函數(shù)逼近方法［3-4］、線性調(diào)頻譜消除方法［5］、基于互相關理論的單頻干擾壓制方法［6］、時間域單頻干擾波壓制方法［7］、Chirp-Z變換譜分析壓制方法［8］、獨立分量分析壓制方法［9］等。這些方法均對已知含有工業(yè)電干擾的地震道進行識別與壓制，且應用效果較明顯，但是對于輸入的任意地震道是否含有工業(yè)電干擾卻無法進行自動識別，也無法達到工業(yè)化推廣應用的目的。目前，常用的工業(yè)電干擾識別方法具有一定的局限性，且發(fā)展較為緩慢，主要思路是在共檢波點域計算均方根振幅能量基礎上［10-11］，再利用人工干預的方法設置門檻加以識別，但這種含工業(yè)電干擾地震道的識別與壓制的計算量較大，且計算效率較低。因此，現(xiàn)行工業(yè)化應用軟件中，仍沒有一種方法能快速簡便地對工業(yè)電干擾進行自動識別與壓制。

在前人研究的基礎上，筆者為了解決常規(guī)方法計算量大、計算精度較低的問題，探索形成一種基于歸一化互相關系數(shù)的含工業(yè)電干擾地震道的自動識別與自適應的工業(yè)電干擾壓制算法［12］，其處理效果優(yōu)于現(xiàn)行工業(yè)化應用軟件。在此基礎上，又將其進一步發(fā)展為適應多種頻率成分的工業(yè)電干擾壓制的新算法，在理論模型和地震資料處理中應用效果均較明顯。因此，新算法可為高精度、高保真疊前地震資料預處理提供重要的技術支撐。

1　方法原理

工業(yè)電干擾通常采用正弦函數(shù)或余弦函數(shù)表示，筆者以余弦函數(shù)為例［12］

式中：C，f和τ分別為工業(yè)電干擾的振幅（無量綱）、頻率（Hz）和時間延遲（無量綱）；Δt為時間采樣率，ms；i為時間采樣序號。

筆者將式（1）推廣到更一般的表達形式（由多種頻率成分工業(yè)電干擾組成），即

式中：k為工業(yè)電干擾成分的序號；K為工業(yè)電干擾的頻率成分總數(shù)量。

為了降低式（2）計算的復雜性，將該式表示為余弦函數(shù)和正弦函數(shù)的組合形式，即

式中：Ak和Bk分別為工業(yè)電干擾頻率fk的余弦函數(shù)和正弦函數(shù)的振幅，無量綱。

于是，建立最小二乘目標函數(shù)

式中：Si為原始地震數(shù)據(jù)；N為分析時窗內(nèi)地震數(shù)據(jù)采樣點個數(shù)。

當實際地震數(shù)據(jù)與預測的工業(yè)電干擾之間的誤差達到最小時，可以根據(jù)式（4）計算得到不同工業(yè)電干擾頻率成分經(jīng)最優(yōu)化逼近后的振幅Ak和Bk以及頻率fk，再根據(jù)正弦函數(shù)和余弦函數(shù)的換算關系，得到每種頻率成分工業(yè)電干擾余弦函數(shù)形式［式（2）］的振幅Ck和時間延遲τk以及相位。

筆者根據(jù)歸一化互相關系數(shù)的值域范圍［0，1］，設置工業(yè)電干擾的自動識別門檻值［12］（通常可取為0.2），從而自動識別地震資料中某地震數(shù)據(jù)道是否含有某種頻率成分的工業(yè)電干擾，這里的歸一化互相關系數(shù)的計算公式如下：

對含有某種頻率成分的工業(yè)電干擾記為yi，k′，將其從地震資料Si中減掉，得到多種頻率成分工業(yè)電干擾壓制后的有效地震信號，即

聯(lián)合式（1）～（6），最終計算得到高精度保幅處理的工業(yè)電干擾壓制結果。由于新算法是逐道實現(xiàn)工業(yè)電干擾的自動識別與自適應壓制處理的，因此，與疊前道集的域無關，同時不需要人工參與識別，輸入的數(shù)據(jù)要求為未做預處理的原始采集數(shù)據(jù)，這就大大簡化了常規(guī)工業(yè)化軟件中工業(yè)電干擾的壓制流程［13］。

2　理論模型測試

為了驗證筆者所提新方法的準確性和有效性，以含有不同振幅、相位和頻率的工業(yè)電干擾理論合成數(shù)據(jù)為例。模型共有336道記錄，每道記錄共有2種頻率成分的工業(yè)電干擾。圖1（a）為理論合成記錄中工業(yè)電干擾頻率在50 Hz附近的分布情況，從左到右其頻率從49.5 Hz增加到50.5 Hz，遞增間隔為1/6Hz；圖1（b）為理論合成記錄中工業(yè)電干擾頻率在100 Hz附近的分布情況，從左到右其頻率從99Hz增加到101Hz，遞增間隔為1/3Hz；圖1（c）和圖1（d）分別為理論合成記錄中2種頻率成分工業(yè)電干擾的振幅和相位分布情況。

圖1　理論合成數(shù)據(jù)工業(yè)電干擾壓制前后結果對比Fig.1 Com parison of theoreticalsynthetic data beforeand after industrial interference suppression

圖1（e）為2種頻率成分余弦函數(shù)組合得到的工業(yè)電干擾記錄，其工業(yè)電干擾的頻率、振幅及相位數(shù)據(jù)如圖1（a）～1（d）所示。圖1（f）和圖1（g）分別為新方法自動識別和自適應壓制掉的50 Hz和100 Hz附近處的工業(yè)電干擾記錄。圖1（h）為壓制工業(yè)電干擾后的有效信號部分，且不含有殘余工業(yè)電干擾。表1為新方法自動識別出部分工業(yè)電干擾振幅、頻率和相位與理論值的對比，誤差幾乎可以忽略。分析圖1和表1可知，筆者提出的新算法具有高精度識別和保幅壓制工業(yè)電干擾的優(yōu)點，驗證了該算法的準確性和有效性。

表1　工業(yè)電干擾頻率、振幅和相位的自動識別結果Table1 The automatic identification resultsof frequency，amplitude and phase of industrial interference

3　實際資料應用

以松遼盆地CD地區(qū)地震資料為例，工區(qū)內(nèi)高壓電線覆蓋范圍較廣，原始地震資料中疊加的工業(yè)電干擾較為嚴重［14-17］。該工區(qū)面積為656.26 km2，總炮數(shù)為29 375炮，地震數(shù)據(jù)總量為1.25 TB。經(jīng)統(tǒng)計，原始疊前炮集數(shù)據(jù)的總地震道數(shù)約為0.68億道，并存儲于一個SGY格式的文件中。由于該疊前地震數(shù)據(jù)量較大，常規(guī)的工業(yè)化軟件需要人工干預和多種處理模塊的串聯(lián)組合，所需磁盤空間較大，計算效率較低，因而難以實現(xiàn)高精度、高效率地開展對含工業(yè)電干擾地震道的識別和壓制處理。為此，筆者開發(fā)了一種新算法（并行計算程序）［13］，從而實現(xiàn)了對該工區(qū)含工業(yè)電干擾地震道的自適應識別與壓制處理。

圖2多種頻率成分工業(yè)電干擾壓制前和壓制后的疊前道集Fig.2 Pre-stack gather before and after multi-frequency com ponent industrial interference suppression

圖2為多種頻率成分工業(yè)電干擾壓制前、后的疊前道集及去掉的工業(yè)電干擾。圖2（a）為CD地區(qū)疊前地震道集中的一個原始采集共檢波點道集。分析可知，該疊前道集含有較強能量的工業(yè)電干擾，信噪比較低，掩蓋了有效的地震反射信號。采用筆者的新算法進行多種頻率成分工業(yè)電干擾的自動識別與壓制處理，在處理后的結果中［圖2（b）］，工業(yè)電干擾得到了有效壓制，有效反射信號更加清晰，信噪比明顯提高，并得到了50 Hz［圖2（c）］、150Hz［圖2（d）］、250Hz［圖2（e）］和350Hz［圖2（f）］左右的工業(yè)電干擾。

根據(jù)頻譜分析［圖3（a）］可知，圖2（a）不僅含有50Hz左右的工業(yè)電干擾頻率成分，還含有150Hz，250Hz和350Hz左右的工業(yè)電干擾頻率成分，常規(guī)工業(yè)電干擾處理方法無法高精度識別和壓制50Hz左右的工業(yè)電干擾頻率成分，對較高頻率成分的工業(yè)電干擾更難以實現(xiàn)自動識別與壓制，而采用新算法處理后的頻譜，在工業(yè)電干擾位置處過渡自然［圖3（b）］，無畸變；在壓制掉的工業(yè)電干擾頻譜中［圖3（c）］，清晰可見從50 Hz開始，以100 Hz為倍頻程的工業(yè)電干擾頻譜能量。

圖3　多種頻率成分工業(yè)電干擾壓制前（a）和壓制后（b）以及壓制掉的工業(yè)電干擾（c）的頻譜Fig.3 Spectrum before（a）and after（b）multi-frequency component industrial interference suppression and the suppressed industrial interference（c）

分析可知，工業(yè)電干擾的頻率并非以整數(shù)50Hz出現(xiàn)，而是在50 Hz左右的一定范圍內(nèi)變化。由此可見，按照常規(guī)工業(yè)電干擾整數(shù)頻率進行壓制的做法是錯誤的。

圖4為新算法在研究區(qū)50 Hz附近自動識別出的工業(yè)電干擾頻率分布結果，紅色代表工業(yè)電干擾頻率小于50 Hz，藍色代表工業(yè)電干擾頻率大于50Hz。圖4中展示的工業(yè)電干擾位置與該工區(qū)近地表的高壓電線分布較為吻合，分析可知，在較強的高壓電線網(wǎng)絡附近，地震資料受工業(yè)電干擾的影響程度較大。由此可見，新算法具有高精度保幅處理的工業(yè)電干擾壓制優(yōu)點。

圖4　50Hz附近自動識別出的工業(yè)電干擾頻率分布Fig.4 Frequency distribution of industrial interference identified automatically around 50Hz

4　結論

（1）筆者提出了一種能夠實現(xiàn)多種頻率成分工業(yè)電干擾壓制的新算法，推導了數(shù)值計算公式。與此同時，根據(jù)最小二乘最優(yōu)逼近和歸一化互相關系數(shù)，實現(xiàn)對含多種頻率成分工業(yè)電干擾地震道的自動識別與工業(yè)電干擾的自適應壓制。

（2）新算法是逐道實現(xiàn)對多種頻率成分工業(yè)電干擾的壓制，因此與疊前道集的域無關，且不需要人為干預，適合于海量疊前道集數(shù)據(jù)的工業(yè)電干擾處理。

（3）理論模型和地震資料應用均表明，新算法具有高精度保幅處理的工業(yè)電干擾壓制的優(yōu)點，因此，能夠滿足當前高精度、高保真疊前地震資料預處理的要求，具有較好的推廣應用價值。

（References）：

［1］凌云.地震數(shù)據(jù)采集·處理·解釋·一體化實踐與探索［M］.北京：石油工業(yè)出版社，2007. Ling Yun.Integrated practiceand exploration of seismic data acquisition，processingand interpretation［M］.Beijing：Petroleum Industry Press，2007.

［2］Linville A F，Meek R A.Canceling stationary sinusoidal noise［J］. Geophysics，1992，57（11）：1493-1501.

［3］胡偉，呂小偉.壓制地震資料中工業(yè)電干擾的余弦逼近法的改進及應用［J］.物探與化探，2005，29（6）：537-540. HuWei，LüXiaowei.The improvement of cosine approachmethod and its application to remove industrial noise in seismic data［J］. Geophysicaland GeochemicalExploration，2005，29（6）：537-540.

［4］高少武，趙波，祝樹云，等.余弦函數(shù)自適應法識別與消除單頻干擾［J］.石油地球物理勘探，2011，46（1）：64-69. Gao Shaowu，Zhao Bo，Zhu Shuyun，etal.Cosine function adaptive method to identify and eliminate single frequency interference［J］. OilGeophysicalProspecting，2011，46（1）：64-69.

［5］高少武，羅國安，趙波，等.利用線性調(diào)頻譜法識別與消除單頻干擾［J］.石油地球物理勘探，2010，45（6）：861-867. Gao Shaowu，Luo Guoan，Zhao Bo，etal.Single frequency interference identification and elimination with linear adjustable spectrum method［J］.OilGeophysicalProspecting，2010，45（6）：861-867.

［6］高少武，趙波，祝樹云，等.自相關法單頻干擾識別與消除方法［J］.地球物理學報，2011，54（3）：854-861. GaoShaowu，ZhaoBo，Zhu Shuyun，etal.Single frequency interference identification and eliminationmethod with autocorrelation［J］.Chinese journalofgeophysics，2011，54（3）：854-861.

［7］高少武，賀振華，趙波，等.時間域單頻干擾波消除方法的改進［J］.石油地球物理勘探，2008，43（3）：270-274. Gao Shaowu，He Zhenhua，Zhao Bo，etal.Improvement single frequency interferencewave elimination method in time domain［J］. OilGeophysicalProspecting，2008，43（3）：270-274.

［8］羅國安，高少武，魏庚雨，等.Chirp-Z變換譜分析壓制地震記錄單頻干擾［J］.石油地球物理勘探，2009，44（2）：166-172. Luo Guoan，Gao Shaowu，WeiGengyu，etal.Suppression of seismic record single frequency interferencewith Chirp Z transform spectrum analysismethod［J］.Oil Geophysical Prospecting，2009，44（2）：166-172.

［9］吳小培，詹長安，周荷琴，等.采用獨立分量分析方法消除信號中的工頻干擾［J］.中國科學技術大學學報，2000，30（6）：671-676. Wu Xiaopei，Yan Changan，Zhou Heqing，etal.Using independent componentanalysismethod to eliminate the power frequency interference signal［J］.Journalof University of Science and Technology ofChina，2000，30（6）：671-676.

［10］步長城，趙志萍，王廣勇，等.在檢波點域分離50Hz工業(yè)干擾［J］.油氣地球物理，2005，3（3）：20-21. Bu Changcheng，Zhao Zhiping，WangGuangyong，etal.50 Hz industrial interference separation in common receiver point domain［J］.Petroleum Geophysics，2005，3（3）：20-21.

［11］李慧，王依川，張春堂，等.地震記錄中工業(yè)電干擾噪聲識別與壓制［J］.大慶石油地質(zhì)與開發(fā)，2012，31（5）：170-174. Li Hui，Wang Yichuang，Zhang Chuntang，et al.Recognition and depression of the interfering noise of the industrial electricity in seismic record［J］.Petroleum Geology&Oilfield Development in Daqing，2012，31（5）：170-174.

［12］陳可洋，陳樹民，李來林，等.單頻干擾的高精度自動識別和自適應壓制方法［J］.巖性油氣藏，2014，26（3）：109-113. Chen Keyang，Chen Shumin，Li Lailin，etal.High-accuracy automatic identification and adaptive suppressionmethod for single frequency interference［J］.Lithologic Reservoirs，2014，26（3）：109-113.

［13］陳可洋.壓制工業(yè)電干擾并行處理模塊的開發(fā)與應用［J］.油氣地球物理，2015，13（2）：1-6. Chen Keyang.Developmentand application ofparallelprocessing moduleon suppressing industrial interference［J］.Petroleum Geophysics.2015，13（2）：1-6.

［14］陳可洋，吳清嶺，李來林，等.松遼盆地三維地震資料連片處理關鍵技術及其應用效果分析［J］.巖性油氣藏，2012，24（2）：87-91. Chen Keyang，Wu Qingling，Li Lailin，etal.3D seismic datamultisurvey joint processing key technology and its application effect analysis in SongLiao Basin［J］.Lithologic Reservoirs，2012，24（2）：87-91.

［15］陳可洋，吳沛熹，楊微.擴散濾波方法在地震資料處理中的應用研究［J］.巖性油氣藏，2014，26（1）：117-122. Chen Keyang，Wu Peixi，YangWei.Application study on the diffusion filteringmethod in the seismic data processing［J］.Lithologic Reservoirs，2014，26（1）：117-122.

［16］陳可洋，楊微，吳清嶺，等.地震反射波與散射波波場分離方法初探［J］.巖性油氣藏，2013，25（2）：76-81. Chen Keyang，YangWei，Wu Qingling，etal.Preliminary study on seismic reflection wave and scattered wave wavefield separating scheme［J］.Lithologic Reservoirs，2013，25（2）：76-81.

［17］陳可洋，吳清嶺，林春華，等.含壓實效應的沙丘曲線靜校正方法及其應用［J］.礦業(yè)工程研究，2011，26（4）：64-68. Chen Keyang，Wu Qingling，Lin Chunhua，et al.Sand dune curve staticmethod considering compaction effectand itsapplication［J］. Mineral AngineeringResearch，2011，26（4）：64-68.

（本文編輯：楊琦）

Suppressing schemeofmulti-frequency com ponent industrial interferenceand itsapplication

CHEN Keyang
（Research Institute of Exploration and Development，PetroChina Daqing Oilfield Company Ltd.，Daqing 163712，Heilongjiang，China）

In order to solve the problemsof the low calculating efficiency and low processing precision of the current industrial interference suppression process in industrial software，this paper presented and derived a new scheme to suppressmulti-frequency component industrial interferences simultaneously.The algorithm uses a combination of cosine functions toapproximate the industrialinterference，and achievesautomatic identification and adaptive suppression ofmulti-frequency component industrial interference according to the optimal approximation of least square method and the normalized cross-correlation coefficient.Itwas applied to theoretical syntheticmodel and practical seismic data in CD area of Songliao Basin.the result shows that the proposed scheme is with perfect amplitudepreserved and high signal to noise ratio，therefore，the new algorithm can be an important technicalguidance for high precision amplitudepreserved seismic datapreprocessing inmatureexploration area.

seismic amplitude-preserved preprocessing；industrial interference；multi-frequency component；automatic identification and adaptive suppression；theoreticalmodelanalysis

P631

A

1673-8926（2015）02-0092-06

2014-10-06；

2014-12-25

國家重點基礎研究發(fā)展計劃（973）項目“火山巖油氣藏的形成機制與分布規(guī)律”（編號：2009CB219307）資助

陳可洋（1983-），男，碩士，工程師，主要從事高精度地震波傳播模擬與逆時成像、高性能集群并行程序開發(fā)、實際地震資料數(shù)字處理方法研究與應用等。地址：（163712）黑龍江省大慶市讓胡路區(qū)大慶石油勘探開發(fā)研究院地震處理二室。電話：（0459）5508524。E-m ail：keyangchen@163.com。