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

    基于改進(jìn)NSGA?Ⅱ算法的航班戰(zhàn)略沖突解脫研究

    2022-12-25 07:47:02胡明華
    關(guān)鍵詞:高度層父代航空器

    徐 滿,胡明華,張 穎,江 灝

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

    隨著經(jīng)濟(jì)的不斷發(fā)展,亞太地區(qū)將成為航空器客運(yùn)量增長的主要?jiǎng)恿?,中國將?022 年成為全球第一大航空運(yùn)輸市場[1]。由于空中交通需求的不斷增長,當(dāng)前的空域會(huì)變得越來越擁擠,空中交通的擁堵情況越來越嚴(yán)重。

    航空器的沖突解脫問題是一個(gè)大規(guī)模組合優(yōu)化問題,且由于航空器之間沖突的存在使得航空器之間相互耦合[2]。有大量的文獻(xiàn)研究航空器在戰(zhàn)術(shù)及預(yù)戰(zhàn)術(shù)階段的沖突解脫問題,主要通過調(diào)整航空器的起飛時(shí)刻、調(diào)整航空器的飛行速度、修改航空器的水平航跡和分配高度層等方式解決相應(yīng)的沖突[3?4]。但是,采用這些方式的短期內(nèi)的沖突解脫存在一定的缺點(diǎn),即沒有從全局角度考慮航空器的沖突解脫。由于航空器之間的相互關(guān)系,短期內(nèi)的沖突解脫可能會(huì)產(chǎn)生“多米諾骨牌”效應(yīng),使得沖突解脫一直處于被動(dòng),危及空域安全。例如,采用地面等待策略解決沖突時(shí),已經(jīng)延誤了的航空器仍然需要繼續(xù)等待其他航空器以滿足空域的容量要求[5]。而戰(zhàn)略沖突解脫研究,旨在根據(jù)航空器的飛行計(jì)劃在起飛前相當(dāng)長一段時(shí)間內(nèi)解決潛在的飛行沖突,是確保航空器安全運(yùn)行的一項(xiàng)關(guān)鍵技術(shù)。目前,戰(zhàn)略沖突解脫問題引起了大量的關(guān)注。

    四維航跡的運(yùn)用以及基于航跡運(yùn)行概念的提出被認(rèn)為是未來空中交通管理的運(yùn)行基礎(chǔ),在此基礎(chǔ)上戰(zhàn)略沖突解脫問題得到了進(jìn)一步的研究。四維航跡將航空器航跡定義為三維空間加上時(shí)間,隨著運(yùn)行技術(shù)的不斷發(fā)展,四維航跡可以有效降低航空器運(yùn)行過程中的不確定性。這就為在一個(gè)較長時(shí)間范圍內(nèi)解決航空器之間潛在沖突提供了技術(shù)支持。

    已有大量的文獻(xiàn)研究了航空器四維航跡戰(zhàn)略沖突解脫問題。沖突解脫由于問題本身的復(fù)雜度較高,常采用啟發(fā)式算法求解,并且在提高算法的求解效率上做了大量的研究。2019 年,Dai 等[6]提出了一種改進(jìn)的模擬退火算法,將求解過程分成兩個(gè)階段:第一階段用于減少?zèng)_突數(shù)量;第二階段用于優(yōu)化航空器的飛行時(shí)間。同時(shí),在沖突解脫時(shí)設(shè)計(jì)了一個(gè)均勻高度層策略,引導(dǎo)航空器均勻分布至各個(gè)可用高度層,引入滑動(dòng)時(shí)間窗來降低問題復(fù)雜度以提高求解效率。但是其將沖突數(shù)量與航空器航跡偏移量進(jìn)行線性組合,解決的仍是單目標(biāo)優(yōu)化問題。除此之外,還有大量運(yùn)用模擬退火算法及其改進(jìn)算法的研究[7?9]。合作協(xié)同進(jìn)化算法[10]在求解航空器沖突解脫問題也擁有良好的性能,合作協(xié)同進(jìn)化算法采用分而治之的思想,將一個(gè)大規(guī)模組合優(yōu)化問題分解成若干個(gè)子問題,從而降低問題的維數(shù),各個(gè)子問題采用特定的進(jìn)化算法優(yōu)化得到部分可行解,通過各個(gè)小組之間的合作進(jìn)化獲取全局可行解。在合作協(xié)同進(jìn)化算法的框架下,采用遺傳算法進(jìn)行優(yōu)化可以極大提高算法的優(yōu)化效率,其研究的仍是單目標(biāo)優(yōu)化問題。

    在實(shí)際運(yùn)行過程中,空中交通管制員往往會(huì)在沖突數(shù)量和航空器航跡調(diào)整量之間尋求一個(gè)很好的平衡,將沖突解脫問題建模為多目標(biāo)優(yōu)化問題更加 符 合 實(shí) 際。Su 等[11]于2017 年 采 用 多 目 標(biāo) 文 化基因算法求解沖突解脫問題。算法采用經(jīng)典非支配排序遺傳算法(Non?dominated sorting genetic al?gorithm Ⅱ,NSGA?Ⅱ)進(jìn)行全局搜索獲取可行解,同時(shí)結(jié)合局部搜索算法加強(qiáng)對當(dāng)前最優(yōu)解的搜索,獲得更優(yōu)可行解,但是算法的復(fù)雜度較高。多目標(biāo)合作協(xié)同進(jìn)化算法[12]同樣被應(yīng)用于沖突解脫問題,在協(xié)同進(jìn)化算法的框架下,采用基于分解的多目標(biāo)進(jìn)化算法(Multi?objective evolutionary algo?rithm based on decomposition,MOEA/D)[13]進(jìn) 行各個(gè)小組的內(nèi)部優(yōu)化。雖然采用了合作協(xié)同進(jìn)化算法的框架,但是算法的求解效率較低。

    針對以上問題,本文采用一種改進(jìn)的NSGA?Ⅱ算法研究航空器沖突解脫中的多目標(biāo)優(yōu)化問題,通過對算法內(nèi)部遺傳算子的設(shè)計(jì),提高算法優(yōu)化效率,同時(shí)控制算法的復(fù)雜度。在解決航空器之間潛在沖突的同時(shí)考慮與航空器航跡調(diào)整量所引起的航跡調(diào)整成本。本文采用調(diào)整航空器起飛時(shí)刻和分配高度層的方式在戰(zhàn)略階段解決航空器之間的潛在沖突,沖突成本來自航班沖突解脫后的飛行計(jì)劃與初始飛行計(jì)劃之間的偏差。主要?jiǎng)?chuàng)新之處在于,本文為染色體中的基因賦予相應(yīng)的局部適應(yīng)度,設(shè)計(jì)了一種自適應(yīng)交叉算子與變異算子,結(jié)合NSGA?Ⅱ框架進(jìn)行全局優(yōu)化。

    1 問題建模

    1.1 沖突探測模型

    本文研究巡航狀態(tài)下的航空器沖突解脫問題,因此假設(shè)航跡起終點(diǎn)為爬升和下降頂點(diǎn)。假設(shè)空域中有n架航空器,F(xiàn)=(f1,f2,f3,…,fn)按照航班計(jì)劃沿著特定的航路飛行,航空器處于巡航狀態(tài),按照固定的高度層以恒定的速度運(yùn)行。航空器的安全間隔在垂直方向上為Nv,水平方向?yàn)镹h,當(dāng)航空器在空域中飛行時(shí),它們的航跡相互影響。假設(shè)航班嚴(yán)格地按照計(jì)劃航線飛行,那么潛在的沖突會(huì)發(fā)生在兩條航路的交叉點(diǎn)附近,如圖1 所示。當(dāng)兩架航空器的保護(hù)區(qū)發(fā)生重疊時(shí),則認(rèn)為兩架航空器存在潛在沖突,這種情況并不意味著兩架航空器一定會(huì)沖突,但是在運(yùn)行過程中是需要避免出現(xiàn)的。當(dāng)航空器fi和航空器fj的水平距離小于Nh且垂直距 離 小 于Nv時(shí),‖Pfi(t)-Pfj(t)‖2

    圖1 航空器沖突圖Fig.1 Aircraft conflict

    為了表示航空器之間的沖突關(guān)系,本文定義了一個(gè)矩陣C來表示航空器之間是否有沖突,即

    1.2 數(shù)學(xué)優(yōu)化模型

    航路網(wǎng)絡(luò)范圍內(nèi)的戰(zhàn)略沖突解脫問題可以描述為:給定一個(gè)航路網(wǎng)絡(luò)范圍,對于航路網(wǎng)絡(luò)中的每一架航空器進(jìn)行四維航跡規(guī)劃,在較低的航空器沖突解脫成本下盡可能減少航空器之間的沖突數(shù)量。在進(jìn)行沖突解脫的過程中,應(yīng)滿足空域容量限制、航空器機(jī)動(dòng)性能限制等。航空器f的航跡由一組四維航跡采樣點(diǎn)組成,即有

    1.2.2 目標(biāo)函數(shù)

    為保證航空器運(yùn)行的安全性,必須盡可能地減少航空器之間的潛在沖突,由式(1)中的矩陣C可得航班之間的沖突關(guān)系。那么,航班fi的沖突數(shù)量可以由式(3)計(jì)算。

    每一個(gè)沖突都被一對沖突航班重復(fù)計(jì)算兩次,所以航班間的沖突數(shù)量為

    在求解過程中,為了盡可能避免與其他航班之間的潛在沖突,可能對航班的初始航跡做出重大修改,從而導(dǎo)致大量偏移其初始航班計(jì)劃。但是,這樣的方式可能會(huì)產(chǎn)生額外的成本。因此,沖突解脫還應(yīng)從經(jīng)濟(jì)角度進(jìn)行評估,即航班應(yīng)盡可能小地偏離初始計(jì)劃,目標(biāo)2 重點(diǎn)關(guān)注最小化航跡調(diào)整成本(Total trajectory amendment cost,TTAC)。本文采用調(diào)整航空器起飛時(shí)刻以及高度層分配解決航空器之間的潛在沖突,TTAC 由地面延誤成本(CostGH)和高度層分配成本(CostFL)構(gòu)成,那么目標(biāo)2 可由式(5)計(jì)算。

    式中:λGH為地面延誤成本系數(shù);λFL為高度層分配成本系數(shù)。由于地面延誤成本與高度層分配成本單位不同,無法直接相加,因此需要進(jìn)行歸一化處理,即將每個(gè)航班的地面延誤和高度調(diào)整量除以最大允許值。

    (1)地面延誤成本

    航空器實(shí)際離場時(shí)間與初始飛行計(jì)劃中離場時(shí)間之差為航空器的地面延誤,即

    式中l(wèi)orig為航空器計(jì)劃飛行高度層。

    1.2.3 約束條件

    考慮到安全性、公平性和航空器的性能,軌跡中的每個(gè)元素應(yīng)盡可能接近其計(jì)劃值。模型的約束集如下所示

    ?f∈F,e∈E,l∈L,t∈T均 應(yīng) 滿 足 上 述 約 束,E是航段集合。其中,約束(8)保證航空器沿著單一可行航線飛行;約束(9)保證航空器的地面延誤不超出最大允許延誤,TfGH表示單個(gè)航空器最大允許延誤;約束(10)保證航空器的飛行高度層只能在允許的高度層內(nèi),NfFL為最大高度調(diào)整量。

    2 多目標(biāo)優(yōu)化算法設(shè)計(jì)

    2.1 改進(jìn)的NSGA?Ⅱ算法

    多目標(biāo)優(yōu)化可以定義為尋找滿足約束條件的決策變量向量,并同時(shí)優(yōu)化多個(gè)目標(biāo)函數(shù)的問題,目標(biāo)函數(shù)之間往往是對立的。NSGA?Ⅱ根據(jù)目標(biāo)函數(shù)特點(diǎn),對可行解進(jìn)行非支配排序,最終得到Pareto 前沿。與經(jīng)典的NSGA?Ⅱ算法及MOEA/D 相比,本文采用的算法根據(jù)問題特點(diǎn)設(shè)計(jì)了自適應(yīng)交叉算子與變異算子[15]加快進(jìn)化過程。

    2.1.1 初始化種群

    本文采用調(diào)整航空器起飛時(shí)刻與高度層分配兩種方式解決航空器間的潛在沖突,在生成初始種群時(shí),根據(jù)當(dāng)前航班計(jì)劃中各個(gè)航班沖突數(shù)量,采用輪盤賭隨機(jī)選擇一條航跡,沖突數(shù)量越多,航跡被選中的概率越大。以概率形式隨機(jī)修改起飛時(shí)刻與飛行高度層,獨(dú)立運(yùn)行N次,生成初始種群,具體過程如圖2 所示。

    圖2 種群初始化流程圖Fig.2 Flow chart of population initialization

    2.1.2 選擇

    本文采用錦標(biāo)賽選擇法。從父代種群中隨機(jī)選擇兩個(gè)父代,首先比較兩個(gè)父代的非支配序,優(yōu)先選擇非支配序較低的父代個(gè)體;若父代的非支配序相同,則比較父代個(gè)體的擁擠距離,優(yōu)先選擇擁擠距離較大的父代個(gè)體。

    2.1.3 交叉

    本文采用的自適應(yīng)交叉算子,在進(jìn)行交叉操作時(shí),對于選擇的父代個(gè)體,比較父代個(gè)體相同位置基因的適應(yīng)度,即在兩個(gè)不同的航班計(jì)劃中比較同一架航班在該計(jì)劃下的沖突解脫成本。個(gè)體中每架航班的局部適應(yīng)度由式(11)計(jì)算。

    式中:AC和BC分別為線性重組生成的兩個(gè)子代,floor 表示向下取整,k為線形重組系數(shù),交叉概率為Pc。自適應(yīng)交叉算子如圖3 所示。

    圖3 自適應(yīng)交叉算子Fig.3 Adaptive crossover operator

    2.1.4 變異

    較高的局部適應(yīng)度意味著該航班的飛行計(jì)劃產(chǎn)生較少的沖突和較低的飛行成本,所以當(dāng)染色體某個(gè)位置基因適應(yīng)度值小于ε時(shí),才發(fā)生變異,其發(fā)生變異的概率為Pm。自適應(yīng)變異算子如圖4所示。

    圖4 自適應(yīng)變異算子Fig.4 Adaptive mutation operator

    2.2 多目標(biāo)優(yōu)化解的評價(jià)

    如何評價(jià)沖突解脫優(yōu)化解集的優(yōu)劣,在使航班盡可能小地偏離初始飛行計(jì)劃的情況下得到較少?zèng)_突數(shù)量的航班計(jì)劃是一個(gè)值得深究的問題。本文擬從Pareto 最優(yōu)解的收斂性和多樣性兩個(gè)方面進(jìn)行評估,采用以下4 個(gè)指標(biāo)[16]來評估本文提出的算法。

    2.2.1 覆蓋率

    覆蓋率(C?Metric,CM)是衡量Pareto 前沿的收斂性的指標(biāo),即

    式中PFU、PFV分別為兩組Pareto 最優(yōu)解集。分子表示解集V中被U中至少一個(gè)解支配的解的數(shù)目,分母表示解集V中所有解的個(gè)數(shù)。 若CCM(PFU,PFV)=1,說明解集V完全被解集U支配,即解集V收斂性比解集U差。

    2.2.2 空間分布

    空間分布(Spacing,SP)是衡量Pareto 前沿的廣泛程度的指標(biāo),即

    式中:up表示各個(gè)解到其他解的最小歐氏距離;-u為均值。CSP(PF)值越小,說明空間分布越均勻。

    2.2.3 超體積

    超體積(Hyper volume,HV)是衡量Pareto 前沿收斂性和多樣性的綜合指標(biāo),即

    式中:δ代表Lebesgue 測度,用來衡量體積;Vp表示第p個(gè)點(diǎn)與參考點(diǎn)圍成區(qū)域的超體積。超體積值越大,表示解集的效果越好。

    2.2.4 平均理想距離

    平均理想距離(Mean ideal distance,MID)是衡量Pareto 前沿和理想點(diǎn)間距離的指標(biāo),即

    式中表示第p個(gè)非支配解的第t個(gè)目標(biāo)值。本文研究雙目標(biāo)優(yōu)化問題,因?yàn)榍蠼庾钚』瘑栴},所以理想點(diǎn)可設(shè)置為坐標(biāo)原點(diǎn)。

    3 仿真驗(yàn)證

    3.1 運(yùn)行場景

    本文重點(diǎn)研究巡航階段沖突解脫問題,數(shù)據(jù)采用2019 年6 月8 日9:00 至12:00 中 國 航 路 網(wǎng) 絡(luò) 中1 472 架航班進(jìn)行仿真驗(yàn)證,在MATLAB R2019a平臺(tái)上運(yùn)行,航班飛過的航路點(diǎn)如圖5 所示,涉及2 136 個(gè)航路點(diǎn)和186 個(gè)起降機(jī)場。本文對每條航跡采用等時(shí)間間隔采樣生成相應(yīng)的航跡采樣點(diǎn),采樣間隔為15 s,利用墨卡托投影將航跡采樣點(diǎn)從經(jīng)緯度坐標(biāo)轉(zhuǎn)化為(x,y)坐標(biāo)。允許的單個(gè)航班最大延誤TfGH為60 min,航班離場時(shí)隙為1 min;單個(gè)航班飛行高度最大調(diào)整量NfFL為1 200 m,遵循“東單西雙”的原則,允許的飛行高度范圍為8 900 m至11 600 m,以300 m 的間隔劃分高度層。采用網(wǎng)格法探測航空器之間的潛在沖突,網(wǎng)格在水平方向上大小為Nh,在垂直方向上為Nv,時(shí)間維度的網(wǎng)格大小與采樣時(shí)間間隔相等。本文采用的改進(jìn)NS?GA?Ⅱ算法各項(xiàng)參數(shù)如表1 所示,其中,種群規(guī)模與進(jìn)化代數(shù)等參數(shù)均通過多次對比實(shí)驗(yàn)得到的最優(yōu)結(jié)果確定。

    圖5 航路點(diǎn)圖Fig.5 Distribution of waypoints

    表1 算法參數(shù)設(shè)置Table 1 Parameter settings of algorithms

    3.2 結(jié)果分析

    3.2.1 算法性能評價(jià)

    本文將改進(jìn)的NSGA?Ⅱ算法與經(jīng)典NSGA?Ⅱ算法以及MOEA/D 進(jìn)行對比驗(yàn)證算法的有效性,同時(shí)利用多目標(biāo)優(yōu)化解的指標(biāo)進(jìn)行評價(jià),具體結(jié)果如表2 所示。啟發(fā)式算法的優(yōu)化效果具有隨機(jī)性,因此,表2 中的結(jié)果均基于15 次獨(dú)立重復(fù)試驗(yàn)得到。

    表2 算法性能對比Table 2 Performance comparison between algorithms

    由表2 可知,改進(jìn)的NSGA?Ⅱ算法在各項(xiàng)性能上均優(yōu)于經(jīng)典的NSGA?Ⅱ算法及MOEA/D。覆蓋率指標(biāo),改進(jìn)的NSGA?Ⅱ算法得到的Pareto前沿完全支配另外兩種算法。改進(jìn)的NSGA?Ⅱ算法的Pareto 最優(yōu)解相對于經(jīng)典算法及MOEA/D而言,解集更占有支配地位(CM),解集的分布更均勻(SP),解集的收斂性更好(HV),解得質(zhì)量也更高(MID)。

    算法的收斂性和多樣性可由反轉(zhuǎn)世代距離(Inverted generational distance,IGD)反映,本文將15 次獨(dú)立運(yùn)行得到的Pareto 最優(yōu)解集P*作為參考集,計(jì)算參考點(diǎn)到最近的解的距離的平均值。IGD值越小,說明算法綜合性越好,即

    IGD 值與進(jìn)化代數(shù)的關(guān)系如圖6 所示,由于IGD 值在前100 代平穩(wěn),所以圖中只選取了前100代的結(jié)果并做歸一化處理,以使變化趨勢更加明顯。

    圖6 IGD 變化趨勢Fig.6 IGD changing trend

    由圖6 可知,在當(dāng)前參數(shù)條件下,3 種算法均可收斂。在3 種算法中,采用了自適應(yīng)遺傳算子的改進(jìn)NSGA?Ⅱ算法收斂速度最快。表3 給出了3 種算法在達(dá)到收斂代數(shù)和達(dá)到最大進(jìn)化代數(shù)時(shí)的運(yùn)行時(shí)間。對比MOEA/D 算法、經(jīng)典NSGA?Ⅱ算法與改進(jìn)的NSGA?Ⅱ算法,本文所提算法在運(yùn)算速度上提高了約8%。

    表3 算法運(yùn)行時(shí)間對比Table 3 Running time comparison between algorithms

    3.2.2 算法優(yōu)化結(jié)果

    圖7 是3 種算法的優(yōu)化結(jié)果對比圖。由圖可知,改進(jìn)的NSGA?Ⅱ算法不僅在沖突解脫的效果上優(yōu)于兩種算法(Obj1),而且產(chǎn)生的航跡偏移量(Obj2)也更小。在最終的Pareto 最優(yōu)解集中,本文的改進(jìn)算法目標(biāo)1 最優(yōu)時(shí)可以將沖突數(shù)量減少至272,而經(jīng)典的算法沖突數(shù)量則多達(dá)422,MOEA/D沖突數(shù)量最小也為324 個(gè),這反映了在引入改進(jìn)了的交叉算子和變異算子后算法在沖突解脫方面的優(yōu)異性能。目標(biāo)2 則反映出改進(jìn)后的NSGA?Ⅱ算法,在解決沖突的時(shí)候產(chǎn)生的航跡調(diào)整量更少,使航空器在沖突解脫時(shí)產(chǎn)生更小的調(diào)整成本,符合航空公司的利益。

    圖7 Pareto 前沿對比圖Fig.7 Comparison chart of Pareto front

    4 結(jié) 論

    本文主要研究了沖突解脫中的多目標(biāo)優(yōu)化問題。針對問題的特點(diǎn),本文將航空器沖突解脫建模為以沖突數(shù)量和航空器航跡調(diào)整量為目標(biāo)的雙目標(biāo)優(yōu)化問題。針對所求解的沖突解脫問題的特點(diǎn),采用了一種改進(jìn)的NSGA?Ⅱ算法,設(shè)計(jì)了自適應(yīng)交叉算子與變異算子,旨在提高算法搜索最優(yōu)解的效率。與經(jīng)典的NSGA?Ⅱ算法以及MOEA/D 相比,改進(jìn)了的算法具有沖突解脫的有效性和求解的高效性,不僅可以同時(shí)兼顧航跡調(diào)整成本解決大量沖突,而且能更快收斂,有很強(qiáng)的參考價(jià)值。

    下一步的工作主要集中對更復(fù)雜沖突解脫方式的建模以及進(jìn)一步提高算法的效率上:沖突解脫方式除了考慮起飛時(shí)間和高度層優(yōu)化分配,還可以考慮航跡的水平調(diào)整以及高度層的分航段優(yōu)化;算法在求解大規(guī)模航班時(shí)缺乏時(shí)效性,可能需要大量的計(jì)算時(shí)間。采用更加多樣化的沖突方式和將算法與局部搜索算法相結(jié)合以提高搜索效率預(yù)期將有助于提高算法的效率。

    猜你喜歡
    高度層父代航空器
    農(nóng)村家庭父代在家庭現(xiàn)代性轉(zhuǎn)型中的作用研究
    中國高等教育的代際傳遞及其內(nèi)在機(jī)制:“學(xué)二代”現(xiàn)象存在嗎?
    延遲退休決策對居民家庭代際收入流動(dòng)性的影響分析
    ——基于人力資本傳遞機(jī)制
    基于高度層的航路短時(shí)利用率模型研究
    男孩偏好激勵(lì)父代掙取更多收入了嗎?
    ——基于子女?dāng)?shù)量基本確定的情形
    論航空器融資租賃出租人的違約取回權(quán)
    航空器的順風(fēng)耳——機(jī)載衛(wèi)星通信
    太空探索(2016年3期)2016-07-12 09:58:45
    火星航空器何時(shí)才能首飛
    太空探索(2016年8期)2016-07-10 09:21:58
    MSG-3在小型航空器系統(tǒng)/動(dòng)力裝置維修要求制訂中的應(yīng)用
    飛機(jī)最佳航路爬升時(shí)機(jī)研究
    国产蜜桃级精品一区二区三区| 亚洲七黄色美女视频| 精品不卡国产一区二区三区| 少妇的逼水好多| 成人二区视频| 偷拍熟女少妇极品色| 丰满的人妻完整版| 国内少妇人妻偷人精品xxx网站| 亚洲精品一卡2卡三卡4卡5卡| 精品久久久久久久人妻蜜臀av| 97热精品久久久久久| av国产免费在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲一区高清亚洲精品| 午夜亚洲福利在线播放| 国产精品人妻久久久影院| 2021天堂中文幕一二区在线观| 国产视频一区二区在线看| 成人国产综合亚洲| 亚洲在线自拍视频| 校园春色视频在线观看| 校园春色视频在线观看| 九九爱精品视频在线观看| 赤兔流量卡办理| 国产伦在线观看视频一区| 国产伦一二天堂av在线观看| 一个人看的www免费观看视频| 99久久九九国产精品国产免费| 天堂av国产一区二区熟女人妻| 日韩欧美在线二视频| 久久久色成人| 免费搜索国产男女视频| 中文字幕熟女人妻在线| 国产伦精品一区二区三区视频9| 无遮挡黄片免费观看| 国内少妇人妻偷人精品xxx网站| .国产精品久久| 不卡视频在线观看欧美| 日日摸夜夜添夜夜添小说| 成人亚洲精品av一区二区| 亚洲欧美激情综合另类| 亚洲精华国产精华精| 深夜精品福利| 又紧又爽又黄一区二区| 亚洲国产日韩欧美精品在线观看| 内地一区二区视频在线| 美女高潮的动态| 好男人在线观看高清免费视频| 淫秽高清视频在线观看| 九九久久精品国产亚洲av麻豆| 深夜a级毛片| 亚洲成人精品中文字幕电影| 大型黄色视频在线免费观看| 国产精品综合久久久久久久免费| 久久香蕉精品热| 听说在线观看完整版免费高清| 三级国产精品欧美在线观看| 亚洲经典国产精华液单| 在线国产一区二区在线| 亚洲第一电影网av| 美女高潮的动态| 国产美女午夜福利| 亚洲国产色片| 最近最新免费中文字幕在线| 亚洲欧美日韩卡通动漫| a级毛片免费高清观看在线播放| 亚洲av成人av| 男女做爰动态图高潮gif福利片| 日韩中字成人| 精品人妻偷拍中文字幕| 网址你懂的国产日韩在线| 伊人久久精品亚洲午夜| 国产亚洲精品av在线| 国产又黄又爽又无遮挡在线| 日本三级黄在线观看| 欧美中文日本在线观看视频| 老女人水多毛片| 久久精品国产亚洲av天美| 干丝袜人妻中文字幕| 国产美女午夜福利| 狂野欧美激情性xxxx在线观看| 九色成人免费人妻av| 日韩欧美三级三区| 久久99热这里只有精品18| 亚州av有码| 天堂动漫精品| 天堂动漫精品| 男女边吃奶边做爰视频| 精品久久国产蜜桃| 国产精品人妻久久久影院| 久久精品夜夜夜夜夜久久蜜豆| 久久久久久伊人网av| 日韩欧美三级三区| 国产精品日韩av在线免费观看| 国产精品一区二区三区四区久久| 国产国拍精品亚洲av在线观看| 久久精品夜夜夜夜夜久久蜜豆| 热99在线观看视频| 色吧在线观看| 在线免费观看的www视频| 99久久精品热视频| 日日撸夜夜添| 亚洲av免费在线观看| 久久精品国产清高在天天线| 美女黄网站色视频| 琪琪午夜伦伦电影理论片6080| 国产欧美日韩精品亚洲av| 国产一区二区亚洲精品在线观看| 国产精品免费一区二区三区在线| 一级毛片久久久久久久久女| 久久九九热精品免费| 久久精品国产亚洲网站| 国产综合懂色| 91狼人影院| 黄色配什么色好看| 国产一级毛片七仙女欲春2| 亚洲第一区二区三区不卡| 亚洲成人久久性| 亚洲专区国产一区二区| 一个人看的www免费观看视频| .国产精品久久| 日本黄色视频三级网站网址| 一区二区三区激情视频| 精品一区二区三区人妻视频| 国产高潮美女av| 亚洲 国产 在线| 精品久久久久久,| 在线观看一区二区三区| 91午夜精品亚洲一区二区三区 | 999久久久精品免费观看国产| 亚洲电影在线观看av| 亚洲中文字幕一区二区三区有码在线看| 在线天堂最新版资源| 国产女主播在线喷水免费视频网站 | 国产精品一区二区三区四区免费观看 | 男女视频在线观看网站免费| 欧美zozozo另类| 色播亚洲综合网| videossex国产| 嫩草影院新地址| 亚洲va在线va天堂va国产| 成人无遮挡网站| 毛片一级片免费看久久久久 | 又爽又黄无遮挡网站| 国产日本99.免费观看| 国产视频内射| 桃色一区二区三区在线观看| 麻豆精品久久久久久蜜桃| 国产日本99.免费观看| 午夜爱爱视频在线播放| 禁无遮挡网站| 毛片女人毛片| 一本久久中文字幕| 在线观看午夜福利视频| 俺也久久电影网| 久久国内精品自在自线图片| av视频在线观看入口| 亚洲性夜色夜夜综合| 成年人黄色毛片网站| 久9热在线精品视频| 午夜激情欧美在线| 狂野欧美白嫩少妇大欣赏| 麻豆成人午夜福利视频| 三级毛片av免费| 狂野欧美白嫩少妇大欣赏| 午夜激情福利司机影院| 国产精品一区www在线观看 | 欧美激情久久久久久爽电影| 国产毛片a区久久久久| 婷婷亚洲欧美| 国产精品国产高清国产av| 成人高潮视频无遮挡免费网站| 亚洲最大成人中文| 婷婷精品国产亚洲av| 免费av不卡在线播放| 成年版毛片免费区| 国产精品1区2区在线观看.| 国产高清不卡午夜福利| 精品不卡国产一区二区三区| 夜夜看夜夜爽夜夜摸| 欧美区成人在线视频| av福利片在线观看| 国产精品98久久久久久宅男小说| 在线天堂最新版资源| 欧美日韩黄片免| 日本免费一区二区三区高清不卡| eeuss影院久久| 3wmmmm亚洲av在线观看| 搡女人真爽免费视频火全软件 | 国产伦在线观看视频一区| 欧美绝顶高潮抽搐喷水| 婷婷亚洲欧美| 国产高清不卡午夜福利| 亚洲国产高清在线一区二区三| 欧美日韩黄片免| 欧美+亚洲+日韩+国产| 日韩一本色道免费dvd| 亚洲美女黄片视频| 久久国产乱子免费精品| 国产精品久久久久久精品电影| 99热只有精品国产| 国产精品日韩av在线免费观看| 国产一区二区三区在线臀色熟女| 欧美日韩乱码在线| 男女视频在线观看网站免费| 性色avwww在线观看| 精品一区二区三区视频在线| 一级黄片播放器| 一级av片app| 亚洲久久久久久中文字幕| 最近最新中文字幕大全电影3| 久久久精品大字幕| 在线播放无遮挡| 欧美日韩瑟瑟在线播放| 国产熟女欧美一区二区| 亚洲成人久久性| av中文乱码字幕在线| 亚洲av.av天堂| 成人无遮挡网站| 欧美中文日本在线观看视频| av黄色大香蕉| 国产av麻豆久久久久久久| 人妻制服诱惑在线中文字幕| 欧美日韩黄片免| 精品久久久久久,| 国产高清视频在线观看网站| 午夜亚洲福利在线播放| 99riav亚洲国产免费| 女人被狂操c到高潮| 亚洲性久久影院| 欧美一区二区精品小视频在线| 久久这里只有精品中国| 女生性感内裤真人,穿戴方法视频| 亚洲va在线va天堂va国产| 久久精品国产亚洲av天美| 亚洲无线观看免费| 午夜福利高清视频| 久久99热这里只有精品18| 丰满乱子伦码专区| 最后的刺客免费高清国语| 久久精品国产亚洲网站| 久久久久精品国产欧美久久久| 亚洲中文日韩欧美视频| 一个人看的www免费观看视频| 久久久久久九九精品二区国产| 久久精品国产自在天天线| 国产亚洲91精品色在线| 欧美成人性av电影在线观看| 99精品久久久久人妻精品| 美女高潮的动态| 精品久久久久久久久久久久久| 中文字幕高清在线视频| 中文亚洲av片在线观看爽| 婷婷六月久久综合丁香| 内地一区二区视频在线| 精品一区二区三区av网在线观看| 欧美激情久久久久久爽电影| 22中文网久久字幕| 日韩欧美一区二区三区在线观看| 亚洲无线在线观看| 国产精华一区二区三区| 99久国产av精品| 大又大粗又爽又黄少妇毛片口| 黄色配什么色好看| 精品人妻1区二区| 黄色视频,在线免费观看| 男人的好看免费观看在线视频| 亚洲天堂国产精品一区在线| 日本-黄色视频高清免费观看| 香蕉av资源在线| 欧美性猛交黑人性爽| 俺也久久电影网| av黄色大香蕉| 伦精品一区二区三区| 悠悠久久av| 久久精品国产清高在天天线| 一级毛片久久久久久久久女| 久9热在线精品视频| 精品久久久久久久久久久久久| 日韩国内少妇激情av| 日韩欧美 国产精品| 成人国产麻豆网| 免费看美女性在线毛片视频| 嫁个100分男人电影在线观看| 亚洲精品456在线播放app | 国产高清激情床上av| av在线天堂中文字幕| 精品一区二区三区视频在线| 精品国内亚洲2022精品成人| 此物有八面人人有两片| 亚洲最大成人av| 国产视频内射| 国产精品一区二区三区四区免费观看 | 国产色爽女视频免费观看| av福利片在线观看| 久久99热这里只有精品18| 色综合站精品国产| 韩国av在线不卡| 亚洲最大成人中文| 久久精品国产自在天天线| 午夜老司机福利剧场| 国产欧美日韩精品一区二区| 亚洲 国产 在线| 深爱激情五月婷婷| 午夜福利成人在线免费观看| 国产探花极品一区二区| eeuss影院久久| 91久久精品国产一区二区三区| 成人无遮挡网站| 啦啦啦韩国在线观看视频| 99久久久亚洲精品蜜臀av| 久久人人爽人人爽人人片va| 国产伦人伦偷精品视频| 午夜福利在线观看吧| 精品久久久久久成人av| 中国美白少妇内射xxxbb| 色噜噜av男人的天堂激情| 好男人在线观看高清免费视频| 国产精品99久久久久久久久| 国产高清激情床上av| 日本欧美国产在线视频| 高清在线国产一区| 国产亚洲精品久久久久久毛片| 亚洲av成人av| 18禁黄网站禁片午夜丰满| 两性午夜刺激爽爽歪歪视频在线观看| 热99在线观看视频| 俺也久久电影网| 久久久久精品国产欧美久久久| 亚洲真实伦在线观看| 最近中文字幕高清免费大全6 | 亚洲在线自拍视频| 99国产极品粉嫩在线观看| 欧美黑人巨大hd| 真人做人爱边吃奶动态| 国产高清激情床上av| 能在线免费观看的黄片| 内射极品少妇av片p| 99热网站在线观看| 午夜激情欧美在线| 国产视频一区二区在线看| a级毛片免费高清观看在线播放| 一卡2卡三卡四卡精品乱码亚洲| 中文亚洲av片在线观看爽| 男女边吃奶边做爰视频| 麻豆av噜噜一区二区三区| 欧美不卡视频在线免费观看| 精品午夜福利在线看| 午夜精品久久久久久毛片777| 色5月婷婷丁香| 免费不卡的大黄色大毛片视频在线观看 | 欧美人与善性xxx| 草草在线视频免费看| 亚洲va在线va天堂va国产| 天堂av国产一区二区熟女人妻| 亚洲精品一区av在线观看| a在线观看视频网站| 天天一区二区日本电影三级| 欧美黑人欧美精品刺激| 非洲黑人性xxxx精品又粗又长| 少妇熟女aⅴ在线视频| 久久久久久久久大av| 99精品久久久久人妻精品| 婷婷亚洲欧美| 亚洲一级一片aⅴ在线观看| 国产亚洲欧美98| av在线亚洲专区| 午夜福利欧美成人| 午夜免费激情av| 国产精品野战在线观看| 美女免费视频网站| 久久久久九九精品影院| 麻豆一二三区av精品| 国产一区二区三区av在线 | 白带黄色成豆腐渣| 亚洲一级一片aⅴ在线观看| 久久久国产成人免费| 国产探花极品一区二区| 亚洲自偷自拍三级| 久久精品国产亚洲网站| 国语自产精品视频在线第100页| 亚洲熟妇中文字幕五十中出| 日韩欧美精品免费久久| 精品国内亚洲2022精品成人| 婷婷精品国产亚洲av在线| 亚洲aⅴ乱码一区二区在线播放| 99久久成人亚洲精品观看| 免费观看在线日韩| 国产一区二区三区在线臀色熟女| 女的被弄到高潮叫床怎么办 | 一个人观看的视频www高清免费观看| 国产激情偷乱视频一区二区| 婷婷色综合大香蕉| 免费搜索国产男女视频| 亚洲av中文字字幕乱码综合| 美女免费视频网站| 俺也久久电影网| 国产午夜精品论理片| 赤兔流量卡办理| 久久这里只有精品中国| 99久久精品国产国产毛片| 赤兔流量卡办理| 久久婷婷人人爽人人干人人爱| 日韩亚洲欧美综合| 中国美白少妇内射xxxbb| 美女免费视频网站| 国产人妻一区二区三区在| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久久久久大av| 国产中年淑女户外野战色| 又紧又爽又黄一区二区| 国产一区二区激情短视频| 亚洲欧美精品综合久久99| 欧美日韩综合久久久久久 | 亚洲成人精品中文字幕电影| 亚洲无线在线观看| 午夜福利视频1000在线观看| 欧美另类亚洲清纯唯美| 久久久久免费精品人妻一区二区| 人妻夜夜爽99麻豆av| 国产私拍福利视频在线观看| 亚洲七黄色美女视频| 欧美色欧美亚洲另类二区| 大又大粗又爽又黄少妇毛片口| 又紧又爽又黄一区二区| 3wmmmm亚洲av在线观看| 超碰av人人做人人爽久久| 熟女人妻精品中文字幕| 大又大粗又爽又黄少妇毛片口| 日本黄色视频三级网站网址| 色综合亚洲欧美另类图片| 岛国在线免费视频观看| 国产视频内射| 亚洲国产日韩欧美精品在线观看| 色哟哟哟哟哟哟| 88av欧美| 1024手机看黄色片| 小蜜桃在线观看免费完整版高清| 国产精品一及| 国产成年人精品一区二区| 女生性感内裤真人,穿戴方法视频| 蜜桃亚洲精品一区二区三区| 一区二区三区激情视频| 偷拍熟女少妇极品色| 国产色爽女视频免费观看| 深夜a级毛片| 热99re8久久精品国产| 不卡一级毛片| 国产亚洲av嫩草精品影院| 欧美三级亚洲精品| 超碰av人人做人人爽久久| 亚洲av中文字字幕乱码综合| 男女下面进入的视频免费午夜| 真人一进一出gif抽搐免费| 免费看光身美女| 少妇的逼好多水| 99久国产av精品| 免费在线观看影片大全网站| 午夜亚洲福利在线播放| 日日夜夜操网爽| 午夜a级毛片| 热99re8久久精品国产| 国产精品1区2区在线观看.| 丰满乱子伦码专区| 成人美女网站在线观看视频| 国产精品无大码| 久久久色成人| 亚洲av五月六月丁香网| 国国产精品蜜臀av免费| 亚洲美女搞黄在线观看 | 夜夜看夜夜爽夜夜摸| 天美传媒精品一区二区| 人妻丰满熟妇av一区二区三区| 日韩欧美精品v在线| 亚洲人成伊人成综合网2020| 国产爱豆传媒在线观看| 最近在线观看免费完整版| 国产亚洲精品久久久久久毛片| 搡老妇女老女人老熟妇| 免费观看精品视频网站| 久久亚洲精品不卡| 欧美色欧美亚洲另类二区| 亚洲电影在线观看av| 91在线精品国自产拍蜜月| a级一级毛片免费在线观看| 色综合亚洲欧美另类图片| 久久久国产成人精品二区| 亚洲成人久久性| 欧美黑人巨大hd| 91狼人影院| 国产成人福利小说| 又爽又黄a免费视频| 精品人妻视频免费看| 简卡轻食公司| 久久香蕉精品热| 欧美成人性av电影在线观看| 看片在线看免费视频| 国产91精品成人一区二区三区| 最新在线观看一区二区三区| 永久网站在线| 久久亚洲真实| 伦精品一区二区三区| 亚洲熟妇熟女久久| 国产av一区在线观看免费| 日本与韩国留学比较| 在现免费观看毛片| av中文乱码字幕在线| 91在线精品国自产拍蜜月| 赤兔流量卡办理| 国产高清视频在线观看网站| 欧美高清性xxxxhd video| 午夜影院日韩av| 乱系列少妇在线播放| 日韩中文字幕欧美一区二区| 亚洲aⅴ乱码一区二区在线播放| 日韩在线高清观看一区二区三区 | 女人十人毛片免费观看3o分钟| 国产亚洲精品av在线| 欧美黑人欧美精品刺激| 日本在线视频免费播放| 国产精品福利在线免费观看| 亚洲av免费高清在线观看| 波野结衣二区三区在线| 日日干狠狠操夜夜爽| 久久久久免费精品人妻一区二区| АⅤ资源中文在线天堂| 精品人妻一区二区三区麻豆 | 国产探花极品一区二区| 日韩强制内射视频| 成人欧美大片| 五月玫瑰六月丁香| 久久精品国产鲁丝片午夜精品 | 精品无人区乱码1区二区| 成人av一区二区三区在线看| 国模一区二区三区四区视频| 天天一区二区日本电影三级| 黄色视频,在线免费观看| 99久久精品一区二区三区| 五月玫瑰六月丁香| 尤物成人国产欧美一区二区三区| eeuss影院久久| 女人被狂操c到高潮| 国产av在哪里看| 狂野欧美白嫩少妇大欣赏| 两个人的视频大全免费| 男插女下体视频免费在线播放| 国产精品一区二区三区四区免费观看 | 九九在线视频观看精品| 国产精品久久久久久av不卡| 又黄又爽又免费观看的视频| 岛国在线免费视频观看| 熟妇人妻久久中文字幕3abv| 偷拍熟女少妇极品色| 中国美白少妇内射xxxbb| 亚洲av一区综合| 中国美白少妇内射xxxbb| 噜噜噜噜噜久久久久久91| 国产欧美日韩精品一区二区| 十八禁网站免费在线| 人妻丰满熟妇av一区二区三区| 十八禁网站免费在线| 午夜视频国产福利| 国产精品三级大全| 日韩在线高清观看一区二区三区 | 日本撒尿小便嘘嘘汇集6| av中文乱码字幕在线| 一本一本综合久久| 日韩精品中文字幕看吧| 国产黄色小视频在线观看| 国产午夜精品论理片| 搡老妇女老女人老熟妇| 少妇高潮的动态图| 成人无遮挡网站| 色哟哟哟哟哟哟| 亚洲成a人片在线一区二区| 搡老岳熟女国产| 亚洲久久久久久中文字幕| 深夜精品福利| 亚洲av中文字字幕乱码综合| 日韩强制内射视频| 欧美日本视频| 亚洲性夜色夜夜综合| 欧美日韩国产亚洲二区| 欧美成人性av电影在线观看| 国产精品一区二区三区四区免费观看 | 欧美又色又爽又黄视频| 精品无人区乱码1区二区| 亚洲精品在线观看二区| 97超级碰碰碰精品色视频在线观看| 久久人人精品亚洲av| 欧美色欧美亚洲另类二区| 成人三级黄色视频| 校园春色视频在线观看| 国产精品免费一区二区三区在线| 国产色婷婷99| 最近视频中文字幕2019在线8| 国产伦在线观看视频一区| 日本a在线网址| 少妇高潮的动态图| 最近最新免费中文字幕在线| 国产老妇女一区| 久久国内精品自在自线图片| 99久久久亚洲精品蜜臀av| 可以在线观看的亚洲视频| 国产三级中文精品| 一级av片app| 日韩精品有码人妻一区| 尤物成人国产欧美一区二区三区| 91麻豆精品激情在线观看国产| 国产成人a区在线观看| 免费高清视频大片| 欧美潮喷喷水| 欧美性感艳星| 少妇高潮的动态图| 男插女下体视频免费在线播放|