Compressed air energy storage cavern rubber sealing layer leakage rate calculation method

Abstract

The invention provides a method for calculating the leakage rate of a rubber sealing layer of a compressed air energy storage cavern, and the method comprises the steps: firstly, arranging infrared temperature detectors at an arch crown and an inverted arch of the cavern in a staggered manner, so as to detect the temperature of the inner surface of the rubber sealing layer, and arranging optical fiber temperature sensors in a crossed manner in the longitudinal direction and the radial direction of the cavern, so as to detect the leakage rate of the rubber sealing layer; through the arrangement mode, the detection range of the infrared temperature detector can cover the inner surface of the rubber sealing layer, and the detection range of the optical fiber temperature sensor can cover the outer surface of the rubber sealing layer. Determining the leakage position of the rubber sealing layer according to the temperature change of the inner and outer surfaces of the rubber sealing layer; finally, high-pressure air is injected into the cavern, the leakage amount and the leakage rate of the rubber sealing layer can be calculated according to a heat conduction formula and temperature detection data at the leakage position, and then the operation efficiency of the air in the cavern in the compressed air energy storage process is evaluated.

Description

Technical Field

[0001] The invention belongs to the field of compressed air energy storage in energy engineering, and in particular relates to a method for calculating the leakage rate of a rubber sealing layer of a compressed air energy storage cavern. Background Art

[0002] Compressed air energy storage technology effectively solves the problems of unstable renewable energy supply and low energy density. In this technology, the selection of gas storage devices is extremely important. Artificially excavated underground gas storage caverns have become an important form of underground gas storage due to their high safety, strong pressure resistance, and flexible construction site selection. The sealing layer is an important guarantee for the good sealing performance of the cavern. At present, many scholars have proposed the use of rubber materials as the sealing layer of artificially excavated underground caverns.

[0003] Rubber materials have good sealing properties, but long-term heat-pressure can damage the rubber sealing layer. Damaged rubber materials will leak, which will affect the long-term safe operation of compressed air energy storage caverns. However, how to directly detect the leakage position and amount after the rubber sealing layer is damaged is still a problem that needs to be solved urgently.

[0004] Therefore, how to provide a method for calculating the leakage rate of the rubber sealing layer of a compressed air energy storage cavern, accurately determine the leakage position of the rubber sealing layer through infrared detectors and optical fiber temperature sensors, and calculate the leakage rate of the rubber sealing layer in combination with the material heat conduction formula, and then evaluate the operating efficiency of the air in the cavern during compressed air energy storage, is a problem that technical personnel in this field urgently need to solve. Summary of the invention

[0005] The purpose of the present invention is to provide a method for calculating the leakage rate of the rubber sealing layer of a compressed air energy storage cavern, so as to solve at least one of the above-mentioned technical problems.

[0006] To achieve the above-mentioned purpose, the first aspect of the present invention provides a method for calculating the leakage rate of a rubber sealing layer in a compressed air energy storage cavern, the calculation method comprising: staggeredly arranging infrared temperature detectors at the vault and invert of the cavern respectively, for detecting the temperature of the inner surface of the rubber sealing layer; cross-arranging optical fiber temperature sensors along the longitudinal direction and radial direction of the cavern respectively, for detecting the temperature of the outer surface of the rubber sealing layer; determining the leakage position of the rubber sealing layer by the temperature change of the inner surface of the rubber sealing layer or the temperature change of the outer surface of the rubber sealing layer; injecting high-pressure air into the cavern, and calculating the leakage amount and leakage rate of the rubber sealing layer according to the heat conduction formula and the temperature detection data at the leakage position.

[0007] In the first aspect, the infrared temperature detectors are staggered at the vault and invert of the cavern, respectively, and specifically include: the cavern includes, from the inside to the outside, a rubber sealing layer, a concrete lining and surrounding rock, and a plurality of infrared temperature detectors are arranged on the surface of the rubber sealing layer, and the plurality of infrared temperature detectors are staggered at the vault and invert of the cavern, and the detection range of the plurality of infrared temperature detectors covers the inner surface of the rubber sealing layer.

