蒙卡計(jì)算和神經(jīng)插值法構(gòu)建退役核設(shè)施輻射場(chǎng)的應(yīng)用機(jī)理

中核四川環(huán)保工程有限責(zé)任公司 張睿超 詹戈輝 張維安

核設(shè)施退役是一項(xiàng)復(fù)雜的工程，其重要工作之一是對(duì)核設(shè)施內(nèi)輻射場(chǎng)情況的掌握，以便根據(jù)設(shè)施內(nèi)輻射場(chǎng)特點(diǎn)設(shè)計(jì)不同的退役處置方案和指導(dǎo)現(xiàn)場(chǎng)人員進(jìn)入實(shí)施。本文圍繞數(shù)字化核退役仿真技術(shù)，基于蒙特卡羅和神經(jīng)差值法提出一種輻射場(chǎng)估算方法，并利用Unity 引擎構(gòu)建輻射場(chǎng)可視化熱力圖，實(shí)現(xiàn)對(duì)退役核設(shè)施內(nèi)部輻射場(chǎng)的仿真。從實(shí)踐效果看，該項(xiàng)研究中提出的輻射場(chǎng)仿真具有一定實(shí)用性。該項(xiàng)研究也對(duì)后續(xù)數(shù)字化核退役的探索提供了數(shù)據(jù)支持。

1 背景

1.1 核設(shè)施退役介紹

當(dāng)前，國(guó)內(nèi)外有相當(dāng)數(shù)量的核設(shè)施處于退役或準(zhǔn)備退役。核設(shè)施退役工作包含了設(shè)計(jì)、實(shí)施、安全管理、退役評(píng)估等環(huán)節(jié)，在設(shè)計(jì)階段的一項(xiàng)重要工作即是掌握設(shè)施內(nèi)輻射場(chǎng)情況，以便根據(jù)設(shè)施內(nèi)輻射場(chǎng)特點(diǎn)設(shè)計(jì)不同的退役處置方案和指導(dǎo)現(xiàn)場(chǎng)人員進(jìn)入實(shí)施[1]。以某單位的廢液存儲(chǔ)廠房為例，該核設(shè)施建造于20世紀(jì)60年代，起初用于貯存該單位其他核設(shè)施運(yùn)行、維護(hù)活動(dòng)中所產(chǎn)生的低放射性液體廢物，現(xiàn)已經(jīng)停止使用，進(jìn)入退役設(shè)計(jì)和申報(bào)實(shí)施階段。在退役設(shè)計(jì)過(guò)程中，對(duì)該設(shè)施選取若干點(diǎn)進(jìn)行輻射場(chǎng)測(cè)量，形成初步輻射場(chǎng)數(shù)據(jù)供參考，但直接使用數(shù)據(jù)點(diǎn)缺乏直觀性，因此有必要將輻射場(chǎng)可視化，便于設(shè)計(jì)和實(shí)施人員使用。

同時(shí)，為了貫徹落實(shí)國(guó)家對(duì)退役治理、環(huán)境保護(hù)、生態(tài)文明建設(shè)的指示批示精神，提高政治站位，增強(qiáng)危機(jī)意識(shí)和極限思維，加快推進(jìn)重大安全隱患消除，堅(jiān)定不移踐行總體國(guó)家安全觀，推進(jìn)“理性、協(xié)調(diào)、并進(jìn)”的核安全觀，做到“充分準(zhǔn)備，一絲不茍，萬(wàn)無(wú)一失，一次成功”，全面推進(jìn)信息化與業(yè)務(wù)工作深度融合。在核設(shè)施退役領(lǐng)域，使用數(shù)字化技術(shù)構(gòu)建輻射場(chǎng)可視化，也是加快退役治理、提高核安全管理能力的必要手段。

1.2 基本模型介紹

