CN111638546B - 一种海水放射性核素能谱峰重叠的判断方法 - Google Patents

一种海水放射性核素能谱峰重叠的判断方法 Download PDF

Info

Publication number
CN111638546B
CN111638546B CN202010405920.7A CN202010405920A CN111638546B CN 111638546 B CN111638546 B CN 111638546B CN 202010405920 A CN202010405920 A CN 202010405920A CN 111638546 B CN111638546 B CN 111638546B
Authority
CN
China
Prior art keywords
peak
channel
peaks
energy spectrum
nuclide
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
CN202010405920.7A
Other languages
English (en)
Other versions
CN111638546A (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.)
Institute of Oceanographic Instrumentation Shandong Academy of Sciences
Original Assignee
Institute of Oceanographic Instrumentation Shandong Academy of Sciences
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 Institute of Oceanographic Instrumentation Shandong Academy of Sciences filed Critical Institute of Oceanographic Instrumentation Shandong Academy of Sciences
Priority to CN202010405920.7A priority Critical patent/CN111638546B/zh
Publication of CN111638546A publication Critical patent/CN111638546A/zh
Application granted granted Critical
Publication of CN111638546B publication Critical patent/CN111638546B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/36Measuring spectral distribution of X-rays or of nuclear radiation spectrometry

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

本发明属于海水检测技术领域,涉及一种海水放射性核素能谱峰重叠的判断方法,包括:(1)采用两种不同的核素对海水进行放射性检测,得到两个核素的峰;(2)针对两个核素的峰,分别建立拟合函数;(3)从左至右,遍历两个核素峰之间的通道;(4)根据当前遍历的通道Xi,分别用两个拟合函数计算每个通道对应的拟合信号值
Figure DDA0002491260460000013
(5)从左到右遍历,如果
Figure DDA0002491260460000011
则结束判断;(6)如果
Figure DDA0002491260460000012
进一步判断:若
Figure DDA0002491260460000015
继续遍历;若
Figure DDA0002491260460000014
则记录该通道chi;(7)若记录的通道满足:(chi‑ch1)≤2.5×bgk1,且(ch2‑chi)≤2.5×bgk2,则判定两种放射性核素的能谱峰发生重叠。本发明的方法不受累积时间和海洋环境干扰的限制,能够识别出两个海水放射性核素能谱的峰是否发生重叠,可以提高计算的准确度。

Description

