李林澤, 常 鳴, 李宏杰, 胡錦鳳, 徐恒志, 劉沛源
(1. 成都理工大學(xué) 地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護(hù)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 四川 成都 610059;2. 山東正元數(shù)字城市建設(shè)有限公司, 山東 煙臺(tái) 264670)
我國(guó)是世界上震后地質(zhì)災(zāi)害活動(dòng)最頻繁、損失最嚴(yán)重的國(guó)家之一,其中泥石流災(zāi)害是震后山區(qū)地質(zhì)災(zāi)害中危險(xiǎn)性最大的災(zāi)害之一。強(qiáng)烈的地震作用在震后短期內(nèi)具有誘發(fā)大量山體滑坡并產(chǎn)生松散堆積物分布在溝道內(nèi)及溝道兩岸的能量,對(duì)山體造成的內(nèi)部損傷影響亦可達(dá)數(shù)百年之久[1-2]。與此同時(shí),地震也大大降低了誘發(fā)泥石流災(zāi)害的降雨閾值,導(dǎo)致災(zāi)難性的泥石流事件頻繁發(fā)生,規(guī)模加大,人民的生命財(cái)產(chǎn)安全受到極大地危害和威脅[3-4]。如2008年5月12日汶川特大地震后在地震的影響范圍內(nèi)接連暴發(fā)了大規(guī)模泥石流,該地震10年后雖然仍有大規(guī)模泥石流暴發(fā),但在震后短期內(nèi)泥石流啟動(dòng)降雨閾值普遍降低的情況下雨季暴發(fā)泥石流的頻次更高,規(guī)模更大,受災(zāi)更嚴(yán)重。因此,對(duì)震后短期內(nèi)可能因強(qiáng)降雨引發(fā)泥石流的運(yùn)動(dòng)過(guò)程和沖出規(guī)模需要更為深入的研究。
近二十年來(lái),國(guó)內(nèi)外大量學(xué)者對(duì)泥石流運(yùn)動(dòng)過(guò)程及沖出規(guī)模開(kāi)展的研究可分為三類(lèi):基于體積方量的預(yù)測(cè)模型、動(dòng)態(tài)預(yù)測(cè)模型和基于地形參數(shù)的預(yù)測(cè)模型[5-9]。隨著科技水平的提高,AI、數(shù)值計(jì)算等方法不斷出現(xiàn)并更新迭代,其中泥石流數(shù)值計(jì)算是理論基礎(chǔ)和計(jì)算數(shù)學(xué)的結(jié)合,將現(xiàn)實(shí)中的泥石流運(yùn)動(dòng)全過(guò)程簡(jiǎn)化為離散方程后進(jìn)行求解,加之?dāng)?shù)值計(jì)算平臺(tái)搭建完善后使用成本低且操作簡(jiǎn)捷,近年來(lái)被大量運(yùn)用到泥石流研究中[10-13]。但如FLO-2D、FLO-3D、VolcFlow、PFC等數(shù)值計(jì)算模型均是以獲得泥石流流深及流速為目的的單一沖出模擬,不能很好地耦合降雨、啟動(dòng)、侵蝕、運(yùn)移、堆積整個(gè)動(dòng)態(tài)過(guò)程,并且忽略了泥石流運(yùn)動(dòng)過(guò)程中尤為重要的侵蝕沖刷。
OpenLISEM數(shù)值計(jì)算模型是一個(gè)集降雨輸入、徑流、啟動(dòng)、侵蝕、運(yùn)移、沖出、堆積于一體的綜合水文模型,其原理是耦合分布式流域水文模型與二相動(dòng)量守恒定律,特點(diǎn)是把流域劃分成若干網(wǎng)格區(qū)域,對(duì)每個(gè)網(wǎng)格區(qū)域分別輸入不同植被、土壤、巖土參數(shù)和高程等屬性,輸入實(shí)時(shí)降雨條件后分別計(jì)算產(chǎn)流量;通過(guò)比較相鄰網(wǎng)格的高程確定各網(wǎng)格的流向,演算降雨過(guò)程中泥石流的啟動(dòng)、運(yùn)移、堆積整個(gè)過(guò)程[14]。該模型不但考慮了地形、植被、物源和水文過(guò)程,同時(shí)也避免了以往常規(guī)模擬軟件需要的泥石流流量過(guò)程線(xiàn)計(jì)算與輸入、啟動(dòng)點(diǎn)位置選取等問(wèn)題的不準(zhǔn)確性,能更加真實(shí)地重現(xiàn)泥石流暴發(fā)時(shí)的場(chǎng)景。
受瀘定6.8級(jí)地震影響,四川省瀘定縣落井溝流域新增大量物源,未來(lái)雨季,在泥石流激發(fā)降雨閾值降低的條件下研究區(qū)極有可能暴發(fā)泥石流,威脅溝口聯(lián)合村隧道、在建的瀘定—石棉高速與溝口居民生命財(cái)產(chǎn)安全。本文以落井溝為研究區(qū),通過(guò)分析其物源分布特征,選取可以較好地模擬泥石流運(yùn)動(dòng)侵蝕全過(guò)程的OpenLISEM數(shù)值計(jì)算模型對(duì)不同降雨強(qiáng)度下的落井溝泥石流運(yùn)動(dòng)全過(guò)程進(jìn)行模擬預(yù)測(cè),為該區(qū)域的泥石流危險(xiǎn)評(píng)價(jià)、預(yù)警預(yù)報(bào)及工程治理提供具有一定適用性和參考性的方法。
落井溝位于四川省瀘定縣得妥鎮(zhèn)中部,溝口坐標(biāo)102°10′42″E,29°33′33″N,其形成區(qū)距2022年9月5日瀘定6.8級(jí)地震震中約1.2 km,位處Ⅸ級(jí)烈度區(qū)(地理位置如圖1)。流域較小,總面積約3.4 km2,地勢(shì)復(fù)雜,東高西低,地形起伏范圍為1 089~2 572 m,主溝長(zhǎng)約2.7 km,流向自東南向西北,平均縱坡降約為390‰,溝谷狹長(zhǎng),兩岸坡度較陡,為明顯的“V”字型,屬典型的高山峽谷區(qū)。流域內(nèi)發(fā)育有兩條小型支溝,于主溝中部交匯并最終流入與溝口正交且流向?yàn)闁|北至西南的大渡河。
圖1 研究區(qū)地理位置Fig.1 Geographical location of the study area
圖2 瀘定縣2022-09-05—2023-04-15地震概況Fig.2 General situation of earthquakes in Luding County during 2022-09-05—2023-04-15
圖3 落井溝坡度分級(jí)與統(tǒng)計(jì)Fig.3 Slope grading and statistics in Luojing Gully
流域年平均氣溫15.7 ℃,年最低氣溫-4.5 ℃,年平均降雨量664.4 mm。從山巔至谷底可觀(guān)察到明顯的氣候垂直分帶現(xiàn)象,海拔低處氣候較溫暖濕潤(rùn),高處較嚴(yán)寒干燥,屬典型的亞熱帶季風(fēng)氣候。研究區(qū)屬華南地層,形成區(qū)主要發(fā)育花崗巖、砂巖;流通區(qū)及堆積區(qū)主要發(fā)育二疊、三疊系變質(zhì)巖,以板巖為主,同時(shí)整個(gè)流域發(fā)育較多沉積巖,第四系堆積物雜亂且分布不均。得妥鄉(xiāng)位于川滇南北向構(gòu)造帶北段,為多組構(gòu)造的交叉部位,構(gòu)造環(huán)境復(fù)雜,新構(gòu)造運(yùn)動(dòng)強(qiáng)烈,主要斷裂為鮮水河斷裂東南段(也稱(chēng)磨西斷裂)及大渡河斷裂,且大渡河斷裂橫穿落井溝堆積區(qū)[15]。
2022年9月5日12時(shí)52分,四川省甘孜州瀘定縣發(fā)生6.8級(jí)地震,震源深度16 km,震中位于29.59°N,102.08°E,震源機(jī)制解以走滑型為主。對(duì)四川省地震局提供的震訊進(jìn)行統(tǒng)計(jì),發(fā)現(xiàn)從2022年9月5日6.8級(jí)地震截止到2023年4月15日,瀘定縣共發(fā)生32次地震,其中包括3次5級(jí)及以上地震,即2022年9月5日6.8級(jí)、2022年10月22日5.0級(jí)及2023年1月26日5.6級(jí)地震,為典型的“主震-余震型”地震。
經(jīng)現(xiàn)場(chǎng)調(diào)查,“9·5”瀘定地震直接造成落井溝內(nèi)產(chǎn)生大量滑坡,溝道內(nèi)堆積較多的松散堆積物,后續(xù)接連暴發(fā)的31次地震更使得溝道兩側(cè)山體內(nèi)部損傷程度及溝道堆積物松散程度逐步增加。泥石流物源在震后得到大量補(bǔ)給,一旦暴發(fā)將對(duì)溝口聯(lián)合村隧道及在建的瀘定—石棉高速的使用和建造造成巨大影響,泥石流若發(fā)生堵江也可能對(duì)溝口東南部約11 km處正在運(yùn)營(yíng)的大崗山水電站造成難以估量的損失。
陡峻的溝道地形條件是形成泥石流至關(guān)重要的一環(huán)[16]。落井溝流域發(fā)育有兩條小型支溝,主溝頂部至溝口平均縱坡降為390‰,地勢(shì)陡峭,為明顯的“V”字型,屬高山峽谷深度切割區(qū)域。其溝谷狹長(zhǎng),兩岸坡度較陡,統(tǒng)計(jì)發(fā)現(xiàn)近70%的流域面積分布在坡度為30°~50°范圍內(nèi),滑坡松散物質(zhì)運(yùn)移到溝道的地形條件極好。這樣的地形條件對(duì)溝道與坡面物源穩(wěn)定程度及降水形成徑流的水動(dòng)力強(qiáng)度均具有正向的交互效應(yīng),利于泥石流發(fā)生。
近年來(lái)對(duì)泥石流物源的識(shí)別技術(shù)愈加豐富,如通過(guò)遙感高精度影像進(jìn)行滑坡的人工解譯[17]與通過(guò)合成孔徑雷達(dá)干涉測(cè)量(Interferometric Synthetic Aperture Radar,ISAR)技術(shù)進(jìn)行潛在滑坡的隱患識(shí)別[18],這些識(shí)別技術(shù)無(wú)一例外都需通過(guò)現(xiàn)場(chǎng)復(fù)核以檢驗(yàn)其適用性及準(zhǔn)確性,而一些極為困難且危險(xiǎn)的野外調(diào)查及測(cè)繪工作難以正常進(jìn)行,因此無(wú)人機(jī)低空數(shù)字?jǐn)z影野外測(cè)量技術(shù)開(kāi)始廣泛應(yīng)用[19]。
瀘定縣6.8級(jí)地震后,通過(guò)遙感解譯技術(shù)獲取落井溝流域新增滑坡37處,面積約0.125 km2。震后新增滑坡位置分布見(jiàn)圖4(b)。
圖4 滑坡震前震后影像對(duì)比Fig.4 Comparison between landslide images before and after the earthquake
震后不久對(duì)其進(jìn)行現(xiàn)場(chǎng)調(diào)查并使用搭載長(zhǎng)測(cè)程激光雷達(dá)的無(wú)人機(jī)進(jìn)行仿地飛行采集數(shù)據(jù),通過(guò)大量點(diǎn)云數(shù)據(jù)處理等方式獲得了分辨率為0.2 m的高精度機(jī)載激光雷達(dá)(Light Detection and Ranging,LiDAR)數(shù)據(jù),包括數(shù)字高程模型(Digital Elevation Model,DEL)及數(shù)字正射影像(Digital Orthophoto,DO)數(shù)據(jù),用于輔助識(shí)別落井溝流域內(nèi)因植被覆蓋未能識(shí)別的崩滑體,并作為基礎(chǔ)數(shù)據(jù)進(jìn)行后期泥石流運(yùn)動(dòng)全過(guò)程數(shù)值計(jì)算。通過(guò)圖5中b1與b2對(duì)比,發(fā)現(xiàn)去植被后的LiDAR數(shù)字高程模型能識(shí)別出影像不易察覺(jué)的滑坡堆積物,同時(shí),由于較多滑坡碎屑流僅是淺層物質(zhì)剝落,故通過(guò)LiDAR航測(cè)處理得到的山體陰影數(shù)據(jù)進(jìn)行精準(zhǔn)識(shí)別,剔除未產(chǎn)生與產(chǎn)生極少物源的滑坡,得到實(shí)際滑坡堆積面積約7.798 萬(wàn)m2(初始滑坡分布見(jiàn)圖1,處理后的滑坡分布見(jiàn)圖5)。
圖5 LiDAR數(shù)字正射影像及去植被的數(shù)字高程模型(圖a1與b1分別是落井溝兩處去植被后 的溝道與坡面滑坡LiDAR數(shù)據(jù);圖a2與b2分別是落井溝兩處未去植被的LiDAR數(shù)字 正射影像數(shù)據(jù);c為滑坡取樣點(diǎn))Fig.5 LiDAR digital orthophoto and digital elevation model for vegetation removal
結(jié)合滑坡體積-面積經(jīng)驗(yàn)公式[20][式(1)]得到震后落井溝溝道及坡面新增的實(shí)際滑坡松散堆積物方量為15.52 萬(wàn)m3,且主要集中在主溝中下段(如圖5)。溝道堆積物易被強(qiáng)降雨作用下形成的地表徑流沖刷侵蝕,溝岸坡腳堆積物被沖刷帶走后形成臨空面又會(huì)不斷失穩(wěn)直至運(yùn)移到溝道,增加溝道物源堆積,且隨海拔下降物源侵蝕會(huì)不斷增加,進(jìn)而暴發(fā)泥石流災(zāi)害。
v=1.315*A1.201 8
(1)
式中:A為同震滑坡面積(m2);v為同震滑坡體積(m3)。
震后大量崩滑體在強(qiáng)降雨作用下極易轉(zhuǎn)化為泥石流[21]。瀘定縣年平均降雨量664.4 mm,雨季降雨充沛,最高月降雨量可達(dá)167 mm,受“V”型峽谷地貌影響,該流域匯水形成徑流的能力較強(qiáng),加之大渡河從溝口流經(jīng),地下水豐富,流域內(nèi)主要為松散堆積層孔隙水和基巖裂隙水,補(bǔ)給方式通常為大氣降水。
泥石流災(zāi)害的發(fā)生是地形、物源及降雨等多因素共同作用的結(jié)果[22]。通過(guò)對(duì)落井溝地形、物源及降雨三個(gè)影響因素分析發(fā)現(xiàn),近70%的流域面積分布在坡度為30°~50°范圍內(nèi),“9·5”瀘定震后溝道及坡面的實(shí)際新增滑坡松散堆積物方量為15.52 萬(wàn)m3,且主要集中在主溝中下段,加之研究區(qū)雨季降雨充沛,地下水豐富,在未來(lái)的強(qiáng)降雨條件下具備暴發(fā)泥石流的初始條件,故對(duì)落井溝泥石流啟動(dòng)、物質(zhì)運(yùn)移及堆積規(guī)模借助數(shù)值計(jì)算工具進(jìn)行深入分析。
OpenLISEM數(shù)值計(jì)算模型通過(guò)建立柵格單元,融合SWATRE滲透模型、理查茲方程、圣維南二維非恒定流方程、二相動(dòng)量守恒定律及無(wú)限邊坡理論,實(shí)現(xiàn)災(zāi)害啟動(dòng)、侵蝕堆積和控制流體的微分方程求解,以此模擬地表水沙平衡,其具有很明顯的時(shí)空差異,是一個(gè)離散的二維數(shù)值模型。該模型可用于模擬流域在降雨期間或降雨后短時(shí)間內(nèi)的物源變化、土地利用類(lèi)型變化、防治工程等制約因素對(duì)泥石流的沖出趨勢(shì)、物源侵蝕、洪水等動(dòng)態(tài)流動(dòng)趨勢(shì)的影響。
模擬所需數(shù)據(jù)主要包括物源數(shù)據(jù)、地形數(shù)據(jù)和降雨數(shù)據(jù),其中地形數(shù)據(jù)利用QGIS平臺(tái)的轉(zhuǎn)換功能將0.2 m的機(jī)載LiDAR數(shù)字高程模型(DEM)轉(zhuǎn)換成模擬所需map格式文件,其他地形數(shù)據(jù)利用PCraster 軟件和Arcgis平臺(tái)制作成map格式文件,包括流向網(wǎng)、地勢(shì)梯度及表面粗糙度等。物源的物理力學(xué)特性指標(biāo)包括滲透系數(shù)、初始含水率、內(nèi)摩擦角、孔隙率、內(nèi)部黏聚力、平均粒徑及特征粒徑D50和D90等,需通過(guò)野外調(diào)查及現(xiàn)場(chǎng)取樣進(jìn)行室內(nèi)試驗(yàn)。降雨數(shù)據(jù)僅需輸入泥石流暴發(fā)期間的降雨過(guò)程,但由于震后落井溝尚未進(jìn)入雨季且未暴發(fā)泥石流,無(wú)法獲取降雨啟動(dòng)泥石流后運(yùn)動(dòng)過(guò)程中的實(shí)際降雨歷時(shí)數(shù)據(jù),故結(jié)合5%、2%、1%降雨強(qiáng)度得到落井溝洪水流量過(guò)程線(xiàn),輸入到OpenLISEM數(shù)值計(jì)算模型中進(jìn)行模擬。
(1) 物源數(shù)據(jù)
通過(guò)0.2 m分辨率的LiDAR航拍數(shù)據(jù)處理得到的落井溝溝道與坡面滑坡松散堆積物體積方量約15.52 萬(wàn)m3。此外,還需獲取物源的土體級(jí)配特征D50及D90進(jìn)行建模,通過(guò)篩除大于20 mm的現(xiàn)場(chǎng)樣品并進(jìn)行室內(nèi)顆粒分析試驗(yàn)最終得到土體級(jí)配曲線(xiàn),獲取了D50及D90分別為5 mm、11 mm(圖6)。取樣點(diǎn)在流域中的位置分布見(jiàn)圖5中c點(diǎn)。
圖6 落井溝滑坡堆積物級(jí)配曲線(xiàn)(<20 mm)Fig.6 Grading curve of landslide deposit in Luojing Gully (<20 mm)
通過(guò)不同土體的初始含水率、孔隙度及滲透系數(shù)等參考值[23]加上現(xiàn)場(chǎng)調(diào)查與室內(nèi)試驗(yàn)綜合得到如表1所列的物源指標(biāo)。
表1 物源特征指標(biāo)及取值
(2) 洪峰流量
OpenLISEM進(jìn)行泥石流數(shù)值計(jì)算的降雨數(shù)據(jù)不需要通過(guò)計(jì)算獲取泥石流流量過(guò)程,僅需提供實(shí)際降雨過(guò)程便能結(jié)合物源及地形參數(shù)進(jìn)行自動(dòng)化模擬,得出最符合實(shí)際的泥石流流量過(guò)程。但由于落井溝尚未暴發(fā)泥石流,故采用四川省中小流域暴雨洪水計(jì)算手冊(cè)進(jìn)行計(jì)算,獲得落井溝泥石流20年、50年、100年一遇暴雨重現(xiàn)期的洪水流量水文圖。采用以下公式(2)~(9)進(jìn)行計(jì)算:
(2)
(3)
(4)
(5)
m=0.221θ0.204
(6)
(7)
(8)
(9)
式中:QB為暴雨流量(m3/s);ψ為洪峰徑流系數(shù);i為最大平均暴雨強(qiáng)度(mm/h);F為流域面積(km2);s為暴雨雨力(mm/h);n為暴雨公式指數(shù);τ為流域匯水時(shí)間(h)。K6、K24為不同重現(xiàn)期的模比系數(shù),可從皮爾遜Ⅲ型曲線(xiàn)模比系數(shù)表中獲取;θ為流域特征參數(shù);m為匯流參數(shù);τ0為當(dāng)ψ=1時(shí)的流域匯流時(shí)間(h);μ為產(chǎn)流參數(shù)(mm/h)。相關(guān)取值列于表2,洪水流量過(guò)程線(xiàn)見(jiàn)圖7。
表2 清水洪峰流量相關(guān)參數(shù)
圖7 落井溝洪水流量過(guò)程線(xiàn)Fig.7 Flood flow hydrograph in Luojing Gully
模擬過(guò)程需提前設(shè)置四個(gè)基礎(chǔ)全英文命名的文件,包含地形和物源的map文件、模擬運(yùn)行的命令流run文件、降雨數(shù)據(jù)rain文件與存儲(chǔ)輸出結(jié)果的res文件,設(shè)置完畢后對(duì)落井溝泥石流基于OpenLISEM模型進(jìn)行運(yùn)動(dòng)全過(guò)程模擬。
結(jié)果顯示:5%、2%、1%降雨強(qiáng)度下,落井溝泥石流運(yùn)動(dòng)過(guò)程的最大泥深分別為6、8、11 m,最大流速分別為5、7、10 m/s,整個(gè)運(yùn)動(dòng)過(guò)程的泥沙堆積顯示溝道內(nèi)淤積較多泥石流固體物質(zhì),推測(cè)受落井溝流域地形影響,于拐彎處發(fā)生彎道超高,降低了泥石流流速,導(dǎo)致泥石流淤積溝道的固體物質(zhì)較多(圖8)。泥沙隨空間的變化關(guān)系呈現(xiàn)為逐步從主溝上端運(yùn)移至溝口(圖9)。
通過(guò)分析落井溝泥石流模擬的最終沖出規(guī)模發(fā)現(xiàn):5%、2%降雨頻率下落井溝泥石流雖沖出溝口但對(duì)大渡河的堵江程度較小,而1%降雨頻率下落井溝泥石流造成大渡河半堵江狀態(tài)。
通過(guò)遙感解譯、現(xiàn)場(chǎng)調(diào)查及機(jī)載LiDAR航拍,對(duì)2022年9月5日瀘定6.8級(jí)地震后落井溝泥石流新增滑坡的發(fā)育特征進(jìn)行分析,結(jié)合OpenLISEM數(shù)值計(jì)算模型對(duì)落井溝泥石流在5%、2%、1%降雨強(qiáng)度下的運(yùn)動(dòng)全過(guò)程進(jìn)行模擬,得到運(yùn)動(dòng)過(guò)程中的最大堆積深度、最大流動(dòng)速度及最終沖出范圍,并給出1%降雨強(qiáng)度下泥石流的泥沙演化全過(guò)程,綜合反映出落井溝泥石流的運(yùn)動(dòng)侵蝕堆積情況,并對(duì)大渡河的堵塞程度給出建議。主要結(jié)論或建議如下:
(1) 受瀘定6.8級(jí)地震直接影響,流域新增滑坡37處,滑坡面積約0.125 km2。通過(guò)現(xiàn)場(chǎng)調(diào)查及0.2 m高精度分辨率機(jī)載LiDAR數(shù)據(jù)對(duì)滑坡進(jìn)行優(yōu)化識(shí)別,并結(jié)合滑坡體積-面積經(jīng)驗(yàn)公式得到落井溝溝道及坡面的實(shí)際滑坡松散堆積物方量為15.52萬(wàn)m3,且主要集中在主溝中下段。
(2) 基于OpenLISEM數(shù)值模擬平臺(tái),結(jié)合精細(xì)化地形、精準(zhǔn)物源及5%、2%、1%降雨強(qiáng)度對(duì)落井溝泥石流進(jìn)行運(yùn)動(dòng)全過(guò)程模擬,獲得運(yùn)動(dòng)過(guò)程最大泥深分別為6、8、11 m,最大流速分別為5、7、10 m/s。通過(guò)對(duì)溝口最終的堆積狀況分析,5%、2%降雨頻率下落井溝泥石流雖沖出溝口但堵江程度極小,1%降雨頻率下落井溝泥石流造成大渡河半堵江狀態(tài)。
(3) 落井溝泥石流溝域面積較小,即便在雨季暴發(fā)對(duì)溝口聯(lián)合村隧道和在建的瀘定—石棉高速的使用和建造在大渡河下游庫(kù)岸相連的大崗山水電站的正常運(yùn)營(yíng)及沿線(xiàn)居民的生命財(cái)產(chǎn)安全造成的威脅仍然較小。但值得注意的是,自2022年9月5日瀘定地震以來(lái)共發(fā)生大小地震32起,包括3次5級(jí)以上地震,這種小而多的地震對(duì)溝道及坡面滑坡物源的各方面影響(尤其是潛在滑坡臨界失穩(wěn)狀態(tài))有著肉眼不易察覺(jué)的特征,雨季豐沛的降雨及形成的徑流對(duì)該流域的潛在物源也有著較大影響,對(duì)瀘定震區(qū)此類(lèi)物源豐富、地勢(shì)陡峻的溝道仍需密切關(guān)注。