基于稀疏分布特征的井下微地震信号识别与提取方法

井下微震监测获得的地震记录往往包含大量的噪声,记录信噪比很低.有效地震信号的识别与提取是进行后续地震定位等工作之前需要优先解决的问题.经过研究发现,井下水压裂微地震信号具有稀疏分布的特征,而井下环境噪声则具有更多的Gaussian分布特征.为此,本文提出将图像处理领域适宜于稀疏分布信号降噪处理的稀疏码收缩方法应用于井下微震监测数据处理.为解决需要利用与待处理数据中有效信号成分具有相似分布特征的无噪信号序列估算正交基以及计算效率等问题,将原方法与小波变换理论相结合.即通过优选小波基函数作为正交基进行小波变换将信号分解为不同级的小波系数,利用稀疏码收缩方法中对稀疏编码施加的非线性收缩方式作为阈值准则对小波系数进行改造.通过多方面的数值实验证明了该方法在处理地震子波及井下微地震信号方面准确可靠.含噪记录经过处理后有效地震信号的到时、波形、时频谱特征等均能得到良好的识别和恢复.并且该方法具有很强的抗噪能力,当信噪比低至-20~-30db时,仍然能够发挥作用.在处理大量实际井下微震监测数据的过程中,面对多种复杂情况,本方法展现出了计算效率高、计算结果可靠、应用简单等优势,证明了其本身具有实际应用价值,值得进一步的研究和推广.     

Automatic event detection in low SNR microseismic signals based on multi-scale permutation entropy and a support vector machine

Microseismic monitoring is an effective means for providing early warning of rock or coal dynamical disasters, and its first step is microseismic event detection, although low SNR microseismic signals often cannot effectively be detected by routine methods. To solve this problem, this paper presents permutation entropy and a support vector machine to detect low SNR microseismic events. First, an extraction method of signal features based on multi-scale permutation entropy is proposed by studying the influence of the scale factor on the signal permutation entropy. Second, the detection model of low SNR microseismic events based on the least squares support vector machine is built by performing a multi-scale permutation entropy calculation for the collected vibration signals, constructing a feature vector set of signals. Finally, a comparative analysis of the microseismic events and noise signals in the experiment proves that the different characteristics of the two can be fully expressed by using multi-scale permutation entropy. The detection model of microseismic events combined with the support vector machine, which has the features of high classification accuracy and fast real-time algorithms, can meet the requirements of online, real-time extractions of microseismic events.     

First arrival time picking for microseismic data based on DWSW algorithm

The first arrival time picking is a crucial step in microseismic data processing. When the signal-to-noise ratio (SNR) is low, however, it is difficult to get the first arrival time accurately with traditional methods. In this paper, we propose the double-sliding-window SW (DWSW) method based on the Shapiro-Wilk (SW) test. The DWSW method is used to detect the first arrival time by making full use of the differences between background noise and effective signals in the statistical properties. Specifically speaking, we obtain the moment corresponding to the maximum as the first arrival time of microseismic data when the statistic of our method reaches its maximum. Hence, in our method, there is no need to select the threshold, which makes the algorithm more facile when the SNR of microseismic data is low. To verify the reliability of the proposed method, a series of experiments is performed on both synthetic and field microseismic data. Our method is compared with the traditional short-time and long-time average (STA/LTA) method, the Akaike information criterion, and the kurtosis method. Analysis results indicate that the accuracy rate of the proposed method is superior to that of the other three methods when the SNR is as low as ??10 dB.     

复小波包域的微震信号分析提取方法

提出了一种在复小波包域分析提取微震信号的新方法。此方法先采用复小波包分解同时获得信号的幅值信息和相位信息,再采用一种称为“比值加权”的复小波包重构法,重构出信噪比大大增强的微震信号。实际应用表明这种新方法在提取微震信号,提高信噪比方面效果显著。     

基于小波变换的大尺度岩体结构微震监测信号去噪方法研究

为将小波去噪方法应用于大尺度岩体结构微震监测信号的去噪研究,首先在MATLAB环境下进行仿真,验证了使用Symlet6小波进行小波去噪的可行性;利用4种自适应阈值规则对含噪信号进行去噪对比,结果表明4种阈值去噪后的信号在均方差较小的情况下都极大地提高了信号的信噪比,有效地去除了噪声,对不同的含噪信号,无偏似然原则阈值去...     

