自动化学报  2018, Vol. 44 Issue (8): 1405-1412   PDF    
基于粒子滤波的工业控制网络态势感知建模
陆耿虹, 冯冬芹     
浙江大学 智能系统与控制研究所 工业控制技术国家重点实验室 杭州 310027
摘要: 粒子滤波(Particle filtering,PF)算法能有效地对工控系统这一类非线性、非高斯噪声系统进行状态估计,但在实际采用经典粒子滤波状态估计检测攻击时,实验结果显示该方法存在很高的漏检率,无法保障系统安全.因此改进经典算法,提出了基于粒子滤波输入估计的态势理解算法.该算法在考虑系统输入与输出关系的同时,结合蒙特卡洛思想,提取工控系统态势特征,计算态势指标,最终实现态势理解.实验结果表明,该算法能有效地感知持续性攻击,并判断系统态势.
关键词: 工控系统     态势感知     粒子滤波     态势理解    
Modeling of Industrial Control Network Situation Awareness With Particle Filtering
LU Geng-Hong, FENG Dong-Qin     
State Key Laboratory of Industrial Control Technology, Institute of Cyber-Systems and Control, Zhejiang University, Hangzhou 310027
Manuscript received : December 19, 2016, accepted: May 22, 2017.
Foundation Item: Supported by National Natural Science Foundation of China (61433006)
Author brief: LU Geng-Hong Ph.D. candidate at the Institute of Cyber-Systems and Control, Zhejiang University. Her research interest covers industrial control system network security situation awareness
Corresponding author. FENG Dong-Qin Professor at the State Key Laboratory of Industrial Control Technology, Institute of Cyber-Systems and Control, Zhejiang University. His research interest covers field bus, real-time ethernet, industrial wireless communication technology, security of industrial control system, and network control system. Corresponding author of this paper
Recommended by Associate Editor GAO Hui-Jun
Abstract: Particle filtering (PF) algorithm can estimate the states of industrial control systems, which are non-linear and have non-Gaussian noises. However, when using classical particle filtering state estimation to detect continuous attacks, it is shown that the false negative rate is too high to ensure the security of the system. Therefore, a situation perception algorithm by means of particle filtering input estimation is proposed to improve the effectiveness of the classical algorithm. Considering the relationship between system input and output and combining Monte-Carlo simulation, the proposed algorithm can extract industrial control system situation features, calculate situation metrics and realize the situation perception. Experimental results indicate that the proposed algorithm can recognize continuous attacks and judge the system situation effectively.
Key words: Industrial control system     situation awareness     particle filtering (PF)     situation perception    

工控网络处于快速发展阶段, 由于工业通信协议中存在不可避免的漏洞, 工控网络容易遭受攻击者的恶意攻击, 给工控网络的安全带来巨大威胁[1], 例如伊朗的震网病毒事件[2], 就是攻击者借助Stuxnet病毒对可编程逻辑控制器(Programmable logic controller, PLC)代码进行篡改, 实现攻击, 从而达到破坏离心机正常运行的恶意目的, 并造成不可逆转的严重事故.此外, 由于工控网络自身存在的复杂性增加了攻击判断与检测的难度, 尤其当管理员处于高压紧张的环境下, 更容易发生判断失误[3].

因此, 如何提高对工控系统网络整体态势的准确判断成为当务之急. 1999年, Bass[4]首次将态势感知(Situation awareness, SA)与网络安全技术相结合, 以期能准确全面地掌握系统的安全状态, 预防事故发生, 为安全管理员提供可靠有效的决策依据. Naderpour等[5]提出新型异常态势模型(Abnormal situation modeling, ASM), 构建贝叶斯网络对多种态势进行分析, 该方法针对安全性要求极高的系统, 利用风险指标判断系统在出现异常态势(Abnormal situation)情况下的危险态势(Hazardous situation)风险等级, 以此确定系统态势. Kim等[6]提出一个基于贝叶斯推论的核电厂态势评估解析模型, 该方法是对核电厂操作员在面对事故发生时的态势评估思考模型(Mental model)建模.

以上两种方法, 均是在系统出现异常或事故的情况下对系统态势进行分析, 其前提是异常或事故报警为可信的; 但是, 若系统遭受到欺骗性攻击(例如假数据注入攻击[7]), 由于系统内的报警系统被攻击者蒙蔽, 此类态势感知技术将无法对系统真实态势进行感知.

针对以上问题, 结合系统在遭受到攻击的情况下系统内部状态值会发生相应改变, 本文提出基于粒子滤波的工业控制网络态势感知建模方法.粒子滤波(Particle filtering, PF)[8-9]是采用蒙特卡罗仿真完成递推贝叶斯滤波过程, 核心是采用一组粒子近似表示系统的后验概率分布, 然后使用近似表示估计非线性系统的状态[10]. Arulampalam等[11]对PF的非线性/非高斯噪声系统状态的估计性能进行了考察, 并将其与扩展卡尔曼滤波进行比较, 证明PF能有效解决非线性/非高斯噪声系统的状态估计问题.

由上述文献的研究结果可知, 传统的PF状态估计, 可以预测非线性非高斯系统的态势变化.但是在实际应用PF状态估计算法对工控系统受到的持续攻击进行检测时, 实验结果显示漏检率高达96%, 这是因为在控制器的作用下, 系统从临界状态趋于稳态, 此时, 由于PF状态估计算法自身的跟踪能力可跟踪到系统的攻击状态导致无法检测出异常, 从而使得该算法在检测持续性攻击时, 存在较高的漏检率, 这一不足也为工控系统安全带来威胁.

因此, 为了实现对持续性攻击的检测, 降低PF状态估计算法对该类攻击的漏检率, 本文将工控系统内在特征与攻击特征相结合, 提出基于PF输入估计的态势理解算法.该方法考虑了工控系统遭受到网络攻击时(以传感器参数篡改为例), 虚假的传感器参数导致系统的输入值与输出值之间的非线性关系发生变化, 利用这一特点, 结合Monte-Carlo思想, 对系统输入的先验概率进行随机取样, 依据相似性对样本粒子分配权值, 并获取系统输入的估计值将估计值与实际值之间的差值, 作为系统态势特征, 判断工控系统是否处于危险态势, 避免了由于PF状态估计引起的漏检, 为针对工控系统网络的持续性攻击检测提出新思路.在本文最后, 对经典PF状态估计算法和态势理解算法进行仿真验证.需要说明, 不同于文献[12]对网络安全态势预测算法的精度进行考察(算法精度越高, 对复杂网络环境的预测结果越准确), 本文考虑到工控系统遭受攻击的后果严重性, 系统中出现的漏报与误报将会误导操作员对工控系统的态势判断, 从而引发灾难, 因此参考Salerno等[13]提出的态势指标, 对本文提出算法报警结果的漏检率与错误率进行分析, 而不再考虑算法精度.实验结果表明, 本文提出的算法能有效感知系统中的不同态势(正常态势及危险态势), 漏检率与错误率均处于较低水平.

1 工控网络态势感知模型 1.1 网络安全态势感知模型简介

Endsley[14]在1988年提出了态势感知的定义:在一定的时空条件下, 对环境因素的获取、理解以及对未来状态的预测.

网络安全态势感知(Network security situation awareness, NSSA)模型分为三个层次, 通过对来自系统的数据进行处理后, 获取系统当前态势, 并对未来态势进行预测.从下至上依次为(如图 1所示):

图 1 态势感知模型 Figure 1 Situation awareness model

1) 态势要素获取:态势要素获取层是NSSA模型的基础, 主要包括对数据的预处理和特征提取.其目的主要在于对工控系统中的海量数据进行缩减, 保留关键信息, 并从中提取特征.

