在20世纪50年代，多道能谱分析技术就在科学、工业、医疗领域占有重要地位。随着数字信号处理与高速模数转换(Analog-to-digital Converter,ADC)技术的发展，数字化多道能谱分析技术已获得广泛应用[1-2]。由于受电源噪声、温度、探测器漏电流和较差的极零相消处理等因素的影响，探测器输出核脉冲信号的基线会发生漂移；另外，高放射性条件下产生的脉冲堆积也会引起核脉冲信号的基线漂移。基线漂移使得有效信号叠加在可能不为零的基线上，干扰真实脉冲幅度的提取，导致测量系统的能量分辨率损失。因此，必须采用有效的方法恢复核脉冲信号的基线，以改善测量系统的能量分辨率[3-10]。

从模拟方法到数字方法，众多机构和学者对基线恢复技术进行了大量理论与方法研究。模拟多道能谱仪采用的基线恢复器虽能有效消除基线漂移和涨落，但易受电子线路和不确定因素的影响。当前广泛使用的数字基线恢复方法有 Kalman滤波器法[5]、IIR滤波器法、最小噪声计权法[6]、零面积滤波恢复法[7]等。如，李彦波等[8]针对相邻采样点取平均值算法的不足，提出了一种基于补偿最小二乘法的基线估计算法；Kafaee[9]提出了双极尖状成形来减轻基线漂移和缓解堆积方法。这些方法运算量较大，不适合现场可编程门阵列(Field-Programmable Gate Array, FPGA)编程；李伟男等[10]提出的动态实时求平均估计基线的方法，易于在FPGA中实现，在高计数率测量条件下，该方法能够改善测量系统能量分辨率，且具有较高的准确性、实时性和灵敏度。以上各基线恢复方法为纯数字算法，没有实际解决探测器输出核脉冲信号基线漂移对前端电路的影响，以及在高放射性测量条件下因核脉冲信号基线漂移引起的谱线漂移和能量分辨率损失问题。本文设计的基于基线自动恢复的数字化多道能谱仪可以很好地改善这些问题。

1 探测器输出信号的基线漂移

闪烁体-光电倍增管系统在核辐射探测器中有广泛应用。它的工作原理是：入射射线与闪烁体原子发生相互作用，使其发生电离或激发，闪烁体原子在退激过程中发射出波长近似可见光的闪烁光子；当闪烁光子到达光阴极后，与光阴极发生光电效应，闪烁光子转换为光电子；聚焦极作用使光阴极产生的光电子尽可能多地打到光电倍增管的第一打拿极上，以产生更多的电子；这些电子在打拿极的作用下加速、倍增，最终在阳极形成电信号，该信号的幅度与入射射线的能量成正比。由于光电倍增管的最末三级打拿极中的次级电子数量已经很大，在相应分压电阻上产生的脉冲电流也很大。因此，通常分压电阻并联旁路电容来改善脉冲电流过大时对极间电压的影响。另外，为了避免空间电荷效应，最末两级打拿极间的电位差需设置高一些。对于高计数率场合，可使用稳定性更好的晶体管射极跟随器分压网络[11-15]，来减小分压器效应的不利影响。

对核能谱测量系统中的基线漂移现象分析发现，后级处理电路引入的基线漂移较小，出现基线漂移的源头在探测器光电倍增管的分压网络中，探测器输出信号的基线漂移量与被测对象的放射性强度有关。为观察在不同放射性条件下探测器输出核脉冲信号基线漂移现象，采用法国圣戈班公司的ø75 mm×75mm LaBr3(Ce)探测器，分别在无放射源、仅137Cs源、137Cs+60Co源和137Cs+60Co+22Na源条件下记录探测输出核脉冲信号的基线道址，实验结果见图1。由图1可以看出，无放射源时，探测器输出核脉冲信号的基线道址集中在 21道附近；仅137Cs源时，基线道址在20道附近，与环境本底测量时相比漂移了1道；当137Cs+60Co源时，基线道址在 18道附近，与环境本底测量时相比漂移了 3道；当137Cs+60Co+22Na源时，基线道址在12道附近，与环境本底测量时相比漂移了9道。可以看出，随着测量对象放射性强度的增加，探测器输出核脉冲信号的基线漂移变大。

3.2.2 鼠和狗的主动脉腔静脉瘘模型 模型优点是该模型适用于所有动物，具有高度可复制性；缺点是这类疾病在人类身上发生率很低[23-26]。

  

图1 不同放射性强度下探测器输出信号的基线道址Fig.1 Baseline channel of detector output signal underdifferent radioactive intensities

2 数字化多道能谱仪设计

2.1 系统原理

前端电路的原理如图4所示。第一级（图4第Ⅰ区）为放大倍数在10倍以内的固定增益电路；第二级（图4第Ⅱ区）为偏置调节电路，反向端BL_ADJ信号来自后级带基线调整量的 DAC输出，用于基线调节；第三级（图4第Ⅲ区）为程控增益放大电路，由MCU控制实现核脉冲信号幅度的线性放大；第四级（图4第Ⅳ区）选用低失真差分ADC驱动器AD8138来抑制共模噪声及干扰，将单端模拟核信号转换成差分输出，使后级ADC性能达到最佳；第五级为高速 ADC转换电路，选用美国 Analog Devices公司的2×107 s-1 12位多级差分流水线结构的模数转换器AD9235（图4第Ⅴ区）。

随着国家对节能工作的越来越重视，随着电气设备生产技术的改进和节能减排意识的提高，电气工程的节能工作有了一定的成效。为了使电气工程进一步提高节能的效果，以适应社会需要，满足国家可持续发展的要求，应该加大节能力度，将节能的理念贯穿于整个用电过程。从生产到管理，从输送到使用，都应该以节能为目的，为合理利用资源做贡献。

  

图2 数字化多道能谱仪框图Fig.2 Block diagram of digital multi-channel spectrometer

2.2 最小平均值基线估计算法

MCU选用STM32F103处理器，它的主要任务是建立数据通信和控制基线调节。MCU中设有基线漂移阈值寄存器，允许基线在一定范围内波动。MCU获取FPGA中的基线估计值后，先将估计值与基线漂移阈值寄存器内的设定值进行比较，当新基线值的漂移量大于设定阈值时，触发后续的基线调整进程；MCU将基线估计值运算后得到带调节量的调整值，该调整值受调整精度阈值的限定；当调整值大于调整精度阈值时，MCU将调整值通过DAC输出至信号调理电路的基线调整端 BL_ADJ（图4第Ⅲ区）进行核脉冲信号的基线调节；当调整值小于调整精度阈值时，此时的基线值认为是0mV，不做基线调节。图6为基线恢复流程图。调整精度阈值由高速ADC决定，本文选用的12位高速ADC，输入信号满量程为2V，则最小刻度单位为 2/212≈0.5mV，故本文中的调整精度阈值为0.5mV。短时间内的脉冲基线漂移量较小，为了较准确地估计一段时间内的基线水平，减小最小平均基线估计误差，应合理设置最小值采样点数M和增加最小值个数N，即以增加采样时间的方式获取基线估计值。本设计中采用的基线恢复调节周期为1s。

 

式中：为最小平均基线估计值；Y为最小值的个数，为便于FPGA移位计算，Y通常取2的n次幂。

由图9可以看出，仅137Cs源时，特征峰峰位在435道，计算得到能量分辨率为3.03%；在137Cs+60Co+22Na源时，137Cs特征峰位于420道，此时的能量分辨率为 3.27%。与仅 137Cs源的测量结果对比，在高放射性条件下，137Cs特征峰向左漂移了15道。

  

图3 基线估计流程图Fig.3 Flow chart of baseline estimation

2.3 前端电路

数字化多道能谱仪(Digital Multi-Channel Analyzer, DMCA)主要由探测器、前端电路、ADC、数字信号处理单元、MCU和计算机组成[2,12]，系统原理框图如图2所示。前端电路的主要功能是对探测器输出核脉冲信号进行线性放大，使其幅度在ADC的采样电压范围；另外，前端电路还需完成核脉冲信号的基线调节工作。ADC采集后的数字信号被送入数字脉冲处理单元进行数字处理，完成基线估计、滤波成形、堆积识别、幅值提取等任务。MCU一方面通过SPI总线获取FPGA中的数字基线估计值，经过基线漂移判别算法计算出调节量，控制数模转换器(Digital-to-analog Converter, DAC)输出模拟调整信号，送至前级偏置调整电路恢复基线；另一方面通过 CAN通讯端口将前级数字脉冲处理后的数据上传至计算机进行数据分析与显示。计算机系统的主要功能有：获取测量数据、完成能量刻度、谱数据分析与处理、测量数据存储显示等。

  

