CN112597631B - 轴向功率偏差的超实时计算方法、电子设备及存储介质 - Google Patents

轴向功率偏差的超实时计算方法、电子设备及存储介质 Download PDF

Info

Publication number
CN112597631B
CN112597631B CN202011396054.6A CN202011396054A CN112597631B CN 112597631 B CN112597631 B CN 112597631B CN 202011396054 A CN202011396054 A CN 202011396054A CN 112597631 B CN112597631 B CN 112597631B
Authority
CN
China
Prior art keywords
real
time
neutron
super
physical model
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
CN202011396054.6A
Other languages
English (en)
Other versions
CN112597631A (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 General Nuclear Power Corp
China Nuclear Power Engineering Co Ltd
CGN Power Co Ltd
Shenzhen China Guangdong Nuclear Engineering Design Co Ltd
Original Assignee
China General Nuclear Power Corp
China Nuclear Power Engineering Co Ltd
CGN Power Co Ltd
Shenzhen China Guangdong Nuclear Engineering Design 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 General Nuclear Power Corp, China Nuclear Power Engineering Co Ltd, CGN Power Co Ltd, Shenzhen China Guangdong Nuclear Engineering Design Co Ltd filed Critical China General Nuclear Power Corp
Priority to CN202011396054.6A priority Critical patent/CN112597631B/zh
Publication of CN112597631A publication Critical patent/CN112597631A/zh
Priority to PCT/CN2021/103398 priority patent/WO2021244667A1/zh
Application granted granted Critical
Publication of CN112597631B publication Critical patent/CN112597631B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply
    • 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

Abstract

本发明涉及轴向功率偏差的超实时计算方法、电子设备及存储介质，包括：设置边界条件；根据所述边界条件确定接口数据；基于所述接口数据并结合实时数据对仿真堆芯物理模型进行初始化；在所述仿真堆芯物理模型完成初始化后，利用所述接口数据在所述仿真堆芯物理模型中进行超实时计算。本发明可预测未来半小时趋势，为核电厂操作员提供了事故情况下的轴向功率调节应对方案的预演验证平台，能够较为准确地预测事故工况下电厂反应性及对中子通量分布的影响的动态特性，预先模拟各种事故工况发展过程，预测后期存在的各种风险，让运行人员有更充足的时间降低甚至消除风险。

Description

轴向功率偏差的超实时计算方法、电子设备及存储介质
技术领域
本发明涉及核电厂堆芯轴向功率偏差的技术领域，更具体地说，涉及一种轴向功率偏差的超实时计算方法、电子设备及存储介质。
背景技术
随着计算机技术的发展及核电增长的需求，超实时仿真技术的研发已成为一种新的核电技术领域，可提供操作预警以保障核电安全运营。核电实时仿真已在培训及工程验证方面应用多年。实时仿真是指每一时间步进行仿真模型计算的时间小于或等于实时时钟，可以向仿真机用户实时显示仿真结果。实时仿真是在时间和表现上复制电站的实际过程，并使用得仿真机的“外观和感觉”与电厂相同。其关键特征是复制，为正常和异常操作以及瞬态事件重新产生电厂表面，甚至与电厂未发生事件的分析结果相匹配。
在实时仿真的基础上，结合电站信息系统的实时监测数据，已在近年开展在线仿真的研发工作。将实时监测数据作为实时仿真的高精度模型的边界条件，与电站运行状态做同步计算，如此可应用于电站系统或设备的表现做监测，并对设备寿命做预期分析，也可应用于电站瞬态或不正常操作的回放及分析，以减少未来操作错误，也可做操作建议。
超实时仿真预警技术是在线仿真的延伸，涉及一系列的技术创新，利用实时采集实际装置获取的现场数据作为初始状态输入，通过加速仿真计算，获得核电厂在未来的状态变化趋势预测，为操作员的操作决策提供有价值的参考信息，保证机组的安全稳定运行，对提升核电厂安全性与经济性有着重要的实际意义。
轴向功率偏差是核电堆芯的重要安全操作参数，必须在其安全范围内操作时，才能保证燃料的完整性。
现有在线仿真是将实时仿真模型以电站测量数据实时校准，仿真模型跟踪(和匹配)仿真范围内的真实电站工况。仿真机与电站信息系统连接，实时获取测量数据，从而与实际电站同步，进行仿真，应用于分析电站事故以及推算操作安全限值。
然而，现有的在线仿真方案的轴向功率偏差基本上在时间上只能平等于电站实时的表现，并依据统计的方式，有限的推测轴向功率偏差相关操作参数的趋势，推测精度有限。
发明内容
本发明要解决的技术问题在于，针对现有技术的上述缺陷，提供一种轴向功率偏差的超实时计算方法、电子设备及存储介质。
本发明解决其技术问题所采用的技术方案是：构造一种轴向功率偏差的超实时计算方法，包括：
设置边界条件；
根据所述边界条件确定接口数据；
基于所述接口数据并结合实时数据对仿真堆芯物理模型进行初始化；
在所述仿真堆芯物理模型完成初始化后，利用所述接口数据在所述仿真堆芯物理模型中进行超实时计算。
其中，所述利用所述接口数据在所述仿真堆芯物理模型中进行超实时计算之后还包括：
对超实时计算结果进行输出和显示。
其中，所述设置边界条件之前包括：
获取所述仿真堆芯物理模型。
其中，所述仿真堆芯物理模型通过以下方式建立：
根据预设三维方程并基于堆芯设计数据中的中子截面库，确定堆芯中子物理特性；
根据所述堆芯中子物理特性，采用三维二群的中子物理动力学模型建立所述仿真堆芯物理模型。
其中，所述预设三维方程为：X、Y和Z方向的网格化的堆芯扩散方程。
其中，所述接口数据包括：一回路热工水力模型接收冷却液的物性参数和燃料温度。
其中，所述基于所述接口数据并结合实时数据对仿真堆芯物理模型进行初始化包括：
在实时计算的条件下，基于所述一回路热工水力模型接收冷却液的物性参数和燃料温度对所述仿真堆芯物理模型进行计算；
在计算完成后，采用所述实时数据对计算结果进行校验；
所述实时数据包括中子量程测量值。
其中，所述在所述仿真堆芯物理模型完成初始化后，利用所述接口数据在所述仿真堆芯物理模型中进行超实时计算包括：
利用所述接口数据在所述仿真堆芯物理模型中计算出形状函数；
根据所述形状函数获得反应性和延迟中子；
基于所述反应性和延迟中子计算出振幅函数；
通过所述形状函数和所述振幅函数，获得中子通量；
根据所述中子通量获得衰变热热功率和中子量程。
其中，所述根据所述中子通量获得衰变热热功率和中子量程包括：
根据所述中子通量，基于与时间关联的中子扩散方程和六群缓发中子平衡方程计算所述衰变热热功率和中子量程。
其中，所述方法还包括：
根据所述中子通量计算中子毒物浓度；
依据所述中子毒物浓度更新中子截面库中的中子截面数据。
其中，所述中子毒物浓度包括：氙的浓度和钐的浓度。
其中，所述方法还包括：
在利用所述接口数据在所述仿真堆芯物理模型中进行超实时计算时：
将所述形状函数单独按照超实时仿真的倍数进行加速计算，且所述形状函数超实时仿真的时间步长与实时仿真的时间步长相同。
其中，所述实时仿真的时间步长为：250毫秒。
其中，所述方法还包括：
在利用所述接口数据在所述仿真堆芯物理模型中进行超实时计算时：
所述反应性、所述延迟中子、所述振幅函数、以及所述中子通量的计算频率为4次/秒。
其中，所述方法还包括：
所述反应性、所述延迟中子、所述振幅函数、以及所述中子通量根据超实时仿真的倍数增加时间步长。
其中，所述堆芯中子物理特性包括：堆芯几何及燃料组态正常运行的堆芯中子物理特性、堆芯几何及燃料组态预期瞬态的堆芯中子物理特性、以及堆芯几何及燃料组态事故瞬态事件的堆芯中子物理特性中的任意一种或者多种。
其中，所述边界条件包括：控制棒位置和堆芯热工水力参数。
本发明还提供一种轴向功率偏差的超实时计算装置，包括：
设置模块，用于设置边界条件；
确定模块，用于根据所述边界条件确定接口数据；
初始化模块，用于基于所述接口数据并结合实时数据对仿真堆芯物理模型进行初始化；
超实时计算模块，用于在所述仿真堆芯物理模型完成初始化后，利用所述接口数据在所述仿真堆芯物理模型中进行超实时计算。
本发明还提供一种电子设备，包括处理器和存储器，所述存储器用于存储计算机程序，所述处理器用于执行所述存储器存储的计算机程序以实现如上所述的轴向功率偏差的超实时计算方法。
本发明还提供一种存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如上所述轴向功率偏差的超实时计算方法的步骤。
实施本发明的轴向功率偏差的超实时计算方法，具有以下有益效果：包括：设置边界条件；根据所述边界条件确定接口数据；基于所述接口数据并结合实时数据对仿真堆芯物理模型进行初始化；在所述仿真堆芯物理模型完成初始化后，利用所述接口数据在所述仿真堆芯物理模型中进行超实时计算。本发明可预测未来半小时趋势，为核电厂操作员提供了事故情况下的轴向功率调节应对方案的预演验证平台，能够较为准确地预测事故工况下电厂反应性及对中子通量分布的影响的动态特性，预先模拟各种事故工况发展过程，预测后期存在的各种风险，让运行人员有更充足的时间降低甚至消除风险。
附图说明
下面将结合附图及实施例对本发明作进一步说明，附图中：
图1是本发明实施例提供的轴向功率偏差的超实时计算方法的流程示意图；
图2是本发明实施例提供的仿真堆芯物理模型计算流程图；
图3是本发明实施例提供的轴向功率偏差的超实时计算装置的结构示意图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。
为了解决目前的仿真轴向功率偏差基本上在时间上只能平等于电站实时的表现、预测精度有限的问题，本发明提供了一种轴功率偏差的超实时计算方法，该方法可以通过超实时计算达到对未来堆芯三维功率分布情况的预测，从而为核电厂操作员提供了事故情况下的轴向功率调节应对方案的预演验证平台，能够较为准确地预测事故工况下核电厂反应性及对中子通量分布的影响的动态特性，预先模拟各种事故工况发展过程，预测后期存在的各种风险，让运行人员有更充足的时间降低甚至消除风险。
具体的，参考图1，图1为本发明实施例提供的轴向功率偏差的超实时计算方法的流程示意图。
如图1所示，该轴向功率偏差的超实时计算方法包括：
步骤S101、设置边界条件。
一些实施例中，边界条件包括：控制棒位置和堆芯热工水力参数。其中，堆芯热工水力参数包括但不限于慢化剂温度、密度、以及燃料温度。
进一步地，一些实施例中，在设置边界条件之前包括：获取仿真堆芯物理模型。
其中，仿真堆芯物理模型通过以下方式建立：
根据预设三维方程并基于堆芯设计数据中的中子截面库，确定堆芯中子物理特性；根据堆芯中子物理特性，采用三维二群的中子物理动力学模型建立仿真堆芯物理模型。可选的，预设三维方程为：X、Y和Z方向的网格化的堆芯扩散方程。
一些实施例中，堆芯中子物理特性包括：堆芯几何及燃料组态正常运行的堆芯中子物理特性、堆芯几何及燃料组态预期瞬态的堆芯中子物理特性、以及堆芯几何及燃料组态事故瞬态事件的堆芯中子物理特性中的任意一种或者多种。
步骤S102、根据边界条件确定接口数据。
其中，接口数据包括：一回路热工水力模型接收冷却液的物性参数和燃料温度。具体的，在确定仿真的边界条件后，即可根据所设定的边界条件确定接口数据，并将所确定的接口数据输入仿真堆芯物理模型中。
步骤S103、基于接口数据并结合实时数据对仿真堆芯物理模型进行初始化。
一些实施例中，基于接口数据并结合实时数据对仿真堆芯物理模型进行初始化包括：在实时计算的条件下，基于一回路热工水力模型接收冷却液的物性参数和燃料温度对仿真堆芯物理模型进行计算；在计算完成后，采用实时数据对计算结果进行校验；实时数据包括中子量程测量值。
具体的，在实时计算的条件下，先利用一回路热工水力模型接收冷却液的物性参数九燃料温度对仿真堆芯物理模型进行计算，再利用实时数据对计算结果进行校验，从而达到对仿真堆芯物理模型的计算结果校验的目的，保证了仿真堆芯物理模型的准确性。
可选的，一些实施例中，一回路热工水力模型接收冷却液的物性参数包括但不限于温度、密度、空隙率、以及硼浓度等。
步骤S104、在仿真堆芯物理模型完成初始化后，利用接口数据在仿真堆芯物理模型中进行超实时计算。
一些实施例中，在仿真堆芯物理模型完成初始化后，利用接口数据在仿真堆芯物理模型中进行超实时计算包括：利用接口数据在仿真堆芯物理模型中计算出形状函数；根据形状函数获得反应性和延迟中子；基于反应性和延迟中子计算出振幅函数；通过形状函数和振幅函数，获得中子通量；根据中子通量获得衰变热热功率和中子量程。其中，仿真堆芯物理模型的计算流程图如图2所示。
具体的，在完成对仿真堆芯物理模型的初始化之后，直接调用所输入的接口数据，在仿真堆芯物理模型中逐步计算出形状函数，进而根据所计算出的形状函数得出反应性和延迟中子，接着通过反应性和延迟中子计算出振幅函数，再通过所计算出的形状函数和振幅函数得到中子通量，进而根据所得到的中子能量计算出衰变热热功率和中子能量的变化情况(即中子量程)。
可选的，一些实施例中，通过中子能量计算出衰变热热功率和中子通量的变化情况可以基于与时间关联的中子扩散方程和六群缓发中子平衡方程计算中子量程及裂变产生的热功率(衰变热热功率)。
进一步地，一些实施例中，在获得中子能量后，还可以通过中子计算中子毒物浓度，并依据中子毒物浓度更新中子截面库中的中子截面数据，将将所计算出的中子毒物浓度更新至中子截面库中的中子截面数据中，以作为下一次仿真计算。可选的，一些实施例中，中子毒物浓度包括：氙的浓度和钐的浓度。
进一步地，一些实施例中，在利用接口数据在仿真堆芯物理模型中进行超实时计算时：将形状函数单独按照超实时仿真的倍数进行加速计算，且形状函数超实时仿真的时间步长与实时仿真的时间步长相同。其中，实时仿真的时间步长为：250毫秒。在利用接口数据在仿真堆芯物理模型中进行超实时计算时：反应性、延迟中子、振幅函数、以及中子通量的计算频率为4次/秒。其中，反应性、延迟中子、振幅函数、以及中子通量根据超实时仿真的倍数增加时间步长。
具体的，在利用仿真堆芯模型进行超实时计算时，将形状函数单独按照超实时仿真的倍数(其中，超实时仿真的倍数可以根据实际计算需求选择确定，例如，可以是2倍、3倍、5倍等)做加速计算，并保持与实时仿真同样的250毫秒的时间步长。而反应性、延迟中子、振幅函数、以及中子能量则保持在每秒4次的计算频率，并依据超实时仿真的倍数增加时间步长，从而达到了实现2－10倍的超实时仿真的时间预测。
进一步地，一些实施例中，在步骤S104之后还可包括：
步骤S105、对超实时计算结果进行输出和显示。
具体的，超实时计算结果包括：衰变热功率和中子量程。其中，对超实时计算结果进行输出和显示可以采用本地方式显示，也可以采用远程方式显示。例如，可以通过本地显示终端进行输出和显示，或者也可以通过远程显示终端进行输出和显示，以达到方便查看的目的。
本发明实施例公开的轴向功率偏差的超实时计算方法可以预测未来半小时或者以上趋势，即可加速模型2-10倍或者以上，为核电厂操作员提供了事故情况下的轴向功率调节应对方案的预演验证平台，能够较为准确地预测事故工况下电厂反应性及对中子通量分布的影响的动态特性，预先模拟各种事故工况发展过程，预测后期存在的各种风险，让运行人员有更充足的时间降低甚至消除风险，同时也能保证仿真堆芯物理模型的精度和实时性。
参考图3，本发明还提供一种轴向功率偏差的超实时计算装置，其可用于实现本发明实施例公开的轴向功率偏差的超实时计算方法。
具体的，如图3所示，该轴向功率偏差的超实时计算装置可包括：
设置模块301，用于设置边界条件。
确定模块302，用于根据边界条件确定接口数据。
初始化模块303，用于基于接口数据并结合实时数据对仿真堆芯物理模型进行初始化。
超实时计算模块304，用于在仿真堆芯物理模型完成初始化后，利用接口数据在仿真堆芯物理模型中进行超实时计算。
本发明还提供一种电子设备，包括处理器和存储器，存储器用于存储计算机程序，处理器用于执行存储器存储的计算机程序以实现轴向功率偏差的超实时计算方法。
本发明还提供一种存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现轴向功率偏差的超实时计算方法的步骤。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
专业人员还可以进一步意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
以上实施例只为说明本发明的技术构思及特点，其目的在于让熟悉此项技术的人士能够了解本发明的内容并据此实施，并不能限制本发明的保护范围。凡跟本发明权利要求范围所做的均等变化与修饰，均应属于本发明权利要求的涵盖范围。

Claims (20)

1.一种轴向功率偏差的超实时计算方法，其特征在于，包括：
设置边界条件；
根据所述边界条件确定接口数据；
基于所述接口数据并结合实时数据对仿真堆芯物理模型进行初始化；
在所述仿真堆芯物理模型完成初始化后，利用所述接口数据在所述仿真堆芯物理模型中进行超实时计算。
2.根据权利要求1所述的轴向功率偏差的超实时计算方法，其特征在于，所述利用所述接口数据在所述仿真堆芯物理模型中进行超实时计算之后还包括：
对超实时计算结果进行输出和显示。
3.根据权利要求1所述的轴向功率偏差的超实时计算方法，其特征在于，所述设置边界条件之前包括：
获取所述仿真堆芯物理模型。
4.根据权利要求1所述的轴向功率偏差的超实时计算方法，其特征在于，所述仿真堆芯物理模型通过以下方式建立：
根据预设三维方程并基于堆芯设计数据中的中子截面库，确定堆芯中子物理特性；
根据所述堆芯中子物理特性，采用三维二群的中子物理动力学模型建立所述仿真堆芯物理模型。
5.根据权利要求4所述的轴向功率偏差的超实时计算方法，其特征在于，所述预设三维方程为：X、Y和Z方向的网格化的堆芯扩散方程。
6.根据权利要求1所述的轴向功率偏差的超实时计算方法，其特征在于，所述接口数据包括：一回路热工水力模型接收冷却液的物性参数和燃料温度。
7.根据权利要求6所述的轴向功率偏差的超实时计算方法，其特征在于，所述基于所述接口数据并结合实时数据对仿真堆芯物理模型进行初始化包括：
在实时计算的条件下，基于所述一回路热工水力模型接收冷却液的物性参数和燃料温度对所述仿真堆芯物理模型进行计算；
在计算完成后，采用所述实时数据对计算结果进行校验；
所述实时数据包括中子量程测量值。
8.根据权利要求1所述的轴向功率偏差的超实时计算方法，其特征在于，所述在所述仿真堆芯物理模型完成初始化后，利用所述接口数据在所述仿真堆芯物理模型中进行超实时计算包括：
利用所述接口数据在所述仿真堆芯物理模型中计算出形状函数；
根据所述形状函数获得反应性和延迟中子；
基于所述反应性和延迟中子计算出振幅函数；
通过所述形状函数和所述振幅函数，获得中子通量；
根据所述中子通量获得衰变热热功率和中子量程。
9.根据权利要求8所述的轴向功率偏差的超实时计算方法，其特征在于，所述根据所述中子通量获得衰变热热功率和中子量程包括：
根据所述中子通量，基于与时间关联的中子扩散方程和六群缓发中子平衡方程计算所述衰变热热功率和中子量程。
10.根据权利要求8所述的轴向功率偏差的超实时计算方法，其特征在于，所述方法还包括：
根据所述中子通量计算中子毒物浓度；
依据所述中子毒物浓度更新中子截面库中的中子截面数据。
11.根据权利要求10所述的轴向功率偏差的超实时计算方法，其特征在于，所述中子毒物浓度包括：氙的浓度和钐的浓度。
12.根据权利要求8所述的轴向功率偏差的超实时计算方法，其特征在于，所述方法还包括：
在利用所述接口数据在所述仿真堆芯物理模型中进行超实时计算时：
将所述形状函数单独按照超实时仿真的倍数进行加速计算，且所述形状函数超实时仿真的时间步长与实时仿真的时间步长相同。
13.根据权利要求12所述的轴向功率偏差的超实时计算方法，其特征在于，所述实时仿真的时间步长为：250毫秒。
14.根据权利要求12所述的轴向功率偏差的超实时计算方法，其特征在于，所述方法还包括：
在利用所述接口数据在所述仿真堆芯物理模型中进行超实时计算时：
所述反应性、所述延迟中子、所述振幅函数、以及所述中子通量的计算频率为4次/秒。
15.根据权利要求14所述的轴向功率偏差的超实时计算方法，其特征在于，所述方法还包括：
所述反应性、所述延迟中子、所述振幅函数、以及所述中子通量根据超实时仿真的倍数增加时间步长。
16.根据权利要求4所述的轴向功率偏差的超实时计算方法，其特征在于，所述堆芯中子物理特性包括：堆芯几何及燃料组态正常运行的堆芯中子物理特性、堆芯几何及燃料组态预期瞬态的堆芯中子物理特性、以及堆芯几何及燃料组态事故瞬态事件的堆芯中子物理特性中的任意一种或者多种。
17.根据权利要求1-16任一项所述的轴向功率偏差的超实时计算方法，其特征在于，所述边界条件包括：控制棒位置和堆芯热工水力参数。
18.一种轴向功率偏差的超实时计算装置，其特征在于，包括：
设置模块，用于设置边界条件；
确定模块，用于根据所述边界条件确定接口数据；
初始化模块，用于基于所述接口数据并结合实时数据对仿真堆芯物理模型进行初始化；
超实时计算模块，用于在所述仿真堆芯物理模型完成初始化后，利用所述接口数据在所述仿真堆芯物理模型中进行超实时计算。
19.一种电子设备，其特征在于，包括处理器和存储器，所述存储器用于存储计算机程序，所述处理器用于执行所述存储器存储的计算机程序以实现如权利要求1-17任一项所述的轴向功率偏差的超实时计算方法。
20.一种存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1-17中任意一项所述轴向功率偏差的超实时计算方法的步骤。
CN202011396054.6A 2020-12-03 2020-12-03 轴向功率偏差的超实时计算方法、电子设备及存储介质 Active CN112597631B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202011396054.6A CN112597631B (zh) 2020-12-03 2020-12-03 轴向功率偏差的超实时计算方法、电子设备及存储介质
PCT/CN2021/103398 WO2021244667A1 (zh) 2020-12-03 2021-06-30 轴向功率偏差的超实时计算方法、电子设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011396054.6A CN112597631B (zh) 2020-12-03 2020-12-03 轴向功率偏差的超实时计算方法、电子设备及存储介质

Publications (2)

Publication Number Publication Date
CN112597631A CN112597631A (zh) 2021-04-02
CN112597631B true CN112597631B (zh) 2023-08-01

Family

ID=75188013

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011396054.6A Active CN112597631B (zh) 2020-12-03 2020-12-03 轴向功率偏差的超实时计算方法、电子设备及存储介质

Country Status (2)

Country Link
CN (1) CN112597631B (zh)
WO (1) WO2021244667A1 (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112597631B (zh) * 2020-12-03 2023-08-01 中广核工程有限公司 轴向功率偏差的超实时计算方法、电子设备及存储介质
CN114330169B (zh) * 2021-12-30 2024-04-09 西安交通大学 一种小型模块化超级安全气冷堆物理热工耦合分析方法
CN116702472B (zh) * 2023-06-07 2024-02-20 西安交通大学 一种针对热管堆的堆芯核热特性数值分析方法
CN116933566B (zh) * 2023-09-14 2023-12-15 中核武汉核电运行技术股份有限公司 基于仿真系统的数据记录、数据回放与数据分析方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4774050A (en) * 1986-04-10 1988-09-27 Westinghouse Electric Corp. Axial power distribution monitor and display using outputs from ex-core detectors and thermocouples
CN104915768A (zh) * 2015-06-02 2015-09-16 中科华核电技术研究院有限公司 一种用于核电厂严重事故诊断及响应支持的方法及系统
CN107301884A (zh) * 2017-07-24 2017-10-27 哈尔滨工程大学 一种混合式核电站故障诊断方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4977529A (en) * 1973-02-23 1990-12-11 Westinghouse Electric Corp. Training simulator for a nuclear power plant
US6061412A (en) * 1995-10-05 2000-05-09 Westinghouse Electric Company Llc Nuclear reaction protection system
CN107342112B (zh) * 2016-12-01 2019-09-10 中广核工程有限公司 蒸汽发生器水位预警方法及装置
CN111897240B (zh) * 2020-08-05 2023-03-31 哈尔滨工程大学 一种基于核动力系统运行的仿真方法及系统
CN112597631B (zh) * 2020-12-03 2023-08-01 中广核工程有限公司 轴向功率偏差的超实时计算方法、电子设备及存储介质

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4774050A (en) * 1986-04-10 1988-09-27 Westinghouse Electric Corp. Axial power distribution monitor and display using outputs from ex-core detectors and thermocouples
CN104915768A (zh) * 2015-06-02 2015-09-16 中科华核电技术研究院有限公司 一种用于核电厂严重事故诊断及响应支持的方法及系统
CN107301884A (zh) * 2017-07-24 2017-10-27 哈尔滨工程大学 一种混合式核电站故障诊断方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Bayesian Long Short-Term Memory Model for Fault Early Warning of Nuclear Power Turbine";Liu Gaojun 等;《IEEE》;全文 *
"核电厂反应堆堆芯物理在线仿真系统研究";赵强;《中国博士学位论文全文数据库（电子期刊）》;全文 *

Also Published As

Publication number Publication date
CN112597631A (zh) 2021-04-02
WO2021244667A1 (zh) 2021-12-09

Similar Documents

Publication Publication Date Title
CN112597631B (zh) 轴向功率偏差的超实时计算方法、电子设备及存储介质
CN107122546B (zh) 一种压水堆稳态计算的多物理耦合方法
Prošek et al. Review of quantitative accuracy assessments with fast Fourier transform based method (FFTBM)
CN109409022B (zh) 一种用于核电堆芯物理仿真可视化建模调试测试方法
MXPA06015255A (es) Metodo para determinar imagenes para limites de operación para una operacion de reactor nuclear.
CN105975713A (zh) 一种采用源重建方法预测电子电路电磁辐射的方法
CN112800683B (zh) 基于卷积神经网络的系统短路电流水平评估方法及系统
Dykin et al. Simulation of in-core neutron noise measurements for axial void profile reconstruction in boiling water reactors
CN110109374B (zh) 液体火箭发动机推力调节系统的半实物仿真方法及装置
Kan et al. Bayesian zero-failure reliability modeling and assessment method for multiple numerical control (NC) machine tools
Yum et al. Uncertainty analyses of neutron noise simulations in a zero-power reactor
Kwon et al. The real-time functional test facility for advanced instrumentation and control in nuclear power plants
Tucker et al. The impact of fueling operations on full core uncertainty analysis in CANDU reactors
Goryunov et al. Optimization and diagnostics code for technological processes: radiochemical production simulator
Bajs et al. Development of the qualified plant nodalization for safety and operational transient analysis
CN111814295A (zh) 热仿真模型的验证方法、系统、计算机设备及存储介质
Ayo-Imoru et al. Hybrid nuclear plant simulator design requirements to enable dynamic diagnostics of plant operations
JP2696049B2 (ja) 原子炉炉心核特性模擬装置
CN114253248B (zh) 自动测试验证核电站dcs控制系统的系统及方法
Pizarro-Vallejos Use Of The Fission Matrix For Steady State And Quasi-static Kinetic Modeling Of The Treat Reactor With Temperature Feedback
KR20140023676A (ko) 노내계측기 신호 기반의 반경방향 첨두계수를 이용한 노심운전제한치감시계통의 Ｐｓｅｕｄｏ Ｈｏｔ Ｐｉｎ 출력분포 구성 방법
Haapanen et al. Validation of programmable industrial automation systems for safety critical applications in NPP's; dynamic testing
JPH04265899A (ja) 原子炉シミュレータ
Zhang et al. Burnup adaptation model in STREAM/RASTK code system
JPH0743437B2 (ja) 核燃料試験体の未臨界度決定方法

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