CN114117732B - Method and device for simulating explosion after gas pipeline leakage - Google Patents
Method and device for simulating explosion after gas pipeline leakage Download PDFInfo
- Publication number
- CN114117732B CN114117732B CN202111218885.9A CN202111218885A CN114117732B CN 114117732 B CN114117732 B CN 114117732B CN 202111218885 A CN202111218885 A CN 202111218885A CN 114117732 B CN114117732 B CN 114117732B
- Authority
- CN
- China
- Prior art keywords
- gas pipeline
- leakage
- evaluation unit
- explosion
- gas
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004880 explosion Methods 0.000 title claims abstract description 132
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000011156 evaluation Methods 0.000 claims abstract description 88
- 238000009792 diffusion process Methods 0.000 claims abstract description 36
- 230000008859 change Effects 0.000 claims abstract description 19
- 238000004088 simulation Methods 0.000 claims description 27
- 238000002485 combustion reaction Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 7
- 238000004134 energy conservation Methods 0.000 claims description 3
- 238000012502 risk assessment Methods 0.000 claims description 3
- 239000007789 gas Substances 0.000 description 192
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 238000009826 distribution Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 230000006378 damage Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 239000003345 natural gas Substances 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000035939 shock Effects 0.000 description 4
- 230000009194 climbing Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 239000002817 coal dust Substances 0.000 description 1
- 239000003034 coal gas Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Examining Or Testing Airtightness (AREA)
Abstract
The invention discloses a method and a device for simulating explosion after gas pipeline leakage, wherein the method comprises the following steps: performing evaluation unit grid division on the gas pipeline in the target area, establishing a three-dimensional model of the gas pipeline in each evaluation unit, and setting observation points for the three-dimensional model of the gas pipeline in each evaluation unit; establishing a gas pipeline leakage control equation and a gas pipeline explosion control equation; constructing a gas pipeline leakage diffusion state model; setting a leakage scene, and calculating leakage parameters based on the leakage scene, the basic parameters of the gas pipeline in each evaluation unit, a gas pipeline leakage control equation, a gas pipeline explosion control equation and a gas pipeline leakage diffusion state model; and respectively simulating the explosion pressure and temperature change rules of the gas pipeline at each observation point at different ignition positions under different leakage scenes according to the divided grids of each evaluation unit, the three-dimensional model of the gas pipeline in each evaluation unit and the calculated leakage parameters.
Description
Technical Field
The invention relates to the technical field of gas pipelines, in particular to a method and a device for simulating explosion after gas pipeline leakage.
Background
The gas pipeline is used as an infrastructure with wide distribution, and plays an important role in resident life and industrial production. However, because of the flammability and explosiveness of the gas itself, there is a great risk of safe operation of the gas pipeline. Especially in recent years, urban infrastructure of China is quickened, and third party damage becomes a main factor of damage to the gas pipeline. Therefore, the research of explosion after leakage aiming at the third party damage of the urban buried gas pipeline is of great significance in guaranteeing the safe operation of the gas pipeline.
Natural gas easily reaches the explosion limit after leakage and diffusion, and can burn and explode when encountering open fire.
Research on numerical simulation of combustion and explosion of gas has been started for a long time abroad, CFD software such as EXSIM, FLACS and AUTOREAGAs capable of dynamically displaying the results of combustion and explosion of gas has been developed, and has been verified by field experiments. Numerical simulation verification is carried out on the experimental results of Ibrahim et al by Namansen et al Cambridge university, and the explosion overpressure in the simulation results is consistent with the experimental results. A, R, green and R.w. Upofford have studied the propagation process of coal dust and gas explosion, and a numerical calculation model is provided. The transition of the gas explosion state (transition from laminar flow to turbulent flow) in a column-like closed space was investigated numerically by Cotese et al. Hansen et al performed extensive verification simulations simulating 30 large-scale explosion tests to evaluate how FLACS predicts far-field pressure. It is concluded that FLACS predicts accident outcome well when the overpressure is below 500 mbar. Bleyer et al have studied the propagation velocity of flames and overpressure in a closed vessel and the effect of obstacles on the propagation of flames and overpressure in order to study the destructive effect of hydrogen generated in a reactor on the reactor after explosion, and have used FLACS software to simulate the situation of explosion of hydrogen in high-pressure closed equipment.
Under the large environment of rapid development of numerical simulation, domestic scholars have also conducted related researches. Lin Baiquan et al simulated the effect of obstacles on gas explosion propagation using Phoenics software, concluded that obstacles would accelerate flame propagation, and analyzed this acceleration because the obstacle-induced turbulent regions could accelerate flame combustion propagation. Qian Xinming and the like use FLACS software to carry out modeling analysis on natural gas explosion accidents occurring in the sea lake area in Beijing city, simulate the explosion results in the room after gas leakage, and study the influence of building doors and windows on explosion overpressure results. Feng Changgen and the like simulate explosion results of different ignition positions in a single-head roadway by adopting AutoReaGas, and summarize the influence of the ignition positions on explosion overpressure values. The experimental method is characterized in that a single-head roadway gas explosion experiment of Chongqing division of a high and new class of reference coal science research institute is implemented, a 200m multiplied by 2.7m multiplied by 2.7m single-head roadway model is constructed by adopting AutoReaGas, the explosion process of gas in a roadway is simulated, and the simulation result is basically consistent with the experimental result.
Therefore, the method has important significance for explosion simulation after gas pipeline leakage.
Disclosure of Invention
In view of the above, the invention provides a method and a device for simulating explosion after gas pipeline leakage, which are used for establishing a leakage control equation and an explosion control equation, calculating leakage parameters and explosion parameters required under different scenes according to the leakage control equation and the explosion control equation, further establishing a three-dimensional model, and simulating the change rule of explosion overpressure values and temperatures at different ignition positions under the condition of whether a wind field exists or not.
The first aspect of the invention provides a method for simulating explosion after leakage of a gas pipeline, which comprises the following steps: performing evaluation unit grid division on the gas pipeline in the target area, establishing a three-dimensional model of the gas pipeline in each evaluation unit, and setting observation points for the three-dimensional model of the gas pipeline in each evaluation unit; establishing a gas pipeline leakage control equation and a gas pipeline explosion control equation; constructing a gas pipeline leakage diffusion state model; setting a leakage scene, and calculating leakage parameters based on the leakage scene, the basic parameters of the gas pipeline in each evaluation unit, a gas pipeline leakage control equation, a gas pipeline explosion control equation and a gas pipeline leakage diffusion state model; and respectively simulating the explosion pressure and temperature change rules of the gas pipeline at each observation point at different ignition positions under different leakage scenes according to the divided grids of each evaluation unit, the three-dimensional model of the gas pipeline in each evaluation unit and the calculated leakage parameters.
A second aspect of the present invention provides a post-leak explosion simulation apparatus for a gas pipeline, the apparatus comprising:
the three-dimensional model building module is used for carrying out evaluation unit grid division on the gas pipeline in the target area, building a three-dimensional model of the gas pipeline in each evaluation unit and setting observation points for the three-dimensional model of the gas pipeline in each evaluation unit; the leakage explosion control equation building module is used for building a gas pipeline leakage control equation and a gas pipeline explosion control equation; the leakage diffusion state model building module is used for building a gas pipeline leakage diffusion state model; the parameter calculation module is used for setting a leakage scene and calculating leakage parameters based on the leakage scene, the basic parameters of the gas pipeline in each evaluation unit, a gas pipeline leakage control equation, a gas pipeline explosion control equation and a gas pipeline leakage diffusion state model; the simulation module is used for respectively simulating the explosion pressure and temperature change rules of the gas pipeline at each observation point at different ignition positions under different leakage scenes according to the divided grids of each evaluation unit, the three-dimensional model of the gas pipeline in each evaluation unit and the calculated leakage parameters.
A third aspect of the present invention provides a post-leak explosion simulation apparatus for a gas pipeline, the apparatus comprising: a memory for storing a computer program; a processor for implementing the following steps when executing the computer program:
performing evaluation unit grid division on the gas pipeline in the target area, establishing a three-dimensional model of the gas pipeline in each evaluation unit, and setting observation points for the three-dimensional model of the gas pipeline in each evaluation unit; establishing a gas pipeline leakage control equation and a gas pipeline explosion control equation; constructing a gas pipeline leakage diffusion state model; setting a leakage scene, and calculating leakage parameters based on the leakage scene, the basic parameters of the gas pipeline in each evaluation unit, a gas pipeline leakage control equation, a gas pipeline explosion control equation and a gas pipeline leakage diffusion state model; and respectively simulating the explosion pressure and temperature change rules of the gas pipeline at each observation point at different ignition positions under different leakage scenes according to the divided grids of each evaluation unit, the three-dimensional model of the gas pipeline in each evaluation unit and the calculated leakage parameters.
A fourth aspect of the present invention provides a computer readable storage medium having stored thereon a computer program which when executed by a processor performs the steps of:
performing evaluation unit grid division on the gas pipeline in the target area, establishing a three-dimensional model of the gas pipeline in each evaluation unit, and setting observation points for the three-dimensional model of the gas pipeline in each evaluation unit; establishing a gas pipeline leakage control equation and a gas pipeline explosion control equation; constructing a gas pipeline leakage diffusion state model; setting a leakage scene, and calculating leakage parameters based on the leakage scene, the basic parameters of the gas pipeline in each evaluation unit, a gas pipeline leakage control equation, a gas pipeline explosion control equation and a gas pipeline leakage diffusion state model; and respectively simulating the explosion pressure and temperature change rules of the gas pipeline at each observation point at different ignition positions under different leakage scenes according to the divided grids of each evaluation unit, the three-dimensional model of the gas pipeline in each evaluation unit and the calculated leakage parameters.
Drawings
For purposes of illustration and not limitation, the invention will now be described in accordance with its preferred embodiments, particularly with reference to the accompanying drawings, in which:
FIG. 1 is a flow chart of a method for simulating explosion after leakage of a gas pipeline according to an embodiment of the invention;
FIG. 2 is a schematic diagram of a grid of partitioned evaluation units;
FIG. 3 is a schematic diagram of a three-dimensional model of a gas pipeline;
FIG. 4 (a) is a graph showing the explosion pressure distribution at ignition point 1 after ignition for 0.025 s;
FIG. 4 (b) is a graph showing the explosion pressure distribution at ignition point 1 after ignition for 0.035 s;
FIG. 5 is a graph of the explosive pressure at each observation point in the target area over time;
FIG. 6 (a) is a graph showing the explosion temperature distribution at ignition point 1 after ignition for 0.025 s;
FIG. 6 (b) is an explosion temperature distribution diagram at ignition point 1 after ignition for 0.035 s;
FIG. 7 is a graph of the change in explosion temperature with time for various observation points within a target area;
FIG. 8 is a schematic diagram of a device for simulating explosion after leakage of a gas pipeline according to another embodiment of the present invention;
fig. 9 is a schematic structural diagram of an explosion simulation device after leakage of a gas pipeline according to another embodiment of the present invention.
Detailed Description
In order that the above-recited objects, features and advantages of the present invention will be more clearly understood, a more particular description of the invention will be rendered by reference to specific embodiments thereof which are illustrated in the appended drawings. It should be noted that, without conflict, the embodiments of the present invention and features in the embodiments may be combined with each other.
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention, and the described embodiments are merely some, rather than all, embodiments of the present invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The terminology used herein in the description of the invention is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention.
Fig. 1 is a flowchart of a method for simulating explosion after leakage of a gas pipeline according to an embodiment of the present application. Referring to fig. 1, the method for simulating explosion after gas pipeline leakage comprises the following steps:
and S100, carrying out grid division on the gas pipeline in the target area, establishing a three-dimensional model of the gas pipeline in each evaluation unit, and setting observation points for the three-dimensional model of the gas pipeline in each evaluation unit.
S101, setting the positions of the leakage points on the gas pipeline in the target area, and performing evaluation unit grid division on the gas pipeline in the target area based on the positions of the leakage points on the gas pipeline in the target area.
The embodiment mainly analyzes the influence of wind, no wind and different ignition positions on the leakage explosion result, the position of the leakage point is not changed as a variable, and the position of the leakage point is set as (64.05,69.2,4) according to the gas pipeline in the target area. The division of the grids directly influences the convergence of the simulation result and the time consumption of the calculation process, the grids are too loose, the calculation result is not accurate enough, even the calculation fails, the grids are too compact, the calculation difficulty is increased, and the calculation time is too long. The three-dimensional modeling and the grid division are coordinated with each other, so that objects are in the grid area, and the smooth operation process is ensured. In general, in order to make the calculation result more accurate, a fine grid is adopted, and the parameter change at the position far from the leakage port is smaller, the flow field is more stable, and a loose grid is generally adopted.
In the embodiment, the 'outer climbing pipe' of the gas pipeline in a certain district is selected to simulate the leakage and explosion of the gas pipeline. In order to reflect the leakage condition of the medium in the pipeline and the subsequent explosion accident more truly, the grid is finely divided around the leakage points, and sparse division is adopted at the position far away from the leakage points. The core area is respectively set as min (64,69,3.8) and max (64.5,69.5,4.3) in the direction of X, Y, Z, the grids in the area are subdivided into 0.06m, the grids are stretched along the min (0, 1) and max (160,118,20) directions at a position far from the leakage point, and the grid division result is shown in fig. 2.
S102, constructing a three-dimensional model of the gas pipeline in each evaluation unit by using FLACS.
Taking a target area as a certain cell as an example, modeling the 6# building, the 7# building, the 8# building and the 10# building in the cell as actual scene reference standards, and assuming that the annual average wind speed of the target area is 3m/s, southeast wind is prevailing. And establishing a three-dimensional model of the gas pipeline in each evaluation unit of the target area by FLACS software according to the collected data such as the outer climbing pipe, the surrounding buildings, the topography and the like of the gas pipeline in the district.
And a certain simplification process is performed when the modeling of the actual scene is performed. In order to sufficiently observe the influence range of the leakage explosion and the influence of the surrounding buildings on the wind field, the ground size was set to 160m×118m. The model comprises the following specific contents: four residential buildings, corresponding roads, greening trees, high-voltage wires, external air conditioners and the like are included in the whole observation area. The four residential buildings are 6# building, 7# building, 8# building and 10# building, the 6# building, the 7# building and the 8# building are on the same side, and the building spacing is 28m and 24m respectively; the building # 10 is positioned on the right side of the building # 7. Six layers of households are arranged in each resident building with the height of 18 meters, and a gas pipeline 'climbing outside' is arranged between the first layer and the second layer and enters the household at the floor connecting position. The three-dimensional model of the gas pipeline in the target area is shown in fig. 3.
S103, setting a plurality of different observation points at the positions of the leakage points of the three-dimensional model of the gas pipeline in each evaluation unit in the target area according to the positions of the leakage points on the gas pipeline.
According to the position (64.05,69.2,4) of the leakage point of the gas pipeline, in order to observe the gas injection condition at the left side (the right east) of the leakage point, observation points P1 (65,69.5,4), P2 (67,69.5,4)), P3 (68,69.5,4) and P4 (64.05,69.5,4) are arranged at different positions in the forward east direction of the leakage point, and the influence of the leakage gas cloud on the nearby high-voltage lines, trees and external air conditioners is observed.
S200, establishing a gas pipeline leakage control equation and a gas pipeline explosion control equation.
When the gas leakage diffusion in the gas pipeline is simulated, the leakage gas diffusion process meets a gas pipeline leakage control equation, and the gas pipeline leakage control equation comprises a mass conservation equation, a momentum conservation equation and an energy conservation equation. Wherein:
(1) Mass conservation equation:
wherein m is mass, kg; v is the gas volume, m 3 ；u j The velocity components in the x, y and z directions are represented by u, v and w respectively; beta j Cell face porosity in grid cells and direction j,%; beta v Cell face porosity in the grid cell and v direction,%.
(2) Momentum conservation equation:
wherein,
R i ＝-f i A i ρ|u i |u i (3)
wherein R is i Resistance in the i direction, N; beta i Cell face porosity in grid cells and direction i,%; f (f) i Is a dimensionless constant that depends on the type and direction of resistance; sigma (sigma) ij Is the stress tensor; mu (mu) eff The effective viscosity is Pa.S.
(3) Energy conservation (enthalpy conservation) equation:
wherein h is enthalpy and J.
In the simulation of gas leakage explosion of a gas pipeline, a gas pipeline explosion control equation is also followed, wherein the gas pipeline explosion control equation comprises a combustion equation and a combustion quality conveying equation,
(4) The combustion mass transfer equation is:
wherein m is the mass fraction of fuel; sigma (sigma) m For the turbulent Prandtl-Schmidt number, 0.7 was taken; r is R m Is the fuel burn rate.
(5) The combustion equation is:
wherein S is L Flame rate, m/s;is the fire source speed, m/s under certain conditions; l is the flame length, m.
S300, constructing a gas pipeline leakage diffusion state model.
And constructing a gas pipeline leakage diffusion state model based on the relation between the mass flow rate and the mass flow state of the natural gas in the gas pipeline, wherein the expression of the gas pipeline leakage diffusion state model is as follows:
when (when)When the gas is defined as subsonic flow, the mass flow calculation formula for leakage is:
A＝πr 2 (2-20)
wherein P is 0 Is ambient pressure, pa; p is the pressure of the medium in the pipeline, pa; k is the adiabatic index of the gas, for polyatomic gas natural gas, k=1.3; c (C) dg As the gas leakage coefficient, 1.0 was taken in this example, depending on the shape of the leakage point; m is the molar mass of the gas, and 0.016kg/mol is taken; r is a gas constant, 8.3144J/(mol.K); t is the gas temperature, and 293K is taken; a is the area of the leakage orifice, m 2 The method comprises the steps of carrying out a first treatment on the surface of the r is the leakage hole radius, m.
When (when)When the gas is defined as sonic flow, the mass flow calculation formula of leakage is:
according to the embodiment, a leakage control equation and an explosion equation are constructed, grid division and three-dimensional modeling are carried out on the gas pipeline in the target area, and proper observation points are set, so that preparation is made for simulation calculation.
S400, setting a leakage scene, and calculating leakage parameters based on the leakage scene, the basic parameters of the gas pipeline in each evaluation unit, a gas pipeline leakage control equation, a gas pipeline explosion control equation and a gas pipeline leakage diffusion state model.
The leakage scene includes wind field conditions, ignition location, leakage point diameter, and ignition time. The leakage scenario includes: scene one is windless, ignition point 1 (65.08,69.1,4.08), ignition point 2 (66.5,69.1,4.08), scene two is southeast wind 3m/s, ignition point 3 (65.08,69.1,4.08), ignition point 4 (66.5,69.1,4.08), leak point diameter is: 40mm, the ignition time was set to 180s.
Leakage parameters include leakage energy, leakage momentum, leakage energy, combustion quality, and mass flow of gas conduit leakage.
The method for calculating the mass flow of the gas pipeline leakage comprises the following steps:
acquiring the pressure and the gas temperature of a gas pipeline;
setting a leakage point position, acquiring the radius of the leakage point, and calculating the area of the leakage point;
substituting the pressure, the gas temperature and the leakage point area of the gas pipeline into a gas pipeline leakage diffusion state model, and calculating the mass flow of gas pipeline leakage.
The position of the leakage point is 64.05,69.2,4, the atmospheric temperature is 20 ℃, and the diameter of the leakage point is: 40mm, the ignition time was 180s, and the area of the leakage hole was 0.001256m 2 The leakage process was continued for 180 seconds, and the mass flow rate of the leaked gas reached 0.144kg/s.
Based on leakage scenes, according to a gas pipeline leakage control equation and a gas pipeline post-leakage explosion control equation, the leakage energy, the leakage momentum, the leakage energy and the combustion quality under different leakage scenes are calculated.
S500, simulating explosion pressure and temperature change rules of the gas pipeline at each observation point at different ignition positions under different leakage scenes according to the divided grids of each evaluation unit, the three-dimensional model of the gas pipeline in each evaluation unit and the calculated leakage parameters.
According to the method, in FLACS, according to the divided grids of each evaluation unit, the three-dimensional model of the gas pipeline in each evaluation unit and the calculated leakage parameters, explosion simulation is carried out on the leakage of the gas pipeline in the target area, and the damage degree of the explosion pressure and the explosion temperature of several scenes to human bodies and buildings is compared and analyzed.
Taking a scene one as an example, under the windless condition, after 180 seconds of gas leakage, the gas is ignited to generate explosion accidents after gas leakage. The ignition is carried out at the ignition point 1 (explosion is caused by the breakage of residential and civil electric wires), so that the explosion is caused, and the explosion result is analyzed by observing the overpressure distribution of high shock waves which are 4m away from the ground at different time after the ignition. The pressure shock wave of the gas explosion lasts for about 0.1 s; the pressure range of 0.025s after ignition is shown in fig. 4 (a), the pressure shock wave spreads around in an elliptical shape, the pressure of the explosion center is 0.5kPa, and the pressure sequentially decreases around the explosion center; as shown in fig. 4 (b), the explosion pressure shock wave range of 0.035s after ignition is widened.
Fig. 5 is a graph showing the change of the overpressure of an explosion of an observation point in a target area with time. The overpressure values of the observation points P1 to P4 change with time. As can be seen from the graph, the explosion overpressure of P1 to P4 alternately fluctuates with time, and the explosion overpressure tends to be stable after 0.6 seconds, and the maximum overpressure generated at the point P2 is 0.55kPa.
Under the condition of no wind in the scene, the explosion result is analyzed by observing the high temperature distribution of 4m from the ground at different time after ignition. The temperature change trend at the ignition point 1 is shown in fig. 6 (a) and 6 (b), the gas is ignited after leakage under the windless condition, the violent combustion at the initial explosion stage generates high temperature, and the temperature of the explosion center is instantaneously raised to 2000K. The high temperature region spreads in a long-strip shape along the leakage trend of the gas as the explosion occurrence time increases. The high temperature region reached a maximum range 0.17s after ignition.
Fig. 7 is a graph showing the change of the explosion temperature of the observation point in the target area with time. With the progress of explosion, the temperatures P1 to P4 are gradually reduced after being rapidly increased, the highest temperature reaches 2000K, and then the temperature is reduced to normal temperature, so that the glass window, the electric wires and the like can be damaged by the high temperature.
According to the method for simulating the explosion after the gas pipeline leaks, the leakage control equation and the explosion control equation are established, and the leakage parameters and the explosion parameters required in different scenes are calculated according to the leakage control equation and the explosion control equation. And a three-dimensional model is built, and the change rule of explosion overpressure values and temperatures at different ignition positions under the condition of wind fields is simulated. When the leaked gas explodes, the leaked gas is rapidly diffused under the action of a wind field, and the concentration is low, so that the explosion intensity is smaller than that of a windless condition. Under the influence of the same wind field, when explosion occurs at different ignition positions, the farther the ignition position is away from the leakage port, the smaller the explosion overpressure value and the impact wave overpressure influence range are relatively. The diffusion direction and the diffusion range of the gas leakage are basically consistent with the wind direction under the influence of the wind direction, when a leakage accident of a target area occurs, people are required to be seriously evacuated and the evacuation problem of residents in the downwind direction is considered, the personal safety of the target area is improved, and the personnel injury caused by the accident is avoided.
Fig. 8 is a block diagram of a device 500 for simulating explosion after leakage of a gas pipeline according to another embodiment of the present invention.
In this embodiment, the gas pipeline leakage explosion risk evaluation device 600 may be applied to a computer device, and the gas pipeline leakage explosion risk evaluation device 600 may include a plurality of functional modules composed of program code segments. Program code for each program segment in the gas pipeline leakage explosion risk assessment device 600 may be stored in a memory of a computer device and executed by at least one processor of the computer device to implement (see fig. 1 for details) a gas pipeline leakage explosion risk assessment function.
In the present embodiment, the gas pipeline leakage explosion risk evaluation apparatus 600 may be divided into a plurality of functional modules according to the functions performed thereby. The gas pipeline leakage explosion risk evaluation apparatus 600 may include: a three-dimensional model building module 601, a leakage explosion control equation building module 602, a leakage diffusion state model building module 603, a parameter calculating module 604 and a simulation module 605. The module referred to in the present invention refers to a series of computer program segments capable of being executed by at least one processor and of performing a fixed function, stored in a memory. In the present embodiment, the functions of the respective modules will be described in detail in the following embodiments.
The three-dimensional model building module 601 is configured to perform evaluation unit grid division on the gas pipeline in the target area, build a three-dimensional model of the gas pipeline in each evaluation unit, and set an observation point for the three-dimensional model of the gas pipeline in each evaluation unit;
the leakage explosion control equation building module 602 is configured to build a gas pipeline leakage control equation, a gas pipeline explosion control equation and a turbulence control equation;
a leakage diffusion state model construction module 603, configured to construct a gas pipeline leakage diffusion state model;
the parameter calculation module 604 is configured to set a leakage scene, and calculate leakage parameters based on the leakage scene, the basic parameters of the gas pipeline in each evaluation unit, a gas pipeline leakage control equation, a gas pipeline explosion control equation, and a gas pipeline leakage diffusion state model;
the simulation module 605 is configured to simulate the explosion pressure and temperature change rules of the gas pipeline at different ignition positions in different leakage scenes according to the divided evaluation grids, the three-dimensional model of the gas pipeline in each evaluation unit and the calculated leakage parameters.
Referring to fig. 9, fig. 9 is a schematic diagram of a post-leak explosion simulation apparatus for a gas pipeline according to another embodiment of the present invention, where the apparatus 700 may include:
a memory 701 for storing a computer program;
the processor 702 may be configured to execute the computer program stored in the memory 701 by performing the following steps:
performing evaluation unit grid division on the gas pipeline in the target area, establishing a three-dimensional model of the gas pipeline in each evaluation unit, and setting observation points for the three-dimensional model of the gas pipeline in each evaluation unit;
establishing a gas pipeline leakage control equation and a gas pipeline explosion control equation;
constructing a gas pipeline leakage diffusion state model;
setting a leakage scene, and calculating leakage parameters based on the leakage scene, the basic parameters of the gas pipeline in each evaluation unit, a gas pipeline leakage control equation, a gas pipeline explosion control equation and a gas pipeline leakage diffusion state model;
and respectively simulating the explosion pressure and temperature change rules of the gas pipeline at each observation point at different ignition positions under different leakage scenes according to the divided grids of each evaluation unit, the three-dimensional model of the gas pipeline in each evaluation unit and the calculated leakage parameters.
For the description of the apparatus provided by the present invention, please refer to the above method embodiment, and the description of the present invention is omitted herein.
Corresponding to the above method embodiments, the present invention also provides a computer readable storage medium having a computer program stored thereon, which when executed by a processor, performs the steps of:
performing evaluation unit grid division on the gas pipeline in the target area, establishing a three-dimensional model of the gas pipeline in each evaluation unit, and setting observation points for the three-dimensional model of the gas pipeline in each evaluation unit;
establishing a gas pipeline leakage control equation and a gas pipeline explosion control equation;
constructing a gas pipeline leakage diffusion state model;
setting a leakage scene, and calculating leakage parameters based on the leakage scene, the basic parameters of the gas pipeline in each evaluation unit, a gas pipeline leakage control equation, a gas pipeline explosion control equation and a gas pipeline leakage diffusion state model;
and respectively simulating the explosion pressure and temperature change rules of the gas pipeline at each observation point at different ignition positions under different leakage scenes according to the divided grids of each evaluation unit, the three-dimensional model of the gas pipeline in each evaluation unit and the calculated leakage parameters.
The computer readable storage medium may include: a U-disk, a removable hard disk, a Read-only memory (ROM), a random access memory (Random Access Memory, RAM), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
For the description of the computer-readable storage medium provided by the present invention, refer to the above method embodiments, and the disclosure is not repeated here.
In this specification, each embodiment is described in a progressive manner, and each embodiment is mainly described in a different point from other embodiments, so that the same or similar parts between the embodiments are referred to each other. The apparatus, device and computer readable storage medium of the embodiments are described more simply because they correspond to the methods of the embodiments, and the description thereof will be given with reference to the method section.
The above embodiments do not limit the scope of the present invention. It will be apparent to those skilled in the art that various modifications, combinations, sub-combinations and alternatives can occur depending upon design requirements and other factors. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the scope of the present invention.
Claims (9)
1. A method for simulating a post-leak explosion of a gas pipeline, comprising:
performing evaluation unit grid division on the gas pipeline in the target area, establishing a three-dimensional model of the gas pipeline in each evaluation unit, and setting observation points for the three-dimensional model of the gas pipeline in each evaluation unit;
establishing a gas pipeline leakage control equation and a gas pipeline explosion control equation;
constructing a gas pipeline leakage diffusion state model;
setting a leakage scene, and calculating leakage parameters based on the leakage scene, the basic parameters of the gas pipeline in each evaluation unit, a gas pipeline leakage control equation, a gas pipeline explosion control equation and a gas pipeline leakage diffusion state model;
according to the divided grids of each evaluation unit, the three-dimensional model of the gas pipeline in each evaluation unit and the calculated leakage parameters, respectively simulating the explosion pressure and temperature change rules of the gas pipeline at each observation point at different ignition positions under different leakage scenes;
wherein:
the step of performing evaluation unit grid division on the gas pipeline in the target area, establishing a three-dimensional model of the gas pipeline in each evaluation unit, and setting observation points on the three-dimensional model of the gas pipeline in each evaluation unit comprises the following steps:
setting the position of a leakage point on a gas pipeline in a target area;
performing evaluation unit grid division on the gas pipeline in the target area based on the positions of the leakage points on the gas pipeline in the target area;
constructing a three-dimensional model of the gas pipeline in each evaluation unit;
according to the positions of the leakage points on the gas pipeline, a plurality of different observation points are arranged at the positions of the set directions of the leakage points of the gas pipeline three-dimensional model in each evaluation unit.
2. The post-leak explosion simulation method for a gas pipeline according to claim 1, wherein the gas pipeline leak control equation includes a mass conservation equation, a momentum conservation equation, and an energy conservation equation.
3. The gas pipeline leakage explosion risk assessment method according to claim 1, wherein the gas pipeline explosion control equation comprises a combustion equation and a combustion quality transportation equation.
4. The post-leak explosion simulation method for a gas pipeline according to claim 1, wherein the gas pipeline leak diffusion state model is:
or alternatively, the first and second heat exchangers may be,
wherein P is 0 Is the ambient pressure; p is the pressure of the medium in the pipeline; k is the adiabatic index of the gas; c (C) dg Is the gas leakage coefficient; m is gasMass molar mass of the body; r is a gas constant; t is the gas temperature; a is the area of the leakage port; r is the leakage hole radius.
5. The post-leak explosion simulation method of a gas pipeline according to claim 1, wherein the leak parameter comprises a mass flow rate of the gas pipeline leak, and the mass flow rate of the gas pipeline leak is calculated by:
acquiring the pressure, the gas temperature and the radius of a leakage point of a gas pipeline;
calculating a leakage point area based on the leakage point radius;
substituting the pressure, the gas temperature and the leakage point area of the gas pipeline into a gas pipeline leakage diffusion state model, and calculating the mass flow of gas pipeline leakage.
6. The post-leak explosion simulation method for a gas pipeline according to claim 1, wherein the leak parameters include leak energy, leak momentum, leak energy and combustion quality; the calculation method of the leakage energy, the leakage momentum, the leakage energy and the combustion quality comprises the following steps:
based on leakage scenes, according to a gas pipeline leakage control equation and a gas pipeline post-leakage explosion control equation, the leakage energy, the leakage momentum, the leakage energy and the combustion quality under different leakage scenes are calculated.
7. A post-leak explosion simulation device for a gas pipeline, comprising:
the three-dimensional model building module is used for carrying out evaluation unit grid division on the gas pipeline in the target area, building a three-dimensional model of the gas pipeline in each evaluation unit and setting observation points for the three-dimensional model of the gas pipeline in each evaluation unit;
the leakage explosion control equation building module is used for building a gas pipeline leakage control equation and a gas pipeline explosion control equation;
the leakage diffusion state model building module is used for building a gas pipeline leakage diffusion state model;
the parameter calculation module is used for setting a leakage scene and calculating leakage parameters based on the leakage scene, the basic parameters of the gas pipeline in each evaluation unit, a gas pipeline leakage control equation, a gas pipeline explosion control equation and a gas pipeline leakage diffusion state model;
the simulation module is used for respectively simulating the explosion pressure and temperature change rules of the gas pipeline at each observation point at different ignition positions under different leakage scenes according to the divided grids of each evaluation unit, the three-dimensional model of the gas pipeline in each evaluation unit and the calculated leakage parameters;
wherein:
the three-dimensional model building module performs evaluation unit grid division on the gas pipeline in the target area, builds a three-dimensional model of the gas pipeline in each evaluation unit, and sets observation points for the three-dimensional model of the gas pipeline in each evaluation unit by the following mode:
setting the position of a leakage point on a gas pipeline in a target area;
performing evaluation unit grid division on the gas pipeline in the target area based on the positions of the leakage points on the gas pipeline in the target area;
constructing a three-dimensional model of the gas pipeline in each evaluation unit;
according to the positions of the leakage points on the gas pipeline, a plurality of different observation points are arranged at the positions of the set directions of the leakage points of the gas pipeline three-dimensional model in each evaluation unit.
8. A post-leak explosion simulation apparatus for a gas pipeline, comprising:
a memory for storing a computer program;
a processor for implementing the steps of the post-leak explosion simulation method for a gas pipeline according to any one of claims 1 to 6 when executing the computer program.
9. A computer-readable storage medium, characterized in that the computer-readable storage medium has stored thereon a computer program which, when executed by a processor, implements the steps of the post-gas pipeline leakage explosion simulation method according to any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111218885.9A CN114117732B (en) | 2021-10-20 | 2021-10-20 | Method and device for simulating explosion after gas pipeline leakage |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111218885.9A CN114117732B (en) | 2021-10-20 | 2021-10-20 | Method and device for simulating explosion after gas pipeline leakage |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114117732A CN114117732A (en) | 2022-03-01 |
| CN114117732B true CN114117732B (en) | 2024-04-05 |
Family
ID=80376255
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111218885.9A Active CN114117732B (en) | 2021-10-20 | 2021-10-20 | Method and device for simulating explosion after gas pipeline leakage |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114117732B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114595648B (en) * | 2022-05-11 | 2022-08-26 | 浙江浙能航天氢能技术有限公司 | Design method and system for inert gas explosion suppression scheme of hydrogenation station |
| CN116629152B (en) * | 2023-05-12 | 2024-05-03 | 中国石油大学(华东) | Energy storage cabin explosion venting modeling method for lithium battery thermal runaway exhaust explosion |
| CN117078489B (en) * | 2023-10-13 | 2024-01-16 | 深圳市胜思达科技有限公司 | Intelligent explosion venting monitoring and control system of energy storage device |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107590336A (en) * | 2017-09-13 | 2018-01-16 | 哈尔滨理工大学 | Method for numerical simulation of the fuel gas pipeline leakage to internal flow |
| CN110630909A (en) * | 2019-09-04 | 2019-12-31 | 常州大学 | Method for predicting action range of jet flame of gas pipeline |
| CN112070273A (en) * | 2020-08-11 | 2020-12-11 | 深圳市燃气集团股份有限公司 | Gas diffusion prediction method, intelligent terminal and storage medium |
| CN113821908A (en) * | 2021-08-19 | 2021-12-21 | 常州大学 | Numerical simulation method in explosion venting process of gas pipeline |
-
2021
- 2021-10-20 CN CN202111218885.9A patent/CN114117732B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107590336A (en) * | 2017-09-13 | 2018-01-16 | 哈尔滨理工大学 | Method for numerical simulation of the fuel gas pipeline leakage to internal flow |
| CN110630909A (en) * | 2019-09-04 | 2019-12-31 | 常州大学 | Method for predicting action range of jet flame of gas pipeline |
| CN112070273A (en) * | 2020-08-11 | 2020-12-11 | 深圳市燃气集团股份有限公司 | Gas diffusion prediction method, intelligent terminal and storage medium |
| CN113821908A (en) * | 2021-08-19 | 2021-12-21 | 常州大学 | Numerical simulation method in explosion venting process of gas pipeline |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114117732A (en) | 2022-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114117732B (en) | Method and device for simulating explosion after gas pipeline leakage | |
| Mao et al. | Simulation and analysis of hydrogen leakage and explosion behaviors in various compartments on a hydrogen fuel cell ship | |
| Zhang et al. | Explosion hazards of LPG-air mixtures in vented enclosure with obstacles | |
| Li et al. | Safety analysis of hydrogen leakage accident with a mobile hydrogen refueling station | |
| Pang et al. | Characteristics of external explosions induced by vented hydrogen deflagration | |
| Wang et al. | Numerical simulation analysis of explosion process and destructive effect by gas explosion accident in buildings | |
| Bie et al. | Simulation analysis on the risk of hydrogen releases and combustion in subsea tunnels | |
| CN107301495A (en) | A kind of appraisal procedure of organic rankine cycle system mixed working fluid leakage security | |
| Makarov et al. | Deflagrations of localised homogeneous and inhomogeneous hydrogen-air mixtures in enclosures | |
| Houf et al. | Evaluation of barrier walls for mitigation of unintended releases of hydrogen | |
| CN116629152B (en) | Energy storage cabin explosion venting modeling method for lithium battery thermal runaway exhaust explosion | |
| Jiang et al. | Multiparameter acceleration characteristics of premixed methane/air explosion in a semi-confined pipe | |
| Wang et al. | Numerical simulation of overpressure loads generated by gas explosions in utility tunnels | |
| CN115496003A (en) | Overpressure injury assessment method for leakage explosion accident of hydrogenation station | |
| Zhang et al. | Influence of geometrical shapes on unconfined vapor cloud explosion | |
| CN113988556A (en) | Method, device and equipment for evaluating explosion risk of gas pipeline leakage and computer readable storage medium | |
| Jiang et al. | Numerical investigation and analysis of indoor gas explosion: A case study of “6· 13” major gas explosion accident in Hubei Province, China | |
| Wang et al. | Simulation and risk assessment of hydrogen leakage in hydrogen production container | |
| Xiu et al. | Effects of combined obstacles on deflagration characteristics of hydrogen-air premixed gas | |
| CN113191067B (en) | Aluminum powder explosion overpressure prediction method and system in constrained space | |
| CN110879919B (en) | Sectional simulation method for poison diffusion under explosion effect | |
| CN106844820B (en) | Carbon dioxide inflow parameter determination method based on deep space exploration test | |
| CN112989717A (en) | Hydrogenation station fire accident thermal radiation assessment method and system | |
| Wang et al. | Research on natural gas leakage diffusion of urban underground pipeline and its explosion hazard | |
| Kang et al. | Study on leakage and explosion law of buried gas pipeline based on scenario construction |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |