核工程與核技術畢業(yè)設計（論文）閃爍探測器探測效率與γ射線能量關系的數(shù)值模擬

1、 畢業(yè)設計（論文） 題 目 閃爍探測器探測效率與射線能量關系的數(shù)值模擬英文題the numerical simulation relationship of scintillation detector detection efficiency and -ray energy學生姓名: 專 業(yè): 核工程與核技術班 級: 指導教師: 二零一一年六月摘 要閃爍能譜儀的探測效率與射線能量關系之間的關系，通過使用蒙特卡羅的數(shù)值模擬方法進行研究。nai閃爍探測器由于其探測效率高、在常溫下可以使用的特點，越來越受到人們的重視。因此研究其探測效率對進一步的開發(fā)和研究閃爍探測器是很有必要的，特別是閃爍探測器與

2、能量之間的關系決定了今后的研究方向。射線主要分為兩大類：一類是，自然射線；另一類是，人工射線，用nai探測器了對射線進行探測。然后用mcnp4c軟件進行模擬實驗，得到了不同射線能量下的能譜。通過能譜圖我們可以根據(jù)全能峰內的相對計數(shù)率，得出相對探測效率。再得出相對探測效率以后，根據(jù)相對應的射線能量，擬合出圖形，得到探測效率與射線能量關系的數(shù)值模擬。通過對上述內容的研究，本課題基本上達到了預期研究目的，取得較好的階段性成果，對閃爍探測器探測效率與能量關系的研究起到一定的科學指導作用，為實驗提供了一定的依據(jù)。關鍵詞：探測效率； 射線； 蒙特卡羅模擬abstractthe relationship b

3、etween detection efficiency of the -ray scintillation spectrometer and the -ray energy, we are using monte carlo simulation method to study. because of its high detection efficiency, can be used at room temperature characteristics, nai scintillation detector more attention has been paid. 顯示對應的拉丁字符的拼

4、音therefore, it is very necessary to study the detection efficiency for further development and research, especially the relationship between scintillation detectors and -ray energy direction of our future research. -ray can be divided into two categories: one is the natural -ray; the other is artifi

5、cial -ray, with the nai detector for -ray detection. then simulation experiments by mcnp4c software. obtained under different -ray energy spectrum of -ray. we obtained the relative detection efficiency by the energy spectra we can according to all-powerful within the relative peak count rate. and th

6、en come after the relative detection efficiency, according to the corresponding -ray energy, fitting out graphics, and -ray detection efficiency has been the relationship between energy simulation. through the above research, this research project is basically to achieve the desired objectives and a

7、chieve good initial results, the efficiency of the scintillation detector and -energy relationship of scientific research play a guiding role for the experiment provided some basis.key words: detection efficiency; -rays; monte carlo simulation目 錄緒 論11.1 選題意義及目的11.2 實驗背景12 模擬實驗基礎32.1 射線與物質的相互作用32.1.1

8、 光電效應32.1.2 康普頓效應32.1.3 光電子對效應42.2 閃爍能譜儀52.2.1 閃爍探測器的原理62.2.2 閃爍探測效率73 構建模型83.1 蒙特卡羅方法及mcnp程序簡介83.1.1 蒙特卡羅計算方法及其應用范圍83.1.2 mcnp程序簡介83.2 建立數(shù)學模型113.2.1 點源模擬113.2.2 探測器柵元卡片填寫113.2.3 模擬模型143.2.4 源定義143.2.5 輸入卡片154 計算結果與處理174.1 模擬輸出結果174.1.1 自然射線的模擬輸出結果174.1.2 人工射線的模擬輸出結果204.2 數(shù)據(jù)處理與結果23結論26致 謝27參考文獻28緒 論

9、1.1 選題意義及目的隨著我國經(jīng)濟的發(fā)展，國民生活水平不斷提高，國內需求不斷增長。核輻射探測器也逐漸被使用到各行各業(yè)，在各類探測器中，閃爍體探測器由于其對射線的探測分辨時間短、探測效率高、能測量射線的能量的優(yōu)點，所以是目前應用的最廣的射線探測器。因此研究閃爍探測器探測效率與射線能量關系對我們提高閃爍探測器探測效率有著重要的意義。本課題主要是應用蒙特卡羅的方法對探測器探測效率與射線能量關系進行數(shù)值模擬，通過這個數(shù)值關系，我們可以找到最佳的射線能量對應著最高的探測效率。這樣就對我們以后的實驗、研究提供了一定的依據(jù)。1.2 實驗背景蒙特卡羅方法又名隨機模擬法或統(tǒng)計實驗法，它是一種具有獨特風格的數(shù)值計

