CN-120706300-B - Method and system for analyzing thermodynamic characteristics of compressed wet air of gas storage

CN120706300BCN 120706300 BCN120706300 BCN 120706300BCN-120706300-B

Abstract

The invention relates to a thermodynamic characteristic analysis method and a thermodynamic characteristic analysis system for compressed humid air of a gas storage, which are used for establishing a thermodynamic theoretical model comprising a flow model, a non-isothermal flow model, a damp-heat coupling model and an actual gas state equation, further establishing a three-dimensional gas storage numerical model comprising a fluid domain, a sealing layer lining domain and a surrounding rock domain, endowing the fluid domain with parameters of humid air materials, endowing the solid domain with parameters of concrete and rock materials, setting initial and boundary conditions, and finally solving the model to obtain a spatial variation process of the temperature, the pressure and the humidity distribution of the compressed humid air in the gas storage. The method can effectively simulate the space change process of thermodynamic characteristics of compressed air in the power station gas storage in consideration of the influence of moisture content in the air, and can be used for design calculation of the underground gas storage of the compressed air energy storage power station and design specification standard formulation of the underground gas storage of the compressed air energy storage power station.

Inventors

Wan Fa
YAO ZHICHAO
HUANG RUIYI
LI HAIFENG
WANG WEIGUANG
LI JIALEI
Wu Gongzheng
LIN JIAYU
Jian Chaofeng
LIU LI

Assignees

珠江水利委员会珠江水利科学研究院

Dates

Publication Date: 20260508
Application Date: 20250612

Claims (8)

