删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

基于虚拟源原理的源边界参数蒙特卡罗反演技术

本站小编 Free考研考试/2021-12-29

摘要:在就地γ谱仪搜索扫描测量“热粒子”、“放射性汇集点”、“放射性汇集区”过程中, 只能给出污染源的大概位置, 不能给出源的污染深度等边界参数. 本文主要对虚拟技术在就地γ谱仪搜索扫描测量细化污染源边界中的应用进行了研究. 将就地γ谱仪测量对象简化成衰减层 + 放射性热区(测量目标源) + 衰减层 + 干扰源的四层理论模型, 运用虚拟技术将源项层虚拟成点源, 进一步简化了理论模型, 使用蒙特卡罗方法模拟计算探测效率与峰谷比等参数, 最后使用最小二乘法使模拟计算结果反演逼近源项实际参数, 从而建立了源边界参数反演计算的理论方法及步骤. 理论研究和实验结果一致, 验证了所建立的计算模型和技术方法是正确可靠的. 目前, 对于均匀分布的放射性核素, 该技术已经能够准确确定污染区域深度分布等边界参数, 从而在治理时达到废物处置减容的目的. 同时, 该技术对于禁核试核查目标核弹头惰层厚度参数的确定也具有重大的参考价值.
关键词: 虚拟点源/
参数表征/
蒙特卡罗模拟/
源边界

English Abstract


--> --> -->
“热粒子”、“放射性汇集点”、“放射性汇集区”的监测是目前国际上核辐射探测领域重点关注的参数, 其中放射性源边界是最为关键的参数. 目前国内外现有的就地γ谱仪搜索扫描测量及伽玛相机成像技术, 只能给出污染源的大概位置, 而对于源边界参数的确定目前研究较少, 尤其是无法给出核素深度方向的参数; 通过出壳γ谱反推核弹头结构材料厚度是一件比较复杂和难度较大的工作, 需要的参数多, 解的不确定性大[1]. 因此, 亟待发展新的测量方法和实验技术来解决这些问题. 放射性源边界测量的准确性最为关键的是探测效率因子这个参数的准确刻度[2,3], 而国内外针对这个问题开展的研究也非常多, 在蒙特卡罗模拟刻度中主要使用点源实验表征晶体参数[4-10], 提高了测量的准确性; 在虚拟刻度方面, 国内外开展了虚拟点探测器理论研究[11-15], 虚拟点探测器模型由Noteal[15]最早提出, 在处理探测器效率刻度问题中该模型被许多研究者使用并进行了大量的实验验证. 在测量源与探测器之间存在一定距离时可用虚拟的点探测器来等效大体积的探测器, 此时探测器峰效率将较精确地满足平方反比定律, 这样, 对于一确定尺寸的探头, 虚拟点探测器位置仅随光子能量的变化而变化. 根据平方反比定律, 可求得空间任意一点处点源的峰效率, 然后可对点源效率积分求解得到体源的探测效率, 即使用点源完成了复杂体源的刻度, 简化和方便了仪器刻度工作. 本文反向使用蒙特卡罗模拟表征晶体参数的技术, 结合虚拟点探测器理论及虚拟源刻度方法等[16,17], 根据测量能谱分析, 就可以确定放射性污染区域深度方向分布参数, 同时对反解核弹头惰层厚度具有参考意义.
HPGe探测器晶体一般为圆柱状, 射线与晶体的作用可以认为是射线与晶体内一虚拟点探测器作用, 即将晶体等效为一个虚拟的点. 这个虚拟的点位置可能不止一个, 但是为了简化模型和计算方便, 一般是在探测器中心对称轴线找到一个且唯一一个代表性的点(代表点), 该点距离晶体上表面的距离又称有效作用深度. 虚拟点探测器: 对于辐射探测器, 在探测器中心对称轴线上存在有一个唯一的代表点位置, 其点探测器对参考点位置点源的探测效率与探测器对点源的探测效率相等且满足平方反比定律, 即该点探测器叫探测器的虚拟点探测器. 虚拟点探测器刻度方法首先用于点源探测效率的刻度, 主要借用参考点位置的标准点源探测效率值通过平方反比定律来计算空间不同位置处点源的探测效率.
根据虚拟点探测器的理论, 将探测器测量对象也进行虚拟化, 定义虚拟点源概念. 虚拟点源: 对于辐射探测器测量对象为单个体源或单个面源的情况下, 在探测器中心对称轴线上有一个唯一的代表点位置, 该位置的放射性点源全能峰探测效率与所对应的放射性体源或面源的探测效率相等, 该点源叫做体源或面源的虚拟点源. 在多数情况下, 测量对象可能不止一个源, 有可能是多个源(体源、面源和点源)组成, 为了区分单个源的情况, 在这里引入等效虚拟点源. 等效虚拟点源: 对于辐射探测器测量对象为混合源(多个体源、面源和点源组成)情况下, 在探测器中心对称轴线上有一个唯一的代表点位置, 该位置的放射性点源全能峰探测效率与混合源探测效率相等, 该点源叫作混合源的等效虚拟点源. 它和虚拟点源的区别是测量对象为多个源, 其他定义都一致.
在无源效率刻度中, 蒙特卡罗模拟表征晶体参数的研究工作非常多, 技术方法也非常成熟. 其基本原理是: 通过标准点源或体源(源的形状大小、材料和活度参数等完全已知)固定位置实验, 获取探测器的实验效率刻度值, 再使用蒙特卡罗方法建立理论计算模型, 调整理论计算模型中的晶体参数使其计算效率值和实验值一致, 从而达到表征晶体参数的目的. 在这一表征的过程中测量对象源的参数是已知的, 测量仪器晶体参数是未知的, 但在实际测量中经常会遇到大量测量对象未知的情况. 针对这种情况, 本文采用逆向思维, 提出源边界参数表征技术, 即晶体参数已知, 源边界参数未知. 其基本原理和晶体参数表征技术类似: 即首先通过实验方法获得测量对象的峰谷比及全能峰探测效率的实验值, 再通过蒙特卡罗方法建立理论计算模型, 为了优化理论计算模型, 使用虚拟技术将测量对象虚拟成点源, 计算调整理论模型中的源边界参数并借助最小二乘方法使计算得到的峰谷比及效率值不断逼近实验值. 具体过程是: 污染沙土的放射性源项层可能会很复杂, 可能有很多各种各样的源组成, 各个源之间又有沙土衰减层, 但不管怎么复杂, 都将它们统一简化成如图1所示的四层理论模型.
图 1 源边界参数反演理论模型
Figure1. The inversion theory model of source boundary parameters.

关心的测量对象目标源一般为体源(目标体源), 又称放射性“热区”, 若直接裸露于地表, 就没有图1中最上一层沙土衰减层. 同时, 热区下方可能存在其他放射源, 我们称之为干扰源, 干扰源是一个相对概念, 当关心它的参数时, 也可以定为目标源. 热区和干扰源之间有衰减层, 干扰源在一些源项简单的测量中可能就不存在. 本文的目标就是在有干扰源存在的情况下反演出热区位置及厚度参数. 具体分三步反演.
第一步: 等效虚拟点位置的实验反演计算. 使用虚拟技术将热区和深层干扰源虚拟成一个等效虚拟点源, 该虚拟点源包含了所有源项的辐射信息. 对于钚污染土, 241Am和239Pu的单位体积活度比一般为常数, 设为$Q ={{{A_{239}}}}/{{{A_{241}}}}$, 其中A表示活度(Bq), 其与探测效率的关系推导为:
${A_{241}} = \frac{{{N_{241}}}}{{t \times {\varepsilon _{241}} \times {P_{{\rm{\gamma }}241}}}},$
${A_{239}} = \frac{{{N_{239}}}}{{t \times {\varepsilon _{239}} \times {P_{{\rm{\gamma }}239}}}},$
$Q = \frac{{{\varepsilon _{241}}}}{{{\varepsilon _{239}}}} \times \frac{{{P_{{\rm{\gamma }}241}} \times {N_{239}}}}{{{P_{{\rm{\gamma }}239}} \times {N_{241}}}},$
其中t为测量时长(s), N表示特征射线全能峰计数, ε表示探测效率, 241表示241Am, 239表示239Pu, ${P_{\rm{\gamma }}}$表示射线分支比. 241Am核素的59.54 keV (分支比0.359, 26.4 keV分支比0.024)和239Pu核素的51.62 keV (分支比0.000271, 129 keV分支比0.000062)射线全能峰计数一般通过实际测量能谱得到, 分支比也是常数, 因此从(3)式可以看出活度比和探测效率比是一一对应的, 因此, 在图1中对称轴上(使用MCNP程序模拟计算对称轴上不同位置处点源探测效率)总能找到一个等效虚拟点位置, 该处的射线全能峰探测效率比算得的239Pu/241Am活度比应与给出的一致. 该步反演除了得到等效虚拟点位置外还得到三个参数: 实测59.54 keV射线的峰谷比以及等效虚拟点处59.54 keV和26.4 keV或129 keV射线全能峰探测效率, 这些参数是下一步反演计算的基础. 其中59.54 keV射线的峰谷比设为X1, 59.54 keV射线全能峰探测效率参数设为X2, 26.4 keV或129 keV射线全能峰探测效率参数设为X3, E表示根据实验获得的数学期望值, 则它们的数学期望分别表示为:
$\frac{{{N_{\rm m}}}}{{{N_{\rm {vm}}}}} = E({X_1}),$
$\varepsilon ({h_{\rm{e}}}) = E({X_2}),$
$\varepsilon {({h_{\rm{e}}})^*} = E({X_3}),$
v表示谷区, 对于59.54 keV射线一般取54—57 keV这个区间, h表示虚拟点位置(虚拟点位置一般在源对称轴线上变化), m表示实测谱值, e表示等效虚拟点, 上标*用于区分虚拟点位置处不同能量射线全能峰探测效率. (5)和(6)式左边的虚拟点源位置he处的探测效率值是通过核素活度比值(3)式定出, 峰谷比通过能谱获得.
第二步: 热区和深层干扰源的虚拟点源位置反演计算. 将热区和深层干扰源分别虚拟成其虚拟点源, 如图1中所示, 用a表示图1中上虚拟点位置(热区虚拟点位置), b表示下虚拟点位置(深层干扰源虚拟点位置), 等效虚拟点是由这两个虚拟点按不同的权重组合而成的. 设59.54 keV和26.4 keV或129 keV射线在对称轴上不同位置按不同权重组合模拟计算结果设为Xi, 表示为:
$\frac{N}{{{N_{\rm v}}}}({h_{\rm{a}}}){w_{\rm{a}}} + \frac{N}{{{N_{\rm v}}}}({h_{\rm{b}}})(1 - {w_{\rm{a}}}) = {X_1},$
$\varepsilon ({h_{\rm{a}}}){w_{\rm{a}}} + \varepsilon ({h_{\rm{b}}})(1 - {w_{\rm{a}}}) = {X_2},$
$\varepsilon {({h_{\rm{a}}})^*}{w_{\rm{a}}} + \varepsilon {({h_{\rm{b}}})^*}(1 - {w_{\rm{a}}}) = {X_3},$
w表示权重因子; 然后结合第一步获取的E(Xi)参数使用最小二乘法计算均方偏差$\sigma (X)$, 如(10)式所示, 找出均方偏差$\sigma (X)$最小的组合(其对应的热区和深层干扰源的虚拟点源位置就是所求)即可.
$\sigma (X) = \sqrt {\frac{1}{3}\sum\limits_{i = 1}^3 {{{\big[({X_i} - E({X_i}))/E({X_i})\big]}^2}} } .$
该步反演结果是分别得到热区和深层干扰源的虚拟点源位置以及它们的射线全能峰探测效率及峰谷比信息.
第三步: 根据第二反演得到热区虚拟点源位置反演计算热区参数. 使用MCNP程序模拟计算热区源项不同参数组合下(不同衰减层厚、不同源厚两两组合)的射线全能峰探测效率和能谱的峰谷比参数, 利用(10)式计算, 找出均方偏差ó最小的组合(其对应的衰减层厚度和热区厚度就是所求)即可.
根据图1的理论计算模型, 设计了两种实验探测模式. 探测模式1: 75 mm × 25 mm 239Pu体源下方放置1.6 cm厚沙土衰减层然后放置在距离平板型高纯锗探测器铍窗0.953 cm位置处, 源与衰减层密度都为1.34 g/cm3, 如图2所示.
图 2 探测模式1
Figure2. The detection mode 1.

