龍澤宇, 宋佩濤, 馮青靚, 李華, 韋加富, 劉瑞文, ?，挷?/p>

(1. 中國輻射防護研究院 輻射安全與防護山西省重點實驗室, 太原, 030006;2. 中國輻射防護研究院 國家環(huán)境保護輻射環(huán)境與健康重點實驗室, 太原, 030006;3. 中國輻射防護研究院 核與輻射前沿技術(shù)研究中心, 太原, 030006)

0 引言

在核設(shè)施退役工作中,工作人員通常需要直接進入到高劑量照射的輻射區(qū)域中進行作業(yè),準確獲取區(qū)域內(nèi)的輻射場分布,不僅可以精準有效地指導工作方案的具體規(guī)劃,還能通過輻射劑量可視化等手段,消除工作人員對輻射“看不見,摸不著”的恐慌心理,對工作人員的輻射防護最優(yōu)化有現(xiàn)實意義。

通常,獲取輻射場的方法可分為基于輻射源項的正演計算方法和基于輻射測量數(shù)據(jù)的插值計算方法兩大類。正演計算是指根據(jù)放射源及場景信息,直接計算輻射場的劑量率分布[1];插值計算方法,則是在利用固定或移動式測量設(shè)備獲取的離散輻射測量數(shù)據(jù)基礎(chǔ)上,利用插值算法實現(xiàn)整個輻射區(qū)域劑量場的重構(gòu)。基于輻射測量數(shù)據(jù)的插值計算方法僅依賴于實測的輻射劑量率數(shù)據(jù),不需要掌握具體的輻射源項信息,適用性要優(yōu)于正演計算方法。

中國工程物理研究院賽雪利用MQ徑向基函數(shù)散亂數(shù)據(jù)插值法,實現(xiàn)了只有少量不規(guī)則分布的劑量數(shù)據(jù)下的輻射場重構(gòu)及可視化[2],但其所建立的模型還較為簡單,僅適用于兩種簡單情況下的輻射場重構(gòu)。哈爾濱工程大學的李夢堃提出了一種核設(shè)施退役輻射場劑量分布仿真方法[3],該方法根據(jù)樣本數(shù)據(jù),先構(gòu)建徑向基神經(jīng)網(wǎng)絡(luò)模型,然后通過反距離權(quán)重法計算核設(shè)施退役劑量場中任意一點的劑量值,最后計算出輻射場的劑量率分布,但實驗結(jié)果表明,只有對于均勻分布的數(shù)據(jù),插值結(jié)果才較好。南華大學張敏提出了一種利用輻射劑量值占優(yōu)構(gòu)造還原算法對已知探測點的劑量進行外插和內(nèi)推的方法[4],重新構(gòu)造出了完整的平面輻射場,但是重構(gòu)得到的輻射場的準確度主要取決于探測點的數(shù)量,當測點數(shù)量較少時,重構(gòu)的準確度較低。

地統(tǒng)計插值技術(shù)起初是從研究礦產(chǎn)儲量計算問題而產(chǎn)生并發(fā)展起來的,但凡是要研究空間分布數(shù)據(jù)的自相關(guān)性,并對這批數(shù)據(jù)進行最優(yōu)無偏插值估計時,均可應(yīng)用地統(tǒng)計插值。輻射場本質(zhì)上由是放射源發(fā)出的射線產(chǎn)生的,同樣具有空間自相關(guān)性,距離越近的點劑量率就越接近,因此將地統(tǒng)計插值用于輻射場重構(gòu)是非常合適的。地統(tǒng)計空間插值盡管都以克里金方法為核心,但其種類繁多,考慮到輻射場還具有變化梯度大的特點,本文嘗試將普通克里金方法與對數(shù)變換相結(jié)合,在離散測量點的基礎(chǔ)上進行輻射場重構(gòu),并利用實測輻射場對其進行驗證,為源項未知情況下的三維輻射場獲取提供一種新方法。

1 方法原理

1.1 普通克里金方法

克里金法的實質(zhì)是利用區(qū)域化變量的觀測點數(shù)據(jù),對插值點的數(shù)據(jù)進行無偏估計,其中應(yīng)用最為廣泛的為普通克里金方法[5-6]。在普通克里金方法中,對于一個插值點的數(shù)據(jù)采用觀測點數(shù)據(jù)的線性組合來估計:

(1)

式中,Z*(x0)為插值點數(shù)據(jù)的估計值;Z(xi)為觀測點的數(shù)據(jù);λi為權(quán)重系數(shù),表示各觀測點對插值點估計值的貢獻程度,其求解必須滿足2個條件:線性無偏和使估計方差最小,分別見式(2)和式(3):

E[Z*(x0)-Z(x0)]=

(2)

(3)

式中,X=Z*(x0)-Z(x0)。

對于式(2),由于普通克里金方法假設(shè)空間任意點都有同樣的期望,即E[Z(x)]=c,可得:

(4)

對于式(3),要使估計方差最小,根據(jù)拉格朗日乘數(shù)法,可得:

(5)

式中,μ為拉格朗日乘數(shù);ri,j為半方差函數(shù)(見1.2節(jié))。對式(4)和式(5)進行求解,獲得權(quán)重系數(shù)λi,代入式(1),即可獲得插值點數(shù)據(jù)的估計值。

1.2 半方差函數(shù)及對數(shù)變換

半方差函數(shù)ri,j反映了區(qū)域化變量的變異程度,其計算公式見式(6):

(6)

普通克里金方法假設(shè)ri,j與兩點之間的距離di,j存在著函數(shù)關(guān)系,因此只需計算出觀測數(shù)據(jù)集中各點間的距離di,j和半方差ri,j,再尋找一個半方差擬合模型來擬合出兩者的關(guān)系,即可根據(jù)插值點和觀測點間的距離,計算出式(5)中的ri,0。

常用半方差擬合模型包括Linear、Power、Gaussian、 Spherical、Exponential等5種,分別見式(7)、式(8)、式(9)、式(10)、式(11):

r(d)=sd+n

(7)

r(d)=sdω+n

(8)

(9)

(10)

(11)

用普通克里金方法進行插值估計時,一般要求數(shù)據(jù)服從正態(tài)分布。然而由于輻射劑量遵循平方反比定律,各點處的劑量率值會在遠離源項的方向迅速下降,因此整個輻射場的劑量率值分布通常呈正偏態(tài)分布,直接使用普通克里金法來進行插值,結(jié)果的偏差可能會較大。

對此,本文在插值前先對觀測點數(shù)據(jù)進行對數(shù)變換,以使其更接近正態(tài)分布:

Z′(xi)=lnZ(xi)

(12)

在對各觀測點數(shù)據(jù)進行對數(shù)變換,并利用普通克里金方法進行插值后,將各點插值結(jié)果按式(13)再進行一次逆變換,作為最后的插值點估計值:

Z*(x0)=eZ*′(x0)

(13)

式中,Z*′(x0)為直接使用經(jīng)對數(shù)變換處理的觀測數(shù)據(jù)獲得的插值點估計值;Z*(x0)為經(jīng)過逆變換后最終的插值點估計值。

需要說明的是,考慮到輻射場具有各向同性的特點,因此本文不對各軸向進行拉伸處理。另外當測量點數(shù)量較多,且分布范圍較大時,為了避免離插值點距離過遠的觀測點對插值結(jié)果產(chǎn)生不利影響,可只采用離插值點距離最近的部分觀測點參與計算。

2 隨機選點情況下的測試驗證

2.1 測試驗證一

為驗證方法在三維輻射場重構(gòu)上的有效性和可靠性,在某廠房中選取了一放置有廢物桶的房間進行現(xiàn)場測量實驗。使用全站儀以及長桿式劑量率儀對廢物桶周圍的輻射場進行測量,劑量率儀測量輻射場中某一點的劑量率值,全站儀記錄該點對應(yīng)的坐標,使坐標與劑量率值一一對應(yīng),得到相應(yīng)的輻射場數(shù)據(jù),總數(shù)為150個點,其三維建模如圖1所示。

圖1 測試驗證一的三維建模圖

對所有的已知點進行隨機選取,選取數(shù)量分別為總數(shù)的15%(23個點)、20%(30個點)和25%(38個點),重復(fù)次數(shù)為5 000次,將選取到的點作為測量點,剩余的點作為輻射場重構(gòu)的驗證點,在不同半方差擬合模型下重構(gòu)結(jié)果的平均偏差分別列于表1、表2、表3(此處的平均偏差為單次隨機選點在各驗證點處重構(gòu)結(jié)果偏差的平均值,而表中的最大值、最小值、平均值為5 000次隨機選點的平均偏差中的最大值、最小值、平均值,下同)。

表1 測試驗證一中測量點數(shù)量為23情況下的重構(gòu)結(jié)果平均偏差

表2 測試驗證一中測量點數(shù)量為30情況下的重構(gòu)結(jié)果平均偏差

表3 測試驗證一中測量點數(shù)量為38情況下的重構(gòu)結(jié)果平均偏差

