Search

We develop a numerical algorithm for simulation of wave propagation in linear nonisothermal poroelastic media, based on Biot theory and a generalized Fourier law of heat transport in analogy with Maxwell model of viscoelasticity. A plane wave analysis indicates the presence of the classical P and S waves and two slow waves, namely, the Biot and the thermal slow modes of propagation, which present diffusive behavior under certain conditions, depending on viscosity, frequency, and the thermoelastic constants. The wavefield is computed with a direct meshing method using the Fourier differential operator to calculate the spatial derivatives. We propose two alternative time‐stepping algorithms, namely, a first‐order explicit Crank‐Nicolson method and a second‐order splitting method. The Fourier differential operator provides spectral accuracy in the calculation of the spatial derivatives. Modeling the thermal diffusive mode is relevant for high‐temperature high‐pressure fields and since it leads to mesoscopic attenuation by mode conversion of the fast waves to the thermal waves. Key Points: The differential equations of dynamic thermoporoelasticity are solved with a direct grid methodThe solver is based on the Fourier pseudospectral method to compute the spatial derivativesThe simulation reveals the behavior of four wave modes, namely, the fast P, slow P, S, and the thermal wave-diffusion field [ABSTRACT FROM AUTHOR]

We establish a generalization of the thermoelasticity wave equation to the porous case, including the Lord–Shulman (LS) and Green–Lindsay (GL) theories that involve a set of relaxation times ( τ i , i = 1 , ... , 4 ). The dynamical equations predict four propagation modes, namely, a fast P wave, a Biot slow wave, a thermal wave, and a shear wave. The plane-wave analysis shows that the GL theory predicts a higher attenuation of the fast P wave, and consequently a higher velocity dispersion than the LS theory if τ 1 = τ 2 > τ 3 , whereas both models predict the same anelasticity for τ 1 = τ 2 = τ 3 . We also propose a generalization of the LS theory by applying two different Maxwell–Vernotte–Cattaneo relaxation times related to the temperature increment ( τ 3 ) and solid/fluid strain components ( τ 4 ), respectively. The generalization predicts positive quality factors when τ 4 ≥ τ 3 , and increasing τ 4 further enhances the attenuation. The wavefields are computed with a direct meshing algorithm using the Fourier pseudospectral method to calculate the spatial derivatives and a first-order explicit Crank–Nicolson time-stepping method. The propagation illustrated with snapshots and waveforms at low and high frequencies is in agreement with the dispersion analysis. The study can be useful for a comprehensive understanding of wave propagation in high-temperature high-pressure fields. [ABSTRACT FROM AUTHOR]