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.

A Generalized Reflection-Transmission Coefficient Matrix And Discrete Wavenumber Method For Synthetic Seismograms .2. For Multiple Sources At Different Depths

Migration of generalized random synthesis of seismic data

An Efficient Numerical Method for Computing Synthetic Seismograms for a Layered Half-space with Sources and Receivers at Close or Same Depths

A New Numerical Technique for Simulating the Coupled Seismic and Electromagnetic Waves in Layered Porous Media

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

Physical Interpretation of the Generalized Reflection and Transmission Coefficients

A Comparative Study On Seismic Wave Methods For Multi-Layered Media With Irregular Interfaces: Irregular Topography Problem

Generalized Synthesis of Seismic Data and Its Migration

Synthetic Seismograms by Normal-Mode Summation: a New Derivation and Numerical Examples

AN EFFICIENT METHOD TO MODEL SEISMIC PROPAGATION IN DIFFUSIVE-VISCOUS MEDIA WITH DIPPING INTERFACES

Numerical simulation of cross-hole seismic wave field

Determining the threshold value in upper limit of wavenumber integration by reflection-transmission coefficient in theoretical seismograms calculation

Complete Synthetic Seismograms Based on A Spherical Self-Gravitating Earth Model with an Atmosphere-Ocean-Mantle-Core Structure

New Developments In Analysis Of High-Frequency Rayleigh Waves

Generation of 6-C Synthetic Seismograms in Stratified Vertically Transversely Isotropic Media Using a Generalized Reflection and Transmission Coefficient Method

Wave Propagation in Multi‐layered Viscoelastic Media with Dipping Interfaces

A simple and fast method to calculate dispersion curves of Rayleigh waves due to models with a low-velocity half-space

A new method for estimation of Rayleigh wave group velocity spectrum

Generalized Transfer Matrix Method for Propagation of Surface Waves in Layered Azimuthally Anisotropic Half-Space

Surface-wave observations after integrating active and passive source data

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