Air kerma conventional true value determining method

Abstract

An air kerma conventional true value determining method is provided, which addresses the problem of on-site and in-situ verification or calibration of radiation protection with existing standard reference radiation, which is large in spatial volume and unable or difficult to be moved. The method includes establishing a minitype reference radiation, selecting a proper radiation source and source intensity for providing incident rays for a shielding box, selecting a plurality of gamma ray dosimeters as samples for training a prediction model to obtain the prediction model of the air kerma conventional true value of a point of test, putting a probe of a dosimeter being verified at the point of test, recording scattering gamma spectrum measured by a gamma spectrometer, with the spectrum applied as input to the prediction model to obtain the air kerma conventional true value. The results are accurate and the reference radiation is small in size.

Claims

1. An air kerma conventional true value determining method, comprising the following steps: Step 1, establishing a minitype reference radiation, which comprises a shielding box comprising a side length not more than 1.5 meters, the shielding box being positioned horizontally and an incident hole being provided on a side thereof for incidence of incident rays, a point of test being arranged in a direction of the incident rays in the shielding box, the shielding box being further provided with a test hole on an upper surface through which a probe of an dosimeter being verified can be put into the shielding box, a reference point on the probe coinciding with the point of test, a monitor point being also arranged in the shielding box, the shielding box being segmented into two parts by one plane perpendicular to a connecting line of the incident hole and the point of test, the monitor point being located at a part adjacent to the incident hole in the shielding box and at a position not directly irradiated by the incident rays, a gamma spectrometer being arranged in the shielding box, a reference point on a probe thereof coinciding with the monitor point and the probe being fixed in the shielding box; Step 2, selecting a proper radiation source and source intensity to provide incident rays for the shielding box; Step 3, selecting a plurality of gamma ray dosimeters as samples for straining a prediction model to obtain the prediction model of the air kerma conventional true value at the point of test; and Step 4, putting the probe of the dosimeter being verified at the point of test, measuring a scattering gamma spectrum by the gamma spectrometer, with the prediction model and the gamma spectrum as input, the air kerma conventional true value at the point of the test being obtained.

2. The air kerma conventional true value determining method of claim 1, wherein Step 3 comprises the following specific steps: Step 31, selecting a plurality of gamma ray dosimeters as samples for straining a prediction model; Step 32, measuring the air kerma conventional true value at the point of test when no sample dosimeter is put, putting a reference point of a probe of an sample dosimeter on the point of test, measuring the air kerma conventional true value at the point of test by adopting an instrument transfer method, and measuring the gamma spectrum of the monitor point via the gamma spectrometer; step 33, extracting a dose feature value by adopting a principal component analysis method according to the gamma energy spectra; and Step 34, training a prediction model of the air kerma conventional true value at the point of test by adopting a support vector machine regression method.

3. The air kerma conventional true value measuring method of claim 2, wherein step 32 comprises the following specific steps: Step 32A, putting a standard graphite cavity ionization chamber at the point of test, and measuring the air kerma conventional true value K.sub.j at the point of test when the intensity of an incident ray beam is V.sub.j; Step 32B, putting the reference point of the probe of the i.sup.th sample dosimeter at the point of test, setting the intensity of the incident ray beam as V.sub.j, recording the indicate value of the sample dosimeter R.sub.ij and measuring the gamma spectrum S.sub.ij of the monitor point at the moment via the gamma spectrometer; Step 32C, putting the sample dosimeter in the standard reference radiation to search a point having the indicate value equal to R.sub.ij, the corresponding air kerma conventional true value of the point being the air kerma conventional true value K.sub.ij at the point of test; Step 32D, sequentially putting the reference points of the probes of the x sample dosimeters in the point of test, and repeating steps 32A to 32C under y source intensity conditions to obtain xy sets of K.sub.ij, S.sub.ij and K.sub.j data for straining a model of the function relationship K.sub.ij=f.sub.1(S.sub.ij,K.sub.j).

4. The air kerma conventional true value measuring method of claim 3, wherein step 33 comprises the following steps: Step 33A, dispersing each acquired S.sub.ij according to a proper energy interval E to obtain a counting rate .sub.ijn array corresponding to energies of the scattering gamma rays, and constructing n-dimensional vectors a.sub.ij of the counting rates using the energies of the scattering gamma rays as a research object; Step 33B, constructing a scattering gamma spectrum counting rate data matrix sample .sub.(xy)n via experiments of the probes of the x sample dosimeters under the y source intensity conditions in step 32D; Step 33C, extracting feature components of the n-dimensional vectors a.sub.ij by adopting a principal component analysis method to obtain principal component vectors .sub.ij=T.sub.nm.sup.T.Math.a.sub.ij of the n-dimensional vectors a.sub.ij, mn, T.sub.nm.sup.T being a transposition of T.sub.nm, and T.sub.nm referring to obtaining a covariance matrix .sub.nn from .sub.(xy)n; and solving a score matrix composed of m first feature vectors of the covariance matrix .sub.nn; and Step 33D, obtaining a function relationship .sub.ij=f.sub.2(S.sub.ij) between .sub.ij and S.sub.ij, thus obtaining the model K.sub.ij=f.sub.3(.sub.ij,K.sub.j)=f.sub.3[f.sub.2(S.sub.ij),K.sub.j].

5. The air kerma conventional true value determining method of claim 4, wherein in step 33A, the proper energy interval E refers to: E=1500/(1282.sup.z) keV, 0z4, z being an integer.

6. The air kerma conventional true value determining method of claim 4, wherein step 33C comprises the following specific steps: Step 33C1, obtaining a covariance matrix .sub.nn from .sub.(xy)n, and solving n feature values .sub.1.sub.2 . . . .sub.n0 of the covariance matrix .sub.nn and corresponding feature vectors t.sub.1, . . . , t.sub.m, . . . t.sub.n; Step 33C2, obtaining a score matrix T.sub.nm=(t.sub.1, . . . , t.sub.m) of the principal components, wherein m is determined by formula .sub.k=1.sup.m.sub.k/.sub.k=1.sup.n.sub.k.sub.m, .sub.m85%; and Step 33C3, obtaining the principal component vectors .sub.ij=T.sub.nm.sup.T.Math.a.sub.ij of the n-dimensional vectors a.sub.ij, wherein mn, and T.sub.nm.sup.T is a transposition of T.sub.nm.

7. The air kerma conventional true value determining method of claim 6, wherein step 34 comprises the following steps: Step 34A, obtaining a data matrix sample (K.sub.ij,.sub.ij,K.sub.j).sub.(xy)(m+2) via the experiments of the probes of the x sample dosimeters under the y kinds of V.sub.j conditions in step 32D, and obtaining a prediction model K.sub.ij=f.sub.3(.sub.ij,K.sub.j) of K.sub.ij by adopting a support vector machine regression method.

8. The air kerma conventional true value determining method of claim 7, wherein in step 34A, the specific method of obtaining a prediction model K.sub.ij=f.sub.3(.sub.ij,K.sub.j) of K.sub.ij by adopting a support vector machine regression method comprises the steps of: selecting radial basis function as a kernel function in the model training process, parameters of the kernel function are determined by a cross validation method; when the model is constructed, allocating the sample data (K.sub.ij,.sub.ij,K.sub.j).sub.(xy)(m+2) to training sets and testing sets according to a certain proportion; and when the test error is not more than 5%, ending the training, and determining the prediction model K.sub.ij=f.sub.3(.sub.ij,K.sub.j).

9. The air kerma conventional true value determining method of claim 8, wherein the certain proportion refers to that the proportion of the training set to the testing set is more than or equal to 1:1.

10. The air kerma conventional true value determining method of claim 9, wherein step 4 comprises the following steps: Step 41, putting the reference point of the probe of the dosimeter being verified at the point of test; Step 42, selecting a proper radiation source and putting it into an isotope radiation source accommodating device, and adjusting an attenuator to obtain the proper intensity V.sub.j of the incident ray beam; and Step 43, measuring the scattering gamma spectrum by using the gamma spectrometer, the scattering gamma spectrum is applied as input of the prediction model K.sub.ij=f.sub.3(.sub.ij,K.sub.j)=f.sub.3[f.sub.2(S.sub.ij),K.sub.j]=f.sub.1(S.sub.ij,K.sub.j) to obtain an air kerma conventional true value at the point of test in the MRR.

11. The air kerma conventional true value determining method of claim 8, wherein step 4 comprises the following steps: Step 41, putting the reference point of the probe of the dosimeter being verified at the point of test; Step 42, selecting a proper radiation source and putting it into an isotope radiation source accommodating device, and adjusting an attenuator to obtain the proper intensity V.sub.j of the incident ray beam; and Step 43, measuring the scattering gamma spectrum by using the gamma spectrometer, the scattering gamma spectrum is applied as input of the prediction model K.sub.ij=f.sub.3(.sub.ij,K.sub.j)=f.sub.3[f.sub.2(S.sub.ij),K.sub.j]=f.sub.1 (S.sub.ij,K.sub.j) to obtain an air kerma conventional true value at the point of test in the MRR.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a structural schematic diagram of a minitype reference radiation in an embodiment of the present invention;

(2) FIG. 2 is a scattering gamma spectrum at a monitor point when dosimeter BH3103A in MMR in the embodiment of the present invention;

(3) FIG. 3 is a schematic diagram of a linear combination coefficient of principle components extracted from scattering gamma spectrum at the monitor point in the embodiment of the present invention.

(4) Among them, 1 is a shielding box, 2 is a dosimeter being verified, 3 is incident rays, 4 is a radiation source, 5 is a test hole, 6 is a point of test, 7 is a monitor point, 8 is an incident hole, and 9 is a gamma spectrometer.

DETAILED DESCRIPTION OF THE EMBODIMENTS

(5) The technical solution of the present invention will be described in detail below in combination with the embodiment and the accompanying drawings.

(6) An air kerma conventional true value determining method of the present invention is that: firstly, establishing a minitype reference radiation (MRR), which comprises a shielding box and a gamma spectrometer. The shielding box is positioned horizontally and an incident hole is set on the side thereof for incidence of incident rays. The point of test is set in the direction of the incident rays in the shielding box. There is also a test hole on the upper surface of the shielding box, by which the probe of a dosimeter being verified can be put into the shielding box. The reference point of the probe should be coincided with the point of test. Further a monitor point is set in the shielding box. The shielding box is segmented into two parts by one plane perpendicular to the connecting line of the incident hole and the point of test. The monitor point is located at the part adjacent to the incident hole in the shielding box where directly irradiate incident rays are avoided. A gamma spectrometer is disposed in the shielding box. The reference point on the probe of the spectrometer is coincided with the monitor point and the probe is fixed in the shielding box. Next, a proper radiation source and source intensity being selected to provide incident rays for the shielding box. And then a plurality of gamma ray dosimeters to be selected as samples for training a prediction model of the air kerma conventional true value at the point of test. Lastly, the probes of dosimeters being verified are arranged at the point of test, then measuring the scattering gamma spectrum with a gamma spectrometer, with the prediction model and the gamma spectrum as input, the air kerma conventional true value at the point of test is obtained.

The Embodiment

(7) In this embodiment, the structural schematic diagram of the minitype reference radiation (MRR) is shown as FIG. 1, the minitype reference radiation comprises a shielding box 1 with a dimension not more than 1.5 meters and a gamma spectrometer 9. The shielding box 1 is positioned horizontally and provided with an incident hole 8 on the side thereof for the incidence of incident rays 3. A point of test 6 is arranged in the direction of the incident rays 3 in the shielding box 1, a test hole 5 is further provided on the upper surface of the shielding box 1, by which a probe of an dosimeter 2 being verified can be put into the shielding box 1 and can make the reference point of the probe being coincided with the point of test 6. A monitor point 7 is further arranged in the shielding box 1, the shielding box 1 is segmented into two parts by one plane that is perpendicular to the connecting line of the incident hole 8 and the point of test 6, and the monitor point 7 is located at the part adjacent to the incident hole 8 in the shielding box 1 and at the position not directly irradiated by the incident rays 3. The gamma spectrometer 9 is disposed in the shielding box 1. The reference point on the probe thereof is coincided with the monitor point 7, and is fixed in the shielding box 1.

