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

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

A first arrival detection method for low SNR microseismic signal

Most of the microseismic signals have low signal-to-noise ratio (SNR) due to the strong background noise, which makes it difficult to locate the first arrival time. Both accuracy and stability of conventional methods are poor in this situation. To overcome this problem, here we proposed a new method based on the adaptive Morlet wavelet and principal component analysis process in wavelet coefficients matrix. The three components of microseismic signal make it possible to extract the features in wavelet coefficients domain. Then the reconstructed signal from weighted features presents an obvious first arrival. Tests on synthetic signals and real data provide a solid evidence for its feasibility in low SNR microseismic signal.     

Full‐waveform based microseismic event detection and signal enhancement: An application of the subspace approach

Microseismic monitoring has proven invaluable for optimizing hydraulic fracturing stimulations and monitoring reservoir changes. The signal to noise ratio of the recorded microseismic data varies enormously from one dataset to another, and it can often be very low, especially for surface monitoring scenarios. Moreover, the data are often contaminated by correlated noises such as borehole waves in the downhole monitoring case. These issues pose a significant challenge for microseismic event detection. In addition, for downhole monitoring, the location of microseismic events relies on the accurate polarization analysis of the often weak P‐wave to determine the event azimuth. Therefore, enhancing the microseismic signal, especially the low signal to noise ratio P‐wave data, has become an important task. In this study, a statistical approach based on the binary hypothesis test is developed to detect the weak events embedded in high noise. The method constructs a vector space, known as the signal subspace, from previously detected events to represent similar, yet significantly variable microseismic signals from specific source regions. Empirical procedures are presented for building the signal subspace from clusters of events. The distribution of the detection statistics is analysed to determine the parameters of the subspace detector including the signal subspace dimension and detection threshold. The effect of correlated noise is corrected in the statistical analysis. The subspace design and detection approach is illustrated on a dual‐array hydrofracture monitoring dataset. The comparison between the subspace approach, array correlation method, and array short‐time average/long‐time average detector is performed on the data from the far monitoring well. It is shown that, at the same expected false alarm rate, the subspace detector gives fewer false alarms than the array short‐time average/long‐time average detector and more event detections than the array correlation detector. The additionally detected events from the subspace detector are further validated using the data from the nearby monitoring well. The comparison demonstrates the potential benefit of using the subspace approach to improve the microseismic viewing distance. Following event detection, a novel method based on subspace projection is proposed to enhance weak microseismic signals. Examples on field data are presented, indicating the effectiveness of this subspace‐projection‐based signal enhancement procedure.     

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

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

基于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.     

基于超宽带雷达及支持向量机的灾后人体呼吸信号识别方法与试验研究

与常规雷达相比,超宽带雷达具有距离分辨力高、近距离盲区小、穿透性强、目标识别率高等特点,已被广泛应用于灾后搜寻、救援工作中,以对受困生命体征目标进行生命探测。为实现使用超宽带雷达对受困生命体征目标的识别定位,本研究提出基于信号多特征提取技术及支持向量机模型的人体呼吸信号识别方法。首先,使用经验模态分解、变分模态分解及希尔伯特变换提取雷达探测信号的微多普勒特征,使用傅里叶变换提取宏观频谱特征,使用相关分析获取相关性特征;然后,以提取的信号特征为输入,使用支持向量机模型对信号进行分类,进而对人体呼吸信号进行识别,对人体位置进行定位。不同障碍物场景下的试验结果表明,本方法可有效识别砖墙、建筑楼板等遮挡物下的受困生命体征目标,并提供其位置信息。     

Convolutional neural networks for automated microseismic detection in downhole distributed acoustic sensing data and comparison to a surface geophone array

Distributed acoustic sensing is a growing technology that enables affordable downhole recording of strain wavefields from microseismic events with spatial sampling down to ∼1 m. Exploiting this high spatial information density motivates different detection approaches than typically used for downhole geophones. A new machine learning method using convolutional neural networks is described that operates on the full strain wavefield. The method is tested using data recorded in a horizontal observation well during hydraulic fracturing in the Eagle Ford Shale, Texas, and the results are compared to a surface geophone array that simultaneously recorded microseismic activity. The neural network was trained using synthetic microseismic events injected into real ambient noise, and it was applied to detect events in the remaining data. There were 535 detections found and no false positives. In general, the signal-to-noise ratio of events recorded by distributed acoustic sensing was lower than the surface array and 368 of 933 surface array events were found. Despite this, 167 new events were found in distributed acoustic sensing data that had no detected counterpart in the surface array. These differences can be attributed to the different detection threshold that depends on both magnitude and distance to the optical fibre. As distributed acoustic sensing data quality continues to improve, neural networks offer many advantages for automated, real-time microseismic event detection, including low computational cost, minimal data pre-processing, low false trigger rates and continuous performance improvement as more training data are acquired.     

Microseismic event denoising via adaptive directional vector median filters

We present a novel denoising scheme via Radon transform-based adaptive vector directional median filters named adaptive directional vector median filter (AD-VMF) to suppress noise for microseismic downhole dataset. AD-VMF contains three major steps for microseismic downhole data processing: (i) applying Radon transform on the microseismic data to obtain the parameters of the waves, (ii) performing S-transform to determine the parameters for filters, and (iii) applying the parameters for vector median filter (VMF) to denoise the data. The steps (i) and (ii) can realize the automatic direction detection. The proposed algorithm is tested with synthetic and field datasets that were recorded with a vertical array of receivers. The P-wave and S-wave direct arrivals are properly denoised for poor signal-to-noise ratio (SNR) records. In the simulation case, we also evaluate the performance with mean square error (MSE) in terms of signal-to-noise ratio (SNR). The result shows that the distortion of the proposed method is very low; the SNR is even less than 0 dB.     

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.     

基于混沌算法的地震信息网络入侵检测研究

为改善传统的基于机器学习的网络入侵检测方法只能检测已知入侵行为,对于未知入侵行为的检测存在误警率高、时效性差的不足,提出一种基于混沌算法的地震信息网络入侵检测方法。创建候选地震信息网络特征-混沌变量映射模型,实现变量之间的转化;采用混沌变量迭代演化算法进行地震信息网络特征选择;使用支持向量机对最优特征进行学习,为提高地震信息网络入侵检测精度,利用柯西蜂群算法对支持向量机参数进行寻优,建立网络入侵检测优化模型。仿真实验证明,基于混沌算法的地震信息网络入侵检测方法能有效实现高检测率、低误报率的入侵检测,具有很高的应用优势。     

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.     

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

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

基于多尺度卷积神经网络的结构损伤识别研究

为了利用结构振动响应的时间多尺度特征来提升卷积神经网络识别结构损伤的能力,给出了两种用于结构损伤识别的多尺度卷积神经网络,即多尺度输入和多尺度卷积核卷积神经网络.对于多尺度输入卷积神经网络,将通过下采样和滑动平均获取的具有不同时间尺度特征的振动信号输入固定尺寸卷积核的分支卷积神经网络;对于多尺度卷积核卷积神经网络,则将...     

Automatic signal detection based on support vector machine

Introduction The automatic processing of continuous seismic data is important for monitoring earthquake, in which real data recorded by field stations located in different regions is transmitted to data cen- tre through internet or satellite communication systems. Automatic processing will run firstly on data, afterwards these automatic processing results will be reviewed and modified. The load of interactive analysis would be increase if there were more false events or missed events after run…     

Comparison of migration‐based location and detection methods for microseismic events