1.A method for analyzing thermodynamic characteristics of compressed humid air in a gas storage, comprising: S1, establishing a thermodynamic theoretical model of compressed wet air of a gas storage, wherein the thermodynamic theoretical model comprises the following steps: the flow model adopts a turbulence k-epsilon model and is used for simulating the turbulence flow characteristics of the gas in the gas storage; and the non-isothermal flow model is used for coupling temperature change in the compression and expansion processes of the gas and is used for representing the flow behavior of the gas under the non-isothermal condition, and the expression of the non-isothermal flow model is as follows: Wherein Q is a heat conduction term, cp is a specific heat capacity, T is a thermodynamic temperature, k is a thermal conductivity, Q is a heat source term, Q p is an energy source term caused by pressure change, and Q vd is viscous dissipation energy; Wherein, the expression of the energy source term Q p caused by the pressure change is: the expression of viscous dissipated energy Q vd is: The boundary wall heat exchange-n.q has the expression: Wherein alpha p is the expansion coefficient, p A is the pressure, tau is the stress tensor, n is the normal vector of the wall surface and points to the inside of the fluid, u τ is the friction speed, tw is the wall surface temperature, and T + is the dimensionless temperature; The wet-heat coupling model is used for describing the humidity diffusion, phase change and latent heat exchange processes and comprises an air moisture transport equation, a turbulent water flow diversion equation and a wet air phase change and temperature coupling equation; Wherein, the expression of the air moisture transport equation is: Wherein M V is the molar mass of the gas phase substance, c v is the gas phase concentration, G w is the gas phase substance transmission quantity per unit time passing unit area, and comprises diffusion and convection contributions; the expression of the water flow equation in turbulence is: Wherein n is the normal vector of an interface unit and points to a gas phase region, C v,w is the vapor equilibrium concentration in a liquid water phase, C v is the actual vapor concentration in the gas phase, and phi w + is a dimensionless concentration parameter; the expression of the wet air phase transition and temperature coupling equation is: Wherein Q evap is the evaporation latent heat flux, which represents the heat taken away in unit time and unit area in the evaporation process, L v is the vaporization latent heat, which represents the heat required for converting unit mass of liquid into gas, Q evap is the phase change mass flux, Q is the total heat flux, which comprises the comprehensive effects of the phase change latent heat and sensible heat transfer, C p,v is the constant pressure specific heat capacity of water vapor, C p,a is the constant pressure specific heat capacity of dry air, C p,l is the constant pressure specific heat capacity of liquid water, and g lc is the liquid water diffusion flux vector; the actual gas state equation is used for representing thermodynamic characteristics of the humid air; S2, constructing a three-dimensional gas storage numerical model comprising a fluid domain, a sealing layer lining domain and a surrounding rock domain based on the thermodynamic theoretical model; s3, wet air material parameters are given to the fluid domain of the three-dimensional gas storage numerical model, concrete and rock material parameters are given to the solid domain, and initial conditions and boundary conditions of the three-dimensional gas storage numerical model are set; And S4, solving the three-dimensional gas storage numerical model to obtain a spatial variation process of temperature, pressure and humidity distribution of compressed humid air in the gas storage.
2. The method for analyzing thermodynamic characteristics of compressed humid air of a gas storage according to claim 1, wherein constructing a three-dimensional gas storage numerical model including a fluid domain, a sealing layer lining domain and a surrounding rock domain based on the thermodynamic theoretical model comprises: According to the gas storage structure, the gas filling and discharging port and the gas storage cavity are used as fluid domains for containing compressed wet air, the sealing layer and the lining structure are used as sealing layer lining domains for sealing and protecting, surrounding rock is used as surrounding rock domains for providing external support, and the domains are numbered; The method comprises the steps of simulating turbulent flow characteristics of gas in a fluid domain through the flow model, coupling temperature change through a non-isothermal flow model to represent flow behavior under a non-isothermal condition, describing humidity diffusion, phase change and latent heat exchange process through a damp-heat coupling model, and representing thermodynamic characteristics of humid air in the fluid domain through an actual gas state equation.
3. The method for analyzing thermodynamic characteristics of compressed humid air of a gas storage according to claim 1, wherein before solving the three-dimensional gas storage numerical model, the method further comprises: and carrying out grid division on the three-dimensional gas storage numerical model, wherein a fluid domain adopts a high-precision boundary layer grid.
4. The method for analyzing thermodynamic characteristics of compressed humid air of air storage according to claim 1, wherein the setting of initial conditions and boundary conditions of the three-dimensional air storage numerical model comprises: setting the initial pressure of the fluid domain to be the ambient atmospheric pressure or the design pressure according to whether the initial cushioning process is distinguished or not; And dynamically adjusting pressure, flow speed and heat flux according to the operation charging and discharging process of the power station, and setting corresponding boundary conditions according to the flow model, the non-isothermal flow model and the air moisture transportation equation.
5. The method for analyzing thermodynamic characteristics of compressed humid air in a gas storage according to claim 1, wherein the flow model comprises a fluid continuity equation, a turbulent kinetic energy k equation, and a turbulent dissipation ratio ε equation, wherein: the expression of the fluid continuity equation is: Wherein ρ is the gas density, t is the time, v is the hamiltonian, u is the velocity vector matrix; the expression of the turbulent kinetic energy k equation is: Wherein k is turbulence kinetic energy, mu is dynamic viscosity, mu T is turbulence viscosity coefficient, and p K is turbulence kinetic energy generation term; The expression of the turbulent dissipation ratio ε equation is: Where σ ε 、C ε1 and C ε2 are both turbulence model constants.
6. A system for analyzing thermodynamic properties of compressed wet air in a gas storage, the system being adapted to perform the method for analyzing thermodynamic properties of compressed wet air in a gas storage according to any one of claims 1 to 5, the system comprising: the thermodynamic theory modeling module is used for establishing a thermodynamic theory model of compressed humid air of the air storage, wherein the thermodynamic theory model comprises the following components: the flow model adopts a turbulence k-epsilon model and is used for simulating the turbulence flow characteristics of the gas in the gas storage; A non-isothermal flow model for coupling temperature changes in the compression and expansion processes of the gas and for characterizing the flow behavior of the gas under non-isothermal conditions; the damp-heat coupling model is used for coupling an air moisture transport equation and describing the humidity diffusion, phase change and latent heat exchange process; the actual gas state equation is used for representing thermodynamic characteristics of the humid air; the three-dimensional model building module is used for building a three-dimensional gas storage numerical model comprising a fluid domain, a sealing layer lining domain and a surrounding rock domain based on the thermodynamic theoretical model; The parameter configuration module is used for giving wet air material parameters to the fluid domain of the three-dimensional gas storage numerical model, giving concrete and rock material parameters to the solid domain, and setting initial conditions and boundary conditions of the three-dimensional gas storage numerical model; and the model solving module is used for solving the three-dimensional gas storage numerical model to obtain a spatial variation process of temperature, pressure and humidity distribution of compressed humid air in the gas storage.
7. An electronic device is provided, which comprises a first electronic device, Characterized by comprising a memory and a processor, wherein, The memory is used for storing programs; the processor is coupled to the memory for executing the program stored in the memory to implement the steps in the method for analyzing compressed humid air thermodynamic characteristics of a gas storage according to any one of the preceding claims 1 to 5.
8. A computer readable storage medium storing a computer readable program or instructions which when executed by a processor is capable of carrying out the steps of the method for analyzing compressed humid air thermodynamic properties of a gas storage according to any one of the preceding claims 1 to 5.

