場(chǎng)區(qū)飛機(jī)滑行效率評(píng)價(jià)模型的設(shè)計(jì)與仿真

耿 宏，朱 浩

(中國(guó)民航大學(xué)電子信息與自動(dòng)化學(xué)院，天津 300300)

1 引言

隨著經(jīng)濟(jì)的高速發(fā)展，航班量日益增長(zhǎng)，我國(guó)大型民用機(jī)場(chǎng)的航班延誤情況也開(kāi)始越來(lái)越多。雖然導(dǎo)致航班延誤的因素是多種多樣的，但是場(chǎng)區(qū)飛機(jī)滑行效率較低是影響航班延誤的主要原因[1]。分析影響場(chǎng)區(qū)飛機(jī)滑行效率的各種因素，準(zhǔn)確評(píng)價(jià)滑行效率，對(duì)減少航班延誤，改善飛機(jī)滑行階段的薄弱環(huán)節(jié)具有重要的價(jià)值和意義。

由于飛機(jī)在滑行過(guò)程中各種不確定性[2]，導(dǎo)致航班延誤問(wèn)題日益嚴(yán)重，相關(guān)專(zhuān)家也開(kāi)始將研究核心轉(zhuǎn)移到優(yōu)化飛機(jī)滑行路徑方面，確保場(chǎng)區(qū)飛機(jī)滑行效率得到有效提升。國(guó)內(nèi)外相關(guān)專(zhuān)家重點(diǎn)針對(duì)場(chǎng)區(qū)飛機(jī)滑行效率評(píng)價(jià)方面的內(nèi)容進(jìn)行了大量研究，并且給出了一些較好的研究成果。例如姜雨等人[3]將航天器在滑行系統(tǒng)中的最短滑行距離作為目標(biāo)函數(shù)，同時(shí)構(gòu)建滑行時(shí)空協(xié)同優(yōu)化模型，最終完成滑行效率評(píng)價(jià)。谷潤(rùn)平等人[4]分析飛機(jī)的導(dǎo)航能力，構(gòu)建飛行側(cè)向偏差模型和飛行高度偏差模型，計(jì)算實(shí)際飛行軌跡和定義航徑之間的偏差，進(jìn)而對(duì)偏差進(jìn)行統(tǒng)計(jì)，根據(jù)統(tǒng)計(jì)結(jié)果進(jìn)行評(píng)價(jià)。由于上述兩種評(píng)價(jià)模型未能在具體評(píng)價(jià)的過(guò)程中加入ISM-AHP-HFS聯(lián)合方法，導(dǎo)致評(píng)價(jià)結(jié)果的誤差和費(fèi)用的大幅度增加。為了更好地解決上述問(wèn)題，提出一種基于ISM-AHP-HFS的場(chǎng)區(qū)飛機(jī)滑行效率評(píng)價(jià)模型。測(cè)試結(jié)果表明，所提模型能夠有效降低評(píng)價(jià)費(fèi)用，獲取理想的評(píng)價(jià)結(jié)果。

2 場(chǎng)區(qū)飛機(jī)滑行效率評(píng)價(jià)模型

2.1 滑行效率影響因素識(shí)別和篩選

為了更好地構(gòu)建場(chǎng)區(qū)飛機(jī)滑行效率評(píng)價(jià)模型，需要優(yōu)先了解飛機(jī)進(jìn)場(chǎng)和離場(chǎng)的滑行過(guò)程。對(duì)于單個(gè)航班，即一架飛機(jī)來(lái)說(shuō)，滑行效率即飛機(jī)的滑行時(shí)間，從飛機(jī)停機(jī)位推出到離開(kāi)跑道起飛的時(shí)間總和。當(dāng)確定滑行路徑后，滑行時(shí)間=滑行距離/滑行速度。對(duì)于連續(xù)時(shí)間內(nèi)的航班而言，滑行效率主要是指設(shè)定時(shí)間內(nèi)所有航班的平均滑行時(shí)間。另外，滑行時(shí)間還包括飛機(jī)在停機(jī)位的等待和在跑道頭等待的總時(shí)間。部分歐美國(guó)家采用暢通滑行時(shí)間來(lái)衡量滑行效率[5，6]，對(duì)于整個(gè)機(jī)場(chǎng)來(lái)說(shuō)，滑行效率指的是航班整個(gè)過(guò)程的參與者，包括空管、航空公司、機(jī)場(chǎng)、旅客對(duì)滑行時(shí)間所做的“貢獻(xiàn)”，包括空管放行，地面服務(wù)保障，旅客安檢，滑行路徑規(guī)劃等。

