張文杰
我國是世界上遭受自然災(zāi)害最為嚴(yán)重的國家之一。自20世紀(jì)80年代以來,隨著經(jīng)濟(jì)的快速發(fā)展,人類工程活動日益頻繁,在多種自然因素及人為活動的影響下,滑坡災(zāi)害成為各類災(zāi)害中頻率最高、危害性最強(qiáng)的地質(zhì)災(zāi)害類型之一,且有逐年加劇趨勢。因此,如何精準(zhǔn)有效地預(yù)防滑坡災(zāi)害已成為地質(zhì)災(zāi)害研究領(lǐng)域的熱點。
近年來,滑坡位移預(yù)測逐漸成為滑坡預(yù)警的有效手段之一,它允許研究人員在滑坡位移監(jiān)測數(shù)據(jù)分析的基礎(chǔ)上,采用相關(guān)數(shù)學(xué)模型進(jìn)行擬合,從而更加精準(zhǔn)地進(jìn)行滑坡位移預(yù)測。這種方法易于實現(xiàn)且精度較高,因此發(fā)展迅速。目前,已有很多關(guān)于滑坡位移預(yù)測的智能算法研究,如BP神經(jīng)網(wǎng)絡(luò)、支持向量機(jī)、灰色模型等。
鑒于僅通過單個模型進(jìn)行滑坡預(yù)測具有局限性,有學(xué)者嘗試將多個模型組合起來共同進(jìn)行滑坡預(yù)測,提出了混合預(yù)測法。本文則采用模擬退火算法(Simulated Annealing,SA)結(jié)合BP神經(jīng)網(wǎng)絡(luò)的混合模型對實例邊坡位移進(jìn)行預(yù)測,算例結(jié)果證明了方法的有效性。
模擬退火算法最早于1953年由N. Metropolis等人提出,是一種以蒙特卡洛迭代求解策略為基礎(chǔ)的隨機(jī)尋優(yōu)算法。該方法由兩部分組成:Metropolis算法和退火過程。它的基本原理是模擬物理學(xué)中固體物質(zhì)的退火過程,當(dāng)溫度升高時,固體中的分子內(nèi)能增加,并做無序運動;當(dāng)溫度降低時,分子內(nèi)能隨溫度下降而減小,并逐漸趨向于有序運動;在任意給定的溫度條件下,固體材料內(nèi)部都能夠?qū)崿F(xiàn)熱平衡,此時可將其看作是一個局部最優(yōu)解;在最低溫度時,固體中的分子內(nèi)能最小,此時可將熱平衡視為全局最佳。
BP神經(jīng)網(wǎng)絡(luò)是1986年由Rumelhart等人提出的一種多層前向反饋神經(jīng)網(wǎng)絡(luò),通過模擬人類的大腦結(jié)構(gòu)與功能來處理多變量問題,具有良好的學(xué)習(xí)能力,能夠較好地處理非線性問題?;镜腂P神經(jīng)網(wǎng)絡(luò)是由三層拓?fù)浣Y(jié)構(gòu)組成:輸入層、隱藏層和輸出層,其結(jié)構(gòu)形式如圖1所示。
BP神經(jīng)網(wǎng)絡(luò)調(diào)整權(quán)值和閾值都是采用梯度下降法,容易使其陷入局部最優(yōu)且收斂速度較慢。模擬退火算法是一種新的隨機(jī)搜索方法,既能實現(xiàn)全局搜索又能進(jìn)行局部優(yōu)化。借助模擬退火算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò),將梯度下降法與隨機(jī)擾動相結(jié)合,采用一種新的方式調(diào)節(jié)權(quán)值和閾值,以達(dá)到解決BP神經(jīng)網(wǎng)絡(luò)容易陷入局部最優(yōu)、收斂慢等問題的目的,從而實現(xiàn)滑坡位移的有效預(yù)測。其計算流程見圖2。
>圖1 BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)圖
> 圖2 SA-BP模型流程圖
>圖3 白水河滑坡GPS監(jiān)測點布置平面圖
> 圖4 不同模型的預(yù)測值與監(jiān)測值對比
白水河滑坡位于三峽庫區(qū)秭歸縣境內(nèi),距三峽大壩約56km處?;w平均厚度約30m,南北方向長約600m,東西方向?qū)捈s700m,體積約為1260×104m3。該滑坡體為特大型老滑坡松散堆積體,出露地層為下侏羅統(tǒng)香溪群中厚層狀砂巖夾、薄層狀泥巖;滑坡后沿高程為410m,以巖土邊界為界,其前沿高程約為70m,現(xiàn)已沒入庫水位以下,東西兩側(cè)以基巖山脊為界,總體坡度約30°;滑坡滑體主要由崩、坡積物及滑坡堆積物組成,滑帶以含碎石或含角礫粉質(zhì)黏土為主。自2003年三峽工程蓄水以來,相關(guān)單位對其進(jìn)行了專業(yè)監(jiān)測,GPS監(jiān)測點布置如圖3所示。
2007年1月至2012年12月,ZG118站點監(jiān)測到的水平累計位移、水位和降雨量數(shù)據(jù),水庫水位和降雨量變化是導(dǎo)致白水河滑坡的主要破壞因素。
選擇監(jiān)測點ZG118的水平累計位移為研究對象,取2007年1月至2012年1月的數(shù)據(jù)作為預(yù)測模型的訓(xùn)練樣本集,2012年2月至2012年12月的數(shù)據(jù)作為預(yù)測檢驗集。
取單月平均降雨量為降雨因子;月均水位為水位因子;時間函數(shù)t為時效因子。將3個影響因子作為輸入,累計位移作為輸出,進(jìn)行基于SA-BP神經(jīng)網(wǎng)絡(luò)的滑坡位移預(yù)測分析。
將SA-BP神經(jīng)網(wǎng)絡(luò)模型與BP神經(jīng)網(wǎng)絡(luò)模型作對比,滑坡位移預(yù)測與監(jiān)測值如圖4所示。
綜合上述分析可知,SA-BP神經(jīng)網(wǎng)絡(luò)模型對白水河滑坡的位移預(yù)測效果比BP模型要好,其誤差值更低,精度更高。
本文結(jié)合SA算法和BP神經(jīng)網(wǎng)絡(luò),建立基于SA算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的滑坡位移預(yù)測模型。通過SA算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的初始連接權(quán)值和閾值,尋找BP神經(jīng)網(wǎng)絡(luò)的全局最優(yōu)權(quán)值和閾值。將所建立的滑坡位移預(yù)測模型和傳統(tǒng)BP神經(jīng)網(wǎng)絡(luò)模型分別應(yīng)用于工程實例中,結(jié)果表明SA-BP神經(jīng)網(wǎng)絡(luò)預(yù)測邊坡位移誤差較小,具有更高的精度。