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

    交通發(fā)生吸引源集聚影響下路網(wǎng)宏觀基本特性*

    2020-04-30 08:33:28丁恒周靜文鄭小燕柏海艦張衛(wèi)華
    物理學(xué)報(bào) 2020年7期
    關(guān)鍵詞:影響實(shí)驗(yàn)

    丁恒 周靜文 鄭小燕 柏海艦 張衛(wèi)華

    (合肥工業(yè)大學(xué)汽車與交通工程學(xué)院, 合肥 230009)

    宏觀基本圖不受路網(wǎng)外部交通發(fā)生吸引源的影響, 描述了勻質(zhì)性路網(wǎng)內(nèi)累積車輛數(shù)與路網(wǎng)旅行完成率之間的關(guān)系, 可直觀表達(dá)路網(wǎng)宏觀基本特性.然而, 當(dāng)路網(wǎng)內(nèi)交通發(fā)生吸引源發(fā)生變化時(shí), 會(huì)影響路網(wǎng)交通密度, 從而對(duì)宏觀基本特性產(chǎn)生影響.為了分析交通發(fā)生吸引源不同集聚狀態(tài)對(duì)宏觀基本圖的影響規(guī)律, 以交通發(fā)生吸引源發(fā)生吸引量和路段阻抗為參數(shù), 建立交通發(fā)生吸引源聚集度模型.以方格式路網(wǎng)為分析對(duì)象,設(shè)計(jì)9組不同交通發(fā)生吸引源方案, 通過對(duì)比交通發(fā)生吸引源聚集度與路網(wǎng)宏觀基本圖曲線, 發(fā)現(xiàn)只有當(dāng)路網(wǎng)處于擁擠流狀態(tài)時(shí), 交通發(fā)生吸引源聚集現(xiàn)象才會(huì)對(duì)宏觀基本圖產(chǎn)生影響, 并且路網(wǎng)中交通發(fā)生吸引源聚集度越低(交通發(fā)生吸引量分布越均衡), 路網(wǎng)旅行完成率就會(huì)越高.

    1 引 言

    宏觀基本圖(macroscopic fundamental diagram, MFD)描述了均質(zhì)性區(qū)域路網(wǎng)內(nèi)車輛密度 (veh·km–1)與空間平均流量 (veh·h–1)之間單峰低散射關(guān)系, 即該區(qū)域內(nèi)累積車輛數(shù) n (t) 與路網(wǎng)旅行完成率 G (n(t)) 之間的固有物理特性, 如圖1(a)所示.在該特性中, 路網(wǎng)旅行完成率 G (n(t)) 隨著累積車輛數(shù) n (t) 變化而變化.當(dāng) n (t)ncr時(shí), 路網(wǎng)旅行完成率 G (n(t)) 又隨著累積車輛數(shù) n (t) 增大而減小, 其中 ncr是臨界累積車輛數(shù).Gerolimins和 Daganzo[1]及 Gonzales等[2]根據(jù)路網(wǎng)宏觀交通流的運(yùn)行狀態(tài), 將MFD劃分為自由流、臨界流和擁擠流3種狀態(tài), 如圖1(b)的Ⅰ、Ⅱ和Ⅲ.由于MFD可直觀地描述區(qū)域路網(wǎng)的宏觀交通狀態(tài), 簡(jiǎn)化了采用復(fù)雜OD數(shù)據(jù)對(duì)大規(guī)模交通網(wǎng)絡(luò)建模過程, 被廣泛應(yīng)用于分析大范圍交通網(wǎng)絡(luò).例如 Haddad、Yang、Ding和 Kim 等根據(jù)MFD描述的路網(wǎng)狀態(tài), 提出了單個(gè)區(qū)域[3?4]、兩個(gè)區(qū)域[5?6]和多個(gè)區(qū)域[7?9]的邊界限流控制方法, 有效地解決了大范圍交通網(wǎng)絡(luò)擁堵區(qū)交通控制問題.

    MFD的概念由Godfrey[10]在1969年首次提出, 但限于當(dāng)時(shí)的交通流檢測(cè)手段, 其存在性并未得到驗(yàn)證.近些年隨著數(shù)據(jù)采集及仿真手段的進(jìn)步, 2008年Geroliminis和 Daganzo[1]根據(jù)日本橫濱路網(wǎng)數(shù)據(jù)、Gonzales等[2]根據(jù)肯尼亞內(nèi)羅畢路網(wǎng)數(shù)據(jù)驗(yàn)證了MFD的存在性, 此外Gao[11]運(yùn)用斯德哥爾摩的交通數(shù)據(jù)發(fā)現(xiàn)了該城市路網(wǎng)也存在MFD.除了上述通過數(shù)據(jù)分析研究MFD存在外,也有部分學(xué)者通過解析模型分析MFD.如Leclercq等[12]利用變分法對(duì)不同路徑上的MFD進(jìn)行估計(jì),然后將其聚合成一個(gè)統(tǒng)一的MFD, 進(jìn)而確定與系統(tǒng)最優(yōu)值相對(duì)應(yīng)的MFD, 該方法對(duì)許多“非同質(zhì)性”網(wǎng)絡(luò)也適用; Daganzo等[13]利用舊金山路網(wǎng)數(shù)據(jù)提出了MFD的解析近似表達(dá)式; Courbon等[14]對(duì)比分析了獲取MFD的3種方法: 分析法、基于路徑的方法、檢測(cè)器法, 并對(duì)該三種方法進(jìn)行了仿真評(píng)價(jià).

    圖1 MFD (a) MFD 的基本特征; (b) 路網(wǎng)狀態(tài)劃分Fig.1.MFD: (a) Basic characteristics of MFD; (b) state classification of MFD curve.

    雖然MFD是特定路網(wǎng)的固有屬性, 但其分布仍然受到路網(wǎng)一些特征參數(shù)的影響.相關(guān)學(xué)者對(duì)影響MFD因素的研究主要從交通條件、道路條件、管控條件以及選擇行為4個(gè)方面進(jìn)行.首先, 關(guān)于交通條件對(duì)MFD的影響[15?17], 典型的有Gayah和Daganzo[18]提出交通擁堵分布不均勻會(huì)導(dǎo)致MFD的“滯回現(xiàn)象”; 許菲菲等[19]采用不同的交通管控措施, 得出交通需求的劇烈變化、公交專用道的設(shè)置、車道禁行都會(huì)不同程度的影響MFD; 朱琳等[20]通過對(duì)比分析路網(wǎng)平均流量-平均密度關(guān)系的時(shí)段分布差異, 證明了交通密度分布的不均衡性是影響路網(wǎng)宏觀交通狀態(tài)的根本因素; Geroliminis和Sun[21]提出路網(wǎng)中車輛密度的空間分布是影響MFD散射特性及其形狀的關(guān)鍵因素之一; Geroliminis等[22]發(fā)現(xiàn)交叉口處左轉(zhuǎn)交通流的存在降低了MFD的最大值.其次道路條件對(duì)MFD的影響研究, 如Buisson等[23]分析了城市路網(wǎng)、穿越城市的高速路網(wǎng)與環(huán)城高速路網(wǎng)對(duì)MFD的影響, 并認(rèn)為高速網(wǎng)絡(luò)不存在MFD, 并且外部OD需求對(duì)MFD是沒有影響的.再次, 不同的交通管控條件也會(huì)影響MFD, 其中管控措施和參數(shù)的研究受到重點(diǎn)關(guān)注.Jin等[24]通過元胞傳輸模型發(fā)現(xiàn)在自由流狀態(tài)下, MFD比較穩(wěn)定, 而間斷流狀態(tài)下信號(hào)控制使得網(wǎng)絡(luò)交通流不穩(wěn)定, 進(jìn)而影響 MFD形狀;Alonso等[25]基于真實(shí)數(shù)據(jù)的統(tǒng)計(jì)分析, 得出交通管制影響鏈路和路網(wǎng)級(jí)別的交通流變化.最后, 針對(duì)出行者不同的選擇行為, Leclercq和Geroliminis[26]通過研究局部非均質(zhì)性交通流對(duì)MFD形狀的影響, 可知路徑選擇和異構(gòu)的局部交通容量可能會(huì)導(dǎo)致交通擁堵的時(shí)空分布不一致, 進(jìn)而影響路網(wǎng)MFD的形狀和分布; Mazloumian等[27]基于對(duì)MFD空間分布不均勻的研究, 可知路徑選擇和交通需求是影響MFD散射的重要因素.總體而言,不同的影響因素歸根結(jié)底都是影響路網(wǎng)密度, 進(jìn)而對(duì)MFD產(chǎn)生影響.在交通發(fā)生吸引源分析方面,雖然Buisson等[23]證明了外部OD需求對(duì)區(qū)域路網(wǎng)的MFD是沒有影響的, 但并沒有說明路網(wǎng)內(nèi)部的交通發(fā)生吸引源分布是否對(duì)MFD產(chǎn)生影響.而實(shí)際上, 路網(wǎng)內(nèi)交通發(fā)生吸引源分布不均勻時(shí), 其產(chǎn)生的交通量會(huì)改變局部路網(wǎng)的交通密度, 進(jìn)而可能影響MFD分布.

    為了分析交通發(fā)生吸引源對(duì)路網(wǎng)MFD產(chǎn)生的影響規(guī)律, 本文主要開展以下兩個(gè)方面工作:1)根據(jù)區(qū)域路網(wǎng)內(nèi)交通發(fā)生吸引源的時(shí)空分布,以交通發(fā)生吸引源發(fā)生吸引量和路段阻抗為動(dòng)態(tài)參數(shù), 提出交通發(fā)生吸引源聚集度分析模型; 2)以4 × 4方格網(wǎng)狀的區(qū)域路網(wǎng)為研究對(duì)象, 根據(jù)不同發(fā)生吸引源的分布, 通過數(shù)據(jù)驗(yàn)證不同聚集度條件下交通發(fā)生吸引源對(duì)MFD的影響規(guī)律.

    2 交通發(fā)生吸引源聚集度模型

    聚集度的概念在經(jīng)濟(jì)學(xué)和地理學(xué)等領(lǐng)域中應(yīng)用比較廣泛, 一般用來形容某個(gè)區(qū)域內(nèi)某個(gè)經(jīng)濟(jì)產(chǎn)業(yè)或某個(gè)地理指標(biāo)的聚集程度.與產(chǎn)業(yè)在空間和發(fā)生吸引強(qiáng)度上分布類似, 路網(wǎng)中交通發(fā)生吸引源也存在著一定的聚集形態(tài).因此將聚集度概念引入交通領(lǐng)域, 用來描述交通發(fā)生吸引源的集聚狀態(tài).然而, 與產(chǎn)業(yè)分布不同在于, 路網(wǎng)中各個(gè)路段的阻抗是實(shí)時(shí)變化的, 因此交通發(fā)生吸引源聚集程度也變化的較為頻繁.因此, 在分析交通發(fā)生吸引源聚集形態(tài)時(shí), 必須考慮到時(shí)間因素.

    假設(shè)路網(wǎng)內(nèi)有N個(gè)交通發(fā)生吸引源, 路網(wǎng)可劃分為M個(gè)區(qū)域.根據(jù)發(fā)生吸引源數(shù)量、發(fā)生吸引量規(guī)模、路網(wǎng)容量, 建立交通發(fā)生吸引源聚集度評(píng)估模型, 步驟如下:

    首先, 根據(jù)N個(gè)交通發(fā)生吸引源空間位置坐標(biāo), 通過M個(gè)區(qū)域劃分對(duì)交通發(fā)生吸引源進(jìn)行分類.分類方法是根據(jù)路段密度和路段速度, 通過路網(wǎng)交通狀態(tài)識(shí)別[28]獲得交通發(fā)生吸引源影響下的路段交通狀態(tài), 分為非常暢通、暢通、輕度擁堵、中度擁堵以及嚴(yán)重?fù)矶?種狀態(tài)等級(jí), 通過最小生成樹法初始劃分、歸一化割法深度劃分以及反向動(dòng)態(tài)合并的過程得到對(duì)交通發(fā)生吸引源分類的M個(gè)區(qū)域[29].

    其次, 計(jì)算交通發(fā)生吸引源空間分布基尼系數(shù):

    式中, Pi表示劃分的M個(gè)區(qū)域中第i個(gè)區(qū)域交通發(fā)生吸引源容量(停車需求)占整個(gè)路網(wǎng)所有發(fā)生吸引源容量(停車需求)的比重, xi表示第i個(gè)區(qū)域路網(wǎng)交通容量占整個(gè)路網(wǎng)交通容量的比重.基尼系數(shù) G =0 時(shí), 表明停車需求在空間分布是均勻的,G越大(最大值為1), 表明停車需求在空間分布越不均勻.

    再次, 采用 Ellision和 Glaeser[30]提出的 EG集聚指數(shù)來測(cè)定交通發(fā)生吸引源聚集度指數(shù), 即評(píng)估模型為

    式中: γ 為交通發(fā)生吸引源聚集度指數(shù); H為赫芬達(dá)爾指數(shù),表示 N 個(gè)交通吸引源中第j個(gè)區(qū)域交通發(fā)生吸引源總?cè)萘?停車需求)占整個(gè)路網(wǎng)所有發(fā)生吸引源容量(停車需求)的比重.

    考慮到交通流的動(dòng)態(tài)變化, (2)式中的空間分布基尼系數(shù)G以及赫芬達(dá)爾指數(shù)H也會(huì)隨著路網(wǎng)交通狀態(tài)即路網(wǎng)交通阻抗的變化而變化.交通越擁堵的地方, 交通發(fā)生吸引強(qiáng)度越大.為了更好地描述路網(wǎng)聚集狀態(tài), 假設(shè)交通發(fā)生吸引源的聚集狀態(tài)與道路交通阻抗呈正相關(guān), 即交通發(fā)生吸引源配置越聚集, 道路交通阻抗越大, 且道路交通阻抗在聚集度模型中屬于正向因子, 則考慮交通阻抗的集聚度評(píng)估改進(jìn)模型為

    考慮路網(wǎng)實(shí)時(shí)狀態(tài), (3)式中引入了路網(wǎng)交通阻抗 ta, 其計(jì)算需要考慮路段不同交通參數(shù).根據(jù)四兵鋒等[31]提出的道路阻抗模型, 建立考慮多種路面的道路阻抗函數(shù)模型:

    3 聚集度指數(shù)計(jì)算過程

    為滿足實(shí)際交通管控的需要, 在交通發(fā)生吸引源聚集度模型基礎(chǔ)上, 計(jì)算路網(wǎng)交通發(fā)生吸引源聚集指數(shù), 并以此分析路網(wǎng)交通發(fā)生吸引源的動(dòng)態(tài)集聚狀態(tài).步驟如下:

    Step1選擇一個(gè)路網(wǎng), 確定路網(wǎng)內(nèi)交通發(fā)生吸引源的分布、規(guī)模.

    Step2劃分子區(qū).在保證空間連續(xù)和子區(qū)內(nèi)密度異質(zhì)性最小化的前提下, 根據(jù)交通流數(shù)據(jù), 然后考慮區(qū)域路網(wǎng)內(nèi)的交通發(fā)生吸引源、交通擁堵狀態(tài)以及交通流大小等因素將分布不均勻的區(qū)域路網(wǎng)劃分為幾個(gè)密度均勻的子區(qū).

    Step3計(jì)算路段交通阻抗.路段阻抗與斷面形式有關(guān).路段的阻抗函數(shù)的計(jì)算步驟基本類似,由于單幅路面的阻抗函數(shù)最復(fù)雜, 本文以單幅路為例給出段阻抗函數(shù)計(jì)算步驟如下:

    首先, 由(4)式道路交通阻抗模型的定義, 需要確定路段零流量條件下的阻抗.根據(jù)1976年Branston[32]路段阻抗函數(shù)的理論特性: 1)當(dāng)流量充分小時(shí), 路段阻抗接近于“零流”阻抗; 2)在流量遠(yuǎn)小于道路通行能力時(shí), 路段阻抗隨流量變化而緩慢變化; 3)在“穩(wěn)態(tài)”系統(tǒng)狀態(tài)下, 阻抗函數(shù)曲線變成飽和流量縱坐標(biāo)的漸近線.可以通過交通調(diào)查或者仿真的手段獲得不同路段上一般機(jī)動(dòng)車的平均零流阻抗.

    然后, 標(biāo)定參數(shù).對(duì)于單幅路面的路段, 假定共觀測(cè)到Z組樣本值, 其中第k組樣本值為以及 Ca, Ca,是確定的.為便于描述作如下變換:

    確定α1,α2,α3,α4,β1,β2,β3,β4參數(shù), 滿足:

    根據(jù)一階極值條件可得關(guān)于參數(shù)a1, a2, a3,a4, b1, b2, b3, b4的非線性方程組.按照常規(guī)的方法求解非線性方程組的解析解難度較大, 本文采用Broyden求解方法[33].為了降低求解誤差的影響, 提高參數(shù)數(shù)值解的準(zhǔn)確性, 設(shè)置求解精度為10–8, 可得滿足精度要求的參數(shù)數(shù)值解.由于該非線性方程組存在多組解, 所以在計(jì)算迭代時(shí), 對(duì)初始迭代值的選取有一定要求.因此,

    最后, 獲得每條路段的交通阻抗, 進(jìn)而得到路網(wǎng)整體交通阻抗.

    Step4由Step2獲得各子區(qū)的交通發(fā)生吸引源的相關(guān)信息, 然后根據(jù)(1)式計(jì)算得到空間分布基尼系數(shù)G, 并計(jì)算赫芬達(dá)爾指數(shù)H, 將Step3得到的路網(wǎng)交通阻抗帶入(3)式, 進(jìn)而可以獲得實(shí)時(shí)的路網(wǎng)交通發(fā)生吸引源的聚集度指數(shù).

    Step5結(jié)束.

    4 交通發(fā)生吸引源聚集效應(yīng)對(duì)路網(wǎng)MFD的影響分析

    根據(jù)Gerolimins等[1]的研究結(jié)果, 存在理想MFD的子區(qū)范圍一般在4—8 km2.在該區(qū)域范圍內(nèi), 路網(wǎng)結(jié)構(gòu)較為單一, 一般接近于方格網(wǎng)狀路網(wǎng).因此, 為了分析交通發(fā)生吸引源聚集度指數(shù)變化對(duì)路網(wǎng)宏觀基本特性的影響, 本文采用4 × 4方格形路網(wǎng)進(jìn)行交通仿真, 如圖2.該路網(wǎng)有16個(gè)交叉口, 96 條路段, 且每條路段均為三車道.將該路網(wǎng)外圍的16個(gè)路口作為路網(wǎng)背景交通量輸入點(diǎn), 然后在路網(wǎng)內(nèi)部設(shè)置8個(gè)車流輸入點(diǎn)作為交通發(fā)生吸引源, 每 300 s 采集交通流數(shù)據(jù)一次, 并得到路網(wǎng)MFD數(shù)據(jù)以及相應(yīng)的聚集度指數(shù).

    圖2 仿真路網(wǎng)Fig.2.Simulation road network.

    4.1 交通發(fā)生吸引源聚集效應(yīng)對(duì)路網(wǎng)密度影響分析

    路網(wǎng)內(nèi)交通發(fā)生吸引源的聚集狀態(tài)一定程度上可以由路網(wǎng)中交通密度分布狀態(tài)體現(xiàn).在上述路網(wǎng)中設(shè)置不同配置的交通發(fā)生吸引源, 如表1所示, 在相同的背景交通流的影響下獲得某一個(gè)信號(hào)周期(150 s)內(nèi)的路網(wǎng)交通密度分布, 如圖3所示.

    圖3中 (a1)、(a2)、(a3)、(a4)和 (a5)是相同的, 均表示在未加載交通發(fā)生吸引源情況下路網(wǎng)交通密度分布.圖3中 (b1)、(b2)、(b3)、(b4)和(b5)分別表示5種總規(guī)模相同但分布不同的交通發(fā)生吸引源配置情況.圖3 中 (c1)、(c2)、(c3)、(c4)和(c5)分別表示在路網(wǎng)背景交通密度和不同配置的交通發(fā)生吸引源條件聯(lián)合影響下, 某一個(gè)信號(hào)周期(150 s)內(nèi)發(fā)生吸引交通量影響下路網(wǎng)交通密度分布.可以發(fā)現(xiàn), 在未加載交通發(fā)生吸引源的情況下, 路網(wǎng)交通密度分布是比較均勻的, 最大的單車道交通密度大約是 5 pcu·km–1; 當(dāng)路網(wǎng)中加載交通發(fā)生吸引源之后, 路網(wǎng)交通密度分布明顯發(fā)生變化, 并且局部最大單車道交通密度達(dá)到134.7 pcu·km–1, 遠(yuǎn)超過 5 pcu·km–1, 并且交通密度較大的路段都是與交通發(fā)生吸引源相連或鄰近的路段.因此, 交通發(fā)生吸引源的聚集效應(yīng)會(huì)改變路網(wǎng)交通密度的分布.

    表1 交通發(fā)生吸引源配置參數(shù)Table 1.Traffic generation and attraction source configuration parameters.

    4.2 交通發(fā)生吸引源聚集效應(yīng)對(duì)路網(wǎng)MFD動(dòng)態(tài)影響分析

    4.2.1 參數(shù)標(biāo)定

    由4.1節(jié)可知, 交通發(fā)生吸引源的聚集效應(yīng)會(huì)影響路網(wǎng)交通密度的分布, 所以進(jìn)而會(huì)影響路網(wǎng)MFD的基本特性.根據(jù)(3)式, 在不同的交通發(fā)生吸引源分布條件下, 分析路網(wǎng)內(nèi)交通發(fā)生吸引源的聚集狀態(tài)對(duì)路網(wǎng)MFD的動(dòng)態(tài)影響.將路網(wǎng)內(nèi)的8個(gè)交通發(fā)生吸引源的總規(guī)模分別設(shè)置為800,1200, 2000 pcu·h–13 種情況, 獲得路網(wǎng)累積車輛數(shù)(n)與路網(wǎng)旅行完成率(G)之間的關(guān)系, 即路網(wǎng)MFD, 如圖4所示.

    根據(jù) Gonzales和 Chavis[34]的研究發(fā)現(xiàn), 盡管MFD的形狀依賴于OD需求, 然而此關(guān)系并不隨時(shí)間變化, 在一定交通范圍內(nèi)可以采用如G(n(t))=a ·(n(t))3+b·(n(t))2+c·n(t)+d的3次方程式近似表示.由實(shí)驗(yàn)數(shù)據(jù)可得, 路網(wǎng)MFD擬合曲線參數(shù)標(biāo)定如表2所示.

    根據(jù)不同交通發(fā)生吸引源條件下MFD數(shù)據(jù)在Ⅰ、Ⅱ和Ⅲ 3種狀態(tài)分布可知, 在Ⅰ、Ⅱ狀態(tài),MFD 并沒有明顯的波動(dòng); 在狀態(tài)Ⅲ, 800 pcu·h–1和 1200 pcu·h–1兩組數(shù)據(jù)的 MFD 有明顯的較大波動(dòng), 而 2000 pcu·h–1組仿真得到的 MFD 波動(dòng)幅度較小.因此, 一定程度上路網(wǎng)宏觀基本圖會(huì)受到交通發(fā)生吸引源的影響.

    圖3 不同交通發(fā)生吸引源集聚影響下交通密度分布Fig.3.Traffic density distribution under different traffic generation and attraction sources.

    圖4 不同交通發(fā)生吸引源規(guī)模下路網(wǎng) MFDs (a) 800 pcu·h–1; (b) 1200 pcu·h–1; (c) 2000 pcu·h–1Fig.4.MFDs under different traffic generation and attraction source scale: (a) 800 pcu·h–1; (b) 1200 pcu·h–1; (c) 2000 pcu·h–1.

    表2 不同交通發(fā)生吸引源規(guī)模下MFD參數(shù)Table 2.MFD parameters under different traffic generation and attraction sourcescale.

    4.2.2 動(dòng)態(tài)分析

    為了進(jìn)一步研究交通發(fā)生吸引源聚集度模型與路網(wǎng)宏觀基本圖波動(dòng)的動(dòng)態(tài)關(guān)系, 通過動(dòng)態(tài)加載交通流量, 并且每 300 s采集交通流數(shù)據(jù)一次, 開展實(shí)驗(yàn)分析.

    在初始狀態(tài), 路網(wǎng)交通流為自由流, 路網(wǎng)交通阻抗可認(rèn)為零流阻抗.再由聚集度模型的定義, 將路網(wǎng)劃分為4個(gè)子區(qū), 每個(gè)子區(qū)包含2個(gè)交通發(fā)生吸引源, 如圖2所示.設(shè)計(jì)的9組仿真實(shí)驗(yàn)按照交通發(fā)生吸引源總規(guī)模分為3大組, 分別為A組、B組和C組, 每大組包含三個(gè)實(shí)驗(yàn)方案, 然后依據(jù)交通發(fā)生吸引源集聚度模型計(jì)算9組仿真實(shí)驗(yàn)的初始聚集度指數(shù), 且每大組的3個(gè)實(shí)驗(yàn)方案的初始聚集度指數(shù)依次遞減, 如表3所示.3大組實(shí)驗(yàn)的MFD及動(dòng)態(tài)聚集度指數(shù)曲線如圖5所示.

    根據(jù)路網(wǎng)的聚集度曲線的變化趨勢(shì)可知, 隨著路網(wǎng)密度的變化, 聚集度曲線有明顯的3個(gè)變化階段, 如圖5所示.參考圖1(b)中MFD的狀態(tài)劃分可知, 聚集度曲線在狀態(tài)Ⅰ時(shí)基本沒有變化, 在狀態(tài)Ⅱ時(shí)開始有小幅的變化, 在狀態(tài)Ⅲ時(shí)開始急劇變化并逐步平穩(wěn)(路網(wǎng)全部阻塞).

    為了更好地對(duì)比分析交通發(fā)生吸引源不同集聚影響下的MFD波動(dòng)特性, 將上述9組仿真實(shí)驗(yàn)MFD數(shù)據(jù)進(jìn)行三次多項(xiàng)式擬合, 相關(guān)參數(shù)標(biāo)定如表4所示.從A組、B組和C組仿真實(shí)驗(yàn)中每一份實(shí)驗(yàn)數(shù)據(jù)的擬合效果來看, R2值都大于0.9, 表明每份數(shù)據(jù)的擬合效果都比較好.

    A組、B組和C組的MFD擬合曲線與相應(yīng)的聚集度曲線如圖5所示.仿真實(shí)驗(yàn)中基本參數(shù)都保持一致, 交通發(fā)生吸引源相同規(guī)模下不同分布是唯一變量, 即交通發(fā)生吸引源不同的聚集度是唯一變量.因此, 路網(wǎng)MFD的變化就是由交通發(fā)生吸引源的聚集效應(yīng)導(dǎo)致的, 而交通發(fā)生吸引源聚集度體現(xiàn)了交通發(fā)生吸引源的聚集效應(yīng), 且仿真過程中每300 s采集一次路網(wǎng)數(shù)據(jù), 可獲得相應(yīng)的MFD數(shù)據(jù)點(diǎn)以及聚集度指數(shù), 則保證了MFD曲線與聚集度曲線基本實(shí)時(shí)對(duì)應(yīng).

    圖5 不同交通發(fā)生吸引源配置條件下MFDs和聚集度曲線 (a1) A組仿真實(shí)驗(yàn)MFD; (a2) A組仿真實(shí)驗(yàn)聚集度曲線;(b1) B 組仿真實(shí)驗(yàn) MFD; (b2) B 組仿真實(shí)驗(yàn)聚集度曲線; (c1) C 組仿真實(shí)驗(yàn) MFD; (c2) C 組仿真實(shí)驗(yàn)聚集度曲線Fig.5.MFDs and aggregation degree curves under different traffic generation and attraction source configuration: (a1) MFDof group A simulation scheme; (a2) aggregation degree curve of group A simulation scheme; (b1) MFDof group B simulation scheme;(b2) aggregation degree curve of group B simulation scheme; (b1) MFDof group C simulation scheme; (b2) aggregation degree curve of group C simulation scheme.

    綜合圖5和表5可知, 在狀態(tài)Ⅰ時(shí), 路網(wǎng)累積車輛數(shù)在0—600 veh之間, A3組平均旅行完成率相比 A1組增加了 0.77%, 相比于 A2組增加了0.61%, A2組平均旅行完成率相比A1組增加了0.16%, A組、B組和C組仿真實(shí)驗(yàn)的MFD擬合曲線都幾乎重疊, 沒有明顯差異, 說明交通發(fā)生吸引源的集聚效應(yīng)對(duì)此狀態(tài)下的MFD幾乎沒有影響; 在狀態(tài)Ⅱ時(shí), 路網(wǎng)累積車輛數(shù)在 600—2000 veh之間, A組、B組和C組仿真實(shí)驗(yàn)的MFD擬合曲線產(chǎn)生一定的差異, A3組平均旅行完成率相比A1組增加了3.13%, 相比于A2組增加了1.14%,A2組平均旅行完成率相比A1組增加了1.96%,說明交通發(fā)生吸引源的集聚效應(yīng)對(duì)此狀態(tài)下的MFD 已經(jīng)產(chǎn)生影響, 但影響有限; 在狀態(tài)Ⅲ時(shí), 路網(wǎng)累積車輛數(shù)在2000—8000 veh之間, A3組平均旅行完成率相比A1組增加了30.96%, 相比于A2組增加了3.29%, A2組平均旅行完成率相比A1組增加了26.79%, 不同交通發(fā)生吸引源條件下MFD擬合曲線之間有明顯的差異, 且A1組MFD擬合曲線偏低、A2組MFD擬合曲線居中、A3組MFD擬合曲線偏高, 對(duì)應(yīng)狀態(tài)Ⅲ的聚集度曲線可知, A1組聚集度曲線值最大、A2組聚集度曲線值相對(duì)次之、A3組聚集度曲線值最小.由表5可知, B組、C組結(jié)果類似, 說明交通發(fā)生吸引源的集聚效應(yīng)對(duì)狀態(tài)Ⅲ下的MFD影響較大, 且交通發(fā)生吸引源聚集度指數(shù)越小, 路網(wǎng)旅行完成率越高, MFD波動(dòng)值越高.

    表3 仿真實(shí)驗(yàn)基本參數(shù)Table 3.Basic parameters of simulation experiment.

    表4 不同交通發(fā)生吸引源配置條件下MFD參數(shù)Table 4.MFD parameters under differenttraffic generation and attraction sourceconfiguration.

    表5 各組仿真實(shí)驗(yàn)平均旅行完成率評(píng)估參數(shù)Table 5.Evaluation parameters of the average trip completionflow of simulation experiment.

    5 結(jié) 論

    本文在MFD波動(dòng)特性的基礎(chǔ)上, 從交通發(fā)生吸引源集聚現(xiàn)象出發(fā), 提出了一種考慮交通阻抗的交通發(fā)生吸引源聚集度模型; 并以9組仿真實(shí)驗(yàn)為例, 獲得路網(wǎng)在不同交通發(fā)生吸引源聚集狀態(tài)下的MFD和聚集度曲線變化特征.對(duì)比不同交通發(fā)生吸引源條件下MFD變化規(guī)律, 有以下結(jié)論:

    1)當(dāng)路網(wǎng)交通流處于自由流狀態(tài)時(shí), 路網(wǎng)交通發(fā)生吸引源集聚不會(huì)對(duì)MFD產(chǎn)生影響; 在臨界流狀態(tài)時(shí), 交通發(fā)生吸引源集聚對(duì)MFD有一定的影響; 當(dāng)路網(wǎng)處于擁擠流狀態(tài)時(shí), 交通發(fā)生吸引源集聚對(duì)MFD有顯著的影響.

    2)在相同的路網(wǎng)流量條件下, 當(dāng)路網(wǎng)處于擁擠流狀態(tài)時(shí), 交通發(fā)生吸引源聚集狀態(tài)越均衡, 即聚集度指數(shù)越小, 路網(wǎng)旅行完成率越高(MFD曲線向上波動(dòng)).

    3)交通發(fā)生吸引源聚集度指數(shù)隨著路網(wǎng)內(nèi)交通發(fā)生吸引源的配置以及路網(wǎng)交通流的變化而變化; 在路網(wǎng)交通發(fā)生吸引源配置不變的情況下, 路網(wǎng)交通發(fā)生吸引源聚集度指數(shù)在自由流狀態(tài)下沒有較大的變化, 在臨界流狀態(tài)下開始有一定的變化, 在擁擠流狀態(tài)下開始急劇下降并逐步平穩(wěn).

    通過本文的研究, 建立了路網(wǎng)交通發(fā)生吸引源與MFD變化之間的宏觀動(dòng)態(tài)變化關(guān)系.依據(jù)該關(guān)系, 可通過調(diào)節(jié)路網(wǎng)內(nèi)各個(gè)停車場(chǎng)交通發(fā)生吸引強(qiáng)度, 改變路網(wǎng)交通發(fā)生吸引源聚集度指數(shù)提高路網(wǎng)旅行完成率, 從而提高路網(wǎng)運(yùn)行效率, 為城市路網(wǎng)內(nèi)大型交通發(fā)生吸引源(停車場(chǎng)、醫(yī)院、學(xué)校、商場(chǎng)等)的配置, 以及交通管控(停車收費(fèi)、路段管制、信號(hào)燈控制等)提供一定的理論支持.由于本文是在一個(gè)規(guī)則的方格形路網(wǎng)分析交通發(fā)生吸引源聚集規(guī)律的, 而實(shí)際路網(wǎng)是比較復(fù)雜的.實(shí)際路網(wǎng)特征的相關(guān)MFD特性還有待進(jìn)一步研究.

    猜你喜歡
    影響實(shí)驗(yàn)
    記一次有趣的實(shí)驗(yàn)
    是什么影響了滑動(dòng)摩擦力的大小
    微型實(shí)驗(yàn)里看“燃燒”
    哪些顧慮影響擔(dān)當(dāng)?
    做個(gè)怪怪長(zhǎng)實(shí)驗(yàn)
    沒錯(cuò),痛經(jīng)有時(shí)也會(huì)影響懷孕
    媽媽寶寶(2017年3期)2017-02-21 01:22:28
    NO與NO2相互轉(zhuǎn)化實(shí)驗(yàn)的改進(jìn)
    實(shí)踐十號(hào)上的19項(xiàng)實(shí)驗(yàn)
    太空探索(2016年5期)2016-07-12 15:17:55
    擴(kuò)鏈劑聯(lián)用對(duì)PETG擴(kuò)鏈反應(yīng)與流變性能的影響
    基于Simulink的跟蹤干擾對(duì)跳頻通信的影響
    午夜福利视频1000在线观看 | 69av精品久久久久久| 国产午夜精品久久久久久| 一进一出好大好爽视频| av视频免费观看在线观看| 亚洲自拍偷在线| 国产av在哪里看| 亚洲av成人不卡在线观看播放网| 一进一出好大好爽视频| 亚洲黑人精品在线| 午夜老司机福利片| 黄色a级毛片大全视频| 又黄又粗又硬又大视频| 亚洲自拍偷在线| 亚洲精品粉嫩美女一区| 高清毛片免费观看视频网站| 欧美日韩亚洲国产一区二区在线观看| 国产免费男女视频| 最好的美女福利视频网| 欧美日本中文国产一区发布| 两性夫妻黄色片| 人妻久久中文字幕网| netflix在线观看网站| 亚洲国产高清在线一区二区三 | 亚洲第一av免费看| 国产成人精品久久二区二区91| 久久这里只有精品19| 免费观看精品视频网站| 欧美成人午夜精品| 18禁国产床啪视频网站| 男人的好看免费观看在线视频 | 久久午夜亚洲精品久久| √禁漫天堂资源中文www| 亚洲在线自拍视频| 国产真人三级小视频在线观看| 99久久久亚洲精品蜜臀av| 亚洲一区二区三区色噜噜| 美女免费视频网站| 他把我摸到了高潮在线观看| 亚洲精品美女久久av网站| 久久香蕉国产精品| 十八禁网站免费在线| 国产精品久久久久久亚洲av鲁大| 999久久久精品免费观看国产| 日韩免费av在线播放| 美女大奶头视频| 在线视频色国产色| 国产精品99久久99久久久不卡| 亚洲精品中文字幕一二三四区| 夜夜爽天天搞| cao死你这个sao货| av超薄肉色丝袜交足视频| 久久香蕉激情| 人人妻人人澡人人看| bbb黄色大片| 午夜亚洲福利在线播放| www.999成人在线观看| 亚洲精品粉嫩美女一区| 免费高清视频大片| 视频区欧美日本亚洲| 少妇裸体淫交视频免费看高清 | 99久久综合精品五月天人人| 亚洲男人的天堂狠狠| 最新美女视频免费是黄的| 黑人欧美特级aaaaaa片| 在线国产一区二区在线| 亚洲一区中文字幕在线| 老司机在亚洲福利影院| 久久国产精品男人的天堂亚洲| 国内毛片毛片毛片毛片毛片| 色综合欧美亚洲国产小说| 老司机深夜福利视频在线观看| 波多野结衣高清无吗| 国产精品九九99| 真人一进一出gif抽搐免费| 成人18禁高潮啪啪吃奶动态图| 国产精品99久久99久久久不卡| 97超级碰碰碰精品色视频在线观看| 亚洲人成电影免费在线| 热re99久久国产66热| 久久精品亚洲精品国产色婷小说| 18禁观看日本| 国产精品自产拍在线观看55亚洲| 可以免费在线观看a视频的电影网站| 韩国av一区二区三区四区| 波多野结衣高清无吗| or卡值多少钱| 欧美色欧美亚洲另类二区 | 亚洲av第一区精品v没综合| 狂野欧美激情性xxxx| 好看av亚洲va欧美ⅴa在| 12—13女人毛片做爰片一| 天堂动漫精品| 国产麻豆69| 黑人操中国人逼视频| 免费看a级黄色片| 叶爱在线成人免费视频播放| 一卡2卡三卡四卡精品乱码亚洲| 激情在线观看视频在线高清| 国产精品久久久av美女十八| 久久久久久久久久久久大奶| 国产精品久久电影中文字幕| 日韩中文字幕欧美一区二区| 国产一区二区在线av高清观看| 久久午夜亚洲精品久久| 日韩大尺度精品在线看网址 | 夜夜夜夜夜久久久久| 日日爽夜夜爽网站| 亚洲精品一区av在线观看| 91成年电影在线观看| 欧美日韩一级在线毛片| netflix在线观看网站| 国产av在哪里看| 丝袜美腿诱惑在线| 国产精品 欧美亚洲| 国产亚洲av高清不卡| 18禁国产床啪视频网站| 国产精品美女特级片免费视频播放器 | 看片在线看免费视频| 亚洲成人久久性| 天堂√8在线中文| 亚洲色图 男人天堂 中文字幕| 欧美一级a爱片免费观看看 | 不卡一级毛片| 午夜福利免费观看在线| 国产单亲对白刺激| 欧美精品亚洲一区二区| 黄色女人牲交| 欧美色视频一区免费| 午夜免费鲁丝| 午夜影院日韩av| 一夜夜www| 婷婷丁香在线五月| 国产精品一区二区精品视频观看| 免费在线观看完整版高清| tocl精华| 无限看片的www在线观看| 亚洲成人精品中文字幕电影| 久久久久久久久久久久大奶| 怎么达到女性高潮| 日本五十路高清| 天天添夜夜摸| 国产精品秋霞免费鲁丝片| 九色国产91popny在线| 国产熟女xx| 国产精品久久久av美女十八| 欧美黑人精品巨大| 国产成人系列免费观看| 日韩精品免费视频一区二区三区| 亚洲无线在线观看| 午夜福利在线观看吧| 久久九九热精品免费| 欧美中文日本在线观看视频| 国产成+人综合+亚洲专区| 成人三级黄色视频| 欧美激情久久久久久爽电影 | 精品日产1卡2卡| 欧美日韩亚洲国产一区二区在线观看| 黑人巨大精品欧美一区二区mp4| 国产私拍福利视频在线观看| 免费高清视频大片| 亚洲中文av在线| 51午夜福利影视在线观看| 一级,二级,三级黄色视频| 国产av又大| 亚洲五月婷婷丁香| 久99久视频精品免费| 精品午夜福利视频在线观看一区| 脱女人内裤的视频| 免费在线观看视频国产中文字幕亚洲| 老司机深夜福利视频在线观看| 欧美日韩黄片免| 国产精品香港三级国产av潘金莲| videosex国产| 国产欧美日韩一区二区精品| 免费高清视频大片| 男女之事视频高清在线观看| 99久久国产精品久久久| 国产精品国产高清国产av| 欧美在线黄色| 中国美女看黄片| 成人国语在线视频| 欧美日韩精品网址| 亚洲成人国产一区在线观看| 国产精品 国内视频| 亚洲国产中文字幕在线视频| 欧美+亚洲+日韩+国产| 亚洲欧美一区二区三区黑人| 中文字幕久久专区| 久久久久久人人人人人| 大型av网站在线播放| 黄色丝袜av网址大全| 老汉色∧v一级毛片| 亚洲人成77777在线视频| 可以免费在线观看a视频的电影网站| 久热爱精品视频在线9| 长腿黑丝高跟| 成人亚洲精品av一区二区| 动漫黄色视频在线观看| 国产亚洲精品久久久久5区| 91麻豆av在线| 日本免费一区二区三区高清不卡 | 国产1区2区3区精品| 999精品在线视频| 国产精品二区激情视频| 国内毛片毛片毛片毛片毛片| 操出白浆在线播放| 黄色毛片三级朝国网站| 夜夜看夜夜爽夜夜摸| 男女午夜视频在线观看| 两个人免费观看高清视频| 亚洲无线在线观看| 青草久久国产| 操美女的视频在线观看| 窝窝影院91人妻| 国产精品98久久久久久宅男小说| 欧美日本亚洲视频在线播放| 男女床上黄色一级片免费看| 制服丝袜大香蕉在线| 国产色视频综合| 婷婷丁香在线五月| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利高清视频| 国产精品久久久久久人妻精品电影| 国产aⅴ精品一区二区三区波| 日本免费一区二区三区高清不卡 | av免费在线观看网站| 国产欧美日韩精品亚洲av| 97超级碰碰碰精品色视频在线观看| 欧美日韩精品网址| 免费高清在线观看日韩| 亚洲欧美日韩高清在线视频| 亚洲美女黄片视频| 欧美日韩一级在线毛片| 在线观看66精品国产| 在线免费观看的www视频| 欧美人与性动交α欧美精品济南到| 国产1区2区3区精品| 成人三级黄色视频| 一区二区三区精品91| 高潮久久久久久久久久久不卡| 黑人巨大精品欧美一区二区蜜桃| 母亲3免费完整高清在线观看| 国产激情欧美一区二区| 在线观看免费视频网站a站| 香蕉国产在线看| 丝袜美足系列| 18禁国产床啪视频网站| 欧美一级毛片孕妇| 成人国产一区最新在线观看| 久久精品91无色码中文字幕| 美国免费a级毛片| 午夜福利一区二区在线看| 午夜精品久久久久久毛片777| 欧美性长视频在线观看| 老汉色av国产亚洲站长工具| 国产成人啪精品午夜网站| 丝袜在线中文字幕| av网站免费在线观看视频| 中文字幕久久专区| 18禁国产床啪视频网站| 精品久久久久久成人av| 精品不卡国产一区二区三区| 欧美日韩亚洲综合一区二区三区_| 高清在线国产一区| 天天躁狠狠躁夜夜躁狠狠躁| 在线观看免费午夜福利视频| 精品不卡国产一区二区三区| 亚洲七黄色美女视频| 亚洲欧美日韩高清在线视频| 久久国产精品男人的天堂亚洲| 99香蕉大伊视频| 日韩精品中文字幕看吧| 久久久久久大精品| 欧美大码av| 国产亚洲精品第一综合不卡| 国产成人影院久久av| 亚洲av成人一区二区三| 亚洲av美国av| www国产在线视频色| 身体一侧抽搐| 亚洲国产欧美一区二区综合| 国产成+人综合+亚洲专区| 亚洲精品美女久久久久99蜜臀| 女人被狂操c到高潮| 精品国产亚洲在线| 亚洲五月天丁香| 91精品国产国语对白视频| 日韩精品免费视频一区二区三区| 精品国产一区二区久久| 久久精品91蜜桃| 亚洲七黄色美女视频| 波多野结衣巨乳人妻| av视频免费观看在线观看| 国产97色在线日韩免费| 每晚都被弄得嗷嗷叫到高潮| 日韩高清综合在线| 亚洲欧美精品综合一区二区三区| 欧美日韩乱码在线| 一级毛片精品| 老汉色av国产亚洲站长工具| 国产精品一区二区免费欧美| 1024香蕉在线观看| 国产精品久久电影中文字幕| 免费少妇av软件| 在线观看免费午夜福利视频| 麻豆久久精品国产亚洲av| 亚洲国产欧美一区二区综合| 欧美激情久久久久久爽电影 | 岛国在线观看网站| 丝袜美腿诱惑在线| 亚洲国产欧美网| 日韩欧美三级三区| 国产人伦9x9x在线观看| 午夜免费观看网址| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产精品999在线| 可以在线观看的亚洲视频| 少妇 在线观看| 大陆偷拍与自拍| 国产激情欧美一区二区| 一进一出好大好爽视频| 国内久久婷婷六月综合欲色啪| 老司机深夜福利视频在线观看| 女性生殖器流出的白浆| 免费在线观看日本一区| 亚洲熟妇中文字幕五十中出| av天堂久久9| 久久狼人影院| 亚洲免费av在线视频| 国产欧美日韩一区二区三区在线| 香蕉国产在线看| 一级黄色大片毛片| 午夜福利在线观看吧| 亚洲国产日韩欧美精品在线观看 | 女生性感内裤真人,穿戴方法视频| 麻豆av在线久日| 伦理电影免费视频| 久久精品91蜜桃| 国产精品免费一区二区三区在线| 国产亚洲av嫩草精品影院| 看免费av毛片| 国产成人欧美| 久久久久久国产a免费观看| 日本一区二区免费在线视频| 亚洲av片天天在线观看| 国产男靠女视频免费网站| 免费人成视频x8x8入口观看| av天堂久久9| 亚洲专区国产一区二区| 男女做爰动态图高潮gif福利片 | 91老司机精品| 国产99白浆流出| 99国产精品免费福利视频| 欧美乱妇无乱码| 国内精品久久久久精免费| 亚洲伊人色综图| 日日干狠狠操夜夜爽| 欧美中文综合在线视频| 少妇熟女aⅴ在线视频| 国产精品,欧美在线| 色av中文字幕| 国产精品香港三级国产av潘金莲| 可以免费在线观看a视频的电影网站| 久久久久亚洲av毛片大全| 99香蕉大伊视频| 亚洲精品中文字幕在线视频| 久9热在线精品视频| 变态另类成人亚洲欧美熟女 | 999精品在线视频| 在线观看免费日韩欧美大片| 国产精品1区2区在线观看.| 少妇裸体淫交视频免费看高清 | 欧美色欧美亚洲另类二区 | 桃色一区二区三区在线观看| 亚洲美女黄片视频| 成年女人毛片免费观看观看9| 久久久久亚洲av毛片大全| 操出白浆在线播放| 97碰自拍视频| 成人手机av| 国产精品亚洲一级av第二区| 亚洲人成电影观看| 中文亚洲av片在线观看爽| 免费看美女性在线毛片视频| 9色porny在线观看| 亚洲五月色婷婷综合| 日韩欧美国产在线观看| 亚洲中文字幕日韩| 免费高清在线观看日韩| 色尼玛亚洲综合影院| 日韩免费av在线播放| 中国美女看黄片| 热re99久久国产66热| 久热这里只有精品99| tocl精华| 日韩国内少妇激情av| 国产精品免费一区二区三区在线| 桃红色精品国产亚洲av| 中文字幕高清在线视频| 久久国产乱子伦精品免费另类| 国产伦人伦偷精品视频| 正在播放国产对白刺激| 亚洲电影在线观看av| 在线永久观看黄色视频| 欧美日韩一级在线毛片| 久久精品亚洲熟妇少妇任你| 色综合亚洲欧美另类图片| 亚洲一卡2卡3卡4卡5卡精品中文| 熟妇人妻久久中文字幕3abv| 国产精品,欧美在线| 女人被狂操c到高潮| 亚洲欧美激情在线| 午夜福利在线观看吧| 亚洲五月色婷婷综合| a在线观看视频网站| 日本黄色视频三级网站网址| 69精品国产乱码久久久| 亚洲专区中文字幕在线| 一进一出抽搐动态| 波多野结衣高清无吗| 两人在一起打扑克的视频| 日本五十路高清| 十分钟在线观看高清视频www| 黄网站色视频无遮挡免费观看| 久久人妻熟女aⅴ| 欧美黄色淫秽网站| 麻豆国产av国片精品| 黄色毛片三级朝国网站| 在线视频色国产色| 亚洲精华国产精华精| 国产熟女xx| 亚洲精品美女久久av网站| 69精品国产乱码久久久| 亚洲自偷自拍图片 自拍| 两个人免费观看高清视频| 日韩av在线大香蕉| 精品欧美一区二区三区在线| 9热在线视频观看99| 99re在线观看精品视频| 精品免费久久久久久久清纯| 日韩一卡2卡3卡4卡2021年| 一边摸一边做爽爽视频免费| 男女床上黄色一级片免费看| 99香蕉大伊视频| 不卡av一区二区三区| 国产精品日韩av在线免费观看 | av免费在线观看网站| 免费女性裸体啪啪无遮挡网站| 看黄色毛片网站| 日本精品一区二区三区蜜桃| 午夜成年电影在线免费观看| 婷婷精品国产亚洲av在线| 午夜精品久久久久久毛片777| 久久香蕉精品热| 最近最新免费中文字幕在线| 又黄又粗又硬又大视频| 国产三级黄色录像| 50天的宝宝边吃奶边哭怎么回事| 免费av毛片视频| 国产激情久久老熟女| 成人三级黄色视频| 操美女的视频在线观看| 天天躁夜夜躁狠狠躁躁| 99re在线观看精品视频| 国产麻豆69| 一进一出抽搐gif免费好疼| 极品教师在线免费播放| 操美女的视频在线观看| 国产一区二区三区视频了| 国产精品久久久久久人妻精品电影| 制服人妻中文乱码| av在线播放免费不卡| 午夜福利视频1000在线观看 | 久久精品影院6| 亚洲av电影不卡..在线观看| 一边摸一边做爽爽视频免费| 啦啦啦 在线观看视频| 99久久国产精品久久久| 欧美乱色亚洲激情| 国产一卡二卡三卡精品| 女人精品久久久久毛片| 嫩草影视91久久| 99在线视频只有这里精品首页| 国产99久久九九免费精品| 无遮挡黄片免费观看| 如日韩欧美国产精品一区二区三区| 久久青草综合色| 激情视频va一区二区三区| 精品国产乱码久久久久久男人| √禁漫天堂资源中文www| 亚洲精品av麻豆狂野| 国产一区在线观看成人免费| 可以在线观看的亚洲视频| 欧美日韩瑟瑟在线播放| av在线播放免费不卡| 国产精华一区二区三区| 又紧又爽又黄一区二区| 国产伦人伦偷精品视频| 日本黄色视频三级网站网址| 亚洲成人久久性| 在线观看午夜福利视频| 999久久久国产精品视频| 欧美乱妇无乱码| 男女午夜视频在线观看| 国产单亲对白刺激| 亚洲男人天堂网一区| 黑人欧美特级aaaaaa片| 免费看a级黄色片| 免费在线观看黄色视频的| 在线观看舔阴道视频| 国产午夜福利久久久久久| 国产成人精品在线电影| 日本撒尿小便嘘嘘汇集6| 国产精品一区二区在线不卡| 每晚都被弄得嗷嗷叫到高潮| 亚洲成人国产一区在线观看| 国产一区二区三区在线臀色熟女| av在线天堂中文字幕| 国产成人欧美在线观看| 免费无遮挡裸体视频| 黄频高清免费视频| 国产人伦9x9x在线观看| 久久亚洲精品不卡| 亚洲一卡2卡3卡4卡5卡精品中文| 黄色成人免费大全| 18禁国产床啪视频网站| svipshipincom国产片| 亚洲精品中文字幕在线视频| 91国产中文字幕| 黑人巨大精品欧美一区二区蜜桃| 法律面前人人平等表现在哪些方面| 一a级毛片在线观看| 脱女人内裤的视频| 精品乱码久久久久久99久播| 国产亚洲精品综合一区在线观看 | 色尼玛亚洲综合影院| 国产精品九九99| 手机成人av网站| av视频在线观看入口| 999精品在线视频| 麻豆av在线久日| 亚洲全国av大片| 女人被狂操c到高潮| 国产野战对白在线观看| 久久人妻av系列| 精品一品国产午夜福利视频| 91老司机精品| 亚洲av第一区精品v没综合| 手机成人av网站| 国产极品粉嫩免费观看在线| 国产精品免费视频内射| 制服诱惑二区| 亚洲国产精品sss在线观看| av福利片在线| 国产熟女xx| 女生性感内裤真人,穿戴方法视频| 久久精品国产99精品国产亚洲性色 | 久久性视频一级片| 国产日韩一区二区三区精品不卡| 亚洲欧美日韩另类电影网站| 91av网站免费观看| 老司机午夜福利在线观看视频| 亚洲国产日韩欧美精品在线观看 | 999久久久国产精品视频| 免费女性裸体啪啪无遮挡网站| 欧美乱码精品一区二区三区| 国产主播在线观看一区二区| 欧美日韩亚洲国产一区二区在线观看| 久久婷婷成人综合色麻豆| 女生性感内裤真人,穿戴方法视频| 99精品在免费线老司机午夜| 久久久国产欧美日韩av| 大码成人一级视频| 妹子高潮喷水视频| 后天国语完整版免费观看| 中文字幕高清在线视频| 久久久久久大精品| 亚洲情色 制服丝袜| 国产精品一区二区三区四区久久 | 精品国产国语对白av| 久久精品人人爽人人爽视色| 国产亚洲av高清不卡| 两人在一起打扑克的视频| 在线观看免费视频网站a站| 国产成人精品久久二区二区免费| 在线天堂中文资源库| 亚洲午夜精品一区,二区,三区| 国产成人一区二区三区免费视频网站| 天天躁狠狠躁夜夜躁狠狠躁| 脱女人内裤的视频| 日韩av在线大香蕉| 国产精品98久久久久久宅男小说| 日本撒尿小便嘘嘘汇集6| 宅男免费午夜| 可以在线观看毛片的网站| cao死你这个sao货| 在线十欧美十亚洲十日本专区| 欧美激情高清一区二区三区| 一级毛片精品| 欧美黄色片欧美黄色片| 色综合婷婷激情| 国产午夜福利久久久久久| 一区二区日韩欧美中文字幕| 日本 欧美在线| 最新美女视频免费是黄的| 每晚都被弄得嗷嗷叫到高潮| 午夜福利免费观看在线| 一区在线观看完整版| 在线国产一区二区在线| 日韩中文字幕欧美一区二区| av在线播放免费不卡| 非洲黑人性xxxx精品又粗又长| 一本大道久久a久久精品| 日本 av在线|