滑坡蠕變時(shí)效穩(wěn)定性理論與滑坡預(yù)測(cè)(報(bào))研究的基本構(gòu)想

楊人光

(中國(guó)四海控股有限公司生態(tài)災(zāi)害測(cè)控中心,北京 100081)

滑坡蠕變時(shí)效穩(wěn)定性理論與滑坡預(yù)測(cè)(報(bào))研究的基本構(gòu)想

楊人光

(中國(guó)四海控股有限公司生態(tài)災(zāi)害測(cè)控中心,北京 100081)

滑坡的形成過(guò)程顯示出其時(shí)效特征,已有研究顯示滑坡的時(shí)效特征有助于滑坡的精確預(yù)測(cè)預(yù)報(bào)研究。本文在馬克斯韋爾、開(kāi)爾文、齋滕迪孝等人的研究基礎(chǔ)上引入了巖土體粘-彈性滯后系數(shù),以巖土材料試件在各種受力狀態(tài)下的實(shí)際監(jiān)測(cè)軌跡曲線(xiàn)為依據(jù),采用擬合法,遵循曲線(xiàn)拐點(diǎn)形變量和形變率連續(xù)性導(dǎo)出關(guān)于時(shí)間的二元非線(xiàn)性超越代數(shù)方程組,并求得時(shí)效共軛點(diǎn)時(shí)間和劇滑失穩(wěn)時(shí)間,最后得到滑坡蠕變時(shí)效可預(yù)測(cè)性的結(jié)論。

滑坡;蠕變;粘-彈性滯后系數(shù);時(shí)效共軛點(diǎn);時(shí)效穩(wěn)定性

1 前言

滑坡是斜坡巖土體沿著貫通的剪切破壞面所發(fā)生的滑移現(xiàn)象,是地表的地質(zhì)災(zāi)害之一。國(guó)內(nèi)外眾多學(xué)者與工程技術(shù)人員進(jìn)行了大量室內(nèi)外監(jiān)測(cè)試驗(yàn)和探索性研究,積累大量資料,并在滑坡穩(wěn)定性與監(jiān)測(cè)預(yù)報(bào)的研究中形成了許多科學(xué)的概念和思路。

近半個(gè)世紀(jì)以來(lái),就滑坡體穩(wěn)定性研究,多數(shù)的國(guó)內(nèi)外學(xué)者及工程技術(shù)人員遵循經(jīng)典力學(xué)思維框架,從滑坡體受力角度或地質(zhì)力學(xué)角度,應(yīng)用巖土力學(xué)中的有限元、離散元等數(shù)值分析法,試圖揭示出非均質(zhì)各向異性,非規(guī)則分布結(jié)構(gòu)面的巖土體的應(yīng)力應(yīng)變狀態(tài),并以微元體主應(yīng)力分布積分及強(qiáng)度指標(biāo),再結(jié)合極限平衡方法確定滑坡滑動(dòng)優(yōu)勢(shì)面,導(dǎo)出滑坡穩(wěn)定性系數(shù)Ksf。它無(wú)疑是滑坡穩(wěn)定性判據(jù)的重要依據(jù),但缺乏時(shí)間的內(nèi)涵,作為滑坡的判據(jù)準(zhǔn)則,就有商榷的余地,正如崔政權(quán)教授(系統(tǒng)工程地質(zhì)導(dǎo)論著者)所述在確定論分析中,指標(biāo)Ksf≥1,表示邊坡處于穩(wěn)定,或者說(shuō)是安全的,問(wèn)題是這一安全究竟有多大把握性的問(wèn)題。為此,又引用以概率論為基礎(chǔ)的置信度或可靠性分析等,這種以純靜態(tài)力學(xué)分析法去研究滑體蠕變運(yùn)動(dòng),類(lèi)同“靜不定”問(wèn)題當(dāng)作靜定處理豈能奏效?