[0008] In the first aspect, the cross-arranging of optical fiber temperature sensors along the longitudinal direction and radial direction of the cavern specifically includes: arranging a plurality of optical fiber temperature sensors on the surface of the rubber sealing layer, and each of the optical fiber temperature sensors is located between the rubber sealing layer and the concrete lining; the plurality of optical fiber temperature sensors are arranged at intervals along the longitudinal direction and radial direction of the cavern, and the detection range of the plurality of optical fiber temperature sensors covers the outer surface of the rubber sealing layer.

[0009] In the first aspect, the heat conduction equation includes:

[0010] (1)

[0011] Where T is the temperature, k is the thermal conductivity, q is the thermal conductivity of the internal heat source, ρ is the material density, and C p is the specific heat of the material, t is the time, and x, y, and z are the three directions of the spatial coordinate system;

[0012] When used to calculate the heat conduction of the complete rubber sealing layer, there is no internal heat source in the rubber sealing layer and it is in a thermally stable state, q = 0, Then the heat conduction formula (1) is simplified to formula (2), which is as follows:

[0013] (2)

[0014] When used to calculate the heat conduction of the damaged rubber sealing layer, there is no heat conduction in the x-axis and y-axis directions, and the formula (2) is further simplified to the formula (3), which is as follows:

[0015] (3)

[0016] In the first aspect, the heat conduction equation further includes: (4) E in +E g -E out =E st ; Among them, E in is the energy entering the system, E g is the internal energy source, E outis the energy flowing out of the system, E st is the energy stored in the system.

[0017] In the first aspect, an envelope line is randomly selected in the region of the leakage position, and formula (5) is obtained according to formula (4), and formula (5) is as follows:

[0018] (5)Q 1 +Q le -Q 2 -Q ed =0;

[0019] Among them, Q 1 Q is the heat transferred from the hot air in the cave to the rubber sealing layer. 2 is the heat transferred from the rubber sealing layer to other structures in the cavern, Q le is the heat transferred from the high-temperature leaked gas to the rubber sealing layer, Q ed is the amount of heat transferred through the envelope.

[0020] In the first aspect, under thermal stability conditions, the amount of heat entering the rubber sealing layer is equal to the amount of heat flowing out of the rubber sealing layer, that is: q 1 -q 1 =0; where q 1 is the heat flux density entering the inner surface of the rubber sealing layer, q 2 is the heat flux density flowing out from the outer surface of the rubber sealing layer;

[0021] According to Fourier's law:

[0022] (6)

[0023] Among them, T 0 is the hot air temperature in the cave, T 0* is the temperature of the inner surface of the rubber sealing layer, T 1* is the temperature of the outer surface of the rubber sealing layer, T 1 is the temperature of the concrete lining, d 1 is the thickness of the inner surface of the rubber sealing layer, d 2 is the thickness of the outer surface of the rubber sealing layer;

[0024] The heat flux calculation formula in the leakage location area is:

[0025] (7)

[0026] Among them, q 1* is the heat flux density entering the inner surface of the rubber sealing layer in the leakage area, q 2* is the heat flux density flowing out from the outer surface of the rubber sealing layer in the leakage location area, T0** is the temperature of the inner surface of the rubber sealing layer in the leakage area, T 1** is the temperature of the outer surface of the rubber sealing layer in the leakage location area;

[0027] Assume that the outer surface temperature of the rubber sealing layer is the same as the inner surface temperature, that is, T 0* =T 1* , T 0** =T 1** ;

[0028] Dividing formula (6) by formula (7) yields:

[0029] (8)

[0030] By performing area integration on the selected envelope line using formula (8), we obtain:

[0031] (9)

[0032] Q ed represents the heat transferred through the envelope, that is:

[0033] (10)

[0034] Among them, k 1 is the thermal conductivity of the rubber sealing layer, dS is the integral area of ​​the leakage location, is the temperature gradient of the rubber sealing layer on the envelope, calculated by the infrared temperature detector;

[0035] Q le Indicates the heat transferred to the rubber sealing layer by the high-temperature leaking gas, that is:

[0036] (11)

[0037] Where C is the specific heat of compressed air, ρ g is the density of compressed air, m is the leakage rate, T v is the average temperature in the leak location area.

[0038] In the first aspect, the leakage rate m of the rubber sealing layer can be calculated by substituting the formulas (9)-(11) into the formula (5):

[0039]

[0040] A second aspect of the present invention provides an electronic device comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor is used to implement a method for calculating the leakage rate of a rubber sealing layer of a compressed air energy storage chamber as described in any one of the first aspects when executing the computer program stored in the memory.

[0041] A third aspect of the present invention provides a computer-readable storage medium having a computer program stored thereon: when the computer program is executed by a processor, the method for calculating the leakage rate of the rubber sealing layer of a compressed air energy storage chamber as described in any one of the first aspects is implemented.

[0042] Beneficial effects:

[0043] The present invention provides a method for calculating the leakage rate of a rubber sealing layer in a compressed air energy storage cavern. Firstly, infrared temperature detectors are staggeredly arranged at the vault and the invert of the cavern to detect the temperature of the inner surface of the rubber sealing layer, and optical fiber temperature sensors are cross-arranged along the longitudinal direction and radial direction of the cavern to detect the temperature of the outer surface of the rubber sealing layer. Through this arrangement, the detection range of the infrared temperature detector can cover the inner surface of the rubber sealing layer, and the detection range of the optical fiber temperature sensor can cover the outer surface of the rubber sealing layer; then, according to the temperature anomaly detected by the infrared detector or the optical fiber temperature sensor or the change of the inner and outer surface temperature of the rubber sealing layer, the leakage position of the rubber sealing layer is determined; finally, high-pressure air is injected into the cavern, and the leakage amount and leakage rate of the rubber sealing layer can be calculated according to the heat conduction formula and the temperature detection data at the leakage position, so as to evaluate the operation efficiency of the air in the cavern during the compressed air energy storage process. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] In order to more clearly illustrate the embodiments of this specification or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0045] Figure 1 A schematic flow chart of a method for calculating the leakage rate of a rubber sealing layer of a compressed air energy storage chamber provided by the present invention;

[0046] Figure 2 A schematic diagram of the layout of the infrared temperature detector and the optical fiber temperature sensor in the cave provided by the present invention;

[0047] Figure 3 A one-dimensional heat transfer schematic diagram provided by the present invention;

[0048] Figure 4A schematic diagram of energy flow in a leaky and damaged rubber sealing layer provided by the present invention;

[0049] Figure 5 A schematic diagram of the structure of an electronic device provided by the present invention;

[0050] Figure 6 A schematic diagram of the structure of a computer-readable storage medium provided by the present invention. DETAILED DESCRIPTION

[0051] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0052] The terms "including" and "having" and any variations thereof in the specification and claims of this application and the above-mentioned drawings are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or modules is not necessarily limited to those steps or modules that are clearly listed, but may include other steps or modules that are not clearly listed or inherent to these processes, methods, products or devices. The naming or numbering of steps in this application does not mean that the steps in the method flow must be executed in the time / logical sequence indicated by the naming or numbering. The named or numbered process steps can change the execution order according to the technical purpose to be achieved, as long as the same or similar technical effects can be achieved.

[0053] The division of modules in this application is a logical division. There may be other division methods when it is implemented in actual applications. For example, multiple modules can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, and the indirect coupling or communication connection between modules can be electrical or other similar forms, which are not limited in this application. In addition, the modules or submodules described as separate components may or may not be physically separated, may or may not be physical modules, or may be distributed in multiple circuit modules, and some or all of the modules may be selected according to actual needs to achieve the purpose of the present application.

[0054] The general idea of ​​the present invention is as follows:

[0055] When air passes through narrow cracks or pores, pressure changes and velocity increases, resulting in a slight increase in temperature. Therefore, when high-pressure air leaks in a cave, a local temperature increase usually occurs. Based on this, the present invention sets temperature detectors on both the inner and outer surfaces of the rubber sealing layer in the cave, uses the temperature abnormality value detected by the temperature detector to determine the leakage position of the rubber sealing layer, and combines the heat conduction formula and Fourier's law to calculate the leakage amount and leakage rate at the leakage position.

