一種空中交通管制系統的飛行衝突預警方法

2023-06-18 09:30:16

一種空中交通管制系統的飛行衝突預警方法
【專利摘要】本發明涉及一種空中交通管制系統的飛行衝突預警方法，所述空中交通管制系統包括數據通信模塊、監視數據融合模塊、機載終端模塊、管制終端模塊，其中監視數據融合模塊用於實現空管雷達監視數據與自動相關監視數據的融合，為管制終端模塊提供實時航跡信息；管制終端模塊包括飛行前無衝突4D航跡生成、飛行中短期4D航跡生成、實時飛行衝突監控與告警、這3個子模塊；上述系統的飛行衝突預警方法，依靠管制終端模塊，處理飛行計劃數據並利用隱馬爾科夫模型生成4D航跡，實現空域交通狀況潛在的交通衝突的分析。本發明可有效預警飛行衝突，提高空中交通的安全性。
【專利說明】一種空中交通管制系統的飛行衝突預警方法

【技術領域】
[0001] 本發明涉及一種空中交通管制系統及方法，尤其涉及一種基於4D航跡運行的空 中交通管制系統對航空器軌跡進行預測並對飛行衝突預警的方法。

【背景技術】
[0002] 隨著全球航空運輸業快速發展與空域資源有限矛盾的日益突出，在空中交通流密 集的複雜空域，仍然採用飛行計劃結合間隔調配的空中交通管理方式逐漸顯示出其落後 性，具體表現在：（1)飛行計劃並未為航空器配置精確的空管間隔，容易造成交通流戰術管 理中的擁擠，降低空域安全性；(2)以飛行計劃為中心的空管自動化系統對飛行剖面的推 算和航跡預測精度差，造成衝突化解能力差；(3)空中交通管制工作仍然側重於保持單個 航空器之間的安全間隔，很難上升到對交通流進行戰略性管理。對於航空器軌跡的預測並 從而並對飛行衝突預警顯得尤為重要。
[0003] 4D航跡是以空間和時間形式，對某一航空器航跡中的各點空間位置（經度、煒度 和高度）和時間的精確描述，基於航跡的運行是指在4D航跡的航路點上使用"控制到達 時間"，即控制航空器通過特定航路點的"時間窗"。在高密度空域把基於4D航跡的運行 (TrajectorybasedOperation)作為基本運行機制之一，是未來對大流量、高密度、小間隔 條件下空域實施管理的一種有效手段，可以顯著地減少航空器航跡的不確定性，提高空域 和機場資源的安全性與利用率。
[0004] 基於航跡運行的空中交通運行方式需要在戰略層面上對單航空器飛行航跡進行 推算和優化，對多航空器構成的交通流實施協同和調整；在預戰術層面上通過修正交通流 中個別航空器的航跡以解決擁塞問題，並保證該交通流中所有航空器的運行效率；而在戰 術層面上預測衝突和優化解脫方案，則非常依賴於能否準確地對航空器的軌跡進行預測並 對飛行衝突進行預警，目前均不能準確實時地對航空器的軌跡進行預測，實時性上做的尤 為的差。


【發明內容】

[0005] 本發明要解決的技術問題是在於克服現有技術的不足，提供一種基於4D航跡運 行的空中交通管制系統的飛行衝突預警方法，可有效、準確、實時地預測航空器的軌跡並預 測飛行衝突。
[0006] 實現本發明目的的技術方案是提供一種空中交通管制系統的飛行衝突預警方法， 所述空中交通管制系統包括機載終端模塊、數據通信模塊、監視數據融合模塊以及管制終 端模塊；監視數據融合模塊用於實現空管雷達監視數據與自動相關監視數據的融合，為管 制終端模塊提供實時航跡信息；
[0007] 所述管制終端模塊包括以下子模塊：
[0008] 飛行前無衝突4D航跡生成模塊，根據飛行計劃和世界區域預報系統的預報數據， 建立航空器動力學模型，然後依據飛行衝突耦合點建立航跡衝突預調配理論模型，生成航 空器無衝突4D航跡；
[0009] 飛行中短期4D航跡生成模塊，依據監視數據融合模塊提供的實時航跡信息，利用 隱馬爾科夫模型，推測未來一定時間窗內的航空器4D軌跡；
[0010] 實時飛行衝突監控與告警模塊，用於建立從航空器的連續動態到離散衝突邏輯的 觀測器，將空中交通系統的連續動態映射為離散觀測值表達的衝突狀態；當系統有可能違 反空中交通管制規則時，對空中交通混雜系統的混雜動態行為實施監控，為管制員提供及 時的告警信息；
[0011] 所述空中交通管制系統的飛行衝突預警方法包括如下幾個步驟：
[0012] 步驟A、飛行前無衝突4D航跡生成模塊根據飛行計劃和世界區域預報系統的預報 數據，建立航空器動力學模型，並依據飛行衝突耦合點建立航跡衝突預調配理論模型，生成 航空器無衝突4D航跡；
[0013] 步驟B、監視數據融合模塊將空管雷達監視數據與自動相關監視數據進行融合， 生成航空器實時航跡信息並提供給管制終端模塊；管制終端模塊中的飛行中短期4D航跡 生成模塊依據航空器實時航跡信息和歷史航跡信息推測未來一定時間窗內的航空器4D軌 跡；所述依據航空器實時航跡信息和歷史航跡信息推測未來一定時間窗內的航空器4D軌 跡的具體實施過程如下：
[0014] 步驟B6、對航空器軌跡數據預處理，依據所獲取的航空器原始離散二維位置序列 X = [X1, X2,…，Xn]和y=Iiy1,y2，…，yn]，採用一階差分方法對其進行處理獲取新的航空器 離散位置序列ΛX= [ΛX1,Λχ2,…，ΛXn-J和Λy= [Λy"Λy2,…，ΛYn-J,其中 Δxb=χb+「xb，Δyb=yb+「yb(b= 1，2，…，n-1);
[0015] 步驟B7、對航空器軌跡數據聚類，對處理後新的航空器離散二維位置序列Λx和 Λy，通過設定聚類個數Μ'，採用遺傳聚類算法分別對其進行聚類；
[0016] 步驟B8、對聚類後的航空器軌跡數據利用隱馬爾科夫模型進行參數訓練，通過將 處理後的航空器運行軌跡數據ΛX和Λy視為隱馬爾科夫過程的顯觀測值，通過設定隱狀 態數目Ν'和參數更新時段ζ'，依據最近的Τ'個位置觀測值並採用B-W算法滾動獲取最新 隱馬爾科夫模型參數λ' ；
[0017] 步驟Β9、依據隱馬爾科夫模型參數，採用Viterbi算法獲取當前時刻觀測值所對 應的隱狀態q;
[0018] 步驟B10、通過設定預測時域h'，基於航空器當前時刻的隱狀態q，獲取未來時段 航空器的位置預測值O;
[0019] 步驟C、實時飛行衝突監控與告警模塊建立從航空器的連續動態到離散衝突邏輯 的觀測器，將空中交通系統的連續動態映射為離散觀測值表達的衝突狀態；當系統有可能 違反空中交通管制規則時，對空中交通混雜系統的混雜動態行為實施監控，為管制員提供 及時的告警信息。
[0020] 進一步的，步驟B中，所述聚類個數M'的值為4,隱狀態數目Ν'的值為3,參數更 新時段ζ'為30秒，Τ'為10,預測時域h'為300秒。
[0021] 進一步的，步驟B的B8具體是指：由於所獲得的航跡序列數據長度是動態變化 的，為了實時跟蹤航空器航跡的狀態變化，有必要在初始航跡隱馬爾科夫模型參數λ' = (π，Α，Β)的基礎上對其重新調整，以便更精確地推測航空器在未來某時刻的位置；每隔時 段ζ'，依據最新獲得的τ'個觀測值（〇1，〇2，…，〇τ，）對航跡隱馬爾科夫模型參數λ' =O，A,B)進行重新估計。
[0022]步驟B的BlO具體是指：每隔時段Φ，根據最新獲得的隱馬爾科夫模型參數λ' = (π，Α，Β)和最近H個歷史觀測值（〇1，〇2，…，〇Η)，基於航空器當前時刻的隱狀態q，通過設定 預測時域h'，在時刻t獲取航空器在未來時段h'的位置預測值0。
[0023] 更進一步的，時段Φ為4秒。
[0024] 進一步的，所述步驟A的航空器無衝突4D航跡按照以下方法生成：
[0025] 步驟Α1、進行航空器狀態轉移建模，根據飛行計劃中航空器的飛行高度剖面，建立 單個航空器在不同航段轉移的Petri網模型：E= (g，G，Pre，Post，m)為航空器階段轉移模 型，其中g表示飛行航段，G表示垂直剖面中飛行狀態參數的轉換點，Pre和Post分別表示 航段和航路點的前後向連接關係，^Z+表示航空器所處的飛行階段；
[0026] 步驟A2、建立航空器全飛行剖面混雜系統模型如下，
[0027]vH=k(vcas,Mach,hp,tL0C),
[0028]vGS=μ(vCAS,Mach,hp,tL0C,vffS,a),
[0029] 其中veAS為校正空速，Mach為馬赫數，hp為氣壓高度，a為風向預報與航路的夾 角，Vws為風速預報值，t為溫度預報值，VA高度變化率，V<^為地速；
[0030] 步驟A3、採用混雜系統仿真的方式推測求解航跡：採用將時間細分的方法， 利用狀態連續變化的特性遞推求解任意時刻航空器在某一飛行階段距參考點的航程 J(X) =J。+ 和高度= 其中Jtl為初始時刻航空器距參考點 的航程，Λτ為時間窗的數值，J(T)為τ時刻航空器距參考點的航程，Iitl為初始時刻航 空器距參考點的高度，1ι(τ)為τ時刻航空器距參考點的高度，由此可以推測得到單航空 器的4D航跡；
[0031] 步驟Α4、對多航空器耦合模型實施無衝突調配：根據兩航空器預達交叉點的時 間，按照空中交通管制原則，對交叉點附近不滿足間隔要求的航空器4D航跡進行二次規 劃，得到無衝突4D航跡。
[0032] 進一步的，所述步驟B中監視數據融合模塊將空管雷達監視數據與自動相關監視 數據進行融合，生成航空器實時航跡信息，具體按照以下方法：
[0033] 步驟Β1、將坐標單位和時間統一；
[0034] 步驟Β2、採用最鄰近數據關聯算法將屬於同一個目標的點相關聯，提取目標航跡； 步驟Β3、將分別從自動相關監視系統和空管雷達提取的航跡數據從不同的時空參
[0035] 考坐標系統變換、對準到管制終端統一的時空參考坐標系統；
[0036] 步驟M、計算兩條航跡的相關係數，若相關係數小於某一預設閾值，則認為兩條航 跡不相關；否則該兩條航跡相關，可以進行融合；
[0037] 步驟Β5、對相關的航跡進行融合。
[0038] 更進一步的，所述步驟Β5中對相關的航跡進行融合，採用基於採樣周期的加權平 均算法，其加權係數根據採樣周期和信息精度確定，再利用加權平均算法將與之相關的自 動相關監視航跡和空管雷達航跡融合為系統航跡。
[0039] 進一步的，所述步驟C的具體實施過程如下：
[0040] 步驟CU構造基於管制規則的衝突超曲面函數集：建立超曲面函數集用以反映系 統的衝突狀況，其中，衝突超曲面中與單一航空器相關的連續函數為第I型超曲 面，與兩架航空器相關的連續函數\ :YXX^R為第II型超曲面；
[0041] 步驟C2、建立由航空器連續狀態至離散衝突狀態的觀測器：需要根據管制規範建 立觀測器，觀測系統系統穿越超曲面而產生的衝突事件，以便控制器做出相應的控制決策 指令；觀測器ξ用於觀測系統中航空器位置的連續變化而產生衝突事件，稱為 第I型觀測器，4S為第II型觀測器；
[0042] 步驟C3、設計從衝突到衝突解脫手段的離散監控器，該離散監控器可描述為函數 βD，其中S是觀測器觀測向量展成的空間，D是所有決策向量d展成的空間；當觀測 器的離散觀測向量表明某一非期望的狀態出現時，立刻發出相應的告警。
[0043] 本發明具有積極的效果：（1)本發明的一種空中交通管制系統的飛行衝突預警方 法在航空器實時軌跡推測過程中，融入了隨機因素的影響，所採用的滾動軌跡推測方案能 夠及時提取外界隨機因素的變化狀況，提高了航空器軌跡推測的準確性。
[0044] (2)本發明的一種空中交通管制系統的飛行衝突預警方法對飛行衝突的預警效果 較好，可有效、準確、實時地預測航空器的軌跡並預測飛行衝突。
[0045] (3)本發明的一種空中交通管制系統的飛行衝突預警方法對飛行剖面的推算和航 跡預測精度高，進而使得衝突化解能力和自動化水平提高，降低了管制員的工作負荷。

【專利附圖】

【附圖說明】
[0046] 圖1為飛行前無衝突4D航跡生成方法流程示意圖；
[0047] 圖2為飛行中短期4D航跡推測方法流程示意圖；
[0048] 圖3為航空器航跡衝突監控與告警方法流程示意圖。

【具體實施方式】
[0049] (實施例1)
[0050] 本實施例的基於4D航跡運行的空中交通管制系統，包括機載終端模塊101、數據 通信模塊102、監視數據融合模塊103以及管制終端模塊104。以下對各部分的具體實施方 式分別進行詳細描述。
[0051] 1.機載終端模塊
[0052] 機載終端模塊101是飛行員獲取地面管制指令、參考4D航跡，以及輸入飛行意圖 的界面，同時還是採集當前航空器位置數據的接口。
[0053]其具體實施方案如下：
[0054] 機載終端模塊101接收如下的信息輸入：（I)ADS-B信息採集單元201通過機載 GPS採集的航空器位置向量、速度向量，以及本航空器的呼號，編碼後通過信息及數據傳遞 給機載數據通信模塊102 ; (2)航空器駕駛員需要將與地面管制指令不一致的飛行意圖，通 過人機輸入界面，以及約定的地面管制員可以識別的形式通過信息及數據傳遞給機載數據 通信模塊102。另外機載終端模塊101實現如下的信息輸出：（1)通過終端顯示屏幕，接收 和顯示飛行員可以識別的飛行管制指令；（2)接收和顯示地面管制終端飛行前生成的無衝 突4D航跡，以及當地面管制終端探測到衝突後計算的最優解脫4D航跡。
[0055] 2.數據通信模塊
[0056] 數據通信模塊102可實現空地雙向數據通信，實現機載實時位置數據和飛行意圖 數據單元202的下行傳輸和地面管制指令單元203,以及參考4D航跡單元204的上行傳輸。
[0057]其具體實施方案如下：
[0058] 下行數據通信：機載終端101通過機載二次雷達應答機將航空器識別標誌和4D 位置信息，以及其他附加數據，如飛行意圖、飛行速度、氣象等信息傳輸給地面二次雷達 (SSR)，二次雷達接收後對數據報文進行解析，並傳輸給中央數據處理組件301解碼，通過 指令航跡數據接口傳輸到管制終端104;上行數據通信：地面管制終端104通過指令航跡數 據接口，經中央數據處理組件301編碼後，地面二次雷達的詢問機將將地面管制指令或參 考4D航跡信息傳遞並顯示在機載終端101。
[0059] 3.監視數據融合模塊
[0060] 監視數據融合模塊103實現空管雷達監視與自動相關監視ADS-B數據的融合，為 管制終端模塊104中的飛行中短期4D航跡生成子模塊和實時飛行衝突監控與告警子模塊 提供實時航跡信息。
[0061] 其具體實施方案如下：
[0062] (1)在預處理階段將坐標單位和時間統一，假設分別從ADS-B和空管雷達中提取 的數據是一系列離散點的坐標（如經度、煒度、海拔高度）、各點對應採集時間；（2)採用最 鄰近數據關聯算法將屬於同一個目標的點相關聯，提取目標航跡；(3)將分別從ADS-B和空 管雷達中提取的航跡數據從不同的時空參考坐標系統變換、對準到管制終端統一的時空參 考坐標系統；(4)計算兩條航跡的相關係數，若相關係數小於某一預設閾值，則認為兩條航 跡不相關，否則該兩條航跡相關，可以進行融合；（5)對相關的航跡進行融合。由於ADS-B和 空管雷達的精度和採樣周期不同，本系統採用基於採樣周期的加權平均算法，其加權係數 根據採樣周期和信息精度確定，再利用加權平均算法將與之相關的ADS-B航跡和空管雷達 航跡融合為系統航跡。
[0063] 4.管制終端模塊
[0064] 管制終端模塊104包括飛行前無衝突4D航跡生成、飛行中短期4D航跡生成、實時 飛行衝突監控與告警這3個子模塊。
[0065] (1)飛行前無衝突4D航跡生成
[0066] 根據飛行數據處理系統（FDP)得到的飛行計劃和世界區域預報系統（WAFS)發布 的風、溫度的GRIB格點預報數據，對空中交通系統建立層次化的混雜系統模型，通過系統 在安全狀態的演化，描述狀態演化的時間軌跡，生成航空器航跡。
[0067] 如圖1所示，其具體實施過程如下：
[0068] 首先，進行航空器狀態轉移建模。航空器沿航跡飛行的過程表現為在航段之間動 態切換過程，根據飛行計劃中航空器的飛行高度剖面，建立單個航空器在不同航段轉移的 Petri網模型：E= (g，G,Pre,Post,m)為航空器階段轉移模型，其中g表示飛行航段，G表 示垂直剖面中飛行狀態參數（包括空速、高度、構型）的轉換點，Pre和Post分別表示航段 和航路點的前後向連接關係，《I^Z1表示航空器所處的飛行階段。
[0069]其次，建立航空器全飛行剖面混雜系統模型。航空器在單個航段內的飛行視為連 續過程，依據質點能量模型，推導航空器在不同的運行階段同氣象條件下的航空器動力學 方程，νΗ=κ(VCAS,Mach,hp,tLQC)，Vgs=μ(VCAS,Mach,hp,tLQC,vws,α)，其中VcAS為校正空速， Mach為馬赫數，hp為氣壓高度，α為風向預報與航路的夾角，Vws為風速預報值，為溫度 預報值，vH為高度變化率，V<^為地速。
[0070] 然後，採用混雜系統仿真的方式推測求解航跡。採用將時間細分的方法，利 用狀態連續變化的特性遞推求解任意時刻航空器在某一飛行階段距參考點的航程 = +f)Λ和高度Α(τ) =Α。+fn?⑴也，其中Jtl為初始時刻航空器距參考點 JO JO 的航程，ΛT為時間窗的數值，J(T)為τ時刻航空器距參考點的航程，Iitl為初始時刻航 空器距參考點的高度，1ι(τ)為τ時刻航空器距參考點的高度，由此可以推測得到單航空 器的4D航跡。
[0071] 最後，對多航空器耦合模型實施無衝突調配。根據兩航空器預達交叉點的時間，按 照空中交通管制原則，對交叉點附近不滿足間隔要求的航空器4D航跡進行二次規劃，得到 無衝突4D航跡。
[0072] (2)飛行中短期4D航跡生成
[0073] 依據管制雷達和自動相關監視系統ADS-B實施融合後獲得航空器實時航跡數據， 利用隱馬爾科夫模型，推測未來5分鐘時間窗內的航空器4D軌跡。
[0074] 如圖2所示，其具體實施過程如下：
[0075] 首先，對航空器軌跡數據預處理，依據所獲取的航空器原始離散二維位置序列X =[XpX2,…，Xn]和y=Iiy^y2，…，yn]，採用一階差分方法對其進行處理獲取新的航空器 離散位置序列ΛX= [ΛX1,Λχ2,…，ΛXn-J和Λy= [Λy"Λy2,…，ΛYn-J,其中 AXb-Xb+l_Xb，AYb-yb+l_yb(b- 1，2，...，n_l)。
[0076] 其次，對航空器軌跡數據聚類。對處理後新的航空器離散二維位置序列ΛX和 Λy，通過設定聚類個數Μ'，採用遺傳聚類算法分別對其進行聚類。
[0077] 然後，對聚類後的航空器軌跡數據利用隱馬爾科夫模型進行參數訓練。通過將處 理後的航空器運行軌跡數據ΛX和Λy視為隱馬爾科夫過程的顯觀測值，通過設定隱狀態 數目Ν'和參數更新時段ζ'，依據最近的Τ'個位置觀測值並採用B-W算法滾動獲取最新 隱馬爾科夫模型參數λ'：由於所獲得的航跡序列數據長度是動態變化的，為了實時跟蹤 航空器航跡的狀態變化，有必要在初始航跡隱馬爾科夫模型參數λ' = (π，Α，Β)的基礎上 對其重新調整，以便更精確地推測航空器在未來某時刻的位置。每隔時段ζ'，依據最新獲 得的Τ'個觀測值（〇1，〇2，…，〇Τ，）對航跡隱馬爾科夫模型參數λ' = (π，Α，Β)進行重新估 計。
[0078] 再而，依據隱馬爾科夫模型參數，採用Viterbi算法獲取當前時刻觀測值所對應 的隱狀態q。
[0079] 最後，每隔時段tP，根據最新獲得的隱馬爾科夫模型參數λ' = (π，A，B)和最近 H個歷史觀測值（〇1，O2，…，〇H)，基於航空器當前時刻的隱狀態q，通過設定預測時域h'，在 時刻t獲取航空器在未來時段h'的位置預測值0。
[0080] 所述聚類個數M'的值為4,隱狀態數目Ν'的值為3,參數更新時段ζ'為30秒， Τ'為10,預測時域h'為300秒，時段φ為4秒。
[0081] (3)實時飛行衝突監控與告警
[0082] 當系統有可能出現違反安全狀態集的狀態時，通過控制器實施狀態監控，對航空 器實施有效的管制措施，避免飛行衝突的發生。
[0083] 如圖3所示，其具體實施過程如下：
[0084] 首先，構造基於管制規則的衝突超曲面函數集。空中交通管制約束的違反都可以 視為被控對象（管制空域飛行的多架航空器）構成系統穿越超曲面而產生的事件，建立超 曲面函數集用以反映系統的衝突狀況。其中，衝突超曲面中與單一航空器相關的連續函數 /I1:Z4i?為第〗型超曲面，而將與兩架航空器相關的連續函數/?:ZxB為第π型超 曲面。
[0085] 然後，建立由航空器連續狀態至離散衝突狀態的觀測器。需要根據管制規範建立 觀測器，觀測系統系統穿越超曲面而產生的衝突事件，以便控制器做出相應的控制決策指 令。觀測器ξ用於觀測系統中航空器位置的連續變化而產生衝突事件，稱為第 I型觀測器，ξ2 ?為第II型觀測器。
[0086] 最後，設計從衝突到衝突解脫手段的離散監控器。當觀測器的離散觀測向量表明 某一非期望的狀態出現時，立刻發出相應的告警。該離散監控器可描述為函數βJ 其中S是觀測器觀測向量展成的空間，D是所有決策向量d展成的空間。
[0087] 顯然，上述實施例僅僅是為清楚地說明本發明所作的舉例，而並非是對本發明的 實施方式的限定。對於所屬領域的普通技術人員來說，在上述說明的基礎上還可以做出其 它不同形式的變化或變動。這裡無需也無法對所有的實施方式予以窮舉。而這些屬於本發 明的精神所引伸出的顯而易見的變化或變動仍處於本發明的保護範圍之中。
【權利要求】
1. 一種空中交通管制系統的航空器軌跡預測方法，所述空中交通管制系統包括機載終 端模塊、數據通信模塊、監視數據融合模塊以及管制終端模塊；監視數據融合模塊用於實現 空管雷達監視數據與自動相關監視數據的融合，為管制終端模塊提供實時航跡信息；其特 徵在於： 所述管制終端模塊包括以下子模塊： 飛行前無衝突4D航跡生成模塊，根據飛行計劃和世界區域預報系統的預報數據，建立 航空器動力學模型，然後依據飛行衝突耦合點建立航跡衝突預調配理論模型，生成航空器 無衝突4D航跡； 飛行中短期4D航跡生成模塊，依據監視數據融合模塊提供的實時航跡信息，利用隱馬 爾科夫模型，推測未來一定時間窗內的航空器4D軌跡； 實時飛行衝突監控與告警模塊，用於建立從航空器的連續動態到離散衝突邏輯的觀測 器，將空中交通系統的連續動態映射為離散觀測值表達的衝突狀態；當系統有可能違反空 中交通管制規則時，對空中交通混雜系統的混雜動態行為實施監控，為管制員提供及時的 告警信息；所述空中交通管制系統的飛行衝突預警方法包括如下幾個步驟： 步驟A、飛行前無衝突4D航跡生成模塊根據飛行計劃和世界區域預報系統的預報數 據，建立航空器動力學模型，並依據飛行衝突耦合點建立航跡衝突預調配理論模型，生成航 空器無衝突4D航跡； 步驟B、監視數據融合模塊將空管雷達監視數據與自動相關監視數據進行融合，生成航 空器實時航跡信息並提供給管制終端模塊；管制終端模塊中的飛行中短期4D航跡生成模 塊依據航空器實時航跡信息和歷史航跡信息推測未來一定時間窗內的航空器4D軌跡；所 述依據航空器實時航跡信息和歷史航跡信息推測未來一定時間窗內的航空器4D軌跡的具 體實施過程如下： 步驟B6、對航空器軌跡數據預處理，依據所獲取的航空器原始離散二維位置序列X= [X1,x2, . . .，xn]和y=Iiy1,y2, . . .，yn]，採用一階差分方法對其進行處理獲取新的航空器離 散位置序列AX= [AX1,Ax2, ? --，Axn_J和Ay= [Ay"Ay2, ? --，Ayn_J,其中 Axb=xb+1-xb,Ayb=yb+1_yb(b=I, 2, ? ? ? ,n-l); 步驟B7、對航空器軌跡數據聚類，對處理後新的航空器離散二維位置序列Ax和Ay， 通過設定聚類個數M'，採用遺傳聚類算法分別對其進行聚類； 步驟B8、對聚類後的航空器軌跡數據利用隱馬爾科夫模型進行參數訓練，通過將處理 後的航空器運行軌跡數據AX和Ay視為隱馬爾科夫過程的顯觀測值，通過設定隱狀態數 目N'和參數更新時段G'，依據最近的T'個位置觀測值並採用B-W算法滾動獲取最新隱馬 爾科夫模型參數A' ； 步驟B9、依據隱馬爾科夫模型參數，採用Viterbi算法獲取當前時刻觀測值所對應的 隱狀態q; 步驟B10、通過設定預測時域h'，基於航空器當前時刻的隱狀態q，獲取未來時段航空 器的位置預測值〇 ; 步驟C、實時飛行衝突監控與告警模塊建立從航空器的連續動態到離散衝突邏輯的觀 測器，將空中交通系統的連續動態映射為離散觀測值表達的衝突狀態；當系統有可能違反 空中交通管制規則時，對空中交通混雜系統的混雜動態行為實施監控，為管制員提供及時 的告警信息。
2. 根據權利要求1所述的一種空中交通管制系統的飛行衝突預警方法，其特徵在於： 步驟B中，所述聚類個數M'的值為4,隱狀態數目N'的值為3,參數更新時段C為30秒， T'為10,預測時域h'為300秒。
3. 根據權利要求1或2所述的一種空中交通管制系統的飛行衝突預警方法，其特徵在 於：步驟B的B8具體是指：由於所獲得的航跡序列數據長度是動態變化的，為了實時跟蹤 航空器航跡的狀態變化，有必要在初始航跡隱馬爾科夫模型參數X' =(IA，B)的基礎上 對其重新調整，以便更精確地推測航空器在未來某時刻的位置；每隔時段G'，依據最新獲 得的T'個觀測值（〇1，〇2,. . .，〇T，）對航跡隱馬爾科夫模型參數A' = (31，A，B)進行重新估 計； 步驟B的BlO具體是指：每隔時段？，根據最新獲得的隱馬爾科夫模型參數X' = (，A，B)和最近H個歷史觀測值（〇1，〇2, . . .，〇H)，基於航空器當前時刻的隱狀態q，通過設 定預測時域h'，在時刻t獲取航空器在未來時段h'的位置預測值0。
4. 根據權利要求1至3之一所述的一種空中交通管制系統的飛行衝突預警方法，其特 徵在於：時段9為4秒。
5. 根據權利要求1至4之一所述的一種空中交通管制系統的飛行衝突預警方法，其特 徵在於：所述步驟A的航空器無衝突4D航跡按照以下方法生成： 步驟A1、進行航空器狀態轉移建模，根據飛行計劃中航空器的飛行高度剖面，建立單個 航空器在不同航段轉移的Petri網模型：E= (g，G，Pre，Post，m)為航空器階段轉移模型， 其中g表示飛行航段，G表示垂直剖面中飛行狀態參數的轉換點，Pre和Post分別表示航段 和航路點的前後向連接關係，Z+表示航空器所處的飛行階段； 步驟A2、建立航空器全飛行剖面混雜系統模型如下， vH-KCAS；Mach,hp,tL0C), vgs-t1CAS)Mach,hp,tL0C,vffS,a), 其中Vffls為校IH空速，Mach為馬赫數，hp為氣壓高度，a為風向預報與航路的夾角，vws 為風速預報值，為溫度預報值，V高度變化率，V<^為地速； 步驟A3、採用混雜系統仿真的方式推測求解航跡：採用將時間細分的方法，利 用狀態連續變化的特性遞推求解任意時刻航空器在某一飛行階段距參考點的航程 卜人+ 和高度_) = &+JTb(T)?^其中Jci為初始時刻航空器距參考點 的航程，AT為時間窗的數值，J(T)為T時刻航空器距參考點的航程，Iltl為初始時刻航 空器距參考點的高度，h(T)為T時刻航空器距參考點的高度，由此可以推測得到單航空 器的4D航跡； 步驟A4、對多航空器耦合模型實施無衝突調配：根據兩航空器預達交叉點的時間，按 照空中交通管制原則，對交叉點附近不滿足間隔要求的航空器4D航跡進行二次規劃，得到 無衝突4D航跡。
6. 根據權利要求1至5之一所述的一種空中交通管制系統的飛行衝突預警方法，其特 徵在於：所述步驟B中監視數據融合模塊將空管雷達監視數據與自動相關監視數據進行融 合，生成航空器實時航跡信息，具體按照以下方法： 步驟Bl、將坐標單位和時間統一； 步驟B2、採用最鄰近數據關聯算法將屬於同一個目標的點相關聯，提取目標航跡； 步驟B3、將分別從自動相關監視系統和空管雷達提取的航跡數據從不同的時空參考坐 標系統變換、對準到管制終端統一的時空參考坐標系統； 步驟B4、計算兩條航跡的相關係數，若相關係數小於某一預設閾值，則認為兩條航跡不 相關；否則該兩條航跡相關，可以進行融合； 步驟B5、對相關的航跡進行融合。
7. 根據權利要求1至6之一所述的一種空中交通管制系統的飛行衝突預警方法，其特 徵在於：所述步驟B5中對相關的航跡進行融合，採用基於採樣周期的加權平均算法，其加 權係數根據採樣周期和信息精度確定，再利用加權平均算法將與之相關的自動相關監視航 跡和空管雷達航跡融合為系統航跡。
8. 根據權利要求1至7之一所述的一種空中交通管制系統的飛行衝突預警方法，其特 徵在於：所述步驟C的具體實施過程如下： 步驟C1、構造基於管制規則的衝突超曲面函數集：建立超曲面函數集用以反映系統的 衝突狀況，其中，衝突超曲面中與單一航空器相關的連續函數hB&為第I型超曲面， 與兩架航空器相關的連續函數心^〃為第II型超曲面； 步驟C2、建立由航空器連續狀態至離散衝突狀態的觀測器：需要根據管制規範建立觀 測器，觀測系統系統穿越超曲面而產生的衝突事件，以便控制器做出相應的控制決策指令； 觀測器I用於觀測系統中航空器位置的連續變化而產生衝突事件，稱S為第I型 觀測器，^ : ZBS為第II型觀測器； 步驟C3、設計從衝突到衝突解脫手段的離散監控器，該離散監控器可描述為函數pD，其中S是觀測器觀測向量展成的空間，D是所有決策向量d展成的空間；當觀測 器的離散觀測向量表明某一非期望的狀態出現時，立刻發出相應的告警。
【文檔編號】G08G5/04GK104504942SQ201510007919
【公開日】2015年4月8日 申請日期:2015年1月7日 優先權日:2015年1月7日
【發明者】韓雲祥, 趙景波, 李廣軍 申請人:江蘇理工學院