基于三阶累积量一维切片的微地震信号时延估计
摘 要：由于微地震有效信号易受环境噪声的干扰，严重影响时延估计值获取的准确性.基于此，提出了基于三阶累积量一维切片的方法进行微地震信号时延估计.首先计算信号的三阶自累积量和互累积量，为减少运算量满足野外实时要求取其一维切片，最终通过准则函数求取时延值.相对于传统的互相关方法，该方法对相关性未知的高斯噪声具有不敏感的特点，对相关高斯噪声有更强的抑制能力.仿真结果表明，在添加信噪比为-5～15 dB的随机高斯噪声和相关高斯噪声时，其时延检测概率和均方根误差优于互相关方法.在实际资料应用中，对实际微地震噪声进行了概率分布统计分析，并且基于井地联合监测方式，选取井中信号作为参考道对地面资料进行时延值的求取.仿真和实际资料处理验证了该方法的有效性和实用性.
关键词：参数估计；噪声抑制；三阶累积量；微地震监测；一维切片
中图分类号：TE319 文献标识码：A
Abstract：The effective microseismic signals are disturbed by noises， so the precision of time delay estimation is affected. Because time delay estimation is the key to determine location accuracy， one dimensional slice method based on third-order cumulants was proposed for microseismic signal time delay estimation. Firstly， the three-order self accumulation and mutual accumulation of the signals were calculated， then， one dimensional slice was obtained in order to reduce the computation and to satisfy the field real-time requirements. Finally， the time delay value was obtained by the criterion function. Compared with the traditional cross correlation method， this method has less sensitivity to the Gauss noise with unknown correlation and has stronger inhibition ability to the correlated Gauss noise. The simulation results show that the detection probability and root mean square error are better than the traditional method under the influence of Gauss random noise and correlated noise with SNR=-5 dB to SNR=15 dB. In the implementation of real data， the real microseismic noise probability statistic characteristics were analyzed and the data in the well were selected as references， and then the time delay value to the ground data was calculated on the basis of the surface and borehole observation. The validity and the practicability of this method were verified.
Key words：parameter estimation； noise suppression； three-order cumulants； microseismic monitoring； one dimensional slice
微地震监测技术是一门交叉学科的新技术，它是通过观测、分析生产活动中所产生的微小地震事件来监测生产活动中产生的影响、效果及地下状态的地球物理技术[1].微地震监测技术在油田水力压裂监测、煤层气压裂监测、矿山与地质灾害监测等领域发挥着越来越重要的作用[2-4].微地震监测分为井中监测和地面监测，地面监测受到的环境干扰很大，相比于地面监测方式井中监测获得的微地震信号信噪比要高，基于井中地面联合监测的方式近些年来得到广泛应用[5-6].
与常规地震勘探信号相比，微地震信号十分微弱，震级大约在-3级到+1级，有效地压制噪声获得时延估计值对后续的微地震信号识别和震源定位等工作至关重要.Eisner等[7-8]提出了一种基于“主事件”（信号明显）与“从事件”（信号不明显）互相关运算的方法来压制噪声并进行震源定位.但是，相关运算的前提是噪声与噪声之间是不相关的，而实际的微地震资料中噪声的相关性未知，且可能具有很强的相关性，这样，相关方法获得的时延估计值是不准确的.
随着高阶统计量理论的逐步完善，由于高阶累积量对高斯噪声不敏感等特性，在地球物理领域尤其在地震信号方面获得了成功的应用[9-10].冯智慧等[11]将四阶累积量应用到地震信号初至拾取中，拾取精度不亚于人工方法，并且效率大大提高；熊晓军等[12]将高阶统计量应用在油气检测中；唐斌等[13]在地震子波恢复也获得了成功应用；宋维琪等[14]研究了基于贝叶斯框架的四阶累积量的自适应算法能够有效地提取出微地震的有效信号. 本文基于微地震井中―地面联合监测方式，以井中信号为参考道，可以降低时延估计过程中的累积误差；运用三阶累积量一维切片对微地震信号进行时延估计，能够降低其估计方差，减少运算量，并且能够有效压制微地震高斯噪声，且获得准确的时延值.理论模型证明了该方法的可行性并且优于互相关方法，实际资料也获得了较好的应用.
2.3 参考道的选取
参考道的选取一般有3种途径：1）选取固定的一道作为参考道，然后进行高阶累积量的时延估计运算；2）每次都选取前一道作为参考道，然后运用高阶累积量对相邻道进行时延估计运算；3）把相似道先进行分组，再进行高阶累积量时延估计运算[11].
实际的微地震信号很微弱，信噪比低，尤其是地面微地震资料.相对于地面微地震资料，井中资料的信噪比较高，容易识别出有效信号.基于此，采用了基于井中地面联合微地震监测方式，以井中资料作为参考道，对地面微地震资料进行时延估计，压制高斯噪声，获得了更准确的时延值.
2.4 基于三阶累积量一维切片方法的微地震信号
时延估计步骤
1）首先截取井中与地面相同时间段和相同时窗长度的数据，其中包含微地震有效信号.这样既可以减少计算量，也避免了无效微震信号噪声的干扰.
2）对截取的数据段进行归一化、去均值、带通滤波等预处理，使其满足高阶累积量运算假设.
3）以井中数据作为参考道，运用上述方法与地面各道信号进行时延估计运算，将得到的值保存.
4）得到的时差值可用于接下来的初至拾取和震源定位处理中.
3 微地震资料噪声分析
微地震监测环境噪声复杂，产生的干扰波可能与微地震传播机理和环境等因素有关.总体可以分为相干噪声和随机噪声.在理想情况下，无限范围的地震噪声根据中心极限定理总体呈现高斯分布[20].为了适用于本文所提方法，首先分析了实际噪声的概率分布情况.图1为2014年在山西进行煤层气压裂监测压裂结束后采集的噪声干扰，截取其一段数据中6道噪声数据进行概率分布分析得到图2（柱状图为实际噪声的概率分布）.
由图2可以看出，有限数据的噪声概率密度分布大部分图2（a）-图2（e）呈“钟形”，基本满足高斯分布特征.图2（e）由于大振幅短时瞬态信号的存在，出现了一个峰度值.图2（f）出现了双峰度，这是由于现场存在大振幅的周期性噪声信号.从对噪声的概率统计特性分析可以得出：当噪声中存在一些窄带强噪声或大振幅的周期性噪声时呈现非高斯分布，与猝发性环境有关，主要是野外施工中机械振动和人为振动产生的干扰，它是少样本随机分布的；而高斯噪声是在全范围内存在的，具有持久保持性.因此，实际监测现场环境噪声信号的分析结果表明该应用场合满足本文所述方法所需的理论基础.
4 仿真实验
仿真中采用地震勘探领域常用的Ricker子波模拟有效信号，将其合成微地震记录来验证提出算法的可靠性.Ricker子波的主频设为50 Hz，合成记
录道集为30道，每道间的时间延迟为10 ms，如图3（a）所示.为了验证基于三阶累积量一维切片的方法在去除高斯噪声方面优于相关方法，分别在原始记录上添加随机高斯噪声和相关高斯噪声，分别如图3（b）和图3（c）所示.可以看出，相关噪声产生了一些假同相轴.
以第1道为参考道，用后续的道与参考道进行本文方法的时延估计运算.选取第2道为例与参考道进行运算，其理论时延值为10 ms.图4为分别添加随机噪声和相关噪声的情况下，本文方法和传统互相关方法通过1 000次Monte-Carlo实验获得的不同信噪比下的时延检测概率.由图4可知，无论是在随机噪声还是相关噪声的影响下，本文方法都优于互相关方法，尤其在相关噪声干扰下，互相关方法在低于0 dB是基本失效.图5为分别添加随机噪声和相关噪声情况下，本文方法和传统互相关方法通过1 000次Monte-Carlo实验获得的不同信噪比下的均方根误差.由图5可知，本文方法均方根误差值小于互相关方法结果，同样通过定量分析证明了本文方法的优势.图6为含有相关噪声的数据在信噪比为-3 dB时运用互相关方法和本文方法进行的一次时延估计结果.由图6中圆圈处可知，基于互相关方法的最大值的时延值在零点左侧，这是相关噪声的影响，干扰了实际时延值的拾取，而基于三阶累积量一维切片方法可以正确计算出时延值，受相关高斯噪声影响很小，并且三阶累积量一维切片方法获得的结果的旁瓣要比互相关方法小，提高了检测分辨率.在实际资料处理中，信噪比极低的情况下也要对资料进行滤波等预处理，尽量提高信噪比.
5 实际资料处理
资料处理采用2014年在山西某地进行煤层气压裂监测的实际数据，仪器系统均为本单位自主研发，其中，井中监测采用五支高精度三分量采集短节，地面采集利用无缆遥测地震仪，外接与井中匹配的高精度三分量检波器.截取了相同时间段的井中和地面数据，图7（a）为井中采集的五级三分量原始数据，信噪比较高；图7（b）为16个地面采集站的原始数据，明显看出其受到了严重的噪声干扰.应用上述方法，选取井中第2道Y分量的数据和地面的第2道Y分量数据，分别对其进行去均值、归一化、带通滤波等预处理，结果如图8（a）和8（b）所示.运用本文提出的三阶累积量一维切片方法对井中和地面数据进行时延估计，所得结果如图8（c）所示.从图中可以看出三阶累积量一维切片方法能够相对准确地求取出时延估计值.图8（d）为运用互相关方法获得的时延值，通过对比可以看出本文方法获得时延值的旁瓣更小，分辨率更高，时延值获取更准确.本文方法得到的每个时延估计值都是地面资料相对于井中资料的各个分量的时间差值，可以进一步求出地面资料各个分量的时延值，根据同行波时差定位法就可以反演出震源位置.
6 结 论
微地震监测资料受噪声干扰严重，本文提出的基于三阶累积量一维切片的方法能够有效地去除高斯噪声，获得时延估计值.通过仿真分析可以知道在添加随机高斯噪声和相关高斯噪声情况下，本文方法在时延拾取概率和均方根误差方面优于常规的相关方法，尤其在相关高斯噪声存在的情况下更能体现其优势.通过对实际微地震噪声的统计特性分析可知，实际噪声具有很强的高斯性.在实际的微地震资料处理中，选取了井中数据作为参考道，避免了选取地面数据作为参考道而受到噪声干扰大而获得的时延值不准确的问题.在信噪比很低的情况下，通过滤波等预处理可以更好地获得时延估计值.获得的时延估计值可以进一步得到地面不同道之间有效信号的时间差值，是后续微地震震源定位运算中的重要参数. 参考文献
[1] 宋维琪， 陈泽东， 毛中华. 水力压裂裂缝微地震监测 技术[M]. 东营： 中国石油大学出版社， 2008： 1-3.
SONG Wei-qi， CHEN Ze-dong，MAO Zhong-hua. Hydraulic fracturing cracks microseismic monitoring technology[M]. Dongying： China University of Petroleum Press， 2008： 1-3.（In Chinese）
[2] 梁兵， 朱广生. 油气田勘探开发中的微地震监测方法[M]. 北京： 石油工业出版社， 2004： 81-100.
LIANG Bing， ZHU Guang-sheng. Microseismic monitoring methods in oil and gas field exploration and development[M]. Beijing： Petroleum Industry Press， 2004： 81-100.（In Chinese）
[3] 姜福兴， 杨淑华， 成云海， 等. 煤矿冲击地压的微地震监测研究[J]. 地球物理学报， 2006， 49（5）： 1511-1516.
JIANG Fu-xing， YANG Shu-hua， CHENG Yun-hai，et al. A study on microseismic monitoring of rock burst in coal mine[J]. Chinese Journal of Geophyics， 2006， 49（5）：1511-1516.（In Chinese）
[4] 柳云龙， 田有， 冯噎， 等. 微震技术与应用研究综述[J].地球物理学进展， 2013， 28（4）： 1801-1808.
LIU Yun-long， TIAN You， FENG Xuan， et al. Review of micmseism technology and its application[J]. Progress in Geophysics， 2013， 28（4）： 1801-1808.（In Chinese）
[5] 李振春， 盛冠群， 王维波， 等. 井地联合观测多分量微地震逆时干涉定位[J]. 石油地球物理勘探， 2014， 49（4） ： 661-666.
LI Zhen-chun， SHENG Guan-qun， WANG Wei-bo， et al.Time-reverse microseismic hypocenter location with interferometric imaging condition based on surface and downhole multi-components[J]. Oil Geophysical Prospecting，2014， 49（4）：661-666.（In Chinese）
[6] 王晨龙， 程玖兵， 尹陈， 等.地面与井中观测条件下的微地震干涉逆时定位算法[J]. 地球物理学报， 2013， 56（9）：3184-3196.
WANG Chen-long， CHENG Jiu-bing， YIN Chen， et al. Microseimic events location of surface and borehole observation with reverse-tine focusing using interferometry technique[J]. Chinese Journal of Geophyics， 2013， 56（9）： 3184-3196.（In Chinese）
[7] EISNER L，ABBOTT D，BARKER W B，et al.Noise suppression for detection and location of microseismic events using a matched filter[C]//Society of Exploration Geophysicists. Las Vegas， Nevada：SEG Annual Meeting， 2008：1431-1435.
[8] GIBBONS S J， RINGDAL F. The detection of low magnitude seismic events using array-based waveform correlation[J]. Geophysical Journal International， 2006，165：149-166.
[9] 李希亮， 刘希强， 董晓娜， 等. 高阶统计量方法在地球物理学中的应用与展望[J].西北地震学报， 2010， 32（2）： 201-205.
LI Xi-liang， LIU Xi-qiang， DONG Xiao-na， et al. Applicationand expectation of high-oder statistics in geophysics[J]. Northwestern Seismological Journal， 2010， 32（2）： 201-205.（In Chinese）
[10]尹成， 伍志明， 邓怀群. 高阶统计量方法在地震勘探中的应用[J]. 地球物理学进展， 2003， 18（3）： 543-550.
YIN Cheng， WU Zhi-ming， DENG Huai-qun. Application of high-order statistics in seismic exploration[J]. Progress in Geophysics， 2003， 18（3）： 543-550.（In Chinese） [11]冯智慧， 刘财， 冯t， 等. 基于互四阶累积量一维切片的地震层位自动拾取方法[J].石油地球物理勘探， 2011， 46（1）： 58-63.
FENG Zhi-hui， LIU Cai， FENG Xuan， et al. An automatic seismic horizon picking-up method based on one dimensional slice of cross fourth order cumulant[J]. Oil Geophysical Prospecting， 2011， 46（1）： 58-63.（In Chinese）
[12]熊晓军， 尹成， 张白林， 等. 高阶统计量油气检测方法研究[J]. 地球物理学报， 2004， 47（5）： 920-927.
XIONG Xiao-jun， YIN Cheng， ZHANG Bai-lin， et al. Method research of forecasting oil and gas using higher-order statistics[J]. Chinese Journal of Geophyics， 2004， 47（5）：920-927.（In Chinese）
[13]唐斌， 尹成. 基于高阶统计的非最小相位地震子波恢复[J]. 地球物理学报， 2011， 54（3）： 404-410.
TANG Bin， YIN Cheng. Non-minimum phase seismic wavelet reconstruction based on higher order statistics[J].Chinese Journal of Geophyics， 2011， 54（3）： 404-410.（In Chinese）
[14]宋维琪， 杨勤勇， 郭全仕，等. 地面微地震资料弱信号提取方法研究[J]. 石油物探， 2013， 52（2）： 131-135.
SONG Wei-qi， YANG Qin-yong， GUO Quan-shi， et al. Weak signal extraction method for surface microseismic monitoring data[J]. Geophysical Prospecting for Petroleum， 2013， 52（2）： 131-135.（In Chinese）
[15]张贤达. 时间序列分析―高阶统计量方法[M]. 北京： 清华大学出版社， 1996： 328-335.
ZHANG Xian-da. Time sequence analysis[M].Beijing ： Tsinghua University Press ， 1996： 328-335.（In Chinese）
[16]张贤达. 现代信号处理[M]. 北京： 清华大学出版社， 2002： 1-20，407-410.
ZHANG Xian-da. Modern signal processing[M]. Beijing： Tsinghua University Press， 2002： 1-20，407-410.（In Chinese）
[17]蒋忠进， 林君， 陈祖斌. 三阶累积量在可控震源勘探信号处理中的应用[J]. 系统工程与电子技术， 2003， 25（11）： 1431-1434.
JIANG Zhong-jin， LIN Jun， CHEN Zu-bin. The application of three-order cumulant in the processing of vibrator data[J].Systems Engineering and Electronics， 2003， 25（11）： 1431-1434.
[18]NIKIAS C L. Time delay estimation in unknown gaussian spatially correlated noise[J]. IEEE Transactions on Acoustics Speech and Signal Processing， 1988， 36（11）： 1706-1714.
[19]柯宏发， 陈永光， 张哲峰. 基于四阶累积量对角切片的短波自适应通信信号检测[J]. 电子学报， 2006， 34（3）： 419-423.
KE Hong-fa， CHEN Yong-guang， ZHANG Zhe-feng. HF adaptive radio communication signal detection based on fourth-order cumulants diagonal slice[J]. Acta Electronica Sinica， 2006， 34（3）： 419-423.（In Chinese）
[20]GROOS J C， RITTER J R R. Time domain classification and quantification of seismic noise in an urban environment[J]. Geophysical Journal International， 2009，179（2）：1213-1231.