一种海水放射性核素能谱峰重叠的判断方法
技术领域
本发明属于海水检测技术领域,涉及一种海水放射性核素能谱峰重叠的判断方法。
背景技术
在海洋放射性物质综合测量过程中,如果存在相应的放射性物质,在所对应的能量区间会出现相应的峰。由于检测到放射性核素信号的峰发生无规律的漂移、重叠,所以二种海水放射性核素的能谱峰,可能存在能谱峰的重叠,影响海水放射性核素的检测结果。检测海水放射性核素能谱峰的判断方法,是海洋检测放射性核素测量所必须的,是目前国内海洋放射性测量发展的核心。现有的海水放射性环境监测方法通常是以假定海水中待测放射性物质所对应的检测信号是平稳的或是时不变的为前提,并不适用于复杂的海洋环境。
但是,实际的海洋现场检测环境是复杂多变的,干扰因素很多。在实际作业过程中发现,海洋环境的变化对海洋物质的测量存在彼此干扰,其现象是检测到的能谱会出现无规律的变化,为得到真正海水放射性核素的检测结果,及时发现海水放射性核素的能量范围。
基于以上原因,如果用两种或两种以上不同的核素进行海水放射性检测,由于在放射性核素检测过程中,会出现很多的峰,有假峰,有叠加峰,也有干扰数据引起的误差峰,漂移峰等,两种核素的峰可能会因为干扰素叠加在一起,影响检测结果的判断。因此,现场测量必须找到一种方法来准确判断不同海水放射性核素能谱峰是否发生重叠。
发明内容
为解决上述技术问题,本发明提供了一种海水放射性核素能谱峰重叠的判断方法,能够识别出二个海水放射性核素能谱的峰,是否发生重叠,可以提高计算的准确度,以达到提高检测放射性核素的效率和准确度的目的。
为达到上述目的,本发明的技术方案如下:
一种海水放射性核素能谱峰重叠的判断方法,包括以下步骤:
(1)采用两种不同的核素对海水进行放射性检测,得到两个核素的峰;
(2)针对两个核素的峰,分别建立拟合函数;
(3)从左至右,遍历两个核素峰之间的通道;
(4)根据当前遍历的通道Xi,分别用两个拟合函数计算每个通道对应的拟合信号值
Figure BDA0002491260440000011
Figure BDA0002491260440000012
(5)从左到右遍历,如果
Figure BDA0002491260440000021
则结束判断;
(6)如果
Figure BDA0002491260440000022
进一步判断:若
Figure BDA0002491260440000023
继续遍历;若
Figure BDA0002491260440000024
则记录该通道chi
(7)若记录的通道满足:(chi-ch1)≤2.5×bgk1,且(ch2-chi)≤2.5×bgk2,则判定两种放射性核素的能谱峰发生重叠,重叠的通道为该遍历点的通道;其中,ch1为第一个核素峰的通道,ch2为第二个核素峰的通道;kgk1为第一个核素峰的半高宽,kgk2为第二个核素峰的半高宽。
作为本发明的一种优选方式,步骤(1)中,判断两个核素的峰是否发生漂移,若发生漂移,则对两个放射性核素的峰进行重新标定,以标定后的峰为准。
进一步优选地,所述步骤(2),采用高斯拟合函数分别对两种放射性核素的能谱峰,建立各自的拟合函数:
Figure BDA0002491260440000025
其中,Ymax为能谱峰的峰信号值,Xmax为能谱峰的峰通道,S/2为能谱峰拟合半高宽的平方值,Xi为拟合曲线的通道值,Yi为拟合曲线通道对应的信号值。
进一步优选地,所述步骤(4)中,左边的峰为第一个峰,右边的峰为第二个峰;第一个峰的峰通道小于第二个峰的峰通道。
通过上述技术方案,本发明提供了一种海水放射性核素能谱峰重叠的判断方法,能够识别出二个海水放射性核素能谱的峰,是否发生重叠,可以提高计算的准确度,以达到提高检测放射性核素的效率和准确度的目的。
附图说明
图1为本发明实施实例所公开的海水放射性核素能谱峰重叠的判断方法流程示意图;
图2为本发明实施实例所公开的所有通道的能谱数据图;
图3为本发明实施实例所公开的放大的能谱曲线;
图4为本发明实施实例所公开的拟合后的能谱曲线;
图5为本发明实施实例放大拟合后两条拟合曲线交汇点。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。
本发明提供了一种海水放射性核素干扰频率的判断方法,如图1所示,具体流程如下:
S101、采用两种不同的核素对海水进行放射性检测,得到两个核素的峰,对于该两个峰,分别判断是否发生漂移,若没有发生漂移,则进行下一步;若发生漂移,则找出漂移峰,对原峰进行重新标定,以标定后的峰为准,再进行下一步。
其中,漂移峰的判断方法为:
(1)、从放射性核素的峰的预设通道,向左右分别各遍历50个通道,设为通道CHn;预设通道为仪器标定的通道ch,通道CHn与预设通道ch的差值分别为-50,-49,-48,-47,……,+49,+50;遍历后的通道CHn为ch-50,ch-49,ch-48,ch-47,……,ch+49,ch+50。
(2)、判断该通道CHn的信号值是否与放射性核素峰的预设通道信号值相差5%以内;预设通道信号值为仪器标定的信号值v;如果通道CHn的信号信号V满足:v*0.95≤V≤v*1.05,就进入下一步骤。标定为累积3小时的数据,在判断的时候也必须要用累积3小时的数据。
(3)、对于上述步骤中,满足v*0.95≤V≤v*1.05的通道CHn,假设该通道为漂移后的峰,分别向左和向右查找峰值一半的点。通道CHn的信号值为V,从CHn向左查找第一个信号值小于半高信号值V’=V/2的点的通道为
Figure BDA0002491260440000031
从CHn向右查找第一个信号值小于半高信号值V’=V/2的点的通道为
Figure BDA0002491260440000032
(4)、分别计算左点、右点与假设的漂移后的峰通道所相差的通道数,如果该差值在预设放射性核素峰的半高宽的±5%之内,记录该通道。
预设放射性核素峰的半高宽为仪器标定的预设放射性核素峰的半高宽bgk,若
Figure BDA0002491260440000033
Figure BDA0002491260440000034
则记录该通道CHn
(5)左、右半高宽差值最小的通道,为漂移后的峰通道,如果在前后各50个点所计算的左右半高宽差值最小的通道不只一个,就取与预设通道最接近的,判为漂移后的峰通道。放射性核素峰发生了漂移,漂移后通道对应的峰为真实放射性核素峰。
在本实施例中,经过判断确定有两个核素的峰通道分别为645、732通道;峰信号值分别为1218、940;半高宽分别为16、17;所有通道的能谱曲线,如图2、图3所示。
S102、采用高斯拟合函数分别对两个放射性核素的能谱峰,建立各自的拟合函数,如图4所示;
根据能谱峰的通道,信号值,和半高宽,计算出能谱峰的高斯函数。高斯函数如下:
Figure BDA0002491260440000041
其中,Ymax为能谱峰的峰信号值,Xmax为能谱峰的峰通道,S/2为能谱峰拟合半高宽的平方值,Xi为拟合曲线的通道值,Yi为拟合曲线通道对应的信号值。
在本实施例中,公式如下:
Figure BDA0002491260440000042
S1=2×162
Figure BDA0002491260440000043
S2=2×172
S103、从第一个峰的通道开始遍历,每一个通道,直到第二个峰的通道;第一个峰的为左边的峰ch1,第二个峰为右边的峰ch2,也就是第一个峰的峰通道小于第二个峰的峰通道;遍历两个峰通道之间的每个通道;
在本实施例中,ch1=645,ch2=732。遍历645通道到732通道之间的每一个通道。
S104、根据当前遍历的通道,分别用两个拟合的函数计算对应的拟合信号值;遍历的通道为Xi,拟合信号值为
Figure BDA0002491260440000044
S105、如果
Figure BDA0002491260440000045
Figure BDA0002491260440000046
继续遍历;如果
Figure BDA0002491260440000047
则退出遍历,判定两种放射性核素的能谱峰不发生重叠;
S106、如果
Figure BDA0002491260440000048
Figure BDA0002491260440000049
记录该通道chi
在本实施实例中,如图5所示,在通道688处,
Figure BDA00024912604400000410
满足
Figure BDA00024912604400000411
Figure BDA00024912604400000412
Figure BDA00024912604400000413
记录该通道ch=688。
S107、若该通道满足:(chi-ch1)≤2.5×bgk1,且(ch2-chi)≤2.5×bgk2,则判定两个放射性核素的能谱峰发生重叠,重叠的通道为该遍历点的通道。
否则这两个海水放射性核素能谱峰没有发生重叠。其中,ch1为第一个核素峰的通道,ch为第二个核素峰的通道;bgk1为第一个核素峰的半高宽,bgk2为第二个核素峰的半高宽。
在本实施实例中,ch-ch1=688-645=43;2.5×bgk1=2.5×16=40;
ch2-ch=732-688=44;2.5×bgk2=2.5×17=42.5;
因为43>40,而且44>42.5;也就是该通道不在第一个峰2.5倍的半高宽之内,也不在第二个峰2.5倍半高宽之内,所以判断这两个核素的能谱峰没有发生重叠。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (2)