探测模式2: 241Am点源通过支架放置在距离2.0 cm处的0.5 cm厚沙土衰减层上再将其放置在75 mm × 25 mm厚241Am体源上, 然后将它们整个放置在平板型高纯锗探测器铍窗进行放射性测量, 源与沙土地密度都为1.43, 如图3所示.
图 3 探测模式2
Figure3. The detection mode 2.

谱获取时间足够长, 保证使其统计涨落控制在1.0%以内, 实验数据见表1表2 (表1表2中分别给出了体源和点源单独测量数据, 只是为了获得两种探测模式的实验参数, 因此不涉及计算参数的数据, 在此未提供).
测量对象测量时长t/105 sN241 (54—57 keV)N241 (59.54 keV)N239 (51.62, 129 keV)A/104 Bq
241Am239Pu
探测模式13.16432513625339979
239Pu体源2.0075248200239711, 527174.5618.7


表1探测模式1实验能谱峰计数及处理结果
Table1.Energy peak count of experimental spectrum and process results for detection mode 1.

测量对象测量时长 t/105 sN241 (26.4 keV)N241 (54—57 keV)N241 (59.54 keV)A/× 104 Bq
26.4 keV59.54 keV
探测模式24.002400509531180659645368.168.91
241Am点源0.5654160024653314567.748.38
241Am体源0.8004534627457730.4230.532


表2探测模式2实验能谱峰计数及处理结果
Table2.Energy peak count of experimental spectrum and process results for detection mode 2.

