基于GIS 的青海玉樹地區(qū)礦山地質(zhì)災(zāi)害風(fēng)險區(qū)劃研究

解立棟，劉 歡

（1.青海省水文地質(zhì)及地?zé)岬刭|(zhì)重點實驗室，青海 西寧 810008；2.青海省水文地質(zhì)工程地質(zhì)環(huán)境地質(zhì)調(diào)查院，青海 西寧 810008）

礦產(chǎn)資源作為人類社會中不可或缺的能源，長期以來一直是人類賴以生存和社會發(fā)展的重要物質(zhì)基礎(chǔ)[1]。青海省作為礦產(chǎn)資源大省，其有著十分豐富的自然資源，迄今為止在青海發(fā)現(xiàn)的礦產(chǎn)地超過一千五百處，其中玉樹為青海省主要礦產(chǎn)地之一，礦種多，門類全，儲量豐富。然而，由于受到認識的局限性以及青海省的地質(zhì)條件等原因的影響，在對礦產(chǎn)資源進行大規(guī)模的開發(fā)利用的同時，也為礦山環(huán)境帶來了大量的問題，對該地區(qū)礦業(yè)的可持續(xù)發(fā)展造成了一定程度上的制約。受本區(qū)自然地質(zhì)條件以及礦山開采產(chǎn)生的地質(zhì)環(huán)境問題帶來的影響，礦區(qū)地質(zhì)災(zāi)害的預(yù)測及防治工作是礦山開發(fā)工作中的一大難題。本文應(yīng)用GIS 技術(shù)對青海玉樹地區(qū)礦山地質(zhì)災(zāi)害進行風(fēng)險區(qū)劃研究，通過對青海玉樹地區(qū)能源礦產(chǎn)資源集中開采區(qū)的礦山地質(zhì)環(huán)境的調(diào)查，探清區(qū)內(nèi)礦山地質(zhì)環(huán)境問題及其成因，進而進一步推測其可能造成的危害，分析區(qū)域地質(zhì)環(huán)境問題可能帶來的地質(zhì)災(zāi)害，為礦山地質(zhì)環(huán)境的管理提供理論依據(jù)，為礦山地質(zhì)災(zāi)害的避免提供參考，研究結(jié)果對礦產(chǎn)資源的合理開發(fā)利用以及礦山地質(zhì)災(zāi)害的防治有著一定的理論意義和實用價值。

1 青海玉樹礦山地質(zhì)環(huán)境背景

青海玉樹地區(qū)位于我國西南“三江”構(gòu)造帶北段，處于青藏高原碰撞造山帶東北部，是滇藏地區(qū)的一部分，屬于東特提斯地質(zhì)構(gòu)造系統(tǒng)[2]。研究區(qū)處于唐古拉—昌都陸塊內(nèi)部幾個三級構(gòu)造單元的交界部位，不同構(gòu)造單元在此交切，各級分界斷裂在此交會，造成區(qū)內(nèi)斷裂構(gòu)造發(fā)育。青海玉樹地區(qū)位于青藏高原的高原面上，整體以高山、峽谷和山原為主，其地勢西北和中部高，東南和東部低，平均海拔較高，境內(nèi)高海拔山峰眾多，大部分終年積雪。區(qū)內(nèi)地層由于經(jīng)受了復(fù)雜的造山運動的改造，形成一系列緊密線型褶皺、表層寬展型褶皺。礦區(qū)的斷裂構(gòu)造主要為北西-南東向及北北東向兩組。其中前者與區(qū)域主構(gòu)造方向一致，為一條順層韌性剪切帶，帶寬300m～500m，傾向北東，控制火山巖帶的空間展布。后者則為后期構(gòu)造，切割火山巖帶，并對礦體起到一定的破壞作用，但由于斷層規(guī)模不大，破壞影響較小。

2 研究區(qū)礦山地質(zhì)問題及其影響因素

2.1 礦山地質(zhì)主要存在的問題

由于礦山開發(fā)產(chǎn)生的地質(zhì)環(huán)境問題主要表現(xiàn)在四個方面。一是礦山開發(fā)引發(fā)地質(zhì)災(zāi)害，使山體產(chǎn)生開裂、滑坡、泥石流等問題；二是礦山開發(fā)破壞地下含水層，使土壤和水環(huán)境遭到污染；三是地表植被破壞嚴重，水土流失加劇，造成山體破損、基巖裸露；四是煤礦的開采引起地面塌陷，破壞土地資源。

2.2 影響礦山地質(zhì)的主要因素分析

導(dǎo)致礦山地質(zhì)災(zāi)害的自然因素主要包括地層巖性、地質(zhì)構(gòu)造特點、礦產(chǎn)資源的分布情況、工程地質(zhì)特性等。研究區(qū)內(nèi)由于受到長期開采的影響，自然生態(tài)環(huán)境條件較脆弱。礦山地質(zhì)主要存在的災(zāi)害包括礦區(qū)積存廢礦渣數(shù)量巨大，礦渣型泥石流災(zāi)害隱患嚴重；采空區(qū)引發(fā)的山體崩塌、滑坡等地質(zhì)災(zāi)害隱患嚴重；含水層疏干破壞嚴重，導(dǎo)致地下水外泄，水位下降嚴重，水土流失嚴重。

煤礦地下開采會引起的地面塌陷等礦山地質(zhì)災(zāi)害，其中以采空塌陷尤為嚴重。從調(diào)查統(tǒng)計分析得知，影響青海玉樹地區(qū)礦山地質(zhì)災(zāi)害的主要問題是采空區(qū)塌陷。

3 基于GIS的研究區(qū)地質(zhì)災(zāi)害風(fēng)險區(qū)劃設(shè)計

GIS 技術(shù)能夠很好地展現(xiàn)地理空間分析算法[3]，為災(zāi)害風(fēng)險區(qū)劃分析提供較好的技術(shù)支持。本研究以礦區(qū)地質(zhì)災(zāi)害風(fēng)險區(qū)劃分析為核心，引入開源GIS 軟件作為技術(shù)支撐。

3.1 設(shè)計礦山地質(zhì)災(zāi)害風(fēng)險區(qū)域劃分模型

影響礦山地質(zhì)災(zāi)害風(fēng)險區(qū)劃結(jié)果的因素眾多，在對礦山地質(zhì)環(huán)境進行風(fēng)險區(qū)劃研究時，建立風(fēng)險區(qū)域劃分模型，選用正確的劃分方法，遵循實用性、相對獨立性等原則，保證結(jié)果的科學(xué)性與客觀性。

由前文分析得知，青海玉樹地區(qū)礦山地質(zhì)災(zāi)害問題的特征和強度與礦山地質(zhì)環(huán)境背景和礦山開發(fā)強度是密切相關(guān)的。因此根據(jù)上述指標(biāo)選取的原則，結(jié)合青海玉樹地區(qū)的區(qū)域地質(zhì)環(huán)境背景特點及調(diào)查統(tǒng)計情況，結(jié)合造成青海玉樹地區(qū)礦山地質(zhì)災(zāi)害的原因，建立青海玉樹地區(qū)礦山地質(zhì)災(zāi)害風(fēng)險評價模型。

采用綜合指數(shù)法評價青海玉樹地區(qū)地質(zhì)災(zāi)害易發(fā)分區(qū)。首先為制作地質(zhì)災(zāi)害因素圖，再加權(quán)求和，按一定的分割方法對綜合指數(shù)進行分割，以達到分區(qū)目的。其計算公式如下：

式中：S—綜合指數(shù)；Fij—影響因素及分類，其中i 表示影響因子的數(shù)量，j 表示影響因子的分類數(shù)量。通過該模型可以完成對青海玉樹地區(qū)的風(fēng)險區(qū)域劃分。

3.2 礦山地質(zhì)災(zāi)害風(fēng)險區(qū)域劃分

青海玉樹地區(qū)礦山地質(zhì)災(zāi)害主要包括崩塌-滑坡、泥石流、巖溶塌陷、地面沉降、采空塌陷和地震液化等，采礦活動、地形地貌、道路建設(shè)、降水量、地質(zhì)環(huán)境條件等因子是主要影響因素。

根據(jù)上述青海玉樹地區(qū)礦山地質(zhì)災(zāi)害的主要影響因素，對該地區(qū)風(fēng)險區(qū)域進行劃分計算。除地面沉降外，其它地質(zhì)災(zāi)害土地資源風(fēng)險區(qū)劃圖是以地質(zhì)災(zāi)害危險性區(qū)劃圖、地質(zhì)災(zāi)害土地資源綜合易損性區(qū)劃圖與土地資源價值密度區(qū)劃圖為基礎(chǔ)，采用風(fēng)險評價原理，結(jié)合青海玉樹地區(qū)的實際地質(zhì)狀況，疊加生成青海玉樹地區(qū)各個單災(zāi)種地質(zhì)災(zāi)害土地資源風(fēng)險區(qū)劃圖。其疊加公式如下：

