Se presenta un nuevo método para resolver problemas de conjugación térmica en modelos estratificados con un salto de conductividad térmica en sus límites. El propósito de este método es aproximar el coeficiente de conductividad térmica inversa, que presenta discontinuidades, utilizando una combinación de funciones escalón. Se construye un operador continuo generalizado en un espacio continuo de medios homogéneos por partes. Obtuvimos una solución analítica para el problema estacionario de conjugación de calor en el modelo estratificado con un grosor finito y con condiciones de Dirichlet-Neumann en los límites externos. Se construyó un algoritmo para la continuación descendente del flujo de calor a profundidades que corresponden a la parte superior de la capa del manto. Las ventajas de este método se ilustran mediante pruebas con datos sísmicos, de gravedad y geotérmicos de una zona de estudio en los Urales y regiones vecinas de Rusia. Examinamos las relaciones estadísticas entre la densidad y los parámetros térmicos y determinamos los componentes del flujo de calor para la corteza y el manto. El método permite una continuación descendente del flujo de calor hasta la base del manto superior y nos permite determinar los efectos térmicos de las características laterales y verticales de las estructuras tectónicas profundas.