The stability of poro-elastic wave equations in saturated porous media

Poro-elastic wave equations are one of the fundamental problems in seismic wave exploration and applied mathematics. In the past few decades, elastic wave theory and numerical method of porous media have developed rapidly. However, the mathematical stability of such wave equations have not been fully studied, which may lead to numerical divergence in the wave propagation simulation in complex porous media. In this paper, we focus on the stability of the wave equation derived from Tuncay’s model and volume averaging method. By analyzing the stability of the first-order hyperbolic relaxation system, the mathematical stability of the wave equation is proved for the first time. Compared with existing poro-elastic wave equations (such as Biot’s theory), the advantage of newly derived equations is that it is not necessary to assume uniform distribution of pores. Such wave equations can spontaneously incorporate complex microscale pore/fracture structures into large-scale media, which is critical for unconventional oil and gas exploration. The process of proof and numerical examples shows that the wave equations are mathematically stable. These results can be applied to numerical simulation of wave field in reservoirs with pore/fracture networks, which is of great significance for unconventional oil and gas exploration.