• <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
    高清在线国产一区| tocl精华| 青草久久国产| 操美女的视频在线观看| 午夜免费鲁丝| 成年人午夜在线观看视频| 亚洲va日本ⅴa欧美va伊人久久 | 亚洲成人免费av在线播放| 一本一本久久a久久精品综合妖精| 精品久久久久久久毛片微露脸 | 人人妻人人添人人爽欧美一区卜| 日韩 亚洲 欧美在线| 久久国产精品男人的天堂亚洲| 亚洲av电影在线进入| 久9热在线精品视频| 亚洲精品国产色婷婷电影| 色老头精品视频在线观看| 国产成人系列免费观看| 日韩欧美免费精品| 夜夜夜夜夜久久久久| 热re99久久国产66热| 12—13女人毛片做爰片一| 好男人电影高清在线观看| 久久亚洲国产成人精品v| 久久久国产成人免费| 青春草亚洲视频在线观看| 夜夜骑夜夜射夜夜干| 自线自在国产av| 亚洲精品国产精品久久久不卡| 国产欧美日韩一区二区精品| 亚洲色图综合在线观看| 不卡一级毛片| 国产欧美亚洲国产| 精品一区二区三区四区五区乱码| 国产成+人综合+亚洲专区| av线在线观看网站| 免费观看人在逋| 女人精品久久久久毛片| 在线 av 中文字幕| avwww免费| 天天操日日干夜夜撸| 精品欧美一区二区三区在线| 精品福利观看| 宅男免费午夜| 国产av又大| 久久热在线av| 亚洲精品中文字幕在线视频| videosex国产| 性少妇av在线| 午夜日韩欧美国产| tube8黄色片| 国内毛片毛片毛片毛片毛片| 欧美另类亚洲清纯唯美| 国产亚洲av高清不卡| 久久精品熟女亚洲av麻豆精品| 狠狠狠狠99中文字幕| 久久久久国产一级毛片高清牌| 国产成+人综合+亚洲专区| 国产亚洲av片在线观看秒播厂| 大香蕉久久网| 首页视频小说图片口味搜索| 水蜜桃什么品种好| 国产又色又爽无遮挡免| 日本av免费视频播放| 久久久久久久久免费视频了| 亚洲国产欧美日韩在线播放| 精品久久久久久久毛片微露脸 | 亚洲精品中文字幕在线视频| 日韩视频一区二区在线观看| 亚洲第一欧美日韩一区二区三区 | 水蜜桃什么品种好| 免费在线观看完整版高清| 久久性视频一级片| 少妇猛男粗大的猛烈进出视频| 成人国产一区最新在线观看| 亚洲欧美精品自产自拍| www.精华液| 一级a爱视频在线免费观看| 国产主播在线观看一区二区| 久久性视频一级片| 丁香六月欧美| 男女午夜视频在线观看| 欧美精品一区二区大全| 9色porny在线观看| 久久久久久久大尺度免费视频| 欧美日本中文国产一区发布| 两人在一起打扑克的视频| 午夜福利影视在线免费观看| 日韩视频在线欧美| 欧美亚洲 丝袜 人妻 在线| 久久人人爽av亚洲精品天堂| 欧美亚洲日本最大视频资源| 欧美日韩中文字幕国产精品一区二区三区 | 国产av精品麻豆| 在线av久久热| 97人妻天天添夜夜摸| 考比视频在线观看| 欧美日韩一级在线毛片| 黑人操中国人逼视频| 新久久久久国产一级毛片| 51午夜福利影视在线观看| 高清av免费在线| 午夜激情av网站| 国产免费福利视频在线观看| 国产一区二区 视频在线| 欧美亚洲日本最大视频资源| 黑丝袜美女国产一区| 亚洲成国产人片在线观看| 一边摸一边抽搐一进一出视频| 午夜福利乱码中文字幕| 别揉我奶头~嗯~啊~动态视频 | 永久免费av网站大全| 国产成人啪精品午夜网站| 久久久久国产精品人妻一区二区| 国产在线免费精品| 国产激情久久老熟女| 又紧又爽又黄一区二区| 国产亚洲欧美在线一区二区| 中文精品一卡2卡3卡4更新| 老司机在亚洲福利影院| 99精品欧美一区二区三区四区| 午夜精品久久久久久毛片777| 91精品伊人久久大香线蕉| av在线app专区| 狠狠狠狠99中文字幕| 亚洲精品国产色婷婷电影| 999精品在线视频| 国产精品一区二区在线观看99| 桃花免费在线播放| 亚洲中文字幕日韩| videosex国产| av一本久久久久| 欧美日韩黄片免| 成人免费观看视频高清| 夫妻午夜视频| 性色av一级| 十八禁网站网址无遮挡| 日韩中文字幕视频在线看片| 黑人巨大精品欧美一区二区mp4| 久久久水蜜桃国产精品网| 精品国产一区二区久久| 精品福利观看| 在线精品无人区一区二区三| 精品久久蜜臀av无| 亚洲精品中文字幕在线视频| 国产成人精品无人区| 伦理电影免费视频| 亚洲少妇的诱惑av| 一区二区三区精品91| 美女福利国产在线| 久久毛片免费看一区二区三区| 麻豆国产av国片精品| 亚洲欧美一区二区三区久久| 国产一卡二卡三卡精品| 各种免费的搞黄视频| 黑人巨大精品欧美一区二区mp4| 国产欧美亚洲国产| 成年女人毛片免费观看观看9 | 伦理电影免费视频| 涩涩av久久男人的天堂| 大香蕉久久成人网| 欧美 亚洲 国产 日韩一| 欧美变态另类bdsm刘玥| 三级毛片av免费| 五月天丁香电影| 亚洲成人国产一区在线观看| 国产精品亚洲av一区麻豆| 777久久人妻少妇嫩草av网站| 深夜精品福利| 国产区一区二久久| 欧美精品亚洲一区二区| 国产男女超爽视频在线观看| 制服人妻中文乱码| 天天躁狠狠躁夜夜躁狠狠躁| 青草久久国产| 欧美性长视频在线观看| 国产一区二区激情短视频 | 欧美激情 高清一区二区三区| 热99re8久久精品国产| 久久人妻福利社区极品人妻图片| 狂野欧美激情性xxxx| 欧美成人午夜精品| 免费高清在线观看日韩| 国产日韩一区二区三区精品不卡| 欧美+亚洲+日韩+国产| 亚洲一码二码三码区别大吗| 国产1区2区3区精品| 天天影视国产精品| 久久久精品国产亚洲av高清涩受| 国产亚洲一区二区精品| 亚洲九九香蕉| 天堂俺去俺来也www色官网| 精品久久久久久久毛片微露脸 | 丝袜喷水一区| 中文字幕人妻熟女乱码| 亚洲专区国产一区二区| 国产精品av久久久久免费| 久久久精品94久久精品| 精品亚洲成a人片在线观看| 久久免费观看电影| 两性午夜刺激爽爽歪歪视频在线观看 | 两人在一起打扑克的视频| 亚洲欧美日韩另类电影网站| 国产亚洲精品久久久久5区| 亚洲国产精品一区三区| videos熟女内射| 亚洲精品日韩在线中文字幕| 美国免费a级毛片| 1024香蕉在线观看| 久久久国产一区二区| 老司机靠b影院| cao死你这个sao货| 亚洲第一电影网av| 亚洲av成人不卡在线观看播放网| 丝袜人妻中文字幕| 欧美高清成人免费视频www| 男女做爰动态图高潮gif福利片| 亚洲人成电影免费在线| 夜夜看夜夜爽夜夜摸| 叶爱在线成人免费视频播放| 国产精品电影一区二区三区| 免费电影在线观看免费观看| 色综合欧美亚洲国产小说| 身体一侧抽搐| 国模一区二区三区四区视频 | 亚洲乱码一区二区免费版| 欧美成人一区二区免费高清观看 | 国产男靠女视频免费网站| 国模一区二区三区四区视频 | 日韩欧美精品v在线| 亚洲成av人片免费观看| av福利片在线观看| 亚洲av熟女| 啦啦啦观看免费观看视频高清| 中文字幕人成人乱码亚洲影| av片东京热男人的天堂| 99久久久亚洲精品蜜臀av| av视频在线观看入口| 欧美最黄视频在线播放免费| 免费电影在线观看免费观看| 国产精品 欧美亚洲| 国产精品影院久久| 天天添夜夜摸| 91老司机精品| 久久精品国产清高在天天线| 久久草成人影院| 欧美一区二区精品小视频在线| 99久久无色码亚洲精品果冻| 亚洲中文日韩欧美视频| 久久精品国产清高在天天线| 国产成人啪精品午夜网站| 国产精品影院久久| 免费在线观看黄色视频的| 亚洲欧美精品综合一区二区三区| 性色av乱码一区二区三区2| 精品国产乱子伦一区二区三区| www日本黄色视频网| 99精品久久久久人妻精品| 国产精品98久久久久久宅男小说| 国产成人精品久久二区二区91| 丁香六月欧美| АⅤ资源中文在线天堂| 久久热在线av| 一进一出好大好爽视频| 国产久久久一区二区三区| 人人妻,人人澡人人爽秒播| 久久久久九九精品影院| 久热爱精品视频在线9| 舔av片在线| 中文字幕久久专区| 欧美午夜高清在线| 美女扒开内裤让男人捅视频| 亚洲成人精品中文字幕电影| 老司机午夜十八禁免费视频| 久久久久久九九精品二区国产 | 97碰自拍视频| 亚洲一卡2卡3卡4卡5卡精品中文| 黄色成人免费大全| 老鸭窝网址在线观看| 在线观看舔阴道视频| 十八禁网站免费在线| 成人18禁在线播放| 亚洲欧美日韩无卡精品| 亚洲一卡2卡3卡4卡5卡精品中文| 999久久久精品免费观看国产| av超薄肉色丝袜交足视频| 久久久久免费精品人妻一区二区| 国产av不卡久久| 高清在线国产一区| 久久精品夜夜夜夜夜久久蜜豆 | 久久久久久人人人人人| 久久午夜亚洲精品久久| 国语自产精品视频在线第100页| 丝袜美腿诱惑在线| 亚洲午夜理论影院| or卡值多少钱| 亚洲国产欧美网| www.999成人在线观看| 国产成人啪精品午夜网站| 亚洲欧美日韩无卡精品| 女人被狂操c到高潮| 日韩免费av在线播放| 欧美成狂野欧美在线观看| av福利片在线观看| 亚洲精华国产精华精| 日韩国内少妇激情av| 久久久久久国产a免费观看| 最新美女视频免费是黄的| 国内揄拍国产精品人妻在线| 91麻豆av在线| 99re在线观看精品视频| 老司机福利观看| 亚洲成人精品中文字幕电影| 欧美丝袜亚洲另类 | 欧美日韩中文字幕国产精品一区二区三区| 国产激情欧美一区二区| 国产欧美日韩精品亚洲av| 制服人妻中文乱码| 精品乱码久久久久久99久播| 一本精品99久久精品77| 丰满人妻一区二区三区视频av | 丝袜人妻中文字幕| 欧美成人午夜精品| 黄片大片在线免费观看| 久久香蕉精品热| a级毛片在线看网站| 两性午夜刺激爽爽歪歪视频在线观看 | 九色国产91popny在线| 国产av一区二区精品久久| 亚洲av电影不卡..在线观看| 国产精品香港三级国产av潘金莲| 亚洲一区高清亚洲精品| 天天躁狠狠躁夜夜躁狠狠躁| 国产v大片淫在线免费观看| 亚洲美女黄片视频| 可以免费在线观看a视频的电影网站| av天堂在线播放| 非洲黑人性xxxx精品又粗又长| 国产精品综合久久久久久久免费| 亚洲人成电影免费在线| 91字幕亚洲| 麻豆成人午夜福利视频| 少妇粗大呻吟视频| 一进一出抽搐动态| 亚洲av成人不卡在线观看播放网| aaaaa片日本免费| 天天躁夜夜躁狠狠躁躁| 狂野欧美白嫩少妇大欣赏| 亚洲成人久久性| 后天国语完整版免费观看| 亚洲中文字幕日韩| 精品欧美国产一区二区三| 日韩中文字幕欧美一区二区| 国产精华一区二区三区| 亚洲精品在线美女| 欧美人与性动交α欧美精品济南到| 精品一区二区三区视频在线观看免费| 制服诱惑二区| 日韩免费av在线播放| 久99久视频精品免费| 亚洲一区二区三区色噜噜| 国产亚洲av嫩草精品影院| 亚洲国产高清在线一区二区三| 午夜免费激情av| 亚洲成人免费电影在线观看| 禁无遮挡网站| 999精品在线视频| 1024手机看黄色片| 好男人电影高清在线观看| 国产麻豆成人av免费视频| 亚洲专区中文字幕在线| 午夜激情福利司机影院| 久久久国产欧美日韩av| 一卡2卡三卡四卡精品乱码亚洲| 在线观看一区二区三区| 国内久久婷婷六月综合欲色啪| 久久婷婷人人爽人人干人人爱| 国产久久久一区二区三区| 国产av一区二区精品久久| 丰满人妻熟妇乱又伦精品不卡| 久久这里只有精品中国| 天堂影院成人在线观看| 亚洲av成人av| 色综合婷婷激情| 亚洲性夜色夜夜综合| 国产午夜精品论理片| 午夜亚洲福利在线播放| 精品福利观看| 男插女下体视频免费在线播放| 精品一区二区三区视频在线观看免费| 欧美日韩一级在线毛片| 99国产极品粉嫩在线观看| 亚洲欧美激情综合另类| 少妇人妻一区二区三区视频| 人妻久久中文字幕网| 亚洲精品久久成人aⅴ小说| 久久伊人香网站| 少妇被粗大的猛进出69影院| 精品第一国产精品| 亚洲精品粉嫩美女一区| 亚洲av片天天在线观看| 制服丝袜大香蕉在线| 99国产精品一区二区三区| 大型av网站在线播放| 国产精品一区二区三区四区免费观看 | 国产不卡一卡二| 午夜成年电影在线免费观看| 国产私拍福利视频在线观看| 精品国内亚洲2022精品成人| 欧美性猛交黑人性爽| 搡老熟女国产l中国老女人| 叶爱在线成人免费视频播放| 婷婷亚洲欧美| 桃红色精品国产亚洲av| 日本免费一区二区三区高清不卡| 国内精品久久久久久久电影| 搡老岳熟女国产| 法律面前人人平等表现在哪些方面| 亚洲激情在线av| 50天的宝宝边吃奶边哭怎么回事| 视频区欧美日本亚洲| 搡老熟女国产l中国老女人| 一夜夜www| 欧美精品亚洲一区二区| 香蕉国产在线看| 亚洲国产精品久久男人天堂| 宅男免费午夜| 国产激情偷乱视频一区二区| 制服丝袜大香蕉在线| 亚洲午夜精品一区,二区,三区| 男女视频在线观看网站免费 | 免费在线观看成人毛片| 美女午夜性视频免费| 国产精品永久免费网站| 久久精品影院6| 精品高清国产在线一区| 狠狠狠狠99中文字幕| 国产精品免费视频内射| 真人一进一出gif抽搐免费| ponron亚洲| xxxwww97欧美| 久久精品aⅴ一区二区三区四区| 国产av一区在线观看免费| 亚洲18禁久久av| 男男h啪啪无遮挡| 制服丝袜大香蕉在线| 岛国在线观看网站| 一级片免费观看大全| 精品国产亚洲在线| 免费在线观看亚洲国产| 亚洲国产看品久久| 最近最新中文字幕大全电影3| 国产精品 国内视频| 麻豆成人av在线观看| 久久九九热精品免费| 嫩草影视91久久| 女人被狂操c到高潮| 婷婷六月久久综合丁香| 精品一区二区三区av网在线观看| 国产精品自产拍在线观看55亚洲| 老司机在亚洲福利影院| 国产熟女午夜一区二区三区| 亚洲美女黄片视频| 久久久久精品国产欧美久久久| 国产精品一区二区三区四区久久| 在线国产一区二区在线| 免费人成视频x8x8入口观看| 三级国产精品欧美在线观看 | 特级一级黄色大片| 亚洲国产欧美一区二区综合| 国产高清视频在线播放一区| 可以免费在线观看a视频的电影网站| 国产免费av片在线观看野外av| 国内精品久久久久精免费| 日韩欧美三级三区| 亚洲人成网站高清观看| 麻豆国产97在线/欧美 | 久久99热这里只有精品18| 欧美日韩亚洲国产一区二区在线观看| 美女黄网站色视频| 亚洲国产欧美网| 日韩av在线大香蕉| 三级国产精品欧美在线观看 | 老司机深夜福利视频在线观看| 一级片免费观看大全| 老司机靠b影院| 性欧美人与动物交配| 草草在线视频免费看| 精品国内亚洲2022精品成人| 搡老妇女老女人老熟妇| 精品国产亚洲在线| 国产69精品久久久久777片 | 成年人黄色毛片网站| 精品人妻1区二区| 亚洲国产中文字幕在线视频| 久久久水蜜桃国产精品网| 一本大道久久a久久精品| www.www免费av| www.自偷自拍.com| 熟女少妇亚洲综合色aaa.| 亚洲精品美女久久久久99蜜臀| 国产午夜精品论理片| 亚洲精品久久国产高清桃花| 波多野结衣高清无吗| 熟女少妇亚洲综合色aaa.| 免费高清视频大片| 亚洲专区字幕在线| 99riav亚洲国产免费| 亚洲全国av大片| 精品午夜福利视频在线观看一区| 中文字幕人妻丝袜一区二区| 熟妇人妻久久中文字幕3abv| 给我免费播放毛片高清在线观看| 免费看a级黄色片| 国产人伦9x9x在线观看| 天天添夜夜摸| 亚洲人成网站在线播放欧美日韩| 日本撒尿小便嘘嘘汇集6| 国产高清激情床上av| 免费一级毛片在线播放高清视频| 免费搜索国产男女视频| 精品电影一区二区在线| 法律面前人人平等表现在哪些方面| 亚洲精品国产精品久久久不卡| 岛国在线观看网站| 欧美日韩一级在线毛片| 日韩精品免费视频一区二区三区| 日日夜夜操网爽| 成人午夜高清在线视频| 久久婷婷成人综合色麻豆| 亚洲第一欧美日韩一区二区三区| 国产欧美日韩精品亚洲av| 欧美精品啪啪一区二区三区| 国产真实乱freesex| 18禁黄网站禁片免费观看直播| 50天的宝宝边吃奶边哭怎么回事| 日本a在线网址| 国产午夜福利久久久久久| 国产单亲对白刺激| 国产欧美日韩精品亚洲av| 757午夜福利合集在线观看| 久久久国产成人免费| 日韩 欧美 亚洲 中文字幕| 高潮久久久久久久久久久不卡| 我要搜黄色片| 三级国产精品欧美在线观看 | 亚洲色图av天堂| 日本一区二区免费在线视频| 国产蜜桃级精品一区二区三区| 在线国产一区二区在线| 99riav亚洲国产免费| 村上凉子中文字幕在线| xxx96com| 日日干狠狠操夜夜爽| 国产高清视频在线播放一区| 日韩欧美在线二视频| 淫秽高清视频在线观看| 久久香蕉精品热| 亚洲一区高清亚洲精品| 天堂动漫精品| 手机成人av网站| 热99re8久久精品国产| 香蕉av资源在线| 国产成+人综合+亚洲专区| tocl精华| 美女午夜性视频免费| 免费一级毛片在线播放高清视频| 午夜福利在线观看吧| 国产三级在线视频| 高潮久久久久久久久久久不卡| 亚洲av熟女| 成人国产一区最新在线观看| 欧美大码av| 亚洲性夜色夜夜综合| 国产午夜福利久久久久久| 99久久综合精品五月天人人| 国产精品九九99| 亚洲专区国产一区二区| av片东京热男人的天堂| 天堂动漫精品| 免费无遮挡裸体视频| 亚洲欧美激情综合另类| 日本精品一区二区三区蜜桃| 曰老女人黄片| 久久久久性生活片| 免费在线观看完整版高清| 欧美精品啪啪一区二区三区| 狂野欧美白嫩少妇大欣赏| 亚洲成人国产一区在线观看| 天天一区二区日本电影三级| 日本一区二区免费在线视频| 亚洲成人国产一区在线观看| 亚洲欧美日韩东京热| 国产成人av教育| 久久天堂一区二区三区四区| 搞女人的毛片| 一级作爱视频免费观看| 一进一出好大好爽视频| 超碰成人久久| 在线观看免费午夜福利视频| 午夜免费观看网址| 日本熟妇午夜| 午夜激情福利司机影院| 欧美黑人欧美精品刺激| 国产精品国产高清国产av| 免费在线观看成人毛片| 天堂动漫精品| 亚洲精品一卡2卡三卡4卡5卡| 国产精品亚洲av一区麻豆| 麻豆国产97在线/欧美 | 久久久久久久精品吃奶| 人人妻,人人澡人人爽秒播| 精品国产乱子伦一区二区三区| a级毛片在线看网站| 少妇人妻一区二区三区视频|