Explosion accident source item simulation analysis method and device

Abstract

The invention discloses an explosion accident source item simulation analysis method and device, and relates to the technical field of explosion accident source item analysis, and the method establishes a set of relatively complete explosion accident source item analysis model, and performs comprehensive and complete analysis simulation on explosion to obtain an explosion accident source item. The device can simulate the distribution conditions of light radiation (light impulse), shock waves and stable clouds after an explosion accident occurs, and provides a basis for evaluating and simulating propagation, influence and prevention after the explosion accident occurs. Therefore, a decision maker can set corresponding protection measures according to the obtained source item parameters, and a scientific basis is provided for emergency response, pollution control and public protection after an explosion accident occurs. Compared with traditional single source item simulation, the method has the advantages that rich and diversified source item data can be obtained through one model, and the model setting cost and the time cost needed by calculation are reduced.

Description

Technical Field

[0001] The present invention relates to the technical field of explosion accident source term analysis, and in particular to an explosion accident source term simulation analysis method and device capable of performing comprehensive and complete analysis. Background Art

[0002] Explosions are a frequent occurrence, and improper handling can lead to immeasurable consequences. Explosion source information is crucial for determining the severity of explosions, assessing their consequences, and developing emergency response measures. However, in the event of a severe explosion, timely and accurate information on the source information may not be available, directly impacting the assessment of consequences and, in particular, the emergency response process, where it serves as a crucial basis for decision-making.

[0003] The core goal of explosion source analysis is to quantitatively assess the scale of material releases. Source analysis of explosions is complex, with source parameters primarily including light radiation, shock waves, and smoke clouds. However, only a few studies have examined a small number of explosion source parameters, and a comprehensive and complete explosion source analysis model has yet to be established.

[0004] Therefore, how to establish a complete explosion accident source analysis model to provide a scientific basis for emergency response, pollution control and public protection after the explosion accident is a technical problem that urgently needs to be solved by technical personnel in this field. Summary of the Invention

[0005] In view of the above problems, the present invention provides a method and device for simulating and analyzing explosion accident source items for overcoming the above problems or at least partially solving the above problems.

[0006] The present invention provides the following solutions:

[0007] A method for simulating and analyzing explosion accident sources, comprising:

[0008] Obtaining basic parameter data of the explosion accident, wherein the basic parameter data includes at least explosion equivalent, explosion mode, explosion type, distance from the explosion center, and explosion height;

[0009] Inputting the basic parameter data into an explosion accident source term analysis model to obtain explosion accident source term simulation analysis results, wherein the explosion accident source term simulation analysis results at least include light radiation damage characteristic quantities, shock wave damage characteristic quantities, and stable cloud distribution;

[0010] The explosion accident source analysis model includes at least a light impulse calculation module, a shock wave damage calculation module and a stable cloud computing module;

[0011] The light impulse calculation module is used to calculate the light impulse and the influence of terrain and objects on light radiation by using the explosion equivalent, the explosion mode and the distance from the explosion center to obtain the light radiation damage characteristic value;

[0012] The shock wave damage calculation module is used to perform air shock wave overpressure calculation, ground explosion shock wave overpressure calculation, and dynamic pressure calculation using the explosion equivalent, the explosion mode, and the distance from the explosion center to obtain the shock wave damage characteristic value;

[0013] The stable cloud computing module is used to describe the stable cloud shape parameters and the sedimentation ash mass parameters, and is used to obtain the distribution of the stable cloud.

[0014] Preferably, the light impulse calculation is used to perform the following operations:

[0015] determining whether to consider attenuation caused by atmospheric absorption and scattering according to the distance from the explosion center;

[0016] If only geometric attenuation is considered and the attenuation caused by atmospheric absorption and scattering is ignored, and the explosion mode is determined to be airburst, the light impulse is calculated using the following formula:

[0017]

[0018] Where: f G represents the light equivalent coefficient; Q represents the explosion equivalent; r represents the distance from the explosion center;

[0019] If only geometric attenuation is considered and the attenuation caused by atmospheric absorption and scattering is ignored, and the explosion method is determined to be a ground explosion, the light impulse is calculated using the following formula:

[0020]

[0021] Where: h B Indicates relative height, f d Indicates the effective light equivalent coefficient;

[0022] Considering the attenuation caused by atmospheric absorption and scattering, and determining that the explosion mode is airburst, the light impulse is calculated using the following formula:

[0023]

[0024] Where: τ represents the atmospheric transmittance;

[0025] Considering the attenuation caused by atmospheric absorption and scattering, and determining that the explosion method is a ground explosion, the light impulse is calculated using the following formula:

[0026]

[0027] Preferably, the atmospheric transmittance τ is expressed by the following formula:

[0028] τ=e -μr +0.23(μr) 1.6 e -0.65μr

[0029] Where: e -μr Indicates the transmittance of direct light; 0.23 (μr) 1.6 e -0.65μr represents the enhancement term of multiple scattering, and μ represents the average atmospheric attenuation coefficient.

[0030] Preferably, the air shock wave overpressure calculation is expressed by the following formula:

[0031]

[0032] The calculation of the ground explosion shock wave overpressure is expressed by the following formula:

[0033]

[0034] The dynamic pressure calculation is expressed by the following formula:

[0035]

[0036] Preferably, the ground explosion shock wave overpressure calculation is performed after the following judgment condition is met:

[0037]

[0038] Preferably, the parameters describing the stable cloud shape include:

[0039] According to the observation data, the height of the stable cloud cap is fitted, H B and top height H T It is expressed by the following formula:

[0040] H B =aQ b

[0041] H T =cQ d

[0042] Where Q represents the explosion equivalent, H B and H T The unit is meter, and its parameters are:

[0043] a=2228,b=0.3463; Q≤4.07kt

[0044] a=2661,b=0.2198; Q>4.07kt

[0045] c=3597,d=0.2553; Q<2.29kt

[0046] c=3170, d=0.4077; 2.29kt≤Q<19kt

[0047] c=6474,d=0.1650; Q≥19kt

[0048] Radius R of the stable cloud cap c The fitting formula is expressed as:

[0049] R C =exp[6.7553+0.32055ln Q+0.01137478(lnQ) 2 ].

[0050] Preferably: the settling ash mass parameter includes the mass of soil plus shell fragments;

[0051] If it is determined to be an underground explosion, the mass parameter of the settled ash is expressed by the following formula:

[0052]

[0053] R=112.5+0.755d-9.6×10 -6 d 3 -9.11×10 -12 d 5

[0054] S=32.7+0.851d-2.52×10 -5 d 3 -1.78×10 -10 d 5

[0055] Where: d represents the specific depth;

[0056] If it is determined to be a ground explosion, the mass parameter of the settled ash is expressed by the following formula:

[0057]

[0058] Where: λ represents the specific height.

[0059] A device for simulating and analyzing an explosion accident source term, used to perform the above-mentioned method for simulating and analyzing an explosion accident source term, comprising:

[0060] A basic parameter acquisition unit is used to acquire basic parameter data of the explosion accident, wherein the basic parameter data at least includes explosion equivalent, explosion mode, explosion type, distance from the explosion center, and explosion height;

[0061] an analysis result acquisition unit, configured to input the basic parameter data into an explosion accident source term analysis model to obtain an explosion accident source term simulation analysis result, wherein the explosion accident source term simulation analysis result at least includes a light radiation damage characteristic quantity, a shock wave damage characteristic quantity, and a distribution of stable clouds;

[0062] The explosion accident source analysis model includes at least a light impulse calculation module, a shock wave damage calculation module and a stable cloud computing module;

[0063] The light impulse calculation module is used to calculate the light impulse and the influence of terrain and objects on light radiation by using the explosion equivalent, the explosion mode and the distance from the explosion center to obtain the light radiation damage characteristic value;

[0064] The shock wave damage calculation module is used to perform air shock wave overpressure calculation, ground explosion shock wave overpressure calculation, and dynamic pressure calculation using the explosion equivalent, the explosion mode, and the distance from the explosion center to obtain the shock wave damage characteristic value;

[0065] The stable cloud computing module is used to describe the stable cloud shape parameters and the sedimentation ash mass parameters, and is used to obtain the distribution of the stable cloud.

[0066] According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects:

[0067] The embodiment of the present application provides a method and device for simulating and analyzing the source term of an explosion accident, which establishes a relatively complete set of explosion accident source term analysis models. By conducting a comprehensive and complete analysis and simulation of the explosion, it is possible to simulate the distribution of light radiation (light impulse), shock waves, and stable clouds after the explosion accident occurs, and provide a basis for evaluating and simulating the propagation, impact, and prevention of explosion accidents. This allows decision makers to set corresponding protective measures based on the various source term parameters obtained, providing a scientific basis for emergency response, pollution control, and public protection after the explosion accident occurs. Compared with the traditional single source term simulation, this method can obtain rich and diverse source term data with one model, reducing the model setting cost and the time cost required for calculation.

[0068] Of course, any product implementing the present invention does not necessarily need to achieve all of the advantages described above at the same time. BRIEF DESCRIPTION OF THE DRAWINGS

[0069] To more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be derived from these drawings without inventive effort.

[0070] Figure 1 This is a flow chart of a method for simulating and analyzing explosion accident sources provided by an embodiment of the present invention;

[0071] Figure 2 Schematic diagram of the shielding effect of terrain on light radiation provided by an embodiment of the present invention;

[0072] Figure 3 This is a flow chart for calculating the characteristic radiation damage quantity, light impulse U, provided by an embodiment of the present invention;

