A Method for Quantitative Interpretation of Stationary Thermal Fields for Layered Media

Abstract

A new method to solve thermal conjugacy problems is presented for layered models with a thermal conductivity jump at their boundaries. The purpose of this method is to approximate the inverse thermal conductivity coefficient, which has breaks, by using a combination of step functions. A generalized continuous operator is constructed in a continuous space of piecewise–homogeneous media. We obtained an analytical solution for the stationary problem of heat conjugacy in the layered model with finite thickness and with Dirichlet–Neumann conditions at the external boundaries. An algorithm was constructed for downward continuation of the heat flux to depths that correspond to the top of the mantle layer. The advantages of this method are illustrated by testing the crustal seismic, gravity and geothermal data of a study area in the Urals and neighboring regions of Russia. We examined statistical relations between density and thermal parameters and determined heat flux components for the crust and the mantle. The method enables a downward continuation of the heat flux to the base of the upper mantle and allows us to determine the thermal effects of the lateral and vertical features of deep tectonic structures.