Method and system for analyzing thermodynamic properties of compressed wet air in gas storage

Abstract

The present application relates to a kind of gas storage compression wet air thermodynamic property analysis method and system, establish the thermodynamic theory model containing flow model, non-isothermal flow model, wet heat coupling model and actual gas state equation;Further construct the three-dimensional gas storage numerical model containing fluid domain, sealing layer lining area and surrounding rock area;For fluid domain, wet air material parameters are assigned, solid domain is assigned concrete and rock material parameters, and initial and boundary conditions are set;Finally, the model is solved, and the spatial variation process of the temperature, pressure and humidity distribution of compressed wet air in the gas storage is obtained.The present application can effectively simulate the spatial variation process of the thermodynamic properties of compressed air in the power station gas storage considering the influence of moisture content in air, and the method can be used for the design calculation of compressed air energy storage power station underground gas storage, and also can be used for the formulation of compressed air energy storage power station underground gas storage design specification standard.

Description

Technical Field

[0001] This invention relates to the field of compressed air energy storage technology, and in particular to a method and system for analyzing the thermodynamic properties of compressed humid air in a gas storage tank. Background Technology

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

[0003] A compressed air energy storage power station consists of six parts: a compressor, an expander turbine, a generator, an air storage tank, a control system, and a thermal storage system. Large-scale compressed air energy storage power stations utilize underground cavities as air storage tanks, and the thermodynamic properties of the compressed air stored within are crucial to the safety, stability, and economic feasibility of the entire power station. Therefore, the changes in the thermodynamic parameters of the compressed air during operation must be clearly predicted during the design and planning stages of the compressed air energy storage power station, as well as in the selection of equipment.

[0004] Current methods for calculating compressed air thermodynamic processes mainly include analytical methods and numerical methods. Analytical methods calculate temperature and pressure by fixing gas storage parameters (such as storage capacity and heat exchange area) and combining the analytical solution of the filling and releasing process; numerical methods predict changes in thermodynamic parameters by simulating the gas flow state within the gas storage tank.

[0005] However, existing technologies do not consider the impact of humidity content in natural air or incoming gas on thermodynamic properties. This neglect of humidity leads to biased calculations. The latent heat of phase change and changes in 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 compression / expansion, threatening equipment lifespan and system stability. Furthermore, the lack of optimization of humidity-related parameters (such as convective heat transfer coefficients) 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 properties of compressed humid air in a gas storage facility, in order to solve the problems of low detection accuracy, background complexity and uneven cell distribution in existing target detection technologies for nasopharyngeal carcinoma organoid images.

[0007] To address the aforementioned problems, in a first aspect, embodiments of the present invention provide a method for analyzing the thermodynamic properties of compressed humid air in a gas storage tank, comprising:

[0008] S1, Establish a thermodynamic theoretical model for compressing humid air in a gas storage tank; 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 tank;

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

[0011] Humidity-heat coupling model: Coupled with the air moisture transport equation, used to describe humidity diffusion, phase change and latent heat exchange processes;

[0012] Real gas law: used to characterize the thermodynamic properties of moist air;

[0013] S2. Based on the aforementioned thermodynamic theoretical model, a three-dimensional numerical model of the gas storage facility is constructed, comprising the fluid domain, the sealing layer lining domain, and the surrounding rock domain.

[0014] S3, Assign humid 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 the initial and boundary conditions of the three-dimensional gas storage numerical model.

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

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

[0017] Based on the gas storage structure, the filling and discharging ports and the gas storage chamber are designated as the fluid domain, used to contain compressed humid air; the sealing layer and lining structure are designated as the sealing layer and lining domain, serving a sealing and protective function; the surrounding rock is designated as the surrounding rock domain, providing external support; and each domain is numbered.

[0018] The flow model simulates the turbulent flow characteristics of gas within the fluid domain; the non-isothermal flow model couples temperature changes to characterize the flow behavior under non-isothermal conditions; the humid-thermal coupling model describes the processes of humidity diffusion, phase change, and latent heat exchange; and the actual gas equation of state characterizes the thermodynamic properties of humid air within the fluid domain.

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

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

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

[0022] Depending on whether the initial cushion gas process is distinguished, the initial pressure of the fluid domain is set to either ambient atmospheric pressure or design pressure.

