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

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

Info

Publication number
CN110990783B
CN110990783B CN201911112539.5A CN201911112539A CN110990783B CN 110990783 B CN110990783 B CN 110990783B CN 201911112539 A CN201911112539 A CN 201911112539A CN 110990783 B CN110990783 B CN 110990783B
Authority
CN
China
Prior art keywords
relative
content
gamma rays
characteristic gamma
characteristic
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
CN201911112539.5A
Other languages
English (en)
Other versions
CN110990783A (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.)
China Nuclear Power Engineering Co Ltd
Original Assignee
China Nuclear Power Engineering Co Ltd
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 China Nuclear Power Engineering Co Ltd filed Critical China Nuclear Power Engineering Co Ltd
Priority to CN201911112539.5A priority Critical patent/CN110990783B/zh
Publication of CN110990783A publication Critical patent/CN110990783A/zh
Application granted granted Critical
Publication of CN110990783B publication Critical patent/CN110990783B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • Y02E30/30Nuclear fission reactors

Abstract

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

Description

一种Pu-242核素相对含量的计算方法
技术领域
本发明属于核辐射与探测技术领域,涉及一种与燃料组件辐照历史无关的乏燃料后处理中242Pu核素相对含量的计算方法。
背景技术
乏燃料组件中钚同位素主要包括238Pu、239Pu、240Pu、241Pu和242Pu。通常,钚核素的识别主要通过测量中低能区特征伽马射线实现的。特征伽马射线法要求核素有能明显区别于其它核素的、特定能量的、可分辨的伽马射线。核素238Pu、239Pu、240Pu、241Pu都有能够用于核素识别的特征伽马射线,为了缓解对探测器定量标定的需求,降低对能量分辨率以及电子学系统性能的过高要求,通常采用相对探测效率法获得同位素相对含量(组成比例),再结合中子计数获得绝对含量。在获得同位素相对含量(组成比例)时,首先选择一种能够在大的能量区间内发射多条具有特征能量的、且能量分布均匀的γ射线的核素作为基准核素,求得探测系统对不同能量的伽马射线的相对探测效率。基于该探测效率,获得其他核素相对于基准核素的核子数之比。
在钚同位素中,可以选择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的含量。
乏燃料组件中的钚源自于238U在压水堆内辐照中俘获中子转化过来的,因此,可以基于238U到钚主要同位素的燃耗链,研究得到乏燃料组件中钚的同位素含量之间的关系,由Pu元素其他核素的相对含量来获得242Pu的相对含量。
目前有通过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核素含量与实际值在一个数量级以内,实用性较低。
发明内容
本发明的目的在于针对现有技术的缺陷,提供一种与燃料组件辐照历史无关的242Pu核素相对含量的计算方法,这种方法处理简单,准确性高,能够有效得到242Pu核素的相对含量。
本发明的技术方案如下:一种242Pu核素相对含量的计算方法,包括:
利用γ谱仪测量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+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
公式中,
x1为238Pu/239Pu相对含量;
x2为240Pu/239Pu相对含量;
x3为241Pu/239Pu相对含量;
y为242Pu/239Pu相对含量;
a1~a20为待定系数。
进一步,如上所述的242Pu核素相对含量的计算方法,其中,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,j238Pu或240Pu或241Pu的特征γ射线k、239Pu的特征γ射线j的强度;λx、λ239238Pu或240Pu或241Pu、239Pu的衰变常数。
进一步,如上所述的242Pu核素相对含量的计算方法,其中,242Pu相对239Pu含量的计算公式是基于典型核电站,通过汇总不同燃料组件、不同燃耗深度下、不同冷却时间下钚同位素的含量原始数据,统计分析后得到的。
进一步,242Pu相对239Pu含量的计算公式中a1~a20待定系数是基于典型核电站不同燃料组件、不同燃耗深度下、不同冷却时间下钚同位素的含量原始数据,利用LM(Levenberg-Marquard)最小二乘法非线性拟合最优化算法进行拟合计算分析后得到的。
更进一步,242Pu相对239Pu含量的计算公式中a1~a20待定系数是根据Pu已知核素相对含量(238Pu/239Pu)的范围进行分组拟合分析得到,进一步提高计算精度。
本发明的有益效果如下:本发明根据多种堆型辐照历史、多种初始富集度、多种燃耗深度、多种冷却时间的乏燃料组件中钚各个核素含量的原始数据,研究得到242Pu/239Pu相对含量随238Pu/239Pu、240Pu/239Pu、241Pu/239Pu变化的通用计算公式,通过探测238Pu、239Pu、240Pu、241Pu的特征γ射线,结合通用计算公式,得到242Pu的相对含量。该方法考虑了238Pu/239Pu、240Pu/239Pu、241Pu/239Pu各个参量之间相互耦合的关系,应用LM(Levenberg-Marquard)最小二乘法非线性拟合最优化算法,收敛快,对冗余的参量不敏感,保证了计算方法的准确性。本发明有效解决了242Pu的特征γ射线与钚的其他同位素的γ射线重叠、难以分辨、无法直接探测的难题。在乏燃料组件的钚浓度测量中,本计算方法具有很好的通用性、适用性、准确性。计算所得242Pu含量的相对误差基本上在1%左右,大部分在1%以下。在钚浓度监测方面,242Pu并没有有效的测量手段,为准确地获得钚各核素的含量,242Pu含量的准确推算是必须的。通过本发明,可以有效、准确地测量出乏燃料中钚各种核素的含量,在乏燃料后处理领域可以有广泛的应用。
附图说明
图1为本发明一种242Pu核素相对含量的计算方法流程图。
图2为本发明242Pu核素相对含量拟合公式与国内所使用公式: 所得到的242Pu拟合值与原始值间的相对误差散点图。
图3是本发明一种242Pu核素相对含量的计算方法按照238Pu/239Pu核素相对含量分组拟合和不分组拟合,242Pu拟合值与原始值间的相对误差散点图。
具体实施方式
下面结合所附的图表、公式对本发明进行详细的描述。
如图1所示本发明的计算流程,使用γ谱仪对含钚系统进行测量,根据钚各同位素的特征γ射线全能峰的计数,识别并确定238Pu、239Pu、240Pu、241Pu的相对含量。具体实施如下:239Pu的特征γ射线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范围内的变化:
其中:εi为相对探测效率,pi为核素衰变时放出相应能量的特征γ射线的强度。
考虑到所得到的相对探测效率曲线的能量范围129~451keV,同时保证所选取的各核素特征γ峰不重叠、有一定强度、易于分辨,因此238Pu的特征γ射线可选用152.700keV,240Pu的特征γ射线可选用212.460keV,241Pu的特征γ射线可选用148.567keV,利用公式可得到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,j238Pu或240Pu或241Pu的特征γ射线k、239Pu的特征γ射线j的强度;λx、λ239238Pu或240Pu或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+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
公式中:
x1,238Pu/239Pu相对含量;
x2,240Pu/239Pu相对含量;
x3,241Pu/239Pu相对含量;
y,242Pu/239Pu相对含量;
a1~a20为待定系数。a1~a20是基于典型核电站不同燃料组件、不同燃耗深度下、不同冷却时间下钚同位素的含量原始数据,利用LM(Levenberg-Marquard)最小二乘法非线性拟合最优化算法进行拟合计算分析后得到的。
根据238Pu/239Pu相对含量的不同,为提高计算精度,根据实际拟合的结果进行分析得到以0.01为界限(对238Pu/239Pu比值在0.01以下的数据进行分组拟合,可以明显改善计算结果的精度,减小误差),表1分别给出计算公式中的待定系数。该方法根据238Pu/239Pu相对含量的范围进行分组处理,可以进一步提高计算精度,如图3所示。
表2是238Pu/239Pu相对含量小于0.01时,利用计算公式对测试组的242Pu/239Pu进行验证计算的误差分析,误差基本在1%以下,个别误差在1%~2%。
表3是238Pu/239Pu相对含量大于0.01时,利用计算公式对测试组的242Pu/239Pu进行验证计算的误差分析,误差基本在1%以下,个别测试组误差在1%~2%。
表1
表2
表3
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若对本发明的这些修改和变型属于本发明权利要求及其同等技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (5)

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+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
公式中,
x1为238Pu/239Pu相对含量;
x2为240Pu/239Pu相对含量;
x3为241Pu/239Pu相对含量;
y为242Pu/239Pu相对含量;
a1~a20为待定系数;
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,j238Pu或240Pu或241Pu的特征γ射线k、239Pu的特征γ射线j的强度;λx、λ239238Pu或240Pu或241Pu、239Pu的衰变常数。
2.如权利要求1所述的Pu-242核素相对含量的计算方法,其特征在于,242Pu相对239Pu含量的计算公式是基于典型核电站,通过汇总不同燃料组件、不同燃耗深度下、不同冷却时间下钚同位素的含量原始数据,统计分析后得到的。
3.如权利要求2所述的Pu-242核素相对含量的计算方法,其特征在于,242Pu相对239Pu含量的计算公式中a1~a20待定系数是基于典型核电站不同燃料组件、不同燃耗深度下、不同冷却时间下钚同位素的含量原始数据,利用LM最小二乘法非线性拟合最优化算法进行拟合计算分析后得到的。
4.如权利要求3所述的Pu-242核素相对含量的计算方法,其特征在于,242Pu相对239Pu含量的计算公式中a1~a20待定系数,根据Pu已知核素相对含量的范围进行分组拟合计算分析得到。
5.如权利要求4所述的Pu-242核素相对含量的计算方法,其特征在于,a1~a20待定系数是根据238Pu/239Pu核素相对含量进行分组拟合计算分析得到。
CN201911112539.5A 2019-11-14 2019-11-14 一种Pu-242核素相对含量的计算方法 Active CN110990783B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911112539.5A CN110990783B (zh) 2019-11-14 2019-11-14 一种Pu-242核素相对含量的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911112539.5A CN110990783B (zh) 2019-11-14 2019-11-14 一种Pu-242核素相对含量的计算方法

Publications (2)

Publication Number Publication Date
CN110990783A CN110990783A (zh) 2020-04-10
CN110990783B true CN110990783B (zh) 2024-03-19

Family

ID=70084352

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911112539.5A Active CN110990783B (zh) 2019-11-14 2019-11-14 一种Pu-242核素相对含量的计算方法

Country Status (1)

Country Link
CN (1) CN110990783B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB8512068D0 (en) * 1984-05-23 1985-06-19 British Nuclear Fuels Plc Measuring plutonium concentration in process liquid streams
US5629525A (en) * 1993-08-12 1997-05-13 British Nuclear Fuels Plc Plutonium assaying
JP2001235547A (ja) * 1999-12-13 2001-08-31 Japan Atom Energy Res Inst 多重ガンマ線検出による高感度核種分析方法
JP2005043173A (ja) * 2003-07-28 2005-02-17 Toshiba Corp 燃料物質内プルトニウム238の非破壊測定方法
JP2012225663A (ja) * 2011-04-15 2012-11-15 Toshiba Corp Pu定量分析方法
CN107656305A (zh) * 2017-08-15 2018-02-02 中国原子能科学研究院 一种测定铀富集度的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8194813B2 (en) * 2004-10-19 2012-06-05 Lawrence Livermore National Security, Llc Absolute nuclear material assay using count distribution (LAMBDA) space

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB8512068D0 (en) * 1984-05-23 1985-06-19 British Nuclear Fuels Plc Measuring plutonium concentration in process liquid streams
US5629525A (en) * 1993-08-12 1997-05-13 British Nuclear Fuels Plc Plutonium assaying
JP2001235547A (ja) * 1999-12-13 2001-08-31 Japan Atom Energy Res Inst 多重ガンマ線検出による高感度核種分析方法
JP2005043173A (ja) * 2003-07-28 2005-02-17 Toshiba Corp 燃料物質内プルトニウム238の非破壊測定方法
JP2012225663A (ja) * 2011-04-15 2012-11-15 Toshiba Corp Pu定量分析方法
CN107656305A (zh) * 2017-08-15 2018-02-02 中国原子能科学研究院 一种测定铀富集度的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
一种钚同位素分析方法及实现;林源根;郭春营;段占元;;核电子学与探测技术(04);第757-760页 *
基于多核素的相对效率法;胡广春;向永春;伍怀龙;刘晓亚;;核技术(03);第030503-1~030503-5页 *
多道γ能谱的精细分析;王思广;冒亚军;唐培家;李泽;;核技术(07);第495-498页 *

