高计数率下核脉冲信号容易发生堆积，将堆积脉冲直接舍弃的方法会导致能谱计数率的下降。本文用硅漂移探测器对X光管激发的铜合金样品进行测量，采用14位20 MHz的ADC采样获取数据，在Matlab上对输出信号进行冲激成形以识别堆积脉冲，应用双指数表达式对堆积脉冲进行重构并进行幅值提取计入能谱中。结果表明，与传统的直接舍弃堆积脉冲相比，该方法提高了能谱计数，提升率约为21.37%。

双指数信号，冲激成形，梯形成形，堆积脉冲识别，脉冲重构

Houhong Zhou, Yi Liu, Chao Li

The College of Nuclear Technology and Automation Engineering, Chengdu University of Technology, Chengdu Sichuan

Received: Dec. 13^th, 2023; accepted: Dec. 25^th, 2023; published: Jan. 31^st, 2024

Nuclear pulse tend to pile up at high count rate. The traditional method is to discard the pile-up pulse, which will decrease the energy spectrum count rate. In this paper, silicon drift detector is used to measure the copper alloy excited by X-ray tube, 14-bit 20 MHzADC is used to obtain the data, impulse pulse-shaping is performed on Matlab to identify the pile-up pulse of the output signal, and the pile-up pulse is reconstructed by double-exponential expression and the amplitude is extracted and included in the energy spectrum. The results show that, compared with directly rejecting the pile-up pulse, this method improves the energy spectrum count by 21.37%.

Keywords:Double-Exponential Signal, Impulse Shaping, Trapezoidal Shaping, Pile-Up Pulse Discrimination, Pulse Reconstruction

Copyright © 2024 by author(s) and beplay安卓登录

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

硅漂移探测器(Silicon Drift Detector, SDD)是由高阻硅制成的X射线探测器，相比于Si-PIN等探测器，SDD器件采用侧向耗尽结构设计，电荷输运电场与探测器耗尽电场相对独立，电荷收集极电容与探测器面积无关，使收集极电容非常小，因而具备优异的能量分辨率和计数率性能 [ 1 ] 。

由于探测器输出的核脉冲具有一定的宽度与较长的下降沿，一个脉冲信号宽度内可能相继探测到多次核事件，这会导致信号堆积，引起脉冲幅度的畸变，影响到脉冲信号幅度的有效提取，进而影响到核谱仪的能量分辨率 [ 2 ] [ 3 ] 。堆积脉冲分离常用的方法是将具有长拖尾的双指数信号通过卷积成形为具有一定宽度的信号，减小或者消除前一个脉冲拖尾对后一个脉冲幅值的影响。例如高斯成形、梯形成形(三角成形)、尖顶成形等方法 [ 4 ] [ 5 ] [ 6 ] [ 7 ] [ 8 ] ，成形之后可对脉冲直接进行幅值提取，其中梯形成形算法具有处理过程简单，成形速度快，成形后梯形脉冲宽度窄，顶部平坦，成形参数可调，线性良好等特点被广泛应用 [ 9 ] ，但是对于成形之后仍然堆积的脉冲信号，上述的成形方法无法对其进行分离，因此需要使用其他方法对该种堆积脉冲进行处理。

本文使用Matlab对ADC采集的硅漂移探测器输出信号数据进行处理，采用冲激成形识别堆积脉冲、梯形成形进行幅值提取，应用核信号的双指数模型对堆积脉冲进行分离并重构，以提高能谱的计数率和分辨率。

硅漂移探测器(SDD)的输出信号近似为脉冲信号或单时间常数的指数衰减信号 [ 10 ] 。随着数字系统的发展，高速ADC的应用降低了半导体探测器系统中的噪声，提高了采样精度。因此核脉冲的上升时间可以被精确地探测和描述，SDD的输出信号应该近似为具有一定上升沿的双指数信号 [ 11 ] 。

双指数信号数学表达式：

