炕金賓 高桂洋 中交（長(zhǎng)沙）建設(shè)有限公司

滑坡預(yù)警預(yù)測(cè)屬于邊坡穩(wěn)定性研究中的一個(gè)重要分支。因坡體自身結(jié)構(gòu)特性較復(fù)雜，故對(duì)滑坡體滑動(dòng)位置的預(yù)測(cè)具有一定的難度。通過(guò)選取合適的模型，對(duì)邊坡變形進(jìn)行預(yù)測(cè)，這種研究手段具有較高的理論及工程應(yīng)用價(jià)值?；疑到y(tǒng)其相關(guān)理論的提出較早，鄧聚龍首次提出并對(duì)其進(jìn)行了較為詳盡的闡述?；疑到y(tǒng)理論為解決信息完整度較差這一問(wèn)題，以數(shù)學(xué)方法為基礎(chǔ)結(jié)合系統(tǒng)論、信息論及控制論觀點(diǎn)對(duì)信息進(jìn)行處理。該理論的本質(zhì)是選取已有的有效信息，通過(guò)一系列的分析實(shí)現(xiàn)對(duì)信息演化趨勢(shì)的精確預(yù)測(cè)，由于對(duì)觀測(cè)數(shù)據(jù)沒(méi)有特殊的限制，其適用范圍較廣運(yùn)用程度較高。本文使用灰色系統(tǒng)理論并結(jié)合現(xiàn)場(chǎng)變形監(jiān)測(cè)數(shù)據(jù)建立炭質(zhì)頁(yè)巖路塹邊坡變形預(yù)測(cè)模型，對(duì)邊坡的變形進(jìn)行動(dòng)態(tài)的預(yù)測(cè)，能較準(zhǔn)確地對(duì)炭質(zhì)頁(yè)巖的變形趨勢(shì)進(jìn)行預(yù)報(bào)。

1.GM（1，1）預(yù)測(cè)模型的建立

灰色模型在本質(zhì)上是運(yùn)用生成函數(shù)和灰色微分方程，通過(guò)一系列的分析、推導(dǎo)建立模型。在模型建立時(shí)可供選擇的原始數(shù)據(jù)樣本較多，但建模實(shí)際運(yùn)用的數(shù)據(jù)有限。不同的數(shù)據(jù)所得到的結(jié)果及誤差差異性較大。為了最大可能保證預(yù)測(cè)的準(zhǔn)確性，本文采用新陳代謝GM（1，1）模型進(jìn)行建模分析。新陳代謝GM（1，1）模型的一般過(guò)程如下：

（1）用X=[x(1),x(2),......,x(n)]建立全信息GM（1，1）模型；

（2）?＞1，令X=[x(k),x(k+1),......,x(n)]建立的模型為部分信息GM（1，1）模型；

（3）若x(n+1)為最新數(shù)據(jù)，把x(n+1)加入X，且去掉老信息x(1)，然后由新的X=[x(2)，x(3)，......，x(n+1)]）建立的模型稱為新信息GM（1，1）模型，即新陳代謝GM（1，1）模型；

（4）對(duì)模型精度進(jìn)行檢驗(yàn)。

2.工程實(shí)例研究

2.1 工程概況

本文將某公路炭質(zhì)頁(yè)巖路塹邊坡作為研究邊坡。該邊坡最大挖高41.0米。坡體巖層傾向與坡向一致，對(duì)坡體穩(wěn)定性不利。地層巖性為：1）覆蓋層，主要為強(qiáng)風(fēng)化炭質(zhì)頁(yè)巖，呈棱角狀，粒徑10cm左右；2）基巖，主要是強(qiáng)中風(fēng)化的炭質(zhì)頁(yè)巖，裂隙發(fā)育，巖芯多呈碎塊狀，巖層產(chǎn)狀為65°∠20°，呈黑色。巖體工程等級(jí)為IV級(jí)軟巖，節(jié)理裂隙較發(fā)育，對(duì)邊坡的穩(wěn)定性影響不大。從現(xiàn)場(chǎng)勘探可知：坡腳處較陡峭；邊坡按正常坡率放坡，其中第一級(jí)邊坡坡率為l:0.75，第二～三級(jí)邊坡坡率為1:1.00，平臺(tái)寬均為2米；因?yàn)檫吰滤诘沫h(huán)境十分復(fù)雜且所處位置較高，所以該邊坡存在很大的安全隱患。如圖1所示。

圖1 某公路炭質(zhì)頁(yè)巖路塹邊坡示意圖

邊坡區(qū)地表水發(fā)育程度較低，地下水以基巖間裂隙水為主，主要接受大氣降水補(bǔ)給。邊坡地處亞熱帶季風(fēng)氣候區(qū)，7、8月份的時(shí)候高溫且多雨，其最高溫度能夠達(dá)到41.4℃左右。由于該地區(qū)為低山丘陵區(qū)，因此在氣流作用下容易形成降雨，所以邊坡所處區(qū)域雨量充沛。通過(guò)實(shí)地探訪并查閱相關(guān)資料，對(duì)該地區(qū)近十年的降雨數(shù)據(jù)進(jìn)行整理，結(jié)果如表1所示。

表1 某地區(qū)2004年-2014年降雨數(shù)據(jù)匯總表

2.2 實(shí)例分析

將灰色新陳代謝GM（1，1）模型運(yùn)用于實(shí)際的工程，并對(duì)其預(yù)測(cè)準(zhǔn)確性和適應(yīng)性進(jìn)行進(jìn)一步分析。

先整理了研究邊坡的變形監(jiān)測(cè)數(shù)據(jù)，然后基于灰色新陳代謝GM（1，1）模型進(jìn)行建模分析。共整理了測(cè)點(diǎn)A1的34組變形數(shù)據(jù)，具體數(shù)據(jù)如表2所示。

表2 邊坡A1測(cè)點(diǎn)變形監(jiān)測(cè)匯總表

在得到上述結(jié)果之后，隨機(jī)抽取總變形中的30組數(shù)據(jù)作為預(yù)測(cè)的原始數(shù)據(jù)。根據(jù)新陳代謝GM（1，1）模型的建模過(guò)程，每7個(gè)為一組不斷進(jìn)行新老替換直到30組數(shù)據(jù)用完。為了使計(jì)算方便，本文借助Matlab軟件進(jìn)行編程分析計(jì)算，在Matlab軟件上運(yùn)行源代碼，得到預(yù)測(cè)結(jié)果。由于篇幅有限，這里只列出最后得到的新陳代謝GM（1，1）模型的時(shí)間響應(yīng)函數(shù)：

得到的最終結(jié)果匯總?cè)绫?、圖2所示。

精度檢驗(yàn)作為評(píng)價(jià)預(yù)測(cè)模型可行性的重要標(biāo)準(zhǔn)，本文在得到上述相關(guān)數(shù)據(jù)后，按擬定精度檢驗(yàn)流程、方法，進(jìn)行驗(yàn)算，具體結(jié)果如下所示：

2）計(jì)算可知，模型的灰色絕對(duì)關(guān)聯(lián)度ε=0.999＞0.900，達(dá)到一級(jí)標(biāo)準(zhǔn)；

從圖2、表3可知：實(shí)測(cè)值和預(yù)測(cè)值較為一致，最大相對(duì)誤差為1.86mm，擬合結(jié)果較好；通過(guò)計(jì)算可知模型預(yù)測(cè)精度滿足要求。所以，該模型可用于預(yù)測(cè)邊坡的后期變形。

圖2 邊坡A1測(cè)點(diǎn)實(shí)測(cè)值與預(yù)測(cè)值對(duì)比圖

表3 邊坡A1測(cè)點(diǎn)累計(jì)變形實(shí)測(cè)值與預(yù)測(cè)值對(duì)比表 單位(mm)

完成上述分析之后，利用后面的4組數(shù)據(jù)對(duì)模型進(jìn)行檢驗(yàn)分析。同樣利用Matlab軟件進(jìn)行計(jì)算，最后得到31～34號(hào)的預(yù)測(cè)結(jié)果如表4及圖3所示。

圖3 邊坡A1測(cè)點(diǎn)實(shí)測(cè)值與預(yù)測(cè)值對(duì)比圖

表4 邊坡A1測(cè)點(diǎn)累計(jì)變形值實(shí)測(cè)值與預(yù)測(cè)值對(duì)比表 單位(mm)

分析可知，序號(hào)31測(cè)點(diǎn)變形實(shí)測(cè)值與預(yù)測(cè)值誤差為0.18mm，序號(hào)32測(cè)點(diǎn)變形實(shí)測(cè)值與預(yù)測(cè)值誤差為0.28mm，序號(hào)33測(cè)點(diǎn)變形實(shí)測(cè)值與預(yù)測(cè)值誤差為0.09mm，序號(hào)34測(cè)點(diǎn)變形實(shí)測(cè)值與預(yù)測(cè)值誤差為0.51mm，實(shí)測(cè)值和預(yù)測(cè)值最大誤差為0.51mm，最大相對(duì)誤差0.39%，模型和實(shí)測(cè)得到的結(jié)果擬合較好，兩者誤差較小且精度滿足要求。

所以，通過(guò)新陳代謝GM（1，1）模型，對(duì)研究邊坡整理好的34組變形數(shù)據(jù)進(jìn)行分析。隨機(jī)抽取4組數(shù)據(jù)進(jìn)行模型驗(yàn)證，剩下的30組數(shù)據(jù)當(dāng)作建模的原始數(shù)據(jù)。在預(yù)測(cè)過(guò)程中以7組數(shù)據(jù)為一個(gè)系列進(jìn)行建模，不斷剔除老數(shù)據(jù)加入新數(shù)據(jù)。得到預(yù)測(cè)結(jié)果之后利用預(yù)留的4組數(shù)據(jù)進(jìn)行驗(yàn)證分析。對(duì)所得結(jié)果分析可知新陳代謝GM（1，1）模型和實(shí)測(cè)得到的結(jié)果擬合較好，最大相對(duì)誤差為1.86mm，兩者誤差較小且精度滿足要求，可以較好的對(duì)邊坡的變形進(jìn)行預(yù)測(cè)，而且效果理想。

3.結(jié)語(yǔ)

本文基于灰色系統(tǒng)理論對(duì)邊坡位變形進(jìn)行分析，得到如下結(jié)論：

（1）新陳代謝GM（1，1）模型是在傳統(tǒng)GM（1，1）模型基礎(chǔ)上的改進(jìn)模型，其利用率高、精度高、結(jié)果更準(zhǔn)確。

（2）基于新陳代謝GM（1，1）模型，對(duì)所研究邊坡變形數(shù)據(jù)進(jìn)行一系列處理分析，所得結(jié)果和實(shí)測(cè)結(jié)果擬合較好，因此該模型可以用作該邊坡變形預(yù)測(cè)模型。