基于小波包和峰度赤池信息量准则的P波震相自动识别方法

基于小波包变换和峰度赤池信息量准则(AIC), 提出了一种新的自动识别P波震相的综合方法, 即小波包-峰度AIC方法． 首先对由加权长短时窗平均比(STA/LTA)法粗略确定的P波到时前后3 s的记录进行小波包三尺度的分解与重构, 分别计算每个尺度重构信号的峰度AIC曲线并将其叠加, 叠加曲线的最小值则为P波震相到时; 然后对原始地震记录进行有限冲激响应自适应滤波以提高信噪比和识别精度; 最后将小波包-峰度AIC方法应用到合成理论地震图及实际地震记录的P波初至自动识别中. 结果表明: 初至清晰度对识别精度的影响比信噪比对其影响更大; 与单独使用加权STA/LTA方法和峰度AIC法相比, 小波包-峰度AIC法具有更强的抗噪能力, 识别精度更高; 当初至清晰时, 小波包-峰度AIC法自动识别与人工识别的P波到时平均绝对差值为(0.077±0.075) s．      

Method for Identifying Micro-seismic <Emphasis Type="Italic">P</Emphasis>-Arrival by Time-frequency Analysis Using Intrinsic Time-Scale Decomposition

A method to identify the P-arrival of microseismic signals is proposed in this work, based on the algorithm of intrinsic timescale decomposition (ITD). Using the results of ITD decomposition of observed data, information of instantaneous amplitude and frequency can be determined. The improved ratio function of short-time average over long-time average and the information of instantaneous frequency are applied to the time-frequency-energy denoised signal for picking the P-arrival of the microseismic signal. We compared the proposed method with the wavelet transform method based on the denoised signal resulting from the best basis wavelet packet transform and the single-scale reconstruction of the wavelet transform. The comparison results showed that the new method is more effective and reliable for identifying P-arrivals of microseismic signals.     

A first arrival picking method of microseismic data based on single time window with window length independent

First arrival picking is a key factor which affects the precision of microseismic data analysis. Here, we propose a new method, which employs the maximum eigenvalue to constraint the Maeda-Akaike Information Criterion (Maeda-AIC) algorithm. First, aims at addressing the pick result affected by signal-to-noise ratio (SNR) of microseismic data, maximum eigenvalue method based on polarization analysis is applied, and the maximum eigenvalue is calculated firstly, as for three component (3C) microseismic data, the maximum eigenvalue is calculated with corresponding covariance matrix, a time window need to be set in the process of building the covariance matrix, and it is the only time window set in the method proposed in this paper, so the method is called single window Maeda-AIC (SWM-AIC), to the single component (1C) microseismic data, the variance of the data is taken as the maximum eigenvalue. Then, to reduce the effect of time window and increase the automation of the algorithm, Maeda-AIC method which is a non-window-based first arrival picking method is applied. Maeda-AIC values in preliminary window are calculated, and the preliminary window is the sequence before the largest eigenvalue of the 3C or 1C data. We validate the developed method with both synthetic and field microseismic data, using a range of signal-to-noise ratios. The developed method is compared with some basic methods, specifically STA/LTA, Maeda-AIC, and the maximum eigenvalue method. The results demonstrate that the new method is much better at identifying first arrival times than basic methods when the data have a low signal-to-noise ratio, and is even faster than the STA/LTA method with 1C data. In contrast to other improved methods, threshold value is not required for this method, and the only time window used in this method is just for maximum eigenvalue calculation, through test in the paper, its length has almost no effect on the first arrival picking.     

基于改进长短时窗比值及优化变分模态分解的微震初至拾取算法

针对低信噪比条件下微震初至拾取准确度低的问题,基于信号幅度变化引入权重因子,对传统长短时窗比值（STA/LTA）算法进行改进,提高初次拾取精度。为了进一步降低拾取误差,对变分模态分解（VMD）算法进行优化,基于互相关系数和排列熵准则自适应确定VMD分解层数,对初次拾取结果前后2—3 s的记录进行优化VMD,并计算分解后各本征模函数（IMF）的峰度赤池信息准则值,得到各IMF的到时,以各IMF的拾取结果及能量比综合加权得到二次拾取到时。仿真实验表明:改进后的STA/LTA在较低信噪比下可降低初次拾取误差约0.01 s以上;相比经验模态分解（EMD)和小波包分解,自适应VMD分解后能再次降低误差,最终与人工拾取结果平均误差在0.023 s以内。实际微震信号初至拾取结果表明,本算法能快速有效地识别初至P波,与人工拾取结果相比误差小,准确率高。      

