亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于BADA模型的飛機(jī)持續(xù)下降進(jìn)近節(jié)能減排研究*

        2017-09-11 13:54:41吳文潔張軍峰
        關(guān)鍵詞:飛機(jī)模型研究

        吳文潔 胡 榮 張軍峰 陳 琳

        (南京航空航天大學(xué)民航學(xué)院 南京 210016)

        基于BADA模型的飛機(jī)持續(xù)下降進(jìn)近節(jié)能減排研究*

        吳文潔 胡 榮 張軍峰 陳 琳

        (南京航空航天大學(xué)民航學(xué)院 南京 210016)

        為了促進(jìn)我國民航業(yè)的節(jié)能減排,提高飛機(jī)運(yùn)行的燃油經(jīng)濟(jì)性和環(huán)保性,考慮實(shí)際氣象條件,提出了飛機(jī)的燃油消耗量和污染物排放量計(jì)算模型.結(jié)合飛機(jī)性能數(shù)據(jù),采用BADA模型仿真,開展了飛機(jī)進(jìn)近階段的節(jié)能減排研究.以廣州白云機(jī)場B767-300飛機(jī)GYA方向進(jìn)場過程為例,利用MATLAB編程計(jì)算飛機(jī)分別采用梯級進(jìn)近和持續(xù)下降進(jìn)近飛行的油耗與排放情況,結(jié)果表明,持續(xù)下降進(jìn)近可節(jié)約超過1/4的燃油消耗量,節(jié)能效果顯著;持續(xù)下降進(jìn)近可以有效減少HC,CO,NOx,CO2和SO2五種污染物的排放量,且對NOx減排效果最佳.推廣持續(xù)下降進(jìn)近的應(yīng)用有助于實(shí)現(xiàn)我國民航的節(jié)能減排.

        持續(xù)下降進(jìn)近;節(jié)能;減排;BADA模型;燃油消耗;航空尾氣排放

        0 引 言

        節(jié)能減排是我國的基本國策,隨著民航運(yùn)輸業(yè)迅速發(fā)展,飛機(jī)起降架次逐漸增多,燃油消耗與尾氣排放逐步上升,燃油價(jià)格波動(dòng)使得民航運(yùn)營成本不斷變化,飛機(jī)尾氣排放導(dǎo)致大氣環(huán)境污染日益加劇,當(dāng)今社會(huì)發(fā)展倡導(dǎo)節(jié)能減排,因此經(jīng)濟(jì)性和環(huán)保性逐漸成為民航發(fā)展的更高追求.

        飛機(jī)是民航重要組成部分,保證民航業(yè)經(jīng)濟(jì)、環(huán)保、可持續(xù)發(fā)展的基礎(chǔ)是實(shí)現(xiàn)飛機(jī)的節(jié)能減排.針對飛機(jī)節(jié)能減排問題,研究學(xué)者已開展深入研究并取得豐碩成果,基于其研究內(nèi)容,主要分為兩個(gè)方面.

        1) 飛機(jī)節(jié)能研究 由于航油價(jià)格不斷波動(dòng),合理控制飛機(jī)航油使用成為省油節(jié)能的關(guān)鍵,而燃油控制的前提則是建立精確的油耗模型.劉婧[1]提出采用飛行數(shù)據(jù)訓(xùn)練BP網(wǎng)絡(luò)計(jì)算燃油流量的模型,簡化油耗估算過程.劉芳[2]根據(jù)不同飛行環(huán)境歷史油耗數(shù)據(jù)建立下降階段的分段匹配燃油消耗模型.Collins[3]考慮飛機(jī)飛行中的動(dòng)能和勢能變化,建立基于能量平衡原理的燃油消耗估計(jì)模型.Nikoleris等[4]利用真實(shí)飛行數(shù)據(jù),對達(dá)拉斯-沃斯堡國際機(jī)場飛機(jī)的燃油消耗總量進(jìn)行估算.Singh等[5]則從宏觀方面對飛機(jī)油耗估算方法進(jìn)行了分類,分析和評價(jià).

        2) 飛機(jī)減排研究 為了維護(hù)民航經(jīng)濟(jì)發(fā)展和環(huán)境保護(hù)之間的平衡,實(shí)現(xiàn)飛機(jī)污染物排放準(zhǔn)確評估成為減排的重要基礎(chǔ).黃勇等[6-7]根據(jù)發(fā)動(dòng)機(jī)海平面靜態(tài)排放數(shù)據(jù),對中國民航飛機(jī)的NOx排放量及其分布進(jìn)行了分析和評估.孫見忠等[8]結(jié)合飛機(jī)飛行參數(shù)和性能數(shù)據(jù),利用相對法模型估算了航班污染物排放總量.夏卿等[9]基于我國民航機(jī)隊(duì)資料和飛行數(shù)據(jù),對我國123個(gè)機(jī)場1周內(nèi)飛機(jī)LTO循環(huán)的氣體污染物排放量進(jìn)行了計(jì)算.魏志強(qiáng)等[10-11]利用飛機(jī)性能軟件建立污染物排放量計(jì)算模型,實(shí)現(xiàn)不同飛行階段污染物排放量的計(jì)算.Owen等[12]則利用歷史排放數(shù)據(jù),對2050年全球航空活動(dòng)所造成的污染物排放數(shù)據(jù)進(jìn)行了大致估算.

        在追求“建設(shè)資源節(jié)約型,環(huán)境友好型社會(huì)”的形勢下,我國民航局明確提出未來五年要全面推進(jìn)民航節(jié)能減排,其中指出在安全運(yùn)行前提內(nèi),采用經(jīng)濟(jì)方式實(shí)施飛行,最大限度地減少燃油消耗和環(huán)境污染[13].持續(xù)下降進(jìn)近就是在環(huán)保要求日益升溫的背景下提出來的.它是指使飛機(jī)保持怠速狀態(tài)實(shí)現(xiàn)恒定下滑角的下降進(jìn)近,即以預(yù)設(shè)的下降梯度,從起始下降高度持續(xù)下降到最低下降高度,中間無平飛航段[14].持續(xù)下降進(jìn)近有下滑航線高、推力小、時(shí)間短和噪聲低等優(yōu)點(diǎn),不僅可以有效減少燃油消耗并降低環(huán)境污染,而且可以減輕飛行員負(fù)擔(dān)并降低事故發(fā)生率[15-16].

        盡管目前針對飛機(jī)節(jié)能減排的研究已取得一定成果,但仍略有不足:①現(xiàn)有油耗和排放的理論研究,較少考慮飛行過程中的氣象條件,預(yù)測精度相對有限,難以滿足民航運(yùn)輸?shù)膶?shí)際需要;②國內(nèi)關(guān)于持續(xù)下降進(jìn)近的研究,大多從宏觀層面開展,由具體案例分析開展的研究成果比較少,與美歐等發(fā)達(dá)國家的差距明顯.

        基于上述考慮,本文以廣州白云機(jī)場GYA方向B767-300飛機(jī)進(jìn)場過程為例,基于BADA模型,考慮氣象條件,利用MATLAB編程,研究梯級進(jìn)近和持續(xù)下降進(jìn)近的飛機(jī)油耗與排放問題,通過模型仿真和理論計(jì)算分析持續(xù)下降進(jìn)近對飛機(jī)燃油消耗和污染物排放的影響.

        1 BADA模型

        飛機(jī)性能數(shù)據(jù)庫(base of aircraft data,BADA)模型是由歐洲航空安全組織開發(fā)的分析飛行動(dòng)力學(xué)的模型,主要應(yīng)用于飛行仿真、航跡預(yù)測、排放評估和油耗計(jì)算四個(gè)方面,且在飛行仿真度、復(fù)雜性和精確度上具有較大優(yōu)勢[17-18].

        1.1 全能量模型

        反映飛機(jī)運(yùn)動(dòng)過程中速度、推力和高度等飛行參數(shù)的變化規(guī)律及約束關(guān)系,核心內(nèi)容是合外力做功等于飛機(jī)動(dòng)能和勢能的變化量.

        (1)

        (3)

        式中:T為發(fā)動(dòng)機(jī)推力;D為飛機(jī)所受阻力;VTAS為飛機(jī)真空速;m為飛機(jī)質(zhì)量;g為重力加速度;h為飛行高度;f(Ma)為能量分配系數(shù),表示沿著選定速度剖面爬升或下降時(shí),用于爬升或下降的能量占所有可用能量的比值;Ma為飛行馬赫數(shù).

        1.2 氣動(dòng)及推力模型

        飛機(jī)在飛行中會(huì)受到升力、重力、阻力及推力作用,直接影響其飛行速度,根據(jù)各飛行階段的性能系數(shù),可計(jì)算飛機(jī)的推力和阻力.

        (4)

        (5)

        (6)

        式中:γ為下降航跡角;CD0,CD1為與機(jī)型有關(guān)的阻力系數(shù);ρ為空氣密度;S為機(jī)翼面積;V為飛行速度.

        1.3 風(fēng)速修正模型

        在大氣邊界層內(nèi),風(fēng)速和高度有關(guān).氣象部門一般會(huì)提供距離地面10 m高度定時(shí)觀測的風(fēng)速資料,不同高度上的風(fēng)速采用乘冪律公式修正.

        (7)

        V=VTAScosφ+VWcosω

        (8)

        式中:VW為距地高度為h處的風(fēng)速;V0為距地高度為10 m的觀測風(fēng)速;h0為距地參考高度,取值為10 m;λ為穩(wěn)定度參數(shù),考慮機(jī)場周邊建筑特征,λ一般取值為0.10;φ為飛機(jī)的航向角;ω為風(fēng)向角.

        1.4 油耗模型

        單位推力的燃油消耗率

        (9)

        (10)

        燃油流量f=ηCfT

        (11)

        式中:η為單位推力的燃油消耗率;Cf1,Cf2為油耗參數(shù);f為燃油流量;Cf為推力等級.

        2 油耗與排放定量計(jì)算

        2.1 燃油消耗量計(jì)算模型

        考慮實(shí)際氣象因素,不同的氣溫、氣壓和風(fēng)速條件都會(huì)對飛機(jī)燃油消耗產(chǎn)生影響,可利用參數(shù)修正法對實(shí)際燃油消耗進(jìn)行修正.

        標(biāo)準(zhǔn)燃油流量模型

        (12)

        (13)

        2.2 污染物排放量計(jì)算模型

        根據(jù)飛機(jī)-發(fā)動(dòng)機(jī)型號組合,參照海平面靜態(tài)排放數(shù)據(jù),結(jié)合氣象因素對各污染物排放指數(shù)進(jìn)行修正,可計(jì)算航空氣體污染物的排放量.

        (14)

        (15)

        (16)

        3 油耗與排放結(jié)果分析

        本文以廣州白云機(jī)場GYA方向進(jìn)場過程為例,開展梯級進(jìn)近和持續(xù)下降進(jìn)近的飛機(jī)油耗與排放評估研究.根據(jù)飛行數(shù)據(jù)可知,GYA方向的進(jìn)場高度是從3 900 m下降至1 800 m.圖1為GYA方向三維飛行軌跡圖,即梯級進(jìn)近與持續(xù)下降進(jìn)近的空間航跡示意圖,其中X坐標(biāo)指“經(jīng)度”,Y坐標(biāo)指“緯度”,Z坐標(biāo)指“飛行高度”.

        圖1 GYA方向三維飛行軌跡圖

        考慮白云機(jī)場繁忙程度且持續(xù)下降進(jìn)近安排在夜間運(yùn)行,對白云機(jī)場2015年8月份00:00—03:00所有航班運(yùn)行情況進(jìn)行統(tǒng)計(jì),結(jié)果表明,B767-300飛機(jī)的運(yùn)行頻次相對較高,故選用 B767-300為本文研究機(jī)型.B767-300飛機(jī)配備的發(fā)動(dòng)機(jī)型號為PW4056,由《Aircraft Engine Exhaust Emissions Data Bank》可知其污染物排放指數(shù),見表1.

        表1 PW4056的污染物排放數(shù)據(jù)

        對于污染物CO2與SO2而言,其排放指數(shù)和發(fā)動(dòng)機(jī)型號及所使用的航空燃油有關(guān),而與飛行階段無關(guān),查閱污染物排放標(biāo)準(zhǔn)可知,PW4056發(fā)動(dòng)機(jī)的CO2排放指數(shù)為3.104 kg/kg,SO2排放指數(shù)為1 g/kg.

        3.1 持續(xù)下降進(jìn)近節(jié)能效果評估

        節(jié)約燃油是在確保飛行安全前提下成本控制首要考慮的因素,基于BADA模型仿真,利用MATLAB計(jì)算B767-300飛機(jī)分別采用梯級進(jìn)近和持續(xù)下降進(jìn)近的燃油消耗量為202.80 kg和149.60 kg,相比梯級進(jìn)近,持續(xù)下降進(jìn)近可減少燃油消耗53.20 kg,節(jié)約航空燃油26.23%.

        3.2 持續(xù)下降進(jìn)近減排效果評估

        一般來說節(jié)能意味著減排,減少燃油消耗的同時(shí)可在一定程度上減輕飛機(jī)排放污染.飛機(jī)排放的污染物主要包括HC,CO,NOx,CO2和SO2等,這些污染物不僅會(huì)對機(jī)場周邊環(huán)境與公共健康造成嚴(yán)重危害,而且會(huì)導(dǎo)致全球溫室效應(yīng)和氣候變化.根據(jù)修正污染物排放計(jì)算模型,利用MATLAB分別估算B767-300飛機(jī)采用梯級進(jìn)近和持續(xù)下降進(jìn)近的五種污染物排放量,具體排放情況見表2.

        表2 五種污染物排放量

        注:減排量=梯級進(jìn)近排放量-持續(xù)下降進(jìn)近排放量;減排比例=減排量/梯級進(jìn)近排放量.

        由表2可知,相比梯級進(jìn)近,持續(xù)下降進(jìn)近能在不同程度上減少飛機(jī)污染物排放并減輕大氣環(huán)境污染,減排效果明顯.

        3.3 持續(xù)下降進(jìn)近節(jié)能減排原因分析

        3.3.1 節(jié)能原因分析

        由燃油消耗量的計(jì)算原理可知,影響油耗的兩大因素是燃油流量和飛行時(shí)間,圖2為燃油流量和飛行時(shí)間對比圖.

        圖2 燃油流量和飛行時(shí)間對比圖

        由圖2a)可知,在GYA方向進(jìn)場飛行過程中,采用持續(xù)下降進(jìn)近飛行的燃油流量始終低于梯級進(jìn)近.由圖2b)可知,相比梯級進(jìn)近,持續(xù)下降進(jìn)近可減少飛行時(shí)間9 s,可節(jié)省時(shí)間2.13%.

        因此綜合燃油流量和飛行時(shí)間考慮,持續(xù)下降進(jìn)近能明顯減少燃油消耗.

        3.3.2 減排原因分析

        由污染物排放量計(jì)算模型可知,影響飛機(jī)排放的兩大因素是燃油消耗量和排放指數(shù).其中,燃油消耗量主要和燃油流量和飛行時(shí)間有關(guān).

        而排放指數(shù)則受外界氣象條件如濕度因子、風(fēng)速、風(fēng)向、氣溫和氣壓等影響,這些影響因素又與飛行高度有關(guān),所以影響飛機(jī)排放的一個(gè)關(guān)鍵因素就是飛行高度.根據(jù)飛行數(shù)據(jù)可知,GYA方向進(jìn)場距離為62 km,圖3為GYA方向飛行剖面示意圖,其中橫坐標(biāo)“飛行距離”是指飛機(jī)到最后進(jìn)近定位點(diǎn)的水平距離.

        圖3 GYA方向飛行剖面圖

        由圖3可知,持續(xù)下降進(jìn)近的飛行高度整體高于梯級進(jìn)近.一般情況下,海拔高度每上升100 m,氣溫平均下降約0.6 ℃,氣壓平均下降約1 000 Pa.相比梯級進(jìn)近,持續(xù)下降進(jìn)近的飛行高度升高,濕度因子增大,外界氣溫和氣壓降低,導(dǎo)致HC,CO和NOx的排放指數(shù)減小,雖然CO2和SO2的排放指數(shù)與飛行高度無關(guān),但是持續(xù)下降進(jìn)近的燃油消耗量明顯低于梯級進(jìn)近,因此綜合燃油消耗量和排放指數(shù)考慮,持續(xù)下降進(jìn)近能有效降低各種污染物排放.

        4 結(jié) 論

        1) 相比梯級進(jìn)近,持續(xù)下降進(jìn)近的發(fā)動(dòng)機(jī)推力小,燃油流量低,飛行時(shí)間短,可節(jié)約超過1/4的燃油.

        2) 相比梯級進(jìn)近,持續(xù)下降進(jìn)近的飛行高度高,污染物排放指數(shù)小,因此持續(xù)下降進(jìn)近對HC,CO和NOx的減排效果明顯,且對NOx的減排效果尤為顯著.

        3) 由于CO2和SO2排放指數(shù)不變,其排放量與油耗成正比,因此持續(xù)下降進(jìn)近對CO2和SO2的減排效果相同.

        4) 持續(xù)下降進(jìn)近節(jié)能減排的主要原因是其推力較小、飛行高度較高和飛行時(shí)間較短.

        然而本文僅從理論計(jì)算和模型仿真的層面研究持續(xù)下降進(jìn)近節(jié)能減排效果,對真機(jī)試驗(yàn)的實(shí)證分析仍需進(jìn)一步研究.此外本文僅研究飛機(jī)自身的節(jié)能減排問題,對于民航業(yè)的節(jié)能減排,未來可以考慮從技術(shù)、運(yùn)營和市場三個(gè)宏觀角度著手研究,比如采用新型發(fā)動(dòng)機(jī)和開發(fā)可替代航空燃油、提高空管和機(jī)場運(yùn)營保障能力以及采用碳抵消和碳稅費(fèi)等經(jīng)濟(jì)措施來倡導(dǎo)民航節(jié)能減排,實(shí)現(xiàn)我國綠色民航發(fā)展.

        [1]劉婧.基于飛行數(shù)據(jù)分析的飛機(jī)燃油估計(jì)模型[D].南京:南京航空航天大學(xué),2010.

        [2]劉芳.下降階段的飛機(jī)油耗優(yōu)化建模方法研究[D].南京:南京航空航天大學(xué),2012.

        [3]COLLINS B P. Estimation of aircraft fuel consumption[J]. Journal of Aircraft, 2015, 19(11):969-975.

        [4]NIKOLERIS T, GUPTA G, KISTLER M. Detailed estimation of fuel consumption and emissions during aircraft taxi operations at Dallas/Fort worth international airport[J]. Transportation Research Part D Transport & Environment, 2011,16(4):302-308.

        [5]SINGH V, SHARMA S K. Fuel consumption optimization in air transport: a review, classification, critique, simple meta-analysis, and future research implications[J]. European Transport Research Review, 2015,7(2):1-24.

        [6]黃勇,周桂林,吳壽生.中國上空民航飛機(jī)NOx排放量及其分布初探[J].環(huán)境科學(xué)學(xué)報(bào),2000,20(2):179-182.

        [7]黃勇,吳冬鶯,王金濤,等.中國上空民航飛機(jī)NOx排放分布再探[J].北京航空航天大學(xué)學(xué)報(bào),2001,27(3):289-292.

        [8]孫見忠,左洪福,劉鵬鵬,等.航空發(fā)動(dòng)機(jī)污染物排放量估算方法[J].交通運(yùn)輸工程學(xué)報(bào),2012,12(2):53-61.

        [9]夏卿,左洪福,楊軍利.中國民航機(jī)場飛機(jī)起飛著陸(LTO)循環(huán)排放量估算[J].環(huán)境科學(xué)學(xué)報(bào),2008,28(7):1469-1474.

        [10]魏志強(qiáng),刁華智,韓博.民用飛機(jī)巡航階段污染物排放量計(jì)算研究[J].科學(xué)技術(shù)與工程,2014,14(19):122-127.

        [11] 魏志強(qiáng),王超.航班飛行各階段污染物排放量估算方法[J].交通運(yùn)輸工程學(xué)報(bào),2010,10(6):48-52.

        [12]OWEN B, LEE D S, LIM L. Flying into the future: aviation emissions scenarios to 2050[J]. Environmental Science & Technology, 2010, 44(7):255-260.

        [13]劉平.全面推進(jìn)民航節(jié)能減排[J].中國民用航空,2011,12(6):1-9.

        [14]孫鵬.民用飛機(jī)連續(xù)下降技術(shù)初探[J].民用飛機(jī)設(shè)計(jì)與研究,2012,22(1):6-9.

        [15]JIN L, CAO Y, SUN D F. Investigation of potential fuel savings due to continuous-descent approach[J]. Journal of Aircraft, 2013,50(3):807-816.

        [16]ALAM S, NGUYEN M H, ABBASS H A. Multi-aircraft dynamic continuous descent approach methodology for low-noise and emission guidance[J]. Journal of Aircraft, 2011, 48(4):1225-1237.

        [17]張軍峰,蔣海行,武曉光,等.基于BADA及航空器意圖的四維航跡預(yù)測[J].西南交通大學(xué)學(xué)報(bào),2014,49(3):553-558.

        [18]閆國華,魏娜,楊曉軍.基于BADA模型的我國民用飛機(jī)持續(xù)進(jìn)近燃油經(jīng)濟(jì)性研究[C].第六屆國際節(jié)能與新能源汽車創(chuàng)新發(fā)展論壇,深圳,2009.

        Research on Aircraft Energy-saving and Emission-reduction of Continuous Descent Approach Based on BADA model

        WU Wenjie HU Rong ZHANG Junfeng CHEN Lin

        (CollegeofCivilAviation,NanjingUniversityofAeronauticsandAstronautics,Nanjing210016,China)

        In order to promote the energy-saving and emission-reduction in China’s civil aviation and improve the fuel economy and environmental protection of aircraft operation, the fuel consumption model and pollutants emission model are presented according to the actual weather conditions. Combined with the aircraft performance data and based on the BADA simulation model, the research on aircraft energy-saving and emission-reduction in approach stage is conducted. Using the example of GYA arrival point of Guangzhou-Baiyun airport and B767-300 aircraft, a MATLAB programming software is used to estimate the fuel consumption and pollutants emission by step-down approach and Continuous Descent Approach, respectively. The results show that: compared with the step-down approach, the Continuous Descent Approach can save over a quarter of the fuel consumption per flight and it has significant energy-saving effect; Continuous Descent Approach also has different emission-reduction effect on five pollutants including HC, CO, NOx, CO2and SO2, and the effect for NOxis more significant. The promotion of Continuous Descent Approach is found to be helpful to reaching the objective of energy-saving and emission-reduction in China’s civil aviation.

        continuous descent approach; energy-saving; emission-reduction; BADA model; fuel consumption; exhaust gas emission

        2017-02-08

        *國家自然科學(xué)基金項(xiàng)目資助(71401072,71201082)

        V249.1

        10.3963/j.issn.2095-3844.2017.04.027

        吳文潔(1993—):女,碩士生,主要研究領(lǐng)域?yàn)榫G色民航,交通運(yùn)輸規(guī)劃與管理

        猜你喜歡
        飛機(jī)模型研究
        一半模型
        飛機(jī)失蹤
        FMS與YBT相關(guān)性的實(shí)證研究
        遼代千人邑研究述論
        重要模型『一線三等角』
        重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
        視錯(cuò)覺在平面設(shè)計(jì)中的應(yīng)用與研究
        科技傳播(2019年22期)2020-01-14 03:06:54
        “拼座飛機(jī)”迎風(fēng)飛揚(yáng)
        EMA伺服控制系統(tǒng)研究
        乘坐飛機(jī)
        一本一道波多野结衣一区| 街拍丝袜美腿美女一区| 蜜桃视频羞羞在线观看| 蜜臀一区二区三区精品| 永久免费人禽av在线观看| 中文字幕在线观看亚洲日韩| 亚洲综合av在线在线播放| 国产福利片无码区在线观看 | 白白色发布在线播放国产| 亚洲人av毛片一区二区| 中文字幕日本av网站| 在线视频观看国产色网| 无码人妻少妇久久中文字幕蜜桃| 男女啪啪免费体验区| 精品国产一级毛片大全| 日韩久久无码免费看A| 清纯唯美亚洲经典中文字幕| 男女裸体做爰视频高清| 色先锋av影音先锋在线| 亚洲色无码播放| 欧美中出在线| 亚洲福利二区三区四区| 国产午夜福利片在线观看| 精品国产乱码久久久久久影片| 国内精品久久久久久久久齐齐| 91精品国产综合久久青草| 日韩中文字幕乱码在线| 国产毛片av一区二区| 妇女bbbb插插插视频| 国产免费又色又爽又黄软件| 国产精品白浆视频免费观看| 少妇高潮太爽了免费网站| 中文字幕女优av在线| 人妻夜夜爽天天爽三区 | 久久se精品一区精品二区国产| av网站入口在线免费观看| 日本不卡视频一区二区三区| 国产精品 无码专区| 久久av高潮av无码av喷吹| 亚洲春色AV无码专区在线播放| 久久夜色精品国产九色|