湯凌志 尹陳艷 陽小華 劉紫靜

Preliminary Study of Radiation Field Reconstruction Based on Radial Basis Function Interpolation Method

TANG Ling-zhi YANG Xiao-hua;

LIU Zi-jing YIN Chen-yan

摘要：目前，我國有許多早期核設(shè)施面臨退役，輻射場的重構(gòu)技術(shù)研究對(duì)于輻射防護(hù)具有重要意義。本文結(jié)合我國早期核設(shè)施退役工程的特點(diǎn)，介紹和探討了基于徑向基函數(shù)插值方法對(duì)退役核設(shè)施進(jìn)行輻射場重構(gòu)的研究現(xiàn)狀和展望。

Abstract： At present， many early nuclear facilities in China are facing decommissioning， and the research on reconstruction technology of radiation fields is of great significance for radiation protection. Based on the characteristics of China's early nuclear facility decommissioning project， this paper introduces and discusses the research status and prospects of radial field reconstruction based on radial basis function interpolation method for decommissioning nuclear facilities.

關(guān)鍵詞：徑向基函數(shù)插值方法;復(fù)雜輻射場重構(gòu);核設(shè)施退役

Key words： radial basis function interpolation method;radiation field reconstruction;nuclear facilities decommissioning

中圖分類號(hào)：TP393.01? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)識(shí)碼：A? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文章編號(hào)：1006-4311（2019）32-0247-02

0? 引言

經(jīng)過60多年的發(fā)展，我國的核工業(yè)已經(jīng)形成了一套完備核工業(yè)體系。目前，為了避免因?yàn)楸O(jiān)管不嚴(yán)或者處置不當(dāng)而使得遺留的放射性物質(zhì)威脅環(huán)境和公眾健康，有許多超期服役、設(shè)施老化的核設(shè)施面臨退役[1]。我國于1990年啟動(dòng)了一批核設(shè)施退役及放射性廢物治理的項(xiàng)目。這些項(xiàng)目主要針對(duì)我國早期的核設(shè)施（如404廠、821廠）的退役以及歷史遺留的放射性廢物的治理[1]。2017年3月，國務(wù)院批復(fù)的《核安全與放射性污染防治“十三五”規(guī)劃及2025年遠(yuǎn)景目標(biāo)》[2]中明確提出了“早期核設(shè)施退役取得明顯成效”的規(guī)劃目標(biāo)。把“統(tǒng)籌推進(jìn)，加快早期核設(shè)施及放射性廢物處理處置”作為重點(diǎn)任務(wù)之一?！昂嗽O(shè)施退役及放射性廢物治理工程”被列為“十三五”重點(diǎn)工程之一。

將虛擬仿真技術(shù)應(yīng)用到核設(shè)施退役工程中，可以將輻射場的分布以可視化的方式顯示在三維虛擬場景中，能夠在實(shí)際拆除之前評(píng)估施工人員所接受的劑量，從而為核設(shè)施退役方案的優(yōu)化與決策提供科學(xué)依據(jù)，以減少人員受到不必要的輻射[3]。而這一切工作的前提是必須在核設(shè)施退役虛擬仿真系統(tǒng)中建立人員作業(yè)區(qū)域內(nèi)的復(fù)雜輻射場。

1? 研究意義

我國早期核設(shè)施退役工程的特點(diǎn)：年代久遠(yuǎn)，設(shè)計(jì)圖紙、源項(xiàng)等數(shù)據(jù)不明、源項(xiàng)分布及場景結(jié)構(gòu)復(fù)雜。而傳統(tǒng)的正演計(jì)算輻射場分布的蒙特卡羅方法和點(diǎn)核積分方法都需要利用源項(xiàng)數(shù)據(jù)和場景結(jié)構(gòu)信息編寫劑量計(jì)算使用的輸入卡，且蒙特卡羅方法對(duì)大空間厚屏蔽問題無法在合理的時(shí)間內(nèi)給出可靠的結(jié)果，而點(diǎn)核方法又不適合復(fù)雜源項(xiàng)和幾何區(qū)域，累積因子估算會(huì)帶來較大誤差[4]。在核設(shè)施退役源項(xiàng)調(diào)查獲得的數(shù)據(jù)的基礎(chǔ)上，可以利用徑向基函數(shù)插值方法對(duì)實(shí)測產(chǎn)生的離散的輻射場測量數(shù)據(jù)進(jìn)行空間重構(gòu)和分析，進(jìn)而反演出完整的輻射場的分布。

