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

    基于換乘導(dǎo)向的大型客運(yùn)樞紐高鐵列車(chē)接續(xù)優(yōu)化

    2019-10-31 09:21:33喬俊孟學(xué)雷王東先湯霖
    計(jì)算機(jī)應(yīng)用 2019年9期

    喬俊 孟學(xué)雷 王東先 湯霖

    摘 要:針對(duì)高速鐵路成網(wǎng)條件下的客運(yùn)樞紐高鐵列車(chē)接續(xù)優(yōu)化問(wèn)題,分析了樞紐內(nèi)的旅客換乘過(guò)程,提出了中長(zhǎng)途客流的換乘滿(mǎn)意度概念;以平均換乘滿(mǎn)意度和樞紐車(chē)站列車(chē)到發(fā)均衡性為優(yōu)化目標(biāo),以大站合理發(fā)車(chē)時(shí)間、合理終到時(shí)間、車(chē)站作業(yè)間隔時(shí)間、旅客換乘時(shí)間、車(chē)站到發(fā)線(xiàn)能力等為約束條件,建立了基于換乘協(xié)同的大型客運(yùn)樞紐高速列車(chē)接續(xù)優(yōu)化模型。設(shè)計(jì)了改進(jìn)染色體編碼方式和選擇策略的遺傳算法對(duì)算例進(jìn)行了求解。改進(jìn)后的遺傳算法同基本遺傳算法、基本模擬退火算法相比,目標(biāo)函數(shù)中所求的平均換乘滿(mǎn)意度分別增加了5.10%、2.93%,樞紐車(chē)站列車(chē)到發(fā)均衡性分別提高了0.27%、2.31%,算例結(jié)果驗(yàn)證了改進(jìn)遺傳算法的有效性和穩(wěn)定性,表明所提方法可以有效地提高大型樞紐高鐵列車(chē)的接續(xù)質(zhì)量。

    關(guān)鍵詞:換乘滿(mǎn)意度;同站換乘;異站換乘;列車(chē)接續(xù);改進(jìn)遺傳算法

    中圖分類(lèi)號(hào):TP301.6

    文獻(xiàn)標(biāo)志碼:A

    High-speed train connection optimization for large passenger transport hub based on transfer orientation

    QIAO Jun1, MENG Xuelei1*, WANG Dongxian1, TANG Lin2

    1.School of Traffic and Transportation, Lanzhou Jiaotong University, Lanzhou Gansu 730070, China;

    2.Lanzhou Railway Station, China Railway Lanzhou Bureau Group Company Limited, Lanzhou Gansu 730070, China

    Abstract:

    In view of the optimization of high-speed train connection in passenger transport hub under the condition of high-speed railway network, the concept of transfer satisfaction of medium and long distance passenger flow was proposed by analyzing the passenger transfer process in hub, and a high-speed train connection optimization model for large passenger transport hub based on transfer orientation was proposed with the average transfer satisfaction and the arrival and departure equilibrium of trains at hub stations as the optimization objective and with the constraint conditions of reasonable originating time of large stations, reasonable terminating time, station operation interval time, passenger transfer time and station arrival and departure line capacity. A genetic algorithm with improved chromosome coding mode and selection strategy was designed to solve the example. Compared with the basic genetic algorithm and the basic simulated annealing algorithm, the improved genetic algorithm increases the average transfer satisfaction in the objective function by 5.10% and 2.93% respectively, and raises the equilibrium of arrival and departure of trains at hub stations by 0.27% and 2.31% respectively. The results of the example verify the effectiveness and stability of the improved genetic algorithm, which indicates that the proposed method can effectively optimize the quality of the high-speed train connection in large passenger transport hub.

    Key words:

    transfer satisfaction; transfer inside one station; transfer between different stations; train connection; improved genetic algorithm

    0 引言

    我國(guó)高速鐵路的快速發(fā)展使得高鐵客流量和換乘客流量逐年攀升,客流種類(lèi)和換乘方式也越來(lái)越多樣化。大型高速鐵路客運(yùn)樞紐作為城市對(duì)內(nèi)對(duì)外交通銜接的重要節(jié)點(diǎn),樞紐內(nèi)的換乘問(wèn)題受到越來(lái)越多的重視??茖W(xué)地分析樞紐的換乘條件,制定合理的換乘接續(xù)方案,優(yōu)化樞紐內(nèi)換乘接續(xù)列車(chē)到發(fā)時(shí)刻的銜接問(wèn)題,對(duì)于保障客流高效地?fù)Q乘具有重要意義。

    近年來(lái),大型鐵路客運(yùn)樞紐的換乘協(xié)同問(wèn)題受到了不少?lài)?guó)內(nèi)外學(xué)者的關(guān)注。國(guó)外學(xué)者多通過(guò)調(diào)整時(shí)刻表的方法研究列車(chē)接續(xù)優(yōu)化問(wèn)題,文獻(xiàn)[1]研究了鐵路與城市公交間的時(shí)刻協(xié)同問(wèn)題,以換乘時(shí)間和公交運(yùn)營(yíng)成本最小為目標(biāo),建立了非線(xiàn)性模型并利用遺傳算法求解。文獻(xiàn)[2]構(gòu)建了鐵路網(wǎng)絡(luò)時(shí)刻表,以換乘時(shí)間最小為目標(biāo)建立了非線(xiàn)性時(shí)刻表優(yōu)化模型,利用遺傳算法求解模型,但模型中沒(méi)有考慮車(chē)站到發(fā)線(xiàn)能力約束。文獻(xiàn)[3]提出通過(guò)調(diào)整鐵路非周期性時(shí)刻表減少乘客換乘時(shí)間,以列車(chē)運(yùn)行時(shí)間、停站時(shí)間等為變量,乘客總的換乘等待時(shí)間最小為目標(biāo)建立了混合整數(shù)規(guī)劃模型,利用改進(jìn)的啟發(fā)式算法求解了模型,但只對(duì)同站換乘進(jìn)行了研究。文獻(xiàn)[4]研究了城市軌道交通換乘可靠性,通過(guò)數(shù)學(xué)推導(dǎo)分析了發(fā)車(chē)間隔時(shí)間、運(yùn)行時(shí)間對(duì)換乘等待時(shí)間的影響,但并未建立模型對(duì)發(fā)車(chē)間隔時(shí)間和運(yùn)行時(shí)間進(jìn)行優(yōu)化。文獻(xiàn)[5]研究了鐵路周期性時(shí)刻表下的旅客換乘問(wèn)題,以換乘旅客的時(shí)間損失最少為目標(biāo),建立了周期性時(shí)刻表優(yōu)化模型,并通過(guò)啟發(fā)式算法求解了模型,但模型并未考慮不同方向列車(chē)之間的需求換乘客流量。文獻(xiàn)[6]研究了以換乘需求為導(dǎo)向的軌道交通換乘優(yōu)化問(wèn)題,以乘客換乘等待時(shí)間最小為目標(biāo),構(gòu)建了線(xiàn)性和改進(jìn)的非線(xiàn)性?xún)?yōu)化模型,并利用粒子群算法求解了改進(jìn)的非線(xiàn)性模型。

    目前國(guó)內(nèi)對(duì)鐵路客運(yùn)樞紐旅客換乘和列車(chē)接續(xù)優(yōu)化問(wèn)題的研究相對(duì)較少。文獻(xiàn)[7]研究了城市軌道交通的換乘,以各條線(xiàn)路首列車(chē)的發(fā)車(chē)時(shí)刻為變量,乘客總的換乘等待時(shí)間最小為目標(biāo),建立了城市軌道交通網(wǎng)列車(chē)時(shí)刻表優(yōu)化模型,運(yùn)用CPLEX工具進(jìn)行了求解,但并未考慮需求換乘客流量約束以及發(fā)車(chē)間隔較大線(xiàn)路間的列車(chē)接續(xù)問(wèn)題。文獻(xiàn)[8]研究了鐵路樞紐與城市公交間的接駁問(wèn)題,以乘客總的換乘時(shí)間最小為目標(biāo),公交發(fā)車(chē)間隔、載客量等為約束建立了樞紐換乘時(shí)間調(diào)度模型。文獻(xiàn)[9]分析了鐵路樞紐內(nèi)的同站和異站換乘過(guò)程,提出了換乘條件評(píng)價(jià)指標(biāo),但并未利用這些評(píng)價(jià)指標(biāo)對(duì)樞紐列車(chē)接續(xù)優(yōu)化進(jìn)行深入研究。文獻(xiàn)[10]分析了歐洲鐵路旅客換乘組織和列車(chē)銜接方案,提出我國(guó)鐵路應(yīng)當(dāng)將列車(chē)銜接方案與開(kāi)行方案、運(yùn)行圖的編制過(guò)程結(jié)合起來(lái),為優(yōu)化樞紐列車(chē)接續(xù)問(wèn)題提供了解決思路。文獻(xiàn)[11]以有效換乘客流量最大和換乘無(wú)效等待時(shí)間最短為目標(biāo)建立了列車(chē)接續(xù)優(yōu)化模型,但模型沒(méi)有考慮樞紐車(chē)站列車(chē)的到發(fā)均衡性和異站換乘。總的來(lái)說(shuō),上述文獻(xiàn)[1-3,6-8,11]研究目標(biāo)多為換乘等待時(shí)間最小,約束條件[1-3,5,8,11]多為乘客必要換乘時(shí)間、列車(chē)停站時(shí)間、發(fā)車(chē)間隔時(shí)間等,且大多通過(guò)調(diào)整時(shí)刻表[1-3,5,7,11]優(yōu)化了列車(chē)接續(xù)問(wèn)題,這些研究成果為本文提供了重要參考,但研究問(wèn)題多限于樞紐內(nèi)單一客運(yùn)站不同交通方式之間的換乘和軌道交通換乘,異站換乘的研究較少,且算法設(shè)計(jì)時(shí)大多未將車(chē)次與到發(fā)時(shí)刻相結(jié)合。鮮有文獻(xiàn)在已知開(kāi)行方案的基礎(chǔ)上,通過(guò)優(yōu)化與編排一定周期內(nèi)換乘樞紐車(chē)站的所有列車(chē)到發(fā)順序和時(shí)間的方法來(lái)研究列車(chē)接續(xù)問(wèn)題。

    針對(duì)這些問(wèn)題,對(duì)不同線(xiàn)路及方向的中長(zhǎng)途客流在樞紐內(nèi)的同站、異站換乘優(yōu)化問(wèn)題進(jìn)行研究,以樞紐內(nèi)需求換乘客流的平均換乘滿(mǎn)意度最大化和樞紐車(chē)站列車(chē)到發(fā)均衡性最大化為目標(biāo),構(gòu)建了基于換乘導(dǎo)向的樞紐高鐵列車(chē)接續(xù)優(yōu)化模型,并設(shè)計(jì)改進(jìn)染色體編碼方式和選擇策略的遺傳算法進(jìn)行求解。

    1 問(wèn)題描述

    1.1 同站換乘過(guò)程

    通常,大型客運(yùn)樞紐高鐵站均設(shè)有專(zhuān)門(mén)的中轉(zhuǎn)換乘通道,同站換乘旅客可以實(shí)現(xiàn)不出站換乘,整個(gè)換乘過(guò)程為旅客下車(chē)步行至接續(xù)列車(chē)??空九_(tái)候車(chē)(或至專(zhuān)用候車(chē)室候車(chē))然后檢票完成上車(chē)。

    如圖1,定義從換出列車(chē)x到達(dá)車(chē)站h時(shí)刻ax起至接續(xù)列車(chē)x從同一車(chē)站h出發(fā)時(shí)刻dx止的時(shí)間為同站換乘時(shí)間,包括兩部分:一是旅客從換出列車(chē)下車(chē)步行至接續(xù)列車(chē)??空九_(tái)(或至專(zhuān)用候車(chē)室)時(shí)止的必要走行時(shí)間tb;二是旅客從到達(dá)接續(xù)列車(chē)站臺(tái)或?qū)S煤蜍?chē)室時(shí)起至接續(xù)列車(chē)出發(fā)時(shí)止的等待時(shí)間tw。

    tb按照實(shí)際查定取值,等待時(shí)間tw受兩個(gè)因素影響,一是接續(xù)列車(chē)開(kāi)行頻率和開(kāi)行均衡性,二是換乘銜接方案。理論上來(lái)講,等待時(shí)間tw最小值為0,最大值為接續(xù)列車(chē)的開(kāi)行間隔時(shí)間。

    1.2 異站換乘過(guò)程

    由于樞紐內(nèi)不同高鐵站銜接的干線(xiàn)方向不同,部分起點(diǎn)和終點(diǎn)不在一個(gè)方向和線(xiàn)路的中長(zhǎng)途旅客為了完成其旅行過(guò)程,通常需要在一個(gè)換乘樞紐內(nèi)的不同車(chē)站間進(jìn)行中轉(zhuǎn),而從一個(gè)車(chē)站到另一車(chē)站,中間需要乘坐城市公共交通。

    如圖2,定義從換出列車(chē)x到達(dá)車(chē)站h時(shí)刻ax起至接續(xù)列車(chē)x從另一車(chē)站h出發(fā)時(shí)刻dx止的時(shí)間為異站換乘時(shí)間,包括5個(gè)部分:一是旅客從換出列車(chē)下車(chē)步行至公共交通站臺(tái)的必要走行時(shí)間tb′;二是旅客從到達(dá)公共交通站臺(tái)到等待公共交通上車(chē)的公共交通候車(chē)時(shí)間tw′;三是旅客乘坐城市公共交通從h車(chē)站到另一接續(xù)車(chē)站h的公共交通旅行時(shí)間時(shí)間t公;四是旅客從另一公共交通站臺(tái)步行至接續(xù)車(chē)站檢票至候車(chē)室的必要走行時(shí)間tb″;五是旅客從接續(xù)車(chē)站候車(chē)室時(shí)起至完成上車(chē),接續(xù)列車(chē)出發(fā)時(shí)止的等待時(shí)間tw。

    其中tb′、tw′、t公按照不同城市公共交通的異站換乘客流分擔(dān)率、發(fā)車(chē)間隔、旅行速度等結(jié)合實(shí)際查定取值,參考文獻(xiàn)[12-15]可知,當(dāng)樞紐采用以軌道交通為主的集疏運(yùn)模式時(shí),軌道交通、公交、出租車(chē)和社會(huì)車(chē)輛的換乘客流分擔(dān)率一般分別為60%、30%、10%;tb″的取值要考慮旅客進(jìn)站排隊(duì)檢票和安檢過(guò)程;tw取值原則與同站換乘一致。

    2 樞紐內(nèi)換乘接續(xù)優(yōu)化模型的構(gòu)建

    2.1 模型假設(shè)

    本文的研究基于以下假設(shè):1)已知研究周期內(nèi)樞紐各個(gè)高鐵站的高速列車(chē)開(kāi)行方案,且不存在不停車(chē)通過(guò)樞紐的高速列車(chē);2)已知研究周期內(nèi)樞紐不同高速列車(chē)之間的需求換乘乘客流量預(yù)測(cè)值;3)不考慮高速列車(chē)運(yùn)行晚點(diǎn)情況;4)同一線(xiàn)路同一方向、同一線(xiàn)路相反方向間不存在高速列車(chē)接續(xù);5)研究周期內(nèi)的換乘客流在整個(gè)旅行過(guò)程中只進(jìn)行一次換乘,且在所研究的樞紐內(nèi)完成換乘;6)本文只考慮同方向列車(chē)的換乘接續(xù)情況。

    14Jmn(dmnj+2x″h-dmnjxh-2Emndh)2](4)

    其中:x,x′,x″∈Xmnah;x,x′,x″∈Xmndh;Emnah=(T2-T1)/Iahmn,Emnah表示樞紐h站連接的m線(xiàn)路n方向上的列車(chē)平均到達(dá)間隔時(shí)間;Emndh=(T2-T1)/Idhmn,Emndh表示樞紐h站連接的m線(xiàn)路n方向上的列車(chē)平均出發(fā)間隔時(shí)間。

    2.4 約束條件

    本文主要從高鐵站列車(chē)合理到發(fā)時(shí)間、換乘客流量、旅客換乘時(shí)間、車(chē)站間隔時(shí)間等方面對(duì)模型進(jìn)行約束:

    1)高速列車(chē)合理到發(fā)時(shí)間約束。

    amnixh-Tmnxh≥DmnE(5)

    dmnjxh-Tmnxh≤AmnL(6)

    T1≤amnixh≤T2(7)

    T1≤dmnjxh≤T2(8)

    式(5)表示到達(dá)樞紐車(chē)站的高速列車(chē)在其始發(fā)站的始發(fā)時(shí)刻不早于始發(fā)車(chē)站的最早發(fā)車(chē)時(shí)刻;式(6)表示離開(kāi)樞紐車(chē)站的高速列車(chē)到達(dá)其終點(diǎn)站的時(shí)刻不晚于終到車(chē)站的最晚終到時(shí)刻;式(7)、(8)表示樞紐車(chē)站的列車(chē)到發(fā)時(shí)刻應(yīng)在研究周期內(nèi)。

    2)換乘客流量約束。

    pqxjmnxi≥f(Qpqx*jmnxi,Q); x≠x*且m≠p(9)

    f(μ1, μ2)=1, μ1≥μ2

    0,μ1<μ2(10)

    表示m線(xiàn)路n方向上第i列到達(dá)樞紐且車(chē)次為x的列車(chē)與p線(xiàn)路q方向上第j列離開(kāi)樞紐且車(chē)次為x*的列車(chē)之間的接續(xù)關(guān)系由兩列車(chē)之間的換乘客流量決定。

    3)換乘時(shí)間約束。

    Tsmin≤pqx*jmnxi(dpqjx*h-amnixh)≤Tsmax;m≠p,x≠x*(11)

    Tdmin≤pqx*jmnxi(dpqjx*h-amnixh*)≤Tdmax;

    h≠h*, m≠p, x≠x*(12)

    式(11)、(12)表示若列車(chē)x與x*構(gòu)成接續(xù)關(guān)系,則兩列車(chē)在樞紐車(chē)站的出發(fā)和到達(dá)時(shí)刻差要滿(mǎn)足換乘時(shí)間約束。

    4)通過(guò)列車(chē)停站時(shí)間約束。

    thmin≤dpqjxh-amnixh≤thmax; m=p時(shí),q≠n(13)

    5)車(chē)站間隔時(shí)間約束。

    amni+1x′h-amnixh≥Iaa(14)

    dmnj+1x′*h-dmnjxh≥Idd(15)

    dmnjxh-amnixh≥Iad(16)

    amnixh-dmnjxh≥Ida(17)

    式(14)~(17)表示樞紐車(chē)站同方向兩列車(chē)到發(fā)的時(shí)間間隔應(yīng)滿(mǎn)足車(chē)站的列車(chē)到達(dá)、出發(fā)、到發(fā)、發(fā)到時(shí)間間隔約束。

    6)樞紐車(chē)站到發(fā)線(xiàn)能力約束。

    對(duì)于任意時(shí)刻t,有如下約束:

    G(ψv,t)=1, ψv≤t

    0,ψv>t; v∈{1,2}(18)

    ψ1(dmnjx*h,amnixh,whxx*)=

    whxx*(dmnjx*h-Tb),? whxx*=-1

    whxx*+1(amnixh-Tap),whxx*=0

    whxx*(amnixh-Tap), whxx*=1(19)

    ψ2(dmnjxh,amnixh,whxx)=

    whxx(dmnjxh-Tde),? whxx=-1

    whxx+1(dmnjxh-Tde),whxx=0

    whxx(amnixh-Ta),whxx=1(20)

    ∑2m=1∑2n=1∑Jmnj=1∑Imni=1[G(ψ1,t)-G(ψ2,t)]≤dfh(21)

    式(19)、(20)中引入?yún)?shù)whxx表示列車(chē)x或x在樞紐內(nèi)h站的停站方式,whxx取-1表示始發(fā)列車(chē),whxx取0表示通過(guò)列車(chē),whxx取1表示終到列車(chē)。ψ1(dmnjxh,amnixh,whxx)表示列車(chē)開(kāi)始占用到發(fā)線(xiàn)的時(shí)間函數(shù),ψ2(dmnjxh,amnixh,whxx)表示列車(chē)結(jié)束占用到發(fā)線(xiàn)的時(shí)間函數(shù)。

    7)區(qū)間通過(guò)能力約束。

    cra+crd≤CrT(22)

    3 算法設(shè)計(jì)及算例分析

    3.1 改進(jìn)遺傳算法設(shè)計(jì)

    遺傳算法在求解一些復(fù)雜的非線(xiàn)性問(wèn)題時(shí)能發(fā)揮其全局搜索優(yōu)勢(shì),且在鐵路運(yùn)輸組織優(yōu)化方面取得了廣泛應(yīng)用。由于模型為多目標(biāo)非線(xiàn)性整數(shù)規(guī)劃模型,為了便于求解,得到研究周期內(nèi)的樞紐車(chē)站所有列車(chē)到發(fā)時(shí)刻,因而將多目標(biāo)轉(zhuǎn)換成單目標(biāo)模型求解。由于兩個(gè)目標(biāo)函數(shù)的量綱和量級(jí)不同,首先應(yīng)統(tǒng)一量綱和量級(jí),使得目標(biāo)函數(shù)規(guī)范化,轉(zhuǎn)換后的目標(biāo)函數(shù)如式(23),式中Z2,min、Z2,max表示到發(fā)均衡性最小值和最大值,λ為加權(quán)系數(shù)。

    min Z=λ(-Z1)+(1-λ)Z2-Z2,minZ2,max-Z2,min(23)

    采用改進(jìn)染色體編碼和選擇策略的遺傳算法,設(shè)計(jì)的算法步驟如下:

    1)構(gòu)造基因編碼和染色體。

    本文關(guān)于換乘樞紐車(chē)站列車(chē)到發(fā)時(shí)間編排問(wèn)題的研究對(duì)象包括車(chē)次、線(xiàn)路、方向、到達(dá)(出發(fā))車(chē)站、到達(dá)(出發(fā))順序和時(shí)間。編碼時(shí)將車(chē)次和車(chē)次對(duì)應(yīng)的到達(dá)(出發(fā))時(shí)刻構(gòu)成一個(gè)元組,作為本文遺傳算法的基因,并采用兩段式染色體編碼方式,以換乘樞紐內(nèi)存在兩個(gè)高鐵站,樞紐連接兩條線(xiàn)路,且兩車(chē)站之間存在聯(lián)絡(luò)線(xiàn)為例,構(gòu)造如圖3所示的染色體編碼結(jié)構(gòu):第一部分表示樞紐內(nèi)兩個(gè)車(chē)站的各個(gè)線(xiàn)路各個(gè)方向的列車(chē)依次到達(dá)(出發(fā))時(shí)刻編排,依次為h站1線(xiàn)路1方向到達(dá)、h站1線(xiàn)路1方向出發(fā)、h站1線(xiàn)路2方向到達(dá)(出發(fā))、h站2線(xiàn)路1方向到達(dá)(出發(fā))、h站2線(xiàn)路2方向到達(dá)(出發(fā))以及h站各個(gè)線(xiàn)路和方向的列車(chē)到發(fā)時(shí)間編排;第二部分表示不同車(chē)站不同線(xiàn)路和方向的列車(chē)依次到達(dá)和出發(fā)時(shí)刻之間的斷點(diǎn)位置。

    2)初始種群的生成。

    由于研究周期[T1,T2]內(nèi)每個(gè)個(gè)體染色體基因中的實(shí)數(shù)編碼、整數(shù)編碼分別表示到達(dá)(離開(kāi))樞紐車(chē)站的列車(chē)車(chē)次和對(duì)應(yīng)列車(chē)在樞紐車(chē)站的到達(dá)(出發(fā))時(shí)刻,首先隨機(jī)生成一個(gè)樞紐車(chē)站所有線(xiàn)路和方向不重復(fù)的到達(dá)(出發(fā))車(chē)次序列,其次生成該車(chē)次序列對(duì)應(yīng)的列車(chē)到達(dá)(出發(fā))時(shí)刻,由于列車(chē)到達(dá)(出發(fā))時(shí)刻是按照時(shí)間序排列的,因此每個(gè)個(gè)體的染色體基因段中的列車(chē)到達(dá)(出發(fā))時(shí)刻應(yīng)該是Xmnah(Xmndh)這些車(chē)次的到達(dá)(出發(fā))時(shí)刻的升序排列。因此可以按照以下方法得到兩個(gè)斷點(diǎn)之間列車(chē)出發(fā)時(shí)刻升序排列的染色體片段:

    ①在[T1,T2]內(nèi)生成一個(gè)隨機(jī)數(shù)t1,作為dmn1xh;

    ②在[t1,T2]內(nèi)生成一個(gè)隨機(jī)數(shù)t2,作為dmn2x′h;

    ③依此類(lèi)推,在[tj-1,T2]內(nèi)隨機(jī)生成tj作為h站發(fā)往m線(xiàn)路n方向第j列車(chē)的出發(fā)時(shí)刻。

    同理可以生成染色體中任意兩個(gè)斷點(diǎn)之間的列車(chē)到達(dá)時(shí)刻升序排列片段,并將這些片段組成完整的染色體個(gè)體。

    3)適應(yīng)度函數(shù)。

    設(shè)K為種群規(guī)模,Zk是當(dāng)前種群中的第k個(gè)染色體的目標(biāo)函數(shù)值,Zmax和Zmin分別為當(dāng)前種群的最大、最小目標(biāo)函數(shù)值,gk為適應(yīng)度函數(shù)。

    gk=Zmax-Zk+γZmax-Zmin+γ; k∈{1,2,…,K}(24)

    其中:γ是(0,1)內(nèi)的隨機(jī)數(shù),引入γ使得當(dāng)染色體間的適應(yīng)值差距較大時(shí),采用適應(yīng)值比例選擇,當(dāng)染色體間的適應(yīng)值差距較小時(shí),則在相互競(jìng)爭(zhēng)的染色體中進(jìn)行純隨機(jī)的選擇。這樣使表現(xiàn)較優(yōu)的個(gè)體適應(yīng)度更大,進(jìn)入下一代的幾率也更大。

    4)改進(jìn)選擇算子。

    文中應(yīng)用如下混合改進(jìn)了的選擇方案,將選擇操作分為以下3個(gè)步驟:

    ①將計(jì)算得到的個(gè)體適應(yīng)度按照從大到小進(jìn)行排序,保留父代種群中的前1/4個(gè)體[16];

    ②隨機(jī)在父代種群中選取1/2K個(gè)個(gè)體,K為種群規(guī)模(排除①中選中的個(gè)體)按照下文的5)交叉算子和6)變異算子生成新子代[16];

    ③將①中選擇出的個(gè)體與②中生成的新子代進(jìn)行模擬退火選擇,其次將模擬退火選擇之后得到的個(gè)體與原父代種群合并組成臨時(shí)種群,設(shè)臨時(shí)種群規(guī)模為K′,臨時(shí)種群中每個(gè)個(gè)體為k′,則保證了隨著算法迭代次數(shù)的下降,改進(jìn)算法的臨時(shí)種群規(guī)模K′始終大于種群規(guī)模K。依據(jù)改進(jìn)方法,借鑒模擬退火算法,將各個(gè)染色體的適應(yīng)度函數(shù)值進(jìn)行拉伸[17],運(yùn)用式(25)計(jì)算臨時(shí)種群中每個(gè)個(gè)體k′被選擇的概率,從臨時(shí)種群中選出K個(gè)個(gè)體組成新的一代進(jìn)行下一輪遺傳算法操作:

    P(k′)=eg(k′)/Tv′/∑K′k′=0eg(k′)/Tv′; k′∈{1,2,…,K′}(25)

    Tv′=1ln(v′/T0+1)(26)

    其中:Tv′為漸趨于0的退火溫度(v′=1,2,…),T0為初始溫度。改進(jìn)的選擇算子借鑒了模擬退火算法中的拉伸方法,可以改善基本遺傳算法的早熟問(wèn)題,同時(shí)使得在遺傳算法后期過(guò)程中,當(dāng)個(gè)體染色體適應(yīng)度趨于一致時(shí),改進(jìn)的算法具有更好的選擇效果,同時(shí)這一選擇算子不僅包含了父代種群中的最佳個(gè)體,保證了算法的收斂性,又保留了最優(yōu)個(gè)體以外的優(yōu)秀個(gè)體,增加了種群的多樣性,避免算法陷入局部最優(yōu)解。

    5)交叉算子。

    文中的染色體是由不同片段的車(chē)次實(shí)數(shù)編碼和列車(chē)有序到(發(fā))時(shí)刻整數(shù)編碼元組構(gòu)成,采用傳統(tǒng)的交叉方法不能得到理想的結(jié)果,因此本文采用區(qū)域交叉,將交叉區(qū)域限定在兩個(gè)相鄰的斷點(diǎn)之間,并將交叉區(qū)域的到(發(fā))時(shí)刻值與另一個(gè)群體基因位的時(shí)刻值進(jìn)行比較,交叉時(shí)車(chē)次不隨時(shí)刻交叉,具體操作步驟如下:

    ①在父代種群中,隨機(jī)選擇兩個(gè)個(gè)體作為父代1和父代2,并隨機(jī)選取兩個(gè)相鄰斷點(diǎn)之間的兩個(gè)交叉點(diǎn),規(guī)定父代1和父代2的兩個(gè)相鄰斷點(diǎn)的選擇相同。將這兩個(gè)交叉點(diǎn)之間的所有時(shí)刻基因值作為交叉對(duì)象進(jìn)行交叉操作。

    ②將父代2交叉區(qū)域內(nèi)的時(shí)刻基因值分別與父代1中所有的時(shí)刻基因值之差的絕對(duì)值進(jìn)行比較。若差值為0,則父代1中的基因值不變;若差值不為0,則將交叉區(qū)域內(nèi)的時(shí)刻基因值與父代1中絕對(duì)值最小的基因位上的時(shí)刻值進(jìn)行替換,若存在兩個(gè)基因時(shí)刻值的絕對(duì)值相等,則將父代1中對(duì)應(yīng)的時(shí)刻基因值較大的時(shí)刻值替換,從而得到交叉后的子代1,同理可以得到交叉后的子代2。例如某兩個(gè)相鄰斷點(diǎn)之間的交叉操作如圖4和圖5。

    6)變異算子。

    在染色體個(gè)體變異之前,為了避免算法陷入局部最優(yōu),同時(shí)得到研究周期內(nèi)樞紐車(chē)站最優(yōu)的列車(chē)到發(fā)順序,應(yīng)隨機(jī)生成K個(gè)不同的樞紐車(chē)站所有線(xiàn)路方向不重復(fù)的到發(fā)車(chē)次序列替換上述交叉后的子代染色體基因中的車(chē)次序列。

    采用均勻變異,將每條染色體上的一個(gè)或多個(gè)基因點(diǎn)設(shè)置為變異點(diǎn),以一定的變異概率Pm確定一個(gè)隨機(jī)數(shù),并確定變異點(diǎn)的取值范圍,以變異后的值取代變異點(diǎn)先前的取值,變異時(shí)基因元組中的車(chē)次不變。例如,隨機(jī)選中的一個(gè)個(gè)體k,個(gè)體k中的amnixh為變異點(diǎn),amnixh的取值范圍為[amnixh-σ,amnixh+σ]。變異點(diǎn)的新基因值為:

    a′mnixh=amnixh-σ+β[(amnixh+σ)-(amnixh-σ)](27)

    其中β∈(0,1),是一個(gè)隨機(jī)數(shù)。

    7)終止條件。

    設(shè)定改進(jìn)的遺傳算法最大進(jìn)化代數(shù),當(dāng)遺傳算法運(yùn)行到最大進(jìn)化代數(shù)之后,停止運(yùn)行,輸出結(jié)果;否則轉(zhuǎn)步驟3)。

    3.2 算例分析

    選擇高速鐵路網(wǎng)中某一大型樞紐,樞紐銜接兩條高速鐵路線(xiàn)路,存在兩個(gè)高鐵站且兩車(chē)站之間存在聯(lián)絡(luò)線(xiàn),該樞紐采用以軌道交通為主的客流集疏運(yùn)模式,圖6是研究客運(yùn)樞紐的換乘接續(xù)示意圖。

    該樞紐詳細(xì)概況如下:dfh1、dfh2為14條、13條;thmin、thmax分別為3min、20min;Iaa、Idd、Ida、Iad分別為2min、2min、3min、3min;tap、tde、tb、ta分別為5min、2min、23min、20min;D11E、D12E、D21E、D22E分別為360、390、330、420;A11L、A12L、A21L、A22L分別為1380、1380、1440、1380;C1T、C2T、C3T、C4T分別為16、7、12、10。選

    取該樞紐的中長(zhǎng)途旅客換乘時(shí)段13:00—14:00作為研究周期,樞紐內(nèi)兩個(gè)高鐵站的列車(chē)開(kāi)行方案如表1所示。由于文獻(xiàn)[11-12,18]研究樞紐與文中所研究樞紐概況、換乘客流量相當(dāng),因而推算得出研究周期內(nèi)不同列車(chē)間的換乘客流量預(yù)測(cè)值,如表2所示。

    根據(jù)文獻(xiàn)[11,18],結(jié)合所研究樞紐的概況、換乘客流量、鐵路部門(mén)的成本和收益等因素綜合確定換乘接續(xù)列車(chē)判定接續(xù)旅客人數(shù)Q為50人。同時(shí)為了保證旅客能夠在一定時(shí)間內(nèi)順利完成換乘,Tsmin、Tsmax按照樞紐內(nèi)客流高峰時(shí)段的乘客必要走行時(shí)間、檢票時(shí)間以及接續(xù)列車(chē)開(kāi)行間隔時(shí)間來(lái)取值,分別取15min、45min;Tdmin、Tdmax結(jié)合樞紐內(nèi)客流高峰時(shí)段的乘客必要走行時(shí)間、檢票時(shí)間、接續(xù)列車(chē)開(kāi)行間隔時(shí)間、城市公共交通的運(yùn)行時(shí)間和換乘客流的分擔(dān)率取值,分別取40min、70min。

    遺傳算法的參數(shù)設(shè)置為:種群規(guī)模K取300,交叉概率pc取0.75,變異概率pm取0.06,最大迭代次數(shù)1000,T0初始溫度取為100,λ取0.6。依據(jù)上述模型和算法,利用Matlab R2014a軟件工具編程求解算例,分別運(yùn)用基本遺傳算法和改進(jìn)的遺傳算法對(duì)算例求解,得到改進(jìn)后遺傳算法的優(yōu)化結(jié)果如表3和表4,目標(biāo)函數(shù)為-0.1257,其中目標(biāo)值Z1滿(mǎn)意度為0.6367,目標(biāo)值Z2列車(chē)到發(fā)均衡性值為640.5826,以及目標(biāo)函數(shù)隨迭代次數(shù)變化示意圖8。

    鑒于本文設(shè)計(jì)改進(jìn)的遺傳算法結(jié)合了模擬退火算法的思想,因此,本文又以基本模擬退火算法作為對(duì)比算法,在同一臺(tái)計(jì)算機(jī)上利用Matlab編程實(shí)現(xiàn)基本模擬退火算法求解模型,通過(guò)進(jìn)行多次實(shí)驗(yàn),得到基本模擬退火算法計(jì)算的結(jié)果為:算法迭代了696次之后,目標(biāo)函數(shù)值趨于穩(wěn)定且到達(dá)最小,此時(shí)求得的目標(biāo)函數(shù)值為-0.1088,其中目標(biāo)值Z1滿(mǎn)意度為0.6186,目標(biāo)值Z2列車(chē)到發(fā)均衡性值為655.7341,進(jìn)一步將基本模擬退火算法、基本遺傳算法以及本文改進(jìn)的遺傳算法進(jìn)行比較,如表5所示。

    在樞紐客流換乘滿(mǎn)意度方面,分析表3和表4可以看出,12對(duì)列車(chē)實(shí)現(xiàn)了有效接續(xù),其中7對(duì)列車(chē)的接續(xù)時(shí)間使得換乘客流滿(mǎn)意度達(dá)到60%以上,1對(duì)列車(chē)的接續(xù)時(shí)間使得換乘客流滿(mǎn)意度最低。同時(shí)由表4求得該樞紐換乘客流的平均同站換乘時(shí)間為27min,平均異站換乘時(shí)間為49.83min。

    在樞紐列車(chē)到發(fā)均衡性方面,h2站連接的2線(xiàn)路2方向列車(chē)到達(dá)和2線(xiàn)路1方向列車(chē)出發(fā)、h1站連接的1線(xiàn)路2方向列車(chē)出發(fā)和2線(xiàn)路1方向列車(chē)到達(dá),這幾個(gè)方向的列車(chē)到發(fā)均

    衡性較優(yōu),而h1站連接的1線(xiàn)路列車(chē)到發(fā)均衡性較差??偟膩?lái)說(shuō),h2站的列車(chē)到發(fā)均衡性較h1站更優(yōu)。

    為了進(jìn)一步驗(yàn)證改進(jìn)后遺傳算法的穩(wěn)定性,分別運(yùn)用基本遺傳算法和改進(jìn)遺傳算法對(duì)算例求解,圖7是隨機(jī)抽取的50次實(shí)驗(yàn)結(jié)果的收斂代數(shù)情況??梢钥闯?,在種群規(guī)模和迭代次數(shù)相同的情況下,改進(jìn)后的遺傳算法穩(wěn)定性更好。

    利用本文模型和改進(jìn)的遺傳算法對(duì)高速鐵路網(wǎng)中某一大型客運(yùn)樞紐高速列車(chē)接續(xù)問(wèn)題進(jìn)行優(yōu)化,得到目標(biāo)函數(shù)值隨迭代次數(shù)變化情況,如圖8所示。

    由于本文設(shè)計(jì)的遺傳算法結(jié)合所研究問(wèn)題設(shè)計(jì)了兩段式染色體編碼方式,將車(chē)次和車(chē)次對(duì)應(yīng)的到達(dá)(出發(fā))時(shí)刻構(gòu)成一個(gè)元組即設(shè)計(jì)的算法基因,并在選擇算子中融入了模擬退火算法的拉伸方法和退火選擇,將遺傳算法全局尋優(yōu)能力強(qiáng)的特點(diǎn)和模擬退火算法局部尋優(yōu)能力強(qiáng)的特點(diǎn)相結(jié)合,使得改進(jìn)的遺傳算法具有更好的選擇效果,有效地提高了本文設(shè)計(jì)算法的解的質(zhì)量。由表5可知,對(duì)比基本模擬退火算法,本文改進(jìn)的遺傳算法,目標(biāo)函數(shù)值更優(yōu),其中旅客平均換乘滿(mǎn)意度增加約2.93%,到發(fā)均衡性增加約2.31%,旅客無(wú)效等待時(shí)間更少,h1站和h2站兩個(gè)車(chē)站的列車(chē)到發(fā)均衡性指標(biāo)值的差異更小;對(duì)比基本遺傳算法,本文改進(jìn)的遺傳算法使得旅客平均換乘滿(mǎn)意度增加約5.10%,到發(fā)均衡性增加約0.27%,且求得的目標(biāo)函數(shù)值、平均無(wú)效等待時(shí)間、h1站和h2站兩個(gè)車(chē)站的列車(chē)到發(fā)均衡性指標(biāo)差異均比基本遺傳算法求得的解更優(yōu)。

    由以上分析結(jié)果可知,本文設(shè)計(jì)改進(jìn)的遺傳算法更加適合所研究問(wèn)題的模型特征,具有更快的收斂速度和更優(yōu)的求解性能,因此,文中基于換乘的高速鐵路樞紐列車(chē)接續(xù)優(yōu)化模型和改進(jìn)的遺傳算法能夠?qū)Q乘樞紐的列車(chē)接續(xù)優(yōu)化問(wèn)題進(jìn)行有效求解。

    利用本文的模型和改進(jìn)遺傳算法對(duì)算例中的大型客運(yùn)樞紐高鐵列車(chē)接續(xù)進(jìn)行優(yōu)化,算法迭代了623次之后,目標(biāo)函數(shù)值趨于穩(wěn)定且達(dá)到最小,得到目標(biāo)函數(shù)值隨迭代次數(shù)變化情況如圖8所示。

    4 結(jié)語(yǔ)

    1) 通過(guò)優(yōu)化與編排周期內(nèi)換乘樞紐車(chē)站的所有列車(chē)到發(fā)順序的方法來(lái)研究換乘樞紐列車(chē)接續(xù)問(wèn)題。建立基于換乘導(dǎo)向的大型高速鐵路樞紐列車(chē)接續(xù)優(yōu)化模型,以換乘客流的平均換乘滿(mǎn)意度和樞紐車(chē)站列車(chē)到發(fā)均衡性為目標(biāo),考慮大站合理發(fā)車(chē)時(shí)間、換乘時(shí)間、車(chē)站間隔時(shí)間等約束,為解決換乘樞紐列車(chē)接續(xù)問(wèn)題提供了理論基礎(chǔ)和方法借鑒。

    2) 定義了換乘滿(mǎn)意度函數(shù),在求解模型時(shí)將多目標(biāo)問(wèn)題轉(zhuǎn)化成單目標(biāo)問(wèn)題,借助于改進(jìn)染色體編碼方式和選擇策略的遺傳算法,有效地求解了模型。研究結(jié)果優(yōu)化了換乘樞紐車(chē)站的列車(chē)接續(xù)問(wèn)題,驗(yàn)證了模型和算法的有效性。

    3) 通過(guò)對(duì)改進(jìn)的遺傳算法與基本遺傳算法、基本模擬退火算法進(jìn)行比較,可以看出本文設(shè)計(jì)的改進(jìn)遺傳算法不僅穩(wěn)定性和求解效率更高,且求得的解的質(zhì)量更優(yōu),對(duì)文中問(wèn)題求解具有較好的適應(yīng)性。

    4) 模型未考慮換出列車(chē)或接續(xù)列車(chē)發(fā)生晚點(diǎn)時(shí)的列車(chē)接續(xù)問(wèn)題,對(duì)于這一問(wèn)題的研究將是后續(xù)工作的重點(diǎn)。

    參考文獻(xiàn)

    [1]SHRIVASTAVA P, DHINGRA S L. Development of coordinated schedules using genetic algorithms [J]. Journal of Transportation Engineering, 2002, 128(1): 89-96.

    [2]YANG Z, KOSEKI T, SONE S. Railway network timetabling for reducing transfer time using genetic algorithm [J]. IEEJ Transactions on Sensors and Micromachines, 2008, 120(6): 802-809.

    [3]WONG R C W, YUEN T W Y, FUNG K W, et al. Optimizing timetable synchronization for rail mass transit [J]. Transportation Science, 2008, 42(1): 57-69.

    [4]MULLER T H J, FURTH P G. Transfer scheduling and control to reduce passenger waiting time [J]. Transportation Research Record: Journal of the Transportation Research Board, 2009, 3(2112): 111-118.

    [5]MESA J A, ORTEGA F A, POZO M A, et al. Rescheduling railway timetables in presence of passenger transfers between lines within a transportation network [M]// de SOUSA J F, ROSSI R. Computer-based Modelling and Optimization in Transportation, AISC 262. Berlin: Springer, 2014: 347-360.

    [6]HASSANNAYEBI E, ZEGORDI S H, YAGHINI M, et al. Timetable optimization models and methods for minimizing passenger waiting time at public transit terminals [J]. Transportation Planning and Technology, 2017, 40(3): 278-304.

    [7]禹丹丹,韓寶明,董寶田,等.基于換乘協(xié)同的軌道交通網(wǎng)列車(chē)時(shí)刻表優(yōu)化模型[J]. 中國(guó)鐵道科學(xué),2015,36(4):129-135.(YU D D, HAN B M, DONG B T, et al. Optimization model of train timetable for rail transit network based on transfer synchronization [J]. China Railway Science, 2015, 36(4): 129-135.)

    [8]姚鳳金,楊浩.旅客樞紐換乘時(shí)間調(diào)度研究[J].交通運(yùn)輸系統(tǒng)工程與信息,2008,8(1):133-137.(YAO F J, YANG H. On transferring time scheduling of passenger hubs [J]. Journal of Transportation Systems Engineering and Information Technology, 2008, 8(1): 133-137.)

    [9]郭瑾.高速鐵路綜合樞紐換乘條件評(píng)價(jià)及應(yīng)用研究[D].北京:北京交通大學(xué),2011:33-59.(GUO J. Research on evaluation and application of transfer condition for integrated hub of high-speed railway [D]. Beijing: Beijing Jiaotong University, 2011: 33-59.)

    [10]郭根材,聶磊,賀振歡,等.基于歐洲列車(chē)時(shí)刻表的高速列車(chē)銜接方案探討[J]. 北京交通大學(xué)學(xué)報(bào),2012,36(6): 21-26. (GUO G C, NIE L, HE Z H, et al. Train connection plan for high-speed railway based on European rail timetable [J]. Journal of Beijing Jiaotong University, 2012, 36(6): 21-26.)

    [11]黃曲.大型換乘站旅客列車(chē)接續(xù)方案優(yōu)化研究[D].成都:西南交通大學(xué),2018:38-39.(HUANG Q. Research on optimization of passenger train connection scheme in large transfer station [D]. Chengdu: Southwest Jiaotong University, 2018: 38-39.)

    [12]王建聰.城市客運(yùn)樞紐換乘組織關(guān)鍵問(wèn)題研究[D].北京:北京交通大學(xué),2006: 30-31.(WANG J C. Key issues on transfer organization of city passenger traffic hubs [D]. Beijing: Beijing Jiaotong University, 2006: 30-31.)

    [13]謝仲磊.高鐵樞紐與城市路網(wǎng)的銜接分析[J]. 交通與運(yùn)輸,2018,34(4):12-14.(XIE Z L. Analysis of connection between high-speed railway hub and urban road network [J]. Traffic and Transportation, 2018, 34(4): 12-14.)

    [14]張發(fā)才.鐵路客運(yùn)與城市交通運(yùn)營(yíng)組織銜接研究[D].南京:東南大學(xué),2006:24-28. (ZHANG F C. Study on passenger transfer of operation between railway station and urban traffic [D]. Nanjing: Southeast University, 2006: 24-28.)

    [15]周軍.高鐵客運(yùn)樞紐與城市交通銜接方式優(yōu)化研究[D].武漢:武漢理工大學(xué),2013: 44-46.(ZHOU J. Study on the optimization of high-speed passenger transport hubs urban transportation mode [D]. Wuhan: Wuhan University of Technology, 2013: 44-46.)

    [16]李政偉,譚國(guó)俊.改進(jìn)的退火遺傳優(yōu)化策略應(yīng)用研究[J].計(jì)算機(jī)工程與應(yīng)用,2010,46(4):245-248. (LI Z W, TAN G J. Research and application of improved simulated annealing genetic strategy [J]. Computer Engineering and Applications, 2010, 46(4): 245-248.)

    [17]劉月娥,何東健,李崢嶸.一種用于BP網(wǎng)絡(luò)優(yōu)化的并行模擬退火遺傳算法[J].計(jì)算機(jī)應(yīng)用,2006,26(1):204-206.(LIU Y E, HE D J, LI Z R. Parallel simulated anneling genetic algorithm for optimizing BP neural network [J]. Journal of Computer Applications, 2006, 26(1): 204-206.)

    [18]王瑩.面向超長(zhǎng)距離客流的高速鐵路網(wǎng)絡(luò)化列車(chē)開(kāi)行方案研究[D].北京:北京交通大學(xué),2018:138-140.(WANG Y. Study on high-speed railway network line planning for extra-long distance passenger flow [D]. Beijing: Beijing Jiaotong University, 2018: 138-140.)

    This work is partially supported by the National Natural Science Foundation of China (71861022, 61563028).

    QIAO Jun, born in 1993, M. S. candidate. Her research interests include operation management and decision optimization of rail transit.

    MENG Xuelei, born in 1979, Ph. D., professor. His research interests include operation management and decision optimization of rail transit.

    WANG Dongxian, born in 1992, M. S. candidate. His research interests include operation management and decision optimization of rail transit.

    TANG Lin, born in 1989, M. S. His research interests include operation management and decision optimization of rail transit.

    色哟哟·www| 国产精品爽爽va在线观看网站| 1000部很黄的大片| 免费观看在线日韩| 欧美一级a爱片免费观看看| 免费在线观看成人毛片| 男人舔女人下体高潮全视频| 亚洲久久久久久中文字幕| 综合色丁香网| 少妇的逼好多水| 国产老妇女一区| 国产高清三级在线| 免费看av在线观看网站| 九九热线精品视视频播放| 久久99精品国语久久久| 久久欧美精品欧美久久欧美| 九九在线视频观看精品| 91精品一卡2卡3卡4卡| 国产精品电影一区二区三区| 成年免费大片在线观看| 18禁动态无遮挡网站| 一边摸一边抽搐一进一小说| 99在线人妻在线中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 亚洲成人中文字幕在线播放| 午夜日本视频在线| 国产亚洲精品av在线| 色视频www国产| 国产亚洲av嫩草精品影院| 最近中文字幕高清免费大全6| 天堂√8在线中文| 国产男人的电影天堂91| 国产极品精品免费视频能看的| 国产综合懂色| 午夜激情福利司机影院| 久久欧美精品欧美久久欧美| 亚洲成人精品中文字幕电影| 男人舔女人下体高潮全视频| 午夜a级毛片| 少妇被粗大猛烈的视频| 少妇被粗大猛烈的视频| 国产一级毛片七仙女欲春2| 伦精品一区二区三区| 中文在线观看免费www的网站| 国产淫片久久久久久久久| 亚洲美女视频黄频| 汤姆久久久久久久影院中文字幕 | 美女内射精品一级片tv| 丰满乱子伦码专区| 高清av免费在线| 男女下面进入的视频免费午夜| 国产一区二区三区av在线| av免费观看日本| 亚洲国产精品成人久久小说| 成人性生交大片免费视频hd| 亚洲精品乱码久久久v下载方式| 亚洲丝袜综合中文字幕| 亚洲成色77777| 久久99精品国语久久久| 色吧在线观看| 美女xxoo啪啪120秒动态图| 2021少妇久久久久久久久久久| 亚洲久久久久久中文字幕| 免费看光身美女| 日本免费一区二区三区高清不卡| 国产亚洲最大av| 最近的中文字幕免费完整| 一级毛片aaaaaa免费看小| 成人毛片a级毛片在线播放| 最近最新中文字幕免费大全7| 成年女人永久免费观看视频| 国产精品综合久久久久久久免费| 最近2019中文字幕mv第一页| 一夜夜www| 综合色丁香网| av卡一久久| 美女内射精品一级片tv| 亚洲欧美日韩高清专用| 欧美成人a在线观看| 国产精品久久久久久精品电影小说 | 欧美成人免费av一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲欧洲国产日韩| 美女内射精品一级片tv| 免费观看精品视频网站| 午夜福利网站1000一区二区三区| 97在线视频观看| 午夜激情福利司机影院| 黄色配什么色好看| 亚洲精品乱久久久久久| 成人性生交大片免费视频hd| 特级一级黄色大片| www.色视频.com| 午夜老司机福利剧场| 丝袜喷水一区| 日韩精品青青久久久久久| 乱人视频在线观看| 国产91av在线免费观看| 日本黄色视频三级网站网址| 天堂√8在线中文| 最近2019中文字幕mv第一页| 亚洲性久久影院| av国产免费在线观看| 一二三四中文在线观看免费高清| 午夜久久久久精精品| 日本免费在线观看一区| 久久精品综合一区二区三区| 91精品国产九色| 欧美不卡视频在线免费观看| 亚洲国产欧美在线一区| 久99久视频精品免费| 精品久久久久久久人妻蜜臀av| 男人舔奶头视频| 深爱激情五月婷婷| 国产精品永久免费网站| 国产免费一级a男人的天堂| 男女边吃奶边做爰视频| 99久国产av精品| 伦理电影大哥的女人| 亚洲国产精品专区欧美| 国产亚洲精品久久久com| 欧美97在线视频| 久久精品人妻少妇| 久久综合国产亚洲精品| 亚洲不卡免费看| 夜夜看夜夜爽夜夜摸| 欧美精品国产亚洲| 欧美3d第一页| 国产av码专区亚洲av| 亚洲av免费在线观看| 国产人妻一区二区三区在| 日韩欧美在线乱码| 男女下面进入的视频免费午夜| 能在线免费看毛片的网站| 国产成人a区在线观看| ponron亚洲| 最近视频中文字幕2019在线8| 免费看av在线观看网站| 亚洲国产高清在线一区二区三| 国产黄a三级三级三级人| 少妇熟女欧美另类| 欧美最新免费一区二区三区| av又黄又爽大尺度在线免费看 | 国产综合懂色| 永久免费av网站大全| 久久久久久国产a免费观看| 2021少妇久久久久久久久久久| 亚洲欧美精品综合久久99| 免费av不卡在线播放| 亚洲精品乱码久久久久久按摩| 久久人人爽人人爽人人片va| 久久草成人影院| 一级黄色大片毛片| a级毛片免费高清观看在线播放| a级毛色黄片| 成人高潮视频无遮挡免费网站| 精品国产三级普通话版| 国产精品久久久久久精品电影小说 | 久久国产乱子免费精品| 欧美xxxx黑人xx丫x性爽| 天堂网av新在线| 亚洲精品日韩在线中文字幕| 久久午夜福利片| 国产色婷婷99| 日韩中字成人| 亚洲欧美日韩高清专用| 欧美激情久久久久久爽电影| 国产国拍精品亚洲av在线观看| 午夜激情福利司机影院| 美女脱内裤让男人舔精品视频| 爱豆传媒免费全集在线观看| 婷婷色av中文字幕| 国产成人精品久久久久久| 久久久久久大精品| 高清在线视频一区二区三区 | videossex国产| 亚洲三级黄色毛片| 久久精品久久久久久噜噜老黄 | 国产精品国产三级专区第一集| av在线播放精品| 国产精品麻豆人妻色哟哟久久 | 欧美性猛交黑人性爽| 免费观看人在逋| 精品久久久久久成人av| 日韩一本色道免费dvd| 亚洲av熟女| 男女那种视频在线观看| 久久国产乱子免费精品| 毛片女人毛片| 日韩在线高清观看一区二区三区| 最新中文字幕久久久久| 亚洲欧美一区二区三区国产| 亚洲色图av天堂| 国产精品一及| 久久精品久久久久久久性| 亚洲在久久综合| 精品久久久久久电影网 | 久久精品91蜜桃| 亚洲内射少妇av| 成人性生交大片免费视频hd| 天堂影院成人在线观看| 18+在线观看网站| 久久久久九九精品影院| 精品久久久久久久末码| 国产真实乱freesex| 51国产日韩欧美| 国产精品久久电影中文字幕| 日韩欧美国产在线观看| 国产精品野战在线观看| 赤兔流量卡办理| 亚洲国产高清在线一区二区三| 美女大奶头视频| 亚洲国产精品专区欧美| 精品久久久久久久久av| 欧美成人a在线观看| 成人av在线播放网站| 韩国av在线不卡| 精品久久久久久成人av| 国产淫语在线视频| 五月伊人婷婷丁香| 欧美色视频一区免费| 久久久久免费精品人妻一区二区| 国产精品一及| 啦啦啦啦在线视频资源| 亚洲欧美精品自产自拍| 日本av手机在线免费观看| 高清在线视频一区二区三区 | 久久综合国产亚洲精品| 啦啦啦韩国在线观看视频| 欧美日本亚洲视频在线播放| 国产精品熟女久久久久浪| 日韩欧美在线乱码| 国产成人免费观看mmmm| 啦啦啦啦在线视频资源| 三级国产精品欧美在线观看| 久久6这里有精品| 大又大粗又爽又黄少妇毛片口| 久久99蜜桃精品久久| 毛片一级片免费看久久久久| 亚洲精华国产精华液的使用体验| 日日撸夜夜添| 国产69精品久久久久777片| 九九在线视频观看精品| av在线播放精品| 99热精品在线国产| 国产成人免费观看mmmm| 亚洲人成网站在线播| 中文天堂在线官网| 亚洲国产最新在线播放| 国产毛片a区久久久久| 国产成人精品婷婷| av线在线观看网站| 亚洲乱码一区二区免费版| 国产亚洲精品久久久com| 久久鲁丝午夜福利片| 国产黄a三级三级三级人| 欧美成人午夜免费资源| 国产一区二区三区av在线| 久久久久网色| 一级黄色大片毛片| 日日干狠狠操夜夜爽| 国产精品野战在线观看| 亚洲国产欧美在线一区| 亚洲天堂国产精品一区在线| 黄片wwwwww| 成人av在线播放网站| 亚洲成色77777| 听说在线观看完整版免费高清| 久久人妻av系列| 久久国产乱子免费精品| 日本黄色片子视频| 日韩av在线免费看完整版不卡| 一级毛片我不卡| 国产成人freesex在线| 精品久久久久久久久av| 国产免费又黄又爽又色| 最近最新中文字幕免费大全7| 91精品国产九色| 久久精品人妻少妇| 国产人妻一区二区三区在| 欧美区成人在线视频| 日韩成人av中文字幕在线观看| 久久精品熟女亚洲av麻豆精品 | 亚洲内射少妇av| 亚洲国产精品国产精品| 亚洲欧洲国产日韩| 韩国高清视频一区二区三区| 亚洲精品日韩在线中文字幕| 精品欧美国产一区二区三| 国产精品一区二区在线观看99 | 欧美成人a在线观看| 久久久久久久久中文| 国产成人免费观看mmmm| 亚洲人成网站在线播| 汤姆久久久久久久影院中文字幕 | 国产黄色视频一区二区在线观看 | 中文字幕精品亚洲无线码一区| 高清午夜精品一区二区三区| 日韩av在线大香蕉| 在线播放无遮挡| 高清av免费在线| 在线播放国产精品三级| 久久国产乱子免费精品| 免费看av在线观看网站| 日产精品乱码卡一卡2卡三| 国产精品嫩草影院av在线观看| 三级经典国产精品| 成人午夜高清在线视频| 久久久久免费精品人妻一区二区| 久久久精品欧美日韩精品| 99热精品在线国产| 长腿黑丝高跟| 欧美一级a爱片免费观看看| 亚洲伊人久久精品综合 | 日韩一本色道免费dvd| 毛片女人毛片| 午夜福利在线观看免费完整高清在| 天天躁日日操中文字幕| 国产白丝娇喘喷水9色精品| 一级二级三级毛片免费看| 可以在线观看毛片的网站| 国产精品,欧美在线| 免费黄网站久久成人精品| 2021天堂中文幕一二区在线观| 麻豆成人午夜福利视频| 欧美三级亚洲精品| 美女cb高潮喷水在线观看| 少妇高潮的动态图| 久久精品久久精品一区二区三区| 搞女人的毛片| 校园人妻丝袜中文字幕| 日本熟妇午夜| 精品一区二区三区人妻视频| 亚洲va在线va天堂va国产| 国语自产精品视频在线第100页| 91在线精品国自产拍蜜月| 亚洲欧美成人精品一区二区| 久久久精品欧美日韩精品| 国产精品野战在线观看| 长腿黑丝高跟| 秋霞伦理黄片| 麻豆国产97在线/欧美| 黄色欧美视频在线观看| 久久99热这里只有精品18| 如何舔出高潮| av免费观看日本| 99热全是精品| 亚洲精品日韩在线中文字幕| 日韩一区二区三区影片| 亚洲欧美精品专区久久| 美女高潮的动态| 日本五十路高清| 51国产日韩欧美| 男的添女的下面高潮视频| 男人和女人高潮做爰伦理| 精品不卡国产一区二区三区| 一级毛片我不卡| 一级黄色大片毛片| 狠狠狠狠99中文字幕| 国产精品麻豆人妻色哟哟久久 | 国产免费福利视频在线观看| 好男人视频免费观看在线| 哪个播放器可以免费观看大片| 一夜夜www| 亚洲av不卡在线观看| 欧美区成人在线视频| 亚洲国产精品专区欧美| 亚洲国产色片| 国产精品99久久久久久久久| 亚洲欧美精品专区久久| 精品免费久久久久久久清纯| 中文字幕人妻熟人妻熟丝袜美| 男女啪啪激烈高潮av片| 麻豆国产97在线/欧美| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利高清视频| 三级毛片av免费| 久久热精品热| 99热这里只有是精品在线观看| 2022亚洲国产成人精品| 国产视频内射| 91av网一区二区| 免费看日本二区| 日本午夜av视频| 国产v大片淫在线免费观看| 亚洲欧美日韩无卡精品| 99久国产av精品| 精品久久久久久久久亚洲| 午夜精品国产一区二区电影 | 久久精品国产亚洲网站| 亚洲成人av在线免费| 26uuu在线亚洲综合色| 日韩欧美精品v在线| 精品一区二区三区视频在线| 成年版毛片免费区| 九色成人免费人妻av| 国产亚洲一区二区精品| 在线观看免费高清a一片| 少妇熟女欧美另类| 国产免费福利视频在线观看| 性色avwww在线观看| 日韩成人伦理影院| 亚洲精品色激情综合| 精品一区在线观看国产| 99久久精品国产国产毛片| 香蕉精品网在线| 国产女主播在线喷水免费视频网站| 日本av手机在线免费观看| 日本欧美视频一区| 国产一区二区在线观看av| 妹子高潮喷水视频| 一级片免费观看大全| 国产在视频线精品| 国产成人精品福利久久| 男女午夜视频在线观看 | 久久精品久久久久久噜噜老黄| 一级黄片播放器| 精品福利永久在线观看| 免费大片18禁| 又黄又爽又刺激的免费视频.| av网站免费在线观看视频| 制服人妻中文乱码| 亚洲欧美色中文字幕在线| 中文天堂在线官网| 十分钟在线观看高清视频www| 午夜激情av网站| 免费看光身美女| 欧美国产精品一级二级三级| 国产乱人偷精品视频| 99国产综合亚洲精品| 高清av免费在线| 18禁观看日本| 人人妻人人澡人人爽人人夜夜| 亚洲伊人久久精品综合| 欧美成人午夜精品| 中文字幕另类日韩欧美亚洲嫩草| 久久人人爽av亚洲精品天堂| 亚洲少妇的诱惑av| 1024视频免费在线观看| 国产精品久久久久成人av| 免费大片黄手机在线观看| 男女边吃奶边做爰视频| 又黄又爽又刺激的免费视频.| 2018国产大陆天天弄谢| 午夜福利,免费看| 高清不卡的av网站| 免费观看av网站的网址| 亚洲精品第二区| 亚洲熟女精品中文字幕| 亚洲国产精品国产精品| 桃花免费在线播放| 久久99热6这里只有精品| 成人黄色视频免费在线看| 久久久久久人妻| 成人国语在线视频| 午夜视频国产福利| 久久青草综合色| 亚洲经典国产精华液单| 欧美国产精品va在线观看不卡| 男女啪啪激烈高潮av片| 香蕉精品网在线| 夫妻性生交免费视频一级片| 中文乱码字字幕精品一区二区三区| 婷婷色麻豆天堂久久| 国产成人欧美| 我的女老师完整版在线观看| av视频免费观看在线观看| 日本av手机在线免费观看| 亚洲人成77777在线视频| 女性被躁到高潮视频| 丝袜人妻中文字幕| 中文字幕制服av| 十分钟在线观看高清视频www| 成年女人在线观看亚洲视频| 在线观看美女被高潮喷水网站| 亚洲成色77777| 在现免费观看毛片| 欧美精品一区二区大全| 啦啦啦啦在线视频资源| 七月丁香在线播放| 亚洲人成网站在线观看播放| 啦啦啦视频在线资源免费观看| av.在线天堂| 免费黄色在线免费观看| 又大又黄又爽视频免费| 日韩电影二区| 人人澡人人妻人| 成人亚洲精品一区在线观看| 美女国产高潮福利片在线看| 香蕉国产在线看| 久久精品国产鲁丝片午夜精品| 亚洲av成人精品一二三区| 晚上一个人看的免费电影| 日本-黄色视频高清免费观看| 人成视频在线观看免费观看| 9191精品国产免费久久| 成年人免费黄色播放视频| 成人影院久久| 国产色婷婷99| 18禁裸乳无遮挡动漫免费视频| 亚洲伊人色综图| 免费在线观看黄色视频的| 国产精品蜜桃在线观看| 在线观看免费日韩欧美大片| 久久99热6这里只有精品| 国产精品国产三级专区第一集| 熟女人妻精品中文字幕| 人成视频在线观看免费观看| 男女边吃奶边做爰视频| 侵犯人妻中文字幕一二三四区| 久久久久久久久久久久大奶| 最近的中文字幕免费完整| 日韩一本色道免费dvd| 日韩一区二区三区影片| 亚洲色图 男人天堂 中文字幕 | 欧美变态另类bdsm刘玥| 亚洲精品国产色婷婷电影| 国产不卡av网站在线观看| 人人澡人人妻人| 亚洲中文av在线| 成人无遮挡网站| 插逼视频在线观看| 夜夜骑夜夜射夜夜干| 免费女性裸体啪啪无遮挡网站| 国产亚洲精品久久久com| 日本-黄色视频高清免费观看| 又黄又粗又硬又大视频| 国产福利在线免费观看视频| 亚洲欧美日韩另类电影网站| 汤姆久久久久久久影院中文字幕| 久久久久国产精品人妻一区二区| 视频在线观看一区二区三区| 欧美日韩av久久| 亚洲精品乱码久久久久久按摩| 亚洲国产精品成人久久小说| 午夜福利网站1000一区二区三区| a 毛片基地| 中文欧美无线码| 国产一级毛片在线| 午夜av观看不卡| 国产精品久久久久久久久免| 欧美精品一区二区免费开放| 国产亚洲一区二区精品| 久久99蜜桃精品久久| 蜜桃国产av成人99| 成人黄色视频免费在线看| 如何舔出高潮| 又大又黄又爽视频免费| 国产深夜福利视频在线观看| 日韩制服骚丝袜av| 精品亚洲成a人片在线观看| 日韩三级伦理在线观看| 久久国产精品男人的天堂亚洲 | 国产色爽女视频免费观看| 欧美精品国产亚洲| 人妻 亚洲 视频| 国产视频首页在线观看| 欧美激情国产日韩精品一区| 久久精品久久精品一区二区三区| 女性生殖器流出的白浆| 精品国产国语对白av| 狠狠婷婷综合久久久久久88av| 老司机亚洲免费影院| 纵有疾风起免费观看全集完整版| 美女国产高潮福利片在线看| 国产黄色免费在线视频| 黄色 视频免费看| 极品人妻少妇av视频| 丁香六月天网| 中文字幕最新亚洲高清| 插逼视频在线观看| 国产一级毛片在线| 2018国产大陆天天弄谢| 性高湖久久久久久久久免费观看| 最近2019中文字幕mv第一页| av片东京热男人的天堂| 深夜精品福利| 曰老女人黄片| 精品国产一区二区久久| 欧美变态另类bdsm刘玥| 午夜福利影视在线免费观看| 香蕉国产在线看| 国产淫语在线视频| 99久久综合免费| 亚洲精品日韩在线中文字幕| 久久久a久久爽久久v久久| 久久99热这里只频精品6学生| 边亲边吃奶的免费视频| 亚洲伊人色综图| 两个人看的免费小视频| 日本免费在线观看一区| 欧美成人精品欧美一级黄| √禁漫天堂资源中文www| 90打野战视频偷拍视频| 丝瓜视频免费看黄片| 欧美精品av麻豆av| 插逼视频在线观看| 飞空精品影院首页| 久久久精品94久久精品| 国产高清不卡午夜福利| 性色av一级| 九色成人免费人妻av| 亚洲国产日韩一区二区| 九色亚洲精品在线播放| 女人精品久久久久毛片| 不卡视频在线观看欧美| 熟女人妻精品中文字幕| 色吧在线观看| 捣出白浆h1v1| 久久久久国产精品人妻一区二区| 亚洲一级一片aⅴ在线观看| 久久久久精品人妻al黑| 水蜜桃什么品种好| 日韩大片免费观看网站| 国产在线视频一区二区|