CN112614604B - Pulse extraction column organic phase ignition accident source item estimation method - Google Patents

Pulse extraction column organic phase ignition accident source item estimation method Download PDF

Info

Publication number
CN112614604B
CN112614604B CN202011346187.2A CN202011346187A CN112614604B CN 112614604 B CN112614604 B CN 112614604B CN 202011346187 A CN202011346187 A CN 202011346187A CN 112614604 B CN112614604 B CN 112614604B
Authority
CN
China
Prior art keywords
mass
smoke
time step
calculating
current time
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
Application number
CN202011346187.2A
Other languages
Chinese (zh)
Other versions
CN112614604A (en
Inventor
梁博宁
张建岗
杨亚鹏
王任泽
冯宗洋
贾林胜
王宁
刘一宁
闫峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Institute for Radiation Protection
Original Assignee
China Institute for Radiation Protection
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Institute for Radiation Protection filed Critical China Institute for Radiation Protection
Priority to CN202011346187.2A priority Critical patent/CN112614604B/en
Publication of CN112614604A publication Critical patent/CN112614604A/en
Application granted granted Critical
Publication of CN112614604B publication Critical patent/CN112614604B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G21NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21FPROTECTION AGAINST X-RADIATION, GAMMA RADIATION, CORPUSCULAR RADIATION OR PARTICLE BOMBARDMENT; TREATING RADIOACTIVELY CONTAMINATED MATERIAL; DECONTAMINATION ARRANGEMENTS THEREFOR
    • G21F9/00Treating radioactively contaminated material; Decontamination arrangements therefor
    • G21F9/007Recovery of isotopes from radioactive waste, e.g. fission products
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • Y02E30/30Nuclear fission reactors

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

The invention relates to a method for estimating an organic phase fire accident source item of a pulse extraction column, which comprises the following steps: (1) starting; (2) acquiring process parameters; (3) Estimating the mass of the solvent and the surface area of the combustion pool of the current leakage; (4) calculating the x type smoke mass generated by fuel combustion; (5) Calculating the residual mass of the xth flue gas after the flue gas is filtered by a filter and is settled by a pipeline; (6) Calculating the activity of the released xth smoke and the total release activity of all kinds of smoke after the smoke is filtered by a filter and is settled by a pipeline; (7) Estimating the residual oxygen content, and returning to the step (4) for re-estimation; (8) if the fuel is burnt out, the calculation is ended. According to the method, a reasonable estimation method is provided by analyzing the environmental characteristics and available state parameters of a pulse extraction column in a spent fuel post-treatment facility and based on the ignition rule of an organic phase (tributyl phosphate-kerosene mixed solution) of the pulse extraction column, and input conditions are provided for emergency response and emergency decision.

Description

Pulse extraction column organic phase ignition accident source item estimation method
Technical Field
The invention belongs to the technical field of nuclear and radiation emergency, and particularly relates to a pulse extraction column organic phase fire accident source item estimation method.
Background
When nuclear fuel is used in a reactor, the consumption of fissile nuclides and the generation of fission products and heavy nuclides cause the change of fuel reactivity, and finally the reactor can not maintain criticality any more, so that the nuclear fuel needs to be replaced to a certain extent. The fuel discharged after irradiation by the reactor is also called spent fuel or irradiated fuel. Spent nuclear fuel must be handled properly because it contains a large amount of radionuclides and is therefore highly radioactive. The spent fuel treatment mainly comprises the processes of storage, transportation, post-treatment, deep geological disposal and the like.
Spent fuel reprocessing is an important link for realizing the recycling and proper disposal of nuclear fuel. The Purex process (Purex) is a chemical process for recovering uranium and plutonium from irradiated nuclear fuel by tributyl phosphate extraction, and is the most effective and successful post-treatment process nowadays. A solvent extraction system, a solvent purification system and an equipment room in a spent fuel post-treatment facility contain a large amount of organic solvents, and potential fire accident potential of organic phase solvents exists. In the uranium and plutonium co-decontamination system, the position with the highest radioactive activity concentration of fission products is a pulse extraction column, the quantity of radioactive nuclides released from a burning solvent is the largest when solvent leakage catches fire in the device, and the main burning substance is a mixed solution of tributyl phosphate and kerosene. In the nuclear emergency condition, the source item estimation needs to be carried out on the accident so as to provide input conditions for emergency response and emergency decision.
Disclosure of Invention
Aiming at the problems in the prior art, the invention aims to provide a pulse extraction column organic phase ignition accident source item estimation method, which is a reasonable estimation method based on the ignition rule of a pulse extraction column organic phase (tributyl phosphate-kerosene mixed solution) by analyzing the environmental characteristics and available state parameters of the pulse extraction column in a spent fuel post-treatment facility and providing input conditions for emergency response and emergency decision.
In order to achieve the above purposes, the invention adopts a technical scheme that: a method for estimating the source term of an organic phase fire accident of a pulse extraction column comprises the following steps:
(1) Judging whether the pulse extraction column has a fire accident or not, and starting calculation if the fire accident occurs;
(2) Acquiring process parameters;
(3) According to pulseExtraction column quality on-line monitoring device and sump liquid level alarm device data estimation current leaked solvent quality M f,0 And the surface area S of the combustion chamber pool
(4) Calculating the x-th smoke mass m generated by fuel combustion at the current time step x,a,t
(5) Calculating the residual mass m of the xth smoke after the xth smoke is filtered by the filter of the exhaust pipeline and subsides in the pipeline and the xth smoke at the current time step x,a,t,J
(6) Calculating the activity A of the released x type of smoke at the current time step after the smoke is filtered by a filter of the exhaust pipeline and is settled by the pipeline x,t And the total release activity A of all kinds of smoke at the current time step t
(7) Judging whether the fuel is burnt out or not, if not, calculating the mass of the residual available oxygen in the equipment room at the current time step
Figure BDA0002799996210000021
And initial available oxygen mass
Figure BDA0002799996210000022
Returning to the step (4) for recalculation;
(8) And if the fuel is burnt out, calculating the total activity A of the organic phase ignition release at all time steps, and finishing the calculation.
Further, one embodiment of the step (1) is to manually start the calculation.
Further, another embodiment in step (1) is to automatically start the calculation by detecting the occurrence of combustion in the pulsed extraction column through the drift behavior of the process monitoring instrument and the increase of the radioactivity of the aerosol in the equipment room.
Further, the process parameters in the step (2) include fixed process parameters, real-time process parameters and real-time monitoring parameters.
Further, the fixed process parameters comprise the height of the equipment room and the volume of the equipment room, the real-time process parameters comprise the volume fraction of tributyl phosphate and the initial mass of the organic solvent of the pulse extraction column in the current process flow, and the real-time monitoring parameters comprise the inlet air flow and the outlet air flow of the equipment room.
Further, the mass M of the solvent currently leaked in the step (3) f,0 And the surface area S of the combustion chamber pool The calculating method comprises the following steps:
M f,0 =M o -M l
Figure BDA0002799996210000031
wherein M is 0 Is the initial mass of the organic solvent of the pulsed extraction column of the current process flow, M 1 Is the current residual organic solvent mass of the pulse extraction column, rho f Is the density of the organic solvent, H pool Is the height of the combustion chamber.
Further, the specific steps of the step (4) include:
(41) Calculating the fuel mass loss rate m' t
Figure BDA0002799996210000032
Where t is the current time step, h c For combustion convection heat flux density, h f For combustion of radiant heat flux, h r Is the surface radiant heat flux density, h l For latent heat of vaporization of fuel, T t For the solution temperature at the current time step, c f (T t ) Is heat capacity of fuel, T b Is the boiling point of the fuel, T In order to be at an infinite space temperature,
Figure BDA0002799996210000033
is the mass ratio of nitrate solution;
(42) Calculating the fuel combustion mass m at the current time step t With the mass of oxygen consumed by combustion of the fuel
Figure BDA0002799996210000034
m t =m" t ·Δt·S pool
Figure BDA0002799996210000035
Wherein, Δ t is the current time step duration,
Figure BDA0002799996210000036
mass of oxygen consumed for combustion of a unit mass of fuel;
(43) Calculating the total quantity of heat Q generated by combustion t And the solution temperature T at the next time step t+1
Q t =am t h t
Figure BDA0002799996210000037
Wherein a is the heat release efficiency, h t As heat of combustion of fuel, M room All the equipment quality (including iron on the ground, wall surface and the like) of the equipment room, c s (T t ) The heat capacity of steel is adopted;
(44) Calculating mass fraction omega of fuel combustion at current time step t
Figure BDA0002799996210000041
Wherein M is f,t The initial remaining fuel mass for the current time step,
mass of oxygen consumed by fuel combustion if current time step
Figure BDA0002799996210000042
Greater than the current remaining available oxygen mass
Figure BDA0002799996210000043
The mass of oxygen consumed by the fuel burning at the current time step
Figure BDA0002799996210000044
Is equal to the currentMass of remaining available oxygen
Figure BDA0002799996210000045
Namely, it is
Figure BDA0002799996210000046
(45) Calculating the x-th smoke mass m generated at the current time step x,a,t
m x,a,t =M x,t ·ω t ·γ x
Wherein M is x,t For the remaining mass of the x-th radionuclide, gamma, at the current time step x The smoke generation rate for the x-th radionuclide.
Further, the specific steps of the step (5) include:
(51) Calculating the residual smoke mass of the smoke before the smoke passes through the jth filter
Figure BDA0002799996210000047
Figure BDA0002799996210000048
Wherein m is x,a,t,j-1 The mass of the residual smoke, zeta, after the xth smoke passes through the jth-1 filter x,j The deposition rate per unit path of the x type flue gas in the j section of the pipeline, l j Is the length of the pipe before the jth filter;
(52) Calculating the x-th smoke mass Deltam filtered by the jth filter x,a,t,j
Figure BDA0002799996210000049
Wherein,
Figure BDA00027999962100000410
the filtering efficiency of the jth filter on the xth type of smoke;
(53) Calculating the residual smoke mass m after passing through the jth filter x,a,t,j
Figure BDA0002799996210000051
Wherein, ω is x Is the mass fraction of the x-th radionuclide in the oxide fume;
(54) Calculating the remaining effective filter mass m for the jth filter p,t,j
Figure BDA0002799996210000052
Wherein m is p,t-1,j For the effective filtering quality remaining one time step before the jth filter,
Figure BDA0002799996210000053
all the filtered smoke mass for the jth filter;
if the filter has a remaining effective filter mass m p,t,j The mass of the filtered smoke is less than or equal to the mass of the filtered smoke, and the mass of the filtered smoke is set as the residual effective filtering mass m of the filter p,t,j And the filter is failed;
(55) After total J filter filtration and pipeline settlement, the residual mass of the x-th smoke at the current time step is m x,a,t,J
Further, the activity A released by the xth smoke in the current time step in the step (6) x,t And the total release activity A of all kinds of smoke at the current time step t The calculating method comprises the following steps:
A x,t =a x m x,a,t,J
Figure BDA0002799996210000054
wherein, a x The specific activity of the X-th radionuclide is more than or equal to 1 and less than or equal to X.
Further, the residual available oxygen in the equipment room at the current time step in the step (7)Quality of
Figure BDA0002799996210000055
And initial available oxygen mass
Figure BDA0002799996210000056
The calculation method comprises the following steps:
for organic solvents, the organic solvents cannot be combusted when the oxygen content in the air reaches eta, so the oxygen content of the exhausted air is supposed to be eta,
Figure BDA0002799996210000057
Figure BDA0002799996210000058
wherein,
Figure BDA0002799996210000061
for the quality of the residual available oxygen in the equipment room at the previous time step, V in Is the inlet air flow, ρ in Is the charge density, V room Is the equipment room air volume.
Further, the method for calculating the total activity a of the organic phase ignition release in the step (8) comprises the following steps:
Figure BDA0002799996210000062
wherein t is more than 0 and less than or equal to t t ,t t The total time for the organic phase to ignite and burn is shown.
The invention has the beneficial effects that: the method provides a reasonable estimation method by analyzing the environmental characteristics and available state parameters of a pulse extraction column in the spent fuel post-treatment facility and based on the ignition rule of an organic phase (tributyl phosphate-kerosene mixed solution) of the pulse extraction column, and provides input conditions for emergency response and emergency decision.
Drawings
FIG. 1 is a flow chart of a method for estimating the source of an organic phase fire accident in a pulsed extraction column according to the present invention.
Detailed Description
In order to make the technical problems solved, the technical solutions adopted, and the technical effects achieved by the present invention clearer, the technical solutions of the embodiments of the present invention will be described in further detail below with reference to the accompanying drawings. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, are within the scope of the present invention.
FIG. 1 is a flow chart of a method for estimating the source of an organic phase fire accident in a pulsed extraction column according to the present invention, the method comprising the following steps:
(1) And judging whether the pulse extraction column has a fire accident or not, and starting calculation if the fire accident occurs.
In the starting module, no interface function is used, and the starting module is started by manual starting or by detecting that the pulse extraction column is burnt through the drifting behavior (pulse extraction column quality and sump liquid level alarm device) of a process monitoring instrument and the radioactive rising phenomenon of the aerosol in the equipment room.
(2) Obtaining process parameters
The acquired process parameters mainly comprise fixed process parameters, real-time process parameters and real-time monitoring parameters.
The fixed process parameters comprise the height of the equipment room, the volume of the equipment room and the like, the real-time process parameters comprise the volume fraction of tributyl phosphate, the initial mass of the organic solvent of the pulse extraction column in the current process flow and the like, and the real-time monitoring parameters comprise the inlet air flow, the outlet air flow and the like of the equipment room.
(3) Estimating the mass M of the solvent leaked at present according to the data of the pulse extraction column quality on-line monitoring device and the water sump liquid level alarm device f,0 (i.e., initial mass of fuel) and combustion bowl surface area S pool The calculation method comprises the following steps:
M f,0 =M o -M l
Figure BDA0002799996210000071
wherein,
M 0 the unit is Kg which is the initial mass of the organic solvent of the pulse extraction column in the current process flow;
M 1 the unit is Kg, which is the mass of the current residual organic solvent of the pulse extraction column;
ρ f is the density of the organic solvent and has the unit of Kg/m 3
H pool Is the height of the combustion chamber in m.
(4) Calculating the x-th smoke mass m generated by fuel combustion at the current time step x,a,t The method comprises the following specific steps:
(41) Firstly, a heat balance formula and a combustion rate empirical formula are combined, the influence of the nitric acid solution on the combustion rate is approximately considered, and the fuel mass loss rate m is calculated " t The calculation method comprises the following steps:
Figure BDA0002799996210000072
wherein,
t is the current time step and is dimensionless;
h c for combustion convective heat flow density, in Kw/m 2
h f The density of heat flow of combustion radiation is given in Kw/m 2
h r Is the surface radiant heat flux density in Kw/m 2
h l The unit is Kj/Kg for the latent heat of fuel evaporation;
T t the solution temperature at the current time step is shown in unit;
c f (T t ) The heat capacity of the fuel is Kj/(Kg. DEG C);
T b is the fuel boiling point, in units of ℃;
T is an infinite space temperature in units of;
Figure BDA0002799996210000081
is the mass ratio of the nitrate solution and has no dimension.
(42) The fuel combustion quality m at the current time step can then be obtained t With the mass of oxygen consumed by combustion of the fuel
Figure BDA0002799996210000082
The calculation method comprises the following steps:
m t =m t "·Δt·S pool
Figure BDA0002799996210000083
wherein,
delta t is the duration of the current time step, and the unit is s;
Figure BDA0002799996210000084
the mass of oxygen consumed by the combustion of fuel of unit mass is dimensionless;
(43) The total heat quantity Q generated by combustion t And the solution temperature T at the next time step t+1 Comprises the following steps:
Q t =am t h t
Figure BDA0002799996210000085
wherein,
a is the heat release efficiency, and is dimensionless;
h t the unit is Kj/Kg as the combustion heat of the fuel;
M room the unit is Kg, which is the mass (including iron on the ground, wall and the like) of all equipment in the equipment room;
c s (T t ) The heat capacity of the steel is expressed in Kj/(Kg. DEG C));
(44) The mass fraction omega of the fuel combustion at the current time step can be obtained t Comprises the following steps:
Figure BDA0002799996210000086
wherein M is f,t The initial residual fuel mass in Kg at the current time step;
mass of oxygen consumed by fuel combustion if current time step
Figure BDA0002799996210000087
Greater than the current remaining available oxygen mass
Figure BDA0002799996210000088
The mass of oxygen consumed by the fuel burning at the current time step
Figure BDA0002799996210000089
Equal to the current remaining available oxygen mass
Figure BDA00027999962100000810
Namely, it is
Figure BDA00027999962100000811
(45) Calculating the x-th smoke mass m generated at the current time step x,a,t
The smoke generated by X radioactive substances and the smoke generated by an organic solvent are shared, wherein the X +1 term is the smoke generated by solvent combustion, the smoke has no radioactivity, and the mass m of the X smoke generated by the current time step x,a,t Comprises the following steps:
m x,a,t =M x,t ·ω t ·γ x
wherein,
M x,t the residual mass of the x-th radionuclide in the current time step is Kg;
γ x the smoke generation rate of the x-th radionuclide is dimensionless.
(5) Calculating the residual mass m of the x-th smoke filtered by the filter of the exhaust pipeline and settled by the pipeline at the current time step x,a,t,J
Since a large amount of flue gas generated by combustion is discharged from the exhaust duct, the duct has a filter, and a part of the flue gas is deposited on the filter, it is necessary to estimate the deposition in the duct and the filter. In addition, the flue gas may cause clogging of the filter, and therefore an evaluation of the filter is also required. The method comprises the following specific steps:
(51) Firstly, the residual smoke mass of the smoke before the smoke passes through the jth filter is calculated
Figure BDA0002799996210000091
Figure BDA0002799996210000092
Wherein,
m x,a,t,j-1 the unit is Kg, which is the mass of the smoke left after the xth smoke passes through the jth-1 filter;
ζ x,j the deposition rate of the xth section of pipeline per unit path of the xth flue gas is m -1
l j Is the length of the pipeline before the jth filter, and is expressed in m;
(52) Calculating the mass delta m of the xth type smoke filtered by the jth filter x,a,t,j
Figure BDA0002799996210000093
Wherein,
Figure BDA0002799996210000094
the filter efficiency of the jth filter to the xth type of smoke is dimensionless;
(53) Calculating the residual smoke mass m after passing through the jth filter x,a,t,j
Figure BDA0002799996210000095
Wherein, ω is x Is the mass fraction of the x-th radionuclide in the oxide smoke, and is dimensionless;
(54) Calculating the effective filter mass m remaining for the jth filter p,t,j
Figure BDA0002799996210000101
Wherein m is p,t-1,j Is the remaining effective filtration mass in Kg for the previous time step of the jth filter;
Figure BDA0002799996210000102
the unit is Kg for all the smoke filtered by the jth filter;
if the remaining effective filter mass m of the filter p,t,j Setting the mass of the filtered smoke as the residual effective filtering mass m of the filter when the mass of the filtered smoke is less than or equal to the mass of the filtered smoke p,t,j And the filter fails.
After filtering by J filters and pipeline sedimentation, the residual mass of the x-th smoke at the current time step is m x,a,t,J
(6) Calculating the activity A of the released x type of smoke at the current time step after the smoke is filtered by a filter of the exhaust pipeline and is settled by the pipeline x,t And the total release activity A of all kinds of smoke at the current time step t The calculation method comprises the following steps:
A x,t =a x m x,a,t,J
Figure BDA0002799996210000103
wherein X is more than or equal to 1 and less than or equal to X, a x The radioactivity of the x type radionuclide is expressed in Bq/Kg;
(7) Judging whether the fuel is burnt out, if not, calculating the current time step equipment roomMass of available oxygen remaining in the reactor
Figure BDA0002799996210000104
And initial available oxygen mass
Figure BDA0002799996210000105
And (5) returning to the step (4) for recalculation.
For organic solvents, the organic solvent cannot be combusted when the oxygen content in the air reaches eta, so that the mass of the residual available oxygen in the equipment room at the current time step is the mass of the residual available oxygen in the equipment room at the current time step on the assumption that the oxygen content of the exhausted air is eta
Figure BDA0002799996210000106
And initial available oxygen mass
Figure BDA0002799996210000107
Comprises the following steps:
Figure BDA0002799996210000108
Figure BDA0002799996210000109
wherein,
Figure BDA0002799996210000111
the mass of the residual available oxygen in the equipment room at the previous time step is Kg;
V in is the intake air flow rate, and has the unit of m 3 /s;
ρ in The unit is Kg/m for the intake air density 3
V room Is the air volume of the equipment room, and has a unit of m 3
(8) And if the fuel is burnt out, calculating the total activity A released by the ignition of the organic phase, and outputting a calculation result.
Total activity A of each radionuclide released by organic phase ignition x Comprises the following steps:
Figure BDA0002799996210000112
the total activity a of the organic phase on fire release was:
Figure BDA0002799996210000113
wherein X is more than or equal to 1 and less than or equal to X, t is more than 0 and less than or equal to t t ,t t The total time of the organic phase ignition combustion is expressed in s.
It will be appreciated by persons skilled in the art that the apparatus and method of the present invention are not limited to the embodiments described in the detailed description, and the detailed description is for the purpose of explanation and not limitation of the invention. Other embodiments will be apparent to those skilled in the art from the following detailed description, which are also included in the scope of the invention as defined in the appended claims.

Claims (7)

1. A method for estimating the source term of an organic phase fire accident of a pulse extraction column is characterized by comprising the following steps:
(1) Judging whether the pulse extraction column has a fire accident or not, and starting calculation if the fire accident occurs;
(2) Acquiring process parameters;
(3) Estimating the mass M of the solvent leaked at present according to the data of the pulse extraction column quality on-line monitoring device and the water sump liquid level alarm device f,0 And the surface area S of the combustion chamber pool
(4) Calculating the mass m of the x-th smoke generated by fuel combustion at the current time step x,a,t The method comprises the following specific steps:
(41) Calculating the fuel mass loss rate m' t
Figure FDA0003727544810000011
Wherein t is the current time step, h c To be burnedConvective heat flux density, h f For combustion of radiant heat flux density, h r Is the surface radiant heat flux density, h l For latent heat of vaporization of fuel, T t For the solution temperature at the current time step, c f (T t ) Is the heat capacity of the fuel, T b Is the boiling point of the fuel, T In order to achieve an infinite space temperature,
Figure FDA0003727544810000012
is the mass ratio of nitrate solution;
(42) Calculating the fuel combustion mass m at the current time step t With the mass of oxygen consumed by combustion of the fuel
Figure FDA0003727544810000013
m t =m″ t ·Δt·S pool
Figure FDA0003727544810000014
Wherein, Δ t is the current time step duration,
Figure FDA0003727544810000015
mass of oxygen consumed for combustion of a unit mass of fuel;
(43) Calculating the total quantity of heat Q generated by combustion t And the solution temperature T of the next time step t+1
Q t =am t h t
Figure FDA0003727544810000016
Wherein a is the heat release efficiency, h t As heat of combustion of fuel, M room For all equipment masses of the equipment room, c s (T t ) Is the heat capacity of steel, M f,t An initial remaining fuel mass for the current time step;
(44) Calculating the current time step fuelMass fraction of combustion omega t
Figure FDA0003727544810000021
Mass of oxygen consumed by fuel combustion if current time step
Figure FDA0003727544810000022
Greater than the current remaining available oxygen mass
Figure FDA0003727544810000023
The mass of oxygen consumed by the fuel burning at the current time step
Figure FDA0003727544810000024
Equal to the current remaining available oxygen mass
Figure FDA0003727544810000025
Namely, it is
Figure FDA0003727544810000026
(45) Calculating the x-th smoke mass m generated at the current time step x,a,t
m x,a,t =M x,t ·ω t ·γ x
Wherein M is x,t For the remaining mass of the x-th radionuclide, gamma, at the current time step x Smoke generation rate for the xth radionuclide;
(5) Calculating the residual mass m of the x-th smoke at the current time step after the x-th smoke is filtered by a filter of the exhaust pipeline and is settled by the pipeline x,a,t,J The method comprises the following specific steps:
(51) Calculating the residual smoke mass of the smoke before the smoke passes through the jth filter
Figure FDA0003727544810000027
Figure FDA0003727544810000028
Wherein m is x,a,t,j-1 The mass of the residual smoke, zeta, after the xth smoke passes through the jth-1 filter x,j The deposition rate per unit path of the x type flue gas in the j section of the pipeline, l j Is the length of the pipe before the jth filter;
(52) Calculating the x-th smoke mass Deltam filtered by the jth filter x,a,t,j
Figure FDA0003727544810000029
Wherein,
Figure FDA00037275448100000210
the filtering efficiency of the jth filter on the xth type of smoke;
(53) Calculating the residual smoke mass m after passing through the jth filter x,a,t,j
Figure FDA00037275448100000211
Wherein, ω is x Is the mass fraction of the x-th radionuclide in the oxide fume;
(54) Calculating the effective filter mass m remaining for the jth filter p,t,j
Figure FDA00037275448100000212
Wherein m is p,t-1,j To provide the effective filtering quality remaining for the time step preceding the jth filter,
Figure FDA00037275448100000213
for the jth filterThe filtered smoke quality is generated by X radionuclides;
if the filter has a remaining effective filter mass m p,t,j Setting the mass of the filtered smoke as the residual effective filtering mass m of the filter when the mass of the filtered smoke is less than or equal to the mass of the filtered smoke p,t,j And the filter is failed;
(55) After filtering by J filters and pipeline sedimentation, the residual mass of the x-th smoke at the current time step is m x,a,t,J
(6) According to the residual mass m of the x-th smoke at the current time step x,a,t,J Calculating the activity A of the x-th smoke released at the current time step x,t And the total release activity A of all kinds of smoke at the current time step t
(7) Judging whether the fuel is burnt out or not, if not, calculating the mass of the residual available oxygen in the equipment room at the current time step
Figure FDA0003727544810000031
And initial available oxygen mass
Figure FDA0003727544810000032
Returning to the step (4) to recalculate,
Figure FDA0003727544810000033
and
Figure FDA0003727544810000034
the calculating method comprises the following steps:
for organic solvents, when the oxygen content in the air reaches eta, the organic solvents cannot be combusted, and the oxygen content of the discharged air is eta,
Figure FDA0003727544810000035
Figure FDA0003727544810000036
wherein,
Figure FDA0003727544810000037
for the remaining available oxygen mass, V, in the equipment room at the previous time step in Is the inlet air flow, ρ in Is the intake density, V room Is the equipment room air volume;
(8) And if the fuel is burnt out, calculating the total activity A of the organic phase ignition release at all time steps, and finishing the calculation.
2. The method as claimed in claim 1, wherein the start-up modes in step (1) include a manual start-up mode and an automatic start-up mode in which the start-up mode is automatically started up when the combustion of the pulsed extraction column is detected by drift behavior of a process monitoring instrument and a rise in activity of an aerosol in an equipment room.
3. The method as claimed in claim 1, wherein the process parameters in step (2) include fixed process parameters, real-time process parameters and real-time monitoring parameters.
4. The method as claimed in claim 3, wherein the fixed process parameters include height of the equipment room and volume of the equipment room, the real-time process parameters include volume fraction of tributyl phosphate in the pulsed extraction column and initial mass of organic solvent in the pulsed extraction column in the current process flow, and the real-time monitoring parameters include inlet air flow and outlet air flow of the equipment room.
5. The method for estimating the source term of organic phase fire accident of pulse extraction column as claimed in claim 1, wherein the mass M of solvent currently leaked in step (3) f,0 And the surface area S of the combustion chamber pool The calculating method comprises the following steps:
M f,0 =M o -M l
Figure FDA0003727544810000041
wherein M is 0 Is the initial mass of the organic solvent of the pulsed extraction column of the current process flow, M 1 Is the current residual organic solvent mass of the pulse extraction column, rho f Density of organic solvent, H pool Is the height of the combustion pool.
6. The method for estimating the source of the organic phase fire accident of the pulse extraction column according to claim 1, wherein the activity A of the released activity of the xth smoke in the current time step in the step (6) x,t And the total release activity A of all kinds of smoke at the current time step t The calculation method comprises the following steps:
A x,t =a x m x,a,t,J
Figure FDA0003727544810000042
wherein, a x The specific activity of the X-th radionuclide is more than or equal to 1 and less than or equal to X.
7. The method for estimating the source term of an organic phase fire accident of a pulse extraction column according to claim 1, wherein the total activity A of the organic phase fire release in all time steps in the step (8) is calculated by:
Figure FDA0003727544810000043
wherein t is more than 0 and less than or equal to t t ,t t The total time for the organic phase to ignite and burn is shown.
CN202011346187.2A 2020-11-26 2020-11-26 Pulse extraction column organic phase ignition accident source item estimation method Active CN112614604B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011346187.2A CN112614604B (en) 2020-11-26 2020-11-26 Pulse extraction column organic phase ignition accident source item estimation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011346187.2A CN112614604B (en) 2020-11-26 2020-11-26 Pulse extraction column organic phase ignition accident source item estimation method

Publications (2)

Publication Number Publication Date
CN112614604A CN112614604A (en) 2021-04-06
CN112614604B true CN112614604B (en) 2022-10-21

Family

ID=75225306

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011346187.2A Active CN112614604B (en) 2020-11-26 2020-11-26 Pulse extraction column organic phase ignition accident source item estimation method

Country Status (1)

Country Link
CN (1) CN112614604B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104061588A (en) * 2014-07-17 2014-09-24 烟台龙源电力技术股份有限公司 Low-nitrogen combustion control method and system based on secondary air door air regulation control
CN110059347A (en) * 2019-03-11 2019-07-26 中国辐射防护研究院 The experiment research platform and method of fire incident in a kind of simulation reprocessing facility
CN110953041A (en) * 2019-12-17 2020-04-03 凯龙高科技股份有限公司 Non-road diesel engine combustor DPF system control method

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NL295851A (en) * 1962-07-26 1900-01-01
NL301862A (en) * 1962-12-26 1900-01-01
US3359078A (en) * 1965-02-10 1967-12-19 Gen Electric Irradaiated nuclear fuel recovery
FR2913971B1 (en) * 2007-03-20 2009-04-24 Saint Gobain GLASS FUSION DEVICE COMPRISING TWO OVENS
US8999168B2 (en) * 2010-10-29 2015-04-07 Ut-Battelle, Llc Supported liquid inorganic membranes for nuclear waste separation
CN102841600B (en) * 2012-08-24 2014-12-17 中国核电工程有限公司 Accidence safety analysis method for nuclear fuel reprocessing plant
KR101736168B1 (en) * 2016-07-28 2017-05-17 한전원자력연료 주식회사 Pulsed column having apparatus for supplying pulse
CN107290990B (en) * 2017-05-03 2021-03-19 中国核电工程有限公司 Measurement and control method of pulse extraction column for nuclear chemical industry
CN108766612A (en) * 2018-04-17 2018-11-06 中国核电工程有限公司 A kind of novel nuclear fuel reprocessing plant's evaporator process detection and debugging system
CN109521053B (en) * 2018-09-21 2021-11-16 中国辐射防护研究院 Experimental device and method for researching post-treatment organic phase fire accident secondary deflagration
CN109509567A (en) * 2018-09-21 2019-03-22 中国辐射防护研究院 The device and method of simulation post-processing organic phase combustion emissions aerosol release
CN109243655A (en) * 2018-09-21 2019-01-18 中国辐射防护研究院 A kind of device and method of research post-processing organic phase ignition behavior
CN110070269B (en) * 2019-03-28 2022-12-09 中国辐射防护研究院 Management optimization method for equipment in post-processing facility based on risk guidance
CN111666531B (en) * 2020-04-23 2023-08-25 中冶华天工程技术有限公司 Fuel combustion calculation method under condition of high content of residual CO in flue gas
CN111423338A (en) * 2020-04-30 2020-07-17 山东微井化工科技股份有限公司 Industrial production system of diazomethane

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104061588A (en) * 2014-07-17 2014-09-24 烟台龙源电力技术股份有限公司 Low-nitrogen combustion control method and system based on secondary air door air regulation control
CN110059347A (en) * 2019-03-11 2019-07-26 中国辐射防护研究院 The experiment research platform and method of fire incident in a kind of simulation reprocessing facility
CN110953041A (en) * 2019-12-17 2020-04-03 凯龙高科技股份有限公司 Non-road diesel engine combustor DPF system control method

