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

    基于MCR-ANN-CA模型的包頭市生態(tài)用地演變模擬

    2019-03-06 09:18:22劉建華張啟斌YANGDi岳德鵬
    農(nóng)業(yè)機械學報 2019年2期
    關鍵詞:元胞包頭市鄰域

    劉建華 張啟斌 YANG Di 岳德鵬 于 強 楊 斕

    (1.北京林業(yè)大學精準林業(yè)北京市重點實驗室, 北京 100083; 2.弗羅里達大學地理系, 蓋恩斯維爾 FL 32611)

    0 引言

    《全國生態(tài)環(huán)境保護綱要(國發(fā)[2000]38號)》中提到了生態(tài)用地一詞,明確指出了生態(tài)用地的重要性[1]。眾多學者基于此不斷發(fā)展生態(tài)用地的概念和內(nèi)涵。一種觀點認為應該從土地類型上來劃分生態(tài)用地,即提供自然生態(tài)系統(tǒng)服務價值的土地都可以被視為生態(tài)用地。另一種觀點認為應該從土地的主體功能來區(qū)分,主張農(nóng)業(yè)用地不應作為生態(tài)用地[2-3]。根據(jù)前人理論成果,結合包頭市實地情況,本文將林地、草地、水體作為生態(tài)用地范圍。

    生態(tài)用地作為自然生態(tài)系統(tǒng)服務的基本載體,是解決城市建設用地擴張與生態(tài)保護矛盾的綜合途徑[4]。在我國西北干旱地區(qū),生態(tài)用地在防止土地沙漠化、水土保持等方面的生態(tài)功能更加突出,對于維持區(qū)域生態(tài)系統(tǒng)健康穩(wěn)定具有至關重要的作用[5-6]。對生態(tài)用地的空間變化進行模擬,可以為干旱地區(qū)生態(tài)用地科學規(guī)劃和生態(tài)安全格局構建提供重要參考價值。

    生態(tài)用地的演變模擬是土地利用演變模擬的一部分,當前土地利用格局的演變模擬主要包括兩種類型,一種是宏觀回顧模型,如邏輯回歸(Logistic regression)模型[7]、馬爾可夫(Markov)模型[8]、系統(tǒng)動力學(System dynamic)模型等[9-10],這些模型根據(jù)土地利用的時空數(shù)據(jù),量化土地利用在宏觀時空尺度上的變化規(guī)律,對未來情景進行預測。另一種是微觀預測模型,包括多智能體模型(Multi-agent)[11]、元胞自動機(Cellular automata, CA)[12]等,這種模型從微觀尺度入手,通過研究土地利用與生態(tài)、人類活動等因子的相互作用,預測宏觀上土地利用的未來演變。將兩種模型結合,可得到二者的混合模型,既可以考慮全局宏觀效應,又可以考慮多個因子對土地利用的微觀作用機制,相比單一模型具有較大優(yōu)勢。將宏觀回顧模型與元胞自動機結合,是這種模型的一種常見形式,如CA-Markov模型[13]。在這種模型中,有兩個關鍵點對其模擬精度有較大影響,一是元胞鄰域內(nèi)不同土地利用類型間轉變規(guī)則,即鄰域規(guī)則;二是地形、生態(tài)環(huán)境等因子對土地利用演變影響,即適宜性規(guī)則[14]。當前傳統(tǒng)元胞自動機鄰域規(guī)則制定方法較多,如支持向量機、多目標決策、Markov模型等[15],用戶往往難以選擇合適的規(guī)則定義方法,部分線性模型也難以根據(jù)鄰域特征精確模擬生態(tài)用地演化等非線性過程;基于邏輯回歸、多標準評價模型及層次分析法等制定的適宜性規(guī)則雖然考慮了多種因素對土地利用變化的影響,但多基于單個像元利用空間疊加手段獨立評價其適宜性,不能根據(jù)生態(tài)用地演化的過程量化研究區(qū)內(nèi)各景觀單元演化為生態(tài)用地的適宜性。

    針對上述問題,選取西北干旱地區(qū)典型城市包頭市為研究區(qū),基于源匯理論,利用最小累積耗費阻力(Minimum cumulative resistance, MCR)模型量化生態(tài)用地演化過程中,從“生態(tài)源地”到其他土地利用類型的適宜性,同時利用人工神經(jīng)網(wǎng)絡(Artificial neural network, ANN)提取CA鄰域轉換規(guī)則,構建MCR-ANN-CA模型對包頭市生態(tài)用地的演變過程進行模擬,以期為區(qū)域生態(tài)用地規(guī)劃及生態(tài)建設提供理論與方法支持。

    1 材料與方法

    1.1 研究區(qū)概況

    包頭市位于內(nèi)蒙古自治區(qū)西部(東經(jīng)109°13′~111°26′,北緯40°13′~42°44′),面積27 768 km2。包頭市深處內(nèi)陸,氣候為典型的溫帶干旱、半干旱大陸性氣候,冬季寒冷干燥、夏季炎熱多雨[16],年平均氣溫2.0~7.7℃,年均降水量175~400 mm,年均蒸發(fā)量為2 100~2 700 mm。包頭市可利用地表水總量為9×108m3,地下水補給量為8.6×109m3。黃河流經(jīng)包頭境內(nèi)214 km,水面寬130~458 m,最大流量6 400 m3/s,年平均徑流量為2.6×1011m3,是包頭市主要用水來源[17]。包頭市海拔976~2 317 m,地勢中間高南北低,北部丘陵、中部山地、南部平原分別占總面積的14.49%、75.51%和10%[18]。干旱的氣候條件與起伏的地貌特征使得包頭市生態(tài)環(huán)境較為脆弱,面臨較高的水土流失與土地沙漠化、荒漠化風險。近年來,包頭市城市建設規(guī)模擴張迅速,房地產(chǎn)開發(fā)項目、工業(yè)園區(qū)等建設項目不斷涌現(xiàn),導致草地、林地、濕地等生態(tài)用地遭到占用與破壞,生態(tài)風險有所升高[19]。

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

    本研究主要數(shù)據(jù)源包括包頭市土地利用數(shù)據(jù)(2006、2011、2016年)、歸一化植被指數(shù)(Normalized difference vegetation index, NDVI)、數(shù)字地面高程模型(Digital elevation model, DEM)、各區(qū)縣旗人口密度及工業(yè)園區(qū)分布數(shù)據(jù)。其中土地利用數(shù)據(jù)由Landsat-8遙感數(shù)據(jù)解譯得到,工業(yè)園區(qū)分布數(shù)據(jù)來自《包頭市土地利用總體規(guī)劃》和《包頭市城市總體規(guī)劃》;歸一化植被指數(shù)(NDVI)、數(shù)字地面高程模型(DEM)來自地理空間數(shù)據(jù)云(http:∥www.gscloud.cn);各區(qū)人口密度數(shù)據(jù)來自包頭市統(tǒng)計年鑒。

    根據(jù)研究需要,從遙感影像中提取耕地、林地、草地、建設用地、水體、未利用地并通過實地驗證確保其精度,其中林地、草地、水體為本文中生態(tài)用地范圍;利用DEM數(shù)據(jù)提取高程、坡度數(shù)據(jù)。

    1.3 最小累積耗費阻力模型

    MCR模型由KANPPEN提出,最早應用于物種遷徙過程研究,之后在物種保護、景觀格局分析等方面取得了廣泛應用[20]。該模型主要考慮3個因素,即“源”、阻力和累積代價,通過對3個因素的分析,對“源”克服阻力向外傳播所耗費的代價或者所做的功進行描述[21]。MCR模型的一般形式為

    R=f∑DijRi

    (1)

    式中R——最小累積阻力

    f——未知負函數(shù),表示最小累積阻力與生態(tài)適宜性的負相關關系

    Dij——從源j到景觀單元i的空間距離

    Ri——景觀單元i處的阻力

    本文將生態(tài)用地的演化過程看作生態(tài)用地對其他景觀的競爭性控制過程,且這種演化必須通過克服阻力實現(xiàn),這樣生態(tài)用地的演化過程就可以抽象為從源(現(xiàn)有生態(tài)用地斑塊)到匯(其他景觀單元)克服阻力做功的水平過程[22]。由于區(qū)域下墊面差異,不同空間位置的土地演化為生態(tài)用地的阻力是不同的,通過由“源”到當前像元的累積阻力可量化當前像元演化為生態(tài)用地的概率,即阻力越大,該像元演化為生態(tài)用地的概率越小,演化為其他用地類型的概率越大。本文利用MCR模型,綜合考慮土地利用、NDVI、坡度、政府規(guī)劃工業(yè)園區(qū)、人口密度、水體距離、高程因子構建累積耗費阻力面,利用該累積耗費阻力面構建CA模型的適宜性規(guī)則,對上述過程進行模擬,提高模型預測精度。

    1.4 元胞自動機模型

    CA模型是一種時間、空間、狀態(tài)都離散的動力學模型,具有明顯的時間與空間特征,CA模型的一般形式為[23]

    Si,t+1=f(Si,t,SN,t)

    (2)

    式中Si,t+1——元胞i在t+1時刻的狀態(tài)

    Si,t——元胞i在t時刻的狀態(tài)

    SN,t——元胞i的鄰域集合在t時刻的狀態(tài)

    f——轉換規(guī)則

    CA模型中,中心元胞在下一時刻的狀態(tài)是其在當前時刻狀態(tài)及其鄰域集合狀態(tài)的函數(shù),準確定義該函數(shù)對于CA模型的模擬精度具有關鍵作用[24-27]。由于生態(tài)用地的演化過程是一種非線性的復雜動力學過程,當前元胞及其鄰域狀態(tài)影響中心元胞過程很難用簡單的規(guī)則定義,因此本文采用人工神經(jīng)網(wǎng)絡模型提取CA模型轉換規(guī)則,對土地利用模擬過程中的現(xiàn)有規(guī)則進行改進[28-29],以提高CA模型的模擬精度。

    本研究中,元胞形狀為正方形的柵格像元,尺寸為30 m×30 m,采用3×3經(jīng)典摩爾鄰域定義鄰域空間,模型初始柵格數(shù)據(jù)中的元胞狀態(tài)定義為研究區(qū)景觀格局類型,依次為耕地、林地、草地、建設用地、水體、未利用地,根據(jù)中心元胞及鄰域元胞狀態(tài),利用ANN模型提取元胞轉換規(guī)則,利用該規(guī)則計算元胞演化為生態(tài)用地的概率。

    1.5 人工神經(jīng)網(wǎng)絡模型

    ANN模型是20世紀80年代以來人工智能領域興起的研究熱點[30]。ANN模型在復雜的非線性系統(tǒng)的模擬方面具有明顯優(yōu)勢,其自組織、自學習、聯(lián)想以及記憶的優(yōu)勢能夠有效簡化CA模型,從原始訓練數(shù)據(jù)中提取CA轉換規(guī)則,避免主觀因素影響,提高模擬精度[31]。

    由于ANN模型在解決此類非線性問題中的明顯優(yōu)勢,國內(nèi)外學者很早就開始了相關研究,嘗試利用ANN模型提取CA模型的轉換規(guī)則[32]。此類研究往往把鄰域規(guī)則和適宜性規(guī)則統(tǒng)一放入ANN模型進行規(guī)則提取,然而在模擬生態(tài)用地演化的CA模型中,相比適宜性規(guī)則,鄰域規(guī)則部分更難定義也更難被人類理解,因此本文利用ANN模型提取CA模型鄰域部分轉換規(guī)則,適宜性規(guī)則采用MCR模型構建。

    采用經(jīng)典BP神經(jīng)網(wǎng)絡算法,網(wǎng)絡結構共3層,分別是輸入層、隱含層、輸出層。模型的輸入層共7個節(jié)點,隱含層共5個節(jié)點,采用tansing激勵函數(shù),輸出層共1個節(jié)點,采用sigmoid激勵函數(shù),其輸出為中心元胞轉換為生態(tài)用地的概率,值域為0~1。輸入變量及其取值范圍如表1所示。

    表1 ANN模型輸入變量及取值范圍Tab.1 Input variables and their ranges of ANN model

    1.6 MCR-ANN-CA模型

    以CA模型為基礎框架,耦合MCR模型與ANN模型,構建MCR-ANN-CA模型用于包頭市生態(tài)用地的演化模擬,技術路線如圖1所示。模型模擬的基本步驟為:

    圖1 技術路線Fig.1 Technical roadmap

    (1)利用MCR模型,綜合考慮多種土地適宜性因子,構建研究區(qū)范圍內(nèi)各像元演化為生態(tài)用地的累積阻力面。

    (2)利用歷史數(shù)據(jù)的轉換情況及鄰域特征,構建訓練數(shù)據(jù),訓練ANN模型,提取CA模型鄰域規(guī)則。

    (3)以當前土地利用為CA模型的輸入數(shù)據(jù),針對任一元胞,利用訓練好的ANN模型,根據(jù)其鄰域狀態(tài)判斷其演化為生態(tài)用地的概率P。

    圖2 阻力因子的阻力分布Fig.2 Spatial distributions of resistance values of resistance factors

    (4)將P與累積阻力面中像元最大值相乘,其結果為S(取值范圍為0到累積阻力面最大值),若S大于當前元胞所對應的累積阻力,則該元胞演化為生態(tài)用地,否則,該元胞演化為非生態(tài)用地。

    (5)將模擬結果轉換為柵格圖像并輸出結果。

    2 結果與分析

    2.1 基于MCR模型的適宜性規(guī)則構建

    對土地利用、NDVI、坡度、政府規(guī)劃工業(yè)園區(qū)、人口密度、水體距離、高程等多種因子進行空間化與標準化處理,結果如圖2所示。其中將NDVI、坡度、人口密度、高程因子數(shù)據(jù)原始值進行歸一化處理作為其阻力。通過對不同用地類型阻力賦值,并對賦值結果歸一化處理得到土地利用數(shù)據(jù)。政府規(guī)劃工業(yè)園區(qū)阻力賦值為1,其他區(qū)域賦值為0,將數(shù)據(jù)柵格化后作為一個阻力因子。

    根據(jù)包頭市2016年土地利用數(shù)據(jù),提取包頭市生態(tài)用地范圍,如圖3a所示,由提取結果可知,包頭市生態(tài)源地主要分布在北部草原帶、大青山一帶及黃河沿岸,其他區(qū)域也有零星生態(tài)用地分布,大致形成了三屏多點的格局。

    利用GIS疊加分析,將圖2中各因子進行空間疊加,得到包頭市各個空間位置演化為生態(tài)用地的阻力面,結果如圖3b所示。基于該阻力面,采用MCR模型,對源地進行累積耗費阻力計算,量化由近及遠的空間范圍內(nèi)各像元演化為生態(tài)用地面臨的阻力,得到如圖3c所示的累積生態(tài)阻力面。由圖3可知,包頭市南北部生態(tài)阻力較低,中部累積生態(tài)阻力較高,尤其圖3c中A區(qū)域及B區(qū)域,累積生態(tài)阻力達到包頭市最高水平。

    圖3 包頭市生態(tài)用地演化累積生態(tài)阻力面構建Fig.3 Construction of ecological land transition minimum accumulative resistance surface of Baotou City

    2.2 基于ANN模型的CA鄰域規(guī)則提取

    基于研究區(qū)2006、2011年景觀格局分布數(shù)據(jù),通過Matlab隨機抽取3 000個景觀像元,統(tǒng)計2006年景觀格局分布數(shù)據(jù)中,以3 000個像元為中心的3×3鄰域內(nèi)的各輸入變量大小,形成ANN模型訓練數(shù)據(jù)庫。選其中70%進行模型訓練,15%作為驗證數(shù)據(jù)集,15%作為測試數(shù)據(jù)集。

    采用量化共軛梯度法對ANN進行訓練,模型經(jīng)57次迭代后均方根誤差達到最小,為0.12,如圖4所示,圖中圓圈處為最優(yōu)擬合點,此時驗證集的決定系數(shù)R2為0.90,模型擬合精度較好地滿足了本研究要求。

    圖5 MCR-ANN-CA與CA-Markov模擬結果對比Fig.5 Comparison of simulation results of MCR-ANN-CA model and CA-Markov model

    圖4 ANN模型均方根誤差Fig.4 Root mean square error of ANN model

    2.3 MCR-ANN-CA生態(tài)用地模擬

    將2011年景觀格局數(shù)據(jù)輸入訓練后的ANN模型,判斷各像元在鄰域規(guī)則影響下是否會演變?yōu)樯鷳B(tài)用地,進而根據(jù)MCR模型生成的累積耗費阻力面,判斷其最終演變方向,將模型模擬結果輸出為柵格圖像。包頭市2016年實際生態(tài)用地8 409.32 km2,MCR-ANN-CA模型的模擬結果中包頭市生態(tài)用地面積為8 670.01 km2,如圖5a、5b所示。

    2.4 模型模擬精度對比分析

    為驗證模型模擬精度,利用CA-Markov模型對包頭市2016年生態(tài)用地進行模擬,結果如圖5c所示。以2016年生態(tài)用地實際分布為參照,利用Idrisi Selva軟件中的Cross Tab模塊對2個模型的模擬結果進行分析,定量分析模型的模擬精度,結果如表2所示。表中KIA指數(shù)為卡帕一致性指數(shù)。

    表2 模型模擬精度評價Tab.2 Accuracy assessment of simulation result

    由表2可知,MCR-ANN-CA模型與CA-Markov模型模擬結果中,生態(tài)用地的面積均大于實際值,但總體上與2016年生態(tài)用地的實際分布保持了較高的一致性,其中MCR-ANN-CA模型的模擬精度略高于CA-Markov模型,二者的相對誤差分別為3.10%與5.31%,KIA指數(shù)分別為0.89和0.87。

    相比CA-Markov模型,MCR-ANN-CA模型通過ANN模型提取了元胞自動機鄰域內(nèi)的轉換規(guī)則,同時利用MCR模型構建累積耗費阻力面,對不同空間位置元胞演化為生態(tài)用地的阻力進行了量化,因此模擬精度得到了進一步提高。

    3 結論

    (1)利用ANN模型提取了元胞自動機的鄰域規(guī)則,同時利用MCR模型構建累積耗費阻力面,基于MCR-ANN-CA模型對包頭市生態(tài)用地演化情況進行模擬,模擬精度較高。

    (2)將MCR-ANN-CA模型與CA-Markov模型模擬結果進行對比,KIA指數(shù)分別為0.89和0.87,相對誤差分別為3.10%和5.31%,MCR-ANN-CA模型對包頭市生態(tài)用地的演化過程具有更高的模擬精度。

    猜你喜歡
    元胞包頭市鄰域
    內(nèi)蒙古包頭市第一實驗小學六(3)中隊
    少先隊活動(2022年5期)2022-06-06 03:45:16
    稀疏圖平方圖的染色數(shù)上界
    基于鄰域競賽的多目標優(yōu)化算法
    自動化學報(2018年7期)2018-08-20 02:59:04
    基于元胞自動機下的交通事故路段仿真
    智富時代(2018年5期)2018-07-18 17:52:04
    關于-型鄰域空間
    包頭市東河區(qū)城市綠化現(xiàn)狀評價
    基于元胞數(shù)據(jù)的多維數(shù)據(jù)傳遞機制
    北京測繪(2016年2期)2016-01-24 02:28:28
    包頭市“健康水工程”建設實踐與思考
    中國水利(2015年4期)2015-02-28 15:12:22
    基于AIS的航道移動瓶頸元胞自動機模型
    中國航海(2014年1期)2014-05-09 07:54:25
    基于時序擴展的鄰域保持嵌入算法及其在故障檢測中的應用
    美女国产视频在线观看| 久久久久人妻精品一区果冻| 久久精品久久精品一区二区三区| 97人妻天天添夜夜摸| 成人二区视频| 三上悠亚av全集在线观看| 久久久久久人人人人人| 男人舔女人的私密视频| 伊人久久国产一区二区| 丁香六月天网| 国产精品不卡视频一区二区| 人成视频在线观看免费观看| 色婷婷av一区二区三区视频| 精品一区二区三区四区五区乱码 | 亚洲精品久久久久久婷婷小说| 国产男女超爽视频在线观看| 黑丝袜美女国产一区| www.av在线官网国产| 亚洲成av片中文字幕在线观看 | 亚洲成色77777| 亚洲色图 男人天堂 中文字幕 | tube8黄色片| 久久久久久久久久久免费av| 最新中文字幕久久久久| 精品少妇内射三级| 亚洲精品一区蜜桃| 国产熟女欧美一区二区| 少妇的丰满在线观看| 99久久综合免费| 18禁动态无遮挡网站| 国产毛片在线视频| av免费观看日本| 美女脱内裤让男人舔精品视频| 在线观看免费日韩欧美大片| 亚洲精品久久午夜乱码| 下体分泌物呈黄色| a级片在线免费高清观看视频| 欧美精品人与动牲交sv欧美| 777米奇影视久久| 又黄又粗又硬又大视频| 国产色婷婷99| 久久久亚洲精品成人影院| 少妇的逼好多水| 99久国产av精品国产电影| 精品亚洲成a人片在线观看| videosex国产| 日本vs欧美在线观看视频| 国产成人免费无遮挡视频| 999精品在线视频| 国产伦理片在线播放av一区| 一区在线观看完整版| 欧美日韩视频高清一区二区三区二| 少妇人妻精品综合一区二区| 亚洲性久久影院| 国产黄色免费在线视频| 校园人妻丝袜中文字幕| 亚洲av在线观看美女高潮| 成人免费观看视频高清| 成人国语在线视频| a 毛片基地| 男女无遮挡免费网站观看| 男女下面插进去视频免费观看 | 午夜福利,免费看| 2022亚洲国产成人精品| 亚洲精品久久午夜乱码| 亚洲精品国产av蜜桃| 欧美日本中文国产一区发布| 桃花免费在线播放| 欧美日韩国产mv在线观看视频| 亚洲av在线观看美女高潮| 国产男女超爽视频在线观看| 国产男人的电影天堂91| 亚洲五月色婷婷综合| 少妇的逼水好多| 22中文网久久字幕| 久久午夜综合久久蜜桃| 亚洲精品国产av蜜桃| 久久久国产精品麻豆| av免费观看日本| 久久 成人 亚洲| 久久人人97超碰香蕉20202| 五月开心婷婷网| 搡女人真爽免费视频火全软件| 亚洲高清免费不卡视频| 熟妇人妻不卡中文字幕| 欧美xxⅹ黑人| 国产亚洲欧美精品永久| 久久99热6这里只有精品| 日韩成人伦理影院| 免费观看a级毛片全部| 日本av手机在线免费观看| 另类精品久久| 免费看光身美女| 最黄视频免费看| 91精品国产国语对白视频| 自线自在国产av| 亚洲av电影在线进入| 国产麻豆69| 视频中文字幕在线观看| 亚洲一区二区三区欧美精品| 亚洲,欧美,日韩| 亚洲精品美女久久av网站| 黄色一级大片看看| 男的添女的下面高潮视频| 欧美日韩精品成人综合77777| 2018国产大陆天天弄谢| 午夜免费鲁丝| 国产69精品久久久久777片| 亚洲三级黄色毛片| 免费高清在线观看日韩| 久久99蜜桃精品久久| 免费不卡的大黄色大毛片视频在线观看| 国产成人av激情在线播放| a 毛片基地| 久久久国产欧美日韩av| 看十八女毛片水多多多| freevideosex欧美| av片东京热男人的天堂| 中国国产av一级| 黄色毛片三级朝国网站| 97人妻天天添夜夜摸| 晚上一个人看的免费电影| 日本91视频免费播放| 赤兔流量卡办理| 99精国产麻豆久久婷婷| 一级,二级,三级黄色视频| 哪个播放器可以免费观看大片| 亚洲精品乱久久久久久| 黄色怎么调成土黄色| 国产亚洲精品久久久com| 婷婷色av中文字幕| 99国产综合亚洲精品| av又黄又爽大尺度在线免费看| 9热在线视频观看99| 欧美老熟妇乱子伦牲交| 少妇被粗大猛烈的视频| 久久久久精品性色| 国产成人一区二区在线| 十分钟在线观看高清视频www| 男人爽女人下面视频在线观看| 赤兔流量卡办理| 日本色播在线视频| 老司机亚洲免费影院| 在线观看免费视频网站a站| 汤姆久久久久久久影院中文字幕| av国产精品久久久久影院| 久久精品国产a三级三级三级| 欧美日韩一区二区视频在线观看视频在线| 一级a做视频免费观看| 三级国产精品片| 国产在线免费精品| 日韩大片免费观看网站| 另类精品久久| 欧美xxⅹ黑人| av国产精品久久久久影院| 丁香六月天网| 国产精品成人在线| 又黄又粗又硬又大视频| 国产在线视频一区二区| 国产国拍精品亚洲av在线观看| 久久久久国产精品人妻一区二区| 黑丝袜美女国产一区| 国产激情久久老熟女| 国产69精品久久久久777片| 国产免费现黄频在线看| 日韩av不卡免费在线播放| 久久久国产一区二区| 亚洲人与动物交配视频| 国产白丝娇喘喷水9色精品| 成人国语在线视频| 飞空精品影院首页| 国产成人精品无人区| 精品第一国产精品| 一级爰片在线观看| 日产精品乱码卡一卡2卡三| 亚洲 欧美一区二区三区| 十八禁高潮呻吟视频| 欧美精品一区二区大全| 亚洲欧美日韩卡通动漫| 欧美精品高潮呻吟av久久| 丝袜人妻中文字幕| av卡一久久| 国产成人精品一,二区| 91午夜精品亚洲一区二区三区| 最近中文字幕高清免费大全6| 晚上一个人看的免费电影| 国产精品国产av在线观看| 欧美精品一区二区免费开放| www.熟女人妻精品国产 | 肉色欧美久久久久久久蜜桃| 精品少妇久久久久久888优播| 涩涩av久久男人的天堂| 男人操女人黄网站| 最近手机中文字幕大全| 深夜精品福利| 久久久国产一区二区| 蜜桃国产av成人99| 久久久亚洲精品成人影院| 久久国产亚洲av麻豆专区| 天天影视国产精品| 大话2 男鬼变身卡| av又黄又爽大尺度在线免费看| 国产精品熟女久久久久浪| 日韩欧美一区视频在线观看| 国产成人精品福利久久| 国产午夜精品一二区理论片| 熟女电影av网| 18禁裸乳无遮挡动漫免费视频| av电影中文网址| 亚洲国产精品国产精品| 人人妻人人爽人人添夜夜欢视频| 欧美精品一区二区大全| 免费av中文字幕在线| 午夜老司机福利剧场| 久久久久久人人人人人| 欧美人与善性xxx| 视频在线观看一区二区三区| 高清av免费在线| 我要看黄色一级片免费的| 女人被躁到高潮嗷嗷叫费观| 亚洲精品一二三| 大香蕉久久成人网| 国产精品一区二区在线不卡| 久久国产精品男人的天堂亚洲 | 一二三四在线观看免费中文在 | 久久午夜综合久久蜜桃| 久久狼人影院| 美女国产视频在线观看| 飞空精品影院首页| 18在线观看网站| 好男人视频免费观看在线| 97超碰精品成人国产| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品成人av观看孕妇| 老司机影院毛片| av天堂久久9| tube8黄色片| 观看av在线不卡| 亚洲精品日本国产第一区| 九草在线视频观看| 黑人欧美特级aaaaaa片| 成年动漫av网址| 26uuu在线亚洲综合色| 狂野欧美激情性bbbbbb| 午夜福利在线观看免费完整高清在| 亚洲精品久久午夜乱码| 日韩成人伦理影院| 丰满乱子伦码专区| 国产亚洲av片在线观看秒播厂| 国产精品偷伦视频观看了| 老司机亚洲免费影院| 免费大片黄手机在线观看| 国产激情久久老熟女| 内地一区二区视频在线| 超色免费av| 99久久精品国产国产毛片| 久久久久国产网址| 欧美精品av麻豆av| 国产又爽黄色视频| 日日撸夜夜添| 免费黄网站久久成人精品| 日韩av不卡免费在线播放| 亚洲综合精品二区| 精品人妻熟女毛片av久久网站| 一级片'在线观看视频| 香蕉精品网在线| 人成视频在线观看免费观看| 2021少妇久久久久久久久久久| 欧美+日韩+精品| 九色成人免费人妻av| 中文字幕亚洲精品专区| 欧美亚洲 丝袜 人妻 在线| 99久久精品国产国产毛片| 大香蕉久久网| 黄色怎么调成土黄色| 亚洲精品久久成人aⅴ小说| 亚洲综合色惰| 久久精品国产亚洲av天美| 国产黄色免费在线视频| 国产黄色视频一区二区在线观看| 亚洲欧美精品自产自拍| 一级毛片我不卡| 成年人免费黄色播放视频| 国产亚洲午夜精品一区二区久久| 国产日韩欧美在线精品| 18禁裸乳无遮挡动漫免费视频| 亚洲国产成人一精品久久久| 2022亚洲国产成人精品| 侵犯人妻中文字幕一二三四区| 精品熟女少妇av免费看| 人妻人人澡人人爽人人| 亚洲成人av在线免费| 校园人妻丝袜中文字幕| 欧美激情国产日韩精品一区| 两个人看的免费小视频| 久久毛片免费看一区二区三区| 久久久久视频综合| av在线app专区| 美女主播在线视频| 久热久热在线精品观看| 波多野结衣一区麻豆| 欧美xxxx性猛交bbbb| 狠狠婷婷综合久久久久久88av| 久久97久久精品| 久久av网站| 国国产精品蜜臀av免费| 一区二区日韩欧美中文字幕 | 国产精品蜜桃在线观看| 日本色播在线视频| 久热久热在线精品观看| av免费观看日本| 建设人人有责人人尽责人人享有的| 久久99一区二区三区| av国产久精品久网站免费入址| 在线精品无人区一区二区三| 久久综合国产亚洲精品| 宅男免费午夜| 免费久久久久久久精品成人欧美视频 | 涩涩av久久男人的天堂| 观看av在线不卡| 国产精品99久久99久久久不卡 | 精品少妇黑人巨大在线播放| 久久热在线av| 亚洲国产精品999| 国产亚洲午夜精品一区二区久久| h视频一区二区三区| 国产精品人妻久久久影院| av黄色大香蕉| av片东京热男人的天堂| 乱码一卡2卡4卡精品| 亚洲av成人精品一二三区| 欧美成人精品欧美一级黄| 只有这里有精品99| 久久精品久久久久久久性| 这个男人来自地球电影免费观看 | 热re99久久国产66热| 在线精品无人区一区二区三| 亚洲精品久久午夜乱码| 国产欧美亚洲国产| 91精品国产国语对白视频| 欧美日本中文国产一区发布| 高清不卡的av网站| 五月开心婷婷网| 久久97久久精品| 十八禁网站网址无遮挡| 亚洲欧美日韩卡通动漫| videosex国产| 9191精品国产免费久久| 在线天堂最新版资源| 精品卡一卡二卡四卡免费| 黄色配什么色好看| 男男h啪啪无遮挡| 国产一区二区三区综合在线观看 | 美女脱内裤让男人舔精品视频| 91精品伊人久久大香线蕉| 看免费成人av毛片| 久久国产亚洲av麻豆专区| 免费大片黄手机在线观看| xxxhd国产人妻xxx| 国产成人精品福利久久| 亚洲婷婷狠狠爱综合网| 在线天堂最新版资源| 一区二区三区乱码不卡18| 国产精品一区二区在线观看99| 少妇猛男粗大的猛烈进出视频| 日韩,欧美,国产一区二区三区| 日本av手机在线免费观看| 91精品三级在线观看| 色视频在线一区二区三区| 欧美激情国产日韩精品一区| av国产精品久久久久影院| 在线精品无人区一区二区三| 亚洲精品av麻豆狂野| 欧美激情极品国产一区二区三区 | 中文字幕人妻丝袜制服| 国产精品久久久久成人av| 免费女性裸体啪啪无遮挡网站| 春色校园在线视频观看| 最近最新中文字幕大全免费视频 | 母亲3免费完整高清在线观看 | 青春草国产在线视频| 汤姆久久久久久久影院中文字幕| 婷婷色av中文字幕| 免费播放大片免费观看视频在线观看| 精品亚洲乱码少妇综合久久| 色吧在线观看| 精品视频人人做人人爽| 久久ye,这里只有精品| 国产又色又爽无遮挡免| 黄片播放在线免费| 亚洲国产日韩一区二区| 亚洲高清免费不卡视频| 免费观看在线日韩| 欧美国产精品va在线观看不卡| 最近最新中文字幕大全免费视频 | 国产黄频视频在线观看| 国产无遮挡羞羞视频在线观看| 日本黄大片高清| 激情五月婷婷亚洲| 有码 亚洲区| 婷婷色av中文字幕| 男女下面插进去视频免费观看 | 丁香六月天网| 亚洲激情五月婷婷啪啪| 亚洲精品av麻豆狂野| 日韩欧美一区视频在线观看| 69精品国产乱码久久久| 纯流量卡能插随身wifi吗| 亚洲美女视频黄频| 国产精品麻豆人妻色哟哟久久| 中国三级夫妇交换| 午夜免费观看性视频| 狂野欧美激情性xxxx在线观看| 午夜福利,免费看| 亚洲精品久久午夜乱码| 欧美最新免费一区二区三区| 老司机亚洲免费影院| 久久久久久久久久久久大奶| 天堂俺去俺来也www色官网| 日本色播在线视频| 伊人久久国产一区二区| 精品久久蜜臀av无| 国产午夜精品一二区理论片| 亚洲精品久久午夜乱码| 久久久久久久久久人人人人人人| 大香蕉97超碰在线| a级毛片在线看网站| 99热这里只有是精品在线观看| 亚洲天堂av无毛| 韩国精品一区二区三区 | 久热这里只有精品99| 激情五月婷婷亚洲| 日韩,欧美,国产一区二区三区| 99久久综合免费| 精品熟女少妇av免费看| 久久久久国产精品人妻一区二区| 亚洲天堂av无毛| 亚洲av在线观看美女高潮| 黑人高潮一二区| 久久精品aⅴ一区二区三区四区 | 国产免费一区二区三区四区乱码| 大话2 男鬼变身卡| 青春草国产在线视频| av福利片在线| 午夜福利视频在线观看免费| 国产精品三级大全| 亚洲久久久国产精品| 一级毛片 在线播放| 九色亚洲精品在线播放| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品av麻豆狂野| 中国美白少妇内射xxxbb| √禁漫天堂资源中文www| 又黄又爽又刺激的免费视频.| 中文字幕人妻熟女乱码| 国产乱来视频区| 精品福利永久在线观看| 我要看黄色一级片免费的| 亚洲成色77777| 久久精品国产综合久久久 | 韩国高清视频一区二区三区| 亚洲欧洲国产日韩| 成人综合一区亚洲| 青青草视频在线视频观看| 久久久久久久久久久免费av| 女性被躁到高潮视频| 国产极品天堂在线| 丝瓜视频免费看黄片| av又黄又爽大尺度在线免费看| 亚洲人成网站在线观看播放| 高清不卡的av网站| 男女免费视频国产| 人妻系列 视频| 免费在线观看黄色视频的| 久久久国产欧美日韩av| 国产亚洲最大av| 免费看光身美女| 一区二区三区乱码不卡18| 美女大奶头黄色视频| 欧美xxxx性猛交bbbb| 精品福利永久在线观看| 精品少妇内射三级| 九色亚洲精品在线播放| 日本-黄色视频高清免费观看| 亚洲,一卡二卡三卡| 色5月婷婷丁香| 国产视频首页在线观看| 日韩欧美一区视频在线观看| 777米奇影视久久| 少妇精品久久久久久久| 免费在线观看黄色视频的| 99国产综合亚洲精品| 亚洲精品av麻豆狂野| 精品人妻熟女毛片av久久网站| 观看美女的网站| kizo精华| 色哟哟·www| 精品国产一区二区三区四区第35| 免费不卡的大黄色大毛片视频在线观看| 久久精品国产a三级三级三级| 精品国产乱码久久久久久小说| 欧美亚洲日本最大视频资源| 成人手机av| 性色av一级| 91国产中文字幕| 国产亚洲一区二区精品| 视频在线观看一区二区三区| 99热6这里只有精品| 婷婷色麻豆天堂久久| 亚洲五月色婷婷综合| 亚洲精品日韩在线中文字幕| 97精品久久久久久久久久精品| 亚洲激情五月婷婷啪啪| 久久久久国产精品人妻一区二区| 亚洲av中文av极速乱| 亚洲精品成人av观看孕妇| 国产黄色免费在线视频| 一级黄片播放器| 哪个播放器可以免费观看大片| tube8黄色片| 两个人看的免费小视频| 亚洲中文av在线| 亚洲国产精品成人久久小说| 69精品国产乱码久久久| 精品国产一区二区三区四区第35| 最近中文字幕高清免费大全6| 香蕉丝袜av| 黄网站色视频无遮挡免费观看| 亚洲av成人精品一二三区| 久久这里有精品视频免费| 日日摸夜夜添夜夜爱| 边亲边吃奶的免费视频| 亚洲精品国产av蜜桃| 美女国产视频在线观看| 免费少妇av软件| 一边亲一边摸免费视频| 丝袜在线中文字幕| 男女啪啪激烈高潮av片| 少妇的逼好多水| 亚洲丝袜综合中文字幕| 夫妻午夜视频| 宅男免费午夜| 成人18禁高潮啪啪吃奶动态图| 伦精品一区二区三区| 亚洲av电影在线观看一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一级爰片在线观看| av一本久久久久| 免费看光身美女| 丰满少妇做爰视频| 男女下面插进去视频免费观看 | 蜜臀久久99精品久久宅男| 国产精品国产三级国产专区5o| 久久99蜜桃精品久久| 女的被弄到高潮叫床怎么办| 少妇被粗大猛烈的视频| 精品亚洲乱码少妇综合久久| 日日撸夜夜添| 久久国产精品男人的天堂亚洲 | 久热这里只有精品99| 久久久久人妻精品一区果冻| 99国产精品免费福利视频| 日日摸夜夜添夜夜爱| 亚洲成人av在线免费| 亚洲av免费高清在线观看| 男的添女的下面高潮视频| 韩国av在线不卡| 欧美日韩精品成人综合77777| 黑人巨大精品欧美一区二区蜜桃 | 久久久久久久亚洲中文字幕| 一边亲一边摸免费视频| 97超碰精品成人国产| 看免费av毛片| 久久久国产精品麻豆| 咕卡用的链子| 纵有疾风起免费观看全集完整版| 99热全是精品| 婷婷色综合大香蕉| av黄色大香蕉| xxx大片免费视频| 久久久欧美国产精品| 久久99热这里只频精品6学生| 亚洲av电影在线进入| 在现免费观看毛片| 午夜日本视频在线| 欧美日韩亚洲高清精品| av片东京热男人的天堂| 成年女人在线观看亚洲视频| av一本久久久久| 久久鲁丝午夜福利片| 国产精品久久久久久av不卡| 永久网站在线| 欧美成人午夜精品| 国产精品久久久久久久电影| 国产成人aa在线观看| 亚洲精品,欧美精品| 久久人人97超碰香蕉20202| 久久久久视频综合| 日韩欧美一区视频在线观看| 国产成人av激情在线播放| 久久女婷五月综合色啪小说| 老司机亚洲免费影院| av在线观看视频网站免费| 久久久久久久国产电影| 各种免费的搞黄视频| 成人午夜精彩视频在线观看| 欧美日韩视频高清一区二区三区二| 国产精品.久久久| 国内精品宾馆在线| 美女大奶头黄色视频| 欧美亚洲日本最大视频资源| 菩萨蛮人人尽说江南好唐韦庄| 欧美97在线视频| 黑人欧美特级aaaaaa片| 夫妻性生交免费视频一级片|