Geometric effect on a laboratory-scale wavefield inferred from a three-dimensional numerical simulation

Abstract The coda part of a waveform transmitted through a laboratory sample should be examined for the high-resolution monitoring of the sample characteristics in detail. However, the origin and propagation process of the later phases in a finite-sized small sample are very complicated with the overlap of multiple unknown reflections and conversions. In this study, we investigated the three-dimensional (3D) geometric effect of a finite-sized cylindrical sample to understand the development of these later phases. This study used 3D finite difference method simulation employing a free-surface boundary condition over a curved model surface and a realistic circular shape of the source model. The simulated waveforms and the visualized 3D wavefield in a stainless steel sample clearly demonstrated the process of multiple reflections and the conversions of the P and S waves at the side surface as well as at the top and bottom of the sample. Rayleigh wave propagation along the curved side boundary was also confirmed, and these waves dominate in the later portion of the simulated waveform with much larger amplitudes than the P and S wave reflections. The feature of the simulated waveforms showed good agreement with laboratory observed waveforms. For the simulation, an introduction of an absorbing boundary condition at the top and bottom of the sample made it possible to efficiently separate the contribution of the vertical and horizontal boundary effects in the simulated wavefield. This procedure helped to confirm the additional finding of vertically propagating multiple surface waves and their conversion at the corner of the sample. This new laboratory-scale 3D simulation enabled the appearance of a variety of geometric effects that constitute the later phases of the transmitted waves.

[1]  Ulrich Wegler,et al.  Radiative transfer of elastic waves versus finite difference simulations in two‐dimensional random media , 2006 .

[2]  Richard D. Miller,et al.  An improved vacuum formulation for 2D finite-difference modeling of Rayleigh waves including surface topography and internal discontinuities , 2012 .

[3]  A. Levander Fourth-order finite-difference P-SV seismograms , 1988 .

[4]  Hiroshi Takenaka,et al.  FDM simulation of seismic-wave propagation for an aftershock of the 2009 Suruga bay earthquake: Effects of ocean-bottom topography and seawater layer , 2012 .

[5]  Xinglin Lei,et al.  Laboratory studies of seismic wave propagation in inhomogeneous media using a laser doppler vibrometer , 1997, Bulletin of the Seismological Society of America.

[6]  J. Hazzard,et al.  Numerical investigation of induced cracking and seismic velocity changes in brittle rock , 2004 .

[7]  T. Blum,et al.  Full-wavefield modeling and reverse time migration of laser ultrasound data: A feasibility study , 2015 .

[8]  Y. Fukushima,et al.  Relationship between Fluctuations of Arrival Time and Energy of Seismic Waves and Scale Length of Heterogeneity: An Inference from Experimental Study , 2001 .

[9]  Robert W. Graves,et al.  Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences , 1996, Bulletin of the Seismological Society of America.

[10]  J. B. Walsh,et al.  Changes in seismic velocity and attenuation during deformation of granite , 1977 .

[11]  Takashi Furumura,et al.  Seismic‐ and Tsunami‐Wave Propagation of the 2011 Off the Pacific Coast of Tohoku Earthquake as Inferred from the Tsunami‐Coupled Finite‐Difference Simulation , 2013 .

[12]  R. Snieder,et al.  Scattering amplitude of a single fracture under uniaxial stress , 2014 .

[13]  Surface waves on cylindrical solids: numerical and experimental study. , 2013, Ultrasonics.

[14]  M. Möllhoff,et al.  Rock fracture compliance derived from time delays of elastic waves , 2010 .

[15]  M. Nafi Toksöz,et al.  Discontinuous-Grid Finite-Difference Seismic Modeling Including Surface Topography , 2001 .

[16]  S. Vialle,et al.  Laboratory measurements of elastic properties of carbonate rocks during injection of reactive CO2‐saturated water , 2011 .

[17]  Bernard A. Chouet,et al.  A free-surface boundary condition for including 3D topography in the finite-difference method , 1997, Bulletin of the Seismological Society of America.

[18]  J. Scales,et al.  Laser characterization of ultrasonic wave propagation in random media. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[19]  Takashi Furumura,et al.  Scattering of high-frequency seismic waves caused by irregular surface topography and small-scale velocity inhomogeneity , 2015 .

[20]  Osamu Nishizawa,et al.  A performance study of a laser Doppler vibrometer for measuring waveforms from piezoelectric transducers , 2009, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[21]  Moshe Reshef,et al.  A nonreflecting boundary condition for discrete acoustic and elastic wave equations , 1985 .

[22]  T. Furumura,et al.  FDM Simulation of Seismic Waves, Ocean Acoustic Waves, and Tsunamis Based on Tsunami-Coupled Equations of Motion , 2011, Pure and Applied Geophysics.

[23]  H. Yukutake,et al.  Fracturing process of granite inferred from measurements of spatial and temporal variations in velocity during triaxial deformations , 1989 .

[24]  M. Toksöz,et al.  Radiation Patterns Of Compressional And Shear Transducers At The Surface Of An Elastic Half Space , 1994 .

[25]  M. Nakatani,et al.  Monitoring frictional strength with acoustic wave transmission , 2008 .

[26]  Wei Wei,et al.  Finite difference modeling of ultrasonic propagation (coda waves) in digital porous cores with un-split convolutional PML and rotated staggered grid , 2014 .

[27]  竹中 博士,et al.  速度・応力型差分法での固体・流体境界の扱いについて , 2005 .

[28]  S. Shapiro,et al.  Modeling the propagation of elastic waves using a modified finite-difference grid , 2000 .