張 明，張 明*，劉 凱，趙 敏

（1.南京航空航天大學(xué)民航學(xué)院，南京 211106；2.中國(guó)民用航空江蘇安全監(jiān)督管理局，南京 210022）

基于WITI的終端空域短期延誤預(yù)測(cè)

張 明1，張 明1*，劉 凱1，趙 敏2

（1.南京航空航天大學(xué)民航學(xué)院，南京 211106；2.中國(guó)民用航空江蘇安全監(jiān)督管理局，南京 210022）

針對(duì)終端空域內(nèi)短期天氣條件變化引起的航班延誤問(wèn)題，首先利用雷達(dá)回波平面位置顯示（PPI）和距離高度顯示（RHI）對(duì)降水等級(jí)分類，通過(guò)圖像識(shí)別方法對(duì)不同降水等級(jí)的云團(tuán)進(jìn)行輪廓提??；其次將對(duì)應(yīng)時(shí)間的交通流信息與降水云團(tuán)的空間范圍匹配，提出了不同降水等級(jí)條件下的三維天氣影響交通指數(shù)（WITI）計(jì)算方法；最后結(jié)合實(shí)際管制規(guī)則，將WITI指數(shù)與管制間隔調(diào)整和改航條件下的延誤時(shí)間作指數(shù)回歸分析，得到短期延誤預(yù)測(cè)模型。通過(guò)算例仿真驗(yàn)證了該模型的擬合優(yōu)度相比于線性回歸模型提高了1.6%～4.3%，且估計(jì)標(biāo)準(zhǔn)誤差趨近于0，說(shuō)明指數(shù)回歸延誤預(yù)測(cè)模型更能反映天氣影響下終端空域的實(shí)際延誤狀況。

終端空域；降水等級(jí)；交通流；WITI；延誤預(yù)測(cè)

近年來(lái)，天氣影響航班延誤的問(wèn)題一直受到國(guó)內(nèi)外學(xué)者的廣泛關(guān)注，對(duì)此展開(kāi)了大量的研究工作。Callaham等[1]提出天氣影響交通指數(shù)WITI概念，WITI指數(shù)可以反映天氣對(duì)交通流的影響程度。Chatterji[2]和Sridhar[3]擴(kuò)展了WITI的概念，并通過(guò)線性回歸方法對(duì)美國(guó)空域系統(tǒng)的延誤與WITI之間的關(guān)系建模。Klein等[4]結(jié)合航路WITI、終端區(qū)WITI和機(jī)場(chǎng)排隊(duì)延誤，提出NAS天氣影響指數(shù)NWX，考慮天氣、交通需求和機(jī)場(chǎng)排隊(duì)等因素分析延誤狀況。Hansen等[5]對(duì)航班延誤與交通流、機(jī)場(chǎng)天氣、航路對(duì)流天氣、天氣預(yù)報(bào)的準(zhǔn)確性等因素之間的關(guān)系建立延誤分析模型。Klein等[6]基于協(xié)同對(duì)流天氣預(yù)測(cè)產(chǎn)品CCFP建立了WITI預(yù)測(cè)模型，由于CCFP每2 h更新預(yù)報(bào)，更新頻率比較低，而麻省理工學(xué)院林肯實(shí)驗(yàn)室開(kāi)發(fā)的走廊綜合氣象系統(tǒng)（CIWS），每5 min的時(shí)間間隔可進(jìn)行一次天氣預(yù)報(bào)，更適用于短期延誤預(yù)測(cè)。Sridhar等[8]根據(jù)CIWS提供的天氣數(shù)據(jù)和ASPM提供的航班現(xiàn)在和過(guò)去的延誤數(shù)據(jù)，利用未來(lái)空中交通管理概念評(píng)估工具FACET計(jì)算出WITI和預(yù)測(cè)WITI（P-WITI），利用線性回歸模型預(yù)測(cè)美國(guó)空域系統(tǒng)未來(lái)2 h的延誤。Chen等[9]利用CIWS提供的降水云團(tuán)的位置和回波頂高等數(shù)據(jù)，結(jié)合航空器高度數(shù)據(jù)計(jì)算WITI指數(shù)，并利用自回歸分析得到延誤預(yù)測(cè)模型。Diana等[10]在終端空域天氣影響條件下，利用整數(shù)規(guī)劃模型動(dòng)態(tài)調(diào)整終端空域進(jìn)場(chǎng)路徑，緩解空域擁塞，并驗(yàn)證了在實(shí)際天氣影響下的適用性。Scot等[11]指出由于終端空域內(nèi)的航空器間隔相對(duì)較小、穩(wěn)定性不高，極易受到對(duì)流天氣的影響，提出基于對(duì)流天氣規(guī)避模型（CWAM）的對(duì)流天氣影響預(yù)測(cè)模型，并提供5 h的影響預(yù)測(cè)。Ison[12]通過(guò)調(diào)研發(fā)現(xiàn)影響終端區(qū)延誤的主要因素有天氣條件、交通流需求以及航空器的進(jìn)場(chǎng)路徑等，且終端區(qū)應(yīng)根據(jù)實(shí)際的運(yùn)行條件和管制規(guī)則建立獨(dú)有的延誤預(yù)測(cè)模型。以往對(duì)天氣影響下航班延誤預(yù)測(cè)的研究大多利用天氣影響交通指數(shù)WITI建立延誤預(yù)測(cè)模型，而WITI指數(shù)與降水等級(jí)的劃分閾值有關(guān)，并沒(méi)有全面考慮不同降水等級(jí)條件下，結(jié)合不同的管制規(guī)則分析延誤狀況。此外，國(guó)外關(guān)于天氣影響交通指數(shù)WITI與航班延誤的關(guān)系問(wèn)題大都針對(duì)國(guó)家空域網(wǎng)絡(luò)層面，而針對(duì)終端空域的復(fù)雜天氣引起延誤的研究較少。

本文主要討論終端空域在短期降水條件影響下的延誤時(shí)間預(yù)測(cè)。首先通過(guò)雷達(dá)回波平面位置顯示（PPI）和距離高度顯示（RHI）圖像的顏色劃分降水等級(jí)，并識(shí)別不同降水等級(jí)的圖像，提取輪廓以確定二維影響范圍和降水云團(tuán)的高度范圍；其次提出三維天氣影響交通指數(shù)WITI，運(yùn)用WITI的計(jì)算方法統(tǒng)計(jì)不同降水等級(jí)條件下受影響航空器數(shù)量和管制規(guī)則不同產(chǎn)生的延誤，利用指數(shù)回歸分析建立短期延誤預(yù)測(cè)模型；最后通過(guò)算例仿真驗(yàn)證在不同降水等級(jí)條件下的三維WITI指數(shù)計(jì)算方法以及延誤預(yù)測(cè)模型的正確性。

1 降水等級(jí)分類

為了解惡劣天氣對(duì)實(shí)際管制調(diào)配措施的影響，對(duì)管制員進(jìn)行調(diào)研訪問(wèn)。當(dāng)終端空域發(fā)生危險(xiǎn)天氣時(shí)，管制員將雷達(dá)回波圖像與終端空域結(jié)構(gòu)匹配，根據(jù)雷達(dá)回波圖像呈現(xiàn)的顏色，判斷惡劣天氣對(duì)進(jìn)場(chǎng)航段影響程度，調(diào)整管制間隔或飛行員申請(qǐng)改航，根據(jù)管制間隔以及改航路徑的不同，航班延誤時(shí)間也有所差異。

目前，國(guó)內(nèi)終端區(qū)的天氣預(yù)報(bào)產(chǎn)品主要是多普勒氣象雷達(dá)，利用多普勒效應(yīng)來(lái)測(cè)量云和降水粒子相對(duì)于雷達(dá)徑向運(yùn)動(dòng)速度，顯示方式有PPI和RHI。根據(jù)目標(biāo)對(duì)雷達(dá)波的反射率對(duì)應(yīng)降水云團(tuán)中水滴的直徑，水滴直徑越大，該降水云團(tuán)的雷達(dá)回波就越強(qiáng)，從而顯示出不同的色調(diào)[13]。根據(jù)雷達(dá)回波圖像中不同顏色的分布來(lái)判斷降水強(qiáng)度的大小。降水等級(jí)、顏色和回波強(qiáng)度的對(duì)應(yīng)關(guān)系[14]，如表1所示。

表1 降水等級(jí)、顏色和回波強(qiáng)度關(guān)系Tab.1 Precipitation grade，colors and echo intensity

