420 likes | 438 Vues
On the Physics and Simulation of Waves at Fluid-Solid Interfaces: Application to NDT, Seismic Exploration and Earthquake Seismology by José M. Carcione (OGS, Italy). The 2D modeling algorithm. 2-D Equations of Motion. Euler-Newton’s Equations:. Constitutive Equations:. Memory Variables:.
E N D
On the Physics and Simulation of Waves at Fluid-Solid Interfaces:Application to NDT, Seismic Exploration and Earthquake SeismologybyJosé M. Carcione (OGS, Italy)
2-D Equations of Motion Euler-Newton’s Equations: Constitutive Equations: Memory Variables:
Scholte wave dispersion equation Relevant roots: Scholte wave Leaky Rayleigh wave
Inhomogeneous waves Elliptical polarization Plane wave
Numerical algorithm Two grids (domain decomposition): ocean and oceanic crust Spatial derivatives Fourier method in the horizontal direction Chebyshev method in the vertical direction Time integration 4th-order Runge-Kutta
AVA analysis Elastic case Anelastic case
Water/plexiglass (soft bottom) No leaky Rayleigh wave
Test with analytical solution Water/plexiglass interface
Test with analytical solution Water/glass interface
Dispersive Scholte waves North Sea. 70 m water depth. Airgun source. Elastic case Anelastic case
Ocean overlying the crust Phase velocity
Ocean overlying the crust Group velocity Dissipation factor
Ocean overlying the crust Ben_Menahem and Singh (1981) Experimental data (Fig. 10.3) Attenuation coefficient
Ocean overlying the crust Phase/group velocities
Ocean overlying the crust High-frequency case Elastic and anelastic solutions
Ocean overlying the crust Low-frequency case Elastic Anelastic
Sediment layer overlying the crust Low-frequency case Elastic Anelastic
Geological model From CRUST 5.1
The Kelvin-Voigt stress-strain relation s =stress components e= strain components u=displacements l, m=Lamé constants l’, m’=damping Lamé constants
Input damping parameters w0=reference frequency QP0 = reference P-wave quality factor QS0 = reference S-wave quality factor
The equations of motion v = particle velocity r= density f=body forces
Tests with analytical solutions Rayleigh waves -- Cagniard-de Hoop solution Pekeris (1955) solution -- unbounded media
Simulation of Rayleigh waves. Seismograms. Lossless case
Simulation of Rayleigh waves. Seismograms. Lossy case
Simulation of Love waves. Seismograms. Lossless case Lossy case
Conclusions Effects of anelastic attenuation Inhomogeneous viscoelastic waves Differences at critical and post-critical angles Rayleigh-window effect Pseudospectral numerical method Verified for reflection/transmission and interface waves Effective tool for seismic exploration studies, NDT and earthquake seismology