Gas storage compressed humid air thermodynamic property analysis method and system

Abstract

The invention relates to a gas storage compressed humid air thermodynamic property analysis method and system. The method comprises the following steps: establishing a thermodynamic theoretical model containing a flow model, a non-isothermal flow model, a humid heat coupling model and an actual gas state equation; further constructing a three-dimensional gas storage numerical model containing a fluid domain, a sealing layer lining domain and a surrounding rock domain; moist air material parameters are assigned to the fluid domain, concrete and rock material parameters are assigned to the solid domain, and initial and boundary conditions are set; and finally, solving the model to obtain the space change process of the temperature, pressure and humidity distribution of the compressed humid air in the gas storage. The method can effectively simulate the spatial change process of the thermodynamic characteristics of the compressed air in the power station gas storage in consideration of the influence of the moisture content in the air, can be used for design calculation of the underground gas storage of the compressed air energy storage power station, and can also be used for formulating design specifications and standards of the underground gas storage of the compressed air energy storage power station.

Description

Technical Field

[0001] The present invention relates to the technical field of compressed air energy storage, and in particular to a method and system for analyzing the thermodynamic characteristics of compressed wet air in a gas storage reservoir. Background Art

[0002] Global demand for renewable energy and energy storage technologies is growing to address climate change and energy security challenges. Compressed air energy storage (CAES) power plants utilize electricity to compress and store air, then release the compressed air when needed to drive a generator for power generation, achieving efficient energy storage and utilization. This technology, with its low cost, environmental friendliness, and flexibility, is gaining increasing attention and support from governments, businesses, and investors. The development of CAES provides a new solution for promoting clean energy and addressing power system challenges.

[0003] A compressed air energy storage power station consists of six components: a compressor, an expansion turbine, a generator, a gas storage reservoir, a control system, and a heat storage system. Large-scale compressed air energy storage plants utilize underground cavities as gas storage reservoirs. The thermodynamic properties of the stored compressed air are crucial to the safety, stability, and economic viability of the entire plant. During the design and planning phase of a compressed air energy storage plant, as well as during equipment selection, it is crucial to clearly predict how the thermodynamic parameters of the compressed air will change during operation.

[0004] Current methods for calculating compressed air thermodynamic processes primarily include analytical and numerical methods. Analytical methods calculate temperature and pressure by fixing gas storage parameters (such as storage capacity and heat exchange area) and combining analytical solutions for the filling and degassing processes. Numerical methods simulate the flow of gas within the storage reservoir to predict changes in thermodynamic parameters.

[0005] However, existing technologies do not consider the impact of the moisture content of natural air or incoming gas on thermodynamic properties. Existing technologies ignore the influence of humidity in the air, resulting in deviations in calculation results. The changes in latent heat of phase change and heat transfer characteristics caused by humidity are not incorporated into the model, resulting in errors in temperature and pressure predictions; humid air may condense or freeze during the compression / expansion process, threatening equipment life and system stability; at the same time, the failure to optimize humidity-related parameters (such as the convective heat transfer coefficient) reduces the overall efficiency of the energy storage system. Summary of the Invention

[0006] In view of this, it is necessary to provide a method and system for analyzing the thermodynamic characteristics of compressed humid air in a gas storage reservoir to solve the problems of low detection accuracy, background complexity and uneven cell distribution in nasopharyngeal carcinoma organoid images caused by existing target detection technologies.

[0007] 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 reservoir, comprising:

[0008] S1, establishing a thermodynamic theoretical model for compressing moist air in a gas storage; wherein the thermodynamic theoretical model includes:

[0009] Flow model: The turbulent k-ε model is used to simulate the turbulent flow characteristics of gas in the gas storage reservoir;

[0010] Non-isothermal flow model: couples the temperature changes during gas compression and expansion to characterize the flow behavior of gas under non-isothermal conditions;

[0011] Hygrothermal coupling model: Coupled air moisture transport equations to describe humidity diffusion, phase change, and latent heat exchange processes;

[0012] Real gas state equation: used to characterize the thermodynamic properties of humid air;

[0013] S2, constructing a three-dimensional gas storage reservoir numerical model including a fluid domain, a sealing layer lining domain, and a surrounding rock domain based on the thermodynamic theoretical model;

[0014] S3, assigning wet air material parameters to the fluid domain of the three-dimensional gas storage numerical model, assigning concrete and rock material parameters to the solid domain, and setting initial conditions and boundary conditions of the three-dimensional gas storage numerical model;

[0015] S4, solving the three-dimensional gas storage numerical model to obtain the spatial variation process of the temperature, pressure and humidity distribution of the compressed wet air in the gas storage.

[0016] Preferably, based on the thermodynamic theoretical model, a three-dimensional gas storage reservoir numerical model including a fluid domain, a sealing layer lining domain, and a surrounding rock domain is constructed, including:

[0017] According to the gas storage structure, the gas filling and discharging ports and the gas storage chamber are regarded as fluid domains to accommodate compressed wet air; the sealing layer and lining structure are regarded as sealing layer lining domains to perform sealing and protection functions; the surrounding rock is regarded as the surrounding rock domain to provide external support; and each domain is numbered;

[0018] The flow model is used to simulate the turbulent flow characteristics of gas in the fluid domain; the non-isothermal flow model is used to couple temperature changes to characterize the flow behavior under non-isothermal conditions; the moisture diffusion, phase change and latent heat exchange processes are described by the wet-heat coupling model; and the actual gas state equation characterizes the thermodynamic characteristics of humid air in the fluid domain.

[0019] Preferably, before solving the three-dimensional gas storage numerical model, the method further includes:

[0020] The three-dimensional gas storage numerical model is meshed, wherein the fluid domain adopts a high-precision boundary layer mesh.

[0021] Preferably, the setting of initial conditions and boundary conditions of the three-dimensional gas storage numerical model includes:

[0022] Depending on whether to distinguish the initial cushioning process, the initial pressure of the fluid domain is set to the ambient atmospheric pressure or the design pressure;

[0023] The pressure, flow rate and heat flux are dynamically adjusted according to the charging and discharging process of the power station operation, and the corresponding boundary conditions are set according to the flow model, non-isothermal flow model and air moisture transport equation.

[0024] Preferably, the flow model includes the fluid continuity equation, the turbulent kinetic energy k equation and the turbulent dissipation rate ε equation; wherein:

[0025] The expression of the fluid continuity equation is:

[0026]

[0027] Where ρ is the gas density; t is the time, ▽ is the Hamiltonian operator; u is the velocity vector matrix;

[0028] The expression of the turbulent kinetic energy k equation is:

[0029]

[0030] Where k is the turbulent kinetic energy; μ is the dynamic viscosity; μ T is the turbulent viscosity coefficient; p K Generate terms for turbulent kinetic energy;

[0031] The expression of the turbulent dissipation rate ε equation is:

[0032]

[0033] Where σ ε 、C ε1 and C ε2 are turbulence model constants.

[0034] Preferably, the expression of the non-isothermal flow model is:

[0035]

[0036] Where q is the heat conduction term, Cp is the specific heat capacity, T is the thermodynamic temperature, k is the thermal conductivity, Q is the heat source term, Q p is the energy source term caused by pressure change, Q vd is the viscous dissipation energy;

[0037] Among them, the energy source term Q caused by pressure change p The expression is:

[0038]

[0039] Viscous dissipated energy Q vd The expression is:

[0040]

[0041] The expression of boundary wall heat transfer -n·q is:

[0042]

[0043] Where, α p is the volume expansion coefficient, p A is the pressure, τ is the stress tensor, n is the wall normal vector pointing to the inside of the fluid; u τ is the friction velocity, Tw is the wall temperature, T + is the dimensionless temperature.

[0044] Preferably, the hygrothermal coupling model includes an air moisture transport equation, a moisture flow equation in turbulence, and a humid air phase change and temperature coupling equation;

[0045] Among them, the expression of the air moisture transport equation is:

[0046]

[0047] Where M V is the molar mass of the gas phase substance, c v is the gas phase concentration, g w is the amount of gaseous material transported per unit time through a unit area, including contributions from diffusion and convection; G is the rate of generation or consumption of gaseous material, and D is the diffusion coefficient;

[0048] The expression of the water flow equation in turbulent flow is:

[0049]

[0050] Where n is the unit normal vector of the interface, pointing to the gas phase region; C v，w is the vapor equilibrium concentration in the liquid water phase; C v is the actual vapor concentration in the gas phase; is the dimensionless concentration parameter;

[0051] The expression of the coupled equation of moist air phase change and temperature is:

[0052]

[0053] Where q evap is the latent heat flux of evaporation, which indicates the amount of heat carried away per unit time and per unit area during the evaporation process; L vis the latent heat of vaporization, which represents the heat required to convert unit mass of liquid into gas; q evap is the phase change mass flux; Q is the total heat flux, including the combined effect of 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 specific heat capacity of liquid water at constant pressure, g lc is the liquid water diffusion flux vector.

[0054] In a second aspect, an embodiment of the present invention provides a system for analyzing thermodynamic properties of compressed humid air in a gas storage reservoir, comprising:

[0055] Thermodynamic theory modeling module, used to establish a thermodynamic theory model of compressed wet air in a gas storage reservoir; wherein the thermodynamic theory model includes:

[0056] Flow model: The turbulent k-ε model is used to simulate the turbulent flow characteristics of gas in the gas storage reservoir;

[0057] Non-isothermal flow model: couples the temperature changes during gas compression and expansion to characterize the flow behavior of gas under non-isothermal conditions;

[0058] Hygrothermal coupling model: Coupled air moisture transport equations to describe humidity diffusion, phase change, and latent heat exchange processes;

[0059] Real gas state equation: used to characterize the thermodynamic properties of humid air;

