CN113936823A - 获取堆外探测器刻度系数的方法、系统、装置及存储介质 - Google Patents

获取堆外探测器刻度系数的方法、系统、装置及存储介质 Download PDF

Info

Publication number
CN113936823A
CN113936823A CN202111192021.4A CN202111192021A CN113936823A CN 113936823 A CN113936823 A CN 113936823A CN 202111192021 A CN202111192021 A CN 202111192021A CN 113936823 A CN113936823 A CN 113936823A
Authority
CN
China
Prior art keywords
power
detector
power distribution
core
normalized
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
CN202111192021.4A
Other languages
English (en)
Other versions
CN113936823B (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.)
Nuclear Power Institute of China
Original Assignee
Nuclear Power Institute of China
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 Nuclear Power Institute of China filed Critical Nuclear Power Institute of China
Priority to CN202111192021.4A priority Critical patent/CN113936823B/zh
Publication of CN113936823A publication Critical patent/CN113936823A/zh
Application granted granted Critical
Publication of CN113936823B publication Critical patent/CN113936823B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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/108—Measuring reactor flux
    • 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

本发明公开了一种获取堆外探测器刻度系数的方法、系统、装置及存储介质，方法包括以下步骤：获取修正系数；基于当前堆芯状态，改变堆芯功率、控制棒棒位或/和轴向氙分布以获得多组其他堆芯状态，并分别计算多组其他堆芯状态下的归一化理论功率分布
Figure DDA0003301565830000012
和堆芯功率Pr；根据修正系数、多组归一化理论功率分布
Figure DDA0003301565830000011
以及多组堆芯功率Pr计算堆外探测器刻度系数。本发明的目的在于提供一种获取堆外探测器刻度系数的方法、系统、装置及存储介质，利用一次堆内通量测量结果来修正核设计理论计算偏差，最终通过求解一个最优化问题获得堆外探测器的刻度系数，从而使获得堆外探测器刻度系数更加精确。

Description

