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

    基于EnKF-MCRP模型的生態(tài)用地?cái)U(kuò)張模擬研究

    2016-10-27 02:04:09岳德鵬DiYang張啟斌李宇彤
    關(guān)鍵詞:源地檢測點(diǎn)阻力

    于 強(qiáng) 岳德鵬 Di Yang 張啟斌 馬 歡 李宇彤

    (1.北京林業(yè)大學(xué)精準(zhǔn)林業(yè)北京市重點(diǎn)實(shí)驗(yàn)室, 北京 100083; 2.佛羅里達(dá)大學(xué)地理系, 蓋恩斯維爾 32611;3.國土資源部土地整治中心, 北京 100035)

    ?

    基于EnKF-MCRP模型的生態(tài)用地?cái)U(kuò)張模擬研究

    于強(qiáng)1岳德鵬1Di Yang2張啟斌1馬歡1李宇彤3

    (1.北京林業(yè)大學(xué)精準(zhǔn)林業(yè)北京市重點(diǎn)實(shí)驗(yàn)室, 北京 100083; 2.佛羅里達(dá)大學(xué)地理系, 蓋恩斯維爾 32611;3.國土資源部土地整治中心, 北京 100035)

    以生態(tài)脆弱區(qū)典型縣域磴口縣為研究區(qū),基于2002、2007、2012、2015年4期遙感影像解譯數(shù)據(jù),將數(shù)據(jù)同化技術(shù)引入生態(tài)源地的變化模擬中,考慮生態(tài)障礙和生態(tài)阻力構(gòu)建EnKF-MCRP模型,進(jìn)行磴口縣生態(tài)用地的變化模擬。結(jié)果表明,引入數(shù)據(jù)同化技術(shù)的EnKF-CA/Markov模型的模擬總精度達(dá)到82.4%,數(shù)據(jù)同化能夠減少誤差的積累。根據(jù)擴(kuò)展能力磴口縣生態(tài)源地共分為5個(gè)等級(jí),其中3、4、5級(jí)生態(tài)源地的空間布局形成東北-西南、西北-西南的沙漠化防護(hù)格局。引入生態(tài)源地變化的EnKF-MCRP模型的生態(tài)用地?cái)U(kuò)張模擬精度最高,生態(tài)用地面積與空間布局最接近實(shí)際情況,方差達(dá)到0.4。此研究可為當(dāng)前以及未來的生態(tài)用地規(guī)劃和管理提供科學(xué)根據(jù)。

    生態(tài)用地; 數(shù)據(jù)同化; EnKF-MCRP模型; 磴口縣

    引言

    生態(tài)用地是具有重要的生態(tài)系統(tǒng)服務(wù)功能或生態(tài)敏感性較高、生態(tài)環(huán)境較為脆弱、對(duì)維護(hù)區(qū)域關(guān)鍵生態(tài)過程發(fā)揮重要作用的土地單元[1]。在西北荒漠綠洲過渡區(qū),覆蓋區(qū)域的生態(tài)系統(tǒng)極其脆弱,生態(tài)用地作為干旱區(qū)防治和減緩?fù)恋鼗哪铀贁U(kuò)展的“緩沖劑”[2],還具有維護(hù)綠洲生態(tài)系統(tǒng)健康和持續(xù)發(fā)展、提供生態(tài)服務(wù)等多重功能,因此荒漠綠洲生態(tài)脆弱區(qū)生態(tài)用地變化模擬已成為全世界學(xué)者們的研究熱點(diǎn)。生態(tài)用地的重建、模擬及預(yù)測可以為荒漠綠洲區(qū)生態(tài)用地資源的合理利用、科學(xué)規(guī)劃以及改善區(qū)域生態(tài)環(huán)境提供有效示范和支撐。

    最小累計(jì)阻力模型(Minimal cumulative resistance,MCR)被廣泛應(yīng)用于物種保護(hù)和景觀格局分析等生態(tài)領(lǐng)域[3-4],眾多學(xué)者利用MCR模型進(jìn)行了旅游用地規(guī)劃、城市用地?cái)U(kuò)張模擬、生態(tài)保護(hù)用地?cái)U(kuò)張趨勢分析、城市土地生態(tài)適宜性評(píng)價(jià)、山區(qū)城市增長邊界劃定等研究[5-6],在構(gòu)建生態(tài)阻力模型時(shí)均考慮了源、距離和阻力3個(gè)因素,目前考慮源地等級(jí)研究較少,考慮源地的變化則基本沒有[7]。作為源地的生態(tài)用地意義重大,其隨著時(shí)間及外界條件不斷發(fā)生變化,數(shù)據(jù)同化能夠同化不同來源的觀測數(shù)據(jù)以修正動(dòng)態(tài)過程模型,目前眾多學(xué)者應(yīng)用不同的數(shù)據(jù)同化算法進(jìn)行了城市擴(kuò)張模擬[8]、土地利用模擬[9]、陸面系統(tǒng)模擬[10]等研究,研究表明數(shù)據(jù)同化能更新(修正)模型并提高模型的模擬精度[8]。本文將目前的研究熱點(diǎn)數(shù)據(jù)同化技術(shù)引進(jìn)到生態(tài)源地的變化模擬中[11],利用同化值修正源地變化的模擬結(jié)果,并將源地的變化引入到生態(tài)阻力面模型中,另外考慮生態(tài)源地?cái)U(kuò)張能力因子修正MCR模型,構(gòu)建適合荒漠綠洲區(qū)的生態(tài)阻力模型。

    本文選擇荒漠綠洲區(qū)典型城市磴口縣為研究區(qū),采用數(shù)據(jù)同化技術(shù)中的集合卡爾曼濾波算法模擬修正生態(tài)源地的變化,構(gòu)建考慮源地變化、源地等級(jí)、距離、生態(tài)阻力/障礙4方面的生態(tài)阻力模型,對(duì)磴口縣的生態(tài)用地?cái)U(kuò)張進(jìn)行模擬。

    1 材料與方法

    1.1研究區(qū)概況

    磴口縣地處中國西北部(東經(jīng)107°05′,北緯40°13′),位于內(nèi)蒙古河套平原源頭,黃河中上游,背靠狼山山脈,西鄰烏蘭布和沙漠。磴口縣氣候干旱少雨,水資源較為短缺,土地沙漠化嚴(yán)重,土地退化嚴(yán)重,區(qū)域蒸發(fā)量大,導(dǎo)致土地鹽漬化程度深,境內(nèi)海拔高度1 030~2 046 m,整個(gè)地形除山區(qū)外,呈東南高西北低,逐漸傾斜。屬中溫帶大陸性季風(fēng)氣候,歷年平均風(fēng)速3.0 m/s,瞬間最大風(fēng)速28 m/s,多年平均降水143.9 mm,多年平均蒸發(fā)量2 327 mm,多年平均氣溫7.6℃,無霜期136 d。全縣有6個(gè)土類,10個(gè)亞類,31個(gè)土屬,258個(gè)土種。黃河流經(jīng)磴口縣52 km,年徑流量3.1×1010m3,共有水域面積約24.07 km2。河套地區(qū)地下水埋深0.5~3 m,沙區(qū)地下水埋深3~10 m,山前洪積扇地下水埋深3~30 m,相對(duì)豐富的地表水與地下水對(duì)磴口縣發(fā)展以及生態(tài)用地?cái)U(kuò)張?zhí)峁┝擞辛ΡU稀?/p>

    1.2數(shù)據(jù)來源與處理

    選取磴口縣夏季且少云的TM影像(2002、2007、2012年)和OLI影像(2015年)為研究對(duì)象,同時(shí)利用研究區(qū)2012年1∶5 萬植被分布圖、2015年1∶5 萬土地利用圖、空間分辨率30 m的DEM等作為輔助數(shù)據(jù)。利用ENVI 5.1軟件對(duì)影像進(jìn)行波段合成、圖像增強(qiáng)和幾何校正處理,選擇最大似然監(jiān)督分類法對(duì)4期遙感影像進(jìn)行目視解譯[12],提取磴口縣的生態(tài)用地類型信息,所提取的磴口縣生態(tài)用地為基礎(chǔ)性生態(tài)用地,其在改善環(huán)境、維持區(qū)域生態(tài)平衡方面具有不可替代的重要作用,包括湖泊、灘地、水庫坑塘等水域和自然荒漠灌林地、人工綠地等林草地。使用ArcMap 10.2進(jìn)行細(xì)碎斑塊處理,運(yùn)用疊加分析工具進(jìn)行空間數(shù)據(jù)分析,最終在ArcInfo 10.2環(huán)境下完成拓?fù)浜透腻e(cuò)處理。

    1.3技術(shù)路線

    利用所構(gòu)建的EnKF-MCRP模型進(jìn)行磴口縣生態(tài)用地?cái)U(kuò)張模擬,總體技術(shù)路線如圖1所示。

    圖1 總體技術(shù)路線Fig.1 Technical roadmap

    1.4基于EnKF的源地變化模擬

    磴口縣的生態(tài)用地類型包括人工型生態(tài)用地和自然型生態(tài)用地,其中能作為生態(tài)源地的生態(tài)用地類型為濕地型生態(tài)用地和荒漠植被型生態(tài)用地[2]。生態(tài)源地是隨著時(shí)間不斷變化的,且生態(tài)源地的變化會(huì)影響區(qū)域生態(tài)用地格局。通過引入集合卡爾曼濾波的數(shù)據(jù)同化方法,使用CA/Markov模型為數(shù)據(jù)同化模型算子模擬源地的變化,并將源地變化引入到生態(tài)用地?cái)U(kuò)張模擬中[13]。

    1.4.1集合卡爾曼濾波

    集合卡爾曼濾波(Ensemble Kalman filter, EnKF)將模型狀態(tài)預(yù)報(bào)看成近似隨機(jī)動(dòng)態(tài)預(yù)報(bào),用一個(gè)狀態(tài)總體(設(shè)數(shù)目為N,即集合數(shù))代表隨機(jī)動(dòng)態(tài)預(yù)報(bào)中的概率密度函數(shù)[14],通過向前積分計(jì)算下一時(shí)刻狀態(tài)總體的概率密度函數(shù),并得到該時(shí)刻的統(tǒng)計(jì)特性(如均值與協(xié)方差)[15-16]。

    集合卡爾曼濾波算法計(jì)算步驟如下:

    (1)模型狀態(tài)集合初始化為

    (1)

    式中n——模型狀態(tài)變量的個(gè)數(shù)

    k——時(shí)間刻度

    (2)計(jì)算每個(gè)狀態(tài)變量在第k+1時(shí)刻的預(yù)報(bào)值

    (2)

    ωk——期望為0、方差為Wk的高斯白噪聲

    Wk——模型的誤差方差矩陣

    M()為模型算子,本文指源的演變CA/Markov模型,在源地變化研究中,將源地變化過程看作Markov過程[17],如利用2002年和2007年的源地?cái)?shù)據(jù)預(yù)測2012年的源地空間分布,并與2012年源地實(shí)際數(shù)據(jù)進(jìn)行對(duì)比,即2007年的源地斑塊對(duì)應(yīng)于Markov過程中的可能狀態(tài),它只與2002年的源地斑塊狀態(tài)相關(guān),不同等級(jí)源地相互轉(zhuǎn)換的面積數(shù)量或比例即為狀態(tài)轉(zhuǎn)移概率[18-19]。

    (3)計(jì)算k+1時(shí)刻的卡爾曼增益矩陣Kk+1為

    (3)

    其中

    H——觀測算子

    R——預(yù)測誤差協(xié)方差矩陣

    (νk+1~N(0,R))

    (4)

    (5)

    其中

    (6)

    式中Yk+1——k+1時(shí)刻的狀態(tài)變量的觀測值

    νk+1——期望為0、方差為R的高斯白噪聲

    (5)判斷是否到結(jié)束時(shí)刻,即目標(biāo)年份。如果未到結(jié)束時(shí)刻,則返回步驟(2),否則結(jié)束。

    1.4.2EnKF與CA/Markov模型實(shí)現(xiàn)

    利用集合卡爾曼濾波數(shù)據(jù)同化算法融合高分辨率遙感影像源地的觀測數(shù)據(jù)得出同化結(jié)果,其輸入數(shù)據(jù)為非布爾型,而CA/Markov模型模擬元胞的結(jié)果狀態(tài)只有兩種,即源地與非源地[13-20]。故將整個(gè)研究區(qū)劃分為若干個(gè)正方格區(qū)域,每個(gè)方格里有若干個(gè)像元(元胞),得出每一個(gè)方格內(nèi)的源地元胞擴(kuò)張強(qiáng)度。

    (7)

    式中con(sij=source)——區(qū)域內(nèi)源地元胞的數(shù)量

    利用ArcGIS軟件的Fishnet模塊將全區(qū)分為513個(gè)規(guī)則方塊,從中選出66個(gè)觀測點(diǎn)和56個(gè)檢測點(diǎn)(圖2)。

    圖2 觀測點(diǎn)與檢測點(diǎn)空間分布圖 Fig.2 Spatial distribution of observation points and detection points

    誤差統(tǒng)計(jì)方式采用方差(VSE),其計(jì)算式為

    (8)

    式中y——真實(shí)值

    x——同化值或預(yù)測值,y和x維數(shù)一致

    1.5MCRP模型

    生態(tài)用地?cái)U(kuò)展可以看作是生態(tài)綠色空間對(duì)其他用地空間的競爭性控制過程,且這種控制和覆蓋必須通過克服阻力來實(shí)現(xiàn)[5,21-22],這樣生態(tài)用地?cái)U(kuò)展就可以模擬為從源到匯克服阻力做功的水平過程[23]?;旧鷳B(tài)阻力模型,即MCR模型的基本公式為

    (9)

    式中VMCR——最小累積阻力面值

    fmin——某土地單元對(duì)不同的生態(tài)用地取累積阻力最小值Dij——從生態(tài)用地j到土地單元i的空間距離

    Ri——用地單元i對(duì)運(yùn)動(dòng)過程的阻力系數(shù)

    考慮不同生態(tài)源地有不同的生態(tài)擴(kuò)張能力,即不同生態(tài)源地的影響力是不同的,引入不同等級(jí)生態(tài)源地的擴(kuò)張能力因子Pj,修正MCR模型,為簡化計(jì)算,按照生態(tài)源地斑塊的面積對(duì)生態(tài)源地進(jìn)行分級(jí),并賦予擴(kuò)張能力值。

    利用EnKF-CA/Markov模型預(yù)測2017年的生態(tài)源地,利用EnKF算法將預(yù)測的2017年的生態(tài)源地與2012年的現(xiàn)狀生態(tài)源地進(jìn)行數(shù)據(jù)同化,得到同化后的生態(tài)源地空間分布值,該值既包含2012年的生態(tài)源地信息,又包含生態(tài)源地的未來變化趨勢,即將源地的變化引入到生態(tài)阻力模型中。

    修正后建立源地?cái)U(kuò)張能力約束下的生態(tài)阻力模型——MCRP模型(Minimal cumulative resistance power)考慮了生態(tài)用地?cái)U(kuò)張的4方面因素,即生態(tài)源地變化、生態(tài)源地等級(jí)、距離和基面阻力特征,修正后公式為

    (10)

    式中VMCRP——生態(tài)用地?cái)U(kuò)展最小生態(tài)累積阻力面值

    Pj——生態(tài)用地j所屬等級(jí)的擴(kuò)張能力因子

    生態(tài)用地的等級(jí)越高,擴(kuò)張能力越強(qiáng)。

    2 結(jié)果與分析

    2.1源地變化模擬與分級(jí)

    根據(jù)2002、2007、2012年影像解譯數(shù)據(jù)結(jié)合磴口縣實(shí)際情況提取生態(tài)源地,在荒漠綠洲區(qū),濕地的生態(tài)作用巨大,故將所有水域提取出來均作為源地,在磴口縣荒漠灌林地的防風(fēng)固沙作用顯著,生態(tài)作用十分重要,故將面積大于0.1 km2的荒漠灌林地提取出來作為源地。

    利用EnKF以CA/Markov模型為數(shù)據(jù)同化算子模擬源地的時(shí)空演變,集合大小設(shè)為30,數(shù)據(jù)同化算子參數(shù)主要包括距離因子、地形因子、水文因子,如圖3所示。

    圖3 各類適宜性因子評(píng)價(jià)圖Fig.3 Evaluation of various types of suitability factors

    分別利用CA/Markov模型和EnKF-CA/Markov模型,基于2002、2007年的源地?cái)?shù)據(jù),預(yù)測2012年的源地空間分布(圖4b、4c),與2012年的生態(tài)用地實(shí)際情況進(jìn)行對(duì)比驗(yàn)證(圖4a)。

    圖4 2012年生態(tài)源地實(shí)際與模型模擬空間分布圖Fig.4 Actual and model simulated spatial distribution map of ecological source in 2012

    利用IDRISI軟件進(jìn)行逐像元對(duì)比分析得出EnKF-CA/Markov模型的模擬總精度為82.4%,模擬結(jié)果方差為0.6,Kappa系數(shù)為0.513。傳統(tǒng)CA/Markov模型的模擬總精度為65.4%,模擬結(jié)果方差為1.3,Kappa系數(shù)為0.342。在模擬過程中,利用數(shù)據(jù)同化技術(shù)修正的模擬結(jié)果在精度與方差上都有改善,本研究通過引入集合卡爾曼濾波,模型能夠減少誤差的積累,提高模擬精度。

    圖5 生態(tài)源地及等級(jí)劃分空間分布圖Fig.5 Distribution map of ecological sources and classification of space

    利用EnKF-CA/Markov模型,基于2007、2012年的源地預(yù)測2017年的源地空間分布情況,得到2017年源地空間分布如圖5a所示。將2017年的生態(tài)源地作為EnKF模型的觀測值,與2012年的實(shí)際生態(tài)源地進(jìn)行數(shù)據(jù)同化,得到2017、2012年的生態(tài)源地同化數(shù)據(jù)(圖5b),其中既包括2012年的生態(tài)源地信息,也包含生態(tài)用地變化至2017年的變化趨勢。

    進(jìn)而根據(jù)源地?cái)U(kuò)展能力,對(duì)源地進(jìn)行等級(jí)劃分,共分為5個(gè)等級(jí)(圖5c),級(jí)別越高表明源地的擴(kuò)張能力因子越大,5級(jí)源地主要分布在磴口縣包勒浩特嘎查、巴音溫都爾嘎查、溫都爾毛道嘎查、烏蘭布和農(nóng)場耕地外圍,自磴口縣東北貫穿至磴口縣西南,與奈倫湖相連。4級(jí)生態(tài)源地主要分布在巴音博日那嘎查、哈騰套海農(nóng)場、納林套海農(nóng)場的耕地外圍。3級(jí)生態(tài)源地主要分布在磴口縣防沙林場、奈倫湖、包爾蓋農(nóng)場附近。3、4、5級(jí)生態(tài)源地的空間布局形成東北-西南、西北-西南的沙漠化防護(hù)格局,對(duì)于防治烏蘭布和沙漠的擴(kuò)張具有重要意義。

    2.2生態(tài)阻力面構(gòu)建

    由于基面特性差異,生態(tài)用地在擴(kuò)張過程中所受的阻力是不同的,阻力系數(shù)反映了生態(tài)用地控制和覆蓋其他用地類型難易程度。將基面阻力分為2類:生態(tài)障礙和生態(tài)阻力。生態(tài)障礙是生態(tài)用地?zé)o法擴(kuò)展到的地塊,如建設(shè)用地、基本農(nóng)田等,構(gòu)成生態(tài)用地?cái)U(kuò)張的剛性約束,阻力系數(shù)無窮大。生態(tài)阻力是對(duì)生態(tài)用地?cái)U(kuò)張有不同等級(jí)阻礙作用的地塊,考慮地塊表面的固有屬性和附加屬性構(gòu)建了生態(tài)阻力評(píng)價(jià)體系,將生態(tài)阻力分為5個(gè)等級(jí),1級(jí)表示阻力最小,5級(jí)表示阻力最大。固有屬性包括用地類型和地質(zhì)地貌特征。附加屬性包括生態(tài)系統(tǒng)服務(wù)價(jià)值和NDVI-水體距離。所構(gòu)建的阻力評(píng)價(jià)體系如表1、2所示。

    首先確定生態(tài)用地空間擴(kuò)張障礙,第①類是基本農(nóng)田,是保證磴口縣區(qū)域經(jīng)濟(jì)、社會(huì)穩(wěn)定的保障,面積為796.28 km2。第②類是建設(shè)用地,是生態(tài)用地?cái)U(kuò)張的剛性限制。生態(tài)障礙分布圖如圖6a所示。

    進(jìn)而確定生態(tài)阻力面,根據(jù)生態(tài)阻力評(píng)價(jià)體系,使用MCRP模型算法,利用ArcGIS軟件中的cost-distance模塊分別生成5個(gè)等級(jí)源地生態(tài)阻力面,進(jìn)行疊加計(jì)算生成生態(tài)用地?cái)U(kuò)張生態(tài)阻力面,如圖6b所示,累計(jì)阻力值最小為0,位于生態(tài)用地內(nèi)部,累計(jì)阻力值最大為77 780.8,位于烏蘭布和沙漠內(nèi)部。由累計(jì)阻力面可以看出,生態(tài)用地外圍阻力均較大,形成明顯的累計(jì)阻力“山脊線”,這是由于在荒漠灌林地外圍是大面積的沙漠,累計(jì)阻力低的區(qū)域較為破碎,整體性差,嚴(yán)重影響景觀流的流動(dòng)。干旱區(qū)由于生態(tài)環(huán)境脆弱,荒漠化嚴(yán)重使得景觀生態(tài)流流動(dòng)不暢,固有屬性導(dǎo)致其景觀流阻力較大,而城市化的加劇導(dǎo)致了人工不透水表面的增加,進(jìn)一步加劇了景觀破碎,增加了景觀流動(dòng)阻力。

    表1 生態(tài)阻力評(píng)價(jià)體系Tab.1 Evaluation system of ecological resistance

    表2 NDVI與水體距離生態(tài)阻力評(píng)價(jià)體系Tab.2 Evaluation system of ecological resistance between NDVI and water distance

    圖6 生態(tài)障礙、阻力面及2015年生態(tài)用地空間分布圖Fig.6 Ecological barrier, resistance surface and spatial distribution of ecological land in 2015

    2.3生態(tài)用地?cái)U(kuò)張模擬

    利用ArcGIS 10.2軟件,使用python腳本語言編寫程序?qū)崿F(xiàn)模型的模擬。以2012年為基期數(shù)據(jù),設(shè)置時(shí)間閾值為3 a,分別使用2012年實(shí)際生態(tài)源地、經(jīng)過數(shù)據(jù)同化后的生態(tài)源地?cái)?shù)據(jù),利用生成的阻力面使用MCR模型和MCRP模型分別進(jìn)行生態(tài)用地?cái)U(kuò)張模擬,即模擬得出至2015年時(shí)的生態(tài)用地情況(圖6c),為便于驗(yàn)證,以2015年的生態(tài)用地面積為檢驗(yàn)標(biāo)準(zhǔn),磴口縣2015年的濕地型生態(tài)用地共137.92 km2。荒漠植被型、天然綠洲型、城鎮(zhèn)綠洲型、夾荒地型生態(tài)用地總計(jì)為1 042.96 km2。

    圖7 生態(tài)用地模型模擬空間分布圖Fig.7 Spatial distribution map of ecological land use through model simulation

    經(jīng)過迭代計(jì)算,利用EnKF-MCRP模型模擬的2015年生態(tài)用地面積共1 260.95 km2(圖7a), MCRP模型模擬的生態(tài)用地面積為983.57 km2(圖7b),EnKF-MCR模型模擬的生態(tài)用地面積為1 631.35 km2(圖7c)。與2015年的實(shí)際生態(tài)用地面積相對(duì)比,利用數(shù)據(jù)同化技術(shù)加入多時(shí)期生態(tài)源地變化的MCRP模型精度最高,逐像元對(duì)比方差為0.4,且整體的生態(tài)用地空間布局接近實(shí)際情況,利用MCRP模型的模擬精度其次,方差為1.7??紤]生態(tài)源地的變化后,MCRP模型模擬結(jié)果比2015年生態(tài)用地真實(shí)面積略大,不使用數(shù)據(jù)同化技術(shù)的MCRP模型模擬結(jié)果要比實(shí)際偏小,說明數(shù)據(jù)同化技術(shù)能夠提高模型模擬精度,將生態(tài)源地的變化作為模型因素十分必要。與EnKF-MCRP模型相比,EnKF-MCR模型的模擬結(jié)果精度較低,逐像元對(duì)比方差為3.1,模擬面積比真實(shí)值高450.37 km2,表明在構(gòu)建生態(tài)阻力面時(shí)一方面要考慮不同等級(jí)源地的生態(tài)阻力,另一方面要考慮生態(tài)障礙的剛性限制,由模擬結(jié)果圖可知,未考慮生態(tài)障礙的生態(tài)用地?cái)U(kuò)張中明顯侵占了基本農(nóng)田和建設(shè)用地。

    對(duì)56個(gè)檢測點(diǎn)的用地類型與實(shí)際用地類型進(jìn)行對(duì)比驗(yàn)證,其中相一致的賦值為1,不一致的賦值為2,56個(gè)檢測點(diǎn)的用地類型一致性對(duì)比發(fā)現(xiàn),EnKF-MCRP模型模擬結(jié)果與實(shí)際情況相符程度最高,MCRP模擬結(jié)果有11個(gè)點(diǎn)與實(shí)際不一致,EnKF-MCR模擬結(jié)果有19個(gè)點(diǎn)與實(shí)際不一致。

    在56個(gè)檢測點(diǎn)中選出16個(gè)具有代表性的檢測點(diǎn)并對(duì)其模擬結(jié)果進(jìn)行重點(diǎn)分析(表3),其中檢測點(diǎn)1、3、4、5、6、11、14、15、16處EnKF-MCRP、MCRP、EnKF-MCR模型的模擬結(jié)果均與實(shí)際相符合,經(jīng)過同化作用后,同化的結(jié)果較明顯,特別是檢測點(diǎn)2,2015年檢測點(diǎn)2處為生態(tài)用地,EnKF-MCRP、EnKF-MCR模型的模擬結(jié)果與實(shí)際相符,但是MCRP模型模擬結(jié)果與實(shí)際不符合,表明數(shù)據(jù)同化對(duì)于模擬精度的提高作用明顯,檢測點(diǎn)10、12位于生態(tài)用地的邊緣,結(jié)合數(shù)據(jù)同化技術(shù)的模型模擬結(jié)果均與實(shí)際相符合,未利用數(shù)據(jù)同化技術(shù)的MCRP模型模擬結(jié)果則精度較低,若迭代年份改變可能會(huì)導(dǎo)致模擬結(jié)果錯(cuò)誤,總體來看,整個(gè)同化過程中各點(diǎn)均能從觀測點(diǎn)中獲取正確的增益信息,因而能夠得出一個(gè)較好的誤差相關(guān)矩陣,作用到增益矩陣中,從而得出較好的同化值。所以數(shù)據(jù)同化能夠根據(jù)觀測點(diǎn)的信息得出比模擬值更接近真實(shí)值的同化值。故引入集合卡爾曼濾波的數(shù)據(jù)同化方法到模型中,能取得比較好的模擬結(jié)果。檢測點(diǎn)7、8、9處于生態(tài)障礙中,檢測點(diǎn)13位于生態(tài)障礙邊緣,EnKF-MCRP、MCRP模型的模擬結(jié)果均與實(shí)際相符,EnKF-MCR模型的模擬結(jié)果為生態(tài)用地,表明生態(tài)用地發(fā)生了過度擴(kuò)張,在構(gòu)建生態(tài)阻力面模型時(shí)考慮不同等級(jí)生態(tài)源地的生態(tài)障礙十分必要。

    表3 部分檢測點(diǎn)用地類型Tab.3 Land type for some detection points

    注:1為生態(tài)用地,2為生態(tài)障礙,3為沙漠,4為非生態(tài)用地。

    3 結(jié)論

    (1) 磴口縣生態(tài)源地包括濕地型生態(tài)用地和荒漠植被型生態(tài)用地,生態(tài)源地對(duì)于維持區(qū)域生態(tài)用地穩(wěn)定發(fā)展有重要作用。本文將數(shù)據(jù)同化技術(shù)應(yīng)用于生態(tài)源地的變化模擬中,將生態(tài)源地的變化趨勢量化,EnKF-CA/Markov模型的模擬總精度達(dá)到82.4%,通過引入集合卡爾曼濾波,能夠減少誤差的積累,提高模擬的精度。

    (2) 根據(jù)生態(tài)源地的擴(kuò)展能力,將磴口縣生態(tài)源地劃分為5個(gè)等級(jí),研究表明:磴口縣的3、4、5級(jí)生態(tài)源地的空間布局形成東北-西南、西北-西南的沙漠化防護(hù)格局,對(duì)于防治烏蘭布和沙漠的擴(kuò)張具有重要意義,對(duì)于維持區(qū)域生態(tài)環(huán)境穩(wěn)定具有重大作用。

    (3) 在構(gòu)建生態(tài)阻力面時(shí)考慮生態(tài)障礙和生態(tài)阻力兩個(gè)方面,將生態(tài)障礙引入阻力面構(gòu)建中能夠防止生態(tài)用地的過度擴(kuò)張,考慮生態(tài)源地變化、生態(tài)源地等級(jí)、距離、基面阻力4方面構(gòu)建MCRP模型,累計(jì)阻力值最小為0,最大為77 780.8,生態(tài)障礙阻力為無窮大。結(jié)果表明引入數(shù)據(jù)同化技術(shù)的EnKF-MCRP模型的生態(tài)用地?cái)U(kuò)張模擬精度最高,生態(tài)用地面積與空間布局最接近實(shí)際情況,數(shù)據(jù)同化對(duì)于模擬精度的提高作用明顯。

    1周銳,王新軍,蘇海龍,等.平頂山新區(qū)生態(tài)用地的識(shí)別與安全格局構(gòu)建[J].生態(tài)學(xué)報(bào),2015,35(6):2003-2012.

    ZHOU Rui, WANG Xinjun, SU Hailong, et al. Identification and security pattern of ecological land in Pingdingshan newly developed area[J]. Acta Ecologica Sinica,2015,35(6):2003-2012.(in Chinese)

    2張紅旗,王立新,賈寶全.西北干旱區(qū)生態(tài)用地概念及其功能分類研究[J].中國生態(tài)農(nóng)業(yè)學(xué)報(bào),2004,12(2):5-8.

    ZHANG Hongqi, WANG Lixin, JIA Baoquan. A conception of ecological land use and its function classification in arid area in Northwest China[J]. Chinese Journal of Eco-Agriculture,2004,12(2):5-8.(in Chinese)

    3岳德鵬,王計(jì)平,劉永兵,等.京郊西北地區(qū)農(nóng)地利用與景觀格局時(shí)空特征[J].農(nóng)業(yè)工程學(xué)報(bào),2008,24(6):89-95.

    YUE Depeng, WANG Jiping, LIU Yongbing, et al. Spatial and temporal characteristics of agricultural land-use and landscape pattern changes in northwest of Beijing[J]. Transactions of the CSAE,2008,24(6):89-95.(in Chinese)

    4梁曼,孫璇,楊振.土地景觀生態(tài)流空間通達(dá)性評(píng)價(jià)——以荊州市為例[J].河南科學(xué),2014,32(8):1654-1658.

    LIANG Man, SUN Xuan, YANG Zhen. Evaluation on spatial accessibility of ecological circulation of landscape—a case of Jingzhou City[J]. Henan Science,2014,32(8):1654-1658.(in Chinese)

    5葉玉瑤,蘇泳嫻,張虹鷗,等.生態(tài)阻力面模型構(gòu)建及其在城市擴(kuò)展模擬中的應(yīng)用[J].地理學(xué)報(bào),2014,69(4):485-496.

    YE Yuyao, SU Yongxian, ZHANG Hong’ou, et al. Ecological resistance surface model and its application in urban expansion simulations[J]. Acta Geographica Sinica,2014,69(4):485-496.(in Chinese)

    6YU Q, YUE D P, WANG J P, et al. The optimization of urban ecological infrastructure network based on the changes of county landscape patterns: a typical case study of ecological fragile zone located at Deng Kou (Inner Mongolia)[J]. Journal of Cleaner Production, doi: 10.1016/j.jclepro.2016.05.014.

    7LI F, Ye Y P, SONG B W. Evaluation of urban suitable ecological land based on the minimum cumulative resistance model: a case study from Changzhou, China[J]. Ecological Modelling,2015,318:194-203.

    8張亦漢,黎夏,劉小平,等. 基于數(shù)據(jù)同化的元胞自動(dòng)機(jī)[J].遙感學(xué)報(bào),2011,15(3):475-491.

    ZHANG Yihan, LI Xia, LIU Xiaoping, et al. The CA model based on data assimilation[J]. Journal of Remote Sensing,2011,15(3):475-491.(in Chinese)

    9桂預(yù)風(fēng),紀(jì)元,楊慧. 基于粒子濾波的隨機(jī)元胞自動(dòng)機(jī)城市用地?cái)U(kuò)張模型研究[J].數(shù)學(xué)的實(shí)踐與認(rèn)識(shí),2015,45(17):222-227.

    GUI Yufeng, JI Yuan, YANG Hui. The research of urban sprawl with the stochastic cellular automata model based on particle filter[J]. Mathematics in Practice and Theory,2015,45(17):222-227.(in Chinese)

    10張生雷,陳良富,蘇林,等.一種基于數(shù)據(jù)同化優(yōu)化陸面過程模式參數(shù)化方案的方法[J].中國科學(xué):地球科學(xué),2015,45(12):1907-1922.

    ZHANG Shenglei, CHEN Liangfu, SU Lin, et al. A data assimilation-based method for optimizing parameterization schemes in a land surface process model[J]. Science China: Earth Sciences,2015,45(12):1907-1922.(in Chinese)

    11黃健熙,李昕璐,劉帝佑,等. 順序同化不同時(shí)空分辨率LAI的冬小麥估產(chǎn)對(duì)比研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2015,46(1):240-248.

    HUANG Jianxi, LI Xinlu, LIU Diyou,et al. Comparison of winter wheat yield estimation by sequential assimilation of three different spatial-temporal resolution remotely sensed LAI datasets[J]. Transactions of the Chinese Society for Agricultural Machinery,2015,46(1):240-248.(in Chinese)

    12李保杰.礦區(qū)土地景觀格局演變及其生態(tài)效應(yīng)研究——以徐州市賈汪礦區(qū)為例[D].徐州:中國礦業(yè)大學(xué),2014.

    LI Baojie. Study on landscape pattern changes and ecological effects in mining area—A case study in Jiawang mining area [D]. Xuzhou: China University of Mining and Technology,2014.(in Chinese)

    13張亦漢,喬紀(jì)綱,艾彬.基于集合卡爾曼濾波動(dòng)態(tài)優(yōu)化CA模型參數(shù)的方法[J].測繪學(xué)報(bào),2013,42(1):123-130.

    ZHANG Yihan, QIAO Jigang, AI Bin. Parameter optimization for CA model using ensemble Kalman filter[J]. Acta Geodaetica et Cartographica Sinica,2013,42(1):123-130.(in Chinese)

    14王文,劉永偉,寇小華,等.基于集合卡爾曼濾波和HYDRUS-1D模型的土壤剖面含水量同化試驗(yàn)[J].水利學(xué)報(bào),2012,43(11):1302-1311.

    WANG Wen, LIU Yongwei, KOU Xiaohua,et al. EnKF and HYDRUS-1D based data assimilation experiments for improving soil moisture profile prediction[J]. Journal of Hydraulic Engineering,2012,43(11):1302-1311.(in Chinese)

    15李淵,基于數(shù)據(jù)同化的太湖葉綠素濃度遙感估算[D].南京:南京師范大學(xué),2014.

    LI Yuan. Chlorophyll a concentration remote evaluation in Lake Taihu based on data assimilation[D]. Nanjing: Nanjing Normal University,2014.(in Chinese)

    16高斌斌.集合Kalman濾波在T106L19中期數(shù)值預(yù)報(bào)譜模式中的應(yīng)用研究[D].南京:南京信息工程大學(xué),2008.

    GAO Binbin. Application research on the ensemble Kalman filter (EnKF) with a medium-range numerical weather prediction (NWP) spectral model at a T106L19 resolution[D]. Nanjing: Nanjing University of Information Science & Technology, 2008.(in Chinese)

    17HALMY M W, GESSLER P E, HICKE J A, et al. Land use/land cover change detection and prediction in the north-western coastal desert of Egypt using Markov-CA[J]. Applied Geography,2015,63:101-112.

    18張瀅,丁建麗.綠洲土地利用變化未來趨勢預(yù)測及其調(diào)控研究[J].干旱區(qū)資源與環(huán)境,2006,20(6):29-35.

    ZHANG Ying, DING Jianli. Study on forecasting and regulation of future trend of oasis land use[J]. Journal of Arid Land Resources and Environment,2006,20(6):29-35.(in Chinese)

    19王友生,余新曉,賀康寧,等.基于CA-Markov模型的藉河流域土地利用變化動(dòng)態(tài)模擬[J].農(nóng)業(yè)工程學(xué)報(bào),2011,27(12):330-336.

    WANG Yousheng, YU Xinxiao, HE Kangning, et al. Dynamic simulation of land use change in Ji He watershed based on CA-Markov model[J]. Transactions of the CSAE,2011,27(12):330-336.(in Chinese)

    20趙冬玲,杜萌,楊建宇,等. 基于CA-Markov模型的土地利用演化模擬預(yù)測研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2016,47(3):278-285.

    ZHAO Dongling, DU Meng, YANG Jianyu, et al. Simulation and forecast study of land use change based on CA-Markov model[J]. Transactions of the Chinese Society for Agricultural Machinery,2016,47(3):278-285.(in Chinese)

    21尹發(fā)能.江漢平原四湖流域景觀生態(tài)規(guī)劃與流域生態(tài)管理研究[D].上海:華東師范大學(xué),2008.

    YIN Faneng. Study of landscape ecological Planning and drainage basin ecological management in the Four-Lake basin[D].Shanghai:East China Normal University,2008.(in Chinese)

    22RAY N, LEHMANN A, JOLY P. Modeling spatial distribution of amphibian populations: a GIS approach based on habitat matrix permeability[J]. Biodiversity & Conservation,2002,11(12):2143-2165.

    23HERZOG F, LAUSCH A, MULLER E, et al. Landscape metrics for assessment of landscape destruction and rehabilitation[J]. Environmental Management,2001,27(1):91-107.

    Simulation on Ecological Land Use Expansion Based on EnKF-MCRP Model

    Yu Qiang1Yue Depeng1Di Yang2Zhang Qibin1Ma Huan1Li Yutong3

    (1.BeijingKeyLaboratoryofPrecisionForestry,BeijingForestryUniversity,Beijing100083,China2.DepartmentofGeographyUniversityofFlorida,GainesvilleFL32611,USA3.CenterforLandConsolidation,MinistryofLandandResources,Beijing100035,China)

    In the ecologically vulnerable area which locates in northwest arid and semiarid areas, ecological land is the important guarantee to maintain the security and stability of regional ecological environment. Studying the simulation of ecological land use change is of great significance. Accordingly, based on the typical ecologically vulnerable area—Dengkou County, this paper simulates ecological land use evolution of Dengkou County, using the four remote sensing image interpretation data of 2002, 2007, 2012 and 2015. Taking ecological eco-resistance barriers and EnKF-MCRP model into account, an ecological resistance surface was built. The EnKF-MCRP model was constructed to simulate the evolution of ecological sources considering the ecological sources change, ecological source level, distance and base surface resistance characteristics. The results showed that the combination of data assimilation and EnKF-CA/Markov model made a total accuracy of 82.4%, by using ensemble Kalman filter, the model can accumulate less errors and improve the accuracy of simulation, i.e., data assimilation can reduce the accumulation of errors. According to the expandability, ecological sources of Dengkou County were divided into five grades, of which the spatial layout of 3, 4, 5 grades formed the northeast-southwest and northwest-southwest pattern of desertification prevention. The building of EnKF-MCRP which takes the ecological source evolution into consideration made the highest precision. And the area of ecological sources and spatial distribution were the closest to the reality, of which the variance met 0.4. Different levels of ecological sources and ecological barriers were used to modify the model, which can improve the accuracy of the simulation results. This study could provide a scientific basis for the current and future ecological land planning and management.

    ecological land; data assimilation; EnKF-MCRP model; Dengkou County

    10.6041/j.issn.1000-1298.2016.09.039

    2016-06-16

    2016-06-30

    國家自然科學(xué)基金項(xiàng)目(41371189)和“十二五”國家科技支撐計(jì)劃項(xiàng)目(2012BAD16B00)

    于強(qiáng)(1987—),男,博士生,主要從事3S技術(shù)在生態(tài)環(huán)境中的應(yīng)用研究,E-mail: yuqiang@bjfu.edu.cn

    岳德鵬(1963—),男,教授,博士生導(dǎo)師,主要從事景觀生態(tài)學(xué)和土地評(píng)價(jià)研究,E-mail: yuedepeng@126.com

    K903

    A

    1000-1298(2016)09-0285-09

    猜你喜歡
    源地檢測點(diǎn)阻力
    核酸檢測點(diǎn)上,有最可愛的平江人
    騎馬做核酸
    鼻阻力測定在兒童OSA診療中的臨床作用
    零阻力
    英語文摘(2020年10期)2020-11-26 08:12:12
    移民與文化認(rèn)同:土家族民歌《吳幺姑》探析
    別讓摩擦成為學(xué)習(xí)的阻力
    飛行器FPGA檢測點(diǎn)優(yōu)化設(shè)置方法
    發(fā)源地
    氣候變化對(duì)渭河源地水文環(huán)境影響分析與探討
    不同種源地漆樹種子生物學(xué)特性研究
    少妇的逼水好多| 国产午夜福利久久久久久| 不卡一级毛片| 精品国产美女av久久久久小说| 乱人视频在线观看| 成人av在线播放网站| 狂野欧美激情性xxxx| 精品欧美国产一区二区三| 女警被强在线播放| 久久久久性生活片| 99久久综合精品五月天人人| 狂野欧美激情性xxxx| 成人亚洲精品av一区二区| 操出白浆在线播放| 国产精品美女特级片免费视频播放器| a在线观看视频网站| 日韩亚洲欧美综合| 国产中年淑女户外野战色| 国产成人欧美在线观看| 亚洲成a人片在线一区二区| www国产在线视频色| 国产成人啪精品午夜网站| 亚洲天堂国产精品一区在线| 免费在线观看日本一区| 窝窝影院91人妻| 欧美日本视频| 亚洲一区二区三区不卡视频| 欧美+日韩+精品| 一级毛片女人18水好多| 亚洲第一欧美日韩一区二区三区| 精品免费久久久久久久清纯| 一二三四社区在线视频社区8| 欧美性猛交黑人性爽| 91在线精品国自产拍蜜月 | svipshipincom国产片| 久久久久精品国产欧美久久久| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美国产在线观看| 看免费av毛片| 特级一级黄色大片| 欧美最新免费一区二区三区 | 窝窝影院91人妻| 午夜福利在线观看免费完整高清在 | 熟女电影av网| 热99在线观看视频| 国产精品久久久久久精品电影| 欧美av亚洲av综合av国产av| 久久精品人妻少妇| 一进一出抽搐gif免费好疼| 91久久精品电影网| 亚洲中文字幕一区二区三区有码在线看| 麻豆久久精品国产亚洲av| 午夜福利在线在线| 久久精品国产亚洲av香蕉五月| 亚洲av免费在线观看| 极品教师在线免费播放| 亚洲精华国产精华精| 精品一区二区三区视频在线 | 午夜福利在线观看免费完整高清在 | 首页视频小说图片口味搜索| 两性午夜刺激爽爽歪歪视频在线观看| 男女下面进入的视频免费午夜| а√天堂www在线а√下载| 久久草成人影院| 国产免费男女视频| 九九久久精品国产亚洲av麻豆| 夜夜爽天天搞| 黑人欧美特级aaaaaa片| 老熟妇仑乱视频hdxx| 国产成人福利小说| xxxwww97欧美| 又黄又爽又免费观看的视频| 高清日韩中文字幕在线| 成人av在线播放网站| 美女被艹到高潮喷水动态| 日本免费一区二区三区高清不卡| 天天添夜夜摸| 色av中文字幕| 午夜两性在线视频| 有码 亚洲区| 51国产日韩欧美| 中文字幕高清在线视频| 在线视频色国产色| 国产成人福利小说| 99精品久久久久人妻精品| 亚洲不卡免费看| 国产探花在线观看一区二区| 女人高潮潮喷娇喘18禁视频| 99久久无色码亚洲精品果冻| 淫妇啪啪啪对白视频| 国产免费一级a男人的天堂| 床上黄色一级片| 国产欧美日韩一区二区精品| 韩国av一区二区三区四区| 亚洲av日韩精品久久久久久密| 午夜福利在线观看吧| 给我免费播放毛片高清在线观看| 久久久久九九精品影院| 1024手机看黄色片| 在线免费观看不下载黄p国产 | 国产精品久久久久久久电影 | 婷婷丁香在线五月| 国产精品野战在线观看| 日本成人三级电影网站| 国产午夜精品久久久久久一区二区三区 | 亚洲电影在线观看av| 欧美日本亚洲视频在线播放| 国产黄a三级三级三级人| 少妇裸体淫交视频免费看高清| 国产v大片淫在线免费观看| 国产毛片a区久久久久| 久久婷婷人人爽人人干人人爱| 嫩草影院入口| 国产激情欧美一区二区| 村上凉子中文字幕在线| 成人无遮挡网站| 亚洲人成网站高清观看| 婷婷精品国产亚洲av在线| 精品久久久久久久毛片微露脸| 少妇的丰满在线观看| 又黄又爽又免费观看的视频| 长腿黑丝高跟| 99久国产av精品| 亚洲av二区三区四区| 一区福利在线观看| 欧美日韩一级在线毛片| 又黄又爽又免费观看的视频| 波多野结衣高清作品| 很黄的视频免费| 99国产精品一区二区三区| 在线观看一区二区三区| 免费看十八禁软件| 99久久精品一区二区三区| 最好的美女福利视频网| 欧美成人一区二区免费高清观看| 麻豆一二三区av精品| 无遮挡黄片免费观看| 最近最新中文字幕大全免费视频| 久久国产精品人妻蜜桃| 亚洲自拍偷在线| 久久6这里有精品| 在线播放国产精品三级| 欧美日韩综合久久久久久 | 国内少妇人妻偷人精品xxx网站| 露出奶头的视频| 女人被狂操c到高潮| 亚洲国产精品合色在线| 女生性感内裤真人,穿戴方法视频| xxxwww97欧美| 三级男女做爰猛烈吃奶摸视频| 国产精品一区二区三区四区久久| 两人在一起打扑克的视频| 日韩精品中文字幕看吧| 国产成人福利小说| 亚洲国产精品久久男人天堂| 好男人电影高清在线观看| 制服丝袜大香蕉在线| 国产久久久一区二区三区| 成年女人看的毛片在线观看| 免费观看人在逋| 久久久久精品国产欧美久久久| 不卡一级毛片| 两个人视频免费观看高清| 精品久久久久久久久久久久久| 精品一区二区三区人妻视频| 亚洲aⅴ乱码一区二区在线播放| 国产 一区 欧美 日韩| 久久精品人妻少妇| 亚洲精品影视一区二区三区av| 国产在线精品亚洲第一网站| 成人三级黄色视频| 一本久久中文字幕| 国产精品99久久久久久久久| 九九久久精品国产亚洲av麻豆| 露出奶头的视频| 欧美黑人巨大hd| 午夜福利视频1000在线观看| 可以在线观看的亚洲视频| 欧美极品一区二区三区四区| 可以在线观看毛片的网站| 欧美三级亚洲精品| av欧美777| 岛国在线免费视频观看| 首页视频小说图片口味搜索| 亚洲无线观看免费| 女人高潮潮喷娇喘18禁视频| 国产一区二区三区在线臀色熟女| av女优亚洲男人天堂| 久久国产精品人妻蜜桃| 午夜两性在线视频| 黄片小视频在线播放| 18禁裸乳无遮挡免费网站照片| 国产综合懂色| 免费av不卡在线播放| 国产精品三级大全| 观看免费一级毛片| 日韩 欧美 亚洲 中文字幕| 男人的好看免费观看在线视频| 成人一区二区视频在线观看| 国产精品香港三级国产av潘金莲| 亚洲欧美日韩高清专用| 岛国在线观看网站| 最近视频中文字幕2019在线8| 色尼玛亚洲综合影院| 伊人久久大香线蕉亚洲五| 91在线精品国自产拍蜜月 | 久久久久久大精品| 国产在视频线在精品| 成人国产一区最新在线观看| 黄片大片在线免费观看| 国产毛片a区久久久久| 19禁男女啪啪无遮挡网站| 欧美日韩亚洲国产一区二区在线观看| 免费在线观看影片大全网站| 国产精品久久视频播放| 午夜免费男女啪啪视频观看 | 99热这里只有是精品50| 色噜噜av男人的天堂激情| 99精品久久久久人妻精品| 亚洲av五月六月丁香网| 亚洲性夜色夜夜综合| 最新美女视频免费是黄的| 又粗又爽又猛毛片免费看| 欧美黑人巨大hd| 真人做人爱边吃奶动态| 久久久久九九精品影院| 国产v大片淫在线免费观看| 身体一侧抽搐| 久久天躁狠狠躁夜夜2o2o| 午夜免费男女啪啪视频观看 | 国产激情偷乱视频一区二区| 久久婷婷人人爽人人干人人爱| 无人区码免费观看不卡| 日本黄大片高清| 女人十人毛片免费观看3o分钟| 精品国内亚洲2022精品成人| 男人舔女人下体高潮全视频| 欧美丝袜亚洲另类 | 亚洲精品乱码久久久v下载方式 | 中文字幕人妻熟人妻熟丝袜美 | 国产高清videossex| 国产av在哪里看| 亚洲成av人片免费观看| 小说图片视频综合网站| 久久久久久久精品吃奶| 国产午夜精品论理片| 97超视频在线观看视频| 成人18禁在线播放| 五月伊人婷婷丁香| 日本撒尿小便嘘嘘汇集6| 午夜两性在线视频| 中亚洲国语对白在线视频| 国产亚洲精品综合一区在线观看| 国产av不卡久久| svipshipincom国产片| 亚洲专区中文字幕在线| 叶爱在线成人免费视频播放| 黄色片一级片一级黄色片| 可以在线观看毛片的网站| 日本免费一区二区三区高清不卡| 19禁男女啪啪无遮挡网站| 精品熟女少妇八av免费久了| 9191精品国产免费久久| 午夜福利18| 欧美zozozo另类| 黄色丝袜av网址大全| 动漫黄色视频在线观看| 波多野结衣高清作品| 一进一出抽搐动态| 欧美乱码精品一区二区三区| 亚洲精品成人久久久久久| 国产精品久久久久久人妻精品电影| 日韩欧美国产在线观看| 国内精品美女久久久久久| 19禁男女啪啪无遮挡网站| 国内精品久久久久久久电影| 国产熟女xx| 久久这里只有精品中国| 成年人黄色毛片网站| 亚洲国产日韩欧美精品在线观看 | 欧美丝袜亚洲另类 | 精品福利观看| 91在线精品国自产拍蜜月 | 男女视频在线观看网站免费| 99国产精品一区二区三区| 九九热线精品视视频播放| 久久人妻av系列| 少妇丰满av| 夜夜夜夜夜久久久久| 久久久久久人人人人人| 一级作爱视频免费观看| 午夜福利在线观看免费完整高清在 | 美女cb高潮喷水在线观看| 亚洲国产精品合色在线| 狠狠狠狠99中文字幕| 日韩大尺度精品在线看网址| 在线免费观看的www视频| 精品99又大又爽又粗少妇毛片 | ponron亚洲| 男插女下体视频免费在线播放| 岛国在线观看网站| 成人午夜高清在线视频| 757午夜福利合集在线观看| 综合色av麻豆| 真人做人爱边吃奶动态| 国产一区二区亚洲精品在线观看| 精品午夜福利视频在线观看一区| 欧美三级亚洲精品| 久久精品91无色码中文字幕| 欧美最新免费一区二区三区 | 禁无遮挡网站| 动漫黄色视频在线观看| 亚洲五月婷婷丁香| 九九在线视频观看精品| 成人无遮挡网站| 欧美不卡视频在线免费观看| 日日干狠狠操夜夜爽| 91九色精品人成在线观看| 丰满乱子伦码专区| 中文字幕高清在线视频| 露出奶头的视频| 免费在线观看日本一区| 免费高清视频大片| 国产精品电影一区二区三区| 国产一区二区激情短视频| 黑人欧美特级aaaaaa片| 国产精品一区二区三区四区久久| 午夜精品一区二区三区免费看| 高清日韩中文字幕在线| 国产高清三级在线| 天堂动漫精品| 看免费av毛片| 变态另类丝袜制服| 成人一区二区视频在线观看| 亚洲精品国产精品久久久不卡| 久久久精品大字幕| АⅤ资源中文在线天堂| 日日夜夜操网爽| 亚洲精品美女久久久久99蜜臀| 亚洲天堂国产精品一区在线| 国产精品免费一区二区三区在线| 九色国产91popny在线| 国产免费av片在线观看野外av| 精品99又大又爽又粗少妇毛片 | 女人高潮潮喷娇喘18禁视频| 美女高潮的动态| 精品人妻1区二区| 亚洲成人久久爱视频| 婷婷精品国产亚洲av| 国产高清视频在线观看网站| 国产高清三级在线| 国产91精品成人一区二区三区| 国产美女午夜福利| 不卡一级毛片| 久99久视频精品免费| 人妻夜夜爽99麻豆av| 俺也久久电影网| 国产视频一区二区在线看| 欧美成人一区二区免费高清观看| 全区人妻精品视频| 真人一进一出gif抽搐免费| 精品一区二区三区视频在线观看免费| 久久午夜亚洲精品久久| 91久久精品电影网| 国产成人aa在线观看| 国产乱人视频| 两人在一起打扑克的视频| av中文乱码字幕在线| 国产视频内射| 国产精品av视频在线免费观看| 美女高潮的动态| 无遮挡黄片免费观看| 久久久国产成人精品二区| 2021天堂中文幕一二区在线观| av片东京热男人的天堂| 最好的美女福利视频网| 婷婷丁香在线五月| 国产精品av视频在线免费观看| 日本三级黄在线观看| 国产午夜精品久久久久久一区二区三区 | 成人精品一区二区免费| 99国产综合亚洲精品| 久久久久免费精品人妻一区二区| 欧美日韩乱码在线| 成年版毛片免费区| 18禁裸乳无遮挡免费网站照片| 国产亚洲精品久久久久久毛片| 欧美成人一区二区免费高清观看| 亚洲精品在线美女| 一进一出抽搐动态| 国产野战对白在线观看| 老司机福利观看| 午夜a级毛片| 久久伊人香网站| www.色视频.com| 丰满人妻熟妇乱又伦精品不卡| 香蕉av资源在线| 久久精品国产99精品国产亚洲性色| 我的老师免费观看完整版| 国产精品精品国产色婷婷| 狂野欧美白嫩少妇大欣赏| 白带黄色成豆腐渣| 午夜免费观看网址| 在线免费观看不下载黄p国产 | 麻豆成人午夜福利视频| 日韩精品中文字幕看吧| 久久久久久久午夜电影| 精品人妻一区二区三区麻豆 | 热99在线观看视频| 中文字幕精品亚洲无线码一区| 变态另类成人亚洲欧美熟女| 18禁在线播放成人免费| 国产色爽女视频免费观看| 久久国产乱子伦精品免费另类| 别揉我奶头~嗯~啊~动态视频| 国产乱人伦免费视频| 国产精品久久视频播放| av在线蜜桃| 国产成人aa在线观看| 精品人妻1区二区| 狂野欧美激情性xxxx| 五月伊人婷婷丁香| 一本精品99久久精品77| 亚洲熟妇熟女久久| 亚洲午夜理论影院| 久久久久久久精品吃奶| 亚洲国产高清在线一区二区三| 国产欧美日韩一区二区精品| 无人区码免费观看不卡| 久久精品91蜜桃| 最近最新免费中文字幕在线| 不卡一级毛片| 国产日本99.免费观看| 91av网一区二区| 亚洲国产精品久久男人天堂| 日本精品一区二区三区蜜桃| 一级毛片高清免费大全| 精品免费久久久久久久清纯| 小说图片视频综合网站| 亚洲精华国产精华精| 欧美精品啪啪一区二区三区| 亚洲熟妇熟女久久| 久久人妻av系列| 国产三级在线视频| 毛片女人毛片| av福利片在线观看| 日日干狠狠操夜夜爽| 91九色精品人成在线观看| 国产一区二区在线av高清观看| 最新美女视频免费是黄的| 日本 av在线| 国产v大片淫在线免费观看| 99国产极品粉嫩在线观看| 操出白浆在线播放| 99久久久亚洲精品蜜臀av| 国产成人影院久久av| 国产精品99久久99久久久不卡| 男女之事视频高清在线观看| 在线播放无遮挡| 亚洲内射少妇av| 久久国产精品人妻蜜桃| 久久婷婷人人爽人人干人人爱| 精品日产1卡2卡| 精品久久久久久,| 国产亚洲欧美98| 一级a爱片免费观看的视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品美女久久久久99蜜臀| 亚洲成av人片在线播放无| 亚洲七黄色美女视频| 中文亚洲av片在线观看爽| 可以在线观看的亚洲视频| 国产成人系列免费观看| 观看免费一级毛片| 亚洲狠狠婷婷综合久久图片| 亚洲av二区三区四区| 亚洲欧美日韩高清专用| 久久国产乱子伦精品免费另类| 成人18禁在线播放| 亚洲精品粉嫩美女一区| 内射极品少妇av片p| 18美女黄网站色大片免费观看| 老汉色∧v一级毛片| 亚洲欧美激情综合另类| 久久精品夜夜夜夜夜久久蜜豆| 99热这里只有是精品50| av在线天堂中文字幕| 婷婷精品国产亚洲av| 别揉我奶头~嗯~啊~动态视频| 午夜福利18| 天堂影院成人在线观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲最大成人中文| 国产亚洲av嫩草精品影院| 成人国产一区最新在线观看| 亚洲人成网站在线播| 欧美黄色片欧美黄色片| 韩国av一区二区三区四区| 男女视频在线观看网站免费| 国产爱豆传媒在线观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲av成人av| 欧美精品啪啪一区二区三区| 国产精品嫩草影院av在线观看 | 国产伦一二天堂av在线观看| 俺也久久电影网| 国产精品影院久久| 日本在线视频免费播放| 国产亚洲欧美在线一区二区| www日本在线高清视频| 桃红色精品国产亚洲av| 久久人人精品亚洲av| 丰满人妻熟妇乱又伦精品不卡| 久久香蕉国产精品| 成年版毛片免费区| 女警被强在线播放| 在线观看一区二区三区| 国产av麻豆久久久久久久| 亚洲av熟女| 三级国产精品欧美在线观看| 女人被狂操c到高潮| 国产精品,欧美在线| 亚洲国产日韩欧美精品在线观看 | 在线观看av片永久免费下载| 色综合欧美亚洲国产小说| 国产成年人精品一区二区| 亚洲一区二区三区不卡视频| 欧美xxxx黑人xx丫x性爽| 午夜福利高清视频| 日韩欧美在线二视频| 琪琪午夜伦伦电影理论片6080| 成年版毛片免费区| 三级毛片av免费| 99久久精品热视频| 99久久99久久久精品蜜桃| 成人无遮挡网站| 久久久久亚洲av毛片大全| 色尼玛亚洲综合影院| 欧美性猛交黑人性爽| 三级男女做爰猛烈吃奶摸视频| 狂野欧美白嫩少妇大欣赏| 男女做爰动态图高潮gif福利片| 波多野结衣巨乳人妻| 亚洲精华国产精华精| 欧美又色又爽又黄视频| 一本精品99久久精品77| 日韩欧美免费精品| 精品一区二区三区人妻视频| 成人亚洲精品av一区二区| 国产一区二区在线av高清观看| 手机成人av网站| 在线观看免费午夜福利视频| 国产伦精品一区二区三区视频9 | 人人妻人人看人人澡| 一个人观看的视频www高清免费观看| 国内少妇人妻偷人精品xxx网站| 午夜福利成人在线免费观看| 岛国视频午夜一区免费看| 国产成+人综合+亚洲专区| 毛片女人毛片| 午夜福利欧美成人| 天堂网av新在线| 亚洲av熟女| 伊人久久大香线蕉亚洲五| 99久久成人亚洲精品观看| 综合色av麻豆| 欧美区成人在线视频| 首页视频小说图片口味搜索| 长腿黑丝高跟| 黄色日韩在线| 亚洲欧美激情综合另类| 88av欧美| www日本在线高清视频| 操出白浆在线播放| 日韩欧美一区二区三区在线观看| 搡老熟女国产l中国老女人| 久久九九热精品免费| 美女免费视频网站| 美女被艹到高潮喷水动态| 久久天躁狠狠躁夜夜2o2o| 欧美日本亚洲视频在线播放| 人妻丰满熟妇av一区二区三区| 亚洲黑人精品在线| 亚洲真实伦在线观看| 黑人欧美特级aaaaaa片| 女人十人毛片免费观看3o分钟| 欧美+亚洲+日韩+国产| 国产 一区 欧美 日韩| 亚洲欧美一区二区三区黑人| 日韩人妻高清精品专区| 制服人妻中文乱码| 国产黄色小视频在线观看| 久久久久久人人人人人| 久久久久精品国产欧美久久久| 在线观看美女被高潮喷水网站 | 97人妻精品一区二区三区麻豆| 看免费av毛片| 最近最新免费中文字幕在线| 村上凉子中文字幕在线| 国产欧美日韩精品一区二区| 在线观看午夜福利视频| 国产成人系列免费观看| 亚洲av成人精品一区久久| 免费人成视频x8x8入口观看| 国产精品综合久久久久久久免费| 欧美一级毛片孕妇| 综合色av麻豆| 亚洲一区二区三区色噜噜| 欧美另类亚洲清纯唯美| 欧美性猛交黑人性爽| 精品一区二区三区视频在线观看免费| 欧美日韩亚洲国产一区二区在线观看| 中文字幕熟女人妻在线| 麻豆久久精品国产亚洲av| 女人被狂操c到高潮|