楊華澤，郭 琳，王騰輝

(1.安徽理工大學(xué)地球與環(huán)境學(xué)院，安徽 淮南 232001;2.江西地礦局贛西地質(zhì)調(diào)查大隊(duì)，江西 南昌 330030)

煤炭是我國(guó)的主體能源，即使隨著新能源的應(yīng)用和發(fā)展，煤炭依舊占據(jù)一次能源生產(chǎn)和消費(fèi)結(jié)構(gòu)中的70%以上。由于淺部煤炭資源開采的日益枯竭，深部開采已成為我國(guó)煤炭生產(chǎn)的必然趨勢(shì)。隨著開采深度和強(qiáng)度的增加，地質(zhì)條件越來越復(fù)雜，地壓、水壓、地溫亦將不斷增加，導(dǎo)致煤層開采過程中受底板高承壓裂隙巖溶水的危害頻率愈加增多?；幢泵禾锾幱谌A北聚煤期南部，是典型的華北型煤田。在其主要含煤地層石炭—二疊系中的下組煤開采過程中，煤層底板廣泛遭到太原組灰?guī)r水和奧陶系灰?guī)r水的威脅，發(fā)生突水事故。譬如，楊莊礦的Ⅱ617工作面發(fā)生了突水，造成礦區(qū)內(nèi)的的4個(gè)生產(chǎn)采區(qū)被淹沒;桃園煤礦二采區(qū)的四個(gè)工作面和南三采區(qū)的1035工作面等發(fā)生了重大突水事故，造成巨大的經(jīng)濟(jì)損失和人員傷亡?？梢?，對(duì)礦井水害的預(yù)測(cè)和防治已經(jīng)是煤礦開采過程中的重中之重，其中，灰?guī)r巖溶類突水是礦井水害防治的重點(diǎn)[1]。

1 采區(qū)概況

楊柳煤礦位于安徽省淮北市孫疃鎮(zhèn)楊柳村，距宿州市23 km，面積約為60.4 km2。楊柳煤礦太原群灰?guī)r含水層為 10煤底板直接充水含水層，該含水層厚度大，水量豐富，水壓高，若與礦井形成直接溝通，后果將會(huì)很嚴(yán)重。

根據(jù)鉆孔數(shù)據(jù)統(tǒng)計(jì)結(jié)果可知，井田內(nèi)共有44個(gè)見太灰的孔，其中，僅揭露一灰的鉆孔達(dá)到31個(gè);揭露到四灰的孔有7個(gè);打穿太原組并揭露奧灰的孔有4個(gè)。其主要是由深灰～灰色的石灰?guī)r、泥巖、粉砂巖和薄煤層組成。

2 底板突水影響因素分析

影響底板突水的因素可歸納為兩點(diǎn):承壓水壓力和有效隔水層厚度及其阻水性能[2]。10煤層是楊柳煤礦的主采煤層，其承壓水上開采的成功與否直接關(guān)系的整個(gè)煤礦的開采安全，意義重大。本次采用灰色模糊綜合評(píng)價(jià)法對(duì)煤層底板的突水危險(xiǎn)性進(jìn)行分區(qū)和評(píng)價(jià)，為10煤層的安全開采提供有利的地質(zhì)依據(jù)。

2.1 巖溶裂隙

灰?guī)r巖溶裂隙的評(píng)價(jià)方法有很多，如直接統(tǒng)計(jì)鉆孔巖心中太灰段裂隙發(fā)育的數(shù)量;通過鉆孔在鉆進(jìn)過程中的漿液消耗量;通過鉆孔的太灰段的巖心采取率。其中，巖心的采取率是指鉆孔中的某一孔段的巖心的長(zhǎng)度與此段的實(shí)際進(jìn)尺的比值，它能夠直接的反映出巖層的完整性。一般來說，巖心采取率低，巖層較破碎，裂隙的發(fā)育程度高，富水性好，底板突水的可能性高;反之，若巖心的采取率高，則巖層的裂隙發(fā)育程度低且富水性弱，底板突水可能性比較低。由于楊柳煤礦中見灰的47個(gè)鉆孔中有35個(gè)孔僅打到一灰，因此選取一灰的巖心采取率作為巖溶裂隙發(fā)育程度的指標(biāo)來評(píng)價(jià)底板的突水性，利用插值法并根據(jù)一灰的巖心采取率的等值線圖來得出每個(gè)塊段的巖心采取率的值。

2.2 斷裂構(gòu)造

