錳浴測量裝置參數(shù)研究

楊竣凱 張健 王平全 張輝 張亞楠

錳浴測量裝置參數(shù)研究

楊竣凱 張健 王平全 張輝 張亞楠

（中國計量科學研究院 北京 100029）

錳浴法是測量中子源強度普遍采用的方法之一。根據(jù)錳浴測量裝置參數(shù)，采用模擬方法開展錳浴測量裝置參數(shù)的研究，結果表明：MnSO4溶液濃度在17.3%～38.6%時，溶液總中子吸收份額和錳的熱中子俘獲份額較高，減少中子逃逸和截面比的不確定度；球形錳池半徑在45～60 cm時，可獲得較高的溶液體積比活度和較低的中子逃逸份額。同時利用蒙特卡羅模擬程序MCNP（Monte Carlo N–Particle Transport Code System）計算時發(fā)現(xiàn)，當粒子數(shù)達到10萬時，修正因子數(shù)值基本穩(wěn)定，由輸入粒子數(shù)帶來的統(tǒng)計漲落不確定度可以忽略。根據(jù)計算結果建立了一套錳浴測量基準裝置，并進行中子源國際比對和12枚中子源測量，測量結果與參考值均在不確定度范圍內(nèi)一致，為以后錳浴測量裝置的建立提供參考。

錳浴法，蒙特卡羅模擬，溶液濃度，半徑，參數(shù)

放射性同位素中子源廣泛應用于核技術應用、核物理實驗和中子計量等領域［1?3］，其中子源強度（也稱為中子發(fā)射率）的準確性直接影響其他儀器設備參數(shù)的測量。目前，國際上普遍采用錳浴法測量放射性同位素中子源的強度，該方法也是國際計量局電離輻射咨詢委員會中子源發(fā)射率K9關鍵比對（Bureau International des Poids et Mesures Consultative Committee for Ionizing Radiation Section III-K9， BIPM CCRI （III）-K9）中最主要方法之一，它具有設備簡單、操作方便和測量不確定度小等優(yōu)點。

不同錳浴測量裝置的參數(shù)，尤其是溶液濃度和錳池半徑，會影響中子逃逸份額和錳的中子吸收份額，從而影響探測器測量結果的統(tǒng)計漲落和準確性；同時可能造成較大的屏蔽負擔，增加裝置以外的建設成本。盡管《計量測試技術手冊》電離輻射部分對于錳池設計給出部分定性要求，如錳池較大、溶液濃度較高［4］，但并未給出具體的參數(shù)范圍。因此，研究錳浴裝置參數(shù)的最佳范圍，對錳浴測量裝置的合理設計并提高測量效率具有重要的參考價值。

1 錳浴法

如圖1所示，被測中子源發(fā)射的中子并非全部被55Mn俘獲，部分中子還會與MnSO4溶液中的其他成分發(fā)生不同反應，同時由于錳池尺寸有限，部分中子會從錳池邊界逃逸［6］，因此需要對計數(shù)率的結果進行修正，主要包括錳俘獲份額、氧硫快中子俘獲份額、源自吸收份額和中子逃逸份額，將計數(shù)率和各修正因子代入式（1）即可得到中子源強度［7］：

圖1　中子進入錳池后所發(fā)生的事件類型

式中：為中子源強度，s–1；為計數(shù)率，s–1；為錳的熱中子俘獲份額占溶液俘獲份額中子的比例；為氧硫快中子俘獲份額；源自吸收份額；為中子逃逸份額。

上述修正因子除了與被測中子源的能譜、包殼和尺寸等相關外，還與錳池尺寸和MnSO4溶液濃度兩項錳浴裝置的參數(shù)有關。

2 分析方法

通過改變MnSO4溶液濃度和錳池尺寸，可以獲得這兩個參數(shù)與各項修正因子的關系，進而獲得這兩個參數(shù)的最佳范圍。早期修正因子主要通過實驗方法確定，每次改變錳池尺寸和濃度都需要重新刻度γ測量系統(tǒng)的探測效率，需要另外制備56Mn放射源，根據(jù)制備方法和放射源管理規(guī)定，制備56Mn放射源比較復雜，降低了可操作性。此外，逃逸中子、氧硫俘獲和源自吸收等修正因子通過實驗測量難度較大，在不同物理模型和錳池參數(shù)條件下利用經(jīng)驗公式計算會引入較大誤差，從而影響測量結果的準確性，因此實驗方法局限性較大。

