翁品迪 陳 博 俞 立
1.浙江工業(yè)大學(xué)網(wǎng)絡(luò)空間安全研究院 杭州 310023 2.浙江工業(yè)大學(xué)信息工程學(xué)院 杭州 310023
信息物理系統(tǒng)(Cyber-physical systems,CPSs)是一個(gè)綜合了計(jì)算、網(wǎng)絡(luò)和物理環(huán)境的多維智能化復(fù)雜系統(tǒng),它借助有線或無(wú)線通信網(wǎng)絡(luò)將各個(gè)關(guān)鍵設(shè)施整合在一起,使得人機(jī)和物理進(jìn)程的交互更加便捷.隨著網(wǎng)絡(luò)通信、嵌入式系統(tǒng)、計(jì)算機(jī)控制及相關(guān)硬件技術(shù)的不斷發(fā)展,信息物理系統(tǒng)已引起了工業(yè)界的廣泛關(guān)注[1-3].然而在此框架下,原有系統(tǒng)的封閉性被打破,使其面臨著來(lái)自網(wǎng)絡(luò)攻擊的安全威脅.例如:2010年,伊朗的首座核電站-布什核電站被震網(wǎng)病毒攻擊,使得伊朗的第一座核設(shè)施推遲發(fā)電,嚴(yán)重?fù)p害了伊朗的工業(yè)設(shè)施[4];2011年美國(guó)伊利諾伊州一處水利控制系統(tǒng)遭到網(wǎng)絡(luò)攻擊,險(xiǎn)些造成大面積的供水中斷.當(dāng)信息物理系統(tǒng)遭受網(wǎng)絡(luò)攻擊時(shí),準(zhǔn)確、及時(shí)地監(jiān)測(cè)到攻擊信號(hào)對(duì)于監(jiān)測(cè)中心采取高效的防御策略是至關(guān)重要的[5].
從控制角度來(lái)看,信息物理系統(tǒng)是傳統(tǒng)數(shù)字控制系統(tǒng)融合通信技術(shù)的下一代網(wǎng)絡(luò)化控制系統(tǒng)[6].信息物理系統(tǒng)中被控對(duì)象的量測(cè)信號(hào)和控制信號(hào)均通過(guò)網(wǎng)絡(luò)進(jìn)行傳輸,因此無(wú)論是傳感器到控制器端,還是控制器到執(zhí)行器端,均有受到惡意網(wǎng)絡(luò)攻擊的可能,從而影響系統(tǒng)的穩(wěn)定性,甚至引發(fā)系統(tǒng)崩潰,造成嚴(yán)重的生產(chǎn)事故與經(jīng)濟(jì)損失[7-8].典型的網(wǎng)絡(luò)攻擊有三種,分別是拒絕服務(wù)攻擊(Denial of service,DoS)、欺騙攻擊和重放攻擊.針對(duì)拒絕服務(wù)攻擊,文獻(xiàn)[9]將DoS 攻擊建模為一類能量約束問(wèn)題,針對(duì)攻擊者與防御者設(shè)計(jì)遞歸分布式卡爾曼估計(jì)器進(jìn)行雙邊優(yōu)化;文獻(xiàn)[10]假定DoS 攻擊有界的情況下,通過(guò)構(gòu)建嵌套切換模型得到了CPS 在基于包控制方法下的穩(wěn)定性條件.針對(duì)重放攻擊,文獻(xiàn)[11]從防御者的角度提出一種帶補(bǔ)償策略的數(shù)學(xué)模型來(lái)描述重放攻擊和帶寬約束,并在線性最小方差意義下設(shè)計(jì)了遞歸分布式卡爾曼融合估計(jì)器.欺騙攻擊又稱為假數(shù)據(jù)注入(False data injection,FDI)攻擊,它通過(guò)向系統(tǒng)注入錯(cuò)誤的控制信號(hào)或測(cè)量信號(hào)影響信號(hào)數(shù)據(jù)的準(zhǔn)確性.注意到FDI 攻擊可以通過(guò)欺騙攻擊檢測(cè)機(jī)制來(lái)影響信息物理系統(tǒng),從而造成攻擊檢測(cè)器的漏報(bào)或虛警[12].文獻(xiàn)[13]中研究了電力系統(tǒng)中針對(duì)FDI 攻擊的狀態(tài)估計(jì)問(wèn)題,揭示了現(xiàn)有錯(cuò)誤測(cè)量檢測(cè)算法中存在的脆弱性.該研究表明,即使攻擊者的資源受限,依舊可以改變狀態(tài)估計(jì)的結(jié)果;而在假定攻擊者可以獲取系統(tǒng)參數(shù)及所有數(shù)據(jù)流的基礎(chǔ)上,此類攻擊則可以對(duì)系統(tǒng)產(chǎn)生一定影響的同時(shí)卻不被檢測(cè)到.目前信息物理系統(tǒng)中的FDI 攻擊主要采用異常信號(hào)檢測(cè)方法進(jìn)行檢測(cè)[14],如基于二元假設(shè)的貝葉斯檢測(cè)[15-18]、基于卡爾曼濾波器的χ2檢測(cè)[19-21]和加權(quán)最小平方檢測(cè)[22-24].后兩種方法均是利用觀測(cè)值構(gòu)造新息殘差,然后通過(guò)與一個(gè)給定的閾值比較來(lái)判決是否受到攻擊.攻擊檢測(cè)的閾值對(duì)檢測(cè)精度起著重要作用,然而閾值通常在檢測(cè)前根據(jù)已有經(jīng)驗(yàn)進(jìn)行選擇,這無(wú)疑會(huì)降低檢測(cè)精度.而且對(duì)于CPSs 而言,特別是電力、醫(yī)療系統(tǒng),僅僅檢測(cè)或識(shí)別到攻擊是不夠的,因?yàn)檫@些基礎(chǔ)設(shè)施系統(tǒng)無(wú)法快速對(duì)其關(guān)閉、重啟、恢復(fù),以對(duì)抗攻擊.為此,不同于文獻(xiàn)[15-24]中基于閾值的FDI 攻擊信號(hào)檢測(cè)方法,本文將對(duì)CPSs 中控制器與執(zhí)行器間的通信網(wǎng)絡(luò)遭受的FDI 攻擊信號(hào)進(jìn)行實(shí)時(shí)估計(jì),不僅可以檢測(cè)攻擊信號(hào)是否存在,而且還可以掌握攻擊信號(hào)的基本特征,從而使防御方采取簡(jiǎn)單有效的補(bǔ)救措施以降低系統(tǒng)性能的損失程度.
為了對(duì)FDI 攻擊信號(hào)進(jìn)行實(shí)時(shí)估計(jì),本文將攻擊信號(hào)建模為狀態(tài)動(dòng)態(tài)方程中的一個(gè)未知參數(shù),然后在每個(gè)傳感器端利用自適應(yīng)卡爾曼濾波設(shè)計(jì)FDI 攻擊信號(hào)的局部估計(jì).由于基于遞歸最小二乘的自適應(yīng)Kalman 濾波方法適應(yīng)于被估計(jì)的參數(shù)是時(shí)不變或者變化緩慢的,而惡意的攻擊信號(hào)往往是時(shí)變的.因此,本文充分利用多傳感器融合所提供的冗余信息,通過(guò)引入補(bǔ)償因子和設(shè)計(jì)分布式融合準(zhǔn)則,來(lái)提高時(shí)變攻擊信號(hào)的估計(jì)精度.注意到補(bǔ)償因子是通過(guò)影響局部互協(xié)方差的信息來(lái)提高估計(jì)性能,而局部互協(xié)方差矩陣存在于分布式融合結(jié)構(gòu)中.最后,通過(guò)兩個(gè)仿真驗(yàn)證所設(shè)計(jì)的融合算法對(duì)FDI 攻擊信號(hào)估計(jì)的有效性.特別地,當(dāng)攻擊信號(hào)是時(shí)變時(shí),仿真結(jié)果表明所引入的補(bǔ)償因子可以明顯提高攻擊信號(hào)的融合估計(jì)性能.
符號(hào)說(shuō)明.E 表示數(shù)學(xué)期望,diag{N1,···,Nn}表示由N1,···,Nn組成的塊對(duì)角矩陣,I表示單位矩陣,⊥表示正交.
















圖1 算例1:情況1 中攻擊信號(hào)和融合估計(jì)的軌跡Fig.1 Example 1:The trajectories of attack signal and its fusion estimation under Case 1

圖2 算例1:情況1 中攻擊信號(hào)的局部估計(jì)器與融合估計(jì)器之間的性能比較Fig.2 Example 1:The performance comparison between local estimators and fusion estimators under Case 1

圖3 算例1:情況2 中攻擊信號(hào)和融合估計(jì)軌跡Fig.3 Example 1:The trajectories of attack signal and its fusion estimation under Case 2

圖4 算例1:情況2 中攻擊信號(hào)的局部估計(jì)器與融合估計(jì)器之間的性能比較
Fig.4 Example 1:The performance comparison between local estimators and fusion estimators under Case 2



圖5 算例1:情況2 中不同補(bǔ)償因子下攻擊信號(hào)融合估計(jì)性能的比較Fig.5 Example 1:The comparison of fusion estimation performance of attack signal under different compensation factors under Case 2

圖6 算例1:情況2 中不同補(bǔ)償因子下攻擊信號(hào)融合估計(jì)的軌跡Fig.6 Example 1:The trajectories of fusion estimation of attack signal under different compensation factors under Case 2


圖7 算例2:不同補(bǔ)償因子下攻擊信號(hào)融合估計(jì)的軌跡Fig.7 Example 2:The trajectories of fusion estimation of attack signal under different compensation factors

圖8 算例2:不同補(bǔ)償因子下攻擊信號(hào)融合估計(jì)性能的比較Fig.8 Example 2:The comparison of fusion estimation performance of attack signal under different compensation factors
的情況(無(wú)補(bǔ)償?shù)那闆r).圖9 基于Monte Carlo 實(shí)驗(yàn)給出了局部估計(jì)器和融合估計(jì)器的MSE 曲線.圖9 表明融合性能明顯優(yōu)于局部估計(jì)器.仿真結(jié)果驗(yàn)證了算法的有效性.

圖9 算例2:攻擊信號(hào)的局部估計(jì)器與融合估計(jì)器之間的性能比較Fig.9 Example 2:The performance comparison of attack signal between local estimators and fusion estimators
本文研究了信息物理系統(tǒng)中遭受的FDI 攻擊信號(hào)的分布式融合估計(jì)問(wèn)題.首先,根據(jù)FDI 攻擊策略對(duì)系統(tǒng)的影響,將攻擊信號(hào)建模為系統(tǒng)狀態(tài)方程中的一個(gè)未知輸入,然后基于自適應(yīng)卡爾曼濾波,利用量測(cè)信息對(duì)系統(tǒng)狀態(tài)和攻擊參數(shù)進(jìn)行聯(lián)合估計(jì),從而得到攻擊參數(shù)的局部變化軌跡.在此基礎(chǔ)上,在分布式融合結(jié)構(gòu)下,引入補(bǔ)償因子并設(shè)計(jì)了最優(yōu)融合估計(jì)方法,提高了攻擊信號(hào)估計(jì)的精度.最后,通過(guò)例子驗(yàn)證了所提算法的有效性.今后的研究方向包括如下幾點(diǎn):首先,選擇合適的補(bǔ)償因子對(duì)本文算法估計(jì)精度的提高有著顯著的影響,如何選擇補(bǔ)償因子將是未來(lái)研究的重要方向之一;此外,許多實(shí)際系統(tǒng)無(wú)法簡(jiǎn)單地利用線性模型描述,因此將該算法推廣至非線性情況的攻擊信號(hào)估計(jì)具有重要的意義;最后,研究控制器與攻擊信號(hào)融合估計(jì)器及攻擊信號(hào)實(shí)時(shí)補(bǔ)償策略的協(xié)同設(shè)計(jì),以降低攻擊所帶來(lái)的損失是另一個(gè)重要的研究方向.