由表1、表2、表3可見,若以5 000次隨機選點測試的重構(gòu)結(jié)果平均偏差的平均值作為評判標準,普通克里金方法在結(jié)合對數(shù)變換后,在除Gaussian之外的模型上,重構(gòu)結(jié)果的準確度均有較大提升。其中,Exponential模型和Spherical模型的準確度要相對更優(yōu),在測量點數(shù)量為23、30、38時,分別為41.23%、31.89%、25.74%和34.98%、26.70%、21.51%。

2.2 測試驗證二

測試驗證二的場景相對更為復(fù)雜,場景中存在某大型設(shè)施作為屏蔽物,劑量率已知點散布在該設(shè)施周圍,總數(shù)為46個,其三維建模如圖2所示。

圖2 測試驗證二的三維建模圖

同樣對所有的已知點進行隨機選取,選取數(shù)量分別為總數(shù)的15%(7個點)、20%(10個點)和25%(12個點),重復(fù)次數(shù)為5 000次。將選取到的點作為測量點,剩余的點作為輻射場重構(gòu)的驗證點,在不同半方差擬合模型下重構(gòu)結(jié)果的平均偏差分別如表4、表5、表6所示。

表4 測試驗證二中測量點數(shù)量為7情況下的重構(gòu)結(jié)果平均偏差

表5 測試驗證二中測量點數(shù)量為10情況下的重構(gòu)結(jié)果平均偏差

表6 測試驗證二中測量點數(shù)量為12情況下的重構(gòu)結(jié)果平均偏差

由表4、表5、表6可見,在測試驗證二中,普通克里金方法在結(jié)合對數(shù)變換后,在除Gaussian外的模型上,重構(gòu)結(jié)果的準確度同樣均有提升。但與測試驗證一不同的是,在測試驗證二中除Gaussian外的4種模型的重構(gòu)準確度都較為接近,在測量點數(shù)量為7時,平均偏差的平均值為70.07%～74.01%;在測點數(shù)量為10時,平均偏差的平均值為63.57%～68.70%;在測點數(shù)量為12時,平均偏差的平均值為61.58%～66.01%。

2.3 綜合結(jié)果

綜合測試驗證一和測試驗證二的結(jié)果,對于結(jié)合對數(shù)變換的普通克里金方法,半方差擬合模型可選用在測試驗證一和測試驗證二中都表現(xiàn)較好的Exponential模型或Spherical模型,在隨機選點的情況下,其重構(gòu)結(jié)果平均偏差的平均值范圍為21.51%～74.01%。

此外,對于同樣數(shù)量的測量點和同樣的半方差擬合模型,重復(fù)選點測試中的最大平均偏差和最小平均偏差存在數(shù)量級的差異,說明測量點選取情況對重構(gòu)結(jié)果準確度的影響較大。因此,為了給實際工程應(yīng)用中的測量點選取工作提供指導,確保結(jié)合對數(shù)變換的普通克里金方法在輻射場重構(gòu)上的準確性,還必須對測量點的較優(yōu)選取方案進行研究。

3 測量點選取方案

考慮到克里金方法本身的空間插值特性,本文從測量點分布的均勻性出發(fā),分析其和輻射場重構(gòu)準確性之間的關(guān)系,從而確定出一種合適的測量點選取方案。為對測量點分布的均勻性進行量化,本文采用文獻[7]中的公式對其進行計算:

(14)

測量點數(shù)目越多,其均勻性的量化就越準確。本文選取測試一中測量點數(shù)目為38,且半方差擬合模型為Exponential和Spherical的實驗數(shù)據(jù)為研究對象,按照式(14)進行計算,繪制的散點圖如圖3所示。

(a) Exponential模型

由圖3可見,當采用結(jié)合對數(shù)變換的普通克里金方法進行三維輻射場重構(gòu)時,無論是采用Exponential模型還是Spherical模型,測量點分布的均勻性都和輻射場重構(gòu)的準確性存在一定的正相關(guān)性。因此為了確保輻射場重構(gòu)的準確性,應(yīng)用中測量點分布應(yīng)盡可能的均勻。

4 結(jié)語

本文對地統(tǒng)計空間插值中結(jié)合對數(shù)變換的普通克里金方法進行了介紹,并通過2個實測的輻射場數(shù)據(jù)對該方法進行了測試驗證。測試驗證結(jié)果表明,該方法在選取Exponential模型和Spherical模型為半方差擬合模型,且測量點為隨機選取的情況下,其重構(gòu)結(jié)果平均偏差的平均值處于21.51%～74.01%的范圍,顯示出該方法在實際應(yīng)用中具備一定的可行性。同時對測量點選取與重構(gòu)準確性之間的關(guān)系進行了初步分析研究,確定了測量點分布應(yīng)盡可能均勻的選取方案。