AIR KERMA CONVENTIONAL TRUE VALUE MEASURING METHOD

Abstract

A measurement method for an air kerma conventional true value comprises: building a small-scale reference radiation field, then selecting a proper radiation source (4) and a source intensity for providing incident rays for a shielding box (1), subsequently selecting a plurality of gamma ray dose measurement instruments as experiment samples for building a prediction model to obtain a prediction model of the air kerma conventional true value of a check point, fmally placing a probe of an instrument to be detected on the check point (6), recording a scattering gamma spectrum detected by a gamma-ray spectrometer (9), and importing the prediction model to obtain the air kerma conventional true value. The method relates to the field of radiation protection detection or calibration, and has the beneficial effects that the result is accurate, the reference radiation field used is small in size, and the method is applied to measurement of the air kerma conventional true value. The method solves the problem that site and in-situ detection or calibration is unlikely to be implemented as the existing standard reference radiation field is too large in space and volume to move or is difficult to move.

Claims

1. An air kerma conventional true value measuring method, comprising the following steps: step 1, establishing a small-scale reference radiation field, the small-scale reference radiation field 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 check point 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 instrument to be detected can be put into the shielding box, a reference point on the probe being superposed with the check point, a dose feature 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 check point and containing the check point, the dose feature point being located at a part close 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 being superposed with the dose feature point and the probe being fixed in the shielding box; step 2, selecting a proper radiation source and source strength to provide incident rays for the shielding box; step 3, selecting a plurality of gamma ray dose measurement instruments as experimental samples for constructing a prediction model to obtain the prediction model of the air kerma conventional true value at the check point; and step 4, putting the probe of the instrument to be detected at the check point, recording scattering gamma spectra measured by the gamma spectrometer, and introducing the scattering gamma spectra to the prediction model to obtain an air kerma conventional true value.

2. The air kerma conventional true value measuring method of claim 1, wherein step 3 comprises the following specific steps: step 31, selecting a plurality of gamma ray dose measurement instruments as experimental samples for constructing a prediction model; step 32, measuring the air kerma conventional true value at the check point when no experimental sample is put, putting a reference point of a probe of an experimental sample on the check point, measuring the air kerma conventional true value at the check point by adopting an instrument transfer method, and acquiring gamma energy spectra of the dose feature point via the gamma spectrometer; step 33, acquiring a dose feature value by adopting a principal component analysis method according to the gamma energy spectra; and step 34, obtaining a prediction model of the air kerma conventional true value at the check point 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 check point, and measuring the air kerma conventional true value K.sub.j′ at the check point when a strength of an incident ray beam is V.sub.j; step 32B, putting the reference point of the probe of the i.sup.th experimental sample at the check point, setting the strength of the incident ray beam as V.sub.j, recording a reading of the experimental sample R.sub.ij and obtaining the gamma energy spectrum S.sub.ij of the dose feature point at the moment via the gamma spectrometer; step 32C, putting the experimental sample in the standard reference radiation field to look for a point having the reading 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 check point; and step 32D, sequentially putting the reference points of the probes of the x experimental samples in the check point, and repeating steps 32A to 32C under y source strength conditions to obtain x×y groups of K.sub.ij, S.sub.ij and K.sub.j′ data for constructing a model of a 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, scattering each acquired S.sub.ij according to a certain energy interval ΔE to obtain a counting rate η.sub.ijn array corresponding to the energy of the scattering gamma ray, and constructing n-dimensional vectors a.sub.ij of counting rates using the energy of the scattering gamma ray as a research object; step 33B, constructing a scattering gamma energy spectrum counting rate data matrix sample Φ.sub.(x×y)×n via experiments of the probes of the x experimental samples under the y source strength conditions in step 32D; step 33C, solving principal 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.n×m.sup.T.Math.a.sub.ij of the n-dimensional vectors a.sub.ij, m≦n, T.sub.n×m.sup.T being a transposition of T.sub.n×m, and T.sub.n×m referring to obtaining a covariance matrix from ξ.sub.n×n from Φ.sub.(x×y)×n; and solving a score matrix composed of m first feature vectors of the covariance matrix ξ.sub.n×n; and step 33D, obtaining a 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′].

5. The air kerma conventional true value measuring method of claim 4, wherein in step 33A, the certain energy interval ΔE refers to: ΔE=1500/(128×2.sup.z)keV, 0≦z≦4, z being an integer.

