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

    古交市礦區(qū)土地覆蓋時(shí)空動(dòng)態(tài)變化研究

    2016-01-18 02:53:53李宏宏
    山西焦煤科技 2015年1期
    關(guān)鍵詞:時(shí)空變化

    (E-mail)lihonghong26@163.com

    古交市礦區(qū)土地覆蓋時(shí)空動(dòng)態(tài)變化研究

    李宏宏

    (山西焦煤集團(tuán) 投資有限公司資源地質(zhì)部,山西太原030021)

    摘要為了研究礦區(qū)土地覆蓋的動(dòng)態(tài)變化規(guī)律,以古交礦區(qū)為研究區(qū),通過(guò)監(jiān)督分類(lèi)解譯獲取礦區(qū)1990年、2000年、2010年相同季節(jié)的土地覆蓋類(lèi)型圖。運(yùn)用轉(zhuǎn)移矩陣和CA-Markov模型方法分別對(duì)研究區(qū)土地覆蓋類(lèi)型動(dòng)態(tài)變化進(jìn)行了研究和預(yù)測(cè),結(jié)果表明:1) 1990—2010年林地的面積持續(xù)增加,年增長(zhǎng)率達(dá)到了5.49%,耕地和草地的面積大幅度減少,尤其是草地的面積,此外城鎮(zhèn)建設(shè)用地和水域輕微減少。2) 2010—2020年耕地和草地的面積進(jìn)一步減少,城鎮(zhèn)和工礦用地以及水域的面積小幅度增加,林地的面積繼續(xù)大幅度增長(zhǎng),年增長(zhǎng)率達(dá)到了2.18%。

    關(guān)鍵詞CA-MARKOV;時(shí)空變化;土地覆蓋類(lèi)型;古交礦區(qū)

    收稿日期:2014-11-23

    作者簡(jiǎn)介:李宏宏(1985—),男,山西代縣人,2010年畢業(yè)于山西農(nóng)業(yè)大學(xué),助理工程師,主要從事煤礦地質(zhì)測(cè)量與防治水等相關(guān)工作

    中圖分類(lèi)號(hào):TD167

    文獻(xiàn)標(biāo)識(shí)碼:B

    文章編號(hào):1672-0652(2015)01-0046-05

    AbstractIn order to study the dynamic change rule of land cover in mining area, takes Gujiao mining area as the study area. Through supervised classification and interpretation obtains the land cover type maps of mining area of the same season in 1990, 2000, 2010. By using the transfer matrix method and CA-Markov model to study and predict the dynamic change of land cover type, the results show that 1) from 1990 to 2010 the areas of forest land continue to increase, the annual growth rate is 5.49%, the areas of farmland and grassland reduce drastically, especially the grassland area, in addition, urban construction land and waters decrease slightly. 2) from 2010 to 2020, the areas of arable land and grassland decrease further. Towns, industrial and mining land and water areas increase small. Forest area continues to increase drastically, the annual growth rate reaches 2.18%.

    近年來(lái),因土地覆蓋/土地覆蓋變化 (LUCC) 導(dǎo)致的全球環(huán)境惡化的問(wèn)題日益突出,目前已成為國(guó)內(nèi)外相關(guān)組織和學(xué)者研究的首要內(nèi)容[1-7]。針對(duì)礦區(qū)的主要研究成果有胡召玲等用轉(zhuǎn)移矩陣對(duì)徐州煤礦區(qū)的土地覆蓋變化進(jìn)行的分析[8],李琳等對(duì)平朔露天礦區(qū)土地覆蓋/覆蓋變化的分析研究等[9]。山西省作為全國(guó)最主要的煤礦基地, 煤礦資源十分豐富,但近年來(lái)在為我國(guó)經(jīng)濟(jì)發(fā)展的增長(zhǎng)做出貢獻(xiàn)的同時(shí),由于過(guò)度開(kāi)采導(dǎo)致了山西省整體土地資源結(jié)構(gòu)不合理和生態(tài)環(huán)境狀況的惡化, 造成了礦區(qū)人地關(guān)系的緊張。本文以TM遙感影像數(shù)據(jù)為基礎(chǔ), 采用轉(zhuǎn)移矩陣定量分析的方法研究了古交礦區(qū) 1990—2010年土地覆蓋的動(dòng)態(tài)變化, 并通CA-Markov模型對(duì)礦區(qū)的土地覆蓋動(dòng)態(tài)變化進(jìn)行預(yù)測(cè),以期為古交市礦區(qū)的土地管理和生態(tài)環(huán)境的綜合整治與恢復(fù)提供決策參考和數(shù)據(jù)支持, 實(shí)現(xiàn)研究區(qū)土地的長(zhǎng)久可持續(xù)利用。

    1 研究區(qū)概括

    古交市作為全國(guó)最大的主焦煤生產(chǎn)基地,地處山西省太原市西北部,作為連接省城太原和山西省西北部的現(xiàn)代化工礦城市,已經(jīng)發(fā)展成為呂梁山東麓的交通樞紐和商品集散地。古交市地處北緯37°40′6〃至38°8′9〃,東經(jīng)111°43′8〃至112°21′5",東西寬50 km,南北長(zhǎng)53 km,總面積1 551 km2. 全市轄10個(gè)鄉(xiāng)鎮(zhèn)和4個(gè)街道辦事處,總?cè)丝?1萬(wàn)。

    礦區(qū)地形復(fù)雜,外圍地形較高,多在1 500 m以上。東部廟前山標(biāo)高1 865 m,西南部孤?tīng)斏綐?biāo)高2 202 m,區(qū)內(nèi)標(biāo)高一般多在1 000~1 350 m.礦區(qū)最低處寨上村汾河河谷,標(biāo)高958.40 m.礦區(qū)屬暖溫帶半干旱大陸性季風(fēng)氣候,雨量集中于7—9月,近年最大降雨量為632 mm.氣溫晝夜變化大,最高氣溫達(dá)40 ℃,最低為-20 ℃,自10月至翌年4月為結(jié)冰期,凍土深度為0.5~0.8 m.此外礦區(qū)交通發(fā)達(dá),太(原)古(交)鐵路從太原市汾河站至古交西曲站41 km,往西穿過(guò)礦區(qū)經(jīng)鎮(zhèn)城底至嵐縣,嵐縣經(jīng)方山至呂梁段正在建設(shè)中。從太原市區(qū)通往礦區(qū)的盤(pán)山公路長(zhǎng)56 km,古交至清徐也有公路相通。研究區(qū)概況見(jiàn)圖1.

    2數(shù)據(jù)源和土地覆蓋類(lèi)型提取

    本研究從地理空間數(shù)據(jù)云上獲取了 1990年、2000年和 2010 年3期landsat4-5的多光譜遙感影像,在ENVI5.0下將3期TM遙感影像以不同波段組合,并按照古交市的行政界線對(duì)TM影像進(jìn)行圖像裁剪;同時(shí)進(jìn)行幾何精校正、輻射定標(biāo)、大氣糾正等處理;最后進(jìn)行圖像拉伸處理,以增強(qiáng)影像判讀,利用監(jiān)督分類(lèi)的最大似然法,結(jié)合 GoogleEarth 工具進(jìn)行人機(jī)交互目視解譯,獲取研究區(qū)1990年、2000年和 2010 年3期的礦區(qū)土地覆蓋圖,見(jiàn)圖2,kappa系數(shù)達(dá)到了0.8以上,分類(lèi)精度超過(guò)了85%,參照研究區(qū)第二次土地調(diào)查圖,表明解譯結(jié)果精度良好。

    圖1 研究區(qū)概況圖

    圖2 1990年、2000年和 2010 年礦區(qū)土地覆蓋類(lèi)型圖

    3研究方法

    為了對(duì)研究區(qū)的土地覆蓋類(lèi)型動(dòng)態(tài)變化進(jìn)行研究,本研究以各地類(lèi)的面積為研究對(duì)象,進(jìn)而綜合反映區(qū)域土地覆蓋變化。面積變化能夠直觀地反應(yīng)不同土地覆蓋類(lèi)型的總量變化,進(jìn)而可掌握土地覆蓋變化總的規(guī)律和土地覆蓋結(jié)構(gòu)的變化態(tài)勢(shì)。首先根據(jù)1990年、2010年兩期土地覆蓋類(lèi)型圖,在ArcGis 10.0中運(yùn)用GIS的空間分析功能獲取研究區(qū)的轉(zhuǎn)移矩陣,并以此來(lái)分析研究區(qū)土地利用類(lèi)型動(dòng)態(tài)變化,獲取動(dòng)態(tài)變化趨勢(shì),然后利用CA-Markov模型對(duì)研究區(qū)的土地利用類(lèi)型進(jìn)行預(yù)測(cè),獲取2020年研究區(qū)土地利用類(lèi)型。

    馬爾科夫(Markov)模型俗稱(chēng)空間概率轉(zhuǎn)換模型,它是以馬爾科夫鏈為基礎(chǔ),研究某一事物或者類(lèi)型的狀態(tài)以及狀態(tài)之間轉(zhuǎn)移規(guī)律的一種隨機(jī)過(guò)程,其主要原理是通過(guò)對(duì)事件某一時(shí)刻T不同狀態(tài)的初始概率以及與T時(shí)刻狀態(tài)之間的轉(zhuǎn)移關(guān)系來(lái)確定T時(shí)刻狀態(tài)變化趨勢(shì)的方法[10]。元胞自動(dòng)機(jī)(Cellular Automata),一般簡(jiǎn)稱(chēng)CA是一種具有時(shí)空計(jì)算特征的動(dòng)力學(xué)模型,它是由單元、狀態(tài)、領(lǐng)域及轉(zhuǎn)換4部分組成,主要的特點(diǎn)是時(shí)間、空間、狀態(tài)都非常離散。分布在規(guī)則格網(wǎng)中的任何一個(gè)單元均取有限空間的離散狀態(tài),遵循同樣的作用規(guī)則,并依據(jù)確定的局部規(guī)則來(lái)做出同步的更新。由于CA模型的狀態(tài)變量是由空間的位置緊密相連,因此,能夠彌補(bǔ)Markov的預(yù)測(cè),不包括空間分布的缺點(diǎn),通過(guò)Idrisi軟件中的Markov模塊和CA-Markov模塊將兩者結(jié)合起來(lái)對(duì)土地覆蓋類(lèi)型的變化過(guò)程進(jìn)行預(yù)測(cè),具體操作步驟如下:

    1) 將研究區(qū)土地覆蓋類(lèi)型圖在ArcGis中轉(zhuǎn)換為ASCII柵格圖,然后在Idrisi中導(dǎo)入,并基于Markov模塊,獲取不同時(shí)間段內(nèi)的土地覆蓋類(lèi)型的轉(zhuǎn)移矩陣。2) 其次確定土地轉(zhuǎn)變適宜性圖集和CA濾波器。土地轉(zhuǎn)變適宜性圖集是約束土地覆蓋的影響因子,選擇地形坡度和道路作為影響耕地、林地和居民地的影響因子,在MCE中獲取土地轉(zhuǎn)變適宜性圖集。CA濾波器的作用是用于產(chǎn)生具有顯著空間特征的權(quán)重因子,本文采用5×5即150 m×150 m范圍內(nèi)的濾波器。3) 最后確定起始時(shí)刻和CA循環(huán)次數(shù)。研究起始時(shí)間是2010年古交市土地覆蓋類(lèi)型分布,設(shè)置CA循環(huán)次數(shù)為10,模擬2020年古交市土地覆蓋類(lèi)型分布狀況。

    4結(jié)果與分析

    4.1土地利用動(dòng)態(tài)變化特征

    根據(jù)1990年和2010年古交市的土地覆蓋評(píng)價(jià)結(jié)果,統(tǒng)計(jì)20多年來(lái)古交市不同土地覆蓋類(lèi)型的土地覆蓋面積,進(jìn)而反映區(qū)域土地覆蓋變化。面積變化首先反映在不同土地覆蓋類(lèi)型的總量變化上,通過(guò)研究分析土地覆蓋類(lèi)型的總量變化,可了解土地覆蓋變化的總體態(tài)勢(shì)和土地覆蓋類(lèi)型結(jié)構(gòu)的改變。根據(jù)1990年、2010年兩期土地覆蓋類(lèi)型圖,在ArcGIS10中運(yùn)用Spatial Analysis的功能,獲取古交市20多年來(lái)的土地覆蓋類(lèi)型面積轉(zhuǎn)移矩陣,見(jiàn)表1.

    為便于更好地研究土地覆蓋動(dòng)態(tài)變化過(guò)程的主導(dǎo)因素,本研究在土地覆蓋面積轉(zhuǎn)移矩陣的基礎(chǔ)上,建立土地覆蓋轉(zhuǎn)換模型,獲取土地覆蓋轉(zhuǎn)移的概率矩陣,不僅可以從定量的角度來(lái)描述古交市礦區(qū)土地覆蓋變化的速度,還能比較土地覆蓋變化的區(qū)域性差異,

    表11990—2010年土地覆蓋狀況轉(zhuǎn)移面積表km2

    并以此作為基礎(chǔ)來(lái)預(yù)測(cè)未來(lái)土地覆蓋時(shí)空動(dòng)態(tài)變化的趨勢(shì)。具體公式如下:

    (1)

    式中:

    S—土地覆蓋類(lèi)型;

    i—轉(zhuǎn)變?yōu)橥恋馗采w類(lèi)型j的轉(zhuǎn)換概率;

    Lij—土地覆蓋類(lèi)型i轉(zhuǎn)變?yōu)橥恋馗采w類(lèi)型j的面積。

    根據(jù)式(1)計(jì)算出土地覆蓋轉(zhuǎn)換概率矩陣,見(jiàn)表2.

    表21990—2010年土地覆蓋轉(zhuǎn)換概率矩陣表%

    為了更好地了解土地覆蓋的年變化規(guī)律,引入單一土地覆蓋類(lèi)型動(dòng)態(tài)度,其表達(dá)式為:

    (2)

    式中:

    M—單一土地覆蓋類(lèi)型動(dòng)態(tài)度;

    T—某一研究時(shí)段;

    Aa—研究期初的某一種土地覆蓋類(lèi)型面積;

    Ab—研究期末的某一種土地覆蓋類(lèi)型面積。

    根據(jù)公式(2)計(jì)算出古交市土地覆蓋類(lèi)型的年變化率表,見(jiàn)表3.

    表31990年、2010年土地覆蓋類(lèi)型年變化率表

    由表1,表2,表3可以看出,古交市近20多年來(lái)的土地利用變化過(guò)程,只有林地的面積大幅度增加,草地面積大幅度減少,此外耕地、城鎮(zhèn)及工礦用地和水域面積是逐漸減少的。其中草地的年均減少率為1.98%,面積由1990年的366.95 km2減少到206.82 km2,凈減少160.13 km2,減少幅度為43.64%.草地主要轉(zhuǎn)變?yōu)楦睾土值?,轉(zhuǎn)移概率為28.79%和30.71%,轉(zhuǎn)移面積105.66 km2和112.70 km2,主要原因是土地整理等政策所致。耕地年變化率0.91%,整體而言耕地的面積變化不大,凈減少了151.93 km2.從轉(zhuǎn)移矩陣類(lèi)型來(lái)看,耕地主要轉(zhuǎn)出類(lèi)型為林地,轉(zhuǎn)移概率為30.21%,主要原因是由于國(guó)家倡導(dǎo)的退耕還林和土地整治,礦區(qū)復(fù)墾等政策的影響。城鎮(zhèn)建設(shè)用地和水域變化不大。城鎮(zhèn)建設(shè)用地年變化率-0.30%,凈減少了5.54 km2,主要是受?chē)?guó)家政策土地整治和礦區(qū)恢復(fù)治理等政策的引導(dǎo),轉(zhuǎn)換為了耕地。水域減少了6.87 km2,年減少率為0.74%,主要轉(zhuǎn)換為了耕地。林地的面積受到國(guó)家退耕還林和草地開(kāi)荒,水土保持等政策的引導(dǎo),近20多年來(lái)增加了325.47 km2,年增長(zhǎng)率達(dá)到了5.49%,主要的來(lái)源是草地和耕地,轉(zhuǎn)換率達(dá)到了30.21%和30.71%,轉(zhuǎn)換了228.66 km2和112.70 km2.以上變化表明,古交市20多年來(lái)的土地利用變化主要是林地面積的增加及城鎮(zhèn)建設(shè)用地和水域輕微減少,此外是草地和耕地的大量減少。

    4.2基于CA-Markov土地覆蓋動(dòng)態(tài)預(yù)測(cè)

    4.2.1模擬精度檢驗(yàn)

    以古交市1990年、2000年兩年的土地覆蓋類(lèi)型數(shù)據(jù)為基礎(chǔ),運(yùn)用CA-Markov模型獲取2010年土地覆蓋類(lèi)型預(yù)測(cè)圖,與實(shí)際解譯的2010年土地覆蓋類(lèi)型分布圖進(jìn)行面積對(duì)比,獲取誤差對(duì)比表見(jiàn)表4.由表4可以看出,面積誤差最大為11.81%,最小為0.16%,精度接近90%.因此,可以用CA-Markov模型來(lái)預(yù)測(cè)古交市土地利用類(lèi)型變化。

    表42010年土地覆蓋類(lèi)型模擬精度評(píng)價(jià)表

    4.2.2研究區(qū)2020年土地覆蓋類(lèi)型模擬結(jié)果

    根據(jù)上述的研究,在IDRISI 軟件中的 CA-Markov模塊下,依次輸入2010年古交市土地覆蓋類(lèi)型分布圖、2000—2010年古交市土地覆蓋類(lèi)型轉(zhuǎn)移面積矩陣,以及2010年土地轉(zhuǎn)變適宜性圖集,對(duì)古交市2020年的土地覆蓋類(lèi)型狀況進(jìn)行模擬,結(jié)果見(jiàn)圖3和表5.

    圖3 古交市礦區(qū)2020年土地覆蓋類(lèi)型圖

    根據(jù)預(yù)測(cè)結(jié)果可以看出,到2020年受到國(guó)家政策的影響,礦區(qū)耕地和草地的面積進(jìn)一步大幅度減少,年減少率為1.74%和2.27%,凈減少面積115.74 km2和51.76 km2. 林地面積繼續(xù)增加,年增長(zhǎng)率達(dá)到了2.18%,凈增長(zhǎng)面積是142.45 km2,主要是由于退耕還林等政策的影響,轉(zhuǎn)入的類(lèi)型主要是耕地和草地。城鎮(zhèn)和工礦用地以及水域,從1990年到2020年的30多年內(nèi)呈現(xiàn)了一種先減少后增加的趨勢(shì),其中城鎮(zhèn)增長(zhǎng)幅度最大,達(dá)到了2.67%,凈增長(zhǎng)面積是23.37 km2,主要增長(zhǎng)區(qū)域位于礦區(qū)的西南部;水域變化不大,增長(zhǎng)率僅為0.38%.

    表52010年和2020年礦區(qū)土地覆蓋變化狀況表

    5結(jié)論

    1) 以古交市礦區(qū)為例,通過(guò)監(jiān)督分類(lèi)獲取1990年、2000年、2010年3期土地覆蓋類(lèi)型圖,通過(guò)轉(zhuǎn)移矩陣和CA-Markov模型分別從研究區(qū)的面積、速度和空間分布3方面對(duì)土地覆蓋類(lèi)型的時(shí)空變化特征進(jìn)行分析;同時(shí)通過(guò)預(yù)測(cè)精度驗(yàn)證獲取CA-Markov模型面積,精度接近90%,表明CA-Markov模型可以反應(yīng)研究區(qū)土地覆蓋類(lèi)型的動(dòng)態(tài)變化狀況,預(yù)測(cè)結(jié)果可為礦區(qū)土地覆蓋類(lèi)型整治及相關(guān)規(guī)劃的編制提供依據(jù)。

    2) 在上述基礎(chǔ)上,利用動(dòng)態(tài)模型對(duì)研究區(qū)1990—2010年的土地覆蓋類(lèi)型進(jìn)行分析,結(jié)果表明:受到國(guó)家退耕還林和土地整理等政策的影響,林地的面積逐年增加,年增長(zhǎng)率達(dá)到了5.49%,主要的轉(zhuǎn)入類(lèi)型為耕地和草地,受此影響耕地和草地的面積大幅度減少,尤其是草地的面積,此外城鎮(zhèn)建設(shè)用地和水域輕微減少。

    3) 利用CA-Markov模型對(duì)研究區(qū)2020年的土地覆蓋類(lèi)型進(jìn)行預(yù)測(cè),結(jié)果表明耕地和草地的面積進(jìn)一步減少,年減少率為1.74%和2.27%.城鎮(zhèn)和工礦用地以及水域的面積小幅度增加,林地的面積繼續(xù)大幅度增長(zhǎng),年增長(zhǎng)率達(dá)到了2.18%,凈增長(zhǎng)面積是142.45 km2.

    參考文獻(xiàn)

    [1]劉紀(jì)遠(yuǎn),張?jiān)鱿?,徐新良,?21世紀(jì)初中國(guó)土地利用變化的空間格局與驅(qū)動(dòng)力分析[J].地理學(xué)報(bào),2009,64(12):1416-1420.

    [2]何英彬,姚艷敏,唐華俊,等.土地利用/覆蓋變化驅(qū)動(dòng)力機(jī)制研究新進(jìn)展[J].中國(guó)農(nóng)學(xué)通報(bào),2013,29(2):190-195.

    [3]毛慧琴,延曉冬,熊喆.土地利用變化對(duì)氣候影響的研究進(jìn)展[J].氣候與環(huán)境研究,2011,16(4):517-524.

    [4]黎治華,高志強(qiáng),高煒,等.中國(guó)1999—2009年土地覆蓋動(dòng)態(tài)變化的時(shí)空特點(diǎn)[J].農(nóng)業(yè)工程學(xué)報(bào),2011,27(2):318-323.

    [5]馮異星,羅格平,周德成,等.近50a土地利用變化對(duì)干旱區(qū)典型流域景觀格局的影響—以新疆瑪納斯河流域?yàn)槔齕J].生態(tài)學(xué)報(bào),2010,30(16):4298-4305.

    [6]楊志平.徐州市土地利用變化及其驅(qū)動(dòng)機(jī)制分析[J].國(guó)土與自然資源研究,2012(5):1-4.

    [7]吳莉,侯西勇,徐新良,等.山東沿海地區(qū)土地利用和景觀格局變化[J].農(nóng)業(yè)工程學(xué)報(bào),2013,29(5):207-212.

    [8]胡召玲,杜培軍,趙昕.徐州煤礦區(qū)土地利用變化分析[J].地理學(xué)報(bào),2007(11):1209-1214.

    [9]李琳,袁春,周偉,等.平朔露天礦區(qū)土地利用/覆蓋變化分析[J].資源與產(chǎn)業(yè),2007(3):5-9.

    [10]楊娟,王昌全,夏建國(guó),等.基于元胞自動(dòng)機(jī)的土地利用空間規(guī)劃輔助研究—以眉山市東坡區(qū)為例[J].土壤學(xué)報(bào),2010(5):847-852.

    Research on Space-time Dynamic Change of

    Land Coverage in Gujiao Mining Area

    LI Honghong

    Key words CA-MARKOV; Temporal and Spatial variation; Land cover types; Gujiao mining area

    猜你喜歡
    時(shí)空變化
    遼西北地區(qū)干旱時(shí)空變化特征分析
    中國(guó)畜禽溫室氣體排放時(shí)空變化及影響因素研究
    基于主成分分析法及熵值法結(jié)合的耕地集約利用時(shí)空變化
    商(2016年29期)2016-10-29 11:10:35
    國(guó)家級(jí)流通節(jié)點(diǎn)城市流通業(yè)效率的時(shí)空變化研究
    江蘇經(jīng)濟(jì)發(fā)展區(qū)域差異時(shí)空變化研究
    南海聲速主躍層的時(shí)空分布
    基于RS與GIS的黃岡市土地利用時(shí)空變化分析
    基于DSSAT和GIS的中國(guó)水稻產(chǎn)量時(shí)空變化模擬研究
    鳳城市耕地土壤pH值的時(shí)空變化
    鳳城市耕地土壤有機(jī)質(zhì)時(shí)空變化研究
    美女cb高潮喷水在线观看| 18禁黄网站禁片午夜丰满| 成人性生交大片免费视频hd| 亚洲人成电影免费在线| 97超视频在线观看视频| 听说在线观看完整版免费高清| 不卡一级毛片| 黄色日韩在线| 在线a可以看的网站| 永久网站在线| 人人妻,人人澡人人爽秒播| 最好的美女福利视频网| 美女免费视频网站| 久久精品夜夜夜夜夜久久蜜豆| 国产精品一及| 欧美午夜高清在线| 男女床上黄色一级片免费看| 欧美日韩中文字幕国产精品一区二区三区| 国产黄色小视频在线观看| 一个人观看的视频www高清免费观看| a在线观看视频网站| av在线天堂中文字幕| 亚洲国产色片| 国产欧美日韩一区二区三| 色噜噜av男人的天堂激情| 欧美午夜高清在线| 久久久久久大精品| 亚洲aⅴ乱码一区二区在线播放| 亚洲欧美日韩高清专用| 欧美中文日本在线观看视频| 欧洲精品卡2卡3卡4卡5卡区| 在线播放国产精品三级| 成人精品一区二区免费| 一个人看的www免费观看视频| av在线蜜桃| 午夜福利18| 欧美日本亚洲视频在线播放| 黄色女人牲交| 深夜精品福利| 午夜亚洲福利在线播放| 青草久久国产| 国产69精品久久久久777片| 国产在视频线在精品| 可以在线观看毛片的网站| 九色国产91popny在线| 老司机深夜福利视频在线观看| 在线a可以看的网站| 日本在线视频免费播放| 性欧美人与动物交配| 免费看日本二区| 国产精品电影一区二区三区| 搞女人的毛片| 一个人看的www免费观看视频| 色哟哟哟哟哟哟| 999久久久精品免费观看国产| 国产男靠女视频免费网站| 久久精品国产自在天天线| 人人妻人人澡欧美一区二区| 久久久久久大精品| 午夜福利视频1000在线观看| 午夜精品在线福利| 国产精品亚洲av一区麻豆| 国内久久婷婷六月综合欲色啪| 美女被艹到高潮喷水动态| 国产乱人伦免费视频| 永久网站在线| 色在线成人网| 少妇人妻一区二区三区视频| 国产私拍福利视频在线观看| av在线天堂中文字幕| 国产欧美日韩一区二区精品| bbb黄色大片| 精品久久国产蜜桃| 观看美女的网站| 日本精品一区二区三区蜜桃| 露出奶头的视频| 日本 欧美在线| 亚洲国产精品sss在线观看| 国产欧美日韩精品亚洲av| 国产成年人精品一区二区| a级毛片免费高清观看在线播放| 91久久精品国产一区二区成人| 真人一进一出gif抽搐免费| 美女高潮的动态| 午夜久久久久精精品| netflix在线观看网站| 麻豆成人av在线观看| 欧美+亚洲+日韩+国产| 婷婷六月久久综合丁香| 舔av片在线| 毛片女人毛片| 国产精品亚洲av一区麻豆| 高潮久久久久久久久久久不卡| 老司机深夜福利视频在线观看| 国产视频一区二区在线看| 国产精品久久久久久人妻精品电影| 日韩国内少妇激情av| 伦理电影大哥的女人| 九色国产91popny在线| 欧美高清性xxxxhd video| 久久久国产成人免费| bbb黄色大片| 国产成人啪精品午夜网站| 老司机福利观看| 国产单亲对白刺激| 欧美一级a爱片免费观看看| 亚洲av成人精品一区久久| 亚洲av五月六月丁香网| 人妻丰满熟妇av一区二区三区| 18+在线观看网站| 亚洲中文字幕日韩| 日韩成人在线观看一区二区三区| 亚洲18禁久久av| 午夜a级毛片| 免费在线观看影片大全网站| 亚洲一区二区三区色噜噜| .国产精品久久| 久久久久久九九精品二区国产| 国产主播在线观看一区二区| 亚洲狠狠婷婷综合久久图片| 午夜老司机福利剧场| 最近中文字幕高清免费大全6 | 亚洲精品日韩av片在线观看| 国产精品久久久久久人妻精品电影| 在线播放国产精品三级| 嫩草影视91久久| 我要看日韩黄色一级片| 亚洲不卡免费看| 亚洲国产高清在线一区二区三| 午夜a级毛片| 性色avwww在线观看| 欧美一级a爱片免费观看看| 国产高潮美女av| 亚洲av日韩精品久久久久久密| 女人被狂操c到高潮| 成人国产综合亚洲| 亚洲国产欧洲综合997久久,| 亚洲欧美激情综合另类| 黄色日韩在线| 综合色av麻豆| 亚洲美女搞黄在线观看 | 中文资源天堂在线| 99久久精品热视频| 在现免费观看毛片| 午夜免费男女啪啪视频观看 | 国产精品综合久久久久久久免费| 免费大片18禁| 国产色爽女视频免费观看| 每晚都被弄得嗷嗷叫到高潮| 观看美女的网站| 在线播放国产精品三级| 婷婷精品国产亚洲av在线| 亚洲七黄色美女视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲中文字幕一区二区三区有码在线看| 久久久久久久午夜电影| 熟女人妻精品中文字幕| 亚洲欧美日韩高清在线视频| 欧美国产日韩亚洲一区| 淫秽高清视频在线观看| 91久久精品国产一区二区成人| 日韩有码中文字幕| 精品日产1卡2卡| 久久久久精品国产欧美久久久| 亚洲人成网站在线播| 亚洲在线观看片| 十八禁国产超污无遮挡网站| 欧美黄色淫秽网站| 九色国产91popny在线| 成人特级黄色片久久久久久久| 日本与韩国留学比较| 国产高清视频在线观看网站| 日韩av在线大香蕉| 国内精品久久久久久久电影| 99精品在免费线老司机午夜| 三级男女做爰猛烈吃奶摸视频| 久久久久性生活片| 97超级碰碰碰精品色视频在线观看| 国产又黄又爽又无遮挡在线| 久久久久久久久中文| 日韩欧美在线乱码| 国产精品乱码一区二三区的特点| 精品久久久久久久久久免费视频| 国产熟女xx| 在线观看av片永久免费下载| 日本五十路高清| 国产av不卡久久| 成年人黄色毛片网站| 老鸭窝网址在线观看| 一区二区三区高清视频在线| www.熟女人妻精品国产| 精品乱码久久久久久99久播| 在线观看午夜福利视频| 别揉我奶头 嗯啊视频| 直男gayav资源| 别揉我奶头 嗯啊视频| 欧美日韩乱码在线| 国产高清视频在线观看网站| 亚洲成人久久爱视频| 国产探花极品一区二区| 欧美+亚洲+日韩+国产| 国产av在哪里看| 国产精品女同一区二区软件 | 日本 av在线| av福利片在线观看| 亚洲av成人不卡在线观看播放网| 亚洲性夜色夜夜综合| 久久久久性生活片| h日本视频在线播放| 宅男免费午夜| 日韩欧美精品免费久久 | 亚洲国产高清在线一区二区三| 亚洲成a人片在线一区二区| 真人一进一出gif抽搐免费| 男女之事视频高清在线观看| www.色视频.com| 毛片女人毛片| 91麻豆av在线| 午夜两性在线视频| 他把我摸到了高潮在线观看| 国产人妻一区二区三区在| 日日摸夜夜添夜夜添小说| 国产高潮美女av| 淫秽高清视频在线观看| aaaaa片日本免费| netflix在线观看网站| 久久国产乱子免费精品| 69人妻影院| av黄色大香蕉| 久久香蕉精品热| 一a级毛片在线观看| 夜夜爽天天搞| 日本一本二区三区精品| 日本免费一区二区三区高清不卡| 欧美性猛交╳xxx乱大交人| 成年女人毛片免费观看观看9| 激情在线观看视频在线高清| 久久精品国产清高在天天线| 最近最新免费中文字幕在线| 国产在线精品亚洲第一网站| 老熟妇仑乱视频hdxx| 最近视频中文字幕2019在线8| 一夜夜www| 又爽又黄无遮挡网站| 久久久久免费精品人妻一区二区| 3wmmmm亚洲av在线观看| 69av精品久久久久久| 国产成人a区在线观看| 日韩精品青青久久久久久| 欧美一区二区亚洲| 色av中文字幕| 国产精品不卡视频一区二区 | 久久午夜福利片| 国产成人啪精品午夜网站| ponron亚洲| 久久久成人免费电影| 成人三级黄色视频| 最后的刺客免费高清国语| 午夜精品一区二区三区免费看| 欧美日韩瑟瑟在线播放| 国产乱人伦免费视频| 精品一区二区免费观看| 精品免费久久久久久久清纯| 亚洲avbb在线观看| 99视频精品全部免费 在线| 99久久精品热视频| 性色avwww在线观看| 免费人成在线观看视频色| 国产精品影院久久| 国产国拍精品亚洲av在线观看| 高清毛片免费观看视频网站| 99久国产av精品| 免费在线观看日本一区| 国产成年人精品一区二区| 亚洲五月婷婷丁香| 真实男女啪啪啪动态图| 欧美不卡视频在线免费观看| 亚洲片人在线观看| 亚洲欧美清纯卡通| 美女被艹到高潮喷水动态| 3wmmmm亚洲av在线观看| 观看美女的网站| 国产 一区 欧美 日韩| 波多野结衣高清无吗| 天美传媒精品一区二区| 国内久久婷婷六月综合欲色啪| 亚洲精品在线美女| 国产三级黄色录像| xxxwww97欧美| 亚洲人成电影免费在线| 国产色婷婷99| 国产真实伦视频高清在线观看 | 精品乱码久久久久久99久播| 亚洲人与动物交配视频| 黄色配什么色好看| 黄片小视频在线播放| av女优亚洲男人天堂| 日日夜夜操网爽| 欧美+亚洲+日韩+国产| 一级毛片久久久久久久久女| 又爽又黄无遮挡网站| 韩国av一区二区三区四区| 久久国产乱子免费精品| 久久6这里有精品| 最新在线观看一区二区三区| 97超级碰碰碰精品色视频在线观看| 1024手机看黄色片| 波多野结衣巨乳人妻| 日韩欧美 国产精品| 国产精品一及| 淫秽高清视频在线观看| 国产亚洲欧美98| 国产精品一区二区三区四区免费观看 | 无人区码免费观看不卡| 免费高清视频大片| 亚洲成av人片免费观看| 精品国内亚洲2022精品成人| 亚洲中文字幕一区二区三区有码在线看| 丰满乱子伦码专区| 青草久久国产| 欧美成狂野欧美在线观看| 蜜桃久久精品国产亚洲av| 99热精品在线国产| 人妻夜夜爽99麻豆av| 欧美日韩乱码在线| 性欧美人与动物交配| 国产探花在线观看一区二区| 国产男靠女视频免费网站| 欧美一区二区国产精品久久精品| 制服丝袜大香蕉在线| 成人av在线播放网站| 欧洲精品卡2卡3卡4卡5卡区| 两性午夜刺激爽爽歪歪视频在线观看| 一a级毛片在线观看| a在线观看视频网站| 亚洲不卡免费看| а√天堂www在线а√下载| 99久久无色码亚洲精品果冻| 在线观看午夜福利视频| 看黄色毛片网站| 国产成人av教育| 特大巨黑吊av在线直播| 国产精品乱码一区二三区的特点| 亚洲激情在线av| 精品欧美国产一区二区三| 国产亚洲精品综合一区在线观看| 99久久精品热视频| 久久久国产成人免费| 啦啦啦韩国在线观看视频| 99在线视频只有这里精品首页| 国产精品电影一区二区三区| 亚洲片人在线观看| 亚洲成a人片在线一区二区| 在线天堂最新版资源| 九九热线精品视视频播放| 国产成人啪精品午夜网站| 日韩人妻高清精品专区| 成人无遮挡网站| 国产精品不卡视频一区二区 | 高清毛片免费观看视频网站| 国产午夜精品论理片| 国产精品av视频在线免费观看| 午夜免费激情av| 啦啦啦观看免费观看视频高清| 精品国产亚洲在线| 久久久久久国产a免费观看| 欧洲精品卡2卡3卡4卡5卡区| 99久久久亚洲精品蜜臀av| 亚洲av电影在线进入| 中文字幕免费在线视频6| 国产极品精品免费视频能看的| 国产三级黄色录像| 成熟少妇高潮喷水视频| 免费在线观看影片大全网站| 国产精品自产拍在线观看55亚洲| 国产单亲对白刺激| 精品国产亚洲在线| 欧美一区二区亚洲| 又爽又黄无遮挡网站| 国产精品一区二区三区四区久久| 岛国在线免费视频观看| 非洲黑人性xxxx精品又粗又长| 国产精品乱码一区二三区的特点| 91麻豆精品激情在线观看国产| 日本一二三区视频观看| 国产探花在线观看一区二区| 日韩欧美精品v在线| 俺也久久电影网| 欧美日韩瑟瑟在线播放| 麻豆国产97在线/欧美| 在线观看午夜福利视频| 三级男女做爰猛烈吃奶摸视频| 欧美激情在线99| 国产精品一区二区三区四区免费观看 | 亚洲不卡免费看| 午夜福利欧美成人| 99热这里只有精品一区| 99久久精品热视频| 久久国产精品人妻蜜桃| 97超视频在线观看视频| 亚洲av成人精品一区久久| 免费电影在线观看免费观看| 高潮久久久久久久久久久不卡| 十八禁网站免费在线| 国产大屁股一区二区在线视频| 久久精品久久久久久噜噜老黄 | 首页视频小说图片口味搜索| 一级a爱片免费观看的视频| 久久久久性生活片| 久久国产乱子免费精品| 老司机福利观看| 亚洲在线自拍视频| 国产三级黄色录像| 国产成+人综合+亚洲专区| 51午夜福利影视在线观看| 久久久久免费精品人妻一区二区| 听说在线观看完整版免费高清| 久久久久免费精品人妻一区二区| 超碰av人人做人人爽久久| 国产精品不卡视频一区二区 | 欧美不卡视频在线免费观看| 男女视频在线观看网站免费| .国产精品久久| www.www免费av| 午夜两性在线视频| 一本一本综合久久| ponron亚洲| 精品人妻熟女av久视频| 日韩欧美国产一区二区入口| 亚洲成a人片在线一区二区| 久久国产乱子伦精品免费另类| 首页视频小说图片口味搜索| 无遮挡黄片免费观看| 亚洲人成伊人成综合网2020| 亚洲美女黄片视频| 麻豆av噜噜一区二区三区| 精品久久久久久,| 欧美3d第一页| 欧美bdsm另类| 高清在线国产一区| 88av欧美| 精品一区二区三区av网在线观看| 丰满的人妻完整版| 亚洲在线自拍视频| 国产一区二区亚洲精品在线观看| 色5月婷婷丁香| 欧美一区二区亚洲| 欧美最黄视频在线播放免费| 国产乱人视频| 波多野结衣巨乳人妻| 日韩中字成人| 亚洲欧美激情综合另类| 国产极品精品免费视频能看的| 女人十人毛片免费观看3o分钟| 亚洲国产日韩欧美精品在线观看| 亚州av有码| 国产成人啪精品午夜网站| 黄片小视频在线播放| 丁香六月欧美| 亚洲七黄色美女视频| 小说图片视频综合网站| 国产精品美女特级片免费视频播放器| 国产午夜精品久久久久久一区二区三区 | 国内精品美女久久久久久| 亚洲人与动物交配视频| 五月伊人婷婷丁香| 男女那种视频在线观看| 长腿黑丝高跟| 色综合婷婷激情| 亚洲七黄色美女视频| 国产精品亚洲一级av第二区| 三级毛片av免费| 亚洲精品成人久久久久久| 中文字幕人妻熟人妻熟丝袜美| 亚洲内射少妇av| 久久久久久久久久成人| 99riav亚洲国产免费| 一区二区三区高清视频在线| 美女高潮的动态| 国产精品av视频在线免费观看| 淫妇啪啪啪对白视频| 舔av片在线| 老熟妇乱子伦视频在线观看| 亚洲av免费高清在线观看| 国产av在哪里看| 国产成人啪精品午夜网站| 免费高清视频大片| 欧美绝顶高潮抽搐喷水| 老鸭窝网址在线观看| 悠悠久久av| 国产单亲对白刺激| 精品乱码久久久久久99久播| 亚洲成av人片在线播放无| 色播亚洲综合网| 亚洲 欧美 日韩 在线 免费| 国产日本99.免费观看| 91午夜精品亚洲一区二区三区 | 脱女人内裤的视频| 男女之事视频高清在线观看| 动漫黄色视频在线观看| 日本黄色片子视频| 国产精品一及| 亚洲国产日韩欧美精品在线观看| 一区二区三区免费毛片| 成人国产一区最新在线观看| 欧洲精品卡2卡3卡4卡5卡区| 国产三级黄色录像| 亚洲色图av天堂| 国产精品嫩草影院av在线观看 | 宅男免费午夜| 韩国av一区二区三区四区| 久久天躁狠狠躁夜夜2o2o| 不卡一级毛片| 中出人妻视频一区二区| 亚洲18禁久久av| 久久久久久久久久成人| 青草久久国产| 日韩成人在线观看一区二区三区| 国产色爽女视频免费观看| 国产精品自产拍在线观看55亚洲| 在线a可以看的网站| 三级毛片av免费| 国产高清激情床上av| 99热精品在线国产| 亚洲欧美日韩卡通动漫| 我的老师免费观看完整版| 性欧美人与动物交配| 在线免费观看的www视频| 搡老熟女国产l中国老女人| 日韩高清综合在线| 人人妻人人澡欧美一区二区| 国产 一区 欧美 日韩| 听说在线观看完整版免费高清| 麻豆成人av在线观看| 很黄的视频免费| 国产免费av片在线观看野外av| 美女被艹到高潮喷水动态| 国产三级黄色录像| 99久久99久久久精品蜜桃| 真实男女啪啪啪动态图| 一本一本综合久久| 国产老妇女一区| 国产亚洲av嫩草精品影院| 一个人看的www免费观看视频| 97热精品久久久久久| 亚洲国产高清在线一区二区三| 欧美成狂野欧美在线观看| 国内少妇人妻偷人精品xxx网站| 青草久久国产| 免费在线观看影片大全网站| 国内精品久久久久久久电影| 97人妻精品一区二区三区麻豆| xxxwww97欧美| 亚洲成人精品中文字幕电影| 欧美极品一区二区三区四区| 91av网一区二区| 日本在线视频免费播放| 深夜精品福利| 国产激情偷乱视频一区二区| 露出奶头的视频| 亚洲,欧美,日韩| 99国产精品一区二区蜜桃av| 免费av不卡在线播放| 免费观看的影片在线观看| 1000部很黄的大片| 热99re8久久精品国产| 午夜影院日韩av| 精品99又大又爽又粗少妇毛片 | 中文字幕人成人乱码亚洲影| 高清在线国产一区| 男女床上黄色一级片免费看| 真人做人爱边吃奶动态| 综合色av麻豆| 天天一区二区日本电影三级| 18美女黄网站色大片免费观看| 国产一区二区三区视频了| 51国产日韩欧美| 亚洲精品日韩av片在线观看| 欧美国产日韩亚洲一区| 美女高潮喷水抽搐中文字幕| 99热精品在线国产| www.999成人在线观看| 亚洲美女搞黄在线观看 | av天堂在线播放| 国产在线精品亚洲第一网站| 亚洲欧美日韩卡通动漫| 欧美+亚洲+日韩+国产| 国产综合懂色| 观看美女的网站| 欧美一区二区精品小视频在线| 直男gayav资源| 我要搜黄色片| 国产高潮美女av| 99久久成人亚洲精品观看| 国产精品女同一区二区软件 | 国产黄a三级三级三级人| 国产黄片美女视频| 夜夜夜夜夜久久久久| 99热只有精品国产| 国产精品国产高清国产av| 99热6这里只有精品| 一进一出好大好爽视频| 一个人看视频在线观看www免费| 久99久视频精品免费| 丰满的人妻完整版| 久久精品国产亚洲av涩爱 | 一a级毛片在线观看| 99久久99久久久精品蜜桃| 欧美成狂野欧美在线观看| 一区福利在线观看| 久久久久免费精品人妻一区二区| 久9热在线精品视频| www.色视频.com| 日本熟妇午夜| 国产中年淑女户外野战色|