CN111538067B - 一种数字化核脉冲直线成形方法 - Google Patents
一种数字化核脉冲直线成形方法 Download PDFInfo
- Publication number
- CN111538067B CN111538067B CN202010372393.4A CN202010372393A CN111538067B CN 111538067 B CN111538067 B CN 111538067B CN 202010372393 A CN202010372393 A CN 202010372393A CN 111538067 B CN111538067 B CN 111538067B
- Authority
- CN
- China
- Prior art keywords
- digital
- nuclear pulse
- nuclear
- pulse
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/36—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T7/00—Details of radiation-measuring instruments
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明公开了一种数字化核脉冲直线成形方法,其特征在于,构建出数字化核脉冲的直线成形方法传递公式;同时利用传递公式内的分量实现重叠核脉冲的堆积判弃;与基线提取结果相减获得数字化核脉冲的准确幅度值;最后给出了传递公式中各个参数的获取方法。与梯形成形相比,该方法算法更简单、无需快成形通道可同时实现核脉冲的定位、只要核脉冲的宽度大于可分辨间隔即可有效准确提取其幅度值;同时成形后核脉冲平顶区更平滑、无需下降沿,这些优势使得本发明算法能在保证能量分辨率的同时提高了能谱计数率。
Description
技术领域
本发明涉及辐射环境测量领域,应用于数字化核能谱测量系统,具体是指一种数字化核脉冲直线成形方法。
背景技术
辐射射线与探测器发生相互作用后被探测器外加电场收集形成的电信号,被前置放大器转变成幅度与入射射线的能量成正比的核脉冲,信号调理电路将核脉冲放大并调节参考电压以匹配高速ADC(模数转换器)采样需求,经高速ADC等间隔采样获得数字化核脉冲,将其存储于FIFO(先进先出)等存储器中等待FPGA(现场可编程门阵列)等可编程器件处理。在FPGA内需实现基线提取、滤波成形、幅度提取、基线恢复、堆积判弃等过程以获得各个数字化核脉冲的幅度值并在相应道记录单元中加1,累积一段时间后按道址大小将相应道记录单元内数值排列显示即为能谱。
被滤波成形方法成形后的核脉冲宽度对能谱的计数率影响很大。目前滤波成形中多采用高斯成形方法(数字化Sallen–Key成形和数字化CR-RCm成形)和梯形成形方法,相比较而言梯形成形方法核脉冲宽度最短、重叠核脉冲的分辨能力最强。
本发明提出一种数字化核脉冲直线成形方法,与梯形成形相比,既能保证能谱的分辨率又能更好地分辨重叠核脉冲。
迄今为止,尚未见到数字化核脉冲直线成形方法的报道。
发明内容
本发明的目的在于克服上述背景技术的不足之处,提出一种数字化核脉冲直线成形方法,核心在于保证能谱的分辨率的同时提高重叠核脉冲的分辨能力。
本发明的目的是通过如下技术方案实现的,一种数字化核脉冲直线成形方法,构建出数字化核脉冲的直线成形方法传递公式;同时利用传递公式内的分量实现重叠核脉冲的堆积判弃;与基线提取结果相减获得数字化核脉冲的准确幅度值;最后给出了传递公式中各个参数的获取方法。
进一步,所述构建出数字化核脉冲的直线成形方法传递公式为:
(1)核脉冲数学模型:
当τ2趋近于0时,yn=fi,即直线成形方法;
(2)离散之前先对xn进行三点平滑,所以离散化后传递公式为:
yn=A·{yn-1+[xn+1-xn-2]/3+[xn+1+xn+xn-1]/3M}。
更进一步,(1)、先确定数字化能谱仪内ADC的采样频率后,采集系列非重叠核脉冲和参考电压离散值,用于确定参数值,包括Nos、M、S和A的值;
(2)、构建数字化核脉冲直线成形方法传递公式:
yn=A·{yn-1+[xn+1-xn-2]/3+[xn+1+xn+xn-1]/3M};
(3)将参数Nos、M、S和A代入上述构建后数字化核脉冲直线成形方法传递公式,并于可编程器件内按如下步骤实现:
①依公式gn=[xn+1-xn-2]/3计算gn并存储在硬件缓冲单元内;
②如果gn+1≥Nos且gn<Nos说明第n点处于新i个核脉冲的起始点Ui,依公式yn=A·{yn-1+gn+[xn+1+xn+xn-1]/3M}计算数字化核脉冲直线成形后的值;否则依公式yn=A·{gn+[xn+1+xn+xn-1]/3M}计算数字化核脉冲直线成形后的值;同时将yn存储在硬件缓冲单元内;
(4)参数Nos、M、S和A的获取方法,是由以下方法得到的:
①参数Nos的确定:通过所用数字化能谱仪采集一段不存在数字化核脉冲的参考电压离散值(不少于五千个点),将其最大值赋于Nos;
②参数M的确定:通过所用数字化能谱仪采集1个非重叠数字化核脉冲,其中非重叠数字化核脉冲幅度采样值大于能谱最大道址的一半;采用公式对上述数字化核脉冲进行拟合,取M≤τ1均可实现数字化核脉冲直线成形,取M=τ1时重叠核脉冲甄别能力最优,能谱计数率最高;
③参数S的确定:通过所用数字化能谱仪采集1个非重叠数字化核脉冲,其中非重叠数字化核脉冲幅度采样值大于能谱最大道址的一半;在A=1时按前述成形方案对上述数字化核脉冲实施直线成形,用公式yn=a-b·e-c·n拟合当前核脉冲的起始点后的直线成形值yn,则取
④参数A的确定:通过所用数字化能谱仪采集12个不同幅度的非重叠数字化核脉冲,其中非重叠数字化核脉冲幅度采样值从20到能谱最大道址;在A=1时按前述方案提取上述12个数字化核脉冲的幅度值Fi(i=1,……,12),并用公式拟合上述12个数字化核脉冲得到幅度拟合值fi(i=1,……,12);运用公式fi=A·Fi拟合12个数据点(Fi,fi)得到参数A的值。
本发明有益效果及创造性分析:
从图1可以看出,经本发明成形后核脉冲中上升沿和平顶与梯形成形结果相似。与梯形成形相比,本发明还存在以下优点:
①成形后核脉冲平顶区更平滑(如图1),证实本发明去除电子学噪声的效果更佳,能量分辨力亦越好;
②无需下降沿(如图2),证实本发明成形后核脉冲可以更窄,重叠核脉冲分辨能力更强。即使重叠核脉冲的起始点间隔小于S+5导致第一个核脉冲幅度无法准确提取(如图2内第4个核脉冲),该方法也不会影响后续核脉冲幅度的准确提取(如图2内第5个核脉冲);
③算法更简单,无需快成形通道可同时实现核脉冲的定位,说明在相同的测量时间内本发明可以处理更多的核脉冲,能谱计数率更高。
附图说明
图1为本发明非重叠数字化核脉冲采样值的数字化直线成形与梯形成形结果;
图2为本发明25.6微秒核脉冲采样值的数字化直线成形与梯形成形结果;
图3为本发明参考电压ADC长时间采样值;
图4为本发明非重叠数字化核脉冲采样值及其拟合分布;
图5为本发明非重叠数字化核脉冲直线成形结果及其拟合分布;
图6为本发明参数A的拟合结果;
图7为本发明实施例1的FPGA内处理流程图;
图8为本发明于硬件中实现数字化直线成形结果后采集得到的137Cs点源的能谱。
具体实施方式
下面结合附图和具体实施例对本发明做进一步的描述,然而本发明的范围并不限于下述实施例。
实施例1:下面通过附图和具体实施方式对本发明作更为详细的描述。
实施例的数字化能谱仪硬件由滨松公司CH249-01型NaI探头、信号调理电路、AD9649、EP3C10E144C8N型FPGA芯片、CY7C68015A型单片机及相应的辅助电路组成。AD9649的采样率设置为40MSPS。
(1)依据核脉冲指数衰减规律,利用衰减补偿原理,建立数字化核脉冲直线成形方法离散传递公式:
①依公式gn=[xn+1-xn-2]/3计算gn并存储在硬件缓冲单元内;
②如果gn+1≥Nos且gn<Nos说明第n点处于新i个核脉冲的起始点Ui,依公式yn=A·{yn-1+gn+[xn+1+xn+xn-1]/3M}计算数字化核脉冲直线成形后的值;否则依公式yn=A·{gn+[xn+1+xn+xn-1]/3M}计算数字化核脉冲直线成形后的值。同时将yn存储在硬件缓冲单元内。
(2)参数的获取
①参数Nos的确定:采用本实施例中数字化能谱仪采集一段不存在数字化核脉冲的参考电压离散值,如图3所示。可以看出其最大值为6,则取Nos=6。
②参数M的确定:采用本实施例中数字化能谱仪采集到的1个幅度较高(幅度采样值大于能谱最大道址的一半)非重叠数字化核脉冲如图4所示。于MATLAB平台采用公式对图4种数字化核脉冲进行拟合(拟合结果如图4所示),得到拟合值τ1=16.87,为实现计数率最大化,取M=τ1=16.87。
④参数A的确定:采用本实施例中数字化能谱仪采集12个不同幅度(幅度采样值从20到能谱最大道址)的非重叠数字化核脉冲,在A=1时按步骤(1)提取上述12个数字化核脉冲的幅度值Fi(i=1,……,12),并用公式拟合上述12个数字化核脉冲得到幅度拟合值fi(i=1,……,12)。运用公式fi=A·Fi拟合12个数据点(Fi,fi),结果如图6所示。得到参数A=0.5176。
(3)将参数Nos、M、S和A代入步骤(1)中②和③即可获得适用于本实施例中数字化能谱仪的数字化直线成形方法的传递公式。于MATLAB平台实现步骤(1),对本实施例中数字化能谱仪获得的25.6微秒连续采样值进行处理,得到数字化核脉冲直线成形方法的输出结果如图2所示。可以看出:在Ui+1-Ui≥41时,数字化直线成形方法能准确提取每个幅度值;即使第4个核脉冲不满足上述条件导致无法提取核脉冲幅度值,也不会对第5和6个核脉冲幅度准确提取造成影响,证实可有效提升能谱仪的计数率。
在本实施例中数字化能谱仪内FPGA芯片上实施数字化核脉冲直线成形方法,详细过程如图7所示,其中基线Bi利用最小平均值法提取;gn由步骤(1)中①确定;yn由步骤(1)中②确定;幅度提取、堆积判弃和基线恢复由步骤(1)中③确定。对每一个可提取幅度的核脉冲而言,能谱记录即将第道计数加1。获得的137Cs点源的能谱如图8所示,分析发现其能量分辨力为7.58%,略高于梯形成形方法的7.61%;计数率较梯形成形提高了4.4%。
本发明所涉及的所有参数符号的具体含义为:
n:数字化能谱仪内ADC连续采样值的序号;
xn:数字化能谱仪内ADC的第n次采样值;
yn:数字化直线成形方法的第n个成形结果;
gn:数字化直线成形方法内分量,可用于确定核脉冲的起始位置Ui;
i:第i个核脉冲;
Ui:第i个核脉冲的起始位置;
A:经数字化直线成形方法成形后核脉冲幅度校正系数;
M:数字化直线成形方法内参数,取值等于核脉冲下降沿衰减系数;
S:经数字化直线成形方法成形后核脉冲上升沿持续时间;
Nos:数字化能谱仪参考电压的ADC采样最大值;
Fi:经数字化直线成形方法成形后第i个核脉冲的幅度值;
fi:采用双指数函数拟合得到的第i个核脉冲的幅度值;
Bi:利用最小平均值法提取的第i个核脉冲的基线值;
τ1:数字化核脉冲下降沿衰减时间常数;
τ2:数字化核脉冲上升沿衰减时间常数;
a:经数字化直线成形方法成形后核脉冲变化规律拟合函数内参数;
b:经数字化直线成形方法成形后核脉冲变化规律拟合函数内参数;
c:经数字化直线成形方法成形后核脉冲变化规律拟合函数内参数。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本实用新型的保护范围之内。
Claims (4)
1.一种数字化核脉冲直线成形方法,其特征在于,构建出数字化核脉冲的直线成形方法传递公式;同时利用传递公式内的分量实现重叠核脉冲的堆积判弃;与基线提取结果相减获得数字化核脉冲的准确幅度值;最后给出了传递公式中各个参数的获取方法;
(1)、先确定数字化能谱仪内ADC的采样频率后,采集系列非重叠核脉冲和参考电压离散值,用于确定参数值,包括Nos、M、S和A的值;
(2)、构建数字化核脉冲直线成形方法传递公式:
yn=A·{yn-1+[xn+1-xn-2]/3+[xn+1+xn+xn-1]/3M};
(3)将参数Nos、M、S和A代入上述构建后数字化核脉冲直线成形方法传递公式,并于可编程器件内按如下步骤实现:
①依公式gn=[xn+1-xn-2]/3计算gn并存储在硬件缓冲单元内;
②如果gn+1≥Nos且gn<Nos说明第n点处于新i个核脉冲的起始点Ui,依公式yn=A·{yn-1+gn+[xn+1+xn+xn-1]/3M}计算数字化核脉冲直线成形后的值;否则依公式yn=A·{gn+[xn+1+xn+xn-1]/3M}计算数字化核脉冲直线成形后的值;同时将yn存储在硬件缓冲单元内;
(4)参数Nos、M、S和A的获取方法,是由以下方法得到的:
①参数Nos的确定:通过所用数字化能谱仪采集一段不存在数字化核脉冲的参考电压离散值,不少于五千个点,将其最大值赋于Nos;
②参数M的确定:通过所用数字化能谱仪采集1个非重叠数字化核脉冲,其中非重叠数字化核脉冲幅度采样值大于能谱最大道址的一半;采用公式对上述数字化核脉冲进行拟合,取M≤τ1均可实现数字化核脉冲直线成形,取M=τ1时重叠核脉冲甄别能力最优,能谱计数率最高;
③参数S的确定:通过所用数字化能谱仪采集1个非重叠数字化核脉冲,其中非重叠数字化核脉冲幅度采样值大于能谱最大道址的一半;在A=1时按前述成形方案对上述数字化核脉冲实施直线成形,用公式yn=a-b·e-c·n拟合当前核脉冲的起始点后的直线成形值yn,则取
3.根据权利要求2所述的一种数字化核脉冲直线成形方法,其特征在于,当Ui+1-UipS+5时,丢弃第i个核脉冲;否则第i个核脉冲的幅度提取值为Fi=(yS+1+yS+2+yS+3+yS+4+yS+5)/5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010372393.4A CN111538067B (zh) | 2020-05-06 | 2020-05-06 | 一种数字化核脉冲直线成形方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010372393.4A CN111538067B (zh) | 2020-05-06 | 2020-05-06 | 一种数字化核脉冲直线成形方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111538067A CN111538067A (zh) | 2020-08-14 |
CN111538067B true CN111538067B (zh) | 2022-09-06 |
Family
ID=71968675
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010372393.4A Active CN111538067B (zh) | 2020-05-06 | 2020-05-06 | 一种数字化核脉冲直线成形方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111538067B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112134545B (zh) * | 2020-09-25 | 2022-04-22 | 中国核动力研究设计院 | 一种基于最佳滤波器的梯形成形方法、系统、终端及介质 |
CN112764082B (zh) * | 2020-12-08 | 2023-05-23 | 武汉第二船舶设计研究所(中国船舶重工集团公司第七一九研究所) | 一种基于fpga的核脉冲数字化成形采样方法 |
CN113835114B (zh) * | 2021-08-25 | 2024-04-26 | 吉林大学 | 一种紧凑式高能伽玛射线反符合叠层探测器 |
CN115980821B (zh) * | 2023-01-16 | 2023-09-05 | 中国科学院近代物理研究所 | 一种基于fpga的自适应能谱测量方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5067090A (en) * | 1989-05-01 | 1991-11-19 | Schlumberger Technology Corporation | Nuclear spectroscopy method and apparatus for digital pulse height analysis |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1003157A (en) * | 1962-03-29 | 1965-09-02 | Decca Ltd | Improvements in or relating to radio guidance systems for aircraft |
JPS6045833B2 (ja) * | 1978-09-26 | 1985-10-12 | 株式会社東芝 | 放射線測定装置 |
US4687622A (en) * | 1985-10-29 | 1987-08-18 | Irt Corporation | Nuclear event detector |
US6609075B1 (en) * | 2001-06-04 | 2003-08-19 | William K. Warburton | Method and apparatus for baseline correction in x-ray and nuclear spectroscopy systems |
CN100396242C (zh) * | 2003-01-20 | 2008-06-25 | 株式会社东芝 | 具有x-射线荧光检查功能的x-射线诊断设备 |
IL225219A (en) * | 2013-03-14 | 2017-08-31 | Micro Tag Temed Ltd | A magnetic resonance facility for the detection and detection of tiny quantities of specific materials quickly and efficiently |
CN103941279A (zh) * | 2014-05-08 | 2014-07-23 | 东华理工大学 | 基于脉宽和幅度双参数的核脉冲信号数字高斯成形方法 |
KR101621366B1 (ko) * | 2015-03-30 | 2016-05-17 | 한양대학교 산학협력단 | 테라헤르츠 전달함수 모델을 이용한 비전도 매질의 테라헤르츠 영역에서의 유전상수 및 감쇠비 계측 장치 및 방법 |
CN105212954B (zh) * | 2015-11-05 | 2018-03-16 | 苏州瑞派宁科技有限公司 | 一种脉冲堆积事件实时处理方法与系统 |
CN106772545B (zh) * | 2017-02-09 | 2019-04-16 | 成都理工大学 | 一种采用尖顶成形算法的数字脉冲幅度分析器 |
CN107193036B (zh) * | 2017-06-26 | 2018-12-18 | 成都理工大学 | 一种改进型核信号梯形脉冲成形方法和装置 |
CN110941007B (zh) * | 2018-09-25 | 2022-02-08 | 北京超分科技有限公司 | 核能谱的处理方法和设备 |
-
2020
- 2020-05-06 CN CN202010372393.4A patent/CN111538067B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5067090A (en) * | 1989-05-01 | 1991-11-19 | Schlumberger Technology Corporation | Nuclear spectroscopy method and apparatus for digital pulse height analysis |
Also Published As
Publication number | Publication date |
---|---|
CN111538067A (zh) | 2020-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111538067B (zh) | 一种数字化核脉冲直线成形方法 | |
WO2014121548A1 (zh) | 数字化闪烁脉冲的基线校正方法及系统 | |
US10175369B2 (en) | Signal data processing for radiation detection | |
JP2008501954A5 (zh) | ||
WO2016110141A1 (zh) | 一种闪烁脉冲的数字化方法 | |
JP2008501954A (ja) | X線及び核分光システムにおける検出限界を向上するための方法及び装置 | |
CN109557573B (zh) | 一种数字化的α/β射线甄别方法 | |
US9767258B2 (en) | System and method of extending the linear dynamic range of event counting | |
EP3073292B1 (en) | Method and system for recovering scintillation pulse information | |
US10795031B2 (en) | Radiation detection system and signal processor for radiation detection | |
CN111380957B (zh) | 一种基于频率-波数域滤波的超声导波检测方法 | |
CN110308476A (zh) | 一种粒子辐射探测方法及探测装置 | |
CN112764082A (zh) | 一种基于fpga的核脉冲数字化成形采样方法 | |
Tang et al. | Counting-loss correction for X-ray spectra using the pulse-repairing method | |
JP7079748B2 (ja) | 放射線検出装置及び放射線検出用信号処理装置 | |
Wu et al. | Implementation of a straight-line shaper on an FPGA for a UAVAGS | |
Paul et al. | Implementation of FPGA based real-time digital DAQ for high resolution, and high count rate nuclear spectroscopy application | |
JP6999593B2 (ja) | 検出装置及びスペクトル生成装置 | |
CN114114375B (zh) | 一种等离子热电子测量系统及测量方法 | |
Garcia-Belmonte et al. | Digital implementation of filters for nuclear applications using the discrete wavelet transform | |
Richards et al. | Real-time pulse shape discrimination for neutron scattering science using EJ-270 | |
Uslenghi et al. | A flexible readout electronics for CZT arrays with real time digital signal processing | |
Bjerge et al. | A system on programmable chip design of a digitizer with improved trapezoidal filter validation | |
JP2003215251A (ja) | 陽電子寿命測定方法及び装置 | |
Peng et al. | A Digital Method of Noise Measurement of Charge-sensitive Preamplifiers Used with Detectors of Ionizing Radiation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |