CN1103053C - 核谱学系统中用于调整和甄别脉冲形状的方法 - Google Patents
核谱学系统中用于调整和甄别脉冲形状的方法 Download PDFInfo
- Publication number
- CN1103053C CN1103053C CN97121430.1A CN97121430A CN1103053C CN 1103053 C CN1103053 C CN 1103053C CN 97121430 A CN97121430 A CN 97121430A CN 1103053 C CN1103053 C CN 1103053C
- Authority
- CN
- China
- Prior art keywords
- pulse
- value
- sampled value
- threshold
- sampled
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 53
- 238000004958 nuclear spectroscopy Methods 0.000 title description 8
- 230000001105 regulatory effect Effects 0.000 title description 2
- 238000005070 sampling Methods 0.000 claims abstract description 42
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 230000005251 gamma ray Effects 0.000 claims description 22
- 239000002245 particle Substances 0.000 claims description 14
- 230000000052 comparative effect Effects 0.000 claims 2
- 230000003993 interaction Effects 0.000 claims 2
- 238000001514 detection method Methods 0.000 abstract description 11
- 230000010354 integration Effects 0.000 abstract description 8
- 238000005516 engineering process Methods 0.000 description 14
- 230000035508 accumulation Effects 0.000 description 13
- 238000009825 accumulation Methods 0.000 description 13
- 230000008569 process Effects 0.000 description 10
- 230000006870 function Effects 0.000 description 7
- 238000001228 spectrum Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000006872 improvement Effects 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000007493 shaping process Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- LBDSXVIYZYSRII-IGMARMGPSA-N alpha-particle Chemical compound [4He+2] LBDSXVIYZYSRII-IGMARMGPSA-N 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 230000003278 mimic effect Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000004904 shortening Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- MYMOFIZGZYHOMD-UHFFFAOYSA-N Dioxygen Chemical compound O=O MYMOFIZGZYHOMD-UHFFFAOYSA-N 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- TVFDJXOCXUVLDH-RNFDNDRNSA-N cesium-137 Chemical compound [137Cs] TVFDJXOCXUVLDH-RNFDNDRNSA-N 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005253 cladding Methods 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000012850 discrimination method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 description 1
- 229910052721 tungsten Inorganic materials 0.000 description 1
- 239000010937 tungsten Substances 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V5/00—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
- G01V5/04—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging
- G01V5/08—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays
- G01V5/12—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays using gamma or X-ray sources
- G01V5/125—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays using gamma or X-ray sources and detecting the secondary gamma- or X-rays in different places along the bore hole
-
- 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/17—Circuit arrangements not adapted to a particular type of detector
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Measurement Of Radiation (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
一种检测、处理、调整和识别核谱系统中信号脉冲形状以区分各种长度脉冲和堆积脉冲从而避免数据丢失和核谱失真的方法，该方法基于数字脉冲采样和检测技术，以与随时出现的实际脉冲异步的固定频率对一个核脉冲群采样。分析表示核粒子或伽马射线与核探测器之间相互作用的脉冲形状从而确定脉冲长度。探测核粒子或伽马射线，形成含有表示所说粒子或伽马射线的脉冲的一个信号；以固定频率对脉冲采样，以形成脉冲群的数字图象；通过对大于一个阈值的连续采样的第一数目的脉冲数计数和通过对大于该阈值的连续采样的第二数目的脉冲数计数来确定脉冲长度；求得连续采样的第一和第二数目的比值；用该比值来调整脉冲长度。
Description
本发明涉及一种探测和处理脉冲信号的方法，具体地说，本发明涉及在核谱学系统中用于调整和甄别信号脉冲形状的方法。
在核谱学应用中，需要测量入射的核粒子能量。在许多情况下是通过测量这些粒子在核探测器中作用的能量来完成这种测量的。这种能量是以一个连续的脉冲群形式体现的，这个脉冲群的长度在理论上是无限长的。为了得到表示作用在探测器中的总能量的一个信号，通常需要对一个电流信号进行积分。一般来说，在电荷积分时间和对处理高脉冲速率的需要之间必须找到一种折衷方案。通常，对于一组不同脉冲以低计数率进行探测很少产生问题，所以可以使用常规的探测装置。但是，有时脉冲的频率或者重复率会在很宽的范围内变化，使得连续脉冲之间的间隔有时变得非常短。因此，这些信号的随机特征和产生的高速率在放大器中必然产生一组“堆积”或者一序列重叠的脉冲。这种堆积常常形成一个由两个或多个经过放大的单一探测信号构成的脉冲，其中所说的每个探测信号表示探测到的伽马射线、中子或者其它核辐射。这种堆积现象导致数据丢失和/或能谱失真。所以，重要的是，首先要从堆积脉冲中识别出各个脉冲，其次当已经探测到堆积脉冲时，适当地处理这些脉冲以恢复原来不同的脉冲，或者至少剔除那些堆积而成并且不表示单个粒子作用在核探测器中的能量的脉冲。
这类脉冲是利用核谱分析仪或脉冲高度分析仪进行分析的。一种核谱分析仪可以包括一个闪烁体探测器、一个光电倍增管、一个耦合器(通常是一个电容器)、一个前置放大器、一个脉冲整形电路和一个脉冲高度分析仪。已知的脉冲高度分析仪包括相互连接的一个脉冲探测器(或者是一个堆积探测器和一个堆积处理电路)、一个模数转换器(ADC)和一个存储器，其不同的信道对应于所探测脉冲的一个给定幅值；脉冲高度分析仪还可以包括设置在ADC上游的一个输入门，用于当ADC忙时，即当该ADC正在处理一个探测脉冲时，阻止其它脉冲进入ADC。
如上所述，如果两个脉冲在同一积分间隔内到达，就会产生一个偏置的“加和”脉冲(堆积脉冲)。这种堆积不仅导致能谱失真，而且使得无法探测到其中任何一个脉冲。脉冲形状也会影响到堆积脉冲的探测。所以，利用一个脉冲整形装置来兼顾积分时间和计数率性能，该装置设置在脉冲数字化和脉冲在能谱(存储器)中累加之前。新型高速电子电路提供了高速率核谱探测性能。但是，这种系统对于脉冲形状的变化更加敏感，并且要求电子信号的形状在较大范围的工作状态下保持不变。能够影响脉冲形状的主要因素是电子元器件特性的改变和核传感器随时间和/或温度的变化。
有一种能谱分析方法就是在授予Bronislow Seeman的美国专利US-5067090中介绍的数字积分技术。按照这种核谱分析技术，对于以其幅度作为粒子，如利用所说辐射探测器捕获的例如伽马射线等粒子能量的一种量度的脉冲的脉冲高度分析是通过以下步骤进行的：(1)以给定速率连续地(异步地)将所探测到的信号转换成数字采样信号；和(2)处理每个数字采样信号以形成每个所探测到的核事件的数字图象。该方法包括探测脉冲到达的步骤，即通过将每个输入采样信号与一个阈值比较，从而确定该采样信号是否表示一个脉冲。更具体地说，是计算出输入采样值和与核事件无关的基本信号之间的差值，并与所说阈值比较。记录任一探测到的脉冲的到达时间。此外，连续地测算基本信号值(基线)，并且在每一采样时间到时予以更新，从而产生一个当前的基本信号值；可取的是，该更新值为输入采样值与以前的采样值的加权平均值。一旦探测到一个脉冲，就通过加和每个采样值(表示该脉冲)与当前的基本信号值之间的差值计算出所说脉冲的能量，所得加和值连续地累计在一个寄存器中。在能量计算过程中，将累计的加和值与一个预置值比较，在溢出的情况下重新设置所说寄存器。可取的是，对于每个脉冲，在所说的脉冲能量计算过程中，将恰好位于表示该脉冲的第一个采样值之前的那个采样值，以及恰好位于表示该脉冲的最后一个采样值之后的那个采样值都考虑在内。
Seeman的专利还包括探测彼此之间在时间上接近的连续脉冲的堆积的步骤。这个探测步骤包括对表示所探测的核事件的采样值的数目计数和将该计数与一个预定的最大计数值比较的步骤。该采样步骤可以利用一个快速模数转换器(ADC)进行。
在Seeman专利中所述的数字积分技术，由于它实质上是计算脉冲下面的面积，所以对于脉冲形状是十分敏感的。通过改变电子元器件的温度或使用时间以及改变核探测器的脉冲响应特性可以影向到脉冲形状。大部分闪烁体探测器在高于室温的高温下会输出快得多的脉冲。
稳定脉冲形状的一种可行的技术记载在授予R.Adolph和B.A.Roscoe的美国专利US-5132540中。这种方法通过按照给定的关系保持堆积脉冲数目和将这个数目与总计数比较来调整脉冲形状。这种方法具有在低计数率情况下比较慢和不可靠的缺点，其中脉冲堆积的总计数是非常小的。
其它核测量应用需要通过辨别从核探测器输出的脉冲形状来区分相互作用的种类。当前，大部分脉冲形状甄别工作依赖于复杂而通常并不稳定的模拟电路。不幸的是，这些技术的使用需要相当多的专业知识。因此，仍然需要有一种能够对不同长度的脉冲和堆积脉冲进行甄别的简单的方法。
本发明的一个目的是开发一种能够调整核谱学系统中脉冲形状的方法。
本发明的第二个目的是开发一种能够区分正常脉冲和堆积脉冲的方法。
本发明的再一个目的是提供一种能够高速率地和以高精度分析信号脉冲形状的方法。
本发明基于数字采样技术，该方法以恒定速率对核脉冲信号群进行采样。对于每个脉冲能够获得许多采样值的事实使得许多过去通常是利用复杂的模拟和数字电路进行的脉冲群分析操作，例如：(1)脉冲积分(确定脉冲下面的面积)，(2)脉冲形状甄别，(3)剔除堆积脉冲和(4)确定信号到达时间(多标度分析)，可以利用一种紧凑的数字装置来进行。本发明涉及数字脉冲群分析的前三个方面。
数字采样和积分技术是获取核谱学数据的一种新型方法。这技术的主要优点是具有非常高的速度和取消了核模数转换器(ADCs)中传统的采样和保持电路。图2表示了数字积分器的基本原理。按照给定的采样频率对从探测器输出的信号进行异步采样。如果一个脉冲超过了预定的数字阈值，就会触发脉冲探测逻辑电路。这种方法的一种可能的实施方式是将大于所说阈值的所有采样值与在该脉冲之前的第一个采样值和在该脉冲之后的第一个采样值相加。将大于阈值的采样值数目多于预置值(例如6)的脉冲剔除，并作为堆积脉冲计数。还根据在一个脉冲期间内所测得的采样值数目进行脉冲形状补偿(延长或缩短脉冲)。
剔除堆积脉冲是以类似于美国专利US-5067090中所述的一种原理为基础的，这种方法利用脉冲形状甄别采区分正常脉冲和堆积脉冲。按照这种方法，正常脉冲具有2-6个大于阈值的采样值。任何具有多于6个大于阈值的采样值的脉冲都被认为是一个堆积脉冲。在本发明中，利用包含其它信息的算法扩展了上述方法。本发明的算法能够完成下述任务，例如：(1)使指示堆积产生的脉冲长度作为脉冲幅度的函数，(2)使堆积剔除阈值(超过阈值的采样值数目)与脉冲峰值建立联系，(3)利用从阈值到峰值的时间确定脉冲的存在，(4)以预定的方式利用一种归纳方法，这种方法假定：对于一个给定脉冲，采样值与以前的采样值相关，和(5)利用上述任务的结合。本发明提供一种处理脉冲的改进方法和调整脉冲形状的直接方法。本发明还能够区分由中子和(不希望的)伽马射线产生的信号(脉冲形状甄别)或者能够实现在核谱学中需要区分不同脉冲形状的其它任务。由探测器产生的中子和伽马射线信号的脉冲形状是彼此不同的。两种成分的幅值之间的比率依赖于闪烁体材料中粒子的电离密度。因此，按照一种数字方法，可以通过观察最大脉冲电压与后面一个时间点的电压的比率区别这些信号。
图1为本发明使用的密度记录仪的示意图。
图2表示一个采样脉冲。
图3a表示补偿不足的采样脉冲。
图3b表示具有正确和所需形状的的采样脉冲。
图3c表示过补偿的采样脉冲。
图4为调整脉冲形状的算法的流程图。
图5表示低计数率(RR)情况下调整比值的一个实例，该调整比值是在脉冲形状补偿恒定设置时计数率的函数。
图6表示指示堆积的脉冲长度，该脉冲长度是脉冲幅度的函数。
图7表示本发明使堆积剔除阈值与脉冲峰值建立关系。
图8表示本发明利用从阈值到峰值的时间确定脉冲存在。
图9a表示本发明中，对于一个给定脉冲形状，采样值与以前的采样值相关。
图9b为本发明的一个流程图，其中对于一个给定的脉冲形状，采样值与以前的采样值相关。
图10表示在脉冲形状甄别基础上剔除堆积的一种普通方法的流程图。
图11表示在利用液体闪烁体的中子谱仪中脉冲形状甄别的三种响应。
本发明提供一种对由一种仪器诸如图1所示的密度记录仪探测的信号进行脉冲形状分析的快速方法。这种仪器包含一个伽马射线源10，一个近距离设置的伽马射线探测器ND，一个中距离设置的伽马射线探测器MD和一个远距离设置的伽马射线探测器FD。这三个伽马射线探测器是衰减型的，在现有技术中这种伽马射线探测器都设置在距离伽马射线源足够远的位置(远大于伽马射线从射线源到探测器的平均自由路径长度，这个距离通常为5厘米)使得计数率随着构成和与之相邻的覆盖材料的密度增加而降低。探测器ND、MD、和FD分别位于距射线源10逐渐增大的距离处。射线源10可以是任何适合的伽马射线源例如铯137。在探测器周围围绕着一层高密度屏蔽材料11例如钨。信号分析仪NA、MA和FA设置在探测器上面。每个分析仪分析来自相应探测器的信号。
参见图2，以预定的采样频率在每个脉冲期间从信号15中异步地选取采样值16。由于信号是异步采样的，需要有探测核探测器中脉冲出现的步骤和确保该脉冲既不是由于电子部分或者探测器噪声产生，也不是由于在短时间间隔内产生的两个或多个脉冲的堆积而产生的步骤。如果采样信号大于一个预定(数字)阈值电平17，并且有2个或多个连续的采样值保持在这个阈值之上，则认为存在一个脉冲。如果有6个以上的连续采样值保持大于阈值，则认为该脉冲是堆积脉冲。
上述脉冲探测和甄别方法是针对具有固定形状的脉冲的。但是，脉冲形状并不总是固定的，所以在某些情况下需要调整脉冲波形。能够引起脉冲形状产生变化的两种原因是电子元器件随时间和温度的变化，更主要的是探测器特性的变化，例如Nal闪烁晶体发射具有一定时间延迟的光。在大约21℃(70°F)的室温下该时间常数大约为240Ns。随着温度的增加，这个时间常数变短。这种时间缩短导致产生更短的“过补偿”脉冲。图3表示各种形状的脉冲。在图3a中，由于脉冲补偿不足，该脉冲持续时间太长。在图3b中，脉冲经过补偿使得该脉冲具有如图2所示的所需形状。图3c表示一个过短的脉冲(过补偿的)。图3a和图3c中的脉冲形状必须予以调整以获得正确形状。通过调节脉冲整形放大器中的脉冲整形滤波器可以实现这种修正。
根据脉冲长度进行脉冲形状补偿
脉冲形状可以根据脉冲长度进行补偿。对于给定的脉冲高度，在没有噪声的情况下，大于阈值的采样值数目总是恒定的。由于采样过程的异步特性存在少许偏差。但是，如果有噪声重叠在该脉冲上，脉冲长度的分布就会变得较长。这个结果意味着即使对于固定形状的脉冲，也会存在具有不同测量长度的脉冲。但是，对于给定的脉冲高度，大于阈值的平均采样值数目是脉冲长度的最佳量度。于是，可以将该平均脉冲长度用于一个调整回路中以保持脉冲形状固定。此外，有可能在一定的脉冲幅度范围内(而不是在一个幅度上)进行脉冲长度测量以增加高频脉动量和增加计数统计量。
图4表示这种脉冲长度补偿方法的流程图。在这个流程图中，是从探测信号20中得到光谱数据的。方框21对这个信号采样，并对采样的数目计数以确定大于预定阈值幅值的连续采样值的数目。确定五个采样脉冲的数目和六个采样脉冲的数目。方框22确定五个采样脉冲与六个采样脉冲的比值。在方框23，所确定的比值与一个预先确定的调整比值进行比较。如果在步骤24判断出五个采样脉冲与六个采样脉冲的比值大于该调整比值，那么在步骤25应将这些脉冲拉长。如果在步骤24判断出五个采样脉冲与六个采样脉冲的比值不大于该调整比值，则在步骤26缩短这些脉冲。在将脉冲缩短或拉长之后，在步骤27清除缓存器，并利用所获取的数据重复该方法。
参见步骤23，可以通过实验或通过理论计算和模拟确定调整比值(RR)。在本发明中，调整比值不需要是一个常数，但是应当是总计数率的函数。这是由于随着计数率的增加，堆积脉冲的数目增加，进而导致未探测的堆积脉冲数目增加。
堆积脉冲是由在时间上没有完全分开，从而被探测为一个脉冲的两个独立脉冲重叠构成的。堆积脉冲平均来说长度大于正常脉冲。这意味着随着计数率的增加，六个采样长脉冲的数目也增大。对于固定的脉冲形状，这种现象使调整比值减小。于是固定脉冲形状的调整比值还成为计数率的函数(对光谱形状的依赖程度更小)。因此在一个采样的情况下，调整比值变为：
RR＝RR0+c*TCR
其中TCR为总计数率，RR0为低计数率的调整比值，c是一个常数，RR为采用计数率修正的调整比值。图5表示对于脉冲形状补偿的一个恒定设置比值RR作为计数率的函数的一个实例。在图5中，垂直轴为调整比值，水平轴为总采样计数率。用标号28和29表示的调整比值曲线是实际测量得出的。它表明在这种情况下，调整比值是计数率的线性函数。
脉冲形状甄别
核谱学中的某些应用需要根据在核探测器中产生的电信号的形状区分这些信号。通常这是利用复杂的模拟电路完成的，这种电路需要进行大量的细微调节，并且易于发生漂移。数字脉冲采样的原理使得这种应用可以采用数字技术和数字处理。下面的一节内容给出脉冲形状甄别的一些例子，更概括地说是介绍甄别技术。
根据脉冲形状甄别剔除堆积脉冲
如上所述，在美国专利US-5067090中所述的数字积分器采用了一种脉冲形状甄别技术以区分正常脉冲和堆积脉冲。这种甄别技术的简单原理如下所述：
-正常脉冲具有2-6个大于阈值的采样值
-任何具有多于6个大于阈值的采样值的脉冲就是一个堆积脉冲。
这种已获成功的技术，可以扩展到那些要求根据比脉冲长度更复杂的判据进行形状甄别和不仅仅包含脉冲长度判据的应用中。
如图6所示的一种应用将表示堆积产生的脉冲长度作为脉冲幅值的函数。在这个实施例中，确定一个具体脉冲是否为堆积脉冲所需的采样数目依赖于脉冲29a的幅值29。如图所示，脉冲30具有幅值10，脉冲31具有幅值6。脉冲30可以具有用于判断脉冲堆积的六个采样限值，脉冲31可以具有用于判断脉冲堆积的四个采样限值。这反映出这样的事实，即具有较小幅值的脉冲具有较少的大于一个预定阈值的采样值。在实施这个实施例之前，首先确定一组预定幅值和相应的采样限值。在上述过程中存取这组幅值和相应的采样限值。在将阈值设置为一个相当高的值以便能够不受噪声的影响进行甄别的情况下，这是十分重要的。
图7表示将堆积脉冲剔除阈值与脉冲峰值建立联系的技术的应用结果。在这个实施例中，检测脉冲采样值的阈值电平正比于该脉冲幅值。在这一应用中首先利用一个预定的阈值检测该脉冲。但是，当确定该脉冲幅值后，又将该阈值设置为所检测幅值的一部分。这样确保了对于较小幅值的脉冲能够剔除适合的堆积脉冲。在图7中，阈值34适合于脉冲幅值36。对于具有该幅值一半的脉冲，阈值应当减小到一个较低的电平值35。这个较低的电平值35必须有所限制以避免由于噪声影响造成错误剔除。这要求将脉冲检测(阈值与之相交)之前的几个采样值存储起来，从而可以利用新的阈值重新分析该脉冲群。
在本发明的另一实施例中，利用从阈值到峰值的时间确定单个脉冲的存在。脉冲堆积常常使到达脉冲峰值的时间较长。于是还可以根据脉冲检测获得的到达峰值时间过长的结果执行堆积脉冲剔除。图8表示两个脉冲，一个具有正确的形状38，另一个37具有不正确的脉冲形状。到达峰值的时间可以用不同的方法测量：
1、可以确定具有最大幅值的采样值的量级。这是一种非常粗略的到达峰值时间测量方法，
2、可以确定峰值的矩心：
矩心＝(求和(Hi *i))/n
其中i为大于阈值的采样值的数目。H为该采样值的高度，n是大于阈值的采样值的总数目。按照一种不同的方法，在阈值计算中可以包括在阈值之前和之后的预定数目的采样值。那些矩心处于预定范围中的脉冲被接受，并且将它们的面积和/或它们的到达峰值时间存储在不同的直方图中。剔除其它脉冲并将它们的面积和/或到达峰值的时间存储在各直方图中可以利用平均的到达峰值时间调整脉冲形状。
图9a和图9b表示一种更通用脉冲形状甄别方法。按照这种方法，将一个脉冲的采样值与由前面的采样值预测的值连续比较。在该实施例中检测大于阈值的第一个采样值。然后预测第二个采样值的允许值的范围。将第二个采样值与该预测值比较。如果该采样值不在预测范围内、就对其作上标记(堆积)。重复这一过程直到一个采样值低于该阈值。如果任何一个采样值都在预测范围以外，则将该脉冲标记为一个堆积脉冲。根据脉冲的类型，使“正常”或“堆积”直方图累计增加。在某些应用中只对堆积脉冲计数，而不将它们存储在直方图中就足够了。
图10表示了用于脉冲形状甄别的一种更复杂的方法。按照这种方法根据许多条件区分脉冲，并进行分类。这使得在核粒子探测器中能区分不同事件的不同脉冲形状。
在另一个不同的实施例中将所获得的脉冲形状与一种预定形状进行比较。这可以在比较两个脉冲之前将它们的面积进行归一化处理而实现。如果平均偏差、均方差或其它适合的量度处于预定的限值内，则认为这些脉冲是一致的。通过构成预定脉冲形状与被检测形状之间的相互关联可以进行更加复杂的分析。利用一组脉冲形状能够区别许多脉冲形状从而剔除堆积脉冲和/或进行脉冲形状甄别。
甄别脉冲形状以区别核信号
在利用液体闪烁体进行中子探测时，常常需要区分由中子产生的信号和由伽马射线产生的不希望有的信号。由探测器产生的这些信号在脉冲形状上是彼此不同的。图11示意性表示了对α粒子46、快中子(质子)47和伽马射线48的三种测量结果。垂直轴表示光发射强度，水平轴表示时间。所有三种脉冲都由至少两种成分组成。两种成分幅度之间的比值依赖于这些粒子在闪烁体材料中的电离密度。所以，按照一种数字方法，可以通过观察最大的脉冲电压与在一个较晚时间点(例如300ns之后)的电压之间的比值采区分这些信号。与此同时，在不影响甄别中子和伽马射线或α粒子的前提下可以对该脉冲进行适合的积分，即汇集闪烁体发出的大部分光。
尽管针对一种具体的探测器描述了本发明的装置和方法，但是本发明的原理可以应用于许多其它的记录仪器，例如中子密度仪、碳/氧探测仪以及许多其它非核探测器。本发明的装置和方法与现有技术相比具有明显的优点。前面已经结合优选实施例介绍了本发明。但是，本发明并不局限于此。在不脱离本发明的发明构思的前提下可以对本发明的基本设计进行改变、改进和改善。此外，对于本领域技术人员来说在得益于上述的教导之后显然可以作出这些改变、改进和改善。所有这些改变、改进和改善都在本发明的范围之内，本发明的范围由所附的权利要求书限定。
Claims (7)
1、一种分析包含在脉冲群中的电脉冲形状的方法，所说脉冲表示核粒子或伽马射线与核探测器之间的相互作用，该方法包括以下步骤：
a)探测所说粒子或伽马射线，产生含有表示所说粒子或伽马射线的脉冲的一个信号；
b)以固定频率对在步骤a)中得到的脉冲采样，以形成该脉冲群的数字图象；
c)通过对大于一个预定阈值电平的连续采样值的数目进行计数来测量该脉冲长度；
d)计算大于所说阈值电平的脉冲的平均数目；和
e)将所说脉冲的平均数目与一个预定的数目比较，并根据测量的脉冲长度与所需脉冲长度之间的比较结果调整该脉冲长度。
2、如权利要求1所述的方法，其特征在于：步骤c)还包括通过对具有第一预定数目的大于所说阈值电平的连续采样值的脉冲数目进行计数来测量该脉冲长度，该方法还包括以下步骤：
f)对具有第二预定数目的大于所说阈值电平的连续采样值的脉冲的第二数目进行计数；
g)确定所说脉冲的第一数目与所说第二数目的比值；和
h)将该比值与一个期望比值比较，并根据比较结果调整脉冲长度。
3、一种分析包含在脉冲群中的电脉冲形状的方法，所说脉冲表示核粒子或伽马射线与核探测器之间的相互作用，该方法包括以下步骤：
a)探测所说粒子或伽马射线，产生含有表示所说粒子或伽马射线的脉冲的一个信号；
b)以固定频率对在步骤a)中得到的脉冲采样，以形成该脉冲群信号的数字图象；
c)确定在一组大于一个预定阈值的连续采样值中最大采样值的幅值；
d)将对应于一个有效采样值的阈值调节为等于在步骤c)中确定的幅值的预定一部分；
e)将具有多于预定数目的幅值大于所说阈值的采样值的脉冲识别为一个堆积脉冲；和
f)剔除所述难积脉冲。
4、如权利要求3所述的方法，还包括以下步骤：
g)确定该脉冲之下的面积，该面积值是所有大于所说阈值的连续采样值与在该阈值之前和之后的预定数目的采样值之和；
h)将所测量的面积与一个或多个预定脉冲形状进行比较；和
i)如果所测量面积与所说的一个或多个预定脉冲形状基本不一致，则将该脉冲识别为一个堆积脉冲。
5、如权利要求3所述的方法，还包括以下步骤：
g)确定该脉冲之下的面积，该面积值是所有大于所说阈值的连续采样值与在该阈值之前的预定数目的采样值之和；
h)对于一个采样值的预定数目n，预测当包含第n+1个脉冲时所得到的面积值；
i)将包含第n+1个采样值时测得的面积与预测的包含第n+1个采样值时的面积比较；和
j)如果所测量的面积不在预测面积的预定范围之内，则将该脉冲识别为一个堆积脉冲。
6、如权利要求3所述的方法，其特征在于：步骤(c)还包括确定在大于预定阈值的第二组连续采样值中的一个第二采样值的幅值，并求出最大采样值的幅值与所说第二采样值的幅值的比值的步骤；以及，步骤(e)还包括如果该比值小于一个预定数值则将该脉冲识别为一个堆积脉冲的步骤。
7、如权利要求3所述的方法，其特征在于：根据大于阈值的采样值的一个预定数目，步骤(c)还包括以下步骤：测量该脉冲达到峰值的时间，并判定所测得的达到峰值时间是否落入达到峰值时间的一个预定范围内；以及，步骤(e)还包括如果该达到峰值时间测量值不在所说的预定范围内则将该脉冲识别为一个堆积脉冲的步骤。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US705,389 | 1991-05-24 | ||
| US705389 | 1996-08-30 | ||
| US08/705,389 US5884234A (en) | 1996-08-30 | 1996-08-30 | Method for pulse shape regulation and discrimination in a nuclear spectroscopy system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1182883A CN1182883A (zh) | 1998-05-27 |
| CN1103053C true CN1103053C (zh) | 2003-03-12 |
Family
ID=24833247
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN97121430.1A Expired - Fee Related CN1103053C (zh) | 1996-08-30 | 1997-08-30 | 核谱学系统中用于调整和甄别脉冲形状的方法 |
Country Status (10)
| Country | Link |
|---|---|
| US (1) | US5884234A (zh) |
| CN (1) | CN1103053C (zh) |
| AU (1) | AU727032B2 (zh) |
| BR (1) | BR9704564A (zh) |
| CA (1) | CA2214231C (zh) |
| EG (1) | EG21438A (zh) |
| GB (1) | GB2316744B (zh) |
| ID (1) | ID18147A (zh) |
| MX (1) | MX9706597A (zh) |
| NO (1) | NO973984L (zh) |
Families Citing this family (41)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6369393B1 (en) * | 1998-05-09 | 2002-04-09 | Canberra Industries, Inc. | Digital pulse de-randomization for radiation spectroscopy |
| EP1183553A1 (en) * | 1999-05-27 | 2002-03-06 | Analyser Systems AG | Method and apparatus for multi-parameter digital pulse pile-up rejection |
| US6327549B1 (en) * | 1999-10-25 | 2001-12-04 | Perkinelmer, Inc. | Differential correction method and apparatus |
| US6609075B1 (en) * | 2001-06-04 | 2003-08-19 | William K. Warburton | Method and apparatus for baseline correction in x-ray and nuclear spectroscopy systems |
| US6953937B2 (en) * | 2003-06-26 | 2005-10-11 | Battelle Energy Alliance, Llc | Method and apparatus for the detection of neutrons and gamma rays |
| US7081626B2 (en) * | 2004-06-02 | 2006-07-25 | The Regents Of The University Of California | Apparatus and method for temperature correction and expanded count rate of inorganic scintillation detectors |
| WO2006029475A1 (en) | 2004-09-16 | 2006-03-23 | Southern Innovation International Pty Ltd | Method and apparatus for resolving individual signals in detector output data. |
| FR2879305B1 (fr) * | 2004-12-15 | 2007-06-29 | Commissariat Energie Atomique | Traitement d'un signal representatif de rayonnement |
| FR2888644B1 (fr) * | 2005-07-12 | 2007-09-14 | Commissariat Energie Atomique | Traitement ameliore d'un signal representatif de rayonnement |
| EP1958004A1 (de) * | 2005-12-06 | 2008-08-20 | Fotec Forschungs- und Technologietransfer Gmbh | Verfahren und vorrichtung zum detektieren und zählen elektrisch geladener teilchen sowie verwendung hiefür |
| WO2007146350A2 (en) * | 2006-06-14 | 2007-12-21 | Baker Hughes Incorporated | Pileup rejection |
| US8069197B2 (en) * | 2006-11-15 | 2011-11-29 | Mirmar Sensor, L.L.C. | Method and apparatus for periodic noise suppression in transient signals |
| MX2010001461A (es) | 2007-08-10 | 2010-03-15 | Schlumberger Technology Bv | Escudos reforzados contra neutrones. |
| US7626178B2 (en) * | 2007-12-03 | 2009-12-01 | General Electric Company | Integrated neutron-gamma radiation detector with adaptively selected gamma threshold |
| US9008969B2 (en) * | 2007-12-10 | 2015-04-14 | Schlumberger Technology Corporation | Pulsed x-ray signal processing |
| JP2011519415A (ja) * | 2008-03-31 | 2011-07-07 | サザン イノヴェーション インターナショナル プロプライアトリー リミテッド | 個別信号解像度を用いた放射線画像形成法 |
| CN104020500A (zh) | 2008-03-31 | 2014-09-03 | 南方创新国际股份有限公司 | 用于钻孔测井的方法和装置 |
| US8954300B2 (en) * | 2008-03-31 | 2015-02-10 | Southern Innovation International Pty Ltd. | Screening method and apparatus |
| JP5581333B2 (ja) | 2008-12-18 | 2014-08-27 | サザン イノヴェーション インターナショナル プロプライアトリー リミテッド | 数学的変換を使用することによってパイルアップされたパルスを分離する方法および装置 |
| CN101561507B (zh) * | 2009-05-21 | 2012-07-04 | 清华大学 | 射线探测器的射线能量检测方法 |
| WO2011002452A1 (en) * | 2009-06-30 | 2011-01-06 | Analogic Corporation | Enhanced photon detection for scanner |
| US8492705B2 (en) * | 2009-10-23 | 2013-07-23 | Saint-Gobain Ceramics & Plastics, Inc. | Radiation detector and method of using a radiation detector |
| US8686348B2 (en) | 2011-02-08 | 2014-04-01 | Schlumberger Technology Corporation | High voltage insulating sleeve for nuclear well logging |
| EP2710413B1 (en) * | 2011-05-17 | 2016-12-14 | Services Pétroliers Schlumberger | High throughput pulse height analyzer |
| CN102445446B (zh) * | 2011-09-28 | 2015-09-02 | 沈阳美加生物医药有限公司 | 一种化学发光免疫检测方法及免疫检测装置 |
| JP5859668B2 (ja) | 2011-11-01 | 2016-02-10 | メリル コーポレーション | 中性子スペクトロメータ |
| JP5392437B2 (ja) * | 2011-11-15 | 2014-01-22 | 富士電機株式会社 | パルス処理装置および放射線分析装置 |
| GB2505897A (en) * | 2012-09-13 | 2014-03-19 | Bae Systems Plc | Discriminating Radiation Pulses |
| CN102916683B (zh) * | 2012-10-18 | 2016-09-14 | 成都理工大学 | 一种可调参数核脉冲模拟方法 |
| CN103226204B (zh) * | 2013-04-01 | 2015-07-15 | 中国科学院高能物理研究所 | 探测器信号的测试系统及方法 |
| GB2514418A (en) * | 2013-05-24 | 2014-11-26 | Oxford Instr Nanotechnology Tools Ltd | Method of processing a particle spectrum |
| US20150090888A1 (en) * | 2013-09-27 | 2015-04-02 | Saint-Gobain Ceramics & Plastics, Inc. | Analyzer Device for Compensating a Scintillator and Method of Using the Same |
| CN104849767B (zh) * | 2014-02-18 | 2017-09-01 | 中国石油天然气集团公司 | 一种核测井探测器信号的处理方法和装置 |
| EP2980540A1 (de) * | 2014-08-01 | 2016-02-03 | VEGA Grieshaber KG | Signalabtastung bei einem radiometrischen Messsystem |
| CN104155674B (zh) * | 2014-08-14 | 2017-06-09 | 西北核技术研究所 | 基于数字化波形实时甄别α/γ混合辐射场粒子的方法 |
| CN104536032B (zh) * | 2014-12-23 | 2017-05-17 | 四川大学 | 基于支持向量机的液体闪烁体探测器n/γ射线甄别方法 |
| CN106483551B (zh) * | 2015-08-28 | 2018-11-09 | 成都理工大学 | 一种仿核信号发生器及其工作方法 |
| CN110133708B (zh) * | 2019-04-24 | 2021-03-09 | 中国科学技术大学 | 一种核脉冲信号的测量装置及测量方法 |
| CN111060982A (zh) * | 2019-12-27 | 2020-04-24 | 中广核久源（成都）科技有限公司 | X光机前辐射探测装置 |
| CN113126140B (zh) * | 2019-12-30 | 2024-08-09 | 中核控制系统工程有限公司 | 一种高速数字n/γ波形实时甄别的系统及方法 |
| CN112597923B (zh) * | 2020-12-28 | 2021-09-21 | 成都大学 | 一种基于形态学与优化灰色模型的脉冲堆积校正方法 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4605854A (en) * | 1984-07-16 | 1986-08-12 | Halliburton Company | Measurement of formation porosity using fast neutron spectroscopy |
| DE3635477C2 (de) * | 1986-10-18 | 1994-02-24 | Gsf Forschungszentrum Umwelt | Verfahren und Schaltungsanordnung zur Reduzierung des Zähluntergrundes der Messung mit einer Koinzidenz-Antikoinzidenz-Auswertung und Verwendung des Verfahrens |
| US5067090A (en) * | 1989-05-01 | 1991-11-19 | Schlumberger Technology Corporation | Nuclear spectroscopy method and apparatus for digital pulse height analysis |
| EP0410828A3 (en) * | 1989-06-28 | 1993-03-17 | Schlumberger Limited | Dead time correction and regulation method and apparatus for nuclear spectroscopy |
| US5210423A (en) * | 1989-09-22 | 1993-05-11 | Siemens Gammasonics, Inc. | Method and apparatus for unpiling pulses generated by piled-up scintillation events |
| US5304808A (en) * | 1993-02-08 | 1994-04-19 | The United States Of America As Represented By The United States Department Of Energy | Method and apparatus for data sampling |
| US5347129A (en) * | 1993-04-14 | 1994-09-13 | University Of Missouri-Columbia | System for determining the type of nuclear radiation from detector output pulse shape |
| US5532944A (en) * | 1994-07-28 | 1996-07-02 | Sorrento Electronics, Inc. | Multi-channel analysis system and method using digital signal processing |
-
1996
- 1996-08-30 US US08/705,389 patent/US5884234A/en not_active Expired - Lifetime
-
1997
- 1997-08-28 ID IDP973002A patent/ID18147A/id unknown
- 1997-08-28 GB GB9718234A patent/GB2316744B/en not_active Expired - Fee Related
- 1997-08-28 BR BR9704564A patent/BR9704564A/pt not_active Application Discontinuation
- 1997-08-29 CA CA002214231A patent/CA2214231C/en not_active Expired - Fee Related
- 1997-08-29 NO NO973984A patent/NO973984L/no not_active Application Discontinuation
- 1997-08-29 AU AU36089/97A patent/AU727032B2/en not_active Ceased
- 1997-08-29 MX MX9706597A patent/MX9706597A/es not_active IP Right Cessation
- 1997-08-30 CN CN97121430.1A patent/CN1103053C/zh not_active Expired - Fee Related
- 1997-08-30 EG EG87197A patent/EG21438A/xx active
Also Published As
| Publication number | Publication date |
|---|---|
| EG21438A (en) | 2001-10-31 |
| US5884234A (en) | 1999-03-16 |
| CN1182883A (zh) | 1998-05-27 |
| GB9718234D0 (en) | 1997-11-05 |
| AU3608997A (en) | 1998-03-05 |
| ID18147A (id) | 1998-03-05 |
| BR9704564A (pt) | 1998-11-03 |
| AU727032B2 (en) | 2000-11-30 |
| MX9706597A (es) | 1998-02-28 |
| CA2214231A1 (en) | 1998-02-28 |
| GB2316744A (en) | 1998-03-04 |
| GB2316744B (en) | 1999-08-11 |
| NO973984D0 (no) | 1997-08-29 |
| CA2214231C (en) | 2000-11-21 |
| NO973984L (no) | 1998-03-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN1103053C (zh) | 核谱学系统中用于调整和甄别脉冲形状的方法 | |
| US7683334B2 (en) | Simultaneous beta and gamma spectroscopy | |
| MXPA97006597A (en) | Method for the regulation and discrimination of the form of the impulse in a nucl spectroscopy system | |
| EP0167119B1 (en) | Semiconductor radiation detector | |
| US4482808A (en) | Apparatus for measuring neutrons and gamma rays | |
| US6732059B2 (en) | Ultra-low background gas-filled alpha counter | |
| US4271356A (en) | Method and apparatus for natural gamma ray spectrum analysis using microprocessor stripper | |
| US10048393B2 (en) | Gain stabilization of photomultipliers | |
| EP0117299B1 (en) | Radiation imaging apparatus | |
| US6822235B2 (en) | Environmental radioactivity monitor | |
| JP7239921B2 (ja) | 陽電子消滅特性測定装置 | |
| GB2332512A (en) | Nuclear spectroscopy signal processing system | |
| US7161153B2 (en) | Apparatus and method for detecting α-ray | |
| JP3103047B2 (ja) | 核分光システムにおけるパルス波形調整及び弁別方法 | |
| JP2972936B2 (ja) | シンチレーションパルス波高データの組合わせ方法および装置 | |
| US3896296A (en) | Live time correction system | |
| Martín et al. | Wilcoxon signed-rank-based technique for the pulse-shape analysis of HPGe detectors | |
| GB2482047A (en) | Continuous air monitor | |
| Bell et al. | Wavelength-based neutron/gamma ray discrimination in CLYC | |
| US5324943A (en) | Method for scintillation counting and a scintillation counter with adjustable coincidence resolving time | |
| Johansson | A coincidence scintillation spectrometer | |
| Hennig et al. | Evaluation of phoswich well detectors for radioxenon monitoring | |
| Bristow | IDA. GAMMA RAY SPECTROMETRIC METHODS IN URANIUM EXPLORAnON-AIRBORNE INSTRUMENTAnON | |
| Hubbard et al. | Digital x‐ray processing electronics for fluorescence EXAFS and spectroscopy | |
| WO1991008497A1 (en) | Method for quench determination and compensation in scintillation counting utilizing pulse shape analysis and a scintillation counter |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C19 | Lapse of patent right due to non-payment of the annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |