算法构建了一个核反应堆的数字孪生监控方法，通过点堆动力学方程与热传递方程耦合的物理模型精确模拟反应堆功率、先驱核浓度和温度的三维状态演化；在正常工况下采集包含测量噪声的模拟数据，训练支持向量机单类分类器学习物理模型预测与传感器测量之间残差的正常分布模式；在实际运行监控中，数字孪生持续进行物理模型前向预测并与实时测量比较生成残差特征，通过预训练的SVM模型检测异常残差模式，同时应用持续性滤波器避免瞬时噪声误报；该系统实现了物理知识驱动与数据驱动融合的混合监控，能够区分正常噪声、缓变故障（冷却能力下降）和突发事故（反应性引入）等多种异常模式。

详细算法步骤：

第一步：核反应堆物理模型构建与求解

基于点堆动力学六群缓发中子模型，建立功率与先驱核浓度的耦合微分方程组

引入热传递方程描述反应堆温度动态变化，考虑功率产热与冷却散热平衡

采用四阶龙格-库塔数值积分方法求解非线性微分方程组，确保数值稳定性

设定合理的数学极限值防止计算溢出，保证长时间模拟的数值可靠性

第二步：数字孪生系统架构设计

设计混合监控系统类，整合物理模型与机器学习组件

物理引擎部分实现基于当前状态的前向预测功能

安全监控模块实现硬阈值检查，当功率超过安全限值时触发停堆

残差分析模块计算传感器测量与物理预测之间的差异特征

第三步：支持向量机分类器训练与校准

在纯正常工况下运行物理模型生成仿真数据，模拟传感器测量噪声

计算每个时间步的功率残差和温度残差，构建正常状态的特征空间

对残差特征进行标准化预处理，消除量纲影响

配置支持向量机单类分类器参数，采用小nu值提高异常检测灵敏度

使用正常残差数据集训练SVM模型，学习正常操作空间的边界

第四步：多场景测试数据生成

设计3600秒长时间测试场景，分为三个阶段模拟不同运行状态

第一阶段保持正常工况，验证基础监控性能

第二阶段模拟冷却系统缓慢退化，冷却系数降低20%

第三阶段模拟反应性引入事故，外部反应性阶跃增加

每个阶段保持一致的传感器噪声水平，确保测试真实性

第五步：实时监控与异常检测流程

初始化数字孪生状态与实测值同步，消除初始偏差

逐时间步执行物理模型预测，生成下一时刻状态估计

应用同步滤波器融合物理预测与当前测量，更新内部状态估计

计算功率和温度残差，构建二维特征向量

将标准化后的残差特征输入训练好的SVM模型进行分类

应用持续性滤波逻辑，要求连续多个时间步检测到异常才触发警报

并行运行硬阈值安全监控，当功率超过65MW时立即触发停堆

第六步：可视化与性能评估

绘制三面板综合展示图表，直观呈现系统运行状态

上层面板显示实测功率曲线，用颜色区域标注不同事件阶段

中间面板展示残差时间序列，可视化模型预测误差的演变

下层面板叠加显示AI预警信号和安全系统停堆决策

通过时间对齐分析，评估混合系统对缓变故障和突发事故的检测能力

验证持续性滤波对减少误报、提高检测可靠性的效果

参考文章：

基于点堆动力学-热传递耦合物理模型与支持向量机残差分析的核反应堆数字孪生混合异常检测算法（以模拟信号为例，Python） - 哥廷根数学学派的文章
https://zhuanlan.zhihu.com/p/1998375229928527716