Microseismic monitoring in the oil and gas industry commonly uses migration‐based methods to locate very weak microseismic events. The objective of this study is to compare the most popular migration‐based methods on a synthetic dataset that simulates a strike‐slip source mechanism event with a low signal‐to‐noise ratio recorded by surface receivers (vertical components). The results show the significance of accounting for the known source mechanism in the event detection and location procedures. For detection and location without such a correction, the ability to detect weak events is reduced. We show both numerically and theoretically that neglecting the source mechanism by using only absolute values of the amplitudes reduces noise suppression during stacking and, consequently, limits the possibility to retrieve weak microseismic events. On the other hand, even a simple correction to the data polarization used with otherwise ineffective methods can significantly improve detections and locations. A simple stacking of the data with a polarization correction provided clear event detection and location, but even better results were obtained for those data combined with methods that are based on semblance and cross‐correlation.     

A robust method for microseismic event detection based on automatic phase pickers

We present a robust method for the automatic detection and picking of microseismic events that consists of two steps. The first step provides accurate single-trace picks using three automatic phase pickers adapted from earthquake seismology. In the second step, a multi-channel strategy is implemented to associate (or not) the previous picks with actual microseismic signals by taking into account their expected alignment in all the available channels, thus reducing the false positive rate. As a result, the method provides the number of declared microseismic events, a confidence indicator associated with each of them, and the corresponding traveltime picks. Results using two field noisy data records demonstrate that the automatic detection and picking of microseismic events can be carried out with a relatively high confidence level and accuracy.     

基于改进EWT和LogitBoost集成分类器的地震事件分类识别算法

快速、准确地识别天然地震和人工爆破事件是地震台网监测的重要工作之一,也是提高地震观测记录质量、开展地震研究工作的重要基础.针对反向传播神经网络、支持向量机等主流分类识别方法在地震事件分类识别应用上的不足,提出一种基于改进 EWT 和 LogitBoost集成分类器的地震事件分类识别算法.首先,基于S谱能量曲线对传统经验小波变换进行改进,将信号自适应分解为按频率和能量分布的本征模函数;其次,提取 P波与S波最大振幅比,前4个本征模函数的香农熵、对数能量熵,以及去噪后重构信号主频等特征;最后,采用基于集成学习 LogitBoost的决策树集成分类器进行分类.实验结果表明,所提算法具有较高的鲁棒性,能有效解决样本不足的问题,识别准确率达93．1%以上,比集成学习 AdaBoost、反向传播神经网络和支持向量机等方法提高了1%以上,且分类识别效果好.     

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

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

Application of Kalman Filtering Techniques for Microseismic Event Detection

—?Microseismic monitoring systems are generally installed in areas of induced seismicity caused by human activity. Induced seismicity results from changes in the state of stress which may occur as a result of excavation within the rock mass in mining (i.e., rockbursts), and changes in hydrostatic pressures and rock temperatures (e.g., during fluid injection or extraction) in oil exploitation, dam construction or fluid disposal. Microseismic monitoring systems determine event locations and important source parameters such as attenuation, seismic moment, source radius, static stress drop, peak particle velocity and seismic energy. An essential part of the operation of a microseismic monitoring system is the reliable detection of microseismic events. In the absence of reliable, automated picking techniques, operators rely upon manual picking. This is time-consuming, costly and, in the presence of background noise, very prone to error. The techniques described in this paper not only permit the reliable identification of events in cluttered signal environments they have also enabled the authors to develop reliable automated event picking procedures. This opens the way to use microseismic monitoring as a cost-effective production/operations procedure. It has been the experience of the authors that in certain noisy environments, the seismic monitoring system may trigger on and subsequently acquire substantial quantities of erroneous data, due to the high energy content of the ambient noise. Digital filtering techniques need to be applied on the microseismic data so that the ambient noise is removed and event detection simplified. The monitoring of seismic acoustic emissions is a continuous, real-time process and it is desirable to implement digital filters which can also be designed in the time domain and in real-time such as the Kalman Filter. This paper presents a real-time Kalman Filter which removes the statistically describable background noise from the recorded seismic traces.     

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.