将计数通过计算转化为活度, 活度计算结果见表1, 其中探测效率数据是由蒙特卡罗模拟计算给出. 体源效率: 51.62 keV为2.03 × 10–2, 59.54 keV为2.3 × 10–2, 129 keV为2.72 × 10–2, 其中活度计算结果18.7 × 104 Bq是取了239Pu核素两条射线计算结果的平均值, A239/A241 = 18.7/4.56 = 4.1 (对于未知源项目标该比值一般是通过取样测量得到, 或根据已有资料查询得到). 241Am源谷区选择54—57 keV能区, 探测模式1下峰谷比为25339979/4325136 = 5.9.
表2中, 241Am点源效率26.4 keV为3.66 × 10–2, 59.54 keV为4.16 × 10–2; 241Am体源效率26.4 keV为4.44 × 10–3, 59.54 keV为2.26 × 10–2. 在探测模式2下, 通过241Am两个射线能量计算其活度, 其总活度为点源活度加体源活度, 如26.4 keV射线能量计算总活度为8.16 × 104 (Bq) = (7.74 + 0.423) × 104 (Bq), 59.54 keV射线能量计算总活度结果为8.91 × 104 (Bq) = (8.38 + 0.532) × 104 (Bq), 且A241 (26.4 keV)/A241 (59.54 keV)比值理论上应该是1.0, 但由于存在系统误差, 241Am (26.4 keV)/241Am (59.54 keV)比值计算结果为8.16/8.91 = 0.92, 对于这个系统误差可以不用校正, 因为后面计算都采用该比值. 探测模式2下, 峰谷比为65964536/9531180 = 6.9.
根据实验设计及探测器参数, 建立了如图4图5所示的蒙特卡罗计算模型(图中的数字及位置坐标表示尺寸参数). 探测器参数: 晶体角死层被两个球及一个圆柱所切得到, 球心在Z轴上坐标分别为1.90和3.44, 球半径分别为1.90和2.06. 图中单位为cm. 铍窗密度1.85, 沙土成分由下列元素(符号)组成: O(52.2%), Si(30.6%), Al(5.49%), Fe(2.21%), Ca(8.04%), Mg(1.40%), 括号内为质量分数.
图 4 探测模式1的MCNP程序计算模型
Figure4. Calculation model of MCNP procedure for detection mode 1.

图 5 探测模式2的MCNP程序计算模型
Figure5. Calculation model of MCNP procedure for detection mode 2.

对于探测器模式1, 只有图1中沙土衰减层和体源, 没有干扰源的存在. 但模拟的过程中都是假设有干扰源存在的, 下面介绍其模拟计算过程.
2
4.1.等效虚拟点源
-->在对称轴上, 计算不同位置处点源探测效率, 241Am主要计算其59.54 keV射线探测效率, 239Pu主要计算其129 keV射线探测效率, 其计算结果如表3中数据所示.
h/cm${\varepsilon _{241}}(h)$/10–3${\varepsilon _{239}}(h)$/10–3A241/104 BqA239/105 BqQ${N}/{ { {N_{\rm v}} } }(h)$
–1.2524.324.60.9210.4975.48.05
–1.6017.919.51.240.6295.06.99
–2.0012.915.11.730.8104.76.12
–2.409.4511.92.361.034.45.44
–2.807.009.453.191.304.14.92
–3.205.247.604.261.613.84.51
–3.603.976.165.621.983.54.18
–3.803.475.566.452.203.44.04
–4.003.035.047.372.433.33.90


表3等效虚拟点源探测效率及峰谷比
Table3.The detection efficiency and peak/valley of equivalent virtual point source.

根据表3中探测效率值, 结合表1中探测模式1的能峰数据计算不同位置坐标处的Q活度比值, 可以看出, 在坐标为–2.80处时Q比值为4.1, 和表1中活度比值一致, 因此, 可以认为等效虚拟点源(如图4中)的位置坐标为(0 0 –2.80). 等效虚拟点源包含了两部分的等效, 一部分是测量对象体源的贡献, 一部分是干扰源的贡献(干扰源可以是点源、面源和体源, 统一将其转化成虚拟点源), 将这两部分各自虚拟成点源, 可以认为等效虚拟点源由体源虚拟点源和干扰源虚拟点源按不同的权重组合而成, 因此下一步的目标就是找出这个权重和体源的虚拟点源位置坐标.
2
4.2.体源虚拟点源
-->等效虚拟点源的探测效率值7.0 × 10–3和9.45 × 10–3即为(5)式和(6)式的数学期望值, 峰谷比5.9为(4)式的数学期望值. 根据表3中不同点位置处的探测效率, 计算不同虚拟点位置和不同权重组合下的等效探测效率及等效峰谷比, 计算结果见表4. 根据这些值按(7)—(9)式计算其理论值Xi, 最后按理论计算公式(10)计算均方偏差, 如表4所列.
wh/cm$\varepsilon (h)$/10–3${\varepsilon ^*}(h)$/10–3${N}/{ { {N_{\rm v}} } }(h)$X2/10–3X3/10–3X1$\sigma (X)$
0.10–1.2524.324.68.05
0.90–3.205.247.604.517.159.304.870.177
0.90–3.603.976.164.186.008.014.560.308
0.90–3.803.475.564.045.557.474.440.385
0.90–4.003.035.043.905.166.994.310.458
0.10–1.6017.919.56.99
0.90–3.205.247.604.516.518.794.760.217
0.90–3.603.976.164.185.367.504.460.396
0.90–3.803.475.564.044.916.964.330.479
0.90–4.003.035.043.904.526.484.210.554
0.10–2.0012.915.16.12
0.90–3.205.247.604.516.018.354.670.277
0.90–3.603.976.164.184.867.064.370.474
0.90–3.803.475.564.044.416.524.240.559
0.90–4.003.035.043.904.026.044.120.635
0.10–2.409.4511.95.44
0.90–3.205.247.604.515.668.034.610.327
0.90–3.603.976.164.184.526.744.300.530
0.90–3.803.475.564.044.066.204.180.616
0.90–4.003.035.043.903.675.724.050.693


