CN111695762B - Correction method and device for nuclear accident diffusion result and result evaluation method and system - Google Patents
Correction method and device for nuclear accident diffusion result and result evaluation method and system Download PDFInfo
- Publication number
- CN111695762B CN111695762B CN202010355824.6A CN202010355824A CN111695762B CN 111695762 B CN111695762 B CN 111695762B CN 202010355824 A CN202010355824 A CN 202010355824A CN 111695762 B CN111695762 B CN 111695762B
- Authority
- CN
- China
- Prior art keywords
- decay
- sub
- release
- nth
- module
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 217
- 238000009792 diffusion process Methods 0.000 title claims abstract description 173
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000011156 evaluation Methods 0.000 title claims abstract description 27
- 230000002285 radioactive effect Effects 0.000 claims abstract description 71
- 230000008021 deposition Effects 0.000 claims description 155
- 230000004907 flux Effects 0.000 claims description 89
- 238000004364 calculation method Methods 0.000 claims description 50
- 239000000941 radioactive substance Substances 0.000 claims description 23
- 230000009471 action Effects 0.000 claims description 16
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000013507 mapping Methods 0.000 claims description 7
- 238000007654 immersion Methods 0.000 claims description 5
- 230000001681 protective effect Effects 0.000 claims description 5
- 239000012857 radioactive material Substances 0.000 claims description 5
- 101100134058 Caenorhabditis elegans nth-1 gene Proteins 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 13
- 238000004458 analytical method Methods 0.000 description 4
- 239000003344 environmental pollutant Substances 0.000 description 4
- 231100000719 pollutant Toxicity 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000005658 nuclear physics Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 239000003102 growth factor Substances 0.000 description 1
- 239000011261 inert gas Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Measurement Of Radiation (AREA)
Abstract
The invention provides a method for correcting a nuclear accident diffusion result, which comprises the following steps: acquiring radioactive release source item data; dividing the source item into M sub-release sections; obtaining diffusion results of M sub-release sections calculated by an atmospheric diffusion model, and enabling each sub-release section to have N groups of diffusion results; calculating decay multiples of the nth moment corresponding to N groups of diffusion results of each sub-release section respectively; multiplying the air concentration field at the nth moment of each sub-release section by the decay multiple of the moment to obtain an air concentration field correction value at the nth moment of the sub-release section; summing the air concentration field correction values at the nth moment of the M sub-release sections to obtain an air concentration field correction value at the nth moment; and traversing N times to obtain the air concentration field correction value of the radioactive release source item. The invention also provides a correction device, a nuclear accident result evaluation method and a system, which are used for solving the problem that the result calculated by the current nuclear accident result evaluation method is deviated from the actual physical process.
Description
Technical Field
The invention belongs to the field of analysis of nuclear facility accident radioactivity results, and particularly relates to a method and a device for correcting a nuclear accident radioactive substance atmospheric diffusion result and a method and a system for evaluating the nuclear accident result.
Background
Under the condition that the nuclear facility has accidents, the space-time distribution of the concentration of the radioactivity and the radiation dose around the nuclear facility is analyzed, predicted and evaluated according to the radioactive release source item and the meteorological condition, and the method is an important component in the processes of nuclear accident emergency preparation and response, post-accident evaluation, environmental impact evaluation and the like. In the nuclear engineering design work, the nuclear facility safety analysis, the environmental impact evaluation, the emergency plan division, the regional nuclear emergency scheme establishment, the three-level PSA analysis and the like all need to take reliable accident result evaluation results as data bases.
The existing result analysis and evaluation models at home and abroad simplify the processing of the decay process due to the complexity and complexity of the radionuclide decay process, so that the calculation result has deviation. Due to the difference of disciplines, the atmospheric diffusion model adopted by the current international mainstream mostly does not consider or only considers the radionuclide decay process in a limited way, meteorological data dynamically changes along with time in actual conditions, the diffusion tracks of radioactive substances released at different moments are different, the concentration at a certain moment is generated by superposition of a plurality of different release sections, and the decay time of the diffusion process is difficult to determine. The main atmospheric diffusion model is mostly built completely, cannot give intermediate results and is difficult to develop secondarily, and the method is introduced and applied to nuclear accident result evaluation and has limitation and difficulty.
Disclosure of Invention
The invention provides a correction method and device of a nuclear accident radioactive substance atmospheric diffusion result and a result evaluation method and system, which are used for solving the problem that the result calculated by the current nuclear accident result evaluation method is deviated from an actual physical process.
In order to achieve the above object, in a first aspect, an embodiment of the present invention provides a method for correcting an atmospheric diffusion result of a nuclear accident radioactive substance, including the steps of:
acquiring radioactive release source item data;
dividing the radioactive release source item data into M sub-release sections according to a first preset duration, wherein the first preset duration is positive integer multiple of delta t, and delta t is the minimum time interval of the calculation result of the atmospheric diffusion model;
obtaining diffusion results of M sub-release sections obtained through calculation of meteorological data and an atmospheric diffusion model, and enabling each sub-release section to have N groups of diffusion results with time intervals of a second preset duration, wherein the second preset duration is positive integer multiple of deltat, N is integer multiple of M, and the diffusion results comprise an air concentration field;
calculating decay multiples of the nth moment corresponding to N groups of diffusion results of each sub-release section respectively, wherein N is a positive integer, and N is more than or equal to 1 and less than or equal to N;
Performing a correction operation on the air concentration field at the nth time of each sub-release section: multiplying the obtained air concentration field in the diffusion result at the nth moment of each sub-release section by the decay multiple corresponding to the moment to obtain an air concentration field correction value at the nth moment of each sub-release section;
performing a summation operation on the air concentration field at the nth time of the M sub-release sections: summing the air concentration field correction values at the nth moment of the M sub-release sections to obtain the air concentration field correction value at the nth moment of the radioactive release source item;
and respectively executing correction operation and summation operation on the air concentration fields at N moments to obtain the air concentration field correction value of the radioactive release source item.
Preferably, the diffusion result further comprises a deposition flux field, and after obtaining the air concentration field correction value of the radioactive emission source term, the method further comprises:
performing a correction operation on the deposition flux field at the nth time of each sub-release segment: dividing the obtained deposition flux field in the diffusion result at the nth time of each sub-release section by the air concentration field at the time, and multiplying the obtained deposition flux field by the air concentration field correction value at the time to obtain the deposition flux field correction value at the nth time of each sub-release section;
Performing a summation operation on the deposited flux fields at the nth time of the M sub-release segments: summing the deposition flux field correction values at the nth time of the M sub-release sections to obtain the deposition flux field correction value at the nth time of the radioactive release source item;
and respectively performing correction operation and summation operation on the deposition flux fields at N moments to obtain the deposition flux field correction value of the radioactive release source item.
Preferably, after obtaining the deposition flux field correction value for the radioactive emission source term, the method further comprises:
performing decay correction with decay time of the ground deposition concentration field at the N-1 time as a second preset duration to obtain contribution from the ground deposition concentration field at the N-1 time in the ground deposition concentration field at the N-1 time, wherein N is a positive integer, and N is more than or equal to 1 and less than or equal to N;
multiplying the obtained deposition flux field correction value at the nth time of the radioactive release source item by a second preset time length to obtain the contribution of air deposition to the ground in the [ n-1, n ] period, and summing the contribution of the air deposition to the ground from the nth time to obtain the ground deposition concentration field at the nth time.
Preferably, the calculating the decay multiple of the nth time corresponding to the N groups of diffusion results of each sub-release section specifically includes:
According to decay chain information of all nuclides in the radioactive release source item data, each nuclide is regarded as a parent nucleus and decay trees are respectively constructed;
according to the decay tree of each parent nucleus, depth-first search is adopted, the decay times of the parent nucleus and the child nucleus at the nth moment corresponding to N groups of diffusion results of each child release section are calculated in parallel, and the decay times are stored in the decay tree.
Preferably, each nuclide is regarded as a parent nucleus and decay trees are respectively constructed, and specifically comprises:
obtaining multi-generation decay chain information of a parent nucleus, and storing each nuclide in the decay chain information as a node of a decay tree, wherein the parent nucleus is a root node of the decay tree, and the child nucleus is a child node of the decay tree;
storing all nodes in the decay tree as nodes in a unidirectional linked list, wherein each node stores node information, and the node information comprises physical properties such as the name of the nuclide, decay constant, decay branch ratio and the like, the position of the nuclide in the unidirectional linked list and the position of each sub-core of the nuclide in the unidirectional linked list;
and constructing the directed relation of all nodes according to the positions of the nuclides in the unidirectionally linked list and the positions of the child cores of the nuclides in the unidirectionally linked list so as to form the decay tree.
In a second aspect, an embodiment of the present invention provides a method for evaluating a nuclear accident outcome, including the steps of:
calculating an effective dose to be received by the public according to the air concentration field correction value, the ground deposition concentration field and the dose conversion factors corresponding to the air concentration field, the ground deposition concentration field and the inhalation amount of human body from air of the radioactive release source item obtained in the first aspect, wherein the effective dose comprises an air immersed external irradiation dose, a ground deposition external irradiation dose and an inhalation internal irradiation dose in an irradiated time period;
comparing the effective dose with the general optimized intervention level, and determining the protective action according to the comparison result.
Preferably, the comparing the effective dose to the universal optimized intervention level, and determining the protective action based on the comparison result specifically comprises:
if T 2d >10mSv, adopting concealing action; if T 7d In (2) with an evacuation action, wherein T2d represents the total effective dose for 2 days and T7d represents the total effective dose for 7 days.
In a third aspect, an embodiment of the present invention provides a device for correcting an atmospheric diffusion result of a nuclear accident radioactive substance, including a splitting module, an obtaining module, a calculating module, and a first correcting module.
The splitting module is used for acquiring the radioactive release source item data and splitting the radioactive release source item data into M sub-release sections according to a first preset duration, wherein the first preset duration is positive integer multiple of Deltat, and Deltat is the minimum time interval of the calculation result of the atmospheric diffusion model;
The acquisition module is connected with the splitting module and is used for acquiring diffusion results of M sub-release sections obtained through calculation of meteorological data and an atmospheric diffusion model, and enabling each sub-release section to have N groups of diffusion results with time intervals of a second preset duration, wherein the second preset duration is positive integer multiple of deltat, N is integer multiple of M, and the diffusion results comprise an air concentration field;
the calculation module is connected with the acquisition module and is used for calculating decay multiples of the nth moment corresponding to N groups of diffusion results of each sub-release section respectively, wherein N is a positive integer, and N is more than or equal to 1 and less than or equal to N;
the first correction module is connected with the acquisition module and the calculation module and is used for executing correction operation on the air concentration field at the nth moment of each sub-release section: the air concentration field in the diffusion result of the nth moment of each sub-release section obtained by the obtaining module is multiplied by the decay multiple corresponding to the moment calculated by the calculating module, so that an air concentration field correction value of the nth moment of each sub-release section is obtained; and performing a summation operation on the air concentration field at the nth time of the M sub-release sections: the air concentration field correction value at the nth moment of the M sub-release sections is summed to obtain the air concentration field correction value at the nth moment of the radioactive release source item; and the air concentration field correction value of the radioactive release source item is obtained by respectively performing correction operation and summation operation on the air concentration fields at N moments.
Preferably, the correction device further comprises a second correction module, the diffusion result further comprises a deposition flux field, and the second correction module is connected with the acquisition module, the calculation module and the first correction module and is used for performing correction operation on the deposition flux field at the nth time of each sub-release section: the method comprises the steps of dividing a deposition flux field in a diffusion result at the nth moment of each sub-release section obtained by an obtaining module by an air concentration field at the moment, and multiplying the deposition flux field by an air concentration field correction value at the nth moment obtained by a first correction module to obtain a deposition flux field correction value at the nth moment of each sub-release section; and performing a summation operation on the deposited flux fields at the nth time of the M sub-release segments: the method comprises the steps of summing the deposition flux field correction values at the nth moment of the M sub-release sections to obtain the deposition flux field correction value at the nth moment of the radioactive release source item; and the correction operation and the summation operation are respectively carried out on the deposition flux fields at N moments, so as to obtain the deposition flux field correction value of the radioactive release source item.
Preferably, the correction device further comprises a third correction module, wherein the third correction module is connected with the acquisition module and the second correction module and is used for carrying out decay correction on the ground deposition concentration field at the N-1 time as decay time of a second preset duration to obtain the contribution quantity from the ground deposition concentration field at the N-1 time in the ground deposition concentration field at the N-1 time, wherein N is a positive integer, and N is more than or equal to 1 and less than or equal to N; and the deposition flux field correction value at the nth moment of the radioactive release source item obtained by the second correction module is multiplied by a second preset duration and then summed with the contribution from the ground deposition concentration field at the nth-1 moment to obtain the ground deposition concentration field at the nth moment.
Preferably, the calculation module comprises a decay multiple unit, the decay multiple unit is connected with the splitting module and the acquisition module, and is used for acquiring decay chain information according to all nuclides in the release source item acquired by the splitting module, and regarding each nuclide as a parent nucleus according to the decay chain information and respectively constructing a decay tree; and the decay multiples of the parent nucleus and the child nucleus at the nth moment, which are respectively corresponding to N groups of diffusion results of each child release section, are calculated in parallel by adopting depth-first search according to the decay tree of each parent nucleus, and are stored in the decay tree.
In a fourth aspect, an embodiment of the present invention provides a nuclear accident outcome evaluation system, including: the correction device for the atmospheric diffusion result of the nuclear accident radioactive substance in the third aspect further comprises a dose module and an evaluation module, wherein the dose module is connected with the correction device for the atmospheric diffusion result of the nuclear accident radioactive substance and is used for calculating an effective dose received by the public according to dose conversion factors corresponding to the air concentration of each nuclide, the ground deposition concentration and the inhalation quantity of human body from air, and the air concentration field correction value and the ground deposition concentration field of the radioactive release source item obtained by the correction device, wherein the effective dose comprises an air immersed external irradiation dose, a ground deposition external irradiation dose and an inhalation internal irradiation dose in an illuminated time period; the evaluation module stores the general optimized intervention level, and pre-stores a mapping table of comparison results and protection actions, and the mapping table is connected with the dosage module and is used for comparing the general optimized intervention level with the effective dosage calculated by the dosage module and determining the corresponding protection actions according to the mapping table according to the comparison results.
According to the correction method for the atmospheric diffusion result of the nuclear accident radioactive substance, provided by the embodiment of the invention, the radioactive release source item is split into M sub-release sections according to the first preset time length, wherein the first preset time length is positive integer multiple of Deltat, deltat is the minimum time interval of the calculation result of the atmospheric diffusion model, so that the diffusion process is integrated into zero, and the simulation of the real diffusion process that the diffusion concentration at a certain place is generated by superposition of a plurality of different release sections is facilitated. And then obtaining diffusion results (the diffusion results comprise an air concentration field and a deposition flux field) of M sub-release sections obtained through calculation of meteorological data and an atmospheric diffusion model, and enabling each sub-release section to have N groups of diffusion results with time intervals of a second preset duration, wherein the second preset duration is a positive integer multiple of deltat, and N is an integer multiple of M. And calculating decay multiples of the nth moment corresponding to N groups of diffusion results of each sub-release section, wherein N is a positive integer, and N is more than or equal to 1 and less than or equal to N. Performing a correction operation on the air concentration field at the nth time of each sub-release section: multiplying the air concentration field in the obtained diffusion result at the nth moment of each sub-release section by the decay multiple corresponding to the moment to obtain an air concentration field correction value at the nth moment of each sub-release section, and then executing summation operation on the air concentration fields at the nth moment of the M sub-release sections: and summing the air concentration field correction values at the nth moment of the M sub-release sections to obtain the air concentration field correction value at the nth moment of the radioactive release source item, and respectively executing correction operation and summation operation on the air concentration fields at the nth moment to obtain the air concentration field correction value of the radioactive release source item. The deposition flux field correction value is then calculated from the air concentration field correction value. Compared with the diffusion result of the currently adopted atmospheric diffusion model, the diffusion result correction value calculated by the correction method is closer to the actual diffusion result, the calculation result is more accurate, and the calculation process is convenient and rapid only by carrying out corresponding correction on the existing calculation method. In addition, the method can directly carry out complete decay correction on the air concentration field and the deposition flux field output by the atmospheric diffusion model, does not need to change the atmospheric diffusion model, and has strong applicability.
Drawings
Fig. 1: the invention relates to a method for correcting the atmospheric diffusion result of nuclear accident radioactive substances, which is a flow chart of the method in the embodiment 1;
fig. 2: the embodiment 2 of the invention discloses a single linked list structure schematic diagram of a nuclear accident result evaluation method;
fig. 3: the decay tree structure schematic diagram of the nuclear accident consequence evaluation method of the embodiment 2 of the invention;
fig. 4: the decay time of the diffusion result of the sub-release section of the nuclear accident consequence evaluation method of the embodiment 2 of the invention is shown in the schematic diagram;
fig. 5: the invention relates to a flow chart of a nuclear accident consequence evaluation method in an embodiment 2;
fig. 6: the device for correcting the atmospheric diffusion result of a nuclear accident radioactive substance in embodiment 3 of the present invention.
Detailed Description
The present invention will be described in further detail below with reference to the drawings and examples for better understanding of the technical scheme of the present invention to those skilled in the art.
The decay process in the atmosphere diffusion of the radionuclide in reality is more complex, and the specific expression is as follows: the meteorological data dynamically changes along with time, the diffusion track of the radioactive substance released at different moments is different, and the concentration at a certain moment is generated by superposition of a plurality of different release sections, so that the decay time of the diffusion process and the like are difficult to determine. According to the method for correcting the atmospheric diffusion result of the nuclear accident radioactive substance, the diffusion process of the radioactive release source item is integrated into zero, the radioactive release source item is split into a plurality of sub-release sections, the time interval between the release end time and the diffusion result of each time obtained by independent diffusion is taken as the decay time, and then the diffusion result is corrected, wherein the diffusion result comprises an air concentration field and a deposition flux field.
Furthermore, it is known from the laws of nuclear physics that the decay times of a parent nucleus and its child nuclei (for child nuclei that do not exist before decay, which may also be referred to as growth factors, hereinafter collectively referred to as decay times) are only related to decay time. The inventor of the present application further obtains that the corrected concentration of a certain nuclide is only related to the concentration before correction and the decay multiple, and the sum of the corrected concentrations of the nuclide is generated for all parent nuclides, so that the correction algorithm provided by the embodiment of the present invention performs correction calculation of the air concentration field by multiplying the concentration before correction by the decay multiple at the corresponding moment.
Example 1:
as shown in fig. 1, the present embodiment provides a method for correcting the atmospheric diffusion result of nuclear accident radioactive substances. The correction method can be applied to a nuclear accident consequence evaluation method, and can directly carry out complete decay correction and calculation on air and a deposition concentration field output by an atmospheric diffusion model. The correction method comprises the following steps:
And 102, splitting the data of the radioactive release source item into M sub-release segments according to a first preset duration, wherein the first preset duration is positive integer multiple of delta t, delta t is the minimum time interval of the calculation result of the atmospheric diffusion model, and M is a positive integer.
Specifically, for a release source item with a total release duration D, taking a first preset duration D, splitting the release source item into M sub-release segments, where m=d/D (i.e., the value of M is the total release duration divided by the first preset duration), where the M sub-release segment has an actual release rate value only in the [ (M-1) D, md ] time period, and no pollutant is released in other time periods, or the release rate is expressed as 0 in other time periods, where M is a positive integer, and 1 is less than or equal to M. It should be noted that, when the release source item contains a plurality of radionuclides (nuclides for short), each nuclide needs to be split according to the method, and the numerical value of each sub-release section after the splitting is irrelevant to the nuclide type. When the value of the first preset duration d is larger, the total number M of the sub-release segments is smaller, the calculation cost is lower, and the accuracy of decay correction is correspondingly reduced.
In this embodiment, the diffusion result of the M sub-release sections obtained by calculating the meteorological data and the atmospheric diffusion model may be obtained in two ways, one is that the meteorological data and the atmospheric diffusion model are set first, and then the diffusion result is obtained by calculating the meteorological data and the atmospheric diffusion model, and the other is that the diffusion result obtained by calculating the external meteorological data and the atmospheric diffusion model is directly obtained. The atmospheric diffusion model in the two acquisition modes can be various pollutant diffusion models established and tested by the atmospheric physics profession, and does not comprise the decay process of the radioactive pollutant, because the decay process of the radioactive pollutant is overlapped with the correction method provided by the embodiment, and the correction result is deviated. Wherein the diffusion result includes an air concentration field and a deposition flux field. Each sub-release segment has a diffusion result N group, where n=d/a second preset duration (i.e., N is the total release duration D divided by the second preset duration). If the cores of the computer processor are enough in calculation, when the values of the first preset duration d and the second preset duration d are deltat, n=m, and the error of decay correction is minimum, namely the accuracy of decay correction is highest. In addition, in order to improve the calculation efficiency, the diffusion results of all the sub-release sections are calculated in parallel, the total number of the parallel sub-release sections is M, and each sub-release section is provided with N groups of diffusion results with intervals of a second preset duration.
In this embodiment, since the diffusion results of the N groups of intervals output by the atmospheric diffusion model are time-varying and recorded at the nth time, in order to obtain the correction value of the air concentration field in the diffusion results, the decay multiple corresponding to the recorded time is required to be obtained, and therefore the decay multiple of the nth time corresponding to the N groups of diffusion results of each sub-release section is required to be calculated. For example, when n=1, the decay multiple of the last time of group 1 (i.e., time 1) needs to be calculated corresponding to the diffusion result of group 1, and when n=5, the decay multiple of the last time of group 5 (i.e., time 5) needs to be calculated corresponding to the diffusion result of group 5.
Optionally, the method of calculating the decay times in this embodiment is as follows: and obtaining the decay chain information of all nuclides according to the nuclide types in the release source item data, and regarding each nuclide as a parent nucleus according to the decay chain information and respectively constructing a decay tree. And then according to the decay tree of each parent nucleus, adopting depth-first search, parallelly calculating the decay times of each parent nucleus and each child nucleus at the nth moment corresponding to N groups of diffusion results of each child release section respectively, and storing the decay times into the decay tree. Wherein the parent nucleus is the radionuclide contained in the source item, and the child nucleus is the nuclide newly generated by the decay of the parent nucleus. Of course, other known ways of calculating decay times may be used.
The decay tree data structure is provided in the embodiment, and the multi-generation decay chain information of the radionuclide can be completely stored in a single linked list mode without any decay chain interception or simplification, so that the problems of numerous nuclides and complex and tedious decay chain structures are solved, and support is provided for decay multiple calculation. Meanwhile, the decay tree structure can simplify a complex decay chain structure, and the decay chain can be quickly traversed by depth-first search, so that the efficiency of multi-generation decay calculation can be improved. In addition, the decay times of the parent nucleus and the child nucleus at the nth moment corresponding to N groups of diffusion results of each child release segment are calculated in parallel, so that the calculation efficiency can be further improved.
Specifically, each nuclide is regarded as a parent nucleus and decay trees are respectively constructed, and the method specifically comprises the following steps:
step 1041, obtaining multi-generation decay chain information of a parent nucleus, and storing each nuclide in the decay chain information as a node of a decay tree, wherein the parent nucleus is a root node of the decay tree, and the child nucleus is a child node of the decay tree.
Step 1042, storing all nodes in the decay tree as nodes in a single-linked list, wherein each node stores node information, and the node information comprises the name of the nuclide, the decay constant, the decay branch ratio and other physical properties, the position of the nuclide in the single-linked list, and the position of each child core of the nuclide in the single-linked list.
Step 1043, constructing a directed relationship of all nodes according to the positions of the nuclides in the unidirectionally linked list and the positions of the sub-cores of the nuclides in the unidirectionally linked list, so as to form a decay tree.
In this embodiment, N groups of diffusion results of M sub-release sections are read in parallel, and then the product of the air concentration field in the diffusion result at the nth time of each sub-release section and the decay multiple corresponding to the nth time is calculated in parallel, where the total number of parallel calculations is M, and the calculation efficiency can be improved by parallel reading the diffusion results and calculating the product of the air concentration field and the decay multiple in parallel.
The radioactive release source item is split into M sub-release sections, so that the complete decay process is supplemented, the concentration of parent nuclei is reduced, and the concentration of the child nuclei is added, so that the air concentration field correction value at the nth moment of the M sub-release sections is required to be summed to obtain the air concentration field correction value at the nth moment of the release source item. The air concentration field correction value of the radioactive emission source item is finally obtained by performing the correction operation and the summation operation on the air concentration fields at the N times respectively for the step 105 and the step 106.
Optionally, after obtaining the air concentration field correction value of the radioactive emission source item, the correction method of the atmospheric diffusion result of the nuclear accident radioactive substance further comprises calculating a deposition flux field correction value, wherein the specific calculation steps are as follows:
performing a correction operation on the deposition flux field at the nth time of each sub-release segment: dividing the obtained deposition flux field in the diffusion result of the nth time of each sub-release section by the air concentration field of the time, and multiplying the obtained deposition flux field by the air concentration field correction value of the time to obtain the deposition flux field correction value of the nth time of each sub-release section. A summation operation is then performed on the deposited flux fields at the nth time of the M sub-release segments: and summing the deposition flux field correction values at the nth time of the M sub-release segments to obtain the deposition flux field correction value at the nth time of the radioactive release source item. By performing the above calculation process N times, that is, performing correction operation and summation operation on the deposition flux fields at N times, respectively, N sets of deposition flux field correction values of the radioactive emission source term can be obtained.
Optionally, after obtaining the correction value of the deposition flux field of the radioactive emission source item, the correction method of the atmospheric diffusion result of the nuclear accident radioactive substance further comprises calculating a ground deposition concentration field:
performing decay correction with decay time of the ground deposition concentration field at the N-1 time being a second preset duration to obtain contribution from the ground deposition concentration field at the N-1 time in the ground deposition concentration field at the N-1 time, wherein N is a positive integer, and N is more than or equal to 1 and less than or equal to N; multiplying the obtained corrected value of the deposition flux field at the nth moment of the radioactive release source item by a second preset time length to obtain the contribution of air deposition to the ground in the [ n-1, n ] period, and summing the contribution of the air deposition to the ground from the ground deposition concentration field at the nth moment and the contribution of the air deposition from the ground deposition concentration field at the nth moment to obtain the ground deposition concentration field at the nth moment.
In this embodiment, when n=1, that is, the ground deposition concentration field at time 0 is 0, since the ground deposition concentration field at time 1 is the product of the deposition flux field correction value at time 1 of the radioactive emission source item and the second preset time period because the ground deposition concentration field at time 0 has not yet been diffused.
Example 2:
the embodiment provides a nuclear accident consequence evaluation method, which comprises the following steps:
Step 201, calculating an effective dose to be applied to the public according to the air concentration field correction value, the ground deposition concentration field, and the dose conversion factors corresponding to the air concentration of each nuclide, the ground deposition concentration, and the inhalation amount of human body from air of the radioactive release source item obtained in the embodiment 1, wherein the effective dose includes an air immersion external irradiation dose, a ground deposition external irradiation dose, and an inhalation internal irradiation dose in an irradiated period. The air immersion external irradiation dose technology considers dose irradiation of a three-dimensional space to a certain point on the ground, and the ground deposition external irradiation dose does not consider inert gas.
Step 202, comparing the effective dose with the general optimized intervention level, and determining the protection action according to the comparison result. Since the corrected diffusion result is more accurate, the result evaluation result calculated by the correction method in embodiment 1 is more accurate, thereby providing a safety protection suggestion for the public.
In this embodiment, the irradiation time period is 2 days, 7 days, 1 month, 1 year, etc. corresponding to the general optimized intervention level in the national standard GB18871-2002, and the calculated irradiation dose is the cumulative dose in the irradiation time period, including the effective dose and the organ dose.
Specifically, if T 2d >10mSv, the concealing action should be adopted; if T 7d Is 50mSv, evacuation action should be employed, where T 2d Represents the total effective dose for 2 days, T 7d Representing the total effective dose for 7 days.
The beneficial effects of the embodiment are as follows:
(1) The radioactive release source item is split into M sub-release sections according to a first preset duration, wherein the first preset duration is positive integer times of delta t, and delta t is the minimum time interval of the calculation result of the atmospheric diffusion model. To achieve the integration of the diffusion process into zero, so as to simulate the real diffusion process that the diffusion concentration at a certain place is generated by superposition of a plurality of different release sections.
(2) Multiplying the obtained air concentration field at the nth time of each sub-release section by the decay multiple corresponding to the time to obtain an air concentration field correction value at the nth time of each sub-release section, summing the air concentration field correction values at the nth time of the M sub-release sections to obtain an air concentration field correction value at the nth time of the radioactive release source item, and calculating a correction value of a deposition flux field according to the air concentration field correction value. And executing the calculation for N times to finish correction of the diffusion result of the existing atmospheric diffusion model, so that the diffusion result is more accurate.
(3) The correction method of the nuclear accident radioactive substance atmospheric diffusion result can directly carry out complete decay correction on the air concentration field and the deposition flux field output by the atmospheric diffusion model, and secondary development on the existing atmospheric diffusion model is not needed, so that the correction method has strong applicability.
(4) And taking each nuclide releasing the source item as a parent nucleus, respectively constructing a decay tree for each parent nucleus, traversing a decay chain by adopting depth-first search, and improving the efficiency of multi-generation decay calculation. And the calculation efficiency can be improved by parallel calculating the diffusion result, parallel calculating the decay multiple and parallel reading the diffusion result of the sub-release section and parallel correcting the diffusion result of the sub-release section.
(5) Because the corrected diffusion result is more accurate, the effective dosage is calculated through the correction value of the air concentration field and the ground deposition concentration field, and the effective protection action is obtained according to the effective dosage, so that the public safety can be ensured.
As shown in fig. 2-5, the method for evaluating the consequences of a nuclear accident in this embodiment is implemented by a specific example, in which the release source item includes only one nuclide I-135, and the release duration is 2 hours = 7200 seconds, and the method includes the following steps:
In step 301, a first preset duration d=1200 seconds is set, and d is taken as a time interval, so that the release source item is split into M sub-release segments, i.e. m=2h/1200 s=6.
In step 303, the decay times of the parent and child nuclei are calculated.
First, the source nuclide I-135 is regarded as a parent nucleus, the information of the previous generation decay chain is selected, and the decay tree of the parent nucleus I-135 is constructed according to the information of the decay chain.
The primary core I-135 includes the child cores Xe-135m, cs-135 and Xe-135 in the information of the third generation decay chain. As shown in FIG. 2, in the structure of the uni-directional linked list of the storage decay tree, each nuclide is stored as a node, wherein the number before the nuclide name is the position of the nuclide node in the uni-directional linked list, for example, the number "0" before the I-135 indicates that the position of the I-135 in the uni-directional linked list is the start bit, in each node, an indefinite length array of the positions of the child cores in the uni-directional linked list is also stored after the nuclide name for quick positioning and searching, for example, the numbers "1" and "3" in the array of the 0 th node I-135 indicate that the child core of the parent core I-135 comprises Xe-135m of the 1 st node and Xe-135 of the 3 rd node respectively.
The decay tree structure which is completely stored by adopting the single-linked list structure is shown in figure 3, wherein the 0 th node of the single-linked list represents a parent nucleus, the parent nucleus is the root node of the decay tree, the rest nodes of the single-linked list represent child cores, the child cores are child nodes of the decay tree, and the directed relation of all the nodes is constructed according to the positions of nuclide nodes in the single-linked list and the positions of all the child cores of the nuclide in the single-linked list so as to form the decay tree, so that the decay chain can be traversed quickly by using a depth-first search method, and the calculation efficiency of decay multiples is improved.
Second, since the first preset duration d=2Δt=1200 seconds and the second preset duration is equal to Δt, which is 600 seconds, the decay time of each diffusion result of each sub-release section is shown in fig. 4. Wherein the decay times are all integer multiples of Δt, N is 1 to N-2, and in addition, since the m=1st sub-release segment is released within the initial 2Δt time, and overlaps with the time period 2Δt of the 1 st time (n=1) and the 2 nd time (n=2), the same conservative treatment is performed on the other sub-release segments in this embodiment, irrespective of the decay of the sub-release segment during the discharge.
Along the decay tree of parent core I-135, decay times of N Δt are calculated and stored into the decay tree, where n=1, 2. Specifically, from half-life and decay branch ratio, the concentration change times (decay times) of I-135 itself and its daughter nuclei after decay for 600s, 1200s, 1800s, … …, 6000s were calculated. By adopting depth-first search, decay trees shown in fig. 3 are searched downwards along the left side, and after searching for the leftmost branch, the decay trees are searched rightwards step by step, wherein the specific sequence is as follows:
(1) Checking the parent nucleus I-135, and respectively calculating decay multiples after the decay time to be calculated according to the half-life period;
(2) Finding the first child nucleus Xe-135m on the left side, and respectively calculating decay times after decay time to be calculated according to the half-life of the parent nucleus I-135, the half-life of the parent nucleus and the branching ratio of the decay of the parent nucleus I-135 into the Xe-135 m;
(3) Finding the left extreme nuclide Cs-135, and calculating decay times from the half-lives of the parent nucleus I-135, the parent nucleus Xe-135m and the parent nucleus thereof and the decay branch ratio between the parent nucleus and the parent nucleus; since Cs-135 decays without producing new daughter nuclei, the leftmost branch of the decay tree structure reaches the end;
(4) After returning Xe-135m, finding the other child core Xe-135 to the right side, and calculating the decay times along I-135- > Xe-135m- > Xe-135;
(5) Finding out the child core Cs-135 of Xe-135 to reach the tail end of the second branch of the decay tree structure; calculating the decay factor of Cs-135 on the branch along the line I-135- > Xe-135m- > Xe-135- > Cs-135, and superposing the decay factor calculated in the step (3) to update the decay factor of Cs-135;
(6) After returning Xe-135, only the branch detected in the step (5) is found, the branch returns upwards to Xe-135m, and both branches are found in the step (3) and the step (5), so that the branch continues to return upwards to the tree root I-135, and the other child core Xe-135 is found to the right; calculating the decay times of the Xe-135 on the branch along the length I-135- > Xe-135, and superposing the decay times calculated in the step (4) to update the decay times of the Xe-135;
(7) Finding out the child core Cs-135 of Xe-135 to reach the tail end of the third branch of the decay tree structure; calculating decay times of the Cs-135 on the branch along the I-135- > Xe-135- > Cs-135 and superposing the result of the step (5) for updating; at this point all branches of the decay tree structure are traversed.
For example, for a decay time of 600s, the decay times of I-135, xe-135m, cs-135 and Xe-135 are respectively: 0.98, 0.01, 3.18E-14 and 0.06; for decay times of 6000s, the decay times of I-135, xe-135m, cs-135 and Xe-135 are respectively: 0.84, 0.105, 3.14E-12 and 0.14.
And 304, carrying out parallel correction, and circularly carrying out N=12 times to calculate to obtain 12 groups of diffusion result correction values of the release source item and a ground deposition concentration field.
Firstly, reading diffusion results of each sub-release section output by an atmospheric diffusion model for M times in parallel, correcting an air concentration field in an nth set of results of each sub-release section for M times in parallel according to decay time shown in fig. 4 to obtain an air concentration field correction value at an nth time of each sub-release section (wherein the nth time is an end time of the nth set of diffusion results), and then superposing the air concentration field correction values at the nth time of all M sub-release sections to obtain an air concentration field correction value at the nth time of a release source item. And then carrying out decay correction with the decay time delta t on the ground deposition concentration field, and accumulating the deposition flux field correction value at the nth moment to the ground deposition concentration field.
For example, as for n=1, as can be seen from fig. 4, the release of the m=1 th sub-release segment is not ended, the conservative consideration does not carry out decay correction on the diffusion result, the m=2 to m=6 th sub-release segments have no radioactivity release yet, and the diffusion result is all 0. Since the air concentration field correction value at the n=1 th time is the sum of the air concentration field correction values of the m=1 to m=6 th sub-release sections, the value thereof is equal to the air concentration field of the m=1 th sub-release section. At this time, the obtained deposition flux field in the diffusion result at the n=1th time in the m=1 th sub-release section is divided by the air concentration field at the time, and then multiplied by the air concentration field correction value at the n=1th time to obtain the deposition flux field correction value at the n=1th time in the m=1 th sub-release section, and the deposition flux field correction value is multiplied by Δt to obtain the concentration field accumulated and deposited on the ground in the period, namely the ground deposition concentration field at the n=1th time.
For n=n=12, as can be seen from fig. 4, the diffusion results of the other sub-release sections except the m=m=6 sub-release section need to be subjected to decay correction, and the decay times are respectively 10 times, 8 times, 6 times, 4 times and 2 times Δt (6000 s, 4800s, 3600s, 2400s, 1200 s). The air concentration field correction value at the n=12 th time is: the air concentration fields of the m=1 to m=5 sub-release sections are multiplied by corresponding decay factors and added, and the air concentration field of the m=6 sub-release section is added to obtain an air concentration field correction value of the nth group, namely the 12 th moment. And (3) carrying out decay correction with the decay time delta t on the n=11 ground deposition concentration field, and multiplying the deposition flux field correction value at the n=12 moment by delta t to accumulate the deposition flux field correction value to the ground deposition concentration field to obtain the ground deposition concentration field at the 12 th moment. The deposition flux field correction value at the n=12 moment is obtained by dividing the deposition flux field at the moment in each sub-release section by the air concentration field at the moment, multiplying the divided deposition flux field with the air concentration field correction value at the moment to obtain the deposition flux field correction value at the n=12 moment in each sub-release section, and summing the deposition flux field correction values at the n=12 moment in the M sub-release sections.
For n=2, 3 …, the calculation process is similar to that described above, and will not be repeated here.
And 305, multiplying the air concentration field correction values of the I-135, the Xe-135m, the Cs-135 and the Xe-135 and the ground deposition concentration field by the respective dose conversion factors DCF (Dose Conversion Factor) respectively and superposing to obtain an effective dose field. For example, suppose I-135, xe-135m, cs-135 and Xe-135 ground floor air concentration c air,1 、c air,2 、c air,3 And c air,4 The ground deposition concentration is c sur,1 、c sur,2 、c sur,3 And c sur,4 The effective dose of DCF irradiated outside air immersion is d respectively air,1 、d air,2 、d air,3 And d air,4 The effective dose of DCF of external irradiation for ground deposition is d respectively sur,1 、d sur,2 、d sur,3 And d sur,4 The effective dose of DCF for inhalation irradiation is d respectively inh,1 、d inh,2 、d inh,3 And d inh,4 When the respiration rate is h, the calculation formula of the effective dose of the inhaled internal radiation of the public at the grid point is shown as formula (1): (c) air,1 d air,1 +c air,2 d air,2 +c air, 3 d air,3 +c air,4 d air,4 ) Δt. (1) The calculation formula of the inhaled internal irradiation effective dose on the grid points is shown as formula (2): (c) sur,1 d sur,1 +c sur,2 d sur,2 +c sur,3 d sur,3 +c sur,4 d sur,4 ) Δt. (2) The calculation formula of the inhaled internal irradiation effective dose on the grid points is shown as formula (3): (c) air,1 d inh,1 +c air,2 d inh,2 +c air,3 d inh,3 +c air, 4 d inh,4 )*h*Δt。 (3)
If the effective dosages of the three routes for 2 days are A, B, C respectively, the total effective dosage T for 2 days 2d ＝A+B+C。
According to the national standard, if T 2d >10mSv, where T should suggest that personnel at the grid point take concealing action 2d Represents a total effective dose for 2 days; if T 7d >50mSv, where T should suggest that personnel at the grid point take evacuation action 7d Representing the total effective dose for 7 days.
Example 3:
as shown in fig. 6, the present embodiment provides a device for correcting the atmospheric diffusion result of a nuclear accident radioactive substance, which includes a splitting module 41, an obtaining module 42, a calculating module 43 and a first correcting module 44.
The splitting module 41 is configured to acquire radioactive release source item data, and split the radioactive release source item data into M sub-release segments according to a first preset duration, where the first preset duration is a positive integer multiple of Δt, and Δt is a minimum time interval of an atmospheric diffusion model calculation result.
The obtaining module 42 is connected with the splitting module 41, and is configured to obtain diffusion results of M sub-release sections obtained by calculating meteorological data and an atmospheric diffusion model, and make each sub-release section have N groups of diffusion results with a time interval being a second preset duration, where the second preset duration is a positive integer multiple of Δt, and N is an integer multiple of M; the diffusion result includes an air concentration field.
The calculating module 43 is connected to the obtaining module 42, and is configured to calculate decay times of N-th moments corresponding to N groups of diffusion results of each sub-release segment, where N is a positive integer, and N is greater than or equal to 1 and less than or equal to N.
Specifically, the calculation module 43 includes a decay multiple unit, where the decay multiple unit is connected to the splitting module 41 and the obtaining module 42, and is configured to obtain decay chain information according to all nuclides in the release source item obtained by the splitting module 41, and treat each nuclide as a parent nucleus according to the decay chain information and respectively construct a decay tree; and the method is used for parallelly calculating the decay times of the parent nucleus and the child nucleus at the nth moment respectively corresponding to N groups of diffusion results of each child release section by adopting depth-first search according to the decay tree of each parent nucleus, and storing the decay times into the decay tree.
The first correction module 44 is connected to the acquisition module 42 and the calculation module 43, and is configured to perform a correction operation on the air concentration field at the nth time of each sub-release section: the air concentration field correction value is obtained by multiplying the air concentration field in the diffusion result of the nth moment of each sub-release section obtained by the obtaining module 42 by the decay multiple corresponding to the moment calculated by the calculating module 43; and performing a summation operation on the air concentration field at the nth time of the M sub-release sections: the method comprises the steps of summing air concentration field correction values at the nth moment of M sub-release sections to obtain air concentration field correction values at the nth moment of a radioactive release source item; and the air concentration field correction value of the radioactive release source item is obtained by respectively performing correction operation and summation operation on the air concentration fields at N moments.
Optionally, the correction device for the atmospheric diffusion result of the nuclear accident radioactive material further comprises a second correction module, and the diffusion result further comprises a deposition flux field, and the second correction module is connected to the acquisition module 42, the calculation module 43, and the first correction module 44, and is configured to perform a correction operation on the deposition flux field at the nth time of each sub-release segment: the method comprises the steps of dividing a deposition flux field in a diffusion result of the nth moment of each sub-release section obtained by the obtaining module 42 by an air concentration field of the moment, and multiplying the deposition flux field by an air concentration field correction value of the nth moment obtained by the first correction module 44 to obtain a deposition flux field correction value of the nth moment of each sub-release section; and performing a summation operation on the deposited flux fields at the nth time of the M sub-release segments: the method comprises the steps of summing the deposition flux field correction values at the nth moment of the M sub-release sections to obtain the deposition flux field correction value at the nth moment of the radioactive release source item; and the correction operation and the summation operation are respectively carried out on the deposition flux fields at N moments, so as to obtain the deposition flux field correction value of the radioactive release source item.
Optionally, the correction device for the atmospheric diffusion result of the nuclear accident radioactive substance further comprises a third correction module, wherein the third correction module is connected with the acquisition module 42 and the second correction module and is used for carrying out decay correction on the ground deposition concentration field at the N-1 time as a second preset duration to obtain the contribution quantity from the ground deposition concentration field at the N-1 time in the ground deposition concentration field at the N-1 time, wherein N is a positive integer, and N is more than or equal to 1 and less than or equal to N; and the correction value of the deposition flux field at the nth time of the radioactive release source item obtained by the second correction module is multiplied by the second preset duration of the acquisition module 42, and then summed with the contribution from the ground deposition concentration field at the (n-1) th time in the ground deposition concentration field at the nth time to obtain the ground deposition concentration field at the nth time.
Example 4:
the embodiment provides a nuclear accident consequence evaluation system, including: the device for correcting the atmospheric diffusion result of the nuclear accident radioactive material in embodiment 3 further comprises a dosage module and an evaluation module.
The dose module is connected with the correction device of the nuclear accident radioactive substance atmospheric diffusion result and is used for calculating the effective dose received by the public according to the dose conversion factors corresponding to the air concentration of each nuclide, the ground deposition concentration and the inhalation quantity of human body from air respectively, and the air concentration field correction value and the ground deposition concentration field of the radioactive release source item obtained by the correction device, wherein the effective dose comprises the air immersion external irradiation dose, the ground deposition external irradiation dose and the inhalation internal irradiation dose in the illuminated time period.
And the evaluation module is used for storing the general optimized intervention level, pre-storing a mapping table of the comparison result and the protection action, connecting the comparison result and the dose module, comparing the general optimized intervention level with the effective dose calculated by the dose module, and determining the corresponding protection action according to the mapping table according to the comparison result.
It is to be understood that the above embodiments are merely illustrative of the application of the principles of the present invention, but not in limitation thereof. Various modifications and improvements may be made by those skilled in the art without departing from the spirit and substance of the invention, and are also considered to be within the scope of the invention.
Claims (9)
1. A method for correcting the atmospheric diffusion result of a nuclear accident radioactive substance, comprising:
acquiring radioactive release source item data;
dividing the radioactive release source item data into M sub-release sections according to a first preset duration, wherein the first preset duration is positive integer multiple of delta t, and delta t is the minimum time interval of the calculation result of the atmospheric diffusion model;
obtaining diffusion results of M sub-release sections obtained through calculation of meteorological data and an atmospheric diffusion model, and enabling each sub-release section to have N groups of diffusion results with time intervals of a second preset duration, wherein the second preset duration is positive integer multiple of deltat, N is integer multiple of M, and the diffusion results comprise an air concentration field;
calculating decay multiples of the nth moment corresponding to N groups of diffusion results of each sub-release section respectively, wherein N is a positive integer, and N is more than or equal to 1 and less than or equal to N;
performing a correction operation on the air concentration field at the nth time of each sub-release section: multiplying the obtained air concentration field in the diffusion result at the nth moment of each sub-release section by the decay multiple corresponding to the moment to obtain an air concentration field correction value at the nth moment of each sub-release section;
performing a summation operation on the air concentration field at the nth time of the M sub-release sections: summing the air concentration field correction values at the nth moment of the M sub-release sections to obtain the air concentration field correction value at the nth moment of the radioactive release source item;
Respectively carrying out correction operation and summation operation on the air concentration fields at N moments to obtain an air concentration field correction value of the radioactive release source item,
the calculating the decay multiple of the nth time corresponding to the N groups of diffusion results of each sub-release section specifically comprises the following steps:
according to decay chain information of all nuclides in the radioactive release source item data, each nuclide is regarded as a parent nucleus and decay trees are respectively constructed;
according to the decay tree of each parent nucleus, depth-first search is adopted, the decay times of each parent nucleus and child nucleus at the nth moment corresponding to N groups of diffusion results of each child release segment are calculated in parallel and stored into the decay tree,
the method for constructing the decay tree by using each nuclide as a parent nucleus comprises the following steps of:
obtaining multi-generation decay chain information of a parent nucleus, and storing each nuclide in the decay chain information as a node of a decay tree, wherein the parent nucleus is a root node of the decay tree, and the child nucleus is a child node of the decay tree;
storing all nodes in the decay tree as nodes in a single-linked list, wherein each node stores node information, and the node information comprises the name of the nuclide, decay constant and decay branch ratio, the position of the nuclide in the single-linked list and the position of each sub-core of the nuclide in the single-linked list;
And constructing the directed relation of all nodes according to the positions of the nuclides in the unidirectionally linked list and the positions of the child cores of the nuclides in the unidirectionally linked list so as to form the decay tree.
2. The method for modifying the atmospheric diffusion result of a nuclear accident radioactive substance according to claim 1, wherein the diffusion result further comprises a deposition flux field,
after obtaining the air concentration field correction value for the radioactive emission source term, the method further comprises:
performing a correction operation on the deposition flux field at the nth time of each sub-release segment: dividing the obtained deposition flux field in the diffusion result at the nth time of each sub-release section by the air concentration field at the time, and multiplying the obtained deposition flux field by the air concentration field correction value at the time to obtain the deposition flux field correction value at the nth time of each sub-release section;
performing a summation operation on the deposited flux fields at the nth time of the M sub-release segments: summing the deposition flux field correction values at the nth time of the M sub-release sections to obtain the deposition flux field correction value at the nth time of the radioactive release source item;
and respectively performing correction operation and summation operation on the deposition flux fields at N moments to obtain the deposition flux field correction value of the radioactive release source item.
3. The method of claim 2, wherein after obtaining the deposition flux field correction value for the radioactive emission source term, the method further comprises:
performing decay correction with decay time of the ground deposition concentration field at the N-1 time being a second preset duration to obtain contribution from the ground deposition concentration field at the N-1 time in the ground deposition concentration field at the N-1 time, wherein N is a positive integer, and N is more than or equal to 1 and less than or equal to N;
multiplying the obtained deposition flux field correction value at the nth time of the radioactive release source item by a second preset time length to obtain the contribution of air deposition to the ground in the [ n-1, n ] period, and summing the contribution of the air deposition to the ground from the nth time to obtain the ground deposition concentration field at the nth time.
4. A method for evaluating the consequences of a nuclear accident, comprising:
calculating an effective dose to be received by the public according to the air concentration field correction value, the ground deposition concentration field and dose conversion factors corresponding to the air concentration of each nuclide, the ground deposition concentration and the inhalation amount of human body from air of the radioactive release source item obtained in claim 3, wherein the effective dose comprises an air immersed external irradiation dose, a ground deposition external irradiation dose and an inhalation internal irradiation dose in an irradiated time period;
Comparing the effective dose with the general optimized intervention level, and determining the protective action according to the comparison result.
5. The method of claim 4, wherein comparing the effective dose to the general optimized intervention level and determining the protective action based on the comparison result comprises:
if T 2d >10mSv, adopting concealing action; if T 7d 50mSv, with evacuation action, where T 2d Represents the total effective dose for 2 days, T 7d Representing the total effective dose for 7 days.
6. The device for correcting the atmospheric diffusion result of the nuclear accident radioactive substance is characterized by comprising a splitting module, an acquisition module, a calculation module and a first correction module,
the splitting module is used for acquiring the radioactive release source item data and splitting the radioactive release source item data into M sub-release sections according to a first preset duration, wherein the first preset duration is positive integer multiple of Deltat, and Deltat is the minimum time interval of the calculation result of the atmospheric diffusion model;
the acquisition module is connected with the splitting module and is used for acquiring diffusion results of M sub-release sections obtained through calculation of meteorological data and an atmospheric diffusion model, and enabling each sub-release section to have N groups of diffusion results with time intervals of a second preset duration, wherein the second preset duration is positive integer multiple of deltat, N is integer multiple of M, and the diffusion results comprise an air concentration field;
The calculation module is connected with the acquisition module and is used for calculating decay multiples of the nth moment corresponding to N groups of diffusion results of each sub-release section respectively, wherein N is a positive integer, and N is more than or equal to 1 and less than or equal to N;
the first correction module is connected with the acquisition module and the calculation module and is used for executing correction operation on the air concentration field at the nth moment of each sub-release section: the air concentration field in the diffusion result of the nth moment of each sub-release section obtained by the obtaining module is multiplied by the decay multiple corresponding to the moment calculated by the calculating module, so that an air concentration field correction value of the nth moment of each sub-release section is obtained;
and performing a summation operation on the air concentration field at the nth time of the M sub-release sections: the air concentration field correction value at the nth moment of the M sub-release sections is summed to obtain the air concentration field correction value at the nth moment of the radioactive release source item;
and the air concentration field correction value of the radioactive release source item is obtained by respectively carrying out correction operation and summation operation on the air concentration fields at N moments,
wherein the calculation module comprises a decay multiple unit,
the decay multiple unit is connected with the splitting module and the acquisition module and is used for acquiring decay chain information according to all nuclides in the release source item acquired by the splitting module, treating each nuclide as a parent nucleus according to the decay chain information and respectively constructing a decay tree;
And the decay multiples of the parent nucleus and the child nucleus at the nth moment corresponding to N groups of diffusion results of each child release section are calculated in parallel by adopting depth-first search according to the decay tree of each parent nucleus and storing the decay multiples to the decay tree,
specifically, the decay multiple unit is used for acquiring multi-generation decay chain information of a parent nucleus, and storing each nuclide in the decay chain information as a node of a decay tree, wherein the parent nucleus is a root node of the decay tree, and the child nucleus is a child node of the decay tree; the method is also used for storing all node arrangements in the decay tree as nodes in the unidirectional linked list, and each node stores node information, wherein the node information comprises the name of the nuclide, decay constant and decay branch ratio, the position of the nuclide in the unidirectional linked list and the position of each child core of the nuclide in the unidirectional linked list; and constructing a directed relation of all nodes according to the positions of the nuclides in the unidirectionally linked list and the positions of the sub-cores of the nuclides in the unidirectionally linked list so as to form the decay tree.
7. The apparatus for modifying the atmospheric diffusion of nuclear accident radioactive materials according to claim 6, further comprising a second modifying module, the diffusion result further comprising a deposition flux field,
The second correction module is connected with the acquisition module, the calculation module and the first correction module and is used for performing correction operation on the deposition flux field at the nth moment of each sub-release section: the method comprises the steps of dividing a deposition flux field in a diffusion result at the nth moment of each sub-release section obtained by an obtaining module by an air concentration field at the moment, and multiplying the deposition flux field by an air concentration field correction value at the nth moment obtained by a first correction module to obtain a deposition flux field correction value at the nth moment of each sub-release section;
and performing a summation operation on the deposited flux fields at the nth time of the M sub-release segments: the method comprises the steps of summing the deposition flux field correction values at the nth moment of the M sub-release sections to obtain the deposition flux field correction value at the nth moment of the radioactive release source item;
and the correction operation and the summation operation are respectively carried out on the deposition flux fields at N moments, so as to obtain the deposition flux field correction value of the radioactive release source item.
8. The apparatus for correcting atmospheric diffusion results of nuclear accident radioactive materials according to claim 7, further comprising a third correction module,
the third correction module is connected with the acquisition module and the second correction module and is used for carrying out decay correction on the ground deposition concentration field at the N-1 time as a second preset duration to obtain the contribution from the ground deposition concentration field at the N-1 time in the ground deposition concentration field at the N-1 time, wherein N is a positive integer, and N is more than or equal to 1 and less than or equal to N;
And the deposition flux field correction value at the nth moment of the radioactive release source item obtained by the second correction module is multiplied by a second preset duration and then summed with the contribution from the ground deposition concentration field at the nth-1 moment to obtain the ground deposition concentration field at the nth moment.
9. A nuclear accident outcome evaluation system, comprising: the nuclear accident radioactive material atmospheric diffusion result correction apparatus of claim 8, further comprising a dosage module and an evaluation module:
the dose module is connected with the correction device of the nuclear accident radioactive substance atmospheric diffusion result and is used for calculating the effective dose received by the public according to the dose conversion factors corresponding to the air concentration of each nuclide, the ground deposition concentration and the inhalation quantity of human body from air respectively, and the air concentration field correction value and the ground deposition concentration field of the radioactive release source item obtained by the correction device, wherein the effective dose comprises the air immersion external irradiation dose, the ground deposition external irradiation dose and the inhalation internal irradiation dose in the illuminated time period;
and the evaluation module is used for storing the general optimized intervention level, pre-storing a mapping table of the comparison result and the protection action, connecting the comparison result and the dose module, comparing the general optimized intervention level with the effective dose calculated by the dose module, and determining the corresponding protection action according to the mapping table according to the comparison result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010355824.6A CN111695762B (en) | 2020-04-29 | 2020-04-29 | Correction method and device for nuclear accident diffusion result and result evaluation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010355824.6A CN111695762B (en) | 2020-04-29 | 2020-04-29 | Correction method and device for nuclear accident diffusion result and result evaluation method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111695762A CN111695762A (en) | 2020-09-22 |
| CN111695762B true CN111695762B (en) | 2023-05-05 |
Family
ID=72476804
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010355824.6A Active CN111695762B (en) | 2020-04-29 | 2020-04-29 | Correction method and device for nuclear accident diffusion result and result evaluation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111695762B (en) |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20080040352A (en) * | 2006-11-03 | 2008-05-08 | 한국전력공사 | Risk-informed population hazard assessment |
| WO2010041192A2 (en) * | 2008-10-10 | 2010-04-15 | Koninklijke Philips Electronics N.V. | Practical spect calibration method for quantification of nuclides with high-energy contributions |
| JP2013088206A (en) * | 2011-10-14 | 2013-05-13 | Mitsubishi Heavy Ind Ltd | Diffusion situation prediction system |
| CN103413062A (en) * | 2013-08-29 | 2013-11-27 | 中国测绘科学研究院 | Computing method of diffusion of radionuclides |
| CN103903106A (en) * | 2014-04-21 | 2014-07-02 | 苏州热工研究院有限公司 | Nuclear emergency decision supporting platform and method based on operation disturbance level |
| CN104217384A (en) * | 2013-05-29 | 2014-12-17 | 江苏省核应急办公室 | A nuclear accident emergency processing and auxiliary decision support system |
| CN104915744A (en) * | 2014-03-14 | 2015-09-16 | 江苏达科智能科技有限公司 | Nuclear emergency command method and command system |
| CN107145700A (en) * | 2016-03-01 | 2017-09-08 | 中国辐射防护研究院 | Core and radiation accident consequence airborne radioactivity dosages of substance evaluation method |
| CN108151967A (en) * | 2017-11-16 | 2018-06-12 | 中国核电工程有限公司 | A kind of nuclear power plant nuclear island workshop liquid leakage rate measuring device |
| CN108510114A (en) * | 2018-03-27 | 2018-09-07 | 环境保护部核与辐射安全中心 | Nucleic Population Doses From Medical prediction technique under a kind of nuclear power plant's future weather scene |
| CN110457829A (en) * | 2019-08-15 | 2019-11-15 | 王博 | A kind of source item release inverting and DIFFUSION PREDICTION method based on integrated model of atmospheric diffusion |
| CN110489789A (en) * | 2019-07-10 | 2019-11-22 | 哈尔滨工程大学 | Radgas spread evaluating method in a kind of elimination of nuclear facilities environment |
| CN110991809A (en) * | 2019-11-06 | 2020-04-10 | 中国辐射防护研究院 | Reactor core inventory real-time estimation method based on Hualong I |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7792774B2 (en) * | 2007-02-26 | 2010-09-07 | International Business Machines Corporation | System and method for deriving a hierarchical event based database optimized for analysis of chaotic events |
| FR2922667A1 (en) * | 2007-10-22 | 2009-04-24 | Commissariat Energie Atomique | METHOD FOR MANAGING A TIME-EVOLVING ACCIDENT |
| TW201220323A (en) * | 2010-11-02 | 2012-05-16 | Inst Nuclear Energy Res Atomic Energy Council | Parameter identification method for severe accidents |
| US9208911B2 (en) * | 2010-11-23 | 2015-12-08 | Westinghouse Electric Company Llc | Full spectrum LOCA evaluation model and analysis methodology |
-
2020
- 2020-04-29 CN CN202010355824.6A patent/CN111695762B/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20080040352A (en) * | 2006-11-03 | 2008-05-08 | 한국전력공사 | Risk-informed population hazard assessment |
| WO2010041192A2 (en) * | 2008-10-10 | 2010-04-15 | Koninklijke Philips Electronics N.V. | Practical spect calibration method for quantification of nuclides with high-energy contributions |
| JP2013088206A (en) * | 2011-10-14 | 2013-05-13 | Mitsubishi Heavy Ind Ltd | Diffusion situation prediction system |
| CN104217384A (en) * | 2013-05-29 | 2014-12-17 | 江苏省核应急办公室 | A nuclear accident emergency processing and auxiliary decision support system |
| CN103413062A (en) * | 2013-08-29 | 2013-11-27 | 中国测绘科学研究院 | Computing method of diffusion of radionuclides |
| CN104915744A (en) * | 2014-03-14 | 2015-09-16 | 江苏达科智能科技有限公司 | Nuclear emergency command method and command system |
| CN103903106A (en) * | 2014-04-21 | 2014-07-02 | 苏州热工研究院有限公司 | Nuclear emergency decision supporting platform and method based on operation disturbance level |
| CN107145700A (en) * | 2016-03-01 | 2017-09-08 | 中国辐射防护研究院 | Core and radiation accident consequence airborne radioactivity dosages of substance evaluation method |
| CN108151967A (en) * | 2017-11-16 | 2018-06-12 | 中国核电工程有限公司 | A kind of nuclear power plant nuclear island workshop liquid leakage rate measuring device |
| CN108510114A (en) * | 2018-03-27 | 2018-09-07 | 环境保护部核与辐射安全中心 | Nucleic Population Doses From Medical prediction technique under a kind of nuclear power plant's future weather scene |
| CN110489789A (en) * | 2019-07-10 | 2019-11-22 | 哈尔滨工程大学 | Radgas spread evaluating method in a kind of elimination of nuclear facilities environment |
| CN110457829A (en) * | 2019-08-15 | 2019-11-15 | 王博 | A kind of source item release inverting and DIFFUSION PREDICTION method based on integrated model of atmospheric diffusion |
| CN110991809A (en) * | 2019-11-06 | 2020-04-10 | 中国辐射防护研究院 | Reactor core inventory real-time estimation method based on Hualong I |
Non-Patent Citations (4)
| Title |
|---|
| 侯姗姗 ; .GIS及大气传输模型在核应急管理中的研究与系统实现.北京测绘.2016,(06),全文. * |
| 刘蕴 ; 方晟 ; 李红 ; 曲静原 ; 姚仁太 ; 范丹 ; .基于四维变分资料同化的核事故源项反演.清华大学学报(自然科学版).2015,(01),全文. * |
| 胡啸峰 ; 陈鹏 ; 曾昭龙 ; .基于WRF的放射性物质大气扩散模拟的环境辐射剂量估算方法.安全与环境工程.2016,(02),全文. * |
| 郭瑞萍 ; 张琼 ; 陈海英 ; 潘昕怿 ; 杨春林 ; .核电厂放射性气态流出物大气弥散模型现状与展望.核安全.2012,(02),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111695762A (en) | 2020-09-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107145698A (en) | Core and radiation accident consequence large scale analogy method and simulation system | |
| Rochman et al. | Uncertainties for Swiss LWR spent nuclear fuels due to nuclear data | |
| Huang et al. | Improvements to the Transmutation Trajectory Analysis of depletion evaluation | |
| Brown | A review of Monte Carlo criticality calculations–convergence, bias, statistics | |
| Leal et al. | Resonance parameter and covariance evaluation for $^{16} $ O up to 6 MeV | |
| CN111695762B (en) | Correction method and device for nuclear accident diffusion result and result evaluation method and system | |
| CN116028010A (en) | Radionuclide hazard sequencing method | |
| Abdurashitov et al. | The BNO–LNGS joint measurement of the solar neutrino capture rate in 71Ga | |
| Mlakar et al. | Relative doses instead of relative concentrations for the determination of the consequences of the radiological atmospheric releases | |
| Pane et al. | Preliminary Development of Radionuclides Release of Individual Dose Code Program for Radiation Monitoring Purposes | |
| Lloyd et al. | Individual external exposures from Nevada Test Site fallout for Utah leukemia cases and controls | |
| Mamikhin et al. | Mathematical model of 137Cs dynamics in the deciduous forest | |
| CN105373665A (en) | Method for obtaining multi-nuclide equivalent parameter used for nuclear station radiation simulation system | |
| Perevolotskaya et al. | Retrospective assessment of the formation of the radiation situation in pine plantations in the first year after the Chernobyl accident | |
| Jang et al. | Verification and Validation of Back-end Cycle Source Term Calculation of the Nodal Code RAST-K | |
| Williams et al. | Efficient Monte Carlo simulation of coincidence effects in radioisotope decays including γ–γ angular correlations | |
| Wu et al. | Three-dimensional dose evaluation system using real-time wind field information for nuclear accidents in Taiwan | |
| Zherdev et al. | ROCOCO: A constants supply system for Monte-Carlo reactor calculation | |
| Shaat et al. | Radiological impact due to atmospheric releases of the source term for FD, Unit 1, Using HOTSPOT code | |
| Perevolotsky et al. | Verification of a Migration–Dosimetric Model Describing the Formation of the Radiation Situation in Pine Stands in the First Year after Catastrophic Fallout | |
| Rovni et al. | High accuracy tritium measurement for the verification of the tritium production rate calculations with MCNPX | |
| Rao et al. | Parity dependence of the nuclear level density at high excitation | |
| Wilson et al. | NESHAPs CY19 Report: Radionuclide Air Emission Report for 2019 | |
| Grogan et al. | Reconstructing Reactor Operating Histories Using the INDEPTH Code | |
| Shanahan et al. | An Interdatabase Comparison of Nuclear Decay and Structure Data Utilized in the Calculation of Dose Coefficients for Radionuclides Produced in a Spallation Neutron Source |
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 |