陳 楠 孔建國

（中國民航飛行學院 交通運輸管理學院，四川 德陽618399）

我國民航業(yè)發(fā)展正從規(guī)模速度向質(zhì)量效率轉(zhuǎn)變，空域資源仍舊緊缺，空域扇區(qū)規(guī)劃不合理，空中交通擁堵問題依舊嚴重，2019 年平均航班正常率為81.43%[1]。某個空域扇區(qū)容量下降或者失效可能會導致大范圍的空域擁堵，進而造成嚴重的航班延誤。因此，從空域扇區(qū)入手，構(gòu)造空域扇區(qū)網(wǎng)絡，利用復雜網(wǎng)絡理論，尋找導致空域擁堵的傳播方式是十分有必要的[2]。復雜網(wǎng)絡已經(jīng)在民航領域被廣泛使用，具有完整的體系和豐富的理論分析。Amaral 最早利用復雜網(wǎng)絡來解決航空領域的問題，分析一定范圍內(nèi)的機場網(wǎng)絡特征，發(fā)現(xiàn)其符合小世界網(wǎng)絡的模型特征[3]。中國學者蔡開泉將導航臺或航向信標臺作為節(jié)點，航路段為邊建立中國航路網(wǎng)絡，分析網(wǎng)絡特性[4]。學者們利用復雜網(wǎng)絡理論研究的對象主要集中于機場網(wǎng)絡和航路網(wǎng)絡，較少從空中管制員角度出發(fā)，從系統(tǒng)層面研究空域扇區(qū)網(wǎng)絡的結(jié)構(gòu)特性及其抗毀性。本文結(jié)合管制員的實際工作定義扇區(qū)網(wǎng)絡的統(tǒng)計特征指標，利用復雜網(wǎng)絡理論分析其拓撲特性，以華北管制扇區(qū)為例分析其抗毀特性。

1 扇區(qū)網(wǎng)絡建模及特征參數(shù)

表1 特征指標定義

圖1 華北管制空域扇區(qū)結(jié)構(gòu)圖

圖2 華北管制空域扇區(qū)網(wǎng)絡圖

2 扇區(qū)網(wǎng)絡實證與特性分析

華北空管局負責我國北京、天津、河北、山西、內(nèi)蒙古五個地區(qū)的民航事務管理，隨著北京大興機場的投入使用，華北區(qū)域管制中心所轄空域?qū)⒏臃泵?。如圖1 所示，為華北區(qū)域管制中心的結(jié)構(gòu)圖。我們?nèi)?019年5 月某日高峰時段6 小時的飛行航跡數(shù)據(jù)，通過航跡重演，得到各扇區(qū)間航班流量情況。以空域各扇區(qū)中心位置為點集，扇區(qū)間的航班流為邊集，建立華北空域管制扇區(qū)網(wǎng)絡圖，如圖2 所示。

各扇區(qū)間的靜態(tài)統(tǒng)計特征指標在一定程度上可以看出各個扇區(qū)的繁忙程度、扇區(qū)的重要性以及扇區(qū)網(wǎng)絡的拓撲結(jié)構(gòu)等，由各統(tǒng)計特征參數(shù)的定義，編程計算得到華北管制區(qū)23 個扇區(qū)的具體參數(shù)數(shù)據(jù)。

其中，度值最大的為北京01 扇，位于華北管制空域的中心地理位置，周圍與其有直接航班聯(lián)系的扇區(qū)眾多，扇區(qū)繁忙。強度最大的扇區(qū)為北京02 扇，該扇區(qū)是終端扇區(qū)，覆蓋最繁忙的機場——北京首都機場，航班流量極大，高峰時段的航班流量超過1200 架次，管制員工作負荷大。北京01 扇區(qū)同樣是介數(shù)最大的扇區(qū)，是多條最短路徑都匯集處，具有較強的中樞性，對扇區(qū)網(wǎng)絡的整體運行狀況有重要作用，在有關扇區(qū)抗毀性研究時應特別關注此扇區(qū)。聚集系數(shù)最大的扇區(qū)為北京07 扇、北京12扇、北京14 扇、北京17 扇，周邊的扇區(qū)大都有直接航班聯(lián)系，空域管理者間溝通密切，便于協(xié)同處理空域擁堵問題。

3 空域扇區(qū)網(wǎng)絡抗毀性研究

為評估扇區(qū)網(wǎng)絡對隨機攻擊和選擇性攻擊的抗毀性，尋找扇區(qū)網(wǎng)絡中的關鍵扇區(qū)，在此選用網(wǎng)絡效率和最大連通子圖相對值作為測度指標。首先，選擇隨機連續(xù)攻擊，用Matlab 依順序移除華北管制扇區(qū)網(wǎng)絡相應序號扇區(qū)。開始移除扇區(qū)較少時，兩個測度指標隨著移除的數(shù)目緩慢下降，整個網(wǎng)絡能夠維持原有特性，正常運轉(zhuǎn)。當移除到第9 個扇區(qū)時，兩個測度指標都急劇下降，整個空域扇區(qū)網(wǎng)絡的不能維持原有的功能。之后網(wǎng)絡效率下降到0.17 以下，網(wǎng)絡最大連通子圖中的扇區(qū)數(shù)下降到7個，扇區(qū)網(wǎng)絡完全崩潰。接著采用選擇性攻擊來探究扇區(qū)網(wǎng)絡的抗毀性，在此又分為度優(yōu)先攻擊和介數(shù)優(yōu)先攻擊。將網(wǎng)絡中各扇區(qū)分別按照度值和介數(shù)降序排列，分別攻擊度值和介數(shù)大的扇區(qū)。抗毀性測度指標兩種優(yōu)先攻擊策略下，抗毀性測度指標從一開始就急劇下降，網(wǎng)絡提前進入崩潰狀態(tài)。介數(shù)優(yōu)先攻擊策略比度值優(yōu)先攻擊，測度指標下降的更快，對整個扇區(qū)的破壞性更大，從而說明介數(shù)是扇區(qū)重要性排序最關鍵的指標，測度指標對整個空域網(wǎng)絡的影響最大。

4 關鍵扇區(qū)的影響力

由于介數(shù)是扇區(qū)重要性排序的關鍵指標，選擇攻擊介數(shù)值最大的扇區(qū)，分析其對其余扇區(qū)以及整個空域扇區(qū)網(wǎng)絡的影響是研究空域扇區(qū)網(wǎng)絡抗毀性的有效方法。在此，將扇區(qū)按介數(shù)降序排序，選擇最大的4 個扇區(qū)：北京01 扇、北京04 扇、北京03 扇、太原02 扇。按順序移除相應的扇區(qū)，并計算其它扇區(qū)的網(wǎng)絡效率E 以及介數(shù)的變化率。尋找關鍵扇區(qū)失效后，影響最大的扇區(qū)即為潛在關鍵扇區(qū)[7]。當北京03 扇和太原02 扇失效后，北京05 扇和太原01 扇的介數(shù)變化率激增，關鍵扇區(qū)內(nèi)的負載外流，導致這些潛在關鍵扇區(qū)內(nèi)的航班量激增，管制員需要提前采取措施疏散這些扇區(qū)內(nèi)的航班，如果沒有提前控制航班流量，將會造成這些扇區(qū)的擁堵，并產(chǎn)生雪崩效應，造成更多扇區(qū)擁堵，出現(xiàn)“級聯(lián)失效”（表2）。

表2 攻擊關鍵扇區(qū)的影響力

5 結(jié)論

本文通過對華北空域扇區(qū)網(wǎng)絡的實證分析，利用復雜網(wǎng)絡理論及各網(wǎng)絡拓撲特性指標的定義，計算出各個扇區(qū)的度、強度、介數(shù)、平均路徑長度等值。通過兩大抗毀性測度指標以及兩種攻擊策略，對華北空域扇區(qū)網(wǎng)絡進行攻擊，得到對扇區(qū)抗毀性影響最大的指標——介數(shù)，從而找到網(wǎng)絡中的關鍵扇區(qū)。之后，繼續(xù)對關鍵扇區(qū)進行選擇性破壞，尋找受影響最大的扇區(qū)。在關鍵扇區(qū)發(fā)生擁堵或失效時，需要對這些扇區(qū)采取保護措施，防止扇區(qū)擁堵進一步擴散。本文首次從空中管制員的視角，構(gòu)建扇區(qū)網(wǎng)絡模型，尋找空域擁堵的傳播途徑和關鍵扇區(qū)，對扇區(qū)網(wǎng)絡的抗毀性進行實證分析，之后的研究工作可以考慮擴大研究對象的范圍以及對扇區(qū)網(wǎng)絡的動態(tài)抗毀性進行深入探索。