劉玉海，潘樹林，胡永泉，高 磊，沈 姣

（1.西南石油大學(xué) 資源與環(huán)境學(xué)院，四川 成都 610500； 2.西南石油大學(xué) 計算機科學(xué)學(xué)院，四川 成都 610500）

0 引言

油氣井射孔過程破壞地下圍巖應(yīng)力平衡狀態(tài)，引起巖石破裂并震動而形成彈性聲波，被埋置于地面的檢波器接收記錄為地震信號，其能量很弱，震級為-1級左右，稱為微地震信號［1-3］.微地震信號大多掩蓋在環(huán)境噪聲下，該環(huán)境噪聲中不僅含有隨機干擾噪聲，還包含多種規(guī)則干擾噪聲，糾合在一起，改變微地震信號形態(tài)，難以準確地辨別微地震剖面上有效信號的同相軸［4-6］.維納方程能夠從隨機干擾噪聲中提取有用的微地震信號而不損失有效信號的頻率成分.針對微地震相波場的復(fù)雜性，文中采用基于維納方程的相關(guān)函數(shù)迭加法檢測微地震信號，該方法能夠有效地衰減微地震信號中的隨機干擾噪聲，突出微地震信號的同相軸，對地面采集的微地震信號有較好的去噪效果，處理結(jié)果達到較高信息保真度.

1 方法原理

1.1 微地震道權(quán)系數(shù)計算

常規(guī)地震信號鄰域加權(quán)檢測法將所有鄰域采樣點權(quán)系數(shù)作相同的假設(shè)不完全合理.一定時窗內(nèi)相鄰道微地震信號相關(guān)性最好［7］，但直接使用相關(guān)函數(shù)法［8］計算權(quán)系數(shù)時，常受噪聲相關(guān)函數(shù)的影響.微地震信號干擾噪聲壓制方法的差異主要表現(xiàn)為權(quán)系數(shù)的計算.

為避免噪聲相關(guān)性影響，可采用遠離零延遲的相關(guān)函數(shù)計算微地震信號檢測的權(quán)系數(shù)［9］，但結(jié)果易造成微地震信號部分能量損失，難以達到提高微地震信號分辨率和突出微地震信號同相軸的目的［10-11］.維納方程法能夠把功率譜全部或部分重疊的微地震信號從隨機干擾噪聲中檢測出來，保留地震信號有效成分的同時使噪聲得到衰減［12-13］.文中利用維納方程的這一特性計算微地震信號的權(quán)系數(shù).

維納方程法應(yīng)用最小平方原理求權(quán)系數(shù)w（m）.為使真實值與估計值的平方誤差最小，對式

求w（m）的偏導(dǎo)數(shù)，并令其等于0，獲得權(quán)系數(shù)的線性方程組：

式中：Rxx為x（n）的自相關(guān)矩陣，Rsx為x（n）和s（n）的互相關(guān)矩陣.即

已知Rxx和Rsx時，由式（2）和式（1）求出微地震道的權(quán)系數(shù)W.

1.2 多道微地震信號分段小時窗處理

式（2）中權(quán)系數(shù)的計算基于統(tǒng)計規(guī)律，只有微地震信號采樣點足夠多時，才能充分衰減微地震記錄中隨機干擾噪聲，清晰突出微地震信號的同相軸.為增加微地震信號的計算采樣點數(shù)，可利用維納方程法計算相鄰的多道微地震信號權(quán)系數(shù)，利用權(quán)系數(shù)對微地震信號進行褶積運算和迭加濾波處理，表示為

式中：k為迭加濾波的時移值；x（n-k）為微地震信號的采集值.

為使得到的權(quán)系數(shù)具有統(tǒng)計特性，采樣點數(shù)N0選取應(yīng)盡可能大，可以利用整道微地震信號進行維納方程權(quán)系數(shù)計算；但該方法易導(dǎo)致弊端：（1）隨時間變化，整個微地震剖面上可能存在2個或2個以上的同相軸，使最終檢測結(jié)果為多個同相軸的共同反映，隨權(quán)系數(shù)變化可能具有平緩或多個極值，難以根據(jù)權(quán)系數(shù)運算結(jié)果判斷同相軸走向，降低微地震信號同相軸的分辨率；（2）最終檢測結(jié)果不能指示同相軸的位置和相對強度，不能顯示不同同相軸的形態(tài).

