基于252Cf源驅(qū)動(dòng)噪聲分析法的鈾部件裂變中子時(shí)間關(guān)聯(lián)符合計(jì)數(shù)模擬研究

周 密 魏 彪 楊 帆 馮 鵬 譙 梁

1（重慶大學(xué)光電技術(shù)及系統(tǒng)教育部重點(diǎn)實(shí)驗(yàn)室 重慶400044）2（重慶理工大學(xué)光電信息學(xué)院 重慶 400054）

基于252Cf源驅(qū)動(dòng)噪聲分析法的鈾部件裂變中子時(shí)間關(guān)聯(lián)符合計(jì)數(shù)模擬研究

周 密1,2魏 彪1楊 帆1馮 鵬1譙 梁1

1（重慶大學(xué)光電技術(shù)及系統(tǒng)教育部重點(diǎn)實(shí)驗(yàn)室 重慶400044）2（重慶理工大學(xué)光電信息學(xué)院 重慶 400054）

依據(jù)中子源驅(qū)動(dòng)噪聲分析法原理，針對(duì)鈾材料或鈾部件的質(zhì)量屬性測(cè)量問題，采用蒙特卡羅(Monte Carlo)模擬方法，模擬研究了幾何形狀相同、富集度相同而質(zhì)量不同的球形金屬鈾部件的質(zhì)量測(cè)量，獲得了與鈾部件質(zhì)量測(cè)量相對(duì)應(yīng)的時(shí)間關(guān)聯(lián)符合計(jì)數(shù)分布。研究結(jié)果表明，中子源注入鈾部件后所構(gòu)成的核信息系統(tǒng)，可根據(jù)源與探測(cè)器間的中子時(shí)間關(guān)聯(lián)符合計(jì)數(shù)分布獲得鈾部件的質(zhì)量等特性。通過對(duì)半徑4.67 cm、質(zhì)量8.00 kg的球形鈾部件樣品的模擬研究與分析，獲得源-探測(cè)器之間的時(shí)間關(guān)聯(lián)符合計(jì)數(shù)，基此，模擬研究了幾何形狀相同、富集度相同而質(zhì)量不同的球形鈾部件，實(shí)現(xiàn)了對(duì)鈾部件質(zhì)量屬性的測(cè)量。

源驅(qū)動(dòng)噪聲分析法，鈾部件，中子信號(hào)，Monte Carlo模擬，時(shí)間關(guān)聯(lián)符合計(jì)數(shù)

鈾材料或鈾部件是核武器最重要的組成之一，其質(zhì)量大小及形狀等基本屬性，不僅直接攸關(guān)核武器的性能，而且是核軍控核查認(rèn)證的基礎(chǔ)[1–4]。鈾材料質(zhì)量屬性的測(cè)量包括實(shí)驗(yàn)測(cè)量研究和計(jì)算機(jī)蒙特卡羅模擬(Monte Carlo，MC)研究。美國(guó)能源部橡樹嶺國(guó)家實(shí)驗(yàn)室(Oak Ridge National Laboratory) Mihalczo等[5–10]采用中子源驅(qū)動(dòng)時(shí)間關(guān)聯(lián)符合測(cè)量法進(jìn)行鈾材料或部件質(zhì)量屬性的測(cè)量，獲得了飛速發(fā)展。此方法的驅(qū)動(dòng)外源需要選擇定時(shí)中子源，如252Cf自發(fā)裂變中子源[11–13]或定時(shí)D-T中子源[14]，基于此，還研制成了核材料識(shí)別系統(tǒng)(Nuclear material identification system, NMIS)，且已成為國(guó)際核軍控核查領(lǐng)域的研究熱點(diǎn)之一。

然而，由于核軍控核查的特殊性，國(guó)外該領(lǐng)域有關(guān)研究工作采取了技術(shù)封鎖的策略，更由于核材料使用的敏感性及強(qiáng)放射性輻射的危害性等因素，致使對(duì)核材料識(shí)別的實(shí)際研究較為困難。為此，本文對(duì)252Cf源驅(qū)動(dòng)時(shí)間關(guān)聯(lián)符合法測(cè)量鈾材料質(zhì)量屬性進(jìn)行了數(shù)值模擬研究。模擬研究了252Cf源注入被測(cè)的鈾部件后測(cè)量系統(tǒng)內(nèi)的中子輸運(yùn)模擬過程、探測(cè)器對(duì)粒子的響應(yīng)過程以及探測(cè)器間中子時(shí)間關(guān)聯(lián)符合響應(yīng)計(jì)數(shù)分布等，研究了幾何形狀相同、富集度相同而質(zhì)量不同的球形鈾部件的中子時(shí)間關(guān)聯(lián)信號(hào)，以此達(dá)到測(cè)量鈾材料質(zhì)量屬性的目的。

1 原理及方法

鈾部件質(zhì)量屬性的測(cè)量是依據(jù)252Cf源驅(qū)動(dòng)噪聲分析法獲得測(cè)量鈾部件的標(biāo)簽參數(shù)，借以進(jìn)行幾何形狀相同、富集度相同而質(zhì)量不同的球形鈾部件的分析與識(shí)別，其原理如圖1所示?？梢姡瑴y(cè)量系統(tǒng)由兩路組成。第一路是252Cf源電離室，為測(cè)量系統(tǒng)輸入驅(qū)動(dòng)中子噪聲源，并作為一個(gè)探測(cè)器，提供與252Cf電離室時(shí)間相關(guān)的脈沖時(shí)間序列信號(hào)。第二路探測(cè)次臨界系統(tǒng)的泄漏中子，作為測(cè)量系統(tǒng)的第二路輸出信號(hào)。此外，該測(cè)量系統(tǒng)可擴(kuò)展為單路輸入、三路或更多路輸出隨機(jī)信號(hào)的結(jié)構(gòu)，可實(shí)現(xiàn)多路信號(hào)的同步采集。第二路探測(cè)器信號(hào)的采集，由第一路252Cf源信號(hào)觸發(fā)或由記錄系統(tǒng)內(nèi)部控制，統(tǒng)計(jì)結(jié)果與觸發(fā)方式無(wú)關(guān)。源信號(hào)與探測(cè)器信號(hào)相互關(guān)聯(lián)，最后根據(jù)采集到的數(shù)據(jù)分析，獲得源與探測(cè)器間的中子時(shí)間關(guān)聯(lián)符合計(jì)數(shù)分布12()Cτ。它反映了源探測(cè)信號(hào)1()Xt與探測(cè)器探測(cè)信號(hào)2()Xt之間的時(shí)間關(guān)聯(lián)符合特性，其定義為：

其中，T為觀測(cè)時(shí)間，τ為延遲時(shí)間。

為定量研究鈾部件的質(zhì)量屬性，定義時(shí)間關(guān)聯(lián)符合積分計(jì)數(shù)為：

其中，τ1為積分時(shí)間起點(diǎn)，τ2為積分時(shí)間終點(diǎn)。選擇合適的積分時(shí)間區(qū)間，R12(τ1,τ2)能夠比較靈敏地反映出幾何形狀相同、富集度相同而質(zhì)量不同鈾部件的質(zhì)量屬性。

圖1 252Cf源驅(qū)動(dòng)噪聲分析法測(cè)量原理圖Fig.1 Principle of 252Cf source-driven noise analysis measurement.

圖2 質(zhì)量8.00 kg鈾部件的中子時(shí)間關(guān)聯(lián)符合計(jì)數(shù)模擬計(jì)算結(jié)果Fig.2 Simulation results of time-correlated pairs neutron of uranium component of 8.00 kg.

2 計(jì)算結(jié)果與討論

模擬研究計(jì)算的鈾部件富集度為93.15%，它由235U、238U、234U及236U四種核素組成，對(duì)應(yīng)的質(zhì)量分別為93.15%、5.64%、0.97%及0.24%。模擬計(jì)算時(shí)，假設(shè)探測(cè)器探測(cè)能量大于0.1 MeV的所有中子信息，探測(cè)效率為100%。在模擬過程中，為了提高計(jì)算效率，不考慮地面、墻壁等周圍環(huán)境對(duì)測(cè)量結(jié)果的影響。圖2給出了半徑4.67 cm、質(zhì)量8.00 kg的球形鈾部件全能區(qū)(0–20 MeV)的12()Cτ隨時(shí)間分布的曲線(見圖2a)和全能區(qū)中裂變中子與散射中子的12()Cτ隨時(shí)間分布的曲線(見圖2b)。

由圖2 (b)可見，全能區(qū)中包括裂變中子和散射中子，且散射中子的比例或份額還很大，這不利于對(duì)鈾部件質(zhì)量的分辨。為此，必須去除散射中子以提高質(zhì)量分辨的效果。由于裂變中子與散射中子的能量不盡相同，所以其產(chǎn)生時(shí)間域也有所不同[15]。

針對(duì)半徑4.67 cm、質(zhì)量8.00 kg的球形鈾部件，計(jì)算出了不同能區(qū)裂變中子信號(hào)與散射中子信號(hào)的能譜對(duì)比圖，如圖3所示。

由圖3可見，在0–20 MeV能區(qū)中，散射中子的份額比較大。當(dāng)能量下閾值增大，散射中子的份額逐漸減少，說明低能量的散射中子占的比重大。為了盡可能地減小散射中子的影響，模擬研究計(jì)算時(shí)，選擇在1.5–20 MeV能區(qū)的裂變中子信號(hào)進(jìn)行鈾部件的質(zhì)量屬性的能譜分析。基此，根據(jù)中子時(shí)間關(guān)聯(lián)符合法對(duì)三個(gè)球形鈾部件開展模擬研究。

三個(gè)球形鈾部件分別為(1) 半徑5.63 cm、質(zhì)量14.02 kg；(2) 半徑5.884 cm、質(zhì)量16.01 kg；(3) 半徑6.13 cm、質(zhì)量18.10 kg。通過模擬計(jì)算，得到了形狀相同、質(zhì)量不同的三個(gè)鈾部件在1.5–20 MeV能區(qū)裂變信號(hào)的12()Cτ隨時(shí)間分布的曲線，如圖4所示。可見，由12()Cτ隨時(shí)間分布的曲線，可以實(shí)現(xiàn)對(duì)鈾部件質(zhì)量的分辨。

根據(jù)時(shí)間關(guān)聯(lián)符合積分計(jì)數(shù)R12(τ1,τ2)的定義，選擇合適的積分時(shí)間區(qū)間，R12(τ1,τ2)能夠比較靈敏地反映出幾何形狀相同、富集度相同而質(zhì)量不同鈾部件的質(zhì)量屬性。分別計(jì)算裂變中子與散射中子的時(shí)間關(guān)聯(lián)符合積分計(jì)數(shù)，定義鈾部件積分計(jì)數(shù)的比值作為選擇合適的積分時(shí)間區(qū)間的依據(jù)，用R表示為：

其中，R12,fission為裂變中子的時(shí)間關(guān)聯(lián)符合積分計(jì)數(shù)，R12,scatter為散射中子的時(shí)間關(guān)聯(lián)符合積分計(jì)數(shù)。

源與探測(cè)器間中子的時(shí)間關(guān)聯(lián)符合計(jì)數(shù)分布，有的是通過裂變鏈相關(guān)聯(lián)的，而散射中子實(shí)際上毫無(wú)關(guān)系。通過對(duì)積分計(jì)數(shù)比值R的計(jì)算可以為R12(τ1,τ2)選擇合適的積分時(shí)間區(qū)間。對(duì)樣本半徑4.67 cm、質(zhì)量8.00 kg的球形鈾部件，在不同的能區(qū)，考慮了不同的積分時(shí)間區(qū)間，得到鈾部件積分計(jì)數(shù)的比值如表1所示。由表1可見，在能量大于1.0 MeV、時(shí)間大于20 ns后，裂變中子明顯大于散射中子。由此，積分時(shí)間起點(diǎn)1τ取為20 ns，積分時(shí)間終點(diǎn)2τ取為100 ns，可以得到R12(20,100)的模擬統(tǒng)計(jì)結(jié)果，如圖5所示。

