Wavefield simulation and analysis with finite difference method in partially saturated double-porosity media
Shi Zhi-Qi,He Xiao,Liu Lin,Chen De-Hua,,,
DOI: https://doi.org/10.7498/aps.73.20240227
IF: 0.906
2024-01-01
Acta Physica Sinica
Abstract:Double-porosity poroelastic models, which account for the effect of mesoscopic flow in heterogeneous rocks on wave dispersion and attenuation, are useful for quantitative seismic interpretation. Wavefield simulation based on double-porosity models not only helps visualize the propagation characteristics of the elastic waves but also lays the foundation for seismic imaging. In this paper, we perform wavefield simulation and analysis based on the Santos-Rayleigh model which incorporates mesoscopic and global flow in partially-saturated double-porosity media. Specifically, the mesoscopic flow mechanism is represented with a Zener viscoelastic model. The comparison shows that the Zener model can accurately capture the propagation characteristics of fast P-wave, but fails to represent the attenuation characteristics of slow P3 wave in the low-frequency band. It suggests that Zener viscoelastic model and slow wave modes follow different mechanisms. Then staggered grid finite-difference method is used to simulate wave propagation in double-porosity media, and the stiff problem is solved with a time-splitting algorithm, which can significantly improve computational efficiency. Based on above methods, the correctness of our algorithm is verified with derived analytical solution for a P-wave source in a uniform partially saturated poroelastic media. Analytical and numerical solutions are in good agreement and mean error is 0.33%. We provide some examples of wavefield snapshots and seismograms in homogeneous and layered heterogeneous media at seismic and ultrasonic frequencies. Simulation results demonstrate the strong attenuation of fast P-wave and no change of S-wave in the seismic band due to mesoscopic flow mechanism, which is consistent with the theoretical predictions of double-porosity model. Moreover, energy of fast P-wave is concentrated in solid phase while slow waves are stronger in fluid phases. This work contributes to the understanding of broadband elastic wave propagation in heterogeneous partially saturated porous media and can be applied in the reservoir imaging with broadband geophysical data.
physics, multidisciplinary