Calculation of temperature distribution in two-dimensional geothermal profile

Summary The temperature-depth distribution along the 250 km long Czechoslovak part of DSS profile VII was calculated to a depth of 40 km. The two-dimensional solution of the steadystate equation of heat conduction in an inhomogeneous medium, obtained by means of the finite difference method, was used. By resolving the observed surface heat flow into two components, one of which is generated by heat sources inside the model and the other is connected with the heat flow from deeper parts of the Earth, the heat flow at the lower boundary of the model was estimated and used here as a boundary condition in the numerical solution. To reduce the differences between the calculated and observed surface heat flow, the heat sources in the upper 10 km of model were optimized. A regularization process, in which the variation of the sought function was limited, was used both to estimate the heat flow and to optimize the sources.

---

u n¶rt;u n ¶rt; (¶rt;u 250 ) nu . . VII ¶rt;u 40 . un ¶rt; u u nn¶rt;u ¶rt; u ¶rt;¶rt; ¶rt; ¶rt; . n u u n n n ¶rt; u, ¶rt; u u uuu n u ¶rt;u, a u n, n uu ¶rt;u u un ¶rt; u u nu u uu u. u na¶rt;u ¶rt; u u u nu nu ¶rt; nuuu uu n u 10 ¶rt;u. nu , u nu nuuuu uu n un n uuu, uuu uau u uu .