CN116187107A - Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium - Google Patents

Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium Download PDF

Info

Publication number
CN116187107A
CN116187107A CN202310469070.0A CN202310469070A CN116187107A CN 116187107 A CN116187107 A CN 116187107A CN 202310469070 A CN202310469070 A CN 202310469070A CN 116187107 A CN116187107 A CN 116187107A
Authority
CN
China
Prior art keywords
field
temperature
abnormal
space
domain
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
Application number
CN202310469070.0A
Other languages
Chinese (zh)
Other versions
CN116187107B (en
Inventor
戴世坤
贾金荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Central South University
Original Assignee
Central South University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Central South University filed Critical Central South University
Priority to CN202310469070.0A priority Critical patent/CN116187107B/en
Publication of CN116187107A publication Critical patent/CN116187107A/en
Application granted granted Critical
Publication of CN116187107B publication Critical patent/CN116187107B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/10—Geothermal energy

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)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium, which are used for carrying out space discretization and time discretization on a target area, and iteratively calculating the total field temperature of a space domain by combining a background field control equation and an abnormal field control equation in space; the space-wave number domain abnormal field temperature of the next time node can be obtained by combining the ground temperature field recursion formula of the explicit differential format according to the space-wave number domain abnormal field temperature of the current time node, and then the space domain total field temperature of the next time node is calculated in an iterative mode; and so on to get the total field temperature of the spatial domain for all time nodes. The three-dimensional ground temperature field dynamic numerical simulation is realized, the dynamic change of the ground temperature field can be truly reflected, the background field control equation, the abnormal field control equation and the ground temperature field recurrence technology are combined, the calculation and storage cost is reduced while the fine simulation of the temperature field dynamic numerical value is ensured, and the calculation efficiency is improved.

Description

Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium
Technical Field
The invention relates to the technical field of ground temperature field numerical simulation, in particular to a three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium.
Background
Geothermal is a renewable energy source with abundant reserves, wide distribution and green and clean, and plays an important role in the sustainable development of society. The research of a fine and efficient three-dimensional ground temperature field numerical simulation method is of great significance in grasping the temperature distribution in the earth, evaluating the geothermal resource quantity and developing geothermal energy.
At present, most of numerical simulation methods for three-dimensional ground temperature fields are steady-state numerical simulation, for example, chinese patent application CN202111423020.6 discloses a three-dimensional ground temperature field numerical simulation method based on a heat conductivity anisotropic medium, and Chinese patent application CN202111423009.X discloses a forward modeling method, device, equipment and medium of a three-dimensional steady-state heat conduction ground temperature field, which are steady-state numerical simulation methods. However, the ground temperature field is dynamically changing in nature, i.e., transient, and dynamic simulation of the ground temperature field is of great practical significance to geothermal research. At present, although researches on transient ground temperature fields are carried out, the problems of large calculated amount and high storage requirement in the dynamic simulation of the three-dimensional ground temperature field numerical values under large-scale complex geological conditions generally exist in the technology, and the cost is high. Aiming at the problem, the invention provides a fine and efficient three-dimensional ground temperature field dynamic numerical simulation technology, which provides important technical support for fine inversion of large-scale ground temperature data.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide a three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium, which can realize fine and efficient simulation of the three-dimensional ground temperature field dynamic numerical value.
In a first aspect, a method for dynamically and numerically simulating a three-dimensional geothermal field is provided, including:
s1: selecting a target area containing an abnormal body, and constructing a target model;
s2: performing spatial discretization and time discretization on the target model to obtain a series of sampling nodes and a series of time nodes;
s3: giving a thermophysical parameter value on each sampling node, wherein the thermophysical parameter value comprises a background thermophysical parameter value at a sampling node without an abnormal body and an abnormal background thermophysical parameter value at a sampling node with an abnormal body;
s4: loading initial conditions and boundary conditions, and taking the background field temperature of the space domain as the total field temperature of the initial space domain based on a background field control equation and a background thermophysical parameter value;
s5: performing horizontal two-dimensional Fourier transform on the abnormal field control equation to obtain a space-wave number domain abnormal field control equation;
s6: according to the total field temperature of the current round space domain and the abnormal background thermophysical parameter value, calculating to obtain the abnormal field temperature of the space-wave number domain by combining a space-wave number domain abnormal field control equation;
s7: performing two-dimensional Fourier inverse transformation on the abnormal field temperature in the space-wave number domain to obtain the abnormal field temperature in the space domain;
s8: obtaining a new spatial domain total field temperature based on the spatial domain background field temperature and the spatial domain abnormal field temperature;
s9: judging whether the new space domain total field temperature meets the iteration convergence condition, if so, outputting the new space domain total field temperature as the ground temperature of the current time node; otherwise, taking the new space domain total field temperature as the current round space domain total field temperature in the next round of iteration, and returning to the step S6;
s10: combining a ground temperature field recursion formula in an explicit differential format, obtaining space-wave number domain abnormal field temperature of the next time node from space-wave number domain abnormal field temperature in the last iteration of the current time node, performing two-dimensional Fourier inverse transformation on the space-wave number domain abnormal field temperature, combining space domain background field temperature to obtain space domain total field temperature as the current round space domain total field temperature of the next time node, and returning to the step S6; and the like until the ground temperature field temperature of all time nodes is obtained.
According to the first aspect, in a possible implementation manner, in the step S2, when the spatial discretization is performed, the target area is spatially subjected to three-dimensional mesh subdivision, and the x direction and the y direction are uniformly subdivided, and the z direction is uniformly or non-uniformly subdivided.
In a possible implementation manner according to the first aspect, the thermophysical parameter values include thermal conductivity, thermal capacity, heat generation rate.
In one possible implementation form according to the first aspect, the background field control equation is expressed as follows:
Figure SMS_1
in the method, in the process of the invention,
Figure SMS_2
for background field temperature, +.>
Figure SMS_3
For background thermal conductivity, +.>
Figure SMS_4
For background heat generation rate, < >>
Figure SMS_5
Is background heat capacity, t is time.
In one possible implementation manner according to the first aspect, the step S5 includes:
the abnormal field control equation is expressed as follows:
Figure SMS_6
in the method, in the process of the invention,
Figure SMS_7
for abnormal field temperature, +.>
Figure SMS_8
For background thermal conductivity, +.>
Figure SMS_9
For background heat capacity, +.>
Figure SMS_10
For an abnormal thermal conductivity,
Figure SMS_11
for abnormal heat generation rate, < >>
Figure SMS_12
Is the abnormal heat capacity, T is the total field temperature, < >>
Figure SMS_13
Is the heat capacity of the fluid, v is the flow rate of the fluid, and t is the time;
performing horizontal two-dimensional Fourier transform on the abnormal field control equation to obtain a space-wave number domain abnormal field control equation, wherein the space-wave number domain abnormal field control equation is expressed as follows:
Figure SMS_14
in the method, in the process of the invention,
Figure SMS_15
abnormal field temperature for the space-wavenumber domain, +.>
Figure SMS_16
Abnormal heat generation rate for space-wave number domain, < >>
Figure SMS_17
Respectively->
Figure SMS_18
Wave number in direction, +.>
Figure SMS_19
Is a two-dimensional Fourier transform symbol, i is an imaginary unit, v x 、v y 、v z The flow rates of the fluid in the x, y, and z directions, respectively.
In one possible implementation manner, according to the first aspect, the iterative convergence condition is that an error is smaller than a preset value; the expression of the error is as follows:
Figure SMS_20
wherein e is an error,
Figure SMS_21
the number of sampling nodes in the x, y and z directions is +.>
Figure SMS_22
、/>
Figure SMS_23
The spatial domain total field temperature for the nth and n+1th iterations, respectively.
According to the first aspect, in a possible implementation manner, in the step S10, a ground temperature field recurrence formula is expressed as follows:
Figure SMS_24
in the method, in the process of the invention,
Figure SMS_25
for the last time node->
Figure SMS_26
Is a spatially-wavenumber domain abnormal field temperature, +.>
Figure SMS_27
For the next time node
Figure SMS_28
K is a coefficient matrix related to the abnormal field of the space-wave number domain, G is a coefficient matrix related to the time derivative term of the abnormal field of the space-wave number domain, and P is a source term.
According to the first aspect, in a possible implementation manner, in the step S10, a catch-up method is used to solve the abnormal field temperature in the space-wave number domain of the next time node.
In a second aspect, there is provided an electronic device comprising:
a memory storing a computer program;
and the processor is used for loading and executing the computer program to realize the steps of the three-dimensional ground temperature field dynamic numerical simulation method.
In a third aspect, a computer readable storage medium is provided, storing a computer program which, when executed by a processor, implements the steps of the three-dimensional earth temperature field dynamic numerical simulation method as described above.
The invention provides a three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium, which are used for carrying out space discretization and time discretization on a target area, and iteratively calculating the total field temperature of a space domain by combining a background field control equation and an abnormal field control equation in space; the space-wave number domain abnormal field temperature of the next time node can be obtained by combining the ground temperature field recursion formula of the explicit differential format according to the space-wave number domain abnormal field temperature of the current time node, and then the space domain total field temperature of the next time node is calculated in an iterative mode; in this way, the total field temperature of the spatial domain of all time nodes can be obtained. The invention realizes three-dimensional ground temperature field dynamic numerical simulation, is ground temperature field transient simulation, can truly reflect the dynamic change of a ground temperature field, combines a background field control equation, an abnormal field control equation and a ground temperature field recurrence technology, reduces calculation and storage costs while guaranteeing the fine simulation of the temperature field dynamic numerical value, improves calculation efficiency, and provides important technical support for fine inversion of large-scale ground temperature data.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for simulating dynamic numerical values of a three-dimensional ground temperature field provided by an embodiment of the invention;
FIG. 2 is a cross-sectional view of a target model according to an embodiment of the present invention, wherein (a) is a cross-sectional view of the target model XOY, and (b) is a cross-sectional view of the target model XOZ;
FIG. 3 is a schematic diagram of the calculation results of the method and COMSOL Multiphysics software at a sampling node according to the embodiment of the present invention;
FIG. 4 is a graph showing the relative error between the calculation results of the method and COMSOL Multiphysics software according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be described in detail below. It will be apparent that the described embodiments are only some, but not all, embodiments of the invention. All other embodiments, based on the examples herein, which are within the scope of the invention as defined by the claims, will be within the scope of the invention as defined by the claims.
The invention simulates the dynamic value of the ground temperature field of each time node based on the dynamic simulation equation of the ground temperature field, and the dynamic simulation equation of the ground temperature field is expressed as follows:
Figure SMS_29
(1)
in the method, in the process of the invention,
Figure SMS_30
for heat conductivity, < >>
Figure SMS_31
Is the total field temperature, +.>
Figure SMS_32
Is the heat capacity of the fluid, v is the flow rate of the fluid, < >>
Figure SMS_33
The heat-generating rate of the heat-generating material,
Figure SMS_34
is the heat capacity of the medium, t is the time.
Based on the superposition principle, the formula (1) is split into a background field control equation and an abnormal field control equation, wherein the background field is set as a uniform layered medium, and the background field control equation is expressed as follows:
Figure SMS_35
(2)
in the method, in the process of the invention,
Figure SMS_36
for background field temperature, +.>
Figure SMS_37
For background thermal conductivity, +.>
Figure SMS_38
For background heat generation rate, < >>
Figure SMS_39
Is background heat capacity.
The abnormal field control equation is expressed as follows:
Figure SMS_40
(3)
in the method, in the process of the invention,
Figure SMS_41
for abnormal field temperature, +.>
Figure SMS_42
For abnormal thermal conductivity +.>
Figure SMS_43
For abnormal heat generation rate, < >>
Figure SMS_44
Is an abnormal heat capacity.
Based on the background field control equation and the abnormal field control equation, the embodiment of the invention provides a three-dimensional ground temperature field dynamic numerical simulation method, as shown in fig. 1, comprising the following steps:
s1: and selecting a target area containing an abnormal body, and constructing a target model.
S2: and performing spatial discretization and time discretization on the target model to obtain a series of sampling nodes and a series of time nodes.
In this embodiment, when performing spatial discretization, the target region is spatially three-dimensionally meshed, with uniform splitting in the x and y directions and uniform or non-uniform z directionEvenly dividing; the number of sampling nodes in the x direction is
Figure SMS_45
The number of sampling nodes in the y direction is +.>
Figure SMS_46
The number of sampling nodes in the z direction is +.>
Figure SMS_47
. When time discretization is carried out, a series of time nodes are selected in time, and the number of the time nodes is +.>
Figure SMS_48
。
S3: the thermophysical parameter values on each sampling node are given, including a background thermophysical parameter value at a non-outlier sampling node and an outlier background thermophysical parameter value at an outlier sampling node.
In this embodiment, the thermophysical parameter values include thermal conductivity
Figure SMS_49
Heat capacity->
Figure SMS_53
And a heat generation rate A, wherein the abnormal heat physical parameter value is different from the surrounding heat physical parameter value. The heat conductivity, the heat capacity and the heat generation rate at the sampling node without abnormal body are background heat conductivity +.>
Figure SMS_54
Background Heat capacity->
Figure SMS_51
Background Heat production Rate->
Figure SMS_52
The background thermophysical parameter value is unchanged in the horizontal direction and only in the vertical direction. The thermal conductivity, the heat capacity and the heat generation rate at the abnormal body sampling node are abnormal thermal conductivity +.>
Figure SMS_55
Abnormal heat capacity->
Figure SMS_56
Abnormal heat generation rate->
Figure SMS_50
。
S4: loading initial conditions and boundary conditions, and taking the background field temperature of the space domain as the total field temperature of the initial space domain based on a background field control equation and a background thermophysical parameter value.
In particular, the boundary condition refers to a known temperature value or a known heat flux density value or a known heat exchange coefficient at the boundary. The initial condition is given initial time
Figure SMS_57
The initial temperature of each sampling node in the target area can be given as original data, and can be calculated based on boundary conditions and model parameters. On the premise of loading initial conditions and boundary conditions, solving based on the background field control equation of the formula (2) and the background thermophysical parameter values of all sampling nodes to obtain the background field temperature of the space domain.
S5: performing horizontal two-dimensional Fourier transform on the abnormal field control equation to obtain a space-wave number domain abnormal field control equation (namely a one-dimensional space domain abnormal field control equation under a two-dimensional wave number domain), wherein the method is expressed as follows:
Figure SMS_58
in the method, in the process of the invention,
Figure SMS_59
abnormal field temperature for the space-wavenumber domain, +.>
Figure SMS_60
Abnormal heat generation rate for space-wave number domain, < >>
Figure SMS_61
Respectively->
Figure SMS_62
Wave number in direction, +.>
Figure SMS_63
Is a two-dimensional Fourier transform symbol, i is an imaginary unit, v x 、v y 、v z The flow rates of the fluid in the x, y, and z directions, respectively.
S6: and solving by adopting a one-dimensional finite element method based on quadratic interpolation according to the total field temperature of the current round space domain and the abnormal background thermophysical parameter value and by combining a space-wave number domain abnormal field control equation to obtain the space-wave number domain abnormal field temperature.
S7: and performing two-dimensional Fourier inverse transformation on the abnormal field temperature in the space-wave number domain to obtain the abnormal field temperature in the space domain.
S8: and adding the space domain background field temperature and the space domain abnormal field temperature to obtain a new space domain total field temperature.
S9: judging whether the new space domain total field temperature meets the iteration convergence condition, if so, outputting the new space domain total field temperature as the ground temperature of the current time node; otherwise, the new spatial domain total field temperature is taken as the current round of spatial domain total field temperature in the next round of iteration, and the step S6 is returned.
In this embodiment, the iteration convergence condition is that an error is smaller than a preset value; the expression of the error is as follows:
Figure SMS_64
wherein e is an error,
Figure SMS_65
the number of sampling nodes in the x, y and z directions is +.>
Figure SMS_66
、/>
Figure SMS_67
The spatial domain total field temperature for the nth and n+1th iterations, respectively. If the error e between the space domain total field temperature obtained by the current iteration round and the space domain total field temperature obtained by the previous iteration round is smaller than a set value, meeting a convergence condition, ending the iteration and outputting the space domain total field temperature; otherwise, if the convergence condition is not satisfied, repeating the steps S6-S9 to continue the iterative computation until convergence.
S10: combining a ground temperature field recursion formula in an explicit differential format, obtaining space-wave number domain abnormal field temperature of the next time node from space-wave number domain abnormal field temperature in the last iteration of the current time node, performing two-dimensional Fourier inverse transformation on the space-wave number domain abnormal field temperature, combining space domain background field temperature to obtain space domain total field temperature as the current round space domain total field temperature of the next time node, and returning to the step S6; and the like until the ground temperature field temperature of all time nodes is obtained.
Specifically, after the final spatial domain total field temperature of the current time node is calculated by utilizing the steps S6-S9, extracting the spatial-wave number domain abnormal field temperature in the last round of iterative calculation of the spatial domain total field temperature of the current time node, then solving the spatial-wave number domain abnormal field temperature of the next time node by combining a ground temperature field recurrence formula of an explicit differential format, performing two-dimensional Fourier inverse transformation on the spatial-wave number domain abnormal field temperature of the next time node to obtain the spatial domain abnormal field temperature of the next time node, adding the spatial domain abnormal field temperature of the next time node and the spatial domain background field temperature to obtain the spatial domain total field temperature of the current round of spatial domain total field temperature of the next time node, returning to the step S6, and executing the iterative calculation process of S6-S9 to obtain the final spatial domain total field temperature (namely the ground temperature field temperature) of the next time node. Step S10 is then performed and so on until the ground temperature field temperature of all time nodes is obtained.
Wherein, the ground temperature field recurrence formula is expressed as follows:
Figure SMS_68
(5)
in the method, in the process of the invention,
Figure SMS_69
for the last time node->
Figure SMS_70
Is a spatially-wavenumber domain abnormal field temperature, +.>
Figure SMS_71
For the next time node
Figure SMS_72
K is a coefficient matrix related to the abnormal field of the space-wave number domain, G is a coefficient matrix related to a time derivative term of the abnormal field of the space-wave number domain, and P is a source term; in the process of solving the abnormal field control equation of the space-wave number domain by adopting the finite element method to obtain the abnormal field temperature of the space-wave number domain, the value of K, G, P can be obtained by carrying out unit analysis and overall synthesis. Equation (5) can be reduced to a diagonal system of equations +.>
Figure SMS_73
Wherein Q is a five-diagonal matrix, X is the unknown quantity to be solved, and B is the right-hand term. The diagonal equation set is quickly solved by adopting the catch-up method, and the ordinary differential equation is solved by adopting the catch-up method, so that the method has high parallelism, and the calculation time and the memory occupation are greatly reduced.
The embodiment of the invention also provides electronic equipment, which comprises:
a memory storing a computer program;
and the processor is used for loading and executing the computer program to realize the steps of the three-dimensional ground temperature field dynamic numerical simulation method according to the embodiment.
The embodiment of the invention also provides a computer readable storage medium storing a computer program which, when executed by a processor, implements the steps of the three-dimensional earth temperature field dynamic numerical simulation method as described in the above embodiment.
The effect of the solution provided by the invention is examined in connection with a specific example.
The computer tested was configured as an Intel (R) Core (TM) i7-11800H, main frequency 2.30GHz, memory 16GB,64 bit win11 system.
The projection of the target model XOY plane and XOZ cross section is shown in fig. 2, the background is a full space medium, the background thermal conductivity is 2W/(m.deg.c), and the thermal capacity is 5×10 5 J/(m 3 DEG C), a fluid abnormal body is arranged, the abnormal heat conductivity is 1.0W/(m DEG C), the fluid flow direction is vertical upwards, and the flow velocity is 1 multiplied by 10 -9 m/s, abnormal heat capacity of 1X 10 6 J/(m 3 at..degree.C.) no internal heat was generated in the background, and the abnormal heat generation rate was 4X 10 -6 W/m 3 . Calculating the range of x-5000 m, y-5000 m, z-0-5000 m, fluid flow range of x-1000 m, y-1000 m and z-3000-700 m. The upper boundary adopts a first type boundary condition (the temperature value is known), the initial temperature value of the boundary is 10 ℃, the lower boundary adopts a second type boundary condition (the heat flux density value is known), and the heat flux density value is 41.86 mW/m 2 . The spatial grid subdivision is 51 x 51, the number of time nodes is 11. In the single-threaded case, a standard Fourier transform is adopted, and the iterative convergence accuracy (preset value) is 10 -4 . Calculation by the method and COMSOL Multiphysics software of the invention respectively, and FIG. 3 is
Figure SMS_74
The calculation results of the method and COMSOL Multiphysics software of the invention at the measuring points are shown in fig. 4, which is a graph of the relative errors of the two algorithms, and the relative errors of the two calculation results are less than 0.05%, thus illustrating the correctness of the method of the invention. Under the same grid condition, the method takes 23 seconds and occupies 0.14GB of memory; COMSOL Multiphysics takes 426 seconds and occupies 10.48GB of memory, which illustrates the high efficiency of the method of the invention.
It is to be understood that the same or similar parts in the above embodiments may be referred to each other, and that in some embodiments, the same or similar parts in other embodiments may be referred to.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
While embodiments of the present invention have been shown and described above, it will be understood that the above embodiments are illustrative and not to be construed as limiting the invention, and that variations, modifications, alternatives and variations may be made to the above embodiments by one of ordinary skill in the art within the scope of the invention.

