彭程輝 王常明 王鋼城

0 引言

冰湖是一種在常年積雪的冰川上發(fā)育的蓄水湖泊。近年來(lái)，隨著全球性氣候變暖加劇[1]，冰川融化速度加快，冰湖的數(shù)量和面積都有增大的趨勢(shì)。發(fā)生潰決災(zāi)害的冰湖類(lèi)型主要為兩種:冰川阻塞湖和冰磧阻塞湖(冰磧湖)[2]。據(jù)統(tǒng)計(jì)，喜馬拉雅山地區(qū)近50年來(lái)至少發(fā)生過(guò)20余次較大的冰湖潰決災(zāi)害[3]。發(fā)生災(zāi)害的冰湖類(lèi)型中，冰磧湖幾乎占了全部。以下為了方便，將冰磧湖簡(jiǎn)稱(chēng)為冰湖。冰湖潰決的機(jī)制一般分成兩類(lèi):1)終磧壟內(nèi)部埋藏的冰核融化引起終磧堤管涌破壞;2)冰川崩塌或者滑坡導(dǎo)致湖水位上漲并引起涌浪，造成終磧堤溢流型潰決[4]。這兩種機(jī)制中，后者發(fā)生的次數(shù)多、災(zāi)害大，是研究的重點(diǎn)。對(duì)于冰湖潰決機(jī)制的研究，很多學(xué)者[4-8]從氣候條件、水文條件和數(shù)值模擬等做了大量的研究，但是冰川滑坡的研究相對(duì)較少。劉晶晶[5]分析了冰舌的穩(wěn)定性，但是沒(méi)有考慮冰川融水和后緣冰川對(duì)冰舌穩(wěn)定性的影響。本文重點(diǎn)研究了冰川滑動(dòng)的力學(xué)條件，并采用潰決危險(xiǎn)性指數(shù) Idl[9]判斷冰湖穩(wěn)定性。

1 冰川滑坡的力學(xué)分析

由于長(zhǎng)期的風(fēng)化作用和應(yīng)力釋放，冰川前緣拉張裂縫密集發(fā)育。裂縫上部寬度可達(dá)0.5 m，下部寬度小于0.2 m，深入冰川內(nèi)部數(shù)米至十余米[10]，使冰川被切割成不連續(xù)的塊體。冰川前緣厚度較薄，容易發(fā)生滑動(dòng)和崩塌，冰崩體可達(dá)數(shù)百萬(wàn)立方米，相當(dāng)于冰湖蓄水容積的1/4～1/2[3]。下面對(duì)冰川前緣進(jìn)行力學(xué)分析。

1.1 冰川前緣滑動(dòng)力學(xué)模型

冰川與下覆基巖層的接觸面為滑動(dòng)面，平面滑動(dòng)。冰川前緣劃分為若干個(gè)相鄰的冰體來(lái)研究，文中取單位厚度的冰川進(jìn)行研究。沿冰川潛在滑動(dòng)方向?qū)Ρw編號(hào)，各冰體對(duì)應(yīng)的長(zhǎng)度為li(i為冰體編號(hào))，冰川總長(zhǎng)度為L(zhǎng);冰川厚度為h;冰川坡角為θ。冰川前緣滑動(dòng)力學(xué)模型如圖1所示。

圖1 冰川前緣滑動(dòng)力學(xué)模型圖

1.2 冰體1受力穩(wěn)定性分析

冰體1位于模型的最上緣，裂縫發(fā)育密集，裂縫寬度較大，延伸到冰川深部，將冰體1與其上緣冰川分割開(kāi)。冰體1受力條件如圖2所示。冰體1受力情況如下:自身重力G，基巖支持力N和抗滑力τl1。由于冰川融水沿著裂縫通道滲入到冰川內(nèi)部，對(duì)冰川滑動(dòng)有一定的潤(rùn)滑作用，因此需要對(duì)抗剪強(qiáng)度τ進(jìn)行折減。設(shè)折減系數(shù)為m且0＜m＜1，則抗剪強(qiáng)度τ為:

其中，l1為冰體1的長(zhǎng)度，m;θ為冰川的坡度，(°);c為內(nèi)聚力，kPa;φ 為內(nèi)摩擦角，(°)。下滑力 T滑=Gsinθ，抗滑力 T抗= τl1。冰體1穩(wěn)定平衡條件為:

將式(1)代入式(2)中，冰體1臨界平衡條件為Gsinθ=τl1即ρghl1sinθ=m(c+ ρghl1cosθ/l1tanφ)l1。其中，ρ為冰川的密度，kg/m3，得臨界冰川厚度h2:

圖2 冰體1受力示意圖

1.3 冰舌受力穩(wěn)定性分析

由于冰體不連續(xù)，不考慮下鄰冰體的作用，則冰體2～n受力條件相同，因此僅對(duì)冰舌進(jìn)行穩(wěn)定性分析。冰舌受力條件如圖3所示?？够?T抗= τln，下滑力 T滑=Pn－1+Gsinθ。其中，Pn－1為冰體n－1在冰舌上的作用力，分兩種情況討論:1)Pn－1=0時(shí)，冰體n－1穩(wěn)定，冰舌不受上部冰體的作用力;2)Pn－1＞0時(shí)，冰體n－1滑動(dòng)，對(duì)冰舌有向下的推動(dòng)力。此時(shí)臨界冰川厚度hn為:

由式(4)可見(jiàn)，當(dāng) Pn－1=0 時(shí)，冰川厚度和式(3)相同;當(dāng) Pn－1＞0時(shí)，臨界冰川厚度減小。

最不利的情況下，冰舌上緣整體滑動(dòng)，此時(shí)下滑力T滑=(G1+G2+…+Gn)sinθ=ρghLsinθ，抗滑力不變。此時(shí)臨界冰川厚度hn'為:

計(jì)算結(jié)果hn'＜hn＜h1，說(shuō)明冰川后緣滑動(dòng)，需要較大的厚度，才能克服抗剪強(qiáng)度;在冰舌部位，冰川厚度小，容易發(fā)生滑動(dòng)。

2 冰湖潰決危險(xiǎn)性評(píng)價(jià)

2.1 臨界冰川厚度計(jì)算

取臨界冰川厚度最大值h1作為冰川滑動(dòng)的判別依據(jù)。當(dāng)冰湖后緣冰川平均厚度大于h1時(shí)，冰川發(fā)生滑動(dòng)。為了方便計(jì)算，對(duì)式(3)進(jìn)行簡(jiǎn)化:取折減系數(shù)m=1，冰川與基巖接觸面上的內(nèi)聚力 c=100 kPa[11]，內(nèi)摩擦角 φ =0°。則式(3)化簡(jiǎn)為:

其中，ρ為冰川密度，取900 kg/m3;g為重力加速度，取10 m/s2。阿亞錯(cuò)等五個(gè)冰湖野外調(diào)查統(tǒng)計(jì)資料如表1所示。

圖3 冰舌受力示意圖

表1 冰湖野外調(diào)查統(tǒng)計(jì)資料

2.2 冰湖潰決危險(xiǎn)性指數(shù)計(jì)算

冰湖潰決危險(xiǎn)性指數(shù)Idl，即危險(xiǎn)冰川的體積與冰湖蓄水體積比值的倒數(shù)1/R[9]，R值越大冰湖潰決的可能性越小[12]。研究區(qū)冰川類(lèi)型屬于亞大陸型冰川[13]，因此危險(xiǎn)冰川的平均厚度H計(jì)算公式如下[14]:

其中，H為危險(xiǎn)冰川平均厚度，m;F為危險(xiǎn)冰川面積，km2。冰湖潰決危險(xiǎn)性指數(shù)計(jì)算結(jié)果如表2所示。

表2 冰湖潰決危險(xiǎn)性指數(shù)計(jì)算結(jié)果

2.3 冰湖潰決危險(xiǎn)性分析

阿亞錯(cuò)危險(xiǎn)冰川平均厚度39.2 m遠(yuǎn)小于臨界冰川厚度68 m，冰川較穩(wěn)定，不會(huì)發(fā)生滑動(dòng)。桑旺錯(cuò)危險(xiǎn)冰川平均厚度71.8 m和臨界冰川厚度73.5 m，兩者相差很小，冰川滑動(dòng)的可能性大。野外調(diào)查得知，該冰湖曾于1954年潰決過(guò)一次，終磧壟厚度超過(guò)600 m，且湖水位與下游河道沒(méi)有集中落差，因此冰湖的穩(wěn)定性較好。郭洛錯(cuò)危險(xiǎn)冰川平均厚度57 m遠(yuǎn)大于臨界冰川厚度27.8 m，冰川容易發(fā)生滑動(dòng)，冰湖不安全。郭洛強(qiáng)錯(cuò)臨界冰川厚度37 m小于危險(xiǎn)冰川平均厚度58.7 m，冰川可能發(fā)生滑動(dòng)。印達(dá)普錯(cuò)臨界冰川厚度74 m大于危險(xiǎn)冰川平均厚度49.9 m，因此冰川不大可能發(fā)生滑動(dòng)。潰決危險(xiǎn)性指數(shù)0.26較小，說(shuō)明冰湖穩(wěn)定，發(fā)生潰決可能性小。但現(xiàn)場(chǎng)調(diào)查顯示，終磧堤容易產(chǎn)生滲透變形破壞，引起冰湖潰決。

