黃 玲，王 霄，張 譯

（貴州大學(xué)電氣工程學(xué)院，貴陽(yáng) 550025）

1 引言

我國(guó)疆域遼闊，地質(zhì)環(huán)境復(fù)雜多樣[1]，地質(zhì)災(zāi)害分布廣、類型多、頻度高、強(qiáng)度大[2]，滑坡、崩塌及泥石流等已經(jīng)成為危害最大的地質(zhì)災(zāi)害[3]，一旦發(fā)生，常給人民的生命財(cái)產(chǎn)安全和工業(yè)生產(chǎn)活動(dòng)造成巨大的損失。針對(duì)礦山開(kāi)采誘發(fā)滑坡的廣泛性、邊坡巖體復(fù)雜性等實(shí)際情況[4]，國(guó)內(nèi)外學(xué)者和技術(shù)人員對(duì)滑坡監(jiān)測(cè)系統(tǒng)開(kāi)發(fā)進(jìn)行了深入研究，許多文獻(xiàn)對(duì)滑坡災(zāi)害監(jiān)測(cè)及預(yù)警預(yù)報(bào)技術(shù)進(jìn)行了研究和探討[5]，但仍然存在許多不足，如監(jiān)測(cè)設(shè)備不夠智能、預(yù)警模型無(wú)法滿足需求等，使專業(yè)監(jiān)測(cè)預(yù)警的作用沒(méi)有完全發(fā)揮出來(lái)。

據(jù)此，從貴州地質(zhì)構(gòu)造的實(shí)踐需要出發(fā)，設(shè)計(jì)并實(shí)現(xiàn)一套適合貴州地質(zhì)特色的集數(shù)據(jù)存儲(chǔ)、展示、分析山體滑坡穩(wěn)定性為一體的綜合山體監(jiān)測(cè)預(yù)警系統(tǒng)。通過(guò)數(shù)學(xué)建模預(yù)測(cè)山體狀態(tài)，進(jìn)行危險(xiǎn)等級(jí)的劃分，再利用數(shù)據(jù)可視化技術(shù)在后端對(duì)所設(shè)計(jì)的數(shù)據(jù)庫(kù)中的數(shù)據(jù)進(jìn)行可視化處理并傳輸?shù)絎eb端[6-7]進(jìn)行展示，通過(guò)訪問(wèn)Web端，實(shí)時(shí)查看山體滑坡的相關(guān)數(shù)據(jù)，從而對(duì)山體滑坡的歷史數(shù)據(jù)進(jìn)行分析評(píng)價(jià)，利用存儲(chǔ)的歷史工作數(shù)據(jù)對(duì)山體滑坡進(jìn)行更好的預(yù)測(cè)?；赪eb的可視化監(jiān)測(cè)預(yù)警系統(tǒng)可實(shí)現(xiàn)數(shù)據(jù)監(jiān)測(cè)、可視化展示以及對(duì)山體滑坡進(jìn)行預(yù)測(cè)等功能，對(duì)避免由山體滑坡帶來(lái)的巨大財(cái)產(chǎn)損失和人傷亡極具現(xiàn)實(shí)意義[8]。

2 系統(tǒng)框架設(shè)計(jì)

系統(tǒng)通過(guò)切向角傳感器、降雨量傳感器、土壤溫濕度傳感器、環(huán)境溫濕度傳感器將采集到的數(shù)據(jù)通過(guò)MQTT傳輸數(shù)據(jù)，再利用阿里云服務(wù)器作為中轉(zhuǎn)站，負(fù)責(zé)下位機(jī)網(wǎng)絡(luò)交互和應(yīng)用信息發(fā)布。

在Python 3.7環(huán)境下配置MQTT通信協(xié)議，實(shí)現(xiàn)數(shù)據(jù)的訂閱、處理與數(shù)據(jù)可視化。其中利用Django框架開(kāi)發(fā)監(jiān)控界面；使用Python語(yǔ)言實(shí)現(xiàn)系統(tǒng)后臺(tái)功能邏輯和滑坡預(yù)測(cè)；使用HTML語(yǔ)言和相關(guān)腳本語(yǔ)言實(shí)現(xiàn)前端模板頁(yè)面；選擇Mysql保存山體滑坡等數(shù)據(jù)并使用Navicat工具對(duì)數(shù)據(jù)庫(kù)進(jìn)行管理，隨后利用數(shù)據(jù)可視化技術(shù)在后端對(duì)數(shù)據(jù)庫(kù)中數(shù)據(jù)進(jìn)行可視化處理并傳輸?shù)絎eb端進(jìn)行展示。

同時(shí)，在系統(tǒng)中融入模糊綜合評(píng)價(jià)算法實(shí)現(xiàn)對(duì)山體滑坡的危險(xiǎn)等級(jí)劃分并對(duì)山體滑坡進(jìn)行預(yù)測(cè)，最后將整個(gè)監(jiān)測(cè)預(yù)警系統(tǒng)進(jìn)行服務(wù)器部署以適應(yīng)實(shí)際工業(yè)環(huán)境的應(yīng)用。整體系統(tǒng)框架如圖1所示。

圖1 山體滑坡預(yù)警系統(tǒng)框架圖

3 層次分析模糊綜合評(píng)價(jià)法

滑坡是非線性系統(tǒng)，等級(jí)狀態(tài)往往具有模糊性和不確定性，模糊綜合評(píng)價(jià)法常作非線性判斷滑坡穩(wěn)定性方法之一。其模型簡(jiǎn)單、易于掌握，因此模糊綜合評(píng)價(jià)法被廣泛用于滑坡的預(yù)測(cè)預(yù)報(bào)中，但仍存在滑坡預(yù)警等級(jí)無(wú)法進(jìn)行精確劃分和評(píng)價(jià)指標(biāo)權(quán)重賦值問(wèn)題?；诖耍诖瞬捎媚：C合評(píng)價(jià)法對(duì)山體滑坡等級(jí)進(jìn)行劃分，利用該方法對(duì)用環(huán)境溫濕度傳感器、土壤濕度傳感器、降雨量傳感器、傾角位移傳感器等采集到的數(shù)據(jù)進(jìn)行融合，從而得到山體滑坡預(yù)警等級(jí)的評(píng)價(jià)結(jié)果。該方法以數(shù)學(xué)為基礎(chǔ)，構(gòu)造出評(píng)價(jià)指標(biāo)體系并根據(jù)層次分析法確定各個(gè)因素的權(quán)重，通過(guò)模糊變換原理，確認(rèn)各因素隸屬度函數(shù)，將構(gòu)造得到的模糊矩陣確定山體滑坡的穩(wěn)定性。方法發(fā)揮了專家經(jīng)驗(yàn)的優(yōu)點(diǎn)，減少人為臆斷導(dǎo)致的誤差，其評(píng)價(jià)結(jié)果具有客觀準(zhǔn)確性。算法的流程圖如圖2所示。

圖2 層次分析模糊綜合評(píng)價(jià)算法流程圖

3.1 建立預(yù)警模型因素和評(píng)價(jià)集

因素集是由影響評(píng)價(jià)對(duì)象因素構(gòu)成的集合，通過(guò)查閱國(guó)內(nèi)外相關(guān)文獻(xiàn)和專家經(jīng)驗(yàn)，此山體滑坡評(píng)價(jià)模型選用環(huán)境溫濕度、土壤濕度、降雨量、位移量作為影響因素，確定山體滑坡監(jiān)測(cè)預(yù)警模型的因素集如下式：

式中，ui（i=1,2,3,4,5）分別對(duì)應(yīng)環(huán)境溫度、環(huán)境濕度、土壤濕度、降雨量和位移量，它們具有不同程度的模糊性。

該山體滑坡監(jiān)測(cè)預(yù)警模型評(píng)價(jià)集為對(duì)山體滑坡穩(wěn)定性評(píng)價(jià)做出的評(píng)價(jià)結(jié)果的集合，劃分為3個(gè)等級(jí)，記為：

式中，vi（i=1,2,3）為山體滑坡監(jiān)測(cè)預(yù)警模型的評(píng)價(jià)指標(biāo)，分別表示安全級(jí)、注意級(jí)和警報(bào)級(jí)，詳細(xì)參數(shù)指標(biāo)如表1所示。

表1 預(yù)警影響因子等級(jí)評(píng)價(jià)指標(biāo)

3.2 隸屬度函數(shù)和因素權(quán)重集

評(píng)價(jià)指標(biāo)一般是用模糊語(yǔ)言來(lái)描述的，使用隸屬函數(shù)數(shù)量化。為消除其模糊性，需構(gòu)建隸屬函數(shù)進(jìn)行數(shù)量化評(píng)價(jià)，因素集到評(píng)語(yǔ)集的模糊映射即為隸屬程度。用(0,1)的數(shù)值表示隸屬程度，其值越高，即隸屬程度越高，模糊性就越高。使用線性關(guān)系得到隸屬函數(shù)值。由此得到單滑坡影響因素的綜合模糊關(guān)系矩陣R如下式：

采用層次分析法確定山體滑坡各因素權(quán)重，建立遞階層次結(jié)構(gòu)模型，將其分為三層，分別為目標(biāo)層、準(zhǔn)則層和方案層。目標(biāo)層需對(duì)山體滑坡監(jiān)測(cè)預(yù)警模型做出穩(wěn)定性評(píng)價(jià)；準(zhǔn)則層為各影響因子；方案層為山體滑坡評(píng)價(jià)結(jié)果，將其分為三個(gè)等級(jí)：安全級(jí)、注意集和警報(bào)集。最后構(gòu)造判斷矩陣，計(jì)算權(quán)重并進(jìn)行一致性檢驗(yàn)。

3.3 預(yù)警模糊綜合評(píng)價(jià)模型

將滑坡因素集上的集合通過(guò)模糊變換為滑坡評(píng)價(jià)集上的模糊集合，即為山體滑坡監(jiān)測(cè)預(yù)警模糊綜合評(píng)價(jià)，其變換式如下式：

上式中，C為滑坡評(píng)價(jià)集上的三個(gè)等級(jí)模糊子集；cj（j=1,2,3）為相應(yīng)等級(jí)子集的隸屬度，cj由ai和ri共同決定。比較各個(gè)影響因素的隸屬度大小情況，再根據(jù)最大隸屬度原則，輸出結(jié)果則為隸屬度最大所對(duì)應(yīng)的那個(gè)滑坡監(jiān)測(cè)預(yù)警等級(jí)。

4 軟件架構(gòu)及模塊設(shè)計(jì)

整個(gè)山體滑坡監(jiān)測(cè)預(yù)警系統(tǒng)基于Web服務(wù)，采用B/S工作方式，通過(guò)瀏覽器訪問(wèn)山體滑坡監(jiān)測(cè)預(yù)警系統(tǒng)。用戶通過(guò)瀏覽器連接到該山體滑坡監(jiān)測(cè)預(yù)警網(wǎng)站時(shí)，處于邏輯層的Web服務(wù)器從文件系統(tǒng)中加載Web頁(yè)面腳本，并將其傳遞給腳本引擎。腳本引擎負(fù)責(zé)解析并執(zhí)行腳本，腳本使用數(shù)據(jù)庫(kù)連接器打開(kāi)存儲(chǔ)層連接并對(duì)數(shù)據(jù)庫(kù)執(zhí)行Sql語(yǔ)句。此山體滑坡監(jiān)測(cè)預(yù)警系統(tǒng)采用Django框架進(jìn)行Web開(kāi)發(fā)，前端使用HTML語(yǔ)言、CSS、JS等實(shí)現(xiàn)網(wǎng)頁(yè)展示，后端主要使用Django自帶的Web框架進(jìn)行一些邏輯編寫(xiě)，融合模糊綜合評(píng)價(jià)法實(shí)現(xiàn)山體滑坡的數(shù)據(jù)處理、等級(jí)判斷和結(jié)果預(yù)測(cè)等。同時(shí)，使用Ajax技術(shù)實(shí)現(xiàn)前后端交互。

4.1 用戶信息和致災(zāi)因子管理模塊

在用戶信息管理模塊，用戶可填寫(xiě)自己的基本信息并能對(duì)自己的信息進(jìn)行修改。用戶登陸以后可以進(jìn)入自己的個(gè)人中心，查看和修改個(gè)人信息。在滑坡致災(zāi)因子管理模塊中，主要是顯示山體滑坡致災(zāi)因子的一些數(shù)據(jù)。經(jīng)授權(quán)登陸后的用戶可以點(diǎn)擊查看山體滑坡影響因子的一些簡(jiǎn)要介紹以及它們對(duì)滑坡的詳細(xì)影響情況。管理員可以在數(shù)據(jù)庫(kù)中添加或刪除山體滑坡影響因子，用戶則可以在前端界面進(jìn)行影響因子查詢。系統(tǒng)顯示界面如圖3所示。

圖3 滑坡影響因子顯示界面

4.2 山體滑坡穩(wěn)定性評(píng)價(jià)模塊

這一模塊是模糊綜合評(píng)價(jià)算法對(duì)山體滑坡做出穩(wěn)定性評(píng)價(jià)的可視化實(shí)現(xiàn)，既可對(duì)山體滑坡進(jìn)行穩(wěn)定性評(píng)價(jià)，還可查看山體滑坡歷史信息。

用戶登陸成功之后進(jìn)入首頁(yè)，當(dāng)用戶點(diǎn)擊“山體滑坡歷史信息查看”這一欄時(shí)，就會(huì)進(jìn)入另一個(gè)界面，通過(guò)MQTT將部署在現(xiàn)場(chǎng)的傳感器采集到的數(shù)據(jù)實(shí)時(shí)傳來(lái)，顯示在此部分，用戶即可實(shí)時(shí)監(jiān)測(cè)當(dāng)?shù)鼗碌木唧w情況，歷史數(shù)據(jù)也會(huì)保留。若點(diǎn)擊山體滑坡穩(wěn)定性可進(jìn)入另一個(gè)界面，如圖4所示。

圖4 滑坡穩(wěn)定性評(píng)價(jià)界面

為讓用戶更直觀地了解當(dāng)?shù)鼗碌攸c(diǎn)及穩(wěn)定性情況，本系統(tǒng)以發(fā)布地圖的形式讓用戶查看玉屏縣的滑坡數(shù)據(jù)以及穩(wěn)定情況。通過(guò)天地圖，調(diào)用API接口，制作出玉屏縣地圖。這里以地圖或者衛(wèi)星的形式進(jìn)行發(fā)布，將玉屏縣的地圖數(shù)據(jù)、滑坡點(diǎn)數(shù)據(jù)發(fā)布成地圖服務(wù)形式載入系統(tǒng)之中，融合算法模型，點(diǎn)擊地圖上的滑坡點(diǎn)，就會(huì)彈出該滑坡點(diǎn)的影響因子情況以及該點(diǎn)的滑坡穩(wěn)定性情況。

5 系統(tǒng)測(cè)試與應(yīng)用

為讓該系統(tǒng)更好為當(dāng)?shù)鼐用穹?wù)，設(shè)計(jì)選取一個(gè)臨災(zāi)階段來(lái)作為研究現(xiàn)場(chǎng)來(lái)進(jìn)行測(cè)試，將硬件部署在所選取的研究區(qū)域?qū)崿F(xiàn)實(shí)時(shí)監(jiān)測(cè)。將系統(tǒng)提前部署在貴州省玉屏縣木弄村的野雞坡水庫(kù)。通過(guò)網(wǎng)頁(yè)訪問(wèn)可得到由傳感器傳來(lái)的數(shù)據(jù)：環(huán)境濕度72.6%、環(huán)境溫度26.3℃、土壤濕度76.2%、切線角α度、降雨量66 mm/h，滑坡發(fā)生時(shí)間為2020年7月9日的13:23，網(wǎng)頁(yè)中發(fā)布預(yù)警時(shí)間為同年7月9日12:01。后經(jīng)實(shí)地勘察，確實(shí)發(fā)現(xiàn)該地點(diǎn)出現(xiàn)塌方現(xiàn)象，后聯(lián)系當(dāng)?shù)叵嚓P(guān)部門到事發(fā)地點(diǎn)進(jìn)行相應(yīng)處理?，F(xiàn)場(chǎng)清理圖如圖5所示。

圖5 滑坡地點(diǎn)現(xiàn)場(chǎng)清理圖

6 結(jié)束語(yǔ)

基于Web技術(shù)，根據(jù)系統(tǒng)設(shè)計(jì)標(biāo)準(zhǔn)流程設(shè)計(jì)并實(shí)現(xiàn)了山體滑坡的監(jiān)測(cè)和預(yù)警功能。用戶通過(guò)簡(jiǎn)單的網(wǎng)站訪問(wèn)即可實(shí)現(xiàn)對(duì)山體滑坡數(shù)據(jù)的實(shí)時(shí)監(jiān)測(cè)并能接收滑坡的預(yù)警信息，從而提前做好預(yù)防措施。有效實(shí)現(xiàn)山體滑坡的遠(yuǎn)程監(jiān)測(cè)改變了以往的村頭播報(bào)滑坡信息方式，提高了遠(yuǎn)程監(jiān)測(cè)管理的效率和居民應(yīng)對(duì)災(zāi)害的能力。將該系統(tǒng)應(yīng)用到實(shí)際案例中，結(jié)果顯示該系統(tǒng)的各個(gè)功能模塊正常運(yùn)行，能夠?qū)逻M(jìn)行有效的監(jiān)測(cè)預(yù)警。