Claims (10)

1. The method for simulating the dynamic numerical value of the three-dimensional ground temperature field is characterized by comprising the following steps of:
s1: selecting a target area containing an abnormal body, and constructing a target model;
s2: performing spatial discretization and time discretization on the target model to obtain a series of sampling nodes and a series of time nodes;
s3: giving a thermophysical parameter value on each sampling node, wherein the thermophysical parameter value comprises a background thermophysical parameter value at a sampling node without an abnormal body and an abnormal background thermophysical parameter value at a sampling node with an abnormal body;
s4: loading initial conditions and boundary conditions, and taking the background field temperature of the space domain as the total field temperature of the initial space domain based on a background field control equation and a background thermophysical parameter value;
s5: performing horizontal two-dimensional Fourier transform on the abnormal field control equation to obtain a space-wave number domain abnormal field control equation;
s6: according to the total field temperature of the current round space domain and the abnormal background thermophysical parameter value, calculating to obtain the abnormal field temperature of the space-wave number domain by combining a space-wave number domain abnormal field control equation;
s7: performing two-dimensional Fourier inverse transformation on the abnormal field temperature in the space-wave number domain to obtain the abnormal field temperature in the space domain;
s8: obtaining a new spatial domain total field temperature based on the spatial domain background field temperature and the spatial domain abnormal field temperature;
s9: judging whether the new space domain total field temperature meets the iteration convergence condition, if so, outputting the new space domain total field temperature as the ground temperature of the current time node; otherwise, taking the new space domain total field temperature as the current round space domain total field temperature in the next round of iteration, and returning to the step S6;
s10: combining a ground temperature field recursion formula in an explicit differential format, obtaining space-wave number domain abnormal field temperature of the next time node from space-wave number domain abnormal field temperature in the last iteration of the current time node, performing two-dimensional Fourier inverse transformation on the space-wave number domain abnormal field temperature, combining space domain background field temperature to obtain space domain total field temperature as the current round space domain total field temperature of the next time node, and returning to the step S6; and the like until the ground temperature field temperature of all time nodes is obtained.
2. The method according to claim 1, wherein in the step S2, the target region is spatially three-dimensionally mesh-divided, the x-direction and the y-direction are uniformly divided, and the z-direction is uniformly or non-uniformly divided.
3. The method of claim 1, wherein the thermophysical parameter values include thermal conductivity, thermal capacity, and heat generation rate.
4. The method of claim 1, wherein the background field control equation is expressed as follows:
Figure QLYQS_1
the method comprises the steps of carrying out a first treatment on the surface of the In (1) the->
Figure QLYQS_2
For background field temperature, +.>
Figure QLYQS_3
For background thermal conductivity, +.>
Figure QLYQS_4
For background heat generation rate, < >>
Figure QLYQS_5
Is background heat capacity, t is time.
5. The method of three-dimensional geothermal field dynamic numerical simulation according to claim 1, wherein the step S5 comprises:
the abnormal field control equation is expressed as follows:
Figure QLYQS_7
the method comprises the steps of carrying out a first treatment on the surface of the In (1) the->
Figure QLYQS_9
For abnormal field temperature, +.>
Figure QLYQS_11
For background thermal conductivity, +.>
Figure QLYQS_8
For background heat capacity, +.>
Figure QLYQS_10
For abnormal thermal conductivity +.>
Figure QLYQS_12
For abnormal heat generation rate, < >>
Figure QLYQS_13
Is the abnormal heat capacity, T is the total field temperature, < >>
Figure QLYQS_6
Is the heat capacity of the fluid, v is the flow rate of the fluid, and t is the time;
performing horizontal two-dimensional Fourier transform on the abnormal field control equation to obtain a space-wave number domain abnormal field control equation, wherein the space-wave number domain abnormal field control equation is expressed as follows:
Figure QLYQS_14
；
in the method, in the process of the invention,
Figure QLYQS_15
abnormal field temperature for the space-wavenumber domain, +.>
Figure QLYQS_16
Abnormal heat generation rate for space-wave number domain, < >>
Figure QLYQS_17
Respectively is
Figure QLYQS_18
Wave number in direction, +.>
Figure QLYQS_19
Is a two-dimensional Fourier transform symbol, i is an imaginary unit, v x 、v y 、v z The flow rates of the fluid in the x, y, and z directions, respectively.
6. The method for simulating the dynamic numerical value of the three-dimensional ground temperature field according to claim 1, wherein the iterative convergence condition is that the error is smaller than a preset value; the expression of the error is as follows:
Figure QLYQS_20
the method comprises the steps of carrying out a first treatment on the surface of the Wherein e is error, ">
Figure QLYQS_21
The number of sampling nodes in the x, y and z directions is +.>
Figure QLYQS_22
、/>
Figure QLYQS_23
The spatial domain total field temperature for the nth and n+1th iterations, respectively.
7. The method according to any one of claims 1 to 6, wherein in the step S10, the ground temperature field recurrence formula is as follows:
Figure QLYQS_24
the method comprises the steps of carrying out a first treatment on the surface of the In (1) the->
Figure QLYQS_25
For the last time node->
Figure QLYQS_26
Is a spatially-wavenumber domain abnormal field temperature, +.>
Figure QLYQS_27
For next time node->
Figure QLYQS_28
K is a coefficient matrix related to the abnormal field of the space-wave number domain, G is a coefficient matrix related to the time derivative term of the abnormal field of the space-wave number domain, and P is a source term.
8. The method according to claim 7, wherein in the step S10, the space-wave number domain abnormal field temperature of the next time node is solved by a catch-up method.
9. An electronic device, comprising:
a memory storing a computer program;
a processor for loading and executing the computer program to implement the steps of the three-dimensional earth temperature field dynamic numerical simulation method as claimed in any one of claims 1 to 8.
10. A computer readable storage medium storing a computer program, characterized in that the computer program when executed by a processor implements the steps of the three-dimensional earth temperature field dynamic numerical simulation method according to any one of claims 1 to 8.
CN202310469070.0A 2023-04-27 2023-04-27 Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium Active CN116187107B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310469070.0A CN116187107B (en) 2023-04-27 2023-04-27 Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310469070.0A CN116187107B (en) 2023-04-27 2023-04-27 Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium

Publications (2)

Publication Number Publication Date
CN116187107A true CN116187107A (en) 2023-05-30
CN116187107B CN116187107B (en) 2023-08-11

Family

ID=86452646

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310469070.0A Active CN116187107B (en) 2023-04-27 2023-04-27 Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium

Country Status (1)

Country Link
CN (1) CN116187107B (en)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012027848A1 (en) * 2010-09-03 2012-03-08 Geosoft Inc. Method and system for modeling anomalous density zones in geophysical exploration
WO2015110599A1 (en) * 2014-01-24 2015-07-30 Ledaflow Technologies Da Method for transient quasi three-dimensional simulation of multiphase fluid flow in pipelines
CN108595762A (en) * 2018-03-26 2018-09-28 中国科学院计算技术研究所 A kind of three-dimensional dividing method and system of general circulation model dynamical frame
CN113792445A (en) * 2021-11-15 2021-12-14 中南大学 Three-dimensional magnetotelluric numerical simulation method based on integral equation method
WO2022000132A1 (en) * 2020-06-28 2022-01-06 大连理工大学 Method for designing and optimizing three-dimensional porous heat dissipation structure on basis of three-cycle minimal curved surface
CN114036805A (en) * 2021-11-26 2022-02-11 中南大学 Forward modeling method, device, equipment and medium for three-dimensional steady-state heat conduction geothermal field
CN114036806A (en) * 2021-11-26 2022-02-11 中南大学 Three-dimensional geothermal field numerical simulation method based on thermal conductivity anisotropic medium
CN114707385A (en) * 2022-04-12 2022-07-05 中国地质大学（武汉） Deep stratum heat conductivity coefficient three-dimensional prediction method and device based on Krylov subspace
US20220261517A1 (en) * 2020-10-29 2022-08-18 Jilin University Simulation method for electromagnetic multi-scale diffusion
CN115659579A (en) * 2022-09-05 2023-01-31 中南大学 Three-dimensional gravity field numerical simulation method and system based on 3D AS-FT

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012027848A1 (en) * 2010-09-03 2012-03-08 Geosoft Inc. Method and system for modeling anomalous density zones in geophysical exploration
WO2015110599A1 (en) * 2014-01-24 2015-07-30 Ledaflow Technologies Da Method for transient quasi three-dimensional simulation of multiphase fluid flow in pipelines
CN108595762A (en) * 2018-03-26 2018-09-28 中国科学院计算技术研究所 A kind of three-dimensional dividing method and system of general circulation model dynamical frame
WO2022000132A1 (en) * 2020-06-28 2022-01-06 大连理工大学 Method for designing and optimizing three-dimensional porous heat dissipation structure on basis of three-cycle minimal curved surface
US20220261517A1 (en) * 2020-10-29 2022-08-18 Jilin University Simulation method for electromagnetic multi-scale diffusion
CN113792445A (en) * 2021-11-15 2021-12-14 中南大学 Three-dimensional magnetotelluric numerical simulation method based on integral equation method
CN114036805A (en) * 2021-11-26 2022-02-11 中南大学 Forward modeling method, device, equipment and medium for three-dimensional steady-state heat conduction geothermal field
CN114036806A (en) * 2021-11-26 2022-02-11 中南大学 Three-dimensional geothermal field numerical simulation method based on thermal conductivity anisotropic medium
CN114707385A (en) * 2022-04-12 2022-07-05 中国地质大学（武汉） Deep stratum heat conductivity coefficient three-dimensional prediction method and device based on Krylov subspace
CN115659579A (en) * 2022-09-05 2023-01-31 中南大学 Three-dimensional gravity field numerical simulation method and system based on 3D AS-FT

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李矗东;魏强;李玉中;黄宏虎;: "基于有限差分法的大线能量焊接热模拟试件的温度分布", 工业加热, no. 01 *

