基于模糊數(shù)學綜合評價法的礦山地質(zhì)環(huán)境影響評價

尹俊凱 汪繼學 張承斌

（山東泰山地質(zhì)勘查公司，山東 泰安 271000）

日照市共開發(fā)利用鐵、銅、鉛等21種礦產(chǎn)，有各類礦山546個，多數(shù)礦山企業(yè)規(guī)模屬中小型。主要開采礦種為飾面用花崗巖、建筑石料用花崗巖、石灰?guī)r（建筑石料用灰?guī)r和水泥用灰?guī)r）等，金屬礦產(chǎn)主要為金礦、鐵礦、鉛鋅礦，其他還有冶金用白云巖、熔劑用蛇紋巖、磚瓦用頁巖、砂巖、鉀長石、粘土礦等非金屬礦產(chǎn)。日照市最主要的礦山地質(zhì)環(huán)境問題是破損山體及露天采坑，對原有地形地貌造成很大破壞并占用大量土地，由于露天開采造成地表植被破壞、水土流失等，使生態(tài)環(huán)境惡化，個別處于“三區(qū)兩線”可視范圍內(nèi)，造成嚴重視覺污染。

日照市礦山產(chǎn)生地質(zhì)環(huán)境問題數(shù)量多，種類雜，且分布區(qū)域廣，為更好地認識礦山開采對地質(zhì)環(huán)境的影響，進一步指導礦山地質(zhì)環(huán)境保護與恢復治理，有必要對礦山地質(zhì)環(huán)境影響進行評價分區(qū)。

1 評價原則

為了規(guī)范分區(qū)，全面反映礦山地質(zhì)環(huán)境影響等級，分區(qū)評價遵循以下原則。

1.1 堅持“以人為本”的原則

選取與人民群眾生活息息相關(guān)的占損土地面積、地質(zhì)災(zāi)害作為主要評價因子，評價其影響程度。當存在多種問題時，或同一問題分散分布時，對人們的影響大于單一問題集中分布，因此增加問題數(shù)量評價因子。以可視距離1km作為礦山的影響范圍邊界，分區(qū)以此為界。

1.2 規(guī)避“三區(qū)兩線”的原則

以“三區(qū)兩線”作為天然的分區(qū)邊界線?！叭齾^(qū)”指城市規(guī)劃區(qū)、自然保護區(qū)、風景名勝區(qū)，“兩線”為重要交通干線和主要水系線，重要交通干線主要考慮高速公路、國道和鐵路，水系線主要考慮河流線和水庫，海岸線、湖泊為當然的分區(qū)邊界線。當“兩線”穿越礦區(qū)時，不再作為天然分界線。

1.3 集中連片的原則

多個礦山集中分布（大于2個），單個礦山相距小于2km，作為一個集中連片區(qū)進行評價；當兩個集中連片區(qū)相距小于2倍可視范圍，且主要礦山地質(zhì)環(huán)境問題相似時，將兩個集中連片區(qū)合并為同一個區(qū)進行評價；單一礦山獨立分布時，且面積小于1km2時，不進行分區(qū)評價，將其歸為一般區(qū)。

1.4 突出主要礦山地質(zhì)環(huán)境問題的原則

當兩個集中連片區(qū)相距小于2倍可視范圍，但主要礦山地質(zhì)環(huán)境問題不同，將兩個集中連片區(qū)分別進行評價，不再合并，以突出主要礦山地質(zhì)環(huán)境問題。

2 評價方法

2.1 評價因子

本次礦山地質(zhì)環(huán)境評價，選擇3個因子，分別是土地占損率A、礦山地質(zhì)環(huán)境問題點分布密度B、地質(zhì)災(zāi)害C。通過對這3個因子進行量化賦值，進行評價分區(qū)。

土地占損率=占損土地面積/分區(qū)面積×100，只統(tǒng)計已有占損（反映現(xiàn)狀）。

礦山地質(zhì)環(huán)境問題點分布密度=區(qū)內(nèi)礦山地質(zhì)環(huán)境問題點數(shù)量/面積×100，即100km2的礦山地質(zhì)環(huán)境問題數(shù)量，便于不同面積分區(qū)之間的對比。（只統(tǒng)計已存在的）

地質(zhì)災(zāi)害分為兩個部分，財產(chǎn)損失和災(zāi)害點數(shù)量，隱患點的影響與已發(fā)生的災(zāi)害點同等重要，險情威脅損失同直接經(jīng)濟損失。該因子同第2個因子，以每百平方千米進行統(tǒng)計。財產(chǎn)損失和災(zāi)害點數(shù)量的權(quán)重系數(shù)按相同處理。

2.2 分區(qū)

（1）將礦山中心點坐標和礦區(qū)范圍投影到分區(qū)圖上，未調(diào)查到礦區(qū)范圍的，以現(xiàn)狀影響范圍為礦區(qū)范圍。

（2）將圖按2km×2km進行網(wǎng)格劃分，查看每個網(wǎng)格的礦山分布情況，將相鄰網(wǎng)格分布有礦山且不屬于同一礦山的劃為同一區(qū)，當相鄰網(wǎng)格無礦山，次一相鄰有礦山的也劃分為同一區(qū)。連續(xù)兩個相鄰網(wǎng)格無礦山則歸為不同的區(qū)。將大量無礦山的地區(qū)歸為一般區(qū)。當相鄰網(wǎng)格有礦山且屬同一礦山的，礦區(qū)面積可局限于一個網(wǎng)格之內(nèi)，且連續(xù)兩個相鄰網(wǎng)格無礦山的，做為一個礦山點分布于圖上，歸為一般區(qū)進行分區(qū)評價，在分區(qū)評價完成后，進行單點評價，因子同上，面積以1個單元網(wǎng)格面積計算。