表4不同组合下等效虚拟点的均方偏差计算数据
Table4.The mean square deviation calculation data of equivalent virtual point at different combination.

表4中热区虚拟点源位置坐标为–1.25, –1.60, –2.00和–2.40 cm, 干扰源的虚拟点源位置坐标为–3.20, –3.60, –3.80和–4.00 cm. 省去中间过程, 直接计算得到虚拟点坐标、权重及均方偏差的一一对应的表5.
h/cmw$\sigma (X)$w$\sigma (X)$w$\sigma (X)$w$\sigma (X)$w$\sigma (X)$
–1.250.100.200.300.400.50
–3.200.900.1770.800.3550.700.6640.600.9880.501.32
–3.600.900.3080.800.2230.700.5080.600.8480.501.20
–3.800.900.3850.800.2010.700.4480.600.7920.501.15
–4.000.900.4580.800.2100.700.3990.600.7440.501.11
–1.600.100.200.300.400.50
–3.200.900.2170.800.1950.700.3600.600.5680.500.784
–3.600.900.3960.800.2100.700.2340.600.4350.500.669
–3.800.900.4790.800.2620.700.2030.600.3840.500.622
–4.000.900.5540.800.3180.700.1960.600.3430.500.583
–2.000.100.200.300.400.50
–3.200.900.2770.800.1870.700.1770.600.2560.500.371
–3.600.900.4740.800.3280.700.2100.600.1810.500.272
–3.800.900.5590.800.3990.700.2570.600.1780.500.238
–4.000.900.6350.800.4650.700.3070.600.1930.500.215
–2.400.100.200.300.400.50
–3.200.900.3270.800.2640.700.2090.600.1730.500.167
–3.600.900.5300.800.4350.700.3440.600.2610.500.195
–3.800.900.6160.800.5100.700.4070.600.3090.500.225
–4.000.900.6930.800.5770.700.4640.600.3560.500.257


表5不同组合下等效虚拟点的均方偏差计算数据
Table5.The mean square deviation calculation data of equivalent virtual point at different combination.

h/cmw$\sigma (X)$w$\sigma (X)$w$\sigma (X)$w$\sigma (X)$
–1.250.600.700.800.9
–3.200.401.650.301.980.202.310.102.64
–3.600.401.550.301.900.202.260.102.61
–3.800.401.510.301.880.202.240.102.60
–4.000.401.480.301.850.202.220.102.60
–1.600.600.700.800.90
–3.200.401.000.301.230.201.450.101.67
–3.600.400.9100.301.160.201.400.101.65
–3.800.400.8720.301.130.201.380.101.64
–4.000.400.8390.301.100.201.370.101.63
–2.000.600.700.800.90
–3.200.400.4980.300.6300.200.7630.100.898
–3.600.400.4100.300.5610.200.7160.100.875
–3.800.400.3750.300.5330.200.6980.100.865
–4.000.400.3470.300.5090.200.6810.100.857
–2.400.600.700.800.90
–3.200.400.1940.300.2440.200.3050.100.372
–3.600.400.1690.300.1990.200.2670.100.351
–3.800.400.1740.300.1860.200.2520.100.342
–4.000.400.1860.300.1780.200.2400.100.335




表5中可以看出, 均方偏差最小的值为0.167, 其对应的体源虚拟点位置坐标为–2.40, 干扰源对应的虚拟点位置坐标–3.20, 而其对应的权重都为0.50. 由于两个点位置很近权重又一样, 因此, 可以认为没有干扰源的存在, 实际上等效虚拟点位置坐标就是体源的虚拟点位置坐标.
2
4.3.体源参数反演
-->表3中坐标–2.80对应的探测效率及峰谷比作为(4)—(6)式的数学期望值, 然后计算不同体源中心坐标不同体源厚度下的探测效率及峰谷比值, 将这些值作为理论计算值, 并计算其与数学期望之间的均方偏差, 计算结果如表6所列.
hV/cm体源厚度/cm$\varepsilon ({h_{\rm{V}}})$
/10–3
${\varepsilon ^*}({h_{\rm{V}}})$
/10–2
${N}/{ { {N_{\rm v}} } }({h_{\rm{V} } })$$\sigma (X)$
–2.800.804.540.6894.690.444
–2.801.24.600.6964.740.431
–2.801.64.700.7054.810.414
–2.450.805.680.8185.030.231
–2.451.25.760.8265.070.217
–2.451.65.890.8375.160.200
–2.452.06.050.8515.260.181
–2.452.56.310.8745.410.159
–2.453.06.650.9035.630.163
–2.454.07.570.9816.240.289
–2.454.98.781.087.120.539
–2.000.807.621.035.580.188
–2.001.27.751.045.640.212
–2.001.67.931.055.750.248
–2.002.08.171.075.880.295
–2.002.58.551.106.110.373
–2.003.09.031.146.410.474
–2.004.010.41.257.350.769
–1.500.8010.71.336.430.744
–1.501.210.91.356.540.781
–1.501.611.21.376.680.834
–1.502.011.51.406.900.905
–1.503.012.91.507.791.18
–0.500.8022.02.3410.12.81


表6体源参数的反演计算数据
Table6.The inversion data of volume source parameters.

体源中心坐标用hV表示. 从表6中可以看出, 均方偏差最小的为0.159, 其对应的体源中心坐标为–2.45, 体源厚度为2.5 cm, 和实际情况完全一致, 充分说明了该技术方法的可行性和准确性. 对于探测模式2, 探测器参数和探测模式1一致, 但是沙土密度变为1.43, 整个测量是贴近探测器进行的, 129 keV射线能量变成26.4 keV. 其计算过程如表7表9所列.
h/cm${\varepsilon _{26.4\;{\rm{keV}}}}(h)$/10–3${\varepsilon _{59.54\;{\rm{keV}}}}(h)$/10–2A26.4 keV/104 BqA59.54 keV/104 BqA59.54 keV/A26.4 keV${N}/{ { {N_{\rm v}} } }(h)$
0.8048.20014.40.05190.3186.1017.0
0.4013.5008.790.18500.5232.8012.0
04.0605.560.61500.8261.309.3
–0.202.2604.481.11001.0300.908.4
–0.401.2703.641.97001.2600.607.7
–0.600.7172.983.49001.5400.407.1
–0.800.4092.456.12001.8800.306.6


表7等效虚拟点源探测效率、峰谷比及活度比
Table7.The detection efficiency, peak/valley and acvitiy ratio of equivalent virtual point source.

hV/cm体源
厚度/cm
${\varepsilon ^*}({h_{\rm{V}}})$/
10–2
$\varepsilon ({h_{\rm{V}}})$/
10–3
${N}/{ { {N_{\rm v}} } }({h_{\rm{V} } })$$\sigma (X)$ hV/cm体源
厚度/cm
${\varepsilon ^*}({h_{\rm{V}}})$/
10–2
$\varepsilon ({h_{\rm{V}}})$/
10–3
${N}/{ { {N_{\rm v}} } }({h_{\rm{V} } })$$\sigma (X)$
0.751.003.5410.411.901.5900 0.152.202.455.009.340.2310
0.750.603.478.0611.401.0100 0.151.602.312.698.540.3470
0.750.303.447.2311.200.8060 0.151.002.211.728.070.5920
0.252.002.595.499.610.3530 0.150.402.171.327.860.6920
0.251.902.564.919.440.2110 02.502.264.409.000.0900
0.251.802.544.429.290.0890 02.452.244.148.900.0470
0.251.702.514.009.140.0220 02.402.233.908.830.0640
0.251.602.493.639.030.1090 02.002.132.528.300.3950
0.251.502.473.328.920.1870 01.502.041.607.830.6260
0.251.302.432.818.720.3130 01.001.971.137.540.7470
0.250.802.372.048.380.5080 –0.251.101.650.6026.860.8910
–0.250.501.610.4576.700.9300


表9体源参数的反演计算数据
Table9.The inversion data of volume source parameters.

