谷志鳴，高文明，魏瀟龍，姚登凱

（1.解放軍91899部隊(duì)，遼寧葫蘆島125001；2.空軍工程大學(xué)空管領(lǐng)航學(xué)院，西安710051）

基于蒙特卡洛法的無(wú)人機(jī)飛行沖突解脫安全評(píng)估*

谷志鳴1，高文明2，魏瀟龍2，姚登凱2

（1.解放軍91899部隊(duì)，遼寧葫蘆島125001；2.空軍工程大學(xué)空管領(lǐng)航學(xué)院，西安710051）

沖突解脫技術(shù)是無(wú)人機(jī)與有人機(jī)共域飛行的前提保障。針對(duì)目前缺少對(duì)沖突解脫技術(shù)安全評(píng)估的方法，提出運(yùn)用蒙特卡洛法對(duì)飛行沖突進(jìn)行仿真模擬，得出沖突解脫算法的安全效率。提出基于改進(jìn)蟻群算法的沖突解脫技術(shù)，實(shí)現(xiàn)無(wú)人機(jī)的飛行沖突解脫；利用蒙特卡洛法，對(duì)飛行沖突過(guò)程進(jìn)行大規(guī)模仿真，計(jì)算沖突解脫效率，與基本蟻群算法進(jìn)行對(duì)比，并從安全角度對(duì)結(jié)果進(jìn)行分析。提出的思路可為未來(lái)無(wú)人機(jī)空域運(yùn)行安全評(píng)估體系的構(gòu)建提供一種理論參考。

無(wú)人機(jī)，蟻群算法，蒙特卡洛，沖突解脫

0 引言

無(wú)人駕駛航空器即無(wú)人機(jī)，因其體積小、使用風(fēng)險(xiǎn)小，成本低、對(duì)環(huán)境要求低等諸多優(yōu)點(diǎn)，在軍事和民用領(lǐng)域受到青睞［1］。隨著無(wú)人機(jī)技術(shù)的發(fā)展，無(wú)人機(jī)早已具備了全空域飛行的能力，但由于“感知-避讓”功能的不完善，國(guó)內(nèi)外無(wú)人機(jī)的管制一般都采用劃設(shè)隔離空域的方式運(yùn)行［2］，限制了無(wú)人機(jī)在各領(lǐng)域的應(yīng)用，為使無(wú)人機(jī)能充分發(fā)揮它的作用，找到一種可靠的飛行沖突解脫技術(shù)成為目前無(wú)人機(jī)系統(tǒng)設(shè)計(jì)的關(guān)鍵技術(shù)之一［3］。

目前，國(guó)內(nèi)外對(duì)無(wú)人機(jī)沖突解脫技術(shù)的研究有很多，但是絕大多數(shù)研究都集中在沖突解脫算法的有效性上，從算法的可靠性以及無(wú)人機(jī)空域飛行安全的角度檢驗(yàn)沖突解脫算法性能的研究很少。在國(guó)外，M.Garcia［4］，Eva Crück［5］在沖突解脫算法上都提出了自己的觀點(diǎn)，但研究中并未對(duì)算法的可靠性上未加證明，很大程度上降低了算法應(yīng)用于實(shí)際工程的可能性。Jose Asmat［6］設(shè)計(jì)了一種利用TCAS（Traffic Alert Collision Avoidance System）技術(shù)來(lái)進(jìn)行無(wú)人機(jī)協(xié)同方式下的沖突解脫技術(shù)，并用Monte Carlo的仿真驗(yàn)證了沖突解脫方案的安全性能；在國(guó)內(nèi)方面對(duì)無(wú)人機(jī)評(píng)估主要集中于任務(wù)執(zhí)行效率方面，王淵［7］、崔莉薇［8］等分別應(yīng)用不同的智能算法，實(shí)現(xiàn)了飛行沖突解脫，但僅僅比較了解脫航跡的長(zhǎng)度和算法運(yùn)算速度，未對(duì)算法的可靠性進(jìn)行說(shuō)明；也有使用風(fēng)險(xiǎn)評(píng)價(jià)指數(shù)矩陣法［9］、層次分析法［10］、多級(jí)模糊綜合評(píng)判法［11］等方法進(jìn)行無(wú)人機(jī)安全評(píng)估，但相對(duì)數(shù)量較少。

鑒于此，本文基于蟻群算法提出了一種無(wú)人機(jī)飛行沖突解脫技術(shù)，再利用蒙特卡洛法（MC，Monte Carlo）對(duì)改進(jìn)的算法進(jìn)行可靠性分析，從安全的角度對(duì)算法的性能進(jìn)行比較。

1 蒙特卡洛方法簡(jiǎn)介

隨著現(xiàn)代計(jì)算機(jī)技術(shù)的飛速發(fā)展，蒙特卡洛法已經(jīng)成為人們解決復(fù)雜統(tǒng)計(jì)模型和隨機(jī)事件模擬的主要工具。蒙特卡洛法又稱(chēng)為隨機(jī)模擬法或統(tǒng)計(jì)模擬法，屬于計(jì)算數(shù)學(xué)的一個(gè)分支，其實(shí)質(zhì)是按一定概率分布產(chǎn)生隨機(jī)數(shù)的方法來(lái)模擬實(shí)際中可能出現(xiàn)的隨機(jī)現(xiàn)象，通過(guò)利用計(jì)算機(jī)軟件進(jìn)行大規(guī)模仿真實(shí)驗(yàn)，得到模擬的實(shí)驗(yàn)數(shù)據(jù)，在此實(shí)驗(yàn)的基礎(chǔ)上進(jìn)行分析判斷，得到實(shí)際問(wèn)題中可能出現(xiàn)的規(guī)律或求解方法。在現(xiàn)實(shí)生活中進(jìn)行工程實(shí)驗(yàn)時(shí)，試驗(yàn)有可能高風(fēng)險(xiǎn)，耗費(fèi)高，耗時(shí)長(zhǎng)甚至具有破壞性，不可能通過(guò)反復(fù)試驗(yàn)進(jìn)行實(shí)驗(yàn)數(shù)據(jù)的采集，此時(shí)蒙特卡洛法不失為最簡(jiǎn)單、經(jīng)濟(jì)、實(shí)用的方法。

本文研究的是無(wú)人機(jī)飛行沖突解脫問(wèn)題，不可能在安全性未達(dá)標(biāo)的情況下在現(xiàn)實(shí)中進(jìn)行沖突解脫效率測(cè)試，因而采取應(yīng)用計(jì)算機(jī)軟件進(jìn)行仿真模擬。飛機(jī)在飛行時(shí)由于設(shè)備精度，飛行員操作誤差，環(huán)境影響等因素，飛行軌跡具有隨機(jī)性，具有蒙特卡洛方法的特征，故而應(yīng)用蒙特卡洛法進(jìn)行仿真，因涉及安全的因素十分復(fù)雜，本文將所有因素，在綜合作用下可以假定為高斯分布［13］。進(jìn)行大規(guī)模的仿真，得出無(wú)人機(jī)沖突解脫算法成功避險(xiǎn)的概率，并對(duì)結(jié)果進(jìn)行分析。

2 改進(jìn)蟻群算法的沖突解脫技術(shù)

蟻群算法作為一種最新提出的仿生進(jìn)化算法具有:魯棒性強(qiáng)、并行計(jì)算能力、全局搜索能力、正反饋機(jī)制、易于與其他算法混合等優(yōu)點(diǎn)［14］。目前為止，蟻群算法已經(jīng)成功解決了旅行商問(wèn)題、指派問(wèn)題、調(diào)度問(wèn)題、車(chē)輛路由問(wèn)題、網(wǎng)絡(luò)路由問(wèn)題等組合優(yōu)化問(wèn)題，且與旅行商問(wèn)題具有一定的相似性，但由于蟻群算法具有搜索時(shí)間長(zhǎng)、易陷入局部解兩大缺點(diǎn)，必須要加以改善才能適應(yīng)無(wú)人機(jī)沖突解脫的復(fù)雜性與對(duì)時(shí)效性的要求。

在傳統(tǒng)蟻群算法中，一般將節(jié)點(diǎn)和節(jié)點(diǎn)之間的距離的倒數(shù)作為啟發(fā)因子，如式（1）:

式中:dij為節(jié)點(diǎn)i到節(jié)點(diǎn)j的距離。螞蟻k由狀態(tài)i向j轉(zhuǎn)移的概率選擇公式為:

式中:Pij（kt）為t時(shí)刻螞蟻k從節(jié)點(diǎn)i轉(zhuǎn)移至j的概率為t時(shí)刻節(jié)點(diǎn)i到j(luò)路徑上的信息素濃度；為節(jié)點(diǎn)j的啟發(fā)函數(shù)；allowedk為螞蟻k下一步可能到達(dá)的路徑點(diǎn)集；α，β分別表示的作用程度。

在所有螞蟻進(jìn)行完一次路徑規(guī)劃后，則需要進(jìn)行信息素的更新，信息素更新公式為:

在無(wú)人機(jī)路徑規(guī)劃問(wèn)題中，有時(shí)候最優(yōu)解與次優(yōu)解之間會(huì)相差較大距離，甚至不在同一方向上。最大最小螞蟻系統(tǒng)中只有最優(yōu)螞蟻才釋放信息素，將所有螞蟻趨向最優(yōu)螞蟻附近，當(dāng)最優(yōu)螞蟻附近的解與真實(shí)的最優(yōu)解之間差距較大時(shí)，最優(yōu)解就有可能被錯(cuò)過(guò)，對(duì)此，引入基于排序的螞蟻系統(tǒng)。每次迭代完成，對(duì)所有螞蟻進(jìn)行排序，除最優(yōu)螞蟻釋放信息素外，適應(yīng)度為前γ的螞蟻都進(jìn)行信息素的釋放，信息素更新公式為:

在完成一次信息素更新之后，一次循環(huán)結(jié)束，對(duì)形成的路徑進(jìn)行適應(yīng)性評(píng)價(jià)，適應(yīng)性函數(shù)為:

其中:w1，w2，w3為權(quán)重；J為綜合代價(jià)；Js為矩形空域威脅代價(jià)；Jr為圓形空域威脅代價(jià)；Ja為動(dòng)態(tài)航空器的綜合威脅代價(jià)；Jl為無(wú)人機(jī)的航程約束代價(jià)。

無(wú)人機(jī)的飛行沖突解脫運(yùn)用了動(dòng)態(tài)路徑規(guī)劃的原理，在每完成一次路徑規(guī)劃之后，進(jìn)行一次狀態(tài)轉(zhuǎn)移，并對(duì)空域態(tài)勢(shì)進(jìn)行再次探測(cè)與更新，基于最新的態(tài)勢(shì)信息，無(wú)人機(jī)進(jìn)行航路的重規(guī)劃，不斷滾動(dòng)前進(jìn)，直至到達(dá)終點(diǎn)。

本文構(gòu)建的沖突場(chǎng)景為航線交叉相遇下的飛行沖突，結(jié)合我國(guó)當(dāng)前航空器運(yùn)行安全間隔的劃分標(biāo)準(zhǔn)［15］，根據(jù)無(wú)人機(jī)每次飛行與其他航空器的最小間隔距離，對(duì)無(wú)人機(jī)運(yùn)行的安全狀況進(jìn)行分類(lèi)。當(dāng)無(wú)人機(jī)與有人機(jī)之間的飛行間隔小于1.76 km，則認(rèn)為無(wú)人機(jī)存在極高的相撞風(fēng)險(xiǎn)，記為一次嚴(yán)重事故。當(dāng)航空器間的距離在10 km～15 km之間時(shí)，則認(rèn)為無(wú)人機(jī)存在發(fā)生飛行沖突的潛在風(fēng)險(xiǎn)。在10 km～15 km的間隔范圍內(nèi)雖然不算一次事故，但對(duì)其進(jìn)行統(tǒng)計(jì)有利于無(wú)人機(jī)沖突解脫反應(yīng)距離的設(shè)置。具體的事故風(fēng)險(xiǎn)等級(jí)劃分如表1所示。

表1 事故風(fēng)險(xiǎn)等級(jí)劃分

改進(jìn)算法的飛行沖突避讓流程如圖1所示。

3 仿真結(jié)果

本文利用Matlab2008進(jìn)行仿真，在100km×100km的空域范圍內(nèi)分別構(gòu)造了單機(jī)沖突場(chǎng)景和雙機(jī)沖突場(chǎng)景，單機(jī)沖突分為45°交叉角和90°交叉角兩種情況，雙機(jī)沖突為兩種角度的綜合。有人機(jī)駕駛的航空器按照飛行計(jì)劃飛行。無(wú)人機(jī)在空域運(yùn)行過(guò)程中需要避開(kāi)空域中運(yùn)行的航空器及各類(lèi)不可進(jìn)入的空域，T代表訓(xùn)練空域，D代表空中危險(xiǎn)區(qū)，F(xiàn)代表空中禁區(qū)，Start為無(wú)人機(jī)某一航段起始點(diǎn)，Goal為該航段終點(diǎn)。無(wú)人機(jī)沖突解脫安全性評(píng)估的運(yùn)行環(huán)境及飛行路徑如圖2所示。

無(wú)人機(jī)沖突解脫過(guò)程中各類(lèi)航空器的綜合軌跡如圖3所示。

圖1 基于改進(jìn)算法的飛行沖突避讓流程圖

圖2 無(wú)人機(jī)空域運(yùn)行環(huán)境

圖3 無(wú)人機(jī)沖突解脫示意圖

本文評(píng)估對(duì)象為無(wú)人機(jī)自主沖突解脫技術(shù)，重點(diǎn)是評(píng)估其安全性水平，因此，管制員與無(wú)人機(jī)操作員的影響因素被去除。而有人機(jī)飛行員并不與無(wú)人機(jī)進(jìn)行飛行沖突解脫的協(xié)同，將會(huì)駕駛飛機(jī)按計(jì)劃飛行，航空器在運(yùn)行過(guò)程中可能會(huì)受到風(fēng)、人為因素影響、設(shè)備運(yùn)行誤差的擾動(dòng)，在綜合作用下可假定為高斯分布。

本文分別對(duì)基本蟻群算法和改進(jìn)蟻群算法的沖突解脫技術(shù)進(jìn)行模擬，分為單機(jī)沖突和雙機(jī)沖突，以檢驗(yàn)各算法對(duì)空域復(fù)雜度增加的情況下的沖突解脫安全性水平，其中單機(jī)沖突按交叉角分為45°交叉角和90°交叉角兩種情況，針對(duì)每一種情況分別進(jìn)行10萬(wàn)次的蒙特卡洛仿真實(shí)驗(yàn)，測(cè)試沖突解脫算法的安全性水平，得出的數(shù)據(jù)如下頁(yè)表2、表3。

對(duì)比表2與表3可以發(fā)現(xiàn)，基本蟻群算法與改進(jìn)蟻群算法在90°交叉角上很好地防止了嚴(yán)重事故和危險(xiǎn)接近的發(fā)生，但是改進(jìn)算法在飛行沖突和存在風(fēng)險(xiǎn)兩個(gè)危險(xiǎn)等級(jí)上改善不明顯。在單機(jī)45°交叉角和雙機(jī)沖突上，改進(jìn)蟻群算法對(duì)比基本算法在事故率上有了很明顯的減少，但按照美軍提出的無(wú)人機(jī)運(yùn)行安全目標(biāo)水平105 h內(nèi)20次事故的標(biāo)準(zhǔn)還有較大的差距［15］。對(duì)比ICAO提出的10-7次事故/h的標(biāo)準(zhǔn)則可判定無(wú)人機(jī)的安全水平無(wú)法達(dá)到有人機(jī)的運(yùn)行標(biāo)準(zhǔn)。說(shuō)明改進(jìn)的算法還需針對(duì)這一情況進(jìn)行進(jìn)一步的改善，來(lái)降低無(wú)人機(jī)沖突解脫過(guò)程中的風(fēng)險(xiǎn)。

表2 基于基本蟻群算法事故率

表3 改進(jìn)算法下的事故率

4 結(jié)論

當(dāng)前國(guó)內(nèi)有關(guān)無(wú)人機(jī)沖突解脫技術(shù)的文獻(xiàn)較多，但大部分研究都集中在沖突解脫技術(shù)的有效性上，很少有針對(duì)沖突解脫安全性方面的驗(yàn)證，并且目前國(guó)內(nèi)針對(duì)無(wú)人機(jī)空域運(yùn)行安全評(píng)估體系還很缺乏。因此，利用蒙特卡洛方法針對(duì)無(wú)人機(jī)飛行沖突解脫技術(shù)的安全性進(jìn)行了模擬仿真，仿真結(jié)果有較大改善，但安全性仍然無(wú)法滿足空域運(yùn)行的要求。得出無(wú)人機(jī)對(duì)45°交叉相遇下的沖突解脫能力存在嚴(yán)重的不足，應(yīng)當(dāng)作為安全瓶頸加以重點(diǎn)研究。本文在構(gòu)建無(wú)人機(jī)空域運(yùn)行環(huán)境的過(guò)程中作了很多簡(jiǎn)化，與實(shí)際環(huán)境具有較大差異，但重點(diǎn)在于引入蒙特卡洛法針對(duì)無(wú)人機(jī)飛行沖突解脫技術(shù)的安全性進(jìn)行模擬，從安全的角度考慮沖突解脫算法的性能，可以為未來(lái)無(wú)人機(jī)空域運(yùn)行安全評(píng)估體系的構(gòu)建提供一種理論參考。

［1］陳金良，高文明，史際剛.無(wú)人機(jī)飛行管理［M］.西安:西北工業(yè)大學(xué)出版社，2014:1-2.

［2］李春錦，文涇.無(wú)人機(jī)系統(tǒng)的運(yùn)行管理［M］.北京:北京航空航天大學(xué)出版社，2011.5.

［3］SHIN H S，ANTONIOS T，WHITE B A，et al.UAV Conflict Detection and Resolution forStatic and Dynamic Obstacles［C］//Guidance，Navigation and Control Conference and Exhibit.Honolulu，Hawaii:American Institute of Aeronautics and Astronautics，2008.

［4］GARCIA M，VIGURIA A，OLLERO A.Dynamic graph-search algorithm for global path planning in presence of hazardous weather［J］.J Intell Robot Syst，2013，69:285-295.

［5］Eva Crück，JOHN L.Sense and avoid system for a MALE UAV［C］//Guidance，Navigation and Control Conference and Exhibit.Hilton Head，South Carolina:American Institute of Aeronautics and Astronautics，2007.

［6］JOSE A，BRETT R，JESICA U，et al.Unmanned aerial collision avoidance system［R］.George Mason University:Department of Systems Engineering and Operations Research，2006.

［7］王淵，孫秀霞，劉樹(shù)光，等.基于改進(jìn)人工蜂群算法的多機(jī)飛行沖突解脫策略［J］.空軍工程大學(xué)學(xué)報(bào)（自然科學(xué)版），2014，15（3）:10-14.

［8］崔莉薇，石為人，劉祥明，等.基于遺傳粒子群算法的飛行沖突解脫［J］.計(jì)算機(jī)工程與應(yīng)用，2013，49（7），263-266.

［9］陽(yáng)再清，劉志強(qiáng).艦載無(wú)人機(jī)飛行安全性評(píng)估方法研究［J］.宇航計(jì)測(cè)技術(shù)，2014，34（6）:84-89.

［10］魏法杰，盧珊，張羽.某實(shí)際型號(hào)無(wú)人機(jī)試飛風(fēng)險(xiǎn)評(píng)估案例［J］.中國(guó)管理科學(xué)，2014，22（11）:199-204.

［11］吳忠杰，張耀中，王強(qiáng).無(wú)人機(jī)復(fù)雜氣象威脅建模及評(píng)估方法［J］.計(jì)算機(jī)應(yīng)用，2013，33（4）:1179-1182.

［12］魏瀟龍，姚登凱，戴喆，等.基于路徑規(guī)劃的無(wú)人機(jī)飛行沖突解脫技術(shù)［J］.火力與指揮控制，2016，41（10）: 48-54.

［13］閆麗萍，夏元清，楊毅.隨機(jī)過(guò)程理論及其在自動(dòng)控制中的應(yīng)用［M］.北京:國(guó)防工業(yè)出版社，2012.4

［14］宋錦娟.一種改進(jìn)的蟻群算法及其在最短路徑問(wèn)題中的應(yīng)用［D］.太原:中北大學(xué)，2013.

［15］張兆寧，王莉莉，李冬賓.飛行間隔安全評(píng)估引論［M］.北京:科學(xué)出版社，2009.

UAV Flight Conflict Resolution Safety Assessment Technique Based on Monte Carlo Method

GU Zhi-ming1，GAO Wen-ming2，WEI Xiao-long2，YAO Deng-kai2
（1.Unit 91899 of PLA，Huludao 125001，China；2.School of Air Traffic Control and Navigation，Air Force Engineering University，Xi’an 710051，China）

Conflict resolution technology is a prerequisite for UAV.Aiming at the lack of the safety assessment method of conflict resolution technology，this paper uses Monte Carlo method to simulate the flight conflict and reached the efficiency of conflict resolution algorithm.In this paper，an improved ant colony algorithm is proposed to achieve the flight conflict resolution；using the Monte Carlo method to simulate the flight conflict process and then compare the efficiency with basic ant colony algorithm.The ideas presented in this paper can provide a theoretical reference for the construction of safety assessment system in the UAV airspace operation in the future.

UAV，ant colony algorithm，Monte Carlo，conflict resolution

V271.4

A

1002-0640（2017）04-0158-04

2016-03-16

2016-04-27

省部級(jí)自然科學(xué)基金資助項(xiàng)目（2015JM7364）

谷志鳴（1988-），男，河北無(wú)極人，碩士研究生。研究方向：空域規(guī)劃與空中交通流量管理。