基于克里金理论的辐射剂量场插值重构与可视化研究

辐射防护 ›› 2019, Vol. 39 ›› Issue (6): 475-482.

基于克里金理论的辐射剂量场插值重构与可视化研究

李华, 赵原, 刘立业, 肖运实, 曹勤剑

中国辐射防护研究院,太原 030006

收稿日期:2018-09-07 出版日期:2019-11-20 发布日期:2020-05-15
通讯作者: 刘立业。E-mail:liuliye@cirp.org.cn
作者简介:李华(1989—),2011年毕业于清华大学工程物理系核工程与核技术专业,2017年博士毕业于清华大学工程物理系核科学与技术专业,副研究员。E-mail:lihua_7559@126.com
基金资助:
国家自然科学基金资助项目(项目编号:11305146、11275174)。

Research on interpolation reconstruction and visualization of radiation dose field based on Kriging theory

LI Hua, ZHAO Yuan, LIU Liye, XIAO Yunshi, CAO Qinjian

China Institute for Radiation Protection, Taiyuan 030006

Received:2018-09-07 Online:2019-11-20 Published:2020-05-15

摘要/Abstract

摘要： 使用基于克里金理论的插值计算方法对均匀网格与随机排列等不同分布类型的辐射剂量场进行了数据重构与可视化研究,并将重构结果与常用的反距离加权插值法进行了对比,分析了不同插值方法对辐射剂量场重构精度与可视化显示效果的影响。结果表明:对于多数情况,克里金插值法在辐射剂量场重构的精度与可视化显示的效果上均明显优于反距离加权插值法,算例中克里金插值法的结果与原始数据之间的偏差在10%以内,而反距离加权插值法的结果与原始数据之间的偏差基本在15%以上;但在辐射剂量场梯度变化较大且已知数据极少的区域,克里金插值法的计算结果偏差会增大,此时可使用反距离加权插值作为其补充,综合使用。本文的研究工作,验证了克里金插值法在辐射剂量场重构与可视化方面的应用潜力,可为推动辐射剂量场可视化仿真技术的发展提供技术支持。

关键词: 剂量场重构, 剂量计算, 克里金理论, 可视化, 插值算法

Abstract: The data reconstruction and visualization of different radiation field distribution types were carried out by interpolation calculation method based on Kriging theory, such as uniform mesh, random arrangement, and so on. The reconstruction results were compared with the commonly used inverse distance weighted interpolation method. The different interpolation methods were used to analyze the effect of reconstruction accuracy and visual display for the radiation field. The results show that for most cases, the Kriging interpolation method is superior to the inverse distance weighted interpolation method in the accuracy of the radiation field reconstruction and the visual display. The deviation between original data and the results of the Kriging interpolation method in the example is within 10%, and the deviation between original data and the results of the IDW is basically above 15%. However, in the region where the gradient of the radiation dose field is large and the known data is very few, the result of the Kriging interpolation method will be biased greatly. In the case where the source position coincides with the interpolation point, the calculation result of the Kriging interpolation method near the source will be biased. In this case, IDW interpolation can be used as a supplement and used comprehensively. The research work in this paper verifies the application potential of Kriging interpolation in radiation field reconstruction and visualization, which provides a technical support for promoting the development of radiation field and dose visualization simulation technology.

Key words: dose field reconstruction, dose calculation, Kriging theory, visualization, interpolation algorithm

中图分类号:

TL72

李华, 赵原, 刘立业, 肖运实, 曹勤剑. 基于克里金理论的辐射剂量场插值重构与可视化研究[J]. 辐射防护, 2019, 39(6): 475-482.

LI Hua, ZHAO Yuan, LIU Liye, XIAO Yunshi, CAO Qinjian. Research on interpolation reconstruction and visualization of radiation dose field based on Kriging theory[J]. RADIATION PROTECTION, 2019, 39(6): 475-482.

参考文献 15