2) 态势理解:对获取的特征进行进一步处理, 包括数据关联、特征检测、识别与分类, 并对多个分类结果进行决策融合, 获取最终决策即整体系统态势.

3) 态势预测:利用预测算法对工控网络态势的趋势进行预测.

1.2 工控网络态势

出于保护工控系统安全的目的, 将工控网络态势分为安全态势和危险态势[15].安全态势指系统中的过程参数均处于系统既定的安全值范围内; 危险态势指系统遭受到攻击, 过程参数超出安全临界值$\delta$的工控系统状态.

一般的PLC系统如图 2所示, 传感器在将传感数据发送至PLC的过程中, 可能会遭受到攻击者的攻击, 真实传感数据被篡改[16], 导致控制系统不稳定.本文假定系统的危险态势是由攻击者施加的假数据注入攻击[17]产生.

$ \begin{align} \widetilde{y}(t)=y(t)+\alpha(t), \quad t\in T_{\rm atc} \end{align} $ (1)
图 2 PLC系统示意图 Figure 2 A diagram of PLC implementation

其中, $y(t)$$\widetilde{y}(t)$分别是在安全态势和危险态势下, PLC接收到的数值. $\alpha(t)$为攻击者施加的攻击参数, $T_{\rm atc}=[t_{{\rm atc}_s}, t_{{\rm atc}_e}]$为从攻击开始$t_{{\rm atc}_s}$到攻击结束$t_{{\rm atc}_e}$的攻击时间.

1.3 工控网络态势感知建模

态势理解层作为态势感知过程中承下启上的中心环节, 利用态势要素对获取的数据进行态势理解, 并将数据处理结果输送至态势预测层.态势理解的优劣, 将直接影响态势感知的结果以及态势预测的性能.因此高效准确的态势理解过程在态势感知中, 显得尤为重要.

本文创新性地提出了工控网络态势感知模型, 如图 3所示, 该模型从下至上共包含三个部分:态势要素获取层、态势理解层以及后续态势感知过程, 本文主要对底部两层进行研究与分析.

图 3 工控网络态势感知模型 Figure 3 Industrial control network situation\\ awareness model

1) 态势要素获取层:采集来自各传感器节点的数据, 将其存储在数据集$D$中.

2) 态势理解层:利用PF算法估计系统输出, 并提取特征, 使用态势指标对要素理解的效果进行衡量, 最后将态势理解结果输出至态势评估与预测中.

定义1. 数据集合$D^{i}_{0:k}=\{d_{j}, j=0, \cdots, k\}$表示从$0\sim k$时刻, 第$i$个节点采集的数据集合.

定义2. 特征提取$E$作为PF输入估计的结果之一, 记录PF输入估计算法的输出结果, 是数据集$D^{i}_{0:k}$在经PF输入估计算法处理后, 得到的特征序列.

$ \begin{align} E=\{e^{i}_{\rm atc}\mid i=1, \cdots, m\} \end{align} $ (2)

定义3.  态势指标$M=\{MR, FR\}$包含两个元素.表示在进行PF状态估计时出现的错误率和漏检率.

$ MR=\frac{\sum\limits_{j=1}^{k}I(m_{j})}{\sum\limits_{j=1}^{k}S(m_{j})} \\[2mm] $ (3)
$ FR=1-\frac{\sum\limits_{j=1}^{k}C(m_{j})}{\sum\limits_{j=1}^{k}W(m_{j})} $ (4)

其中, $I$为误测特征个数, $S$为测得特征个数, $C$为正确测得的危险态势特征个数, $W$为期望得到的危险态势特征个数.

由于在特征提取过程中, 可能存在感知错误、遗漏等情况, 因此提出对态势特征获取过程中PF出现的错误特征数量和丢失的特征数量进行衡量.

本文提出的态势感知模型中的态势理解层不仅能对来自工控系统众多节点的海量数据进行简化, 给出系统态势特征, 而且可以利用态势指标对态势理解过程的质量进行计算, 为后续阶段的态势感知结果可靠性提供依据.

1.3.1 粒子滤波状态估计

为了获取非线性、非高斯噪声的工控系统状态特征, 采用PF算法进行状态估计, 是一种很有效的非线性滤波技术[18], 适用于任何能用状态空间模型表示的非高斯背景的非线性随机系统, 精度可以逼近最优估计.该算法的实质是由粒子及其权重组成的离散随机测度近似相关的概率分布, 并根据算法递推更新离散随机测度[19].

对于一个非线性、非高斯过程建模如下:

$ x_{k}=g(x_{k-1}, w_{k}) $ (5)
$ y_{k}=h(x_{k}, v_{k}) $ (6)

其中, $x_{k}$为待估计的状态量, 假设状态$x_{k}$的先验分布$p(x_{0})$已知; $y_{k}$为已知的观测量; $g(\cdot)$为状态转移方程, $h(\cdot)$为观测方程; $w_{k}$$v_{k}$都为独立的噪声, 分别称为状态噪声和观测噪声.

步骤1. 初始化.设$k=0$, 采样: $x^{i}_{0}\sim p(x_{0})$, 根据$p(x_{0})$分布采样得到$N$$x^{i}_{0}$, $i=1, 2, \cdots, N$.

步骤2. 重要性权值计算.采样$x^{i}_{k} \sim q(x_{k}|$ $x^{i}_{0:k-1}, y_{0:k})$, $i=1, 2, \cdots, N$, 利用下式计算重要性权值.

$ \begin{align} w^{i}_{k}=&\ w^{i}_{k-1}p(y_{k}|x^{i}_{k-1})=\notag\\[2mm] &\ w^{i}_{k-1}\frac{p(y_{k}|x^{i}_{k})p(x_{k}|x^{i}_{k-1})}{ q(x_{k}|x^{i}_{0:k-1}, y_{k})} \end{align} $ (7)

其中, $p(y_{k}|x^{i}_{k})$为似然函数概率密度.

$ \begin{align} p(y_{k}|x^{i}_{k})=p_{e_{k}}(y_{k}-h(x^{i}_{k})) \end{align} $ (8)

式(7)中的$q(x_{k}|x^{i}_{0:k-1}, y_{k})$是重要性概率密度, 在本算法中, 选取使权系数方差最小的最优重要密度函数.

$ \begin{align} q(x_{k}|x^{i}_{0:k-1}, y_{k})=p(x_{k}|x^{i}_{k-1}) \end{align} $ (9)

在计算重要性权值后, 进行归一化.

$ \begin{align} \bar{w}^{i}_{k}=\dfrac{w^{i}_{k}}{\sum\limits^{N}_{i=1}w^{i}_{k}} \end{align} $ (10)

步骤3. 重采样.设定有效样本数$N_{\rm threshold}$, 并计算退化因子$N_{\rm eff}$

$ \begin{align} N_{\rm eff}=\dfrac{1}{\sum\limits^{N}_{i=1}(\bar{w}^{i}_{k})^{2}} \end{align} $ (11)

$N_{\rm eff}$越小, 意味着退化现象越严重.

$N_{\rm eff} < N_{\rm threshold}$, 则进行重采样, 将原来的带权样本$\{{x}^{i}_{0:k}, \bar{w}^i_k\}^N_{i=1}$映射为等权样本$\{\tilde{x}^{i}_{0:k}$, ${1}/{N}\}^{N}_{i=1}$.其中, $\tilde{x}^{i}_{0:k}$为重采样后的粒子.

步骤4. 预测.获取状态预测$\hat{x}_{k}$, 据此计算观测值估计$\hat{y}_{k}$, 并计算估计值与真实值的偏差$e_{k}$(特征).

$ \hat{x}_{k}=\sum\limits^{N}_{i=1}\bar{w}^{i}_{k}\tilde{x}^{i}_{k} $ (12)
$ \hat{y}_{k}=h(\hat{x}_{k}, v_{k}) $ (13)
$ e_{k}=y_{k}-\hat{y}_{k} $ (14)

步骤5. 输出.将$e_{k}$与安全阈值$\delta$进行比较, 当$|e_{k}|> \delta$时, 记录系统第$i$次出现危险态势的时刻$k^{i}_{\rm atc}$以及对应偏差$e^{i}_{k}$.

1.3.2 基于粒子滤波输入估计的态势理解算法

在实际应用PF状态估计对系统受到的攻击进行检测时, 实验结果显示, 对于系统中施加持续时间为100小时的攻击, PF状态估计算法可以在系统被攻击(即状态值发生突然变化)后的10分钟内检测到攻击, 即偏差值$e$将会超出阈值; 但是由于PF估计算法自身具备的状态变化跟踪能力, 在系统遭受到攻击的2小时后, 系统状态将趋于稳定, 即此时的偏差值$e$将会接近于零, 从而导致PF状态估计算法无法持续检测出系统受到的攻击.

考虑到对于输入输出之间存在函数关系$y=f(u)$的系统, 在遭受到攻击时, 系统的输入$u$也会发生变化, 因此, 本文提出基于PF输入估计的态势感知算法, 以实现对长时间持续攻击的态势感知, 降低系统漏报率.

基于PF输入估计的态势感知算法, 不再考虑状态的估计, 利用输出残差进行报警; 而是利用Monte-Carlo思想, 在已获得的输出$y_{k}$基础上, 估计输入$\hat{u}$.利用实际输入(控制器输出) $u$$\hat{u}$的差值, 对系统态势进行感知.

算法1.基于PF输入估计的态势感知算法

步骤1.  获取$t=0, \cdots, k$时刻的输出$y_{k}$;

步骤2. 初始化.设$k=0$, 在$[u_{\min}, u_{\max}]$内随机采样$N$次, 获得$N$个随机样本(粒子), 构成序列$\{u^{*}_{k}(i)|i=1, \cdots, N\}$;

步骤3. 权值分配.对每个粒子$u^{*}_{k}(i)$分配相应的权重$q^{*}_{k}(i)$, $i=1, \cdots, N$, 计算方式如下:

$ \begin{align} q^{*}_{k}(i)=\frac{p(y_{k}|u^{*}_{k}(i))}{\sum^{N}\limits_{j=1}p(y_{k}|u^{*}_{k}(j))} \end{align} $ (15)

其中, $p(y_{k}|u^{*}_{k}(i))$为相似性, 即

$ \begin{align} p(y_{k}|u^{*}_{k}(i))=p_{e^{*}_{k}}(y_{k}-f(u^{*}_{k}(i))) \end{align} $ (16)

步骤4. 重采样.采用多项式重采样法, 将原来的带权样本$\{{u}^{*}_{0:k}(i), {q}^{*}_{k}(i)\}^N_{i=1}$映射为等权样本$\{{\tilde{u}}^*_{0:k}, {1}/{N}\}^N_{i=1}$.其中, $\tilde{u}^{*}_{0:k}(i)$为重采样后的粒子, 获得新粒子样本.

步骤5. 输出估计.获取输出估计值预测$\hat{u}_{k}$, 并计算估计值与真实值的偏差$e^{*}_{k}$ (特征).

$ \hat{u}_{k}=\sum\limits^{N}_{i=1}q^{*}_{k}(i)\tilde{u}^{*}_{k}(i) $ (17)
$ e^{*}_{k}=u_{k}-\hat{u}_{k} $ (18)

步骤6. 计算态势特征.将$e^{*}_{k}$与安全阈值$\delta$进行比较, 当$|e^{*}_{k}|>\delta$时, 记录系统第$j$次出现危险态势的时刻$k^{*}_{\rm atc}(j)$以及对应偏差$e^{*}_{k}(j)$.

步骤7. 计算态势指标.利用式(3)和式(4)计算该系统的态势指标, 态势指标$M=\{MR, FR\}$.

步骤8. 态势理解结果输出.依据步骤6和步骤7的计算结果, 整理得到态势理解结果并作为算法输出$E$, $M$.

2 实验验证 2.1 仿真对象

某精馏塔提馏段温度单回路控制方案[20]图 4所示, 蒸馏塔提馏段某块板的温度为主变量, 控制器TC 21通过控制信号$u$控制蒸汽控制阀对温度进行控制, 温度传感器TT 21能对提馏段的温度$y$进行检测(系统输出), $y_{sp}$为设定值.

2.2 仿真过程

对PF状态估计算法和态势理解算法进行仿真, 通过计算两种算法的漏报率, 对算法有效性进行分析.

阶段1.  仿真模型建立.依据图 4的控制方案, 建立相应的控制系统方框图, 如图 5所示.对图 5中各环节进行参数设置: $G_{Tm}$为温度测量环节, $G_{Tm}$ $=$ ${1}/(s+1)$; 控制阀$G_{v}$为近似线性阀, $G_{v}=1$; 蒸汽流量对象$G_{p2}={0.1}/{(1.5s+1)}$; 提馏段温度对象的控制通道与扰动通道动态特性的参数设置分别为$G_{p1}$ $=$ ${5}/{(4s^{2}+5s+1)}$, $G_{d}={-0.5}/{(3s+1)}$; 单回路控制器TC的PID参数$K_{c}=2.4$, $T_i=8.8$, $T_{d}$ $=$ $2.2$.

图 4 提馏段温度单回路控制方案 Figure 4 Temperature single loop control scheme of distillation
图 5 提馏段温度单回路控制系统方框图 Figure 5 Block diagram of temperature single loop control scheme of distillation

阶段2. 态势设置.对安全态势及遭受不同时长攻击的两种危险态势进行仿真, 设定系统运行总时长$T_{\rm run}=500$.

1) 安全态势:系统正常运行, 温度设定值$y_{sp}=20$;

2) 危险态势1:在$k_{1}=200$$k_{2}=300$时刻, 攻击者篡改温度传感器值, 分别施加攻击强度为$\alpha(k_{1})$ $=30$, $\alpha(k_{2})=-30$的攻击, 每次攻击持续时长为$T_{\rm atc}=5$;

3) 危险态势2:在$k=200$时, 攻击者开始篡改温度传感器值, 攻击强度为: $\alpha(k)=30$, 攻击持续时长为$T_{\rm atc}=100$.

图 6为系统中可能出现的三种态势仿真结果示意图.

图 6 三种不同态势情况 Figure 6 The three different situations

阶段3.  算法验证.算法验证部分选择$t=70$以后的数据进行分析, 即对系统处于稳定后的数据进行仿真处理.

1) PF状态估计. 图 7为利用PF状态估计对三种情况进行态势检测的结果, 假设安全阈值$\delta=2$.从图 7可以看出, a) $M_{a}=\{0, 0\}$, 此时算法没有检测到攻击, 判断系统处于正常运行状态; b) $M_{b}=$ $\{0.0200, 0.1667\}$, 此时算法检测到两次攻击, 存在较低的错误率与漏检率; c) $M_{c}=\{0.2040, 0.9604\}$, 由态势指标可得, 该算法在面对持续时间较长的攻击时, 漏检率高达0.9604.但是能检测到攻击开始以及结束的时刻.

图 7 PF状态估计算法仿真结果 Figure 7 The simulation results related to PF state estimation algorithm

2) 基于PF输入估计的态势理解. 图 8为利用态势理解算法对三种情况进行态势检测的结果, 假设安全阈值$\delta=5$.从图 8可以看出, a) $M_{a}=$ $\{0, $ $0\}$, 此时算法没有检测到攻击, 判断系统处于正常运行状态; b) $M_{b}=\{0.0820, 0.1667\}$, 此时算法检测到两次攻击, 错误率与漏检率均较低; c) $M_{c}=$ $\{0.048, $ $0.0396\}$; 检测到长持续时间下的攻击, 错误率与漏检率均低于$5$%.

图 8 态势理解算法仿真结果 Figure 8 The simulation results related to situation awareness algorithm
2.3 仿真结果分析

对比两种算法的仿真结果可知:

1) PF状态估计算法

由a)和b)的检测结果可知, 当系统存在突然发生的变化时, 该算法能进行有效的检测, 特征值变化明显, 检测的错误率与漏检率较低; 由c)的仿真结果可知, 该算法无法感知到系统中出现的长持续性攻击, 存在很高的漏报率, 对系统安全产生严重威胁.

2) 基于PF输入估计的态势理解算法

a) 和b)的仿真结果证明, 该算法能有效跟踪和预测输入的变化趋势, 两种算法的漏检率相同, 尽管态势理解算法的错误率相对PF状态估计算法偏高$6$%, 但是对于系统安全性能而言, 该算法依然有效检测出了b)的危险态势; 由c)的仿真结果可得, 态势理解算法能有效检测到系统中存在的长时间持续的攻击, 错误率与漏检率均小于$5$%, 远低于PF状态估计算法.

对比仿真结果可知, 本文提出的基于PF输入估计的态势理解算法能有效判断系统中出现的危险态势, 为工控网络安全态势感知提供可靠的感知判据.

3 结论

本文介绍了基于粒子滤波的工业控制网络态势感知建模, 将文中构造的态势感知模型中的态势理解层作为研究重点, 提出了基于PF输入估计的态势理解算法.该算法是对经典PF状态估计的改进, 在Monte-Carlo思想的基础上, 考虑到系统输入与输出之间的联系, 通过输出值校准输入的估计值, 利用实际值与估计值之间存在的差值, 判断系统是处于安全态势还是危险态势.该算法具有较好的输入值估计能力, 可准确获取系统态势特征, 给出态势指标, 对特征的可靠性进行定量衡量.该算法可为后续的态势预测过程提供科学可靠的数据信息, 提升态势感知结果的准确性.

实验结果表明, 本文提出的态势理解算法, 能有效地利用数据源中的信息对系统内存在的危险态势进行感知.仿真结果表明该算法具有较高的可靠性与准确性.

参考文献
1
Huang Jia-Hui, Feng Dong-Qin, Wang Hong-Jian. A method for quantifying vulnerability of industrial control system based on attack graph. Acta Automatica Sinica, 2016, 42(5): 792-798.
( 黄家辉, 冯冬芹, 王虹鉴. 基于攻击图的工控系统脆弱性量化方法. 自动化学报, 2016, 42(5): 792-798.)
2
Genge B, Nai Fovino I, Siaterlis C, Masera M. Analyzing cyber-physical attacks on networked industrial control systems. Berlin Heidelberg, Germany: Springer, 2011, 167-183.
3
Lu J, Yang X W, Zhang G Q. Support vector machine-based multi-source multi-attribute information integration for situation assessment. Expert Systems with Applications, 2008, 34(2): 1333-1340.
4
Bass T. Multisensor data fusion for next generation distributed intrusion detection systems. In: Proceedings of the 1999 IRIS National Symposium on Sensor and Data Fusion. Washington, USA: IRIS, 1999. 24-27
5
Naderpour M, Lu J, Zhang G Q. An abnormal situation modeling method to assist operators in safety-critical systems. Reliability Engineering and System Safety, 2015, 133: 33-47. DOI:10.1016/j.ress.2014.08.003
6
Kim M C, Seong P H. An analytic model for situation assessment of nuclear power plant operators based on Bayesian inference. Reliability Engineering and System Safety, 2006, 91(3): 270-282. DOI:10.1016/j.ress.2005.01.012
7
Jia Chi-Qian, Feng Dong-Qin. Industrial control system devices security assessment with multi-objective decision. Acta Automatica Sinica, 2016, 42(5): 706-714.
( 贾驰千, 冯冬芹. 基于多目标决策的工控系统设备安全评估方法研究. 自动化学报, 2016, 42(5): 706-714.)
8
Doucet A, de Freitas N, Gordon N. Sequential Monte Carlo Methods in Practice. New York, USA: Springer, 2001.
9
Carpenter J, Clifford P, Fearnhead P. Improved particle filter for nonlinear problems. IEE Proceedings-Radar, Sonar and Navigation, 1999, 146(1): 2-7. DOI:10.1049/ip-rsn:19990255
10
Kadirkamanathan V, Li P, Jaward M H, Fabri S G. Particle filtering-based fault detection in non-linear stochastic systems. International Journal of Systems Science, 2002, 33(4): 259-265. DOI:10.1080/00207720110102566
11
Arulampalam S, Maskell S, Gordon N, Clapp T. A tutorial on particle filters for online nonlinear/non-Gaussian Bayesian tracking. IEEE Transactions on Signal Processing, 2002, 50(2): 174-188. DOI:10.1109/78.978374
12
Tang Yong-Li, Li Wei-Jie, Yu Jin-Xia, Yan Xi-Xi. Research on a prediction method of network security situation based on particle filter. Computer Applications and Software, 2017, 34(1): 293-297.
( 汤永利, 李伟杰, 于金霞, 闫玺玺. 基于粒子滤波的网络安全态势预测方法研究. 计算机应用与软件, 2017, 34(1): 293-297. DOI:10.3969/j.issn.1000-386x.2017.01.053)
13
Salerno J J, Blasch E P, Hinman M, Boulware D M. Evaluating algorithmic techniques in supporting situation awareness. In: Proceedings of the 2005 Multisensor, Multisource Information Fusion: Architectures, Algorithms, and Applications. Orlando, Florida, USA: SPIE, 2005. 96-104
14
Endsley M R. Design and evaluation for situation awareness enhancement. In: Proceedings of the 32nd Human Factors Society Annual Meeting. Santa Monica, USA: SAGE, 1988. 97-101
15
Naderpour M, Lu J, Zhang G Q. A situation risk awareness approach for process systems safety. Safety Science, 2014, 64: 173-189. DOI:10.1016/j.ssci.2013.12.005
16
Gonzalez C A, Hinton A. Detecting malicious software execution in programmable logic controllers using power fingerprinting. Critical Infrastructure Protection Ⅷ. Berlin Heidelberg, Germany: Springer, 2014, 15-27.
17
Lu Geng-Hong, Feng Dong-Qin. Industrial control system network security situation awareness modeling and algorithm implementation. Control Theory and Applications, 2016, 33(8): 1054-1060.
( 陆耿虹, 冯冬芹. 工控网络安全态势感知算法实现. 控制理论与应用, 2016, 33(8): 1054-1060.)
18
Cheng Q, Varshney P K, Michels J, Belcastro C M. Distributed fault detection via particle filtering and decision fusion. In: Proceedings of the 8th International Conference on Information Fusion. Philadelphia, PA, USA: IEEE, 2005. 1239-1246
19
Li Tian-Cheng, Fan Hong-Qi, Sun Shu-Dong. Particle filtering:theory, approach, and application for multitarget tracking. Acta Automatica Sinica, 2015, 41(12): 1981-2002.
( 李天成, 范红旗, 孙树栋. 粒子滤波理论、方法及其在多目标跟踪中的应用. 自动化学报, 2015, 41(12): 1981-2002.)
20
Dai Lian-Kui, Yu Ling, Tian Xue-Min, Wang Shu-Qing. Process Control Engineering. 3rd edition. Beijing: Chemical Industry Press, 2012.
( 戴连奎, 于玲, 田学民, 王树青. 过程控制工程. 第3版. 北京: 化学工业出版社, 2012.)