[0056] During the compressed gas energy storage process, the temperature in the cavern will change due to the extraction and release of gas. In order to avoid the influence of the extraction and release process on the temperature field of the rubber sealing layer, when evaluating and calculating the air leakage, it is necessary to inject a certain amount of high-pressure air into the cavern and perform a period of detection to obtain the leakage location and leakage amount.

[0057] In order to further calculate the leakage amount at the leakage position, the formula was derived by combining the heat conduction formula and Fourier's law. When high-temperature gas leaks, the high-temperature gas will transfer heat to the rubber sealing layer, making the temperature of the sealing layer in the leakage position area higher than other parts. The leakage position of the sealing layer can be located by measuring the temperature field of the sealing layer, and the air leakage rate can be calculated in combination with the temperature detector.

[0058] Embodiment 1

[0059] See also Figure 1 The first embodiment provides a method for calculating the leakage rate of a rubber sealing layer in a compressed air energy storage cavern, the calculation method comprising: staggeredly arranging infrared temperature detectors at the vault and the invert of the cavern, respectively, for detecting the temperature of the inner surface of the rubber sealing layer; cross-arranging optical fiber temperature sensors along the longitudinal direction and radial direction of the cavern, respectively, for detecting the temperature of the outer surface of the rubber sealing layer; determining the leakage position of the rubber sealing layer by the temperature change of the inner surface of the rubber sealing layer or the temperature change of the outer surface of the rubber sealing layer; injecting high-pressure air into the cavern, and calculating the leakage amount and leakage rate of the rubber sealing layer according to the heat conduction formula and the temperature detection data at the leakage position.

[0060] Specifically, the present invention provides a method for calculating the leakage rate of a rubber sealing layer in a compressed air energy storage cavern. First, infrared temperature detectors are staggered at the vault and invert of the cavern to detect the temperature of the inner surface of the rubber sealing layer, and optical fiber temperature sensors are cross-arranged along the longitudinal direction and radial direction of the cavern to detect the temperature of the outer surface of the rubber sealing layer. Through this arrangement, the detection range of the infrared temperature detector can cover the inner surface of the rubber sealing layer, and the detection range of the optical fiber temperature sensor can cover the outer surface of the rubber sealing layer; then, according to the temperature anomaly detected by the infrared detector or the optical fiber temperature sensor or the change in the temperature of the inner and outer surfaces of the rubber sealing layer, the leakage position of the rubber sealing layer is determined; finally, high-pressure air is injected into the cavern, and the leakage amount and leakage rate of the rubber sealing layer can be calculated according to the heat conduction formula and the temperature detection data at the leakage position, thereby evaluating the operating efficiency of the air in the cavern during the compressed air energy storage process.

[0061] In some possible implementations, such as Figure 2 As shown, the infrared temperature detectors are staggered at the vault and the invert of the cavern respectively, and specifically include: the cavern includes: a rubber sealing layer, a concrete lining and surrounding rock from the inside to the outside, and a plurality of infrared temperature detectors are arranged on the surface of the rubber sealing layer, and the plurality of infrared temperature detectors are staggered at the vault and the invert of the cavern in sequence, and the detection range of the plurality of infrared temperature detectors covers the inner surface of the rubber sealing layer.

[0062] Those skilled in the art can understand that infrared temperature detectors can sense infrared radiation emitted by objects. When high-pressure air leaks, the air in the leakage area will generate infrared radiation slightly higher than the surrounding environment due to the increase in temperature. Therefore, a number of infrared temperature detectors are arranged on the inner wall of the cave chamber, i.e., the inner surface of the rubber sealing layer, so that the detection range covers the inner wall of the cave chamber. In order to improve the accuracy of temperature detection, a number of infrared temperature detectors are staggered in sequence at the vault and invert of the cave chamber.

[0063] In some possible implementations, such as Figure 2 As shown, the optical fiber temperature sensors are cross-arranged along the longitudinal direction and radial direction of the cavern respectively, specifically comprising: a plurality of optical fiber temperature sensors are arranged on the surface of the rubber sealing layer, and each of the optical fiber temperature sensors is located between the rubber sealing layer and the concrete lining; a plurality of the optical fiber temperature sensors are arranged at intervals along the longitudinal direction and radial direction of the cavern, and the detection range of the plurality of the optical fiber temperature sensors covers the outer surface of the rubber sealing layer.

[0064] Furthermore, a plurality of optical fiber temperature sensors are arranged between the rubber sealing layer and the concrete lining, and the plurality of optical fiber temperature sensors are cross-arranged along the longitudinal direction and radial direction of the cavern to form a grid distributed optical fiber temperature sensor to detect the temperature of the outer surface of the rubber sealing layer.

[0065] In some possible implementations, the heat conduction formula includes:

[0066] (1)

[0067] Where T is the temperature, k is the thermal conductivity, q is the thermal conductivity of the internal heat source, ρ is the material density, and C p is the specific heat of the material, t is the time, and x, y, and z are the three directions of the spatial coordinate system;

[0068] When used to calculate the heat conduction of the complete rubber sealing layer, there is no internal heat source in the rubber sealing layer and it is in a thermally stable state, q = 0, Then the heat conduction formula (1) is simplified to formula (2), which is as follows:

[0069] (2)

[0070] When used to calculate the heat conduction of the damaged rubber sealing layer, there is no heat conduction in the x-axis and y-axis directions, and the formula (2) is further simplified to the formula (3), which is as follows: (3)

[0071] Specifically, by applying the heat conduction formula (1) to the present invention, the complex three-dimensional heat transfer can be simplified to one-dimensional heat transfer, such as Figure 3 As shown, in Figure 3 Where L is the thickness of the sealing layer, T 0 is the air temperature in the cave, T 0* is the inner surface temperature of the sealing layer, which can be obtained by an infrared temperature detector, T 1* is the outer surface temperature of the sealing layer, T 1 is the concrete lining temperature, which can be obtained through the optical fiber temperature sensor. Generally, T 0 >T 0* >T 1* >T 1 .

[0072] In some possible implementations, the heat conduction formula further includes: (4) E in +E g -E out =E st ; Among them, E in is the energy entering the system, E gis the internal energy source, E out is the energy flowing out of the system, E st is the energy stored in the system.

[0073] In some possible implementations, an envelope line is randomly selected in the region of the leakage position, and formula (5) is obtained according to formula (4), and formula (5) is as follows:

[0074] (5)Q 1 +Q le -Q 2 -Q ed =0;

[0075] Among them, Q 1 Q is the heat transferred from the hot air in the cave to the rubber sealing layer. 2 is the heat transferred from the rubber sealing layer to other structures in the cavern, Q le is the heat transferred from the high-temperature leaked gas to the rubber sealing layer, Q ed is the amount of heat transferred through the envelope.

[0076] Figure 4 Figure 5 is a schematic diagram of energy flow in a leaky and damaged rubber sealing layer. During the compressed air energy storage process, the heat transfer at any envelope line can be expressed by formula (5).

[0077] In some possible implementations, under thermal stability conditions, the amount of heat entering the rubber sealing layer is equal to the amount of heat flowing out of the rubber sealing layer, that is: 1 -q 1 =0; where q 1 is the heat flux density entering the inner surface of the rubber sealing layer, q 2 is the heat flux density flowing out from the outer surface of the rubber sealing layer;

[0078] According to Fourier's law:

[0079] (6)

[0080] Among them, T 0 is the hot air temperature in the cave, T 0* is the temperature of the inner surface of the rubber sealing layer, T 1* is the temperature of the outer surface of the rubber sealing layer, T 1 is the temperature of the concrete lining, d 1 is the thickness of the inner surface of the rubber sealing layer, d 2 is the thickness of the outer surface of the rubber sealing layer;

[0081] The heat flux calculation formula in the leakage location area is:

[0082] (7)

