CN109471147A - A kind of the gamma dose rate calculation method and system discrete based on variable weight - Google Patents
A kind of the gamma dose rate calculation method and system discrete based on variable weight Download PDFInfo
- Publication number
- CN109471147A CN109471147A CN201811115395.4A CN201811115395A CN109471147A CN 109471147 A CN109471147 A CN 109471147A CN 201811115395 A CN201811115395 A CN 201811115395A CN 109471147 A CN109471147 A CN 109471147A
- Authority
- CN
- China
- Prior art keywords
- point
- source
- calculated
- discrete
- dose rate
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/02—Dosimeters
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Measurement Of Radiation (AREA)
- Nuclear Medicine (AREA)
Abstract
The gamma dose rate calculation method and system discrete based on variable weight that the invention discloses a kind of, method includes: S1, using the geometric center in body source or face source as origin, establishes X, Y, Z coordinate axis respectively；S2, discrete weight parameter is arranged according to full-size of the distance and body source or face source of point to be calculated to each reference axis on each change in coordinate axis direction；S3, by body source or face source it is discrete be several point sources；S4, the gamma dose rate contribution for calculating each point source；S5, the gamma dose rate contribution of each point source is accumulated, obtains the gamma dose rate of body source or face source at point to be calculated；S6, the gamma dose rate of every individual sources and each face source at point to be calculated is accumulated, obtains the final gamma dose rate of point to be calculated.Method and system provided by the present invention, it is different to distance between source according to point, using different weights, keep the corresponding dispersed number in proper source and face source different, not only ensure that the precision that γ external irradiation dose rate calculates, additionally it is possible to greatly improve the efficiency of calculating.
Description
Technical field
The present invention relates to radiation protection technique fields, and in particular to a kind of gamma dose rate calculating side discrete based on variable weight
Method and system.
Background technique
Point- kernel integral method is to calculate the common method of γ external irradiation dose, is widely used in SHIELDING CALCULATION and external exposure
In Rapid Dose Calculation, compared with Monte Carlo method, it is not limited not only by bulk and shielding body thickness, but also is measured when machine
It is small, there is faster calculating speed, disadvantage is that this method cannot carry out accurate consideration to KPT Scatter problem, needs to introduce
Build up factor is modified KPT Scatter, and calculated result is often relatively conservative.
This method uses discrete method, conversion based on the calculating of isotopic point source, for body source or face source
For several isotopic point sources, the calculated result of each point source is added up, to obtain whole individual sources or face source
Calculated result.Body source or face source are separated into the number of point source, it is not only related with the precision of calculated result, also with required meter
Evaluation time is related.On the one hand traditional Point- kernel integral algorithm does not optimize computational efficiency at present, and calculating speed is already not
It is able to satisfy current demand, on the other hand as advantage of the Point- kernel integral method in gamma radiation field calculating shows, while with
The needs of radiation field visualization display, the calculative dot matrix quantity for representing three dimensional radiation field is more and more huger, this is to calculating
Speed propose stern challenge.
Summary of the invention
In view of the deficiencies in the prior art, the γ agent discrete based on variable weight that the purpose of the present invention is to provide a kind of
Dose rate calculation method and system, according to calculating point to the distance between source difference, using different weights to body source or face source into
Row is discrete, guarantees to improve the efficiency that gamma radiation field calculates under the premise of computational accuracy.
To achieve the above object, The technical solution adopted by the invention is as follows:
A kind of gamma dose rate calculation method discrete based on variable weight, comprising:
S1, local coordinate system is established inside integrated source or one side source, with the geometric center in the body source or the face source
For origin, X, Y, Z coordinate axis are established respectively；
S2, calculate a point to be calculated arrive each reference axis distance, according to the distance of the point to be calculated to each reference axis with
Discrete weight parameter is arranged in the full-size of the body source or the face source on each change in coordinate axis direction；
S3, according to the point to be calculated to the distance of each reference axis, the full-size and the discrete weight parameter,
The body source or the face source is discrete for several point sources；
S4, the gamma dose rate contribution that each point source is calculated according to the distance of each point source to the point to be calculated；
S5, the gamma dose rate contribution of each point source is accumulated, obtains the body source or the face source described wait count
Calculate the gamma dose rate at point；
S6, the gamma dose rate of every individual sources and each face source at the point to be calculated is accumulated, obtain it is described to
Calculate the final gamma dose rate of point.
Further, a kind of gamma dose rate calculation method discrete based on variable weight as described above, step S3 include:
S31, according to the point to be calculated to the distance of each reference axis, the full-size and the discrete weight parameter,
Dispersion number of the point to be calculated on each change in coordinate axis direction is calculated；
S32, from the maximum dispersion number on preset each change in coordinate axis direction and the point to be calculated is calculated in each coordinate
Smaller is taken in dispersion number in axis direction, as final dispersion number of the point to be calculated on each change in coordinate axis direction；
S33, by the body source or the face source it is discrete be several point sources, point source number be the point to be calculated in each coordinate
The product of final dispersion number in axis direction.
Further, a kind of gamma dose rate calculation method discrete based on variable weight as described above, step S31 are specifically wrapped
It includes:
If Ri≠ 0, then the point to be calculated is obtained after calculating and round up according to the following formula on each change in coordinate axis direction
Dispersion number Ni；
Wherein, n is the discrete weight parameter, DiFor the maximum of the body source or the face source on each change in coordinate axis direction
Size, RiDistance for the point to be calculated to each reference axis, NiFor dispersion number of the point to be calculated on i change in coordinate axis direction,
NiFor positive integer, i=i | x, y, z }；
If Ri=0, i.e., the described point to be calculated is located in reference axis, then institute is obtained after calculating and round up according to the following formula
State dispersion number N of the point to be calculated on the change in coordinate axis direction；
R '=R-D/2,
Wherein, n is the discrete weight parameter, and R is the point to be calculated to the distance of origin, and D is the body source or institute
Full-size of the face source on the change in coordinate axis direction is stated, N is positive integer.
Further, a kind of gamma dose rate calculation method discrete based on variable weight as described above, in step S2, according to institute
State point to be calculated to each reference axis distance and the body source or the face source on each change in coordinate axis direction full-size setting
Discrete weight parameter n, specifically includes:
IfThen the value of n is set as
IfThen the value of n is set as the arbitrary number between 10-20 according to the actual situation；
IfThen the value of n is set as the arbitrary number between 1-10 according to the actual situation；
Wherein, DiFor the full-size of the body source or the face source on each change in coordinate axis direction, RiFor the point to be calculated
To the distance of each reference axis, i=i | x, y, z }.
Further, a kind of gamma dose rate calculation method discrete based on variable weight as described above, step S4 are specifically included:
The gamma dose rate contribution of each point source is calculated according to the following formula；
Wherein, F (E) is transfer function, and E is the energy of γ photon, AiFor the activity of i-th of point source, B is build up factor, e-μtFor exponential damping item, μ is line attenuation coefficient of the material to γ photon, and t is the distance that γ photon passes through within this material, and r is
Distance of i-th of point source to the point to be calculated, HiFor the gamma dose rate contribution of i-th of point source, i is positive integer.
A kind of gamma dose rate computing system discrete based on variable weight is additionally provided in the embodiment of the present invention, comprising:
Module is established, for establishing local coordinate system inside integrated source or one side source, with the body source or the face source
Geometric center be origin, establish X, Y, Z coordinate axis respectively；
Setup module, the distance for calculating a point to be calculated to each reference axis, according to the point to be calculated to each coordinate
Discrete weight parameter is arranged in the full-size of the distance of axis and the body source or the face source on each change in coordinate axis direction；
Discrete block, for according to the point to be calculated to the distance of each reference axis, the full-size and it is described from
Weight parameter is dissipated, the body source or the face source is discrete for several point sources；
Computing module calculates the gamma dose rate tribute of each point source for the distance according to each point source to the point to be calculated
It offers；
First accumulation module obtains the body source or described for accumulating the gamma dose rate contribution of each point source
Gamma dose rate of the face source at the point to be calculated；
Second accumulation module, it is tired for carrying out the gamma dose rate of every individual sources and each face source at the point to be calculated
Product, obtains the final gamma dose rate of the point to be calculated.
Further, a kind of gamma dose rate computing system discrete based on variable weight as described above, the discrete block packet
It includes:
Computational submodule, for according to the point to be calculated to the distance of each reference axis, the full-size and it is described from
Weight parameter is dissipated, dispersion number of the point to be calculated on each change in coordinate axis direction is calculated；
Determine submodule, for from preset each change in coordinate axis direction maximum dispersion number and be calculated it is described to be calculated
Point smaller is taken in the dispersion number on each change in coordinate axis direction, as the point to be calculated on each change in coordinate axis direction it is final from
Dissipate number；
Discrete submodule, for by the body source or the face source it is discrete be several point sources, point source number is described wait count
Calculate the product of final dispersion number of the point on each change in coordinate axis direction.
Further, a kind of gamma dose rate computing system discrete based on variable weight as described above, the computational submodule
It is specifically used for:
If Ri≠ 0, then the point to be calculated is obtained after calculating and round up according to the following formula on each change in coordinate axis direction
Dispersion number Ni；
Wherein, n is the discrete weight parameter, DiFor the maximum of the body source or the face source on each change in coordinate axis direction
Size, RiDistance for the point to be calculated to each reference axis, NiFor dispersion number of the point to be calculated on i change in coordinate axis direction,
NiFor positive integer, i=i | x, y, z }；
If Ri=0, i.e., the described point to be calculated is located in reference axis, then institute is obtained after calculating and round up according to the following formula
State dispersion number N of the point to be calculated on the change in coordinate axis direction；
R '=R-D/2,
Wherein, n is the discrete weight parameter, and R is the point to be calculated to the distance of origin, and D is the body source or institute
Full-size of the face source on the change in coordinate axis direction is stated, N is positive integer.
Further, a kind of gamma dose rate computing system discrete based on variable weight as described above, the setup module tool
Body is used for:
IfThen the value of n is set as
IfThen the value of n is set as the arbitrary number between 10-20 according to the actual situation；
IfThen the value of n is set as the arbitrary number between 1-10 according to the actual situation；
Wherein, DiFor the full-size of the body source or the face source on each change in coordinate axis direction, RiFor the point to be calculated
To the distance of each reference axis, i=i | x, y, z }.
Further, a kind of gamma dose rate computing system discrete based on variable weight as described above, the first accumulation mould
Block is specifically used for:
The gamma dose rate contribution of each point source is calculated according to the following formula；
Wherein, F (E) is transfer function, and E is the energy of γ photon, AiFor the activity of i-th of point source, B is build up factor, e-μtFor exponential damping item, μ is line attenuation coefficient of the material to γ photon, and t is the distance that γ photon passes through within this material, and r is
Distance of i-th of point source to the point to be calculated, HiFor the gamma dose rate contribution of i-th of point source, i is positive integer.
The beneficial effects of the present invention are: method and system provided by the present invention, no according to the distance between point to source
Together, using different weights, keep the corresponding dispersed number in proper source and face source different, the calculating for distance sources compared with far region
Point, less dispersed number can reach very high computational accuracy, and the closer region of distance sources, need biggish dispersed number
It can be only achieved corresponding computational accuracy.It not only ensure that the precision that γ external irradiation dose rate calculates, additionally it is possible to substantially increase meter
The efficiency of calculation.
Detailed description of the invention
Fig. 1 is that a kind of process of the gamma dose rate calculation method discrete based on variable weight provided in the embodiment of the present invention is shown
It is intended to；
Fig. 2 is that a kind of structure of the gamma dose rate computing system discrete based on variable weight provided in the embodiment of the present invention is shown
It is intended to.
Specific embodiment
The present invention is described in further detail with specific embodiment with reference to the accompanying drawings of the specification.
As shown in Figure 1, a kind of gamma dose rate calculation method discrete based on variable weight, comprising:
S1, local coordinate system is established inside integrated source or one side source, using the geometric center in body source or face source as origin, point
X, Y, Z coordinate axis are not established；
The distance of S2,1 point to be calculated of calculating to each reference axis, according to the distance of point to be calculated to each reference axis and body source
Or discrete weight parameter is arranged in full-size of the face source on each change in coordinate axis direction；
S3, according to point to be calculated to the distance of each reference axis, full-size and discrete weight parameter, by body source or face source
Discrete is several point sources；
S4, the gamma dose rate contribution that each point source is calculated according to the distance of each point source to point to be calculated；
S5, the gamma dose rate contribution of each point source is accumulated, obtains the γ agent of body source or face source at point to be calculated
Dose rate；
S6, the gamma dose rate of every individual sources and each face source at point to be calculated is accumulated, obtains point to be calculated
Final gamma dose rate.
Step S3 includes:
S31, according to point to be calculated to the distance of each reference axis, full-size and discrete weight parameter, be calculated wait count
Calculate dispersion number of the point on each change in coordinate axis direction；
S32, from the maximum dispersion number on preset each change in coordinate axis direction and point to be calculated is calculated in each reference axis side
Smaller is taken in upward dispersion number, as final dispersion number of the point to be calculated on each change in coordinate axis direction；
S33, by body source or face source it is discrete be several point sources, point source number be point to be calculated on each change in coordinate axis direction most
The product of whole dispersion number.
Step S31 is specifically included:
If Ri≠ 0, then obtained after calculating and round up according to the following formula point to be calculated on each change in coordinate axis direction from
Dissipate number Ni；
Wherein, n is discrete weight parameter, DiFor the full-size of body source or face source on each change in coordinate axis direction, RiFor wait count
Calculate the distance that point arrives each reference axis, NiFor dispersion number of the point to be calculated on i change in coordinate axis direction, NiFor positive integer, i=i | x, y,
z}；
If Ri=0, i.e., point to be calculated is located in reference axis, then obtains after calculating and round up according to the following formula to be calculated
Dispersion number N of the point on the change in coordinate axis direction；
R '=R-D/2,
Wherein, n is discrete weight parameter, and R is distance of the point to be calculated to origin, and D is body source or face source in the reference axis
Full-size on direction, N are positive integer.
In step S2, according to the distance of point to be calculated to each reference axis and body source or face source on each change in coordinate axis direction most
Discrete weight parameter n is arranged in large scale, specifically includes:
IfThen the value of n is set as
IfThen the value of n is set as the arbitrary number between 10-20 according to the actual situation；
IfThen the value of n is set as the arbitrary number between 1-10 according to the actual situation；
Wherein, DiFor the full-size of body source or face source on each change in coordinate axis direction, RiFor point to be calculated to each reference axis
Distance, i=i | x, y, z }.
Step S4 is specifically included:
The gamma dose rate contribution of each point source is calculated according to the following formula；
Wherein, F (E) is transfer function, and E is the energy of γ photon, AiFor the activity of i-th of point source, B is build up factor, e-μtFor exponential damping item, μ is line attenuation coefficient of the material to γ photon, and t is the distance that γ photon passes through within this material, and r is
Distance of i-th of point source to point to be calculated, HiFor the gamma dose rate contribution of i-th of point source, i is positive integer.
Body source or face source dispersed number are more, and obtained radiation field data is also more accurate.With the increasing of dispersed number
More, the error between calculated value and analytic solutions can be smaller and smaller, and when dispersed number increase to a certain extent after, calculated value with
Error between analytic solutions tends to 0, further increases dispersed number, and computational accuracy will not get a promotion, instead can be significantly
Increase calculate needed for time, while with calculate point and source between distance increase, when calculated value tends towards stability it is required from
Scattered number can also gradually become smaller, therefore comprehensively consider computational accuracy and time, use weight discrete method herein, i.e., for distance
The calculating point of body source or face source immediate area, using biggish dispersed number；For the meter apart from body source or face source compared with far region
Point is calculated, using lesser dispersed number.It is discrete using weight for body source herein in computation of radiation field, it can be according to practical meter
It calculates and needs, comprehensively consider calculating point to the distance on source surface and the full-size in source itself, select suitable discrete weight parameter
n。
Embodiment one
The geometric dimension of one staircase is 6.3m*4.25m*9.3m, and will be drawn at equal intervals in the place by length, width and height simultaneously
It is divided into 20*15*31 totally 9300 gamma dose rates point to be calculated.The body source is evenly distributed in pipeline there are a pipe surface source
Wall, activity 109Bq, energy 1.173MeV, pipeline outer radius are 20cm, inside radius 19cm, are highly 200cm, are located at
On second layer step, ignore influence of wall thickness.It is discrete using variable weight to pipe surface source item after method of the invention,
Discrete weight parameter n in short transverse takes 20, on same computer, the calculating of calculation method and conventional method after optimization
Time comparison is as shown in the table.After variable weight discrete optimizing method, the calculating time is reduced to about original 1/3.
Meanwhile in order to compare the calculated result of the two, in source, near zone has chosen 600 calculating points, between the two
Average deviation be 1.32%, maximum deviation 4.02%；Distance sources compared with far region choose 900 calculating points, the two it
Between average deviation be 0.61%, maximum deviation 3.55%.The calculating deviation of this degree meets scene γ radiation enough
The calculating requirement of field.
As can be seen that the method for computational efficiency optimization proposed in this paper, can not only effectively improve the meter of gamma radiation field
Efficiency is calculated, and can correspondingly guarantee the computational accuracy of result.
As shown in Fig. 2, a kind of gamma dose rate computing system discrete based on variable weight is additionally provided in the embodiment of the present invention,
Include:
Module 1 is established, for establishing local coordinate system inside integrated source or one side source, in the geometry in body source or face source
The heart is origin, establishes X, Y, Z coordinate axis respectively；
Setup module 2, the distance for calculating a point to be calculated to each reference axis, according to point to be calculated to each reference axis
Discrete weight parameter is arranged in the full-size of distance and body source or face source on each change in coordinate axis direction；
Discrete block 3, for according to point to be calculated to the distance of each reference axis, full-size and discrete weight parameter,
Body source or face source is discrete for several point sources；
Computing module 4 calculates the gamma dose rate contribution of each point source for the distance according to each point source to point to be calculated；
First accumulation module 5, for the contribution of the gamma dose rate of each point source to be accumulated, obtain body source or face source to
Calculate the gamma dose rate at point；
Second accumulation module 6, for the gamma dose rate of every individual sources and each face source at point to be calculated to be accumulated,
Obtain the final gamma dose rate of point to be calculated.
Discrete block 3 includes:
Computational submodule, for according to point to be calculated to the distance of each reference axis, full-size and discrete weight parameter, meter
Calculation obtains dispersion number of the point to be calculated on each change in coordinate axis direction；
Submodule is determined, for from the maximum dispersion number on preset each change in coordinate axis direction and point to be calculated being calculated existing
Smaller is taken in dispersion number on each change in coordinate axis direction, as final dispersion number of the point to be calculated on each change in coordinate axis direction；
Discrete submodule, for by body source or face source it is discrete be several point sources, point source number be point to be calculated in each coordinate
The product of final dispersion number in axis direction.
Computational submodule is specifically used for:
If Ri≠ 0, then obtained after calculating and round up according to the following formula point to be calculated on each change in coordinate axis direction from
Dissipate number Ni；
Wherein, n is discrete weight parameter, DiFor the full-size of body source or face source on each change in coordinate axis direction, RiFor wait count
Calculate the distance that point arrives each reference axis, NiFor dispersion number of the point to be calculated on i change in coordinate axis direction, NiFor positive integer, i=i | x, y,
z}；
If Ri=0, i.e., point to be calculated is located in reference axis, then obtains after calculating and round up according to the following formula to be calculated
Dispersion number N of the point on the change in coordinate axis direction；
R '=R-D/2,
Wherein, n is discrete weight parameter, and R is distance of the point to be calculated to origin, and D is body source or face source in the reference axis
Full-size on direction, N are positive integer.
Setup module 2 is specifically used for:
IfThen the value of n is set as
IfThen the value of n is set as the arbitrary number between 10-20 according to the actual situation；
IfThen the value of n is set as the arbitrary number between 1-10 according to the actual situation；
Wherein, DiFor the full-size of body source or face source on each change in coordinate axis direction, RiFor point to be calculated to each reference axis
Distance, i=i | x, y, z }.
First accumulation module 5 is specifically used for:
The gamma dose rate contribution of each point source is calculated according to the following formula；
Wherein, F (E) is transfer function, and E is the energy of γ photon, AiFor the activity of i-th of point source, B is build up factor, e-μtFor exponential damping item, μ is line attenuation coefficient of the material to γ photon, and t is the distance that γ photon passes through within this material, and r is
Distance of i-th of point source to point to be calculated, HiFor the gamma dose rate contribution of i-th of point source, i is positive integer.
Obviously, various changes and modifications can be made to the invention without departing from essence of the invention by those skilled in the art
Mind and range.In this way, if these modifications and changes of the present invention belongs to the range of the claims in the present invention and its equivalent technology
Within, then the present invention is also intended to include these modifications and variations.
Claims (10)
1. a kind of gamma dose rate calculation method discrete based on variable weight characterized by comprising
S1, local coordinate system is established inside integrated source or one side source, is original with the geometric center in the body source or the face source
Point establishes X, Y, Z coordinate axis respectively；
S2, calculate a point to be calculated to each reference axis distance, according to the distance of the point to be calculated to each reference axis and described
Discrete weight parameter is arranged in the full-size of body source or the face source on each change in coordinate axis direction；
S3, according to the point to be calculated to the distance of each reference axis, the full-size and the discrete weight parameter, by institute
It is discrete for several point sources to state body source or the face source；
S4, the gamma dose rate contribution that each point source is calculated according to the distance of each point source to the point to be calculated；
S5, the gamma dose rate contribution of each point source is accumulated, obtains the body source or the face source in the point to be calculated
The gamma dose rate at place；
S6, the gamma dose rate of every individual sources and each face source at the point to be calculated is accumulated, is obtained described to be calculated
The final gamma dose rate of point.
2. a kind of gamma dose rate calculation method discrete based on variable weight according to claim 1, which is characterized in that step
S3 includes:
S31, according to the point to be calculated to the distance of each reference axis, the full-size and the discrete weight parameter, calculate
Obtain dispersion number of the point to be calculated on each change in coordinate axis direction；
S32, from the maximum dispersion number on preset each change in coordinate axis direction and the point to be calculated is calculated in each reference axis side
Smaller is taken in upward dispersion number, as final dispersion number of the point to be calculated on each change in coordinate axis direction；
S33, by the body source or the face source it is discrete be several point sources, point source number be the point to be calculated in each reference axis side
The product of upward final dispersion number.
3. a kind of gamma dose rate calculation method discrete based on variable weight according to claim 2, which is characterized in that step
S31 is specifically included:
If Ri≠ 0, then it is discrete on each change in coordinate axis direction that the point to be calculated is obtained after calculating and round up according to the following formula
Number Ni；
Wherein, n is the discrete weight parameter, DiFor the full-size of the body source or the face source on each change in coordinate axis direction,
RiDistance for the point to be calculated to each reference axis, NiFor dispersion number of the point to be calculated on i change in coordinate axis direction, NiFor
Positive integer, i=i | x, y, z }；
If Ri=0, i.e., the described point to be calculated is located in reference axis, then obtained after calculating and round up according to the following formula it is described to
Calculate dispersion number N of the point on the change in coordinate axis direction；
Wherein, n is the discrete weight parameter, and R is distance of the point to be calculated to origin, and D is the body source or the face
Full-size of the source on the change in coordinate axis direction, N are positive integer.
4. a kind of gamma dose rate calculation method discrete based on variable weight according to claim 3, which is characterized in that step
In S2, according to the distance of the point to be calculated to each reference axis and the body source or the face source on each change in coordinate axis direction most
Discrete weight parameter n is arranged in large scale, specifically includes:
IfThen the value of n is set as
IfThen the value of n is set as the arbitrary number between 10-20 according to the actual situation；
IfThen the value of n is set as the arbitrary number between 1-10 according to the actual situation；
Wherein, DiFor the full-size of the body source or the face source on each change in coordinate axis direction, RiIt is the point to be calculated to respectively
The distance of reference axis, i=i | x, y, z }.
5. a kind of gamma dose rate calculation method discrete based on variable weight according to claim 1-4, feature
It is, step S4 is specifically included:
The gamma dose rate contribution of each point source is calculated according to the following formula；
Wherein, F (E) is transfer function, and E is the energy of γ photon, AiFor the activity of i-th of point source, B is build up factor, e-μtFor
Exponential damping item, the distance that μ passes through the line attenuation coefficient of γ photon, t for γ photon for material within this material, i-th of r
Distance of the point source to the point to be calculated, HiFor the gamma dose rate contribution of i-th of point source, i is positive integer.
6. a kind of gamma dose rate computing system discrete based on variable weight characterized by comprising
Module is established, for establishing local coordinate system inside integrated source or one side source, with the several of the body source or the face source
What center is origin, establishes X, Y, Z coordinate axis respectively；
Setup module, the distance for calculating a point to be calculated to each reference axis, according to the point to be calculated to each reference axis
Discrete weight parameter is arranged in the full-size of distance and the body source or the face source on each change in coordinate axis direction；
Discrete block, for according to the point to be calculated to the distance of each reference axis, the full-size and the discrete power
Weight parameter, the body source or the face source is discrete for several point sources；
Computing module calculates the gamma dose rate contribution of each point source for the distance according to each point source to the point to be calculated；
First accumulation module obtains the body source or the face source for accumulating the gamma dose rate contribution of each point source
Gamma dose rate at the point to be calculated；
Second accumulation module, for the gamma dose rate of every individual sources and each face source at the point to be calculated to be accumulated,
Obtain the final gamma dose rate of the point to be calculated.
7. a kind of gamma dose rate computing system discrete based on variable weight according to claim 6, which is characterized in that described
Discrete block includes:
Computational submodule, for according to the point to be calculated to the distance of each reference axis, the full-size and the discrete power
Weight parameter, is calculated dispersion number of the point to be calculated on each change in coordinate axis direction；
Submodule is determined, for from the maximum dispersion number on preset each change in coordinate axis direction and the point to be calculated being calculated existing
Smaller is taken in dispersion number on each change in coordinate axis direction, it is final discrete on each change in coordinate axis direction as the point to be calculated
Number；
Discrete submodule, for by the body source or the face source it is discrete be several point sources, point source number be the point to be calculated
The product of final dispersion number on each change in coordinate axis direction.
8. a kind of gamma dose rate computing system discrete based on variable weight according to claim 7, which is characterized in that described
Computational submodule is specifically used for:
If Ri≠ 0, then it is discrete on each change in coordinate axis direction that the point to be calculated is obtained after calculating and round up according to the following formula
Number Ni；
Wherein, n is the discrete weight parameter, DiFor the full-size of the body source or the face source on each change in coordinate axis direction,
RiDistance for the point to be calculated to each reference axis, NiFor dispersion number of the point to be calculated on i change in coordinate axis direction, NiFor
Positive integer, i=i | x, y, z }；
If Ri=0, i.e., the described point to be calculated is located in reference axis, then obtained after calculating and round up according to the following formula it is described to
Calculate dispersion number N of the point on the change in coordinate axis direction；
Wherein, n is the discrete weight parameter, and R is distance of the point to be calculated to origin, and D is the body source or the face
Full-size of the source on the change in coordinate axis direction, N are positive integer.
9. a kind of gamma dose rate computing system discrete based on variable weight according to claim 8, which is characterized in that described
Setup module is specifically used for:
IfThen the value of n is set as
IfThen the value of n is set as the arbitrary number between 10-20 according to the actual situation；
IfThen the value of n is set as the arbitrary number between 1-10 according to the actual situation；
Wherein, DiFor the full-size of the body source or the face source on each change in coordinate axis direction, RiIt is the point to be calculated to respectively
The distance of reference axis, i=i | x, y, z }.
10. according to a kind of described in any item gamma dose rate computing systems discrete based on variable weight of claim 6-9, feature
It is, the first accumulation module is specifically used for:
The gamma dose rate contribution of each point source is calculated according to the following formula；
Wherein, F (E) is transfer function, and E is the energy of γ photon, AiFor the activity of i-th of point source, B is build up factor, e-μtFor
Exponential damping item, the distance that μ passes through the line attenuation coefficient of γ photon, t for γ photon for material within this material, i-th of r
Distance of the point source to the point to be calculated, HiFor the gamma dose rate contribution of i-th of point source, i is positive integer.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811115395.4A CN109471147B (en) | 2018-09-25 | 2018-09-25 | Gamma dose rate calculation method and system based on variable weight dispersion |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811115395.4A CN109471147B (en) | 2018-09-25 | 2018-09-25 | Gamma dose rate calculation method and system based on variable weight dispersion |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109471147A true CN109471147A (en) | 2019-03-15 |
| CN109471147B CN109471147B (en) | 2022-10-18 |
Family
ID=65663195
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811115395.4A Active CN109471147B (en) | 2018-09-25 | 2018-09-25 | Gamma dose rate calculation method and system based on variable weight dispersion |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109471147B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110221338A (en) * | 2019-05-17 | 2019-09-10 | 华南理工大学 | A method of reconstruct shield neighboring area radiation field |
Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0824697A1 (en) * | 1995-04-18 | 1998-02-25 | Wyntek Diagnostics, Inc. | One step immunochromatographic device and method of use |
| WO2003046611A1 (en) * | 2001-11-27 | 2003-06-05 | Bnfl (Ip) Ltd | Environmental radiation detector |
| US20080230705A1 (en) * | 2004-11-09 | 2008-09-25 | Spectrum Dynamics Llc | Radioimaging |
| EP2167971A2 (en) * | 2007-07-19 | 2010-03-31 | bioMérieux | Ezrin assay method for the in vitro diagnosis of colorectal cancer |
| WO2011019303A1 (en) * | 2009-08-13 | 2011-02-17 | Закрытое Акционерное Общество "Hayчнo-Пpoизвoдcтвeннaя Коммерческая Фирма "Элtah Лtд" | Multi-element x-ray radiation detector, rare earth x-ray luminophore therefor, and method for forming a multi-element scintillator and detector as a whole |
| EP2556389A1 (en) * | 2010-04-09 | 2013-02-13 | Isp Investments Inc. | Radiation dosimetry method |
| US20130173220A1 (en) * | 2006-03-17 | 2013-07-04 | Canberra Industries, Inc. | Radiation Analysis System and Method |
| CN104483693A (en) * | 2014-12-24 | 2015-04-01 | 西北核技术研究所 | Non-uniform distributed source detection efficiency calculation and simulation device and method |
| CN105092612A (en) * | 2014-05-06 | 2015-11-25 | 阿斯特菲公司 | Computed tomography system for cargo and transported containers |
| WO2016012620A1 (en) * | 2014-07-24 | 2016-01-28 | Dosevue Nv | Direct surface radiation dose measurement system with quantitative optical read-out |
| CN106061541A (en) * | 2014-02-27 | 2016-10-26 | 皇家飞利浦有限公司 | System for applying radiation to a target region within a subject |
| CN106199676A (en) * | 2015-04-30 | 2016-12-07 | 北京中智核安科技有限公司 | A kind of gamma detector passive efficiency scale new method |
| CN106814384A (en) * | 2015-11-27 | 2017-06-09 | 华北电力大学 | Nuclear power plant's point source radiation source strength backstepping method and point source radiation source strength backstepping system |
| CN106932810A (en) * | 2017-04-01 | 2017-07-07 | 西安体医疗科技有限公司 | A kind of convolutional calculation method of gamma rays dosage |
| CN107290769A (en) * | 2016-04-12 | 2017-10-24 | 华北电力大学 | The recombination radiation source strength backstepping method and system of nuclear power plant's point source body source combination |
| CN108287357A (en) * | 2018-01-15 | 2018-07-17 | 东华理工大学 | A kind of source peak detection efficient acquisition methods of cylinder bulk detector to point source |
| CN108549753A (en) * | 2018-03-28 | 2018-09-18 | 中国船舶重工集团公司第七〇九研究所 | A kind of radiation shield computational methods that Point- kernel integral method is coupled with Monte Carlo method |
| CA3076763A1 (en) * | 2017-10-20 | 2019-04-25 | Australian Nuclear Science And Technology Organisation | Compressive imaging method and system |
-
2018
- 2018-09-25 CN CN201811115395.4A patent/CN109471147B/en active Active
Patent Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0824697A1 (en) * | 1995-04-18 | 1998-02-25 | Wyntek Diagnostics, Inc. | One step immunochromatographic device and method of use |
| WO2003046611A1 (en) * | 2001-11-27 | 2003-06-05 | Bnfl (Ip) Ltd | Environmental radiation detector |
| US20080230705A1 (en) * | 2004-11-09 | 2008-09-25 | Spectrum Dynamics Llc | Radioimaging |
| US20130173220A1 (en) * | 2006-03-17 | 2013-07-04 | Canberra Industries, Inc. | Radiation Analysis System and Method |
| EP2167971A2 (en) * | 2007-07-19 | 2010-03-31 | bioMérieux | Ezrin assay method for the in vitro diagnosis of colorectal cancer |
| WO2011019303A1 (en) * | 2009-08-13 | 2011-02-17 | Закрытое Акционерное Общество "Hayчнo-Пpoизвoдcтвeннaя Коммерческая Фирма "Элtah Лtд" | Multi-element x-ray radiation detector, rare earth x-ray luminophore therefor, and method for forming a multi-element scintillator and detector as a whole |
| EP2556389A1 (en) * | 2010-04-09 | 2013-02-13 | Isp Investments Inc. | Radiation dosimetry method |
| CN106061541A (en) * | 2014-02-27 | 2016-10-26 | 皇家飞利浦有限公司 | System for applying radiation to a target region within a subject |
| CN105092612A (en) * | 2014-05-06 | 2015-11-25 | 阿斯特菲公司 | Computed tomography system for cargo and transported containers |
| WO2016012620A1 (en) * | 2014-07-24 | 2016-01-28 | Dosevue Nv | Direct surface radiation dose measurement system with quantitative optical read-out |
| CN104483693A (en) * | 2014-12-24 | 2015-04-01 | 西北核技术研究所 | Non-uniform distributed source detection efficiency calculation and simulation device and method |
| CN106199676A (en) * | 2015-04-30 | 2016-12-07 | 北京中智核安科技有限公司 | A kind of gamma detector passive efficiency scale new method |
| CN106814384A (en) * | 2015-11-27 | 2017-06-09 | 华北电力大学 | Nuclear power plant's point source radiation source strength backstepping method and point source radiation source strength backstepping system |
| CN107290769A (en) * | 2016-04-12 | 2017-10-24 | 华北电力大学 | The recombination radiation source strength backstepping method and system of nuclear power plant's point source body source combination |
| CN106932810A (en) * | 2017-04-01 | 2017-07-07 | 西安体医疗科技有限公司 | A kind of convolutional calculation method of gamma rays dosage |
| CA3076763A1 (en) * | 2017-10-20 | 2019-04-25 | Australian Nuclear Science And Technology Organisation | Compressive imaging method and system |
| CN108287357A (en) * | 2018-01-15 | 2018-07-17 | 东华理工大学 | A kind of source peak detection efficient acquisition methods of cylinder bulk detector to point source |
| CN108549753A (en) * | 2018-03-28 | 2018-09-18 | 中国船舶重工集团公司第七〇九研究所 | A kind of radiation shield computational methods that Point- kernel integral method is coupled with Monte Carlo method |
Non-Patent Citations (10)
| Title |
|---|
| C. M. HOPPER.ET: "Slide Rule for Rapid Response Estimation of Radiological Dose from Criticality Accidents", 《COMPUTATIONAL PHYSICS AND ENGINEERING DIVISION》 * |
| JOHN CANNY: "GaP: A Factor Model for Discrete Data", 《ACM》 * |
| MICHALIS TZORTZIS.ET: "Gamma radiation measurements and dose rates in commercially-used natural tiling rocks (granites)", 《UCY−PHY》 * |
| RAJESH ACHARYA.ET: "Software Application for Gamma Ray Computed Tomography Data Acquisition with Discrete Detectors", 《SOFTWARE APPLICATION FOR GAMMA RAY COMPUTED TOMOGRAPHY DATA》 * |
| S. THYKIER-NIELSEN: "Calculation Method for Gamma Dose Rates From Gaussian Puffs", 《RISO》 * |
| XUEJUN GU.ET: "GPU-based fast gamma index calculation", 《PHYS MED BIOL》 * |
| 张普忠等: "SCALE5.1程序系统中蒙特卡罗方法模块与离散纵坐标模块在乏燃料运输容器屏蔽计算中的比较分析", 《辐射防护》 * |
| 徐立鹏等: "FFT-BP神经网络模型对车载γ能谱辐射剂量率的预测分析", 《光谱学与光谱分析》 * |
| 李华等: "基于MCNP对γ射线吸收剂量累积因子的计算与研究", 《辐射防护》 * |
| 赛雪等: "Multiquadric散乱数据插值方法在γ辐射场可视化中的应用初探", 《核技术》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110221338A (en) * | 2019-05-17 | 2019-09-10 | 华南理工大学 | A method of reconstruct shield neighboring area radiation field |
| CN110221338B (en) * | 2019-05-17 | 2021-02-12 | 华南理工大学 | Method for reconstructing radiation field in peripheral area of shield |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109471147B (en) | 2022-10-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Flyer et al. | A guide to RBF-generated finite differences for nonlinear transport: shallow water simulations on a sphere | |
| O’Rourke et al. | An improved collision damping time for MP-PIC calculations of dense particle flows with applications to polydisperse sedimenting beds and colliding particle jets | |
| Desrayaud et al. | Benchmark solutions for natural convection flows in vertical channels submitted to different open boundary conditions | |
| US12118279B2 (en) | Lattice Boltzmann based solver for high speed flows | |
| Tang et al. | Convergence analysis for stochastic collocation methods to scalar hyperbolic equations with a random wave speed | |
| CN110489709A (en) | The method for numerical simulation of parsing Wall-function based on compressible flows | |
| CN108549753A (en) | A kind of radiation shield computational methods that Point- kernel integral method is coupled with Monte Carlo method | |
| CN101192307A (en) | Point cloud triangular topological relations construction method | |
| CN109471147A (en) | A kind of the gamma dose rate calculation method and system discrete based on variable weight | |
| Yildizeli et al. | Multi-objective optimization of multiple impinging jet system through genetic algorithm | |
| CN103405852A (en) | Rapid calculation method of stereotactic radiotherapy dosage field distribution | |
| CN102789547A (en) | Stay cable force calculation method taking actions of vibration attenuation damper into account | |
| Sharma et al. | hybrid: a CPU–GPU Monte Carlo method for modeling indirect x-ray detectors with columnar scintillators | |
| CN106354552A (en) | Parallel Computer Task Distribution Method and Device | |
| Liu et al. | Helical swimming in Stokes flow using a novel boundary-element method | |
| Puchalski et al. | Applications of four-body exponentially correlated functions | |
| Gotoh et al. | Spectrum of passive scalar at very high Schmidt number in turbulence | |
| Talukdar et al. | Investigation on optimization of the thermal performance for compressible laminar natural convection flow in open-ended vertical channel | |
| CN111650630B (en) | Arrangement method and system of radioactive sources | |
| Ghiyasvand et al. | A simple approximation algorithm for computing Arrow-Debreu prices | |
| Chopade et al. | Estimation of power of heaters in a radiant furnace for uniform thermal conditions on 3-D irregular shaped objects | |
| Kim et al. | Numerical simulation of air flow around multiple objects using the discrete vortex method | |
| Takeda et al. | A comparison of four viscous models for the discrete vortex method | |
| Yang et al. | MOCUM solutions and sensitivity study for C5G7 benchmark | |
| Wang et al. | Parallel 3D deterministic particle transport on Intel MIC architecture |
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 |