亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        展寬段排隊(duì)溢出狀態(tài)下區(qū)段延誤模型研究

        2016-06-30 01:52:06蔡艷鋒孫大志唐貴濤
        公路交通科技 2016年6期

        蔡艷鋒,孫大志,唐貴濤

        (1.西南交通大學(xué) 交通運(yùn)輸與物流學(xué)院,四川 成都 610031;2.昆明理工大學(xué) 交通工程學(xué)院,云南 昆明 650093)

        展寬段排隊(duì)溢出狀態(tài)下區(qū)段延誤模型研究

        蔡艷鋒1,孫大志1,唐貴濤2

        (1.西南交通大學(xué)交通運(yùn)輸與物流學(xué)院,四川成都610031;2.昆明理工大學(xué)交通工程學(xué)院,云南昆明650093)

        摘要:為研究過飽和信號交叉口排隊(duì)溢出狀態(tài)下單周期通過展寬段車輛的區(qū)段延誤模型,首先介紹了3大經(jīng)典延誤模型,提出了3大延誤模型均未考慮展寬段長度和飽和流率的波動(dòng)問題。其次,以展寬段為研究對象,在經(jīng)典延誤模型理論的基礎(chǔ)上,推導(dǎo)出了展寬段區(qū)段延誤模型。最后,從公式推導(dǎo)計(jì)算和VISSIM仿真兩個(gè)方面對優(yōu)化模型和定數(shù)理論延誤模型進(jìn)行了對比。結(jié)果表明:該模型能夠有效地反映展寬段溢出狀態(tài)下產(chǎn)生的區(qū)段延誤,可為展寬段的科學(xué)設(shè)計(jì)提供理論基礎(chǔ)。

        關(guān)鍵詞:交通工程;展寬段;VISSIM仿真;區(qū)段延誤模型;排隊(duì)溢出

        0引言

        信號交叉口是城市路網(wǎng)的重要組成部分,如何提高交通流在交叉口處的通行效率一直是城市交通管控的重點(diǎn)和難點(diǎn)。拓寬進(jìn)口道可以增加進(jìn)口道車道數(shù)并提高進(jìn)口道的通行能力,其長度由展寬漸變段和展寬段構(gòu)成[1]。紅燈期間到達(dá)的車輛在展寬段內(nèi)壓縮存儲(chǔ),綠燈期間高效釋放。早在1967年,Leisch[2]將進(jìn)出口道拓寬部分相對于路段增加的車道稱為短車道,并提出了若干短車道長度建議值;Lee[3]通過模擬手段分析了拓寬進(jìn)口道的通行能力,但并未考慮溢出擁堵問題;Tian[4]、楊曉光[5]、程國柱[6]從概率論角度研究了展寬段長度受限下的交叉口通行能力。

        延誤是指行駛在路段上的車輛由于受到道路環(huán)境、交通管理與控制或其他車輛的干擾等影響而損失的時(shí)間[7]?,F(xiàn)有延誤模型多以整個(gè)進(jìn)口道上的所有車輛為研究對象,較少對空間某路段通過車輛延誤進(jìn)行研究。隨著展寬段的普遍設(shè)置,為了更好地發(fā)揮展寬段的功能,研究展寬段單周期內(nèi)通過車輛所產(chǎn)生的區(qū)段延誤極有必要。區(qū)段延誤是指單位時(shí)間內(nèi)通過某固定路段車輛所產(chǎn)生的延誤。本文在3大經(jīng)典延誤理論的基礎(chǔ)上,以整個(gè)展寬段為研究對象,在考慮展寬段最大排隊(duì)長度和飽和流率波動(dòng)的前提下,對展寬段排隊(duì)溢出時(shí)所產(chǎn)生的區(qū)段延誤進(jìn)行分析研究。

        1經(jīng)典延誤模型

        多年來,國內(nèi)外對交叉口延誤的研究積累了豐富的成果,比較經(jīng)典的理論方法有穩(wěn)態(tài)理論、定數(shù)理論和過渡函數(shù)理論。

        (1)穩(wěn)態(tài)理論

        穩(wěn)態(tài)理論研究的是低飽和度穩(wěn)態(tài)交通流延誤時(shí)間和停車次數(shù)的數(shù)學(xué)表達(dá)式。車輛平均延誤時(shí)間可以簡化為均衡相位平均延誤時(shí)間與隨機(jī)平均延誤時(shí)間兩部分,較為經(jīng)典的Webster延誤模型[8]見式(1)、式(2)。另外,Miller和Akcelik也提出了相似的延誤計(jì)算模型[9-10]。

        穩(wěn)態(tài)理論總延誤:

        (1)

        穩(wěn)態(tài)理論平均延誤:

        (2)

        式中,D為總延誤;d為平均延誤;C為信號周期時(shí)長;q為車輛到達(dá)率;g為綠燈時(shí)長;S為飽和流率;x為飽和度。

        穩(wěn)態(tài)理論在低飽和度的情況下是比較切合實(shí)際的,但隨著飽和度的增高,穩(wěn)態(tài)理論計(jì)算結(jié)果與實(shí)際情況有較大出入,因此穩(wěn)態(tài)理論不適用于計(jì)算過飽和狀態(tài)下的延誤。另外從式(1)、式(2)可以看出,穩(wěn)態(tài)理論在計(jì)算延誤時(shí)沒有考慮研究路段的長度及飽和流率的波動(dòng)問題。

        (2)定數(shù)理論

        定數(shù)理論是針對過飽和交叉口排隊(duì)車輛隨周期累積的現(xiàn)象提出的,其中典型的研究成果是May在1965年編著的《交通流理論》和其與Heller在1976年合作編著的《定數(shù)排隊(duì)模式》[11]。定數(shù)理論假設(shè)車輛到達(dá)率q大于通行能力Q,且在研究時(shí)段內(nèi)為定值。定數(shù)理論把過飽和交通流作為一種確定的情況考慮,交通流排隊(duì)長度、延誤時(shí)間隨時(shí)間線性增加。

        定數(shù)理論總延誤:

        (3)

        定數(shù)理論平均延誤:

        (4)

        式中,r為紅燈時(shí)長;N為平均滯留車輛數(shù)。

        但定數(shù)理論同樣沒有考慮研究路段的長度及飽和流率的波動(dòng)問題,因此計(jì)算出來的理論排隊(duì)長度及延誤時(shí)間是無限制的,這與本文要研究的展寬段內(nèi)車輛周期區(qū)段延誤相矛盾。

        (3)過渡函數(shù)理論

        過渡函數(shù)理論首先由Whiting提出,主要用于解決仿真程序中關(guān)于隨機(jī)延誤的模擬。后來,Kimber和Hollis深入研究了穩(wěn)態(tài)理論與定數(shù)理論的結(jié)合方法,提出了過渡函數(shù)理論的推導(dǎo)模型[12]。

        過渡函數(shù)總延誤:

        (5)

        過渡函數(shù)平均延誤:

        (6)

        過渡函數(shù)理論的建立不僅解決了近飽和交通狀態(tài)下車輛受阻滯程度的定量分析,而且彌補(bǔ)了過飽和交通狀態(tài)下被定數(shù)理論所忽略的隨機(jī)阻滯。但是同穩(wěn)態(tài)理論和定數(shù)理論一樣,過渡函數(shù)也沒有考慮路段長度及飽和流率的波動(dòng)問題,所以對過飽和狀態(tài)下展寬段區(qū)段延誤的計(jì)算不合適。

        以上3種經(jīng)典延誤模型在計(jì)算信號交叉口延誤時(shí)均沒有考慮到路段的長度及飽和流率的波動(dòng)問題,因此對于展寬段排隊(duì)溢出狀態(tài)下區(qū)段延誤的計(jì)算結(jié)果將會(huì)與實(shí)際出現(xiàn)偏差。因此,本文以3種經(jīng)典延誤模型為基礎(chǔ),推導(dǎo)適合展寬段排隊(duì)溢出狀態(tài)的區(qū)段延誤計(jì)算模型。

        2展寬段溢出狀態(tài)下區(qū)段延誤模型的建立

        本文以展寬段為研究對象,考慮在展寬段排隊(duì)溢出時(shí),單位周期通過展寬段車輛的區(qū)段總延誤及車均延誤。在這種情況下,展寬段內(nèi)排隊(duì)長度是固定的,不會(huì)像定數(shù)理論那樣隨著時(shí)間的延續(xù)一直增加。當(dāng)綠燈啟亮?xí)r,隨著隊(duì)首車輛的駛離,展寬段后排隊(duì)車輛陸續(xù)進(jìn)入展寬段。

        (1)以圖1所示交叉口為研究對象,各個(gè)方向設(shè)計(jì)參數(shù)相同,采用單口放行;

        (2)將展寬漸變段等同于展寬段進(jìn)行分析,陰影部分為研究區(qū)域;

        (3)西進(jìn)口道由于過飽和引起排隊(duì)溢出現(xiàn)象的發(fā)生;

        (5)由于設(shè)有專用車道且不受信號控制,不對右轉(zhuǎn)車輛進(jìn)行研究。

        圖1 信號交叉口Fig.1 Signalized intersection

        2.1展寬段駛?cè)肓髀蕅s和消散飽和流率S不恒定

        實(shí)際情況中,展寬段駛?cè)肓髀蕅s和消散飽和流率S為變量。根據(jù)研究區(qū)段內(nèi)交通流的特性,繪制出展寬段內(nèi)單車道排隊(duì)車輛數(shù)Qu與延誤時(shí)間t的關(guān)系曲線,如圖2所示。

        要實(shí)現(xiàn)自動(dòng)氣象站的正常運(yùn)轉(zhuǎn),在對其進(jìn)行安裝以及相關(guān)系統(tǒng)設(shè)置上要嚴(yán)格遵循相關(guān)要求進(jìn)行,以保證自動(dòng)氣象站在進(jìn)行實(shí)際測量時(shí)不會(huì)受到雷電的干擾。自動(dòng)氣象站防雷技術(shù)在使用時(shí)要結(jié)合當(dāng)?shù)氐膶?shí)際情況靈活運(yùn)用,以保證防雷效果的提高,避免自動(dòng)氣象站受損。

        圖2 qs與S不恒定時(shí)的排隊(duì)車輛數(shù)與延誤時(shí)間Fig.2 Number of queuing vehicles and delay time when qs and S are unsteady

        圖2中,n為展寬段內(nèi)單車道所能容納的最大車輛數(shù);δ為綠燈啟亮后啟動(dòng)波傳播至展寬段末端所用的時(shí)間;m為綠燈期間駛離展寬段的車輛數(shù),也是后續(xù)駛?cè)胝箤挾蔚能囕v數(shù);T為后續(xù)駛?cè)胝箤挾蝝輛車所用的時(shí)間[13]。求解單周期內(nèi)單車道通過車輛的區(qū)段延誤問題可以轉(zhuǎn)化為求解不規(guī)則圖形ABCDEFG的面積問題,即為:

        (7)

        式中,D為單車道區(qū)間延誤;Δ1為圖形ABJ的面積;Δ2為圖形IEF的面積。

        2.2展寬段駛?cè)肓髀蕅s和消散飽和流率S恒定

        過飽和引起排隊(duì)溢出的情況,單周期展寬段內(nèi)車輛無法完全釋放,釋放車流流率經(jīng)歷了綠燈初期流率快速增長階段和以飽和流率消散兩個(gè)階段。引入當(dāng)量飽和流率Se的概念[14],將飽和流率轉(zhuǎn)化為定值,計(jì)算公式為:

        (8)

        式中,Se為車流當(dāng)量飽和流率;S(t)為綠燈初期車流流率隨時(shí)間的變化函數(shù);S1為車流飽和流率;g為相位有效綠燈時(shí)間;g0為綠燈初期流率增長階段時(shí)長。

        同樣將展寬段車流駛?cè)肓髀兽D(zhuǎn)化為定值,計(jì)算公式為:

        (9)

        式中,qe為當(dāng)量車流駛?cè)肓髀?;qs(t)為車流駛?cè)肓髀孰S時(shí)間的變化函數(shù)。

        從而繪制出展寬段駛?cè)肓髀屎拖柡土髀屎愣〞r(shí)的排隊(duì)車輛數(shù)Qu與延誤時(shí)間t的關(guān)系曲線,如圖3所示。

        圖3 qs與S恒定時(shí)的排隊(duì)車輛數(shù)與延誤時(shí)間Fig.3 Number of queuing vehicles and delay time when qs and S are steady

        從而得出各參數(shù)的計(jì)算公式如下:

        (10)

        (11)

        (12)

        m=Seg,

        (13)

        T=Seg/qe,

        (14)

        這種情況下,單周期內(nèi)單車道延誤即可簡化為:

        D=n(r+g)+1/2mg-1/2mT-mδ。

        (15)

        把式(10)~(14)代入式(15)得:

        (16)

        整個(gè)展寬段區(qū)段延誤Dall的計(jì)算公式即為:

        (17)

        在一個(gè)信號周期內(nèi),展寬段內(nèi)共進(jìn)入與通過的車輛數(shù)為3(n+m),所以車輛平均延誤為:

        (18)

        3本文延誤模型與定數(shù)理論的比較

        3大理論中,定數(shù)理論與本文研究模型的假設(shè)條件較為接近,為驗(yàn)證本文研究模型的優(yōu)越性,以圖1所示交叉口為例,分別依照定數(shù)理論和本文提出的區(qū)段延誤模型繪制單車道排隊(duì)車輛數(shù)Qu與延誤時(shí)間t隨周期的變化趨勢,如圖4所示。

        圖4 定數(shù)理論與本文延誤模型的比較Fig.4 Comparison of fixed number theory and proposed delay model

        假設(shè)展寬段后路段車流到達(dá)率恒定,從圖4可以看出,定數(shù)理論單車道延誤即為不規(guī)則圖形ABCHG的面積,定數(shù)理論每周期排隊(duì)車輛數(shù)和車輛總延誤隨著周期的延續(xù)逐漸增加,而本文提出的區(qū)段延誤模型所得的排隊(duì)車輛數(shù)與車輛總延誤在定值周圍浮動(dòng)。假設(shè)展寬段內(nèi)單車道排隊(duì)車輛數(shù)到達(dá)n輛排隊(duì)溢出時(shí)為研究的第1周期,則利用定數(shù)理論計(jì)算的該周期總延誤為:

        (19)

        式中,n1=3n,ql為展寬段后路段車流到達(dá)率。

        利用本文延誤模型計(jì)算的總延誤為:

        (20)

        前者比后者多出的總延誤為:

        (21)

        利用定數(shù)理論計(jì)算的第2周期總延誤為:

        (22)

        式中n2=n1+qlC-Seg。

        利用本文區(qū)段延誤模型計(jì)算的總延誤為:

        (23)

        第2周期前者比后者多出的總延誤為:

        (24)

        溢出車流消散之前的第i周期,定數(shù)理論比本文延誤模型多出的總延誤為:

        (25)

        總延誤增加量的斜率為:

        (26)

        由以上計(jì)算過程可以得出,過飽和狀態(tài)下定數(shù)理論計(jì)算出的總延誤量與周期呈線性正相關(guān),延誤量增長速率與每周期內(nèi)滯留的車輛數(shù)及信號周期有關(guān),用定數(shù)理論延誤模型計(jì)算出的總延誤大于實(shí)際展寬段每周期通過車輛所產(chǎn)生的延誤。相反,本文提出的展寬段區(qū)段延誤模型能夠很好地反映展寬段在排隊(duì)溢出狀態(tài)下進(jìn)口道總延誤和排隊(duì)長度的特性。

        4仿真驗(yàn)證

        以圖1所示交叉口為研究對象,應(yīng)用VISSIM微觀仿真軟件進(jìn)行區(qū)段延誤模型的模擬驗(yàn)證,仿真界面如圖5所示。設(shè)該信號交叉口采用單口放行,各相位綠燈時(shí)長25 s,黃閃時(shí)長3 s,周期總計(jì)時(shí)長112 s,西進(jìn)口路段車流到達(dá)率取1 200 veh/h,展寬段單車道車流駛?cè)肓髀嗜? 500 veh/h,單車道車流消散飽和流率取1 600 veh/h,展寬漸變段長度為20 m,展寬段長度為60 m,車頭間距取7 m,啟動(dòng)波波速取5 m/s,車流運(yùn)行數(shù)據(jù)自第3周期開始檢測。分別應(yīng)用定數(shù)理論模型和本文區(qū)段延誤模型進(jìn)行單周期通過車輛總延誤計(jì)算,并將計(jì)算結(jié)果與仿真所得區(qū)段總延誤通過圖表形式展現(xiàn)出來,具體如表1和圖6所示。

        圖5 VISSIM仿真界面Fig.5 VISSIM simulation interface

        第3周期第4周期第5周期第6周期第7周期仿真延誤值/s30243056304730633054優(yōu)化模型計(jì)算值/s31353135313531353135定數(shù)理論計(jì)算值/s53705818626667147162

        圖6 三種所得延誤值對比圖(單位:s)Fig.6 Contrast figure of 3 kinds of delay value(unit:s)

        由表1和圖6可以看出,仿真延誤值以3 049 s為中心上下輕微波動(dòng)??傮w來說,仿真延誤值略小于區(qū)段延誤模型計(jì)算值,較為穩(wěn)定且不隨周期的延續(xù)發(fā)生較大的波動(dòng)。區(qū)段延誤模型計(jì)算值為定值3 135 s,不隨周期的延續(xù)而發(fā)生變化。定數(shù)理論計(jì)算值在第3周期起始已經(jīng)大于仿真延誤值和區(qū)段延誤模型計(jì)算值,且隨著周期的延續(xù),定數(shù)理論計(jì)算值將以448 s為公差逐漸增加,直至無窮大。因此可以看出,與定數(shù)理論計(jì)算方法比較,本文提出的區(qū)段延誤模型能夠更好地反映信號交叉口單周期通過展寬段車輛的區(qū)段延誤。

        5結(jié)論

        本文在定數(shù)理論的基礎(chǔ)上,考慮展寬段長度和飽和流率的波動(dòng)性,建立了展寬段排隊(duì)溢出狀態(tài)下的區(qū)段延誤模型。該模型能較好地反映溢出狀態(tài)下單周期通過展寬段車輛所產(chǎn)生的延誤,為今后展寬段的科學(xué)設(shè)計(jì)提出參考依據(jù)。但模型的建立與研究過程中仍有一些不足之處,如未考慮右轉(zhuǎn)車輛對直行和左轉(zhuǎn)車輛的干擾以及右轉(zhuǎn)車輛產(chǎn)生的延誤。對常規(guī)四相位設(shè)置情況下區(qū)段延誤模型的建立、區(qū)段延誤模型與展寬段設(shè)置形式之間的內(nèi)在邏輯關(guān)系均有待進(jìn)一步的研究。

        參考文獻(xiàn):

        References:

        [1]GB 50647—2011, 城市道路交叉口規(guī)劃規(guī)范[S].

        GB 50647—2011, Code for Planning of Intersections on Urban Roads[S].

        [2]LEISCH J E. Capacity Analysis Techniques for Design of Signalized Intersections: Installment No. 2 Part IV: High Type Facilities and Interchanges[J]. Public Roads, 1967, 34(10): 171-210.

        [3]LEE J J. Impact of Lane Drops on Intersection Capacity[D]. North Caroline: North Carolina State University, 2005.

        [4]TIAN Z Z, WU N. Probabilistic Model for Signalized Intersection Capacity with a Short Right-turn Lane[J]. Journal of Transportation Engineering, 2006, 132(3): 205-212.

        [5]楊曉光, 趙靖, 曾瀅, 等. 短車道對信號交叉口通行能力影響研究[J]. 公路交通科技, 2008, 25(12): 151-156.

        YANG Xiao-guang, ZHAO Jing, ZENG Ying, et al. Research on Impact of Short Lane on Signalized Intersection Capacity[J]. Journal of Highway and Transportation Research and Development, 2008, 25(12): 151-156.

        [6]程國柱. 路口拓寬條件下信號交叉口通行能力研究[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2002.

        CHENG Guo-zhu. Research on Capacity of Widened Signalized Intersection[D]. Harbin: Harbin Institute of Technology, 2002.

        [7]李江. 交通工程學(xué)[M]. 北京: 人民交通出版社, 2002: 66-67.

        LI Jiang. Traffic Engineering[M]. Beijing: China Communications Press, 2002: 66-67.

        [8]WEBSTER F V. Traffic Signal Settings, Road Research Technical Paper No. 39[R]. London: Road Research Laboratory, 1958.

        [9]MILLER A J. Settings for Fixed-cycle Traffic Signals[J]. Journal of the Operational Research Society, 1963, 14(4): 373-386.

        [10]AKCELIK R. Comments on the Application of Queueing Theory to Delays at Signals[J]. Australian Road Research, 1990, 20(3): 53-61.

        [11]MAY A D, HELLER H E M. A Deterministic Queueing Model[J]. Transportation Research, 1967, 1(2): 117-128.

        [12]KIMBER R M, HOLLIS E M. Peak Period Traffic Delays at Road Junctions and other Bottlenecks[J]. Traffic Engineering and Control, 1978, 19(10): 442-446.

        [13]劉巖, 王殿海. 基于短連線的過飽和信號交叉口最大延誤模型[J].中國公路學(xué)報(bào),2011,24(6):91-95.LIU Yan, WANG Dian-hai. Maximum Delay Model in Short Length Signalized Intersection under Supersaturated Traffic Flow[J]. China Journal of Highway and Transport, 2011, 24(6): 91-95.

        [14]尚斌. 基于信號交叉口渠化段長度的配時(shí)參數(shù)優(yōu)化方法[D]. 長春: 吉林大學(xué), 2007.

        SHANG Bin. Optimal Method of Signal Parameters Based on Channelization Length at Signalized Intersections[D]. Changchun: Jilin University, 2007.

        Study on Section Delay Model under Widening Zone Queue Overflow Condition

        CAI Yan-feng1, SUN Da-zhi1, TANG Gui-tao2

        (1. School of Transportation and Logistics, Southwest Jiaotong University, Chengdu Sichuan 610031, China;2. School of Traffic Engineering, Kunming University of Science and Technology,Kunming Yunnan 650093, China)

        Abstract:In order to study the section delay model for vehicles passing widening zone in single cycle at oversaturated signalized intersection under queue overflow condition, first, we introduced 3 classical delay models, and put forward none of them considered the widening zone length and the fluctuation of saturation flow rate. Second, taking the widening zone as the research object, we derived the widening zone delay model on the basis of the classical delay model theory. Finally, we compared the optimization model and the fixed number theory based delay model by formula deriving/calculating and VISSIM simulation. The result shows that the proposed model can effectively reflect the section delay in widening zone under queue overflow condition, it can provide a theoretical basis for scientific design of widening zone.

        Key words:traffic engineering; widening zone; VISSIM simulation; section delay model; queue overflow

        收稿日期:2015-10-08

        作者簡介:蔡艷鋒(1986-),男,湖南益陽人,博士研究生. (cyfsecret@126.com)

        doi:10.3969/j.issn.1002-0268.2016.06.020

        中圖分類號:U491.2+3

        文獻(xiàn)標(biāo)識(shí)碼:A

        文章編號:1002-0268(2016)06-0123-05

        怡春院欧美一区二区三区免费| 成年人视频在线观看麻豆| 久久精品亚州中文字幕| 成人欧美一区二区三区在线观看 | 亚洲精品国产二区三区在线| 99热婷婷一区二区三区| 欧美怡春院一区二区三区| 三级特黄60分钟在线观看| 精品少妇ay一区二区三区| 亚洲国产毛片| 国产精品亚洲av国产| 日本一区二区三区熟女俱乐部 | 亚洲岛国一区二区三区| 香港aa三级久久三级| 黄又色又污又爽又高潮动态图| 97色噜噜| 在线观看精品国产福利片87| 米奇亚洲国产精品思久久| 国产自拍精品在线免费观看| 国产精品理论片| 在线亚洲午夜理论av大片| 视频一区二区在线播放| 国产一区二区精品久久凹凸| 91热久久免费频精品99| 国产精品国产三级久久| 狠狠躁天天躁中文字幕| 国产精选免在线观看| 国产三级三级精品久久| 久久综合伊人77777麻豆| 久久久久久九九99精品| 亚洲五月天综合| 日韩精品电影在线观看| 白白白色视频在线观看播放| 好吊妞视频这里有精品| 亚洲av无码精品色午夜果冻不卡 | 免费a级毛片又大又粗又黑| 亚洲精品夜夜夜妓女网| 亚洲欧美日韩综合久久| 中文字幕亚洲综合久久| 亚洲美女毛片在线视频| 无码吃奶揉捏奶头高潮视频|