崔 堯,張向陽,何高魁
(中國原子能科學研究院,北京 102413)
核燃料組件無損檢測探測系統(tǒng)設計
崔 堯,張向陽,何高魁
(中國原子能科學研究院,北京 102413)
燃料組件是反應堆的核心部分,在高溫、高壓及強中子輻射場等復雜環(huán)境條件下,燃料棒中芯塊會出現(xiàn)腫脹、變形甚至包殼破裂,嚴重威脅反應堆的安全運行。為了更好地了解燃料組件在反應堆內(nèi)的變化,研究高燃耗的燃料組件中燃料棒的中心空洞形成和燃料棒的變形情況,高能X射線無損檢測是一種有效的技術手段。由于輻照后核燃料組件自身具有強放射性,探測系統(tǒng)設計中必須考慮減弱燃料組件自身輻射對探測采集的影響,因此組件探測系統(tǒng)中探測器陣列及準直器的優(yōu)化設計十分必要。經(jīng)過建模及相關模擬計算,得到了探測器單元最佳尺寸,優(yōu)化了后準直器的結構設計,為提高燃料組件無損檢測系統(tǒng)重建圖像的質(zhì)量提供幫助。
高能X射線無損檢測;烏酸鎘晶體;能量沉積率;串擾率;后準直器
燃料組件是反應堆的核心部分,燃料組件的安全性對于整個反應堆的安全運行有著重要的意義。在高溫、高壓及強中子輻射場等復雜環(huán)境條件下,燃料棒中芯塊會出現(xiàn)腫脹、變形甚至包殼破裂,嚴重威脅反應堆的安全運行。為了研究核反應堆燃料組件的安全性和完整性,了解其熱力學性能和機械性能非常重要。通過傳統(tǒng)的檢測方法,如金相學,來檢測燃料棒中心空洞和裂紋,需要花費大量的時間,且需要破壞燃料組件自身的結構,無法檢測燃料棒的形變。X射線無損檢測技術可以簡易且快速的檢測燃料組件內(nèi)部的結構,并且費用低,不產(chǎn)生廢物[3]。
探測系統(tǒng)是整個燃料組件無損檢測系統(tǒng)中的核心部分,其性能對于重建圖像的質(zhì)量有著重要的影響。為了提高高能X射線在探測器中的能量沉積率,一般選用鎢酸鎘(CdWO4)晶體。鎢酸鎘晶體具有密度大(ρ=7.9 g/cm3)、原子序數(shù)高(Z=64.2)、熒光轉換率(1.2×104~1.5×104光子/MeV)高等優(yōu)點,且其熒光波長(470/540 nm)與一般的光電二極管敏感波長相匹配,在高能X射線無損檢測中得到了廣泛的應用[5]。
針對核燃料組件探測的特殊性,本文通過蒙特卡羅方法,研究了在9 MeV直線加速器下,不同尺寸的CdWO4晶體的能量沉積率和串擾率,得到了最佳的晶體尺寸;模擬了燃料組件自身強輻射對數(shù)據(jù)采集的影響,優(yōu)化后準直器的設計從而將燃料組件自身輻射的影響降至最低。經(jīng)過優(yōu)化設計的探測系統(tǒng)將為改善重建圖像的質(zhì)量奠定基礎。
1.1 蒙特卡羅方法
蒙特卡羅方法是以概率統(tǒng)計為理論基礎的一種數(shù)學計算方法。蒙卡方法可以真實地模擬實際粒子運動的物理過程,其仿真結果與實際結果比較吻合,可以得到滿意的結果。MCNP可以解決電子、光子以及中子的聯(lián)合運輸問題,光子和電子的能量范圍從1 keV到1 000 MeV[4],選用MCNP5版本,抽樣次數(shù)為一百萬次,準直器孔徑選用0.4 mm。
1.2 CdWO4晶體能量沉積率的估算
在以往的文獻中,采用蒙特卡羅方法研究探測器的能量沉積率,大多采用單能的X射線。對于9 MeV直線加速器,其產(chǎn)生的X射線為連續(xù)譜,為了更好地接近實際,采用麥克斯韋譜模擬加速器所產(chǎn)生的連續(xù)譜(如圖1),并在此基礎上估算晶體的能量沉積率。所采用的模型如圖2,CdWO4晶體高為5 mm,X射線沿x軸入射。
圖1 模擬加速器所產(chǎn)生的連續(xù)譜Fig.1 Continuous spectrum generated by accelerator
圖2 模擬模型Fig.2 Simulation model
改變CdWO4晶體截面寬度和晶體長度進行模擬,可以得到其能量沉積率與截面寬度和晶體長度之間的關系,如圖3所示。從圖3中可以看出,隨著晶體截面寬度和晶體長度的變大,能量沉積率隨之增大。從圖中可以得出,在晶體尺寸為2 mm×5 mm×30 mm時,能量沉積率為38%。
圖3 能量沉積率與晶體尺寸的關系Fig.3 The relationship between energy deposition rate and the size of scintillation
當進一步增加晶體的橫截面積時發(fā)現(xiàn),CdWO4晶體截面尺寸對能量沉積率的影響逐步減小。因此,進一步模擬當長度一定時(30 mm),晶體能量沉積率與晶體截面寬度之間的關系,如圖5所示。
圖4 增加晶體尺寸Fig.4 Increasing the size of the scintillation
圖5 能量沉積率與截面寬度的關系Fig. 5 The relationship between energy deposition rate and the width of cross section
當截面寬度較小時,隨著截面寬度的增長,能量沉積率快速的增長,當截面寬度大于3 mm時,增長變緩。截面寬度從1 mm增加到2 mm時,能量沉積率增長了19.1%;截面寬度從3 mm增加倒4 mm時,能量沉積率僅增長了2.6%。
因此,綜合考慮能量沉積率以及系統(tǒng)空間分辨率,最佳晶體尺寸選擇為3 mm×5 mm×30 mm。
1.3 串擾率
通過輻照后核燃料組件的探測環(huán)境十分復雜,各種背景干擾都將對探測成像產(chǎn)生嚴重的影響,其中探測單元之間的串擾也是不容忽視的問題。如圖6所示,經(jīng)過準直器入射到CdWO4晶體內(nèi)的X射線,其中的一部分散射光子、光電子以及X射線等會被相鄰的探測器采集到,形成串擾信號,對重建圖像的質(zhì)量產(chǎn)生很大影響。
改變相鄰兩塊CdWO4晶體之間的隔離層的厚度,模擬串擾率與隔離層厚度之間的關系,隔離層采用鉛作為屏蔽材料。通過蒙特卡羅模擬,選用寬度為0.4 mm的準直器,隨著隔離層厚度的改變,串擾率隨之改變,結果示于圖7。
圖6 串擾模擬Fig.6 Crosstalk simulation
圖7 串擾率與隔離層厚度的關系Fig.7 The relationship between crosstalk and the thickness of isolation layer
由圖7可以看出,在沒有隔離層的情況下,串擾率為11.5%,當隔離層為0.8 mm時,串擾率下降為1.6%。
圖8 不同隔離材料的串擾率Fig.8 The crosstalk ratio of different materials
改變隔離層的材料,用Wu代替Pb進行模擬計算,計算結果如圖8所示。在隔離層厚度處于0.2 mm到0.4 mm之間時,由于Wu的密度大于鉛,可以更好地隔離X射線,減小串擾率。當隔離層厚度大于0.6 mm時,兩者的隔離作用趨同,0.6 mm時僅相差0.41%。
因此,隔離層最佳厚度選為0.6 mm,隔離材料選用鉛。
由于輻照后的燃料組件具有很強的放射性,乏燃料經(jīng)冷卻后,大部分放射性源自銫-137和鍶-90。鍶-90發(fā)生β衰變,其半衰期為28.1 a,銫-137的半衰期為30.17 a。銫-137衰變放出的能量為662 keV的光子是輻照后燃料組件強輻射的主要來源。燃料組件每秒約產(chǎn)生1013個光子,其輻射出的光子進入探測器,形成背景輻射,對信號采集和圖像重建產(chǎn)生嚴重的影響。采用的模型中,燃料組件活性區(qū)長度為1.5 m,其輻射出的一部分光子可以不經(jīng)過后準直器而直接進入探測器,從而產(chǎn)生很大的噪聲信號。因此,需要對后準直器進行優(yōu)化設計,減弱燃料組件自身輻射帶來的影響。
通常情況下,后準直器,前縫板縫寬5 mm,中縫板垂直縫寬為0.4 mm,若采用此種設計,經(jīng)模擬,每秒鐘燃料組件輻射出的光子進入探測器的數(shù)量約為3×107,若不加以優(yōu)化設計,將對圖像重建產(chǎn)生嚴重影響。
為了降低燃料組件輻射的影響,加速器工作在脈沖狀態(tài),探測器與加速器同步工作。當加速器脈沖頻率為250 Hz,脈沖持續(xù)時間為2.5 μs時,中心軸線上距離靶1 m處X射線劑量率為3 000 cGy/min,加速器距離探測器的1.4 m,晶體尺寸為4 mm×7 mm×30 mm。經(jīng)計算,單個脈沖時間內(nèi),單個探測單元接受到的由加速器產(chǎn)生的X射線光子數(shù)為2.7×106。由于燃料組件自身的強輻射,單個探測單元接受到的來自燃料組件的光子數(shù)約為75。
優(yōu)化后準直器的設計模型示于圖9。在中縫板的后方,增加后縫板,即在探測模塊上下兩側增加屏蔽板,以屏蔽燃料組件輻射出的直接被探測器接收的光子。屏蔽板采用鎢鐵鎳合金,厚度為15 mm。在此基礎上進行模擬,一個脈沖時間內(nèi),進入探測器的光子數(shù)量減小為7.6。減小了背景輻射在探測器中的貢獻,有利于圖像重建質(zhì)量改善。
圖9 優(yōu)化后的模型Fig.9 The optimal model
通過對燃料組件無損檢測探測系統(tǒng)的模擬計算,實現(xiàn)了對探測系統(tǒng)的優(yōu)化設計。通過模擬計算,綜合考慮CdWO4晶體截面寬度、晶體厚度、晶體之間隔離材料、隔離材料厚度等因素,可以得出,在晶體截面寬度為3 mm,晶體厚度為30 mm,隔離材料選用鉛,且隔離層厚度為0.6 mm時,由CdWO4晶體構成的線性陣列探測器的性能達到最佳狀態(tài)。通過模擬,對后準直器進行優(yōu)化,增加了后縫板,減小了背景輻射在探測器中的貢獻。
[1] 樓琪. 用于高能X射線成像的CdWO4閃爍體探測器靈敏度的研究[J]. 核電子學與探測器,2002,22(4):331-336.
[2] 陳志強. 高能射線工業(yè)CT最新進展[J]. CT理論與應用研究. 2005. 14(4):1-4.
[3] Kozo Katsuyama. High energy X-ray CT study on the central void formations and the fuel pin deformations of FBR fuel assemblies[J]. Nuclear Instruments and Methods Physics Research B, 2007, 255: 365-372.
[4] 田麗. 工業(yè)CT圖像的探測器串擾偽影研究[D]. 重慶:重慶大學,2001,4:1-71.
[5] 張朝宗. 工業(yè)CT技術和原理[M]. 北京:科學出版社,2009:110.
Detection System for the Fuel Assembly Nondestructive Testing
CUI Yao, ZHANG Xiang-yang, HE Gao-kui
(ChinaInstituteofAtomicEnergy,Beijing102413,China)
Fuel assemblies are the central components of a reactor. The core fuel pellets in the fuel pins will swell and deform and the fuel cladding may even break under the complex environment of high temperature, high pressure and intense neutron radiation field, which threats the safety of the reactor. To better understand the changes in the behavior of the fuel assembly in the reactor and study the central void formations and deformations of fuel pins in fuel assemblies to high burn-up, high-energy X-ray non-destructive testing is an effective technical means. Irradiated nuclear fuel assembly has a strong radioactivity, it is necessary to optimize the design of the detector system and the collimator to reduce the effect from gamma rays emitted from the irradiated fuel assembly during detection system designing phase. Through modeling, estimating and optimization, the optimal size of the detector unit is obtained and the collimator design is optimized which can lay the foundation to improve the quality of the reconstructed images of the fuel assembly nondestructive system.
high energy X-ray nondestructive testing; CdWO4scintillation; energy deposition rate; ratio of crosstalk; collimator
10.7538/tws.2015.28.03.0167
2015-04-14;
2015-05-17
崔 堯(1990—),男,山西長治人,碩士研究生,核技術及應用專業(yè)
TL364+2
A
1000-7512(2015)03-0167-04