應(yīng)用蒙特卡羅方法表征HPGe探測器

牟婉君，李 梅，蹇 源

（中國工程物理研究院 核物理與化學(xué)研究所，四川 綿陽 621900）

探測器的探測效率是表征γ射線照射量率與探測器輸出脈沖之間的一重要物理量，在實際測量中，需針對不同的測量條件，如源的形式、幾何位置、環(huán)境測量等，通過對標(biāo)準(zhǔn)樣品的測量得到探測器在此探測條件下的探測效率，這種方法常常昂貴又費時。隨著計算機(jī)技術(shù)的不斷發(fā)展，蒙特卡羅模擬在計算探測器對光子的探測效率方面大量使用，該方法可實現(xiàn)樣品的無源效率刻度，具有節(jié)約、簡便及準(zhǔn)確度高等特點。但探測器的生產(chǎn)廠家難以給出精確的晶體幾何參數(shù)，同一批產(chǎn)品各探測器間的幾何參數(shù)差異也較大，這些差異對測量結(jié)果的分析有較大影響。為了使蒙特卡羅計算的結(jié)果準(zhǔn)確，需先對計算模型的關(guān)鍵尺寸進(jìn)行調(diào)整，一般廠家提供的晶體直徑與晶體長度宏觀可測，因此，認(rèn)為這兩項數(shù)值準(zhǔn)確合理，不需調(diào)整［1－5］。本工作通過實驗測量效率與蒙特卡羅方法模擬效率相結(jié)合的方法，準(zhǔn)確確定HPGe探測器晶體的死層厚度及冷指尺寸，并進(jìn)行實驗驗證。

1 材料與方法

1.1 探測器與計算模型

美國ORTEC公司生產(chǎn)的同軸型高純鍺γ譜儀，相對探測效率70%，分辨率（FHWM）1.85keV（60Co，1.33MeV），廠家提供的探測器幾何結(jié)構(gòu)及尺寸如下：探測器的晶體尺寸為φ2.96cm×6.77cm，死層厚度為0.09cm，測量的γ射線點源為241Am、152Eu、60Co、137Cs，點源與探測器的軸向距離為20cm，可有效減少符合相加的影響及源位置引起的誤差。

1.2 蒙特卡羅方法

實驗中采用蒙特卡羅計算程序MCNP 4C，該軟件是一多功能的蒙特卡羅計算程序，主要可計算中子、電子和光子3種粒子的運輸問題，同時包括相互之間的耦合輸運等，它具有較強(qiáng)的通用性和幾何處理能力、使用精確的點截面參數(shù)、適用多種問題的計算及豐富的降低方差技巧等特點，可很好用于跟蹤計算、物理實驗?zāi)M及輻射防護(hù)。

1.3 效率刻度

將4種標(biāo)準(zhǔn)源分別在γ譜儀上進(jìn)行測量（計數(shù)統(tǒng)計誤差小于0.5%），全能峰的凈面積由GrammaVision32軟件獲得。實驗測量中，樣品測量的幾何條件固定，但受實驗條件的限制，未做樣品的自吸收效正。按能量分布，采用下式［6］計算全能峰的探測效率：

式中：εp（Ei）為譜儀對能量為Ei的γ射線的探測效率；Na（Ei）是能量為Ei的γ射線在譜儀上觀測到的全能峰面積（凈計數(shù)率，s－1）；A0（Ei）為加入標(biāo)準(zhǔn)源能量為Ei的γ核素活度；P（Ei）為分支比。

2 實驗結(jié)果

2.1 實驗效率與探測效率的比較

采用HPGe探測器對標(biāo)準(zhǔn)241Am、152Eu、60Co、137Cs源在距探測器軸向距離為20cm處分別進(jìn)行測量，利用式（1）得到軸向距離的全能峰探測效率，源的實驗全能峰效率與計算效率的對比示于圖1。由圖1可知，在MCNP模擬計算中確定源與探測器的軸向距離20cm，按廠家提供的的死層參數(shù)0.09cm、冷指半徑0.28cm和冷指長度1.75cm進(jìn)行模擬計算，計算效率低于實測效率，這除了與計算中所取的放射源尺寸以及源到探測器距離和實際值之間存在偏差有關(guān)外，源本身的吸收和散射更是造成計算值與實測值出現(xiàn)偏差的原因，而產(chǎn)生這種偏差的主要原因來自探測器晶體的死層厚度和冷指尺寸。

圖1 20cm處計算效率與實驗效率的對比Fig.1 Comparison of simulated and measured counting efficiencies in 20cm detection height

2.2 死層厚度的調(diào)整

探測器晶體死層厚度對高、低能部分的γ射線的探測效率均有影響，但對低能γ射線的影響顯著，所以，實驗中主要通過低能部分計算效率與實測效率的對比和分析來找到兩者相符合的死層厚度。當(dāng)計算效率大于實測效率時，需增加死層厚度；反之，則需減小死層厚度。

由圖1知，源的計算效率低于實測效率，則應(yīng)在模擬計算中減小死層厚度。模擬了探測器對能量為59.54、81.01及121.78keV 的低能光子的探測效率，并與實測效率進(jìn)行了比較，從而對晶體死層厚度進(jìn)行調(diào)整。計算的探測器探測效率隨死層厚度的變化情況及與實測結(jié)果的相對偏差列于表1。由表1可知，死層厚度設(shè)在0.07cm時，低能部分計算效率與實測效率相對偏差在5%以內(nèi)，符合效果很好。圖2為死層厚度為0.07cm時源的計算效率與實測效率的對比，從圖中可看出，在低能區(qū)計算效率與實測效率符合較好，但在高能區(qū)兩者有一定的偏差。

