柏 磊,王仲奇,邵婕文,何麗霞
(中國原子能科學研究院 放射化學研究所,北京 102413)
核燃料生產過程中,在化工轉化和產品加工時會產生大量含鈾殘余物料,目前采用回收利用后再裝桶暫存的處理方式,這些殘余物料經多次循環(huán)處理之后,累積形成一定量的中低密度非均勻核廢料,根據我國核材料管制及放射性廢物管理相關規(guī)定,其中的核材料量應準確測定。這類核廢料介質裝填密度較小,介質和核材料分布不均勻,對其中核材料進行定量分析時,很難取得代表性的樣品,部分樣品甚至無法采用化學分析方法精確測定其中的核材料量,只能采用非破壞性測量技術進行分析。通過采集分析放射性物質發(fā)射的中子、γ信號或熱功率數據反推其中的核材料量,通常需要采用相近的工作標準樣品完成系統(tǒng)探測效率刻度而進行測量數據溯源。對于中低密度非均勻核廢料,源于其介質和放射性分布不均勻的特征,很難制備具有代表性的標準樣品,結合實際情況,宜采用無源探測效率刻度技術進行儀器設備的探測效率刻度。
本研究擬建立非均勻含鈾物料測量系統(tǒng)的無源效率刻度方法,在實驗室條件下進行標準點源和硝酸鈾酰溶液工作樣品的實驗驗證,并針對實際的“酸不溶渣”樣品,通過增量法原理進行測量結果評價的嘗試性研究,對于采用該方法進行非均勻物料測量的結果評價具有一定的指導意義。
圖1 非均勻含鈾物料測量裝置Fig.1 Heterogeneous uranium-bearing material measurement equipment
非均勻含鈾物料測量裝置示于圖1。該專用系統(tǒng)采用分段γ射線掃描(segmented gamma scanning, SGS)測量技術[1-3]原理和設計,根據核燃料生產設施核廢料的幾何特征,進行技術集成后達到定制開發(fā)目標。系統(tǒng)由透射源、樣品旋轉平臺、探測器、自動控制和數據分析軟件等部分構成,透射源用于獲取基體的密度分布;樣品旋轉平臺用于樣品放置和分層級均勻化處理;探測器為高分辨率γ射線探測器,自帶準直器和屏蔽體,用于采集透射源和自發(fā)射γ射線能譜。
采用SGS測量技術進行含鈾物料分析時,通過樣品徑向旋轉、軸向分段的方式實現樣品分層和層內介質及放射性近似均勻化處理。若采用標準樣品進行效率刻度和數據校正,即為過質量刻度的方式,則核廢料與標準樣品的測量數據之間滿足如下關系:
M=TCR/(K0×K)
(1)
式中,K0為單位質量核廢料發(fā)射的某特定能量γ射線的發(fā)射率;K為刻度因子,通常由相關標準樣品刻度測量得到;TCR為經過校正后樣品的特定能量γ射線的總計數率。
在得到某特定能量γ射線的發(fā)射率之后,利用系統(tǒng)的探測效率進行運算,可推導出其中的放射性物質或者核材料的質量,計算關系如式(2):
(2)
其中,ε為探測器的探測效率,C(E)為每層樣品中特征γ射線的實際發(fā)射率,C為探測器測量得到的每層物料特定射線的計數率。
式(2)中探測效率ε為關鍵數據,可采用有源或者無源方式獲得,有源即采用實物標準樣品,無源即采用模擬計算的方式。無源探測效率刻度過程中,利用蒙特卡羅粒子輸運工具模擬特征γ射線在樣品、探測器、準直器及屏蔽體中的輸運過程,記錄其在探測器靈敏體積中能量沉積譜,獲取特定幾何條件下特征γ射線的探測效率[4],本研究過程中,通過探測器表征、點源驗證、條件驗證等步驟實現。
探測器表征,指基于探測器出廠參數,采用表征實驗與模擬計算結合確定探測器的結構參數。表征實驗選用標準點放射源,精確定位后測定不同源探距條件下對應的特征γ射線的絕對探測效率,模擬計算則基于探測器幾何模型。通常,準確的探測器外部尺寸獲取較為容易,工藝相關參數難于得到,兩種參數同時影響測量和計算效率[5-6]。探測器表征過程中,根據其外部尺寸和工藝參數典型值進行初始化,建立含外殼、支架、晶體、冷指、冷孔多參數的探測器幾何模型,如圖2所示。
圖2 蒙特卡羅模擬計算的HPGe探測器幾何模型Fig.2 HPGe detector geometry figure used for Monte Carlo simulation
探測器表征過程中,逐步調整蒙特卡羅計算所用探測器模型的幾何參數,使絕對效率計算值逐步逼近最終與實驗值保持一致。非均勻含鈾物料測量系統(tǒng)所用探測器的表征結果列于表1,探測器的絕對探測效率計算值與實驗值之間偏差小于2%。
探測器表征之后,在非均勻含鈾物料測量系統(tǒng)平臺上進行無源探測效率刻度模型驗證。探測器外加鉛屏蔽體和準直器,放射源采用支架定位放置在樣品旋轉平臺中心、探測器晶體軸的延長線上;體放射源采用SGS測量工作模式分別完成透射和自發(fā)射γ射線能譜測量分析。
點源驗證采用137Cs和60Co標準點放射源,實驗現場如圖3所示,樣品旋轉平臺固定,能譜測量時間根據特征γ射線的全能峰面積而定,保證統(tǒng)計漲落小于0.5%。
表1 探測器絕對探測效率實驗測量與模擬計算結果比對Table 1 Comparison of absolute detection efficiency between measurement and simulation
圖3 探測器表征驗證實驗現場Fig.3 The confirmatory experiment of detector characterisation
經測量與計算,探測器的絕對探測效率計算值與實驗值數據見表2,二者之間偏差小于2%。
采用SGS工作模式測量經過質譜分析方法定值的硝酸鈾酰溶液樣品,進行無源探測效率刻度模型驗證。實驗現場如圖4所示,硝酸鈾酰溶液樣品放置在旋轉平臺中心位置,每個樣品分為4層進行連續(xù)掃描透射測量結合固定自發(fā)射方式測量,透射和自發(fā)射測量時間分別設定為100 s和2 000 s。
表2 裝置絕對探測效率的蒙特卡羅模擬計算與實驗數據比對Table 2 Comparison of absolute detection efficiency between measurement and simulation
圖4 硝酸鈾酰溶液樣品驗證實驗現場Fig.4 The confirmatory experiment with uranium nitrate
透射源選用152Eu,標稱活度10 mCi,選取其能量為122、244、344、779、964、1 408 keV的特征γ射線作為透射效率擬合射線。硝酸鈾酰溶液盛裝在直徑為9.6 cm、壁厚0.3 cm、高度15 cm聚乙烯桶中,分別編號為1、2、3號,其中含235U分別為8.35、4.39、2.16 g。
首先,打開透射源獲取有樣品和無樣品條件下透射源的特征γ射線能譜,通過解譜計算得到各特征γ射線能峰的計數率,兩者相比即可獲得各特征γ射線的透射率T,透射率與能量之間經擬合后再插值計算介質密度,透射率T與特征γ射線能量之間的關系見下式[7]:
(4)
式中,T為測量對象對能量為E的特征峰的透射率, a、b、c、d為擬合參數。
結合透射率及探測器表征模型,建立非均勻含鈾物料體源測量蒙特卡羅計算模型[8],計算得出體源驗證系統(tǒng)絕對探測效率矩陣,數據見表3,由表3數據可以看出,三個樣品的絕對探測效率較為一致。
表3 絕對探測效率矩陣Table 3 The matrix of absolute detection efficiency
自發(fā)射測量獲取硝酸鈾酰的γ射線能譜,經解譜計算其中235U發(fā)射的能量為186 keV的特征γ射線計數率,利用表3所列探測效率,結合公式(1)和式(5)計算每層樣品中所含235U的質量(g)和樣品所含總235U的質量:
(5)
表4 硝酸鈾酰樣品235U質量分析結果Table 4 Resluts of 235U mass in uranium nitrate sample
利用經過標準點源和具有化學方法定值的均勻硝酸鈾酰溶液樣品驗證的SGS裝置,對需要進行測量分析的實際樣品進行應用嘗試。
根據核材料衡算需求,選用“酸不溶渣”(圖5)作為測量分析應用對象。所謂酸不溶渣,是鈾生產工藝線殘留物質,由大小不一的結塊組成。這種物質不僅在介質密度和核材料分布方面具有明顯非均勻特性,同時還難以完全溶解于酸溶液。這決定了酸不溶渣無法采用化學方法對其中的核材料(235U)含量進行測定,只能采用非破壞性分析方法進行測量分析。
采用SGS測量分析工作模式,可以給出酸不溶渣中核材料含量的估計值。考慮酸不溶渣的非均勻分布和難溶于酸的兩個主要特性,本文采用增量法對酸不溶渣的SGS測量分析結果進行評價。
圖5 “酸不溶渣”樣品形態(tài)及包裝Fig.5 The form and package of sample
增量法源于化學分析中被廣泛用于檢驗儀器準確度的標準加入法,通常被用在不具備刻度標準曲線的場景,特點就是將樣品和標準品混在一起同時測定。在非破壞性測量分析中(non-destruction analysis, NDA),特別是在SGS測量分析工作模式中,尚未見到采用增量法進行分析結果準確度評價的報道。
非破壞性測量分析中標準樣品的加入方法與化學分析中均勻式加入方式不同,因此如何將非均勻性加入等價于均勻性加入,是采用增量法評價測量分析結果的關鍵。
基于非對等性刻度方法研究[10],選取合適的徑向位置r(r/r0=0.75,其中r為點源距離樣品軸線的距離,r0為樣品的半徑),將工作樣品放入該處,可以近似等價于加入均勻體樣品。
基于4個鈾工作標準樣品,參照標準加入法,將未加入工作標準樣品的酸不溶渣樣品視為基礎樣品,利用經過標準點源和具有化學方法定制的均勻硝酸鈾酰溶液工作標準樣品驗證的無源效率刻度SGS裝置,測量其中所含的235U質量,采用增量法評價其測量分析結果的準確性。
在經過測量分析的基礎樣品中,依次將4個點源逐個埋入酸不溶渣的確定位置,并依次測量分析總235U質量,根據測量數據進行增量分析并判斷其準確性,進而評估基礎樣品的測量分析結果準確性。
基于SGS方法中的各層均勻性的假設,結合非對等刻度方法,根據樣品各層的透射率數據,開展模擬計算,計算探測器對于該層的探測效率與上述點源的探測效率,并給出兩者之間的相對偏差,見表5。相對偏差在2%以內,與上文表征偏差數據一致。
表5 體源與點源探測效率偏差Table 5 The relative deviation of point source and body source
按照上述流程,酸不溶渣測量數據見表6。樣品的多次測量相對標準偏差在1%左右,其主要來源與能譜測量的統(tǒng)計性偏差(該裝置采用高精度伺服電機控制探測器上下移動和樣品旋轉,因電機走位帶來的重復性誤差影響非常小),測量得到的增量235U與工作樣品標稱值之間的偏差在8.5%以內。相對于均勻樣品的測量偏差,對于采用增量法的“酸不溶渣”測量偏差相對較大,其主要原因包括SGS方法中采用每層均勻的假設與實際情況不符合帶來的誤差,上述非均勻刻度方法引入的模擬誤差,非理想點源引入的體積誤差,以及點源放置位置(r/r0=0.75)帶來的方法誤差。此位置主要由不同的射線能量、不同的介質確定,工業(yè)應用上選擇固定的0.75必然與真實的結果存在誤差,依據參考文獻[10]和[11],可通過實驗結果(介質透射率)和模擬計算給出此時的位置假設誤差為1.45%。
表6 235U質量分析結果Table 6 Results of 235U mass
根據測量分析結果,將235U的測量值與標準樣品數量(測量基礎樣品時,數量按照0代入)作線性擬合處理,效果見圖6,兩者之間的相關系數R2=0.995 4,即標準樣品中235U的質量分析結果與標稱值間較為吻合。
圖6 測量235U質量與顆粒個數曲線Fig.6 The curve between the mass of 235U and number of samples
由表5和圖6數據可得出,利用經過標準點源和具有化學定值的均勻硝酸鈾酰溶液樣品驗證的SGS裝置,對非均勻含鈾物料酸不溶渣的測量分析結果為0.75 g。采用增量法擬合外推方式,推算酸不溶渣中235U的質量為0.74 g,與實測值偏差為-1.87%。
本研究建立了非均勻含鈾物料測量系統(tǒng)的無源效率刻度技術,并在實驗室條件下進行了標準點源和具有化學方法定值的均勻硝酸鈾酰溶液工作樣品實驗測量分析驗證。驗證測量分析實驗數據的相對偏差均小于2%。
針對非均勻、難溶于酸溶液的“酸不溶渣” 典型實際樣品,測量分析了其中核材料的含量,并嘗試利用增量法評價測量分析結果的準確性,針對增加的工作樣品測量相對偏差在8.5%以內。原則上講,對于SGS這類一維體素劃分方法測量分析結果的評價需要在更高的維度上進行[8,12]??紤]到實際情況,直接沿用化學分析增量法進行實際樣品測量分析結果評價嘗試,得到的評估結果可能與實際情況存在一定的偏差,針對不均勻程度更大的測量對象,SGS方法中層內均勻的假設帶來的誤差影響必然更大。
計劃利用TGS[12]和SGS、SRGS[13]測量分析方法,結合非對等刻度方法[10],從樣品各層透射測量的分析結果出發(fā),通過模擬計算的方式,得到探測器對于加入樣品與基礎樣品的探測效率的差異,針對無源刻度SGS裝置建立更為有效可靠的測量分析結果準確性評價方法。同時,今后將陸續(xù)開展不同含量的放射性物料和工作點源組合測量研究,針對測量方法與誤差關系開展下一步的應用研究工作。