[0083] Among them, q 1* is the heat flux density entering the inner surface of the rubber sealing layer in the leakage area, q 2* is the heat flux density flowing out from the outer surface of the rubber sealing layer in the leakage location area, T 0** is the temperature of the inner surface of the rubber sealing layer in the leakage area, T 1** is the temperature of the outer surface of the rubber sealing layer in the leakage location area;

[0084] Assume that the outer surface temperature of the rubber sealing layer is the same as the inner surface temperature, that is, T 0* =T 1* , T 0** =T 1** ;

[0085] Dividing formula (6) by formula (7) yields:

[0086] (8)

[0087] By performing area integration on the selected envelope line using formula (8), we obtain:

[0088] (9)

[0089] Q ed represents the heat transferred through the envelope, that is:

[0090] (10)

[0091] Among them, k 1 is the thermal conductivity of the rubber sealing layer, dS is the integral area of ​​the leakage location, is the temperature gradient of the rubber sealing layer on the envelope, calculated by the infrared temperature detector;

[0092] Q le Indicates the heat transferred to the rubber sealing layer by the high-temperature leaking gas, that is:

[0093] (11)Q le =∫Cρ g (T 0 -T(V))dV=Cm(T 0 -T v );

[0094] Where C is the specific heat of compressed air, ρ g is the density of compressed air, m is the leakage rate, T v is the average temperature in the leak location area.

[0095] Specifically, since the leaking gas in the rubber sealing layer leakage area will transfer heat to the rubber sealing layer, the temperature of the rubber sealing layer at the damaged area is higher than that of other areas. Therefore, the heat transfer formula (6) of the rubber sealing layer in the thermally stable state and the heat transfer formula (7) at the leakage location can be obtained. Furthermore, the thickness of the rubber sealing layer is much smaller than the size of the cavern energy storage structure, so the temperature difference between the inner and outer surfaces of the rubber sealing layer can be ignored, that is, T 0* =T 1* , T 0** =T 1** On this basis, divide formula (6) by formula (7) to obtain formula (8), and perform area integration on formula (8) to calculate Q in formula (4): 1 -Q 2 The value of T 0** (S) and T 0* It can be obtained based on the detection data of the infrared temperature detector.

[0096] In addition, in formula (4), Q le The heat transferred from the high-temperature leaked gas to the rubber sealing layer can be calculated by the amount of leaked gas and the temperature drop, and the volume of the leaked gas is what needs to be calculated in the present invention. It is set as an unknown number and the volume of the leaked gas is integrated. The temperature drop can be obtained by subtracting the final temperature from the initial temperature in the cavern, where the initial temperature is the temperature in the cavern at the beginning of the test, and the final temperature is the average temperature in the leakage location area at the end of the test.

[0097] In some possible implementations, the leakage rate m of the rubber sealing layer can be calculated by substituting the formulas (9)-(11) into the formula (5):

[0098]

[0099] In summary, the method for calculating the leakage rate of the rubber sealing layer of the compressed air energy storage cavern provided by the present invention can accurately locate the leakage position of the rubber sealing layer according to the temperature change of the detection data by staggeredly setting a number of infrared temperature detectors on the inner surface of the rubber sealing layer and cross-setting a number of optical fiber temperature sensors on the inner surface of the rubber sealing layer, and gradually derive the calculation formula of the leakage rate m of the rubber sealing layer in combination with the heat conduction formula and Fourier's law, and then evaluate the operating efficiency of the air in the cavern during the compressed air energy storage process.

[0100] Embodiment 2

[0101] See also Figure 5 , Figure 5 FIG. 1 is a schematic diagram of an embodiment of an electronic device 100 according to an embodiment of the present invention, comprising:

[0102] A memory 101, a processor 102, and a computer program 103 stored in the memory and executable on the processor, wherein the processor implements the above-mentioned method for calculating the leakage rate of the rubber sealing layer of the compressed air energy storage chamber when executing the computer program 103 stored in the memory.

[0103] For the sake of convenience, only the parts related to the embodiment of the present invention are shown. For specific technical details not disclosed, please refer to the part of the method for calculating the leakage rate of the rubber sealing layer of the compressed air energy storage cavern in the first embodiment of the present invention. The memory 101 can be used to store a computer program 103, which includes software programs, modules and data. The processor 102 executes the computer program 103 stored in the memory 101, thereby executing various functional applications and data processing of the electronic device.

