肖陽(yáng)，葉唐進(jìn)

拉薩市郊地質(zhì)災(zāi)害預(yù)測(cè)及防治輔助決策

肖陽(yáng)，葉唐進(jìn)

(西藏大學(xué)工學(xué)院，拉薩 850000）

拉薩市郊區(qū)的自然環(huán)境和地質(zhì)環(huán)境較為脆弱，在雨量充沛的6~9月易形成崩塌、滑坡和泥石流等地質(zhì)災(zāi)害。為了對(duì)地質(zhì)災(zāi)害進(jìn)行有效的預(yù)測(cè)和防治輔助決策，根據(jù)野外調(diào)查資料，結(jié)合模糊數(shù)學(xué)、小波人工神經(jīng)網(wǎng)絡(luò)和灰色決策等理論進(jìn)行降雨量預(yù)測(cè)、危險(xiǎn)性計(jì)算以及防治輔助決策，并運(yùn)用Java SE Development Kit 8平臺(tái)開(kāi)發(fā)了地質(zhì)災(zāi)害預(yù)測(cè)及防治輔助決策系統(tǒng)。

地質(zhì)災(zāi)害；預(yù)測(cè)；輔助決策；拉薩

拉薩市位于喜馬拉雅山脈北側(cè)，屬高原溫帶半干旱季風(fēng)氣候，日照充足，空氣干燥，晝夜溫差大，降水量小，蒸發(fā)量大。由于拉薩河谷兩側(cè)海拔較高，植被覆蓋率低，山體寒凍風(fēng)化較為嚴(yán)重，其自然環(huán)境和地質(zhì)環(huán)境較為脆弱，遇大暴雨容易形成泥石流、崩塌、滑坡等地質(zhì)災(zāi)害[1]。

拉薩市郊區(qū)地質(zhì)災(zāi)害的研究，陳之強(qiáng)[2]等對(duì)色拉寺后山危巖的穩(wěn)定性進(jìn)行了分析并提出防治措施，王連宏[3]等對(duì)哲蚌寺泥石流的危害性進(jìn)行了計(jì)算分析并提出了相應(yīng)的治理措施，朱坤[4]等對(duì)墨竹工卡縣中學(xué)后山的滑坡穩(wěn)定性進(jìn)行了分析并進(jìn)行了抗滑樁的設(shè)計(jì)，張永昇[5]對(duì)甲瑪?shù)V區(qū)的融雪型碎屑流進(jìn)行了穩(wěn)定性分析并進(jìn)行數(shù)值模擬加以驗(yàn)證。但是對(duì)于拉薩市郊區(qū)地質(zhì)災(zāi)害的預(yù)測(cè)和防治輔助決策研究目前還沒(méi)有學(xué)者涉及。因此通過(guò)野外實(shí)地勘察，并結(jié)合模糊數(shù)學(xué)、小波人工神經(jīng)網(wǎng)絡(luò)和灰色決策等理論，運(yùn)用計(jì)算機(jī)編程和數(shù)據(jù)庫(kù)等技術(shù)，開(kāi)發(fā)了地質(zhì)災(zāi)害預(yù)測(cè)及防治輔助決策系統(tǒng)，進(jìn)而對(duì)地質(zhì)災(zāi)害發(fā)生的時(shí)間和空間進(jìn)行預(yù)測(cè)，同時(shí)提出防治輔助決策。

1 地質(zhì)背景

1.1 地形地貌與地層巖性

研究區(qū)屬岡底斯山脈東延部分的藏南高山寬谷地貌單元，總體為高山挾持河谷的地貌形態(tài)。受地質(zhì)構(gòu)造的控制以及拉薩河、堆龍河等水系侵蝕切割的影響，區(qū)內(nèi)形成了梁谷相間的高山河谷地形，地勢(shì)總體呈現(xiàn)南北高中部低、嶺高谷深的特點(diǎn)[6]。該區(qū)屬于岡底斯—騰沖地層區(qū)、拉薩沃卡分區(qū)，出露的地層主要有上侏羅系、白堊系、古近系及第四系。除第四系，其余地層分布面積約401.40km2，多呈不整合接觸，與巖漿巖呈侵入接觸。第四系分布面積363.81km2，成因類型有沖積、洪積、沖洪積、坡洪積、坡積、沼澤堆積、風(fēng)積等[7]。

1.2 地質(zhì)構(gòu)造與水文地質(zhì)

該區(qū)位于雅魯藏布江以北東西向構(gòu)造中部地帶，即拉薩—林芝巖帶，南側(cè)以雅魯藏布江深大斷裂為界與喜馬拉雅弧形構(gòu)造相鄰。處于當(dāng)雄羊八井-尼木與桑日-錯(cuò)那活動(dòng)構(gòu)造帶間的“安全島”，距離雅魯藏布江深大斷裂和當(dāng)雄羊八井-尼木地震活躍帶60~100km。該區(qū)域河流、溪溝發(fā)育，呈樹(shù)枝狀展布，屬雅魯藏布江水系。常年性河流主要有拉薩河、堆龍河和流沙河，其余皆為季節(jié)性溪流。地表水資源豐富，但時(shí)空分布極為不均。地下水按賦存空間與水力性質(zhì)可劃分為松散巖類孔隙水、基巖裂隙水、碳酸鹽巖類裂隙溶洞水三大類型[7]。

1.3 氣象與不良地質(zhì)

該區(qū)屬內(nèi)陸高原溫帶半干旱氣候[7]，區(qū)內(nèi)多年平均氣溫8℃，最高氣溫29.4℃，最低氣溫-16.5℃，多年平均降水量424.56mm，最大年降水量613.8mm，最小年降水量229.6mm，5~9月集中了全年94.8%的降水量，多為歷時(shí)較短的中到大雨，最大日降水量為61.7mm。區(qū)內(nèi)蒸發(fā)量為2 306.70mm，寒凍風(fēng)化較為強(qiáng)烈。

研究區(qū)的不良地質(zhì)主要為崩塌、滑坡、泥石流、碎屑流和地面塌陷等幾種類型，不良地質(zhì)的發(fā)生與地形地貌、地層巖性、大氣降水等因素密切相關(guān)。其中，泥石流主要分布在拉薩河、堆龍河兩岸的支溝與山前坡麓地帶。崩塌主要分布在風(fēng)化較為嚴(yán)重的花崗巖邊坡地段，部分發(fā)育于溝谷的自然陡坡區(qū)?；铝阈欠植加诶_河兩岸山前斜坡地段。碎屑流主要分布于海拔較高、寒凍風(fēng)化較嚴(yán)重的山頂[7]。

2 預(yù)測(cè)及防治輔助決策的計(jì)算方法

對(duì)于拉薩市地質(zhì)災(zāi)害的預(yù)測(cè)及防治輔助決策，本文通過(guò)野外調(diào)查分析，再結(jié)合小波人工神經(jīng)網(wǎng)絡(luò)計(jì)算出2016年拉薩市降雨的規(guī)律，然后利用模糊數(shù)學(xué)計(jì)算，得到其時(shí)間序列預(yù)測(cè)結(jié)果，最后利用灰色決策對(duì)地質(zhì)災(zāi)害提出防治輔助決策。

2.1 小波人工神經(jīng)網(wǎng)絡(luò)

運(yùn)用小波人工神經(jīng)網(wǎng)絡(luò)組合模型[8]，預(yù)測(cè)2016年拉薩市的降雨量，其原始數(shù)據(jù)來(lái)源于歐洲中期天氣預(yù)報(bào)中心（ECMWF）拉薩市1979~2015年的降雨數(shù)據(jù)。然后利用Matlab進(jìn)行小波快速分解，采用小波為bior 2.4，分解級(jí)數(shù)為二級(jí)，得到小波分解序列D2，D1，A2，見(jiàn)圖1。

圖1 小波分解

圖2 時(shí)間序列響應(yīng)

經(jīng)分析，WANN結(jié)構(gòu)為3-15-1，即輸入為[D2,D1,A2]，由37年降雨數(shù)據(jù)（1979~2015年）建模，訓(xùn)練1 000次，得出2016年拉薩市的年降雨量為474.15mm。其時(shí)間序列響應(yīng)見(jiàn)圖2，誤差在允許范圍內(nèi)，可用來(lái)進(jìn)行預(yù)測(cè)。

2016年實(shí)際降雨量與預(yù)測(cè)降雨量對(duì)比統(tǒng)計(jì)見(jiàn)表1，因今年的數(shù)據(jù)只有1~5月，因此只對(duì)比了前五個(gè)月的降雨量，從表中可以看出，預(yù)測(cè)降雨量與實(shí)際降雨量之間的誤差都在2mm之內(nèi)，誤差在允許范圍，可用來(lái)進(jìn)行計(jì)算。

2.2 模糊數(shù)學(xué)

表1 拉薩市2016年1~5月實(shí)際降雨和預(yù)測(cè)降雨對(duì)比結(jié)果（單位：mm）

采用模糊數(shù)學(xué)中的模糊綜合評(píng)判[9]，以野外調(diào)查的拉薩市典型的六處地質(zhì)災(zāi)害點(diǎn)為例建立因素集U1，因素集U1={地層巖性，坡度，坡高，結(jié)構(gòu)面發(fā)育，地質(zhì)構(gòu)造，新構(gòu)造運(yùn)動(dòng)與地震，人類工程活動(dòng)，月降雨量}。在模糊綜合評(píng)判過(guò)程中，因崩塌與滑坡共同點(diǎn)較多，都選用以上的8個(gè)評(píng)價(jià)因子，泥石流選用另外的8個(gè)因子，分別為：巖性因素、流域面積、主溝平均比降、溝床彎曲系數(shù)、流域切割密度、單位面積物源量、不良地質(zhì)現(xiàn)象、月降雨量。

將地質(zhì)災(zāi)害穩(wěn)定狀態(tài)分為4個(gè)等級(jí)，評(píng)價(jià)集V為：V={穩(wěn)定、基本穩(wěn)定、欠穩(wěn)定、不穩(wěn)定}。采用半梯形模糊分布與專家打分法確定隸屬函數(shù)，計(jì)算出單因素評(píng)判矩陣，采用層次分析法[10]確定各因子權(quán)重，再使用上文的月降雨量進(jìn)行預(yù)測(cè)，選擇崩塌1、崩塌2、崩塌3、滑坡1、滑坡2、泥石流分別編號(hào)為1#、2#、3#、4#、5#、6#。然后進(jìn)行模糊綜合評(píng)判，得到評(píng)價(jià)集中“不穩(wěn)定”狀態(tài)的隸屬度值的時(shí)間序列預(yù)測(cè)結(jié)果（圖3）。

圖3 典型地質(zhì)災(zāi)害點(diǎn)時(shí)間序列預(yù)測(cè)結(jié)果折線圖

從圖3可知，在12個(gè)月中，“不穩(wěn)定”隸屬度值的峰值一般出現(xiàn)在7月和8月，兩側(cè)分別呈遞減趨勢(shì)，因此危險(xiǎn)月份為6~9月，7~8月最危險(xiǎn)，其余月份則比較安全。

2.3 灰色決策

采用灰色決策中的多目標(biāo)加權(quán)灰靶決策[11]，以拉薩市典型的6處地質(zhì)災(zāi)害點(diǎn)的時(shí)間序列預(yù)測(cè)結(jié)果進(jìn)行決策。首先建立事件集、對(duì)策集及決策方案集，將最危險(xiǎn)的災(zāi)害點(diǎn)作為事件a1，事件集A={a1}。選擇崩塌1、崩塌2、崩塌3、滑坡1、滑坡2、泥石流建立對(duì)策集B={b1,b2,b3,b4,b5,b6}，由事件A和對(duì)策B構(gòu)造決策方案集：

表2 各個(gè)目標(biāo)及相應(yīng)指標(biāo)的決策權(quán)重

然后確定三個(gè)決策目標(biāo)：最大值、變化率1、變化率2，其中最大值與變化率1為效益型指標(biāo)，變化率2屬于成本型指標(biāo)。最后確定個(gè)目標(biāo)的決策權(quán)（表2）。

求得各目標(biāo)的效果樣本向量為：

U(1)=（0.45,0.30,0.40,0.29,0.29,0.32）

U(2)=（0.13,0.06,0.17,0.11,0.05,0.10）

U(3)=（0.13,0.05,0.16,0.11,0.06,0.09）

根據(jù)計(jì)算得出的多組數(shù)據(jù)，設(shè)定最大值的臨界值為0.2，變化率1的臨界值為0.05，設(shè)定變化率2的臨界值為0.1，求得一致效果測(cè)度向量為：

R(1)=（1.0,1.0,0.80,0.36,0.36,0.48）

R(2)=（1.0,0.89,0.22,0.67,0.00,0.56）

R(3)=（-1.00,-0.75,1.0,-0.25,1.0,0.25）

3 軟件系統(tǒng)實(shí)現(xiàn)

拉薩市地質(zhì)災(zāi)害預(yù)測(cè)及防治輔助決策系統(tǒng)必須具備專業(yè)性、實(shí)用性與通用性，使其具有數(shù)據(jù)儲(chǔ)存管理、實(shí)時(shí)更新預(yù)測(cè)和為政府相關(guān)部門提出防治輔助決策的功能。

圖4 系統(tǒng)登錄界面

圖5 系統(tǒng)主頁(yè)界面

3.1 系統(tǒng)組成與系統(tǒng)功能

拉薩市地質(zhì)災(zāi)害預(yù)測(cè)及防治輔助決策系統(tǒng)主要包括三大模塊：危險(xiǎn)性計(jì)算模塊、預(yù)測(cè)及防治輔助決策模塊、數(shù)據(jù)庫(kù)管理模塊。系統(tǒng)基于Java SE Development Kit 8平臺(tái)開(kāi)發(fā)，可運(yùn)行在具備Java SE 8環(huán)境的Windows、Linux、Mac等個(gè)人計(jì)算機(jī)上。系統(tǒng)采用NetBeans 8.1 IDE開(kāi)發(fā)工具，采用Apache Derby開(kāi)源數(shù)據(jù)庫(kù)。登錄系統(tǒng)作為系統(tǒng)的入口，輸入正確的用戶名與密碼方可進(jìn)入主頁(yè)，有控制用戶權(quán)限、記錄用戶行為、保護(hù)操作安全的作用，圖4為登錄界面。其數(shù)據(jù)庫(kù)管理是對(duì)已有數(shù)據(jù)進(jìn)行管理的功能，方便錄入、查看、修改災(zāi)害點(diǎn)數(shù)據(jù)。在主頁(yè)部分，除菜單欄外展示了拉薩市郊區(qū)地質(zhì)災(zāi)害的簡(jiǎn)介、照片、小波分析結(jié)果、時(shí)間序列預(yù)測(cè)結(jié)果，還具有快速引導(dǎo)操作的作用，見(jiàn)圖5。

3.2 危險(xiǎn)性計(jì)算

采用小波人工神經(jīng)網(wǎng)絡(luò)組合模型，使用Matlab等軟件進(jìn)行計(jì)算，預(yù)測(cè)出拉薩市2016年的降雨量，再進(jìn)行權(quán)重計(jì)算，然后進(jìn)行地質(zhì)災(zāi)害的隸屬度計(jì)算及綜合評(píng)判，最后得出地質(zhì)災(zāi)害的危險(xiǎn)性評(píng)價(jià)結(jié)果。崩塌滑坡權(quán)重計(jì)算、隸屬度計(jì)算以及綜合評(píng)判界面詳見(jiàn)圖6、圖7。

圖6 崩塌滑坡權(quán)重計(jì)算

圖7 崩塌滑坡模糊綜合評(píng)判

3.3 時(shí)間序列預(yù)測(cè)及防治輔助決策

時(shí)間序列預(yù)測(cè)方法是將2016年的月降雨序列帶入模糊綜合評(píng)判中進(jìn)行計(jì)算，得出12個(gè)評(píng)價(jià)集中“不穩(wěn)定”狀態(tài)的隸屬度值，然后繪圖，并為輔助決策作準(zhǔn)備（圖8）。防治輔助決策采用灰色決策模型中的多目標(biāo)加權(quán)灰靶決策模型，數(shù)據(jù)來(lái)源于上文中的結(jié)果，將其進(jìn)行計(jì)算，并得出最終輔助決策結(jié)果（圖9）。

圖8 時(shí)間序列預(yù)測(cè)

圖9 輔助決策界面

3.4 數(shù)據(jù)庫(kù)管理及赤平投影

為了更好地實(shí)現(xiàn)系統(tǒng)功能，本系統(tǒng)中加入了數(shù)據(jù)庫(kù)管理功能，采用Apache Derby數(shù)據(jù)庫(kù)，并將其內(nèi)嵌至軟件中，可較全面地錄入地質(zhì)災(zāi)害點(diǎn)的名稱、地點(diǎn)、編號(hào)、地理位置、邊坡環(huán)境、邊坡基本特征、失穩(wěn)情況等數(shù)據(jù)，便于查詢與使用。數(shù)據(jù)管理界面見(jiàn)圖10。系統(tǒng)中極射赤平投影主要用來(lái)表示結(jié)構(gòu)面的穩(wěn)定情況，可以對(duì)單個(gè)巖體的穩(wěn)定性進(jìn)行判別，詳見(jiàn)圖11。

圖10 數(shù)據(jù)管理界面

圖11 赤平投影界面

4 結(jié)論

通過(guò)對(duì)拉薩市郊區(qū)幾處典型地質(zhì)災(zāi)害點(diǎn)的野外調(diào)查，并結(jié)合模糊數(shù)學(xué)、小波人工神經(jīng)網(wǎng)絡(luò)、灰色決策理論和計(jì)算機(jī)技術(shù)分析研究，得出以下幾點(diǎn)結(jié)論。

1）通過(guò)野外調(diào)查和理論分析得出拉薩市地質(zhì)災(zāi)害主要影響因素為降雨、地震等，地質(zhì)災(zāi)害發(fā)生主要集中在雨季的6~9月。

2）通過(guò)防治輔助決策分析，得出先治理比較危險(xiǎn)崩塌3，其次是崩塌2、崩塌1和泥石流，最后治理滑坡2和滑坡1。

3）運(yùn)用Java SE Development Kit 8平臺(tái)開(kāi)發(fā)了地質(zhì)災(zāi)害預(yù)測(cè)及防治輔助決策系統(tǒng)軟件。

[1] 王曉軍，程紹敏．西藏主要?dú)夂蛱卣鞣治鯷J]．高原山地氣象研究，2009，04：81-84．

[2] 陳之強(qiáng)，王宏，姜福堂，等．拉薩色拉寺后山危巖防治措施研究[J]．市場(chǎng)與技術(shù)，2013，20(6)：5-7．

[3] 王連宏，葉唐進(jìn)，姜福堂，等．哲蚌寺泥石流的防治措施研究[J]．四川建材，2014，40(179)：147-149．

[4] 朱坤，葉唐進(jìn)，陳貴波，等．墨竹工卡縣中學(xué)后山滑坡體穩(wěn)定性分析研究[J]．四川建材，2014，40（182）：94-96．

[5] 張永昇，王萍淋，喬晨光，等．融雪型滑坡—碎屑流的穩(wěn)定性分析[J]．四川建材，2015，12（6）：234-236．

[6] 張豐述．拉薩市城市環(huán)境地質(zhì)評(píng)價(jià)[Z]．成都理工大學(xué)，2011．

[7] 地礦部九一五水文地質(zhì)工程地質(zhì)大隊(duì). 西藏拉薩河河谷平原綜合水文地質(zhì)勘查報(bào)告（1∶10萬(wàn)）[R]. 1991．

[8] 王文圣，丁晶，李躍清．水文小波分析[M]．北京：化學(xué)工業(yè)出版社，2005．

[9] 梁保送，曹殿立．模糊數(shù)學(xué)及其應(yīng)用[M]．北京：科學(xué)出版社，2007．

[10] 汪應(yīng)洛．系統(tǒng)工程簡(jiǎn)明教程[M]．北京：高等教育出版社，2009．

[11] 劉思峰，楊英杰，吳利豐．灰色系統(tǒng)理論及其應(yīng)用[M]．北京：科學(xué)出版社，2014．

Prediction and Control of Geohazards in the Suburban District of Lhasa

XIAO Yang YE Tang-jin

(School of Engineering, Tibet University, Lhasa 850000)

The natural and geological environment in the outskirts of Lhasa is much fragile which makes easy of geohazards such as collapse, landslide, debris flow, etc. from June to September. This study makes precipitation prediction, risk calculation and prevention and control of auxiliary decision-making based on the fuzzy mathematics and develops geological hazard prediction and prevention and control of auxiliary decision-making system by the Java SE Development Kit 8 platform.

geohazard prediction; prevention and control; auxiliary decision-making; fuzzy mathematics; Lhasa

P642.2

A

1006-0995（2017）03-0449-05

10.3969/j.issn.1006-0995.2017.03.022

2016-11-24

西藏自治區(qū)自然科學(xué)基金項(xiàng)目（2016ZR-15-7）和西藏自治區(qū)大學(xué)生創(chuàng)新性實(shí)驗(yàn)訓(xùn)練計(jì)劃項(xiàng)目（2015QCX026）資助

肖陽(yáng)（1994-），男，四川大竹人，本科學(xué)生，主要從事地質(zhì)災(zāi)害預(yù)測(cè)預(yù)報(bào)研究

葉唐進(jìn)（1981-），男，貴州銅仁人，副教授，主要從事地質(zhì)災(zāi)害相關(guān)方面的研究