时窗能量特征法拾取微地震波初始到时的可行性研究

微地震（MS）波初始到时的自动拾取是MS监测数据处理的关键技术之一,也是实现MS震源自动定位的技术难点.本文在MS震源定位结果反演与推断的研究基础上,对不同类型MS波的到时点特征进行了分析与描述,并对不同时窗长度下能量特征值的变化规律进行了研究,提出了控制时窗移动范围和确定时窗长度自适应参数的具体方法,利用建立的MS波初始到时点特征的模式识别库,对拾取的到时进行模式归类、定量评价和匹配,提高了自动拾取结果的可靠性.研究结果表明,对典型的信噪比高的MS波,到时自动拾取的结果与手工拾取的结果基本一致;对无量纲大振幅的MS波,到时自动拾取结果的可靠性要高于手工拾取,对信噪比低和到时点不清晰的MS波自动拾取的可靠性较低.     

基于时频稀疏性分析法的低信噪比微震事件识别与恢复

微震监测是直观评价压裂过程和压裂效果的有效手段.微震事件识别是微震监测的首要步骤.然而对于低信噪比微震监测数据,常规识别方法很难取得满意效果.基于微震事件在时频域中的稀疏性,本文提出利用Renyi熵值表示微震监测数据的时频稀疏程度,并以时频距离为约束条件,建立以低熵值的道数为判别阈值的目标函数.本文方法能在识别出微震事件的同时,恢复出较为清晰的微震事件.通过数值计算和对实际监测数据的测试,表明该方法对低信噪比的微震监测数据有较好的处理效果.     

基于Context模型的Shearlet变换地面微地震随机噪声压制

地面微地震监测采集到的微地震信号通常能量微弱,信噪比低,如何提高微震数据的信噪比是数据处理的难题.Shearlet变换是一种新型的多尺度几何分析方法,具有敏感的方向性和较强的稀疏表示特性,能起到很好的随机噪声压制效果.由于地面微震数据的有效信号大多被淹没在噪声中,基于传统阈值的Shearlet变换（the traditional threshold-based Shearlet transform TST）只考虑到尺度或方向的阈值,在去噪过程中会过度扼制有效信号系数,造成有效信号能量损失.因而,本文建立Context模型,得到基于Context模型的Shearlet变换（the Context-model-based Shearlet transform CMST）方法,改进传统Shearlet阈值方法的不足.我们通过所建立的Context模型将能量相近的各方向系数划分为同一组,并分组估计阈值,分别处理各部分系数,达到微弱同相轴有效恢复的目的.通过TST及CMST的模拟实验与实际地面微震记录处理结果对比可知,本文方法在低信噪比条件下比对比方法更加有效地恢复地面微震数据的微弱信号,随机噪声压制效果明显,在-10 dB条件下,提升信噪比18.3741 dB.     

微地震事件初至拾取SLPEA算法

微地震事件初至拾取是微地震数据处理的关键步骤之一.实际微地震监测资料中存在大量低信噪比事件,而传统方法对这些事件的应用效果并不理想.为了克服传统方法抗噪性弱的缺点,本文通过综合地震信号与环境噪声在振幅、偏振以及统计特征等方面的存在的差异,设计了一种针对低信噪比微地震事件的初至拾取方法——SLPEA算法.为了检验本文方法的可行性和有效性,分别对模型数据和实际资料进行了处理,并将处理结果与传统方法及手工拾取的结果进行了对比.分析表明,利用本文方法得到的初至到时与手工拾取结果的绝对误差平均值仅为1.33×10~(-3)s,小于3个采样点;方差为3.21×10~(-6)s~2;初至到时在手工拾取结果±0.005s误差范围内的个数占总数的95.8%.这些参数值均优于传统方法的同类参数,证明了本文方法的可靠性.     

Weak Seismic Signal Extraction Based on the Curvelet Transform