1.一种海水放射性核素能谱峰重叠的判断方法,其特征在于:包括以下步骤:
(1)采用两种不同的核素对海水进行放射性检测,得到两个核素的峰;
(2)针对两个核素的峰,采用高斯拟合函数分别对两种放射性核素的能谱峰,建立各自的拟合函数:
Figure FDA0003656338600000011
Figure FDA0003656338600000012
其中,
Figure FDA0003656338600000013
为第一个核素能谱峰的峰信号值,
Figure FDA0003656338600000014
为第二个核素能谱峰的峰信号值;
Figure FDA0003656338600000015
分别为第一个核素能谱峰、第二个核素能谱峰的峰通道;S1、S2分别为第一个核素能谱峰、第二个核素能谱峰拟合半高宽的平方值×2;Xi为拟合曲线的通道值;
Figure FDA0003656338600000016
为采用第一个核素峰拟合曲线计算的每个通道对应拟合信号值;
Figure FDA0003656338600000017
为采用第二个核素峰拟合曲线计算的每个通道对应拟合信号值;
(3)从左至右,遍历两个核素峰之间的通道;左边的峰为第一个峰,右边的峰为第二个峰;第一个峰的峰通道小于第二个峰的峰通道;
(4)根据当前遍历的通道Xi,分别用两个拟合函数计算每个通道对应的拟合信号值
Figure FDA0003656338600000018
(5)从左到右遍历,如果
Figure FDA0003656338600000019
则结束判断;
(6)如果
Figure FDA00036563386000000110
进一步判断:若
Figure FDA00036563386000000111
继续遍历;若
Figure FDA00036563386000000112
则记录该通道chi
(7)若记录的通道满足:(chi-ch1)≤2.5×bgk1,且(ch2-chi)≤2.5×bgk2,则判定两种放射性核素的能谱峰发生重叠,重叠的通道为该遍历点的通道;其中,ch1为第一个核素峰的通道,ch2为第二个核素峰的通道;bgk1为第一个核素峰的半高宽,bgk2为第二个核素峰的半高宽。
2.根据权利要求1所述的一种海水放射性核素能谱峰重叠的判断方法,其特征在于:所述步骤(1)中,判断两个核素的峰是否发生漂移,若发生漂移,则对两个放射性核素的峰进行重新标定,以标定后的峰为准。
CN202010405920.7A 2020-05-14 2020-05-14 一种海水放射性核素能谱峰重叠的判断方法 Active CN111638546B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010405920.7A CN111638546B (zh) 2020-05-14 2020-05-14 一种海水放射性核素能谱峰重叠的判断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010405920.7A CN111638546B (zh) 2020-05-14 2020-05-14 一种海水放射性核素能谱峰重叠的判断方法

