999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于旅游資源整合的高速公路網(wǎng)絡(luò)雙層規(guī)劃模型及實例研究

2015-11-04 10:28:36王文輝朱曉英
城市道橋與防洪 2015年1期
關(guān)鍵詞:高速公路規(guī)劃旅游

王文輝,朱曉英

(1.同濟(jì)大學(xué)土木工程學(xué)院,上海市200092;2.天津城建設(shè)計院有限公司,天津市300121)

基于旅游資源整合的高速公路網(wǎng)絡(luò)雙層規(guī)劃模型及實例研究

王文輝1,朱曉英2

(1.同濟(jì)大學(xué)土木工程學(xué)院,上海市200092;2.天津城建設(shè)計院有限公司,天津市300121)

將高速公路資源整合到旅游交通網(wǎng)絡(luò)中,是當(dāng)前高速公路網(wǎng)絡(luò)規(guī)劃的趨勢和方向?;诖?,從旅游資源整合的角度出發(fā),建立了基于旅游資源的高速公路網(wǎng)絡(luò)雙層規(guī)劃模型;然后,研究了模型的求解方法,采用遺傳算法,并用Mat1ab編寫了求解的算法程序;最后,給出黑龍江省高速公路網(wǎng)絡(luò)規(guī)劃的實例研究,論證了模型的可行性和合理性。

旅游資源;雙層規(guī)劃;遺傳算法;黑龍江

0 引言

隨著旅游產(chǎn)業(yè)發(fā)展的不斷深入,公路作為旅游活動的主要載體備受關(guān)注,成為旅游研究的重要方面。旅游交通作為旅游系統(tǒng)的重要組成部分,是區(qū)域旅游業(yè)發(fā)展的重要支撐和保障[1]。其服務(wù)質(zhì)量和水平影響著游客的旅游決策、旅游目的地的選擇以及旅游資源的開發(fā)建設(shè),進(jìn)而影響到區(qū)域的旅游開發(fā)。一方面,交通為旅游的發(fā)展提供了必要條件;另一方面,旅游對交通的發(fā)展也起到了相當(dāng)大的促進(jìn)作用[2]。

高速公路作為社會經(jīng)濟(jì)發(fā)展到一定階段的產(chǎn)物,具有縮短行程時間、節(jié)約花費、便捷舒適、安全等優(yōu)點。隨著居民生活水平的提高、私家車保有量的增長以及人們對出行路網(wǎng)環(huán)境需求的提升,將高速公路資源整合到旅游交通網(wǎng)絡(luò)中,建立合理、高效、便捷的交通系統(tǒng),是當(dāng)前區(qū)域路網(wǎng)規(guī)劃的趨勢和方向,也是促進(jìn)旅游景區(qū)的開發(fā)與資源整合的必然需求。實踐證明,旅游景點附近高速公路交通網(wǎng)絡(luò)的建設(shè),能使依靠公路交通運輸?shù)穆糜物L(fēng)景區(qū)輻射半徑擴(kuò)大320~500 km[3]。

近年來,公路網(wǎng)絡(luò)優(yōu)化的主要方法是實證調(diào)查、數(shù)理分析方法、系統(tǒng)分析、四階段方法和網(wǎng)絡(luò)設(shè)計問題(NDP)。網(wǎng)絡(luò)設(shè)計問題是指在約束條件下選擇最優(yōu)的投資決策,通過新建路段或改擴(kuò)建路段來使整個公路網(wǎng)絡(luò)的性能達(dá)到最優(yōu),滿足遠(yuǎn)期交通增長的需求[4-5]。對于研究區(qū)域的公路網(wǎng),NDP能使網(wǎng)絡(luò)布局工作更加定量化和科學(xué)化,布局優(yōu)化方案很多,在方案選擇上也更趨于多元化,如何最大限度地滿足公路交通需求成為路線優(yōu)化的重點[6-8]。網(wǎng)絡(luò)設(shè)計中的決策變量不同決定了其模型也會不一樣,決策變量主要有新建路段和改建路段提升能力兩種,同時NDP也被分為三類:離散性、連續(xù)性和混合性三種,分別對應(yīng)不同的決策變量[9]。由于雙層規(guī)劃能夠充分解釋網(wǎng)絡(luò)設(shè)計問題的決策過程,因此Gwo-Hshiung Tzeng利用多目標(biāo)數(shù)學(xué)規(guī)劃來設(shè)計連續(xù)性網(wǎng)絡(luò)設(shè)計模型,利用雙層規(guī)劃的思想解釋并討論帶有多目標(biāo)決策的網(wǎng)絡(luò)改進(jìn)問題,建立了非線性雙層規(guī)劃模型[10]。

在以往的省域高速公路網(wǎng)絡(luò)規(guī)劃時,一般只考慮主要城市的連通性,未將旅游景區(qū)的連通性納入其中。隨著旅游業(yè)的不斷發(fā)展以及旅游景區(qū)服務(wù)需求的提高,高速旅游公路的規(guī)劃將成為必然?;诖?,本文從高速公路網(wǎng)與旅游資源整合的角度進(jìn)行研究,建立基于旅游資源的高速公路網(wǎng)絡(luò)雙層規(guī)劃模型,并給出了以黑龍江省高速公路網(wǎng)絡(luò)為例的實證研究。

1 問題描述

在進(jìn)行基于旅游資源的公路網(wǎng)絡(luò)規(guī)劃時,需要考慮多方面的因素,既要考慮到政府部門的規(guī)劃決策,又要考慮到出行者的出行需求。從網(wǎng)絡(luò)規(guī)劃者(政府部門)的角度考慮,希望整個網(wǎng)絡(luò)系統(tǒng)的出行總費用最少,同時,在將高速公路網(wǎng)與旅游資源整合時,希望規(guī)劃后的路網(wǎng)旅游資源需求增加最大。從網(wǎng)絡(luò)用戶(出行者)的角度考慮,每個用戶都希望出行的時間或者費用最少。綜合考慮,為了使模型建立更接近于真實情況,單一的優(yōu)化目標(biāo)不能滿足此規(guī)劃的要求,因此,本文采用雙層規(guī)劃模型。雙層規(guī)劃模型中上層從網(wǎng)絡(luò)規(guī)劃者(政府部門)的角度考慮,下層規(guī)劃模型從網(wǎng)絡(luò)用戶(出行者)的角度考慮。

2 模型建立

2.1優(yōu)化目標(biāo)

基于旅游資源的高速公路網(wǎng)絡(luò)規(guī)劃從微觀分析入手,建立雙層規(guī)劃模型。上層模型的優(yōu)化目標(biāo)為:整個網(wǎng)絡(luò)系統(tǒng)的出行總費用(時間)最少,路網(wǎng)旅游資源需求增加最大。下層模型的優(yōu)化目標(biāo)為:單個用戶的出行費用(時間)最少。

其中路網(wǎng)旅游資源需求如何量化是目標(biāo)函數(shù)的關(guān)鍵。通常計算旅游資源需求時,會將旅游景區(qū)的現(xiàn)有規(guī)模、開發(fā)規(guī)劃、服務(wù)水平以及旅游者的消費水平等作為考慮因素,計算出的需求增加是由于旅游景區(qū)的開發(fā)等因素導(dǎo)致的。本文衡量的是規(guī)劃路網(wǎng)建成后對于旅游資源需求的影響,可以通過交通量的增加來反映,在此提出旅游交通誘增量的概念。計算分析規(guī)劃路網(wǎng)建成后由于路網(wǎng)結(jié)構(gòu)的改善等而增加的旅游交通量。

在計算交通誘增量時,一般引用重力模型進(jìn)行計算[11]。本文提出的旅游交通誘增量為交通誘增量的一個分支,因此重力模型也適用于旅游交通誘增量的計算。借鑒張航提出的基于重力模型預(yù)測交通量的方法[12],給出旅游交通誘增量模型。

式中:qNij為i區(qū)到j(luò)區(qū)間的旅游交通出行量;Qij為 i區(qū)到j(luò)區(qū)的旅游趨勢交通量分布,veh/d;Q′ij為i區(qū)到j(luò)區(qū)的旅游誘增交通量,veh/d;ta為無此項目時(現(xiàn)狀)i區(qū)到j(luò)區(qū)的交通阻抗(時間);ta′為有此項目時(將來),i區(qū)到j(luò)區(qū)的交通阻抗(時間);K、a、β、γ為重力模型參數(shù)。

2.2符號定義

其中,xa,ya為決策變量,是0-1變量,xa為新建路段時,其值為1,否則為0;ya為擴(kuò)建路段時,其值為1,否則為0。va(xa,ya)為路段a上的交通量;ta(va(xa,ya),xa,ya)為路段a上的出行時間;minz(xa,ya)為公路網(wǎng)的總出行時間;manFxa,ya)為基于旅游需求的總誘增交通量。

在約束條件中,需要考慮路網(wǎng)造價、規(guī)劃里程、服務(wù)水平等約束。其中,m1為新建路段a的單車道造價,萬元/km;m2為擴(kuò)建路段a的造價,萬元/km;M為預(yù)期建設(shè)資金的最大值,常數(shù);la為新建路段a的里程,km;lyy為規(guī)劃前(即原有)公路網(wǎng)里程,km;La為公路網(wǎng)規(guī)劃總里程;ωa為規(guī)劃路段a的服務(wù)水平,即V/C(路段交通量與路段通行能力之比);ωmin為最小允許服務(wù)水平;ωman為最大允許服務(wù)水平。

2.3雙層規(guī)劃模型形式

由上述可以建立公路網(wǎng)絡(luò)路線布局的雙層規(guī)劃模型。

其中,優(yōu)化目標(biāo)①、②為上層優(yōu)化目標(biāo)函數(shù),①表示基于旅游需求的總誘增交通量最大,②表示路網(wǎng)系統(tǒng)的總時間最少。優(yōu)化目標(biāo)③為下層優(yōu)化目標(biāo)函數(shù),③表示單個用戶出行時間最少。其實質(zhì)是路網(wǎng)的平衡配流模型。本文采用Beckmannde用戶平衡分配模型。此模型較為成熟并且應(yīng)用頗多,在此不再贅述。

約束條件中①表示高速公路網(wǎng)絡(luò)規(guī)模在最佳規(guī)模內(nèi);②表示路段服務(wù)水平(V/C)在一定范圍內(nèi);③表示規(guī)劃總投資在預(yù)期建設(shè)資金內(nèi);④表示交通流守恒;⑤表示路段流量守恒。

3 模型求解算法

3.1上層優(yōu)化模型目標(biāo)函數(shù)的轉(zhuǎn)化

在本文上層優(yōu)化模型中,有兩個目標(biāo)函數(shù),應(yīng)先轉(zhuǎn)化為單目標(biāo)函數(shù)。將多目標(biāo)函數(shù)轉(zhuǎn)化為單目標(biāo)函數(shù)一般可采用加法或乘除法。本文中上層優(yōu)化模型中的兩個目標(biāo)為:基于旅游需求的總誘增交通量最大,路網(wǎng)系統(tǒng)的運行時間最小,兩目標(biāo)處于沖突狀態(tài)。因此,本文適合用乘除法進(jìn)行目標(biāo)函數(shù)的轉(zhuǎn)化,轉(zhuǎn)化后的單目標(biāo)函數(shù)為:

3.2遺傳算法求解步驟

類似上述交通網(wǎng)絡(luò)的雙層規(guī)劃模型常用的解法通常有Bender分解法、分枝定界法、極點算法、互補旋轉(zhuǎn)算法、模擬退火算法以及一些啟發(fā)式算法求解[13]。由于模型的上層和下層均為0-1混合整數(shù)規(guī)劃問題,比較適合采用遺傳算法進(jìn)行求解。

遺傳算法的應(yīng)用目前比較成熟且較為廣泛,因此有關(guān)算法的原理、迭代步驟等在此不再贅述。本文對遺傳算法的求解步驟進(jìn)行改進(jìn)。

考慮現(xiàn)狀路網(wǎng)結(jié)構(gòu)、規(guī)劃年0D交通量以及專家意見等因素,生成初始可能連接圖。

首先,進(jìn)行交通流分配,計算各個體對應(yīng)的目標(biāo)函數(shù)值,對個體按目標(biāo)函數(shù)值進(jìn)行排序。將目標(biāo)函數(shù)值不大于當(dāng)前最優(yōu)值的個體轉(zhuǎn)入下層規(guī)劃進(jìn)行交通流分配,計算路段交通量和負(fù)荷度V/C。如果Gen>Maxgen或stop>Maxstop,結(jié)束計算;否則,Gen=Gen+1,根據(jù)交通流分配結(jié)果檢查上層約束條件,將不滿足條件的個體省去。然后,計算種群的個體適應(yīng)度函數(shù)值,利用“無回放隨機(jī)選擇法”進(jìn)行選擇操作,依概率進(jìn)行交叉、變異操作,交叉采用兩點交叉操作,轉(zhuǎn)向第1步。結(jié)束,輸出結(jié)果,得到最終的路網(wǎng)優(yōu)化方案。

其中,設(shè)置初始進(jìn)化代數(shù)Gen=0,最優(yōu)值維持步長stop=0;設(shè)置初始目標(biāo)值Zmax為一相當(dāng)大的正數(shù),隨機(jī)設(shè)定一個由初始可能連接圖的子圖構(gòu)成的種群作為上層規(guī)劃的初始方案。群體規(guī)模取10,交叉概率和變異概率分別取0.9和0.1,最大進(jìn)化代數(shù)取100。利用Mat1ab軟件及其提供的遺傳算法工具箱編寫程序。

4 算例分析

為了方便對雙層規(guī)劃模型的計算結(jié)果進(jìn)行分析,以黑龍江省現(xiàn)有高速公路網(wǎng)絡(luò)為例,來驗證模型的可行性和合理性。

4.1方案描述

為了更好地反映實際規(guī)劃需求,本文給出了4個優(yōu)化方案分別進(jìn)行模型的求解,進(jìn)行對比分析,選取最優(yōu)的優(yōu)化結(jié)果。

現(xiàn)階段黑龍江省大部分旅游景區(qū)并未在高速公路網(wǎng)絡(luò)節(jié)點內(nèi),即旅游景區(qū)與高速公路網(wǎng)絡(luò)是不連通的,這會對模型的求解造成困難。因此,在設(shè)計方案之初,首先要連接一些旅游高速公路,從而保證路網(wǎng)的連通性。本文將5個旅游景區(qū)(亞布力、五大連池、杜爾伯特、鏡泊湖、饒河縣)納入路網(wǎng)節(jié)點,除亞布力外,其他景區(qū)均與高速公路網(wǎng)絡(luò)不連通,因此,需要在現(xiàn)有路網(wǎng)基礎(chǔ)上增加4條旅游高速公路作為設(shè)計方案。

將高速公路網(wǎng)絡(luò)與5個旅游景區(qū)連接可以有很多種方案,考慮到設(shè)計旅游高速公路線路的人為因素,僅做出一種設(shè)計方案進(jìn)行路網(wǎng)優(yōu)化的計算是不合理的、片面的。由于篇幅的限制,本文不能給出所有旅游高速公路連接的可能性,在綜合專家建議、選線合理性等情況的基礎(chǔ)上,本文給出了3個設(shè)計方案。依據(jù)上文提出的模型,將3個設(shè)計方案進(jìn)行網(wǎng)絡(luò)優(yōu)化,選擇出最優(yōu)的路網(wǎng)優(yōu)化方案,以期為實際的網(wǎng)絡(luò)設(shè)計提供參考。

