CN113515863A - Method for calculating heat quantity of middle-deep sleeve type heat exchanger based on numerical inversion - Google Patents
Method for calculating heat quantity of middle-deep sleeve type heat exchanger based on numerical inversion Download PDFInfo
- Publication number
- CN113515863A CN113515863A CN202110793294.8A CN202110793294A CN113515863A CN 113515863 A CN113515863 A CN 113515863A CN 202110793294 A CN202110793294 A CN 202110793294A CN 113515863 A CN113515863 A CN 113515863A
- Authority
- CN
- China
- Prior art keywords
- heat exchanger
- temperature
- heat
- formula
- pipe
- 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
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
- 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)
- Heat-Exchange Devices With Radiators And Conduit Assemblies (AREA)
Abstract
The invention provides a method for calculating the heat quantity of a middle-deep sleeve type heat exchanger based on numerical inversion, which comprises the following steps: step one, establishing a middle-deep sleeve type heat exchanger model, wherein the middle-deep sleeve type heat exchanger model is divided into two parts by taking a borehole wall as a boundary, one part is a heat exchanger and a well cementation cement sheath in a borehole, and the other part is a stratum outside the borehole; step two, solving a formation heat conduction equation: step three, solving the temperature of the circulating water in the double-pipe heat exchanger: step four, solving the steps from step two to step three to obtain the water temperature in the double-pipe heat exchangerAndthe solution in the frequency domain. The invention is in the given placeUnder the conditions of the layer initial temperature distribution function, the heat exchanger inlet water temperature and the related physical parameters of the heat exchanger and the stratum, the distribution rule of the outlet water temperature of the double-pipe heat exchanger and the water temperature in the heat exchanger along the depth direction at any moment can be obtained, and then the change of the heat taken by the double-pipe heat exchanger along with the time can be obtained.
Description
Technical Field
The invention belongs to the technical field of development and utilization of deep geothermal energy, relates to heat extraction, and particularly relates to a method for calculating heat extraction of a deep sleeve type heat exchanger based on numerical inversion.
Background
A mid-deep ground source heat pump system is a technology for indirectly utilizing geothermal energy, which obtains heat from the ground through a double pipe heat exchanger installed in a borehole. The double pipe heat exchanger structure is shown in figure 1. The double-pipe heat exchanger penetrates into a high-temperature rock stratum below 2000 m underground through a drill hole, cold water flows in from an annular area between the inner sleeve pipe and the outer sleeve pipe, heat is absorbed from the surrounding stratum when the cold water flows downwards, and the heated water is pumped out of the ground through the inner pipe for heating. The initial temperature of the formation rock mass around the geothermal well can reach more than 100K at most, the temperature of the formation rock mass around the drill hole is gradually reduced along with the heat extraction process of the double-pipe heat exchanger, and meanwhile, the heat extraction quantity of the double-pipe heat exchanger can also be continuously reduced. The method has the advantages that the water temperature distribution in the heat exchanger and the outlet water temperature are accurately calculated, and the method is particularly important for evaluating the heat extraction capacity of the geothermal well and optimizing the design of the geothermal well.
The heat transfer process of the intermediate-deep double-pipe heat exchanger is very complex, the involved time span is long, and the space area is large. Compared with a shallow layer buried pipe heat exchanger, the vertical variation range of the formation temperature is large, and the heat exchange process of the shallow layer buried pipe heat exchanger relates to a very complex mathematical physics problem. At present, researches on a method for calculating heat of a shallow geothermal well heat exchanger are numerous, because the depth of the shallow geothermal well is generally not more than 200m, the vertical variation of the formation temperature is very small, the formation temperature is assumed to be a constant value when the analytic solution of the shallow geothermal well heat exchanger is solved, and the vertical variation of the formation temperature is not considered. Therefore, the method for calculating the heat quantity of the shallow geothermal well heat exchanger cannot be applied to the calculation of the heat quantity of the intermediate and deep geothermal well heat exchanger.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide a method for calculating the heat extraction amount of a middle-deep sleeve type heat exchanger based on numerical inversion, and solve the technical problem that the method for calculating the heat extraction amount of a shallow geothermal well heat exchanger in the prior art cannot be applied to the calculation of the heat extraction amount of the middle-deep sleeve type heat exchanger.
In order to solve the technical problems, the invention adopts the following technical scheme:
a method for calculating heat quantity of a middle-deep sleeve type heat exchanger based on numerical inversion comprises the following steps:
step one, establishing a middle-deep sleeve type heat exchanger model, wherein the middle-deep sleeve type heat exchanger model is divided into two parts by taking a borehole wall as a boundary, one part is a heat exchanger and a well cementation cement sheath in a borehole, and the other part is a stratum outside the borehole;
on the basis of the change of the formation temperature along with the depth, the heat transfer processes of the two parts are respectively calculated, and the two parts are coupled through the temperature of the wall of the drill hole;
step two, solving a formation heat conduction equation:
assuming that the heat conduction of the formation in the depth direction is zero, at any depth z along the axial direction of the double-pipe heat exchanger, the formation heat conduction is a one-dimensional process under a cylindrical coordinate, and can be described by the following partial differential equation:
initial conditions: t when T is 0e=az+Tair;
Solving to obtain:
in the formula:
q is the radial heat loss of the stratum at a certain depth and is expressed by W;
ρeis the density of the stratum rock mass with the unit of kg/m3;
ceThe specific heat capacity of the stratum rock mass is expressed as J/(kg.K);
Tethe temperature of the stratum rock mass at the depth z is represented by K;
Tezis the formation initial temperature in K;
t is time in units of s;
t0is a time factor of the formation temperature, dimensionless;
Tairis the surface temperature in K;
λethe heat conductivity coefficient of the stratum rock mass is W/(m.K);
r is the distance from the axis of the borehole in m;
rbis the borehole radius in m;
dQ is the radial heat loss of the formation at the borehole wall in units of W;
Tbis the borehole wall temperature in K;
z is the distance from the earth's surface in m;
a is the earth temperature gradient with the unit of K/m;
step three, solving the temperature of the circulating water in the double-pipe heat exchanger:
when circulating water flows in from the outer pipe and flows out from the inner pipe, the heat balance equation along the flow direction is expressed by a formula VII and a formula VIII;
inner tube heat balance equation:
outer tube heat balance equation:
wherein:
the conditions for the solution were as follows:
T2(0,t)=Tinjformula XI;
T1(H,t)=T2(H, t) formula XII;
in the formula:
A1is the water passing cross-sectional area of the inner pipe, and the unit is m2;
A2Is the cross-sectional area of water passing between the inner and outer tubes, and has unit of m2;
ρfIs the density of water in kg/m3;
c is the specific heat capacity of water and the unit is J/(kg. K);
T1is the temperature of water in the inner tube in unitsIs K;
T2is the temperature of the water in the outer tube, in K;
Tbis the borehole wall temperature in K;
w is the flow of water in the sleeve, and the unit is kg/s;
r1is the inner radius of the inner tube, and the unit is m;
r2is the outer radius of the inner tube, and is expressed in m;
r3is the inner radius of the outer tube, and is expressed in m;
r4is the outer radius of the outer tube, in m;
rbis the borehole radius in m;
h is the geothermal well depth, and the unit is m;
h1is the convective heat transfer coefficient of the inner wall of the inner tube and has the unit of W/(m)2·K);
h2Is the convective heat transfer coefficient of the outer wall of the inner pipe and has the unit of W/(m)2·K);
h3Is the convective heat transfer coefficient of the outer tube wall, and has the unit of W/(m)2·K);
kiThe heat conductivity coefficient of the inner tube wall is W/(m.K);
kothe thermal conductivity coefficient of the outer tube wall is W/(m.K);
kgthe unit is W/(m.K) for the heat conductivity coefficient of the cementing material;
Tinjwater temperature at the inlet, in units of K;
setting the heat quantity flowing out of the borehole wall to be equal to the heat quantity obtained by the outer pipe from the stratum, the heat transfer coupling condition between the sleeve and the stratum is as follows:
from the formula XIII:
order:
in the formula:
z is a geothermal well depth factor and is dimensionless;
a is the total heat transfer coefficient factor of the inner pipe wall, and is dimensionless;
b is a factor of the total heat transfer coefficient between the outer pipe wall and the well wall, and is dimensionless;
φ1is a water temperature factor in the inner pipe, and has no dimension;
φ2is the water temperature factor in the outer tubeNo dimension;
tDthe water temperature time factor in the double-pipe heat exchanger is dimensionless;
M1、M2all are drilling aperture factors and have no dimension;
substituting formula xiv into formula viii, dimensionless the formula vii, formula viii and the conditions of the solution:
wherein:
φ2(0,tD)=1;
φ1(1,tD)=φ2(1,tD);
in the formula:
beta is a factor of the geothermal gradient and is dimensionless;
laplace transformations of the formulae XXIII and XXIV:
in the formula:
s is a complex variable and is dimensionless;
solving the linear differential equations of the formulae XXIII and XXIV to obtain:
wherein:
a1、a2、k1、k2、C1and C2All are intermediate quantities;
a1=M1s+(1+k1)A;
a2=M1s+(1+k2)A;
C1、C2is obtained by dissolving after substituting the formula XXV and the formula XXVI under the definite dissolving condition;
step four, solving the steps from step two to step three to obtain the water temperature in the double-pipe heat exchangerAndthe solution in the frequency domain.
Compared with the prior art, the invention has the following technical effects:
according to the invention, under the condition of giving a stratum initial temperature distribution function, the heat exchanger inlet water temperature and the related physical parameters of the heat exchanger and the stratum, the distribution rule of the outlet water temperature of the double-pipe heat exchanger and the water temperature in the heat exchanger along the depth direction at any moment can be obtained, and then the change of the heat extraction quantity of the double-pipe heat exchanger along with the time can be obtained. The calculation result can be used as the basis for evaluating the heat extraction capacity of the medium-deep geothermal well and optimizing the design scheme of the geothermal well.
Drawings
Fig. 1(a) is a schematic view of a vertical sectional structure of a double pipe heat exchanger.
Fig. 1(b) is a schematic view of a horizontal sectional structure of the double pipe heat exchanger.
Fig. 2 is a graph of outlet water temperature versus time.
The meaning of the individual reference symbols in the figures is: 1-stratum, 2-drilling wall, 3-cementing cement sheath, 4-outer well, 5-inner well, 6-heat-taking circulation water inlet, and 7-heat-taking circulation water outlet.
The present invention will be explained in further detail with reference to examples.
Detailed Description
For a middle-deep geothermal well, the vertical variation of the formation temperature is large, and the assumption of constant geothermal temperature inevitably causes that the calculated value of the heat exchange medium temperature in the double-pipe heat exchanger deviates from the true value in the flowing direction thereof to be overlarge. Therefore, the influence of vertical variation of the ground temperature is necessarily considered when the heat extraction process of the middle-deep sleeve type heat exchanger is researched.
The invention aims to provide a method for calculating heat extraction quantity of a casing pipe type heat exchanger of a middle-deep geothermal well, which is characterized by comprising the following steps: vertical variation of the formation temperature is considered in the solving process, but vertical heat conduction of the formation is not considered; simplifying the convective heat transfer process in the double-pipe heat exchanger into a one-dimensional linear problem and processing according to an unsteady state; in the solving process, firstly, a heat conduction equation of a stratum is solved to obtain the temperature of the wall of the drill hole, then, the solution of a differential equation set on a frequency domain is obtained by using laplace transformation solving, and finally, the water temperature distribution in the double-pipe heat exchanger and the water temperature at an outlet are obtained by using a numerical inversion method to obtain the heat extraction quantity of the double-pipe heat exchanger. The invention provides a technical basis for simplifying the calculation method of the geothermal well heat extraction quantity, estimating the geothermal well heat extraction capacity and optimizing the geothermal well design.
It is to be understood that all parts, devices and apparatus of the present invention, unless otherwise specified, are intended to be included within the scope of the present invention as defined in the appended claims.
The following embodiments of the present invention are provided, and it should be noted that the present invention is not limited to the following embodiments, and all equivalent changes based on the technical solutions of the present invention are within the protection scope of the present invention.
Example (b):
the embodiment provides a method for calculating the heat extraction quantity of a middle-deep double-pipe heat exchanger based on numerical inversion, which comprises the following steps:
step one, establishing a middle-deep sleeve type heat exchanger model, wherein the middle-deep sleeve type heat exchanger model is divided into two parts by taking a borehole wall as a boundary, one part is a heat exchanger and a well cementation cement sheath in a borehole, and the other part is a stratum outside the borehole;
on the basis of the change of the formation temperature along with the depth, the heat transfer processes of the two parts are respectively calculated, and the two parts are coupled through the temperature of the wall of the drill hole;
step two, solving a formation heat conduction equation:
assuming that the heat conduction of the formation in the depth direction is zero, at any depth z along the axial direction of the double-pipe heat exchanger, the formation heat conduction is a one-dimensional process under a cylindrical coordinate, and can be described by the following partial differential equation:
initial conditions: t when T is 0e=az+Tair;
Solving to obtain:
in the formula:
q is the radial heat loss of the stratum at a certain depth and is expressed by W;
ρeis the density of the stratum rock mass with the unit of kg/m3;
ceThe specific heat capacity of the stratum rock mass is expressed as J/(kg.K);
Tethe temperature of the stratum rock mass at the depth z is represented by K;
Tezis the formation initial temperature in K;
t is time in units of s;
t0is a time factor of the formation temperature, dimensionless;
Tairis the surface temperature in K;
λethe heat conductivity coefficient of the stratum rock mass is W/(m.K);
r is the distance from the axis of the borehole in m;
rbis the borehole radius in m;
dQ is the radial heat loss of the formation at the borehole wall in units of W;
Tbis the borehole wall temperature in K;
z is the distance from the earth's surface in m;
a is the earth temperature gradient with the unit of K/m;
step three, solving the temperature of the circulating water in the double-pipe heat exchanger:
when circulating water flows in from the outer pipe and flows out from the inner pipe, the heat balance equation along the flow direction is expressed by a formula VII and a formula VIII;
inner tube heat balance equation:
outer tube heat balance equation:
wherein:
the conditions for the solution were as follows:
T2(0,t)=Tinjformula XI;
T1(H,t)=T2(H, t) formula XII;
in the formula:
A1is the water passing cross-sectional area of the inner pipe, and the unit is m2;
A2Is the cross-sectional area of water passing between the inner and outer tubes, and has unit of m2;
ρfIs the density of water in kg/m3;
c is the specific heat capacity of water and the unit is J/(kg. K);
T1the temperature of water in the inner tube is expressed in K;
T2is the temperature of the water in the outer tube, in K;
Tbis the borehole wall temperature in K;
w is the flow of water in the sleeve, and the unit is kg/s;
r1is the inner radius of the inner tube, and the unit is m;
r2is the outer radius of the inner tube, and is expressed in m;
r3is the inner radius of the outer tube, and is expressed in m;
r4is the outer radius of the outer tube, in m;
rbis the borehole radius in m;
h is the geothermal well depth, and the unit is m;
h1is the convective heat transfer coefficient of the inner wall of the inner tube and has the unit of W/(m)2·K);
h2Is the convective heat transfer coefficient of the outer wall of the inner pipe and has the unit of W/(m)2·K);
h3Is the convective heat transfer coefficient of the outer tube wall, and has the unit of W/(m)2·K);
kiThe heat conductivity coefficient of the inner tube wall is W/(m.K);
kothe thermal conductivity coefficient of the outer tube wall is W/(m.K);
kgthe unit is W/(m.K) for the heat conductivity coefficient of the cementing material;
Tinjwater temperature at the inlet, in units of K;
setting the heat quantity flowing out of the borehole wall to be equal to the heat quantity obtained by the outer pipe from the stratum, the heat transfer coupling condition between the sleeve and the stratum is as follows:
from the formula XIII:
order:
in the formula:
z is a geothermal well depth factor and is dimensionless;
a is the total heat transfer coefficient factor of the inner pipe wall, and is dimensionless;
b is a factor of the total heat transfer coefficient between the outer pipe wall and the well wall, and is dimensionless;
φ1is a water temperature factor in the inner pipe, and has no dimension;
φ2is a water temperature factor in the outer pipe, and has no dimension;
tDthe water temperature time factor in the double-pipe heat exchanger is dimensionless;
M1、M2all are drilling aperture factors and have no dimension;
substituting formula xiv into formula viii, dimensionless the formula vii, formula viii and the conditions of the solution:
wherein:
φ2(0,tD)=1;
φ1(1,tD)=φ2(1,tD);
in the formula:
beta is a factor of the geothermal gradient and is dimensionless;
laplace transformations of the formulae XXIII and XXIV:
in the formula:
s is a complex variable and is dimensionless;
solving the linear differential equations of the formulae XXIII and XXIV to obtain:
wherein:
a1、a2、k1、k2、C1and C2All are intermediate quantities;
a1=M1s+(1+k1)A;
a2=M1s+(1+k2)A;
C1、C2is obtained by dissolving after substituting the formula XXV and the formula XXVI under the definite dissolving condition;
step four, solving the steps from step two to step three to obtain the water temperature in the double-pipe heat exchangerAndthe solution in the frequency domain.
In the invention, a time domain solution phi of the water temperature in the double-pipe heat exchanger is obtained1And phi2Need to be aligned withAndperforming inverse laplace transform becauseAndthe form of the double-pipe heat exchanger is complex, laplace inverse change cannot be directly carried out, and a numerical inversion method is required to obtain a solution phi of the water temperature in the double-pipe heat exchanger on a time domain1And phi2And then T is obtained1(T, z) and T2(t,z)。
Application example:
in this application example, the change of the outlet water temperature of the double pipe heat exchanger with time is calculated by using the analytic solution of the method for calculating the heat quantity of the deep double pipe heat exchanger based on numerical inversion in the above embodiment 1, and the parameters required for calculation are shown in table 1. In addition, in the application example, the surface temperature is a constant value of 16K, the ground temperature gradient is 0.03K/m, and the water temperature at the inlet of the heat exchanger is 15K.
TABLE 1 analysis and calculation parameter table for double-pipe heat exchanger
Taking the first heating season as an example, the time-dependent change relationship of the outlet water temperature when the formation thermal conductivity is 2.5W/(m.K) is obtained, as shown in FIG. 2. The outlet water temperature is the same as the inlet water temperature at time zero. The outlet water temperature rose rapidly in the initial phase and was at its maximum at 0.2 days, with a maximum temperature of about 60K, at which time the heat removal was 1500 kW. Then, the outlet water temperature gradually decreases again as the heat extraction progresses. After continuously taking heat for 120 days, the outlet water temperature is 27K, and the heat taking amount is 420 kW. Compared with the analytic solution obtained by the steady-state treatment of the convection heat transfer in the double-pipe heat exchanger, the analytic solution has the trend that the water temperature at the outlet at the initial stage of calculation is increased firstly and then decreased, and the analytic solution is more in line with the actual situation.
Claims (1)
1. A method for calculating heat quantity of a middle-deep sleeve type heat exchanger based on numerical inversion is characterized by comprising the following steps:
step one, establishing a middle-deep sleeve type heat exchanger model, wherein the middle-deep sleeve type heat exchanger model is divided into two parts by taking a borehole wall as a boundary, one part is a heat exchanger and a well cementation cement sheath in a borehole, and the other part is a stratum outside the borehole;
on the basis of the change of the formation temperature along with the depth, the heat transfer processes of the two parts are respectively calculated, and the two parts are coupled through the temperature of the wall of the drill hole;
step two, solving a formation heat conduction equation:
assuming that the heat conduction of the formation in the depth direction is zero, at any depth z along the axial direction of the double-pipe heat exchanger, the formation heat conduction is a one-dimensional process under a cylindrical coordinate, and can be described by the following partial differential equation:
initial conditions: t when T is 0e=az+Tair;
Solving to obtain:
in the formula:
q is the radial heat loss of the stratum at a certain depth and is expressed by W;
ρeis the density of the stratum rock mass with the unit of kg/m3;
ceThe specific heat capacity of the stratum rock mass is expressed as J/(kg.K);
Tethe temperature of the stratum rock mass at the depth z is represented by K;
Tezis the formation initial temperature in K;
t is time in units of s;
t0is a time factor of the formation temperature, dimensionless;
Tairis the surface temperature in K;
λethe heat conductivity coefficient of the stratum rock mass is W/(m.K);
r is the distance from the axis of the borehole in m;
rbis the borehole radius in m;
dQ is the radial heat loss of the formation at the borehole wall in units of W;
Tbis the borehole wall temperature in K;
z is the distance from the earth's surface in m;
a is the earth temperature gradient with the unit of K/m;
step three, solving the temperature of the circulating water in the double-pipe heat exchanger:
when circulating water flows in from the outer pipe and flows out from the inner pipe, the heat balance equation along the flow direction is expressed by a formula VII and a formula VIII;
inner tube heat balance equation:
outer tube heat balance equation:
wherein:
the conditions for the solution were as follows:
T2(0,t)=Tinjformula XI;
T1(H,t)=T2(H, t) formula XII;
in the formula:
A1is the water passing cross-sectional area of the inner pipe, and the unit is m2;
A2Is the cross-sectional area of water passing between the inner and outer tubes, and has unit of m2;
ρfIs the density of water in kg/m3;
c is the specific heat capacity of water and the unit is J/(kg. K);
T1the temperature of water in the inner tube is expressed in K;
T2is the temperature of the water in the outer tube, in K;
Tbis the borehole wall temperature in K;
w is the flow of water in the sleeve, and the unit is kg/s;
r1is the inner radius of the inner tube, and the unit is m;
r2is the outer radius of the inner tube, and is expressed in m;
r3is the inner radius of the outer tube, and is expressed in m;
r4is the outer radius of the outer tube, in m;
rbis the borehole radius in m;
h is the geothermal well depth, and the unit is m;
h1is the convective heat transfer coefficient of the inner wall of the inner tube and has the unit of W/(m)2·K);
h2Is the convective heat transfer coefficient of the outer wall of the inner pipe and has the unit of W/(m)2·K);
h3Is the convective heat transfer coefficient of the outer tube wall, and has the unit of W/(m)2·K);
kiThe heat conductivity coefficient of the inner tube wall is W/(m.K);
kothe thermal conductivity coefficient of the outer tube wall is W/(m.K);
kgthe unit is W/(m.K) for the heat conductivity coefficient of the cementing material;
Tinjwater temperature at the inlet, in units of K;
setting the heat quantity flowing out of the borehole wall to be equal to the heat quantity obtained by the outer pipe from the stratum, the heat transfer coupling condition between the sleeve and the stratum is as follows:
from the formula XIII:
order:
in the formula:
z is a geothermal well depth factor and is dimensionless;
a is the total heat transfer coefficient factor of the inner pipe wall, and is dimensionless;
b is a factor of the total heat transfer coefficient between the outer pipe wall and the well wall, and is dimensionless;
φ1is a water temperature factor in the inner pipe, and has no dimension;
φ2is a water temperature factor in the outer pipe, and has no dimension;
tDthe water temperature time factor in the double-pipe heat exchanger is dimensionless;
M1、M2all are drilling aperture factors and have no dimension;
substituting formula xiv into formula viii, dimensionless the formula vii, formula viii and the conditions of the solution:
wherein:
φ2(0,tD)=1;
φ1(1,tD)=φ2(1,tD);
in the formula:
beta is a factor of the geothermal gradient and is dimensionless;
laplace transformations of the formulae XXIII and XXIV:
in the formula:
s is a complex variable and is dimensionless;
solving the linear differential equations of the formulae XXIII and XXIV to obtain:
wherein:
a1、a2、k1、k2、C1and C2All are intermediate quantities;
a1=M1s+(1+k1)A;
a2=M1s+(1+k2)A;
C1、C2is obtained by dissolving after substituting the formula XXV and the formula XXVI under the definite dissolving condition;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110793294.8A CN113515863B (en) | 2021-07-14 | 2021-07-14 | Method for calculating heat quantity of middle-deep sleeve type heat exchanger based on numerical inversion |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110793294.8A CN113515863B (en) | 2021-07-14 | 2021-07-14 | Method for calculating heat quantity of middle-deep sleeve type heat exchanger based on numerical inversion |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113515863A true CN113515863A (en) | 2021-10-19 |
CN113515863B CN113515863B (en) | 2022-11-29 |
Family
ID=78066888
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110793294.8A Active CN113515863B (en) | 2021-07-14 | 2021-07-14 | Method for calculating heat quantity of middle-deep sleeve type heat exchanger based on numerical inversion |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113515863B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115270512A (en) * | 2022-08-24 | 2022-11-01 | 西南石油大学 | Optimal flow design method of closed-cycle geothermal utilization system |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3739749A1 (en) * | 1987-11-24 | 1989-06-08 | Mueller Christa Dr | Selective treatment, at the same time and place, of tumours by means of heat and radioactive emissions by two-zone combination needles |
CN109800527A (en) * | 2019-01-30 | 2019-05-24 | 西安交通大学 | A kind of calculation method for self circulation heat exchange amount in the buried double-tube heat exchanger of mid-deep strata |
CN109946103A (en) * | 2019-04-23 | 2019-06-28 | 山东建筑大学 | A kind of underground heat parameter test system and method based on mid-deep strata ground heat exchanger |
CN111551288A (en) * | 2020-05-13 | 2020-08-18 | 山东建筑大学 | Prediction method for circulating liquid temperature of middle-deep layer U-shaped heat exchange well |
CN112597639A (en) * | 2020-12-11 | 2021-04-02 | 山东中瑞新能源科技有限公司 | Heat transfer dimension reduction analysis method and system for porous ground heat exchanger |
-
2021
- 2021-07-14 CN CN202110793294.8A patent/CN113515863B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3739749A1 (en) * | 1987-11-24 | 1989-06-08 | Mueller Christa Dr | Selective treatment, at the same time and place, of tumours by means of heat and radioactive emissions by two-zone combination needles |
CN109800527A (en) * | 2019-01-30 | 2019-05-24 | 西安交通大学 | A kind of calculation method for self circulation heat exchange amount in the buried double-tube heat exchanger of mid-deep strata |
CN109946103A (en) * | 2019-04-23 | 2019-06-28 | 山东建筑大学 | A kind of underground heat parameter test system and method based on mid-deep strata ground heat exchanger |
CN111551288A (en) * | 2020-05-13 | 2020-08-18 | 山东建筑大学 | Prediction method for circulating liquid temperature of middle-deep layer U-shaped heat exchange well |
CN112597639A (en) * | 2020-12-11 | 2021-04-02 | 山东中瑞新能源科技有限公司 | Heat transfer dimension reduction analysis method and system for porous ground heat exchanger |
Non-Patent Citations (2)
Title |
---|
A. SHTYM; I. ZHURMILOVA: ""Ground heat exchangers of geothermal heat pumps and analysis of their constructive features and types"", 《2017 INTERNATIONAL CONFERENCE ON INDUSTRIAL ENGINEERING, APPLICATIONS AND MANUFACTURING (ICIEAM)》 * |
杜甜甜等: "套管式中深层地埋管换热器传热建模及取热分析", 《可再生能源》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115270512A (en) * | 2022-08-24 | 2022-11-01 | 西南石油大学 | Optimal flow design method of closed-cycle geothermal utilization system |
CN115270512B (en) * | 2022-08-24 | 2024-04-16 | 西南石油大学 | Optimum flow design method for closed circulation geothermal utilization system |
Also Published As
Publication number | Publication date |
---|---|
CN113515863B (en) | 2022-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104806230B (en) | Method for calculating wellbore temperature field and pressure field distribution of supercritical steam injection well | |
CN103775058B (en) | Method for determining heat loss of shaft | |
CN110826225B (en) | Vertical single-U-shaped buried pipe heat exchanger temperature field analysis method | |
CN106770439A (en) | Rock-soil layer is layered Determination of conductive coefficients method | |
Wang et al. | A new analytical solution for horizontal geothermal heat exchangers with vertical spiral coils | |
CN111400893A (en) | Method for analyzing fluid temperature field of sleeve type buried pipe heat exchanger | |
CN105370255B (en) | Method and device for determining temperature distribution of gas injection and electric ignition shaft of in-situ combustion cage | |
CN113515863B (en) | Method for calculating heat quantity of middle-deep sleeve type heat exchanger based on numerical inversion | |
CN111539130A (en) | Design and calculation method for drilling depth of geothermal heat exchanger of middle-deep buried pipe | |
Jiashu et al. | The influence of thermal boundary conditions of wellbore on the heat extraction performance of deep borehole heat exchangers | |
Kurevija et al. | Steady-state heat rejection rates for a coaxial borehole heat exchanger during passive and active cooling determined with the novel step thermal response test method | |
CN105422084A (en) | Method and device for acquiring shaft temperature field of hot water circulation heating viscosity reduction process | |
CN107145705A (en) | A kind of method and device for obtaining circulating temperature | |
CN111625967B (en) | Simplified calculation method for heat exchange performance of U-shaped medium-deep buried pipe heat exchanger | |
Bao et al. | An analytical heat transfer model for the mid-deep U-shaped borehole heat exchanger considering groundwater seepage | |
CN206235584U (en) | A kind of experimental system for verifying energy stake Calculation of Heat Transfer model under the conditions of seepage action of ground water | |
CN109033012B (en) | Method and device for determining temperature field of hollow sucker rod hot water injection circulation shaft | |
Xu et al. | A comprehensive investigation on U-tube ground heat exchanger performance considering insulation layer in deep geothermal system | |
CN106934106B (en) | Method and device for acquiring well cementation circulating temperature | |
Liao et al. | A novel 2-D equivalent numerical model of helix energy pile based on heat transfer characteristics of internal heat convection | |
CN113505454B (en) | Method for calculating heat quantity of middle-deep geothermal well casing type heat exchanger | |
Benim et al. | A computational investigation of the thermohydraulics of an EGS project | |
CN113468743B (en) | Medium-deep buried pipe fluid temperature field analysis method considering groundwater seepage | |
CN116628919A (en) | Operation design method for efficient sustainable heat extraction of medium-deep coaxial cased well | |
CN115270512B (en) | Optimum flow design method for closed circulation geothermal utilization system |
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 |