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.

No abstract

No abstract

Summary We analyze the reflection coefficient of an inhomogeneous plane wave incident on the thermally insulated surface of a thermo-poroelastic medium. The theory, which includes the classic Lord-Shulman (LS) and Green-Lindsay (GL) theories as well as a generalization of the LS model, predicts three inhomogeneous longitudinal waves and one transverse wave, described by potential functions specified by the propagation direction and inhomogeneity angle. The GL model can give a stronger P1-wave thermal attenuation and consequently a stronger velocity dispersion than the LS model. We investigate the influence of inhomogeneity angle, type of incident wave, frequency and surface boundary conditions. The generalized LS model exhibits increased P1-wave thermal attenuation with increasing Maxwell-Vernotte-Cattaneo (MVC) relaxation time and consequently predicts more interference energy, irrespective if the surface is open or sealed. The inhomogeneity angle affects the energy partitions particularly near the grazing incidence, with a significant interference energy, which must be taken into account to satisfy the energy conservation. The thermal dispersion occurs at frequencies around the thermal relaxation peak, which moves to low frequencies when the conductivity increases.