10、算方法，它既能求解定性的數(shù)學問題，也能求解隨即問題。半個多世紀以來，由于科學技術的發(fā)展和電子計算機的發(fā)明，這種方法作為一種獨立的方法被提出來，但與一般數(shù)值計算方法有很大區(qū)別，它是以概率統(tǒng)計理論為基礎的一種方法。由于蒙特卡羅方法能夠比較逼真地描述事物的特點及物理實驗過程，解決一些數(shù)值方法難以解決的問題，因而該方法的應用領域日趨廣泛1。蒙特卡洛方法由于其簡單性、靈活性和普遍性在獲得廣泛應用。尤其在核物理上更是得到很大應用，20世紀80年代，schlumberger、lanl和美國北卡羅來納州立大學開始將蒙特卡羅方法用于核測井的模擬，取得很大的成功2。以后的發(fā)展里，隨著計算機技術的不斷進步，現(xiàn)代計算

11、機技術使隨即采樣成為解決復雜問題的實用方法。到了20世紀90年代，蒙特卡羅模擬已經(jīng)成為核測井科研和生產中的重要工具。近一些年來，一些大型的蒙特卡羅方法應用軟件先后生產并發(fā)展，例如，美國的橡樹嶺國家實驗室的輻射屏蔽信息中心rsic、阿貢國家實驗的軟件中心nesc和計算機物理庫cpc等所收集的大量應用軟件，就有相當數(shù)量的蒙特卡洛方法應用軟件。具體的用在核技術領域的軟件有美國橡樹嶺國家實驗室的mose程序、美國stanford linear accelerator center的egs程序、美國生地國家實驗室的sandyl和tiger程序、美國洛斯阿拉莫斯國家實驗室的mcnp程序、歐洲核子研究中心的

12、gernt程序、英國的mcbend程序以及法國的tripoli程序等。這些程序在粒子輸運模擬中發(fā)揮了重要作用，但這些程序一般為通用程序，并且操作復雜不便計算和探測問題。因而，研究核探測問題的蒙特卡羅計算方法及軟件成為當今輻射輸運領域的熱門課題2,3。在運用蒙特卡羅方法對射線的響應函數(shù)進行計算以及探測效率刻度方面，國內外重20世紀60年代就已經(jīng)有很多人開始研究工作，而從20世紀70年代以來，國內外有關這方面的研究工作可分為：對裸晶體，隨著光子能量越來越高，考慮的反應機制日趨完善；對低能光子，在反應機制較簡單的情況下，研究了體源、反射層及源與探測器的幾何關系對響應函數(shù)及探測效率的影響2,3。進入2

13、0世紀90年代以后，國內外持續(xù)有學者對蒙特卡洛方法計算響應函數(shù)和探測效率的方法進行改進完善，像h.vnguyen等利用ngrc程序、fayezh.h等利用egs4程序模擬了nai(tl)閃爍探測器對射線的響應函數(shù)，清華大學的hu-xia-shi等人利用mcnp模擬了nai(tl) 閃爍探測器對射線的響應函數(shù)3,4?？v觀前人的研究成果，可以發(fā)現(xiàn)蒙特卡洛方法在計算射線在探測器中與探測器效率刻度方面，為我們研究提供了極大的便利，并可以減少大量的實驗工作，但是在大多停留在實用國外的通用軟件上，而國內還沒有比較成功的自主研發(fā)陳鼓工我們方便使用，并且在探測效率刻度大多是正對源與探測器的空間幾何位置來進行，

14、很少有人對射線能量與探測效率之間的關系進行研究3,4。2 模擬實驗基礎2.1 射線與物質的相互作用2.1.1 光電效應當射線與原子殼層軌道上的內層電子碰撞時，將所有能量交給殼層的一個內層電子，電子克服了電離能，脫離原子而運動，產生光電子，而射線被完全吸收，這種作用稱為光電效應。為使光電效應保持動量守恒，必須有第三者，即反沖核。由于自由電子不能吸收射線而形成光電子，因此光電效應產生的幾率將隨著電子結合能的增加而增加，靠近k層、l層產生光電效應的幾率最大1,2。圖21 光電效應示意圖當入射射線的能量大于原子的電離能，才能產生該物質的光電效應。光電子的動能是入射射線能量與電離能之差。一般認為，光電效

