CN117972266A - Calculation method and device for atmospheric dispersion ground concentration and dose evaluation method - Google Patents
Calculation method and device for atmospheric dispersion ground concentration and dose evaluation method Download PDFInfo
- Publication number
- CN117972266A CN117972266A CN202410008427.XA CN202410008427A CN117972266A CN 117972266 A CN117972266 A CN 117972266A CN 202410008427 A CN202410008427 A CN 202410008427A CN 117972266 A CN117972266 A CN 117972266A
- Authority
- CN
- China
- Prior art keywords
- evaluation
- emission source
- emission
- distance
- subarea
- 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.)
- Pending
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 212
- 239000006185 dispersion Substances 0.000 title claims abstract description 53
- 238000004364 calculation method Methods 0.000 title claims abstract description 19
- 238000000034 method Methods 0.000 claims abstract description 27
- 230000005855 radiation Effects 0.000 claims abstract description 19
- 238000009792 diffusion process Methods 0.000 claims description 16
- 230000000694 effects Effects 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 6
- 239000000941 radioactive substance Substances 0.000 claims description 6
- 238000000469 dry deposition Methods 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 8
- 238000012545 processing Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000007774 longterm Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 239000002901 radioactive waste Substances 0.000 description 1
- 239000000779 smoke Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Measurement Of Radiation (AREA)
Abstract
The invention discloses a calculation method and a device for the concentration of a multipoint source exhaust amplified gas dispersion ground surface and a radiation dose evaluation method, and belongs to the technical field of nuclear power sources. The method comprises the following steps: acquiring the distance from the evaluation subarea to the emission source S; acquiring the wind direction of the evaluation subarea blown by the emission source S according to the distance from the evaluation subarea to the emission source S; and calculating the concentration of the atmosphere dispersion ground of the evaluation subarea under the influence of the emission sources S according to the distance from the evaluation subarea to the emission sources S and the wind direction blown by the emission sources S, wherein the number of the emission sources S is a plurality of. The method can solve the problem of low accuracy of the result of the atmospheric dispersion ground concentration calculated by taking the factory site center as a single point source based on the Gaussian plume model in the related technology.
Description
Technical Field
The invention belongs to the technical field of nuclear power energy, and particularly relates to a calculation method and device for the concentration of ground with amplified gas dispersion of a multi-point source row and an evaluation method for radiation dose.
Background
During the selection of a nuclear power plant site, the radiation influence of the gas-liquid effluent of the nuclear power plant on the environment needs to be evaluated, namely, the dosage evaluation of radioactive effluent to the public and organisms is performed along with the dispersion of the atmosphere and water body into the environment. The radiation environment influence evaluation work of the nuclear power plant is taken as an important component part of the radiation protection optimization task of the newly built nuclear power plant, is a content which must be considered in the design, and also relates to various aspects of emission monitoring, operation management, radioactive waste minimization and the like after operation.
The radionuclide in the loop can enter the factory air through leakage during normal operation of the nuclear power plant, and then is discharged from the chimney through the filtering function of the factory ventilation system. The traditional gaussian plume model is used to calculate the atmospheric dispersion of a single point source, but in practice a nuclear plant site is typically composed of multiple units, each unit having a respective exhaust stack, and the exhaust conditions, exhaust concentrations, and exhaust location coordinates are all different. If the distances between a plurality of point sources of one site are far, the atmospheric dispersion is calculated by taking the site center as a single point source, so that the error of the calculation result is large.
In the environmental impact evaluation, the evaluation range is usually a range of 80 km radius centered on the factory site, 16 azimuth, 12 subareas per azimuth, and 192 fan-shaped subareas in total. The wind direction angle of each fan-shaped subarea relative to each emission source is not matched with 16 standard wind direction angles of meteorological data, and the emission distances are different, so that the accuracy of the result of the atmospheric dispersion ground concentration calculated by taking a factory site center as a single-point source based on a Gaussian smoke plume model is low, and the safety of a nuclear power plant is further influenced.
Disclosure of Invention
Aiming at the defects of the prior art, the technical problem to be solved by the invention is to provide a calculation method and a device for the concentration of the amplified gas dispersion ground of the multi-point source row and an evaluation method for radiation dosage, which can ensure that the wind direction angle of each fan-shaped evaluation subarea relative to each emission source is matched with the standard wind direction of meteorological data, and realize the accurate calculation of the concentration of the amplified gas dispersion ground of the multi-point source row.
In a first aspect, the invention provides a method for calculating the concentration of a multipoint source row amplified gas dispersion ground, which comprises the following steps: acquiring the distance from the evaluation subarea to the emission source S; acquiring the wind direction of the evaluation subarea blown by the emission source S according to the distance from the evaluation subarea to the emission source S; and calculating the concentration of the atmosphere dispersion ground of the evaluation subarea under the influence of the emission sources S according to the distance from the evaluation subarea to the emission sources S and the wind direction blown by the emission sources S, wherein the number of the emission sources S is a plurality of.
Preferably, before the acquiring the distance of the evaluation sub-region from the emission source S, the calculation method further includes: discretizing an evaluation area of the ground to obtain an azimuth index and a distance index of an evaluation subarea A (theta, gamma) in Cartesian coordinates with a factory site center O as an origin, wherein theta is the azimuth index and gamma is the distance index.
Preferably, the acquiring the distance from the evaluation sub-area to the emission source S specifically includes:
The angle between the evaluation subregion a (θ, γ) and the cartesian coordinate x-axis with the factory site center O as origin is calculated according to the following model:
the distance of the evaluation subregion a (θ, γ) to the emission source S is calculated according to the following model:
Wherein,
Alpha A is the angle between the evaluation subarea A (theta, gamma) and the X-axis of Cartesian coordinates taking the center O of the factory address as the origin;
to evaluate the x-axis component of the sub-zone a (θ, γ) from the emission source S;
to evaluate the y-axis component of the sub-zone a (θ, γ) from the emission source S;
d (gamma) is the distance from the index gamma to the center O of the plant site, wherein the index gamma corresponds to the evaluation subarea;
x S is the Cartesian coordinate X-axis component of the discharge source S with the factory site center O as the origin;
y S is the Cartesian coordinate Y-axis component of the discharge source S with the factory site center O as the origin;
to evaluate the distance of the sub-zone a (θ, γ) to the emission source S.
Preferably, the step of acquiring the wind direction of the evaluation subregion blown by the emission source S according to the distance from the evaluation subregion to the emission source S specifically includes:
The subscript of the angle between the evaluation subregion a (θ, γ) and the cartesian x-axis with the emission source S as origin is calculated according to the following model:
Converting the included angle subscript of the evaluation subarea A (theta, gamma) and the Cartesian coordinate x-axis taking the emission source S as an origin into the azimuth of the evaluation subarea A (theta, gamma) in a polar coordinate system taking the emission source S as the origin so as to match the standard azimuth in the meteorological data;
calculating the wind direction of the evaluation subarea A (theta, gamma) blown by the emission source S according to the azimuth in the polar coordinate system, wherein, An included angle subscript for evaluating the sub-area A (theta, gamma) and the Cartesian coordinate x-axis taking the emission source S as an origin; /(I)To evaluate the y-axis component of the sub-zone a (θ, γ) from the emission source S; /(I)To evaluate the distance of the sub-zone a (θ, γ) to the emission source S.
Preferably, the converting the subscript of the included angle between the evaluation sub-area a (θ, γ) and the cartesian coordinate x-axis with the emission source S as the origin to the azimuth of the evaluation sub-area a (θ, γ) in the polar coordinate system with the emission source S as the origin specifically includes:
the subscript of the angle of the evaluation subregion a (θ, γ) with the cartesian x-axis with the emission source S as the origin is converted into the orientation of the evaluation subregion a (θ, γ) in the polar coordinate system with the emission source S as the origin according to the following model:
Where θ S is the azimuth of the evaluation sub-region a (θ, γ) in a polar coordinate system with the emission source S as the origin.
Preferably, the calculating the wind direction of the evaluation subarea a (θ, γ) blown by the emission source S according to the azimuth in the polar coordinate system specifically includes:
The wind direction of the evaluation zone blown by the emission source S is calculated according to the following model:
Wherein, To evaluate the wind direction of the sub-zone a (θ, γ) when blown by the discharge source S.
Preferably, the calculating the concentration of the evaluation subarea on the atmosphere dispersion ground affected by the emission source S according to the distance from the evaluation subarea to the emission source S and the wind direction blown by the emission source S specifically includes:
The concentration of the evaluation sub-zone a (θ, γ) at the surface of the atmospheric diffusion affected by the single emission source was calculated according to the following model:
The concentration of the evaluation subregion a (θ, γ) at the surface of the atmospheric diffusion affected by the emission source S was calculated according to the following model:
wherein the number of the emission sources S is a plurality of emission sources;
To evaluate the concentration of sub-zone a (θ, γ) at the surface of the atmosphere diffusion affected by a single emission source;
x is the downwind distance;
i is a wind direction subscript, i=0 to 15, and represents 16 directions respectively;
j is an atmospheric stability category subscript, j=0 to 5, and represents a to F respectively;
k is a wind speed grade index, and k=0 to 5;
q S is the emission rate of the species emitted by emission source S;
f ijk is the combined frequency of wind direction i, stability class j and wind speed class k in the average range;
u jk is the average wind speed of the stability class j, the wind speed class k in the average range;
h is the equivalent height of the single emission source;
q A is the concentration of the evaluation subarea a (θ, γ) on the ground of the atmospheric dispersion affected by several emission sources;
attenuation factors caused by nuclide decay;
attenuation factors resulting from dry deposition of the species;
to evaluate the direction of the wind blown by the discharge source S in the sub-zone A (θ, γ);
to evaluate the distance of the sub-zone a (θ, γ) to the emission source S.
In a second aspect, the invention further provides a computing device for amplifying the gas dispersion ground concentration by the multi-point source row, which comprises a first acquisition module, a second acquisition module and a computing module.
And a first acquisition module for acquiring the distance from the evaluation zone to the emission source S. And the second acquisition module is connected with the first acquisition module and is used for acquiring the wind direction blown by the emission source S in the evaluation subarea according to the distance from the evaluation subarea to the emission source S. The calculation module is connected with the first acquisition module and the second acquisition module and is used for calculating the atmosphere dispersion ground concentration of the evaluation subarea under the influence of the emission sources S according to the distance from the evaluation subarea to the emission sources S and the wind direction blown by the emission sources S, wherein the number of the emission sources S is a plurality of.
In a third aspect, the invention also provides a computing device comprising a memory having stored therein a computer program and a processor arranged to run the computer program to implement the method of computing a multi-point source row amplified gas-dispersed ground concentration according to the first aspect.
In a fourth aspect, the present invention also provides a method for evaluating radiation dose, including: according to the method for calculating the concentration of the amplified gas dispersion ground of the multi-point source row in the first aspect, the activity of radioactive substances discharged to the environment by the nuclear power plant is obtained; evaluating the effect of radiation dose on living things based on the activity of radioactive substances emitted to the environment by nuclear power plant annual
According to the calculation method and device for the amplified gas dispersion ground concentration of the multi-point source row and the evaluation method for the radiation dose, the distances from the evaluation subarea to the emission sources in the evaluation area are obtained, the wind directions blown by the corresponding emission sources from the evaluation subarea are obtained according to the distances from the evaluation subarea to the emission sources, the directions are standard wind directions in meteorological data, and the atmospheric dispersion ground concentration of the evaluation subarea under the influence of the emission sources is calculated according to the distances from the evaluation subarea to the emission sources and the corresponding wind directions. The atmospheric dispersion ground concentration is determined according to the actual distance from the evaluation subarea to each emission source and the corresponding wind direction, and the wind direction is matched with the standard wind direction in the meteorological data, so that the accuracy of the obtained atmospheric dispersion ground concentration result is high, and the accuracy of the radiation environment evaluation result is improved.
Drawings
FIG. 1 is a schematic diagram of a method for calculating the concentration of a multipoint source row amplified gas dispersion ground in accordance with embodiment 1 of the present invention;
FIG. 2 is a schematic diagram of discrete identification of an evaluation sub-region according to embodiment 1 of the present invention;
FIG. 3 is a schematic view of a Gaussian plume model of example 1 of the present invention;
Fig. 4 is a schematic diagram of a calculation method of a multi-point source row amplified gas dispersion ground concentration according to embodiment 2 of the present invention.
Fig. 5 is a schematic structural diagram of a computing device for amplifying gas dispersion ground concentration by a multi-point source row according to embodiment 3 of the present invention.
Detailed Description
In order to make the technical scheme of the present invention better understood by those skilled in the art, the following detailed description of the embodiments of the present invention will be given with reference to the accompanying drawings.
It is to be understood that the specific embodiments and figures described herein are merely illustrative of the invention, and are not limiting of the invention.
It is to be understood that the various embodiments of the invention and the features of the embodiments may be combined with each other without conflict.
It is to be understood that only the portions relevant to the present invention are shown in the drawings for convenience of description, and the portions irrelevant to the present invention are not shown in the drawings.
It should be understood that each unit and module in the embodiments of the present invention may correspond to only one physical structure, may be formed by a plurality of physical structures, or may be integrated into one physical structure.
It will be appreciated that, without conflict, the functions and steps noted in the flowcharts and block diagrams of the present invention may occur out of the order noted in the figures.
It is to be understood that the flowcharts and block diagrams of the present invention illustrate the architecture, functionality, and operation of possible implementations of systems, apparatuses, devices, methods according to various embodiments of the present invention. Where each block in the flowchart or block diagrams may represent a unit, module, segment, code, or the like, which comprises executable instructions for implementing the specified functions. Moreover, each block or combination of blocks in the block diagrams and flowchart illustrations can be implemented by hardware-based systems that perform the specified functions, or by combinations of hardware and computer instructions.
It should be understood that the units and modules related in the embodiments of the present invention may be implemented by software, or may be implemented by hardware, for example, the units and modules may be located in a processor.
Example 1:
According to the emission characteristics of the nuclear power plant emission chimney under the normal operation condition, the embodiment assumes that the emission rate is constant, the concentrations of the evaluation subareas in the horizontal direction are equal, and the concentrations outside the evaluation subareas are zero. And according to the relative positions of the target evaluation subarea and each emission source, a key algorithm of coordinate conversion is given, and finally the atmosphere dispersion ground concentration of the multipoint sources of the evaluation subarea is given. And because the relative angles of the target evaluation subarea and each emission source are not matched with 16 orientations of the meteorological data, in order to consider conservation, the concentration of the evaluation fan-shaped subarea along an arc line is assumed to be equal, and the concentration outside the target fan-shaped subarea is zero, so that all the diffused concentrations are concentrated in a uniformly distributed fan-shaped area. Based on the assumption, the angles of the target evaluation subareas relative to all emission sources are rounded relative to multiples of 22.5 degrees, so that corresponding azimuth subscripts can be obtained, and the actual wind direction is matched with the standard wind direction angle in meteorological data.
As shown in fig. 1, the embodiment provides a method for calculating the concentration of the amplified gas dispersion ground of the multi-point source row, which includes:
step 101, the distance from the evaluation subregion to the emission source S is acquired.
In this embodiment, a plurality of units are generally installed at a plant site of a nuclear power plant, each unit is provided with an emission chimney, the types of nuclides emitted by each chimney are different, the activity of emitted nuclides is also different, and even the emission conditions (such as emission height, emission coordinates, and nearby building heights) of each chimney are different. The emission sources S refer to emission stacks, and the number of the emission sources S in this embodiment is a plurality. In the radiation environmental impact evaluation process, it is necessary to evaluate the radiation impact of the public in 192 subregions around the factory site center, and thus it is necessary to give the atmospheric dispersion factor of 192 evaluation subregions. For the target evaluation zone, the angle and distance with respect to each discharge point source need to be considered.
Optionally, before acquiring the distance from the evaluation sub-region to the emission source S, the method for calculating the point source emission atmospheric dispersion ground concentration further includes: discretizing an evaluation area of the ground to obtain an azimuth index and a distance index of an evaluation subarea A (theta, gamma) in Cartesian coordinates with a factory site center O as an origin, wherein theta is the azimuth index and gamma is the distance index.
In this embodiment, discretization processing is performed on the evaluation area of the ground in accordance with 16 orientations, 12 sub-areas for each orientation, and a total of 192 evaluation sub-areas. As shown in fig. 2, the center of the circle is the factory site center O, if the target evaluation sub-area is the evaluation sub-area a (θ, γ), θ is the azimuth index, and γ is the distance index of the evaluation sub-area. The value range of θ is 0-15, representing 16 orientations of Cartesian coordinates with the factory site center O as the origin, and the value of γ is 0-11, representing 12 sub-areas under each orientation. The discretization processing of the embodiment is suitable for other numbers of orientations and division conditions of subareas, and the discretization processing is beneficial to coordinate conversion of the evaluation subareas.
Optionally, step 101: the distance from the evaluation subarea to the emission source S is acquired, and the method specifically comprises the following steps:
The angle between the evaluation subregion a (θ, γ) and the cartesian coordinate x-axis with the factory site center O as origin is calculated according to the following model:
the distance of the evaluation subregion a (θ, γ) to the emission source S is calculated according to the following model:
Wherein,
Alpha A is the angle between the evaluation subarea A (theta, gamma) and the X-axis of Cartesian coordinates taking the center O of the factory address as the origin;
to evaluate the x-axis component of the sub-zone a (θ, γ) from the emission source S;
to evaluate the y-axis component of the sub-zone a (θ, γ) from the emission source S;
d (gamma) is the distance from the index gamma to the center O of the plant site, wherein the index gamma corresponds to the evaluation subarea;
x S is the Cartesian coordinate X-axis component of the discharge source S with the factory site center O as the origin;
y S is the Cartesian coordinate Y-axis component of the discharge source S with the factory site center O as the origin;
to evaluate the distance of the sub-zone a (θ, γ) to the emission source S.
In this embodiment, an included angle model of the evaluation subregion and the cartesian coordinate x-axis using the factory site center O as the origin is constructed according to the discretization result of the evaluation subregion. When the theta value of the included angle model is 4, alpha A is calculated to be 0, namely the evaluation subarea is positioned on the X axis of Cartesian coordinates taking the center O of the factory address as the origin. And then calculating the distance from the evaluation subarea to the emission source according to the calculated included angle between the evaluation subarea and the X-axis of the Cartesian coordinate taking the center O of the factory address as the origin.
And 102, acquiring the wind direction blown by the emission source S in the evaluation subarea according to the distance from the evaluation subarea to the emission source S.
Specifically, the method for acquiring the wind direction of the evaluation subregion blown by the emission source S according to the distance from the evaluation subregion to the emission source S comprises the following steps:
The subscript of the angle between the evaluation subregion a (θ, γ) and the cartesian x-axis with the emission source S as origin is calculated according to the following model:
Converting the included angle subscript of the evaluation subarea A (theta, gamma) and the Cartesian coordinate x-axis taking the emission source S as an origin into the azimuth of the evaluation subarea A (theta, gamma) in a polar coordinate system taking the emission source S as the origin so as to match the standard azimuth in the meteorological data;
calculating the wind direction of the evaluation subarea A (theta, gamma) blown by the emission source S according to the azimuth in the polar coordinate system, wherein, An included angle subscript for evaluating the sub-area A (theta, gamma) and the Cartesian coordinate x-axis taking the emission source S as an origin; /(I)To evaluate the y-axis component of the sub-zone a (θ, γ) from the emission source S; /(I)To evaluate the distance of the sub-zone a (θ, γ) to the emission source S.
In this embodiment, since the relative angles of the evaluation sub-area and each emission source are not matched with the 16 orientations of the meteorological data, and in order to consider conservation, it is assumed that the concentrations of the evaluation fan-shaped sub-areas along the arc line shown in fig. 2 are equal, and the concentrations outside the target fan-shaped sub-areas are zero, so that all the diffused concentrations are concentrated in the uniformly distributed fan-shaped area. Based on the above assumption, rounding the angle of the evaluation subregion relative to each emission source relative to a multiple of 22.5 ° gives the corresponding azimuth index of the evaluation subregion to simplify modeling. Therefore, an included angle subscript model of the evaluation subarea A (theta, gamma) and the Cartesian coordinate x-axis taking the emission source S as an origin is constructed based on the assumption, and the included angle subscript of the evaluation subarea A (theta, gamma) and the Cartesian coordinate x-axis taking the emission source S as the origin is obtained according to the included angle subscript model and is converted into the azimuth of the evaluation subarea A (theta, gamma) in the polar coordinate system taking the emission source S as the origin, so that the standard azimuth in the meteorological data is matched, and the accuracy of the concentration calculation result is prevented from being influenced by the mismatching of the azimuth in the meteorological data.
Optionally, converting the subscript of the included angle between the evaluation subarea a (θ, γ) and the cartesian coordinate x-axis with the emission source S as the origin to the azimuth of the evaluation subarea a (θ, γ) in the polar coordinate system with the emission source S as the origin specifically includes: the subscript of the angle of the evaluation subregion a (θ, γ) with the cartesian x-axis with the emission source S as the origin is converted into the orientation of the evaluation subregion a (θ, γ) in the polar coordinate system with the emission source S as the origin according to the following model:
Where θ s is the azimuth of the evaluation sub-region a (θ, γ) in a polar coordinate system with the emission source S as the origin.
In this embodiment, a transformation model θ S is constructed to implement matching between the index of the included angle between the evaluation sub-region a (θ, γ) and the x-axis of the cartesian coordinate with the emission source S as the origin, and the azimuth in the meteorological data, so as to complete coordinate transformation. In other words, the present embodiment calculates the azimuth counterclockwise starting from the x-axis of the cartesian coordinate with the source S as the origin, and matches the azimuth clockwise starting from the north direction in the polar coordinate system with the source S as the origin.
Optionally, the wind direction of the evaluation subarea a (θ, γ) blown by the emission source S is calculated according to the azimuth in the polar coordinate system, which specifically includes:
The wind direction of the evaluation zone blown by the emission source S is calculated according to the following model:
Wherein, To evaluate the wind direction of the sub-zone a (θ, γ) when blown by the discharge source S.
In this embodiment, the wind direction blown by the discharge source S is calculated according to the azimuth of the evaluation subregion in the polar coordinate system, so as to realize the conversion between the azimuth of the evaluation subregion and the wind direction.
And 103, calculating the concentration of the atmosphere dispersion ground of the evaluation subarea under the influence of the emission sources S according to the distance from the evaluation subarea to the emission sources S and the wind direction blown by the emission sources S, wherein the number of the emission sources S is a plurality of.
Specifically, according to the distance from the evaluation subarea to the emission source S and the wind direction blown by the emission source S, the method for calculating the concentration of the evaluation subarea on the atmosphere dispersion ground influenced by the emission source S comprises the following steps:
The concentration of the evaluation sub-zone a (θ, γ) at the surface of the atmospheric diffusion affected by the single emission source was calculated according to the following model:
The concentration of the evaluation subregion a (θ, γ) at the surface of the atmospheric diffusion affected by the emission source S was calculated according to the following model:
wherein the number of the emission sources S is a plurality of emission sources;
To evaluate the concentration of sub-zone a (θ, γ) at the surface of the atmosphere diffusion affected by a single emission source;
x is the downwind distance;
i is a wind direction subscript, i=0 to 15, and represents 16 directions respectively;
j is an atmospheric stability category subscript, j=0 to 5, and represents a to F respectively;
k is a wind speed grade index, and k=0 to 5;
q S is the emission rate of the species emitted by emission source S;
f ijk is the combined frequency of wind direction i, stability class j and wind speed class k in the average range;
u jk is the average wind speed of the stability class j, the wind speed class k in the average range;
h is the equivalent height of the single emission source;
q A is the concentration of the evaluation subarea a (θ, γ) on the ground of the atmospheric dispersion affected by several emission sources;
attenuation factors caused by nuclide decay;
attenuation factors resulting from dry deposition of the species;
to evaluate the direction of the wind blown by the discharge source S in the sub-zone A (θ, γ);
to evaluate the distance of the sub-zone a (θ, γ) to the emission source S.
In this embodiment, as shown in fig. 3, when the distance from the evaluation sub-area to the emission source S and the wind direction blown by the emission source S are obtained, the ground concentration with the leeward distance x is:
Wherein q (x) is the ground concentration with the leeward distance x;
q (x, y, 0) is the ground concentration at a downwind distance x and a horizontal distance y;
q is the emission rate of nuclear species emitted by the nuclear power plant;
u is the average wind speed;
σ z is the vertical atmospheric diffusion parameter;
σ y is the horizontal atmospheric diffusion parameter;
H is the equivalent height of the point source.
In this embodiment, the concentration of the evaluation subregion a (θ, γ) on the atmosphere dispersion ground affected by the single emission source is calculated from the model, and then the effects of the individual emission sources are summed up, the concentration of the evaluation subregion a (θ, γ) on the atmosphere dispersion ground affected by the emission source S is calculated, and the process takes into account the decay and deposition effects of the nuclide.
According to the calculation method of the multi-point source row amplified gas dispersion ground concentration, the distances from the evaluation subarea to the emission sources in the evaluation area are obtained, the wind directions blown by the corresponding emission sources from the evaluation subarea are obtained according to the distances from the evaluation subarea to the emission sources, the directions are standard wind directions in meteorological data, and the atmospheric dispersion ground concentration of the evaluation subarea under the influence of the emission sources is calculated according to the distances from the evaluation subarea to the emission sources and the corresponding wind directions. The atmospheric dispersion ground concentration is determined according to the actual distance from the evaluation subarea to each emission source and the corresponding wind direction, and the wind direction is matched with the standard wind direction in the meteorological data, so that the accuracy of the obtained atmospheric dispersion ground concentration result is high, and the accuracy of the radiation environment evaluation result is improved. Specifically, the evaluation region is processed by discretization to facilitate the subsequent coordinate conversion of the evaluation subregion. Rounding the angle of the evaluation subregion relative to the respective emission source relative to a multiple of 22.5 ° yields the corresponding orientation subscript for the evaluation subregion, thereby simplifying modeling. And (3) starting with the X axis of the Cartesian coordinate taking the discharge source S as an origin, calculating the azimuth anticlockwise, matching to the north direction in the polar coordinate system taking the discharge source S as the origin, calculating the azimuth clockwise, and completing matching of the azimuth to the standard wind direction in the meteorological data so as to improve the accuracy of the concentration calculation result. In summary, the specific calculation method of the air dispersion ground concentration of the multi-point source row is provided, the ground concentration of the target evaluation subarea is provided by considering the difference between the discharge condition, the discharge concentration and the discharge position coordinates of the multi-point source row, and a design basis is provided for plant site evaluation and radiation environment influence evaluation. .
Example 2:
As shown in fig. 4, the present embodiment provides a calculation method of a multipoint source row amplified gas dispersion ground concentration, which establishes a mathematical model of a multipoint source on the basis of a gaussian plume model, the calculation method includes:
Step 1: discretization of the evaluation subregions, for 16 orientations, 12 subregions per orientation, and a total of 192 evaluation subregions. The target subarea is an evaluation subarea A (theta, gamma), theta is a subscript of the azimuth, and gamma is a subscript of the subarea distance. The basis of the nuclear power station radiation environment influence evaluation is to give the atmospheric dispersion factors of the 192 subareas.
Step 2: the included angle between the evaluation subarea A (theta, gamma) and the Cartesian coordinate x-axis taking the factory site center O as the origin is as follows:
Alpha A -the angle of the evaluation subregion a (θ, γ) with the x-axis of the cartesian coordinates with the origin of the factory site center O.
Step 3: the distance of the evaluation subregion a (θ, γ) to the emission source S is:
-evaluating the x-axis component of the sub-zone a (θ, γ) and the emission source S;
-evaluating the y-axis component of the sub-area a (θ, γ) and the emission source S;
d (gamma) -subscript gamma corresponds to the distance of the sub-zone from the center of the site;
X S -Cartesian coordinate X-axis component of the emission source S with the factory site center O as the origin;
y S -a cartesian Y-axis component of the emission source S with the factory site center O as the origin;
the relative distance of the sub-zone a (θ, γ) from the emission source S is evaluated.
Step 4: the angle between the evaluation subarea a (θ, γ) and the cartesian coordinate x-axis with the discharge source S as the origin is given by the subscript:
-subscript of the angle of the evaluation sub-zone a (θ, γ) to the cartesian coordinate x-axis with the emission source S as origin.
Step 5: the orientation of the sub-zone a (θ, γ) in polar coordinates with the emission source S as origin is evaluated:
θ S —evaluate the orientation of the subregion a (θ, γ) in polar coordinates with the emission source S as origin.
Step 6: evaluation sub-zone A (θ, γ) in wind direction contaminated by emission source S
-Evaluating the direction of the sub-zone a (θ, γ) when it is exposed to the wind blown by the emission source a.
Step 7: since wind direction is usually given in 16 orientations each corresponding to a 22.5 sector in meteorological data, it is assumed that each angle of the same sector has the same wind frequency, i.e. is equal in the y-direction. The beneficial effects of this assumption are: although the azimuth angle of 192 subregions is an integer multiple of 22.5 °, in the evaluation of the multipoint sources, the angle of each emission source relative to the evaluation subregion is not an integer multiple of 22.5 °, so this assumption can be used to give a subscript to the azimuth by rounding off the multiple of 22.5 °. Based on this assumption, the ground concentration for a leeway distance x is:
q (x) is the ground concentration with a leeward distance x;
q (x, y, 0) is the ground concentration with the leeward distance x and the horizontal distance y;
Q-is the emission rate of nuclear species emitted by the nuclear power plant;
u-average wind speed;
sigma z -vertical atmospheric diffusion parameter;
Sigma y -horizontal atmospheric diffusion parameter;
h-equivalent height of the point source.
The long-term average concentration of the evaluation sub-zone a (θ, γ) under the influence of the emission source S is:
-evaluating the long-term average concentration of the sub-zone a (θ, γ) at the time affected by the emission source S;
x-downwind distance;
i-wind direction subscripts, i=0 to 15, representing 16 orientations respectively;
j-atmospheric stability category subscripts, j=0 to 5, respectively representing a to F;
k-wind speed grade subscript, k=0 to 5;
q S -is the emission rate of the nuclide emitted by the emission source S in the nuclear power plant;
f ijk -the combined frequency of wind direction i, stability class j and wind speed class k in the average range;
u jk -average wind speed of stability class j, wind speed class k in average range;
h-equivalent height of the point source.
Step 8: the long-term average concentration of the evaluation sub-zone a (θ, γ) at the influence of a plurality of point sources is:
q A -evaluate the long-term average concentration of sub-zone a (θ, γ) over time affected by multiple point sources.
-A decay factor caused by nuclide decay;
-a decay factor caused by dry deposition of the species.
Example 3:
As shown in fig. 5, the present embodiment provides a computing device for amplifying the concentration of the gas dispersion ground by a multi-point source row, which includes a first obtaining module 51, a second obtaining module 52, and a computing module 53.
A first acquisition module 51 for acquiring the distance of the evaluation zone to the emission source S.
The second acquisition module 52 is connected to the first acquisition module 51, and is configured to acquire a wind direction blown by the emission source S in the evaluation subregion according to a distance from the evaluation subregion to the emission source S.
The calculating module 53 is connected to the first acquiring module 51 and the second acquiring module 52, and is configured to calculate, according to a distance from the evaluation subarea to the emission source S and a wind direction blown by the emission source S, an atmospheric dispersion ground concentration of the evaluation subarea affected by the emission source S, where the number of the emission sources S is a plurality.
Optionally, the computing device further comprises a processing module. The processing module is connected with the first acquisition module and is used for discretizing an evaluation area of the ground to obtain an azimuth index and a distance index of the evaluation subarea A (theta, gamma) in Cartesian coordinates with the factory site center O as an origin, wherein theta is the azimuth index and gamma is the distance index.
Optionally, the first obtaining module is configured to calculate an angle between the evaluation sub-area a (θ, γ) and a cartesian coordinate x-axis with the factory site center O as an origin according to the following model:
the distance of the evaluation subregion a (θ, γ) to the emission source S is calculated according to the following model:
Wherein,
Alpha A is the angle between the evaluation subarea A (theta, gamma) and the X-axis of Cartesian coordinates taking the center O of the factory address as the origin;
to evaluate the x-axis component of the sub-zone a (θ, γ) from the emission source S;
to evaluate the y-axis component of the sub-zone a (θ, γ) from the emission source S;
d (gamma) is the distance from the index gamma to the center O of the plant site, wherein the index gamma corresponds to the evaluation subarea;
x S is the Cartesian coordinate X-axis component of the discharge source S with the factory site center O as the origin;
y S is the Cartesian coordinate Y-axis component of the discharge source S with the factory site center O as the origin;
to evaluate the distance of the sub-zone a (θ, γ) to the emission source S.
Optionally, the second obtaining module is configured to calculate an angle subscript between the evaluation sub-area a (θ, γ) and a cartesian coordinate x-axis with the emission source S as an origin according to the following model:
Converting the included angle subscript of the evaluation subarea A (theta, gamma) and the Cartesian coordinate x-axis taking the emission source as an origin into the azimuth of the evaluation subarea A (theta, gamma) in a polar coordinate system taking the emission source S as the origin so as to match the standard azimuth in the meteorological data;
calculating the wind direction of the evaluation subarea A (theta, gamma) blown by the emission source S according to the azimuth in the polar coordinate system, wherein, An included angle subscript for evaluating the sub-area A (theta, gamma) and the Cartesian coordinate x-axis taking the emission source S as an origin; /(I)To evaluate the y-axis component of the sub-zone a (θ, γ) from the emission source S; /(I)To evaluate the distance of the sub-zone a (θ, γ) to the emission source S.
Optionally, the second obtaining module is further configured to convert an angle subscript between the evaluation sub-area a (θ, γ) and a cartesian coordinate x-axis with the emission source S as an origin into an orientation of the evaluation sub-area a (θ, γ) in a polar coordinate system with the emission source S as an origin according to the following model:
Where θ S is the azimuth of the evaluation sub-region a (θ, γ) in a polar coordinate system with the emission source S as the origin.
Optionally, the second obtaining module is further configured to calculate the wind direction of the evaluation zone blown by the emission source S according to the following model:
Wherein, To evaluate the wind direction of the sub-zone a (θ, γ) when blown by the discharge source S.
Optionally, the calculation module is configured to calculate the concentration of the evaluation subregion a (θ, γ) at the surface of the atmospheric diffusion affected by the single emission source according to the following model:
The concentration of the evaluation subregion a (θ, γ) at the surface of the atmospheric diffusion affected by the emission source S was calculated according to the following model:
wherein the number of the emission sources S is a plurality of emission sources;
To evaluate the concentration of sub-zone a (θ, γ) at the surface of the atmosphere diffusion affected by a single emission source;
x is the downwind distance;
i is a wind direction subscript, i=0 to 15, and represents 16 directions respectively;
j is an atmospheric stability category subscript, j=0 to 5, and represents a to F respectively;
k is a wind speed grade index, and k=0 to 5;
q S is the emission rate of the species emitted by emission source S;
f ijk is the combined frequency of wind direction i, stability class j and wind speed class k in the average range;
u jk is the average wind speed of the stability class j, the wind speed class k in the average range;
h is the equivalent height of the single emission source;
q A is the concentration of the evaluation subarea a (θ, γ) on the ground of the atmospheric dispersion affected by several emission sources;
attenuation factors caused by nuclide decay;
Is the decay factor caused by dry deposition of the species.
Example 4:
The present embodiment provides an electronic device comprising a memory and a processor, the memory having stored therein a computer program, the processor being arranged to run the computer program to implement the method of calculating a multi-point source row amplified gas-dispersed ground concentration as described in embodiment 1 or embodiment 2.
Example 5:
The embodiment provides a method for evaluating radiation dose, which comprises the following steps: the method for calculating the concentration of the ground surface of the amplified gas dispersion of the multi-point source row according to the embodiment 1 or the embodiment 2 obtains the activity of the radioactive substances discharged to the environment by the nuclear power plant; the effect of radiation dose on organisms is evaluated based on the activity of the annual emitted radioactive substances to the environment of the nuclear power plant.
It is to be understood that the above embodiments are merely illustrative of the application of the principles of the present invention, but not in limitation thereof. Various modifications and improvements may be made by those skilled in the art without departing from the spirit and substance of the invention, and are also considered to be within the scope of the invention.
Claims (10)
1. The method for calculating the concentration of the amplified gas dispersion ground of the multi-point source row is characterized by comprising the following steps of:
acquiring the distance from the evaluation subarea to the emission source S;
Acquiring the wind direction of the evaluation subarea blown by the emission source S according to the distance from the evaluation subarea to the emission source S;
And calculating the concentration of the atmosphere dispersion ground of the evaluation subarea under the influence of the emission sources S according to the distance from the evaluation subarea to the emission sources S and the wind direction blown by the emission sources S, wherein the number of the emission sources S is a plurality of.
2. The method according to claim 1, further comprising, prior to said obtaining the distance of the evaluation sub-zone from the emission source S:
discretizing an evaluation area of the ground to obtain an azimuth index and a distance index of an evaluation subarea A (theta, gamma) in Cartesian coordinates with a factory site center O as an origin, wherein theta is the azimuth index and gamma is the distance index.
3. The method according to claim 2, characterized in that said obtaining the distance of the evaluation zone from the emission source S comprises in particular:
The angle between the evaluation subregion a (θ, γ) and the cartesian coordinate x-axis with the factory site center O as origin is calculated according to the following model:
the distance of the evaluation subregion a (θ, γ) to the emission source S is calculated according to the following model:
Wherein,
Alpha A is the angle between the evaluation subarea A (theta, gamma) and the X-axis of Cartesian coordinates taking the center O of the factory address as the origin;
to evaluate the x-axis component of the sub-zone a (θ, γ) from the emission source S;
to evaluate the y-axis component of the sub-zone a (θ, γ) from the emission source S;
d (gamma) is the distance from the index gamma to the center O of the plant site, wherein the index gamma corresponds to the evaluation subarea;
x S is the Cartesian coordinate X-axis component of the discharge source S with the factory site center O as the origin;
y S is the Cartesian coordinate Y-axis component of the discharge source S with the factory site center O as the origin;
to evaluate the distance of the sub-zone a (θ, γ) to the emission source S.
4. The method according to claim 2, wherein the step of obtaining the wind direction of the evaluation subregion when blown by the emission source S according to the distance from the evaluation subregion to the emission source S specifically comprises:
The subscript of the angle between the evaluation subregion a (θ, γ) and the cartesian x-axis with the emission source S as origin is calculated according to the following model:
Converting the included angle subscript of the evaluation subarea A (theta, gamma) and the Cartesian coordinate x-axis taking the emission source S as an origin into the azimuth of the evaluation subarea A (theta, gamma) in a polar coordinate system taking the emission source S as the origin so as to match the standard azimuth in the meteorological data;
Calculating the wind direction of the evaluation subarea A (theta, gamma) blown by the emission source S according to the azimuth in the polar coordinate system,
Wherein,An included angle subscript for evaluating the sub-area A (theta, gamma) and the Cartesian coordinate x-axis taking the emission source S as an origin;
to evaluate the y-axis component of the sub-zone a (θ, γ) from the emission source S;
to evaluate the distance of the sub-zone a (θ, γ) to the emission source S.
5. The method according to claim 4, wherein the converting the subscript of the included angle of the evaluation subregion a (θ, γ) with respect to the cartesian x-axis having the emission source S as the origin into the orientation of the evaluation subregion a (θ, γ) in the polar coordinate system having the emission source S as the origin specifically comprises:
the subscript of the angle of the evaluation subregion a (θ, γ) with the cartesian x-axis with the emission source S as the origin is converted into the orientation of the evaluation subregion a (θ, γ) in the polar coordinate system with the emission source S as the origin according to the following model:
Where θ s is the azimuth of the evaluation sub-region a (θ, γ) in a polar coordinate system with the emission source S as the origin.
6. The method according to claim 5, wherein said calculating the wind direction of the evaluation sub-zone a (θ, γ) when blown by the emission source S from the azimuth in the polar coordinate system, specifically comprises:
The wind direction of the evaluation zone blown by the emission source S is calculated according to the following model:
Wherein, To evaluate the wind direction of the sub-zone a (θ, γ) when blown by the discharge source S.
7. The method according to claim 2, wherein the calculating the concentration of the evaluation subarea on the atmosphere dispersion ground affected by the emission source S according to the distance from the evaluation subarea to the emission source S and the wind direction blown by the emission source S specifically comprises:
The concentration of the evaluation sub-zone a (θ, γ) at the surface of the atmospheric diffusion affected by the single emission source was calculated according to the following model:
The concentration of the evaluation subregion a (θ, γ) at the surface of the atmospheric diffusion affected by the emission source S was calculated according to the following model:
wherein the number of the emission sources S is a plurality of emission sources;
To evaluate the concentration of sub-zone a (θ, γ) at the surface of the atmosphere diffusion affected by a single emission source;
x is the downwind distance;
i is a wind direction subscript, i=0 to 15, and represents 16 directions respectively;
j is an atmospheric stability category subscript, j=0 to 5, and represents a to F respectively;
k is a wind speed grade index, and k=0 to 5;
q S is the emission rate of the species emitted by emission source S;
f ijk is the combined frequency of wind direction i, stability class j and wind speed class k in the average range;
u jk is the average wind speed of the stability class j, the wind speed class k in the average range;
h is the equivalent height of the single emission source;
q A is the concentration of the evaluation subarea a (θ, γ) on the ground of the atmospheric dispersion affected by several emission sources;
attenuation factors caused by nuclide decay;
attenuation factors resulting from dry deposition of the species;
to evaluate the direction of the wind blown by the discharge source S in the sub-zone A (θ, γ);
to evaluate the distance of the sub-zone a (θ, γ) to the emission source S.
8. A computing device for amplifying gas dispersion ground concentration by a multi-point source row is characterized by comprising a first acquisition module, a second acquisition module and a computing module,
A first acquisition module for acquiring the distance from the evaluation zone to the emission source S,
A second acquisition module connected with the first acquisition module and used for acquiring the wind direction of the evaluation subarea blown by the emission source S according to the distance from the evaluation subarea to the emission source S,
The calculation module is connected with the first acquisition module and the second acquisition module and is used for calculating the atmosphere dispersion ground concentration of the evaluation subarea under the influence of the emission sources S according to the distance from the evaluation subarea to the emission sources S and the wind direction blown by the emission sources S, wherein the number of the emission sources S is a plurality of.
9. A computing device comprising a memory and a processor, the memory having stored therein a computer program, the processor being arranged to run the computer program to implement the method of computing a multi-point source-drain amplified gas-dispersed ground concentration as claimed in any one of claims 1 to 7.
10. A method of evaluating radiation dose, comprising:
The method for calculating the concentration of the amplified gas dispersion ground of the multi-point source row according to any one of claims 1 to 7, wherein the activity of radioactive substances discharged to the environment by the nuclear power plant is obtained;
The effect of radiation dose on organisms is evaluated based on the activity of the annual emitted radioactive substances to the environment of the nuclear power plant.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410008427.XA CN117972266A (en) | 2024-01-03 | 2024-01-03 | Calculation method and device for atmospheric dispersion ground concentration and dose evaluation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410008427.XA CN117972266A (en) | 2024-01-03 | 2024-01-03 | Calculation method and device for atmospheric dispersion ground concentration and dose evaluation method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117972266A true CN117972266A (en) | 2024-05-03 |
Family
ID=90863620
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410008427.XA Pending CN117972266A (en) | 2024-01-03 | 2024-01-03 | Calculation method and device for atmospheric dispersion ground concentration and dose evaluation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117972266A (en) |
-
2024
- 2024-01-03 CN CN202410008427.XA patent/CN117972266A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102819661A (en) | New algorithm for atmospheric environment capacity by using region air quality model | |
| CN103926569B (en) | Three-dimensional radar net is based on the associated centralization of cross bearing point and interferes discrimination method with distributed compacting | |
| CN114841395A (en) | Real-time power prediction system and method based on regional building types | |
| CN117972266A (en) | Calculation method and device for atmospheric dispersion ground concentration and dose evaluation method | |
| CN111999754B (en) | Evaluation system based on nuclear facility airborne effluent monitoring data | |
| CN113236487A (en) | Wind power plant noise control method, system, device and readable storage medium | |
| CN103425875A (en) | Fast obtaining method of dose distribution of photon beam in target body | |
| CN108827279A (en) | A kind of analog sun sensor anti-interference method based on earth light | |
| CN115618718A (en) | Landfill methane emission estimation method based on surface monitoring and genetic algorithm optimization | |
| Wu et al. | On the origin of surface ozone episode in Shanghai over Yangtze River Delta during a prolonged heat wave | |
| CN107423541B (en) | Improved atmospheric dispersion calculation method for potential consequences of design benchmark accidents of nuclear power plant | |
| CN115790705A (en) | Industry cluster VOCs pollution tracing analysis method based on navigation monitoring | |
| Li et al. | Characteristics of secondary inorganic aerosols and contributions to PM2. 5 pollution based on machine learning approach in Shandong Province | |
| CN107783174B (en) | Method for calibrating detection efficiency of nuclear facility liquid effluent on-line monitor | |
| CN114722635A (en) | Layout method of atmospheric pollutant monitoring sensors in chemical industrial park | |
| CN113393014A (en) | Method and device for dividing nuclear power plant smoke plume emergency plan zone | |
| CN114441714A (en) | Power plant carbon emission intensity monitoring method and system based on vehicle-mounted measuring system | |
| CN113378367A (en) | Chemical simulation environment construction method based on complex terrain | |
| CN112990643B (en) | Design method of dosage calculation system under accident condition | |
| JP2767088B2 (en) | Estimation method of gaseous substance leakage point | |
| Miller et al. | A Mathematical Structure for Refinement of Sound Ranging Estimates. | |
| CN117273308A (en) | Regional planning method, device, system and medium for photovoltaic module | |
| Brykin et al. | Search for the Critical Point in the Calculation of the Diffusion of Radioactive Substances Emitted into the Atmosphere | |
| CN116773745B (en) | Method and device for positioning pollutant gas emission source of industrial park and electronic equipment | |
| CN116644884A (en) | Method and device for linking regional environment tracking evaluation and pollution discharge permission |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination |