CN115048809A - 一种用于多物理多尺度耦合系统的不确定性量化方法 - Google Patents

Abstract

本发明涉及反应堆热工水力技术领域，具体涉及一种用于多物理多尺度耦合系统的不确定性量化方法，包括以下步骤：收集用于多物理多尺度耦合系统不确定性评价和验证的实验数据，并建立实验数据库；针对目标分析对象装置，结合多物理多尺度耦合系统不确定性源识别装置在目标工况中存在的不确定性源及其程序表示，并将不确定性根据来源进行分类；对不同种类的不确定性输入源进行评估和量化；对目标分析装置和工况，识别和量化多物理多尺度耦合系统中各输入不确定性源后，使用量化得到的各输入参数的不确定性分布执行定量敏感性分析；执行不确定性传播计算。本发明能够识别和量化多物理多尺度耦合系统中存在的不确定性。

Description

一种用于多物理多尺度耦合系统的不确定性量化方法

技术领域

本发明涉及反应堆热工水力技术领域，具体涉及一种用于多物理多尺度耦合系统的不确定性量化方法。

背景技术

核反应堆的安全分析依赖于计算机程序，由于核反应堆系统十分复杂，不同分析工况下的现象进程以及分析对象均存在差异，因此存在诸多不同种类的安全分析软件，如中子物理程序、燃料行为程序、热工水力系统程序、安全壳热工水力程序等。这些分析程序的功能各异，各程序具有特定的使用范围，且适用的物理场和尺度均存在一定区别。如中子物理类程序主要计算反应堆中子动力学，其分析尺度涵盖微观中子截面至宏观堆芯；燃料行为类程序主要计算燃料芯块及包壳的温度场及应力，其涵盖介观材料蠕变和宏观包壳变形、爆裂；而热工水力系统程序主要计算反应堆整体冷却剂系统的瞬态宏观表现。

在对核反应堆进行事故工况安全分析时，需使用所有种类的安全分析程序进行综合分析。传统方法中采用“数据接口”的方式耦合各类程序，即把一类程序计算结束后得到的输出作为其他程序的输入，而该过程不是实时的，故这种程序耦合方法也被称为弱耦合。弱耦合无法真实模拟核反应堆的真实状态，存在过多的冗余保守性，因此有必要发展强耦合技术。相比之下，强耦合能够构建一个实时分析的耦合系统，将系统内所有程序的数据串通，实现不同功能程序间的数据实时交互，大大提高模拟结果的精度。

其中，基于强耦合技术构建的多物理多尺度核反应堆安全分析系统平台也被称为数字堆，具有最佳估算的能力。最佳估算计算中要求考虑分析中的不确定性，因此，有必要开发一种方法用于量化多物理多尺度耦合系统中存在的各类不确定性。

发明内容

本发明提供了一种用于多物理多尺度耦合系统的不确定性量化方法，能够识别和量化多物理多尺度耦合系统中存在的各类不确定性、量化耦合系统中各功能程序的计算不确定性，并能通过不确定性传播计算量化整个耦合系统的计算不确定性，实现最佳估算加不确定性分析的目标。

本发明通过下述技术方案实现：

一种用于多物理多尺度耦合系统的不确定性量化方法，包括以下步骤：

S10、收集用于多物理多尺度耦合系统不确定性评价和验证的实验数据，并建立实验数据库；

S20、针对目标分析对象装置，结合多物理多尺度耦合系统不确定性源识别装置在目标工况中存在的不确定性源及其程序表示，并将不确定性根据来源进行分类；

S30、对不同种类的不确定性输入源进行评估和量化；

S40、对目标分析装置和工况，识别和量化多物理多尺度耦合系统中各输入不确定性源后，使用量化得到的各输入参数的不确定性分布执行定量敏感性分析；

S50、执行不确定性传播计算，以量化多物理多尺度耦合系统的计算不确定性。

本发明提供的用于多物理多尺度耦合系统的不确定性量化方法，先构建用于评估和量化耦合系统不确定性所需的基准题和数据库，并识别多物理多尺度耦合系统中存在的不确定性源和对其进行分类，后对不同种类的不确定性源进行量化和评估，最后量化多物理多尺度耦合系统的计算不确定性，能够识别和量化多物理多尺度耦合系统中存在的各类不确定性、量化耦合系统中各功能程序的计算不确定性，并能通过不确定性传播计算量化整个耦合系统的计算不确定性，实现最佳估算加不确定性分析的目标。

具体而言，步骤S10中建立实验数据库包括以下步骤：

S11、将用于多物理多尺度耦合系统不确定性评价和验证的实验数据，按照不同物理场进行整理，并初步分类为分离效应实验和整体效应实验；

S12、基于多准则决策分析和试验源分配技术开发数据库分类方法，将实验数据按照工况、组件、现象以及参数范围进行分类，形成规范化实验数据库；

S13、基于得到的完善数据库建立评估和验证用基准题，建立标准化数据库使用规范。

具体而言，步骤S20中不确定性分类方法包括以下步骤：

S21、针对目标分析对象装置，使用层次分析法对其进行尺度由大至小的结构化分析；

S22、针对局部部件的不同物理场使用层次分析法进行分析；

S23、结合多物理多尺度耦合系统，识别目标装置在目标工况中存在的不确定性源及其程序表示，并将不确定性根据来源进行分类。

具体而言，步骤S22中的物理场包括中子物理场、温度场、流场、应力场、水化学场。

具体而言，步骤S23中不确定性的分类包括核数据库、几何结构、设计测量、材料物性、初始/边界条件、本构模型以及数值离散。

具体而言，步骤S30中的评估和量化的方法为：

核数据库不确定性，使用微扰分析法量化其不确定性；

几何结构、设计测量以及材料物性的不确定性，基于设计、施工或实验过程的测量不确定度给定；

初始/边界条件，基于建立的数据库，结合参数统计法评估其不确定性；

程序中本构模型的不确定性，使用参数/非参数统计法并结合开发的实验数据库对齐进行量化；

程序数值离散的不确定性，使用证据理论结方法量化其不确定性。

具体而言，所述参数统计法包括假设检验和频数分析法，所述非参数统计法包括非参数曲线估计法、基于贝叶斯理论的马尔科夫链蒙特卡罗以及基于快速傅里叶变换的方法。

具体而言，步骤S40中，对于目标分析装置和工况，识别和量化多物理多尺度耦合系统中各输入不确定性源后，使用量化得到的各输入参数的不确定性分布执行采用矩独立全局敏感性进行定量敏感性分析计算，并基于敏感性分析结果筛选对目标输出次要或无影响的输入参数和确定影响参数。

具体而言，步骤S50中采用逐层耦合的方式传播和量化。

具体而言，步骤S50中传播和量化的方法包括以下步骤：

S51、针对微观尺度程序，使用步骤S30中量化得到的重要输入参数不确定性分布，结合步骤S10中开发得到的基准题执行不确定性随机抽样计算，得到微观尺度程序对应输出的不确定性传播抽样计算结果，通过比较不确定性抽样计算结果与基准题中对应的实验测量数据，量化微观尺度程序的计算不确定性；

S52、将所有微观程序视为独立的黑匣子，并将微观程序的输出不确定性分布作为介观尺度程序的输入不确定性，执行不确定性随机抽样传播计算以量化介观程序的计算结果不确定性分布；

S53、以步骤S52相同的原理将介观分析程序视为黑匣子作为宏观分析程序的不确定性本构模块，量化宏观分析程序的输出不确定性分布。

本发明具有如下的优点和有益效果：

1、本发明提供的用于多物理多尺度耦合系统的不确定性量化方法，先构建用于评估和量化耦合系统不确定性所需的基准题和数据库，并识别多物理多尺度耦合系统中存在的不确定性源和对其进行分类，后对不同种类的不确定性源进行量化和评估，最后量化多物理多尺度耦合系统的计算不确定性，在使用多物理多尺度耦合系统对目标核反应堆进行工况分析时，能够识别和量化多物理多尺度耦合系统中存在的各类不确定性，能量化耦合系统中各功能程序的计算不确定性，并能通过不确定性传播计算量化整个耦合系统的计算不确定性，实现最佳估算加不确定性分析的目标。

2、本发明提供的用于多物理多尺度耦合系统的不确定性量化方法，能够解决传统耦合系统中存在的冗余保守性的问题，由于多物理多尺度耦合系统具有最佳估算能力，在考虑了不确定性的情况下其能够真实的模拟事故进程，大大减小了程序模拟现象失真的可能性，实现核反应堆瞬态工况“现实性”模拟的安全分析要求，增加了安全分析结果的可靠性，并能够提供更多的安全裕量，因此，本发明对于核反应堆安全分析方法的完善以及高精度数字堆系统的构建具有重要促进作用。

附图说明

此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部分，并不构成对本发明实施例的限定。

在附图中：

图1为本发明实施例用于多物理多尺度耦合系统的不确定性量化方法的流程示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。

实施例1

结合图1，本实施例提供了一种用于多物理多尺度耦合系统的不确定性量化方法，包括以下步骤：

S10、收集用于多物理多尺度耦合系统不确定性评价和验证的实验数据，并建立实验数据库。

具体而的，建立实验数据库包括以下步骤：

S11、收集用于多物理多尺度耦合系统不确定性评价和验证的实验数据后缘，将其按照不同物理场进行整理，并初步分类为分离效应实验和整体效应实验；

S12、基于多准则决策分析和试验源分配技术开发数据库分类方法，将实验数据按照工况、组件、现象以及参数范围进行分类，形成规范化实验数据库；

S13、基于得到的完善数据库建立评估和验证用基准题，建立标准化数据库使用规范，以减小不确定性量化过程中人为因素引入的影响。

S20、针对目标分析对象装置，结合多物理多尺度耦合系统不确定性源识别装置在目标工况中存在的不确定性源及其程序表示，并将不确定性根据来源进行分类。

具体来说，不确定性分类方法包括以下步骤：

S21、针对目标分析对象装置，使用层次分析法对其进行尺度由大至小的结构化分析，顺序为整体核反应堆装置-系统-组件-部件；

S22、针对局部部件的不同物理场使用层次分析法进行分析，其中物理场包括中子物理场、温度场、流场、应力场、水化学场；

S23、结合多物理多尺度耦合系统，识别目标装置在目标工况中存在的不确定性源及其程序表示，并将不确定性根据来源进行分类，其中不确定性的分类包括核数据库、几何结构、设计测量、材料物性、初始/边界条件、本构模型以及数值离散。

S30、对不同种类的不确定性输入源进行评估和量化。

具体来讲，评估和量化的方法为：

针对核数据库不确定性，使用微扰分析法量化其不确定性；

对于几何结构、设计测量以及材料物性的不确定性，基于设计、施工或实验过程的测量不确定度给定；

对于初始/边界条件，基于建立的数据库，结合参数统计法评估其不确定性，最常使用的参数统计法包括假设检验和频数分析法；

对于程序中本构模型的不确定性，使用参数/非参数统计法并结合开发的实验数据库对齐进行量化，最常使用的参数统计法包括假设检验和频数分析法，最常使用的非参数统计法包括非参数曲线估计法、基于贝叶斯理论的马尔科夫链蒙特卡罗以及基于快速傅里叶变换的方法；

对于程序数值离散的不确定性，使用证据理论结或同时结合经验的方法量化其不确定性。

S40、对目标分析装置和工况，识别和量化多物理多尺度耦合系统中各输入不确定性源后，使用量化得到的各输入参数的不确定性分布执行定量敏感性分析。

具体的，对于目标分析装置和工况，识别和量化多物理多尺度耦合系统中各输入不确定性源后，使用量化得到的各输入参数的不确定性分布执行采用矩独立全局敏感性进行定量敏感性分析计算，并基于敏感性分析结果筛选对目标输出次要或无影响的输入参数和确定影响参数(影响较大的重要参数)。

S50、执行不确定性传播计算，以量化多物理多尺度耦合系统的计算不确定性。

即，确定目标分析工况的重要输入参数及其不确定性分布后，执行不确定性传播计算以量化多物理多尺度耦合系统的计算不确定性。其中，考虑到多物理多尺度耦合系统的复杂性，因此使用逐层耦合的方式传播和量化。

具体而言，传播和量化的方法包括以下步骤：

S51、针对微观尺度程序，如中子截面类程序，使用步骤S30中量化得到的重要输入参数不确定性分布，结合步骤S10中开发得到的基准题执行不确定性随机抽样计算，得到微观尺度程序对应输出的不确定性传播抽样计算结果，通过比较不确定性抽样计算结果与基准题中对应的实验测量数据，量化微观尺度程序的计算不确定性；

S52、将所有微观程序视为独立的黑匣子，其作为介观尺度程序的一个本构模块，在使用介观尺度程序，如燃料行为分析程序，进行计算时不再考虑微观程序的输入不确定性，直接将微观程序的输出不确定性分布作为介观尺度程序的输入不确定性，进而执行不确定性随机抽样传播计算以量化介观程序的计算结果不确定性分布；

S53、以步骤S52相同的原理将介观分析程序视为黑匣子作为宏观分析程序的不确定性本构模块，量化宏观分析程序的输出不确定性分布，以逐层耦合方式按照微观-介观-宏观的顺序传播和量化耦合系统中的不确定性。

因此，步骤S50通过逐层耦合评价不确定性的方法，不仅实现了量化多物理多尺度耦合系统中不同尺度单个功能程序的计算不确定性，并且实现了不确定性在耦合系统中的传播，能得到最终分析所需的目标输出的不确定性分布。

综上，本实施例提供的方法，并识别多物理多尺度耦合系统中存在的不确定性源和对其进行分类，后对不同种类的不确定性源进行量化和评估，最后量化多物理多尺度耦合系统的计算不确定性，在使用多物理多尺度耦合系统对目标核反应堆进行工况分析时，能够识别和量化多物理多尺度耦合系统中存在的各类不确定性、量化耦合系统中各功能程序的计算不确定性，并能通过不确定性传播计算量化整个耦合系统的计算不确定性，实现最佳估算加不确定性分析的目标。

以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种用于多物理多尺度耦合系统的不确定性量化方法，其特征在于，包括以下步骤：

S10、收集用于多物理多尺度耦合系统不确定性评价和验证的实验数据，并建立实验数据库；

S20、针对目标分析对象装置，结合多物理多尺度耦合系统不确定性源识别装置在目标工况中存在的不确定性源及其程序表示，并将不确定性根据来源进行分类；

S30、对不同种类的不确定性输入源进行评估和量化；

S40、对目标分析装置和工况，识别和量化多物理多尺度耦合系统中各输入不确定性源后，使用量化得到的各输入参数的不确定性分布执行定量敏感性分析；

S50、执行不确定性传播计算，以量化多物理多尺度耦合系统的计算不确定性。

2.根据权利要求1所述的一种用于多物理多尺度耦合系统的不确定性量化方法，其特征在于，步骤S10中建立实验数据库包括以下步骤：

S11、将用于多物理多尺度耦合系统不确定性评价和验证的实验数据，按照不同物理场进行整理，并初步分类为分离效应实验和整体效应实验；

S12、基于多准则决策分析和试验源分配技术开发数据库分类方法，将实验数据按照工况、组件、现象以及参数范围进行分类，形成规范化实验数据库；

S13、基于得到的完善数据库建立评估和验证用基准题，建立标准化数据库使用规范。

3.根据权利要求1所述的一种用于多物理多尺度耦合系统的不确定性量化方法，其特征在于，步骤S20中不确定性分类方法包括以下步骤：

S21、针对目标分析对象装置，使用层次分析法对其进行尺度由大至小的结构化分析；

S22、针对局部部件的不同物理场使用层次分析法进行分析；

S23、结合多物理多尺度耦合系统，识别目标装置在目标工况中存在的不确定性源及其程序表示，并将不确定性根据来源进行分类。

4.根据权利要求3所述的一种用于多物理多尺度耦合系统的不确定性量化方法，其特征在于，步骤S22中的物理场包括中子物理场、温度场、流场、应力场、水化学场。

5.根据权利要求3所述的一种用于多物理多尺度耦合系统的不确定性量化方法，其特征在于，步骤S23中不确定性的分类包括核数据库、几何结构、设计测量、材料物性、初始/边界条件、本构模型以及数值离散。

6.根据权利要求5所述的一种用于多物理多尺度耦合系统的不确定性量化方法，其特征在于，步骤S30中的评估和量化的方法为：

核数据库不确定性，使用微扰分析法量化其不确定性；

几何结构、设计测量以及材料物性的不确定性，基于设计、施工或实验过程的测量不确定度给定；

初始/边界条件，基于建立的数据库，结合参数统计法评估其不确定性；

程序中本构模型的不确定性，使用参数/非参数统计法并结合开发的实验数据库对齐进行量化；

程序数值离散的不确定性，使用证据理论结方法量化其不确定性。

7.根据权利要求6所述的一种用于多物理多尺度耦合系统的不确定性量化方法，其特征在于，所述参数统计法包括假设检验和频数分析法，所述非参数统计法包括非参数曲线估计法、基于贝叶斯理论的马尔科夫链蒙特卡罗以及基于快速傅里叶变换的方法。

8.根据权利要求1所述的一种用于多物理多尺度耦合系统的不确定性量化方法，其特征在于，步骤S40中，对于目标分析装置和工况，识别和量化多物理多尺度耦合系统中各输入不确定性源后，使用量化得到的各输入参数的不确定性分布执行采用矩独立全局敏感性进行定量敏感性分析计算，并基于敏感性分析结果筛选对目标输出次要或无影响的输入参数和确定影响参数。

9.根据权利要求7所述的一种用于多物理多尺度耦合系统的不确定性量化方法，其特征在于，步骤S50中采用逐层耦合的方式传播和量化。

10.根据权利要求9所述的一种用于多物理多尺度耦合系统的不确定性量化方法，其特征在于，步骤S50中传播和量化的方法包括以下步骤：

S51、针对微观尺度程序，使用步骤S30中量化得到的重要输入参数不确定性分布，结合步骤S10中开发得到的基准题执行不确定性随机抽样计算，得到微观尺度程序对应输出的不确定性传播抽样计算结果，通过比较不确定性抽样计算结果与基准题中对应的实验测量数据，量化微观尺度程序的计算不确定性；

S52、将所有微观程序视为独立的黑匣子，并将微观程序的输出不确定性分布作为介观尺度程序的输入不确定性，执行不确定性随机抽样传播计算以量化介观程序的计算结果不确定性分布；

S53、以步骤S52相同的原理将介观分析程序视为黑匣子作为宏观分析程序的不确定性本构模块，量化宏观分析程序的输出不确定性分布。

Images