6. The air kerma conventional true value measuring method of claim 4, wherein step 33C comprises the following specific steps: step 33C1, obtaining a covariance matrix ξ.sub.n×n from Φ.sub.(x×y)×n, and solving n feature values λ.sub.1≧λ.sub.2≧ . . . ≧λ.sub.n≧0 of the covariance matrix ξ.sub.n×n and corresponding feature vectors t.sub.1, . . . , t.sub.m, . . . t.sub.n; step 33C2, obtaining a score matrix T.sub.n×m=(t.sub.1, . . . , t.sub.m) of the principal components, wherein m is determined by formula Σ.sub.k=1.sup.nλ.sub.k≧δ.sub.m, δ.sub.m≧85%; and step 33C3, obtaining the principal component vectors ψ.sub.ij=T.sub.n×m.sup.T.Math.a.sub.ij of the n-dimensional vectors a.sub.ij, wherein m≦n, and T.sub.n×m.sup.T is a transposition of T.sub.n×m.

7. The air kerma conventional true value measuring 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.(x×y)×(m+2) via the experiments of the probes of the x experimental samples 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 measuring 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: 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 constructed, allocating the data samples (K.sub.ij,ψ.sub.ij,K.sub.j′).sub.(x×y)×(m+2) to a training set and a test set 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 measuring method of claim 8, wherein the certain proportion refers to that the proportion of the training set to the test set is more than or equal to 1:1.

10. The air kerma conventional true value measuring method of claim 9, wherein step 4 comprises the following steps: step 41, putting the reference point of the probe of the instrument to be detected at the check point; step 42, selecting a proper radiation source and putting it into an isotope radiation source accommodating device, and adjusting an attenuator to obtain the measured strength V.sub.j of the incident ray beam; and step 43, measuring the scattering gamma spectrum at the moment by using the gamma spectrometer, and introducing the scattering gamma spectrum data into the prediction model constructed 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 check point.

11. The air kerma conventional true value measuring method of claim 8, wherein step 4 comprises the following steps: step 41, putting the reference point of the probe of the instrument to be detected at the check point; step 42, selecting a proper radiation source and putting it into an isotope radiation source accommodating device, and adjusting an attenuator to obtain the measured strength V.sub.j of the incident ray beam; and step 43, measuring the scattering gamma spectrum at the moment by using the gamma spectrometer, and introducing the scattering gamma spectrum data into the prediction model constructed 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 check point.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0042] FIG. 1 is a structural schematic diagram of a small-scale reference radiation field in an embodiment of the present invention;

[0043] FIG. 2 is a scattering gamma energy spectrum diagram of a gamma dose instrument BH3103A in MMR at a feature dose point in the embodiment of the present invention;

[0044] FIG. 3 is a schematic diagram of a linear combination coefficient of principle components of scattering gamma energy spectra at the feature dose point in the embodiment of the present invention.

[0045] Among them, 1 is a shielding box, 2 is an instrument to be detected, 3 is incident ray, 4 is a radiation source, 5 is a test hole, 6 is a check point, 7 is a dose feature point, 8 is an incident hole, and 9 is a gamma spectrometer.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0046] The technical solution of the present invention will be described in detail below in combination with an embodiment and the accompanying drawings.

[0047] An air kerma conventional true value measuring method of the present invention is that: establishing a small-scale reference radiation field (MRR) first, the small-scale reference radiation field comprising a shielding box and a gamma spectrometer, the shielding box being positioned horizontally and provided with an incident hole on the side the shielding box for the incidence of incident rays, a check point being arranged in the direction of the incident rays in the shielding box, the shielding box being further provided with a test hole on the upper surface through which a reference point on a probe of an instrumentation to be detected can be superposed with the check point, a dose feature point being further arranged in the shielding box, the shielding box being segmented into two parts by one plane that is perpendicular to the connecting line of the incident hole and the check point and that contains the check point, the dose feature point being located at the part close to the incident hole in the shielding box and at the position not directly irradiated by the incident rays, a gamma spectrometer being arranged in the shielding box, the reference point on a probe of the gamma spectrometer being superposed with the dose feature point and the probe being fixed in the shielding box; then selecting a proper radiation source for providing incident rays for the shielding box, selecting a plurality of gamma ray dose measurement instruments as experimental samples for establishing a prediction model to obtain a prediction model of the air kerma conventional true value at the check point, and finally during a test: putting the probe of the instrument to be detected at the check point, recording the scattering gamma spectra measured by the gamma spectrometer, and introducing the prediction model to obtain an air kerma conventional true value.