[0073] Figure 4 This is a flow chart for calculating the shock wave damage characteristic quantities overpressure Δp and dynamic pressure q provided by an embodiment of the present invention;

[0074] Figure 5 Schematic diagram of the geometric dimensions of a stable smoke cloud provided by an embodiment of the present invention;

[0075] Figure 6 Schematic diagram of an explosion accident source term simulation and analysis device provided by an embodiment of the present invention;

[0076] Figure 7 It is a schematic diagram of an explosion accident source term simulation and analysis device provided by an embodiment of the present invention. DETAILED DESCRIPTION

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

[0078] See also Figure 1 , is a method for simulating and analyzing explosion accident sources provided by an embodiment of the present invention, such as Figure 1 As shown, the method may include:

[0079] S101: Obtain basic parameter data of the explosion accident, wherein the basic parameter data at least includes explosion equivalent, explosion mode, explosion type, distance from the explosion center, and explosion height;

[0080] S102: Inputting the basic parameter data into an explosion accident source term analysis model to obtain explosion accident source term simulation analysis results, wherein the explosion accident source term simulation analysis results at least include light radiation damage characteristic quantities, shock wave damage characteristic quantities, and stable cloud distribution;

[0081] The explosion accident source analysis model includes at least a light impulse calculation module, a shock wave damage calculation module and a stable cloud computing module;

[0082] The light impulse calculation module is configured to use the explosion equivalent, the explosion mode, and the distance from the explosion center to perform light impulse calculation and calculate the influence of terrain and objects on light radiation to obtain the light radiation damage characteristic value. In specific implementation, the embodiment of the present application can provide the light impulse calculation module to perform the following operations:

[0083] determining whether to consider attenuation caused by atmospheric absorption and scattering according to the distance from the explosion center;

[0084] If only geometric attenuation is considered and the attenuation caused by atmospheric absorption and scattering is ignored, and the explosion mode is determined to be airburst, the light impulse is calculated using the following formula:

[0085]

[0086] Where: f G represents the light equivalent coefficient; Q represents the explosion equivalent; r represents the distance from the explosion center;

[0087] If only geometric attenuation is considered and the attenuation caused by atmospheric absorption and scattering is ignored, and the explosion method is determined to be a ground explosion, the light impulse is calculated using the following formula:

[0088]

[0089] Where h B Indicates relative height, f d Indicates the effective light equivalent coefficient;

[0090] Considering the attenuation caused by atmospheric absorption and scattering, and determining that the explosion mode is airburst, the light impulse is calculated using the following formula:

[0091]

[0092] Where: τ represents the atmospheric transmittance;

[0093] Considering the attenuation caused by atmospheric absorption and scattering, and determining that the explosion method is a ground explosion, the light impulse is calculated using the following formula:

[0094]

[0095] Furthermore, the atmospheric transmittance τ is expressed by the following formula:

[0096] τ=e -μr +0.23(μr) 1.6 e -0.65μr

[0097] Where: e -μr Indicates the transmittance of direct light; 0.23 (μr)1.6 e -0.65μr represents the enhancement term of multiple scattering, and μ represents the average atmospheric attenuation coefficient.

[0098] The shock wave damage calculation module is used to perform air shock wave overpressure calculation, ground explosion shock wave overpressure calculation, and dynamic pressure calculation using the explosion equivalent, the explosion mode, and the distance from the explosion center to obtain the shock wave damage characteristic quantity; the shock wave calculation can be divided into two types: air shock wave and ground explosion shock wave. In specific implementation, the embodiment of the present application can provide that the air shock wave overpressure calculation is expressed by the following formula:

[0099]

[0100] The calculation of the ground explosion shock wave overpressure is expressed by the following formula:

[0101]

[0102] The dynamic pressure calculation is expressed by the following formula:

[0103]

[0104] The calculation of ground explosion shock wave overpressure is performed after the following judgment conditions are met:

[0105]

[0106] The stable cloud computing module is used to describe the stable cloud shape parameters and the sedimentation ash mass parameters, and is used to obtain the distribution of the stable cloud. In specific implementation, the embodiment of the present application can provide that the description of the stable cloud shape parameters includes:

[0107] According to the observation data, the height of the stable cloud cap is fitted, H B and top height H T It is expressed by the following formula:

[0108] H B =aQ b

[0109] H T =cQ d

[0110] Where Q represents the explosion equivalent, H B and H T The unit is meter, and its parameters are:

[0111] a=2228,b=0.3463; Q≤4.07kt

[0112] a=2661,b=0.2198; Q>4.07kt

[0113] c=3597,d=0.2553; Q<2.29kt

[0114] c=3170, d=0.4077; 2.29kt≤Q<19kt

[0115] c=6474,d=0.1650; Q≥19kt

[0116] Radius R of the stable cloud cap c The fitting formula is expressed as:

[0117] R C =exp[6.7553+0.32055lnQ+0.01137478(ln Q) 2 ].

[0118] The settling ash mass parameter includes the mass of soil plus shell fragments;

[0119] If it is determined to be an underground explosion, the mass parameter of the settled ash is expressed by the following formula:

[0120]

[0121] R=112.5+0.755d-9.6×10 -6 d 3 -9.11×10 -12 d 5

[0122] S=32.7+0.851d-2.52×10 -5 d 3 -1.78×10 -10 d 5

[0123] Where: d represents the specific depth;

[0124] If it is determined to be a ground explosion, the mass parameter of the settled ash is expressed by the following formula:

[0125]

[0126] Where: λ represents the specific height.

[0127] The explosion accident source term simulation and analysis method provided in the embodiments of this application, through comprehensive and complete analysis and simulation of large-yield explosions, can obtain source term parameters such as light radiation, shock waves, and smoke clouds. This allows decision makers to set corresponding protective measures based on the obtained source term parameters, providing a scientific basis for emergency response, pollution control, and public protection after the explosion accident. Compared with traditional single source term simulation, this method can obtain rich and diverse source term data with a single model, reducing model setup costs and computational time costs.

[0128] The following is a detailed introduction to the large-yield explosion accident source term simulation analysis method provided in the embodiments of the present application.

[0129] The embodiment of this application provides a complete explosion accident source term analysis model, which provides a scientific basis for emergency response, pollution control and public protection after the explosion accident. In specific implementation, a variety of different source terms can be simulated and analyzed.

[0130] (1) Light radiation (light impulse).

[0131] The main characteristic quantity of light radiation damage is the light impulse. Light impulse U refers to the total energy projected by the fireball onto a unit area of an object perpendicular to the direction of light propagation during the entire luminous time. Its international unit is J / m2.

[0132] Calculation of light impulse at close range:

[0133] The attenuation effect of the atmosphere on visible light and infrared rays is not obvious at close range, so when the distance from the explosion center is not too far, the calculation of the light impulse only considers the geometric attenuation factor. e It is the light radiation energy received per unit time on a unit area perpendicular to the radiation receiving surface at different distances. Since it indicates the intensity of the light and the degree to which the surface of the object is illuminated, the light impulse can be determined by integrating the illumination over the luminous time of the fireball, that is,

[0134]

[0135] Equivalent to 0 to t e The integral of the fireball's radiation power P divided by the surface area of the sphere with radius r during the time (light radiation action time) is:

[0136]

[0137] That is, the radiation energy E of the explosion accident G , so the light impulse can be expressed as follows:

[0138]

[0139] Light radiation equivalent Q G The radiation energy E of the explosion accident G The following relationship is satisfied:

[0140] Q G =E G / 4.19×10 12 (Formula 4)

[0141] Combining (3) and (4), E GIt can be expressed by the following formula:

[0142] E G =f G Q×4.19×10 12 (J) (Formula 5)

[0143] Therefore, when only the geometric attenuation is considered and the attenuation caused by atmospheric absorption and scattering is ignored, the optical impulse U can be calculated using the following formula 7:

[0144]

[0145] Among them, f G represents the light equivalent coefficient; Q represents the explosion accident equivalent; r represents the distance from the explosion center.

[0146] If it is an airburst, the light equivalent coefficient can be used to calculate the light impulse by referring to Table 1.

[0147] Table 1 Light equivalent coefficients of different equivalents during airburst

[0148]

[0149] If it is a ground explosion or surface explosion accident, f G f d Instead, the calculation formula for the light impulse of ground explosion accident is as follows:

[0150]

[0151] Where h B Indicates relative height, f d Indicates the effective light equivalent coefficient;

[0152] Calculation of light impulse taking into account atmospheric attenuation.

[0153] If the target is far from the explosion center and atmospheric visibility is poor, the calculation of the optical impulse must account for the attenuation of light radiation by atmospheric absorption and scattering. The transmission of optical radiation through the atmosphere is affected by a variety of factors. For example, varying atmospheric composition significantly attenuates certain wavelengths; uneven atmospheric density affects the propagation of optical radiation; and scattering attenuates optical radiation at close range but enhances it at greater distances, to the point where scattered light can dominate. In these cases, atmospheric transmittance must be considered when calculating the optical impulse.

[0154] Introducing the atmospheric transmittance τ, if it is an airburst, the calculation of the light impulse satisfies the following formula 9:

[0155]

[0156] If it is a ground explosion, the calculation of the light impulse satisfies the following formula 10:

[0157]

[0158] The atmospheric transmittance τ is the ratio of the transmitted light intensity to the incident light intensity after a parallel light beam passes through the atmosphere of a certain thickness. The expression formula is:

[0159]

