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

        ?

        基于遺傳算法的平面交叉口信號配時優(yōu)化

        2020-12-14 04:26:52陸穎林麗
        森林工程 2020年6期
        關(guān)鍵詞:遺傳算法

        陸穎 林麗

        摘 要:為了提高平面交叉口的運行效率,以平面交叉口為研究對象,建立交叉口信號配時參數(shù)優(yōu)化的非線性模型,該模型以交叉口的平均延誤值最小為目標,約束條件為有效綠燈時間和周期時長,并采用浮點數(shù)編碼遺傳算法在Matlab軟件中對信號配時參數(shù)優(yōu)化模型進行求解。為驗證優(yōu)化模型的正確性和遺傳算法的有效性,選取江蘇省鹽城市某一交叉口為算例,采用遺傳算法對該模型進行求解。結(jié)果顯示,優(yōu)化后的平面交叉口的周期時長從123 s下降到97 s,下降率為21.1%;平均延誤從32.79 s下降到14.91 s,下降率為54.5%;平面交叉口的服務(wù)水平從C等級上升到B等級。最后通過Vissim仿真對算例進行進一步驗證,結(jié)果表明平面交叉口平均延誤值下降幅度為24.1%。故遺傳算法適用于平面交叉口信號配時方案優(yōu)化,并且具有一定的可靠性與優(yōu)越性。

        關(guān)鍵詞:平面交叉口;信號配時;平均延誤;非線性模型;遺傳算法

        中圖分類號:U491.51??? 文獻標識碼:A?? 文章編號:1006-8023(2020)06-0103-07

        Signal Timing Optimization of Intersections Based on Genetic Algorithm

        LU Ying,LIN Li*

        (College of Automobile and Traffic Engineering, Nanjing Forestry University, Nanjing 210037, China)

        Abstract:In order to improve the operation efficiency of a plane intersection, a nonlinear model for the optimization of the signal timing parameters of a plane intersection was established with the intersection as the research object. The goal of the model was to minimize the average delay value of the intersection, and the constraints are the effective green light time and cycle time. The model was solved by using the floating-point coding genetic algorithm in Matrix Laboratory software. In order to verify the correctness of the optimization model and the effectiveness of genetic algorithm, an intersection in Yancheng city of Jiangsu Province was selected as an example, and the model was solved by genetic algorithm. The results showed that the cycle time of the intersection after optimization decreased from 123 s to 97 s, with a decrease rate of 21.1%; the average delay decreased from 32.79 s to 14.91 s, the decline rate was 54.5%; the service level of intersections increased from C level to B level. At last, the example was further verified by VISSIM simulation. The result showed that the average delay value of the intersection decreased by 24.1%. Therefore, the genetic algorithm was suitable for the optimization of the signal timing scheme of the intersection, and had certain reliability and superiority.

        Keywords:Plane intersection; signal timing; average delay; nonlinear model; genetic algorithm

        收稿日期:2020-05-10

        基金項目:住房城鄉(xiāng)建設(shè)部2014年科學(xué)技術(shù)計劃項目(2014-k5-007)

        第一作者簡介:陸穎,碩士研究生。研究方向:交通運輸規(guī)劃與管理,E-mail:1124065472@qq.com

        通信作者:林麗,碩士,副教授。研究方向:交通運輸規(guī)劃與管理,E-mail:LinLi401@njfu.com.cn

        引文格式:陸穎,林麗.基于遺傳算法的平面交叉口信號配時優(yōu)化[J].森林工程,2020,36(6):103-109.

        LU Y, LIN L. Signal timing optimization of intersections based on genetic algorithm[J]. Forest Engineering,2020,36(6):103-109.

        0 引言

        國家經(jīng)濟的迅猛發(fā)展導(dǎo)致城市用地規(guī)模不斷擴張,城市道路是城市的筋脈,而平面交叉口是城市道路交通的重要組成部分,同時也是城市道路交通的“堵點”所在[1]。因此,對平面交叉口的研究與治理將一直持續(xù),其中對平面交叉口交通控制方面進行改善不失為一種有效方法[2]。其中對信號配時進行優(yōu)化,不需要對平面交叉口的空間資源進行大規(guī)模改擴建,是一種成本低、見效快的有效方式。但是,目前許多城市道路信號配時方案還是根據(jù)經(jīng)典Webster法進行求解,計算所求得的有效綠燈時間運用到交通流中產(chǎn)生的延誤值往往較大,并會出現(xiàn)進口道的時間損失,浪費平面交叉口的時間和空間資源。

        相比較經(jīng)典的Webster法,智能優(yōu)化算法具有一定的先進性,國內(nèi)外大量學(xué)者將各種優(yōu)化算法應(yīng)用到智能交通控制領(lǐng)域中,在理論研究方面取得了較多的成果。Rouphail等[3]提出基于信號配時優(yōu)化的CORSIM模型,通過遺傳算法對模型進行求解,以網(wǎng)絡(luò)延誤和排隊時間為評價指標。Kim等[4]提出利用遺傳算法自動確定模糊控制交通系統(tǒng)中的最優(yōu)參數(shù)和隸屬函數(shù),并且通過多個平面交叉口仿真驗證了該方法的有效性。高云峰等[5]根據(jù)延誤、停車率和排隊長度3個優(yōu)化指標,建立了基于非支配排序遺傳算法的交叉口信號配時多目標模型,但是算例結(jié)果顯示平均停車率比Webster算法得出的數(shù)值要高,3個優(yōu)化指標之間并不能達到平衡。伍尚昆等[6]將改進的蟻群算法運用到多目標信號配時優(yōu)化模型中,并且通過交叉口平峰和高峰兩個時段進行模型驗證,但是在平峰階段,經(jīng)典的Webster算法得出的總延誤時間最小。李振龍等[7]從車輛延誤、排隊長度和尾氣排放量3個性能指標,建立了以平面交叉口的信號控制優(yōu)化模型,并采用遺傳算法進行求解。劉暢等[8]以各相位有效綠燈時間為自變量建立多目標信號配時優(yōu)化模型,其中優(yōu)化目標為平面交叉口人均延誤和人均 CO 排放,并通過粒子群算法進行求解。姚志洪等[9]以累計排隊車輛延誤為基礎(chǔ),提出平面交叉口信號配時雙目標模型即通行能力最大、平均延誤最小,但是算例結(jié)果表明兩者并不能同時達到預(yù)期效果,模型效益不能最大化。

        總體來看,粒子群算法、蟻群算法和遺傳算法在信號配時參數(shù)計算中運用較為廣泛,但是對于研究非線性規(guī)劃模型求解問題,粒子群算法[10]在收斂性分析方面較薄弱,而遺傳算法不僅擁有成熟的收斂性分析方法,并且能夠估計收斂速度,蟻群算法[11]適合求解圖中搜索路徑問題,計算開銷大,而遺傳算法在選擇操作時將會淘汰一些較差的種群,從而減少之后的計算開銷。相比較于傳統(tǒng)的遺傳算法在編碼階段通常使用二進制編碼,而二進制編碼并不適合求解多維、高精度的連續(xù)函數(shù)優(yōu)化問題,而浮點數(shù)編碼不僅適合求解高精度函數(shù)問題,并且能夠提高運算效率。因此,本文以平面交叉口平均延誤最小為目標,建立非線性規(guī)劃模型,采用浮點數(shù)編碼遺傳算法進行求解。該算法可通過常規(guī)軟件進行操作,可廣泛應(yīng)用于交叉口信號配時優(yōu)化工作,實用性強。

        1 建立信號配時優(yōu)化模型

        平均延誤、通行能力、飽和度、最大排隊長度、平均停車次數(shù)、停車率以及油耗等通常作為交通效益評價的指標[12]。其中平均延誤的大小可以最直觀的體現(xiàn)平面交叉口的運行效率,故本文選取平均延誤最小值作為目標函數(shù)。

        1.1平均延誤

        平面交叉口的延誤時間是指車輛在擁堵情況下,通過交叉口所需時間與暢通行駛同樣距離所需時間的差值[13]。本文依據(jù)HCM延誤公式[12],將延誤分為均勻延誤和隨機附加延誤進行計算,其計算公式為:

        2 遺傳算法

        遺傳算法(Genetic Algorithms,GA)是一種全局自適應(yīng)概率搜索算法,以自然選擇和基因遺傳學(xué)原理作為基礎(chǔ),以生物進化優(yōu)勝劣汰的自然選擇機理和生物界繁衍進化的遺傳機制作為指導(dǎo)思想[14]。

        本文建立的信號配時優(yōu)化非線性模型具有不規(guī)則、不連續(xù)的特點,遺傳算法定義的適應(yīng)度函數(shù)有很強的適應(yīng)性,即便是非規(guī)則、非連續(xù)或是有噪音的情況下,也能找到全局最優(yōu)解,因此遺傳算法適合求解本文的信號配時優(yōu)化模型。遺傳算法的具體步驟如下。

        (1)初始參數(shù)的設(shè)定。種群規(guī)模n、交叉概率Pc、變異概率Pm、進化代數(shù)t。

        (2)染色體編碼。本文選取浮點數(shù)編碼方法,該方法可以克服二進制編碼因優(yōu)化問題的精度要求而導(dǎo)致編碼串很長、搜索空間急劇擴大、計算量大且時間長的一系列問題。其中,個體的編碼長度與其決策變量的位數(shù)相等。

        (3)適應(yīng)度函數(shù)。本文以各相位有效綠燈時間ge作為優(yōu)化變量,以平均延誤為優(yōu)化目標函數(shù)dI,有必要將目標函數(shù)轉(zhuǎn)換成最大化問題形式,且為函數(shù)值非負的適應(yīng)度函數(shù)。本文要求解平面交叉口平均延誤的最小值,為了滿足適應(yīng)度函數(shù)的條件,適應(yīng)度函數(shù)DI和目標函數(shù)dI之間的對應(yīng)關(guān)系需要被建立,本文以目標函數(shù)的倒數(shù)作為適應(yīng)度函數(shù),即:

        DI=1dI。? (13)

        (4)選擇。選取輪盤賭法作為選擇算子,該方法是一種有退還的隨機選擇策略。群體中的每個個體都是圓盤中扇形的一部分,扇形面積越大,個體的適應(yīng)值就越高,該個體進入下一輪選擇的概率也就越大。假設(shè)種群個體數(shù)目為N,個體i的適應(yīng)度為Fi,則它被選中的概率為:

        Pi=Fi∑Nj=1Fj。? (14)

        (5)交叉。選擇實數(shù)交叉法[15],該方法是指兩個基因個體通過線性組合產(chǎn)生兩個新的個體,如第m個染色體am和第n個染色體an在j位的交叉操作方法為:

        amj=amj(1-b)+anjb

        anj=anj(1-b)+amjb。(15)

        式中:b是區(qū)間[0,1]之間的隨機數(shù)。

        (6)變異。采用非均勻變異。第i個個體的第j個基因aij進行變異的操作方法為:

        aij=aij+(aij-amax)f(t),r>0.5

        aij+(amin-aij)f(t),r≤0.5。(16)

        式中:amax是基因aij的上界;amin是基因aij的下界;f(t)=r2(1-t/T)2,r2是隨機數(shù);t為當前迭代次數(shù);T為最大進化次數(shù),r為[0,1]區(qū)間的隨機數(shù)。

        (7)進化終止條件。設(shè)最大進化代數(shù)為T,當?shù)螖?shù)滿足t=T條件時,基于遺傳算法的最優(yōu)解為種群中適應(yīng)度值最高的個體,最后進行個體解碼工作,得到平面交叉口的信號配時優(yōu)化方案[16]。

        3 實例分析

        為了驗證本文信號配時優(yōu)化模型的正確性以及遺傳算法求解的優(yōu)越性,選取江蘇省鹽城市某一平面交叉口進行實地調(diào)查。該平面交叉口的幾何構(gòu)造和渠化方案如圖1所示。該平面交叉口的交通控制方案為三相位,第一相位為東西進口道直行與右轉(zhuǎn)車輛放行;第二相位為東西進口道左轉(zhuǎn)車輛放行;第三相位為南北進口道左轉(zhuǎn)、直行和右轉(zhuǎn)車輛一起放行。經(jīng)過多次調(diào)查發(fā)現(xiàn),該平面交叉口的高峰小時出現(xiàn)在早上8:00—9:00,選取一天中的高峰小時對其進行高峰小時交通調(diào)查,并進行當量交通量的換算,最后調(diào)查結(jié)果見表1。

        首先將現(xiàn)狀交通量數(shù)據(jù)和配時參數(shù)運用Webster法重新進行配時計算,其次將本文建立的信號配時優(yōu)化模型在Matlab環(huán)境下進行遺傳算法求解。根據(jù)高峰小時周期時長推薦值[10]設(shè)置最小信號周期為40 s,最大信號周期為120 s;為了滿足行人過街時長,設(shè)置第一相位最小有效綠燈時間為15 s,最大有效綠燈時間為60 s;第二相位最小有效綠燈時間為10 s,最大綠燈時間為50 s;第三相位最小有效綠燈時間為30 s,最大綠燈時間為60 s。

        遺傳算法中控制參數(shù)的選擇也非常關(guān)鍵,通常這些參數(shù)包括種群規(guī)模n、交叉概率Pc、變異概率Pm和進化代數(shù)t。種群規(guī)模的大小會直接影響遺傳算法的收斂性,推薦值為10~200,實例中取值為n=50;交叉概率控制著交叉操作被使用的頻度,一般建議Pc的取值范圍為0.4~0.99,實例中取值為Pc=0.6;變異概率控制著變異操作被使用的概率,一般建議值為0.000 1~0.1,實例中取值為Pm=0.01;進化代數(shù)t的一般取值為100~400,實例中取值為t=300。

        確定各個控制參數(shù)后,將目標函數(shù)文件,選擇、交叉及變異文件與算法主函數(shù)文件在同一文件夾下運行,在300次運行優(yōu)化過程中,目標函數(shù)從120代開始收斂,即目標函數(shù)達到最優(yōu),運行結(jié)果見表2。從表2可知,周期時長從現(xiàn)狀123 s下降到97 s,下降幅度為21.1%;平面交叉口的平均延誤從現(xiàn)狀32.79 s/pcu下降到14.91 s/pcu,下降率達54.5%;由遺傳算法求解得出的信號配時如圖2所示,第一相位的綠燈時間為46 s,第二相位的綠燈時間為11 s,第三相位綠燈時間為31 s。平面交叉口機動車延誤與服務(wù)水平對應(yīng)關(guān)系[17]見表3。由表3標準可知,現(xiàn)狀交叉口平均延誤為32.79 s/pcu,對應(yīng)交叉口服務(wù)水平為C;Webster算法計算得出的平均延誤對照表3,其交叉口服務(wù)水平依然為C;遺傳算法得出的結(jié)果服務(wù)水平上升到B。Webster算法和遺傳算法在一定程度上都使得現(xiàn)狀交叉口的平均延誤值得到下降,但從下降幅度來看遺傳算法對信號配時優(yōu)化模型的求解要優(yōu)于現(xiàn)狀以及經(jīng)典的Webster算法,使得平面交叉口的平均延誤進一步降低,服務(wù)水平從C等級上升到B等級,平面交叉口的運行效率隨之提升。

        為了進一步驗證遺傳算法的優(yōu)越性,將現(xiàn)狀、Webster算法以及遺傳算法得出的配時方案在Vissim4.3中運行。根據(jù)表1平面交叉口當日高峰小時數(shù)據(jù),按換算后的當量小汽車數(shù)據(jù)進行車輛輸入,建軍東路(即東西向)小汽車期望速度為40~60 km/h,富康路(即南北向)小汽車期望速度為20~40 km/h,分別對應(yīng)相應(yīng)的期望車速分布曲線;對平面交叉口各個進口道分別設(shè)置檢測器對三項方案分別進行延誤指標評價,仿真時長為3 600 s。仿真結(jié)果見表4。

        根據(jù)表4中的仿真延誤數(shù)據(jù),總體來看大小順序為:遺傳算法得出的交叉口平均延誤值、Webster算法延誤值、現(xiàn)狀延誤值。從進口道延誤來看,南、北進口增減量變化幅度最大,與現(xiàn)狀平均延誤值相比分別下降了17.5? s/pcu和18.5 s/pcu;東、西進口延誤增減量幅度較小。從交叉口總體平均延誤來看,基于遺傳算法的平均延誤值相比于現(xiàn)狀值下降了5.6 s/pcu,下降幅度近似為24.1%;相比于Webster算法,下降了1.6 s/pcu,下降率幅度近似為8.3%。因此,基于遺傳算法信號配時優(yōu)化模型求解得出的平面交叉口平均延誤值要低于Webster算法,可使平面交叉口平均延誤值進一步降低,提高交叉口的運行效率。

        4 結(jié)束語

        信號配時參數(shù)的合理設(shè)置可以影響整個平面交叉口車流的暢通度,本文將平面交叉口的交通效益做出了量化評價,以平面交叉口平均延誤值最小為目標函數(shù),建立了以有效綠燈時間和周期時長作為約束條件的信號配時優(yōu)化模型,通過浮點數(shù)編碼遺傳算法對模型進行求解,得出的結(jié)論如下。

        (1)相比于現(xiàn)狀和Webster算法,遺傳算法求解得出的配時方案運用于平面交叉口,可使得交叉口的平均延誤進一步降低,并且使得平面交叉口的服務(wù)水平上升到B等級,同時通過Vissim4.3仿真實驗進一步驗證了遺傳算法對平面交叉口信號配時算法的適用性。

        (2)本文的信號配時優(yōu)化模型廣泛適用于城市道路的平面交叉口,并且遺傳算法的求解過程也可通過常規(guī)軟件進行操作,實用性強,效果明顯。

        對于本文建立的信號配時優(yōu)化模型,后續(xù)研究中可加入最大排隊長度、通行能力等指標,建立多目標信號配時優(yōu)化模型,從而更全面地評價平面交叉口的交通效益。

        【參 考 文 獻】

        [1]成衛(wèi),袁滿榮,姬利娜.基于粒子群算法的道路交叉口信號配時優(yōu)化模型:以昆明市為例[J].城市交通,2013,11(3):27-33.

        CHENG W, YUAN M R, JI L N. Signal timing model based on particle swarm optimization: a case study of Kunming[J]. Urban Transport of China, 2013, 11(3): 27-33.

        [2]田豐,邊婷婷.基于自適應(yīng)遺傳算法的交通信號配時優(yōu)化[J].計算機仿真,2010,27(6):305-308.

        TIAN F, BIAN T T. Traffic signal timing optimization based on adaptive genetic algorithm[J]. Computer Simulation, 2010, 27(6): 305-308.

        [3]ROUPHAIL N M, PARK B B, SACKS J. Direct signal timing optimization: Strategy development and results[C]// In XI Pan American Conference in Traffic and Transportation Engineering, 2000.

        [4]KIM J, KIM B M, HUH N C. Genetic algorithm approach to generate rules and membership functions of fuzzy traffic controller[C]// 10th IEEE International Conference on Fuzzy Systems.(Cat. No. 01CH37297), IEEE, Melbourne, Australia, 2001, 1: 525-528.

        [5]高云峰,徐立鴻,胡華,等.交叉口定周期信號控制多目標優(yōu)化方法[J].中國公路學(xué)報,2011,24(5):82-88.

        GAO Y F, XU L H, HU H, et al. Multi-objective optimization method for fixed-time signal control at intersection[J]. China Journal of Highway and Transport, 2011, 24(5): 82-88.

        [6]伍尚昆,陳翠宜,祝勝林.基于多種群蟻群算法的交叉路口信號配時優(yōu)化[J].計算機應(yīng)用與軟件,2014,31(5):83-88.

        WU S K, CHEN C Y, ZHU S L. Timing optimisation for intersection signal based on multi-colony ant algorithm[J]. Computer Applications and Software, 2014, 31(5): 83-88.

        [7]李振龍,董文會,韓建龍,等.基于遺傳算法的交叉口信號控制多目標優(yōu)化[J].計算機應(yīng)用,2016,36(S2):82-84.

        LI Z L, DONG W H, HAN J L, et al. Multi-object optimization for intersection signal control based on genetic algorithm[J]. Journal of Computer Applications, 2016, 36(S2): 82-84.

        [8]劉暢,魏麗英.考慮人均延誤和人均排放的信號配時優(yōu)化模型[J].哈爾濱工業(yè)大學(xué)學(xué)報,2018,50(9):83-88.

        LIU C, WEI L Y. Signal timing optimization model considering per capita delay and per capita emissions[J]. Journal of Harbin Institute of Technology, 2018, 50(9): 83-88.

        [9]姚志洪,蔣陽升,吳北川.基于休假排隊理論的信號配時仿真優(yōu)化模型[J].計算機仿真,2018,35(6):156-159.

        YAO Z H, JIANG Y S, WU B C. Traffic signal timing simulation optimization based on vacation queueing theory[J]. Computer Simulation, 2018, 35(6): 156-159.

        [10]首艷芳,徐建閩.信號交叉口多目標動態(tài)決策模型及其優(yōu)化方法[J].公路交通科技,2012,29(11):92-97.

        SHOU Y F, XU J M. Multi-objective dynamic decision-making model of signalized intersection and its optimization method[J]. Journal of Highway and Transportation Research and Development, 2012, 29(11): 92-97.

        [11]肖業(yè)偉,黃輝先,王宸昊.基于蟻群算法的交叉路口多相位信號配時優(yōu)化[J].計算機工程與應(yīng)用,2008,44(19):241-244.

        XIAO Y W, HUANG H X, WANG C H. Multiphase traffic signal timing optimization of intersection based on ant colony algorithm[J]. Computer Engineering and Applications, 2008, 44(19): 241-244.

        [12]吳兵.交通管理與控制(5版)[M].北京:人民交通出版社,2015.

        WU B. Traffic management and control(Fifth edition)[M]. Beijing: Peoples Traffic Press, 2015.

        [13]慕飛飛,張惠珍.基于遺傳算法的單點交叉口信號配時優(yōu)化[J].上海理工大學(xué)學(xué)報,2015,37(6):600-604.

        MU F F, ZHANG H Z. Signal timing optimization at single-point intersection based on genetic algorithm[J]. Journal of University of Shanghai for Science and Technology, 2015, 37(6): 600-604.

        [14]卓金武, 王鴻鈞. MATLAB數(shù)學(xué)建模方法與實踐[M]. 北京: 北京航空航天大學(xué)出版社, 2018.

        ZHUO J W, WANG H J. Mathematical modeling method and practice of MATLAB [M]. Beijing: Beijing University of Aeronautics & Astronautics Press, 2018.

        [15]李碩,付珊,賀文,等.基于遺傳算法的左轉(zhuǎn)待行區(qū)交叉口信號配時優(yōu)化研究[J].公路工程,2018,43(5):131-137.

        LI S, FU S, HE W, et al. Study on signal timing optimization of the intersection with left-turn waiting area based on genetic algorithm[J]. Highway Engineering, 2018, 43(5): 131-137.

        [16]王逸,姚志洪,蔣陽升,等.基于自適應(yīng)遺傳算法的雙環(huán)信號配時優(yōu)化模型[J].工業(yè)工程,2018,21(5):72-80.

        WANG Y, YAO Z H, JIANG Y S, et al. The dual-phase signal timing optimization model based on adaptive genetic algorithm[J]. Industrial Engineering Journal, 2018, 21(5): 72-80.

        [17]王秋平,譚學(xué)龍,張生瑞.城市單點交叉口信號配時優(yōu)化[J].交通運輸工程學(xué)報,2006,6(2):60-64.

        WANG Q P, TAN X L, ZHANG S R. Signal timing optimization of urban single-point intersections[J]. Journal of Traffic and Transportation Engineering, 2006, 6(2): 60-64.

        猜你喜歡
        遺傳算法
        基于遺傳算法的模糊控制在過熱汽溫控制系統(tǒng)優(yōu)化中的應(yīng)用
        電子制作(2019年16期)2019-09-27 09:34:44
        遺傳算法對CMAC與PID并行勵磁控制的優(yōu)化
        基于自適應(yīng)遺傳算法的CSAMT一維反演
        基于遺傳算法的建筑物沉降回歸分析
        一種基于遺傳算法的聚類分析方法在DNA序列比較中的應(yīng)用
        基于遺傳算法和LS-SVM的財務(wù)危機預(yù)測
        遺傳算法識別模型在水污染源辨識中的應(yīng)用
        協(xié)同進化在遺傳算法中的應(yīng)用研究
        軟件發(fā)布規(guī)劃的遺傳算法實現(xiàn)與解釋
        基于改進的遺傳算法的模糊聚類算法
        女同性恋一区二区三区四区| 高h小月被几个老头调教| 狠狠躁天天躁无码中文字幕图| 久激情内射婷内射蜜桃| 人与嘼交av免费| 日本加勒比东京热日韩| 国产美女黄性色av网站| 91九色播放在线观看| 97在线视频免费人妻| 成人a级视频在线观看| 亚洲av区无码字幕中文色| 国产丝袜一区二区三区在线不卡| 丰满人妻AV无码一区二区三区| 色av色婷婷18人妻久久久| 久久久精品视频网站在线观看| 亚洲av无码专区首页| 小12箩利洗澡无码视频网站| 中文字幕一区二区在线| 国产精品熟女一区二区三区| 亚洲va无码va在线va天堂| 国产精品无套内射迪丽热巴| 精品的一区二区三区| 久久精品国产熟女亚洲av麻豆| 天天摸夜夜摸摸到高潮| 黑人巨大白妞出浆| 官网A级毛片| 亚洲不卡在线免费视频| 成 人片 黄 色 大 片| 久久综合精品国产丝袜长腿| 亚洲国产精品久久久久婷婷软件| 国产特黄a三级三级三中国| 久久中文字幕人妻淑女| 欧洲精品免费一区二区三区| 亚洲日韩精品久久久久久| av国产免费在线播放| 久久精品国产亚洲av香蕉| 亚洲影院天堂中文av色| 久久久亚洲精品一区二区| 国产日本精品视频一区二区| 无码成人aaaaa毛片| 国产精品久久1024|