• <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í)空變化研究
    久久99精品国语久久久| 亚洲久久久国产精品| 免费黄色在线免费观看| 99热全是精品| 久久久久国产精品人妻一区二区| 欧美xxxx性猛交bbbb| 国产精品一二三区在线看| 欧美成人午夜精品| 晚上一个人看的免费电影| 另类亚洲欧美激情| 美女中出高潮动态图| 国产一区有黄有色的免费视频| 久热这里只有精品99| 亚洲中文av在线| 成年美女黄网站色视频大全免费| 国产极品天堂在线| 亚洲av欧美aⅴ国产| 伦理电影免费视频| 日本av免费视频播放| av免费观看日本| 免费女性裸体啪啪无遮挡网站| 伦精品一区二区三区| 久久狼人影院| 亚洲国产看品久久| 两个人免费观看高清视频| 在线观看一区二区三区激情| 免费人成在线观看视频色| 国产精品蜜桃在线观看| 成年人午夜在线观看视频| 韩国高清视频一区二区三区| 国产乱来视频区| 乱人伦中国视频| 成人免费观看视频高清| 日日撸夜夜添| 国产精品久久久久成人av| 丝瓜视频免费看黄片| 岛国毛片在线播放| 美女国产视频在线观看| 9191精品国产免费久久| 亚洲欧美色中文字幕在线| 一级毛片黄色毛片免费观看视频| 又粗又硬又长又爽又黄的视频| 精品卡一卡二卡四卡免费| 亚洲精品国产av蜜桃| 国产欧美亚洲国产| 亚洲伊人久久精品综合| 热99国产精品久久久久久7| 亚洲国产毛片av蜜桃av| 国产永久视频网站| 亚洲av电影在线观看一区二区三区| 九色亚洲精品在线播放| 久久精品国产亚洲av天美| 99国产精品免费福利视频| 欧美亚洲 丝袜 人妻 在线| 国产亚洲最大av| 97精品久久久久久久久久精品| 精品一品国产午夜福利视频| 国产精品国产三级专区第一集| 欧美另类一区| 妹子高潮喷水视频| 黄色一级大片看看| 成人手机av| 精品第一国产精品| 亚洲精品成人av观看孕妇| 高清在线视频一区二区三区| 美女视频免费永久观看网站| 亚洲国产成人一精品久久久| 侵犯人妻中文字幕一二三四区| 久久午夜福利片| freevideosex欧美| 蜜臀久久99精品久久宅男| www.色视频.com| 国产成人av激情在线播放| 久久99热这里只频精品6学生| 亚洲国产色片| 国产乱人偷精品视频| 赤兔流量卡办理| 免费观看av网站的网址| 日本av免费视频播放| 日本av免费视频播放| 一级片'在线观看视频| 好男人视频免费观看在线| 日韩 亚洲 欧美在线| 久热久热在线精品观看| 亚洲少妇的诱惑av| 午夜视频国产福利| 国产高清不卡午夜福利| 一级毛片黄色毛片免费观看视频| 如日韩欧美国产精品一区二区三区| 国产精品秋霞免费鲁丝片| 爱豆传媒免费全集在线观看| 男女边吃奶边做爰视频| 国产片内射在线| av国产精品久久久久影院| 亚洲,一卡二卡三卡| 亚洲高清免费不卡视频| 精品久久久精品久久久| 日本爱情动作片www.在线观看| 国产精品女同一区二区软件| 久久精品久久久久久噜噜老黄| 伊人亚洲综合成人网| 18在线观看网站| 国产熟女午夜一区二区三区| 日韩三级伦理在线观看| 午夜久久久在线观看| 一区二区三区四区激情视频| 有码 亚洲区| 国产成人a∨麻豆精品| 亚洲欧美日韩另类电影网站| 欧美亚洲日本最大视频资源| 亚洲国产毛片av蜜桃av| 久久毛片免费看一区二区三区| 大陆偷拍与自拍| 超色免费av| 99久久中文字幕三级久久日本| 久久久久久久精品精品| 国产黄色视频一区二区在线观看| 久久精品国产a三级三级三级| 99国产精品免费福利视频| 蜜臀久久99精品久久宅男| 观看av在线不卡| 自拍欧美九色日韩亚洲蝌蚪91| 搡老乐熟女国产| 国产高清国产精品国产三级| 亚洲国产精品成人久久小说| 黑人巨大精品欧美一区二区蜜桃 | 少妇人妻 视频| 国产视频首页在线观看| 精品卡一卡二卡四卡免费| 亚洲精品中文字幕在线视频| 一区二区日韩欧美中文字幕 | 国产永久视频网站| 日本色播在线视频| 在线观看三级黄色| 侵犯人妻中文字幕一二三四区| 亚洲成av片中文字幕在线观看 | 国产麻豆69| 日韩欧美精品免费久久| 久久久欧美国产精品| 亚洲天堂av无毛| 国产免费一区二区三区四区乱码| 欧美变态另类bdsm刘玥| 亚洲美女视频黄频| 色婷婷久久久亚洲欧美| 成年人免费黄色播放视频| 日本wwww免费看| 大香蕉97超碰在线| 成人午夜精彩视频在线观看| 国产国语露脸激情在线看| 日韩免费高清中文字幕av| 精品视频人人做人人爽| 亚洲激情五月婷婷啪啪| 久久 成人 亚洲| 国产精品麻豆人妻色哟哟久久| 1024视频免费在线观看| 国产1区2区3区精品| 涩涩av久久男人的天堂| av播播在线观看一区| 免费观看性生交大片5| 亚洲美女搞黄在线观看| 日韩人妻精品一区2区三区| 亚洲av电影在线观看一区二区三区| 天天躁夜夜躁狠狠躁躁| 国产成人精品久久久久久| 一级毛片电影观看| 自线自在国产av| 免费观看在线日韩| 校园人妻丝袜中文字幕| 国产精品久久久av美女十八| 国产精品国产三级国产专区5o| 热99久久久久精品小说推荐| 大话2 男鬼变身卡| 草草在线视频免费看| 美女xxoo啪啪120秒动态图| 美女内射精品一级片tv| 麻豆精品久久久久久蜜桃| 一级,二级,三级黄色视频| 日韩制服骚丝袜av| 狂野欧美激情性xxxx在线观看| 国语对白做爰xxxⅹ性视频网站| 欧美日本中文国产一区发布| 亚洲av中文av极速乱| 日本色播在线视频| 一本色道久久久久久精品综合| 两个人看的免费小视频| 最近2019中文字幕mv第一页| 亚洲国产毛片av蜜桃av| 亚洲经典国产精华液单| 国产欧美亚洲国产| 精品一区二区三卡| 久久久久精品人妻al黑| 亚洲欧美精品自产自拍| 亚洲国产av新网站| 亚洲国产精品一区三区| 人妻少妇偷人精品九色| tube8黄色片| 老司机影院毛片| 国产精品国产三级专区第一集| 免费不卡的大黄色大毛片视频在线观看| 国产成人av激情在线播放| 在线 av 中文字幕| 香蕉精品网在线| 日韩欧美一区视频在线观看| 捣出白浆h1v1| 欧美日韩视频精品一区| 成人漫画全彩无遮挡| 中文字幕精品免费在线观看视频 | 狠狠婷婷综合久久久久久88av| 亚洲精品久久久久久婷婷小说| 免费在线观看黄色视频的| 777米奇影视久久| 亚洲 欧美一区二区三区| 亚洲欧美清纯卡通| 韩国av在线不卡| av线在线观看网站| 日韩电影二区| 黄色配什么色好看| 日本免费在线观看一区| 制服诱惑二区| 精品久久蜜臀av无| 91在线精品国自产拍蜜月| 久久 成人 亚洲| 侵犯人妻中文字幕一二三四区| 飞空精品影院首页| 精品人妻一区二区三区麻豆| 国产精品麻豆人妻色哟哟久久| 久久精品国产亚洲av涩爱| 最近中文字幕2019免费版| 久久久久久久久久久久大奶| 久久久久精品久久久久真实原创| 十分钟在线观看高清视频www| 国产激情久久老熟女| 免费观看无遮挡的男女| 国产 精品1| 男女午夜视频在线观看 | 亚洲色图 男人天堂 中文字幕 | 久久人人爽av亚洲精品天堂| 久久精品国产综合久久久 | 国产av码专区亚洲av| 欧美精品人与动牲交sv欧美| 蜜臀久久99精品久久宅男| 视频区图区小说| 国产一区亚洲一区在线观看| 亚洲国产精品一区三区| 免费观看无遮挡的男女| 咕卡用的链子| a 毛片基地| 欧美 日韩 精品 国产| 日本91视频免费播放| 人妻 亚洲 视频| 日韩熟女老妇一区二区性免费视频| 久久久久久久久久久久大奶| 久久精品久久久久久噜噜老黄| 人人妻人人澡人人看| 色吧在线观看| 成人18禁高潮啪啪吃奶动态图| av国产久精品久网站免费入址| 最近2019中文字幕mv第一页| 亚洲精品国产av蜜桃| 久久鲁丝午夜福利片| 国产精品秋霞免费鲁丝片| 两个人看的免费小视频| 视频中文字幕在线观看| 中文字幕av电影在线播放| 美女国产高潮福利片在线看| 激情视频va一区二区三区| 亚洲精品一区蜜桃| 一区二区日韩欧美中文字幕 | 熟女电影av网| 国产成人精品无人区| 久久久亚洲精品成人影院| 制服诱惑二区| 有码 亚洲区| 亚洲精品av麻豆狂野| 高清欧美精品videossex| 一边摸一边做爽爽视频免费| 欧美3d第一页| av在线老鸭窝| 国产精品麻豆人妻色哟哟久久| 国产不卡av网站在线观看| 精品人妻熟女毛片av久久网站| 成人免费观看视频高清| 大话2 男鬼变身卡| 日产精品乱码卡一卡2卡三| 精品国产乱码久久久久久小说| 99久久中文字幕三级久久日本| 少妇的丰满在线观看| 黄色怎么调成土黄色| 18禁在线无遮挡免费观看视频| 欧美日韩av久久| 性色av一级| 日韩制服丝袜自拍偷拍| 精品一品国产午夜福利视频| 国产毛片在线视频| 国产永久视频网站| 亚洲美女视频黄频| 久久久久久久精品精品| 热re99久久精品国产66热6| 免费大片黄手机在线观看| 日日撸夜夜添| 99久久综合免费| 亚洲欧美一区二区三区国产| 自线自在国产av| 色5月婷婷丁香| 日本免费在线观看一区| 91aial.com中文字幕在线观看| 视频区图区小说| 午夜激情久久久久久久| 高清不卡的av网站| 色婷婷久久久亚洲欧美| 亚洲五月色婷婷综合| av黄色大香蕉| 欧美老熟妇乱子伦牲交| 永久免费av网站大全| 婷婷色综合大香蕉| 亚洲欧洲精品一区二区精品久久久 | 久久鲁丝午夜福利片| 国产黄色免费在线视频| 女人精品久久久久毛片| 黄片播放在线免费| 日韩欧美精品免费久久| 久久久久久久国产电影| 亚洲精品久久成人aⅴ小说| 精品人妻熟女毛片av久久网站| 亚洲精品自拍成人| 久久人人爽人人爽人人片va| 性色avwww在线观看| 免费黄网站久久成人精品| 国产探花极品一区二区| 国产精品国产三级专区第一集| √禁漫天堂资源中文www| 一级毛片我不卡| 久久99热6这里只有精品| 国产日韩欧美在线精品| 在线天堂中文资源库| av播播在线观看一区| 亚洲成国产人片在线观看| av国产精品久久久久影院| 国产 精品1| 亚洲三级黄色毛片| 最近中文字幕2019免费版| 日本av免费视频播放| 国产片内射在线| 男人爽女人下面视频在线观看| 久久精品国产鲁丝片午夜精品| 免费大片黄手机在线观看| 日产精品乱码卡一卡2卡三| 狂野欧美激情性xxxx在线观看| 国产精品熟女久久久久浪| 涩涩av久久男人的天堂| 99视频精品全部免费 在线| 成人黄色视频免费在线看| 搡老乐熟女国产| 国产黄色视频一区二区在线观看| 在线观看一区二区三区激情| 九九在线视频观看精品| 在线观看三级黄色| 美女内射精品一级片tv| 欧美 日韩 精品 国产| 少妇被粗大的猛进出69影院 | 亚洲综合色惰| 美国免费a级毛片| 国产色婷婷99| 寂寞人妻少妇视频99o| 亚洲欧洲国产日韩| 最近中文字幕高清免费大全6| 晚上一个人看的免费电影| 高清视频免费观看一区二区| 两性夫妻黄色片 | 老司机影院成人| 欧美97在线视频| 久久免费观看电影| 国国产精品蜜臀av免费| 亚洲精品日本国产第一区| 久久毛片免费看一区二区三区| 男女边摸边吃奶| 在线观看国产h片| 黄网站色视频无遮挡免费观看| 日本欧美视频一区| 五月开心婷婷网| 91aial.com中文字幕在线观看| 亚洲情色 制服丝袜| 亚洲精品成人av观看孕妇| av视频免费观看在线观看| 免费观看av网站的网址| 亚洲国产av影院在线观看| 国产免费视频播放在线视频| 婷婷色av中文字幕| 亚洲,欧美精品.| 伊人亚洲综合成人网| 在线观看三级黄色| 日韩 亚洲 欧美在线| 国产一区有黄有色的免费视频| 成人免费观看视频高清| 免费播放大片免费观看视频在线观看| av播播在线观看一区| 丰满迷人的少妇在线观看| 99久久精品国产国产毛片| 久久精品国产鲁丝片午夜精品| 国产在线免费精品| 婷婷色综合大香蕉| 国产综合精华液| 亚洲成人一二三区av| 免费看光身美女| 日本午夜av视频| 亚洲精品久久午夜乱码| 午夜精品国产一区二区电影| 亚洲在久久综合| 亚洲熟女精品中文字幕| 少妇猛男粗大的猛烈进出视频| 久久久久久久国产电影| 久久精品国产亚洲av涩爱| 永久免费av网站大全| 久久久精品94久久精品| 欧美人与性动交α欧美软件 | 1024视频免费在线观看| 丝袜美足系列| 久久久久久人人人人人| 美女脱内裤让男人舔精品视频| 亚洲精品乱久久久久久| av不卡在线播放| 国产欧美日韩一区二区三区在线| 亚洲中文av在线| 免费人妻精品一区二区三区视频| 成人影院久久| 人妻系列 视频| 日本黄色日本黄色录像| 一区在线观看完整版| 美女内射精品一级片tv| 女的被弄到高潮叫床怎么办| 亚洲国产精品999| freevideosex欧美| 人人澡人人妻人| 男人舔女人的私密视频| 精品国产一区二区三区久久久樱花| 成人国产av品久久久| 满18在线观看网站| 亚洲av福利一区| 最近2019中文字幕mv第一页| 久久久久久久久久成人| 观看av在线不卡| 好男人视频免费观看在线| 十八禁高潮呻吟视频| 午夜福利在线观看免费完整高清在| 免费大片18禁| 99热这里只有是精品在线观看| 午夜免费男女啪啪视频观看| 成人毛片60女人毛片免费| 少妇的丰满在线观看| 中文乱码字字幕精品一区二区三区| 亚洲丝袜综合中文字幕| 蜜臀久久99精品久久宅男| 中文字幕av电影在线播放| 亚洲少妇的诱惑av| 国产成人精品久久久久久| 国产亚洲午夜精品一区二区久久| 国产片特级美女逼逼视频| 18+在线观看网站| 日韩中文字幕视频在线看片| 免费黄色在线免费观看| 日本黄色日本黄色录像| 久久人人爽人人爽人人片va| 综合色丁香网| 韩国高清视频一区二区三区| 免费在线观看完整版高清| freevideosex欧美| 人人澡人人妻人| 国产成人欧美| 亚洲成色77777| 亚洲国产欧美在线一区| 亚洲欧美中文字幕日韩二区| 久久狼人影院| 如日韩欧美国产精品一区二区三区| 日日爽夜夜爽网站| 国产日韩欧美亚洲二区| 久久99一区二区三区| 丰满乱子伦码专区| 欧美成人精品欧美一级黄| 99视频精品全部免费 在线| 国产精品欧美亚洲77777| 亚洲精品国产av蜜桃| 制服丝袜香蕉在线| 人人澡人人妻人| 高清在线视频一区二区三区| 男女边摸边吃奶| 精品卡一卡二卡四卡免费| 久久毛片免费看一区二区三区| 免费播放大片免费观看视频在线观看| 国产免费一区二区三区四区乱码| 国产免费又黄又爽又色| 国产极品天堂在线| 丝袜美足系列| 少妇的逼水好多| 少妇人妻 视频| 美国免费a级毛片| 18+在线观看网站| 国产 一区精品| 国产片内射在线| 国产精品久久久久久av不卡| 国产在视频线精品| 毛片一级片免费看久久久久| 在线观看美女被高潮喷水网站| 午夜久久久在线观看| 9色porny在线观看| av在线老鸭窝| 中文欧美无线码| 老司机影院成人| 国产成人精品福利久久| 一边亲一边摸免费视频| 国产精品国产三级国产专区5o| 99精国产麻豆久久婷婷| 青春草亚洲视频在线观看| 伊人久久国产一区二区| 亚洲av欧美aⅴ国产| 国产在视频线精品| 赤兔流量卡办理| 99热网站在线观看| 丰满饥渴人妻一区二区三| 日本免费在线观看一区| 中国国产av一级| 新久久久久国产一级毛片| 超色免费av| 中文字幕精品免费在线观看视频 | 免费在线观看黄色视频的| 老司机亚洲免费影院| 亚洲在久久综合| 久久久精品94久久精品| 国产乱人偷精品视频| 欧美日韩成人在线一区二区| 亚洲欧美清纯卡通| 黄色配什么色好看| 国产男人的电影天堂91| 亚洲美女搞黄在线观看| 1024视频免费在线观看| videossex国产| 欧美激情 高清一区二区三区| 亚洲精品乱久久久久久| 久热久热在线精品观看| 中文天堂在线官网| 久久精品国产综合久久久 | 男女边吃奶边做爰视频| 青春草国产在线视频| 男女啪啪激烈高潮av片| 最后的刺客免费高清国语| kizo精华| 日韩中文字幕视频在线看片| 国产爽快片一区二区三区| 国产精品 国内视频| 日韩欧美精品免费久久| 久久久久久久久久成人| 九色成人免费人妻av| 韩国精品一区二区三区 | 男人爽女人下面视频在线观看| 亚洲av电影在线进入| 国产视频首页在线观看| 国产成人精品一,二区| 极品人妻少妇av视频| 成年av动漫网址| 免费观看a级毛片全部| 欧美成人午夜精品| 纵有疾风起免费观看全集完整版| 美女福利国产在线| 久久精品国产亚洲av涩爱| 多毛熟女@视频| 精品一区二区三卡| 一级毛片 在线播放| 中文字幕免费在线视频6| av播播在线观看一区| 最近的中文字幕免费完整| 国产av码专区亚洲av| 久久综合国产亚洲精品| 久久久精品免费免费高清| 少妇的丰满在线观看| 在线亚洲精品国产二区图片欧美| freevideosex欧美| 在线 av 中文字幕| 国产免费现黄频在线看| 久久青草综合色| 少妇被粗大的猛进出69影院 | 人人澡人人妻人| 成年动漫av网址| 久久久久久伊人网av| 精品99又大又爽又粗少妇毛片| 久久99蜜桃精品久久| 99热全是精品| 亚洲av中文av极速乱| 精品久久久精品久久久| 久久久久精品人妻al黑| 在线观看国产h片| 人妻一区二区av| 国产精品欧美亚洲77777| 男人添女人高潮全过程视频| 国产一区亚洲一区在线观看| 天天影视国产精品| 2018国产大陆天天弄谢| 久久这里有精品视频免费| 亚洲精品国产色婷婷电影| 宅男免费午夜| 天堂8中文在线网| 亚洲在久久综合| 日日爽夜夜爽网站| 国产高清不卡午夜福利| 最近手机中文字幕大全| 寂寞人妻少妇视频99o| www日本在线高清视频| 人成视频在线观看免费观看| 丝袜脚勾引网站| 18在线观看网站| 亚洲精品成人av观看孕妇| 如日韩欧美国产精品一区二区三区| 飞空精品影院首页| 国产成人免费观看mmmm| 欧美97在线视频| 中文字幕最新亚洲高清| 久久久久久久久久成人|