CN114418362A - Control room anti-knock analysis method and system - Google Patents

Control room anti-knock analysis method and system Download PDF

Info

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
Application number
CN202111675046.XA
Other languages
Chinese (zh)
Inventor
张明广
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing South Industrial Safety Technology Co ltd
Original Assignee
Nanjing South Industrial Safety Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing South Industrial Safety Technology Co ltd filed Critical Nanjing South Industrial Safety Technology Co ltd
Priority to CN202111675046.XA priority Critical patent/CN114418362A/en
Publication of CN114418362A publication Critical patent/CN114418362A/en
Pending legal-status Critical Current

Links

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

Control room anti-knock analysis method and system
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:
Figure BDA0003450878650000021
the pipeline leakage model is as follows:
Figure BDA0003450878650000022
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:
Figure BDA0003450878650000023
the subsonic model is:
Figure BDA0003450878650000024
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;
Figure BDA0003450878650000025
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
Figure BDA0003450878650000026
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:
Figure BDA0003450878650000031
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:
Figure BDA0003450878650000032
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:
Figure BDA0003450878650000041
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:
Figure BDA0003450878650000042
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:
Figure BDA0003450878650000043
in the formula, teThe positive acting time of overpressure on the control chamber is s;
Figure BDA0003450878650000044
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:
Figure BDA0003450878650000051
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
Figure BDA0003450878650000061
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:
Figure BDA0003450878650000062
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:
Figure BDA0003450878650000063
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:
Figure BDA0003450878650000064
the subsonic model is:
Figure BDA0003450878650000071
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.
Figure BDA0003450878650000072
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:
Figure BDA0003450878650000073
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.
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:
Figure BDA0003450878650000081
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:
Figure BDA0003450878650000091
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:
Figure BDA0003450878650000092
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:
Figure BDA0003450878650000101
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:
Figure BDA0003450878650000102
in the formula, teThe positive acting time of overpressure on the control chamber is s;
Figure BDA0003450878650000103
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:
Figure FDA0003450878640000011
the pipeline leakage model is as follows:
Figure FDA0003450878640000012
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:
Figure FDA0003450878640000013
the subsonic model is:
Figure FDA0003450878640000014
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;
Figure FDA0003450878640000021
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
Figure FDA0003450878640000022
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:
Figure FDA0003450878640000023
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:
Figure FDA0003450878640000031
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:
Figure FDA0003450878640000032
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:
Figure FDA0003450878640000033
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:
Figure FDA0003450878640000041
in the formula, teThe positive acting time of overpressure on the control chamber is s;
Figure FDA0003450878640000042
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.
CN202111675046.XA 2021-12-31 2021-12-31 Control room anti-knock analysis method and system Pending CN114418362A (en)

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)

* Cited by examiner, † Cited by third party
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

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103605822A (en) * 2013-10-08 2014-02-26 交通运输部天津水运工程科学研究所 Dynamic classification method of risk of liquid chemical pipeline leakage of petrochemical wharf
CN111539126A (en) * 2020-05-14 2020-08-14 山东慧泰智能科技有限公司 Chemical storage tank steam cloud explosion accident simulation analysis method and system
CN112668752A (en) * 2020-11-27 2021-04-16 合肥泽众城市智能科技有限公司 Safe evacuation method and system for gas pipeline leakage emergency

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103605822A (en) * 2013-10-08 2014-02-26 交通运输部天津水运工程科学研究所 Dynamic classification method of risk of liquid chemical pipeline leakage of petrochemical wharf
CN111539126A (en) * 2020-05-14 2020-08-14 山东慧泰智能科技有限公司 Chemical storage tank steam cloud explosion accident simulation analysis method and system
CN112668752A (en) * 2020-11-27 2021-04-16 合肥泽众城市智能科技有限公司 Safe evacuation method and system for gas pipeline leakage emergency

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
伍东;宋文华;张茹;王鹏;: "火电厂氢气储罐火灾爆炸危险性分析", 消防科学与技术, no. 11, 15 November 2008 (2008-11-15), pages 847 - 851 *
吴月华;杨林青;刘迁;王天惠;李林;: "泄漏与蒸气云爆炸模型在氢气管道定量评价中的应用", 山东化工, no. 07, 8 April 2015 (2015-04-08), pages 183 - 186 *
夏兰生;: "中心控制室抗爆模型选择与计算", 石油化工自动化, no. 02, 20 April 2015 (2015-04-20), pages 51 - 55 *
张少鹏;刘秀琴;: "石化企业中心控制室的抗爆需求探讨", 石油化工自动化, no. 05, 20 October 2013 (2013-10-20), pages 11 - 15 *
王德国;: "基于管道爆炸数值模拟的架空天然气管道并行间距研究", 中国石油大学学报(自然科学版), no. 05, 20 October 2013 (2013-10-20), pages 183 - 186 *

Cited By (2)

* Cited by examiner, † Cited by third party
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
Le et al. Safety investigation of hydrogen energy storage systems using quantitative risk assessment
Rigas et al. Evaluation of hazards associated with hydrogen storage facilities
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
Pang et al. Characteristics of external explosions induced by vented hydrogen deflagration
CN114418362A (en) Control room anti-knock analysis method and system
CN102096768B (en) Dangerous chemical leakage pollution evaluation method
CN107301495A (en) A kind of appraisal procedure of organic rankine cycle system mixed working fluid leakage security
Qi et al. Flexible alarming mechanism of a general GDS deployment for explosive accidents caused by gas leakage
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
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
Çetinyokuş et al. Consequence Analysis of An Industrial Accident at a Fuel Station
Büyükkıdan et al. The Risk Calculation of Hazardous Zones Created By Flammable And Explosive Chemicals, LPG Tank Example
CN113821765A (en) Risk-based building anti-explosion value evaluation method and system
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
Le et al. Simulation of the LPG leakage and dispersion process to the factory environment using computational fluid dynamics (CFD)
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
Bull Review of large‐scale explosion experiments
Pourgol-Mohammad et al. Simulation of fuel storage hazards associated with petroleum refineries tanks

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.