CN112882082B - 基于脉冲形状甄别的中子-γ周围剂量当量率仪的测量方法 - Google Patents
基于脉冲形状甄别的中子-γ周围剂量当量率仪的测量方法 Download PDFInfo
- Publication number
- CN112882082B CN112882082B CN202110034832.5A CN202110034832A CN112882082B CN 112882082 B CN112882082 B CN 112882082B CN 202110034832 A CN202110034832 A CN 202110034832A CN 112882082 B CN112882082 B CN 112882082B
- Authority
- CN
- China
- Prior art keywords
- gamma
- pulse
- neutron
- charge
- signal
- 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/16—Measuring radiation intensity
- G01T1/20—Measuring radiation intensity with scintillation detectors
- G01T1/208—Circuits specially adapted for scintillation detectors, e.g. for the photo-multiplier section
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
本发明公开了基于脉冲形状甄别的中子‑γ周围剂量当量率仪,包括有机闪烁体探测器、光电转换器件、信号处理电路、脉冲形状甄别电路和仪器控制电路;有机闪烁体探测器在中子或γ射线的照射下发出荧光,光电转换器件将荧光转换成电流信号,信号预处理电路将光电转换器件产生的电流信号转换成时间常数有差异的电压脉冲信号,脉冲形状甄别电路可根据电压信号的脉冲形状差异进行中子和γ甄别,并分别记录中子和γ产生的电荷积分谱,仪器控制电路将单位时间内中子和γ产生的脉冲的电荷积分分别换算成中子和γ周围剂量当量率。本发明可以用同一个探测器实现中子和γ周围剂量当量率的测量,可以大大减轻中子周围剂量当量率仪的重量。
Description
技术领域
本发明涉及核辐射探测技术领域,具体涉及一种基于脉冲形状甄别的中子-γ周围剂量当量率仪。
背景技术
中子和γ周围剂量当量率通常由不同的两种仪器进行测量。
中子周围剂量当量率仪一般由探头和信号采集处理系统两部分组成。探头一般由探测器和慢化体组成。探测器一般使用3He正比计数管、BF3正比计数管或6LiI闪烁体探测器等。中子与探测器内的灵敏物质作用后产生脉冲信号,信号采集处理系统将这些脉冲信号进行适当处理后计数,并将计数率换算成中子周围剂量当量率。慢化体用来对能量较高的中子进行慢化,使之能够被探测器探测到,起到优化中子周围剂量当量率仪能量响应的目的。现行的中子周围剂量当量测量装置的探头多为上述结构,其缺点是能量响应性能不尽如人意,且重量重、体积大。
γ周围剂量当量率仪也一般由探测器和信号处理系统组成,探测器一般是GM计数管、闪烁体探测器或半导体探测器中的一种。使用 GM计数管作探测器时一般需在GM计数管外增加能量补偿层,信号处理系统对GM计数管输出的脉冲进行计数或时间-计数处理得到剂量率结果。使用闪烁体探测器或半导体探测器作探测器时一般需对探测器输出的信号进行脉冲幅度分析,再利用G(E)函数转换成剂量率结果。
发明内容
针对目前中子和γ周围剂量当量测量需要两个不同仪器完成,且中子周围剂量当量测量装置体积大重量重的问题,本发明提出一种基于脉冲形状甄别的中子-γ周围剂量当量率仪。
为实现上述目的,本发明所设计的基于脉冲形状甄别的中子-γ周围剂量当量率仪,其特殊之处在于,所述剂量当量率仪包括:
有机闪烁体探测器:用于对于中子和γ的作用发出在时间特性上存在差异的荧光信号;
光电转换器件:用于将有机闪烁体探测器发出的荧光信号转换成电压或电流脉冲信号;
信号处理电路:用于对光电转换器件输出的电压或电流脉冲信号进行处理,转换为信号幅度和阻抗与脉冲形状甄别电路相匹配的电压脉冲信号,并保留中子和γ信号的时间特性的差异;
脉冲形状甄别电路:用于将信号处理电路输出的电压脉冲信号进行数字转换,根据中子和γ信号的时间特性的差异对其甄别,完成信号的电荷积分计算,分别得到中子和γ脉冲信号的电荷积分;
仪器控制电路:用于将脉冲形状甄别电路处理得到的中子和γ脉冲信号的电荷积分转换成中子和γ周围剂量当量结果。
进一步地,所述有机闪烁体探测器采用掺杂10B元素的闪烁体探测器。
更进一步地,所述脉冲形状甄别电路包括高速AD转换电路和 FPGA系统,所述高速AD转换电路用于对信号处理电路输出的电压脉冲信号进行数字化处理,所述FPGA系统用于采用脉冲形状甄别 PSD算法甄别中子和γ信号的时间特性的差异,脉冲形状甄别PSD算法为上升时间法、电荷比较法、脉冲梯度法和频率梯度法中的一种或若干种。
更进一步地,FPGA系统对经过数字化处理的脉冲信号判断信号幅度,当幅度超过设定触发阈值,确定为一次有效的粒子脉冲事件,并将该事件的数据缓存到FPGA内部FIFO存储器中,同时启动运算,根据设定的脉冲起始点确定方法自动寻找脉冲起始点位置和结束点位置,然后根据设定的长门时间窗Tlong、短门时间窗Tshort对时间窗内的电荷进行积分,分别得出两个时间窗内的电荷积分Qlong和Qshort;同时,对每个脉冲信号的起始点位置到结束点位置时间段内的电荷进行积分,得出每个脉冲信号的电荷积分Q;根据每一个脉冲的Qlong和Qshort、Q计算脉冲形状甄别PSD值,通过PSD值的取值区间分别对中子和γ计数,FPGA将计算结果发送给仪器控制电路。
更进一步地,所述仪器控制电路:在每秒内确定不同电荷积分Q 的γ计数分布,根据该分布和不同电荷积分Q所对应的γ电荷积分- 剂量转换函数Gγ(Q),计算出γ的剂量当量率;在每秒内确定不同电荷积分Q的中子计数分布,根据该分布和不同电荷积分Q所对应的中子电荷积分-剂量转换函数Gn(Q),计算出中子的剂量当量率。
更进一步地,所述脉冲形状甄别PSD值的计算方法为:PSD= (Qlong--Qshort)/Qlong;其中,Qlong为脉冲信号长门的电荷积分,时间选取脉冲信号的长度;Qshort为脉冲信号短门的电荷积分,时间选取脉冲信号的前沿。
更进一步地,计算γ的剂量当量率的公式为:
其中,i=1,2,3,…,n;n为有机闪烁体探测器探测到γ能谱的道数; Qi为有机闪烁体探测器探测到的一个电荷积分值;表示γ在有机闪烁体探测器探测到电荷积分值为Qi所对应的计数率,Gγ(Qi)表示该道所对应的伽马电荷积分-剂量转换函数Gγ(Q)的权重值。
更进一步地,计算中子的剂量当量率的公式为:
更进一步地,伽马电荷积分-剂量转换函数Gγ(Q)的函数表达式为:
其中k为指数函数关系式的项数,Ak为权重系数。
更进一步地,中子电荷积分-剂量转换函数Gn(Q)的函数表达式为:
Gn(Q)=aQ3+bQ2+cQ+d
其中a,b,c,d均为常数系数。
与现有技术相比,本发明具有的有益效果为:
(1)本发明可同时测量中子和γ周围剂量当量,根据电压信号的脉冲形状差异进行中子和γ甄别,计算中子和γ周围剂量当量率;
(2)本发明提出的电荷积分-剂量转换函数Gn(Q)和Gγ(Q)将单位时间内中子和γ产生的脉冲的电荷积分分别换算成中子和γ周围剂量当量率;
(3)本发明采用的有机闪烁体探测器的闪烁体材料中掺杂10B元素,以提高低能中子的响应;
(4)本发明可在保证测量的能量范围较大的情况下,体积和重量均有较大幅度的减小。
附图说明
图1是本发明的整体结构示意图;
图2是经过ADC采样的指数衰减脉冲信号示意图;
图3是中子-γ光脉冲信号形状区别示意图;
图4是电荷比较法示意图;
图5是本仪器实验中对某Am-Be中子源的甄别效果;
图6为MCNP单能中子的脉冲分布示意图;
图7为剂量转换系数曲线示意图。
具体实施方式
为了使本发明技术方案及有益效果更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。
如图1所示,本发明提出的基于脉冲形状甄别的中子-γ周围剂量当量率仪,包括有机闪烁体探测器、光电转换器件、信号处理电路、脉冲形状甄别电路、仪器控制电路。其中
有机闪烁体探测器:用于对于中子和γ的作用发出在时间特性上存在差异的荧光信号;
光电转换器件:用于将有机闪烁体探测器发出的荧光信号转换成电压或电流脉冲信号;
信号处理电路:用于对光电转换器件输出的电压或电流脉冲信号进行处理,转换为信号幅度和阻抗与脉冲形状甄别电路相匹配的电压脉冲信号,并保留中子和γ信号的时间特性的差异;
脉冲形状甄别电路:用于将信号处理电路输出的电压脉冲信号进行数字转换,根据中子和γ信号的时间特性的差异对其甄别,完成信号的电荷积分计算,分别得到中子和γ脉冲信号的电荷积分;
仪器控制电路:用于将脉冲形状甄别电路处理得到的中子和γ脉冲信号的电荷积分转换成中子和γ周围剂量当量结果。
有机闪烁体探测器为化学合成材料,可以是液态的或固态的,可对于中子和γ的作用发出在时间特性上存在差异的荧光信号。有机闪烁体探测器的不同成分的含量对于中子剂量测量性能有一定的影响,主要包括探测器的能量响应性能和n-γ甄别性能。常见的用于n-γ脉冲形状甄别的有机闪烁体有BC系列、NE系列、EJ系列等,不同型号的有机闪烁体探测器在性能上有所差异,但是在化学成份上是类似的,本实施例优选EJ339A型闪烁体探测器,其在有机闪烁体材料中掺杂10B元素,以提高低能中子的响应。
闪烁探测器的工作原理是:(1)放射性粒子打到探测器上,并在闪烁体内产生闪光;(2)所产生的闪烁光子经光导传输到光电倍增管的光阴极上转换成光电子;(3)光阴极上产生的电子经光电倍增管系统的多次倍增产生足够大小的输出信号;(4)光电倍增管输出电信号经配套电子仪器放大、处理而提供所需要的信息。
在本发明的一个实施方式中,光电转换器件可以是光电倍增管或 SiPIN二极管,本实施例采用光电倍增管作为光电转换器件,将光电倍增管与EJ339A型闪烁体通过硅脂耦合后封装在铝合金外壳内,外壳内配套有磁屏蔽套、内置高压电源模块和分压器电路。
信号处理电路用于对光电转换器件输出的电压或电流脉冲信号进行初步处理,使之成为在信号幅度和阻抗方面与脉冲形状甄别电路相匹配的电压脉冲信号,且保留中子和γ信号的时间特性的差异。本实施例的信号处理电路为一个程控放大电路,可在不改变中子和γ信号时间特性差异的同时对光电倍增管输出的信号进行适当放大,使电压脉冲信号的幅度在0~1V之间,满足后续电路AD转换的要求。
脉冲形状甄别电路由高速AD转换部分和数字化脉冲信息处理部分组成。高速AD转换部分将经信号处理电路处理的电压脉冲进行数字转换,使之成为数字信号。数字化脉冲信息处理部分是一种数字化信号处理系统,可以根据中子和γ产生的脉冲的电压脉冲信号的时间特性的差异对它们进行甄别,同时可完成信号的电荷积分计算,分别得到中子和γ脉冲产生信号的电荷积分。
本实施例的脉冲形状甄别电路由高速AD转换电路和FPGA系统组成。高速AD转换电路采用AD9680芯片,采样率1GSps,采样精度14位,采用JESD204B高速数字接口与FPGA系统通讯。FPGA系统采用XC7K325T芯片,用于处理AD转换后的数字信号,最大数据率为12.5Gb/s。本实施例利用FPGA中的算法对中子和γ的脉冲信号进行脉冲形状甄别,可采用上升时间法、电荷比较法、脉冲梯度法和频率梯度法等不同的PSD算法,本实施例优选电荷比较法作为FPGA 中运行的PSD算法。光电倍增管输出的信号通过信号处理电路后,产生快速的指数衰减脉冲信号,信号经过高速ADC数字化后,由FPGA 实时获取,经过ADC采样的指数衰减脉冲信号如图2所示(图示横轴单位ns,纵轴单位mV)。
FPGA获取到ADC实时数据之后,首先判断信号幅度,当幅度超过设定触发阈值,系统认为这是一次有效的粒子脉冲事件,并将该事件的数据缓存到FPGA内部FIFO存储器中,FPGA内部的脉冲形状甄别模块同时启动运算,脉冲形状甄别模块根据设定的脉冲起始点确定方法自动寻找脉冲起始点位置和结束点位置,然后根据设定好的长门时间窗Tlong,短门时间窗Tshort对时间窗内的电荷进行积分,分别得出两个时间窗内的电荷积分Qlong和Qshort。同时,对脉冲信号的起始点位置到结束点位置时间段内的电荷进行积分,得出每个脉冲信号的电荷积分Q;再根据每一个脉冲的Qlong和Qshort、Q计算脉冲形状甄别 PSD值,通过PSD的取值区间分别对中子和γ计数,FPGA将计算结果发送给仪器控制电路。
仪器控制电路是一种运行有嵌入式软件的微处理器系统,可将脉冲形状甄别电路处理得到的中子和γ脉冲的电荷积分转换成中子和γ周围剂量当量结果。本实施例的仪器控制电路是一个单片机系统,仪器控制电路在接收到一个脉冲的计算结果之后,利用单片机内的嵌入式程序计算出该脉冲的PSD值,仪器控制电路根据该脉冲的PSD 值确定该脉冲是γ还是中子作用产生的。
每秒(s)内确定的不同电荷积分Q的伽马计数分布,根据该分布和不同电荷积分Q所对应的伽马电荷积分-剂量转换函数Gγ(Q),单片机计算出伽马的剂量当量率。
每秒(s)内确定的不同电荷积分Q的中子计数分布,根据该分布和不同电荷积分Q所对应的中子电荷积分-剂量转换函数Gn(Q),单片机计算出中子的剂量当量率。
仪器控制电路采用脉冲形状甄别方法区分中子和γ射线。脉冲形状甄别是利用n-γ射线在闪烁体中产生的闪烁光脉冲信号形状存在差异这一特征实现中子和γ射线的甄别。对同一闪烁体而言,闪烁体内入射粒子激发的荧光脉冲分为快成分、慢成分,不同的入射粒子,快、慢成分的强度比也不同,中子的快成分份额小、慢成分份额大,而γ射线快成分份额大、慢成分份额小,因而n-γ射线通过闪烁体探测器探测到的光脉冲信号形状就有所区别,如图3所示,通过对这种光脉冲信号形状区别的处理分析,就可以实现n-γ射线的甄别。
本实施例采用电荷比较法的脉冲形状甄别方法。本方法中,提到的电荷积分是单个脉冲信号幅值对时间的积分。
电荷比较法将脉冲信号中瞬发光子持续所产生的前沿的积分值定义为快成分,由缓发光子持续所产生的后沿的积分值定义为慢成分,两者之和称为总电荷积分。如图4所示,电荷比较法使用长门、短门来确定信号积分的位置,本实施例脉冲形状甄别PSD值作为甄别依据,PSD定义如下:
PSD=(Qlong--Qshort)/Qlong
其中Qlong为脉冲信号长门的电荷积分,时间选取完整脉冲信号的长度;Qshort为脉冲信号短门的电荷积分,短门的积分时间选取信号波形的前沿。显然,Qlong即为总电荷积分,Qlong-Qshort即为慢成分电荷积分。由于中子比γ射线产生的脉冲信号的衰减速度慢,因而中子的脉冲形状甄别PSD值要比γ射线大。通过这样的定量区别即可实现 n-γ脉冲信号甄别。
电荷比较法实现过程中,根据脉冲信号起始点位置,设定好时间窗Tlong和Tshort,然后对时间窗内的电荷进行积分,分别得出两个时间窗内的电荷积分Qlong和Qshort。其中,Tlong和Tshort分别为长门、短门的时间域。
Tlong和Tshort的选取过程如下:当脉冲信号积分值第一次大于较小的数值(本仪器经实验摸索定义为30)时,定义该位置为起始位置时间t0,则长门时间结束位置tlong=t0+Δtlong,短门时间结束位置 tshort=t0+Δtshort;Δtlong和Δtshort由操作时设定、调整甄别效果最佳的数值,则可得到Tlong时间窗为[t0,tlong]、Tshort时间窗为[t0,tshort]。通过脉冲信号对时间各个时间窗的积分获得Qlong、Qshort。由PSD的计算公式可得每一个脉冲信号的脉冲形状甄别PSD值,通过PSD值大小分布的不同区别开中子和γ。
电荷比较法实现的Matlab程序语言。本仪器实验中对某Am-Be 中子源的甄别效果如图5所示。
区别中子、γ粒子后,仪器控制电路根据不同粒子产生的脉冲的电荷积分值Q,以及电荷积分-剂量转换函数Gn(Q)和Gγ(Q)函数可计算出中子和γ所产生的剂量当量率,计算方法如下:
脉冲的电荷积分-剂量转换函数Gγ(Q)
其中,i=1,2,3,…,n;n为探测器探测到γ能谱的道数;Qi为探测器内一个电荷积分值;N(Q)为探测器测得的γ能谱分布;表示γ在探测器探测到电荷积分值为Qi所对应的计数率;Gγ(Qi)表示该道所对应的Gγ(Q)函数权重值。由于实验上标准点源的个数有限,所以采用最小二乘法求解Gγ(Q)函数。使Gγ(Q)表示为指数函数关系式,如下:
k为指数函数关系式的项数,k=1,2,3,…;Ak为权重系数;则探测器根据Gγ(Q)函数计算得出的剂量率应为:
对于第j个标准γ源(j=1,2,3,…)其根据Gγ(Q)函数计算对应的剂量率为:
根据此等式,求出权重系数Ak各值,即可得出γ脉冲的电荷积分 -剂量转换函数Gγ(Q)。则探测器探测到的空间某点γ吸收剂量率为:
中子电荷积分-剂量转换函数Gn(Q)
测量中子周围剂量当量率时,通过使用反冲质子法测量中子。一般情况下,中子进入闪烁体后,与闪烁体内的氢原子发生弹性碰撞反应,设中子的能量为En,而反冲质子的能量为Ep,反冲角为θ,则存在以下关系式:
Ep=Encosθ
由于反冲质子的能量与反冲角有关,所以即便入射的是单能中子,反冲质子的能量仍然是连续分布的,且最大值等于En。反冲质子法是利用中子同氢核的弹性散射产生反冲质子,在含氢物质靶相当薄的条件下,某一能量入射中子的注量率同产生的反冲质子数成正比,测定反冲质子数和其他相关的量,就绝对测定了中子注量率。
本EJ339A液体闪烁体探测器中,产生的反冲质子引起闪烁体发光,若发光后经过电子学线路探测得到的脉冲信号电荷量为Q,则Q 与入射中子的能量En成正比。中子通过闪烁体探测器探测后得到的脉冲高度分布为P(Q),则探测器测量到的中子周围剂量当量率为:
实际中若使用多个单能中子源进行Gn(Qi)计算,难度较大。通过基于蒙特卡罗方法的MCNP程序F8卡模拟多个单能中子源对空间某点闪烁体探测器进行照射,获得不同能量为En(n=1,2,…,n;En-1<En) 单能中子源的脉冲高度分布矩阵P(Ep)(实际是也是反冲质子的脉冲高度分布矩阵)为:
P(Epi,En)为单能中子能量为En时探测器反冲质子能量为Ep(i-1)到Epi的概率。本实施中,将n和m最大值相等,则P(Ep)为方矩阵。通过MCNP计算EJ339A不同单能中子入射后的脉冲高度分布,如图 6所示。
设通过MCNP单能中子源计算出的脉冲高度分布对应的剂量转换函数为G(Ep)。矩阵G(Ep):
[G(Ep1) G(Ep2) G(Ep3) ...... G(Epm)]T
G(Epi)是脉冲高度大小为Ep(i-1)到Epi的光谱权重值。
通过MCNP模拟计算的单能中子源在探测器位置中子周围剂量当量标准值为H,标准值由单能中子注量和ICRP74号数据给出的注量- 剂量转换系数乘积得到。不同单能中子源的中子周围剂量当量标准值矩阵H为:
那么探测器根据G(Ep)计算的中子周围剂量当量和MCNP模拟计算的中子周围剂量当量标准值相等的,则有下列等式:
表示为矩阵方程:H=P(Ep)·G(Ep),即:G(Ep)=P-1(Ep)·H
通过此方程可求出G(Ep)各个元素值。剂量转换系数曲线如图7 所示。
对G(Ep)和Ep拟合为四次曲线,可得G(Ep)函数为:
G(Ep)=xEp 3+yEp 2+zEp+t
其中x,y,z,t均为常数系数。本实施例中,最佳拟合曲线的x,y,z,t 数值分别为:4.253,10.9593,-1.273,0.0731。
当然求出G(Ep)对于实测还需要进行试验标定。由于探测器实际探测到的脉冲信号电荷Q和反冲质子能量、入射中子能量都是正比关系,可设比例系数为K,那么就存在:Q=K˙Ep这样的关系,K为待定系数。显然,中子的剂量当量值通过MCNP和实验测得的两种结果也存在这样的比例关系。实际中,通过实验测量某一标准中子辐射场,根据不同位置剂量的不同,确定此待定系数。确定系数后,根据G(Ep) 函数拟合曲线的公式,将Ep值替换为Q即可得Gn(Q)函数式为:
G(Q)=aQ3+bQ2+cQ+d
其中a,b,c,d均为常数系数。由此即可实现探测器在实际测量中子时获得中子周围剂量当量率:
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
本说明书未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (5)
1.基于脉冲形状甄别的中子-γ周围剂量当量率仪的测量方法,其特征在于:所述剂量当量率仪包括:
有机闪烁体探测器:用于对于中子和γ的作用发出在时间特性上存在差异的荧光信号;
光电转换器件:用于将有机闪烁体探测器发出的荧光信号转换成电压或电流脉冲信号;
信号处理电路:用于对光电转换器件输出的电压或电流脉冲信号进行处理,转换为信号幅度和阻抗与脉冲形状甄别电路相匹配的电压脉冲信号,并保留中子和γ信号的时间特性的差异;
脉冲形状甄别电路:用于将信号处理电路输出的电压脉冲信号进行数字转换,根据中子和γ信号的时间特性的差异对其甄别,完成信号的电荷积分计算,分别得到中子和γ脉冲信号的电荷积分;
仪器控制电路:用于将脉冲形状甄别电路处理得到的中子和γ脉冲信号的电荷积分转换成中子和γ周围剂量当量结果;所述仪器控制电路:在每秒内确定不同电荷积分Q的γ计数分布,根据该分布和不同电荷积分Q所对应的伽马电荷积分-剂量转换函数Gγ(Q),计算出γ的剂量当量率;在每秒内确定不同电荷积分Q的中子计数分布,根据该分布和不同电荷积分Q所对应的中子电荷积分-剂量转换函数Gn(Q),计算出中子的剂量当量率;
其中,i=1,2,3,…,n;n为探测器探测到γ能谱的道数;Qi为探测器内一个电荷积分值;N(Q)为探测器测得的γ能谱分布;表示γ在探测器探测到电荷积分值为Qi所对应的计数率;Gγ(Qi)表示该道所对应的Gγ(Q)函数权重值,采用最小二乘法求解Gγ(Q)函数;使Gγ(Q)表示为指数函数关系式:
k为指数函数关系式的项数,k=1,2,3,…;Ak为权重系数;则探测器根据Gγ(Q)函数计算得出的剂量率应为:
对于第j个标准γ源(j=1,2,3,…)其根据Gγ(Q)函数计算对应的剂量率为:
根据此等式,求出权重系数Ak各值,得出γ脉冲的电荷积分-剂量转换函数Gγ(Q);
中子电荷积分-剂量转换函数Gn(Q)的计算方法为:对G(Ep)和Ep拟合为四次曲线,得G(Ep)函数为:G(Ep)=xEp 3+yEp 2+zEp+t,其中G(Ep)为脉冲高度分布对应的剂量转换函数,Ep为反冲质子的能量,x,y,z,t均为常数系数;进行试验标定,探测器实际探测到的脉冲信号电荷Q和反冲质子能量、入射中子能量都是正比关系,设比例系数为K,Q=K˙Ep,K为待定系数,通过实验测量标准中子辐射场,根据不同位置剂量的不同,确定此待定系数K后,根据G(Ep)函数拟合曲线的公式,将Ep值替换为Q得到Gn(Q)函数式为G(Q)=aQ3+bQ2+cQ+d其中a,b,c,d均为常数系数。
2.根据权利要求1所述的基于脉冲形状甄别的中子-γ周围剂量当量率仪的测量方法,其特征在于:所述有机闪烁体探测器采用掺杂10B元素的闪烁体探测器。
3.根据权利要求1所述的基于脉冲形状甄别的中子-γ周围剂量当量率仪的测量方法,其特征在于:所述脉冲形状甄别电路包括高速AD转换电路和FPGA系统,所述高速AD转换电路用于对信号处理电路输出的电压脉冲信号进行数字化处理,所述FPGA系统用于采用脉冲形状甄别PSD算法甄别中子和γ信号的时间特性的差异,脉冲形状甄别PSD算法为上升时间法、电荷比较法、脉冲梯度法和频率梯度法中的一种或若干种。
4.根据权利要求3所述的基于脉冲形状甄别的中子-γ周围剂量当量率仪的测量方法,其特征在于:FPGA系统对经过数字化处理的脉冲信号判断信号幅度,当幅度超过设定触发阈值,确定为一次有效的粒子脉冲事件,并将该事件的数据缓存到FPGA内部FIFO存储器中,同时启动运算,根据设定的脉冲起始点确定方法自动寻找脉冲起始点位置和结束点位置,然后根据设定的长门时间窗Tlong、短门时间窗Tshort对时间窗内的电荷进行积分,分别得出两个时间窗内的电荷积分Qlong和Qshort;同时,对每个脉冲信号的起始点位置到结束点位置时间段内的电荷进行积分,得出每个脉冲信号的电荷积分Q;根据每一个脉冲的Qlong和Qshort、Q计算脉冲形状甄别PSD值,通过PSD值的取值区间分别对中子和γ计数,FPGA将计算结果发送给仪器控制电路。
5.根据权利要求4所述的基于脉冲形状甄别的中子-γ周围剂量当量率仪的测量方法,其特征在于:所述脉冲形状甄别PSD值的计算方法为:PSD=(Qlong--Qshort)/Qlong;其中,Qlong为脉冲信号长门的电荷积分,时间选取脉冲信号的长度;Qshort为脉冲信号短门的电荷积分,时间选取脉冲信号的前沿。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110034832.5A CN112882082B (zh) | 2021-01-12 | 2021-01-12 | 基于脉冲形状甄别的中子-γ周围剂量当量率仪的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110034832.5A CN112882082B (zh) | 2021-01-12 | 2021-01-12 | 基于脉冲形状甄别的中子-γ周围剂量当量率仪的测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112882082A CN112882082A (zh) | 2021-06-01 |
CN112882082B true CN112882082B (zh) | 2022-03-29 |
Family
ID=76045273
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110034832.5A Active CN112882082B (zh) | 2021-01-12 | 2021-01-12 | 基于脉冲形状甄别的中子-γ周围剂量当量率仪的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112882082B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113848580A (zh) * | 2021-09-26 | 2021-12-28 | 中国科学院高能物理研究所 | 一种X/γ射线辐射量探测系统及探测方法 |
CN114035220A (zh) * | 2021-11-30 | 2022-02-11 | 成都天核科技有限公司 | 一种便携式能谱型n-γ个人剂量当量测量仪及其方法 |
CN114325809B (zh) * | 2022-01-05 | 2023-05-12 | 中国科学院高能物理研究所 | 基于电流积分型电子学系统的中子剂量仪 |
CN114839662A (zh) * | 2022-05-16 | 2022-08-02 | 哈尔滨工程大学 | 一种宽量程伽马剂量仪 |
CN115144886B (zh) * | 2022-09-05 | 2023-01-06 | 中国人民解放军军事科学院军事医学研究院 | 高通量中子剂量评估方法及装置 |
CN116660970B (zh) * | 2023-07-07 | 2024-03-22 | 安徽中能聚控科技有限公司 | 一种基于中子伽马混合辐射的脉冲甄别方法及系统 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI497528B (zh) * | 2012-08-07 | 2015-08-21 | Iner Aec Executive Yuan | 輻射監測器遊校系統及方法 |
CN104536032B (zh) * | 2014-12-23 | 2017-05-17 | 四川大学 | 基于支持向量机的液体闪烁体探测器n/γ射线甄别方法 |
CN105093261B (zh) * | 2015-08-20 | 2018-08-24 | 中国原子能科学研究院 | 一种数字化中子谱仪的中子、γ甄别系统 |
CN106873019B (zh) * | 2017-01-06 | 2019-04-05 | 中国科学院高能物理研究所 | 一种辐射剂量测量方法 |
CN107422359A (zh) * | 2017-05-16 | 2017-12-01 | 中国工程物理研究院材料研究所 | 一种基于液体闪烁体探测器的中子伽马剂量的测量方法 |
CN107907901B (zh) * | 2017-11-06 | 2020-03-17 | 中广核工程有限公司 | 核电站中子、伽马能谱与剂量的测量方法和系统 |
CN108646284B (zh) * | 2018-05-15 | 2020-11-13 | 张金钊 | 一种γ能谱组合探测系统及γ能谱测量方法 |
CN111221030A (zh) * | 2019-12-03 | 2020-06-02 | 北京卫星环境工程研究所 | 基于物理集成的中子伽马探测器及中子伽马在线甄别方法 |
CN113466913A (zh) * | 2021-05-27 | 2021-10-01 | 中国工程物理研究院材料研究所 | 基于液体闪烁体探测器的混合辐射场剂量测量系统 |
CN113466915A (zh) * | 2021-08-18 | 2021-10-01 | 内蒙古民族大学 | 一种基于高纯锗探测器的辐射剂量率测量系统及方法 |
-
2021
- 2021-01-12 CN CN202110034832.5A patent/CN112882082B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112882082A (zh) | 2021-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112882082B (zh) | 基于脉冲形状甄别的中子-γ周围剂量当量率仪的测量方法 | |
US4217496A (en) | Portable instrument for measuring neutron energy spectra and neutron dose in a mixed n-γ field | |
EP2773982B1 (en) | Neutron spectrometer | |
US5298756A (en) | Scintillator fiber optic long counter | |
JP5832892B2 (ja) | 孔内検層のための方法および装置 | |
US4217497A (en) | Portable instrument for measuring neutron energy spectra and neutron dose in a mixed n-γ field | |
Dietze | Energy calibration of NE-213 scintillation counters by δ-rays | |
Parks | The Compton effect-Compton scattering and gamma ray spectroscopy | |
Jie et al. | Energy calibration of a BC501A liquid scintillator using a γ-γ coincidence technique | |
Luo et al. | Test of digital neutron–gamma discrimination with four different photomultiplier tubes for the NEutron Detector Array (NEDA) | |
JP2001235546A (ja) | 放射性ガス測定装置と燃料破損検出システム | |
CN108535758B (zh) | 一种脉冲形状甄别算法 | |
Bourne et al. | Digital pile-up rejection for plutonium experiments with solution-grown stilbene | |
CN116381772A (zh) | 一种用于剂量率仪的实时能量响应校正方法、系统及终端 | |
Archer et al. | An algorithmic approach to single‐probe Cherenkov removal in pulsed x‐ray beams | |
Pywell et al. | Photon flux monitor for a mono-energetic gamma ray source | |
Farsoni et al. | A system for simultaneous beta and gamma spectroscopy | |
Xu et al. | A well typed phoswich detector consisting of CsI and plastic scintillators for low level radioactivity measurements | |
CN114740520A (zh) | 一种放射性惰性气体活度测量装置及方法 | |
Ren et al. | Introduction of a C6D6 detector system on the Back-n of CSNS | |
Potapov et al. | A combined spectrometric detector of fast neutrons | |
Li et al. | Determination method of high fluence rate for DT neutron source with long counter | |
Williams et al. | The CsI ball ancillary detector array for TIP and TIGRESS at TRIUMF | |
Korun | Measurement of the total-to-peak ratio of a low-energy germanium gamma-ray detector | |
Maspero et al. | A real time scintillating fiber time of flight spectrometer for linac photoproduced neutrons |
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 |