在滑行的過(guò)程中，可能會(huì)出現(xiàn)以下幾方面的沖突：交叉沖突；對(duì)頭沖突和追尾沖突。為了更加全面掌握影響場(chǎng)區(qū)飛機(jī)滑行效率的因素，需要盡可能獲取影響飛機(jī)滑行的全部關(guān)鍵點(diǎn)。由于涉及的內(nèi)容比較廣泛，將全部影響滑行效率的因素進(jìn)行合并，獲取的客觀指標(biāo)和主觀指標(biāo)如表1和表2所示。

表1 滑行效率客觀指標(biāo)

通過(guò)影響指標(biāo)的相關(guān)定義，將上述指標(biāo)分別劃分為客觀指標(biāo)和主觀指標(biāo)，通過(guò)ISM分析各個(gè)指標(biāo)之間的關(guān)聯(lián)性[7，8]。ISM模型的主要目的就是將復(fù)雜的問(wèn)題簡(jiǎn)單化，盡量選取核心體系。設(shè)定aij代表布爾矩陣中的元素，則各個(gè)元素之間的關(guān)聯(lián)性能夠表示為

(1)

其中，可達(dá)矩陣是通過(guò)矩陣的形式表示在較長(zhǎng)一段時(shí)間內(nèi)一定長(zhǎng)度路徑后達(dá)到有向圖中各個(gè)節(jié)點(diǎn)的程度。優(yōu)先利用可達(dá)矩陣M計(jì)算可達(dá)集合P(Si)和先行集合Q(Sj)，則有

(2)

其中，可達(dá)集合中元素Si就是可達(dá)矩陣M上第i行對(duì)應(yīng)的影響因素取值；先行集合中元素Sj就是矩陣M中第j行對(duì)應(yīng)的影響因素取值。將上述兩種集合相結(jié)合，獲取共同集合L(Sij)：

L(Sij)={Si|P(Si)}∩Q(Sj)=P(Sij)

(3)

接下來(lái)，使用優(yōu)先的層級(jí)抽取規(guī)則方法進(jìn)行層級(jí)劃分，同時(shí)根據(jù)層級(jí)劃分結(jié)果將元素放置到對(duì)應(yīng)的層次中，同時(shí)將可達(dá)矩陣轉(zhuǎn)換為多個(gè)不同的矩陣，利用有向分支將各階要素進(jìn)行連接，構(gòu)建低階層次結(jié)構(gòu)圖[9]，進(jìn)而得到各個(gè)影響因素之間的主要關(guān)聯(lián)。

通過(guò)主成分分析法對(duì)各項(xiàng)指標(biāo)進(jìn)行識(shí)別和篩選，最終的指標(biāo)篩選結(jié)果如表3所示。

表3 評(píng)價(jià)指標(biāo)篩選結(jié)果

為了確保指標(biāo)的科學(xué)性，需要篩選出主要的場(chǎng)區(qū)飛機(jī)滑行效率指標(biāo)，對(duì)民航領(lǐng)域相關(guān)的專(zhuān)業(yè)人員進(jìn)行問(wèn)卷調(diào)查。分析最終的統(tǒng)計(jì)結(jié)果，刪除得到最低的四項(xiàng)，獲取有效的指標(biāo)，對(duì)指標(biāo)進(jìn)行分類(lèi)，得到如圖1所示的指標(biāo)體系。

圖1 滑行效率評(píng)價(jià)指標(biāo)體系

2.2 基于ISM-AHP-HFS的場(chǎng)區(qū)飛機(jī)滑行效率評(píng)價(jià)模型的建立

結(jié)合2.1小節(jié)對(duì)場(chǎng)區(qū)飛機(jī)滑行效率影響因素之間的關(guān)聯(lián)性分析，將分析結(jié)果設(shè)定為層次結(jié)構(gòu)模型的輸入，主要通過(guò)ISM-AHP-HFS三者相結(jié)合的方式構(gòu)建廠區(qū)滑行效率評(píng)價(jià)模型，具體的操作步驟如下所示：

通過(guò)對(duì)研究對(duì)象深入地調(diào)查和分析，獲取評(píng)價(jià)對(duì)象的關(guān)鍵要素[10]，同時(shí)將各個(gè)重要的影響因素通過(guò)層析結(jié)構(gòu)進(jìn)行分解，其中層次結(jié)構(gòu)模型如圖2所示。

圖2 層次結(jié)構(gòu)模型

對(duì)各個(gè)研究對(duì)象周邊的環(huán)境和自然因素進(jìn)行深入分析，進(jìn)而從代表性以及系統(tǒng)性等多個(gè)角度出發(fā)，組建場(chǎng)區(qū)飛機(jī)滑行效率評(píng)價(jià)指標(biāo)體系。為了更好觀察各個(gè)指標(biāo)的變化規(guī)律，優(yōu)先對(duì)全部評(píng)價(jià)指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理，如式(4)所示

(4)

其中，r(i，j)代表經(jīng)過(guò)標(biāo)準(zhǔn)化后的評(píng)價(jià)指標(biāo)取值；x(i，j)代表未進(jìn)行評(píng)價(jià)的樣本數(shù)據(jù)；xmax(i)代表樣本集中第i個(gè)評(píng)價(jià)指標(biāo)的最大值；xmin(i)代表樣本集中第i個(gè)評(píng)價(jià)指標(biāo)的最小值。

通過(guò)式(5)對(duì)越小越優(yōu)型指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理

(5)

通過(guò)式(6)對(duì)越中越優(yōu)型指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理

(6)

根據(jù)模糊評(píng)價(jià)矩陣R組建判斷矩陣B，主要目的是篩選具有價(jià)值的評(píng)價(jià)指標(biāo)。通過(guò)模糊綜合評(píng)價(jià)方法進(jìn)行評(píng)價(jià)的過(guò)程即選取最優(yōu)評(píng)價(jià)指標(biāo)的過(guò)程。接下來(lái)從綜合評(píng)價(jià)角度出發(fā)[11]，設(shè)定評(píng)價(jià)指標(biāo)1比評(píng)價(jià)指標(biāo)2中的樣本變化趨勢(shì)更加明顯，證明評(píng)價(jià)指標(biāo)1中傳遞的信息更多且更有價(jià)值。結(jié)合上述分析，對(duì)各個(gè)評(píng)價(jià)指標(biāo)的樣本差進(jìn)行計(jì)算，具體如式(7)所示

(7)

通過(guò)式(7)全面描述評(píng)價(jià)指標(biāo)對(duì)綜合評(píng)價(jià)的影響程度，同時(shí)應(yīng)用其構(gòu)建判斷矩陣B

(8)

式(8)代表各個(gè)評(píng)價(jià)指標(biāo)下樣本系列的平均取值。通過(guò)式(9)獲取1～9級(jí)的判斷矩陣，即

(9)

其中，smin和smax分別代表s(i)的最小值和最大值；bm代表相對(duì)重要性程度參數(shù)。

對(duì)判斷矩陣進(jìn)行一致性檢驗(yàn)以及修正等相關(guān)計(jì)算，要求判斷矩陣B必須滿(mǎn)足以下的約束條件

(10)

當(dāng)判斷矩陣的維數(shù)偏高時(shí)，通過(guò)相關(guān)定義計(jì)算判斷矩陣的特征值和特征向量。利用算法平均方法計(jì)算各個(gè)指標(biāo)的加權(quán)平均列向量，如式(11)所示

(11)

一致性比率計(jì)算需要優(yōu)先計(jì)算判斷矩陣的最大特征根，以下選取方根方法計(jì)算最大特征根。

優(yōu)先計(jì)算判斷矩陣M中各行的乘積幾何平均值di

(12)

將全部向量進(jìn)行歸一化處理，獲取如式(13)所示的特征向量

(13)

通過(guò)特征向量計(jì)算得到最大特征根，即

(14)

對(duì)判斷矩陣進(jìn)行一致性檢驗(yàn)，假設(shè)通過(guò)一致性檢驗(yàn)，則將全部向量進(jìn)行歸一化處理，得到對(duì)應(yīng)的權(quán)重向量；反之，則重新組建判斷矩陣。

在上述分析的基礎(chǔ)上，結(jié)合滑行效率評(píng)價(jià)指標(biāo)體系，邀請(qǐng)具有經(jīng)驗(yàn)的專(zhuān)家對(duì)各個(gè)大型機(jī)場(chǎng)中評(píng)價(jià)指標(biāo)對(duì)滑行效率的影響進(jìn)行評(píng)價(jià)打分，將重復(fù)數(shù)據(jù)刪除，同時(shí)利用專(zhuān)家的可信度，構(gòu)建具有可信度的猶豫模糊評(píng)價(jià)矩陣[12]。

通過(guò)具有可信度誘導(dǎo)猶豫模糊加權(quán)算子對(duì)猶豫模糊作為數(shù)據(jù)進(jìn)行集成，獲取信息集成值構(gòu)建信息集成矩陣，具體結(jié)果如表4所示。

表4 信息集成矩陣

通過(guò)猶豫模糊距離公式，分別計(jì)算各個(gè)指標(biāo)中的猶豫模糊基本元素和對(duì)應(yīng)屬性集成值之間的測(cè)度，獲取距離矩陣

D=(dij)10×4

(15)