表8中热区虚拟点源位置坐标为0.80, 0.40, 和0 cm, 干扰源的虚拟点源位置坐标为–0.40, –0.60和–0.80 cm. 从表8中可以看出, 均方偏差最小是0.112, 其对应的上下虚拟点源权重都为0.50, 上下虚拟点源坐标分别为0和–0.8 cm, 等效虚拟点源坐标为–0.20 cm, 比较靠近上虚拟点位置, 两虚拟点位置较远, 权重也占了一半, 因此, 可以认为有一个强的干扰源的存在. 干扰源只需证明它存在或者不存在即可, 它的存在及其设置只是为了更为准确地解出体源边界参数, 在本文中它的边界参数是不必求解的.
h/cmw$\sigma (X)$w$\sigma (X)$w$\sigma (X)$w$\sigma (X)$w$\sigma (X)$
0.800.100.200.300.400.50
–0.400.902.4140.801.7830.702.8070.603.8320.504.86
–0.600.902.1150.801.5120.702.5690.603.6270.504.69
–0.800.901.9510.801.3620.702.4330.603.5100.504.59
0.400.100.200.300.400.50
–0.400.900.4660.800.2910.700.5290.600.7800.501.03
–0.600.900.1770.800.1190.700.2860.600.5670.500.857
–0.800.900.1850.800.2450.700.1750.600.4470.500.753
00.100.200.300.400.50
–0.400.900.1900.800.1600.700.1680.600.2350.500.327
–0.600.900.4310.800.3460.700.2140.600.1230.500.165
–0.800.900.6200.800.5120.700.3500.600.1970.500.112


表8不同组合下等效虚拟点的均方偏差计算数据
Table8.The mean square deviation calculation data of equivalent virtual point at different combination.

h/cmw$\sigma (X)$w$\sigma (X)$w$\sigma (X)$w$\sigma (X)$
0.800.600.700.800.90
–0.400.405.880.306.910.207.930.108.96
–0.600.405.750.306.810.207.870.108.92
–0.800.405.670.306.750.207.830.108.90
0.400.600.700.800.90
–0.400.401.290.301.550.201.800.102.06
–0.600.401.150.301.440.201.730.102.03
–0.800.401.060.301.380.201.690.102.01
00.600.700.800.90
–0.400.400.4270.300.5320.200.6390.100.747
–0.600.400.2870.300.4250.200.5670.100.711
–0.800.400.2070.300.3610.200.5230.100.689