當(dāng)天氣條件處于降水等級(jí)1時(shí)，即進(jìn)場(chǎng)航段受藍(lán)色降水天氣影響時(shí)，管制員會(huì)增大航空器的水平間隔以保證航空器安全飛行。根據(jù)雷達(dá)管制最低水平間隔標(biāo)準(zhǔn)，進(jìn)近管制不得小于6 km[15]，但同時(shí)考慮終端空域容量限制，實(shí)際管制間隔應(yīng)小于進(jìn)近管制間隔標(biāo)準(zhǔn)的2倍。設(shè)雷達(dá)管制最低水平間隔標(biāo)準(zhǔn)為d，管制間隔系數(shù)為α，受藍(lán)色降水天氣影響后的管制間隔為d′，則d′=αd。顯然管制間隔越大，航班的延誤時(shí)間越長(zhǎng)。

當(dāng)天氣條件處于降水等級(jí)2及以上時(shí)，即進(jìn)場(chǎng)航段受綠色、黃色、橙色、淡紅、紅色降水天氣影響，航空器將進(jìn)行改航、備降或返航。對(duì)于航空器改航情況，設(shè)受降水影響的進(jìn)場(chǎng)航段長(zhǎng)度為L(zhǎng)，改航距離系數(shù)為β，改航之后的路徑長(zhǎng)度為L(zhǎng)′，則L′=βL，β＞1。由于改航路徑是由飛行員根據(jù)機(jī)載雷達(dá)設(shè)備所提供的天氣信息進(jìn)行機(jī)動(dòng)飛行，所以改航路徑的距離對(duì)航班延誤時(shí)間有直接影響，改航距離越長(zhǎng)，航班的延誤時(shí)間越長(zhǎng)。

2 三維WITI計(jì)算模型

WITI指數(shù)是指受降水天氣影響的航空器數(shù)量，航空器數(shù)量越大，天氣影響程度越嚴(yán)重。本文將WITI的概念應(yīng)用于終端空域短期天氣條件變化的背景下，并結(jié)合多普勒氣象雷達(dá)提供的PPI和RHI圖像，以及管制部門(mén)提供的雷達(dá)軌跡數(shù)據(jù)，提出三維WITI指數(shù)的計(jì)算模型，其計(jì)算公式如下

其中：Pi和hi分別為第i架航空器在終端空域的平面坐標(biāo)和高度坐標(biāo)；W（t）為t時(shí)段的WITI指數(shù)；m（t）為t時(shí)段終端空域內(nèi)航空器數(shù)量；n（t）為t時(shí)段終端空域降水云團(tuán)數(shù)量；為航空器ai的平面位置判定參數(shù)；為航空器ai的高度位置判定參數(shù)；Wj為第j個(gè)降水云團(tuán)的二維平面輪廓內(nèi)的點(diǎn)集合；為第j個(gè)降水云團(tuán)的云底高；為第j個(gè)降水云團(tuán)的云頂高。

根據(jù)多普勒雷達(dá)提供的平面位置顯示PPI圖像，可利用Matlab圖像處理的顏色識(shí)別和輪廓提取方法獲得不同等級(jí)的降水云團(tuán)輪廓線，結(jié)合航空器的飛行軌跡，判斷航空器是否位于降水云團(tuán)所影響的范圍內(nèi)，若是；反之

航空器ai所在的高度層根據(jù)t時(shí)刻航空器所在的終端區(qū)進(jìn)場(chǎng)航段確定，根據(jù)機(jī)場(chǎng)標(biāo)準(zhǔn)儀表進(jìn)場(chǎng)程序，并結(jié)合實(shí)際管制員的行為習(xí)慣，可得到不同進(jìn)場(chǎng)航段所處的不同高度層。降水云團(tuán)的高度數(shù)據(jù)包括云底高、云頂高，根據(jù)多普勒氣象雷達(dá)提供的距離高度顯示RHI圖像獲取云高數(shù)據(jù)。若航空器的高度位于云底高與云頂高之間，則；反之

根據(jù)雷達(dá)回波圖像反映出的不同降水等級(jí)，將會(huì)有不同的管制間隔和改航路徑。通過(guò)Matlab識(shí)別出雷達(dá)回波圖像中藍(lán)色區(qū)域和綠色區(qū)域的邊界，并匹配交通流數(shù)據(jù)后，利用上述WITI指數(shù)計(jì)算方法可得到不同降水等級(jí)條件下受影響的航空器數(shù)量。

綜上，WITI指數(shù)是由終端空域內(nèi)的天氣條件和交通流狀況決定的，而終端空域內(nèi)的天氣條件和交通流狀況隨時(shí)間變化而變化，呈現(xiàn)動(dòng)態(tài)性，故WITI指數(shù)在不同時(shí)間段的數(shù)值將會(huì)不同。

3 終端空域短期延誤預(yù)測(cè)模型

根據(jù)降水等級(jí)分為兩種情景模式：

情景1 天氣條件為降水等級(jí)1，管制間隔系數(shù)為α。在情景1中引入管制間隔系數(shù)α，其意義為該天氣條件下2架航空器的管制間隔標(biāo)準(zhǔn)與正常天氣條件下管制間隔標(biāo)準(zhǔn)之間的比值，α取值在1～2之間。

情景2 天氣條件為降水等級(jí)2及以上，改航距離系數(shù)為β。在情景2條件下引入改航距離系數(shù)β，其意義在于航空器改航路徑的距離與原進(jìn)場(chǎng)航段距離的比值，β＞1。

文獻(xiàn)[10]對(duì)航路階段的延誤時(shí)間與天氣影響交通指數(shù)WITI建立了線性回歸模型，其研究結(jié)果表明延誤時(shí)間與WITI呈現(xiàn)明顯的線性關(guān)系；而針對(duì)終端空域的復(fù)雜天氣條件和交通流狀況，兩者有可能線性相關(guān)關(guān)系，更有可能是非線性關(guān)系。當(dāng)終端空域達(dá)到飽和容量時(shí)，其管制負(fù)荷也急劇增大，產(chǎn)生的航班延誤急劇增大，故對(duì)延誤時(shí)間與WITI建立指數(shù)回歸模型。

假設(shè)延誤時(shí)間與WITI指數(shù)之間的關(guān)系為

其中：D為延誤時(shí)間；W為WITI指數(shù)；a、b為相關(guān)系數(shù)。

4 算例仿真

本文以上海終端空域2016年04月02日12：00—20：30的降水過(guò)程為例，選取終端區(qū)氣象雷達(dá)回波圖像以及交通流數(shù)據(jù)，利用Visual Basic對(duì)上海終端區(qū)的交通流數(shù)據(jù)進(jìn)行仿真，采用Matlab的圖像處理工具對(duì)氣象雷達(dá)回波圖像進(jìn)行顏色識(shí)別與輪廓提取，得到天氣與交通流的匹配數(shù)據(jù)，計(jì)算三維WITI指數(shù)，并得到不同情景下的延誤預(yù)測(cè)模型，如圖1所示。

4.1 不同降水等級(jí)的WITI指數(shù)

第1步 不同降水等級(jí)下受影響的進(jìn)場(chǎng)航段。

情景1 針對(duì)圖1（a）中雷達(dá)回波PPI圖提取藍(lán)色區(qū)域輪廓，與上海終端空域結(jié)構(gòu)匹配，得到受影響的進(jìn)場(chǎng)航段為AND-DADAT-BELOP-XSY航段和AND-IDNIK-IGLIT-XSY航段，如圖1（c）所示。

情景2 針對(duì)圖1（a）中雷達(dá)回波PPI圖提取綠色區(qū)域輪廓，與上海終端空域結(jié)構(gòu)匹配，得到受影響的進(jìn)場(chǎng)航段為AND-IDNIK-IGLIT-XSY航段，如圖1（d）所示。

由于AND-IDNIK-IGLIT-XSY航段同時(shí)受到藍(lán)色區(qū)域和綠色區(qū)域覆蓋，故以AND-IDNIK-IGLITXSY進(jìn)場(chǎng)航段為例對(duì)不同天氣情景條件下的延誤狀況進(jìn)行分析。

為了精確得到影響AND-IDNIK-IGLIT-XSY航段降水云團(tuán)的高度范圍，對(duì)圖1（b）雷達(dá)回波RHI圖提取藍(lán)色輪廓和綠色輪廓得到圖1（e）和圖1（f）。由于RHI像素為600×600，而RHI圖像中坐標(biāo)軸原點(diǎn)的像素為（41，481），高度20.0 km的像素點(diǎn)為（41，118），故RHI圖像中高度范圍內(nèi)每個(gè)像素代表55 m。提取藍(lán)色輪廓圖中最高點(diǎn)的像素為（92，337），云頂高為（481-337）×55=7 920 m；提取綠色輪廓圖中最高點(diǎn)的像素為（65，395），其云頂高為（481-395）×55=4 730 m。

圖1 雷達(dá)回波PPI和RHI圖Fig.1 Radar echo PPI and RHI

由以上分析可知，AND-IDNIK-IGLIT-XSY進(jìn)場(chǎng)航段共同受藍(lán)色區(qū)域和綠色區(qū)域影響，且識(shí)別的藍(lán)色區(qū)域云頂高為7920m，綠色區(qū)域云頂高為4730m。而庵東進(jìn)場(chǎng)高度為5 100 m，DADAT報(bào)告點(diǎn)高度為3 500 m。

第2步 不同天氣情景條件下計(jì)算WITI指數(shù)。

將圖2中7架航空器自左向右分別編號(hào)i=1，…，7，2個(gè)降水云團(tuán)分別編號(hào)j=1，2。第1個(gè)降水云團(tuán)的云底高，云頂高；第2個(gè)降水云團(tuán)的云底高，云頂高。對(duì)執(zhí)行航班CCA198的第1架航空器a1，其平面位置位于綠色輪廓內(nèi)，故；航空器a1所在高度為4 620 m，在云高與之間，故。對(duì)其余6架航空器分別計(jì)算，結(jié)果如表2所示。

圖2 PPI綠色輪廓與交通流示意圖Fig.2 PPI green contour and air traffic flow

表2 航空器1～7的位置參數(shù)和高度參數(shù)Tab.2 Position and altitude parameters of Aircraft No.1～No.7

該時(shí)刻終端空域內(nèi)WITI指數(shù)W（t）計(jì)算如下

該空域內(nèi)有7架航空器，而WITI指數(shù)的計(jì)算結(jié)果為3，表示該空域內(nèi)只有a1、a4、a6這3架航空器受天氣影響；a2、a5為飛越航空器，飛行高度大于降水云團(tuán)云頂高，故不受天氣影響；航空器a3、a7在降水云團(tuán)輪廓外，故不受天氣影響。

由于12：00—20：30的降水過(guò)程中降水云團(tuán)的位置高度、降水強(qiáng)度均在實(shí)時(shí)動(dòng)態(tài)變化，AND-IDNIKIGLIT-XSY進(jìn)場(chǎng)航段受降水天氣影響的范圍和程度也在實(shí)時(shí)變化，故選取每0.5 h的雷達(dá)回波圖像與交通流數(shù)據(jù)進(jìn)行匹配，仿真統(tǒng)計(jì)WITI指數(shù)的周期為0.5 h，一共得到12 h的WITI變化趨勢(shì)，如圖3所示。其中B-WITI表示情景1天氣條件影響下的WITI指數(shù)，GWITI表示情景2天氣條件影響下的WITI指數(shù)。

圖3 AND-IDNIK-IGLIT-XSY航段WITI指數(shù)Fig.3WITI index of AND-IDNIK-IGLIT-XSY

4.2 不同情景的延誤預(yù)測(cè)模型

對(duì)情景1條件下設(shè)定不同的管制間隔系數(shù)α，情景2條件下設(shè)定不同的改航距離系數(shù)β，仿真得到的延誤時(shí)間分析如下：

情景1 在原有的前后機(jī)間隔為60 s的基礎(chǔ)上，分別設(shè)定管制間隔系數(shù)α=1.2、1.5、1.8。為了分析延誤時(shí)間與B-WITI指數(shù)的變化趨勢(shì)，將延誤時(shí)間縮小至原來(lái)的1/100，使兩者取得相同的數(shù)量級(jí)，如圖4所示。

圖4 α=1.2、1.5、1.8時(shí)延誤時(shí)間Fig.4 Delays when α=1.2，1.5，1.8

從圖4可看出，在12：00—20：30受降水條件影響的過(guò)程中，WITI指數(shù)和延誤時(shí)間的變化趨勢(shì)均相同，表明管制間隔影響下的延誤時(shí)間與WITI正相關(guān)，且管制間隔越大，延誤時(shí)間越大。

對(duì)延誤時(shí)間和WITI建立指數(shù)回歸模型，與文獻(xiàn)[10]的線性回歸模型相比，兩者的擬合曲線圖如圖5所示，模型相關(guān)系數(shù)、擬合優(yōu)度和估計(jì)標(biāo)準(zhǔn)誤差，如表3所示，對(duì)比分析得到指數(shù)模型的擬合優(yōu)度相對(duì)更優(yōu)，且指數(shù)模型的估計(jì)標(biāo)準(zhǔn)誤差SEE均趨于0。在管制間隔系數(shù)α=1.2、1.5、1.8的條件下，得到3個(gè)短期延誤預(yù)測(cè)模型，即

圖5 α=1.2、1.5、1.8時(shí)延誤時(shí)間與WITI指數(shù)的關(guān)系Fig.5 Relationship between delay and WITI when α=1.2，1.5，1.8

表3 α=1.2、1.5、1.8時(shí)模型匯總Tab.3 Model summary when α=1.2,1.5,1.8

情景2 原AND-IDNIK-IGLIT-XSY進(jìn)場(chǎng)航段距離為105.43 km，由于受降水天氣影響現(xiàn)規(guī)劃3條改航路徑，距離分別為126.62km、137.10km和159.25km，與原有進(jìn)場(chǎng)航段距離相比分別是其1.2、1.3、1.5倍。故將仿真程序中改航距離系數(shù)β分別設(shè)定為1.2、1.3、1.5。為了分析延誤時(shí)間與G-WITI指數(shù)的變化趨勢(shì)，將延誤時(shí)間縮小至原來(lái)的1/300，使兩者取得相同的數(shù)量級(jí)，如圖6所示。

從圖6可看出，在12：00—20：30受降水條件影響的過(guò)程中，WITI指數(shù)和延誤時(shí)間的變化趨勢(shì)均相同，表明改航影響下的延誤時(shí)間與WITI正相關(guān)，且改航距離越大，延誤時(shí)間越大。

對(duì)延誤時(shí)間和WITI建立指數(shù)回歸模型，與線性回歸模型相比，指數(shù)模型的擬合優(yōu)度相對(duì)更優(yōu)，如圖7所示，且指數(shù)模型的估計(jì)標(biāo)準(zhǔn)誤差SEE均趨于0，如表4所示。在改航距離系數(shù)β=1.2、1.3、1.5的條件下，得到3個(gè)短期延誤預(yù)測(cè)模型，即

綜上分析，對(duì)情景1和情景2條件下的WITI指數(shù)進(jìn)行比較，情景1條件下的WITI指數(shù)比情景2的大，表明雷達(dá)回波PPI圖像中藍(lán)色區(qū)域的影響范圍比綠色區(qū)域的大，故受降水等級(jí)1影響的航空器數(shù)量大；對(duì)情景1和情景2條件下的延誤時(shí)間進(jìn)行比較，改航情況下的延誤時(shí)間比增大管制間隔所引起的延誤時(shí)間更大，相應(yīng)的延誤成本也更大。

圖7 β=1.2、1.3、1.5時(shí)延誤時(shí)間與WITI的關(guān)系Fig.7 Relationship between delay and WITI when β=1.2,1.5,1.8

表4 β=1.2、1.3、1.5時(shí)模型匯總Tab.4 Model summary when β=1.2,1.3,1.5

5 結(jié)語(yǔ)

本文提出在不同降水等級(jí)條件下的三維終端空域WITI指數(shù)的計(jì)算方法，并利用線性回歸方法建立短期延誤預(yù)測(cè)模型，通過(guò)案例仿真驗(yàn)證了延誤預(yù)測(cè)模型的正確性和有效性。具體包括：①提出三維WITI指數(shù)計(jì)算模型，即考慮了降水云團(tuán)云頂高數(shù)據(jù)，與航空器高度進(jìn)行比較，確定航空器是否受到影響，更加精確統(tǒng)計(jì)了受降水天氣影響的航空器數(shù)量；②考慮在不同降水等級(jí)條件下不同管制措施對(duì)延誤時(shí)間的影響，即通過(guò)指數(shù)回歸分析得到短期延誤預(yù)測(cè)模型，相比線性回歸模型預(yù)測(cè)的延誤時(shí)間更加符合實(shí)際，分析得到延誤時(shí)間與天氣條件和管制員管制調(diào)配策略有關(guān)。

[1]CALLAHAM M B,DEARMON J S,COOPER A M,et al.Assessing NAS Performance：Normalizing for the Effects of Weather[C]//4th USA/Europe Air Traffic Management R&D Symposium,Santa Fe,NM,2001.

[2]CHATTERJIG,SRIDHARB.NationalAirspaceSystemDelayEstimation Using Weather Weighted Traffic Counts[C]//AIAA Guidance Navigation and Control Conference,San Francisco,2005:1-17.

[3]SRIDHAR B,SWEI S.Classification and Computation of Aggregate DelayUsingCenter-BasedWeatherImpactedTraffic Index[C]//7th AIAA AviationTechnology,IntegrationandOperations Conference(ATIO),Belfast,Northern Ireland,2007.

[4]KLEIN A,JEHLEN R,LIANG D.Weather Index With Queuing Component For National Airspace System Performance Assessment[C]//7th USA/Europe ATM R&D Seminar,Barcelona,Spain,2007.

[5]HANSEN M,XIONG J.Weather Normalization for Evaluating National Airspace System（NAS）Performance[C]//7th USA/Europe ATM R&D Seminar,Barcelona,Spain,2007.

[6]KLEIN A,KAVOUSSI S,HICKMAN D,et al.Predicting Weather Impact on Air Traffic[C]//ICNS Conference,Herndon,VA,2007.

[7]KLINGLE WILSON D,EVANS J.Description of the Corridor Integrated Weather System(CIWS)Weather Products,Project Report ATC-317[R].Lexington:MIT Lincoln Laboratory,2005.

[8]SRIDHAR B,CHEN N.Short term national airspace system delay prediction using weatherimpactedtrafficindex[J].Journalof Guidance Control and Dynamics,2009,32(2):657-662.

[9]CHEN N,SRIDHAR B.Estimation of Air Traffic Delay Using Three Dimensional Weather Information[C]//The 26th Congress of International CounciloftheAeronauticalSciences,Anchorage,Alaska,2008:8916.

[10]DIANA M,HAMSA B.Dynamic Reconfiguration of Terminal Airspace During Convective Weather[C]//49th IEEE Conference on Decision and Control,Atlanta,GA,2010:4875-4881.

[11]SCOT C,RICH D.Convective Weather Impact Forecasting in the Terminal Area[C]//14th AIAA Aviation Technology Integration and Operations Conference,Atlanta,GA,2014：2715.

[12]ISON D,WEILAND L,MCANDREW I,et al.Identification of air traffic management principles influential in the development of an airport arrival delay prediction model[J].Journal of Aviation/Aerospace Education&Research,2015,24(2):39-53.

[13]黃儀方.航空氣象[M].2版.成都：西南交通大學(xué)出版社,2011.

[14]LINDHOLM T,SHARMAN R,KROZEL J,et al.Translating Weather into Traffic Flow Management Impacts for NextGen[C]//14th Conference on Aviation Range and Aerospace Meteorology American Meteorological Society,Atlanta,GA,2010

[15]董襄寧.空中交通管理基礎(chǔ)[M].北京：科學(xué)出版社,2011.

（責(zé)任編輯：孟 欣）

Short-term terminal airspace delay prediction based on weather impacted traffic index

ZHANG Ming1,ZHANG Ming1*,LIU Kai1,ZHAO Min2
（1.College of Civil Aviation,Nanjing University of Aeronautics and Astronautics,Nanjing 211106,China; 2.Jiangsu Province Safety Oversight Bureau,CAAC,Nanjing 210022,China）

Aiming at the flight delay in terminal airspace caused by the changes of short-term weather conditions,radar echo PPI（plane-position indication）and RHI（range-height indication）are firstly used to classify the precipitation grades,then the contour of different grades of convective clouds is extracted by image recognition method.Secondly,the corresponding time of traffic flow information is matched to the convective clouds,and the calculation method of 3D WITI（weather impacted traffic index）based on different levels of precipitation is put forward.Finally,according to air traffic control rules,a short-term delay prediction model is proposed through exponential regression analysis between WITI and delay time under different control intervals or rerouting path distance.Simulation results show that under the condition of control interval adjustment and rerouting,compared with the linear regression model,the fit goodness of exponential regression model increases 1.6%～4.3%,and the standard error of estimate is closed to 0,indicating that the exponential regression delay prediction model can better reflect the actual delay situation of terminal airspace under the influence of precipitation.

terminal airspace;precipitation grades;air traffic flow;WITI;delay prediction

V355.1

A

1674－5590（2017）02－0005－06

2016-06-14；

2016-07-13

國(guó)家自然科學(xué)基金項(xiàng)目（U1233101，71271113）；中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)（NS2016062）

張明（1991—），男，山東臨沂人，碩士研究生，研究方向?yàn)榭罩薪煌ㄒ?guī)劃與管理.

張明*（1975—），男，江蘇南京人，副教授，博士，研究方向?yàn)榭罩薪煌ü芾?