Description

Method and system for analyzing thermodynamic characteristics of compressed wet air of gas storage Technical Field The invention relates to the technical field of compressed air energy storage, in particular to a method and a system for analyzing thermodynamic characteristics of compressed wet air of a gas storage. Background Worldwide there is an increasing demand for renewable energy and energy storage technologies to address the challenges of climate change and energy security. The compressed air energy storage power station compresses and stores air by utilizing electric power, and then releases compressed air to drive the generator to generate power when needed, so that the high-efficiency storage and utilization of energy are realized. The technology has the advantages of low cost, environmental protection, strong flexibility and the like, and is increasingly concerned and supported by governments, enterprises and investors. The development of compressed air energy storage power stations provides a new solution to the popularization of clean energy sources and the challenges of coping with power systems. The compressed air energy storage power station consists of six parts, namely a compressor, an expansion turbine, a generator, a gas storage, a control system and a heat storage system. Large-scale compressed air energy storage power stations use underground cavities as reservoirs, wherein the thermodynamic properties of the stored compressed air relate to the safety stability and economic viability of the whole power station. In the design planning stage of the compressed air energy storage power station and the selection of equipment, the thermodynamic parameter change process of the compressed air in the operation process must be clearly and definitely predicted. The current calculation method of the thermodynamic process of the compressed air mainly comprises an analysis method and a numerical method. The analysis method calculates the temperature and the pressure by fixing the parameters (such as the storage capacity and the heat exchange area) of the gas storage and combining the analysis solution of the gas charging and discharging process, and the numerical method predicts the change of thermodynamic parameters by simulating the gas flowing state in the gas storage. However, none of the prior art considers the effect of natural air or warehouse entry gas moisture content on thermodynamic properties. In the prior art, the influence of humidity in the air is ignored, so that the calculation result has deviation. The change of the latent heat of phase change and heat transfer characteristic caused by humidity is not included in the model, so that temperature and pressure prediction errors are caused, wet air can be condensed or frozen in the compression/expansion process, the service life of equipment and the stability of a system are threatened, and meanwhile, the overall efficiency of the energy storage system is reduced due to the fact that humidity related parameters (such as convection heat transfer coefficient) are not optimized. Disclosure of Invention In view of the foregoing, it is necessary to provide a method and a system for analyzing thermodynamic characteristics of compressed humid air in a gas storage, so as to solve the problems of low detection accuracy, background complexity and uneven cell distribution in nasopharyngeal carcinoma organoid images in the existing target detection technology. In order to solve the above problems, in a first aspect, an embodiment of the present invention provides a method for analyzing thermodynamic characteristics of compressed humid air in a gas storage, including: S1, establishing a thermodynamic theoretical model of compressed wet air of a gas storage, wherein the thermodynamic theoretical model comprises the following steps: the flow model adopts a turbulence k-epsilon model and is used for simulating the turbulence flow characteristics of the gas in the gas storage; A non-isothermal flow model for coupling temperature changes in the compression and expansion processes of the gas and for characterizing the flow behavior of the gas under non-isothermal conditions; the damp-heat coupling model is used for coupling an air moisture transport equation and describing the humidity diffusion, phase change and latent heat exchange process; the actual gas state equation is used for representing thermodynamic characteristics of the humid air; S2, constructing a three-dimensional gas storage numerical model comprising a fluid domain, a sealing layer lining domain and a surrounding rock domain based on the thermodynamic theoretical model; s3, wet air material parameters are given to the fluid domain of the three-dimensional gas storage numerical model, concrete and rock material parameters are given to the solid domain, and initial conditions and boundary conditions of the three-dimensional gas storage numerical model are set; And S4, solving the three-dim