Seismic signal denoising is a key step in seismic data processing. Airgun signals are easy to be interfered with by noise when it travels a long distance due to the weak energy of active source signal of the airgun. Aiming to solve this problem, and considering that the conventional Curvelet transform threshold processing method does not use the seismic spectrum information, we independently process the Curvelet scale layer corresponding to valid data based on the characteristics of the Curvelet transform of multi-scale, multi-direction and capable of expressing the sparse seismic signals in order to fully excavate the information features. Combined with the Curvelet adaptive threshold denoising the algorithm, we apply the Curvelet transform to denoising seismic signals while retaining the weak information in the signal as much as possible. The simulation experiments show that the improved threshold denoising method based on Curvelet transform is superior to the frequency domain filtering, wavelet denoising and traditional Curvelet denoising method in detailed information extraction and signal denoising of low SNR signals. The calculation accuracy of the relative wave velocity variation of underground medium is improved.     

Identification of blasting vibration and coal-rock fracturing microseismic signals

Α new method based on variational mode decomposition (VMD) is proposed to distinguish between coal-rock fracturing and blasting vibration microseismic signals. First, the signals are decomposed to obtain the variational mode components, which are ranked by frequency in descending order. Second, each mode component is extracted to form the eigenvector of the energy of the original signal and calculate the center of gravity coefficient of the energy distribution plane. Finally, the coal-rock fracturing and blasting vibration signals are classified using a decision tree stump. Experimental results suggest that VMD can effectively separate the signal components into coal-rock fracturing and blasting vibration signals based on frequency. The contrast in the energy distribution center coefficient after the dimension reduction of the energy distribution eigenvector accurately identifies the two types of microseismic signals. The method is verified by comparing it to EMD and wavelet packet decomposition.     

HOMOMORPHIC FILTERING — THEORY AND PRACTICE*

The application of homomorphic filtering in marine seismic reflection work is investigated with the aims to achieve the estimation of the basic wavelet, the wavelet deconvolution and the elimination of multiples. Each of these deconvolution problems can be subdivided into two parts: The first problem is the detection of those parts in the cepstrum which ought to be suppressed in processing. The second part includes the actual filtering process and the problem of minimizing the random noise which generally is enhanced during the homomorphic procedure. The application of homomorphic filters to synthetic seismograms and air-gun measurements shows the possibilities for the practical application of the method as well as the critical parameters which determine the quality of the results. These parameters are:

• a) the signal-to-noise ratio (SNR) of the input data
• b) the window width and the cepstrum components for the separation of the individual parts
• c) the time invariance of the signal in the trace.

In the presence of random noise the power cepstrum is most efficient for the detection of wavelet arrival times. For wavelet estimation, overlapping signals can be detected with the power cepstrum up to a SNR of three. In comparison with this, the detection of long period multiples is much more complicated. While the exact determination of the water reverberation arrival times can be realized with the power cepstrum up to a multiples-to-primaries ratio of three to five, the detection of the internal multiples is generally not possible, since for these multiples this threshold value of detectibility and arrival time determination is generally not realized. For wavelet estimation, comb filtering of the complex cepstrum is most valuable. The wavelet estimation gives no problems up to a SNR of ten. Even in the presence of larger noise a reasonable estimation can be obtained up to a SNR of five by filtering the phase spectrum during the computation of the complex cepstrum. In contrast to this, the successful application of the method for the multiple reduction is confined to a SNR of ten, since the filtering of the phase spectrum for noise reduction cannot be applied. Even if the threshold results are empirical, they show the limits fór the successful application of the method.     

基于重构激发信号的电磁式可控震源相关检测参考信号方法研究