滑坡的形成有其特有的過(guò)程,該過(guò)程極其復(fù)雜且受多方面的因素影響,如:局部地質(zhì)結(jié)構(gòu)破壞、地應(yīng)力釋放、水的浸泡和蝕化作用、周邊環(huán)境(坡腳沖刷、重復(fù)性微震動(dòng)等)以及氣候狀況影響等諸多因素。根據(jù)滑坡類(lèi)型,探求滑坡的穩(wěn)定性系數(shù),在求解過(guò)程中的分析方法,便涉及到諸多的物理力學(xué)性質(zhì)參數(shù),參數(shù)的正確選取必然要利用多種手段和方法,但“未能取得令人滿(mǎn)意的效果”。

總之,滑坡的穩(wěn)定性問(wèn)題,是巖土體結(jié)構(gòu)及模型材料的基礎(chǔ)性理論課題,并非僅依賴(lài)于傳統(tǒng)經(jīng)典力學(xué)構(gòu)架范疇。至于是否系與現(xiàn)代邊緣性學(xué)科的交叉呢?也不得而知。作者以為任何基礎(chǔ)性學(xué)科應(yīng)遵循著學(xué)科固有的抽象模型和運(yùn)動(dòng)規(guī)律,這就驅(qū)使作者去追逐新思路,重新確立新的理念。

2 研究現(xiàn)狀

可喜的是,國(guó)內(nèi)外有許多學(xué)者走出經(jīng)典力學(xué)的思維框架,從滑坡監(jiān)測(cè)預(yù)測(cè)入手,建立巖土體流變(包括蠕變)的物理數(shù)學(xué)模型,如馬克斯韋爾(C. Maxwell)體或沃伊特(Voigt)和開(kāi)爾文(Kelvin)體。尤其是馬克斯韋爾粘-彈性滯后物理模型和開(kāi)爾文的數(shù)學(xué)模式,兩者相結(jié)合即為粘-彈性滯后模型。在恒定載荷力作用下,作者把它改寫(xiě)為二階微分方程式,即:引進(jìn)了巖土體粘-彈性滯后系數(shù)E、η分別為巖土彈性模量和粘滯性系數(shù)。它盡可擬合巖土穩(wěn)定蠕變時(shí)段,然而,卻不能推廣于蠕變?nèi)^(guò)程,因?yàn)榭繂我粷u變線(xiàn)性方程是無(wú)法預(yù)測(cè)(報(bào))失穩(wěn)突變時(shí)間的。故它被擱淺了,因而諸多科技人員從臨滑實(shí)測(cè)蠕變軌跡曲線(xiàn),采用回歸多項(xiàng)式方程,梯度-正弦或者應(yīng)用流變方程等擬合外推,即擬合法,有人也簡(jiǎn)稱(chēng)為“灰色理論”的組成部分,實(shí)際上往往把滑坡體初始萌發(fā)的流變時(shí)段加以擬合外推,就帶有經(jīng)驗(yàn)性和揣測(cè)性,同樣是不能奏效。

但近半個(gè)世紀(jì)來(lái),諸多學(xué)者作出了大量在拉、壓、彎曲、扭轉(zhuǎn)的受力狀態(tài)下,施加不同載荷比的巖土試件的滯后形變軌跡曲線(xiàn)。并且指出“滑坡劇滑時(shí)間就是在該曲線(xiàn)上尋求時(shí)間增量Δt→0位移增量Δs→∞的那一點(diǎn)”,這猶如在灰暗荒原燃起一垛智慧的篝火。

特別是1965年日本齋滕迪孝等人在室內(nèi)沙土盤(pán)蠕變時(shí)效試驗(yàn)基礎(chǔ)上提出第(Ⅲ)階段(作者稱(chēng)之為失穩(wěn)蠕變)數(shù)學(xué)模式,即

式中,A、α為待定系數(shù)。

在此基礎(chǔ)上提出第(Ⅲ)蠕變階段滑坡時(shí)間預(yù)報(bào)的經(jīng)驗(yàn)公式,因限于試驗(yàn)條件,他們忽略了初期的穩(wěn)定時(shí)段或者寓于線(xiàn)性穩(wěn)定蠕變時(shí)段如圖1。導(dǎo)致其誤差幅值就相差十幾倍(即10±0.59)?？上藗冞B同齋滕迪孝先生提出失穩(wěn)時(shí)段的數(shù)學(xué)模式也缺乏應(yīng)有的關(guān)注。

圖1 土質(zhì)斜坡等速蠕變速率與達(dá)到破壞時(shí)間關(guān)系圖解Fig.1 Illustrate the relation between the constant rate creep and the time to failure of soil slope

這就是當(dāng)前巖土蠕變學(xué)理論,滑坡穩(wěn)定性以及預(yù)測(cè)預(yù)報(bào)理論模式研究、探索的基本成果和認(rèn)同狀況。

3 巖土體蠕變時(shí)效穩(wěn)定性理論與滑坡預(yù)測(cè)預(yù)報(bào)研究的基本構(gòu)想

作者以為馬克斯韋爾、開(kāi)爾文、齋滕迪孝等人所探索出上述基本成果是卓越的,是開(kāi)拓性的,是鑄造出一塊塊堆放著的構(gòu)筑蠕變學(xué)或者說(shuō)巖土結(jié)構(gòu)穩(wěn)定性理論大廈的基石。

3.1 已有研究成果的啟示

(1)齋滕迪孝先生以沙土盤(pán)室內(nèi)蠕變時(shí)效試驗(yàn)為依據(jù)歸納出二階微分方程式,作者從中受到極大的啟迪。因?yàn)樵撐⒎址匠虜?shù)學(xué)模式是非線(xiàn)性的,量化了滑坡體從漸變到突變的演化,并且不含有受力狀態(tài)的顯式項(xiàng)。故而該數(shù)學(xué)模式有望適合于巖土體在各種復(fù)合受力狀態(tài)下失穩(wěn)蠕變時(shí)效性,賦予它普遍性的特征。作者根據(jù)該方程的解,擬合蠕變后期實(shí)際監(jiān)測(cè)賦值點(diǎn),便預(yù)測(cè)(報(bào))出滑坡體短期或臨滑的劇滑時(shí)間tf,因而它不失為失穩(wěn)時(shí)段的可預(yù)測(cè)預(yù)報(bào)性。

(2)馬克斯韋爾、開(kāi)爾文的粘-彈性滯后體的數(shù)學(xué)模式,恒定載荷力對(duì)時(shí)間的導(dǎo)數(shù)后亦不含有載荷項(xiàng),作者把它改寫(xiě)為二階微分方程式。如式(1)所示,它描述前期穩(wěn)態(tài)時(shí)段形變時(shí)效的規(guī)律,引出一個(gè)尚待共認(rèn)的巖體材料的粘-彈性滯后體的物性參量,以它來(lái)度量巖土材料的時(shí)效性的強(qiáng)弱,它是巖土脆性材料所特有的,即呈現(xiàn)出粘-彈性滯后特性。為擬合實(shí)際監(jiān)測(cè)曲線(xiàn)求解,反饋出巖土材料參量ξ,這就是規(guī)律的可認(rèn)識(shí)性與和諧性。因?yàn)閷?duì)巖土體而言,嚴(yán)格來(lái)說(shuō),參量ξ完全相同的巖土塊體是不存在的,更何況參與滑坡體蠕變時(shí)效并非滑塊自身而是滑動(dòng)面的巖土介質(zhì)。

(3)參量ξ制約著蠕變?nèi)^(guò)程,所謂可預(yù)測(cè)性是指穩(wěn)定時(shí)段終端,即失穩(wěn)點(diǎn),與它貼近的失穩(wěn)伊始時(shí)段就預(yù)測(cè)出滑體劇滑時(shí)間tf,因而穩(wěn)定時(shí)段是可預(yù)測(cè)預(yù)報(bào)性不可或缺的組成部分。

3.2 基本研究構(gòu)想

作者的工作就是將那一塊塊堆放著的基石構(gòu)筑起一座大廈。首先要論證以下兩道關(guān)鍵性的問(wèn)題。

(1)原本孤立量化出的巖土材料蠕變時(shí)效形變兩個(gè)二階微分方程(1)、(2)式,作者把它分別確認(rèn)為穩(wěn)定和失穩(wěn)蠕變時(shí)段的數(shù)學(xué)模式,它們能否擬合在同一條蠕變時(shí)效監(jiān)測(cè)軌跡曲線(xiàn)上呢?其聯(lián)立求解能否相互耦合,是否存在一個(gè)唯一的共軛點(diǎn)呢?即在曲線(xiàn)連續(xù)條件下求解出軌跡曲線(xiàn)的拐點(diǎn),即時(shí)效共軛點(diǎn)時(shí)間tp。

(2)滑體一旦進(jìn)入穩(wěn)定轉(zhuǎn)向失穩(wěn)伊始時(shí)段賦值點(diǎn)(即緊貼著拐點(diǎn))能否就預(yù)測(cè)預(yù)報(bào)出其劇滑時(shí)間tf,即賦予微分方程組的可預(yù)測(cè)預(yù)報(bào)性呢?

為此,作者以巖土材料試件在各種受力狀態(tài)下的實(shí)際監(jiān)測(cè)軌跡曲線(xiàn)為依據(jù),采用(1)、(2)式微分方程解的擬合法,遵循曲線(xiàn)拐點(diǎn)形變量和形變率連續(xù)性導(dǎo)出關(guān)于以tp、tf為未知量的二元非線(xiàn)性超越代數(shù)方程組,論證其函數(shù)方程在tp、tf的時(shí)域范圍內(nèi)函數(shù)值發(fā)生變化,則函數(shù)有零解。并應(yīng)用多重二分法數(shù)值程序求解得出實(shí)數(shù)根tp、tf,最終才成功地鎖定(1)、(2)式微分方程聯(lián)立確認(rèn)為巖土蠕變時(shí)效微分方程組,即巖土材料蠕變學(xué)理論的基本方程,以揭示出巖土材料在恒定施加載荷作用下,隨時(shí)間推移而發(fā)生形變的時(shí)間效應(yīng),暫稱(chēng)為蠕變時(shí)效規(guī)律。作者抓住馬克斯韋爾、開(kāi)爾文粘-彈性時(shí)效性模式和齋滕迪孝非線(xiàn)性微分方程,才導(dǎo)出巖土蠕變的可預(yù)測(cè)性。

同時(shí)應(yīng)指出,在以上兩道關(guān)鍵性問(wèn)題論證時(shí),蠕變時(shí)效載體的規(guī)律是指巖土整體所呈現(xiàn)出宏觀規(guī)律。如同巖土材料試件是指試件整體的形變量或者整體位移量的宏觀規(guī)律。這是巖土體自身物理力學(xué)性質(zhì)參數(shù),以及各點(diǎn)位微元體形變的無(wú)規(guī)律性,卻呈現(xiàn)出整體宏觀的規(guī)律性。

至此,讀者將會(huì)提出一個(gè)深層次的問(wèn)題:(1)、(2)式二階微分方程均不含有載荷項(xiàng),如何反映出施加于試件上的載荷比β,或者它的倒數(shù)即滑坡穩(wěn)定性系數(shù)Ksf在巖土材料蠕變時(shí)效規(guī)律中內(nèi)在關(guān)系,或者說(shuō)能否賦予Ksf更深刻的內(nèi)涵。