表9中, 均方偏差的计算只用了26.4 keV射线能量探测效率和59.54 keV射线峰谷比这两项, 没有使用59.54 keV射线能量探测效率. 从表9中可以看出, 最小均方偏差为0.0220, 其对应的体源中心坐标为0.25 cm, 体源厚度为1.7 cm, 计算得沙土衰减层厚度为1.53 mm, 由于这一值相对体源厚度太小, 在表层不可能存在这么薄的沙土衰减层, 因此可以将其直接忽略, 并入体源厚度, 可认为体源厚度为2.0 cm, 相对真实值2.5 cm, 大概有20.0%左右的相对偏差, 其偏差稍微偏大的原因是干扰点源使用低能26.4 keV射线穿透3.0 cm左右沙土衰减层后在探测器晶体内产生的有效计数太少, 统计涨落过大, 造成了计算误差也大. 因此, 不能选用能量过低的低能射线, 比如本文选用的26.4 keV能量射线就有些偏低, 会造成结果一定的偏差.
本文提供一种确定放射性可疑区源边界参数的虚拟点源方法, 解决了现有探测方法无法准确判定放射性污染区深度边界参数的技术问题. 该技术方法可应用于多种核素或发射多条特征能量射线核素的源项目标, 其过程是: 首先根据核素活度浓度比值(一种核素该比值可认为是1.0, 但要测量其两条及以上特征能量射线)找出等效虚拟点源, 该虚拟点源特点是由该点算出的活度浓度比值与给出的一致, 它包含了目标体源及其他干扰源的信息; 其次再将该等效虚拟点按一定的权重分成上下两个虚拟点源, 若上下两个虚拟点源与等效虚拟点源距离过近且权重相等, 则不需要将等效虚拟点源分成上下两个虚拟点源, 它本身就是体源的虚拟点源, 不存在干扰源, 若上下两个虚拟点源距离等效点源距离较远或权重不一致, 则可认为上虚拟点源就是目标体源的虚拟点源, 该虚拟点源的探测效率及峰谷比和目标体源的一致; 最后就是根据目标体源的虚拟点源的探测效率及峰谷比计算并反演得到体源边界参数. 该技术的优点是: 1)利用虚拟点探测器理论及虚拟点源技术, 成功地解决了方程组法遇到深层干扰源情况下存在病态或不确定大的难题, 使就地及航测γ谱仪发展为一项更为有效的探测手段; 2)使用了新原理及测量方式, 不增加新的物理探测设备, 大大减少了经费开支, 减少了工作量及时间.
相关话题/计算 测量 干扰 实验 技术

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 锶原子光晶格钟碰撞频移的测量
    摘要:中性原子光晶格钟的系统不确定度评估中,碰撞频移引起的频移修正量和不确定度是其中重要的一项,且其评估结果将直接影响交流斯塔克频移的评估.碰撞频移来源于囚禁在同一个格点里面原子间的相互作用,其大小与原子的密度有关.本文实验测量了国家授时中心87Sr光晶格钟的碰撞频移.利用水平方向的一维光晶格囚禁数 ...
    本站小编 Free考研考试 2021-12-29
  • 基于波长调制技术的高温高压流场温度测量方法
    摘要:温度是衡量燃烧效率的重要参数之一,温度的测量对工业燃烧过程的节能减排控制和发动机状态诊断等都具有重要意义.可调谐半导体吸收光谱技术是一种非侵入式测量技术,具有较强的环境适应性,可实现快速、原位检测.本文基于H2O在7185.6,6807.8以及7444.35/37cm–1三条吸收线集成测量系统 ...
    本站小编 Free考研考试 2021-12-29
  • 宣布式单光子源宣布效率的宣布测量基相关性
    摘要:对非线性晶体参量下转换过程所产生的纠缠光子对的一路光子进行测量,“宣布”式地在另一路获得单光子,是制备单光子源的一种重要途径.其中,核心的问题是如何得到较高的宣布效率.本文以I类相位匹配偏硼酸钡(β-BaB2O4,BBO)非线性晶体参量下转换过程中所产生的偏振纠缠光子对为例,通过研究其纠缠特性 ...
    本站小编 Free考研考试 2021-12-29
  • 磁路和天线位置对2 cm电子回旋共振离子推力器性能影响的实验研究
    摘要:磁路和天线位置对电子回旋共振离子推力器(electroncyclotronresonanceionthruster,ECRIT)的电子加热、等离子体约束和输运过程都有很大的影响,进而影响离子束流引出和中和器耦合电压.本文开展不同磁路和天线位置下2cmECRIT离子束流和耦合电压变化规律的实验研 ...
    本站小编 Free考研考试 2021-12-29
  • 碱金属和碱土金属掺杂二维GaN材料电磁特性的第一性原理计算
    摘要:基于密度泛函理论和投影缀加波赝势方法,采用广义梯度近似算法研究了碱金属(Li,Na,K和Rb)和碱土金属(Be,Mg和Sr)掺杂二维GaN单层的电子结构和磁学性质.研究表明,除Be原子位于GaN单层平面内之外,其余掺杂原子均略微隆起于平面.通过比较七种掺杂体系在不同环境下的形成能,发现在富N环 ...
    本站小编 Free考研考试 2021-12-29
  • 基于紧束缚模型的拓扑物理微波实验验证平台的开发
    摘要:拓扑光子学、拓扑物理与光学的结合,为凝聚态理论的验证以及新型光学器件的构建提供了新的视角.紧束缚模型是凝聚态物理的重要研究手段.我们发现,将传统光子晶体的背景材料由通常的空气改为有效介电常数为负数的材料之后,这样的光子晶体和紧束缚模型有一一对应的关系,可以用于相关理论的验证.通过数值仿真实验, ...
    本站小编 Free考研考试 2021-12-29
  • 有源光纤中稀土离子激光上能级寿命测量的研究
    摘要:提出了一种基于脉冲光纤激光放大器能量特性测量有源光纤中稀土离子上能级寿命的方法.根据光纤激光器速率方程,能够确定有源光纤中反转粒子数储能随抽运功率和时间的变化关系;实验测量不同种子光脉冲重复频率下放大器输出单脉冲能量的变化,可以反映出反转粒子数随时间的变化情况,进而根据理论模型得到激活离子的激 ...
    本站小编 Free考研考试 2021-12-29
  • 半导体激光器储备池计算系统的工作点选取方法
    摘要:半导体激光器储备池计算系统的性能受很多因素的影响,如虚节点间隔、激光器的偏置电流和反馈强度等.对于光注入信号方式,注入强度和频率失谐的大小也会影响系统的性能,使得工作点更难确定.为此,本文以10阶非线性自回归移动平均任务为基础,提出一种选取半导体激光器储备池计算系统的最佳反馈强度与注入强度的方 ...
    本站小编 Free考研考试 2021-12-29
  • 忆容器多谐振荡器及其实验
    摘要:忆容器是一种具有记忆性的非线性电容,为研究忆容器的电路特性,提出了一种压控型忆容器的二次曲线模型,利用电流反馈型运放等器件构建了能够动态模拟忆容器q-v特性的仿真器.通过仿真和实验观测到忆容器的滞回曲线,以及随外加激励频率增加而收缩的特性.分析了周期性激励信号的参数对忆容值取值范围的影响,并对 ...
    本站小编 Free考研考试 2021-12-29
  • 结构光照明技术在二维激光诱导荧光成像去杂散光中的应用
    摘要:背景杂散光信号的干扰制约了激光片光成像技术的发展,本文将结构光照明技术应用到激光片光成像测量中来消除杂散光的干扰.介绍了基于结构光照明技术的工作原理和实验测量系统,基于Matlab软件理论分析了相位移动结构光照明技术具有完全消除杂散光的作用,并针对稳定罗丹明B溶液进行了二维激光诱导荧光成像实验 ...
    本站小编 Free考研考试 2021-12-29