[0104] Embodiment 3

[0105] The present invention also provides a computer-readable storage medium. Figure 6 , Figure 6 The present invention is a schematic diagram of an embodiment of a computer-readable storage medium, wherein the computer-readable storage medium may store a computer program, which, when executed, includes part or all of the steps of the method for calculating the leakage rate of the rubber sealing layer of the compressed air energy storage chamber recorded in the first embodiment of the method.

[0106] Technicians in the relevant field can clearly understand that, for the convenience and simplicity of description, the specific working process of the above-described device, electronic device and computer-readable storage medium can refer to the corresponding process of the method for calculating the leakage rate of the rubber sealing layer of the compressed air energy storage chamber described in the aforementioned method embodiment, and will not be repeated here.

[0107] In the several embodiments provided in the present application, it should be understood that the disclosed systems and methods can be implemented in other ways. For example, the system embodiments described above are only schematic, for example, the division of units is only a logical function division, and there may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of systems or units, which can be electrical, mechanical or other forms.

[0108] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0109] In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0110] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art or all or part of the technical solution can be embodied in the form of a software product. The computer software product is stored in a storage medium, including a number of instructions to enable a computer device (which can be a personal computer, server, or network device, etc.) to perform all or part of the steps of the method for calculating the leakage rate of the rubber sealing layer of the compressed air energy storage chamber in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk and other media that can store program codes.

[0111] As described above, the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit the same. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that the technical solutions described in the aforementioned embodiments may still be modified, or some of the technical features thereof may be replaced by equivalents. 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 embodiments of the present invention.

Claims

1. A method for calculating the leakage rate of the rubber sealing layer of a compressed air energy storage cavern, characterized in that: The calculation method includes: Infrared temperature detectors are staggeredly arranged at the vault and invert of the cavern to detect the temperature of the inner surface of the rubber sealing layer; Optical fiber temperature sensors are cross-arranged along the longitudinal direction and radial direction of the chamber to detect the temperature of the outer surface of the rubber sealing layer; Determining the leakage position of the rubber sealing layer by the temperature change of the inner surface of the rubber sealing layer or the temperature change of the outer surface of the rubber sealing layer; High-pressure air is injected into the cavern, and the leakage amount and leakage rate of the rubber sealing layer are calculated according to the heat conduction formula and the temperature detection data at the leakage position.

2. The method for calculating the leakage rate of the rubber sealing layer of the compressed air energy storage cavern according to claim 1 is characterized in that: The infrared temperature detectors are staggeredly arranged at the vault and the invert of the cavern respectively, specifically comprising: The cavern includes, from the inside to the outside, a rubber sealing layer, a concrete lining and surrounding rock. A plurality of infrared temperature detectors are arranged on the surface of the rubber sealing layer. The plurality of infrared temperature detectors are alternately arranged at the vault and the invert of the cavern, and the detection range of the plurality of infrared temperature detectors covers the inner surface of the rubber sealing layer.

3. The method for calculating the leakage rate of the rubber sealing layer of the compressed air energy storage cavern according to claim 2 is characterized in that: The optical fiber temperature sensors are arranged crosswise along the longitudinal direction and radial direction of the cavern respectively, specifically comprising: A plurality of optical fiber temperature sensors are arranged on the surface of the rubber sealing layer, and each of the optical fiber temperature sensors is located between the rubber sealing layer and the concrete lining; the plurality of optical fiber temperature sensors are arranged at intervals along the longitudinal direction and radial direction of the cavern, and the detection range of the plurality of optical fiber temperature sensors covers the outer surface of the rubber sealing layer.