同時,給出方案4作為基本方案,僅將15個主要城市作為路網(wǎng)節(jié)點,不考慮5個旅游景區(qū)的情況,與方案1~3進(jìn)行對比分析。

4.2模型參數(shù)取值

根據(jù)《黑龍江省骨架路網(wǎng)規(guī)劃》的公路遠(yuǎn)景規(guī)劃,2020年公路建設(shè)總投資約600億元,高速公路總里程將達(dá)到7 000 km。單位里程造價新建道路2 000萬元/km,擴(kuò)建道路800萬元/km。上層模型中路段負(fù)荷度約束上限取為1,下限取0。

4.3數(shù)據(jù)信息

(1)路網(wǎng)節(jié)點編碼

將黑龍江省內(nèi)的15個主要城市(包括地級市)與5個旅游景區(qū)進(jìn)行編號,節(jié)點編號如表1所示。

表1 路網(wǎng)節(jié)點編碼

(2)黑龍江省高速公路網(wǎng)絡(luò)距離矩陣

根據(jù)《黑龍江公路水路交通“十一五”規(guī)劃》和《黑龍江省骨架路網(wǎng)規(guī)劃》以及專家意見,在現(xiàn)狀公路網(wǎng)的基礎(chǔ)上,給出每種方案的路網(wǎng)距離矩陣。矩陣中包括現(xiàn)有路段和備選路段,現(xiàn)有路段的距離為已知數(shù)據(jù),備選路段的距離根據(jù)專家選線的方案得出。

(3)黑龍江省內(nèi)小區(qū)劃分及編碼

為了進(jìn)行交通流分配,以各市政府所在地為形心,所管轄范圍為小區(qū)范圍,將路網(wǎng)劃分為8個交通小區(qū),各小區(qū)編號如表2所示。圖1為黑龍江省0D交通小區(qū)劃分示意圖。

表2 黑龍江省OD小區(qū)編號

(4)黑龍江省高速公路網(wǎng)絡(luò)機(jī)動車0D矩陣

根據(jù)《黑龍江省骨架路網(wǎng)規(guī)劃》及《黑龍江公路水路交通“十一五”規(guī)劃》中的公路運輸量需求分析以及主要公路通道交通量預(yù)測結(jié)果進(jìn)行反推,預(yù)測各節(jié)點交通生成量后進(jìn)行交通分布,得到2020年各小區(qū)之間的機(jī)動車0D交換量的預(yù)測結(jié)果,如表3所示。但是由于調(diào)查方式及數(shù)據(jù)的局限性,所計算出的機(jī)動車0D交換量存在一定的誤差,僅為此案例做以計算分析使用。

圖1 黑龍江省OD交通小區(qū)劃分示意圖

表3 黑龍江省高速公路網(wǎng)絡(luò)機(jī)動車OD表(單位:100pcu)

4.4模型優(yōu)化結(jié)果

將所有的數(shù)據(jù)信息錄入,運用Mat1ab程序進(jìn)行運算,得出4種方案的優(yōu)化結(jié)果。將4種方案的優(yōu)化結(jié)果進(jìn)行對比分析,如表4所示。

表4 四種優(yōu)化方案結(jié)果比較

方案4的總建設(shè)費用較少,但是由于其只考慮了主要城市的連通性,而沒有將主要旅游景區(qū)納入其內(nèi),不能更好地發(fā)揮高速公路網(wǎng)在旅游資源優(yōu)化上的優(yōu)勢,只作為對比方案使用,不能作為最優(yōu)方案。

方案1、2、3均為納入旅游景區(qū)為節(jié)點的高速公路網(wǎng)絡(luò)優(yōu)化,雖然總建設(shè)費用較多,但都在其約束條件的范圍內(nèi)。不僅能夠滿足整個路網(wǎng)總通行時間最小的優(yōu)化目標(biāo),還能滿足提高旅游誘增交通量的目標(biāo)。方案1共新建了8條連接旅游景區(qū)的高速公路,分別為:杜爾伯特—大慶、五大連池—杜爾伯特、五大連池—北安、饒河縣—佳木斯、亞布力—佳木斯、饒河縣—同江、亞布力—鏡泊湖、鏡泊湖—牡丹江。方案2共新建7條連接旅游景區(qū)的高速公路,分別為:五大連池—北安、杜爾伯特—齊齊哈爾、鏡泊湖—牡丹江、饒河縣—佳木斯、五大連池—杜爾伯特、亞布力—佳木斯、亞布力—鏡泊湖。方案3共新建8條連接旅游景區(qū)的高速公路,分別為:五大連池—北安、杜爾伯特—大慶、鏡泊湖—牡丹江、饒河縣—同江、杜爾伯特—齊齊哈爾、五大連池—杜爾伯特、亞布力—佳木斯、亞布力—鏡泊湖,且方案3的總建設(shè)費用比方案2少27.36億元,比方案1少42.64億元,大大減少了投資費用。綜合比較,選擇方案3作為最終的優(yōu)化方案。

5 結(jié)語

高速公路網(wǎng)絡(luò)與旅游資源整合的發(fā)展策略將成為未來公路網(wǎng)規(guī)劃的一個趨勢和方向。本文建立了基于旅游資源的高速公路網(wǎng)絡(luò)雙層規(guī)化模型,采用遺傳算法進(jìn)行求解,并用Mat1ab軟件及其工具箱編寫了算法程序。以黑龍江省高速公路網(wǎng)絡(luò)優(yōu)化作為實例分析,驗證了模型的可行性。

[1]黃琳, 金海龍, 包瑞. 新疆旅游交通現(xiàn)狀及發(fā)展研究[J]. 新疆師范大學(xué)學(xué)報(自然科學(xué)版),2008,27(1):103-106.

[2]孫有望, 李云清.論旅游交通與交通旅游[J].上海鐵道大學(xué)學(xué)報,1999,20(10):65-69.

[3]歐越男, 謝光秋. 高速公路建設(shè)與發(fā)展旅游經(jīng)濟(jì)的可行性探討[J]. 交通與運輸,2005(12):60-62.

[4]Luis Gouveia, Jose Manue1 Pires. Mode1s for a Steiner ring network design Prob1em with revenues[J]. EuroPean Jouma1 of 0peration Researeh,2001,(133): 21-31.

[5]Me1kote Daskin M. S. An intenratedmode1 of faci1ity 1oeation and transportation network design[J]. TransPortation Research PartA,2001, 35(6): 515-538.

[6]Magnanti T. L. and Wong R. T. Network design and transportation p1anning:mode1s and a1gorithms[J]. TransPortation Science,1984,(18): 1-55.

[7]Daganzo C. Queue spi11overs in transportation networks with a route choice[J]. Transportation Science,1998(32): 3-11.

[8]Chen,M.and A1fa,A.S.A network design a1gorithm using astochsatic incrementa1 traffic assignment approach[J]. Transportation Science,1991,(25): 215-224.

[9]JohnsonD.S,Lenstra J.K.and Rinnooy Kan A.H.G.The comp1exity of the network design Prob1em[J]. Network,1978(8): 279-285.

[10]Chen Y W, Tzeng G H. A fuzzymu1ti-objectivemode1 for reconstructing the post-quake road-network by genetic a1gorithm[J]. Internationa1 Journa1 of Fuzzy Systems, 1999, 1(2): 85-95.