[0160] In the formula: r is the distance (km) from the measurement point to the explosion center; μ is the average atmospheric attenuation coefficient (km -1 )

[0161] In practice, in most cases, the incident light is not a parallel light but has a certain solid angle. The empirical formula for the atmospheric transmittance summarized according to the test site is:

[0162] τ = e -μr +0.23(μr) 1.6 e -0.65μr (Equation 11) [[ID=?]] [[ID=?]]

[0163] In the above formula: e -μr is the transmittance of the direct light; 0.23(μr) 1.6 e -0.65μr is the enhancement term for multiple scattering. In the case of long distances, the latter plays a major role.

[0164] (1) When 0 < H ≤ 3 - Z0, the expression of the average atmospheric attenuation coefficient μ is:

[0165]

[0166] (2) When 3 - Z0 < H ≤ 30 - Z0, the expression of the average atmospheric attenuation coefficient μ is:

[0167]

[0168] In the formula: D is the ground horizontal atmospheric visibility (km); H is the explosion height (km); Z0 is the altitude of the projection point of the explosion center.

[0169]

[0170] Among them, the ground horizontal atmospheric visibility can be considered as the maximum distance (which needs to be input by the user or obtained from the system platform) to clearly see the outlines of distant gray targets such as hillocks and buildings against the background of the sky near the horizon with normal vision.

[0171] The influence of terrain and ground objects on light radiation: <000Any terrain or landform that can block light, such as hills, mountains, earth piles, houses, deep trenches, large craters, etc., can affect the propagation of light radiation and thus affect the distribution of light impulse. Figure 2 shown.

[0173] Calculation process of optical radiation damage characteristic quantity:

[0174] The calculation process of the characteristic quantity of light radiation damage, light impulse U, is as follows: Figure 3 shown.

[0175] (2) Shock wave.

[0176] When an explosion occurs in the atmosphere, a huge amount of energy is released instantly, forming a high-temperature and high-pressure fireball. As the fireball expands violently, it rapidly compresses the surrounding air and spreads at an extremely high speed, forming a shock wave.

[0177] The main characteristic quantities of shock wave damage are overpressure and dynamic pressure. Overpressure refers to the pressure exceeding atmospheric pressure, represented by Δp, with the unit of Pa. Dynamic pressure refers to the impact pressure generated by the high-speed airflow within the shock wave in its direction of motion, represented by q, with the unit of Pa.

[0178] Calculation of shock wave overpressure in air explosion accidents:

[0179] The calculation of shock wave overpressure usually adopts the uniform atmosphere point explosion theory. This theory assumes that after the shock wave of an explosion accident is formed, the volume and mass of the air it contains greatly exceed the volume and mass of the projectile. The projectile has no effect on the propagation of the shock wave. The explosion accident can be regarded as a point explosion with no mass and no volume, which releases huge energy in an instant.

[0180] If the explosion occurs under standard atmospheric conditions, that is, p0=1.01325×10 5 Pa, ρ0=1.225kg / m 3 , in an infinite uniform atmospheric medium with T0 = 288.16K, the calculation formula for the shock wave overpressure is:

[0181]

[0182] Where Q is the explosion equivalent, in kt; r is the distance from the explosion center, in m.

[0183] Calculation of overpressure of ground explosion shock wave:

[0184] When an explosion occurs near the ground or touches the ground, the shock wave is approximately hemispherical or semi-spherical. The air shock wave front propagating in the air layer near the ground is always perpendicular to the ground surface. The conditions for the shock wave to interact with the ground surface are fixed, and the energy consumed in the deformation of the ground soil usually only accounts for a few percent of the total explosion energy and can be ignored. In the volume enclosed by the shock wave, the energy density of the air is approximately twice that of an air explosion of the same equivalent. Therefore, in the case of a ground explosion, the shock wave overpressure is calculated as:

[0185]

[0186] Where Q is the explosion equivalent, in kt; r is the distance from the explosion center, in m.

[0187] Note that the condition for using formula (16) is

[0188] Dynamic pressure calculation:

[0189] Dynamic pressure is characterized by the kinetic energy of air per unit volume, namely:

[0190]

[0191] In formula (11), V i is the velocity of the shock wave airflow (m / s); ρ i is the air density on the shock wave front, which is related to the overpressure on the wave front and the undisturbed atmospheric pressure p0 (i.e. normal atmospheric pressure) as follows:

[0192]

[0193] Assuming that the explosion accident occurs under standard atmospheric conditions, ρ0, γ and V i Substituting the numerical values and expressions of the parameters into formula (18), the dynamic pressure calculation formula of the shock wave is obtained as follows:

[0194]

[0195] Calculation process of shock wave damage characteristic quantity:

[0196] The calculation process of shock wave damage characteristic quantities overpressure Δp and dynamic pressure q is as follows: Figure 4 shown.

[0197] (3) Stable cloud.

[0198] like Figure 5 As shown in Figure 2, the source term model is used to describe the calculation model of parameters such as the shape, particle distribution, and activity distribution of the "stable cloud".

[0199] Shape parameters:

[0200] The explosion cloud is assumed to be a flat ellipsoid. According to the observed data, the height of the fitted stable cloud cap is H B and top height H T It can be expressed as:

[0201] H B =aQ b (Equation 19)

[0202] H T =cQ d (Equation 20)

[0203] Where Q represents the explosive equivalent (kt), H B and H T The unit is meter, and its parameters are:

[0204] a=2228, b=0.3463; Q≤4.07kt (Formula 21)

[0205] a=2661, b=0.2198; Q>4.07kt (Formula 22)

[0206] c=3597,d=0.2553;Q<2.29kt (Equation 23)

[0207] c=3170,d=0.4077;2.29kt≤Q<19kt (Equation 24)

[0208] c=6474,d=0.1650;Q≥19kt (Equation 25) Radius R of the stable cloud cap c The fitting formula is:

[0209] R C =exp[6.7553+0.32055ln Q+0.01137478(ln Q) 2 ](Equation 26)

[0210] Where W represents the explosion energy equivalent (kt).

[0211] Settled ash mass:

[0212] The mass of fallen ash is the mass of soil plus weapon fragments.

[0213] Underground Explosion:

[0214]

[0215] Among them, Q is the explosion accident equivalent (kt),

[0216] R=112.5+0.755d-9.6×10 -6 d 3-9.11×10 -12 d 5 (Equation 28)

[0217] S=32.7+0.851d-2.52×10 -5 d 3 -1.78×10 -10 d 5 (Equation 29)

[0218] d is the relative depth

[0219] Ground explosion:

[0220]

[0221] Among them, λ is the specific height

[0222] Air or high altitude explosion: soil is not considered, and the shell fragments of the explosion accident are 90.7kg.

[0223] In summary, the explosion accident source term simulation and analysis method provided in this application has established a relatively complete explosion accident source term analysis model, which can simulate the distribution of light radiation (light impulse) and shock wave stable cloud after the explosion accident occurs, and provide a basis for evaluating and simulating the propagation, impact and prevention of explosion accidents after the occurrence of the event.

[0224] See also Figure 6 , the embodiment of the present application can also provide an explosion accident source term simulation analysis device, such as Figure 6 As shown, the device for executing the above-mentioned explosion accident source term simulation analysis method may include:

[0225] A basic parameter acquisition unit 601 is used to acquire basic parameter data of an explosion accident, wherein the basic parameter data includes at least explosion equivalent, explosion mode, explosion type, distance from the explosion center, and explosion height;

[0226] An analysis result acquisition unit 602 is configured to input the basic parameter data into an explosion accident source term analysis model to obtain an explosion accident source term simulation analysis result, wherein the explosion accident source term simulation analysis result includes at least a light radiation damage characteristic quantity, a shock wave damage characteristic quantity, and a distribution of stable clouds;

[0227] The explosion accident source analysis model includes at least a light impulse calculation module, a shock wave damage calculation module and a stable cloud computing module;

[0228] The light impulse calculation module is used to calculate the light impulse and the influence of terrain and objects on light radiation by using the explosion equivalent, the explosion mode and the distance from the explosion center to obtain the light radiation damage characteristic value;

[0229] The shock wave damage calculation module is used to perform air shock wave overpressure calculation, ground explosion shock wave overpressure calculation, and dynamic pressure calculation using the explosion equivalent, the explosion mode, and the distance from the explosion center to obtain the shock wave damage characteristic value;

[0230] The stable cloud computing module is used to describe the stable cloud shape parameters and the sedimentation ash mass parameters, and is used to obtain the distribution of the stable cloud.

[0231] The embodiment of the present application may further provide an explosion accident source term simulation analysis device, the device comprising a processor and a memory:

[0232] The memory is used to store program code and transmit the program code to the processor;

[0233] The processor is used to execute the steps of the above-mentioned explosion accident source term simulation analysis method according to the instructions in the program code.

[0234] like Figure 7 As shown, an explosion accident source term simulation and analysis device provided by an embodiment of the present application may include: a processor 10, a memory 11, a communication interface 12, and a communication bus 13. The processor 10, the memory 11, and the communication interface 12 all communicate with each other via the communication bus 13.

[0235] In the embodiment of the present application, the processor 10 may be a central processing unit (CPU), a graphics processing unit (GPU), an application-specific integrated circuit, a digital signal processor, a field programmable gate array, or other programmable logic devices.

[0236] The processor 10 may call a program stored in the memory 11 . Specifically, the processor 10 may execute operations in an embodiment of the explosion accident source term simulation analysis method.

[0237] The memory 11 is used to store one or more programs. The program may include program code, and the program code includes computer operating instructions. In the embodiment of the present application, the memory 11 stores at least a program for implementing the following functions:

[0238] Obtaining basic parameter data of the explosion accident, wherein the basic parameter data includes at least explosion equivalent, explosion mode, explosion type, distance from the explosion center, and explosion height;

[0239] Inputting the basic parameter data into an explosion accident source term analysis model to obtain explosion accident source term simulation analysis results, wherein the explosion accident source term simulation analysis results at least include light radiation damage characteristic quantities, shock wave damage characteristic quantities, and stable cloud distribution;

[0240] The explosion accident source analysis model includes at least a light impulse calculation module, a shock wave damage calculation module and a stable cloud computing module;

[0241] The light impulse calculation module is used to calculate the light impulse and the influence of terrain and objects on light radiation by using the explosion equivalent, the explosion mode and the distance from the explosion center to obtain the light radiation damage characteristic value;

[0242] The shock wave damage calculation module is used to perform air shock wave overpressure calculation, ground explosion shock wave overpressure calculation, and dynamic pressure calculation using the explosion equivalent, the explosion mode, and the distance from the explosion center to obtain the shock wave damage characteristic value;

[0243] The stable cloud computing module is used to describe the stable cloud shape parameters and the sedimentation ash mass parameters, and is used to obtain the distribution of the stable cloud.

[0244] In one possible implementation, the memory 11 may include a program storage area and a data storage area, wherein the program storage area can store an operating system and application programs required for at least one function (such as a file creation function, a data reading and writing function), etc.; the data storage area can store data created during use, such as initialization data, etc.

[0245] In addition, the memory 11 may include a high-speed random access memory and may also include a non-volatile memory, such as at least one disk storage device or other volatile solid-state storage device.

[0246] The communication interface 12 may be an interface of a communication module, and is used to connect to other devices or systems.

[0247] Of course, it needs to be explained that Figure 7 The structure shown does not constitute a limitation on the explosion accident source term simulation and analysis device in the embodiment of the present application. In actual application, the explosion accident source term simulation and analysis device may include Figure 7 More or fewer components than shown, or combinations of certain components.

[0248] An embodiment of the present application may also provide a computer-readable storage medium, wherein the computer-readable storage medium is used to store program code, and the program code is used to execute the steps of the above-mentioned explosion accident source term simulation analysis method.

[0249] It should be noted that, in this document, relational terms such as first and second, etc., are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply the existence of any such actual relationship or order between these entities or operations. Moreover, the terms "comprises," "comprising," or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article, or device comprising a series of elements includes not only those elements, but also other elements not explicitly listed, or elements inherent to such process, method, article, or device. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of other identical elements in the process, method, article, or device comprising the element.

[0250] Through the description of the above embodiments, it can be seen that those skilled in the art can clearly understand that the present application can be implemented by means of software plus a necessary general hardware platform. Based on this understanding, the technical solution of the present application, or the part that contributes to the prior art, can be embodied in the form of a software product, which can be stored in a storage medium such as ROM / RAM, a magnetic disk, an optical disk, etc., and includes a number of instructions for enabling a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in various embodiments or certain parts of the embodiments of the present application.

[0251] Each embodiment in this specification is described in a progressive manner. The same or similar parts between the embodiments can be referred to each other. Each embodiment focuses on the differences from other embodiments. In particular, for system or system embodiments, since they are basically similar to method embodiments, the description is relatively simple. For relevant parts, refer to the partial description of the method embodiment. The system and system embodiments described above are merely schematic, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment. A person of ordinary skill in the art can understand and implement it without expending creative work.

[0252] The above description is only a preferred embodiment of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention are included in the scope of protection of the present invention.

Claims

1. A method for simulating and analyzing explosion accident sources, characterized in that: include: Obtaining basic parameter data of the explosion accident, wherein the basic parameter data includes at least explosion equivalent, explosion mode, explosion type, distance from the explosion center, and explosion height; Inputting the basic parameter data into an explosion accident source term analysis model to obtain explosion accident source term simulation analysis results, wherein the explosion accident source term simulation analysis results at least include light radiation damage characteristic quantities, shock wave damage characteristic quantities, and stable cloud distribution; The explosion accident source analysis model includes at least a light impulse calculation module, a shock wave damage calculation module and a stable cloud computing module; The light impulse calculation module is used to calculate the light impulse and the influence of terrain and objects on light radiation by using the explosion equivalent, the explosion mode and the distance from the explosion center to obtain the light radiation damage characteristic value; The shock wave damage calculation module is used to perform air shock wave overpressure calculation, ground explosion shock wave overpressure calculation, and dynamic pressure calculation using the explosion equivalent, the explosion mode, and the distance from the explosion center to obtain the shock wave damage characteristic value; The stable cloud computing module is used to describe the stable cloud shape parameters and the sedimentation ash mass parameters, and is used to obtain the distribution of the stable cloud.

