Simulation of thermoelastic waves based on the Lord-Shulman theory

Thermoelasticity is important in seismic propagation due to the effects related to wave attenuation and velocity dispersion. We have applied a novel finite-difference (FD) solver of the Lord Shulman thermoelasticity equations to compute synthetic seismograms that include the effects of the thermal properties (expansion coefficient, thermal conductivity, and specific heat) compared with the classic forward-modeling codes. We use a time splitting method because the presence of a slow quasistatic mode (the thermal mode) makes the differential equations stiff and unstable for explicit time-stepping methods. The spatial derivatives are computed with a rotated staggered-grid FD method, and an unsplit convolutional perfectly matched layer is used to absorb the waves at the boundaries, with an optimal performance at the grazing incidence. The stability condition of the modeling algorithm is examined. The numerical experiments illustrate the effects of the thermoelasticity properties on the attenuation of the fast P-wave (or E-wave) and the slow thermal P-wave (or T-wave). These propagation modes have characteristics similar to the fast and slow P-waves of poroelasticity, respectively. The thermal expansion coefficient has a significant effect on the velocity dispersion and attenuation of the elastic waves, and the thermal conductivity affects the relaxation time of the thermal diffusion process, with the T mode becoming wave-like at high thermal conductivities and high frequencies.