15、應視整個原子看成一個統(tǒng)一的體系與入射光子發(fā)生作用，吸收能量，發(fā)射光電子。沒有足夠的能量克服軌道電子的結合能，電子就不可能脫離原子的束縛而成為自由電子1,2。2.1.2 康普頓效應 當入射的射線與一個自由電子發(fā)生一次碰撞時，射線將部分能量傳給電子，使它與射線的初始運動方向呈角射出，而射線則與初始運動方向呈角散射。這種現(xiàn)象稱為康普頓吳有訓效應，簡稱康普頓效應1,2。將入射光子看成是具有一定能量和動量的粒子與電子發(fā)生剛性碰撞。在這一碰撞過程中，入射光子將一部分能量傳遞給電子。而光子發(fā)生能量和方向的改變，形成能量為，動量為的光量子。散射角為。電子在獲得一定能量之后，克服結合能，沿方向出射。稱為反沖電子

16、或康普頓電子。其動能取決與電子獲得能量與電子結合能之差；并具有相應的動量，作用過程見圖22所示。根據(jù)能量和動量守恒定律，可以得到沿方向散射的射線能量： （21）沿方向出射的反沖電子能量： （22） 式中：電子靜止質量; 光速。圖22 康普頓效應示意圖2.1.3 光電子對效應隨著入射光子能量的增高，光電效應的吸收作用很寬減弱，康普頓效應也逐漸減弱。當光子能量大于1.02mev時，就存在形成電子對效應的幾率，即光子完全被吸收而產生一正、負電子對。電子對的動能為 (2-3)式中，入射光仔能量； 正、負電子的靜止質量； 光速。正電子與負電子不同，它很不穩(wěn)定，其壽命大約為。正電子在物理中由于電力而損失能

17、量后，將和電子結合，通常轉化為兩個能來那個為0.51mev且方向相反的光子，這種現(xiàn)象即正電子湮沒1,2。當能量大于1.02mev的光子通過原子序數(shù)較大的物質時，將容易測到這種射線。圖23 電子對效應示意圖2.2 閃爍能譜儀閃爍體是透明的絕緣體。射線與閃爍體相互作用所產生的次級電子可以再和它作用，使其原子分子電離和激發(fā)。顯然，產生的電子離子是無法從物質中引出來的，但其過程中產生的閃光可從物質中透射出來，被收集轉換成點脈沖放大后，可供分析記錄。因此，閃爍體可作射線探測器的探測介質。閃爍體可分為有機和無機兩大類，有機閃爍體又分為液體閃爍體、有機晶體和塑料閃爍體三種。這些閃爍體對射線有中等大小的吸收本

18、領。有機閃爍體的探測器有比較高的探測效率。盡管其能量分辨率不好，但響應快、時間分辨好。因此，常用這類探測器用作快符合裝置上的探測器。無機閃爍體的典型代表是鉈激發(fā)的碘化鈉晶體和碘化銫晶體。他們原子序數(shù)高、密度大，對射線有較高的吸收本領。其時間響應快，能量分辨率也好，是20世紀50年代以來探測射線的主要探測器2,3。2.2.1 閃爍探測器的原理閃爍探測器主要由閃爍體、光電倍增管和相應的電子儀器贊歌主要部分組成。圖2-4是閃爍能譜的結構框圖，探頭部分主要包括閃爍體和光電倍增管和前置放大器；主機部分包括主放大器，多道脈沖幅度分析器，接口電路、電源電路以及計算機和相應的應用軟件等。圖2-4 閃爍能譜儀原

19、理框圖目前，閃爍晶體主要使用nai(tl)晶體，它的密度大（=3.67g/cm3）,而且高原子系數(shù)的碘（z=53）占重要的85%，所以對射線探測效率特別高，相對發(fā)光效率大，它的光譜最強波長為415nm左右，能與光電倍增管的光譜響應較好的一種。但是nai(tl)晶體容易潮解，吸收空氣水分會變質失效，所以一般nai(tl)晶體使用時都是裝在密封的金屬盒中。射線與nai(tl)晶體相互作用，會使其電離、激發(fā)而發(fā)射熒光，從閃爍體出來的光子通過光導射向光電倍增管的光陰極，由于光電效應，在光陰極上打出電子。光電子經(jīng)過各級之間的電壓加速，聚焦后射向打拿極。每個光子在打拿極上擊出幾個電子，這些電子射向后續(xù)的打

20、拿極。所以，最后射向陽極的電子數(shù)目是很多的，陰極把這些電子收集起來，轉變成電信號并輸出。經(jīng)過nai(tl)探測器和光電倍增管輸出的信號往往較弱，需要加前置放大器。其主要作用是：提高系統(tǒng)的信噪比；減小信號經(jīng)電纜傳送時外界干擾的影響；阻抗匹配，功率放大。探頭的信號經(jīng)前置放大器輸出，但前置放大器輸出脈沖幅度和波形并不適合后面的分析測量設備的要求，所以對信號還需要進一步放大和成形，在放大和成形的過程中必須嚴格保持探測器輸出的有用信號（如射線的能量和時間信息），可能減少它們的失真。這樣的放大和成形就由放大器來完成。放大器輸出的信號要適應分析測量設備的要求，必須解決兩個問題。一個是把小信號放大到需要的幅度

21、；另一個就是改造信號形狀，通常稱為濾波成形。目的是放大有用的信號，降低造成，提高信噪比，適合于后續(xù)電路的測量。在這個過程中盡可能不損失有用的信息。主放大器輸出的信號就可以給多道脈沖幅度分析器（mca）進行脈沖的幅度分析。多道譜儀是按照脈沖的幅度對輻射脈沖信號進行分類計數(shù)，以及對后續(xù)輻射信息保存、處理以及獲取有用信息的儀器。它的核心部件就是多道脈沖幅度分析器，其基本功能就是按輸入脈沖的幅度分類計數(shù)。多道脈沖幅度分析器的任務是將被測量的脈沖幅度范圍平均分成2n個幅度間隔，然后測量脈沖幅度在每一個“幅度間隔”內的輸入脈沖個數(shù)，最后得到輸入信號的脈沖幅度分布曲線。其測量方法是用計算機技術中的a/d轉換

22、以及數(shù)據(jù)存儲技術。它將能夠分析的脈沖幅度范圍分成多個幅度間隔，這些幅度間隔的個數(shù)就是脈沖幅度分析器的道數(shù)；而幅度間隔的寬度就是道寬。道數(shù)越多，幅度分布分析得越精細，各個道的計數(shù)也相應減少，則需要測量的時間也要加長，當然硬件電路也就越復雜，所以不應盲目地追求道數(shù)。通常，要求在幅度峰的半寬度范圍內應有510道，對于采用nai(tl)探測器的多道能譜儀，由于它的能量分辨率比較差，128道至512道就能滿足測量要求2,3。2.2.2 閃爍探測效率要知道探測效率與射線能量之間的數(shù)值關系，我們首先要探測效率的定義和如何計算探測效率。目前，作為探測效率的定義是多樣的，例如集合效率、本征效率、峰-總比效率、本

23、征全能效率、絕對全能峰效率、絕對雙逃逸峰效率和絕對單逃逸峰效率等3。在本文中，我們是用絕對全能峰效率來進行對探測器探測效率的計算。如果只測量全能峰內的計數(shù)率，則探測裝置的噪聲和干擾輻射對計數(shù)率的影響都可減少。全能峰很容易辨認，從測得的能譜中求得全能峰內的凈計數(shù)率都是比較容易的。根據(jù)全能峰內的凈計數(shù)率與全能峰探測效率，求樣品射線發(fā)射率的測量方法稱為全能峰法。利用全能峰法我們就可以測出全能峰探測效率。例如，在一定條件，測量發(fā)射率為n0的標準源，得到全能峰內的凈計數(shù)率為np則全能峰探測效率為5： （2-1） 3 構建模型3.1 蒙特卡羅方法及mcnp程序簡介蒙特卡羅方法又稱隨機抽樣技巧或統(tǒng)計試驗方法

24、。半個多世紀以來，由于科學技術的發(fā)展和電子計算機的出現(xiàn)與發(fā)展，這種方法作為一種獨立的方法被提出來，并首先在核武器的試驗與研制中得到應用。3.1.1 蒙特卡羅計算方法及其應用范圍蒙特卡羅方法是一種計算方法，但與一般數(shù)值計算方法有很大區(qū)別。它是以概率統(tǒng)計理論為基礎的一種方法。由于蒙特卡羅方法能夠比較逼真地描述事物的特點及物理實驗過程，解決一些數(shù)值方法難以解決的問題，因而該方法的應用領域日趨廣泛4。蒙特卡羅方法的基本思想，當所求問題的解是某個事件的概率，或者是某個隨機變量的數(shù)學期望，或者是與概率、數(shù)學期望有關的量時，通過某種試驗的方法，得出該事件發(fā)生的頻率，或者該隨機變量若干個具體觀察值的算術平均值

25、，通過它得到問題的解4,5。蒙特卡羅方法的優(yōu)點，能夠比較逼真地描述具有隨機性質的事物的特點及物理實驗過程；受幾何條件限制??；收斂速度與問題的維數(shù)無關；具有同時計算多個方案與多個未知量的能力；誤差容易確定；程序結構簡單，易于實現(xiàn)。同時它的缺點，收斂速度慢；誤差具有概率性；在粒子輸運問題中，計算結果與系統(tǒng)大小有關。故使用蒙特卡羅方法需“揚長避短”，發(fā)揮蒙特卡羅的特長，使其應用范圍更加廣泛4,5,6。蒙特卡羅方法應用范圍包括：粒子輸運問題，統(tǒng)計物理，典型數(shù)學問題，真空技術，激光技術以及醫(yī)學，生物，探礦等方面。其在粒子輸運問題中的應用范圍主要包括：實驗核物理，反應堆物理，高能物理等方面。在實驗核物理中

26、的應用范圍主要包括：通量及反應率，中子探測效率，光子探測效率，光子能量沉積譜及響應函數(shù)，氣體正比計數(shù)管反沖質子譜，多次散射與通量衰減修正等方面6。3.1.2 mcnp程序簡介mcnp程序全名為monte carlo neutron and photo transport code。它是美國los alamos國家實驗室開發(fā)的大型多功能通用蒙特卡羅程序，可以計算中子、光子和電子的聯(lián)合輸運問題以及臨界問題，中子能量范圍從10-11mev至20mev，光子和電子的能量范圍從1kev至1000mev。程序采用獨特的曲面組合幾何結構，使用點截面數(shù)據(jù)，程序通用性較強，與其它程序相比，mcnp程序中的減方差

27、技巧是比較多而全的4,5,6。mcnp程序主要特點，程序中的幾何是三維任意組態(tài)；可以使用精細的點截面數(shù)據(jù)；該程序功能齊全；在減小方差技巧方面，內容十分豐富；具有很強的通用性7。mcnp的輸入包括幾個文件，但主要的一個是由用戶編寫的inp文件，該文件包括描述問題所必須的全部輸入信息。文件采用卡片結構，每行代表一張卡片，文件由一系列卡片組成，對于任一特定的問題，只需用到inp全部輸入卡片的一小部分。例如對不同的源，只須修改源的描述卡；不同的計數(shù)類型，有對應的計數(shù)描述卡8。mcnp輸入數(shù)據(jù)中，物理量的單位見表3-1：表3-1 mcnp輸入文件中物理量的單位長度厘米能量mev時間10-8秒溫度mev(

28、kt)原子密度1024個原子/厘米3截面10-24厘米2原子量中子質量的1.008664967質量密度克/厘米3阿伏伽德羅常數(shù)6.0231023mcnp程序輸入文件的基本形式，使用這種形式建立一個蒙特卡羅計算問題，對問題的幾何結構、材料、記數(shù)要求等等給以描述，如果需要，可以直接運行。該文件的格式如下：表3-2 mcnp初始運行的輸入文件信息塊空行分隔符選擇項標題卡僅一行，占用地1-80列。作為輸入標題。柵元卡空行分隔符定義構成整個系統(tǒng)的各個基本戒指單元以及相應的物理信息。曲面卡空行分隔符定義組成柵元的曲面信息。數(shù)據(jù)卡空行分隔符其它數(shù)據(jù)，包括問題類型、源描述、材料描述、計數(shù)描述，問題截斷條件等。

