基于成本指數(shù)的民機(jī)四維航跡預(yù)測優(yōu)化算法

樊劉仡 鄭智明 薛 飛 王 青

(上海飛機(jī)設(shè)計(jì)研究院，上海 201210)

0 引言

隨著社會(huì)經(jīng)濟(jì)的快速發(fā)展，空中交通運(yùn)輸需求與日俱增，空域擁堵、航班延誤等問題日益嚴(yán)重，給民航安全、高效運(yùn)行帶來了嚴(yán)峻的壓力。為了應(yīng)對目前空中交通運(yùn)輸上存在的問題，國際民航組織提出了全球空中交通管理運(yùn)行戰(zhàn)略[1-3]。該戰(zhàn)略指導(dǎo)未來空中交通管制從基于指令運(yùn)行(Command-Based Operation，簡稱CBO)向基于軌跡運(yùn)行(Trajectory-Based Operation，簡稱TBO)轉(zhuǎn)變，而基于軌跡運(yùn)行的理念在時(shí)間與空間上對航空器飛行軌跡的預(yù)測提出了更高精度的要求。

軌跡預(yù)測是按照飛行計(jì)劃中的控制方案驅(qū)動(dòng)飛機(jī)性能模型進(jìn)行飛行全過程的快速數(shù)值仿真[4]。它不僅能夠計(jì)算得到飛機(jī)飛行的三維軌跡，而且能夠預(yù)測出飛機(jī)在飛行軌跡上每一點(diǎn)的時(shí)間、空速、地速、質(zhì)量、油耗、升力和阻力等全狀態(tài)信息，得到飛行全過程的高維預(yù)測數(shù)據(jù)。為減少空氣污染，節(jié)省燃料消耗，除了在飛機(jī)設(shè)計(jì)階段采取多種優(yōu)化技術(shù)以外，也可以在運(yùn)營階段利用FMS計(jì)算最優(yōu)高度速度計(jì)劃表和最優(yōu)軌跡對飛機(jī)的性能進(jìn)行優(yōu)化[5-6]。

優(yōu)化參考軌跡主要有兩種手段，一是考慮飛行器運(yùn)動(dòng)方程(Equation of Motion，簡稱EOM)，用最優(yōu)控制的方法來得到最優(yōu)軌跡[7-8]；二是動(dòng)態(tài)查詢性能數(shù)據(jù)庫[9-10]，考慮在風(fēng)或者高度變化因素影響下的性能優(yōu)化問題。由于FMS計(jì)算效率不高，難以使用EOM進(jìn)行大量實(shí)時(shí)運(yùn)算，因此通常使用查詢性能數(shù)據(jù)庫(Performance Database，簡稱PDB)中的經(jīng)驗(yàn)表格的方法來完成計(jì)算。PDB中的數(shù)據(jù)以離散形式存儲，因此可以將軌跡優(yōu)化問題視為一個(gè)組合優(yōu)化問題。對于優(yōu)化指標(biāo)的選擇一般有燃油消耗和飛行時(shí)間兩種，飛行成本指數(shù)能夠?qū)烧哌M(jìn)行有效的權(quán)衡，因此本文提出了一種基于成本指數(shù)的方法對民用飛機(jī)的四維航跡預(yù)測進(jìn)行優(yōu)化。

1 成本指數(shù)與代價(jià)函數(shù)

飛行成本是飛行航跡優(yōu)化過程中必須考慮的關(guān)鍵指標(biāo)。一般來說，飛行成本可以分為固定成本和變動(dòng)成本。其中固定成本在一定范圍內(nèi)與航空公司業(yè)務(wù)量變化關(guān)系不大，譬如發(fā)動(dòng)機(jī)大修費(fèi)用、機(jī)組固定工資、月固定租賃費(fèi)用、飛行訓(xùn)練費(fèi)用等。變動(dòng)成本包括燃油成本和其他與飛行時(shí)間相關(guān)的成本(以下稱為時(shí)間成本)，譬如機(jī)組的小時(shí)薪水提成、小時(shí)租賃費(fèi)用、飛機(jī)維護(hù)成本等。

飛行航跡的優(yōu)化通常只考慮變動(dòng)成本，即燃油成本和時(shí)間成本，而不將固定成本考慮在內(nèi)。記單位燃油成本為Cf元/kg，單位時(shí)間成本為Ct元/min，則可以定義成本指數(shù)(Cost Index，簡稱CI)為：

(1)

單位為：kg/min。

成本指數(shù)一般由機(jī)組人員在飛機(jī)起飛前通過顯示控制系統(tǒng)的人機(jī)界面輸入到飛行管理計(jì)算機(jī)中。成本指數(shù)的取值反映了時(shí)間成本和燃油成本在飛行總成本中的構(gòu)成比例。當(dāng)成本指數(shù)為0時(shí)，表征時(shí)間成本可以忽略不計(jì)，飛行成本的全部權(quán)重集中于燃油成本，此時(shí)優(yōu)化的結(jié)果為最省油的飛行模式，即最大燃油里程模式。當(dāng)成本指數(shù)取最大值時(shí)，表征燃油成本可以忽略不計(jì)，飛行成本的全部權(quán)重集中于時(shí)間成本，此時(shí)優(yōu)化的結(jié)果為最短飛行時(shí)間。一般來說，按照最大燃油里程模式飛行速度太慢。當(dāng)飛機(jī)處于最大燃油里程狀態(tài)飛行時(shí)，以多消耗1%燃油為代價(jià)就可以換取3%～5%的速度提升。因此，在實(shí)際工程中，一般不會(huì)將成本指數(shù)設(shè)為0。同理，也不會(huì)將成本指數(shù)設(shè)為最大值。對于不同的飛機(jī)和不同的航空公司，成本指數(shù)的取值范圍也會(huì)有很大差別。以波音系列飛機(jī)為例，表1展示了部分型號飛機(jī)飛行過程中選用的成本指數(shù)的范圍。

表1 波音系列飛機(jī)成本指數(shù)取值范圍

在實(shí)際運(yùn)行過程中，通常會(huì)將成本指數(shù)設(shè)置為兩者之間較為合適的值。表2顯示了波音系列飛機(jī)成本指數(shù)推薦的取值。

表2 波音系列飛機(jī)成本指數(shù)推薦取值

通過成本指數(shù)，可以計(jì)算飛行的直接運(yùn)行成本(Direct Operating Cost, 簡稱DOC)：

(2)

其中,t0為起飛時(shí)刻，tf為降落時(shí)刻，f為燃油流量(單位時(shí)間內(nèi)消耗的燃油量)。

飛行成本的變化與巡航速度的選取有密切的關(guān)系，如圖1所示。時(shí)間成本隨巡航速度的增大而遞減，燃油成本隨巡航速度的變化會(huì)出現(xiàn)極小值，極小值處取得最大燃油里程速度。當(dāng)速度小于最大燃油里程速度時(shí)，燃油成本隨巡航速度增大而減小；當(dāng)速度大于最大燃油里程速度時(shí)，燃油成本隨巡航速度增大而增大。總成本隨巡航速度的變化也會(huì)出現(xiàn)極小值，極小值處取得最經(jīng)濟(jì)巡航速度。

成本指數(shù)對垂直剖面的影響如圖2所示。成本指數(shù)增加會(huì)導(dǎo)致爬升頂點(diǎn)和下降頂點(diǎn)在水平方向上的投影均向前運(yùn)動(dòng)。

圖2 成本指數(shù)變化對上升頂點(diǎn)和下降頂點(diǎn)位置的影響

針對軌跡優(yōu)化需要定下優(yōu)化的目標(biāo)，即獲取計(jì)算航程代價(jià)的函數(shù)，根據(jù)需求的不同，航程代價(jià)一般指航行中的燃油消耗、時(shí)間、成本指數(shù)CI的綜合代價(jià)。

定義航程成本如式(3)所示：

Cost=FBurned+CI×TFlight

(3)