（3）分區(qū)邊界，以可視范圍1km做為分區(qū)邊界，穿過村莊、城鎮(zhèn)、湖泊等不宜穿過的地區(qū)時，向內(nèi)側(cè)偏移盡可能避過。

圖1 礦山地質(zhì)環(huán)境評價分區(qū)示意圖

2.3 分區(qū)評價方法

2.3.1 土地占損率模糊判別

土地占損率

式中：

ai-第i區(qū)土地占損率；

Fij-第i區(qū)第j項礦山地質(zhì)環(huán)境問題占損土地面積（km2）；

Si-第i區(qū)面積（km2）。

以函數(shù)2進行模糊判別

該方程表達的主要思想是以土地占損大于20歸為嚴重，其他小于20的按該公式判別其模糊歸屬度。

2.3.2 礦山地質(zhì)環(huán)境問題點分布密度模糊判別

礦山地質(zhì)環(huán)境問題點分布密度

式中：

bi-第i區(qū)礦山地質(zhì)環(huán)境問題點分布密度(個/百平方千米)；

n-第i區(qū)礦山地質(zhì)環(huán)境問題點數(shù)量。

以函數(shù)4進行模糊判別

2.3.3 地質(zhì)災(zāi)害

地質(zhì)災(zāi)害判別因子C分為單位財產(chǎn)損失cc和災(zāi)害點分布密度cm兩部分，地質(zhì)災(zāi)害隱患點的影響與已發(fā)生的災(zāi)害點同等重要，險情威脅損失同直接經(jīng)濟損失。同樣按百平方千米進行統(tǒng)計。

單位面積財產(chǎn)損失

災(zāi)害點分布密度

式中：

CCSij-第i區(qū)第j個地質(zhì)災(zāi)害財產(chǎn)損失；

m-第i區(qū)地質(zhì)災(zāi)害點數(shù)量。

參考地質(zhì)災(zāi)害危險性評估規(guī)范（DZ/T 0286-2015）表—地質(zhì)災(zāi)害危害程度分級表，以100萬和500萬作為分級標準，≥500萬元的危害程度大，≤100萬元的地質(zhì)災(zāi)害危害程度小。本文以單位面積財產(chǎn)損失≥500萬元/km2嚴重歸屬度為1，≤100萬元/km2嚴重歸屬度為0.1，建立模糊判別函數(shù)（7）如下：

本文以5個小型災(zāi)害可比一個大型災(zāi)害，為增加其收斂性，將災(zāi)害點分布密度擴大10倍，建立災(zāi)害點分布密度模糊判別函數(shù)（8）。

2.3.4 建立模糊判別矩陣Rij

根據(jù)上述函數(shù)，可分別計算其礦山地質(zhì)環(huán)境影響嚴重程度歸屬度，建立模糊判別矩陣Rij。

式中：

2.3.5 截矩陣進行判別

引入截矩陣進行礦山地質(zhì)環(huán)境影響嚴重程度判別。

取λ=0.9，0.5，0.1分別判別，取R0.9（rij≥0.9）歸為嚴重區(qū)，R0.5不含R0.9（0.5≤rij＜0.9）歸為較嚴重區(qū)，R0.1不含R0.5（0.1≤rij＜0.5）歸為輕微區(qū)，R0.1（rij＜0.1）歸為一般區(qū)。

2.3.6 分區(qū)賦值計算

根據(jù)上述評價原則及方法對日照市礦山地質(zhì)環(huán)境進行分區(qū)評價，根據(jù)調(diào)查數(shù)據(jù)進行賦值計算，通過2.3.5所列截矩陣進行礦山地質(zhì)環(huán)境影響嚴重程度判別。各典型分區(qū)計算結(jié)果見表1。

表1 典型分區(qū)計算結(jié)果表

3 評價結(jié)果

根據(jù)日照市礦山地質(zhì)環(huán)境調(diào)查數(shù)據(jù)，對日照市礦山地質(zhì)環(huán)境影響程度現(xiàn)狀進行了評價分區(qū)，將日照市劃分為40個小區(qū)，其中礦山地質(zhì)環(huán)境影響嚴重區(qū)12個，礦山地質(zhì)環(huán)境影響較嚴重區(qū)16個，礦山地質(zhì)環(huán)境影響輕微區(qū)11個，礦山地質(zhì)環(huán)境影響一般區(qū)（即空白區(qū)）1個（見圖2）。

3.1 礦山地質(zhì)環(huán)境影響嚴重區(qū)

本次調(diào)查劃分礦山地質(zhì)環(huán)境影響嚴重區(qū)12個，面積299.88km2，共包括298個礦山，473處礦山地質(zhì)環(huán)境問題，占損土地面積約930公頃。

3.2 礦山地質(zhì)環(huán)境影響較嚴重區(qū)

本次調(diào)查劃分礦山地質(zhì)環(huán)境影響較嚴重區(qū)16個，面積440.73km2，共包含礦山192個，245處地質(zhì)環(huán)境問題，占損土地面積約為660公頃。

3.3 礦山地質(zhì)環(huán)境影響輕微區(qū)

本次調(diào)查劃分礦山地質(zhì)環(huán)境影響輕微區(qū)11個，面積192.89km2，共包括41個礦山，50處礦山地質(zhì)環(huán)境問題，占損土地面積約180公頃。

圖2 日照市礦山地質(zhì)環(huán)境影響評價分區(qū)圖

對比日照市礦山地質(zhì)環(huán)境實際情況，利用基于模糊數(shù)學的綜合評價法對礦山地質(zhì)環(huán)境影響進行的評價分區(qū)，其應(yīng)用方便，結(jié)果合理，對下一步的日照市礦山地質(zhì)環(huán)境保護與治理規(guī)劃有重要的指導與借鑒意義。