陈 炼 梁福田 王宝琛 金 革

基于球谱仪的中子剂量率仪设计

陈 炼 梁福田 王宝琛 金 革

（核探测与核电子学国家重点实验室 中国科学技术大学 合肥 230026）

研制了一套基于球谱仪的中子剂量率仪。采用单个慢化球及位置灵敏计数器的探测结构，通过单次测量即可计算出中子剂量等信息。在基于电荷分配法的读出电路中，通过对电荷灵敏放大器及滤波成型电路的参数进行优化，使探测器系统的平均位置分辨达到6mm。利用硬件寻峰等数字信号处理技术，使数据获取系统的最高计数率达到200 k·s−1，满足剂量率仪实时性的要求。

中子剂量率仪，位置灵敏计数器，电荷分配法，数字脉冲处理

20世纪60年代，Bramblett等[1]提出了利用多球系统测量中子能谱的方法。多球谱仪由一系列中心放置热中子探测器的不同直径慢化球组成，慢化球的最小直径增量为12.7mm。根据测量计数和能谱响应函数，就可以通过解谱算法求解出中子能谱。由于中子能量响应峰值与慢化球直径之间并没有线性对应关系[2]，因此在不同辐射场中使用时，需要的慢化球数量及大小也不尽相同，这使得多球谱仪使用起来十分不便。

Toyokawa[3]发展了基于位置灵敏计数器的单球中子谱仪。该谱仪通过位置灵敏探测器将一个慢化球分成不同的计数区，仅需一次测量就可以获得与BSS谱仪近似的能谱结果，大大简化了谱仪的使用复杂度。由于基于球型能谱仪的中子剂量率仪具有近似各向同性的能量响应及较宽的能量测量范围，其给出的剂量当量值也更为可信，因而在中子辐射防护领域得到了广泛应用[4]。

在能谱型中子剂量率仪中，通过测量慢化球不同位置处的中子注量，求解出辐射场的中子能谱，再利用注量剂量转换系数就可以计算出中子周围剂量当量，计算结果的准确度依赖于能量区间的划分和几何模型的精确描述。为此我们设计了基于单球中子谱仪的能谱型中子剂量率仪系统，利用高位置分辨实现对慢化球的准确区域划分，以使中子剂量的测量更为准确。

1 基于球谱仪的中子剂量率仪结构

基于球谱仪的中子剂量率仪结构如图1所示。

图1 中子剂量率仪结构示意图Fig.1 Schematic diagram of the spectrometric neutron dosimeter.

中子剂量率仪由探测器系统和电子学系统两部分构成。其中探测器系统采用直径30cm的聚乙烯球作为慢化体，并按互相垂直的方式将三根位置灵敏计数器安装在慢化球内，计数器相互间的轴心间距为7mm。计数器的有效长度为250mm，外径为12.7mm，管壁材料为304不锈钢。计数气体为He-3，气压为8×105Pa。电子学系统包括电荷灵敏放大器及数据获取系统等。

2 电荷分配法原理

电荷分配法原理如图2所示。其中Cb为探测器的隔直电容，Cf、Rf为反馈电容与反馈电阻。

反应产生的带电粒子经电离、倍增后被两端的阳极丝收集。根据电荷分配原理，两端收集的电荷与入射点到两端的电阻成反比。这些电荷经电荷灵敏放大器放大后，输出一定幅度的电压信号，则中子入射位置与两端输出电压信号关系可表示为：

式中，L为探测器的有效长度；x为中子入射点位置；R为探测器阳极丝的总电阻；r1、r2为入射点到阳极丝两端的电阻；q1、q2为阳极丝两端收集到的电荷；V1、V2为前端电路输出电压信号的幅度。考虑在非理想情况下，电荷灵敏放大器的输入阻抗、隔直电容的等效阻抗以及连接线缆等都会对电荷分配造成影响，假设这些因素带来的额外阻抗为ε，则此时两端输出电压的关系为：

即两端输出电压与中子入射位置仍具有线性关系，但会受到ε的影响。ε越大，两者线性关系的斜率就越小，即位置分辨的动态范围就越小；而R越大，ε所带来的影响就越小。为获得较好的位置分辨，在选择探测器时，阳极丝应具有较高的电阻率；而当探测器选定时，则应优化前端电子学尤其是电荷灵敏放大器的设计。

3 电子学系统设计

3.1 电荷灵敏放大器设计

电荷灵敏放大器原理如图3所示，其中Cd、Cs、Cin分别为探测器的结电容、分布电容和放大器的输入电容，A0为放大器的开环增益。

图3 电荷灵敏放大器结构示意图Fig.3 Schematic diagram of charge-sensitive amplifier.

当探测器输出电荷Qin时，在放大器输入端形成的电压信号Uin可表示为：

当(1+A0)Cf＞＞Cd+Cb+Cs+Cin满足时，放大器输出信号幅度为：

其反馈系数为：

在电荷灵敏前放中，隔直电容Cb会带入电荷平衡效应，从而影响探测器位置分辨的动态范围，Cb越小影响越明显[5]；Cb选择过大，则会使反馈系数变小。反馈电容Cf决定了电荷放大的灵敏度，其容值越稳定，则输出信号的幅度也能保持恒定，同时小容值的Cf能获得更大的输出信号幅度；但Cf选择过小，则反馈系数就越小，电路的稳定性就会变差。反馈电容Cf通过反馈电阻Rf进行放电，同时放电的时间常数也决定了系统的最大计数率。为减小弹道亏损的影响，Rf应选择较大的电阻，但Rf越大，放电时间常数就越大，会导致系统的最大计数率变低。放大器的开环增益决定了电路的工作状态，运放的输入阻抗越大、增益带宽积越宽，电路的稳定性就越好。

综合考虑信噪比、反馈深度、最大计数率等因素，我们选择了100pF的隔直电容、1pF的反馈电容、10MΩ的反馈电阻以及带宽为1.6GHz的OPA657实现电荷灵敏放大器设计。

3.2 滤波成形电路设计

从电荷灵敏放大器输出的信号通常混杂着较多的噪声，因此需要通过滤波成形电路来提高信噪比。根据最佳滤波原理，可以先用白化滤波器将噪声变为白噪声，再利用匹配滤波器获得最佳信噪比。由于在实际电路中白化滤波器后的匹配滤波基本是无法实现的，因此滤波成型电路采用设计如图4所示，使用CR高通滤波器来实现白化，并使用有源滤波器实现多级RC积分电路来近似匹配滤波器。该电路时间常数选择为1μs，在获得较好的信噪比的同时，输出波形也接近于高斯波形，方便幅度分析使用。

图4 滤波成型电路结构示意图Fig.4 Schematic diagram of filter and shaping circuit.

3.3 数据获取系统设计

为了及时了解中子周围剂量当量的变化，需要中子剂量率仪能实时给出中子能谱、总注量及总剂量等相应信息。而基于传统多道模式的采集系统由于死时间较长，使系统计数率受到限制，为此我们设计了基于FPGA (Field Programmable Gate Array)的高计数率数据获取系统，以满足剂量率仪数据处理实时性的要求。

数据获取系统的结构如图5所示。由滤波成形电路输出的信号，经主放放大后由模拟数字转换器(Analog-to-digital converter, ADC)进行数字化，数字化后的数据直接送入FPGA。当中子击中位置灵敏探测器时，会在探测器两端同时产生信号，因此通过时间符合逻辑就可以将偶然触发的单端信号剔除。在完成事例选择后，利用快速硬件寻峰算法找出对应的两个信号峰值，并计算出粒子入射位置。然后根据事先确定的区域划分方法，对不同区域的粒子入射数做统计。统计完成的数据由USB总线传输至上位机，利用解谱算法求解出中子能谱，进而给出中子剂量等信息。利用FPGA内部的硬件资源，使信号处理过程与信号采集同步完成，数据获取系统并不产生额外的死时间，从而保证了剂量率仪实时性的要求。同时系统还通过FPGA控制高压模块，为探测器提供工作电压。

图5 数据获取系统结构示意图Fig.5 Block diagram of data acquisition system.

4 电子学性能测试

4.1 计数率测试

使用Agilent 33250A型信号发生器进行计数率测试。测试结果表明，在成型滤波电路的时间常数选择为1μs时，数据获取系统可以达到200 k·s−1的最高计数率。由于电荷灵敏放大器的放电时间常数为10μs，探测器的最大计数率小于100 k·s−1，因此系统完全能满足数据实时获取的需求。

4.2 一致性和信噪比测试

我们使用信号发生器对电子学系统的一致性和信噪比进行了测试。由信号发生器输出固定幅度的方波经电容微分后转换为稳定幅值的电流脉冲作为测试信号，同时输入两个前置放大器后，由数据获取系统采集两路信号峰值，并进行统计。测试结果如图6所示，经高斯拟合后两通道峰值位置分别为7110和7073道，其相对偏差为0.5%，半高宽分别为46.61和37.56，计算得到信噪比分别为51.11dB、52.94dB。

图6 固定输入下双通道输出幅度谱Fig.6 Amplitude spectra of two channels with fixed input.

4.3 位置分辨测试

我们使用一台小型加速器中子源对系统进行了位置分辨率测试。测试时，14MeV的中子束通过100mm厚的聚乙烯慢化至热中子能区。He-3探测器放置在两块包裹镉套的含硼聚乙烯及铅砖屏蔽体之后。两块屏蔽体之间留有约1mm的缝隙。测量时缝隙距离He-3管的左端距离分别为15mm、35mm、60mm、90mm、125mm、160mm、190mm、215mm及235mm。测量结果如图7所示，在探测器的有效长度内，中子入射位置与测量结果具有较好的线性度。探测器中心的位置分辨达到了3.5mm，由于电荷平衡等效应的影响，两端的位置分辨慢慢变差，其平均分辨为6.1mm，能满足对应BBS谱仪最小12.7mm增量的区域划分需求。

图7 位置分辨测试结果 (a) 相对位置-计数分布图，(b)中子入射位置与相对位置关系图，(c) 相对位置分辨示意图Fig.7 Test result of position resolution. (a) Relative position and counter, (b) Incident position and relative position, (c) Incident position and position resolution

5 结语

本文设计了基于位置灵敏探测器的能谱型中子剂量率仪。通过优化电荷灵敏放大器及滤波成型电路的参数选择，使系统具有较高的信噪比及较好的位置分辨。测试的结果表明电荷灵敏放大器具有超过50dB的信噪比，且不同通道间的一致性差异小于0.5%。位置灵敏探测器的位置分辨具有较好的线性度，平均分辨达到6.1mm，满足了中子剂量计算中对慢化球空间分区最小12.7mm增量的需求。而基于FPGA的数据获取系统可以达到200K的最高计数率，也满足了剂量率仪数据处理实时性的需求。

1 Bramblett R L, Ewing R I, Bonner T W. A new type of neutron spectrometer[J]. Nuclear Instruments and Methods, 1960, 9(1): 1-12

2 Mares V, Schraube G, Schraube H, et al. Calculated neutron response of a Bonner sphere spectrometer with3He counter[J]. Nuclear Instruments and Methods in Physics Research A, 1991, 307: 398-412

3 Toyokawa H, Uritani A, Mori C, et al. Neutron spectrometer with position-sensitive proportional counters[J]. Nuclear Instruments and Methods in Physics Research A, 1996, 381: 481-487

4 Thomas D J, Alevra A V. Bonner sphere spectrometers-a critical review[J]. Nuclear Instruments and Methods in Physics Research A, 2002, 476: 12-20

5 黄土琛, 宫辉, 邵贝贝. 位置灵敏3He管中子探测器前端电子学设计[J]. 清华大学学报(自然科学版), 2012, 52(12): 1731-1735

HUANG Tuchen, GONG Hui, SHAO Beibei. Design of the front-end electronics for linear position-sensitive3He neutron detector[J]. Journal of Tsinghua University (Natural Science Edition), 2012, 52(12): 1731-1735

CLC TL82

Design of the neutron dosimeter based on neutron sphere spectrometer

CHEN Lian LIANG Futian WANG Baochen JIN Ge
(State Key Laboratory of Particle Detection and Electronics, University of Science and Technology of China, Hefei 230026, China)

Background: In neutron radiation protection applications, the spectrometric neutron dosimeter has almost isotropic response and can cover the energy from thermal to GeV neutrons, and the radiation filed information such as neutron spectrum and neutron dose can be calculated more accurately. Purpose: A spectrometric neutron dosimeter which has high precision and high counting rate is developed. Methods: By optimizing the parameters of the readout circuit, the position resolution of the detector system is improved and the geometric division of the spherical moderate is more accurate. Using digital pulse processing technology, the dead-time of the system is reduced. Results: The position resolution of the neutron dosimeter has good linearity, and the average resolution is about 6 mm. The high counting rate of the data acquisition system is 200 k·s−1. Conclusion: We developed a neutron dosimeter based on neutron sphere spectrometer which has high position resolution and high counting rate. The performance of the dosimeter is expected to be improved for further optimization design after the test in the radiation field.

Neutron dosimeter, Position-sensitive counter, Charge division, Digital pulse processing

TL82

10.11889/j.0253-3219.2015.hjs.38.010402

项目(No.11375179)资助

陈炼，男，1980年出生，2012年于中国科学技术大学获博士学位

金革，E-mail: goldjin@ustc.edu.cn

2014-09-04，

2014-10-30