The application of the unequal interval difference method in the terrestrial heat field distribution simulation

Aiming at simulating the unsteady, heterogeneous terrestrial heat field distribution with local anomalies, this article adopted unequal interval difference method to deduce the numerical solution of two-dimensional heat conduction equation. Based on the deduction, we dynamically simulate a terrestrial heat field distribution with local anomalies, where heat conductivity, density and specific heat are 1.5(W/M·°C), 2500(kg/m³) and 1500(J/K·°C). The simulation shows a conclusion, that it spends 0.4 million years that the distribution recover to initial equilibrium state, which has important significance in hydrocarbon pool forming dynamics simulation and evaluation.