最終評(píng)價(jià)結(jié)果:阿亞錯(cuò)和桑旺錯(cuò)穩(wěn)定性較高，發(fā)生潰決的可能性較小;郭洛錯(cuò)和郭洛強(qiáng)錯(cuò)穩(wěn)定性較低，發(fā)生潰決的可能性較大;印達(dá)普錯(cuò)穩(wěn)定性低，潰決可能性大。

3 結(jié)論和建議

冰湖地處高原地區(qū)，不可能做到實(shí)時(shí)監(jiān)測(cè)，而冰湖潰決造成的災(zāi)害往往是難以估量的，因此對(duì)冰湖潰決的危險(xiǎn)性評(píng)價(jià)顯得尤為重要。通過(guò)文中研究得到了以下結(jié)論:1)冰川滑動(dòng)的部位主要集中在冰川前緣，前緣裂縫發(fā)育。2)臨界冰川厚度可以作為冰川滑動(dòng)的判別依據(jù)，對(duì)冰川的穩(wěn)定性進(jìn)行定量分析。3)危險(xiǎn)性指數(shù)適用于冰川滑坡引起的冰湖潰決危險(xiǎn)性評(píng)價(jià)，而對(duì)于因終磧堤滲透變形破壞引起的冰湖潰決評(píng)價(jià)結(jié)果不甚理想。4)危險(xiǎn)冰川的范圍界定是未來(lái)研究的一個(gè)重要方向。

[1]Houghton J T，Ding Y，Griggs D J，et al.Climate Change 2000:The Scientific Basis[M].Cambridge，UK:Cambridge University Press，2001:1-880.

[2]徐道明，馮清華.冰川泥石流與冰湖潰決災(zāi)害研究[J].冰川凍土，1988，10(3):284-289.

[3]徐道明，馮清華.西藏喜馬拉雅山區(qū)危險(xiǎn)冰湖及其潰決特征[J].地理學(xué)報(bào)，1989，44(3):343-352.

[4]崔 鵬，馬東濤，陳寧生，等.冰湖潰決泥石流的形成、演化與減災(zāi)對(duì)策[J].第四紀(jì)研究，2003，23(6):621-628.

[5]劉晶晶，程尊蘭，李 泳，等.西藏終磧湖潰決形式研究[J].地學(xué)前緣(中國(guó)地質(zhì)大學(xué)(北京)北京大學(xué))，2009，16(4):372-380.

[6]蔣忠信，崔 鵬，蔣良濰.冰磧湖漫溢型潰決臨界水文條件[J].鐵道工程學(xué)報(bào)，2004，4(84):21-26.

[7]陳儲(chǔ)軍，劉 明，張 幟.西藏年楚河冰川終磧湖潰決條件及洪水估算[J].冰川凍土，1996，18(4):347-352.

[8]岳志遠(yuǎn)，曹志先，車(chē) 濤，等.冰湖潰決洪水的二維水動(dòng)力學(xué)數(shù)值模擬[J].冰川凍土，2007，29(5):756-763.

[9]呂儒仁，唐邦興，朱平一，等.西藏泥石流與環(huán)境[M].成都:四川科技大學(xué)出版社，1999:69-105.

[10]李德基，游 勇.西藏波密米堆冰湖潰決淺議[J].山地研究，1992，10(4):219-224.

[11]施雅風(fēng)，黃茂桓，任炳輝，等.中國(guó)冰川概論[M].北京:科學(xué)出版社，1988:92.

[12]陳曉青，崔 鵬，楊 忠，等.喜馬拉雅山中段波曲流域近期冰湖潰決危險(xiǎn)性分析與評(píng)估[J].冰川凍土，2007，29(4):509-516.

[13]施雅風(fēng)，黃茂桓，姚檀棟，等.中國(guó)冰川與環(huán)境——現(xiàn)在、過(guò)去和未來(lái)[M].北京:科學(xué)出版社，2000:26-27.

[14]施雅風(fēng)，劉潮海，王宗太，等.簡(jiǎn)明中國(guó)冰川目錄[M].上海:上?？茖W(xué)普及出版社，2005:15.