張衛(wèi)，徐均，李奇，賈婧，劉紅義，岳鐵鑄

（北京京航計(jì)算通訊研究所，北京100074）

引 言

在強(qiáng)降雨的誘發(fā)下，滑坡、崩塌、泥石流、塌陷等地質(zhì)災(zāi)害經(jīng)常發(fā)生，資料表明，由暴雨誘發(fā)的地質(zhì)災(zāi)害的概率約占85%，其中近95%的滑坡等地質(zhì)災(zāi)害發(fā)生在暴雨開(kāi)始的10小時(shí)以后，86%的滑坡等地質(zhì)災(zāi)害發(fā)生在強(qiáng)降雨開(kāi)始的24小時(shí)以后［1］。無(wú)論在暴雨過(guò)后多長(zhǎng)時(shí)間發(fā)生地質(zhì)災(zāi)害，往往和監(jiān)測(cè)區(qū)域的有效降水量、河道水流情況、地質(zhì)條件等密切相關(guān)。由于滑坡等地質(zhì)災(zāi)害具有突發(fā)性、隨機(jī)性，以及短時(shí)間內(nèi)能造成巨大損失的特點(diǎn)，傳統(tǒng)模型并不能實(shí)時(shí)準(zhǔn)確地預(yù)測(cè)監(jiān)測(cè)區(qū)域的危險(xiǎn)系數(shù)。本系統(tǒng)能根據(jù)降雨型滑坡形成機(jī)理，將無(wú)線傳感網(wǎng)、物聯(lián)網(wǎng)等技術(shù)應(yīng)用到滑坡等地質(zhì)災(zāi)害預(yù)測(cè)模型中，從而全面地感知監(jiān)測(cè)區(qū)域的多維數(shù)據(jù)，實(shí)現(xiàn)動(dòng)態(tài)預(yù)警的同時(shí)提高了預(yù)警準(zhǔn)確性。

1 系統(tǒng)總體設(shè)計(jì)

針對(duì)山區(qū)復(fù)雜地形，結(jié)合物聯(lián)網(wǎng)、無(wú)線傳感網(wǎng)等技術(shù)的應(yīng)用特點(diǎn)，特提出基于物聯(lián)網(wǎng)的滑坡地質(zhì)災(zāi)害預(yù)警監(jiān)控系統(tǒng)，其系統(tǒng)結(jié)構(gòu)如圖1所示。該系統(tǒng)主要由前端無(wú)線感知網(wǎng)絡(luò)、中間傳輸網(wǎng)絡(luò)、數(shù)據(jù)處理中心等組成。其中前端感知網(wǎng)絡(luò)負(fù)責(zé)實(shí)時(shí)感知監(jiān)測(cè)區(qū)域的動(dòng)態(tài)變化，為合理決策提供科學(xué)依據(jù)；中間傳輸網(wǎng)絡(luò)部分主要包括GPRS無(wú)線傳輸和IP網(wǎng)絡(luò)通信等，負(fù)責(zé)把前端無(wú)線感知網(wǎng)絡(luò)采集的數(shù)據(jù)包傳送至遠(yuǎn)端的數(shù)據(jù)處理中心；數(shù)據(jù)處理中心具有數(shù)據(jù)分析統(tǒng)計(jì)、數(shù)據(jù)建模、模糊判斷、數(shù)據(jù)共享等功能。同時(shí)，系統(tǒng)可以通過(guò)網(wǎng)站信息、應(yīng)急短信、電話、LED 顯示屏發(fā)布相關(guān)預(yù)警信息。

圖1 預(yù)警系統(tǒng)結(jié)構(gòu)圖

根據(jù)資料分析與專家意見(jiàn)，將滑坡等地質(zhì)災(zāi)害的影響因子主要?jiǎng)澐譃檎T發(fā)因子與地質(zhì)條件兩大類，具體的影響因子如表1所列［2］。其中，外在誘發(fā)因子主要指有效降水量和水流情況，而內(nèi)在地質(zhì)條件主要包括斷層裂縫、靠近水域、地質(zhì)巖性、地貌起伏度、高程等。根據(jù)影響因子與相關(guān)技術(shù)建立有效的預(yù)警模型，實(shí)現(xiàn)對(duì)監(jiān)測(cè)區(qū)域的滑坡等地質(zhì)災(zāi)害的實(shí)時(shí)預(yù)警。

表1 滑坡等地質(zhì)災(zāi)害影響因子

2 網(wǎng)絡(luò)節(jié)點(diǎn)設(shè)計(jì)

2.1 利用ZigBee技術(shù)互聯(lián)通信

某一監(jiān)測(cè)點(diǎn)內(nèi)部采用ZigBee技術(shù)進(jìn)行組網(wǎng)通信，實(shí)現(xiàn)對(duì)監(jiān)測(cè)區(qū)域的全方面感知與控制。ZigBee WiFi、Bluetooth、GPRS／GSM 等常用無(wú)線通信技術(shù)的具體參數(shù)如表2所列。ZigBee具有成本低、體積小、功耗低、易于擴(kuò)展、感應(yīng)性強(qiáng)等特點(diǎn)［3］。考慮到降雨型滑坡系統(tǒng)的監(jiān)測(cè)網(wǎng)絡(luò)節(jié)點(diǎn)較多、數(shù)據(jù)傳輸量不大，但覆蓋面相對(duì)較廣的特點(diǎn)，感知區(qū)域適合采用ZigBee技術(shù)進(jìn)行互聯(lián)通信。

表2 常用無(wú)線通信技術(shù)參數(shù)

2.2 網(wǎng)絡(luò)節(jié)點(diǎn)的硬件結(jié)構(gòu)設(shè)計(jì)

網(wǎng)絡(luò)節(jié)點(diǎn)是一個(gè)微型嵌入式系統(tǒng)，是構(gòu)成數(shù)據(jù)采集與信息傳輸系統(tǒng)的基本單元，它能在較小的體積內(nèi)集成多種功能，如數(shù)據(jù)采集、信息處理和無(wú)線數(shù)據(jù)收發(fā)等功能。網(wǎng)絡(luò)節(jié)點(diǎn)一般由處理器模塊、ZigBee通信模塊、能量供應(yīng)模塊、傳感器模塊（協(xié)調(diào)節(jié)點(diǎn)可沒(méi)有傳感器模塊）等模塊構(gòu)成。其中傳感器模塊作為整個(gè)系統(tǒng)的“感知器官”，負(fù)責(zé)監(jiān)測(cè)區(qū)域內(nèi)信息的采集和數(shù)據(jù)轉(zhuǎn)換；處理器模塊作為整個(gè)節(jié)點(diǎn)的“大腦”，主要負(fù)責(zé)對(duì)傳感器模塊、ZigBee通信模塊的控制及整個(gè)節(jié)點(diǎn)的協(xié)調(diào)等；ZigBee通信模塊作為網(wǎng)絡(luò)節(jié)點(diǎn)的“耳朵＋嘴巴”，負(fù)責(zé)與其他網(wǎng)絡(luò)節(jié)點(diǎn)進(jìn)行無(wú)線數(shù)據(jù)通信。其中終端節(jié)點(diǎn)僅包含Zigbee通信模塊，協(xié)調(diào)器節(jié)點(diǎn)包含ZigBee通信模塊和GPRS數(shù)據(jù)收發(fā)模塊；能量供應(yīng)模塊則是節(jié)點(diǎn)的“心臟”，為網(wǎng)絡(luò)節(jié)點(diǎn)的各個(gè)組成部分提供足夠的動(dòng)力和能量。

根據(jù)實(shí)際功能需求及部署特點(diǎn)，采用星型的無(wú)線網(wǎng)絡(luò)分布模型，具有組網(wǎng)簡(jiǎn)單快捷、數(shù)據(jù)傳輸快等特點(diǎn)。在該模型下只需設(shè)計(jì)兩類節(jié)點(diǎn)：終端節(jié)點(diǎn)和協(xié)調(diào)器節(jié)點(diǎn)，其結(jié)構(gòu)分別如圖2和圖3所示。其中，協(xié)調(diào)器節(jié)點(diǎn)負(fù)責(zé)收集小范圍內(nèi)終端節(jié)點(diǎn)所采集的數(shù)據(jù)信息，而終端節(jié)點(diǎn)主要負(fù)責(zé)信息采集并實(shí)時(shí)發(fā)送至協(xié)調(diào)器節(jié)點(diǎn)。在網(wǎng)絡(luò)節(jié)點(diǎn)的能量供應(yīng)設(shè)計(jì)方面，由于協(xié)調(diào)器節(jié)點(diǎn)既要負(fù)責(zé)收集整個(gè)Zig-Bee網(wǎng)絡(luò)的信息，又要將數(shù)據(jù)包通過(guò)GPRS 模塊轉(zhuǎn)發(fā)出去，需要較大的功率，可采用太陽(yáng)能或線路供電等模式。而終端節(jié)點(diǎn)只需要定時(shí)將采集的數(shù)據(jù)包發(fā)送至協(xié)調(diào)器節(jié)點(diǎn)，然后就迅速進(jìn)入睡眠狀態(tài)，其功率較小，可采用電池單獨(dú)供電，或者通過(guò)線路與協(xié)調(diào)器節(jié)點(diǎn)共用供電系統(tǒng)。

圖2 終端節(jié)點(diǎn)結(jié)構(gòu)圖

圖3 協(xié)調(diào)器節(jié)點(diǎn)結(jié)構(gòu)圖

在網(wǎng)絡(luò)節(jié)點(diǎn)的操作系統(tǒng)選擇方面，考慮到節(jié)點(diǎn)的任務(wù)量相對(duì)少，而數(shù)據(jù)傳輸?shù)膶?shí)時(shí)性要求較高，系統(tǒng)可選用一些實(shí)時(shí)性較強(qiáng)、可移植、可固化、可裁剪、搶先式多任務(wù)內(nèi)核的操作系統(tǒng)，如Tiny OS、μC／OS-II等。

2.3 傳感器的全面感知

在監(jiān)測(cè)區(qū)域選擇性部署電子自動(dòng)雨量計(jì)，并在已發(fā)生（潛在發(fā)生）滑坡、泥石流所經(jīng)河道的水壩前設(shè)置多個(gè)孔洞，每個(gè)孔洞下端部署一個(gè)液位傳感器，在不同深度部署數(shù)個(gè)液體流速傳感器，實(shí)時(shí)采集監(jiān)測(cè)區(qū)域降雨及河道水流信息。其中，測(cè)量的有效降雨量作為山體滑坡危險(xiǎn)度的第一指標(biāo)，河道水位深度及流速作為輔助指標(biāo)。同時(shí)，在易發(fā)斷層裂縫處部署多個(gè)縫距傳感設(shè)備，測(cè)量地表裂縫分割體之間相對(duì)張開(kāi)、閉合、位錯(cuò)及垂直向升降的變化量，從而全面感知監(jiān)測(cè)區(qū)域的動(dòng)態(tài)變化。

3 預(yù)警模型

3.1 影響因子分析

在地形較為復(fù)雜的山區(qū)，降雨型滑坡等地質(zhì)災(zāi)害一般由一系列外在誘發(fā)因子和內(nèi)在地質(zhì)條件等影響因子共同作用誘發(fā)產(chǎn)生。由于各因子在誘發(fā)災(zāi)害過(guò)程中的作用程度不同，內(nèi)在則表現(xiàn)為各影響因子在數(shù)據(jù)模型中所占權(quán)重不同。由于很難獲取大量的歷史滑坡統(tǒng)計(jì)數(shù)據(jù)，本文結(jié)合相關(guān)資料，相關(guān)影響因子權(quán)值采用主觀賦權(quán)法（也稱為專家賦權(quán)法）來(lái)確定，即通過(guò)一定方法綜合各位專家對(duì)各指標(biāo)給出的權(quán)重進(jìn)行的賦權(quán)。經(jīng)統(tǒng)計(jì)，各影響因子的權(quán)重調(diào)整范圍及初始值分配如表3所列。

表3 各級(jí)影響因子權(quán)重分配

（1）影響因子

研究表明降雨對(duì)山體滑坡的誘發(fā)作用，不僅取決于當(dāng)日降雨量，還和近期降水量、河道水流等情況密切相關(guān)。系統(tǒng)利用部署的電子雨量計(jì)實(shí)時(shí)測(cè)量降水量r，在河道中安裝液位傳感器，測(cè)量數(shù)據(jù)為h。近期滑坡等地質(zhì)災(zāi)害有效降水量Re是指近期各日降雨量（M 天內(nèi)）與其影響系數(shù)乘積之和，而有效估計(jì)降水量R為有效降水量Re、近期（次日）預(yù)計(jì)降水量Rf和各自影響系數(shù)乘積之和，具體計(jì)算公式為：

其中i表示雨期距今天數(shù)，i＝0表示當(dāng)天，i＝1表示前一天，Ri為某天降水量（當(dāng)天降水量可以分時(shí)統(tǒng)計(jì)）；αi為該天的降水影響系數(shù)，考慮到持續(xù)型降雨誘發(fā)地質(zhì)災(zāi)害的滯后性，隨著計(jì)算雨量日期的前移，系數(shù)αi逐漸減??；β1和β2分別為有效降水量、近期預(yù)計(jì)降水量的影響因子。結(jié)合近期發(fā)生泥石流時(shí)的平均警戒有效降水量Rav、河壩的警戒水位高度H，參照公式（3）聯(lián)合估計(jì)有效降雨所產(chǎn)生的不穩(wěn)定分值f。

其中，參數(shù)λ1、λ2、K1、K2均為可調(diào)整參數(shù)，K1、K2初始值為1。

（2）地質(zhì)巖性

根據(jù)當(dāng)?shù)氐牡刭|(zhì)測(cè)繪圖，并結(jié)合實(shí)地考察，對(duì)各個(gè)監(jiān)測(cè)點(diǎn)的地質(zhì)巖性進(jìn)行量化，初步劃分為4 個(gè)等級(jí)，即V1松散巖組、V2粘軟巖、V3中硬巖組、V4堅(jiān)硬巖組［1］。

（3）到水系距離

滑坡等地質(zhì)災(zāi)害與到周邊水系（如湖泊、池塘、水庫(kù)、較大河流等）的有效距離有著比較緊密關(guān)系。據(jù)統(tǒng)計(jì)，距水系的有效距離越近，發(fā)生地質(zhì)災(zāi)害的頻率越高。系統(tǒng)可根據(jù)監(jiān)測(cè)區(qū)域地質(zhì)地圖等資料，計(jì)算監(jiān)測(cè)點(diǎn)距周邊水系的有效距離。

（4）高 程

可以利用監(jiān)測(cè)區(qū)域的GIS地圖，對(duì)前段監(jiān)測(cè)區(qū)域的高度進(jìn)行統(tǒng)計(jì)。根據(jù)高程對(duì)滑坡等地質(zhì)災(zāi)害影響的程度和實(shí)際情況，初步將高程劃分為3個(gè)等級(jí)，即500 m 以下、500～1000m、1000m 以上。

（5）地貌起伏度

地貌起伏度反映地表起伏變化，常指某一確定面積（監(jiān)測(cè)區(qū)域）內(nèi)最高點(diǎn)和最低點(diǎn)海拔高度之差。可以利用GIS等技術(shù)建立數(shù)字高程模型，提取地貌起伏度、坡度和坡向等地形參數(shù)。將監(jiān)測(cè)區(qū)域地貌分為：H1平坦起伏（0～50m）、H2小起伏（50～100m）、H3中起伏（100～300 m）、H4山地起伏（300～600 m）、H5 高山起伏（600～1500m）等。由于極高起伏地區(qū)的人類活動(dòng)稀少，不在統(tǒng)計(jì)范圍內(nèi)。

此外，要在易發(fā)斷層觸發(fā)帶安裝裂縫器以測(cè)量裂縫的情況，全面感知潛在誘發(fā)因子。由于裂縫和地質(zhì)附著力等地質(zhì)環(huán)境密切相關(guān)，此部分要結(jié)合實(shí)際情況而定。

3.2 多維權(quán)重組合模型

在組合型數(shù)據(jù)模型中，需要在一定范圍內(nèi)調(diào)整一階因子權(quán)值，使其和為1。任取二級(jí)因子網(wǎng)格單元k，其所屬的一級(jí)因子的權(quán)值為wi，針對(duì)實(shí)際測(cè)量情況，在一級(jí)因子確定的情況下，量化其二級(jí)因子，確定歸于或者相近于哪一具體子集，從而確定二級(jí)權(quán)值，即確定j值。例如，在一級(jí)影響因子地質(zhì)巖性的范圍內(nèi)，經(jīng)過(guò)查閱GIS等勘測(cè)資料和實(shí)地考察等方式，確定最接近于二級(jí)子集v1、v2、v3、v4中的哪一種，從而選擇合適的二級(jí)權(quán)值。對(duì)于組合了各二級(jí)因子的網(wǎng)格單元k，定義其不穩(wěn)定分值Si＝wi·，然后根據(jù)各級(jí)不穩(wěn)定分值確定滑坡等地質(zhì)災(zāi)害的易發(fā)性指數(shù)（LandsLide Susceptibility Index，LSI）。

LSI是評(píng)估監(jiān)測(cè)區(qū)域的最敏感性指標(biāo)，值越大，發(fā)生滑坡等地質(zhì)災(zāi)害的概率越大，當(dāng)其超過(guò)某具體閥值，就會(huì)提前預(yù)警。

3.3 數(shù)據(jù)分析

按照上述公式（1）～（4），計(jì)算各個(gè)監(jiān)測(cè)點(diǎn)的LSI，根據(jù)地質(zhì)災(zāi)害的易發(fā)性指數(shù)將檢測(cè)區(qū)域劃分為極輕微、輕微、中等、較嚴(yán)重、嚴(yán)重、極嚴(yán)重6級(jí)，并用不同顏色進(jìn)行標(biāo)記，然后采用等間距、均值-標(biāo)準(zhǔn)差等方法，實(shí)時(shí)（間隔固定時(shí)間，如半小時(shí)繪制一次）繪制整個(gè)監(jiān)測(cè)區(qū)域的滑坡等地質(zhì)災(zāi)害危險(xiǎn)系數(shù)預(yù)警圖，從而實(shí)現(xiàn)滑坡等地質(zhì)災(zāi)害預(yù)報(bào)的精細(xì)化管理。

在數(shù)據(jù)建模方面，采用了自調(diào)整的動(dòng)態(tài)神經(jīng)網(wǎng)絡(luò)模型，利用歷史數(shù)據(jù)或者實(shí)施過(guò)程中逐步獲取的數(shù)據(jù)，對(duì)模型的自身結(jié)構(gòu)及學(xué)習(xí)規(guī)則進(jìn)行了動(dòng)態(tài)優(yōu)化，調(diào)整各級(jí)權(quán)重，使模型達(dá)到最優(yōu)，從而使得系統(tǒng)具有更強(qiáng)的適應(yīng)性和準(zhǔn)確性。

結(jié) 語(yǔ)