[0023] The pressure, flow rate, and heat flux are dynamically adjusted according to the gas charging and discharging process during power plant operation, and 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 for the fluid continuity equation is:

[0026]

[0027] In the formula, ρ is the gas density; t is time; ▽ is the Hamiltonian operator; and u is the velocity vector matrix.

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

[0029]

[0030] In the formula, k is the turbulent kinetic energy; μ is the dynamic viscosity; μ T p is the turbulent viscosity coefficient. K This is the turbulent kinetic energy generation term;

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

[0032]

[0033] In the formula, σ ε C ε1 and C ε2 All of these are constants for the turbulence model.

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

[0035]

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

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

[0038]

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

[0040]

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

[0042]

[0043] In the formula, α p p is the coefficient of volume expansion. A Let τ be the pressure, τ be the stress tensor, and n be the wall normal vector pointing into the fluid; u τ Let T be the friction velocity, Tw be the wall temperature, and T be the friction velocity. + The value is a dimensionless temperature.

[0044] Preferably, the humid-thermal coupling model includes an air moisture transport equation, a moisture flow equation in turbulent flow, and a humid air phase change and temperature coupling equation.

[0045] The equation for air moisture transport is as follows:

[0046]

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

[0048] The equation for water flow in turbulent flow is:

[0049]

[0050] In the formula, n is the unit normal vector of the interface, pointing to the gas phase region; C v，w C is the vapor equilibrium concentration in the liquid aqueous phase. v This represents the actual vapor concentration in the gas phase. This is a dimensionless concentration parameter;

[0051] The expression for the coupling equation between the phase transition and temperature of moist air is:

[0052]

[0053] In the formula, q evap L represents the latent heat flux of vaporization, which indicates the amount of heat removed per unit time and per unit area during the evaporation process. vThe latent heat of vaporization (q) represents the amount of heat required to convert a unit mass of liquid into a gas. evap C is the phase change mass flux; Q is the total heat flux, which includes the combined effects of latent heat of phase change and sensible heat transfer; p，v C is the isobaric specific heat capacity of water vapor. p，a C is the specific heat capacity of dry air at constant pressure. p，l The specific heat capacity at constant pressure of liquid water, g lc This is the diffusion flux vector of liquid water.

[0054] Secondly, embodiments of the present invention provide a system for analyzing the thermodynamic properties of compressed humid air in a gas storage tank, comprising:

[0055] A thermodynamic theoretical modeling module is used to establish a thermodynamic theoretical model of compressed humid air in a gas storage tank; wherein, the thermodynamic theoretical model includes:

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

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

[0058] Humidity-heat coupling model: Coupled with the air moisture transport equation, used to describe humidity diffusion, phase change and latent heat exchange processes;

[0059] Real gas law: used to characterize the thermodynamic properties of moist air;

[0060] The three-dimensional model building module is used to construct a three-dimensional numerical model of the gas storage tank, including the fluid domain, the sealing layer lining domain, and the surrounding rock domain, based on the thermodynamic theoretical model.

[0061] The parameter configuration module is used to assign moist 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 the 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 temperature, pressure and humidity distribution of compressed humid air in the gas storage.

[0063] Thirdly, the present invention also provides an electronic device, including a memory and a processor, wherein,

[0064] The memory is used to store programs;

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

[0066] Fourthly, the present invention also provides a computer-readable storage medium for storing a computer-readable program or instruction, which, when executed by a processor, enables the implementation of the steps in the method for analyzing the thermodynamic properties of compressed humid air in a gas storage tank as described in the first aspect embodiment of the present invention.

[0067] The method and system for analyzing the thermodynamic properties of compressed humid air in a gas storage tank provided by this invention have the following advantages compared with the prior art:

[0068] (1) Existing technologies neglect the influence of air humidity on thermodynamic properties, leading to deviations in calculation results. This invention establishes a humidity-heat coupling model, coupled with air moisture transport equations, and considers the latent heat of phase change and changes in heat transfer characteristics caused by humidity, which can more accurately predict parameters such as temperature and pressure and reduce prediction errors.

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

[0070] (3) This invention can realistically and effectively simulate the spatial change process of the thermodynamic properties of compressed air in the gas storage tank of a power station, taking into account the influence of the humidity in the air. This method can be used for the design calculation of underground gas storage tanks of compressed air energy storage power stations, and can also be used for the formulation of design specifications and standards for underground gas storage tanks of compressed air energy storage power stations. Attached Figure Description

[0071] Figure 1 Flowchart of the method for analyzing the thermodynamic properties of compressed humid air in a gas storage tank provided by the present invention;

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

[0073] Figure 3 The structural block diagram of the thermodynamic property analysis system for compressed humid air in a gas storage tank provided by the present invention is shown below.

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

[0075] Preferred embodiments of the present invention will now be described in detail with reference to the accompanying drawings, which form part of this application and are used together with the embodiments of the present invention to illustrate the principles of the present invention, but are not intended to limit the scope of the present invention.

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

[0077] Existing technologies for thermodynamic analysis of compressed air energy storage mainly suffer from the following problems: The influence of humidity content in natural air or incoming gas on thermodynamic properties is not considered, leading to biased calculation results. Changes in latent heat of phase change and heat transfer characteristics caused by humidity are not included in the model, resulting in errors in temperature and pressure predictions. Humid air may condense or freeze during compression / expansion, threatening equipment lifespan and system stability. The lack of optimization of humidity-related parameters reduces the overall efficiency of the energy storage system, failing to meet the precise requirements for the design, planning, and equipment selection of compressed air energy storage power stations.

[0078] In view of this, the present invention provides a method for analyzing the thermodynamic characteristics of compressed humid air in a gas storage facility. This method considers three-dimensional numerical calculations of the thermodynamic characteristics of humid air and, by introducing the actual equation of state for humid gases, accurately simulates the temperature, pressure, and humidity distribution of compressed air within the storage facility, thereby improving the design accuracy and operational safety of the energy storage system. The following will elaborate and describe this method through several embodiments.

[0079] Figure 1 The flowchart of the method for analyzing the thermodynamic properties of compressed humid air in a gas storage tank provided by this invention is shown below. Figure 1 The method for analyzing the thermodynamic properties of compressed humid air in a gas storage tank provided by this invention includes at least the following steps:

[0080] Step S1: Establish a thermodynamic theoretical model for compressing humid air in a gas storage tank; wherein, the thermodynamic theoretical model includes a flow model, a non-isothermal flow model, a hygrothermal coupling model, and an actual gas equation of state;

[0081] Specifically, the flow model adopts the turbulent k-ε model to simulate the turbulent flow characteristics of gas within the gas storage tank; the flow model includes the fluid continuity equation, the turbulent kinetic energy k equation, and the turbulent dissipation rate ε equation; among which:

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

[0083]

[0084] In the formula, ρ is the gas density; t is time. Here, is the Hamiltonian operator; u is the velocity vector matrix;

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

[0086]

[0087] In the formula, k is the turbulent kinetic energy; μ is the dynamic viscosity; μ T p is the turbulent viscosity coefficient. K This is the turbulent kinetic energy generation term;

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

[0089]

[0090] In the formula, σ ε C ε1 and C ε2 All of these are constants for the turbulence model.

