CN114418362A - Control room anti-knock analysis method and system - Google Patents
Control room anti-knock analysis method and system Download PDFInfo
- Publication number
- CN114418362A CN114418362A CN202111675046.XA CN202111675046A CN114418362A CN 114418362 A CN114418362 A CN 114418362A CN 202111675046 A CN202111675046 A CN 202111675046A CN 114418362 A CN114418362 A CN 114418362A
- Authority
- CN
- China
- Prior art keywords
- explosion
- leakage
- calculating
- gas
- control room
- 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.)
- Pending
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 21
- 238000004880 explosion Methods 0.000 claims abstract description 106
- 239000000463 material Substances 0.000 claims abstract description 41
- 239000007788 liquid Substances 0.000 claims abstract description 37
- 239000000126 substance Substances 0.000 claims abstract description 23
- 238000000034 method Methods 0.000 claims abstract description 15
- 230000035939 shock Effects 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims description 39
- 239000002360 explosive Substances 0.000 claims description 20
- 238000003860 storage Methods 0.000 claims description 17
- 238000002485 combustion reaction Methods 0.000 claims description 10
- 238000005422 blasting Methods 0.000 claims description 9
- 239000000203 mixture Substances 0.000 claims description 9
- 239000000446 fuel Substances 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims description 4
- 230000009471 action Effects 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000009413 insulation Methods 0.000 claims description 3
- 230000002787 reinforcement Effects 0.000 abstract description 3
- 230000000903 blocking effect Effects 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 231100000331 toxic Toxicity 0.000 description 2
- 230000002588 toxic effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009435 building construction Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000383 hazardous chemical Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- 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/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/08—Construction
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- General Business, Economics & Management (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Game Theory and Decision Science (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Examining Or Testing Airtightness (AREA)
Abstract
The invention provides a control room antiknock analysis method and a system, comprising the following steps: acquiring dangerous source data in a safe distance around a control room; calculating mass flow rate of the leaked materials according to the dangerous source data; if the leakage is liquid, calculating explosion data when chemical explosion occurs according to the mass flow rate of the leakage material, and if the leakage is gas, calculating explosion data when physical explosion occurs; and calculating the positive acting time of the steam cloud explosion shock wave received by the control room according to the explosion data. According to the method and the device, the explosion result can be calculated according to the dangerous source and the dangerous scene around the laboratory, and the result analysis is obtained according to the explosion, so that the reinforcement of the periphery of the laboratory is facilitated to reduce the influence after the explosion occurs.
Description
Technical Field
The invention belongs to the technical field of hazardous chemical management, and particularly relates to a control room antiknock analysis method and system.
Background
With the continuous refinement of petrochemical engineering building standard planning, many chemical enterprises with long-term construction age refer to the plant layout designed according to building design specifications, namely, the difference exists between the building construction specification standard and the petrochemical engineering building specification standard, and the difference causes a plurality of hidden danger problems of building spacing, chemical devices, devices of control rooms and the like. Forcing these chemical enterprises to perform according to new petrochemical specifications requires that the enterprises demolish buildings, rebuild significant equipment facilities, such as storage tanks, production areas, control rooms, etc., which are unacceptable solutions for the enterprises. In order to solve the actual requirements of the enterprises, safety analysis of various laboratories under the emergency situation facing peripheral hazard sources is analyzed, so that heat radiation, explosion intensity, explosion range and the like under the conditions of actual explosion, leakage, combustion and the like need to be scientifically calculated, and an effective anti-explosion reinforcement design is made according to the analysis result.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the control room anti-explosion analysis method and system can calculate explosion risks and consequences according to the surrounding danger sources of the laboratory, and therefore protection work of the laboratory is guided.
The invention is realized by the following steps: a control room antiknock analysis method comprises the following steps:
acquiring dangerous source data in a safe distance around a control room;
calculating mass flow rate of the leaked materials according to the dangerous source data;
if the leakage is liquid, calculating explosion data when chemical explosion occurs according to the mass flow rate of the leakage material, and if the leakage is gas, calculating explosion data when physical explosion occurs;
and calculating the positive acting time of the steam cloud explosion shock wave received by the control room according to the explosion data.
Further, said calculating a mass flow rate of the leaking material from the hazard source data comprises:
if the leakage object is liquid, judging whether the leakage scene is storage tank leakage or pipeline leakage, if so, determining a storage tank leakage model, and if so, determining a pipeline leakage model;
and if the leakage object is gas, judging whether the leakage scene is sonic flow or subsonic flow, if so, determining a sonic flow model, and if so, determining a subsonic flow model.
And calculating the mass flow rate of the leaked materials according to the leakage model and the dangerous source data.
Further, the tank leakage model is:
the pipeline leakage model is as follows:
in the formula, QmIs the mass flow rate in kg/s; p1Is the liquid pressure in the storage tank, and the unit is Pa; p0Is ambient pressure in Pa; c1The leakage coefficient of the liquid in the storage tank; a is the area of the leakage hole and is given in m2；P2Is the liquid pressure in the pipeline, and the unit is Pa; ρ is the density of the leaking liquid in kg/m3(ii) a g is the acceleration of gravity; h isLThe height of the liquid above the leakage hole is m; c2Is the leakage coefficient of the liquid in the pipe.
The sonic flow model is as follows:
the subsonic model is:
in the formula, QmIs the mass flow rate in kg/s; cdIs the gas leakage coefficient; a is the area of the leakage hole in m2；P3Is the gas pressure in the container, and the unit is Pa; m is the molecular weight of the leaking gas or vapor; rgIs an ideal gas constant with the unit of J/(mol.K); t is the gas temperature in K; y is the efflux coefficient;
in the formula, P0Is ambient pressure in Pa; p3Is the gas pressure in the container, and the unit is Pa; gamma is adiabatic index, gamma is cp/cvCp is constant-pressure heat capacity, and cv is constant-capacity heat capacity.
Further, the determining whether the leakage scene is sonic flow or subsonic flow is:
if it is
Then sonic flow, otherwise subsonic flow.
Further, the calculating explosion data when a chemical explosion occurs according to the mass flow rate of the leaked material comprises:
and (3) calculating the TNT equivalent according to the mass flow rate of the leaked material, wherein the calculation formula is as follows:
WTNT＝aA1WfQf/QTNT，
in the formula, a is the ground explosion coefficient; a. the1TNT equivalent coefficient for vapor cloud; wfIs the mass of combustible gas in the vapor cloud, and the unit is kg, Wf＝Qm·tDrain device，tDrain deviceThe time from leakage occurrence to emergency plugging completion is shortened; qfThe unit is kJ/kg of combustion heat of combustible gas; qTNTIs the explosive heat of TNT in kJ/kg;
and calculating the overpressure value of the control chamber according to the TNT equivalent, wherein the calculation formula is as follows:
in the formula, PControlThe unit is 0.1MPa for controlling the overpressure value of the chamber; wTNTTNT equivalent of explosive material in kg; r is the distance of the explosion point from the control room in m.
Further, the calculating explosion data when a chemical explosion occurs according to the mass flow rate of the leaked material comprises:
determining the range of the combustible gas cloud by using a gas diffusion model, wherein the calculation formula is as follows:
Vc＝Q/(ρ×cs)，
in the formula, VcIs the volume of the combustible gas cloud in m3(ii) a Q is the amount of leakage substance in kg, Q ═ Qm·tDrain device，tDrain deviceThe time from leakage occurrence to emergency plugging completion is shortened; rho is the vapor density in kg/m3；csIs volume fraction of fuel mixed with air, and the unit is;
calculating the combustion energy of the fuel-air mixture in the area according to the range of the combustible gas cloud, wherein the calculation formula is as follows:
E＝Vc×3.5×106，
wherein E is the combustion energy of the fuel-air mixture in the explosive source, and is given in J;
calculating a dimensionless, comparable distance for a single equivalent fuel-air mixture by the formula:
wherein R is the distance between an explosion point and a control room and is in the unit of m; pBig (a)Is the local atmospheric pressure, and has the unit of 0.1 MPa;
determining the dimensionless analog maximum lateral overpressure according to the dimensionless analog distance and the intensity of the explosive source;
calculating the overpressure applied to the control chamber according to the dimensionless simulation maximum lateral overpressure, wherein the calculation formula is as follows:
Pcontrol＝ΔPsPBig (a)，
In the formula, PControlThe unit is 0.1MPa for the overpressure value to which the control chamber is subjected.
Further, the calculating explosion data when a physical explosion occurs includes:
and calculating the blasting energy by the following formula:
in the formula, EgIs the blasting energy of the gas, in kJ; pInsulation boardIs the absolute pressure of the gas in the container, in MPa; v is the volume of the container, in m3(ii) a k is the adiabatic index of the gas, namely the ratio of the constant-pressure specific heat to the constant-volume specific heat of the gas;
and (3) calculating the TNT equivalent according to the blasting energy, wherein the calculation formula is as follows:
WTNT＝Eg/QTNT，
in the formula, QTNTIs the explosive heat of TNT in kJ/kg;
and calculating the overpressure value of the control chamber according to the TNT equivalent, wherein the calculation formula is as follows:
in the formula, PControlThe unit is 0.1MPa for controlling the overpressure value of the chamber; wTNTTNT equivalent of explosive material in kg; r is the distance of the explosion point from the control room in m.
Further, a calculation formula for calculating the positive acting time of the steam cloud explosion shock wave received by the control room according to the explosion data is as follows:
in the formula, teThe positive acting time of overpressure on the control chamber is s;
the dimensionless positive pressure action time; e is total energy of gas cloud explosion, and the unit is MJ; pBig (a)Is the local atmospheric pressure, and has the unit of 0.1 MPa; c. C0Is the speed of sound propagation in air.
The invention also provides a control room antiknock analysis system, which comprises:
the danger source data acquisition module is used for acquiring danger source data in a safe distance around the control room;
the leakage material mass flow rate calculation module is used for calculating the leakage material mass flow rate according to the dangerous source data;
the explosion data calculation module is used for calculating explosion data when chemical explosion occurs according to the mass flow rate of the leaked materials when the leaked materials are liquid, and calculating the explosion data when physical explosion occurs when the leaked materials are gas;
and the time calculation module is used for calculating the positive acting time of the steam cloud explosion shock wave received by the control room according to the explosion data.
The invention has the following beneficial effects: according to the method and the device, the explosion result can be calculated according to the dangerous source and the dangerous scene around the laboratory, and the result analysis is obtained according to the explosion, so that the reinforcement of the periphery of the laboratory is facilitated to reduce the influence after the explosion occurs.
Drawings
FIG. 1 is a flow chart of a method of the present invention;
FIG. 2 is a comparison of the comparative overpressure of the TNO model Sachs in the present invention;
FIG. 3 is a graph of the relationship between dimensionless normal operating time and dimensionless distance in the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The following examples are only for illustrating the technical solutions of the present invention more clearly, and the protection scope of the present invention is not limited thereby.
As shown in fig. 1, the present invention provides a control room antiknock analysis method, including the following steps:
s1, acquiring hazard source data including physical and chemical characteristics, whether the hazard source data belong to flammable, explosive, toxic and extremely toxic articles and the like within a safe distance around a control room, simulating the process and consequence simulation of accident occurrence, wherein the cause of the accident may be leakage of a liquid storage tank, dynamic leakage superposition condition of leaked amount and leaked substances, dynamic superposition condition of leakage through a pipeline, dynamic superposition condition of leakage of a gas substance storage tank, dynamic leakage superposition condition of leaked amount and leaked substances, and data acquisition through pipeline leakage, dynamic superposition of leakage and the like.
And S2, calculating the mass flow rate of the leaked materials according to the danger source data. The method specifically comprises the following steps:
and S21, determining a leakage model according to the leakage object and the leakage scene. The size of the leakage pore can be divided into two categories, namely complete fracture and pore leakage. When the diameter of the device (facility) is less than 150mm, the case of hole leakage and the case of complete rupture are taken for less than the diameter of the device (facility), see table below:
if the leakage object is liquid, judging whether the leakage scene is storage tank leakage or pipeline leakage, if so, determining a storage tank leakage model, and if so, determining a pipeline leakage model; if the leakage object is gas, judging whether the leakage scene is sonic flow or subsonic flow, and if the leakage scene is subsonic flow
Then the flow model of the sonic speed is determined for calculation, otherwise, the flow model of the subsonic speed is determined for calculation.
And S22, calculating the mass flow rate of the leaked materials according to the leakage model and the dangerous source data.
The storage tank leakage model is as follows:
in the formula, QmIs the mass flow rate of liquid leaking from the holes in the storage tank, and the unit is kg/s; p1Is the liquid pressure in the storage tank, and the unit is Pa; p0Is ambient pressure in Pa; g is the gravity acceleration and takes 9.8m/s2；hLThe height of the liquid above the leakage hole; c1The values for the leakage coefficient of the liquid in the tank are given in the following table:
| reynolds number | Circular shape | Triangle shape | Rectangle |
| ＞100 | 0.65 | 0.6 | 0.55 |
| ≤100 | 0.5 | 0.45 | 0.4 |
The pipeline leakage model is as follows:
in the formula, QmThe mass flow rate of liquid leaking through the holes in the pipeline is kg/s; p0Is ambient pressure in Pa; c1The leakage coefficient of the liquid in the storage tank; a is the area of the leakage hole and is given in m2；P2Is the liquid pressure in the pipeline, and the unit is Pa; ρ is the density of the leaking liquid in kg/m3(ii) a The unit is m; c2For liquids in pipesThe leakage coefficient is obtained by referring to the following experience:
(1) for sharp holes and reynolds numbers greater than 30000, the liquid leakage coefficient is approximately 0.61, for which case the outflow rate of liquid is independent of the size of the slit;
(2) for a smooth nozzle, the liquid leakage coefficient may be taken to be approximately 1;
(3) for a nozzle attached to the container (i.e., a length to diameter ratio of not less than 3), the liquid leakage coefficient is approximately 0.81;
(4) when the liquid leakage coefficient is unknown or cannot be determined, taking 1.0 maximizes the calculated flow rate.
The sonic flow model is as follows:
the subsonic model is:
in the formula, QmIs the mass flow rate in kg/s; cdThe gas leakage coefficient is related to the shape of a leakage hole, wherein the shape of the leakage hole is 1.00 when the leakage hole is circular, 0.95 when the leakage hole is triangular, and 0.90 when the leakage hole is rectangular; a is the area of the leakage hole in m2；P3Is the gas pressure in the container, and the unit is Pa; m is the molecular weight of the leaking gas or vapor; rgIs an ideal gas constant with the unit of J/(mol.K); t is the gas temperature in K; and Y is the outflow coefficient.
In the formula, P0Is ambient pressure in Pa; p3Is the gas pressure in the container, and the unit is Pa; gamma is adiabatic index, gamma is cp/cvCp is constant-pressure heat capacity, and cv is constant-capacity heat capacity.
And S3, if the leakage is liquid, calculating explosion data when chemical explosion occurs according to the mass flow rate of the leakage material, and if the leakage is gas, calculating the explosion data when physical explosion occurs.
The TNT equivalent method, S31, is a simple method of equating a fuel of known energy to TNT, based on the assumption that the fuel explosion behaves like a TNT explosion with equal energy. The method comprises the following steps:
s311, calculating the TNT equivalent according to the mass flow rate of the leaked material, wherein the calculation formula is as follows:
WTNT＝aA1WfQf/QTNT，
in the formula, a is the ground explosion coefficient and is 1.8; a. the1Taking the TNT equivalent coefficient of the vapor cloud to be 0.04; wfIs the mass of combustible gas in the vapor cloud, and the unit is kg, Wf＝Qm·tDrain device，tDrain deviceGenerally taking 10min for the time from leakage occurrence to emergency plugging completion; qfThe unit is kJ/kg of combustion heat of combustible gas; qTNTThe explosive heat of TNT is expressed in kJ/kg, and 4500kJ/kg is taken.
S312, calculating the overpressure value of the control chamber according to the TNT equivalent, wherein the calculation formula is as follows:
The multi-energy method of S32, TNO (the Netherlands organization) was proposed and gradually perfected by the Dutch TNO laboratory in 1985 on the basis of a large number of experiments and numerical studies. The TNO multipotency method belongs to a one-dimensional numerical model, is improved on the basis of the TNO model, comprehensively considers various factors such as turbulence acceleration, local constraint, gas activity and the like, and is a method for analyzing the vapor cloud explosion principle. It assumes a vapor cloud with a hemispherical shape, a central ignition, and a flame propagating outward at a constant velocity. Defining different explosion intensities according to the difference of the propagation speeds of the explosion intensities, and obtaining corresponding data through calculation to obtain a correlation curve of dimensionless peak value lateral overpressure and dimensionless distance, which is called an explosion intensity curve, and the following steps of:
s321, determining the range of the combustible gas cloud according to the gas diffusion model, wherein the calculation formula is as follows:
Vc＝Q/(ρ×cs)，
in the formula, VcIs the volume of the combustible gas cloud in m3；
Q is the amount of leakage substance in kg, Q ═ Qm·tDrain device，tDrain deviceGenerally taking 10min for the time from leakage occurrence to emergency plugging completion; rho is the vapor density in kg/m3；csIs the volume fraction of the fuel mixed with air in%.
S322, calculating the combustion energy of the fuel-air mixture in the area according to the range of the combustible gas cloud, wherein the calculation formula is as follows:
E＝Vc×3.5×106，
where E is the combustion energy of the fuel-air mixture in the explosive source and is given in J.
S323, calculating a dimensionless analog distance of the single equivalent fuel-air mixture, wherein the calculation formula is as follows:
wherein R is the distance between an explosion point and a control room and is in the unit of m;
Pbig (a)Is the local atmospheric pressure and has the unit of 0.1 MPa.
And S324, determining the dimensionless comparative maximum lateral overpressure according to the dimensionless comparative distance and the intensity of the explosive source. Intensity of explosive source R0And the value is any integer between 1 and 10. Each explosion intensity level was determined using the Kinsella method, which was considered based on the ignition energy of each event, the degree of obstruction in the source of the explosion, and the degree of confinement, as shown in the following table:
when the ignition energy of the ignition source is less than 100MJ, the ignition source is called weak ignition energy, and conversely, the ignition source is strong ignition energy; the strong blocking degree represents that the volume of the obstacles is more than 30 percent of the volume of the whole blocking area and the distance between the obstacles is less than 3m, and the weak blocking degree represents that the obstacles exist in the area but cannot simultaneously include the two conditions and does not represent that the area has no obstacles; the existence of the constraint represents that the combustible gas cloud is limited by 2 or 3 solid surfaces, and the existence of the constraint represents that the combustible gas cloud is limited by only one surface on the earth surface.
S325, calculating the overpressure applied to the control chamber according to the dimensionless simulation maximum lateral overpressure, wherein the calculation formula is as follows:
Pcontrol＝ΔPsPBig (a)，
In the formula, PControlFor controlling the overpressure value applied to the chamber, the unit is 0.1MPa, delta PsFor dimensionless comparison of the peak overpressure, the dimensionless comparison of the maximum lateral overpressure Δ P is found from FIG. 2s。
And S33, if the leakage is gas, calculating the overpressure after the chemical explosion and calculating explosion data when the physical explosion occurs. The physical explosion is formed by the sudden change of the state or pressure of a substance, and no chemical reaction occurs in the process of the physical explosion. When a physical explosion occurs in the container, the energy released by the expansion of the gas (i.e., the energy of the explosion) is related not only to the gas pressure and the volume of the container but also to the physical phase of the medium within the container. The medium contained in the container in which the physical explosion occurs is non-thermal gas without enthalpy value and entropy value; the compressed gas is called under the pressure-bearing state. And calculating the overpressure of the physical explosion by using a compressed gas explosion energy calculation model. The method comprises the following steps:
s331, calculating the blasting energy, wherein the calculation formula is as follows:
in the formula, EgIs the blasting energy of the gas, in kJ; pInsulation boardIs the absolute pressure of the gas in the container, in MPa; v is the volume of the container, in m3(ii) a k is the adiabatic index of the gas, namely the ratio of the constant-pressure specific heat to the constant-volume specific heat of the gas, and is 1.4.
S332, calculating the TNT equivalent according to the blasting energy, wherein the calculation formula is as follows:
WTNT＝Eg/QTNT，
in the formula, QTNTIs the explosive heat of TNT in kJ/kg.
S333, calculating an overpressure value applied to the control chamber according to the TNT equivalent, wherein the calculation formula is as follows:
in the formula, PControlThe unit is 0.1MPa for controlling the overpressure value of the chamber; wTNTTNT equivalent of explosive material in kg; r is the distance of the explosion point from the control room in m.
And S4, calculating the positive acting time of the steam cloud explosion shock wave received by the control room according to the explosion data. The calculation formula is as follows:
in the formula, teThe positive acting time of overpressure on the control chamber is s;
the dimensionless positive pressure action time is looked up from fig. 3; e is total energy of gas cloud explosion, and the unit is MJ; pBig (a)Is the local atmospheric pressure, and has the unit of 0.1 MPa; c. C0340m/s is taken as the propagation speed of sound in air.
The invention also provides a control room antiknock analysis system, which comprises:
the danger source data acquisition module is used for acquiring danger source data in a safe distance around the control room;
the leakage material mass flow rate calculation module is used for calculating the leakage material mass flow rate according to the dangerous source data;
the explosion data calculation module is used for calculating explosion data when chemical explosion occurs according to the mass flow rate of the leaked materials when the leaked materials are liquid, and calculating the explosion data when physical explosion occurs when the leaked materials are gas;
and the time calculation module is used for calculating the positive acting time of the steam cloud explosion shock wave received by the control room according to the explosion data.
In the description of the present invention, it is to be understood that the terms "upper", "lower", "left", "right", and the like, indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, are only for convenience in describing the present invention and simplifying the description, and do not indicate or imply that the referred device or element must have a specific orientation and a specific orientation configuration and operation, and thus, should not be construed as limiting the present invention. Furthermore, "first" and "second" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature.
In the description of the present invention, "a plurality" means two or more unless otherwise specified. In the description of the present invention, it should be noted that, unless otherwise explicitly specified or limited, the terms "mounted," "connected," and the like are to be construed broadly and may be, for example, fixedly connected, detachably connected, or integrally connected; can be mechanically or electrically connected; they may be directly connected or indirectly connected through an intermediate member, or they may be connected through two or more elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
The above is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, several modifications and variations can be made without departing from the technical principle of the present invention, and these modifications and variations should also be regarded as the protection scope of the present invention.
Claims (9)
1. A control room antiknock analysis method is characterized by comprising the following steps:
acquiring dangerous source data in a safe distance around a control room;
calculating mass flow rate of the leaked materials according to the dangerous source data;
if the leakage is liquid, calculating explosion data when chemical explosion occurs according to the mass flow rate of the leakage material, and if the leakage is gas, calculating explosion data when physical explosion occurs;
and calculating the positive acting time of the steam cloud explosion shock wave received by the control room according to the explosion data.
2. The control room antiknock analysis method of claim 1, wherein said calculating a mass flow rate of leaked material from said hazard source data comprises:
if the leakage object is liquid, judging whether the leakage scene is storage tank leakage or pipeline leakage, if so, determining a storage tank leakage model, and if so, determining a pipeline leakage model;
and if the leakage object is gas, judging whether the leakage scene is sonic flow or subsonic flow, if so, determining a sonic flow model, and if so, determining a subsonic flow model.
And calculating the mass flow rate of the leaked materials according to the leakage model and the dangerous source data.
3. The control room antiknock analysis method of claim 2, wherein the tank leak model is:
the pipeline leakage model is as follows:
in the formula, QmIs the mass flow rate in kg/s; p1Is the liquid pressure in the storage tank, and the unit is Pa; p0Is ambient pressure in Pa; c1The leakage coefficient of the liquid in the storage tank; a is the area of the leakage hole and is given in m2；P2Is the liquid pressure in the pipeline, and the unit is Pa; ρ is the density of the leaking liquid in kg/m3(ii) a g is the acceleration of gravity; h isLThe height of the liquid above the leakage hole is m; c2Is the leakage coefficient of the liquid in the pipeline;
the sonic flow model is as follows:
the subsonic model is:
in the formula, QmIs the mass flow rate in kg/s; cdIs the gas leakage coefficient; a is the area of the leakage hole in m2；P3Is the gas pressure in the container, and the unit is Pa; m is the molecular weight of the leaking gas or vapor; rgIs an ideal gas constant with the unit of J/(mol.K); t is the gas temperature in K; y is the efflux coefficient;
in the formula, P0To ambient pressure, singlyThe bit is Pa; p3Is the gas pressure in the container, and the unit is Pa; gamma is adiabatic index, gamma is cp/cvCp is constant-pressure heat capacity, and cv is constant-capacity heat capacity.
4. The control room antiknock analysis method according to claim 3, wherein the determining whether the leakage scene is sonic flow or subsonic flow is:
if it is
Then sonic flow, otherwise subsonic flow.
5. The control room antiknock analysis method of claim 3, wherein said calculating explosion data at the time of a chemical explosion based on said mass flow rate of leaking material comprises:
and (3) calculating the TNT equivalent according to the mass flow rate of the leaked material, wherein the calculation formula is as follows:
WTNT＝aA1WfQf/QTNT，
in the formula, a is the ground explosion coefficient; a. the1TNT equivalent coefficient for vapor cloud; wfIs the mass of combustible gas in the vapor cloud, and the unit is kg, Wf＝Qm·tDrain device，tDrain deviceThe time from leakage occurrence to emergency plugging completion is shortened; qfThe unit is kJ/kg of combustion heat of combustible gas; qTNTIs the explosive heat of TNT in kJ/kg;
and calculating the overpressure value of the control chamber according to the TNT equivalent, wherein the calculation formula is as follows:
in the formula, PControlThe unit is 0.1MPa for controlling the overpressure value of the chamber; wTNTTNT equivalent of explosive material in kg; r is the distance of the explosion point from the control room in m.
6. The control room antiknock analysis method of claim 3, wherein said calculating explosion data at the time of a chemical explosion based on said mass flow rate of leaking material comprises:
determining the range of the combustible gas cloud by using a gas diffusion model, wherein the calculation formula is as follows:
Vc＝Q/(ρ×cs)，
in the formula, VcIs the volume of the combustible gas cloud in m3(ii) a Q is the amount of leakage substance in kg, Q ═ Qm·tDrain device，tDrain deviceThe time from leakage occurrence to emergency plugging completion is shortened; rho is the vapor density in kg/m3；csIs volume fraction of fuel mixed with air, and the unit is;
calculating the combustion energy of the fuel-air mixture in the area according to the range of the combustible gas cloud, wherein the calculation formula is as follows:
E＝Vc×3.5×106，
wherein E is the combustion energy of the fuel-air mixture in the explosive source, and is given in J;
calculating a dimensionless, comparable distance for a single equivalent fuel-air mixture by the formula:
wherein R is the distance between an explosion point and a control room and is in the unit of m; pBig (a)Is the local atmospheric pressure, and has the unit of 0.1 MPa;
determining the dimensionless analog maximum lateral overpressure according to the dimensionless analog distance and the intensity of the explosive source;
calculating the overpressure applied to the control chamber according to the dimensionless simulation maximum lateral overpressure, wherein the calculation formula is as follows:
Pcontrol＝ΔPsPBig (a)，
In the formula, PControlThe unit is 0.1MPa for the overpressure value to which the control chamber is subjected.
7. The method according to claim 3, wherein the calculating explosion data when a physical explosion occurs comprises:
and calculating the blasting energy by the following formula:
in the formula, EgIs the blasting energy of the gas, in kJ; pInsulation boardIs the absolute pressure of the gas in the container, in MPa; v is the volume of the container, in m3(ii) a k is the adiabatic index of the gas, namely the ratio of the constant-pressure specific heat to the constant-volume specific heat of the gas;
and (3) calculating the TNT equivalent according to the blasting energy, wherein the calculation formula is as follows:
WTNT＝Eg/QTNT，
in the formula, QTNTIs the explosive heat of TNT in kJ/kg;
and calculating the overpressure value of the control chamber according to the TNT equivalent, wherein the calculation formula is as follows:
in the formula, PControlThe unit is 0.1MPa for controlling the overpressure value of the chamber; wTNTTNT equivalent of explosive material in kg; r is the distance of the explosion point from the control room in m.
8. The control room antiknock analysis method according to claim 5 or 6, wherein the calculation formula for calculating the normal acting time of the steam cloud explosion shock wave received by the control room according to the explosion data is as follows:
in the formula, teThe positive acting time of overpressure on the control chamber is s;
the dimensionless positive pressure action time; e is total energy of gas cloud explosion, and the unit is MJ; pBig (a)Is the local atmospheric pressure, and has the unit of 0.1 MPa; c. C0Is the speed of sound propagation in air.
9. A control room antiknock analysis system, comprising:
the danger source data acquisition module is used for acquiring danger source data in a safe distance around the control room;
the leakage material mass flow rate calculation module is used for calculating the leakage material mass flow rate according to the dangerous source data;
the explosion data calculation module is used for calculating explosion data when chemical explosion occurs according to the mass flow rate of the leaked materials when the leaked materials are liquid, and calculating the explosion data when physical explosion occurs when the leaked materials are gas;
and the time calculation module is used for calculating the positive acting time of the steam cloud explosion shock wave received by the control room according to the explosion data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111675046.XA CN114418362A (en) | 2021-12-31 | 2021-12-31 | Control room anti-knock analysis method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111675046.XA CN114418362A (en) | 2021-12-31 | 2021-12-31 | Control room anti-knock analysis method and system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114418362A true CN114418362A (en) | 2022-04-29 |
Family
ID=81270764
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111675046.XA Pending CN114418362A (en) | 2021-12-31 | 2021-12-31 | Control room anti-knock analysis method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114418362A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116703170A (en) * | 2023-05-24 | 2023-09-05 | 北京市燃气集团有限责任公司 | Investigation analysis method and device for reasons of indoor natural gas explosion accidents |
-
2021
- 2021-12-31 CN CN202111675046.XA patent/CN114418362A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116703170A (en) * | 2023-05-24 | 2023-09-05 | 北京市燃气集团有限责任公司 | Investigation analysis method and device for reasons of indoor natural gas explosion accidents |
| CN116703170B (en) * | 2023-05-24 | 2024-04-12 | 北京市燃气集团有限责任公司 | Investigation analysis method and device for reasons of indoor natural gas explosion accidents |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Shen et al. | Recent application of Computational Fluid Dynamics (CFD) in process safety and loss prevention: A review | |
| Baraldi et al. | Development of a model evaluation protocol for CFD analysis of hydrogen safety issues the SUSANA project | |
| Rigas et al. | Evaluation of hazards associated with hydrogen storage facilities | |
| Le et al. | Safety investigation of hydrogen energy storage systems using quantitative risk assessment | |
| Salzano et al. | The analysis of domino accidents triggered by vapor cloud explosions | |
| CN103605822B (en) | The dangerous dynamic classification method of Petrochemical Wharf liquid chemicals pipe leakage | |
| Khan et al. | MAXCRED–a new software package for rapid risk assessment in chemical process industries | |
| CN114418362A (en) | Control room anti-knock analysis method and system | |
| CN107301495A (en) | A kind of appraisal procedure of organic rankine cycle system mixed working fluid leakage security | |
| Ustolin et al. | Computational fluid dynamics modeling of liquid hydrogen release and dispersion in gas refuelling stations | |
| Gu et al. | Experimental and numerical study on deflagration characteristics of large-scale propane–air mixture | |
| Tulach et al. | CFD simulation of vented explosion and turbulent flame propagation | |
| Tan et al. | Enhanced mitigation of fire and explosion risks due to hydrogen leakage using targeted nitrogen nozzle spray approach | |
| Kim et al. | On the application of CFD codes for natural gas dispersion and explosion in gas fuelled ship | |
| CN113821765A (en) | Risk-based building anti-explosion value evaluation method and system | |
| Büyükkıdan et al. | The Risk Calculation of Hazardous Zones Created By Flammable And Explosive Chemicals, LPG Tank Example | |
| Baraldi et al. | Gap analysis of CFD modelling of accidental hydrogen release and combustion | |
| Patel et al. | Hazard evaluation using aloha tool in storage area of an oil refinery | |
| Miller et al. | New criteria for safety distances during pneumatic pressure testing of vessels and pipes | |
| Çetinyokuş et al. | Consequence Analysis of An Industrial Accident at a Fuel Station | |
| Blaylock et al. | Analyses in Support of Risk-Informed Natural Gas Vehicle Maintenance Facility Codes and Standards: Phase II | |
| Sulaiman | Gas Explosion Characteristics in Confined Straight and 90 Degree Bend Pipes | |
| Czujko | Consequences of explosions in various industries | |
| Pourgol-Mohammad et al. | Simulation of fuel storage hazards associated with petroleum refineries tanks | |
| Baraldi et al. | Prioritisation of Research and Development for modelling the safe production, storage, delivery and use of hydrogen |
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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 210000 Room 401, building e, No. 606, ningliu Road, Jiangbei new area, Nanjing, Jiangsu Province Applicant after: Nanjing Nangong emergency Technology Co.,Ltd. Address before: 210000 Room 401, building e, No. 606, ningliu Road, Jiangbei new area, Nanjing, Jiangsu Province Applicant before: Nanjing South Industrial Safety Technology Co.,Ltd. |