李春娟 張偉華 肖雪夫 王志強(qiáng) 劉毅娜

(中國(guó)原子能科學(xué)研究院,計(jì)量與校準(zhǔn)技術(shù)重點(diǎn)實(shí)驗(yàn)室,北京102413)

GEM-TEPC電場(chǎng)模擬及增益分析方法研究

李春娟 張偉華 肖雪夫 王志強(qiáng) 劉毅娜

(中國(guó)原子能科學(xué)研究院,計(jì)量與校準(zhǔn)技術(shù)重點(diǎn)實(shí)驗(yàn)室,北京102413)

為對(duì)現(xiàn)有組織等效正比計(jì)數(shù)器(TEPC)進(jìn)行改進(jìn),使其滿足高注量率輻射場(chǎng)測(cè)量,本實(shí)驗(yàn)室研制一套基于氣體電子倍增技術(shù)(GEM)的寬量程TEPC。對(duì)GEM-TEPC來說,探測(cè)器的增益是其關(guān)鍵指標(biāo),施加在GEM電極上的電壓、漂移區(qū)場(chǎng)強(qiáng)、收集區(qū)場(chǎng)強(qiáng)、靈敏區(qū)充入氣體的成分、氣壓均影響探測(cè)器的有效增益。本工作通過ANSYS程序?qū)EM-TEPC的電場(chǎng)進(jìn)行模擬計(jì)算,并對(duì)利用Garfield程序?qū)ζ湓鲆嫣匦赃M(jìn)行模擬的方法以及增益測(cè)量方法進(jìn)行研究。

GEM-TEPC 電場(chǎng) 增益 ANSYS Garfield

1 引 言

微劑量學(xué)實(shí)驗(yàn)研究方法的基本出發(fā)點(diǎn)是使用一個(gè)組織等效氣體體積來模擬一個(gè)正常密度的小生物組織體積。H.H.Rossi最早提出并使用正比計(jì)數(shù)器測(cè)量微劑量參數(shù)[1]。用于微劑量學(xué)量測(cè)量時(shí),正比計(jì)數(shù)器的壁材料常選用組織等效塑料(A-150),空腔中充入組織等效氣體,這一類正比計(jì)數(shù)器被稱為組織等效正比計(jì)數(shù)器(TEPC)。

氣體電子倍增技術(shù)(GEM)是20世紀(jì)90年代,由歐洲核子中心的F.Sauli等人發(fā)明的一種微結(jié)構(gòu)氣體探測(cè)器[2],該探測(cè)器具有高計(jì)數(shù)率、抗輻照、信號(hào)讀出模式多樣等優(yōu)點(diǎn)。由于在高注量率輻射場(chǎng)測(cè)量時(shí),常規(guī)TEPC,存在信號(hào)堆積現(xiàn)象,為解決這一問題,拓展現(xiàn)有TEPC的注量率測(cè)量范圍,原子能院計(jì)量測(cè)試部張偉華等人將GEM應(yīng)用于TEPC,制作成GEM-TEPC。

對(duì)GEM-TEPC來說,探測(cè)器的增益是其關(guān)鍵指標(biāo)。GEM膜上下面的電壓差對(duì)形成GEM的有效增益很重要,故需要深入研究GEM的電場(chǎng)特性。另外、漂移區(qū)場(chǎng)強(qiáng)、收集區(qū)場(chǎng)強(qiáng),靈敏區(qū)充入氣體的成分、氣壓均影響探測(cè)器的有效增益。本工作通過ANSYS[3]程序?qū)EM-TEPC的電場(chǎng)進(jìn)行模擬計(jì)算,并對(duì)利用Garfield[4]程序?qū)ζ湓鲆嫣匦赃M(jìn)行模擬的方法以及增益測(cè)量方法進(jìn)行研究。

2 裝置介紹

基于GEM的氣體探測(cè)器一般主要由漂移電極、GEM膜和PCB(印刷電路板)讀出電極三層組成,如圖1所示。標(biāo)準(zhǔn)GEM膜是雙面覆銅且雙錐形通孔呈均勻分布的聚酰亞胺(Kapton)膜,如圖2所示。

其中Kapton膜厚度為50μm,銅箔厚度為5μm,雙錐形孔間距140μm,孔表面直徑70μm以及孔心直徑50μm。目前使用的GEM膜大多是從歐洲核子中心(CERN)購(gòu)買,也有個(gè)別實(shí)驗(yàn)室自制厚層GEM膜。當(dāng)在上下兩層銅箔上加一定電壓,會(huì)在孔中形成很強(qiáng)的電場(chǎng)(＞100kV/cm),電子穿過小孔時(shí)在強(qiáng)電場(chǎng)的驅(qū)動(dòng)下會(huì)在工作氣體中倍增放大,放大后的電子在收集區(qū)的電場(chǎng)作用下繼續(xù)向下漂移,最后被PCB收集。

將GEM應(yīng)用于TEPC,其結(jié)構(gòu)主要包括A150陰極板、GEM膜和,陽極板,用不同厚度的絕緣框隔出漂移區(qū)、收集區(qū)的距離。可通過電阻鏈分壓供電的方式在漂移電極、GEM上下銅箔和讀出電極上分別加上不同的電壓,如圖3所示。分壓電阻值的設(shè)置需考慮使探測(cè)器獲得合適的增益所需要的漂移區(qū)場(chǎng)強(qiáng)、收集區(qū)場(chǎng)強(qiáng)及GEM膜壓差。

為保證探測(cè)器的組織等效性,漂移電極選用國(guó)際通用的A150組織等效塑料,壁厚滿足帶電粒子平衡要求。探測(cè)器內(nèi)充組織等效氣體,組織等效氣體的選擇須考慮到組織成分的性質(zhì)、氣體增益性質(zhì),以及氣體成分和壁材料成分是否匹配。目前組織等效氣體主要選用兩種混合氣體,一種混合氣體以甲烷為主(64.4%CH4,32.5%CO2,3.1%N2);另外一種混合氣體以丙烷為主(55%C3H8,39.6%CO2, 5.4%N2)。經(jīng)過前期實(shí)驗(yàn),基于丙烷的混合氣體增益更高,選用基于丙烷的混合氣體。

3 模擬軟件平臺(tái)介紹

由CERN開發(fā)研制的Garfield程序主要用來模擬多絲正比室、電離室等傳統(tǒng)型氣體探測(cè)器中的自由電子、離子與氣體分子的相互作用。其主要運(yùn)用內(nèi)部的各種程序包模擬氣體探測(cè)器的性能,每個(gè)程序包模塊都有其特定的功能,如:CELL確定探測(cè)器的幾何結(jié)構(gòu),并且可導(dǎo)入ANSYS軟件模擬計(jì)算的三維電場(chǎng)值;GAS計(jì)算氣體的屬性,包括湯姆森系數(shù)、吸附系數(shù)、橫向和縱向擴(kuò)散系數(shù)、漂移速度隨電場(chǎng)強(qiáng)度變化等情況;FIELD顯示探測(cè)器的電場(chǎng)強(qiáng)度分布,而電勢(shì)與電場(chǎng)圖與探測(cè)器的參數(shù)設(shè)置有關(guān);DRIFT模擬自由電荷的漂移雪崩;SIGNAL計(jì)算電極的感應(yīng)電荷的大小等等。由于Garfield程序本身無法進(jìn)行探測(cè)器的三維結(jié)構(gòu)的構(gòu)建,不能進(jìn)行電場(chǎng)強(qiáng)度的模擬計(jì)算,在模擬電子漂移雪崩過程中,需要借助第三方軟件來計(jì)算探測(cè)器內(nèi)部電場(chǎng)強(qiáng)度的分布,常用的模擬計(jì)算電場(chǎng)強(qiáng)度的軟件有ANSYS與Maxwell兩種。本工作的電場(chǎng)場(chǎng)強(qiáng)均是由ANSYS軟件計(jì)算。

ANSYS軟件是融結(jié)構(gòu)、流體、電磁場(chǎng)、聲場(chǎng)和耦合場(chǎng)于一體的大型通用有限元分析軟件,由世界上最大的有限元分析軟件公司之一的美國(guó)ANSYS公司開發(fā),它能與多數(shù)CAD軟件接口,實(shí)現(xiàn)數(shù)據(jù)的共享和交換。軟件主要包括三個(gè)部分:前處理模塊,分析計(jì)算模塊和后處理模塊。

4 電場(chǎng)計(jì)算方法及結(jié)果

利用ANSYS軟件對(duì)電場(chǎng)進(jìn)行模擬計(jì)算主要包括三個(gè)步驟:

1)通過前處理模塊創(chuàng)建物理環(huán)境,主要包括過濾圖形界面、定義工作標(biāo)題、指定工作名、定義單元類型及材料屬性等;

2)通過分析計(jì)算模塊建立三維幾何模型,對(duì)每個(gè)幾何單元賦予單元類型及材料特性,劃分網(wǎng)格,加載電壓及周期性邊界條件,并對(duì)電場(chǎng)進(jìn)行分析計(jì)算;

3)通過后處理模塊進(jìn)行畫圖,并導(dǎo)出數(shù)據(jù)結(jié)果。

由于GEM膜是周期性的重復(fù)單元,在建立三維幾何模型時(shí)可以選取一個(gè)最小單元來構(gòu)建。由于標(biāo)準(zhǔn)的GEM膜孔是雙錐形,在構(gòu)建過程中,需要四個(gè)部件才能構(gòu)建一個(gè)GEM膜孔,即兩個(gè)圓錐體以及兩個(gè)圓筒,根據(jù)所選擇的最小單元模塊,用同樣的方式構(gòu)建其他的四個(gè)1/4孔,即可以在ANSYS軟件中構(gòu)建出最小單元模型。

計(jì)算中建立的三維幾何模型如圖4所示,網(wǎng)格劃分情況如圖5所示,在漂移區(qū)、GEM膜以及收集區(qū)的電壓分別為100V,400V,100V時(shí),電場(chǎng)及電勢(shì)分布計(jì)算結(jié)果如圖6至圖7所示。可以看出,由于GEM微孔通道直徑很小,漂移電極和讀出電極之間的電力線在GEM微孔中密集產(chǎn)生高強(qiáng)度雙極電場(chǎng)。

5 增益計(jì)算方法

將ANSYS計(jì)算的最小單元的電場(chǎng)強(qiáng)度計(jì)算結(jié)果導(dǎo)入Garfield軟件,通過x-y平面周期平鋪的方法構(gòu)建整個(gè)標(biāo)準(zhǔn)GEM膜的電場(chǎng)及電勢(shì)分布,在Garfield軟件平臺(tái)上依據(jù)CELL,GAS,FIELD以及DRIFT這幾個(gè)程序包來進(jìn)一步計(jì)算氣體的性質(zhì),模擬帶電粒子在探測(cè)器的電場(chǎng)作用下產(chǎn)生的電離電荷并追蹤其雪崩和漂移的過程,模擬計(jì)算GEM-TEPC的增益特性。

GEM探測(cè)器的增益主要是由入孔系數(shù)、絕對(duì)增益、出孔系數(shù)三者共同決定的,而GEM膜的絕對(duì)增益隨著GEM膜的上下兩端電壓差的增加而增大。GEM膜的電子入孔系數(shù)是由漂移電場(chǎng)與GEM膜壓差的比值決定的,在滿足GEM膜絕對(duì)增益的條件下,可以固定GEM膜的壓差,進(jìn)而根據(jù)電子的入孔系數(shù)得到漂移電場(chǎng)的優(yōu)化值。

GEM膜的出孔系數(shù)與收集區(qū)電場(chǎng)和GEM膜壓差有關(guān),并且在很大程度上影響了孔中雪崩電子被陽極吸收的概率。在GEM孔內(nèi)的一定范圍內(nèi)隨機(jī)給定電子,統(tǒng)計(jì)不同收集區(qū)電場(chǎng)和GEM膜壓差下電子的出孔系數(shù)。電子出孔系數(shù)通常不能達(dá)到100%,這主要是因?yàn)橛胁糠蛛娮颖籊EM膜下層的銅膜俘獲。

靈敏區(qū)充入氣體的成分、氣壓也會(huì)影響探測(cè)器的增益。探測(cè)器內(nèi)充以丙烷為主(55%C3H8, 39.6%CO2,5.4%N2)的組織等效氣體。實(shí)際使用時(shí),氣壓由擬模擬組織的尺度以及探測(cè)器的靈敏體積決定[8],一般為0.1atm～0.5atm,模擬計(jì)算不同氣壓(0.1atm～0.4atm)時(shí),GEM有效增益隨GEM各電極上的電壓的變化。

6 探測(cè)器增益測(cè)量方法

探測(cè)器漂移電極上設(shè)置的校準(zhǔn)腔室置有241Am α源,利用其對(duì)探測(cè)器的增益進(jìn)行測(cè)量,具體測(cè)量方法如下:設(shè)漂移區(qū)原初電離產(chǎn)生電子個(gè)數(shù)為Q1,到達(dá)陽極板收集到的電子個(gè)數(shù)為Q2,經(jīng)過前放放大后測(cè)量到的脈沖信號(hào)幅度為V,Eα為α粒子在漂移區(qū)內(nèi)的最大能量沉積,可由SRIM計(jì)算得到(在2μm軟組織中的能量沉積為266keV),α粒子在組織等效氣體中的電離功W約為(29～31)eV,計(jì)算中選取30eV,根據(jù)前放142PC的使用說明書,其反饋電容為0.15pF,則Q1由公式(1)計(jì)算得到,Q2由公式(2)計(jì)算得到,增益G由公式(3)計(jì)算得到。

通過實(shí)驗(yàn)測(cè)量及模擬計(jì)算相結(jié)合的方法,對(duì)探測(cè)器增益進(jìn)行分析,可給出漂移區(qū)與收集區(qū)場(chǎng)強(qiáng)的參考值范圍,為探測(cè)器分壓電路的設(shè)置提供優(yōu)化參數(shù)。

7 結(jié)束語

綜上,通過ANSYS和Garfield程序相結(jié)合的方法對(duì)GEM-TEPC的電場(chǎng)及增益特性進(jìn)行模擬研究是切實(shí)可行的,與通過實(shí)驗(yàn)對(duì)增益進(jìn)行測(cè)量相結(jié)合,可給出漂移區(qū)與收集區(qū)場(chǎng)強(qiáng)的參考值范圍,為探測(cè)器分壓電路的設(shè)置提供優(yōu)化參數(shù),給出探測(cè)器合理的高壓供電范圍。本工作通過模擬計(jì)算給出了探測(cè)器電場(chǎng)分布,并對(duì)通過Garfield程序?qū)μ綔y(cè)器增益進(jìn)行計(jì)算的方法進(jìn)行研究,對(duì)增益的實(shí)驗(yàn)測(cè)量方法進(jìn)行研究,為下一步對(duì)探測(cè)器增益進(jìn)一步精細(xì)設(shè)計(jì)計(jì)算及實(shí)驗(yàn)測(cè)量,對(duì)探測(cè)器工作參數(shù)進(jìn)行優(yōu)化分析奠定了基礎(chǔ)。

[1] Rossi H H,Zaider M.Microdosimetry and its applications [M].Springer,1996.

[2] SAULI F.GEM:Nucl Instr Meth A[M].1997:386～531.

[3] http://www.Ansys.com/[DB/OL].

[4] http://garfield.web.cern.ch/garfield/[DB/OL].

[5] 張偉華,王志強(qiáng).組織等效正比計(jì)數(shù)器的測(cè)量原理和方法[J].劑量研究,2008(3):45～54.

Simulation of Three Dimension Distribufion of Electric Field in GEM-TEPC and Methods to Study of its Effective Gain

LI Chun-juan ZHANG Wei-hua XIAO Xue-fu WANG Zhi-qiang LIU Yi-na
(National Key Laboratory for Metrology and Calibration Techniques,China Institute of Atomic Energy,Beijing 102413,China)

To optimize the design of tissue-equivalent proportional counter(TEPC)and make it satisfy the measurement of radiation field with high fluence rate,we set up a TEPC based on a Gas Electron Multiplier(GEM).The effective gain is very important for GEM-TEPC,and the key factors affecting the effective gain are the voltage of GEM electrode,the drift area field,collected area field,the gas components filling into sensitive area and the gas pressure.Three dimension distribufion of electric field in GEM-TEPC were simulated with code ANSYS,and the way how to simulate the maximum safe gain of GEM-TEPC with code Garfield and the way how to measure the maximum safe gain were studied.

GEM-TEPC Electric field Effective gain ANSYS Garfield

1000-7202(2017)01-0072-04

O571.5

A

2017-02-17,

2017-02-24

李春娟(1981-),女,副研究員,碩士,主要研究方向:電離輻射計(jì)量技術(shù)。