因此，在已有相關(guān)研究的基礎(chǔ)上，利用多學(xué)科交叉技術(shù)，開展復(fù)雜輻射場重構(gòu)方法研究，可以解決我國早期核設(shè)施退役工程的退役虛擬仿真系統(tǒng)的關(guān)鍵技術(shù)問題，為進(jìn)一步開展退役仿真技術(shù)研究和工程應(yīng)用系統(tǒng)的開發(fā)奠定基礎(chǔ)，在核設(shè)施運(yùn)行、檢修、退役等輻射領(lǐng)域，有廣闊的應(yīng)用前景[5]和較高的實(shí)用價(jià)值。

2? 國內(nèi)外研究現(xiàn)狀

在歐、美、日等發(fā)達(dá)國家，核電發(fā)展相對(duì)比較早，在核與輻射安全領(lǐng)域研究并開發(fā)了相應(yīng)的軟件系統(tǒng)[6]，代表性的有法國的NARVEOS系統(tǒng)、比利時(shí)核能研究中心的VISIPLAN、西班牙的SICOMORO、日本的DEXUS和VRdose等。這些核設(shè)施虛擬仿真軟件系統(tǒng)主要包括作業(yè)任務(wù)的模擬、工作人員受照射劑量的估算，退役方案的制定及優(yōu)化等方面，都是應(yīng)用于核電廠，其輻射場分布計(jì)算方法大部分使用的是點(diǎn)核積分法。

目前，國內(nèi)僅少量科研院所和大學(xué)在開展輻射作業(yè)劑量虛擬仿真方面軟件的研發(fā)工作。中國科學(xué)院核能安全技術(shù)研究所 FDS 團(tuán)隊(duì)開發(fā)的Super MC/RVIS。其核心計(jì)算所使用的是蒙特卡羅方法SuperMC。中國核動(dòng)力研究設(shè)計(jì)院張永領(lǐng)等人進(jìn)行了反應(yīng)堆退役三維仿真技術(shù)的研究，研究出了基于點(diǎn)核積分算法的三維輻射場計(jì)算算法，開發(fā)了基于DELMIA和VIRTOOLS平臺(tái)的反應(yīng)堆退役三維仿真原型系統(tǒng)。南華大學(xué)宋英明、梁燁等人針對(duì)核設(shè)施退役過程中多源項(xiàng)的三維輻射場，采用點(diǎn)核積分的方法進(jìn)行重構(gòu)，設(shè)計(jì)開發(fā)了面向核設(shè)施退役過程的輻射場重構(gòu)與拆除路徑優(yōu)化功能軟件[7]。

對(duì)于復(fù)雜幾何空間中的輻射場計(jì)算，目前常用的方法有蒙特卡羅方法、離散縱標(biāo)法和點(diǎn)核積分方法，這些方法實(shí)際上是對(duì)輻射傳輸方程（Boltzmann方程）的求解，可歸于“正演”方法。然而，要對(duì)輻射傳輸方程進(jìn)行求解計(jì)算，不論是用確定論方法還是隨機(jī)性方法，都需要在清楚地了解放射源系統(tǒng)的構(gòu)造信息的前提下，方可構(gòu)建準(zhǔn)確的系統(tǒng)模型進(jìn)行粒子輸運(yùn)模擬，對(duì)于我國早期核設(shè)施退役工程這種源項(xiàng)數(shù)據(jù)不清楚、源項(xiàng)分布及場景結(jié)構(gòu)復(fù)雜的情況不適用[8][9]。