隨著中子輸運計算技術的發(fā)展，目前可利用蒙特卡羅模擬方法（Monte Carlo N–Particle Transport Code System，MCNP）針對不同的錳浴測量裝置參數(shù)，精確計算各項修正因子，解決實驗方法存在的問題［8?10］。通常，基于模擬方法獲得的修正因子，錳浴法中子源強度測量結果的合成不確定度基本在0.8%（=1）左右，近兩次的中子發(fā)射率國際比對也充分證明模擬方法的可靠性，因此本工作利用MCNP研究了錳浴測量裝置參數(shù)的最佳范圍。

3 參數(shù)研究

以一枚241Am-Be中子源為例，中子源活性區(qū)半徑0.95 cm，高1.8 cm，其中241AmO2質(zhì)量分數(shù)5.67%，9Be質(zhì)量分數(shù)94.33%，密度1.1 g?cm–3，中子源包殼內(nèi)的活性區(qū)外有鋁墊片。測量使用的不銹鋼中子源定位容器內(nèi)半徑1.6 cm，厚度1 mm，其他結構的材料參數(shù)和物理模型見表1和圖2。

被測241Am-Be源包殼外的中子能譜采用ISO 8529-1—2001給出的能譜［11］。根據(jù)實際測量條件，中子源包殼內(nèi)的鋁墊片、包殼及定位容器會影響中子出射能譜，導致活性區(qū)發(fā)射的中子數(shù)與進入錳池的中子數(shù)不一致，因此需要對中子源活性區(qū)能譜和發(fā)射權重進行修正。首先利用MCNP將包殼外的中子能譜多次反迭代至源活性區(qū)內(nèi)，使其發(fā)射至包殼表面的源能譜與ISO 8529近似；由于源材料和源包殼的（n，2n）反應，將源發(fā)射中子權重設為0.984，以對源包殼外中子進行歸一化處理［8］。

表1 基于MCNP的錳浴測量裝置結構及材料參數(shù)

圖2　MCNP模擬的錳池幾何模型示意圖

默認參數(shù)下，錳池半徑為55 cm，MnSO4溶液氫錳原子數(shù)比H/Mn為60，分析過程中，當某一參數(shù)發(fā)生改變時，其他參數(shù)均保持不變。

MnSO4溶液吸收的中子份額abs、各項修正因子子項及標準偏差需要被嚴格記錄。修正因子主要包括：1）錳池邊界泄漏；2）、、的快中子俘獲n和熱中子俘獲n，包括氧原子的（n，γ）和（n，α）反應3、O3，的（n，γ）、（n，α）、（n，p）反應2、2、1和的（n，γ）、4；3）中子源本身和邊界材料對中子的俘獲n，包括中子源材料俘獲1、源包殼俘獲2、源容器俘獲3以及鋁墊片俘獲4；4）Mn的（n，γ）反應俘獲熱中子1［6，12–14］。修正因子為各子項之和，源自吸收和包殼材料修正因子=1+2+3+4，、俘獲的修正因子=1+2+3，錳俘獲熱中子份額占溶液俘獲中子比例=1/（1+2+3+4）。

3.1　輸入粒子數(shù)

圖3　各個記錄量的吸收份額比和標準偏差比

3.2　 MnSO4溶液濃度

圖4　MnSO4溶液、H和Mn的中子吸收份額

表2 各國家計量院錳浴測量裝置NH/NMn[16]

3.3　錳池半徑

MnSO4溶液對中子慢化、吸收能力及裝置的中子泄漏受錳池半徑影響，且溶液比活度隨55Mn的中子吸收份額而改變，從而影響γ測量裝置的計數(shù)率和統(tǒng)計漲落。由于測量時，僅對一定體積的溶液進行測量，因此該處定義體積比活度v，即單位體積放射源的活度；而溶液v與55Mn的中子俘獲份額成正比，可將溶液v轉換成體積比吸收份額r，即單位體積內(nèi)55Mn中子俘獲份額（式（2））。

式中：r為體積比活度份額，cm–3；為錳池半徑，cm。

如圖5所示，假設H/Mn=60，當錳池半徑增大時，MnSO4溶液體積增多，導致溶液的慢化和吸收能力增強，Mn的中子吸收份額4上升，錳池的中子泄漏份額減少，當錳池半徑超過60 cm之后，Mn的中子吸收份額幾乎不再隨半徑的增大而增大，基本接近飽和，但v和r會下降，導致γ測量裝置的計數(shù)率下降和統(tǒng)計漲落升高，因此錳池半徑不宜超過60 cm；當錳池半徑減少時，其中子泄漏份額不斷增大，=60 cm時僅有0.5%，而=15 cm時中子逃逸份額高達49.6%，接近一半的中子泄漏，不僅增大了錳浴測量裝置的屏蔽負擔，泄漏出來的中子可能通過墻面、地面以及其他物體反射再次進入錳池中造成中子的二次慢化和吸收，影響中子源強測量結果，并使結果的不確定度增大。因此，為有效平衡溶液的比活度和中子泄漏份額的關系，保證測量結果準確性，錳池半徑宜在45～60 cm較為合適，此時溶液體積比活度和中子泄漏份額分別為9.7×10–7～4.3×10–7cm–3和2.3×10–2～5×10–3。如表3所示，多數(shù)國家實驗室的錳池半徑基本處于此范圍內(nèi)。

圖5　錳浴球半徑對溶液體積比活度、中子泄漏份額和Mn吸收份額的影響

表3 NIM測量結果不確定度清單

4 實驗驗證

4.1　錳浴測量裝置

根據(jù)上述結論，中國計量科學研究院建立的錳浴測量裝置（中子源強度基準）（圖6）主要由球形錳池、中子源定位容器、MnSO4溶液循環(huán)管路、鉛屏蔽室、γ測量裝置（雙路NaI探測器）及相關電子學設備組成，其中錳池半徑為55 cm，溶液濃度為21.9%，H/Mn為59.704。

圖6　錳浴測量裝置 (a) 球形錳池，(b) 鉛室中的NaI探測器Fig.6　Mn bath measurement device (a) Spherical Mn bath, (b) NaI detector in lead chamber

4.2　國際比對

中子源強度比對是國際計量局電離輻射咨詢委員會組織的中子計量關鍵比對之一，最近完成的國際比對中，主導實驗室為英國國家物理實驗室（National Physical Laboratory，NPL），參比實驗室為中國計量科學研究院（National Institute of Metrology，NIM），比對結果分別是NPL=（8.082±0.102）×106s–1，NIM=（8.167+0.127）×106s–1，結果相差1.05%，中國計量科學研究院不確定分析見表3，NPL測量結果不確定度NPL=0.63%（=1），將結果和不確定度代入式（3），得到歸一化誤差n=0.135，表明依據(jù)上述研究結論建立的錳浴測量裝置測量結果準確，實現(xiàn)了國際等效。

式中：為中子源強；為不確定度；n為歸一化誤差，在0～1范圍內(nèi)表示參與實驗室的測量值之間的差異小于或等于兩個實驗室的合成擴展不確定度（=2）。

4.3　 12枚中子源測量

中國計量科學研究院新建的熱中子參考輻射裝置（熱中子注量率基準）內(nèi)的12枚241Am-Be中子源與S026/19中子源為同一批次生產(chǎn)的中子源，S026/19由英國NPL實驗室測量，結果為（1.164±0.015）×107s–1，12枚中子源由計量院測量，結果如表4所示，在不確定度范圍內(nèi)與英國NPL的測量結果一致［17］。

5 結語

表4 12枚241Am-Be中子源測量結果

注：參考時間為2019.05.22

Note: Reference time was 22/05/2019

致謝 感謝中國原子能科學研究院李春娟老師對MCNP模擬提供的指導和幫助。

作者貢獻聲明 楊竣凱、王平全：直接參與、文章撰寫；張健、張輝：直接參與、文章撰寫、工作支持；張亞楠：直接參與、工作支持。

1 丁大釗, 葉春堂, 趙志祥. 中子物理學: 原理、方法與應用[M]. 北京: 原子能出版社, 2001.

DING Dazhao, YE Chuntang, ZHAO Zhixiang. Neutron physics: principles, methods and applications[M]. Beijing: Atomic Energy Press, 2001.

2 張鋒, 袁超, 黃隆基. 中國同位素中子源測井技術與應用進展[J]. 同位素, 2011, 24(1): 12–19.

ZHANG Feng, YUAN Chao, HUANG Longji. Technology and application development of isotopic neutron source well logging in China[J]. Journal of Isotopes, 2011, 24(1): 12–19.

3 劉鎮(zhèn)洲, 陳金象. Am-Be中子源輻射場周圍劑量當量與吸收劑量的計算[J]. 原子能科學技術, 2009, 43(1): 11–15.

LIU Zhenzhou, CHEN Jinxiang. Calculations of dose and ambient dose equivalent in radiation field for Am-Be neutron source[J]. Atomic Energy Science and Technology, 2009, 43(1): 11–15.

4 《計量測試技術手冊》委員會. 計量測試技術手冊第12卷: 電離輻射計量測試技術手冊[M]. 北京: 中國計量出版社, 1997.

Technical Manual of Measurement and Testing Committee. Measurement and testing technical manual Volume 12: technical manual for measurement and testing of ionizing radiation[M]. Beijing: China Metrology Press, 1997.

5 盧希庭. 原子核物理[M]. 第2版. 北京: 原子能出版社, 2000.

LU Xiting. Nuclear physics[M]. 2nd Ed. Beijing: Atomic Press, 2000.

6 李春娟, 劉毅娜, 張偉華, 等. 錳浴法測量中子源發(fā)射率中各修正因子的蒙特卡羅計算[J]. 原子能科學技術, 2014, 48(10): 1876–1881. DOI: 10.7538/yzk.2014.48.10.1876.

LI Chunjuan, LIU Yina, ZHANG Weihua,. Monte Carlo calculation of correction factors for radionuclide neutron source emission rate measurement by manganese bath method[J]. Atomic Energy Science and Technology, 2014, 48(10): 1876–1881. DOI: 10.7538/yzk.2014.48.10.1876.

7 Croft S, Favalli A. How the choice of data reduction can strongly influence uncertainty assessment: a re-analysis of Mn-bath experiments[J]. Radiation Measurements, 2012, 47(7): 481–485. DOI: 10.1016/j.radmeas.2012.05.008.

8 李春娟, 陳軍, 吳海成, 等.238U裂變電離室修正因子及探測響應的蒙特卡羅模擬計算[J]. 原子能科學技術, 2009, 43(3): 208–214.

LI Chunjuan, CHEN Jun, WU Haicheng,. Monte-Carlo simulation for correction factors and detection response of238U fission chamber[J]. Atomic Energy Science and Technology, 2009, 43(3): 208–214.

9 Khabaz R. Appraisement of the correction factors for neutron reaction in the manganese bath using Monte Carlo calculation[J]. Journal of Radioanalytical and Nuclear Chemistry, 2012, 293(2): 455–462. DOI: 10.1007/s10967-012-1792-0.

10 許淑艷. 蒙特卡羅方法在實驗核物理中的應用[M]. 第2版. 北京: 原子能出版社, 2006.

XU Shuyan. Application of Monte Carlo method in experimental nuclear physics[M]. 2nd Ed. Beijing: Atomic Press, 2006.

11 ISO 8529-1—2001: Reference neutron radiations-Part 1:characteristics and method of production[S]. 2001.

12 Murphey W M. The correction factor for fast neutron reactions on sulfur and oxygen in the manganous-sulfate-bath calibration of neutron sources[J]. Nuclear Instruments and Methods, 1965, 37: 13–21. DOI: 10.1016/0029-554X(65)90330-7.

13 凌玉德. 錳浴法測中子源強度的修正項及241Am(α,n)Be中子源能譜的計算[J]. 原子能科學技術, 1984, 18(6): 661–666.

