1、核截面引起积分参数keff不确定度的一维分析程序开发报告人： 刚直2014年11月1(1) 选题背景及意义；(2) 国内外研究现状 ；(3) 基本原理和方法；(4) 程序开发；(5) 程序验证及初步应用；(6) 结论；(7) 总结和展望；报告内容：2一. 背景及意义图1-1 堆物理计算系统示意图核数据库（ENDF.）问题模型（输入卡）计算程序（确定论/ M-C ）计算结果计算过程1.反应堆物理计算过程（一）引言32.计算过程误差来源图1-2 堆物理计算误差来源示意图核数据（ENDF.）问题模型（输入卡）计算程序（确定论/ M-C ）计算结果误差来源：核反应截面理论计算参数共振参数角分布能谱数据
2、分类：统计误差系统误差误差来源：几何模型- （1D/2D/3D）几何尺寸材料成分及密度源分布描述 误差来源：理论模型： 输运模型 扩散模型 1维/2维/3维数学方法： M-C统计 Sn离散,散射阶数 收敛准则计算结果= RR（一）引言4 1）随着计算机技术的迅猛发展和计算方法的不断完善，目前发展成熟的先进计算程序本身的系统误差对计算结果的影响逐渐减小，而由核数据引入的计算误差对结果的影响越来越不可忽视。 2） 国外从20世纪70年代开始进行核截面数据对堆参数计算的不确定度的研究，到目前已发展了各类比较成熟的计算程序，并广泛用于核工程设计计算中。（二）背景5 3）在提高反应堆安全和经济性方面有重
3、要作用，如： 提高结果可信度：确定积分参数的误差使结果更合 理可信，有助于进行安全分析。 设计改进：在工程设计中确立主要参数的误差情 况，根据设计限值对其进行调整和优化。 反应堆压力容器延长寿期：通过计算压力容器辐照 剂量等技术参数的不确定度来确定安全裕度，进而预计或延长 其服役期限。 核数据调整：根据对基准实验的计算值和实验值比 较，在不确定度范围内对占主导地位的核截面数据进行适当调 整，使利用调整后的截面数据得到的计算结果更接近实验值， 完成对核数据的有效调整。（二）背景6 4）近年来随着ENDF/B-6.8及JENDL-3.3等最新微观评价数据库相继释放，这些库中包含大多数重要核素核数据
4、的误差信息协方差数据。数据中心已为CEFR制作了包含15个主要核素的6/12群截面协方差矩阵。（二）背景7（一）国外发展状况二.国内外研究现状 1.发展成熟了利用微扰技术进行核参数对截面的灵敏 度和不确定性分析方法。 目前对核装置采用一阶或二阶微扰技术就足够精确了，而在空间上逐步从一维发展到三维以进一步提高计算精度。2.应用于核工程的许多领域。 考虑核截面本身的不确定对核参数计算结果的影响来提高结果的精确度和可靠程度从而达到提高安全性和经济性的目标，国外开展此项工作涉及如裂变堆芯物理参数分析、燃料循环分析、PWR压力容器延寿、聚变层研究、积分实验前后分析、剂量学及医学应用，油井勘探等方面。本论
5、文关注的裂变堆芯物理参数的误差分析方面，包括如临界keff、控制棒价值、安全棒价值、空泡效应、燃耗反应性损失、有效份额以及多普勒效应等参数。8（1）确定论方法 1）SWANLAKE程序 该程序是美国橡树岭国家实验室（ORNL）在20世纪七十年代末开发的进行屏蔽灵敏度计算的一维程序。它接受ANISN程序提供的通量和共轭通量以及截面数据，可以计算如剂量率或反应率等屏蔽参量对于截面数据的灵敏度。（二）国外一些主要计算工具 二.国内外研究现状 9（1 1）确定论方法）确定论方法 2）TSUNAMI (SCALE-5)程序 SCALE程序包是美国ORNL开发的为核电站许可申请开展相关计算评估的模块程序系
6、统，可以进行包括临界、屏蔽、源项、燃耗、衰变热和传热等项的计算分析。该程序包目前已发展到SCALE5.0版本。在新版本里加入了灵敏度和不确定度分析程序TSUNAMI(Tools for Sensitivity and Uncertainty Analysis Methodology Implementation)可以进行反应堆装置积分参数Keff对相关截面的灵敏度分析和不确定度分析。（二）国外一些主要计算工具 二.国内外研究现状 10（1 1）确定论方法）确定论方法 3）SUSD3D程序 该程序是由斯洛文尼亚的Jozef Stefan研究院通过欧洲“经济合作与发展组织”(OECD)中的原子能机
7、构（NEA）发布的用于核装置对于截面数据灵敏度和不确定分析的计算程序，从初始的SUSD-1D，SUSD-2D已发展到现在的SUSD-3D版本。 该程序采用了一阶微扰理论和计算方法求得灵敏度函数，再利用协方差数据进行探测响应或堆设计参数等积分量对截面数据的不确定度计算。（二）国外一些主要计算工具 二.国内外研究现状 11（1 1）确定论方法）确定论方法 3）SUSD3D程序 新版本中能考虑引起积分参数计算不确定的多种不确定源项，如中子/截面、能量相关的响应函数以及聚变堆中涉及的二次角分布和二次能量分布等不确定项。 SUSD3D通过二进制的交互文件格式能接受大多程序计算输出的通量和共轭通量，这些程
8、序包括DOORS程序包中的ANISN、DORT、 TORT、 ONEDANT、 TWODANT、和THREEDANT等。程序基本结构见图2-1。 （二）国外一些主要计算工具 12图2-1 SUSD3D灵敏度/不确定度分析程序结构示意图13（1 1）确定论方法）确定论方法 4）俄罗斯灵敏度/不确定度分析程序。 俄罗斯在反应堆积分参数对截面的灵敏度和不确定分析方面进行了卓有成效的研究工作，独立开发了一批计算程序，形成了自己的计算系统。其中，包括了不确定计算模块CORE、灵敏度计算程序包、协方差数据处理模块和宏观实验评价等，具体计算系统结构见图2-2，主要灵敏度计算程序见表2-1。积分参数Keff对
9、相关截面的灵敏度分析和不确定度分析。（二）国外一些主要计算工具 14图2-2 俄罗斯灵敏度/不确定度分析程序系统协方差数据灵敏度系数不确定度实验评价CORECORELUNDLUNDLSENSLSENSLEMEXLEMEX15表2-1 俄罗斯灵敏度函数主要计算程序几何确定论方法（近似方法）蒙特卡罗方法子群扩散方法0DNULGEO1DCRABCRAB-12DTVK-2DTWODANT3DTRIGEXMMK-KENONS16（2）蒙特卡罗方法 M-C方法: 利用M-C方法进行灵敏度计算，可灵活处理各种复杂几何结构，因此通常在进行三维灵敏度计算时被采用。但同时为提高计算精度要耗费大量时间，对有些很微小
10、的扰动，其结果可能和计算统计偏差相当而变得没有意义。 1）MCNP4C程序 MCNP4C作为著名的粒子输运计算程序，功能非常强大，它采用了微分算子方法将微扰量转换为响应量的径迹估计而实现随机模拟。该程序可通过1阶微扰或2阶微扰进行积分参数对材料成分、密度、几何结构和截面数据的灵敏度分析。 2）其他 美国KENO程序俄罗斯开发的MMK-KENO都是可用于截面数据对积分量灵敏度和不确定度计算的三维蒙特卡罗程序。（二）国外一些主要计算工具 二.国内外研究现状 17（三）国内发展情况国内发展情况 二.国内外研究现状 在国内有关核截面不确定性引起的积分参数误差的分析只有核数据中心在核截面库协方差矩阵制作
11、和灵敏度计算方面有了一定基础，但把二者相结合用于误差分析还没有开展太多工作。总体说来我们在这方面的研究才刚刚起步，因此有必要尽早投入人力物力开展此项研究工作。18（四）核截面不确定性引起反应堆积分参数误差举例 二.国内外研究现状 BN-600目前的UO2堆芯按计划要过渡到MOX堆芯，为了设计验证需要，俄罗斯IPPE和OKBM研究院做了零功率模拟实验和理论计算，利用TRIGEX扩散程序、MMKKENO蒙特卡罗程序和ABBN-93数据库进行了计算分析，其中利用TRIGEX程序和CORE程序联合计算出了由截面数据不确定性引起的主要堆芯核参数误差数据，参见表2-2。BN-600混合堆芯核参数计算不确定
12、度利用积分实验调整后不确定性调整偏差Keff，%1.50.3+0.1控制棒价值，%5.21.8+0.5安全棒价值，%6.31.9+0.3钠空泡效应，%k/k0.320.11-0.02燃耗，%k/k0.200.12-0.0619表2-2 BN-600混合堆芯主要核参数由截面不确定引起的误差情况（四）核截面不确定性引起反应堆积分参数误差举例核截面不确定性引起反应堆积分参数误差举例 二.国内外研究现状 20三.基本原理和方法图3-1 核装置积分参数不确定度计算流程（评价）（用NJOY加工处理）灵敏度系数计算群截面协方差矩阵（PUFF-2）不确定度计算（如SUSD3D），*计算 ANISN. 微观评价
13、数据库（n：ENDF/B-VI.8, JEFF-3）测量和理论处理（EXFOR）（一）基本流程21 定义：对一个核装置定义一个宏观特 征参数R,它是若干变量如反应截面、裂变谱、权重谱 等的函数，其中一个变量q的单位改变引起R的单位变化量的比值就被定义为特征参数R相应于参数q的灵敏度函数，表示为： 其中 为相空间位置矢量和微分体元。1.灵敏度函数（二）基本概念三.基本原理和方法22（1）协方差基本概念： 随机变量的数学期望：方差为则随机变量 和 之间的协方差为相对协方差可以表示为（二）基本概念2.核截面的协方差数据 23（2）微观点截面协方差矩阵 在微观核截面数据的测量和评价中，各个反应的截面数
14、据也是由更基本的若干直接测量数据根据一定关系导出，这样直接测量数据带有的误差通过函数关系传递给导出量，导致不同点截面之间存在着关联，这种关联就形成了点截面数据的协方差矩阵。举例： 由多个测量值扣除一个共同的常数本底属于这种类型这样 的误差将在导出量 ， 间引入关联，由于由定义得到其中（二）基本概念2.核截面的协方差数据 24（2）微观点截面协方差矩阵 利用 ， 的表达式，可构造出 的协方差矩阵 为（由于协方差矩阵为对称阵，只写出下三角部分）： 2.核截面的协方差数据 （二）基本概念25（3）微观群截面协方差矩阵 从ENDF/B库出发，经过共振重造（RECONR模块）、群平均截面计算（GROUP
15、R模块）后，即把共振参数还原成点截面，再根据具体的能群结构计算出全能区的相应的群平均截面，然后再通过ERROR模块计算出群到群的相对协方差数据。 （二）基本概念2.核截面的协方差数据 26 考虑函数，其中 为直接测量值，误差为 。那么随机变量 是一个导出量， 测量值 的误差 要通过函数关系传递给 ，得到 。一般而言，函数关系可能为线性，也可能为非线性；但在推导误差传递公式时总是将 表述成或近似表述成的 线性函数。实验测量值分布为正态分布，根据正态变量的再现性，可知导出量 必然也是一个正态变量，因此在区间 包含真值的概率是0.683。（二）基本概念3.误差传递 27 举例：则数学期望为于是导出量
16、的测量误差为根据定义式中 为的 方差， 为 ， 之间的协方差。上式就是线性函数的误差传递公式，而推导非线性函数误差传递公式的基本思想是设法将其近似表成线性函数，如用台劳级数展开，再进行线性函数误差传递。（二）基本概念3.误差传递 28（二）基本概念4.积分量不确定度 基本核数据的测量和评价误差会通过函数关系传递给所要计算的反应堆积分参数，形成积分参数计算不确定度。要计算反应堆积分参数的不确定度，首先要根据基础核数据库给出的误差信息库（核数据协方差数据）通过处理程序计算出不同反应的群截面协方差数据库。然后根据堆的结构及几何布置计算出群截面的灵敏度函数，再根据灵敏度函数与群截面协方差数据计算出反应
17、堆积分参数的误差。 29设反应堆的某个积分参数为 ， 而与 计算有关的多群截面数据为 。 可以是堆内某区 某群中子 的某种反应 的平均截面 。函数关系为 。 取一阶扰动则有 ： （二）基本概念4.积分量不确定度 其中30由此得到的相对不确定度其中（二）基本概念4.积分量不确定度 31 提出问题： 通过灵敏度和不确定分析程序可计算得到由截面数据本身的误差带给堆积分参数的不确定度。各类计算表明由核数据不确定引入的误差数值可观，已占积分参数总误差的较大比例。因此，有效降低核数据误差造成的积分参数不确定度是提高反应堆设计精确化的重要内容。 （三）利用积分实验减小不确定度方法 32 解决方法： 目前，比
18、较有效的方法是充分利用各类实验装置开展的积分实验，通过Bayes统计方法将核装置上实验数据和反应堆的计算数据进行相关计算处理，利用实验结果的小不确定度来降低计算结果的大不确定度。 开展这项工作有两个基本前提：一方面实验装置和反应堆上考察的积分参数对主要核数据有近似的灵敏度系数，即二者有较好的相关性；另一方面装置上积分参数的实验测量不确定度必须小于堆积分参数计算得到的不确定度。 具体方法和公式参见论文。（三）利用积分实验减小不确定度方法 331）方法1 确定论方法 方法描述：方法描述： 根据问题的物理性质所建立的数学模型可以用一个或一组确定的数学物理方程来表示，而后对这些方程可以采用数学方法求出
19、其精确或近似的解。堆物理模型求解主要考虑包括能量，位置和方向三个量的数学处理，相应发展了各种近似解法，见表3-1。（四）数值离散方法34自变量E E近似方法分群近似单群双群多群有限差分法粗网格（节块）方法 节点展开法 解析节块法 格林函数节块法有限元方法降维方法综合通量法 近似离散坐标方法积分输运方法 碰撞概率法 面流方法球谐近似方法 扩散近似( )表3-1 堆物理数值求解确定论方法352）方法2 微分算子方法（蒙特卡罗方法） 该方法称为试验统计或“非确定论方法”，它是基于统计（或概率）理论的数值方法，对所要研究问题构造一随机模拟模型，通过计算机进行抽样试验来求得问题的近似解。蒙特卡罗方法特别
20、适合求解本身带有随机性的物理问题如粒输运等。尤其随着计算机的飞速发展，蒙特卡罗方法已在各领域得到广泛应用。 在求解微扰问题上，特别开发了微分算子方法，可以进行材料成分、密度和截面的灵敏度计算。 （四）数值离散方法363）本论文采用方法 SN方法,求解灵敏度函数和不确定度函数 本论文采用Sn方法计算时所需要的通量、共轭通量等参数可由ANISN求解得到，因此论文中进行灵敏度函数和误差函数的求解也采用了SN方法离散形式。大多数积分参数对于截面数据的灵敏度计算的一维、二维程序中均采用了离散坐标方法进行求解，主要考虑了计算时间少和获得较具体的通量分布等优点。在数值离散中，对能量E分成6群（或12群），空
21、间位置采用一维有限差分，方向采用SN离散坐标方法。具体离散和推导见附录二。（四）数值离散方法37 本程序采用Fortran90语言编写,选用了微软公司开发的Fortran PowerStation4.0 版本编译系统进行编译。FortranPowerStation4.0采用32位计算，工作在Windows平台，不仅完全支持Fortran 90，也支持与VC+等的混合编程，完全满足本程序的编译要求。四.程序开发（一）程序开发环境38程序调试灵敏度验证计算keff相应各截面灵敏度函数推导keff灵敏度函数数值离散keff灵敏度函数子程序编写MCNP灵敏度计算不确定度子程序编写灵敏度/不确定度程序调
22、试灵敏度/不确定度程序初步验算ANISN计算（ ）群截面协方差数据接口子程序编写四.程序开发（二）程序开发流程图4.1 SUCA1D一维程序开发流程图39四.程序开发（三）程序主要公式 (1)keff对于各反应截面的灵敏度函数 0 BA40四.程序开发（三）程序主要公式 (1)keff对于各反应截面的灵敏度函数 中子消失项（泄漏和移出）中子散射到本相空间项41四.程序开发（三）程序主要公式 (1)keff对于各反应截面的灵敏度函数 中子裂变产生项 系统内中子裂变积分项 42x序号类型公式内容备注1裂变截面2俘获截面 34 5散射截面 等散射6总截面 DTSizgakizga, 1,DTTSiz
23、gfizgfkizgf, 3, 1,DTSizgkizg, 3,DTSizgkizg, 3,DTTSizgsizgskizgs, 2, 1,DTTTSizgfizgsizgtkizgt, 3,2, 1,31,TT为中 ),(fn1T为),(),(),(anpnn3T为中中x3T中x为，与f互换31,TT中x),(),(nnnn321,TTT分别为fst,表4.1 keff对各反应核数据的灵敏度计算公式43四.程序开发（三）程序主要公式 (3)keff相关截面数据的不确定度44 本程序SUCA1D是基于ANISN计算提供的通量、共轭通量等基本参数和NJOY处理得到的协方差数据来进行堆积分参数ke
24、ff对各截面的灵敏度和不确定度一维计算程序，主要包含了基本参数读入模块、宏观截面处理模块、灵敏度计算模块、不确定计算模块和输出模块等，目前可以计算包括裂变截面、吸收截面、裂变中子、裂变谱、散射截面及总截面对应keff的灵敏度系数和不确定度。 SUCA1D程序主要包含了12个子程序，共约1200行。程序的计算框图见图4.2，程序主要模块功能介绍见表4.2。四.程序开发（四）程序系统及特点45图4.2 SUCA1D程序计算流程46表4.2 SUCA1D程序主要模块介绍序号模块名称内容介绍1FREADFREAD程序主要从三个文件中读入输入数据。其中从ANISN计算后分别输出的通量和共轭通量两个文件中
25、直接读入群通量、角通量、微观截面、差分节点体积、密度、求积权重系数、有效增殖因子、裂变谱等主要参数。另外一个为编写的输入文件，其主要包括相对协方差数据及输出控制参数。2MACROSSMACROSS由于本程序计算需要单独的各个核素在不同区域的宏观截面，而ANISN一般输出中都是混合宏观截面，因此需要将ANISN输出中的各材料区中每个核素的微观截面计算处理为宏观截面，这样避免了直接调用微观截面的繁琐并节省了内存。3SPERTSPERT调用FREAD中读入的基本数据和MACROSS处理后的宏观截面，依据上面给出的灵敏度计算公式，编写了通用的灵敏度计算模块，供具体计算各不同反应截面灵敏度时调用。4SR
26、ESULTSRESULT根据各个不同反应类型的灵敏度计算要求确定输入参数，调用SPERT模块，得到灵敏度计算结果。5ERRORERROR调用SRESULT计算得到的灵敏度系数和FREAD中读入的协方差数据，依据误差公式计算得到不同核素基于区域、能群和截面的keff的误差。6FOUTFOUT按照输入卡的控制要求，可输出按区域、能群和核素等微分或积分的灵敏度系数和keff的误差数据。47 中国实验快堆（CEFR）是一座热功率为65 MW，电功率为25MW的池式钠冷快中子反应堆，计划于2009年达到首次临界。CEFR运行初期采用高富集度的二氧化铀燃料，以后逐步过渡到MOX燃料。本文计算是针对CEFR
27、平衡态堆芯进行的。堆芯装载图见图5.1，主要参数见表5.1。五.程序初步验证（一）CEFR堆芯介绍48图5.1 CEFR平衡态堆芯布置图（一）CEFR堆芯介绍49表5.1 CEFR平衡态额定功率下堆芯主要参数表名 称燃料类型235U重量含量，燃料组件数，盒UO264.481组件对边距，cm组件栅元面积，cm2堆芯等效半径，cm 30.2活性段，cm 45.6过渡区，cm 8.6 上/下转换区，cm 10.1/25.1上屏蔽区/下气腔，cm 506.13732.62内 容（一）CEFR堆芯介绍50五.程序初步验证（二）程序主要输入参数1）一维模型（参考图5.2 R-Z 模型）图5.2 CEFR堆
28、芯R-Z模型结构示意图51五.程序初步验证（二）程序主要输入参数2）能群结构表5.3 程序计算采用的能群结构能群序号能量上限，MeV能群中点，MeV1 158.177E+002 1.3539.255E-013 4.979E-12.694E-014 4.087E-22.499E-025 9.119E-34.787E-036 4.540E-42.270E-04能群序号计算通量共轭通量ANISNMCNPANISN10.1600.1790.19620.2620.2540.19630.4800.4750.18040.0650.0740.15050.0270.0170.13860.0060.0010.14
29、152五.程序初步验证（二）程序主要输入参数3）通量和共轭通量表5.4 CEFR堆芯6群归一化能谱53五.程序初步验证（二）程序主要输入参数3）通量和共轭通量1E-41E-30.010.11100.00.10.20.30.40.5Relative FluxE, Mev ANISN MCNP图5.3 CEFR一维堆芯归一化能谱10-410-310-210-11001010.130.140.150.160.170.180.190.20Relative AdfluxE, Mev ANISN图5.4 CEFR一维堆芯归一化共轭通量谱54五.程序初步验证（二）程序主要输入参数4）协方差数据 核数据中心可
30、通过NJOY程序系统将微观评价库ENDF/B中的协方差数据制作成供反应堆积分参数不确定度分析使用的多群协方差矩阵，现已为CEFR制作了包含了15个重要核素的6群和12群相对协方差数据。 本论文结合灵敏度计算情况和已有的协方差数据，选取了235U、238U、Fe和Na核素裂变、俘获等重要核素反应截面的相对协方差矩阵 。55五.程序初步验证（三）计算结果1）灵敏度系数 主要计算得到了CEFR堆芯中235U、238U、Fe和Na核素的相应裂变或辐射俘获反应灵敏度系数，并绘制了灵敏度曲线，分别参见表5.10-5.12和图5.5-5.10。56五.程序初步验证（三）计算结果1）灵敏度系数1 235U裂变
31、、辐射俘获表5.10 keff对235U主要反应截面灵敏度能群序号SUCA1DMCNPSUCA1DMCNP11.09E-011.34E-01-2.20E-03-2.30E-0321.31E-011.69E-01-1.01E-02-8.60E-0332.77E-013.80E-01-4.02E-02-4.42E-0245.73E-029.20E-02-9.10E-03-1.85E-0252.32E-024.74E-02-3.80E-03-1.14E-0261.66E-039.70E-03-3.00E-04-2.90E-03总灵敏度5.99E-018.33E-01-6.57E-02-8.78E-0
32、2235U235U57五.程序初步验证（三）计算结果10-410-310-210-11001010.00.10.20.30.4Sensitivity coefficientE, Mev SUCA1D MCNP图5.5 keff对235U裂变反应截面灵敏度曲线10-410-310-210-1100101-0.05-0.04-0.03-0.02-0.010.00Sensitivity coefficientE, Mev MCNP SUCA1D 图5.6 keff对235U 辐射俘获截面灵敏度曲线1）灵敏度系数1 235U裂变、辐射俘获58五.程序初步验证（三）计算结果1）灵敏度系数2 238U裂变
33、、辐射俘获能群序号SUCA1DMCNPSUCA1DMCNP12.49E-023.45E-02-7.68E-04-9.60E-0427.54E-041.20E-03-4.92E-03-4.03E-0331.06E-052.00E-05-1.10E-02-1.17E-0241.26E-060.00E+00-3.31E-03-6.18E-0354.05E-070.00E+00-1.37E-03-3.01E-0369.74E-100.00E+00-2.16E-04-3.80E-04总灵敏度2.56E-023.57E-02-2.16E-02-2.63E-02238U238U表5.11 keff对238U
34、主要反应的灵敏度10-410-310-210-11001010.000.010.020.030.04Sensitivity coefficientE, Mev SUCA1D MCNP59五.程序初步验证（三）计算结果1）灵敏度系数2 238U裂变、辐射俘获10-410-310-210-1100101-1.5x10-2-1.0 x10-2-5.0 x10-30.0Sensitivity coefficientE, Mev MCNP SUCA1D 图5.8 keff对238U 辐射俘获截面灵敏度曲线图5.7 keff对238U 裂变反应截面灵敏度曲线60五.程序初步验证（三）计算结果1）灵敏度系数
35、3 keff对Fe 、Na辐射俘获反应表5.11 keff对Fe 和Na核素辐射俘获截面的灵敏度 能群序号SUCA1DMCNPSUCA1DMCNP1-1.45E-03-2.00E-04-2.40E-04-1.00E-052-1.25E-03-1.29E-03-3.93E-05-3.00E-053-4.35E-03-3.59E-03-2.06E-04-1.70E-044-1.10E-03-1.25E-03-3.08E-06-6.00E-055-2.61E-03-1.18E-03-1.04E-04-1.90E-046-5.51E-04-5.90E-04-1.40E-05-3.00E-05总灵敏度-
36、1.13E-02-8.10E-03-6.06E-04-4.90E-0456Fe23Na61五.程序初步验证（三）计算结果1）灵敏度系数3 keff对Fe 、Na辐射俘获反应10-410-310-210-1100101-5.0 x10-3-4.0 x10-3-3.0 x10-3-2.0 x10-3-1.0 x10-30.0Sensitivity coefficientE, Mev MCNP SUCA1D 图5.9 keff对Fe辐射俘获截面灵敏度曲线10-410-310-210-1100101-4.0 x10-4-3.0 x10-4-2.0 x10-4-1.0 x10-40.0Sensitivi
37、ty coefficientE, Mev MCNP SUCA1D 图5.10 keff对Na辐射俘获截面灵敏度曲线62五.程序初步验证（三）计算结果2）不确定度计算结果 序号类 型数值，11.2722.203 0.1514 0.73150.126 0.0047总不确定度2.658参考总不确定度* 1.902235U-238U-238U-Fe-Na-235U- 表5.13 CEFR主要核素反应截面引起的keff不确定度 63五.程序初步验证（三）结果分析 （1）从灵敏度计算公式可以看出，积分参数keff对某个截面的灵敏程度主要与核素截面、反应堆通量和共轭通量等相关。其中共轭通量的物理意义为中子价
38、值，从5.2节中ANISN计算得到的共轭通量谱分布来看，各群中子对反应堆keff的价值相差不大。因此灵敏度大小主要由截面和通量来确定。 64五.程序初步验证（三）结果分析 （2）235U裂变反应的积分灵敏度最大，反应了其维持反应堆临界的主导地位。而相比之下，辐射俘获反应的灵敏度要小于裂变反应的。 65五.程序初步验证（三）结果分析 （3）比较各核素截面的群灵敏度可以看出，235U裂变、俘获反应和238U、Fe及Na的俘获反应的第3群灵敏度最大，是积分灵敏度中的主要部分。这主要是因为对于各群截面相近的反应，群通量大小决定了其灵敏度结果，而CEFR堆芯中子通量主要集中在平均能量为100keV的第三
39、能群，使得上述核素反应的第三群灵敏度最大。66五.程序初步验证（三）结果分析 （4）238U裂变反应灵敏度计算结果表明，灵敏度主要贡献来自第1群。这是因为238U裂变反应是主要发生在1MeV以上的阈能反应，其反应截面占了主导因素，群通量为次要因素，所以代表高能区的第1群的灵敏度数值最大。67五.程序初步验证（三）结果分析 （5）从SUCA1D和MCNP程序对几个主要核素截面灵敏度结果比较来看，二者计算结果趋势基本一致，数据接近。其差别可能主要来自MCNP采用了ENDF/B5点截面库和三维模型进行蒙特卡罗模拟计算，而SUCA1D是基于NVITAMIN-C的6群库和一维模型进行计算。68五.程序初
40、步验证（三）结果分析 （6）235U裂变和辐射俘获反应截面是引起的积分参数keff不确定度的主要来源。积分参数不确定度大小由截面灵敏度和协方差数据两方面因素来确定，从灵敏度分析得到，235U的裂变反应灵敏度是keff对各截面中的最大灵敏度，但由于辐射俘获反应的协方差数据远大于裂变反应的，因此最终辐射俘获反应截面引起的keff不确定度大于裂变反应的贡献。69五.程序初步验证（三）结果分析 （7）238U裂变和辐射俘获反应截面引起积分参数keff不确定度的近似计算结果表明，其不确定度影响是处于次要地位的。由于缺乏对应的协方差数据，238U裂变和辐射俘获反应截面是引起的积分参数keff不确定度是近似
41、用235U的协方差数据代替，计算结果表明238U裂变和辐射俘获反应截面对应的keff不确定度要比235U的小将近一个量级。这主要是由238U裂变和辐射俘获反应截面较小的灵敏度系数决定的。70五.程序初步验证（三）结果分析 （8）本论文程序计算出CEFR堆芯物理计算系统中由核截面误差引起的keff不确定度约为2.65，而俄罗斯程序系统得到的不确定度为1.9，二者基本接近，初步反映了两套程序系统具有一定的可比性，这也通过两个系统并行完成的CEFR设计计算得到了验证。71六.结论结论一： 本论文在调研了相关文献基础上验证推导了堆系统积分参数keff由核数据引起的灵敏度和不确定度计算公式，并根据国内物
42、理屏蔽计算程序系统的特点（输运计算、方法和ANISN格式等），对灵敏度公式进行了相应的离散，获得了可靠的理论公式，为编写程序打下基础。同时通过调研学习，掌握了对反应堆系统积分参数利用微扰技术进行核数据灵敏度和不确定度分析的基本方法。 72六.结论结论二： 本论文在国内首次建立了专门用于核数据灵敏度和不确定度分析的一维程序，并利用现有的协方差数据库，首次给出了CEFR堆芯keff由235U裂变等重要反应截面引起的不确定度。73六.结论结论三： SUCA1D程序目前主要考虑了由裂变截面、俘获截面、(n,p)、(n,a)、裂变谱分布、裂变中子产额以及总截面等参数对keff的灵敏度和不确定度的计算。程
43、序可以直接从ANISN输出的通量和共轭通量两个文件以及编写的协方差输入卡中读入相关参数，进行上述计算内容的计算，并且根据输出要求可分别按能量、空间和核素进行微分或积分结果的输出。 74六.结论结论四： 通过和MCNP程序一同对“CEFR平衡态装载堆芯keff对应核数据的灵敏度”的对比计算，得到了比较相近的计算结果，并且结果也符合理论规律。而对几个重要核素的主要截面引起keff总误差和俄罗斯计算结果的比较，可看出二者基本接近。综合上述两方面的初步比对计算，可以认为SUCA1D程序基本正确，计算结果合理可信。75七.展望展望一： 程序功能和界面的完善。目前SUCA1D程序只采用了简单的用户交互方式
44、,程序中也没有绘图功能。为了更方便使用该程序，在后续工作中可以改善交互界面，并在程序中增加灵敏度和不确定度对于能群、空间节点及核素的函数曲线绘图功能，这样可方便直观地观察计算结果。同时这些程序模块的建立也为以后开发2维、3维程序提供基本的技术支持。76七. 展望展望二： 开展其他堆积分参数的灵敏度和不确定度分析研究。要全面对反应堆设计中堆芯物理参数进行由核数据引起的误差分析，除了计算keff误差外，还有必要考虑如控制棒价值、钠空泡反应性效应、燃耗反应性效应等参数的误差情况。77七.展望展望三： 开发堆积分参数对核数据灵敏度和不确定度的二维、三维计算程序。为了进一步提高计算的准确度，有必要利用现
45、有程序系统，开发二维和三维的灵敏度计算程序。二维程序系统中可以考虑利用现有的DORT二维Sn计算程序来进行通量和共轭通量的计算。而三维程序开发可同时探索用TORT等确定论程序和MCNP等非确定论程序进行主要积分参数的灵敏度分析。78七.展望展望四： 程序系统的建立。开展反应堆系统积分参数对核数据的灵敏度和不确定度计算主要涉及协方差数据制作、灵敏度和不确定程序开发和通量等基本参数计算三方面。因此要开展好这方面的工作，必须要根据条件许可选择合适的技术路线，处理好三方面的接口衔接，组成一个方便、高效和可靠的计算系统。尤其是在协方差数据制作上，目前从基础数据库中提取的截面协方差数据还较少，仍缺乏关键截面的协方差数据。79七.展望展望五： 开展利用相关实验结果减小积分误差的研究。利用灵敏度和不确定计算系统，可以获得核数据不确定性传递到堆系统积分参数所带来的误差情况，由于核数据误差是由实验和数据处理过程中的统计误差、系统误差等组成，要单独减小这些微观数据的误差并非易事。而利用在相近核装置上进行的一些相关实验的结果可以对计算获得的堆系统积分参数的不确定度进行修正处理，从而降低了该积分参数由核数据不确定性带来的误差。因此有必要开展这方面的技术研究工作，为以后可能利用CEFR、东风六号快堆临界装置或俄罗斯的堆芯零功率模拟装置的相关实验数据进行积分参数的误差降低处理打下基础。