29、其它選擇項信息塊空行分隔符選擇項continue寫在1-8列數(shù)據(jù)卡空行分隔符只允許部分數(shù)據(jù)卡。（fq,dd,nps,ctme,idum,rdum,prdmp,lost,dbcn,print,kcode,mplot,za,zb和zc）其它選擇項卡片填寫格式注意問題，inp輸入文件的每一行（稱之為一張卡片）都限于使用第180列并構成卡片映象。大部分輸入卡片按行填寫；然而，對數(shù)據(jù)卡允許按列填寫。$符號為它所在那行數(shù)據(jù)的結束符，在$符號后面的內容作為注釋，它可從$符號后面的任一列開始。標題卡只占一行，整行都可填入用戶需要的信息，也可以是空行。但要注意在其它地方使用空行是作為結束符或者分隔符。輸入文件中

30、，在標題卡之后及最后的空行結束卡之前的任何地方都可插入注釋卡。注釋卡必須是字母“c”寫在15列中的任意位置，且至少用一個空格隔開后面的注釋內容9。其中行輸入格式，柵元卡、曲面卡和數(shù)據(jù)卡的書寫格式是相同的。必須從15列開始填寫這些卡片相應的名字(或編號)和粒子標識符，后面填寫用空格分隔的數(shù)據(jù)項。如果15列為空，則表示它是前一張卡片的繼續(xù)卡。如果在一行的末尾有一個用空格隔開的符號“&”，則表示下一行是該行的繼續(xù)卡，數(shù)據(jù)可填寫在180列。一個數(shù)據(jù)項必須在一張卡片上寫完，不得跨到下一張卡片上。完全空白的一行則為兩組卡片的分隔符。對任何給定的帶有粒子標識符的類型卡只能有一張。需要整數(shù)的數(shù)據(jù)項必須填寫整數(shù)

31、，其它數(shù)據(jù)可填寫為整數(shù)或浮點數(shù)以及mcnp能讀的數(shù)據(jù)10。列輸入格式對輸入柵元參數(shù)及源描述是非常有用的。對于柵元重要性及體積等參數(shù)，輸入按行排列時其可讀性差，且在增加或刪除一些柵元時容易出錯。用列輸入格式，一個柵元的所有柵元參數(shù)是放在標有該柵元名字的那行上。如果刪掉一個柵元，用戶只需刪除該柵元參數(shù)行，而不需要在每一個柵元參數(shù)卡上尋找該柵元所對應的數(shù)據(jù)項。對于源描述也有類似的情況。用列格式，卡片名字逐個放在一個輸入行上，并且在這些卡片名字下面按列列出數(shù)據(jù)項。后續(xù)各行為各個柵元的數(shù)據(jù)。如果填寫某個柵元名字，則必須填寫全部柵元名字，且柵元的順序可以任意排列；如果沒有指定柵元名字，則按柵元卡描述的順序

32、排列。在一個輸入文件中，允許有多個列數(shù)據(jù)塊11。3.2 建立數(shù)學模型3.2.1 點源模擬本實驗，模擬探測器使用nai（tl）多道能譜儀（圖3-1）替代閃爍探測器。為了得到探測器探測效率和射線能量之間的數(shù)值關系，我們使用點源進行模擬，所選用用的點源包括自然射線和人工射線兩大類，最低能量值為0.352mev，最大的能量值2.62mev(詳見表3-5) 7,8。 圖3-1 nai（tl）多道能譜儀及其探頭部分3.2.2 探測器柵元卡片填寫模擬nai(tl)探測器內部結構示意見圖3-2，探測器分6個柵元，使用15個界面組成。下圖中阿拉伯數(shù)字表示該曲面的曲面號，圓圈數(shù)字表示各柵元號。表3-3為探測器各部

33、分的尺寸大小8。 圖3-2 nai(tl)探測器內部結構示意圖表3-3 探測器各部分尺寸大小單位（cm）內徑外徑厚度高nai-3.8100-7.062mgo3.81003.99500.18507.850sio2-3.9950-0.300fe-3.9950-0.300air4.04504.36460.3196-fe4.36464.44060.07609.705柵元卡填寫：1 1 -3.67 1 -2 -3 imp:p=1 imp:e=12 2 -2.00 (5 -6 -3)#1 imp:p=1 imp:e=13 3 -2.65 3 -6 -4 imp:p=1 imp:e=14 4 -7.93 (