V ( t ) = { A ⋅ ( e − t τ 1 − e − t τ 2 ) , t ≥ 0 0 , t < 0 (1)

式中， A 为脉冲振幅， τ 1 和 τ 2 分别为信号的慢时间常数和快时间常数。

双指数梯形成形算法 [ 12 ] 是一种将ADC采集得到的双指数脉冲信号转换成梯形脉冲信号的算法，原理如图1所示。理想的梯形脉冲输出信号表达为：

V t r ( t ) = ( A t r ) ⋅ [ t ⋅ u ( t ) − ( t − t r ) ⋅ u ( t − t r ) − ( t − t t ) ⋅ u ( t − t t ) + ( t − t d ) ⋅ u ( t − t d ) ] (2)

V t r ( t ) 为等腰梯形脉冲； u ( t ) 为单位阶跃函数； t r 为梯形脉冲上升结束点； t t 为梯形脉冲平顶结束点； t d 为梯形脉冲下降结束点。

探测器输出信号为双指数信号，其表达式为：

V i ( n ) = A ⋅ ( e − n ⋅ T S τ 1 − e − n ⋅ T S τ 2 ) (3)

式中， A 为脉冲振幅， τ 1 和 τ 2 分别为脉冲的下降部分和上升部分， T S 为ADC的采样周期， n 为采样点序列。

图1. 理想双指数信号与理想梯形成形信号

分别将 V i ， V t r 离散后进行 Z 变换得到：

V i ( z ) = n r ⋅ ( α − β ) ⋅ ( 1 − z − 1 ) 2 (4)

V t r ( z ) = ( 1 − z − n r − z − n t + z − n d ) ⋅ ( 1 − α ⋅ z − 1 ) ⋅ ( 1 − β ⋅ z − 1 ) (5)

式中 n r = t r T S 决定梯形脉冲上升边长度； n t = t t T S 决定梯形脉冲平顶宽度； n d = t d T S 决定梯形脉冲下降长度； α = e − T S τ 1 ， β = e − T S τ 2 。

将上式相除后再作Z反变换即可得出双指数梯形成形算法时域递推式：

V o ( n ) = 2 ⋅ V o ( n − 1 ) − V o ( n − 2 ) + [ ( V i ( n ) − V i ( n − n r ) − V i ( n − n t ) + V i ( n − n d ) )     − ( α + β ) ⋅ ( V i ( n − 1 ) − V i ( n − n r − 1 ) − V i ( n − n t − 1 ) + V i ( n − n d − 1 ) )     + α ⋅ β ⋅ ( V i ( n − 2 ) − V i ( n − n r − 2 ) − V i ( n − n t − 2 ) + V i ( n − n d − 2 ) ) ] / [ n r ⋅ ( α − β ) ] (6)

式中 V i ( n ) 输入双指数脉冲序列；n为采样序列号； V 为输出梯形成形后的序列。

本文用AmpTek公司生产的型号为XR-100SDD半导体探测器对铜(Cu)标准样品进行测量，采用型号为KYW2000A-Ag X光管作为激发源，光管电压设为20 kV，光管电流设置为78.4 uA，测量时间为100 s，采用14位20 MHz的ADC对输出信号进行采集，输出核信号最大幅度为16,384，时间分辨率为50 ns，将所得的数据导入到Matlab中进行处理，对多个非堆积脉冲的时间常数进行拟合求平均值，求得双指数信号的快慢衰减时间常数 τ 1 = 974.595 n s o ， τ 2 = 121.755 n s 。采用梯形成形算法对输出信号进行成形并提取幅值，梯形信号的上升时间设置为1.25 us，平顶时间设置为0.5 us，脉冲宽度为3 us，采用冲激成形对脉冲堆积进行判断并记录堆积脉冲位置，最后对堆积脉冲应用双指数信号模型进行重构，将堆积脉冲重新进行幅值提取后计入到能谱中，具体过程如下所示。

本文根据SDD半导体探测器的双指数输出信号模型，在Matlab上对ADC采集到的数据进行处理，对信号进行梯形成形幅值提取、冲激成形识别堆积脉冲，应用双指数信号模型对被丢弃的堆积脉冲进行重构。结果显示，冲激成形能够准确识别堆积脉冲的位置，重构之后的脉冲信号能够较好的还原前一个脉冲的形状与幅值，并且对第二个脉冲信号的幅值也能较好地还原；与传统的直接舍弃堆积脉冲相比，该方法提高了脉冲计数率。当脉冲未发生堆积时，无需对脉冲进行重构；当脉冲发生前沿堆积时，舍弃该脉冲；当脉冲发尾堆积时，对脉冲进行重构。经过计算，该方法使得能谱计数率提高了21.37%，可以用于堆积脉冲的处理。

国家自然科学基金(12075038)。

周厚泓,刘 易,李 超. 基于双指数拟合的核信号堆积脉冲重构方法Nuclear Pile-Up Pulse Reconstruction Method Based on Double-Exponential Fitting[J]. 核科学与技术, 2024, 12(01): 63-69. https://doi.org/10.12677/NST.2024.121007