為了增進在清華開放式水池反應器(THOR)硼中子捕獲治療(BNCT)設施的中子能譜測量精準度，本研究將傳統的中子活化偵檢器測量法做進一步改良，將傳統的裸活化箔片嵌在球型緩速體及吸收體中，有三項優點：(1)可減少中子射源的角分布對偵檢器的影響；(2)利用緩速體及吸收體來改變球心的箔片對中子的響應，可產生多種能涵蓋不同能量範圍的響應函數；(3)使用球型偵檢器可以改善箔片對超熱中子的偵測效率，減少整體實驗量測的時間，有利於精進能譜測量的結果。
本研究會介紹這種新型的球型中子活化偵檢器，偵檢器的響應函數計算以及利用SAND-Ⅱ能譜演算法為主的能譜反解程式SAND-EX。對於球型活化檢器的測量效果，會先以理論測試來評估能譜反解的結果，並實際在THOR BNCT設施進行實驗量測，傳統的裸活化箔片會做為反解結果比較的基準，不論理論或實驗，結果都證明使用球型偵檢器能有效改善中子能譜測量結果。
本研究會探討反解參數不同對反解結果的影響，包含使用不同初始能譜對反解結果影響，射源的角分布不同對偵檢器響應函數的影響，還有在活化箔片中加入包覆鎘片的活化箔片對反解結果的影響，另外本研究也利用MAXED反解程式以探討使用不同的反解演算法的結果。總結上述反解參數的靈敏度影響的結果，在BNCT中子能譜測量上，球型偵檢器還是比傳統的活化箔片能得到較穩定較好的結果。

A set of spherical-type activation detectors was developed aiming to provide better determination of the neutron spectrum at the Tsing Hua Open-pool Reactor (THOR) BNCT facility. The activation foil embedded in a spherically shaped holder exhibits three advantages: (1) minimizing the effect of neutron angular dependence, (2) creating response functions with broadened coverage of neutron energies by introducing additional moderators or absorbers to the central activation foil, and (3) reducing overall detection time by increasing efficiencies of most activation foils under the irradiation of epithermal neutrons. This paper presents the design concept and the calculated response functions of new detectors. Theoretical and experimental demonstrations of the performance of the detectors are provided through comparisons of the unfolded neutron spectra determined using this method and conventional multiple-foil activation techniques.

摘要 i
Abstract ii
致謝 iii
目錄 iv
表目錄 vi
圖目錄 viii
第一章、緒論 1
1.1 硼中子捕獲治療 1
1.2 中子能譜量測 2
1.2.1 中子束特性測量 2
1.2.2 球型活化偵檢器 3
第二章、中子活化分析 4
2.1 裸活化箔片(bare activation foil) 4
2.1.1 活化箔片介紹 4
2.1.2 文獻回顧 8
2.2 球型中子活化偵檢器 11
2.3 THOR 實驗設施 12
2.3.1 硼中子捕獲治療設施 12
2.3.2 中子即時監測系統 14
2.3.3 高純度鍺偵檢器 14
2.3.4 能譜分析軟體Gamma Vision-32 19
2.4 飽和活度 21
第三章、中子能譜反解 24
3.1 SAND-Ⅱ能譜反解演算 24
3.2 能譜反解程式 26
3.2.1 SAND-EX能譜反解程式 26
3.2.2能譜平滑修正 29
3.3 初始能譜 30
3.4偵檢器響應函數 31
3.4.1蒙地卡羅法計算響應函數 31
3.4.2響應函數計算結果 32
第四章、理論測試與實驗驗證 38
4.1 偵檢器選擇 38
4.2 結果與討論 41
4.2.1 理論測試 41
4.2.2 實驗量測 45
第五章、靈敏度分析 51
5.1 初始能譜 51
5.2 射源角分布對響應函數影響 55
5.3 活化箔片與鎘片 67
5.3.1 響應函數 67
5.3.2 結果與討論 70
5.4 MAXED反解程式 76
5.4.1 MAXED反解演算 76
5.4.2 初始能譜 78
5.4.3 結果與討論 79
5.4.4 SAND-EX與MAXE程式結果比較 83
第六章、結論與建議 89
參考文獻 91

