婁新剛，康志勇，劉 超

(遼寧省冶金地質(zhì)勘查局地質(zhì)勘查研究院，遼寧 鞍山 114038)

遼寧省鞍山市是我國重要的工業(yè)城市，素有鋼都之稱，城區(qū)人口146萬?，F(xiàn)已發(fā)展成為以冶金工業(yè)為主、機械、化工、電子紡織等工業(yè)門類較齊全的工業(yè)城市，隨著改革開放以來城市經(jīng)濟發(fā)展，旅游等第三產(chǎn)業(yè)得到了飛速發(fā)展。城市建設(shè)日新月異，已初見現(xiàn)代化大都市的雛形。

鞍山市區(qū)的地質(zhì)環(huán)境復雜，地質(zhì)災害發(fā)育。城區(qū)的周邊分布著多處以鐵礦為主礦山，截止到2009年底，全市有各類礦山企業(yè)70家(生產(chǎn)17家、閉坑礦山企業(yè)53家)。這些企業(yè)構(gòu)成了鞍山市的產(chǎn)業(yè)基礎(chǔ)，為鞍山市的經(jīng)濟建設(shè)做出了巨大的貢獻，但同時也帶來了較多的環(huán)境問題，如破壞或占用大量土地，引起土地荒漠化、水土流失、生態(tài)環(huán)境退化。一些礦區(qū)還發(fā)生了滑坡、崩塌、地面塌陷等地質(zhì)災害，給人民生命財產(chǎn)造成威脅，制約經(jīng)濟的發(fā)展。

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

1.1 礦業(yè)開發(fā)對土地資源的占用及破壞

由于鞍山市是礦業(yè)型城市，礦產(chǎn)資源開發(fā)利用歷史悠久，礦山數(shù)量較多，開發(fā)規(guī)模較大。礦產(chǎn)資源開發(fā)利用在為國家經(jīng)濟建設(shè)作出巨大貢獻的同時，也給鞍山地區(qū)土地和生態(tài)環(huán)境帶來了巨大壓力和破壞。以四大鐵礦山(西鞍山未列)統(tǒng)計數(shù)據(jù)為例：現(xiàn)有大型采坑4個，占地11.6km2；5個大型排巖場，堆積巖渣總量超過16億t，占地22.3km2；6個大型尾礦庫，堆積尾礦約5億t，占地7.3km2。這些采坑、排巖場、尾礦庫總面積達41.2km2，約占市區(qū)(含市郊)總面積的6.6%。其中東鞍山鐵礦占地面積1.75km2，有排石場一個，占地面積4.5km2，尾礦庫3個，占地面積2.46km2；大孤山鐵礦占地2.38km2，有排石場一個，占地面積7.54km2，尾礦庫一個，占地面積2.45km2；齊大山鐵礦占地4.81km2，排石場一個，占地4.41km2，尾礦庫兩個，占地2.39km2；眼前山鐵礦占地2.66km2，排石場一個占地5.84km2，尾礦庫與大孤山礦山共用。在一些采坑周圍、排巖場區(qū)、尾礦庫區(qū)現(xiàn)已成裸地，植被度在2度以下，荒漠化十分嚴重。

1.2 礦業(yè)開發(fā)對地下水系統(tǒng)的影響

礦山開采時進行疏干排水，導致水均衡系統(tǒng)破壞，造成地表河流干涸斷流，地下水位下降形成降落漏斗。開口面積1～2km以至幾平方千米，降深達200m左右的四大采坑。由于生產(chǎn)所需人工疏干，在千山這一系統(tǒng)的周圍形成幾個巨大的地下水漏斗，每個漏斗的影響面積都在十幾平方公里左右，礦山周圍的民用井大部分都已廢棄。影響范圍直接深入千山風景區(qū)的核心部位，地下水位的嚴重下降，直接影響千山風景區(qū)的植被生長。

1.3 形成一個礦區(qū)生產(chǎn)活動帶

已有的地質(zhì)災害大都發(fā)生這些大型鐵礦的排石場、尾礦庫、露天邊坡，此外，部分日偽時期留下的采空區(qū)發(fā)生過地面塌陷。

2 地質(zhì)環(huán)境質(zhì)量評價

2.1 空間數(shù)據(jù)庫的建立

收集研究區(qū)內(nèi)的所有地質(zhì)環(huán)境資料，并對資料篩選分類，數(shù)據(jù)按照實體描述類型分為空問數(shù)據(jù)和屬性數(shù)據(jù)兩類。采用 MAPGIS對研究區(qū)基礎(chǔ)資料進行數(shù)字化處理，然后對基礎(chǔ)圖件進行分離圖層，進行分層式管理。按照空間數(shù)據(jù)各要素的特點，將空問數(shù)據(jù)分解為單一性質(zhì)的基本要素圖層，各圖層主要為略陽縣基本地理圖層、地形等高線圖層、地貌圖層、地層巖性圖層、地質(zhì)構(gòu)造圖層、礦產(chǎn)分布圖層等。同時，對需要進行投影變換、圖幅拼接的地圖數(shù)據(jù)，進行投影變換、圖幅拼接，以實現(xiàn)不同比例尺地圖的合成，以及跨分度帶地圖的拼接。屬性數(shù)據(jù)包括對圖形進行描述或說明的數(shù)據(jù)及用于進行評價的基礎(chǔ)數(shù)據(jù)，如高程信息、地理標注、地層巖性特征、植被發(fā)育情況等。屬性數(shù)據(jù)通過Access數(shù)據(jù)庫進行錄入、管理，MAPGIS通過接口進行調(diào)用。

2.2 評價指標體系

在對資料進行分析后 ，通過兩兩比較法對指標進行篩選優(yōu)化，選取地形坡度、工程地質(zhì)巖組、巖土體結(jié)構(gòu)類型、地下水位、植被發(fā)育情況、月平均降雨量 、人類工程活動強度等為主要評價指標 。地質(zhì)環(huán)境質(zhì)量分為優(yōu)、良、中、差 4個等級，相應地將評價指標也用4級來描述，建立評價指標體系(表1)。

表1 地質(zhì)環(huán)境質(zhì)量評價指標分級標準

計算中，采用層次分析法來確定評價指標的權(quán)重，具體步驟如下：

1)構(gòu)造判斷矩陣。對于z1，z2，…，zn個評價指標，由專家組運用兩兩比較得判斷矩陣Z。

2)將判斷矩陣按行求

(1)

3)歸一化

(2)

4)求最大特征值

(3)

式中:P=(w1w2……wn)T。

5)進行一致性和隨機性檢驗

C.I=(λmax-n)/(n-1)

(4)

C.R=C.I/R.I

(5)

式中：C.I為一致性指標；n為矩陣階數(shù)；R.I為平均隨機一致性指標 ；C.R為隨機一致性比率。只有當C.R<0.10時，求出的權(quán)值才比較合理。

通過計算 ，得出各指標的權(quán)重如表2。

表2 評價指標權(quán)重

2.3 評價模型

評價數(shù)學模型采用由陳守煜提出的多級模糊模式識別模型，其以相對級別特征值作為判斷或識別的依據(jù)，克服了最大隸屬度原則所不適用的地方，而且以相對隸屬度、隸屬函數(shù)為基礎(chǔ)理論，使隸屬度、隸屬函數(shù)的計算更容易。實現(xiàn)步驟為：

1)規(guī)格化處理。設(shè)某樣本需要對模糊子集 A進行識別，樣本的m個指標組成的特征向量為

X=(x1x2…xixm)T

式中：xi為樣本指標i的特征值；i=1，2，…，m。

如樣本依據(jù)m個指標c個級別的已知指標標準特征值進行識別，則有指標標準特征值矩陣

式中：Yih為指標i級別h的標準特征值；h=1，2…，c；i=1，2，…，m。

則遞減型指標對A的相對隸屬度

(6)

h級指標標準特征值對A的相對隸屬度

(7)

遞增型指標對A的相對隸屬度

(8)

h級指標標準特征值，對 A的相對隸屬度

(9)

2)求最優(yōu)相對隸屬度向量。樣本對A的各個級別相對隸屬度矩陣

U=(u1u2…un)T=uh

式中 ：uh為樣本對A的級別h的相對隸屬度，按下式求解

(10)

式中：a、b為樣本級別區(qū)間范圍；wi為樣本指標i的權(quán)重；p為距離參數(shù) ，p= 1為海明距離，p=2為歐氏距離；dh為廣義權(quán)距離。

3)求級別特征值。相對級別特征值為級別變量h以對應的相對隸屬度為權(quán)重的總和，即

(11)

式中：Hu為相對級別特征值，是描述狀態(tài)或級別的無量綱數(shù)，且1≤Hu≤c。由于其利用了全部相對隸屬度信息，其判斷結(jié)論比最大隸屬度原則更符合實際情況。因此，可應用相對級別特征值來判定此樣本屬于何種級別。

2.4 評價結(jié)果分析

用 C++Builder 6．O開發(fā)的評價模型，對鞍山市區(qū)進行了地質(zhì)環(huán)境的質(zhì)量評價，將評價分為四大區(qū)八個亞區(qū)。分區(qū)結(jié)果見鞍山市區(qū)地質(zhì)環(huán)境質(zhì)量分區(qū)表(表3)。

表3 鞍山市區(qū)地質(zhì)環(huán)境質(zhì)量分區(qū)表

續(xù)表

3 結(jié)語

運用 GIS技術(shù)，建立了鞍山市區(qū)空間數(shù)據(jù)庫，綜合考慮地形坡度、工程地質(zhì)巖組 、巖土體結(jié)構(gòu)類型地下水位 、植被發(fā)育的情況、月平均降雨量、人類工程活動強度等指標 ，運用層次分析法確定各權(quán)值；采用多級模糊模式識別模型，對鞍山市區(qū)進行了地質(zhì)環(huán)境質(zhì)量評價，把鞍山市區(qū)地質(zhì)環(huán)境劃分為優(yōu)、良、中、差4個級別。評價結(jié)果，為鞍山市區(qū)地質(zhì)環(huán)境的開發(fā)、利用和保護以及地質(zhì)災害的預防，提供了一定的參考依據(jù)。

[1]王景明，王春梅，劉科.地裂縫及其災害研究的新進展[J].地球科學進展，2001，16(3)：303-313.

[2]田大佑.地質(zhì)災害防治與制度建設(shè)[J].湖北地礦，2000，14(3-4)：64-66.

[3]張開平.滑坡泥石流治理探討[J].水土保持研究，2000，7(3)：174-176.

[4]鹿愛莉，曹清華，薛全全.建立地質(zhì)災害調(diào)查項目社會評價指標體系[J].資源·產(chǎn)業(yè)，2003，5(4)：16-18.

[5]柳源.建立具有中國特色的地質(zhì)災害防災預警系統(tǒng)[J].資源·產(chǎn)業(yè)，2003，5(2)：56-57.

[6]袁紹武.塌陷防治與可持續(xù)發(fā)展[J].中國地質(zhì)災害與防治學報，1997，8(增刊)，29-33.

[7]李毅，李蘅，張靜.我國礦山地質(zhì)災害主要類型和勘查防治方法[J].礦產(chǎn)與地質(zhì)，2004，18(1)：62-64.

[8]楊裕云，楊紅剛，吳有才.與地下水作用有關(guān)的地質(zhì)災害[J].水文地質(zhì)工程地質(zhì)，2004(增刊)：1-7.

[9]康志勇.鞍山市“九五”時期地礦成果匯編(1996～2000年)[C].鞍山，遼寧冶金地質(zhì)勘查局，2001.