CN113094908A - Combustible gas pipeline explosion intensity deduction method based on section trace measurement - Google Patents

Combustible gas pipeline explosion intensity deduction method based on section trace measurement Download PDF

Info

Publication number
CN113094908A
CN113094908A CN202110409224.8A CN202110409224A CN113094908A CN 113094908 A CN113094908 A CN 113094908A CN 202110409224 A CN202110409224 A CN 202110409224A CN 113094908 A CN113094908 A CN 113094908A
Authority
CN
China
Prior art keywords
pipeline
explosion
gas
combustible gas
method based
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.)
Granted
Application number
CN202110409224.8A
Other languages
Chinese (zh)
Other versions
CN113094908B (en
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.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
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 China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN202110409224.8A priority Critical patent/CN113094908B/en
Publication of CN113094908A publication Critical patent/CN113094908A/en
Application granted granted Critical
Publication of CN113094908B publication Critical patent/CN113094908B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/14—Pipes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Abstract

The invention discloses a combustible gas pipeline explosion intensity deduction method based on section trace measurement, which comprises the following steps: analyzing the stepped traces formed on the cross section of the combustible gas pipeline after explosion, and determining the crack initiation position according to the opposite convex directions of the stepped traces on the two sides; measuring the spacing values of at least three groups of stepped traces from the crack initiation position to two sides respectively, and calculating the average spacing of the stepped traces at two sides respectively
Figure DDA0003023509770000011
And calculating the detonation pressure of the mixed gas in the pipeline according to a formula, namely the explosion intensity of the pipeline. The method can realize the deduction of the explosion intensity without finite element analysis, has high accuracy of the deduction result and can realize quantitative calculation.

Description

Combustible gas pipeline explosion intensity deduction method based on section trace measurement
Technical Field
The invention relates to the field of inversion of pipeline explosion accidents, in particular to a combustible gas pipeline explosion intensity deduction method based on section trace measurement.
Background
The pipeline is widely applied to the process industries of chemical industry, nuclear power, petroleum and the like, and combustible and explosive dangerous liquid or gas is often conveyed inside the pipeline. For the combustible gas pipeline, because air or oxygen may be mixed in the pipeline due to careless operation or other external factors, combustible mixed gas is formed in the pipeline, the mixed gas in the pipeline can be ignited under certain conditions, and the flame forms deflagration and detonation after accelerating in the pipeline, so that the pipeline explosion accident is caused, and serious harm is caused. For example, in the pipeline explosion accident of the Hamaoka nuclear power station in japan, the investigation result shows that the pipeline is ignited by hydrogen and oxygen to explode.
The important point in the analysis and investigation of pipeline explosion accidents is the inversion of explosion intensity, and then the cause and development process of the accidents are deduced. For the explosion of a combustible gas pipeline, the inversion of the explosion intensity is mainly the inversion of the detonation pressure because the interior tends to form detonation. At present, the inversion of the explosion intensity mainly adopts an empirical or semi-empirical method, and the deduction result is often not accurate enough and is easily influenced by the professional degree and experience of personnel according to the damage degree of pipelines, surrounding buildings and other structures and lack of clear and quantitative data basis.
Recent research shows that a 'step' -shaped trace often exists on the section of the pipeline after explosion, but the corresponding relation between the relevant characteristics of the trace and the detonation pressure or the explosion intensity is not clear, and an effective method for deducing the explosion intensity of the pipeline according to the trace characteristics of the section of the pipeline is not available at present.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a combustible gas pipeline explosion intensity deduction method based on section trace measurement, the deduction result is high in accuracy, and quantitative calculation can be realized.
A combustible gas pipeline explosion intensity deduction method based on section trace measurement comprises the following steps:
(1) the protrusions of the stepped traces on two sides of the crack initiation position on the pipeline section are opposite in direction, so that the stepped traces formed on the exploded combustible gas pipeline section are analyzed, and the crack initiation position is determined according to the opposite protrusion directions of the stepped traces on the two sides;
(2) because the initiation position and the crack initiation position are often not at the same axial position, the expansion rates of the two crack tips are different due to the asymmetry of the load, and the distance between the formed step-shaped traces is always large on one side and small on the other side. Therefore, the distance values of at least three groups of stepped traces are measured from the crack initiation position to two sides respectively, and the distance between the stepped traces on the two sides of the crack initiation position is defined as d1、d2Respectively calculating the average distance between the stepped traces on both sides
Figure BDA0003023509750000021
And is
Figure BDA0003023509750000022
(3) Looking up the technological parameters of pipeline operation to obtain the type and composition of medium/gas in pipeline and determine the absolute pressure p of gas in pipeline before explosion0And a thermodynamic temperature T. Calculating the average molar mass of the mixed gas according to the type and the volume fraction of the mixed gas in the pipeline
Figure BDA0003023509750000023
Figure BDA0003023509750000024
In the formula, MiIs the molar mass of the i-th gas, niIs the volume fraction or mole fraction of the ith gas.
Calculating the detonation pressure p of the mixed gas in the pipeline according to the following formula, namely the explosion intensity of the pipeline:
Figure BDA0003023509750000025
wherein p is0The absolute pressure of gas in the pipeline before explosion, T is thermodynamic temperature, E, rho and v are respectively elastic modulus, density and Poisson ratio of the pipe, and the absolute pressure is obtained by looking up a material mechanical property manual or a database or sampling the pipeline after explosion and carrying out a material mechanical property experiment; r is an ideal gas constant which is 8.314J/(mol.K); gamma is the adiabatic index of the mixed gas in the pipeline before explosion, the value of the mixed gas can be 1.4 for common gases such as hydrogen, methane, air and the like, and other types or mixed gases can be obtained by searching data; d is the middle diameter of the pipeline, and is obtained by looking up the design data of the pipeline or measuring the pipeline after explosion; alpha is the ratio of the mean rates of axial expansion of the two split tips, alpha<1, the value is usually in the range of 0.70 to 0.98, and in most cases, α may take a value of 0.95.
The invention has the following beneficial effects:
(1) the deduction method has clear steps, relevant physical parameters are clear and easy to obtain, the deduction of the explosion intensity depends on quantitative calculation, and the accuracy of the deduction result is high;
(2) the deduction method has low requirement on the professional degree of the implementer, only needs to measure the step-shaped trace distance of the pipeline and further calculates according to a formula without complex measuring equipment or tools and developing finite element simulation analysis.
Drawings
FIG. 1 is a schematic view of a combustible gas conduit of an embodiment of the invention;
FIG. 2 is a schematic illustration of a cross-sectional stepped trace of a combustible gas conduit according to an embodiment of the invention.
Detailed Description
The present invention will be described in detail below with reference to the accompanying drawings and preferred embodiments, and the objects and effects of the present invention will become more apparent, it being understood that the specific embodiments described herein are merely illustrative of the present invention and are not intended to limit the present invention.
In the embodiment, the 6061-T6 aluminum alloy pipeline shown in FIG. 1 is taken as an object, and the mixed gas of ethylene and oxygen inside the pipeline is detonated, so that the pipeline is broken. According to the implementation steps of the invention, the explosion intensity is deduced.
1. Determining geometric parameters and mechanical properties of material
The pipe was measured to have an outside diameter of 41.28mm and a wall thickness of 0.89mm, where the diameter D was 41.28-0.89-40.39 mm. Looking up the mechanical property parameters of 6061-T6 aluminum alloy material to obtain the elastic modulus E of the material as 6.90 multiplied by 1010Pa, density rho 2780kg/m3And poisson's ratio v ═ 0.33.
2. Determining internal gas thermodynamic state parameters before pipeline explosion
Inquiring the experimental conditions of the pipeline to obtain the stoichiometric ratio of ethylene to oxygen in the pipeline as 1: 3 initial pressure p of mixed gas0180kPa, temperature T298K, the average molar mass of the mixed gas is calculated as:
Figure BDA0003023509750000031
3. determining the average distance between the stepped traces on two sides of the crack initiation position of the pipeline section
Figure BDA0003023509750000032
Observing the stepped traces on the section of the pipeline, respectively measuring the distances between the stepped traces on two sides of the 3 groups of crack initiation positions and averaging to obtain the average distance between the stepped traces on two sides
Figure BDA0003023509750000033
As shown in fig. 2.
4. Calculating the explosion intensity of the pipeline
And substituting the data into the following formula to calculate the detonation pressure of the pipeline, wherein the adiabatic index gamma of the mixed gas is 1.4, and the adiabatic index alpha of the mixed gas is 0.95.
Figure BDA0003023509750000034
The detonation pressure of the gas in the pipeline measured by an actual experiment is 6.1MPa, and the error between the deduction result and the actual value is-8.2%. The invention provides the result with better accuracy, considering that the precise inversion of the explosion intensity of the pipeline is always extremely difficult work.
It will be understood by those skilled in the art that the foregoing is only a preferred embodiment of the present invention, and is not intended to limit the invention, and although the invention has been described in detail with reference to the foregoing examples, it will be apparent to those skilled in the art that various changes in the form and details of the embodiments may be made and equivalents may be substituted for elements thereof. All modifications, equivalents and the like which come within the spirit and principle of the invention are intended to be included within the scope of the invention.

Claims (4)

1. A combustible gas pipeline explosion intensity deduction method based on section trace measurement is characterized by comprising the following steps:
(1) analyzing the stepped traces formed on the cross section of the combustible gas pipeline after explosion, and determining the crack initiation position according to the opposite convex directions of the stepped traces on the two sides;
(2) measuring the spacing values of at least three groups of stepped traces from the crack initiation position to two sides respectively, and calculating the average spacing of the stepped traces at two sides respectively
Figure FDA0003023509740000011
And is
Figure FDA0003023509740000012
(3) Calculating the detonation pressure p of the mixed gas in the pipeline according to the following formula, namely the explosion intensity of the pipeline:
Figure FDA0003023509740000013
Figure FDA0003023509740000014
wherein p is0Is the absolute pressure of the gas in the pipeline before explosion, T is the thermodynamic temperature,
Figure FDA0003023509740000015
is the average molar mass, M, of the gas mixture in the conduitiIs the molar mass of the i-th gas, niIs the volume fraction or mole fraction of the ith gas; E. rho and ν are respectively the elastic modulus, density and Poisson's ratio of the pipe, R is an ideal gas constant which is 8.314J/(mol. K); gamma is the adiabatic index of the mixed gas in the pipeline before explosion; alpha is the ratio of the mean rates of axial expansion of the two split tips, alpha<1; d is the middle diameter of the pipeline.
2. The combustible gas pipeline explosion intensity deduction method based on section mark measurement as claimed in claim 1, wherein the value range of α is 0.70-0.98.
3. The combustible gas pipeline explosion intensity deduction method based on section mark measurement according to claim 1, wherein α is 0.95.
4. The combustible gas pipeline explosion intensity deduction method based on section trace measurement according to claim 1, characterized in that the values of gamma for hydrogen, methane and air are 1.4.
CN202110409224.8A 2021-04-16 2021-04-16 Combustible gas pipeline explosion intensity deduction method based on section trace measurement Active CN113094908B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110409224.8A CN113094908B (en) 2021-04-16 2021-04-16 Combustible gas pipeline explosion intensity deduction method based on section trace measurement

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110409224.8A CN113094908B (en) 2021-04-16 2021-04-16 Combustible gas pipeline explosion intensity deduction method based on section trace measurement

Publications (2)

Publication Number Publication Date
CN113094908A true CN113094908A (en) 2021-07-09
CN113094908B CN113094908B (en) 2022-03-08

Family

ID=76678465

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110409224.8A Active CN113094908B (en) 2021-04-16 2021-04-16 Combustible gas pipeline explosion intensity deduction method based on section trace measurement

Country Status (1)

Country Link
CN (1) CN113094908B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013104461A1 (en) * 2012-01-13 2013-07-18 L'air Liquide,Societe Anonyme Pour L'etude Et L'exploitation Des Procedes Georges Claude Method for the preparation of compressed oxidant-fuel gas mixtures
CN110334465A (en) * 2019-07-15 2019-10-15 中国石油大学(华东) It is a kind of consider damage evolution gas burst under pipeline Dynamic Fracture prediction technique
CN111191382A (en) * 2020-01-09 2020-05-22 中国石油大学(华东) Method for calculating progressive crack propagation length of metal pipeline in front and back directions under internal explosion
CN111751074A (en) * 2020-07-01 2020-10-09 中国科学院力学研究所 Detonation-driven high-enthalpy shock tunnel automatic inflation control system
WO2021012136A1 (en) * 2019-07-22 2021-01-28 中国石油大学(华东) Method for calculating internal explosion load rate based on progressive extension distance of pipeline cracks

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013104461A1 (en) * 2012-01-13 2013-07-18 L'air Liquide,Societe Anonyme Pour L'etude Et L'exploitation Des Procedes Georges Claude Method for the preparation of compressed oxidant-fuel gas mixtures
CN110334465A (en) * 2019-07-15 2019-10-15 中国石油大学(华东) It is a kind of consider damage evolution gas burst under pipeline Dynamic Fracture prediction technique
WO2021012136A1 (en) * 2019-07-22 2021-01-28 中国石油大学(华东) Method for calculating internal explosion load rate based on progressive extension distance of pipeline cracks
CN111191382A (en) * 2020-01-09 2020-05-22 中国石油大学(华东) Method for calculating progressive crack propagation length of metal pipeline in front and back directions under internal explosion
CN111751074A (en) * 2020-07-01 2020-10-09 中国科学院力学研究所 Detonation-driven high-enthalpy shock tunnel automatic inflation control system

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
BO LI .ETC: ""Effect of steam explosion on dietary fiber, polysaccharide, protein and physicochemical properties of okara"", 《FOOD HYDROCOLLOIDS》 *
LIUTAO SUN .ETC: ""Effects of changes in pipe cross-section on the explosion-proof distance and the propagation characteristics of gas explosions"", 《JOURNAL OF NATURAL GAS SCIENCE AND ENGINEERING》 *
YANG DU .ETC: ""Consequence analysis of premixed flammable gas explosion occurring in pipe using a coupled fluid-structure-fracture approach"", 《JOURNAL OF LOSS PREVENTION IN THE PROCESS INDUSTRIES》 *
YANG DU .ETC: ""Suppressions of gasoline-air mixture explosion by non-premixed nitrogen in a closed tunnel"", 《JOURNAL OF LOSS PREVENTION IN THE PROCESS INDUSTRIES》 *
付传雄 等: ""基于弹性模量缩减法的钢衬钢筋混凝土压力管道极限承载力计算"", 《水利水电科技进展》 *
张宏彬: ""管道气压试验安全技术与管理"", 《海洋工程装备与技术》 *
杜洋 等: ""考虑流固耦合的管道爆炸后果预测与分析"", 《浙江大学学报（工学版）》 *

Also Published As

Publication number Publication date
CN113094908B (en) 2022-03-08

Similar Documents

Publication Publication Date Title
Movileanu et al. Explosion of gaseous ethylene–air mixtures in closed cylindrical vessels with central ignition
CN107209164B (en) Method for calculating methane value and apparatus for measuring methane value
Wang et al. Effects of flammable gases on the explosion characteristics of CH4 in air
Mitu et al. Influence of initial pressure and vessel’s geometry on deflagration of stoichiometric methane–air mixture in small-scale closed vessels
Jones Instrument Technology: Measurement of pressure, level, flow and temperature
CN105628752A (en) Calibration method of electrochemical sensor
CN102589809A (en) Portable leak detector calibration system and method
CN113094908B (en) Combustible gas pipeline explosion intensity deduction method based on section trace measurement
Liao et al. Approximation of flammability region for natural gas–air–diluent mixture
McElroy et al. Compression-factor measurements on methane, carbon dioxide, and (methane+ carbon dioxide) using a weighing method
Budenholzer et al. Phase Equlibria in Hydrocarbons Systems
KR102603253B1 (en) Direct combustion calorimetry system of flare gas
CN110414822B (en) Method for calculating internal explosion load rate according to pipeline crack progressive propagation distance
CN101726441A (en) Method for detecting low-temperature brittleness of zinc anode ribbon
WO2021012136A1 (en) Method for calculating internal explosion load rate based on progressive extension distance of pipeline cracks
CN110378075A (en) A kind of acid and nonacid natural gas water content amount prediction technique
US11835440B2 (en) Multi-capillary force curve averaging method based on multi-sample overall virtual measurement
US10371678B2 (en) Method and measuring apparatus for determining gas properties by correlation
Shu et al. Experimental and modeling studies on the correlation between auto-ignition delays and the methane number of Liquefied Natural Gas (LNG) and Liquefied Biogas (LBG)
Schildberg et al. Experimental Determination of the Static Equivalent Pressures of Detonative Explosions of Cyclohexane/O 2/N 2-Mixtures in Long and Short Pipes (part 1of 3).
CN205808782U (en) A kind of vacuum withdraw device of short neck nozzle
Dielenschneider et al. Some guidelines for the experimental characterization of turbocharger compressors
Gan Flammability Characteristics of Light Hydrocarbons and Their Mixtures at Elevated Conditions
CN113218816B (en) Method for testing density of unknown gas at different temperatures by using thermogravimetric method
US6094993A (en) Method for measuring flow rate of a fluid in a conduit

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant