Forward modeling method, device, equipment and medium for three-dimensional steady-state heat conduction geothermal field

Abstract

The forward modeling method, device, equipment and medium of the three-dimensional steady-state heat conduction geothermal field, construct a three-dimensional prism model of the target area containing a three-dimensional anomaly body and mesh it, assign the thermophysical parameters of each node in the three-dimensional prism model according to the distribution of the thermophysical parameters of the target area, load the boundary conditions, calculate the space domain background field temperature, and use it as the initial space domain temperature total field; calculate the space wave number mixed domain abnormal temperature field based on the thermophysical parameters and the current space domain temperature total field; perform Fourier inverse transformation on the space wave number mixed domain abnormal temperature field to obtain the space domain abnormal field temperature; obtain a new space domain temperature total field based on the space domain background field temperature and the space domain abnormal field temperature; and continuously iterate until the new space domain temperature total field obtained by iteration meets the iteration convergence condition. The present invention realizes the efficient and high-precision numerical simulation of the three-dimensional steady-state heat conduction geothermal field.

Description

Technical Field

[0001] The present invention belongs to the technical field of geothermal field numerical simulation, and in particular relates to a forward modeling method, device, equipment and medium for a three-dimensional steady-state heat conduction geothermal field. Background Art

[0002] With the development of social and economic environment, the demand for green energy is increasing day by day. Geothermal energy as a clean renewable energy is gaining more and more attention. However, due to the complexity of underground geological conditions, the exploration, mining and utilization of geothermal resources are facing challenges, so the numerical simulation technology of geothermal field is in urgent need of development.

[0003] The current defects of numerical simulation of geothermal fields are: conventional numerical simulation methods of geothermal fields (such as finite element method, finite difference method, boundary element method, finite volume method) are all calculated in the spatial domain. In order to achieve higher accuracy requirements in large-scale complex media, the model needs to be finely divided. The number of divided units is large, resulting in large storage volume and long calculation time. Summary of the invention

[0004] In view of the current situation that conventional methods such as the finite element method, the finite difference method, the boundary element method, etc. have large computational complexity and high storage requirements when processing large-scale numerical simulation of temperature fields, the present invention aims to propose a forward modeling method, device, equipment and medium for three-dimensional steady-state heat conduction geothermal field to meet the needs of refined forward modeling calculations for undulating terrain.

[0005] In order to achieve the above technical objectives, the technical solution proposed by the present invention is:

[0006] In one aspect, the present invention provides a forward modeling method for a three-dimensional steady-state heat conduction geothermal field, comprising:

[0007] S1. Determine a three-dimensional anomaly, determine a target area containing the three-dimensional anomaly, and construct a three-dimensional prism model of the target area;

[0008] S2. Meshing the three-dimensional prism model along the x, y, and z directions to obtain a series of mesh nodes and meshing parameters of the three-dimensional prism model;

[0009] S3. assigning values ​​of the thermophysical parameters of each subdivision node in the three-dimensional prism model according to the distribution of the thermophysical parameters of the target area;

[0010] S4. Load boundary conditions, calculate the space domain background field temperature, and use it as the initial space domain temperature total field;

[0011] S5. Calculate the abnormal temperature field in the spatial wavenumber mixed domain based on the thermophysical parameters and the current total temperature field in the spatial domain;

[0012] S6. Perform inverse Fourier transform on the abnormal temperature field in the spatial wavenumber mixed domain to obtain the abnormal field temperature in the spatial domain;

[0013] S7. based on the space domain background field temperature and the space domain abnormal field temperature, obtain a new space domain total temperature field;

[0014] S8. Set the iteration convergence condition. If the new spatial domain temperature total field meets the iteration convergence condition, output the new spatial domain temperature total field. Otherwise, use the new spatial domain temperature total field as the current spatial domain temperature total field in the next iteration and return to S5.

[0015] In another aspect, the present invention provides a forward modeling device for a three-dimensional steady-state heat conduction geothermal field, comprising:

[0016] The first module is used to determine a three-dimensional anomaly, determine a target area containing the three-dimensional anomaly, and construct a three-dimensional prism model of the target area;

[0017] The second module is used to mesh the three-dimensional prism model along the x, y, and z directions to obtain a series of mesh nodes and mesh parameters of the three-dimensional prism model;

[0018] The third module is used to assign values ​​of the thermophysical parameters of each segmentation node in the three-dimensional prism model according to the distribution of the thermophysical parameters of the target area;

[0019] The fourth module is used to load boundary conditions and calculate the background field temperature of the space domain, and use it as the total temperature field of the initial space domain;

[0020] The fifth module is used to calculate the abnormal temperature field in the spatial wavenumber mixed domain based on the thermophysical parameters and the total temperature field in the current spatial domain;

[0021] The sixth module is used to perform inverse Fourier transform on the abnormal temperature field in the spatial wave number mixed domain to obtain the abnormal field temperature in the spatial domain;

[0022] The seventh module is used to obtain a new total temperature field in the space domain based on the background field temperature in the space domain and the abnormal field temperature in the space domain;

[0023] The eighth module is used to set the iterative convergence conditions and determine whether the new spatial domain temperature total field meets the iterative convergence conditions. If so, the new spatial domain temperature total field is output; otherwise, the new spatial domain temperature total field is used as the current spatial domain temperature total field in the next iteration and input into the fifth module.

[0024] On the other hand, the present invention provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the steps in the forward modeling method of the three-dimensional steady-state heat conduction geothermal field are implemented.

[0025] On the other hand, the present invention further provides a computer-readable storage medium having a computer program stored thereon, which implements the steps in the forward modeling method of the three-dimensional steady-state heat conduction geothermal field when the computer program is executed by a processor.

[0026] Compared with the prior art, the advantages of the present invention are:

[0027] The present invention determines the model by setting the geothermal field calculation range and the shape, size and position of the anomaly, and gives the thermophysical property parameters of the model; uses two-dimensional Fourier transform to transform the three-dimensional problem into a one-dimensional problem, thereby realizing efficient and high-precision numerical simulation of the three-dimensional steady-state heat conduction geothermal field.

[0028] The present invention solves the problems of large computational complexity and high storage requirements of conventional methods in current large-scale numerical simulation of three-dimensional geothermal fields, and provides an effective forward modeling solution for realizing the fine inversion and interpretation of geothermal exploration measured data under complex field terrain or underground complex structures. BRIEF DESCRIPTION OF THE DRAWINGS

[0029] Figure 1 is a flow chart of an embodiment of the present invention;

[0030] Figure 2 is a schematic diagram of a calculation area and a model in one embodiment of the present invention;

[0031] Figure 3 A comparison diagram of calculation results using the method provided by the present invention and the traditional COMSOL Multiphysics method according to an embodiment of the present invention;

[0032] Figure 4 A relative error diagram of calculation results using the method provided by the present invention and the traditional COMSOL Multiphysics method according to an embodiment of the present invention;

[0033] Figure 5 It is a structural schematic diagram of an embodiment of the present invention. DETAILED DESCRIPTION

[0034] In order to make the purpose, technical scheme and advantages of the embodiments of the present invention more clearly understood, the following will be used with drawings and detailed descriptions to clearly illustrate the spirit of the contents disclosed by the present invention. After understanding the embodiments of the contents of the present invention, any person skilled in the art can change and modify the techniques taught by the contents of the present invention without departing from the spirit and scope of the contents of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention, but are not intended to limit the present invention.

[0035] Reference Figure 1 In one embodiment of the present invention, a forward modeling method for a three-dimensional steady-state heat conduction geothermal field is provided, comprising:

[0036] S1. Determine a three-dimensional anomaly, determine a target area containing the three-dimensional anomaly, and construct a three-dimensional prism model of the target area;

[0037] S2. Meshing the three-dimensional prism model along the x, y, and z directions to obtain a series of mesh nodes and meshing parameters of the three-dimensional prism model;

[0038] S3. assigning values ​​of the thermophysical parameters of each subdivision node in the three-dimensional prism model according to the distribution of the thermophysical parameters of the target area;

[0039] S4. Load boundary conditions, calculate the space domain background field temperature, and use it as the initial space domain temperature total field;

[0040] S5. Calculate the abnormal temperature field in the spatial wavenumber mixed domain based on the thermophysical parameters and the current total temperature field in the spatial domain;

[0041] S6. Perform inverse Fourier transform on the abnormal temperature field in the spatial wavenumber mixed domain to obtain the abnormal field temperature in the spatial domain;

[0042] S7. based on the space domain background field temperature and the space domain abnormal field temperature, obtain a new space domain total temperature field;

[0043] S8. Set the iteration convergence condition. If the new spatial domain temperature total field meets the iteration convergence condition, output the new spatial domain temperature total field. Otherwise, use the new spatial domain temperature total field as the current spatial domain temperature total field in the next iteration and return to S5.

[0044] The thermophysical parameters of the present invention include background thermal conductivity, background heat generation rate, abnormal thermal conductivity, and abnormal heat generation rate. In step S3, the background thermal conductivity, background heat generation rate, abnormal thermal conductivity, and abnormal heat generation rate of each dissection node are assigned values.

[0045] The boundary conditions described in the present invention are boundary conditions of the geothermal field, which are generally classified into three categories, namely, boundary conditions of given temperature, boundary conditions of given heat flow, and boundary conditions of given heat exchange coefficient. In S4 of the present invention, any of the boundary conditions of the geothermal field can be loaded. The boundary condition of given temperature means that the temperature of the subdivided nodes distributed on the boundary of the target area is known, the boundary condition of given heat flow means that the heat flux density of the subdivided nodes distributed on the boundary of the target area is known, and the boundary condition of given heat exchange coefficient means that the temperature of the heat source and the heat exchange coefficient on the boundary of the target area are known.

[0046] In S4 of the present invention, a one-dimensional finite element method based on quadratic function interpolation is used to calculate the background field temperature in the space domain, and the calculation formula is:

[0047]

[0048] In the formula, λ0(x i ,y j ,z k ) is the coordinate position (x i ,y j ,z k ) of the mesh node, T0(x i ,y j ,z k ) is the coordinate position (x i ,y j ,z k ) of the spatial background field temperature of the node, Q0(x i ,y j ,z k ) is the coordinate position (x i ,y j ,z k ) is the background heat generation rate of the meshed nodes. In each element, the function is assumed to be a quadratic function.

[0049] In S5 of the present invention, the abnormal field temperature in the spatial domain, the abnormal heat flux density in the spatial domain, and the abnormal heat generation rate in the spatial domain are subjected to a two-dimensional Fourier transform to obtain the abnormal field temperature in the spatial wave number mixing domain Abnormal heat flux in spatial wavenumber mixing domain Abnormal heat generation rate in spatial wave number mixing domain

[0050]

[0051]

[0052] Where: i=1,2,...N x ,j=1,2,...N y ,k=1,2,...N k , N x 、N y 、N k are the number of nodes obtained after meshing the three-dimensional prism model along the x, y, and z directions; T a (x i ,y j ,z k ) is the coordinate position (x i ,y j ,z k ) of the spatial domain abnormal field temperature of the node, q a (x i ,y j ,z k ) is the coordinate position (xi ,y j ,z k ) of the node in the spatial domain, Q a (x i ,y j ,z k ) is the coordinate position (x i ,y j ,z k ) of the node in the spatial domain, k x , k y is the discrete offset wave number in the x and y directions, Δx and Δy are the sampling intervals when the three-dimensional prism model is meshed along the x and y directions, e represents a natural constant, represents the imaginary unit, (k x ,k y ,z k ) is the coordinate position (x i ,y j ,z k )'s corresponding coordinate position in the spatial wave number mixing domain;

[0053] Calculate the temperature of the abnormal field in the spatial wavenumber mixed domain When , the finite element method based on quadratic function interpolation is used to solve, and the calculation formula is:

[0054]

[0055]

[0056]

[0057]

[0058] Where, T (n-1) It represents the total temperature field of the current spatial domain corresponding to the nth iteration. In the first iteration, the spatial domain background field temperature T0 is used as the total temperature field T of the current spatial domain corresponding to the first iteration. (0) , i is the imaginary unit, They are the anomalous heat flux in the spatial wavenumber mixing domain The three components in the x, y, and z directions, λ a (x i ,y j ,z k ) is the coordinate position (x i ,y j ,z k )’s spatial domain anomalous thermal conductivity of the meshed nodes.

[0059] Next, the temperature of the anomalous field in the spatial wavenumber mixing domain is analyzed. Perform two-dimensional inverse Fourier transform to obtain the abnormal temperature field T in the spatial domain a (x i ,y j ,z k ).

[0060] In the nth iteration, the sum of the spatial domain background field temperature and the spatial domain abnormal field temperature calculated in the nth iteration is the new spatial domain total temperature field:

[0061] T (n) (x i ,y j ,z k )=T0(x i ,y j ,z k )+T a (x i ,y j ,z k )

[0062] In another embodiment of the present invention, k x , k y The calculation method is as follows:

[0063] The number of Gaussian points is selected as 4, and the node coordinate t is:

[0064]

[0065] Then the discrete offset wave number k in the x and y directions is x , k y for:

[0066]

[0067]

[0068] In the formula, Δk x ,Δk y are the fundamental wave numbers in the x and y directions respectively,

[0069] When N x 、N y When is an even number,

[0070]

[0071]

[0072] When N x 、N y When is an odd number,

[0073]

[0074]

[0075] It can be understood that the preset iteration termination condition refers to a preset model calculation constraint condition, which is used to constrain the entire model to converge in the process of performance calculation so that the model can output a result that meets the condition. In the present invention, the iteration termination condition in (S8) can be set to:

[0076]

[0077] Where: T (n-1) represents the total temperature field in the current spatial domain corresponding to the nth iteration, T (n) It represents the new total temperature field in the spatial domain obtained by the nth iteration calculation, and ε0 is the set numerical accuracy.

[0078] When the above iterative convergence conditions are met, the iteration stops. Of course, in practical applications, those skilled in the art may also set other iterative termination conditions based on the prior art, conventional technical means in the field or common knowledge, which are not limited to the iterative termination conditions set in the above preferred embodiments of the present application.

[0079] The accuracy and efficiency of the forward modeling method of the three-dimensional steady-state heat conduction geothermal field provided by the present invention are tested below.

[0080] The method of this embodiment is implemented by programming in Fortran language, and the test computer configuration is i5-4590, main frequency is 3.30GHz, and memory is 16GB.

[0081] A three-dimensional thermal conductivity model is designed with a size of 10km×10km×10km and a thermal conductivity of 2W / (m·℃). An anomaly is set in the center of the model with a size of 2km×2km×2km and a thermal conductivity of 6W / (m·℃). The projection of the center of the anomaly on the ground is taken as the coordinate origin. The model schematic diagram is shown in the figure below. Figure 2 The grid nodes are 51×51×51, the upper boundary adopts the first type of boundary conditions, the temperature value is 10℃, the lower boundary adopts the second type of boundary conditions, and the heat flux value is 41.86mW / m 2 In the single-threaded case, the Fourier transform of four Gaussian points is taken, and the expected numerical accuracy of the iteration is 10 -4 . The algorithm of the present invention and COMSOL Multiphysics software are used for calculation respectively. Figure 3 This is a comparison chart of the calculation results of this algorithm and COMSOL Multiphysics on the x=0, y=0 measurement line. Figure 4The relative error diagram between the two algorithms shows that the relative error of the calculation results of the two algorithms is less than 0.2%, which shows the correctness of the method of this embodiment. Under the same grid conditions, the algorithm of the present invention takes 0.5s and occupies 0.05GB of memory, while COMSOL Multiphysics takes 31s and occupies 4.2GB of memory, which fully shows the efficiency of the method of this embodiment.

[0082] In one embodiment of the present invention, a forward modeling device for a three-dimensional steady-state heat conduction geothermal field is provided, comprising:

[0083] The first module is used to determine a three-dimensional anomaly, determine a target area containing the three-dimensional anomaly, and construct a three-dimensional prism model of the target area;

[0084] The second module is used to mesh the three-dimensional prism model along the x, y, and z directions to obtain a series of mesh nodes and mesh parameters of the three-dimensional prism model;

[0085] The third module is used to assign values ​​of the thermophysical parameters of each segmentation node in the three-dimensional prism model according to the distribution of the thermophysical parameters of the target area;

[0086] The fourth module is used to load boundary conditions and calculate the background field temperature of the space domain, and use it as the total temperature field of the initial space domain;

[0087] The fifth module is used to calculate the abnormal temperature field in the spatial wavenumber mixed domain based on the thermophysical parameters and the total temperature field in the current spatial domain;

[0088] The sixth module is used to perform inverse Fourier transform on the abnormal temperature field in the spatial wave number mixed domain to obtain the abnormal field temperature in the spatial domain;

[0089] The seventh module is used to obtain a new total temperature field in the space domain based on the background field temperature in the space domain and the abnormal field temperature in the space domain;