[0060] A three-dimensional model building module is used to build a three-dimensional gas storage reservoir numerical model including a fluid domain, a sealing layer lining domain and a surrounding rock domain based on the thermodynamic theoretical model;

[0061] a parameter configuration module, configured to assign wet air material parameters to the fluid domain of the three-dimensional gas storage numerical model, assign concrete and rock material parameters to the solid domain, and set initial conditions and boundary conditions of the three-dimensional gas storage numerical model;

[0062] The model solving module is used to solve the three-dimensional gas storage numerical model to obtain the spatial variation process of the temperature, pressure and humidity distribution of the compressed wet air in the gas storage.

[0063] In a third aspect, the present invention further provides an electronic device comprising a memory and a processor, wherein:

[0064] The memory is used to store programs;

[0065] The processor is coupled to the memory and is used to execute the program stored in the memory to implement the steps in the method for analyzing the thermodynamic characteristics of compressed wet air in a gas storage reservoir as described in the embodiment of the first aspect of the present invention.

[0066] In a fourth aspect, the present invention further provides a computer-readable storage medium for storing computer-readable programs or instructions, which, when executed by a processor, can implement the steps in the method for analyzing the thermodynamic characteristics of compressed humid air in a gas storage reservoir as described in the embodiment of the first aspect of the present invention.

[0067] The method and system for analyzing the thermodynamic characteristics of compressed humid air in a gas storage provided by the present invention have the following beneficial effects compared with the prior art:

[0068] (1) Existing technologies ignore the effect of humidity on thermodynamic properties, leading to deviations in calculation results. The present invention establishes a coupled hygrothermal model, coupled with the air moisture transport equation, etc., and considers the changes in phase change latent heat and heat transfer characteristics caused by humidity. This allows for more accurate predictions of parameters such as temperature and pressure, reducing prediction errors.

[0069] (2) Wet air may condense or freeze during the compression / expansion process, threatening equipment life and system stability. This invention takes the impact of humidity into account and can predict these potential problems in advance, providing a basis for equipment selection and protective measures design, helping to ensure the safe and stable operation of compressed air energy storage power stations.

[0070] (3) The present invention can realistically and effectively simulate the spatial variation process of the thermodynamic characteristics of compressed air in the gas storage reservoir of a power station taking into account the influence of the moisture content in the air. This method can be used for the design calculation of underground gas storage reservoirs of compressed air energy storage power stations, and can also be used for the formulation of design specifications and standards for underground gas storage reservoirs of compressed air energy storage power stations. BRIEF DESCRIPTION OF THE DRAWINGS

[0071] Figure 1 Flow chart of the method for analyzing the thermodynamic characteristics of compressed wet air in a gas storage provided by the present invention;

[0072] Figure 2 A schematic structural diagram of a three-dimensional gas storage numerical model provided by the present invention;

[0073] Figure 3 A structural block diagram of the thermodynamic characteristics analysis system for compressed wet air in a gas storage provided by the present invention;

[0074] Figure 4 This is a structural block diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0075] The preferred embodiments of the present invention will be described in detail below in conjunction with the accompanying drawings, wherein the accompanying drawings constitute a part of this application and are used together with the embodiments of the present invention to illustrate the principles of the present invention, and are not used to limit the scope of the present invention.

[0076] References herein to "embodiments" mean that a particular feature, structure, or characteristic described in connection with the embodiments may be included in at least one embodiment of the present application. The appearance of this phrase in various places in the specification does not necessarily refer to the same embodiment, nor does it constitute an independent or alternative embodiment that is mutually exclusive of other embodiments. It is understood, both explicitly and implicitly, by those skilled in the art that the embodiments described herein may be combined with other embodiments.

[0077] Existing technologies for thermodynamic analysis of compressed air energy storage (CAES) suffer from the following major issues: The impact of the moisture content of natural air or incoming gas on thermodynamic properties is not considered, leading to biased calculation results. Humidity-induced changes in latent heat of phase change and heat transfer properties are not incorporated into the model, resulting in errors in temperature and pressure predictions. Wet air can condense or freeze during compression and expansion, threatening equipment life and system stability. Humidity-related parameters are not optimized, reducing the overall efficiency of the energy storage system and failing to meet the precise requirements for design, planning, and equipment selection of CAESE power plants.

[0078] In light of this, the present invention provides a method for analyzing the thermodynamic properties of compressed moist air in a gas storage facility. This method incorporates three-dimensional numerical calculations of the thermodynamic properties of moist air and, by incorporating the equation of state for actual humid gases, accurately simulates the temperature, pressure, and humidity distribution of compressed air within the facility, thereby improving the design accuracy and operational safety of energy storage systems. This method will be described and illustrated below through multiple examples.

[0079] Figure 1 The flow chart of the method for analyzing the thermodynamic characteristics of compressed wet air in the gas storage provided by the present invention is as follows: Figure 1 The method for analyzing the thermodynamic characteristics of compressed humid air in a gas storage provided by the present invention comprises at least the following steps:

[0080] Step S1, establishing a thermodynamic theoretical model of compressed humid air in a gas storage reservoir; wherein the thermodynamic theoretical model includes a flow model, a non-isothermal flow model, a hygrothermal coupling model, and an actual gas state equation;

[0081] Specifically, the flow model adopts the turbulent k-ε model to simulate the turbulent flow characteristics of gas in the gas storage reservoir. The flow model includes the fluid continuity equation, the turbulent kinetic energy k equation, and the turbulent dissipation rate ε equation. Among them:

[0082] The expression of the fluid continuity equation is:

[0083]

[0084] Where ρ is the gas density; t is the time, is the Hamiltonian operator; u is the velocity vector matrix;

[0085] The expression of the turbulent kinetic energy k equation is:

[0086]

[0087] Where k is the turbulent kinetic energy; μ is the dynamic viscosity; μ T is the turbulent viscosity coefficient; p K Generate terms for turbulent kinetic energy;

[0088] The expression of the turbulent dissipation rate ε equation is:

[0089]

[0090] Where σ ε 、C ε1 and C ε2 are turbulence model constants.

[0091] The non-isothermal flow model couples the temperature changes during gas compression and expansion to characterize the flow behavior of gas under non-isothermal conditions. The expression of the non-isothermal flow model is:

[0092]

[0093] Where q is the heat conduction term, Cp is the specific heat capacity, T is the thermodynamic temperature, k is the thermal conductivity, Q is the heat source term, Q p is the energy source term caused by pressure change, Q vd is the viscous dissipation energy;

[0094] Among them, the energy source term Q caused by pressure change p The expression is:

[0095]

[0096] Viscous dissipated energy Q vd The expression is:

[0097]

[0098] The expression of boundary wall heat transfer -n·q is:

[0099]

[0100] Where, α p is the volume expansion coefficient, p A is the pressure, τ is the stress tensor, n is the wall normal vector pointing to the inside of the fluid; u τ is the friction velocity, Tw is the wall temperature, T+ is the dimensionless temperature.

[0101] The coupled hygrothermal model couples the air moisture transport equation to describe humidity diffusion, phase change, and latent heat exchange. The coupled hygrothermal model includes the air moisture transport equation, the moisture flow equation in turbulent flow, and the coupled equations for moist air phase change and temperature.

[0102] Among them, the expression of the air moisture transport equation is:

[0103]

[0104] Where M V is the molar mass of the gas phase substance, c v is the gas phase concentration, g w is the amount of gaseous material transported per unit time through a unit area, including contributions from diffusion and convection; G is the rate of generation or consumption of gaseous material, and D is the diffusion coefficient;

[0105] The expression of the water flow equation in turbulent flow is:

[0106]

[0107] Where n is the unit normal vector of the interface, pointing to the gas phase region; C v，w is the vapor equilibrium concentration in the liquid water phase; C v is the actual vapor concentration in the gas phase; is the dimensionless concentration parameter;

[0108] The expression of the coupled equation of moist air phase change and temperature is:

[0109] q evap =-L v g evap

[0110]

[0111] Where q evap is the latent heat flux of evaporation, which indicates the amount of heat carried away per unit time and per unit area during the evaporation process; L v is the latent heat of vaporization, which represents the heat required to convert unit mass of liquid into gas; q evap is the phase change mass flux; Q is the total heat flux, including the combined effect of 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 specific heat capacity of liquid water at constant pressure, g lc is the liquid water diffusion flux vector.

[0112] The real gas state equation is used to characterize the thermodynamic properties of humid air. In this embodiment, the real gas state equation is used instead of the ideal gas equation. The real gas state equation that includes humidity (such as the Redlich-Kwong equation) can be used.

[0113] In step S1, a thermodynamic theory model consisting of multiple sub-models is established to comprehensively describe the characteristics of moist air compressed in a gas storage facility. The turbulent k-ε model is used to simulate the turbulent flow characteristics of gas within the gas storage facility to accurately grasp the complex flow state of the gas. The non-isothermal flow model couples the temperature changes during gas compression and expansion to accurately characterize the flow behavior of gas under temperature changes. The hygrothermal coupling model couples the air moisture transport equation to describe the humidity diffusion, phase change, and latent heat exchange processes, taking into account the influence of humidity on the thermodynamic properties of the gas. The actual gas equation of state is used to accurately characterize the thermodynamic properties of moist air. Combining these models allows for a comprehensive and precise thermodynamic theory system for moist air compressed in a gas storage facility.

[0114] Step S2: constructing a three-dimensional gas storage reservoir numerical model including a fluid domain, a sealing layer lining domain, and a surrounding rock domain based on the thermodynamic theoretical model.

[0115] Figure 2 This is a schematic diagram of the structure of the three-dimensional gas storage numerical model provided by the present invention, referring to Figure 2 The gas storage reservoir consists of a gas inlet and outlet port 1, a gas storage chamber 2, a sealing layer and lining structure 3, and surrounding rock 4. The gas inlet and outlet port 1 connects to the gas storage chamber 2, forming the core flow channel of the fluid domain; the sealing layer and lining structure 3 tightly wrap the gas storage chamber to prevent gas leakage and provide mechanical support; the surrounding rock 4 serves as the external geological environment.

[0116] In step S2, based on the gas storage structure, the gas filling and discharging ports and the gas storage chamber are used as fluid domains to accommodate compressed wet air; the sealing layer and lining structure are used as the sealing layer lining domain to perform sealing and protection functions; the surrounding rock is used as the surrounding rock domain to provide external support; and each domain is numbered.

[0117] Within each constructed domain, flow models are used to analyze gas turbulence within the fluid domain, understanding characteristics such as airflow velocity and energy dissipation. Non-isothermal flow models are used to account for the effects of temperature changes on gas flow during compression and expansion. Hygrothermal coupling models can analyze humidity diffusion, water vapor phase transitions, and latent heat exchange within the fluid domain. The real gas equation of state accurately describes the thermodynamic relationships between pressure, volume, and temperature of moist air within the fluid domain. This enables the numerical simulation of complex physical processes within the gas storage reservoir.

[0118] Step S3: assigning wet air material parameters to the fluid domain of the three-dimensional gas storage numerical model, assigning concrete and rock material parameters to the solid domain, and setting the initial conditions and boundary conditions of the three-dimensional gas storage numerical model.

[0119] Specifically, the fluid domain is assigned humid air material parameters. These material parameters include the convective heat transfer coefficient, specific heat capacity, molar mass, dynamic viscosity, density, thermal conductivity, and specific heat rate. These parameters are key indicators for describing the thermodynamic and transport properties of humid air. For example, the convective heat transfer coefficient affects the efficiency of heat transfer between the gas and the surrounding environment; the specific heat capacity determines the degree of temperature change when the gas absorbs or releases heat; and the dynamic viscosity affects the flow resistance of the gas. By accurately assigning these parameters, processes such as the flow, heat transfer, and humidity changes of humid air within the fluid domain of the gas storage reservoir can be more accurately simulated.

[0120] The solid domain consists of the sealing layer lining domain and the surrounding rock domain, which are assigned relevant material parameters for concrete and rock, respectively. These material parameters include elastic modulus, Poisson's ratio, density, and thermal conductivity. These parameters describe the mechanical and thermal properties of the solid domain. For example, the elastic modulus and Poisson's ratio determine the deformation characteristics of the solid structure under load, while thermal conductivity affects the heat conduction within the solid.

[0121] Furthermore, the initial and boundary conditions of the 3D gas storage numerical model are set. The initial pressure of the fluid domain is set based on whether the initial cushioning process is distinguished. If the initial cushioning process is distinguished, the initial pressure may be set to the design pressure; if not, it is set to the ambient atmospheric pressure.

[0122] In terms of boundary condition settings, this embodiment dynamically adjusts pressure, flow rate, and heat flux according to the power plant's operating charging and discharging processes. By combining theoretical models such as flow models, non-isothermal flow models, and air moisture transport equations, the model can simulate complex physical processes within the gas storage reservoir under different operating conditions, such as the high-speed inflow of gas during charging and the pressure drop and accompanying temperature changes during discharging.

[0123] Step S4: solving the three-dimensional gas storage numerical model to obtain the spatial variation process of the temperature, pressure and humidity distribution of the compressed wet air in the gas storage.

[0124] Specifically, before solving the three-dimensional gas storage numerical model, the three-dimensional gas storage numerical model is meshed, wherein the fluid domain adopts a high-precision boundary layer mesh. In the three-dimensional gas storage numerical model, the complex gas storage structure can be decomposed into many small, computable units through meshing, so that numerical methods (such as finite element method, finite volume method, etc.) can be applied to these units to solve the physical processes such as gas flow, heat transfer, humidity change, etc. in the gas storage. The fluid domain adopts a high-precision boundary layer mesh, which can more accurately capture the drastic changes in physical quantities in the boundary layer. Compared with ordinary grids, it can more accurately describe the flow behavior of the fluid near the wall, such as the development of turbulence on the wall, the characteristics of heat conduction in the boundary layer, etc.

[0125] In step S4, an appropriate solver is selected based on the actual operation of the gas storage facility. In actual operation, the gas storage facility involves complex physical processes such as gas flow, heat transfer, and humidity changes. For the three-dimensional gas storage numerical model used in this invention, a solver capable of handling multi-physics coupling problems can be selected to simultaneously solve the interrelated equations for flow, heat transfer, and humidity diffusion. The solver performs numerical calculations on the three-dimensional gas storage numerical model based on previously set initial conditions, boundary conditions, and assigned material parameters.

[0126] The calculated results are post-processed to extract parameters related to the thermodynamic properties of the compressed moist air within the gas storage, such as temperature, pressure, and humidity, from the massive amount of calculated data. Data analysis reveals the distribution of these parameters within the gas storage space and how they change over time.

[0127] In this embodiment, the analysis results can be presented in intuitive visualizations, such as contour plots and cloud maps of temperature, pressure, and humidity. These visualizations clearly demonstrate the spatial variations of various parameters within the gas storage, helping researchers intuitively understand the thermodynamic properties of compressed humid air within the gas storage, and providing strong support for the design, operation, and optimization of the gas storage.

[0128] Figure 3 This is a structural diagram of the thermodynamic characteristics analysis system for compressed wet air in a gas storage provided by the present invention, refer to Figure 3 The gas storage compressed wet air thermodynamic characteristics analysis system 300 includes:

[0129] Thermodynamic theory modeling module 301 is used to establish a thermodynamic theory model of compressed wet air in a gas storage reservoir; wherein the thermodynamic theory model includes:

[0130] Flow model: The turbulent k-ε model is used to simulate the turbulent flow characteristics of gas in the gas storage reservoir;

[0131] Non-isothermal flow model: couples the temperature changes during gas compression and expansion to characterize the flow behavior of gas under non-isothermal conditions;

[0132] Hygrothermal coupling model: Coupled air moisture transport equations to describe humidity diffusion, phase change, and latent heat exchange processes;

[0133] Real gas state equation: used to characterize the thermodynamic properties of humid air;

[0134] A three-dimensional model building module 302 is used to build a three-dimensional gas storage reservoir numerical model including a fluid domain, a sealing layer lining domain, and a surrounding rock domain based on the thermodynamic theoretical model;

[0135] A parameter configuration module 303 is used to assign wet air material parameters to the fluid domain of the three-dimensional gas storage numerical model, assign concrete and rock material parameters to the solid domain, and set initial conditions and boundary conditions of the three-dimensional gas storage numerical model;

[0136] The model solving module 304 is used to solve the three-dimensional gas storage numerical model to obtain the spatial variation process of the temperature, pressure and humidity distribution of the compressed wet air in the gas storage.

[0137] The system for analyzing the thermodynamic characteristics of compressed humid air in a gas storage reservoir provided by the present invention is used to execute the method for analyzing the thermodynamic characteristics of compressed humid air in a gas storage reservoir provided by the aforementioned embodiments. The method for analyzing the thermodynamic characteristics of compressed humid air in a gas storage reservoir has been described in detail in the aforementioned embodiments and will not be repeated in this embodiment.

[0138] The method and system for analyzing the thermodynamic characteristics of compressed humid air in a gas storage provided by the present invention have the following beneficial effects compared with the prior art:

[0139] (1) Existing technologies ignore the effect of humidity on thermodynamic properties, leading to deviations in calculation results. The present invention establishes a coupled hygrothermal model, coupled with the air moisture transport equation, etc., and considers the changes in phase change latent heat and heat transfer characteristics caused by humidity. This allows for more accurate predictions of parameters such as temperature and pressure, reducing prediction errors.

[0140] (2) Wet air may condense or freeze during the compression / expansion process, threatening equipment life and system stability. This invention takes the impact of humidity into account and can predict these potential problems in advance, providing a basis for equipment selection and protective measures design, helping to ensure the safe and stable operation of compressed air energy storage power stations.

[0141] (3) The present invention can realistically and effectively simulate the spatial variation process of the thermodynamic characteristics of compressed air in the gas storage reservoir of a power station taking into account the influence of the moisture content in the air. This method can be used for the design calculation of underground gas storage reservoirs of compressed air energy storage power stations, and can also be used for the formulation of design specifications and standards for underground gas storage reservoirs of compressed air energy storage power stations.

[0142] Figure 4 The structural block diagram of the electronic device provided by the present invention is as follows: Figure 4 As shown, the present invention further provides an electronic device, wherein the electronic device 400 can be a computing device such as a mobile terminal, a desktop computer, a notebook, a PDA, or a server. The electronic device 400 includes a processor 401 and a memory 402, wherein the memory 402 stores a program 403 for analyzing the thermodynamic properties of compressed humid air in a gas storage facility.

[0143] In some embodiments, the memory 402 may be an internal storage unit of a computer device, such as a hard disk or memory of the computer device. In other embodiments, the memory 402 may also be an external storage device of the computer device, such as a plug-in hard disk equipped on the computer device, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, a flash card (Flash Card), etc. Furthermore, the memory 402 may also include both an internal storage unit of the computer device and an external storage device. The memory 402 is used to store application software and various types of data installed on the computer device, such as program codes for installing the computer device. The memory 402 may also be used to temporarily store data that has been output or is to be output. In one embodiment, when the gas storage compressed humid air thermodynamic characteristics analysis program 403 is executed by the processor 401, the following steps are implemented:

[0144] S1, establishing a thermodynamic theoretical model for compressed moist air in a gas storage reservoir; wherein the thermodynamic theoretical model includes a flow model, a non-isothermal flow model, a moist-heat coupling model, and an actual gas state equation;

[0145] S2, constructing a three-dimensional gas storage reservoir numerical model including a fluid domain, a sealing layer lining domain, and a surrounding rock domain based on the thermodynamic theoretical model;

[0146] S3, assigning wet air material parameters to the fluid domain of the three-dimensional gas storage numerical model, assigning concrete and rock material parameters to the solid domain, and setting initial conditions and boundary conditions of the three-dimensional gas storage numerical model;

[0147] S4, solving the three-dimensional gas storage numerical model to obtain the spatial variation process of the temperature, pressure and humidity distribution of the compressed wet air in the gas storage.

[0148] In some embodiments, the processor 401 may be a central processing unit (CPU), a microprocessor, or other data processing chip, configured to execute program codes or process data stored in the memory 402, such as executing a thermodynamic properties analysis program for compressed humid air in a gas storage facility.

[0149] This embodiment further provides a computer-readable storage medium storing a program for analyzing the thermodynamic characteristics of compressed humid air in a gas storage facility. When the program is executed by a processor, the following steps are implemented:

[0150] S1, establishing a thermodynamic theoretical model for compressed moist air in a gas storage reservoir; wherein the thermodynamic theoretical model includes a flow model, a non-isothermal flow model, a moist-heat coupling model, and an actual gas state equation;

[0151] S2, constructing a three-dimensional gas storage reservoir numerical model including a fluid domain, a sealing layer lining domain, and a surrounding rock domain based on the thermodynamic theoretical model;

[0152] S3, assigning wet air material parameters to the fluid domain of the three-dimensional gas storage numerical model, assigning concrete and rock material parameters to the solid domain, and setting initial conditions and boundary conditions of the three-dimensional gas storage numerical model;

[0153] S4, solving the three-dimensional gas storage numerical model to obtain the spatial variation process of the temperature, pressure and humidity distribution of the compressed wet air in the gas storage.

[0154] The above-described embodiments merely illustrate several embodiments of the present invention, and while the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the present invention. It should be noted that a person skilled in the art would be able to make various modifications and improvements without departing from the spirit of the present invention, all of which fall within the scope of protection of the present invention. Therefore, the scope of protection of the present invention shall be determined by the appended claims.

[0155] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the various embodiments of the present invention.

Claims

1. A method for analyzing the thermodynamic characteristics of compressed moist air in a gas storage facility, characterized in that: include: S1, establishing a thermodynamic theoretical model for compressing moist air in a gas storage; wherein the thermodynamic theoretical model includes: Flow model: The turbulent k-ε model is used to simulate the turbulent flow characteristics of gas in the gas storage reservoir; Non-isothermal flow model: couples the temperature changes during gas compression and expansion to characterize the flow behavior of gas under non-isothermal conditions; Hygrothermal coupling model: Couples the air moisture transport equation to describe the humidity diffusion, phase change and latent heat exchange processes; Real gas state equation: used to characterize the thermodynamic properties of humid air; S2, constructing a three-dimensional gas storage reservoir numerical model including a fluid domain, a sealing layer lining domain, and a surrounding rock domain based on the thermodynamic theoretical model; S3, assigning wet air material parameters to the fluid domain of the three-dimensional gas storage numerical model, assigning concrete and rock material parameters to the solid domain, and setting initial conditions and boundary conditions of the three-dimensional gas storage numerical model; S4, solving the three-dimensional gas storage numerical model to obtain the spatial variation process of the temperature, pressure and humidity distribution of the compressed wet air in the gas storage.

2. The method for analyzing thermodynamic characteristics of compressed humid air in a gas storage facility according to claim 1, characterized in that: Based on the thermodynamic theoretical model, a three-dimensional gas storage reservoir numerical model including the fluid domain, the sealing layer lining domain, and the surrounding rock domain is constructed, including: According to the gas storage structure, the gas filling and discharging ports and the gas storage chamber are regarded as fluid domains to accommodate compressed wet air; the sealing layer and lining structure are regarded as sealing layer lining domains to perform sealing and protection functions; the surrounding rock is regarded as the surrounding rock domain to provide external support; and each domain is numbered; The flow model is used to simulate the turbulent flow characteristics of gas in the fluid domain; the non-isothermal flow model is used to couple temperature changes to characterize the flow behavior under non-isothermal conditions; the moisture diffusion, phase change and latent heat exchange processes are described by the wet-heat coupling model; and the actual gas state equation characterizes the thermodynamic characteristics of humid air in the fluid domain.

3. The method for analyzing thermodynamic characteristics of compressed humid air in a gas storage facility according to claim 1, characterized in that: Before solving the three-dimensional gas storage numerical model, the method further includes: The three-dimensional gas storage numerical model is meshed, wherein the fluid domain adopts a high-precision boundary layer mesh.

4. The method for analyzing thermodynamic characteristics of compressed humid air in a gas storage facility according to claim 1, characterized in that: The setting of the initial conditions and boundary conditions of the three-dimensional gas storage numerical model includes: Depending on whether to distinguish the initial cushioning process, the initial pressure of the fluid domain is set to the ambient atmospheric pressure or the design pressure; The pressure, flow rate and heat flux are dynamically adjusted according to the charging and discharging process of the power station operation, and the corresponding boundary conditions are set according to the flow model, non-isothermal flow model and air moisture transport equation.

5. The method for analyzing thermodynamic characteristics of compressed humid air in a gas storage facility according to claim 1, characterized in that: The flow model includes the fluid continuity equation, the turbulent kinetic energy k equation and the turbulent dissipation rate ε equation; wherein: The expression of the fluid continuity equation is: Where ρ is the gas density; t is the time, is the Hamiltonian operator; u is the velocity vector matrix; The expression of the turbulent kinetic energy K equation is: Where k is the turbulent kinetic energy; μ is the dynamic viscosity; μ T is the turbulent viscosity coefficient; p K Generate terms for turbulent kinetic energy; The expression of the turbulent dissipation rate ε equation is: Where σ ε 、C ε1 and C ε2 are turbulence model constants.

6. The method for analyzing thermodynamic characteristics of compressed humid air in a gas storage facility according to claim 5, characterized in that: The expression of the non-isothermal flow model is: Where q is the heat conduction term, Cp is the specific heat capacity, T is the thermodynamic temperature, k is the thermal conductivity, Q is the heat source term, Q p is the energy source term caused by pressure change, Q vd is the viscous dissipation energy; Among them, the energy source term Q caused by pressure change p The expression is: Viscous dissipated energy Q vd The expression is: The expression of boundary wall heat transfer -n·q is: Where, α p is the volume expansion coefficient, p A is the pressure, τ is the stress tensor, n is the wall normal vector pointing to the inside of the fluid; u τ is the friction velocity, Tw is the wall temperature, T + is the dimensionless temperature.

7. The method for analyzing thermodynamic characteristics of compressed humid air in a gas storage facility according to claim 6, characterized in that: The hygrothermal coupling model includes the air moisture transport equation, the moisture flow equation in turbulence, and the humid air phase change and temperature coupling equation; Among them, the expression of the air moisture transport equation is: Where M V is the molar mass of the gas phase substance, c v is the gas phase concentration, g w is the amount of gaseous material transported per unit time through a unit area, including contributions from diffusion and convection; G is the rate of generation or consumption of gaseous material, and D is the diffusion coefficient; The expression of the water flow equation in turbulent flow is: Where n is the unit normal vector of the interface, pointing to the gas phase region; C v，w is the vapor equilibrium concentration in the liquid water phase; C v is the actual vapor concentration in the gas phase; is the dimensionless concentration parameter; The expression of the coupled equation of moist air phase change and temperature is: q evap =-L v g evap Where q evap is the latent heat flux of evaporation, which indicates the amount of heat carried away per unit time and per unit area during the evaporation process; L v is the latent heat of vaporization, which represents the heat required to convert unit mass of liquid into gas; q evap is the phase change mass flux; Q is the total heat flux, including the combined effect of 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 specific heat capacity of liquid water at constant pressure, g lc is the liquid water diffusion flux vector.

8. A system for analyzing thermodynamic characteristics of compressed humid air in a gas storage facility, characterized in that: include: Thermodynamic theory modeling module, used to establish a thermodynamic theory model of compressed wet air in a gas storage reservoir; wherein the thermodynamic theory model includes: Flow model: The turbulent k-ε model is used to simulate the turbulent flow characteristics of gas in the gas storage reservoir; Non-isothermal flow model: couples the temperature changes during gas compression and expansion to characterize the flow behavior of gas under non-isothermal conditions; Hygrothermal coupling model: Couples the air moisture transport equation to describe the humidity diffusion, phase change and latent heat exchange processes; Real gas state equation: used to characterize the thermodynamic properties of humid air; A three-dimensional model building module is used to build a three-dimensional gas storage reservoir numerical model including a fluid domain, a sealing layer lining domain and a surrounding rock domain based on the thermodynamic theoretical model; a parameter configuration module, configured to assign wet air material parameters to the fluid domain of the three-dimensional gas storage numerical model, assign concrete and rock material parameters to the solid domain, and set initial conditions and boundary conditions of the three-dimensional gas storage numerical model; The model solving module is used to solve the three-dimensional gas storage numerical model to obtain the spatial variation process of the temperature, pressure and humidity distribution of the compressed wet air in the gas storage.

9. An electronic device, It is characterized in that comprising a memory and a processor, wherein, The memory is used to store programs; The processor is coupled to the memory and is used to execute the program stored in the memory to implement the steps of the method for analyzing the thermodynamic characteristics of compressed wet air in a gas storage reservoir as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that Used to store computer-readable programs or instructions, which, when executed by a processor, can implement the steps of the method for analyzing the thermodynamic characteristics of compressed wet air in a gas storage reservoir as described in any one of claims 1 to 7.

Images