為此作者在論證(1)、(2)方程組彼此耦合、并求解蠕變時(shí)效共軛點(diǎn)唯一性的同時(shí),就演繹出方程組各參量數(shù)和施加載荷比β獨(dú)立的相關(guān)函數(shù)關(guān)系的曲線(xiàn)族。并且發(fā)現(xiàn)時(shí)效穩(wěn)定指數(shù)α曲線(xiàn)族卻重合出一條曲線(xiàn)。α-β相關(guān)函數(shù)曲線(xiàn)與巖土試件材料物性ξ以及受力狀態(tài)(拉、壓、彎曲、剪切、扭轉(zhuǎn))無(wú)關(guān),作者稱(chēng)之謂α-β規(guī)律關(guān)系曲線(xiàn)。其實(shí)α為方程(2)式的無(wú)因次指數(shù),故之α-β規(guī)律關(guān)系與方程(2)式是相和諧的。因而從任何滑坡監(jiān)測(cè)出軌跡曲線(xiàn),演釋出α值,則可由α-β規(guī)律關(guān)系曲線(xiàn)圖反查出β值。它的倒數(shù)即為滑坡穩(wěn)定性系數(shù)Ksf。因此蠕變時(shí)效方程組才賦予Ksf的穩(wěn)定時(shí)間(年限)的內(nèi)涵。對(duì)于滑坡體而言,嚴(yán)格來(lái)說(shuō),唯有滑體監(jiān)測(cè)出軌跡曲線(xiàn)方可給出有可信度的Ksf。

滑坡穩(wěn)定性系數(shù)Ksf定義和清晰概念的確立,不難確認(rèn)所謂滑坡理論即為滑坡穩(wěn)定性理論。那么如何更確切賦予蠕變時(shí)效方程組(1)、(2)式以及時(shí)效共軛點(diǎn)物理力學(xué)性的含義呢?

作者應(yīng)用量綱分析法,從(1)、(2)二階微分方程式出發(fā),導(dǎo)出蠕變時(shí)效系統(tǒng)的無(wú)因次數(shù),其中唯時(shí)效共軛點(diǎn)(tp、yp)組合參數(shù)

當(dāng)ve=1,巖土結(jié)構(gòu)處于永久穩(wěn)定;

當(dāng)ve>1,巖土結(jié)構(gòu)從穩(wěn)定轉(zhuǎn)向失穩(wěn)蠕變。

故ve具有判斷巖土結(jié)構(gòu)從穩(wěn)定轉(zhuǎn)向失穩(wěn)的特征,故稱(chēng)ve為結(jié)構(gòu)失穩(wěn)判據(jù),則其共軛點(diǎn)即為失穩(wěn)點(diǎn)。于是巖土材料粘-彈性時(shí)效模型、穩(wěn)定時(shí)段方程式(1)和失穩(wěn)非線(xiàn)性時(shí)效方程式(2),它們構(gòu)成巖土蠕變時(shí)效穩(wěn)定性理論的基本方程組?；路€(wěn)定性系為巖土結(jié)構(gòu)穩(wěn)定性一個(gè)最典型的特例。其實(shí),凡在恒定重力場(chǎng)驅(qū)動(dòng)下地殼板塊、區(qū)域性地塊以及地表上下的建筑物和構(gòu)筑物等均涉及蠕變時(shí)效穩(wěn)定性問(wèn)題,并揭示出蠕變時(shí)效穩(wěn)定性一般性原理。

通過(guò)分別論述蠕變時(shí)效可預(yù)測(cè)性、低值載荷比β相關(guān)函數(shù)分析,導(dǎo)出失穩(wěn)時(shí)段“自擬合”規(guī)律與巖土抗形變極限強(qiáng)度時(shí)效性,蠕變時(shí)效方程組的積分導(dǎo)出蠕變能量累積、貯存和瞬息釋放機(jī)制。并應(yīng)用現(xiàn)代數(shù)值程序模擬試驗(yàn)分析法,探求區(qū)域性地塊超巨大能量參量數(shù)的內(nèi)在關(guān)系。

