楊 帆 許 強 范宣梅 葉 微

(①地質(zhì)災害防治與地質(zhì)環(huán)境保護國家重點實驗室(成都理工大學) 成都 610059)(②四川大學商學院 成都 610065)

0 引 言

滑坡的預測預警是國內(nèi)外學者研究的熱點問題。自20世紀60年代日本學者齋藤(Satio， 1965)提出滑坡位移經(jīng)驗公式預報法以來，眾多學者對斜坡變形破壞及滑坡預報方法進行不斷的探索(張倬元等， 1994)。目前，滑坡預報手段通常有兩種方法：一是基于降雨閾值方法進行預報，通常在特定地區(qū)范圍較為適用(伍宇明等， 2014; Palladino et al.， 2017)； 二是基于現(xiàn)場實測數(shù)據(jù)進行預報(包括位移變形、孔隙水壓、含水率及微震監(jiān)測等方法)(Busslinger， 2009; Li et al.， 2012)。由于后者監(jiān)測數(shù)據(jù)具有實時動態(tài)的特點，在一定程度上具有更高的預報精度。針對滑坡位移預報方法多種多樣總體可分為確定性模型、統(tǒng)計預報模型和非線性預報模型3類(蘇愛軍等， 1990； 廖野瀾等， 1996； 黃潤秋等， 1997； 蘭恒星等， 2000； 李秀珍等， 2003； 汪洋等， 2004； 尹光志等， 2007； 許強等， 2008； 杜鵑等， 2009； Li et al.， 2012； 董秀軍等， 2015； 苗發(fā)盛等， 2016； 郭子正等， 2018； Li et al.， 2018； 楊背背等， 2018； Zhu et al.， 2018)。

隨著現(xiàn)代科學技術的不斷發(fā)展，滑坡位移監(jiān)測及預測預報手段越來越豐富，監(jiān)測精度得到顯著提升?；骂A測預報研究已從之前的單因素預報方法，逐漸轉(zhuǎn)變?yōu)榛趯崟r監(jiān)測數(shù)據(jù)，結(jié)合定性、定量及現(xiàn)代數(shù)值預報技術等多種手段的綜合動態(tài)預報，并開發(fā)自動化預警系統(tǒng)，實現(xiàn)主動防災減災的新階段(許強等， 2004， 2014； 蘇白燕等， 2018)。選用適合的方法對位移-時間曲線的監(jiān)測數(shù)據(jù)進行分析，將其應用于滑坡監(jiān)測預警系統(tǒng)作為預警判據(jù)，對防災減災具有十分重要的意義。

隨著模式識別與智能計算的發(fā)展，一些數(shù)據(jù)挖掘的方法(如神經(jīng)網(wǎng)絡模型、灰色模型、時間序列分析模型、小波分析模型、極限學習機模型、支持向量機、混沌序列模型等)被應用到滑坡位移預測研究中(劉曉等， 2005； 杜鵑等， 2009； 徐峰等， 2011； 張俊等， 2015； 周超等， 2015； 闕金聲等， 2016； 鄧冬梅等， 2017； 郭子正等， 2018； Li et al.， 2018； Zhu et al.， 2018)。劉曉等(2005)用BP神經(jīng)網(wǎng)絡方法提取位移趨勢項，結(jié)合時間序列ARMA進行滾動預測； 杜娟等(2009)結(jié)合BP神經(jīng)網(wǎng)絡分別對滑坡趨勢項和周期項位移進行預測，根據(jù)時間序列提取趨勢項并預測，以BP神經(jīng)網(wǎng)絡方法對周期項位移進行預測。徐峰等(2011)結(jié)合灰色模型和自回歸AR模型的滑坡位移預測方法，采用灰色模型和自回歸AR模型進行預測研究； 張俊等(2015)采用多項式及PSO-SVR對趨勢項和周期項位移進行預測； Li et al.(2018)用小波函數(shù)聯(lián)合ELM、OS-ELM建立滑坡位移預測模型。均取得較好的結(jié)果。鄧冬梅等(2017)在張俊等(2015)成果基礎上，加入集合經(jīng)驗模態(tài)分解(EEMD)的粒子群優(yōu)化-支持向量機回歸方法進行位移預測研究； 在前人研究成果基礎上，郭子正等(2018)基于小波分析的三角函數(shù)模型進行周期項位移預測研究； 楊背背等(2018)運用基于時間序列與長短時基于網(wǎng)絡(LSTM)進行位移預測研究。

上述方法中，杜娟等(2009)存在研究樣本年份過短、單個多項式對趨勢項位移進行擬合精度不高、灰色系統(tǒng)對于非線性序列研究誤差較大等問題(徐峰等， 2011)； 人工神經(jīng)網(wǎng)絡缺乏完備的理論基礎，存在模型結(jié)構(gòu)難確定、精度難保證及泛化能力不高等缺點(張俊等， 2015)。支持向量機是Corinna Cortes和Vapnik于1995年首次提出的一種基于統(tǒng)計學習模型的監(jiān)督式機器分類技術(Vapnik， 1995)，有堅實的數(shù)學理論基礎，可以較好地解決小樣本、非線性和局部最小等問題(楊帆等， 2017)。在模型訓練過程中特征子集的選取以及參數(shù)的設定對分類結(jié)果有顯著的影響(Zhao et al.， 2011)。常用支持向量機參數(shù)優(yōu)化算法有遺傳算法(GA)、蟻群算法(ACO)、粒子群算法(PSO)及模擬退火算法等(Shelokar et al.， 2004; 楊成祥等， 2005； Huang et al.， 2006; Niknam et al.， 2010; Zhao et al.， 2011; 劉愛軍等， 2013； 吳華鋒等， 2013； 陳健飛等， 2016)。但上述傳統(tǒng)算法的性能很大程度上依賴其選用的參數(shù)，同時，在解決多峰值問題的尋優(yōu)過程中，結(jié)果容易陷入局部最優(yōu)問題，從而直接影響預測結(jié)果的準確性(陳健飛等， 2016)。人工蜂群算法具有較好的全局搜索能力，魯棒性強，可較好地解決上述問題，避免陷入局部最優(yōu)。

針對上述問題，本文提出了一種基于時間序列與人工蜂群算法支持向量回歸機(ABC-SVR)的滑坡位移預測模型。以三峽庫區(qū)白水河滑坡為例，用時間序列加法模型結(jié)合平滑法中的移動平均法，將監(jiān)測點位移數(shù)據(jù)分解為兩項：滑坡位移趨勢項和滑坡位移周期項。用多項式最小二乘法對位移趨勢項數(shù)據(jù)進行分段擬合，由人工蜂群支持向量機基于位移周期項數(shù)據(jù)進行訓練及預測。同時選取長江庫水位、本月累積降雨、前一月累積雨量、前兩月累積雨量、前一月庫水位變化量、前兩月庫水位變化量與監(jiān)測點年累計位移量7項因子結(jié)合灰色系統(tǒng)關聯(lián)分析法研究各因素與周期項位移間的關聯(lián)性。將趨勢項、周期項位移預測值疊加得到總位移預測值并同實際監(jiān)測位移進行對比。通過與BP神經(jīng)網(wǎng)絡和PSO-SVR模型預測結(jié)果的對比，該模型表現(xiàn)出較好的預測結(jié)果。

1 滑坡位移預測方法介紹

1.1 時間序列加法模型

地殼表層是個復雜的開放系統(tǒng)，滑坡在其發(fā)展演化的過程中會受到多種因素的影響與控制，其位移時間序列理論可概括為3部分(楊叔子， 1991)：包括受地形、地質(zhì)構(gòu)造等內(nèi)部地質(zhì)因素影響的趨勢項； 受降雨、地下水及庫水位變動等外部作用影響的周期項； 以及一個隨機項。 由于在現(xiàn)有技術條件下，難以對隨機項因素進行準確評估，因此，在本文的滑坡位移預測研究中暫不考慮隨機項因素，時間加法序列模型可以簡化為：

X(t)=φ(t)+η(t)

(1)

式中，X(t)為時間序列函數(shù)，φ(t)為位移趨勢項，η(t)為位移周期項。

1.1.1 趨勢項位移提取

滑坡位移趨勢項代表了滑坡位移長期發(fā)展的趨勢。采用移動平均法提取趨勢項位移，可削弱或消除時間序列中季節(jié)和不規(guī)則周期變動的影響。計算式如下：

(2)

1.1.2 周期項位移提取

由時間序列加法模型公式(式1)可知，用滑坡總位移減去滑坡位移趨勢項即可得到滑坡位移的周期項。

1.2 支持向量機

支持向量機，是1995年由Corinna Cortes和Vapnik首次提出的一種基于統(tǒng)計學習VC維(Vapnik-Chervonenkisdimension)和結(jié)構(gòu)風險最小化理論的監(jiān)督學習分類技術。該模型采用核映射思想，在解決小樣本、非線性、高維數(shù)的問題時表現(xiàn)出很好的優(yōu)勢(Vapnik et al.， 1995)，能較好的克服神經(jīng)網(wǎng)絡等方法的固有缺陷，在處理非線性問題時具有一定的優(yōu)越性(楊曉偉等， 2013； 楊帆等， 2017)。

支持向量機(SVM)可用于解決線性可分問題，而對于線性不可分問題，利用核函數(shù)將低維數(shù)據(jù)映射到高維空間，從而低維空間上線性不可分問題便轉(zhuǎn)化為高維空間上的線性可分問題，求解后將其映射到低維空間。假設數(shù)據(jù)集為S，原始空間集Rn，從低維到高維映射關系為：φ：x→φ(x)，具體運算步驟為：

(1)尋找一個核函數(shù)K(s，t)使得

K(xi，xj)=(φ(xi)，φ(xj))，

(3)

(2)構(gòu)造并求解凸二次規(guī)劃問題：

(4)

(3)計算b

(5)

(4)構(gòu)造非線性函數(shù)

(6)

本文核函數(shù)選用徑向基核函數(shù)(RBF)：

K(xi，xj)=exp(-|xi-xj|2/2γ2)

(7)

1.3 人工蜂群算法

人工蜂群算法(Artificial Bee Colony)是一種受蜂群行為啟發(fā)的代數(shù)優(yōu)化算法，由Karaboga等人在2005年提出(楊淑瑩等， 2015)。該算法具有較強的發(fā)現(xiàn)最優(yōu)解、信息正反饋、收斂速度快、魯棒性強及易于和其他方法結(jié)合等優(yōu)點，能較大程度上避免陷入局部最優(yōu)的問題。算法步驟如下：

1.3.1 蜜源初始化

初始化種群參數(shù)，隨機生成SN個可行解(與雇傭蜂數(shù)量相同)，計算適應度(Fitness)。隨機產(chǎn)生可行解公式：

(8)

式中，xi(i=1， 2，…，SN)為D維向量，D是待優(yōu)化參數(shù)的數(shù)目。j∈(1， 2，…，D)。

1.3.2 蜜源更新公式

引領蜂與跟隨蜂依據(jù)下列公式記錄截止目前的最優(yōu)值，且在當前蜜源鄰域內(nèi)搜索最優(yōu)值：

vi， j=xi， j+φi， j(xi， j-xk， j)

(9)

式中，j∈(1， 2，…，D)，k∈(1， 2，…，SN)，k隨機生成且k≠i，φi， k∈[-1， 1]的隨機數(shù)。

1.3.3 選擇雇傭蜂概率公式

(10)

式中，fit(xi)為第i個Fitness解所對應的富源度。

1.3.4 偵查蜂的產(chǎn)生

若有蜜源連續(xù)未被更新的次數(shù)大于limit時，則對該蜜源進行重新初始化。算法流程圖(圖 1)：

圖 1 程序流程圖Fig. 1 Flow chart of the procedure

2 三峽庫區(qū)白水河滑坡位移監(jiān)測分析

在我國三峽庫區(qū)，由于季節(jié)性降水的分配不均及水庫水位的顯著波動等因素，水庫沿岸滑坡災害頻發(fā)，據(jù)統(tǒng)計，截至2014年長江三峽庫區(qū)崩滑地質(zhì)災害點多達5000余處(許強等， 2014)，對長江沿岸居民生命財產(chǎn)安全造成嚴重威脅。因此對持續(xù)變形的滑坡進行位移預測被認為是減少或避免滑坡災害造成損失的重要且經(jīng)濟的方法(殷躍平， 2003； 許強等， 2004； 許強，2012; Zhu et al.，2018; Ma et al.，2018)。

2.1 滑坡工程地質(zhì)概況

白水河滑坡位于三峽庫區(qū)寬河谷段，長江主干道南岸(凸岸)，距三峽大壩約56km，屬秭歸縣沙鎮(zhèn)溪鎮(zhèn)樂豐村。秭歸縣地處亞熱帶季風氣候區(qū)，降雨連續(xù)集中且雨季多暴雨。地勢南高北低，為侏羅系下統(tǒng)香溪群砂巖組成的順向岸坡，巖層產(chǎn)狀15°∠36°，滑坡前緣高程約70m，后緣高程約410m，東西兩側(cè)為近南北向山脊，北側(cè)為白水河單面山體。該滑坡南北向約600m，東西向約700m，滑坡總體坡度約30°，滑坡總面積0.42km2，滑坡體積126×105￣￣m3(圖 2)。

圖 2 白水河滑坡平面示意圖Fig. 2 Plane sketch map of Baishuihe Landslide

白水河滑坡歷史上就曾發(fā)生過滑動，多年來，隨著庫水位的不斷波動及汛期強降雨，該滑坡變形特征較為強烈。2003年6月隨著三峽水庫蓄水至135m后，該滑坡出現(xiàn)明顯的位移增長趨勢，特別是2007年6月后，受強降雨和庫水位等聯(lián)合作用，滑坡變形驟然增加(圖 3)。ZG93和ZG118監(jiān)測點從2003年6月開始監(jiān)測。2005年5月和10月又陸續(xù)加設了XD1、XD2和XD3、XD4監(jiān)測點。

本文選用ZG93及ZG118監(jiān)測點數(shù)據(jù)進行研究并建立模型。因為相比于其他監(jiān)測點，ZG93和ZG118持續(xù)監(jiān)測時間長，且均位于白水河滑坡體中部區(qū)域，可以較好的記錄和反映滑坡位移運動的整體趨勢。

圖 3 2003～2013年白水河滑坡月降雨量及長江水位與累計位移關系圖Fig. 3 Displacements， rainfall and reservoir level during the period 2003～2013

3 預測模型的應用

滑坡變形的產(chǎn)生受內(nèi)部地質(zhì)因素和外部誘發(fā)因素共同影響，本文以三峽庫區(qū)白水河滑坡為研究對象，選用ZG93和ZG118位移監(jiān)測點自2003年6月至2013年3月間的監(jiān)測數(shù)據(jù)進行分析，將總位移分解為趨勢項和周期項位移，分別用多項式最小二乘法、人工蜂群支持向量機模型進行擬合訓練及預測分析。

3.1 趨勢項位移分析

移動平均可以削弱或消除原時間序列中季節(jié)和不規(guī)則周期性的變動，從而分離出趨勢項。移動平均法中取周期長度M=12個月，根據(jù)移動平均法公式(式(2))提取出滑坡位移趨勢項，利用多項式最小二乘法對其進行分段擬合。

由于趨勢項提取值曲線呈現(xiàn)明顯的分段特征，故本文將其分為3個階段(2004年6月～2007年5月， 2007年6月～2008年7月， 2008年8月～2011年6月)進行研究。ZG93和ZG118滑坡監(jiān)測點趨勢項位移的提取值和預測值(圖 4)。

圖 4 趨勢項位移的提取及預測(ZG93、ZG118)Fig. 4 Extraction and forecasting the displacement of trend term(ZG93、ZG118)

利用Matlab2010年對趨勢項位移提取值進行多項式最小二乘法擬合(圖 4)，擬合方程及精度結(jié)果如下：

監(jiān)測點ZG93：

監(jiān)測點ZG118：

3.2 周期項位移分析

3.2.1 周期項位移提取

根據(jù)時間序列加法模型公式(式1)，周期項位移等于總位移量減去趨勢項位移。因此，監(jiān)測點ZG93與ZG118的周期項位移提取值如圖 5所示。

圖 5 監(jiān)測點周期項位移提取值Fig. 5 Periodic displacement value of monitoring points

3.2.2 周期項位移預測

降雨和庫水位的變化是滑坡位移的兩個重要的誘發(fā)因素(杜鵑等， 2009； 楊背背等， 2018)。通過選取合適的影響因子，研究分析各誘發(fā)因素與滑坡位移之間的響應關系，可以更好地對滑坡位移進行研究和準確評價。本文在比對了前人研究資料的基礎上(杜鵑等， 2009； 張俊等， 2015)，選取長江庫水位、本月累積降雨、前一月累積雨量、前兩月累積雨量、前一月庫水位變化量、前兩月庫水位變化量與監(jiān)測點年累計位移7項因素，結(jié)合灰色系統(tǒng)關聯(lián)度分析法研究各因子同滑坡監(jiān)測點周期項位移間的關聯(lián)度，取分辨系數(shù)為0.5，從而得到監(jiān)測點周期性位移與各影響因子間的關聯(lián)度rk如表 1所示。其中，長江庫水位、前兩月降雨、前一月降雨、本月累積降雨同監(jiān)測點ZG93周期項位移監(jiān)測數(shù)據(jù)之間的關聯(lián)度分別為0.8525、0.8203、0.7984、0.7912，并且后3項數(shù)值比較接近； 與監(jiān)測點ZG118表現(xiàn)出的現(xiàn)象近似。反映出監(jiān)測點“階躍式”累計位移同長江庫區(qū)水位有密切的關系； 同時，降雨因素在滑坡位移中稍有滯后性。

表 1 滑坡位移周期項與各因子的相關性Table 1 Correlation between periodic terms of landslide displacement and various factors

3.3 模型預測及對比分析

用可優(yōu)化選參的人工蜂群算法(ABC)與支持向量回歸機(SVR)結(jié)合，構(gòu)建ABC-SVR滑坡位移預測模型，ABC算法可尋求最優(yōu)的懲罰因子C及核函數(shù)參數(shù)g，可提高模型的預測精度(許國根， 2012； 楊景明等， 2014； 卓金武， 2014； 高雷阜等， 2016)。以白水河滑坡2004年6月～2011年6月間的滑坡周期項位移數(shù)據(jù)為訓練樣本，釆用ABC-SVR模型對2011年7月～2013年6月間的周期性位移進行預測。詳細過程如下：

(1)樣本數(shù)據(jù)預處理。以離差標準化(Min-Max Normalization)對原始數(shù)據(jù)做線性變換，使結(jié)果映射到[0，1]之間。

(2)模型參數(shù)初始化。規(guī)定蜂群規(guī)模為NP=20，蜜源數(shù)量為NP/2，最大搜索次數(shù)limit，最大迭代次數(shù)maxcycle，核函數(shù)使用RBF核函數(shù)，支持向量回歸機待優(yōu)化參數(shù)的取值區(qū)間。

(3)訓練ABC-SVR模型，用ABC算法搜尋SVR懲罰因子C及核函數(shù)參數(shù)g的最優(yōu)值，用得到的最優(yōu)參數(shù)值帶入支持向量回歸機中進行學習訓練和預測。并且用ABC-SVR模型與文獻中的BP神經(jīng)網(wǎng)絡模型、PSO-SVR模型的預測結(jié)果進行對比，對比結(jié)果見表 2。

表 2 不同模型預測精度對比Table 2 Comparisons of prediction accuracy of different models

預測結(jié)果顯示，監(jiān)測點ZG93的均方根誤差MSE及平均絕對百分誤差MAPE分別為0.0311， 0.0768，相關系數(shù)(Correlation coefficient)為0.96，監(jiān)測點ZG118的均方根誤差MSE及平均絕對百分誤差MAPE分別為0.0418， 0.0828，相關系數(shù)(Correlation coefficient)為0.95，預測精度較高。

從上表可知，PSO-SVR在時間上優(yōu)于ABC-SVR，但預測效果低于ABC-SVR。BP神經(jīng)網(wǎng)絡預測模型運行時間較長，預測結(jié)果不如其他兩個模型。綜上所述，ABC-SVR預測模型預測效果最優(yōu)，說明其具有較好的尋優(yōu)性能及預測效果。

3.4 累計位移預測

根據(jù)時間序列加法模型公式(式(1))可知，監(jiān)測點趨勢項與周期項位移預測值的數(shù)量和即為監(jiān)測點滑坡累計總位移的預測值。將滑坡累計位移預測值與實測值進行對比(圖 6)。結(jié)果顯示，預測結(jié)果與實測值變化趨勢較為一致，預測結(jié)果較好。ZG93和ZG118累計位移預測擬合優(yōu)度值R2分別為0.97、0.96，平均絕對百分誤差MAPE分別為0.083和0.078，說明基于時間序列人工蜂群支持向量機模型具有較好的預測精度，可以用于滑坡位移的預測研究。

圖 6 監(jiān)測點累積位移的預測Fig. 6 Prediction of cumulative displacement of monitoring points

4 結(jié) 論

(1)本文提出了一種基于時間序列分析的結(jié)合人工蜂群算法優(yōu)化支持向量機(ABC-SVR)的滑坡位移預測模型。用移動平均法從滑坡累計位移監(jiān)測數(shù)據(jù)中提取滑坡趨勢項位移，用多項式最小二乘法對趨勢項位移進行分段擬合和預測。該模型的提出解決了灰色系統(tǒng)誤差較大問題，避免了神經(jīng)網(wǎng)絡存在的模型結(jié)構(gòu)難確定、精度難保證及泛化能力不高的缺點，為非線性滑坡預測模型在尋優(yōu)過程中易陷入局部最優(yōu)的情況提供了新的思路。

(2)根據(jù)灰色系統(tǒng)關聯(lián)分析法，分析7項影響因子與滑坡周期性位移數(shù)據(jù)的關聯(lián)度。結(jié)果顯示，長江庫水位、前兩個月累積降雨和前1個月累積降雨為相關性rk最大的3項，分別達到(ZG93監(jiān)測點)0.8525、0.8203、0.7987，表明庫水位變動在滑坡位移中的重要性，并且降雨因素對滑坡位移影響具有一定的滯后性。

(3)將人工蜂群算法優(yōu)化支持向量機(ABC-SVR)模型所得結(jié)果同BP神經(jīng)網(wǎng)絡及PSO-SVR方法的結(jié)果進行對比，ABC-SVR模型展現(xiàn)出較好的位移預測結(jié)果。表明基于時間序列的人工蜂群支持向量機位移預測模型可以較好地對滑坡位移進行預測研究，有較好的實用性及應用前景。