(8) In this embodiment, the shielding box 1 could be a cube with a sectional size of 1 meter, e.g., a 1 m1 m1 m sized cube, and could also be a cuboid or in other shape, the specific size being determined by the total weight of the MRR which allowed by the intended use. The incident hole 8 could be located in the center position of the side of the shielding box, the point of test 6 could also be located in the geometrical center of the shielding box, and the monitor point 7 is generally located at the inner bottom of the shielding box 1.

(9) In use, the specific method includes the following steps:

(10) Step 1, establishing the aforesaid minitype reference radiation device.

(11) Step 2, selecting a proper radiation source and source intensity for providing incident rays for the shielding box.

(12) Step 3, selecting a plurality of gamma ray dosimeters as samples for training a prediction model of the air kerma conventional true value at the point of test in the MRR.

(13) This step includes the following specific steps:

(14) Step 31, selecting a plurality of gamma ray dosimeters as samples for training a prediction model of the air conventional true value, and the types of the plurality of selected gamma ray dosimeters can be BH3103A, FJ317E, SSM-1, FD-3013B, CIT-2000FX.Math., Inspector1000 and Canberra Radiagem2000;

(15) Step 32, measuring the air kerma conventional true value at the point of test when no sample dosimeter is arranged, coinciding the reference point of the probe of a sample dosimeter with the point of test, then measuring the air kerma conventional true value at the point of test by adopting an instrument transfer method, measuring the scattering gamma spectrum of the monitor point by a gamma spectrometer, The specific method is as follows:

(16) Step 32A, putting a standard graphite cavity ionization chamber at the point of test, and measuring the air kerma conventional true value K.sub.j at the point of test, when the intensity of an incident ray beam is V.sub.j;

(17) Step 32B, putting the reference point of the probe of the i.sup.th sample dosimeter at the point of test, setting the intensity of the incident ray beam as V.sub.j, recording the indicate value of the sample dosimeter R.sub.ij and measuring the gamma spectrum S.sub.ij at the monitor point with a gamma spectrometer;

(18) Step 32C, putting the sample dosimeter in the standard reference radiation to search a point where the indicate value of the dosimeter is R.sub.ij, the air kerma conventional true value corresponding to the point being the air kerma conventional true value at the point of test K.sub.ij;

(19) Step 32D, sequentially putting the reference points of the probes of the x sample dosimeters in the point of test, repeating steps 32A to 32C under y kinds of Vj conditions to obtain xy groups of K.sub.ij, S.sub.ij and K.sub.j data, and obtaining the function relationship K.sub.ij=f.sub.1(S.sub.ij,K.sub.j) thereof;

(20) Step 33, acquiring a dose feature value by adopting a principal component analysis method according to the gamma spectrum. The specific method is as follows:

(21) Step 33A, dispersing each acquired S.sub.ij according to a certain energy interval E to obtain a counting rate .sub.ijn array corresponding to the energies of the scattering gamma rays, and constructing n-dimensional vectors a.sub.ij of the counting rates using the energies of the scattering gamma rays as research objects; here, the certain energy interval E refers to E=1500/(1282.sup.z) keV, and z is an integer more than or equal to 0 and less than or equal to 4;

(22) Step 33B, constructing a scattering gamma spectrum counting rate data matrix sample .sub.(xy)n via the experiments of the probes of the x sample dosimeters under the y kinds of V.sub.j conditions in step 32D;