图4 前端电路原理图Fig.4 Schematic diagram of the front-end circuit

2.4 数字脉冲处理

图5 是FPGA中的数字处理单元功能框图，图中细实线箭头表示控制信号，粗实线箭头表示总线数据。采样控制模块控制高速ADC，将调理后的核信号转换为数字信号，并送至先进先出(First In First Out, FIFO)缓存。基线估计模块从FIFO读取数据，运行最小平均基线估计算法，该模块与后级 MCU控制的 DAC和前级偏置调整电路共同完成基线恢复。梯形成形分为快成形和慢成形两个结构相同的模块，快成形模块将脉冲成形为宽度较小的梯形脉冲，用作为慢成形的控制信号；快成形后的数据进入阈值去噪模块，剔除噪声后送入堆积判弃模块，一方面用作慢成形时的幅值提取的控制信号，另一方面产生死时间信号供死时间记录。慢成形后的数据，进入幅值提取模块，然后根据幅度大小存入谱线RAM中形成谱数据。图中的FIFO、波形RAM、谱线RAM均为双端口，输入输出端口可以独立进行读写操作。借助多路选择器，上位机可自由选择读取 ADC采样原始脉冲数据、梯形成形后脉冲数据、多道谱数据，方便对核信号的分析与处理。逻辑控制模块用于协调各模块之间的数据同步，参数寄存器用于存储系统的配置参数，如启停状态、脉冲幅值的放大倍数、阈值去噪模块的阈值、梯形成形参数等。

  

图5 FPGA中的实时处理功能框图Fig.5 Block diagram of the real time processing unit in FPGA

2.5 基线恢复的实现

在数字化多道能谱仪中，核脉冲信号的基线估计在FPGA中完成。由于FPGA不擅长除法运算，为提高系统的运行效率，应尽量采用简单算法。实际应用中，短时间内脉冲基线的漂移量较小，因此设计了一种最小平均值算法来求取一段时间的基线估计值，方法如下：高速 ADC对脉冲信号连续采样，FPGA取X个连续采样值并求出最小值，记为Ai（i=1, 2, 3, …, Y），然后重复Y次上述操作，得到最小值数列，求该数列的平均值并将其作为基线估计值，算法如下：

图1为文献[11]采用的实验模型, 是由680 mm×380 mm的基础大平板和可改变角度的三角楔体两部分组成. 在表面沿纵向布有3条测量线, 分别是I线: 沿模型的对称线布置; II线: 距模型对称线58 mm; III线: 距离模型对称线80 mm, 由于三维楔体的横向总宽度为140 mm, 所以III线实际上位于楔体侧面平板上, 距离楔体侧面10 mm.

  

图6 基线自动恢复流程图Fig.6 Flow chart of auto-baseline restoration

3 实验与测试

3.1 基线实时调整测试

为验证算法的准确性，分别在137Cs源和137Cs+60Co+22Na源两种情况下进行能谱测量，137Cs+60Co+22Na是为了模拟高放射性条件。在不采用基线自动恢复算法情况下，不同放射性强度测量结果如图 9所示。

  

图7 137Cs源能谱Fig.7 Energy spectrum of 137Cs source

  

图8 初始偏置值小于(a)和大于(b)实际值时的基线自动调整测试Fig.8 Auto-baseline regulation test with the set value more (a) and less (b) than the actual values

3.2 基线自动恢复测试

对 137Cs进行多次测量，系统的能量分辨率和特征峰所在道址均稳定，图7为任选出的一组谱线。特征峰661keV所在的峰位在435道，能量分辨率为 2.92%。测量条件不变，通过改变探测器输出核脉冲信号的偏置电压，模拟核脉冲信号的基线漂移。图 8(a)为设定的基线值小于实际基线值时测得的能谱图。由图 8(a)可以看出，Ⅰ区谱线低能区域的高计数率异常谱线产生的原因为初始基线值较低，噪声被作为有效信号被记录；Ⅱ区所示谱线区域，谱线向低能方向自动调整；图8(b)为设定的基线值大于实际基线值时测得的能谱图，由图8(b)可以看出，Ⅲ区所示谱线区域，特征峰向高能方向自动调整。综合分析图8(a)、(b)的谱线可得，Ⅱ、Ⅲ区谱线的异常是由于初始基线值与实际基线值相差较大，影响了脉冲信号的幅值，进而使峰位发生漂移，两种测试方式的137Cs特征峰最终均稳定在435道。调整后的基线值可作为后续同等测量条件下的基线初始值。

X为采样点数，为能获取到堆积脉冲的最低幅值，X的取值应大于两个脉冲宽度所对应的采样点数。探测器输出信号基线估计流程如图3所示。

  

图9 未采用基线恢复时不同放射性强度测量能谱Fig.9 Measurement spectrums under different radioactive intensities without baseline restoration

采用基线自动恢复算法后，在相同测量条件下，不同放射性强度测量结果如图10所示。可以看出，在高计数率条件下137Cs的特征峰峰位在434道，计算得到能量分辨率为3.17%，与仅137Cs源时测量结果对比，峰位仅左移了1道，能量分辨率指标下降了0.25%。

许多人问，低压高，但是高压正常是怎么一回事？这种现象算不算高血压呢？需要怎么处理才好呢？今天就来说一说关于低压高的事。

  

图10 采用基线恢复时不同放射性强度测量能谱Fig.10 Measurement spectrums under different radioactive intensities with baseline restoration

3.3 能量线性测试

采用137Cs+60Co+22Na源，测量得到混合能谱，各个放射源的特征峰能量以及在能谱中的道址如表1所示。对这些能量与道址进行线性拟合，得到二者之间的线性关系，如图11所示，其中线性拟合参数R2=0.99999，结果表明：测量系统具有较好的能量线性。

从口岸运行62年的整体来看，二连浩特口岸的进出口贸易工作完成得非常出色，成为了我国北方贸易繁荣发展的重要门户。之所以能够如此出色，与口岸边防力量提供的和平稳定环境密不可分。

 

表1 特征峰峰位与对应道址Table 1 Peak positions of nuclide and their corresponding channels

  

核素Nuclide 能量Energy / keV 道址Channel 22Na 511.00 333.64 137Cs 661.65 434.54 60Co 1173.2 782.05 22Na 1274.5 850.97 137Co 1332.5 888.57

  

图11 道址与能量的线性拟合结果Fig.11 Liner fitting result between channel and energy

4 结语

针对不同放射性强度的测量应用中，闪烁体-光电倍增管探测器输出核脉冲信号的基线发生漂移的问题，提出了数字基线估计与硬件基线恢复相结合的方法，并基于LaBr3(Ce)探测器，设计了具有基线自动恢复功能的数字化多道能谱仪。测试表明：设计系统在 511～1332.5keV 能量范围，能量-道址线性拟合参数 R2=0.99999；在 137Cs+60Co+22Na源条件下，采用基线自动恢复方法可将 137Cs的661keV特征峰稳定在434道，与仅137Cs源时相比，峰位相差1道；同时，采用基线自动恢复可以使系统能量分辨率保持稳定。测量系统能够解决核脉冲信号基线漂移引起的谱线漂移问题。

参考文献

1 Simões J B, Correia C. Pulse processing architectures[J].Nuclear Instruments & Methods in Physics Research,1999, 422(1-3): 405-410. DOI: 10.1016/S0168-9002(98)00992-9.

2 陈亮. 核素识别算法及数字化能谱采集系统研究[D].北京: 清华大学, 2009.CHEN Liang. Nuclide identification algorithm and digital spectra acquisition system[D]. Beijing: Tsinghua University, 2009.

3 Abbiati R, Gatti E, Geraci A, et al. A new digital estimation technique for baseline restoration[J]. Nuclear Instruments & Methods in Physics Research A, 2005,548(3): 507-516. DOI: 10.1016/j.nima.2005.04.072.

4 Geraci A, Rech I, Gatti E, et al. Shared baseline restoration at minimum noise for high resolution spectroscopy[J]. Nuclear Instruments and Methods in Physics Research A, 2002, 482(1-2): 441-448. DOI:10.1016/S0168-9002(01)01509-1.

5 张同锋, 方方, 王敏, 等. 核辐射脉冲幅度分析的基线自适应 Kalman滤波估计[J]. 核电子学与探测技术,2012, 32(2): 203-206. DOI: 10.3969/j.issn.0258-0934.2012.02.020.ZHANG Tongfeng, FANG Fang, WANG Min, et al.Baseline adaptive Kalman filter estimation method for nuclear radiation pulse height analysis[J]. Nuclear Electronics & Detection Technology, 2012, 32(2):203-206. DOI: 10.3969/j.issn.0258-0934.2012.02. 020.

6 肖无云, 魏义祥, 艾宪芸. 多道脉冲幅度分析中的数字基线估计方法[J]. 核电子学与探测技术, 2005, 25(6):601-604.XIAO Wuyun, WEI Yixiang, AI Xianyun. Digital baseline estimation method for multi-channel pulse height analyzing[J]. Nuclear Electronics & Detection Technology, 2005, 25(6): 601-604.

7 Geraci A, Pullia A, Ripamonti G. Quasi-optimum gamma and X spectroscopy based on real-time digital techniques[J]. Nuclear Instruments & Methods in Physics Research A, 2000, 439(2): 378-384. DOI: 10.1016/S0168-9002(99)00897-9.

8 李彦波, 李朝海, 何子述, 等. DMCA 中基于补偿最小二乘法的基线估计算法[J]. 核电子学与探测技术, 2014,34(8): 991-994. DOI: 10.3969/j.issn.0258-0934.2014.08.017.LI Yanbo, LI Chaohai, HE Zishu, et al. A baseline estimation algorithm based on penalized least squares in DMCA[J]. Nuclear Electronics & Detection Technology,2014, 34(8): 991-994. DOI: 10.3969/j.issn. 0258-0934.2014.08.017.

9 Kafaee M, Moussavi-Zarandi A. Baseline restoration and pile-up correction based on bipolar cusp-like shaping for high-resolution radiation spectroscopy[J]. Journal of the Korean Physical Societ, 2016, 68(8): 960-964. DOI:10.3938/jkps.68.960.

10 李伟男, 杨朝文, 周荣. 基于FPGA脉冲幅度分析器的数字化基线估计方法[J]. 核技术, 2015, 38(6): 060403.DOI: 10.11889/j.0253-3219.2015.hjs.38.060403.LI Weinan, YANG Chaowen, ZHOU Rong. Baseline estimation method of digital multi-channel pulse height analyzer based on FPGA[J]. Nuclear Techniques, 2015,38(6): 060403. DOI: 10.11889/j.0253-3219.2015.hjs.38.060403.

11 梁卫平, 胡颖睿, 肖无云, 等. 数字化多道脉冲幅度分析器调理电路设计[J]. 核电子学与探测技术, 2012,32(4): 462-465. DOI: 10.3969/j.issn.0258-0934.2012.04.021.LIANG Weiping, HU Yingrui, XIAO Wuyun, et al.Digital multi-channel pulse height analyze pre-filter circuit design[J]. Nuclear Electronics & Detection Technology, 2012, 32(4): 462-465. DOI: 10.3969/j.issn.0258-0934.2012.04.021.

12 钱建复, 沈庭云. 核辐射剂量学[M]. 北京: 国防工业出版社, 2009: 161-193.QIAN Jianfu, SHEN Tingyun. Nuclear radiation dosimetry[M]. Beijing: National Defense Industry Press,2009: 161-193.

13 曹忠, 黄正德, 陈群. 高计数率下闪烁体-光电倍增管系统增益的稳定[J]. 核技术, 1990, 13(2): 85-90.CAO Zhong, HUANG Zhengde, CHEN Qun. Gain stabilization of the scintillator-photomultiplier counting rates[J]. Nuclear Techniques, 1990, 13(2): 85-90.

14 张思玲, 董成富, 原振东, 等. 低功耗和快速补偿的光电倍增管偏置电路[J]. 核电子学与探测技术, 2000,20(3): 236-239.ZHANG Siling, DONG Chengfu, YUAN Zhendong, et al.Low consumption and fast recovery photomultiplier bases[J]. Nuclear Electronics & Detection Technology,2000, 20(3): 236-239.

15 魏福利, 王培伟, 袁媛, 等. 光电倍增管线性电流扩展研究[J]. 核电子学与探测技术, 2011, 31(9): 986-990.DOI: 10.3969/j.issn.0258-0934.2011.09.012.WEI Fuli, WANG Peiwei, YUAN Yuan, et al. The study of linear current extention of PMT[J]. Nuclear Electronics& Detection Technology, 2011, 31(9): 986-990. DOI:10.3969/j.issn.0258-0934.2011.09.012.