Numerical Analysis of Seismoelectric Conversion in Stratified Low-Permeability Porous Rocks

A new set of complex dynamic localized response coefficients (LRC-NEW), related to the fast P waves, slow P waves, S waves, and electromagnetic (EM) waves, respectively, is derived to overcome the numerical inaccuracy occurring while computing the seismoelectric (SE) wave fields in stratified low-permeability porous media based on the reflectivity method. In addition, the dynamic electrokinetic coupling coefficient is modified to make it applicable for the electrical double layer of arbitrary thickness. After a thorough validation, we find that the EM waves simulated by using the known LRC expressions differ significantly with the analytical solutions when the static permeability is below 10 –15 m 2 , while those simulated by using LRC-NEW achieve excellent fits with the analytical solutions even when the static permeability is as low as 10 –20 m 2 . A sensitivity study demonstrates that an optimum static permeability exists where the amplitude of evanescent SE responses reaches its maximum. The normalized root-mean-square amplitude of the evanescent SE signals in low-permeability rocks such as shales is more sensitive to permeability than seismic waves, even when the interdependence of porosity and static permeability is considered. Sensitivity studies also manifest that evanescent SE responses are sensitive to salinity and water saturation even in low-permeability shales. Notably, the amplitude of the evanescent SE conversions increases monotonically with water saturation, which is much simpler than the relationship between the seismic amplitude and water saturation. Therefore, the SE method has the potential to characterize the permeability, salinity, and water saturation of low-permeability reservoirs.