[0090] The eighth module is used to set the iterative convergence conditions and determine whether the new spatial domain temperature total field meets the iterative convergence conditions. If so, the new spatial domain temperature total field is output; otherwise, the new spatial domain temperature total field is used as the current spatial domain temperature total field in the next iteration and input into the fifth module.

[0091] The methods for implementing the functions of the above modules can be implemented using the same methods as in the above embodiments, which will not be described in detail here.

[0092] In this embodiment, a computer device is provided. The computer device may be a server, and its internal structure diagram may be as shown in FIG. Figure 5As shown. The computer device includes a processor, a memory, a network interface and a database connected through a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store sample data. The network interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, the forward modeling method of the above-mentioned three-dimensional steady-state heat conduction geothermal field is implemented.

[0093] Those skilled in the art will understand that Figure 5 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0094] In one embodiment, a computer device is provided, including a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps of the forward modeling method of the three-dimensional steady-state heat conduction geothermal field in the above embodiment when executing the computer program.

[0095] In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the steps of the forward modeling method of the three-dimensional steady-state heat conduction geothermal field in the above embodiment are implemented.

[0096] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).

[0097] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0098] The above-mentioned embodiments only express several implementation methods of the present application, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the invention patent. It should be pointed out that, for a person of ordinary skill in the art, several variations and improvements can be made without departing from the concept of the present application, and these all belong to the protection scope of the present application. Therefore, the protection scope of the patent of the present application shall be subject to the attached claims.

Claims

1. A forward modeling method for three-dimensional steady-state heat conduction geothermal field, characterized in that: include: S1. Determine a three-dimensional anomaly, determine a target area containing the three-dimensional anomaly, and construct a three-dimensional prism model of the target area; S2. Meshing the three-dimensional prism model along the x, y, and z directions to obtain a series of mesh nodes and meshing parameters of the three-dimensional prism model; S3. assigning values ​​of the thermophysical parameters of each subdivision node in the three-dimensional prism model according to the distribution of the thermophysical parameters of the target area; S4. Load boundary conditions, calculate the space domain background field temperature, and use it as the initial space domain temperature total field; S5. Calculate the abnormal temperature field in the spatial wavenumber mixed domain based on the thermophysical parameters and the total temperature field in the current spatial domain; S6. Perform inverse Fourier transform on the abnormal temperature field in the spatial wavenumber mixed domain to obtain the abnormal field temperature in the spatial domain; S7. based on the space domain background field temperature and the space domain abnormal field temperature, obtain a new space domain total temperature field; S8. Set the iteration convergence condition. If the new spatial domain temperature total field meets the iteration convergence condition, output the new spatial domain temperature total field. Otherwise, use the new spatial domain temperature total field as the current spatial domain temperature total field in the next iteration and return to S5.

2. The forward modeling method of three-dimensional steady-state heat conduction geothermal field according to claim 1 is characterized in that: Thermophysical parameters include background thermal conductivity, background heat generation rate, abnormal thermal conductivity and abnormal heat generation rate. The background thermal conductivity, background heat generation rate, abnormal thermal conductivity and abnormal heat generation rate are assigned to each dissection node.

3. The forward modeling method of three-dimensional steady-state heat conduction geothermal field according to claim 1 or 2, characterized in that: The boundary condition is any one of a boundary condition of a given temperature, a boundary condition of a given heat flow, and a boundary condition of a given heat exchange coefficient, wherein the boundary condition of a given temperature means that the temperature of the segmentation nodes distributed on the boundary of the target area is known, the boundary condition of a given heat flow means that the heat flux density of the segmentation nodes distributed on the boundary of the target area is known, and the boundary condition of a given heat exchange coefficient means that the heat source temperature and the heat exchange coefficient on the boundary of the target area are known.

4. The forward modeling method of three-dimensional steady-state heat conduction geothermal field according to claim 3 is characterized in that: The one-dimensional finite element method based on quadratic function interpolation is used to calculate the background field temperature in the space domain. The calculation formula is: In the formula, λ0(x i ,y j ,z k ) is the coordinate position (x i ,y j ,z k ) of the mesh node, T0(x i ,y j ,z k ) is the coordinate position (x i ,y j ,z k ) of the spatial background field temperature of the node, Q0(x i ,y j ,z k ) is the coordinate position (x i ,y j ,z k )’s background heat generation rate at the meshed nodes.

5. The forward modeling method of three-dimensional steady-state heat conduction geothermal field according to claim 3 is characterized in that: S5 includes: The abnormal field temperature in the spatial domain, the abnormal heat flux density in the spatial domain, and the abnormal heat generation rate in the spatial domain are transformed into the abnormal field temperature in the spatial wavenumber mixing domain. Abnormal heat flux in the spatial wavenumber mixing domain Abnormal heat generation rate in spatial wave number mixing domain Where: i=1,2,...N x ,j=1,2,...N y ,k=1,2,...N k , N x 、N y 、N k are the number of nodes obtained after meshing the three-dimensional prism model along the x, y, and z directions; T a (x i ,y j ,z k ) is the coordinate position (x i ,y j ,z k ) of the spatial domain abnormal field temperature of the node, q a (x i ,y j ,z k ) is the coordinate position (x i ,y j ,z k ) of the node in the spatial domain, Q a (x i ,y j ,z k ) is the coordinate position (x i ,y j ,z k ) of the node in the spatial domain, k x , k y is the discrete offset wave number in the x and y directions, Δx and Δy are the sampling intervals when the three-dimensional prism model is meshed along the x and y directions, respectively, e represents a natural constant, represents the imaginary unit, (k x ,k y ,z k ) is the coordinate position (x i ,y j ,z k )'s corresponding coordinate position in the spatial wavenumber mixing domain; Calculate the temperature of the abnormal field in the spatial wavenumber mixed domain When , the finite element method based on quadratic function interpolation is used to solve it, and the calculation formula is: Where, T (n-1) It represents the total temperature field of the current spatial domain corresponding to the nth iteration. In the first iteration, the spatial domain background field temperature T0 is used as the total temperature field T of the current spatial domain corresponding to the first iteration. (0) , i is the imaginary unit, They are the anomalous heat flux in the spatial wavenumber mixing domain The three components in the x, y, and z directions, λ a (x i ,yj,z k ) is the coordinate position (x i ,y j ,z k )’s spatial domain anomalous thermal conductivity of the meshed nodes.

6. The forward modeling method of three-dimensional steady-state heat conduction geothermal field according to claim 5 is characterized in that: k x , k y The calculation method is as follows: The number of Gaussian points is selected as 4, and the node coordinate t is: Then the discrete offset wave number k in the x and y directions is x , k y for: In the formula, Δk x ,Δk y are the fundamental wave numbers in the x and y directions respectively, When N x 、N y When is an even number, When N x 、N y When is an odd number, 7. The forward modeling method of three-dimensional steady-state heat conduction geothermal field according to claim 1, 2, 4, 5 or 6, characterized in that: The set iterative convergence condition is: Where: T (n-1) represents the total temperature field in the current spatial domain corresponding to the nth iteration, T (n) It represents the new total temperature field in the spatial domain obtained by the nth iteration calculation, and ε0 is the set numerical accuracy.

8. A forward modeling device for three-dimensional steady-state heat conduction geothermal field, characterized in that: include: The first module is used to determine a three-dimensional anomaly, determine a target area containing the three-dimensional anomaly, and construct a three-dimensional prism model of the target area; The second module is used to mesh the three-dimensional prism model along the x, y, and z directions to obtain a series of mesh nodes and mesh parameters of the three-dimensional prism model; The third module is used to assign values ​​of the thermophysical parameters of each segmentation node in the three-dimensional prism model according to the distribution of the thermophysical parameters of the target area; The fourth module is used to load boundary conditions and calculate the background field temperature of the space domain, and use it as the total temperature field of the initial space domain; The fifth module is used to calculate the abnormal temperature field in the spatial wavenumber mixed domain based on the thermophysical parameters and the total temperature field in the current spatial domain; The sixth module is used to perform inverse Fourier transform on the abnormal temperature field in the spatial wave number mixed domain to obtain the abnormal field temperature in the spatial domain; The seventh module is used to obtain a new total temperature field in the space domain based on the background field temperature in the space domain and the abnormal field temperature in the space domain; The eighth module is used to set the iterative convergence conditions and determine whether the new spatial domain temperature total field meets the iterative convergence conditions. If so, the new spatial domain temperature total field is output; otherwise, the new spatial domain temperature total field is used as the current spatial domain temperature total field in the next iteration and input into the fifth module.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps in the forward modeling method of the three-dimensional steady-state heat conduction geothermal field described in claim 1 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps in the forward modeling method of the three-dimensional steady-state heat conduction geothermal field described in claim 1 are implemented.

Images