• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    考慮飛機排放因素的飛機巡航性能參數優(yōu)化方法

    2016-11-20 07:21:33魏志強張文秀韓博
    航空學報 2016年11期
    關鍵詞:燃油污染物飛機

    魏志強, 張文秀, 韓博

    1.中國民航大學 空中交通管理學院, 天津 300300

    2.中國民航大學 天津市空管運行規(guī)劃與安全技術重點實驗室, 天津 300300

    考慮飛機排放因素的飛機巡航性能參數優(yōu)化方法

    魏志強1,2*, 張文秀1, 韓博1

    1.中國民航大學 空中交通管理學院, 天津 300300

    2.中國民航大學 天津市空管運行規(guī)劃與安全技術重點實驗室, 天津 300300

    巡航是民航飛機主要的飛行階段,航班飛行中絕大多數的燃油消耗、飛行時間和污染物排放均發(fā)生在該階段?;趪H民航組織(ICAO)基準排放數據和BM2方法,建立了污染物排放量和排放成本計算模型;提出污染物排放價格權重、成本指數排放因子和綜合成本指數概念,以改進飛行成本計算模型,考慮污染物排放對飛行成本優(yōu)化的影響;通過計算分析輸入成本指數對性能參數和飛行成本的影響,建立了基于搜索方法的綜合成本指數優(yōu)化流程,并采用Visual Studio予以開發(fā)實現;然后計算了綜合成本指數飛行的經濟效益,以及飛機質量、成本指數和巡航高度對飛行成本的影響。結果表明,選擇合適的綜合成本指數和巡航高度,可使總飛行成本達到最小,提高運行經濟性。與傳統(tǒng)巡航相比,在典型情況下采用綜合成本指數巡航可以降低相同條件下的總飛行成本的0.15%。

    環(huán)境污染; 巡航性能; 綜合成本指數; 排放指數; 巡航高度

    根據民用飛機在飛行中廢氣排放的時空分布情況,飛機發(fā)動機對環(huán)境影響可分為兩大領域:全球范圍內的溫室效應和機場區(qū)域內的空氣質量問題[1]。為控制和減少航空業(yè)對溫室效應的影響,歐盟力爭將航空公司的飛機碳排放納入了碳排放交易體系[2]。根據發(fā)動機廠家審定數據,國際民航組織(International Civil Aviation Organization,ICAO)建立了機型排放數據庫,采用標準的起飛著陸循環(huán)(Landing and takeoff,LTO)來計算在機場范圍內1 000 m以下飛機的滑行、起飛、爬升、進近著陸等階段的CO、HC、SO2和NOx排放量,但對于耗時近80%的航路飛行排放,則沒有涉及[3]。

    為考慮航路飛行階段的排放影響,波音公司提出了BM2方法,實現對ICAO模型的修正。采用該模型,Tamara和Robert根據空中交通流量數據,對英國的航空排放問題進行了總體估算[4];基于澳大利亞2008年的空管航跡數據,Pham等研究建立了精度為1°×1°×1 000 ft (1 ft=30.48 cm)的四維排放數據庫清單[5];Owen等對2050年全球航空活動所造成的污染物排放數據進行了估算[6]。魏志強等依據建立了排放指數修正模型,實現對各個飛行階段的污染物排放量的估算分析[7-8]。上述研究實現了對污染物排放清單的估算,但都沒有涉及到如何通過優(yōu)化飛行參數來降低飛機污染物的排放量。

    巡航是民航飛機主要的飛行階段,航班飛行中絕大多數的燃油消耗、飛行時間和污染物物排放量均發(fā)生在該階段。傳統(tǒng)的研究主要是通過選擇合適的巡航方式、速度和高度,以達到節(jié)油或降低成本的目的。Singh引入成本指數的概念,從降低運行成本的角度研究了混合動力飛機垂直飛行剖面參數的優(yōu)化問題[9]。Miyazawa等構建了自由飛行概念下的飛機巡航參數計算模型,研究發(fā)動機特性參數和大氣因素對巡航性能及參數優(yōu)化的影響[10]。Erzberger和Lee研究了在推力受限條件下的巡航參數優(yōu)化問題。研究表明,對推力的約束會影響到最小運行成本條件下的最佳航程、巡航高度和速度[11]。Franco和Rivas在給定高度巡航參數優(yōu)化中,引入了航路隨機風的干擾,使得優(yōu)化出的巡航速度參數具有更高的魯棒性[12]。Valenzuela和Rivas研究了基于空中交通規(guī)則的約束的巡航參數優(yōu)化,以降低直接運營成本[13]。Sridhar和Grabbe研究了高空風對最佳巡航高度的影響,表明在最佳風速航跡飛行中,通過優(yōu)化垂直剖面可以降低1%~3%的飛機油耗與時間[14]。針對軍用大型運輸機巡航段航跡優(yōu)化問題,楊杰等提出了一種近似工程化的計算方法,實現對巡航參數的優(yōu)化計算[15];魏志強等研究了基于基本性能數據(Basic Aircraft Data, BADA)模型的遠程巡航方式(Long Range Cruise, LRC)下的巡航參數計算方法[16];王超等研究了基于BADA模型的飛行軌跡預測與優(yōu)化算法[17]。上述研究在對飛行參數進行優(yōu)化時,沒有考慮飛機排放的影響。

    通過優(yōu)化飛行參數以降低航班飛行對環(huán)境的影響,是國外研究的新動態(tài)之一。Schumann等建立了基于尾跡云和燃油消耗的航路優(yōu)化方法,以盡量減少航班運行對環(huán)境的影響[18]。Kaiser等建立了飛行航跡優(yōu)化模型,通過優(yōu)化飛行氣壓高度參數,以降低發(fā)動機排放對尾跡云形成和發(fā)展的影響[19]。Crewe將排放成本引入航班運行成本中,實現對飛行軌跡的優(yōu)化[20]。Cook等研究了不同污染物之間的權重關系,提出涵蓋排放成本的動態(tài)成本指數概念,并用于航班延誤分析[21]。

    通過對巡航參數的優(yōu)化,可以降低飛行活動對環(huán)境的影響?;贗CAO基準排放數據和BM2方法,本文建立了污染物排放成本計算模型;提出了污染物排放價格權重、成本指數排放因子和綜合成本指數概念,以改進飛行成本計算模型,考慮污染物排放對飛行成本優(yōu)化的影響;通過計算分析輸入成本指數對性能參數和飛行成本的影響,建立了基于搜索方法的綜合成本指數優(yōu)化流程,并采用Visual Studio予以編程開發(fā);最后分析了綜合成本指數飛行的經濟效益,以及飛機質量、成本指數和巡航高度對飛行成本的影響。

    1 污染物排放計算模型

    1.1 污染物排放成本計算模型

    飛機發(fā)動機的污染物排放種類主要包括CO2、NOx、HC、CO和SO2。其排放成本(CE)與污染物排放價格、排放量有關[19]。而排放量又取決于飛行中的燃油消耗量和發(fā)動機污染物排放指數(即單位油耗下的污染物排放量)。即

    (1)

    式中:下標i為污染物的類別;CEi為i類污染物的排放成本(¥),具體計算公式為

    CEi=PPiEPi=PPiFfEi

    (2)

    式中:PPi為i類污染物的排放價格(¥/kg);EPi為i污染物的排放量(kg);Ff為燃油流量(kg/h);Ei為i類污染物的排放指數(kg/kg)。為考慮污染物排放的影響,本文將污染物的排放價格(PPi,¥/kg)與燃油價格(PF,¥/kg)之比定義為污染物的排放價格權重(Wi),即

    (3)

    在此基礎上,提出無量綱的成本指數排放因子K,以反映污染物排放價格權重和排放指數的影響,具體定義為

    (4)

    將式(2)、式(3)和式(4)代入式(1),可以推導出給定巡航距離(R, km)下的污染物排放成本計算公式為

    (5)

    式中:VG為地速大小(m/s)。

    1.2 污染物排放量計算模型

    ICAO定義了標準的起飛降落循環(huán),包括進近著陸、滑行、起飛和爬升4個飛行階段。在此基礎上,發(fā)動機廠家基于試飛數據,提供了海平面、國際標準大氣模型(International Standard Atmosphere Model,ISAM)條件下的排放指數數據。本文以B737-800飛機(CFM56-7B24發(fā)動機)為例進行巡航參數的優(yōu)化研究,所用到的排放指數如表2所示[3]。

    由于飛機真實飛行條件不同于適航審定狀態(tài),需要根據實際飛行條件修正排放指數,以得到更為精確的污染物排放成本。在航空燃油種類一定的情況下,CO2和SO2的排放指數與其他因素無關,無需修正。由于飛機的實際巡航狀態(tài)與基準飛行階段所對應的飛行條件不同,需要依據ICAO的BM2方法對其他污染物的排放指數進行修正,以考慮巡航參數(燃油流量、飛行氣壓高度、溫度偏差等)的影響,具體如文獻[7]所示。

    表1 燃油價格與污染物排放價格信息[21]Table 1 Information about fuel price and pollution emission price[21]

    表2 CFM56-7B24發(fā)動機基準排放數據[3]Table 2 Basic emission data of CFM56-7B24 engine[3]

    2 民用飛機總飛行成本計算模型與實現

    2.1 傳統(tǒng)的飛行成本計算模型

    飛行成本是經濟巡航方式下的優(yōu)化目標。通常情況下的飛行成本包括燃油成本、時間成本和固定成本3部分。計算公式為

    C=CT+CF+CC=

    (6)

    式中:C為飛行成本(¥);CT為時間成本(¥);CF為燃油成本(¥);CC為固定成本(¥),與航空公司財務狀況有關,在優(yōu)化時可假定為固定常數;CI成本指數(kg/h),即時間價格與燃油價格之比。根據航空公司的運行狀況,財務部門計算出與上一年度的飛行時間成本(包括維修費、租賃費、飛行小時費等)和時間價格,然后用時間價格除以航油價格即可得到所需要的成本指數大小。

    在飛行前,飛行員可在飛行管理系統(tǒng)(FMS)中輸入CI,即可使FMS計算出優(yōu)化后巡航速度,并控制飛機按此最佳速度進行巡航,實現巡航成本的最優(yōu)化。對于B737-800而言,CI的取值范圍為0~999。CI越大優(yōu)化出的飛行速度就越大;當CI取到999時,FMS將忽略燃油價格,此時將會使飛行按最大限制速度巡航;反之當CI取0時對應于理論上最節(jié)油的最大航程速度。

    2.2 考慮排放影響的總飛行成本計算模型

    為考慮發(fā)動機污染物排放成本的影響,在式(6)中增加污染物排放成本的影響,進而得到總飛行成本的計算公式為

    C=CT+CF+CE+CC

    (7)

    再將式(5)代入式(7),可得

    (8)

    式中:CIH為綜合成本指數(kg/h),反映了污染物排放價格對經濟巡航的影響,其定義為

    (9)

    飛行成本是飛機巡航參數優(yōu)化的目標。文中式(6)為機載FMS在進行巡航性能優(yōu)化時的優(yōu)化目標計算公式,而式(8)為考慮污染物排放影響后的優(yōu)化目標計算公式。由于飛機上的FMS在進行巡航速度優(yōu)化時尚未考慮污染物的影響(即只能輸入成本指數),為此可以在航空公司財務部門提供的CI基礎上,根據式(9)計算出CIH后輸入到FMS系統(tǒng)中。從式(9)中可以看出,考慮污染物排放成本影響時的CIH要小于傳統(tǒng)的成本指數,即可以適度降低CI以使考慮污染物排放成本影響時總飛行成本達到最低。

    2.3 輸入成本指數對飛行性能參數及成本的影響

    從理論上講,CI是由公司財務部門確定的。但在實際飛行中,飛行員可以通過輸入不同的成本指數,滿足飛行性能優(yōu)化需要。在性能參數優(yōu)化中,由于燃油流量與飛機速度之間的關系很難顯性表示,致使Ei也無法直接根據飛行速度來計算。飛機廠家提供的性能軟件(如波音公司的BPS軟件)可以計算給定質量、巡航高度、溫度偏差和成本指數條件下的經濟巡航馬赫數、燃油流量,但沒有考慮污染物排放量的影響。

    本文采用BPS軟件計算不同輸入成本指數(CII)下的燃油流量數據,然后按照上述計算模型編程計算出在給定燃油價格、成本指數和污染物價格的情況下的總飛行成本(C),以考察輸入成本指數(CII)對飛行性能參數和飛行成本的影響。

    計算條件如下:燃油價格為6 000 ¥/t,成本指數為(CI)30,輸入成本指數(CII)計算范圍為0~60(即實際飛行中的常用取值范圍);燃油價格見表1,B737-800飛機,初始巡航質量68 t,飛行距離為100 km,巡航高度為7 800 m,溫度偏差為ISA+10,靜風。具體計算結果如表3和表4所示。

    從表3看出,在飛機發(fā)動機的污染物排放量中,CO2占絕大優(yōu)勢,是民航業(yè)節(jié)能減排的關注重點。但從表4的成本對比中可以看出,盡管NOx的排放量不及CO2的1%,但由于NOx的排放價格較高(見表1),使得其排放成本接近CO2的2/3,需要引起業(yè)界的高度重視。在其他參數一定情況下,改變輸入成本指數可以實現對飛機性能參數的調整,進而影響總飛行成本的大小及構成情況。飛機總飛行成本隨輸入成本指數增加會先減后增,因此可以優(yōu)化方法搜索到最小飛行成本所對應的輸入成本指數即為綜合成本指數。

    表3 污染物的排放量與輸入成本指數的關系Table 3 Relations between pollution emissions and inputted cost index

    表4 污染物排放成本與輸入成本指數的關系Table 4 Relations between pollution emission cost and inputted cost index

    2.4 綜合成本指數計算流程設計及軟件實現

    本文采取優(yōu)化遍歷搜索方法計算出給定大氣相對濕度(φ)、飛行氣壓高度(HP)、溫度偏差(ΔT)、飛機質量(m)、燃油價格(PF)、成本指數(CI)、航程(R)等條件下的綜合成本指數(CIH)。具體優(yōu)化計算流程圖如圖1所示。

    依據圖1所示的綜合成本指數計算流程圖和文中的相關計算公式,采用Visual Studio C#開發(fā)了民用飛機綜合巡航性能計算分析軟件,主要包括BPS軟件計算結果數據加載、大氣參數計算、排放指數修正計算、成本計算、最小成本條件判斷及綜合成本指數計算等功能模塊,能根據輸入的飛機質量、巡航高度、溫度偏差及成本指數等參數,實現對綜合成本指數的計算。

    圖1 綜合成本指數優(yōu)化計算流程圖Fig.1 Flowchart of optimizing calculation of integrated cost index

    3 基于綜合成本指數的巡航性能計算分析

    3.1 經濟效益分析

    傳統(tǒng)的巡航方式是按照給定的成本指數由機載飛行管理系統(tǒng)來對飛行參數進行優(yōu)化控制,沒有考慮污染物排放的成本影響。依據本文開發(fā)的民用飛機綜合巡航性能計算軟件,以B737-800飛機為例選取典型條件計算出綜合成本指數,并對巡航參數進行計算分析。具體計算條件為:巡航初始質量為65 t;巡航高度為8 900 m;巡航距離為1 000 km;標準大氣溫度環(huán)境;靜風;航空公司財務部門直接提供的成本指數(CI)為30(不考慮污染物排放影響);燃油價格6 000 ¥/t。表5的第2行為飛機按給定成本指數(CI=30)飛行的性能數據;第3行為按照優(yōu)化出的綜合成本指數(CIH=22)飛行時的性能數據。表6則為上述兩種飛行方式下的飛行成本對比。

    表5 成本指數巡航與綜合成本指數巡航性能數據的對比Table 5 Comparison of cruise performance data between cost index and integrated cost index

    表6 成本指數巡航與綜合成本指數巡航的對比Table 6 Comparison of cruise cost data between cost index and integrated cost index

    從表5、表6可以看出,計算出的綜合成本指數小于航空公司財務部門直接提供的成本指數大小,采用綜合成本指數飛行可降低相同條件下的總飛行成本的 0.15%。由于綜合成本指數偏小會導致飛行速度減緩,因此飛行時間略有增加,導致時間成本增加 2.6%。同時,燃油消耗會有所減少,燃油成本降低 1.1%。除CO外的其他污染物排放量也都會降低,其中絕對降幅和相對降幅最大的是NOx。因此,采用該飛行方式將會更有利于降低民航運輸飛行對大氣環(huán)境的影響。

    3.2 飛行成本的影響因素分析

    在飛機型號、污染物價格一定的情況下,飛機初始質量、成本指數和巡航高度是飛行成本的主要影響因素。圖2給出了飛機初始巡航質量對飛行成本的影響。計算條件:巡航高度9 500 m,巡航距離1 200 km,靜風,ISA+10,成本指數(CI)為35,燃油價格6 000 ¥/t。

    從圖2可以看出,隨著飛機質量增加,飛機的巡航時間基本不變,但燃油消耗有較大增加,使得燃油成本和污染物排放成本也有較大增加。圖3給出了不同巡航方式(成本指數)對飛行參數和飛行成本的影響。計算條件:初始巡航質量68 t,巡航距離1 500 km,靜風,ISA,巡航高度8 900 m,燃油價格6 000 ¥/t。

    從圖3可以看出,燃油成本和污染物成本基本與公司的成本指數無關,但由于時間成本隨成本指數增加而增加,使得總飛行成本也逐漸減增大。圖4給出了巡航高度對飛行參數和飛行成本的影響。計算條件:初始巡航質量70 t,巡航距離800 km,靜風,ISA,成本指數(CI)為30,燃油價格6 000 ¥/t。

    圖2 巡航初始質量對飛行成本的影響Fig.2 Influence of initial cruise mass on flight cost

    圖3 成本指數對飛行成本的影響Fig.3 Influence of cost index on flight cost

    圖4 巡航高度對飛行成本的影響Fig.4 Influence of cruise altitude on flight cost

    從圖4可以看出,在常用巡航高度范圍內,總飛行成本隨高度升高而先降后略增,其最低值對應于飛機的最佳巡航高度。

    通過對不同條件下總飛行成本的計算和對比(見圖2~圖4),在運行上可通過如下舉措來優(yōu)化飛機性能參數,降低航班運行成本。

    1) 在控制飛機質量方面,盡可能減少不必要的載重。

    2) 在成本指數方面,盡可能減少成本指數大小。為此需要航空公司進一步壓縮飛機維修費用和機組人員的小時費,通過降低時間價格來減小成本指數。

    3) 在巡航高度方面,可計算出不同飛行氣壓高度下的總飛行成本,然后盡可能地選擇最小成本對應的高度來巡航,降低總飛行成本。

    4 結 論

    基于ICAO基準排放數據和BM2方法,本文建立了污染物排放成本計算模型,改進了飛行成本計算模型,建立了基于搜索方法的綜合成本指數優(yōu)化流程,開發(fā)相應的計算軟件,以考慮污染物排放對飛行成本優(yōu)化的影響。在其他條件一定時,選擇合適的綜合成本指數和巡航高度,可使總飛行成本達到最小,提高運行經濟性。

    通過改變污染物的排放價格權重,還可以計算不同污染物種類關注度下的最佳巡航性能參數,實現飛行成本、燃油成本和排放成本的總體最優(yōu)化飛行。與傳統(tǒng)巡航相比,在典型情況下采用文中的綜合成本指數巡航可以降低相同條件下的總飛行成本的0.15%,其中飛行時間成本增加2.6%、燃油成本降低1.1%、污染物總體排放成本降低3.6%、各種污染物排放成本中降幅最大的是NOx。研究表明,采用該飛行方式將會更有利于降低民航運輸飛行對大氣環(huán)境的影響。研究結果可用于民航飛機的日常巡航性能優(yōu)化運行,提高航班運行的經濟性和環(huán)境友好性。

    [1] 張弛, 林宇震, 徐華勝, 等. 民用航空發(fā)動機低排放燃燒室技術發(fā)展現狀及水平[J]. 航空學報, 2014, 35(2): 332-350.

    ZHANG C, LIN Y Z, XU H S, et al. Development status and level of low emissions combustor technologies for civil aero-engine[J]. Acta Aeronautica et Astronautica Sinica, 2014, 35(2): 332-350 (in Chinese).

    [2] MALINA R, MCCONNACHIE D, WINCHESTER N, et al. The impact of the european union emissions trading scheme on US aviation[J]. Journal of Air Transport Management, 2012, 19(3): 36-41.

    [3] ICAO. Engine exhaust emissions data bank: Doc 9646-AN/943[R]. Montreal: ICAO, 2009.

    [4] TAMARA P, ROBERT B. Estimates of UK CO2emissions from aviation using air traffic data[J]. Climatic Change, 2007, 88(3-4): 367-384.

    [5] PHAM V, TANG J, ALAM S, et al. Aviation emission inventory development and analysis[J]. Environmental Modelling & Software, 2010, 25(12): 1738-1753.

    [6] OWEN B, LEE D S, LIM L. Flying into the future: Aviation emissions scenarios to 2050[J]. Environmental Science Technology, 2010, 44(10): 2255-2260.

    [7] 魏志強, 王超. 航班飛行各階段污染物排放量估算方法[J]. 交通運輸工程學報, 2010, 10(6): 48-52.

    WEI Z Q, WANG C. Estimating method of pollution emissions for scheduled flight in different phases[J]. Journal of Traffic and Transportation Engineering, 2010, 10(6): 48-52 (in Chinese).

    [8] 魏志強, 刁華智, 韓博. 民用飛機巡航階段污染物排放量計算研究[J]. 科學技術與工程, 2014, 14(19): 1617-1815.

    WEI Z Q, DIAO H Z, HAN B. Research on calculating of aircraft pollution emissions in cruise phase[J]. Science Technology and Engineering, 2014, 14(19): 1617-1815 (in Chinese).

    [9] SINGH R. Cost-based flight technique optimization for hybrid energy aircraft[J]. Aircraft Engineering & Aerospace Technology, 2014, 86(6): 591-598.

    [10] MIYAZAWA Y, HARADA A, WICKRAMASINGHE N K, et al. Effect of jet passenger aircraft performance model on the optimal periodic cruise maneuver[J]. Space Technology Japan the Japan Society for Aeronautical & Spaceences, 2013, 12: 99-105.

    [11] ERZBERGER H, LEE H. Constrained optimum trajectories with specified range[J]. Journal of Guidance Control & Dynamics, 2012, 3(1): 78-85.

    [12] FRANCO A, RIVAS D. Minimum-cost cruise at constant altitude of commercial aircraft including wind effects[J]. Journal of Guidance Control & Dynamics, 2011, 34(4): 1253-1260.

    [13] VALENZUELA A, RIVAS D. Optimization of aircraft cruise procedures using discrete trajectory patterns[J]. Journal of Aircraft, 2014, 51(51): 1632-1640.

    [14] SRIDHAR B, GRABBE S. Optimizing aircraft trajectories with multiple cruise altitudes in the presence of winds[J]. Journal of Aerospace Information Systems, 2014, 11(1): 35-47.

    [15] 楊杰, 薛建平, 王發(fā)威, 等. 大型運輸機巡航航跡優(yōu)化方案建模與分析[J]. 飛行力學, 2012, 30(4): 314-317.

    YANG J, XUE J P, WANG F W, et al. Cruise trajectory optimization modeling and analysis for large transport[J]. Flight Dynamics, 2012, 30(4): 314-317 (in Chinese).

    [16] WEI Z Q, WANG C, LIU F. Calculation of long range cruise performance based on the BADA model[C]//ICCTP 2010: Integrated Transportation Systems: Green, Intelligent, Reliable, 2010: 1097-1103.

    [17] 王超, 郭九霞, 沈志鵬. 基于基本飛行模型的4D航跡預測方法[J]. 西南交通大學學報, 2009, 44(2): 295-300.

    WANG C, GUO J X, SHEN Z P. Prediction of 4D trajectory based on basic flight models[J]. Journal of Southwest Jiaotong University, 2009, 44(2): 295-300 (in Chinese).

    [18] SCHUMANN U, GRAF K, MANNSTEIN H. Potential to reduce the climate impact of aviation by flight level changes[C] //3rd AIAA Atmosphere Space Environments Conference. Reston: AIAA, 2011.

    [19] KAISER M, ROSENOW J, FRICKE H, et al. Tradeoff between optimum altitude and contrail layer to ensure maximum ecological en-route performance using the enhanced trajectory prediction model (ETPM)[C]//International Conference on Application and Theory of Automation in Command and Control Systems, 2012: 127-134.

    [20] GREWE V, FR?MMING C, J?CKEL P, et al. Climate cost functions as a basis for climate optimized flight trajectories[C] //Tenth USA/Europe Air Traffic Management Research and Development Seminar (ATM2013), 2013.

    [21] COOK A, TANNER G, WILLIAMS V, et al. Dynamic cost indexing managing airline delay costs[J]. Journal of Air Transport Management, 2010: 26-35.

    魏志強男, 碩士, 副教授, 碩士生導師。主要研究方向: 飛機性能。

    Tel.: 022-24092905

    E-mail: weizhiqia@sina.com

    張文秀女, 碩士研究生。主要研究方向: 飛機性能。

    Tel.: 022-24092434

    E-mail: wxzhang@cauc.edu.cn

    韓博男, 博士, 講師。主要研究方向: 航空氣象與環(huán)境。

    Tel.: 022-24092905

    E-mail: bohan@cauc.edu.cn

    *Correspondingauthor.Tel.:022-24092905E-mail:weizhiqia@sina.com

    Optimizationmethodofaircraftcruiseperformanceparametersconsideringpollutionemissions

    WEIZhiqiang1,2,*,ZHANGWenxiu1,HANBo1

    1.CollegeofAirTrafficManagement,CivilAviationUniversityofChina,Tianjin300300,China2.TianjinKeyLaboratoryofOperationProgrammingandSafetyTechnologyofAirtrafficManagement,CivilAviationUniversityofChina,Tianjin300300,China

    Cruiseisthemainflightphaseforcivilaircraft,andmostofthefuelconsumption,flighttimeandpollutantsareconsumedorproducedinthisphase.ThepollutantemissioncalculationmodelandemissioncostmodelareestablishedbasedontheInternationalCivilAviationOrganization(ICAO)basicemissiondataandBM2methods.Inordertoconsidertheeffectofpollutionemission,theconceptsofpriceweightforpollutionemission,costindexemissionfactorandintegratedcostindexareproposedinthispapertoimprovethetotalflightcostmodel.Duetocalculatingandanalyzingtheinfluenceofinputtedcostindexonflightcost,theoptimumflowchartforcalculatingintegratedcostindexisestablishedbasedonoptimumsearchmethod.TheemissioncalculationsoftwareisdevelopedbyVisualStudioforcalculatingtheeconomicincomeforcruiseflightusingintegratedcostindex.Meanwhile,basedonresultdatacalculatedbythesoftware,theimpactsofinitialcruisemass,costindexandcruisealtitudeonflightcostarealsoanalyzed.Theresultsshowthattheflightcostcanachievetheminimumbychoosingadequateintegratedcostindexandcruisealtitude,asotherconditionsareconstant.Comparedwithtraditionalflightmode,thetotalflightcostcandecrease0.15%atgivenflightcondition.

    environmentalpollution;cruiseperformance;integratedcostindex;emissionindex;cruisealtitude

    2015-11-20;Revised2016-03-14;Accepted2016-04-11;Publishedonline2016-05-041435

    URL:www.cnki.net/kcms/detail/11.1929.V.20160504.1435.012.html

    s:NationalHigh-techResearchandDevelopmentProgramofChina(2014AA110501);NationalNaturalScienceFoundationofChina(U1533116,21407174);AeronauticalScienceFoundationofChina(20140267002);TianjinResearchProgramofApplicationFoundationandAdvancedTechnology(14JCQNJC08100)

    2015-11-20;退修日期2016-03-14;錄用日期2016-04-11; < class="emphasis_bold">網絡出版時間

    時間:2016-05-041435

    www.cnki.net/kcms/detail/11.1929.V.20160504.1435.012.html

    國家“863”計劃 (2014AA110501); 國家自然科學基金 (U1533116,21407174); 航空科學基金 (20140267002); 天津市應用基礎與前沿技術研究計劃 (14JCQNJC08100)

    *

    .Tel.:022-24092905E-mailweizhiqia@sina.com

    魏志強, 張文秀, 韓博. 考慮飛機排放因素的飛機巡航性能參數優(yōu)化方法J. 航空學報,2016,37(11):3485-3493.WEIZQ,ZHANGWX,HANB.OptimizationmethodofaircraftcruiseperformanceparametersconsideringpollutionemissionsJ.ActaAeronauticaetAstronauticaSinica,2016,37(11):3485-3493.

    http://hkxb.buaa.edu.cnhkxb@buaa.edu.cn

    10.7527/S1000-6893.2016.0119

    V212.13+3

    A

    1000-6893(2016)11-3485-09

    猜你喜歡
    燃油污染物飛機
    飛機失蹤
    菌株出馬讓畜禽污染物變廢為寶
    燃油泄漏闖了禍
    《新污染物治理》??鞲鍐⑹?/a>
    《新污染物治理》??鞲鍐⑹?/a>
    你能找出污染物嗎?
    “拼座飛機”迎風飛揚
    當代陜西(2019年11期)2019-06-24 03:40:28
    乘坐飛機
    神奇飛機變變變
    奔馳S500車燃油表不準
    一二三四在线观看免费中文在| 欧美黄色淫秽网站| 热99re8久久精品国产| 国产成人av教育| 性欧美人与动物交配| 国产成人啪精品午夜网站| 国产成人av激情在线播放| 看片在线看免费视频| 中文资源天堂在线| 国语自产精品视频在线第100页| 两性午夜刺激爽爽歪歪视频在线观看 | 中亚洲国语对白在线视频| 免费搜索国产男女视频| 国产极品粉嫩免费观看在线| 婷婷丁香在线五月| 久久午夜综合久久蜜桃| 国产成人精品久久二区二区免费| 久久中文字幕一级| 黄色女人牲交| 午夜免费激情av| 国内久久婷婷六月综合欲色啪| 欧美精品啪啪一区二区三区| 欧美在线黄色| 国产成人精品无人区| 久久久国产成人精品二区| 精品卡一卡二卡四卡免费| 久久久久久九九精品二区国产 | 1024视频免费在线观看| 色综合欧美亚洲国产小说| 韩国av一区二区三区四区| 999久久久国产精品视频| 亚洲精品美女久久久久99蜜臀| 91国产中文字幕| 精品久久久久久久久久久久久 | 精品久久久久久久毛片微露脸| 给我免费播放毛片高清在线观看| 国产麻豆成人av免费视频| 久久中文看片网| 一区二区日韩欧美中文字幕| 国产精品二区激情视频| 久久久久久国产a免费观看| 久久久久久大精品| 亚洲七黄色美女视频| 日韩av在线大香蕉| 禁无遮挡网站| 亚洲人成伊人成综合网2020| 中出人妻视频一区二区| 最新美女视频免费是黄的| 人人澡人人妻人| 精品高清国产在线一区| 国产伦在线观看视频一区| 老司机在亚洲福利影院| 一进一出抽搐gif免费好疼| 他把我摸到了高潮在线观看| 最近在线观看免费完整版| 91老司机精品| 男人舔奶头视频| 国产三级黄色录像| 成年免费大片在线观看| 曰老女人黄片| 国产精品久久久人人做人人爽| 午夜福利在线观看吧| 日本在线视频免费播放| 精品少妇一区二区三区视频日本电影| 国产私拍福利视频在线观看| 亚洲五月婷婷丁香| 50天的宝宝边吃奶边哭怎么回事| 91字幕亚洲| 亚洲av日韩精品久久久久久密| 美女 人体艺术 gogo| 一级毛片精品| 成人三级黄色视频| 国产精品久久久久久人妻精品电影| 老司机福利观看| 美女国产高潮福利片在线看| 欧美午夜高清在线| 操出白浆在线播放| 一本大道久久a久久精品| 亚洲国产毛片av蜜桃av| 亚洲精华国产精华精| 欧美国产日韩亚洲一区| 欧美亚洲日本最大视频资源| 淫秽高清视频在线观看| 91麻豆av在线| 国产亚洲精品久久久久久毛片| 亚洲av美国av| 精品无人区乱码1区二区| 黑人操中国人逼视频| 国产av在哪里看| 亚洲性夜色夜夜综合| 成人国语在线视频| 亚洲最大成人中文| 精品国产美女av久久久久小说| 国产成人影院久久av| 国产aⅴ精品一区二区三区波| 国产亚洲精品久久久久久毛片| 国产精品二区激情视频| 欧美日韩一级在线毛片| 亚洲专区字幕在线| 国产欧美日韩一区二区三| 桃红色精品国产亚洲av| 看免费av毛片| 97人妻精品一区二区三区麻豆 | 亚洲男人的天堂狠狠| 国产精品99久久99久久久不卡| 草草在线视频免费看| 国产精品久久视频播放| 免费观看精品视频网站| 久久久国产欧美日韩av| 久久久久九九精品影院| 窝窝影院91人妻| 婷婷精品国产亚洲av在线| 精品免费久久久久久久清纯| 天堂影院成人在线观看| 波多野结衣高清无吗| av在线天堂中文字幕| 可以免费在线观看a视频的电影网站| 欧美国产精品va在线观看不卡| 久久精品aⅴ一区二区三区四区| 这个男人来自地球电影免费观看| 麻豆一二三区av精品| 久久人人精品亚洲av| 亚洲七黄色美女视频| 欧美国产日韩亚洲一区| 男人舔女人的私密视频| 亚洲黑人精品在线| 免费在线观看日本一区| 天天躁狠狠躁夜夜躁狠狠躁| www.www免费av| 欧美一级a爱片免费观看看 | 亚洲中文av在线| 亚洲成a人片在线一区二区| www日本在线高清视频| 国产主播在线观看一区二区| av福利片在线| 色婷婷久久久亚洲欧美| 国产成年人精品一区二区| bbb黄色大片| 一本大道久久a久久精品| 男女做爰动态图高潮gif福利片| 大香蕉久久成人网| 搡老熟女国产l中国老女人| 精品国产一区二区三区四区第35| 免费看十八禁软件| 黄色毛片三级朝国网站| 欧美黑人精品巨大| bbb黄色大片| 脱女人内裤的视频| 国产爱豆传媒在线观看 | 夜夜夜夜夜久久久久| 亚洲一区中文字幕在线| 女生性感内裤真人,穿戴方法视频| 欧美色欧美亚洲另类二区| 色婷婷久久久亚洲欧美| 在线天堂中文资源库| e午夜精品久久久久久久| 女性被躁到高潮视频| 久久精品国产清高在天天线| 美女 人体艺术 gogo| 又黄又爽又免费观看的视频| 国产一区二区激情短视频| aaaaa片日本免费| 在线国产一区二区在线| 激情在线观看视频在线高清| 18禁观看日本| 中文字幕人妻丝袜一区二区| 免费搜索国产男女视频| av片东京热男人的天堂| 看免费av毛片| 在线观看免费视频日本深夜| 99re在线观看精品视频| videosex国产| 给我免费播放毛片高清在线观看| 成人三级做爰电影| 国产免费男女视频| 此物有八面人人有两片| 变态另类丝袜制服| 日韩大码丰满熟妇| 国产精品亚洲av一区麻豆| 亚洲中文日韩欧美视频| 激情在线观看视频在线高清| 亚洲国产精品sss在线观看| or卡值多少钱| 在线观看免费日韩欧美大片| 亚洲精品国产一区二区精华液| 亚洲一区二区三区不卡视频| 夜夜看夜夜爽夜夜摸| 亚洲片人在线观看| 国产精品av久久久久免费| 国产日本99.免费观看| 久久精品国产综合久久久| 久久久久久九九精品二区国产 | 美女午夜性视频免费| 女性生殖器流出的白浆| a在线观看视频网站| 最新在线观看一区二区三区| 日日摸夜夜添夜夜添小说| 亚洲成人免费电影在线观看| 日韩精品免费视频一区二区三区| 亚洲狠狠婷婷综合久久图片| 中文字幕人妻熟女乱码| 国产精品一区二区免费欧美| 日本成人三级电影网站| 丁香欧美五月| 亚洲人成77777在线视频| bbb黄色大片| 亚洲精品av麻豆狂野| 天天躁夜夜躁狠狠躁躁| 女性生殖器流出的白浆| 很黄的视频免费| 国产精品电影一区二区三区| 亚洲人成网站在线播放欧美日韩| 国产真人三级小视频在线观看| 丝袜美腿诱惑在线| 国产真实乱freesex| 中文字幕久久专区| 日韩中文字幕欧美一区二区| 国产精品,欧美在线| 亚洲黑人精品在线| 老司机福利观看| 在线观看舔阴道视频| 波多野结衣高清无吗| 欧美激情极品国产一区二区三区| 一本一本综合久久| 国产成人av激情在线播放| 麻豆成人午夜福利视频| 午夜久久久久精精品| 99国产精品一区二区三区| 亚洲精品在线观看二区| or卡值多少钱| 欧美黄色淫秽网站| 一级黄色大片毛片| 波多野结衣高清作品| 黄色片一级片一级黄色片| 美国免费a级毛片| 亚洲国产看品久久| 国产高清有码在线观看视频 | 岛国在线观看网站| 黄片播放在线免费| 久久婷婷成人综合色麻豆| 99在线人妻在线中文字幕| 精品一区二区三区四区五区乱码| 99国产精品一区二区三区| av欧美777| 亚洲精华国产精华精| 国产精品一区二区免费欧美| 欧美又色又爽又黄视频| 欧美性猛交╳xxx乱大交人| 给我免费播放毛片高清在线观看| 国产激情偷乱视频一区二区| 免费在线观看影片大全网站| 国产91精品成人一区二区三区| 男女之事视频高清在线观看| 老司机午夜十八禁免费视频| 中文字幕av电影在线播放| 老汉色∧v一级毛片| 国产一卡二卡三卡精品| 久久性视频一级片| 久久国产精品人妻蜜桃| 久久久久久免费高清国产稀缺| 他把我摸到了高潮在线观看| 嫁个100分男人电影在线观看| 日韩欧美国产在线观看| 国产精品1区2区在线观看.| 久热这里只有精品99| 日韩高清综合在线| 黄色a级毛片大全视频| 在线观看66精品国产| 99精品欧美一区二区三区四区| 女性生殖器流出的白浆| 日韩 欧美 亚洲 中文字幕| 亚洲七黄色美女视频| 中文在线观看免费www的网站 | 亚洲精品一卡2卡三卡4卡5卡| 国产精品一区二区三区四区久久 | 99精品在免费线老司机午夜| 精品国产超薄肉色丝袜足j| 日韩精品免费视频一区二区三区| 99国产精品一区二区蜜桃av| 欧美丝袜亚洲另类 | 久久狼人影院| 少妇裸体淫交视频免费看高清 | 亚洲精华国产精华精| 久久午夜综合久久蜜桃| 亚洲 欧美 日韩 在线 免费| 啦啦啦 在线观看视频| 亚洲成a人片在线一区二区| 亚洲精品中文字幕在线视频| 日韩av在线大香蕉| 亚洲一区二区三区不卡视频| 久久久久久国产a免费观看| 亚洲天堂国产精品一区在线| 午夜福利一区二区在线看| 女人被狂操c到高潮| 国产爱豆传媒在线观看 | 欧美大码av| 在线观看一区二区三区| 亚洲欧美激情综合另类| 神马国产精品三级电影在线观看 | 国产成人av激情在线播放| 久久精品国产亚洲av香蕉五月| 桃色一区二区三区在线观看| 亚洲精品在线观看二区| 美女 人体艺术 gogo| 亚洲精品中文字幕在线视频| 欧美乱妇无乱码| 亚洲精品一卡2卡三卡4卡5卡| 男女之事视频高清在线观看| 日韩国内少妇激情av| 嫩草影院精品99| 窝窝影院91人妻| 久久久久久久久久黄片| 中文字幕最新亚洲高清| 一区二区日韩欧美中文字幕| svipshipincom国产片| 黄色成人免费大全| 成人午夜高清在线视频 | 18美女黄网站色大片免费观看| 又黄又粗又硬又大视频| 91av网站免费观看| 国产成+人综合+亚洲专区| 久久久国产精品麻豆| 50天的宝宝边吃奶边哭怎么回事| 免费电影在线观看免费观看| 观看免费一级毛片| 久久久久久久久久黄片| 精品久久久久久久久久免费视频| 88av欧美| 日韩大尺度精品在线看网址| 99在线人妻在线中文字幕| 欧美国产日韩亚洲一区| 国产午夜福利久久久久久| www.999成人在线观看| 中文字幕人成人乱码亚洲影| 嫁个100分男人电影在线观看| 精品卡一卡二卡四卡免费| 久久精品aⅴ一区二区三区四区| 国产av又大| 在线十欧美十亚洲十日本专区| 午夜福利一区二区在线看| bbb黄色大片| 欧美日韩瑟瑟在线播放| 国产精品,欧美在线| 午夜精品久久久久久毛片777| 午夜日韩欧美国产| 亚洲精品一区av在线观看| 母亲3免费完整高清在线观看| 怎么达到女性高潮| 色在线成人网| 午夜福利在线观看吧| 可以在线观看毛片的网站| 成人国产一区最新在线观看| 大型av网站在线播放| 黄色a级毛片大全视频| 国产三级在线视频| 色精品久久人妻99蜜桃| 一级毛片女人18水好多| 村上凉子中文字幕在线| 大香蕉久久成人网| 少妇熟女aⅴ在线视频| 国产精品av久久久久免费| 中文资源天堂在线| 国产精品av久久久久免费| 亚洲国产日韩欧美精品在线观看 | 我的亚洲天堂| 免费高清视频大片| 久久国产亚洲av麻豆专区| 99在线人妻在线中文字幕| 看免费av毛片| 亚洲无线在线观看| 国产伦一二天堂av在线观看| 国产亚洲精品久久久久5区| 中文字幕精品免费在线观看视频| 亚洲七黄色美女视频| 天天一区二区日本电影三级| 国产黄色小视频在线观看| 超碰成人久久| 亚洲一码二码三码区别大吗| 欧美精品亚洲一区二区| 一本久久中文字幕| 老司机靠b影院| 午夜精品久久久久久毛片777| 日韩欧美 国产精品| 1024手机看黄色片| 嫩草影视91久久| 午夜老司机福利片| 亚洲美女黄片视频| 熟女少妇亚洲综合色aaa.| 亚洲片人在线观看| 久久人妻av系列| 深夜精品福利| 亚洲av中文字字幕乱码综合 | 欧美绝顶高潮抽搐喷水| 两个人视频免费观看高清| 99国产精品一区二区三区| 成人午夜高清在线视频 | 国产熟女午夜一区二区三区| 日韩精品青青久久久久久| 热99re8久久精品国产| 午夜免费成人在线视频| 久久精品国产99精品国产亚洲性色| 亚洲国产精品sss在线观看| 欧美成人免费av一区二区三区| 国产又爽黄色视频| 久久久久精品国产欧美久久久| 免费av毛片视频| 日本一本二区三区精品| 久久久久久国产a免费观看| 午夜日韩欧美国产| 狠狠狠狠99中文字幕| 老司机在亚洲福利影院| 亚洲五月婷婷丁香| 亚洲 国产 在线| aaaaa片日本免费| 精品久久久久久,| 成熟少妇高潮喷水视频| 久久天堂一区二区三区四区| 精品一区二区三区av网在线观看| av有码第一页| 国产亚洲精品一区二区www| 亚洲av第一区精品v没综合| 成年版毛片免费区| 国产色视频综合| 国产爱豆传媒在线观看 | 韩国av一区二区三区四区| 制服丝袜大香蕉在线| 欧美日韩瑟瑟在线播放| 99riav亚洲国产免费| 国产99白浆流出| 国产熟女午夜一区二区三区| 亚洲精品一区av在线观看| 久久久久国产一级毛片高清牌| 亚洲片人在线观看| 淫妇啪啪啪对白视频| 亚洲精品粉嫩美女一区| 女人爽到高潮嗷嗷叫在线视频| 黄频高清免费视频| 亚洲精品一卡2卡三卡4卡5卡| 在线观看免费日韩欧美大片| 久久精品人妻少妇| 可以免费在线观看a视频的电影网站| 亚洲片人在线观看| 精品久久久久久久久久免费视频| 露出奶头的视频| 欧美性长视频在线观看| 丝袜人妻中文字幕| 老熟妇仑乱视频hdxx| 97超级碰碰碰精品色视频在线观看| 国产又爽黄色视频| 久久人妻av系列| 国产三级黄色录像| 亚洲va日本ⅴa欧美va伊人久久| 成人三级做爰电影| 亚洲片人在线观看| 午夜免费鲁丝| 久久久精品欧美日韩精品| 天天一区二区日本电影三级| 国产成人影院久久av| 一区二区三区高清视频在线| 亚洲成国产人片在线观看| 一区福利在线观看| 日韩大尺度精品在线看网址| 丝袜美腿诱惑在线| 一级黄色大片毛片| 女警被强在线播放| 美女大奶头视频| 精品久久久久久久久久久久久 | 亚洲国产精品久久男人天堂| 一二三四社区在线视频社区8| 国产99久久九九免费精品| 亚洲久久久国产精品| 中文字幕最新亚洲高清| 男女做爰动态图高潮gif福利片| 成熟少妇高潮喷水视频| 2021天堂中文幕一二区在线观 | 亚洲一区二区三区色噜噜| 老鸭窝网址在线观看| 免费在线观看视频国产中文字幕亚洲| 一级片免费观看大全| 999精品在线视频| 亚洲成a人片在线一区二区| 午夜a级毛片| 男女那种视频在线观看| 琪琪午夜伦伦电影理论片6080| 欧美性长视频在线观看| 久久性视频一级片| 亚洲真实伦在线观看| 欧洲精品卡2卡3卡4卡5卡区| 日韩三级视频一区二区三区| 香蕉久久夜色| 在线观看午夜福利视频| 亚洲av中文字字幕乱码综合 | 白带黄色成豆腐渣| 最近在线观看免费完整版| 亚洲真实伦在线观看| 欧洲精品卡2卡3卡4卡5卡区| 日韩三级视频一区二区三区| 国产精品乱码一区二三区的特点| 国产成人影院久久av| 久久精品人妻少妇| 国产成人啪精品午夜网站| АⅤ资源中文在线天堂| 日本熟妇午夜| 香蕉国产在线看| 午夜影院日韩av| 天堂影院成人在线观看| 免费高清在线观看日韩| 国产精品免费视频内射| 午夜老司机福利片| 国产高清激情床上av| 看免费av毛片| 最近最新中文字幕大全免费视频| 亚洲一区二区三区色噜噜| 成人国语在线视频| 岛国视频午夜一区免费看| 亚洲av成人一区二区三| 亚洲精华国产精华精| 亚洲专区中文字幕在线| 美女扒开内裤让男人捅视频| 国产一区二区在线av高清观看| 亚洲自偷自拍图片 自拍| 色播亚洲综合网| 国产高清视频在线播放一区| 国产亚洲欧美精品永久| 波多野结衣高清作品| 国产免费男女视频| 一本一本综合久久| 啦啦啦免费观看视频1| 黄色a级毛片大全视频| 国产亚洲精品久久久久5区| 亚洲国产高清在线一区二区三 | 成人av一区二区三区在线看| 91麻豆av在线| 日韩三级视频一区二区三区| 美国免费a级毛片| 最近最新中文字幕大全免费视频| 老熟妇乱子伦视频在线观看| 熟女电影av网| 男女视频在线观看网站免费 | 人人妻人人看人人澡| 19禁男女啪啪无遮挡网站| 一级毛片女人18水好多| 亚洲熟妇中文字幕五十中出| 国产在线精品亚洲第一网站| 中文字幕人妻丝袜一区二区| 黄片播放在线免费| 亚洲国产欧洲综合997久久, | 亚洲精品色激情综合| 淫妇啪啪啪对白视频| 一二三四社区在线视频社区8| 可以在线观看毛片的网站| 听说在线观看完整版免费高清| 国产熟女xx| av视频在线观看入口| 少妇 在线观看| 国产亚洲av嫩草精品影院| 悠悠久久av| 国产成人av教育| 久久久国产欧美日韩av| 老司机福利观看| 欧美午夜高清在线| 国产在线精品亚洲第一网站| 19禁男女啪啪无遮挡网站| 午夜日韩欧美国产| 中文字幕另类日韩欧美亚洲嫩草| 日韩精品中文字幕看吧| 日本 欧美在线| 不卡一级毛片| 亚洲av电影不卡..在线观看| 女生性感内裤真人,穿戴方法视频| 777久久人妻少妇嫩草av网站| 美女 人体艺术 gogo| 香蕉av资源在线| 无限看片的www在线观看| 1024手机看黄色片| 亚洲精品久久成人aⅴ小说| 看黄色毛片网站| 1024视频免费在线观看| 久久久国产欧美日韩av| 免费在线观看日本一区| 亚洲国产欧美一区二区综合| 亚洲精品一区av在线观看| 91老司机精品| 免费人成视频x8x8入口观看| 此物有八面人人有两片| 九色国产91popny在线| 性色av乱码一区二区三区2| 悠悠久久av| 欧美国产日韩亚洲一区| 俄罗斯特黄特色一大片| 国产一区二区三区视频了| 美女扒开内裤让男人捅视频| 亚洲狠狠婷婷综合久久图片| 一区二区三区激情视频| 成人亚洲精品av一区二区| 欧美日韩中文字幕国产精品一区二区三区| 少妇 在线观看| 少妇的丰满在线观看| 久久香蕉国产精品| 国产亚洲av高清不卡| 欧美一级毛片孕妇| 亚洲中文av在线| 熟女电影av网| 人妻久久中文字幕网| 一个人免费在线观看的高清视频| 精品国产亚洲在线| 校园春色视频在线观看| 欧美性猛交黑人性爽| 黄网站色视频无遮挡免费观看| 人妻久久中文字幕网| 亚洲狠狠婷婷综合久久图片| 色哟哟哟哟哟哟| 欧美中文综合在线视频| 久久中文字幕一级| 久久久久久九九精品二区国产 |