Analytical regularization of the high-frequency instability problem in numerical simulation of seismoelectric wave-fields in multi-layered porous media

Abstract

A new method of numerical simulation of seismoelectric wave-fields in multi-layered porous media has been developed by extending Chen′s technique of computing synthetic seismograms to coupled seismic and electromagnetic waves. However, the existence of an exponential growth factor leads to numerical instability for high-frequencies. A natural regularization approach eliminating the high-frequency instability is to create a fictitious thin flat layer including the source point in the numerical calculation, which requires an adequate thickness of the fictitious source layer as well as some extra computation cost. In this article, an analytical regularization approach is developed to deal with the high-frequency instability problem in numerical simulations of seismoelectric wave-fields in multi-layered porous media. The results show that the analytical regularization approach is much more effective in solving the above high-frequency instability problem than the natural regularization approach.