A computational model is presented for simulations small amplitude wavefields in a deformable porous medium saturated with a compressible fluid under temperature variations. The model is based on the governing equations derived with the use of the theory of Symmetric Hyperbolic Thermodynamically Compatible (SHTC) systems in conjunction with the efficient staggered grid finite difference numerical method. It is numerically shown that the characteristics of wavefields in a saturated porous medium strongly depend on the porosity, which varies with temperature.