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

    基于XGBoost模型的地鐵列車運(yùn)行狀態(tài)仿真方法

    2022-04-11 10:36:52門元昊劉曉雙秦梟喃羅森林
    城市軌道交通研究 2022年3期
    關(guān)鍵詞:列車運(yùn)行加速度列車

    門元昊 吳 亮 劉曉雙 秦梟喃 羅森林

    (1. 北京理工大學(xué)信息系統(tǒng)及安全對(duì)抗實(shí)驗(yàn)中心, 100081, 北京; 2. 通號(hào)城市軌道交通技術(shù)有限公司, 100081, 北京∥第一作者, 碩士研究生)

    列車運(yùn)行狀態(tài)仿真即通過建立真實(shí)運(yùn)行環(huán)境的近似模型,模擬地鐵列車在ATC(列車自動(dòng)控制)系統(tǒng)下的運(yùn)行表現(xiàn),可結(jié)合電子地圖實(shí)現(xiàn)列車位置、速度的估計(jì)反饋,是測(cè)試列車控制系統(tǒng)的重要工具。仿真的目的是為ATC系統(tǒng)的研究提供模擬試驗(yàn)平臺(tái)。由于ATC在設(shè)計(jì)、開發(fā)、測(cè)試等各個(gè)環(huán)節(jié)所需投入的時(shí)間長(zhǎng)且系統(tǒng)不斷更新迭代,在其研制過程中持續(xù)地在真實(shí)的生產(chǎn)環(huán)境中進(jìn)行開發(fā)與調(diào)試是不現(xiàn)實(shí)的,因此,通常需要先建立列車運(yùn)行狀態(tài)的仿真模型,在仿真模型上開展相關(guān)系統(tǒng)的研發(fā)工作。目前,針對(duì)列車運(yùn)行狀態(tài)的仿真仍以傳統(tǒng)的物理模型或改進(jìn)的物理模型為主[1-3],這些方法均對(duì)列車的受力及運(yùn)動(dòng)學(xué)模型進(jìn)行了簡(jiǎn)化與近似,且面對(duì)不同的列車與線路環(huán)境時(shí),需對(duì)照真實(shí)的運(yùn)行數(shù)據(jù)采用人工方式來調(diào)整模型參數(shù),存在著列車適配性差、仿真誤差大等問題,難以高效、準(zhǔn)確地模擬ATC在真實(shí)環(huán)境中的實(shí)際表現(xiàn),不利于系統(tǒng)的評(píng)估與改進(jìn)。

    針對(duì)上述問題,基于機(jī)器學(xué)習(xí)的理論與方法,本文通過挖掘列車在目標(biāo)環(huán)境下的運(yùn)行表現(xiàn),結(jié)合特征工程提高算法對(duì)ATC系統(tǒng)延時(shí)等特性的學(xué)習(xí)能力,基于XGBoost(極端梯度提升)算法建立具備更好適應(yīng)性與準(zhǔn)確性的列車運(yùn)行狀態(tài)仿真模型,為ATC系統(tǒng)的研發(fā)提供基礎(chǔ)技術(shù)支撐。

    1 現(xiàn)有的列車運(yùn)行狀態(tài)仿真方法概述

    現(xiàn)有的列車運(yùn)行狀態(tài)仿真系統(tǒng)一般以物理模型為基礎(chǔ),結(jié)合人工或求解式的參數(shù)調(diào)整方法,實(shí)現(xiàn)列車控制輸出與實(shí)際表現(xiàn)間函數(shù)關(guān)系的擬合,進(jìn)而對(duì)加速度、速度、位置等描述列車運(yùn)行狀態(tài)的屬性進(jìn)行估計(jì)與仿真,常采用包括單質(zhì)點(diǎn)模型、多質(zhì)點(diǎn)模型、繩體模型等描述方法擬合列車的動(dòng)力學(xué)表現(xiàn)。文獻(xiàn)[4]提出了一種針對(duì)高速列車運(yùn)行的仿真方法,使用繩體模型對(duì)列車進(jìn)行動(dòng)力學(xué)建模,能夠較好地模擬出列車質(zhì)量分布均衡、動(dòng)力分布分散的特點(diǎn),準(zhǔn)確模擬列車運(yùn)行狀態(tài),但由于不同列車的質(zhì)量與動(dòng)力分布特點(diǎn)差異較大等原因,該方法難以適配不同型號(hào)的列車和復(fù)雜運(yùn)行的線路。

    綜上,隨著仿真技術(shù)的不斷發(fā)展,基于物理模型的方法對(duì)輸出動(dòng)力的作用效果進(jìn)行估算的準(zhǔn)確性在持續(xù)提升,但仍缺乏適應(yīng)ATC系統(tǒng)輸出延時(shí)、適配不同列車性能特點(diǎn)的能力。因此,引入機(jī)器學(xué)習(xí)理論,研發(fā)可自動(dòng)化利用真實(shí)運(yùn)行數(shù)據(jù),構(gòu)建可準(zhǔn)確擬合控制命令與列車表現(xiàn)間復(fù)雜函數(shù)關(guān)系的仿真模型,以提高對(duì)列車運(yùn)行狀態(tài)的仿真能力。

    2 基于XGBoost算法的列車運(yùn)行狀態(tài)仿真方法

    2.1 算法原理

    本文提出一種基于XGBoost集成學(xué)習(xí)算法的列車運(yùn)行狀態(tài)仿真方法,其原理如圖1所示。該方法基于真實(shí)的運(yùn)行數(shù)據(jù),首先針對(duì)列車運(yùn)行過程中在不同速度、坡度及軌道狀態(tài)等條件下施加相同的牽引/制動(dòng)級(jí)位所產(chǎn)生不同加速度的作用效果,以及特定列車具有獨(dú)特且伴有一定隨機(jī)性的系統(tǒng)延時(shí)等問題,采用特征工程擴(kuò)展優(yōu)化原始數(shù)據(jù)維度,篩選并建立軌道信息、速度狀態(tài)以及級(jí)位序列的特征子集,合并構(gòu)成新的特征空間;在此基礎(chǔ)上,利用XGBoost算法,學(xué)習(xí)任意采樣時(shí)刻不同狀態(tài)下控制命令(包括歷史命令信息)與實(shí)際加速度間的關(guān)系,從而實(shí)現(xiàn)對(duì)下1個(gè)采樣時(shí)刻列車運(yùn)行狀態(tài)(借由速度、位置進(jìn)行描述)的估計(jì)。由此建立可與ATC在電子地圖輔助下進(jìn)行交互反饋的列車運(yùn)行狀態(tài)仿真模型。

    圖1 基于XGBoost模型的地鐵列車運(yùn)行狀態(tài)仿真方法原理Fig.1 Principle of metro train running state simulation method based on XGBoost model

    2.2 特征工程

    訓(xùn)練數(shù)據(jù)的樣本質(zhì)量與特征空間是決定機(jī)器學(xué)習(xí)效果的關(guān)鍵要素之一,因此首先對(duì)原始數(shù)據(jù)中存在的不完整數(shù)據(jù)(如采樣遺漏等)、異常值(如日志錯(cuò)誤等)進(jìn)行清洗,提取出完整的站間運(yùn)行數(shù)據(jù),然后針對(duì)列車運(yùn)行狀態(tài)仿真任務(wù),采用特征工程方法引導(dǎo)機(jī)器學(xué)習(xí)算法準(zhǔn)確挖掘目標(biāo)函數(shù)關(guān)系。

    列車運(yùn)行狀態(tài)仿真模型的核心任務(wù)是:觀察當(dāng)前采樣時(shí)刻的軌道坡度、列車速度、控制級(jí)位等數(shù)據(jù),基于隱含了列車運(yùn)動(dòng)模型、控制延時(shí)特性、級(jí)位作用表現(xiàn)等信息的作用效果估計(jì)模型,預(yù)測(cè)列車在下1個(gè)采樣時(shí)刻的速度表現(xiàn)。因此,特征空間應(yīng)包含軌道信息、速度狀態(tài)和級(jí)位序列,而訓(xùn)練標(biāo)簽可通過計(jì)算采樣時(shí)刻間的列車速度變化獲取。

    綜上,為建立可面向不同列車和不同線路環(huán)境、訓(xùn)練成本低、特征空間具有一定地鐵環(huán)境數(shù)據(jù)普適性的列車運(yùn)行狀態(tài)仿真模型,挖掘列車運(yùn)動(dòng)模型、控制延時(shí)特性、級(jí)位作用的表現(xiàn)特征,本文所述方法將當(dāng)前時(shí)刻的坡度數(shù)據(jù)、速度信息、歷史級(jí)位作為特征空間的主要組成。設(shè)t時(shí)刻列車的運(yùn)行特征為xt={et,vt,lt},其中:et為t時(shí)刻列車所處位置的坡度;vt為t時(shí)刻列車運(yùn)行速度;lt為t時(shí)刻列車輸出的級(jí)位序列{lt,0,lt,1,…,lt,k},lt,k為t-k時(shí)刻列車運(yùn)行系統(tǒng)輸出的級(jí)位。依據(jù)目標(biāo)列車的經(jīng)驗(yàn)最大延時(shí)周期數(shù),將k設(shè)為11,歷史時(shí)刻不足的采樣點(diǎn)則實(shí)施補(bǔ)零操作。同時(shí),將列車t+1時(shí)刻與t時(shí)刻的速度差at(即t時(shí)刻的加速度)設(shè)為樣本標(biāo)簽。最終構(gòu)建出13維的特征向量及1維的樣本標(biāo)簽,由此可依照有監(jiān)督學(xué)習(xí)方法進(jìn)行后續(xù)的模型訓(xùn)練及驗(yàn)證。

    2.3 模型構(gòu)建

    利用預(yù)處理后的數(shù)據(jù)集訓(xùn)練XGBoost模型,以實(shí)現(xiàn)對(duì)采樣點(diǎn)列車瞬時(shí)加速度的估計(jì)。XGBoost算法是文獻(xiàn)[5]提出的一種可擴(kuò)展的端到端基于樹的Boosting(提升方法)家族算法,可通過構(gòu)建多個(gè)弱學(xué)習(xí)器,將這些弱學(xué)習(xí)器組合形成強(qiáng)學(xué)習(xí)器。通常選取決策樹等子模型作為此算法的弱分類器。XGBoost算法因其優(yōu)異的泛化性能被廣泛應(yīng)用于有監(jiān)督學(xué)習(xí)領(lǐng)域,適合用于列車運(yùn)行狀態(tài)的仿真任務(wù),相比GBDT(梯度提升決策樹)方法,XGBoost算法在代價(jià)函數(shù)中加入了正則項(xiàng),可有效控制模型的復(fù)雜度,同時(shí)支持多線程訓(xùn)練。采用XGBoost算法后模型的計(jì)算速度更快,可滿足不同規(guī)模列車的運(yùn)行數(shù)據(jù)仿真要求。

    在模型訓(xùn)練中,XGBoost算法利用坡度數(shù)據(jù)、速度信息、歷史級(jí)位等參數(shù)自適應(yīng)地挖掘列車行駛過程中的加速度表現(xiàn)和控制延時(shí)特點(diǎn),以有效模擬列車在ATC系統(tǒng)下的真實(shí)運(yùn)行表現(xiàn)。與基于列車輸出級(jí)位和加速度間的對(duì)應(yīng)關(guān)系表(以下簡(jiǎn)稱“加速度表”)的傳統(tǒng)物理模型相比,XGBoost算法能夠適應(yīng)性地降低因環(huán)境因素干擾或列車實(shí)際參數(shù)與出廠參數(shù)不完全匹配等問題引發(fā)的仿真誤差。XGBoost算法的基本原理如下:

    對(duì)于樣本集N={(x1,y1),(x2,y2),…,(xn,yn)},訓(xùn)練優(yōu)化函數(shù)參數(shù)空間F,使得全局上目標(biāo)函數(shù)O(F)最小。目標(biāo)函數(shù)可表示為:

    O(F)=l(F)+Ω(F)

    (1)

    式中:

    Ω(F)——正則化項(xiàng),用以約束F的復(fù)雜性,控制模型的復(fù)雜度。

    XGBoost算法采用同Boosting算法家族一致的殘差修正思想。針對(duì)每個(gè)樣本,將現(xiàn)有估計(jì)與真實(shí)標(biāo)簽間的殘差作為即將加入的新弱學(xué)習(xí)器的訓(xùn)練標(biāo)簽,由此可得到加入第t棵樹后的估計(jì)結(jié)果為:

    (2)

    式中:

    xi——模型輸入的第i條樣本數(shù)據(jù);

    fi(xi)、ft(xi)——分別為第i棵樹和第t棵樹的分類結(jié)果;

    將損失函數(shù)和正則化項(xiàng)代入式(1)目標(biāo)函數(shù)中,可得第t次迭代時(shí)的目標(biāo)函數(shù)為:

    (3)

    式中:

    n——樣本數(shù)量;

    yi——第i條樣本數(shù)據(jù)對(duì)應(yīng)的真實(shí)標(biāo)簽;

    Ω(ft)——正則化項(xiàng),用以約束第t棵樹ft的復(fù)雜性,控制模型的復(fù)雜度。

    對(duì)目標(biāo)函數(shù)中的損失函數(shù)項(xiàng)進(jìn)行二階泰勒展開,并代入單個(gè)弱學(xué)習(xí)器的函數(shù)式與模型復(fù)雜度的函數(shù)式表達(dá),能夠得到展開的目標(biāo)函數(shù)如下:

    (4)

    式中:

    T——葉子節(jié)點(diǎn)的個(gè)數(shù);

    wj——葉子節(jié)點(diǎn)權(quán)重;

    K——構(gòu)建樹的總棵數(shù);

    γ、λ——均為調(diào)整參數(shù),用以防止過擬合;

    Gj——損失函數(shù)一階偏導(dǎo)的累加值之和;

    Hj——損失函數(shù)二階偏導(dǎo)的累加值之和。

    基于式(4)可計(jì)算每個(gè)葉子節(jié)點(diǎn)的最優(yōu)權(quán)重和模型最小損失,從而在迭代過程中使用貪婪算法逐步訓(xùn)練K棵樹,以組成具備強(qiáng)估計(jì)能力的集成模型。

    2.4 狀態(tài)仿真

    地鐵列車在運(yùn)行過程中會(huì)依據(jù)地圖信息和當(dāng)前運(yùn)行數(shù)據(jù)不斷切換其行車狀態(tài),并輸出相應(yīng)的制動(dòng)級(jí)位或牽引級(jí)位。在使用物理模型仿真時(shí),常常會(huì)依據(jù)列車出廠時(shí)提供的加速度表和級(jí)位傳輸延時(shí)信息來預(yù)測(cè)計(jì)算當(dāng)前時(shí)刻列車的加速度。具體的計(jì)算流程為:首先將列車系統(tǒng)輸出的級(jí)位依次存到級(jí)位盒中,然后依據(jù)列車延時(shí)表獲得當(dāng)前作用在列車上的級(jí)位,接著根據(jù)加速度表和當(dāng)前的列車速度找到對(duì)應(yīng)的理論加速度,最后依據(jù)式(5)更新下1個(gè)周期的列車速度,并依據(jù)式(6)更新該周期內(nèi)列車行駛距離。

    vt+1=vt+atT0

    (5)

    (6)

    式中:

    vt——t時(shí)刻的列車速度;

    T0——列車運(yùn)行的1個(gè)周期時(shí)長(zhǎng);

    st——t到t+1時(shí)刻內(nèi)列車的運(yùn)行距離。

    由于列車型號(hào)的不一致使得在模擬仿真過程中需要不斷地更換參數(shù),所以難以用一套固定的模型參數(shù)實(shí)現(xiàn)準(zhǔn)確的模擬仿真。此外,由于環(huán)境因素的干擾,使得列車在實(shí)際行駛過程出現(xiàn)與出廠參數(shù)不符的偏差。

    綜上,本文采用上述基于XGBoost算法構(gòu)建的列車速度估計(jì)模型,替換傳統(tǒng)物理模型的加速度計(jì)算模塊。仿真時(shí)先對(duì)列車的線路信息及狀態(tài)進(jìn)行初始化,并進(jìn)行與訓(xùn)練數(shù)據(jù)一致的特征工程處理;然后將采樣點(diǎn)輸入至XGBoost模型中,獲得當(dāng)前時(shí)刻的列車加速度,同理依據(jù)式(5)~(6)更新列車的速度和行駛距離;接著利用地圖信息、ATC系統(tǒng)輸出更新的坡度信息和級(jí)位序列。上述步驟依次迭代進(jìn)行,直到列車到站停車。最后,將列車在線路上的運(yùn)行仿真數(shù)據(jù)進(jìn)行記錄,并分析其仿真效果。

    3 基于XGBoost算法的列車運(yùn)行狀態(tài)仿真試驗(yàn)

    3.1 試驗(yàn)?zāi)康?/h3>

    為驗(yàn)證基于XGBoost算法的列車運(yùn)行狀態(tài)仿真方法的有效性,本文對(duì)3列來自某地鐵線路的真實(shí)列車數(shù)據(jù)進(jìn)行學(xué)習(xí)及仿真測(cè)試,并與基于性能參數(shù)表的物理模型仿真方法進(jìn)行比較。數(shù)據(jù)特征參數(shù)的具體說明如表1所示。

    表1 模型中采用的列車特征參數(shù)及其單位

    3.2 評(píng)價(jià)方法

    采用每個(gè)采樣周期內(nèi)加速度的均方根誤差R和每站速度曲線誤差百分比E兩個(gè)指標(biāo),其計(jì)算式如式(7)、式(8)所示。R和E的值越小,表示模型的準(zhǔn)確性越高,與實(shí)際列車運(yùn)行情況的符合程度越高,即仿真效果越好。

    (7)

    (8)

    式中:

    N——數(shù)據(jù)集中的數(shù)據(jù)總條數(shù);

    ai——列車的實(shí)際加速度;

    m——線路上的數(shù)據(jù)采樣數(shù);

    vj——j時(shí)刻列車的實(shí)際速度;

    3.3 試驗(yàn)過程

    首先對(duì)原始數(shù)據(jù)進(jìn)行清洗,對(duì)加速度異常、缺失的數(shù)據(jù)點(diǎn)進(jìn)行刪減處理,然后依據(jù)算法原理中的特征工程方法逐條完成特征空間的建立。接著預(yù)留1個(gè)車站的數(shù)據(jù)(該站不參與模型訓(xùn)練)并用作單站列車運(yùn)行速度曲線仿真,以驗(yàn)證模型的實(shí)用性。將其余站點(diǎn)的采樣數(shù)據(jù)按照7∶2∶1的比例分別拆分至訓(xùn)練集、驗(yàn)證集和測(cè)試集內(nèi)。

    在進(jìn)行列車運(yùn)行狀態(tài)仿真模型訓(xùn)練時(shí),首先將訓(xùn)練集輸入至XGBoost模型中,并用驗(yàn)證集進(jìn)行超參數(shù)調(diào)整,最后使用測(cè)試集對(duì)仿真模型的性能進(jìn)行分析評(píng)價(jià)。此外,將基于XGBoost模型的仿真結(jié)果與基于列車性能參數(shù)表物理模型得到的仿真結(jié)果進(jìn)行對(duì)比。

    基于驗(yàn)證集的表現(xiàn),確定試驗(yàn)中采用的XGBoost模型超參數(shù)如下:最大樹深度為8,學(xué)習(xí)率為0.6,基分類器數(shù)量為500。在實(shí)用性檢測(cè)時(shí),將基于XGBoost算法的仿真模型應(yīng)用到?jīng)]參加模型訓(xùn)練的預(yù)留站上并進(jìn)行模擬仿真,將預(yù)留站列車的實(shí)際運(yùn)行速度、基于XGBoost模型得到的仿真速度、基于物理模型得到的仿真速度進(jìn)行對(duì)比。

    3.4 試驗(yàn)結(jié)果分析

    參與試驗(yàn)的3列車的編號(hào)為T1、T2、T3,數(shù)據(jù)量分別為460 511條、795 88條、472 68條。將基于物理模型得到的仿真結(jié)果與基于XGBoost模型得到的結(jié)果進(jìn)行對(duì)比,如表2所示。從表2中可以看出,與基于物理模型得到的結(jié)果相比,基于XGBoost模型的列車運(yùn)行狀態(tài)仿真模型的R和E在多列車上的測(cè)試結(jié)果均較低,且在不同列車上的仿真性能表現(xiàn)更為穩(wěn)定,這說明了基于XGBoost模型的仿真方法具有良好的環(huán)境適應(yīng)性,能與不同的列車車型相匹配。

    表2 基于物理模型和基于XGBoost模型的仿真結(jié)果對(duì)比

    將預(yù)留站列車的實(shí)際運(yùn)行速度、基于XGBoost模型得到的仿真速度、基于物理模型得到的仿真速度進(jìn)行對(duì)比,得到的結(jié)果如圖2所示。從圖2可以看出,基于物理模型得到的仿真速度在模擬控制級(jí)位變換較為頻繁,且在復(fù)雜的目標(biāo)速度調(diào)整(采樣周期為(200~400)×0.2 s)及巡航速度控制(采樣周期為(400~1 300)×0.2 s)階段均難以實(shí)現(xiàn)準(zhǔn)確的速度估計(jì),從而影響了整個(gè)車站的仿真效果。而基于XGBoost模型得到的仿真速度曲線則與列車的實(shí)際運(yùn)行曲線更為貼合。

    圖2 3種不同測(cè)試工況下列車在單站運(yùn)行時(shí)的 速度曲線仿真對(duì)比

    綜上可知,本文提出的基于XGBoost的列車運(yùn)行速度仿真模型構(gòu)建方法可從實(shí)際運(yùn)行數(shù)據(jù)中有效挖掘列車運(yùn)動(dòng)模型、控制延時(shí)特性、級(jí)位作用表現(xiàn)等特征并加以利用,進(jìn)而建立列車運(yùn)行狀態(tài)的仿真模型。該模型具有良好的環(huán)境適應(yīng)性,能與不同的列車車型相匹配,可生成與實(shí)際情況接近的列車運(yùn)行速度曲線,較為準(zhǔn)確地實(shí)現(xiàn)對(duì)地鐵列車運(yùn)行狀態(tài)的仿真。

    4 結(jié)語

    列車運(yùn)行狀態(tài)仿真是研發(fā)與優(yōu)化ATC的重要工具,具備自動(dòng)化適應(yīng)能力與高精度仿真性能的狀態(tài)仿真模型將助力軌道交通自動(dòng)化技術(shù)的不斷發(fā)展。本文提出的基于XGBoost模型的仿真方法可從實(shí)際運(yùn)行數(shù)據(jù)中有效挖掘列車運(yùn)動(dòng)模型、控制延時(shí)特性、級(jí)位作用表現(xiàn)等特征,實(shí)現(xiàn)對(duì)ATC系統(tǒng)輸出與實(shí)際作用效果間復(fù)雜函數(shù)關(guān)系的學(xué)習(xí),進(jìn)而建立具備一定適應(yīng)性且模擬精度較高的列車運(yùn)行狀態(tài)仿真模型。試驗(yàn)采用真實(shí)數(shù)據(jù)驗(yàn)證了使用該方法所建立的仿真模型的有效性,且該方法適用于針對(duì)不同列車車型的運(yùn)行狀態(tài)仿真場(chǎng)景。未來應(yīng)進(jìn)一步結(jié)合其他先進(jìn)的序列數(shù)據(jù)建模方法,不斷改進(jìn)、優(yōu)化該模型的訓(xùn)練效果,為列車自動(dòng)運(yùn)行技術(shù)的發(fā)展提供更可靠、有效的技術(shù)支撐。

    猜你喜歡
    列車運(yùn)行加速度列車
    “鱉”不住了!從26元/斤飆至38元/斤,2022年甲魚能否再跑出“加速度”?
    登上末日列車
    關(guān)愛向列車下延伸
    改善地鐵列車運(yùn)行舒適度方案探討
    穿越時(shí)空的列車
    天際加速度
    汽車觀察(2018年12期)2018-12-26 01:05:42
    創(chuàng)新,動(dòng)能轉(zhuǎn)換的“加速度”
    金橋(2018年4期)2018-09-26 02:24:46
    死亡加速度
    列車運(yùn)行控制系統(tǒng)技術(shù)發(fā)展趨勢(shì)分析
    相同徑路的高速列車運(yùn)行圖編制方法
    亚洲人成电影观看| 永久免费av网站大全| 精品人妻一区二区三区麻豆| 黄色 视频免费看| 亚洲男人天堂网一区| 欧美日本中文国产一区发布| www.精华液| 三级毛片av免费| 大陆偷拍与自拍| 91成人精品电影| 成人18禁高潮啪啪吃奶动态图| 欧美日韩成人在线一区二区| 伊人亚洲综合成人网| 免费高清在线观看视频在线观看| 免费一级毛片在线播放高清视频 | 亚洲中文日韩欧美视频| 国产不卡av网站在线观看| 美女高潮喷水抽搐中文字幕| 欧美黑人精品巨大| 欧美精品高潮呻吟av久久| 最近中文字幕2019免费版| 91国产中文字幕| 亚洲熟女精品中文字幕| 亚洲av日韩在线播放| 天天添夜夜摸| 国产精品免费视频内射| 满18在线观看网站| 国产区一区二久久| 亚洲成av片中文字幕在线观看| 亚洲欧美色中文字幕在线| 嫩草影视91久久| 午夜精品国产一区二区电影| 国产无遮挡羞羞视频在线观看| 国产男女内射视频| 性色av乱码一区二区三区2| 麻豆av在线久日| 亚洲精华国产精华精| 国产熟女午夜一区二区三区| 人妻一区二区av| 国产又色又爽无遮挡免| 又紧又爽又黄一区二区| 国产欧美日韩综合在线一区二区| 日本av免费视频播放| 精品免费久久久久久久清纯 | 日韩视频一区二区在线观看| 久久精品国产亚洲av香蕉五月 | kizo精华| 99国产精品一区二区三区| 亚洲一区中文字幕在线| 自线自在国产av| 亚洲性夜色夜夜综合| 黄色a级毛片大全视频| 50天的宝宝边吃奶边哭怎么回事| 搡老熟女国产l中国老女人| a级毛片黄视频| 精品国产超薄肉色丝袜足j| 在线观看www视频免费| 免费在线观看影片大全网站| 亚洲少妇的诱惑av| 一区福利在线观看| 日本黄色日本黄色录像| 天堂中文最新版在线下载| 亚洲av欧美aⅴ国产| 国产成人一区二区三区免费视频网站| 免费av中文字幕在线| 啦啦啦在线免费观看视频4| 午夜激情av网站| 欧美性长视频在线观看| 丰满迷人的少妇在线观看| 日韩免费高清中文字幕av| 久久精品国产a三级三级三级| 久久亚洲国产成人精品v| 中亚洲国语对白在线视频| 大片电影免费在线观看免费| 国产精品熟女久久久久浪| 亚洲伊人色综图| 人人妻人人澡人人爽人人夜夜| 菩萨蛮人人尽说江南好唐韦庄| xxxhd国产人妻xxx| 国产一级毛片在线| 久久久水蜜桃国产精品网| 国产精品成人在线| 亚洲 欧美一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 久久精品熟女亚洲av麻豆精品| 男女下面插进去视频免费观看| 久9热在线精品视频| 久久 成人 亚洲| 精品国产超薄肉色丝袜足j| 亚洲精品美女久久久久99蜜臀| 久久九九热精品免费| 人人澡人人妻人| 久久久久久免费高清国产稀缺| 十八禁网站网址无遮挡| 久9热在线精品视频| 真人做人爱边吃奶动态| 亚洲专区国产一区二区| 狠狠精品人妻久久久久久综合| 国产成人精品在线电影| 久久久欧美国产精品| 悠悠久久av| 人妻 亚洲 视频| 性高湖久久久久久久久免费观看| 91精品三级在线观看| 亚洲精品美女久久久久99蜜臀| 国产有黄有色有爽视频| 色精品久久人妻99蜜桃| e午夜精品久久久久久久| 日本av手机在线免费观看| 欧美性长视频在线观看| 日本五十路高清| 国产在线视频一区二区| 丝袜人妻中文字幕| 精品欧美一区二区三区在线| 少妇猛男粗大的猛烈进出视频| 一级黄色大片毛片| 亚洲欧美色中文字幕在线| videosex国产| 亚洲人成电影观看| 女警被强在线播放| 熟女少妇亚洲综合色aaa.| 精品高清国产在线一区| 婷婷色av中文字幕| 考比视频在线观看| 日韩欧美免费精品| 香蕉国产在线看| 青青草视频在线视频观看| 午夜精品国产一区二区电影| 欧美黑人欧美精品刺激| 亚洲色图 男人天堂 中文字幕| 啦啦啦中文免费视频观看日本| 97人妻天天添夜夜摸| 精品人妻在线不人妻| 欧美日韩亚洲综合一区二区三区_| 最新的欧美精品一区二区| 免费av中文字幕在线| 国产无遮挡羞羞视频在线观看| 女人久久www免费人成看片| 操美女的视频在线观看| 黑人欧美特级aaaaaa片| 婷婷色av中文字幕| 天天添夜夜摸| 亚洲国产中文字幕在线视频| 国产精品免费视频内射| 捣出白浆h1v1| 国产精品秋霞免费鲁丝片| 一本久久精品| 日韩欧美一区二区三区在线观看 | 97精品久久久久久久久久精品| 国产免费av片在线观看野外av| 欧美日韩视频精品一区| 大片免费播放器 马上看| 久久国产精品男人的天堂亚洲| 国产精品 国内视频| 国产精品欧美亚洲77777| 丝袜在线中文字幕| 视频区欧美日本亚洲| 最新在线观看一区二区三区| 国产精品久久久久久精品电影小说| 欧美黄色片欧美黄色片| 亚洲熟女精品中文字幕| 亚洲欧美清纯卡通| 国产主播在线观看一区二区| 中亚洲国语对白在线视频| 亚洲国产精品一区三区| 12—13女人毛片做爰片一| 亚洲av电影在线进入| 成年女人毛片免费观看观看9 | 多毛熟女@视频| 国产一级毛片在线| 免费观看人在逋| 午夜激情久久久久久久| 亚洲色图综合在线观看| 国产精品国产三级国产专区5o| 捣出白浆h1v1| 精品欧美一区二区三区在线| 激情视频va一区二区三区| 亚洲av日韩精品久久久久久密| 一本色道久久久久久精品综合| 亚洲欧美精品综合一区二区三区| 99精品欧美一区二区三区四区| 三级毛片av免费| 成年人午夜在线观看视频| 啦啦啦免费观看视频1| a在线观看视频网站| 丁香六月欧美| 日韩电影二区| 国产av又大| 大片免费播放器 马上看| 在线观看免费视频网站a站| 国产精品成人在线| 夜夜骑夜夜射夜夜干| 精品人妻在线不人妻| 99热国产这里只有精品6| 亚洲自偷自拍图片 自拍| 国产欧美日韩一区二区三 | 亚洲五月色婷婷综合| 色94色欧美一区二区| av网站在线播放免费| 在线观看www视频免费| 国产有黄有色有爽视频| 国产人伦9x9x在线观看| 精品人妻1区二区| 亚洲精品久久久久久婷婷小说| 啦啦啦视频在线资源免费观看| 人人妻,人人澡人人爽秒播| 人成视频在线观看免费观看| 美女扒开内裤让男人捅视频| 一边摸一边抽搐一进一出视频| 宅男免费午夜| 亚洲欧美成人综合另类久久久| 国产老妇伦熟女老妇高清| 欧美黄色片欧美黄色片| 精品久久久久久电影网| 9热在线视频观看99| 国产精品国产av在线观看| 国产色视频综合| 日韩免费高清中文字幕av| 欧美黄色片欧美黄色片| 亚洲第一欧美日韩一区二区三区 | 黑人巨大精品欧美一区二区蜜桃| 视频在线观看一区二区三区| av电影中文网址| 午夜福利,免费看| 免费不卡黄色视频| 亚洲精品一卡2卡三卡4卡5卡 | 成人18禁高潮啪啪吃奶动态图| 亚洲专区中文字幕在线| 免费女性裸体啪啪无遮挡网站| 黄色a级毛片大全视频| 国产日韩欧美在线精品| 亚洲天堂av无毛| 又黄又粗又硬又大视频| 中文字幕人妻丝袜一区二区| 一区在线观看完整版| 真人做人爱边吃奶动态| 久久久久久免费高清国产稀缺| 69精品国产乱码久久久| 精品第一国产精品| 欧美亚洲日本最大视频资源| 久久久精品免费免费高清| 成在线人永久免费视频| 精品熟女少妇八av免费久了| av在线播放精品| 亚洲少妇的诱惑av| 国产91精品成人一区二区三区 | 黄色视频,在线免费观看| 99久久国产精品久久久| netflix在线观看网站| 久久国产精品男人的天堂亚洲| 欧美在线一区亚洲| 国产精品欧美亚洲77777| 少妇粗大呻吟视频| 亚洲av电影在线进入| 高清黄色对白视频在线免费看| 一区二区三区精品91| 黑人欧美特级aaaaaa片| 国产精品国产三级国产专区5o| 国产免费av片在线观看野外av| 国产一区二区三区av在线| 窝窝影院91人妻| 人人妻,人人澡人人爽秒播| 免费观看av网站的网址| 自拍欧美九色日韩亚洲蝌蚪91| 日韩 亚洲 欧美在线| 在线精品无人区一区二区三| 别揉我奶头~嗯~啊~动态视频 | 国产精品一区二区在线观看99| 18在线观看网站| 欧美久久黑人一区二区| 12—13女人毛片做爰片一| 五月开心婷婷网| 精品久久久久久电影网| 五月天丁香电影| 国产一卡二卡三卡精品| 日韩人妻精品一区2区三区| 午夜久久久在线观看| 免费在线观看影片大全网站| 亚洲一码二码三码区别大吗| 中文字幕高清在线视频| 免费在线观看完整版高清| 各种免费的搞黄视频| 另类精品久久| 高清黄色对白视频在线免费看| 色播在线永久视频| 一本一本久久a久久精品综合妖精| 久久久国产一区二区| 男人舔女人的私密视频| 亚洲人成电影观看| 欧美日韩av久久| 天天添夜夜摸| 久久国产亚洲av麻豆专区| a级片在线免费高清观看视频| 国产成人欧美在线观看 | 亚洲欧美日韩另类电影网站| 久久人妻福利社区极品人妻图片| 两个人看的免费小视频| 日韩精品免费视频一区二区三区| 国产精品亚洲av一区麻豆| 一级黄色大片毛片| 人成视频在线观看免费观看| 无限看片的www在线观看| www.999成人在线观看| 国产高清videossex| 精品国产一区二区三区四区第35| 极品人妻少妇av视频| 妹子高潮喷水视频| 男女边摸边吃奶| 国产老妇伦熟女老妇高清| 人人妻人人澡人人爽人人夜夜| 一级毛片女人18水好多| 精品人妻在线不人妻| av又黄又爽大尺度在线免费看| 国产欧美日韩一区二区三区在线| 亚洲精品久久久久久婷婷小说| 一本综合久久免费| 十八禁网站网址无遮挡| 久久综合国产亚洲精品| 日韩电影二区| 国产男人的电影天堂91| 亚洲国产毛片av蜜桃av| 夜夜夜夜夜久久久久| 热re99久久精品国产66热6| 国产欧美亚洲国产| 国产激情久久老熟女| 国产野战对白在线观看| 后天国语完整版免费观看| 岛国在线观看网站| 精品免费久久久久久久清纯 | 国产成人av激情在线播放| 久久天堂一区二区三区四区| 午夜激情久久久久久久| 天天影视国产精品| 可以免费在线观看a视频的电影网站| 不卡一级毛片| 纵有疾风起免费观看全集完整版| 一本久久精品| 性高湖久久久久久久久免费观看| 亚洲精品av麻豆狂野| 一二三四社区在线视频社区8| 精品第一国产精品| 亚洲第一青青草原| 男女国产视频网站| 啦啦啦免费观看视频1| 少妇的丰满在线观看| 精品亚洲成a人片在线观看| 日韩有码中文字幕| 黑人操中国人逼视频| www.av在线官网国产| 免费在线观看日本一区| 久久中文看片网| 国产亚洲欧美精品永久| 国产又爽黄色视频| 一本—道久久a久久精品蜜桃钙片| 人人妻,人人澡人人爽秒播| 黄色片一级片一级黄色片| 在线亚洲精品国产二区图片欧美| 国产精品偷伦视频观看了| 亚洲午夜精品一区,二区,三区| 美女脱内裤让男人舔精品视频| 亚洲欧美一区二区三区久久| 99国产综合亚洲精品| 精品高清国产在线一区| 一个人免费在线观看的高清视频 | 亚洲精品美女久久久久99蜜臀| 久久久久久久久免费视频了| 男女床上黄色一级片免费看| 亚洲欧美日韩高清在线视频 | 波多野结衣一区麻豆| 一区二区三区激情视频| 欧美精品一区二区免费开放| 国产视频一区二区在线看| 丝袜人妻中文字幕| 激情视频va一区二区三区| 中文欧美无线码| 999久久久精品免费观看国产| 日韩大码丰满熟妇| 国产精品免费视频内射| 中亚洲国语对白在线视频| 狂野欧美激情性bbbbbb| 最新的欧美精品一区二区| 99热全是精品| 超色免费av| 国产视频一区二区在线看| 超色免费av| 999久久久精品免费观看国产| 如日韩欧美国产精品一区二区三区| 亚洲中文字幕日韩| 亚洲熟女精品中文字幕| 亚洲精品国产av蜜桃| 99久久国产精品久久久| 女人高潮潮喷娇喘18禁视频| 夜夜夜夜夜久久久久| 无遮挡黄片免费观看| 久久久久国产一级毛片高清牌| 国产精品麻豆人妻色哟哟久久| 精品欧美一区二区三区在线| 久久亚洲精品不卡| 免费av中文字幕在线| 亚洲视频免费观看视频| 两个人看的免费小视频| 啦啦啦 在线观看视频| 久久国产亚洲av麻豆专区| 91成年电影在线观看| 丝袜美足系列| 国产亚洲精品久久久久5区| www日本在线高清视频| 欧美精品人与动牲交sv欧美| 亚洲欧美激情在线| 国产一区二区三区av在线| 成人av一区二区三区在线看 | 亚洲国产毛片av蜜桃av| 亚洲国产中文字幕在线视频| 国产91精品成人一区二区三区 | 美女视频免费永久观看网站| 国产成人av教育| av在线app专区| 成人国产av品久久久| 91麻豆精品激情在线观看国产 | 国内毛片毛片毛片毛片毛片| 国产在视频线精品| 国产男女超爽视频在线观看| 99国产精品99久久久久| 国产伦理片在线播放av一区| 黄色a级毛片大全视频| 日韩三级视频一区二区三区| 精品卡一卡二卡四卡免费| 国产av精品麻豆| 丝袜美足系列| 成年美女黄网站色视频大全免费| 成年人午夜在线观看视频| 伊人亚洲综合成人网| 97精品久久久久久久久久精品| 中文欧美无线码| 午夜精品国产一区二区电影| 18禁裸乳无遮挡动漫免费视频| 精品视频人人做人人爽| 大码成人一级视频| 美女午夜性视频免费| 18禁裸乳无遮挡动漫免费视频| 一区福利在线观看| 50天的宝宝边吃奶边哭怎么回事| 视频区欧美日本亚洲| 又黄又粗又硬又大视频| 久久久久久久久久久久大奶| 国产亚洲av高清不卡| 69av精品久久久久久 | 亚洲精品久久成人aⅴ小说| 国产精品免费大片| 黄片小视频在线播放| 久久99热这里只频精品6学生| 国产在线免费精品| 欧美乱码精品一区二区三区| 亚洲精品中文字幕一二三四区 | 欧美午夜高清在线| 中文字幕av电影在线播放| 久久人妻福利社区极品人妻图片| 一二三四在线观看免费中文在| 嫁个100分男人电影在线观看| 999精品在线视频| 99久久国产精品久久久| 精品国产一区二区三区四区第35| 久久久久久久久免费视频了| 一个人免费在线观看的高清视频 | 亚洲性夜色夜夜综合| 成年人午夜在线观看视频| 嫁个100分男人电影在线观看| 亚洲精品国产av蜜桃| 各种免费的搞黄视频| 国产又色又爽无遮挡免| 午夜福利视频精品| 欧美xxⅹ黑人| 在线观看www视频免费| 国产麻豆69| 日韩电影二区| 久久99热这里只频精品6学生| 99精国产麻豆久久婷婷| 中文字幕精品免费在线观看视频| 18在线观看网站| 999久久久国产精品视频| 日韩一区二区三区影片| 国产精品秋霞免费鲁丝片| 国产精品久久久av美女十八| 两个人看的免费小视频| 国产亚洲欧美在线一区二区| av不卡在线播放| 中文字幕人妻丝袜制服| 亚洲午夜精品一区,二区,三区| 欧美日韩成人在线一区二区| 欧美激情久久久久久爽电影 | 精品一区二区三区四区五区乱码| 精品亚洲成a人片在线观看| 一本色道久久久久久精品综合| 久久久欧美国产精品| 欧美日韩成人在线一区二区| 精品人妻1区二区| 多毛熟女@视频| 19禁男女啪啪无遮挡网站| av不卡在线播放| 大香蕉久久成人网| 悠悠久久av| 精品久久蜜臀av无| 日韩视频在线欧美| 亚洲精品久久午夜乱码| 国产亚洲欧美在线一区二区| 国产免费一区二区三区四区乱码| 日本a在线网址| 少妇猛男粗大的猛烈进出视频| 国产成人精品无人区| 老汉色∧v一级毛片| 99久久精品国产亚洲精品| 成人国产一区最新在线观看| 亚洲综合色网址| 十八禁网站免费在线| 亚洲av成人一区二区三| av网站免费在线观看视频| 国产高清视频在线播放一区 | 国产主播在线观看一区二区| 黄色怎么调成土黄色| 美女脱内裤让男人舔精品视频| 国产视频一区二区在线看| 国产一区二区在线观看av| av视频免费观看在线观看| 亚洲精品av麻豆狂野| 69av精品久久久久久 | 久久精品人人爽人人爽视色| 欧美精品啪啪一区二区三区 | 99热全是精品| 国产成人免费观看mmmm| 日韩 欧美 亚洲 中文字幕| 免费在线观看日本一区| 亚洲国产精品一区二区三区在线| 悠悠久久av| 一级片'在线观看视频| 久久亚洲国产成人精品v| 淫妇啪啪啪对白视频 | 日本av手机在线免费观看| 大陆偷拍与自拍| 国产日韩欧美视频二区| 欧美黄色淫秽网站| 在线永久观看黄色视频| 亚洲欧美一区二区三区久久| 久久精品国产亚洲av高清一级| 成年人午夜在线观看视频| a级毛片在线看网站| 亚洲欧美一区二区三区久久| 亚洲av电影在线观看一区二区三区| 国产精品国产av在线观看| 淫妇啪啪啪对白视频 | 亚洲欧美成人综合另类久久久| 天堂俺去俺来也www色官网| av福利片在线| 一个人免费看片子| 久久久久精品国产欧美久久久 | 婷婷丁香在线五月| avwww免费| 成人影院久久| 国产一区二区三区综合在线观看| 美女福利国产在线| 性高湖久久久久久久久免费观看| 色播在线永久视频| 99国产精品一区二区蜜桃av | 热re99久久国产66热| 国产精品一区二区在线观看99| 99国产精品一区二区三区| 亚洲精品日韩在线中文字幕| 999久久久国产精品视频| 欧美一级毛片孕妇| 考比视频在线观看| 日日夜夜操网爽| 大陆偷拍与自拍| 日韩 亚洲 欧美在线| 99香蕉大伊视频| 老熟女久久久| 亚洲精品日韩在线中文字幕| 每晚都被弄得嗷嗷叫到高潮| 老司机午夜十八禁免费视频| 午夜精品国产一区二区电影| 亚洲国产av新网站| 国产成人系列免费观看| 一级片'在线观看视频| 黄片小视频在线播放| 日韩中文字幕视频在线看片| 午夜福利一区二区在线看| 日日爽夜夜爽网站| 亚洲七黄色美女视频| 亚洲精品第二区| 午夜老司机福利片| 999久久久精品免费观看国产| 精品久久久久久电影网| 日韩精品免费视频一区二区三区| 色视频在线一区二区三区| 大陆偷拍与自拍| 一边摸一边抽搐一进一出视频| 永久免费av网站大全| 亚洲五月色婷婷综合| 99热国产这里只有精品6| 午夜免费成人在线视频| 久久久久国产一级毛片高清牌| 另类亚洲欧美激情| 日日夜夜操网爽| 中文字幕人妻丝袜制服| 侵犯人妻中文字幕一二三四区| 久久精品亚洲av国产电影网| 欧美成人午夜精品| 女警被强在线播放| 啦啦啦免费观看视频1| 久久国产精品人妻蜜桃| 91成年电影在线观看| 亚洲国产精品成人久久小说| 久久天躁狠狠躁夜夜2o2o| 亚洲熟女毛片儿| 亚洲av男天堂| av视频免费观看在线观看| 亚洲精品av麻豆狂野|