The Embodiment

[0048] In this embodiment, the structural schematic diagram of the small-scale reference radiation field (MRR) is shown as FIG. 1, and the small-scale reference radiation field comprises a shielding box 1 having a side length 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 check point 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, a probe of an instrumentation 2 to be detected can be put into the shielding box 1 via the test hole 5 to make a reference point on the probe being superposed with the check point 6. A dose feature 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 check point 6 and that contains the check point 6, and the dose feature point 7 is located at the part close 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 arranged in the shielding box 1, a reference point on a probe thereof is superposed with the dose feature point 7, and the gamma spectrometer 9 is fixed in the shielding box 1.

[0049] In this embodiment, the shielding box 1 may be a cube having a sectional size of 1 meter, e.g., a 1 m×1 m×1 m sized cube, and may also be a cuboid or in other shape, the specific size being determined by the total weight of the MRR allowed by the intended use. The incident hole 8 may be located in the center position of the side of the shielding box, the check point 6 may also be located in the geometrical center of the shielding box, and the dose feature point 7 is generally located on the inner bottom of the shielding box 1.

[0050] In use, the specific method includes the following steps:

[0051] Step 1, constructing the aforesaid small-scale reference radiation field device.

[0052] Step 2, selecting a proper radiation source and source strength for providing incident rays for the shielding box.

[0053] Step 3, selecting a plurality of gamma ray dose measurement instruments as experimental samples for establishing a prediction model to obtain the prediction model of the air kerma conventional true value at the check point.

[0054] This step includes the following specific steps:

[0055] step 31, selecting a plurality of gamma ray dose measurement instruments as experimental samples for establishing a prediction model, and the plurality of selected gamma ray dose measurement instruments may be the gamma ray dose measurement instruments of types BH3103A, FJ317E, SSM-1, FD-3013B, CIT-2000FX.γ, Inspector1000 and Canberra Radiagem2000;

[0056] step 32, measuring the air kerma conventional true value at the check point when no experimental sample is arranged, putting a reference point of a probe of an experimental sample on the check point, measuring the air kerma conventional true value at the check point by adopting an instrument transfer method, and acquiring gamma energy spectra of the dose feature point via the gamma spectrometer. The specific method is as follows:

[0057] step 32A, putting a standard graphite cavity ionization chamber at the check point, and measuring the air kerma conventional true value at the check point, denoted as K.sub.j′, when the strength of an incident ray beam is V.sub.j;

[0058] step 32B, putting the reference point of the probe of the i.sup.th experimental sample at the check point, setting the strength of the incident ray beam as V.sub.j, recording the reading of the experimental sample R.sub.ij and obtaining the gamma energy spectrum of the dose feature point at the moment via the gamma spectrometer S.sub.ij;

[0059] step 32C, putting the experimental sample in the standard reference radiation field to look for a point having the reading equal to R.sub.ij, the air kerma conventional true value corresponding to the point being the air kerma conventional true value at the check point K.sub.ij;

[0060] step 32D, sequentially putting the reference points of the probes of the x experimental samples in the check point, repeating steps 32A to 32C under y kinds of Vj conditions to obtain x×y groups of K.sub.ij, S.sub.ij and K.sub.j′ data, and obtaining a function relationship K.sub.ij=f.sub.1(S.sub.ij, K.sub.j′) thereof;

[0061] step 33, acquiring a dose feature value by adopting a principal component analysis method according to the gamma energy spectrum. The specific method is as follows:

[0062] step 33A, scattering each acquired S.sub.ij according to a certain energy interval ΔE to obtain a counting rate η.sub.ijn array corresponding to the energy of the scattering gamma ray, and constructing n-dimensional vectors a.sub.ij of the counting rates using the energy of the scattering gamma ray as a research object; here, the certain energy interval ΔE refers to ΔE=1500/(128×2.sup.z)keV, and z is an integer more than or equal to 0 and less than or equal to 4;

[0063] step 33B, constructing a scattering gamma energy spectrum counting rate data matrix sample Φ.sub.(x×y)×n via the experiments of the probes of the x experimental samples under the y kinds of V.sub.j conditions in step 32D;

[0064] step 33C, solving the principal 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.n×m.sup.T.Math.a.sub.ij of the n-dimensional vectors a.sub.ij, wherein m≦n, T.sub.n×m.sup.T is a transposition of T.sub.n×m, and T.sub.n×m refers to obtaining a covariance matrix ξ.sub.n×n from Φ.sub.(x×y)×n; and solving a score matrix composed of m first feature vectors of the covariance matrix ξ.sub.n×n. The specific method is as follows:

[0065] step 33C1, obtaining a covariance matrix ξ.sub.n×n from Φ.sub.(x×y)×n, and solving n feature values λ.sub.1≧λ.sub.2≧ . . . ≧λ.sub.n≧0 of the covariance matrix ξ.sub.n×n and corresponding feature vectors t.sub.1, . . . , t.sub.m, . . . t.sub.n;

[0066] step 33C2, a score matrix of the principal components is T.sub.n×m=(t.sub.1, . . . , t.sub.m), wherein m is determined by the formula Σ.sub.k=1.sup.mλ.sub.k/Σ.sub.k=1.sup.nλ.sub.k≧δ.sub.m, and δ.sub.m≧85%;

[0067] step 33C3, the principal component vectors of the n-dimensional vectors a.sub.ij is ψ.sub.ij=T.sub.n×m.sup.T.Math.a.sub.ij, wherein m≦n, and T.sub.n×m.sup.T is a transposition of T.sub.n×m;

[0068] step 33D, obtaining a 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′];

[0069] step 34, obtaining regressionally a prediction model of the air kerma conventional true value at the check point by adopting a support vector machine method. The specific method is as follows:

[0070] step 34A, obtaining a data matrix sample (K.sub.ij,ψ.sub.ij,K.sub.j′).sub.(x×y)×(m+2) via the experiments of the probes of the x experimental samples 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.(x×y)×(m+2) to a training set and a test 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 test set is more than or equal to 1:1.

[0071] Step 4, putting the probe of the instrument to be detected at the check point, recording scattering gamma spectra measured by the gamma spectrometer, and introducing the scattering gamma spectra to the prediction model to obtain an air kerma conventional true value.

[0072] Step 4 includes the following specific steps:

[0073] step 41, putting the reference point of the probe of the instrument to be detected at the check point;

[0074] step 42, selecting a proper radiation source and putting it into an isotope radiation source accommodating device, and adjusting an attenuator to obtain the strength V.sub.j of the measured incident ray beam; and

[0075] step 43, measuring the scattering gamma spectrum at the moment by using the gamma spectrometer, and introducing the scattering gamma spectrum data into the established prediction model K.sub.ij=f.sub.3(ψ.sub.ij,K.sub.j′)=f.sub.1(S.sub.ij,K.sub.j′) to obtain an air kerma conventional true value at the check point.

[0076] The methohd may further include the following step:

[0077] obtaining a correction factor ω=K.sub.ij/R.sub.ij in combination with the reading of the instrument to be detected R.sub.ij.

[0078] The energy response characteristic of the instrument to be detected can also be obtained via the above method by switching the radiation sources with different energy; or the angle response data of the instrument to be detected can be obtained by rotating the probe of the instrument to be detected, and other verification items stipulated in JJG393-2003 can also be realized.

[0079] A specific example is as follows:

[0080] A radioisotope source .sup.137Cs is selected as the radiation source 4 in this embodiment to provide a radiation ray source for the small-scale 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 having a proper thickness, thus ensuring personnel safety when the device is used.

[0081] 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 check point 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 check point 6 and that contains the check point 6, and the dose feature point 7 is located at the part close 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 check point 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 to be detected; scattering gamma ray spectra in the box are measured by using an Inspector1000 portable gamma spectrometer of Canberra company, the reference point of the probe of the gamma spectrometer 9 is aligned with the dose feature point 7 on the bottom of the shielding box 1, and the probe of the gamma spectrometer 9 is fixed.

[0082] 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 strength of the radiation source and the range of the instrument to be detected. The times of an attenuator is adjusted according to the range of the common gamma ray dose (rate) and dose equivalent (rate) instrument to obtain five experiment source strengths V.sub.j,(j=1, 2, . . . , 5), and the range of the dose rate is 65 μGy/h−1.25 mGy/h.

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

[0084] Step A

[0085] A PTW-32005 standard graphite cavity ionization chamber is arranged at the check point 6 of the shielding box, and the air kerma value at the check point 6 K.sub.j′ when the source strength is V.sub.j is measured.