LING Yude. Correction term of neutron source intensity measured by manganese bath method and calculation of energy spectrum of241Am(α,n)Be neutron source[J]. Atomic Energy Science and Technology, 1984, 18(6): 661–666.

14 劉毅娜, 王志強, 李春娟, 等. 錳浴法絕對測量中子源發(fā)射率[J]. 原子能科學技術, 2013, 47(6): 1044–1047. DOI: 10.7538/yzk.2013.47.06.1044.

LIU Yina, WANG Zhiqiang, LI Chunjuan,. Absolute measurement of neutron source emission rate with manganese bath method[J]. Atomic Energy Science and Technology, 2013, 47(6): 1044–1047. DOI: 10.7538/yzk.2013.47.06.1044.

15 倪育才. 實用測量不確定度評定[M]. 第2版. 北京: 中國計量出版社, 2007.

NI Yucai. Evaluation of practical measurement uncertainty[M]. 2nd Ed. Beijing: China Metrology Publishing House, 2007.

16 Roberts N J, Jones L N, Wang Z,. International key comparison of measurements of neutron source emission rate (1999–2005): CCRI(III)-K9.AmBe[J]. Metrologia, 2011, 48(1A): 06018. DOI: 10.1088/0026-1394/48/1a/06018.

17 楊竣凱, 王平全, 張輝, 等. 熱中子參考輻射裝置參數(shù)的實驗測量[J]. 核技術, 2021, 44(11): 110501. DOI: 10.11889/j.0253-3219.2021.hjs.44.110501.

YANG Junkai, WANG Pingquan, ZHANG Hui,. Experimental measurement of parameters of thermal neutron reference field[J]. Nuclear Techniques, 2021, 44(11): 110501. DOI: 10.11889/j.0253-3219.2021.hjs.44.110501.

Parameters of manganese bath measurement device

YANG Junkai ZHANG Jian WANG Pingquan ZHANG Hui ZHANG Yanan

()

Mn bath method is one of the most commonly used methods to measure the intensity of radionuclide neutron sources in the world.This study aims to guarantee the accuracy of numerical transmission and give a reasonable parameter range during the design of manganese bath device. [Methods] First of all, Monte Carlo N–Particle Transport Code System (MCNP) was employed to calculate the correction factors, solution concentrations and Mn-bath sphere. Then, a reference device of Mn bath measurement was setup according to the MCNP calculation results, and applied to the measurement of 12 neutron sources located in the standard of thermal neutron fluence rate with International Comparison of neutron sources.Simulation results show higher ratio of total neutron absorption and manganese neutron capture is obtained when the solution concentration is 17.3%～38.6% in which the uncertainty of neutron-escaping and section ratio is reduced, and the radius of Mn-bath sphere is 45～60 cm in which the volume specific activity of solution is increased and neutron-escaping is decreased. When the number of input particles is up to 100 000, the values of correction factors remain stable, hence the statistical fluctuation can be ignored. The international comparison of neutron sources for the measurement of 12 neutron sources indicates consistent results with the reference values within the uncertainty range.The parameters range of this study provides an effective reference for the establishment of Mn bath device in the future.

Mn bath method, Monte Carlo simulation, Solution concentration, Radius, Parameters

Supported by Fundamental Research Funds of National Institute of Metrology, China (No.AKYZZ2113), National Key Research and Development Program of China (No.2017YFF0206205)

YANG Junkai, male, born in 1994, graduated from China Institute of Atomic Energy with a master's degree in 2020, focusing on particle physics and nuclear physics

ZHANG Hui, E-mail: zhanghui@nim.ac.cn

2021-12-29,

2022-03-03

TL99

10.11889/j.0253-3219.2022.hjs.45.060501

中國計量科學研究院基本業(yè)務費(No.AKYZZ2113)、國家重點研發(fā)計劃(No.2017YFF0206205)資助

楊竣凱，男，1994年出生，2020年于中國原子能科學研究院獲碩士學位，研究領域為粒子物理與原子核物理

張輝，E-mail：zhanghui@nim.ac.cn

2021-12-29，

2022-03-03