(23) Step 33C, Extracting the principal components of the n-dimensional vectors a.sub.ij by adopting a principal component analysis method to obtain the principal component vectors .sub.ij=T.sub.nm.sup.T.Math.a.sub.ij of the n-dimensional vectors a.sub.ij, wherein mn, T.sub.nm.sup.T is a transposition of T.sub.nm, and T.sub.nm refers to a covariance matrix .sub.nn obtained from .sub.(xy)n; and solving a score matrix composed of m first feature vectors of the covariance matrix .sub.nn. The specific method is as follows:

(24) Step 33C1, obtaining a covariance matrix from .sub.(xy)n, and solving n feature values .sub.1.sub.2 . . . .sub.n0 of the covariance matrix .sub.nn and corresponding feature vectors t.sub.1, . . . , t.sub.m, . . . t.sub.n;

(25) Step 33C2, a score matrix of the principal components is T.sub.nm=(t.sub.1, . . . , t.sub.m), wherein m is determined by formula .sub.k=1.sup.m.sub.k/.sub.k=1.sup.n.sub.k.sub.m, and .sub.m85%;

(26) Step 33C3, the principal component vector of the n-dimensional vectors a.sub.ij is .sub.ij=T.sub.nm.sup.T.Math.a.sub.n, wherein mn, and T.sub.nm.sup.T is a transposition of T.sub.nm;

(27) Step 33D, obtaining the function relationship .sub.ij=f.sub.2(S.sub.ij) between .sub.ij and S.sub.ij, thus obtaining K.sub.ij=f.sub.3(.sub.ij,K.sub.j)=f.sub.3[f.sub.2(S.sub.ij),K.sub.j];

(28) Step 34, obtaining a prediction model of the air kerma conventional true value at the point of test by adopting a support vector machine method. The specific method is as follows:

(29) Step 34A, obtaining a data matrix sample (K.sub.ij,.sub.ij,K.sub.j).sub.(xy)(m+2) via the experiments of the probes of the x sample dosimeters under the y kinds of V.sub.j conditions in step 32D, and obtaining a prediction model K.sub.ij=f.sub.3(.sub.ij,K.sub.j) of K.sub.ij by adopting a support vector machine regression method. The specific method is as follows: training a kernel function selected by a regression prediction model as a radial basis kernel, the parameter of the kernel function is determined by a cross validation method; when the model is established, allocating the data samples (K.sub.ij,.sub.ij,K.sub.j).sub.(xy)(m+2) to a training set and a testing set according to a certain proportion; and when the test error is not more than 5%, ending the training, and determining the prediction model as K.sub.ij=f.sub.3(.sub.ij,K.sub.j). Herein, the certain proportion refers to that the proportion of the training set to the testing set is more than or equal to 1:1.

(30) Step 4, putting the probe of the dosimeter being verified at the point of test, then monitoring the scattering gamma spectrum with a gamma spectrometer, according to the prediction model, the air kerma conventional true value at the point of test is obtained;

(31) Step 4 includes the following specific steps:

(32) Step 41, putting the reference point of the probe of the dosimeter being verified at the point of test of test;

(33) Step 42, selecting a proper radiation source and putting it into an isotope radiation source accommodating device, and adjusting attenuation rate to obtain the proper intensity V.sub.j of the incident ray beam;

(34) Step 43, measuring the scattering gamma spectrum with a gamma spectrometer, inputting the gamma spectrum to the prediction model K.sub.ij=f.sub.3(.sub.ij,K.sub.j)=f.sub.3[f.sub.2(S.sub.ij),k.sub.j]=f.sub.1(S.sub.ij,K.sub.j), and the air kerma conventional true value at the point of test in the MRR is obtained.

(35) The method may further include the following step:

(36) Obtaining a correction factor =K.sub.ij/R.sub.ij in combination with the indicate value R.sub.ij of the dosimeter being verified.

(37) The energy response characteristic of the dosimeter being verified can also be obtained via the above method by switching the radiation sources with different energies; or the angle response data of dosimeter being verified can be obtained by rotating the probe thereof, and other verification items stipulated in JJG393-2003 can also be realized.

(38) A specific example is as follows:

(39) A radioisotope source .sup.137Cs is selected as the radiation source 4 in this embodiment to provide a radiation ray source for the minitype reference radiation (MRR), and a calibration device for calibration of a gamma ray radiation protection instrument is constructed, the structure thereof being shown in FIG. 1. According to the requirement of radiation protection, the shielding box is made of a material such as lead, tungsten alloy or the like with a proper thickness, thus ensuring personnel security when the device is used.

(40) The shielding box 1 is in the shape of a cube having a side length of 1 m, and the geometric center thereof is set as a point of test 12. The incident hole 8 having the diameter of 120 mm and used for the incidence of gamma rays is arranged in the geometric center of the side which is close to the radiator (the radiation source 4), of the shielding box 1. The shielding box 1 is segmented into two parts by one plane that is perpendicular to the connecting line of the incident hole 8 and the point of test 6, and the monitor point 7 is located at the part adjacent to the incident hole 8 on the bottom center line of the shielding box 1 and spaced 100 mm from the projection point at the point of test 6 on the bottom; the test hole 5 having the diameter of 200 mm is arranged at the top of the shielding box 1, and is used for putting the probe of the instrument 2 being verified; scattering gamma ray spectrum in the box are measured by an Inspector1000 portable gamma spectrometer of Canberra company, the reference point of the probe of the gamma spectrometer 9 is aligned with the monitor point 7 on the bottom of the shielding box 1, and the probe of the gamma spectrometer 9 is fixed.

(41) The activity degree of the .sup.137Cs radioisotope source is 1 Ci, the incident ray beam 3 is provided for the shielding box 1 via the device 4 such as a radiator or the like, and the attenuation times of the incident ray beam 3 is adjusted according to the source intensity of the radiation source and the range of the dosimeter being verified. The times of an attenuator is adjusted according to the range of common gamma ray dose (rate) meter and dose equivalent (rate) meter to obtain five experiment source intensities V.sub.j, (j=1, 2, . . . , 5), and the range of the dose rate is 65 Gy/h-1.25 mGy/h.

(42) An experiment is carried out according to the method of this embodiment, a prediction model of the air kerma conventional true value at the point of test 6 is obtained, and the specific implementation steps are as follows:

(43) Step A

(44) A PTW-32005 standard graphite cavity ionization chamber is arranged at the point of test 6 of the shielding box, and the air kerma value at the point of test 6 K.sub.j is measured when the source intensity is V.sub.j.

(45) Step B

(46) Totally nine different types of common gamma ray dose (rate) meters BH3103A, FJ317E, SSM-1, FD-3013B, CIT-2000FX.Math., Inspector1000 (containing two kinds of probes: an IPRON-3 probe and an IPROS-2 probe) and Canberra Radiagem2000 are selected as samples for the experiment, and consecutively numbered as 1, 2 . . . 9, i.e., i=1 . . . 9.

(47) Step C

(48) The probe of the above No. 1 instrument is vertically arranged into the shielding box. And the reference point on the probe is coincided with the point of test 6. The source intensity V.sub.j is sequentially switched for measurement, and recording the indicate value R.sub.1j of No. 1 instrument and the gamma spectrum S.sub.1j recorded by the Inspector1000. FIG. 2 shows the scattering gamma spectrum of No. 1 instrument S.sub.1j in the shielding box under five different intensities V.sub.j.

(49) Step D

(50) The No. 1 instrument is arranged in the standard radiation of -ray air kerma (protection level) measurement standards of an ionizing radiation metrology station of China Academy of Engineering Physics to search a point P.sub.1j having the indicate value R.sub.1j. The air kerma conventional true value at the point P.sub.1j is obtained according to the existing parameters of the standard radiation. And the value is the air kerma conventional true value K.sub.1j at the point of test 6 of the shielding box when the probe of No. 1 instrument is arranged at the point of test 6 of the shielding box under the source intensity V.sub.j in step C.

(51) Step E

(52) The probes of No. 2 to No. 9 instruments are respectively arranged at the point of test 6 in the shielding box, and 45 groups of K.sub.ij, S.sub.ij and K.sub.j data can be obtained by repeating steps C and D under five V.sub.j conditions. The data has a function relationship K.sub.ij=f.sub.1(S.sub.ij,K.sub.j), which is a mathematic prediction model for predicting the air kerma conventional true value of the probe of the instrument being verified at the point of test 6 of the shielding box in the method of the present invention.

(53) Step F

(54) The S.sub.ij is dispersed according to an energy interval 3 keV to obtain an array of 512 counting rates .sub.ijn corresponding to the energies of the scattering gamma rays. According to S.sub.ij features, in order to reduce the dimension of calculating data, the first 250 counting rates having obvious features are selected as valid data, and 250-dimensional vectors a.sub.ij of the counting rates using the energies of the scattering gamma rays as a research object are constructed. Then a scattering gamma spectrum counting rate data matrix sample .sub.45250 is constructed via the experimental data of the probes of the nine gamma dose (rate) meters under the five source intensity V.sub.j conditions. The principal components of 45 pieces of 250-dimensional vectors a.sub.ij are solved by adopting a principal component analysis (PCA) method, i.e., first a covariance matrix .sub.250250 is obtained from .sub.45250, and then 250 feature values .sub.1.sub.2 . . . .sub.2500 of the covariance matrix .sub.250250 and the corresponding feature vectors t.sub.1, . . . , t.sub.m, . . . t.sub.n are solved. The score matrix of the principal components is T.sub.250m=(t.sub.1, . . . , t.sub.m), wherein m is determined by formula .sub.k=1.sup.m.sub.k/.sub.k=1.sup.n.sub.k.sub.m. The principal component score matrix of the 250-dimensional vectors a.sub.ij is T.sub.nm(t.sub.1, . . . , t.sub.m), mn. When .sub.m is 90%, m=2. The linear combination coefficient of the score vectors t.sub.1 and t.sub.2 of two principal components is shown as FIG. 3.

(55) Step G

(56) According to step F, a function relationship .sub.ij=f.sub.2(S.sub.ij) between .sub.ij and S.sub.ij can be obtained. The simulated prediction model K.sub.ij=f.sub.1(S.sub.ij,K.sub.j) in step E can be simplified into K.sub.ij=f.sub.3(.sub.ij,K.sub.j) by replacing S.sub.ij with .sub.ij. Moreover, a data matrix sample (K.sub.ij,.sub.ij,K.sub.j).sub.45(m+2) can be obtained via experiments by using the nine gamma ray dose (rate) meters under the conditions of five different radiation source intensities V.sub.j.

(57) Step H

(58) Based on the data matrix sample (K.sub.ij,.sub.ij,K.sub.j).sub.45(m+2) obtained via experiments, a prediction model K.sub.ij=f.sub.3(.sub.ij,K.sub.j) of K.sub.ij is obtained by adopting a least squares support vector machine (LS-SVM, an improved form of SVM) regression method in this embodiment.

(59) The prediction model is trained on a Matlab software platform for the Windows7

(60) system. And the version of the Matlab software is 2012a. A radial basis function