2. The explosion accident source term simulation analysis method according to claim 1, characterized in that: The light impulse calculation is used to perform the following operations: determining whether to consider attenuation caused by atmospheric absorption and scattering according to the distance from the explosion center; If only geometric attenuation is considered and the attenuation caused by atmospheric absorption and scattering is ignored, and the explosion mode is determined to be airburst, the light impulse is calculated using the following formula: Where: f G represents the light equivalent coefficient; Q represents the explosion equivalent; r represents the distance from the explosion center; If only geometric attenuation is considered and the attenuation caused by atmospheric absorption and scattering is ignored, and the explosion method is determined to be a ground explosion, the light impulse is calculated using the following formula: Where: h B Indicates relative height, f d Indicates the effective light equivalent coefficient; Considering the attenuation caused by atmospheric absorption and scattering, and determining that the explosion mode is airburst, the light impulse is calculated using the following formula: Where: τ represents the atmospheric transmittance; Considering the attenuation caused by atmospheric absorption and scattering, and determining that the explosion method is a ground explosion, the light impulse is calculated using the following formula:

3. The explosion accident source term simulation analysis method according to claim 2, characterized in that: The atmospheric transmittance τ is expressed by the following formula: τ=e -μr +0.23(μr) 1.6 in -0.65μr Where: e -μr Indicates the transmittance of direct light; 0.23 (μr) 1.6 e -0.65μr represents the enhancement term of multiple scattering, and μ represents the average atmospheric attenuation coefficient.

4. The explosion accident source term simulation analysis method according to claim 1, characterized in that: The airborne shock wave overpressure calculation is expressed by the following formula: The calculation of the ground explosion shock wave overpressure is expressed by the following formula: The dynamic pressure calculation is expressed by the following formula:

5. The explosion accident source term simulation analysis method according to claim 4, characterized in that: The calculation of ground explosion shock wave overpressure is performed after the following judgment conditions are met:

6. The explosion accident source term simulation analysis method according to claim 1, characterized in that: The parameters describing the stable cloud shape include: According to the observation data, the height of the stable cloud cap is fitted, H B and top height H T It is expressed by the following formula: H B =aQ b H T =cQ d Where Q represents the explosion equivalent, H B and H T The unit is meter, and its parameters are: a=2228,b=0.3463; Q≤4.07kt a=2661,b=0.2198; Q>4.07kt c=3597,d=0.2553; Q<2.29kt c=3170, d=0.4077; 2.29kt≤Q<19kt c=6474,d=0.1650; Q≥19kt Radius R of the stable cloud cap c The fitting formula is expressed as: R C =exp[6.7553+0.32055ln Q+0.01137478(lnQ) 2 ]。 7. The explosion accident source term simulation analysis method according to claim 1, characterized in that: The settling ash mass parameter includes the mass of soil plus shell fragments; If it is determined to be an underground explosion, the mass parameter of the settled ash is expressed by the following formula: R=112.5+0.755d-9.6×10 -6 d 3 -9.11×10 -12 d 5 S=32.7+0.851d-2.52×10 -5 d 3 -1.78×10 -10 d 5 Where: d represents the specific depth; If it is determined to be a ground explosion, the mass parameter of the settled ash is expressed by the following formula: Where: λ represents the specific height.

8. An explosion accident source term simulation and analysis device, characterized in that: The device is used to perform the explosion accident source term simulation analysis method according to any one of claims 1 to 7, comprising: A basic parameter acquisition unit is used to acquire basic parameter data of the explosion accident, wherein the basic parameter data at least includes explosion equivalent, explosion mode, explosion type, distance from the explosion center, and explosion height; an analysis result acquisition unit, configured to input the basic parameter data into an explosion accident source term analysis model to obtain an explosion accident source term simulation analysis result, wherein the explosion accident source term simulation analysis result at least includes a light radiation damage characteristic quantity, a shock wave damage characteristic quantity, and a distribution of stable clouds; The explosion accident source analysis model includes at least a light impulse calculation module, a shock wave damage calculation module and a stable cloud computing module; The light impulse calculation module is used to calculate the light impulse and the influence of terrain and objects on light radiation by using the explosion equivalent, the explosion mode and the distance from the explosion center to obtain the light radiation damage characteristic value; The shock wave damage calculation module is used to perform air shock wave overpressure calculation, ground explosion shock wave overpressure calculation, and dynamic pressure calculation using the explosion equivalent, the explosion mode, and the distance from the explosion center to obtain the shock wave damage characteristic value; The stable cloud computing module is used to describe the stable cloud shape parameters and the sedimentation ash mass parameters, and is used to obtain the distribution of the stable cloud.

Images