CN110568468B - 一种辐射脉冲计数突变算法 - Google Patents

一种辐射脉冲计数突变算法 Download PDF

Info

Publication number
CN110568468B
CN110568468B CN201910797059.0A CN201910797059A CN110568468B CN 110568468 B CN110568468 B CN 110568468B CN 201910797059 A CN201910797059 A CN 201910797059A CN 110568468 B CN110568468 B CN 110568468B
Authority
CN
China
Prior art keywords
pulses
counting
pulse
mutation
data
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
Application number
CN201910797059.0A
Other languages
English (en)
Other versions
CN110568468A (zh
Inventor
丁炳亮
张华�
魏鸿林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fuzhou Gem Instrument Equipment Co ltd
Original Assignee
Fuzhou Gem Instrument Equipment Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fuzhou Gem Instrument Equipment Co ltd filed Critical Fuzhou Gem Instrument Equipment Co ltd
Priority to CN201910797059.0A priority Critical patent/CN110568468B/zh
Publication of CN110568468A publication Critical patent/CN110568468A/zh
Application granted granted Critical
Publication of CN110568468B publication Critical patent/CN110568468B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/15—Instruments in which pulses generated by a radiation detector are integrated, e.g. by a diode pump circuit

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Measurement Of Radiation (AREA)

Abstract

本发明提供了脉冲计数技术领域的一种辐射脉冲计数突变算法，包括如下步骤：步骤S1、设定一动态阈值σ以及一触发数值n；步骤S2、获取脉冲个数m以及各脉冲的计数数据；步骤S3、m≥n？，若是，进入步骤S4；若否，结束流程；步骤S40、设定一初始计算数k以及一当前计算数i；步骤S5、i＜m/2？，若是，进入步骤S6；若否，结束流程；步骤S6、计算从0至m‑i个脉冲的计数数据的平均值A1，计算从m‑i至m个脉冲的计数数据的平均值A2，|A1‑A2|＞σ？若是，进入步骤S7；若否，则i+1，进入步骤S5；步骤S7、记录辐射脉冲计数突变1次，此时的i对应辐射脉冲计数突变位置。本发明的优点在于：实现提升脉冲型辐射检测仪器对脉冲突变计数的检测效果，提高响应速度，且不增加硬件成本。

Description

一种辐射脉冲计数突变算法
技术领域
本发明涉及脉冲计数技术领域，特别指一种辐射脉冲计数突变算法。
背景技术
脉冲型辐射检测仪器的核心是对脉冲计数的处理，传统的脉冲型辐射检测仪器大多使用低灵敏度脉冲型辐射传感器，并采用固定的单阈值对脉冲计数进行判断。由于脉冲型辐射传感器输出的脉冲时间分布具有随机性，使得在相同间隔的统计时间内脉冲计数有涨落，但脉冲计数服从泊松分布，可以求得脉冲计数的均方根
Figure BDA0002181260480000011
进而求得相对统计误差
Figure BDA0002181260480000012
其中N表示测量时间内的记脉冲计数；为了减小相对统计误差需要增大N，传统上采用提高脉冲型辐射传感器的灵敏度以及延长测量时间的方法。
但是，传统的方法存在有如下缺点：1、采用固定的单阈值对脉冲计数进行判断，在低辐射强度下检测计数突变效果差，且在超过阈值的辐射强度后，若辐射强度继续变高就无法再次检测到计数的突变，即使采用多个固定的单阈值算法，本质还是固定阈值，检测效果没有得到本质上的改善；2、灵敏度高的脉冲型辐射传感器价格昂贵；3、延长测量时间会牺牲脉冲型辐射检测仪器的响应速度。
因此，如何提供一种辐射脉冲计数突变算法，实现提升脉冲型辐射检测仪器对脉冲突变计数的检测效果，提高响应速度，且不增加硬件成本，成为一个亟待解决的问题。
发明内容
本发明要解决的技术问题，在于提供一种辐射脉冲计数突变算法，实现提升脉冲型辐射检测仪器对脉冲突变计数的检测效果，提高响应速度，且不增加硬件成本。
本发明是这样实现的：一种辐射脉冲计数突变算法，包括如下步骤：
步骤S10、设定一动态阈值σ以及一触发数值n；
步骤S20、获取脉冲个数m以及各脉冲的计数数据；
步骤S30、判断脉冲个数m是否大于等于所述触发数值n，若是，则进入步骤S40；若否，则结束流程；
步骤S40、设定一初始计算数k以及一当前计算数i，且当前计算数i从初始计算数k开始进行计算；其中n、m、k、i均为正整数，且n≤m，k≤i；
步骤S50、判断当前计算数i是否小于脉冲个数m的一半，若是，则进入步骤S60；若否，则结束流程；
步骤S60、计算从0至m-i个脉冲的计数数据的平均值A1，计算从m-i至m个脉冲的计数数据的平均值A2，判断A2与A1差值的绝对值是否大于所述动态阈值σ，若是，则进入步骤S70；若否，则将当前计算数i加1，并进入步骤S50；
步骤S70、记录辐射脉冲计数突变1次，此时的i对应辐射脉冲计数突变位置；删除从0至m-i个脉冲的计数数据，保留从m-i至m个脉冲的计数数据。
进一步地，所述步骤S10中，
所述动态阈值
Figure BDA0002181260480000021
其中N1表示0至m-i个脉冲的计数数据之和，N2表示从m-i至m个脉冲的计数数据之和，
Figure BDA0002181260480000022
表示用m-i至m个脉冲的计数数据之和计算动态阈值σ，
Figure BDA0002181260480000023
为使动态阈值σ在低本底下有更好表现的系数。
进一步地，所述步骤S60中，所述A2与A1差值的绝对值具体为：
Figure BDA0002181260480000031
本发明的优点在于：
1、通过延长测量时间获得更多的脉冲个数来减小相对统计误差，同时每秒对各脉冲计数数据的涨落的误差进行累积(对时间积分)，判断某一时长内误差累积值是否大于所述动态阈值σ，即判断A2与A1差值的绝对值是否大于所述动态阈值σ，极大的提升了脉冲型辐射检测仪器对脉冲突变计数的检测效果。
2、通过设置动态阈值σ，实现更快跟踪各脉冲的计数数据的变化量，极大的提高了脉冲型辐射检测仪器检测的响应速度。
3、通过延长测量时间的方式，利用灵敏度低的脉冲型辐射传感器也能实现很好的检测效果，即不增加硬件成本。
附图说明
下面参照附图结合实施例对本发明作进一步的说明。
图1是本发明一种辐射脉冲计数突变算法的流程图。
图2是本发明动态阈值系数项的仿真图。
具体实施方式
请参照图1至图2所示，本发明一种辐射脉冲计数突变算法的较佳实施例，包括如下步骤：
步骤S10、设定一动态阈值σ以及一触发数值n；通过设置动态阈值σ，实现更快跟踪各脉冲的计数数据的变化量，极大的提高了脉冲型辐射检测仪器检测的响应速度；通过设置触发数值n，当脉冲个数m大于触发数值n时才启用突变算法，有效预防由于低灵敏度的GM管在本底下的计数值太小，而意外检测到数据突变的情况；
步骤S20、获取脉冲个数m以及各脉冲的计数数据；各脉冲的计数数据为一个数组，数组长度根据需要的相对统计误差而自行设定；
步骤S30、判断脉冲个数m是否大于等于所述触发数值n，若是，则进入步骤S40；若否，则结束流程；
步骤S40、设定一初始计算数k以及一当前计算数i，且当前计算数i从初始计算数k开始进行计算；其中n、m、k、i均为正整数，且n≤m，k≤i；如果传感器灵敏度不大且响应时间要求不严格，初始计算数k建议取大于1的值；
步骤S50、判断当前计算数i是否小于脉冲个数m的一半，若是，则进入步骤S60；若否，则结束流程；
步骤S60、计算从0至m-i个脉冲的计数数据的平均值A1，计算从m-i至m个脉冲的计数数据的平均值A2，判断A2与A1差值的绝对值是否大于所述动态阈值σ，若是，则进入步骤S70；若否，则将当前计算数i加1，并进入步骤S50；
步骤S70、记录辐射脉冲计数突变1次，此时的i对应辐射脉冲计数突变位置；删除从0至m-i个脉冲的计数数据，保留从m-i至m个脉冲的计数数据。保留从m-i至m个脉冲的计数数据，即保留突变位置后的数据是为了更好的稳定数据。
所述步骤S10中，
所述动态阈值
Figure BDA0002181260480000041
其中N1表示0至m-i个脉冲的计数数据之和，N2表示从m-i至m个脉冲的计数数据之和，
Figure BDA0002181260480000042
表示用m-i至m个脉冲的计数数据之和计算动态阈值σ，
Figure BDA0002181260480000043
为使动态阈值σ在低本底下有更好表现的系数。
所述步骤S60中，所述A2与A1差值的绝对值具体为：
Figure BDA0002181260480000044
实际求取A2与A1差值的绝对值时，为了减少低本底下经常检测到突变的情况，需要缩小i倍：
Figure BDA0002181260480000051
A2与A1差值的绝对值大于所述动态阈值σ的表达式如下：
Figure BDA0002181260480000052
进而推导出：
Figure BDA0002181260480000053
即动态阈值包括一系数项
Figure BDA0002181260480000054
当m＝60，i＝0…29时的仿真图参见图2，从图2可知：
随着i的增大，动态阈值的系数项也会随着增大，这样可以避免本底下经常超过阈值的情况，当i＞22时动态阈值的系数项开始下降，这是考虑到i增大时，m-i至m个脉冲的计数数据之和也增大，计数数据的涨落将减小，即动态阈值的系数项也应对应减小。
本发明在计算动态阈值时需要进行频繁的开方运算，对于低配的MCU是非常占用CPU资源的，采用如下开方算法可以提高计算效率：
Figure BDA0002181260480000055
本发明的试验验证：
实验采用型号为J304盖革计数管(能量补偿型)，灵敏度为1uSv≥2CPS，测量统计时间为60S，每秒读取一次计数值，使用Am241作为测试辐射源。靠近Am241约2-3S立即检测到计数值突变，剂量值显示2uSv/h，响应效果非常好。放在辐射源上5分钟没有出现意外检测到计数突变，更长时间(半小时以上)的持续观察偶尔会出现意外检测到计数突变，说明数据稳定性不错，和市场常用的同灵敏度辐射检测仪比较在计数突变和数据稳定性方面有了非常大的改善。
综上所述，本发明的优点在于：
1、通过延长测量时间获得更多的脉冲个数来减小相对统计误差，同时每秒对各脉冲计数数据的涨落的误差进行累积(对时间积分)，判断某一时长内误差累积值是否大于所述动态阈值σ，即判断A2与A1差值的绝对值是否大于所述动态阈值σ，极大的提升了脉冲型辐射检测仪器对脉冲突变计数的检测效果。
2、通过设置动态阈值σ，实现更快跟踪各脉冲的计数数据的变化量，极大的提高了脉冲型辐射检测仪器检测的响应速度。
3、通过延长测量时间的方式，利用灵敏度低的脉冲型辐射传感器也能实现很好的检测效果，即不增加硬件成本。
虽然以上描述了本发明的具体实施方式，但是熟悉本技术领域的技术人员应当理解，我们所描述的具体的实施例只是说明性的，而不是用于对本发明的范围的限定，熟悉本领域的技术人员在依照本发明的精神所作的等效的修饰以及变化，都应当涵盖在本发明的权利要求所保护的范围内。

Claims (1)

1.一种辐射脉冲计数突变算法，其特征在于：包括如下步骤：
步骤S10、设定一动态阈值σ以及一触发数值n；
步骤S20、获取脉冲个数m以及各脉冲的计数数据；
步骤S30、判断脉冲个数m是否大于等于所述触发数值n，若是，则进入步骤S40；若否，则结束流程；
步骤S40、设定一初始计算数k以及一当前计算数i，且当前计算数i从初始计算数k开始进行计算；其中n、m、k、i均为正整数，且n≤m，k≤i；
步骤S50、判断当前计算数i是否小于脉冲个数m的一半，若是，则进入步骤S60；若否，则结束流程；
步骤S60、计算从0至m-i个脉冲的计数数据的平均值A1，计算从m-i至m个脉冲的计数数据的平均值A2，判断A2与A1差值的绝对值是否大于所述动态阈值σ，若是，则进入步骤S70；若否，则将当前计算数i加1，并进入步骤S50；
步骤S70、记录辐射脉冲计数突变1次，此时的i对应辐射脉冲计数突变位置；删除从0至m-i个脉冲的计数数据，保留从m-i至m个脉冲的计数数据；
所述步骤S10中，
所述动态阈值
Figure FDA0003946167120000011
其中N1表示0至m-i个脉冲的计数数据之和，N2表示从m-i至m个脉冲的计数数据之和，
Figure FDA0003946167120000012
表示用m-i至m个脉冲的计数数据之和计算动态阈值σ，
Figure FDA0003946167120000013
为使动态阈值σ在低本底下有更好表现的系数；
所述步骤S60中，所述A2与A1差值的绝对值具体为：
Figure FDA0003946167120000021
CN201910797059.0A 2019-08-27 2019-08-27 一种辐射脉冲计数突变算法 Active CN110568468B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910797059.0A CN110568468B (zh) 2019-08-27 2019-08-27 一种辐射脉冲计数突变算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910797059.0A CN110568468B (zh) 2019-08-27 2019-08-27 一种辐射脉冲计数突变算法

Publications (2)

Publication Number Publication Date
CN110568468A CN110568468A (zh) 2019-12-13
CN110568468B true CN110568468B (zh) 2023-01-10

Family

ID=68776296

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910797059.0A Active CN110568468B (zh) 2019-08-27 2019-08-27 一种辐射脉冲计数突变算法

Country Status (1)

Country Link
CN (1) CN110568468B (zh)

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3586982A (en) * 1969-02-19 1971-06-22 Westinghouse Electric Corp Pulse count detector having fail-safe features
US3696399A (en) * 1970-09-11 1972-10-03 Coulter Electronics Range expansion method and apparatus for multichannel pulse analysis
US4767929A (en) * 1986-10-06 1988-08-30 The United States Of America As Represented By The United State Department Of Energy Extended range radiation dose-rate monitor
US5258926A (en) * 1988-08-08 1993-11-02 Osterreichesches Forschungszentrum Seibersdorf Gmbh Method of measuring radiation for a radiation measuring device
US5377282A (en) * 1991-09-19 1994-12-27 International Business Machines Corporation Optical inspection system utilizing dynamic analog-to-digital thresholding
US5640109A (en) * 1995-10-27 1997-06-17 Mts Systems Corporation Pulse detector
US5928150A (en) * 1997-10-04 1999-07-27 Neoprobe Corporation System for locating and detecting a source of photon emissions
JP2004333326A (ja) * 2003-05-08 2004-11-25 Nippon Kankyo Chosa Kenkyusho:Kk 放射線測定器および放射線測定方法
EP1906211A1 (en) * 2005-07-22 2008-04-02 ICX Radiation GmbH Detector for the measurement of Ionizing radiation
CN102338880A (zh) * 2011-06-02 2012-02-01 中国科学技术大学 核脉冲幅度数字化方法及系统
CN102859391A (zh) * 2011-04-21 2013-01-02 株式会社东芝 电子设备以及输出方法
JP2014085183A (ja) * 2012-10-22 2014-05-12 Panasonic Corp 放射線測定装置、放射線測定方法及び放射線測定プログラム
CN108886594A (zh) * 2016-04-12 2018-11-23 雷神公司 具有模拟计数器元件的数字单位单元

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8178848B2 (en) * 2009-12-21 2012-05-15 General Electric Company Systems and methods for filtering noise in pixelated image detectors
EP4342521A3 (en) * 2017-07-11 2024-05-08 RefleXion Medical Inc. Methods for pet detector afterglow management

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3586982A (en) * 1969-02-19 1971-06-22 Westinghouse Electric Corp Pulse count detector having fail-safe features
US3696399A (en) * 1970-09-11 1972-10-03 Coulter Electronics Range expansion method and apparatus for multichannel pulse analysis
US4767929A (en) * 1986-10-06 1988-08-30 The United States Of America As Represented By The United State Department Of Energy Extended range radiation dose-rate monitor
US5258926A (en) * 1988-08-08 1993-11-02 Osterreichesches Forschungszentrum Seibersdorf Gmbh Method of measuring radiation for a radiation measuring device
US5377282A (en) * 1991-09-19 1994-12-27 International Business Machines Corporation Optical inspection system utilizing dynamic analog-to-digital thresholding
US5640109A (en) * 1995-10-27 1997-06-17 Mts Systems Corporation Pulse detector
US5928150A (en) * 1997-10-04 1999-07-27 Neoprobe Corporation System for locating and detecting a source of photon emissions
JP2004333326A (ja) * 2003-05-08 2004-11-25 Nippon Kankyo Chosa Kenkyusho:Kk 放射線測定器および放射線測定方法
EP1906211A1 (en) * 2005-07-22 2008-04-02 ICX Radiation GmbH Detector for the measurement of Ionizing radiation
CN102859391A (zh) * 2011-04-21 2013-01-02 株式会社东芝 电子设备以及输出方法
CN102338880A (zh) * 2011-06-02 2012-02-01 中国科学技术大学 核脉冲幅度数字化方法及系统
JP2014085183A (ja) * 2012-10-22 2014-05-12 Panasonic Corp 放射線測定装置、放射線測定方法及び放射線測定プログラム
CN108886594A (zh) * 2016-04-12 2018-11-23 雷神公司 具有模拟计数器元件的数字单位单元

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于雪崩二极管阵列的光子计数器研制;夏海峰等;《渭南师范学院学报》;20110215;第26卷(第02期);第40-43页 *
基于高速电流脉冲前放的数字式电荷积分型能谱仪系统设计;曾国强等;《湘潭大学自然科学学报》;20180815(第04期);第32-38页 *

Also Published As

Publication number Publication date
CN110568468A (zh) 2019-12-13

Similar Documents

Publication Publication Date Title
CN109558295B (zh) 一种性能指标异常检测方法及装置
JP4160275B2 (ja) エネルギー測定方法及び測定装置
WO2016110141A1 (zh) 一种闪烁脉冲的数字化方法
WO2017020747A1 (zh) 一种检测慢盘的方法及装置
JP7265791B2 (ja) パルス信号をサンプリングするための方法及びデバイス並びにコンピュータプログラム媒体
US8468161B2 (en) Determining a seasonal effect in temporal data
JP2010151843A (ja) 波形特徴化方法
US20140268140A1 (en) Pulse discriminator for particle counter
CN110568468B (zh) 一种辐射脉冲计数突变算法
CN109900382B (zh) 光纤温度信号解调方法、装置及光纤温度解调仪
JP6355491B2 (ja) 日射計性能低下状態推定装置、日射計性能低下状態推定システム及び日射計性能低下状態推定方法
CN114966470A (zh) 电力扰动检测方法、装置、计算机设备、存储介质
US4830020A (en) Measurement signal interference elimination
US20210080436A1 (en) Liquid chromatograph
US10819368B2 (en) Method for compressing and restoring time series data
CN109214447B (zh) 磁盘寿命预测方法及装置
TWI618003B (zh) 基於感測資料之估測方法及基於感測資料之估測系統
Kotov et al. CCD base line subtraction algorithms
JP5146226B2 (ja) 有義波高算出装置、プログラム、及び有義波高算出方法
WO2020196257A1 (ja) 測距方法、測距装置、及び、プログラム
CN108174196A (zh) 基于距离加权的成像系统调制传递函数测量方法
CN110196446B (zh) 一种高能粒子探测单元脉冲频率的快速平滑方法
KR101426753B1 (ko) 지문 영상 취득 장치 및 방법
CN113670369A (zh) 基于移动终端的风速测量及风噪声检测方法及装置
CN112016045A (zh) 纳秒级脉冲功率计数据处理方法

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