Also Published As

Publication number Publication date
CN110990783A (zh) 2020-04-10

Similar Documents

Publication Publication Date Title
Kim et al. Analysis of high burnup pressurized water reactor fuel using uranium, plutonium, neodymium, and cesium isotope correlations with burnup
CN110990783B (zh) 一种Pu-242核素相对含量的计算方法
Roudil et al. CETAMA contribution to safeguards and nuclear forensic analysis based on nuclear reference materials
Gauld et al. Investigation of passive gamma spectroscopy to verify spent nuclear fuel content
Gauld et al. In-field performance testing of the fork detector for quantitative spent fuel verification
Lee et al. Determination of the plutonium mass and curium ratio of spent fuel assemblies for input nuclear material accountancy of pyroprocessing, and analysis of their errors
Phillips et al. Neutron measurement techniques for the nondestructive analysis of irradiated fuel assemblies
CN209879001U (zh) 一种中子多重性测量装置
EP0042099A2 (en) Self-regulating neutron coincidence counter
Sung et al. Experimental characterization of the accuracy of multidetector boron meter for operational safety of reactors
Koskelo et al. Sustainability of gamma-ray isotopic evaluation codes
Gunnink Plutonium isotopic analysis of nondescript samples by gamma-ray spectrometry
Abousahl et al. The MGA code for the determination of the isotopic composition of plutonium and MOX by gamma spectrometry—A performance study
Cheatham et al. Determining Spent Nuclear Fuel's Plutonium Content, Initial Enrichment, Burnup, and Cooling Time
Gunnink et al. TASTEX gamma spectrometer system for measuring isotopic and total plutonium concentrations in solutions
Walton Nondestructive Assay
JPH04249797A (ja) 照射燃料集合体の燃焼度測定方法
McNeece et al. The Graphite Isotope Ratio Method (GIRM): A Plutonium Production Verification Tool
Ilas et al. New Radiochemical Assay Experiments to Improve the Validation Basis for High Burnup Fuel
Grogan et al. Verifying Safeguards Declarations with INDEPTH: A Sensitivity Study
Gunnink et al. PAT, a Program for Determining the Category of Plutonium Samples
Chaudhury et al. Non-destructive assay of plutonium in absence of gamma-ray spectrometry
Lin et al. Attribution of nuclear material by non-destructive radiation measurement methods
Croft et al. Considerations on the use of an Anti-Compton Suppression Technique for the Measurement of Low Energy Gamma Rays in Dense Waste Materials-18144
Villani et al. Optimal Differential Die-Away Early Gate Study for Difficult Assays

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