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

    泥石流啟動(dòng)臨界土體含水量及其預(yù)警應(yīng)用

    2014-09-21 19:33:29胡凱衡馬超
    關(guān)鍵詞:滲透系數(shù)泥石流降雨

    胡凱衡+馬超

    基金項(xiàng)目:中國(guó)科學(xué)院重點(diǎn)部署項(xiàng)目(KZZDEW0501);國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目(41030742);中國(guó)科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所青年百人團(tuán)隊(duì)項(xiàng)目(110900K235)

    摘要:傳統(tǒng)的泥石流預(yù)警方法多基于前期和實(shí)時(shí)降雨量等間接指標(biāo),但實(shí)際上直接影響泥石流啟動(dòng)的關(guān)鍵物理參數(shù)是土體含水量,通過(guò)分析土體含水量的變化來(lái)判斷泥石流啟動(dòng)更為直接可靠。首先定義了泥石流啟動(dòng)的臨界土體含水量的概念,然后基于國(guó)內(nèi)外泥石流啟動(dòng)的觀測(cè)試驗(yàn)數(shù)據(jù),采用逐步回歸分析方法,建立了臨界土體含水量與土體滲透系數(shù)、孔隙度和顆粒曲率系數(shù)的經(jīng)驗(yàn)關(guān)系,進(jìn)而提出一種基于臨界土體含水量和實(shí)時(shí)降雨的泥石流預(yù)警方法。最后,以云南東川蔣家溝1999年7月16日發(fā)生的一場(chǎng)泥石流為實(shí)例進(jìn)行演算和驗(yàn)證。結(jié)果表明:該方法在可靠性和準(zhǔn)確性上優(yōu)于傳統(tǒng)利用臨界線(xiàn)和暴發(fā)線(xiàn)判別泥石流的預(yù)測(cè)模型。

    關(guān)鍵詞:泥石流;預(yù)警;臨界土體含水量;孔隙度;滲透系數(shù);曲率系數(shù);降雨

    中圖分類(lèi)號(hào):P642.23文獻(xiàn)標(biāo)志碼:A

    Critical Soil Moisture for Debris Flow Initiation and Its Application in Forecasting

    HU Kaiheng1,2, MA Chao1,2,3

    (1. Key Laboratory of Mountain Hazards and Earth Surface Processes, Chinese Academy of Sciences, Chengdu 610041,

    Sichuan, China; 2. Institute of Mountain Hazards and Environment, Chinese Academy of Sciences

    and Ministry of Water Resources, Chengdu 610041, Sichuan, China; 3. University of

    Chinese Academy of Sciences, Beijing 100049, China)

    Abstract: Most of traditional debris flow forecasting methods are based on indirect variables such as antecedent and realtime rainfalls. But the key factor influencing directly the debris flow initiation is soil moisture, which is more reliable for debris flow forecasting. Firstly, the concept of critical soil moisture for debris flow initiation was defined; secondly, based on the experimental and observation data of debris flow initiation at home and abroad, the empirical relationships between critical soil moisture and permeability coefficient, porosity, coefficient of paritcle curvature were obtained by the means of stepwise regression analysis; finally, the forcasting method for debris flow based on critical soil moisture and realtime rainfall was proposed. The method was tested and verified by the debris flow happened at Jiangjiagou of Dongchuan, Yunnan on 16 July, 1999. The results show that the reliability and accuracy of new forcasting method is better than the traditional forecast method for debris flow based on critical line and occurrence line.

    Key words: debris flow; forecasting; critical soil moisture; porosity; permeability coefficient; coefficient of curvature; rainfall

    0引言

    自然界中,泥石流往往由上游溝岸或坡面的固體物質(zhì)失穩(wěn)進(jìn)入溝道,并在溝道水流的動(dòng)力作用下形成。初始形成的泥石流規(guī)模不大,在運(yùn)動(dòng)過(guò)程中通過(guò)侵蝕、裹挾松散物質(zhì),規(guī)模才逐漸增長(zhǎng)[12]。

    一般而言,坡面土體要達(dá)到一定的含水量(體積比,下同)才能啟動(dòng)形成泥石流。坡面土體的含水量是反映泥石流形成的直接參數(shù)。目前常用的前期降雨和當(dāng)期降雨量預(yù)警指標(biāo),間接體現(xiàn)了降雨使土體含水量增加,抗剪強(qiáng)度和土體穩(wěn)定性降低,從而使泥石流發(fā)生的可能性增大[36]。國(guó)內(nèi)外許多泥石流預(yù)警模型多以“前期有效雨量雨強(qiáng)/雨強(qiáng)歷時(shí)”為基本模式[710]。但是前期有效雨量與流域下墊面特性、巖土體物理特性等有關(guān),其計(jì)算過(guò)程中涉及的遞減系數(shù)、前期降水衰減系數(shù)等重要常數(shù)的確定需要長(zhǎng)時(shí)間的土體含水量和降雨量觀測(cè),且不同學(xué)者使用的計(jì)算方法以及得到的結(jié)果可能不一樣[1114]。不同泥石流溝的有效前期雨量衰減規(guī)律也不同。因此,前期降雨量預(yù)警模式存在局地性強(qiáng)、預(yù)警精度不高、參數(shù)不易確定等缺點(diǎn)。

    針對(duì)“有效前期雨量雨強(qiáng)/雨強(qiáng)歷時(shí)”預(yù)警模型所存在的不足,筆者通過(guò)分析國(guó)內(nèi)外大量泥石流形成過(guò)程中的土體含水量變化數(shù)據(jù),直接以土體含水量為預(yù)警指標(biāo),提出泥石流啟動(dòng)臨界土體含水量的概念和經(jīng)驗(yàn)計(jì)算公式;經(jīng)過(guò)云南東川蔣家溝泥石流觀測(cè)數(shù)據(jù)的驗(yàn)證,所得到的方法和公式可應(yīng)用于泥石流預(yù)警工作中。

    1泥石流臨界土體含水量

    1.1泥石流形成過(guò)程中土體含水量變化

    Cannon等對(duì)美國(guó)加利福利亞、科羅拉多火災(zāi)地區(qū)大量野外土體含水量和泥石流暴發(fā)時(shí)間監(jiān)測(cè)表明,不同深度坡面土體(一般監(jiān)測(cè)最大深度為60 cm)在泥石流暴發(fā)時(shí)的含水量都未達(dá)到飽和[1517]。以加利福尼亞南部桑加布里埃爾山區(qū)2009年發(fā)生的泥石流為例,167 mm的累積雨量共造成7條流域產(chǎn)生泥石流,而流域中的土體含水量卻始終維持在22%,遠(yuǎn)小于飽和度40%[18]。這些地區(qū)的泥石流在較好的地表植被條件下(比如火災(zāi)過(guò)后植被逐漸恢復(fù)后),可由徑流觸發(fā)轉(zhuǎn)變?yōu)榛录ぐl(fā),1 h激發(fā)雨量也逐漸增大[19];蔣家溝暴雨泥石流的多年監(jiān)測(cè)也表明:降雨過(guò)程中,雨水在坡面松散土體最大可以下滲60~80 cm,斜坡坡腳位置的土體含水量最大,但是都未超過(guò)孔隙度 [2021];一些模型試驗(yàn)結(jié)果也表明,在土體啟動(dòng)時(shí)大部分土體處于未飽和狀態(tài)[2224]。這與傳統(tǒng)認(rèn)識(shí)的泥石流形成理論并不一致,尤其從非飽和土力學(xué)角度來(lái)看,泥石流啟動(dòng)的原因應(yīng)是土體飽和后短歷時(shí)雨強(qiáng)造成孔隙水壓力劇增、有效應(yīng)力下降或喪失。事實(shí)上,降雨過(guò)程中土體含水量空間分布不均勻,垂直深度上土體含水量是非線(xiàn)性關(guān)系。對(duì)孔隙度大、強(qiáng)度低、顆粒間黏結(jié)程度差的坡面松散物源體來(lái)說(shuō),土體失穩(wěn)不僅僅是由于含水量增加、抗剪強(qiáng)度減小、下滑力增加導(dǎo)致的,土體失穩(wěn)在泥石流形成過(guò)程中只是初步階段也是必不可少的階段,降雨下滲在土體中形成壤中流以及表面形成地表徑流也是非常重要的原因。在一些地區(qū)甚至徑流作用在泥石流形成作用中更是起主導(dǎo)作用。例如,根據(jù)蔣家溝流域坡面徑流和溝道匯流匯集過(guò)程的觀測(cè),大多是清水變成泥石流體的過(guò)程,在很多情況下,在坡面上就初步形成了小泥石流體或漿體[25]。美國(guó)加利福尼亞南部、科羅拉多地區(qū)火災(zāi)后的泥石流多是流域坡面土體含水量達(dá)到一定值后,由峰值降雨時(shí)段在引發(fā)一定水力條件的坡面徑流激發(fā)形成的[26]。

    1.2臨界土體含水量定義和經(jīng)驗(yàn)計(jì)算公式

    無(wú)論是土力類(lèi)還是水力類(lèi)的泥石流,暴發(fā)時(shí)刻大多在峰值降雨時(shí)段附近。隨著雨水入滲,土體含水量逐漸增加,當(dāng)含水量達(dá)到某個(gè)臨界值時(shí),源區(qū)坡面土體達(dá)到極限平衡狀態(tài),或者坡面土體入滲和失水達(dá)到動(dòng)態(tài)平衡。后續(xù)的降雨強(qiáng)度如果能夠達(dá)到或者超過(guò)滲透率,坡面產(chǎn)生地表徑流而且土體失穩(wěn)下滑,形成最初的泥石流。這一臨界含水量稱(chēng)為泥石流啟動(dòng)的臨界土體含水量。國(guó)內(nèi)外研究中有關(guān)臨界土體含水量的概念已有所涉及。例如,在Brocca等提出的一些降雨徑流模型中,土體水分平衡主要考慮了降雨入滲率、雨強(qiáng)、蒸發(fā)率以及因壤中流和深層土體滲透的土體排水率,臨界土體含水量與飽和度并不一致[27]。楊大文等將遂川江流域?qū)崪y(cè)土體含水量與不同時(shí)間的雨量結(jié)合起來(lái),建立流域內(nèi)不同點(diǎn)土體飽和度和警戒雨量的關(guān)系[28]。該結(jié)果間接反映了含水量越大,警戒雨量越小。實(shí)測(cè)的土體飽和度中也反映了臨界土體含水量的概念。

    通過(guò)對(duì)國(guó)內(nèi)外大量降雨激發(fā)淺表層滑坡、泥石流過(guò)程中土體含水量變化與土體物理參數(shù)的分析(表1),發(fā)現(xiàn)泥石流暴發(fā)時(shí)的土體含水量與土體滲透系數(shù)、孔隙度以及顆粒級(jí)配存在正相關(guān)關(guān)系,用Matlab的Stepwise函數(shù)作交互式逐步回歸分析,得到多重線(xiàn)性關(guān)系式

    Wa =-1.12K+0.46n+0.12Cc-0.165(1)

    式中:Wa為臨界土體含水量;K為土體滲透系數(shù);n為土體孔隙度;Cc為土體顆粒的曲率系數(shù),Cc=d230/(d10d60),其中d10、d30和d60分別是土體顆粒質(zhì)量累計(jì)含量(質(zhì)量分?jǐn)?shù),下同)為10%、30%和60%時(shí)的粒徑。

    式(1)擬合所用數(shù)據(jù)如表1,判定系數(shù)為0.896 5,調(diào)整判定系數(shù)為0.883,均方根誤差為0.055 6。對(duì)該方程顯著性進(jìn)行F檢驗(yàn),查表得到F0.05(3,23)值為3028,遠(yuǎn)小于統(tǒng)計(jì)量F值(66.4),因此,回歸公式顯著。

    從式(1)可以看出,臨界土體含水量與滲透系數(shù)成負(fù)相關(guān)關(guān)系,與孔隙度和土體顆粒曲率系數(shù)成正相關(guān)關(guān)系。滲透系數(shù)越大,說(shuō)明土體在一定水力梯度下可以很快下滲并在短時(shí)間內(nèi)轉(zhuǎn)化為壤中流。因內(nèi)部壤中流流動(dòng)和深層滲漏的土體排水速度與雨強(qiáng)大小接近時(shí),土體含水量才達(dá)到臨界值,并短時(shí)間維持在相對(duì)穩(wěn)定的水平。孔隙度越大,說(shuō)明水、氣兩相占據(jù)土體內(nèi)部空間越大。一般而言,在非飽和階段,土體內(nèi)部顆粒之間具有一定的黏聚力和咬合力并使土體具備一定的抗剪強(qiáng)度,從而保持穩(wěn)定。國(guó)內(nèi)外許多研究證明土體啟動(dòng)時(shí)是非飽和的,土體內(nèi)部需要更多的水才能使抗剪強(qiáng)度下降,下滑力增加。因此,臨界土體含水量與孔隙度成正相關(guān)關(guān)系說(shuō)明水分在松散碎屑坡面體穩(wěn)定性中具有重要作用。土體顆粒的曲率系數(shù)反映了土顆粒粒徑分布曲線(xiàn)形態(tài)。從表1中曲率系數(shù)可以看出,泥石流源區(qū)粒徑級(jí)配累計(jì)曲線(xiàn)斜率比較連續(xù),細(xì)顆粒在泥石流源區(qū)土體仍占很大部分。盡管級(jí)配連續(xù),但是泥石流源區(qū)土體是不均勻的,存在不連續(xù)粒徑,且均勻系數(shù)變化較大。由此說(shuō)明臨界土體含水量與曲率系數(shù)成正相關(guān)關(guān)系,曲率系數(shù)越大,土體內(nèi)部存在不連續(xù)粒徑,土顆粒間的存水空間越多,臨界土體含水量越大。

    另外,式(1)中的3個(gè)變量綜合反映了流域內(nèi)土體的平均最大蓄水量。臨界土體含水量實(shí)際上與平均最大蓄水量物理涵義一致。超過(guò)臨界土體含水量或者平均最大蓄水量的降水將從地表流走,形成地表徑流,并聚集形成足夠水動(dòng)力條件的溝道水流。坡面匯流而來(lái)的初步小規(guī)模泥石流漿體與溝道水流

    表1泥石流形成時(shí)的土體含水量、滲透系數(shù)、孔隙度以及曲率系數(shù)

    Tab.1Soil Moisture, Permeability Coefficient, Porosity and

    Coefficient of Curvature when Debris Flow Formed

    編號(hào)WaK/(mm·s-1)nCc

    1

    2

    30.500 00.023 5100.753.333 333

    0.590 00.014 4930.763.333 333

    0.660 00.004 6670.693.333 333

    0.540 00.007 6190.723.333 333

    0.480 00.003 7040.763.333 333

    0.540 00.005 0720.763.333 333

    0.240 00.005 6260.542.000 000

    0.133 00.005 6260.450.408 333

    0.120 00.010 2600.580.888 889

    0.113 00.087 7600.580.888 889

    0.124 00.010 2600.580.888 889

    0.103 60.087 7600.580.888 889

    0.161 80.010 2600.580.888 889

    0.152 60.005 6260.450.408 333

    0.117 40.010 2600.580.888 889

    0.152 90.005 6260.451.125 000

    0.326 00.003 8190.691.481 481

    0.310 00.003 8190.691.481 481

    0.357 00.003 8190.691.481 481

    0.237 00.003 8190.691.481 481

    0.340 00.003 8190.691.481 481

    0.340 00.003 8190.691.481 481

    0.352 00.003 8190.691.481 481

    0.378 00.003 8190.691.481 481

    0.313 00.003 8190.691.481 481

    0.347 00.003 8190.691.481 481

    0.347 00.003 8190.691.481 481

    注:編號(hào)1的數(shù)據(jù)引自文獻(xiàn)[29];編號(hào)2的數(shù)據(jù)引自文獻(xiàn)[30];編號(hào)3的數(shù)據(jù)引自文獻(xiàn)[24]。

    匯集增大了其中的含沙量和攜帶固體物質(zhì)能力。溝道水流在流經(jīng)泥石流動(dòng)床時(shí)強(qiáng)烈侵蝕固體物質(zhì)并形成泥石流。由于式(1)中的3個(gè)參數(shù)是泥石流源區(qū)土體物理特征值,所以該公式主要針對(duì)泥石流形成區(qū)的土體,尤其是細(xì)顆粒含量多的泥石流物源體。

    2基于臨界土體含水量和實(shí)時(shí)降雨的預(yù)警方法

    根據(jù)臨界土體含水量的定義以及泥石流預(yù)警模式中有效前期降雨難以精確計(jì)算的問(wèn)題,筆者提出一種基于臨界土體含水量和實(shí)時(shí)降雨的泥石流預(yù)警方法(圖1),該方法的具體流程如下。

    圖1基于土體含水量和實(shí)時(shí)降雨的泥石流預(yù)警方法實(shí)施過(guò)程

    Fig.1Flow Chart of Debris Flow Forecasting System Based on Critical Soil Moisture and Realtime Rainfall

    (1)計(jì)算泥石流啟動(dòng)的臨界土體含水量。根據(jù)泥石流形成區(qū)的巖土體特征,利用式(1)計(jì)算對(duì)應(yīng)臨界土體含水量Wa。

    (2)計(jì)算當(dāng)前土體含水量與臨界土體含水量的差值。在得到臨界土體含水量的情況下,通過(guò)土壤含水量傳感器,以太陽(yáng)能板和蓄電池作為電源,運(yùn)用GPRS無(wú)線(xiàn)網(wǎng)絡(luò)信息傳輸和室內(nèi)數(shù)據(jù)接收終端,得到降雨過(guò)程中t時(shí)刻的土體含水量Wt。通過(guò)室內(nèi)數(shù)據(jù)處理系統(tǒng)反復(fù)計(jì)算差值含水量ΔW,其表達(dá)式為

    ΔW=Wa-Wt(2)

    (3)根據(jù)當(dāng)前雨強(qiáng)計(jì)算達(dá)到臨界土體含水量所需要的時(shí)間T。具體計(jì)算方法為:通過(guò)雨量傳感器,以無(wú)線(xiàn)網(wǎng)絡(luò)方式發(fā)出和接收雨量值,通過(guò)室內(nèi)系統(tǒng)判定t時(shí)刻雨強(qiáng)Rt和土體滲透系數(shù)K的大??;并選擇大于滲透系數(shù)的計(jì)算公式和小于滲透系數(shù)的計(jì)算公式對(duì)達(dá)到臨界土體含水量所需要的時(shí)間進(jìn)行計(jì)算。

    T=(Wa-Wt)/KRt>K

    (Wa-Wt)/RtRt≤K(4)

    (4)根據(jù)計(jì)算的T值進(jìn)行預(yù)警,并每隔一定時(shí)間重復(fù)第(3)、(4)步直至達(dá)到臨界土體含水量或降雨結(jié)束。

    如果T>0,則多通道數(shù)據(jù)反復(fù)確認(rèn)降雨過(guò)程中土體含水量是否達(dá)到臨界土體含水量,并每隔一定時(shí)間(比如10 s)重復(fù)第(2)、(3)步操作內(nèi)容。如果降雨一直持續(xù)且T≤0,那么安排現(xiàn)場(chǎng)查看或發(fā)出泥石流即將發(fā)生的警報(bào)。如果降雨結(jié)束,則停止泥石流預(yù)警。

    3實(shí)例分析

    筆者利用提出的基于臨界土體含水量和實(shí)時(shí)降雨的泥石流預(yù)警方法,對(duì)云南東川蔣家溝1999年7月16日泥石流進(jìn)行演算。

    3.1臨界土體含水量確定

    蔣家溝角礫土密度為1.954 g·cm-3,孔隙度n為0.381 8,滲透系數(shù)k為0.008 07 mm·s-1,d10為0.01 mm,d30為0.25 mm,d60為3 mm,角礫土顆粒級(jí)配曲線(xiàn)的曲率系數(shù)Cc為3.125(圖2)。由此得到蔣家溝泥石流暴發(fā)的臨界土體含水量為401%。陳曉清等在蔣家溝人工降雨激發(fā)滑坡失穩(wěn)形成泥石流的試驗(yàn)中發(fā)現(xiàn),盡管不同深度的土體含水量不同,但每一組試驗(yàn)中土體含水量最大值都介于45%和35%之間,且土體破壞前的土體含水量大致在40%上下劇烈波動(dòng)[23]。土體在降雨作用下達(dá)到破壞前狀態(tài),其含水量劇烈變動(dòng),但是變動(dòng)的幅度基本在40%左右。這與通過(guò)式(1)計(jì)算得到的土體含水量值基本相同。因此,蔣家溝源區(qū)土體含水量接近越該值,泥石流暴發(fā)的可能性就越大。

    圖2云南東川蔣家溝角礫土的顆粒級(jí)配曲線(xiàn)

    Fig.2Particle Grading Curve of Breccia Soil in Jiangjiagou of Dongchuan, Yunnan

    3.2演算過(guò)程

    根據(jù)東川泥石流觀測(cè)站提供的降雨數(shù)據(jù),激發(fā)蔣家溝1999年7月16日泥石流的降雨過(guò)程見(jiàn)圖3。

    圖31999年7月16日泥石流暴發(fā)前的降雨過(guò)程和前20 d的雨量過(guò)程

    Fig.3Rainfall Processes Before Occurrence of

    Debris Flow on 16 July, 1999 and During the 20 Days Before the Debris Flow

    (1)計(jì)算差值含水量。初始時(shí)土體含水量為4%,未達(dá)到臨界土體含水量(40.1%),因此,在此時(shí)沒(méi)有泥石流發(fā)生,差值含水量ΔW為36.1%。

    (2)比較雨強(qiáng)、滲透系數(shù)大小并計(jì)算每單位時(shí)間的土體含水量。蔣家溝角礫土的滲透系數(shù)為0008 07 mm·s-1,相當(dāng)于每10 min降雨484 mm,與該溝的始發(fā)雨強(qiáng)(5 mm)非常接近。滲透系數(shù)0008 07 mm·s-1是經(jīng)過(guò)原位滲透試驗(yàn)得到的參數(shù),即土體達(dá)到穩(wěn)定滲透階段時(shí)的滲透系數(shù)。滲透系數(shù)是隨時(shí)間和土體含水量變化的,但無(wú)論滲透系數(shù)在降雨過(guò)程中如何變化,穩(wěn)定滲透階段的滲透系數(shù)在不同測(cè)試手段下差異不是很大。比如陳寧生等經(jīng)人工降雨試驗(yàn)測(cè)得降雨開(kāi)始時(shí)的土體初始滲透系數(shù)為0009 2 mm·s-1 [31],相當(dāng)于每10 min降雨552 mm的等效雨強(qiáng)。

    通過(guò)實(shí)際10 min雨量(R10)與滲透系數(shù)的對(duì)比,以單位面積和單位垂直深度的土體作為分析對(duì)象得到每單位10 min末的含水量

    Wt=WtΔt+ΔtRAHA(5)

    式中:WtΔt為tΔt時(shí)刻(對(duì)應(yīng)t時(shí)刻之初)的土體含水量;t為實(shí)際降雨記錄時(shí)刻;Δt為單位雨量記錄持續(xù)時(shí)間(這里為10 min);R為雨強(qiáng);

    A為土體面積;ΔtRA為時(shí)間Δt內(nèi)進(jìn)入土體的降水體積;HA為土體的總體積,H為一般含水量傳感器的探針長(zhǎng)度(60 mm)。

    (3)根據(jù)筆者提出的預(yù)警方法演算。由于蔣家溝沒(méi)有實(shí)測(cè)的泥石流暴發(fā)過(guò)程土體含水量變化,這里在得到臨界土體含水量、確定基本土體含水量和實(shí)際雨量過(guò)程后,根據(jù)提出的方法具體流程進(jìn)行演算。演算結(jié)果見(jiàn)表2。

    表2基于臨界土體含水量和實(shí)時(shí)降雨的泥石流預(yù)警方法演算過(guò)程

    Tab.2Calculation Process by the Forcasting

    Method for Debris Flow Based on Critical Soil Moisture and Realtime Rainfall

    編號(hào)時(shí)間段R10/mmWtΔtWa-WtΔtT/min

    121:50~22:000.50.040 00.361 0433.200 0

    222:00~22:101.70.048 30.352 7124.482 4

    322:10~22:201.60.076 70.324 3121.612 5

    422:20~22:301.10.103 30.297 7162.381 8

    522:30~22:400.80.121 70.279 3209.475 0

    622:40~22:500.90.135 00.266 0177.333 3

    722:50~23:001.30.150 00.251 0115.846 2

    823:00~23:102.00.171 70.229 368.790 0

    923:10~23:204.40.205 00.196 026.727 3

    1023:20~23:301.40.278 30.122 752.585 7

    1123:30~23:401.10.301 70.099 354.163 6

    1223:40~23:503.00.320 00.081 016.200 0

    1323:50~次日00:002.30.370 00.031 08.087 0

    14次日00:00~次日00:101.30.408 3<0.000 0<0.000 0

    注:次日01:12:34時(shí)刻,監(jiān)測(cè)到泥石流;R10值始終小于滲透系數(shù)。

    3.3精度比較

    將本文提出的預(yù)警方法與蔣家溝泥石流預(yù)報(bào)臨界線(xiàn)和暴發(fā)線(xiàn)判別式進(jìn)行結(jié)果對(duì)比,得到以下結(jié)果

    R10=5.5-0.098(Pa0+h)>0.5 mm(6)

    R10=6.9-0.123(Pa0+h)>1.0 mm(7)

    式中:Pa0為泥石流暴發(fā)前某一天的指數(shù);h為泥石流暴發(fā)前的當(dāng)日降雨量。

    臨界線(xiàn)式(6)的物理意義是:在10 min雨強(qiáng)大于0.5 mm的降水過(guò)程中,某10 min降水量只要等于5.5-0.098(Pa0+h),則蔣家溝泥石流就可能暴發(fā);暴發(fā)線(xiàn)式(7)的物理意義為:在10 min雨強(qiáng)大于1 mm的降水過(guò)程中,某10 min降水量只要等于69-0.123(Pa0+h),則蔣家溝就會(huì)暴發(fā)泥石流。

    利用該次泥石流過(guò)程之前的降雨量過(guò)程和蔣家溝前期雨量計(jì)算公式,得到臨界線(xiàn)和暴發(fā)線(xiàn)10 min雨量分別為2.26、2.86 mm(圖3)。從該次降雨過(guò)程來(lái)看,23:10~23:20時(shí)段4.4 mm的降雨是造成該次泥石流的主要原因。但泥石流暴發(fā)時(shí)并不與該次降雨過(guò)程的峰值雨量時(shí)段重合,滯后近1 h。表1演算結(jié)果表明,該時(shí)段土體含水量并未達(dá)到臨界土體含水量。而泥石流暴發(fā)時(shí)段雨強(qiáng)僅0.9 mm,大于臨界線(xiàn)而小于暴發(fā)線(xiàn)。從達(dá)到臨界土體含水量和泥石流暴發(fā)時(shí)間上來(lái)看,達(dá)到臨界土體含水量和泥石流暴發(fā)之間相差約1 h,而原方法可提前預(yù)警17~200 min[8]。因此,從臨界土體含水量結(jié)合實(shí)時(shí)降雨過(guò)程來(lái)判別泥石流的發(fā)生更準(zhǔn)確。筆者提出的基于臨界土體含水量和實(shí)時(shí)降雨的預(yù)警方法比傳統(tǒng)利用臨界線(xiàn)和暴發(fā)線(xiàn)判別泥石流的物理意義更明確,方法更可靠。

    4結(jié)語(yǔ)

    (1)在泥石流形成過(guò)程中,一般理論認(rèn)為前期降雨使源區(qū)土體飽和,短歷時(shí)雨強(qiáng)造成飽和后的土體產(chǎn)生高孔隙水壓力使土體失穩(wěn)并轉(zhuǎn)化為泥石流。但是國(guó)內(nèi)外大量泥石流形成過(guò)程監(jiān)測(cè)表明,源區(qū)坡面土體降雨過(guò)程和泥石流形成過(guò)程中土體都未達(dá)到飽和,且土體含水量存在一個(gè)臨界值。由此,本文提出了臨界土體含水量的概念,并通過(guò)擬合國(guó)內(nèi)外野外監(jiān)測(cè)數(shù)據(jù)得到計(jì)算臨界土體含水量的經(jīng)驗(yàn)公式。

    (2)基于臨界土體含水量的概念和經(jīng)驗(yàn)公式,發(fā)展了一種基于臨界土體含水量和實(shí)時(shí)降雨的泥石流預(yù)警方法。通過(guò)云南東川蔣家溝1999年7月16日暴發(fā)的泥石流的實(shí)際觀測(cè)資料,對(duì)該方法進(jìn)行了實(shí)例演算。演算結(jié)果表明,該場(chǎng)泥石流暴發(fā)時(shí)刻并未與峰值降雨時(shí)段重合,而是在達(dá)到臨界土體含水量后約1 h。

    (3)由于臨界土體含水量計(jì)算公式是經(jīng)驗(yàn)性的,在后續(xù)研究中有必要從土體降雨入滲以及激發(fā)坡面土體失穩(wěn)的物理過(guò)程并結(jié)合泥石流形成區(qū)監(jiān)測(cè)進(jìn)行深入研究,提出更具有物理意義的臨界土體含水量概念,建立以水文學(xué)、水力學(xué)、泥沙運(yùn)動(dòng)學(xué)等為基礎(chǔ)的預(yù)警模型。

    參考文獻(xiàn):

    References:

    [1]崔鵬.泥石流起動(dòng)條件及機(jī)理的實(shí)驗(yàn)研究[J].科學(xué)通報(bào),1991,36(21):16501652.

    CUI Peng.Experiment Study on the Mechanism and Condition of Starting Up of Debris Flow[J].Chinese Science Bulletin,1991,36(21):16501652.

    [2]崔鵬.中國(guó)山地災(zāi)害研究進(jìn)展與未來(lái)應(yīng)關(guān)注的科學(xué)問(wèn)題[J].地理科學(xué)進(jìn)展,2014,33(2):145152.

    CUI Peng.Progress and Prospects in Research on Mountain Hazards in China[J].Progress in Geography,2014,33(2):145152.

    [3]CAINE N.The Rainfall Intensity:Duration Control of Shallow Landslides and Debris Flows[J].Geografiska Annaler.Series A,Physical Geography,1980,62(1/2):2327.

    [4]WIECZOREK G F,GLADE T.Climaticc Factors Influencing Occurrence of Debris Flows[M]∥JAKOB M,HUNGR O.Debrisflow Hazards and Related Phenomena.Berlin:Springer,2005:325362.

    [5]馬超,胡凱衡,宋國(guó)虎,等.汶川地震災(zāi)區(qū)帽殼子滑坡形成泥石流的過(guò)程和特征[J].地球科學(xué)與環(huán)境學(xué)報(bào),2013,35(4):98103.

    MA Chao,HU Kaiheng,SONG Guohu,et al.Processes and Characteristics of Debris Flows Induced by Maoqiaozi Landslide in Wenchuan Earthquake Stricken Area[J].Journal of Earth Sciences and Environment,2013,35(4):98103.

    [6]謝洪,劉維明,趙晉恒,等.四川石棉2012年“7·14”唐家溝泥石流特征[J].地球科學(xué)與環(huán)境學(xué)報(bào),2013,35(4):9097.

    XIE Hong,LIU Weiming,ZHAO Jinheng,et al.Characteristics of Tangjiagou Debris Flow in Shimian of Sichuan in July 14,2012[J].Journal of Earth Sciences and Environment,2013,35(4):9097.

    [7]譚萬(wàn)沛,王成華,姚令侃,等.暴雨泥石流滑坡的區(qū)域預(yù)測(cè)與預(yù)報(bào):以攀西地區(qū)為例[M].成都:四川科學(xué)技術(shù)出版社,1994.

    TAN Wanpei,WANG Chenghua,YAO Lingkan,et al.Regional Forecasting and Predicting of Rainfall Induced Debris Flows and Landslides:Take the Western Panzhihua as an Example[M].Chengdu:Sichuan Science and Technology Press,1994.

    [8]陳景武.降雨預(yù)報(bào)泥石流的原理及方法[C]∥中國(guó)科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所.第二屆全國(guó)泥石流學(xué)術(shù)會(huì)議論文集.北京:科學(xué)出版社,1989:8490.

    CHEN Jingwu.The Principle and Method of Debris Flow Forecast Based on Rainfall[C]∥Institute of Mountain Hazards and Environment,Chinese Academy of Sciences and Ministry of Water Resources.Proceeding of the Second National Debris Flow Conference.Beijing:Science Press,1989:8490.

    [9]GUZZETTI F,PERUCCACCI S,ROSSI M,et al.The Rainfall Intensityduration Control of Shallow Landslides and Debris Flows:An Update[J].Landslides,2008,5(1):317.

    [10]戚國(guó)慶,黃潤(rùn)秋.泥石流成因機(jī)理的非飽和土力學(xué)理論研究[J].中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào),2003,14(3):1215.

    QI Guoqing,HUANG Runqiu.Study on Genetic and Mechanical Anlysis of Debris Flow Based on Unsaturated Soils Mechanics[J].The Chinese Journal of Geological Hazard and Control,2003,14(3):1215.

    [11]崔鵬,楊坤,陳杰.前期降雨對(duì)泥石流形成的貢獻(xiàn)——以蔣家溝泥石流形成為例[J].中國(guó)水土保持科學(xué),2003,1(1):1115.

    CUI Peng,YANG Kun,CHEN Jie.Relationship Between Occurrence of Debris Flow and Antecedent Precipitation—Taking the Jiangjia Gully as an Example[J].Science of Soil and Water Conservation,2003,1(1):1115.

    [12]韋方強(qiáng),胡凱衡,陳杰.泥石流預(yù)報(bào)中前期有效降水量的確定[J].山地學(xué)報(bào),2005,23(4):453457.

    WEI Fangqiang,HU Kaiheng,CHEN Jie.Determination of Effective Antecedent Rainfall for Debris Flow Forecast[J].Journal of Mountain Science,2005,23(4):453457.

    [13]李鐵鋒,叢威青.基于Logistic回歸及前期有效雨量的降雨誘發(fā)型滑坡預(yù)測(cè)方法[J].中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào),2006,17(1):3335.

    LI Tiefeng,CONG Weiqing.A Method for Rainfallinduced Landslides Prediction Based on Logistic Regression and Effective Antecedent Rainfall[J].The Chinese Journal of Geological Hazard and Control,2006,17(1):3335.

    [14]王裕宜,鄒仁元,李昌志.泥石流土體侵蝕與始發(fā)雨量的相關(guān)性研究[J].土壤侵蝕與水土保持學(xué)報(bào),1999,5(6):3438.

    WANG Yuyi,ZOU Renyuan,LI Changzhi.Study on Relationship Between Erosion of Debris Flows and Critical Rain Quantity[J].Journal ofSoil Erosion and Soil and Water Conservation,1999,5(6):3438.

    [15]CANNON S H,KIRKHAM R M,PARISE M.Wildfirerelated Debrisflow Initiation Processes,Storm King Mountain,Colorado[J].Geomorphology,2001,39(3/4):171188.

    [16]SANTI P M,DEWOLFE V G,HIGGINS J D,et al.Sources of Debris Flow Material in Burned Areas[J].Geomorphology,2008,96(3/4):310321.

    [17]CANNON S H,GARTNER J E,WILSON R C,et al.Storm Rainfall Conditions for Floods and Debris Flows from Recently Burned Areas in Southwestern Colorado and Southern California[J].Geomorphology,2008,96(3/4):250269.

    [18]KEAN J W,STALEY D M,CANNON S H.In Situ Measurements of Postfire Debris Flows in Southern California:Comparisons ofthe Timing and Magnitude of 24 Debrisflow Events with Rainfall and Soil Moisture Conditions[J].Journal of Geophysical Research:Earth Surface,2011,116(F4):121.

    [19]MEYER G A,PIERCE J L,WOOD S H,et al.Fire,Storms,and Erosional Events in the Idaho Batholith[J].Hydrological Processes,2001,15(15):30253038.

    [20]吳積善,康志成,田連權(quán),等.云南蔣家溝泥石流觀測(cè)研究[M].北京:科學(xué)出版社,1990.

    WU Jishan,KANG Zhicheng,TIAN Lianquan,et al.Observation Investigations on Debris Flows in Jiangjia Ravine,Yunnan Province[M].Beijing:Science Press,1990.

    [21]王裕宜,詹錢(qián)登,嚴(yán)壁玉.泥石流體結(jié)構(gòu)和流變特性[M].長(zhǎng)沙:湖南科學(xué)技術(shù)出版社,2001.

    WANG Yuyi,ZHAN Qiandeng,YAN Biyu.Debris flow Structure and Rheology[M].Changsha:Hunan Science and Technology Press,2001.

    [22]PONZIANI F,PANDOLFO C,STELLUTI M,et al.Assessment of Rainfall Thresholds and Soil Moisture Modeling for Operational Hydrogeological Risk Prevention in the Umbria Region(Central Italy)[J].Landslides,2012,9(2):229237.

    [23]陳曉清,崔鵬,馮自立,等.滑坡轉(zhuǎn)化泥石流起動(dòng)的人工降雨試驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào),2006,25(1):106116.

    CHEN Xiaoqing,CUI Peng,FENG Zili,et al.Artificial Rainfall Experimental Study on Landslide Transition to Debris Flow[J].Chinese Journal of Rock Mechanics and Engineering,2006,

    25(1):106116.

    [24]CHAE B G,KIM M I.Suggestion of a Method for Landslide Early Warning Using the Change in the Volumetric Water Content Gradient Due to Rainfall Infiltration[J].Environmental Earth Sciences,2012,66(7):19731986.

    [25]李椷,吳濟(jì)難.云南東川蔣家溝泥石流形成條件的初步分析[C]∥中國(guó)科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所.泥石流論文集(1).重慶:科學(xué)技術(shù)出版社重慶分社,1981:8792.

    LI Jian,WU Jinan.Preliminary Analysis on the Formation Conditions of Debris Flows in Jiangjia Ravine,Dongchuan,Yunnan Province[C]∥Institute of Mountain Hazards and Environment,Chinese Academy of Sciences and Ministry of Water Resources.Proceedings of Debris Flows:The First Volume.Chongqing:Chongqing Branches of Science and Technology Press,1981:8792.

    [26]LARSEN I J,MACDONALD L H,BROWN E,et al.Causes of Postfire Runoff and Erosion:Water Repellency,Cover,or Soil Sealing?[J].Soil Science Society of America Journal,2007,73(4):13931407.

    [27]BROCCA L,BARBETTA S,MELONE F,et al.A Continuous Rainfallrunoff Model Derived from Investigations in a Small Experimental Basin[C]∥SCHUMANN S A,HOLKO L.Status and Perspectives of Hydrology in Small Basins.GoslarHahnenklee:IAHS Press,2010:179185.

    [28]楊大文,龔偉,劉志雨,等.基于分布式模型土壤含水量評(píng)估的山洪預(yù)警指標(biāo)體系[C]∥中國(guó)水利學(xué)會(huì).中國(guó)水利學(xué)會(huì)2010學(xué)術(shù)年會(huì)論文集:上冊(cè).鄭州:黃河水利出版社,2010:464473.

    YANG Dawen,GONG Wei,LIU Zhiyu,et al.An Fooding Forecasting Method Based Distributed Model of Estimating Soil Moisutre[C]∥Chinese Hydraulic Engineering Society.2010 Annual Conference Proceedings of Chinese Hydraulic Engineering Society:The First Volume.Zhengzhou:Yellow River Water Conservancy Press,2010:464473.

    [29]GRECO R,GUIDA A,DAMIANO E,et al.Soil Water Content and Suction Monitoring in Model Slopes for Shallow Flowslides Early Warning Applications[J].Physics and Chemistry of the Earth,Parts A/B/C,2010,35(3/4/5):127136.

    [30]COE J A,KINNER D A,GODT J W.Initiation Conditions for Debris Flows Generated by Runoff at Chalk Cliffs,Central Colorado[J].Geomorphology,2008,96(3/4):270297.

    [31]陳寧生,張軍.泥石流源區(qū)弱固結(jié)礫石土的滲透規(guī)律[J].山地學(xué)報(bào),2001,19(1):169171.

    CHEN Ningsheng,ZHANG Jun.The Research of Permeability on Lose Gravelly Soil in Debris Flow Original Area[J].Journal of Mountain Science,2001,19(1):169171.

    CHEN Xiaoqing,CUI Peng,FENG Zili,et al.Artificial Rainfall Experimental Study on Landslide Transition to Debris Flow[J].Chinese Journal of Rock Mechanics and Engineering,2006,

    25(1):106116.

    [24]CHAE B G,KIM M I.Suggestion of a Method for Landslide Early Warning Using the Change in the Volumetric Water Content Gradient Due to Rainfall Infiltration[J].Environmental Earth Sciences,2012,66(7):19731986.

    [25]李椷,吳濟(jì)難.云南東川蔣家溝泥石流形成條件的初步分析[C]∥中國(guó)科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所.泥石流論文集(1).重慶:科學(xué)技術(shù)出版社重慶分社,1981:8792.

    LI Jian,WU Jinan.Preliminary Analysis on the Formation Conditions of Debris Flows in Jiangjia Ravine,Dongchuan,Yunnan Province[C]∥Institute of Mountain Hazards and Environment,Chinese Academy of Sciences and Ministry of Water Resources.Proceedings of Debris Flows:The First Volume.Chongqing:Chongqing Branches of Science and Technology Press,1981:8792.

    [26]LARSEN I J,MACDONALD L H,BROWN E,et al.Causes of Postfire Runoff and Erosion:Water Repellency,Cover,or Soil Sealing?[J].Soil Science Society of America Journal,2007,73(4):13931407.

    [27]BROCCA L,BARBETTA S,MELONE F,et al.A Continuous Rainfallrunoff Model Derived from Investigations in a Small Experimental Basin[C]∥SCHUMANN S A,HOLKO L.Status and Perspectives of Hydrology in Small Basins.GoslarHahnenklee:IAHS Press,2010:179185.

    [28]楊大文,龔偉,劉志雨,等.基于分布式模型土壤含水量評(píng)估的山洪預(yù)警指標(biāo)體系[C]∥中國(guó)水利學(xué)會(huì).中國(guó)水利學(xué)會(huì)2010學(xué)術(shù)年會(huì)論文集:上冊(cè).鄭州:黃河水利出版社,2010:464473.

    YANG Dawen,GONG Wei,LIU Zhiyu,et al.An Fooding Forecasting Method Based Distributed Model of Estimating Soil Moisutre[C]∥Chinese Hydraulic Engineering Society.2010 Annual Conference Proceedings of Chinese Hydraulic Engineering Society:The First Volume.Zhengzhou:Yellow River Water Conservancy Press,2010:464473.

    [29]GRECO R,GUIDA A,DAMIANO E,et al.Soil Water Content and Suction Monitoring in Model Slopes for Shallow Flowslides Early Warning Applications[J].Physics and Chemistry of the Earth,Parts A/B/C,2010,35(3/4/5):127136.

    [30]COE J A,KINNER D A,GODT J W.Initiation Conditions for Debris Flows Generated by Runoff at Chalk Cliffs,Central Colorado[J].Geomorphology,2008,96(3/4):270297.

    [31]陳寧生,張軍.泥石流源區(qū)弱固結(jié)礫石土的滲透規(guī)律[J].山地學(xué)報(bào),2001,19(1):169171.

    CHEN Ningsheng,ZHANG Jun.The Research of Permeability on Lose Gravelly Soil in Debris Flow Original Area[J].Journal of Mountain Science,2001,19(1):169171.

    CHEN Xiaoqing,CUI Peng,FENG Zili,et al.Artificial Rainfall Experimental Study on Landslide Transition to Debris Flow[J].Chinese Journal of Rock Mechanics and Engineering,2006,

    25(1):106116.

    [24]CHAE B G,KIM M I.Suggestion of a Method for Landslide Early Warning Using the Change in the Volumetric Water Content Gradient Due to Rainfall Infiltration[J].Environmental Earth Sciences,2012,66(7):19731986.

    [25]李椷,吳濟(jì)難.云南東川蔣家溝泥石流形成條件的初步分析[C]∥中國(guó)科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所.泥石流論文集(1).重慶:科學(xué)技術(shù)出版社重慶分社,1981:8792.

    LI Jian,WU Jinan.Preliminary Analysis on the Formation Conditions of Debris Flows in Jiangjia Ravine,Dongchuan,Yunnan Province[C]∥Institute of Mountain Hazards and Environment,Chinese Academy of Sciences and Ministry of Water Resources.Proceedings of Debris Flows:The First Volume.Chongqing:Chongqing Branches of Science and Technology Press,1981:8792.

    [26]LARSEN I J,MACDONALD L H,BROWN E,et al.Causes of Postfire Runoff and Erosion:Water Repellency,Cover,or Soil Sealing?[J].Soil Science Society of America Journal,2007,73(4):13931407.

    [27]BROCCA L,BARBETTA S,MELONE F,et al.A Continuous Rainfallrunoff Model Derived from Investigations in a Small Experimental Basin[C]∥SCHUMANN S A,HOLKO L.Status and Perspectives of Hydrology in Small Basins.GoslarHahnenklee:IAHS Press,2010:179185.

    [28]楊大文,龔偉,劉志雨,等.基于分布式模型土壤含水量評(píng)估的山洪預(yù)警指標(biāo)體系[C]∥中國(guó)水利學(xué)會(huì).中國(guó)水利學(xué)會(huì)2010學(xué)術(shù)年會(huì)論文集:上冊(cè).鄭州:黃河水利出版社,2010:464473.

    YANG Dawen,GONG Wei,LIU Zhiyu,et al.An Fooding Forecasting Method Based Distributed Model of Estimating Soil Moisutre[C]∥Chinese Hydraulic Engineering Society.2010 Annual Conference Proceedings of Chinese Hydraulic Engineering Society:The First Volume.Zhengzhou:Yellow River Water Conservancy Press,2010:464473.

    [29]GRECO R,GUIDA A,DAMIANO E,et al.Soil Water Content and Suction Monitoring in Model Slopes for Shallow Flowslides Early Warning Applications[J].Physics and Chemistry of the Earth,Parts A/B/C,2010,35(3/4/5):127136.

    [30]COE J A,KINNER D A,GODT J W.Initiation Conditions for Debris Flows Generated by Runoff at Chalk Cliffs,Central Colorado[J].Geomorphology,2008,96(3/4):270297.

    [31]陳寧生,張軍.泥石流源區(qū)弱固結(jié)礫石土的滲透規(guī)律[J].山地學(xué)報(bào),2001,19(1):169171.

    CHEN Ningsheng,ZHANG Jun.The Research of Permeability on Lose Gravelly Soil in Debris Flow Original Area[J].Journal of Mountain Science,2001,19(1):169171.

    猜你喜歡
    滲透系數(shù)泥石流降雨
    基于Origin的滲透系數(shù)衰減方程在地?zé)崴毓嘀械膽?yīng)用
    泥石流
    雜文月刊(2018年21期)2019-01-05 05:55:28
    多孔材料水滲透系數(shù)預(yù)測(cè)的隨機(jī)行走法
    輸水渠防滲墻及基巖滲透系數(shù)敏感性分析
    “民謠泥石流”花粥:唱出自己
    海峽姐妹(2017年6期)2017-06-24 09:37:36
    泥石流
    滄州市2016年“7.19~7.22”與“8.24~8.25”降雨對(duì)比研究
    紅黏土降雨入滲的定量分析
    河北平原新近系熱儲(chǔ)層滲透系數(shù)規(guī)律性分析
    機(jī)械班長(zhǎng)
    啦啦啦啦在线视频资源| 18禁观看日本| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | av女优亚洲男人天堂| 久久久久久久国产电影| 国产精品.久久久| 91国产中文字幕| 欧美精品国产亚洲| 日韩免费高清中文字幕av| 免费看av在线观看网站| 亚洲国产精品999| 伊人亚洲综合成人网| 搡女人真爽免费视频火全软件| 青青草视频在线视频观看| 亚洲人与动物交配视频| kizo精华| 亚洲av国产av综合av卡| 亚洲精品国产色婷婷电影| 老司机影院成人| 国产色爽女视频免费观看| 人人妻人人澡人人爽人人夜夜| 男人操女人黄网站| 黑人欧美特级aaaaaa片| a级毛色黄片| 赤兔流量卡办理| 亚洲精品aⅴ在线观看| 久久久久视频综合| 男人爽女人下面视频在线观看| 丝袜喷水一区| 国产女主播在线喷水免费视频网站| 亚洲怡红院男人天堂| 免费观看性生交大片5| 性色av一级| 人妻 亚洲 视频| 国产成人精品无人区| 亚州av有码| 日韩欧美精品免费久久| 久久精品国产a三级三级三级| 五月开心婷婷网| 日本午夜av视频| 人人妻人人添人人爽欧美一区卜| 中文字幕制服av| 国产在线视频一区二区| 哪个播放器可以免费观看大片| 久久久精品94久久精品| 亚洲人与动物交配视频| 99国产精品免费福利视频| av在线观看视频网站免费| 男女国产视频网站| a 毛片基地| 男女无遮挡免费网站观看| 永久网站在线| 多毛熟女@视频| 熟女人妻精品中文字幕| 国产不卡av网站在线观看| 亚洲精品视频女| 黑人猛操日本美女一级片| 亚洲欧美中文字幕日韩二区| 国产精品人妻久久久久久| 久久热精品热| 久久人人爽av亚洲精品天堂| 另类精品久久| 少妇高潮的动态图| 国产精品秋霞免费鲁丝片| 一二三四中文在线观看免费高清| 丁香六月天网| 男的添女的下面高潮视频| 午夜久久久在线观看| 下体分泌物呈黄色| 亚洲久久久国产精品| 男人操女人黄网站| 啦啦啦在线免费观看视频4| 中文字幕制服av| 女性生殖器流出的白浆| 视频区图区小说| 中文字幕高清在线视频| 国产在线观看jvid| 99国产精品一区二区蜜桃av | 激情在线观看视频在线高清 | 亚洲精华国产精华精| 色尼玛亚洲综合影院| 黄色成人免费大全| 他把我摸到了高潮在线观看 | 99久久人妻综合| 91成年电影在线观看| 午夜老司机福利片| 成人精品一区二区免费| av线在线观看网站| 亚洲一码二码三码区别大吗| 久久国产精品人妻蜜桃| 在线亚洲精品国产二区图片欧美| 免费观看人在逋| 人人妻人人爽人人添夜夜欢视频| 中文字幕人妻熟女乱码| 无限看片的www在线观看| 少妇精品久久久久久久| 亚洲一区二区三区欧美精品| 日日爽夜夜爽网站| 搡老岳熟女国产| 啪啪无遮挡十八禁网站| 日本av手机在线免费观看| 丝袜在线中文字幕| 男女之事视频高清在线观看| 成人影院久久| 国产深夜福利视频在线观看| 51午夜福利影视在线观看| 夜夜夜夜夜久久久久| 窝窝影院91人妻| 国产亚洲精品一区二区www | 久久影院123| 9色porny在线观看| 少妇 在线观看| 久久亚洲精品不卡| 狠狠婷婷综合久久久久久88av| 久久久国产成人免费| 黄色丝袜av网址大全| 少妇猛男粗大的猛烈进出视频| 麻豆成人av在线观看| 国产淫语在线视频| 性高湖久久久久久久久免费观看| 免费看十八禁软件| 久久午夜综合久久蜜桃| 激情在线观看视频在线高清 | 国产深夜福利视频在线观看| 天天影视国产精品| a级片在线免费高清观看视频| 男女免费视频国产| 日韩中文字幕视频在线看片| netflix在线观看网站| 亚洲全国av大片| 18禁国产床啪视频网站| 国内毛片毛片毛片毛片毛片| 久久久精品国产亚洲av高清涩受| 啦啦啦在线免费观看视频4| 国产亚洲一区二区精品| 国产欧美日韩一区二区三| 老鸭窝网址在线观看| 国产成人av教育| 久久热在线av| 欧美激情久久久久久爽电影 | 国产精品二区激情视频| 欧美精品人与动牲交sv欧美| 老司机午夜福利在线观看视频 | 精品一区二区三卡| 午夜福利在线免费观看网站| 亚洲熟女精品中文字幕| 少妇的丰满在线观看| 男男h啪啪无遮挡| 激情视频va一区二区三区| 亚洲精品久久成人aⅴ小说| 日韩欧美免费精品| 美女主播在线视频| 久久午夜综合久久蜜桃| 欧美精品亚洲一区二区| 巨乳人妻的诱惑在线观看| 99精国产麻豆久久婷婷| 精品国产乱码久久久久久小说| 国产激情久久老熟女| 女人爽到高潮嗷嗷叫在线视频| 建设人人有责人人尽责人人享有的| av欧美777| 久久精品国产亚洲av高清一级| 99热网站在线观看| 亚洲av成人一区二区三| 三上悠亚av全集在线观看| 91精品三级在线观看| 亚洲av国产av综合av卡| 91麻豆av在线| 欧美乱码精品一区二区三区| 飞空精品影院首页| 1024香蕉在线观看| 日韩人妻精品一区2区三区| av网站在线播放免费| 国产精品熟女久久久久浪| 涩涩av久久男人的天堂| 午夜福利在线免费观看网站| 视频区图区小说| 国产av国产精品国产| 欧美日韩亚洲高清精品| 免费在线观看完整版高清| 久久精品亚洲av国产电影网| 人人妻人人添人人爽欧美一区卜| 国产无遮挡羞羞视频在线观看| 亚洲国产毛片av蜜桃av| 成人亚洲精品一区在线观看| 如日韩欧美国产精品一区二区三区| 成人18禁在线播放| 免费在线观看视频国产中文字幕亚洲| 99国产极品粉嫩在线观看| 久久久国产精品麻豆| 电影成人av| 亚洲天堂av无毛| 久久人人爽av亚洲精品天堂| 亚洲精品成人av观看孕妇| 精品国产一区二区三区四区第35| 精品国产乱码久久久久久男人| 三级毛片av免费| 窝窝影院91人妻| 下体分泌物呈黄色| 丁香六月天网| 午夜福利免费观看在线| cao死你这个sao货| 可以免费在线观看a视频的电影网站| 91精品国产国语对白视频| 欧美日韩亚洲高清精品| 国产精品二区激情视频| 久久精品亚洲av国产电影网| 正在播放国产对白刺激| 国产精品国产高清国产av | 国产无遮挡羞羞视频在线观看| 日韩欧美三级三区| 欧美精品一区二区免费开放| 国产在线精品亚洲第一网站| 一二三四在线观看免费中文在| 在线观看人妻少妇| 国产精品av久久久久免费| 夜夜夜夜夜久久久久| 一本综合久久免费| 丁香六月天网| 大陆偷拍与自拍| 视频区图区小说| 亚洲国产欧美网| 极品少妇高潮喷水抽搐| 精品国产一区二区三区久久久樱花| 日本vs欧美在线观看视频| 丰满少妇做爰视频| 免费看a级黄色片| 少妇裸体淫交视频免费看高清 | 啦啦啦在线免费观看视频4| 999精品在线视频| www.999成人在线观看| 欧美成狂野欧美在线观看| 黄色a级毛片大全视频| 国产高清视频在线播放一区| 十八禁网站免费在线| 久9热在线精品视频| 啦啦啦视频在线资源免费观看| 久久久久久久国产电影| 日韩熟女老妇一区二区性免费视频| 国产国语露脸激情在线看| 啦啦啦在线免费观看视频4| 欧美+亚洲+日韩+国产| 亚洲人成77777在线视频| 三级毛片av免费| 国产成人免费无遮挡视频| 日韩一区二区三区影片| 他把我摸到了高潮在线观看 | 成年人午夜在线观看视频| 一个人免费看片子| 久久中文字幕人妻熟女| 久久久久视频综合| 老司机深夜福利视频在线观看| 国产深夜福利视频在线观看| 国产精品九九99| 精品一区二区三区四区五区乱码| 国产成人免费无遮挡视频| 亚洲欧美激情在线| 80岁老熟妇乱子伦牲交| 中文欧美无线码| 免费在线观看日本一区| 老熟妇仑乱视频hdxx| 女人久久www免费人成看片| 99精品久久久久人妻精品| 国产精品电影一区二区三区 | 麻豆乱淫一区二区| 考比视频在线观看| 国产精品自产拍在线观看55亚洲 | 久久久国产成人免费| 国产精品一区二区在线不卡| 国产日韩欧美视频二区| 欧美激情 高清一区二区三区| 久久亚洲真实| 欧美日韩亚洲国产一区二区在线观看 | 国产精品影院久久| 久久久水蜜桃国产精品网| 国产黄色免费在线视频| 精品一区二区三区av网在线观看 | 性高湖久久久久久久久免费观看| 女人精品久久久久毛片| 亚洲一区中文字幕在线| 久久国产精品男人的天堂亚洲| a在线观看视频网站| 久久久久久久久免费视频了| 三上悠亚av全集在线观看| 国内毛片毛片毛片毛片毛片| 高清av免费在线| 在线观看66精品国产| 欧美激情 高清一区二区三区| 免费人妻精品一区二区三区视频| 制服人妻中文乱码| 动漫黄色视频在线观看| 捣出白浆h1v1| 99久久精品国产亚洲精品| 又黄又粗又硬又大视频| 一夜夜www| 大香蕉久久成人网| 黄色a级毛片大全视频| 久久久久久久国产电影| 叶爱在线成人免费视频播放| 老司机福利观看| 老熟妇乱子伦视频在线观看| 国产精品久久久人人做人人爽| 99riav亚洲国产免费| 亚洲精品久久午夜乱码| 91成人精品电影| 51午夜福利影视在线观看| 中文字幕最新亚洲高清| 少妇粗大呻吟视频| 国产精品自产拍在线观看55亚洲 | av福利片在线| xxxhd国产人妻xxx| 亚洲全国av大片| 叶爱在线成人免费视频播放| √禁漫天堂资源中文www| 久久久久久免费高清国产稀缺| 色94色欧美一区二区| 人妻久久中文字幕网| 黄网站色视频无遮挡免费观看| 国产不卡av网站在线观看| 少妇裸体淫交视频免费看高清 | 首页视频小说图片口味搜索| 757午夜福利合集在线观看| 国产免费av片在线观看野外av| 国产欧美日韩一区二区精品| 热re99久久国产66热| 免费黄频网站在线观看国产| 老鸭窝网址在线观看| av国产精品久久久久影院| 国产在线一区二区三区精| 男女下面插进去视频免费观看| 午夜福利视频精品| 久久久精品免费免费高清| 久久这里只有精品19| 一二三四在线观看免费中文在| 免费高清在线观看日韩| 色尼玛亚洲综合影院| 亚洲伊人色综图| 两性夫妻黄色片| 成年人午夜在线观看视频| 十八禁网站免费在线| 亚洲精品美女久久av网站| 国产成人精品无人区| 又黄又粗又硬又大视频| 亚洲av第一区精品v没综合| 9191精品国产免费久久| 久久中文字幕人妻熟女| 我要看黄色一级片免费的| 精品国产一区二区三区久久久樱花| 亚洲成人手机| 叶爱在线成人免费视频播放| 超碰成人久久| 精品福利永久在线观看| 少妇 在线观看| 最近最新免费中文字幕在线| 国产老妇伦熟女老妇高清| 狠狠婷婷综合久久久久久88av| 久久精品国产综合久久久| 一个人免费在线观看的高清视频| 黑人操中国人逼视频| 国产欧美亚洲国产| 啦啦啦在线免费观看视频4| 亚洲熟妇熟女久久| 啦啦啦视频在线资源免费观看| 黄色视频不卡| 国产免费视频播放在线视频| 日韩视频一区二区在线观看| 中文字幕精品免费在线观看视频| 亚洲精品中文字幕一二三四区 | 99国产精品一区二区蜜桃av | 精品视频人人做人人爽| 黑人猛操日本美女一级片| 午夜福利一区二区在线看| 国产亚洲午夜精品一区二区久久| 热99久久久久精品小说推荐| 亚洲欧洲日产国产| 亚洲精品久久成人aⅴ小说| 午夜久久久在线观看| 久久精品aⅴ一区二区三区四区| 久久久久视频综合| 国产又爽黄色视频| 国产成人啪精品午夜网站| 国产主播在线观看一区二区| 中文字幕最新亚洲高清| 91精品三级在线观看| 脱女人内裤的视频| 91成年电影在线观看| 午夜日韩欧美国产| 免费观看人在逋| 亚洲av成人一区二区三| 黄色视频,在线免费观看| 成人精品一区二区免费| 国产精品国产av在线观看| 午夜成年电影在线免费观看| h视频一区二区三区| 色尼玛亚洲综合影院| av在线播放免费不卡| 久久久国产一区二区| 精品人妻在线不人妻| 亚洲欧洲精品一区二区精品久久久| 欧美午夜高清在线| 精品国产乱子伦一区二区三区| 捣出白浆h1v1| 精品国产国语对白av| 国产精品av久久久久免费| 欧美人与性动交α欧美软件| av网站免费在线观看视频| 久久国产精品影院| 亚洲色图 男人天堂 中文字幕| 国产一区二区三区视频了| 波多野结衣一区麻豆| 18禁观看日本| 国产单亲对白刺激| 在线观看免费视频网站a站| 成人18禁在线播放| 亚洲av美国av| 国产99久久九九免费精品| 亚洲av日韩精品久久久久久密| 午夜91福利影院| 久久中文字幕人妻熟女| 9热在线视频观看99| 最新的欧美精品一区二区| 999久久久精品免费观看国产| 99精品久久久久人妻精品| 99久久人妻综合| 中文字幕色久视频| av视频免费观看在线观看| 久久午夜综合久久蜜桃| 国产一区二区激情短视频| 99九九在线精品视频| 亚洲精品中文字幕一二三四区 | 久久久久国内视频| 天堂俺去俺来也www色官网| 男人操女人黄网站| 男女午夜视频在线观看| 国产男女内射视频| 国产精品二区激情视频| 亚洲午夜理论影院| 人成视频在线观看免费观看| 亚洲中文字幕日韩| 热99re8久久精品国产| 丝袜喷水一区| 欧美精品av麻豆av| 国产精品九九99| 美国免费a级毛片| 天天操日日干夜夜撸| 女同久久另类99精品国产91| 午夜精品久久久久久毛片777| 99久久人妻综合| 大码成人一级视频| 精品乱码久久久久久99久播| 男男h啪啪无遮挡| 天堂中文最新版在线下载| 成年版毛片免费区| 男人操女人黄网站| 亚洲美女黄片视频| 高清在线国产一区| 美女高潮到喷水免费观看| 国产成+人综合+亚洲专区| 国产精品久久久久久人妻精品电影 | 露出奶头的视频| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美精品高潮呻吟av久久| 18禁国产床啪视频网站| 少妇猛男粗大的猛烈进出视频| 少妇裸体淫交视频免费看高清 | 欧美精品av麻豆av| 亚洲精品粉嫩美女一区| 亚洲自偷自拍图片 自拍| 男女边摸边吃奶| 久久婷婷成人综合色麻豆| 亚洲人成伊人成综合网2020| av欧美777| 色视频在线一区二区三区| 久久国产精品大桥未久av| 国产精品国产高清国产av | 怎么达到女性高潮| 香蕉丝袜av| av免费在线观看网站| 两性夫妻黄色片| 午夜免费成人在线视频| 午夜久久久在线观看| 99国产精品免费福利视频| 亚洲欧美一区二区三区久久| 一级片免费观看大全| 老司机靠b影院| 两人在一起打扑克的视频| 人妻 亚洲 视频| 国产精品电影一区二区三区 | 性少妇av在线| 欧美变态另类bdsm刘玥| 窝窝影院91人妻| 免费观看人在逋| 国产视频一区二区在线看| 国产国语露脸激情在线看| 精品一区二区三卡| 看免费av毛片| 狠狠狠狠99中文字幕| 麻豆乱淫一区二区| 性少妇av在线| 久久天堂一区二区三区四区| 99国产精品99久久久久| 夜夜爽天天搞| 天天躁夜夜躁狠狠躁躁| www.精华液| 99在线人妻在线中文字幕 | 中文字幕制服av| 精品午夜福利视频在线观看一区 | 亚洲精品国产一区二区精华液| 91精品三级在线观看| 夜夜骑夜夜射夜夜干| 亚洲精品自拍成人| 乱人伦中国视频| 成人18禁在线播放| 90打野战视频偷拍视频| 欧美激情久久久久久爽电影 | 无遮挡黄片免费观看| 成人18禁在线播放| 亚洲av日韩精品久久久久久密| 美女国产高潮福利片在线看| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美av亚洲av综合av国产av| 欧美久久黑人一区二区| 99国产综合亚洲精品| 麻豆成人av在线观看| 日本wwww免费看| 80岁老熟妇乱子伦牲交| 大片免费播放器 马上看| 久久久久久久久久久久大奶| 亚洲国产成人一精品久久久| av网站在线播放免费| 成年人免费黄色播放视频| 757午夜福利合集在线观看| tocl精华| 99re在线观看精品视频| 日日爽夜夜爽网站| 电影成人av| 韩国精品一区二区三区| 正在播放国产对白刺激| 精品国产一区二区三区久久久樱花| 国产精品亚洲一级av第二区| 精品一区二区三区四区五区乱码| 欧美午夜高清在线| 天天躁日日躁夜夜躁夜夜| 黄片大片在线免费观看| 丁香六月天网| 人人妻人人添人人爽欧美一区卜| 亚洲国产精品一区二区三区在线| 老汉色∧v一级毛片| 老司机午夜福利在线观看视频 | 亚洲精品国产色婷婷电影| 成在线人永久免费视频| 天堂俺去俺来也www色官网| 亚洲成国产人片在线观看| 国产91精品成人一区二区三区 | 1024视频免费在线观看| 亚洲少妇的诱惑av| 高清在线国产一区| 黄色视频,在线免费观看| 国产日韩欧美视频二区| 亚洲精品乱久久久久久| 亚洲熟妇熟女久久| 国产精品一区二区精品视频观看| 久久这里只有精品19| 国产91精品成人一区二区三区 | 中国美女看黄片| 男女午夜视频在线观看| 国产精品一区二区免费欧美| 男女免费视频国产| 久久久久久人人人人人| 丰满人妻熟妇乱又伦精品不卡| 亚洲性夜色夜夜综合| 亚洲成人免费电影在线观看| av又黄又爽大尺度在线免费看| 亚洲精品粉嫩美女一区| 黄色成人免费大全| 久久久久久亚洲精品国产蜜桃av| 欧美黄色片欧美黄色片| 超碰成人久久| 国产极品粉嫩免费观看在线| 久久精品国产亚洲av香蕉五月 | 国产在线一区二区三区精| 久久午夜综合久久蜜桃| 精品免费久久久久久久清纯 | 国产欧美亚洲国产| 女性被躁到高潮视频| 国产精品美女特级片免费视频播放器 | 日本撒尿小便嘘嘘汇集6| 人妻一区二区av| 美女国产高潮福利片在线看| 精品国产一区二区久久| 人妻一区二区av| 曰老女人黄片| 亚洲伊人久久精品综合| 天堂动漫精品| 成年版毛片免费区| 夜夜夜夜夜久久久久| 国产在线视频一区二区| 黄色成人免费大全| 亚洲成国产人片在线观看| 国产精品九九99| 日韩一区二区三区影片| 色综合欧美亚洲国产小说| 90打野战视频偷拍视频| 动漫黄色视频在线观看| 国产激情久久老熟女| 99国产综合亚洲精品| 欧美日韩亚洲高清精品| 欧美精品av麻豆av| 狠狠狠狠99中文字幕| 韩国精品一区二区三区| 91老司机精品| 久久亚洲真实| 波多野结衣av一区二区av| 国产成人免费无遮挡视频| 国产黄频视频在线观看|