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

    考慮換道壓力增益的高速公路分流區(qū)跟馳模型

    2025-08-15 00:00:00楊龍海車婷婷馬文奎章錫俏陳霖
    重慶大學(xué)學(xué)報(bào) 2025年8期
    關(guān)鍵詞:交通流匝道標(biāo)準(zhǔn)差

    中圖分類號(hào):U491.265 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):1000-582X(2025)08-040-14

    doi:10.11835/j.issn.1000-582X.2024.005

    Research on car-following model for freeway diversion area considering lane-changing pressure gain

    YANG Longhai,CHE Tingting,MA Wenkui, ZHANG Xiqiao,CHEN Lin (School ofTransportation Science and Engineering,Harbin Institute ofTechnology,Harbin150o90,P.R.China)

    Abstract: To more accurately capture vehicle-following behavior in freeway diversion areas,this study proposes anenhanced car-following model by incorporating a lane-changing pressure gain factor into the full velocity diference model (FVDM).The proposed model accounts for both lane-changing pressure and lane-changing behavior. First,a linear stability analysis reveals that the stability region of the freeway diversion area diminishes as the lane-changing pressure gain factor increases.Second, using trajectory data from 92 lane-changing vehicles extracted from the NGSIM dataset,the improved car-folowing model is calibrated andvalidated.Simulation results demonstrate that the improved model more accurately reproduces vehicle speed and position.Compared with the original FVDM,the simulation error of the proposed pressure-based FVDM (P-FVDM) is reduced by 16% .Compared with the lane pressre FVDM(LP-FVDM),the proposed pressure-lane pressure FVDM(PLPFVDM) reduces error by 12% .Finally, the improved model is used to simulate trafc oscillation in the diversion area.Results show that lane changing behavior can trigger traffic oscilations.Higher driving speeds attenuate oscillation amplitudes,while lower speeds result in a concave growth patrn of speed standard deviation along the vehicle platoon in the upstream direction.Furthermore,lane changes occurring closer to the exit ramp exacerbate oscillations, whereas shorter lane-changing duration help suppress them.

    Keywords: transportation engineering; car-following model; traffic oscilation; lane-changing pressure; freeway diversionarea;numerical simulation

    “交通震蕩\"是一種車輛“走走停停\"的現(xiàn)象[2],對(duì)交通流的運(yùn)行有嚴(yán)重的負(fù)面影響,如:燃料消耗加劇、有害氣體排放量大、安全風(fēng)險(xiǎn)劇增、增加行車延誤等。換道行為是高速公路中非常普遍的現(xiàn)象,是引發(fā)交通震蕩的一個(gè)重要誘因[3-4],具有非常高的研究?jī)r(jià)值,由于其過(guò)程的復(fù)雜性,目前的定量研究相對(duì)較少。

    換道行為研究主要是從駕駛員行為和道路交通狀況的角度研究城市和高速公路的變道特征和影響因素。Moridpour等[5描述了車輛換道的過(guò)程,主要包括換道動(dòng)機(jī)的產(chǎn)生、換道時(shí)機(jī)的選擇及換道的執(zhí)行過(guò)程。根據(jù)研究動(dòng)機(jī),可將變道行為分為強(qiáng)制性變道和任意性變道。張穎達(dá)等利用Next Generation Simulation(NGSIM)數(shù)據(jù)集中的軌跡數(shù)據(jù),提出半對(duì)數(shù)和多項(xiàng)式模型來(lái)研究強(qiáng)制型和自由型換道行為的換道時(shí)間差異。換道行為中換道位置距離交叉口的距離、行駛速度和換道持續(xù)時(shí)間對(duì)交通流的穩(wěn)定性有一定的影響8-]。Toledo等[12分析了換道時(shí)間的分布范圍,發(fā)現(xiàn)了對(duì)換道時(shí)間產(chǎn)生影響的最主要因素是具有換道行為車輛與當(dāng)前車道前后相鄰車輛的速度差值和與目標(biāo)車道前后相鄰車輛的間距值。

    影響交通震蕩的因素包括縱向相互作用的不穩(wěn)定性、跟馳和換道行為、加減速不對(duì)稱特性、駕駛員類型以及道路線形等[13-14]。高速公路分流區(qū)域頻發(fā)的換道行為比交通流縱向不穩(wěn)定性更易誘發(fā)交通震蕩[15]。Mauch等通過(guò)研究高速公路交通流軌跡數(shù)據(jù)發(fā)現(xiàn)換道行為對(duì)交通震蕩的生成和發(fā)展有著非常大的影響。張娟采用累積移動(dòng)序列法研究上海市快速路匝道交通震蕩的特性發(fā)現(xiàn),震蕩傳播超越駛?cè)朐训篮?,振幅被削弱,穿過(guò)駛出匝道后被放大。黃永現(xiàn)[8為探究交通震蕩的發(fā)展特點(diǎn),進(jìn)行了3次實(shí)車實(shí)驗(yàn),大致計(jì)算出車隊(duì)臨界速度的區(qū)間為 30~40km/h ,并將其用來(lái)度量交通流的穩(wěn)定性。Yao等利用現(xiàn)場(chǎng)試驗(yàn)發(fā)現(xiàn),2輛連續(xù)車輛之間的時(shí)間間隔與前車的速度標(biāo)準(zhǔn)差和后車的初始速度呈負(fù)相關(guān),且與前車的平均速度和初始間距呈正相關(guān)。這些發(fā)現(xiàn)為更深入的理解交通震蕩現(xiàn)象提供了基礎(chǔ)。

    1 數(shù)據(jù)來(lái)源與篩選

    1.1 數(shù)據(jù)來(lái)源

    NGSIM數(shù)據(jù)集是對(duì)攝像機(jī)獲得的視頻數(shù)據(jù)應(yīng)用技術(shù)手段提取而得到的,收集了詳細(xì)、高質(zhì)量的交通流數(shù)據(jù),精度較高,多次被研究者們應(yīng)用于微觀駕駛行為分析、模型參數(shù)標(biāo)定等研究。NGSIM數(shù)據(jù)集包括Peachtree、Lankershim Boulevard、I-80和US-101共4條路段的數(shù)據(jù),其中,前面2條是城市主干道,后面2條是高速公路。

    為研究高速公路分流區(qū)的交通震蕩現(xiàn)象,選取的數(shù)據(jù)采集區(qū)域位于加利福尼亞州的US-101高速公路路段,道路渠化信息如圖1所示。該區(qū)域道路長(zhǎng)約2100ft(約 640.08m ),共有8條車道(主路車道5條,減速車道1條,匝道2條),每條車道的寬度為 3.66m ,最高限速 112.65km/h 。數(shù)據(jù)采集于2005年6月15日早晨07:50~08:05、08:05~08:20和08:20~08:35共3個(gè)時(shí)間段,采集對(duì)象為固定時(shí)間段和固定路段上的全部車輛,數(shù)據(jù)采集間隔0.1s,即1s產(chǎn)生10條數(shù)據(jù)。

    研究選取分流區(qū)內(nèi)換道車輛的駛出路線為從主路最右側(cè)車道5換道到減速車道6最后再到出口匝道車道8,換道駛出車輛必須在車道5進(jìn)行強(qiáng)制換道,主要以車道5、車道6和車道8的車輛數(shù)據(jù)為基礎(chǔ)進(jìn)行研究。

    Fig.1 Schematic diagram ofUS-101collection section

    1.2 數(shù)據(jù)處理

    1.2.1 數(shù)據(jù)篩選

    選擇的數(shù)據(jù)處理工具為python編程語(yǔ)言,對(duì)數(shù)據(jù)進(jìn)行初步預(yù)處理,具體步驟及原則如下。

    1)以從車道5換道駛?cè)胲嚨?時(shí)車輛跨越車道線的時(shí)刻為依據(jù),提取此時(shí)刻前后15s(150條)的數(shù)據(jù)作為換道車輛行為分析的原始數(shù)據(jù)。

    2)提取具有換道行為的車輛軌跡數(shù)據(jù)。車輛的換道行為必須從車道5開(kāi)始,換道至車道6,最后從出口匝道車道8駛出。剔除既沒(méi)有換道行為,也沒(méi)有對(duì)車輛的換道行為產(chǎn)生影響的車輛。

    3)剔除發(fā)生連續(xù)多次換道和換道失敗等行為的數(shù)據(jù),提取從車道5至車道6時(shí)只發(fā)生一次換道行為的車輛數(shù)據(jù)。

    4)換道開(kāi)始時(shí)須保證車道6的上游匝道無(wú)車輛,否則將會(huì)對(duì)換道車輛的換道行為產(chǎn)生影響,進(jìn)而會(huì)對(duì)研究分流區(qū)交通震蕩現(xiàn)象產(chǎn)生影響。

    經(jīng)上述篩選后,共提取出符合研究需求的換道車輛92輛。為研究換道駛出車輛的換道位置距出口匝道鼻端的距離和換道時(shí)間長(zhǎng)度對(duì)分流區(qū)交通震蕩的影響,需提取換道車輛的換道時(shí)間和換道位置進(jìn)行分析,如圖2所示。由圖2可知,大部分車輛的換道行為都在3~12 s完成,只有極少數(shù)車輛的換道時(shí)間大于12s,同時(shí)可以發(fā)現(xiàn)當(dāng)換道起始點(diǎn)距出口匝道鼻端較近時(shí),駕駛員完成換道行為的時(shí)間較短;當(dāng)距鼻端較遠(yuǎn)時(shí),換道時(shí)間較長(zhǎng)。大部分駕駛員的換道時(shí)間大于4s,只有極少數(shù)的駕駛員能夠在較短時(shí)間內(nèi)完成換道。因此,可以得出,大多數(shù)駕駛員會(huì)盡可能地選擇距出口匝道鼻端較遠(yuǎn)的位置進(jìn)行換道。綜上,出口匝道換道可能會(huì)對(duì)分流區(qū)內(nèi)駕駛員的心理產(chǎn)生影響,與普通路段的換道行為有明顯的差異。

    圖1 US-101采集路段示意圖

    圖2車輛換道時(shí)間-位置圖

    Fig.2Time-locationdiagramofvehicle lanechange

    1.2.2 數(shù)據(jù)修正

    NGSIM數(shù)據(jù)通過(guò)視頻拍攝的方式獲取,在拍攝過(guò)程中存在被前車遮擋、光線強(qiáng)度變化、周圍環(huán)境顏色突變以及圖像識(shí)別不準(zhǔn)確(黑色車和陰影、白色車輛和標(biāo)志標(biāo)線等現(xiàn)象無(wú)法區(qū)分)等現(xiàn)象,造成數(shù)據(jù)中有部分缺失。此外,其速度值是由縱向位置數(shù)據(jù)求一階倒數(shù)獲得,存在突變誤差,研究采取一階差分修正數(shù)據(jù)中的異常值。

    根據(jù)已有研究,滿足車輛性能和人體所能承受極限的加速度值為 -8~5m/s2[20] ,同時(shí),軌跡數(shù)據(jù)集中位移、速度和加速度相鄰幀之間的數(shù)據(jù)差值應(yīng)符合運(yùn)動(dòng)學(xué)的規(guī)律,相鄰幀的差值大約相等,即

    xt–xt-1≈xt-1-xt-2.

    式中: xt 為當(dāng)前幀軌跡數(shù)據(jù); xt-1 為前1幀軌跡數(shù)據(jù); xt-2 為前2幀軌跡數(shù)據(jù)。

    利用前1幀和前2幀的數(shù)據(jù)差值來(lái)估算當(dāng)前幀的軌跡數(shù)據(jù)值,計(jì)算方式如式(2)所示。

    xt*≈xt-1+(xt-1-xt-2),

    式中, xt* 為當(dāng)前幀軌跡數(shù)據(jù)的預(yù)測(cè)值。

    根據(jù)確定的加速度范圍確定誤差范圍 Δx ,當(dāng)預(yù)測(cè)值和采樣值的差值超出誤差范圍 Δx 時(shí),確定當(dāng)前時(shí)間幀的軌跡數(shù)據(jù)采樣值為異常值,然后將 xt* 重新替換為 x 要規(guī)避異常值為預(yù)測(cè)起始點(diǎn)的情況,即要滿足式(3)。

    當(dāng)連續(xù)出現(xiàn)3個(gè)異常值時(shí),應(yīng)重新選取起始點(diǎn)進(jìn)行預(yù)測(cè)。以104號(hào)車的速度數(shù)據(jù)為例演示異常值修正效果,如圖3所示。

    圖3104號(hào)車速度修正效果圖

    Fig.3Speed correction effectofcar104

    圖3中,藍(lán)色線條表示的是車輛的原始軌跡數(shù)據(jù),紅色線條為車輛修正突變誤差之后的軌跡數(shù)據(jù)。從圖中可以看到被黃色圓圈標(biāo)注的異常值已經(jīng)被剔除并替換,正常值未發(fā)生變化,修復(fù)后的速度數(shù)據(jù)與前后速度變化趨勢(shì)一致。數(shù)據(jù)從視頻中提取和誤差修正時(shí),技術(shù)處理會(huì)帶來(lái)一定的誤差,造成與真實(shí)數(shù)據(jù)差距較大,因此,還需將測(cè)量誤差和數(shù)據(jù)處理誤差應(yīng)用濾波方法進(jìn)行去噪處理。

    2考慮換道壓力增益的跟馳模型

    2.1 換道壓力增益因子構(gòu)建

    換道壓力[211是指外部因素對(duì)車輛換道行為的影響程度,與離最遲換道位置的距離有關(guān),越靠近最遲換道位置壓力越大,車輛換道的概率越大。

    為研究高速公路分流區(qū)出口匝道對(duì)換道行為的影響,提取US-101路段中92輛車換道起始點(diǎn)距出口匝道鼻端的距離。繪制換道起始點(diǎn)距出口匝道鼻端距離 L 的分布頻率圖,如圖4所示。經(jīng)擬合發(fā)現(xiàn)換道起始點(diǎn)距出口匝道鼻端距離 L 的頻率分布圖符合正態(tài)分布曲線,取置信區(qū)間為 95% ,可得出 L 的分布區(qū)間為 120~344m 。定義換道壓力增益因子 K(Ln-1) 衡量換道壓力對(duì)駕駛員心理的影響,即

    式中: μ 為換道壓力敏感系數(shù),取值范圍為 [0,1];Ln-1 為換道車輛換道位置距出口匝道鼻端的距離; Lmax 和 Lmin 為換道車輛換道位置距出口匝道鼻端的最遠(yuǎn)距離和最近距離。如圖5所示, Lmax 為 344m,Lmin 為 120m ,后文將 K(Ln-1) 簡(jiǎn)化為 K(?) 。

    圖4換道起始點(diǎn)距離出口匝道鼻端距離分布頻率直方圖

    Fig.4 Frequency of the distance from the start of the lane change to the nose end of the exit ramp

    圖5換道車輛距出口匝道鼻端距離

    Fig.5Distance from the nose of the exit ramp for lane changevehicles

    2.2 改進(jìn)跟馳模型

    全速度差模型(fullvelocity difference model,F(xiàn)VDM)[22]是在廣義力模型的基礎(chǔ)上,考慮了前后車正、負(fù)速度差對(duì)駕駛員的影響所建立的應(yīng)用廣泛的跟馳模型之一。換道車輛到達(dá)出口匝道影響范圍之前車輛的行駛采用FVDM,根據(jù)前后2輛車的速度差和間距差調(diào)整自身加速度,改變自身的駕駛狀態(tài)。

    當(dāng)前車行駛到出口匝道影響范圍時(shí),可能發(fā)生換道駛出行為,會(huì)刺激后車集中注意力,有意識(shí)的去調(diào)整駕駛狀態(tài),跟馳車輛的駕駛行為受前車和后車的間距差 Δxn,n-1(t), 第 n 輛車在 t 時(shí)刻與前車的速度差Δνn,n-1(t) 及 Ln-1 的影響。結(jié)合式FVDM的基本形式,建立考慮換道壓力增益的全速度差跟馳模型(pressure-full velocity difference model,P-FVDM)

    α=1/τ,

    κ=λ/τ

    速度優(yōu)化函數(shù)可表示為

    式中: 為跟馳車輛 n 雞在難 t 時(shí)刻的加速度, m/s2;a 為間距差敏感系數(shù);τ為駕駛員反應(yīng)時(shí)間, s;νn(t) 為第 n 輛車在 t 時(shí)刻的速度, m/s;κ 為速度差敏感系數(shù);λ為敏感系數(shù)常量; lc 為P-FVDM中車輛的最小安全間距, m ;C1 和 C2 為常數(shù)項(xiàng)。

    前車行駛一定距離后,開(kāi)始發(fā)生換道駛出行為,此時(shí)轉(zhuǎn)向燈亮起,開(kāi)始發(fā)生橫向偏移,如圖5所示。此時(shí)跟馳車輛的駕駛行為受 橫向偏移 ΔWn,n+1 及后車距出口匝道鼻端的距離 Ln+1 的影響。應(yīng)用雙車道跟馳模型理論,在FVDM的基礎(chǔ)之上,考慮前車駛離當(dāng)前車道時(shí)產(chǎn)生的橫向偏移對(duì)后車的跟馳行為的影響,王玄金[23]建立換道全速度差模型(lane pass-fullvelocity diference model,LP-FVDM)。結(jié)合LP-FVDM的基本形式,建立考慮換道壓力增益和換道駛出間距增益的全速度差模型(pressure lane pass-full velocitydifference model,PLP-FVDM)

    速度優(yōu)化函數(shù)可表示為

    式中: ρ 為橫向間距敏感系數(shù); G(?) 為橫向間距增益因子; ΔWmax 為車道寬度, m 。

    圖6橫向偏移示意圖

    Fig.6Schematicdiagramoflateraloffset

    2.3 線性穩(wěn)定性分析

    假設(shè)在初始狀態(tài)下,車隊(duì)的車頭間距都相同為 p,m ,車輛運(yùn)行速度為優(yōu)化速度 V((1+G(?)+K(?))p) 則得到各車輛的初始位置為

    假設(shè)車輛 n 在 t 時(shí)刻由于小的擾動(dòng),產(chǎn)生偏離,定義為 yn(t) ,有

    xn({})=xn0(t)+yn(t)

    將式(11)和式(12)帶人到式(10)中,可得

    式中

    Δyn(t)=yn-1(t)-yn(t)

    將式(17)帶人到式(15)中,并將其進(jìn)行展開(kāi),得到關(guān)于 z 的方程

    z2=α[(1+G(?)+K(?))V((1+G(?)+K(?))(e-1)-z]+κ(e-1)z0.

    令:

    z=z1(iκ)+z2(iκ)2+…

    將式(19)帶入式(18)中,計(jì)算得 iκ 的一階、二階項(xiàng)的系數(shù)分別為

    當(dāng) z2gt;0 時(shí),均勻交通流會(huì)逐漸演變?yōu)榉€(wěn)定狀態(tài);當(dāng) z2lt;0 時(shí),均勻交通流會(huì)逐漸演變?yōu)椴环€(wěn)定狀態(tài)。當(dāng)z2=0 時(shí)擾動(dòng)衰減,此時(shí)可得模型的中性穩(wěn)定條件為

    為更好的分析 K(?) 對(duì)模型穩(wěn)定性的影響,分別令 λ=0.5,G(?)=0.3 做出壓力增益因子取不同值時(shí)的臨界穩(wěn)定曲線圖,如圖7所示。

    圖7 λ=0.5,G(?)=0.3 的穩(wěn)定性臨界圖

    Fig.7Stabilitycriticalitydiagramfor λ=0.5 and δG(?)=0.3

    圖7中,隨著 K(?) 的增大,臨界曲線向左上方移動(dòng),反應(yīng)系數(shù)的最大值升高,穩(wěn)定區(qū)域范圍變小,此時(shí)交通流能夠在相對(duì)較小的車頭間距下穩(wěn)定運(yùn)行,但要求駕駛員的反應(yīng)時(shí)間足夠小。在臨界曲線附近更容易出現(xiàn)交通流不穩(wěn)定的現(xiàn)象,表明出口匝道對(duì)交通流的運(yùn)行有較大影響,在分流區(qū)內(nèi)交通流的穩(wěn)定性下降,換道駛出的車輛對(duì)主路車流的影響較大,極易引發(fā)交通震蕩不斷向上游傳播。

    當(dāng) G(?)=0 時(shí),式(19)退化為P-FVDM的穩(wěn)定性條件,同樣可以得出隨著 K(?) 值的增大,模型的穩(wěn)定區(qū)域變小,對(duì)駕駛員反應(yīng)能力的要求更高,更易引發(fā)交通震蕩。

    3改進(jìn)跟馳模型精度驗(yàn)證

    3.1模型參數(shù)標(biāo)定

    利用前文得到的換道車輛和其跟馳車輛的軌跡數(shù)據(jù)對(duì)改進(jìn)模型的參數(shù)進(jìn)行標(biāo)定,結(jié)合基礎(chǔ)模型的參數(shù)值確定標(biāo)定參數(shù)的取值范圍,見(jiàn)表1所示。標(biāo)定方法選取遺傳算法,遺傳算法的參數(shù)設(shè)置見(jiàn)表2。

    表1參數(shù)取值范圍Table1 Parametersrangeofvalues

    表2遺傳算法參數(shù)

    Table2 Parametersof thegeneticalgorithm

    研究選取的性能指標(biāo)為跟馳車輛的加速度,使用RMSE作為擬合優(yōu)度函數(shù)

    式中: Xobs,i 為實(shí)際值; Xmodel,i 為模型計(jì)算值; N 為樣本量。

    經(jīng)GA參數(shù)標(biāo)定后,P-FVDM和PLP-FVDM的模型參數(shù)的均值、中值、最大值、最小值和標(biāo)準(zhǔn)差分別如表3和表4所示。

    Table3 CalibrationresultsoftheP-FVDMmodel

    表4PLP-FVDM標(biāo)定結(jié)果

    表3P-FVDM標(biāo)定結(jié)果

    Table4 CalibrationresultsofPLP-FVDMmodel

    3.2 模型仿真

    為了驗(yàn)證P-FVDM和PLP-FVDM的效果,選取車輛ID為13、1226共2輛車的軌跡數(shù)據(jù),利用P-FVDM和PLP-FVDM對(duì)其進(jìn)行速度和位置的仿真,與實(shí)際軌跡數(shù)據(jù)對(duì)比分析。用原始軌跡中的初始車輛速度和位置,利用模型計(jì)算出加速度,由此可以計(jì)算出車輛的仿真速度和位置。速度仿真效果和位置仿真效果分別如圖8和9所示。

    圖8車輛速度仿真效果圖Fig.8Vehicle speed simulation

    由圖8可知,P-FVDM與PLP-FVDM對(duì)車輛速度的模擬效果良好,其中1226號(hào)車的速度仿真誤差在 1m/s 范圍內(nèi),13號(hào)車的速度仿真誤差在 1.5m/s 范圍內(nèi)。綜上,P-FVDM和PLP-FVDM速度仿真誤差都在 1.5m/s 以內(nèi),誤差范圍較小,雖然當(dāng)車輛速度突變的程度較大時(shí)仿真效果相對(duì)較差,但總體與車輛的速度變化趨勢(shì)是相符的,能夠再現(xiàn)分流區(qū)真實(shí)交通流的速度信息。

    圖9車輛位置仿真效果圖

    Fig.9 Vehiclepositionsimulation

    由圖9可知,1226號(hào)車位置仿真誤差范圍在 3m 內(nèi);13號(hào)車的位置仿真誤差范圍在 4m 內(nèi)。因此,模型對(duì)車輛位置的仿真效果總體較好,總體位置仿真誤差范圍在 4m 內(nèi),P-FVDM和PLP-FVDM能夠較好的預(yù)測(cè)車輛的位置,能夠再現(xiàn)真實(shí)交通流中車輛的位置信息,進(jìn)一步驗(yàn)證了模型的準(zhǔn)確性。

    3.3 模型有效性驗(yàn)證

    為檢驗(yàn)?zāi)P偷挠行?,?duì)比P-FVDM和FVDM、PLP-FVDM和LP-FVDM的速度仿真效果,分別計(jì)算13、48、100、359、703和1226號(hào)車輛改進(jìn)前和改進(jìn)后模型的仿真速度的均方根百分比誤差值,即RMSPE,計(jì)算公式為式(24),結(jié)果分別如圖10和11所示。

    由圖10可知,與FVDM相比,在換道車輛進(jìn)入到出口匝道影響區(qū)域時(shí),P-FVDM所控制的跟馳車輛的仿真誤差比FVDM的誤差減小,其中703號(hào)車的仿真誤差減小的最多(約為 32% ),359號(hào)車的仿真誤差減小的最少(約為 8% ),模型的總體提升約 16% 。因此,P-FVDM能夠更加精準(zhǔn)的反映換道車輛進(jìn)入出口匝道影響區(qū)發(fā)生換道行為之前主路跟馳車輛的駕駛情況。由圖11可知,與LP-FVDM相比,在換道車輛發(fā)生換道行為過(guò)程中,PLP-FVDM所控制的跟馳車輛的速度仿真誤差與真實(shí)值的誤差有所減小,模型的精度總體提升約12% ,能夠更加精準(zhǔn)的反映換道車輛換道過(guò)程中主路跟馳車輛的駕駛情況。

    圖10P-FVDM效果對(duì)比

    Fig.10 Comparison of the effect of P-FVDM model

    圖11PLP-FVDM效果對(duì)比 Fig.11 Comparison of the effect of PLP-FVDM model

    綜上,提出的模型對(duì)分流區(qū)交通流的模擬效果更好。

    4交通震蕩仿真研究

    4.1 仿真條件

    采用城市交通模擬系統(tǒng)(simulation of urban mobility,SUMO)[24]仿真軟件模擬換道車輛駛出時(shí)的分流區(qū)內(nèi)交通震蕩,初始狀態(tài)設(shè)定換道車輛后面由20輛車組成的車隊(duì)跟馳行駛,所有車輛之間的跟馳模型為FVDM,給定初始速度和初始間距使車隊(duì)行駛至穩(wěn)定狀態(tài)。

    t=0 s時(shí),換道駛出車輛到達(dá)出口匝道影響區(qū)域的起點(diǎn),換道駛出車輛和車隊(duì)第一輛車之間的控制模型改為P-FVDM,車隊(duì)內(nèi)部的車輛控制模型始終為FVDM,保持不變。在 t=21 s時(shí),換道駛出車輛開(kāi)始換道,發(fā)生橫向偏移,換道時(shí)間為 6s 。此時(shí)換道駛出車輛和車隊(duì)第一輛車之間的控制模型改變?yōu)镻LP-FVDM。換道車輛沿道路縱向的行駛速度不變,橫向偏移的方式由車輛的換道時(shí)間和車道寬度( 3.66m 計(jì)算得知,在前2.5s內(nèi)車輛勻加速,速度達(dá)到最大時(shí)車輛橫向位置到達(dá)主路車道和減速車道的車道線,在后 2.5s 內(nèi)車輛橫向速度開(kāi)始均勻減速,當(dāng)橫向速度減為 0m/s 時(shí),車輛到達(dá)減速車道的中心線,此時(shí)完成換道行為。在 時(shí),換道駛出車輛完成換道,此時(shí)對(duì)車隊(duì)的第一輛車不產(chǎn)生影響,車隊(duì)第一輛車以當(dāng)前速度繼續(xù)行駛。

    選取速度序列的標(biāo)準(zhǔn)差作為衡量震蕩振幅的指標(biāo)

    式中: u(t) 為當(dāng)前車輛的速度采集值; um 為速度序列的平均值; n 為一輛車的仿真數(shù)據(jù)點(diǎn)數(shù); In 為仿真的時(shí)間序列。

    4.2交通流速度對(duì)交通震蕩的影響

    為研究不同交通流速度對(duì)交通震蕩的影響,結(jié)合高速公路出口匝道附近車輛運(yùn)行特性,設(shè)定換道駛出車輛和車隊(duì)的運(yùn)行速度為 30,40,50km/h ,初始車頭間距由FVDM計(jì)算得出,如表5所示。為了使交通震蕩現(xiàn)象更加便于觀察,繪制出了車隊(duì)中車輛的速度標(biāo)準(zhǔn)差分布圖,仿真結(jié)果如圖12所示。

    圖12車隊(duì)速度標(biāo)準(zhǔn)差圖

    Fig.12Diagramof standarddeviationofvehiclefleetspeed

    由圖12可知,車隊(duì)以 30km/h 勻速行駛時(shí),車流內(nèi)部產(chǎn)生擾動(dòng),隨著震蕩沿車隊(duì)向上游傳播,車隊(duì)中車輛的速度標(biāo)準(zhǔn)差呈現(xiàn)凹增長(zhǎng)的模式。位于車隊(duì)上游的車輛速度標(biāo)準(zhǔn)差較大,即速度變化較大,震蕩振幅增大,震蕩現(xiàn)象更明顯,且隨著傳播過(guò)程快速增長(zhǎng)。隨著震蕩不斷向上游傳播,車輛的速度標(biāo)準(zhǔn)差的增長(zhǎng)幅度逐漸減緩,甚至呈現(xiàn)降低趨勢(shì),即此時(shí)震蕩的增長(zhǎng)開(kāi)始趨于飽和。當(dāng)車隊(duì)以 40km/h 勻速行駛時(shí),所有車輛的速度標(biāo)準(zhǔn)差非常小,均小于0.5,且沿車隊(duì)向上游傳播過(guò)程中迅速減小,即交通震蕩在傳播過(guò)程中逐漸消散。當(dāng)車隊(duì)以 50km/h 勻速行駛時(shí),車輛速度標(biāo)準(zhǔn)差更小,均小于0.15,且降低速度更快,即交通震蕩消散的更快。

    車隊(duì)以 30km/h 勻速行駛時(shí),交通流內(nèi)部不穩(wěn)定和換道車輛的換道行為共同作用,產(chǎn)生交通震蕩。車隊(duì)分別以 40,50km/h 勻速行駛時(shí),換道車輛的換道行為引起了交通震蕩,但此時(shí)車輛的車頭間距較大,前方出現(xiàn)擾動(dòng)時(shí)后方駕駛員有更多的時(shí)間對(duì)車輛進(jìn)行調(diào)控,因此,交通震蕩在向上游傳播的過(guò)程中逐漸消散。

    綜上所述,車輛換道行為能夠引發(fā)交通震蕩,但是當(dāng)行駛速度較高時(shí)震蕩會(huì)衰減,當(dāng)行駛速度較低時(shí)速度標(biāo)準(zhǔn)差沿車隊(duì)向上游呈凹增長(zhǎng)趨勢(shì)增長(zhǎng)。

    4.3換道位置對(duì)交通震蕩的影響

    初始條件下,讓換道駛出車輛和跟馳車隊(duì)均以 50km/h 的速度勻速行駛,車頭間距為 28.5m ,所有車輛均以FVDM跟馳行駛,在到達(dá)出口匝道影響區(qū)域之前使交通流達(dá)到穩(wěn)定狀態(tài)。 t=0s 時(shí),換道駛出車輛到達(dá)匝道影響區(qū)域的起點(diǎn),設(shè)定換道駛出車輛的加速度的值為 -0.2m/s2 ,分別在 s和 t=25 s開(kāi)始換道,換道時(shí)間均為6s,其余仿真設(shè)置與4.2節(jié)設(shè)置相同。

    車隊(duì)車輛的速度標(biāo)準(zhǔn)差仿真結(jié)果如圖13所示。換道時(shí)刻越晚,表示換道車輛的換道起始點(diǎn)的位置距出口匝道鼻端的距離越近。由圖13可以看出,換道車輛在不同的位置換道駛出時(shí),均發(fā)生了明顯的交通震蕩現(xiàn)象,并向車隊(duì)上游不斷傳播。分析發(fā)現(xiàn),換道車輛在不同位置換道駛出時(shí),車輛速度標(biāo)準(zhǔn)差都以凹增長(zhǎng)的方式增加,與以往研究相符。隨著換道起始點(diǎn)距出口匝道鼻端距離減小,換道車輛的換道行為對(duì)后方跟馳車輛形成的影響越大,同一輛車的速度標(biāo)準(zhǔn)差增大,沿車隊(duì)向上游傳播的速度標(biāo)準(zhǔn)差的增加幅度也隨之增大。即換道車輛越早發(fā)生換道駛出行為對(duì)交通流的穩(wěn)定越有利。

    圖13車隊(duì)速度標(biāo)準(zhǔn)差圖

    Fig.13Diagramofstandarddeviationofvehiclefleetspeed

    圖14車隊(duì)速度標(biāo)準(zhǔn)差圖

    Fig.14Diagramof standard deviationofvehicle fleetspeed

    4.4換道時(shí)長(zhǎng)對(duì)交通震蕩的影響分析

    初始條件下,讓換道駛出車輛和跟馳車隊(duì)均以 50km/h 的速度勻速行駛,車頭間距為 28.5m ,所有車輛均以FVDM跟馳,達(dá)到穩(wěn)定狀態(tài)。

    時(shí),換道駛出車輛到達(dá)出口匝道影響區(qū)域的起點(diǎn),設(shè)定換道駛出車輛的加速度的值為 -0.2m/s2 ,在 時(shí)開(kāi)始換道,產(chǎn)生橫向位移,分別在 s和 時(shí)完成換道,即換道時(shí)間分別為5s,7s和9s,其余仿真設(shè)置與文中設(shè)置相同。車隊(duì)中車輛的速度標(biāo)準(zhǔn)差仿真結(jié)果如圖14所示。

    從圖14中可以發(fā)現(xiàn),換道時(shí)間越長(zhǎng),換道行為對(duì)后方車輛的影響越大,速度標(biāo)準(zhǔn)差沿車隊(duì)向上游傳播增長(zhǎng)的越快,同時(shí)單個(gè)車輛的速度標(biāo)準(zhǔn)差更大,震蕩更加明顯。因此,換道時(shí)間越短對(duì)交通流的穩(wěn)定性越有利。

    5結(jié)束語(yǔ)

    以NGSIM數(shù)據(jù)集的US-101路段數(shù)據(jù)為數(shù)據(jù)源,對(duì)數(shù)據(jù)處理后發(fā)現(xiàn),分流區(qū)內(nèi)換道駛出車輛的換道時(shí)間范圍為3~12s,換道位置距出口匝道鼻端的距離為 50~350m ,換道位置和換道時(shí)間之間具有一定的關(guān)系。

    通過(guò)構(gòu)建高速公路分流區(qū)換道壓力增益因子對(duì)模型進(jìn)行改進(jìn),并進(jìn)行穩(wěn)定性分析,發(fā)現(xiàn)換道壓力敏感系數(shù)越高,模型的穩(wěn)定區(qū)域越小。利用重構(gòu)后的軌跡數(shù)據(jù)對(duì)模型進(jìn)行參數(shù)標(biāo)定,仿真分析發(fā)現(xiàn),改進(jìn)后的模型能夠更好地模擬車輛的速度和位置,P-FVDM的仿真精度提升了 16% ,PLP-FVDM的仿真精度提升了 12% 。利用改進(jìn)后的模型對(duì)高速公路分流區(qū)內(nèi)的交通震蕩進(jìn)行仿真,研究結(jié)果表明,車速為 30km/h 時(shí)產(chǎn)生明顯的交通震蕩現(xiàn)象,并呈凹增長(zhǎng)趨勢(shì),車速為40和 50km/h 時(shí)交通震蕩在向后傳播過(guò)程中逐漸消散;隨著換道位置距離出口匝道鼻端距離變小、換道時(shí)間變長(zhǎng),速度標(biāo)準(zhǔn)差更大,即交通震蕩振幅更大,交通流的穩(wěn)定性更差。

    研究采用的數(shù)據(jù)集來(lái)自美國(guó),不同國(guó)家的交通環(huán)境和人文因素相差較大,因此,交通流也存在一定差異,提出的模型仍需要中國(guó)的交通流數(shù)據(jù)進(jìn)行進(jìn)一步的驗(yàn)證。再者,研究只考慮了分流區(qū)內(nèi)一輛車換道駛出對(duì)主路交通流的影響,后續(xù)可考慮多輛車連續(xù)或間斷駛出時(shí)對(duì)主路交通流的影響。

    參考文獻(xiàn)

    1」Euie LU.Car-onowmg anu steauy-state meory ior noncongested trac[J].Uperauons Kesearcn,1yol,y(1): o-/0.

    [2]YeoH,SkabardonisA.Understandingstop-and-gotraficinviewofasymmetrictraffic theoryM]/TransportationandTraffic Theory 2009:Golden Jubilee.Boston,MA: Springer US,2009:99-115.

    [3]KernerBS,Rehborn H.Experimentalfeaturesandcharacteristicsof traicjams[J].Physical ReviewE,StatisticalPhysics, Plasmas,F(xiàn)luids,and Related Interdisciplinary Topics,1996,53(2): R1297-R1300.

    [4]Ahn S,CassidyMJ.Freeway traicoscilationsandvehiclelane-change maneuversC]//rocedingsofthe7thIntetioal SymposiumonTraficandTransportationTheory.July23-25,20o7,London,UK.Amsterdam:Elsevier,20o7:691-710.

    [5]Moridpour S,Rose G,SarviM.Modelling the heavy vehicle drivers’lanechanging decisionunder heavy traficconditios[J]. Road amp; Transport Research,2009,18(4):49-57.

    [6]張穎達(dá),邵春福,李慧軒,等.基于NGSIM軌跡數(shù)據(jù)的換道行為微觀特性分析[J].交通信息與安全,2015,33(6):19-24,32. Zhang YD,ShaoCF,LiHX,et al. Microscopic characteristicsof lane-change maneuvers based onNGSIM[J].Joualof Transport Information and Safety, 2015,33(6): 19-24,32.(in Chinese)

    [7]U.S.DepartmentofTransportation,F(xiàn)ederalHighwayAdministration(FHWA).Nextgenerationsimulationfactsheet[B/OL]. (2011-05-27)[2023-12-06]. htp://ps.fhwa.dot.gov/trafficanalysistools/ngsim.htm.

    [8]SusilawatiS,TanKY,SamadKMA,etal.Thestudyoflocationspecificlanechangeimpactsotraficdelayusingetended cell transmissonmodel[C]/Proceedingsof the2018EUREKA International ConferenceonEnginering,Technologyand Construction (EUREKA2018).April24-26,2018,Kuching,Sarawak,Malaysia.Les Ulis:EDP Sciences,2018,203:05008. DOI: 10.1051/matecconf/201820305008.

    [9]LiXP,CuiJX,AnS,etal.Stop-and-gotrafcanalysis:theoreticalproperties,environmentalimpactsandoilation mitigation[J].Transportation Research Part B:Methodological, 2014,70: 319-339.

    [10]鄭淑欣,李振龍,趙曉華.霧天高速公路換道行為特性研究[J].交通信息與安全,2020,38(1):35-42,51. Zheng SX,LiZL,Zhao XH.Behavioral characteristicsof lane-changingon highwayin foggyweather[J].Journalof Transport Information and Safety, 2020,38(1): 35-42, 51.(in Chinese)

    [11]Zheng ST,JiangR,TianJF,etal.Empiricalandexperimentalsudyontegowthpaternof taicoscilationsupstramof fixed bottleneck and model test[J].Transportation Research Part C:Emerging Technologies,2022,140:103729.

    [12]ToledoT,ZoharD.Modelingdurationof lanechanges[J].Transportation Research Record:Journalof theTransportation Research Board,2007,1999(1): 71-78.

    [13]Jiang R,Hu MB,Zhang HM,et al.Onsome experimental featuresof car-folowing behaviorand how to model them[J]. Transportation Research Part B:Methodological, 2015,80:338-354.

    [14]JiangR,JinCJ,ZhangHM,etal.ExperimentalandmpiricalinvestigatiooftraicfloistabilityCProeedingoft 6th International Symposium of the Transport and Trafic ScienceForum (HEART2016).May17-20,2016,Xi'an,China. Amsterdam: Elsevier,2017,23:157-173.DOI: 10.1016/j.trpro.2017.05.011.

    [15]賀玉龍,劉磊,遲佳欣.高速公路車輛換道行為風(fēng)險(xiǎn)研究[J].重慶交通大學(xué)學(xué)報(bào)(自然科學(xué)版),2021,40(4):26-33. HeYL,LiuL,ChiJX.Risk of vehiclelanechange behavioronhighway[J].JournalofChongqing Jiaotong University (Natural Science),2021,40(4): 26-33.(in Chinese)

    [16]Mauch ,CassidyM.Freeway traffc:oscilltions observations and predictionsM]/Transportationand TraffcTheorynthe 21stCentury.Amsterdam:Elsevier,2002: 653-673.

    [17]張娟.城市快速路瓶頸交通震蕩特性研究[J].交通信息與安全,2015,33(5):43-48. ZhangJ.Astudyoftraffcoscillationsat thebottlenecksonurbanexpressways[J].JouralofTransportInformationandSafety 2015,33(5): 43-48.(in Chinese)

    [18]黃永現(xiàn).交通震蕩演化特性的實(shí)驗(yàn)研究和建模仿真[D].合肥:中國(guó)科學(xué)技術(shù)大學(xué),2019. Huang YX.Experimental research and modelingontheevolution of trafic oscilation[D]. Hefei: UniversityofScience and Technology of China,2019.(in Chinese)

    [19]YaoHD,LiQWLiXP.Astudyofrelationshipsintraicoscilltonfeaturesbasdonfieldexperiments[J].Trasprtation ResearchPartA:Policyand Practice,2020,141:339-355.

    [20]馬艷鋒.基于NGSIM數(shù)據(jù)的高速公路出入口強(qiáng)制換道車輛執(zhí)行及后調(diào)整行為研究[D].西安:長(zhǎng)安大學(xué),2020. MaYF.Research onvehicles’mandatorylane change behavior in executionand post-adjustment periodonmulti-lane freeway based on NGSIM trajectory data[D]. Xi'an; Chang'an University,2020.(in Chinese)

    [21]劉有軍,曹珊.基于元胞自動(dòng)機(jī)的強(qiáng)制換道模型研究[J].交通信息與安全,2009,27(3):78-80. LiuYJ,Cao S.Compulsorylane-changing traficmodel based on celularautomaton[J].Journalof Tansport Informationand Safety,2009,27(3):78-80.(in Chinese)

    [22]JiangR,WuQS,ZhuZJ.Fullvelocitydierence modelforacar-following heory[J].PhysicalReviewE,2001,64017101.

    [23]王玄金.考慮前車換道過(guò)程的跟馳模型研究[D].重慶:重慶大學(xué),2017. Wang X J.Research on car following model considering preceding vehicle’slane-changing processD].Chongqing: Chongqing University,2017. (in Chinese)

    [24]Barcel6J.FundamentalsoftraficsimulationM].NewYork:Springer,20l0.(InternationalSeriesiOperations Researchamp; Management Science,vol 145).

    (編輯 呂建斌)

    猜你喜歡
    交通流匝道標(biāo)準(zhǔn)差
    基于車聯(lián)網(wǎng)的中心導(dǎo)航云性能分析模型
    中國(guó)綠色化發(fā)展水平的時(shí)空格局及優(yōu)化策略
    智慧收費(fèi)系統(tǒng)在全椒西收費(fèi)站的應(yīng)用
    袋栽鹿茸菇菌棒物理力學(xué)特性研究
    芻議教材“二次開(kāi)發(fā)”的若干視角
    免费久久久久久久精品成人欧美视频 | 久久精品熟女亚洲av麻豆精品| 一本久久精品| 嫩草影院入口| 日本av手机在线免费观看| 看非洲黑人一级黄片| 哪个播放器可以免费观看大片| 亚洲精品日韩在线中文字幕| 两个人的视频大全免费| 免费av不卡在线播放| 亚洲国产日韩一区二区| 久久久a久久爽久久v久久| 中文字幕人妻熟人妻熟丝袜美| 中文字幕制服av| 九九爱精品视频在线观看| 国产白丝娇喘喷水9色精品| 亚洲成人av在线免费| 免费人妻精品一区二区三区视频| 亚洲内射少妇av| 成人无遮挡网站| 亚洲欧美一区二区三区国产| 国产一区有黄有色的免费视频| 高清视频免费观看一区二区| 精品人妻视频免费看| 国产男人的电影天堂91| 五月天丁香电影| 中文字幕人妻熟人妻熟丝袜美| 中文资源天堂在线| 国产精品福利在线免费观看| 午夜免费鲁丝| 特大巨黑吊av在线直播| www.av在线官网国产| 亚洲精品日韩在线中文字幕| 欧美少妇被猛烈插入视频| 在现免费观看毛片| 精品久久久噜噜| 国产av码专区亚洲av| 欧美bdsm另类| 女人十人毛片免费观看3o分钟| 国产精品一区www在线观看| 亚洲国产精品国产精品| 欧美最新免费一区二区三区| 男的添女的下面高潮视频| 久久久久久九九精品二区国产| 在线天堂最新版资源| 你懂的网址亚洲精品在线观看| 欧美日韩综合久久久久久| 欧美区成人在线视频| 九九在线视频观看精品| 亚洲精品自拍成人| 日韩精品有码人妻一区| 国产亚洲欧美精品永久| 91久久精品国产一区二区三区| 亚洲欧洲日产国产| 激情 狠狠 欧美| 一级二级三级毛片免费看| 精品午夜福利在线看| 亚洲国产高清在线一区二区三| 午夜福利在线观看免费完整高清在| 亚洲丝袜综合中文字幕| 少妇精品久久久久久久| 亚洲精品日本国产第一区| 一区二区三区四区激情视频| 欧美日韩国产mv在线观看视频 | 青青草视频在线视频观看| 激情五月婷婷亚洲| 精品久久国产蜜桃| 久久久久人妻精品一区果冻| 天堂8中文在线网| 亚洲精品国产成人久久av| 久久国产亚洲av麻豆专区| 一级毛片我不卡| 简卡轻食公司| 久久久久久久国产电影| 精品一区二区三卡| 在线精品无人区一区二区三 | 国产免费视频播放在线视频| 老师上课跳d突然被开到最大视频| 精品少妇黑人巨大在线播放| 26uuu在线亚洲综合色| 在线观看免费高清a一片| 成年av动漫网址| 三级经典国产精品| 色婷婷av一区二区三区视频| a 毛片基地| 久热这里只有精品99| 22中文网久久字幕| 校园人妻丝袜中文字幕| 亚洲美女搞黄在线观看| 制服丝袜香蕉在线| 久久久精品94久久精品| 国产精品久久久久久久久免| 久久女婷五月综合色啪小说| 最近最新中文字幕免费大全7| 深爱激情五月婷婷| 亚洲精品乱码久久久v下载方式| 精品酒店卫生间| 亚洲怡红院男人天堂| videossex国产| 又黄又爽又刺激的免费视频.| 九草在线视频观看| 亚洲欧美一区二区三区黑人 | 最近最新中文字幕大全电影3| 国产无遮挡羞羞视频在线观看| 51国产日韩欧美| 18禁动态无遮挡网站| 亚洲图色成人| 日本欧美视频一区| 99国产精品免费福利视频| 国产熟女欧美一区二区| 免费高清在线观看视频在线观看| 国产亚洲5aaaaa淫片| 亚洲综合色惰| 狂野欧美激情性xxxx在线观看| 校园人妻丝袜中文字幕| 国产精品99久久99久久久不卡 | 久久99蜜桃精品久久| 免费久久久久久久精品成人欧美视频 | 只有这里有精品99| 亚洲欧洲日产国产| 精品一区二区三卡| 国产成人a区在线观看| 国产精品不卡视频一区二区| 国产免费一级a男人的天堂| 久久久成人免费电影| 校园人妻丝袜中文字幕| 人妻少妇偷人精品九色| 国产深夜福利视频在线观看| 日韩一区二区视频免费看| 亚洲性久久影院| 身体一侧抽搐| 3wmmmm亚洲av在线观看| 欧美少妇被猛烈插入视频| 我要看黄色一级片免费的| 内地一区二区视频在线| 亚洲人与动物交配视频| av国产久精品久网站免费入址| 亚洲四区av| 国产爽快片一区二区三区| kizo精华| 婷婷色麻豆天堂久久| 久久久国产一区二区| 欧美zozozo另类| 99热6这里只有精品| 国产精品麻豆人妻色哟哟久久| 在线播放无遮挡| 久久热精品热| 成人国产av品久久久| 亚洲精品久久久久久婷婷小说| 中文字幕av成人在线电影| 精品久久久久久电影网| 国产精品一区二区在线不卡| 国产在视频线精品| 干丝袜人妻中文字幕| 中文字幕人妻熟人妻熟丝袜美| 国产av一区二区精品久久 | 大话2 男鬼变身卡| 国产真实伦视频高清在线观看| 婷婷色综合大香蕉| 国产国拍精品亚洲av在线观看| 观看美女的网站| 一区二区三区免费毛片| 韩国高清视频一区二区三区| 国产成人精品一,二区| 欧美激情极品国产一区二区三区 | 大话2 男鬼变身卡| 欧美日韩视频精品一区| 中文精品一卡2卡3卡4更新| 日韩成人av中文字幕在线观看| 王馨瑶露胸无遮挡在线观看| 干丝袜人妻中文字幕| 赤兔流量卡办理| 天堂中文最新版在线下载| av国产免费在线观看| 色吧在线观看| 久久精品久久精品一区二区三区| 久久久久网色| 欧美成人午夜免费资源| 国产伦精品一区二区三区四那| 97超碰精品成人国产| 久久久久久久久久人人人人人人| 天堂中文最新版在线下载| 精品人妻视频免费看| 啦啦啦在线观看免费高清www| 深爱激情五月婷婷| 久久久久精品性色| 亚洲图色成人| 中文资源天堂在线| 永久免费av网站大全| 最新中文字幕久久久久| 国产亚洲午夜精品一区二区久久| 亚洲国产精品一区三区| 中文乱码字字幕精品一区二区三区| 亚洲欧美日韩另类电影网站 | 不卡视频在线观看欧美| 精品熟女少妇av免费看| 欧美日韩精品成人综合77777| 日韩伦理黄色片| 有码 亚洲区| 一本久久精品| 在线亚洲精品国产二区图片欧美 | 激情五月婷婷亚洲| 在现免费观看毛片| 99久久人妻综合| 不卡视频在线观看欧美| 亚洲av欧美aⅴ国产| 七月丁香在线播放| 亚洲色图综合在线观看| 最近2019中文字幕mv第一页| 在线亚洲精品国产二区图片欧美 | 啦啦啦中文免费视频观看日本| 夜夜骑夜夜射夜夜干| 91精品国产九色| 午夜免费男女啪啪视频观看| 亚洲电影在线观看av| 久久久久精品久久久久真实原创| 亚洲国产欧美人成| 男人和女人高潮做爰伦理| 亚洲人成网站在线播| 看非洲黑人一级黄片| 久久久久人妻精品一区果冻| 亚洲国产精品成人久久小说| 国产大屁股一区二区在线视频| 久久毛片免费看一区二区三区| 国产精品一区二区在线观看99| 精品人妻偷拍中文字幕| 精品一区二区免费观看| 91在线精品国自产拍蜜月| 国产人妻一区二区三区在| 国产亚洲av片在线观看秒播厂| 日本黄色日本黄色录像| 日产精品乱码卡一卡2卡三| 免费av中文字幕在线| 我要看黄色一级片免费的| 99久久人妻综合| 人人妻人人爽人人添夜夜欢视频 | 亚洲无线观看免费| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲综合色惰| 日本一二三区视频观看| 又黄又爽又刺激的免费视频.| av.在线天堂| 日本与韩国留学比较| 亚洲经典国产精华液单| 日本午夜av视频| 欧美性感艳星| 亚洲欧美一区二区三区国产| 日韩一区二区视频免费看| 成年美女黄网站色视频大全免费 | 天堂俺去俺来也www色官网| 亚洲国产精品999| 国产精品福利在线免费观看| 日本猛色少妇xxxxx猛交久久| 亚洲高清免费不卡视频| 国产精品久久久久久久久免| 91狼人影院| av播播在线观看一区| 久久久久人妻精品一区果冻| 少妇人妻精品综合一区二区| 秋霞伦理黄片| 久久久久久九九精品二区国产| 爱豆传媒免费全集在线观看| 王馨瑶露胸无遮挡在线观看| 国产在线一区二区三区精| 亚洲中文av在线| 日本vs欧美在线观看视频 | 亚洲,一卡二卡三卡| 如何舔出高潮| 亚洲,欧美,日韩| av免费观看日本| 在线观看免费高清a一片| 国产91av在线免费观看| 高清视频免费观看一区二区| 国产真实伦视频高清在线观看| 噜噜噜噜噜久久久久久91| a级毛色黄片| 亚洲性久久影院| 又大又黄又爽视频免费| 亚洲成人手机| 内地一区二区视频在线| 久久久久久久久久久免费av| 亚洲精品日韩av片在线观看| 久久久久久久精品精品| 99热这里只有精品一区| h日本视频在线播放| 99久久精品热视频| 国产成人免费无遮挡视频| 中文字幕亚洲精品专区| 人妻制服诱惑在线中文字幕| 黄片无遮挡物在线观看| 一级毛片久久久久久久久女| 嘟嘟电影网在线观看| 亚洲最大成人中文| 欧美日韩精品成人综合77777| 亚洲精品国产成人久久av| 午夜激情福利司机影院| 色哟哟·www| 黑人高潮一二区| 久久久久网色| 国产av一区二区精品久久 | 最新中文字幕久久久久| av国产精品久久久久影院| 亚洲人成网站在线播| 国产精品.久久久| 免费播放大片免费观看视频在线观看| 色视频在线一区二区三区| 欧美日韩亚洲高清精品| 麻豆成人av视频| 中文欧美无线码| 亚洲欧美精品专区久久| 天堂8中文在线网| 香蕉精品网在线| h日本视频在线播放| 色吧在线观看| 久久久久精品性色| 久久毛片免费看一区二区三区| 1000部很黄的大片| 男女边吃奶边做爰视频| 久久久久国产精品人妻一区二区| 免费少妇av软件| 国产极品天堂在线| 观看美女的网站| 国产91av在线免费观看| 高清午夜精品一区二区三区| 精品久久国产蜜桃| 99热全是精品| 亚洲人成网站在线观看播放| av在线观看视频网站免费| 91在线精品国自产拍蜜月| 亚洲美女搞黄在线观看| 亚洲国产最新在线播放| 精品人妻熟女av久视频| 国产精品人妻久久久久久| 免费人妻精品一区二区三区视频| 成人特级av手机在线观看| 妹子高潮喷水视频| 亚洲成色77777| av又黄又爽大尺度在线免费看| 91久久精品国产一区二区成人| 亚洲精品一区蜜桃| 在线看a的网站| 久久ye,这里只有精品| 如何舔出高潮| 人人妻人人澡人人爽人人夜夜| 精品国产乱码久久久久久小说| 六月丁香七月| 国产色婷婷99| 亚洲精品自拍成人| 嫩草影院入口| 亚洲国产毛片av蜜桃av| 中文天堂在线官网| 我要看日韩黄色一级片| 久热久热在线精品观看| 你懂的网址亚洲精品在线观看| 日韩中文字幕视频在线看片 | 五月开心婷婷网| 日韩大片免费观看网站| 韩国av在线不卡| 大香蕉久久网| 免费不卡的大黄色大毛片视频在线观看| 国产亚洲精品久久久com| 欧美日本视频| 草草在线视频免费看| 下体分泌物呈黄色| 国产探花极品一区二区| 国产精品秋霞免费鲁丝片| 久久久久网色| 日韩,欧美,国产一区二区三区| 黄片wwwwww| 国产精品熟女久久久久浪| 91狼人影院| 免费不卡的大黄色大毛片视频在线观看| 大片免费播放器 马上看| 欧美97在线视频| av播播在线观看一区| 成人无遮挡网站| 99久久精品一区二区三区| 精品国产三级普通话版| 91精品国产国语对白视频| 中文字幕久久专区| 日本av免费视频播放| av不卡在线播放| 尾随美女入室| 国产精品国产三级国产av玫瑰| 美女cb高潮喷水在线观看| 国产成人一区二区在线| 七月丁香在线播放| 日韩一区二区三区影片| av天堂中文字幕网| av一本久久久久| av线在线观看网站| 人妻夜夜爽99麻豆av| 久久精品久久精品一区二区三区| 纯流量卡能插随身wifi吗| 高清黄色对白视频在线免费看 | 成人特级av手机在线观看| 99久久人妻综合| 成人毛片a级毛片在线播放| 久久6这里有精品| av专区在线播放| 少妇人妻 视频| 国产 精品1| 久久鲁丝午夜福利片| 国产亚洲91精品色在线| 亚洲精品456在线播放app| 亚洲第一区二区三区不卡| 久久精品久久精品一区二区三区| tube8黄色片| 自拍偷自拍亚洲精品老妇| 美女cb高潮喷水在线观看| 麻豆成人av视频| 18禁裸乳无遮挡动漫免费视频| 99热6这里只有精品| 美女脱内裤让男人舔精品视频| 中国国产av一级| 国产日韩欧美在线精品| 国产 一区精品| 下体分泌物呈黄色| 亚洲av二区三区四区| 亚洲国产精品一区三区| 精品久久久噜噜| 99热国产这里只有精品6| 久久精品夜色国产| 欧美日韩一区二区视频在线观看视频在线| 亚洲内射少妇av| 国产乱人偷精品视频| 欧美精品人与动牲交sv欧美| 91精品国产国语对白视频| 另类亚洲欧美激情| 我的女老师完整版在线观看| 亚洲av综合色区一区| 欧美精品人与动牲交sv欧美| av国产久精品久网站免费入址| 成年女人在线观看亚洲视频| 日日撸夜夜添| 99热6这里只有精品| 欧美日韩一区二区视频在线观看视频在线| 精品久久国产蜜桃| 亚洲图色成人| 欧美精品一区二区免费开放| 国产人妻一区二区三区在| 美女脱内裤让男人舔精品视频| 青春草亚洲视频在线观看| 国产亚洲av片在线观看秒播厂| 国产精品偷伦视频观看了| 亚洲精品中文字幕在线视频 | 日韩欧美精品免费久久| 久久av网站| 在线观看三级黄色| 亚洲精品乱码久久久v下载方式| 久久人妻熟女aⅴ| 久久韩国三级中文字幕| 在线播放无遮挡| 久久人人爽人人爽人人片va| 成年av动漫网址| 只有这里有精品99| 下体分泌物呈黄色| 高清午夜精品一区二区三区| 国产成人aa在线观看| 美女脱内裤让男人舔精品视频| 亚洲欧洲日产国产| 九九爱精品视频在线观看| av网站免费在线观看视频| 日本爱情动作片www.在线观看| 中文在线观看免费www的网站| 亚洲精品视频女| 国语对白做爰xxxⅹ性视频网站| 中文在线观看免费www的网站| 日韩人妻高清精品专区| 18禁在线无遮挡免费观看视频| 国产伦精品一区二区三区四那| 日产精品乱码卡一卡2卡三| av不卡在线播放| 性高湖久久久久久久久免费观看| 一级毛片 在线播放| 日本黄色片子视频| 天堂中文最新版在线下载| 国产av国产精品国产| 久久这里有精品视频免费| 亚洲国产av新网站| 国产精品一二三区在线看| 国产高潮美女av| 国产免费视频播放在线视频| 亚洲欧美成人精品一区二区| 爱豆传媒免费全集在线观看| 国产69精品久久久久777片| 人体艺术视频欧美日本| 在线观看美女被高潮喷水网站| 免费黄频网站在线观看国产| 只有这里有精品99| 美女cb高潮喷水在线观看| 十八禁网站网址无遮挡 | 婷婷色麻豆天堂久久| 日本欧美国产在线视频| 青春草视频在线免费观看| 日韩一本色道免费dvd| 亚洲av电影在线观看一区二区三区| 人妻少妇偷人精品九色| 秋霞伦理黄片| 日本黄大片高清| 国产精品欧美亚洲77777| 久久精品国产亚洲网站| 亚洲中文av在线| 亚洲欧洲国产日韩| av国产久精品久网站免费入址| 亚洲精品456在线播放app| 久久女婷五月综合色啪小说| 日韩中文字幕视频在线看片 | 熟女av电影| 精品久久久久久电影网| 国产片特级美女逼逼视频| 欧美日本视频| 日本wwww免费看| 国产精品一区www在线观看| 国产伦在线观看视频一区| 免费播放大片免费观看视频在线观看| 午夜福利视频精品| 日本欧美视频一区| 国产黄色免费在线视频| 午夜福利视频精品| 一级毛片aaaaaa免费看小| 街头女战士在线观看网站| 亚洲成人手机| 国产成人a区在线观看| 免费人妻精品一区二区三区视频| 久久久久久久精品精品| 国产 一区 欧美 日韩| freevideosex欧美| 少妇猛男粗大的猛烈进出视频| 国产在线一区二区三区精| 男女边摸边吃奶| 在线免费十八禁| 亚洲欧美一区二区三区黑人 | 大话2 男鬼变身卡| 国语对白做爰xxxⅹ性视频网站| 国产成人freesex在线| 免费大片18禁| 国产精品不卡视频一区二区| 男女国产视频网站| a级一级毛片免费在线观看| 国产视频内射| 纵有疾风起免费观看全集完整版| 亚洲成人手机| 黄色一级大片看看| 卡戴珊不雅视频在线播放| 伦理电影大哥的女人| 卡戴珊不雅视频在线播放| 亚洲欧美日韩东京热| 欧美日韩视频高清一区二区三区二| 日本欧美视频一区| 五月玫瑰六月丁香| 国产免费又黄又爽又色| 久久久国产一区二区| 男人添女人高潮全过程视频| 久久久国产一区二区| 大又大粗又爽又黄少妇毛片口| 亚洲怡红院男人天堂| 精品人妻一区二区三区麻豆| 亚洲av男天堂| 精品人妻一区二区三区麻豆| 七月丁香在线播放| 免费黄频网站在线观看国产| 国产男人的电影天堂91| 老熟女久久久| 久久久久久伊人网av| 国产日韩欧美在线精品| 九色成人免费人妻av| 欧美一级a爱片免费观看看| 在线观看三级黄色| 青春草视频在线免费观看| 久久人人爽人人爽人人片va| 成人午夜精彩视频在线观看| a级毛片免费高清观看在线播放| 免费看日本二区| 欧美日韩精品成人综合77777| 欧美成人精品欧美一级黄| 在线观看一区二区三区| 直男gayav资源| 十分钟在线观看高清视频www | 女人十人毛片免费观看3o分钟| 国产探花极品一区二区| av在线播放精品| 久久久久久久久久人人人人人人| 亚洲,一卡二卡三卡| 人人妻人人澡人人爽人人夜夜| 欧美zozozo另类| 亚洲av二区三区四区| 亚洲精品国产av成人精品| 身体一侧抽搐| 亚洲精品成人av观看孕妇| 精品人妻视频免费看| 蜜桃久久精品国产亚洲av| 高清日韩中文字幕在线| 国产成人a区在线观看| 精品久久久久久久久av| 国产精品人妻久久久久久| 久久久久久久精品精品| 国产精品.久久久| 亚洲国产最新在线播放| 亚洲欧美一区二区三区黑人 | 日韩制服骚丝袜av| 久久久a久久爽久久v久久| 国产男女内射视频| 人人妻人人看人人澡| 少妇裸体淫交视频免费看高清| 亚洲在久久综合| 91精品伊人久久大香线蕉| 一级毛片黄色毛片免费观看视频| 亚洲综合精品二区| av网站免费在线观看视频| 干丝袜人妻中文字幕| 日韩制服骚丝袜av| 免费少妇av软件| 人妻一区二区av| 青青草视频在线视频观看|