CN109740279A - 一种特征频率谱库的电可控中子伽马能谱解析方法 - Google Patents

一种特征频率谱库的电可控中子伽马能谱解析方法 Download PDF

Info

Publication number
CN109740279A
CN109740279A CN201910031291.3A CN201910031291A CN109740279A CN 109740279 A CN109740279 A CN 109740279A CN 201910031291 A CN201910031291 A CN 201910031291A CN 109740279 A CN109740279 A CN 109740279A
Authority
CN
China
Prior art keywords
gamma
spectrum
library
function
energy
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.)
Granted
Application number
CN201910031291.3A
Other languages
English (en)
Other versions
CN109740279B (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.)
East China Institute of Technology
Original Assignee
East China Institute of Technology
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 East China Institute of Technology filed Critical East China Institute of Technology
Priority to CN201910031291.3A priority Critical patent/CN109740279B/zh
Publication of CN109740279A publication Critical patent/CN109740279A/zh
Application granted granted Critical
Publication of CN109740279B publication Critical patent/CN109740279B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Radiation (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明公开一种特征频率谱库的电可控中子伽马能谱解析方法，其步骤：探测器响应函数的建立，单元素特征伽马谱库的计算，单元素特征频率谱库的建立，电可控中子伽马能谱最小二乘拟合。本发明的优点在于：采用模拟计算大大降低了单元素特征伽马谱库的实验成本。采用快速傅里叶变化将伽马谱库转化成单元素特征频率谱库，利用幅度与频率的函数提高了元素特征信息的识别度，能够显著提高电可控中子伽马能谱解析精度，克服传统伽马能谱统计性差、各元素特征峰信息不明显、传统能窗谱分析方法准确性差等难题。

Description

一种特征频率谱库的电可控中子伽马能谱解析方法
技术领域
本发明涉及核测井中电可控中子伽马分析技术领域，涉及核技术，尤其是 一种伽马能谱解析方法。
背景技术
核测井中，采用电可控中子源已成为一种趋势，它有其它分析技术所不具 有的很多优点，由于中子不带电，具有极强的穿透能力，尤其是中子与物质发 生核反应时，在非常短的时间内(10-14秒内)发出特征伽马射线，从而来地层 的主要成分含量。近年来，基于电可控中子伽马测井技术的发展。而对伽马能 谱解析是其中关键技术，目前所采取的主要分析方法为能窗法。
核测井中电可控中子伽马能谱是由电可控中子管发射的中子与样品中元 素发生核反应后释放的伽马射线被能谱仪收集、放大及数字化后得到，目前电 可控中子管产生的中子通量相对较低，在用于测井得到的伽马能谱统计性差、 各元素特征峰信息不明显、传统能窗谱分析方法准确性差。因而，现如今急需 一种快速、准确的新型能谱解析方法。
发明内容
本发明的目的在于提供一种特征频率谱库的电可控中子伽马能谱解析方 法，其采用模拟计算大大降低了单元素特征伽马谱库的实验成本，采用快速傅 里叶变化将伽马谱库转化成单元素特征频率谱库，利用幅度与频率的函数提高 了元素特征信息的识别度，能够显著提高电可控中子伽马能谱解析精度，克服 传统伽马能谱统计性差、各元素特征峰信息不明显、传统能窗谱分析方法准确 性差等难题。
本发明的技术方案为：一种特征频率谱库的电可控中子伽马能谱解析方 法，其步骤：
步骤一、探测器响应函数的建立
根据同位素伽马源释放的伽马射线与探测器反应后的能谱，采用半经验公 式法，拟合探测器的能量刻度曲线、探测效率、能量分辨率参数，形成探测器 响应函数；其中探测器响应函数包含5部分：(1)洛伦兹分布特征峰函数，L(n)； (2)畸变的高斯函数，D(n)；(3)阶梯函数，S(n)；(4)本底函数，B(n)；(5)Si逃 逸峰高斯函数，ESi(n)；因此，探测器响应函数DRF(n)和各部分函数的表达式 如下：
DRF(n)＝L(n)+D(n)+S(n)+B(n)+ESi(n)
式中n为多道分析器输出的能谱的道数，HG、HS、HD、HB和HK分别为全 能峰的高斯函数、阶梯函数、高斯畸变函数、本底的短阶梯函数、Si逃逸峰的 高斯函数的幅度值，Γ为洛伦兹宽度，no为全能峰的道址，nT为畸变高斯峰的 道址，nk为Si逃逸峰的道址，σ为全能峰的标准偏差，β为指数部分的斜率， k＝1-8，erfc为反误差函数；HD是全能峰的净峰面积；HS、HB和Hk分别是阶梯 函数、本底和Si逃逸峰的幅度值；
步骤二、单元素特征伽马谱库的计算
通过蒙特卡罗模拟计算软件MCNP得到样品中每个元素不同含量时的特 征伽马谱库，利用粒子信息追踪Ptrac卡，区分元素i的俘获伽马和非弹伽马 能量的强度，元素种类不少于地层中常用的13种元素；在利用俘获伽马能量 的强度Ii与步骤一中得到的探测器响应函数DRF(n)卷积，得到单元素的特征伽 马谱P(n)；即：
式中，P(n)为计算得到的元素i的单元素谱特征伽马谱，Ii为元素i的俘获 伽马能量的强度，DRF为权利1中求得的探测器响应函数，n为多道分析器输 出的能谱的道数,N为待卷积的道数；
采用模拟计算大大降低了单元素特征伽马谱库的实验成本；
步骤三、单元素特征频率谱库的建立
离散的单元素特征伽马谱库P(n)通过快速傅里叶变化，即：
式中，Y(w)为单元素特征伽马谱库P(n)经过傅里叶变化后的频率谱，w 指频率，F为快速傅里叶变化，i指虚数，n为能谱道数；
将计数与能量的关系转化成幅度与频率的函数，建立单元素特征频率谱 库，元素种类不少于13种元素；采用快速傅里叶变化将伽马谱库转化成单元 素特征频率谱库，利用幅度与频率的函数提高了元素特征信息的识别度，能够 显著提高电可控中子伽马能谱解析精度，克服传统伽马能谱统计性差、各元素 特征峰信息不明显、传统能窗谱分析方法准确性差等难题；
步骤四、电可控中子伽马能谱最小二乘拟合
本方法认为实验中总谱是由每个元素的单元素谱累加组成的，因此，
yn是待测样品在第n道址上的总计数，xnj为频率谱库中j元素在第n道的计数， 为第n道上的随机误差，通过上述步骤后，
aj为需要计算得到的元素的含量值，通过χ2最小值得到：
式中σn 2为第n道上计数yn的方差值，通过循环迭代得到最小的χ2值即为频率 谱库最小二乘计算得到的值；利用不同元素特征频率谱的差异，采用最小二乘 拟合实验中的电可控中子伽马能谱，求得样品中多个元素的含量信息，质量百 分含量为100％，实现能谱解析。
所形成的探测器响应函数适用于碘化钠、锗酸铋、溴化镧探测器。
单元素特征伽马谱库包含俘获伽马核非弹散射能谱。
单元素特征频率谱库包含实数部分。
电可控中子伽马能谱最小二乘拟合适用于迭代、差值等方法。
该方法同样适用于瞬发瞬发伽马射线中子活化分析技术的能谱分析。
本发明的优点在于：采用模拟计算大大降低了单元素特征伽马谱库的实 验成本。采用快速傅里叶变化将伽马谱库转化成单元素特征频率谱库，利用 幅度与频率的函数提高了元素特征信息的识别度，能够显著提高电可控中子 伽马能谱解析精度，克服传统伽马能谱统计性差、各元素特征峰信息不明 显、传统能窗谱分析方法准确性差等难题。
具体实施方式
本实施例中，一种特征频率谱库的电可控中子伽马能谱解析方法，该能谱 解析方法是基于样品中各元素的伽马能谱的特征频率谱的差异性进行的。该方 法包括以下步骤：
步骤一、探测器响应函数的建立
本实施例中，探测器选用圣戈班的锗酸铋(BGO)探测器，利用60Co，137Cs 等同位素伽马源释放的1.17MeV、1.33MeV、0.667MeV伽马射线与探测器反 应后的能谱，采用半经验公式法，拟合探测器的能量刻度系数、探测效率、能 量分辨率等关键参数分比为0.0018、0.70、9.6％，形成探测器响应函数。
步骤二、单元素特征伽马谱库的计算
本实施例中，通过模拟软件MCNP得到样品中每个元素不同含量时的特 征伽马谱库，地层种类采用Si、Ca、Mg、H、O、C、N、Fe、Al、S、Cl、Ti、K、U等14种，其含量分别为：5％、28％、1％、24％、10％、5％、6％、2％、 1％、0.5％、2％、10％、5％、0.05％。
步骤三、单元素特征频率谱库的建立
本实施例中，Si、Ca、Mg、H、O、C、N、Fe、Al、S、Cl、Ti、K、U等 14种单元素特征伽马谱库通过快速傅里叶变化，将计数与能量的关系转化成幅 度与频率的函数，建立单元素特征频率谱库。
步骤四、电可控中子伽马能谱最小二乘拟合
本实施例中，利用不同元素特征频率谱的差异，采用最小二乘拟合实验中 的电可控中子伽马能谱，判断公式如下：
式中σn 2为第n道上计数的方差值，通过循环迭代得到最小的χ2值即为最 小二乘计算得到的值。yn是待测样品在第n道址上的总计数，xnj为谱库中j元 素在第n道的计数，为第n道上的随机误差，i、m是元素，aj为需要计算 得到的元素的含量值。通过上面的公式计算得到样品中Si、Ca、Mg、H、O、 C、N、Fe、Al、S、Cl、Ti、K、U等14种元素的含量分别为：5.12％、28.27％、 0.91％、23.81％、10.05％、5.06％、6.04％、2.11％、1.09％、0.49％、1.89％、10.07％、 5.06％、0.049％，其与理论值的偏差如表1所示，从表1中可以看出通过该方 法计算得到的元素含量与理论值偏差小于0.27％，实现能谱准确解析。
表1特征频率谱库最小二乘法的计算结果与理论值比较
以上所述，仅是本发明的叫价实施例，并非对本发明作任何限制，根据本 发明技术实质可以有各种更改和变化。凡是在本发明的精神和原则之内，所作 的任何修改、等同提花、改进等，均应包含在本发明的保护范围之内。

Claims (6)

1.一种特征频率谱库的电可控中子伽马能谱解析方法，其步骤：
步骤一、探测器响应函数的建立
根据同位素伽马源释放的伽马射线与探测器反应后的能谱，采用半经验公式法，拟合探测器的能量刻度曲线、探测效率、能量分辨率参数，形成探测器响应函数；其中探测器响应函数包含5部分：(1)洛伦兹分布特征峰函数，L(n)；(2)畸变的高斯函数，D(n)；(3)阶梯函数，S(n)；(4)本底函数，B(n)；(5)Si逃逸峰高斯函数，ESi(n)；因此，探测器响应函数DRF(n)和各部分函数的表达式如下：
DRF(n)＝L(n)+D(n)+S(n)+B(n)+ESi(n)
式中n为多道分析器输出的能谱的道数，HG、HS、HD、HB和HK分别为全能峰的高斯函数、阶梯函数、高斯畸变函数、本底的短阶梯函数、Si逃逸峰的高斯函数的幅度值，Γ为洛伦兹宽度，no为全能峰的道址，nT为畸变高斯峰的道址，nk为Si逃逸峰的道址，σ为全能峰的标准偏差，β为指数部分的斜率，k＝1-8，erfc为反误差函数；HD是全能峰的净峰面积；HS、HB和Hk分别是阶梯函数、本底和Si逃逸峰的幅度值；
步骤二、单元素特征伽马谱库的计算
通过蒙特卡罗模拟计算软件MCNP得到样品中每个元素不同含量时的特征伽马谱库，利用粒子信息追踪Ptrac卡，区分元素i的俘获伽马和非弹伽马能量的强度，元素种类不少于地层中常用的13种元素；在利用俘获伽马能量的强度Ii与步骤一中得到的探测器响应函数DRF(n)卷积，得到单元素的特征伽马谱P(n)；即：
式中，P(n)为计算得到的元素i的单元素谱特征伽马谱，Ii为元素i的俘获伽马能量的强度，DRF为权利1中求得的探测器响应函数，n为多道分析器输出的能谱的道数，N为待卷积的道数；
采用模拟计算大大降低了单元素特征伽马谱库的实验成本；
步骤三、单元素特征频率谱库的建立
离散的单元素特征伽马谱库P(n)通过快速傅里叶变化，即：
式中，Y(w)为单元素特征伽马谱库P(n)经过傅里叶变化后的频率谱，w指频率，F为快速傅里叶变化，i指虚数，n为能谱道数；
将计数与能量的关系转化成幅度与频率的函数，建立单元素特征频率谱库，元素种类不少于13种元素；采用快速傅里叶变化将伽马谱库转化成单元素特征频率谱库，利用幅度与频率的函数提高了元素特征信息的识别度，能够显著提高电可控中子伽马能谱解析精度，克服传统伽马能谱统计性差、各元素特征峰信息不明显、传统能窗谱分析方法准确性差等难题；
步骤四、电可控中子伽马能谱最小二乘拟合
本方法认为实验中总谱是由每个元素的单元素谱累加组成的，因此，
yn是待测样品在第n道址上的总计数，xnj为频率谱库中j元素在第n道的计数，为第n道上的随机误差，通过上述步骤后，
aj为需要计算得到的元素的含量值，通过χ2最小值得到：
式中σn 2为第n道上计数yn的方差值，通过循环迭代得到最小的χ2值即为频率谱库最小二乘计算得到的值；利用不同元素特征频率谱的差异，采用最小二乘拟合实验中的电可控中子伽马能谱，求得样品中多个元素的含量信息，质量百分含量为100％，实现能谱解析。
2.根据权利要求1所述特征频率谱库的电可控中子伽马能谱解析方法，其特征在于：所形成的探测器响应函数适用于碘化钠、锗酸铋、溴化镧探测器。
3.根据权利要求1所述特征频率谱库的电可控中子伽马能谱解析方法，其特征在于：单元素特征伽马谱库包含俘获伽马核非弹散射能谱。
4.根据权利要求1所述特征频率谱库的电可控中子伽马能谱解析方法，其特征在于：单元素特征频率谱库包含实数部分。
5.根据权利要求1所述特征频率谱库的电可控中子伽马能谱解析方法，其特征在于：电可控中子伽马能谱最小二乘拟合适用于迭代、差值等方法。
6.根据权利要求1所述特征频率谱库的电可控中子伽马能谱解析方法，其特征在于：该方法同样适用于瞬发瞬发伽马射线中子活化分析技术的能谱分析。
CN201910031291.3A 2019-01-14 2019-01-14 一种特征频率谱库的电可控中子伽马能谱解析方法 Active CN109740279B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910031291.3A CN109740279B (zh) 2019-01-14 2019-01-14 一种特征频率谱库的电可控中子伽马能谱解析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910031291.3A CN109740279B (zh) 2019-01-14 2019-01-14 一种特征频率谱库的电可控中子伽马能谱解析方法

Publications (2)

Publication Number Publication Date
CN109740279A true CN109740279A (zh) 2019-05-10
CN109740279B CN109740279B (zh) 2022-06-14

Family

ID=66364659

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910031291.3A Active CN109740279B (zh) 2019-01-14 2019-01-14 一种特征频率谱库的电可控中子伽马能谱解析方法

Country Status (1)

Country Link
CN (1) CN109740279B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112258432A (zh) * 2020-10-23 2021-01-22 成都大学 基于复合型结构元素数学形态学的中子-伽马甄别方法
CN112326708A (zh) * 2020-10-20 2021-02-05 汪金龙 基于瞬发伽马光谱的人体组织密度和元素组成重建算法
CN113671588A (zh) * 2021-07-15 2021-11-19 中国石油大学(北京) 元素测井解谱方法、装置、电子设备及存储介质
CN115508879A (zh) * 2022-11-21 2022-12-23 山东省科学院海洋仪器仪表研究所 一种海水放射性检测的伽马能谱全谱分析方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4503328A (en) * 1981-09-14 1985-03-05 Halliburton Company Evaluation of the interaction with radiant energy of substances traversed by a borehole
US20080230707A1 (en) * 2007-03-23 2008-09-25 Verista Imaging, Inc. High resolution near-field imaging method and apparatus
CN102449505A (zh) * 2009-04-30 2012-05-09 R·P·塔里雅克汗 测定辐射方向性的组合物及方法
CN103135125A (zh) * 2011-12-01 2013-06-05 中国辐射防护研究院 一种过滤x射线参考辐射的能谱分析方法
CN103377746A (zh) * 2012-04-14 2013-10-30 中国科学技术大学 实现显微镜系统超分辨成像的方法
CN104656117A (zh) * 2015-01-29 2015-05-27 中国人民解放军国防科学技术大学 一种利用全耗尽型硅探测器甄别质子和α粒子的方法
CN107907901A (zh) * 2017-11-06 2018-04-13 中广核工程有限公司 核电站中子、伽马能谱与剂量的测量方法和系统
CN108333637A (zh) * 2017-12-14 2018-07-27 中国矿业大学 一种提高元素测井技术确定元素含量准确度的方法
CN108825219A (zh) * 2018-05-25 2018-11-16 东华理工大学 融合自然γ能谱与中子时间谱的铀矿测井刻度参数求法
CN108918565A (zh) * 2018-05-11 2018-11-30 南京航空航天大学 一种基于瞬发伽马射线中子活化分析技术的样品元素分布测量装置及方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4503328A (en) * 1981-09-14 1985-03-05 Halliburton Company Evaluation of the interaction with radiant energy of substances traversed by a borehole
US20080230707A1 (en) * 2007-03-23 2008-09-25 Verista Imaging, Inc. High resolution near-field imaging method and apparatus
CN102449505A (zh) * 2009-04-30 2012-05-09 R·P·塔里雅克汗 测定辐射方向性的组合物及方法
CN103135125A (zh) * 2011-12-01 2013-06-05 中国辐射防护研究院 一种过滤x射线参考辐射的能谱分析方法
CN103377746A (zh) * 2012-04-14 2013-10-30 中国科学技术大学 实现显微镜系统超分辨成像的方法
CN104656117A (zh) * 2015-01-29 2015-05-27 中国人民解放军国防科学技术大学 一种利用全耗尽型硅探测器甄别质子和α粒子的方法
CN107907901A (zh) * 2017-11-06 2018-04-13 中广核工程有限公司 核电站中子、伽马能谱与剂量的测量方法和系统
CN108333637A (zh) * 2017-12-14 2018-07-27 中国矿业大学 一种提高元素测井技术确定元素含量准确度的方法
CN108918565A (zh) * 2018-05-11 2018-11-30 南京航空航天大学 一种基于瞬发伽马射线中子活化分析技术的样品元素分布测量装置及方法
CN108825219A (zh) * 2018-05-25 2018-11-16 东华理工大学 融合自然γ能谱与中子时间谱的铀矿测井刻度参数求法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张鹏飞 等: "傅里叶变换法扣除中子活化分析能谱本底", 《核电子学与探测技术》 *
李飞: "基于蒙特卡罗模拟和人工神经算法的伽马射线能谱数据处理技术", 《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅱ辑》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112326708A (zh) * 2020-10-20 2021-02-05 汪金龙 基于瞬发伽马光谱的人体组织密度和元素组成重建算法
CN112326708B (zh) * 2020-10-20 2023-10-24 汪金龙 基于瞬发伽马光谱的人体组织密度和元素组成重建算法
CN112258432A (zh) * 2020-10-23 2021-01-22 成都大学 基于复合型结构元素数学形态学的中子-伽马甄别方法
CN112258432B (zh) * 2020-10-23 2021-09-21 成都大学 基于复合型结构元素数学形态学的中子-伽马甄别方法
CN113671588A (zh) * 2021-07-15 2021-11-19 中国石油大学(北京) 元素测井解谱方法、装置、电子设备及存储介质
CN113671588B (zh) * 2021-07-15 2022-10-14 中国石油大学(北京) 元素测井解谱方法、装置、电子设备及存储介质
CN115508879A (zh) * 2022-11-21 2022-12-23 山东省科学院海洋仪器仪表研究所 一种海水放射性检测的伽马能谱全谱分析方法
CN115508879B (zh) * 2022-11-21 2023-02-24 山东省科学院海洋仪器仪表研究所 一种海水放射性检测的伽马能谱全谱分析方法

Also Published As

Publication number Publication date
CN109740279B (zh) 2022-06-14

Similar Documents

Publication Publication Date Title
CN109740279A (zh) 一种特征频率谱库的电可控中子伽马能谱解析方法
JP6255552B2 (ja) 中性子検出器
CN107526113B (zh) 俘获伽马能谱的确定方法和装置
CN103135125A (zh) 一种过滤x射线参考辐射的能谱分析方法
Mijatović et al. Study of the cross section determination with the PRISMA spectrometer: The 40 Ar+ 208 Pb case
Schrack et al. Neutral meson photoproduction from complex nuclei
Pechousek et al. Software emulator of nuclear pulse generation with different pulse shapes and pile-up
Boardman et al. A gamma-ray identification algorithm based on Fisher linear discriminant analysis
Lawrence Neutron Spectrum Unfolding with Organic Scintillators for Arms-control Verification.
CN107238856B (zh) 一种高通量氘-氚中子发生器中子平均能量的确定方法
Bushnin et al. J/ψ particle production in π−-N interactions at 27 and 40 GeV/c
Prestwich et al. Gamma-Ray Spectra of Co 60 and Mn 56 Following Resonance-Neutron Capture in Co 59 and Mn 55
Luo et al. Application research on neutron-gamma discrimination based on BC501A liquid scintillator
Schneid et al. Study of Inelastic Scattering of Protons from Sn 116, Sn 122, and Sn 124 at Isobaric Analog Resonances
Kim et al. Monte Carlo simulation for the analysis of various solid samples using handheld X-ray fluorescence spectrometer and evaluation of the effect by environmental interferences
Chichester et al. Measurement of the neutron spectrum of a DD electronic neutron generator
CN113340978B (zh) 非天然同位素丰度的元素浓度测量方法及系统
Wen et al. Evaluation of a modified HPGe preamplifier for high-rate spectroscopy measurements in a pulsed photonuclear environment
Foley Short-lived Photofission Product Yields and Analytical Methods for Nuclear Forensic Application
Watson et al. Distribution of K X Rays as a Function of Mass and Atomic Number in the Spontaneous Fission of Cf 252
Papp et al. Radiation Measurements without Posing the Inversion Problem
Kulisek et al. Lead Slowing-Down Spectrometry Time Spectral Analysis for Spent Fuel Assay: FY12 Status Report
Yi et al. Discrimination of piled-up neutron-gamma pulses using charge comparison method and neural network for CLYC detectors
Gadey Gamma Spectrum Enhancement in a Metal Loaded Plastic Scintillator using PSD
Kelley et al. Measurement of the Fast Neutron Response for $^ 4 {\rm He} $ Scintillation Detectors Using a Coincidence Scattering Method

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