[11]王兵, 董志宏,曹中杰.誘增交通量理論研究[J].河北工業(yè)大學(xué)學(xué)報,2002,31(6):102-106.

[12]張航, 張玲.基于重力模型預(yù)測誘增交通量方法研究.公路交通技術(shù),2006,2(1):111-113.

[13]雷英杰, 張善文,等.MATLAB遺傳算法工具箱及應(yīng)用[M].西安電子科技大學(xué)出版社,2005.

U412.36+6

A

1009-7716(2015)01-0133-05

2014-11-20

王文輝(1971-),男,河北石家莊人,高級工程師,在職研究生,研究方向為建筑與土木工程,現(xiàn)從事市政能源、碼頭航道等工作。

猜你喜歡
高速公路規(guī)劃旅游
規(guī)劃引領(lǐng)把握未來
快遞業(yè)十三五規(guī)劃發(fā)布
商周刊(2017年5期)2017-08-22 03:35:26
旅游
多管齊下落實規(guī)劃
高速公路與PPP
迎接“十三五”規(guī)劃
旅游的最后一天
高速公路上的狗
小說月刊(2014年4期)2014-04-23 08:52:20
GPS在高速公路中的應(yīng)用
河南科技(2014年18期)2014-02-27 14:15:06
出國旅游的42個表達(dá)
海外英語(2013年8期)2013-11-22 09:16:04
主站蜘蛛池模板: 成人午夜视频网站| 欧美成人第一页| 九九久久99精品| 久久99国产乱子伦精品免| 久久精品丝袜| 久久久精品无码一二三区| 欧美综合成人| 亚洲第一成人在线| 色综合天天操| 呦系列视频一区二区三区| 欧美亚洲日韩中文| 青青操视频在线| 婷婷亚洲最大| 国产成人区在线观看视频| 日韩a级毛片| 在线综合亚洲欧美网站| 亚洲精品你懂的| 午夜福利视频一区| 666精品国产精品亚洲| 国产99视频在线| 波多野结衣在线一区二区| 国产一级视频在线观看网站| 国产靠逼视频| 欧美日韩中文国产| 91久久性奴调教国产免费| 91精品视频网站| 韩国v欧美v亚洲v日本v| 亚洲 欧美 偷自乱 图片| 亚洲另类色| 欧美日韩中文字幕在线| 永久天堂网Av| 免费一极毛片| 国产成人免费手机在线观看视频| 在线人成精品免费视频| 九九九久久国产精品| 亚洲无码熟妇人妻AV在线| 手机精品视频在线观看免费| 欧美在线精品怡红院| 国产免费怡红院视频| 亚洲成a人片在线观看88| 中文字幕第4页| 永久在线播放| a天堂视频在线| 欧美成人在线免费| 久久国语对白| 国产视频a| 国产精品午夜福利麻豆| 久久婷婷色综合老司机| 久久国产精品麻豆系列| 高清久久精品亚洲日韩Av| 强奷白丝美女在线观看| 国产精品2| 色婷婷亚洲综合五月| 欧美精品啪啪| 国产小视频免费观看| 久久永久免费人妻精品| 亚洲国产天堂久久综合226114| 激情综合网激情综合| 精品免费在线视频| 亚洲一区二区在线无码| 在线人成精品免费视频| 一本大道香蕉高清久久| 日韩欧美国产成人| 国产乱码精品一区二区三区中文| 欧美视频二区| 黄色网页在线观看| 亚洲免费福利视频| 欧美亚洲另类在线观看| 中文字幕在线视频免费| 成人精品亚洲| 在线免费看片a| 久久国产精品娇妻素人| 亚洲日韩Av中文字幕无码| 国产精选小视频在线观看| 国内精品视频| 91外围女在线观看| 國產尤物AV尤物在線觀看| 精品五夜婷香蕉国产线看观看| 久久久久亚洲精品成人网 | 男女男免费视频网站国产| 国产精品无码久久久久AV| 中文字幕va|