CN116467849A - 核仪表系统指示信号的量化方法、装置、设备和存储介质 - Google Patents

核仪表系统指示信号的量化方法、装置、设备和存储介质 Download PDF

Info

Publication number
CN116467849A
CN116467849A CN202310283228.5A CN202310283228A CN116467849A CN 116467849 A CN116467849 A CN 116467849A CN 202310283228 A CN202310283228 A CN 202310283228A CN 116467849 A CN116467849 A CN 116467849A
Authority
CN
China
Prior art keywords
reactor core
power distribution
parameters
condition
conditions
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.)
Pending
Application number
CN202310283228.5A
Other languages
English (en)
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 Technology Research Institute Co Ltd
CGN Power Co Ltd
Original Assignee
China General Nuclear Power Corp
China Nuclear Power Technology Research Institute Co Ltd
CGN Power 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 Technology Research Institute Co Ltd, CGN Power Co Ltd filed Critical China General Nuclear Power Corp
Priority to CN202310283228.5A priority Critical patent/CN116467849A/zh
Publication of CN116467849A publication Critical patent/CN116467849A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21C—NUCLEAR REACTORS
    • G21C17/00—Monitoring; Testing ; Maintaining
    • G21C17/001—Mechanical simulators
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D21/00—Measuring or testing not otherwise provided for
    • G01D21/02—Measuring two or more variables by means not covered by a single other subclass
    • 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/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • 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
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/08—Probabilistic or stochastic CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/10—Noise analysis or noise optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/12—Timing analysis or timing optimisation
    • 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

核仪表系统指示信号的量化方法、装置、设备和存储介质
技术领域
本申请涉及核仪表测量技术领域，特别是涉及一种核仪表系统指示信号的量化方法、装置、设备和存储介质。
背景技术
核仪表系统(Nuclear Instrumentation，RPN)是核电站的关键测量系统，其测量的主要原理为堆芯核功率、泄漏至堆外探测器的中子注量率、堆外探测器响应之间呈线性关系，而在机组运行期间，RPN堆芯核功率和堆外探测器响应之间并不是严格的线性关系，这会导致RPN的指示信号与实际信号之间存在偏差，即存在指示信号的不确定性，为使指示信号的不确定性减小，需要对指示信号的不确定性进行评估。
目前，主要是通过对RPN在长期的运行过程中进行运行数据的统计，并根据统计的运行数据来评估RPN指示信号的不确定性，或者根据检修人员或维护人员的经验值来评估指示信号的不确定性。
但是，上述评估方法存在准确度低的问题。
发明内容
基于此，有必要针对上述技术问题，提供一种能够提高核仪表系统指示信号评估准确度的核仪表系统指示信号的量化方法、装置、设备和存储介质。
第一方面，本申请提供了一种核仪表系统指示信号的量化方法。所述方法包括：
通过模拟堆芯在不同运行工况下的功率分布情况，得到堆芯在不同运行工况下的功率分布参数；
根据堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。
在其中一个实施例中，通过模拟堆芯在不同运行工况下的功率分布情况，得到堆芯在不同运行工况下的功率分布参数，包括：
构建堆芯的实际物理模型；
获取不同运行工况下的特征参数；
将不同运行工况下的特征参数输入至堆芯的实际物理模型，得到堆芯在不同运行工况下的功率分布参数。
在其中一个实施例中，不同运行工况包括正常运行工况和事故运行工况，通过模拟堆芯在不同运行工况下的功率分布情况，得到堆芯在不同运行工况下的功率分布参数，包括：
通过模拟堆芯在正常运行工况下的功率分布情况、以及堆芯在事故运行工况下的功率分布情况，得到堆芯在正常运行工况下的功率分布参数、、以及堆芯在事故运行工况下的功率分布参数。
在其中一个实施例中，根据堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化，包括：
根据堆芯在不同运行工况下的功率分布参数，确定外部探测器堆芯在不同运行工况下的指示信号；
根据堆芯外部探测器在不同运行工况下的指示信号和基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。
在其中一个实施例中，不同运行工况包括正常运行工况和事故运行工况，根据堆芯外部探测器在不同运行工况下的功率分布参数，确定堆芯在不同运行工况下的指示信号，包括：
根据堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数，确定堆芯堆芯外部探测器在正常运行工况下和事故工况下的探测电流；
根据堆芯在标定过程中的功率分布参数，确定标定参数；
根据堆芯的标定参数和堆芯外部探测器在正常运行工况下和事故工况下的探测电流，确定堆芯外部探测器在正常运行工况下和事故工况下的指示信号。
在其中一个实施例中，根据堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数，确定堆芯外部探测器在正常运行工况下和事故工况下的探测电流，包括：
获取堆芯外部探测器的响应；
根据堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数、以及堆芯外部探测器的响应，确定堆芯在正常运行工况下和事故工况下的探测电流。
在其中一个实施例中，根据堆芯在标定过程中的功率分布参数，确定标定参数，包括：
获取堆芯外部探测器的响应，以及堆芯在标定过程中的基准功率分布参数；
根据堆芯在标定过程中的基准功率分布参数和堆芯外部探测器的响应，确定堆芯在标定过程中的探测电流；
根据堆芯在标定过程中的探测电流和堆芯在标定过程中的基准物理参数，确定堆芯的标定参数。
第二方面，本申请还提供了一种核仪表系统指示信号的量化装置。所述装置包括：
获取模块，用于通过模拟堆芯在不同运行工况下的功率分布情况，得到堆芯在不同运行工况下的功率分布参数；
量化模块，用于根据堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。
第三方面，本申请还提供了一种计算机设备。所述计算机设备包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：
通过模拟堆芯在不同运行工况下的功率分布情况，得到所述堆芯在不同运行工况下的功率分布参数；
根据所述堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对所述堆芯外部探测器的指示信号进行量化。
第四方面，本申请还提供了一种计算机可读存储介质。所述计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
通过模拟堆芯在不同运行工况下的功率分布情况，得到所述堆芯在不同运行工况下的功率分布参数；
根据所述堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对所述堆芯外部探测器的指示信号进行量化。
第五方面，本申请还提供了一种计算机程序产品。所述计算机程序产品，包括计算机程序，该计算机程序被处理器执行时实现以下步骤：
通过模拟堆芯在不同运行工况下的功率分布情况，得到所述堆芯在不同运行工况下的功率分布参数；
根据所述堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对所述堆芯外部探测器的指示信号进行量化。
上述核仪表系统指示信号的量化方法、装置、设备和存储介质，通过模拟堆芯在不同运行工况下的功率分布情况，得到堆芯在不同运行工况下的功率分布参数，并根据堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。上述方法考虑了堆芯在不同运行工况下的功率分布参数，并根据堆芯在不同运行工况下的功率分布参数和基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化，与现有技术相比，本方案无需采集实际工况下的实测数据，适用于缺少实测数据的新型堆芯控制及保护信号定值设计，因此本方案的应用性较强；另外，本方案通过理论分析对堆芯外部探测器的指示信号的不确定性进行了量化，可以得到更加准确的指示信号。
附图说明
图1为一个实施例中核仪表系统指示信号的量化方法的应用环境图；
图2为一个实施例中核仪表系统指示信号的量化方法的流程示意图；
图3为图2对应的实施例中S201步骤的流程示意图；
图4为图2对应的实施例中S202步骤的流程示意图；
图5为图4对应的实施例中S401步骤的流程示意图；
图6为图5对应的实施例中S4011步骤的流程示意图；
图7为图6对应的实施例中S4012步骤的流程示意图；
图8为另一个实施例中信号的量化方法的流程示意图；
图9为一个实施例中核仪表系统指示信号的量化装置的结构框图；
图10为一个实施例中核仪表系统指示信号的量化装置的结构框图；
图11为一个实施例中核仪表系统指示信号的量化装置的结构框图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。
核仪表系统(ReactorProtection Neutronic，RPN)是核电站的关键测量系统，其测量的主要原理为堆芯核功率、泄漏至堆外探测器的中子注量率、堆外探测器响应之间呈线性关系。然而，在机组运行期间，RPN堆芯核功率和堆外探测器响应之间并不是严格的线性关系，这会导致RPN的指示信号与实际监测物理量之间存在偏差，即存在指示信号的不确定性。
另外，由于堆芯独特的物理设计以及堆外探测器标定方式，实际堆外探测器响应会受到功率分布的影响，而机组在运行期间，燃耗效应、氙振荡及控制棒棒位变化等因素会导致堆内功率分布发生变化，这一变化将使得堆外探测器响应-堆芯指示功率之间不是严格的线性关系，这也会造成RPN指示信号与堆芯实际信号出现偏差。
而RPN的指示信号通常用于堆芯保护及堆芯控制，若RPN指示信号与真实信号的偏差评估不合理，可能导致机组误触发控制信号及保护信号，影响机组的可用率；也可能造成事故分析考虑的保守性不足，发生事故情况下无法及时触发相应的保护。
目前，主要是通过对RPN在长期的运行过程中进行运行数据的统计，并根据统计的运行数据来评估RPN指示信号的不确定性，或者根据检修人员或维护人员的经验值来评估指示信号的不确定性。然而，上述评估RPN指示信号不确定性的准确度较低。本申请旨在解决该问题。
在上述介绍完本申请实施例提供的核仪表系统指示信号的量化方法的背景技术之后，下面，将对本申请实施例提供的核仪表系统指示信号的量化方法所涉及到的实施环境进行简要说明。本申请实施例提供的核仪表系统指示信号的量化方法，可以应用于如图1所示的计算机设备中。该计算机设备包括通过系统总线连接的处理器、存储器，该存储器中存储有计算机程序，处理器执行该计算机程序时可以执行下述方法实施例的步骤。可选的，该计算机设备还可以包括输入/输出接口。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器，该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的输入/输出接口用于与外部的终端通过网络连接通信。可选的，该计算机设备可以是服务器，可以是个人计算机，还可以是个人数字助理，还可以是其他的终端设备，例如平板电脑、手机等等，还可以是云端或者远程服务器，本申请实施例对计算机设备的具体形式并不做限定。
本领域技术人员可以理解，图1中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
上述介绍了本申请实施例提供的核仪表系统指示信号的量化方法的应用场景后，下面重点介绍信号的量化过程。
在一个实施例中，如图2所示，提供了一种核仪表系统指示信号的量化方法，以该方法应用于图1中的计算机设备为例进行说明，包括以下步骤：
S201、通过模拟堆芯在不同运行工况下的功率分布情况，得到堆芯在不同运行工况下的功率分布参数。
其中，不同运行工况可以包括但不限于正常运行工况和事故运行工况，以及理想运行工况等。不同运行工况下的功率分布情况可以包括但不限于不同运行工况下的扰动装载情况、燃耗、棒位、氙振荡时间等。不同运行工况下的功率分布情况不同。
其中，物理堆芯模型是针对数值反应堆开发的具备高保真的稳态计算程序。不同运行工况下的功率分布参数不同。
本申请实施例中，可以预先设计一个实验使用的物理堆芯模型，并使该物理堆芯模型在不同运行工况下运行，从而分析该物理堆芯模型在不同运行工况下的功率分布情况，从而得到该物理堆芯模型在不同运行工况下的功率分布参数；可选的，也可以预先在仿真应用上构建仿真堆芯模型，并使该仿真堆芯模型在不同运行工况下运行，从而分析该仿真堆芯模型在不同运行工况下的功率分布情况，从而得到该仿真堆芯模型在不同运行工况下的功率分布参数。
在需要对核仪表系统指示信号的不确定性进行评估时，计算机设备可以通过上述物理堆芯模型或仿真堆芯模型分析堆芯在不同运行工况下的功率分布情况，得到堆芯在不同运行工况下的功率分布参数。需要说明的是，物理堆芯模型或仿真堆芯模型均是实际堆芯的等效模型，能够表征实际的堆芯状态。
S202、根据堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。
其中，基准功率分布参数可以是计算机设备中预设的堆芯在正常工况下的标准功率分布参数和事故工况下的标准功率分布参数。堆芯外部探测器的指示信号可以是轴向功率偏差参数、核功率等。
本申请实施例中，在上述S201获取到堆芯在不同运行工况下的功率分布参数后，计算机设备可以对堆芯在不同运行工况下的功率分布参数和计算机设备中预设的对应的基准功率分布参数进行比较，得到比较结果，该比较结果可以反应堆芯在不同运行工况下的功率分布参数与计算机设备中预设的对应的基准功率分布参数的接近程度。比如，若堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数比较接近，说明堆芯外部探测器的指示信号符合当前工况下的真实信号；若堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数相差较大，说明堆芯外部探测器的指示信号不符合当前工况下的真实信号。因此，计算机设备可以根据比较结果基于非参数统计分析方法确定堆芯外部探测器在不同运行工况下的功率分布参数与对应的基准功率分布参数不确定性程度，即基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。
本申请实施例提供的核仪表系统指示信号的量化方法，通过模拟堆芯在不同运行工况下的功率分布情况，得到堆芯在不同运行工况下的功率分布参数，并根据堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。上述方法考虑了堆芯在不同运行工况下的功率分布参数，并根据堆芯在不同运行工况下的功率分布参数和基准功率分布参数对堆芯外部探测器的指示信号进行量化，与现有技术相比，本方案无需采集实际工况下的实测数据，适用于缺少实测数据的新型堆芯控制及保护信号定值设计，因此本方案的应用性较强；另外，本方案通过理论分析对堆芯外部探测器的指示信号的不确定性进行了量化，可以得到更加准确的指示信号。
在一个实施例中，在图2所示实施例的基础上，可以对通过模拟堆芯在不同运行工况下的功率分布情况，得到堆芯在不同运行工况下的功率分布参数的过程进行描述，如图3所示，上述S201“通过模拟堆芯在不同运行工况下的功率分布情况，得到堆芯在不同运行工况下的功率分布参数”可以包括以下步骤：
S301、构建堆芯的实际物理模型。
其中，堆芯的实际物理模型可以是反应堆模型，或者，也可以是上述实施例中提到的仿真堆芯模型，堆芯的实际物理模型用于根据不同运行工况下的特征参数得到堆芯在不同工况下的功率分布参数。可选的，堆芯的实际物理模型可以基于实际堆芯在不同运行工况下的数据样本和预先构建的初始堆芯模型训练得到。
本申请实施例中，在需要对堆芯外部探测器指示信号的不确定性进行评估时，计算机设备可以根据堆芯的几何参数或物理参数构建相应的堆芯的实际物理模型；可选的，也可以预先构建初始堆芯模型，再获取实际堆芯在不同运行工况下的数据样本，并基于这些样本数据训练上述初始堆芯模型得到堆芯的物理模型。
S302、获取不同运行工况下的特征参数。
其中，运行工况下的特征参数的类别可以分为扰动不同循环燃耗、不同堆芯装载方案、不同控制棒棒况、不同功率水平，不同氙振荡时间；其中，扰动不同燃耗覆盖从寿期初到寿期末典型燃耗；不同装载方案包括：电厂运行可能出现的燃料管理典型循环方案；不同插棒工况包括：控制棒在顶部、最佳控制范围的典型棒况及插入限制位置；不同氙振荡时间包括：瞬时氙、氙峰、平衡氙等典型氙特征工况。
本申请实施例中，可以通过对上述各类别中提取上述各类别中的任意一项并进行排列组合，得到一种工况下的特征参数；进一步地，遍历提取各类别中所有条件进行组合，并进行运行限制条件的筛选，得到覆盖运行限制条件范围内不同运行工况下的特征参数。示例性的，对扰动BLX、C01循环装载方案、R棒在典型棒位ARO，以及GN棒在典型棒位ARO进行组合，得到一种工况下的特征参数；进一步地，对扰动6000MWd/tU、C01循环装载方案、R棒在典型棒位ARO，以及GN棒在典型棒位ARO进行组合，得到另一种工况下的特征参数；依次遍历所有特征参数的类别，得到各种工况下的特征参数，再对运行限制条件的筛选，得到覆盖运行限制条件范围内不同运行工况下的特征参数
S303、将不同运行工况下的特征参数输入至堆芯的实际物理模型，得到堆芯在不同运行工况下的功率分布参数。
本申请实施例中，计算机设备可以将上述S302中不同运行工况下的特征参数输入至构建好的堆芯的实际物理模型进行计算，得到计算结果，该计算结果即可表征堆芯在不同运行工况下的功率分布参数。
进一步地，通过模拟堆芯在正常运行工况下的功率分布情况、以及堆芯在事故运行工况下的功率分布情况，得到堆芯在正常运行工况下的功率分布参数、以及堆芯在事故运行工况下的功率分布参数。
示例性的，堆芯在正常运行工况下的功率分布情况可以是扰动BLX、C01循环装载方案、R棒在典型棒位ARO，以及GN棒在典型棒位ARO进行组合得到的特征参数，将这些参数输入至堆芯的实际物理模型中进行计算，得到堆芯在正常运行工况下的功率分布参数。之后，以堆芯在正常运行工况下的功率分布情况为基础，模拟事故进程即得到堆芯在事故运行工况下的功率分布参数。
本申请实施例提供的获取不同运行工况下的功率分布参数的方法，考虑了正常运行工况下、事故运行工况下以及标定过程中的功率分布情况，并基于不同工况下的功率分布情况得到不同运行工况下的功率分布参数，使得考虑的运行工况更加全面，进而使得后续对堆芯外部探测器的指示信号不确定性的量化更加准确。
在一个实施例中，在图2所示实施例的基础上，可以对根据堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化的过程进行描述，如图4所示，上述S202“根据堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化”可以包括以下步骤：
S401、根据堆芯在不同运行工况下的功率分布参数，确定堆芯外部探测器在不同运行工况下的指示信号。
其中，堆芯在不同运行工况下的指示信号可以包括但不限于核功率Pr和轴向功率偏差信号ΔI。
本申请实施例中，在计算机设备得到堆芯在不同运行工况下的功率分布参数后，可以将得到的堆芯在不同运行工况下的功率分布参数分别进行运算，得到运算结果，该运算结果可以表征堆芯外部探测器在不同运行工况下的指示信号。示例性的，可以将得到的堆芯在不同运行工况下的功率分布参数代入预设的物理模型中进行数值运算，得到堆芯外部探测器在不同运行工况下的指示信号；又一示例性的，还可以将得到的堆芯在不同运行工况下的功率分布参数与堆芯的其他参数进行数值运算，得到堆芯外部探测器在不同运行工况下的指示信号。
本申请实施例中，可以根据堆芯在不同运行工况下的功率分布参数与堆芯相关的物理模型，来确定堆芯外部探测器在不同运行工况下的指示信号。下面提供一种确定堆芯外部探测器在不同运行工况下的指示信号的方法。
示例性地，如图5所示，不同运行工况包括正常运行工况和事故运行工况，上述S401“根据堆芯在不同运行工况下的功率分布参数，确定堆芯外部探测器在不同运行工况下的指示信号”，包括：
S4011、根据堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数，确定堆芯外部探测器在正常运行工况下和事故工况下的探测电流。
其中，探测电流可以是探测器的响应电流。
本申请实施例中，在计算机设备得到堆芯在正常运行工况下和事故运行工况下的功率分布参数后，可以将得到的堆芯在正常运行工况下和事故运行工况下的功率分布参数分别进行运算，得到运算结果，该运算结果可以分别表征堆芯外部探测器在正常运行工况下和事故工况下的探测电流。示例性的，可以将得到的堆芯在正常运行工况下的功率分布参数和事故工况下的功率分布参数代入预设的电流模型中进行数值运算，分别得到堆芯外部探测器在正常运行工况下和事故工况下的探测电流；又一示例性的，还可以将得到的堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数与堆芯的其他参数进行数值运算，得到堆芯外部探测器在正常运行工况下和事故工况下的探测电流。
本申请实施例中，可以根据堆芯在不同运行工况下的功率分布参数与堆芯外部探测器的响应，来确定堆芯外部探测器在正常运行工况下和事故工况下的探测电流。下面提供一种确定堆芯外部探测器在正常运行工况下和事故工况下的探测电流的方法。
示例性地，如图6所示，上述S4011“根据堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数，确定堆芯外部探测器在正常运行工况下和事故工况下的探测电流”，包括：
S40111、获取堆芯外部探测器的响应。
其中，堆芯外部探测器的响应可以是堆芯中的一个中子散经过散射碰撞后到达堆芯外探测器处的概率，一般情况下，堆芯(x,y,z)节块的堆芯外部探测器的响应可以用W(x,y,z)表示。
本申请实施例中，计算机设备可以在预设的模型中直接获取堆芯外部探测器的响应。
S40112、根据堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数、以及堆芯外部探测器的响应，确定堆芯外部探测器在正常运行工况下和事故工况下的探测电流。
本申请实施例中，计算机设备可以对上述步骤中获取的堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数、以及堆芯外部探测器的响应进行数值运算，得到数值运算结果，该数值运算结果可以表征堆芯外部探测器在正常运行工况下和事故工况下的探测电流。示例性的，堆芯外部探测器在不同运行工况下的探测电流可以用以下公式(1)表示：
I＝C×∑x＝1,2,…,nR(x,y,z) (1)；
y＝1,2,…,n
z＝1,2,…,n
其中，I为堆芯外部探测器在不同运行工况下的探测电流，C为与堆芯外部探测器灵敏度相关的系数，该系数一般为常数，R(x,y,z)为堆芯中任一位置(x,y,z)节块内裂变产生的中子在探测器内产生的反应率。需要说明的是，不同工况下的探测电流根据节块内裂变产生的中子在探测器内产生的反应率的变化而变化。
进一步地，公式(1)中的堆芯中任一位置(x,y,z)节块内裂变产生的中子在探测器内产生的反应率R(x,y,z)可以用以下公式(2)表示：
R(x,y,z)＝w(x,y,z)×[v(x,y,z)×∑f(x,y,z)×Φ(x,y,z)] (2)；
其中，W(x,y,z)是堆芯中任意位置(x,y,z)节块的中子散射概率，表征一个中子经过散射碰撞后到达堆芯外探测器处的概率，可通过计算机设备在反应堆模型中模拟得到。[v(x,y,z)×∑f(x,y,z)×Φ(x,y,z)]是堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数，表征堆芯任一位置(x,y,z)节块处的裂变中子产生率，与堆芯状态及功率分布相关，其中的Φ(x,y,z)可通过确定论程序计算得出的。
S4012、根据堆芯在标定过程中的功率分布参数，确定标定参数。
其中，堆芯在标定过程中的功率分布参数可以是堆芯在试运行过程中理想运行(正常运行)的功率分布参数；标定参数可以是标定工况下(满功率工况下)测得的堆芯内的指示信号(Pr和ΔI)和堆外探测器响应计算出的信号(Pr'和ΔI')之间的偏差。
探测器电流的标定系数，一般情况下包括：α、上部探测器电流的标定系数KU及下部探测器电流的标定系数KL。
本申请实施例中，计算机设备可以对上述步骤中确定的堆芯在标定过程中的功率分布参数进行数值运算，得到标定参数。示例性的，可以先将上述步骤中确定的堆芯在标定过程中的功率分布参数代入上述公式(1)和(2)中，得到堆芯在标定情况下的探测电流，再对堆芯在标定情况下的探测电流进行数值运算，得到标定参数。
本申请实施例中，可以根据堆芯在不同运行工况下的功率分布参数与堆芯相关的散射概率，来确定标定参数。下面提供一种确定标定参数的方法。
示例性地，如图7所示，上述S4012“根据堆芯在标定过程中的功率分布参数，确定标定参数”，包括：
S40121、获取堆芯外部探测器的响应，以及堆芯在标定过程中的基准功率分布参数。
本申请实施例中，计算机设备可以在预设的模型中直接获取堆芯外部探测器的响应，其中预设的模型可以是反应堆模型；也可以获取计算机设备中预设的堆芯在标定过程中的基准功率分布参数，即堆芯在试运行过程中理想的功率分布参数。
S40122、根据堆芯在标定过程中的基准功率分布参数和堆芯外部探测器的响应，确定堆芯在标定过程中的探测电流。
其中，堆芯在标定过程中的探测电流可以是堆芯在试运行过程中理想的探测电流。
本申请实施例中，可以对堆芯在标定过程中的基准功率分布参数和堆芯外部探测器的响应进行运算，得到运算结果，该运算结果表征堆芯在标定过程中的探测电流。示例性的，可以将堆芯在标定过程中的基准功率分布参数和堆芯外部探测器的响应代入至上述公式(1)和(2)中，得到堆芯在标定过程中的探测电流。
S40123、根据堆芯在标定过程中的探测电流和堆芯在标定过程中的基准物理参数，确定标定参数。
本申请实施例中，计算机设备可以先获取堆芯在标定过程中的指示信号，该指示信号表征堆芯在试运行过程中理想的核功率以及轴向功率偏差参数，再根据堆芯在标定过程中的探测电流和堆芯在标定过程中的基准物理参数，确定标定参数。
示例性的，堆芯的标定参数(核功率)可以用以下公式(3)表示：
Pr＝KU×IU+KL×IL (3)；
其中，Pr是堆芯的核功率，KU、KL分别为上部探测器电流的标定系数和下部探测器电流的标定系数，IU、IL分别为上部探测电流和下部探测电流。
示例性的，堆芯的标定参数(轴向功率偏差参数)可以用以下公式(4)表示：
ΔI＝α(KU×IU-KL×IL) (4)；
其中，ΔI是堆芯的轴向功率偏差参数，KU、KL分别为上部探测器电流的标定系数和下部探测器电流的标定系数，IU、IL分别为上部探测电流和下部探测电流，α为探测器电流的一般标定系数。
S4013、根据标定参数和堆芯外部探测器在正常运行工况下和事故工况下的探测电流，确定堆芯外部探测器在正常运行工况下和事故工况下的指示信号。
本申请实施例中，在上述得到标定参数，以及堆芯在正常运行工况下和事故工况下的探测电流后，可以将标定参数，以及堆芯在正常运行工况下和事故工况下的探测电流代入预设的核物理模型中，得到堆芯外部探测器在正常运行工况下和事故工况下的指示信号。示例性的，可以将标定参数，以及堆芯在正常运行工况下和事故工况下的探测电流代入上述公式(3)和(4)中，得到堆芯外部探测器在正常运行工况下和事故工况下的指示信号。
S402、根据堆芯外部探测器在不同运行工况下的指示信号和基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。
其中，基准功率分布参数可以通过是计算机设备中预设确定论程序的反应堆模型获得，功率分布参数与探测器响应关系模型可以通过蒙卡程序获得。
本申请实施例中，可以先将堆芯在不同运行工况下的基准功率分布参数分别代入上述公式(1)、(2)、(3)和(4)中，得到堆芯在不同运行工况下的基准指示信号，再对堆芯外部探测器在不同运行工况下的指示信号和堆芯在不同运行工况下的基准指示信号进行数值运算，得到数值运算结果，该数值运算结果的大小表征堆芯外部探测器的指示信号的准确度。可选的，可以先对堆芯外部探测器在不同运行工况下的指示信号和堆芯在不同运行工况下的基准指示信号作差，再将堆芯外部探测器在不同运行工况下的指示信号与上述差作比，得到第一运算结果，同时对堆芯外部探测器在不同运行工况下的指示信号和堆芯在不同运行工况下的基准指示信号作差，再将堆芯在不同运行工况下的基准指示信号与上述差作比，得到第二运算结果，再将第一运算结果与第二运算结果进行数值比较，数值越小，说明堆芯外部探测器指示信号的准确度较高，数值越大，则说明堆芯外部探测器指示信号的准确度较低。
本申请实施例中，提供了一种对堆芯外部探测器的指示信号进行量化的方法，即对堆芯外部探测器的指示信号采用WILKS非参数统计分析方法进行不确定性统计。示例性的，根据得到的堆芯外部探测器指示信号的Pr{Pr1,PrI2,…PrN}/ΔI{ΔI1,ΔI2,…ΔIN}，与堆芯的基准指示信号Pr/ΔI做差，得到堆芯外部探测器的指示信号与堆芯的基准指示信号的偏差{x1,x2,…xN}。其中，{x1,x2,…xN}可以认为是N次独立的、与堆芯状态相关的随机函数x的输出，其概率分布函数f(x)连续。再将{x1,x2,…xN}按照从小到大的方式进行排列，形成抽样结果，同时定义x0＝-∞，xN+1＝+∞。
在WILKS非参数统计分析方法中，对于概率份额γ，由Wilks非参数统计公式可以确定容忍区间(L,U)，使得在容忍区间(L,U)内概率的置信水平β用以下公式(5)表示为：
其中，N为抽样样本数；r和s为阶数，分别代表了样本置信区间的下限和上限，0≤r＜s≤N，L＝y(r)，U＝y(s)；
若堆芯外部探测器指示信号的轴向功率偏差参数ΔI落在所要求的范围内，就认为是堆芯运行是安全的，不会触发反应堆保护系统的保护信号。表1给出了由式(5)计算得出的满足95％置信度下置信区间所需计算次数。在给定堆芯外部探测器指示信号测量偏差的样本数n、95％置信水平和95％概率份额的情况下，找到样本数n在表1右侧数据栏中对应的位置N(r)<n<N(r+1)，根据表1确定该组样本点满足两个95％的取值阶数r，那么该样本点满足95％置信度下的置信区间为(r+1,N-r-1)。该区间表明，在考虑95％的置信区间后，堆芯外部探测器指示信号的不确定性范围为(r+1,N-r-1)。示例性地，在三代核电机组典型循环的分析中，显示信号核功率的最大正偏差为5.08％，最大负偏差为-1.87％；轴向功率偏差信号的最大正偏差为2.64％，最大负偏差为-1.23％。考虑95％的置信区间后，核功率的不确定性范围为(-1.05，3.85)，轴向功率偏差信号的不确定性范围为(-0.94，1.97)。
表1满足95％置信度下置信区间所需的计算次数
本申请实施例提供的核仪表系统指示信号的量化方法，考虑了堆芯在不同运行工况下的功率分布参数，并根据堆芯在不同运行工况下的功率分布参数和基准功率分布参数对堆芯外部探测器的指示信号进行量化，与现有技术相比，本方案无需采集实际工况下的实测数据，即可通过理论分析对堆芯外部探测器指示信号的不确定性进行量化，以适用于缺少实测数据的新型堆芯控制及保护信号定值设计；另外，本方法采用WILKS非参数分布统计方法，对数据的整体分布类型没有要求，也不依赖指示信号与基准指示信号的偏差，只需要根据样本数和置信水平即可确认置信区间，从而使得确定出的堆芯外部探测器指示信号进不确定性的准确度更高。
综合上述所有实施例所述方法，图8为本申请实施例提供的一种核仪表系统指示信号的量化方法的流程图，如图8所示，该方法可以包括以下步骤：
S10、构建堆芯的实际物理模型；
S11、获取不同运行工况下的特征参数；
S12、将不同运行工况下的特征参数输入至堆芯的实际物理模型，得到堆芯在不同运行工况下的功率分布参数；
S13、获取堆芯外部探测器的响应；
S14、根据堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数、以及堆芯外部探测器的响应，确定堆芯外部探测器在正常运行工况下和事故工况下的探测电流；
S15、获取堆芯外部探测器的响应，以及堆芯在标定过程中的基准功率分布参数；
S16、根据堆芯在标定过程中的基准功率分布参数和堆芯外部探测器的响应，确定堆芯在标定过程中的探测电流；
S17、根据堆芯在标定过程中的探测电流和堆芯在标定过程中的基准物理参数，确定标定参数；
S18、根据标定参数和堆芯外部探测器在正常运行工况下和事故工况下的探测电流，确定堆芯外部探测器在正常运行工况下和事故工况下的指示信号；
S19、根据堆芯外部探测器在不同运行工况下的指示信号和基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。
本申请实施例提供的核仪表系统指示信号的量化方法，通过模拟堆芯在不同运行工况下的功率分布情况，得到堆芯在不同运行工况下的功率分布参数，并根据堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。上述方法考虑了堆芯在不同运行工况下的功率分布参数，并根据堆芯在不同运行工况下的功率分布参数和基准功率分布参数对堆芯外部探测器的指示信号进行量化，与现有技术相比，本方案无需采集实际工况下的实测数据，适用于缺少实测数据的新型堆芯控制及保护信号定值设计，因此本方案的应用性较强；另外，本方案通过理论分析对堆芯外部探测器的指示信号的不确定性进行了量化，可以得到更加准确的指示信号。
应该理解的是，虽然如上所述的各实施例所涉及的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，如上所述的各实施例所涉及的流程图中的至少一部分步骤可以包括多个步骤或者多个阶段，这些步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤中的步骤或者阶段的至少一部分轮流或者交替地执行。
基于同样的发明构思，本申请实施例还提供了一种用于实现上述所涉及的核仪表系统指示信号的量化方法的信号的量化装置。该装置所提供的解决问题的实现方案与上述方法中所记载的实现方案相似，故下面所提供的一个或多个信号的量化装置实施例中的具体限定可以参见上文中对于核仪表系统指示信号的量化方法的限定，在此不再赘述。
在一个实施例中，如图9所示，提供了一种核仪表系统指示信号的量化装置，包括：模拟模块10、量化模块20，其中：
模拟模块10，用于通过模拟堆芯在不同运行工况下的功率分布情况，得到堆芯在不同运行工况下的功率分布参数。
量化模块20，用于根据堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。
在一个实施例中，上述模拟模块10，包括构建单元100、获取单元101和输入单元102，其中：
构建单元100，具体用于构建堆芯的实际物理模型。
获取单元101，具体用于获取不同运行工况下的特征参数。
输入单元102，具体用于将不同运行工况下的特征参数输入至堆芯的实际物理模型，得到堆芯在不同运行工况下的功率分布参数。
在一个实施例中，不同运行工况包括正常运行工况和事故运行工况，上述模拟模块10，包括：通过模拟堆芯在正常运行工况下的功率分布情况、以及堆芯在事故运行工况下的功率分布情况，得到堆芯在正常运行工况下的功率分布参数、以及堆芯在事故运行工况下的功率分布参数。
在一个实施例中，如图11所示，量化模块20，包括：确定单元200和量化单元201，其中：
确定单元200，具体用于根据堆芯在不同运行工况下的功率分布参数，确定堆芯外部探测器在不同运行工况下的指示信号。
量化单元201，具体用于根据堆芯外部探测器在不同运行工况下的指示信号和基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。
在一个实施例中，不同运行工况包括正常运行工况和事故运行工况，确定单元200，具体用于根据堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数，确定堆芯外部探测器在正常运行工况下和事故工况下的探测电流；根据堆芯在标定过程中的功率分布参数，确定标定参数；根据标定参数和堆芯外部探测器在正常运行工况下和事故工况下的探测电流，确定外部探测器堆芯在正常运行工况下和事故工况下的指示信号。
在一个实施例中，确定单元200，具体用于获取堆芯外部探测器的响应；根据堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数、以及堆芯外部探测器的响应，确定堆芯外部探测器在正常运行工况下和事故工况下的探测电流。
在一个实施例中，确定单元200，具体用于获取堆芯外部探测器的响应，以及堆芯在标定过程中的基准功率分布参数；根据堆芯在标定过程中的基准功率分布参数和堆芯外部探测器的响应，确定堆芯在标定过程中的探测电流；根据堆芯在标定过程中的探测电流和堆芯在标定过程中的基准物理参数，确定标定参数。
上述核仪表系统指示信号的量化装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中，提供了一种计算机设备，该计算机设备可以是服务器，其内部结构图可以如图1所示。该计算机设备包括处理器、存储器、输入/输出接口(Input/Output，简称I/O)和通信接口。其中，处理器、存储器和输入/输出接口通过系统总线连接，通信接口通过输入/输出接口连接到系统总线。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质和内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的数据库用于存储指示信号数据。该计算机设备的输入/输出接口用于处理器与外部设备之间交换信息。该计算机设备的通信接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种核仪表系统指示信号的量化方法。
本领域技术人员可以理解，图1中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：
通过模拟堆芯在不同运行工况下的功率分布情况，得到堆芯在不同运行工况下的功率分布参数；
根据堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
构建堆芯的实际物理模型；
获取不同运行工况下的特征参数；
将不同运行工况下的特征参数输入至堆芯的实际物理模型，得到堆芯在不同运行工况下的功率分布参数。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
通过模拟堆芯在正常运行工况下的功率分布情况、以及堆芯在事故运行工况下的功率分布情况，得到堆芯在正常运行工况下的功率分布参数、以及堆芯在事故运行工况下的功率分布参数。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
根据堆芯在不同运行工况下的功率分布参数，确定堆芯外部探测器在不同运行工况下的指示信号；
根据外部探测器堆芯在不同运行工况下的指示信号和基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
根据堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数，确定堆芯外部探测器在正常运行工况下和事故工况下的探测电流；
根据堆芯在标定过程中的功率分布参数，确定标定参数；
根据堆芯的标定参数和堆芯外部探测器在正常运行工况下和事故工况下的探测电流，确定堆芯在正常运行工况下和事故工况下的指示信号。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
获取堆芯外部探测器的响应；
根据堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数、以及堆芯外部探测器的响应，确定堆芯外部探测器在正常运行工况下和事故工况下的探测电流。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
获取堆芯外部探测器的响应，以及堆芯在标定过程中的基准功率分布参数；
根据堆芯在标定过程中的基准功率分布参数和堆芯外部探测器的响应，确定堆芯在标定过程中的探测电流；
根据堆芯在标定过程中的探测电流和堆芯在标定过程中的基准物理参数，确定标定参数。
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
通过模拟堆芯在不同运行工况下的功率分布情况，得到堆芯在不同运行工况下的功率分布参数；
根据堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
构建堆芯的实际物理模型；
获取不同运行工况下的特征参数；
将不同运行工况下的特征参数输入至堆芯的实际物理模型，得到堆芯在不同运行工况下的功率分布参数。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
通过模拟堆芯在正常运行工况下的功率分布情况、以及堆芯在事故运行工况下的功率分布情况，得到堆芯在正常运行工况下的功率分布参数、以及堆芯在事故运行工况下的功率分布参数。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
根据堆芯在不同运行工况下的功率分布参数，确定堆芯外部探测器在不同运行工况下的指示信号；
根据堆芯外部探测器在不同运行工况下的指示信号和基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
根据堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数，确定堆芯外部探测器在正常运行工况下和事故工况下的探测电流；
根据堆芯在标定过程中的功率分布参数，确定标定参数；
根据堆芯的标定参数和堆芯外部探测器在正常运行工况下和事故工况下的探测电流，确定堆芯外部探测器在正常运行工况下和事故工况下的指示信号。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
获取堆芯外部探测器的响应；
根据堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数、以及堆芯外部探测器的响应，确定堆芯外部探测器在正常运行工况下和事故工况下的探测电流。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
获取堆芯外部探测器的响应，以及堆芯在标定过程中的基准功率分布参数；
根据堆芯在标定过程中的基准功率分布参数和堆芯外部探测器的响应，确定堆芯在标定过程中的探测电流；
根据堆芯在标定过程中的探测电流和堆芯在标定过程中的基准物理参数，确定标定参数。
在一个实施例中，提供了一种计算机程序产品，包括计算机程序，该计算机程序被处理器执行时实现以下步骤：
通过模拟堆芯在不同运行工况下的功率分布情况，得到堆芯在不同运行工况下的功率分布参数；
根据堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
构建堆芯的实际物理模型；
获取不同运行工况下的特征参数；
将不同运行工况下的特征参数输入至堆芯的实际物理模型，得到堆芯在不同运行工况下的功率分布参数。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
通过模拟堆芯在正常运行工况下的功率分布情况、以及堆芯在事故运行工况下的功率分布情况，得到堆芯在正常运行工况下的功率分布参数、以及堆芯在事故运行工况下的功率分布参数。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
根据堆芯在不同运行工况下的功率分布参数，确定堆芯外部探测器在不同运行工况下的指示信号；
根据堆芯外部探测器在不同运行工况下的指示信号和基准功率分布参数，基于非参数统计分析方法对堆芯外部探测器的指示信号进行量化。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
根据堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数，确定堆芯外部探测器在正常运行工况下和事故工况下的探测电流；
根据堆芯在标定过程中的功率分布参数，确定标定参数；
根据堆芯的标定参数和堆芯外部探测器在正常运行工况下和事故工况下的探测电流，确定堆芯外部探测器在正常运行工况下和事故工况下的指示信号。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
获取堆芯外部探测器的响应；
根据堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数、以及堆芯外部探测器的响应，确定堆芯外部探测器在正常运行工况下和事故工况下的探测电流。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
获取堆芯外部探测器的响应，以及堆芯在标定过程中的基准功率分布参数；
根据堆芯在标定过程中的基准功率分布参数和堆芯外部探测器的响应，确定堆芯在标定过程中的探测电流；
根据堆芯在标定过程中的探测电流和堆芯在标定过程中的基准物理参数，确定标定参数。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、数据库或其它介质的任何引用，均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(Read-OnlyMemory，ROM)、磁带、软盘、闪存、光存储器、高密度嵌入式非易失性存储器、阻变存储器(ReRAM)、磁变存储器(Magnetoresistive RandomAccess Memory，MRAM)、铁电存储器(Ferroelectric RandomAccess Memory，FRAM)、相变存储器(Phase Change Memory，PCM)、石墨烯存储器等。易失性存储器可包括随机存取存储器(RandomAccess Memory，RAM)或外部高速缓冲存储器等。作为说明而非局限，RAM可以是多种形式，比如静态随机存取存储器(Static RandomAccess Memory，SRAM)或动态随机存取存储器(Dynamic RandomAccessMemory，DRAM)等。本申请所提供的各实施例中所涉及的数据库可包括关系型数据库和非关系型数据库中至少一种。非关系型数据库可包括基于区块链的分布式数据库等，不限于此。本申请所提供的各实施例中所涉及的处理器可为通用处理器、中央处理器、图形处理器、数字信号处理器、可编程逻辑器、基于量子计算的数据处理逻辑器等，不限于此。
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本申请专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请的保护范围应以所附权利要求为准。

Claims (10)

1.一种核仪表系统指示信号的量化方法，其特征在于，所述方法包括：
通过模拟堆芯在不同运行工况下的功率分布情况，得到所述堆芯在不同运行工况下的功率分布参数；
根据所述堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对所述堆芯外部探测器的指示信号进行量化。
2.根据权利要求1所述的方法，其特征在于，所述通过模拟堆芯在不同运行工况下的功率分布情况，得到所述堆芯在不同运行工况下的功率分布参数，包括：
构建所述堆芯的实际物理模型；
获取所述不同运行工况下的特征参数；
将所述不同运行工况下的特征参数输入至所述堆芯的实际物理模型，得到所述堆芯在不同运行工况下的功率分布参数。
3.根据权利要求1或2所述的方法，其特征在于，所述不同运行工况包括正常运行工况和事故运行工况，所述通过模拟堆芯在不同运行工况下的功率分布情况，得到所述堆芯在不同运行工况下的功率分布参数，包括：
通过模拟所述堆芯在正常运行工况下的功率分布情况、以及所述堆芯在事故运行工况下的功率分布情况，得到所述堆芯在正常运行工况下的功率分布参数、以及所述堆芯在事故运行工况下的功率分布参数。
4.根据权利要求1所述的方法，其特征在于，所述根据所述堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对所述堆芯外部探测器的指示信号进行量化，包括：
根据所述堆芯在不同运行工况下的功率分布参数，确定所述堆芯外部探测器在不同运行工况下的指示信号；
根据所述堆芯外部探测器在不同运行工况下的指示信号和所述基准功率分布参数，基于所述非参数统计分析方法对所述堆芯外部探测器的指示信号进行量化。
5.根据权利要求4所述的方法，其特征在于，所述不同运行工况包括正常运行工况和事故运行工况，所述根据所述堆芯在不同运行工况下的功率分布参数，确定所述堆芯外部探测器在不同运行工况下的指示信号，包括：
根据所述堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数，确定所述堆芯外部探测器在正常运行工况下和事故工况下的探测电流；
根据所述堆芯在标定过程中的功率分布参数，确定标定参数；
根据所述标定参数和所述堆芯外部探测器在正常运行工况下和事故工况下的探测电流，确定所述堆芯外部探测器在正常运行工况下和事故工况下的指示信号。
6.根据权利要求5所述的方法，其特征在于，所述根据所述堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数，确定所述堆芯外部探测器在正常运行工况下和事故工况下的探测电流，包括：
获取所述堆芯外部探测器的响应；
根据所述堆芯在正常运行工况下的功率分布参数和在事故运行工况下的功率分布参数、以及所述堆芯外部探测器的响应，确定所述堆芯外部探测器在正常运行工况下和事故工况下的探测电流。
7.根据权利要求5所述的方法，其特征在于，所述根据所述堆芯在标定过程中的功率分布参数，确定标定参数，包括：
获取所述堆芯外部探测器的响应，以及所述堆芯在标定过程中的基准功率分布参数；
根据所述堆芯在标定过程中的基准功率分布参数和所述堆芯外部探测器的响应，确定所述堆芯在标定过程中的探测电流；
根据所述堆芯在标定过程中的探测电流和所述堆芯在标定过程中的基准物理参数，确定所述标定参数。
8.一种核仪表系统指示信号的量化装置，其特征在于，所述包括：
获取模块，用于通过模拟堆芯在不同运行工况下的功率分布情况，得到所述堆芯在不同运行工况下的功率分布参数；
量化模块，用于根据所述堆芯在不同运行工况下的功率分布参数和对应的基准功率分布参数，基于非参数统计分析方法对所述堆芯外部探测器的指示信号进行量化。
9.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至7中任一项所述的方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至7中任一项所述的方法的步骤。
CN202310283228.5A 2023-03-22 2023-03-22 核仪表系统指示信号的量化方法、装置、设备和存储介质 Pending CN116467849A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310283228.5A CN116467849A (zh) 2023-03-22 2023-03-22 核仪表系统指示信号的量化方法、装置、设备和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310283228.5A CN116467849A (zh) 2023-03-22 2023-03-22 核仪表系统指示信号的量化方法、装置、设备和存储介质

