李延河 劉寶敏

(中國平煤神馬集團平禹一礦，河南省禹州市，450000)

我國華北石炭—二疊紀(jì)巖溶型煤田普遍存在著底板灰?guī)r承壓水突水災(zāi)害問題，隨著開采深度的加深和生產(chǎn)規(guī)模的擴大，作用在煤層底板上的水壓越來越高，突水危險性越來越大。準(zhǔn)確評價和預(yù)測底板承壓水突水危險性是做好煤礦防治水工作的重要前提，《煤礦防治水規(guī)定》提出了安全隔水層厚度法和突水系數(shù)法，公式簡單實用，已廣泛用于巷道和采煤工作面的底板突水危險性評價。武強院士將確定底板突水多種主控因素權(quán)重系數(shù)的信息融合方法與具有強大空間信息分析處理功能的GIS耦合于一體，提出了評價底板突水危險的脆弱性指數(shù)法。這種方法綜合考慮了各種因素對底板承壓水突水的“貢獻(xiàn)”，評價結(jié)果比突水系數(shù)法更符合實際，從而得到廣泛認(rèn)同。但是，GIS軟件在基層煤礦普及程度十分有限，同時脆弱性指數(shù)缺乏全國統(tǒng)一的標(biāo)準(zhǔn)，從而影響了脆弱性指數(shù)法廣泛應(yīng)用。在地質(zhì)測量等領(lǐng)域廣泛使用的Surfer軟件具有多元信息復(fù)合疊加功能，最突出優(yōu)勢是其擁有十幾種數(shù)據(jù)網(wǎng)格化方法，能夠方便快捷地繪制各種等值線圖和分區(qū)圖，而這些正是脆弱性指數(shù)法常選擇GIS平臺的重要理由所在。目前，以Surfer為平臺評價煤礦底板突水脆弱性的實例較少，但已有成功的案例。為利于脆弱性指數(shù)法的推廣應(yīng)用，本文以平禹一礦為研究對象，以Surfer為平臺，聯(lián)合運用脆弱性指數(shù)和突水系數(shù)評價礦井深部煤層底板巖溶承壓水的突水危險性。

1 平禹一礦灰?guī)r水害背景

平禹一礦位于河南省禹州市北部9 km的古城鎮(zhèn)，井田東西長8 km，南北寬0.35～2.4 km，面積13.5 km2。主采二疊系山西組二1煤層和二3煤層，厚度4～7 m，生產(chǎn)能力達(dá)150萬t/a，現(xiàn)開采標(biāo)高為-415 m。平禹一礦處在禹州巖溶水系統(tǒng)最下游的承壓自流排泄區(qū)，煤層底板太灰和寒灰?guī)r溶發(fā)育，富水性和導(dǎo)水性強，承壓水補給充沛。自1972年以來，已經(jīng)發(fā)生49次突水，造成2次淹井和2次淹采區(qū)。尤其是2005年在東大巷發(fā)生了特大型寒灰水突水，最大突水量達(dá)到38056 m3/h。礦井轉(zhuǎn)入到-415 m水平生產(chǎn)后，作用在煤層底板上的水壓超過3 MPa，煤層底板巖溶承壓水威脅著礦井安全生產(chǎn)，為安全高效生產(chǎn)，需要評價底板突水危險性，并進(jìn)行水害分區(qū)和綜合治理。

2 煤層底板突水影響因素分析

通過綜合分析平禹一礦水文地質(zhì)條件、充水條件及歷史上底板巖溶水突水事故，確定了影響平禹一礦底板寒灰?guī)r溶承壓水突水的4個主控因素，包括寒灰水壓、有效隔水層厚度、斷層和灰?guī)r巖溶發(fā)育程度及富水性。

2.1 水壓

