CN109949953A - 核燃料棒富集度无源检测系统及富集度重构方法 - Google Patents
核燃料棒富集度无源检测系统及富集度重构方法 Download PDFInfo
- Publication number
- CN109949953A CN109949953A CN201910264940.4A CN201910264940A CN109949953A CN 109949953 A CN109949953 A CN 109949953A CN 201910264940 A CN201910264940 A CN 201910264940A CN 109949953 A CN109949953 A CN 109949953A
- Authority
- CN
- China
- Prior art keywords
- detector
- nuclear fuel
- fuel rod
- energy peak
- pellet
- 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
Links
- 239000003758 nuclear fuel Substances 0.000 title claims abstract description 107
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000001514 detection method Methods 0.000 title claims abstract description 38
- 239000008188 pellet Substances 0.000 claims abstract description 87
- 230000004044 response Effects 0.000 claims abstract description 27
- 238000005259 measurement Methods 0.000 claims description 19
- JFALSRSLKYAFGM-OIOBTWANSA-N uranium-235 Chemical compound [235U] JFALSRSLKYAFGM-OIOBTWANSA-N 0.000 claims description 10
- 230000005251 gamma ray Effects 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 6
- 241000287196 Asthenes Species 0.000 claims 1
- 239000000446 fuel Substances 0.000 description 10
- 230000008569 process Effects 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000004992 fission Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 2
- 230000002285 radioactive effect Effects 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000005255 beta decay Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- OOAWCECZEHPMBX-UHFFFAOYSA-N oxygen(2-);uranium(4+) Chemical compound [O-2].[O-2].[U+4] OOAWCECZEHPMBX-UHFFFAOYSA-N 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000008646 thermal stress Effects 0.000 description 1
- FCTBKIHDJGHPPO-UHFFFAOYSA-N uranium dioxide Inorganic materials O=[U]=O FCTBKIHDJGHPPO-UHFFFAOYSA-N 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Monitoring And Testing Of Nuclear Reactors (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明提供一种核燃料棒富集度无源检测系统及富集度重构方法。系统包括多个依次串行等间距排列的中间有孔洞的探测器。方法包括获取探测器对核燃料棒中所有芯块的响应分布；将各探测器所探测到的目标芯块全能峰计数率的平均值作为迭代初始值；通过EM迭代算法的正投影，依据响应分布、迭代的初始值以及所有探测器在核燃料棒运动时间t处对所有芯块的测量值得到估计值；将目标芯块经过各个探测器时各个探测器相应的测量值叠加得到真实值；通过EM迭代算法的反投影，依据估计值与真实值的比值对全能峰计数率进行修正；以及将进行修正后的全能峰计数率作为新的迭代的初始值重复依次进行EM迭代算法的正投影和反投影直至得到达到预期标准的全能峰计数率。
Description
技术领域
本发明涉及核安全领域和核探测数据处理领域，具体涉及一种核燃料棒富集度无源检测系统及富集度重构方法。
背景技术
压水堆是我国核电站常用的反应堆。核燃料棒作为反应堆的关键组成，由多个二氧化铀(UO2)芯块组成。在燃料芯块中，将核燃料中235U(铀235)的质量分数称为核燃料芯块的富集度，压水堆核燃料芯块的富集度一般为2％～5％。在反应堆运行时，燃料棒中芯块富集度的不均匀(非设计的整体偏高/偏低)或存在异常单个芯块(不同富集度的燃料芯块混入燃料棒)，均会导致燃料棒发热不均匀，产生热应力，进而导致燃料棒破裂，停堆，引发安全事故和重大损失。因此，对每根燃料棒的芯块富集度都要进行100％检测。
核燃料富集度的检测方法常见的为有源检测法。有源检测法是将慢化后的252Cf中子照射燃料棒，热中子与238U发生反应238U(n，γ)→239Np→239Pu生成239Pu，与235U发生反应235U(n，f)→裂变产物+中子，产生的裂变产物不稳定，伴随一系列β衰变，释放出大量缓发伽马射线，计算某种裂变产物核素和239Pu的特征伽马射线全能峰净面积可以推算出235U与238U的比值。这种方法需要使用中子源对核燃料棒进行活化，因而检测成本高，同时也给核燃料棒的存储、运输、操作带来了安全隐患。
因此，需要一种新的核燃料棒富集度检测系统及方法。
在所述背景技术部分公开的上述信息仅用于加强对本发明的背景的理解，因此它可以包括不构成对本领域普通技术人员已知的现有技术的信息。
发明内容
本发明的目的在于提供一种核燃料棒富集度无源检测系统及富集度重构方法，进而至少在一定程度上克服由于相关技术的限制和缺陷而导致的一个或者多个问题。
本发明的其他特性和优点将通过下面的详细描述变得显然，或部分地通过本发明的实践而习得。
根据本发明的第一方面，公开一种核燃料棒富集度无源检测的系统，其特征在于，包括：
多个依次串行等间距排列的探测器，其中探测器为多边形伽马射线探测器，且探测器的中间具有允许核燃料棒通过的孔洞；以及
处理单元，用于接收和处理探测器的测量结果；
其中核燃料棒由多个大小相同且沿核燃料棒的轴向分布的芯块构成，测量时核燃料棒匀速经由孔洞依次经过多个探测器，在一个单位时间内核燃料棒移动与探测器的间距相同的距离且探测器进行一次测量。
根据本发明的一实施方式，其中多边形伽马射线探测器为八边形；以及多边形伽马射线探测器的数量为128个。
根据本发明的第二方面，公开一种基于EM迭代算法以及如任一前述的系统的核燃料棒富集度重构方法，其特征在于，包括：
获取探测器对核燃料棒中所有芯块的响应分布；
将各个探测器所探测到的目标芯块全能峰计数率的平均值作为迭代的初始值；
通过EM迭代算法的正投影，依据响应分布、迭代的初始值以及所有探测器在核燃料棒运动时间t处对所有芯块的测量值得到估计值；
将目标芯块经过各个探测器时各个探测器相应的测量值叠加得到真实值，记为countt；
通过EM迭代算法的反投影，依据估计值与真实值的比值对全能峰计数率进行修正；以及
将进行修正后的全能峰计数率作为新的迭代的初始值重复依次进行EM迭代算法的正投影和反投影直至得到达到预期标准的全能峰计数率。
根据本发明的一实施方式，所述方法还包括：依据实验标定的全能峰计数率-富集度线性关系以及达到预期标准的全能峰计数率，进行核燃料棒中各芯块的核燃料富集度的重构。
根据本发明的一实施方式，其中获取探测器对核燃料棒中所有芯块的响应分布包括：
测量一个单个的芯块移动经过所有探测器，所有探测器各个时刻的全能峰计数率，并将测量结果进行归一化处理以作为核燃料棒中第一个芯块的响应分布，记为则探测器对核燃料棒中所有芯块的响应分布为：
其中j为探测器标号，i为芯块标号，j、i及t均为整数，j∈[1，J]，J为探测器的总个数，i∈[1，I]，I为芯块的总个数。
根据本发明的一实施方式，其中将各个探测器所探测到的目标芯块全能峰计数率的平均值作为迭代的初始值包括：
将目标芯块经过各个探测器时各个探测器相应的测量值叠加后取平均值作为迭代的初始值。
根据本发明的一实施方式，其中估计值为：
其中n∈[0，N]，N为得到达到预期标准的全能峰计数率所需的迭代次数，delta[j]为核燃料棒器由第1个探测器移动至第j个探测器所需的时间，ei n为第i个芯块经过第n次迭代修正后的全能峰计数率。
根据本发明的一实施方式，其中进行修正后的全能峰计数率表示为；
其中
根据本发明的一实施方式，其中核燃料棒中的核燃料为铀235。
根据本发明的一实施方式，其中核燃料棒中的核燃料为铀235，以及所述全能峰计数率富集度线性关系为铀235的全能峰计数率富集度线性关系。
根据本发明的一些实施方式，通过考虑在相同时刻，核燃料棒中多个芯块均会放射特征伽马射线并被同一探测器探测的实际情况，从而使富集度的重构结果更加精确。
根据本发明的一些实施方式，仅需对一个芯块移动通过探测器时的计数进行测量，便可得到探测器对核燃料棒中所有芯块的响应分布，操作简便快捷。
应当理解的是，以上的一般描述和后文的细节描述仅是示例性的，并不能限制本发明。
附图说明
通过参照附图详细描述其示例实施例，本发明的上述和其它目标、特征及优点将变得更加显而易见。
图1示出根据本发明一示例实施方式的一核燃料棒富集度无源检测系统的示意图。
图2示出核燃料棒富集度无源检测系统的实体实拍图。
图3示出根据本发明一示例实施方式的一核燃料棒富集度重构方法的流程图。
图4示出根据本发明另一示例实施方式的一核燃料棒富集度重构方法的流程图。
具体实施方式
现在将参考附图更全面地描述示例实施方式。然而，示例实施方式能够以多种形式实施，且不应被理解为限于在此阐述的范例；相反，提供这些实施方式使得本发明的描述将更加全面和完整，并将示例实施方式的构思全面地传达给本领域的技术人员。附图仅为本发明的示意性图解，并非一定是按比例绘制。图中相同的附图标记表示相同或类似的部分，因而将省略对它们的重复描述。
此外，所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施方式中。在下面的描述中，提供许多具体细节从而给出对本发明的实施方式的充分理解。然而，本领域技术人员将意识到，可以实践本发明的技术方案而省略所述特定细节中的一个或更多，或者可以采用其它的方法、组元、步骤等。在其它情况下，不详细示出或描述公知结构、方法、实现或者操作以避免喧宾夺主而使得本发明的各方面变得模糊。
附图中所示的一些方框图是功能实体，不一定必须与物理或逻辑上独立的实体相对应。可以采用软件形式来实现这些功能实体，或在一个或多个硬件模块或集成电路中实现这些功能实体，或在不同网络和/或处理器装置和/或微控制器装置中实现这些功能实体。
本发明公开一种核燃料棒富集度无源检测系统及富集度重构方法。检测系统包括多个依次串行等间距排列的探测器，其中探测器为多边形伽马射线，且探测器的中间具有允许核燃料棒通过的孔洞。核燃料棒富集度重构方法包括获取探测器对核燃料棒中所有芯块的响应分布；初始化目标芯块所发射伽马特征射线的全能峰计数率作为迭代的初始值；通过EM迭代算法的正投影，依据响应分布、迭代的初始值以及所有探测器在核燃料棒运动时间t处对所有芯块的测量值得到估计值；将目标芯块经过各个探测器时各个探测器相应的测量值叠加得到真实值；通过EM迭代算法的反投影，依据估计值与真实值的比值对全能峰计数率进行修正；以及将进行修正后的全能峰计数率作为新的迭代的初始值重复依次进行EM迭代算法的正投影和反投影直至得到达到预期标准的全能峰计数率。本发明通过考虑在相同时刻，核燃料棒中多个芯块均会放射特征伽马射线并被同一探测器探测的实际情况，从而使富集度的重构结果更加精确。同时仅需对一个芯块移动通过探测器时的计数进行测量，便可得到探测器对核燃料棒中所有芯块的响应分布，操作简便快捷。
与有源检测法不同，无源检测法是直接测量235U的伽马特征峰计数，由于核燃料富集度的不同，235U的97keV与185.7keV天然特征峰的计数和会发生显著变化，借助这种关系，推算出核燃料棒的富集度。相比于有源检测法，无源检测无需额外放射源，成本低，操作安全，不需要处理放射源废源，是核燃料富集度检测的发展趋势。
中国科学院高能物理研究所协助北方核燃料公司研制的核燃料富集度快速无源检测系统，其由32*4个8边形探测器及其后端电子学和配套软件构成，每个探测器中间开有小孔，测量过程中，核燃料棒从中间小孔通过，探测器获取计数。每个核燃料棒由多个核燃料芯块组成，在富集度无源检测过程中，核燃料棒通过探测器单元并被探测。一种核燃料富集度重构方式可以如下所述：假设每个探测器单位时间探测到的计数仅来自于燃料棒中正对于该探测器的芯块。考虑到核燃料棒匀速运动，这种重构方法通过获取各个探测器每一时刻伽马特征射线能谱的全能峰计数，将每个探测器开始产生计数的时间点作为该探测器的计数时间零点，把所有探测器相同的计数时间点的全能峰计数率进行叠加后取平均，可得到所有探测器对各个芯块产生的特征伽马射线探测的平均探测结果，以此作为每个芯块所放射出的伽马射线的强度，经过对比获取各个芯块的富集度。以上重构方式可称为平均值重构方式。
实际上探测器获取计数时，除了记录正对于其的芯块产生的伽马射线，也会记录其他位置芯块产生的伽马射线。而平均值重构方式忽略了其他位置芯块产生的伽马射线对探测器的影响，所重构得到的核燃料富集度不够准确。
为了对核燃料富集度进行更精确的计算和重构，本发明提出如图1-3所示的核燃料棒富集度无源检测系统及富集度重构方法。
下面结合图1-4对本发明的核燃料棒富集度无源检测系统及富集度重构方法进行详细说明，其中，图1示出根据本发明一示例实施方式的一核燃料棒富集度无源检测系统的示意图；图2示出核燃料棒富集度无源检测系统的实体实拍图；图3示出根据本发明一示例实施方式的一核燃料棒富集度重构方法的流程图；图4示出根据本发明另一示例实施方式的一核燃料棒富集度重构方法的流程图。
首先结合图1-2对本发明一示例实施方式的一核燃料棒富集度无源检测系统进行具体说明。
如图1所示，核燃料棒富集度无源检测系统包括：多个依次串行等间距排列的探测器1，其中探测器为多边形伽马射线，且探测器的中间具有允许核燃料棒2通过的孔洞11；以及处理单元(未图示)，用于接收和处理探测器的测量结果；其中核燃料棒由多个大小相同且沿核燃料棒的轴向分布的芯块构成，测量时核燃料棒匀速经由孔洞依次经过多个探测器，在一个单位时间内核燃料棒移动与探测器的间距相同的距离且探测器进行一次测量。图2示出核燃料棒富集度无源检测系统的实体实拍图。
根据本发明的一实施方式，其中多边形伽马射线探测器为八边形；以及多边形伽马射线探测器的数量为128个。
下面结合图3对本发明的基于EM迭代算法以及如任一前述的核燃料棒富集度无源检测系统的核燃料棒富集度重构方法进行说明。其中EM迭代算法，指的是最大期望算法(Expectation Maximization Algorithm，又译期望最大化算法)，是一种迭代算法，在统计学中被用于寻找，依赖于不可观察的隐性变量的概率模型中，参数的最大似然估计。
图3示出根据本发明一示例实施方式的一核燃料棒富集度重构方法的流程图。
如图3所示，在S31，获取探测器对核燃料棒中所有芯块的响应分布。响应分布是指核燃料棒每移动的一个单位距离，各个探测器对每个芯块的响应，即各个探测器探测到芯块所发射的伽马特征射线的概率。它是由探测器标号j、芯块标号i以及核燃料棒运动时间t三个信息索引的三维矩阵，记为其中j、i及t均为整数，j∈[1，J]，J为探测器的总个数，i∈[1，I]，I为芯块的总个数。为获取响应分布，测量一个芯块移动经过所有探测器，探测器各个时刻的全能峰计数，并将测量结果进行归一化处理，以此作为核燃料棒中第一个芯块的响应分布，记为在某一特定时间t0，第一个芯块正对于探测器j0，由于第二个芯块与第一个芯块相邻，此时探测器j0对第二个芯块的响应与第一个芯块在t0-1时刻的探测器j0的响应相同，即：
依次类推，对于芯块i，它的响应分布可以表示为：
借助这种方法，只需对一个芯块移动通过探测器时的计数进行测量，便可得到探测器对核燃料棒中所有芯块的响应分布。
在S32，初始化目标芯块i所发射的伽马特征射线的全能峰计数率作为迭代的初始值。本发明基于EM迭代算法，该算法对初始值的选定并不敏感。但为尽可能减少的迭代次数，将前述的平均值重构方式中所得的各芯块235U(铀235)的全能峰计数率作为迭代的初始值，记为
也就是说，根据本发明的一实施方式，其中初始化目标芯块i所发射的伽马特征射线的全能峰计数率作为迭代的初始值包括：将目标芯块i经过各个探测器时各个探测器相应的测量值叠加后取平均值作为迭代的初始值。
在S33，通过EM迭代算法的正投影，依据响应分布、迭代的初始值以及所有探测器在核燃料棒运动时间t处对所有芯块的测量值得到估计值。
正投影是实际探测值的估算过程。各个探测器在每个单位时间都会记录一个全能峰计数值，此计数是对应探测器探测到每个芯块中235U所放射出的伽马特征射线计数的叠加，因而该计数值可以估计为：
其中n为迭代次数，n∈[0，N]，N为得到达到预期标准的全能峰计数率所需的迭代次数，ei n为第i个芯块经过第n次迭代修正后的全能峰计数率。
在核探测中计数值n的统计误差，可以表示为：
因而计数值n越大，统计误差越小，为了减小估计结果的统计误差，将每个探测器开始产生计数的时间作为时间零点，把各探测器的时间零点对齐，将每个对应时间的探测器计数值相加，作为最终估计值：
其中delta[j]指核燃料棒器由第1个探测器移动至第j个探测器所需的时间。
在S34，将目标芯块经过各个探测器时各个探测器相应的测量值叠加得到真实值，记为countt。
在S35，通过EM迭代算法的反投影，依据估计值与真实值的比值对全能峰计数率进行修正。
反投影是利用投影值与真实值的比值(6)对核燃料芯块中235U的全能峰计数率进行修正，可表示为：
在S36，将进行修正后的全能峰计数率作为新的迭代的初始值重复依次进行EM迭代算法的正投影和反投影直至得到达到预期标准的全能峰计数率。也就是说，将每次修正/更新后的计数率重新代入到正投影过程，不断迭代，直至获得理想的迭代结果。
根据本发明的一实施方式，所述方法还包括：依据实验标定的全能峰计数率-富集度线性关系以及达到预期标准的全能峰计数率，进行核燃料棒中各芯块的核燃料富集度的重构，如附图4中的S47所示。其中，S41-S46与S31-S36相同，在此不再赘述。。
根据本发明的一实施方式，其中核燃料棒中的核燃料为铀235，以及所述全能峰计数率-富集度线性关系为铀235的全能峰计数率-富集度线性关系。
表1示出平均值重构方式与基于EM算法异常芯块富集度检测的结果对比。
表1
由表1可以看到基于EM算法的异常芯块富集度检测/重构方法检测到的富集度偏差与富集度实际偏差极为接近。
通过以上的详细描述，本领域的技术人员易于理解，根据本发明实施例的系统和方法具有以下优点中的一个或多个。
根据本发明的一些实施方式，通过考虑在相同时刻，核燃料棒中多个芯块均会放射特征伽马射线并被同一探测器探测的实际情况，从而使富集度的重构结果更加精确。
根据本发明的一些实施方式，仅需对一个芯块移动通过探测器时的计数进行测量，便可得到探测器对核燃料棒中所有芯块的响应分布，操作简便快捷。
本领域技术人员在考虑说明书及实践这里公开的发明后，将容易想到本发明的其它实施方案。本申请旨在涵盖本发明的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本发明的真正范围和精神由下面的权利要求指出。
应当理解的是，本发明并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围进行各种修改和改变。本发明的范围仅由所附的权利要求来限制。
Claims (10)
1.一种核燃料棒富集度无源检测的系统，其特征在于，包括：
多个依次串行等间距排列的探测器，其中探测器为多边形伽马射线探测器，且探测器的中间具有允许核燃料棒通过的孔洞；以及
处理单元，用于接收和处理探测器的测量结果；
其中核燃料棒由多个大小相同且沿核燃料棒的轴向分布的芯块构成，测量时核燃料棒匀速经由孔洞依次经过多个探测器，在一个单位时间内核燃料棒移动与探测器的间距相同的距离且探测器进行一次测量。
2.如权利要求1所述的系统，其中多边形伽马射线探测器为八边形；以及多边形伽马射线探测器的数量为128个。
3.一种基于EM迭代算法以及如权利要求1-2任一所述的系统的核燃料棒富集度重构方法，其特征在于，包括：
获取探测器对核燃料棒中所有芯块的响应分布；
将各个探测器所探测到的目标芯块全能峰计数率的平均值作为迭代的初始值；
通过EM迭代算法的正投影，依据响应分布、迭代的初始值以及所有探测器在核燃料棒运动时间t处对所有芯块的测量值得到估计值；
将目标芯块经过各个探测器时各个探测器相应的测量值叠加得到真实值，记为countt；
通过EM迭代算法的反投影，依据估计值与真实值的比值对全能峰计数率进行修正；以及
将进行修正后的全能峰计数率作为新的迭代的初始值重复依次进行EM迭代算法的正投影和反投影直至得到达到预期标准的全能峰计数率。
4.如权利要求3所述的方法，还包括：依据实验标定的全能峰计数率-富集度线性关系以及达到预期标准的全能峰计数率，进行核燃料棒中各芯块的核燃料富集度的重构。
5.如权利要求3或4所述的方法，其中获取探测器对核燃料棒中所有芯块的响应分布包括：
测量一个单个的芯块移动经过所有探测器，所有探测器各个时刻的全能峰计数率，并将测量结果进行归一化处理以作为核燃料棒中第一个芯块的响应分布，记为则探测器对核燃料棒中所有芯块的响应分布为：
其中j为探测器标号，i为芯块标号，j、i及t均为整数，j∈[1，J]，J为探测器的总个数，i∈[1，I]，I为芯块的总个数。
6.如权利要求5所述的方法，其中将各个探测器所探测到的目标芯块全能峰计数率的平均值作为迭代的初始值包括：
将目标芯块经过各个探测器时各个探测器相应的测量值叠加后取平均值作为迭代的初始值。
7.如权利要求6所述的方法，其中估计值为：
其中n∈[0，N]，N为得到达到预期标准的全能峰计数率所需的迭代次数，delta[j]为核燃料棒器由第1个探测器移动至第j个探测器所需的时间，ei n为第i个芯块经过第n次迭代修正后的全能峰计数率。
8.如权利要求7所述的方法，其中进行修正后的全能峰计数率表示为；
其中
9.如权利要求3或8所述的方法，其中核燃料棒中的核燃料为铀235。
10.如权利要求3或8所述的方法，其中核燃料棒中的核燃料为铀235，以及所述全能峰计数率-富集度线性关系为铀235的全能峰计数率-富集度线性关系。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910264940.4A CN109949953B (zh) | 2019-04-02 | 2019-04-02 | 核燃料棒富集度无源检测系统及富集度重构方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910264940.4A CN109949953B (zh) | 2019-04-02 | 2019-04-02 | 核燃料棒富集度无源检测系统及富集度重构方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109949953A true CN109949953A (zh) | 2019-06-28 |
| CN109949953B CN109949953B (zh) | 2021-07-13 |
Family
ID=67012539
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910264940.4A Active CN109949953B (zh) | 2019-04-02 | 2019-04-02 | 核燃料棒富集度无源检测系统及富集度重构方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109949953B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111736201A (zh) * | 2020-05-25 | 2020-10-02 | 中国核电工程有限公司 | 一种核燃料棒有源检测系统及方法 |
| CN113050151A (zh) * | 2021-03-15 | 2021-06-29 | 中国原子能科学研究院 | 一种识别放射性核素的方法以及测量系统 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4822552A (en) * | 1987-02-25 | 1989-04-18 | Westinghouse Electric Corp. | Method and apparatus for passively gamma scanning a nuclear fuel rod |
| CN101283913A (zh) * | 2008-05-30 | 2008-10-15 | 首都师范大学 | Ct图像重建的gpu加速方法 |
| JP2009156724A (ja) * | 2007-12-27 | 2009-07-16 | Toshiba Corp | 燃料棒検査方法、燃料棒検査装置および燃料集合体の製造方法 |
| CN201378438Y (zh) * | 2009-03-04 | 2010-01-06 | 中国核电工程有限公司 | 一种252Cf中子活化核燃料棒232U富集度及均匀性的γ射线测量探测器 |
| CN102565845A (zh) * | 2010-12-31 | 2012-07-11 | 同方威视技术股份有限公司 | 利用多个探测器的伽马能谱核素识别方法 |
| CN103163548A (zh) * | 2013-03-07 | 2013-06-19 | 北京辛耕普华医疗科技有限公司 | 基于伽马相机的放射性物质探测方法及其装置和系统 |
| CN106596496A (zh) * | 2017-01-11 | 2017-04-26 | 中核建中核燃料元件有限公司 | 一种用于燃料棒丰度检测的γ探头及探头系统 |
| CN107133462A (zh) * | 2017-04-27 | 2017-09-05 | 上海联影医疗科技有限公司 | 数据处理方法、装置及设备 |
| WO2018067208A2 (en) * | 2016-06-09 | 2018-04-12 | Phoenix Nuclear Labs Llc | System and method for performing active scanning of a nuclear fuel rod |
-
2019
- 2019-04-02 CN CN201910264940.4A patent/CN109949953B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4822552A (en) * | 1987-02-25 | 1989-04-18 | Westinghouse Electric Corp. | Method and apparatus for passively gamma scanning a nuclear fuel rod |
| JP2009156724A (ja) * | 2007-12-27 | 2009-07-16 | Toshiba Corp | 燃料棒検査方法、燃料棒検査装置および燃料集合体の製造方法 |
| CN101283913A (zh) * | 2008-05-30 | 2008-10-15 | 首都师范大学 | Ct图像重建的gpu加速方法 |
| CN201378438Y (zh) * | 2009-03-04 | 2010-01-06 | 中国核电工程有限公司 | 一种252Cf中子活化核燃料棒232U富集度及均匀性的γ射线测量探测器 |
| CN102565845A (zh) * | 2010-12-31 | 2012-07-11 | 同方威视技术股份有限公司 | 利用多个探测器的伽马能谱核素识别方法 |
| CN103163548A (zh) * | 2013-03-07 | 2013-06-19 | 北京辛耕普华医疗科技有限公司 | 基于伽马相机的放射性物质探测方法及其装置和系统 |
| WO2018067208A2 (en) * | 2016-06-09 | 2018-04-12 | Phoenix Nuclear Labs Llc | System and method for performing active scanning of a nuclear fuel rod |
| CN106596496A (zh) * | 2017-01-11 | 2017-04-26 | 中核建中核燃料元件有限公司 | 一种用于燃料棒丰度检测的γ探头及探头系统 |
| CN107133462A (zh) * | 2017-04-27 | 2017-09-05 | 上海联影医疗科技有限公司 | 数据处理方法、装置及设备 |
Non-Patent Citations (2)
| Title |
|---|
| 张雷 等: "_252_Cf中子活化核燃料棒__省略_35_U富集度检测设备的软件设计", 《核电子学与探测技术》 * |
| 张雷 等: "新型无源法核燃料棒_235_U富集度检测系统研制", 《核电子学与探测技术》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111736201A (zh) * | 2020-05-25 | 2020-10-02 | 中国核电工程有限公司 | 一种核燃料棒有源检测系统及方法 |
| CN111736201B (zh) * | 2020-05-25 | 2023-11-14 | 中国核电工程有限公司 | 一种核燃料棒有源检测系统及方法 |
| CN113050151A (zh) * | 2021-03-15 | 2021-06-29 | 中国原子能科学研究院 | 一种识别放射性核素的方法以及测量系统 |
| CN113050151B (zh) * | 2021-03-15 | 2022-11-25 | 中国原子能科学研究院 | 一种识别放射性核素的方法以及测量系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109949953B (zh) | 2021-07-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106202868B (zh) | 一种获取反应堆多群核数据库中的中间共振因子的方法 | |
| CN109949953A (zh) | 核燃料棒富集度无源检测系统及富集度重构方法 | |
| Verbeke et al. | Stochastic analog neutron transport with TRIPOLI-4 and FREYA: Bayesian uncertainty quantification for neutron multiplicity counting | |
| McSpaden et al. | Preliminary NoMAD results of the MUSiC experiment | |
| US20180137945A1 (en) | Systems and methods for assaying nuclear fuel | |
| Sugaya et al. | Inverse estimation methods of unknown radioactive source for fuel debris search | |
| CN114694863B (zh) | 利用裂变产额和衰变数据计算多群缓发光子产额的方法 | |
| Jacobsson et al. | Outcomes of the JNT 1955 phase I viability study of gamma emission tomography for spent fuel verification | |
| Lee et al. | Development of the ACP safeguards neutron counter for PWR spent fuel rods | |
| Nakajima et al. | Modified conversion ratio measurement in light water-moderated UO2 lattices | |
| Kim et al. | Dual-Particle Imaging Performance of a Cs 2 LiYCl 6: Ce (CLYC)-Based Rotational Modulation Collimator (RMC) System | |
| Seo et al. | Burnup measurement of spent fuel assembly by CZT-based gamma-ray spectroscopy for input nuclear material accountability of pyroprocessing | |
| Wang et al. | A study of beryllium moderator thickness for a fission chamber with fast neutron measurements | |
| Swift | Materials and configuration from NMIS type neutron imaging and gamma spectroscopy | |
| Fang | Multi-Mode Imaging for Triso-Fueled Pebble Identification | |
| Checchia et al. | Muon tomography for spent nuclear fuel control | |
| Georgadze et al. | A remote reactor monitoring with plastic scintillation detector | |
| Swift et al. | Attributes from NMIS time coincidence, fast-neutron imaging, fission mapping, and gamma-ray spectrometry data | |
| Rana et al. | Thermal neutron flux measurements in self-serve and tray rod facilities of Dhruva reactor | |
| Fang et al. | Attenuation and Scattering Correction in Passive Gamma Emission Tomography Reconstruction | |
| Lee et al. | Performance evaluation of a multiple-scattering Compton imager for distribution of prompt gamma-rays in proton therapy | |
| Wang et al. | DRF Methodology for Reactor Pressure Vessel Dosimetry | |
| van Galen | MCNP modelling of the secondary and uncollimated radiation effects for the DIII-D Gamma-Ray Imager | |
| Andersson et al. | Feasibility Study of Using Gamma Emission Tomography for Identification of Leaking Fuel Rods in Commercial Fuel Assemblies | |
| CN116148913A (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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Li Chunmiao Inventor after: Wang Peilin Inventor after: Wang Xiaoming Inventor after: Liu Shuangquan Inventor after: Yu Dongbao Inventor after: Wei Cunfeng Inventor after: Wei Long Inventor after: Tang Hui Inventor after: Sun Xiaoli Inventor after: Liu Yantao Inventor after: Zhu Yongli Inventor before: Li Chunmiao Inventor before: Liu Shuangquan Inventor before: Wei Cunfeng Inventor before: Wei Long Inventor before: Sun Xiaoli Inventor before: Liu Yantao Inventor before: Wang Peilin Inventor before: Wang Xiaoming |
|
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20210527 Address after: 100049 Beijing Shijingshan District 19 Yuquanlu Road No. 2 hospital Applicant after: Institute of High Energy Physics Chinese Academy of Sciences Applicant after: CHINA NORTH NUCLEAR FUEL Co.,Ltd. Address before: 100049 Beijing Shijingshan District 19 Yuquanlu Road No. 2 hospital Applicant before: Institute of High Energy Physics Chinese Academy of Sciences |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |