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

    基于列車備選集的技術(shù)站單組列車編組計(jì)劃優(yōu)化模型

    2021-04-02 05:58:38魏玉光
    鐵道學(xué)報(bào) 2021年10期
    關(guān)鍵詞:優(yōu)化模型

    方 波,魏玉光,楊 浩

    (北京交通大學(xué) 交通運(yùn)輸學(xué)院, 北京 100044)

    技術(shù)站列車編組計(jì)劃作為貨物列車編組計(jì)劃的重要組成部分[1],是一個(gè)復(fù)雜的大規(guī)模組合優(yōu)化問題。而單組列車作為我國技術(shù)站貨物列車的主要組織形式,承擔(dān)了大部分技術(shù)站間的車流運(yùn)輸任務(wù)。因此,技術(shù)站單組列車編組計(jì)劃的優(yōu)化是提高路網(wǎng)車流運(yùn)輸組織水平的關(guān)鍵。

    大規(guī)模實(shí)際路網(wǎng)涉及范圍廣、因素多,使得單組列車編組計(jì)劃成為NP-hard問題,用現(xiàn)有模型和方法幾乎難以得到最優(yōu)解,簡化實(shí)際問題并求解出符合運(yùn)輸需要的編組計(jì)劃成為研究重點(diǎn)。國內(nèi)外學(xué)者圍繞技術(shù)站單組列車編組計(jì)劃開展了大量研究,研究成果主要集中于車流徑路優(yōu)化、單組列車編組計(jì)劃優(yōu)化、單組列車與分組列車編組計(jì)劃協(xié)同優(yōu)化、車流徑路與單組列車編組計(jì)劃的綜合優(yōu)化等方面。

    車流徑路優(yōu)化方面:江南等[2]提出適合鐵路管理方式的車流徑路制定參照模型及算法,并針對(duì)各種車流徑路的求解方法給出詳細(xì)的數(shù)學(xué)模型與算法;溫旭紅等[3]以路網(wǎng)上車流的廣義里程成本最小化為目標(biāo),以車流徑路的樹形結(jié)構(gòu)約束、車流強(qiáng)度守恒以及路段通過能力限制等為約束條件,提出基于樹形結(jié)構(gòu)的鐵路車流徑路優(yōu)化模型。

    單組列車編組計(jì)劃優(yōu)化方面:史峰等[4]將給定編組去向集時(shí)有改編能力約束的編組計(jì)劃問題轉(zhuǎn)化為車流徑路問題,提出交替添加、剔除編組去向的啟發(fā)式算法;許紅等[5]綜合考慮車站編組能力、解體能力、調(diào)車線容車數(shù)等影響因素,構(gòu)建了以技術(shù)站車流總消耗最小及改編能力利用程度最均衡的多目標(biāo)0-1規(guī)劃模型,并以改進(jìn)的遺傳算法求解;林柏梁等[6]基于最遠(yuǎn)站法則,構(gòu)建了列車編組計(jì)劃的雙層模型,上層確定列車到發(fā)站,下層確定車流改編方案,并采用模擬退火算法對(duì)模型進(jìn)行求解;Chen等[7]假設(shè)車流徑路已知,構(gòu)建了基于路徑的單組列車編組計(jì)劃線性模型,以車站最大編組去向數(shù)代替分類線數(shù)量約束來避免模型的非線性,并設(shè)計(jì)了一種啟發(fā)式的樹分解算法進(jìn)行求解。

    單組列車與分組列車編組計(jì)劃協(xié)同優(yōu)化方面:梁棟等[8]分析了分組列車相較于單組列車的耗費(fèi)節(jié)省,在單組列車編組計(jì)劃的基礎(chǔ)上,以分組列車取代單組列車的節(jié)省最大為目標(biāo)求解分組列車編組計(jì)劃模型;肖杰等[9]建立了單組列車和分組列車編組計(jì)劃的綜合優(yōu)化模型,以遺傳算法求解得到的單組列車編組計(jì)劃為基礎(chǔ),采用遍歷算法求解綜合優(yōu)化模型;另外,Xiao等[10]開發(fā)了一種基于遺傳算法和禁忌搜索的混合算法來求解單組列車編組計(jì)劃模型,然后將分組列車替換為單組列車,并提出一種貪婪算法求解綜合優(yōu)化模型。

    車流徑路與單組列車綜合優(yōu)化方面:李德義[11]提出另一種列車備選集方法,構(gòu)建以技術(shù)站車流集結(jié)、改編中轉(zhuǎn)、無改編中轉(zhuǎn)和走行總消耗最小的車流徑路與編組計(jì)劃綜合優(yōu)化模型,并設(shè)計(jì)基于分支定界和列生成算法的分支定價(jià)算法進(jìn)行求解;林柏梁等[12]提出了將車流徑路選擇、裝車地直達(dá)列車編組計(jì)劃及技術(shù)站直達(dá)列車編組計(jì)劃三者統(tǒng)一優(yōu)化的模型;嚴(yán)余松等[13]根據(jù)機(jī)會(huì)約束規(guī)劃理論,建立了基于車流量波動(dòng)的列車編組計(jì)劃與車流徑路綜合優(yōu)化模型,并設(shè)計(jì)基于隨機(jī)模擬的改進(jìn)分枝定界法對(duì)模型進(jìn)行求解。

    文獻(xiàn)[4-10,12-13]可歸納為“車流-車組-列車”的遞歸方法,即先將車流合并為車組(編組去向),再將車組分配給列車,大多基于點(diǎn)-弧或弧-路模型建??蚣埽ㄟ^引入決策變量表示車流是否在某站改編來遞歸求解車流的完整改編過程。這種方法容易導(dǎo)致模型的非線性結(jié)構(gòu),求解難度高,依賴于啟發(fā)式算法,并且車流遞歸難以直觀展示各支車流的改編過程。文獻(xiàn)[11]可概括為“車流-列車”的研究方法,模型基于多商品網(wǎng)絡(luò)流框架,車流改編接續(xù)過程直觀,但冗余變量過多,求解結(jié)果不理想。

    針對(duì)現(xiàn)有研究存在的不足,本文基于多商品流服務(wù)網(wǎng)絡(luò)設(shè)計(jì)原理,提出列車備選集的方法,首先構(gòu)建列車備選集,再進(jìn)行車流分配,把單組列車編組計(jì)劃問題轉(zhuǎn)化為列車服務(wù)網(wǎng)絡(luò)上的車流分配問題,并建立了兩階段線性規(guī)劃模型分別求解最優(yōu)車流徑路和最優(yōu)編組計(jì)劃。第一階段的車流徑路模型用于求解各OD對(duì)的最優(yōu)車流徑路;第二階段為編組計(jì)劃模型,以上一階段所得最優(yōu)車流徑路生成的列車備選集為輸入,求解最優(yōu)列車服務(wù)和車流分配方案。最后分別以小規(guī)模路網(wǎng)和大規(guī)模路網(wǎng)案例驗(yàn)證模型的可行性和結(jié)果的精確性。

    1 列車備選集的構(gòu)建及優(yōu)化

    1.1 列車備選集的構(gòu)建

    列車備選集是指技術(shù)站之間可開行的全部直達(dá)(直通)列車和區(qū)段列車的集合。以圖1所示的直線方向5個(gè)技術(shù)站為例,t為列車,假定兩相鄰技術(shù)站之間必開區(qū)段列車,即t12、t23、t34、t45;方向上所有可開行的直達(dá)(或直通)列車(以下統(tǒng)稱直達(dá)列車)為t13、t14、t15、t24、t25、t35,共6列。此時(shí)的列車備選集為{t12,t23,t34,t45,t13,t14,t15,t24,t25,t35}。同理,當(dāng)直線方向有n個(gè)技術(shù)站時(shí),備選列車數(shù)q為

    圖1 直線方向5個(gè)技術(shù)站

    式中:第一項(xiàng)為直達(dá)列車數(shù);第二項(xiàng)為區(qū)段列車數(shù)。

    對(duì)于直線方向的鐵路線路,相當(dāng)于已經(jīng)給定了任意兩個(gè)站之間的車流徑路,列車備選集的構(gòu)建簡單明了。簡單路網(wǎng)示意見圖2,以O(shè)D對(duì)(1,4)之間的車流為例,車流的走行徑路可能為1→2→4或1→2→3→4,這兩種車流徑路會(huì)生成不同的備選列車{t12,t24,t124}和{t12,t23,t34,t123,t234,t1234}。因此,需要先確定路網(wǎng)條件下每個(gè)OD對(duì)之間的最終車流徑路,從而將路網(wǎng)轉(zhuǎn)化為各OD對(duì)間只有唯一徑路的簡單直線線路,以此構(gòu)造路網(wǎng)條件下的列車備選集。

    圖2 簡單路網(wǎng)示意圖

    假設(shè)圖2路網(wǎng)中所有OD間的所有車流徑路都滿足繞道率且最終徑路已知,圖2所示路網(wǎng)(單方向)列車備選集的構(gòu)建過程見表1,剩余未列舉的OD對(duì)恰好與表1中列舉的對(duì)稱。

    表1 簡單路網(wǎng)單方向的列車備選集構(gòu)建

    需要注意,對(duì)于備選集里的列車,每支車流只能選擇由其車流徑路產(chǎn)生的列車。以O(shè)D對(duì)(1,3)之間的車流為例,該支車流只能選擇徑路1→2→3產(chǎn)生的列車{t12,t23,t123}。該規(guī)定能夠剔除大量冗余方案,進(jìn)一步減少變量規(guī)模。

    由于實(shí)際路網(wǎng)往往比圖2所示情況復(fù)雜,每個(gè)OD對(duì)之間可能存在成百上千條車流徑路,這些徑路一方面難以全部納入考慮范圍,另一方面包含了大量過長的無效徑路。為降低問題的復(fù)雜度,在求解各OD對(duì)間的備選車流徑路時(shí),引入車流徑路繞道率μr[14]來剔除過長的不合理徑路,定義為

    (1)

    式中:Lr為某OD間的任意一條徑路r的長度;Lrmin為該OD間最短徑路rmin的長度。

    當(dāng)某條徑路的繞道率大于規(guī)定閾值(本文取1.35)時(shí),該條徑路將不予以考慮。

    1.2 列車備選集的優(yōu)化

    上述列車備選集中的直達(dá)列車,是理論可開行的所有直達(dá)列車,而實(shí)際生產(chǎn)中,由于技術(shù)站作業(yè)設(shè)備、線路條件、改編時(shí)間消耗和車流量大小等條件的影響,兩技術(shù)站間可能必須開行直達(dá)列車或是禁止開行直達(dá)列車,以下稱為必開條件和禁開條件。

    (1)必開條件

    傳統(tǒng)的分析計(jì)算法指出,若某一單支車流無改編通過其運(yùn)行途中任何一個(gè)技術(shù)站所得車小時(shí)節(jié)省都大于或等于其在編成站的集結(jié)車小時(shí)消耗,則該車流滿足絕對(duì)條件,可直接列入方向最優(yōu)編組方案[1]。對(duì)應(yīng)列車備選集,若某支車流滿足絕對(duì)條件,則該車流起訖點(diǎn)之間必開直達(dá)列車。除此之外,規(guī)定兩相鄰技術(shù)站之間必開區(qū)段列車。

    (2)禁開條件

    編組計(jì)劃需遵守技術(shù)站禁開去向規(guī)定,通過剔除禁開去向的直達(dá)列車,可以進(jìn)一步縮減列車備選集。

    2 單組列車編組計(jì)劃兩階段優(yōu)化模型

    模型中的集合、參數(shù)及變量說明分別見表2~表4。

    表2 集合符號(hào)及定義

    表3 參數(shù)符號(hào)及定義

    表4 變量符號(hào)及定義

    2.1 第一階段——車流徑路模型

    (2)

    (3)

    (4)

    (5)

    (6)

    目標(biāo)函數(shù)式(2)第一項(xiàng)指車流選擇不同徑路時(shí)的總走行車小時(shí)消耗,第二項(xiàng)指車流完全改編車小時(shí)消耗,表示所有車流在其運(yùn)行途中的每個(gè)技術(shù)站均進(jìn)行改編時(shí)的車小時(shí)總和。由于車流總改編車小時(shí)消耗等于車流完全改編消耗減去車流無改編節(jié)省,而每支車流在每條徑路下的完全改編車小時(shí)都是確定的,為了使總的車流改編車小時(shí)盡可能少,車流徑路模型求解出的最終徑路既要使車流走行總消耗少,也要使車流完全改編消耗盡可能少,同時(shí),編組計(jì)劃模型目標(biāo)函數(shù)使車流無改編節(jié)省最大,兩個(gè)模型的目標(biāo)函數(shù)相加將能保證總的車流改編消耗最少。約束式(3)為徑路唯一性約束,表示每支車流最終只能選擇一條車流徑路。約束式(4)為弧通過能力約束,弧通過能力以“列”為單位,每條有向弧代表一個(gè)單向站間區(qū)間,因此,流經(jīng)任一有向弧的車流量不能超過其可用通過能力。

    (7)

    對(duì)于式(4)和式(7),模型任選一式即可。

    2.2 第二階段——編組計(jì)劃模型

    (8)

    (9)

    (10)

    (11)

    ?v∈V?d∈D

    (12)

    (13)

    xt=1 ?t∈Ssuff

    (14)

    xt=0 ?t∈Sforb

    (15)

    xt∈{0,1} ?t∈P

    (16)

    (17)

    wt∈N?t∈P

    (18)

    式中:wt為列車t的車流所占用的分類線數(shù),為非負(fù)整數(shù);k為每條分類線一晝夜可服務(wù)的貨車輛數(shù),取值為200車/條[5-6],并規(guī)定每條分類線只能由一個(gè)編組去向的車組占用;M為一個(gè)很大的正數(shù)。

    (19)

    則每個(gè)編組去向占用的分類線數(shù)量為

    (20)

    因此,技術(shù)站的分類線約束為

    (21)

    式(21)為非線性約束,故引入輔助決策變量wt,因此,式(21)可線性化為

    (22)

    根據(jù)式(20),wt滿足

    k(wt-1)

    (23)

    式(12)為流平衡約束,能夠保證兩技術(shù)站間只有一種列車開行方案以及每支車流只選擇一種輸送方案,同時(shí)還能保證車流從始發(fā)站運(yùn)輸至終到站;式(13)為變量的邏輯約束,只有在列車開行的情況下,車流才能選擇該列車;式(14)、式(15)分別表示滿足必開和禁開條件時(shí)列車開行決策變量的取值;式(16)~式(18)分別表示列車開行決策變量、車流分配決策變量、輔助決策變量的取值范圍。

    經(jīng)線性化處理后,模型整體為線性規(guī)劃,可以通過功能強(qiáng)大的商業(yè)求解器或者傳統(tǒng)精確算法求出精確解。模型的求解流程見圖3。

    圖3 模型整體求解流程

    2.3 模型復(fù)雜度分析

    表3 模型規(guī)模對(duì)比

    3 案例求解

    分別以INFORMS 2019 RAS競賽[14]中的8個(gè)技術(shù)站和32個(gè)技術(shù)站的路網(wǎng)對(duì)本文所提方法和模型進(jìn)行求解驗(yàn)證。

    3.1 小規(guī)模路網(wǎng)案例

    8個(gè)技術(shù)站的小規(guī)模路網(wǎng)見圖4,路網(wǎng)數(shù)據(jù)及參數(shù)均按照文獻(xiàn)[14]要求設(shè)置,不考慮必開去向和禁開去向。在處理器為Inter(R) Core(TM) i7-9750 H 2.60 GHz,內(nèi)存為16 GB,操作系統(tǒng)為Windows 10的Lenovo ThinkPad P53計(jì)算機(jī)上,用Python 3.7語言編程實(shí)現(xiàn),并運(yùn)用新一代大規(guī)模數(shù)學(xué)規(guī)劃優(yōu)化器Gurobi 9.1.0對(duì)提出的模型進(jìn)行求解驗(yàn)證。

    圖4 小規(guī)模路網(wǎng)示意圖

    求解得到車流總的消耗為236 135.1車·h,最終開行34列列車,其中區(qū)段列車18列,直達(dá)列車16列。各技術(shù)站的改編能力及分類線使用情況見表4,各區(qū)間的通過能力利用情況見表5,最終的直達(dá)列車服務(wù)方案見圖5,最終的車流分配方案見表6。

    圖5 最終的直達(dá)列車服務(wù)方案

    表4 各技術(shù)站的改編能力及分類線使用情況

    表5 各區(qū)間通過能力利用情況

    表6 最終的車流分配方案

    文獻(xiàn)[14]求解結(jié)果與本文兩階段線性規(guī)劃模型的求解結(jié)果對(duì)比見表7。其中,文獻(xiàn)[14]基于模擬退火算法求解非線性綜合優(yōu)化模型,而本文采用精確方法求解線性模型,兩者求解結(jié)果相比,本文總消耗減少了488.9車·h,結(jié)果表明本文所提方法和模型有效。

    表7 小規(guī)模路網(wǎng)的求解結(jié)果對(duì)比 車·h

    3.2 大規(guī)模路網(wǎng)案例

    8個(gè)技術(shù)站的小規(guī)模案例已證明本文所提模型的可行性,為了更好驗(yàn)證和展示模型的求解效率和精確性,以文獻(xiàn)[14]中32個(gè)站的大規(guī)模路網(wǎng)為例,路網(wǎng)共包括32個(gè)技術(shù)站,128個(gè)單向站間區(qū)間,992支車流,見圖6。

    圖6 大規(guī)模路網(wǎng)示意圖

    由于該路網(wǎng)各技術(shù)站之間連接的復(fù)雜性和區(qū)間長度的特點(diǎn),導(dǎo)致兩技術(shù)站間滿足繞道率μ=1.35的備選徑路可能存在成百上千條,為了降低求解的復(fù)雜度,在滿足繞道率的條件下,任意兩技術(shù)站間的備選徑路最多取20條,設(shè)置求解GAP為0.1%,其余參數(shù)不變,模型共耗時(shí)91 s得到最優(yōu)解,車流總消耗為5 163 471.8車·h,較文獻(xiàn)[14]所給結(jié)果共節(jié)省1 855車·h。

    為研究模型的穩(wěn)定性以及各OD對(duì)間備選徑路數(shù)量對(duì)求解結(jié)果的影響,本文在求解車流徑路模型時(shí),每個(gè)OD對(duì)間滿足繞道率的K-短路最多分別選5、10、20、30、40、50條。求解結(jié)果見表8,模型求解時(shí)間在91~248 s之間,平均耗時(shí)147 s,證明本文模型求解效率高且穩(wěn)定;6次實(shí)驗(yàn)總消耗最大相差2 953.6車·h,平均總消耗為5 164 922.3車·h,較文獻(xiàn)[14]平均節(jié)省404車·h,證明備選徑路數(shù)量對(duì)本文模型的求解結(jié)果影響較小,同時(shí)也說明,在利用本文模型求解單組列車編組計(jì)劃時(shí),只需為每個(gè)OD對(duì)選取少量的備選徑路便可得到近似最優(yōu)解。

    表8 各OD對(duì)間不同備選徑路數(shù)對(duì)求解結(jié)果的影響

    4 結(jié)論

    技術(shù)站列車編組計(jì)劃作為貨物列車編組計(jì)劃的重要組成部分,規(guī)定了路網(wǎng)上所有技術(shù)站間的車流徑路、列車的編成、到達(dá)站(或是裝卸、解體站)以及列車的編組內(nèi)容和編掛辦法等,是車流組織計(jì)劃,也是路網(wǎng)上各車站分工的戰(zhàn)略部署[1],對(duì)發(fā)揮鐵路運(yùn)輸能力,提高運(yùn)輸效率和經(jīng)濟(jì)效益有著重要作用,準(zhǔn)確、高效地求解出符合運(yùn)輸需求的技術(shù)站列車編組計(jì)劃對(duì)運(yùn)輸生產(chǎn)尤為關(guān)鍵。

    本文基于多商品流服務(wù)網(wǎng)絡(luò)設(shè)計(jì)原理,提出了列車備選集的構(gòu)建方法,并建立了基于列車備選集的技術(shù)站單組列車編組計(jì)劃兩階段線性規(guī)劃模型,得到以下結(jié)論:①較車流遞歸方法,本文模型更加簡單,能夠直觀展示車流的改編過程;②構(gòu)建列車備選集可以限制每支車流只能分配到其最優(yōu)徑路生成的備選列車上,以此剔除不合理的車流分配變量,減小模型規(guī)模,加快求解效率;③將分類線約束有效轉(zhuǎn)化為線性約束,既避免了復(fù)雜的非線性結(jié)構(gòu),也避免將分類線約束轉(zhuǎn)化為最大編組去向數(shù)而造成與運(yùn)輸實(shí)際不符;④小規(guī)模和大規(guī)模路網(wǎng)的求解結(jié)果表明,本文模型能夠快速、準(zhǔn)確地得到最優(yōu)解,求解結(jié)果質(zhì)量優(yōu)于非線性模型;⑤各OD對(duì)間不同備選徑路數(shù)對(duì)本文模型求解結(jié)果影響很小,模型穩(wěn)定可靠。下一階段的研究將著重于車流徑路和編組計(jì)劃的綜合協(xié)同優(yōu)化,同時(shí)將分組列車納入考慮,構(gòu)建更加符合運(yùn)輸實(shí)際的編組計(jì)劃模型。

    猜你喜歡
    優(yōu)化模型
    一半模型
    超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
    民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數(shù)”優(yōu)化運(yùn)算——以2021年解析幾何高考題為例
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    欧美黄色淫秽网站| 久久久久亚洲av毛片大全| 男女做爰动态图高潮gif福利片| 亚洲自偷自拍三级| 午夜福利在线观看吧| 日本黄色片子视频| 看片在线看免费视频| 亚洲国产高清在线一区二区三| 欧美bdsm另类| 欧美激情在线99| 亚洲一区高清亚洲精品| 亚洲国产日韩欧美精品在线观看| 成人国产一区最新在线观看| 亚洲,欧美精品.| 免费大片18禁| 日韩高清综合在线| 级片在线观看| 日本黄色视频三级网站网址| 亚洲综合色惰| 亚洲欧美清纯卡通| 成人美女网站在线观看视频| 国内精品美女久久久久久| 国产主播在线观看一区二区| 精品免费久久久久久久清纯| 免费在线观看亚洲国产| 亚洲中文日韩欧美视频| 精品人妻偷拍中文字幕| 九九热线精品视视频播放| 国产欧美日韩一区二区精品| 啦啦啦观看免费观看视频高清| 99视频精品全部免费 在线| 男人舔女人下体高潮全视频| 欧美激情在线99| 亚洲第一欧美日韩一区二区三区| 人妻制服诱惑在线中文字幕| 美女大奶头视频| av在线天堂中文字幕| 熟女人妻精品中文字幕| 亚洲av中文字字幕乱码综合| 神马国产精品三级电影在线观看| 国产亚洲欧美在线一区二区| 欧美色视频一区免费| 久久人妻av系列| 久久国产精品影院| 国产精品免费一区二区三区在线| 亚洲最大成人手机在线| 欧美最黄视频在线播放免费| 欧美+亚洲+日韩+国产| 成年女人永久免费观看视频| 日本撒尿小便嘘嘘汇集6| 午夜视频国产福利| 成人毛片a级毛片在线播放| 国产欧美日韩一区二区三| 99久久无色码亚洲精品果冻| 直男gayav资源| 99久久九九国产精品国产免费| 美女免费视频网站| 夜夜爽天天搞| 国产又黄又爽又无遮挡在线| 色哟哟哟哟哟哟| 深夜精品福利| 狂野欧美白嫩少妇大欣赏| 脱女人内裤的视频| 一区二区三区高清视频在线| 丰满乱子伦码专区| 麻豆国产av国片精品| 丰满人妻熟妇乱又伦精品不卡| 国产精品98久久久久久宅男小说| 欧美黄色片欧美黄色片| 精品熟女少妇八av免费久了| av专区在线播放| 久久久久久九九精品二区国产| АⅤ资源中文在线天堂| 免费在线观看成人毛片| 国产精品亚洲美女久久久| av天堂在线播放| 中文字幕熟女人妻在线| 免费在线观看亚洲国产| 成人国产综合亚洲| 香蕉av资源在线| 97热精品久久久久久| 丰满人妻熟妇乱又伦精品不卡| 亚洲专区国产一区二区| 国产成年人精品一区二区| 毛片女人毛片| 少妇的逼水好多| 91午夜精品亚洲一区二区三区 | 亚洲精品影视一区二区三区av| 99热这里只有精品一区| 日本成人三级电影网站| 国产国拍精品亚洲av在线观看| 在现免费观看毛片| 18+在线观看网站| av在线观看视频网站免费| 极品教师在线免费播放| 免费黄网站久久成人精品 | 毛片女人毛片| 国产精品人妻久久久久久| av在线观看视频网站免费| 在线观看美女被高潮喷水网站 | 成人国产一区最新在线观看| 一区二区三区四区激情视频 | 成年女人看的毛片在线观看| 少妇的逼好多水| 亚洲人成电影免费在线| 真人做人爱边吃奶动态| 亚洲乱码一区二区免费版| 嫩草影院精品99| 国产日本99.免费观看| 在线播放无遮挡| 午夜影院日韩av| 俄罗斯特黄特色一大片| 高清在线国产一区| 色哟哟哟哟哟哟| 亚洲av日韩精品久久久久久密| 三级毛片av免费| 国产亚洲精品av在线| 亚洲精品日韩av片在线观看| 91在线精品国自产拍蜜月| 99久国产av精品| 久久久久国产精品人妻aⅴ院| 国产主播在线观看一区二区| 真人一进一出gif抽搐免费| 欧美极品一区二区三区四区| 国产精品99久久久久久久久| 久久婷婷人人爽人人干人人爱| 色哟哟·www| 久久性视频一级片| 国产激情偷乱视频一区二区| 午夜福利在线在线| 色精品久久人妻99蜜桃| 国产色婷婷99| 91av网一区二区| 亚洲精品影视一区二区三区av| 午夜久久久久精精品| 永久网站在线| 琪琪午夜伦伦电影理论片6080| 在线看三级毛片| 日韩国内少妇激情av| 亚洲美女视频黄频| 国产欧美日韩一区二区精品| 小蜜桃在线观看免费完整版高清| 最新中文字幕久久久久| 精品熟女少妇八av免费久了| 白带黄色成豆腐渣| 亚洲精品粉嫩美女一区| 欧美极品一区二区三区四区| 日韩人妻高清精品专区| 嫩草影院新地址| 亚洲男人的天堂狠狠| 中文亚洲av片在线观看爽| 亚洲内射少妇av| 久久久久性生活片| 久久久久国产精品人妻aⅴ院| 高清日韩中文字幕在线| 亚洲人成伊人成综合网2020| 又紧又爽又黄一区二区| 亚洲最大成人av| 成人亚洲精品av一区二区| 亚洲人成网站在线播放欧美日韩| 噜噜噜噜噜久久久久久91| 欧美中文日本在线观看视频| 亚洲成人免费电影在线观看| 一级毛片久久久久久久久女| 人人妻,人人澡人人爽秒播| 国产精品电影一区二区三区| eeuss影院久久| 日韩av在线大香蕉| 可以在线观看毛片的网站| 亚洲真实伦在线观看| 91在线精品国自产拍蜜月| 尤物成人国产欧美一区二区三区| 身体一侧抽搐| 男人舔女人下体高潮全视频| 久久国产乱子免费精品| 变态另类丝袜制服| 亚洲一区高清亚洲精品| 精品一区二区三区视频在线观看免费| a级毛片a级免费在线| 极品教师在线视频| 国产aⅴ精品一区二区三区波| 精品人妻视频免费看| 日韩中字成人| 国产欧美日韩一区二区精品| 欧美日韩国产亚洲二区| 丰满的人妻完整版| 99热这里只有是精品在线观看 | 长腿黑丝高跟| 内地一区二区视频在线| 男人舔奶头视频| 欧美在线一区亚洲| 国产探花在线观看一区二区| 久久久久亚洲av毛片大全| 欧美在线一区亚洲| 亚洲成av人片在线播放无| 国产主播在线观看一区二区| xxxwww97欧美| 不卡一级毛片| 岛国在线免费视频观看| 免费一级毛片在线播放高清视频| 日本 av在线| 国产淫片久久久久久久久 | www.999成人在线观看| 国产日本99.免费观看| 国产精品98久久久久久宅男小说| 欧美性感艳星| 可以在线观看毛片的网站| eeuss影院久久| 真人一进一出gif抽搐免费| 欧美乱色亚洲激情| 天堂动漫精品| 国模一区二区三区四区视频| 国产精品久久久久久久久免 | 久久久精品欧美日韩精品| 午夜精品在线福利| 网址你懂的国产日韩在线| 一二三四社区在线视频社区8| 国产69精品久久久久777片| 九九在线视频观看精品| 亚洲av第一区精品v没综合| 国产精品永久免费网站| 人人妻,人人澡人人爽秒播| 日韩av在线大香蕉| 国产午夜精品论理片| 亚洲综合色惰| 色综合欧美亚洲国产小说| 亚洲第一区二区三区不卡| 国产伦人伦偷精品视频| 午夜视频国产福利| 亚洲18禁久久av| 国产69精品久久久久777片| 日韩中文字幕欧美一区二区| 超碰av人人做人人爽久久| 久久国产精品人妻蜜桃| av欧美777| 一进一出抽搐动态| 亚洲第一区二区三区不卡| 午夜亚洲福利在线播放| 赤兔流量卡办理| 国产精品久久久久久精品电影| 美女 人体艺术 gogo| 亚洲精品一卡2卡三卡4卡5卡| 90打野战视频偷拍视频| 午夜精品一区二区三区免费看| 搡老岳熟女国产| 成人av在线播放网站| 黄色一级大片看看| 欧美一区二区国产精品久久精品| 国产久久久一区二区三区| 97碰自拍视频| 久久九九热精品免费| 久久久久久久久久成人| 综合色av麻豆| 亚洲三级黄色毛片| 亚洲综合色惰| 欧美高清性xxxxhd video| 亚洲 国产 在线| 一区二区三区免费毛片| 99国产精品一区二区三区| 精品一区二区三区人妻视频| 中亚洲国语对白在线视频| 在线观看舔阴道视频| 五月玫瑰六月丁香| 亚洲精品乱码久久久v下载方式| 男人舔奶头视频| 一本久久中文字幕| 国产精品1区2区在线观看.| 久久久久久久午夜电影| 淫妇啪啪啪对白视频| 日韩中字成人| 国产单亲对白刺激| 日韩精品青青久久久久久| 欧美成人一区二区免费高清观看| 757午夜福利合集在线观看| 国产一区二区亚洲精品在线观看| 人妻丰满熟妇av一区二区三区| 少妇被粗大猛烈的视频| 亚洲精品456在线播放app | 国产精品日韩av在线免费观看| 久久伊人香网站| 午夜影院日韩av| av女优亚洲男人天堂| 国产精品98久久久久久宅男小说| 亚洲aⅴ乱码一区二区在线播放| 国产激情偷乱视频一区二区| 精品乱码久久久久久99久播| 精品99又大又爽又粗少妇毛片 | 偷拍熟女少妇极品色| 欧美一级a爱片免费观看看| 最新中文字幕久久久久| 老司机福利观看| 久久欧美精品欧美久久欧美| 精品熟女少妇八av免费久了| 两个人的视频大全免费| 美女cb高潮喷水在线观看| 757午夜福利合集在线观看| or卡值多少钱| 波多野结衣巨乳人妻| 亚洲,欧美精品.| 精品久久久久久成人av| 美女免费视频网站| 国产精品久久久久久久久免 | 久9热在线精品视频| 国产精品亚洲一级av第二区| 色精品久久人妻99蜜桃| .国产精品久久| 搡老妇女老女人老熟妇| 午夜免费男女啪啪视频观看 | 久久久成人免费电影| 国产精品久久久久久久久免 | 久久久国产成人免费| 变态另类成人亚洲欧美熟女| 国产精品久久视频播放| h日本视频在线播放| 嫩草影院新地址| 国产亚洲精品av在线| 亚洲内射少妇av| av在线老鸭窝| 88av欧美| 国产在线精品亚洲第一网站| 亚洲精品色激情综合| 丰满人妻熟妇乱又伦精品不卡| 欧美在线黄色| 人妻久久中文字幕网| 国产三级在线视频| 国产日本99.免费观看| 中文字幕av成人在线电影| 1024手机看黄色片| 身体一侧抽搐| 欧美色视频一区免费| 亚洲成人免费电影在线观看| 久久天躁狠狠躁夜夜2o2o| 天天一区二区日本电影三级| 深爱激情五月婷婷| 搡老熟女国产l中国老女人| 国产av不卡久久| 亚洲欧美激情综合另类| 国内毛片毛片毛片毛片毛片| 国产精品不卡视频一区二区 | 国产又黄又爽又无遮挡在线| 99热只有精品国产| 日本在线视频免费播放| 日本 av在线| 亚洲男人的天堂狠狠| 搡老岳熟女国产| 成人鲁丝片一二三区免费| 亚洲精品亚洲一区二区| 免费在线观看亚洲国产| 午夜福利18| 99热6这里只有精品| 一进一出抽搐gif免费好疼| 成年人黄色毛片网站| 在线天堂最新版资源| 18美女黄网站色大片免费观看| 黄色一级大片看看| 精品久久久久久久末码| 亚洲色图av天堂| 亚洲欧美日韩卡通动漫| 激情在线观看视频在线高清| 中亚洲国语对白在线视频| 国产白丝娇喘喷水9色精品| 日日干狠狠操夜夜爽| 亚洲欧美日韩东京热| 一级毛片久久久久久久久女| 成人高潮视频无遮挡免费网站| 免费大片18禁| 欧美精品啪啪一区二区三区| 日本 欧美在线| 欧美最新免费一区二区三区 | 天天一区二区日本电影三级| 免费人成视频x8x8入口观看| 国产老妇女一区| 69av精品久久久久久| 日韩国内少妇激情av| 在线观看美女被高潮喷水网站 | 午夜福利免费观看在线| 精品久久久久久久人妻蜜臀av| 久久精品影院6| 免费在线观看日本一区| 国产欧美日韩精品亚洲av| 国产探花在线观看一区二区| 有码 亚洲区| 噜噜噜噜噜久久久久久91| xxxwww97欧美| 日韩欧美精品v在线| 成人鲁丝片一二三区免费| 少妇人妻一区二区三区视频| 亚洲av不卡在线观看| 精华霜和精华液先用哪个| 久久99热6这里只有精品| 国产精品影院久久| 国产单亲对白刺激| 久久这里只有精品中国| 国产在视频线在精品| 国产精品久久久久久精品电影| 高清在线国产一区| 97超视频在线观看视频| 亚洲在线观看片| 国产高清视频在线观看网站| 成人精品一区二区免费| or卡值多少钱| 天天一区二区日本电影三级| 中亚洲国语对白在线视频| 亚洲欧美日韩东京热| 女人十人毛片免费观看3o分钟| 激情在线观看视频在线高清| 俺也久久电影网| 婷婷精品国产亚洲av| 亚洲精品久久国产高清桃花| 亚洲中文字幕一区二区三区有码在线看| 日韩欧美精品免费久久 | 久99久视频精品免费| 熟女电影av网| 国产亚洲精品久久久com| 国产一级毛片七仙女欲春2| 在现免费观看毛片| 神马国产精品三级电影在线观看| 蜜桃亚洲精品一区二区三区| 久久婷婷人人爽人人干人人爱| 国产aⅴ精品一区二区三区波| 国产成年人精品一区二区| 欧美潮喷喷水| 舔av片在线| 一本综合久久免费| 99久久无色码亚洲精品果冻| 一区二区三区激情视频| bbb黄色大片| 人人妻人人澡欧美一区二区| 中文字幕免费在线视频6| 国产伦精品一区二区三区视频9| 黄色丝袜av网址大全| 他把我摸到了高潮在线观看| 免费人成在线观看视频色| a级毛片免费高清观看在线播放| 色5月婷婷丁香| 一个人免费在线观看电影| 亚洲国产精品合色在线| 99riav亚洲国产免费| 亚洲成人久久性| 久久久国产成人免费| 久久久久亚洲av毛片大全| 蜜桃亚洲精品一区二区三区| 亚洲av.av天堂| 天堂网av新在线| 久久久久久大精品| 人妻制服诱惑在线中文字幕| 午夜免费激情av| 久久热精品热| 久久久国产成人精品二区| 丁香欧美五月| 一本精品99久久精品77| 变态另类丝袜制服| 国产免费男女视频| 色精品久久人妻99蜜桃| 国产高潮美女av| 美女高潮的动态| 亚洲美女视频黄频| 成人美女网站在线观看视频| 国产精品精品国产色婷婷| 少妇熟女aⅴ在线视频| 亚洲美女视频黄频| 久久久精品大字幕| 久久久精品欧美日韩精品| 黄色配什么色好看| 久久午夜福利片| 亚洲一区二区三区色噜噜| 三级毛片av免费| 深夜a级毛片| 精品久久久久久久久久久久久| 男女做爰动态图高潮gif福利片| 国产色爽女视频免费观看| 久久久久久国产a免费观看| 国产精品亚洲美女久久久| 最近最新免费中文字幕在线| 动漫黄色视频在线观看| 免费在线观看成人毛片| 最新中文字幕久久久久| 国产探花在线观看一区二区| 国产伦人伦偷精品视频| 亚洲国产高清在线一区二区三| 欧美xxxx性猛交bbbb| 五月伊人婷婷丁香| 久久久精品大字幕| 日韩欧美免费精品| 在线天堂最新版资源| 午夜激情欧美在线| 中文字幕人成人乱码亚洲影| 男女那种视频在线观看| 欧美性猛交╳xxx乱大交人| 国产黄片美女视频| 99在线视频只有这里精品首页| 精品人妻1区二区| 蜜桃久久精品国产亚洲av| 亚州av有码| 国产熟女xx| 欧美成人一区二区免费高清观看| 日本免费一区二区三区高清不卡| 最近最新中文字幕大全电影3| 成年女人永久免费观看视频| 国产毛片a区久久久久| 日韩欧美精品免费久久 | 日本撒尿小便嘘嘘汇集6| 热99在线观看视频| 欧美最黄视频在线播放免费| 日日干狠狠操夜夜爽| 999久久久精品免费观看国产| 非洲黑人性xxxx精品又粗又长| 久久精品久久久久久噜噜老黄 | 国产高潮美女av| 国产精品久久久久久久久免 | 亚洲不卡免费看| 美女高潮喷水抽搐中文字幕| 国产成+人综合+亚洲专区| 日韩欧美精品免费久久 | 亚洲不卡免费看| 床上黄色一级片| 欧美日韩瑟瑟在线播放| 精品欧美国产一区二区三| 亚洲欧美日韩无卡精品| 国产 一区 欧美 日韩| 成人特级av手机在线观看| 久久久久国产精品人妻aⅴ院| 精品欧美国产一区二区三| 日韩成人在线观看一区二区三区| 国产精品一区二区免费欧美| 极品教师在线免费播放| 中文亚洲av片在线观看爽| 欧美色欧美亚洲另类二区| 国产欧美日韩一区二区三| 免费av观看视频| 国模一区二区三区四区视频| 午夜免费成人在线视频| 国产综合懂色| 听说在线观看完整版免费高清| 国产精品1区2区在线观看.| 欧美xxxx性猛交bbbb| 国产精品av视频在线免费观看| 99热精品在线国产| 免费人成视频x8x8入口观看| 久久久久久九九精品二区国产| 中文字幕免费在线视频6| 色哟哟·www| 成年版毛片免费区| 热99re8久久精品国产| av欧美777| 真人做人爱边吃奶动态| 夜夜爽天天搞| 欧美又色又爽又黄视频| 黄色丝袜av网址大全| 国产精品精品国产色婷婷| 亚洲中文字幕日韩| 国产精品伦人一区二区| 99在线人妻在线中文字幕| 欧美精品啪啪一区二区三区| 简卡轻食公司| 99久久精品热视频| 中文资源天堂在线| 午夜免费激情av| 国产真实伦视频高清在线观看 | 国内精品美女久久久久久| 91麻豆av在线| 天堂√8在线中文| 成年版毛片免费区| 一级黄片播放器| 91麻豆精品激情在线观看国产| 欧美日韩国产亚洲二区| 一本一本综合久久| 久久精品久久久久久噜噜老黄 | 日本五十路高清| 最近最新中文字幕大全电影3| 日本免费a在线| 国产极品精品免费视频能看的| 极品教师在线视频| 国产三级黄色录像| 国产老妇女一区| 性欧美人与动物交配| 黄片小视频在线播放| 熟妇人妻久久中文字幕3abv| 亚洲最大成人手机在线| 午夜福利欧美成人| 中文字幕av在线有码专区| 18美女黄网站色大片免费观看| 亚洲乱码一区二区免费版| 高潮久久久久久久久久久不卡| 亚洲av熟女| 欧美+亚洲+日韩+国产| 99久久成人亚洲精品观看| 超碰av人人做人人爽久久| 99久久久亚洲精品蜜臀av| 中文字幕免费在线视频6| 露出奶头的视频| 国产精品久久久久久久电影| 亚洲第一欧美日韩一区二区三区| 亚洲国产欧美人成| 午夜两性在线视频| 人妻久久中文字幕网| a级毛片a级免费在线| 日韩中文字幕欧美一区二区| 精品日产1卡2卡| 久久久久久久午夜电影| 成人特级黄色片久久久久久久| 久久草成人影院| 久久久久久久午夜电影| 亚洲国产高清在线一区二区三| 老司机福利观看| 亚洲最大成人av| 三级国产精品欧美在线观看| 国产成人啪精品午夜网站| 欧美一区二区亚洲| 中文亚洲av片在线观看爽| 成熟少妇高潮喷水视频| 噜噜噜噜噜久久久久久91| 色哟哟哟哟哟哟| 久久伊人香网站| 日本免费a在线|