获取堆外探测器刻度系数的方法、系统、装置及存储介质
技术领域
本发明涉及核电技术领域，尤其涉及一种获取堆外探测器刻度系数的方法、系统、装置及存储介质。
背景技术
确保核反应堆安全是核反应堆正常运行的先决条件。为了保证核反应堆安全，需要连续不断地监测核反应堆的状态。目前我国核电厂监测的主要手段之一是堆外核仪表系统。该系统能监测反应堆启堆到正常运行全过程，而监测对象正常功率运行的监测信号主要由功率量程探测器(以下简称为“堆外探测器”)给出。堆外探测器信号被处理后，可反映出反应堆的功率和堆芯轴向功率分布，从而实现对核反应堆的监测。
反映堆外探测器电流和堆芯功率与轴向功率分布偏差之间的系数称之为堆外探测器刻度系数。传统确定堆外探测器刻度系数方法是多点法：通过在不同的功率台阶上引入氙震荡，通过堆内可移动式探测器测量和堆外电流信号，从而计算探测器刻度系数。然而，该方式需要人为引入氙震荡，影响了反应堆安全性和经济性。现国内部分核电站采用一点法校刻堆外探测器刻度系数：即用理论模拟氙震荡代替真实的电厂实验。但是现有技术中，仅利用实测电流进行修正，并且在构造模拟计算时仅考虑满功率状态下，因此与实际情况存在很大的误差。
发明内容
本发明的目的在于提供一种获取堆外探测器刻度系数的方法、系统、装置及存储介质，利用一次堆内通量测量结果来修正核设计理论计算偏差，最终通过求解一个最优化问题获得堆外探测器的刻度系数，从而使获得堆外探测器刻度系数更加精确。
本发明通过下述技术方案实现：
在本申请的一个方面中，提供了一种获取堆外探测器刻度系数的方法，包括以下步骤：
S1：获取修正系数；所述修正系数为当前堆芯状态的归一化实测功率分布
Figure BDA0003301565810000011
与归一化理论功率分布
Figure BDA0003301565810000012
的比值；
S2：基于当前堆芯状态，改变堆芯功率、控制棒棒位或/和轴向氙分布以获得多组其他堆芯状态，并分别计算多组所述其他堆芯状态下的归一化理论功率分布
Figure BDA0003301565810000013
和堆芯功率Pr；
S3：根据所述修正系数、多组所述归一化理论功率分布
Figure BDA0003301565810000014
以及多组所述堆芯功率Pr计算堆外探测器刻度系数。
现国内部分核电站采用一点法校刻堆外探测器刻度系数：即用理论模拟氙震荡代替真实的电厂实验。但是现有技术中，仅利用实测电流进行修正，并且在构造模拟计算时仅考虑满功率状态下，因此与实际情况存在很大的偏差。基于此，在本申请中，提供了一种获取堆外探测器刻度系数的方法，通过在线监测的方式建模堆芯理论模型，这样理论计算的堆芯状态和实际堆芯状态将更接近；同时通过多种方式模拟构造堆芯状态，如模拟插棒、模拟氙震荡、模拟低功率到满功率等各个功率水平的状态等更多实际运行的可能状态，这样使得模拟的状态范围更广，能适应于不同的场景；此外，本方案还摒弃了传统方法中的假设条件，在最小假设的情况下，基于求解最优化问题获得最佳的堆外探测器刻度系数，从而使获得的堆外探测器刻度系数更加准确。
优选地，所述S1包括以下子步骤：
S11：获取当前堆芯状态的实测通量图，并根据所述实测通量图获取归一化实测功率分布
Figure BDA0003301565810000021
实测总功率PrM以及功率差值ΔIM，所述功率差值ΔIM为堆外探测器上部功率与堆外探测器下部功率的差值；
获取堆外探测器各段电流值
Figure BDA0003301565810000022
并根据所述各段电流值
Figure BDA0003301565810000023
获取堆外探测器上部电流
Figure BDA0003301565810000024
和堆外探测器下部电流值
Figure BDA0003301565810000025
S12：根据所述归一化实测功率分布
Figure BDA0003301565810000026
所述实测总功率PrM、所述堆外探测器上部电流
Figure BDA0003301565810000027
所述堆外探测器下部电流值
Figure BDA0003301565810000028
以及所述功率差值ΔIM进行建模，得到与当前堆芯状态对应的第一理论模型，并根据所述第一理论模型计算所述归一化理论功率分布
Figure BDA0003301565810000029
S13：根据所述归一化实测功率分布
Figure BDA00033015658100000210
与所述归一化理论功率分布
Figure BDA00033015658100000211
获取所述修正系数：
Figure BDA00033015658100000212
其中，Rj表示修正系数。
优选地，所述S2包括以下子步骤：
S21：改变所述第一理论模型中的堆芯功率、控制棒棒位或/和轴向氙分布以获得多组其他堆芯状态对应的第二理论模型；
S22：根据多组所述第二理论模型分别计算所述归一化理论功率分布
Figure BDA00033015658100000213
和所述堆芯功率Pr。
优选地，所述S3包括以下子步骤：
S31：根据所述修正系数和所述归一化理论功率分布
Figure BDA0003301565810000031
获取修正后的归一化理论功率分布P′j：
Figure BDA0003301565810000032
S32：对所述归一化理论功率分布Pj′进行归一化处理得到归一化理论功率分布
Figure BDA0003301565810000033
S33：根据所述堆芯功率Pr和所述归一化理论功率分布
Figure BDA0003301565810000034
计算修正后的功率差值ΔI；
S34：根据所述堆芯功率Pr、所述归一化理论功率分布
Figure BDA0003301565810000035
以及所述功率差值ΔI计算堆外探测器刻度系数；
Figure BDA0003301565810000036
Figure BDA0003301565810000037
Figure BDA0003301565810000038
Figure BDA0003301565810000039
Figure BDA00033015658100000310
Figure BDA00033015658100000311
其中，Ku、KL和α均为堆外探测器刻度系数，
Figure BDA00033015658100000312
为第i组的堆外探测器上部电流、
Figure BDA00033015658100000313
为第i组的堆外探测器下部电流，w为权重系数，N表示有N组数据，
Figure BDA00033015658100000314
为第i组的预测电流，Pri为第i组的堆芯功率，
Figure BDA00033015658100000315
为第i组的归一化理论功率分布，Vj为堆芯节块体积，wj,k为堆外探测器响应因子，表示节块j对探测器k的贡献。
在本申请的第二个方面中，提供了一种获取堆外探测器刻度系数的系统，包括：
获取模块，用于获取修正系数；所述修正系数为当前堆芯状态的归一化实测功率分布
Figure BDA00033015658100000316
与归一化理论功率分布
Figure BDA0003301565810000041
的比值；
处理模块，基于当前堆芯状态，改变堆芯功率、控制棒棒位或/和轴向氙分布以获得多组其他堆芯状态，并分别计算多组所述其他堆芯状态下的归一化理论功率分布
Figure BDA0003301565810000042
和堆芯功率Pr；
计算模块，用于根据所述修正系数、多组所述归一化理论功率分布
Figure BDA0003301565810000043
以及多组所述堆芯功率Pr计算堆外探测器刻度系数。
现国内部分核电站采用一点法校刻堆外探测器刻度系数：即用理论模拟氙震荡代替真实的电厂实验。但是现有技术中，仅利用实测电流进行修正，并且在构造模拟计算时仅考虑满功率状态下，因此与实际情况存在很大的偏差。基于此，在本申请中，提供了一种获取堆外探测器刻度系数的系统，通过在线监测的方式建模堆芯理论模型，这样理论计算的堆芯状态和实际堆芯状态将更接近；同时通过多种方式模拟构造堆芯状态，如模拟插棒、模拟氙震荡、模拟低功率到满功率等各个功率水平的状态等更多实际运行的可能状态，这样使得模拟的状态范围更广，能适应于不同的场景；此外，本方案还摒弃了传统方法中的假设条件，在最小假设的情况下，基于求解最优化问题获得最佳的堆外探测器刻度系数，从而使获得的堆外探测器刻度系数更加准确。
优选地，所述获取模块包括：
第一获取单元，用于获取当前堆芯状态的实测通量图，并根据所述实测通量图获取归一化实测功率分布
Figure BDA0003301565810000044
实测总功率PrM以及功率差值ΔIM，所述功率差值ΔIM为堆外探测器上部功率与堆外探测器下部功率的差值；
第二获取单元，用于获取堆外探测器各段电流值
Figure BDA0003301565810000045
并根据所述各段电流值
Figure BDA0003301565810000046
获取堆外探测器上部电流
Figure BDA0003301565810000047
和堆外探测器下部电流值
Figure BDA0003301565810000048
建模单元，用于根据所述归一化实测功率分布
Figure BDA0003301565810000049
所述实测总功率PrM、所述堆外探测器上部电流
Figure BDA00033015658100000410
所述堆外探测器下部电流值
Figure BDA00033015658100000411
以及所述功率差值ΔIM进行建模，得到与当前堆芯状态对应的第一理论模型；
第一计算单元，用于根据所述第一理论模型计算所述归一化理论功率分布
Figure BDA00033015658100000412
第二计算单元，用于所述归一化实测功率分布
Figure BDA00033015658100000413
与所述归一化理论功率分布
Figure BDA00033015658100000414
获取所述修正系数：
Figure BDA0003301565810000051
其中，Rj表示修正系数。
优选地，所述处理模块包括：
处理单元，用于改变所述第一理论模型中的堆芯功率、控制棒棒位或/和轴向氙分布以获得多组其他堆芯状态对应的第二理论模型；
第三计算单元，用于根据多组所述第二理论模型分别计算所述归一化理论功率分布
Figure BDA0003301565810000052
和所述堆芯功率Pr。
优选地，所述计算模块包括：
修正单元，用于根据所述修正系数和所述归一化理论功率分布
Figure BDA0003301565810000053
获取修正后的归一化理论功率分布P′j：
Figure BDA0003301565810000054
归一化单元，用于对所述归一化理论功率分布Pj′进行归一化处理得到归一化理论功率分布
Figure BDA0003301565810000055
第四计算单元，用于根据所述堆芯功率Pr和所述归一化理论功率分布
Figure BDA0003301565810000056
计算修正后的功率差值ΔI；
第五计算单元，用于根据所述堆芯功率Pr、所述归一化理论功率分布
Figure BDA0003301565810000057
以及所述功率差值ΔI计算堆外探测器刻度系数；
Figure BDA0003301565810000058
Figure BDA0003301565810000059
Figure BDA00033015658100000510
Figure BDA00033015658100000511
Figure BDA0003301565810000061
Figure BDA0003301565810000062
其中，Ku、KL和α均为堆外探测器刻度系数，
Figure BDA0003301565810000063
为第i组的堆外探测器上部电流、
Figure BDA0003301565810000064
为第i组的堆外探测器下部电流，w为权重系数，N表示有N组数据，
Figure BDA0003301565810000065
为第i组的预测电流，Pri为第i组的堆芯功率，
Figure BDA0003301565810000066
为第i组的归一化理论功率分布，Vj为堆芯节块体积，wj,k为堆外探测器响应因子，表示节块j对探测器k的贡献。
在本申请的第三个方面中，提供了一种电子装置，包括处理器和存储器；
存储器，用于存储处理器可执行指令；
处理器，被配置为执行如上所述的获取堆外探测器刻度系数的方法。
在本申请的第四个方面中，提供了一种计算机可读存储介质，包括存储的计算机程序，所述计算机程序运行时执行如上所述的获取堆外探测器刻度系数的方法。
本发明与现有技术相比，具有如下的优点和有益效果：
1、在堆外探测器刻度系数计算中，只需要一次堆内测量结果，不再需要人为移动控制棒改变堆芯轴向功率分布，减少刻度系数试验的时间，而且避免了堆芯额外的不安全性；
2、减少了可移动式堆内探测器的使用频率，增加了堆内探测器的寿命；
3、通过在线监测方式或采用离线跟踪方式构造与实际堆芯很接近的理论模型；
4、模拟计算了不同功率水平、不同实际运行可能会出现的堆芯状态下的功率分布情况，使得堆外探测器刻度系数适用性广泛。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部分，并不构成对本发明实施例的限定。在附图中：
图1为本发明的流程示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。
实施例1
本实施例提供了一种获取堆外探测器刻度系数的方法，如图1所示，包括以下步骤：
S1：获取修正系数，具体地，包括：
利用在线监测系统或者利用堆芯核设计程序(比如SCIENCE软件包或者中国核动力研究设计院NESTOR软件包等)离线跟踪堆芯运行状态至当前堆芯状态；本实施例所说的当前堆芯状态指的是堆芯的最新状态，从而使得后续建立的理论模型中的各燃耗分布、可燃毒物分布和实际情况尽可能的接近，以减小误差；
处理实测堆内通量测量结果(比如堆内可移动式探测器或固定式探测器测量得到的实测通量图)，得到实测总功率PrM、归一化实测功率分布
Figure BDA0003301565810000071
以及功率差值ΔIM；同时记录对应的堆外探测器各段电流值
Figure BDA0003301565810000072
并根据各段电流值
Figure BDA0003301565810000073
计算得到堆外探测器上部电流
Figure BDA0003301565810000074
和下部电流值
Figure BDA0003301565810000075
根据当前堆芯状态(即获取的归一化实测功率分布
Figure BDA0003301565810000076
实测总功率PrM、堆外探测器上部电流
Figure BDA0003301565810000077
堆外探测器下部电流值
Figure BDA0003301565810000078
以及功率差值ΔIM)，用核设计软件建模当前堆芯(即第一理论模型)，并计算出归一化理论功率分布
Figure BDA0003301565810000079
计算归一化实测功率分布与归一化理论功率分布的比值Rj，
Figure BDA00033015658100000710
得到修正系数。
值得说明的是，本实施例中的上标M表示测量值，下标j表示三维堆芯的节块编号，k表示探测器轴向编号，u表示堆外探测器上部电流之和，L表示堆外探测器下部电流之和，mc表示对应于测量状态下的计算值。
在本实施例中，在堆外探测器刻度系数计算中，只需要一次堆内测量结果，用以获取修正系数，在后续操作过程中，不再需要人为移动控制棒改变堆芯轴向功率分布，而是直接通过模拟的方式获取其他堆芯状态，从而减少操作时间，而且避免了堆芯额外的不安全性。
S2：基于当前堆芯状态，改变堆芯功率、控制棒棒位或/和轴向氙分布以获得多组其他堆芯状态，并分别计算多组其他堆芯状态下的归一化理论功率分布
Figure BDA00033015658100000711
和堆芯功率Pr；具体地：
本实施例提供了两种方式获取多组其他堆芯状态：
第一种：在核设计软件建模的第一理论模型中，改变堆芯功率、控制棒棒位和轴向氙分布的一种或多种，以获得第一组第二理论模型；在获得的第一组第二理论模型的基础上，在改变堆芯功率、控制棒棒位和轴向氙分布的一种或多种，以获得第二组第二理论模型，重复上述步骤，以获得N组第二理论模型；并且每获得一组第二理论模型，便用核设计软件计算当前第二理论模型对应的堆芯状态的归一化理论功率分布
Figure BDA00033015658100000712
和堆芯功率Pr；
第二种：在核设计软件建模的第一理论模型中，通过多次改变堆芯功率、控制棒棒位和轴向氙分布的一种或多种，且每次改变的堆芯功率、控制棒棒位和轴向氙分布的值不完全相同，以获得N组不同的第二理论模型；并且每获得一组第二理论模型，便用核设计软件计算当前第二理论模型对应的堆芯状态的归一化理论功率分布
Figure BDA0003301565810000081
和堆芯功率Pr。
在具体实施时，两种方法可单独实施，也可混合一起实施，在本实施例中，为了简化操作过程，选择第一种方法获取N组第二理论模型。
在本实施例中，在计算堆外探测器刻度系数时，不仅仅只考虑了一种状态下的堆芯状态，而是模拟计算了不同功率水平、不同实际运行可能会出现的堆芯状态下的功率分布情况，使得堆外探测器刻度系数适用性广泛，在面对不同的堆芯状态时，也可以进行测量。
S3：根据修正系数、多组归一化理论功率分布
Figure BDA0003301565810000082
以及多组堆芯功率Pr计算堆外探测器刻度系数；以下步骤为总的处理流程，在具体实施时，每一组第二理论模型都应执行下式步骤：
首先，根据修正系数和归一化理论功率分布
Figure BDA0003301565810000083
获取修正后的归一化理论功率分布P′j，
Figure BDA0003301565810000084
其次，对修正后的归一化理论功率分布P′j重新进行归一化处理得到
Figure BDA0003301565810000085
然后，根据堆芯功率Pr和归一化理论功率分布
Figure BDA0003301565810000086
计算修正后的功率差值ΔI；
每一组第二理论模型执行上述步骤后，可得到每一组的堆芯功率Pr、归一化理论功率分布
Figure BDA0003301565810000087
以及功率差值ΔI；
基于每一组的堆芯功率Pr、归一化理论功率分布
Figure BDA0003301565810000088
以及功率差值ΔI，可以得到堆外探测器刻度系数；
Figure BDA0003301565810000089
Figure BDA00033015658100000810
Figure BDA00033015658100000811
Figure BDA0003301565810000091
Figure BDA0003301565810000092
Figure BDA0003301565810000093
其中，Ku、KL和α均为堆外探测器刻度系数，
Figure BDA0003301565810000094
为第i组的堆外探测器上部电流、
Figure BDA0003301565810000095
为第i组的堆外探测器下部电流，w为权重系数，取值为正值，代表了最小化功率误差和最小化ΔI误差之间的某种平衡，N表示有N组数据，
Figure BDA0003301565810000096
为第i组的预测电流，Pri为第i组的堆芯功率，
Figure BDA0003301565810000097
为第i组的归一化理论功率分布，Vj为堆芯节块体积，wj,k为堆外探测器响应因子，表示节块j对探测器k的贡献。
值得说明的是，在实际处理时，也可以按照下式堆外探测器刻度系数；
Figure BDA0003301565810000098
但是，相较于前一种计算方式，此种方式的计算过程稍显复杂，且计算结果的稳定性也不如第一种，因此在本实施例中，选用第一种方式获取堆外探测器刻度系数。
在本实施例中，摒弃了传统方法中的假设条件(传统方法中，需要假定堆内AO和堆外AO线性)，直接通过求解极值问题获得最佳的堆外探测器刻度系数，即：通过在最小假设的情况下，基于求解最优化问题获得最佳的堆外探测器刻度系数，从而使获得的堆外探测器刻度系数更加准确。
实施例2
本实施例在实施例1的基础上提供了一种获取堆外探测器刻度系数的系统，包括：
获取模块，用于获取修正系数；修正系数为当前堆芯状态的归一化实测功率分布
Figure BDA0003301565810000099
与归一化理论功率分布
Figure BDA00033015658100000910
的比值；
具体地，本实施例中的获取模块包括：
第一获取单元，用于获取当前堆芯状态的实测通量图，并根据实测通量图获取归一化实测功率分布
Figure BDA00033015658100000911
实测总功率PrM以及功率差值ΔIM，功率差值ΔIM为堆外探测器上部功率与堆外探测器下部功率的差值；
第二获取单元，用于获取堆外探测器各段电流值
Figure BDA00033015658100000912
并根据各段电流值
Figure BDA00033015658100000913
获取堆外探测器上部电流
Figure BDA0003301565810000101
和堆外探测器下部电流值
Figure BDA0003301565810000102
建模单元，用于根据归一化实测功率分布
Figure BDA0003301565810000103
实测总功率PrM、堆外探测器上部电流
Figure BDA0003301565810000104
堆外探测器下部电流值
Figure BDA0003301565810000105
以及功率差值ΔIM进行建模，得到与当前堆芯状态对应的第一理论模型；
第一计算单元，用于根据第一理论模型计算归一化理论功率分布
Figure BDA0003301565810000106
第二计算单元，用于归一化实测功率分布
Figure BDA0003301565810000107
与归一化理论功率分布
Figure BDA0003301565810000108
获取修正系数：
Figure BDA0003301565810000109
其中，Rj表示修正系数。
处理模块，基于当前堆芯状态，改变堆芯功率、控制棒棒位或/和轴向氙分布以获得多组其他堆芯状态，并分别计算多组其他堆芯状态下的归一化理论功率分布
Figure BDA00033015658100001010
和堆芯功率Pr；
具体地，本实施例中的处理模块包括：
处理单元，用于改变第一理论模型中的堆芯功率、控制棒棒位或/和轴向氙分布以获得多组其他堆芯状态对应的第二理论模型；
第三计算单元，用于根据多组第二理论模型分别计算归一化理论功率分布
Figure BDA00033015658100001011
和堆芯功率Pr。
计算模块，用于根据修正系数、多组归一化理论功率分布
Figure BDA00033015658100001012
以及多组堆芯功率Pr计算堆外探测器刻度系数。
具体地，本实施例中的计算模块包括：
修正单元，用于根据修正系数和归一化理论功率分布
Figure BDA00033015658100001013
获取修正后的归一化理论功率分布P′j：
Figure BDA00033015658100001014
归一化单元，用于对归一化理论功率分布P′j进行归一化处理得到归一化理论功率分布
Figure BDA00033015658100001015
第四计算单元，用于根据堆芯功率Pr和归一化理论功率分布
Figure BDA00033015658100001016
计算修正后的功率差值ΔI；
第五计算单元，用于根据堆芯功率Pr、归一化理论功率分布
Figure BDA0003301565810000111
以及功率差值ΔI计算堆外探测器刻度系数；
Figure BDA0003301565810000112
Figure BDA0003301565810000113
Figure BDA0003301565810000114
Figure BDA0003301565810000115
Figure BDA0003301565810000116
Figure BDA0003301565810000117
其中，Ku、KL和α均为堆外探测器刻度系数，
Figure BDA0003301565810000118
为第i组的堆外探测器上部电流、
Figure BDA0003301565810000119
为第i组的堆外探测器下部电流，w为权重系数，N表示有N组数据，
Figure BDA00033015658100001110
为第i组的预测电流，Pri为第i组的堆芯功率，
Figure BDA00033015658100001111
为第i组的归一化理论功率分布，Vj为堆芯节块体积，wj,k为堆外探测器响应因子，表示节块j对探测器k的贡献。
实施例3
本实施例在实施例1的基础上还提供了一种电子装置，包括处理器和存储器；
存储器，用于存储处理器可执行指令；
处理器，被配置为执行如实施例1所提供的获取堆外探测器刻度系数的方法。
实施例4
本实施例在实施例1的基础上还提供了一种计算机可读存储介质，包括存储的计算机程序，计算机程序运行时执行如实施例1所提供的获取堆外探测器刻度系数的方法。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.获取堆外探测器刻度系数的方法，其特征在于，包括以下步骤：
S1：获取修正系数；所述修正系数为当前堆芯状态的归一化实测功率分布
Figure FDA0003301565800000011
与归一化理论功率分布
Figure FDA0003301565800000012
的比值；
S2：基于当前堆芯状态，改变堆芯功率、控制棒棒位或/和轴向氙分布以获得多组其他堆芯状态，并分别计算多组所述其他堆芯状态下的归一化理论功率分布
Figure FDA0003301565800000013
和堆芯功率Pr；
S3：根据所述修正系数、多组所述归一化理论功率分布
Figure FDA0003301565800000014
以及多组所述堆芯功率Pr计算堆外探测器刻度系数。
2.根据权利要求1所述的获取堆外探测器刻度系数的方法，其特征在于，所述S1包括以下子步骤：
S11：获取当前堆芯状态的实测通量图，并根据所述实测通量图获取归一化实测功率分布
Figure FDA0003301565800000015
实测总功率PrM以及功率差值ΔIM，所述功率差值ΔIM为堆外探测器上部功率与堆外探测器下部功率的差值；
获取堆外探测器各段电流值
Figure FDA0003301565800000016
并根据所述各段电流值
Figure FDA0003301565800000017
获取堆外探测器上部电流
Figure FDA0003301565800000018
和堆外探测器下部电流值
Figure FDA0003301565800000019
S12：根据所述归一化实测功率分布
Figure FDA00033015658000000110
所述实测总功率PrM、所述堆外探测器上部电流
Figure FDA00033015658000000111
所述堆外探测器下部电流值
Figure FDA00033015658000000112
以及所述功率差值ΔIM进行建模，得到与当前堆芯状态对应的第一理论模型，并根据所述第一理论模型计算所述归一化理论功率分布
Figure FDA00033015658000000113
S13：根据所述归一化实测功率分布
Figure FDA00033015658000000114
与所述归一化理论功率分布
Figure FDA00033015658000000115
获取所述修正系数：
Figure FDA00033015658000000116
其中，Rj表示修正系数。
3.根据权利要求2所述的获取堆外探测器刻度系数的方法，其特征在于，所述S2包括以下子步骤：
S21：改变所述第一理论模型中的堆芯功率、控制棒棒位或/和轴向氙分布以获得多组其他堆芯状态对应的第二理论模型；
S22：根据多组所述第二理论模型分别计算所述归一化理论功率分布
Figure FDA00033015658000000117
和所述堆芯功率Pr。
4.根据权利要求1-3中任意一项所述的获取堆外探测器刻度系数的方法，其特征在于，所述S3包括以下子步骤：
S31：根据所述修正系数和所述归一化理论功率分布
Figure FDA0003301565800000021
获取修正后的归一化理论功率分布P′j：
Figure FDA0003301565800000022
S32：对所述归一化理论功率分布P′j进行归一化处理得到归一化理论功率分布
Figure FDA0003301565800000023
S33：根据所述堆芯功率Pr和所述归一化理论功率分布
Figure FDA0003301565800000024
计算修正后的功率差值ΔI；
S34：根据所述堆芯功率Pr、所述归一化理论功率分布
Figure FDA0003301565800000025
以及所述功率差值ΔI计算堆外探测器刻度系数；
Figure FDA0003301565800000026
Figure FDA0003301565800000027
Figure FDA0003301565800000028
Figure FDA0003301565800000029
Figure FDA00033015658000000210
Figure FDA00033015658000000211
其中，Ku、KL和α均为堆外探测器刻度系数，
Figure FDA00033015658000000212
为第i组的堆外探测器上部电流、
Figure FDA00033015658000000213
为第i组的堆外探测器下部电流，w为权重系数，N表示有N组数据，
Figure FDA00033015658000000214
为第i组的预测电流，Pri为第i组的堆芯功率，
Figure FDA00033015658000000215
为第i组的归一化理论功率分布，Vj为堆芯节块体积，wj,k为堆外探测器响应因子，表示节块j对探测器k的贡献。
5.获取堆外探测器刻度系数的系统，其特征在于，包括：
获取模块，用于获取修正系数；所述修正系数为当前堆芯状态的归一化实测功率分布
Figure FDA0003301565800000031
与归一化理论功率分布
Figure FDA0003301565800000032
的比值；
处理模块，基于当前堆芯状态，改变堆芯功率、控制棒棒位或/和轴向氙分布以获得多组其他堆芯状态，并分别计算多组所述其他堆芯状态下的归一化理论功率分布
Figure FDA0003301565800000033
和堆芯功率Pr；
计算模块，用于根据所述修正系数、多组所述归一化理论功率分布
Figure FDA0003301565800000034
以及多组所述堆芯功率Pr计算堆外探测器刻度系数。
6.根据权利要求5所述的获取堆外探测器刻度系数的系统，其特征在于，所述获取模块包括：
第一获取单元，用于获取当前堆芯状态的实测通量图，并根据所述实测通量图获取归一化实测功率分布
Figure FDA0003301565800000035
实测总功率PrM以及功率差值ΔIM，所述功率差值ΔIM为堆外探测器上部功率与堆外探测器下部功率的差值；
第二获取单元，用于获取堆外探测器各段电流值
Figure FDA0003301565800000036
并根据所述各段电流值
Figure FDA0003301565800000037
获取堆外探测器上部电流
Figure FDA0003301565800000038
和堆外探测器下部电流值
Figure FDA0003301565800000039
建模单元，用于根据所述归一化实测功率分布
Figure FDA00033015658000000310
所述实测总功率PrM、所述堆外探测器上部电流
Figure FDA00033015658000000311
所述堆外探测器下部电流值
Figure FDA00033015658000000312
以及所述功率差值ΔIM进行建模，得到与当前堆芯状态对应的第一理论模型；
第一计算单元，用于根据所述第一理论模型计算所述归一化理论功率分布
Figure FDA00033015658000000313
第二计算单元，用于所述归一化实测功率分布
Figure FDA00033015658000000314
与所述归一化理论功率分布
Figure FDA00033015658000000315
获取所述修正系数：
Figure FDA00033015658000000316
其中，Rj表示修正系数。
7.根据权利要求6所述的获取堆外探测器刻度系数的系统，其特征在于，所述处理模块包括：
处理单元，用于改变所述第一理论模型中的堆芯功率、控制棒棒位或/和轴向氙分布以获得多组其他堆芯状态对应的第二理论模型；
第三计算单元，用于根据多组所述第二理论模型分别计算所述归一化理论功率分布
Figure FDA0003301565800000041
和所述堆芯功率Pr。
8.根据权利要求5-8中任意一项所述的获取堆外探测器刻度系数的系统，其特征在于，所述计算模块包括：
修正单元，用于根据所述修正系数和所述归一化理论功率分布
Figure FDA0003301565800000042
获取修正后的归一化理论功率分布P′j：
Figure FDA0003301565800000043
归一化单元，用于对所述归一化理论功率分布P′j进行归一化处理得到归一化理论功率分布
Figure FDA0003301565800000044
第四计算单元，用于根据所述堆芯功率Pr和所述归一化理论功率分布
Figure FDA0003301565800000045
计算修正后的功率差值ΔI；
第五计算单元，用于根据所述堆芯功率Pr、所述归一化理论功率分布
Figure FDA0003301565800000046
以及所述功率差值ΔI计算堆外探测器刻度系数；
Figure FDA0003301565800000047
Figure FDA0003301565800000048
Figure FDA0003301565800000049
Figure FDA00033015658000000410
Figure FDA00033015658000000411
Figure FDA00033015658000000412
其中，Ku、KL和α均为堆外探测器刻度系数，
Figure FDA00033015658000000414
为第i组的堆外探测器上部电流、
Figure FDA00033015658000000413
为第i组的堆外探测器下部电流，w为权重系数，N表示有N组数据，
Figure FDA0003301565800000051
为第i组的预测电流，Pri为第i组的堆芯功率，
Figure FDA0003301565800000052
为第i组的归一化理论功率分布，Vj为堆芯节块体积，wj,k为堆外探测器响应因子，表示节块j对探测器k的贡献。
9.电子装置，其特征在于，包括处理器和存储器；
所述存储器，用于存储所述处理器可执行指令；
所述处理器，被配置为执行如权利要求1-4中任意一项所述的获取堆外探测器刻度系数的方法。
10.计算机可读存储介质，其特征在于，包括存储的计算机程序，所述计算机程序运行时执行如权利要求1-4中任意一项所述的获取堆外探测器刻度系数的方法。
CN202111192021.4A 2021-10-13 2021-10-13 获取堆外探测器刻度系数的方法、系统、装置及存储介质 Active CN113936823B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111192021.4A CN113936823B (zh) 2021-10-13 2021-10-13 获取堆外探测器刻度系数的方法、系统、装置及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111192021.4A CN113936823B (zh) 2021-10-13 2021-10-13 获取堆外探测器刻度系数的方法、系统、装置及存储介质

Publications (2)

Publication Number Publication Date
CN113936823A true CN113936823A (zh) 2022-01-14
CN113936823B CN113936823B (zh) 2023-11-14

Family

ID=79278880

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111192021.4A Active CN113936823B (zh) 2021-10-13 2021-10-13 获取堆外探测器刻度系数的方法、系统、装置及存储介质

Country Status (1)

Country Link
CN (1) CN113936823B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115050496A (zh) * 2022-08-11 2022-09-13 西安交通大学 一种压水堆堆外探测器单点校刻方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100104059A1 (en) * 2008-10-29 2010-04-29 Korea Electric Power Corporation, A Korean Corporation Renormalization method of excore detector
CN106128529A (zh) * 2016-06-27 2016-11-16 中广核研究院有限公司 一种堆芯三维功率分布的在线测量方法
CN111489842A (zh) * 2020-04-20 2020-08-04 上海核星核电科技有限公司 一种在氙毒尚未平衡时实施压水堆堆芯功率分布测量的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100104059A1 (en) * 2008-10-29 2010-04-29 Korea Electric Power Corporation, A Korean Corporation Renormalization method of excore detector
CN106128529A (zh) * 2016-06-27 2016-11-16 中广核研究院有限公司 一种堆芯三维功率分布的在线测量方法
CN111489842A (zh) * 2020-04-20 2020-08-04 上海核星核电科技有限公司 一种在氙毒尚未平衡时实施压水堆堆芯功率分布测量的方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115050496A (zh) * 2022-08-11 2022-09-13 西安交通大学 一种压水堆堆外探测器单点校刻方法

Also Published As

Publication number Publication date
CN113936823B (zh) 2023-11-14

Similar Documents

Publication Publication Date Title
CN114266157B (zh) 实现任意时刻反应堆功率分布试验的方法、装置和设备
KR20130079340A (ko) 핵 원자로에서 엑스코어 검출기를 캘리브레이팅하는 방법
CN110739091B (zh) 一种核电厂次临界度直接测量系统及硬件平台和测量方法
CN115392141B (zh) 一种自适应的电流互感器误差评估方法
CN113936823A (zh) 获取堆外探测器刻度系数的方法、系统、装置及存储介质
CN115508770A (zh) 一种基于kl-nb算法的电能表运行状态在线评估方法
CN112821420B (zh) 一种基于XGBoost的ASFR模型中动态阻尼因子、多维频率指标的预测方法及系统
JP4918345B2 (ja) 未臨界度測定方法及び未臨界度測定用プログラム、並びに未臨界度測定装置
CN110749919A (zh) 一种标定核反应堆堆外探测器的方法及装置
CN116682997A (zh) 燃料电池发动机供氢系统控制方法、装置和电子设备
CN114420328B (zh) 反应堆次临界度的监测方法及装置
CN111797545B (zh) 一种基于实测数据的风电机组偏航折减系数计算方法
CN116206789B (zh) 一种核反应堆在线次临界度监测方法和系统
JP2000162364A (ja) 原子炉炉心性能計算装置
CN116298572A (zh) 基于新能源场站的硬件在环测试系统及其测试方法、装置
CN114896730B (zh) 一种s3d风管支架恒荷载自动分析计算的系统与方法
CN117198572B (zh) 基于pod和中子探测数据的核燃料包壳破损在线监测方法
KR20140023676A (ko) 노내계측기 신호 기반의 반경방향 첨두계수를 이용한 노심운전제한치감시계통의 Ｐｓｅｕｄｏ Ｈｏｔ Ｐｉｎ 출력분포 구성 방법
CN117153440B (zh) 轴向氙振荡试验及堆外校刻方法、装置、存储介质和终端
JPH1020072A (ja) 炉心性能計算装置
CN107887907B (zh) 一种电力系统动态状态估计多时段滚动优化方法
CN117153437A (zh) 一种cpr机组升功率物理试验方法和系统
CN117691689A (zh) 基于多策略改进秃鹰搜索算法的光伏逆变器参数辨识方法
CN113887014A (zh) 临界热流密度关系式的确定方法、可读存储介质及开发装置
CN115910396A (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