李 奇,樊元慶,王世聯(lián),趙允剛,常印忠,賈懷茂,劉蜀疆,張新軍
(禁核試北京國家數(shù)據(jù)中心和北京放射性核素實(shí)驗(yàn)室,北京 100085)
地下核試驗(yàn)產(chǎn)生的中子與圍巖中的鈣發(fā)生40Ca(n,α)37Ar反應(yīng)產(chǎn)生大量37Ar,根據(jù)周圍介質(zhì)中鈣含量的不同,37Ar產(chǎn)生量約為4×1011~5×1012Bq/kt[1]。氬常溫下為氣態(tài),空氣中含量約為0.93%,化學(xué)性質(zhì)穩(wěn)定,較容易通過巖石原有縫隙和核爆炸過程中產(chǎn)生的裂隙遷移到地表[2],因此,37Ar是全面禁止核試驗(yàn)條約(CTBT)現(xiàn)場(chǎng)視察(OSI)重要監(jiān)測(cè)核素之一。在OSI過程中,通過采集地表土壤中的空氣,分離濃縮得到氬樣品,在現(xiàn)場(chǎng)進(jìn)行測(cè)量或送至指定實(shí)驗(yàn)室進(jìn)行測(cè)量,確定樣品中的37Ar含量。由于37Ar遷移至地表過程復(fù)雜,且隨著時(shí)間推移不斷衰變,對(duì)于封閉較好的地質(zhì)結(jié)構(gòu),樣品中37Ar活度較低,因此,高靈敏地探測(cè)37Ar活度是地下核試驗(yàn)監(jiān)測(cè)的關(guān)鍵技術(shù)之一。
37Ar半衰期為35.04 d,軌道電子俘獲(EC)衰變至37Cl基態(tài),不發(fā)射β或γ射線,僅發(fā)射低能X射線和俄歇電子。其中,K殼層EC發(fā)射的X射線或俄歇電子能量為2.82 keV,分支比為0.902[3],L殼層EC發(fā)射的X射線或俄歇電子能量為0.27 keV,分支比為0.089,其他殼層EC發(fā)射的X射線或俄歇電子能量更低。一般探測(cè)設(shè)備很難探測(cè)這些低能射線,通常采用內(nèi)充氣正比計(jì)數(shù)管系統(tǒng)探測(cè)K殼層EC發(fā)射的X射線和俄歇電子測(cè)量37Ar活度,為提高探測(cè)靈敏度,通常在鉛屏蔽室外裝配塑料閃爍體探測(cè)器,通過反符合降低宇宙射線本底[4-6]。美國西北太平洋國家實(shí)驗(yàn)室在地下實(shí)驗(yàn)室中建立了37Ar活度測(cè)量裝置[7],進(jìn)一步提高了探測(cè)靈敏度。本文將介紹作為CTBT 16個(gè)國際放射性核素實(shí)驗(yàn)室和OSI指定場(chǎng)外分析實(shí)驗(yàn)室之一的北京放射性核素實(shí)驗(yàn)室建立的低本底內(nèi)充氣正比計(jì)數(shù)管(LBPC)37Ar活度測(cè)量裝置和方法。
建立的LBPC系統(tǒng)由內(nèi)充氣正比計(jì)數(shù)管、屏蔽體、塑料閃爍體探測(cè)器及相應(yīng)電子學(xué)系統(tǒng)組成。計(jì)數(shù)管置于屏蔽體內(nèi)部,長(zhǎng)度為590.3 mm,直徑為51 mm,內(nèi)部體積為1 003.8 mL,陰極材料(管壁)為高導(dǎo)電無氧銅;屏蔽體采用分層結(jié)構(gòu),外層為5 cm厚的普通鉛,內(nèi)層為2.5 cm厚的低本底鉛(210Pb比活度<5 Bq/kg),內(nèi)襯2 mm厚的無氧銅和5 mm厚的有機(jī)玻璃;屏蔽體上、下各放置一100 cm×60 cm×5 cm的塑料閃爍體,每個(gè)塑料閃爍體各裝配4支光電倍增管。塑料閃爍體探測(cè)器的8支光電倍增管由8路高壓提供偏壓,輸出信號(hào)經(jīng)各自的下閾甄別后“或”邏輯相加為1個(gè)信號(hào),經(jīng)延遲作為“反符合”信號(hào)屏蔽宇宙射線在計(jì)數(shù)管中產(chǎn)生的本底,內(nèi)充氣正比計(jì)數(shù)管由高壓電源提供偏壓,信號(hào)經(jīng)前置放大器、主放大器分為兩路,一路獲取反宇宙射線能譜,另一路獲取原始能譜。
將氧化鈣粉末用鋁箔包裹置入鋁罐中制成輻照靶,放入反應(yīng)堆堆芯輻照。輻照后的氧化鈣粉末裝入離心管中,離心管頂部用橡膠塞塞住并用蓋子蓋緊。用注射器向離心管中緩慢注入稀硝酸,待氧化鈣溶解后用注射器抽取上層氣體即得到37Ar樣品。
內(nèi)充氣正比計(jì)數(shù)管中間區(qū)域電場(chǎng)是均勻的,37Ar在該靈敏體積內(nèi)的探測(cè)效率可達(dá)100%,而其兩端區(qū)域由于電場(chǎng)畸變存在端效應(yīng),使得探測(cè)效率下降,因此需準(zhǔn)確刻度37Ar的探測(cè)效率。37Ar探測(cè)效率通常采用3H效率過渡[8]或通過長(zhǎng)度補(bǔ)償法得到37Ar活度標(biāo)準(zhǔn)實(shí)驗(yàn)刻度得到[9-10],本文采用長(zhǎng)度補(bǔ)償法準(zhǔn)確測(cè)量37Ar活度刻度得到了37Ar探測(cè)效率,效率計(jì)算公式為:
(1)
式中:ε為L(zhǎng)BPC對(duì)37Ar的探測(cè)效率;n為L(zhǎng)BPC測(cè)量的37Ar的K殼層俘獲俄歇電子和X射線的凈峰計(jì)數(shù)率,s-1;C為37Ar的活度濃度,Bq/mL;V為L(zhǎng)BPC的體積,mL;PK為K殼層EC分支比。
在制樣裝置(圖1)上,用注射器將37Ar氣體樣品注入至制樣裝置的進(jìn)樣口,37Ar和工作氣體(氬甲烷)在混氣室中充分混合后,同時(shí)充入長(zhǎng)度不同而其他結(jié)構(gòu)完全相同的長(zhǎng)、中、短3支內(nèi)充氣正比計(jì)數(shù)管以及本文建立的LBPC中,充氣壓力為200 kPa,此時(shí),充入這些計(jì)數(shù)管中的37Ar活度濃度相同。1次測(cè)量的37Ar能譜示于圖2,高斯擬合得到各計(jì)數(shù)管的37Ar凈峰面積。長(zhǎng)、中、短3支計(jì)數(shù)管通過長(zhǎng)度補(bǔ)償法準(zhǔn)確測(cè)量37Ar活度濃度。
圖1 制樣裝置示意圖Fig.1 Schematic diagram of sample preparation device
圖2 長(zhǎng)、中、短內(nèi)充氣正比計(jì)數(shù)管及LBPC的37Ar樣品譜Fig.2 37Ar spectra of long, medium, and short internal gas proportional counters and LBPC
長(zhǎng)、中、短3支計(jì)數(shù)管直徑均為38 mm,端結(jié)構(gòu)相同,內(nèi)部體積分別為361.3、266.3、171.5 mL。任意兩支計(jì)數(shù)管的37Ar凈峰計(jì)數(shù)率差除以體積差(長(zhǎng)度補(bǔ)償),即可消除端效應(yīng)影響。由于計(jì)數(shù)管端結(jié)構(gòu)不可能完全一致,因此,采用長(zhǎng)度不同的3支計(jì)數(shù)管對(duì)其計(jì)數(shù)率和體積作線性擬合,擬合得到的斜率除以37Ar K殼層EC分支比即得到37Ar活度濃度,斜率的相對(duì)標(biāo)準(zhǔn)偏差為長(zhǎng)度補(bǔ)償引入的不確定度。內(nèi)充氣正比計(jì)數(shù)管測(cè)量氣體活度還存在壁效應(yīng)影響,壁效應(yīng)是充氣壓力的函數(shù),充氣壓力越大壁效應(yīng)越小,文獻(xiàn)[9]研究表明充氣壓力為200 kPa時(shí),壁效應(yīng)影響小于1%,可忽略。37Ar活度濃度測(cè)量結(jié)果列于表1(活度濃度參考時(shí)間為2017-03-08 12:00),共進(jìn)行了7組測(cè)量,測(cè)定的37Ar活度濃度為(2.21±0.02) Bq/mL(k=1)。活度濃度不確定度主要包括長(zhǎng)度補(bǔ)償引入的不確定度0.5%、計(jì)數(shù)管體積不確定度0.5%、統(tǒng)計(jì)不確定度0.5%、K殼層EC分支比的不確定度0.3%。
表1 37Ar活度濃度測(cè)量結(jié)果Table 1 Measurement result of 37Ar activity concentration
根據(jù)測(cè)定的37Ar活度濃度以及LBPC體積和測(cè)量的37Ar凈峰計(jì)數(shù)率(已校正至活度濃度參考時(shí)間),由式(1)計(jì)算得到LBPC對(duì)37Ar的探測(cè)效率為0.899±0.015(k=1),測(cè)量結(jié)果列于表2。探測(cè)效率不確定度主要包括37Ar活度濃度不確定度1%、計(jì)數(shù)管體積不確定度0.5%、統(tǒng)計(jì)不確定度0.5%、K殼層EC分支比的不確定度0.3%。
表2 37Ar探測(cè)效率測(cè)量結(jié)果Table 2 Measurement result of 37Ar detection efficiency
37Ar活度A和最小可探測(cè)活度(MDA)計(jì)算公式分別為:
(2)
(3)
實(shí)驗(yàn)測(cè)量了LBPC的37Ar反宇宙射線能譜和原始能譜(圖3),利用反宇宙射線能譜和原始能譜計(jì)算的37Ar活度列于表3。從能譜和活度測(cè)量結(jié)果看出,反宇宙射線能譜和原始能譜計(jì)算的37Ar活度相差僅3%,反宇宙射線能譜37Ar峰明顯突出,不確定度顯著降低。
圖3 37Ar反宇宙射線能譜和原始能譜Fig.3 Anti-cosmic spectrum and original spectrum of 37Ar
表3 反宇宙射線能譜和原始能譜計(jì)算的37Ar活度Table 3 Activity of 37Ar from anti-cosmic spectrum and original spectrum
實(shí)驗(yàn)獲取了LBPC本底的37Ar反宇宙射線能譜和原始能譜(圖4),測(cè)量時(shí)間為390 000 s。相對(duì)于原始本底能譜,反宇宙射線本底能譜中37Ar峰區(qū)本底計(jì)數(shù)率降至原來的7.3%,僅0.016 s-1。測(cè)量24 h,37Ar的MDA達(dá)2.5 mBq。
研究建立了LBPC37Ar活度測(cè)量裝置和方法。采用長(zhǎng)度補(bǔ)償法測(cè)量37Ar活度準(zhǔn)確刻度了37Ar的探測(cè)效率;采用塑料閃爍體探測(cè)器反符合屏蔽宇宙射線在計(jì)數(shù)管中產(chǎn)生的本底,37Ar峰區(qū)本底計(jì)數(shù)率降至原來的7.3%,測(cè)量24 h,37Ar的MDA達(dá)到了2.5 mBq。
圖4 37Ar反宇宙射線本底能譜和原始本底能譜Fig.4 Anti-cosmic background spectrum and original background spectrum of 37Ar