為避免弊端，采用分段小時窗法處理整個測區(qū)的微地震信號，表示為

式中：wm（k）為第m個時窗內(nèi)微地震信號檢測的權(quán)系數(shù)；xm（n-k）為第m個時窗內(nèi)微地震信號的采集值；（n）為第m個時窗內(nèi)微地震信號的檢測值；M為窗口寬度；N為窗口長度.

檢測微地震信號時，每次將M×N個測點所組成的測線段作為一個時窗，利用維納方程法計算權(quán)系數(shù)，對微地震信號進行檢測處理；然后逐點推移整個時窗的測線段，直至遍及整個測區(qū)［17-18］.

針對不同的窗口參數(shù)M和N，隨著隨機干擾噪聲方差的增大，對應(yīng)的窗口最佳參數(shù)隨之增大.對微地震信號進行小時窗相關(guān)迭加處理時，除考慮窗口長度N隨隨機干擾噪聲方差變化外，還需考慮窗口寬度M的選擇與微地震信號的穩(wěn)定性相關(guān).對于走向長度較大且比較穩(wěn)定的微地震同相軸，較大窗口寬度對隨機干擾噪聲的衰減效果較好.如果同相軸的走向長度較小或走向變化較大，窗口寬度M選太大將影響同相軸清晰度、降低對同相軸走向變化的分辨能力［19］.

多道微地震信號分段小時窗處理法還應(yīng)考慮計算量的影響.時窗選取的越小計算量越大，在不影響微地震信號處理效果前提下，應(yīng)盡量使用較大的時窗［20］.使用該方法對微地震信號進行迭加運算時，有時難以判斷同相軸形態(tài)，可以進行不同時窗的權(quán)系數(shù)運算，選取濾波效果最佳的時窗參數(shù).

2 仿真分析

為驗證維納方程相關(guān)函數(shù)迭加法檢測地面采集的微地震信號效果的有效性，建立模擬微地震信號（見圖1）和加入45%隨機干擾噪聲的微地震信號（見圖2）.

由圖2可見：加入與模擬微地震信號同等強度的隨機干擾噪聲后，微地震信號表現(xiàn)為噪聲能量強，信號信噪比降低，信號清晰度差，信號的同相軸模糊不清.

利用維納方程相關(guān)函數(shù)迭加法對隨機干擾噪聲的去噪效果見圖3.由圖3可見：對加入隨機干擾噪聲的地震信號進行濾波處理后，微地震信號的同相軸被清晰地顯示出來.表明使用維納濾波相關(guān)函數(shù)迭加法能夠衰減微地震信號中的隨機干擾噪聲，提高信號的信噪比，并且保持信號特征不變，證明該方法檢測微地震信號可行.

3 實例應(yīng)用

為驗證維納方程相關(guān)函數(shù)迭加法實際應(yīng)用效果，使用該方法處理地面采集的微地震信號（見圖4）.微地震信號包含200條地震道、1 000個采樣點，采樣間隔為1ms.

由圖4可見：原始微地震信號中隨機干擾噪聲很強，同微地震信號迭加在一起，降低微地震信號的分辨率，難以準確判斷微地震信號同相軸的走向和形態(tài).對原始微地震信號進行分段小時窗剖分后，在每個小時窗內(nèi)利用維納方程相關(guān)函數(shù)迭加法對其進行檢測，檢測結(jié)果見圖5.

由圖5可見：經(jīng)過維納方程相關(guān)函數(shù)迭加法濾波處理后的微地震信號，信噪比有所提高，微地震信號同相軸的走向、形態(tài)和分辨率比圖4的突出和清晰.說明對地面采集的微地震信號采用維納方程相關(guān)函數(shù)迭加法檢測處理，其方法有效.

圖4和圖5微地震信號采用維納方程相關(guān)函數(shù)迭加法檢測前后的平均振幅譜見圖6.

由圖6可見：地面采集的微地震信號經(jīng)過維納方程相關(guān)函數(shù)迭加法濾波后，高能量的隨機干擾噪聲被衰減，保留部分有效的微地震信號.對微地震信號濾波前后的平均振幅譜對比表明，利用維納方程相關(guān)函數(shù)迭加法能夠衰減混合在微地震信號中的隨機干擾噪聲，對微地震信號的信噪比起到很好的改善作用.

