Full Information Numerical Simulation of Two-Dimensional Steady Ground Temperature Field

An efficient and precise numerical simulation method of geothermal field is of great significance for the exploitation, evaluation, and sustainable development of geothermal resources. In the face of complex geological conditions, the conventional numerical simulation algorithm of geothermal field has many problems, such as large memory consumption and long calculation times. To solve these problems, this paper proposes an efficient and high-precision numerical simulation method for a full information geothermal field. The temperature field is divided into background field and anomalous field for calculation from the perspective of the superposition of temperature field. The background field is a uniform layered medium, and the anomalous field is converted into one-dimensional ordinary differential equations under different wavenumbers via one-dimensional full information Fourier transform. Memory consumption and computational time are significantly reduced. The equations under different wavenumbers are independent of each other and highly parallel, while retaining vertical spatial domain, with strict upper and lower boundary conditions, which can adapt to the needs of complex terrain simulation. The finite element method is used to solve the one-dimensional ordinary differential equation, and the Thomas algorithm is used to efficiently solve the linear equations with fixed bandwidth, which improves the computational efficiency and realizes the efficient and high-precision numerical simulation of the ground temperature field. The results of the algorithm were compared with COMSOL Multiphysics finite element simulation software to verify the correctness of the algorithm and compare and analyze the improvement of the algorithm in computation time and memory consumption under the same mesh division.