聶肖虎,謝璐瑤,楊冬,劉新,劉青
(1.中國石化銷售股份有限公司華中分公司, 湖北 武漢 430023; 2.紫光軟件系統(tǒng)有限公司, 北京 100084)
我國是世界上地質災害分布最嚴重的國家之一,據(jù)2019年的地質災害累計數(shù)據(jù)統(tǒng)計,我國共分布約25萬多個地質災害點,包括滑坡、斜坡、崩塌、泥石流、地面塌陷、地面沉降等(李媛等,2004),其中滑坡斜坡類地質災害數(shù)量占比最大,可達60%,其次為崩塌,約22%,我國地質災害類型數(shù)量占比見圖1。廣義的滑坡是指各種類型的巖土體沿斜坡向下運動的過程,通常包括泥石流和崩塌,我國學者多采用狹義滑坡定義,多指斜坡土層或堅硬的巖層,在重力作用下沿一個或多個滑動面向下移動的過程?;碌鹊刭|災害的發(fā)生一方面取決于基礎地質條件,比如地層巖性、坡度、高程、土壤等;另一方面又受到多種外界因素的誘發(fā)影響,如氣象、人類工程活動等(丁繼新等,2004;李長江等,2008)。在各種因素誘發(fā)的滑坡地質災害中,以降雨誘發(fā)的滑坡是分布最廣、頻次最高的地質災害(劉廣潤等,2002;張玉成等,2007)(以下簡稱降雨型滑坡)。據(jù)相關統(tǒng)計,我國已發(fā)生的突發(fā)性地質災害(滑坡、崩塌、泥石流)中,降雨誘發(fā)的災害數(shù)量可達90%。有學者對重慶地區(qū)滑坡災害進行統(tǒng)計,近50年來重慶地區(qū)由降雨誘發(fā)的滑坡數(shù)量約147個,占總滑坡數(shù)量的97%,且多數(shù)滑坡的發(fā)生與24小時內的降雨量有極大相關性(周國兵等,2003),因此開展降雨型滑坡地質災害的研究對于國家開展防災減災工作具有十分重要的意義(文海家等,2018)。
降雨型滑坡一直以來都是國內外地質災害研究的熱點領域,經過近幾十年的發(fā)展,其研究成果十分豐碩。從探索降雨誘發(fā)滑坡的作用機理,到研究與滑坡是否發(fā)生相關的具體臨界閾值,并將相關閾值模型用于滑坡預報預警實踐中,以及近些年來隨著RS、GIS技術的融入,使滑坡的早期識別和區(qū)域分級預警得以實現(xiàn)。降雨型滑坡研究內容見圖2。降雨誘發(fā)滑坡的作用機理主要體現(xiàn)在水體下滲過程中造成斜坡巖土體的物理性質及內部孔隙壓力的變化,使巖土體的抗剪切強度降低(李長江等,2008)。雨水在坡體中的入滲規(guī)律研究多是針對單體滑坡進行模型試驗及數(shù)值模擬,旨在揭示滑坡發(fā)生的內在規(guī)律,為不同滑坡的關鍵監(jiān)測指標提供參考,更偏于理論性質的研究,但其也是制定滑坡監(jiān)測方案的第一步。適用于報警的閾值研究采用的閾值參數(shù)根據(jù)不同區(qū)域的滑坡特點有所不同,暴雨頻次、降雨的周期變化、降雨歷時、降雨量以及降雨雨型與地質災害的發(fā)生存在著密切的聯(lián)系(馬力等,2002),該類研究通常對降雨滑坡歷史數(shù)據(jù)進行統(tǒng)計分析(高華喜和殷坤龍,2007),采用相應的模型進行閾值確定,也有部分學者將地質因素加入到閾值計算的模型中(劉艷輝和劉麗楠,2016),總的來說,閾值研究仍是降雨型滑坡的一個研究難點。但是目前這兩項研究一般是割裂進行的,未能共同作用于最終的預警預報工作中,如何將二者研究結合進行預警是今后降雨型滑坡研究的一個突破點。不論是單體滑坡的誘發(fā)機制研究還是區(qū)域性的降雨閾值研究,兩者在空間上的預測都略顯劣勢,但GIS技術的發(fā)展使得具有空間性質的柵格單元內的滑坡預測研究也隨之加強(彭令等,2010),為了進行更加精細化的空間地質災害預警,GIS是不可缺少的重要工具。
本文首先介紹降雨型滑坡的發(fā)育特征、雨水入滲滑坡的內在規(guī)律和相關模型,進而說明不同的降雨特征與滑坡之間的關系,并在此基礎上對目前的降雨型滑坡閾值研究進行論述和總結,并提出適合長輸管道地質災害預警工作的降雨閾值計算模型,最后指出降雨型滑坡研究目前存在的問題,并對今后的研究方向進行展望。
現(xiàn)有的研究多通過兩種方式探討降雨與滑坡之間的關系,一是研究滑坡內在的規(guī)律及穩(wěn)定性,利用各類監(jiān)測設備(水位計、孔壓計等)獲取滑坡體內各項參數(shù)在降雨條件下的變化情況,以此來建立雨水入滲模型及滑坡發(fā)生的內在機制;二是從外在的降雨特征著手,在收集充足的歷史降雨數(shù)據(jù)基礎上,將滑坡發(fā)生的時間點與之一一對應,研究滑坡發(fā)生時的不同降雨參數(shù)(雨量、雨強、歷時、雨型)的特征及閾值,并以此來對滑坡進行預報。
降雨誘發(fā)滑坡的實質是滲入斜坡的雨水通過與巖土體的相互作用,軟化滑面處的巖土體,降低巖土體的抗剪切強度,從而破壞邊坡的穩(wěn)定性,當巖土體的強度不能抵抗剪切力時,滑坡就發(fā)生了(楊文東,2006;李長江等,2008)。兩者相互作用的結果就是導致巖土體的內聚力和摩擦力下降(徐則民等,2005),具體誘發(fā)機理可分為四種:雨水入滲到斜坡巖土體內,使巖土體的強度減小而重量不斷增加,致使巖土體下滑力隨之增加;降雨入滲降低斜坡非飽和區(qū)域的基質吸力,致使構成斜坡的土體抗剪強度也隨之降低;降雨使得巖土體干濕交替,產生大量裂隙,為雨水入滲提供了更多的通道,加速了滑坡的發(fā)生;降雨抬升了地下水位,增加了對巖土體的浮托力(盛玲,2015)。
降雨觸發(fā)滑坡的內在因素包括構成斜坡體的物質條件(如巖質滑坡與土質滑坡)、斜坡體的坡度以及地質構造等(吳仁銑,2013),而雨水進入邊坡體的入滲率取決于降雨的方式、降雨強度以及斜坡巖土體的滲水物理屬性。不同物質組成的斜坡滲流過程也不相同,土質邊坡表現(xiàn)為由非飽和到飽和的變化過程,形成浸潤區(qū),抗剪強度降低;節(jié)理較發(fā)育、孔隙較大的巖質邊坡因雨水滲流在地下水位以上的區(qū)域形成上層滯水區(qū),孔隙壓力增大,降低抗滑力,降雨入滲下伏基巖類斜坡和土質斜坡示意見圖3。
國內外學者對于降雨入滲過程及機理進行了許多研究,提出了多種不同形式的入滲模型(朱文彬和劉寶琛,2002;謝劍明,2004)、過程模擬方法(楊文東,2006)及斜坡穩(wěn)定性評價方法(陳祖煜,2003)。降雨入滲實質上是水分在土壤飽氣帶中的運動(吳仁銑,2013),雨水滲入邊坡到達潛水面經歷了一個飽和—非飽和的滲流過程,從物理學角度得出邊坡入滲過程的非飽和滲流控制方程如下:
(1)
式(1)中:h在飽和區(qū),為壓力水頭,在非飽和區(qū),為基質吸力(非飽和帶中空氣壓力與水壓力之間的差值)的勢函數(shù),稱為基質勢;k(h)ij為滲透系數(shù),在非飽和區(qū)為基質勢h的函數(shù);S為貯水系數(shù);β在飽和區(qū),β=0,在非飽和區(qū),β=1;R為源匯項,為降雨入滲的供水強度;x1,x2,x3分別表示笛卡爾坐標系中的x(x1)軸,y(x2)軸,z(x3)軸為正向向上的鉛直軸。
Green-Ampt模型是一種簡化的入滲模型,用來確定瞬態(tài)入滲鋒面所到達的位置,可以看作是入滲問題的近似解,應用較為廣泛。但實際發(fā)生的滑坡機制是十分復雜的,Mein and Larson(1973)在該模型上根據(jù)降雨強度與土體入滲能力的關系進行改進,將雨水入滲分為2個階段,即自由入滲階段和積水入滲階段。國內學者的相關試驗表明,雨水入滲過程中還需考慮坡度因子,目前比較符合實際入滲情況的模型由楊文東(2006)提出,該模型在Green-Ampt模型的基礎上增加了入滲過程及斜坡效應。此外,楊文東(2006)除了對入滲機理和穩(wěn)定性進行分析外,還建立了模擬滑坡發(fā)生過程的模型,可以重現(xiàn)邊坡滑坡的動態(tài)破壞過程。
雨水改變了邊坡內的孔隙水壓力,從而降低邊坡穩(wěn)定性,在了解雨水入滲規(guī)律的基礎上,進而研究邊坡在降雨條件下的穩(wěn)定性?;路€(wěn)定性評價的方法包括極限平衡法、有限元強度折減法(銀曉鵬,2008)、模糊綜合評價(楊文東,2006)以及概率分析方法。降雨型滑坡穩(wěn)定性分析中,土體含水量和孔隙水壓力是兩個基本參數(shù),這是因為邊坡里有水的滲流作用,土的含水量增加,抗剪強度降低,致使邊坡極不穩(wěn)定。滲透力因素是否加入到穩(wěn)定性分析中則取決于研究對象的選擇,毛昶熙和段祥寶(2012)對地震和降雨等多個因素作用下的滑坡進行穩(wěn)定性計算時,為了更好地結合有限元法一次完成計算,加入了滲透力概念。陳祖煜(2003)考慮到滑裂面上的孔壓與坡面上的水壓力,認為不需要計算滲透力。在運用現(xiàn)場實測的土體含水量和孔隙水壓力對降雨型滑坡穩(wěn)定性進行評價時,不僅僅需要討論穩(wěn)定性分析物理模型的適用性,還需要對穩(wěn)定性分析輸入?yún)?shù)的合理性及其對研究目標的敏感性進行研究(Bordoni et al.,2015)。
1.2.1 時空分布特性
因降雨導致的滑坡等地質災害在空間分布上具有一定的區(qū)域性和同時性,且往往突發(fā)性強,災害程度大。國內已經有許多學者對降雨與滑坡等地質災害之間的時空分布關系進行了深入研究,并取得了非常多的成果。陳正洪和孟斌(1995)從幾十個典型事例出發(fā),分析了湖北省降雨型滑坡的時空分布特征,并且與降雨的時空分布規(guī)律進行相關分析,該研究促進了對二者之間對應關系的認識。劉艷輝等(2009)通過統(tǒng)計分析2006—2007年汛期全國地質害的實際發(fā)生情況得出地質災害的發(fā)育密度與年均雨量呈正比,在時間分布上與月降雨量、月暴雨日數(shù)的分布總體上具備一定的對應關系,且主要受典型強降雨過程的落區(qū)控制。林孝松和郭躍(2001)用統(tǒng)計分析的方法研究了滑坡與降雨的耦合關系,暴雨與滑坡的發(fā)生在時間上具有很好的一致性,通常也伴隨著一定的滯后效應,即滑坡發(fā)生于降雨之后的一段時間內,這就證實了滑坡的發(fā)生與降雨之間存在不穩(wěn)定周期性,且其規(guī)模與雨量之間存在一定的正相關。中國地質災害分布及中國年平均降雨量見圖4、圖5,年均降雨量大的區(qū)域也正是滑坡等地質災害分布密集區(qū)域。
1.2.2 降雨量、降雨強度、降雨歷時、雨型對滑坡的影響
國內外學者對降雨量與滑坡地質災害之間的關系進行了大量研究,結果表明滑坡地質災害除了受降雨在空間和時間上的分布影響外,降雨的雨量、雨強、歷時、雨型等特征也控制影響著滑坡的變形過程。
鐘蔭乾(1998)發(fā)現(xiàn)滑坡變形的過程與變形階段間的降雨強度變化的相關性較大,雨量越大的區(qū)域滑坡越發(fā)育,以新灘滑坡和黃臘石滑坡為例,說明滑坡的變形和劇滑與降雨的間斷性、時間性相關性較高。也有較多研究通過滑坡位移的監(jiān)測數(shù)據(jù)研究滑坡與降雨量之間的相關性,楊順泉(2002)對福建某滑坡體進行為期三年的監(jiān)測,將水平位移數(shù)據(jù)、垂直位移數(shù)據(jù)與降雨量進行回歸運算分析,結果顯示滑坡體各方位的位移量與降雨量相關性顯著,存在降雨量越大,位移變形越大的趨勢。
降雨強度指在一個給定時段內的累積降雨量或降雨速率(麻土華等,2011),暴雨引發(fā)的滑坡往往比長時間降雨型的臨界降雨量低,以降雨強度較大的深圳為例(高華喜和殷坤龍,2007),將滑坡發(fā)生次數(shù)分別與降雨強度、降雨量、降雨時間進行耦合分析,結果表明,滑坡與1~2天的降雨量(大暴雨)、 1~4天及一次降雨過程的累計降雨量的相關性系數(shù)均在0.8以上;在降雨時間上,暴雨發(fā)生的當天,便有45.3%的滑坡出現(xiàn),僅有13.3%的滑坡滯后于暴雨的3天以后,最晚在4天左右。
張勇等(2020)在分析降雨強度的基礎上,根據(jù)降雨的過程特征和降雨時間,對巴中市2007—2016年間降雨事件進行了細分,共劃分為6類降雨雨型,并分別統(tǒng)計不同雨型下發(fā)生的滑坡次數(shù),結果發(fā)現(xiàn)遞增型降雨下的滑坡占比最大,可達46.9%;單峰型降雨次之,約占25.02%。此外,對這兩類雨型下的滑坡進行穩(wěn)定性分析得知,其滑坡穩(wěn)定性在降雨后出現(xiàn)明顯降低,表明降雨雨型因子也對滑坡有著顯著影響。
降雨引發(fā)滑坡等地質災害過程是含有多個相互作用的復雜非線性關系系統(tǒng),我們將該過程看作一個系統(tǒng),影響滑坡的各類因素(降雨、巖土體、地形等)看作是該系統(tǒng)的輸入變量,那么滑坡的發(fā)生與否(發(fā)生與不發(fā)生)則是該系統(tǒng)的輸出變量,對輸入變量和輸出變量之間的非線性關系進行研究便是滑坡的預報模型(Crozier and Eyles,1980)。目前對于滑坡的預報模型研究主要分為兩類,一是與滑坡發(fā)生有關的閾值研究,分區(qū)域預測滑坡發(fā)生的時間段,偏時間預測;二是根據(jù)已知滑坡點繪制滑坡敏感性圖,一般不考慮誘發(fā)因素的影響,偏空間預測,因此,本文主要對前者進行介紹。根據(jù)有效應力原理,當斜坡巖土體的孔隙壓力到達失穩(wěn)狀態(tài)的臨界條件時,滑坡才能發(fā)生,而孔隙壓力與地下水的狀態(tài)相關,雨水下滲又增加了地下水的補給,從而誘發(fā)滑坡,因此,這是通過對降雨量的研究預報滑坡的理論基礎。通過閾值計算對滑坡進行預報是目前地質災害研究的熱點,現(xiàn)階段的閾值研究通常采用降雨、土壤含水量;或者其他水文地質參數(shù),這些參數(shù)往往具有區(qū)域性,因而降雨型滑坡預報幾乎都是根據(jù)區(qū)域特點計算降雨閾值而建立的。
降雨閾值是研究引發(fā)滑坡的降雨下限或臨界降雨,通過對誘發(fā)滑坡地質災害的歷史降雨事件進行統(tǒng)計,以降雨時間為橫坐標,降雨量為縱坐標,將兩項數(shù)據(jù)繪制在笛卡爾坐標、半對數(shù)坐標或者雙對數(shù)坐標中,然后利用相關的統(tǒng)計分析軟件(如R語言、SPSS軟件),對這些數(shù)據(jù)點的下限水平線進行擬合,即可得到該區(qū)域的降雨閾值示意圖(圖6)。該曲線中前期降雨強度小,從“0點”開始,降雨強度明顯增大直到滑坡發(fā)生。臨界降雨量則是指這段時間內的累計降雨量,閾值研究則是想通過降雨量變化的曲線來推斷滑坡發(fā)生的時間,建立一個經驗公式,可以用來作為滑坡啟動的降雨判據(jù)(陳洪凱等,2012),是目前降雨型滑坡還未得到標準答案的難點問題,通常采用過程降雨閾值模型、前期降雨閾值模型或者其他類型的組合模型進行研究。
2.1.1 過程降雨閾值
過程降雨閾值是對引發(fā)滑坡的降雨強度的閾值進行研究,通過對降雨歷史數(shù)據(jù)進行分析,得到其滑坡發(fā)生時所需要的降雨強度,即一定時間內的降雨速率,該時間段可以是短時,小時降雨量、日降雨量等;也可以是長時間觀測,以該階段的降雨量平均值進行表示。不同的時間尺度上,過程降雨閾值模型的確定也就相對復雜,且在不同的空間尺度上,閾值也不盡相同。目前包括四種模型:降雨強度—歷時(ID);使用年均降雨量(MAP)或其他參數(shù);過程累積雨量—歷時(ED);過程累積雨量—降雨強度(EI)。其中,降雨強度—歷時關系在世界各地降雨預報中使用頻率最高。
(1)降雨強度—歷時(ID)
該模型最早由Caine(1980)提出,可用下式表示:
I=C+∝×Dβ
(2)
式(2)中:I為誘發(fā)滑坡體等地質災害的降雨強度,短歷時取整個降雨事件的峰值數(shù)據(jù),長歷時則取平均值/(mm/h);D為誘發(fā)滑坡的降雨事件歷時/h;α、β均為統(tǒng)計參數(shù);C≥0。
該模型中的降雨歷時D取值范圍為1~100 h,全球各區(qū)域的ID模型閾值可通過意大利的水文地質研究保護所的降雨閾值庫進行查詢。國內外地區(qū)的ID模型參數(shù)統(tǒng)計見表1。
(2)年均降雨量閾值(MAP)或其他參數(shù)
不同區(qū)域之間引發(fā)滑坡的降雨閾值是不同的,一些學者則根據(jù)研究區(qū)的氣候特點采用當?shù)氐慕涤炅刻卣鲾?shù)據(jù)進行閾值計算,較為少用,偏經驗性與區(qū)域性。
(3)累積降雨量—歷時關系閾值(ED)
滑坡的發(fā)生在一定程度上與一個降雨事件的相關性較大,因此有些學者將累計雨量與降雨歷時進行分析,具體表達式如下:
E=C+∝×Dβ
(3)
式(3)中:E為累積降雨量/mm;D、α、β、C意義單位同上。
(4)累積降雨量—平均降雨強度(EI)
使用平均降雨強度I,式(3)中的降雨歷時D:
E=C+∝×Iβ
(4)
2.1.2 前期有效降雨量模型
降雨誘發(fā)滑坡、泥石流等地質災害是一個過程性問題,前期降雨量會直接影響到地下水水位和土壤含水量的變化,這也是誘發(fā)滑坡的重要因素。意大利東北部的滑坡研究表明,當前15天的累計降雨超過200 mm時,滑坡的發(fā)生與2天內降雨量是否大于200 mm有很大的相關性,當大于200 mm時,總會發(fā)生滑坡(Kohler and Linsley,1951),也就是說前期降雨量較大時,即便在較低的降雨強度下,滑坡也有可能發(fā)生,但前期降雨量會隨著時間增大而衰減,前期降雨量示意見圖7。如何界定有效的前期降雨量也是重要的研究問題。
Kohler and Linsley(1951)采用獲取的水文數(shù)據(jù)來測定地表徑流,確定降雨對徑流的影響;Crozier(1999)則將其引入降雨誘發(fā)滑坡預報研究,用于計算前期降雨量指數(shù)API(A ntecedent Precipitation Index),該指數(shù)計算公式如下:
Pa0=KR1+K2R2+…+KnRn,0 (5) 式(5)中:Pa0為相對0d的經校正的前期有效降雨量;R1為滑坡發(fā)生前一天的降雨量;Rn為滑坡發(fā)生前第n天的降雨量,K為土層中水流出的衰減系數(shù)。 許多研究者主要采用上式或者對其略做修改后作為計算前期有效降雨量的方法,根據(jù)研究區(qū)特點選擇合適的n和K。下表為一些研究學者的計算公式,現(xiàn)有各種前期有效降雨量計算方法見表2。 在公式(5)中,滑坡發(fā)生前的每天降雨衰減都獨立的,但事實中并不以天為單位遞減,每天的降雨衰減速度與其前一天的降雨事件相關,在此基礎上,李長江等(2011)學者建立了基于降雨—滑坡分形關系的前期有效降雨量模型: (6) 式(6)中:PEa為相對滑坡發(fā)生時(0 d)的前期有效降雨量;α為標度指數(shù);Pi為0 d前第i天的降雨量 筆者對國內以前期降雨量為參數(shù)的降雨臨界閾值確定進行了統(tǒng)計,國內一些地區(qū)的降雨量閾值研究方法見表3。 除了考慮降雨閾值與滑坡的關系外,降雨導致的土壤含水量和潛在蒸發(fā)量變化也會直接導致滑坡的發(fā)生,因此可以通過建立土壤含水量與滑坡是否發(fā)生之間的關系對滑坡進行預報(廖忠湞等,2019)。土壤含水量一方面可以通過布設監(jiān)測設備獲取,另一方面可以根據(jù)收集到的降雨、氣溫和潛在蒸發(fā)量與土壤屬性進行計算。該模型在已知降雨數(shù)據(jù)和土體含水量數(shù)據(jù)的條件下,計算滑坡發(fā)生的概率,采用土體含水量指數(shù)進行表示。土體水狀態(tài)指數(shù)有兩種表達形式:土體含水量不足量(DS<0)和日降雨過剩量(EP≥0)。 DS0=DS1-(P0-PE0) (7) 式(7)中:DS0是土體含水量不足量;DS1是第0天的前一天土體含水量不足量;P0是第0天的降雨量;PE0是第0天潛在蒸發(fā)量。當DS0≥0時,土體飽和,過剩雨量出現(xiàn),土體具有正的孔隙水壓力的特征,其作用是滑坡發(fā)生的前期因素。 日降雨過剩量EP的表達式為: EP0=(P0-PE0)-DS1 (8) 式(8)中:EP0是第0天的降雨過剩量。將經過衰減計算后的過剩降雨量進行累計后,就構成了正的土壤含水量指數(shù)。 EPa0=KEP1+K2EP2+…+KnEPn (9) 該模型可以通過氣象站提供的日蒸發(fā)量和降雨量對土壤含水量進行校正,結合上式計算出每天的前期土壤含水量,根據(jù)歷時滑坡數(shù)據(jù)建立的回歸方程計算出下一天達到滑坡條件下的土壤含水量閾值的所需的日降雨量,然后進行預報(譚萬沛和韓慶玉,1992)。 前兩個模型均需要收集研究區(qū)歷時滑坡發(fā)生時間及對應的降雨數(shù)據(jù)才能開展閾值計算,對于歷史數(shù)據(jù)不充分的區(qū)域便無法運用。水文學—斜坡穩(wěn)定性則是通過物理模型,考慮滑坡體本身的力學特性、土層厚度等特征,將研究區(qū)按評價單元網格進行劃分,結合DEM數(shù)據(jù)提取坡度等地形數(shù)據(jù),對斜坡體的穩(wěn)定性進行計算,從而確定失穩(wěn)時對應的臨界降雨量。該模型所需的物理參數(shù)較多,因空間差異使得較難大范圍地應用,其表達式如下: (10) 式(10)中:k為滲透系數(shù),由斜坡的巖土性質決定;t為潛在滑體厚度/m;b為斜坡體的過水寬度;θ為潛在滑面傾角,近似取斜坡坡度;A為區(qū)域柵格匯水面積;γ為斜坡坡體物質的重度;γw為地下水的重度;c為有效內聚力;φ為有效內摩擦角。 國內學者崔鵬等(2011)以汶川為研究區(qū),在收集其地形數(shù)據(jù)、巖土體物理屬性的基礎上,采用該模型對其降雨閾值進行計算,研究滑坡啟動的臨界降雨條件,并統(tǒng)計不同海拔、坡度內的閾值柵格數(shù)量。國外學者(Wilkinson et al.,2002)在歷史數(shù)據(jù)不充分的情況下,將坡地水文模型與斜坡穩(wěn)定性相結合,并結合GIS技術,實現(xiàn)區(qū)域的滑坡降雨閾值區(qū)劃,對降雨型滑坡災害預警有重要借鑒意義。 張洪奎(2020)指出,模型研究與監(jiān)測數(shù)據(jù)都是為了最終的預報預警。目前預報預警的研究區(qū)域多為地質地貌、氣象水文環(huán)境相對一致的區(qū)域,而長輸管道線路的研究較少。由于長輸管道沿線所跨區(qū)域的地質環(huán)境及氣候條件差異很大,很難以一個通用的閾值模型對地質災害進行預警預報,因此在管道沿線地質災害預警工作中,必須結合管道沿線的地質環(huán)境、氣象水文進行分區(qū)域的模型研究。 冼國棟等(2019)對蘭成原油管道途經區(qū)域劃分為黃土高原區(qū)、秦嶺大巴山高中山區(qū)和四川盆地區(qū)三個地貌單元,分別構建相應的評價指標體系進行地質災害易損性評估,利用GIS因子疊加技術,采用綜合模型將三段線路的危險性劃分為五個等級,但其研究成果偏向于空間上評估預測,未能從時間上及時對地質災害進行預警預報。也有學者在對管道沿線區(qū)域地質災害易發(fā)性評估的基礎上,根據(jù)實時的氣象信息計算前期有效降雨量,當降雨數(shù)據(jù)達到設定好的閾值時便發(fā)出預警信息(白利平等,2007),該監(jiān)測預警工作首先采用遙感技術對全線地質災害進行普查和識別,結合地質災害歷史數(shù)據(jù)進行管道沿線區(qū)域的易發(fā)性評估,分重點地布設北斗微位移監(jiān)測設備,并且采用前期有效降雨量模型和臨界降雨閾值設定了四個預警等級,對管道類地質災害監(jiān)測起到一定的借鑒參考意義,但其臨界降雨閾值是通用的,并未結合其易發(fā)性分區(qū)的結果進行設定。徐輝和海知(2019)通過北斗微位移的監(jiān)測數(shù)據(jù),與降雨數(shù)據(jù)進行對比,建立了北斗系統(tǒng)與氣象預報融合的監(jiān)測預警系統(tǒng),可以很好的在線路工程中應用。 如何將空間上的評估結果和降雨閾值結合起來進行預報是管道地質災害監(jiān)測預警的研究重點和難點,加之管道一般遠離居住區(qū),其地質災害歷史數(shù)據(jù)在很多情況是缺失的,無法構成降雨與滑坡發(fā)生之間的統(tǒng)計關系,這也說明通過統(tǒng)計分析的方法進行臨界降雨閾值的計算在多數(shù)情況下是不可行的。在上述的預報模型研究介紹中,水文學—斜坡穩(wěn)定性模型可能是管道沿線地質災害監(jiān)測一個比較好的解決方案,該模型通過計算斜坡穩(wěn)定性來確定臨界降雨量閾值,計算公式(10)中滑體厚度t一般需要現(xiàn)場踏勘測量潛在滑坡點進而獲得,其他參數(shù)可通過查閱《工程地質手冊》及利用ArcGIS軟件分析獲得,且其臨界降雨量的計算結果具有空間性,并非整個區(qū)域的通用閾值。綜上,長輸管道等行業(yè)地質災害監(jiān)測可首先采用遙感技術(SAR、高分影像)進行潛在災害點的識別,若歷史災害數(shù)據(jù)充分,可采用降雨閾值的方法計算臨界降雨量,對沿線地質災害進行實時監(jiān)測及預報,若歷史災害數(shù)據(jù)不充分,可對現(xiàn)場進行調查并采集潛在滑坡體厚度數(shù)據(jù),結合斜坡不穩(wěn)定模型計算降雨閾值,從而進行監(jiān)測預報。 降雨與滑坡發(fā)生之間存在著非常緊密的關系,有關的研究也非常之多,但其發(fā)生啟動機理以及預警預報工作仍有待進一步的研究與發(fā)展。首先應立足于滑坡本體的研究,聚焦于坡體結構、地層巖性、水文條件等基礎地質條件,探究滲流、水壓力、抗剪強度等內部參數(shù)對滑坡啟動之間的關系,總結降雨型滑坡的啟動機制理論研究。同時著重于工程應用,對于降雨等氣象條件在滑坡預警預報中的規(guī)律進行研究,總結不同地區(qū)、不同地質條件下降雨與滑坡發(fā)生之間的關系??偨Y降雨與滑坡之間的預警預報模型,提出降雨強度的計算方法以及預警閾值等參數(shù)。 目前利用GPS、GIS、北斗定位等技術,聯(lián)合地面物聯(lián)網設備對滑坡體進行實時監(jiān)測,通過對其巖土體形變位移的研究,評估其所處的狀態(tài)以及預測未來的發(fā)展趨勢?;诨掳l(fā)生發(fā)展的一般規(guī)律,以物聯(lián)網監(jiān)測數(shù)據(jù)作為其發(fā)展階段的標志。并結合降雨模型及閾值的研究,建立二者互相驗證、互相校核的預警預報體系。近些年出現(xiàn)的InSAR監(jiān)測技術已經越來越多的應用于災害監(jiān)測和預警當中,其高精度、高分辨率、全天候、范圍廣等特點,在災害和地表形變監(jiān)測中逐步成為一個新興方向,其良好的應用前景可結合常規(guī)的降雨、含水率、物聯(lián)網等參數(shù)組成一個可靠的驗證體系,搭建一個基于上述多手段、多數(shù)據(jù)的災害監(jiān)測預警系統(tǒng),為工程安全和人民生命安全服務。2.2 前期土壤含水量模型(ASWM)
2.3 水文學—斜坡穩(wěn)定性模型
3 長輸管道地質災害監(jiān)測預警應用
4 總結與展望