一种基于敏感性和不确定性分析的目标精度评估方法技术

一种基于敏感性和不确定性分析的目标精度评估方法，1、针对传统反应堆物理计算两步法，基于微扰理论获取组件计算和堆芯计算目标参数的灵敏度系数；2、采用评价库加工程序NJOY，制作多群截面协方差数据库，结合目标参数灵敏度系数，获取目标参数的不确定度；3、对比目标参数的不确定度限制，对于不确定度大于限制的目标参数，建立带约束条件的最小值问题，采用序列二次规划算法等方法获取指定核数据的不确定度限制，为核数据的优化提供数据支持；能够根据不确定度超过限制的目标参数，重新评估核数据精度，同时该评估具有针对性，对特定核数据的改进具有指导意义。

【技术实现步骤摘要】

本专利技术涉及核反应堆堆芯设计和安全，以及核数据评估领域，具体涉及一种根据核反应堆设计目标精度，对核数据进行精度评估的方法。
技术介绍
随着核反应堆战略地位的不断提高，其安全性受到进一步的重视，尤其是福岛事故之后，国际上对核反应堆的安全性更是提出了前所未有的苛刻要求，这使得人们对核反应堆设计计算精度的要求也不断提高。中子学计算是核反应堆设计中最基础和关键的问题之一，其计算精度决定着核反应堆设计的安全性，长期以来，如何提高中子学计算的精度一直是国际上研究的热点和重点，近年来各种先进的中子学计算方法得到长足的发展，由物理模型、数值方法造成的误差大为降低。而核数据作为中子学计算最基本的输入参数之一，由于其测量值不可避免地存在不确定度，已逐渐成为核反应堆中子学计算不确定度的重要源头。目前，反应堆中子学计算由核数据引入的不确定度已经得到较为广泛的研究，不少数值分析结果表明目前的核数据精度无法满足反应堆设计的精度要求。所以提高核数据的测量精度成为减小计算结果不确定度，满足核反应设计要求的重要手段。然而，海量的核数据都带有自身的不确定度，不可能要求每一种核数据都经过重新测量，这样既没有针对性又浪费大量的人力财力。目前并没有一种有效的方法来对海量的核数据进行评估，因此有必要专利技术一种方法，可以对特定核数据进行敏感性和不确定性分析，有针对性地对特定核数据进行评估，为其改进提供数据支持。
技术实现思路
为了根据核反应堆目标设计参数的不确定度要求，对核数据的精度改进提供数据支持，本专利技术提出一种基于敏感性和不确定性分析的目标精度评估方法，能够有针对性地对某些核数据的改进提供数据支持。为了实现上述目的，本专利技术采用以下技术方案予以实施：一种基于敏感性和不确定性分析的目标精度评估方法，该方法包括以下步骤：步骤1：根据传统的反应堆物理计算两步法，基于微扰理论，分别获取核反应堆组件计算和堆芯计算的计算结果的灵敏度系数；其中组件计算结果包括组件目标参数即有效增殖因子、功率，以及作为堆芯计算的输入的组件参数；堆芯计算结果包括堆芯目标参数即有效增殖因子、归一化功率；为了叙述简便，如无特别说明，以下将组件或堆芯的计算结果统一称为响应；传统的核反应堆物理计算两步法中，组件计算进行中子输运计算，堆芯计算进行中子扩散计算；设反应堆物理组件计算或者堆芯计算某响应为R，表示成含通量或者共轭通量的泛函，设为：式中：H1，H2——依赖于截面数据的算子；Φ——中子角通量密度或标通量密度；Φ*——共轭中子角通量密度或标通量密度；<>——表示在求解空间积分；那么，该响应R对核数据α的相对灵敏度系数的计算式为：其中Γ*和Γ分别称为广义共轭通量和广义通量；Γ*为广义共轭方程的解，广义共轭方程为：式中：M*——输运或扩散算子的共轭算子；Γ为广义方程的解，广义方程为：M——输运算子或扩散算子；采用成熟的组件模块化特征线方法作为二维输运求解方法，对中子输运方程、共轭中子输运方程、广义方程以及广义共轭输运方程进行求解，获取这些方程的通量解，再根据式(5)求得响应R的相对灵敏度系数；上述计算过程获取的是组件计算或者堆芯计算响应对核数据扰动的直接影响，称为显式敏感性；但是，在组件计算的中子输运计算之前，要通过共振计算获取共振核素的有效自屏截面，因此获取组件响应的灵敏度系数的时候，为了使得灵敏度系数的计算更加精确，需要考虑核数据α对共振计算影响，此时灵敏度系数为式中：j——共振核素标识；x——共振反应标识；g——能群标识；——共振核素j的x反应的第g群共振截面；式(8)中右端第一项即为上述计算过程获取的显式灵敏度系数；右端第二项表示核数据α通过影响共振计算对组件计算响应的间接影响，称为隐式敏感性；针对共振计算方法——子群方法，进行式(8)中右端第二项，即隐式灵敏度系数的求解；首先对有效自屏截面σx,g建立广义子群共轭方程：式中：——子群输运算子的共轭算子；——子群广义共轭通量；——子群广义共轭源；其中式中i——子群标识；I——子群总数——第i个子群的源项；σx,g,i——第i个子群截面；其次，采用成熟的模块化特征线方法解得子群广义共轭通量，根据公式(16)，得到共振自屏截面σx,g对核数据α的灵敏度系数式中：L——子群输运算子；——子群截面对核数据α的灵敏度系数；——子群输运方程源项对核数据α的灵敏度系数；采用基于窄共振近似的子群方法中，当α为非共振核素核数据时，源项Qg,i与核数据α无关，因此当核数据α为相关共振核素核数据时，为子群概率pg,i对核数据α的灵敏度系数；此外，式(16)中表示子群截面对核数据α的灵敏度系数子群截面σx,g,i和子群概率pg,i统称为子群参数，统一用表示，子群参数的对核数据α的灵敏度系数采用直接扰动方法求解；对于各共振核素，逐群扰动其连续能量截面，即扰动该共振核素的共振积分表，设扰动百分比为δ，则根据差商代替微分的方法求得子群参数的灵敏度系数，即式中：为未扰动的子群参数，即子群截面σx,g,i或子群概率pg,i；和分别为正向和负向扰动核数据α时的子群参数；δ为α的扰动百分比；至此，便根据式(16)获取共振自屏截面对核数据α的灵敏度系数；结合前面获取的组件响应的显式灵敏度系数，按式(18)计算隐式灵敏度系数：当通过以上步骤计算出组件响应对多群核数据的灵敏度系数，以及堆芯响应对堆芯计算的输入参数，即组件参数的灵敏度系数之后，根据灵敏度系数的传递关系，堆芯响应对多群核数据的灵敏度系数为：式中：a——组件参数标志；c——堆芯响应标志；h——堆芯响应标识；m——组件参数标识；M——组件参数总数；——堆芯响应对组件参数灵敏度系数；——组件参数对多群核数据灵敏度系数；——堆芯响应对多群核数据的灵敏度系数；步骤2：采用核数据库加工程序NJOY，制作核数据的协方差数据库；结合步骤1获取的灵敏度系数，得到组件响应和堆芯响应的不确定度；协方差数据描述了核数据本身的不确定度，保存于核评价数据库中；采用NJOY程序，从评价核数据库出发加工出与组件计算能群结构相一致的多群协方差数据；基于加工出的多群协方差数据和步骤1得到的灵敏度系数，计算出组件响应的协方差数据：式中：m，n——组件响应标识；p，q——核数据标识；——第m个组件响应；——第n个组件响应；——第m个组件响应对第p种核数据的灵敏度系数；——第n个组件响应对第q种核数据的灵敏度系数；——第p种和第q种核数据之间的相对协方差；——第m个组件响应与第n个组件响应之间的相对协方差；组件参数包含于上述组件响应的范畴，是堆芯计算的输入参数，在获取组件参数的不确定度之后，结合步骤1得到的堆芯响应对组件参数的灵敏度系数，即得到堆芯响应的不确定度；式中：h，k——堆芯响应标识；——第h个堆芯响应；——第k个堆芯响应；——第m个组件参数与n个组件参数之间的相对协方差；——第h个堆芯响应与第k个堆芯响应之间的相对协方差；步骤3：步骤2获取了组件响应和堆芯响应的不确定度，其中包括组件目标参数和堆芯目标参数的不确定度，对比该目标参数的不确定度限制，判断是否小于这个限制；如果不满足要求，就需要对核数据进行评估；设目标参数的精度要求为δ2，建立一个带约束条件的最小值问题：使得：式中：p——核数据标识；y——目标参数标识；P——核本文档来自技高网...

【技术保护点】
一种基于敏感性和不确定性分析的目标精度评估方法，其特征在于：该方法包括以下步骤：步骤1：根据传统的反应堆物理计算两步法，基于微扰理论，分别获取核反应堆组件计算和堆芯计算的计算结果的灵敏度系数；其中组件计算结果包括组件目标参数即有效增殖因子、功率，以及作为堆芯计算的输入的组件参数；堆芯计算结果包括堆芯目标参数即有效增殖因子、归一化功率；为了叙述简便，如无特别说明，以下将组件或堆芯的计算结果统一称为响应；传统的核反应堆物理计算两步法中，组件计算进行中子输运计算，堆芯计算进行中子扩散计算；设反应堆物理组件计算或者堆芯计算某响应为R，表示成含通量或者共轭通量的泛函，设为：式中：H1，H2——依赖于截面数据的算子；Φ——中子角通量密度或标通量密度；Φ*——共轭中子角通量密度或标通量密度；<>——表示在求解空间积分；那么，该响应R对核数据α的相对灵敏度系数的计算式为：其中Γ*和Γ分别称为广义共轭通量和广义通量；Γ*为广义共轭方程的解，广义共轭方程为：式中：M*——输运或扩散算子的共轭算子；Γ为广义方程的解，广义方程为：M——输运算子或扩散算子；采用成熟的组件模块化特征线方法作为二维输运求解方法，对中子输运方程、共轭中子输运方程、广义方程、广义方程以及广义共轭输运方程进行求解，获取这些方程的通量解，再根据式(5)求得响应R的相对灵敏度系数；上述计算过程获取的是组件计算或者堆芯计算响应对核数据扰动的直接影响，称为显式敏感性；但是，在组件计算的中子输运计算之前，要通过共振计算获取共振核素的有效自屏截面，因此获取组件响应的灵敏度系数的时候，为了使得灵敏度系数的计算更加精确，需要考虑核数据α对共振计算影响，此时灵敏度系数为式中：j——共振核素标识；x——共振反应标识；g——能群标识；——共振核素j的x反应的第g群共振截面；式(8)中右端第一项即为上述计算过程获取的显式灵敏度系数；右端第二项表示核数据α通过影响共振计算对组件计算响应的间接影响，称为隐式敏感性；针对共振计算方法——子群方法，进行式(8)中右端第二项，即隐式灵敏度系数的求解；首先对有效自屏截面σx,g建立广义子群共轭方程：式中：——子群输运算子的共轭算子；——子群广义共轭通量；——子群广义共轭源；其中式中i——子群标识；I——子群总数——第i个子群的源项；σx,g,i——第i个子群截面；其次，采用成熟的模块化特征线方法解得子群广义共轭通量，根据公式(16)，得到共振自屏截面σx,g对核数据α的灵敏度系数式中：L——子群输运算子；——子群截面对核数据α的灵敏度系数；——子群输运方程源项对核数据α的灵敏度系数；采用基于窄共振近似的子群方法中，当α为非共振核素核数据时，源项Qg,i与核数据α无关，因此当核数据α为相关共振核素核数据时，为子群概率pg,i对核数据α的灵敏度系数；此外，式(16)中表示子群截面对核数据α的灵敏度系数；子群截面σx,g,i和子群概率pg,i统称为子群参数，统一用表示；子群参数的对核数据α的灵敏度系数采用直接扰动方法求解，对于各共振核素，逐群扰动其连续能量截面，即扰动该共振核素的共振积分表，设扰动百分比为δ，则根据差商代替微分的方法求得子群参数的灵敏度系数，即式中：为未扰动的子群参数，即子群截面σx,g,i或子群概率pg,i；和分别为正向和负向扰动核数据α时的子群参数；δ为α的扰动百分比；至此，便根据式(16)获取共振自屏截面对核数据α的灵敏度系数；结合前面获取的组件响应的显式灵敏度系数，按式(18)计算隐式灵敏度系数：当通过以上步骤计算出组件响应对多群核数据的灵敏度系数，以及堆芯响应对堆芯计算的输入参数，即组件参数的灵敏度系数之后，根据灵敏度系数的传递关系，堆芯响应对多群核数据的灵敏度系数为：式中：a——组件参数标志；c——堆芯响应标志；h——堆芯响应标识；m——组件参数标识；M——组件参数总数；——堆芯响应对组件参数灵敏度系数；——组件参数对多群核数据灵敏度系数；——堆芯响应对多群核数据的灵敏度系数；步骤2：采用核数据库加工程序NJOY，制作核数据的协方差数据库；结合步骤1获取的灵敏度系数，得到组件响应和堆芯响应的不确定度；协方差数据描述了核数据本身的不确定度，保存于核评价数据库中；采用NJOY程序，从评价核数据库出发加工出与组件计算能群结构相一致的多群协方差数据；基于加工出的多群协方差数据和步骤1得到的灵敏度系数，计算出组件响应的协方差数据：式中：m，n——组件响应标识；p，q——核数据标识；——第m个组件响应；——第n个组件响应；——第m个组件响应对第p种核数据的灵敏度系数；——第n个组件响应对第q种核数据的灵敏度系数；——第p种和第q种核数据之间的相对协方差；——第m个组件响应与第n个组件响应之间的相对协方差；组件参数包含于上述组件响应的范畴，是堆...

【技术特征摘要】
1.一种基于敏感性和不确定性分析的目标精度评估方法，其特征在于：该方法包括以下步骤：步骤1：根据传统的反应堆物理计算两步法，基于微扰理论，分别获取核反应堆组件计算和堆芯计算的计算结果的灵敏度系数；其中组件计算结果包括组件目标参数即有效增殖因子、功率，以及作为堆芯计算的输入的组件参数；堆芯计算结果包括堆芯目标参数即有效增殖因子、归一化功率；为了叙述简便，如无特别说明，以下将组件或堆芯的计算结果统一称为响应；传统的核反应堆物理计算两步法中，组件计算进行中子输运计算，堆芯计算进行中子扩散计算；设反应堆物理组件计算或者堆芯计算某响应为R，表示成含通量或者共轭通量的泛函，设为：式中：H1，H2——依赖于截面数据的算子；Φ——中子角通量密度或标通量密度；Φ*——共轭中子角通量密度或标通量密度；<>——表示在求解空间积分；那么，该响应R对核数据α的相对灵敏度系数的计算式为：其中Γ*和Γ分别称为广义共轭通量和广义通量；Γ*为广义共轭方程的解，广义共轭方程为：式中：M*——输运或扩散算子的共轭算子；Γ为广义方程的解，广义方程为：M——输运算子或扩散算子；采用成熟的组件模块化特征线方法作为二维输运求解方法，对中子输运方程、共轭中子输运方程、广义方程、广义方程以及广义共轭输运方程进行求解，获取这些方程的通量解，再根据式(5)求得响应R的相对灵敏度系数；上述计算过程获取的是组件计算或者堆芯计算响应对核数据扰动的直接影响，称为显式敏感性；但是，在组件计算的中子输运计算之前，要通过共振计算获取共振核素的有效自屏截面，因此获取组件响应的灵敏度系数的时候，为了使得灵敏度系数的计算更加精确，需要考虑核数据α对共振计算影响，此时灵敏度系数为式中：j——共振核素标识；x——共振反应标识；g——能群标识；——共振核素j的x反应的第g群共振截面；式(8)中右端第一项即为上述计算过程获取的显式灵敏度系数；右端第二项表示核数据α通过影响共振计算对组件计算响应的间接影响，称为隐式敏感性；针对共振计算方法——子群方法，进行式(8)中右端第二项，即隐式灵敏度系数的求解；首先对有效自屏截面σx,g建立广义子群共轭方程：式中：——子群输运算子的共轭算子；——子群广义共轭通量；——子群广义共轭源；其中式中i——子群标识；I——子群总数——第i个子群的源项；σx,g,i——第i个子群截面；其次，采用成熟的模块化特征线方法解得子群广义共轭通量，根据公式(16)，得到共振自屏截面σx,g对核数据α的灵敏度系数式中：L——子群输运算子；——子群截面对核数据α的灵敏度系数；——子群输运方程源项对核数据α的灵敏度系数；采用基于窄共振近似的子群方法中，当α为非共振核素核数据时，源项Qg,i与核数据α无关，因此当核数据α为相关共振核素核数据时，为子群概率pg,i对核数据α的灵敏度系数；此外，式(16)中表示子群截面对核数据α的灵敏度系数；子群截面σx,g,i和子群概率pg,i统称为子群参数，统一用表示；子群参数的对核数据α的灵敏度系数采用...

【专利技术属性】
技术研发人员：曹良志，刘勇，祖铁军，吴宏春，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：陕西;61