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

    云南鎮(zhèn)雄縣滑坡主要致災(zāi)因子篩選及其回歸模型的建立

    2015-03-01 09:13:30談樹成龍叔林楊軍楊煬
    關(guān)鍵詞:鎮(zhèn)雄縣巖組坡向

    談樹成,龍叔林,楊軍,楊煬

    (云南大學(xué)資源環(huán)境與地球科學(xué)學(xué)院,昆明 650091)

    滑坡是一個(gè)復(fù)雜的,由各種因素相互作用及影響構(gòu)成的非線性系統(tǒng)系統(tǒng),是內(nèi)因(環(huán)境因子)和外因(觸發(fā)因子)共同作用的結(jié)果。鎮(zhèn)雄縣滑坡地質(zhì)災(zāi)害頻發(fā),制約著該縣的經(jīng)濟(jì)發(fā)展、威脅公路暢通和人民財(cái)產(chǎn)安全。鎮(zhèn)雄縣滑坡地質(zhì)災(zāi)害發(fā)育主要受地形地貌、地層巖性、地質(zhì)構(gòu)造、水文地質(zhì)等基本地質(zhì)條件及人類工程活動(dòng)、大氣降雨等誘發(fā)因素的控制[1]。

    1 研究區(qū)概況及地質(zhì)災(zāi)害情況

    鎮(zhèn)雄縣位于云南省東北部,縣境內(nèi)重巒疊嶂,峰高谷深,全縣無(wú)壩子,半山區(qū)、山區(qū)、高山區(qū)分別占全縣總面積的29%、46%、25%。部分地段為懸崖陡坡,受地形、構(gòu)造、巖性、氣候的控制,加之,隨著社會(huì)經(jīng)濟(jì)的發(fā)展及人口資源的不斷增加,人類工程活動(dòng)頻繁,坡體植被破壞嚴(yán)重,水土流失較重,山坡物質(zhì)極易向山下轉(zhuǎn)移,屬地質(zhì)災(zāi)害易發(fā)的地區(qū)。鎮(zhèn)雄縣已發(fā)滑坡地質(zhì)災(zāi)害111個(gè),屬于滑坡地質(zhì)災(zāi)害多發(fā)地區(qū),滑坡的頻繁發(fā)生,給鎮(zhèn)雄縣人民的生命財(cái)產(chǎn)帶來(lái)了巨大的損失。

    2 滑坡地質(zhì)災(zāi)害致災(zāi)因子敏感性分析

    本文所采取的數(shù)據(jù)量化及分析方法為滑坡確定性系數(shù)(CF)分析方法[2-4],該方法是一種結(jié)合確定性系數(shù)(CF)用來(lái)分析影響滑坡發(fā)生各個(gè)因素敏感性的方法?;麓_定性系數(shù)(CF)分析方法為一個(gè)概率函數(shù),最早由E.H.Shortliffe和B.G.Buchanan[5]提出,由 D.Hecker man[6]進(jìn)行完善改進(jìn),具體表示為:

    式中,PPa為事件(滑坡)在數(shù)據(jù)類a中發(fā)生的條件概率,在實(shí)際滑坡應(yīng)用時(shí)可以表示為代表數(shù)據(jù)類a的單元中存在的滑坡面積與單元面積的比值:PPs為事件在整個(gè)研究區(qū)A中發(fā)生的先驗(yàn)概率,可以表示為整個(gè)研究區(qū)的滑坡面積與研究區(qū)面積的比值。

    通過式(1)的函數(shù)轉(zhuǎn)化,CF 為[-1,l]。正值代表事件發(fā)生確定性的增長(zhǎng),即滑坡變形失穩(wěn)的確定性高;負(fù)值代表確定性的降低,表示滑坡變形失穩(wěn)的確定性低;接近于0值,代表先驗(yàn)概率與條件概率十分接近,時(shí)間發(fā)生的確定性不可能進(jìn)行確定。通過此方法能確定每個(gè)一級(jí)因子中較為重要的二級(jí)因子范圍,即每個(gè)二級(jí)因子的敏感性。

    鎮(zhèn)雄縣滑坡災(zāi)害分析選取的影響因子主要包括:高程、地形起伏度、巖性、坡度、坡向、與河流距離、與斷層距離、與公路距離等8個(gè)主要因子。對(duì)每個(gè)因子進(jìn)行了分組,具體結(jié)果見表1。

    2.1 高程

    鎮(zhèn)雄縣地處烏蒙山脈北緣斜坡山地,處于金沙江南岸,總體地勢(shì)西南高、東北低,全縣最低高程520 m,最高高程2 416 m,縣境內(nèi)最大高差達(dá)1 896 m。

    表1 致災(zāi)因子分類表

    通過 (1)式的函數(shù)轉(zhuǎn)化,得到表2高程分級(jí)及CF值。

    由表2與圖1可知高程因子中700~1 200 m滑坡易發(fā)性高,即高程因子中700~1 200 m敏感性高。

    2.2 坡度

    坡度在斜坡穩(wěn)定性分析中,是一個(gè)非常重要的要素。在平緩地帶,滑坡是很難形成的。斜坡失穩(wěn)形成滑坡一般發(fā)生在<60°的地帶,超過60°,斜坡失穩(wěn)一般以崩塌的形式出現(xiàn)。

    鎮(zhèn)雄縣屬于山區(qū),坡度變化明顯,高坡度地區(qū)主要集中在鎮(zhèn)雄縣中東部地區(qū),坡度高達(dá)50°,中部及中西部相對(duì)平緩。

    通過(1)式的函數(shù)轉(zhuǎn)化,得到表2坡度分級(jí)及CF值。由表2與圖2可知坡度因子中5°~15°滑坡易發(fā)性高,即高程因子中5°~15°敏感性高。

    2.3 起伏度

    起伏度為特定區(qū)域內(nèi)(30 m×30 m柵格單元內(nèi)),最高高程與最低高程的差值,及相對(duì)高差。本文為了方便后期數(shù)據(jù)分析,引入起伏度的概念對(duì)地形變化進(jìn)行量化。鎮(zhèn)雄縣境內(nèi)山高坡陡,溝谷縱橫,主要為山區(qū),面積占縣域總面積的46%。

    通過上述(1)式的函數(shù)轉(zhuǎn)化,得到表2起伏度分級(jí)及CF值。由表2與圖3可知起伏度因子中30~200 m滑坡易發(fā)性高,即起伏度因子中30~200 m敏感性高。

    2.4 巖性(巖組)

    鎮(zhèn)雄縣地層出露較齊全,除第三系、泥盆系外,其他地層均有出露,以二疊系、三疊系較為發(fā)育,且上二疊統(tǒng)峨眉山玄武巖分布廣泛。

    表2 滑坡致災(zāi)因子CF值表

    圖1 鎮(zhèn)雄縣高程與滑坡分布圖

    圖2 鎮(zhèn)雄縣坡度與滑坡分布圖

    巖體總體以沉積巖為主,主要巖性為碳酸鹽巖、碎屑巖,其次是濱海相、海陸相交互相的砂巖、泥巖、頁(yè)巖和含煤建造,巖漿巖比較單一,僅出露二疊系上統(tǒng)峨眉山組(P2β)玄武巖。第四系沉積類型單一,山間小盆谷堆積、偶見河流階地,厚度小。

    通過上述(1)式的函數(shù)轉(zhuǎn)化,得到表2起伏度分級(jí)及CF值。由表2與圖4可知巖組因子中半堅(jiān)硬巖組滑坡易發(fā)性高,即巖組因子中半堅(jiān)硬巖組敏感性高。

    2.5 坡向

    有學(xué)者研究認(rèn)為,斜坡的朝向可能也與滑坡的易發(fā)性程度有關(guān),但要因地區(qū)而異。本文根據(jù)斜坡朝向?qū)⑵孪蛞蜃臃譃闁|、西、南、北與平地等5個(gè)二級(jí)因子。其中北朝向指角度在0°~45°與315°~360°范圍內(nèi),東朝向?yàn)槠孪蚪嵌仍?5°~135°范圍內(nèi),南朝向?yàn)槠孪蚪嵌仍?35°~225°范圍內(nèi),西朝向?yàn)槠孪蚪嵌仍?25°~315°范圍內(nèi)。各坡向分布情況如圖5所示。通過上述(1)式的函數(shù)轉(zhuǎn)化,得到表2中坡向分級(jí)及CF值。由表2與圖5可知坡向南區(qū)段滑坡易發(fā)性高,即坡向南區(qū)段敏感性高。

    圖3 鎮(zhèn)雄縣起伏度與滑坡分布圖

    圖4 鎮(zhèn)雄縣巖組與滑坡分布圖

    圖5 鎮(zhèn)雄縣坡向與滑坡分布圖

    2.6 與構(gòu)造線距離

    與構(gòu)造線的距離分區(qū)是在由Arc Map環(huán)境下通過對(duì)構(gòu)造線取緩沖區(qū)得到,如圖6所示。

    通過上述(1)式的函數(shù)轉(zhuǎn)化,得到表2與構(gòu)造線距離分級(jí)及CF值。由表2與圖6可知與構(gòu)造線距離在500~100 m區(qū)段滑坡易發(fā)性高,即與構(gòu)造線距離在500~100 m區(qū)段敏感性高。

    2.7 與水系距離

    河谷切割提供了發(fā)生滑坡的空間條件,是滑坡分布的主要孕育場(chǎng),鎮(zhèn)雄縣內(nèi)深切河谷的切割密度明顯地影響著地質(zhì)災(zāi)害分布的密度??傮w上,鎮(zhèn)雄縣大小河流共計(jì)21條,長(zhǎng)度超過10 k m的河流共17條,地表河流天然落差較大,水能資源較豐富,河流河床陡險(xiǎn)多石,所有河流均屬于長(zhǎng)江水系。河流兩側(cè)邊坡前緣常年受水浸潤(rùn),土體之間的摩擦力減小,土體的抗滑力下降,引發(fā)了大量滑坡。

    圖6 鎮(zhèn)雄縣構(gòu)造線緩沖區(qū)與滑坡分布圖

    通過上述(1)式的函數(shù)轉(zhuǎn)化,得到每個(gè)區(qū)段的CF值,見表2。由表2與圖7可知與水系距離在200~400 m區(qū)段滑坡易發(fā)性高,即與水系距離在200~400 m區(qū)段敏感性高。

    2.8 與道路距離

    在道路的修建和運(yùn)營(yíng)過程,對(duì)于斜坡土體有一個(gè)擾動(dòng)作用,可能破壞土體的平衡,尤其是在土體處于極限平衡狀態(tài)時(shí),故本文將公路基礎(chǔ)建設(shè)作為滑坡的一個(gè)誘發(fā)因子。

    圖7 鎮(zhèn)雄縣水系緩沖區(qū)與滑坡分布圖

    通過上述(1)式的函數(shù)轉(zhuǎn)化,得到每個(gè)級(jí)別的CF值,見表2。由表2與圖8可知與道路距離在<500 m區(qū)段滑坡易發(fā)性高,即與道路距離在<500 m區(qū)段敏感性高。

    3 回歸模型的建立

    3.1 主要致災(zāi)因子確定

    圖8 鎮(zhèn)雄縣道路緩沖區(qū)與滑坡分布圖

    計(jì)算出每一數(shù)據(jù)層的CF后,需要將不同數(shù)據(jù)層的CF進(jìn)行合并。假定要合并的兩數(shù)據(jù)層的CF值分別為x和y,合并得到的結(jié)果為Z,合并的公式如下:

    每個(gè)因子數(shù)據(jù)層的CF計(jì)算及合并均在Arcgis軟件中進(jìn)行。首先將每個(gè)因子數(shù)據(jù)層分別賦予公式(1)計(jì)算CF值,再在Arcgis軟件中進(jìn)行各個(gè)數(shù)據(jù)層的數(shù)據(jù)疊加及數(shù)據(jù)運(yùn)算。所涉及的功能模塊為處理柵格數(shù)據(jù)常用的柵格計(jì)算器,兩兩圖層疊加,每個(gè)柵格圖層的柵格單元中的CF值按照(2)式進(jìn)行運(yùn)算得到Z值,按一定的規(guī)則對(duì)合并后得到的Z值進(jìn)行重新劃分,得到滑坡形成條件分區(qū)結(jié)果。通過對(duì)比已有滑坡數(shù)據(jù),可確定每一種致災(zāi)因子對(duì)滑坡發(fā)生的影響程度,從而確定滑坡主要致災(zāi)因子。根據(jù)公式(2)進(jìn)行圖層疊加分析得出鎮(zhèn)雄縣滑坡主要致災(zāi)因子有:巖性(巖組)、坡度、坡向、起伏度、與水系距離、與道路距離。

    將鎮(zhèn)雄縣影響滑坡地質(zhì)災(zāi)害的主要致災(zāi)圖層:巖性(巖組)、坡度、坡向、起伏度、與水系距離、與道路距離,在Arcgis中進(jìn)行疊加分析,得到進(jìn)行統(tǒng)計(jì)計(jì)算的柵格單元,每個(gè)柵格單元均包含每一個(gè)主要致災(zāi)因子的單一條件分組(表3),從而可以認(rèn)定單元之間是相互獨(dú)立的,有利于統(tǒng)計(jì)分析的進(jìn)行。各個(gè)主要致災(zāi)因子疊加情況見圖9。

    3.2 回歸模型的建立

    統(tǒng)計(jì)學(xué)中回歸的含義是指,變量之間定量的依存關(guān)系[7-9]。

    對(duì)于標(biāo)準(zhǔn)的線性回歸方程[10-11]:

    在數(shù)理統(tǒng)計(jì)和數(shù)值分析上,常用最小二乘法對(duì)這樣的函數(shù)進(jìn)行擬合。也就是說已知各組自變量,要求確定各組自變量之間存在的關(guān)系,并且在在已知自變量的情況下,得出因變量Y值。

    將回歸原理用于滑坡形成條件分析中,xi表示各自變量,即各致災(zāi)因子指標(biāo)值,取值范圍[-1,1]。因變量Y就是滑坡形成條件有利級(jí)別,滑坡形成條件有利級(jí)別取值范圍[1,6]。

    隨機(jī)選取局部區(qū)域柵格單元(本文采用隨機(jī)選取200個(gè)柵格作為樣本數(shù)據(jù),樣本數(shù)據(jù)見附表A),采用各因子數(shù)據(jù)類CF值作為回歸模型的自變量,滑坡形成條件有利級(jí)別作為因變量,將數(shù)據(jù)導(dǎo)入Spss統(tǒng)計(jì)分析軟件,進(jìn)行多元線性回歸分析,得到鎮(zhèn)雄縣滑坡地質(zhì)災(zāi)害主要因子的回歸模型,如表4。將最終統(tǒng)計(jì)量帶入(3)式得出下式:

    式中,x1,x2,x3…x6,分別代表巖性(巖組)、坡度、坡向、與水系距離、與道路距離以及起伏度等致災(zāi)因子的CF值。

    表3 主要致災(zāi)因子疊加后部分單元格CF值

    圖9 主要致災(zāi)因子疊加形成條件分區(qū)圖

    表4 最終統(tǒng)計(jì)量值

    根據(jù)Spss得到的結(jié)果,顯著性水平Sig.均接近于零(如表4),說明自變量對(duì)因變量均有顯著影響。

    3.3 回歸模型的檢驗(yàn)

    根據(jù)所得回歸模型,進(jìn)行主要致災(zāi)因子圖層疊加分析,得到最終滑坡地質(zhì)災(zāi)害形成條件分區(qū),結(jié)果見表5。

    由表5可以看出,已發(fā)的滑坡主要集中在4、5、6類中,比例合計(jì)83.17%,所占滑坡總面積的83.17%。這一結(jié)果表明所建模型可以解釋80%左右的已發(fā)滑坡。因此,所建鎮(zhèn)雄縣滑坡地質(zhì)災(zāi)害主要致災(zāi)因子多元線性回歸模型精度較高。

    表5 模型所得分區(qū)結(jié)果分析表

    4 結(jié)論

    本次研究以GIS技術(shù)對(duì)鎮(zhèn)雄縣滑坡主要致災(zāi)因子進(jìn)行篩選分析,并采用CF法建立回歸模型對(duì)鎮(zhèn)雄縣滑坡形成條件進(jìn)行分區(qū)對(duì)比研究,得出了以下結(jié)論:

    (1)通過上述滑坡主要致災(zāi)因子的分析,可以得出,巖性(巖組)、坡度、坡向、起伏度、與水系距離、與道路距離等因子能作為影響鎮(zhèn)雄縣滑坡發(fā)育的主要因子。并結(jié)合致災(zāi)因子敏感性的分析,可得出:1)巖性(巖組):從巖組上看,半堅(jiān)硬巖組為滑坡易發(fā)巖組;2)地形:5°~15°是鎮(zhèn)雄縣滑坡地質(zhì)災(zāi)害易發(fā)坡度范圍,南方向?yàn)殒?zhèn)雄縣滑坡易發(fā)坡向范圍,起伏度在30~200 m之間,為鎮(zhèn)雄縣滑坡易發(fā)起伏度范圍;3)水系:水系緩沖區(qū)在200~400 m范圍為鎮(zhèn)雄縣滑坡易發(fā)區(qū)段;4)道路:道路緩沖區(qū)在<500 m范圍為鎮(zhèn)雄縣滑坡易發(fā)區(qū)段。

    (2)回歸模型研究表明,已發(fā)的滑坡主要集中在潛在不利區(qū)、較不利區(qū)、不利區(qū)3類中,比例合計(jì)83.17%,所占滑坡總面積的83.17%。這一結(jié)果表明所建模型可以解釋80%左右的已發(fā)滑坡。因此,所建鎮(zhèn)雄縣滑坡地質(zhì)災(zāi)害主要致災(zāi)因子多元線性回歸模型精度較高。

    [1]云南省地質(zhì)局第二區(qū)域地質(zhì)測(cè)量大隊(duì).1∶200 000中華人民共和國(guó)區(qū)域地質(zhì)調(diào)查報(bào)告(昭通幅)[R].1976-1978.

    [2]許沖,戴福初,姚鑫,等.基于GIS的汶川地震滑坡災(zāi)害影響因子[J].巖石力學(xué)與工程學(xué)報(bào),2010,29(5):2972-2981.

    [3]蘭恒星,伍法權(quán),王思敬.基于GIS的滑坡CF多元回歸模型及其應(yīng)用[J].山地學(xué)報(bào),2002,20(6):732-737.

    [4]蘭恒星,王苓涓,周成虎.云南小江流域滑坡關(guān)鍵影響因子研究[J].中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào),2003,14(1):100-106.

    [5]徐錫偉,聞學(xué)澤,葉建青.汶川M s8.0地震地表破裂帶及其發(fā)震構(gòu)造[J].地震地質(zhì),2008,04(3):10-15.

    [6]張培震,徐錫偉,聞學(xué)澤.2008年汶川8.0級(jí)地震發(fā)震斷裂的滑動(dòng)速率、復(fù)發(fā)周期和構(gòu)造成因[J].地球物理,2008,4(3):120-125.

    [7]何春雄.應(yīng)用數(shù)理統(tǒng)計(jì)基礎(chǔ)[M].廣州:華南理工大學(xué)出版社,2006.

    [8]周爽.社會(huì)統(tǒng)計(jì)分析:SPSS應(yīng)用教程[M].北京:清華大學(xué)出版社,2006.

    [9]龔有容.應(yīng)用統(tǒng)計(jì)學(xué)[M].上海:立信會(huì)計(jì)出版社,2006.

    [10]李雪平.基于GIS的區(qū)域斜坡穩(wěn)定性評(píng)價(jià)Logistic回歸模型研究[D].武漢:中國(guó)地質(zhì)大學(xué)工程學(xué)院,2005.

    [11]王征征.MAPGIS和SPSS的區(qū)域型滑坡危險(xiǎn)性區(qū)劃[D].北京:中國(guó)地質(zhì)大學(xué)地球科學(xué)與資源學(xué)院,2006.

    猜你喜歡
    鎮(zhèn)雄縣巖組坡向
    銀川市地下水賦存條件及動(dòng)態(tài)特征
    鎮(zhèn)雄縣坡頭鎮(zhèn)海塘歷史遺留硫磺冶煉廢渣綜合處置及生態(tài)恢復(fù)工程通過初步驗(yàn)收
    鎮(zhèn)雄縣總工會(huì):為高考考生助力
    鎮(zhèn)雄縣總工會(huì):開展送溫暖活動(dòng)
    鎮(zhèn)雄縣總工會(huì):開展禁毒宣傳活動(dòng)
    DEM地表坡向變率的向量幾何計(jì)算法
    淮南潘集深部勘查區(qū)15-2孔工程地質(zhì)巖組劃分
    白音華煤田三號(hào)露天礦區(qū)水文地質(zhì)條件分析
    青藏高原東緣高寒草甸坡向梯度上植物光合生理特征研究
    土壤與作物(2015年3期)2015-12-08 00:46:55
    坡向和坡位對(duì)小流域梯田土壤有機(jī)碳、氮變化的影響
    久久性视频一级片| 久热爱精品视频在线9| 一本久久精品| 不卡av一区二区三区| 热99国产精品久久久久久7| 在线观看免费视频网站a站| 18禁观看日本| 亚洲欧美激情在线| 久久青草综合色| 亚洲欧美日韩高清在线视频 | 欧美大码av| 麻豆国产av国片精品| 欧美日韩一级在线毛片| 永久免费av网站大全| 久久免费观看电影| 亚洲精品日韩在线中文字幕| 后天国语完整版免费观看| 久久性视频一级片| 人人妻人人澡人人看| 久久这里只有精品19| 欧美日韩视频精品一区| 99香蕉大伊视频| 在线观看免费午夜福利视频| 国产成人一区二区三区免费视频网站 | 捣出白浆h1v1| 19禁男女啪啪无遮挡网站| 夫妻性生交免费视频一级片| 精品免费久久久久久久清纯 | 免费观看人在逋| 精品国产乱码久久久久久男人| 美女脱内裤让男人舔精品视频| 精品免费久久久久久久清纯 | e午夜精品久久久久久久| 精品欧美一区二区三区在线| 男女高潮啪啪啪动态图| √禁漫天堂资源中文www| 欧美大码av| 精品少妇内射三级| 亚洲国产精品成人久久小说| 国产精品99久久99久久久不卡| 黄色 视频免费看| 精品亚洲成a人片在线观看| av天堂在线播放| 国产xxxxx性猛交| 欧美激情 高清一区二区三区| 九草在线视频观看| 久久久国产精品麻豆| 最近最新中文字幕大全免费视频 | 国产精品99久久99久久久不卡| 亚洲欧美中文字幕日韩二区| 亚洲精品在线美女| 国产亚洲精品第一综合不卡| 免费少妇av软件| 免费看不卡的av| 看免费成人av毛片| 成在线人永久免费视频| 欧美精品人与动牲交sv欧美| 又大又爽又粗| 午夜av观看不卡| 中文字幕人妻熟女乱码| 色婷婷av一区二区三区视频| 国产片内射在线| 久久国产亚洲av麻豆专区| 99久久99久久久精品蜜桃| 下体分泌物呈黄色| 老汉色∧v一级毛片| 亚洲国产精品国产精品| 国产午夜精品一二区理论片| 欧美黑人精品巨大| 19禁男女啪啪无遮挡网站| 大陆偷拍与自拍| 美女大奶头黄色视频| 在线观看免费午夜福利视频| 99香蕉大伊视频| 国产成人欧美在线观看 | 国产在视频线精品| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲情色 制服丝袜| 精品卡一卡二卡四卡免费| 97精品久久久久久久久久精品| 国产一级毛片在线| 狠狠精品人妻久久久久久综合| 无遮挡黄片免费观看| h视频一区二区三区| 国产欧美日韩精品亚洲av| 免费看不卡的av| 成在线人永久免费视频| 色精品久久人妻99蜜桃| 在线观看免费日韩欧美大片| 日本午夜av视频| 又大又爽又粗| 汤姆久久久久久久影院中文字幕| 久久天躁狠狠躁夜夜2o2o | 国产成人欧美| 欧美97在线视频| 日本91视频免费播放| 久热这里只有精品99| 成人免费观看视频高清| 日韩熟女老妇一区二区性免费视频| 高清黄色对白视频在线免费看| 免费在线观看黄色视频的| 69精品国产乱码久久久| 久久久久国产一级毛片高清牌| 高潮久久久久久久久久久不卡| 成年人午夜在线观看视频| 国产精品久久久久成人av| 亚洲精品国产区一区二| tube8黄色片| 午夜影院在线不卡| 一二三四社区在线视频社区8| 51午夜福利影视在线观看| 午夜av观看不卡| 最新在线观看一区二区三区 | 建设人人有责人人尽责人人享有的| 亚洲 国产 在线| 精品人妻在线不人妻| 欧美人与善性xxx| 人人妻人人澡人人爽人人夜夜| 久久免费观看电影| 99热网站在线观看| 人人澡人人妻人| 热re99久久国产66热| 国产黄色视频一区二区在线观看| 亚洲av成人不卡在线观看播放网 | 女性生殖器流出的白浆| 成年女人毛片免费观看观看9 | 午夜激情久久久久久久| 欧美精品av麻豆av| 午夜福利视频精品| 99热全是精品| 国产精品久久久av美女十八| 午夜福利视频在线观看免费| 免费在线观看完整版高清| 国产免费福利视频在线观看| 操出白浆在线播放| 老司机午夜十八禁免费视频| 亚洲精品国产色婷婷电影| 欧美人与性动交α欧美软件| 男人舔女人的私密视频| 欧美精品高潮呻吟av久久| 欧美日本中文国产一区发布| 国产精品 国内视频| 99久久精品国产亚洲精品| 日本av手机在线免费观看| 亚洲国产欧美网| 999精品在线视频| 国产精品久久久人人做人人爽| 久久人妻福利社区极品人妻图片 | 久久久久久久大尺度免费视频| 91成人精品电影| 五月开心婷婷网| 免费人妻精品一区二区三区视频| 色综合欧美亚洲国产小说| 久久热在线av| www日本在线高清视频| www.自偷自拍.com| 亚洲第一青青草原| 国产精品成人在线| 丰满迷人的少妇在线观看| 亚洲欧美一区二区三区国产| 日韩电影二区| 欧美精品一区二区大全| 日韩av免费高清视频| 日韩免费高清中文字幕av| 97在线人人人人妻| 欧美日韩亚洲综合一区二区三区_| 日本一区二区免费在线视频| xxxhd国产人妻xxx| 久久久久网色| av天堂在线播放| 啦啦啦 在线观看视频| 亚洲精品久久久久久婷婷小说| 中文字幕人妻熟女乱码| 人人澡人人妻人| 亚洲熟女毛片儿| 久热这里只有精品99| 天天操日日干夜夜撸| 99热网站在线观看| 国产高清国产精品国产三级| 一边摸一边抽搐一进一出视频| 国产日韩欧美亚洲二区| 女性生殖器流出的白浆| 一级片免费观看大全| 大码成人一级视频| 国语对白做爰xxxⅹ性视频网站| 热re99久久精品国产66热6| 国产又爽黄色视频| 午夜福利免费观看在线| 婷婷色综合大香蕉| 91老司机精品| 日韩一区二区三区影片| 丁香六月天网| 夫妻性生交免费视频一级片| 久久精品久久久久久噜噜老黄| 黄片播放在线免费| 国产视频一区二区在线看| 国产亚洲精品第一综合不卡| 19禁男女啪啪无遮挡网站| 亚洲色图 男人天堂 中文字幕| 国产精品一区二区在线不卡| 99久久综合免费| 亚洲精品国产区一区二| 亚洲国产精品999| 伊人亚洲综合成人网| 一边亲一边摸免费视频| 亚洲精品国产av成人精品| 妹子高潮喷水视频| 亚洲成人国产一区在线观看 | 久久精品国产综合久久久| 国产一级毛片在线| 免费少妇av软件| 天堂中文最新版在线下载| 两人在一起打扑克的视频| 亚洲成人手机| 国产亚洲一区二区精品| 久久久精品免费免费高清| 一级片免费观看大全| 欧美日韩亚洲国产一区二区在线观看 | 欧美人与性动交α欧美精品济南到| 亚洲国产精品一区三区| 少妇 在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品av久久久久免费| 日本wwww免费看| 亚洲色图综合在线观看| 免费在线观看影片大全网站 | 亚洲激情五月婷婷啪啪| 一本色道久久久久久精品综合| 天天躁日日躁夜夜躁夜夜| 新久久久久国产一级毛片| 亚洲国产日韩一区二区| 男人操女人黄网站| 激情五月婷婷亚洲| 一边摸一边做爽爽视频免费| 免费看不卡的av| 国产成人影院久久av| 国产精品免费大片| 久久久久网色| 男人添女人高潮全过程视频| 男女下面插进去视频免费观看| 色综合欧美亚洲国产小说| 亚洲成人手机| 操出白浆在线播放| 国产人伦9x9x在线观看| 水蜜桃什么品种好| 考比视频在线观看| 母亲3免费完整高清在线观看| xxx大片免费视频| 在现免费观看毛片| 97在线人人人人妻| 久久人妻福利社区极品人妻图片 | 自线自在国产av| 五月开心婷婷网| 后天国语完整版免费观看| 制服人妻中文乱码| 黄网站色视频无遮挡免费观看| 免费在线观看完整版高清| 国产精品一国产av| 午夜福利视频精品| 在线观看www视频免费| 亚洲精品美女久久久久99蜜臀 | 999精品在线视频| 大香蕉久久网| 成年人黄色毛片网站| 高潮久久久久久久久久久不卡| 美女视频免费永久观看网站| 成年人免费黄色播放视频| 亚洲伊人久久精品综合| 午夜福利视频精品| 成年av动漫网址| 国产精品av久久久久免费| 精品久久蜜臀av无| 手机成人av网站| 亚洲国产精品999| 久久精品久久久久久噜噜老黄| 亚洲,一卡二卡三卡| 99久久人妻综合| 丝袜美腿诱惑在线| 91麻豆av在线| 日韩 亚洲 欧美在线| 国产精品香港三级国产av潘金莲 | 国产片特级美女逼逼视频| 亚洲欧美色中文字幕在线| 亚洲av欧美aⅴ国产| 国产人伦9x9x在线观看| 啦啦啦啦在线视频资源| 热99久久久久精品小说推荐| 91国产中文字幕| 91精品国产国语对白视频| 精品人妻1区二区| 国产女主播在线喷水免费视频网站| 狂野欧美激情性xxxx| 手机成人av网站| 亚洲男人天堂网一区| 国产欧美日韩综合在线一区二区| 午夜91福利影院| 精品国产一区二区三区四区第35| 18在线观看网站| 国产精品 欧美亚洲| 国产亚洲一区二区精品| 女性生殖器流出的白浆| 日韩中文字幕视频在线看片| 免费日韩欧美在线观看| 亚洲视频免费观看视频| 99香蕉大伊视频| 最近手机中文字幕大全| 国产av一区二区精品久久| 日本av免费视频播放| 免费在线观看日本一区| 亚洲图色成人| 51午夜福利影视在线观看| 亚洲国产精品一区三区| 性色av一级| 色94色欧美一区二区| 狂野欧美激情性xxxx| 亚洲欧美日韩另类电影网站| 高清视频免费观看一区二区| 免费不卡黄色视频| 欧美人与性动交α欧美软件| 国产欧美日韩一区二区三 | 国产免费福利视频在线观看| 久久午夜综合久久蜜桃| 亚洲精品一区蜜桃| 天堂中文最新版在线下载| av天堂久久9| 久久久久久久国产电影| 51午夜福利影视在线观看| 久久ye,这里只有精品| www.自偷自拍.com| 国产av精品麻豆| svipshipincom国产片| 欧美日韩黄片免| 色网站视频免费| 欧美大码av| 黄色 视频免费看| 中文字幕另类日韩欧美亚洲嫩草| 男女高潮啪啪啪动态图| 久久久久久久大尺度免费视频| 国产成人精品久久久久久| 777米奇影视久久| 亚洲欧美日韩另类电影网站| 岛国毛片在线播放| 少妇的丰满在线观看| 亚洲成人手机| 国产精品国产三级国产专区5o| 深夜精品福利| 女人精品久久久久毛片| 国产精品久久久久成人av| 久久人人97超碰香蕉20202| 精品免费久久久久久久清纯 | 桃花免费在线播放| 国产日韩欧美视频二区| 国产免费现黄频在线看| 首页视频小说图片口味搜索 | 久久九九热精品免费| 久久精品国产综合久久久| 赤兔流量卡办理| 免费高清在线观看视频在线观看| 中文字幕高清在线视频| 美女中出高潮动态图| 最黄视频免费看| 男人爽女人下面视频在线观看| 亚洲精品美女久久av网站| 中文字幕色久视频| 又紧又爽又黄一区二区| 性高湖久久久久久久久免费观看| 国产成人啪精品午夜网站| 中文字幕亚洲精品专区| 亚洲欧美精品自产自拍| 欧美中文综合在线视频| 中文字幕制服av| 大香蕉久久网| 夫妻午夜视频| 丝瓜视频免费看黄片| 久久久久久久久久久久大奶| 国精品久久久久久国模美| av在线app专区| 欧美日韩成人在线一区二区| 精品卡一卡二卡四卡免费| 久久九九热精品免费| 国产一区二区 视频在线| 新久久久久国产一级毛片| 亚洲人成电影观看| 亚洲欧美清纯卡通| 久久久久精品人妻al黑| 亚洲国产中文字幕在线视频| 超碰97精品在线观看| 叶爱在线成人免费视频播放| www.熟女人妻精品国产| 亚洲少妇的诱惑av| 亚洲五月婷婷丁香| 在现免费观看毛片| 高清视频免费观看一区二区| 中文字幕色久视频| 欧美 日韩 精品 国产| cao死你这个sao货| 美女福利国产在线| 精品一区二区三卡| 日本av免费视频播放| 亚洲一区中文字幕在线| 成人免费观看视频高清| www.熟女人妻精品国产| 欧美av亚洲av综合av国产av| 最新的欧美精品一区二区| 女人爽到高潮嗷嗷叫在线视频| 深夜精品福利| 熟女少妇亚洲综合色aaa.| 亚洲国产成人一精品久久久| 欧美日韩综合久久久久久| 国产男女内射视频| 亚洲专区中文字幕在线| 亚洲国产欧美一区二区综合| 一级毛片我不卡| 91麻豆精品激情在线观看国产 | 久久人妻熟女aⅴ| 久久久国产一区二区| 久久久久久久国产电影| 最新的欧美精品一区二区| 亚洲自偷自拍图片 自拍| 精品人妻熟女毛片av久久网站| 啦啦啦在线观看免费高清www| 黑人巨大精品欧美一区二区蜜桃| 中文字幕人妻丝袜制服| 高清不卡的av网站| a级片在线免费高清观看视频| 久久久久视频综合| 久久久精品区二区三区| 交换朋友夫妻互换小说| 91精品国产国语对白视频| 99re6热这里在线精品视频| 婷婷丁香在线五月| 亚洲专区中文字幕在线| 国产精品 国内视频| 一本大道久久a久久精品| 日韩,欧美,国产一区二区三区| 日本猛色少妇xxxxx猛交久久| av线在线观看网站| 国产色视频综合| 老司机在亚洲福利影院| 欧美日本中文国产一区发布| 国产人伦9x9x在线观看| 飞空精品影院首页| 亚洲精品久久久久久婷婷小说| 久久久久网色| 大陆偷拍与自拍| 久久国产亚洲av麻豆专区| 高清视频免费观看一区二区| 久久久久视频综合| 午夜视频精品福利| 免费高清在线观看视频在线观看| 亚洲精品日本国产第一区| 国产高清不卡午夜福利| av在线播放精品| 一本一本久久a久久精品综合妖精| 爱豆传媒免费全集在线观看| 成人国产一区最新在线观看 | 老汉色∧v一级毛片| 亚洲精品日韩在线中文字幕| 亚洲精品日本国产第一区| 国产熟女欧美一区二区| xxxhd国产人妻xxx| 欧美变态另类bdsm刘玥| 各种免费的搞黄视频| 日韩制服丝袜自拍偷拍| 男人操女人黄网站| 伦理电影免费视频| 一区二区三区乱码不卡18| 看免费成人av毛片| 中国国产av一级| 日韩精品免费视频一区二区三区| 欧美黄色淫秽网站| 成人国产一区最新在线观看 | 一区二区三区激情视频| 交换朋友夫妻互换小说| 亚洲自偷自拍图片 自拍| 国产精品一区二区在线观看99| 91老司机精品| 国产一卡二卡三卡精品| www.熟女人妻精品国产| 在线观看人妻少妇| 亚洲国产av新网站| 午夜精品国产一区二区电影| 成年人黄色毛片网站| 国产亚洲一区二区精品| 日韩av在线免费看完整版不卡| 曰老女人黄片| 国产在线视频一区二区| 久久午夜综合久久蜜桃| 久久久精品免费免费高清| 久久久国产一区二区| 久久国产精品男人的天堂亚洲| 欧美黄色淫秽网站| 精品亚洲成国产av| 免费人妻精品一区二区三区视频| 97人妻天天添夜夜摸| 国产成人精品无人区| 又黄又粗又硬又大视频| 久久中文字幕一级| 亚洲精品国产av成人精品| 叶爱在线成人免费视频播放| 人人澡人人妻人| kizo精华| 久久国产精品人妻蜜桃| 少妇被粗大的猛进出69影院| 欧美中文综合在线视频| avwww免费| 免费久久久久久久精品成人欧美视频| 亚洲专区中文字幕在线| 99香蕉大伊视频| 中文字幕制服av| 精品熟女少妇八av免费久了| 老汉色∧v一级毛片| 久9热在线精品视频| 欧美日韩亚洲高清精品| 两性夫妻黄色片| 建设人人有责人人尽责人人享有的| 十八禁人妻一区二区| 一区二区三区激情视频| 久久精品亚洲av国产电影网| 一区二区三区激情视频| 大香蕉久久成人网| 国产一级毛片在线| 欧美亚洲 丝袜 人妻 在线| 侵犯人妻中文字幕一二三四区| 成年美女黄网站色视频大全免费| 一本大道久久a久久精品| 一本综合久久免费| 色94色欧美一区二区| 丝袜美足系列| 午夜福利免费观看在线| 久久精品久久精品一区二区三区| 国产不卡av网站在线观看| 中国国产av一级| 黄片小视频在线播放| 欧美大码av| 91精品三级在线观看| 亚洲精品久久午夜乱码| 久热这里只有精品99| www.熟女人妻精品国产| 亚洲精品一二三| 久久影院123| 男女床上黄色一级片免费看| 黄色a级毛片大全视频| 国产一区二区三区综合在线观看| 国产人伦9x9x在线观看| 成年动漫av网址| 精品一区在线观看国产| 免费在线观看完整版高清| 久久青草综合色| 视频区欧美日本亚洲| netflix在线观看网站| 99国产综合亚洲精品| 丰满少妇做爰视频| 国产男女超爽视频在线观看| 一本综合久久免费| 亚洲一码二码三码区别大吗| 久久国产精品影院| 欧美精品一区二区免费开放| 蜜桃国产av成人99| 视频区图区小说| 亚洲成人免费av在线播放| 婷婷色综合大香蕉| 丝瓜视频免费看黄片| 国产97色在线日韩免费| 日韩 欧美 亚洲 中文字幕| svipshipincom国产片| 欧美黄色片欧美黄色片| 啦啦啦在线免费观看视频4| 色网站视频免费| 一本久久精品| 大型av网站在线播放| 秋霞在线观看毛片| 波多野结衣一区麻豆| 韩国精品一区二区三区| 一级毛片 在线播放| 在线观看免费视频网站a站| 亚洲情色 制服丝袜| 国产欧美日韩精品亚洲av| 国产一区二区 视频在线| 午夜免费男女啪啪视频观看| 七月丁香在线播放| 三上悠亚av全集在线观看| 性少妇av在线| 亚洲成人免费av在线播放| 国产亚洲午夜精品一区二区久久| a级片在线免费高清观看视频| 国产成人免费无遮挡视频| 欧美xxⅹ黑人| 成年人午夜在线观看视频| 汤姆久久久久久久影院中文字幕| 亚洲精品中文字幕在线视频| 亚洲国产av影院在线观看| 在线观看一区二区三区激情| 一级毛片电影观看| 亚洲国产精品999| 少妇猛男粗大的猛烈进出视频| 女性生殖器流出的白浆| 最近手机中文字幕大全| 久久久亚洲精品成人影院| 丰满饥渴人妻一区二区三| 99香蕉大伊视频| 国产亚洲精品第一综合不卡| svipshipincom国产片| 久久鲁丝午夜福利片| 丁香六月欧美| 日本a在线网址| 国产亚洲av高清不卡| 久久精品久久久久久噜噜老黄| av在线老鸭窝| 国产高清国产精品国产三级| 啦啦啦 在线观看视频| 国产精品久久久久成人av| 超色免费av| 欧美日韩亚洲高清精品| 一本—道久久a久久精品蜜桃钙片|