[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 for the non-isothermal flow model is:

[0092]

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

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

[0095]

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

[0097]

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

[0099]

[0100] In the formula, α p p is the coefficient of volume expansion. A Let τ be the pressure, τ be the stress tensor, and n be the wall normal vector pointing into the fluid; u τ Let T be the friction velocity, Tw be the wall temperature, and T be the friction velocity.+ The value is a dimensionless temperature.

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

[0102] The equation for air moisture transport is as follows:

[0103]

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

[0105] The equation for water flow in turbulent flow is:

[0106]

[0107] In the formula, n is the unit normal vector of the interface, pointing to the gas phase region; C v，w C is the vapor equilibrium concentration in the liquid aqueous phase. v This represents the actual vapor concentration in the gas phase. This is a dimensionless concentration parameter;

[0108] The expression for the coupling equation between the phase transition and temperature of moist air is:

[0109] q evap =-L v g evap

[0110]

[0111] In the formula, q evap L represents the latent heat flux of vaporization, which indicates the amount of heat removed per unit time and per unit area during the evaporation process. v The latent heat of vaporization (q) represents the amount of heat required to convert a unit mass of liquid into a gas. evap C is the phase change mass flux; Q is the total heat flux, which includes the combined effects of latent heat of phase change and sensible heat transfer; p，v C is the isobaric specific heat capacity of water vapor. p，a C is the specific heat capacity of dry air at constant pressure. p，l The specific heat capacity at constant pressure of liquid water, g lc This is the diffusion flux vector of liquid water.

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

[0113] In step S1, a thermodynamic theoretical model comprising multiple sub-models is established to comprehensively describe the characteristics of compressed humid air in the gas storage tank. A turbulent k-ε model is used to simulate the turbulent flow characteristics of the gas within the storage tank, accurately capturing the complex flow state of the gas. A non-isothermal flow model couples the temperature changes during gas compression and expansion, realistically characterizing the gas flow behavior under temperature variations. A humid-thermal coupling model couples the air moisture transport equation, describing humidity diffusion, phase change, and latent heat exchange processes, considering the influence of humidity on the thermodynamic properties of the gas. The actual gas state equation is used to accurately characterize the thermodynamic properties of humid air. Integrating these models allows for a comprehensive and accurate construction of the thermodynamic theoretical system for compressed humid air in the gas storage tank.

[0114] Step S2: Based on the thermodynamic theoretical model, construct a three-dimensional numerical model of the gas storage tank that includes the fluid domain, the sealing layer lining domain, and the surrounding rock domain.

[0115] Figure 2 This is a schematic diagram of the structure of the three-dimensional gas storage numerical model provided by the present invention, with reference to... Figure 2 The gas storage facility consists of an inlet / outlet 1, a gas storage chamber 2, a sealing layer and lining structure 3, and surrounding rock 4. The inlet / outlet 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 wraps 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 filling and discharging ports and the gas storage chamber are designated as fluid domains to contain compressed humid air; the sealing layer and lining structure are designated as sealing layer lining domains to provide sealing and protection; the surrounding rock is designated as surrounding rock domains to provide external support; and each domain is numbered.

[0117] Within the constructed domains, flow models are used to analyze gas turbulence and understand characteristics such as airflow velocity and energy dissipation. Non-isothermal flow models are used to consider the impact of temperature changes on gas flow during compression and expansion. Humid-thermal coupling models can analyze humidity diffusion, water vapor phase change, and latent heat exchange within the fluid domain. Actual gas equations of state accurately describe the thermodynamic relationships of pressure, volume, and temperature of humid air within the fluid domain. This enables the numerical simulation of complex physical processes within a gas storage facility.

[0118] Step S3: Assign moist 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 the initial and boundary conditions of the three-dimensional gas storage numerical model.

[0119] Specifically, material parameters for humid air are assigned to the fluid domain. These parameters include convective heat transfer coefficient, specific heat capacity, molar mass, dynamic viscosity, density, thermal conductivity, and specific heat rate. These parameters are key indicators describing the thermodynamic and transport characteristics of humid air. For example, the convective heat transfer coefficient affects the efficiency of heat transfer between the gas and its 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, the flow, heat transfer, and humidity changes of humid air within the fluid domain of the gas storage tank can be simulated more precisely.

[0120] The solid domain comprises the sealing lining domain and the surrounding rock domain, which are assigned relevant material parameters to concrete and rock, respectively. The material parameters of concrete and rock include elastic modulus, Poisson's ratio, density, thermal conductivity, etc. These parameters are used to describe the mechanical and thermal properties of the solid domain. For example, elastic modulus and Poisson's ratio determine the deformation characteristics of the solid structure under stress, and thermal conductivity affects the heat conduction inside the solid.

[0121] Furthermore, the initial and boundary conditions of the three-dimensional gas storage numerical model are set. The initial pressure of the fluid domain is set according to whether the initial cushion gas process is distinguished. If the initial cushion gas process is distinguished, the initial pressure may be set to the design pressure; if not, it is set to the ambient atmospheric pressure.

[0122] Regarding the setting of boundary conditions, this embodiment dynamically adjusts the pressure, flow rate, and heat flux according to the gas charging and discharging process during power plant operation. 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 tank under different operating conditions, such as the flow state of high-speed gas inflow during charging, and the pressure reduction and accompanying temperature changes during discharging.

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

[0124] Specifically, before solving the numerical model of the three-dimensional gas storage facility, the model is meshed, with a high-precision boundary layer mesh used for the fluid domain. Meshing decomposes the complex gas storage structure into many small, computable elements, allowing numerical methods (such as the finite element method and the finite volume method) to be applied to these elements, thereby solving for physical processes such as gas flow, heat transfer, and humidity changes within the storage facility. The high-precision boundary layer mesh for the fluid domain more accurately captures the dramatic changes in physical quantities within the boundary layer. Compared to ordinary meshes, it more precisely describes the flow behavior of fluids near the wall, such as the development of turbulence at the wall and the characteristics of heat conduction in the boundary layer.

[0125] In step S4, a suitable solver is selected based on the actual operation of the gas storage facility. In actual operation, complex physical processes are involved within the gas storage facility, such as gas flow, heat transfer, and humidity changes. For the three-dimensional gas storage facility numerical model in this invention, a solver capable of handling multiphysics coupling problems can be selected to simultaneously solve the interrelated equations related to flow, heat transfer, and humidity diffusion. The solver will perform numerical calculations on the three-dimensional gas storage facility numerical model based on the previously set initial conditions, boundary conditions, and assigned material parameters.

[0126] The calculated results are post-processed to extract relevant parameters of the thermodynamic properties of compressed humid air within the gas storage facility, such as temperature, pressure, and humidity, from the large amount of computational data. Data analysis allows us to understand the distribution of these parameters within the gas storage space and their variation over time.

[0127] In this embodiment, the analysis results can be presented in an intuitive visualization format, such as plotting contour maps and cloud maps of temperature, pressure, and humidity. These visualizations can clearly show the spatial variation process of each parameter within the gas storage facility, helping researchers intuitively understand the thermodynamic characteristics of compressed humid air within the gas storage facility, and providing strong support for the design, operation, and optimization of the gas storage facility.

[0128] Figure 3 The structural block diagram of the compressed humid air thermodynamic property analysis system for gas storage provided by the present invention is shown below. Figure 3 The 300-unit compressed humid air thermodynamic property analysis system for gas storage includes:

[0129] Thermodynamic theory modeling module 301 is used to establish a thermodynamic theory model of compressed humid air in a gas storage tank; 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 tank;

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

[0132] Humidity-heat coupling model: Coupled with the air moisture transport equation, used to describe humidity diffusion, phase change and latent heat exchange processes;

[0133] Real gas law: used to characterize the thermodynamic properties of moist air;

[0134] The three-dimensional model building module 302 is used to construct a three-dimensional gas storage numerical model that includes a fluid domain, a sealing layer lining domain, and a surrounding rock domain based on the thermodynamic theoretical model.

[0135] The parameter configuration module 303 is used to assign moist 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 the 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 humid air in the gas storage.

[0137] The gas storage compressed humid air thermodynamic characteristic analysis system provided by the present invention is used to execute the gas storage compressed humid air thermodynamic characteristic analysis method provided in the foregoing embodiments. The gas storage compressed humid air thermodynamic characteristic analysis method has been described in detail in the above embodiments, and will not be repeated here.

[0138] The method and system for analyzing the thermodynamic properties of compressed humid air in a gas storage tank provided by this invention have the following advantages compared with the prior art:

[0139] (1) Existing technologies neglect the influence of air humidity on thermodynamic properties, leading to deviations in calculation results. This invention establishes a humidity-heat coupling model, coupled with air moisture transport equations, and considers the latent heat of phase change and changes in heat transfer characteristics caused by humidity, which can more accurately predict parameters such as temperature and pressure and reduce prediction errors.

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

[0141] (3) This invention can realistically and effectively simulate the spatial change process of the thermodynamic properties of compressed air in the gas storage tank of a power station, taking into account the influence of the humidity in the air. This method can be used for the design calculation of underground gas storage tanks of compressed air energy storage power stations, and can also be used for the formulation of design specifications and standards for underground gas storage tanks of compressed air energy storage power stations.

[0142] Figure 4 A structural block diagram of the electronic device provided by the present invention, such as Figure 4 As shown, the present invention also provides an electronic device, which can be a mobile terminal, desktop computer, laptop, handheld computer, server, or other computing device. 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 tank.

[0143] In some embodiments, memory 402 may be an internal storage unit of a computer device, such as a hard disk or memory. In other embodiments, memory 402 may be an external storage device of a computer device, such as a plug-in hard disk, smart media card (SMC), secure digital card (SD), flash card, etc. Further, memory 402 may include both internal and external storage units of the computer device. Memory 402 is used to store application software and various types of data installed on the computer device, such as program code installed on the computer device. Memory 402 may also be used to temporarily store data that has been output or will be output. In one embodiment, when the compressed humid air thermodynamic characteristic analysis program 403 of the gas storage tank is executed by processor 401, the following steps are implemented:

[0144] S1. Establish a thermodynamic theoretical model for compressing humid air in a gas storage tank; wherein, the thermodynamic theoretical model includes a flow model, a non-isothermal flow model, a hygrothermal coupling model, and an actual gas equation of state;

[0145] S2. Based on the aforementioned thermodynamic theoretical model, a three-dimensional numerical model of the gas storage facility is constructed, comprising the fluid domain, the sealing layer lining domain, and the surrounding rock domain.

[0146] S3, Assign humid 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 the initial and boundary conditions of the three-dimensional gas storage numerical model.

[0147] S4. Solve the three-dimensional numerical model of the gas storage tank to obtain the spatial variation process of the temperature, pressure, and humidity distribution of the compressed humid air inside the gas storage tank.

[0148] In some embodiments, processor 401 may be a central processing unit (CPU), microprocessor or other data processing chip, used to run program code stored in memory 402 or process data, such as executing a program to analyze the thermodynamic properties of compressed humid air in a gas storage tank.

[0149] This embodiment also provides a computer-readable storage medium storing a program for analyzing the thermodynamic properties of compressed humid air in a gas storage tank. When executed by a processor, this program performs the following steps:

[0150] S1. Establish a thermodynamic theoretical model for compressing humid air in a gas storage tank; wherein, the thermodynamic theoretical model includes a flow model, a non-isothermal flow model, a hygrothermal coupling model, and an actual gas equation of state;

[0151] S2. Based on the aforementioned thermodynamic theoretical model, a three-dimensional numerical model of the gas storage facility is constructed, comprising the fluid domain, the sealing layer lining domain, and the surrounding rock domain.

[0152] S3, Assign humid 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 the initial and boundary conditions of the three-dimensional gas storage numerical model.

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

[0154] The embodiments described above are merely illustrative of several implementations of the present invention, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of the invention. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of the present invention, and these modifications and improvements all fall within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should 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, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for analyzing the thermodynamic properties of compressed humid air in a gas storage tank, characterized in that, include: S1, Establish a thermodynamic theoretical model for compressing humid air in a gas storage tank; 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 tank; Non-isothermal flow model: Coupled with temperature changes during gas compression and expansion, it characterizes the flow behavior of gas under non-isothermal conditions; the expression of the non-isothermal flow model is: In the formula, q is the heat conduction term. Cp For specific heat capacity, T Thermodynamic temperature k For thermal conductivity, Q For heat source items, Q p The energy source term is caused by pressure changes. Q vd Viscous energy dissipation; Among them, the energy source term caused by pressure change Q p The expression is: Viscous dissipation energy Q vd The expression is: The expression for boundary wall heat transfer -n·q is: In the formula, α p The coefficient of volume expansion. p A For pressure, τ Let n be the stress tensor, and n be the wall normal vector pointing into the fluid. u τ For friction speed, Tw The wall temperature, T + Dimensionless temperature; Humidity-heat coupling model: Coupled with the air moisture transport equation, it is used to describe the processes of humidity diffusion, phase change and latent heat exchange; the humidity-heat coupling model includes the air moisture transport equation, the moisture flow equation in turbulence, and the coupling equation between the phase change and temperature of moist air; The equation for air moisture transport is as follows: In the formula, M V Where is the molar mass of the gaseous substance. c v For gas phase concentration, g w It represents the amount of gaseous matter transported per unit time through a unit area, including contributions from diffusion and convection. G The rate of formation or consumption of gaseous substances. D The diffusion coefficient is denoted as . The equation for water flow in turbulent flow is: In the formula, n is the unit normal vector of the interface, pointing to the gas phase region; C v，w C is the vapor equilibrium concentration in the liquid aqueous phase. v This represents the actual vapor concentration in the gas phase. φ w + This is a dimensionless concentration parameter; The expression for the coupling equation between the phase transition and temperature of moist air is: In the formula, q evap Evaporative latent heat flux represents the amount of heat removed per unit time and per unit area during the evaporation process; L v Latent heat of vaporization represents the amount of heat required to convert a unit mass of liquid into a gas. q evap This refers to the phase transition mass flux; Q Total heat flux, which includes the combined effects of latent heat of phase change and sensible heat transfer; C p，v The specific heat capacity of water vapor at constant pressure. C p，a The specific heat capacity at constant pressure of dry air. C p，l The specific heat capacity at constant pressure of liquid water, g lc This is the diffusion flux vector of liquid water; Real gas law: used to characterize the thermodynamic properties of moist air; S2. Based on the aforementioned thermodynamic theoretical model, a three-dimensional numerical model of the gas storage facility is constructed, comprising the fluid domain, the sealing layer lining domain, and the surrounding rock domain. S3, Assign humid 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 the initial and boundary conditions of the three-dimensional gas storage numerical model. S4. Solve the three-dimensional numerical model of the gas storage tank to obtain the spatial variation process of the temperature, pressure and humidity distribution of the compressed humid air in the gas storage tank.

2. The method for analyzing the thermodynamic properties of compressed humid air in a gas storage tank according to claim 1, characterized in that, Based on the aforementioned thermodynamic theoretical model, a three-dimensional numerical model of a gas storage facility is constructed, comprising a fluid domain, a sealing layer lining domain, and a surrounding rock domain, including: Based on the gas storage structure, the filling and discharging ports and the gas storage chamber are designated as the fluid domain, used to contain compressed humid air; the sealing layer and lining structure are designated as the sealing layer and lining domain, serving a sealing and protective function; the surrounding rock is designated as the surrounding rock domain, providing external support; and each domain is numbered. The flow model simulates the turbulent flow characteristics of gas within the fluid domain; the non-isothermal flow model couples temperature changes to characterize the flow behavior under non-isothermal conditions; the humid-thermal coupling model describes the processes of humidity diffusion, phase change, and latent heat exchange; and the actual gas equation of state characterizes the thermodynamic properties of humid air within the fluid domain.

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

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

5. The method for analyzing the thermodynamic properties of compressed humid air in a gas storage tank according to claim 1, characterized in that, The flow model includes the fluid continuity equation and turbulent kinetic energy. k Equations and turbulent dissipation rates ε Equation; where: The expression for the fluid continuity equation is: In the formula, ρ The density of the gas; t Let ∠ be time, ▽ be the Hamiltonian operator, and u be the velocity vector matrix; Turbulent kinetic energy k The equation is expressed as: In the formula, k It is turbulent kinetic energy; μ Dynamic viscosity; μ T The viscosity coefficient is the turbulent viscosity coefficient. p K This is the turbulent kinetic energy generation term; Turbulent dissipation rate ε The equation is expressed as: In the formula, σ ε , C ε1 and C ε2 All of these are constants for the turbulence model.

6. A system for analyzing the thermodynamic properties of compressed humid air in a gas storage tank, the system being used to execute the method for analyzing the thermodynamic properties of compressed humid air in a gas storage tank as described in any one of claims 1-5, characterized in that, The system includes: A thermodynamic theoretical modeling module is used to establish a thermodynamic theoretical model of compressed humid air in a gas storage tank; 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 tank; Non-isothermal flow model: Couples the temperature changes during gas compression and expansion processes to characterize the flow behavior of gas under non-isothermal conditions; Humidity-heat coupling model: Coupled with the air moisture transport equation, used to describe humidity diffusion, phase change and latent heat exchange processes; Real gas law: used to characterize the thermodynamic properties of moist air; The three-dimensional model building module is used to construct a three-dimensional numerical model of the gas storage tank, including the fluid domain, the sealing layer lining domain, and the surrounding rock domain, based on the thermodynamic theoretical model. The parameter configuration module is used to assign moist 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 the 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 temperature, pressure and humidity distribution of compressed humid air in the gas storage.

7. An electronic device, Its features are, Including memory and processor, among which, The memory is used to store programs; The processor, coupled to the memory, is used to execute the program stored in the memory to implement the steps in the method for analyzing the thermodynamic properties of compressed humid air in a gas storage tank as described in any one of claims 1 to 5.

8. 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 in the method for analyzing the thermodynamic properties of compressed humid air in a gas storage tank as described in any one of claims 1 to 5.

Images