龔華根，徐 芬，楊 輝

(1.江西省地礦資源勘查開發(fā)中心，江西 南昌 330030;2.江西省勘察設(shè)計(jì)研究院，江西 南昌 330030)

礦山巖溶塌陷預(yù)測(cè)評(píng)估作為礦山環(huán)境評(píng)價(jià)中重要的單環(huán)境問題評(píng)價(jià)之一，是礦山環(huán)境問題研究的一個(gè)重要環(huán)節(jié)[1]，如何科學(xué)合理的對(duì)礦山巖溶塌陷作出評(píng)價(jià)也是礦山環(huán)境研究的一個(gè)重大課題，也是礦山生態(tài)環(huán)境保護(hù)與恢復(fù)治理規(guī)劃的基礎(chǔ)[2]。

礦山巖溶塌陷預(yù)測(cè)評(píng)估一直是礦山地質(zhì)環(huán)境保護(hù)與恢復(fù)治理中的重難點(diǎn)問題。目前，礦山地質(zhì)環(huán)境保護(hù)與恢復(fù)治理中的礦山巖溶塌陷預(yù)測(cè)評(píng)估大部分采用定性描述或根據(jù)量化評(píng)價(jià)標(biāo)準(zhǔn)表進(jìn)行半定量打分進(jìn)行評(píng)價(jià)。然而，先進(jìn)的GIS技術(shù)與科學(xué)的預(yù)測(cè)模型相結(jié)合，是對(duì)巖溶塌陷等進(jìn)行定量預(yù)測(cè)的一種非常行之有效的辦法[3]。因此，本文用此方法對(duì)某大理石礦巖溶塌陷進(jìn)行預(yù)測(cè)評(píng)估。

1 礦區(qū)地質(zhì)環(huán)境背景

礦區(qū)位于江西省中東部，礦種為大理石礦;開采方式為露天開采;礦區(qū)面積1.970 2 km2;開采標(biāo)高為 +280～+508 m。礦區(qū)屬亞熱帶東南季風(fēng)性氣候區(qū)，氣候溫和，四季分明，歷年平均降雨量為1 728.5 mm。礦區(qū)地形變化不大，屬淺切割的低山地形巖溶地貌區(qū)(見圖1)。地形南低北高，中部高，東部及西部低。最高海拔 508.91 m，最低海拔 272.30 m[4]。

礦區(qū)內(nèi)地層出露較簡(jiǎn)單，具體可見圖2(三維地質(zhì)圖)。主要為藤田盆地發(fā)育的石炭系下統(tǒng)上部為大塘組(C1d)、石炭系中統(tǒng)黃龍組(C2h)、上統(tǒng)船山組(C3c)和第四系(Q)地層。礦區(qū)以單斜巖層為特征，地層總體走向?yàn)槟媳毕?，傾向250°～270°，傾角一般為 25°～35°之間。礦區(qū)內(nèi)斷裂構(gòu)造不發(fā)育，未發(fā)現(xiàn)有斷層分布[4]。

2 礦山巖溶塌陷機(jī)理分析

巖溶塌陷產(chǎn)生與否主要由自然因素所決定[5]，這些因素主要有覆蓋層特征、巖溶發(fā)育程度和地下水活動(dòng)。它們共同構(gòu)成了影響巖溶塌陷的基本因素，也就是巖溶塌陷的內(nèi)因，只要具備了塌陷的內(nèi)因，巖溶塌陷就可能形成。而自然因素和人類活動(dòng)的影響，主要起加速和引發(fā)作用，它們共同構(gòu)成了影響巖溶塌陷的外因。巖溶塌陷機(jī)理分析見圖3。

圖1 地形地貌曲面圖

圖2 三維地質(zhì)模型柵欄圖

圖3 巖溶塌陷機(jī)理分析

礦山開采導(dǎo)致巖溶塌陷的外因又具有其專屬特點(diǎn)。礦山基建及開采過程中可能促進(jìn)巖溶塌陷形成的活動(dòng)主要有:抽排地下水、施工機(jī)械及車輛的荷載與振動(dòng)、挖填土層使覆蓋層變薄或增厚。礦山巖溶塌陷的發(fā)生主要是在覆蓋型巖溶礦區(qū)淺部巖溶發(fā)育、第四系厚度較小的地段，當(dāng)進(jìn)行抽排水、荷載、振動(dòng)等人類活動(dòng)時(shí)，或改變了地下水的運(yùn)動(dòng)狀態(tài)，水力坡度增大，使地下水潛蝕、淘空、搬運(yùn)第四系土體及巖溶充填物的能力增強(qiáng)，或改變地下巖溶頂板應(yīng)力狀態(tài)，引起自然條件下的相對(duì)應(yīng)力平衡遭到破壞，因而產(chǎn)生巖溶塌陷。