（1）蒙特卡羅方法：又稱(chēng)隨機(jī)抽樣方法或統(tǒng)計(jì)模擬方法，屬于試驗(yàn)數(shù)學(xué)的一個(gè)分支。該方法是以概率統(tǒng)計(jì)理論為基礎(chǔ)的計(jì)算方法，利用隨機(jī)數(shù)進(jìn)行統(tǒng)計(jì)試驗(yàn)，以求得的統(tǒng)計(jì)特征值(如均值、概率等)作為待解問(wèn)題的數(shù)值解。在信息技術(shù)不斷發(fā)展過(guò)后，蒙特卡羅方法已被廣泛應(yīng)用于核物理、工程等各個(gè)方面。例如在研究中子、光子和電子等粒子輸運(yùn)的領(lǐng)域，就有Los Alamos 國(guó)家實(shí)驗(yàn)室應(yīng)用理論物理部研發(fā)的Monte Carlo Neutron and Photo Transport Code (蒙特卡羅中子-光子輸運(yùn)程序，MCNP)，利用計(jì)算機(jī)程序?qū)崿F(xiàn)蒙特卡羅方法模擬單粒子的徑跡，追蹤粒子全生命周期運(yùn)行流程，通過(guò)重復(fù)實(shí)驗(yàn)推導(dǎo)出粒子輸運(yùn)過(guò)程模擬[2]。

（2）神經(jīng)插值：神經(jīng)插值即是用神經(jīng)網(wǎng)絡(luò)技術(shù)進(jìn)行插值計(jì)算。神經(jīng)網(wǎng)絡(luò)一般是由大量的簡(jiǎn)單計(jì)算單元（即神經(jīng)元）構(gòu)成的非線性系統(tǒng)，該系統(tǒng)可一定程度上對(duì)人腦神經(jīng)系統(tǒng)進(jìn)行模擬，實(shí)現(xiàn)信息處理、存儲(chǔ)及檢索功能，并以此解決復(fù)雜計(jì)算問(wèn)題。在各類(lèi)神經(jīng)網(wǎng)絡(luò)中，使用徑向基函數(shù)（Radial Basis Function Network，RBF）作為激活函數(shù)構(gòu)造的神經(jīng)網(wǎng)絡(luò)被廣泛用于序列預(yù)測(cè)、系統(tǒng)控制和插值計(jì)算等領(lǐng)域。RBF 神經(jīng)網(wǎng)絡(luò)構(gòu)造的基本思想，即是用RBF 作為隱單元的“基”構(gòu)成隱含層空間，以此將輸入矢量直接映射到隱含層空間，而不需要通過(guò)權(quán)連接。當(dāng)RBF 的中心點(diǎn)確定以后，這種映射關(guān)系也就確定了。而隱含層空間到輸出空間的映射是線性的，即網(wǎng)絡(luò)的輸出是隱單元輸出的線性加權(quán)和，此處的權(quán)即為網(wǎng)絡(luò)可調(diào)參數(shù)?；赗BF 神經(jīng)網(wǎng)絡(luò)的插值法具有局部逼近的優(yōu)點(diǎn)，因?yàn)镽BF 神經(jīng)網(wǎng)絡(luò)是一種性能優(yōu)良的前饋型神經(jīng)網(wǎng)絡(luò)，可以任意精度逼近任意的非線性函數(shù)，且具有全局逼近能力，從根本上解決了經(jīng)典神經(jīng)網(wǎng)絡(luò)（如BP 網(wǎng)絡(luò)）的局部最優(yōu)問(wèn)題；而且拓?fù)浣Y(jié)構(gòu)緊湊，結(jié)構(gòu)參數(shù)可實(shí)現(xiàn)分離學(xué)習(xí)，收斂速度快。RBF 網(wǎng)絡(luò)和模糊邏輯能夠?qū)崿F(xiàn)很好的互補(bǔ)，提高神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)泛化能力。RBF網(wǎng)絡(luò)能夠處理系統(tǒng)內(nèi)的難以解析的規(guī)律性，具有良好的泛化能力，并有很快的學(xué)習(xí)收斂速度，已成功應(yīng)用于非線性函數(shù)逼近、時(shí)間序列分析、數(shù)據(jù)分類(lèi)、模式識(shí)別、信息處理、圖像處理、系統(tǒng)建模、控制和故障診斷等，是神經(jīng)網(wǎng)絡(luò)插值計(jì)算的首選方法[3]。

2 探索與應(yīng)用

2.1 基于MCNP5 軟件的蒙特卡羅計(jì)算

MCNP5 軟件的系統(tǒng)架構(gòu)如圖1所示。

圖1 MCNP5 軟件架構(gòu)Fig.1 Software structs of MCNP5

初始輸入文件的結(jié)構(gòu)，對(duì)目標(biāo)的幾何結(jié)構(gòu)、材料、記數(shù)要求等給以描述（如表1所示）。

