CN111695762A - Correction method and device for nuclear accident diffusion result and consequence evaluation method and system - Google Patents

Correction method and device for nuclear accident diffusion result and consequence evaluation method and system Download PDF

Info

Publication number
CN111695762A
CN111695762A CN202010355824.6A CN202010355824A CN111695762A CN 111695762 A CN111695762 A CN 111695762A CN 202010355824 A CN202010355824 A CN 202010355824A CN 111695762 A CN111695762 A CN 111695762A
Authority
CN
China
Prior art keywords
decay
release
sub
nth
moment
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010355824.6A
Other languages
Chinese (zh)
Other versions
CN111695762B (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 Nuclear Power Engineering Co Ltd
Original Assignee
China Nuclear Power Engineering Co Ltd
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 Nuclear Power Engineering Co Ltd filed Critical China Nuclear Power Engineering Co Ltd
Priority to CN202010355824.6A priority Critical patent/CN111695762B/en
Publication of CN111695762A publication Critical patent/CN111695762A/en
Application granted granted Critical
Publication of CN111695762B publication Critical patent/CN111695762B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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; acquiring 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 respectively corresponding to the N groups of diffusion results of each sub-release section; multiplying the nth moment air concentration field of each sub-release section by the decay multiple of the moment to obtain the nth moment air concentration field correction value of the sub-release section; summing the corrected values of the air concentration field at the nth moment of the M sub-release sections to obtain the corrected value of the air concentration field at the nth moment; and traversing N moments to obtain the air concentration field correction value of the radioactive release source item. The invention also provides a correction device, a nuclear accident consequence evaluation method and a nuclear accident consequence evaluation system, which are used for solving the problem that the result calculated by the conventional nuclear accident consequence evaluation method has deviation with the actual physical process.

Description

Correction method and device for nuclear accident diffusion result and consequence evaluation method and system
Technical Field
The invention belongs to the field of nuclear facility accident radioactive consequence analysis, 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 a nuclear accident consequence.
Background
Under the condition that the nuclear facility has an accident, the spatial-temporal distribution of the radioactivity activity concentration and the radiation dose around the nuclear facility is analyzed, predicted and evaluated according to the radioactivity release source item and the meteorological conditions, and the method is an important component of the processes of emergency preparation and response of the nuclear accident, post-accident evaluation, environmental impact evaluation and the like. In the nuclear engineering design work, the safety analysis of nuclear facilities, the environmental impact evaluation, the emergency planning division, the establishment of a regional nuclear emergency scheme, the three-level PSA analysis and the like all need to take reliable accident consequence evaluation results as data bases.
Due to the complexity and the complexity of the radionuclide decay process, the processing of the decay process by the existing consequence analysis and evaluation models at home and abroad is simplified, so that the calculation result has deviation. Due to differences of disciplines, the atmospheric diffusion model adopted in the current international mainstream mostly does not consider or only considers the radionuclide decay process in a limited way, and in the actual situation, meteorological data dynamically changes along with time, diffusion tracks of radioactive substances released at different moments are different, the concentration at a certain position at a certain time is generated by overlapping a plurality of different release sections, and the decay time of the diffusion process is difficult to determine. At present, most mainstream atmospheric diffusion models are completely constructed, intermediate results cannot be given, secondary development is difficult, and limitation and difficulty exist when the atmospheric diffusion models are introduced and applied to nuclear accident consequence evaluation.
Disclosure of Invention
The invention provides a method and a device for correcting an atmospheric diffusion result of a radioactive substance in a nuclear accident, and an outcome evaluation method and an outcome evaluation system, which are used for solving the problem that the result calculated by the conventional nuclear accident outcome evaluation method is deviated from the 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 radioactive substance in a nuclear accident, including the following steps:
acquiring radioactive release source item data;
splitting the radioactive release source item data into M sub-release sections according to a first preset time length, wherein the first preset time length is a positive integral multiple of delta t, and the delta t is a minimum time interval of an atmospheric diffusion model calculation result;
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 second preset time duration, wherein the second preset time duration is positive integral multiple of delta t, N is integral multiple of M, and the diffusion results comprise air concentration fields;
calculating decay multiples of the nth moment respectively corresponding to the N groups of diffusion results of each sub-release section, wherein N is a positive integer and is more than or equal to 1 and less than or equal to N;
and performing a correction operation on the air concentration field at the nth moment of each sub-release segment: multiplying the 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 the air concentration field correction value at the nth moment of each sub-release section;
and performing summation operation on the air concentration field at the nth moment of the M sub-release sections: summing the air concentration field correction values of the M sub-release sections at the nth moment to obtain the air concentration field correction value of the radioactive release source item at the nth moment;
and respectively executing correction operation and summation operation on the air concentration fields at the N moments to obtain the air concentration field correction value of the radioactive release source item.
Preferably, the diffusion result further includes a deposition flux field, and after obtaining the air concentration field correction value of the radioactive emission source term, the method further includes:
performing a correction operation on the deposition flux field at the nth time of each sub-release segment: dividing the deposition flux field in the diffusion result of each sub-release section at the nth moment by the air concentration field at the moment, and multiplying the divided deposition flux field by the correction value of the air concentration field at the moment to obtain the correction value of the deposition flux field at the nth moment of each sub-release section;
performing a summation operation on the deposition flux fields at the nth time of the M sub-release segments: summing the deposition flux field correction values of the M sub-release sections at the nth moment to obtain the deposition flux field correction value of the radioactive release source item at the nth moment;
and respectively executing correction operation and summation operation on the deposition flux fields at the 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 of the radioactive emission source term, the method further comprises:
carrying out decay correction on the ground sediment concentration field at the nth-1 moment by taking the decay time as a second preset time length to obtain the contribution amount from the ground sediment concentration field at the nth-1 moment in the ground sediment concentration field at the nth moment, wherein N is a positive integer, and is more than or equal to 1 and less than or equal to N;
and multiplying the obtained deposition flux field correction value of the radioactive release source item at the nth moment by a second preset time length to obtain the contribution amount of air deposition to the ground in the [ n-1, n ] time period, and then summing the contribution amount of the ground deposition concentration field from the nth-1 moment to obtain the ground deposition concentration field at the nth moment.
Preferably, the calculating the decay multiple of the nth time respectively corresponding to the N groups of diffusion results of each sub-release segment specifically includes:
according to the decay chain information of all nuclides in the radioactive release source item data, taking each nuclide as a parent nucleus and respectively constructing a decay tree;
and according to the decay tree of each parent nucleus, depth-first search is adopted, the decay multiples of each parent nucleus and each child nucleus at the nth moment respectively corresponding to the N groups of diffusion results of each child release segment are calculated in parallel, and the decay multiples are stored in the decay tree.
Preferably, the treating each nuclide as a parent nucleus and respectively constructing a decay tree specifically includes:
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 a child nucleus is a child node of the decay tree;
arranging and storing all nodes in the decay tree as nodes in a single-direction linked list, wherein each node stores node information, and the node information comprises physical attributes such as the name of a nuclide, a decay constant, a decay branch ratio and the like, the position of the nuclide in the single-direction linked list and the position of each sub-core of the nuclide in the single-direction linked list;
and constructing the directed relation of all nodes according to the positions of the nuclides in the single linked list and the positions of all sub-kernels of the nuclides in the single linked list to form the decay tree.
In a second aspect, an embodiment of the present invention provides a method for evaluating a nuclear accident consequence, including the following steps:
calculating the effective dose of the public according to the air concentration field correction value and the ground deposition concentration field of the radioactive release source item obtained in the first aspect and dose conversion factors corresponding to the air concentration field of each nuclide, the ground deposition concentration field and the inhaled quantity of the human body from the air, wherein the effective dose comprises the air immersion external irradiation dose, the ground deposition external irradiation dose and the inhaled internal irradiation dose in the irradiated time period;
and comparing the effective dose with the universal optimization intervention level, and determining a protection action according to the comparison result.
Preferably, the comparing the effective dose with the universal optimal intervention level and determining the protective action according to the comparison result specifically comprises:
if T2d>10mSv, using stealth action; if T7dTg 50mSv, evacuation action is adopted, wherein T2d represents total effective dose for 2 days, and T7d represents 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 radioactive release source item data and splitting the radioactive release source item data into M sub-release sections according to a first preset time length, wherein the first preset time length is a positive integral multiple of delta t, and the delta t is a minimum time interval of an atmospheric diffusion model calculation result;
the acquisition module is connected with the splitting module and 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 second preset time duration, wherein the second preset time duration is positive integral multiple of delta t, N is integral multiple of M, and the diffusion results comprise air concentration fields;
the calculation module is connected with the acquisition module and is used for calculating the decay multiple of the nth moment respectively corresponding to the N groups of diffusion results of each sub-release section, wherein N is a positive integer and 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 segment: the air concentration field in the diffusion result at the nth moment of each sub-release section acquired by the acquisition module is multiplied by the decay multiple corresponding to the moment calculated by the calculation module to obtain the air concentration field correction value at the nth moment of each sub-release section; and the air concentration field at the nth moment of the M sub-release sections is used for performing summation operation: 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 is used for respectively performing correction operation and summation operation on the air concentration fields at the N moments to obtain the air concentration field correction value of the radioactive release source item.
Preferably, the correction device further includes a second correction module, the diffusion result further includes a deposition flux field, and the second correction module is connected to the acquisition module, the calculation module, and the first correction module, and is configured to perform a correction operation on the deposition flux field at the nth time of each sub-release segment: the device is used for dividing the deposition flux field in the diffusion result at the nth moment of each sub-release section acquired by the acquisition module by the air concentration field at the moment, and multiplying the deposition flux field by the air concentration field correction value at the nth moment acquired by the first correction module to acquire the deposition flux field correction value at the nth moment of each sub-release section; and for performing a summation operation on the deposition flux fields at the nth time instant of the M sub-release segments: the deposition flux field correction value at the nth moment of the M sub-release sections is summed 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 executed on the deposition flux fields at the N moments to obtain the deposition flux field correction value of the radioactive release source item.
Preferably, the correction device further comprises 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 sediment concentration field at the nth-1 moment by using the decay time as a second preset time length to obtain the contribution amount from the ground sediment concentration field at the nth-1 moment in the ground sediment concentration field at the nth moment, wherein N is a positive integer, and N is greater than or equal to 1 and less than or equal to N; and the device is used for multiplying the deposition flux field correction value of the radioactive release source item at the nth time obtained by the second correction module by a second preset time length, and then summing the deposition flux field correction value and the contribution amount of the ground deposition concentration field from the nth-1 time to obtain the ground deposition concentration field at the nth time.
Preferably, the calculation module includes a decay multiple unit, and the decay multiple unit is connected to the splitting module and the obtaining module, and is configured to obtain decay chain information according to all nuclides in the release source item obtained by the splitting module, and regard each nuclide as a parent nucleus according to the decay chain information and respectively construct a decay tree; and the system is used for calculating the decay multiple of each mother nucleus and each subnucleus at the nth moment respectively corresponding to the N groups of diffusion results of each sub-release section in parallel by adopting depth-first search according to the decay tree of each mother nucleus and storing the decay multiples into the decay tree.
In a fourth aspect, an embodiment of the present invention provides a nuclear accident consequence 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 the effective dose of the public according to dose conversion factors corresponding to the air concentration of each nuclide, the ground deposition concentration and the inhaled quantity of the human body from the air, 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 inhaled internal irradiation dose in the illuminated time period; the evaluation module is stored with a general optimization intervention level, is prestored with a mapping table of a comparison result and a protection action, is connected with the dosage module, and is used for comparing the general optimization intervention level with the effective dosage calculated by the dosage module and determining the corresponding protection action according to the comparison result and the mapping table.
According to the method for correcting the atmospheric diffusion result of the radioactive substance in the nuclear accident, the radioactive release source item is divided into M sub-release sections according to the first preset time length, wherein the first preset time length is positive integral multiple of delta t, and the delta t is the minimum time interval of the calculation result of the atmospheric diffusion model, so that the whole diffusion process is divided into zero, and the simulation of the real diffusion process that the diffusion concentration at a certain position at a certain time is generated by overlapping 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 second preset time duration, wherein the second preset time duration is positive integral multiple of delta t, and N is integral multiple of M. And calculating the decay multiple of the nth moment respectively corresponding to the N groups of diffusion results of each sub-release section, wherein N is a positive integer and is more than or equal to 1 and less than or equal to N. And performing a correction operation on the air concentration field at the nth moment of each sub-release segment: multiplying the obtained air concentration field in the diffusion result of the nth moment of each sub-release segment by the decay multiple corresponding to the moment to obtain the air concentration field correction value of the nth moment of each sub-release segment, and then performing summation operation on the air concentration fields of the nth moments of the M sub-release segments: and summing the air concentration field correction values of the M sub-release sections at the nth moment to obtain the air concentration field correction value of the radioactive release source item at the nth moment, and respectively performing correction operation and summation operation on the air concentration fields at the N moments to obtain the air concentration field correction value of the radioactive release source item. The deposition flux field correction is then calculated from the air concentration field correction. 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, corresponding correction is only performed on the existing calculation method, and the calculation process is convenient and rapid. 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 embodiment 1 is a flow chart of a correction method of atmospheric diffusion results of nuclear accident radioactive substances;
FIG. 2: the invention embodiment 2 is a one-way linked list structure diagram of a nuclear accident consequence evaluation method;
FIG. 3: the decay tree structure schematic diagram of the nuclear accident consequence evaluation method in embodiment 2 of the invention;
FIG. 4: the decay time schematic diagram of the diffusion result of the sub-release segment of the nuclear accident consequence evaluation method in embodiment 2 of the invention is shown;
FIG. 5: the invention embodiment 2 is a flow chart of a nuclear accident consequence evaluation method;
FIG. 6: the invention provides a device for correcting the atmospheric diffusion result of a nuclear accident radioactive substance in embodiment 3.
Detailed Description
In order to make the technical solutions of the present invention better understood, the present invention is further described in detail below with reference to the accompanying drawings and examples.
The decay process in the atmospheric diffusion of the radionuclide in reality is complex, and the specific expression is as follows: the meteorological data dynamically changes with time, the diffusion tracks of radioactive substances released at different times are different, and the concentration at a certain position at a certain time is generated by overlapping a plurality of different release sections, so that the decay time of the diffusion process is difficult to determine. In the method for correcting the atmospheric diffusion result of the radioactive substance in the nuclear accident provided by this embodiment, the diffusion process of the radioactive release source item is broken into whole parts, and the radioactive release source item is divided into a plurality of sub-release sections, and the diffusion result is corrected by using the time interval between the release end time and the diffusion result at each time obtained by independent diffusion thereof as decay time, where the diffusion result includes an air concentration field and a deposition flux field.
Furthermore, from the laws of atomic physics, it is known that the decay times of a parent nucleus and its child nucleus (for a non-existing child nucleus before decay, it may also be referred to as a growth factor, hereinafter collectively referred to as the decay times) are related only to the decay time. Further, the inventors of the present invention have found that the corrected concentration of a certain nuclide is only related to the pre-corrected concentration and the decay multiple, and the sum of the corrected concentrations of the nuclide is generated for all the parent nuclei, so that the correction algorithm provided by the embodiment of the present invention performs the correction calculation of the air concentration field by multiplying the pre-corrected concentration by the decay multiple at the corresponding time.
Example 1:
as shown in fig. 1, the present embodiment provides a method for correcting the result of atmospheric diffusion of a radioactive substance in a nuclear accident. The correction method can be applied to a nuclear accident consequence evaluation method, and can directly carry out complete decay correction and calculation on the air and deposition concentration field output by the atmospheric diffusion model. The correction method comprises the following steps:
step 101, radioactive release source item data is acquired. Specifically, the radioactive release source item data includes the radionuclide species in the source item, the release rate, the release duration, and the like.
Step 102, splitting the radioactive release source item data into M sub-release sections according to a first preset time length, wherein the first preset time length is a positive integral multiple of delta t, the delta t is a minimum time interval of an atmospheric diffusion model calculation result, and M is a positive integer.
Specifically, for a release source item with a total release duration D, a first preset duration D is taken, and the release source item is split into M sub-release segments, where M is D/D (that is, the value of M is the total release duration divided by the first preset duration), the mth sub-release segment has an actual release rate value only in the [ (M-1) D, md ] time period, no pollutant is released in other time periods, or the release rate in other time periods is equal to 0, M is a positive integer, and M is greater than or equal to 1 and less than or equal to M. It should be noted that, when the release source item contains multiple radionuclides (nuclides for short), the release of each nuclide needs to be split according to the method, and the numerical value of each sub-release segment after splitting is independent of the nuclide species. When the value of the first preset time length d is larger, the total number M of the sub-release sections is smaller, the calculation cost is lower, and the accuracy of decay correction is correspondingly reduced.
And 103, 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 second preset time duration, wherein the second preset time duration is positive integral multiple of delta t, N is integral multiple of M, and N is a positive integer.
In this embodiment, the diffusion results of the M sub-release sections obtained through calculation of the meteorological data and the atmospheric diffusion model may be obtained through two ways, one is to set the meteorological data and the atmospheric diffusion model first, and then obtain the diffusion results through calculation of the meteorological data and the atmospheric diffusion model, and the other is to directly obtain the diffusion results obtained through calculation of the external meteorological data and the atmospheric diffusion model. The atmospheric diffusion models in the two acquisition modes can be various pollutant diffusion models established and tested by the atmospheric physics profession, and do not include the decay process of radioactive pollutants, because the decay process of the radioactive pollutants is superposed with the correction method provided by the embodiment, the correction result is deviated. Wherein the diffusion results include an air concentration field and a deposition flux field. Each sub-release segment has N groups of diffusion results, where N is D/a second preset duration (i.e., the value of N is the total release duration D divided by the second preset duration). If the computer processor core is enough in calculation, and the values of the first preset time length d and the second preset time length are both delta t, when N is M, the decay correction error is minimum, namely the decay correction accuracy is highest. In addition, in order to improve the calculation efficiency, the diffusion results of the sub-release sections are calculated in parallel, the parallel total number is the total number M of the sub-release sections, and each sub-release section has N groups of diffusion results with the interval of a second preset time length.
And 104, calculating the decay multiple of the nth moment respectively corresponding to the N groups of diffusion results of each sub-release section, wherein N is a positive integer and is more than or equal to 1 and less than or equal to N.
In this embodiment, since the diffusion results output by the atmospheric diffusion model and spaced by the second preset time period are time-varying and recorded at the nth time, respectively, to obtain the correction value of the air concentration field in the diffusion results, the decay multiple corresponding to the recording time needs to be obtained, and therefore the decay multiple corresponding to the nth time respectively corresponding to the N diffusion results of each sub-release segment needs to be calculated. For example, when n is 1, the multiple of decay at the end of the 1 st group (i.e., the 1 st time) is calculated corresponding to the 1 st group diffusion result, and when n is 5, the multiple of decay at the end of the 5 th group (i.e., the 5 th time) is calculated corresponding to the 5 th group diffusion result.
Optionally, the method for calculating the decay multiple in this embodiment is as follows: and obtaining attenuation chain information of all nuclides according to the nuclide species in the release source item data, and regarding each nuclide as a parent nucleus according to the attenuation chain information and respectively constructing a decay tree. And then according to the decay tree of each parent nucleus, depth-first search is adopted, the decay multiples of each parent nucleus and each child nucleus at the nth moment respectively corresponding to the N groups of diffusion results of each child release segment are calculated in parallel, and the decay multiples are stored in the decay tree. Wherein, the mother nucleus is the radioactive nuclide contained in the source item, and the daughter nucleus is the nuclide newly generated by the decay of the mother nucleus. Of course, other known ways of calculating decay times may be used.
The decay tree data structure provided in the embodiment can completely store the multi-generation decay chain information of the radioactive nuclide in a form of the single-chain table without any decay chain truncation or simplification, solves the problems of numerous nuclides, complex and fussy decay chain structure, and provides support for decay multiple calculation. Meanwhile, the decay tree structure can simplify a complex decay chain structure, the decay chain can be quickly traversed through depth-first search, and the efficiency of multi-generation decay calculation can be improved. In addition, the decay multiples of the mother nucleus and the daughter nucleus at the nth moment respectively corresponding to the N groups of diffusion results of each sub-release segment are calculated in parallel, so that the calculation efficiency can be further improved.
Specifically, regarding each nuclide as a parent nucleus and respectively constructing a decay tree, the method specifically comprises the following steps:
step 1041, obtaining the multigeneration decay chain information of the mother nucleus, and storing each nuclide in the decay chain information as a node of the decay tree, wherein the mother nucleus is a root node of the decay tree, and the child nucleus is a child node of the decay tree.
Step 1042, arranging and storing all nodes in the decay tree as nodes in a single-direction linked list, wherein each node stores node information, and the node information comprises physical attributes such as the name of a nuclide, a decay constant, a decay branch ratio and the like, the position of the nuclide in the single-direction linked list, and the position of each sub-core of the nuclide in the single-direction linked list.
Step 1043, constructing the directed relationship of all nodes according to the positions of the nuclides in the single linked list and the positions of each sub-core of the nuclides in the single linked list to form a decay tree.
Step 105, performing correction operation on the air concentration field at the nth moment of each sub-release segment: and multiplying the obtained air concentration field at the nth moment of each sub-release section by the decay multiple corresponding to the moment to obtain the air concentration field correction value at the nth moment of each sub-release section. According to the physical law of atomic nuclei, the decay multiple of the parent nucleus and the child nucleus is only related to decay time, so that the corrected concentration of a certain nuclide is only related to the concentration before correction and the decay multiple, and the corrected value of the air concentration field at the nth moment of each child release segment is obtained on the basis of the product of the air concentration field at the moment and the corresponding decay multiple of the moment of each child release segment output by an atmospheric diffusion model.
In this embodiment, N groups of diffusion results of M sub-release segments 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 segment and the decay multiple corresponding to that time is calculated in parallel, where the total number of parallel calculations is M, and the calculation efficiency can be improved by reading the diffusion results in parallel and calculating the product of the air concentration field and the decay multiple in parallel.
Step 106, performing summation operation on the air concentration field at the nth moment of the M sub-release sections: and summing the air concentration field correction values of the M sub-release sections at the nth moment to obtain the air concentration field correction value of the radioactive release source item at the nth moment.
The radioactive release source item is divided into M sub-release sections, so that a complete decay process is supplemented, the concentration of parent nuclei is reduced, and the concentration of sub-nuclei is added, so that the corrected values of the air concentration field at the nth moment of the M sub-release sections are summed to obtain the corrected value of the air concentration field at the nth moment of the release source item. Finally, the air concentration field correction value of the radioactive release source item is obtained by executing N times to the step 105 and the step 106, namely respectively executing correction operation and summation operation to the air concentration field at N times.
Optionally, after obtaining the air concentration field correction value of the radioactive release source item, the method for correcting the atmospheric diffusion result of the radioactive substance in the nuclear accident further includes calculating the deposition flux field correction value, and 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: and dividing the deposition flux field in the diffusion result of the nth moment of each sub-release section by the air concentration field at the moment, and multiplying the division by the air concentration field correction value at the moment to obtain the deposition flux field correction value of the nth moment of each sub-release section. The summation operation is then performed on the deposition flux fields at the nth moment of the M sub-release segments: and summing the deposition flux field correction values of the M sub-release sections at the nth moment to obtain the deposition flux field correction value of the radioactive release source item at the nth moment. By executing the calculation process for N times, namely respectively executing correction operation and summation operation on the deposition flux fields at N moments, N groups of deposition flux field correction values of the radioactive release source term can be obtained.
Optionally, after obtaining the deposition flux field correction value of the radioactive emission source item, the method for correcting the atmospheric diffusion result of the radioactive substance in the nuclear accident further includes calculating a ground deposition concentration field:
carrying out decay correction on the ground sediment concentration field at the nth-1 moment by taking the decay time as a second preset time length to obtain the contribution amount from the ground sediment concentration field at the nth-1 moment in the ground sediment concentration field at the nth moment, wherein N is a positive integer, and is more than or equal to 1 and less than or equal to N; and multiplying the obtained deposition flux field correction value of the radioactive release source item at the nth moment by a second preset time length to obtain the contribution amount of air deposition to the ground in the [ n-1, n ] time period, and summing the contribution amount of the air deposition to the ground from the ground deposition concentration field at the nth moment in the ground deposition concentration field to obtain the ground deposition concentration field at the nth moment.
In the present embodiment, when n is 1, that is, the ground deposition concentration field at the 0 th time is 0, because the 0 th time has not been diffused, the ground deposition concentration field at the 1 st time is the product of the deposition flux field correction value at the 1 st time of the radioactive release source term and the second preset time period.
Example 2:
the embodiment provides a nuclear accident consequence evaluation method, which comprises the following steps:
step 201, calculating an effective dose 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 inhaled quantity of the human body from the air obtained in the embodiment 1, wherein the effective dose includes an air immersion external irradiation dose, a ground deposition external irradiation dose, and an inhaled internal irradiation dose in an irradiated time period. The technology of air immersion external irradiation dose considers the dose irradiation of a certain point on the ground in a three-dimensional space, and the ground deposition external irradiation dose does not consider inert gas.
Step 202, comparing the effective dose with the universal optimization intervention level, and determining a protective action according to the comparison result. The corrected diffusion result is more accurate, so that the result evaluation result calculated by the correction method in the embodiment 1 is more accurate, and a safe protection suggestion is provided for the public.
In the embodiment, the irradiation time period is 2 days, 7 days, 1 month, 1 year and the like corresponding to the universal optimization intervention level in national standard GB18871-2002, and the calculated irradiation dose is the accumulated dose in the irradiation time period and comprises an effective dose and an organ dose.
In particular, if T2d>10mSv, concealed action should be taken; if T7dMore than 50mSv, an evacuation action should be taken, where T2dRepresents the total effective dose, T, for 2 days7dRepresents the total effective dose for 7 days.
The beneficial effects of the above embodiment are as follows:
(1) the radioactive release source item is divided into M sub-release sections according to a first preset time length, wherein the first preset time length is a positive integral multiple of delta t, and the delta t is the minimum time interval of the calculation result of the atmospheric diffusion model. So as to realize the breaking of the diffusion process into whole parts, and facilitate the simulation of the real diffusion process that the diffusion concentration at a certain position at a time is generated by the superposition of a plurality of different release sections.
(2) Multiplying the obtained air concentration field of each sub-release section at the nth moment by the decay multiple corresponding to the nth moment to obtain the air concentration field correction value of each sub-release section at the nth moment, summing the air concentration field correction values of the M sub-release sections at the nth moment to obtain the air concentration field correction value of the radioactive release source item at the nth moment, and calculating the correction value of the deposition flux field according to the air concentration field correction value. And executing the calculation for N times to finish the 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 atmospheric diffusion result of the nuclear accident radioactive substance can directly carry out complete decay correction on the air concentration field and the deposition flux field output by the atmospheric diffusion model, and does not need to carry out secondary development on the existing atmospheric diffusion model, so that the correction method has strong applicability.
(4) Each nuclide releasing a source item is taken as a parent nucleus, a decay tree is respectively constructed for each parent nucleus, and a depth-first search traversal decay chain is adopted, so that the efficiency of multi-generation decay calculation can be improved. And the calculation efficiency can be improved by parallelly calculating the diffusion result, parallelly calculating the decay multiple, parallelly reading the diffusion result of the sub-release section and parallelly correcting the diffusion result of the sub-release section.
(5) Because the corrected diffusion result is more accurate, the effective dose is calculated through the air concentration field correction value and the ground deposition concentration field, and the public safety can be better ensured according to the effective protection action obtained by the effective dose.
As shown in fig. 2 to 5, the method for evaluating the nuclear accident consequence in this embodiment is implemented by a specific example, in this example, the release source item only contains one nuclide I-135, and the release duration is 2 hours 7200 seconds, and the method includes the following steps:
step 301, setting a first preset time duration d to 1200 seconds, and splitting the release source item into M sub-release segments with d as a time interval, that is, M is 2h/1200s is 6.
Step 302, calling meteorological data and an atmospheric diffusion model, and calculating N groups of diffusion results of each sub-release segment within 2 hours, namely calculating an I-135 air concentration field and a deposition flux field of each grid point of each time interval within the research range, wherein the time interval delta t of the calculation results of the atmospheric diffusion model is 600 seconds. And setting the second preset time period equal to the time t, so that each sub-release segment has N groups of diffusion results with the time interval of the second preset time period, wherein N is 2 h/delta t, 7200s/600s is 12.
Step 303, calculating decay times of the mother nucleus and the daughter nucleus.
Firstly, taking a source item nuclide I-135 as a parent nucleus, selecting previous three generations of decay chain information, and constructing a decay tree of the parent nucleus I-135 according to the decay chain information.
The third generation decay chain information of the parent nucleus I-135 comprises sub nuclei Xe-135m, Cs-135 and Xe-135. As shown in FIG. 2, a schematic diagram of a singly linked list structure for storing decay trees is shown, in which each nuclear species is stored as a node, where the number before the name of the nuclear species is the position of the nuclear species node in the singly linked list, such as the number "0" before I-135, which means that the position of I-135 in the singly linked list is the start bit, and in each node, an indefinite length array of the positions of its sub-cores in the singly linked list is also stored after the name of the nuclear species for fast locating and searching, for example, the numbers "1" and "3" in the array of the 0 th node I-135 respectively indicate that the sub-core of the mother core I-135 includes Xe-135m of the 1 st node and Xe-135 of the 3 rd node.
The decay tree structure completely stored by adopting the unidirectional chain table structure is shown in fig. 3, wherein the 0 th node of the unidirectional chain table represents a mother core, the mother core is a root node of the decay tree, the other nodes of the unidirectional chain table represent sub-cores, the sub-cores are sub-nodes of the decay tree, and the directed relation of all the nodes is constructed according to the positions of nuclide nodes in the unidirectional chain table and the positions of the respective sub-cores of the nuclide in the unidirectional chain table to form the decay tree, so that the decay chain is conveniently traversed by using a depth-first search method, and the calculation efficiency of decay multiples is improved.
Secondly, since the first preset time period d is 1200 seconds and the second preset time period is 600 seconds, the decay time of each diffusion result of each sub-release segment is as shown in fig. 4. The decay time is an integral multiple of Δ t, N is 1 to N-2, and since the m-1-th sub-release segment is released within the initial 2 Δ t time and overlaps with the time segment 2 Δ t of the 1 st time (N-1) and the 2 nd time (N-2), the same process is applied to the other sub-release segments in the present embodiment regardless of the decay of the sub-release segments during the discharge process.
Along the decay tree of parent core I-135, the decay time is calculated as a decay multiple of N Δ t, where N ═ 1, 2. Specifically, from the half-life and the decay branch ratio, the concentration of I-135 itself and its daughter nuclei varied in multiples (decay times) after 600s, 1200s, 1800s, … …, and 6000s decay. Depth-first search is adopted, the decay tree shown in fig. 3 is searched downwards along the left side, after the leftmost branch is searched, the right search is performed step by step, and the specific sequence is as follows:
(1) checking the parent nucleus I-135, and respectively calculating decay multiples after the decay time required to be calculated according to the half-life period of the parent nucleus I-135;
(2) checking the first left subnucleus Xe-135m, calculating the decay times after the decay time required to be calculated respectively according to the half-life of the mother nucleus I-135, the self half-life and the branch ratio of I-135 decaying into Xe-135 m;
(3) finding the tail-most nuclide Cs-135 at the left side, and calculating the decay multiple according to the half-life of the parent nucleus I-135, parent nucleus Xe-135m and the parent nucleus itself and the decay branch ratio between the parent nucleus and the parent nucleus; since Cs-135 decays without producing a new daughter nucleus, the leftmost branch of the decay tree structure reaches the end;
(4) returning to Xe-135m, finding another sub-nucleus Xe-135 to the right side, and calculating the decay times of the sub-nucleus Xe-135 along the I-135- > Xe-135m- > Xe-135;
(5) finding the nucleus Cs-135 of Xe-135 to reach the end of the second branch of the decay tree structure; calculating the decay multiple of the Cs-135 on the branch along the I-135- > Xe-135m- > Xe-135- > Cs-135, and overlapping the decay multiple calculated in the step (3) to update the decay multiple of the Cs-135;
(6) after returning to Xe-135, finding that only the branch found in step (5) returns to Xe-135m upwards, finding that both branches have been found in steps (3) and (5), then continuing to return to root I-135 upwards and finding another sub-nucleus Xe-135 to the right; calculating the decay multiple of the Xe-135 on the branch along the I-135- > Xe-135, and adding the decay multiple calculated in the step (4) to update the decay multiple of the Xe-135;
(7) finding the sub-nucleus Cs-135 of Xe-135 to reach the end of the third branch of the decay tree structure; calculating the decay multiple 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 time, all branches of the decay tree structure are traversed.
For example, for a decay time of 600s, the decay times for I-135, Xe-135m, Cs-135 and Xe-135 are: 0.98, 0.01, 3.18E-14 and 0.06; for a decay time 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 step 304, carrying out parallel correction, and calculating 12 groups of diffusion result correction values and ground sediment concentration fields of the release source items by circulating for 12 times of N-times.
Firstly, reading the diffusion result of each sub-release section output by the atmospheric diffusion model M times in parallel, correcting the air concentration field in the nth group of results of each sub-release section M times in parallel according to the decay time shown in FIG. 4 to obtain the air concentration field correction value of the nth time of each sub-release section (wherein the nth time is the last time of the nth group of diffusion results), and then overlapping the air concentration field correction values of the nth time of all M sub-release sections to obtain the air concentration field correction value of the nth time of the release source item. And then carrying out decay correction on the ground sediment concentration field with the decay time delta t, and accumulating the sediment flux field correction value at the nth moment to the ground sediment concentration field.
For example, as shown in fig. 4, for n-1, the release of the m-1 th sub-release segment is not completed, and the decay correction is not performed on the diffusion result under conservative consideration, and the m-2 th to m-6 th sub-release segments have not been released with radioactivity, and all the diffusion results are 0. Since the air concentration field correction value at the time point n-1 is the cumulative sum of the air concentration field correction values of the sub release segments m-1 to m-6, it is numerically equal to the air concentration field of the sub release segment m-1. At this time, the deposition flux field in the diffusion result at the nth-1 moment in the obtained m-1 sub release segment is divided by the air concentration field at the moment, and then multiplied by the air concentration field correction value at the nth-1 moment to obtain the deposition flux field correction value at the nth-1 moment in the m-1 sub release segment, 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 nth-1 moment.
As can be seen from fig. 4, when N is 12, the diffusion results of all the sub-release segments except the M-th sub-release segment, i.e., M is 6, are corrected for decay, and the decay times are Δ t (6000 s, 4800s, 3600s, 2400s, and 1200s) 10 times, 8 times, 6 times, 4 times, and 2 times, respectively. The corrected value of the air concentration field at the n-th time point is as follows: and multiplying the air concentration fields of the m-1 th to m-5 th sub-release sections by corresponding decay multiples respectively, and adding the multiplied air concentration fields to the air concentration fields of the m-6 th sub-release sections to obtain the corrected value of the air concentration field of the Nth group, namely the 12 th moment. And (3) carrying out decay correction with the decay time of delta t on the ground deposition concentration field with the n-11 th moment, and multiplying the deposition flux field correction value with the delta t at the n-12 th moment by the delta t to accumulate the deposition flux field correction value into the ground deposition concentration field to obtain the ground deposition concentration field at the 12 th moment. And the deposition flux field correction value at the nth-12 moment is obtained by dividing the deposition flux field at the moment in each sub-releasing section by the air concentration field at the moment, multiplying the division by the air concentration field correction value at the moment to obtain the deposition flux field correction value at the nth-12 moment in each sub-releasing section, and then summing the deposition flux field correction values at the nth-12 moments in the M sub-releasing sections.
For n-2, 3 …, the calculation process is similar to that described above and will not be described herein.
And 305, multiplying the air concentration field correction values of I-135, Xe-135m, Cs-135 and Xe-135 and the ground deposition concentration field by respective dose Conversion factors DCF (dose Conversion factor) respectively and overlapping to obtain an effective dose field. For example, assume I at a grid pointA ground layer air concentration of-135, Xe-135m, Cs-135 and Xe-135 of cair,1、cair,2、cair,3And cair,4The ground deposition concentration is csur,1、csur,2、csur,3And csur,4The effective dose DCF of the air immersion external irradiation is dair,1、dair,2、dair,3And dair,4The effective doses DCF of the ground deposition external irradiation are respectively dsur,1、dsur,2、dsur,3And dsur,4The effective dose DCF for inhalation is dinh,1、dinh,2、dinh,3And dinh,4And if the respiratory rate is h, the calculation formula of the effective dose of the public inhaled radiation at the grid point is shown as the formula (1): (c)air,1dair,1+cair,2dair,2+cair, 3dair,3+cair,4dair,4) Δ t. (1) The calculation formula of the effective dose of the inhaled radiation on the grid points in the public is shown as the formula (2): (c)sur,1dsur,1+csur,2dsur,2+csur,3dsur,3+csur,4dsur,4) Δ t. (2) The calculation formula of the effective dose of the inhaled radiation on the grid points in the public is shown as the formula (3): (c)air,1dinh,1+cair,2dinh,2+cair,3dinh,3+cair,4dinh,4)*h*Δt。(3)
If the 2-day effective dose of each of the three routes is A, B, C, the 2-day total effective dose T2d＝A+B+C。
And step 306, comparing the calculated effective dose with the universal optimization intervention level, and giving a protective action suggestion.
According to the national standard, if T2d>10mSv, the person at this grid point should be advised to take covert action, where T2dRepresents a total effective dose for 2 days; if T7d>50mSv, at which point persons at the grid point should be advised to take evacuation action, where T7dRepresents 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 material, 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 obtain radioactive release source item data, and split the radioactive release source item data into M sub-release sections according to a first preset time duration, where the first preset time duration is a positive integer multiple of Δ t, and Δ t is a minimum time interval of an atmospheric diffusion model calculation result.
An obtaining module 42, connected to the splitting module 41, configured to obtain diffusion results of M sub-release segments calculated through meteorological data and an atmospheric diffusion model, and enable each sub-release segment to have N groups of diffusion results with a time interval being a second preset time duration, where the second preset time duration is a positive integer multiple of Δ t, and N is an integer multiple of M; the diffusion result includes an air concentration field.
And the calculating module 43 is connected to the obtaining module 42, and is configured to calculate decay multiples of the nth time respectively corresponding to the 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, which 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 released source item obtained by the splitting module 41, and regard each nuclide as a parent nucleus according to the decay chain information and respectively construct a decay tree; and the system is used for calculating the decay multiple of each mother nucleus and each subnucleus at the nth moment respectively corresponding to the N groups of diffusion results of each sub-release section in parallel by adopting depth-first search according to the decay tree of each mother nucleus and storing the decay multiple into the decay tree.
A first correcting module 44, connected to the obtaining module 42 and the calculating module 43, for performing a correcting operation on the air concentration field at the nth time of each sub-release segment: the air concentration field calculation module is used for multiplying the air concentration field in the diffusion result at 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 calculation module 43 to obtain the air concentration field correction value at the nth moment of each sub-release section; and the air concentration field at the nth moment of the M sub-release sections is used for performing summation operation: 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 is used for respectively performing correction operation and summation operation on the air concentration fields at the N moments to obtain the air concentration field correction value of the radioactive release source item.
Optionally, the device for correcting the atmospheric diffusion result of the nuclear accident radioactive substance further includes a second correction module, the diffusion result further includes 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 deposition flux field correction module is used for dividing the deposition flux field in the diffusion result at the nth moment of each sub-release section obtained by the obtaining module 42 by the air concentration field at the moment, and multiplying the division result by the air concentration field correction value at the nth moment obtained by the first correction module 44 to obtain the deposition flux field correction value at the nth moment of each sub-release section; and for performing a summation operation on the deposition flux fields at the nth time instant of the M sub-release segments: the deposition flux field correction values at the nth moment of the M sub-release sections are summed 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 executed on the deposition flux fields at the N moments 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 includes a third correction module, the third correction module is connected to the obtaining module 42 and the second correction module, and is configured to perform decay correction on the ground deposition concentration field at the nth-1 moment, where the decay time is a second preset time duration, so as to obtain a contribution amount from the ground deposition concentration field at the nth-1 moment in the ground deposition concentration field at the nth moment, where N is a positive integer, and N is greater than or equal to 1 and less than or equal to N; and the device is used for multiplying the deposition flux field correction value of the radioactive release source item at the nth moment obtained by the second correction module by the second preset time length of the acquisition module 42, and then summing the deposition flux field correction value and the contribution amount from the ground deposition concentration field at the nth moment in the ground deposition concentration field at the n-1 moment to obtain the ground deposition concentration field at the nth moment.
Example 4:
the embodiment provides a nuclear accident consequence evaluation system, which comprises: the apparatus for correcting the atmospheric diffusion result of the nuclear accident radioactive material according to embodiment 3 further includes a dosing module and an evaluation module.
And the dose module is connected with the correction device of the atmospheric diffusion result of the nuclear accident radioactive substance and is used for calculating the effective dose of the public according to the dose conversion factors corresponding to the air concentration of each nuclide, the ground deposition concentration and the inhaled quantity of the human body from the air, the air concentration field correction value of the radioactive release source item and the ground deposition concentration field, wherein the air concentration field correction value and the ground deposition concentration field are obtained by the correction device, and the effective dose comprises the air immersion external irradiation dose, the ground deposition external irradiation dose and the inhaled internal irradiation dose in the illuminated time period.
And the evaluation module is internally stored with a general optimized intervention level, is prestored with a mapping table of a comparison result and a protection action, is connected with the dosage module, is used for comparing the general optimized intervention level with the effective dosage calculated by the dosage module, and determines the corresponding protection action according to the comparison result and the mapping table.
It will be understood that the above embodiments are merely exemplary embodiments taken to illustrate the principles of the present invention, which is not limited thereto. It will be apparent to those skilled in the art that various modifications and improvements can be made without departing from the spirit and substance of the invention, and these modifications and improvements are also considered to be within the scope of the invention.

Claims (12)

1. A method for correcting an atmospheric diffusion result of a radioactive substance in a nuclear accident is characterized by comprising the following steps:
acquiring radioactive release source item data;
splitting the radioactive release source item data into M sub-release sections according to a first preset time length, wherein the first preset time length is a positive integral multiple of delta t, and the delta t is a minimum time interval of an atmospheric diffusion model calculation result;
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 second preset time duration, wherein the second preset time duration is positive integral multiple of delta t, N is integral multiple of M, and the diffusion results comprise air concentration fields;
calculating decay multiples of the nth moment respectively corresponding to the N groups of diffusion results of each sub-release section, wherein N is a positive integer and is more than or equal to 1 and less than or equal to N;
and performing a correction operation on the air concentration field at the nth moment of each sub-release segment: multiplying the 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 the air concentration field correction value at the nth moment of each sub-release section;
and performing summation operation on the air concentration field at the nth moment of the M sub-release sections: summing the air concentration field correction values of the M sub-release sections at the nth moment to obtain the air concentration field correction value of the radioactive release source item at the nth moment;
and respectively executing correction operation and summation operation on the air concentration fields at the N moments to obtain the air concentration field correction value of the radioactive release source item.
2. The method for correcting atmospheric diffusion result of radioactive substances in nuclear accident according to claim 1, wherein the diffusion result further comprises a sedimentary flux field,
after obtaining the air concentration field correction value of the radioactive emission source item, the method further comprises:
performing a correction operation on the deposition flux field at the nth time of each sub-release segment: dividing the deposition flux field in the diffusion result of each sub-release section at the nth moment by the air concentration field at the moment, and multiplying the divided deposition flux field by the correction value of the air concentration field at the moment to obtain the correction value of the deposition flux field at the nth moment of each sub-release section;
performing a summation operation on the deposition flux fields at the nth time of the M sub-release segments: summing the deposition flux field correction values of the M sub-release sections at the nth moment to obtain the deposition flux field correction value of the radioactive release source item at the nth moment;
and respectively executing correction operation and summation operation on the deposition flux fields at the N moments to obtain the deposition flux field correction value of the radioactive release source item.
3. The method for correcting atmospheric diffusion results of radioactive substances in nuclear accidents according to claim 2, wherein after obtaining the deposition flux field correction value of the radioactive release source item, the method further comprises:
carrying out decay correction on the ground sediment concentration field at the nth-1 moment by taking the decay time as a second preset time length to obtain the contribution amount from the ground sediment concentration field at the nth-1 moment in the ground sediment concentration field at the nth moment, wherein N is a positive integer, and is more than or equal to 1 and less than or equal to N;
and multiplying the obtained deposition flux field correction value of the radioactive release source item at the nth moment by a second preset time length to obtain the contribution amount of air deposition to the ground in the [ n-1, n ] time period, and then summing the contribution amount of the ground deposition concentration field from the nth-1 moment to obtain the ground deposition concentration field at the nth moment.
4. The method for correcting atmospheric diffusion results of radioactive substances in a nuclear accident according to claim 1, wherein the calculating the decay multiple of the nth time corresponding to each of the N groups of diffusion results of each sub-release segment specifically includes:
according to the decay chain information of all nuclides in the radioactive release source item data, taking each nuclide as a parent nucleus and respectively constructing a decay tree;
and according to the decay tree of each parent nucleus, depth-first search is adopted, the decay multiples of each parent nucleus and each child nucleus at the nth moment respectively corresponding to the N groups of diffusion results of each child release segment are calculated in parallel, and the decay multiples are stored in the decay tree.
5. The method for correcting the atmospheric diffusion result of the radioactive substance in the nuclear accident according to claim 4, wherein the step of constructing the decay tree by taking each nuclide as a parent nucleus comprises:
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 a child nucleus is a child node of the decay tree;
arranging and storing all nodes in the decay tree as nodes in a single-direction linked list, wherein each node stores node information, and the node information comprises the name of a nuclide, a decay constant and a decay branch ratio, the position of the nuclide in the single-direction linked list and the position of each sub-core of the nuclide in the single-direction linked list;
and constructing the directed relation of all nodes according to the positions of the nuclides in the single linked list and the positions of all sub-kernels of the nuclides in the single linked list to form the decay tree.
6. A nuclear accident consequence evaluation method is characterized by comprising the following steps:
calculating the effective dose of the public according to the air concentration field correction value of the radioactive release source term, the ground deposition concentration field and the dose conversion factors corresponding to the air concentration of each nuclide, the ground deposition concentration and the inhaled quantity of the human body from the air, wherein the effective dose comprises the air immersion external irradiation dose, the ground deposition external irradiation dose and the inhaled internal irradiation dose in the irradiated time period;
and comparing the effective dose with the universal optimization intervention level, and determining a protection action according to the comparison result.
7. The method for evaluating the outcome of a nuclear accident according to claim 6, wherein the comparing the effective dose with the universal optimal intervention level and the determining the protective action according to the comparison result specifically comprise:
if T2d>10mSv, using stealth action; if T7dTg 50mSv, and an evacuation action, wherein T2dRepresents the total effective dose, T, for 2 days7dRepresents the total effective dose for 7 days.
8. A correction device for atmospheric diffusion results of radioactive substances in nuclear accidents 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 radioactive release source item data and splitting the radioactive release source item data into M sub-release sections according to a first preset time length, wherein the first preset time length is a positive integral multiple of delta t, and the delta t is a minimum time interval of an atmospheric diffusion model calculation result;
the acquisition module is connected with the splitting module and 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 second preset time duration, wherein the second preset time duration is positive integral multiple of delta t, N is integral multiple of M, and the diffusion results comprise air concentration fields;
the calculation module is connected with the acquisition module and is used for calculating the decay multiple of the nth moment respectively corresponding to the N groups of diffusion results of each sub-release section, wherein N is a positive integer and 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 segment: the air concentration field in the diffusion result at the nth moment of each sub-release section acquired by the acquisition module is multiplied by the decay multiple corresponding to the moment calculated by the calculation module to obtain the air concentration field correction value at the nth moment of each sub-release section;
and the air concentration field at the nth moment of the M sub-release sections is used for performing summation operation: 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 is used for respectively performing correction operation and summation operation on the air concentration fields at the N moments to obtain the air concentration field correction value of the radioactive release source item.
9. The device for correcting the atmospheric diffusion result of the radioactive substances in the nuclear accident according to claim 8, further comprising a second correction module, wherein the diffusion result further comprises 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 executing correction operation on the deposition flux field at the nth moment of each sub-release segment: the device is used for dividing the deposition flux field in the diffusion result at the nth moment of each sub-release section acquired by the acquisition module by the air concentration field at the moment, and multiplying the deposition flux field by the air concentration field correction value at the nth moment acquired by the first correction module to acquire the deposition flux field correction value at the nth moment of each sub-release section;
and for performing a summation operation on the deposition flux fields at the nth time instant of the M sub-release segments: the deposition flux field correction value at the nth moment of the M sub-release sections is summed 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 executed on the deposition flux fields at the N moments to obtain the deposition flux field correction value of the radioactive release source item.
10. The apparatus for correcting atmospheric diffusion results of radioactive substances in nuclear accidents according to claim 9, further comprising a third correcting 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 sediment concentration field at the nth-1 moment by taking the decay time as a second preset time length to obtain the contribution amount from the ground sediment concentration field at the nth-1 moment in the ground sediment concentration field at the nth moment, wherein N is a positive integer, and N is more than or equal to 1 and less than or equal to N;
and the device is used for multiplying the deposition flux field correction value of the radioactive release source item at the nth time obtained by the second correction module by a second preset time length, and then summing the deposition flux field correction value and the contribution amount of the ground deposition concentration field from the nth-1 time to obtain the ground deposition concentration field at the nth time.
11. The apparatus for correcting atmospheric diffusion result of nuclear accident radioactive material according to claim 10, wherein the calculation module comprises a decay multiple unit,
the decay multiple unit is connected with the splitting module and the obtaining module and used for obtaining decay chain information according to all nuclides in the release source items obtained by the splitting module, and taking each nuclide as a parent nucleus according to the decay chain information and respectively constructing a decay tree;
and the system is used for calculating the decay multiple of each mother nucleus and each subnucleus at the nth moment respectively corresponding to the N groups of diffusion results of each sub-release section in parallel by adopting depth-first search according to the decay tree of each mother nucleus and storing the decay multiples into the decay tree.
12. A nuclear accident outcome evaluation system, comprising: the apparatus for correcting atmospheric diffusion results of radioactive substances in a nuclear accident according to claim 10, further comprising a dosing module and an evaluation module:
the dose module is connected with the correction device of the atmospheric diffusion result of the nuclear accident radioactive substance and is used for calculating the effective dose of the public according to the dose conversion factors corresponding to the air concentration of each nuclide, the ground deposition concentration and the inhaled quantity of the human body from the air, the air concentration field correction value of the radioactive release source item and the ground deposition concentration field, wherein the air concentration field correction value and the ground deposition concentration field are obtained by the correction device, and the effective dose comprises the air immersion external irradiation dose, the ground deposition external irradiation dose and the inhaled internal irradiation dose in the illuminated time period;
and the evaluation module is internally stored with a general optimized intervention level, is prestored with a mapping table of a comparison result and a protection action, is connected with the dosage module, is used for comparing the general optimized intervention level with the effective dosage calculated by the dosage module, and determines the corresponding protection action according to the comparison result and the mapping table.
CN202010355824.6A 2020-04-29 2020-04-29 Correction method and device for nuclear accident diffusion result and result evaluation method and system Active CN111695762B (en)

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 true CN111695762A (en) 2020-09-22
CN111695762B 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 (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20080040352A (en) * 2006-11-03 2008-05-08 한국전력공사 Risk-informed population hazard assessment
US20080208902A1 (en) * 2007-02-26 2008-08-28 Robert R Friedlander System and method for deriving a hierarchical event based database optimized for analysis of chaotic events
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
US20100324871A1 (en) * 2007-10-22 2010-12-23 Commissariat A L'energie Atomique Et Aux Energies Alternatives Method for determining a time course of an accident occurring in a risk-prone installation
US20120109618A1 (en) * 2010-11-02 2012-05-03 Institute Of Nuclear Energy Research, Atomic Energy Council, Executive Yuan Accident parameter identification method for severe accidents
US20120185222A1 (en) * 2010-11-23 2012-07-19 Westinghouse Electric Company, Llc Full spectrum loca evaluation model and analysis methodology
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

Patent Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20080040352A (en) * 2006-11-03 2008-05-08 한국전력공사 Risk-informed population hazard assessment
US20080208902A1 (en) * 2007-02-26 2008-08-28 Robert R Friedlander System and method for deriving a hierarchical event based database optimized for analysis of chaotic events
US20100324871A1 (en) * 2007-10-22 2010-12-23 Commissariat A L'energie Atomique Et Aux Energies Alternatives Method for determining a time course of an accident occurring in a risk-prone installation
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
US20120109618A1 (en) * 2010-11-02 2012-05-03 Institute Of Nuclear Energy Research, Atomic Energy Council, Executive Yuan Accident parameter identification method for severe accidents
US20120185222A1 (en) * 2010-11-23 2012-07-19 Westinghouse Electric Company, Llc Full spectrum loca evaluation model and analysis methodology
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)

* Cited by examiner, † Cited by third party
Title
侯姗姗;: "GIS及大气传输模型在核应急管理中的研究与系统实现" *
刘蕴;方晟;李红;曲静原;姚仁太;范丹;: "基于四维变分资料同化的核事故源项反演" *
胡啸峰;陈鹏;曾昭龙;: "基于WRF的放射性物质大气扩散模拟的环境辐射剂量估算方法" *
郭瑞萍;张琼;陈海英;潘昕怿;杨春林;: "核电厂放射性气态流出物大气弥散模型现状与展望" *

Also Published As

Publication number Publication date
CN111695762B (en) 2023-05-05

Similar Documents

Publication Publication Date Title
Forrest et al. FISPACT-2007: User manual
CN107145698A (en) Core and radiation accident consequence large scale analogy method and simulation system
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
Fiorito et al. Inventory calculation and nuclear data uncertainty propagation on light water reactor fuel using ALEPH-2 and SCALE 6.2
Abdurashitov et al. The BNO–LNGS joint measurement of the solar neutrino capture rate in 71Ga
Smith et al. An analysis of various solution strategies and perturbations on inputs of the reactor shielding problem
CN116028010A (en) Radionuclide hazard sequencing method
Mickael A fast, automated, semideterministic weight windows generator for MCNP
Jang et al. Verification and Validation of Back-end Cycle Source Term Calculation of the Nodal Code RAST-K
Perevolotskaya et al. Retrospective assessment of the formation of the radiation situation in pine plantations in the first year after the Chernobyl accident
Zherdev et al. ROCOCO: A constants supply system for Monte-Carlo reactor calculation
Wu et al. Three-dimensional dose evaluation system using real-time wind field information for nuclear accidents in Taiwan
Gougar et al. Matrix formulation of pebble circulation in the PEBBED code
Bahran et al. Development of a research reactor protocol for neutron multiplication measurements
Wilson et al. NESHAPs CY19 Report: Radionuclide Air Emission Report for 2019
Dyachkov et al. Methods for Measuring Daughter Products of Radon Decay in the Surface Atmospheric Layer of the Earth
Hofman Assimilation scenario for long-term deposition of 137Cs
Bourret et al. A Perspective on the Successes of the NNSS Underground Test Area (UGTA) Activity-20221
Holcomb et al. 10.1. SCALE Cross Section Libraries
Burn Learning aspects of the Direct Statistical Approach to the optimization of Monte Carlo radiation transport calculations
Omar TRIMON: THEORY AND GUIDELINES
López et al. WWER fuel rod isotopics by MONTEBURNS 1.0—Influence on the multiplication factor and comparison with the CB3 benchmark data
Kelley et al. TUNL contributions in the US nuclear data program and NSDD
Hooper et al. Predictive Fallout Composition Modeling: Improvements and Applications of the Defense Land Fallout Interpretive Code

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