Table of Contents Author Guidelines Submit a Manuscript
Shock and Vibration
Volume 2016, Article ID 4051701, 6 pages
Research Article

Algorithm Indicating Moment of P-Wave Arrival Based on Second-Moment Characteristic

1KGHM CUPRUM Ltd., CBR, Sikorskiego 2-8, 53-659 Wrocław, Poland
2Faculty of Pure and Applied Mathematics, Hugo Steinhaus Center, Wroclaw University of Science and Technology, Janiszewskiego 14a, 50-370 Wrocław, Poland
3KGHM Polska Miedź S.A. O/ZG “Rudna”, Henryka Dąbrowskiego 50, 59-100 Polkowice, Poland

Received 22 April 2016; Accepted 1 August 2016

Academic Editor: Shimin Liu

Copyright © 2016 Jakub Sokolowski et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


The moment of P-wave arrival can provide us with many information about the nature of a seismic event. Without adequate knowledge regarding the onset moment, many properties of the events related to location, polarization of P-wave, and so forth are impossible to receive. In order to save time required to indicate P-wave arrival moment manually, one can benefit from automatic picking algorithms. In this paper two algorithms based on a method finding a regime switch point are applied to seismic event data in order to find P-wave arrival time. The algorithms are based on signals transformed via a basic transform rather than on raw recordings. They involve partitioning the transformed signal into two separate series and fitting logarithm function to the first subset (which corresponds to pure noise and therefore it is considered stationary), exponent or power function to the second subset (which corresponds to nonstationary seismic event), and finding the point at which these functions best fit the statistic in terms of sum of squared errors. Effectiveness of the algorithms is tested on seismic data acquired from O/ZG “Rudna” underground copper ore mine with moments of P-wave arrival initially picked by broadly known STA/LTA algorithm and then corrected by seismic station specialists. The results of proposed algorithms are compared to those obtained using STA/LTA.