4 結(jié)束語

維納方程相關(guān)函數(shù)迭加法能夠衰減地面檢波器接收的微地震信號中的隨機干擾噪聲，提高微地震信號的信噪比，是在各種噪聲干擾背景上突出微地震信號同相軸的有效方法.

根據(jù)地面接收的微地震信號特點，選擇最佳的分段小時窗相關(guān)函數(shù)法迭加參數(shù)，能夠突出微地震信號的同相軸和提高微地震信號的分辨率，使處理后信號保真度不變.結(jié)果表明，維納方程相關(guān)函數(shù)迭加法可衰減噪聲中的隨機干擾噪聲，檢測出微地震信號中的規(guī)則波形信號，是微地震信號處理中有效的方法.

［1］李智敏，茍先太，金煒東，等.微地震信號的頻率特征［J］.地質(zhì)工程學(xué)報，2010，30（6）：106-112.

［2］王利亞，印春生.強噪聲背景中微弱信號檢測的初步研究［J］.分析化學(xué)，1999，12（27）：128-133.

［3］Leo E，Michael T.Challenges for microseismic moniting［J］.Geophysics，2010，75（5）：146-154.

［4］David W E.Q determination，corner frequency and spectral characteristics of microseismicity induced by hydraulic fracturing［C］.San Antonio：SEG Annual Meeting，2011：1666-1673.

［5］王有新.應(yīng)用地震數(shù)據(jù)處理方法［M］.北京：石油工業(yè)出版社，2009：63-70.

［6］劉俊，張斌珍.微弱信號檢測技術(shù)［M］.北京：電子工業(yè)出版社，2005：86-93.

［7］程佩青.數(shù)字信號處理［M］.2版.北京：清華大學(xué)出版社，2003：546-587.

［8］石穎，劉洪.地震信號的復(fù)地震道分析及應(yīng)用［J］.地球物理學(xué)進展，2009，23（5）：134-141.

［9］侯榆青，楊旭朗，韓敏.基于偏微分方程與維納濾波的混合去噪方法［J］.計算機工程，2010，36（10）：193-198.

［10］張?zhí)扈?基于維納濾波的小波圖像去噪算法研究［J］.鄭州大學(xué)學(xué)報：自然科學(xué)版，2011，36（1）：92-96.

［11］崔曉杰.維納濾波的應(yīng)用研究［J］.勘探地球物理進展，2009，32（3）：182-187.

［12］顏君，鄭寶玉.改進的維納濾波在地震資料處理中的應(yīng)用［J］.吉林大學(xué)學(xué)報：自然科學(xué)版，2009，23（4）：108-117.

［13］段少華，顧學(xué)邁，張乃通.通信系統(tǒng)仿真中的信噪比估計［J］.系統(tǒng)工程與電子技術(shù)，2004，21（4）：122-128.

［14］李鏘，滕建輔.小波域 wiener濾波器信號的去噪方法［J］.天津大學(xué)學(xué)報，2003，36（3）：334-337.

［15］牟永光.地震資料數(shù)學(xué)處理方法［M］.2版.北京：石油工業(yè)出版社，2006.122-167.

［16］邱天爽，劉文紅，郭瑩，等.現(xiàn)代數(shù)字信號處理與噪聲降低［M］.3版.北京：電子工業(yè)出版社，2009：78-95.

［17］高晉占.微弱信號檢測［M］.北京：清華大學(xué)出版社，2004：69-76.

［18］蔡希玲.聲波和強能量干擾的分頻自適應(yīng)檢測與壓制方法［J］.石油地球物理勘探，2008：34（8）：90-95.

［19］Farnoush，F(xiàn)orghani-Arani，Mark W.Analysis of passive surface-wave noise in surface microseismic data and its implications［C］.San Antonio：SEG 2011Annual Meeting，2011：1493-1499.

［20］David E D，Nanoseis，Stephen P.Gardner.Comparison of simultaneous downhole and surface microseismic monitoring in the Williston Basin［C］.San Antonio：SEG 2011Annual Meeting，2011：1687-11674.