首页 / 专利库 / 分销网络和设备 / 发电厂 / 核能电厂 / 核反应堆 / 堆芯 / 燃料组件 / 一种Pu-242核素相对含量的计算方法

一种Pu-242核素相对含量的计算方法

阅读:556发布:2020-05-08

专利汇可以提供一种Pu-242核素相对含量的计算方法专利检索,专利查询,专利分析的服务。并且本 发明 属于核 辐射 与探测技术领域,涉及一种Pu-242核素相对含量的计算方法。该方法利用γ谱仪测量238Pu、239Pu、240Pu、241Pu的特征γ射线,利用239Pu的特征γ射线,获得不同 能量 γ射线的相对探测效率曲线;根据238Pu、239Pu、240Pu、241Pu的特征γ射线的计数,结合γ射线探测效率随能量的变化,计算238Pu、240Pu、241Pu相对239Pu的含量,然后计算出242Pu相对239Pu的含量。本发明利用现有的探测仪器仪表,解决了无法直接识别242Pu、探测其含量的问题, 算法 中已考虑多种组件、多种堆型的多种辐照历史、各种燃耗深度、各个乏 燃料 冷却时间下的情况,是有效和准确的。计算的242Pu含量的相对误差基本在1%左右,具有很高的 精度 。,下面是一种Pu-242核素相对含量的计算方法专利的具体信息内容。

1.一种Pu-242核素相对含量的计算方法,其特征在于,包括:
利用γ谱仪测量238Pu、239Pu、240Pu、241Pu的特征γ射线,利用239Pu的特征γ射线,获得不同能量γ射线的相对探测效率曲线;
根据238Pu、239Pu、240Pu、241Pu的特征γ射线的计数,结合γ射线探测效率随能量的变化,计算238Pu、240Pu、241Pu相对239Pu的含量,即238Pu/239Pu、240Pu/239Pu、241Pu/239Pu;
根据如下公式计算出242Pu相对239Pu的含量:
y=a1+a2*x1+a3*x2+a4*x3+a5*x12+a6*x22+a7*x32+a8*x1*x2+a9*x1*x3+a10*x2*x3+
3 3 3 2 2 2 2 2
a11*x1+a12*x2+a13*x3+a14*x1*x2+a15*x1*x3+a16*x2*x1+a17*x2*x3+a18*x3*x1+a19*x32*x2+a20*x1*x2*x3
公式中,
x1为238Pu/239Pu相对含量;
240 239
x2为 Pu/ Pu相对含量;
x3为241Pu/239Pu相对含量;
y为242Pu/239Pu相对含量;
a1~a20为待定系数。
2.如权利要求1所述的Pu-242核素相对含量的计算方法,其特征在于,238Pu、240Pu、241Pu相对239Pu的含量根据如下公式计算:
其中, 为238Pu或240Pu或241Pu相对于239Pu的含量;Ax,k、A239,j分别为238Pu或240Pu或
241Pu的特征γ射线k、239Pu的特征γ射线j的全能峰计数;εj、εk分别为特征γ射线j、k的相对探测效率;px,k、p239,j是238Pu或240Pu或241Pu的特征γ射线k、239Pu的特征γ射线j的强度;λx、λ239是238Pu或240Pu或241Pu、239Pu的衰变常数。
3.如权利要求1所述的Pu-242核素相对含量的计算方法,其特征在于,242Pu相对239Pu含量的计算公式是基于典型核电站,通过汇总不同燃料组件、不同燃耗深度下、不同冷却时间下钚同位素的含量原始数据,统计分析后得到的。
4.如权利要求3所述的Pu-242核素相对含量的计算方法,其特征在于,242Pu相对239Pu含量的计算公式中a1~a20待定系数是基于典型核电站不同燃料组件、不同燃耗深度下、不同冷却时间下钚同位素的含量原始数据,利用LM最小二乘法非线性拟合最优化算法进行拟合计算分析后得到的。
5.如权利要求4所述的Pu-242核素相对含量的计算方法,其特征在于,242Pu相对239Pu含量的计算公式中a1~a20待定系数,根据Pu已知核素相对含量的范围进行分组拟合计算分析得到。
6.如权利要求5所述的Pu-242核素相对含量的计算方法,其特征在于,a1~a20待定系数是根据238Pu/239Pu核素相对含量进行分组拟合计算分析得到。

说明书全文

一种Pu-242核素相对含量的计算方法

技术领域

[0001] 本发明属于核辐射与探测技术领域,涉及一种与燃料组件辐照历史无关的乏燃料后处理中242Pu核素相对含量的计算方法。

背景技术

[0002] 乏燃料组件中钚同位素主要包括238Pu、239Pu、240Pu、241Pu和242Pu。通常,钚核素的识别主要通过测量中低能区特征伽射线实现的。特征伽马射线法要求核素有能明显区别于其它核素的、特定能量的、可分辨的伽马射线。核素238Pu、239Pu、240Pu、241Pu都有能够用于核素识别的特征伽马射线,为了缓解对探测器定量标定的需求,降低对能量分辨率以及电子学系统性能的过高要求,通常采用相对探测效率法获得同位素相对含量(组成比例),再结合中子计数获得绝对含量。在获得同位素相对含量(组成比例)时,首先选择一种能够在大的能量区间内发射多条具有特征能量的、且能量分布均匀的γ射线的核素作为基准核素,求得探测系统对不同能量的伽马射线的相对探测效率。基于该探测效率,获得其他核素相对于基准核素的核子数之比。
[0003] 在钚同位素中,可以选择239Pu作为基准核素,根据γ射线能谱上相应的γ射线峰的数据,利用相对转换系数法可得到其余钚同位素与239Pu含量的比值。但是,242Pu主要的三条特征γ射线44.915keV、103.5keV、159.018keV,均不宜用作特征伽马射线,一方面是因为其强度弱,另一方面与钚的其它同位素的相近能量的伽马射线重叠而难以分辨,44.915keV的伽马射线与240Pu的45.244keV、241Pu的44.2keV的伽马射线难以分辨,103.5keV的伽马射线与239Pu的103.06keV、240Pu的104.234keV、241Pu的103.68keV的伽马射线难以分辨,159.018keV的伽马射线与240Pu的160.308keV、241Pu的159.955keV的伽马射线难以分辨。因此,无法通过直接探测242Pu的特征伽马射线来识别和计算242Pu的含量。
[0004] 乏燃料组件中的钚源自于238U在压堆内辐照中俘获中子转化过来的,因此,可以基于238U到钚主要同位素的燃耗链,研究得到乏燃料组件中钚的同位素含量之间的关系,由Pu元素其他核素的相对含量来获得242Pu的相对含量。
[0005] 目前有通过Pu元素其他核素来推算242Pu核素含量的计算方法,如美国加州大学劳伦斯辐射实验室(UCRL)研发的MGA程序所使用的公式:[242]=a+b[239]+c[240]+d[241]+e[239]2+f[240]2+g[241]2,其中a~g代表系数,中括号内的代表Pu元素各核素的含量;以及国内所使用的公式: 但是现有的计算方法并没有考虑238Pu核素含量的影响,而且计算结果与实际值之间误差较大,相对误差一般在50%左右,只能保证所得到的242Pu核素含量与实际值在一个数量级以内,实用性较低。

发明内容

[0006] 本发明的目的在于针对现有技术缺陷,提供一种与燃料组件辐照历史无关的242Pu核素相对含量的计算方法,这种方法处理简单,准确性高,能够有效得到242Pu核素的相对含量。
[0007] 本发明的技术方案如下:一种242Pu核素相对含量的计算方法,包括:
[0008] 利用γ谱仪测量238Pu、239Pu、240Pu、241Pu的特征γ射线,利用239Pu的特征γ射线,获得不同能量γ射线的相对探测效率曲线;
[0009] 根据238Pu、239Pu、240Pu、241Pu的特征γ射线的计数,结合γ射线探测效率随能量的变化,计算238Pu、240Pu、241Pu相对239Pu的含量,即238Pu/239Pu、240Pu/239Pu、241Pu/239Pu;
[0010] 根据如下公式计算出242Pu相对239Pu的含量:
[0011] y=a1+a2*x1+a3*x2+a4*x3+a5*x12+a6*x22+a7*x32+a8*x1*x2+a9*x1*x3+a10*x2*3 3 3 2 2 2 2 2
x3+a11*x1 +a12*x2+a13*x3+a14*x1 *x2+a15*x1 *x3+a16*x2*x1+a17*x2*x3+a18*x3 *x1+a19*x32*x2+a20*x1*x2*x3
[0012] 公式中,
[0013] x1为238Pu/239Pu相对含量;
[0014] x2为240Pu/239Pu相对含量;
[0015] x3为241Pu/239Pu相对含量;
[0016] y为242Pu/239Pu相对含量;
[0017] a1~a20为待定系数。
[0018] 进一步,如上所述的242Pu核素相对含量的计算方法,其中,238Pu、240Pu、241Pu相对239Pu的含量根据如下公式计算:
[0019]
[0020] 其中, 为238Pu或240Pu或241Pu相对于239Pu的含量;Ax,k、A239,j分别为238Pu或240Pu或241Pu的特征γ射线k、239Pu的特征γ射线j的全能峰计数;εj、εk分别为特征γ射线j、k的相对探测效率;px,k、p239,j是238Pu或240Pu或241Pu的特征γ射线k、239Pu的特征γ射线j的强度;λx、λ239是238Pu或240Pu或241Pu、239Pu的衰变常数。
[0021] 进一步,如上所述的242Pu核素相对含量的计算方法,其中,242Pu相对239Pu含量的计算公式是基于典型核电站,通过汇总不同燃料组件、不同燃耗深度下、不同冷却时间下钚同位素的含量原始数据,统计分析后得到的。
[0022] 进一步,242Pu相对239Pu含量的计算公式中a1~a20待定系数是基于典型核电站不同燃料组件、不同燃耗深度下、不同冷却时间下钚同位素的含量原始数据,利用LM(Levenberg-Marquard)最小二乘法非线性拟合最优化算法进行拟合计算分析后得到的。
[0023] 更进一步,242Pu相对239Pu含量的计算公式中a1~a20待定系数是根据Pu已知核素相对含量(238Pu/239Pu)的范围进行分组拟合分析得到,进一步提高计算精度
[0024] 本发明的有益效果如下:本发明根据多种堆型辐照历史、多种初始富集度、多种燃耗深度、多种冷却时间的乏燃料组件中钚各个核素含量的原始数据,研究得到242Pu/239Pu相238 239 240 239 241 239 238 239
对含量随 Pu/ Pu、 Pu/ Pu、 Pu/ Pu变化的通用计算公式,通过探测 Pu、 Pu、
240Pu、241Pu的特征γ射线,结合通用计算公式,得到242Pu的相对含量。该方法考虑了238Pu/
239Pu、240Pu/239Pu、241Pu/239Pu各个参量之间相互耦合的关系,应用LM(Levenberg-Marquard)最小二乘法非线性拟合最优化算法,收敛快,对冗余的参量不敏感,保证了计算方法的准确性。本发明有效解决了242Pu的特征γ射线与钚的其他同位素的γ射线重叠、难以分辨、无法直接探测的难题。在乏燃料组件的钚浓度测量中,本计算方法具有很好的通用性、适用性、准确性。计算所得242Pu含量的相对误差基本上在1%左右,大部分在1%以下。在钚浓度监测方面,242Pu并没有有效的测量手段,为准确地获得钚各核素的含量,242Pu含量的准确推算是必须的。通过本发明,可以有效、准确地测量出乏燃料中钚各种核素的含量,在乏燃料后处理领域可以有广泛的应用。
附图说明
[0025] 图1为本发明一种242Pu核素相对含量的计算方法流程图
[0026] 图2为本发明242Pu核素相对含量拟合公式与国内所使用公式:所得到的242Pu拟合值与原始值间的相对误差散点图。
[0027] 图3是本发明一种242Pu核素相对含量的计算方法按照238Pu/239Pu核素相对含量分组拟合和不分组拟合,242Pu拟合值与原始值间的相对误差散点图。

具体实施方式