本文在原有模型的基礎(chǔ)上進(jìn)行改進(jìn)，利用傳感器網(wǎng)絡(luò)實(shí)時(shí)感知監(jiān)測(cè)區(qū)域的多種指標(biāo)，運(yùn)用多因素動(dòng)態(tài)聯(lián)合建模，采用定量分析法，對(duì)降雨型滑坡等地質(zhì)災(zāi)害進(jìn)行相對(duì)準(zhǔn)確預(yù)報(bào)。系統(tǒng)根據(jù)每個(gè)監(jiān)測(cè)點(diǎn)的LSI，利用相關(guān)方法繪制整個(gè)監(jiān)測(cè)區(qū)域的危險(xiǎn)系數(shù)預(yù)警圖。對(duì)于超過(guò)閾值的，立即報(bào)警，有效預(yù)防和避免了地質(zhì)災(zāi)害，保護(hù)人民生命財(cái)產(chǎn)安全。

系統(tǒng)具有高集成度、高精度、全天候自動(dòng)化實(shí)時(shí)監(jiān)測(cè)與自動(dòng)預(yù)警的能力，同時(shí)還具有監(jiān)測(cè)范圍大、部署靈活的特點(diǎn)。在指導(dǎo)降雨型滑坡的防災(zāi)與減災(zāi)方面具有較強(qiáng)的實(shí)用價(jià)值。

［1］王衛(wèi)東.基于GIS的區(qū)域公路地質(zhì)災(zāi)害管理與空間決策支持系統(tǒng)研究［D］.長(zhǎng)沙：中南大學(xué)商學(xué)院，2009.

［2］解傳銀.基于權(quán)重模型的滑坡災(zāi)害敏感性評(píng)價(jià)［J］.中南大學(xué)學(xué)報(bào)：自然科學(xué)版，2011，42（6）：1772-1779.

［3］Sena.The comparison of Wi-Fi，Bluetooth and ZigBee［EB／OL］.［2012-11-10］.http：／／www.sena.com／blog／？p＝359.

［4］Brew G，Barazangi M.Tectonic map and geologic evolution of Syria：The role of GIS［J］.The Leading Edge，2000（19）：176-182.

［5］郭芳芳，楊農(nóng).地形起伏度和坡度分析在區(qū)域滑坡災(zāi)害評(píng)價(jià)中的應(yīng)用［J］.中國(guó)地質(zhì)，2008，35（1）：131-141.

［6］國(guó)土資源部.四川雅安地質(zhì)災(zāi)害預(yù)警示范區(qū)［EB／OL］.［2012-11-10］.http：／／www.mlr.gov.cn／dzhj／201003／t20100326＿142901.htm.2012.

［7］張春山，張業(yè)成，馬寅生.黃河上游地區(qū)崩塌、滑坡、泥石流地質(zhì)災(zāi)害區(qū)域危險(xiǎn)性評(píng)價(jià)［J］.地質(zhì)力學(xué)學(xué)報(bào)，2003，9（2）：143-153.

［8］殷坤龍，朱良峰.滑坡災(zāi)害空間區(qū)劃及GIS 應(yīng)用研究［J］.地學(xué)前緣，2001，8（2）：279-284.

［9］劉藝梁，殷坤龍.邏輯回歸和人工神經(jīng)網(wǎng)絡(luò)模型在滑坡災(zāi)害空間預(yù)測(cè)中的應(yīng)用［J］.水文地質(zhì)工程地質(zhì)，2010，37（5）：92-96.

［10］曹炯清.基于IEEE 802.15.4無(wú)線傳感器網(wǎng)絡(luò)的山體滑坡檢測(cè)系統(tǒng)［J］.電腦知識(shí)與技術(shù)，2009，5（8）：1817-1819.

［11］張繼賢.3S支持下的滑坡地質(zhì)災(zāi)害監(jiān)測(cè)評(píng)估與建模［J］.測(cè)繪工程，2005，14（2），1-5.

［12］謝洪，鐘敦倫，矯震，等.2008年汶川地震重災(zāi)區(qū)的泥石流［J］.山地學(xué)報(bào)，2008，27（4）：501 509.