[1]胡海濤,項(xiàng)式均.關(guān)中西部滑坡的結(jié)構(gòu)、構(gòu)造特性及穩(wěn)定性分析[J].地質(zhì)學(xué)報(bào),1965,45(3).

[2]李四光.地質(zhì)力學(xué)[M].科學(xué)出版社,1976.

[3]谷德振.巖體工程地質(zhì)力學(xué)基礎(chǔ)[M].科學(xué)出版社,1976.

[4]西德G,哥德赫(編),張清,張彌(譯).有限元法在巖土力學(xué)中的應(yīng)用[M].中國(guó)鐵道出版社,1983.

[5]楊人光.爆破拋體加速運(yùn)動(dòng)的滲流模型[J].力學(xué)學(xué)報(bào),1984,16 (1).

[6]張倬元.工程地質(zhì)分析原理[M].地質(zhì)出版社,1985.

[7]C C維來(lái)洛夫.土力學(xué)的流變?cè)韀M].科學(xué)出版社,1988.

[8]王思敬,1987全國(guó)滑坡學(xué)術(shù)討論會(huì)滑坡論文集[C].四川科學(xué)技術(shù)出版社,1989.

[9]國(guó)際滑坡與巖土工程學(xué)術(shù)會(huì)議論文集[C].華中理工大學(xué)出版社,1991.

[10]崔政權(quán).系統(tǒng)工程地質(zhì)導(dǎo)論[M].水利電力出版社,1992.

[11]張倬元.國(guó)際滑坡研究進(jìn)展——第六屆國(guó)際滑坡學(xué)術(shù)討論會(huì)述評(píng)[J].中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào),1993,4(3).

[12]晏同珍,等.滑坡文集第十集[C].中國(guó)鐵道出版社,1993.

[13]楊人光,孫善平.人工爆破治理火山災(zāi)害可行性探討[M].地震出版社,1994.

[14]文寶萍,李緩.黃土地區(qū)典型滑坡預(yù)測(cè)預(yù)報(bào)及減災(zāi)對(duì)策研究[M].地質(zhì)出版社,1997.

[15]鄭穎人.邊坡與滑坡工程治理[M].人民交通出版社,2007.

[16]殷躍平.汶川地震與滑坡災(zāi)害概論[M].地質(zhì)出版社,2009.

AGING STABIL ITY THEORY OF LANDSL IDE CREEP AND BASIC CONCEPTIONOF LANDSL IDE PRED ICTION

YANG Ren-guang
(CHina Sihai Holdings Co.,L td,Beijing100081,China)

Formation of landslide show s its aging characteristics w hile it is revealed by the existing researches that the aging characteristics of landslide can help to p redict and fo recast the landslide p recisely.On the basis of research findings achieved by Maxwell,Kelvin and Saitou Susumu Takashi etc.,the viscoelastic lag coefficient of rock-soilmass is introduced in this paper, together w ith the p ractical monito ring geometric locus of test piece of rock and soil materials under different strained conditions,and then the fitting p rocess is adop ted.On the p remiseof follow ing the continuitiesof both defo rmation and degreeof defo rmation of knee of curve,the binary nonlinear transcendental algebraic equations concerning time is deduced,and also the time of aging conjugate point and destabilization time caused by severe sliding are calculated,finally the conclusion in w hich the ageing of landslide creep can be p redicted is obtained.

landslide;creep;viscoelastic lag coefficient;aging conjugate point;aging stability

P642.22

:A

1006-4362(2010)02-0071-03

楊人光(1939- ),男,原中科院力學(xué)所研究員,現(xiàn)任中國(guó)四?？毓捎邢薰旧鷳B(tài)災(zāi)害測(cè)控中心總工程師,長(zhǎng)期致力于巖土理論與實(shí)踐,主編的《控制爆炸技術(shù)研究與應(yīng)用》獲國(guó)家一等獎(jiǎng)。

2010-03-26改回日期:2010-04-29