除了通過計(jì)算獲得輻射場的分布以外，在核設(shè)施退役工程的實(shí)際應(yīng)用中，還可以通過在作業(yè)區(qū)域用實(shí)驗(yàn)測量的方式采集輻射場數(shù)據(jù)。但是由于條件限制，不可能測量整個(gè)輻射場的空間分布?？梢酝ㄟ^插值算法對(duì)源項(xiàng)調(diào)查實(shí)測獲得的有限的、離散的采樣數(shù)據(jù)進(jìn)行分析和空間重構(gòu)，從而獲得完整的輻射場分布數(shù)據(jù)。在從輻射場測量數(shù)據(jù)到輻射場數(shù)據(jù)重構(gòu)（反演）這個(gè)研究方向上，僅有國內(nèi)少數(shù)研究者提出了自己的方法。

華南理工大學(xué)的蔡杰進(jìn)、王壯等人[10][11]提出了一種基于網(wǎng)格插值的三維輻射場劑量的反演方法。但該方法需要預(yù)先測量矩形網(wǎng)格上較多點(diǎn)的輻射劑量值，而在早期核設(shè)施退役實(shí)際應(yīng)用中會(huì)出現(xiàn)探測器因復(fù)雜幾何條件無法布置的情況。

哈爾濱工程大學(xué)的劉永闊等人[12]提出了一種核設(shè)施退役輻射場劑量分布仿真方法。該方法根據(jù)樣本數(shù)據(jù)，先構(gòu)建徑向基神經(jīng)網(wǎng)絡(luò)模型，然后通過反距離權(quán)重法計(jì)算核設(shè)施退役劑量場中任意一點(diǎn)的劑量值，最后計(jì)算出輻射場的劑量分布。反距離權(quán)重法具有計(jì)算快，實(shí)現(xiàn)方便的優(yōu)點(diǎn)，但其局限比較多，性能受采樣點(diǎn)分布影響較大，只有對(duì)于均勻分布的數(shù)據(jù)，得到的結(jié)果較好[13]。

中國工程物理研究院賽雪、韋孟伏的文獻(xiàn)[8][9]中，基于MQ徑向基函數(shù)散亂插值方法，利用少量、離散、不規(guī)則分布的模擬數(shù)據(jù)和實(shí)測數(shù)據(jù)，對(duì)γ輻射場進(jìn)行了重構(gòu)及可視化。證明了將MQ方法應(yīng)用于γ輻射場數(shù)據(jù)重構(gòu)是可行的。

散亂數(shù)據(jù)插值方法包括一系列插值或逼近的算法。目前其在測繪學(xué)、地形學(xué)、地球物理以及產(chǎn)品外形設(shè)計(jì)等領(lǐng)域廣泛應(yīng)用。反距離權(quán)重插值法、有限元法以及徑向基函數(shù)插值法都屬于目前最常用的中、小規(guī)模散亂數(shù)據(jù)插值算法[14]。

徑向基函數(shù)是一類以徑向距離為變量的基函數(shù)集合。徑向基函數(shù)插值表達(dá)形式簡潔，計(jì)算精度優(yōu)異、滿足適定性要求，與維度無關(guān)，是一種精確空間插值方法[15]。

徑向基函數(shù)法同時(shí)也屬于無網(wǎng)格方法。20世紀(jì)90年代以來，這種僅僅基于節(jié)點(diǎn)而不需劃分單元的數(shù)值方法取得了很大的發(fā)展，目前已成為計(jì)算科學(xué)研究的熱點(diǎn)之一。

R.L. Hardy在1971年提出的MQ插值方法[16]是目前應(yīng)用最為成功的徑向基函數(shù)插值算法之一。目前其在測繪學(xué)、水力學(xué)、攝影制圖法、地球物理等方面均有廣泛應(yīng)用。

3? 結(jié)論與展望

賽雪等人的文獻(xiàn)[8][9]以兩種輻射場為例，基于少量采樣點(diǎn)，實(shí)現(xiàn)了輻射場的重構(gòu)及可視化，探究了MQ徑向基函數(shù)插值方法在γ輻射場重構(gòu)的可行性。但是其采用的物理模型還比較簡單，今后的研究中可研究、嘗試將該方法應(yīng)用于分布更復(fù)雜的退役核設(shè)施的輻射場重構(gòu)與分析中;可與神經(jīng)網(wǎng)絡(luò)算法結(jié)合，使輻射場重構(gòu)的結(jié)果逼近精度更高，且具有較好的抗噪、平滑能力;可采用逐點(diǎn)交叉驗(yàn)證方法，能有效獲取插值模型的最佳形態(tài)參數(shù)，提高結(jié)果的精度和可靠性;可與輻射防護(hù)算法相結(jié)合，提出符合退役核設(shè)施輻射場分布特點(diǎn)的基函數(shù)，從而在輻射場重構(gòu)的精度和計(jì)算運(yùn)行的時(shí)間效率上都滿足工程實(shí)踐的需求。

參考文獻(xiàn)：

[1]劉華. 核與輻射設(shè)施退役監(jiān)管[A].中國核學(xué)會(huì)輻射防護(hù)分會(huì).全國核與輻射設(shè)施退役學(xué)術(shù)研討會(huì)論文集[C].中國核學(xué)會(huì)輻射防護(hù)分會(huì)：中國核學(xué)會(huì)，2007：1.

[2]中華人民共和國環(huán)境保護(hù)部.國務(wù)院批復(fù)《核安全與放射性污染防治“十三五”規(guī)劃及2025年遠(yuǎn)景目標(biāo)》[EB/OL].http：//hssaq.mee.gov.cn/dtxx/201703/t20170323_408666.shtml，2017-03-23/2018-11-28.

[3]王俊峰.核設(shè)施退役工程[M].中國原子能出版社，2013.

[4]畢遠(yuǎn)杰，楊宏偉，郭金森，駱志平，李傳龍，陳凌.研究堆退役中的快速輻射傳輸計(jì)算方法研究[J].原子能科學(xué)技術(shù)，2015，49（S1）：470-474.

[5]劉中坤，彭敏俊，朱海山，成守宇，鞏誠.核設(shè)施退役虛擬仿真系統(tǒng)框架研究[J].原子能科學(xué)技術(shù)，2011，45（09）：1080-1086.

[6]李華.γ輻射場快速計(jì)算與源項(xiàng)反演算法研究及初步應(yīng)用[D].清華大學(xué)，2016.

[7]宋英明，梁燁，葉凱萱，張震宇，張秋楠，王寧，顏佳偉，朱志超，鄒樹梁.核設(shè)施退役過程中的輻射場重構(gòu)與拆除路徑優(yōu)化[J].核技術(shù)，2017，40（05）：54-60.

[8]賽雪.基于散亂數(shù)據(jù)插值方法的γ輻射場可視化技術(shù)研究[D].中國工程物理研究院，2017.

[9]賽雪，陳穎，韋孟伏.Multiquadric散亂數(shù)據(jù)插值方法在γ輻射場可視化中的應(yīng)用初探[J].核技術(shù)，2016（10）：55-61.

[10]蔡杰進(jìn)，王壯，劉榮.一種三維輻射場劑量的反演方法[P]. 中國專利：CN107797132A，2018-03-13.

[11]Zhuang Wang， Jiejin Cai，Inversion of radiation field on nuclear facilities A method based on net function interpolation，Radiation Physics and Chemistry，Volume 153，2018，Pages 27-34.

[12]劉永闊，李夢堃，彭敏俊，楊立群.一種核設(shè)施退役輻射場劑量分布仿真方法 [P].中國專利：CN107330187A，2017-11-07.

[13]陳堃.空間散亂分布位場數(shù)據(jù)插值研究[D].浙江大學(xué)，2013.

[14]殷浩，戴光明.散亂數(shù)據(jù)可視化研究綜述[J].微機(jī)發(fā)展，2005（07）：7-10，13.

[15]Buhmann， Martin D . Radial basis functions：theory and implementations[M]. Cambridge University Press， 2003.

[16]Hardy R L . Multiquadric equations of topography and other irregular surfaces[J]. Journal of Geophysical Resarch， 1971， 76（8）：1905-1915.

基金項(xiàng)目：湖南省教育廳科學(xué)研究項(xiàng)目（課題號(hào)：11C1082）。

作者簡介：湯凌志（1976-），男，湖南岳陽人，實(shí)驗(yàn)師，研究生，研究方向?yàn)楹藬?shù)據(jù)處理與計(jì)算。