圖3 裂變中子和散射中子的時(shí)間關(guān)聯(lián)符合計(jì)數(shù)模擬計(jì)算結(jié)果Fig.3 Simulation results of the time-correlated pairs neutron from fission and scattering.

圖4 不同質(zhì)量鈾部件裂變中子的時(shí)間關(guān)聯(lián)符合計(jì)數(shù)的模擬計(jì)算結(jié)果Fig.4 Simulation results of the time-correlated pairs fission neutrons from uranium components of different masses.

表1 裂變中子和散射中子的時(shí)間關(guān)聯(lián)符合積分計(jì)數(shù)比值模擬計(jì)算結(jié)果Table1 Simulation ratios of time-correlated integral pairsneutrons from fission and scattered neutrons.

圖5 不同質(zhì)量鈾部件裂變中子的時(shí)間關(guān)聯(lián)符合積分計(jì)數(shù)模擬計(jì)算結(jié)果Fig.5 Simulation results of the time-correlated integral pairs fission neutrons from uranium components of different masses.

3 結(jié)語(yǔ)

依據(jù)中子源驅(qū)動(dòng)噪聲分析法原理，針對(duì)鈾材料或鈾部件的質(zhì)量屬性測(cè)量問題，采用Monte Carlo方法，模擬研究了幾何形狀相同、富集度相同而質(zhì)量不同的球形金屬鈾部件的質(zhì)量測(cè)量，獲得了與鈾部件質(zhì)量測(cè)量相對(duì)應(yīng)的時(shí)間關(guān)聯(lián)符合計(jì)數(shù)分布。研究結(jié)果表明，中子源注入鈾部件后由此構(gòu)成的核信息系統(tǒng)，根據(jù)源與探測(cè)器間的中子時(shí)間關(guān)聯(lián)符合計(jì)數(shù)分布，可以分析與識(shí)別鈾部件的質(zhì)量等特性。通過對(duì)樣本半徑4.67 cm、質(zhì)量8.00 kg的球形鈾部件的模擬研究與分析，獲得了源-探測(cè)器之間的時(shí)間關(guān)聯(lián)符合計(jì)數(shù)，基此，模擬研究了幾何形狀相同、富集度相同而質(zhì)量不同的球形鈾部件，實(shí)現(xiàn)了對(duì)鈾部件質(zhì)量屬性的測(cè)量。

1 杜祥琓. 核軍備控制的科學(xué)技術(shù)基礎(chǔ)[M]. 北京: 國(guó)防工業(yè)出版社, 1996

DU Xiangwan. An elementary of nuclear arms control technology[M]. Beijing: National Defense Industry Press, 1996

2 劉成安, 伍鈞. 核軍備控制核查技術(shù)概論[M]. 北京:國(guó)防工業(yè)出版社, 2007

LIU Chengan, WU Jun. An introduction of verification technology of nuclear arms control[M]. Beijing: National Defense Industry Press, 2007

3 Reilly D, Ensslin N, Smith H J. Passive nondestructive assay of nuclear of nuclear materials[C]. Nuclear Regulatory Commission, Washington DC, 1991

4 Mihalczo J T, Pare V K. Nuclear weapons identification system[C]. Presented at the proceedings of the 1996 IEEE nuclear science symposium, Anaheim, CA, USA, 1996

5 Mihalczo J T, King W T, Blakeman E D. Subcriticality measurements for coupled uranium metal cylinders using the252Cf-source-driven neutron noise analysis method[J]. Nuclear Science and Engineering, 1987, 95(1): 1–13

6 Chiles M M, Mihalczo J T, Fowler C E. Small, annular, double-contained252Cf fission chamber for source-driven subcriticality measurements[J]. Institute of Electrical and Electronics Engineers Nuclear Science Symposium, Nuclear Science and Engineering, 1993, 40(4): 816–818

7 Mihalczo J T, Valentine T E. Calculational verification and process control applications utilizing the high sensitivity of noise measurement parameters to fissile system configuration[J]. Nuclear Science & Engineering, 1995, 121(2): 286–300

8 Mihalczo J T, Mullens J A, Mattingly J K, et al. Physical description of nuclear materials identification system (NMIS) signatures[J]. Nuclear Instruments and Methods in Physics Research A, 2000, 450(2): 531–555

9 Valentine T E, Mihalczo J T. MCNP-DSP: A neutron and gamma ray Monte Carlo calculation of source-driven noise-measured parameters[J]. Annals of Nuclear Energy, 1996, 23(3): 12–71

10 Marseguerra M, Padovani E, Pozzi S A. Use of the MCNP-POLIMI code for time-correlation safeguards measurements[J]. Progress in Nuclear Energy, 2003, 43(4): 305–311

11 Mihalczo J T, King W T, Blakeman E D.252Cf source-driven neutron noise analysis[C]. Workshop on subciritical reactivity measurements, Albuquerque, New Mexico Conference 8508/05, 1985

12 Mihalczo J T, Blakeman E D, Rogan G E, et al. Dynamic subcriticality measurements using the sup252Cf-sourcedriven noise analysis method[J]. Nuclear Science & Engineering, 1990, 104(4): 314–338

13 李建勝, 張翼, 金宇, 等.252Cf快裂變室研制[J]. 核電子學(xué)與探測(cè)技術(shù), 2001, 21(4): 264–267

LI Jiansheng, ZHANG Yi, JIN Yu, et al. Development of252Cf fast fission chamber[J]. Nuclear Electronics & Detection Technlogy, 2001, 21(4): 264–267

14 Mihalczo J T, Mattingly J K, Valentine T E, et al. Source options for nuclear weapons identification system[R]. ORNL/TM-13025, 1995

15 Mihalczo J T. The use of252Cf as a randomly pulsed neutron source for prompt neutron decay measurements[J]. Nuclear Science and Engineering, 1970, 41(1): 296–300

CLCTL277

Simulation study of neutrons time-correlated coincidence count for uranium components based on252Cf source-driven noise analysis method

ZHOU Mi1,2WEI Biao1YANG Fan1FENG Peng1QIAO Liang1

1(Key Laboratory of Opto-Electronic Technology and System,Ministry of Education,Chongqing University,Chongqing 400044,China)
2(School of Optoelectronic Information,Chongqing University of Technology,Chongqing 400054,China)

Background:For the intensity of high-enriched uranium’s neutrons from fission is weak, the active detection method is commonly adopted. Purpose: The quality attribute of uranium components is analysed. Methods: According to the theory of neutron source-driven noise analysis method and mass measurement problem of uranium material or uranium component, Monte Carlo simulation is used to study the mass measurement of metal uranium sphere components with the same geometry and enrichment and the different masses. The time-correlation coincidence count distributions to different uranium components are obtained. Results: The source-driven time correlation coincidence measurements can provide quantities, time-dependent coincidence distributions between two detectors, which can be related to the mass of uranium components. Conclusions: By studying on the sample (radius: 4.67 cm, mass: 8 kg), it realized mass measurement to metal uranium sphere components with different masses.

Source-driven noise analysis method, Uranium component, Neutron signal, Monte Carlo simulation, Time-correlation coincidence count

TL277

10.11889/j.0253-3219.2013.hjs.36.060202

國(guó)家自然科學(xué)基金(61175005)資助

周密，女，1979年出生，2004年于重慶大學(xué)獲碩士學(xué)位，專業(yè)：核信號(hào)計(jì)算機(jī)仿真研究

譙梁，E-mail: qiaoycwh@126.com

2012-12-11，

2013-04-03