Publications (1)

Publication Number Publication Date
CN116467849A true CN116467849A (zh) 2023-07-21

Family

ID=87181510

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310283228.5A Pending CN116467849A (zh) 2023-03-22 2023-03-22 核仪表系统指示信号的量化方法、装置、设备和存储介质

Country Status (1)

Country Link
CN (1) CN116467849A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116879513A (zh) * 2023-09-07 2023-10-13 中碳实测(北京)科技有限公司 气体分析系统的验证方法、装置、设备及存储介质

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116879513A (zh) * 2023-09-07 2023-10-13 中碳实测(北京)科技有限公司 气体分析系统的验证方法、装置、设备及存储介质
CN116879513B (zh) * 2023-09-07 2023-11-14 中碳实测(北京)科技有限公司 气体分析系统的验证方法、装置、设备及存储介质

Similar Documents

Publication Publication Date Title
Cacuci The second-order adjoint sensitivity analysis methodology
CN116467849A (zh) 核仪表系统指示信号的量化方法、装置、设备和存储介质
CN115147012B (zh) 一种基于神经网络模型的碳排放量核算方法及装置
Loy et al. Diagnostic tools for hierarchical linear models
CN112529296A (zh) 水质预测方法、装置及服务器
CN115841046A (zh) 基于维纳过程的加速退化试验数据处理方法和装置
CN108877969B (zh) 核功率理论模型建立及验证方法、系统及终端设备
CN116227045B (zh) 一种结构试件的局部应力应变场构造方法及系统
CN116484209A (zh) 大坝结构监测异常识别方法及装置
CN114487976B (zh) Mcm电子式互感器校验仪溯源不确定度评定方法及系统
CN114490412A (zh) 基于自减逆向云发生器的三维cad软件性能度量方法及装置
CN114625731A (zh) 核数据截面库生成方法、装置、设备及存储介质
CN115392715A (zh) 一种用电数据风险评估方法、装置、设备及存储介质
Radaideh et al. Sensitivity and uncertainty analysis of the fundamental delayed neutron data in lwrs
CN111027680B (zh) 基于变分自编码器的监控量不确定性预测方法及系统
Faure et al. Fuel performance uncertainty quantification and sensitivity analysis in the presence of epistemic and aleatoric sources of uncertainties
CN114707310A (zh) 基于奇异值分解变换的核反应堆物理仿真方法、装置、计算机设备及存储介质
CN113986970A (zh) 一种基于基线库数据的量费计算结果检测方法
CN112488528A (zh) 数据集的处理方法、装置、设备和存储介质
Burr et al. Assessing confidence in inferring reactor type and fuel burnup: A Markov chain Monte Carlo approach
Pizarro-Vallejos Use Of The Fission Matrix For Steady State And Quasi-static Kinetic Modeling Of The Treat Reactor With Temperature Feedback
CN112379318B (zh) 一种功率量程探测器的系数校准方法、装置、设备和介质
Chin et al. Improvement on Fourier fundamental mode coefficient for Monte Carlo criticality calculations
Druzhaev et al. Method of WWER active core power parameters evaluation by readings of ex-core neutron flux monitoring system using additional information
CN117153438A (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