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

    基于ESGA的純電動需求響應公交路徑規(guī)劃

    2023-01-03 05:55:54郭獻超夏戈松徐士偉陳興博
    甘肅科學學報 2022年6期
    關(guān)鍵詞:停靠站公交車公交

    郭獻超,夏戈松,徐士偉,陳興博

    (1.廣州市交通規(guī)劃研究院,廣東 廣州 510000; 2.華南理工大學土木與交通學院,廣東 廣州 510000; 3.廣東省交通運輸規(guī)劃研究中心,廣東 廣州 510000)

    近年來隨著中國經(jīng)濟的快速發(fā)展和城鎮(zhèn)化的加快,機動車輛迎來了爆發(fā)式增長,大量農(nóng)村人口進入城市,城市的交通壓力越來越大,雖然城市公共交通進入快速發(fā)展時期,但仍無法滿足市民便捷出行的需求。如今,我國的公共交通已經(jīng)進入了高速發(fā)展軌道,出現(xiàn)了各種各樣的交通方式,乘客的出行需求也各不相同,越來越傾向個性化發(fā)展的趨勢。需求響應公交的出現(xiàn)滿足了許多乘客的要求,迎來了快速發(fā)展的契機。通過對純電動需求響應公交路徑規(guī)劃的研究,可以進一步完善我國公共交通體系,增加公共交通的多樣性,擴大公共交通的覆蓋面和提高乘客對公共交通的滿意度。

    我國對需求響應公交的研究開始的較晚,國內(nèi)缺乏需求響應公交的研究內(nèi)容,尤其是純電動騎車在需求響應公交的應用研究。Lajunen[1]模擬了不同公交線路,通過對混合動力公交和純電動公交車能量消耗分析,計算出了它們在運營中的成本效益;李軍等[2]根據(jù)純電動公交車的特點,考慮了充電時間、電池狀態(tài)、充電速率等因素,設(shè)計了以營運車輛數(shù)最小為目標的純電動公交車輛調(diào)度算法;楊揚等[3]應用整數(shù)規(guī)劃、網(wǎng)絡流等基本理論,將純電動公交調(diào)度計劃問題轉(zhuǎn)換為網(wǎng)絡模型,通過使用列生成方法對模型進行求解;Schilde等[4]研究了不同速度下對DRT服務水平的影響;潘述亮等[5]對包含特殊需求的DRT模型進行了研究;盧小林等[6]提出一種以最小化公交車運營時間為目標的DRT路徑優(yōu)化模型;Masmoudi等[7]研究了異質(zhì)車輛的DAR的路徑規(guī)劃問題,提出了一種針對已知DAR問題特性的混合遺傳算法;Masmoudi等[8]建立了以總路徑成本最小化為目標的混合整數(shù)非線性規(guī)劃模型;鄭漢等[9]構(gòu)建需求響應公交服務定制等價分解模型,采用分布式列生成算法對模型進行求解;柳伍生等[10]提出了多需求響應機制下的城市綠色定制公交線網(wǎng)優(yōu)化模型,利用上車線網(wǎng)和下車線網(wǎng)相結(jié)合的分層線網(wǎng)算法對模型進行求解。目前國內(nèi)外研究主要集中在傳統(tǒng)公交的研究中,隨著電動車的推廣,許多城市都將常規(guī)公交換成了電動公交,然而有關(guān)電動公交應用在DRT的研究卻很少。

    純電動需求響應公交路徑規(guī)劃可分為臨時??空疽?guī)劃和路徑規(guī)劃兩階段,首先規(guī)劃純電動公交臨時停靠站點,通常通過K-means算法實現(xiàn)[11]。對公交路徑規(guī)劃可看作是一種特殊形式的帶有容量約束的車輛路徑問題(CVRP,capacitated vehicle routing problem),不同的是公交車需要從同一場站出發(fā)到達同一目的地,目的地與出發(fā)地往往不同。車輛路徑規(guī)劃已被證實為NP(non-deterministic polynomial)完全問題,通常設(shè)計啟發(fā)式算法[12-14]求解。遺傳算法(GA,genetic algorithm)是通過模擬生物在自然環(huán)境中的遺傳和進化過程而設(shè)計,是一種廣泛使用的啟發(fā)式搜索算法,其選擇策略通常是基于輪盤賭(RWSGA,roulette wheel selection genetic algorithm)的概率選擇。為了提高算法收斂速度和搜索結(jié)果,提出以精英選擇遺傳算法[15](ESGA,elitist selection genetic algorithm)求解純電動需求響應公交路徑規(guī)劃問題,并通過實例驗證猜想。

    1 問題描述和模型建立

    1.1 問題描述

    純電動公交車需求響應式公交路徑優(yōu)化問題可描述如下:在研究的區(qū)域內(nèi),通過調(diào)度公交場站若干輛相同的純電動公交車以滿足該區(qū)域乘客的出行需求,所有乘客會在預先規(guī)劃好的區(qū)域上車,并到達同一終點站。通過合理規(guī)劃公交車的行駛路徑,以實現(xiàn)企業(yè)利潤的最大化。該問題屬于帶有軟時間窗的車輛路徑問題(VRPSTW),主要分兩個部分進行處理。

    第一部分為需求響應站點規(guī)劃,根據(jù)乘客的響應點的位置分布,規(guī)劃若干需求響應站點,用于使需求響應公交車停車接收上車的乘客,并根據(jù)站點乘客的期望上車時間分布,確定需求響應站點的時間窗。

    第二部分為純電動公交車路徑規(guī)劃,根據(jù)已規(guī)劃的需求響應站點,規(guī)劃需求響應公交車的行駛路徑,這一部分可以看作帶有軟時間窗的車輛路徑問題,但是由于純電動公交車的特殊性,在路徑規(guī)劃的過程中,需要根據(jù)純電動公交車的續(xù)航狀態(tài),合理安排公交車到充電樁的充電策略。為了方便模型建立和求解,有以下基本假設(shè):

    (1) 假設(shè)純電動公交車不會載送乘客駛?cè)氤潆姌冻潆?

    (2) 假設(shè)配送車輛為相同車型的純電動公交車,最大載客量相同、單次充電行駛里程相同;

    (3) 假設(shè)配送車輛為相同車型的純電動公交車,最大載客量相同、單次充電行駛里程相同;

    (4) 假設(shè)從場站發(fā)出的純電動公交車只有兩種:第1種是純電動公交車發(fā)出前已經(jīng)充滿電(即k≤P0);第2種是純電動公交車需要到場站附近的充電樁(即h=n+1)完成快速充電后再發(fā)車(即k>P0);

    (5) 假設(shè)在所研究的網(wǎng)絡上分布若干充電樁能夠滿足純電動公交車的充電需求;

    (6) 假設(shè)每輛公交車都只從固定的站場出發(fā),完成所有乘客的接收后,到達統(tǒng)一的目的地;

    (7) 假設(shè)公交車只在臨時??空就\?每一區(qū)域的乘客需要到最近的臨時??空境塑?

    (8) 假設(shè)每輛公交車只開一條線路,每個停靠點只能被公交車訪問一次。

    基于以上對純電動需求響應式公交路徑規(guī)劃問題的描述和對模型構(gòu)建的假設(shè),為方便建立純電動需求響應公交路徑規(guī)劃的數(shù)學模型,定義以下數(shù)學符號的含義:

    p0表示公交場站充電樁的總數(shù)量;

    V表示純電動公交車的集合,v={k},k=1,2,…,p,當k≤p0時,表示純電動公交車在公交場站已充滿電,即滿電量發(fā)車;當k>p0時,表示純電動公交車需要去公交場站臨近的充電樁充滿電后再發(fā)車;

    S表示純電動公交車??空镜暮霞?v={i},i=1,2,…,n;

    E表示路網(wǎng)中充電樁的合集,E={j},j=1,2,…,m,j=1表示公交場站附近的充電樁,即k>p0的車輛需要先到該充電樁充電,再去接送乘客;

    A表示路網(wǎng)中所有研究點的合集,即充電樁總數(shù)量、純電動公交車??空究倲?shù)量、出發(fā)地、目的地的四者數(shù)量之和,A={h},h=0,1,2,…,n,n+1,…,n+m,n+m+1,h=0表示純電動公交車出發(fā)地,h=n+m+1表示純電動公交車目的地;

    O表示乘客的集合,O={a},a=1,2,…,b;

    Rk表示所有車輛的單次充電的行駛路徑集;

    Lk表示所有車輛完成響應需求的行駛路徑集;

    Cf表示純電動公交車的啟動成本,其與車輛行駛里程無關(guān),而與車輛開行數(shù)量有關(guān);

    C0表示純電動公交車的運行成本,其與車輛單位行駛里程成正比;

    ETi表示??空军ci的時間窗上界;

    LTi表示??空军ci的時間窗下界;

    C0表示純電動公交車在充電樁充電的時間;

    s表示純電動公交車在行駛過程中的平均行駛速度;

    Ct表示單位時間的懲罰成本;

    ua表示乘客a使用需求響應公交車的意愿,ua=1表示乘客具有使用意愿,否則ua=0,a=1,2,…,b;

    ea表示乘客a支付需求響應公交的意愿(元);

    e0表示企業(yè)響應乘客最低支付意愿(元);

    ra表示企業(yè)對乘客a出行需求的響應意愿,ra=1表示企業(yè)響應乘客需求,否則ra=0;

    dij表示點i到點j的行駛距離(i,j=0,1,2,…,n,n+1,…,n+m,n+m+1);

    d0表示純電動公交車充滿電后的最大行駛里程;

    xijk=1表示車輛k從點i行駛到站點j,否則xijk=0(i,j=0,1,2,…,n,n+1,…,n+m,n+m+1,k=1,2,…,p);

    Qi表示第i個??空军c的上車人數(shù),i=1,2,…,n;

    Q0表示需求響應公交車的最大載客量。

    1.2 模型建立

    純電動需求響應公交路徑規(guī)劃模型的目標函數(shù)包括總利潤和時間懲罰成本兩部分,純電動需求響應公交的成本包括運營車輛的人力成本和可變成本,收入為乘客支付的乘車票價總和,二者之差為企業(yè)的總利潤。將公交車超出時間窗的時間量化為時間函數(shù),并通過懲罰函數(shù)計入到總優(yōu)化目標中。綜上所述,純電動需求響應公交路徑規(guī)劃模型如下:

    (1)

    (2)

    當xijk=1時,

    (3)

    p≤pmax,

    (4)

    (5)

    (6)

    (7)

    (8)

    (9)

    (10)

    (11)

    (12)

    xijk=0或1,

    (13)

    k>p0?x0,n+1,k=1,

    (14)

    式(2)表示純電動公交車到達臨時停靠站i超出時間窗的偏離值;式(3)表示純電動公交車到達臨時??空緄的行程時間;式(4)表示企業(yè)可投入用于響應乘客需求的純電動公交車數(shù)量不大于場站公交車的總數(shù)量;式(5)表示純電動公交車只能在接送乘客之前駛?cè)氤潆姌冻潆?式(6)表示每輛純電動公交車的乘客數(shù)量不大于公交車的載客上限;式(7)表示每輛純電動公交車在每次充電后的行駛里程不超過其最大行駛里程;式(8)表示運營企業(yè)只響應支付意愿大于等于e0的乘客出行需求;式(10)表示每個臨時??空军c只有唯一的一輛純電動公交車駛?cè)?式(11)表示每個臨時??空军c只有唯一的一輛純電動公交車駛出;式(12)表示每輛純電動公交車的起點都是停車場站,終點都是乘客所要到達的目的地;式(13)表示決策變量被限制為0、1變量;式(14)表示對于未在公交場站充電的公交車,需要先去場站附近的充電樁快速充電,再去接送乘客。

    2 模型求解

    2.1 基于K-means算法的??空疽?guī)劃

    首先基于K-means聚類算法規(guī)劃了公交臨時??空尽-means聚類算法可以根據(jù)樣本之間的內(nèi)在聯(lián)系,將相同的樣本劃分為一類,從而將樣本整體劃分為有不同特征的若干類。K-means聚類的基本流程是先隨機選取k個樣本作為初始的類中心點。然后計算每個樣本與各個種類中心點的距離,把每個樣本劃分為與其距離最近的類中心點。當所有樣本對象都劃分完成后,每一類的類中心點會根據(jù)類中現(xiàn)有的樣本重新計算,這個過程一直重復直到?jīng)]有中心再發(fā)生變化為止。K-means算法的基本流程如圖1所示。

    圖1 K-means算法流程Fig.1 K-means algorithm flow

    通過選擇歐式距離作為樣本點到類中心點的判別距離,其計算公式為

    (15)

    選擇點合集的重心作為每一類的中心點,其計算公式為

    (16)

    2.2 基于ESGA的路徑規(guī)劃

    采用精英種群選擇的遺傳算法(ESGA,elitist selection genetic algorithm)來求解純電動需求響應公交路徑規(guī)劃模型,相比于常規(guī)的遺傳算法,ESGA用每一代種群中適應度較高的個體建立精英種群,精英種群不經(jīng)過選擇的操作直接進入子代種群。因此ESGA相對于常規(guī)的遺傳算法收斂速度更快,能夠更快地獲得求解結(jié)果。ESGA流程如下:

    (1) 編碼 編碼是為了將需要解決的實際問題轉(zhuǎn)化為可被遺傳算法直接求解的數(shù)學形式。在現(xiàn)代生物學中,個體的基因代表著不同個體特征。相似的,在遺傳算法中經(jīng)過基因編碼的個體代表優(yōu)化目標的一組可行解。常用的編碼方案主要有0、1編碼和整數(shù)編碼兩種,本文采用整數(shù)編碼方案,并根據(jù)純電動公交車行駛路徑的實際問題設(shè)計了編碼方案。

    為了解釋ESGA中編碼過程,通過實例來具體描述。首先,假設(shè)路網(wǎng)中有10個規(guī)劃的公交車臨時??空竞?個充電樁,此外還有公交車起點和目的地2個點。編碼的具體流程如下:

    步驟1:對路網(wǎng)中的點進行編號,0表示公交車起點,1~10表示臨時??空?11~14表示充電樁;

    步驟2:隨機生成1~10不重復的數(shù)列,表示公交車到達臨時??空镜捻樞?

    步驟3:在上述數(shù)列隨機插入若干個公交車的出發(fā)點0,表示需要安排調(diào)度的公交車數(shù)量;

    步驟4:在上述的列序號0后面隨機插入若干充電樁編碼,表示純電動公交車需要駛?cè)氤潆姌冻潆?

    步驟5:將上述的數(shù)列作為個體的基因編碼;

    步驟6:重復步驟2~步驟5,形成一定規(guī)模的種群。

    上述過程中個體基因編碼的示意圖如圖2所示。

    圖2 染色體編碼Fig.2 Chromosome coding

    (2) 適應度計算 ESGA通過選擇算子淘汰劣勢個體,保留優(yōu)勢個體,適應度為個體選擇的依據(jù),目標函數(shù)以最小化企業(yè)利潤的相反數(shù)為優(yōu)化目標,因此,使得目標函數(shù)值較小的個體具有較高的適應度。以目標函數(shù)的倒數(shù)作為計算個體的適應度函數(shù),即

    (3) 基于精英保留的選擇算子 通過采用了精英個體保留的選擇策略,在每一代種群中選擇一部分適應度最高的個體,不進行任何操作而直接保留到下一代,剩下部分的個體再通過選擇得到。保留的精英個體只占種群的小部分,這是為了避免種群陷入局部最優(yōu)解。精英個體保留的選擇策略可以保證種群的適應度不會增加,加快收斂速度,從而提高最優(yōu)解的搜索速度。

    (4) 交叉算子 遺傳算法使用交叉算子通過模擬種群進化中染色體的交叉重組,從而產(chǎn)生子代個體。ESGA采用的交叉算子操作流程如圖3所示。

    圖3 交叉操作Fig.3 Cross operation

    (5) 變異算子 遺傳算法的變異算子通過模擬種群進化中個體的基因發(fā)生突變,進而產(chǎn)生新的個體。ESGA采用兩點變異,操作流程如圖4所示。

    圖4 變異算子Fig.4 Mutation operator

    (6) 反轉(zhuǎn)算子 反轉(zhuǎn)算子從本質(zhì)上來說也是基因突變的一種類型,通過逆轉(zhuǎn)染色體基因片段的基因位順序,從而實現(xiàn)改變?nèi)旧w的基因。反轉(zhuǎn)算子也可以提高遺傳算法搜索能力。ESGA采用的反轉(zhuǎn)算子操作流程如圖5所示。

    圖5 反轉(zhuǎn)算子Fig.5 Inversion operator

    (7) 進化終止條件 從理論上來說,當搜索到最優(yōu)解時,即算法達到100%收斂時,算法會自動終止,然而在實際情況下,由于算法會陷入局部最優(yōu)很難到達100%收斂。因此,需要設(shè)計終止條件以避免算法進入無限迭代循環(huán)。遺傳算法的終止條件判別主要有以下兩種方法:

    ① 判定相鄰兩代種群適應度的差值小于某一閾值時終止進化。即當種群進化速度緩慢,種群進化進一步搜索到的最優(yōu)解對目標函數(shù)沒有顯著變化,繼續(xù)求解所花費的時間和資源對優(yōu)化的目標沒有顯著意義,目前的求解結(jié)果已經(jīng)可以滿足實際規(guī)劃方案了。

    ② 當種群進化到某一代數(shù)時終止循環(huán)。第二種方法相對第一種更便于理解和操作,可以根據(jù)結(jié)果靈活調(diào)整。因此選擇種群進化最大代數(shù)作為種群進化終止的判定條件。

    3 實例分析

    3.1 實驗算例

    通過對某區(qū)域的公交客流點進行支付意愿調(diào)查,該區(qū)域內(nèi)有5個充電站,該區(qū)域的公交客流需求點、起終點、充電樁如圖6所示。

    圖6 乘客坐標Fig.6 Passenger coordinates

    3.2 求解結(jié)果

    首先根據(jù)坐標信息使用K-means算法將乘客分為32類確定公交??空军c,如圖7所示。然后根據(jù)乘客支付意愿確定每個站點的上車人數(shù),公交??奎c信息如表1所列。

    圖7 公交??奎c規(guī)劃Fig.7 Bus stop planning

    表1 公交??奎c信息Table 1 Bus stop information table

    然后基于ESGA規(guī)劃純電動需求響應公交路徑,其中種群規(guī)模為500,進化代數(shù)為1 000,精英種群比例為15%,變異概率為0.1,反轉(zhuǎn)概率為0.1,站場車輛數(shù)為14,每輛車的限載人數(shù)為30,單位時間懲罰成本為0.5元/min,車輛啟動成本司機工資為150元/輛,可變成本為0.49元/km。經(jīng)過20次的獨立重復試驗,選取最好試驗結(jié)果作為最終結(jié)果,其進化曲線、不同情況下的車輛路徑和計算結(jié)果見圖8~圖11和表2。

    圖8 進化曲線Fig.8 Evolution curve

    圖9 不需充電的車輛行駛路徑Fig.9 The driving path of the vehicle without charging

    圖10 需要充電的車輛行駛路徑Fig.10 The driving path of the vehicle that needs to be charged

    圖11 補充電能車輛行駛路徑Fig.11 Supplementary electric vehicle driving path

    表2 計算結(jié)果Table 2 Calculation result table

    由表2可以看出每輛公交車的實際行駛里程均在30 km左右,約為單次充電行駛里程(60 km)的一半。需求響應公交載客率均在50%以上,平均載客率為86.67%,其中11輛車的載客率在80%以上,企業(yè)的成本包括人力成本、可變成本和時間懲罰成本3部分,分別為2 100.00元、226.35元和361.97元,企業(yè)的總收入全部來源于乘客乘坐響應公交所支付的票價,共7 643元。企業(yè)的總利潤是企業(yè)的總收入減去總成本,為4 954.68元。

    4 結(jié)語

    通過描述純電動需求響應式公交路徑優(yōu)化問題,構(gòu)建了特定假設(shè)下的數(shù)學模型,采用ESGA算法進行問題求解,并通過實際算例驗證了ESGA算法,結(jié)果表明ESGA算法精英種群的規(guī)模設(shè)為15%左右時求解結(jié)果較好。由于研究水平有限,在對純電動需求響應公交進行路徑規(guī)劃時有部分問題沒有考慮到,存在一定的局限性和不足。今后主要從以下幾個方面進行研究和完善:

    (1) 此次研究只考慮了單一的純電動車型,而企業(yè)在實際運營過程中,可以根據(jù)乘客的個性化需求和實際乘客數(shù)量來考慮同時使用多種車型,這樣可以節(jié)約運營成本,提高需求響應公交滿意度。

    (2) 此次研究未考慮動態(tài)響應需求,在實際運營中會出現(xiàn)新的需求,考慮實時的動態(tài)需求不僅可以增加企業(yè)的收入,也可以為更多乘客提供需求響應服務。

    猜你喜歡
    ??空?/a>公交車公交
    你們認識嗎
    一元公交開進太行深處
    寒假像一列火車
    等公交
    等公交
    公交車上
    公交車奇妙日
    幼兒畫刊(2017年5期)2017-06-21 21:17:02
    公交??空驹O(shè)置形式對路段交通運行狀態(tài)影響分析
    城里的公交車
    小布老虎(2016年12期)2016-12-01 05:46:57
    城市公交??空驹O(shè)置問題的探討
    福建建筑(2015年2期)2015-02-18 02:32:22
    久久国产精品大桥未久av| 欧美日韩福利视频一区二区| 欧美av亚洲av综合av国产av| 精品少妇一区二区三区视频日本电影| 黑人猛操日本美女一级片| 少妇的丰满在线观看| 亚洲五月婷婷丁香| 老鸭窝网址在线观看| 亚洲av国产av综合av卡| 国产精品自产拍在线观看55亚洲 | 国精品久久久久久国模美| 中文精品一卡2卡3卡4更新| 69av精品久久久久久 | 亚洲一码二码三码区别大吗| av线在线观看网站| 黄网站色视频无遮挡免费观看| 在线观看免费午夜福利视频| 国产精品九九99| 国产日韩欧美视频二区| 国产一区二区 视频在线| 性少妇av在线| tocl精华| 在线看a的网站| 99精品欧美一区二区三区四区| 精品亚洲成国产av| 美女国产高潮福利片在线看| 精品久久蜜臀av无| 最黄视频免费看| 人人妻人人澡人人看| 日韩一卡2卡3卡4卡2021年| 色婷婷av一区二区三区视频| 男女边摸边吃奶| 最近中文字幕2019免费版| 日韩人妻精品一区2区三区| 捣出白浆h1v1| 亚洲第一青青草原| 国产精品二区激情视频| 99国产精品免费福利视频| 美女中出高潮动态图| 日韩中文字幕欧美一区二区| 国产免费一区二区三区四区乱码| 操美女的视频在线观看| 最近最新中文字幕大全免费视频| 少妇被粗大的猛进出69影院| 韩国精品一区二区三区| 亚洲午夜精品一区,二区,三区| 久久ye,这里只有精品| 美女大奶头黄色视频| 久久亚洲国产成人精品v| 国产又爽黄色视频| 男人舔女人的私密视频| 777久久人妻少妇嫩草av网站| 9色porny在线观看| 国产日韩一区二区三区精品不卡| 国产麻豆69| 国产成人欧美在线观看 | 亚洲欧美精品综合一区二区三区| 丰满少妇做爰视频| 曰老女人黄片| 老熟妇乱子伦视频在线观看 | 亚洲欧洲日产国产| 狠狠婷婷综合久久久久久88av| 精品少妇久久久久久888优播| 在线av久久热| tocl精华| 在线av久久热| 国产高清videossex| 正在播放国产对白刺激| 亚洲国产精品一区三区| svipshipincom国产片| 久久综合国产亚洲精品| 亚洲色图综合在线观看| 99久久综合免费| 欧美激情久久久久久爽电影 | 亚洲欧美精品自产自拍| 纵有疾风起免费观看全集完整版| 亚洲av国产av综合av卡| 午夜福利在线观看吧| 老司机福利观看| 制服人妻中文乱码| 美女福利国产在线| 老熟妇乱子伦视频在线观看 | 人人妻人人爽人人添夜夜欢视频| 免费少妇av软件| 久久天堂一区二区三区四区| 国产成人系列免费观看| 国产主播在线观看一区二区| 黄色视频,在线免费观看| 一级,二级,三级黄色视频| 久久天躁狠狠躁夜夜2o2o| 久久性视频一级片| 免费久久久久久久精品成人欧美视频| 国产精品 欧美亚洲| 一区二区三区四区激情视频| 免费在线观看日本一区| 国产精品免费视频内射| 欧美av亚洲av综合av国产av| 国产精品二区激情视频| 人成视频在线观看免费观看| 狠狠狠狠99中文字幕| 一区二区三区精品91| 国产成人免费观看mmmm| 亚洲熟女精品中文字幕| 中文字幕色久视频| 一级毛片女人18水好多| 9热在线视频观看99| 建设人人有责人人尽责人人享有的| 美女国产高潮福利片在线看| 一区二区三区乱码不卡18| 国产精品自产拍在线观看55亚洲 | 91老司机精品| 高清欧美精品videossex| av欧美777| 女性被躁到高潮视频| 国产不卡av网站在线观看| 午夜免费观看性视频| 视频区欧美日本亚洲| 亚洲一区二区三区欧美精品| 男人操女人黄网站| 国产欧美日韩一区二区精品| 国产一区二区三区在线臀色熟女 | 久久精品成人免费网站| 国产男人的电影天堂91| 王馨瑶露胸无遮挡在线观看| 亚洲中文字幕日韩| 午夜日韩欧美国产| 丝瓜视频免费看黄片| 18禁国产床啪视频网站| 肉色欧美久久久久久久蜜桃| e午夜精品久久久久久久| 夜夜骑夜夜射夜夜干| 一级a爱视频在线免费观看| 日韩视频在线欧美| 精品人妻熟女毛片av久久网站| 黄色视频不卡| 香蕉国产在线看| 男女高潮啪啪啪动态图| 热re99久久国产66热| 不卡av一区二区三区| 天天躁夜夜躁狠狠躁躁| 国产91精品成人一区二区三区 | 一个人免费在线观看的高清视频 | 欧美黄色片欧美黄色片| 十八禁网站免费在线| 永久免费av网站大全| 日韩大码丰满熟妇| 久久久精品国产亚洲av高清涩受| 美女扒开内裤让男人捅视频| 国产成+人综合+亚洲专区| 女人被躁到高潮嗷嗷叫费观| 另类亚洲欧美激情| 十八禁网站网址无遮挡| 久久99一区二区三区| 国产成人影院久久av| 日韩一卡2卡3卡4卡2021年| 国产人伦9x9x在线观看| 欧美乱码精品一区二区三区| 每晚都被弄得嗷嗷叫到高潮| www.999成人在线观看| 日韩人妻精品一区2区三区| 久久影院123| 国产1区2区3区精品| 999久久久国产精品视频| 丝袜美腿诱惑在线| 精品一区二区三区av网在线观看 | 少妇精品久久久久久久| 1024视频免费在线观看| 国产成人免费无遮挡视频| 三上悠亚av全集在线观看| 精品国产一区二区三区久久久樱花| 99香蕉大伊视频| 啦啦啦在线免费观看视频4| 国产亚洲av片在线观看秒播厂| 国产成人精品无人区| 高清在线国产一区| 亚洲成人免费av在线播放| 最近最新免费中文字幕在线| 黄色 视频免费看| 十八禁网站免费在线| 成人国语在线视频| 激情视频va一区二区三区| 日本黄色日本黄色录像| 亚洲国产欧美在线一区| 欧美 日韩 精品 国产| 久久国产精品大桥未久av| 国产亚洲欧美精品永久| 国产欧美日韩一区二区三区在线| 肉色欧美久久久久久久蜜桃| 淫妇啪啪啪对白视频 | 欧美日韩成人在线一区二区| 脱女人内裤的视频| 日韩电影二区| 热99re8久久精品国产| 老司机深夜福利视频在线观看 | 黑人欧美特级aaaaaa片| 少妇精品久久久久久久| 久久精品成人免费网站| 国产成人精品久久二区二区免费| 亚洲精品久久成人aⅴ小说| 后天国语完整版免费观看| 黑人巨大精品欧美一区二区mp4| 91麻豆精品激情在线观看国产 | 天堂俺去俺来也www色官网| 大香蕉久久网| 夫妻午夜视频| 亚洲av成人不卡在线观看播放网 | 97人妻天天添夜夜摸| 欧美黄色片欧美黄色片| 欧美日韩亚洲国产一区二区在线观看 | 日本av免费视频播放| 婷婷色av中文字幕| 国产日韩欧美视频二区| 精品亚洲乱码少妇综合久久| 免费在线观看完整版高清| 天堂中文最新版在线下载| 最黄视频免费看| 国产高清videossex| 下体分泌物呈黄色| 国产亚洲av片在线观看秒播厂| 热99国产精品久久久久久7| 男人舔女人的私密视频| 欧美黄色片欧美黄色片| 又大又爽又粗| 久久久久精品国产欧美久久久 | 日韩一区二区三区影片| 蜜桃国产av成人99| 久久久精品94久久精品| 国产精品二区激情视频| 免费观看av网站的网址| 制服诱惑二区| 久久久久久久久免费视频了| 91字幕亚洲| 成人亚洲精品一区在线观看| 日本一区二区免费在线视频| 精品高清国产在线一区| 999久久久国产精品视频| 精品国产乱码久久久久久男人| 欧美亚洲 丝袜 人妻 在线| 欧美黄色淫秽网站| 亚洲va日本ⅴa欧美va伊人久久 | 黄色视频,在线免费观看| 欧美日韩一级在线毛片| 无限看片的www在线观看| 欧美av亚洲av综合av国产av| 1024视频免费在线观看| 老司机午夜福利在线观看视频 | 亚洲人成电影免费在线| 久久ye,这里只有精品| 女性被躁到高潮视频| 国产在线一区二区三区精| 国产三级黄色录像| 亚洲人成电影免费在线| 精品国内亚洲2022精品成人 | 51午夜福利影视在线观看| 爱豆传媒免费全集在线观看| 日本av手机在线免费观看| 一边摸一边做爽爽视频免费| 热99久久久久精品小说推荐| 视频在线观看一区二区三区| √禁漫天堂资源中文www| 亚洲色图综合在线观看| 两人在一起打扑克的视频| 免费在线观看日本一区| 老司机靠b影院| 19禁男女啪啪无遮挡网站| 后天国语完整版免费观看| 免费在线观看日本一区| 蜜桃国产av成人99| 国产主播在线观看一区二区| 亚洲一区二区三区欧美精品| 人妻久久中文字幕网| 亚洲va日本ⅴa欧美va伊人久久 | 日韩人妻精品一区2区三区| 亚洲伊人久久精品综合| 一区二区三区乱码不卡18| 国产又色又爽无遮挡免| 亚洲精品乱久久久久久| 欧美+亚洲+日韩+国产| 嫁个100分男人电影在线观看| 女人爽到高潮嗷嗷叫在线视频| 亚洲av美国av| 伦理电影免费视频| 男女无遮挡免费网站观看| bbb黄色大片| 黄网站色视频无遮挡免费观看| 精品久久久久久久毛片微露脸 | 国产日韩欧美在线精品| 菩萨蛮人人尽说江南好唐韦庄| 99热全是精品| 国产黄色免费在线视频| 美女国产高潮福利片在线看| 精品国产乱码久久久久久男人| 高潮久久久久久久久久久不卡| 青春草视频在线免费观看| 亚洲av男天堂| 亚洲 国产 在线| 国产片内射在线| 在线十欧美十亚洲十日本专区| 十八禁网站网址无遮挡| 青春草亚洲视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 99久久精品国产亚洲精品| 女人被躁到高潮嗷嗷叫费观| 欧美日韩成人在线一区二区| 日韩 亚洲 欧美在线| 国内毛片毛片毛片毛片毛片| 在线观看舔阴道视频| 国产欧美亚洲国产| 人人妻人人添人人爽欧美一区卜| 精品国产一区二区三区四区第35| 一进一出抽搐动态| 每晚都被弄得嗷嗷叫到高潮| 美女主播在线视频| 亚洲成人免费电影在线观看| 汤姆久久久久久久影院中文字幕| 老司机在亚洲福利影院| 国产成人精品久久二区二区91| 在线观看免费高清a一片| 男人添女人高潮全过程视频| 少妇的丰满在线观看| 18禁观看日本| 免费在线观看日本一区| 欧美变态另类bdsm刘玥| 久久久国产一区二区| 99久久人妻综合| 亚洲黑人精品在线| 老司机深夜福利视频在线观看 | 亚洲专区字幕在线| 免费观看人在逋| 午夜91福利影院| 99国产精品免费福利视频| 久久久精品94久久精品| 欧美午夜高清在线| 伦理电影免费视频| 一区福利在线观看| 欧美日韩亚洲国产一区二区在线观看 | 久9热在线精品视频| 激情视频va一区二区三区| 母亲3免费完整高清在线观看| 国产成+人综合+亚洲专区| 99精国产麻豆久久婷婷| 亚洲av成人一区二区三| 欧美日韩亚洲综合一区二区三区_| 国产人伦9x9x在线观看| 亚洲国产精品999| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品日韩在线中文字幕| 丝袜脚勾引网站| 欧美激情高清一区二区三区| 免费在线观看视频国产中文字幕亚洲 | 多毛熟女@视频| 91成人精品电影| 人人妻人人澡人人看| 欧美亚洲日本最大视频资源| 日本av免费视频播放| 秋霞在线观看毛片| 欧美日韩成人在线一区二区| 欧美成人午夜精品| 亚洲精品自拍成人| 午夜激情av网站| 国产成人精品久久二区二区91| 搡老熟女国产l中国老女人| 婷婷丁香在线五月| 日本五十路高清| 色婷婷av一区二区三区视频| 99久久人妻综合| 亚洲人成电影免费在线| av线在线观看网站| 精品一区二区三区四区五区乱码| 丝瓜视频免费看黄片| 99香蕉大伊视频| 不卡一级毛片| 午夜精品久久久久久毛片777| 丁香六月欧美| 久久精品人人爽人人爽视色| 人人澡人人妻人| 在线永久观看黄色视频| 国产一区二区 视频在线| 亚洲男人天堂网一区| 天堂8中文在线网| 在线永久观看黄色视频| 久久午夜综合久久蜜桃| 1024香蕉在线观看| xxxhd国产人妻xxx| 岛国毛片在线播放| 在线观看免费高清a一片| 在线天堂中文资源库| 久久人人爽av亚洲精品天堂| 热99久久久久精品小说推荐| 伦理电影免费视频| 午夜福利免费观看在线| 久久国产精品人妻蜜桃| 精品人妻熟女毛片av久久网站| 国产激情久久老熟女| 欧美在线黄色| 午夜成年电影在线免费观看| 欧美黄色片欧美黄色片| 啦啦啦在线免费观看视频4| 性高湖久久久久久久久免费观看| 窝窝影院91人妻| 亚洲中文av在线| 国产1区2区3区精品| 大香蕉久久网| 2018国产大陆天天弄谢| 91精品国产国语对白视频| 国产精品二区激情视频| 大片电影免费在线观看免费| 亚洲精品中文字幕在线视频| 美女视频免费永久观看网站| 精品人妻1区二区| 国产黄色免费在线视频| 亚洲av片天天在线观看| 亚洲,欧美精品.| 国产97色在线日韩免费| 每晚都被弄得嗷嗷叫到高潮| 麻豆av在线久日| 一区二区三区乱码不卡18| 人人妻人人爽人人添夜夜欢视频| 老司机午夜十八禁免费视频| 超碰97精品在线观看| 日本a在线网址| 亚洲欧美一区二区三区久久| 91精品三级在线观看| 国产亚洲午夜精品一区二区久久| 好男人电影高清在线观看| 亚洲人成电影观看| 日韩 亚洲 欧美在线| 成年美女黄网站色视频大全免费| 在线观看免费高清a一片| 久久久国产精品麻豆| 午夜日韩欧美国产| 久久99一区二区三区| 美女福利国产在线| 丝袜美腿诱惑在线| 97人妻天天添夜夜摸| 中亚洲国语对白在线视频| 国产亚洲一区二区精品| 国产欧美日韩一区二区精品| 可以免费在线观看a视频的电影网站| 亚洲少妇的诱惑av| 一区二区三区四区激情视频| 国产亚洲欧美在线一区二区| 无限看片的www在线观看| 纵有疾风起免费观看全集完整版| 又紧又爽又黄一区二区| 中文字幕人妻熟女乱码| 欧美+亚洲+日韩+国产| 国产成人免费观看mmmm| 国产在线观看jvid| 高清在线国产一区| 亚洲成av片中文字幕在线观看| 青青草视频在线视频观看| 亚洲精品美女久久av网站| 国产成人免费无遮挡视频| 精品久久久久久久毛片微露脸 | 精品一区二区三卡| 黑人猛操日本美女一级片| 宅男免费午夜| 日日夜夜操网爽| 丝袜美腿诱惑在线| 又紧又爽又黄一区二区| 国产极品粉嫩免费观看在线| 99国产极品粉嫩在线观看| 午夜激情久久久久久久| 国产精品国产av在线观看| 欧美精品av麻豆av| 国产欧美日韩精品亚洲av| 亚洲国产欧美网| 亚洲国产精品一区二区三区在线| 久久久国产一区二区| 午夜精品国产一区二区电影| 日本vs欧美在线观看视频| 亚洲欧美色中文字幕在线| 亚洲欧洲日产国产| 午夜免费观看性视频| 大码成人一级视频| 99久久人妻综合| 国产一级毛片在线| 久久精品国产亚洲av高清一级| 中文精品一卡2卡3卡4更新| 动漫黄色视频在线观看| 亚洲国产欧美一区二区综合| 欧美日韩亚洲国产一区二区在线观看 | 两人在一起打扑克的视频| 亚洲精品国产av蜜桃| 免费女性裸体啪啪无遮挡网站| 中文字幕色久视频| 韩国高清视频一区二区三区| 大型av网站在线播放| 99国产精品一区二区三区| 性色av乱码一区二区三区2| 王馨瑶露胸无遮挡在线观看| 国产欧美亚洲国产| 久久综合国产亚洲精品| 午夜福利,免费看| 黑人猛操日本美女一级片| 青春草视频在线免费观看| 欧美日韩亚洲国产一区二区在线观看 | 老汉色av国产亚洲站长工具| 男女国产视频网站| 亚洲精品久久成人aⅴ小说| 欧美一级毛片孕妇| 日韩精品免费视频一区二区三区| 美女扒开内裤让男人捅视频| 成年人免费黄色播放视频| netflix在线观看网站| a级毛片在线看网站| 国产精品久久久久久精品古装| 久久这里只有精品19| 捣出白浆h1v1| 搡老乐熟女国产| 十八禁网站免费在线| 99九九在线精品视频| 亚洲黑人精品在线| 亚洲精品国产av蜜桃| 桃花免费在线播放| 少妇被粗大的猛进出69影院| 亚洲成人手机| 交换朋友夫妻互换小说| 999久久久精品免费观看国产| 国产精品秋霞免费鲁丝片| 搡老岳熟女国产| 我的亚洲天堂| 亚洲欧美一区二区三区黑人| 亚洲精品中文字幕在线视频| 亚洲免费av在线视频| 蜜桃在线观看..| 91精品国产国语对白视频| 日本av手机在线免费观看| 亚洲精品粉嫩美女一区| 免费日韩欧美在线观看| 一区二区三区精品91| 亚洲av成人不卡在线观看播放网 | 1024香蕉在线观看| 搡老熟女国产l中国老女人| 日韩大码丰满熟妇| videos熟女内射| www.熟女人妻精品国产| 国产男女超爽视频在线观看| 欧美日韩av久久| 一区在线观看完整版| 精品亚洲成a人片在线观看| 国产亚洲午夜精品一区二区久久| 免费在线观看完整版高清| 亚洲av成人一区二区三| 丝袜人妻中文字幕| 乱人伦中国视频| 国产激情久久老熟女| 日韩,欧美,国产一区二区三区| 我要看黄色一级片免费的| 日韩大码丰满熟妇| 久久久国产欧美日韩av| 国产高清国产精品国产三级| 久久久国产精品麻豆| a 毛片基地| 免费观看a级毛片全部| 亚洲欧美色中文字幕在线| 成人18禁高潮啪啪吃奶动态图| 国产亚洲欧美在线一区二区| 一级黄色大片毛片| 国产日韩欧美视频二区| 日韩制服丝袜自拍偷拍| 国产亚洲av高清不卡| 99国产精品一区二区蜜桃av | 人人妻人人澡人人爽人人夜夜| 中文字幕av电影在线播放| 岛国毛片在线播放| 日日摸夜夜添夜夜添小说| 久久国产精品大桥未久av| 欧美乱码精品一区二区三区| 亚洲精品中文字幕一二三四区 | 久久人人爽人人片av| 国产国语露脸激情在线看| 成人av一区二区三区在线看 | 免费在线观看影片大全网站| 精品国产一区二区久久| 精品亚洲乱码少妇综合久久| e午夜精品久久久久久久| 精品人妻在线不人妻| 国产精品一区二区在线不卡| 亚洲成人国产一区在线观看| 精品久久久精品久久久| 超色免费av| 日本五十路高清| 黄色a级毛片大全视频| 国产成人精品久久二区二区91| 成人av一区二区三区在线看 | 精品少妇黑人巨大在线播放| av不卡在线播放| 一区在线观看完整版| 日韩视频在线欧美| 窝窝影院91人妻| 亚洲精品日韩在线中文字幕| 午夜免费成人在线视频| 午夜老司机福利片| 久久精品国产亚洲av香蕉五月 | 日韩熟女老妇一区二区性免费视频| 中文字幕高清在线视频| 大香蕉久久成人网| 丁香六月天网| 黑人欧美特级aaaaaa片| 51午夜福利影视在线观看| 亚洲全国av大片| 欧美国产精品一级二级三级| 国产精品秋霞免费鲁丝片| 久久精品国产亚洲av香蕉五月 | 中国国产av一级| 男女免费视频国产| 亚洲精品国产av蜜桃| 免费av中文字幕在线| 黑丝袜美女国产一区| 午夜两性在线视频| 人妻久久中文字幕网| 老司机影院毛片| 黄网站色视频无遮挡免费观看|