氡釷射氣連續測量裝置的研制及其現場驗證

劉嘯晨,楊明理,王 攀,初旭陽

(核工業北京化工冶金研究院,北京 101149)

連續監測氡與釷射氣在輻射防護領域具有重要意義。氡與釷射氣及其子體是人體所受天然放射源內照射劑量的主要貢獻者,在高釷本底地區或稀土廠房,釷射氣對內照射劑量的貢獻甚至高于氡[1-3]。因此,為了輻射防護與劑量監測工作的開展,需要對現場環境中的氡與釷射氣的放射性濃度進行連續監測。

氡與釷射氣的半衰期不同,在不同材料中的擴散系數以及發射出的α粒子能量也不同,在氡與釷射氣的混合環境中通常利用這三方面的不同進行甄別測量[4-8]?，F有氡釷射氣監測裝置,按照原理主要分為固體徑跡探測器、脈沖電離室探測器、半導體探測器、閃爍室探測器四類。固體徑跡探測器體積小,攜帶方便;但所得到的是一段時間內氡與釷射氣造成的累積效應,無法用于連續監測,無法得到一段時間內氡與釷射氣所致劑量的變化情況。半導體探測器與脈沖電離室探測器可以通過能量分辨識別氡與釷射氣并進行連續監測;但半導體探測器在高濕度環境中受到復合損失的影響較大,脈沖電離室探測器的最小探測限比半導體探測器更低,當環境中氡與釷射氣的濃度之比超出(1∶5)～(5∶1)范圍時,對其中相對較少的核素測量值準確性下降。在所有探測器中,閃爍室探測器具有計數穩定性好、受濕度影響較小,以及抗電磁干擾能力強等特點[9-12];但使用閃爍室探測器進行長期監測時,可能受到前序周期生成的附壁釷子體影響(氡子體半衰期較短,在長期監測中影響遠遠小于釷子體)。因此,排除濕度變化、電磁干擾以及附壁釷子體的影響,準確地對氡與釷射氣進行甄別測量,是氡釷射氣聯合監測的關鍵。

筆者使用閃爍室構造氡與釷射氣連續監測裝置并解決濕度變化影響,建立了氡與釷射氣分辨測量、連續監測方法,研制了氡釷射氣連續測量裝置原型機(簡稱為原型機),并對原型機進行了現場測試。

1 氡釷射氣連續監測原理和裝置研制

1.1 監測原理

利用氡與釷射氣半衰期的差異,在2個探測單元間增加1個延遲器,如果采樣氣流中的釷射氣在延遲單元中接近完全衰變,則可以認為前1個探測器所產生的計數是由氡與釷射氣及其子體衰變發出的α粒子所致,后1個探測器所產生的計數完全由氡及其子體衰變發出的α粒子導致,前后探測單元的計數相減則可近似得到由釷射氣及其子體產生的影響,即通過2個探測器產生的計數反推2種核素的放射性濃度。

1.2 裝置研制

原型機主要由延遲器、閃爍室探測器、光電倍增管、單片機、電源構成(圖1),閃爍室記錄放射性衰變產生的計數,并將結果依次上傳至信號采集單元、單片機和上位機。上位機接收裝置各元件上傳的信號,并通過單片機對高低壓電源、甄別閾值、測量周期、采樣流率等進行控制。

圖1 氡釷射氣連續監測原型機框圖Fig. 1 The diagram of radon-thoron continuous monitoring prototype machine

氡與釷射氣在閃爍室內產生的子體,一部分附在閃爍室壁上(附壁效應),另一部分會隨采樣氣體流出閃爍室,附壁比例主要受濕度與閃爍室結構的影響。目前已有附壁比例不受濕度影響的改進型ZnS閃爍室[13],為排除環境濕度對結果的干擾,本研究使用此改進型閃爍室來構建原型機。

采樣氣體先后流經過濾器、第1閃爍室、延遲器、第2閃爍室。氣體中原有的放射性子體被過濾器全部過濾,在第1閃爍室中氡、釷射氣及其生成的子體衰變產生的α粒子被計數(以下簡稱通道計數);采樣氣體流入延遲器經過較長時間(釷射氣半衰期的8倍以上)后再進入第2閃爍室,在其中僅有氡及其子體衰變產生的α粒子被計數。通過第2閃爍室的計數可以計算出采樣氣體中的氡濃度,而通過2個閃爍室在測量周期內通道計數之差則可計算出采樣氣體中的釷射氣濃度。

2 氡釷射氣連續監測裝置性能測試

2.1 遺留系數標定

氡子體中半衰期最長的RaB的半衰期僅有26.7 min,按6倍半衰期考慮,附壁氡子體對相隔3 h后的測量值的影響已經可以忽略。在長周期(測量周期24 h)取樣監測過程中,由于釷子體中存在半衰期較長的核素,因此前序采樣周期中附壁的釷子體依然會對后續周期的通道計數產生貢獻。將初始采樣周期在一定周期后遺留的附壁釷子體所產生的通道計數與初始采樣周期內由釷射氣及其子體產生的通道計數之比定義為遺留系數。為了排除在長期監測中附壁釷子體的影響,建立了數學模型,并將數學模型計算出的遺留系數與試驗測定出的遺留系數進行比較。

數學模型的建立過程如下。為使推導過程的指代清晰明了,釷射氣(Tn)及其子體核素對應的下標見表1。由于使用了附壁狀態不受濕度影響的改進型閃爍室,可以假設生成的子體全部留在閃爍室內。

表1 釷射氣及其子體核素編號Table 2 The code list for thoron and its progenies

第n周期前后通道有效計數之差受到所有前序周期的影響,可以表述為

(1)

式(1)中,余量系數Rn-i+1表示第i周期的附壁釷子體遺留物對第n周期前后通道有效計數之差的影響,其計算公式為

(2)

式(2)中,Gn(i)表示第i周期遺留物在第n周期產生的計數。接下來需確定R的變化規律,將逐周期分析。

2.1.1 起始周期余量系數計算

起始周期存在關系:

G前(1)-G后(1)=R1C前(1),

(3)

子體末態數量計算公式為

(4)

(5)

周期內總α衰變數量計算公式為

G前(1)-G后(1)=3VC前(1)T-N1(T)-N2(T)。

(6)

通過上述公式可算出R1,下面將計算其他余量系數。

2.1.2 任意周期余量系數計算

后續周期遺留物變化規律(t>T):

(7)

(8)

特別地

(9)

2.1.3 遺留系數計算

定義遺留系數(γi)=Ri/R1,以定量分析第i周期的余量對測量結果的影響。遺留系數(γ)與測量結果(C)的關系為

(10)

顯然,遺留系數僅受測量周期(T)與核素衰變常數(λ)的影響,每一周期的測量結果都可由本周期的計數與前序周期的測量結果計算得到。附壁釷子體的影響經過60 h(6倍半衰期)后可以忽略不計,因此用于長期監測時,若將測量周期設置為24 h,則僅需考慮前3個遺留系數即可。

測量方法:將原型機與釷射氣濃度恒定的釷室(平衡濃度約為127 500 Bq/m3,遠大于實驗室環境中的釷射氣濃度)閉環連接,采樣24 h(周期時長設置為24 h),得到此周期內前1個閃爍室產生的總通道計數。然后將原型機與環境空氣連接,與釷室斷開,從而確保此后通道一的計數完全由采樣周期產生的附壁釷子體產生,測量并計算與空氣連接后第0～24、24～48、48～72 h的遺留系數,結果見表2?？梢钥闯?二者偏差較小,當測量周期設置為24 h時,在原型機的算法中,可采用遺留系數試驗結果的算術平均值(0.189、0.039 6、0.008 58)來扣除前序周期對當前周期的影響。

表2 附壁子體遺留系數的理論計算與測量結果Table 2 The calculation and measurement results of remnant coefficients

表2中理論計算與試驗測量結果的偏差定義為

(11)

式中:σ—理論與實驗結果的偏差,無量綱;CT—理論計算結果;CE—試驗測量結果。

2.2 線性響應

為了檢驗原型機能否在較寬濃度區間內測量氡與釷射氣的濃度,檢驗了原型機的線性響應,標定了刻度系數。

試驗方法:將活度濃度為44 841、37 479、25 393、8 474、1 376 Bq/m3的氡室與活度濃度為106 244、72 579、35 825、14 597、7 411 Bq/m3的釷室(濃度均由Alpha Guard DF2000標定)分別與原型機閉環連接,測得不同濃度下原型機的通道計數率。測量結束后,將通道計數率與Alpha Guard DF2000(以下簡稱DF2000)測得的對應濃度進行線性擬合,通過線性相關系數判斷原型機的線性響應。若原型機線性響應較好,則進一步得到通道對核素的刻度系數。由刻度系數與遺留系數即可得到通過通道計數計算活度濃度的公式為

(12)

(13)

式中:CRn(i)—第i周期測得氡濃度,Bq/m3;CTn(i)—第i周期測得釷射氣濃度,Bq/m3;G1(i)—第i周期通道一的計數,c;G2(i)—第i周期通道二的計數,c;εR1—通道一對氡的刻度系數,c/(Bq/m3);εR2—通道二對氡的刻度系數,c/(Bq/m3);εT1—通道一對釷射氣的刻度系數,c/(Bq/m3);γ1—附壁釷子體對下一周期的遺留系數,無量綱;γ2—附壁釷子體對兩周期后的遺留系數,無量綱;γ3—附壁釷子體對三周期后的遺留系數,無量綱。

DF2000具有濕度修正與能量分辨功能,可以在高濕度環境下分別測量Tn與Rn的放射性濃度,能夠直接將母體核素與放射性子體區分開來并排除濕度變化的干擾,該儀器經中國計量科學研究院檢定合格,具有較高的精度與準確度。試驗結果見表3。

表3 線性響應的試驗結果Table 3 The experimental results of linear response

