Abstract:For a shallow or ground surface source and receiver at the same level or close depth, it is very difficult or computationally inefficient to simulate seismoelectric or electroseismic wave-fields in stratified porous media by current reflectivity methods, such as the Luco-Apsel-Chen generalized reflection and transmission method (LAC GRTM). In this work, the peak-trough averaging method which has been proved effective and efficient in dealing with this kind of computational problem is adopted to update the seismoelectric and electroseismic modeling algorithm based on LAC GRTM. After thoroughly verifying the accuracy and computational efficiency of the updated algorithm, we utilize it to numerically investigate both the electroseismic and seismoelectric couplings. Snapshots of electroseismic wave-fields indicate evanescent electroseismic conversion, a reverse process of evanescent seismoelectric conversion, dominates at relatively larger ratios of seismic wavelength to interface depth, whereas the interfacial radiation electroseismic conversion is more prominent for the opposite situation. Our seismoelectric modeling results demonstrate that electric signals can arrive at the ground surface a few milliseconds earlier than their causative seismic signals due to evanescent seismoelectric conversion. This is the first modeling result considering source-receiver geometries on the surface capable of explaining similar phenomena reported in geophysical field observations of seismically induced electrokinetic effects over a long history. The updated algorithm offers an accurate and efficient tool for forward modeling and will benefit interpretations of field observations as well as future inversion studies. Plain Language Summary For the case of active seismoelectric or electroseismic measurements with source and receivers located on or near the ground surface, it is computationally difficult for current reflectivity methods, such as Luco-Apsel-Chen generalized reflection and transmission method (LAC GRTM), to simulate the wave-fields in stratified porous media. In this study, the peak-trough averaging method is adopted to update the LAC-GRTM-based seismoelectric and electroseismic modeling algorithms to solve this computational problem. We verify the accuracy and efficiency of the updated algorithm and apply it to investigate the electroseismic and seismoelectric couplings. Two kinds of electroseismic conversions were identified by our modellings. When the ratio of seismic wavelength to interface depth is large, evanescent electroseismic conversion will be dominant; otherwise, the interfacial radiation electroseismic conversion will become prominent. In our seismoelectric modellings, it is found electric signals recorded by ground surface receivers could be several milliseconds earlier than their causative seismic signals due to evanescent seismoelectric conversion. This is the first modeling result capable of explaining similar phenomena that have been reported in seismoelectric field measurements over a long history. Generally speaking, the updated algorithm is a useful forward modeling tool, which could benefit interpretations of field data and future inversion studies.

The Doppler Curves of Different Seismic Phases in Acoustic-seismic Coupling Signals

The Seismic Responses and Its Doppler Effects of Moving Aircraft Caused by Acoustic-Seismic Coupling

Theoretical Analysis and Validation of Multiple Modes Doppler Curves

Sound Field Separating on Arbitrary Surfaces Enclosing a Sound Scatterer Based on Combined Integral Equations.

Numerical Investigation of Rayleigh-wave Propagation on Topography Surface

Seismic stress path induced by obliquely incident waves and its simulation

Investigation and Use of Surface-wave Characteristics for Near-surface Applications

Simple equations guide high-frequency surface-wave investigation techniques

Progress and prospect of surface-wave imaging techniques in near-surface applications

Discussion on Some Practical Equations with Implications to High-Frequency Surface-Wave Techniques

Power-law frequency-dependent Q simulations in viscoacoustic media using decoupled fractional Laplacians

Phase-velocity inversion from data-based diffraction kernels: seismic Michelson interferometer

A Three-dimensional Immersed Boundary Method for Accurate Simulation of Acoustic Wavefields with Complex Surface Topography

Seismoelectric and Electroseismic Modeling in Stratified Porous Media with a Shallow or Ground Surface Source

A Comprehensive Comparison Between the Refraction Microtremor and Seismic Interferometry Methods for Phase-Velocity Estimation

Modeling obliquely incident seismic wave propagation in rocks via the FDM-DEM coupling scheme: Algorithms and applications

New Developments In Analysis Of High-Frequency Rayleigh Waves

Power-Law Frequency-Dependent Q Simulation and Reverse-Time Migration in Transversely Isotropic Acoustic Media

The exponentiated phase measurement, and objective-function hybridization for adjoint waveform tomography

Numerical Investigation of Implementation of Air-Earth Boundary by Acoustic-Elastic Boundary Approach

Near-field seismoacoustic wave scattering due to an irregular interface: a unified framework