CN111554419A - 一种基于不确定度分析的核反应堆次临界度测量方法 - Google Patents
一种基于不确定度分析的核反应堆次临界度测量方法 Download PDFInfo
- Publication number
- CN111554419A CN111554419A CN202010418695.0A CN202010418695A CN111554419A CN 111554419 A CN111554419 A CN 111554419A CN 202010418695 A CN202010418695 A CN 202010418695A CN 111554419 A CN111554419 A CN 111554419A
- Authority
- CN
- China
- Prior art keywords
- uncertainty
- analysis
- measurement
- subcritical
- neutron signal
- 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
Classifications
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21C—NUCLEAR REACTORS
- G21C17/00—Monitoring; Testing ; Maintaining
- G21C17/10—Structural combination of fuel element, control rod, reactor core, or moderator structure with sensitive instruments, e.g. for measuring radioactivity, strain
- G21C17/104—Measuring reactivity
-
- 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
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Plasma & Fusion (AREA)
- General Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明公开了一种基于不确定度分析的核反应堆次临界度测量方法,包括以下步骤:对测量中子信号的不确定度分析,基于最优化算法,获得具有最优不确定度的中子信号测量数据组,并计算其不确定度,用于后续次临界度测量分析;构建基于不确定度分析的次临界度修正模型,确定修正因子及其不确定度,用于后续次临界度分析;构建基于不确定度分析的次临界度测量模型,结合步骤1获得的中子信号数据、修正因子数据及其不确定度进行不确定度分析,得出最终次临界度测量结果及其不确定度。本发明采用基于不确定度分析的数学算法,减小中子信号测量不确定度,减小修正因子引入的不确定度,提高次临界度测量效率和准确性。
Description
技术领域
本发明涉及一种反应堆次临界度测量方法,具体涉及一种基于不确定度分析的核反应堆次临界度测量方法。
背景技术
根据法规要求,在反应堆装卸料过程、或物理启动试验前的控制棒调试过程,为确保核安全,需要测量反应堆的次临界度,在反应堆零功率物理试验时,也需要测量反应堆次临界度,为校核理论计算提供数据。
现有次临界度测量方法基于点堆模型,反应堆内有源、次临界状态下,堆内中子增殖m代后(m趋于无穷),堆内中子总数呈以下规律:
其中N为堆内中子总数,Nc为探测器中子信号,ε为与探测器特性相关的常数,S0为中子源强,l为瞬发中子寿命。
由式(2)可得:
通过计算得到参考状态的有效中子倍增因子kref,并通过测量得到参考状态的中子信号Nc,ref以及待测状态的中子信号Nc,1后,由上述两式即可得到待测状态的有效中子倍增因子k1,实现次临界度测量:
在次临界度测量过程中,每一步测量的中子信号是进行次临界度计算的关键参数,该量在测量过程中受统计涨落、高压纹波浮动、放大倍数漂移以及电磁干扰等因素的影响,测量结果存在一定不确定度,将直接作用于临界外推结果。
发明内容
针对上述技术问题,本发明提供了解决上述问题的一种基于不确定度分析的核反应堆次临界度测量方法,采用基于不确定度分析的数学算法,减小中子信号测量不确定度,减小修正因子引入的不确定度,提高次临界度测量效率和准确性。
本发明通过下述技术方案实现:
一种基于不确定度分析的核反应堆次临界度测量方法,包括以下步骤:
步骤1,对测量中子信号的不确定度分析,基于最优化算法,获得具有最优不确定度的中子信号测量数据组,并计算其不确定度,用于后续次临界度测量分析;
步骤2,构建基于不确定度分析的次临界度修正模型,确定修正因子及其不确定度,用于后续次临界度分析;
步骤3,构建基于不确定度分析的次临界度测量模型,结合步骤1获得的中子信号数据、修正因子数据及其不确定度进行不确定度分析,得出最终次临界度测量结果及其不确定度。
进一步优选,所述步骤1中,最优化算法操作包括以下步骤:
对于连续测量的n个中子信号测量数据xi,i=1,2,…,n,从中选取m个连续的中子信号测量数据,基于不确定度分析模型,计算所选测量数据组算术平均值的不确定度;
进一步优选,所述步骤2中:在次临界度测量过程中,若采用了修正算法对测量结果进行修正,则基于步骤1获得的中子信号不确定度,确定修正因子的不确定度限值范围;在计算修正所需的物理量时,控制各物理量的计算不确定度,使得修正因子的合成不确定度满足要求;
进一步优选,所述步骤2中:在计算修正所需的物理量时,各物理量的计算不确定度无法达到设定的最小值,无法使修正因子的合成不确定度满足要求,则直接计算修正因子及其不确定度,用于后续次临界度分析。
进一步优选,所述步骤2中:若开展修正计算时,尚无中子信号测量结果,则直接计算修正因子,并计算其不确定度,用于后续次临界度分析。
本发明具有如下的优点和有益效果:
1、本发明针对次临界度测量过程中子信号统计涨落、高压纹波浮动、放大倍数漂移以及电磁干扰等因素的影响,建立了一种能够减小上述影响的反应堆次临界度测量方法。采用基于不确定度分析的数学算法,减小中子信号测量不确定度,减小修正因子引入的不确定度,针对从浅至深局部或均匀引入反应性的各种次临界度测量,弥补现有反应堆次临界度测量方法的不足,提高反应堆次临界度测量的效率及准确性。
2、本发明通过在次临界度测量过程进行不确定度分析,优化测量数据,降低其不确定度的影响,从而提高测量效率和准确性。目前国内尚无其他单位开展过基于不确定度分析的次临界度测量方法开发或相关专利技术的报道。因此,需要针对自主堆芯与国外堆芯的差异性,掌握自主化的基于不确定度分析的次临界度测量方法,为堆芯次临界度的高效、准确测量提供能力。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例1
本实施例提供了一种基于不确定度分析的核反应堆次临界度测量方法,具体步骤如下所示:
步骤1,测量中子信号
采用单个或多个中子信号管测量反应堆的中子信号。在反应堆处于次临界稳态下,通过前置放大器将中子信号管测量的信号放大,经长距离传输后送到脉冲主放模块,脉冲主放经微积分放大、基线恢复、幅度甄别、整形后由缓冲隔离方式输出要求幅度和宽度的中子脉冲信号到数据采集卡,通过处理终端对数据进行分析处理,获得待测状态以及基准状态下的中子信号数据。其中,中子信号包括中子计数率等。
步骤2,基于最优化算法的中子信号测量数据处理方法
1)在测量过程中,基于不确定度分析模型,对测量的中子信号数据进行不确定度分析,可采用的不确定度分析模型有多种类型,本实施例提供如下分析方法:
在某个待测的次临界状态下,对中子信号进行测量,获得单位时间(如1秒)内的中子信号(如中子计数率),在一段时间内得到n个中子信号测得值xi,其中,i=1,2,…,n。单个测得值xi的测量不确定度s(xi)为:
在测量过程中,基于上述(1)(2)式或其它测量不确定度分析模型,对测量的中子信号数据进行不确定度分析,基于下述选取最优化中子信号测量值的数学算法,获得具有最优不确定度的中子信号测量数据组,并计算其不确定度,用于后续次临界度测量分析。
2)选取最优化中子信号测量值的数学算法:
对于连续测量的n个中子信号测量数据xi,其中i=1,2,…,n,从中选取m个连续的中子信号测量数据,基于(1)(2)式,计算所选测量数据组算术平均值的不确定度。改变所选m个连续中子信号测量数据的起点和终点,获得m=1,2,…,n时的所有测量数据组的算术平均值的不确定度从中选取具有最优不确定度的测量数据组,采用该组数据作为最优的中子信号测量值进行后续次临界度测量分析。
步骤3,构建基于不确定度分析的次临界度修正模型,包括以下模型:
1)在次临界度测量过程中,若采用了修正算法对测量结果进行修正,则基于步骤2获得的最优中子信号测量值及其不确定度,基于法规规定的标准不确定度合成算法或其他近似不确定度合成算法,根据合成不确定度限值范围,确定修正因子的不确定度限值范围。在计算修正所需的物理量时,控制各物理量的计算不确定度,使得修正因子的合成不确定度满足要求。
2)若在计算修正所需的物理量时,各物理量的计算不确定度无法达到足够低的水平,无法使修正因子的合成不确定度满足要求。则直接计算修正因子及其不确定度,用于后续次临界度分析。
3)若开展修正计算时,尚无中子信号测量结果,则直接计算修正因子,并计算其不确定度,用于后续次临界度分析。
步骤4,构建基于不确定度分析的次临界度测量模型
结合步骤1获得的中子信号数据和步骤2获得的修正因子数据及其不确定度,基于法规规定的标准不确定度合成算法或其他近似不确定度合成算法,结合下式(3)或其它的次临界度测量或修正算法,进行不确定度合成,得出最终次临界度测量结果及其不确定度。
ρ0为基准状态的次临界度,ρ1为待测状态的次临界度,N0和N1分别为基准状态和待测状态的中子计数率,C为修正因子。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种基于不确定度分析的核反应堆次临界度测量方法,其特征在于,包括以下步骤:
步骤1,对测量中子信号的不确定度分析,基于最优化算法,获得具有最优不确定度的中子信号测量数据组,并计算其不确定度,用于后续次临界度测量分析;
步骤2,构建基于不确定度分析的次临界度修正模型,确定修正因子及其不确定度,用于后续次临界度分析;
步骤3,构建基于不确定度分析的次临界度测量模型,结合步骤1获得的中子信号数据、修正因子数据及其不确定度进行不确定度分析,得出最终次临界度测量结果及其不确定度。
3.根据权利要求1所述的一种基于不确定度分析的核反应堆次临界度测量方法,其特征在于,所述步骤2中:
在次临界度测量过程中,若采用了修正算法对测量结果进行修正,则基于步骤1获得的中子信号不确定度,确定修正因子的不确定度限值范围;在计算修正所需的物理量时,控制各物理量的计算不确定度,使得修正因子的合成不确定度满足要求。
4.根据权利要求1所述的一种基于不确定度分析的核反应堆次临界度测量方法,其特征在于,所述步骤2中:
在计算修正所需的物理量时,各物理量的计算不确定度无法达到设定的最小值,无法使修正因子的合成不确定度满足要求,则直接计算修正因子及其不确定度,用于后续次临界度分析。
5.根据权利要求1所述的一种基于不确定度分析的核反应堆次临界度测量方法,其特征在于,所述步骤2中:
若开展修正计算时,尚无中子信号测量结果,则直接计算修正因子,并计算其不确定度,用于后续次临界度分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010418695.0A CN111554419B (zh) | 2020-05-18 | 2020-05-18 | 一种基于不确定度分析的核反应堆次临界度测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010418695.0A CN111554419B (zh) | 2020-05-18 | 2020-05-18 | 一种基于不确定度分析的核反应堆次临界度测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111554419A true CN111554419A (zh) | 2020-08-18 |
CN111554419B CN111554419B (zh) | 2022-11-15 |
Family
ID=72001235
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010418695.0A Active CN111554419B (zh) | 2020-05-18 | 2020-05-18 | 一种基于不确定度分析的核反应堆次临界度测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111554419B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112967825A (zh) * | 2021-03-19 | 2021-06-15 | 中国核动力研究设计院 | 一种基于修正信号不确定度分析的反应性测量方法 |
CN113161029A (zh) * | 2021-03-19 | 2021-07-23 | 中国核动力研究设计院 | 一种基于采样信号变频分析的反应性测量方法 |
CN113161028A (zh) * | 2021-03-19 | 2021-07-23 | 中国核动力研究设计院 | 一种基于修正信号优化处理的反应性测量方法 |
CN113192659A (zh) * | 2021-03-19 | 2021-07-30 | 中国核动力研究设计院 | 基于原测中子信号平均采样的核反应堆反应性测量方法 |
CN113987784A (zh) * | 2021-10-26 | 2022-01-28 | 中国核动力研究设计院 | 一种量化压水堆核设计软件包获取不确定度的方法及装置 |
CN114420328A (zh) * | 2022-01-10 | 2022-04-29 | 中国原子能科学研究院 | 反应堆次临界度的监测方法及装置 |
CN116206789A (zh) * | 2022-12-06 | 2023-06-02 | 中国核动力研究设计院 | 一种核反应堆在线次临界度监测方法和系统 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3188470A (en) * | 1963-03-05 | 1965-06-08 | Charles W Ricker | Method of determining the negative reactivity of neutronic multiplying systems |
US6181759B1 (en) * | 1999-07-23 | 2001-01-30 | Westinghouse Electric Company Llc | Method and apparatus for determining nearness to criticality of a nuclear fueled electric power generating unit |
JP2006138795A (ja) * | 2004-11-15 | 2006-06-01 | Toshiba Corp | 中性子検出器感度校正方法および燃料収納体系の未臨界度評価方法 |
CN102246242A (zh) * | 2008-12-11 | 2011-11-16 | 西屋电气有限责任公司 | 次临界反应性测量方法 |
CN102800372A (zh) * | 2012-08-17 | 2012-11-28 | 中国原子能科学研究院 | 一种反应堆反应性测量值的外源修正方法 |
ES2413567A1 (es) * | 2012-10-22 | 2013-07-16 | Centro De Investigaciones Energéticas, Medioambientales Y Tecnológicas (Ciemat) | Procedimiento de determinación de la constante de criticidad en un reactor nuclear subcrítico |
CN104036834A (zh) * | 2014-06-20 | 2014-09-10 | 中国科学院合肥物质科学研究院 | 一种次临界系统次临界度的测量方法 |
CN107689256A (zh) * | 2017-09-19 | 2018-02-13 | 中国核动力研究设计院 | 一种核反应堆次临界度测量方法 |
WO2018084940A2 (en) * | 2016-09-16 | 2018-05-11 | Elysium Industries Ltd. | Reactor control |
CN109903866A (zh) * | 2019-03-18 | 2019-06-18 | 中国原子能科学研究院 | 一种监测次临界反应堆反应性的方法 |
CN109961862A (zh) * | 2017-12-14 | 2019-07-02 | 中国核动力研究设计院 | 核电站反应堆深度次临界状态下控制棒价值动态刻度设备 |
CN109961861A (zh) * | 2017-12-14 | 2019-07-02 | 中国核动力研究设计院 | 一种核电站反应堆深度次临界状态下反应性测量仪器 |
-
2020
- 2020-05-18 CN CN202010418695.0A patent/CN111554419B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3188470A (en) * | 1963-03-05 | 1965-06-08 | Charles W Ricker | Method of determining the negative reactivity of neutronic multiplying systems |
US6181759B1 (en) * | 1999-07-23 | 2001-01-30 | Westinghouse Electric Company Llc | Method and apparatus for determining nearness to criticality of a nuclear fueled electric power generating unit |
JP2006138795A (ja) * | 2004-11-15 | 2006-06-01 | Toshiba Corp | 中性子検出器感度校正方法および燃料収納体系の未臨界度評価方法 |
CN102246242A (zh) * | 2008-12-11 | 2011-11-16 | 西屋电气有限责任公司 | 次临界反应性测量方法 |
CN102800372A (zh) * | 2012-08-17 | 2012-11-28 | 中国原子能科学研究院 | 一种反应堆反应性测量值的外源修正方法 |
ES2413567A1 (es) * | 2012-10-22 | 2013-07-16 | Centro De Investigaciones Energéticas, Medioambientales Y Tecnológicas (Ciemat) | Procedimiento de determinación de la constante de criticidad en un reactor nuclear subcrítico |
CN104036834A (zh) * | 2014-06-20 | 2014-09-10 | 中国科学院合肥物质科学研究院 | 一种次临界系统次临界度的测量方法 |
WO2018084940A2 (en) * | 2016-09-16 | 2018-05-11 | Elysium Industries Ltd. | Reactor control |
CN107689256A (zh) * | 2017-09-19 | 2018-02-13 | 中国核动力研究设计院 | 一种核反应堆次临界度测量方法 |
CN109961862A (zh) * | 2017-12-14 | 2019-07-02 | 中国核动力研究设计院 | 核电站反应堆深度次临界状态下控制棒价值动态刻度设备 |
CN109961861A (zh) * | 2017-12-14 | 2019-07-02 | 中国核动力研究设计院 | 一种核电站反应堆深度次临界状态下反应性测量仪器 |
CN109903866A (zh) * | 2019-03-18 | 2019-06-18 | 中国原子能科学研究院 | 一种监测次临界反应堆反应性的方法 |
Non-Patent Citations (4)
Title |
---|
PALMIOTTI, G. SALVATORES, M.: "Nuclear Data Target Accuracy Requirements for MA Burners", 《TRANSACTIONS OF THE AMERICAN NUCLEAR SOCIETY》 * |
万波 等: "脉冲中子源法测量铅基零功率反应堆Venus-Ⅱ次临界度", 《原子能科学技术》 * |
李泽光 等: "基于微扰的蒙特卡罗反应堆临界搜索研究", 《原子能科学技术》 * |
黄礼渊 等: "YBW临界装置技术改造后的首次临界试验和零功率物理试验", 《核电子学与探测技术》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112967825A (zh) * | 2021-03-19 | 2021-06-15 | 中国核动力研究设计院 | 一种基于修正信号不确定度分析的反应性测量方法 |
CN113161029A (zh) * | 2021-03-19 | 2021-07-23 | 中国核动力研究设计院 | 一种基于采样信号变频分析的反应性测量方法 |
CN113161028A (zh) * | 2021-03-19 | 2021-07-23 | 中国核动力研究设计院 | 一种基于修正信号优化处理的反应性测量方法 |
CN113192659A (zh) * | 2021-03-19 | 2021-07-30 | 中国核动力研究设计院 | 基于原测中子信号平均采样的核反应堆反应性测量方法 |
CN113161028B (zh) * | 2021-03-19 | 2023-01-24 | 中国核动力研究设计院 | 一种基于修正信号优化处理的反应性测量方法 |
CN113987784A (zh) * | 2021-10-26 | 2022-01-28 | 中国核动力研究设计院 | 一种量化压水堆核设计软件包获取不确定度的方法及装置 |
CN114420328A (zh) * | 2022-01-10 | 2022-04-29 | 中国原子能科学研究院 | 反应堆次临界度的监测方法及装置 |
CN114420328B (zh) * | 2022-01-10 | 2024-05-14 | 中国原子能科学研究院 | 反应堆次临界度的监测方法及装置 |
CN116206789A (zh) * | 2022-12-06 | 2023-06-02 | 中国核动力研究设计院 | 一种核反应堆在线次临界度监测方法和系统 |
CN116206789B (zh) * | 2022-12-06 | 2024-02-27 | 中国核动力研究设计院 | 一种核反应堆在线次临界度监测方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN111554419B (zh) | 2022-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111554419B (zh) | 一种基于不确定度分析的核反应堆次临界度测量方法 | |
CN101871767B (zh) | 零件形位公差检测系统及方法 | |
CN101105986B (zh) | 反应堆反应性测量方法 | |
CN106782710B (zh) | 一种反应堆测瞬发中子衰减参数的数据处理方法 | |
CN111554418B (zh) | 一种核反应堆外推临界方法 | |
CN109215822B (zh) | 一种落棒反应性测量方法 | |
CN111812703B (zh) | 一种伽玛能谱背景计数方法 | |
CN114487976B (zh) | Mcm电子式互感器校验仪溯源不确定度评定方法及系统 | |
JPH07119828B2 (ja) | 炉心反応度を測定し、妥当性を検査する方法 | |
CN110097986B (zh) | 一种核反应堆碱性水质条件下硼浓度的监测方法及其应用 | |
CN103235350B (zh) | 放射性测井仪器稳定性检测及刻度方法与装置 | |
CN101252026B (zh) | 反应堆中子源和伽马噪声测量法 | |
CN109192341B (zh) | 基于三维时空动力学的大反应性测量方法 | |
EP1927995A2 (en) | System and method for stabilizing the measurement of radioactivity | |
CN112967825B (zh) | 一种基于修正信号不确定度分析的反应性测量方法 | |
CN114420328B (zh) | 反应堆次临界度的监测方法及装置 | |
JP2018159669A (ja) | 中性子検出器等の信号のみに基づいて核分裂性物質の組成、未臨界度、遅発中性子割合、中性子世代時間、即発中性子寿命を測定する方法。 | |
CN110749919A (zh) | 一种标定核反应堆堆外探测器的方法及装置 | |
CN113936823B (zh) | 获取堆外探测器刻度系数的方法、系统、装置及存储介质 | |
EP3977487B1 (en) | System and method to determine reactivity | |
CN112016045B (zh) | 纳秒级脉冲功率计数据处理方法 | |
CN113192659B (zh) | 基于原测中子信号平均采样的核反应堆反应性测量方法 | |
CN115877425B (zh) | 一种基于ai模型的ssr格网主站确定方法及装置 | |
US9842663B2 (en) | Core monitoring system | |
CN113161029B (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 |