由于釷射氣在延遲器中已接近全部衰變,因此對釷射氣的線性刻度僅考慮釷射氣濃度與通道一計數率的關系。由多組結果所確定的通道線性相關系數為0.999,表明原型機在較寬濃度范圍內的測量結果與DF2000保持一致,線性響應良好。在此濃度范圍內、周期24 h情況下,通道一、二對氡的刻度系數分別為122.4 c/(Bq/m3)和99.4 c/(Bq/m3),通道一對釷射氣的刻度系數為20.1 c/(Bq/m3)。

3 氡釷射氣監測裝置現場測試

3.1 測試方案

南方某新材料有限公司稀土廠房板框過濾車間的環境空氣,具有高Tn濃度、低Rn濃度、高濕度、強酸性、強γ放射性的特點(隨工作狀況變化幅度較大),適于探測器性能的檢驗。

為了驗證原型機對氡與釷射氣長期甄別測量的準確性,將原型機與DF2000布置于板框過濾車間的同一位置,從2023年2月18日上午開始進行96 h的連續測試,每經過24 h,將2臺儀器在此周期內各自測量結果的算術平均值比較一次,評價原型機的現場測量準確性。

3.2 試驗結果與分析

定義原型機與DF2000測量偏差的計算公式為

(14)

ε=|1-σ|×100%,

(15)

式中:σ—原型機與DF2000測量結果的比值,無量綱;ε—原型機與DF2000測量結果的偏差,無量綱;CP—原型機的測量結果,Bq/m3;CD—DF2000的測量結果,Bq/m3。

現場96 h比對測試結果見表4。

表4 現場試驗結果Table 4 The results of on-site experiment

從表4可看出,即使在氡與釷射氣濃度之比極為懸殊的情況下,原型機與DF2000對Tn與Rn的測量結果偏差均小于10%?，F場試驗中,原型機與DF2000對釷射氣測量結果比值的平均值為0.975,對氡測量結果比值的平均值為0.990,可以認為原型機在高濕度、強酸性、高Tn濃度、低Rn濃度環境中,24 h周期下對Tn與Rn的測量準確性符合要求,在長期監測過程中,能夠排除前序周期遺留物的影響。

3.3 分析與討論

為了解原型機對Rn與Tn放射性濃度探測的準確性,將原型機的現場測量結果與其他氡釷射氣聯合探測方法進行比較,由于原型機的現場測試中僅有一個比較對象(DF2000),而DF2000相對真值存在一定的偏差,因此對原型機與真值的偏差進行計算。計算方法為

εP(Rn)=[1-(1-εD(Rn))×σAv(Rn)]×100%,

(16)

εP(Tn)=[1-(1-εD(Tn))×σAv(Tn)]×100%,

(17)

式中:εP(Rn)—原型機對氡測量結果與真值的偏差,無量綱;εP(Tn)—原型機對釷射氣測量結果與真值的偏差,無量綱;εD(Rn)—DF2000對氡測量結果與真值的偏差,無量綱;εD(Tn)—DF2000對釷射氣測量結果與真值的偏差,無量綱;σAv(Rn)—原型機與DF2000對氡測量結果比值的平均值,無量綱;σAv(Tn)—原型機與DF2000對釷射氣測量結果比值的平均值,無量綱。計算結果見表5。

表5 氡釷射氣聯合監測方法比較Table 5 Comparison of radon-thoron simultaneous monitoring methods

2種固體徑跡探測器(GSF、LD-P)對氡與釷射氣的測量值與真值的偏差均略大于原型機。2種半導體探測器(RAD 7、RTM-2100)測量結果受濕度的影響很大,除RAD 7對氡的探測結果(與真值的偏差為3%)比原型機更準確外,其他測量結果與真值的偏差均大于原型機?；趦啥斡嫈捣ǖ膶嶒炑b置,主要用于在氡與釷射氣混合環境中對2種核素的放射性濃度進行快速測量,其對釷射氣的測量結果比原型機更準確但對氡的測量結果偏差更大,不具備對附壁子體的修正能力。Alpha Guard DF2000在純氡、純釷或氡釷濃度之比適中的環境下,均能給出較準確的測量結果,且不受濕度變化的影響。

本研究所構造的原型機,具備在氡與釷射氣混合環境中分別測量2種核素活度濃度的功能,測量結果與真實值的偏差均小于10%;且能夠在長期連續監測中,排除濕度與前序周期釷子體遺留物的干擾。

4 結論

使用改進型ZnS閃爍室,基于延遲計數法構建了氡釷射氣連續監測方法,在2個閃爍室之間增加延遲器使釷射氣充分衰變,分別得到由氡與釷射氣導致的計數,實現了在氡與釷射氣混合環境中甄別測量2種核素放射性濃度的目標。

在氡與釷射氣濃度之比較為懸殊的環境中,原型機依然能夠進行甄別測量。原型機測量準確度受濕度影響小,能夠排除前序周期子體遺留物的影響,且長期計數的穩定性高。