1. IAEA, Current State of Neutron Capture Therapy. 2001, IAEA: Vienna.
2. Voorbraak, W.P., Recommendations for the Dosimetry of Boron Neutron Therapy. 2003, NRG: Petten.
3. 張志豪, 清華水池式反應器硼中子捕獲治療射束中子通率測量, 工程與系統科學系. 2006, 國立清華大學.
4. 蔡碧恩, 利用間接中子造影法進行硼中子捕獲治療中子射束之特性研究, 工程與系統科學系. 2009, 國立清華大學.
5. 劉淵豪, 硼中子捕獲治療超熱中子束之中子特性研究, 工程與系統科學系. 2009, 國立清華大學.
6. Marek, M. and L. Viererbl, Bonner sphere spectrometer for characterization of BNCT beam. Applied Radiation and Isotopes, 2011. 69(12): p. 1918-1920.
7. Ueda, H., et al., The optimization study of Bonner sphere in the epi-thermal neutron radiation field for BNCT. Applied Radiation and Isotopes, 2011. 69: p. 1657-1659.
8. 劉鴻鳴, 先進硼中子捕獲治療劑量系統研發. 2013.
9. Liu, Y.-W.H., et al., Renovation of epithermal neutron beam for BNCT at THOR. Applied Radiation and Isotopes, 2004. 61(5): p. 1039-1043.
10. Liu, Y.H., et al., Neutron Spectra Measurement and Comparison of HFR and THOR BNCT Beams. Applied Radiation and Isotopes, 2009. 67(7-8 Suppl): p. S137-S140.
11. Liu, Y.H., et al., On-line Neutron Monitor System of Epithermal Neutron Beam for BNCT at THOR, in International Congress on Neutron Capture Therapy. 2006, Takamatsu kagawa: Japan. p. 413-416.
12. ORTEC, ORTEC MCB Connections-32 Hardware Property Diologs Manual. 2004, USA.
13. Berg, S. and W.N. McElroy, A computer-automated iterative method for neutron flux spectra determination by foil activation, in SAND-Ⅱ and associate code. 1967, Air Force Weapons Laboratory: USA.
14. 陳韋霖, SAND-EXII使用說明. 2012.
15. Liu, Y.H., et al., Coarse-scailing adjustment of fine-group neutron spectra for epithermal neutron beams in BNCT using multiple activation detectors. Nuclear Instruments and Methods in Physics Research A, 2009. 598: p. 764-773.
16. Liu, Y.H., W.L. Chen, and S.H. Jiang, Improvement of Activation -Detector-Based Spectrum Unfolding Technique for BNCT. IEEE Transactions on Nuclear Science, 2013. 60(2): p. 783-791.
17. OriginLab. Origin: Smoothing using Adjacent Averaging. Available from: http://www.originlab.com/index.aspx?go=Products/Origin/DataAnalysis/SignalProcessing/SmoothingAndFitting&pid=62.
18. Pelowitz, D.B.E., MCNPX user's manual, v2.7.0. 2009, Los Alamos National Laboratory: USA.
19. Reginatto, M., P. Goldhagen, and S. Neumann, Spectrum unfolding, sensitivity analysis and propagation of uncertainty with the maximun entropy deconvalution code MAXED. Nuclear Instruments and Methods in Physics Research A, 2002. 476: p. 242-246.
20. Reginatto, M. and P. Goldhagen, MAXED, A COMPUTER CODE FOR MAXIMUM ENTROPY DECONVOLUTION OF MULTISPHERE NEUTRON SPECTROMETER DATA. Health Physics, 1999. 77(5): p. 579-583.
21. UMG 3.3 - Unfolding with MAXED and GRAVEL. 2004, OAK RIDGE NATIONAL LABORATORY.