表1 MCNP 輸入數(shù)據(jù)格式Tab.1 Input data format of MCNP

在MCNP 使用中，源項(xiàng)和輻射粒子的種類(lèi)是通過(guò)SDEF 命令準(zhǔn)確指定來(lái)的，通過(guò)修改SDEF 參數(shù)可以定義目標(biāo)中源項(xiàng)的各類(lèi)特征。輸入文件中只允許唯一的SDEF 卡。

關(guān)于材料卡，MCNP 運(yùn)算中填充各個(gè)柵元的材料說(shuō)明包含以下部分：（1）定義一個(gè)唯一的材料編號(hào)；（2）元素（或核素）組成；（3）所使用的截面庫(kù)（如表2所示）。

表2 MCNP 材料卡Tab.2 MCNP material cards

計(jì)算模型的建立，包括核素名稱(chēng)、輻射類(lèi)型（光子或中子）、活度濃度、輻射能譜、源項(xiàng)空間分布類(lèi)型、分布參數(shù)等源項(xiàng)數(shù)據(jù)。由于廠區(qū)內(nèi)地下貯罐實(shí)際儲(chǔ)存情況比較復(fù)雜，由于長(zhǎng)時(shí)間貯藏的泥漿廢液多次蒸發(fā)沉淀，部分貯罐形成上層清液下層泥漿的狀態(tài)，而部分貯罐由于長(zhǎng)時(shí)間未使用，已經(jīng)結(jié)板固化。精確性層面上無(wú)法用準(zhǔn)確的質(zhì)量密度活度來(lái)進(jìn)行數(shù)字模型描述。由此客觀條件引起的誤差，同時(shí)輔以源項(xiàng)調(diào)查監(jiān)測(cè)出來(lái)的N 個(gè)定點(diǎn)數(shù)據(jù)為樣本進(jìn)行徑向基函數(shù)(RBF)插值，進(jìn)行大區(qū)域復(fù)雜環(huán)境的輻射場(chǎng)量補(bǔ)償模擬計(jì)算。

2.2 徑向基函數(shù)(RBF)插值補(bǔ)償

（1）模型假設(shè)：設(shè)定地面網(wǎng)格是等距，設(shè)置單位網(wǎng)格面積為1m×1m；鑒于科研嘗試，所涉及的廠區(qū)范圍不大，將該科研區(qū)域暫定為一個(gè)平面區(qū)域；科研區(qū)域的網(wǎng)格點(diǎn)排序，需要和數(shù)據(jù)文件中的數(shù)據(jù)坐標(biāo)精確映射；源項(xiàng)調(diào)查數(shù)據(jù)按照生產(chǎn)管理經(jīng)驗(yàn)必須覆蓋高危險(xiǎn)區(qū)域數(shù)據(jù)，假設(shè)其他未檢測(cè)點(diǎn)為連續(xù)衰減點(diǎn)位數(shù)據(jù)。

估算的準(zhǔn)確率和誤差水平受廠區(qū)現(xiàn)場(chǎng)污染水平程度、以及源項(xiàng)調(diào)查測(cè)量能力、相關(guān)算法的影響較大，以某年進(jìn)行的源項(xiàng)調(diào)查報(bào)告監(jiān)測(cè)出來(lái)的N 個(gè)定點(diǎn)數(shù)據(jù)進(jìn)行徑向基函數(shù)(RBF)插值。數(shù)據(jù)格式如表3所示。

表3 輻射場(chǎng)調(diào)查原始數(shù)據(jù)格式Tab.3 Raw data format of radiation survey

首先還是導(dǎo)入所需要的庫(kù)或模塊（如表4所示）。

表4 導(dǎo)入模塊Tab.4 Import modules

（2）主要參考算法模型：我們建立了一個(gè)RBF 的類(lèi)函數(shù)，之后隨機(jī)選取了RBF 的計(jì)算中心，分別計(jì)算輸入值和RBF 中心的距離，訓(xùn)練之后的數(shù)組用之后的學(xué)習(xí)數(shù)據(jù)再次帶入測(cè)試計(jì)算準(zhǔn)確度（如表5所示）。

表5 算法實(shí)現(xiàn)Tab.5 Algorithm implementation

