韋 薇,何 晶
(上海工程技術(shù)大學(xué)航空運(yùn)輸學(xué)院,上海 201620)
在民航業(yè)高速發(fā)展的背景下我國旅客運(yùn)輸量不斷增加,航班資源供給與需求之間存在的矛盾逐漸凸顯出來。多機(jī)場區(qū)域在國家政策和市場需求的推動下逐漸形成。但也存在一些問題,包括差異化發(fā)展不顯著、資源分配不均勻、機(jī)場利用率低和機(jī)場過度飽和等。在生產(chǎn)運(yùn)營活動中航空公司的核心和基礎(chǔ)是航班計(jì)劃,航空公司的市場競爭力和經(jīng)濟(jì)狀況都會受到航班計(jì)劃的影響。航班頻率是航班計(jì)劃中的重要環(huán)節(jié),因此需要對多機(jī)場群航班頻率進(jìn)行優(yōu)化。
樂美龍等人考慮票價(jià)、旅行時(shí)間、延誤時(shí)間和過站時(shí)間等因素構(gòu)建旅行負(fù)效用函數(shù),并在多行程路徑的基礎(chǔ)上將航空公司收益最大化作為目標(biāo),構(gòu)建目標(biāo)函數(shù),通過旅行負(fù)效用函數(shù)和目標(biāo)函數(shù)構(gòu)成多機(jī)場群航班頻率優(yōu)化模型,完成航班頻率的優(yōu)化,該方法沒有考慮旅客的出行成本,導(dǎo)致方法存在出行成本高的問題。姜思露等人在多機(jī)型多航線競爭環(huán)境下將航班運(yùn)營成本和旅客計(jì)劃延誤成本最小作為目標(biāo)函數(shù),通過建模優(yōu)化方法構(gòu)建航班頻率優(yōu)化模型,該方法優(yōu)化后存在延誤時(shí)間長和延誤航班數(shù)量多的問題。張兆寧等人將相關(guān)運(yùn)行模式下的優(yōu)化目標(biāo)選取最小航班延誤成本,構(gòu)建航班頻率優(yōu)化模型,并通過滾動時(shí)間窗啟示算法進(jìn)行求解,完成航班頻率優(yōu)化,該方法優(yōu)化后的飛機(jī)日利用率較低。
為了解決上述方法中存在的問題,提出考慮出行成本的多機(jī)場群航班頻率優(yōu)化模型。
通過航班頻率的合理制定,可以提高交通資源的利用率,在出行過程中指導(dǎo)游客作出選擇行為,航班頻率的制定會參考游客的出行選擇行為。通過對其進(jìn)行博弈分析,建立多機(jī)場群航班頻率優(yōu)化模型。
R
代表的是出行方式對應(yīng)的社會收益;f
代表的是出行方式j
對應(yīng)的航班頻率,此時(shí)存在下式maxR
=q
(V
-V
+p
)(1)
式中,q
代表的是客流量;V
代表的是最大出行費(fèi)用;p
代表的是第j
種出行方式對應(yīng)的票價(jià);V
代表的是旅客廣義費(fèi)用。通過上述公式可知,q
、V
受航班頻率f
的影響。出行方式在上層規(guī)劃過程中對應(yīng)的社會效應(yīng)如下(2)
f
-描述除了第j
種出行方式外,旅客選擇其他出行方式的頻率,為了最大化社會效益,需要根據(jù)f
-制定出行方式j
的頻率f
,對出行方式航班頻率進(jìn)行博弈分析,構(gòu)建多機(jī)場群航班頻率的上層優(yōu)化模型(3)
式中,Ω
受f
-的影響,是由出行方式構(gòu)成的可行決策集合。用K
描述不同出行方式下對應(yīng)的服務(wù)水平,可通過下式計(jì)算得到K
=U
+αf
(4)
式中,α
為針對航班頻率,旅客對應(yīng)的敏感程度;U
代表的是除航班頻率外,第j
種出行方式的地面交通時(shí)間、準(zhǔn)點(diǎn)率以及出行時(shí)間等綜合因素對應(yīng)的效用值。設(shè)X
代表的是游客選擇出行方式的概率,可通過Logit
分類模型計(jì)算得到,出行方式j
對應(yīng)的旅客流量,可根據(jù)概率X
計(jì)算得到,即q
=X
Q
,其中,Q
描述的是客流量總數(shù)。根據(jù)上述分析可將社會效益R
轉(zhuǎn)化為下式(5)
頻率f
可通過出行方式j
對應(yīng)的流量計(jì)算得到(6)
式中,X
-代表的是除了出行方式j
,其他出行方式對應(yīng)的選擇概率數(shù)值。在目標(biāo)函數(shù)中帶入上式,獲得關(guān)于頻率出行方式j
社會效益對應(yīng)的函數(shù),通過下述公式求其一階導(dǎo)數(shù)和二階導(dǎo)數(shù)(7)
通過上述分析可知,針對q
,社會效益屬于嚴(yán)格凹函數(shù)。設(shè)b
代表的是出行方式j
對應(yīng)的策略,b
-代表的是除了出行方式j
,其余出行方式的策略,因此最大化社會效益對應(yīng)的數(shù)據(jù)模型如下(8)
c
(q
)(9)
式中,δ
代表的是準(zhǔn)點(diǎn)率因素。可通過下式描述旅客選擇模型Z
(q
)(10)
式中,w
代表的是正定矩陣。考慮出行成本的多機(jī)場群航班頻率優(yōu)化模型構(gòu)建方法在Stackelberg博弈模型的基礎(chǔ)上構(gòu)建多機(jī)場群航班頻率優(yōu)化模型。
上層博弈模型為各種出行方式的最大化社會效益
(11)
最小化旅客廣義出行費(fèi)用為下層博弈模型
q
(f
)=arg minZ
(q
,f
)(12)
考慮出行成本的多機(jī)場群航班頻率優(yōu)化模型構(gòu)建方法采用遺傳算法對多機(jī)場群航班頻率優(yōu)化模型進(jìn)行求解,具體過程如下:
M
條跑道和N
架飛機(jī),則可通過以下兩個(gè)部分描述染色體I
(13)
式中,A
=(a
,a
,…,a
)描述的是需要著陸的飛機(jī)構(gòu)成的排列,即飛機(jī)鏈表;R
=(r
,r
,…,r
)代表的是飛機(jī)降落的序列;其中r
∈{1,2,…,N
}代表的是各飛機(jī)在鏈表中對應(yīng)的跑道,即跑道向量。A
中存在的基因?qū)儆陲w機(jī)基因;R
中存在的基屬于跑道基因。P
描述飛機(jī)在相同跑道中降落的順序,用U
={(i
,j
):ipj
}描述飛機(jī)序列構(gòu)成的集合,如果第一架在跑道中降落的是飛機(jī)i
,則用E
描述飛機(jī)i
的降落時(shí)間,否則重新對飛機(jī)i
降落的時(shí)間進(jìn)行計(jì)算。t
,用d
(I
,I
)描述個(gè)體I
和I
在集合HQS中的距離,其計(jì)算公式如下(14)
I
描述交叉操作中的母體;用I
描述交叉操作中的父體,其表達(dá)式分別如下(15)
考慮出行成本的多機(jī)場群航班頻率優(yōu)化模型構(gòu)建方法通過以下兩種變異方式對染色體進(jìn)行變異處理:
1)變異操作1
根據(jù)變異概率染色體I
對應(yīng)的跑道向量和飛機(jī)鏈表中存在的基因依次進(jìn)行變異操作,隨機(jī)獲取整數(shù)i
,j
∈[1,2,…,N
],并將i
、j
位置上的基因在跑道向量和飛機(jī)鏈表中進(jìn)行交換處理,獲得變異處理后的染色體。2)變異操作2
設(shè)置迭代終止條件,滿足條件時(shí)輸出多機(jī)場群航班頻率優(yōu)化模型的最優(yōu)解,完成多機(jī)場群航班頻率的優(yōu)化。
考慮出行成本的多機(jī)場群航班頻率優(yōu)化模型構(gòu)建方法采用遺傳算法求解模型的流程圖如圖1所示。
圖1 模型求解流程
為了驗(yàn)證考慮出行成本的多機(jī)場群航班頻率優(yōu)化模型構(gòu)建方法的整體有效性,需要對考慮出行成本的多機(jī)場群航班頻率優(yōu)化模型構(gòu)建方法進(jìn)行測試。
分別采用考慮出行成本的多機(jī)場群航班頻率優(yōu)化模型構(gòu)建方法、文獻(xiàn)[3]方法和文獻(xiàn)[4]方法進(jìn)行如下測試。
將計(jì)劃延誤時(shí)間作為測試指標(biāo),對所提方法、文獻(xiàn)[3]方法和文獻(xiàn)[4]方法進(jìn)行測試,測試結(jié)果如圖2所示。
圖2 不同方法優(yōu)化后的計(jì)劃延誤時(shí)間
根據(jù)圖2中的數(shù)據(jù)可知,采用所提方法對多機(jī)場群航班頻率進(jìn)行優(yōu)化后,與優(yōu)化前的計(jì)劃延誤時(shí)間相比,優(yōu)化后的計(jì)劃延誤時(shí)間明顯降低;采用文獻(xiàn)[3]方法對多機(jī)場群航班頻率進(jìn)行優(yōu)化后,公司1、公司2和公司3的計(jì)劃延誤時(shí)間有所降低,但對公司4和公司5的優(yōu)化效果不明顯;采用文獻(xiàn)[4]方法對多機(jī)場群航班頻率進(jìn)行優(yōu)化后,計(jì)劃延誤時(shí)間沒有明顯得到降低。對比上述方法的測試結(jié)果可知,所提方法優(yōu)化后的計(jì)劃延誤時(shí)間得到了明顯的改善,驗(yàn)證了所提方法的有效性。
圖3為所提方法、文獻(xiàn)[3]方法和文獻(xiàn)[4]方法優(yōu)化后的航班延誤數(shù)量。
圖3 不同方法的航班延誤數(shù)量
根據(jù)圖3中的數(shù)據(jù)可知,對不同公司的航班頻率優(yōu)化后,所提方法的航班延誤數(shù)量均在5架以內(nèi),文獻(xiàn)[4]方法的航班延誤數(shù)量在10架左右波動,文獻(xiàn)[3]方法的航班延誤數(shù)量最高,在15架左右。對比所提方法、文獻(xiàn)[3]方法和文獻(xiàn)[4]方法的測試結(jié)果可知,所提方法優(yōu)化后的延誤航班最少。
將出行成本作為指標(biāo),對所提方法、文獻(xiàn)[3]方法和文獻(xiàn)[4]方法進(jìn)行測試,測試結(jié)果如圖4所示。
圖4 不同方法優(yōu)化后的出行成本
對圖4進(jìn)行分析可知,所提方法優(yōu)化后的出行成本明顯降低,因?yàn)樗岱椒ㄔ跇?gòu)建多機(jī)場群航班頻率優(yōu)化的雙層模型時(shí),在下層博弈模型中考慮了旅客的出行成本,使其優(yōu)化后的出行成本得到了明顯的降低,文獻(xiàn)[3]方法和文獻(xiàn)[4]方法的出行成本在優(yōu)化后雖然有所降低,但與所提方法相比,出行成本的降低幅度較低,因?yàn)橐陨蟽煞N方法在構(gòu)建航班頻率優(yōu)化模型時(shí)沒有考慮旅客的出行成本,出行成本優(yōu)化效果較差。
將日利用率作為測試指標(biāo),驗(yàn)證所提方法、文獻(xiàn)[3]方法和文獻(xiàn)[4]方法的有效性,測試結(jié)果如表1所示。
表1 不同方法的日利用率
根據(jù)表1中的數(shù)據(jù)可知,所提方法對航班頻率優(yōu)化后,不同機(jī)型的日利用率均在10小時(shí)以上,文獻(xiàn)[4]方法對航班頻率優(yōu)化后,不同機(jī)型的日利用率與所提方法相比略差,文獻(xiàn)[3]方法對航班頻率優(yōu)化后,不同機(jī)型的日利用率最低。通過上述對比驗(yàn)證了所提方法的有效性。
在結(jié)構(gòu)上我國機(jī)場存在發(fā)展不平衡的問題,在航班數(shù)量不斷增加的背景下,我國樞紐機(jī)場已經(jīng)無法滿足要求。導(dǎo)致惡劣天氣下機(jī)場會出現(xiàn)大量航班延誤,機(jī)場擁堵的問題,航空運(yùn)輸具有網(wǎng)絡(luò)型的特點(diǎn),整個(gè)網(wǎng)絡(luò)都會受到航空延誤的影響,給游客和航空公司造成經(jīng)濟(jì)損失的同時(shí)降低了運(yùn)輸系統(tǒng)的效率。為了解決上述問題,亟需對多機(jī)場群航班頻率進(jìn)行優(yōu)化。目前多機(jī)場群航班頻率優(yōu)化模型存在延誤時(shí)間長、延誤航班數(shù)量多、出行成本高和日利用率低的問題。提出考慮出行成本的多機(jī)場群航班頻率優(yōu)化模型構(gòu)建方法,考慮旅客的出行成本,構(gòu)建了多機(jī)場群航班頻率優(yōu)化模型,并采用遺傳算法對多機(jī)場群航班頻率優(yōu)化模型求解,完成多機(jī)場群航班頻率優(yōu)化,解決了目前方法中存在的問題,提高了機(jī)場運(yùn)輸系統(tǒng)的效率。