1. 概述
核电仿真模型是实现孪生核电站的基础。真实核电机组的不同设备、系统对应不同的仿真模型，这就导致仿真建模需要具备设备多样、参数多样、属性多样的复杂关系属性。建立核电数字孪生系统首先要对各种复杂设备的参数、性能、特征、以及设备之间的关系的进行数字表达，满足静态以及动态的设备及系统建模需要(如设备约束关系、设备上下游关系、设备参数动态定义、设备参数关联性等)，以适应不同设备、系统、机组的定义需要 [1] 。
要实现现实核电站与数字孪生系统的交互与协同，就需要使用高效的数据分析和计算技术，譬如利用分布式并行计算技术、多线程技术等，将各类核电设备和系统的性能参数、运行状态、影响因素和运行流程实时反馈在仿真模型中，继而实现针对系统和设备运行过程的真实映射 [2] ，这就是仿真模型状态跟踪技术。
DCS初始化是仿真模型状态跟踪技术中的关键环节。DCS初始化有两种思路：一种是通过对过程控制层(DCS一层)的初始化实现，这种方法虽然精度和效率较高，但要求核电机组存储海量的历史数据，一般核电站出于安全考虑不会存储一层的I/O及中间变量，因此这种方法不能普及；另一种是通过对人机界面层(DCS二层)的初始化实现，这种方法不能保证初始化后的数字孪生样机与机组DCS完全一致，但是能够保证样机所有设备状态与机组一致，满足仿真模型状态跟踪的要求，因此，本文采用基于DCS二层的初始化方法。
2. 原理分析
2.1. 系统组成结构
DCS初始化工具的功能是将数字孪生系统仿真模型中的各个设备在较短时间内自动达到指定的电厂某个时刻实际运行状态。DCS初始化工具的输入是电厂某时刻的实际运行数据，输出是达成和电厂运行状态一致的设备仿真模型。(图1)
Figure 1. DCS initialization system composition structure diagram
图1. DCS初始化系统组成结构图
2.2. DCS初始化工具功能模块
根据电厂实际运行数据分析和模型建模经验，将DCS初始化工具分为创建跟踪点、初态分析以及初始化三大模块，各模块依次按顺序执行。下文介绍三大模块各自的主体功能设计。(图2)
Figure 2. DCS initialization tool module relationship diagram
图2. DCS初始化工具模块关系图
3. 实现方法
3.1. 创建跟踪点
创建跟踪点是指DCS初始化工具接收到电厂DCS某一时刻的运行数据后，对数据进行分析，初步判断当前的数据和已保存的仿真模型IC数据的参数近似性，并将数字孪生系统复位到与电厂该时刻最相近状态的过程。状态确认分为机组数据导入、数据分析和检索、和数字孪生系统状态复位3个步骤。
为驱动状态复位功能，需预存足够的模型IC (初始状态)。按5%功率间隔从满功率到零功率降低过程保存21个IC状态和零功率到40%全功率8个IC状态，共29个预存IC状态点。按功率P、温度T对应IC号N关系建立检索函数关系式：
f(P,T) = N
3.1.1. 机组数据导入
在核电数字孪生系统中，预先将机组数据存储在时序数据库IOTDB中，IOTDB是一款聚焦工业物联网、高性能轻量级的时序数据管理系统，也是一款开源时序数据库，为用户提供数据收集、存储和分析等服务 [3] 。
在DCS初始化工具的客户端界面选择初始化的时间点，工具查询时序数据库获取机组该时刻的电功率、平均温度以及所有设备的状态，并写入数字孪生系统的共享内存中，供后续数据分析和初始化使用。
3.1.2. 数据分析和检索
数据分析过程中首先获取共享内存中当前电功率数据，将电功率和额定功率进行比较计算得到归一化功率。如归一化功率大于40%，则从满功率向下寻找最接近该功率的IC；如归一化功率小于40%，则先依据功率和温度判断该过程是升功率或降功率，确定过程后再依据功率检索最接近IC。整个过程如下图3所示：
3.1.3. 数字孪生系统状态复位
通过数据分析和检索获取数字孪生系统中最接近的IC后，DCS初始化工具通过仿真平台的外部接口复位该IC，此时，数字孪生系统与机组达到较为相近的运行工况。
3.2. 初态分析
以福清核电3号机组为例，DCS二层采用的是ATOS的ADACS_N系统，该系统中设备状态(device state，简称DSTA)分为9类：全开(启/合闸)、全关(停/分闸)、正在开、正在关、中间位置、开故障、关故障、中间位置故障、其他故障；设备基本状态(basic device state，简称BSTA)分为4类：全开(启/合闸)、全关(停/分闸)、中间位置、故障。两者之间的关系如表1所示：
Table 1. The corresponding relationship between equipment status and equipment basic status
表1. 设备状态与设备基本状态对应关系表
ADACS_N系统数据库中存储的设备状态码为“基本设备状态码”，数字孪生系统的计算服务器也计算了设备BSTA，遍历数字孪生系统中所有设备，获取BSTA并与机组当前时刻的设备状态进行比较，即可得到待初始化设备清单，机组设备BSTA为仿真模型目标状态。
3.3. 对象初始化
初始化对象按照PO (即process object：过程对象类型)可分为AIR_OP_DAMPER、BINARY_ACTUATOR、BLOCK_VALVE、BLOCK_STP_VALVE、SOL_VALVE_FO、SOL_VALVE_FC、LOOP_CONTROL、LOOP_VALVE_FO、LOOP_VALVE_FC、SELECTOR、MULTI_COMMANDS共11类 [4] ，为了简化模型，将初始化对象抽象为如下3类(表2)：
Table 2. Initialization object classification table
表2. 初始化对象分类表
3.3.1. 普通开关类设备的初始化
对于普通开关类设备，初始化之前记住数字孪生系统各PO的手/自动状态，然后，发送CIN_5命令把所有PO统一切到手动状态，此时可以通过外部操作命令来控制设备。根据目标状态发送CIN_1或CIN_2命令打开或关闭设备，最后将样机初始状态为自动的设备切回自动状态。初始化流程如图4所示。
Figure 4. The flow chart of the initialization of ordinary switch equipment
图4. 普通开关类设备初始化流程图
3.3.2. 调节阀的初始化
对于调节阀，根据ADACS_N系统历史数据存储机制，分为有目标阀位的调节阀和只有开关状态的调节阀两类 [5] 。在对调节阀进行初始化之前记住样机各PO的开/闭环调节状态，发送CIN_8命令把所有PO统一切到开环调节状态。对于只有开关状态的调节阀，将阀位调到0或100；对于有目标阀位的调节阀，将阀位调到目标阀位，然后仿真模型通过目标阀位计算出闭环调节模式下PID的设定值输入，以免调节阀切回闭环调节模式时PID模块产生扰动。最后，将样机初始状态为闭环调节的设备切回闭环调节状态。初始化流程如图5所示。
Figure 5. Flow chart of regulating valve initialization
图5. 调节阀初始化流程图
3.3.3. 组控开关和选择器的初始化
对于组控开关(KG)和选择器(KC)的初始化，一般用于互为备用的设备的列间切换，不单独对KC、KG进行初始化。由于部分设备的开/关存在一些信号闭锁，需要通过KC、KG解除闭锁信号才能对设备进行操作，对于这种情况设计了配置文件处理此类包含复杂逻辑的PO。配置文件包含操作对象、操作类型、前置条件和后置条件，初始化程序首先按照前置条件操作相应对象；然后对操作对象执行预定义的操作；最后按照后置条件操作相应对象。配置文件样式如表3所示：
Table 3. Device initialization configuration table containing complex logic
表3. 包含复杂逻辑的设备初始化配置表
4. 实验与分析
4.1. 初始化对象状态分析
实验选取福清核电3号机“100% FP工况运行”为初始化目标，福清核电3号机组包含1345个设备，在初始化之前，其中1136个设备状态与数字孪生系统100% FP标准工况一致。在实验过程中，记录机组与数字孪生系统设备状态的一致性，结果如图6所示。
Figure 6. Device state consistency during initialization
图6. 初始化过程中设备状态一致性
4.2. 数字孪生系统状态分析
为了验证DCS初始化后，数字孪生系统的总体运行状况，选取了M310机组中与运行指标紧密关联的29个重要参数，记录经初始化后这些参数与机组100%FP工况下相关参数的偏差值和偏差比例，结果见表4。
Table 4. Important parameter deviation table
表4. 重要参数偏差表
5. 结论
经以上分析，使用本DCS初始化技术能够使初始化后的数字孪生样机设备状态与机组完全一致，DCS初始化后，数字孪生系统保持稳定，且重要运行参数的偏差在合理范围内。未来结合工艺系统、堆芯、热工、电气等子系统的初始化技术，能够最终实现核电数字孪生系统仿真模型状态跟踪。