Integrated processing method for microseismic signal based on deep neural network

Denoising and onset time picking of signals are essential before extracting source information from collected seismic/microseismic data. We proposed an advanced deep dual-tasking network (DDTN) that integrates these two procedures sequentially to achieve the optimal performance. Two homo-structured encoder–decoder networks with specially designed structures and parameters are connected in series for handling the denoising and detection of microseismic signals. Based on the similarity of data types, the output of the denoising network will be imported into the detection network to obtain labels for the signal duration. The procedures of denoising and duration detection can be completed in an integrated way, where the denoised signals can improve the accuracy of onset time picking. Results show that the method has a good performance for the denoising of microseismic signals that contain various types and intensities of noise. Compared with existing methods, DDTN removes the noise with a minor waveform distortion. It is ideal for recovering the microseismic signal while maintaining a good capacity for onset time picking when the signal-to-noise ratio is low. Based on that, the second network can detect a more accurate duration of microseismic signals and thus obtain more accurate onset time. The method has great potential to be extended to the study of exploration seismology and earthquakes

MSSPN: Automatic First Arrival Picking using Multi-Stage Segmentation Picking Network

Picking the first arrival times of prestack gathers is called First Arrival Time (FAT) picking, which is an indispensable step in seismic data processing, and is mainly solved manually in the past. With the current increasing density of seismic data collection, the efficiency of manual picking has been unable to meet the actual needs. Therefore, automatic picking methods have been greatly developed in recent decades, especially those based on deep learning. However, few of the current supervised deep learning-based method can avoid the dependence on labeled samples. Besides, since the gather data is a set of signals which are greatly different from the natural images, it is difficult for the current method to solve the FAT picking problem in case of a low Signal to Noise Ratio (SNR). In this paper, for hard rock seismic gather data, we propose a Multi-Stage Segmentation Pickup Network (MSSPN), which solves the generalization problem across worksites and the picking problem in the case of low SNR. In MSSPN, there are four sub-models to simulate the manually picking processing, which is assumed to four stages from coarse to fine. Experiments on seven field datasets with different qualities show that our MSSPN outperforms benchmarks by a large margin.Particularly, our method can achieve more than 90\% accurate picking across worksites in the case of medium and high SNRs, and even fine-tuned model can achieve 88\% accurate picking of the dataset with low SNR

Earthquake Arrival Association with Backprojection and Graph Theory

The association of seismic wave arrivals with causative earthquakes becomes progressively more challenging as arrival detection methods become more sensitive, and particularly when earthquake rates are high. For instance, seismic waves arriving across a monitoring network from several sources may overlap in time, false arrivals may be detected, and some arrivals may be of unknown phase (e.g., P- or S-waves). We propose an automated method to associate arrivals with earthquake sources and obtain source locations applicable to such situations. To do so we use a pattern detection metric based on the principle of backprojection to reveal candidate sources, followed by graph-theory-based clustering and an integer linear optimization routine to associate arrivals with the minimum number of sources necessary to explain the data. This method solves for all sources and phase assignments simultaneously, rather than in a sequential greedy procedure as is common in other association routines. We demonstrate our method on both synthetic and real data from the Integrated Plate Boundary Observatory Chile (IPOC) seismic network of northern Chile. For the synthetic tests we report results for cases with varying complexity, including rates of 500 earthquakes/day and 500 false arrivals/station/day, for which we measure true positive detection accuracy of > 95%. For the real data we develop a new catalog between January 1, 2010 - December 31, 2017 containing 817,548 earthquakes, with detection rates on average 279 earthquakes/day, and a magnitude-of-completion of ~M1.8. A subset of detections are identified as sources related to quarry and industrial site activity, and we also detect thousands of foreshocks and aftershocks of the April 1, 2014 Mw 8.2 Iquique earthquake. During the highest rates of aftershock activity, > 600 earthquakes/day are detected in the vicinity of the Iquique earthquake rupture zone

Joint detection and classification of rockfalls in a microseismic monitoring network

A rockfall (RF) is a ubiquitous geohazard that is difficult to monitor or predict and poses a significant risk for people and transportation in several hilly and mountainous environments. The seismic signal generated by RF carries abundant physical and mechanical information. Thus, signals can be used by researchers to reconstruct the event location, onset time, volume and trajectory, and develop an efficient early warning system. Therefore, the precise automatic detection and classification of RF events are important objectives for scientists, especially in seismic monitoring arrays. An algorithm called DESTRO (DEtection and STorage of ROckfalls) aimed at combining seismic event automatic detection and classification was implemented ad hoc within the MATLAB environment. In event detection, the STA/LTA (short-time-average through long-time-average) method combined with other parameters, such as the minimum duration of an RF and the minimum interval time between two continuous seismic events is used. Furthermore, nine significant features based on the frequency, amplitude, seismic waveform, duration and multiple station attributes are newly proposed to classify seismic events in a RF environment. In particular, a three-step classification method is proposed for the discrimination of five different source types: RFs, earthquakes (EQs), tremors, multispike events (MSs) and subordinate MS events. Each component (vertical, east–west and north–south) at each station within the monitoring network is analysed, and a three-step classification is performed. At a given time, the event series detected from each component are integrated and reclassified component by component and station by station into a final event-type series as an output result. By this algorithm, a case study of the seven-month-long seismic monitoring of a former quarry in Central Italy was investigated by means of four triaxial velocimeters with continuous acquisition at a sampling rate of 200 Hz. During this monitoring period, a human-induced RF simulation was performed, releasing 95 blocks (in which 90 blocks validated) of different sizes from the benches of the quarry. Consequently, 64.9 per cent of EQs within 100 km were confirmed in a one-month monitoring period, 88 blocks in the RF simulation were classified correctly as RF events and 2 blocks were classified as MSs given their small energy. Finally, an ad hoc section of the algorithm was designed specifically for RF classification combined with EQ recognition. The algorithm could be applied in slope seismic monitoring to monitor the dynamic states of rock masses, as well as in slope instability forecasting and risk evaluation in EQ-prone areas

Use of the shearlet energy entropy and of the support vector machine classifier to process weak microseismic and desert seismic signals

Low-amplitude signal detection is a key procedure in borehole microseismic and desert seismic exploration. Usually, signals are difficult to detect due to their low amplitude and noise contamination. To solve this problem, we propose a method combining shearlet energy entropy with a support vector machine (SVM) to detect low-amplitude signals. In the proposed method, the signal feature is extracted using shearlet energy entropy. The signal is more sparsely represented in the shearlet domain because of the multi-scale and multi-direction characteristic of the shearlet transform, which favours signal feature extraction. Furthermore, in calculating shearlet energy entropy, we use the correlation of shearlet coefficients to enhance the difference between signal and noise in the shearlet domain. Shearlet energy entropy makes the SVM achieve a more accurate classification result compared with other traditional features such as amplitude and energy. The results of synthetic and field data show that our method is more effective than the STA/LTA and the convolutional neural network for low-amplitude microseismic signal and desert seismic signal detection

Use of the shearlet energy entropy and of the support vector machine classifier to process weak microseismic and desert seismic signals

Low-amplitude signal detection is a key procedure in borehole microseismic and desert seismic exploration. Usually, signals are difficult to detect due to their low amplitude and noise contamination. To solve this problem, we propose a method combining shearlet energy entropy with a support vector machine (SVM) to detect low-amplitude signals. In the proposed method, the signal feature is extracted using shearlet energy entropy. The signal is more sparsely represented in the shearlet domain because of the multi-scale and multi-direction characteristic of the shearlet transform, which favours signal feature extraction. Furthermore, in calculating shearlet energy entropy, we use the correlation of shearlet coefficients to enhance the difference between signal and noise in the shearlet domain. Shearlet energy entropy makes the SVM achieve a more accurate classification result compared with other traditional features such as amplitude and energy. The results of synthetic and field data show that our method is more effective than the STA/LTA and the convolutional neural network for low-amplitude microseismic signal and desert seismic signal detection