[0086] Step B

[0087] Totally nine different types of common gamma ray dose (rate) instruments 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.

[0088] Step C

[0089] The probe of the above No. 1 instrument is vertically arranged into the shielding box, and the reference point on the probe is superposed with the check point 6. The source strength V.sub.j is sequentially switched for measurement, and the readings of the No. 1 instrument and R.sub.1j the gamma spectra S.sub.1j recorded by the Inspector1000 are recorded. FIG. 2 shows the scattering gamma spectra of the No. 1 instrument S.sub.1j in the shielding box and under five source strengths V.sub.j.

[0090] Step D

[0091] The No. 1 instrument is arranged in the standard radiation field of “γ-ray air kerma (protection level) measurement standards” of an ionizing radiation metrology station of China Academy of Engineering Physics to look for a point P.sub.ij having the reading R.sub.1j. The air kerma conventional true value at the point P.sub.ij is obtained according to the existing parameters of the standard radiation field, and the value is the air kerma conventional true value K.sub.1j at the check point 6 of the shielding box when the probe of the No. 1 instrument is arranged at the check point 6 of the shielding box under the source strength V.sub.j in step C.

[0092] Step E

[0093] The probes of the No. 2 to No. 9 instruments are respectively arranged at the check point 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 above 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 to be detected at the check point 6 of the shielding box in the method of the present invention.

[0094] Step F

[0095] The S.sub.ij is scattered according to an energy interval 3 keV to obtain an array of 512 counting rates η.sub.ijn corresponding to the energy of the scattering gamma ray. According to S.sub.ij features, in order to reduce the volume 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 energy of the scattering gamma ray as a research object are constructed; then a scattering gamma energy spectrum counting rate data matrix sample Φ.sub.45×250 is constructed via the experimental data of the probes of the nine gamma dose (rate) instruments under the five source strength 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., a covariance matrix ξ.sub.250×250 is obtained first from Φ.sub.45×250, and then 250 feature values λ.sub.1≧λ.sub.2≧ . . . ≧λ.sub.250≧0 of the covariance matrix ξ.sub.250×250 and the corresponding feature vectors are solved. The score matrix of the principal components is T.sub.250×m=(t.sub.1, . . . , t.sub.m), wherein m is determined by formula Σ.sub.k=1.sup.mλ.sub.k≧δ.sub.m. The principal component score matrix of the 250-dimensional vectors a.sub.ij is T.sub.n×m=(t.sub.1, . . . , t.sub.m), m≦n. 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.

[0096] Step G

[0097] 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) instruments under five different radiation source strength V.sub.j conditions.

[0098] Step H

[0099] 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.

[0100] The prediction model is trained on a Matlab software platform for the Windows7 system, and the version of the Matlab software is 2012a. A radial basis function

[00003]$K\left(x,x_i\right)=\exp \left(-\frac{{.Math.x-x_i.Math.}^2}{2.Math..Math.{\sigma }^2}\right)$

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 test set according to a proportion of 6:3; and 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 and K.sub.j′ are respectively the principle component vector of the energy spectrum S.sub.ij when the instrument to be detected is introduced into the shielding box and the air kerma value at the check point of the shielding box when no probe is introduced under the source strength, ψ′ 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 check point 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″)]′×αa+b, i.e., K.sub.ij=f.sub.1(S.sub.ij,K.sub.j′).

[0101] When the BH3103A gamma ray dose rate instrument 2 to be detected is calibrated, a probe of the BH3103A is put into the shielding box, and the reference point of the probe is superposed with the check point 6 of the MRR; a proper radiation source strength 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 reading of the BH3103A is nearby the midpoint of the calibration range, scattering gamma spectra 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 check point 6 of the MRR under such condition is 91.27 μGy/h, the mean R of readings of the five instruments is 89.82 μGy/h, a calibration factor is obtained according to formula

[00004]$\omega =\frac{K_{\mathrm{ij}}}{\stackrel{\_}{R}}-1.016,$

and calibration of the instrument is thus realized.

[0102] 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 the small-scale reference radiation MRR is not limited to a cube, the 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 check point and the dose feature 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 check point in the MRR, the SVM method has multiple forms and is being 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 check point in the MRR.

[0103] Other than one .sup.137Cs cesium source for calibration of the gamma ray dose (rate) instrument 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) instrument. 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) instruments.

[0104] 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.