本文档发表于《欧洲物理期刊A》（Eur. Phys. J. A）2015年第51卷第182期，由瑞士保罗谢勒研究所（PSI）的D. Rochman与荷兰核研究咨询集团（NRG）的A.J. Koning和S.C. van der Marck共同完成。这是一篇关于改进中子学（neutronics）模拟与核数据不确定度量化方法的原创性研究论文。研究领域属于核数据评价与反应堆物理计算。

研究背景与目标 核数据（nuclear data），如中子截面、裂变谱等，是进行反应堆设计、临界安全分析等中子学模拟的基础。核数据评价者的目标是提供与所有可用测量数据（包括微分测量和积分测量）及理论预期最吻合的“最佳估计”值及其不确定度。微分测量（differential measurements）针对单一核数据（如某个能量点的裂变截面），而积分测量（integral measurements），例如大量临界基准实验（criticality benchmarks）的有效增殖因子（keff）测量，则能综合反映一系列核数据的整体质量。通常，通过比较模拟计算值（C）与实验测量值（E）的比值（C/E）或计算χ²值，可以评估不同核数据库（如美国的ENDF/B、OECD/NEA的JEFF、日本的JENDL）的性能。

然而，如图1所示，尽管随着时间的推移，各大核数据库的χ²值总体呈下降趋势（即与积分基准实验的符合度在提高），但近年来的改进已变得有限。此外，传统的“全蒙特卡洛”（Total Monte Carlo, TMC）不确定度量化方法，即使用大量随机核数据文件进行重复计算，虽然能提供保守的不确定度估计，但计算成本高昂，且未充分利用积分基准实验的信息来“约束”或“优选”核数据。

因此，本研究旨在提出一种新颖且相对简单的方法，以同时实现两个目标：1) 通过有限次数的计算，显著改善中子学模拟（如临界计算）与积分基准实验的符合度；2) 通过将积分基准实验的信息整合到模拟方案中，降低因核数据不确定度带来的模拟结果不确定度。研究以²³⁹Pu核数据为例进行演示。

详细研究流程 本研究主要包含两个核心部分：随机核数据文件的生成与基于积分基准实验的筛选和加权组合。

第一部分：随机²³⁹Pu核数据文件的生成 1. 方法：研究团队使用核反应理论计算系统TALYS，生成了10,000个随机版本的²³⁹Pu核数据ENDF文件。这些随机文件的产生并非直接拟合实验数据，而是以性能较好的JEFF-3.2库的评价为参考中心值，围绕其在一定范围内随机扰动模型参数。这种方法既节省了重新分析海量微分实验数据的时间，又能确保生成的随机截面在物理上是合理的，并覆盖了足够的“核相空间”。 2. 处理：使用NJOY程序将所有随机ENDF文件处理成可用于MCNP蒙特卡洛模拟程序的ACE格式截面库。 3. 特点：图3展示了部分随机文件的裂变截面及其不确定度。这些随机文件之间的分散度，既考虑了实验和理论的不确定度，也允许了更大的变化范围，以探索更广的可能性。表1对比了主要数据库与随机文件平均值的若干关键核数据（如热中子裂变截面、俘获截面、平均裂变中子数等），显示了随机文件的统计特性。

第二部分：基于积分基准实验的筛选与加权 1. 基准实验选择：从国际临界安全基准实验手册（ICSBEP）中，选取了35个以²³⁹Pu为主的临界基准实验（其中22个快谱，13个热谱），如表2所示。选择标准基于计算资源的限制（需为每个随机文件计算每个基准）和案例的覆盖面。图4展示了三个主流数据库（ENDF/B-VII.1, JENDL-4.0, JEFF-3.2）对这35个基准的C/E值表现，其χ²值分别为0.8, 2.3和0.6。 2. 单个随机文件的性能评估（“良好拟合”）： * 流程：将每一个随机²³⁹Pu文件（替换掉原始库中的²³⁹Pu数据）分别与三个主流库（用于提供其他核素的数据）结合，使用MCNP计算全部35个基准的keff。 * 评估：为每个“库+随机文件”组合计算其相对于35个基准的χ²值（公式1）。χ²值越低，表示该组合与基准实验的整体符合度越好。 * 结果：对于每个原始库，都能找到一个最佳的随机文件（如ENDF/B-VII.1对应文件#6174，JENDL-4.0对应#912，JEFF-3.2对应#5989），使其χ²值显著低于原始库本身（见表3）。例如，ENDF/B-VII.1的χ²从0.8降至0.26。图5-7直观展示了使用最佳随机文件后，C/E值更集中地分布在1.0附近，改善明显。 * 发现：一个有趣的现象是，针对某个库找到的最佳随机文件，与其他库搭配时也能取得不错的效果（见表3“平均χ²”列），这表明这些优选文件本身具有普适的优良特性。 3. 多个随机文件的加权组合（“更好拟合”）： * 动机：为了寻求比单一最佳文件更好的整体符合度，并为进一步的不确定度量化铺路，研究者探索将多个随机文件的结果进行线性组合。 * 方法：组合后的keff计算值由公式2给出：C_i = (Σ w_j * keff_ij) / (Σ w_j)，其中w_j是第j个随机文件的权重。 * 权重策略探索：研究比较了三种策略：(1) 所有文件等权重（即传统TMC方法）；(2) 权重为exp(-χ²_j)；(3) 选取少量（m个）最佳文件，并通过优化算法寻找一组特定的正权重{w_j}，使得组合结果的χ²最小化。 * 优化流程：对于策略(3)，采用迭代算法：随机尝试添加一个文件，并随机搜索其权重，看是否能降低总体χ²；如果不能则丢弃该文件。重复此过程直至χ²不再显著降低。 * 结果：策略(3)被证明最有效且实用。如图8所示，对于ENDF/B-VII.1库，仅需7个精心挑选的随机文件及其特定权重，即可将χ²降至约0.17（表4），且继续增加文件数对改善χ²贡献甚微。类似地，JENDL-4.0需要6个文件，JEFF-3.2需要9个文件。这实现了以极少的额外计算量（仅需运行m次模拟）获得最优拟合。 4. 方法的外推验证（预测能力）： * 流程：为了测试方法的泛化能力，研究者将上述优选出的“库+文件组合”应用于更大的基准实验集——400个²³⁹Pu临界基准实验（即包含未用于筛选的365个基准）。 * 结果：如图9和表5所示，无论是使用单个最佳文件还是多个文件的加权组合，其χ²值均显著优于原始库。例如，ENDF/B-VII.1的χ²从3.8分别降至2.9（单文件）和2.5（7文件组合）。这证明了该方法不仅对用于筛选的基准有效，也具有良好的预测能力。 5. 不确定度的降低： * 流程：利用优选出的m个加权随机文件，不仅可以得到改进的中心值（keff），还可以根据公式3计算其标准偏差，作为核数据不确定度导致的keff不确定度。公式3确保了中心值与不确定度计算的一致性。 * 结果：图10展示了对于400个基准，使用7个加权文件（策略3）计算的不确定度，与使用10,000个等权重文件（策略1，传统TMC）计算的不确定度对比。前者在几乎所有基准案例中都显著降低了不确定度。这是因为加权组合融入了积分基准实验的信息，过滤掉了那些与基准实验严重不符的随机文件，从而提供了“最佳估计加不确定度”，而非保守的“保守估计加不确定度”。 * 应用示例：图11展示了一个MOX燃料栅元（pincell）燃耗计算示例。使用7个加权随机文件得到的keff不确定度带，远窄于使用10,000个等权重文件得到的不确定度带，清晰地展示了将临界基准信息传递到实际应用中所带来的不确定度缩减效益。 6. 对核数据评价的反馈： * 流程：分析那些能带来更好积分性能的随机文件（或其组合）的微分核数据特征，可以为核数据评价者提供改进方向的线索。 * 结果与推断：图12-14对比了优选文件与原始库的核数据差异。研究指出了一些可能有助于改善积分结果的微分数据调整趋势，但强调这些结论是“指示性而非决定性的”，因为可能存在不同核素数据间的补偿效应。主要发现包括： * 热中子裂变截面可能需要增加约3.2%，但关键共振区（如0.29 eV）变化不大。 * 热中子平均裂变中子数（ν）变化不大，但高于1 keV能量后有略微降低的趋势。 * 热中子入射下的平均裂变中子谱（PFNS）能量可能需增加约80 keV（谱“变硬”），而0.5 MeV入射时则可能需降低约700 keV（谱“变软”）。 * 非弹性散射截面（n, inl）可能需要显著降低。 * 热中子俘获截面未见明显调整需求。 * 局限性：作者明确指出，该方法目前仅针对临界基准实验进行优化，要形成通用性评价，还需结合其他类型应用（如屏蔽、衰变热）和微分实验数据进行综合判断。同时，缺乏灵敏度向量（sensitivity vectors）分析也限制了对特定截面变化重要性的精确解读。

结论与意义 本研究提出了一种通过筛选和加权组合随机核数据文件来改进中子学模拟与核数据不确定度量化的创新方法。其核心价值在于： 1. 提升模拟精度：该方法能有效利用大量积分基准实验的信息，显著改善核数据库对临界基准实验的再现能力（降低χ²），如图15所总结。 2. 降低计算成本与不确定度：仅需运行少数几次（如7次）加权模拟，即可同时获得改进的中心值估计和更具现实意义的（更小的）核数据不确定度，实现了从“保守估计”到“最佳估计加不确定度”的转变。 3. 为评价提供指引：方法产生的优选核数据文件，能为核数据评价者提供基于积分性能的微分数据调整线索。 4. 方法通用性强：该框架可扩展至其他重要核素（如²³⁵U, ²³⁸U），并整合更多类型的积分实验数据，从而有望系统性提升通用核数据库的质量。

研究亮点 1. 方法新颖且高效：将随机核数据生成技术与基于积分实验的统计筛选、加权组合相结合，构思巧妙。以极少的计算代价（数倍于单次模拟）同时实现了精度提升和不确定度量化。 2. 理念先进：明确提出了将积分实验信息作为“约束”融入核数据不确定度量化流程，使最终的不确定度估计既基于微分数据的先天分散性，也受益于积分实验的后验验证，更具实际应用价值。 3. 实证充分：研究以²³⁹Pu为例，从少量基准筛选、到多基准验证、再到实际栅元应用，层层递进，完整展示了方法的可行性和有效性。 4. 兼具实用性与理论价值：所提方法可直接应用于工程设计和安全分析，以降低设计裕量；同时也为核数据评价工作提供了新的分析工具和思路。

其他有价值内容 文章还讨论了该方法的潜在局限性，如可能存在核素间补偿效应、对线性组合假设的依赖、以及需要结合灵敏度分析以获得更可靠结论等，体现了作者严谨的科学态度。同时，也展望了该方法与更传统评价方法相结合、以及应用于更广泛基准实验类型的前景。