式中：R—土地資源風(fēng)險，萬元/a；P—地質(zhì)災(zāi)害發(fā)生概率；V—土地資源綜合易損性；W—土地資源總價值。

由于采空塌陷為青海玉樹地區(qū)礦山地質(zhì)的主要災(zāi)害，為此計算該地區(qū)地面沉降風(fēng)險。地面沉降對土地資源造成的直接損失主要體現(xiàn)在地面標(biāo)高損失，主要是對城區(qū)、重要交通干線、重要工程的標(biāo)高損失，風(fēng)險按修復(fù)成本法計算，即它等于人工填土夯筑使地面恢復(fù)到地面沉降發(fā)生前的標(biāo)高所需的費用，計算公式如下：

式中：S—某鄉(xiāng)鎮(zhèn)城區(qū)、重要交通干線、重要工程的占地面積，m2；Vm—地面沉降平均速率，m/a；Jd：人工填土夯筑單價，元/m3。

最后對地質(zhì)災(zāi)害綜合風(fēng)險進行計算，即崩塌-滑坡、泥石流、巖溶塌陷、地面沉降、采空塌陷和地震液化等風(fēng)險疊加，疊加公式如下：

式中：Rq：地質(zhì)災(zāi)害土地資源綜合風(fēng)險；Rqi：某類地質(zhì)災(zāi)害土地資源風(fēng)險。

按公式（2）、（3）、（4）對青海玉樹地區(qū)地質(zhì)災(zāi)害可能造成的風(fēng)險以鄉(xiāng)鎮(zhèn)為單元進行預(yù)測評估，得到該地區(qū)礦山地質(zhì)災(zāi)害風(fēng)險區(qū)域劃分的結(jié)果。

3.3 與傳統(tǒng)風(fēng)險區(qū)域劃分方式的對比

由于青海玉樹地區(qū)處于地震帶上，且其礦山地質(zhì)的主要災(zāi)害為采空塌陷，因此以采空塌陷以及地震液化作為對比目標(biāo)，以采空塌陷、地震液化造成的經(jīng)濟損失，即物質(zhì)財富風(fēng)險作為區(qū)劃標(biāo)準(zhǔn)，將傳統(tǒng)風(fēng)險區(qū)域劃分的方法與本文風(fēng)險區(qū)域劃分的方法進行對比。采空塌陷及地震液化物質(zhì)財富風(fēng)險區(qū)劃標(biāo)準(zhǔn)如表1 所示。

表1 采空塌陷及地震液化物質(zhì)財富風(fēng)險區(qū)劃標(biāo)準(zhǔn)

根據(jù)物質(zhì)財富風(fēng)險區(qū)劃標(biāo)準(zhǔn)能夠確定青海玉樹地區(qū)采空塌陷及地震液化的物質(zhì)財富風(fēng)險曲線的分布規(guī)律。本文使用的風(fēng)險區(qū)域劃分模型以及風(fēng)險區(qū)域劃分方法不僅能夠完成對于礦山地質(zhì)災(zāi)害風(fēng)險區(qū)劃研究，滿足地質(zhì)風(fēng)險區(qū)劃的要求，而且顯著優(yōu)于傳統(tǒng)對于風(fēng)險區(qū)域的劃分方式，得到的結(jié)果更準(zhǔn)確，其數(shù)據(jù)更有說服力。

4 結(jié)束語

本文以GIS 技術(shù)為基礎(chǔ)，對青海玉樹地區(qū)礦山地質(zhì)災(zāi)害風(fēng)險區(qū)劃進行研究，本研究可以從區(qū)域上對青海玉樹地區(qū)的地質(zhì)災(zāi)害進行預(yù)測、預(yù)防，加強減災(zāi)防災(zāi)工作的開展與推進，能夠在一定程度上保護當(dāng)?shù)厝嗣竦纳敭a(chǎn)安全，促進青海玉樹地區(qū)礦山的開發(fā)與開采。取得的成果可為本地區(qū)礦山的開采工作以及其它地區(qū)地質(zhì)災(zāi)害風(fēng)險研究工作提供一定的參考。