張鈺沁 潘衛(wèi)軍 姜沿強
(中國民用航空飛行學院空中交通管理學院 廣漢 618307)
近年來我國大力新建中小型支線機場,中小型機場是我國航空業(yè)發(fā)展的重要基礎設施。根據(jù)2019 年民航局的統(tǒng)計公報,我國年旅客吞吐量200萬人次以下的中小型機場有165 個,占中國機場總數(shù)量的70%,占全國機場吞吐量的6.8%,從地區(qū)來看,現(xiàn)西南地區(qū)中小型機場占比最大,且目前西南地區(qū)部分支線機場地理距離較近,分布密集,形成了多個支線機場群。例如騰沖、保山、芒市機場群,滄源、瀾滄、西雙版納、思茅機場群等。由于各支線機場因地理位置架構的機場群周邊空域沖突較大,管制的協(xié)調(diào)量增大,管制效率降低,存在一定的安全隱患。
對于存在沖突空域的中小機場群,空中交通服務能力的下降會引發(fā)容需不平衡,不可避免會造成延誤的發(fā)生,需要建立一種機場管制聯(lián)合運行策略對現(xiàn)已有的中小機場進行聯(lián)合網(wǎng)絡節(jié)點調(diào)度。本文提出一種適用于執(zhí)行西南地區(qū)支線機場航班的調(diào)度方案:通過建立機場群聯(lián)合運行的優(yōu)勢,可通過協(xié)調(diào)運行的方式,計算每架航空器運行過程中的累積延誤,并對累積延誤高的航班進行合理的優(yōu)先排序的調(diào)度方案,分攤航班總延誤,有助于在保持吞吐量的同時盡可能減少平均延誤,也可以利用減少累積延誤來避免在空中和地面的過度延誤,從而減少排隊阻塞造成的級聯(lián)擁塞效應。
航班延誤通常定義為計劃起飛/到達時間與實際起飛/到達時間的差值[1],一般將15min 作為劃分航班延誤的依據(jù),作為衡量交通系統(tǒng)運行性能的關鍵指標,將航班晚于計劃15min 起飛或到達視為航班起飛/到達延誤[2]。造成延誤的原因分為幾類:天氣、航空公司、空管等。對于天氣原因:一方面機場容量易受天氣影響,惡劣的氣象條件會造成機場容量下降,導致容需不平衡,引發(fā)航班延誤[3~4]。另一方面,航路上的惡劣氣象條件可能會使飛機繞飛,導致總飛行時間(gate-to-gate time)延長[5]。對于航空公司的原因:航空公司通常在排班時會為航班在預定的時間內(nèi)增加一些填充的延誤時間,以便在預期的延誤情況下保持預期的準點率,因此航班時刻表的時刻并不代表無延誤的行程。對于空管的原因:當空域飛行量超過區(qū)域或終端區(qū)扇區(qū)保障能力時,為保障空中交通的整體流動性,會實施空中交通管制措施[6];地面等待策略(Ground Delay Programs,GDP)和地面保障資源的限制也會造成了部分航班的延誤[7~8]。GDP 是指當航班到達機場容量下降后,不能在指定時段為航班提高到達服務時,通常在出發(fā)機場實施地面等待,從而將空中延誤轉(zhuǎn)化為地面延誤,以減小延誤損失。除上述因素外,航空公司的資源管理問題也會導致出現(xiàn)延誤。
當一架航空器執(zhí)行飛行任務前,通常唯一可參考的是航空公司的排班時刻表,但時間表包含了航空公司對該航班的填充延誤,且忽略了該航空器在之前飛行任務中造成的延誤[9]。每階段航程所積累的延誤會由于航班之間的連接性的特點傳向下游航班,造成持續(xù)的延誤傳播[9]。
本文提出了累積延誤的概念,不同于上文中的延誤的定義,是基于未阻礙的預期時刻UT而不是排班時刻為基準計算的,無障礙時刻是基于歷史統(tǒng)計數(shù)據(jù)或軌跡生成的暢通時刻。以無障礙預期時刻為參考基準,更能挖掘每段行程的潛在延誤,便于通過算法對航班進行合理的調(diào)配[11~12]。累積延誤可用于度量執(zhí)行連續(xù)飛行任務的航空器在當前時刻積累的延誤量。因此要對每架航空器的延誤進行連續(xù)監(jiān)測計量,探測連續(xù)的航班延誤狀態(tài)有助于整合連續(xù)的延誤決策。這種以航空器為導向的視角是對“gate to gate”視角的補充[13]。這種觀點對于新型的高密度、短途飛行模式尤其重要。通常的飛行會依據(jù)FCFS 的算法分配航班順序[10],而累積延誤指標可用于優(yōu)化分配算法,重新制定航空器的起降優(yōu)先級。
后文將分別對基于無阻礙暢行時刻、可行時刻、累積延誤三種航空器的調(diào)度參數(shù)進行對比分析。
1)使用以下的符號來跟蹤表示不同類型的時刻,且都表示為時刻點而不是時間段,給每個符號加上“D”來表示出發(fā)時刻,加上“A”來表示到達時刻,例如“FTD”表示可行的出發(fā)時刻,“UTA”表示無阻礙的到達時刻。
(1)UTA/D(Unimpeded Times):未受阻礙的預期進場/離場時刻。用于表示某航空器在任何階段都不產(chǎn)生延誤,僅根據(jù)航行軌跡生成的暢通時間為準。
(2)FTA/D(Feasible Times):進場/離場的最早可行時刻。用于表示某航空器基于離開上一節(jié)點/離開該節(jié)點的時刻所更新的預計到達該節(jié)點/到達下一節(jié)點的時刻。
(3)STA/D(Schedule Times):依據(jù)調(diào)度程序時間間隔標準分配的連續(xù)進場/離場時刻。
(4)ATA/D(Actual Times):實際到場/離場時刻。
2)使用以下符號來跟蹤表示不同類型的時間,且都表示時間段而不是時刻。
(1)st:連續(xù)進離場時間間隔。
(2)ut1(unimpeded travel time):無阻礙行程時間。
(3)ut2(unimpeded turn time):無阻礙的周轉(zhuǎn)時間。
2.3.1 基于無障礙暢行時刻(UT)進行調(diào)度
通過將無阻礙的航空器行程時與周轉(zhuǎn)時間相加,可根據(jù)起飛時間推算出每架航空器處于每個節(jié)點時的UTs。UTs在模擬的過程中是固定的,因為航班運行存在著不同程度的延誤狀況,該方法無法用于實際航班運行調(diào)度,但可以作為無阻礙時間的標準用于測量航空器的累積延誤。
2.3.2 基于可行時刻(FT)對航班進行調(diào)度
基于航空器在前一節(jié)點累積延誤的可行期望時刻進行調(diào)度的方法。期望時刻是通過計算同一航空器在其行程的前一個節(jié)點上所產(chǎn)生的延誤來估計的到達下一個節(jié)點的時刻。其中FTs被假定由航空器更新,并與節(jié)點調(diào)度程序共享。在每個節(jié)點更新航空器的狀態(tài),計算FT值如下:
使當前節(jié)點為i,航空器離開前一個節(jié)點i-1(FTA)或者到達當前節(jié)點i(FTD):
其中,ut1為無阻礙行程時間,ut2無阻礙的周轉(zhuǎn)時間。
基于可行時刻的調(diào)度需要對每架航空器行駛過程中到達目標節(jié)點的可行時刻進行估算并不斷更新,以反映任何階段產(chǎn)生的延誤。網(wǎng)絡調(diào)度程序會估計每架在任意節(jié)點航空器的FTs,應用最小時間間隔程序重新分配進離場時刻STs,并按照FCFS 原則生成航班時刻表,該方法具有較高的可行性,但其更新頻率會在一定程度上影響航班的吞吐量。
基于累積延誤算法調(diào)度程序是將基于可行時刻的調(diào)度方法進行優(yōu)化,加入了累積延誤對航班排序的影響,累積延誤稱為AD。首先使用FTs調(diào)度程序按照FCFS原則不斷更新航班時刻表,利用更新后的航班時刻表,利用調(diào)度算法在不造成額外延誤的情況下,安排AD高的航班優(yōu)先起飛或著陸,使低延誤航班承擔部分高延誤航班的累積延誤。在航班的交換過程中,總體延誤不變,僅減小航班的延誤平均方差。
圖1 航班調(diào)度算法流程圖
1)累積延誤的計算公式為其中,ST是基于FT分配的到該節(jié)點的時刻,UT為該航空器在無阻礙情況下該節(jié)點的時刻。
因為FT只能反映出當前節(jié)點所產(chǎn)生的延誤,對于還未到達目標節(jié)點的航班,累積延誤量取決于ST與UT的差值,對于已到達目標節(jié)點的航班,累積延誤量取決于AT與UT的差值。
2)根據(jù)計算的AD進行排序。
對于可能存在沖突調(diào)度的航班,對AD高的航班,在可行范圍內(nèi)為了不產(chǎn)生新的延誤,只在AD方差減小的情況下交換航班順序,并為交換后的航班分配新的ST與AD,AD能夠直觀地反映調(diào)度后航班延誤的變化情況。該種調(diào)度方法只減小AD的方差,其均值在交換后保持不變,因為減少高延誤航班的減少與低延誤航班的增加是完全匹配的,目的在于緩解高延誤,平均航班的總體延誤。
在僅僅在出現(xiàn)累積延誤方差減小的情況下對現(xiàn)已有的航班次序交換后,可根據(jù)st計算出每架航班交換后的新的航班時刻表的ST,其算法如下:
首先將各航空器交換后產(chǎn)生的新FT從小到大排列得:對于FT1……FTi,
其中,ST為依據(jù)調(diào)度程序時間間隔標準分配的連續(xù)進場/離場時刻,F(xiàn)T為場/離場的最早可行時刻,st為過渡時間。
建立該仿真模型目的是為了模擬仿真航空器調(diào)度方案對于實際運行的可操作性和緩解航班累積延誤的有效性。在本節(jié)中,將以云南地區(qū)各機場群為例,使用節(jié)點來代替描述機場,節(jié)點網(wǎng)絡調(diào)度程序描述機場群的聯(lián)合運行模式,對機場群進行網(wǎng)絡調(diào)度進行模擬仿真。
跟蹤執(zhí)行中小機場群飛行任務的B737 航空器的實際飛行數(shù)據(jù),為模擬仿真場景提供數(shù)據(jù)參考。主要是以昆明長水機場為樞紐機場,來往于往香格里拉迪慶、麗江三義、保山云瑞、普洱思茅等中小機場群。據(jù)此建立一個以主要節(jié)點A,四個次要節(jié)點B、C、D、E 與四個支線節(jié)點B1、C1、D1、E1所組成的仿真模型。節(jié)點A代表樞紐機場,其余節(jié)點代表中小機場群,設置一組航空器從分別從節(jié)點A 出發(fā),前往B、C、D、E 節(jié)點。節(jié)點之間的每段行程距離不等,參考實際距離取值,且在模擬過程中保持不變。每架航空器的速率參考B737飛機的平均航行速度進行仿真。航空器在機場內(nèi)變更飛行任務時的周轉(zhuǎn)時間在40min~45min范圍內(nèi)取值。
圖2 場景仿真模型
在模擬過程中,四架航空器分別從節(jié)點A 出發(fā),沿著各條航線前往次要節(jié)點與支線節(jié)點,再執(zhí)行返航任務。其中,每段航程的UTs為固定值累加,在模擬的行程過程中產(chǎn)生的延誤隨機。FTs取值根據(jù)航班的隨機延誤情況進行更新,并生成對應的STs值。各架航空器每到達一個節(jié)點后,基于機場群的網(wǎng)絡調(diào)度程序(機場聯(lián)合運行程序),通過計算模擬仿真出航班實際運行中的累積延誤,輸出其對應的UT、FT、ST、AD值。
通過進行完整的場景模擬仿真,輸出每段航程的UT、FT、ST、AD值。經(jīng)過模擬仿真得出多組用于調(diào)配的航班數(shù)據(jù),觀察大量仿真數(shù)據(jù)可得知,對于時刻表上可以進行調(diào)度的延誤航班,還需要具備一定的條件,才能進行調(diào)度。否則兩架航空器在交換航班次序后,不但不能減小平均延誤,會使其累積延誤的方差會增大。
通過應用上文提出的算法,可以篩除無需進行調(diào)度的航班數(shù)據(jù),僅僅留存交換后累積延誤方差減小的航班進行調(diào)度,表1 為其中一組趨于平均的仿真數(shù)據(jù)。
表1 航程輸出數(shù)據(jù)表
由表1 的數(shù)據(jù)得出,當五架航空器匯聚于A 機場時,按照FT先到先服務的對航空器的進場時間的排序為b、a、e、c、d。
由表1 可得知航空器c 的累積延誤最高,且根據(jù)排序規(guī)則,由于航班e 的ST小于航班c 的FT,表明航空器c與航空器e進行交換后不會產(chǎn)生新的延誤,且模擬交換后,所得到的AD方差明顯減小,如表2所示。
表2 調(diào)度數(shù)據(jù)更新表
最終航空器調(diào)度后按順序從小到大排序為航空器b、a、c、e、d。
根據(jù)仿真模擬出場景,對累積延誤調(diào)度的結果進行優(yōu)化分析。在進一步的模擬中,對FT 采用不同的更新頻率,采用每分鐘更新一次航空器的FT,與5min更新一次進行對比。
圖3 FT不同更新頻率下的累積延誤
仿真結果表明對于FT 不同更新頻率下,航空器所產(chǎn)生的累積延誤不同,更新頻率越快,航空器的累積延誤越低,表明通過提高航空器的時刻位置信息實時更新頻率可以有效減少航班的累積延誤,同樣有助于降低航班的整體延誤水平。
通過探究網(wǎng)絡節(jié)點的容量是否會影響調(diào)度的效果。將處于某節(jié)點的航空器數(shù)量設定為固定值,其余條件均不變,分別模擬處于不同節(jié)點容量時,累積延誤調(diào)度方式的效果。
圖4 不同節(jié)點容量網(wǎng)絡調(diào)度效果圖
從模擬出的網(wǎng)絡調(diào)度效果圖來看,對于不同的節(jié)點容量,基于累積延誤的調(diào)度方法同樣能夠有效地減小延誤的方差值。
通過對歷史數(shù)據(jù)的分析發(fā)現(xiàn),同一架航空器從一個航班到另一個航班的傳播延誤運輸是航空公司延誤的主要來源,而航空公司的排班時刻表填充隱藏了這種延誤運輸效應。本文針對西南地區(qū)中小機場群的運行模擬仿真,研究了利用累積延誤度量節(jié)點進行網(wǎng)絡調(diào)度的方法,減輕了該延誤在航空運輸中的累積效應。使用無阻礙預期時刻來代替航空公司的排班表作為測量延誤的基準,并根據(jù)行程中造成的延誤來更新預期的可行時刻,從而提高節(jié)點調(diào)度的靈活性和可行性。模擬仿真結果證明調(diào)度后的平均延誤基本保持不變,僅減小調(diào)度航空器之間累積延誤的方差,能夠縮短高延誤航空器的累積延誤,提高航班的整體運行效率。