楊凱達
(中國人民解放軍95808部隊 酒泉 735006)
目標毀傷效果評估(Battle Damage Assessment,BDA)是目標保障工作的一項重要內(nèi)容,貫穿作戰(zhàn)全流程,是指揮員進行作戰(zhàn)決策、部隊檢驗訓(xùn)練效果的重要評估手段。目前,現(xiàn)有評估主要通過航空、航天偵察圖像,電子偵察及抵近偵察等手段作為評估依據(jù),在不同的戰(zhàn)場環(huán)境中,受限于時效性、安全性及敵方干擾和破壞等因素,制約目標毀傷效果評估工作的開展,難以實現(xiàn)高效準確的評估。隨著武器裝備的發(fā)展,電視制導(dǎo)導(dǎo)彈以其命中精度高、殺傷威力大和作戰(zhàn)效費比高等優(yōu)點,被廣泛應(yīng)用于實現(xiàn)對目標的“點穴”式毀傷。利用電視制導(dǎo)導(dǎo)彈攻擊目標時近實時的視頻圖像,進行目標毀傷效果評估,是一種全新的目標毀傷效果評估手段,便于飛行人員及武器操縱員在攻擊單目標情況下選擇目標易損傷部位,在攻擊目標群的情況下選擇體系中要害目標,進而提升實戰(zhàn)條件下目標選擇和戰(zhàn)場應(yīng)用能力。
基于導(dǎo)彈視頻圖像實現(xiàn)目標毀傷效果評估,相關(guān)研究鮮見于國內(nèi)。美軍從1999年起,就將“武器系統(tǒng)視頻”(Weapon System Video,WSV)作為進行物理毀傷評估所需的數(shù)據(jù)資料來源之一[1]。John提出了從多個層面改進目標毀傷效果評估的技術(shù)解決方案,其中包含加裝彈藥傳感器(Munition Sensors)[2]這一構(gòu)想。美國空軍實驗室(Air Force Research Lab)的研究人員在文獻[3]中,嘗試在精導(dǎo)彈藥上加裝傳感器以實現(xiàn)毀傷數(shù)據(jù)的收集,在安全的前提下,實現(xiàn)大范圍的毀傷信息的反饋。同時,他們還探索延長傳感器收集數(shù)據(jù)的時間,從而進行毀傷信息的收集、傳輸、轉(zhuǎn)發(fā)等研究。但到目前為止,還沒有哪種相關(guān)方法投入實際運用的報道。我方裝備了大量電視制導(dǎo)導(dǎo)彈,在制定作戰(zhàn)計劃和作戰(zhàn)實施階段都需要及時了解被打擊目標的毀傷情況,雖然有一些評估方法,但總體上比較簡單、科學(xué)性不強,需要進一步研究和優(yōu)化[4~6]。
目前,大部分BDA方法和系統(tǒng)都難以對具有大量時序信息的信息源進行信息的處理和運用,本文提出的技術(shù)路線如圖1所示,在實施首輪空地打擊時,利用導(dǎo)彈回傳的視頻圖像,確定導(dǎo)彈在被打擊目標上的落點,通過對導(dǎo)彈戰(zhàn)斗部威力和目標易損性的分析,利用動態(tài)貝葉斯網(wǎng)絡(luò)模型,加入毀傷時間流的演化過程,構(gòu)建目標毀傷效果評估及預(yù)測模型,實現(xiàn)對目標的毀傷效果評估以及預(yù)測目標毀傷趨勢。
圖1 技術(shù)路線
目前,作者已完成了兩方面的研究,一是基于導(dǎo)彈實時傳回的視頻圖像為依據(jù),改進圖像處理算法,實現(xiàn)導(dǎo)彈落點定位[7];二是基于戰(zhàn)斗部威力/目標易損性(V/L)方法,結(jié)合戰(zhàn)斗部威力、雷達目標模型,建立毀傷效果模型,實現(xiàn)了毀傷節(jié)點數(shù)據(jù)模擬[8]。
為了正確地處理、表述、轉(zhuǎn)換和評估目標毀傷情況,BDA模型不僅需要能夠系統(tǒng)地處理戰(zhàn)場的不確定性信息[9],同時,還應(yīng)該實現(xiàn)對目標毀傷效果進行時間流上的推理和評估[10]。BDA模型的構(gòu)建應(yīng)著重思考三個方面的問題[11]:一是評估信息的不確定性,戰(zhàn)時評估信息來源廣泛,包括各類參戰(zhàn)人員、武器平臺、目標整編信息等,其中,存在著大量具有時間變化和因果關(guān)系的信息;二是戰(zhàn)場決策的時效性,戰(zhàn)場決策的緊迫性要求指揮員要在極短的時間內(nèi)完成決策,決策者需要了解當(dāng)次打擊后目標的毀傷程度,以便對是否進行二次打擊進行決策;三是戰(zhàn)場態(tài)勢的動態(tài)變化性。戰(zhàn)場態(tài)勢的變化,對目標毀傷的要求也隨著時間的推移發(fā)生變化。
貝葉斯網(wǎng)絡(luò)(Bayesian Network,BN)[12]主要包括靜態(tài)貝葉斯網(wǎng)絡(luò)(Static Bayesian Network,SBN)和動態(tài)貝葉斯網(wǎng)絡(luò)(Dynamic Bayesian Network,DBN)兩種[13],是一種基于貝葉斯定理和貝葉斯假設(shè),描述各變量之間概率關(guān)系,進行不確定事件分析和推理的工具[14]。戰(zhàn)爭條件下,戰(zhàn)場環(huán)境不斷變化,利用SBN網(wǎng)絡(luò)模型無法處理動態(tài)系統(tǒng)中存在的問題,特征變量在不同時刻下會呈現(xiàn)出不同的狀態(tài)。評估實驗將毀傷時間流引入BN,構(gòu)建可以處理時間序列數(shù)據(jù)的動態(tài)貝葉斯網(wǎng)絡(luò)模型。
利用狀態(tài)空間模型可以實現(xiàn)在時間流上基于先驗知識的預(yù)測,構(gòu)建包括先驗概率P(Χ1) 、狀態(tài)轉(zhuǎn)換模型P(Χt|Χt-1)及觀測函數(shù)P(Yt|Χt) 的模型。設(shè)一個隨時間變化的變量 Zt,用 Ut、Xt、Yt表示狀態(tài)空間的三個部分,即輸入、隱含和輸出變量[15]。設(shè)一個D維實值觀測向量{Y1,Y1,…,YT},在每一個時間步Y(jié)t中,利用Χt進行建模,則:
模型構(gòu)成結(jié)構(gòu)為隱含狀態(tài)變量Χ={x0,x1,…,xT-1}序列的概率分布函數(shù)和觀測變量Y={y0,y1,…,yT-1}的序列,變量數(shù)為T,表示系統(tǒng)發(fā)生的時間界,即:
一般要求模型滿足如下兩個假設(shè):一是時不變假設(shè),即條件概率變化一直平穩(wěn),表示為P(Χt|Χt-1)與時間t無關(guān);二是馬爾科夫假設(shè),即P(Χt+1|Χ1,Χ2,…,Χt)=P(Χt+1|Χt)。
2時間片貝葉斯網(wǎng)是一個貝葉斯網(wǎng)絡(luò)片段,節(jié)點有Χ∪Χ'。其中,Χ的節(jié)點沒有雙親節(jié)點,Χ'的節(jié)點條件概率分布為P(Χ′|Pa(Χ′)),該片段表現(xiàn)的條件概率分布為
DBN由先驗網(wǎng)B1和轉(zhuǎn)移網(wǎng)B→兩部分構(gòu)成,B1為先驗分布P(Χ1)的BN,B→為一個2時間片貝葉斯網(wǎng),通過有向無環(huán)圖對P(Χt|Χt-1)定義為
給定窗口長度T,通過疊加(B1,B→) 形成一個完整的BN,擴展到計算DBN中每一個點的聯(lián)合概率分布。
3.2.1 拓撲結(jié)構(gòu)
BN模型拓撲結(jié)構(gòu)主要體現(xiàn)變量之間的相互關(guān)系,各個變量對于目標毀傷效果評估的決定權(quán)重是不同的,需要將所有因素融合處理,才能降低不確定性對評估結(jié)果的影響[16]。建立DBN模型是一個逐層推理的過程,先根據(jù)評估特征變量確定BN拓撲結(jié)構(gòu),再利用轉(zhuǎn)移BN模型,對先驗BN模型進行擴展,根據(jù)相鄰時刻特征變量之間的相互影響程度,建立因果關(guān)系,充分反映特征變量之間的概率依賴關(guān)系和毀傷時間流下的演化過程。將SBN網(wǎng)絡(luò)沿時間流展開,得到DBN目標毀傷效果評估模型,如圖2所示。
圖2 DBN評估模型結(jié)構(gòu)
3.2.2 評估變量
將{武器威力,武器落點,目標類型,目標作戰(zhàn)能力,目標防御能力,目標修復(fù)能力}作為目標毀傷效果評估的特征變量屬性集,其子變量的節(jié)點狀態(tài)集合如表1所示。
表1 模型各節(jié)點狀態(tài)集合
采用Netica軟件開展仿真,依據(jù)DBN模型,建立推理圖形。模型初始化后,加入初始數(shù)據(jù),通過仿真,得到的模擬導(dǎo)彈落點及V/L結(jié)果,更新模型節(jié)點數(shù)據(jù)。網(wǎng)絡(luò)推理被觸發(fā),各網(wǎng)絡(luò)節(jié)點的概率分布數(shù)據(jù)更新,獲得T時刻的評估結(jié)果。設(shè)定兩種典型作戰(zhàn)想定方案,如表2所示。
表2 兩種典型作戰(zhàn)想定
打擊目標為某雷達永備陣地,配屬維修分隊,維修實力較強,對遭受中度和輕度的打擊,具有一定的抗毀傷能力。相鄰時刻的條件轉(zhuǎn)移概率的取值與相鄰時刻的選取密不可分,間隔時刻不同,條件概率也就隨之變化。根據(jù)維修的成效性評估,設(shè)定條件轉(zhuǎn)移評估時刻的間隔為1天,評估時刻為時刻T1至?xí)r刻T10。模型時刻間的條件轉(zhuǎn)移概率表見表3。
表3 相鄰時刻的條件轉(zhuǎn)移概率表
初始數(shù)據(jù)輸入后,得到兩種想定條件下的評估模型推理結(jié)果(時刻T1),各節(jié)點狀態(tài)概率如圖3、圖4所示。
圖3 想定方案A模型推理結(jié)果(時刻T1)
圖4 想定方案B模型推理結(jié)果(時刻T1)
BDA與V/L(LE、VU)三個關(guān)鍵節(jié)點在各時刻的概率分布表如表4、表5所示。
表4 關(guān)鍵節(jié)點概率分布表(想定方案A)
表5 關(guān)鍵節(jié)點概率分布表(想定方案B)
不考慮時間流對目標毀傷效果評估的影響,只能實現(xiàn)基于輸入數(shù)據(jù)的靜態(tài)評估,初始數(shù)據(jù)輸入后,BDA與V/L(LE、VU)三個關(guān)鍵節(jié)點的概率分布表如表6所示。
表6 關(guān)鍵節(jié)點概率分布表
如圖5所示,想定方案A條件下,隨著時間流的演進,目標處于重度毀傷(T)、中度毀傷(M)、輕度毀傷(L)三個狀態(tài)的概率分布,在10個時刻內(nèi)變化較為明顯。目標處于T、M、L狀態(tài)的概率演化分別呈現(xiàn)遞減、平穩(wěn)、上升的趨勢。T1時刻的目標狀態(tài),表現(xiàn)為其遭受打擊后實時的毀傷情況,目標處于T狀態(tài)的概率較高,達56.7%。T1~T6時刻期間,隨著維修力量的介入,目標整體毀傷情況有所緩解。T6時刻目標處于T狀態(tài)的概率達50.1%,與處于M、L狀態(tài)的概率之和基本持平,之后目標毀傷情況明顯好轉(zhuǎn),其處于T狀態(tài)的下降趨勢加快,相對應(yīng)的,處于L狀態(tài)的趨勢明顯增強。T6~T7時刻之間,目標處于M、L狀態(tài)的概率分布交匯,T9~T10時刻之間,目標處于T、L狀態(tài)的概率分布交匯,目標毀傷狀態(tài)明顯好轉(zhuǎn)。通過分析可知,T6時刻是實施二次打擊的最佳時機。
圖5 想定方案A條件下目標毀傷狀態(tài)分布趨勢
如圖6所示,想定方案B條件較為有利,我方已部分取得制空權(quán),敵方無干擾,導(dǎo)彈制導(dǎo)精度較高,在T8時刻之前,目標處于T狀態(tài)的概率,一直處于80%以上,反映此次毀傷效果較為嚴重,取得了較好的打擊效果。整體上看,目標處于T狀態(tài)的概率經(jīng)歷了小幅度提升,在T4時刻達到峰值,達83.8%,說明有毀傷二次效應(yīng)出現(xiàn);在T7時刻后,概率下降趨勢開始加強,目標處于M狀態(tài)、L狀態(tài)的概率同步有所提升,但都在10%之內(nèi)。通過分析可知,目標在較長時間內(nèi)毀傷較為嚴重,功能恢復(fù)的可能性較小,沒有二次打擊的必要。
圖6 想定方案B條件下目標毀傷狀態(tài)分布趨勢
在同樣觀測數(shù)據(jù)下,兩種模型在T1時刻的毀傷評估等級基本一致,未加入時間流的BN模型只能評估T1時刻目標的毀傷情況,而加入時間流的DBN模型可以對T1~T10時刻的目標毀傷情況進行評估,實現(xiàn)對時刻間毀傷情況的利用與轉(zhuǎn)化,較好地展示了目標毀傷情況在時間流下的變化。
結(jié)合某次實彈演練數(shù)據(jù),對DBN模型的可行性進行驗證。演練主要對某型雷達目標進行打擊,打擊共分為兩個階段:第一階段載機在較遠距離發(fā)射導(dǎo)彈,對雷達目標進行打擊,打擊后由情報人員進行目標的毀傷效果判定;第二階段在第一階段的基礎(chǔ)上,對目標進行二次打擊,打擊后由情報人員進行目標的毀傷效果判定。導(dǎo)彈采用爆破戰(zhàn)斗部,由視頻監(jiān)視設(shè)備確定導(dǎo)彈在目標上的落點。第一階段打擊結(jié)束后,目標受損,經(jīng)專家綜合評估后,判定目標的毀傷為輕度毀傷。通過DBN模型的計算,得到第一階段打擊后雷達處于T、M和L狀態(tài)的概率分別為14%、21%和65%,目標處于L狀態(tài)的概率最大,模型計算結(jié)果與專家評估相一致。第二階段的打擊結(jié)束后,經(jīng)專家綜合評估后,判定雷達的毀傷為中度毀傷。DBN模型計算得到目標處于T、M和L狀態(tài)的毀傷概率分別為37%、48%和15%。目標處于M狀態(tài)的概率最大,模型計算結(jié)果與專家評估相一致,基本符合時間演化過程毀傷評估結(jié)果變化。綜上所述,利用DBN模型開展目標毀傷效果評估的準確性較高,與人工評估結(jié)果基本一致。
構(gòu)建基于DBN的目標毀傷效果評估模型,根據(jù)V/L變量分析結(jié)果,建立拓撲結(jié)構(gòu)并確定節(jié)點參數(shù),實現(xiàn)了典型作戰(zhàn)想定下的目標毀傷效果評估。通過與BN模型的實驗對比,DBN模型能夠?qū)崿F(xiàn)對時間演化過程的描述,評估目標實時毀傷效果及趨勢,模型評估結(jié)果準確性較高,有效克服了傳統(tǒng)評估方法的主觀性和不確定性。