（3）中心：這里trainFunc 是使用K-means 方法，就是傳統(tǒng)的算法，先隨機(jī)選k 個(gè)樣本作為中心，然后按照歐氏距離對(duì)每個(gè)樣本分組，再重新確定聚類(lèi)中心，再不斷重復(fù)上面的步驟，直到最終聚類(lèi)中心變化在一定范圍內(nèi)。

（4）寬度：寬度向量影響著神經(jīng)元對(duì)輸入信息的作用范圍，寬度越小，相應(yīng)隱含層神經(jīng)元作用函數(shù)的形狀越窄，那么處于其他神經(jīng)元中心附近的信息在該神經(jīng)元出的響應(yīng)就越??；就跟高斯函數(shù)圖像兩邊的上升下降區(qū)域的寬度一樣。

（5）訓(xùn)練步驟：先初始化中心、寬度、最后一層權(quán)重；計(jì)算損失，如果在可以接受的范圍內(nèi)，停止訓(xùn)練；利用梯度更新的方法更新中心、寬度、權(quán)重；返回前一步計(jì)算損失。

（6）應(yīng)用限制：神經(jīng)網(wǎng)絡(luò)最限制實(shí)施的問(wèn)題是沒(méi)有業(yè)務(wù)相關(guān)理論來(lái)解釋自己的推理過(guò)程和推理依據(jù)；當(dāng)提供的原始數(shù)據(jù)不充分的情況下，神經(jīng)網(wǎng)絡(luò)模型就無(wú)法進(jìn)行工作；隱層基函數(shù)的中心是在輸入樣本集中選取的，在很多情況下難以反映系統(tǒng)真正的輸入輸出關(guān)系，并且初始中心點(diǎn)數(shù)太多，優(yōu)選過(guò)程會(huì)異常困難。

2.3 Unity3D 中對(duì)數(shù)據(jù)進(jìn)行可視化

針對(duì)輻射場(chǎng)的可視化，可以采用三維建模疊加輻射場(chǎng)數(shù)據(jù)顯示的方式實(shí)現(xiàn)[4]?；赨nity3D 的工業(yè)可視化的項(xiàng)目中常見(jiàn)的需求之一，就是顯示某個(gè)設(shè)備的計(jì)量信息，或者是顯示一間房間的計(jì)量信息。在退役治理領(lǐng)域，常見(jiàn)的就是空間輻射場(chǎng)可視化。在得到計(jì)算后的輻射場(chǎng)數(shù)據(jù)后，使用Unity3D 平臺(tái)及相關(guān)組件繪制熱力圖（如表6所示）。

表6 熱力圖主邏輯實(shí)現(xiàn)Tab.6 Implementation of heatmap

（1）首先從上述預(yù)估方法輸出文檔中解析數(shù)據(jù)，獲取需要的數(shù)據(jù)。就是說(shuō)有n×n 個(gè)數(shù)據(jù)，可以用二維數(shù)組去接收數(shù)據(jù)，然后用一個(gè)List 數(shù)組去接收解析完的數(shù)據(jù)。

（2）利用坐標(biāo)點(diǎn)中獲取的α、β、γ 值進(jìn)行取相應(yīng)的顏色值，將貼圖數(shù)據(jù)賦值給Image 的Sprite。

3 結(jié)語(yǔ)

在數(shù)字化核退役仿真項(xiàng)目中，經(jīng)過(guò)有限源項(xiàng)調(diào)查數(shù)據(jù)繪制退役廠區(qū)的輻射場(chǎng)熱力圖，總結(jié)得出在復(fù)雜多源的建構(gòu)物中進(jìn)行輻射場(chǎng)估算繪制，利用單一的算法無(wú)法達(dá)到科研預(yù)期效果，針對(duì)這種場(chǎng)景條件下的作業(yè)，多種算法耦合相互彌補(bǔ)缺陷是當(dāng)前解決方案的最優(yōu)解，也是一種行之有效的辦法。從本研究結(jié)果看，由于硬件層面的限制，輻射場(chǎng)估算的計(jì)算用時(shí)還很受影響，后續(xù)要實(shí)現(xiàn)動(dòng)態(tài)生成輻射場(chǎng)估算，可以通過(guò)優(yōu)化算法的運(yùn)行時(shí)間提高效率。本研究可作為核設(shè)施退役工程設(shè)計(jì)的參考[5]，為最優(yōu)化人員行動(dòng)路線規(guī)劃、退役過(guò)程人員輻射劑量評(píng)估等提供數(shù)據(jù)支持依據(jù)。