Seismoelectric and Electroseismic Modeling in Stratified Porous Media with a Shallow or Ground Surface Source
Xu-Zhen Zheng,Hengxin Ren,Karl E. Butler,Haiming Zhang,Yao-Chong Sun,Wei Zhang,Qinghua Huang,Xiaofei Chen
DOI: https://doi.org/10.1029/2021jb021950
2021-01-01
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.