34、7 -8 -9)#1#2#3 imp:p=1 imp:e=15 5 -0.001293 (10 -11 -12)#1#2#3#4 imp:p=1 imp:e=16 6 -7.93 (13 -14 -15)#1#2#3#4#5 imp:p=1 imp:e=1第一列表示柵元號（參考圖 3-2）；第二列表示各柵元對應的材料號；第三列表示該材料的密度；第四列表示該柵元組成的曲面號（參考圖 3-2），imp：p=1 表示該柵元中光子的重要性為 1，同樣 imp：e=1 表示該柵元中電子的重要性為 1。對應各個曲面的幾何位置，第一列為曲面號；第二列中 pz 表示垂直 z 軸的平面，cz 表示軸心在 z 軸

35、的圓柱；第三列表示對應第二列方程數(shù)據(jù)項。1 pz 1.8112 cz 3.813 pz 9.4314 pz 9.7315 pz 1.6266 cz 3.9957 pz 1.5768 cz 4.0459 pz 9.78110 pz 0.07611 cz 4.364612 pz 11.28113 pz 014 cz 4.440615 pz 11.357探測器材料組成見表3-3表3-4 探測器組成材料及其各元素含量-物質密度（g/cm3）成分比探測器nai&tli3.67na:-0.1532i:-0.8462tl:-0.0006mgo2.00mg:0.5o:0.5-sio22.65si:0.3333

36、o:0.6667-fe&c7.93fe：-0.99995c:-0.00005-空氣air0.001293n:-0.7562o:-0.2438-注:成分比部分帶負值表示物質中元素的質量比，正值表示物質中元素的原子成分比材料卡描述：m1 11023 -0.15318 53127 -0.846201 81204 -0.000616m2 12000 0.5 8016 0.5m3 14000 0.33333 8016 0.66667m4 26000 -0.99995 6012 -0.00005m5 8016 -0.24377 7014 -0.75623第一列表示表示材料號；第二列表示對應各材料的物質組成

37、比例；以上材料分別對應 nai(tl)晶體、反射體 mgo、有機玻璃 sio2、不銹鋼罐以及空氣12。3.2.3 模擬模型建立點源模擬系統(tǒng)，見圖 3-3，建立三維坐標系，將探測器置于坐標原點，原點正下方 10cm 處放置點源，現(xiàn)以cs137源為例。以原點為中心建立一個半徑為 150cm 的球體屏蔽系統(tǒng)，以便快速結束粒子追蹤，球體填充空氣，密度大小為 0.001293 g/cm3，球體外為真空8。圖3-3 點源模擬系統(tǒng)示意剖面圖3.2.4 源定義我們選取多種射線能量的進行模擬，本文中我們選取自然射線和人工放射性兩大類射線。各個能量和元素見下表8。表3-5 各個核素與能量自然射線核素能量(mev)

38、pb-2140.352bi-2140.609bi-2141.12k-401.46bi-2141.76ti-2082.62人工射線ba-1310.356ba-1310.389cs-1370.662co-601.173co-601.332na-241.369下文中我們用137cs源進行模擬示范。137cs源為各向同性點源，能量為 662kev，位置在原點正下方 10cm 處。發(fā)射粒子權重為 1，發(fā)射源的粒子類型是光子，同時考慮光子作用。描述如下12：mode p esdef erg=0.662 pos=0 0 -10 par=2 wgt=13.2.5 輸入卡片# nai探測器點源137cs模擬81

39、 1 -3.67 1 -2 -3 imp:p=1 imp:e=1 $nai晶體2 2 -2.00 (5 -6 -3)#1 imp:p=1 imp:e=1 $mgo3 3 -2.65 3 -6 -4 imp:p=1 imp:e=1 $sio24 4 -7.93 (7 -8 -9)#1#2#3 imp:p=1 imp:e=1 $不銹鋼殼5 5 -0.001293 (10 -11 -12)#1#2#3#4 imp:p=1 imp:e=1 $空氣6 6 -7.93 (13 -14 -15)#1#2#3#4#5 imp:p=1 imp:e=1 $不銹鋼殼7 7 -0.001293 -16#1#2#3#4

40、#5#6 imp:p=1 imp:e=1 $空氣8 0 16 imp:p=0 imp:e=0 $真空1 pz 1.811 $曲面定義2 cz 3.813 pz 9.4314 pz 9.7315 pz 1.6266 cz 3.9957 pz 1.5768 cz 4.0459 pz 9.78110 pz 0.07611 cz 4.364612 pz 11.28113 pz 014 cz 4.440615 pz 11.35716 so 150mode p e $源定義sdef erg=0.662 pos=0 0 -10 par=2 wgt=1m1 11023 -0.15318 53127 -0.84

41、6201 81204 -0.000616 $材料定義卡m2 12000 0.5 8016 0.5m3 14000 0.33333 8016 0.66667m4 26000 -0.99995 6012 -0.00005m5 8016 -0.24377 7014 -0.75623m6 26000 -0.99995 6012 -0.00005m7 8016 -0.24377 7014 -0.75623f8:p 1 $記數(shù)卡e8 0 1e-5 0.01 400i 1.0 $記數(shù)能量卡ft8 geb 0 0.0252 $記數(shù)特殊處理卡nps 50000000 $計算粒子數(shù)截斷卡print $打印輸出結果

42、4 計算結果與處理4.1 模擬輸出結果4.1.1 自然射線的模擬輸出結果本次模擬粒子數(shù)為50000000pb-214能量為0.352mev的能譜圖如下：圖4-1214pb0.352mev能譜圖bi-214能量為0.609mev的能譜圖如下:圖4-2215bi0.609mev能譜圖 bi-214能量為1.12mev的能譜圖如下:圖4-3 214bi1.12mev能譜圖k-40能量為1.46mev的能譜圖如下:圖4-4 40k1.46mev能譜圖bi-214能量為1.76mev的能譜圖如下:圖4-5 214bi1.76mev能譜圖ti-208能量為2.26mev的能譜圖如下:圖4-6 208ti2

43、.26mev能譜圖對于能量大于1.02mev的射線與探測介質相互作用就有可能產生電子對事件，在此事件中，正電子湮滅產生兩個0.511mev的光子，一個或兩個均逃離探測介質，探測介質吸收余下的能量，則產生單逃逸峰和雙逃逸峰。在實驗結果中這一現(xiàn)象都表現(xiàn)出來了，特別是在1.5mev以上，譜形上又出現(xiàn)單逃逸峰和雙逃逸峰13，見圖4-5和圖4-6中。圖4-1到4-6是所選的自然射線能譜圖，根據(jù)測量各個能譜上全能峰上的相對計數(shù)率，與對應的能量，列出下表：表4-1 各個元素能量與對應的計數(shù)率元素能量（mev）相對計數(shù)率pb-2140.3520.0132bi-2140.6090.00854bi-2141.12

44、0.00513k-401.460.00409bi-2141.760.00350ti-2082.620.002424.1.2 人工射線的模擬輸出結果以下是所選取的人工射線的能譜圖：ba-131能量為0.356mev的能譜圖如下： 圖4-7 131ba0.356mev能譜圖ba-131能量為0.389mev的能譜圖如下： 圖4-8 131ba0.389mev能譜圖cs-137能量為0.662mev的能譜圖如下： 圖4-9 137cs0.662mev能譜圖co-60能量為1.173mev的能譜圖如下： 圖4-10 60co1.173mev能譜圖co-60能量為1.332mev的能譜圖如下： 圖4-1

45、1 60co1.332mev能譜圖na-24能量為1.369mev的能譜圖如下： 圖4-12 24na1.369mev能譜圖圖4-7到4-12是所選的人工射線能譜圖，根據(jù)測量各個能譜上全能峰上的相對計數(shù)率，與對應的能量，列出下表：表4-2 各個元素能量與對應的計數(shù)率元素能量（mev）相對計數(shù)率ba-1310.3560.0131ba-1310.3890.0123cs-1370.6620.00796co-601.1730.00491co-601.3320.00442na-241.3690004314.2 數(shù)據(jù)處理與結果根據(jù)表4-1與4-2可分別模擬出人工射線與自然射線的射線能量和閃爍探測器的探測效率之間的數(shù)值模擬圖和擬合關系式。自然射線能量值和相對計數(shù)率的模擬圖如下： 圖4-13 自然射線能量值和相對計數(shù)率的模擬圖則根據(jù)圖4-13擬合出指數(shù)衰減函數(shù)下面關系：p=0.02052*exp(-e/0.51684)+0.00263 （4-1）人工射線能量值和相對計數(shù)率的模擬圖如下：圖4-14 人工射線能量值和相對計數(shù)率的模擬