根據(jù)標(biāo)準(zhǔn)，直接設(shè)定一致性閾值，將猶豫模糊源內(nèi)距離測(cè)度作為研究基礎(chǔ)，進(jìn)行一致性檢驗(yàn)。

對(duì)全部指標(biāo)進(jìn)行綜合一致性評(píng)分，根據(jù)評(píng)分結(jié)果組建場(chǎng)區(qū)飛機(jī)滑行效率評(píng)價(jià)模型，通過(guò)模型實(shí)現(xiàn)最終的評(píng)價(jià)。

3 仿真研究

為了驗(yàn)證所提基于ISM-AHP-HFS的場(chǎng)區(qū)飛機(jī)滑行效率評(píng)價(jià)模型的有效性，選取場(chǎng)面航空器滑行效率評(píng)價(jià)模型和民用飛機(jī)RNP導(dǎo)航能力的適航評(píng)估模型作為對(duì)比對(duì)象進(jìn)行測(cè)試。

為了驗(yàn)證評(píng)價(jià)模型獲取的評(píng)價(jià)結(jié)果是否可信，將均方根誤差、平均絕對(duì)誤差以及相對(duì)誤差設(shè)定為衡量指標(biāo)，具體的實(shí)驗(yàn)對(duì)比結(jié)果如表5～表7所示。

表5 不同模型的均方根誤差對(duì)比結(jié)果

表6 不同模型的平均絕對(duì)誤差對(duì)比結(jié)果

表7 不同模型的相對(duì)誤差對(duì)比結(jié)果

分析表5至表7中的實(shí)驗(yàn)數(shù)據(jù)可知，隨著實(shí)驗(yàn)次數(shù)的增加，各個(gè)方法的各種評(píng)價(jià)指標(biāo)也發(fā)生了十分明顯的變化。由于場(chǎng)面航空器滑行效率評(píng)價(jià)模型和民用飛機(jī)RNP導(dǎo)航能力的適航評(píng)估模型未能在評(píng)價(jià)的過(guò)程中加入ISM-AHP-HFS聯(lián)合方法，促使兩種模型的均方根誤差、平均絕對(duì)誤差以及相對(duì)誤差明顯更高一些，充分驗(yàn)證了所提模型的優(yōu)越性，同時(shí)也說(shuō)明所提模型能夠獲取更加精準(zhǔn)的評(píng)價(jià)結(jié)果，且評(píng)價(jià)結(jié)果更加可信。

為了更全面地驗(yàn)證所提模型的優(yōu)越性，以下實(shí)驗(yàn)測(cè)試對(duì)比三種模型的評(píng)價(jià)費(fèi)用，具體實(shí)驗(yàn)對(duì)比結(jié)果如圖3所示：

圖3 場(chǎng)區(qū)飛機(jī)滑行效率評(píng)價(jià)對(duì)比結(jié)果

分析圖3中的實(shí)驗(yàn)數(shù)據(jù)可知，所提模型的場(chǎng)區(qū)飛機(jī)滑行效率評(píng)價(jià)費(fèi)用在三種評(píng)價(jià)模型中為最低；場(chǎng)面航空器滑行效率評(píng)價(jià)模型的評(píng)價(jià)費(fèi)用次之；民用飛機(jī)RNP導(dǎo)航能力的適航評(píng)估模型的評(píng)價(jià)費(fèi)用最高。由此可見(jiàn)，所提模型加入ISM-AHP-HFS三種聯(lián)合方法進(jìn)行評(píng)價(jià)是有效且可行的。

4 結(jié)束語(yǔ)

現(xiàn)階段我國(guó)對(duì)于場(chǎng)區(qū)飛機(jī)滑行效率方面的研究?jī)?nèi)容仍然處于起步階段，開(kāi)展相關(guān)對(duì)于提升場(chǎng)區(qū)飛機(jī)滑行效率具有十分重要的意義。提出一種基于ISM-AHP-HFS的場(chǎng)區(qū)飛機(jī)滑行效率評(píng)價(jià)模型，仿真結(jié)果表明，所提模型能夠獲取準(zhǔn)確的評(píng)價(jià)結(jié)果，同時(shí)還能夠有效降低評(píng)價(jià)費(fèi)用。

所提模型由于受到研究時(shí)間和環(huán)境的限制，仍然存在一定的不足，后續(xù)將重點(diǎn)針對(duì)以下幾方面的內(nèi)容進(jìn)行研究：

1)進(jìn)一步研究外界因素對(duì)評(píng)價(jià)結(jié)果產(chǎn)生的影響。

2)綜合分析各個(gè)機(jī)場(chǎng)的實(shí)際情況，進(jìn)而對(duì)滑行路徑進(jìn)行規(guī)劃。