4. The method for calculating the leakage rate of the rubber sealing layer of the compressed air energy storage cavern according to claim 3 is characterized in that: The heat conduction formula includes: Where T is the temperature, k is the thermal conductivity, q is the thermal conductivity of the internal heat source, ρ is the material density, and C p is the specific heat of the material, t is the time, and x, y, and z are the three directions of the spatial coordinate system; When used to calculate the heat conduction of the complete rubber sealing layer, there is no internal heat source in the rubber sealing layer and it is in a thermally stable state, q = 0, Then the heat conduction formula (1) is simplified to formula (2), which is as follows: When used to calculate the heat conduction of the damaged rubber sealing layer, there is no heat conduction in the x-axis and y-axis directions, and the formula (2) is further simplified to the formula (3), which is as follows:

5. The method for calculating the leakage rate of the rubber sealing layer of the compressed air energy storage cavern according to claim 4 is characterized in that: The heat conduction formula also includes: (4)And in +E g -AND out =And st ; Among them, E in is the energy entering the system, E g is the internal energy source, E out is the energy flowing out of the system, E st is the energy stored in the system.

6. The method for calculating the leakage rate of the rubber sealing layer of the compressed air energy storage cavern according to claim 5 is characterized in that: An envelope line is randomly selected in the region of the leakage position, and formula (5) is obtained according to formula (4), and formula (5) is as follows: (5)Q1+Q le -Q2-Q ed =0; Among them, Q1 is the heat transferred from the hot air in the cave to the rubber sealing layer, Q2 is the heat transferred from the rubber sealing layer to other structures in the cave, and Q le is the heat transferred from the high-temperature leaked gas to the rubber sealing layer, Q ed is the amount of heat transferred through the envelope.

7. The method for calculating the leakage rate of the rubber sealing layer of the compressed air energy storage cavern according to claim 6 is characterized in that: In a thermally stable state, the heat entering the rubber sealing layer is equal to the heat flowing out of the rubber sealing layer, that is, q1-q1=0; wherein q1 is the heat flux density entering the inner surface of the rubber sealing layer, and q2 is the heat flux density flowing out of the outer surface of the rubber sealing layer; According to Fourier's law: Among them, T0 is the hot air temperature in the cave, T 0* is the temperature of the inner surface of the rubber sealing layer, T 1* is the temperature of the outer surface of the rubber sealing layer, T1 is the temperature of the concrete lining, d1 is the thickness of the inner surface of the rubber sealing layer, and d2 is the thickness of the outer surface of the rubber sealing layer; The heat flux calculation formula in the leakage location area is: Among them, q 1* is the heat flux density entering the inner surface of the rubber sealing layer in the leakage area, q 2* is the heat flux density flowing out from the outer surface of the rubber sealing layer in the leakage location area, T 0** is the temperature of the inner surface of the rubber sealing layer in the leakage area, T 1** is the temperature of the outer surface of the rubber sealing layer in the leakage location area; Assume that the outer surface temperature of the rubber sealing layer is the same as the inner surface temperature, that is, T 0* =T 1* , T 0** =T 1** ; Dividing formula (6) by formula (7) yields: By performing area integration on the selected envelope line using formula (8), we obtain: Q ed represents the heat transferred through the envelope, that is: Where k1 is the thermal conductivity of the rubber sealing layer, dS is the integral area of ​​the leakage location, is the temperature gradient of the rubber sealing layer on the envelope, calculated by the infrared temperature detector; Q le Indicates the heat transferred to the rubber sealing layer by the high-temperature leaking gas, that is: (11)Q le =∫Cρ g (T0-T(V))dV=Cm(T0-T v ); Where C is the specific heat of compressed air, ρ g is the density of compressed air, m is the leakage rate, T v is the average temperature in the leak location area.

8. The method for calculating the leakage rate of the rubber sealing layer of the compressed air energy storage cavern according to claim 7 is characterized in that: Substituting the above formulas (9)-(11) into the above formula (5), the leakage rate m of the rubber sealing layer can be calculated:

9. An electronic device, comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor is used to implement the method for calculating the leakage rate of the rubber sealing layer of a compressed air energy storage chamber as described in any one of claims 1 to 8 when executing the computer program stored in the memory.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, it implements the method for calculating the leakage rate of the rubber sealing layer of the compressed air energy storage cavern according to any one of claims 1 to 8.

Images