張建萍,史慧君,曹金亮

(1.山西省地質(zhì)環(huán)境監(jiān)測中心,山西 太原 030024; 2.太原碧藍(lán)水利工程設(shè)計有限公司,山西 太原 030024)

礦產(chǎn)資源過度開發(fā)引起的一系列地質(zhì)環(huán)境問題不僅威脅到人民生命財產(chǎn)安全,還制約著經(jīng)濟(jì)的可持續(xù)發(fā)展。對礦山開發(fā)所引起地質(zhì)環(huán)境問題的風(fēng)險性進(jìn)行評價,可以為政府及有關(guān)部門合理規(guī)劃開發(fā)礦產(chǎn)資源、保護(hù)礦區(qū)生態(tài)環(huán)境及后期的恢復(fù)治理規(guī)劃提供基礎(chǔ)資料和科學(xué)依據(jù),具有重要意義。

大同市擁有豐富的煤炭及其它礦產(chǎn)資源,是全國礦產(chǎn)資源開發(fā)程度較高的地區(qū)之一。由于長期礦產(chǎn)開采,導(dǎo)致該地區(qū)地質(zhì)環(huán)境遭到嚴(yán)重破壞,地質(zhì)環(huán)境風(fēng)險持續(xù)加大,因此有必要開展相應(yīng)的地質(zhì)環(huán)境風(fēng)險評估工作,為下一步礦山地質(zhì)環(huán)境問題綜合防治提供依據(jù)。

礦山地質(zhì)環(huán)境風(fēng)險產(chǎn)生的根源一方面是由于自然地質(zhì)環(huán)境的脆弱性,另一方面是由人類礦業(yè)開發(fā)活動所引發(fā)。目前對礦山地質(zhì)環(huán)境的研究多側(cè)重于環(huán)境質(zhì)量綜合評價,定性或定量地評價礦業(yè)開發(fā)對周邊地質(zhì)環(huán)境造成的影響程度[1～3]。突出了人為活動對形成地質(zhì)環(huán)境問題的影響,而忽視了礦山地質(zhì)環(huán)境自身的脆弱性。本文在借鑒各類自然災(zāi)害風(fēng)險研究的基礎(chǔ)上,提出礦山地質(zhì)環(huán)境風(fēng)險評價思路,并以山西省大同市為例,綜合分析大同市礦山地質(zhì)環(huán)境遭受破壞的風(fēng)險性。

1 研究區(qū)概況

1.1 自然地質(zhì)條件

大同市位于山西省最北端,北臨內(nèi)蒙古,東與河北接壤,總面積為14 176 km2。該地區(qū)屬大陸性季風(fēng)氣候,風(fēng)多雨少,溫差較大；境內(nèi)河流眾多,分屬海河和黃河流域。本區(qū)屬構(gòu)造山地的山坡及山麓地帶,屬中低山地貌,地表大部分為第四系黃土覆蓋,溝谷發(fā)育。大同市基底為太古界集寧群變質(zhì)巖。上伏地層為寒武系、奧陶系、石炭系、二疊系、侏羅系、白堊系、新近系、第四系。本區(qū)構(gòu)造為一軸向大致NNE、南部轉(zhuǎn)為NE向的不對稱向斜。西北翼寬緩,地層近水平。東南翼狹窄,邊部地層直立或倒轉(zhuǎn)。較大斷層如青磁窯斷層、煤峪口斷層、白洞斷層、鵝毛口斷層等,見于區(qū)內(nèi)東南邊緣,走向NE—NNE[3]。

1.2 礦山地質(zhì)環(huán)境現(xiàn)狀

大同市礦產(chǎn)資源豐富,其中以煤炭資源儲量占我國1/8而著稱煤都,位居中國產(chǎn)煤城市之首。研究區(qū)現(xiàn)有各類礦山313 座,具體見表1。

表1 研究區(qū)礦山統(tǒng)計表Table 1 Statistics of mines in the study area

研究區(qū)發(fā)育的礦山地質(zhì)環(huán)境問題主要有采礦引發(fā)的地質(zhì)災(zāi)害(地面塌陷、地裂縫、崩塌、滑坡等),地形地貌的破壞,矸石、尾礦堆放,廢水對水環(huán)境的污染等。據(jù)2015年調(diào)查資料,研究區(qū)共發(fā)育各類地質(zhì)災(zāi)害454 處,其中地面塌陷、地裂縫433 處,崩塌17 處,滑坡4 處,造成直接經(jīng)濟(jì)損失達(dá)28 053.68 萬元。矸石堆放125 處,廢石堆放58 處,占用和破壞土地達(dá)623.9 km2(圖1)。

2 礦山地質(zhì)環(huán)境風(fēng)險評價方法

礦山地質(zhì)環(huán)境風(fēng)險理解為研究區(qū)的礦山地質(zhì)環(huán)境在當(dāng)?shù)氐V業(yè)開發(fā)活動影響下破壞的概率或可能性。礦業(yè)開發(fā)是風(fēng)險的施加者,研究區(qū)地質(zhì)環(huán)境是風(fēng)險的承載體。即礦業(yè)開發(fā)對礦山地質(zhì)環(huán)境的影響力度越大,造成的礦山地質(zhì)環(huán)境風(fēng)險越高；研究區(qū)地質(zhì)環(huán)境越脆弱,則抵抗破壞的能力越小,造成的礦山地質(zhì)環(huán)境風(fēng)險越高。礦業(yè)開發(fā)對礦山地質(zhì)環(huán)境的影響和地質(zhì)環(huán)境本身的脆弱性這兩種效應(yīng)相互作用就構(gòu)成了研究區(qū)的礦山地質(zhì)環(huán)境風(fēng)險性[3～7]。

圖1 大同市主要地質(zhì)環(huán)境問題圖Fig.1 Map of the major geo-environment problems in Datong

2.1 評價方法

目前,自然災(zāi)害風(fēng)險評價多采用加權(quán)比較法、綜合評判法(也稱綜合指數(shù)法)、模糊綜合評判法、概率法、神經(jīng)網(wǎng)絡(luò)以及GIS圖層疊加分析法等。本文采用GIS圖層疊加分析法與綜合評判法相結(jié)合對研究區(qū)礦山地質(zhì)環(huán)境風(fēng)險性進(jìn)行評價,用風(fēng)險性指數(shù)來反映研究區(qū)礦山地質(zhì)環(huán)境的風(fēng)險程度。

2.2 評價指標(biāo)體系

本文將礦山地質(zhì)環(huán)境評價系統(tǒng)分為礦業(yè)開發(fā)影響力和地質(zhì)環(huán)境脆弱性兩個子系統(tǒng),并以此確定評價指標(biāo)體系。

礦業(yè)開發(fā)影響力評價指標(biāo)體系包括兩個層面：要素層和指標(biāo)層。要素層包括采礦活動、突發(fā)型地質(zhì)環(huán)境問題和漸變型地質(zhì)環(huán)境問題,反映礦業(yè)開發(fā)造成的主要影響。指標(biāo)層即每一個要素又包括若干個單指標(biāo),來反映其內(nèi)在原因和潛在因素。

地質(zhì)環(huán)境脆弱性指標(biāo)體系要素層包括基礎(chǔ)地質(zhì)和自然地理兩個方面,來反映研究區(qū)礦山地質(zhì)背景。

綜上,構(gòu)建出礦山地質(zhì)環(huán)境風(fēng)險評價指標(biāo)體系(圖2)。

2.3 評價指標(biāo)分級量化

在參考中國地質(zhì)調(diào)查局制定的《礦山地質(zhì)環(huán)境調(diào)查評價規(guī)范》基礎(chǔ)上[8],將礦山地質(zhì)環(huán)境風(fēng)險程度分為高風(fēng)險、中等風(fēng)險、低風(fēng)險三個等級,礦業(yè)開發(fā)影響力和地質(zhì)環(huán)境脆弱性評價分為嚴(yán)重、較嚴(yán)重和一般三個等級。根據(jù)這三個等級的劃分,評價指標(biāo)也分為三個等級。并通過專家咨詢和實地調(diào)查制定了分級量化評分標(biāo)準(zhǔn),具體分級見表2。

圖2 大同市礦山地質(zhì)環(huán)境風(fēng)險評價指標(biāo)體系Fig.2 Index system of mine geo-environment risk assessment of Datong

表2 礦山地質(zhì)環(huán)境風(fēng)險評價指標(biāo)分級打分表Table 2 Index classification table of mine geo-environment risk assessment of Datong

3 大同市礦山地質(zhì)環(huán)境風(fēng)險評價實現(xiàn)

3.1 評價模型

參考災(zāi)害風(fēng)險計算公式：風(fēng)險度(R)=危險性(H)×易損性(V),礦山地質(zhì)環(huán)境風(fēng)險值可用礦業(yè)開發(fā)影響力和地質(zhì)環(huán)境脆弱性值來表達(dá)[6],即用以下數(shù)學(xué)模型來表示：

F=K×D

式中：F——礦山地質(zhì)環(huán)境風(fēng)險值；

K——礦業(yè)開發(fā)影響力值；

D——地質(zhì)環(huán)境脆弱性值。

礦業(yè)開發(fā)影響力和地質(zhì)環(huán)境脆弱性值可用模糊綜合評價模型來表示,以礦業(yè)開發(fā)影響力為例,綜合評價模型為：

式中：Kj——單要素的加權(quán)分值；

Wi——單指標(biāo)的權(quán)重值；

Fi——每一要素層中單指標(biāo)的分值；

n1——各要素中的單指標(biāo)個數(shù)；

K——礦業(yè)開發(fā)影響力值；

Wj——各要素的權(quán)重值；

Fj——各要素的加權(quán)分值；

n——要素個數(shù),此處n為3。

地質(zhì)環(huán)境脆弱性評價模型與礦業(yè)開發(fā)影響力評價模型相似,故不重復(fù)介紹。

3.2 評價指標(biāo)權(quán)重計算

指標(biāo)的權(quán)重反映不同評價因子間重要性程度的差異。對于進(jìn)行礦山地質(zhì)環(huán)境風(fēng)險評價時,選擇的評價因子既有定量指標(biāo),也有非定量指標(biāo),并且指標(biāo)數(shù)量較大,給確定權(quán)重帶來了一定的困難。應(yīng)選擇合適的定權(quán)方法來計算指標(biāo)權(quán)重。

本次礦山地質(zhì)環(huán)境風(fēng)險評價指標(biāo)權(quán)重采用層次分析法結(jié)合專家打分法來確定。

(1)組織專家組

專家組成員是在相關(guān)方面經(jīng)驗豐富的專家,對各因素的相對重要性進(jìn)行打分,構(gòu)成判斷矩陣。

(2)構(gòu)造判斷矩陣并計算權(quán)重

采用1～9標(biāo)度法構(gòu)造判斷矩陣,在統(tǒng)計分析軟件R中求出各判斷矩陣的最大特征向量。并對判斷矩陣進(jìn)行一致性檢驗,經(jīng)檢驗所構(gòu)建的判斷矩陣一致性很好。最后將特征向量歸一化處理,即為權(quán)重值W。

①判斷矩陣F-K和D層(表3)

表3 F-K和D層判斷矩陣及權(quán)重表Table 3 Judgment matrix and weight of layer F-K and D

②判斷矩陣K-K1、K2、K3,D-D1和D2層(表4～5)

表4 K-K1、K2、K3層判斷矩陣及權(quán)重Table 4 Judgment matrix and weight of layer K-K1、K2、K3

表5 D-D1和D2層判斷矩陣及權(quán)重Table 5 Judgment matrix and weight of layer D-D1 and D2

③判斷矩陣K1-K1i,K2-K2i,K3-K3i,D1-D1i,D2-D2i層(表6～10)

表6 K1-K1i層判斷矩陣及權(quán)重Table 6 Judgment matrix and weight of layer K1-K1i

表7 K2-K2i層判斷矩陣及權(quán)重Table 7 Judgment matrix and weight of layer K2-K2i

表8 K3-K3i層判斷矩陣及權(quán)重Table 8 Judgment matrix and weight of layer K3-K3i

表9 D1-D1i層判斷矩陣及權(quán)重Table 9 Judgment matrix and weight of layer D1-D1i

表10 D2-D2i層判斷矩陣及權(quán)重Table 10 Judgment matrix and weight of layer D2-D2i

3.3 評價結(jié)果

由上述評價模型和指標(biāo)權(quán)重,可形成如下計算公式：

K1=0.08×K11+0.11×K12+0.27×K13+0.54×K14；

K2=0.31×K21+0.31×K22+0.10×K23+0.10×K24+0.18×K25；

K3=0.23×K31+0.08×K32+0.69×K33；

K=0.22×K1+0.13×K2+0.65×K3；

D1=0.33×D11+0.67×D12；

D2=0.09×D21+0.07×D22+0.21×D23+0.63×D24；

D=0.4×D1+0.6×D2；

F=0.6×K+0.4×D。

根據(jù)已構(gòu)建的評價指標(biāo)體系,在GIS平臺繪制相應(yīng)的單指標(biāo)圖層。運用ArcGIS9.3軟件的圖層聯(lián)合功能(Union),根據(jù)上述公式,對各圖層屬性進(jìn)行加權(quán)疊加計算,得到一個綜合的圖層即為礦山地質(zhì)環(huán)境風(fēng)險評價分區(qū)圖(圖3)。

圖3 大同市礦山地質(zhì)環(huán)境風(fēng)險評價分區(qū)圖Fig.3 Zoning plan of mine geo-environment risk assessment of Datong

根據(jù)大同市礦山地質(zhì)環(huán)境風(fēng)險分區(qū)圖,結(jié)合大同市主要礦山地質(zhì)環(huán)境問題,可得出：

大同市礦山地質(zhì)環(huán)境高風(fēng)險區(qū)面積為1 277.26 km2,占總面積的9.01%。主要分布于大同市南郊區(qū)、新榮區(qū)、左云縣煤礦集中開采地及渾源南部、靈丘東北部等采石場分布地區(qū)。區(qū)內(nèi)礦山地質(zhì)環(huán)境問題發(fā)育較多,對礦山地質(zhì)環(huán)境破壞嚴(yán)重。主要礦山地質(zhì)環(huán)境問題為采煤引發(fā)的地面塌陷、地裂縫；采煤破壞地下含水層；矸石堆放和采石場對地形地貌的破壞等。

大同市礦山地質(zhì)環(huán)境中風(fēng)險區(qū)面積為1 505.49 km2,占總面積的10.62%。主要分布于大同市天鎮(zhèn)東部、陽高北部、渾源中部、靈丘北、南等金屬礦分布地區(qū)。區(qū)內(nèi)礦山地質(zhì)環(huán)境問題主要為采礦引發(fā)的地裂縫、尾礦堆放等,對區(qū)內(nèi)礦山地質(zhì)環(huán)境影響較嚴(yán)重。

大同市礦山地質(zhì)環(huán)境低風(fēng)險區(qū)面積為11 393.25 km2,占總面積的80.37%。主要分布于除高風(fēng)險區(qū)和中風(fēng)險區(qū)的未進(jìn)行礦業(yè)開發(fā)的區(qū)域,區(qū)內(nèi)礦產(chǎn)資源較少,對礦山地質(zhì)環(huán)境沒有破壞。

4 風(fēng)險防控建議

風(fēng)險防控主要從兩方面進(jìn)行,一是從風(fēng)險源頭對風(fēng)險的施加者進(jìn)行管控；二是對風(fēng)險的承災(zāi)體進(jìn)行保護(hù)提升其防災(zāi)能力。對于大同市地質(zhì)環(huán)境風(fēng)險性,根據(jù)不同的分區(qū)和主要礦山地質(zhì)環(huán)境問題,提出相應(yīng)的風(fēng)險防控建議[9～10]。

高風(fēng)險區(qū)：主要是煤礦集中開采地與露天采石場聚集地,形成的礦山地質(zhì)環(huán)境問題主要為地面塌陷、地裂縫,及對地形地貌的破壞。因此對本區(qū)的防控建議首先是加大治理力度,應(yīng)因地制宜,綜合治理。對地面塌陷、地裂縫,可填埋復(fù)墾；對地形地貌的破壞,可進(jìn)行復(fù)綠；同時加強礦業(yè)開發(fā)的準(zhǔn)入制度,嚴(yán)格礦山地質(zhì)環(huán)境保護(hù)制度,采用新技術(shù)、新的采礦方法,將礦山開采對環(huán)境的破壞控制在最小限度。要加強事中監(jiān)管,對礦山開采對環(huán)境的影響實施跟蹤監(jiān)測,發(fā)現(xiàn)問題及時督促礦山企業(yè)落實保護(hù)治理措施。

中風(fēng)險區(qū)：主要是金屬礦分布區(qū),形成的礦山地質(zhì)環(huán)境問題主要是廢水廢渣的排放,其次是地裂縫等。因此本區(qū)的防控建議,是對廢水、廢液進(jìn)行凈化治理,有污水排放的礦山企業(yè)都建起坑口沉淀池和污水處理設(shè)施,做到達(dá)標(biāo)排放；對于尾礦進(jìn)行專業(yè)的開發(fā)利用和回收處理。

低風(fēng)險區(qū)：建議在維持現(xiàn)狀的基礎(chǔ)上,在今后礦業(yè)開發(fā)活動中注重礦山地質(zhì)環(huán)境的保護(hù)。

5 結(jié)論

(1)在以往礦山地質(zhì)環(huán)境綜合評價的基礎(chǔ)上,借鑒自然災(zāi)害風(fēng)險評價理論,提出了礦山地質(zhì)環(huán)境風(fēng)險評價的新思路和方法。認(rèn)為礦山地質(zhì)環(huán)境風(fēng)險性是風(fēng)險施加者礦業(yè)開發(fā)影響力和風(fēng)險承災(zāi)體地質(zhì)環(huán)境脆弱性共同作用形成的。

(2)大同市礦山地質(zhì)環(huán)境風(fēng)險評價結(jié)果分為高風(fēng)險區(qū)、中風(fēng)險區(qū)、低風(fēng)險區(qū),面積分別為1 277.26 km2、1 505.49 km2、11 393.25 km2,分別占大同市總面積的9.01%,10.62%,80.37%。通過驗證,該評價結(jié)果與實際情況基本相符,評價結(jié)果有效。

(3)為了達(dá)到防災(zāi)減災(zāi)效果,提出了大同市風(fēng)險防控建議,使評價結(jié)果能更好地為礦產(chǎn)資源開發(fā)利用和保護(hù)服務(wù)。

參考文獻(xiàn)：

[1] 唐朝暉, 劉楠, 柴波,等. 合山市礦山地質(zhì)環(huán)境影響評價研究[J]. 水文地質(zhì)工程地質(zhì), 2012, 39(6): 124-129. [ TANG Z H, LIU N, CHAI B,etal. Mining geo-environmental impact assessment of Heshan City[J]. Hydrogeology & Engineering Geology, 2012, 39(6):124-129. (in Chinese)]

[2] 梁偉. 基于RS與GIS的鶴崗礦山地質(zhì)環(huán)境評價研究[D]. 長春:吉林大學(xué), 2012. [ LIANG W. The research of mine geological environment assessment of Hegang based on RS and GIS[D]. Changchun: Jilin University, 2013. (in Chinese)]

[3] 張建萍, 曹金亮. 山西省礦山環(huán)境地質(zhì)問題及其分區(qū)研究[J]. 工程地質(zhì)學(xué)報, 2007, 15(2): 25-34. [ ZHANG J P, CAO J L. Geo-environmental problems in mines in Shanxi and a study on effect zone[J]. Journal of Engineering Geology, 2007, 15(2): 25-34. (in Chinese)]

[4] 常青, 邱瑤, 謝苗苗，等. 基于土地破壞的礦區(qū)生態(tài)風(fēng)險評價[J]. 生態(tài)學(xué)報, 2012, 32(16): 5164-5174. [ CHANG Q, QIU Y, XIE M M,etal. Theory and method of ecological risk assessment for mining areas based on the land destruction[J]. Acta Ecologica Sinica, 2012, 32(16): 5164-5174. (in Chinese)]

[5] 匡恒, 馬傳明, 和澤康. 基于城市規(guī)劃的開封市地質(zhì)環(huán)境風(fēng)險評價[J]. 監(jiān)測與評價, 2016(2): 133-137. [ KUANG H, MA C M, HE Z K,etal. The geo-environmental risk assessment of KaiFeng city based on the urban planning[J]. Environmental Monitoring (Assessment, 2016, (2): 133-137. (in Chinese)]

[6] 劉德玉. 青海省木里煤田江倉礦區(qū)地質(zhì)生態(tài)環(huán)境風(fēng)險評價[D]. 北京:中國地質(zhì)科學(xué)院, 2013. [LIU D Y. Risk assessment of geo-ecological environment in Jiangcang diggings Muli coal field, Qinghai Province[D]. Beijing: Chinese Academy of Geological Sciences, 2013. (in Chinese)]

[7] 羅炳佳,沈誠.貴州礦山地質(zhì)環(huán)境影響評估[J].水文地質(zhì)工程地質(zhì),2013,40(1):134-138.[LUO B J,SHEN C.Impact assessment of mine geological environment of Guizhou Province[J].Hydrogeology & Engineering Geology,2013,40(1):134-138.(in Chinese)]

[8] DD2014—05 礦山地質(zhì)環(huán)境調(diào)查評價規(guī)范[S]. [DD2014—05 The specification for investigation and evaluation of mine geological environment[S]. (in Chinese)]

[9] 孟慶華. 秦嶺山區(qū)地質(zhì)災(zāi)害風(fēng)險評估方法研究[D]. 北京:中國地質(zhì)科學(xué)院, 2011. [ MENG Q H. Study on the methods of geo-hazards risk assessment in Qinling Mountain[D]. Beijing: Chinese Academy of Geological Sciences, 2011. (in Chinese)]

[10] 向喜瓊, 黃潤秋. 地質(zhì)災(zāi)害風(fēng)險評價與風(fēng)險管理[J]. 地質(zhì)災(zāi)害與環(huán)境保護(hù), 2000, 11(1): 92-96. [XIANG X X, HUANG R Q. Risk assessment and risk management for slope geohazards[J]. Journal of Geological Hazards and Environment Preservation, 2000, 11(1): 92-96. (in Chinese)]