搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

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

田自宁 欧阳晓平 陈伟 王雪梅 邓宁 刘文彪 田言杰

引用本文:
Citation:

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

田自宁, 欧阳晓平, 陈伟, 王雪梅, 邓宁, 刘文彪, 田言杰

Source boundary parameter of Monte Carlo inversion technology based on virtual source principle

Tian Zi-Ning, Ouyang Xiao-Ping, Chen Wei, Wang Xue-Mei, Deng Ning, Liu Wen-Biao, Tian Yan-Jie
PDF
HTML
导出引用
  • 在就地γ谱仪搜索扫描测量“热粒子”、“放射性汇集点”、“放射性汇集区”过程中, 只能给出污染源的大概位置, 不能给出源的污染深度等边界参数. 本文主要对虚拟技术在就地γ谱仪搜索扫描测量细化污染源边界中的应用进行了研究. 将就地γ谱仪测量对象简化成衰减层 + 放射性热区(测量目标源) + 衰减层 + 干扰源的四层理论模型, 运用虚拟技术将源项层虚拟成点源, 进一步简化了理论模型, 使用蒙特卡罗方法模拟计算探测效率与峰谷比等参数, 最后使用最小二乘法使模拟计算结果反演逼近源项实际参数, 从而建立了源边界参数反演计算的理论方法及步骤. 理论研究和实验结果一致, 验证了所建立的计算模型和技术方法是正确可靠的. 目前, 对于均匀分布的放射性核素, 该技术已经能够准确确定污染区域深度分布等边界参数, 从而在治理时达到废物处置减容的目的. 同时, 该技术对于禁核试核查目标核弹头惰层厚度参数的确定也具有重大的参考价值.
    In the in situ γ spectrometer based measurement of " hot particular”, " radioactive collection point” and " radioactive collection area”, only the position of the pollution source can be located roughly, but its boundary parameters such as the thickness of pollution source cannot be given. In this paper, the application of virtual technology to the scanning of γ spectrometer is studied. We convert γ spectrometer measurement objects into a four-layer theoretical model, which are attenuation thickness + radioactive hot area + attenuation thickness + disturb source. Then, the source item layer is virtualized into a point source by using virtual technology. So, the theoretical model is further simplified. Then the detection efficiency and peak/valley ratio parameter of source term are simulated by Monte Carlo method. Finally, the source term parameters are retrieved by using the least square method, and thus establishing the theoretical method and procedure of inversion calculation of source boundary parameters. In this paper, the theoretical and experimental results are shown to be consistent with each other. So, this method is verified to be correct and practicable. Currently, the method can accurately determine the depth distribution parameters of radioactive contamination area for uniformly distributed radio nuclides. In conclusion, the technical achievements can be used to accurately determine the boundary range of the radioactive hot zone, thus realizing the purpose of reducing the waste disposal capacity during the treatment. At the same time, the determination of the inert layer thickness parameters of the target nuclear warhead of Nuclear Test Ban Treaty has a significant reference value.
      通信作者: 田自宁, tzn1019@126.com
    • 基金项目: 国家自然科学基金(批准号: 11405134)资助的课题
      Corresponding author: Tian Zi-Ning, tzn1019@126.com
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 11405134)
    [1]

    熊宗华, 亢武, 龚建, 胡广春, 向永春, 裴永全 2003 物理学报 52 1Google Scholar

    Xiong Z H, Kang W, Gong J, Hu G C, Xiang Y C, Pei Y Q 2003 Atca Phys. Sin. 52 1Google Scholar

    [2]

    Kováik A, Sy′kora I, Povinec P P 2013 J. Radioanal. Nucl. Chem. 298 665Google Scholar

    [3]

    Peyres V, García-Torao E 2007 Nucl. Instr. Meth. Phys. Res. A 580 296Google Scholar

    [4]

    田自宁, 欧阳晓平, 殷经鹏, 张洋, 杨文静 2013 原子能科学技术 47 1411Google Scholar

    Tian Z N, Ouyang X P, Yin J P, Zhang Y, Yang W J 2013 Atomic Energy Sci. Technol. 47 1411Google Scholar

    [5]

    Elanique A, Marzocchi O, Leone D, Hegenbart L, Breustedt B, Oufni L 2012 Appl. Radiat. Isot. 70 538Google Scholar

    [6]

    Budjá D, Heisel M, Maneschg W, Simgen H 2009 Appl. Radiat. Isot. 67 706Google Scholar

    [7]

    Gasparro J, Hult M, Johnston P N, Tagziria H 2008 Nucl. Instr. Meth. Phys. Res. A 594 196Google Scholar

    [8]

    Huya N Q, Binhb D Q, An V X 2007 Nucl. Instr. Meth. Phys. Res. A 573 384Google Scholar

    [9]

    Huy N Q 2010 Nucl. Instr. Meth. Phys. Res. A 621 390Google Scholar

    [10]

    Luís R, Bento J, Carvalhal G, Nogueira P, Silva L, Teles P, Vaz P 2010 Nucl. Instr. Meth. Phys. Res. A 623 1014Google Scholar

    [11]

    Mohammadi M A, Abdi M R, Kamali M, Mostajaboddavati M, Zare M R 2011 Appl. Radiat. Isot. 69 521Google Scholar

    [12]

    Presler O, German U, Pelled O, Alfassi Z B 2004 Appl. Radiat. Isot. 60 213Google Scholar

    [13]

    Mahling S, Orion I, Alfassi Z B 2006 Nucl. Instr. Meth. Phys. Res. A 557 544Google Scholar

    [14]

    熊文彬, 仇春华, 段天英, 刘浩杰, 潘君艳, 陈海涛, 刘进辉 2011 原子能科学技术 45 999

    Xiong W B, Qiu C H, Duan T Y, Liu H J, Pan J Y, Chen H T, Liu J H 2011 Atomic Energy Sci. Technol. 45 999

    [15]

    Noteal A 1971 Nucl. Instr. Meth. 91 513Google Scholar

    [16]

    田自宁, 欧阳晓平, 曾鸣, 成智威 2013 物理学报 62 162902Google Scholar

    Tian Z N, Ouyang X P, Zeng M, Cheng Z W 2013 Acta Phys. Sin. 62 162902Google Scholar

    [17]

    田自宁, 陈伟, 韩斌, 田言杰, 刘文彪, 冯天成, 欧阳晓平 2016 物理学报 65 062901Google Scholar

    Tian Z N, Chen W, Han B, Tian Y J, Liu W B, Feng T C, Ouyang X P 2016 Acta Phys. Sin. 65 062901Google Scholar

  • 图 1  源边界参数反演理论模型

    Fig. 1.  The inversion theory model of source boundary parameters.

    图 2  探测模式1

    Fig. 2.  The detection mode 1.

    图 3  探测模式2

    Fig. 3.  The detection mode 2.

    图 4  探测模式1的MCNP程序计算模型

    Fig. 4.  Calculation model of MCNP procedure for detection mode 1.

    图 5  探测模式2的MCNP程序计算模型

    Fig. 5.  Calculation model of MCNP procedure for detection mode 2.

    表 1  探测模式1实验能谱峰计数及处理结果

    Table 1.  Energy peak count of experimental spectrum and process results for detection mode 1.

    测量对象测量时长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
    下载: 导出CSV

    表 2  探测模式2实验能谱峰计数及处理结果

    Table 2.  Energy peak count of experimental spectrum and process results for detection mode 2.

    测量对象测量时长 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
    下载: 导出CSV

    表 3  等效虚拟点源探测效率及峰谷比

    Table 3.  The detection efficiency and peak/valley of equivalent virtual point source.

    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
    下载: 导出CSV

    表 4  不同组合下等效虚拟点的均方偏差计算数据

    Table 4.  The mean square deviation calculation data of equivalent virtual point at different combination.

    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
    下载: 导出CSV

    表 5  不同组合下等效虚拟点的均方偏差计算数据

    Table 5.  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)$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
    下载: 导出CSV
    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
    下载: 导出CSV

    表 6  体源参数的反演计算数据

    Table 6.  The inversion data of volume source parameters.

    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
    下载: 导出CSV

    表 7  等效虚拟点源探测效率、峰谷比及活度比

    Table 7.  The detection efficiency, peak/valley and acvitiy ratio of equivalent virtual point source.

    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
    下载: 导出CSV

    表 9  体源参数的反演计算数据

    Table 9.  The inversion data of volume source parameters.

    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
    下载: 导出CSV

    表 8  不同组合下等效虚拟点的均方偏差计算数据

    Table 8.  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)$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
    下载: 导出CSV
    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
    下载: 导出CSV
  • [1]

    熊宗华, 亢武, 龚建, 胡广春, 向永春, 裴永全 2003 物理学报 52 1Google Scholar

    Xiong Z H, Kang W, Gong J, Hu G C, Xiang Y C, Pei Y Q 2003 Atca Phys. Sin. 52 1Google Scholar

    [2]

    Kováik A, Sy′kora I, Povinec P P 2013 J. Radioanal. Nucl. Chem. 298 665Google Scholar

    [3]

    Peyres V, García-Torao E 2007 Nucl. Instr. Meth. Phys. Res. A 580 296Google Scholar

    [4]

    田自宁, 欧阳晓平, 殷经鹏, 张洋, 杨文静 2013 原子能科学技术 47 1411Google Scholar

    Tian Z N, Ouyang X P, Yin J P, Zhang Y, Yang W J 2013 Atomic Energy Sci. Technol. 47 1411Google Scholar

    [5]

    Elanique A, Marzocchi O, Leone D, Hegenbart L, Breustedt B, Oufni L 2012 Appl. Radiat. Isot. 70 538Google Scholar

    [6]

    Budjá D, Heisel M, Maneschg W, Simgen H 2009 Appl. Radiat. Isot. 67 706Google Scholar

    [7]

    Gasparro J, Hult M, Johnston P N, Tagziria H 2008 Nucl. Instr. Meth. Phys. Res. A 594 196Google Scholar

    [8]

    Huya N Q, Binhb D Q, An V X 2007 Nucl. Instr. Meth. Phys. Res. A 573 384Google Scholar

    [9]

    Huy N Q 2010 Nucl. Instr. Meth. Phys. Res. A 621 390Google Scholar

    [10]

    Luís R, Bento J, Carvalhal G, Nogueira P, Silva L, Teles P, Vaz P 2010 Nucl. Instr. Meth. Phys. Res. A 623 1014Google Scholar

    [11]

    Mohammadi M A, Abdi M R, Kamali M, Mostajaboddavati M, Zare M R 2011 Appl. Radiat. Isot. 69 521Google Scholar

    [12]

    Presler O, German U, Pelled O, Alfassi Z B 2004 Appl. Radiat. Isot. 60 213Google Scholar

    [13]

    Mahling S, Orion I, Alfassi Z B 2006 Nucl. Instr. Meth. Phys. Res. A 557 544Google Scholar

    [14]

    熊文彬, 仇春华, 段天英, 刘浩杰, 潘君艳, 陈海涛, 刘进辉 2011 原子能科学技术 45 999

    Xiong W B, Qiu C H, Duan T Y, Liu H J, Pan J Y, Chen H T, Liu J H 2011 Atomic Energy Sci. Technol. 45 999

    [15]

    Noteal A 1971 Nucl. Instr. Meth. 91 513Google Scholar

    [16]

    田自宁, 欧阳晓平, 曾鸣, 成智威 2013 物理学报 62 162902Google Scholar

    Tian Z N, Ouyang X P, Zeng M, Cheng Z W 2013 Acta Phys. Sin. 62 162902Google Scholar

    [17]

    田自宁, 陈伟, 韩斌, 田言杰, 刘文彪, 冯天成, 欧阳晓平 2016 物理学报 65 062901Google Scholar

    Tian Z N, Chen W, Han B, Tian Y J, Liu W B, Feng T C, Ouyang X P 2016 Acta Phys. Sin. 65 062901Google Scholar

计量
  • 文章访问数:  6012
  • PDF下载量:  38
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-07-16
  • 修回日期:  2019-09-19
  • 上网日期:  2019-11-27
  • 刊出日期:  2019-12-05

/

返回文章
返回