![通讯作者](http://jst.tsinghuajournals.com/CN/rhhtml/images/REcor.gif)
![gonghui02@mails.tsinghua.edu.cn](http://jst.tsinghuajournals.com/CN/rhhtml/images/REemail.gif)
清华大学 工程物理系, 粒子技术与辐射成像教育部重点实验室, 北京 100084
收稿日期:2016-03-22
基金项目:国家自然科学基金资助项目(11475094)
作者简介:黄宇雁(1992—), 女, 博士研究生
通信作者:宫辉, 副教授, E-mail:gonghui02@mails.tsinghua.edu.cn
摘要:针对常见的指数型核脉冲信号,研究了高计数率下的梯形成形算法,特别针对计数率增大所带来的堆积问题进行了研究。提出了一种利用反卷积定位信号到达时刻的算法,有效提升了对堆积脉冲的识别能力,降低了堆积效应对能量分辨率的影响。结合上位机中用LabVIEW编写的采集软件,搭建了一套完整可靠的脉冲幅度分析系统,将梯形成形、堆积判别、能谱构建等一系列功能在现场可编程门阵列(FPGA)中实时实现。将该系统实际应用于使用硅漂移探测器的X射线衍射仪上,获得了良好的结果,在1.8×105 cps(每s计数)的条件下,所得能量分辨率达2.46%。
关键词:实时梯形成形高计数率堆积效应X射线衍射仪
Real time trapezoidal shaping algorithm at high count rates
HUANG Yuyan, GONG Hui
![Corresponding author](http://jst.tsinghuajournals.com/CN/rhhtml/images/REcor.gif)
![gonghui02@mails.tsinghua.edu.cn](http://jst.tsinghuajournals.com/CN/rhhtml/images/REemail.gif)
Key Laboratory of Particle and Radiation Imaging of Ministry of Education, Department of Engineering Physics, Tsinghua University, Beijing 100084, China
Abstract: A trapezoidal shaping algorithm was developed for common exponential nuclear signals and the pile-up effect at high count rate was studied. Deconvolution was used for time-of-arrival estimates to effectively improve the ability to recognize pile-up pulses and reduce the pile-up effect on the energy resolution. A complete pulse height analyzer was built in LabVIEW to realize real time trapezoidal shaping, pile-up recognition and spectrum construction in the field-programmable gate array (FPGA). Measurements with an X-ray diffractometer with a silicon drift detector (SDD) give a 2.46% energy resolution at a count rate of 1.8×105 cps (counts per second).
Key words: real time trapezoidal shapinghigh count ratepile-up effectX-ray diffractometer
梯形成形算法是一种广泛应用于核信号与能谱测量的数字滤波成形方法[1],其本身有一定的噪声抑制能力,在梯形平顶时间大于探测器电荷收集时间的情况下,能够有效避免弹道亏损,且具有脉冲形状对称、脉宽窄、上升时间与平顶时间独立可调等诸多优点[2]。尤其在高计数率情况下,通过减小成形脉冲宽度,梯形成形算法能够达到一定的去除堆积、减小计数率损失的目的。
堆积效应是核信号测量中的普遍问题,也是影响幅度测量和能谱分辨率的关键问题。长期以来,中外****对此问题进行了大量研究,目前提出的方法主要有:减小脉冲宽度以减小堆积概率[3-4]、依据参考波形外推[5-7]、极大似然估计[8-9]、能谱反卷积[10-11]等。在高计数率下,为了在现场可编程门阵列(field-programmable gate array, FPGA)中实时实现梯形成形和堆积判别等功能,减小成形脉冲宽度是一个可行且高效的方法,但此脉宽不能无限减小,需要综合考虑输入信号的上升时间、弹道亏损效应、最终对能谱的能量分辨率要求等。
X射线衍射仪是一种常用的晶体结构分析仪器。本文实验中的射线源为铜的特征X射线,主要存在Kα和Kβ两种特征,对应能量分别为8.05和8.9 keV,由于后续晶体分析实验中只用到Kα峰,故需通过单道上下阈值的方式将Kβ峰去除,防止它对分析结果产生干扰。本文实验系统结构如图 1所示。
![]() |
图 1 X射线衍射实验系统结构 |
图选项 |
本文采用实时梯形成形算法,结合利用反卷积定位信号到达时刻的堆积脉冲识别算法,在FPGA上实现了实时的数字脉冲多道仪的大部分功能,并实际应用于X射线衍射仪上,获得良好的结果。
1 梯形成形算法在FPGA中的实现理想的梯形成形脉冲如图 2所示。
![]() |
图 2 理想梯形脉冲 |
图选项 |
指数输入信号可表示为
${V_i}\left( t \right) = A{{\rm{e}}^{-t/\tau }}u\left( t \right).$ | (1) |
$H\left( z \right) = \frac{{{z^{-1}}\left( {1-d{z^{-1}}} \right)\left( {1 - {z^{ - {n_a}}}} \right)\left( {1 - {z^{ - {n_b}}}} \right)}}{{{n_a}{{\left( {1 - {z^{ - 1}}} \right)}^2}}}.$ | (2) |
$d = {{\rm{e}}^{-{T_{\rm{s}}}/\tau }}, $ | (3) |
${n_a} = {t_a}/{T_{\rm{s}}}, {\rm{ }}{n_b} = {t_b}/{T_{\rm{s}}}, {\rm{ }}{n_c} = {t_c}/{T_{\rm{s}}}.$ | (4) |
观察系统函数H(z)的表达式不难发现,此为一个无限冲激响应(infinite impulse response, IIR)数字滤波系统,具有反馈单元。为了降低系统设计的复杂度,避免量化误差与舍入误差的积累,故选择级联结构来实现。对于线性移不变系统,交换子系统的级联顺序并不会影响最终的输出,但在FPGA中具体实现时,由于存在误差与有限字长效应,系统的分解方式与级联顺序成为必须考虑的问题。
文[13]提出了一种系统的分解方式与级联顺序,并作了较为详细的分析和阐述。本文在此基础上进行了改进,利用更少的资源实现更稳定的梯形成形算法,
H(z)=H1(z)H2(z)H3(z)H4(z).
其中:
$\left\{ {\begin{array}{*{20}{c}}{{H_1}\left( z \right) = 1-{z^{-{n_b}}}, \;\;\;\;\;}\\{{H_2}\left( z \right) = \frac{{1-d{z^{ - 1}}}}{{{n_a}\left( {1 - {z^{ - 1}}} \right)}}, }\\{{H_3}\left( z \right) = \frac{{{z^{ - 1}}}}{{1 - {z^{ - 1}}}}, \;\;\;\;\;\;}\\{{H_4}\left( z \right) = 1 - {z^{ - {n_a}}}.\;\;\;\;\;\;}\end{array}} \right.$ | (6) |
实验表明,此滤波器结构所实现的梯形成形算法,加上简单的空闲时间判定和基线恢复模块,即可准确实现长时测量。
2 堆积识别算法精确提取脉冲幅度依赖于准确判定脉冲的到达时刻。普遍的做法是利用两个通道的梯形成形,分别称为慢通道和快通道[12]。慢通道梯形脉宽较宽,用于幅度提取,以便获得较好的分辨率,但在高计数率下,为了减小堆积效应,此脉宽也不能太大。快通道梯形脉宽很窄,用于定位信号脉冲的到达时刻,取为输入信号上升时间的3~4倍即可。
当检测到相邻两个快通道信号之间间隔小于某一阈值时,即可判定信号产生了堆积。然而,快通道的梯形信号本身也有产生堆积的可能,使得最终获得的能谱产生误差。尤其在高计数率条件下,这种现象不可忽视。因而,本文研究了一种利用反卷积与快通道梯形成形相结合的脉冲定位方法,其原理如下:
归一化输入信号为
${V_i}\left( t \right) = {{\rm{e}}^{-t/\tau }}u(t).$ | (7) |
${V_i}\left( z \right) = \frac{1}{{1-d{z^{-1}}}}.$ | (8) |
$H\left( z \right) = \frac{{X\left( z \right)}}{{{V_i}\left( z \right)}} = 1-d{z^{-1}}.$ | (9) |
$y\left( n \right) = x\left( n \right)-dx(n-1).$ | (10) |
图 3为脉冲到达时刻检测的ModelSim仿真结果,其激励信号为实采探测器输出信号,快通道成形时间为300 ns。可以看出,在输入指数脉冲与快通道梯形脉冲均发生堆积的情况下,利用反卷积算法能够准确地识别出堆积。
![]() |
图 3 堆积识别算法仿真波形 |
图选项 |
3 实验与应用结果本实验中所用晶体为硅粉,探测器为KETEK公司的AXAS-MH150型探测系统,含硅漂移探测器和前放,输出信号的上升时间为100 ns。数字多道部分的ADC型号为AD9255,采样频率80 MHz,FPGA为Xilinx公司的XC6SLX16。
3.1 能谱与分辨率对成形时间从600 ns到1 200 ns、计数率从约25 kcps (counts per second, 每s计数)到约270 kcps进行了能谱采集测试,数字多道设置为4 096道。在成形时间1.2 μs、180 kcps高计数率条件下,所得能谱如图 4所示。对能谱作Gauss拟合可得其均值与半高宽分别为1 429.4与35.164,则分辨率=半高宽/均值=2.46%。能谱均值对应铜的Kα射线能量8.05 keV,则等效的分辨率=8.05 keV×2.46%=198 eV。可见,本实验的能量分辨率达2.46%(198 eV),能够有效分离Kα和Kβ特征。
![]() |
图 4 8.05 keV铜Kα特征X射线能谱 |
图选项 |
不同输入计数率下,不同成形时间条件下所得能谱的能量分辨率曲线如图 5所示。可以看出,同一成形时间下,分辨率随计数率的升高而变差,这是由于计数率的增长使得脉冲堆积效应更加明显,对脉冲幅度的提取误差增大。在相同计数率条件下,梯形成形的上升时间越长,其对高频噪声的抑制越好;平顶时间越长,越能降低由探测器电荷收集时间所带来的弹道亏损[2],故而能量分辨率越好。
![]() |
图 5 不同计数率与成形时间下的能量分辨率 |
图选项 |
3.2 计数率校正核信号通常服从Poisson分布,
$p\left( {x = n} \right) = \frac{{{\lambda ^n}{{\rm{e}}^{-\lambda }}}}{{n!}}.$ | (11) |
$p(\Delta t \ge t) = {{\rm{e}}^{-\lambda }}.$ | (12) |
${p_{{\rm{record}}}} = {{\rm{e}}^{-2\lambda }}, $ | (13) |
$\lambda = {R_{{\rm{real}}}}{t_{{\rm{pile-up}}}}.$ | (14) |
相应的计数率校正公式为
${\rm{ }}{R_{{\rm{real}}}}.{{\rm{e}}^{-2{R_{{\rm{real}}}}{t_{{\rm{pile-up}}}}}} = {R_{{\rm{measure}}}}.$ | (15) |
![]() |
图 6 校正前后脉冲输入计数率与输出计数率的关系 |
图选项 |
4 结论本文探讨了高计数率下的梯形成形算法,重点研究了合理的滤波器结构和脉冲到达时刻定位算法,在FPGA中实现了实时的数字脉冲多道谱仪的大部分功能。在X射线衍射仪系统上配合硅漂移探测器,在成形时间1.2 μs、180 kcps高计数率条件下,所得能量分辨率达2.46%(198 eV)。另外,本文所研究的梯形成形算法能够方便地应用于核领域不同的信号源和不同的探测器。
参考文献
[1] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Radeka V. Trapezoidal filtering of signals from large germanium detectors at high rates[J]. IEEE Transactions on Nuclear Science, 1972, 99(3): 525–539. |
[2] | Journal of Central South University(Science and Technology), 41(2):649-654.-->肖无云, 魏义祥, 艾宪芸. 数字化多道脉冲幅度分析中的梯形成形算法[J]. 清华大学学报:自然科学版, 2005, 45(6): 810–812.XIAO Wuyun, WEI Yixiang, AI Xianyun. Trapezoidal shaping algorithm for digital multi-channel pulse height analysis[J]. Journal of Tsinghua University: Science and Technology, 2005, 45(6): 810–812. (in Chinese) |
[3] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Bolic M, Drndarevic V. Digital gamma-ray spectroscopy based on FPGA technology[J]. Nuclear Instruments & Methods in Physics Research, 2002, 482(3): 761–766. |
[4] | Journal of Central South University(Science and Technology), 41(2):649-654.--> Bolic M, Drndarevic V. Processing architecture for high count rate spectrometry with NaI(Tl) detector [C]//IEEE Instrumentation and Measurement Technology Conference Record. Victoria, BC, Canada: IEEE Press, 2008: 274-278. |
[5] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Cano-Ott D, Tain J L, Gadea A, et al. Pulse pileup correction of large NaI(Tl) total absorption spectra using the true pulse shape[J]. Nuclear Instruments & Methods in Physics Research, 1999, 430(2): 488–497. |
[6] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Wong W H, Li H. A scintillation detector signal processing technique with active pileup prevention for extending scintillation count rates[J]. IEEE Transactions on Nuclear Science, 1998, 45(3): 838–842. DOI:10.1109/23.682647 |
[7] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Haselman M D, Pasko J, Hauck S, et al. FPGA-based pulse pile-up correction with energy and timing recovery[J]. IEEE Transactions on Nuclear Science, 2012, 59(5): 1823–1830. DOI:10.1109/TNS.2012.2207403 |
[8] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Ehrenberg J E, Ewart T E, Morris R D. Signal-processing techniques for resolving individual pulses in a multipath signal[J]. Journal of the Acoustical Society of America, 1978, 63(6): 1861–1865. DOI:10.1121/1.381926 |
[9] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Bolic M, Drndarevic V, Gueaieb W. Pileup correction algorithms for very-high-count-rate gamma-ray spectrometry with NaI(Tl) detectors[J]. IEEE Transactions on Instrumentation & Measurement, 2010, 59(1): 122–130. |
[10] | Journal of Central South University(Science and Technology), 41(2):649-654.-->LIU Zhenzhou, CHEN Jinxiang. A Monte Carlo based technique for analysing gamma-ray spectra[J]. Measurement Science & Technology, 2008, 19(8): 085102. |
[11] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Meng L J, Ramsden D. An inter-comparison of three spectral-deconvolution algorithms for gamma-ray spectroscopy[J]. IEEE Transactions on Nuclear Science, 2000, 47(4): 1329–1336. DOI:10.1109/23.872973 |
[12] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Imperiale C, Imperiale A. On nuclear spectrometry pulses digital shaping and processing[J]. Measurement, 2001, 30(1): 49–73. DOI:10.1016/S0263-2241(00)00057-9 |
[13] | Journal of Central South University(Science and Technology), 41(2):649-654.-->陈新光, 李翔宇, 孙义和. 粒子探测器读出电路数字滤波器设计[J]. 电子产品世界, 2010, 17(10): 29–31.CHEN Xinguang, LI Xiangyu, SUN Yihe. Particle detector readout circuits digital filter design[J]. Electronic Engineering & Product World, 2010, 17(10): 29–31. DOI:10.3969/j.issn.1005-5517.2010.09.006(in Chinese) |