[0028] 下面结合所附的图表、公式对本发明进行详细的描述。
[0029] 如图1所示本发明的计算流程,使用γ谱仪对含钚系统进行测量,根据钚各同位素的特征γ射线全能峰的计数,识别并确定238Pu、239Pu、240Pu、241Pu的相对含量。具体实施如239
下: Pu的特征γ射线129.296keV、203.55keV、332.845keV、375.054keV、413.713keV、
422.598keV和451.481keV的全能峰计数分别为A1、A2、A3、A4、A5、A6,假如以413.713keV的γ射线的全能峰计数为基准(选取其他全能峰计数也可以,如332.845keV),利用公式可得到各特征γ射线相对于基准特征γ射线的相对探测效率,相对探测效率仅与γ射线能量有关,与钚核素种类无关,由此可得到相对探测效率随γ射线能量在129~451keV范围内的变化:
[0030]
[0031] 其中:εi为相对探测效率,pi为核素衰变时放出相应能量的特征γ射线的强度。
[0032] 考虑到所得到的相对探测效率曲线的能量范围129~451keV,同时保证所选取的各核素特征γ峰不重叠、有一定强度、易于分辨,因此238Pu的特征γ射线可选用152.700keV,240Pu的特征γ射线可选用212.460keV,241Pu的特征γ射线可选用148.567keV,利用公式可得到238Pu、240Pu、241Pu相对于239Pu的含量:
[0033]
[0034] 其中, 为238Pu、240Pu、241Pu相对于239Pu的含量;Ax,k、A239,j分别为238Pu或240Pu或241Pu的特征γ射线k、239Pu的特征γ射线j的全能峰计数;εj、εk分别为特征γ射线j、k的相对探测效率;px,k、p239,j是238Pu或240Pu或241Pu的特征γ射线k、239Pu的特征γ射线j的强度;238 240 241 239
λx、λ239是 Pu或 Pu或 Pu、 Pu的衰变常数。
[0035] 基于国内典型核电站,统计汇总了不同燃料组件、不同燃耗深度下、不同冷却时间下钚同位素的含量原始数据,研究得到242Pu/239Pu的计算公式如下:
[0036] y=a1+a2*x1+a3*x2+a4*x3+a5*x12+a6*x22+a7*x32+a8*x1*x2+a9*x1*x3+a10*x2*x3+a11*x13+a12*x23+a13*x33+a14*x12*x2+a15*x12*x3+a16*x22*x1+a17*x22*x3+a18*x32*x1+a19*x32*x2+a20*x1*x2*x3
[0037] 公式中:
[0038] x1,238Pu/239Pu相对含量;
[0039] x2,240Pu/239Pu相对含量;
[0040] x3,241Pu/239Pu相对含量;
[0041] y,242Pu/239Pu相对含量;
[0042] a1~a20为待定系数。a1~a20是基于典型核电站不同燃料组件、不同燃耗深度下、不同冷却时间下钚同位素的含量原始数据,利用LM(Levenberg-Marquard)最小二乘法非线性拟合最优化算法进行拟合计算分析后得到的。
[0043] 根据238Pu/239Pu相对含量的不同,为提高计算精度,根据实际拟合的结果进行分析得到以0.01为界限(对238Pu/239Pu比值在0.01以下的数据进行分组拟合,可以明显改善计算238 239
结果的精度,减小误差),表1分别给出计算公式中的待定系数。该方法根据 Pu/ Pu相对含量的范围进行分组处理,可以进一步提高计算精度,如图3所示。
[0044] 表2是238Pu/239Pu相对含量小于0.01时,利用计算公式对测试组的242Pu/239Pu进行验证计算的误差分析,误差基本在1%以下,个别误差在1%~2%。
[0045] 表3是238Pu/239Pu相对含量大于0.01时,利用计算公式对测试组的242Pu/239Pu进行验证计算的误差分析,误差基本在1%以下,个别测试组误差在1%~2%。
[0046] 表1
[0047]
[0048] 表2
[0049]
[0050]
[0051] 表3
[0052]
[0053] 显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若对本发明的这些修改和变型属于本发明权利要求及其同等技术的范围之内,则本发明也意图包含这些改动和变型在内。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