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

    基于改進(jìn)魚群算法的城軌列車運(yùn)行均衡調(diào)整研究

    2020-12-16 10:10:58楊艷娟徐永能董瑞超景順利姜建萍
    兵器裝備工程學(xué)報(bào) 2020年11期
    關(guān)鍵詞:晚點(diǎn)魚群列車運(yùn)行

    楊艷娟,徐永能,董瑞超,景順利,姜建萍

    (1.南京理工大學(xué) a.自動(dòng)化學(xué)院; b.瞬態(tài)物理國(guó)家重點(diǎn)實(shí)驗(yàn)室, 南京 210094;2.南京恩瑞特實(shí)業(yè)有限公司, 南京 211106)

    隨著基于無(wú)線通信的CBTC(Communication Based Train Control System)系統(tǒng)在城市軌道交通列車上的應(yīng)用,促進(jìn)了列車控制技術(shù)的發(fā)展,在很大程度上縮短了列車的追蹤間隔。但是,列車在實(shí)際運(yùn)行過(guò)程中容易受到隨機(jī)干擾因素的影響,使得列車偏離計(jì)劃運(yùn)行圖。由于列車運(yùn)行間隔短,若前行列車晚點(diǎn)時(shí)間過(guò)長(zhǎng)時(shí),可能導(dǎo)致列車運(yùn)行的大面積紊亂,導(dǎo)致列車群連帶晚點(diǎn)。因此,列車運(yùn)行自動(dòng)控制系統(tǒng)需要采取一定的調(diào)整措施,對(duì)偏離計(jì)劃運(yùn)行圖的列車進(jìn)行實(shí)時(shí)、高效地調(diào)整,使列車盡快恢復(fù)按圖定計(jì)劃運(yùn)行,達(dá)到保障城市軌道交通運(yùn)營(yíng)效率和服務(wù)水平的目的。

    目前,針對(duì)城市軌道交通列車運(yùn)行調(diào)整已經(jīng)取得了豐富的研究成果,國(guó)內(nèi)外學(xué)者大多從構(gòu)建模型的角度出發(fā),將列車調(diào)整問(wèn)題抽象為具體的數(shù)學(xué)模型[1-4],并采用粒子群算法[5-7]、遺傳算法[8-11]、螢火蟲算法[12]等智能優(yōu)化算法對(duì)模型進(jìn)行求解。本研究從綜合考慮乘客出行體驗(yàn)和企業(yè)運(yùn)營(yíng)服務(wù)考核指標(biāo)的角度出發(fā),以列車的總晚點(diǎn)時(shí)間和乘客延誤等待總時(shí)間最少為優(yōu)化目標(biāo),建立城市軌道交通列車運(yùn)行均衡調(diào)整模型,并提出一種改進(jìn)人工魚群算法對(duì)模型進(jìn)行求解。

    1 列車運(yùn)行均衡調(diào)整模型構(gòu)建

    在構(gòu)建城市軌道交通列車運(yùn)行均衡調(diào)整模型時(shí),本研究給出以下假設(shè):

    1) 列車晚點(diǎn)時(shí)間在ATR(Automatic Train Regulation,ATR)閾值內(nèi),無(wú)需人工介入調(diào)整;

    2) 列車晚點(diǎn)不包括改變列車運(yùn)行次序、增加/取消列車車次以及發(fā)生設(shè)備故障的情況;

    3) 延誤只對(duì)發(fā)生初始延誤的后續(xù)列車有影響。

    1.1 約束條件

    城市軌道交通列車調(diào)整的約束條件是復(fù)雜且相對(duì)固定的,主要包括停站時(shí)間約束、區(qū)間運(yùn)行時(shí)間約束、追蹤間隔時(shí)間約束和到發(fā)車時(shí)間約束等四項(xiàng)約束[13-14]。

    1) 停站時(shí)間約束

    停站約束時(shí)間可以表示成:

    (1)

    其中,Ts min和Ts max分別表示列車li在車站sj的最小停站時(shí)間和最大停站時(shí)間。

    2) 區(qū)間運(yùn)行時(shí)間約束

    為了保障乘客出行的安全性和舒適性,列車的運(yùn)行時(shí)間約束不僅對(duì)最小區(qū)間運(yùn)行時(shí)間進(jìn)行約束,也對(duì)最大區(qū)間運(yùn)行時(shí)間也進(jìn)行了約束,具體表示如下:

    (2)

    其中,Tr min和Tr max分別表示列車的最小區(qū)間運(yùn)行時(shí)間和最大區(qū)間運(yùn)行時(shí)間。

    3) 追蹤間隔時(shí)間約束

    為了保障安全性,防止相鄰列車發(fā)生追尾事件,規(guī)定列車追蹤間隔時(shí)間約束如下:

    (3)

    其中:Tg min表示列車間的最小追蹤間隔;Tfd min表示相鄰列車到發(fā)最小間隔。

    4) 到發(fā)車時(shí)間約束

    為防止發(fā)生延誤的列車對(duì)前面的列車產(chǎn)生影響,盡快實(shí)現(xiàn)“按圖行車”,提高調(diào)整效率,對(duì)列車到發(fā)車時(shí)間進(jìn)行約束,保證不會(huì)出現(xiàn)早于計(jì)劃時(shí)間到發(fā)車的情況:

    (4)

    1.2 優(yōu)化目標(biāo)

    列車的總晚點(diǎn)時(shí)間可以表示為各列車到達(dá)各車站的晚點(diǎn)時(shí)分和各列車從各車站出發(fā)的晚點(diǎn)時(shí)間之和受晚點(diǎn)影響的列車數(shù)Nl和到發(fā)時(shí)間線數(shù)Ns的改進(jìn)公式:

    (5)

    (6)

    其中:tdelay為列車晚點(diǎn)時(shí)間;Tyz為列車追蹤間隔冗余時(shí)間;Tyr為區(qū)間運(yùn)行冗余時(shí)間;N和M分別車站總數(shù)和列車總數(shù);n和m表示發(fā)生延誤的車站和列車;a為延誤類型標(biāo)志,a=0表示列車到達(dá)延誤,a=1表示列車發(fā)車延誤;ceil(x)和floor(x) 函數(shù)分別表示x趨近于正無(wú)窮的整數(shù)和趨近于負(fù)無(wú)窮大的整數(shù)。列車總晚點(diǎn)時(shí)間最小函數(shù)可以表示成:

    (7)

    由于存在列車到發(fā)車約束條件,所以上述公式可以簡(jiǎn)化為:

    (8)

    晚點(diǎn)時(shí)域內(nèi)乘客等待總時(shí)間最小函數(shù)表示如下:

    (9)

    (10)

    城市軌道交通列車運(yùn)行均衡調(diào)整優(yōu)化目標(biāo)函數(shù)用列車總晚點(diǎn)時(shí)間最小函數(shù)和晚點(diǎn)時(shí)域內(nèi)乘客等待總時(shí)間最小函數(shù)的加權(quán)和表示如下:

    F=α1×F1+α2×F2=

    (11)

    其中,α1和α2分別為總晚點(diǎn)時(shí)間最小函數(shù)和晚點(diǎn)時(shí)域內(nèi)乘客等待總時(shí)間最小函數(shù)的權(quán)重值,滿足α1+α2=1。

    2 人工魚群優(yōu)化算法

    2003年,李曉磊博士提出了人工魚群算法[15],它是一種模擬自然魚群行為的智能優(yōu)化算法,是一種自上而下的尋優(yōu)結(jié)構(gòu),其基本行為包括覓食、聚群、追尾和隨機(jī),逐步聚集在食物濃度最大的地方,通過(guò)模仿自然魚的行為更快地找到全局最優(yōu)值。為了兼顧列車自動(dòng)調(diào)整優(yōu)化效果與求解時(shí)間,本研究采用改進(jìn)的人工魚群算法對(duì)上述模型進(jìn)行求解。

    2.1 行為描述

    假設(shè)有fishnum條魚,每條魚的狀態(tài)Xi={x1,x2,…,xm},(i=1,2,…,fishnum)為待尋優(yōu)變量,食物濃度f(wàn)(Xi)為適應(yīng)度函數(shù),人工魚之間距離dij=||Xi-Xj||,visual為人工魚的視野,step為人工魚的步長(zhǎng),δ為人工魚擁擠度因子,try_number為嘗試次數(shù)。

    在城市軌道交通列車運(yùn)行調(diào)整模型中,目標(biāo)函數(shù)F>0恒成立,因此可以得到人工魚群的食物濃度函數(shù):

    1) 覓食行為:假設(shè)人工魚在t時(shí)刻的位置為Xi,在視野內(nèi)隨機(jī)選擇狀態(tài)Xj,如果f(Xi)

    在城市軌道交通列車運(yùn)行調(diào)整模型中,對(duì)前進(jìn)條件的要求較為嚴(yán)格,本研究提出改進(jìn)的覓食行為Matlab偽代碼如下:

    functiony=AF_prey(Xi)

    fori=1∶trynumber

    Xj=Xi+Rand*visual;

    iff(Xi)

    if judge(Xnew)==1

    y=Xnew

    break;

    else

    continue;

    end

    else

    continue;

    end

    end

    end

    其中:judge(Xnew)為判斷Xnew是否滿足約束條件的函數(shù),如果滿足約束條件,置1;反之,置0。

    2) 聚群行為:假設(shè)人工魚在t時(shí)刻的位置為Xi,計(jì)算視野內(nèi)滿足dijδf(Xi),則Xi向Xc方向按式(12)得到Xnew,判斷Xnew是否滿足要求,如果滿足則Xnew為下一個(gè)狀態(tài);否則執(zhí)行覓食行為。

    (12)

    3) 追尾行為:假設(shè)人工魚在t時(shí)刻的位置為Xi,搜索視野內(nèi)滿足dijδf(Xi),則按式(13)得到Xnew,判斷Xnew是否滿足要求,如果滿足則Xnew為下一個(gè)狀態(tài);否則執(zhí)行覓食行為。

    (13)

    4) 隨機(jī)行為:假設(shè)人工魚在t時(shí)刻的位置為Xi,當(dāng)覓食行為重復(fù)try_number次后,仍不滿足約束條件則執(zhí)行隨機(jī)行為:

    Xnew=Xi+Rand*step

    if judge(Xnew)==1

    y=Xnew;

    else

    y=Xi;

    2.2 求解流程

    本研究應(yīng)用改進(jìn)的人工魚群算法,對(duì)列車運(yùn)行的出發(fā)時(shí)間和到達(dá)時(shí)間進(jìn)行調(diào)整,求解目標(biāo)函數(shù)值最優(yōu)下的各列列車在各個(gè)站臺(tái)的合理到發(fā)時(shí)問(wèn),具體求解流程如下:

    步驟1參數(shù)初始化。初始化列車運(yùn)行圖XR、列車總數(shù)M、車站總數(shù)N、列車追蹤間隔冗余時(shí)間Tyz、區(qū)間運(yùn)行冗余時(shí)間Tyr和列車各約束條件Ts min、Ts max、Tr min、Tr max、Tg min和Tfd min,輸入列車晚點(diǎn)時(shí)間tdelay;

    其中:

    步驟2人工魚群參數(shù)初始化。初始化人工魚群種群數(shù)目fishnum、人工魚視野visual、人工魚步長(zhǎng)step、人工魚擁擠度因子δ、嘗試次數(shù)try_number和最大迭代次數(shù)Maxgen;

    步驟3根據(jù)式(5)和式(6)計(jì)算晚點(diǎn)時(shí)域參數(shù);

    步驟4根據(jù)晚點(diǎn)時(shí)域參數(shù),初始化人工魚的狀態(tài)Xi(i=1,2,…,fishnum),并計(jì)算對(duì)應(yīng)的食物濃度f(wàn)(Xi);

    其中

    由于城市軌道交通列車運(yùn)行調(diào)整模型約束條件復(fù)雜,因此為了保證初始化的人工魚的狀態(tài)Xi有足夠可行解,提出以下初始化規(guī)則:

    行規(guī)則:

    Xi(ii,jj+1)=Xi(ii,jj)+Tg min+k*rand

    列規(guī)則:

    Xi(ii+1,jj)=Xi(ii,jj)+Ts min+k*rand

    Xi(ii,jj)為到達(dá)線

    Xi(ii+1,jj)=Xi(ii,jj)+Tr min+k*rand

    Xi(ii,jj)為出發(fā)線

    其中,k是初始化參數(shù),根據(jù)列車晚點(diǎn)時(shí)間tdelay不同而調(diào)整。在進(jìn)行初始化時(shí),先初始化人工魚的狀態(tài)的第一行和第一列,然后按照行規(guī)則和列規(guī)則分別計(jì)算下一行的各個(gè)數(shù)據(jù),對(duì)比選擇既滿足行規(guī)則和列規(guī)則的計(jì)算結(jié)果,判斷結(jié)果是否滿足到發(fā)車約束,如果滿足則初始化成功;否則,選擇相應(yīng)的原到發(fā)車時(shí)間代替不滿足約束的數(shù)據(jù)。

    步驟5對(duì)人工魚的狀態(tài)Xi按食物濃度f(wàn)(Xi)降序排列,并設(shè)定全局最優(yōu)pbest(X1)=f(X1);

    步驟6按照人工魚的4種行為對(duì)人工魚的狀態(tài)Xi進(jìn)行更新,并選擇食物濃度最大的行為更新結(jié)果作為人工魚的狀態(tài)Xi;

    步驟7判斷是否滿足結(jié)束條件,如果滿足,進(jìn)行步驟8;否則,轉(zhuǎn)到步驟5繼續(xù)執(zhí)行;

    步驟8對(duì)人工魚的狀態(tài)Xi按食物濃度f(wàn)(Xi)降序排列,并輸出調(diào)整的列車運(yùn)行時(shí)刻以及對(duì)應(yīng)的全局最優(yōu)解pbest(X1);

    步驟9算法結(jié)束。

    算法流程如圖1所示。

    圖1 人工魚群算法求解流程框圖

    3 實(shí)例仿真

    3.1 基本參數(shù)

    南京地鐵3號(hào)線全長(zhǎng)44.9 km,其中地下段42.5 km、高架段2.4 km,共有林場(chǎng)站到秣周東路29座站臺(tái),本研究算法選擇三號(hào)線勝太西路站到新莊站路段早高峰(7∶15∶48—8∶24∶07)期間運(yùn)行的8列車,構(gòu)造列車計(jì)劃時(shí)刻表。列車區(qū)間運(yùn)行參數(shù)以及魚群算法相關(guān)參數(shù)如表1所示。

    表1 列車區(qū)間運(yùn)行參數(shù)以及魚群算法參數(shù)

    3.2 魚群算法仿真

    為了驗(yàn)證算法的適用性和有效性,假設(shè)第4輛車在經(jīng)過(guò)南京南站時(shí)由于隨機(jī)因素發(fā)生晚點(diǎn),根據(jù)式(5)和式(6),受晚點(diǎn)影響列車數(shù)Nl和到發(fā)時(shí)間線數(shù)Ns與追蹤間隔冗余時(shí)間Tyz有關(guān)系,故分別選取小于追蹤間隔冗余時(shí)間的晚點(diǎn)時(shí)間、大于1倍追蹤間隔冗余時(shí)間的晚點(diǎn)時(shí)間和大于2倍追蹤間隔冗余時(shí)間的晚點(diǎn)時(shí)間。

    1) 當(dāng)晚點(diǎn)時(shí)間小于追蹤間隔冗余時(shí)間Tyz時(shí),假設(shè)tdelay=40 s,根據(jù)式(5)和式(6)可以確定晚點(diǎn)時(shí)域受影響列車Nl=1和受影響到發(fā)時(shí)間線數(shù)Ns=5,根據(jù)魚群算法得到的適應(yīng)度收斂圖和列車最優(yōu)調(diào)整方案如圖2和圖3所示。

    圖2 tdelay=40 s時(shí)魚群收斂曲線

    圖3 tdelay=40 s時(shí)列車最優(yōu)調(diào)整方案曲線

    從圖2可以看出,改進(jìn)魚群算法的收斂曲線在前50代收斂速度較快,當(dāng)經(jīng)過(guò)60代之后收斂曲線接近水平,算法優(yōu)化結(jié)果為列車總晚點(diǎn)時(shí)間F1=53.055 6 s、列車晚點(diǎn)時(shí)域內(nèi)乘客總等待時(shí)間F2=84 526 s;從圖3可以看出,列車在南京南站到站晚點(diǎn)40 s時(shí),經(jīng)過(guò)本站發(fā)車時(shí)間調(diào)整,在明發(fā)廣場(chǎng)站就調(diào)整回列車計(jì)劃運(yùn)行圖,調(diào)整效率較高。

    2) 當(dāng)晚點(diǎn)時(shí)間大于追蹤間隔冗余時(shí)間Tyz時(shí),假設(shè)tdelay=80 s,根據(jù)式(5)和式(6)可以確定晚點(diǎn)時(shí)域受影響列車Nl=2和受影響到發(fā)時(shí)間線數(shù)Ns=10,應(yīng)用改進(jìn)魚群算法得到的適應(yīng)度收斂圖和列車最優(yōu)調(diào)整方案如圖4和圖5所示。

    圖4 tdelay=80 s時(shí)魚群收斂曲線

    圖5 tdelay=80 s時(shí)列車最優(yōu)調(diào)整方案曲線

    從圖4可以看出,改進(jìn)魚群收斂曲線在前20代收斂速度較快,當(dāng)經(jīng)過(guò)30代之后收斂曲線接近水平,算法優(yōu)化結(jié)果為列車總晚點(diǎn)時(shí)間F1=389.849 6 s、列車晚點(diǎn)時(shí)域內(nèi)乘客總等待時(shí)間F2=168 180 s;從圖5可以看出,由于列車在南京南站到站晚點(diǎn)80s大于間隔冗余時(shí)間Tyz,所以后續(xù)列車也受到影響,相繼偏離計(jì)劃運(yùn)行圖,在經(jīng)過(guò)改進(jìn)魚群算法調(diào)整后,延誤的車輛均較快調(diào)整到列車計(jì)劃運(yùn)行圖。

    3) 為了驗(yàn)證改進(jìn)魚群算法對(duì)列車發(fā)生程度較大的晚點(diǎn)的調(diào)整能力,假設(shè)tdelay=130 s,同理計(jì)算得到晚點(diǎn)時(shí)域受影響列車Nl=3和受影響到發(fā)時(shí)間線數(shù)Ns=15,根據(jù)魚群算法得到的適應(yīng)度收斂圖和列車最優(yōu)調(diào)整方案如圖6和圖7所示。

    圖6 tdelay=130 s時(shí)魚群收斂曲線

    圖7 tdelay=130 s時(shí)列車最優(yōu)調(diào)整方案曲線

    從圖6可以看出,改進(jìn)魚群算法的收斂曲線在前60代收斂速度較快,當(dāng)經(jīng)過(guò)60代之后收斂曲線接近水平,算法優(yōu)化結(jié)果為列車總晚點(diǎn)時(shí)間F1=1 344.8 s、列車晚點(diǎn)時(shí)域內(nèi)乘客總等待時(shí)間F2=8.884 4×105s;從圖7可以看出,在經(jīng)過(guò)改進(jìn)魚群算法調(diào)整后,發(fā)生較大程度晚點(diǎn)的列車也能較快地調(diào)整到列車計(jì)劃運(yùn)行圖,減少了列車間的延誤傳播。

    4 結(jié)論

    針對(duì)城市軌道交通列車運(yùn)行調(diào)整問(wèn)題,本研究建立了列車運(yùn)行均衡調(diào)整模型;將列車運(yùn)行調(diào)整的多約束條件作為人工魚群的前進(jìn)條件,提出了改進(jìn)人工魚群算法;選取南京地鐵三號(hào)線的列車運(yùn)行數(shù)據(jù),分別對(duì)晚點(diǎn)時(shí)間40 s、80 s和 130 s 進(jìn)行仿真,形成最優(yōu)調(diào)整方案。仿真結(jié)果表明,改進(jìn)的人工魚群算法可以滿足城市軌道交通列車運(yùn)行調(diào)整的實(shí)時(shí)調(diào)整要求,具有優(yōu)良的適應(yīng)性和高效性,減少了列車間的延誤傳播,保障了列車運(yùn)行效率和客運(yùn)服務(wù)質(zhì)量。

    猜你喜歡
    晚點(diǎn)魚群列車運(yùn)行
    基于馬爾科夫鏈的高鐵列車連帶晚點(diǎn)橫向傳播
    晚點(diǎn)的火車(外三首)
    金沙江文藝(2022年4期)2022-04-26 14:14:22
    改善地鐵列車運(yùn)行舒適度方案探討
    高速鐵路初始晚點(diǎn)致因-影響列車數(shù)分布模型
    魚群漩渦
    中外文摘(2017年19期)2017-10-10 08:28:41
    列車運(yùn)行控制系統(tǒng)技術(shù)發(fā)展趨勢(shì)分析
    基于改進(jìn)魚群優(yōu)化支持向量機(jī)的短期風(fēng)電功率預(yù)測(cè)
    基于人工魚群算法的光伏陣列多峰MPPT控制策略
    相同徑路的高速列車運(yùn)行圖編制方法
    多子群并行人工魚群算法的改進(jìn)研究
    自线自在国产av| 色视频在线一区二区三区| 两个人看的免费小视频| 热99国产精品久久久久久7| 亚洲中文字幕日韩| 国产xxxxx性猛交| 久久久久精品人妻al黑| 国产精品成人在线| 曰老女人黄片| 自线自在国产av| 蜜桃在线观看..| 国产成人精品久久二区二区91| 最黄视频免费看| 免费观看a级毛片全部| 亚洲一区中文字幕在线| 热99re8久久精品国产| 国产精品久久久人人做人人爽| 成年人免费黄色播放视频| 成年美女黄网站色视频大全免费| www.自偷自拍.com| 中文精品一卡2卡3卡4更新| 亚洲一码二码三码区别大吗| 国产福利在线免费观看视频| 国产福利在线免费观看视频| 大型av网站在线播放| 狠狠狠狠99中文字幕| 天堂中文最新版在线下载| 性色av乱码一区二区三区2| 国产男女超爽视频在线观看| 亚洲国产日韩一区二区| 纵有疾风起免费观看全集完整版| 欧美另类亚洲清纯唯美| 久久久久网色| 国产成人系列免费观看| 中文字幕人妻丝袜制服| 亚洲熟女精品中文字幕| 久久影院123| 99国产精品一区二区三区| av线在线观看网站| 宅男免费午夜| 欧美日韩成人在线一区二区| xxxhd国产人妻xxx| 嫩草影视91久久| 欧美黄色片欧美黄色片| 久9热在线精品视频| av一本久久久久| 国产精品1区2区在线观看. | 亚洲欧美成人综合另类久久久| 久久久精品94久久精品| 欧美成人午夜精品| 日韩欧美一区二区三区在线观看 | 人人妻人人澡人人看| 91成人精品电影| 少妇粗大呻吟视频| 国产深夜福利视频在线观看| 99久久人妻综合| 国产日韩欧美在线精品| 日韩欧美一区视频在线观看| 亚洲视频免费观看视频| 亚洲av日韩精品久久久久久密| 国产精品一区二区在线观看99| 日韩欧美免费精品| 国产黄色免费在线视频| 精品国产一区二区久久| 91精品国产国语对白视频| 久久久久久人人人人人| 成人免费观看视频高清| 少妇人妻久久综合中文| 在线 av 中文字幕| 亚洲成人国产一区在线观看| 又大又爽又粗| 我的亚洲天堂| 97在线人人人人妻| 亚洲人成77777在线视频| 亚洲欧美成人综合另类久久久| 亚洲国产欧美网| 十八禁高潮呻吟视频| 久久久国产精品麻豆| 黄色片一级片一级黄色片| 精品少妇内射三级| 国产精品麻豆人妻色哟哟久久| 国产伦人伦偷精品视频| 久久久久久亚洲精品国产蜜桃av| 欧美另类一区| 黄色毛片三级朝国网站| 精品人妻1区二区| 国产精品自产拍在线观看55亚洲 | cao死你这个sao货| 欧美精品av麻豆av| 亚洲精品国产av蜜桃| 美女高潮到喷水免费观看| 女警被强在线播放| 在线十欧美十亚洲十日本专区| 国产又爽黄色视频| 18在线观看网站| 高清av免费在线| 一级a爱视频在线免费观看| 国产免费av片在线观看野外av| 亚洲一区二区三区欧美精品| 国产人伦9x9x在线观看| 天堂俺去俺来也www色官网| 亚洲男人天堂网一区| 波多野结衣一区麻豆| 午夜福利视频精品| 无遮挡黄片免费观看| 亚洲七黄色美女视频| 大码成人一级视频| 狂野欧美激情性bbbbbb| 欧美黑人欧美精品刺激| 99热网站在线观看| 97精品久久久久久久久久精品| 免费在线观看完整版高清| 欧美日韩视频精品一区| 精品国产国语对白av| 老熟女久久久| 亚洲激情五月婷婷啪啪| 啦啦啦在线免费观看视频4| 精品亚洲乱码少妇综合久久| 在线精品无人区一区二区三| 亚洲国产av影院在线观看| 91精品国产国语对白视频| 国产野战对白在线观看| 人成视频在线观看免费观看| 色视频在线一区二区三区| 亚洲av成人不卡在线观看播放网 | 久久久久网色| 精品乱码久久久久久99久播| 精品国产乱码久久久久久小说| 亚洲国产欧美网| 无遮挡黄片免费观看| 美女大奶头黄色视频| 国产激情久久老熟女| 亚洲欧美日韩另类电影网站| 亚洲精品美女久久av网站| 久久这里只有精品19| 人妻人人澡人人爽人人| 午夜福利,免费看| 国产极品粉嫩免费观看在线| 欧美激情久久久久久爽电影 | 一区二区三区精品91| 亚洲色图综合在线观看| 国产深夜福利视频在线观看| 久久天躁狠狠躁夜夜2o2o| 国产精品欧美亚洲77777| 亚洲精品中文字幕在线视频| 51午夜福利影视在线观看| √禁漫天堂资源中文www| 日韩中文字幕欧美一区二区| 黄色视频,在线免费观看| 亚洲国产欧美在线一区| 97人妻天天添夜夜摸| 国产伦理片在线播放av一区| 夜夜夜夜夜久久久久| 国产精品久久久久久精品电影小说| 亚洲精品一二三| 别揉我奶头~嗯~啊~动态视频 | 久久精品成人免费网站| 午夜日韩欧美国产| 在线永久观看黄色视频| 精品一区二区三区四区五区乱码| 一级a爱视频在线免费观看| 亚洲国产日韩一区二区| 法律面前人人平等表现在哪些方面 | 国产精品影院久久| 1024视频免费在线观看| 久久精品国产综合久久久| 午夜福利乱码中文字幕| 国产麻豆69| 丝袜脚勾引网站| 一二三四社区在线视频社区8| 日韩制服骚丝袜av| 日韩 亚洲 欧美在线| 亚洲精品粉嫩美女一区| 国产成人啪精品午夜网站| 不卡一级毛片| 91成年电影在线观看| 欧美黄色淫秽网站| 久久久久久久久免费视频了| 无限看片的www在线观看| 国产亚洲午夜精品一区二区久久| 国产精品亚洲av一区麻豆| 成人三级做爰电影| 久久影院123| 曰老女人黄片| 午夜免费成人在线视频| 久久久欧美国产精品| 俄罗斯特黄特色一大片| 国产熟女午夜一区二区三区| 国产亚洲欧美精品永久| 精品国内亚洲2022精品成人 | 欧美日本中文国产一区发布| 国产黄色免费在线视频| 正在播放国产对白刺激| 99热全是精品| 成年美女黄网站色视频大全免费| 日日摸夜夜添夜夜添小说| 免费少妇av软件| 亚洲性夜色夜夜综合| 夜夜夜夜夜久久久久| 国产亚洲av高清不卡| 国产精品九九99| 高清av免费在线| 久久久久久久大尺度免费视频| 国产成人精品无人区| 精品国产乱子伦一区二区三区 | 99国产综合亚洲精品| 亚洲专区中文字幕在线| 国产真人三级小视频在线观看| 国产在线免费精品| 免费在线观看黄色视频的| 日韩中文字幕视频在线看片| 脱女人内裤的视频| 亚洲中文字幕日韩| av视频免费观看在线观看| 男女免费视频国产| av片东京热男人的天堂| 久久99热这里只频精品6学生| 午夜两性在线视频| h视频一区二区三区| 亚洲精品中文字幕一二三四区 | 精品熟女少妇八av免费久了| 免费在线观看日本一区| 亚洲专区国产一区二区| 亚洲国产日韩一区二区| www.熟女人妻精品国产| 女人高潮潮喷娇喘18禁视频| 亚洲国产欧美一区二区综合| 午夜91福利影院| 蜜桃在线观看..| 日本撒尿小便嘘嘘汇集6| 777久久人妻少妇嫩草av网站| 国产欧美亚洲国产| 正在播放国产对白刺激| 国产淫语在线视频| 精品国产乱码久久久久久小说| 97人妻天天添夜夜摸| 精品久久久久久久毛片微露脸 | 国产精品麻豆人妻色哟哟久久| 欧美日韩视频精品一区| 又紧又爽又黄一区二区| 亚洲欧美精品自产自拍| 精品卡一卡二卡四卡免费| 日韩欧美一区二区三区在线观看 | 久久人妻熟女aⅴ| 天堂8中文在线网| 国产日韩欧美在线精品| 国产区一区二久久| 999久久久国产精品视频| 丝袜脚勾引网站| 人妻一区二区av| 久久久久国产一级毛片高清牌| 另类精品久久| 亚洲激情五月婷婷啪啪| a级毛片黄视频| 亚洲精品美女久久久久99蜜臀| 99热网站在线观看| 嫩草影视91久久| 亚洲av欧美aⅴ国产| 久久影院123| 视频区欧美日本亚洲| 各种免费的搞黄视频| 老熟妇仑乱视频hdxx| 国产成人精品久久二区二区91| 亚洲中文日韩欧美视频| 18禁国产床啪视频网站| 亚洲国产中文字幕在线视频| av线在线观看网站| 亚洲天堂av无毛| 精品福利观看| 亚洲欧美一区二区三区黑人| 午夜91福利影院| 久久精品熟女亚洲av麻豆精品| 男女边摸边吃奶| 多毛熟女@视频| 曰老女人黄片| 蜜桃在线观看..| 中文欧美无线码| 50天的宝宝边吃奶边哭怎么回事| 日韩 亚洲 欧美在线| 少妇被粗大的猛进出69影院| av超薄肉色丝袜交足视频| 亚洲精品国产色婷婷电影| 国产视频一区二区在线看| 19禁男女啪啪无遮挡网站| 日本vs欧美在线观看视频| 成年人午夜在线观看视频| 一区二区日韩欧美中文字幕| 汤姆久久久久久久影院中文字幕| 99久久综合免费| 男女床上黄色一级片免费看| 国产高清国产精品国产三级| 久久精品aⅴ一区二区三区四区| 成人国产av品久久久| tube8黄色片| 视频区图区小说| 不卡一级毛片| 国产一区二区三区av在线| 久久精品亚洲熟妇少妇任你| kizo精华| 汤姆久久久久久久影院中文字幕| 日本av手机在线免费观看| 男人舔女人的私密视频| 最近中文字幕2019免费版| 少妇的丰满在线观看| 97人妻天天添夜夜摸| 亚洲精品美女久久久久99蜜臀| 久久这里只有精品19| 久久久久网色| 高清视频免费观看一区二区| 一区二区av电影网| 色婷婷久久久亚洲欧美| 黄片小视频在线播放| 女警被强在线播放| 欧美国产精品va在线观看不卡| 国产日韩欧美在线精品| 亚洲黑人精品在线| 丁香六月欧美| 久久精品熟女亚洲av麻豆精品| 国产97色在线日韩免费| www.自偷自拍.com| 欧美精品亚洲一区二区| a级片在线免费高清观看视频| 亚洲综合色网址| 女人高潮潮喷娇喘18禁视频| 两人在一起打扑克的视频| 中文字幕制服av| 91国产中文字幕| 水蜜桃什么品种好| 亚洲国产精品一区二区三区在线| 老司机深夜福利视频在线观看 | 国产精品久久久久久精品电影小说| 免费观看人在逋| 男男h啪啪无遮挡| 国产欧美亚洲国产| 久久午夜综合久久蜜桃| 黄片小视频在线播放| 亚洲精品国产色婷婷电影| av欧美777| 午夜激情久久久久久久| 在线观看免费视频网站a站| 久久久久国产一级毛片高清牌| 日本91视频免费播放| 熟女少妇亚洲综合色aaa.| 国产精品影院久久| 最新在线观看一区二区三区| 亚洲欧美色中文字幕在线| 精品少妇一区二区三区视频日本电影| 母亲3免费完整高清在线观看| 激情视频va一区二区三区| 欧美激情 高清一区二区三区| 在线观看免费午夜福利视频| 一本久久精品| 三上悠亚av全集在线观看| 肉色欧美久久久久久久蜜桃| 亚洲avbb在线观看| tube8黄色片| av网站免费在线观看视频| 亚洲国产精品一区三区| 久久久久久久国产电影| 亚洲欧美一区二区三区黑人| 黑人欧美特级aaaaaa片| 国产在线一区二区三区精| 国产av精品麻豆| 如日韩欧美国产精品一区二区三区| 国产精品一区二区免费欧美 | 国产精品久久久久久人妻精品电影 | 天天躁夜夜躁狠狠躁躁| 亚洲国产毛片av蜜桃av| 国产av又大| 久久国产精品影院| 久久精品国产亚洲av高清一级| 啦啦啦视频在线资源免费观看| 亚洲精品美女久久久久99蜜臀| 欧美久久黑人一区二区| 欧美变态另类bdsm刘玥| 一级片'在线观看视频| 国产一区有黄有色的免费视频| 两人在一起打扑克的视频| 久久久久久久国产电影| 日本撒尿小便嘘嘘汇集6| 好男人电影高清在线观看| 三级毛片av免费| 国产精品久久久久久精品古装| 最新的欧美精品一区二区| 亚洲自偷自拍图片 自拍| 欧美黑人精品巨大| 欧美精品啪啪一区二区三区 | 久久久久国产一级毛片高清牌| 亚洲av成人不卡在线观看播放网 | 精品亚洲成a人片在线观看| 一区二区av电影网| 新久久久久国产一级毛片| 日本一区二区免费在线视频| 亚洲专区中文字幕在线| 国产成人免费观看mmmm| 亚洲第一av免费看| av网站在线播放免费| 久久狼人影院| 亚洲欧洲日产国产| 1024香蕉在线观看| 人人妻人人添人人爽欧美一区卜| 免费在线观看完整版高清| 国产成人av激情在线播放| 亚洲,欧美精品.| 久久人人爽av亚洲精品天堂| 国产在线视频一区二区| 91精品伊人久久大香线蕉| 久久国产精品大桥未久av| 中文字幕高清在线视频| 国产精品一区二区在线观看99| 丰满少妇做爰视频| 中文字幕人妻丝袜制服| av不卡在线播放| 欧美在线一区亚洲| 手机成人av网站| 天天操日日干夜夜撸| 精品福利观看| 国产黄色免费在线视频| 国产亚洲一区二区精品| 免费在线观看黄色视频的| 国产一级毛片在线| 国产精品秋霞免费鲁丝片| 国产亚洲av片在线观看秒播厂| 精品国产一区二区久久| 婷婷色av中文字幕| 精品久久久久久电影网| 99久久人妻综合| e午夜精品久久久久久久| 男女午夜视频在线观看| 最新在线观看一区二区三区| 精品国产一区二区三区四区第35| 国产成人啪精品午夜网站| 一区二区日韩欧美中文字幕| 两个人看的免费小视频| 十八禁网站网址无遮挡| 咕卡用的链子| 亚洲欧美一区二区三区黑人| 侵犯人妻中文字幕一二三四区| 亚洲av国产av综合av卡| 久久久精品区二区三区| 精品一区二区三区四区五区乱码| 91成人精品电影| 自拍欧美九色日韩亚洲蝌蚪91| 两个人免费观看高清视频| 高清在线国产一区| 19禁男女啪啪无遮挡网站| 91精品国产国语对白视频| 18禁黄网站禁片午夜丰满| 在线看a的网站| 精品少妇久久久久久888优播| 99久久99久久久精品蜜桃| kizo精华| 别揉我奶头~嗯~啊~动态视频 | 中文字幕制服av| 久久国产精品大桥未久av| 青春草亚洲视频在线观看| 三上悠亚av全集在线观看| 这个男人来自地球电影免费观看| 国产免费视频播放在线视频| 欧美老熟妇乱子伦牲交| 国产成人一区二区三区免费视频网站| 97精品久久久久久久久久精品| 肉色欧美久久久久久久蜜桃| 亚洲 欧美一区二区三区| 啦啦啦免费观看视频1| 免费观看a级毛片全部| 久久久久网色| 亚洲国产日韩一区二区| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美色中文字幕在线| 99久久国产精品久久久| 国产亚洲av片在线观看秒播厂| 高清av免费在线| 巨乳人妻的诱惑在线观看| 免费在线观看黄色视频的| 精品国产乱码久久久久久小说| 91九色精品人成在线观看| 十八禁网站免费在线| 高清在线国产一区| tube8黄色片| 久久精品久久久久久噜噜老黄| 侵犯人妻中文字幕一二三四区| 叶爱在线成人免费视频播放| 国产成人精品久久二区二区免费| 日本精品一区二区三区蜜桃| √禁漫天堂资源中文www| 亚洲成人免费av在线播放| 一进一出抽搐动态| 久久久精品区二区三区| 成年美女黄网站色视频大全免费| 天天躁日日躁夜夜躁夜夜| 成人免费观看视频高清| 50天的宝宝边吃奶边哭怎么回事| 亚洲第一青青草原| 欧美国产精品一级二级三级| videos熟女内射| 精品熟女少妇八av免费久了| 国产精品一二三区在线看| 亚洲一码二码三码区别大吗| 无限看片的www在线观看| 精品国产一区二区三区久久久樱花| 淫妇啪啪啪对白视频 | 国产一区二区三区综合在线观看| 精品高清国产在线一区| 国产精品 国内视频| 曰老女人黄片| 欧美成狂野欧美在线观看| 国产精品影院久久| 美女扒开内裤让男人捅视频| 宅男免费午夜| 一区二区三区乱码不卡18| 中文字幕另类日韩欧美亚洲嫩草| 黄色视频在线播放观看不卡| 国产色视频综合| 国产欧美日韩一区二区三区在线| 最新在线观看一区二区三区| avwww免费| 亚洲av电影在线观看一区二区三区| 国产精品.久久久| 夜夜骑夜夜射夜夜干| 一级片免费观看大全| 水蜜桃什么品种好| 亚洲av国产av综合av卡| 悠悠久久av| 少妇的丰满在线观看| 深夜精品福利| 高清欧美精品videossex| 超色免费av| 久久精品国产a三级三级三级| 亚洲国产中文字幕在线视频| 亚洲精品粉嫩美女一区| 成年人午夜在线观看视频| 国产一区二区三区av在线| 亚洲成人免费av在线播放| 欧美久久黑人一区二区| 黄频高清免费视频| 超碰97精品在线观看| 日韩中文字幕视频在线看片| 一本一本久久a久久精品综合妖精| av电影中文网址| 桃红色精品国产亚洲av| 久久精品国产a三级三级三级| 国产精品欧美亚洲77777| 午夜福利免费观看在线| 老司机亚洲免费影院| 69精品国产乱码久久久| 欧美精品高潮呻吟av久久| 黄色怎么调成土黄色| 成年美女黄网站色视频大全免费| 国产一级毛片在线| 亚洲综合色网址| 国产精品九九99| 老司机在亚洲福利影院| 久久精品久久久久久噜噜老黄| 精品视频人人做人人爽| 中文字幕av电影在线播放| 人成视频在线观看免费观看| 三上悠亚av全集在线观看| 黑人猛操日本美女一级片| 成人亚洲精品一区在线观看| 少妇 在线观看| 亚洲伊人久久精品综合| 国产伦理片在线播放av一区| 在线观看舔阴道视频| 亚洲五月婷婷丁香| 久久国产精品男人的天堂亚洲| 久久综合国产亚洲精品| 国产成+人综合+亚洲专区| 亚洲av国产av综合av卡| 日本91视频免费播放| 精品国产乱子伦一区二区三区 | 97人妻天天添夜夜摸| 黄色片一级片一级黄色片| 日韩电影二区| 脱女人内裤的视频| a级毛片在线看网站| 日韩欧美国产一区二区入口| 久久九九热精品免费| 欧美日韩亚洲高清精品| av福利片在线| 黄网站色视频无遮挡免费观看| 精品国产乱码久久久久久男人| 色婷婷久久久亚洲欧美| 久久99一区二区三区| 亚洲五月色婷婷综合| 两性夫妻黄色片| 狠狠精品人妻久久久久久综合| 99久久精品国产亚洲精品| 欧美av亚洲av综合av国产av| 美女扒开内裤让男人捅视频| 亚洲精品一二三| 国产av国产精品国产| 国产男女超爽视频在线观看| 女人爽到高潮嗷嗷叫在线视频| 国产精品自产拍在线观看55亚洲 | 免费观看a级毛片全部| 韩国精品一区二区三区| 爱豆传媒免费全集在线观看| 国产真人三级小视频在线观看| 脱女人内裤的视频| 视频在线观看一区二区三区| 精品亚洲成国产av| 国精品久久久久久国模美| 亚洲av国产av综合av卡| 无遮挡黄片免费观看| 丁香六月欧美| 又大又爽又粗| 国产真人三级小视频在线观看| 午夜日韩欧美国产| 丝袜人妻中文字幕| 欧美在线黄色| 啦啦啦视频在线资源免费观看| a在线观看视频网站| 亚洲欧美色中文字幕在线| 国产成人免费无遮挡视频|