在平禹一礦49次突水中，以太灰L7～8灰?guī)r巖溶水為水源的突水13次，以太灰L1-4灰?guī)r巖溶水為水源的突水8次，以寒武系灰?guī)r巖溶水為水源的突水6次，這也就意味著煤層底板灰?guī)r突水占總突水次數(shù)的55%。歷史上的4次突水淹井和淹采區(qū)事故其水源均為寒灰?guī)r溶水。在27次底板灰?guī)r突水中，21次突水的水壓超過3 MPa，最大水壓為4 MPa，5次突水的水壓在1.8～2.6 MPa之間，1次突水水壓為1.1 MPa。由此可見，水壓是誘發(fā)平禹一礦底板灰?guī)r水突水的關(guān)鍵因素，水壓越高，突水危險性越大。當(dāng)前，平禹一礦寒灰水位為-100 m，-415 m水平承壓3.15 MPa，因此帶壓采煤時仍然存在寒灰水突水的潛在危險。

2.2 有效隔水層厚度

通過對井田內(nèi)300多個鉆孔資料的統(tǒng)計，確定了煤層底板隔水層厚度的空間分布。二1煤層底板至太灰上段灰?guī)r間的隔水層平均厚度為11.3 m，最小厚度為4.3 m，最大厚度為26.2 m。二1煤層底板至太灰下段灰?guī)r間的隔水層平均厚度為57.1 m，最小厚度為44.0 m，最大厚度為73.5 m。二1煤層底板至寒灰含水層間的隔水層平均厚度為87.3 m，最小厚度為72.5 m，最大厚度為108.8 m。煤層底板隔水層起到阻隔其下伏承壓水進(jìn)入礦井的作用，其厚度越大，阻水能力就越強。但受采動影響，煤層底板巖層會產(chǎn)生采動破壞裂隙，使得采動破壞深度范圍內(nèi)的隔水層喪失原有隔水性，真正起阻水作用的是除底板破壞深度后剩余的有效隔水層。

2.3 斷層

大量資料表明，我國80%以上的底板突水發(fā)生在斷裂構(gòu)造附近。平禹一礦以斷層為直接突水通道的底板巖溶水突水雖然只有4次，但因斷層溝通了寒武系灰?guī)r和石炭系薄層灰?guī)r的垂直水力聯(lián)系，又是構(gòu)造薄弱帶和富水帶，因此，斷層引起的突水具有突水強度大和危險性高的特征。突水點位置與構(gòu)造關(guān)系十分密切，斷層及影響帶、小斷層密集帶、地層產(chǎn)狀突變處、斷層轉(zhuǎn)折段或分支處均是易發(fā)生突水的部位。

2.4 灰?guī)r巖溶發(fā)育程度及富水性

在寒武系地層的露頭區(qū)，中統(tǒng)張夏組鮞狀灰?guī)r和上統(tǒng)崮山組白云質(zhì)灰?guī)r巖溶發(fā)育，巖溶形態(tài)主要是巖溶裂隙、溶溝、溶槽和小溶洞。在平禹一礦井田范圍內(nèi)，物探、鉆探和采掘工程揭露證實，寒灰和太灰?guī)r溶發(fā)育，富水性強，富水異常區(qū)分布廣且面積大，沿斷層呈帶狀分布。寒灰?guī)r溶發(fā)育和富水性強的另外一種表現(xiàn)是注漿堵水鉆孔的注漿材料消耗量非常大，如12041工作面的注1孔在L1-3灰?guī)r注漿段共消耗水泥6021.6 t，注2孔在L1-3灰?guī)r注漿段共消耗水泥10701 t，注4孔在寒灰注漿段共消耗水泥14632 t。

3 AHP法確定突水主控因素的權(quán)重

合理確定各主控因素的權(quán)重是運用脆弱性指數(shù)法評價底板突水危險性的關(guān)鍵問題，常用方法有人工神經(jīng)網(wǎng)絡(luò)法(ANN)、證據(jù)權(quán)重法、Logistic回歸法和層次分析法(AHP)、熵權(quán)法-模糊層次法等。為使其與突水系數(shù)的評價結(jié)果能相互對比，選擇了影響平禹一礦底板突水的4個關(guān)鍵因素作為構(gòu)建AHP層次模型和確定權(quán)重的主要依據(jù)，所建立的層次模型結(jié)構(gòu)如圖1所示。

圖1 底板突水危險評價層次結(jié)構(gòu)模型

運用專家評分法對影響突水的因素進(jìn)行評分，構(gòu)造判斷矩陣并計算出各主控因素對突水的貢獻(xiàn)(權(quán)重)。水壓、富水性、有效隔水層厚度和斷層影響指數(shù)共四項主控因素所對應(yīng)的權(quán)重值依次是0.3678、0.1134、0.3571和0.1617。據(jù)此建立了評價平禹一礦煤層底板突水脆弱性的數(shù)學(xué)模型：

式中：VI——脆弱性指數(shù)；

fk(x,y)——單因素影響值函數(shù)；

x,y——地理坐標(biāo)。

4 基于Surfer的突水危險性評價

4.1 脆弱性指數(shù)的計算

GIS具有強大的空間信息分析和疊加功能，所以常用GIS軟件實現(xiàn)煤層底板突水脆弱性的評價。Surfer也具有空間信息疊加功能，而且其插值方法豐富，繪制等值線圖和分區(qū)圖的功能遠(yuǎn)比GIS強，而這些都是脆弱性指數(shù)法評價底板突水危險的重要環(huán)節(jié)。為此，本文基于Surfer軟件評價平禹一礦底板突水危險性。實現(xiàn)步驟如下：對評價因子的數(shù)據(jù)進(jìn)行歸一化處理，然后按網(wǎng)格進(jìn)行剖分量化，繪制各主控因素的專題層圖；按照各因素權(quán)重貢獻(xiàn)值對專題圖進(jìn)行復(fù)合疊加，形成煤層底板突水脆弱性指數(shù)分布圖；結(jié)合突水系數(shù)法及脆弱性法對平禹一礦深部采區(qū)的底板突水危險性進(jìn)行分區(qū)。以下簡要說明繪制水壓專題圖的方法：

(1)對平禹一礦二1煤層底板等高線數(shù)據(jù)進(jìn)行柵格化，創(chuàng)建等高線.GRD〗文件，簡稱A。

(2)對平禹一礦井田寒灰水位數(shù)據(jù)進(jìn)行柵格化，創(chuàng)建水位.GRD〗文件，簡稱B。柵格化時要修改水位數(shù)據(jù)的坐標(biāo)范圍和柵格數(shù)量，使柵格化后的坐標(biāo)范圍和柵格數(shù)量與煤層等高線相同，以便于進(jìn)行相關(guān)的運算。

(3)利用math模塊，對上述兩個柵格數(shù)據(jù)文件進(jìn)行計算，令C1=(B-A)/100，其結(jié)果C1表示巖溶水作用在煤層底板上的水壓，得到的水壓等值線圖就作為水壓專題圖。

采用同樣的方法，可繪制等效隔水層厚度的專題圖C3。由于含水層富水性和斷層是屬性數(shù)據(jù)，要先將二者量化后，才能繪制各自的專題圖，分別用C2和C4表示。

(4)同樣利用math模塊，將4個主控因素的柵格化數(shù)據(jù)分步合并運算，即可得到煤層底板突水脆弱性指數(shù)等值線圖，即：VI=0.3678×C1+0.1134×C2+0.3571×C3+0.1617×C4。

4.2 突水系數(shù)的計算

