The geophone can daily detect a large amount of microseismic signals and it is time consuming for manual picking the arrival times of P-S waves. We propose an automatic P-and S-wave onset picking method based on single-component geophone, which is widely used in mines. Firstly, we filter the original signal with a Butterworth bandpass filter in order to increase the SNR. The mine microseismic signal frequency ranges from 10 Hz to 2000 Hz, so we choose the band-pass filter interval to be 3-2400 Hz. Secondly, we need to locate the P-S wave phases in the signal. An improved characteristic function can strengthen the change of frequency and amplitude. We define two neighboring windows of different sizes, with the short time window lying in the end of the long time window. Then we compute the average value of the characteristic function within the window, which is also called STA/LTA (short-term average value/long-term average value). With the sliding of the long-short time window, we can get the positions of the P-S wave onset windows. The ratio of STA/LTA will have two peaks because the arrivals of P wave and S wave. Through a preset threshold, we can locate two time windows including P wave onset time and S wave onset time, respectively. In the last step, Akaike's information criterion(AR-AIC criterion) is used in the two time windows, the global minimum of AIC in the window is the arriving time. We have tested this method with one thousand single-component microseismic signals, the result has shown that 71% of P-wave arrival time picking is correct and 79% of S-wave picking is correct.
WANG Liguan
,
WANG Hui
,
ZHAO Junjie
. An automatic P-and S-wave onset-picking designed for single-component microseismic signals[J]. Science & Technology Review, 2016
, 34(2)
: 184
-189
.
DOI: 10.3981/j.issn.1000-7857.2016.2.030
[1] 刘劲松, 王赟, 姚振兴. 微地震信号到时自动拾取方法[J]. 地球物理学报, 2013, 56(5): 1660-1666. Liu Jingsong, Wang Yun, Yao Zhenxing. On micro-seismic first arrival identification: A case study[J]. Chinese Journal of Geohpysics, 2013, 56 (5): 1660-1666.
[2] Allen R V. Automatic earthquake recognition and timing from single trac-es[J]. Bulletin of the Seismological Society of America, 1978, 68(5): 1521-1532.
[3] Boschetti F, Dentith M D, List R D. A fractal-based algorithm for detect-ing first arrivals on seismic traces[J]. Geophysics, 1996, 61(4): 1095-1102.
[4] Anant K S, Dowla F U. Wavelet transform methods for phase identifica-tion in three-component seismograms[J]. Bulletin of the Seismological So-ciety of America, 1997, 87(6):1598-1612.
[5] Gibbons S J, Ringdal F, Kvarna T. Detection and characterization of seis-mic phases using continuous spectral estimation on incoherent and partial-ly coherent arrays[J]. Geophysical Journal International, 2008, 172(1): 405-421.
[6] Maeda N. A method for reading and checking phase times in auto-pro-cessing system of seismic wave data[J]. The Seismological Society Japan, 1985, 38(3): 365-379.
[7] Saragiotis C H D, Hadjileontiadis L J, Panas S M. PAI-S/K: A robust auto-matic seismic P phase arrival identification[J]. IEEE Transactions on Geo-science and Remote Sensing, 2002, 40(6): 1395-1404.
[8] 吴治涛, 李仕雄. STA/LTA算法拾取微地震事件P波到时对比研究[J]. 地球物理学进展, 2010, 25(5): 1577-1582. Wu Zhitao, Li Shixiong. Comparison of STA/LTA P-pickers for microseis-mic monitoring[J]. Progress in Geoρhysics, 2010, 25(5): 1577-1582.
[9] Earle P S, Shearer P M. Characterization of global seismograms using an automatic-picking algorithm[J]. Bulletin of the Seismological Society of America, 1994, 84(2): 366-376.
[10] 王彩霞, 白超英, 王馨. 地震震相初至自动检测技术综述[J]. 地球物理学进展, 2013, 28(5): 2363-2375. Wang Caixia, Bai Chaoying, Wang Xin. Review of automatic onset time picking for seismic arrival[J]. Progress in Geophysics, 2013, 28(5): 2363-2375.
[11] Sleeman R, van Eck T. Robust automatic P-phase picking: An on-line implementation in the analysis of broad-band seismogram recording[J]. Physics of the Earth and Planetary Interiors, 1999, 113(1): 265-275.
[12] Akaike H. Markovian representation of stochastic processes and its appli-cation to the analysis of autoregressive moving average processes[J]. An-nals of the Institute of Statistica Mathematics, 1974, 26(1): 363-387.
[13] 贾建科, 韩团军, 朱宁洲. 基于Matlab GUI的模拟带通滤波器的设计 [J]. 现代电子技术, 2010, 10: 58-62. Jia Jianke, Han Tuanjun, Zhun Ningzhou. Design of analog band-pass filter based on Matlab GUI[J]. Modern Electronic Technology, 2010, 10: 58-62.
[14] Glenn S R. S-Transform analysis of gravity wave activity from a small scale network of airglow imagers[D]. Ontario: The University of Western Ontario, 1999.
[15] 刘代志, 王仁明, 李夕海, 等. 基于小波包分解及AR模型的单通道地震波信号初至点检测[J]. 地球物理学报, 2005, 48(5): 1098-1102. Liu Daizhi, Wang Renming, Li Xihai, et al. Onset point identification of single channel seismic signal based on wavelet packet and the AR model [J]. Chinese Journal of Geohpysics, 2005, 48(5): 1098-1102.
