Automatic detection and characterisation of the first P- and S-wave pulse in rocks using ultrasonic transmission method

Abstract In this paper, we propose a new methodology for the automatic picking of P- and S- wave arrivals. The measurement of elastic wave velocities is crucial for characterisation of the elastic properties of rocks, as well as their physico-mechanical properties and durability. P-waves are easy to generate and acquire, whereas the acquisition and subsequent analysis of the S-waves present major difficulties in most investigations of rocks. In our proposed method the recorded signal is first pre-processed in the wavelet domain, removing low-frequency disturbances and filtering noise. The signal is then analysed in the time-domain to detect and characterise the first pulse and, subsequently, estimate the corresponding onset time. The automatic approach analyses all pulses detected in the output signal and selects the first pulse of the P or S wave relative to symmetry, amplitude and duration criteria. This triple check provides greater confidence in the obtained results. We record P- and S- waveforms through the transmission method for a broad range of sedimentary, igneous and metamorphic rocks used in buildings. Results are compared to manual picking, which is considered as a true or reference value. The recorded signals show that microstructural components of rocks have a strong influence on the output signal. Mineralogical composition, porosity and particle size affect the wave velocity, attenuation, wavelength and waveform, which in turn influences the manual picking of the onset time. The proposed methodology provides precise values of the onset time of P- and S- waves. The discrepancies between automatic and manual P- and S- wave picking are 0.7 and 4.4%, respectively. P-wave signals have a high signal-to-noise ratio and their arrival times are clear and easy to determine. However, the arrival time of S-waves presents significant problems, mainly for medium to coarse-grained rocks. In this case, the output signal is contaminated by P-waves and has a lower signal-to-noise ratio. Besides, propagation of the S-wave through the rock affects the frequency and waveform of the first pulse, which complicates manual picking of the onset time. The developed methodology distinguishes between the S-wave and the remaining P-waves contained in the transversal signal, independently of the skill and subjectivity of a human analyst.

[1]  M. Leonard,et al.  Comparison of Manual and Automatic Onset Time Picking , 2000 .

[2]  S. Demirdag The effect of using different polymer and cement based materials in pore filling applications on technical parameters of travertine stone , 2009 .

[3]  N. Magotra,et al.  Single-station seismic event detection and location , 1989 .

[4]  Jihun Choi,et al.  Improved modified energy ratio method using a multi-window approach for accurate arrival picking , 2017 .

[5]  Arye Nehorai,et al.  Analysis of a polarized seismic wave model , 1996, IEEE Trans. Signal Process..

[6]  Yifan Zhao,et al.  Comparison of alternatives to amplitude thresholding for onset detection of acoustic emission signals , 2017 .

[7]  Reinoud Sleeman,et al.  Robust automatic P-phase picking: an on-line implementation in the analysis of broadband seismogram recordings , 1999 .

[8]  Federico Botella,et al.  A real-time earthquake detector with prefiltering by wavelets , 2003 .

[9]  Jerry M. Mendel,et al.  First break refraction event picking using fuzzy logic systems , 1994, IEEE Trans. Fuzzy Syst..

[10]  D. Benavente,et al.  Impact of salt and frost weathering on the physical and durability properties of travertines and carbonate tufas used as building material , 2018, Environmental Earth Sciences.

[11]  Juan J. Galiana-Merino,et al.  Seismic $P$ Phase Picking Using a Kurtosis-Based Criterion in the Stationary Wavelet Domain , 2008, IEEE Transactions on Geoscience and Remote Sensing.

[12]  Sergio I. Molina,et al.  Wavelet Transform Methods for Azimuth Estimation in Local Three-Component Seismograms , 2007 .

[13]  David Benavente,et al.  Spatial attenuation: The most sensitive ultrasonic parameter for detecting petrographic features and decay processes in carbonate rocks , 2011 .

[14]  D. Marcotte,et al.  Correlations between compressional and shear wave velocities and corresponding Poisson's ratios for some common rocks and sulfide ores , 2009 .

[15]  S. Mostafa Mousavi,et al.  Automatic microseismic denoising and onset detection using the synchrosqueezed continuous wavelet transform , 2016 .

[16]  Christopher John Young,et al.  A comparison of select trigger algorithms for automated global seismic phase and event detection , 1998, Bulletin of the Seismological Society of America.

[17]  Hengchang Dai,et al.  Automatic picking of seismic arrivals in local earthquake data using an artificial neural network , 1995 .

[18]  J. J. Giner,et al.  De-Noising of Short-Period Seismograms by Wavelet Packet Transform , 2003 .

[19]  J. G. Cai,et al.  Experimental study of ultrasonic wave attenuation across parallel fractures , 2006 .

[20]  Said Gaci,et al.  The Use of Wavelet-Based Denoising Techniques to Enhance the First-Arrival Picking on Seismic Traces , 2014, IEEE Transactions on Geoscience and Remote Sensing.

[21]  D. Benavente,et al.  Sedimentary structures and physical properties of travertine and carbonate tufa building stone , 2012 .

[22]  W. Yi,et al.  Propagation characteristics of acoustic emission wave in reinforced concrete , 2017 .

[23]  Joël Sarout,et al.  A semi-automatic processing technique for elastic-wave laboratory data. , 2009, Ultrasonics.

[24]  H. Akaike Markovian representation of stochastic processes and its application to the analysis of autoregressive moving average processes , 1974 .

[25]  Miguel Gómez-Heras,et al.  Non-linear decay of building stones during freeze–thaw weathering processes , 2013 .

[26]  Peter M. Shearer,et al.  Characterization of global seismograms using an automatic-picking algorithm , 1994, Bulletin of the Seismological Society of America.

[27]  Alberto Michelini,et al.  Automatic picking of P and S phases using a neural tree , 2006 .

[28]  D. Benavente,et al.  The influence of petrophysical properties on the salt weathering of porous building rocks , 2007 .

[29]  Juan J. Galiana-Merino,et al.  SeismicWaveTool: Continuous and discrete wavelet analysis and filtering for multichannel seismic data , 2013, Comput. Phys. Commun..

[30]  Z. Lafhaj,et al.  Microcracking based rock classification using ultrasonic and porosity parameters and multivariate analysis methods , 2013 .

[31]  J. Vidale Waveform effects of a high‐velocity, subducted slab , 1987 .

[32]  Clifford H. Thurber,et al.  Automatic P-Wave Arrival Detection and Picking with Multiscale Wavelet Analysis for Single-Component Recordings , 2003 .

[33]  M. Gómez-Heras,et al.  Improved correlation between the static and dynamic elastic modulus of different types of rocks , 2016 .

[34]  Jochen H Kurz,et al.  Strategies for reliable automatic onset time picking of acoustic emissions and of ultrasound signals in concrete. , 2005, Ultrasonics.

[35]  Giuseppe Lacidogna,et al.  Reliable onset time determination and source location of acoustic emissions in concrete structures , 2012 .

[36]  R. V. Allen,et al.  Automatic earthquake recognition and timing from single traces , 1978, Bulletin of the Seismological Society of America.

[37]  G. Crosta,et al.  Ultrasonic and X-ray computed tomography characterization of progressive fracture damage in low-porous carbonate rocks , 2016 .

[38]  H. Akaike Markovian Representation of Stochastic Processes and Its Application to the Analysis of Autoregressive Moving Average Processes , 1974 .