用脆弱性指數(shù)來評價和劃分煤層底板突水危險，由于尚未建立全國性或地區(qū)性的脆弱性指數(shù)限制，其評價結(jié)果是相對的。而突水系數(shù)是水壓和隔水層厚度的比值，有全國統(tǒng)一的突水臨界值，即受構(gòu)造破壞塊段臨界突水系數(shù)為0.06 MPa/m，正常構(gòu)造塊段臨界突水系數(shù)為0.1 MPa/m。在2009年頒布實施的《煤礦防治水規(guī)定》中，將突水系數(shù)法作為評價采煤工作面底板突水危險的首選方法。鑒于此，本文以Surfer軟件為平臺，再用突水系數(shù)法評價平禹一礦深部采區(qū)的底板突水危險性。對前面計算得到的巖溶水作用在煤層底板上的水壓柵格數(shù)據(jù)再進(jìn)行以下計算：

(1)將作用在煤層底板上的水壓C1與隔水層厚度M進(jìn)行相加，求得巖溶水作用在隔水層底板上的水壓P(單位為MPa)，即：

P=C1+M×0.01

(2)

(2)將水壓P和隔水層厚度M再進(jìn)行除法運算，就得到突水系數(shù)T的柵格化數(shù)據(jù)，并繪制等值線圖。

(3)

4.3 底板突水危險性的綜合評價

綜合脆弱性指數(shù)和突水系數(shù)的兩種評價結(jié)果，并參照平煤神馬集團關(guān)于采煤工作面底板灰?guī)r承壓水突水危險性的劃分方法，將平禹一礦井田劃分為安全區(qū)、突水威脅區(qū)和突水危險區(qū)3個區(qū)域，如圖2所示。

圖2 平禹一礦底板突水危險性分區(qū)圖

結(jié)果顯示，平禹一礦二1煤層底板標(biāo)高大于-60 m區(qū)域，處在灰?guī)r承壓水位以上，煤層底板不承受水壓，不受底板灰?guī)r水的影響，屬于安全區(qū),如圖2中綠色區(qū)域。標(biāo)高-60～-400 m且不受斷層影響的區(qū)域?qū)儆谕凰{區(qū)，如圖2中橙色區(qū)域，處在灰?guī)r含水層水位以下，但突水系數(shù)小于0.06 MPa/m，正常地段沒有底板突水危險，但斷層或隔水層薄弱地段存在突水的可能。標(biāo)高-400 m以深和標(biāo)高-60～-400 m的斷層影響區(qū)域?qū)儆谕凰ｋU區(qū),如圖2中紅色區(qū)域，作用在煤層底板的水壓大，或突水系數(shù)超過0.06 MPa/m，或存在斷層垂向?qū)ǖ?，開采過程中存在突水危險。

4.4 底板灰?guī)r水害的綜合治理

平禹一礦為保障深部的安全生產(chǎn)，對灰?guī)r水害進(jìn)行了綜合探查和防治，相繼實施了水文地質(zhì)條件的綜合探查、巖溶水補給通道的帷幕注漿截流、大流量疏放寒灰水、工作面煤層底板地層預(yù)注漿改造加固、排水系統(tǒng)改造等一系列的防治水措施，治理效果顯著，使平禹一礦在極其復(fù)雜的水文地質(zhì)背景下，實現(xiàn)了連續(xù)10年的安全生產(chǎn)，礦井產(chǎn)能也由原來的60萬t/a提高到150萬t/a。

5 結(jié)語

(1)Surfer軟件具有GIS軟件的空間信息疊加功能，由于其普及性更強，使用方便，用于煤礦底板突水脆弱性評價也十分方便。

(2)以Surfer軟件為平臺，建立了評價平禹一礦底板突水脆弱性的層次模型，綜合脆弱性指數(shù)和突水系數(shù)的評價結(jié)果，將平禹一礦劃分為安全區(qū)、突水威脅區(qū)和突水危險區(qū)。

(3)平禹一礦對-200 m大巷以下采區(qū)灰?guī)r水害采取了綜合探查和防治措施，自2008年13091工作面突水淹采區(qū)以來，沒有再次發(fā)生底板突水，原本極其嚴(yán)重的底板水害實現(xiàn)了可防、可控和可治的目標(biāo)。