3 層次分析法確定權(quán)重

3.1 建立層次結(jié)構(gòu)

通過對(duì)礦山巖溶塌陷的機(jī)理分析，選取地下水活動(dòng)、覆蓋層特征、巖溶發(fā)育程度和人類活動(dòng)強(qiáng)度這四個(gè)方面，共6個(gè)評(píng)價(jià)指標(biāo)，建立礦山巖溶塌陷預(yù)測(cè)評(píng)估層次結(jié)構(gòu)。見圖4。

圖4 層次分析結(jié)構(gòu)模型

3.2 構(gòu)造判斷矩陣

根據(jù)表 1 重要性標(biāo)度含義表，對(duì) C1、C2、C3、C4、C5、C6六個(gè)因子進(jìn)行兩兩比較，得出重要性標(biāo)度含義表，通過專家咨詢、篩選、打分，確定每一層各個(gè)因子之間的相對(duì)重要性比值，并構(gòu)造判斷矩陣，見表2。

表1 重要性標(biāo)度含義表

表2 判斷矩陣

3.3 計(jì)算權(quán)向量

通過計(jì)算判斷矩陣的權(quán)向量得出每一個(gè)指標(biāo)及每一層次各個(gè)因子的權(quán)重。計(jì)算權(quán)向量有特征根法、和法、根法、冪法等，這里采用較為簡(jiǎn)便的和法。

3.4 一致性檢驗(yàn)

表3 平均隨機(jī)一致性指標(biāo) R.I.表 (1000次正互反矩陣計(jì)算結(jié)果)

經(jīng)計(jì)算，λmax≈6.1，C.R.=0.016 ＜ 0.1。因此，認(rèn)為判斷矩陣的一致性是可以接受的。

各因子 權(quán) 重 為:Wi=[0.077，0.037，0.283，0.037，0.283，0.283]T。

圖5 最低地下水位分區(qū)圖

圖6 巖溶率分區(qū)圖

圖7 下伏巖體巖性分區(qū)圖

4 對(duì)各評(píng)價(jià)因子進(jìn)行等級(jí)劃分

依據(jù)層次分析法數(shù)學(xué)模型，根據(jù)專家評(píng)定及大量實(shí)例統(tǒng)計(jì)相結(jié)合，將礦區(qū)巖溶塌陷各評(píng)價(jià)因子進(jìn)行量化分級(jí)。詳見表4。

表4 巖溶塌陷各評(píng)價(jià)因子等級(jí)劃分表

5 基于GIS對(duì)評(píng)估區(qū)內(nèi)各評(píng)價(jià)因子進(jìn)行分區(qū)

最低地下水水位分區(qū):根據(jù)地形地貌及第四系覆蓋情況，位于溝谷地形較為平緩的第四系覆蓋區(qū)最低地下水水位一般位于基巖頂面上下，而地形較陡的山坡地區(qū)最低地下水位一般位于基巖頂面以下，賦存于地下巖溶通道、暗河及裂隙中。

巖溶率分區(qū):根據(jù)儲(chǔ)量核實(shí)報(bào)告中的鉆孔數(shù)據(jù)、物探數(shù)據(jù)及地表溶蝕洼地、溶洞、落水洞等的分布進(jìn)行圈定。

下伏巖體巖性分區(qū):根據(jù)鉆孔揭露的地層及區(qū)域地質(zhì)構(gòu)造進(jìn)行劃分。

第四系覆蓋層厚度分區(qū):根據(jù)鉆孔揭露的第四系厚度及現(xiàn)場(chǎng)踏勘的結(jié)果進(jìn)行劃分。

抽排水影響范圍分區(qū):根據(jù)開采最低標(biāo)高+280 m及含水層的展布、地下水流向進(jìn)行劃定，低于+280 m的區(qū)域及評(píng)估區(qū)東北側(cè)大塘組隔水層等區(qū)域?yàn)橛绊戄^小區(qū)，再根據(jù)距離抽排水區(qū)域的距離進(jìn)一步劃分為500 m以內(nèi)及500 m以外區(qū)域。

振動(dòng)及荷載分區(qū):根據(jù)礦山開發(fā)利用方案設(shè)計(jì)的開采范圍及廢石場(chǎng)、荒料場(chǎng)場(chǎng)址劃定。

評(píng)估區(qū)內(nèi)各評(píng)價(jià)因子分區(qū)結(jié)果見圖5～圖10。

圖8 第四系覆蓋層厚度分區(qū)圖

圖9 抽排水影響范圍分區(qū)圖

圖10 振動(dòng)及荷載分區(qū)圖

6 對(duì)各評(píng)價(jià)因子分區(qū)圖進(jìn)行空間分析

利用MAPGIS空間分析功能，對(duì)評(píng)估區(qū)內(nèi)各評(píng)價(jià)因子分區(qū)結(jié)果進(jìn)行復(fù)合、求交、疊加運(yùn)算，并對(duì)各評(píng)價(jià)因子的屬性字段進(jìn)行計(jì)算，確定各預(yù)測(cè)單元的隸屬函數(shù)矩陣。按照隸屬度最大原則，最大值則為該預(yù)測(cè)評(píng)判單元所對(duì)應(yīng)的級(jí)別。

以上通過空間分析及編制程序自動(dòng)完成，并將評(píng)判結(jié)果存到各預(yù)測(cè)單元的屬性庫中，最后根據(jù)計(jì)算結(jié)果，把相同隸屬度的單元?jiǎng)澐譃橥患?jí)別，不同級(jí)別賦予不同顏色，得到某大理石礦巖溶塌陷易發(fā)性分區(qū)圖。分區(qū)結(jié)果見圖11。

圖11 巖溶塌陷易發(fā)性分區(qū)圖

7 結(jié)語

通過對(duì)上述預(yù)測(cè)評(píng)估結(jié)果的分析，并結(jié)合評(píng)估區(qū)地質(zhì)環(huán)境條件及巖溶塌陷現(xiàn)狀，可以確定:在自然條件下，礦區(qū)露采范圍、荒料場(chǎng)內(nèi)、北露采區(qū)北側(cè)和東側(cè)及南露采區(qū)南部的部分第四系覆蓋較薄區(qū)域，由于受地下水位波動(dòng)及礦山露采抽排水的影響，屬于巖溶塌陷易發(fā)區(qū);礦區(qū)范圍周邊其他受礦山露采抽排水和振動(dòng)、荷載影響的區(qū)域，由于下伏巖體存在可溶巖體的分布，屬于巖溶塌陷較易發(fā)區(qū);其他區(qū)域由于受礦山開采影響較小、巖溶發(fā)育程度較低、局部地區(qū)不存在可溶巖體的分布，屬于巖溶塌陷不易發(fā)區(qū)。

所建立的預(yù)測(cè)評(píng)估模型能反映礦山巖溶塌陷形成的一般規(guī)律，模型中的評(píng)價(jià)因子基本符合某大理石礦巖溶塌陷的實(shí)際情況，現(xiàn)場(chǎng)調(diào)查發(fā)現(xiàn)的兩個(gè)已有巖溶塌陷點(diǎn)均位于巖溶塌陷的易發(fā)區(qū)。但評(píng)價(jià)因子的選取及級(jí)別的劃分還有待今后進(jìn)一步完善。本文所采用的評(píng)價(jià)方法高效、簡(jiǎn)潔，減少了人為因素在預(yù)測(cè)中的影響作用，且可視化程度較高、便于數(shù)據(jù)管理并不斷完善，可以為礦山巖溶塌陷預(yù)警系統(tǒng)的建立奠定基礎(chǔ)。

[1]武強(qiáng)，薛東，連會(huì)青.礦山環(huán)境評(píng)價(jià)方法綜述[J].水文地質(zhì)工程地質(zhì).2005，03:84-88.

[2]黃敬軍，李向前，陸華.江蘇省露采礦山地質(zhì)環(huán)境綜合評(píng)價(jià)研究[J].水文地質(zhì)工程地質(zhì).2005，05:97-100.

[3]胡成，陳植華，丁國平，等.GIS技術(shù)在巖溶塌陷預(yù)測(cè)中的應(yīng)用[J].桂林工學(xué)院學(xué)報(bào).2000，20(02):117-119.

[4]李文臣，劉發(fā)榮，等.江西省永豐縣石馬鎮(zhèn)上升礦區(qū)飾面大理石礦資源儲(chǔ)量核實(shí)報(bào)告[R].中材地質(zhì)工程勘查研究院.2012.

[5]賓紅衛(wèi).層次因子綜合預(yù)測(cè)法在巖溶地面塌陷地質(zhì)災(zāi)害預(yù)測(cè)分析評(píng)估中的應(yīng)用[J].西部探礦工程.2010，04:111-117.