Also Published As

Publication number Publication date
CN116187107B (en) 2023-08-11

Similar Documents

Publication Publication Date Title
De Boer et al. Mesh deformation based on radial basis function interpolation
Guo et al. Three-dimensional transient heat conduction analysis by Laplace transformation and multiple reciprocity boundary face method
Roderick et al. Polynomial regression approaches using derivative information for uncertainty quantification
Dehghan et al. Numerical study of three-dimensional Turing patterns using a meshless method based on moving Kriging element free Galerkin (EFG) approach
Wang et al. Solving of two-dimensional unsteady-state heat-transfer inverse problem using finite difference method and model prediction control method
WO2021068901A1 (en) Material microstructure evolution simulation solved on the basis of exponential time difference format
Chang et al. Review of computational schemes in inverse heat conduction problems
Fang et al. Isogeometric boundary element analysis for two-dimensional thermoelasticity with variable temperature
Li et al. Meshless local Petrov Galerkin method for 2D/3D nonlinear convection–diffusion equations based on LS-RBF-PUM
Liang et al. A new alternating iteration strategy based on the proper orthogonal decomposition for solving large-scaled transient nonlinear heat conduction problems
Drakoulas et al. FastSVD-ML–ROM: A reduced-order modeling framework based on machine learning for real-time applications
Zang et al. Isogeometric boundary element for analyzing steady-state heat conduction problems under spatially varying conductivity and internal heat source
Moghaderi et al. Mixed two‐grid finite difference methods for solving one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations
Danko Functional or operator representation of numerical heat and mass transport models
CN116187107B (en) Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium
Liu et al. An efficient data-driven optimization framework for designing graded cellular structures
CN113158538A (en) Soft measurement method for heat flux density of complex-structure boiling surface
Gu et al. A fast inversion approach for the identification of highly transient surface heat flux based on the generative adversarial network
Kadantsev et al. The formulation and implementation of analytic energy gradients for periodic density functional calculations with STO/NAO Bloch basis set
Ding et al. A pressure-based segregated solver for incompressible flow on unstructured grids
CN105677995A (en) Method for numerical solution of fuzzy steady state heat conduction problem based on full grid point collocation theory
CN114036805A (en) Forward modeling method, device, equipment and medium for three-dimensional steady-state heat conduction geothermal field
CN105808508A (en) Random orthogonal expansion method for solving uncertain heat conduction problem
Kramer Model reduction of the coupled Burgers equation in conservation form
CN114036806A (en) Three-dimensional geothermal field numerical simulation method based on thermal conductivity anisotropic medium

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