Also Published As

Publication number Publication date
CN112614604A (en) 2021-04-06

Similar Documents

Publication Publication Date Title
CN112614604B (en) Pulse extraction column organic phase ignition accident source item estimation method
CN109313945A (en) Nuclear fuel salt
Cipiti et al. Fusion transmutation of waste: design and analysis of the in-zinerator concept
Ayer et al. Nuclear fuel cycle facility accident analysis handbook
Kim et al. Development of a Korean roadmap for technical issue resolution for fission product behavior during severe accidents
Finney et al. Correlation of radioactive waste treatment costs and the environmental impact of waste effluents in the nuclear fuel cycle for use in establishing''as low as practicable''guides: nuclear fuel reprocessing
CN109300564B (en) Device and method for simulating steam blocking and corrosion of filter
JP2021032750A (en) Method and system for processing spent fuel
Davis Jr et al. Correlation of radioactive waste treatment costs and the environmental impact of waste effluents in the nuclear fuel cycle: reprocessing of high-temperature gas-cooled reactor fuel containing U-233 and thorium
Vargo A Brief History of Nuclear Criticality Accidents in Russia—1953–1997
Lian et al. Research and Development of Solvent Fire in the Reprocessing Plant
Arm Flowsheet Evaluation of Dissolving Used Nuclear Fuel in PUREX Solvent
Raab et al. OPERATING EXPERIENCE USING SILVER REACTORS FOR RADIOIODINE REMOVAL IN THE HANFORD PUREX PLANT.
Bower Chemical Technology Branch, Annual Report, Fiscal Year 1968
Cowley Organic solvent topical report
Hallman HB-Line Dissolver Dilution Flows and Dissolution Capability with Dissolver Charge Chute Cover Off
Pechin et al. Correlation of radioactive waste treatment costs and the environmental impact of waste effluents in the nuclear fuel cycle for use in establishing''as low as practicable''guides: fabrication of light-water reactor fuel from enriched uranium dioxide
Shida et al. A development and an application of MIXSET-X computer code for simulating the PUREX solvent extraction system
GB2257561A (en) Process and apparatus for the continuous elimination of the radioactive iodine contained in irradiation nuclear fuel elements
Bancroft et al. The Volatilization and Collection of Ruthenium and Cesium in a System for Incorporating Fission Products into Glass
Rosenberg et al. Air pollution control process and apparatus
Wavada Application of the extended Pitzer equation to nuclear fuel reprocessing
Sone et al. Reduction of radioactive secondary waste with steam reforming in treatment of waste TBP/dodecane
Zhang et al. Analysis of airborne radioactivity in containment and adjacent buildings under PWR accident conditions
Lackey et al. Assessment of coater size for the fuel refabrication prototype plant

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