CN114036806A - Three-dimensional geothermal field numerical simulation method based on thermal conductivity anisotropic medium - Google Patents
Three-dimensional geothermal field numerical simulation method based on thermal conductivity anisotropic medium Download PDFInfo
- Publication number
- CN114036806A CN114036806A CN202111423020.6A CN202111423020A CN114036806A CN 114036806 A CN114036806 A CN 114036806A CN 202111423020 A CN202111423020 A CN 202111423020A CN 114036806 A CN114036806 A CN 114036806A
- Authority
- CN
- China
- Prior art keywords
- field
- temperature
- domain
- dimensional
- abnormal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000004088 simulation Methods 0.000 title claims abstract description 35
- 230000002159 abnormal effect Effects 0.000 claims abstract description 75
- 238000004364 calculation method Methods 0.000 claims description 32
- 230000020169 heat generation Effects 0.000 claims description 23
- 230000004907 flux Effects 0.000 claims description 9
- 238000012887 quadratic function Methods 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 2
- 238000004590 computer program Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000002547 anomalous effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The three-dimensional geothermal field numerical simulation method based on the thermal conductivity anisotropic medium comprises the steps of constructing a three-dimensional prism model containing a three-dimensional abnormal body inside, carrying out grid subdivision on the three-dimensional prism model, and assigning values to thermal physical parameters of each subdivision node in the three-dimensional prism model according to distribution conditions of the thermal physical parameters, wherein the thermal conductivity is tensor; calculating the ambient field temperature of the spatial domain, and taking the ambient field temperature as an initial spatial domain temperature total field; calculating an abnormal temperature field of the spatial wave number mixed domain based on the thermophysical parameters and the current total spatial domain temperature field; performing two-dimensional inverse Fourier transform on the spatial wave number mixed domain abnormal temperature field to obtain the temperature of the spatial domain abnormal field; obtaining a new total spatial domain temperature field based on the spatial domain background field temperature and the spatial domain abnormal field temperature; and continuously iterating until the new space domain temperature total field obtained by iteration meets the iteration convergence condition. The invention realizes the high-efficiency and high-precision numerical simulation of the three-dimensional geothermal field of the thermal conductivity anisotropic medium.
Description
Technical Field
The invention belongs to the technical field of geothermal field numerical simulation, and particularly relates to a three-dimensional geothermal field numerical simulation method based on a thermal conductivity anisotropic medium.
Background
Geothermal energy is used as a green energy source and has wider and wider application prospect. The research on the numerical simulation of the geothermal field has important economic value and social significance on the development of the geothermal theory and the exploitation and utilization of geothermal energy.
The current numerical simulation of the ground temperature field has the following defects: a conventional geothermal field numerical simulation method (such as a finite element method, a finite difference method, a boundary element method and a finite volume method) is used for calculating in a spatial domain, a model needs to be finely divided when a high precision requirement is met in a large-scale complex medium, and the number of dividing units is large, so that the storage capacity is large and the calculation time is long.
Under the actual complex geological conditions, on one hand, the three-dimensional geothermal exploration method and the geothermal field numerical simulation method used in the exploration system have the problems of large calculation amount and low calculation efficiency, and on the other hand, the thermal conductivity of the rock mostly has anisotropy, and the anisotropy of the thermal conductivity needs to be considered in the geothermal field numerical simulation. Therefore, it is necessary to provide a numerical simulation method of a three-dimensional geothermal field from the anisotropy point of view to solve the problems in the prior art.
Disclosure of Invention
The invention aims to provide a three-dimensional geothermal field numerical simulation method based on a thermal conductivity anisotropic medium, aiming at the current situations of large calculation amount and high storage requirement when the conventional methods such as the current finite element method, the finite difference method, the boundary element method and the like are used for processing large-scale numerical simulation of a geothermal field.
In order to achieve the technical purpose, the technical scheme provided by the invention is as follows:
in one aspect, the invention provides a three-dimensional geothermal field numerical simulation method based on a thermal conductivity anisotropic medium, comprising the following steps:
s1, determining a three-dimensional abnormal body, determining a target area containing the three-dimensional abnormal body, and constructing a three-dimensional prism model of the target area;
s2, carrying out mesh subdivision on the three-dimensional prism model along the directions of x, y and z, and subdividing to obtain a series of subdivision nodes and mesh subdivision parameters of the three-dimensional prism model;
s3, giving the thermal physical parameters of each subdivision node in the three-dimensional prism model according to the distribution condition of the thermal physical parameters of the target area, wherein the thermal physical parameters comprise thermal conductivity and heat generation rate, and the thermal conductivity is tensor;
s4, calculating the background field temperature of the spatial domain, and taking the background field temperature as an initial spatial domain temperature total field;
s5, calculating an abnormal temperature field of the spatial wave number mixed domain based on the thermal physical parameters and the current total spatial domain temperature field;
s6, performing inverse Fourier transform on the spatial wave number mixed domain abnormal temperature field to obtain the spatial domain abnormal field temperature;
s7, obtaining a new total spatial domain temperature field based on the spatial domain background field temperature and the spatial domain abnormal field temperature;
and S8, setting an iteration convergence condition, outputting a new space domain temperature total field if the new space domain temperature total field meets the iteration convergence condition, and returning to S5 if the new space domain temperature total field is not used as the current space domain temperature total field in the next iteration.
In another aspect, the present invention provides a three-dimensional geothermal field numerical simulation apparatus based on a thermal conductivity anisotropic medium, including:
the system comprises a first module, a second module and a third module, wherein the first module is used for determining a three-dimensional abnormal body, determining a target area containing the three-dimensional abnormal body inside and constructing a three-dimensional prism model of the target area;
the second module is used for carrying out mesh subdivision on the three-dimensional prism model along the directions of x, y and z to obtain a series of subdivision nodes and mesh subdivision parameters of the three-dimensional prism model through subdivision;
the third module is used for giving the thermal physical parameters of each subdivision node in the three-dimensional prism model according to the distribution condition of the thermal physical parameters of the target area, wherein the thermal physical parameters comprise thermal conductivity and heat generation rate, and the thermal conductivity is tensor;
the fourth module is used for calculating the ambient field temperature of the spatial domain and taking the ambient field temperature as the total temperature field of the initial spatial domain;
a fifth module, configured to calculate a spatial wave number mixed domain abnormal temperature field based on the thermophysical parameters and the current spatial domain total temperature field;
the sixth module is used for performing inverse Fourier transform on the spatial wave number mixed domain abnormal temperature field to obtain the spatial domain abnormal field temperature;
the seventh module is used for obtaining a new total temperature field of the space domain based on the background field temperature of the space domain and the abnormal field temperature of the space domain;
an eighth module, configured to set an iterative convergence condition, determine whether the new spatial domain total temperature field meets the iterative convergence condition, and if yes, output the new spatial domain total temperature field; and otherwise, taking the new total spatial domain temperature field as the current total spatial domain temperature field in the next iteration, and inputting the current total spatial domain temperature field into the fifth module.
In another aspect, the present invention provides a computer device, which includes a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps in the three-dimensional geothermal field numerical simulation method based on a thermal conductivity anisotropic medium when executing the computer program.
In still another aspect, the present invention further provides a computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, implements the steps in the method for numerical simulation of a three-dimensional geothermal field based on a thermal conductivity anisotropic medium.
Compared with the prior art, the invention has the advantages that:
the invention provides a high-efficiency and high-precision numerical simulation method for a geothermal field with any shape and any thermal conductivity distribution, and the thermal conductivity of rocks has anisotropy in consideration of the actual complex situation underground. For large-scale numerical simulation of a three-dimensional geothermal field of a thermal conductivity anisotropic medium, a method with high calculation speed and small occupied memory is provided, and compared with the conventional common mature software COMSOL Multiphysics, the method has remarkable efficiency advantage.
Drawings
FIG. 1 is a flow chart in one embodiment of the present invention;
FIG. 2 is a schematic diagram of a calculation region and a model according to an embodiment of the invention;
fig. 3 is a comparison graph of the calculation results of the method provided by the present invention and the conventional COMSOL Multiphysics method according to an embodiment of the present invention, wherein (a) is a graph of the calculation results of the conventional COMSOL Multiphysics software on the y-0 plane, and (b) is a graph of the calculation results of the method of the present invention on the y-0 plane;
FIG. 4 is a diagram illustrating relative error of the calculation results of the method of the present invention and the conventional COMSOL Multiphysics method according to an embodiment of the present invention;
fig. 5 is a schematic structural diagram according to an embodiment of the present invention.
Detailed Description
For the purpose of promoting a clear understanding of the objects, aspects and advantages of the embodiments of the invention, reference will now be made to the drawings and detailed description, wherein there are shown in the drawings and described below specific embodiments of the invention, in which modifications and variations can be made by one skilled in the art without departing from the spirit and scope of the invention. The exemplary embodiments of the present invention and the description thereof are provided to explain the present invention and not to limit the present invention.
Referring to fig. 1, in an embodiment of the present invention, a three-dimensional geothermal field numerical simulation method based on a thermal conductivity anisotropic medium is provided, including:
s1, determining a three-dimensional abnormal body, determining a target area containing the three-dimensional abnormal body, and constructing a three-dimensional prism model of the target area;
s2, carrying out mesh subdivision on the three-dimensional prism model along the directions of x, y and z, and subdividing to obtain a series of subdivision nodes and mesh subdivision parameters of the three-dimensional prism model;
s3, giving the thermal physical parameters of each subdivision node in the three-dimensional prism model according to the distribution condition of the thermal physical parameters of the target area, wherein the thermal physical parameters comprise thermal conductivity and heat generation rate, and the thermal conductivity is tensor;
s4, loading boundary conditions, calculating the background field temperature of the spatial domain, and taking the background field temperature as an initial spatial domain temperature total field;
s5, calculating an abnormal temperature field of the spatial wave number mixed domain based on the thermal physical parameters and the current total spatial domain temperature field;
s6, performing inverse Fourier transform on the spatial wave number mixed domain abnormal temperature field to obtain the spatial domain abnormal field temperature;
s7, obtaining a new total spatial domain temperature field based on the spatial domain background field temperature and the spatial domain abnormal field temperature;
and S8, setting an iteration convergence condition, outputting a new space domain temperature total field if the new space domain temperature total field meets the iteration convergence condition, and returning to S5 if the new space domain temperature total field is not used as the current space domain temperature total field in the next iteration.
The thermophysical parameters include thermal conductivity and heat generation, and the thermal conductivity and the heat generation of each subdivision node are assigned in step S3. The heat conductivity of each subdivision node is composed of background heat conductivity and abnormal heat conductivity, the heat generation rate of each subdivision node is composed of background heat generation rate and abnormal heat generation rate, and the background heat conductivity, the background heat generation rate, the abnormal heat conductivity and the abnormal heat generation rate of each subdivision node are assigned. Considering the thermal conductivity anisotropy, the thermal conductivity of each subdivision node is tensor, and the thermal conductivity of each subdivision node has 9 components. Considering that most media have symmetrical properties, the matrix form of the thermal conductivity of each subdivision node can be expressed as:
wherein: lambda [ alpha ]0(xi,yj,zk) Is set as (x) at the coordinate positioni,yj,zk) Background thermal conductivity, λ, of the subdivision node of (1)a(xi,yj,zk) Is set as (x) at the coordinate positioni,yj,zk) The abnormal thermal conductivity of the subdivision node of (x)i,yj,zk) N. N denotes the coordinate position of the subdivision node numbered i, j, k, i ═ 1,2x，j＝1,2,...Ny，k＝1,2,...Nk，Nx、Ny、NzAnd respectively obtaining the number of subdivision nodes after the three-dimensional prism model carries out mesh subdivision along the directions of x, y and z.
The boundary conditions of the present invention are boundary conditions of a ground temperature field, and generally include three types, namely, boundary conditions of a given temperature, boundary conditions of a given heat flow, and boundary conditions of a given heat exchange coefficient. In the present invention S4, any of the above boundary conditions of the earth temperature field may be loaded. The boundary condition of the given temperature means that the temperature of the subdivision nodes distributed on the boundary of the target area is known, the boundary condition of the given heat flow means that the heat flow density of the subdivision nodes distributed on the boundary of the target area is known, and the boundary condition of the given heat exchange coefficient means that the temperature of a heat source and the heat exchange coefficient on the boundary of the target area are known.
In the invention S4, a one-dimensional finite element method based on quadratic function interpolation is adopted to calculate the background field temperature of each subdivision node space domain, and the calculation formula is as follows:
in the formula of0(xi,yj,zk) Is set as (x) at the coordinate positioni,yj,zk) Background thermal conductivity, T, of the subdivision node0(xi,yj,zk) Is set as (x) at the coordinate positioni,yj,zk) The spatial domain background field temperature, Q, of the subdivision node0(xi,yj,zk) Is set as (x) at the coordinate positioni,yj,zk) The background heat generation rate of the subdivision nodes;
for T0(xi,yj,zk) Respectively making derivation in x, y and z directions to obtain three components T in x, y and z directions0x(xi,yj,zk)、T0y(xi,yj,zk)、T0z(xi,yj,zk)。
In an embodiment of the present invention, S5 includes:
s5.1, calculating the spatial domain abnormal heat flow density on each subdivision node based on the abnormal heat conductivity and the current spatial domain temperature total field;
s5.2, performing two-dimensional Fourier transform on the spatial domain abnormal heat flux density on each subdivision node, and calculating to obtain the spatial wave number mixed domain abnormal heat flux density on each subdivision node;
and S5.3, calculating the abnormal field temperature of the spatial wave number mixed domain based on the abnormal heat flow density of the spatial wave number mixed domain.
In an embodiment of the present invention, an implementation method of S5.1 includes:
for the coordinate position is (x)i,yj,zk) The spatial domain abnormal heat flux density q on the subdivision nodea(xi,yj,zk)，qax(xi,yj,zk)、qay(xi,yj,zk)、qaz(xi,yj,zk) Is qa(xi,yj,zk) The three components in the x, y, z directions are:
wherein: t is(n)(xi,yj,zk) Is a total field of the temperature of the current space domain corresponding to the nth iteration, and the background field temperature T of the space domain during the first iteration0(xi,yj,zk) As a current total field of spatial domain temperature, λ, corresponding to the first iterationa(xi,yj,zk) Is set as (x) at the coordinate positioni,yj,zk) The abnormal thermal conductivity of the subdivision nodes.
The total field T of the temperature of the current space domain corresponding to the first iteration(0)(xi,yj,zk) Three components in the x, y, z directions are
Namely, it is
In an embodiment of the present invention, an implementation method of S5.2 includes: the coordinate position of the spatial wave number mixed domain is (k)x,ky,zk) The spatial wave number mixed domain abnormal heat flux density on the subdivision node
The three components in the x, y, z directions are:
wherein: k is a radical ofx，kyDiscrete offset wavenumbers in the x, y directions, e represents a natural constant,
denotes the unit of imaginary number, (k)x,ky,zk) Is set as (x) at the coordinate positioni,yj,zk) And (4) corresponding coordinate positions of the subdivision nodes in the space wave number mixed domain.
In an embodiment of the present invention, an implementation method of S5.3 includes:
spatial wavenumber mixed domain anomalous field temperature
The following governing equation is satisfied:
in the formula,
kx，kydiscrete offset wavenumbers in the x, y directions;
indicates a coordinate position of (x)i,yj,zk) The abnormal heat generation rate of the space wave number mixed domain of the subdivision node is (x) from the coordinate positioni,yj,zk) Abnormal heat generation rate Q of the divided nodea(xi,yj,zk) Performing two-dimensional Fourier transform to obtain the Fourier transform; i is an imaginary unit;
solving the control equation by using a finite element method based on quadratic function interpolation to obtain the abnormal field temperature of the spatial wave number mixed domain
For abnormal field temperature in space wave number domain
Performing two-dimensional Fourier inverse transformation to obtain the temperature T of the spatial domain abnormal fielda(xi,yj,zk) I.e. by
The spatial domain abnormal field temperature T obtained by the current iteration calculationa(xi,yj,zk) And the spatial domain background temperature field T0(xi,yj,zk) And the sum is the new total field of the spatial domain temperature obtained by the current iterative calculation.
In another embodiment of the present invention, kx，kyThe calculation method of (2) is as follows:
kx，kythe calculation method of (2) is as follows:
kx＝p·Δkx
ky＝q·Δky
wherein,
delta x and delta y are unit interval lengths between adjacent subdivision nodes when the three-dimensional prism model carries out grid subdivision along the x direction and the y direction respectively, and Nx、NyRespectively performing mesh subdivision on the three-dimensional prism model and then dividing the node number in the x and y directions;
when N is presentx、NyIn the case of an even number, the number of the first,
when N is presentx、NyIn the case of an odd number of the groups,
it is understood that the preset iteration termination condition refers to a preset model calculation constraint condition for constraining the whole model to converge in the performance calculation process, so that the model can output a result meeting the condition. In the present invention, the iteration termination condition in (S8) may be set to:
wherein: t is(n+1)(xi,yj,zk) Represents the new total temperature field of the space domain, epsilon, obtained by the nth iteration calculation0Abs () represents the absolute value for a set numerical precision;
if the iteration convergence condition is met, stopping iteration, outputting a new space domain temperature total field obtained by the nth iteration calculation, and if the iteration convergence condition is not met, taking the new space domain temperature total field obtained by the nth iteration calculation as a current space domain temperature total field corresponding to the (n + 1) th iteration, and continuing the iteration.
Of course, in practical applications, a person skilled in the art may set other iteration termination conditions based on the prior art, the conventional technical means in the field, or the common general knowledge, and is not limited to the iteration termination conditions set in the above preferred embodiments of the present application.
The accuracy and efficiency of the three-dimensional geothermal field numerical simulation method based on the thermal conductivity anisotropic medium provided by the invention are tested.
The method is realized by Fortran language programming, and the configuration of a testing computer is i9-7980XE, the main frequency is 2.60GHz and the internal memory is 64.0 GB.
Fig. 2 shows a schematic diagram of the calculation region and the three-dimensional prism model in this embodiment. The projection position of the center of the three-dimensional prism model on the ground is set as the origin of coordinates, the range of a calculation region in the x direction is-5 km, the range in the y direction is-5 km, the range in the z direction is 0-10 km, a background field is an isotropic medium, the thermal conductivity is 2W/(m.DEG C), an abnormal body is arranged in the center of the three-dimensional prism model, the abnormal body is an anisotropic medium, and the size is 2km multiplied by 2 km. Coefficient of thermal conductivity of
The first type of boundary conditions are adopted in consideration of the upper and lower boundaries, the upper boundary temperature value is 20 ℃, and the lower boundary temperature value is 100 ℃. The number of split nodes is 101 × 101 × 101. The calculation was performed using a standard FFT with an expected numerical precision of 10 to achieve the convergence criterion-4. The calculation is respectively carried out by using the method of the invention and the traditional COMSOL Multiphysics software, and fig. 3 is a comparison graph of the calculation results of the method of the invention and the traditional COMSOL Multiphysics on a y-0 measuring surface, wherein (a) is a calculation result graph of the traditional COMSOL Multiphysics software on the y-0 measuring surface, and (b) is a calculation result graph of the method of the invention on the y-0 measuring surface. FIG. 4 is a relative error chart between the method of the present invention and conventional COMSOL Multiphysics software, and it can be seen that the relative error of the calculation results of the two is below 1.0%, which illustrates the correctness of the method of the present embodiment. Under the same grid condition, the algorithm consumes 27.8s, occupies 0.76GB, the COMSOL multiprophy consumes 204.2s and occupies 21.8GB, and the high efficiency of the method is fully explained.
The invention provides a three-dimensional geothermal field numerical simulation device based on a thermal conductivity anisotropic medium in an embodiment, which comprises:
the system comprises a first module, a second module and a third module, wherein the first module is used for determining a three-dimensional abnormal body, determining a target area containing the three-dimensional abnormal body inside and constructing a three-dimensional prism model of the target area;
the second module is used for carrying out mesh subdivision on the three-dimensional prism model along the directions of x, y and z to obtain a series of subdivision nodes and mesh subdivision parameters of the three-dimensional prism model through subdivision;
the third module is used for giving the thermal physical parameters of each subdivision node in the three-dimensional prism model according to the distribution condition of the thermal physical parameters of the target area, wherein the thermal physical parameters comprise thermal conductivity and heat generation rate, and the thermal conductivity is tensor;
the fourth module is used for calculating the ambient field temperature of the spatial domain and taking the ambient field temperature as the total temperature field of the initial spatial domain;
a fifth module, configured to calculate a spatial wave number mixed domain abnormal temperature field based on the thermophysical parameters and the current spatial domain total temperature field;
the sixth module is used for performing inverse Fourier transform on the spatial wave number mixed domain abnormal temperature field to obtain the spatial domain abnormal field temperature;
the seventh module is used for obtaining a new total temperature field of the space domain based on the background field temperature of the space domain and the abnormal field temperature of the space domain;
an eighth module, configured to set an iterative convergence condition, determine whether the new spatial domain total temperature field meets the iterative convergence condition, and if yes, output the new spatial domain total temperature field; and otherwise, taking the new total spatial domain temperature field as the current total spatial domain temperature field in the next iteration, and inputting the current total spatial domain temperature field into the fifth module.
The implementation method of the functions of the modules can be implemented by the same method in the foregoing embodiments, and details are not repeated here.
In this embodiment, a computer device is provided, and the computer device may be a server, and its internal structure diagram may be as shown in fig. 5. The computer device includes a processor, a memory, a network interface, and a database connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device comprises a nonvolatile 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 an operating system and computer programs 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 for communicating with an external terminal through a network connection. The computer program is executed by a processor to realize the three-dimensional geothermal field numerical simulation method based on the thermal conductivity anisotropic medium.
Those skilled in the art will appreciate that the architecture shown in fig. 5 is merely a block diagram of some of the structures associated with the disclosed aspects and is not intended to limit the computing devices to which the disclosed aspects apply, as particular computing devices may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
In one embodiment, a computer device is provided, which includes a memory and a processor, the memory stores a computer program, and the processor implements the steps of the three-dimensional geothermal field numerical simulation method based on a thermal conductivity anisotropic medium in the above embodiments when executing the computer program.
In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored, which, when being executed by a processor, implements the steps of the thermal conductivity anisotropic medium-based three-dimensional geothermal field numerical simulation method in the above-described embodiments.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in the embodiments provided herein may include non-volatile and/or volatile memory, among others. 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. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), Dynamic RAM (DRAM), Synchronous DRAM (SDRAM), Double Data Rate SDRAM (DDRSDRAM), Enhanced SDRAM (ESDRAM), Synchronous Link DRAM (SLDRAM), Rambus Direct RAM (RDRAM), direct bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (10)
1. The three-dimensional geothermal field numerical simulation method based on the thermal conductivity anisotropic medium is characterized by comprising the following steps of:
s1, determining a three-dimensional abnormal body, determining a target area containing the three-dimensional abnormal body, and constructing a three-dimensional prism model of the target area;
s2, carrying out mesh subdivision on the three-dimensional prism model along the directions of x, y and z, and subdividing to obtain a series of subdivision nodes and mesh subdivision parameters of the three-dimensional prism model;
s3, giving the thermal physical parameters of each subdivision node in the three-dimensional prism model according to the distribution condition of the thermal physical parameters of the target area, wherein the thermal physical parameters comprise thermal conductivity and heat generation rate, and the thermal conductivity is tensor;
s4, calculating the background field temperature of the spatial domain, and taking the background field temperature as an initial spatial domain temperature total field;
s5, calculating an abnormal temperature field of the spatial wave number mixed domain based on the thermal physical parameters and the current total spatial domain temperature field;
s6, performing inverse Fourier transform on the spatial wave number mixed domain abnormal temperature field to obtain the spatial domain abnormal field temperature;
s7, obtaining a new total spatial domain temperature field based on the spatial domain background field temperature and the spatial domain abnormal field temperature;
and S8, setting an iteration convergence condition, outputting a new space domain temperature total field if the new space domain temperature total field meets the iteration convergence condition, and returning to S5 if the new space domain temperature total field is not used as the current space domain temperature total field in the next iteration.
2. The method for numerically simulating the three-dimensional geothermal field based on the thermal conductivity anisotropic medium according to claim 1, wherein the thermal conductivity of each subdivision node is composed of a background thermal conductivity and an abnormal thermal conductivity, the heat generation rate of each subdivision node is composed of a background heat generation rate and an abnormal heat generation rate, and the background thermal conductivity, the background heat generation rate, the abnormal thermal conductivity and the abnormal heat generation rate of each subdivision node are assigned, wherein the thermal conductivity of each subdivision node is tensor.
3. The method for numerical simulation of a three-dimensional geothermal field based on a thermally conductive anisotropic medium of claim 2, wherein the thermal conductivity of each subdivision node has 9 components, expressed as:
wherein: lambda [ alpha ]0(xi,yj,zk) Is set as (x) at the coordinate positioni,yj,zk) Background thermal conductivity, λ, of the subdivision node of (1)a(xi,yj,zk) Is set as (x) at the coordinate positioni,yj,zk) The abnormal thermal conductivity of the subdivision node of (x)i,yj,zk) N. N denotes the coordinate position of the subdivision node numbered i, j, k, i ═ 1,2x，j＝1,2,...Ny，k＝1,2,...Nk，Nx、Ny、NzAnd respectively obtaining the number of subdivision nodes after the three-dimensional prism model carries out mesh subdivision along the directions of x, y and z.
4. The three-dimensional geothermal field numerical simulation method based on the thermal conductivity anisotropic medium of claim 3, wherein the background field temperature of each subdivision node space domain is calculated by adopting a one-dimensional finite element method based on quadratic function interpolation, and the calculation formula is as follows:
in the formula of0(xi,yj,zk) Is set as (x) at the coordinate positioni,yj,zk) Background thermal conductivity, T, of the subdivision node0(xi,yj,zk) Is set as (x) at the coordinate positioni,yj,zk) The spatial domain background field temperature, Q, of the subdivision node0(xi,yj,zk) Is set as (x) at the coordinate positioni,yj,zk) The background heat generation rate of the subdivision nodes;
for T0(xi,yj,zk) Respectively making derivation in x, y and z directions to obtain three components T in x, y and z directions0x(xi,yj,zk)、T0y(xi,yj,zk)、T0z(xi,yj,zk)。
5. The method for numerical simulation of a three-dimensional geothermal field based on a thermal conductivity anisotropic medium according to claim 1,2, 3 or 4, wherein S5 comprises:
calculating the spatial domain abnormal heat flow density on each subdivision node based on the abnormal heat conductivity and the current spatial domain temperature total field;
performing two-dimensional Fourier transform on the spatial domain abnormal heat flux density on each subdivision node, and calculating to obtain the spatial wave number mixed domain abnormal heat flux density on each subdivision node;
and calculating the abnormal field temperature of the spatial wavenumber mixed domain based on the abnormal heat flux density of the spatial wavenumber mixed domain.
6. The method of numerical simulation of a three-dimensional geothermal field based on a thermally conductive anisotropic medium of claim 5, wherein the coordinate position is (x)i,yj,zk) The spatial domain abnormal heat flux density q on the subdivision nodea(xi,yj,zk)，qax(xi,yj,zk)、qay(xi,yj,zk)、qaz(xi,yj,zk) Is qa(xi,yj,zk) The three components in the x, y, z directions are:
wherein: t is(n)(xi,yj,zk) Is a total field of the temperature of the current space domain corresponding to the nth iteration, and the background field temperature T of the space domain during the first iteration0(xi,yj,zk) As a current total field of spatial domain temperature, λ, corresponding to the first iterationa(xi,yj,zk) Is set as (x) at the coordinate positioni,yj,zk) The abnormal thermal conductivity of the subdivision nodes.
7. The three-dimensional geothermal field numerical simulation method based on a thermal conductivity anisotropic medium of claim 6, wherein the spatial-wavenumber mixed domain coordinate position is (k)x,ky,zk) The spatial wave number mixed domain abnormal heat flux density on the subdivision node
The three components in the x, y, z directions are:
wherein: k is a radical ofx，kyDiscrete offset wavenumbers in the x, y directions, e represents a natural constant,
denotes the unit of imaginary number, (k)x,ky,zk) Is set as (x) at the coordinate positioni,yj,zk) And (4) corresponding coordinate positions of the subdivision nodes in the space wave number mixed domain.
8. The three-dimensional geothermal field numerical simulation method based on thermal conductivity anisotropic medium of claim 7, wherein the spatial wavenumber mixed domain abnormal field temperature
The following governing equation is satisfied:
in the formula,
kx，kydiscrete offset wavenumbers in the x, y directions;
indicates a coordinate position of (x)i,yj,zk) The abnormal heat generation rate of the space wave number mixed domain of the subdivision node is (x) from the coordinate positioni,yj,zk) Abnormal heat generation rate Q of the divided nodea(xi,yj,zk) Performing two-dimensional Fourier transform to obtain the Fourier transform; i is an imaginary unit;
solving the control equation by using a finite element method based on quadratic function interpolation to obtain the abnormal field temperature of the spatial wave number mixed domain
9. The method for numerical simulation of a three-dimensional geothermal field based on a thermally conductive anisotropic medium according to claim 7 or 8, wherein k isx，kyThe calculation method of (2) is as follows:
kx＝p·Δkx
ky＝q·Δky
wherein,
delta x and delta y are unit interval lengths between adjacent subdivision nodes when the three-dimensional prism model carries out grid subdivision along the x direction and the y direction respectively, and Nx、NyRespectively performing mesh subdivision on the three-dimensional prism model and then dividing the node number in the x and y directions;
when N is presentx、NyIn the case of an even number, the number of the first,
when N is presentx、NyIn the case of an odd number of the groups,
10. the method for numerical simulation of a three-dimensional geothermal field based on a thermal conductivity anisotropic medium according to claim 1,2, 3, 4, 6, 7 or 8, wherein the iterative convergence conditions are set as follows:
abs(T(n+1)(xi,yj,zk)-T(n)(xi,yj,zk))/abs(T(n+1)(xi,yj,zk))＜ε0
wherein: t is(n+1)(xi,yj,zk) Represents the new total temperature field of the space domain, epsilon, obtained by the nth iteration calculation0Abs () represents the absolute value for a set numerical precision;
and if the iteration convergence condition is met, outputting a new space domain temperature total field obtained by the nth iteration calculation, and if the iteration convergence condition is not met, taking the new space domain temperature total field obtained by the nth iteration calculation as the current space domain temperature total field corresponding to the (n + 1) th iteration.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111423020.6A CN114036806B (en) | 2021-11-26 | 2021-11-26 | Three-dimensional ground temperature field numerical simulation method based on heat conductivity anisotropic medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111423020.6A CN114036806B (en) | 2021-11-26 | 2021-11-26 | Three-dimensional ground temperature field numerical simulation method based on heat conductivity anisotropic medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114036806A true CN114036806A (en) | 2022-02-11 |
| CN114036806B CN114036806B (en) | 2024-07-09 |
Family
ID=80145735
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111423020.6A Active CN114036806B (en) | 2021-11-26 | 2021-11-26 | Three-dimensional ground temperature field numerical simulation method based on heat conductivity anisotropic medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114036806B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116187107A (en) * | 2023-04-27 | 2023-05-30 | 中南大学 | Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130132050A1 (en) * | 2011-11-21 | 2013-05-23 | Schlumberger Technology Corporation | Dynamic Prediction of Downhole Temperature Distributions |
| JP2017049971A (en) * | 2015-09-03 | 2017-03-09 | 住友重機械工業株式会社 | Simulation method, simulation device and simulation program |
| CN113656750A (en) * | 2021-10-20 | 2021-11-16 | 中南大学 | Magnetic induction intensity calculation method of strong magnetic medium based on space wave number domain |
-
2021
- 2021-11-26 CN CN202111423020.6A patent/CN114036806B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130132050A1 (en) * | 2011-11-21 | 2013-05-23 | Schlumberger Technology Corporation | Dynamic Prediction of Downhole Temperature Distributions |
| JP2017049971A (en) * | 2015-09-03 | 2017-03-09 | 住友重機械工業株式会社 | Simulation method, simulation device and simulation program |
| CN113656750A (en) * | 2021-10-20 | 2021-11-16 | 中南大学 | Magnetic induction intensity calculation method of strong magnetic medium based on space wave number domain |
Non-Patent Citations (2)
| Title |
|---|
| 余正良;赵敏;: "巷道导水构造体围岩温度场-渗流场影响因素分析", 现代矿业, no. 01, 25 January 2018 (2018-01-25), pages 100 - 106 * |
| 张建平;胡慧瑶;王树森;龚曙光;刘庭显;: "正交各向异性结构的三维无网格法稳态传热模型及应用", 材料导报, no. 08, 18 April 2020 (2020-04-18), pages 8036 - 8041 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116187107A (en) * | 2023-04-27 | 2023-05-30 | 中南大学 | Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium |
| CN116187107B (en) * | 2023-04-27 | 2023-08-11 | 中南大学 | Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114036806B (en) | 2024-07-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108345741B (en) | Two-dimensional thermal deformation and thermal stress analysis method based on non-grid RKPM anisotropic material | |
| CN112287534A (en) | NUFFT-based two-dimensional magnetic anomaly fast forward modeling method and device | |
| CN114036805A (en) | Forward modeling method, device, equipment and medium for three-dimensional steady-state heat conduction geothermal field | |
| CN116842813B (en) | Three-dimensional geoelectromagnetic forward modeling method | |
| CN114065511B (en) | Method, device, equipment and medium for simulating two-dimensional forward modeling of magnetotelluric under undulating terrain | |
| Zhang et al. | Coupled models and parallel simulations for three-dimensional full-Stokes ice sheet modeling | |
| CN116522713A (en) | Aviation electromagnetic data partition three-dimensional inversion method based on partition coordinate descent algorithm | |
| CN113420487A (en) | Three-dimensional gravity gradient tensor numerical simulation method, device, equipment and medium | |
| CN114036806B (en) | Three-dimensional ground temperature field numerical simulation method based on heat conductivity anisotropic medium | |
| CN113255230A (en) | Gravity model forward modeling method and system based on MQ radial basis function | |
| CN107247686B (en) | FETD simulation method based on parallel algorithm | |
| Li et al. | Verified error bounds for isolated singular solutions of polynomial systems: case of breadth one | |
| CN111241728A (en) | Intermittent Galerkin finite element numerical solution method of Euler equation | |
| Yang et al. | An upwind CESE scheme for 2D and 3D MHD numerical simulation in general curvilinear coordinates | |
| CN113076678A (en) | Frequency domain two-degree body weight abnormity rapid numerical simulation method and device | |
| CN111339688B (en) | Method for solving rocket simulation model time domain equation based on big data parallel algorithm | |
| CN103530451A (en) | Multi-grid Chebyshev parallel spectral element method with complex medium elastic wave propagation and simulation | |
| CN114970289B (en) | Three-dimensional magnetotelluric anisotropy forward modeling numerical simulation method, equipment and medium | |
| Tanbay et al. | Numerical solution of the multigroup neutron diffusion equation by the meshless RBF collocation method | |
| CN113642189B (en) | Gravity gradient tensor rapid numerical simulation method and device based on product decomposition | |
| CN113779818B (en) | Three-dimensional geologic body electromagnetic field numerical simulation method, device, equipment and medium thereof | |
| CN114462264B (en) | Cable joint state evaluation method and device based on multi-physical field coupling | |
| CN113673163A (en) | Three-dimensional magnetic anisotropy constant fast forward modeling method and device and computer equipment | |
| Ostermann | Three-dimensional modeling of heat transport in deep hydrothermal reservoirs | |
| Ang et al. | Nonlinear heat equation for nonhomogeneous anisotropic materials: A dual‐reciprocity boundary element solution |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |