CN107195342A - A kind of method for measuring nuclear power plant's nucleic activity distribution - Google Patents
A kind of method for measuring nuclear power plant's nucleic activity distribution Download PDFInfo
- Publication number
- CN107195342A CN107195342A CN201710216323.8A CN201710216323A CN107195342A CN 107195342 A CN107195342 A CN 107195342A CN 201710216323 A CN201710216323 A CN 201710216323A CN 107195342 A CN107195342 A CN 107195342A
- Authority
- CN
- China
- Prior art keywords
- nucleic
- point
- impact
- currently assigned
- energy
- 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
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21C—NUCLEAR REACTORS
- G21C17/00—Monitoring; Testing ; Maintaining
- G21C17/017—Inspection or maintenance of pipe-lines or tubes in nuclear installations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/56—Particle system, point based geometry or rendering
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Abstract
The present invention provides a kind of method for measuring nuclear power plant's nucleic activity distribution, nucleic and measurement data are specified including obtaining detector, actual measurement measurement data is entered into preset mode as current continuous data, each single-row vector and its element number for specifying the actual measurement probability right value of nucleic to line up is obtained;Classification using each specified nucleic is mock up flat single-row vectorial element number corresponding with its as period, according to geometrical model circulating analog detector, obtain the measurement data simulated every time and enter preset mode processing, obtain each analog probability weighted value for specifying nucleic when it is modeled every time and the homotype matrix formed by analog probability weighted value;The single-row vectorial and its corresponding homotype matrix of nucleic is specified according to each, each activity distribution for specifying nucleic is calculated.Implement the present invention, existing radiation " focus " detection means can be overcome the shortcomings of, accurately measure out factory building radiation " focus " activity distribution situation.
Description
Technical field
The present invention relates to technical field of nuclear radiation measurement, more particularly to a kind of side for measuring nuclear power plant's nucleic activity distribution
Method.
Background technology
With the gradually aging of active service nuclear power plant of China, in reactor internal pipeline activation corrosion product ad-hoc location by
Year accumulation, forms gamma emitter, and this kind of radiation " focus " is being increased year by year so that during overhaul suffered by staff
Dose of radiation also increases year by year, and the problem of causing control collective's dose of radiation more protrudes.
Whether weaken ray, or the radiation field emerged by calculating " focus " to optimize operation by the method for shielding
Path, all be unable to do without the accurate measurement to " focus ".But, measurement radiation " focus " and common " find and lose radioactive source " are no
Together:First, " look for source " problem generally can be seen as point source, and nucleic species is single, is carried out into ripe gamma camera
Position as just caning be found that radioactive source, but factory building radiation " focus " generally all accumulates in pipeline or device interior, with one
Fixed geometry and size, it is impossible to be considered as point source, and nucleic composition is complicated, and gamma energy scope is wide, to detector
Energy resolution require high;Secondly, factory building radiation " focus " is related to three big key parameters:Position, energyWith
Activity, these three parameters are not three scalars, but the function that three independents variable are constituted, equivalent to obtain byx×y×z ×i×jThe matrix of individual unknown number composition.
In the prior art, generally using gamma ray spectrometer, cellular type γ cameras(Including orifice type and coding cellular type)And Compton
Scattering imaging type gamma camera radiates " focus " activity distribution situation to measure factory building.But, inventor has found, although gamma ray spectrometer
Power spectrum can be surveyed and energy and intensity is provided(I.e.i×jMatrix), but without position resolution, survey is all position radiation
The mixing resultant in source;Cellular type γ cameras can provide two-dimensional position and intensity(I.e.x×z×jMatrix), but on the one hand because generally adopting
With scintillator detector, energy resolution is typically inferior to 7%, it is impossible to effective district pyrene element, and survey be all nucleic mixing it is strong
Degree, the another aspect visual field is generally less than 60 °, just must respectively shot towards all multi-angles when measuring space in a big way, time-consuming
Arduously;Although Compton scatter imaging formula gamma camera can provide two-dimensional position, energy and intensity simultaneously(I.e.x×z×i×j
Matrix), and without collimater, visual angle is up to 360 °, but the result of imaging is only capable of presenting with the two dimensional form at longitude and latitude angle, shoots
Go out the transmitted intensity picture being superimposed on optical photograph, it is impossible to solve radiation " focus " by quantitative measurment source strength during material blocks
The problem of.
Therefore, need a kind of method of measurement factory building radiation " focus " activity distribution badly, existing radiation " focus " can be overcome
The deficiency of detection means, accurately measures out factory building radiation " focus " activity distribution situation.
The content of the invention
Technical problem to be solved of the embodiment of the present invention is that there is provided a kind of side for measuring nuclear power plant's nucleic activity distribution
Method, can overcome the shortcomings of existing radiation " focus " detection means, accurately measure out factory building radiation " focus " activity distribution situation,
Nuclear power plant's radiation " focus " is reflected while related to three-dimensional position, energy and intensity(Give three big key parametersx
×y×z×i×jMatrix).
In order to solve the above-mentioned technical problem, the embodiments of the invention provide a kind of side for measuring nuclear power plant's nucleic activity distribution
Method, methods described includes:
Obtain to be surveyed to specify after nuclear power plant's radioactive source to be measured by detector and switch to hit inside nucleic and gained correspondence gamma events
The measurement data of point, and enter calculating in preset mode using the actual measurement measurement data got as current continuous data,
And using the classification of each specified nucleic as unit of account, respectively obtain each actual measurement probability right value for specifying nucleic, enter
One step obtains each and specifies nucleic single-row vectorial by its correspondence constituted single-row vectorial of actual measurement probability right value and correspondence
Element number.
The geometrical model of space and the detector residing for radioactive source to be measured in the nuclear power plant is set up using illiteracy card program,
And single-row vectorial member is constituted by actual measurement probability right value by mock up flat and its correspondence using the classification of each specified nucleic
Plain number is period, blocks detector detection nuclear power plant radiation to be measured described in the geometrical model circulating analog built by described cover
Behind source, measurement data when gamma events obtained by simulating the detector every time switch to the point of impact is obtained respectively, and by each institute
State the analogue measurement data got and enter calculating in the preset mode respectively as current continuous data, obtain each and refer to
Determine analog probability weighted value of the nucleic when it is modeled every time, same specified nucleic is further modeled what is obtained every time
Analog probability weighted value is that column vector forms its corresponding homotype matrix;Wherein, the actual measurement measurement data and the simulation are surveyed
Measure the data three-dimensional position including the point of impact, hit time and energy deposition.
According to each obtained specify nucleic by its correspondence actual measurement probability right value constituted it is single-row vector and
The homotype matrix formed by its correspondence analog probability weighted value by column vector, passes through equationAx=b, calculate each and specify
The activity distribution of nucleic;Wherein,AFor each specify nucleic by its correspondence analog probability weighted value by column vector formed it is same
Type matrix,bThe single-row vector that nucleic is made up of its correspondence actual measurement probability right value is specified for each,xSpecified for each
The vector that to be solved activity value of the nucleic in each space lattice confined is lined up.
Wherein, the current continuous data enters the specific steps calculated in the preset mode and included:
It is determined that current continuous data enters the classification of corresponding currently assigned nucleic when calculating, and according to the class of currently assigned nucleic
Not, in default nucleic characteristic gamma ray energy scale, the gamma energy of currently assigned nucleic is obtained, and according to described
The gamma energy of the preset energy resolution ratio of detector and the currently assigned nucleic, determines currently assigned nucleic in energy
Energy window scope in spectrum;
According to the hit time of each point of impact in current continuous data, all hits in the current continuous data are clicked through
Row primary screening, is met the point of impact and its combination of the corresponding point of impact of default combination condition, and according to the current finger
Determine each point of impact in the gamma energy and its corresponding energy window scope and the current continuous data of nucleic
Energy is deposited, and the point of impact obtained respectively to primary screening carries out postsearch screening, obtains the point of impact needed for currently assigned nucleic
And its corresponding point of impact combination, and further according to the gamma energy of the currently assigned nucleic and its corresponding hit
Point combination, obtains the Compton scattering point and photoelectric absorption point of currently assigned nucleic;
Deposited according to the three-dimensional position of each point of impact in current continuous data and energy, determine own in currently assigned nucleic
Three-dimensional position and the energy deposition of Compton scattering point and photoelectric absorption point in point of impact combination, and with currently assigned nucleic
Any group of point of impact combination is unit, according to Compton scattering point and photoelectricity in each group of point of impact combination in currently assigned nucleic
The three-dimensional position of point and the gamma energy of energy deposition and currently assigned nucleic are absorbed, to any in currently assigned nucleic
Point of impact combination carries out circular conical surface virtually, obtains combining corresponding virtual circular cone with each point of impact in currently assigned nucleic
Face;Wherein, each described virtual circular conical surface is with Compton scattering point in correspondence point of impact combination in currently assigned nucleic
Three-dimensional position is summit, with Compton scattering point in correspondence point of impact combination in currently assigned nucleic to relative photoelectric absorption point
Vector be axle and using in currently assigned nucleic in the combination of the correspondence point of impact specific Compton scattering angle as semi-cone angle is built
Into;Specific Compton scattering angle is the point of impact group in its currently assigned nucleic of correspondence in each described virtual circular conical surface
The energy deposition of Compton scattering point and the gamma energy of currently assigned nucleic are determined in conjunction;
Discretization traversal processing is carried out to each described virtual circular conical surface obtained by currently assigned nucleic, and according to discrete
Change traversal processing result, in the cube grid that default three dimensions is divided, filtering out storage has currently assigned nucleic
The cube grid of interior virtual circular conical surface and its quantity of the virtual circular conical surface of correspondence storage, and further screened according to each
Cube grid deposits the quantity of virtual circular conical surface, obtains currently assigned nucleic and is surveyed corresponding to current continuous data for the actual measurement
When actual measurement probability right value or currently assigned nucleic when measuring data correspond to current continuous data for the analogue measurement data
Analog probability weighted value;Wherein, the default three dimensions is by building in the nuclear power plant residing for radioactive source to be measured
The three-dimensional cartesian coordinate system that space matches is realized, including multiple cube grids with identical specifiable lattice size.
Wherein, in the current continuous data of the basis each point of impact the hit time, to the current continuous data
In all points of impact carry out primary screenings, be met the point of impact and its combination of the corresponding point of impact of default combination condition, and
According in the gamma energy of the currently assigned nucleic and its corresponding energy window scope and the current continuous data
The energy deposition of each point of impact, the point of impact obtained respectively to primary screening carries out postsearch screening, obtains currently assigned core
The required point of impact of element and its combination of the corresponding point of impact, and further according to the gamma energy of the currently assigned nucleic
And its combination of the corresponding point of impact, specifically wrap the step of the Compton scattering point and photoelectric absorption point that obtain currently assigned nucleic
Include:
Ordered to any combination of two of the point of impact in the current continuous data, and according to each in the current continuous data
The hit time at midpoint, two points of impact of all any combination are screened, retains and to meet the default combobar
The point of impact of part and its combination of formation;
According to the energy deposition of each point of impact in the current continuous data, determine that the energy of the reservation point of impact sinks
Product, and calculate the cumulative sum of energy deposition of every group of point of impact in the reservation point of impact correspondence combination;
According to the energy of every group of point of impact in the energy window scope of the currently assigned nucleic and the reservation combination calculated
The cumulative sum of deposition, is screened to the reservation combination, is obtained currently assigned nucleic and is met the group that preset energy contrasts condition
Close and its required point of impact of correspondence;
According to the gamma energy of the currently assigned nucleic and its corresponding combination, the gamma of currently assigned nucleic is set in
When ray energy is more than or equal to predetermined threshold value, to deposit the big point of impact equal for energy in the corresponding any combination of currently assigned nucleic
It is photoelectric absorption point to deposit the small point of impact for Compton scattering point, energy;Or be set in the gamma of currently assigned nucleic and penetrate
It is light that energy, which deposits the big point of impact, when heat input is less than the predetermined threshold value, in the corresponding any combination of currently assigned nucleic
It is Compton scattering point that electric absorption point, energy, which deposit the small point of impact,.
Wherein, the default combination condition should be less than for the hit time interval of two points of impact among any reservation combination
Two points of impact are respectively provided with uniqueness among the detector response time, and any reservation combination.
Wherein, the preset energy contrast condition is any cumulative sum of energy deposition for retaining two points of impact in combination
It should be located in the range of the energy window of the currently assigned nucleic.
Wherein, specific Compton scattering angle passes through formula in each described described virtual circular conical surfaceTo realize;Wherein,θFor Compton scattering angle, E1To be any in currently assigned nucleic
The energy deposition of Compton scattering point, E in point of impact combination0For the gamma energy of currently assigned nucleic, mec2For photon matter
Amount, it is constant.
Wherein, it is described that the described virtual circular conical surface of each resulting in currently assigned nucleic is carried out at discretization traversal
Reason, and according to discretization traversal processing result, in the cube grid that default three dimensions is divided, filtering out storage has
The quantity of the cube grid of virtual circular conical surface and its virtual circular conical surface of correspondence storage in currently assigned nucleic, and further basis
Each screens the quantity that cube grid deposits virtual circular conical surface, obtains currently assigned nucleic corresponding to current continuous data
Actual measurement probability right value or currently assigned nucleic during for the actual measurement measurement data corresponding to current continuous data are the mould
Intend measurement data when analog probability weighted value the step of specifically include:
All cube grids in the default three dimensions are entered as 0;
Virtual circular conical surface execution first described to each resulting in currently assigned nucleic is recirculated after traversal, is continued to working as
Each described virtual circular conical surface obtained by preceding specified nucleic performs second and recirculated traversal, filters out double searching loop
Storage has in currently assigned nucleic the cube grid of virtual circular conical surface and counts each and screen and deposited in cube grid afterwards
Put the quantity of virtual circular conical surface;Wherein, described first traversal is recirculated using imaginary circles cone element length as amplitude variations, and set
Tentering degree variate-value is incrementally traveled through since 0 with specific amplitude increment, until set amplitude variations value is beyond described
Untill the border of default three dimensions;Described second recirculate traversal when with each amplitude variations value being respectively radius of circle institute it is right
The central angle answered is angle variables, and set angle variate-value incrementally traveled through since 0 ing with special angle increment, up to institute
Untill the angle variables value set is 2 π;Wherein, the specific amplitude increment is set to the 1/ of cube grid specifiable lattice size
10-1/3 one;The special angle increment be set to arc length correspondence cube grid specifiable lattice size 1/10-1/5 its
Central angle for the moment;
Each is screened into quantity that cube grid counts while being multiplied by a set-point, and by resulting product
Corresponding to imparting respectively, each screens cube grid;
When detecting current continuous data for the actual measurement measurement data, each by the imparting screens cube net
The product of lattice as currently assigned nucleic actual measurement probability right value;
When detecting current continuous data for the analogue measurement data, each by the imparting screens cube net
The product of lattice as currently assigned nucleic analog probability weighted value.
Wherein, the three-dimensional position X of the point of impact, Y-direction pass through in the actual measurement measurement data and the analogue measurement
The pixel sequence number of the detector determines, Z-direction be using nuclear electronics record the electronics drift duration of gamma events come
It is determined that.
Wherein, the actual measurement measurement data is by the way that the detector is placed in residing for radioactive source to be measured described in nuclear power plant
Outside space and in the range of 1 to 10 meter of the radioactive source to be measured, open and detect the radioactive source to be measured 1 to 120 minute
And obtain.
Wherein, the classification of the specified nucleic includes58Co、60Co、110mAg、124Sb、137Cs、51Cr、54Mn、59Fe、95Zr
、95Nb、134Cs、131I。
Implement the embodiment of the present invention, have the advantages that:
1st, in embodiments of the present invention, either after detector direct detection radioactive source to be measured gained actual measurement measurement data, also
It is the analogue measurement data that analog prober detects gained after radioactive source to be measured, the two includes the point of impact(I.e. radiation " focus ",
Nucleic)Three-dimensional position(x×y×z), energy deposition(i), and the two is involved in into preset mode calculating and draws nucleic
Activity(j)Distribution situation, so as to set up nuclear power plant's radiation " focus " while related to three-dimensional position, energy and intensity(Give
Three big key parameters are gone outx×y×z×i×jMatrix), overcome the deficiency of existing radiation " focus " detection means;
2nd, in embodiments of the present invention, the corresponding mould of analogue measurement data as obtained by detecting after radioactive source to be measured analog prober
Actual measurement of the homotype matrix that quasi-probability weighted value is formed to correct gained after detector direct detection radioactive source to be measured measures number
According to the single-row vector that is formed of actual measurement probability right value, the activity distribution situation of the nucleic drawn is more accurate, so as to improve
Computational accuracy.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
There is the accompanying drawing used required in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments of invention, for those of ordinary skill in the art, without having to pay creative labor, according to
These accompanying drawings obtain other accompanying drawings and still fall within scope of the invention.
Fig. 1 is a kind of flow chart for the method for measuring nuclear power plant's nucleic activity distribution provided in an embodiment of the present invention.
Embodiment
To make the object, technical solutions and advantages of the present invention clearer, the present invention is made into one below in conjunction with accompanying drawing
It is described in detail on step ground.
As shown in figure 1, in the embodiment of the present invention there is provided a kind of measurement nuclear power plant nucleic activity distribution method, institute
The method of stating includes:
Step S1, acquisition are surveyed gained after nuclear power plant's radioactive source to be measured by detector and specified in nucleic and gained correspondence gamma events
Portion switchs to the measurement data of the point of impact, and enters default mould using the actual measurement measurement data got as current continuous data
Calculated in formula, and using the classification of each specified nucleic as unit of account, respectively obtain each actual measurement probability for specifying nucleic
Weighted value, further obtains each single-row vectorial and correspondence list for specifying nucleic to be made up of its correspondence actual measurement probability right value
The element number of column vector;
Detailed process is, first before step S1, and the selection sensitive gamma detector of three-dimensional position treats measuring to survey nuclear power plant
Penetrate source, this kind of detector is able to record that each hit detector and produces the gamma events of energy deposition there is provided it in detection
The three-dimensional position of the point of impact, hit time and energy deposition inside device.In one embodiment, this kind of detector generally uses tellurium
Zinc cadmium(CdZnTe)Semiconductor array detector, can also be made up of other materials.
The three-dimensional position of this kind of detector record point of impact can use but be not limited to following methods:X, Y-direction position
Determined, i.e., realized by the way that detector sensitive area is carried out into pixelation, typical way by the pixel sequence number of detector
That anode is divided into latticed, common cathode, using anode sequence number be converted into X, Y-direction positional information;Z-direction is profit
The electronics drift duration of gamma events is recorded with nuclear electronics to determine, i.e., the electric float of each event is recorded using nuclear electronics
Duration is moved, the depth information of Z-direction is converted into, so allows for the three-dimensional of the point of impact of each gamma events inside detector
Positional information can be read and be used.
The sensitive gamma detector of three-dimensional position is placed in space residing for nuclear power plant's radioactive source to be measured, i.e., may such as be contained
Radiate " focus "(Radioactive source i.e. to be measured)Factory building in because " focus " is typically deposited on the positions such as pipeline or valve, and distance
In the range of 1 to 10 meter of radioactive source to be measured, open and detect radioactive source to be measured 1 to 120 minute, record each gamma events pair
Three-dimensional position, hit time and the energy deposition of the point of impact should be switched to.These information constitute gamma events data set, Ke Yi
The data that recorded are transmitted to computer at any time in detection process, total data can also be packed again after detection terminates and passed
It is defeated.As can be seen here, the measured data that computer is got includes three-dimensional position, hit time and the energy deposition of the point of impact.
Secondly, before step S1, in addition it is also necessary to set up the feature of each specified nucleic of actual measurement nuclear power plant radioactive source to be measured
Gamma energy table(The characteristic gamma ray energy scale of i.e. pre-set specified nucleic), wherein, specify the classification of nucleic
Including but not limited to58Co、60Co、110mAg、124Sb、137Cs、51Cr、54Mn、59Fe、95Zr、95Nb、134Cs、131I.Selecting each
When planting the characteristic gamma ray of specified nucleic as measurement object, following strategy is mainly followed:
(1)Specifying the energy of nucleic gamma ray will be located in the response range of detector(250-1500keV preferentially is selected, secondly
Select 50-250keV and 1500-2000keV);(2)Preferential, the ray energy for specifying the emission probability of nucleic gamma ray high
Larger is preferential;(3)According to the energy resolution of detector, avoid with falling into oblivion peak(511keV), backscattering peak(~180keV)Cross
In close energy;(4)For specifying nucleic to have the alternative nucleic of a plurality of gamma ray, 1 gamma-rays is generally selected as master
Object is surveyed, reselection 1-3 bar gamma rays are used as verification.
In one embodiment, according to above-mentioned strategy, the characteristic gamma ray of specified nucleic is selected, it is as shown in table 1 below:
Table 1
Nucleic | Main survey gamma energyE 0 | Gamma ray for verification |
58Co | 810.8keV | Nothing |
60Co | 1332.5keV | 1173.2keV |
110mAg | 884.7keV | 657.8keV、937.5keV、1384.3keV |
124Sb | 1691.0keV | 602.7keV |
137Cs | 661.7keV | Nothing |
51Cr | 320.1keV | Nothing |
54Mn | 834.8keV | Nothing |
59Fe | 1099.2keV | 1291.6keV |
95Zr | 724.2keV | 756.7keV |
95Nb | 765.8keV | Nothing |
134Cs | 795.84keV | 604.7keV |
131I | 364.48keV | Nothing |
Finally, each obtained before step S1 is specified into nucleic, actual measurement measurement data and default nucleic characteristic gamma
Radial energy scale is all imported in computer, and will actual measurement measurement data as continuous data by internal preset mode computation, with
Just obtain the actual measurement probability right value of each nucleic and its be correspondingly formed single-row vectorial element number.
In embodiments of the present invention, actual measurement measurement data enters the process tool that preset mode is calculated as current continuous data
Body is as follows:
Step S11, the classification for determining corresponding currently assigned nucleic when current continuous data enters calculating, and according to currently assigned
The classification of nucleic, in default nucleic characteristic gamma ray energy scale, obtains the gamma energy of currently assigned nucleic, and
According to the preset energy resolution ratio of the detector and the gamma energy of the currently assigned nucleic, determine currently assigned
Energy window scope of the nucleic in power spectrum;
Specifically, according to the classification of currently assigned nucleic, corresponding gamma energy is found, according to the preset energy of detector
Resolution ratio(Representative value is 1%-3%), converse energy window of the currently assigned nucleic in power spectrum in upper table 1.Energy window width can
To take 2.54 times of the all-round peak width FWHM of detector, other width can also be taken.Formula is as follows:Energy window bound=E 0±1.27·FWHM = E 0±1.27·R·E 0;
In one embodiment, the currently assigned nucleic into calculating is58Co, its characteristic gamma ray energyE 0In upper table 1
Find as 810.8keV, now detector is the energy resolution of the sensitive gamma detector of three-dimensional positionR=2%, converse to be measured
Nucleic58Energy window scopes of the Co in power spectrum is 790.2keV-831.4keV.
Step S12, the hit time according to each point of impact in current continuous data, in the current continuous data
All points of impact carry out primary screening, are met the point of impact and its combination of the corresponding point of impact of default combination condition, and root
According to every in the gamma energy and its corresponding energy window scope of the currently assigned nucleic and the current continuous data
The energy deposition of one point of impact, the point of impact obtained respectively to primary screening carries out postsearch screening, obtains currently assigned nucleic
The required point of impact and its combination of the corresponding point of impact, and further according to the gamma energy of the currently assigned nucleic and
Its corresponding point of impact combination, obtains the Compton scattering point and photoelectric absorption point of currently assigned nucleic;
Specifically, step S121, to any combination of two of the point of impact in the current continuous data, and according to the current meter
The hit time of each point of impact in data is measured, two points of impact of all any combination are screened, retains and meets
The point of impact of the default combination condition and its combination of formation;Wherein, the default combination condition is any reservation combination
In the hit time intervals of two points of impact should be less than the detector response time, and any retain two hits among combination
Point is respectively provided with uniqueness;
It is understood that step S121 is that " double point of impact events " screening is carried out to current continuous data, method is first to judge
The hit time of each point of impact, it is less than detector response time for the interval that any two point of impact hits the time, all
It is considered as and is triggered by same gammaphoton, thus establishes the association in time of each point of impact, it follows which
The individual point of impact is triggered by same gammaphoton.The event that one gammaphoton produces 2 points of impact just is referred to as " double hits
Point event ", every such event is retained, and remaining is rejected more or less than the event of 2 points of impact.
Step S122, the energy deposition according to each point of impact in the current continuous data, determine the reservation life
The energy deposition at midpoint, and calculate the cumulative sum of energy deposition of every group of point of impact in the reservation point of impact correspondence combination;
It is understood that step S122 is to two points of impact in " double point of impact events " is remained in step S121 combination
Energy deposition summation, so as to postsearch screening, obtain the point of impact that currently assigned nucleic needs.
Every group in step S123, the energy window scope according to the currently assigned nucleic and the reservation combination calculated
The cumulative sum of energy deposition of the point of impact, screens to the reservation combination, obtains currently assigned nucleic and meet preset energy
The combination of contrast condition and its required point of impact of correspondence;Wherein, the preset energy contrast condition retains in combination to be any
The cumulative sum of energy deposition of two points of impact should be located in the range of the energy window of the currently assigned nucleic;
It is understood that the cumulative sum of energy deposition that step S123 is used for judging to retain one by one two points of impact in combination is
It is no in the range of the energy window of currently assigned nucleic;If it is, continuing to remain as can imaging event;If not,
Then abandon.
Step S124, the gamma energy according to the currently assigned nucleic and its corresponding combination, are set in current
When specifying the gamma energy of nucleic to be more than or equal to predetermined threshold value, energy sinks in the corresponding any combination of currently assigned nucleic
The big point of impact of product is that to deposit the small point of impact be photoelectric absorption point for Compton scattering point, energy;Or it is set in current finger
When the gamma energy for determining nucleic is less than the predetermined threshold value, energy deposition is big in the corresponding any combination of currently assigned nucleic
The point of impact be that to deposit the small point of impact be Compton scattering point for photoelectric absorption point, energy.
It is understood that step S124 is used for judging if the gamma energy of currently assigned nucleicE 0>=preset threshold
Value(Such as 400keV), then energy in any combination in currently assigned nucleic is deposited into the larger point of impact and is labeled as Compton scattering
Point, energy deposits the less point of impact and is labeled as photoelectric absorption point, if conversely, then judging the gamma ray of currently assigned nucleic
EnergyE 0<Predetermined threshold value(Such as 400keV), energy in any combination in currently assigned nucleic will then be deposited to the larger point of impact
Labeled as photoelectric absorption point, energy deposits the less point of impact and is labeled as Compton scattering point.
Step S13, deposited according to the three-dimensional position of each point of impact in current continuous data and energy, it is determined that currently finger
Determine three-dimensional position and the energy deposition of Compton scattering point and photoelectric absorption point in all point of impact combinations in nucleic, and with current
It is unit to specify any group of point of impact combination in nucleic, according to Compton in each group of point of impact combination in currently assigned nucleic
Three-dimensional position and the energy deposition and the gamma energy of currently assigned nucleic of scattering point and photoelectric absorption point, to currently referring to
Determine any point of impact combination in nucleic and carry out circular conical surface virtually, obtain combining right with each point of impact in currently assigned nucleic
The virtual circular conical surface answered;Wherein, each described virtual circular conical surface is with health in correspondence point of impact combination in currently assigned nucleic
The three-dimensional position of general scattering point is summit, with Compton scattering point in correspondence point of impact combination in currently assigned nucleic to relative
The vector of photoelectric absorption point be axle and using in currently assigned nucleic in the combination of the correspondence point of impact specific Compton scattering angle as
Semi-cone angle is built-up;Specific Compton scattering angle is by its currently assigned core of correspondence in each described virtual circular conical surface
The energy deposition of Compton scattering point and the gamma energy of currently assigned nucleic are determined in point of impact combination in plain;
Specifically, determining the three-dimensional position of Compton scattering point and photoelectric absorption point in all point of impact combinations in currently assigned nucleic
Put and deposited with energy, and in units of any group of point of impact combination in currently assigned nucleic, with the three-dimensional of Compton scattering point
Position is summit, with the vector of Compton scattering point to relative photoelectric absorption point for axle, with specific Compton scattering angleθFor half
Cone angle, builds virtual circular conical surface in space.
Wherein, specific Compton scattering angleθPass through formulaTo realize;Its
In,θFor Compton scattering angle,E 1Energy for Compton scattering point in any point of impact combination in currently assigned nucleic is deposited,E 0
For the gamma energy of currently assigned nucleic,m e c 2For photon quality, it is constant.
At step S14, virtual circular conical surface progress discretization traversal described to each resulting in currently assigned nucleic
Reason, and according to discretization traversal processing result, in the cube grid that default three dimensions is divided, filtering out storage has
The quantity of the cube grid of virtual circular conical surface and its virtual circular conical surface of correspondence storage in currently assigned nucleic, and further basis
Each screens the quantity that cube grid deposits virtual circular conical surface, obtains currently assigned nucleic corresponding to current continuous data
Actual measurement probability right value or currently assigned nucleic during for the actual measurement measurement data corresponding to current continuous data are the mould
Intend analog probability weighted value during measurement data;Wherein, the default three dimensions is treated by building in the nuclear power plant
The three-dimensional cartesian coordinate system that matches of space residing for radioactive source is surveyed to realize, including it is multiple vertical with identical specifiable lattice size
Cube grid.
Specifically, all cube grids in the default three dimensions are entered as into 0;
Virtual circular conical surface execution first described to each resulting in currently assigned nucleic is recirculated after traversal, is continued to working as
Each described virtual circular conical surface obtained by preceding specified nucleic performs second and recirculated traversal, filters out double searching loop
Storage, which has in currently assigned nucleic the cube grid of virtual circular conical surface and counts each, afterwards screens quilt in cube grid
The traversed quantity of virtual circular conical surface;Wherein, described first traversal is recirculated using imaginary circles cone element length as amplitude variations,
And set amplitude variations value and incrementally traveled through with specific amplitude increment since 0, until set amplitude variations value exceeds
Untill the border of the default three dimensions;Described second recirculates traversal when with each amplitude variations value being respectively radius of circle
Corresponding central angle is angle variables, and set angle variate-value incrementally traveled through since 0 with special angle increment, directly
Untill set angle variables value is 2 π;Wherein, the specific amplitude increment is set to cube grid specifiable lattice size
1/10-1/3 one;The special angle increment is set to arc length correspondence in cube grid specifiable lattice size 1/10-1/5
Its central angle for the moment;
Each is screened into cube grid by the traversed quantity of virtual circular conical surface while be multiplied by a set-point, and will
Resulting product corresponds to imparting respectively, and each screens cube grid.Then, each by the imparting is screened
The product of cube grid as currently assigned nucleic actual measurement probability right value.
It should be noted that default three dimensions is matched by building space residing for radioactive source to be measured in nuclear power plant
Three-dimensional cartesian coordinate system realize, including multiple cube grids with identical specifiable lattice size, the specifiable lattice chi
Very little can be 1-10cm or other sizes.Size of mesh opening determines the space essence of final source strength distribution measuring results
Fineness, therefore the setting of size of mesh opening will be with the typical component in space to be measured(For example, the pipe of " focus " may be deposited in factory building
Road)Geometric scale be adapted.
It is understood that recirculated by above-mentioned two, each virtual circular conical surface translates into computer program can be with
The series of points of traversal.For the point that each is traversed, according to its three-dimensional position judgment, which cube grid it belongs to,
And added up a upper value for the corresponding probability right of the grid(I.e. virtual circular conical surface quantity is multiplied by a value), this value can consolidate
It is set to 1 or other numerical value or the variable relevant with traversal point position.
Survey after the determination of probability right value, although it reflects the probability that its corresponding grid contains the nucleic, this
Probable value is not yet equal to the nucleic activity that the grid contains, but the two has unique conversion relation, it is necessary to by visiting in the measurements
Efficiency calibration is surveyed to determine.Therefore, after step s 14, it is necessary to which the actual measurement probability right value of currently assigned nucleic is converted into
The form of three-dimensional matrice, is singly lined up in single-row vectorial form by element and counts element number, so as to follow-up
The realization of detection efficiency calibration.
Step S2, set up the several of space residing for radioactive source to be measured and the detector in the nuclear power plant using covering card program
What model, and using each specify the classification of nucleic for mock up flat and its correspondingly from actual measurement probability right value constitute it is single-row to
The element number of amount is period, and detecting nuclear power plant by detector described in the geometrical model circulating analog for covering card structure treats
Survey after radioactive source, measurement data when gamma events obtained by simulating the detector every time switch to the point of impact is obtained respectively, and will
The analogue measurement data got enter in the preset mode respectively as current continuous data every time calculates, and obtains every
Same specified nucleic, is further modeled by a kind of analog probability weighted value of specified nucleic when it is modeled every time every time
Obtained analog probability weighted value is that column vector forms its corresponding homotype matrix;Wherein, it is described to survey measurement data and described
Analogue measurement the data three-dimensional position including the point of impact, hit time and energy deposition;
Detailed process is to be converted from common using full energy peak counting unlike the gamma spectrometer of nucleic activity, some is empty
Between the corresponding virtual circular conical surface of Compton scatter imaging of grid other grids in space can all be contributed, such asnFirst vectorbIn
Each element be from owningnThe superposition of the contribution of individual source item grid, it is therefore desirable ton×nMatrix could reflectnIt is individual
Source item activity element withnThe correlation of individual probability space element, remembers thisn×nMatrix be efficiency calibration matrix.Efficiency calibration
The of matrixjShownIndividual element, representing successively ought onlyjWhen individual source item grid has unit activity nucleic, to alln
Contribution of the individual source item grid in probability space matrix.Therefore detection efficiency calibration seeks to all members efficiency calibration matrix
Element is decided.
Accordingly, it would be desirable to set up space residing for radioactive source to be measured and the detector in the nuclear power plant using card program is covered
Geometrical model, and constitute single-row by actual measurement probability right value by mock up flat and its correspondence using the classification of each specified nucleic
The element number of vector is period, is detected by geometrical model circulating analog detector after nuclear power plant's radioactive source to be measured, respectively
Gamma events obtained by obtaining each analog prober switch to the measurement data during point of impact, and by the analogue measurement got every time
Data enter in preset mode respectively as current continuous data to be calculated, and is obtained each and is specified nucleic when it is modeled every time
Analog probability weighted value, it is column vector that same specified nucleic is further modeled to obtained analog probability weighted value every time
Form its corresponding homotype matrix.
It is understood that during geometrical model circulating analog, each can be confined according in step S1 in fact
The cube network in nucleic corresponding three-dimensional space is specified, circulating analog detector is so only needed to the spatial cuboids grid confined
Detection efficient, obtain each and specify analog probability weighted value for circulating each time of nucleic, will further follow each time
The form for the analog probability weighted value defiled vector that ring is obtained, and by each of same specified nucleic after all circulations terminate
Individual single-row vector lines up homotype matrix.
As can be seen here, one cycle is often completed, specifies the efficiency calibration matrix of nucleic just to fill in a row, is circulated more than
It is fully completed, entirelyn×nEfficiency calibration matrix be just filled, that is, form homotype matrixA。
In embodiments of the present invention, in the cyclic process each time of detection efficiency calibration, the analogue measurement number of nucleic is specified
According to continuous data can be taken as to be calculated into preset mode, the process of its each cycle calculations in step S1 with surveying measurement data
The process for entering preset mode calculating as current continuous data is identical, and the result simply finally drawn is different, that is, works as inspection
When measuring current continuous data for analogue measurement data, each of imparting is screened to the product of cube grid as current
The analog probability weighted value of nucleic is specified, specific calculating process refers to step S11 to step S14, and this is no longer going to repeat them.
Step S3, specified according to each obtained nucleic by its correspondence survey that probability right value constituted it is single-row
Vector and the homotype matrix formed by its correspondence analog probability weighted value by column vector, pass through equationAx=b, calculate each
Plant the activity distribution of specified nucleic;Wherein,ANucleic is specified to be column vector institute shape by its correspondence analog probability weighted value for each
Into homotype matrix,bThe single-row vector that nucleic is made up of its correspondence actual measurement probability right value is specified for each,xTo be each
Plant the vector that to be solved activity value of the specified nucleic in each space lattice confined is lined up.
Detailed process is to determine detection efficient matrix, i.e., after homotype matrix A, it is known that probability space vectorbWith not
Relation between the source item activity vector x known, you can be concisely expressed as equation with matrix and vector multiplicationAx=b.Solve this homogeneous
System of linear equations obtains vectorx, that is, obtainnThe nucleic activity of individual source item grid.
The application scenarios of method to measuring nuclear power plant's nucleic activity distribution in the embodiment of the present invention are described further, with
In the factory building with three horizontal pipes, contain in measurement horizontal pipe58Exemplified by Co nucleic activity, it is described as follows:
Before computer disposal, a cadmium-zinc-teiluride semiconductor array of the selection with the sensitive read out function of three-dimensional position is detected
Device, its number of pixels is 11 × 11, and detector is positioned in front of horizontal pipe on ground, is with horizontal pipe minimum distance
4m, start detection 10min, record actual measurement measurement data, i.e., three-dimensional position, hit time and the energy of each point of impact sink
Product, passes through total data on network cable transmission to computer after end to be detected;
Computer is got after data, actual measurement measurement data in carry out " double point of impact events " screening, leave 68433 it is " double
Point of impact event ", determines nucleic58Co characteristic gamma ray is 810.8keV, according to the sensitive gamma detector of three-dimensional position
Energy resolutionR=2%, converse nucleic58Energy window scopes of the Co in power spectrum is 790.2keV-831.4keV;
To each group " double point of impact events "(Retain combination)In two points of impact energy deposition summation, judge energy deposition
Whether sum belongs in the range of 790.2keV-831.4keV energy window;If belonging to, continue to retain, and be included into nucleic58Co
Combination in;If being not belonging to, give up;
To nucleic58Co any combination, because of gamma energyE=810.8keV>400keV, therefore by nucleic58All combinations of Co
In among two points of impact energy deposit the larger point of impact and be labeled as Compton scattering point, energy deposits less point of impact mark
It is designated as photoelectric absorption point;
To nucleic58Co each combination is unit, with the vector from photoelectric absorption point to Compton scattering point for reference direction,
And remember nucleic58The energy of Compton scattering point is deposited as in each combinations of CoE 1, converse nucleic58Compton in each combinations of Co
Angle of scatteringθ;Now, using Compton scattering point coordinate in space as summit, using reference direction as axle, withθFor semi-cone angle,
Virtual circular conical surface is constituted in space, then each combination corresponds to virtual circular conical surface as one;
The factory building for three horizontal pipes sets up three-dimensional cartesian coordinate system in a computer, and factory building size is 5m × 5m × 5m, and is drawn
Discrete cube grid, size of mesh opening is set as 4cm, and storing a probability right in a computer for each cube grid becomes
Amount, initial value is 0;
To all virtual circular conical surfaces, discretization traversal is carried out according to following double round-robin method one by one:First recirculates:With
Bus length is amplitude variations, is incrementally traveled through from 0 to 650cm, amplitude delta is set as 1cm;Second recirculates:For the first weight
Circulate with the girth of bus length correspondence circle, and with the central angle of the circleφFor angle variables, incrementally traveled through from 0 to 2 π, angle
Increment is set as;Each point that above-mentioned double searching loop is arrived, first judges which it belongs to empty according to its coordinate
Between grid, then added up a upper value for the corresponding probability right of the grid, this value is set as;
Complete nucleic58After the traversal of all virtual circular conical surfaces of Co, nucleic is obtained58Co actual measurement probability right value, so that one need to be entered
Step carries out detection efficiency calibration;
Now, space lattice is compareed with the 3-D geometric model of factory building, the part of source item nucleic is there may be in factory building
It is three horizontal pipes, the grid that they are accommodated has 360, and these grids are referred to as source item grid, remembers its numbern=360, that
Each source item grid Radionuclide58Co activity is all unknown number to be solved.According to order from high to low be source item grid according to
It is secondary be numbered with 1 to 360 sequence number, the unknown number vector of activity will be designated asx, it has 360 elements.At the same time, the probability obtained
The corresponding actual measurement probability right value of source item grid sequentially number is filtered out in the array of space, they are formed a line and is designated as vectorb, it
Also there are 360 elements;
Remember efficiency calibration matrix(That is homotype matrix)ForA, it is 360 × 360 square formation, and is covered by the PARTICLE TRANSPORT FROM of computer
Special Monte Carlo Simulation of Ions Inside program is carried out:360 × 360 array is set up with storage matrixA;With cover card program Geant4 set up factory building,
Water, the geometrical models of detector in three horizontal pipes, pipe, are only covering the of snap gauge typejSetting unit activity in individual source item grid
Nucleic, produce gamma ray according to emission probability, then carry out covering all events of card analog prober record, newly-built is general
Rate space carries out probability cone traversal, the probability right of all source item grids in probability space is read according to numbering, and write according to numbering
MatrixAThej1st to the 360th element of row.After the completion of circulating above, whole efficiency calibration matrixAJust it is filled.
Solve system of homogeneous linear equationsAx=b, obtain vectorx, that is, obtained 360 source item grids58Co nucleic activity, from
And the nucleic in three horizontal pipes58Co activity is distributed with regard to this determination.
Implement the embodiment of the present invention, have the advantages that:
1st, in embodiments of the present invention, either after detector direct detection radioactive source to be measured gained actual measurement measurement data, also
It is the analogue measurement data that analog prober detects gained after radioactive source to be measured, the two includes the point of impact(I.e. radiation " focus ",
Nucleic)Three-dimensional position(x×y×z), energy deposition(i), and the two is involved in into preset mode calculating and draws nucleic
Activity(j)Distribution situation, so as to set up nuclear power plant's radiation " focus " while related to three-dimensional position, energy and intensity(Give
Three big key parameters are gone outx×y×z×i×jMatrix), overcome the deficiency of existing radiation " focus " detection means;
2nd, in embodiments of the present invention, the corresponding mould of analogue measurement data as obtained by detecting after radioactive source to be measured analog prober
Actual measurement of the homotype matrix that quasi-probability weighted value is formed to correct gained after detector direct detection radioactive source to be measured measures number
According to the single-row vector that is formed of actual measurement probability right value, the activity distribution situation of the nucleic drawn is more accurate, so as to improve
Computational accuracy.
Can be with one of ordinary skill in the art will appreciate that realizing that all or part of step in above-described embodiment method is
The hardware of correlation is instructed to complete by program, described program can be stored in a computer read/write memory medium,
Described storage medium, such as ROM/RAM, disk, CD.
Above disclosure is only preferred embodiment of present invention, can not limit the right model of the present invention with this certainly
Enclose, therefore the equivalent variations made according to the claims in the present invention, still belong to the scope that the present invention is covered.
Claims (10)
1. a kind of method for measuring nuclear power plant's nucleic activity distribution, it is characterised in that methods described includes:
Obtain to be surveyed to specify after nuclear power plant's radioactive source to be measured by detector and switch to hit inside nucleic and gained correspondence gamma events
The measurement data of point, and enter calculating in preset mode using the actual measurement measurement data got as current continuous data,
And using the classification of each specified nucleic as unit of account, respectively obtain each actual measurement probability right value for specifying nucleic, enter
One step obtains each and specifies nucleic single-row vectorial by its correspondence constituted single-row vectorial of actual measurement probability right value and correspondence
Element number;
The geometrical model of space and the detector residing for radioactive source to be measured in the nuclear power plant is set up using card program is covered, and with
Each specifies the classification of nucleic to constitute single-row vectorial element by actual measurement probability right value by mock up flat and its correspondence
Number is period, and nuclear power plant's radioactive source to be measured is detected by detector described in the geometrical model circulating analog for covering card structure
Afterwards, measurement data when gamma events obtained by simulating the detector every time switch to the point of impact is obtained respectively, and will be described every time
The analogue measurement data got enter in the preset mode respectively as current continuous data to be calculated, and is obtained each and is specified
Same specified nucleic, is further modeled obtained mould by analog probability weighted value of the nucleic when it is modeled every time every time
Quasi-probability weighted value is that column vector forms its corresponding homotype matrix;Wherein, the actual measurement measurement data and the analogue measurement
The data three-dimensional position including the point of impact, hit time and energy deposition;
Nucleic is specified by its single-row vector that correspondingly actual measurement probability right value is constituted and by it according to each obtained
The homotype matrix that correspondence analog probability weighted value is formed by column vector, passes through equationAx=b, calculate each and specify nucleic
Activity distribution;Wherein,AThe homotype square that nucleic is formed by its correspondence analog probability weighted value by column vector is specified for each
Battle array,bThe single-row vector that nucleic is made up of its correspondence actual measurement probability right value is specified for each,xNucleic is specified for each
The vector that activity value to be solved in each space lattice confined is lined up.
2. the method as described in claim 1, it is characterised in that the current continuous data, which enters in the preset mode, to be calculated
Specific steps include:
It is determined that current continuous data enters the classification of corresponding currently assigned nucleic when calculating, and according to the class of currently assigned nucleic
Not, in default nucleic characteristic gamma ray energy scale, the gamma energy of currently assigned nucleic is obtained, and according to described
The gamma energy of the preset energy resolution ratio of detector and the currently assigned nucleic, determines currently assigned nucleic in energy
Energy window scope in spectrum;
According to the hit time of each point of impact in current continuous data, all hits in the current continuous data are clicked through
Row primary screening, is met the point of impact and its combination of the corresponding point of impact of default combination condition, and according to the current finger
Determine each point of impact in the gamma energy and its corresponding energy window scope and the current continuous data of nucleic
Energy is deposited, and the point of impact obtained respectively to primary screening carries out postsearch screening, obtains the point of impact needed for currently assigned nucleic
And its corresponding point of impact combination, and further according to the gamma energy of the currently assigned nucleic and its corresponding hit
Point combination, obtains the Compton scattering point and photoelectric absorption point of currently assigned nucleic;
Deposited according to the three-dimensional position of each point of impact in current continuous data and energy, determine own in currently assigned nucleic
Three-dimensional position and the energy deposition of Compton scattering point and photoelectric absorption point in point of impact combination, and with currently assigned nucleic
Any group of point of impact combination is unit, according to Compton scattering point and photoelectricity in each group of point of impact combination in currently assigned nucleic
The three-dimensional position of point and the gamma energy of energy deposition and currently assigned nucleic are absorbed, to any in currently assigned nucleic
Point of impact combination carries out circular conical surface virtually, obtains combining corresponding virtual circular cone with each point of impact in currently assigned nucleic
Face;Wherein, each described virtual circular conical surface is with Compton scattering point in correspondence point of impact combination in currently assigned nucleic
Three-dimensional position is summit, with Compton scattering point in correspondence point of impact combination in currently assigned nucleic to relative photoelectric absorption point
Vector be axle and using in currently assigned nucleic in the combination of the correspondence point of impact specific Compton scattering angle as semi-cone angle is built
Into;Specific Compton scattering angle is the point of impact group in its currently assigned nucleic of correspondence in each described virtual circular conical surface
The energy deposition of Compton scattering point and the gamma energy of currently assigned nucleic are determined in conjunction;
Discretization traversal processing is carried out to each described virtual circular conical surface obtained by currently assigned nucleic, and according to discrete
Change traversal processing result, in the cube grid that default three dimensions is divided, filtering out storage has currently assigned nucleic
The cube grid of interior virtual circular conical surface and its quantity of the virtual circular conical surface of correspondence storage, and further screened according to each
Cube grid deposits the quantity of virtual circular conical surface, obtains currently assigned nucleic and is surveyed corresponding to current continuous data for the actual measurement
When actual measurement probability right value or currently assigned nucleic when measuring data correspond to current continuous data for the analogue measurement data
Analog probability weighted value;Wherein, the default three dimensions is by building in the nuclear power plant residing for radioactive source to be measured
The three-dimensional cartesian coordinate system that space matches is realized, including multiple cube grids with identical specifiable lattice size.
3. method as claimed in claim 2, it is characterised in that the life of each point of impact in the current continuous data of basis
The middle time, primary screening is carried out to all points of impact in the current continuous data, the hit of default combination condition is met
Point and its corresponding point of impact combination, and gamma energy and its corresponding energy window model according to the currently assigned nucleic
Enclose and the current continuous data in each point of impact energy deposition, the point of impact that is obtained respectively to primary screening carries out
Postsearch screening, obtains the point of impact needed for currently assigned nucleic and its combination of the corresponding point of impact, and further work as according to described
The gamma energy of preceding specified nucleic and its combination of the corresponding point of impact, obtain currently assigned nucleic Compton scattering point and
The step of photoelectric absorption point, specifically includes:
Ordered to any combination of two of the point of impact in the current continuous data, and according to each in the current continuous data
The hit time at midpoint, two points of impact of all any combination are screened, retains and to meet the default combobar
The point of impact of part and its combination of formation;
According to the energy deposition of each point of impact in the current continuous data, determine that the energy of the reservation point of impact sinks
Product, and calculate the cumulative sum of energy deposition of every group of point of impact in the reservation point of impact correspondence combination;
According to the energy of every group of point of impact in the energy window scope of the currently assigned nucleic and the reservation combination calculated
The cumulative sum of deposition, is screened to the reservation combination, is obtained currently assigned nucleic and is met the group that preset energy contrasts condition
Close and its required point of impact of correspondence;
According to the gamma energy of the currently assigned nucleic and its corresponding combination, the gamma of currently assigned nucleic is set in
When ray energy is more than or equal to predetermined threshold value, to deposit the big point of impact equal for energy in the corresponding any combination of currently assigned nucleic
It is photoelectric absorption point to deposit the small point of impact for Compton scattering point, energy;Or be set in the gamma of currently assigned nucleic and penetrate
It is light that energy, which deposits the big point of impact, when heat input is less than the predetermined threshold value, in the corresponding any combination of currently assigned nucleic
It is Compton scattering point that electric absorption point, energy, which deposit the small point of impact,.
4. method as claimed in claim 3, it is characterised in that the default combination condition is two among any reservation combination
Two points of impact have among the hit time interval of the point of impact should be less than the detector response time, and any reservation combination
There is uniqueness.
5. method as claimed in claim 3, it is characterised in that the preset energy contrast condition retains two in combination to be any
The cumulative sum of energy deposition of the individual point of impact should be located in the range of the energy window of the currently assigned nucleic.
6. method as claimed in claim 2, it is characterised in that specific Compton in each described described virtual circular conical surface
Angle of scattering passes through formulaTo realize;Wherein,θFor Compton scattering angle, E1For
The energy of Compton scattering point is deposited in any point of impact combination in currently assigned nucleic, E0Gamma for currently assigned nucleic is penetrated
Heat input, mec2For photon quality, it is constant.
7. method as claimed in claim 2, it is characterised in that described to described in each resulting in currently assigned nucleic
Virtual circular conical surface carries out discretization traversal processing, and according to discretization traversal processing result, is divided in default three dimensions
Cube grid in, filtering out storage has in currently assigned nucleic the cube grid of virtual circular conical surface and its correspondence storage void
Intend the quantity of circular conical surface, and the quantity that cube grid deposits virtual circular conical surface is further screened according to each, worked as
Actual measurement probability right value or currently assigned core when preceding specified nucleic corresponds to current continuous data for the actual measurement measurement data
The step of element corresponds to analog probability weighted value when current continuous data is the analogue measurement data specifically includes:
All cube grids in the default three dimensions are entered as 0;
Virtual circular conical surface execution first described to each resulting in currently assigned nucleic is recirculated after traversal, is continued to working as
Each described virtual circular conical surface obtained by preceding specified nucleic performs second and recirculated traversal, filters out double searching loop
Storage, which has in currently assigned nucleic the cube grid of virtual circular conical surface and counts each, afterwards screens quilt in cube grid
The traversed quantity of virtual circular conical surface;Wherein, described first traversal is recirculated using imaginary circles cone element length as amplitude variations,
And set amplitude variations value and incrementally traveled through with specific amplitude increment since 0, until set amplitude variations value exceeds
Untill the border of the default three dimensions;Described second recirculates traversal when with each amplitude variations value being respectively radius of circle
Corresponding central angle is angle variables, and set angle variate-value incrementally traveled through since 0 with special angle increment, directly
Untill set angle variables value is 2 π;Wherein, the specific amplitude increment is set to cube grid specifiable lattice size
1/10-1/3 one;The special angle increment is set to arc length correspondence in cube grid specifiable lattice size 1/10-1/5
Its central angle for the moment;
Each is screened into cube grid by the traversed quantity of virtual circular conical surface while be multiplied by a set-point, and will
Resulting product corresponds to imparting respectively, and each screens cube grid;
When detecting current continuous data for the actual measurement measurement data, each by the imparting screens cube net
The product of lattice as currently assigned nucleic actual measurement probability right value;
When detecting current continuous data for the analogue measurement data, each by the imparting screens cube net
The product of lattice as currently assigned nucleic analog probability weighted value.
8. the method as any one of claim 1-7, it is characterised in that the actual measurement measurement data and the simulation are surveyed
The three-dimensional position X of the point of impact, Y-direction are the pixel sequence numbers by the detector to determine in amount, and Z-direction is to utilize core
The electronics of electronics record gamma events drifts about duration to determine.
9. the method as any one of claim 1-7, it is characterised in that the actual measurement measurement data is by by described in
Detector is placed in residing for radioactive source to be measured described in nuclear power plant outside space and apart from the scope of 1 to 10 meter of the radioactive source to be measured
It is interior, open and detect the radioactive source to be measured 1 to 120 minute and obtain.
10. the method as any one of claim 1-7, it is characterised in that the classification of the specified nucleic includes58Co、60Co、110mAg、124Sb、137Cs、51Cr、54Mn、59Fe、95Zr、95Nb、134Cs、131I。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710216323.8A CN107195342B (en) | 2017-04-05 | 2017-04-05 | A method of measurement nuclear power plant's nucleic activity distribution |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710216323.8A CN107195342B (en) | 2017-04-05 | 2017-04-05 | A method of measurement nuclear power plant's nucleic activity distribution |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107195342A true CN107195342A (en) | 2017-09-22 |
CN107195342B CN107195342B (en) | 2019-07-30 |
Family
ID=59871088
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710216323.8A Active CN107195342B (en) | 2017-04-05 | 2017-04-05 | A method of measurement nuclear power plant's nucleic activity distribution |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107195342B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108319160A (en) * | 2018-02-01 | 2018-07-24 | 中国核电工程有限公司 | Nuclear power plant master control room simulator system based on virtual reality technology |
CN108535768A (en) * | 2018-03-13 | 2018-09-14 | 江苏超敏科技有限公司 | A kind of gamma camera based on double detector technology |
CN109814147A (en) * | 2019-01-01 | 2019-05-28 | 中国人民解放军63653部队 | Use the method for gamma activity Nuclear analysis around detector plus decaying bulk measurement aperture |
WO2021139020A1 (en) * | 2020-01-12 | 2021-07-15 | 兰州大学 | Multi-parameter detection system and method for intra-drill-hole nuclear facility decommissioning waste radioactivity |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103853929A (en) * | 2014-03-17 | 2014-06-11 | 东华理工大学 | Low-resolution gamma energy spectrum inversion analysis process and method based on Monte Carlo response matrix |
CN104020507A (en) * | 2014-06-13 | 2014-09-03 | 西北核技术研究所 | Radioactivity area determining method based on virtual point detector principle |
CN204009085U (en) * | 2014-07-15 | 2014-12-10 | 中广核工程有限公司 | A kind of nuclear power plant hydrogeneous waste gas radioactivity continuous monitoring device |
CN105355256A (en) * | 2015-09-30 | 2016-02-24 | 湖北省农业科学院农产品加工与核农技术研究所 | Method for arranging radioactive sources |
CN106199676A (en) * | 2015-04-30 | 2016-12-07 | 北京中智核安科技有限公司 | A kind of gamma detector passive efficiency scale new method |
-
2017
- 2017-04-05 CN CN201710216323.8A patent/CN107195342B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103853929A (en) * | 2014-03-17 | 2014-06-11 | 东华理工大学 | Low-resolution gamma energy spectrum inversion analysis process and method based on Monte Carlo response matrix |
CN104020507A (en) * | 2014-06-13 | 2014-09-03 | 西北核技术研究所 | Radioactivity area determining method based on virtual point detector principle |
CN204009085U (en) * | 2014-07-15 | 2014-12-10 | 中广核工程有限公司 | A kind of nuclear power plant hydrogeneous waste gas radioactivity continuous monitoring device |
CN106199676A (en) * | 2015-04-30 | 2016-12-07 | 北京中智核安科技有限公司 | A kind of gamma detector passive efficiency scale new method |
CN105355256A (en) * | 2015-09-30 | 2016-02-24 | 湖北省农业科学院农产品加工与核农技术研究所 | Method for arranging radioactive sources |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108319160A (en) * | 2018-02-01 | 2018-07-24 | 中国核电工程有限公司 | Nuclear power plant master control room simulator system based on virtual reality technology |
CN108319160B (en) * | 2018-02-01 | 2021-04-16 | 中国核电工程有限公司 | Nuclear power station main control room simulator system based on virtual reality technology |
CN108535768A (en) * | 2018-03-13 | 2018-09-14 | 江苏超敏科技有限公司 | A kind of gamma camera based on double detector technology |
CN108535768B (en) * | 2018-03-13 | 2021-08-20 | 江苏超敏科技有限公司 | Gamma camera based on double-detector technology |
CN109814147A (en) * | 2019-01-01 | 2019-05-28 | 中国人民解放军63653部队 | Use the method for gamma activity Nuclear analysis around detector plus decaying bulk measurement aperture |
CN109814147B (en) * | 2019-01-01 | 2021-01-29 | 中国人民解放军63653部队 | Method for measuring gamma radionuclide distribution around small hole by using detector and attenuator |
WO2021139020A1 (en) * | 2020-01-12 | 2021-07-15 | 兰州大学 | Multi-parameter detection system and method for intra-drill-hole nuclear facility decommissioning waste radioactivity |
Also Published As
Publication number | Publication date |
---|---|
CN107195342B (en) | 2019-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107195342B (en) | A method of measurement nuclear power plant's nucleic activity distribution | |
US9958561B2 (en) | Apparatus and method for the determination of one or more neutron source characteristics | |
US10281595B2 (en) | Method and apparatus for distinguishing radionuclide by using plastic scintillator | |
CN103901052B (en) | A kind of SGS and TGS combined measurement device and collimating apparatus optimization method | |
CN103645491B (en) | Radioactive source localization method, Apparatus and system | |
US7279676B2 (en) | Position sensitive radiation spectrometer | |
CN102262237A (en) | Photonic radiation detection device, and methods for dimensioning and operating such device | |
CA3175821A1 (en) | Gamma ray detection system and calibration method thereof | |
CN107688194A (en) | Data processing method and data processing equipment | |
CN108287357B (en) | Method for obtaining source peak detection efficiency of cylinder detector to point source | |
Lévai et al. | Feasibility of gamma emission tomography for partial defect verification of spent LWR fuel assemblies: Summary report on simulation and experimental studies including design options and cost-benefit analysis: Task JNT A1201 of the Support Programmes of Finland (FINSP), Hungary (HUNSP) and Sweden (SWESP) to the IAEA Safeguards | |
Poitrasson-Rivière et al. | Monte Carlo investigation of a high-efficiency, two-plane Compton camera for long-range localization of radioactive materials | |
Olson et al. | Passive detection of small low-emission sources: two-dimensional numerical case studies | |
CN111512184B (en) | Method for determining the amount of radioisotope | |
Jacobsson et al. | Outcomes of the JNT 1955 phase I viability study of gamma emission tomography for spent fuel verification | |
Brabants et al. | Minimal detection time for localization of radioactive hot spots in low and elevated background environments using a czt gamma-ray spectrometer | |
Lee et al. | Optimization of the SPECT systems based on a CdTe pixelated semiconductor detector using novel parallel-hole collimators | |
Virta et al. | Verifying spent nuclear fuel with Passive Gamma Emission Tomography prior to disposal in a geological repository in Finland | |
Stanley et al. | See inside: The development of a cosmic ray muon imaging system to aid the clean up of the UK’s nuclear waste legacy | |
Brown et al. | Thermal neutron source location using a 3-D position-sensitive CdZnTe detector array | |
CN116577819B (en) | Multi-head Compton detection method and system | |
Swift et al. | Attributes from NMIS time coincidence, fast-neutron imaging, fission mapping, and gamma-ray spectrometry data | |
Stockman et al. | Facility scale in-situ source localization and assay via a sparse 3He neutron detector array: enhancing nuclear material control and accounting in nuclear fuel cycle facilities | |
Gagliardi et al. | Novel applications of state-of-the-art gamma-ray imaging technique: from nuclear decommissioning and radioprotection to radiological characterization and safeguards | |
Goldston et al. | Autonomous mobile directionally and spectrally sensitive neutron detectors |
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 |