表1 不同死層厚度下計算效率與實驗效率的相對偏差Table 1 Deviation between simulated and measured counting efficiencies with different dead layers thickness

圖2 死層厚度為0.07cm時計算效率與實驗效率的對比Fig.2 Comparison of simulated and measured counting efficiencies with dead layers thickness of 0.07cm

2.3 冷指尺寸的調(diào)整

冷指處于探測器的后端，對穿透力較強(qiáng)的高能γ射線的影響相對較大，在上述確定死層厚度0.07cm的基礎(chǔ)上，實驗中利用152Eu、60Co源的高能部分的強(qiáng)峰來觀察冷指尺寸變化對模擬計算效率的影響，當(dāng)計算效率大于實測效率時，增大冷指尺寸；反之，減小冷指尺寸。從圖2可看出，調(diào)整死層厚度后，在高能區(qū)計算效率大于實測效率，需通過調(diào)整冷指尺寸來使兩者相符合，實驗中對冷指尺寸的調(diào)整主要針對冷指的長度。

模擬了在不同冷指長度下探測器對能量為964.0、1 173.2、1 332.6及1 408.1keV 的高能光子的探測效率，計算效率與實測效率的相對偏差列于表2。由表2可知，當(dāng)冷指長度設(shè)置為3.25cm時，高能部分計算效率與實測效率相對偏差在5%以內(nèi)，兩者符合較好。圖3為冷指長度為3.25cm時源的計算效率與實測效率的對比。從圖中可看出，通過調(diào)整死層厚度與冷指尺寸，在59.54～1 408.1keV光子能量范圍內(nèi)計算效率與實測效率符合較好。

表2 不同冷指長度下計算效率與實驗效率的相對偏差Table 2 Deviation between simulated and measured counting efficiencies with different sizes of cold refer

圖3 冷指長度為3.25cm時計算效率與實驗效率的對比Fig.3 Comparison of simulated and measured counting efficiencies with cold refer length of 3.25cm

2.4 結(jié)果驗證

通過上述實驗確定的Monte Carlo方法模擬計算光子探測效率的最佳模型參數(shù)為：晶體死層厚度0.07cm、冷指半徑0.28cm、冷指長度3.25cm。基于此參數(shù)，對以上結(jié)論進(jìn)行了驗證，采用MCNP程序模擬了源距探測器軸向距離25cm下不同徑向距離的探測效率，并將計算效率與實驗效率進(jìn)行了對比，結(jié)果如圖4所示。由圖4可知，在軸向距離25cm下，不同徑向距離處計算效率與實驗效率的低能與高能部分符合較好，證明修改參數(shù)后的模型對模擬計算是合理的。

圖4 源距探測器不同徑向距離時計算效率與實驗效率的對比Fig.4 Comparison of simulated and measured counting efficiencies in different radial distances

3 結(jié)論

采用MCNP能準(zhǔn)確計算高純鍺探測器的探測效率，晶體的死層厚度對低能光子的探測效率影響較大，而冷指尺寸對高能光子探測效率影響較大，在實驗中通過反復(fù)實驗驗證獲得了最佳計算模型參數(shù)，計算效率與實測效率相對偏差在5%以內(nèi)，這為以后開展無源刻度奠定了基礎(chǔ)。

在本實驗分析過程中忽略了樣品自吸收的影響，在測量過程中須保持所有點源樣品到探測器軸線距離一致，所得數(shù)據(jù)才具可靠性。在測量過程中的主要誤差來源有全能峰凈計數(shù)率的統(tǒng)計漲落和樣品與探頭距離的精密控制。

［1]許淑艷.蒙特卡羅方法在實驗核物理中的應(yīng)用［M］.北京：原子能出版社，1996.

［2]向東，郭南英.γ能譜的蒙特卡羅［J］.數(shù)理醫(yī)藥雜志，2006，19（3）：77－79.XIANG Dong，GUO Nanying.Monte Carlo simulation of gamma spectrum［J］.Journal of Mathematical Medicine，2006，19（3）：77－79（in Chinese）.

［3]張斌全，馬吉增，程建平，等.蒙特卡羅方法計算用于低能光子測量的高純鍺探測器的效率［J］.核電子學(xué)與探測技術(shù)，2005，25（3）：274－277.ZHANG Binquan，MA Jizeng，CHENG Jianping，et al.Calculation of the detection efficiency of an HPGe detector in low energy photon measurement with Monte Carlo method［J］.Nucl Electron Detect Technol，2005，25（3）：274－277（in Chinese）.

［4]MONES L，DEDONDER J，XILEI L，et al.Calculation of the absolute peak efficiency of gamma－ray detectors for different counting geometries［J］.Nucl Instrum Methods，1981，187：451－472.

［5]MONES L，HOSTE J.Calculation of the peak efficiency of high－purity germanium detectors［J］.Int J Appl Radiat Isot，1983，34 （8）：1 085－1 095.

［6]向長欣.一組實用γ譜分析效率刻度曲線［J］.核電子學(xué)與探測技術(shù)，1994，14（6）：363－366.XIANG Changxin.A group of applied efficiency calibration curves for γ－spectrometric analysis［J］.Nucl Electron Detect Technol，1994，14（6）：363－366（in Chinese）.