(61) $K (x, x_{i}) = \exp (- \frac{{.Math. x - x_{i} .Math.}^{2}}{2^{2}})$
is selected as the kernel function of the model by calling a least squares support vector machine toolbox (LS-SVMlab Toolbox User's Guide version 1.5) in the platform. And the parameter .sup.2 of the kernel function and the regularization parameter c are determined by an L-fold cross validation method. L is set to be equal to 10. And the data sample (K.sub.ij,.sub.ij,K.sub.j).sub.45(m+2) is allocated to a training set and a testing set according to a proportion of 6:3. And the training is ended when the test error is less than or equal to 5%. The prediction model of K.sub.ij is K.sub.ij=F[(.sub.ij,K.sub.j), (,K)]+b is finally acquired, wherein F is the kernel function, a and b are parameters of the model, .sub.ij is the principle component vector of the energy spectrum S.sub.ij when the dosimeter being verified is introduced into the shielding box, And K.sub.j is the air kerma value at the point of test of the shielding box when no probe is introduced under the source intensity, and K are sample data of the principle component vector of the energy spectrum for training the model and air kerma sample data at the point of test in the shielding box when no probe is introduced. In combination with the function .sub.ij=f.sub.2(S.sub.ij), the model can be expressed as K.sub.ij=F[(f.sub.2(S.sub.ij),K.sub.j), (,K)]+b, i.e., K.sub.ij=f.sub.1(S.sub.ij,K.sub.j).

(62) When the BH3103A gamma ray dose rate meter 2 being verified is calibrated, a probe of the BH3103A is put into the shielding box, and the reference point of the probe is coincide with the point of test 6 of the MRR; a proper radiation source intensity V.sub.j is determined according to the range of the BH3103A in a manner of selecting an attenuator or the like so that the indicate value of the BH3103A is nearby the midpoint of the calibration range, scattering gamma spectrum measured by the gamma spectrometer 9 are recorded, the principle component vector .sub.i of the spectrum data is extracted and introduced into the prediction model established K.sub.ij=F[(.sub.ij,K.sub.j), (,K)]+b, the air kerma conventional true value K.sub.ij at the point of test 6 of the MRR under such condition is 91.27 Gy/h, the mean R of indicate values of the five dosimeters is 89.82 Gy/h, a calibration factor is obtained according to formula

(63) $= \frac{K_{ij}}{\overline{R}} = 1.016,$
and calibration of the dosimeter is thus realized.

(64) The aforesaid embodiment is only an example for realizing the present invention, and the present invention can be realized in multiple ways. For example, the shape of minitype reference radiation MRR is not limited to a cube, A MRR in other shape such as a cuboid or the like does not influence the effect of the present invention, and the methods of introducing gamma rays via the shielding box and limiting the gamma rays into a small closed space are all implementations of the present invention; the point of test and the monitor point are not limited to the positions in the embodiment, as long as they are located in the MRR, can fulfill the purposes required by the claims and do not influence the effect of the present invention; as for the SVM method for establishing the prediction model K.sub.ij=f.sub.1(S.sub.ij,K.sub.j) of the air kerma conventional true value at the point of test in the MRR, the SVM method has multiple forms and is rapidly developed, the SVM regression mode is not limited to the least squares support vector machine LS-SVM used in this embodiment, and other modes of SVM, C-SVM, v-SVM and the like adopting an SMO (Sequential Minimal Optimization) algorithm are available for fulfilling the purpose of establishing the prediction model K.sub.ij=f.sub.1(S.sub.ij,K.sub.j) of the air kerma conventional true value at the point of test in the MRR.

(65) Other than one .sup.137Cs cesium source for calibration of the gamma ray dose (rate) dosimeter in this embodiment, .sup.137Cs, .sup.241Am and .sup.60Co sources and the method introduced in the present invention can also be simultaneously adopted to obtain the indicators of energy response, angle response and the like of the gamma ray dose (rate) dosimeter. An X ray machine serving as a ray source and the method of the present invention can also be adopted for verification and calibration of gamma and X ray dose (rate) dosimeters.

(66) Although the content of the present invention has been introduced in detail via the above preferred embodiment, the above introduction shall not be regarded as a limitation to the present invention. It would be obvious for a person having professional knowledge and skills to make various modifications, substitutions and avoidances to the present invention upon reading the above content. Therefore, the protection scope of the present invention should be defined by the appended claims.

Air kerma conventional true value determining method

Assignee

Inventors

Cpc classification

Classification Explorer

G01N2223/05

PHYSICS

Classification Explorer

G01T1/02

PHYSICS

Classification Explorer

G01N2223/638

PHYSICS

Classification Explorer

G01N23/20

PHYSICS

Classification Explorer

G01T7/005

PHYSICS

International classification

Classification Explorer

G01T1/00

PHYSICS

Classification Explorer

G01T1/02

PHYSICS

Classification Explorer

G01N23/20

PHYSICS

Abstract

Claims

Description