朱建國(guó)，叢沛桐，劉 浩

(1.華南農(nóng)業(yè)大學(xué) 水利與土木工程學(xué)院，廣東 廣州 510642；2.廣東省農(nóng)業(yè)水土信息無(wú)人機(jī)采集遙感工程技術(shù)研究中心，廣東 廣州 510642)

1 概述

鐘山縣處于我國(guó)西南地區(qū)，行政區(qū)劃屬?gòu)V西東部偏北的賀州市，北緯N24°31′45.35″，東經(jīng)E111°17′51.26″，北回歸線偏北區(qū)域，氣候?qū)儆趤啛釒Ъ撅L(fēng)氣候是典型的亞熱帶過(guò)渡帶。鐘山縣屬于橫斷山脈邊緣區(qū)的丘陵山地地區(qū)的富鐘賀巖溶盆地，境內(nèi)徑流量充沛，水資源總量豐富，但常發(fā)生局部性旱澇災(zāi)害。

國(guó)外的Varnes、Hansell、Choi等學(xué)者分別提出山洪地質(zhì)災(zāi)害三原理；預(yù)測(cè)了20世紀(jì)80年代2月中旬的暴風(fēng)雨期間發(fā)生的重大滑坡事件的發(fā)生時(shí)間；通過(guò)頻率比模型，利用滑坡發(fā)生因子，分析了滑坡的敏感性及韓國(guó)張黃地區(qū)的滑坡災(zāi)害[1-2]。國(guó)內(nèi)的學(xué)者殷坤龍、劉傳正、賈茜淳、王晶等人也分別建立了浙江省與網(wǎng)絡(luò)連接的預(yù)警預(yù)報(bào)的防災(zāi)減災(zāi)體系，應(yīng)用到地質(zhì)災(zāi)害氣象預(yù)報(bào)預(yù)警實(shí)際工作中，建立廣義回歸神經(jīng)網(wǎng)絡(luò)系統(tǒng)，預(yù)測(cè)出不同降雨量下，不同地區(qū)的山洪地質(zhì)災(zāi)害危險(xiǎn)等級(jí)；引入土壤松散程度因子，對(duì)廣東省博羅縣瀾石河小流域的山洪地質(zhì)災(zāi)害危險(xiǎn)性和易損性進(jìn)行評(píng)價(jià)[3-7]。

山洪歷時(shí)較短，短時(shí)內(nèi)洪峰流量可達(dá)到500 m3/s以上，我國(guó)平均每年山洪地質(zhì)災(zāi)害造成的傷亡率占洪澇災(zāi)害的比例超過(guò)60%[8-10]。

山洪地質(zhì)災(zāi)害往往伴隨形成泥石流、滑坡等，具有突發(fā)性強(qiáng)、影響范圍廣、破壞能力強(qiáng)等特點(diǎn)，其一旦形成，大概率會(huì)造成人員傷亡，財(cái)產(chǎn)(道橋、農(nóng)田、牲畜、房屋等)受到不同程度的破壞，對(duì)國(guó)民經(jīng)濟(jì)和人民生命財(cái)產(chǎn)造成極大的危害。面對(duì)山洪地質(zhì)災(zāi)害防治的需求，提高數(shù)據(jù)接收、處理和服務(wù)能力，增強(qiáng)降雨帶來(lái)的山洪地質(zhì)災(zāi)害精準(zhǔn)預(yù)警，提高長(zhǎng)期山洪地質(zhì)災(zāi)害短時(shí)預(yù)測(cè)能力，為人民生命財(cái)產(chǎn)安全提供保障。

2 鐘山山洪地質(zhì)災(zāi)害影響因子調(diào)查

對(duì)鐘山縣裸地土壤類型研究采取實(shí)地調(diào)查。鐘山縣地質(zhì)構(gòu)造復(fù)雜多樣，表層有強(qiáng)風(fēng)化巖石，土質(zhì)條件較脆弱松散，地形多山區(qū)丘陵，易發(fā)生山洪地質(zhì)災(zāi)害。鐘山縣大部分為泥盆系地層，花山鄉(xiāng)、兩安鄉(xiāng)多為花崗巖體，鳳翔鎮(zhèn)、公安鎮(zhèn)、紅花鎮(zhèn)、回龍鎮(zhèn)、清塘鎮(zhèn)、珊瑚鎮(zhèn)、石龍鎮(zhèn)、同古鎮(zhèn)、燕塘鎮(zhèn)、鐘山鎮(zhèn)的主要巖性為砂質(zhì)頁(yè)巖和灰?guī)r。鐘山裸露巖體多為砂巖、粉質(zhì)粘土，土體疏松，土壤松散程度大，斜坡穩(wěn)定性差，抗剪強(qiáng)度低，土質(zhì)松散地帶極易發(fā)生滑坡、泥石流、崩塌等山洪地質(zhì)災(zāi)害。地質(zhì)災(zāi)害調(diào)查點(diǎn)數(shù)量統(tǒng)計(jì)見(jiàn)表1。

表1 鐘山縣地質(zhì)災(zāi)害調(diào)查點(diǎn)數(shù)量 個(gè)

在強(qiáng)降雨的作用下，容易使松散土體和水的聚集，導(dǎo)致山體表面土體失穩(wěn)，形成滑坡泥石流等地質(zhì)災(zāi)害。有研究表明不同松散程度的土壤發(fā)生山洪地質(zhì)災(zāi)害的種類也不同，如泥石流、滑坡、崩塌、地面塌陷等。

經(jīng)調(diào)查，近3 a鐘山縣發(fā)生地質(zhì)災(zāi)害155起，累計(jì)傷亡約15人，平均每年造成的直接損失300余萬(wàn)元。根據(jù)災(zāi)害大小，把山洪地質(zhì)災(zāi)害造成損失面積小于0.1 km2的定為低危險(xiǎn)類型，造成損失面積在0.1～1.0 km2的定為一般危險(xiǎn)類型，造成損失面積在1.0～5 km2的定為中等危險(xiǎn)類型，造成損失面積大于5 km2的定為高危險(xiǎn)類型。圖1為4等災(zāi)害類型占比示意，其中高危、中等、一般和低危類型山洪地質(zhì)災(zāi)害占比分別為16%、12%、33%、39%。

面對(duì)山洪地質(zhì)災(zāi)害防治的需求，提高數(shù)據(jù)接收、處理和服務(wù)能力，提高山洪地質(zhì)災(zāi)害短時(shí)預(yù)測(cè)能力，增強(qiáng)降雨帶來(lái)的山洪地質(zhì)災(zāi)害短時(shí)精準(zhǔn)預(yù)警顯得更加緊迫。

調(diào)查表明導(dǎo)致山洪地質(zhì)災(zāi)害的直接誘因是突發(fā)性強(qiáng)降水。廣西鐘山縣地質(zhì)災(zāi)害隱患點(diǎn)多面廣，具有群發(fā)性、突發(fā)性、隱蔽性強(qiáng)、危害巨大等特點(diǎn)，切坡建房引發(fā)的地質(zhì)災(zāi)害隱患極為普遍，地質(zhì)災(zāi)害防治工作依然嚴(yán)峻，防災(zāi)減災(zāi)任務(wù)十分艱巨。因此研發(fā)地質(zhì)災(zāi)害預(yù)警APP，可以有效便捷地把山洪地質(zhì)災(zāi)害風(fēng)險(xiǎn)等級(jí)帶給廣大群眾，定期推送鐘山縣政府的防災(zāi)宣傳內(nèi)容，提高民眾危險(xiǎn)意識(shí)，增強(qiáng)民眾對(duì)災(zāi)害的辨識(shí)、避免因?yàn)?zāi)受傷、提高自救能力。

圖1 4等災(zāi)害類型占比示意

3 技術(shù)路線

收集廣西壯族自治區(qū)鐘山縣各鄉(xiāng)鎮(zhèn)歷史地質(zhì)災(zāi)害發(fā)生區(qū)域的資料，從所收集到的鐘山縣歷史調(diào)查資料進(jìn)一步進(jìn)行深入的分析，探究鐘山縣山洪地質(zhì)災(zāi)害發(fā)生的規(guī)律及特點(diǎn)，找出形成災(zāi)害的關(guān)鍵影響因子。依靠定量降水預(yù)報(bào)、結(jié)合ArcGIS不同圖層的數(shù)據(jù)以及調(diào)查得出的土壤條件資料，基于MATLAB神經(jīng)網(wǎng)絡(luò)算法，對(duì)鐘山縣山洪地質(zhì)災(zāi)害風(fēng)險(xiǎn)進(jìn)評(píng)價(jià)。將MATLAB神經(jīng)網(wǎng)絡(luò)工具箱編譯為可被Java類型，前端基于leaflet開(kāi)源地圖并結(jié)合當(dāng)?shù)貧庀箢A(yù)報(bào)信息，把鐘山縣山洪地質(zhì)災(zāi)害預(yù)警APP在手機(jī)端展示，為鐘山縣山洪地質(zhì)災(zāi)害預(yù)警提供技術(shù)的支撐。進(jìn)入危險(xiǎn)區(qū)時(shí)，手機(jī)實(shí)時(shí)定位功能數(shù)據(jù)傳入后臺(tái)，顯示鐘山縣其他位置的災(zāi)害危險(xiǎn)情況，啟動(dòng)語(yǔ)音報(bào)警機(jī)制。圖2為研究技術(shù)路線示意。

3.1 GIS提取特征因子

鐘山縣有豐富的松散體物質(zhì)，地形有溝狀縱橫的坡度和充沛的降雨量，對(duì)于松散土體和水的聚集有利。鐘山縣這種裸露花崗巖強(qiáng)風(fēng)化巖體、碎屑巖和變質(zhì)巖的地質(zhì)條件，地形較為高聳坡陡，縱橫溝壑坡度大，便于水源的大量聚集，成為山洪地質(zhì)災(zāi)害易發(fā)多發(fā)區(qū)。

在山洪地質(zhì)災(zāi)害影響因子中選取地貌、地形、降雨、土壤巖性、植被等進(jìn)行量化，在Arcgis中提取高程、坡度、降雨量、植被覆蓋指數(shù)、土壤松散度作為神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)的特征因子。圖3為鐘山縣地質(zhì)災(zāi)害點(diǎn)示意。

圖2 研究技術(shù)路線示意

圖3 鐘山縣地質(zhì)災(zāi)害點(diǎn)示意

3.2 降雨因子

山洪形成的滑坡泥石流災(zāi)害須同時(shí)具備3個(gè)基本條件：有豐富的松散體物質(zhì)，短時(shí)間內(nèi)有大量來(lái)水和有一定坡度(一般小于25°)的溝狀地形，對(duì)于松散土體和水的聚集有利。裸露花崗巖強(qiáng)風(fēng)化巖體、碎屑巖和變質(zhì)巖的地質(zhì)條件，地形較為高聳坡陡，縱橫溝壑坡度大便于水源的大量聚集，多為形成泥石流的地區(qū)。降雨量對(duì)山洪地質(zhì)災(zāi)害的形成影響因子比重較大，有學(xué)者提出降雨類型的滑坡預(yù)警模型[11-13]。根據(jù)氣象資料結(jié)合實(shí)地調(diào)查，5—9月是廣西鐘山縣突發(fā)地質(zhì)災(zāi)害的重點(diǎn)防范期，也是降雨型山洪地質(zhì)災(zāi)害多發(fā)期。

鐘山縣2018—2019年月均降雨量見(jiàn)圖4所示。通過(guò)2019年的降雨及山洪地質(zhì)災(zāi)害發(fā)生數(shù)據(jù)的分析表明：在5—9月降雨集中時(shí)段是75%以上山洪地質(zhì)災(zāi)害發(fā)生的時(shí)段，且持續(xù)降雨和突發(fā)強(qiáng)降雨對(duì)更容易引發(fā)山洪地質(zhì)災(zāi)害。調(diào)查發(fā)現(xiàn)：對(duì)于持續(xù)性降雨引發(fā)的山洪地質(zhì)災(zāi)害有滯后性，強(qiáng)降雨引發(fā)的山洪地質(zhì)災(zāi)害有即時(shí)性。

圖4 2018-2019年月均降雨量示意

4 結(jié)果分析

4.1 神經(jīng)網(wǎng)絡(luò)模擬

由于影響山洪地質(zhì)災(zāi)害各因子有高程、坡度、植被覆蓋率、土壤松散度和降雨等，這些因子數(shù)據(jù)的單位、量級(jí)的不同，不可以直接作為因子進(jìn)行MATLAB的深度學(xué)習(xí)神經(jīng)網(wǎng)絡(luò)工具箱分析，因此，在使用MATLAB的深度學(xué)習(xí)神經(jīng)網(wǎng)絡(luò)工具箱之前，對(duì)各個(gè)因子進(jìn)行異質(zhì)數(shù)據(jù)歸一化的處理顯得尤為重要。而使用神經(jīng)網(wǎng)絡(luò)工具箱解決山洪地質(zhì)災(zāi)害預(yù)測(cè)時(shí)最大的限制條件是如何使影響因子異質(zhì)數(shù)據(jù)歸一化。

把調(diào)查整理的435組數(shù)據(jù)輸入神經(jīng)網(wǎng)絡(luò)模型，取其中70%進(jìn)行訓(xùn)練，剩下的30%一半進(jìn)行驗(yàn)證一半數(shù)據(jù)用來(lái)預(yù)測(cè)，結(jié)果預(yù)測(cè)精度在83.8%以上。圖5為神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)示意，圖6為神經(jīng)網(wǎng)絡(luò)線性擬合分析示意。

綜合分析與MATLAB的深度學(xué)習(xí)神經(jīng)網(wǎng)絡(luò)工具箱分析結(jié)果可以得出參數(shù)的歸一化處理與條件判斷處理均能達(dá)到滿意的分析精度；分別使用歸一化函數(shù)和條件判斷兩種方法對(duì)MATLAB的深度學(xué)習(xí)神經(jīng)網(wǎng)絡(luò)工具箱數(shù)據(jù)進(jìn)行處理，神經(jīng)網(wǎng)絡(luò)均能達(dá)到令人滿意的準(zhǔn)確度，且兩種方式的準(zhǔn)確度相近。

圖5 神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)示意

圖6 神經(jīng)網(wǎng)絡(luò)線性擬合分析示意

4.2 開(kāi)發(fā)顯示

使用Matlab和Java結(jié)合編程，創(chuàng)建具有高度交互性和可視化效果的虛擬或遠(yuǎn)程應(yīng)用程序。在這種情況下，Matlab通過(guò)工具箱用于計(jì)算和處理數(shù)據(jù)，而Java提供后臺(tái)支持和Java應(yīng)用程序的用戶交互[14-18]。選擇leaflet開(kāi)源地圖的天地圖插件作為顯示地圖，相對(duì)于比一般地圖API，天地圖在解決經(jīng)緯度換算偏移誤差的問(wèn)題上比較有優(yōu)勢(shì)。圖7為手機(jī)端APP顯示示意。

APP完善后將為對(duì)鐘山縣公眾開(kāi)放，界面簡(jiǎn)潔易于普通民眾操作。手機(jī)進(jìn)入危險(xiǎn)區(qū)時(shí)，手機(jī)實(shí)時(shí)定位功能數(shù)據(jù)傳入后臺(tái)，啟動(dòng)語(yǔ)音報(bào)警機(jī)制，也可以通過(guò)點(diǎn)擊地圖了解鐘山縣其他位置的災(zāi)害危險(xiǎn)情況。民眾也可通過(guò)點(diǎn)擊手機(jī)地圖上相應(yīng)位置得到對(duì)應(yīng)的危險(xiǎn)等級(jí)報(bào)送。

圖7 手機(jī)端APP顯示示意

5 結(jié)語(yǔ)

本文通過(guò)對(duì)鐘山縣已發(fā)生山洪地質(zhì)災(zāi)害的隱患點(diǎn)調(diào)查研究，基于山洪地質(zhì)災(zāi)害非線性多場(chǎng)耦合預(yù)警模型，構(gòu)建了適于鐘山縣的山洪地質(zhì)災(zāi)害預(yù)警APP應(yīng)用，實(shí)現(xiàn)實(shí)時(shí)預(yù)警鐘山縣災(zāi)害危險(xiǎn)等級(jí)信息，初步實(shí)現(xiàn)了從使用MATLAB與Java的混合開(kāi)發(fā)方式，實(shí)現(xiàn)了數(shù)據(jù)的在線運(yùn)算，彌補(bǔ)了Java模型運(yùn)算不足的缺點(diǎn)，也解決了MATLAB開(kāi)發(fā)后臺(tái)的缺陷，Leaflet1.6在線開(kāi)源地圖使山洪地質(zhì)災(zāi)害預(yù)警模型的理念轉(zhuǎn)變成為手機(jī)APP實(shí)際應(yīng)用。預(yù)警APP相比其他模型具有提高了實(shí)用性、運(yùn)算速度在毫秒級(jí)別、運(yùn)算能力便捷及其通用性良好，并且與山洪地質(zhì)災(zāi)害預(yù)警模型相結(jié)合，模型的研究得到了實(shí)際運(yùn)用，為鐘山縣的山洪地質(zhì)災(zāi)害預(yù)警提供了技術(shù)支撐。