针对电磁式可控震源地震数据的相关检测,研究发现,在地下结构复杂、基板-大地耦合不佳时,常规方法——基于震源控制信号或基板附近信号作为参考信号检测得到的地震记录中,存在子波到时误差和虚假多次波问题.本文分析了上述问题的理论原因,并提出基于重构激发信号的相关检测参考信号方法(Correlation Detection Reference Signal Based on the Reconstructed Excitation Signal,CDRSBRES).首先,利用直达波与其他地震波到时不一致的特点,从震源基板附近信号中分离、提取直达波.然后,利用直达波重构震源激发信号并作为参考信号对地震数据进行相关检测.最后,应用谱白化技术提高检测结果质量.数值模拟研究表明,重构激发信号与理想激发信号的相关系数为0.9869,达到高度线性相关,CDRSBRES方法检测的地震记录在子波到时和波形特征上均与模型相符.随后,在某金属矿区开展了可控震源对比实验.与液压式可控震源MiniVib T15000检测结果相比,电磁式可控震源PHVS 500的检测结果中:基于震源控制信号的检测结果存在子波到时误差约0.012s,对应垂向精度误差约11.16m;基于基板附近信号的检测结果部分区域出现虚假多次波,信噪比降低;而CDRSBRES方法的检测结果子波到时误差约0.001s,对应垂向精度误差约0.93m,波形特征一致,相同区域无虚假多次波.综上,本方法适用于电磁式可控震源地震数据的高精度检测,尤其对于地下结构复杂区域的高分辨率地震勘探具有重要意义.     

First-break picking for microseismic data based on cascading use of Shearlet and Stockwell transforms

First-break picking of microseismic data is a significant step in microseismic monitoring. There is a great error in conventional first-break picking methods based on time domain analysis in low signal to noise ratio. S-transform may provide a novel approach, it can extract the time–frequency features of the signal and reduce the picking error because of its high time–frequency resolution and good time–frequency clustering; however, the S-transform is not well suited for microseismic data with high noise. For applications to array data where the weak signal has spatial coherency as well as some distinct temporal characteristics, we propose to combine the shearlet transform with a time–frequency transform. In the proposed method, the shearlet transform is used to capture spatial coherency features of the signal. The information of the signal and noise in shearlet domain is represented by shearlet coefficients. We use the correlation of signal coefficients at adjacent fine scales to give prominence to signal features to accurately discriminate the signal from noise. The prominent signal coefficients make the signal better gathered in time–frequency spectrum of the S-transform. Finally, we can get reliable and accurate first breaks based on the change of energy. The performance of the proposed method was tested on synthetic and field microseismic data. The experimental results indicated that our method is outstanding in terms of both picking precision and adaptability to noise.     

煤矿矿震定位中异向波速模型的构建与求解

针对煤矿上覆岩层层状赋存和离层带的特点,构建矿井尺度的微震监测系统异向波速模型,模型中波速向量由地面探头速度与井下探头速度组成.研究了在只有强矿震信号和混有爆破信号两种条件下,以到时残差最小为目标和震源定位误差最小为目标的两种求解模型,模型求解选用具有全局寻优特性的遗传算法与CMEAS算法结合的混合算法.现场实际应用得出,只使用爆破信号的到时残差法最优,混有强矿震信号的到时残差法其次;与爆破信号定位所用的统一简化波速模型相比,震源定位误差大幅度降低.在此基础上进一步减低定位误差,还需从微震台网的优化布设方面解决.     

二氧化碳相变技术应用于新型震源研发的可行性研究

为了研究二氧化碳物理相变技术应用于新型震源研发的可行性,在地下成层性较好的某煤田地震测区,开展了利用二氧化碳相变技术激发地震波的野外人工震源激发-接收实验．并与传统炸药震源进行了对比．地震数据利用Aries2.66型垂直分量反射地震仪和PDS-2型三分量地震仪接收．根据实测地震数据,从野外地震记录震相识别,初至波传播距离分析,震源近场地震信号时频分析,CO_2相变激发震源子波提取和基于CO_2震源子波的地震初至波波形反演实验等多个方面,进行了关于CO_2相变激发技术能否产生地震波信号以及能否将其应用于新型震源研发的可行性研究．研究结果表明CO_2物理相变膨胀能够产生能量集中的地震波信号;在实验区地质条件和激发参量下地震记录中初至波的可识别的传播距离约为1km;震源近场地震信号的主频集中在8～13Hz;利用震源近场数据提取了CO_2震源子波;通过地震初至波波形反演实验认为这种震源子波能够应用于波形反演等方面的研究．因为CO_2相变激发具有绿色、环保、安全等方面的优点,若能进一步在激发能量、激发—延迟时间一致性等方面加以改进,该技术有望在城市隐伏活动断层探测、城市地下空间探测、煤矿高瓦斯环境人工地震勘探等领域发挥重要的作用．