斷裂構(gòu)造是決定煤層底板突水地點(diǎn)的主要因素之一，因?yàn)樵跀嗔褬?gòu)造較發(fā)育、斷層的密度較大的地方，其富水性較強(qiáng)，突水的可能性也比較大。斷裂分維值恰好能綜合的反映出斷裂構(gòu)造的發(fā)育程度，它包括斷裂的長(zhǎng)度、密度以及裂隙相互交叉關(guān)系等信息。一般來說，分維值越大，則表示斷裂的長(zhǎng)度越長(zhǎng)、密度越大、分支越多，其對(duì)巖體的切割越嚴(yán)重，富水性越強(qiáng)，易發(fā)生底板突水;分維值越小，則表示斷裂的長(zhǎng)度比較小、分布比較稀疏、分支也比較少，對(duì)巖體的切割也比較輕微，富水性越弱，不易發(fā)生底板突水。因此，選取斷裂分維值作為斷裂構(gòu)造發(fā)育程度的指標(biāo)來評(píng)價(jià)底板的突水性。由鉆孔資料的統(tǒng)計(jì)知，本礦區(qū)10煤底板的斷裂構(gòu)造分維值為0.831 7～1.549 1，利用插值法并根據(jù)斷裂構(gòu)造分維值的等值線圖來得出每個(gè)塊段的斷裂構(gòu)造分維值。

2.3 砂巖率

底板砂巖率是指煤層底板巖層段中砂巖層的厚度占總體巖層厚度的比值，它是決定巖層富水性的前提，是影響地下水的賦存和巖層隔水性能的主要因素。一般來說，砂巖率高，巖層的儲(chǔ)水量大，富水性好，隔水性能弱，底板的突水可能性大;砂巖率低，巖層的儲(chǔ)水量小，富水性弱 ，隔水性能強(qiáng)，底板的突水可能性小。因此，選取煤層底板隔水層的砂巖率作為指標(biāo)來評(píng)價(jià)底板的突水性。由鉆孔資料的統(tǒng)計(jì)知，本礦區(qū)10煤底板至一灰頂隔水層段主要由粉砂巖和泥巖構(gòu)成，砂巖率的值為0.556 5～1，利用插值法并根據(jù)10煤層底板至一灰頂巖層段的砂巖率等值線圖來得出每個(gè)塊段的砂巖率，如圖 1。

圖1 楊柳煤礦10煤至一灰頂砂巖率圖

2.4 底板隔水層厚度

隔水層厚度是指開采煤層的底板與下伏的含水層的頂面之間的隔水巖層的厚度，其值的大小會(huì)直接影響到煤層底板是否突水。一般來說，隔水層的厚度越大，底板的阻隔水能力越強(qiáng)，突水的可能性越小;反之，隔水層的厚度越小，底板的阻隔水能力越弱，突水可能性就越大。

據(jù)該礦區(qū)鉆孔資料得知，在正常情況下，10煤底板至一灰頂?shù)拈g距為 33.90～69.68 m，平均值為 59.16 m，但由于煤層開采會(huì)造成底板巖層發(fā)生破壞，產(chǎn)生裂隙，使部分隔水層失去隔水能力甚至變成導(dǎo)水層。由底板破壞深度經(jīng)驗(yàn)公式計(jì)算法和FLAC3D數(shù)值模擬結(jié)果的值，綜合得出底板破壞深度為13.15 m。因此，楊柳煤礦的10煤底板有效隔水層厚度為 20.75～56.53 m，平均厚度約為 46.01 m。由于煤層底板的有效隔水層厚度更能準(zhǔn)確的表示出煤層底板隔水層的隔水性能，因此選取10煤底板的有效隔水層厚度指標(biāo)來評(píng)價(jià)底板的突水性，利用插值法并根據(jù)10煤底板有效隔水層厚度的等值線圖來得出每個(gè)塊段的有效隔水層厚度，如圖2。

3 底板突水危險(xiǎn)性預(yù)測(cè)

影響底板突水的因素可歸納為兩點(diǎn):承壓水壓力和有效隔水層厚度及其阻水性能[2]。10煤層是楊柳煤礦的主采煤層，其承壓水上開采的成功與否直接關(guān)系的整個(gè)煤礦的開采安全，意義重大。本次采用灰色模糊綜合評(píng)價(jià)法對(duì)煤層底板的突水危險(xiǎn)性進(jìn)行分區(qū)和評(píng)價(jià)，為10煤層的安全開采提供有利的地質(zhì)依據(jù)。

3.1 灰色模糊綜合評(píng)價(jià)法