Publications (2)

Publication Number Publication Date
CN111638546A CN111638546A (zh) 2020-09-08
CN111638546B true CN111638546B (zh) 2022-07-22

Family

ID=72328480

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010405920.7A Active CN111638546B (zh) 2020-05-14 2020-05-14 一种海水放射性核素能谱峰重叠的判断方法

Country Status (1)

Country Link
CN (1) CN111638546B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09304542A (ja) * 1996-05-20 1997-11-28 Aloka Co Ltd 放射線測定装置
CN102298153A (zh) * 2010-06-23 2011-12-28 成都理工大学 放射性测量中多重谱峰的分解方法
CN105203565A (zh) * 2014-06-11 2015-12-30 成都理工大学 一种能谱重叠峰解析方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI495896B (zh) * 2011-04-19 2015-08-11 Iner Aec Executive Yuan 碘化鈉偵檢器量測所得能譜之活度計算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09304542A (ja) * 1996-05-20 1997-11-28 Aloka Co Ltd 放射線測定装置
CN102298153A (zh) * 2010-06-23 2011-12-28 成都理工大学 放射性测量中多重谱峰的分解方法
CN105203565A (zh) * 2014-06-11 2015-12-30 成都理工大学 一种能谱重叠峰解析方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于模糊逻辑的γ能谱核素识别;王一鸣等;《清华大学学报(自然科学版)》;20121215(第12期);全文 *

Also Published As

Publication number Publication date
CN111638546A (zh) 2020-09-08

Similar Documents

Publication Publication Date Title
US7462837B2 (en) Nuclide identifier system
CN111046833B (zh) 一种海水放射性核素检测的二次寻峰方法
CN108375783B (zh) 一种海水放射性检测的自动寻峰方法
CN109696702B (zh) 一种海水放射性核素k40检测的重叠峰判断方法
CN109901216B (zh) 一种检测海水放射性核素的寻峰方法
CN111638546B (zh) 一种海水放射性核素能谱峰重叠的判断方法
CN111666532B (zh) 一种海水放射性核素峰重叠的分峰计算方法
JPS6375573A (ja) 波形識別方法
CN109669205B (zh) 一种海水放射性核素k40元素的寻峰方法
CN111443371B (zh) 一种海水放射性核素峰漂移的判断方法
CN110082808B (zh) 一种基于复杂背景下核脉冲信号快速探测及识别方法
CN109765601B (zh) 一种海水中放射性核素k40元素的计数率的计算方法
CN111060950B (zh) 一种海水放射性核素检测的重叠峰去除方法
CN109525308B (zh) 一种光模块的检测方法、装置、设备及存储介质
CN113092916B (zh) 一种孔缝端口的电磁环境适应性边界场景快速推演方法
CN113009544A (zh) 一种海水放射性核素峰范围的判断方法
CA3194994A1 (en) Multi-sector data correction method and system for gamma imaging while drilling
CN111045069B (zh) 一种海水放射性核素检测的数据修正方法
US4075486A (en) Method and device for diminishing the background effect in a proportional counter
TW550554B (en) BCA evaluation device and BCA evaluation method for rewritable DVDs
CN112965098A (zh) 一种海水放射性核素检测的能量刻度方法
CN117871998A (zh) 一种串口设备异常波形检测的方法及系统
WO2021133224A1 (ru) Способ диагностирования стабильности работы устройства с коронным счетчиком
JP3494912B2 (ja) 周波数誤差検出装置
CN116087716A (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