CN109580433B - Source term estimation method for diffusion of conventional explosive radioactive aerosol - Google Patents
Source term estimation method for diffusion of conventional explosive radioactive aerosol Download PDFInfo
- Publication number
- CN109580433B CN109580433B CN201811258555.0A CN201811258555A CN109580433B CN 109580433 B CN109580433 B CN 109580433B CN 201811258555 A CN201811258555 A CN 201811258555A CN 109580433 B CN109580433 B CN 109580433B
- Authority
- CN
- China
- Prior art keywords
- explosion
- explosive
- cloud
- aerosol
- particle size
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000002360 explosive Substances 0.000 title claims abstract description 81
- 239000000443 aerosol Substances 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000002285 radioactive effect Effects 0.000 title claims abstract description 30
- 238000009792 diffusion process Methods 0.000 title claims abstract description 27
- 238000004880 explosion Methods 0.000 claims abstract description 95
- 239000000779 smoke Substances 0.000 claims abstract description 75
- 239000002245 particle Substances 0.000 claims abstract description 64
- 230000000694 effects Effects 0.000 claims abstract description 22
- 239000007789 gas Substances 0.000 claims description 35
- 239000000941 radioactive substance Substances 0.000 claims description 16
- 239000000463 material Substances 0.000 claims description 14
- 238000005474 detonation Methods 0.000 claims description 13
- 239000000126 substance Substances 0.000 claims description 12
- 230000005484 gravity Effects 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000000295 complement effect Effects 0.000 claims description 3
- 229910052770 Uranium Inorganic materials 0.000 description 9
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 8
- 239000012857 radioactive material Substances 0.000 description 7
- 239000000428 dust Substances 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000002485 combustion reaction Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 239000011824 nuclear material Substances 0.000 description 4
- 230000005855 radiation Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 230000008021 deposition Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical class [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 1
- WZECUPJJEIXUKY-UHFFFAOYSA-N [O-2].[O-2].[O-2].[U+6] Chemical compound [O-2].[O-2].[O-2].[U+6] WZECUPJJEIXUKY-UHFFFAOYSA-N 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 125000004432 carbon atom Chemical group C* 0.000 description 1
- 229910002090 carbon oxide Inorganic materials 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000005202 decontamination Methods 0.000 description 1
- 230000003588 decontaminative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000037406 food intake Effects 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 230000008821 health effect Effects 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 229910017464 nitrogen compound Inorganic materials 0.000 description 1
- 150000002830 nitrogen compounds Chemical class 0.000 description 1
- 238000005025 nuclear technology Methods 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 241000894007 species Species 0.000 description 1
- XTQHKBHJIVJGKJ-UHFFFAOYSA-N sulfur monoxide Chemical class S=O XTQHKBHJIVJGKJ-UHFFFAOYSA-N 0.000 description 1
- 229910052815 sulfur oxide Inorganic materials 0.000 description 1
- 230000002459 sustained effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 229910000439 uranium oxide Inorganic materials 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/02—Investigating particle size or size distribution
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N2015/0023—Investigating dispersion of liquids
- G01N2015/0026—Investigating dispersion of liquids in gas, e.g. fog
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N2015/0042—Investigating dispersion of solids
- G01N2015/0046—Investigating dispersion of solids in gas, e.g. smoke
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measurement Of Radiation (AREA)
Abstract
The invention discloses a method for estimating the diffusion source item of a conventional explosive radioactive aerosol, which comprises the following steps: aiming at the conventional explosion of the radioactive texture surface, the explosion equivalent is estimated, and the size of the explosion smoke cloud is estimated based on the explosion equivalent, wherein the size of the explosion smoke cloud comprises the following steps: the height of the explosion smoke cloud, the radius of the explosion smoke cloud, the volume of the explosion smoke cloud and the density of the explosion smoke cloud are used for estimating the generation share of the explosion aerosol and the particle size distribution of the explosion aerosol, estimating the influence of heavy gas effect on source items and setting nuclides needed to be concerned about in an explosion event. The method is based on the explosion equivalent, the explosive mass, the local wind speed and the like, comprehensively considers the size of the explosion smoke cloud, the generation portion of the explosion aerosol, the aerosol particle size distribution and the influence of heavy gas effect on the source item, forms a source item estimation method, is more in line with the actual situation, and has the effect superior to a point source diffusion model.
Description
Technical Field
The invention relates to the field of radiation environment influence evaluation, in particular to a source item estimation method for diffusion of conventional explosive radioactive aerosol. Namely a calculation method of radioactive aerosol diffusion source items generated by conventional explosion.
Background
Since the 21 st century, especially after the 9.11 th incident, the unconventional security threats represented by various terrorist explosions have become more serious, and the nuclear materials and facilities have become one of the main targets of terrorists due to their particularity, and the security problems of the nuclear materials and facilities have become a major concern.
IAEA (international atomic energy agency) has divided nuclear and radiation terrorist attacks into 4 main categories:
1. illegally obtaining radioactive materials (radioactive materials), manufacturing Radioactive Distribution Devices (RDDs), and performing radioactive attacks, such as dirty bombs;
2. nuclear facilities such as an attack Nuclear Power Plant (NPP) and the like cause nuclear accidents, and radioactive substances are released into the environment;
3. illegally obtaining special nuclear materials, manufacturing a rough nuclear device (IND), and carrying out nuclear explosion;
4. and illegally acquiring a complete nuclear weapon and carrying out nuclear explosion.
Of these, the 2 nd, 3 rd and 4 th have great influence, but since it is difficult for terrorists to obtain sufficient raw nuclear weapon material, and the manufacturing technique is difficult, and the protection measures for nuclear facilities and nuclear weapons are strict, the probability of occurrence is very low, and the nuclear and radiation terrorist attacks most likely to be implemented by terrorists are the use of "dirty bombs" to make terrorist events.
With the continuous expansion of the application range of nuclear technology, radioactive materials are more widely applied in the fields of industry, agriculture, medicine, military and the like. Terrorists can obtain radioactive materials through various ways such as theft, smuggling, illegal purchase and the like. Although a "dirty bomb" terrorist attack has not occurred to date, the threat of a "dirty bomb" terrorist attack is real.
In the evaluation of radiation environment influence, the source item refers to the type, form and quantity (which can be called as source intensity) of the radioactive nuclide used by, carried by or generated by the source and the physicochemical morphological characteristics of the radioactive nuclide. Currently, analog calculation is mostly performed in a point source form in environmental impact evaluation, and a virtual point source form is also mostly adopted in an explosion event. The point sources have identifiable ranges that can be distinguished from other sources of contamination, simplifying the calculations in the form of one point in the mathematical model.
However, for an explosive event, whether initiated by accident or intentionally, a simple point source diffusion model is difficult to solve for atmospheric diffusion of explosive pollutants. The source terms required by the outcome calculation model of an explosion event (accident) are: after the explosion is over, a stable cloud (explosion cloud) is formed, which begins to diffuse under atmospheric action, with the instantaneous aerosol characteristics.
The size of the explosive smoke cloud and the source characteristics such as the distribution of internal substances affect the concentration of radioactive substances on the ground and in the air from the area near the explosive point to the area of several kilometers or more, and other consequences, so the problem of the explosive source needs to be considered in the evaluation of the consequences. The prior art lacks corresponding research and analysis, and has no corresponding forming technical scheme.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide a conventional explosive radioactive aerosol diffusion source item estimation method, and provides a source item calculation method for establishing an explosive event consequence evaluation model.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows:
a method for estimating the diffusion source term of a conventional explosive radioactive aerosol comprises the following steps:
aiming at the conventional explosion of the radioactive texture surface, the explosion equivalent is estimated,
estimating the size of the explosion cloud based on the explosion equivalent, the size of the explosion cloud comprising: height of the detonation smoke cloud, radius of the detonation smoke cloud, volume of the detonation smoke cloud, and density of the detonation smoke cloud, wherein:
when the height of the explosion smoke cloud is estimated, the explosion equivalent and the local wind speed are integrated,
the explosive equivalent is taken into account when estimating the radius of the explosive cloud,
the explosive quality and the volume of the explosive cloud are integrated when the density of the explosive cloud is estimated,
the explosive aerosol generation fraction and the aerosol particle size distribution were estimated,
the effect of heavy gas effects on the source terms is estimated and the nuclides of interest needed to set the explosive event.
Further, as mentioned above, the method for estimating the source term of the conventional explosive radioactive aerosol diffusion means to convert the explosive energy into TNT equivalent with a conversion relationship of 4.19MJ (megajoule) 1 kgTNT.
Further, the method for estimating the source term of the diffusion of the conventional explosive radioactive aerosol, which estimates the size of the explosive smoke cloud based on the explosive equivalent weight specifically includes:
the height of the explosion smoke cloud is estimated by the following formula,
h(w,t,v)=(6.3±1.0)×w0.29±0.03×t0.5±0.05stability or neutrality
h(w,t,v)=(6.3±1.0)×w0.29±0.013×t(0.7±0.013)-(0.03±0.002)vInstability of the film
Wherein:
h is the effective height that the smoke cloud can reach, and the unit is m,
w is the explosive TNT equivalent, the applicable range of the explosive equivalent is between 0.25kg and 100kg,
t is the time to reach the effective altitude, in units of s,
v is the local wind speed, in m/s,
the radius of the explosion smoke cloud is estimated by the following formula,
R(w)=3.5w0.375
the volume of the explosion smoke cloud is estimated by the following formula,
V=πR2h
the density of the explosive cloud is estimated by the following formula,
ρ=ρa+(Wd+WE)/V
wherein:
rho is the density of the explosive smoke cloud and has the unit of kg/m3,
ρaIs the air density in kg/m3,
WdIs the mass of the explosive, and the unit is kg,
WEis the total mass in kg of other materials contained in the explosive material.
Further, the method for estimating the source term of the diffusion of the conventional explosive radioactive aerosol includes:
calculating the mass fraction of the airborne radioactive substance, using the following formula,
F=2.378(WT/WE)-0.6383 WT≥5WE
F=1.0 WT<5WE
wherein:
WTis the TNT equivalent, in kg,
WEis contained in the exploderThe total mass of the other materials in the fry mass, in kg,
the particle size distribution of the airborne radioactive substance is usually lognormal distributed with the mass, the particle size distribution of the airborne radioactive substance is calculated, the following formula is adopted,
wherein:
FRis the particle diameter d1And d2The mass fraction of the gas-carrying substance in between,
d1and d2Respectively two particle sizes, d1>d2The unit is a number of microns,
m is a mass median diameter in μm, obtained experimentally,
g is the geometric standard deviation, an empirical value,
erfc is a complementary error function.
Further, the conventional method for estimating the source term of the explosive radioactive aerosol diffusion as described above is divided into two groups according to the particle size of more than 15 μm and the particle size of 15 μm or less when calculating the particle size distribution of the airborne radioactive material,
the part with the particle size of more than 15 μm and less than 100 μm is calculated according to the particles and divided into four groups according to the ground reflection coefficient: d is more than 15 and less than or equal to 30 mu m, d is more than 30 and less than or equal to 47 mu m, d is more than 47 and less than or equal to 75 mu m, d is more than 75 and less than or equal to 100 mu m,
fractions having a particle size of 15 μm or less, considered as gases, were divided into three groups: d is more than 0 and less than or equal to 1 mu m, d is more than 1 and less than or equal to 10 mu m, d is more than 10 and less than or equal to 15 mu m,
more than 100 μm are independently used as a group,
the fractions within the different fractions are determined from this particle size fraction.
Further, as described above, the method for estimating the source term of diffusion of a conventional explosive radioactive aerosol, where the method for estimating the influence of Heavy Gas effects on the source term adopts an IIT Heavy Gas Model-I instantaneous release Heavy Gas Model, specifically including:
gravity settling, using the following formula:
wherein:
dR/dt is the change rate of the radius of the smoke cloud along with time, and the unit is m/s,
ρais the density of air, in kg/m3,
K is the gravity settling coefficient of the smoke cloud,
Δ ρ is the density difference between the cloud and air in kg/m3,
Air entrainment, using the following formula:
wherein:
α*to control the parameters of the side draw rate,
ue is the velocity of the top entrainment, which is the Raphson number Ri and the radial turbulence velocity UlAs a function of (c).
Further, as described above in a conventional method for estimating the source term of explosive radioactive aerosol diffusion, the source term estimation is started after a stable smoke cloud is formed after explosion, and the temperature of the smoke cloud and the temperature of the surrounding atmosphere are considered to be equal.
The invention has the beneficial effects that: based on the explosion equivalent, the explosive mass, the local wind speed and the like, the size of the explosion smoke cloud, the generation share of the explosion aerosol, the aerosol particle size distribution and the influence of heavy gas effect on the source item are comprehensively considered to form a source item estimation method, which is more in line with the actual situation and has the effect superior to a point source diffusion model.
Drawings
FIG. 1 is a graph of the effect of heavy gas on the height and radius of the source term in accordance with the present invention. Wherein:
a is a graph of the height h over time.
And b is a graph of the change of the radius R with time.
FIG. 2 is a block diagram of a method of the present invention.
Detailed Description
The invention is described in further detail below with reference to the drawings and the detailed description.
Fig. 2 is a flowchart illustrating a method for estimating a source term of a conventional explosive radioactive aerosol diffusion according to an embodiment of the present invention, where the method mainly includes:
aiming at the conventional explosion of the radioactive texture surface, the explosion equivalent is estimated,
estimating the size of the explosion cloud based on the explosion equivalent, the size of the explosion cloud comprising: height of the detonation smoke cloud, radius of the detonation smoke cloud, volume of the detonation smoke cloud, and density of the detonation smoke cloud, wherein:
when the height of the explosion smoke cloud is estimated, the explosion equivalent and the local wind speed are integrated,
the explosive equivalent is taken into account when estimating the radius of the explosive cloud,
the explosive quality and the volume of the explosive cloud are integrated when the density of the explosive cloud is estimated,
the explosive aerosol generation fraction and the aerosol particle size distribution were estimated,
the effect of heavy gas effects on the source term is estimated and the nuclear species of interest needed to set the explosive event ("dirty bomb" explosive event).
The conventional explosion of the radioactive texture surface refers to the conventional explosion of the radioactive texture surface such as 'dirty bomb' explosion and the like.
Further, the estimation of the explosion equivalent refers to converting the energy of explosion into the TNT equivalent, and the conversion relation is 4.19MJ (megajoule) ═ 1 kgTNT.
The explosives adopted by conventional explosion can be different, and for convenience of calculation, in the dirty bomb explosion scene, the explosive explosion energy is converted into TNT equivalent.
Further, estimating the size of the explosion smoke cloud based on the explosion equivalent specifically comprises:
the height of the explosion smoke cloud is estimated by the following formula,
h(w,t,v)=(6.3±1.0)×w0.29±0.03×t0.5±0.05stability or neutrality
h(w,t,v)=(6.3±1.0)×w0.29±0.013×t(0.7±0.013)-(0.03±0.002)vInstability of the film
Wherein:
h is the effective height that the smoke cloud can reach, and the unit is m,
w is the explosive TNT equivalent, the applicable range of the explosive equivalent is between 0.25kg and 100kg,
t is the time to reach the effective altitude, in units of s,
v is the local wind speed, in m/s,
the formula for t is a function of the explosive equivalent, using the following formula,
t(w)=(28.0±1.1)w(0.33±0.4),
the radius of the explosion smoke cloud is estimated by the following formula,
R(w)=3.5w0.375
as for the radius of the smoke cloud formed after explosion, because the radius is slightly influenced by the environment and is greatly influenced by the explosion equivalent,
the volume of the explosion smoke cloud is estimated by the following formula,
V=πR2h
the density of the explosive cloud is estimated by the following formula,
ρ=ρa+(Wd+WE)/V
wherein:
rho is the density of the explosive smoke cloud and has the unit of kg/m3,
ρaIs the air density in kg/m3,
WdIs the mass of the explosive, and the unit is kg,
WEis the total mass in kg of other materials contained in the explosive material.
The process of 'dirty bomb' explosion is mainly divided into three stages:
stage one: rapid diffusion of explosive substances generated in an explosion. This process heats the air near the explosion point. Typically lasting from a few tenths of a millisecond to a second, is the fastest and least affected by meteorological conditions for the entire process.
And a second stage: the explosive material, hot air and entrained dust create expansion and elevation of the smoke cloud. At this stage, the average temperature within the cloud is significantly higher than the ambient temperature. The size and vertical velocity of the smoke cloud are determined by the average density within the smoke cloud, the ambient environmental conditions, the viscosity and velocity of the air, and the temperature distribution within the smoke cloud. If the explosive volume is relatively large, the heat in the cloud may continue to increase because of the presence of products of sustained combustion after the explosion. At this stage, the effective process of cooling the smoke cloud is radioactive, transformation of the smoke cloud, heat transfer to the environment, and expansion of the smoke cloud. This process lasts from a few seconds to tens of seconds.
And a third stage: the cloud continues to expand and rise. At this stage, the average temperature of the cloud is close to that of the surrounding environment, the cloud continues to increase the inertial force of the internal gas, and the expansion of the cloud is mainly due to diffusion and turbulence. This phase lasts for hundreds of seconds in a particularly large explosion and ends when it reaches the effective altitude.
When studying the size of the smoke cloud at the beginning of the explosion, the assumption is made that the heat rise of the smoke cloud after the explosion is rapid, before most diffusion, i.e. expansion precedes diffusion.
Further, the estimating of the generation fraction of the explosive aerosol and the aerosol particle size distribution specifically includes:
calculating the mass fraction of the airborne radioactive substance, using the following formula,
F=2.378(WT/WE)-0.6383 WT≥5WE
F=1.0 WT<5WE
wherein:
WTis the TNT equivalent, in kg,
WEis contained in an explosive substanceThe total mass of the material is kg,
the particle size distribution of the airborne radioactive substance is usually lognormal distributed with the mass, the particle size distribution of the airborne radioactive substance is calculated, the following formula is adopted,
wherein:
FRis the particle diameter d1And d2The mass fraction of the gas-carrying substance in between,
d1and d2Respectively two particle sizes, d1>d2The unit is a number of microns,
m is the mass median diameter, the value of m is 5, the unit is mum, and the m is obtained according to domestic experiments,
g is the geometric standard deviation, the value of g is 4, which is an empirical value,
erfc is a complementary error function.
The particle size of the aerosol particles affects the diffusion process after the formation of the explosive cloud, the deposition rate of the particles and the gas in the atmospheric environment is different, and the aerosol can generate some characteristics of heavy gas due to the existence of the particles. Evaluating the post-detonation particle scale problem involves two calculations:
(1) calculating the mass fraction of the substance that becomes airborne radioactive due to the explosion;
(2) the mass fraction of all airborne substances in the particle size range which can have a detrimental health effect on humans is calculated.
When calculating the particle size distribution (aerosol particle size distribution) of the airborne radioactive substance, the invention divides the aerosol into two groups and seven groups according to the particle size,
the two groups refer to the particle size of more than 15 μm and the particle size of less than or equal to 15 μm.
The part with the grain diameter more than 15 μm and less than 100 μm is calculated according to the particles and divided into four groups according to the ground reflection coefficient,
the fraction having a particle size of 15 μm or less, considered as a gas, aerosol particles having a particle size of 15 μm or less are divided into three groups,
the seven particle size groups are specifically: d is more than 0 and less than or equal to 1 mu m, d is more than or equal to 1 and less than or equal to 10 mu m, d is more than or equal to 10 and less than or equal to 15 mu m, d is more than 15 and less than or equal to 30 mu m, d is more than 30 and less than or equal to 47 mu m, d is more than 47 and less than or equal to 75 mu m, and d is more than 75 and less than or equal. The particles having a particle size of 100 μm or more (excluding the number) are individually grouped. The invention determines the fraction within the different fractions on the basis of this particle size fraction.
It should be noted that: according to the summary of the relevant documents of the aerosol, the particles with the smallest aerosol deposition rate are about 1 μm in particle size, some books are different from guiding rules, and the aerosol with the particle size below 10 μm is regarded as gas, so the invention takes the two points as the basis to group the aerosols with the particle size of 15 μm or less.
Examples are as follows:
a uranium mass of 4kg and a TNT equivalent weight of 0.5kg were selected and the calculated particle size groups are given in Table 1.
TABLE 1 particle size groups and fraction of each group
Further, the estimating of the influence of the Heavy Gas effect on the source term adopts an IIT Heavy Gas Model-I instantaneous Heavy Gas release Model, which specifically comprises:
gravity settling, using the following formula:
wherein:
dR/dt is the change rate of the radius of the smoke cloud along with time, and the unit is m/s,
ρais the density of air, in kg/m3,
K is the gravity settling coefficient of the smoke cloud, the value is 1,
Δ ρ is the density difference between the cloud and air in kg/m3,
Air entrainment, using the following formula:
wherein:
α*to control the parameters of the side draw rate,
ue is the velocity of the top entrainment, which is the Raphson number Ri and the radial turbulence velocity UlAs a function of (a) or (b),
cloud cluster heating adopts the following formula:
wherein:
t is the temperature of the smoke cloud in units of,
dT/dT is the amount of change in the smoke cloud temperature over time,
delta Ta is the difference between the temperature of the smoke cloud and the temperature of the air in units of,
cp is the specific heat capacity of the released substance, in J/(kg ℃ C.),
m is the mass of the released substance,
qc is the rate of heat exchange of the smoke cloud due to turbulent and forced exchanges, depending on the temperature difference between the ground and the smoke cloud and various thermodynamic properties of the material,
cpa is the specific heat capacity of the entrained air,
ma is the mass of air entrained.
The aerosol formed after the explosion of the dirty bomb mainly comprises products after the explosion and the combustion of the explosive, products after the combustion or the oxidation of radioactive substances, dust and the like which are absorbed into the smoke cloud. The gas products after explosive combustion are mainly carbon oxides, sulfur oxides, nitrogen compounds and the like, and the molecular weights of the compounds are not much different from the molecular weight of air and are gases; the solid particles such as radioactive aerosol particles and dust are not continuous gas like gas, and the distance between the solid particles is far larger than the distance between air molecules, so that the solid particles can be used for preparing the gasConsidering that the density of the aerosol is the sum of the density of air and the density of radioactive substances and dust, so that the density of the aerosol is greater than that of air, it can be regarded as heavy gas, and the heavy gas effect should be considered to affect the explosive source item of the dirty bomb. The standard for judging heavy gas effect is density difference<0.0001kg/m3And also as a criterion for the loss of the heavy gas effect.
In the prior art, three equations of a box model are used for heavy gas effect judgment, and are correlated with each other, so that calculation is troublesome.
The invention adopts an IIT Heavy Gas Model-I instantaneous release Heavy Gas Model, the side entrainment is considered to enter, the entrainment speed Ue depends on UlAnd Ri formula as follows:
Ue=α′UlRi-1
in the above formula, the first and second carbon atoms are,
wherein:
α' is the entrainment constant, the IIT Heavy Gas Model-I Model recommends a value of 0.8,
Ulis that the radial turbulence rate is proportional to the friction speed U, UlParameter a is related to atmospheric stability, as shown in table 2:
TABLE 2 stability vs. parameter a
g is the gravity acceleration, 9.8m/s is taken,
ls is the turbulence length, which is a function of height, and the expression for ls is given by Taylor (1970) as follows:
ls=5.88h0.48。
as mentioned above, the density of the aerosol is the sum of the air density and the density of the radioactive material and the dust, and thus the density of the plume after entrainment of air can be considered as the ratio of the sum of the mass of the entrained air and the mass of the initial plume to the volume after entrainment of air, i.e.:
wherein:
a is the initial mass of the smoke cloud,
and M is the increased volume of entrained air multiplied by the density of the air, M ═ R2h-400). times.1.293, and obtaining the following by taking the derivatives of the two sides of the formula:
as can be seen from the foregoing, the present invention,
thus, it is possible to obtain
The parameters calculated in the previous step are brought in and simplified to obtain
While
And substituting the parameters to obtain the result of R, and solving h by using MATLAB.
In the 'dirty bomb' explosion scene, the source term estimation is started after stable smoke cloud is formed after explosion, the temperature of the smoke cloud is equal to the temperature of the surrounding atmosphere, and therefore the radius R of the explosion smoke cloud and the effective height h which the smoke cloud can reach are obtained without considering cloud cluster heating. For example, fig. 1 is a graph showing the results of the effects of heavy gas on the spread of smoke cloud calculated when the mass of uranium is 4kg and the TNT equivalent amount is 0.5 kg. Wherein: FIG. 1.a is a graph of the height h as a function of time. FIG. 1.b is a graph of the change of radius R with time.
Further, the nuclide needed to be concerned by the set explosion event is selected according to the following principle:
(1) ease of acquisition of radioactive materials;
(2) a greater number of possibilities can be obtained;
(3) the use and the manufacture are simple;
(4) easy shielding;
(5) higher activity levels;
(6) a sufficiently long half-life.
The method specifically comprises the following steps:
60Co、137Cs、192Ir、90Sr、241Am、252Cf、238Pu。
currently, the most commonly used isotopes in the medical, scientific and industrial industries are likely to be the commercial sources for "dirty bombs" and are set as nuclides of interest.
The properties of the isotopes most commonly used in the medical, scientific and industrial industries are listed in table 3.
TABLE 3 commercial sources that may be used for "dirty bombs
In addition, uranium also requires attention in addition to the several nuclides mentioned above because:
(1) the application is wide, the uranium is easy to obtain, the domestic and foreign control of the uranium is relaxed, and the application of the uranium in the civil aspect is also very wide.
(2) Emitters, which are primarily alpha-rays, low-energy beta-rays and gamma-rays, are easier to shield than Co and Cs, and are also convenient for terrorists to operate.
(3) Most of the detectors at present are used for detecting gamma rays, and the gamma ray detector is difficult to detect alpha rays, which causes certain difficulty for detecting uranium oxide after the explosion of a dirty bomb.
(4) Inhalation and ingestion of uranium can cause chemical and radiological hazards to the human body.
(5) Uranium is widely known as the primary material for making nuclear weapons, and if the material of a "dirty bomb" is uranium, it may cause greater panic than other nuclear materials not known.
The invention has the beneficial technical effects that: on the basis of basic conditions such as explosion equivalent, explosive mass, local wind speed and the like, the influences of the size of explosion smoke cloud, the generation portion of explosion aerosol, aerosol particle size distribution and heavy gas effect on source items are fully considered, and nuclides needing attention are listed. Compared with a point source, the source item estimation method provided by the invention is more suitable for the situation of 'dirty bomb' explosion in practice.
The invention can provide a preliminary basis for emergency plan and decontamination work of the place where the dirty bomb explosion event occurs, provides a basis for calculating the diffusion of radioactive substances in city streets, and provides a material and a theoretical basis for further research on the dirty bomb terrorist attack event in the future.
The present invention is not limited to the above specific embodiments, and those skilled in the art can derive other embodiments according to the technical solutions of the present invention, which also belong to the technical innovation scope of the present invention.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is intended to include such modifications and variations.
Claims (3)
1.A method for estimating the diffusion source term of a conventional explosive radioactive aerosol comprises the following steps:
aiming at the conventional explosion of the radioactive texture surface, the explosion equivalent is estimated,
estimating the size of the explosion cloud based on the explosion equivalent, the size of the explosion cloud comprising: height of the detonation smoke cloud, radius of the detonation smoke cloud, volume of the detonation smoke cloud, and density of the detonation smoke cloud, wherein:
when the height of the explosion smoke cloud is estimated, the explosion equivalent and the local wind speed are integrated,
the explosive equivalent is taken into account when estimating the radius of the explosive cloud,
the explosive quality and the volume of the explosive cloud are integrated when the density of the explosive cloud is estimated,
the explosive aerosol generation fraction and the aerosol particle size distribution were estimated,
estimating the influence of heavy gas effect on source items and setting the nuclide of interest required by the explosion event;
the estimation of the explosion equivalent refers to converting the explosion energy into TNT equivalent, and the conversion relation is 4.19MJ (megajoule) which is 1 kgTNT;
the estimating of the size of the explosion smoke cloud based on the explosion equivalent specifically comprises:
the height of the explosion smoke cloud is estimated by the following formula,
h(w,t,v)=(6.3±1.0)×w0.29±0.03×t0.5±0.05stability or neutrality
h(w,t,v)=(6.3±1.0)×w0.29±0.013×t(0.7±0.013)-(0.03±0.002)vInstability of the film
Wherein:
h is the effective height that the smoke cloud can reach, and the unit is m,
w is the explosive TNT equivalent, the applicable range of the explosive equivalent is between 0.25kg and 100kg,
t is the time to reach the effective altitude, in units of s,
v is the local wind speed, in m/s,
the radius of the explosion smoke cloud is estimated by the following formula,
R(w)=3.5w0.375
the volume of the explosion smoke cloud is estimated by the following formula,
V=πR2h
the density of the explosive cloud is estimated by the following formula,
ρ=ρa+(Wd+WE)/V
wherein:
rho is the density of the explosive smoke cloud and has the unit of kg/m3,
ρaIs the air density in kg/m3,
WdIs the mass of the explosive, and the unit is kg,
WEis the total mass of other materials contained in the explosive substance in kg;
the estimating of the generation portion of the explosive aerosol and the aerosol particle size distribution specifically comprises the following steps:
calculating the mass fraction of the airborne radioactive substance, using the following formula,
F=2.378(WT/WE)-0.6383 WT≥5WE
F=1.0 WT<5WE
wherein:
WTis the TNT equivalent, in kg,
WEis the total mass of other materials contained in the explosive substance, in kg,
the particle size distribution of the airborne radioactive substance is usually lognormal distributed with the mass, the particle size distribution of the airborne radioactive substance is calculated, the following formula is adopted,
wherein:
FRis the particle diameter d1And d2The mass fraction of the gas-carrying substance in between,
d1and d2Respectively two particle sizes, d1>d2The unit is a number of microns,
m is a mass median diameter in μm, obtained experimentally,
g is the geometric standard deviation, an empirical value,
erfc is a complementary error function;
the method for estimating the influence of the Heavy Gas effect on the source item adopts an IIT Heavy Gas Model-I instantaneous release Heavy Gas Model, and specifically comprises the following steps:
gravity settling, using the following formula:
wherein:
dR/dt is the change rate of the radius of the smoke cloud along with time, and the unit is m/s,
ρais the density of air, in kg/m3,
K is the gravity settling coefficient of the smoke cloud,
Δ ρ is the density difference between the cloud and air in kg/m3,
Air entrainment, using the following formula:
wherein:
α*to control the parameters of the side draw rate,
ue is the rate of top entrainment, which is the Raphson number Ri and the radial turbulenceSpeed UlAs a function of (c).
2. The method for estimating the source term of a conventional explosive radioactive aerosol diffusion according to claim 1, wherein: when calculating the particle size distribution of the airborne radioactive substance, the method is divided into two groups according to the particle size of more than 15 μm and the particle size of less than or equal to 15 μm,
the part with the particle size of more than 15 μm and less than 100 μm is calculated according to the particles and divided into four groups according to the ground reflection coefficient: d is more than 15 and less than or equal to 30 mu m, d is more than 30 and less than or equal to 47 mu m, d is more than 47 and less than or equal to 75 mu m, d is more than 75 and less than or equal to 100 mu m,
fractions having a particle size of 15 μm or less, considered as gases, were divided into three groups: d is more than 0 and less than or equal to 1 mu m, d is more than 1 and less than or equal to 10 mu m, d is more than 10 and less than or equal to 15 mu m,
more than 100 μm are independently used as a group,
the fractions within the different fractions are determined from this particle size fraction.
3. The method for estimating the source term of a conventional explosive radioactive aerosol diffusion according to claim 1, wherein: the source term estimation begins after a stable cloud of smoke is formed after the explosion, at which point the temperature of the cloud of smoke and the temperature of the surrounding atmosphere are considered to be equal.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811258555.0A CN109580433B (en) | 2018-10-26 | 2018-10-26 | Source term estimation method for diffusion of conventional explosive radioactive aerosol |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811258555.0A CN109580433B (en) | 2018-10-26 | 2018-10-26 | Source term estimation method for diffusion of conventional explosive radioactive aerosol |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109580433A CN109580433A (en) | 2019-04-05 |
CN109580433B true CN109580433B (en) | 2021-05-28 |
Family
ID=65920979
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811258555.0A Active CN109580433B (en) | 2018-10-26 | 2018-10-26 | Source term estimation method for diffusion of conventional explosive radioactive aerosol |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109580433B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109387544B (en) * | 2018-09-03 | 2021-01-15 | 中国辐射防护研究院 | Method for estimating explosion source item of hydrogen mixed gas in high-level radioactive waste liquid storage tank |
CN112651595B (en) * | 2020-12-01 | 2022-05-17 | 中国辐射防护研究院 | Method for formulating emergency action level of post-processing plant |
CN112781942A (en) * | 2020-12-16 | 2021-05-11 | 中国辐射防护研究院 | Explosive aerosol source item acquisition and environment monitoring method |
Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005327214A (en) * | 2004-05-17 | 2005-11-24 | Mitsui Sumitomo Insurance Co Ltd | Disaster occurrence frequency estimation device, method, software and recording medium, and disaster event occurrence frequency estimation system |
CN101201311A (en) * | 2006-12-12 | 2008-06-18 | 中国科学院过程工程研究所 | Method for measuring heterogeneous texture parameter distribution of particle fluid two-phase flow system |
CN102567808B (en) * | 2010-12-31 | 2015-04-29 | 北京工业大学 | Method for forecasting and warning accident consequence of major hazard installation by combining with real-time meteorological information |
CN104933294A (en) * | 2015-05-25 | 2015-09-23 | 苏州热工研究院有限公司 | Method for constructing migration and propagation evaluation model of radioactive substance |
CN105138715A (en) * | 2015-06-29 | 2015-12-09 | 中国人民解放军军事医学科学院生物工程研究所 | Microbial aerosol atmospheric diffusion hazard assessment method |
CN105228761A (en) * | 2013-02-11 | 2016-01-06 | 冲击及清洁装置有限责任公司 | The method and apparatus of the inner space of clean container and equipment |
CN205164420U (en) * | 2015-12-04 | 2016-04-20 | 山东理工大学 | Aerosol remove device |
CN105549065A (en) * | 2015-12-09 | 2016-05-04 | 中国船舶重工集团公司第七一九研究所 | Artificial radioactive aerosol source machine |
CN107145699A (en) * | 2016-03-01 | 2017-09-08 | 中国辐射防护研究院 | Airborne radionuclide long-distance migration Lagrangian Particle Dispersion computational methods |
CN107145698A (en) * | 2016-03-01 | 2017-09-08 | 中国辐射防护研究院 | Core and radiation accident consequence large scale analogy method and simulation system |
CN107145613A (en) * | 2016-03-01 | 2017-09-08 | 中国辐射防护研究院 | Nuclear fission blast source item evaluation method based on TNT equivalents |
CN107391845A (en) * | 2017-07-23 | 2017-11-24 | 南京理工大学 | A kind of method for numerical simulation in tobacco edulcoration device flow field |
CN107526852A (en) * | 2016-06-21 | 2017-12-29 | 中国辐射防护研究院 | A kind of nuclear facilities Accident Off-site Consequence real-time online evaluation method and system |
CN108020448A (en) * | 2016-11-03 | 2018-05-11 | 中国辐射防护研究院 | A kind of method that polonium purifies source processed in aerosol |
CN108563840A (en) * | 2018-03-23 | 2018-09-21 | 西安交通大学 | A kind of nuclear reactor vapour explosion comprehensive analysis method |
CN109633548A (en) * | 2018-12-10 | 2019-04-16 | 禁核试北京国家数据中心 | A kind of underwater sound platform gateway linked method |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
PL3134403T3 (en) * | 2014-04-23 | 2020-07-13 | Incyte Corporation | 1h-pyrrolo[2,3-c]pyridin-7(6h)-ones and pyrazolo[3,4-c]pyridin-7(6h)-ones as inhibitors of bet proteins |
US20170316080A1 (en) * | 2016-04-29 | 2017-11-02 | Quest Software Inc. | Automatically generated employee profiles |
-
2018
- 2018-10-26 CN CN201811258555.0A patent/CN109580433B/en active Active
Patent Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005327214A (en) * | 2004-05-17 | 2005-11-24 | Mitsui Sumitomo Insurance Co Ltd | Disaster occurrence frequency estimation device, method, software and recording medium, and disaster event occurrence frequency estimation system |
CN101201311A (en) * | 2006-12-12 | 2008-06-18 | 中国科学院过程工程研究所 | Method for measuring heterogeneous texture parameter distribution of particle fluid two-phase flow system |
CN102567808B (en) * | 2010-12-31 | 2015-04-29 | 北京工业大学 | Method for forecasting and warning accident consequence of major hazard installation by combining with real-time meteorological information |
CN105228761A (en) * | 2013-02-11 | 2016-01-06 | 冲击及清洁装置有限责任公司 | The method and apparatus of the inner space of clean container and equipment |
CN104933294A (en) * | 2015-05-25 | 2015-09-23 | 苏州热工研究院有限公司 | Method for constructing migration and propagation evaluation model of radioactive substance |
CN105138715A (en) * | 2015-06-29 | 2015-12-09 | 中国人民解放军军事医学科学院生物工程研究所 | Microbial aerosol atmospheric diffusion hazard assessment method |
CN205164420U (en) * | 2015-12-04 | 2016-04-20 | 山东理工大学 | Aerosol remove device |
CN105549065A (en) * | 2015-12-09 | 2016-05-04 | 中国船舶重工集团公司第七一九研究所 | Artificial radioactive aerosol source machine |
CN107145699A (en) * | 2016-03-01 | 2017-09-08 | 中国辐射防护研究院 | Airborne radionuclide long-distance migration Lagrangian Particle Dispersion computational methods |
CN107145698A (en) * | 2016-03-01 | 2017-09-08 | 中国辐射防护研究院 | Core and radiation accident consequence large scale analogy method and simulation system |
CN107145613A (en) * | 2016-03-01 | 2017-09-08 | 中国辐射防护研究院 | Nuclear fission blast source item evaluation method based on TNT equivalents |
CN107526852A (en) * | 2016-06-21 | 2017-12-29 | 中国辐射防护研究院 | A kind of nuclear facilities Accident Off-site Consequence real-time online evaluation method and system |
CN108020448A (en) * | 2016-11-03 | 2018-05-11 | 中国辐射防护研究院 | A kind of method that polonium purifies source processed in aerosol |
CN107391845A (en) * | 2017-07-23 | 2017-11-24 | 南京理工大学 | A kind of method for numerical simulation in tobacco edulcoration device flow field |
CN108563840A (en) * | 2018-03-23 | 2018-09-21 | 西安交通大学 | A kind of nuclear reactor vapour explosion comprehensive analysis method |
CN109633548A (en) * | 2018-12-10 | 2019-04-16 | 禁核试北京国家数据中心 | A kind of underwater sound platform gateway linked method |
Non-Patent Citations (4)
Title |
---|
A review of very large vapour cloud explosions:cloud formation and explosion severity;G Atkinson 等;《Journal of loss prevention in the process industries》;20170731;第48卷;第367-375页 * |
核武器运输******事故气溶胶扩散数值模拟;唐秀欢 等;《安全与环境学报》;20140228;第14卷(第1期);第142-148页 * |
******条件下释放气溶胶扩散研究;谢朝阳 等;《中国安全科学学报》;20081031;第18卷(第10期);第87-91页 * |
脏弹***事件后果评估模型的构建与应用;苏自强 等;《四川环境》;20191231;第38卷(第6期);第65-71页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109580433A (en) | 2019-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109580433B (en) | Source term estimation method for diffusion of conventional explosive radioactive aerosol | |
Harper et al. | Realistic radiological dispersal device hazard boundaries and ramifications for early consequence management decisions | |
Papastefanou | Radioactive aerosols | |
Shultis et al. | Radiation shielding technology | |
Green et al. | Overview of the full-scale radiological dispersal device field trials | |
Krest’yanikova et al. | Stochastic models of hot planetary and satellite coronas: A photochemical source of hot oxygen in the upper atmosphere of Mars | |
Di Giovanni et al. | Two realistic scenarios of intentional release of radionuclides (Cs-137, Sr-90)-the use of the HotSpot code to forecast contamination extent | |
Dvornik et al. | Assessment of 137Cs contamination of combustion products and air pollution during the forest fires in zones of radioactive contamination | |
Stephens | Secondary antiproton spectrum in the atmosphere with the inclusion of Fermi motion in air nuclei | |
CN107145613B (en) | TNT equivalent based nuclear fission explosion source item estimation method | |
Kolbin et al. | The behaviour of Eu, Pu, Am radionuclide at burning radioactive graphite in an oxygen atmosphere. Computer experiments | |
Ouf et al. | Contribution to the study of particle resuspension kinetics during thermal degradation of polymers | |
Mishima | A review of research on plutonium releases during overheating and fires | |
Tian et al. | Experimental Study on the Natural Deposition Characteristics of Erosols in Containment | |
Vasconcelos et al. | Coulomb-nuclear interference in the scattering of 3He and α-particles by 122,124 Te, 124Sn and 114Cd | |
Dwivedi et al. | Optimization of controlling parameters of plasma torch aerosol generator and characteristics of synthesized metal oxide aerosols in context of NAF program | |
Elokhin | On the longitudinal stability of technogenic ionized formations | |
Hale et al. | DETERMINATION OF 239PU AIRBORNE CONCENTRATION ALPHA CORRECTION FACTOR FOR A ZINC SULFIDE DETECTOR VIA AMBIENT 222RN PROGENY AIR SAMPLING | |
Condit et al. | Pseudo-Evaporation of High Specific Activity Alpha-Emitting Materials | |
Skubacz | Transmission of ultrafine particles through separating systems of dust samplers | |
Ogorodnikov et al. | Radon and Its Effect on the Radiation Situation in the Chernobyl NPP Sarcophagus | |
Kogan et al. | Plutonium release fractions from accidental fires | |
Parkhurst et al. | Physicochemical characterization of Capstone depleted uranium aerosols I: uranium concentration in aerosols as a function of time and particle size | |
Suhana et al. | Evaluation of radiological hazards of particulates emissions from a coal fired power plant | |
Jallu et al. | The use of non-destructive passive neutron measurement methods in dismantling and radioactive waste characterization |
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 |