CN113836700B - Cross-season soil heat storage modeling method and device suitable for multi-energy flow system - Google Patents
Cross-season soil heat storage modeling method and device suitable for multi-energy flow system Download PDFInfo
- Publication number
- CN113836700B CN113836700B CN202111028987.4A CN202111028987A CN113836700B CN 113836700 B CN113836700 B CN 113836700B CN 202111028987 A CN202111028987 A CN 202111028987A CN 113836700 B CN113836700 B CN 113836700B
- Authority
- CN
- China
- Prior art keywords
- soil
- vertical
- shaped pipe
- heat storage
- formula
- 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.)
- Active
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
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- 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
Abstract
The invention relates to the technical field of operation and energy management of a multi-energy flow system, and discloses a cross-season soil heat storage modeling method and device suitable for the multi-energy flow system, wherein the method comprises the following steps: constructing a control equation of a three-dimensional temperature field of soil around the vertical U-shaped pipe; determining initial conditions for solving the control equation; determining boundary conditions between fluid in the vertical U-shaped pipe and the inner wall vertical U-shaped pipe; determining boundary conditions of the outer wall of the vertical U-shaped pipe and soil; determining an energy equation of the soil heat storage system; the heat storage/release rate and the soil heat storage system heat storage capacity limit are determined. According to the cross-season soil heat storage modeling method and device suitable for the multi-energy flow system, the built model comprehensively considers the multi-structure layers inside the cross-season soil heat storage system, the internal physical operation state is described, and the method and device can be effectively used for large-scale optimization of the multi-energy flow system.
Description
Technical Field
The invention relates to the technical field of operation and energy management of a multi-energy flow system, in particular to a cross-season soil heat storage modeling method and device suitable for the multi-energy flow system.
Background
With the pollution, environment and ecological problems caused by the traditional energy heating becoming more serious, renewable energy heating is highly concerned by all countries in the world, and the clean transformation of heating is prominent. However, heating by renewable energy is limited by factors such as weather, regions and seasons, and has the characteristics of strong intermittency, instability and seasonality. The cross-season heat storage stores renewable energy sources in seasons with rich heat sources and low heat demands, and releases the renewable energy sources in seasons with poor heat sources and high heat demands, so that the unmatched characteristics of the renewable energy source heating system in time, space and strength can be effectively solved.
For systems where there is a multi-energy flow of electric heat gas, etc., thermal energy storage is a key technology to address renewable energy uncertainty, using an efficient storage medium to store excess generated heat or cooling fluid for later use in useful applications. Generally, thermal energy storage units can be divided into three broad categories: sensible heat, latent heat and thermochemical storage systems. Heat storage technology is relatively mature at present and has been used for a long time in various heat storage devices using water, rock and soil as common storage media. Such a system is cheap and simple and does not change the material phase by increasing the temperature depending on the specific heat capacity of the storage material.
The cross-season heat storage can transfer heat such as solar energy and industrial waste heat from summer or transition seasons to winter, overcomes the defects of instability and low utilization rate of short-term heat storage, and expands the depth and the breadth of renewable energy utilization. The main research directions are hot water heat storage (large heat capacity, small influence by hydrogeological conditions and large heat storage/release power) and buried pipe heat storage (soil is used as a heat storage body, the problem of recharge is avoided, and underground water quality is not damaged).
As is well known, ground Source Heat Pumps (GSHP) are renewable and environmentally friendly and have been widely used. Suitable areas for GSHP are areas where the subsurface soil temperature is between 10 ℃ and 20 ℃, or where more subsurface heat exchangers may be installed in cold climates. A solar-ground coupled heat pump (SGCHP) can solve the above problems. The basic goal of SGCHP is to achieve a higher heating or cooling coefficient of performance (COP) compared to conventional heat pump systems.
At present, modeling of cross-season soil heat storage in a thermodynamic system at home and abroad mainly focuses on the field of thermodynamic research, and a partial differential equation mathematical model is established through a thermodynamic correlation law. The model only considers the dynamic processes inside the soil heat storage and does not relate to the interaction with other elements in the multi-energy flow system and the problem of consumption of renewable power sources. For a multi-energy flow system, due to coupling of energy sources such as electricity, heat and the like, the cross-season soil heat storage is not only applied to a thermodynamic system, but also more used for stabilizing wind and light fluctuation brought by a renewable power source in an electric power system. The problem to be solved urgently is to establish a mathematical model suitable for the unified analysis and optimized scheduling of the multi-energy flow system by considering the coupling relation and interaction between the cross-season soil heat storage and other energy flows.
Disclosure of Invention
The invention provides a cross-season soil heat storage modeling method and device suitable for a multi-energy flow system, which comprehensively consider multi-structure layering inside the cross-season soil heat storage system and depict the internal physical operation state, so that the method and device are efficiently used for large-scale optimization of the multi-energy flow system.
In order to solve the above technical problem, in a first aspect, an embodiment of the present invention provides a cross-season soil heat storage modeling method suitable for a multi-energy flow system, including the following steps:
constructing a control equation of a three-dimensional temperature field of the surrounding soil of the vertical U-shaped pipe:
in the formula (1) (. Rho) S Is the density of the soil; c. C S The specific heat capacity of the soil is shown, a phi function represents the instantaneous state of energy in the soil, and a T function represents the instantaneous temperature of the soil; τ is time; z is the depth of layer; r is the distance from one point in the vertical U-shaped pipe to the shaft; theta is a tangential component; lambda is the conductivity of the soil;
determining initial conditions for solving the control equations:
T f (τ)＝T he (r,θ,τ)＝T s (r,θ,τ)＝T 0 (τ＝0) (2)
in the formula (2), T f Is the temperature, T, of the fluid in the vertical U-tube conduit he Is the temperature, T, of the soil heat exchanger s Is the temperature, T, of the soil surrounding the vertical U-tube 0 Is the initial condition of the temperature;
determining boundary conditions between the fluid in the vertical U-shaped pipe and the inner wall vertical U-shaped pipe:
in formula (3), λ he Is the thermal conductivity of the soil heat exchanger, alpha is the surface heat exchange coefficient of the inner wall of the vertical U-shaped pipe, r in Is vertical to the radius of the inner wall of the U-shaped pipe;
determining boundary conditions of the outer wall of the vertical U-shaped pipe and soil:
in the formula (4), r out Is perpendicular to the outer wall radius of the U-shaped pipe, lambda s Is the thermal conductivity of the soil;
the boundary conditions of the bottom edge ab of the outer wall of the vertical U-shaped pipe and the soil boundary, the side bc of the soil boundary and the side ef of the soil boundary are adiabatic conditions, and are expressed as follows:
the boundary condition of the outer wall of the vertical U-shaped pipe and the top edge cd of the soil boundary is a third type of boundary condition and is expressed as follows:
the boundary condition of the side ad of the vertical U-tube is a second type of boundary condition, which is expressed as:
in the formula (5), z 0 The depth of a layer for embedding the vertical U-shaped pipe; in the formula (6), r 0 The distance from the soil boundary to the center of the vertical U-shaped pipe is obtained; in the formula (7), h a Is the surface heat exchange coefficient of the U-shaped pipe and the air, T oe The ambient temperature; in the formula (8), r out Is vertical to the outer wall radius of the U-shaped pipe;
determining an energy equation of the soil heat storage system:
in the formula (9), the reaction mixture is,
for the available thermal energy stored in the soil thermal storage system at time t,
for the available thermal energy stored in the soil thermal storage system at time t-1,
for the functional rate of the soil thermal storage system at time t,
in order to achieve the heat storage efficiency of the soil heat storage system,
the heat release efficiency of the soil heat storage system; in formula (10), SQ s,0 For the total thermal energy stored in the soil thermal storage system at the initial moment,
total heat after energy storage for soil heat storage systemEnergy is saved;
determining the heat storage/release rate and the heat storage capacity limit of the soil heat storage system:
in the formula (11), the reaction mixture is,
is the maximum rate of heat storage; in the formula (12), the reaction mixture is,
is the maximum rate of heat release; in the formula (13), the reaction mixture is,
the minimum value of the heat energy which can be stored by the soil heat storage system,
the maximum value of the heat energy which can be stored by the soil heat storage system.
In order to solve the above technical problem, in a second aspect, an embodiment of the present invention provides a cross-season soil heat storage modeling apparatus suitable for a multi-energy flow system, including:
the control equation building module is used for building a control equation of a three-dimensional temperature field of the surrounding soil of the vertical U-shaped pipe:
in the formula (1), ρ S Is the density of the soil; c. C S The specific heat capacity of the soil is shown, a phi function represents the instantaneous state of energy in the soil, and T is the instantaneous temperature of the soil; τ is time; z is the depth of layer; r is the distance from one point in the vertical U-shaped pipe to the shaft; theta is a tangential component; lambda is the conductivity of the soil;
an initial condition determining module, configured to determine an initial condition for solving the control equation:
T f (τ)＝T he (r,θ,τ)＝T s (r,θ,τ)＝T 0 (τ＝0) (2)
in the formula (2), T f Is the temperature of the fluid in the vertical U-shaped pipe, T he Is the temperature, T, of the soil heat exchanger s Is the temperature of the soil around the vertical U-shaped pipe, T 0 Is the initial condition of temperature;
a first boundary condition determining module, configured to determine a boundary condition between a fluid in the vertical U-tube pipe and the inner wall vertical U-tube pipe:
in formula (3), λ he Alpha is the surface heat exchange coefficient of the vertical U-shaped pipe inner wall, r is the thermal conductivity of the soil heat exchanger in Is vertical to the radius of the inner wall of the U-shaped pipe;
the second boundary determining module is used for determining the boundary conditions between the outer wall of the vertical U-shaped pipe and the soil:
in the formula (4), r out Is perpendicular to the outer wall radius of the U-shaped pipe, lambda s Is the thermal conductivity of the soil;
the boundary conditions of the bottom edge ab of the outer wall of the vertical U-shaped pipe and the soil boundary, the side bc of the soil boundary and the side ef of the soil boundary are adiabatic conditions, and are expressed as follows:
the boundary condition of the outer wall of the vertical U-shaped pipe and the top edge cd of the soil boundary is a third type of boundary condition and is represented as follows:
the boundary conditions for the side ad of the vertical U-tube are of a second type, and are expressed as:
in the formula (5), z 0 The depth of a layer for embedding the vertical U-shaped pipe; in the formula (6), r 0 The distance from the soil boundary to the center of the vertical U-shaped pipe is obtained; in the formula (7), h a Is the surface heat exchange coefficient of the U-shaped pipe and the air, T oe Is the outside temperature; in the formula (8), r out The radius is vertical to the outer wall of the U-shaped pipe;
the capacity equation determining module is used for determining an energy equation of the soil heat storage system:
in the formula (9), the reaction mixture is,
for the available thermal energy stored in the soil thermal storage system at time t,
for the available thermal energy stored in the soil thermal storage system at time t-1,
for the functional rate of the soil heat storage system at time t,
in order to achieve the heat storage efficiency of the soil heat storage system,
the heat release efficiency of the soil heat storage system; in formula (10), SQ s,0 For the total thermal energy stored in the soil thermal storage system at the initial moment,
storing the total heat energy of the soil heat storage system after the energy storage is finished;
a heat storage capacity limit determination module for determining a heat storage/release rate and a soil heat storage system heat storage capacity limit:
in the formula (11), the reaction mixture is,
is the maximum rate of heat storage; in the formula (12), the reaction mixture is,
is the maximum rate of heat release; in the formula (13), the reaction mixture is,
the minimum value of the heat energy which can be stored by the soil heat storage system,
the maximum value of the heat energy which can be stored by the soil heat storage system.
In order to solve the foregoing technical problem, in a third aspect, an embodiment of the present invention provides a terminal device, including:
a memory for storing a computer program;
a processor for executing the computer program;
wherein the processor, when executing the computer program, implements a cross-season soil thermal storage modeling method suitable for a multi-energy flow system as described in the first aspect.
Compared with the prior art, the cross-season soil heat storage modeling method and device suitable for the multi-energy flow system, provided by the embodiment of the invention, have the beneficial effects that: the built model comprehensively considers multiple structural layers inside the cross-season soil heat storage system, describes the internal physical running state, and can be efficiently used for the large-scale optimization problem of the multi-energy flow system.
Drawings
In order to more clearly illustrate the technical features of the embodiments of the present invention, the drawings required to be used in the embodiments of the present invention will be briefly described below, and it is apparent that the drawings described below are only some embodiments of the present invention, and it is obvious to those skilled in the art that other drawings may be obtained based on the drawings without paying inventive labor.
FIG. 1 is a schematic flow diagram of a preferred embodiment of a cross-season soil heat storage modeling method suitable for a multi-energy flow system;
FIG. 2 is a schematic view of a model of a vertical U-shaped pipe and surrounding soil provided by the present invention;
fig. 3 is a schematic structural diagram of a preferred embodiment of a terminal device provided by the present invention.
Detailed Description
In order to make the technical features, objects and effects of the present invention more clearly understood, the following detailed description of the embodiments of the present invention is made with reference to the accompanying drawings and examples. The following examples are merely illustrative of the present invention and are not intended to limit the scope of the present invention. Other embodiments, which can be derived by a person skilled in the art from the embodiments of the present invention without inventive step, shall fall within the scope of the present invention.
Fig. 1 is a schematic flow chart of a preferred embodiment of a cross-season soil heat storage modeling method suitable for a multi-energy flow system provided by the invention.
As shown in fig. 1, the method comprises the steps of:
s10: constructing a control equation of a three-dimensional temperature field of soil around the vertical U-shaped pipe:
in the formula (1) (. Rho) S Is the density of the soil; c. C S The specific heat capacity of the soil is shown, a phi function represents the instantaneous state of energy in the soil, and T is the instantaneous temperature of the soil; τ is time; z is the depth of layer; r is the distance from one point in the vertical U-shaped pipe to the shaft; theta is a tangential component; λ is the conductivity of the soil.
It should be noted that, a typical seasonal soil heat storage heating air conditioning system is composed of four parts: ceiling solar collection systems, soil heat storage systems (also known as underground embedded heat exchange systems), heat pump systems and floor radiant heating systems (also known as air conditioning terminal units). A planar heat exchanger is adopted for heat exchange between the ceiling solar energy collecting system and the soil heat storage system and between the solar hot water direct heating element and the floor radiation heating system, and the soil heat storage system adopts a vertical U-shaped pipe (also called a vertical U-shaped heat exchanger) buried in deep soil for heat exchange.
The vertical U-shaped pipe is the main component of the system, and the heat transfer effect of the vertical U-shaped pipe directly influences the coefficient of performance (COP). The actual heat transfer between the vertical U-tubes and the surrounding soil is a complex and unstable process. Necessary simplifications must be made to facilitate the analysis of the problem. The physical model is simplified as follows:
1) Conductivity changes due to thermal coupling wet transfer effects and ground water convective heat transfer effects are ignored. The heat transfer between the soil and the embedded exchanger is considered to be purely conductive. The soil was stratified according to depth. The conductivity of each layer is constant.
2) There is no idea of heat transfer between the bottom of the vertical well and the soil. The boundary may be considered adiabatic.
3) The thermophysical parameters of the soil, cement slurry, buried pipe and fluid in the pipe do not change as heat transfer occurs.
Theoretically, the impact of a buried pipe on the surrounding soil can reach a limited area. In fact, the influence of the vertical U-shaped pipe on the soil temperature field is smaller and smaller with the increase of the distance, and the influence can be ignored at a far place.
According to the above assumptions, the control equation of the three-dimensional temperature field of the surrounding soil of the vertical U-shaped pipe can be obtained according to equation (1).
S20: determining initial conditions for solving the control equations:
T f (τ)＝T he (r,θ,τ)＝T s (r,θ,τ)＝T 0 (τ＝0) (2)
in the formula (2), T f Is the temperature of the fluid in the vertical U-shaped pipe, T he Temperature, T, of soil heat exchangers s Is the temperature of the soil around the vertical U-shaped pipe, T 0 Is the initial condition of the temperature.
After the control equation of the three-dimensional temperature field of the soil around the vertical U-shaped pipe is obtained, the initial condition of the soil is determined according to the equation (2).
S30: determining boundary conditions between the fluid in the vertical U-shaped pipe and the inner wall vertical U-shaped pipe:
in formula (3), λ he Is the thermal conductivity of the soil heat exchanger, alpha is the surface heat exchange coefficient of the inner wall of the vertical U-shaped pipe, r in Is vertical to the radius of the inner wall of the U-shaped pipe.
S40: determining boundary conditions of the outer wall of the vertical U-shaped pipe and soil:
in the formula (4), r out Is perpendicular to the outer wall radius of the U-shaped pipe, lambda s Is the thermal conductivity of the soil;
the boundary conditions of the bottom edge ab of the outer wall of the vertical U-shaped pipe and the soil boundary, the side edge bc of the soil boundary and the side edge ef of the soil boundary are thermal insulation conditions, and are represented as follows:
the boundary condition of the outer wall of the vertical U-shaped pipe and the top edge cd of the soil boundary is a third type of boundary condition and is represented as follows:
the boundary conditions for the side ad of the vertical U-tube are of a second type, and are expressed as:
in the formula (5), z 0 The depth of layer is embedded in the vertical U-shaped pipe; in the formula (6), r 0 The distance from the soil boundary to the center of the vertical U-shaped pipe is obtained; in the formula (7), h a Is U-shapedSurface heat transfer coefficient of tube and air, T oe The ambient temperature; in the formula (8), r out Is vertical to the outer wall radius of the U-shaped pipe.
After the initial conditions are determined, the boundary conditions of the vertical U-shaped tube and the surrounding soil need to be determined, wherein a model schematic diagram of the vertical U-shaped tube and the surrounding soil is shown in fig. 2.
Based on the previous analysis and some assumptions, a seasonal soil thermal storage heating air conditioning system control model that can determine when and how to perform thermal storage/release can be proposed according to equations (9) - (10). In this model, the decision variable is the heat energy flowing into the soil thermal storage system per time period. The available thermal energy and the corresponding heat loss in the soil heat storage system are changed simultaneously. The equation represents the thermal energy stored in the storage unit at time t. The equation requires that the total amount of thermal energy for the last moment is given in advance, usually the same as the initial thermal energy storage amount.
S50: determining an energy equation of a soil heat storage system:
in the formula (9), the reaction mixture is,
for the available thermal energy stored in the soil thermal storage system at time t,
for the available thermal energy stored in the soil thermal storage system at time t-1,
for the functional rate of the soil thermal storage system at time t,
for the heat storage efficiency of the soil heat storage system,
the heat release efficiency of the soil heat storage system; in formula (10), SQ s,0 For the total thermal energy stored in the soil thermal storage system at the initial moment,
and storing the total heat energy after the energy storage for the soil heat storage system is finished.
S60: determining the heat storage/release rate and the heat storage capacity limit of the soil heat storage system:
in the formula (11), the reaction mixture is,
is the maximum rate of heat storage; in the formula (12), the reaction mixture is,
is the maximum rate of heat release; in the formula (13), the reaction mixture is,
the minimum value of the heat energy which can be stored by the soil heat storage system,
the maximum value of the heat energy which can be stored by the soil heat storage system.
Wherein, the boundary condition between the outer wall of the vertical U-shaped pipe and the soil determines the heightEfficiency between heat exchange between solar energy collection system and soil heat storage system
And
the energy storage and energy supply efficiency of the soil energy storage system is shown and is related to the equipment characteristics of specific device parameters.
According to the cross-season soil heat storage modeling method suitable for the multi-energy flow system, the built model comprehensively considers multi-structure layering inside the cross-season soil heat storage system, the internal physical operation state is described, and the method can be effectively used for large-scale optimization of the multi-energy flow system.
In a preferred embodiment, the method further comprises:
s70: and determining the optimal heat storage time and the optimal heat release time of the soil heat storage system by taking the established model as constraint and the lowest operation cost as a target.
After the soil heat storage model is constructed, the optimal heat storage time and the optimal heat release time of the soil heat storage system can be rapidly determined through a mature commercial optimization solver, and a reference basis is provided for heat storage/release.
It should be understood that all or part of the processes of the cross-season soil heat storage modeling method for the multi-energy flow system may be implemented by a computer program, which may be stored in a computer readable storage medium and executed by a processor, so as to implement the steps of the cross-season soil heat storage modeling method for the multi-energy flow system. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, an executable file or some intermediate form, etc. The computer readable medium may include: any entity or device capable of carrying the computer program code, recording medium, U disk, removable hard disk, magnetic disk, optical disk, computer Memory, read-Only Memory (ROM), random Access Memory (RAM), electrical carrier wave signal, telecommunications signal, and software distribution medium. It should be noted that the computer-readable medium may contain suitable additions or subtractions depending on the requirements of legislation and patent practice in the jurisdiction, for example in some jurisdictions, in accordance with legislation and patent practice, the computer-readable medium may not include electrical carrier signals or telecommunication signals.
Correspondingly, the invention provides a cross-season soil heat storage and modeling device suitable for a multi-energy flow system, which comprises:
the control equation building module is used for building a control equation of a three-dimensional temperature field of the surrounding soil of the vertical U-shaped pipe:
in the formula (1), ρ S Is the density of the soil; c. C S The specific heat capacity of the soil is shown, a phi function represents the instantaneous state of energy in the soil, and T is the instantaneous temperature of the soil; τ is time; z is the depth of layer; r is the distance from one point in the vertical U-shaped pipe to the shaft; theta is a tangential component; lambda is the conductivity of the soil;
an initial condition determining module, configured to determine an initial condition for solving the control equation:
T f (τ)＝T he (r,θ,τ)＝T s (r,θ,τ)＝T 0 (τ＝0) (2)
in the formula (2), T f Is the temperature, T, of the fluid in the vertical U-tube conduit he Is the temperature, T, of the soil heat exchanger s Is the temperature of the soil around the vertical U-shaped pipe, T 0 Is the initial condition of the temperature;
a first boundary condition determining module, configured to determine a boundary condition between a fluid in the vertical U-tube and the inner-wall vertical U-tube:
in formula (3), λ he Is the thermal conductivity of the soil heat exchanger, alpha is the surface heat exchange coefficient of the inner wall of the vertical U-shaped pipe, r in Is vertical to the radius of the inner wall of the U-shaped pipe;
the second boundary determining module is used for determining the boundary condition between the outer wall of the vertical U-shaped pipe and soil:
in the formula (4), r out Is perpendicular to the outer wall radius of the U-shaped pipe, lambda s Is the thermal conductivity of the soil;
the boundary conditions of the bottom edge ab of the outer wall of the vertical U-shaped pipe and the soil boundary, the side bc of the soil boundary and the side ef of the soil boundary are adiabatic conditions, and are expressed as follows:
the boundary condition of the outer wall of the vertical U-shaped pipe and the top edge cd of the soil boundary is a third type of boundary condition and is represented as follows:
the boundary condition of the side ad of the vertical U-tube is a second type of boundary condition, which is expressed as:
in the formula (5), z 0 The depth of a layer for embedding the vertical U-shaped pipe; in the formula (6), r 0 The distance from the boundary of the soil to the center of the vertical U-shaped pipe; formula (A), (B) and7) In (h) a Is the surface heat exchange coefficient of the U-shaped pipe and the air, T oe Is the outside temperature; in the formula (8), r out Is vertical to the outer wall radius of the U-shaped pipe;
the capacity equation determining module is used for determining an energy equation of the soil heat storage system:
in the formula (9), the reaction mixture is,
for the available thermal energy stored in the soil thermal storage system at time t,
for the available thermal energy stored in the soil thermal storage system at time t-1,
for the functional rate of the soil heat storage system at time t,
in order to achieve the heat storage efficiency of the soil heat storage system,
the heat release efficiency of the soil heat storage system; in formula (10), SQ s,0 For the total thermal energy stored in the soil thermal storage system at the initial moment,
storing the total heat energy of the soil heat storage system after the energy storage is finished;
a heat storage capacity limit determination module for determining a heat storage/release rate and a heat storage capacity limit of the soil heat storage system:
in the formula (11), the reaction mixture is,
is the maximum rate of heat storage; in the formula (12), the reaction mixture is,
is the maximum rate of heat release; in the formula (13), the reaction mixture is,
the minimum value of the heat energy which can be stored by the soil heat storage system,
the maximum value of the heat energy which can be stored by the soil heat storage system.
Fig. 3 is a schematic structural diagram of a preferred embodiment of a terminal device provided by the present invention, where the terminal device is capable of implementing all processes of the cross-season soil heat storage modeling method applicable to a multi-energy flow system and achieving corresponding technical effects.
As shown in fig. 3, the apparatus includes:
a memory 21 for storing a computer program;
a processor 22 for executing the computer program;
when the processor 22 executes the computer program, the method for modeling the trans-seasonal soil heat storage suitable for the multi-energy flow system according to the embodiment is realized.
Illustratively, the computer program may be divided into one or more modules/units, which are stored in the memory 21 and executed by the processor 22 to accomplish the present invention. The one or more modules/units may be a series of computer program instruction segments capable of performing specific functions, which are used for describing the execution process of the computer program in the terminal device.
The Processor 22 may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic device, discrete hardware component, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory 21 can be used for storing the computer programs and/or modules, and the processor 22 implements various functions of the terminal device by running or executing the computer programs and/or modules stored in the memory 21 and calling data stored in the memory 21. The memory 21 may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required by at least one function (such as a sound playing function, an image playing function, etc.), and the like; the storage data area may store data (such as audio data, a phonebook, etc.) created according to the use of the cellular phone, and the like. In addition, the memory 21 may include high speed random access memory, and may also include non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), at least one magnetic disk storage device, a Flash memory device, or other volatile solid state memory device.
It should be noted that the terminal device includes, but is not limited to, a processor and a memory, and those skilled in the art will understand that the structural diagram in fig. 3 is only an example of the terminal device, and does not constitute a limitation to the terminal device, and may include more components than those shown in the figure, or may combine some components, or may include different components.
The above description is only a preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and it should be noted that, for those skilled in the art, it is possible to make several equivalent obvious modifications and/or equivalent substitutions without departing from the technical principle of the present invention, and these obvious modifications and/or equivalent substitutions should also be regarded as the scope of the present invention.
Claims (3)
1. A cross-season soil heat storage modeling method suitable for a multi-energy flow system is characterized by comprising the following steps:
constructing a control equation of a three-dimensional temperature field of the surrounding soil of the vertical U-shaped pipe:
in the formula (1), ρ S Is the density of the soil; c. C S The specific heat capacity of the soil is shown, a phi function represents the instantaneous state of energy in the soil, and a T function represents the instantaneous temperature of the soil; τ is time; z is the depth of layer; r is the distance from one point in the vertical U-shaped pipe to the shaft; theta is a tangential component; lambda is the conductivity of the soil;
determining initial conditions for solving the governing equations:
T f (τ)＝T he (r,θ,τ)＝T s (r,θ,τ)＝T 0 (τ＝0) (2)
in formula (2), T f Is the temperature of the fluid in the vertical U-shaped pipe, T he Is the temperature of the soil heat exchanger, T s Is the temperature of the soil around the vertical U-shaped pipe, T 0 Is the initial condition of temperature;
determining boundary conditions between the fluid in the vertical U-shaped pipe and the inner wall vertical U-shaped pipe:
in formula (3), λ he Is the thermal conductivity of the soil heat exchanger, alpha is the surface heat exchange coefficient of the inner wall of the vertical U-shaped pipe, r in Is vertical to the radius of the inner wall of the U-shaped pipe;
determining boundary conditions of the outer wall of the vertical U-shaped pipe and soil:
in the formula (4), r out Is perpendicular to the outer wall radius of the U-shaped pipe, lambda s Is the thermal conductivity of the soil;
the boundary conditions of the bottom edge ab of the outer wall of the vertical U-shaped pipe and the soil boundary, the side edge bc of the soil boundary and the side edge ef of the soil boundary are thermal insulation conditions, and are represented as follows:
the boundary condition of the outer wall of the vertical U-shaped pipe and the top edge cd of the soil boundary is a third type of boundary condition and is expressed as follows:
the boundary condition of the side ad of the vertical U-tube is a second type of boundary condition, which is expressed as:
in the formula (5), z 0 The depth of layer is embedded in the vertical U-shaped pipe; in the formula (6), r 0 The distance from the soil boundary to the center of the vertical U-shaped pipe; in the formula (7), h a Is the surface heat exchange coefficient of the U-shaped pipe and the air, T oe Is the outside temperature; in the formula (8), r out Is vertical to the radius of the outer wall of the U-shaped pipe;
determining an energy equation of the soil heat storage system:
in the formula (9), the reaction mixture is,
for the available thermal energy stored in the soil thermal storage system at time t,
the available thermal energy stored in the soil thermal storage system for time t-1,
for the functional rate of the soil heat storage system at time t,
in order to achieve the heat storage efficiency of the soil heat storage system,
the heat release efficiency of the soil heat storage system; in formula (10), SQ s,0 For the total thermal energy stored in the soil thermal storage system at the initial moment,
storing the total heat energy of the soil heat storage system after the energy storage is finished;
determining the heat storage/release rate and the heat storage capacity limit of the soil heat storage system:
in the formula (11), the reaction mixture is,
is the maximum rate of heat storage; in the formula (12), the reaction mixture is,
is the maximum rate of heat release; in the formula (13), the reaction mixture is,
the minimum value of the heat energy which can be stored by the soil heat storage system,
the maximum value of the heat energy which can be stored by the soil heat storage system.
2. A cross-season soil heat storage modeling device suitable for a multi-energy flow system, comprising:
the control equation building module is used for building a control equation of a three-dimensional temperature field of the surrounding soil of the vertical U-shaped pipe:
in the formula (1) (. Rho) S Is the density of the soil; c. C S The specific heat capacity of the soil is shown, a phi function represents the instantaneous state of energy in the soil, and a T function represents the instantaneous temperature of the soil; τ is time; z is the depth of layer; r is the distance from one point in the vertical U-shaped pipe to the shaft; theta is a tangential component; lambda is the conductivity of the soil;
an initial condition determining module, configured to determine an initial condition for solving the control equation:
T f (τ)＝T he (r,θ,τ)＝T s (r,θ,τ)＝T 0 (τ＝0) (2)
in formula (2), T f Is the temperature, T, of the fluid in the vertical U-tube conduit he Is the temperature of the soil heat exchanger, T s Is the temperature of the soil around the vertical U-shaped pipe, T 0 Is the initial condition of temperature;
the first boundary condition determining module is used for determining the boundary condition between the fluid in the vertical U-shaped pipe pipeline and the inner wall vertical U-shaped pipe pipeline:
in formula (3), λ he Alpha is the surface heat exchange coefficient of the inner wall of the vertical U-shaped pipe, r is the thermal conductivity of the soil heat exchanger in Is vertical to the radius of the inner wall of the U-shaped pipe;
the second boundary determining module is used for determining the boundary conditions between the outer wall of the vertical U-shaped pipe and the soil:
in the formula (4), r out Is perpendicular to the outer wall radius of the U-shaped pipe, lambda s Is the thermal conductivity of the soil;
the boundary conditions of the bottom edge ab of the outer wall of the vertical U-shaped pipe and the soil boundary, the side bc of the soil boundary and the side ef of the soil boundary are adiabatic conditions, and are expressed as follows:
the boundary condition of the outer wall of the vertical U-shaped pipe and the top edge cd of the soil boundary is a third type of boundary condition and is expressed as follows:
the boundary conditions for the side ad of the vertical U-tube are of a second type, and are expressed as:
in the formula (5), z 0 The depth of a layer for embedding the vertical U-shaped pipe; in the formula (6), r 0 The distance from the soil boundary to the center of the vertical U-shaped pipe; in the formula (7), h a Is the surface heat exchange coefficient of the U-shaped pipe and the air, T oe Is the outside temperature; in the formula (8), r out Is vertical to the radius of the outer wall of the U-shaped pipe;
the capacity equation determining module is used for determining an energy equation of the soil heat storage system:
in the formula (9), the reaction mixture is,
for the available thermal energy stored in the soil thermal storage system at time t,
the available thermal energy stored in the soil thermal storage system for time t-1,
for the functional rate of the soil heat storage system at time t,
in order to achieve the heat storage efficiency of the soil heat storage system,
the heat release efficiency of the soil heat storage system; in formula (10), SQ s,0 For the total thermal energy stored in the soil thermal storage system at the initial moment,
storing the total heat energy of the soil heat storage system after the energy storage is finished;
a heat storage capacity limit determination module for determining a heat storage/release rate and a soil heat storage system heat storage capacity limit:
in the formula (11), the reaction mixture is,
is the maximum rate of heat storage; in the formula (12), the reaction mixture is,
is the maximum rate of heat release; in the formula (13), the reaction mixture is,
the minimum value of the heat energy which can be stored by the soil heat storage system,
the maximum value of the heat energy which can be stored by the soil heat storage system.
3. A terminal device, comprising:
a memory for storing a computer program;
a processor for executing the computer program;
wherein the processor, when executing the computer program, implements the cross-season soil thermal storage modeling method for a multi-energy flow system of claim 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111028987.4A CN113836700B (en) | 2021-09-02 | 2021-09-02 | Cross-season soil heat storage modeling method and device suitable for multi-energy flow system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111028987.4A CN113836700B (en) | 2021-09-02 | 2021-09-02 | Cross-season soil heat storage modeling method and device suitable for multi-energy flow system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113836700A CN113836700A (en) | 2021-12-24 |
| CN113836700B true CN113836700B (en) | 2022-12-06 |
Family
ID=78962106
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111028987.4A Active CN113836700B (en) | 2021-09-02 | 2021-09-02 | Cross-season soil heat storage modeling method and device suitable for multi-energy flow system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113836700B (en) |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110059389B (en) * | 2019-04-10 | 2023-04-07 | 北京石油化工学院 | Solar cross-season soil heat storage POD (POD) rapid prediction method |
| CN110489912B (en) * | 2019-08-27 | 2022-12-13 | 北京石油化工学院 | Method for simulating numerical value of solar cross-season soil heat storage layered slice |
| CN111125921A (en) * | 2019-12-27 | 2020-05-08 | 常州工学院 | Method for rapidly and accurately realizing dynamic simulation of performance of vertical U-shaped ground heat exchanger |
-
2021
- 2021-09-02 CN CN202111028987.4A patent/CN113836700B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN113836700A (en) | 2021-12-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liu et al. | Performance and feasibility study of hybrid ground source heat pump system assisted with cooling tower for one office building based on one Shanghai case | |
| Soni et al. | Ground coupled heat exchangers: A review and applications | |
| Fan et al. | Thermal performance and operation strategy optimization for a practical hybrid ground-source heat-pump system | |
| CN102550340B (en) | Solar regenerative heating device for greenhouse heating | |
| Al-Ameen et al. | Thermal performance of a solar assisted horizontal ground heat exchanger | |
| Wong et al. | Recent inter-seasonal underground thermal energy storage applications in Canada | |
| CN215675924U (en) | Double-heat-storage heat pump heating device for middle-deep geothermal buried pipe | |
| Sun et al. | Research on multi-objective optimization of control strategies and equipment parameters for a combined heating system of geothermal and solar energy in cold and arid regions based on TRNSYS | |
| Cho et al. | Simulation modeling of ground source heat pump systems for the performance analysis of residential buildings | |
| CN113836700B (en) | Cross-season soil heat storage modeling method and device suitable for multi-energy flow system | |
| Aisa et al. | Modelling and simulation of a solar water heating system with thermal storage | |
| Mourad et al. | An energy-efficient smart home for new cities in Egypt | |
| CN206269225U (en) | A kind of many power heating systems | |
| Assad et al. | Heating and cooling by geothermal energy | |
| Hwang et al. | Transient modeling of seasonal borehole thermal energy storage system during heat energy storing process | |
| Emmi et al. | A Management Strategy for Multi-Source Heat Pump Systems | |
| Woo et al. | Performance investigation of the hybrid-renewable energy system with geothermal and solar heat sources for a residential building in South Korea | |
| Matysko | Ground as a source of heat for heat pumps | |
| Liu et al. | Storage periods and heating terminals types influence on the solar heating thermal storage volume in the Tibetan Plateau regions | |
| CN116293866A (en) | Operation control method and system of medium-deep geothermal coupling heat storage and supply system | |
| CN219160452U (en) | Utilize ground to carry out adjusting device who strides middle-low temperature heat accumulation in season | |
| Wang et al. | Energy optimization in the design and well layout of Shengjing Shopping Mall in Tiexi District, Shenyang City | |
| Fang et al. | The Ground-Coupled Heat Pump Technology in China | |
| Varney et al. | Simulations of a photovoltaic-thermal ground source heat pump system | |
| Eslami Nejad et al. | Heat pump capacity effects on peak electricity consumption and total length of self-and solar-assisted shallow ground heat exchanger networks |
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 |