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

    斜坡地質(zhì)災(zāi)害敏感性評(píng)價(jià)中地勢起伏度提取最佳尺度研究

    2020-11-19 02:28:46張明媚薛永安
    關(guān)鍵詞:變點(diǎn)斜坡均值

    張明媚,薛永安

    (1.太原理工大學(xué) 礦業(yè)工程學(xué)院,太原 030024;2.山西能源學(xué)院 地質(zhì)測繪工程系,山西 晉中 030600)

    斜坡地質(zhì)災(zāi)害敏感性是指在特殊地形或某些因素作用下發(fā)生滑坡、崩塌災(zāi)害的可能性[1],地質(zhì)災(zāi)害敏感性評(píng)價(jià)的本質(zhì)就是利用數(shù)學(xué)語言來表達(dá)在給定的地質(zhì)環(huán)境條件下地質(zhì)災(zāi)害在空間上的發(fā)生概率[2-3]。目前,對(duì)于地質(zhì)災(zāi)害敏感性評(píng)價(jià)研究已有許多成果。ACHOUR et al[4]基于層次分析法和信息量模型開展了滑坡敏感性研究。YALCIN et al[5]的對(duì)比研究表明層次分析法對(duì)滑坡敏感性的分區(qū)與實(shí)際情況更加符合。NICU et al[6]采用坡度、高程、曲率、巖性、降水、土地利用、地形潤濕性指數(shù)、地形、坡向、平面曲率和距離河流距離等11個(gè)條件因子對(duì)滑坡的敏感性程度進(jìn)行了統(tǒng)計(jì)分析。HONG et al[7]選取坡度、坡向、海拔、平面曲率、剖面曲率、巖性、距離斷層距離、距離河流距離、距離道路距離、土地利用、歸一化植被指數(shù)等15個(gè)因子作為滑坡災(zāi)害敏感性評(píng)價(jià)因子集合。THAI et al[8]選用坡度、曲率、海拔等15個(gè)影響因子作為滑坡災(zāi)害敏感性評(píng)價(jià)因子集合。而KALANTAR et al[9]對(duì)Dodangeh流域進(jìn)行滑坡災(zāi)害敏感性評(píng)價(jià)研究中采用的因子為曲率、平面曲率、剖面曲率、高程、坡度、坡向、距離斷層距離、距離河流距離、地形濕度指數(shù)、河流水力指數(shù)、地形粗糙度指數(shù)、沉積物搬運(yùn)指數(shù)、巖性及土地利用。在上述這些研究中,以高程、坡度、坡向、曲率等為代表的數(shù)字地形因子是斜坡地質(zhì)災(zāi)害敏感性評(píng)價(jià)的重要因子。

    然而,作為描述一個(gè)區(qū)域地形特征的宏觀性指標(biāo),地勢起伏度指在一個(gè)特定的區(qū)域內(nèi),最高點(diǎn)海拔高度與最低點(diǎn)海拔高度的差值,它能反映地表起伏變化,與斜坡地質(zhì)災(zāi)害的發(fā)生存在很大的相關(guān)性[10]。按照地貌發(fā)育的基本理論,一種地貌類型存在一個(gè)使最大高差達(dá)到相對(duì)穩(wěn)定的最佳分析區(qū)域,傳統(tǒng)方法是利用柵格窗口的遞增方法來尋找最佳的分析窗口[11]。這是地貌分析中地勢起伏度計(jì)算的最佳統(tǒng)計(jì)單元思路,即隨著統(tǒng)計(jì)單元半徑的增大,地勢起伏度隨之變化,到達(dá)某一臨界點(diǎn)后趨于穩(wěn)定,臨界點(diǎn)即為地勢起伏度的最佳統(tǒng)計(jì)單元,也是真實(shí)反映地形起伏的DEM地勢起伏度提取單位面積。針對(duì)地勢起伏度最佳提取單元的研究[12-13],眾多學(xué)者以不同尺度的DEM在不同實(shí)驗(yàn)區(qū)開展了研究工作,結(jié)果各異。涂漢明等[14]通過對(duì)全國600個(gè)樣點(diǎn)和2個(gè)小區(qū)的詳細(xì)研究,運(yùn)用模糊數(shù)學(xué)方法,確定了全國地勢起伏度最佳統(tǒng)計(jì)單元為21 km2.張軍等[15]則得出基于1∶250 000 DEM的新疆地勢起伏度計(jì)算的最佳統(tǒng)計(jì)單元為2.56 km2.趙斌濱等[16]基于90 m分辨率SRTM-DEM數(shù)據(jù)提取得到中國地勢起伏度的最佳統(tǒng)計(jì)單元面積為2.25 km2.王讓虎等[17]基于ASTER GDEM數(shù)據(jù)提取得到中國東北地區(qū)地形起伏度的最佳統(tǒng)計(jì)單元面積為2.62 km2.陳學(xué)兄等[18]以30 m分辨率ASTER GDEM數(shù)據(jù)為基礎(chǔ),運(yùn)用均值變點(diǎn)分析法計(jì)算得到0.260 1 km2為山西地形起伏度提取的最佳統(tǒng)計(jì)單元。楊艷林等[19]基于ASTER GDEM數(shù)據(jù)分析得出長江中游的地形起伏度最佳分析窗口為19×19像元,面積約0.324 9 km2.而郭芳芳等[20]從地勢起伏度分析在區(qū)域斜坡災(zāi)害評(píng)價(jià)中的應(yīng)用出發(fā),尋找反映斜坡點(diǎn)對(duì)應(yīng)的真實(shí)地勢起伏度單元,對(duì)不同統(tǒng)計(jì)單元地形起伏度值與斜坡個(gè)數(shù)統(tǒng)計(jì)曲線峰值分布特征進(jìn)行分析,確定了研究區(qū)地形起伏度最佳統(tǒng)計(jì)單元面積為4 km2.

    基于上述研究,本文以山西省太原市西山地質(zhì)塊體為研究區(qū),以ASTER GDEM V2為地勢起伏度提取數(shù)據(jù)源,以2012年斜坡地質(zhì)災(zāi)害分布信息為斜坡災(zāi)害發(fā)育敏感性評(píng)價(jià)數(shù)據(jù)源,以均值變點(diǎn)法開展斜坡災(zāi)害敏感性評(píng)價(jià)中地勢起伏度因子最佳提取單元研究。

    1 數(shù)據(jù)與方法

    1.1 研究區(qū)與數(shù)據(jù)源

    研究區(qū)位于山西省中部、太原盆地的北端,北接太原市陽曲縣,西鄰太原市古交市、清徐縣,東鄰太原盆地,南靠太原市清徐縣,包括太原市尖草坪區(qū)、萬柏林區(qū)和晉源區(qū),總面積441.063 km2,見圖1.研究區(qū)的地理坐標(biāo)為東經(jīng)112°15′03″-112°29′25″,北緯37°38′54″-38°01′14″.

    圖1 研究區(qū)地理位置圖Fig.1 Geographical map of the studied area

    1.1.1DEM數(shù)據(jù)

    ASTER GDEM(先進(jìn)星載熱輻射和反射儀全球數(shù)字高程模型)是根據(jù)NASA的新一代對(duì)地觀測衛(wèi)星Terra的詳盡觀測結(jié)果制作完成的,其全球空間分辨率約為30 m,垂直分辨率為20 m,空間參考為WGS84/EGM96[21].ASTER GDEM目前共有兩版數(shù)據(jù),本文研究采用第2版,于2011年10月發(fā)布,數(shù)據(jù)來源于中國科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心地理空間數(shù)據(jù)云平臺(tái)。以研究區(qū)ASTER GDEM V2重采樣后30 m分辨率DEM為基準(zhǔn)數(shù)據(jù)(數(shù)據(jù)參數(shù)及精度見表1).

    表1 研究區(qū)主要地形參數(shù)及信息源精度Table 1 Main terrain parameters and information source accuracy of the studied area

    1.1.2斜坡地質(zhì)災(zāi)害數(shù)據(jù)

    從太原市規(guī)劃和自然資源局收集到研究區(qū)2012年斜坡地質(zhì)災(zāi)害分布信息,見圖2.

    圖2 研究區(qū)斜坡地質(zhì)災(zāi)害空間分布圖Fig.2 Spatial distribution of slope geological hazards in the studied area

    1.2 均值變點(diǎn)法

    目前,最佳地勢起伏度統(tǒng)計(jì)單元的尋找主要依賴人工判斷,以平均地勢起伏度隨統(tǒng)計(jì)單元面積變化的轉(zhuǎn)折點(diǎn)作為最佳統(tǒng)計(jì)單元面積,也就是人工判斷曲線上由陡變緩的位置;這種判斷方法顯然主觀性太強(qiáng),不夠準(zhǔn)確。隨著最佳地勢起伏度統(tǒng)計(jì)單元研究工作的深入[22],統(tǒng)計(jì)學(xué)領(lǐng)域的均值變點(diǎn)法逐漸被用作計(jì)算曲線上由陡變緩點(diǎn)位的客觀方法。

    均值變點(diǎn)法的原理及步驟如下[23]:

    1) 檢驗(yàn)是否存在變點(diǎn),即檢驗(yàn)原假設(shè)H0是否存在變點(diǎn)。若H0被接受,則該數(shù)據(jù)序列沒有變點(diǎn);若H0被否定,則假設(shè)該序列中至多存在q個(gè)變點(diǎn),對(duì)m1,m2,…,mq變點(diǎn)進(jìn)行估計(jì)。

    2) 估計(jì)變點(diǎn)個(gè)數(shù)。

    3) 估計(jì)變點(diǎn)處均值的跳躍度。

    在實(shí)際估計(jì)中,如果先驗(yàn)知識(shí)足以肯定序列中變點(diǎn)的存在,則可跳過步驟1)而直接進(jìn)入后幾步。

    均值變點(diǎn)問題的離散數(shù)據(jù)模型如下。

    設(shè)xt=at+et,t=1,2,…,N,則

    a1=…=am1-1=b1,am1=…=am2-1=b2,…,amq=…=aN=bq+1.

    (1)

    式中:1

    H0∶b1=b2=…=bq+1.

    (2)

    此處特別強(qiáng)調(diào)本檢驗(yàn)與多樣本檢驗(yàn)的差別,即此檢驗(yàn)中m1,m2,…,mq為未知。

    原假設(shè)H0的檢驗(yàn)步驟如下。

    1) 令i=2,…,N,每個(gè)i將樣本分為兩段:

    x1,x2,…,xi-1和xi,xi+1,…,xN.

    (3)

    (4)

    2) 計(jì)算統(tǒng)計(jì)量:

    (5)

    3) 計(jì)算期望值:

    E(S-Si),i=2,3,…,N.

    (6)

    4) 求極大值:

    (7)

    在平均意義下認(rèn)為:

    S*=min(S2,…,SN) .

    (8)

    5) 取檢驗(yàn)顯著性水平為α,計(jì)算C值:

    C=σ2(2lg lgN+lg lg lgN-lg π-2lg(-0.5lg(1-α))) .

    (9)

    如果σ2未知,則用下面估計(jì)來替代:

    (10)

    6) 如果S-S*>C,則否定H0,認(rèn)為無變點(diǎn),否則接受H0,該檢驗(yàn)有漸近水平α.

    平均來說,如果序列中存在變點(diǎn),則S和Si的差距會(huì)因變點(diǎn)的存在而增大,所以均值變點(diǎn)法原理明確,操作簡單。同時(shí),均值變點(diǎn)法對(duì)只有一個(gè)變點(diǎn)的檢驗(yàn)最為有效,存在多個(gè)變點(diǎn)時(shí)有可能因?yàn)榫档亩啻紊刀窒薙和Si之間的差距。

    2 結(jié)果與討論

    2.1 平均地勢起伏度與統(tǒng)計(jì)單元之間的關(guān)系

    基于ArcGIS平臺(tái),利用鄰域分析法提取研究區(qū)不同面積單元的地勢起伏度,提取單元分別為30 m分辨率DEM數(shù)據(jù)網(wǎng)格的2×2,3×3,4×4,5×5,…,25×25,對(duì)應(yīng)的提取單元面積(km2)分別為:0.003 6,0.008 1,0.014 4,0.022 5,0.032 4,0.044 1,0.057 6,0.072 9,0.090 0,0.108 9,0.129 6,0.152 1,0.176 4,0.202 5,0.230 4,0.260 1,0.291 6,0.324 9,0.360 0,0.396 9,0.435 6,0.476 1,0.518 4,0.562 5,共24個(gè)面積單元。經(jīng)統(tǒng)計(jì),研究區(qū)不同提取單元大小與平均地勢起伏度之間關(guān)系如表2所示。

    由表2可以看出,隨著提取的地勢起伏度統(tǒng)計(jì)單元的增大,研究區(qū)平均地勢起伏度呈現(xiàn)上升趨勢,從最小的13.53 m上升到173.96 m,這有助于分析研究區(qū)斜坡地質(zhì)災(zāi)害發(fā)育與地形起伏之間的關(guān)系。至于哪個(gè)尺度的分析窗口是最佳窗口,可以通過區(qū)域地勢起伏度與斜坡地質(zhì)災(zāi)害發(fā)育之間的關(guān)系進(jìn)行客觀的分析評(píng)價(jià)。

    表2 研究區(qū)不同統(tǒng)計(jì)單元與平均地勢起伏度之間的關(guān)系Table 2 Relationship between statistical unit and average relief amplitude of the area

    以統(tǒng)計(jì)分析單元面積為橫坐標(biāo),相應(yīng)的平均地勢起伏度計(jì)算值為縱坐標(biāo),對(duì)地勢起伏度隨統(tǒng)計(jì)單元面積的變化情況進(jìn)行擬合,結(jié)果如圖3所示。

    對(duì)于不同統(tǒng)計(jì)單元面積與所提取地勢起伏度之間的關(guān)系,通過對(duì)數(shù)方程進(jìn)行擬合,得到擬合曲線方程。通過統(tǒng)計(jì)學(xué)檢驗(yàn),擬合效果良好,擬合方程為:

    y=34.097lnx-52.50,R2=0.972 9 .

    (11)

    式中:y為平均地勢起伏度,m;x為統(tǒng)計(jì)單元面積,km2;R2為決定系數(shù)。

    從圖3的曲線變化趨勢來看,應(yīng)有1個(gè)變點(diǎn)存在,過了變點(diǎn)后地勢起伏度的增速變緩,而這個(gè)變點(diǎn)所對(duì)應(yīng)的面積就是研究區(qū)提取地勢起伏度的最佳統(tǒng)計(jì)單元面積。

    圖3 地勢起伏度與統(tǒng)計(jì)單元面積對(duì)應(yīng)關(guān)系擬合曲線Fig.3 Fitting curve of the corresponding relation between relief amplitude and area of statistical unit

    2.2 地勢起伏度提取最佳統(tǒng)計(jì)單元分析

    對(duì)不同統(tǒng)計(jì)單元的平均地勢起伏度、不同統(tǒng)計(jì)單元地勢起伏度分級(jí)中的斜坡地質(zhì)災(zāi)害峰值、地勢起伏度峰值進(jìn)行統(tǒng)計(jì),結(jié)果見圖4.

    圖4 平均地勢起伏度、地勢起伏度峰值、斜坡地質(zhì)災(zāi)害峰值與統(tǒng)計(jì)單元關(guān)系圖Fig.4 Relationship between average value of relief amplitude, peak value of relief amplitude, and peak value of slope geological hazards as one side and statistical unit as another side

    根據(jù)圖4統(tǒng)計(jì)結(jié)果分別構(gòu)建樣本序列,利用均值變點(diǎn)法的公式編制Excel程序,計(jì)算樣本序列的統(tǒng)計(jì)量S和Si,求得S-Si(見表3),構(gòu)建S與Si差值的變化擬合曲線。

    表3 均值變點(diǎn)法分析統(tǒng)計(jì)結(jié)果Table 3 Statistical results of mean change point

    平均地勢起伏度樣本序列S值為16.994 5,S與Si差值的變化曲線見圖5.地勢起伏度峰值樣本序列S的值為20.130 4,S與Si差值的變化曲線見圖6.斜坡地質(zhì)災(zāi)害峰值樣本序列S的值為89.359 3,S與Si差值的變化曲線見圖7.

    圖5 基于平均地勢起伏度的S-Si變化曲線Fig.5 Variation curve based on average value of relief amplitude S-Si

    圖6 基于地勢起伏度峰值的S-Si變化曲線Fig.6 Variation curve based on peak value of relief amplitude S-Si

    圖7 基于斜坡地質(zhì)災(zāi)害峰值的S-Si變化曲線Fig.7 Variation curve based on peak value of slope geological hazards S-Si

    2.3 討論

    1) 圖4顯示,每一種統(tǒng)計(jì)單元都存在一個(gè)斜坡地質(zhì)災(zāi)害發(fā)育的峰值區(qū)間,隨統(tǒng)計(jì)單元的增大,峰值逐漸減小,到達(dá)某一個(gè)臨界值后峰值衰減趨于穩(wěn)定。

    2) 圖5顯示,i為11時(shí),S與Si差值達(dá)最大,i=11處即為變點(diǎn),對(duì)應(yīng)統(tǒng)計(jì)單元為12×12網(wǎng)格。因此,由平均地勢起伏度以均值變點(diǎn)法來分析得到研究區(qū)的最佳地勢起伏度提取單元為12×12網(wǎng)格,最佳統(tǒng)計(jì)面積為0.129 6 km2.

    3) 圖6顯示,i為11時(shí),S與Si差值達(dá)最大,i=11處即為變點(diǎn),對(duì)應(yīng)統(tǒng)計(jì)單元為12×12網(wǎng)格。因此,由地勢起伏度峰值以均值變點(diǎn)法來分析得到研究區(qū)的最佳地勢起伏度提取單元為12×12網(wǎng)格,最佳統(tǒng)計(jì)面積為0.129 6 km2.

    4) 圖7顯示,i為8時(shí),S與Si差值達(dá)最大,i=8處即為變點(diǎn),對(duì)應(yīng)統(tǒng)計(jì)單元為9×9網(wǎng)格。因此,通過斜坡地質(zhì)災(zāi)害峰值以均值變點(diǎn)法來分析得到研究區(qū)的最佳地勢起伏度提取單元為9×9網(wǎng)格,最佳統(tǒng)計(jì)面積為0.072 9 km2.

    5) 對(duì)比三種樣本序列所計(jì)算得到的最佳地勢起伏度提取單元,使用平均地勢起伏度和地勢起伏度峰值樣本序列計(jì)算結(jié)果一致,均為12×12網(wǎng)格,而斜坡地質(zhì)災(zāi)害峰值樣本序列計(jì)算結(jié)果則為9×9網(wǎng)格。作為客觀存在的區(qū)域地貌特征因子,斜坡地質(zhì)災(zāi)害演變頻繁,而平均地勢起伏度演變較小。因此,以研究區(qū)數(shù)字地貌角度分析,基于平均地勢起伏度所計(jì)算的最佳提取單元更合理;而以地質(zhì)災(zāi)害發(fā)育角度分析,基于斜坡地質(zhì)災(zāi)害峰值所計(jì)算的最佳提取單元?jiǎng)t更合理。

    6) 針對(duì)斜坡地質(zhì)災(zāi)害敏感性評(píng)價(jià)中地勢起伏度最佳統(tǒng)計(jì)單元開展研究,重點(diǎn)關(guān)注的是研究區(qū)地勢起伏度與斜坡地質(zhì)災(zāi)害發(fā)育之間的關(guān)系,因此所尋找的是能反映地勢起伏度與斜坡地質(zhì)災(zāi)害發(fā)育之間關(guān)系的最佳統(tǒng)計(jì)單元,而不是數(shù)字地貌研究中區(qū)域地勢起伏度計(jì)算的最佳統(tǒng)計(jì)單元。因此,選擇9×9網(wǎng)格作為研究區(qū)斜坡地質(zhì)災(zāi)害敏感性評(píng)價(jià)中地勢起伏度因子最佳提取單元。

    3 結(jié)論

    1) 地勢起伏度與統(tǒng)計(jì)單元面積之間存在顯著的對(duì)數(shù)變化關(guān)系,研究區(qū)內(nèi)地勢起伏度隨統(tǒng)計(jì)單元面積的增大而迅速增大,當(dāng)統(tǒng)計(jì)單元面積達(dá)到一定數(shù)值后,增大速度逐漸變緩,最后近似趨于平穩(wěn),這為開展地勢起伏度最佳提取單元提供了基礎(chǔ)。

    2) 地勢起伏度提取最佳統(tǒng)計(jì)單元與地質(zhì)災(zāi)害分布峰值、平均地勢起伏度和地勢起伏度峰值之間存在明顯的相關(guān)性,目前的研究主要基于平均地勢起伏度,利用均值變點(diǎn)法或主觀方式分析確定最佳提取單元,這給斜坡地質(zhì)災(zāi)害敏感性評(píng)價(jià)帶來一定的不確定性,所以應(yīng)綜合三種情況開展統(tǒng)計(jì)并綜合分析后確定研究區(qū)地勢起伏度最佳統(tǒng)計(jì)尺度。

    3) 以研究區(qū)ASTER GDEM重采樣30 m分辨率DEM和2012年斜坡類地質(zhì)災(zāi)害空間分布信息為數(shù)據(jù)源,以2×2,3×3,…,25×25共24個(gè)分析窗口提取了研究區(qū)的地勢起伏度,以均值變點(diǎn)法對(duì)不同統(tǒng)計(jì)單元的平均地勢起伏度、不同統(tǒng)計(jì)單元地勢起伏度分級(jí)中的斜坡地質(zhì)災(zāi)害峰值、地勢起伏度峰值進(jìn)行統(tǒng)計(jì)分析。綜合三種樣本序列計(jì)算結(jié)果,本文選擇9×9網(wǎng)格作為研究區(qū)斜坡類地質(zhì)災(zāi)害敏感性評(píng)價(jià)時(shí)的最佳地勢起伏度提取單元,這為研究區(qū)斜坡類地質(zhì)災(zāi)害敏感性評(píng)價(jià)研究提供了可靠的地勢起伏度因子,同時(shí)也為同類研究提供了數(shù)字地形因子較為準(zhǔn)確的提取方法。

    猜你喜歡
    變點(diǎn)斜坡均值
    回歸模型參數(shù)的變點(diǎn)檢測方法研究
    正態(tài)分布序列均值變點(diǎn)檢測的貝葉斯方法
    基于二元分割的多變點(diǎn)估計(jì)
    獨(dú)立二項(xiàng)分布序列變點(diǎn)的識(shí)別方法
    信仰的“斜坡”
    夢是長長的斜坡(外一首)
    天津詩人(2017年2期)2017-11-29 01:24:12
    均值不等式失效時(shí)的解決方法
    均值與方差在生活中的應(yīng)用
    關(guān)于均值有界變差函數(shù)的重要不等式
    對(duì)偶均值積分的Marcus-Lopes不等式
    亚洲 国产 在线| 亚洲av中文字字幕乱码综合| 久久久久精品国产欧美久久久| 亚洲人成网站在线播| 人妻丰满熟妇av一区二区三区| 小蜜桃在线观看免费完整版高清| 久久久久久久亚洲中文字幕 | 搞女人的毛片| 搞女人的毛片| 日韩中文字幕欧美一区二区| 久久欧美精品欧美久久欧美| 国产午夜精品久久久久久一区二区三区 | 亚洲男人的天堂狠狠| 天堂影院成人在线观看| 黄色丝袜av网址大全| 在线国产一区二区在线| 日韩精品中文字幕看吧| 又黄又粗又硬又大视频| 中文资源天堂在线| 人妻丰满熟妇av一区二区三区| 欧美大码av| 一卡2卡三卡四卡精品乱码亚洲| 欧美大码av| 国产亚洲精品一区二区www| 亚洲人成网站在线播放欧美日韩| 亚洲一区二区三区色噜噜| 精品人妻偷拍中文字幕| av国产免费在线观看| 色吧在线观看| 乱人视频在线观看| 亚洲欧美一区二区三区黑人| 久久亚洲真实| 麻豆国产av国片精品| 最新在线观看一区二区三区| 欧美另类亚洲清纯唯美| 国产欧美日韩精品一区二区| 韩国av一区二区三区四区| 手机成人av网站| 神马国产精品三级电影在线观看| 男女下面进入的视频免费午夜| 久久久成人免费电影| 岛国视频午夜一区免费看| 免费看美女性在线毛片视频| 成人无遮挡网站| x7x7x7水蜜桃| 中文在线观看免费www的网站| 一区二区三区国产精品乱码| 色老头精品视频在线观看| 身体一侧抽搐| 久久人人精品亚洲av| 18禁国产床啪视频网站| 香蕉丝袜av| 69av精品久久久久久| 精品一区二区三区视频在线 | 久久久久国产精品人妻aⅴ院| 国产成人a区在线观看| 少妇的逼水好多| av天堂中文字幕网| 老熟妇仑乱视频hdxx| 精品国产三级普通话版| 国产精品爽爽va在线观看网站| 嫩草影院入口| bbb黄色大片| 国产v大片淫在线免费观看| 精品久久久久久久久久免费视频| 欧美乱色亚洲激情| 国产精品亚洲一级av第二区| 欧美日本亚洲视频在线播放| 国产高清激情床上av| 最近视频中文字幕2019在线8| 高清毛片免费观看视频网站| 欧美色欧美亚洲另类二区| 国产高清videossex| 在线国产一区二区在线| 国产探花极品一区二区| 久久香蕉精品热| 村上凉子中文字幕在线| 久久中文看片网| or卡值多少钱| 亚洲av免费高清在线观看| 97超级碰碰碰精品色视频在线观看| 天堂网av新在线| 国内精品一区二区在线观看| 精品电影一区二区在线| 搡老岳熟女国产| 欧美性猛交黑人性爽| 母亲3免费完整高清在线观看| 在线免费观看不下载黄p国产 | 99精品欧美一区二区三区四区| 制服丝袜大香蕉在线| 国产精华一区二区三区| 亚洲人成网站在线播| 少妇的逼水好多| 999久久久精品免费观看国产| 国产成人系列免费观看| 有码 亚洲区| 午夜福利在线观看吧| 日本黄色片子视频| 亚洲在线观看片| 国产黄片美女视频| av天堂在线播放| 一二三四社区在线视频社区8| 欧美zozozo另类| 国产成人av教育| 精品人妻一区二区三区麻豆 | 韩国av一区二区三区四区| 99久久久亚洲精品蜜臀av| 久久精品夜夜夜夜夜久久蜜豆| 国产精品三级大全| 国模一区二区三区四区视频| 色综合欧美亚洲国产小说| 精品国内亚洲2022精品成人| 国语自产精品视频在线第100页| 午夜免费成人在线视频| 国产国拍精品亚洲av在线观看 | 日韩欧美国产在线观看| 9191精品国产免费久久| 麻豆久久精品国产亚洲av| 日本 av在线| 综合色av麻豆| 亚洲成av人片免费观看| svipshipincom国产片| 香蕉久久夜色| 欧美一区二区精品小视频在线| 美女黄网站色视频| 欧美在线一区亚洲| 国产中年淑女户外野战色| 精品久久久久久久久久免费视频| 婷婷精品国产亚洲av| 亚洲乱码一区二区免费版| 中文字幕av成人在线电影| 日韩欧美一区二区三区在线观看| 3wmmmm亚洲av在线观看| 亚洲国产欧美人成| 欧美成狂野欧美在线观看| 国产高清videossex| 亚洲专区国产一区二区| 亚洲国产精品sss在线观看| 哪里可以看免费的av片| 亚洲精品久久国产高清桃花| 午夜福利成人在线免费观看| 好看av亚洲va欧美ⅴa在| 每晚都被弄得嗷嗷叫到高潮| 在线免费观看的www视频| 一夜夜www| 午夜激情欧美在线| 国产高清videossex| 老熟妇乱子伦视频在线观看| 人人妻人人澡欧美一区二区| 亚洲内射少妇av| 桃红色精品国产亚洲av| 老司机午夜福利在线观看视频| 成人精品一区二区免费| 国产99白浆流出| 日韩 欧美 亚洲 中文字幕| 欧美一区二区亚洲| 一进一出抽搐动态| 国产aⅴ精品一区二区三区波| 欧美日韩瑟瑟在线播放| 男女午夜视频在线观看| 丰满乱子伦码专区| 精品一区二区三区视频在线 | 天堂动漫精品| 99riav亚洲国产免费| 欧美成人免费av一区二区三区| 久久精品亚洲精品国产色婷小说| 男人舔奶头视频| 亚洲av日韩精品久久久久久密| 国产亚洲精品久久久com| 日韩欧美在线二视频| 成人三级黄色视频| 亚洲熟妇熟女久久| a在线观看视频网站| 色尼玛亚洲综合影院| 色老头精品视频在线观看| 日本五十路高清| 国产三级黄色录像| 久久久国产成人免费| 国产色爽女视频免费观看| 国产视频一区二区在线看| 搡女人真爽免费视频火全软件 | 99riav亚洲国产免费| 国产高清videossex| 欧美日韩综合久久久久久 | svipshipincom国产片| 波野结衣二区三区在线 | 成人av在线播放网站| 久久久精品大字幕| 精品一区二区三区视频在线观看免费| 美女高潮喷水抽搐中文字幕| 国产亚洲av嫩草精品影院| 国产aⅴ精品一区二区三区波| 又紧又爽又黄一区二区| 久久久成人免费电影| 免费人成视频x8x8入口观看| 日本三级黄在线观看| h日本视频在线播放| 久久这里只有精品中国| 亚洲欧美日韩卡通动漫| 国产v大片淫在线免费观看| 岛国在线观看网站| 亚洲成人精品中文字幕电影| 亚洲精品美女久久久久99蜜臀| 国产精品嫩草影院av在线观看 | 亚洲精品美女久久久久99蜜臀| 男人的好看免费观看在线视频| 国产av一区在线观看免费| 日本成人三级电影网站| 精品国产美女av久久久久小说| 最后的刺客免费高清国语| 国产在视频线在精品| 在线观看免费视频日本深夜| 色综合亚洲欧美另类图片| 十八禁人妻一区二区| 国产综合懂色| 动漫黄色视频在线观看| 国产高清视频在线观看网站| 亚洲最大成人手机在线| 国产精品一区二区免费欧美| 亚洲一区二区三区色噜噜| 国产欧美日韩精品亚洲av| 日本一二三区视频观看| 亚洲人成网站高清观看| 久久精品91蜜桃| 国产欧美日韩精品亚洲av| 久久久国产成人免费| 午夜免费男女啪啪视频观看 | av视频在线观看入口| 91麻豆av在线| 国产精品国产高清国产av| 久久婷婷人人爽人人干人人爱| а√天堂www在线а√下载| 亚洲久久久久久中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 草草在线视频免费看| 色精品久久人妻99蜜桃| 精品人妻一区二区三区麻豆 | 午夜影院日韩av| 欧美一区二区亚洲| 蜜桃亚洲精品一区二区三区| 日韩精品中文字幕看吧| 国产蜜桃级精品一区二区三区| 欧美成人一区二区免费高清观看| 真实男女啪啪啪动态图| 全区人妻精品视频| 男女床上黄色一级片免费看| 日韩av在线大香蕉| www.熟女人妻精品国产| 全区人妻精品视频| 99热只有精品国产| 熟女人妻精品中文字幕| 在线观看舔阴道视频| 中文字幕av成人在线电影| 超碰av人人做人人爽久久 | 国产视频一区二区在线看| 长腿黑丝高跟| 别揉我奶头~嗯~啊~动态视频| 色综合婷婷激情| 久久国产乱子伦精品免费另类| 99久久无色码亚洲精品果冻| 午夜福利成人在线免费观看| 在线a可以看的网站| 色吧在线观看| 亚洲人成电影免费在线| 亚洲午夜理论影院| 亚洲国产精品sss在线观看| 性色av乱码一区二区三区2| 老司机福利观看| 免费av不卡在线播放| 男女床上黄色一级片免费看| 国产美女午夜福利| 免费av观看视频| 3wmmmm亚洲av在线观看| 午夜福利在线观看吧| 看免费av毛片| 亚洲精品一区av在线观看| 欧美在线一区亚洲| 制服人妻中文乱码| 每晚都被弄得嗷嗷叫到高潮| 极品教师在线免费播放| 色综合欧美亚洲国产小说| 身体一侧抽搐| 俄罗斯特黄特色一大片| 色尼玛亚洲综合影院| 两个人的视频大全免费| 偷拍熟女少妇极品色| 国产高清videossex| 日韩免费av在线播放| 极品教师在线免费播放| 色尼玛亚洲综合影院| 一进一出抽搐gif免费好疼| 午夜亚洲福利在线播放| 国产高清videossex| 91字幕亚洲| 日本a在线网址| 亚洲最大成人中文| 欧美在线黄色| 嫁个100分男人电影在线观看| 怎么达到女性高潮| 欧美日本亚洲视频在线播放| 中亚洲国语对白在线视频| 国产探花在线观看一区二区| 村上凉子中文字幕在线| 午夜老司机福利剧场| 免费一级毛片在线播放高清视频| 亚洲精品日韩av片在线观看 | 日日摸夜夜添夜夜添小说| 国产三级在线视频| 日本 欧美在线| 国产精品电影一区二区三区| 午夜福利18| 国产精品久久视频播放| 美女被艹到高潮喷水动态| 无遮挡黄片免费观看| 一区二区三区免费毛片| 一级毛片高清免费大全| 久久这里只有精品中国| 18禁黄网站禁片午夜丰满| 亚洲最大成人中文| 国产高清视频在线观看网站| 最新在线观看一区二区三区| 日本成人三级电影网站| a级毛片a级免费在线| 法律面前人人平等表现在哪些方面| 两个人的视频大全免费| 91在线精品国自产拍蜜月 | 女生性感内裤真人,穿戴方法视频| 国产成人欧美在线观看| 亚洲中文字幕一区二区三区有码在线看| 老汉色av国产亚洲站长工具| 亚洲人与动物交配视频| 欧美最新免费一区二区三区 | 欧美午夜高清在线| 亚洲人与动物交配视频| 国产精品1区2区在线观看.| 19禁男女啪啪无遮挡网站| 欧美成人免费av一区二区三区| 国产爱豆传媒在线观看| 女人十人毛片免费观看3o分钟| 中文字幕久久专区| 欧美在线黄色| 久久久国产成人免费| 色老头精品视频在线观看| 国产一区二区激情短视频| 亚洲国产中文字幕在线视频| 91字幕亚洲| 久久精品91无色码中文字幕| 9191精品国产免费久久| 99热只有精品国产| 成人鲁丝片一二三区免费| 一区二区三区激情视频| 亚洲av电影在线进入| 丰满人妻熟妇乱又伦精品不卡| 99国产极品粉嫩在线观看| 亚洲五月天丁香| 99久久精品一区二区三区| 亚洲片人在线观看| 在线观看免费午夜福利视频| 国产黄片美女视频| 天堂动漫精品| 99久国产av精品| 国产免费一级a男人的天堂| 一个人免费在线观看的高清视频| 精品人妻一区二区三区麻豆 | 丰满的人妻完整版| 国产美女午夜福利| 欧美最新免费一区二区三区 | 国产精品一及| 在线看三级毛片| 成年女人看的毛片在线观看| 最近最新中文字幕大全免费视频| 成人一区二区视频在线观看| 欧美bdsm另类| 狂野欧美白嫩少妇大欣赏| 婷婷精品国产亚洲av| 黄色视频,在线免费观看| av黄色大香蕉| 女人高潮潮喷娇喘18禁视频| tocl精华| 亚洲一区二区三区不卡视频| av欧美777| 色av中文字幕| 美女大奶头视频| 色精品久久人妻99蜜桃| 91九色精品人成在线观看| 国产成年人精品一区二区| www.熟女人妻精品国产| av中文乱码字幕在线| 夜夜爽天天搞| 窝窝影院91人妻| 亚洲片人在线观看| x7x7x7水蜜桃| 夜夜看夜夜爽夜夜摸| 国内毛片毛片毛片毛片毛片| 国产在视频线在精品| 亚洲精品日韩av片在线观看 | 亚洲午夜理论影院| 亚洲,欧美精品.| 99久久综合精品五月天人人| 12—13女人毛片做爰片一| 欧美在线一区亚洲| 看黄色毛片网站| 美女黄网站色视频| 亚洲精品456在线播放app | 99热这里只有精品一区| 日韩人妻高清精品专区| 色精品久久人妻99蜜桃| 欧美日本视频| 亚洲av免费在线观看| 欧美成人免费av一区二区三区| 人人妻,人人澡人人爽秒播| 日本一二三区视频观看| 亚洲国产精品合色在线| 噜噜噜噜噜久久久久久91| av黄色大香蕉| 精品电影一区二区在线| 精品一区二区三区视频在线 | 欧美色视频一区免费| 伊人久久精品亚洲午夜| av在线天堂中文字幕| 国产欧美日韩精品一区二区| 国产成人福利小说| 国产在线精品亚洲第一网站| av专区在线播放| 每晚都被弄得嗷嗷叫到高潮| 18美女黄网站色大片免费观看| 欧美日韩瑟瑟在线播放| 精品免费久久久久久久清纯| 亚洲国产欧美网| 国产一区二区亚洲精品在线观看| 亚洲在线观看片| 中文字幕人成人乱码亚洲影| 精品一区二区三区视频在线 | 超碰av人人做人人爽久久 | 亚洲人成网站在线播放欧美日韩| 欧美性猛交╳xxx乱大交人| 国产 一区 欧美 日韩| av天堂在线播放| 99riav亚洲国产免费| 欧美日韩国产亚洲二区| 日韩欧美三级三区| 首页视频小说图片口味搜索| 欧美日韩黄片免| 午夜免费男女啪啪视频观看 | 国产黄色小视频在线观看| 国产欧美日韩精品一区二区| 国产成人福利小说| 国产高清三级在线| 久久99热这里只有精品18| 欧美另类亚洲清纯唯美| 亚洲精品成人久久久久久| 国产av在哪里看| 99精品久久久久人妻精品| 国产99白浆流出| 男插女下体视频免费在线播放| 免费av不卡在线播放| 一个人观看的视频www高清免费观看| 日韩国内少妇激情av| 日韩欧美免费精品| 欧美乱码精品一区二区三区| 90打野战视频偷拍视频| 午夜视频国产福利| 亚洲国产欧美人成| 国产成人欧美在线观看| 久久久久久久亚洲中文字幕 | 亚洲成人免费电影在线观看| 18禁裸乳无遮挡免费网站照片| 国产视频一区二区在线看| 免费看十八禁软件| 三级男女做爰猛烈吃奶摸视频| 亚洲aⅴ乱码一区二区在线播放| 少妇高潮的动态图| 欧美成人a在线观看| 国产精品98久久久久久宅男小说| 国产精品久久视频播放| 亚洲精品一卡2卡三卡4卡5卡| 欧美xxxx黑人xx丫x性爽| 九九在线视频观看精品| 亚洲精品久久国产高清桃花| 亚洲人成网站在线播放欧美日韩| 亚洲五月婷婷丁香| 又黄又粗又硬又大视频| 51午夜福利影视在线观看| 国产伦在线观看视频一区| 叶爱在线成人免费视频播放| 国产午夜精品论理片| 天美传媒精品一区二区| 久久精品国产亚洲av涩爱 | 亚洲精品粉嫩美女一区| 日韩大尺度精品在线看网址| 国产精品免费一区二区三区在线| 成人18禁在线播放| 久久精品国产自在天天线| 人妻久久中文字幕网| www日本黄色视频网| 99久国产av精品| 啦啦啦免费观看视频1| 性欧美人与动物交配| 久久欧美精品欧美久久欧美| 女警被强在线播放| 精品一区二区三区人妻视频| 高清毛片免费观看视频网站| 99视频精品全部免费 在线| 一边摸一边抽搐一进一小说| 女同久久另类99精品国产91| 在线观看午夜福利视频| 国产精品久久久久久久电影 | 久久午夜亚洲精品久久| 岛国在线观看网站| 69人妻影院| 麻豆成人av在线观看| 69av精品久久久久久| 亚洲av五月六月丁香网| 亚洲五月天丁香| 麻豆成人av在线观看| 夜夜爽天天搞| 一区二区三区激情视频| 婷婷六月久久综合丁香| 欧美极品一区二区三区四区| 波多野结衣高清无吗| 日韩有码中文字幕| 精品人妻偷拍中文字幕| 久久久色成人| 人妻夜夜爽99麻豆av| 麻豆成人午夜福利视频| 1000部很黄的大片| 成人av在线播放网站| 精品人妻偷拍中文字幕| 热99re8久久精品国产| 97碰自拍视频| 欧美色欧美亚洲另类二区| 欧美日本亚洲视频在线播放| 国产69精品久久久久777片| 久久人人精品亚洲av| 女人被狂操c到高潮| 日韩成人在线观看一区二区三区| 国产v大片淫在线免费观看| 午夜福利欧美成人| 亚洲avbb在线观看| 亚洲片人在线观看| 亚洲人成网站在线播| 他把我摸到了高潮在线观看| 成熟少妇高潮喷水视频| 又粗又爽又猛毛片免费看| 免费观看的影片在线观看| 好男人电影高清在线观看| 97超级碰碰碰精品色视频在线观看| 日韩成人在线观看一区二区三区| 国产精品久久电影中文字幕| 九九热线精品视视频播放| 久久国产精品影院| 午夜激情欧美在线| 91麻豆精品激情在线观看国产| 人人妻人人看人人澡| 亚洲精品久久国产高清桃花| 午夜影院日韩av| 亚洲在线观看片| 女人十人毛片免费观看3o分钟| 久久久久久人人人人人| 亚洲人成网站在线播| 欧美性猛交黑人性爽| 在线观看免费视频日本深夜| 欧美3d第一页| АⅤ资源中文在线天堂| 国产黄片美女视频| 国产日本99.免费观看| 99热只有精品国产| 在线十欧美十亚洲十日本专区| 成年人黄色毛片网站| 午夜福利欧美成人| 国产欧美日韩精品一区二区| 久久天躁狠狠躁夜夜2o2o| 日韩欧美国产一区二区入口| 九色国产91popny在线| 18禁美女被吸乳视频| 两个人的视频大全免费| 亚洲国产精品久久男人天堂| 久久久精品欧美日韩精品| 精品不卡国产一区二区三区| 欧美黄色淫秽网站| 亚洲熟妇熟女久久| 欧美最新免费一区二区三区 | 午夜激情欧美在线| 宅男免费午夜| www日本在线高清视频| 一进一出抽搐动态| 色吧在线观看| 国内精品一区二区在线观看| 午夜激情福利司机影院| 99热只有精品国产| 激情在线观看视频在线高清| 非洲黑人性xxxx精品又粗又长| 黄片小视频在线播放| 国产精品乱码一区二三区的特点| 18禁黄网站禁片午夜丰满| 色吧在线观看| 欧美日韩瑟瑟在线播放| 99热6这里只有精品| 深夜精品福利| 国产乱人伦免费视频| aaaaa片日本免费| 国产精品一区二区三区四区久久| 亚洲色图av天堂| 免费在线观看亚洲国产| 天堂动漫精品| 波野结衣二区三区在线 | 久久精品人妻少妇| 尤物成人国产欧美一区二区三区| 一夜夜www| 日韩精品青青久久久久久| 亚洲在线自拍视频| 免费人成视频x8x8入口观看| 国产野战对白在线观看| 成人欧美大片| 美女cb高潮喷水在线观看|