模糊數(shù)學(xué)是將模糊信息用數(shù)學(xué)方法來描述的工具，揭露出模糊信息的本質(zhì)和規(guī)律?；疑碚撌怯脭?shù)學(xué)方法處理復(fù)雜系統(tǒng)，了解到系統(tǒng)內(nèi)部的變化趨勢(shì)以及相互關(guān)系等。若只用模糊數(shù)學(xué)法易使信息丟失;若只用灰色理論則不能體現(xiàn)評(píng)價(jià)規(guī)則的模糊性，兩種結(jié)果均會(huì)造成結(jié)果的偏差。因此，將這兩種方法結(jié)合運(yùn)用進(jìn)行綜合評(píng)判，即灰色模糊綜合評(píng)判法。本次選取含水層厚度、巖溶裂隙發(fā)育、斷裂構(gòu)造、10煤至一灰的砂巖率以及底板隔水層的有效厚度作為評(píng)價(jià)指標(biāo)，建立數(shù)學(xué)模型，對(duì)底板的突水危險(xiǎn)性進(jìn)行評(píng)價(jià)。

圖2 楊柳煤礦10煤底板至一灰頂厚度等值線圖

3.2 底板突水危險(xiǎn)性評(píng)價(jià)

在建模前需將各影響因素進(jìn)行量化轉(zhuǎn)變?yōu)槎恳蛩?，以便建立模型，量化值如?所示。在模型中，選取斷裂分維值作為主序列(X0≡ X0(1)，X0(2)，…，X0(n

{}))，含水層厚度(X1)、巖心采取率(X2)、砂巖率(X3)、底板隔水層厚度(X4)作為子序列。首先對(duì)量化后的數(shù)據(jù)進(jìn)行歸一化處理，再根據(jù)子序列(Xi)與主系列(X0)在k點(diǎn)時(shí)的關(guān)聯(lián)度計(jì)算公式，即:

其中，Δi(k)=|X0(k)－Xi(k)|，ρ為分辨率系數(shù)，其值可以決定最大差對(duì)因素?cái)?shù)據(jù)轉(zhuǎn)化的影響程度，本次選取ρ=0.5。則子序列對(duì)主系列的關(guān)聯(lián)度系數(shù) γi為:

其中，n為樣本的個(gè)數(shù)，在本次模擬中劃分為56個(gè)塊段，因此n=56。在計(jì)算出關(guān)聯(lián)度系數(shù)后，對(duì)結(jié)果進(jìn)行單位化處理，得出因素權(quán)重[3](αi)。再選取含水層厚度、巖心采取率、斷裂分維值、砂巖率和底板隔水層厚度五個(gè)因素為底板突水危險(xiǎn)性的評(píng)價(jià)指標(biāo)，則因素集可表示為 U={u1，u2，u3，u4，u5}。在模型的評(píng)價(jià)過程中，由于巖心采取率以及底板隔水層厚度兩個(gè)因素與其余的因素為負(fù)相關(guān)，因此為了便于計(jì)算，取其相對(duì)于的正值。對(duì)評(píng)價(jià)指標(biāo)的數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，確定出突水危險(xiǎn)性等級(jí)的評(píng)價(jià)標(biāo)準(zhǔn)，如表1所示。設(shè)評(píng)語集為 V={v1，v2，v3，v4}，其中，v1為安全區(qū);v2為相對(duì)安全區(qū);v3為威脅區(qū);v4為危險(xiǎn)區(qū)。

表1 評(píng)價(jià)指標(biāo)等級(jí)劃分標(biāo)準(zhǔn)

因素集U與評(píng)語集V之間的模糊關(guān)系可以用矩陣R來表示:

其中，rij表示第i種因素可以被劃分為第 j級(jí)標(biāo)準(zhǔn)的可能性，即i對(duì)j的隸屬度。運(yùn)用聚類權(quán)重的計(jì)算公式:

其中，n指所選取的評(píng)價(jià)因素的個(gè)數(shù)，在此n=5;Wi是其中的第i個(gè)評(píng)價(jià)因素的聚類權(quán)重值;xi是因素歸一化處理后的數(shù)值;Si1是第i個(gè)因素的評(píng)語集所相應(yīng)的1類和2類的分界值。

白化權(quán)函數(shù)可以對(duì)某對(duì)象隨評(píng)價(jià)指標(biāo)數(shù)據(jù)的大小而變化的過程進(jìn)行定量描述，能直觀反映出此評(píng)價(jià)對(duì)象隸屬與劃分的某個(gè)灰類的程度。白化函數(shù)的確定可以說是在評(píng)價(jià)過程中由定性分析到定量建模的關(guān)鍵。

基于模糊數(shù)學(xué)理論和表2中的評(píng)價(jià)指標(biāo)等級(jí)劃分標(biāo)準(zhǔn)，把底板突水危險(xiǎn)性評(píng)價(jià)劃分為“安全區(qū)”、“相對(duì)安全區(qū)”、“威脅區(qū)”和“危險(xiǎn)區(qū)”四個(gè)灰類等級(jí)。以底板的砂巖率為例，其四類白化權(quán)函數(shù)為:

3.3 灰色模糊計(jì)算

由上述白化權(quán)函數(shù)[4]，把各塊段內(nèi)的參數(shù)的值帶入函數(shù)得出其分別對(duì)于四個(gè)等級(jí)的白化權(quán)函數(shù)值。由得出的白化權(quán)函數(shù)值和之前計(jì)算的雙權(quán)重值，根據(jù)公式:

得出在第i個(gè)塊段中，第 j個(gè)因素屬于級(jí)別 k的聚類系數(shù)。其中，n為選取的評(píng)價(jià)因素的個(gè)數(shù)，在此 n=5;fij[k](xij)指第i個(gè)塊段中，第j個(gè)因素屬于級(jí)別 k的白化權(quán)函數(shù)值，ηij為其相應(yīng)的雙權(quán)重值。

以第一個(gè)塊段的突水危險(xiǎn)性評(píng)價(jià)為例，根據(jù)上述計(jì)算公式分別求出其相對(duì)于四類等級(jí)的聚類系數(shù)值 σ1[1]、σ1[2]、σ1[3]、σ1[4]，如:σ1[1]= η11f11[1](x11)+ η12f12[1](x12)+ …，+η15 f15[1](x15)

同理，得出 σ1[2]、σ1[3]、σ1[4]的值。選取其中最大值所對(duì)應(yīng)的等級(jí)作為該塊段的突水危險(xiǎn)性等級(jí)，結(jié)果如表2所示。

表2 楊柳煤礦10煤層底板突水危險(xiǎn)性等級(jí)評(píng)價(jià)值及評(píng)價(jià)結(jié)果

Ⅱ

根據(jù)灰色模糊計(jì)算中得出的聚類系數(shù)值，得出危險(xiǎn)性分區(qū)的預(yù)測(cè)結(jié)果，如表2所示。再根據(jù)之前礦區(qū)所劃分塊段的中心點(diǎn)坐標(biāo)，結(jié)合預(yù)測(cè)結(jié)果作出礦區(qū)底板突水危險(xiǎn)性預(yù)測(cè)的分區(qū)圖，如圖3所示。

圖3 楊柳礦10煤層底板突水危險(xiǎn)性分區(qū)圖

由評(píng)價(jià)結(jié)果顯示，本礦區(qū)的突水危險(xiǎn)性等級(jí)主要為危險(xiǎn)區(qū)、威脅區(qū)和相對(duì)安全區(qū)，其中相對(duì)安全區(qū)和威脅區(qū)約占到整個(gè)礦區(qū)面積的85%左右。在煤層的開采過程中，需重點(diǎn)對(duì)危險(xiǎn)和威脅區(qū)進(jìn)行重點(diǎn)的勘探和預(yù)防。

4 結(jié)語

(1)通過對(duì)楊柳煤礦 10煤層隔水底板穩(wěn)定性分析，將煤層底板突水危險(xiǎn)性影響因素統(tǒng)一無量綱化后進(jìn)行灰色關(guān)聯(lián)度比較，不僅提高數(shù)據(jù)使用的準(zhǔn)確性，而且確定影響該煤層底板突水的主要因素含水層厚度、巖溶裂隙發(fā)育、斷裂構(gòu)造、10煤至一灰的砂巖率以及底板隔水層的有效厚度。

(2)以楊柳煤礦10煤層為例，選用灰色關(guān)聯(lián)分析法得出的“因素權(quán)重”和模糊聚類分析法得出的“程度權(quán)重”相結(jié)合的雙權(quán)重思想，不僅能夠給出突水主要控制因素，而且增加底板突水危險(xiǎn)性評(píng)價(jià)的客觀性與準(zhǔn)確性。

[1]施龍青，韓進(jìn).底板突水機(jī)理及預(yù)測(cè)預(yù)報(bào)[M].徐州:中國(guó)礦業(yè)大學(xué)出版社.2004.

[2]李忠建，等.運(yùn)用突水系數(shù)法和模糊聚類法綜合評(píng)價(jià)煤層底板突水危險(xiǎn)性[J].礦業(yè)安全與環(huán)保.2010，(37):24 －26.

[3]孫芳芳.淺議灰色關(guān)聯(lián)度分析方法及其應(yīng)用[J].科技信息.2010，(17):880 －882.

[4]董奮義，等.灰色系統(tǒng)教學(xué)中白化權(quán)函數(shù)的構(gòu)造方法分析[J].華北水利水電學(xué)院學(xué) 報(bào).2010，31(3):97－99.

[5]嚴(yán)大思 ，奚硯濤 ，潘玲玲.基于突水系數(shù)的煤礦底板突水危險(xiǎn)性評(píng)價(jià)系統(tǒng)－以山西太原西山煤田古交礦區(qū)為例[J].中國(guó)煤炭地質(zhì).2010，22(1):31 －34.