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

    基于GIS 空間分析的城市空間增長(zhǎng)模擬與應(yīng)用

    2022-06-07 02:39:36房浩然
    經(jīng)緯天地 2022年2期
    關(guān)鍵詞:城市化用地預(yù)測(cè)

    房浩然

    (山東科技大學(xué),山東 泰安 271000)

    0.引言

    近年來(lái),我國(guó)城市化發(fā)展進(jìn)入了全新的發(fā)展周期,總體方向從量的無(wú)序擴(kuò)張,轉(zhuǎn)變?yōu)橄蚋哔|(zhì)量發(fā)展。新城區(qū)的發(fā)展需解決舊城區(qū)滯后發(fā)展帶來(lái)的弊端,同時(shí)需創(chuàng)造新的動(dòng)力,為新的經(jīng)濟(jì)發(fā)展提供著力點(diǎn),構(gòu)建全新的、具有多方位綜合性功能的新城區(qū)。新城區(qū)的發(fā)展使土地效益、經(jīng)濟(jì)效益得到顯著的提高。但出現(xiàn)了不少問(wèn)題,如,空間分布泛化、經(jīng)濟(jì)增長(zhǎng)速變緩、增長(zhǎng)極核作用減弱、低效用地較多、產(chǎn)區(qū)與城區(qū)分離問(wèn)題突出。因此,探索城市空間擴(kuò)展因素,預(yù)測(cè)城市空間未來(lái)拓展方向,已成為現(xiàn)階段城市發(fā)展研究的重要課題。

    在城市空間增長(zhǎng)模擬預(yù)測(cè)領(lǐng)域,國(guó)外最早在上世紀(jì)60年代,結(jié)合地區(qū)宗教、市場(chǎng)等客觀環(huán)境提出了理想化的靜態(tài)預(yù)測(cè)模型,后結(jié)合城市空間機(jī)制、外部空間拓展研究對(duì)靜態(tài)模型進(jìn)行了優(yōu)化,模型趨于理性化。在此基礎(chǔ)上進(jìn)行了動(dòng)態(tài)模型的探索,進(jìn)入了新城市化發(fā)展階段,在新城市化的基礎(chǔ)上融入了生態(tài)足跡理論,使動(dòng)態(tài)模型進(jìn)入成熟階段,構(gòu)建了現(xiàn)在廣泛使用的經(jīng)典CA模型,并在CA模型的基礎(chǔ)上進(jìn)行了非柵格的CA模型、社會(huì)經(jīng)濟(jì)CA模型、概率統(tǒng)計(jì)CA模型、ANNCA模型、SLEUTH模型的衍生[1]。國(guó)內(nèi)城市空間研究方向主要以真實(shí)案例分析為主,集中在發(fā)達(dá)城市,在空間結(jié)構(gòu)的定量化、模擬等方面進(jìn)步較為明顯,但在動(dòng)態(tài)模型與應(yīng)用、數(shù)理統(tǒng)計(jì)技術(shù)發(fā)展方面較為欠缺,因此加強(qiáng)動(dòng)態(tài)模型數(shù)理技術(shù)應(yīng)用,是現(xiàn)階段較為重要的發(fā)展方向。為實(shí)現(xiàn)城市空間增長(zhǎng)模擬預(yù)測(cè),將GIS空間分析技術(shù)構(gòu)造空間分析模型,對(duì)城市的擴(kuò)展速度、強(qiáng)度、方向和形態(tài)進(jìn)行分析研究,建立模型,預(yù)測(cè)空間形態(tài)演化趨勢(shì),優(yōu)化城市用地布局。

    1.城市空間增長(zhǎng)分析理論基礎(chǔ)

    1.1 GIS在城市空間增長(zhǎng)模型中的應(yīng)用

    地理信息系統(tǒng)(Geographic Information System,GIS)以計(jì)算機(jī)技術(shù)為基礎(chǔ),融合地理地圖學(xué)、測(cè)繪學(xué),是地球科學(xué)與信息科學(xué)的交叉性學(xué)科。它以計(jì)算機(jī)為硬件基礎(chǔ),利用地理空間數(shù)據(jù)庫(kù),對(duì)空間相關(guān)數(shù)據(jù)進(jìn)行采集、管理、操作、分析、模擬和顯示,并采用地理模型等分析方法,實(shí)時(shí)提供多種空間和動(dòng)態(tài)的地理信息,為地理研究和地理決策提供服務(wù),是現(xiàn)階段規(guī)劃管理、分析決策和研究地理空間數(shù)據(jù)的重要工具[4]。

    空間分析是以地理對(duì)象的時(shí)空位置信息為基礎(chǔ),結(jié)合其空間形態(tài)特征,實(shí)現(xiàn)地理對(duì)象的提取和空間信息傳輸??臻g分析功能是城市規(guī)劃應(yīng)用GIS技術(shù)的核心能力,在城市規(guī)劃應(yīng)用中發(fā)揮了巨大作用。GIS空間分析方法種類(lèi)多,不同的空間分析方法可通過(guò)不同方式定向定量地分析城市空間發(fā)展情況,探索城市擴(kuò)展特征,利用GIS可視化表達(dá)功能,詳細(xì)分析城市演變與更新過(guò)程(如表1所示):

    表1 空間分析類(lèi)型及作用

    1.2 元胞自動(dòng)機(jī)(CA)與SLEUTH模型

    元胞自動(dòng)機(jī)(CA)不是指單一的模型,而是一類(lèi)模型的總稱(chēng)。這類(lèi)模型是在C語(yǔ)言基礎(chǔ)上,結(jié)合特定的參數(shù),設(shè)定增長(zhǎng)規(guī)則,隨著時(shí)間的變化,可實(shí)現(xiàn)模擬數(shù)據(jù)不斷動(dòng)態(tài)更新。元胞自動(dòng)機(jī)模型較多且均滿(mǎn)足上述規(guī)律,適合應(yīng)用在城市空間擴(kuò)展模擬分析中。

    SLEUTH模型又稱(chēng)城市增長(zhǎng)元胞自動(dòng)機(jī)模型,包含了UGM(城市增長(zhǎng)模型)和LCD(土地利用模型)兩個(gè)子模型。該模型利用城市空間真實(shí)增長(zhǎng)現(xiàn)狀,獲取相關(guān)參數(shù),預(yù)測(cè)未來(lái)城市發(fā)展趨勢(shì)。具有六個(gè)輸入圖層:坡度(Slope)、交通(Transportat ion)、城市范圍(Urban)、土地利用(Land use)、排除層(Excluded)、山體陰影(Hillshade)。

    SLEUTH設(shè)計(jì)原理是將城市元胞化,并以此作為擴(kuò)散發(fā)起點(diǎn),在各種約束條件的加持下,分析計(jì)算各元胞發(fā)展可能性,反映城市趨勢(shì)。主要控制條件有:擴(kuò)散系數(shù)(Diffusion)、傳播系數(shù)(Spread)、繁殖系數(shù)(Breed)、道路引力系數(shù)(Road Gravity)和坡度阻礙系數(shù)(Slope)。在上述5個(gè)約束條件下,可以造成4種增長(zhǎng)類(lèi)型:

    (1)自發(fā)增長(zhǎng):表示會(huì)隨機(jī)產(chǎn)生城市化點(diǎn),即每個(gè)非城市化點(diǎn)都存在化為城市點(diǎn)的概率。

    (2)新擴(kuò)展中心增長(zhǎng):在不涉及外界干預(yù)的自然發(fā)展情況下,新建立的城市為城市擴(kuò)展中心。

    (3)邊界增長(zhǎng):從新擴(kuò)展中心增長(zhǎng)的城市擴(kuò)展中心邊界向外發(fā)展,邊界增長(zhǎng)是新擴(kuò)展中心增長(zhǎng)和中心增長(zhǎng)的總和。

    (4)道路影響型增長(zhǎng):在前面三個(gè)步驟中城市化的元胞與控制路邊產(chǎn)生的擴(kuò)展中心沿著交通線模擬新增長(zhǎng)的趨勢(shì)。

    1.3 SLEUTH預(yù)測(cè)模型的應(yīng)用

    SLEUTH模型的基本運(yùn)行單位是增長(zhǎng)環(huán)。首先通過(guò)多個(gè)增長(zhǎng)環(huán)實(shí)現(xiàn)預(yù)測(cè)分析,一般包括各項(xiàng)系數(shù)初始化,其次結(jié)合各類(lèi)增長(zhǎng)規(guī)則進(jìn)行數(shù)據(jù)模擬,最后生成模擬結(jié)果,如圖像數(shù)據(jù)和表格數(shù)據(jù)等[3]。

    SLEUTH模型包含了三個(gè)應(yīng)用模塊:測(cè)試模塊、校準(zhǔn)模塊和預(yù)測(cè)模塊。

    測(cè)試模塊,利用測(cè)試模塊分析相關(guān)數(shù)據(jù)可行性,測(cè)試該模型是否能夠使用。

    校準(zhǔn)模塊,包括粗校準(zhǔn)(Coarse Calibration)、精校準(zhǔn)(Fine Calibration)和終校準(zhǔn)(Final Calibration)三個(gè)階段,主要利用蒙特卡洛迭代方法,計(jì)算最佳控制系數(shù)。

    預(yù)測(cè)模型,在排除圖層中。不同類(lèi)型的用地,城鎮(zhèn)化阻力值不同。根據(jù)不同阻力值,利用校準(zhǔn)模型對(duì)預(yù)測(cè)模型進(jìn)行校準(zhǔn),將基年數(shù)據(jù)帶入校準(zhǔn)后的模型進(jìn)行預(yù)測(cè),可得到不同背景下的城市空間擴(kuò)展情況。具體流程(如圖1所示):

    圖1 SLEUTH模型執(zhí)行過(guò)程

    2.城市空間拓展分析

    城市發(fā)展過(guò)程,內(nèi)在表現(xiàn)為經(jīng)濟(jì)結(jié)構(gòu)和文化水平,外在表現(xiàn)為城市空間結(jié)構(gòu)分布,反映在城市各類(lèi)型用地的幾何形狀中。研究城市空間演變特征,預(yù)測(cè)城市空間拓展速度、方向,分析城市空間結(jié)構(gòu)分布,探索城市演變機(jī)制,可結(jié)合GIS技術(shù)中的空間分析法。利用GIS空間分析法,將表達(dá)城市擴(kuò)展程度的指數(shù)可視化,包括城市擴(kuò)展速度、強(qiáng)度、城市緊湊程度、形狀指數(shù)、分形維數(shù)指數(shù)等,利用這些指數(shù)可從時(shí)間維度、空間維度對(duì)城市空間擴(kuò)展進(jìn)行詳細(xì)分析。

    針對(duì)5組不同時(shí)間節(jié)點(diǎn)的遙感影像數(shù)據(jù),利用ArcGIS最大似然法進(jìn)行監(jiān)督分類(lèi),分析獲取這5組不同時(shí)間節(jié)點(diǎn)的土地利用情況,在此基礎(chǔ)上,利用重分類(lèi)工具對(duì)不同時(shí)期的土地利用圖進(jìn)行處理,結(jié)合GIS空間疊加分析功能進(jìn)行疊加處理,分析面積變化情況。

    2.1 擴(kuò)展速度分析

    擴(kuò)展速度是分析城市空間擴(kuò)展情況的基礎(chǔ)指標(biāo)。結(jié)合不同時(shí)期城市空間在固定時(shí)間范圍內(nèi)的面積變化,可計(jì)算出城市年均擴(kuò)展面積,具體計(jì)算方法如式(1)所示。為方便評(píng)價(jià)擴(kuò)展速度可進(jìn)行等級(jí)劃分(如表2所示),分析統(tǒng)計(jì)結(jié)果(如表3所示):

    表2 城市擴(kuò)展評(píng)價(jià)等級(jí)

    表3 擴(kuò)展分析數(shù)據(jù)統(tǒng)計(jì)

    式(1)中,V為該區(qū)域年均擴(kuò)展速度,單位為km2/a;Ua為該區(qū)域初期用地面積;Ub為該區(qū)域末期用地面積,單位為km2;Δt為初期與末期的時(shí)間差,單位為a。

    2.2 擴(kuò)展強(qiáng)度分析

    擴(kuò)展強(qiáng)度是評(píng)價(jià)城市不同時(shí)間段擴(kuò)散力度強(qiáng)弱的重要指標(biāo),指的是城區(qū)新增擴(kuò)展面積占總體土地面積的比值,擴(kuò)展強(qiáng)度指標(biāo)越高,表示擴(kuò)展力度越大,擴(kuò)展強(qiáng)度指標(biāo)越低,表示擴(kuò)展力度越小,具體計(jì)算式如式(2)所示:

    式(2)中,EⅡ?yàn)閿U(kuò)散強(qiáng)度指數(shù);U為該區(qū)域總體面積,單位為km2;Ua為該區(qū)域初期用地面積,單位為km2;Ub為該區(qū)域末期用地面積,單位為km2;Δt為初期與末期的時(shí)間差,單位為a。

    2.3 外部形態(tài)緊湊度

    城市外圍輪廓形狀緊湊度反映了城市空間形狀的緊湊程度,可通過(guò)計(jì)算城市面積與周長(zhǎng)的關(guān)系,分析城市擴(kuò)展的周期性[7]。緊湊度計(jì)算方法如式(3)所示:

    式(3)中,BCI為城市用地邊界緊湊度;A為用地面積,單位為km2;P為城市邊界周長(zhǎng)。由式(3)可知:BCI取值范圍在0-1之間,越接近1,緊湊性程度越高,形狀越趨于圓形;越接近0,則緊湊程度越低。

    2.4 擴(kuò)展方向分析

    城市擴(kuò)展方向分析通常采用象限方位分析法,該方法通常以原始城區(qū)為原點(diǎn),以城市擴(kuò)展最遠(yuǎn)點(diǎn)的連線為半徑,構(gòu)建圓形,以正北方向?yàn)槠瘘c(diǎn),進(jìn)行象限劃分,象限可劃分8個(gè)區(qū)域,結(jié)合GIS疊加分析得到各區(qū)域的面積,繪制不同時(shí)期城市擴(kuò)展方向象限圖(如圖2所示):

    圖2 擴(kuò)展方向分析圖

    由表3可知:在1996-2000、2001-2005兩個(gè)階段,該研究區(qū)域擴(kuò)展速度為較快發(fā)展階段,相應(yīng)的擴(kuò)展強(qiáng)度也屬于中等強(qiáng)度;在2006-2010、2011-2015兩個(gè)階段,擴(kuò)展速度較快,呈現(xiàn)高速發(fā)展現(xiàn)象,同時(shí)期擴(kuò)展強(qiáng)度也呈現(xiàn)高強(qiáng)度發(fā)展現(xiàn)象。結(jié)合圖2擴(kuò)展方向分析,隨著擴(kuò)展強(qiáng)度和擴(kuò)展速度地不斷加快,緊湊度不斷變小,發(fā)展方向出現(xiàn)不均衡現(xiàn)象。

    3.SLEUTH 模型

    3.1 建立數(shù)據(jù)集

    模型的預(yù)測(cè)過(guò)程是在歷史數(shù)據(jù)的基礎(chǔ)上進(jìn)行的,但隨著數(shù)據(jù)不斷增加,數(shù)據(jù)分辨率不斷提高,預(yù)測(cè)時(shí)間也會(huì)出現(xiàn)指數(shù)增長(zhǎng)現(xiàn)象,同時(shí)數(shù)據(jù)分辨率對(duì)預(yù)測(cè)結(jié)構(gòu)也有著重要影響。由文獻(xiàn)可知,將分辨率為60m的原始數(shù)據(jù),用于預(yù)測(cè)模型基礎(chǔ)數(shù)據(jù)進(jìn)行重采樣,在預(yù)測(cè)精度和耗費(fèi)時(shí)間上都有著比較好的效果。

    SLEUTH模型建立的基礎(chǔ)是原始城區(qū)各類(lèi)型數(shù)據(jù),原始數(shù)據(jù)集的建立對(duì)未來(lái)城市空間增長(zhǎng)邊界的預(yù)測(cè)起決定性作用。收集多個(gè)不同年份的已有數(shù)據(jù)集主要包括衛(wèi)星影像數(shù)據(jù)、城區(qū)總體規(guī)劃圖、電子地圖、DEM數(shù)據(jù)等。

    針對(duì)原始衛(wèi)星影像數(shù)據(jù)進(jìn)行配準(zhǔn)、重采樣與重分類(lèi),以便于準(zhǔn)確提取城市化圖層。針對(duì)總體規(guī)劃圖,利用ArcGIS軟件,結(jié)合實(shí)地調(diào)查數(shù)據(jù)進(jìn)行幾何配準(zhǔn),分別提取出各個(gè)年份交通道路圖,并進(jìn)行柵格處理。利用ArcGIS軟件對(duì)DEM數(shù)據(jù)進(jìn)行重采樣處理。

    3.2 SLEUTH模型構(gòu)建

    輸入數(shù)據(jù)集利用ArcGIS軟件,截取DEM數(shù)字高程圖整數(shù)值建立坡度層、山體陰影層。利用ArcGIS軟件將城區(qū)矢量數(shù)據(jù)轉(zhuǎn)換為分辨率60m的柵格數(shù)據(jù),并進(jìn)行二進(jìn)制重分類(lèi)處理,構(gòu)建城市化層,并以此作為控制層進(jìn)行校準(zhǔn)模型。為提取道路權(quán)重分配的基礎(chǔ)圖層,矢量化處理遙感影像,提取各時(shí)期道路矢量數(shù)據(jù),針對(duì)不同等級(jí)道路賦予不同等級(jí)權(quán)重值,賦值權(quán)重從大到小,依次為高速公路=4、主干路=2、次干路=1、非道路=0。城市空間擴(kuò)展的限制因素主要通過(guò)排除層進(jìn)行限制,利用排除層設(shè)置城市空間擴(kuò)展的允許或限制區(qū)域,同時(shí)根據(jù)對(duì)不同區(qū)域賦值的大小來(lái)限制擴(kuò)展速度、強(qiáng)度以及擴(kuò)展的先后順序。設(shè)置排除層,考慮固有限制條件,如,基本耕地、水域、生態(tài)保護(hù)區(qū)、歷史文化保護(hù)區(qū)等,將此類(lèi)限制因素進(jìn)行賦值,調(diào)整模型,限制城市增長(zhǎng)速度及方向[5]。

    3.3 模型校準(zhǔn)

    在利用模型預(yù)測(cè)前,需對(duì)模型進(jìn)行參數(shù)校準(zhǔn),在對(duì)SLEUTH模型進(jìn)行校準(zhǔn)過(guò)程中,會(huì)產(chǎn)生13個(gè)參數(shù),使用OSM(the Optimal SLEUTH Metric)算法評(píng)估模型校準(zhǔn)結(jié)果,得到最佳控制系數(shù)組合[6]。具體模型校準(zhǔn)過(guò)程如下:

    (1)粗校準(zhǔn)階段

    初始校準(zhǔn)階段,數(shù)據(jù)空間分辨率120m,模擬系數(shù)取值范圍設(shè)定為0-100,取值步長(zhǎng)為25,進(jìn)行10次蒙特卡洛迭代,以O(shè)SM為關(guān)鍵值進(jìn)行降序排列,排列順序?yàn)?-5的值作為取值范圍,產(chǎn)生5個(gè)新的系數(shù)區(qū)段。

    (2)精校準(zhǔn)階段

    精校準(zhǔn)數(shù)據(jù)空間分辨為120m,利用粗校準(zhǔn)階段生成的新的區(qū)段,取值步長(zhǎng)為5再進(jìn)行多次蒙特卡洛迭代。再次重復(fù)選擇5個(gè)新區(qū)段。

    (3)終校準(zhǔn)階段

    終校準(zhǔn)階段空間分辨為60m,利用精校準(zhǔn)階段生成的新的區(qū)段,取步長(zhǎng)為2再進(jìn)行多次蒙特卡洛迭代,取OSM排名第一的作為最優(yōu)控制系數(shù)的值。

    3.4 多情景城市增長(zhǎng)預(yù)測(cè)

    結(jié)合最佳控制系數(shù),針對(duì)不同情景環(huán)境分析城市增長(zhǎng)情況,考慮城市規(guī)劃、現(xiàn)行耕地保護(hù)政策以及生態(tài)環(huán)境,設(shè)計(jì)4種背景,進(jìn)行預(yù)測(cè)分析。

    背景1:排除層設(shè)置,首先考慮已規(guī)劃區(qū)域的各種用地類(lèi)型,包括居住用地、城市中心用地、城鎮(zhèn)產(chǎn)業(yè)發(fā)展用地,統(tǒng)一將此類(lèi)用地灰度賦值為0,未規(guī)劃區(qū)域灰度賦值為80,水體等灰度賦值為100。

    背景2:考慮土地利用分類(lèi)情況,結(jié)合耕地保護(hù)政策,排除層設(shè)置為:在道路1km范圍內(nèi)耕地灰度賦值為80,超1km范圍耕地灰度賦值為50,水體賦值為100。

    背景3:考慮生態(tài)保護(hù)要求,將生態(tài)保護(hù)的區(qū)域和水域賦值為100,其余區(qū)域賦值為0。

    以60m為預(yù)測(cè)分辨率,利用SLEUTH模型的預(yù)測(cè)模塊,設(shè)定蒙特卡洛迭代次數(shù)為100,預(yù)測(cè)2015-2030年的城市空間擴(kuò)展情況。結(jié)果(如表4所示):

    表4 不同情景下2030年模擬用地分類(lèi)與比重

    由表4可知:到2030年,在城市規(guī)劃的影響下,可被城市化的高概率用地占比為3.02%,可被城市化的低概率用地占比5.63%,在耕地保護(hù)政策的影響下,可被城市化的高概率用地占比為8.65%,可被城市化的低概率用地占比4.63%;在生態(tài)環(huán)境保護(hù)的影響下,可被城市化的高概率用地占比為7.26%,可被城市化的低概率用地占比3.15%。在城市規(guī)劃背景的,三個(gè)不同限制條件影響下,可城市化的像元最少約為38152,比重約為8.65%。

    4.結(jié)束語(yǔ)

    為規(guī)劃城市空間結(jié)構(gòu),促進(jìn)城市合理發(fā)展,在GIS空間技術(shù)的基礎(chǔ)上,研究應(yīng)用了SLEUTH模型對(duì)城市空間擴(kuò)展進(jìn)行預(yù)測(cè),具體工作內(nèi)容如下:

    (1)論述了GIS空間分析技術(shù)應(yīng)用細(xì)節(jié)及范圍,針對(duì)SLEUTH模型,從其基本算法原理、數(shù)據(jù)集、控制條件、擴(kuò)散形式以及使用細(xì)節(jié),進(jìn)行了分析研究,為實(shí)現(xiàn)擴(kuò)展預(yù)測(cè)提供了技術(shù)原理。

    (2)針對(duì)歷史數(shù)據(jù)進(jìn)行了擴(kuò)散情況分析,包括城市擴(kuò)展速度、強(qiáng)度、城市緊湊度,并利用象限圖進(jìn)行了擴(kuò)展方向的分析。利用多種原始數(shù)據(jù)集,進(jìn)行了預(yù)測(cè)模型的構(gòu)建,通過(guò)多次模型校準(zhǔn),進(jìn)行了最佳參數(shù)的選取,并以此參數(shù)實(shí)現(xiàn)了對(duì)不同發(fā)展背景下城市擴(kuò)展的預(yù)測(cè)。

    猜你喜歡
    城市化用地預(yù)測(cè)
    無(wú)可預(yù)測(cè)
    黃河之聲(2022年10期)2022-09-27 13:59:46
    選修2-2期中考試預(yù)測(cè)卷(A卷)
    選修2-2期中考試預(yù)測(cè)卷(B卷)
    城鄉(xiāng)建設(shè)用地增減掛鉤研究進(jìn)展綜述
    不必預(yù)測(cè)未來(lái),只需把握現(xiàn)在
    城鄉(xiāng)建設(shè)用地增減掛鉤政策的演變
    城鄉(xiāng)建設(shè)用地增減掛鉤的實(shí)踐與認(rèn)識(shí)
    城鄉(xiāng)建設(shè)用地增減掛鉤的實(shí)踐與認(rèn)識(shí)
    失衡的城市化:現(xiàn)狀與出路
    湖湘論壇(2015年4期)2015-12-01 09:30:08
    “城市化”諸概念辨析
    亚洲欧美日韩东京热| 国产亚洲精品av在线| 午夜福利在线观看吧| 久久国产乱子伦精品免费另类| 熟妇人妻久久中文字幕3abv| 男人的好看免费观看在线视频| 91av网一区二区| 欧美日韩亚洲国产一区二区在线观看| 美女高潮喷水抽搐中文字幕| 国产精品野战在线观看| 90打野战视频偷拍视频| av在线蜜桃| 最新美女视频免费是黄的| 天堂√8在线中文| 国产老妇女一区| 成人亚洲精品av一区二区| 精品乱码久久久久久99久播| 老熟妇乱子伦视频在线观看| 人人妻人人看人人澡| 久久草成人影院| 国产中年淑女户外野战色| 亚洲成av人片免费观看| 极品教师在线免费播放| 此物有八面人人有两片| 真人做人爱边吃奶动态| 别揉我奶头~嗯~啊~动态视频| 亚洲aⅴ乱码一区二区在线播放| 伊人久久大香线蕉亚洲五| 在线天堂最新版资源| av专区在线播放| 欧美丝袜亚洲另类 | 午夜视频国产福利| 9191精品国产免费久久| 国产精品美女特级片免费视频播放器| 亚洲av美国av| 亚洲 国产 在线| 亚洲欧美日韩高清在线视频| 国产精华一区二区三区| 久久久久久久久久黄片| 伊人久久大香线蕉亚洲五| 日韩欧美国产在线观看| 日本黄色片子视频| 在线观看午夜福利视频| 久久久久久久午夜电影| 女生性感内裤真人,穿戴方法视频| 99久久久亚洲精品蜜臀av| 99热这里只有是精品50| 很黄的视频免费| 日韩欧美精品免费久久 | 久久久久久久久中文| 老司机福利观看| 99riav亚洲国产免费| 亚洲内射少妇av| 91麻豆精品激情在线观看国产| 老鸭窝网址在线观看| 国产精品一及| 一个人看的www免费观看视频| 欧美黄色淫秽网站| 久久久久久久精品吃奶| 宅男免费午夜| 久久欧美精品欧美久久欧美| 国产99白浆流出| 夜夜躁狠狠躁天天躁| 桃色一区二区三区在线观看| 久久性视频一级片| 乱人视频在线观看| 丁香欧美五月| 国产99白浆流出| 在线十欧美十亚洲十日本专区| 国产综合懂色| av片东京热男人的天堂| 在线观看日韩欧美| 国产高清视频在线观看网站| 国产不卡一卡二| 久久久久久人人人人人| 免费av观看视频| 国产av一区在线观看免费| 少妇人妻精品综合一区二区 | 变态另类成人亚洲欧美熟女| 狂野欧美激情性xxxx| 午夜福利在线观看吧| 色吧在线观看| 国产精品久久久久久亚洲av鲁大| 少妇高潮的动态图| 啦啦啦观看免费观看视频高清| 国产精品永久免费网站| 亚洲国产色片| av天堂中文字幕网| 成人一区二区视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品久久国产高清桃花| 九色成人免费人妻av| 色噜噜av男人的天堂激情| 白带黄色成豆腐渣| 欧美成人免费av一区二区三区| 亚洲熟妇中文字幕五十中出| 亚洲18禁久久av| 伊人久久大香线蕉亚洲五| 波多野结衣巨乳人妻| 国产三级黄色录像| 一进一出抽搐gif免费好疼| netflix在线观看网站| 中出人妻视频一区二区| 久久久久久久精品吃奶| 国产乱人视频| 好看av亚洲va欧美ⅴa在| 天堂影院成人在线观看| a在线观看视频网站| 亚洲熟妇熟女久久| 欧美极品一区二区三区四区| 日本免费一区二区三区高清不卡| 脱女人内裤的视频| 亚洲精品456在线播放app | 女警被强在线播放| 欧美日韩乱码在线| 日本精品一区二区三区蜜桃| 日韩欧美精品v在线| 黄色日韩在线| 天堂av国产一区二区熟女人妻| 午夜激情欧美在线| 一本久久中文字幕| 亚洲精品粉嫩美女一区| 村上凉子中文字幕在线| 日本黄色视频三级网站网址| av专区在线播放| 美女高潮的动态| 欧美乱妇无乱码| 国产精品精品国产色婷婷| 九九久久精品国产亚洲av麻豆| 成人av一区二区三区在线看| 老司机深夜福利视频在线观看| 日本黄色片子视频| 嫁个100分男人电影在线观看| 岛国视频午夜一区免费看| 日本黄色视频三级网站网址| 91在线观看av| 久久人妻av系列| 欧美黄色淫秽网站| 欧美日韩福利视频一区二区| 91在线精品国自产拍蜜月 | 长腿黑丝高跟| 亚洲欧美日韩无卡精品| 欧美三级亚洲精品| 韩国av一区二区三区四区| 男女那种视频在线观看| xxxwww97欧美| 国产精华一区二区三区| 听说在线观看完整版免费高清| 色噜噜av男人的天堂激情| 网址你懂的国产日韩在线| 最近最新中文字幕大全电影3| 亚洲人与动物交配视频| 1024手机看黄色片| 国产精品久久久人人做人人爽| 国产三级中文精品| 又紧又爽又黄一区二区| 天堂av国产一区二区熟女人妻| 欧美中文综合在线视频| 99国产精品一区二区蜜桃av| АⅤ资源中文在线天堂| 在线观看一区二区三区| 99久久久亚洲精品蜜臀av| 色哟哟哟哟哟哟| 无限看片的www在线观看| 亚洲不卡免费看| 欧美3d第一页| 国产精品美女特级片免费视频播放器| 18美女黄网站色大片免费观看| 欧美日韩中文字幕国产精品一区二区三区| 日本 av在线| 91久久精品电影网| 一个人看视频在线观看www免费 | 欧美激情在线99| 久久国产精品影院| 欧美乱色亚洲激情| 一本综合久久免费| 久久人妻av系列| 51国产日韩欧美| 日日干狠狠操夜夜爽| 国产真人三级小视频在线观看| 99久久99久久久精品蜜桃| 啦啦啦韩国在线观看视频| 日日夜夜操网爽| 国产成人欧美在线观看| 国产中年淑女户外野战色| 久久性视频一级片| 岛国在线免费视频观看| 亚洲内射少妇av| 欧美黑人欧美精品刺激| 亚洲美女视频黄频| 国产欧美日韩一区二区精品| 亚洲欧美日韩高清专用| 午夜精品在线福利| 天天躁日日操中文字幕| 欧美日本亚洲视频在线播放| 欧美高清成人免费视频www| 脱女人内裤的视频| 高清在线国产一区| 精品熟女少妇八av免费久了| 嫩草影院精品99| 久久精品夜夜夜夜夜久久蜜豆| 91av网一区二区| 日韩欧美 国产精品| 天堂av国产一区二区熟女人妻| 久久精品国产亚洲av涩爱 | 国产av麻豆久久久久久久| 校园春色视频在线观看| 91麻豆av在线| 美女 人体艺术 gogo| 久久久国产成人免费| 久99久视频精品免费| 成人无遮挡网站| 又黄又爽又免费观看的视频| 国产伦在线观看视频一区| 国产单亲对白刺激| 看片在线看免费视频| 婷婷精品国产亚洲av| 免费av不卡在线播放| 欧美乱码精品一区二区三区| 国产伦精品一区二区三区四那| 国产在视频线在精品| 尤物成人国产欧美一区二区三区| 亚洲不卡免费看| 欧美极品一区二区三区四区| 亚洲欧美激情综合另类| 一本一本综合久久| 天天躁日日操中文字幕| 亚洲aⅴ乱码一区二区在线播放| 久久6这里有精品| 亚洲国产精品成人综合色| 亚洲黑人精品在线| 日韩高清综合在线| 男插女下体视频免费在线播放| 一本久久中文字幕| 亚洲精品影视一区二区三区av| а√天堂www在线а√下载| 亚洲五月婷婷丁香| 欧美在线一区亚洲| 久久天躁狠狠躁夜夜2o2o| 中国美女看黄片| 麻豆成人午夜福利视频| 美女 人体艺术 gogo| 操出白浆在线播放| 老司机午夜十八禁免费视频| 热99re8久久精品国产| 在线看三级毛片| 在线观看日韩欧美| 99精品久久久久人妻精品| 久久精品91无色码中文字幕| 一夜夜www| 亚洲国产精品成人综合色| 欧美激情在线99| 亚洲aⅴ乱码一区二区在线播放| 久久精品91无色码中文字幕| 欧美成人a在线观看| 两个人视频免费观看高清| 黄片大片在线免费观看| 在线观看午夜福利视频| 亚洲av免费高清在线观看| 非洲黑人性xxxx精品又粗又长| 成人鲁丝片一二三区免费| 亚洲av电影不卡..在线观看| 99热6这里只有精品| 别揉我奶头~嗯~啊~动态视频| 别揉我奶头~嗯~啊~动态视频| bbb黄色大片| 免费在线观看亚洲国产| 亚洲一区高清亚洲精品| 老司机午夜福利在线观看视频| 99久国产av精品| 夜夜躁狠狠躁天天躁| 精品国内亚洲2022精品成人| 欧美一级a爱片免费观看看| 美女被艹到高潮喷水动态| 尤物成人国产欧美一区二区三区| 欧美丝袜亚洲另类 | 日韩有码中文字幕| 美女 人体艺术 gogo| 国产视频内射| 欧美丝袜亚洲另类 | 不卡一级毛片| 欧美激情久久久久久爽电影| 欧美日韩精品网址| 一级黄色大片毛片| 国产一区二区三区视频了| 日韩免费av在线播放| 两个人的视频大全免费| 久久久久久九九精品二区国产| 国内少妇人妻偷人精品xxx网站| 日韩欧美 国产精品| h日本视频在线播放| 麻豆成人av在线观看| 亚洲人成网站在线播放欧美日韩| 午夜福利在线观看免费完整高清在 | 成年女人永久免费观看视频| 淫秽高清视频在线观看| 岛国在线免费视频观看| 亚洲性夜色夜夜综合| 国产精品精品国产色婷婷| 日本免费a在线| 日韩av在线大香蕉| 亚洲黑人精品在线| 午夜老司机福利剧场| 国产淫片久久久久久久久 | 国产激情偷乱视频一区二区| 亚洲最大成人中文| 免费看a级黄色片| 岛国在线免费视频观看| 亚洲aⅴ乱码一区二区在线播放| 国产亚洲av嫩草精品影院| 欧美大码av| 三级男女做爰猛烈吃奶摸视频| 草草在线视频免费看| 免费搜索国产男女视频| 日本三级黄在线观看| 国产探花极品一区二区| 欧美三级亚洲精品| 国产真实伦视频高清在线观看 | АⅤ资源中文在线天堂| 亚洲精品粉嫩美女一区| 欧美性感艳星| 国产精品 国内视频| 欧美色视频一区免费| 国产亚洲欧美98| 少妇的丰满在线观看| 午夜福利成人在线免费观看| 亚洲天堂国产精品一区在线| 国产精品亚洲av一区麻豆| 美女高潮的动态| 国产 一区 欧美 日韩| 国产激情欧美一区二区| 国产成人a区在线观看| 亚洲在线观看片| 在线观看舔阴道视频| 精品福利观看| 成年女人看的毛片在线观看| 91字幕亚洲| 国产真实伦视频高清在线观看 | 一进一出好大好爽视频| 十八禁网站免费在线| 国产伦人伦偷精品视频| a在线观看视频网站| 国产中年淑女户外野战色| 法律面前人人平等表现在哪些方面| 免费在线观看日本一区| 国产亚洲精品一区二区www| 婷婷精品国产亚洲av| av在线蜜桃| 亚洲精品一区av在线观看| 久久亚洲精品不卡| 欧美黑人欧美精品刺激| 欧美日韩乱码在线| 一区二区三区免费毛片| 午夜精品久久久久久毛片777| 人妻久久中文字幕网| 精品久久久久久久久久免费视频| 国产精品99久久久久久久久| 全区人妻精品视频| 女人高潮潮喷娇喘18禁视频| 色吧在线观看| 日本免费a在线| 亚洲av免费高清在线观看| 免费在线观看亚洲国产| 亚洲人成网站在线播| 一个人免费在线观看电影| 男插女下体视频免费在线播放| 色视频www国产| 国产精品 欧美亚洲| 高清在线国产一区| 黑人欧美特级aaaaaa片| 欧美日本视频| 国产乱人伦免费视频| 在线观看av片永久免费下载| 男女做爰动态图高潮gif福利片| 天天添夜夜摸| 午夜激情欧美在线| 看免费av毛片| 免费av观看视频| 国产成人av激情在线播放| 亚洲在线观看片| 亚洲第一电影网av| 一本综合久久免费| 国产极品精品免费视频能看的| 少妇人妻一区二区三区视频| 香蕉av资源在线| 日本a在线网址| 黄色片一级片一级黄色片| 国内精品久久久久久久电影| 国产高清videossex| 在线播放无遮挡| 老汉色∧v一级毛片| 国产爱豆传媒在线观看| 草草在线视频免费看| 俄罗斯特黄特色一大片| 嫁个100分男人电影在线观看| 国产v大片淫在线免费观看| 国产探花在线观看一区二区| 在线a可以看的网站| 亚洲精品亚洲一区二区| 久久久久久久久大av| а√天堂www在线а√下载| 在线观看免费午夜福利视频| 国产v大片淫在线免费观看| 国产伦精品一区二区三区四那| 一区二区三区激情视频| 男人舔女人下体高潮全视频| 乱人视频在线观看| 久久久国产成人免费| 亚洲天堂国产精品一区在线| 色老头精品视频在线观看| 午夜福利在线在线| av欧美777| 国产精品精品国产色婷婷| 欧美乱妇无乱码| 国内精品久久久久久久电影| 国内精品久久久久精免费| 国产黄片美女视频| 午夜影院日韩av| 精品久久久久久久毛片微露脸| 精品欧美国产一区二区三| 久久伊人香网站| 日韩国内少妇激情av| 久久精品人妻少妇| 欧美黑人欧美精品刺激| 国产一区二区三区在线臀色熟女| 真人一进一出gif抽搐免费| 首页视频小说图片口味搜索| 99精品在免费线老司机午夜| 亚洲精品在线美女| 一进一出好大好爽视频| 午夜福利在线观看吧| 亚洲人成网站在线播| 成人精品一区二区免费| 国产免费一级a男人的天堂| 精品国产三级普通话版| 日日干狠狠操夜夜爽| 99久久精品国产亚洲精品| 老司机午夜福利在线观看视频| 18+在线观看网站| 此物有八面人人有两片| 天美传媒精品一区二区| 亚洲男人的天堂狠狠| 人妻夜夜爽99麻豆av| 我要搜黄色片| 国产高清videossex| 中文字幕人成人乱码亚洲影| 久久久久性生活片| 精品一区二区三区视频在线观看免费| 91久久精品国产一区二区成人 | 人人妻人人看人人澡| 精华霜和精华液先用哪个| 久久精品91无色码中文字幕| 精品久久久久久久久久免费视频| 嫩草影视91久久| 亚洲国产精品成人综合色| 制服人妻中文乱码| 国产老妇女一区| 美女黄网站色视频| 日本一二三区视频观看| 久久久精品大字幕| 久久久久久人人人人人| 亚洲国产欧美人成| 亚洲五月婷婷丁香| 欧美黄色淫秽网站| 午夜精品一区二区三区免费看| 热99re8久久精品国产| 欧美大码av| 成人高潮视频无遮挡免费网站| 亚洲精品一区av在线观看| 深夜精品福利| 嫩草影院入口| 男插女下体视频免费在线播放| 国产av麻豆久久久久久久| 婷婷精品国产亚洲av在线| 99在线人妻在线中文字幕| 色视频www国产| 在线天堂最新版资源| 18禁在线播放成人免费| 人妻丰满熟妇av一区二区三区| 国产真实伦视频高清在线观看 | tocl精华| 天天一区二区日本电影三级| 午夜福利欧美成人| 男插女下体视频免费在线播放| 观看免费一级毛片| 国产男靠女视频免费网站| 国产成人影院久久av| 在线播放无遮挡| eeuss影院久久| 成人性生交大片免费视频hd| 9191精品国产免费久久| 亚洲精品成人久久久久久| 麻豆国产av国片精品| 欧美精品啪啪一区二区三区| 99国产精品一区二区三区| 在线a可以看的网站| 国产成人福利小说| 国产淫片久久久久久久久 | 国产一级毛片七仙女欲春2| 99在线人妻在线中文字幕| 别揉我奶头~嗯~啊~动态视频| 国产熟女xx| 午夜免费成人在线视频| 亚洲精品久久国产高清桃花| 国产成+人综合+亚洲专区| 深夜精品福利| 亚洲人成网站在线播| 人人妻人人澡欧美一区二区| 美女cb高潮喷水在线观看| 国产成+人综合+亚洲专区| 色吧在线观看| 很黄的视频免费| 精品福利观看| 精华霜和精华液先用哪个| 国产精品日韩av在线免费观看| 黄色视频,在线免费观看| 99久久成人亚洲精品观看| 国产乱人视频| 免费搜索国产男女视频| x7x7x7水蜜桃| 国产伦在线观看视频一区| 男人舔女人下体高潮全视频| 99在线视频只有这里精品首页| 日韩欧美在线乱码| 国产野战对白在线观看| 色播亚洲综合网| 在线观看日韩欧美| 国产探花在线观看一区二区| 91av网一区二区| 一本久久中文字幕| 在线观看舔阴道视频| 亚洲成人精品中文字幕电影| 亚洲黑人精品在线| svipshipincom国产片| 看免费av毛片| 国产av在哪里看| 亚洲av免费在线观看| 成年免费大片在线观看| 夜夜躁狠狠躁天天躁| 每晚都被弄得嗷嗷叫到高潮| 最近视频中文字幕2019在线8| svipshipincom国产片| 精品人妻1区二区| 麻豆成人午夜福利视频| 国产亚洲精品av在线| 色吧在线观看| 麻豆国产97在线/欧美| 欧美日韩亚洲国产一区二区在线观看| 免费看a级黄色片| 别揉我奶头~嗯~啊~动态视频| 动漫黄色视频在线观看| 免费看美女性在线毛片视频| 两个人看的免费小视频| 五月伊人婷婷丁香| 欧美色视频一区免费| 午夜久久久久精精品| 在线观看av片永久免费下载| 男女视频在线观看网站免费| 欧美一级毛片孕妇| 悠悠久久av| 国产精品嫩草影院av在线观看 | 久久精品综合一区二区三区| 欧美成人免费av一区二区三区| 伊人久久大香线蕉亚洲五| 一本一本综合久久| 亚洲av二区三区四区| 一个人免费在线观看的高清视频| 夜夜躁狠狠躁天天躁| 国产蜜桃级精品一区二区三区| 精品熟女少妇八av免费久了| 精品福利观看| 日韩欧美免费精品| 国产免费男女视频| 日本a在线网址| 脱女人内裤的视频| 欧美乱码精品一区二区三区| 亚洲av成人av| 一本精品99久久精品77| 一区福利在线观看| 日韩欧美国产在线观看| 九九久久精品国产亚洲av麻豆| 亚洲天堂国产精品一区在线| 国产91精品成人一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 国产蜜桃级精品一区二区三区| 成人18禁在线播放| 久久亚洲精品不卡| 99久久成人亚洲精品观看| 久久99热这里只有精品18| 蜜桃久久精品国产亚洲av| av欧美777| 日韩欧美在线二视频| 亚洲美女视频黄频| 日本在线视频免费播放| 久久婷婷人人爽人人干人人爱| 法律面前人人平等表现在哪些方面| 午夜福利欧美成人| 国产v大片淫在线免费观看| 色在线成人网| 免费在线观看亚洲国产| 亚洲乱码一区二区免费版| 色av中文字幕| 国产精品av视频在线免费观看| 国产aⅴ精品一区二区三区波| 男女床上黄色一级片免费看| 观看免费一级毛片| 男女床上黄色一级片免费看| 天堂网av新在线| 国产成人福利小说| 麻豆一二三区av精品| 90打野战视频偷拍视频| 国产亚洲精品综合一区在线观看| 舔av片在线| 欧美日韩国产亚洲二区| 叶爱在线成人免费视频播放|