[1]	万海峡. 辐射剂量场可视化平台剂量分布快速计算程序开发[D]. 华北电力大学, 2012.WAN Haixia. Development of rapid calculation code for dose rate distribution in visualization platform of radiation field[D]. North China Electric Power University, 2012.
[2]	Vermeersch F. ALARA pre-job studies using the VISIPLAN 3D ALARA planning tool [J]. Radiation Protection Dosimetry[J]. 2005, 115(1-4): 294-297.
[3]	Jeong K S, Choi B S, Moon J K, et al. Real-time assessment of exposure dose to workers in radiolohical environments during decommissionig of nuclear facilities[J]. Annals of Nuclear Energy, 2014, 73: 441-445.
[4]	Prokhorets I M. Point-kernel method for radiation fields simulation [J]. Problems of Atomic Science and Technology, 2007, N5: 106-109.
[5]	陈宜正, 王鑫, 武帧, 等.用于迷宫屏蔽计算的蒙特卡罗减方差方法研究 [J]. 原子能科学技术, 2018, 52(3): 562-569.CHEN Yizheng, WANG Xin, WU Zhen, et al. Study on Monte Carlo variance-reduction method used for maze shielding calculation[J]. Atomic Energy Science and Technology, 2018, 52(3): 562-569.
[6]	蔡元菲. 快速三维克里金插值方法研究及实现[D]. 电子科技大学,2013.CAI Yuanfei. The research and implementation of the fast three-dimentional kriging interpolation method[D]. University of Electronic Science and Technology of China, 2013.
[7]	Buhmann M D. Radial basis functions: theory and implementations[M]. Cambridge University Press, 2003.
[8]	张文波, 曹耀钦, 孙伟. 基于有限元插值方法的电磁覆盖范围重建[J]. 计算机仿真, 2013, 30(7): 170-173.ZHANG Wenbo, CAO Yaoqin, SUN Wei. Electromagnetic wave coverage area reconstruction based on finite element interpolation[J]. Computer Simulation, 2013, 30(7): 170-173.
[9]	王家华, 高海余, 周叶. 克里金地质绘图技术——计算机的模型和算法[M]. 北京: 石油工业出版社,1999.WANG Jiahua, GAO Haiyu, ZHOU Ye. Kriging geological mapping technology—Computer models and algorithms[M]. Beijing: Petroleum industry press, 1999.
[10]	Franke R. Scattered data interpolation: test of some methods[J]. Mathematics of Computation, 1982, 38(157): 181-200.
[11]	王彬彬, 石丽红, 卢月明. 一种顾及高程的时空反距离加权插值方法[J]. 测绘与空间地理信息, 2018, 41(10): 82-85.WANG Binbin, SHI Lihong, LU Yueming. A space-time reverse distance weighted interpolation method considering elevation[J]. Geomatics & Spatial Information Technology, 2018, 41(10): 82-85.
[12]	卢月明, 王亮, 仇阿根, 等. 局部加权线性回归模型的PM2.5空间插值方法[J]. 测绘科学, 2018, 43(11): 79-84.LU Yueming, WANG Liang, QIU Agen, et al. PM2.5 spatial interpolation method based on local weighted linear regression model[J]. Science of Surveying and Mapping, 2018, 43(11): 79-84.
[13]	顾军华, 许鹏, 董瑶, 等. 基于克里金插值的自适应VIRE室内定位算法研究[J]. 计算机工程与应用, 2018, 54(12): 100-104.GU Junhua, XU Peng, DONG Yao, et al. Research on adaptive VIRE indoor positioning algorithm based on Kriging interpolation[J]. Computer Engineering and Applications, 2018, 54(12): 100-104.
[14]	Lee S, Wolberg G, Shin S Y. Scattered data interpolation with multilevel b-splines[J]. IEEE Trans on Visualization and Computer Graphics, 1997, 3(3): 228-244.
[15]	Lophaven S N, Nielsrn H B, Sandergaard J. DACE-A matlab kriging toolbox, Version2.0[R]. Denmark: Technical University of Denmark, 2002.