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

    汶川震區(qū)北川縣泥石流流域崩滑體時(shí)空演變特征

    2020-06-02 11:09:32李明威史青云
    水文地質(zhì)工程地質(zhì) 2020年3期
    關(guān)鍵詞:研究

    李明威,唐 川,陳 明,史青云

    (成都理工大學(xué)地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護(hù)國(guó)家重點(diǎn)實(shí)驗(yàn)室,四川 成都 610059)

    強(qiáng)震不僅直接引發(fā)大量同震崩塌、滑坡,還會(huì)在震后較長(zhǎng)時(shí)間內(nèi)提高山體斜坡不穩(wěn)定性,崩塌、滑坡堆積體是后續(xù)降雨誘發(fā)泥石流的主要固體物質(zhì)來(lái)源[1-3]。1999 年臺(tái)灣集集地震發(fā)生后,同震滑坡面積高達(dá)159.77 km2[4],在2001年7月30日的Thphoon Toraji臺(tái)風(fēng)事件中該地區(qū)暴發(fā)了嚴(yán)重的泥石流災(zāi)害[5]。2008年“5.12”汶川地震誘發(fā)的同震崩塌、滑坡所形成的固體松散物總體積達(dá)52.5×108m3[6]。在震后幾年內(nèi),汶川震區(qū)暴發(fā)了大規(guī)模群發(fā)性泥石流災(zāi)害。例如,2008年9月24日,持續(xù)性強(qiáng)降雨導(dǎo)致北川縣城暴發(fā)以魏家溝為主的群發(fā)性泥石流災(zāi)害[7];2010年汶川縣映秀鎮(zhèn)、都江堰市龍池鎮(zhèn)、綿竹市在短時(shí)強(qiáng)降雨作用下暴發(fā)了特大泥石流群災(zāi)害事件[8-11];2013年汶川縣綿虒鎮(zhèn)、佛堂壩村在持續(xù)性暴雨作用下發(fā)生了泥石流災(zāi)害[12-13]。因此,汶川震后滑坡泥石流活動(dòng)趨勢(shì)的預(yù)測(cè)研究極為迫切。

    汶川地震后滑坡泥石流活動(dòng)研究中,眾多學(xué)者主要集中于對(duì)同震滑坡以及泥石流發(fā)育特征、形成機(jī)制、分布規(guī)律的研究[14-16],僅有少數(shù)學(xué)者對(duì)震后滑坡、泥石流活動(dòng)演化趨勢(shì)展開(kāi)研究。常鳴等[14]綜合選取震前及震后幾年,特別是大規(guī)模群發(fā)性泥石流災(zāi)害暴發(fā)前后的多期遙感影像,對(duì)比分析不同時(shí)期崩滑體的空間分布規(guī)律及演變規(guī)律。Tang等[18]以北川縣城為例,選取兩期遙感影像,解譯了汶川地震及隨后的強(qiáng)降雨事件誘發(fā)的滑坡災(zāi)害,并分析了暴雨作用對(duì)于泥石流的影響,預(yù)測(cè)汶川震區(qū)地質(zhì)災(zāi)害在震后10~15 a內(nèi)仍將保持高活動(dòng)性。黃潤(rùn)秋[19]基于多年災(zāi)害頻數(shù)統(tǒng)計(jì)數(shù)據(jù),結(jié)合關(guān)東大地震與臺(tái)灣集集地震的災(zāi)害活動(dòng)演變趨勢(shì),預(yù)測(cè)汶川震區(qū)地質(zhì)災(zāi)害活動(dòng)將持續(xù)20~25 a,并以4~5 a為一個(gè)周期震蕩式衰減。Tang等[20]利用龍池地區(qū)多期遙感影像解譯結(jié)果,對(duì)崩滑體活動(dòng)性進(jìn)行分級(jí),并基于此對(duì)未來(lái)崩滑體演變情況進(jìn)行了預(yù)測(cè)和分析,認(rèn)為震后10 a研究區(qū)滑坡泥石流活動(dòng)已趨于微弱。

    滑坡泥石流的活動(dòng)與流域植被恢復(fù)有著相當(dāng)密切的關(guān)系,植被恢復(fù)情況是影響崩滑體活動(dòng)性的關(guān)鍵因素之一,植被恢復(fù)對(duì)震后崩滑堆積體穩(wěn)定性的影響研究顯得極為關(guān)鍵。早期,一些學(xué)者對(duì)震區(qū)植被動(dòng)態(tài)變化展開(kāi)了少量研究,基于多期NDVI的監(jiān)測(cè)結(jié)果,動(dòng)態(tài)分析了汶川震區(qū)植被的時(shí)空變化特征及恢復(fù)情況[21-22],但他們對(duì)于汶川震區(qū)崩滑體的研究多集中在震后短期時(shí)間內(nèi)的對(duì)比分析上,對(duì)崩滑體長(zhǎng)期演變特征的分析較少,同時(shí),對(duì)于植被與泥石流活動(dòng)性的關(guān)系研究也較為欠缺。

    1 研究區(qū)概況

    本文選擇汶川地震極重災(zāi)區(qū)北川縣城8條泥石流流域?yàn)檠芯繀^(qū)(圖1),總面積為663.9×104m2。北川縣城地處龍門(mén)山前山與后山的交界地帶,以山地為主,海拔700~1 500 m,地勢(shì)整體西高東低,由西北向東南傾斜。研究區(qū)位于映秀—北川斷裂帶上盤(pán),在地震中受影響極大,毀壞嚴(yán)重。研究區(qū)出露地層主要包括寒武系砂巖、砂頁(yè)巖、泥質(zhì)灰?guī)r,志留系板巖、千枚巖、灰?guī)r,泥盆系和石炭系的碳酸鹽巖,以及廣泛分布于河流兩側(cè)和山前溝口地帶的階地和洪積扇地帶的第四系松散堆積物。

    研究區(qū)位于亞熱帶濕潤(rùn)季風(fēng)氣候區(qū),多年平均氣溫15.6 ℃。區(qū)內(nèi)雨量充沛,年平均降水量1 399 mm,主要集中在6—9月,占全年的74%。研究區(qū)各支溝溝口及兩側(cè)山體斜坡多為林地,植被覆蓋率較高,可達(dá)80%以上?!?.12”地震時(shí)地質(zhì)災(zāi)害發(fā)生地段植被遭受?chē)?yán)重破壞。

    研究區(qū)在汶川地震前無(wú)較大范圍滑坡和崩塌分布,僅有零星小規(guī)?;l(fā)育在溝谷兩側(cè)?!?.12”汶川地震后,北川縣城周邊山體發(fā)生大面積山體滑坡,例如王家?guī)r滑坡、新北川中學(xué)崩滑體。2008年9月23—24日,北川縣遭遇汶川地震后的最大一場(chǎng)強(qiáng)降雨過(guò)程,累計(jì)降水量達(dá)331.9 mm[23],導(dǎo)致北川縣城西側(cè)的8條溝道暴發(fā)泥石流。根據(jù)綿陽(yáng)市統(tǒng)計(jì)年鑒查閱得到,北川縣歷年雨季降雨量如表1所示。

    表1 2008—2016年北川縣雨季(6—9月)平均降水量

    圖1 研究區(qū)位置及泥石流流域概況

    2 數(shù)據(jù)來(lái)源與研究方法

    汶川震區(qū)北川縣城泥石流源區(qū)崩滑體空間分布特征分析,主要利用從各種來(lái)源獲得的多期遙感圖像進(jìn)行崩滑體解譯,并利用ArcGIS軟件分析其空間分布特征。

    結(jié)合前人研究資料和相關(guān)數(shù)據(jù),綜合考慮分辨率、云層覆蓋率等相關(guān)因素,共選擇2008年“5.12”地震發(fā)生后、2008年“9.24”泥石流發(fā)生后、2010年、2011年、2013年、2014年、2015年、2016年八期遙感影像數(shù)據(jù),以通過(guò)ALOS影像下載收集得到研究區(qū)12.5 m分辨率的DEM為基準(zhǔn),對(duì)遙感影像進(jìn)行地理配準(zhǔn)與幾何校正。以“5.12”地震后的崩滑體解譯結(jié)果為基準(zhǔn),進(jìn)行震后幾年崩滑體解譯,保持滑坡邊界一致,盡可能減少圖像質(zhì)量差異造成的誤差。

    針對(duì)震后研究區(qū)內(nèi)植被恢復(fù)的研究主要基于像元二分模型[24]計(jì)算植被覆蓋度(VFC)和植被覆蓋度恢復(fù)率(VCRR)[22],其計(jì)算模型如下:

    (1)

    式中:NDVIsoil——完全無(wú)植被覆蓋區(qū)域的NDVI值;

    企業(yè)戰(zhàn)略是指企業(yè)為了實(shí)現(xiàn)長(zhǎng)期生存和發(fā)展,通過(guò)外部環(huán)境與內(nèi)部資源的整合,對(duì)公司的發(fā)展目標(biāo)、業(yè)務(wù)方向、指導(dǎo)方針、實(shí)施步驟等作出的總體性謀劃。這強(qiáng)調(diào)了企業(yè)戰(zhàn)略必須要確保企業(yè)適應(yīng)環(huán)境變化,實(shí)現(xiàn)企業(yè)資源條件與內(nèi)外部環(huán)境因素的匹配,才能運(yùn)用管理的手段為企業(yè)整體和企業(yè)有關(guān)的活動(dòng)制定行動(dòng)綱領(lǐng),從而明確企業(yè)發(fā)展的方向和目的,界定業(yè)務(wù)范圍,構(gòu)筑競(jìng)爭(zhēng)優(yōu)勢(shì)。

    NDVIveg——完全被植被所覆蓋區(qū)域的NDVI值。

    當(dāng)研究區(qū)內(nèi)近似可取VFCmax=100%,VFCmin=0%時(shí),式(1)變?yōu)椋?/p>

    (2)

    (3)

    式中:NDVIi——地震后某一年的植被指數(shù);

    NDVIpost——地震后植被覆蓋度最低的一期植被指數(shù);

    NDVIpre——地震前植被指數(shù)。

    本文選擇研究區(qū)與遙感影像對(duì)應(yīng)的八期Landsat 5/7影像,利用ENVI軟件提取其歸一化植被指數(shù)(NDVI),選取閾值為5%,利用式(2)計(jì)算研究區(qū)的植被覆蓋度,在ArcGIS軟件中利用式(3)計(jì)算植被恢復(fù)率(VCRR),并分析其恢復(fù)情況。

    本文所使用的數(shù)據(jù)見(jiàn)表2。

    3 崩滑體演化與植被恢復(fù)特征

    研究區(qū)的震后崩滑體分布如圖2所示。汶川地震發(fā)生后,研究區(qū)發(fā)生大量崩塌滑坡,導(dǎo)致坡面巖土體極為松動(dòng)。解譯的崩滑體總面積達(dá)145.9×104m2,占研究區(qū)總面積的21.9%。2008年9月24日,北川縣遭遇持續(xù)性暴雨,促使原有不穩(wěn)定巖土體進(jìn)一步滑動(dòng),崩滑體面積進(jìn)一步擴(kuò)大,同時(shí)又產(chǎn)生部分新增崩滑體。因而,“9.24”泥石流災(zāi)害發(fā)生后,研究區(qū)崩滑體總面積達(dá)到峰值,為177.3×104m2,占研究區(qū)總面積的26.7%。并在隨后幾年一直保持在較高水平。隨著時(shí)間的延長(zhǎng),研究區(qū)內(nèi)未再發(fā)生強(qiáng)震、暴雨等強(qiáng)大的外動(dòng)力作用,崩滑體活動(dòng)性下降,面積減少。同時(shí),隨著坡面植被逐漸恢復(fù),進(jìn)一步抑制崩滑體活動(dòng)。2011年,崩滑體面積有了顯著減少,并在其后保持穩(wěn)定恢復(fù)狀態(tài)。研究區(qū)歷年降雨量、植被覆蓋度及崩滑體總面積的變化趨勢(shì)如圖3所示。

    表2 遙感影像數(shù)據(jù)類(lèi)型及其特征參數(shù)

    圖2 研究區(qū)震后不同時(shí)間崩滑體分布

    圖3 震后降雨量、植被覆蓋度、崩滑體面積變化

    3.1 崩滑體在高程上的分布與演化分析

    研究區(qū)流域范圍內(nèi)高程最低為585 m,最高為1 710 m,本文以12.5 m分辨率的DEM為基礎(chǔ)數(shù)據(jù),利用ArcGIS軟件將崩滑體圖層與DEM進(jìn)行疊加分析,并對(duì)崩滑體在高程上的分布面積進(jìn)行統(tǒng)計(jì)。根據(jù)地形條件將研究區(qū)高程分為五個(gè)等值區(qū)間,差值為300 m。參照常鳴等[22]分析綿遠(yuǎn)河流域泥石流形成區(qū)崩滑體特征所使用的密度比計(jì)算方法,定義某段高程區(qū)間上崩滑體總面積與研究區(qū)內(nèi)該段高程區(qū)間總面積的比值為密度比,計(jì)算其對(duì)應(yīng)的密度比,得到不同年份崩滑體在高程上的空間分布特征及演化規(guī)律(圖4)。

    從圖4可以看出,2008年“5.12”汶川地震發(fā)生后,崩滑體在高程700~1 500 m范圍內(nèi)發(fā)育面積較大。其中,900~1 100 m范圍內(nèi)崩滑體最為發(fā)育,總面積為49.3×104m2;其次為1 100~1 300 m范圍崩滑體較為發(fā)育,總面積為35.9×104m2;這與Tang 等[18]的統(tǒng)計(jì)結(jié)果基本一致?!?.24”泥石流災(zāi)害發(fā)生后,高程900~1 100 m范圍內(nèi)崩滑體面積為61×104m2,1 100~1 300 m范圍崩滑體面積為45.1×104m2??傮w來(lái)看,崩滑體在1 300~1 500 m范圍內(nèi)密度比最大。由此說(shuō)明,在海拔較高的地區(qū),巖土體較為松散,在地震、暴雨等外界因素的誘發(fā)下,極易發(fā)生崩塌滑坡。但隨著海拔的進(jìn)一步升高,崩滑體的面積有所減少,但其密度比仍維持在較高水平。隨著時(shí)間的推移,震后幾年內(nèi)各位置的滑坡均有不同程度的恢復(fù),海拔較低地區(qū)恢復(fù)程度較差,海拔大于1 500 m地區(qū)恢復(fù)程度最高,其密度比減少最快,700~1 500 m地區(qū)崩滑體呈穩(wěn)定恢復(fù)狀態(tài)。

    圖4 崩滑體在高程上的空間分布規(guī)律

    3.2 崩滑體在坡度上的分布與演化分析

    利用ArcGIS軟件中的3D-Analysis功能,從分辨率為12.5 m的DEM中提取研究區(qū)地形坡度,根據(jù)地形條件將坡度分為五個(gè)等值區(qū)間,差值為15°。利用ArcGIS的空間分析功能對(duì)崩滑體在坡度上的分布面積進(jìn)行統(tǒng)計(jì),并計(jì)算其對(duì)應(yīng)的密度比,得到不同年份崩滑體在坡度上的空間分布特征及演化規(guī)律(圖5)。

    圖5 崩滑體在坡度上的空間分布規(guī)律

    從圖5可以看出,2008年“5.12”地震發(fā)生后,研究區(qū)崩滑體在坡度30°~45°范圍內(nèi)發(fā)育面積最大,其面積為82.6×104m2。“9.24”泥石流后,研究區(qū)崩滑體在坡度30°~45°范圍內(nèi)面積增加到100.5×104m2。在坡度較高的地方崩滑體的發(fā)育面積有所減少,但其密度比仍維持在較高水平,由此說(shuō)明,隨著坡度的增加,土體的不穩(wěn)定性提高,崩滑體更容易發(fā)生。隨著時(shí)間的推移,在沒(méi)有強(qiáng)烈外界因素的影響下,崩滑體呈穩(wěn)定恢復(fù)狀態(tài),其中在坡度大于60°范圍內(nèi)的崩滑體恢復(fù)程度最高。

    3.3 崩滑體在坡向上的分布與演化分析

    不同坡向上由于山體的朝向不同,導(dǎo)致其日照時(shí)間、降雨因素等都有顯著差異,因而其崩滑體的發(fā)育程度也有所不同。利用DEM生成研究區(qū)山體坡向,以正北方向?yàn)?°,每隔45°為一個(gè)坡向區(qū)間,共分為8個(gè)區(qū)間。運(yùn)用ArcGIS中的空間分析功能,對(duì)崩滑體在不同坡向上的分布面積進(jìn)行統(tǒng)計(jì),并計(jì)算其密度比,得到不同年份崩滑體在坡向區(qū)間內(nèi)的空間分布特征及演化規(guī)律(圖6)。

    圖6 崩滑體在坡向上的空間分布規(guī)律

    從圖6可以看出,2008年“5.12”地震發(fā)生后,崩滑體多集中發(fā)育在坡向?yàn)?°~180°的范圍內(nèi),其面積為128×104m2?!?.24”泥石流后,該區(qū)域的崩滑體面積達(dá)150.7×104m2,占崩滑體總面積的85%。在180°~360°的坡向區(qū)間,崩滑體發(fā)育較少,其中在270°~315°的坡向區(qū)間幾乎無(wú)崩滑體發(fā)育。隨著時(shí)間的推移,由于在0°~180°坡向區(qū)間日照較為充足,植被發(fā)育較好,其崩滑體恢復(fù)程度也維持在較高水平。而西側(cè)坡面的恢復(fù)程度則較低。

    3.4 崩滑體在距溝道距離上的分布與演化分析

    研究區(qū)主要由8條泥石流溝組成,泥石流物源主要來(lái)自于坡面崩滑體。參考楊濤[26]對(duì)綿虒鎮(zhèn)崩滑體在距溝道距離上的統(tǒng)計(jì)分析,結(jié)合本研究區(qū)地形條件,本文以150 m為差值,將距溝道距離等間距分為五級(jí)對(duì)崩滑體在不同面積的分布距離進(jìn)行統(tǒng)計(jì),并計(jì)算其相應(yīng)的密度比,得到不同年份崩滑體在距溝道距離上的空間分布特征和演化規(guī)律(圖7)。

    圖7 崩滑體在距溝道距離上的空間分布規(guī)律

    從圖7可以看出,2008年“5.12”地震后研究區(qū)崩滑體面積最大的區(qū)域是距溝道距離0~150 m的范圍,其面積為61.2×104m2。在“9.24”泥石流災(zāi)害發(fā)生后,該范圍內(nèi)的崩滑體總面積增長(zhǎng)到了74.6×104m2,占崩滑體總面積的42%。但密度比最大的區(qū)域是300~450 m的范圍,為26.1%,暴雨導(dǎo)致崩滑體進(jìn)一步增加后,密度比達(dá)到了32.7%。隨著時(shí)間的推移,除600~750 m范圍外,其余崩滑體均呈現(xiàn)穩(wěn)定的恢復(fù)狀態(tài)。在距離溝道600~750 m的坡面上,崩滑體在2011年到2013年間有所增加,如今已基本恢復(fù),面積很小,密度比極低。其余范圍內(nèi),崩滑體呈穩(wěn)定恢復(fù)狀態(tài)。

    3.5 泥石流源區(qū)植被覆蓋度演化分析

    根據(jù)式(1)計(jì)算得到,研究區(qū)內(nèi)植被覆蓋度(VFC)在“9.24”泥石流發(fā)生后達(dá)到最低值,為0.64(圖3)。從2010年起至2013年初,研究區(qū)內(nèi)植被覆蓋度呈緩慢恢復(fù)狀態(tài),根據(jù)圖3可以看出,2013年降雨量較之前幾年有了顯著增加,雨季過(guò)后,2014年研究區(qū)植被覆蓋度略有降低,但未發(fā)生大規(guī)模泥石流災(zāi)害,崩滑體面積未增加。2014—2016年,研究區(qū)植被恢復(fù)良好,崩滑體面積減少,泥石流活動(dòng)性降低。

    根據(jù)許積層等[21]對(duì)汶川震區(qū)岷江河谷段植被恢復(fù)特征的研究,結(jié)合本文研究區(qū)實(shí)際情況,將植被覆蓋度恢復(fù)率VCRR分為4級(jí):<0為植被恢復(fù)極差;0~0.5為植被恢復(fù)較差;0.5~1為植被恢復(fù)中等;>1為植被恢復(fù)良好。

    研究區(qū)泥石流源區(qū)崩滑體植被覆蓋度恢復(fù)率統(tǒng)計(jì)結(jié)果如圖8所示。

    圖8 植被覆蓋度恢復(fù)率

    從圖8可以看出,2008—2010年,植被恢復(fù)程度依然處于較低水平,植被恢復(fù)極差區(qū)域占比78.4%,由此說(shuō)明,在強(qiáng)烈外力破壞作用發(fā)生后的短期時(shí)間內(nèi),崩滑體活動(dòng)性仍處于較高水平,缺少植被恢復(fù)的穩(wěn)定條件。從2011年起,植被覆蓋度恢復(fù)中等以上的區(qū)域占比逐漸增高,研究區(qū)植被處于緩慢恢復(fù)狀態(tài),崩滑體活動(dòng)性逐年降低,趨于穩(wěn)定。2013年雨季降雨的顯著增加,使得植被覆蓋度恢復(fù)率有所降低,2014年植被恢復(fù)極差的區(qū)域占整個(gè)研究區(qū)面積的56.1%,但植被整體仍處于恢復(fù)狀態(tài),2015年之后植被覆蓋度恢復(fù)率中等以上區(qū)域面積進(jìn)一步增加。

    4 結(jié)論

    (1)汶川地震發(fā)生后,北川縣城崩滑體發(fā)育面積較震前顯著增加,震后不久,9月23—24日短時(shí)強(qiáng)降雨發(fā)生后,進(jìn)一步誘發(fā)大量崩滑體,導(dǎo)致“9.24”群發(fā)性泥石流災(zāi)害的發(fā)生,崩滑體面積達(dá)到峰值,并延續(xù)到2011年崩滑體面積仍處在較高水平。隨著生態(tài)和地質(zhì)環(huán)境的逐漸恢復(fù),2013年起研究區(qū)崩滑體面積顯著降低,至2016年已恢復(fù)至較低水平。

    (2)研究表明,崩滑體在高程900~1 100 m、坡度30°~45°、坡向90°~135°、距溝道150 m范圍內(nèi)發(fā)育面積最大。密度比最高的為高程1 300~1 500 m、坡度大于60°、坡向45°~90°、距溝道300~450 m范圍內(nèi)。

    (3)根據(jù)計(jì)算,研究區(qū)植被覆蓋度(VFC)在“9.24”泥石流發(fā)生后達(dá)到最低值,震后至2010年植被恢復(fù)水平較低,2011—2013年,植被恢復(fù)程度略有提高。受2013年降雨作用影響,2014年植被覆蓋度恢復(fù)率極差區(qū)域占比56.1%,隨后幾年研究區(qū)植被呈穩(wěn)定恢復(fù)狀態(tài)。

    (4)與映秀、龍池等地區(qū)相比,研究區(qū)內(nèi)2010年、2013年降雨強(qiáng)度與往常相比也有所增加,但未發(fā)生大規(guī)模泥石流群發(fā)性災(zāi)害。由此認(rèn)為,除地層巖性、微地貌等因素影響外,植被也對(duì)泥石流活動(dòng)性具有一定的抑制作用。

    猜你喜歡
    研究
    FMS與YBT相關(guān)性的實(shí)證研究
    2020年國(guó)內(nèi)翻譯研究述評(píng)
    遼代千人邑研究述論
    視錯(cuò)覺(jué)在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    關(guān)于遼朝“一國(guó)兩制”研究的回顧與思考
    EMA伺服控制系統(tǒng)研究
    基于聲、光、磁、觸摸多功能控制的研究
    電子制作(2018年11期)2018-08-04 03:26:04
    新版C-NCAP側(cè)面碰撞假人損傷研究
    關(guān)于反傾銷(xiāo)會(huì)計(jì)研究的思考
    焊接膜層脫落的攻關(guān)研究
    電子制作(2017年23期)2017-02-02 07:17:19
    免费观看a级毛片全部| 啦啦啦 在线观看视频| 曰老女人黄片| 我的亚洲天堂| 男女高潮啪啪啪动态图| 一区二区av电影网| 在线看a的网站| 国产免费又黄又爽又色| 成年女人毛片免费观看观看9 | 新久久久久国产一级毛片| 久久久精品94久久精品| 亚洲精品久久成人aⅴ小说| 亚洲精品av麻豆狂野| 天天操日日干夜夜撸| 日韩一区二区三区影片| 麻豆av在线久日| 人人妻,人人澡人人爽秒播 | 亚洲婷婷狠狠爱综合网| 欧美久久黑人一区二区| 国产激情久久老熟女| 亚洲熟女毛片儿| 亚洲四区av| 国产高清国产精品国产三级| 日韩精品有码人妻一区| 欧美日韩亚洲综合一区二区三区_| 国产不卡av网站在线观看| 国产日韩欧美在线精品| 纯流量卡能插随身wifi吗| 一边摸一边抽搐一进一出视频| 日韩免费高清中文字幕av| 黄色怎么调成土黄色| 亚洲情色 制服丝袜| 中文精品一卡2卡3卡4更新| 超色免费av| 国产免费福利视频在线观看| 五月开心婷婷网| 国产在线免费精品| 国产欧美亚洲国产| 精品人妻熟女毛片av久久网站| 美女中出高潮动态图| 日韩视频在线欧美| 日韩一区二区三区影片| 亚洲欧美清纯卡通| 久久 成人 亚洲| 啦啦啦在线免费观看视频4| 青春草亚洲视频在线观看| 男女下面插进去视频免费观看| 人人妻人人澡人人爽人人夜夜| 亚洲国产精品999| 午夜免费观看性视频| 欧美在线黄色| 午夜福利乱码中文字幕| 免费黄网站久久成人精品| 午夜福利影视在线免费观看| av电影中文网址| 在线观看免费午夜福利视频| 在线观看免费日韩欧美大片| 日韩av不卡免费在线播放| 国产 精品1| 啦啦啦中文免费视频观看日本| 十分钟在线观看高清视频www| 国产极品天堂在线| 伦理电影大哥的女人| 国产色婷婷99| 日韩 欧美 亚洲 中文字幕| 97在线人人人人妻| 纯流量卡能插随身wifi吗| 久久精品久久久久久久性| 夜夜骑夜夜射夜夜干| 咕卡用的链子| 黄色怎么调成土黄色| 午夜福利视频精品| 亚洲精品久久午夜乱码| 久久久久久久久免费视频了| 国产成人91sexporn| 亚洲av日韩在线播放| 亚洲第一青青草原| 亚洲,欧美精品.| 蜜桃在线观看..| 熟女av电影| 少妇被粗大的猛进出69影院| 久久婷婷青草| 精品国产一区二区三区四区第35| 久久综合国产亚洲精品| 久久久国产一区二区| 久久ye,这里只有精品| 人人妻人人添人人爽欧美一区卜| 日本午夜av视频| 日韩一本色道免费dvd| 国产av码专区亚洲av| 国产成人一区二区在线| 国产野战对白在线观看| 自线自在国产av| 欧美日韩亚洲高清精品| 999精品在线视频| 悠悠久久av| 自线自在国产av| 嫩草影视91久久| 男人舔女人的私密视频| 国产精品蜜桃在线观看| 精品少妇内射三级| 多毛熟女@视频| 一边亲一边摸免费视频| 久久精品人人爽人人爽视色| 日韩,欧美,国产一区二区三区| 妹子高潮喷水视频| 亚洲国产精品成人久久小说| 亚洲国产中文字幕在线视频| 最近最新中文字幕免费大全7| 免费女性裸体啪啪无遮挡网站| 岛国毛片在线播放| 精品福利永久在线观看| 尾随美女入室| 日本欧美视频一区| 另类精品久久| 天堂俺去俺来也www色官网| 一区二区三区激情视频| 丝袜脚勾引网站| 日韩 亚洲 欧美在线| 777久久人妻少妇嫩草av网站| 天天躁夜夜躁狠狠久久av| 久久精品人人爽人人爽视色| 不卡视频在线观看欧美| 欧美成人精品欧美一级黄| 亚洲国产欧美网| 亚洲精品av麻豆狂野| 丝袜美足系列| 国产精品嫩草影院av在线观看| 亚洲国产中文字幕在线视频| 中文字幕最新亚洲高清| 精品少妇久久久久久888优播| 国产免费一区二区三区四区乱码| 最近2019中文字幕mv第一页| 欧美乱码精品一区二区三区| 黄色视频在线播放观看不卡| 国产精品一二三区在线看| 少妇 在线观看| 国产精品国产三级国产专区5o| 2021少妇久久久久久久久久久| 午夜日韩欧美国产| 国产成人a∨麻豆精品| 成人影院久久| 国产亚洲最大av| 一本—道久久a久久精品蜜桃钙片| 亚洲欧洲国产日韩| 午夜日本视频在线| 国产精品国产三级专区第一集| 你懂的网址亚洲精品在线观看| 精品酒店卫生间| 老汉色av国产亚洲站长工具| 99re6热这里在线精品视频| 免费看不卡的av| 18禁国产床啪视频网站| 中文字幕人妻丝袜制服| 一区在线观看完整版| 久久久精品94久久精品| 美女福利国产在线| 婷婷成人精品国产| 国产精品久久久人人做人人爽| 自线自在国产av| 亚洲欧美成人精品一区二区| 黑丝袜美女国产一区| 免费高清在线观看视频在线观看| 亚洲激情五月婷婷啪啪| 性少妇av在线| 亚洲 欧美一区二区三区| 曰老女人黄片| 高清在线视频一区二区三区| 赤兔流量卡办理| 欧美精品高潮呻吟av久久| 久久精品熟女亚洲av麻豆精品| 青春草视频在线免费观看| 人人妻人人添人人爽欧美一区卜| 亚洲成人手机| 丝袜脚勾引网站| 女的被弄到高潮叫床怎么办| 亚洲精品在线美女| 亚洲精品中文字幕在线视频| 黑人巨大精品欧美一区二区蜜桃| av线在线观看网站| h视频一区二区三区| 不卡视频在线观看欧美| 99精国产麻豆久久婷婷| 国产成人精品久久二区二区91 | 香蕉国产在线看| 亚洲免费av在线视频| 国产日韩一区二区三区精品不卡| 久久人人爽av亚洲精品天堂| 国产免费视频播放在线视频| 高清欧美精品videossex| 韩国av在线不卡| 99热全是精品| 狂野欧美激情性bbbbbb| 国产99久久九九免费精品| 欧美日韩亚洲国产一区二区在线观看 | 十八禁网站网址无遮挡| 日本91视频免费播放| 国产精品蜜桃在线观看| 亚洲自偷自拍图片 自拍| 人人澡人人妻人| 丰满迷人的少妇在线观看| 桃花免费在线播放| 亚洲天堂av无毛| 国产黄色免费在线视频| 国产精品秋霞免费鲁丝片| 久久国产亚洲av麻豆专区| 国产在视频线精品| 日韩免费高清中文字幕av| 婷婷成人精品国产| 成年动漫av网址| 亚洲中文av在线| av在线观看视频网站免费| 国产一区二区在线观看av| 国产精品三级大全| 五月天丁香电影| 国产爽快片一区二区三区| 宅男免费午夜| 日韩av在线免费看完整版不卡| 最近中文字幕2019免费版| 丝袜人妻中文字幕| 交换朋友夫妻互换小说| 少妇人妻 视频| 日韩成人av中文字幕在线观看| 亚洲综合色网址| 天堂8中文在线网| 18禁观看日本| 日本欧美视频一区| 国产精品偷伦视频观看了| 一区二区三区四区激情视频| 美女视频免费永久观看网站| 天天添夜夜摸| 久久久精品94久久精品| 天天躁日日躁夜夜躁夜夜| 美女中出高潮动态图| 大香蕉久久网| 亚洲成色77777| 午夜福利免费观看在线| 操美女的视频在线观看| 国产精品熟女久久久久浪| 久久精品aⅴ一区二区三区四区| 亚洲精品自拍成人| 午夜久久久在线观看| 建设人人有责人人尽责人人享有的| kizo精华| 免费观看a级毛片全部| 91精品伊人久久大香线蕉| 高清在线视频一区二区三区| 大话2 男鬼变身卡| 爱豆传媒免费全集在线观看| 2021少妇久久久久久久久久久| 亚洲国产毛片av蜜桃av| 色网站视频免费| 桃花免费在线播放| 久久久久久久精品精品| 国产精品久久久久成人av| 男女无遮挡免费网站观看| 伊人亚洲综合成人网| 亚洲av福利一区| 国产激情久久老熟女| 侵犯人妻中文字幕一二三四区| 在现免费观看毛片| 久久精品熟女亚洲av麻豆精品| 国产在线免费精品| 十分钟在线观看高清视频www| av国产精品久久久久影院| 一二三四在线观看免费中文在| 午夜免费观看性视频| 成人手机av| 啦啦啦视频在线资源免费观看| 十八禁高潮呻吟视频| 人成视频在线观看免费观看| 一本色道久久久久久精品综合| 亚洲精品国产av成人精品| 久久久久久人妻| 中文字幕另类日韩欧美亚洲嫩草| 成人三级做爰电影| 激情视频va一区二区三区| 90打野战视频偷拍视频| 国产精品免费大片| 女性被躁到高潮视频| 天天影视国产精品| 亚洲色图 男人天堂 中文字幕| 国产免费一区二区三区四区乱码| 免费少妇av软件| 亚洲,欧美,日韩| 亚洲天堂av无毛| 亚洲综合精品二区| 老司机影院成人| 成人午夜精彩视频在线观看| 亚洲精品视频女| 欧美日韩亚洲国产一区二区在线观看 | 男人舔女人的私密视频| 制服丝袜香蕉在线| 成人三级做爰电影| 久久久精品区二区三区| 亚洲成人av在线免费| 欧美中文综合在线视频| 天堂俺去俺来也www色官网| www.av在线官网国产| 女的被弄到高潮叫床怎么办| 欧美人与性动交α欧美软件| 欧美老熟妇乱子伦牲交| 黄片小视频在线播放| 黄片无遮挡物在线观看| 在线观看国产h片| 久久久久久免费高清国产稀缺| 欧美黄色片欧美黄色片| 亚洲免费av在线视频| 水蜜桃什么品种好| 欧美变态另类bdsm刘玥| 一级片免费观看大全| 国产免费视频播放在线视频| 一区福利在线观看| 大香蕉久久成人网| 精品少妇久久久久久888优播| 国产精品蜜桃在线观看| 精品酒店卫生间| 午夜激情av网站| 中文字幕人妻丝袜一区二区 | 亚洲成人一二三区av| 成人午夜精彩视频在线观看| 最近最新中文字幕大全免费视频 | 国产成人精品无人区| 大话2 男鬼变身卡| av又黄又爽大尺度在线免费看| av电影中文网址| 欧美变态另类bdsm刘玥| 日本色播在线视频| 日本av手机在线免费观看| 五月天丁香电影| 久久精品熟女亚洲av麻豆精品| 欧美精品人与动牲交sv欧美| 欧美精品av麻豆av| 热re99久久精品国产66热6| 看十八女毛片水多多多| 可以免费在线观看a视频的电影网站 | av电影中文网址| 国产一区有黄有色的免费视频| 大陆偷拍与自拍| 国产精品 欧美亚洲| 欧美成人午夜精品| 国产一区二区三区综合在线观看| 久久久精品国产亚洲av高清涩受| 国语对白做爰xxxⅹ性视频网站| 日本欧美视频一区| 如日韩欧美国产精品一区二区三区| 成人国语在线视频| 麻豆精品久久久久久蜜桃| 巨乳人妻的诱惑在线观看| 91国产中文字幕| 免费在线观看视频国产中文字幕亚洲 | 伦理电影大哥的女人| 日韩中文字幕视频在线看片| 另类精品久久| 香蕉丝袜av| 国产精品一区二区在线不卡| 亚洲五月色婷婷综合| 成人亚洲精品一区在线观看| √禁漫天堂资源中文www| 狠狠精品人妻久久久久久综合| 久久精品国产亚洲av高清一级| 亚洲欧洲国产日韩| 在线免费观看不下载黄p国产| 亚洲成av片中文字幕在线观看| 波多野结衣一区麻豆| 久久精品亚洲熟妇少妇任你| 19禁男女啪啪无遮挡网站| 色婷婷av一区二区三区视频| 青春草国产在线视频| 国产午夜精品一二区理论片| 制服诱惑二区| 久久热在线av| 国产精品麻豆人妻色哟哟久久| 欧美xxⅹ黑人| 一区二区三区四区激情视频| 亚洲av成人不卡在线观看播放网 | 9191精品国产免费久久| 久久久久国产精品人妻一区二区| 亚洲精品美女久久久久99蜜臀 | 男女边摸边吃奶| 观看美女的网站| 97人妻天天添夜夜摸| 精品亚洲成国产av| 人成视频在线观看免费观看| 在线观看国产h片| 国产乱人偷精品视频| 一区二区日韩欧美中文字幕| 午夜日本视频在线| 久久97久久精品| 国产xxxxx性猛交| 波多野结衣一区麻豆| 成人漫画全彩无遮挡| 侵犯人妻中文字幕一二三四区| 国产又色又爽无遮挡免| 日韩大码丰满熟妇| 美国免费a级毛片| 在线观看免费视频网站a站| 欧美日韩福利视频一区二区| 五月天丁香电影| 老汉色av国产亚洲站长工具| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品久久午夜乱码| 亚洲自偷自拍图片 自拍| 看十八女毛片水多多多| 天美传媒精品一区二区| 国产日韩欧美在线精品| 亚洲熟女精品中文字幕| 国产亚洲av高清不卡| 欧美xxⅹ黑人| 国产成人免费无遮挡视频| 久久天堂一区二区三区四区| 国产精品人妻久久久影院| 亚洲国产欧美日韩在线播放| 日韩av不卡免费在线播放| 亚洲国产欧美网| 亚洲,欧美精品.| 国产成人啪精品午夜网站| 女的被弄到高潮叫床怎么办| 纵有疾风起免费观看全集完整版| 美女扒开内裤让男人捅视频| 乱人伦中国视频| 青草久久国产| 亚洲精品日本国产第一区| 女性生殖器流出的白浆| 久久午夜综合久久蜜桃| 午夜福利乱码中文字幕| 亚洲综合色网址| 精品亚洲成a人片在线观看| 永久免费av网站大全| 亚洲免费av在线视频| av国产精品久久久久影院| 性色av一级| 夫妻性生交免费视频一级片| 18在线观看网站| 国产精品三级大全| 亚洲精品成人av观看孕妇| 国产日韩欧美在线精品| 国产av码专区亚洲av| 亚洲伊人久久精品综合| 亚洲自偷自拍图片 自拍| 日韩欧美一区视频在线观看| 一区二区三区精品91| 啦啦啦中文免费视频观看日本| 欧美精品一区二区大全| 亚洲国产中文字幕在线视频| 精品久久久精品久久久| 晚上一个人看的免费电影| videos熟女内射| 国产成人精品无人区| 伦理电影免费视频| 亚洲精品一区蜜桃| 不卡视频在线观看欧美| 九草在线视频观看| 久久精品久久久久久噜噜老黄| 久久久精品区二区三区| 国产极品天堂在线| 最近中文字幕高清免费大全6| 伦理电影大哥的女人| 色婷婷av一区二区三区视频| 五月天丁香电影| 看免费成人av毛片| 成年女人毛片免费观看观看9 | av不卡在线播放| 亚洲成人国产一区在线观看 | 色精品久久人妻99蜜桃| 色播在线永久视频| 1024香蕉在线观看| 一边亲一边摸免费视频| 亚洲五月色婷婷综合| 国产伦人伦偷精品视频| 青青草视频在线视频观看| 精品亚洲乱码少妇综合久久| 只有这里有精品99| 国语对白做爰xxxⅹ性视频网站| 一区二区日韩欧美中文字幕| 亚洲国产av影院在线观看| 亚洲欧美激情在线| 欧美xxⅹ黑人| 国产成人免费无遮挡视频| 国产又色又爽无遮挡免| 久久天躁狠狠躁夜夜2o2o | 日日爽夜夜爽网站| 午夜免费男女啪啪视频观看| 国产极品粉嫩免费观看在线| 秋霞在线观看毛片| 天天添夜夜摸| 国精品久久久久久国模美| 久久精品aⅴ一区二区三区四区| 国产熟女午夜一区二区三区| 777久久人妻少妇嫩草av网站| 精品酒店卫生间| 国产亚洲午夜精品一区二区久久| 97精品久久久久久久久久精品| 久久久国产精品麻豆| 在线观看国产h片| 国产色婷婷99| 亚洲伊人久久精品综合| 亚洲av综合色区一区| 日韩av不卡免费在线播放| 精品卡一卡二卡四卡免费| 亚洲男人天堂网一区| www.精华液| 国产精品麻豆人妻色哟哟久久| 又大又黄又爽视频免费| 亚洲成人一二三区av| 这个男人来自地球电影免费观看 | 一本一本久久a久久精品综合妖精| 极品人妻少妇av视频| 建设人人有责人人尽责人人享有的| 叶爱在线成人免费视频播放| 日日摸夜夜添夜夜爱| 欧美变态另类bdsm刘玥| 免费不卡黄色视频| 精品国产超薄肉色丝袜足j| 国产欧美亚洲国产| 嫩草影院入口| 美女脱内裤让男人舔精品视频| 免费久久久久久久精品成人欧美视频| 女人高潮潮喷娇喘18禁视频| 一级毛片 在线播放| 波多野结衣av一区二区av| 亚洲免费av在线视频| 欧美av亚洲av综合av国产av | 嫩草影视91久久| 69精品国产乱码久久久| 成年av动漫网址| 国产精品一国产av| 在线观看www视频免费| 国产精品成人在线| 国产一级毛片在线| 99香蕉大伊视频| 伊人久久国产一区二区| 亚洲免费av在线视频| 久久狼人影院| 国产精品欧美亚洲77777| 亚洲,欧美精品.| 最近中文字幕高清免费大全6| 久久婷婷青草| 丰满迷人的少妇在线观看| 国产一区二区三区综合在线观看| 欧美日韩亚洲国产一区二区在线观看 | 久久这里只有精品19| 国产精品麻豆人妻色哟哟久久| 国产成人免费观看mmmm| av线在线观看网站| 日韩伦理黄色片| 黑人巨大精品欧美一区二区蜜桃| 日本色播在线视频| 国产黄色视频一区二区在线观看| 天天操日日干夜夜撸| 亚洲精品成人av观看孕妇| 免费观看a级毛片全部| 国产无遮挡羞羞视频在线观看| 制服诱惑二区| 成年av动漫网址| 性少妇av在线| 91成人精品电影| 美女扒开内裤让男人捅视频| 成人午夜精彩视频在线观看| 少妇人妻 视频| 精品亚洲成国产av| 多毛熟女@视频| 老汉色av国产亚洲站长工具| 一区二区三区四区激情视频| www.av在线官网国产| 亚洲精品av麻豆狂野| 夫妻性生交免费视频一级片| 综合色丁香网| 99国产精品免费福利视频| 卡戴珊不雅视频在线播放| 成人免费观看视频高清| 色播在线永久视频| av.在线天堂| 无限看片的www在线观看| 国产成人精品久久二区二区91 | 建设人人有责人人尽责人人享有的| 国产在线免费精品| 欧美中文综合在线视频| 性少妇av在线| 新久久久久国产一级毛片| 十分钟在线观看高清视频www| 亚洲成人一二三区av| 91国产中文字幕| 妹子高潮喷水视频| 操出白浆在线播放| 亚洲av成人不卡在线观看播放网 | 九九爱精品视频在线观看| 国产 一区精品| 国产乱来视频区| 亚洲av国产av综合av卡| 亚洲精品日本国产第一区| 伊人亚洲综合成人网| 亚洲欧美日韩另类电影网站| 午夜日本视频在线| 久久人人爽人人片av| 99精国产麻豆久久婷婷| 在线观看免费高清a一片| 夫妻性生交免费视频一级片| 夜夜骑夜夜射夜夜干| 亚洲国产欧美一区二区综合| 久久精品久久久久久久性| 青草久久国产| 99久久人妻综合| 爱豆传媒免费全集在线观看| 一区福利在线观看| 丰满迷人的少妇在线观看| 熟女av电影| 日韩中文字幕视频在线看片| 亚洲国产欧美日韩在线播放| 视频在线观看一区二区三区| 亚洲一区中文字幕在线| 免费观看人在逋| 中文字幕高清在线视频| 91精品三级在线观看| 精品一区在线观看国产|