《Annals of Nuclear Energy》：Differential equation for reactivity monitoring using the superposition principle
编辑推荐：
本综述提出了一种基于逆点动力学模型的新型反应性计算方法，该方法无需核功率历史数据即可实现实时反应性估算。通过将总反应性解耦为瞬发中子贡献项和六组缓发中子贡献项之和，并推导出六个一阶常微分方程（ODEs），实现了对缓发中子效应的独立求解。该方法有效克服了传统积分模型在系统重启或计算中断时依赖历史功率数据的局限，提升了反应性监测的连续性与响应性，为新型核反应堆（特别是具有较长中子代时间的反应堆）的实时安全监控提供了新思路。
在核反应堆的安全监控与运行控制中，反应性（ρ）是衡量堆芯偏离临界状态程度的关键参数。传统的反应性估算方法，特别是基于经典逆点动力学模型的积分形式，严重依赖于对核功率历史数据的连续采集与积分计算。这在实时监控系统中存在明显短板：一旦数据采集或计算过程中断，系统的重新启动需要等待积累足够的历史数据，才能恢复准确的计算，从而造成反应性监测的空窗期，不利于反应堆的连续安全监控。此外，高频率的数据采样需求也增加了系统对噪声的敏感性。因此，开发一种不依赖于核功率历史、能够快速响应并适用于更广泛反应堆类型（尤其是新兴反应堆技术）的新型反应性计算方法，成为近年来的研究重点。
积分形式的逆点动力学模型及其局限
经典的逆点动力学模型是大多数反应性计的基础。该模型通过求解一个积分方程，从测量的核功率P(t)中反推出反应性ρ(t)。其表达式综合考虑了瞬发中子的贡献（与功率变化率相关）和六组缓发中子先驱核的贡献（通过积分项体现）。该模型虽然成熟，但其核心缺陷在于积分项∫P(τ)dτ。这项计算意味着当前时刻的反应性估值，强烈依赖于从初始时刻到当前时刻的全部核功率历史。任何计算的中断都意味着历史数据的丢失，重新计算需要时间重新积累数据，导致反应性计在启动或重启后存在一段不可用期。
新型微分方程方法：解耦反应性贡献
为了克服上述局限，本文提出了一种创新的视角和计算方法。其核心思想是将总反应性视为瞬发中子贡献（ρP(t)）与所有缓发中子贡献之和（ρD(t)）。其中，缓发中子贡献又可依据叠加原理，分解为六组不同半衰期和份额的缓发中子先驱核的独立贡献之和（ρD(t) = ΣρDi(t), i=1,...,6）。
基于这一分解，研究团队从传统的积分形式出发，经过一系列数学变换，成功地为每一组缓发中子先驱核的贡献ρDi(t)推导出了一个独立的一阶常微分方程（ODE）。这个ODE的形式为：
dρDi/dt + ( (1/P) dP/dt + λi) ρDi- (βi/P) dP/dt = 0
其中，λi是第i组缓发中子先驱核的衰变常数，βi是其对应的有效缓发中子份额。
这个方程的关键优势在于，它完全消除了对核功率历史积分项的依赖。要计算当前时刻的ρDi，只需要知道当前的核功率P(t)和功率变化率dP/dt，以及反应堆的动力学参数（Λ, βi, λi）。
方法实现与优势
在实际计算中，反应性估算按以下步骤进行：
- 1.
瞬发贡献计算：直接由公式 ρP(t) = (Λ/P(t)) dP/dt 计算。
- 2.
缓发贡献计算：并行或依次求解六个上述的一阶ODE，得到每一组的ρDi(t)。论文中采用了经典的四阶龙格-库塔法（RK4）进行数值求解。然后将六组结果相加，得到总的缓发中子贡献ρD(t)。
- 3.
总反应性合成：总反应性由两者相加得到：ρ(t) = ρP(t) + ρD(t)。
该方法带来了多重显著优势：
- •
无需历史数据：最大的突破是摆脱了对核功率历史数据的依赖。反应性计可以在任何时刻启动或重启，无需等待数据积累，立即开始提供有效的反应性估算，极大地增强了系统的可用性和鲁棒性。
- •
适用性更广：与之前一些基于微分方程的方法相比，本文提出的方法对中子代时间（Λ）没有限制（例如要求Λ ≤ 10-5s），因此尤其适用于中子代时间较长的新型反应堆（如某些加速器驱动次临界系统或熔盐堆）。
- •
采样灵活性：该方法允许以较长的间隔（如1秒）采样核功率数据，降低了对数据采集系统实时性的苛刻要求，也减少了高频采样带来的噪声干扰。
- •
物理图像清晰：将反应性解耦为瞬发和缓发贡献，有助于操作员更清晰地理解瞬变过程中不同中子成分的作用。例如，在缓发中子主导的瞬变中，ρP的贡献很小；而在瞬发中子主导的快速瞬变中，ρD的贡献则占主导。
结论
本文提出的基于逆点动力学微分方程模型的无历史反应性计算方法，是核反应堆实时监控技术的一项重要进展。它通过巧妙的数学重构，将原本依赖全局历史信息的积分问题，转化为一组仅依赖当前状态的微分方程求解问题。这不仅解决了传统反应性计在中断恢复后的数据空窗问题，也提升了对多种堆型的适应性。该方法的成功实施，有望为下一代核反应堆，特别是那些安全运行窗口更小、对监控实时性要求更高的新型反应堆，提供更加可靠、灵活和快速的反应性监测工具，为核能的安全利用增添一道重要的技术保障。