其中，F(xiàn)Burned是總航程的燃油消耗，單位kg，TFlight是總航程時(shí)間，單位min，CI是成本指數(shù)，單位kg/min。CI一般由運(yùn)營方提供，CI越高，表明時(shí)間成本比重越高。

2 四維航跡優(yōu)化方法

根據(jù)ARINC702A中關(guān)于性能計(jì)算模塊的規(guī)定，機(jī)載性能數(shù)據(jù)庫主要由以下三個(gè)方面組成:

1)飛機(jī)氣動(dòng)模型數(shù)據(jù)庫，內(nèi)容包括升力系數(shù)、阻力極曲線等飛機(jī)氣動(dòng)外形數(shù)據(jù)；

2)發(fā)動(dòng)機(jī)數(shù)據(jù)庫，主要由發(fā)動(dòng)機(jī)推力模型和耗油特性數(shù)據(jù)組成，內(nèi)容包括發(fā)動(dòng)機(jī)推力模式、推力限制、推力使用狀態(tài)以及對應(yīng)的耗油特性等數(shù)據(jù)。

3)空中性能數(shù)據(jù)庫，內(nèi)容包括飛機(jī)在起飛、爬升、巡航、下降和著陸階段等飛行模態(tài)下的速度、燃油消耗、時(shí)間、水平距離和重量等數(shù)據(jù)。同時(shí)還包括基于性能數(shù)據(jù)二次計(jì)算得到的與航路相關(guān)的性能數(shù)據(jù)，以便機(jī)載性能管理。

飛機(jī)性能手冊中提供的數(shù)據(jù)包含曲線和表格兩種形式，將數(shù)據(jù)整理構(gòu)建成便于系統(tǒng)調(diào)用的數(shù)據(jù)庫，涵蓋飛行包線數(shù)據(jù)、發(fā)動(dòng)機(jī)數(shù)據(jù)和空中性能數(shù)據(jù)三個(gè)方面以供FMC調(diào)用。由于很多數(shù)據(jù)是離散存儲的，因此可以使用插值方法進(jìn)行制備和查詢。

為了尋找最優(yōu)的軌跡，需要計(jì)算按照不同速度高度計(jì)劃表飛行的航程代價(jià)并進(jìn)行比較，如果采取巡航階段階梯爬升，也需要計(jì)算爬升所帶來的額外代價(jià)。最后，利用窮舉的辦法，計(jì)算所有可行的軌跡，選擇總航程代價(jià)最低者作為最優(yōu)的垂直剖面。

根據(jù)第1節(jié)的描述，以航程代價(jià)最小為目標(biāo)時(shí)，軌跡優(yōu)化的主要環(huán)節(jié)在巡航階段，可以用基于成本指數(shù)的方法進(jìn)行優(yōu)化。

3 最經(jīng)濟(jì)巡航速度的計(jì)算

飛機(jī)在巡航階段飛行時(shí)，高度變化量近似為0，速度變化量近似為0，巡航過程中迎角很小，α也可近似視為0。在飛機(jī)水平航跡確定的情況下，其在垂直剖面內(nèi)的運(yùn)動(dòng)可以近似用質(zhì)心運(yùn)動(dòng)方程來描述：

(4)

其中，x為飛機(jī)在水平方向上的位置，v為飛行真空速，T為發(fā)動(dòng)機(jī)推力，D為阻力，L為飛機(jī)升力，m為飛機(jī)質(zhì)量，g為重力加速度，η為推力單位燃料消耗。

根據(jù)第2節(jié)對飛機(jī)性能模型的描述，可以得到如下關(guān)系：

在垂直剖面的優(yōu)化過程中，優(yōu)化目標(biāo)是飛行總成本：

其中，F(xiàn)為總?cè)加拖牧浚琭為燃油流量，t為飛行時(shí)間，CI為成本指數(shù)，Cf為燃油價(jià)格，Ct為單位時(shí)間成本，η為推力單位燃料消耗，D為飛機(jī)阻力。在優(yōu)化過程中給定了初始飛行時(shí)刻t0，初始飛機(jī)重量m0，初始位置x0。

考慮到燃油價(jià)格Cf是給定的值，因此優(yōu)化目標(biāo)可以寫成：

s.t.

x(t0)=x0,x(tf)=xf,m(t0)=m0,h(t)=h0

(7)

取v為優(yōu)化變量，使用HJB方程求解上述優(yōu)化問題，得到最經(jīng)濟(jì)巡航速度v*，并通過仿真的方法作出v*隨巡航高度的變化關(guān)系如圖3所示。圖中，各條曲線代表不同成本指數(shù)下的情況，巡航重量取240 t，成本指數(shù)變化范圍取0～100，每間隔10取一個(gè)樣本點(diǎn)；橫坐標(biāo)代表巡航高度，變化范圍取0 m～12 000 m，每隔100 m取一個(gè)樣本點(diǎn)；縱坐標(biāo)代表計(jì)算得到的最經(jīng)濟(jì)巡航速度。

圖3 最經(jīng)濟(jì)巡航速度隨高度的變化關(guān)系示意圖(不同成本指數(shù)下)

從圖中可以看出，當(dāng)巡航高度在7 000 m～11 000 m范圍內(nèi)變化時(shí)，計(jì)算出的最經(jīng)濟(jì)巡航速度與高度呈現(xiàn)近似線性關(guān)系。當(dāng)高度達(dá)到對流層頂后，溫度不再隨高度增加而變化，因此圖像上曲線斜率在該點(diǎn)處發(fā)生改變，如圖4所示。

圖4 最經(jīng)濟(jì)巡航速度隨高度的變化關(guān)系示意圖(7 000 m～12 000 m段)

由于v*與h之間存在局部近似線性關(guān)系，因此可以利用最小二乘的方法擬合曲線斜率。這樣，只要知道了給定成本指數(shù)下某一巡航高度對應(yīng)的最經(jīng)濟(jì)巡航速度，很快就能通過擬合公式計(jì)算其他高度對應(yīng)的最經(jīng)濟(jì)巡航速度。

4 最經(jīng)濟(jì)巡航高度的計(jì)算

類似的，以巡航高度為優(yōu)化變量，使用HJB方程求解上述優(yōu)化問題，發(fā)現(xiàn)在巡航速度確定后，巡航高度越高，巡航總成本越小。因此對于高度的優(yōu)化需要利用飛行包線的限制。

受飛機(jī)性能和空中交通管制約束，民航飛機(jī)只能在一定的飛行包線內(nèi)飛行。以波音747飛機(jī)為例，其巡航馬赫數(shù)限制不能超過0.89 Ma，巡航高度不能高于45 000 ft。取最大巡航馬赫數(shù)為0.85 Ma, 假設(shè)巡航高度變化范圍為7 000 m～12 000 m。飛機(jī)在某個(gè)巡航高度上飛行時(shí)，速度不能無限增大。加速過程中，當(dāng)最大推力等于阻力時(shí)，飛機(jī)達(dá)到最大速度。類似地，飛機(jī)飛行高度也不能無限增大。隨著飛行高度的增大，發(fā)動(dòng)機(jī)最大推力減小，直到不足以克服飛行受到的阻力，飛機(jī)的飛行高度下降。

根據(jù)飛機(jī)性能模型，有：

其中，ΔTeff=ΔT-CTc,4,0≤ΔTeff·CTc,5≤0.4,CTc,5≥0,CTc,3=0,CFcr=0.966。

又：

(9)

其中，

令最大巡航推力等于阻力，

(11)

則：

解得：

(13)

此即為最大巡航高度的約束條件，發(fā)現(xiàn)其與飛機(jī)巡航重量相關(guān)。飛機(jī)巡航重量減小時(shí)，最大巡航高度會(huì)隨之增大。對于起飛重量250 t的情況，可以計(jì)算最大巡航高度約為10 700 m，約合35 000 ft。

下面討論在固定飛行高度的情況下巡航速度的約束條件,可得：

在巡航重量240 t的條件下，改變巡航高度，可以得到速度約束如圖5所示。

圖5 最大巡航速度隨巡航高度變化曲線(推力約束下)

可見隨著巡航高度的增大，最大巡航速度限值隨之減小。取速度最大限值時(shí)飛行總成本的變化曲線如圖6所示，其中巡航高度變化范圍為10 000 m～10 600 m，每隔100 m取一個(gè)樣本點(diǎn)，成本指數(shù)取為30。發(fā)現(xiàn)在發(fā)動(dòng)機(jī)最大推力約束下，總成本隨巡航高度的增大呈現(xiàn)出先減小后增大的變化趨勢。

考慮0.85 Ma的最大巡航馬赫數(shù)限值。固定巡航速度為0.85 Ma，改變巡航高度，從7 000 m～10 000 m，每隔500 m取一個(gè)樣本點(diǎn)。得到飛行總成本隨巡航高度的變化曲線如圖7所示。在最大巡航馬赫數(shù)的條件下，隨著巡航高度的增大，飛行總成本逐漸減小。

圖6 取速度最大限值時(shí)飛行總成本隨巡航高度變化曲線(推力約束下)

圖7 飛行總成本隨巡航高度變化曲線(最大巡航馬赫數(shù)約束下)

綜合以上分析，可以得到最經(jīng)濟(jì)巡航高度的求解方法，如圖8所示。

圖8 最經(jīng)濟(jì)巡航高度求解示意圖

其中一族平行的直線表示不同成本指數(shù)下最經(jīng)濟(jì)巡航速度隨巡航高度的變化情況；黑色虛線表示最大巡航馬赫數(shù)的限制，本文中取0.85 Ma；黑色實(shí)線表示發(fā)動(dòng)機(jī)最大推力對巡航速度的限制。

以成本指數(shù)CI取30，起飛重量取250 t的情況為例進(jìn)行分析。圖8中①過程可以分解為④過程加⑤過程。在④過程中巡航速度不變，巡航時(shí)間成本近似不變，當(dāng)巡航高度增大時(shí)，由于大氣密度下降，阻力減小，因此燃油消耗減少，燃油成本降低，導(dǎo)致飛行總成本降低；在⑤過程中，由于是在高度不變的情況下向最經(jīng)濟(jì)巡航速度的方向變化，因此飛行總成本會(huì)降低。綜合對④、⑤過程的分析，可知①過程中飛行總成本是單調(diào)下降的，即隨著巡航高度的增加，最經(jīng)濟(jì)巡航速度對應(yīng)的飛行總成本單調(diào)遞減。②過程沿著最大巡航馬赫數(shù)限制線變化，根據(jù)圖7的實(shí)驗(yàn)結(jié)果知，當(dāng)巡航高度增加時(shí)，飛行總成本呈現(xiàn)下降趨勢。③過程沿著發(fā)動(dòng)機(jī)最大推力限制線變化，根據(jù)圖6的實(shí)驗(yàn)結(jié)果知，當(dāng)巡航高度增加時(shí)，飛行總成本呈現(xiàn)先減小后增大的變化趨勢。

根據(jù)對①、②、③過程的分析可知，在成本指數(shù)為30、起飛重量取250 t的情況下，最優(yōu)巡航條件應(yīng)該在③過程中取到。結(jié)合圖6的實(shí)驗(yàn)結(jié)果，最經(jīng)濟(jì)巡航高度應(yīng)取10 400 m，約合34 000 ft；最經(jīng)濟(jì)巡航速度應(yīng)取243 m/s，約合真空速472節(jié)，校正空速284節(jié)，巡航馬赫數(shù)0.817 Ma。

5 結(jié)論

本文綜合考慮民用飛機(jī)在飛行過程中燃油以及時(shí)間代價(jià)，基于成本指數(shù)提出了一種航跡預(yù)測優(yōu)化方法，以波音747飛機(jī)為例，結(jié)合飛機(jī)性能模型，理論推導(dǎo)了最經(jīng)濟(jì)巡航速度以及最經(jīng)濟(jì)巡航高度的計(jì)算方法，通過仿真驗(yàn)證了本文優(yōu)化方法的可行性，為民用飛機(jī)四維航跡精確預(yù)測的工程實(shí)現(xiàn)提供了理論參考。