• <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
    基于時序擴展的鄰域保持嵌入算法及其在故障檢測中的應用
    亚洲激情在线av| 亚洲精品中文字幕一二三四区| 一区福利在线观看| 99热只有精品国产| 亚洲欧美日韩卡通动漫| 国产亚洲av嫩草精品影院| 成人性生交大片免费视频hd| www国产在线视频色| 日韩欧美国产在线观看| 成人欧美大片| 色老头精品视频在线观看| 麻豆久久精品国产亚洲av| 久久国产乱子伦精品免费另类| 午夜福利欧美成人| 床上黄色一级片| 嫩草影院精品99| av视频在线观看入口| 亚洲国产欧美一区二区综合| 久久天躁狠狠躁夜夜2o2o| 亚洲五月婷婷丁香| 高清毛片免费观看视频网站| 国产高清三级在线| 亚洲中文字幕一区二区三区有码在线看 | 男人舔女人的私密视频| 视频区欧美日本亚洲| 精品欧美国产一区二区三| 日韩欧美国产在线观看| 我的老师免费观看完整版| 亚洲av美国av| 母亲3免费完整高清在线观看| 女人被狂操c到高潮| 亚洲国产欧美网| 身体一侧抽搐| 国产av不卡久久| 给我免费播放毛片高清在线观看| 淫妇啪啪啪对白视频| 亚洲精品一卡2卡三卡4卡5卡| 婷婷精品国产亚洲av在线| 国产精品九九99| 又大又爽又粗| 午夜久久久久精精品| h日本视频在线播放| 亚洲片人在线观看| 性色avwww在线观看| 亚洲av电影不卡..在线观看| 国产精品99久久99久久久不卡| www日本黄色视频网| 国产精品爽爽va在线观看网站| 成人性生交大片免费视频hd| 亚洲国产欧美网| 91字幕亚洲| 亚洲精品粉嫩美女一区| 午夜精品久久久久久毛片777| 久久久久久久久免费视频了| 悠悠久久av| 五月玫瑰六月丁香| 操出白浆在线播放| 91在线观看av| 无遮挡黄片免费观看| 亚洲av电影在线进入| 久99久视频精品免费| 怎么达到女性高潮| 国产爱豆传媒在线观看| 黄色视频,在线免费观看| 亚洲人与动物交配视频| 国产熟女xx| 国产精品久久久久久亚洲av鲁大| 91av网站免费观看| 色综合欧美亚洲国产小说| 又黄又粗又硬又大视频| www日本黄色视频网| 两个人视频免费观看高清| 精品久久蜜臀av无| 精品不卡国产一区二区三区| 热99在线观看视频| 午夜福利视频1000在线观看| 午夜精品一区二区三区免费看| 精品国产超薄肉色丝袜足j| 久久中文看片网| 丝袜人妻中文字幕| 国产av麻豆久久久久久久| 国产一区二区三区视频了| 88av欧美| 精品一区二区三区四区五区乱码| 国产午夜精品论理片| 日韩欧美精品v在线| 在线观看美女被高潮喷水网站 | 日本精品一区二区三区蜜桃| 欧美日韩福利视频一区二区| 9191精品国产免费久久| 亚洲国产精品999在线| 99国产综合亚洲精品| xxxwww97欧美| 日韩国内少妇激情av| 亚洲国产精品999在线| 成人午夜高清在线视频| 国产综合懂色| 日本免费一区二区三区高清不卡| 制服丝袜大香蕉在线| 欧美丝袜亚洲另类 | 成人午夜高清在线视频| 亚洲av成人精品一区久久| 国产午夜福利久久久久久| 久久久久久久久久黄片| 午夜激情福利司机影院| 色老头精品视频在线观看| 国产主播在线观看一区二区| 久久性视频一级片| 久久久精品大字幕| av中文乱码字幕在线| 偷拍熟女少妇极品色| 国产97色在线日韩免费| 国产亚洲精品一区二区www| 国产真人三级小视频在线观看| av视频在线观看入口| 一卡2卡三卡四卡精品乱码亚洲| 岛国在线观看网站| 丁香欧美五月| 欧美xxxx黑人xx丫x性爽| 精品久久蜜臀av无| 久久精品91蜜桃| 亚洲,欧美精品.| 成人午夜高清在线视频| 国产高清三级在线| 嫩草影视91久久| a级毛片a级免费在线| 后天国语完整版免费观看| 日日夜夜操网爽| 999久久久国产精品视频| a级毛片a级免费在线| 国产69精品久久久久777片 | 非洲黑人性xxxx精品又粗又长| 一级毛片女人18水好多| 亚洲 国产 在线| 黄色女人牲交| 久久久国产欧美日韩av| 亚洲真实伦在线观看| 亚洲精品色激情综合| 欧美成人一区二区免费高清观看 | 成年版毛片免费区| ponron亚洲| 亚洲中文av在线| 99热6这里只有精品| 国产熟女xx| 久久欧美精品欧美久久欧美| 国产精品久久久人人做人人爽| 一个人看的www免费观看视频| 1000部很黄的大片| 在线观看免费视频日本深夜| 国产成人系列免费观看| or卡值多少钱| 久久九九热精品免费| 99久久成人亚洲精品观看| 国产毛片a区久久久久| 老汉色∧v一级毛片| 宅男免费午夜| 一本久久中文字幕| 色视频www国产| 法律面前人人平等表现在哪些方面| 欧美日韩精品网址| 亚洲精品美女久久av网站| 黄色片一级片一级黄色片| 亚洲av片天天在线观看| 亚洲美女视频黄频| 国产麻豆成人av免费视频| 不卡av一区二区三区| 国产精品久久久人人做人人爽| 精品一区二区三区视频在线 | 国产精品98久久久久久宅男小说| 天堂√8在线中文| 日本撒尿小便嘘嘘汇集6| 亚洲va日本ⅴa欧美va伊人久久| 中文字幕人成人乱码亚洲影| 99在线视频只有这里精品首页| 亚洲av成人一区二区三| 黄色女人牲交| 国内精品久久久久久久电影| 天天一区二区日本电影三级| 床上黄色一级片| www国产在线视频色| 小说图片视频综合网站| 亚洲性夜色夜夜综合| 久久久久国产一级毛片高清牌| 三级男女做爰猛烈吃奶摸视频| 亚洲精品美女久久久久99蜜臀| 亚洲欧美激情综合另类| 一二三四在线观看免费中文在| 床上黄色一级片| 日韩免费av在线播放| 琪琪午夜伦伦电影理论片6080| 久久草成人影院| 日本黄色视频三级网站网址| 亚洲av电影在线进入| 国产久久久一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 国产精品永久免费网站| 亚洲真实伦在线观看| 国模一区二区三区四区视频 | 日本a在线网址| 国产高清videossex| 久久欧美精品欧美久久欧美| 18美女黄网站色大片免费观看| 99热精品在线国产| 色综合亚洲欧美另类图片| 欧美日韩一级在线毛片| 久久国产精品影院| 亚洲狠狠婷婷综合久久图片| 丰满人妻熟妇乱又伦精品不卡| 日韩人妻高清精品专区| www.www免费av| 少妇的逼水好多| 草草在线视频免费看| 精品一区二区三区视频在线 | 久久婷婷人人爽人人干人人爱| 精品欧美国产一区二区三| 老司机深夜福利视频在线观看| 亚洲激情在线av| 99久久精品国产亚洲精品| 婷婷精品国产亚洲av| 亚洲精品粉嫩美女一区| 99久久成人亚洲精品观看| 国产精品影院久久| 99在线视频只有这里精品首页| 欧美日韩瑟瑟在线播放| 中文字幕熟女人妻在线| 精品国内亚洲2022精品成人| 热99在线观看视频| 亚洲成a人片在线一区二区| 成人特级黄色片久久久久久久| 嫩草影院精品99| 成年免费大片在线观看| 亚洲成av人片在线播放无| 在线永久观看黄色视频| 男插女下体视频免费在线播放| xxx96com| 91av网一区二区| 久久午夜亚洲精品久久| 在线a可以看的网站| 舔av片在线| 久久人妻av系列| 国产三级在线视频| 国产亚洲av嫩草精品影院| 国产精品亚洲美女久久久| 91麻豆精品激情在线观看国产| 黄色片一级片一级黄色片| 黄色日韩在线| а√天堂www在线а√下载| 亚洲在线观看片| 日韩 欧美 亚洲 中文字幕| 久9热在线精品视频| 一个人看视频在线观看www免费 | 无人区码免费观看不卡| 精品一区二区三区视频在线观看免费| 国产精品 国内视频| 中文字幕av在线有码专区| 日本 av在线| bbb黄色大片| 在线观看免费视频日本深夜| 亚洲aⅴ乱码一区二区在线播放| 99久久精品一区二区三区| 国产蜜桃级精品一区二区三区| 亚洲欧美激情综合另类| 国产极品精品免费视频能看的| 观看美女的网站| 精品国产乱子伦一区二区三区| 成人三级黄色视频| 99久久国产精品久久久| 1000部很黄的大片| 亚洲成av人片免费观看| 国产精品99久久99久久久不卡| 欧美性猛交╳xxx乱大交人| 久久久久久久久免费视频了| 热99在线观看视频| 色av中文字幕| 日本与韩国留学比较| 久久午夜亚洲精品久久| 国产精品久久久久久久电影 | 午夜福利成人在线免费观看| 无遮挡黄片免费观看| 亚洲人成电影免费在线| 999久久久精品免费观看国产| 亚洲欧美一区二区三区黑人| 人人妻人人看人人澡| 国产91精品成人一区二区三区| 国产精品亚洲av一区麻豆| 日本免费一区二区三区高清不卡| 亚洲色图 男人天堂 中文字幕| 国产精品国产高清国产av| 舔av片在线| 精品国内亚洲2022精品成人| 欧美黑人欧美精品刺激| 日韩欧美在线二视频| 亚洲自拍偷在线| 国产欧美日韩一区二区精品| 精品国产亚洲在线| 嫁个100分男人电影在线观看| 国产免费男女视频| 综合色av麻豆| 亚洲 欧美一区二区三区| 国产主播在线观看一区二区| 老司机福利观看| 久久久久亚洲av毛片大全| 久久九九热精品免费| 欧美日韩精品网址| 国产成人精品久久二区二区91| 亚洲av中文字字幕乱码综合| 国产精品一区二区三区四区久久| 国产精品一区二区免费欧美| 久久精品夜夜夜夜夜久久蜜豆| 成人鲁丝片一二三区免费| 大型黄色视频在线免费观看| 国产欧美日韩一区二区精品| 欧美日韩黄片免| 国产精品 国内视频| 亚洲最大成人中文| 舔av片在线| 99久久精品一区二区三区| 深夜精品福利| 欧美大码av| 国产主播在线观看一区二区| 亚洲avbb在线观看| 巨乳人妻的诱惑在线观看| 亚洲人成伊人成综合网2020| 我的老师免费观看完整版| 免费看十八禁软件| 欧美日韩精品网址| 夜夜爽天天搞| 在线十欧美十亚洲十日本专区| 亚洲av电影在线进入| 免费在线观看日本一区| 日韩欧美三级三区| 九色国产91popny在线| 一a级毛片在线观看| 首页视频小说图片口味搜索| 日本五十路高清| 亚洲成人久久爱视频| 欧美日韩黄片免| 性色avwww在线观看| 网址你懂的国产日韩在线| 观看美女的网站| netflix在线观看网站| 久久久久久久午夜电影| 成人精品一区二区免费| 亚洲av电影在线进入| 国产激情偷乱视频一区二区| 亚洲中文av在线| 免费观看精品视频网站| 99在线人妻在线中文字幕| xxx96com| 亚洲18禁久久av| 国内久久婷婷六月综合欲色啪| 在线观看一区二区三区| 特大巨黑吊av在线直播| 国产精品久久久av美女十八| 老司机午夜福利在线观看视频| 亚洲专区中文字幕在线| 高潮久久久久久久久久久不卡| 一进一出好大好爽视频| 久久久久久久久免费视频了| 国产av不卡久久| 免费在线观看日本一区| 亚洲熟妇熟女久久| 99精品在免费线老司机午夜| 19禁男女啪啪无遮挡网站| 亚洲 欧美一区二区三区| 亚洲国产欧洲综合997久久,| 午夜精品一区二区三区免费看| 在线免费观看的www视频| 国产综合懂色| 一二三四社区在线视频社区8| 最新美女视频免费是黄的| 久久天堂一区二区三区四区| 九色国产91popny在线| 亚洲激情在线av| 高清在线国产一区| 久久久成人免费电影| 成熟少妇高潮喷水视频| 精品国产三级普通话版| 在线播放国产精品三级| 亚洲av免费在线观看| 精品久久久久久久久久免费视频| 亚洲国产欧洲综合997久久,| 99热精品在线国产| 亚洲熟妇熟女久久| 免费一级毛片在线播放高清视频| 99在线人妻在线中文字幕| 亚洲欧美日韩高清在线视频| 亚洲精品美女久久久久99蜜臀| 男女那种视频在线观看| 深夜精品福利| 美女被艹到高潮喷水动态| 久久久久久久午夜电影| 国产亚洲精品一区二区www| 免费大片18禁| 两性午夜刺激爽爽歪歪视频在线观看| 成人午夜高清在线视频| 婷婷精品国产亚洲av| 最近最新免费中文字幕在线| 一二三四在线观看免费中文在| 亚洲精华国产精华精| 免费av不卡在线播放| 精品久久久久久久末码| 国产极品精品免费视频能看的| 久9热在线精品视频| 99久久久亚洲精品蜜臀av| 亚洲欧美日韩高清专用| 91字幕亚洲| 在线免费观看不下载黄p国产 | 久久久久国产一级毛片高清牌| 少妇人妻一区二区三区视频| 免费看日本二区| 在线视频色国产色| 国产精品久久电影中文字幕| 色吧在线观看| 99热精品在线国产| 成在线人永久免费视频| 无限看片的www在线观看| 搡老妇女老女人老熟妇| 国产黄片美女视频| 男插女下体视频免费在线播放| 国产午夜福利久久久久久| av女优亚洲男人天堂 | 久久久久国产一级毛片高清牌| 啦啦啦免费观看视频1| 欧洲精品卡2卡3卡4卡5卡区| 国产精品久久久久久精品电影| 色尼玛亚洲综合影院| 国产精品久久久人人做人人爽| 亚洲天堂国产精品一区在线| 长腿黑丝高跟| 亚洲成av人片在线播放无| 偷拍熟女少妇极品色| 少妇熟女aⅴ在线视频| 亚洲七黄色美女视频| 亚洲最大成人中文| 久久亚洲真实| 熟女电影av网| 亚洲成av人片免费观看| 久久精品人妻少妇| 中出人妻视频一区二区| 色吧在线观看| 波多野结衣巨乳人妻| 老司机午夜福利在线观看视频| 亚洲欧美日韩无卡精品| 色综合欧美亚洲国产小说| 久久久成人免费电影| 久久久久久久久免费视频了| 综合色av麻豆| 少妇熟女aⅴ在线视频| 久久精品91无色码中文字幕| 亚洲一区二区三区色噜噜| 国产淫片久久久久久久久 | 真人做人爱边吃奶动态| 精品免费久久久久久久清纯| 国产精品一区二区三区四区免费观看 | 在线观看日韩欧美| 日本五十路高清| 校园春色视频在线观看| 国产成年人精品一区二区| 国产av麻豆久久久久久久| 俄罗斯特黄特色一大片| 在线观看66精品国产| 亚洲欧美日韩高清专用| 国产三级在线视频| 欧美国产日韩亚洲一区| 国产一区二区在线av高清观看| 亚洲av成人一区二区三| 国产精品久久久久久人妻精品电影| 在线观看66精品国产| 久久久国产成人免费| 精品人妻1区二区| 人妻久久中文字幕网| 久久精品国产亚洲av香蕉五月| 亚洲欧美日韩高清在线视频| 极品教师在线免费播放| 激情在线观看视频在线高清| 男人和女人高潮做爰伦理| 午夜免费成人在线视频| 色av中文字幕| 中文字幕人妻丝袜一区二区| 欧美日韩国产亚洲二区| 精品免费久久久久久久清纯| 国产精品久久久久久久电影 | 亚洲中文日韩欧美视频| 两个人视频免费观看高清| 精品不卡国产一区二区三区| 久99久视频精品免费| 久久久国产欧美日韩av| www.熟女人妻精品国产| 十八禁网站免费在线| 色综合婷婷激情| 日本a在线网址| 中文亚洲av片在线观看爽| 亚洲国产色片| 欧美av亚洲av综合av国产av| 级片在线观看| a级毛片a级免费在线| 国内精品一区二区在线观看| 久久久久九九精品影院| 伊人久久大香线蕉亚洲五| 久久久久国产一级毛片高清牌| 欧美午夜高清在线| 美女被艹到高潮喷水动态| 免费大片18禁| 国产视频一区二区在线看| 99在线视频只有这里精品首页| 免费看a级黄色片| 国产激情久久老熟女| 天堂动漫精品| 欧美乱色亚洲激情| 日本一二三区视频观看| 久久精品人妻少妇| 国产激情久久老熟女| 波多野结衣高清作品| 色综合婷婷激情| 免费人成视频x8x8入口观看| 亚洲精品在线美女| 美女高潮的动态| 91麻豆精品激情在线观看国产| 国产1区2区3区精品| 亚洲七黄色美女视频| 国产极品精品免费视频能看的| 97超级碰碰碰精品色视频在线观看| 国内揄拍国产精品人妻在线| 国产欧美日韩精品一区二区| 99国产精品一区二区蜜桃av| 成人av一区二区三区在线看| 中出人妻视频一区二区| 99精品欧美一区二区三区四区| 在线观看舔阴道视频| 亚洲国产欧美一区二区综合| 伊人久久大香线蕉亚洲五| www.精华液| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av五月六月丁香网| 亚洲国产色片| 99re在线观看精品视频| 人妻丰满熟妇av一区二区三区| 久久中文字幕人妻熟女| 亚洲人成网站在线播放欧美日韩| 在线看三级毛片| 1024香蕉在线观看| 波多野结衣巨乳人妻| 亚洲国产欧美一区二区综合| 欧美+亚洲+日韩+国产| 久久伊人香网站| 99在线人妻在线中文字幕| 国产精品乱码一区二三区的特点| 一本一本综合久久| 搡老岳熟女国产| 观看免费一级毛片| 免费观看人在逋| 动漫黄色视频在线观看| 成年免费大片在线观看| 国产成年人精品一区二区| 国产精品亚洲美女久久久| 男人的好看免费观看在线视频| www.自偷自拍.com| av黄色大香蕉| 亚洲av成人av| 国产精品亚洲美女久久久| 久久午夜亚洲精品久久| 欧美在线黄色| 脱女人内裤的视频| 日韩av在线大香蕉| 午夜福利在线观看免费完整高清在 | 中文字幕人妻丝袜一区二区| 国产精品乱码一区二三区的特点| 麻豆成人午夜福利视频| 久久精品aⅴ一区二区三区四区| 一区二区三区国产精品乱码| 天天躁日日操中文字幕| 久久精品人妻少妇| 亚洲欧洲精品一区二区精品久久久| 久久久水蜜桃国产精品网| 色老头精品视频在线观看| 国产成人aa在线观看| 在线观看一区二区三区| 中文字幕最新亚洲高清| 国产高清视频在线观看网站| 国产成人精品久久二区二区免费| 亚洲av成人不卡在线观看播放网| 国内揄拍国产精品人妻在线| 在线看三级毛片| 免费看a级黄色片| 免费看日本二区| 一级a爱片免费观看的视频| 99久久综合精品五月天人人| 久久久成人免费电影| 国产aⅴ精品一区二区三区波| 亚洲美女视频黄频| 久久久久精品国产欧美久久久| 婷婷六月久久综合丁香| 久久久久精品国产欧美久久久| 国产三级黄色录像| av天堂中文字幕网| 欧美黄色淫秽网站| 日韩三级视频一区二区三区| 免费无遮挡裸体视频| 免费在线观看成人毛片| 99热只有精品国产| 婷婷亚洲欧美| 久久久久性生活片| 日韩欧美三级三区| 久久久久国产精品人妻aⅴ院| 热99re8久久精品国产| 黄片小视频在线播放| 超碰成人久久| 国产美女午夜福利| 日韩欧美三级三区| 国产精品国产高清国产av| 免费在线观看亚洲国产| 中文字幕人妻丝袜一区二区| 免费无遮挡裸体视频| e午夜精品久久久久久久|