蔣莊煤礦太原組煤層開采底板突水危險性分區分級研究

尹會永，Liliana Lefticariu，魏久傳，謝道雷，肖樂樂，張偉杰，佟文亮

（1.中國礦業大學 （北京）地球科學與測繪工程學院，北京100083；2.山東科技大學地球科學與工程學院 山東省沉積成礦作用與沉積礦產重點實驗室，山東 青島266590；3.山東科技大學 礦山災害預防控制省部共建國家重點實驗室培育基地，山東 青島266590；4.南伊利諾伊大學地質系，伊利諾伊卡本代爾62901；5.棗莊礦業 （集團）有限責任公司蔣莊煤礦，山東 棗莊 277519）

蔣莊煤礦太原組煤層開采底板突水危險性分區分級研究

尹會永1，2，3，Liliana Lefticariu4，魏久傳2，3，謝道雷2，肖樂樂2，張偉杰2，佟文亮5

（1.中國礦業大學 （北京）地球科學與測繪工程學院，北京100083；2.山東科技大學地球科學與工程學院 山東省沉積成礦作用與沉積礦產重點實驗室，山東 青島266590；3.山東科技大學 礦山災害預防控制省部共建國家重點實驗室培育基地，山東 青島266590；4.南伊利諾伊大學地質系，伊利諾伊卡本代爾62901；5.棗莊礦業 （集團）有限責任公司蔣莊煤礦，山東 棗莊 277519）

華北地區越來越多的煤礦將要開采太原組下部煤層，基底奧灰突水危險性隨之增加。為合理評價滕州礦區蔣莊煤礦162、163采區太原組16煤底板奧陶突水危險性，在分析采區地質及水文地質特征的基礎上，利用基于GIS的AHP型“脆弱性指數法”，選取奧灰（富水性與水壓力）、斷裂構造（密度、斷層交點和端點的分布、斷層規模指數）和有效隔水層厚度因素作為評價指標，建立預測模型。應用本模型對研究區奧灰突水危險性進行了分區分級，分為安全、較安全、較危險和危險四個等級，對每個等級的分布進行了分區。經驗證，預測效果符合實際，并與“突水系數法”預測結果相比，更詳細、具體。分區分級結果有助于指導礦井16煤開采防治水工作。

太原組；底板突水；脆弱性指數；分區分級

華北含煤區作為我國主要的產煤區之一，多年來以開采山西組和太原組上部煤層為主，但隨著上部煤層的日漸殆盡，越來越多的煤礦要面臨開采太原組下部煤層的局面。然而，太原組下部煤層距離基底奧陶系灰巖（奧灰）巖溶含水層較近，易發生底板突水災害。如何較為準確的確定底板突水危險性是礦井防治水工作需要解決的主要任務。我國礦井水文地質工作者經過近六十多年的實踐、研究，已形成了許多實用、先進的評價理論和方法，如“突水系數法”［1］、下三帶”理論［2］、“關鍵層”（KS）理論［3］、薄板結構理論［4］、“零位破壞”與“原位張裂”理論［5］、基于 GIS的多源信息復合［6］、突水概率指數法［7］、神 經 網 絡［8］、優 勢 面［9］、信 息 分 析－熵 的 方法［10］、“脆弱性指數法”［11］、“突水系數－單位涌水量法”［12］、模糊聚類法［13］、無量綱信息融合法［14］、基于巖體極限平衡理論［15］、基于模糊物元理論［16］等理論和方法，對預測底板突水起到了良好的指導作用。近年來，“脆弱性指數”系列方法應用越來越廣［17－19］，本文即利用脆弱性指數系列方法中的基于GIS（地理信息系統）的AHP（層次分析）型評價方法對山東省滕州礦區蔣莊煤礦太原組下部煤層開采受底板奧灰突水危險性進行分區分級，并對比突水系數評價結果，以指導礦井水害防治工作。

1 礦井概況

蔣莊煤礦位于濟寧市微山縣歡城鎮和滕州市西崗鎮交界處，北距滕州市城區25km，南距微山縣城16km。井田基本呈長方型，東西寬約3.5km，南北長約10km，面積36.463km2。蔣莊煤礦于1989年6月24日正式投產，主采3上、3下、12下、16煤，核定生產能力為275萬t/a，生產水平為－320m和－430m。

隨著蔣莊煤礦上組煤煤炭資源的日漸殆盡，礦井今后將逐漸進入開采下組煤為主。其中太原組16煤將是今后主要開采對象。16煤的162采區、163采區是礦井下組煤先期開拓開采區域，162采區、163采區自西向東埋藏深度逐漸增大，底板突水威脅逐漸增加，尤其是底板奧灰承壓含水層水壓也越來越大，奧灰巖溶水的突水危險性也隨之增大，因此研究16煤底板突水危險性對于礦井今后開展防治水工作具有重要的指導意義。162采區、163采區即為本文研究對象，162采區、163采區位于井田的中西部，面積約為6.67km2。

2 162、163采區地質與水文地質概況

2.1 162采區、163采區地質概況

該采區地層屬典型的華北型含煤地層，屬于隱蔽型煤田。地層發育有寒武系、奧陶系、石炭系、二疊系、侏羅系、第四系。井田位于滕南煤田的中部，162采區、163采區總體為一傾向南東的單斜構造，在162采區發育一向斜，地層傾角在3～10°。采區目前共揭露及控制斷層19條，斷層密度2.8條/km2，其中：落差（H）≥100m有3條，50≤H＜100m有3條，20≤H＜50m有8條，10≤H＜20m有3條，5≤H＜10m有2條。

采區內16煤底板標高為－400～－560m，埋藏深度整體呈由西向東逐漸加深。16煤厚度0.7～1.91m，平均0.89m，呈現由北向南逐漸變厚的趨勢。采區內及其周邊共有49個鉆孔，其中48個鉆孔穿過16煤，鉆探密度為6.15孔/km2。采區16煤資源儲量約750萬t。163采區西部已回采16301工作面，16303工作面、16306工作面已放出。

2.2 162采區、163采區水文地質概況

2.2.1 含水層

礦井直接充水含水層主要有：山西組3煤頂砂（包括3上、3下煤層頂部砂巖簡稱3煤頂砂）、太原組第三層石灰巖（簡稱三灰）、太原組第十下層灰巖（簡稱十下灰）；間接充水含水層主要有：第四系砂礫層、上侏羅統礫巖、石盒子組頂底部砂巖、本溪組第十四層石灰巖（簡稱十四灰）、奧陶系石灰巖（簡稱奧灰）。本文僅對16煤底板具有突水威脅的奧灰含水層進行簡要敘述。

奧灰為青灰、灰白色厚層狀石灰巖。淺部裂隙溶洞發育，富水性弱至強，為煤系含水層重要補給水源。鉆孔揭露奧灰上部鋁質較多，下部少量方解石、少量泥質充填，含硅質結合，堅硬，完整，縫合線發育，方解石充填裂隙，堅硬完整。單位涌水量0.000125～0.612L/s.m，水質為硫酸鹽、氯化物－鈣鈉型，礦化度3.828～4.9g/L，2014年6月31日水位標高－169.68m（水12－2孔）～－156.95m（2006－1孔）。

2.2.2 隔水層

主要隔水層為：第四系隔水層、太原組煤層間的黏土巖和本溪組黏土巖。其中本溪組黏土巖主要有十二灰至十四灰之間的黏土巖，間夾細砂巖；16煤與奧灰間距平均57.52m。十四灰與奧灰之間為雜色黏土巖、鐵質泥巖、黏土巖、G層黏土巖。對于阻隔奧灰水上升、減少奧灰水對礦井生產的威脅起到重要作用。

3 煤層底板突水威脅性分區分級研究

為了確定研究區煤層開采受底板奧灰突水的威脅程度，并進行威脅程度分區分級，本文利用基于GIS的AHP型“脆弱性指數法”的原理建立各個影響因素專題圖，并確定各個影響因素的權重，最終綜合各影響因素并統一在預測模型中，建立研究區煤層開采受底板奧灰突水威脅程度的分區分級圖。

3.1 影響因素分析及其圖件

在分析蔣莊煤礦已有的生產資料和參考其他文獻的基礎上［20－21］，本文選取影響突水的“矛”（奧灰水壓及其富水性，斷裂構造（密度、交點和端點個數、規模））和“盾”（有效隔水層厚度）等主要影響因素并繪制各個圖件（圖1）。

圖1 各影響因素圖

1）煤層底板奧灰含水層。奧灰是影響162采區、163采區16煤開采最具威脅性的含水層，分析其水壓與富水性特征。①奧灰含水層水壓是發生底板突水的動力因素。利用蔣莊煤礦奧灰水文觀測孔的水位和標高數據以及區內鉆孔的相關數據推算得到奧灰水頭壓力，據此生成研究區奧灰水壓圖（圖1（a），單位為 MPa）。可以看出，奧灰水壓是隨著埋深的增加而增加，突水的概率也隨之增加。奧灰水壓最大在162采區中部，呈現由西向東逐漸變大的趨勢。②奧灰含水層的富水性是發生底板突水物質基礎。井田內奧灰含水層的富水性主要受埋深、斷裂構造、巖溶發育程度等因素控制，具有極大的不均一性。根據收集到的抽水鉆孔單位涌水量大小、鉆孔漏失量、埋深、斷裂構造等因素，將奧灰富水性劃分為弱、中等、富三個級別。據此數據插值生成了富水性等值線圖，進而制作奧灰含水層的富水性圖（圖1（b），圖中數字為相對大小）。總體呈現，采區邊緣富水性較弱，中部富水性較強。

2）斷裂構造。斷裂構造的存在增加了突水概率，斷裂構造的分布及其發育情況，規模大小等對底板突水影響較大。①斷層存在的區域，突水概率增加，因此，斷層分布是影響煤層底板突水的重要因素。在采區16煤底板等高線圖上，以200m×200m的網度，按照文獻［11］的方法繪制斷層密度等值線圖（圖1（c），單位為條數），整體西部復雜，東部較為簡單。②斷層規模主要從斷層落差和延展長度考慮，相同條件下，斷層落差越大、延展長度越長，則斷層效應越大，導水性增強。以200m×200m的網度繪制斷層規模指數等值線圖（圖1（d），單位為落差與長度乘積），整體西部復雜，東部較為簡單。③在斷層相交與斷層的端點處，巖體裂隙發育，導水的可能性增強。在采區16煤底板等高線圖上，同樣以200m×200m的網度繪制斷層點和端點分布圖（圖1（e），單位為個數），整體西部復雜，東部較為簡單。

3）隔水層厚度。根據“下三帶”理論，真正起到阻水作用的是有效隔水層，其數值等于16煤至奧灰間距減去底板破壞帶數值與導升帶高度。根據鉆孔資料，先統計間距值，然后再減去以“三下壓煤規程”［22］中的公式計算出來的底板破壞帶深度（工作面寬度以現有工作面值為參考，為160m）。資料表明，采區原始導升高度和礦壓導升高度為零。16煤至奧灰有效隔水層厚度見圖1（f）（單位為m），可知，16煤底板有效隔水層厚度為34.62～37.32m，平均35.91m，總體呈自西向東逐漸變薄的趨勢。

3.2 基于AHP方法的影響權重確定

通過前述影響因素分析可知，為了確定底板奧灰突水危險性，需對“矛”和“盾”的指標進行權重衡量，如圖2所示，建立三個級別來逐一確定權重大小。

根據對影響蔣莊煤礦162采區、163采區煤層底板突水因素的分析，運用“征集專家評分”的方法［11］，來構建162采區、163采區16煤底板奧灰突水評價的判斷矩陣，并最終獲得各個影響因素的權重大小（表1～4）。

通過表1～4可知，各組矩陣計算出λ，CI與CR，存在的CR值都小于0.1，判斷矩陣通過一致性檢驗。所以，得到的各因素權重值如表4中w列所示。

圖2 底板奧灰突水危險性層次結構示意圖

表1 判斷矩陣A～Bi（i＝1～3）

表2 判斷矩陣B1～Ci（i＝1～2）

表3 判斷矩陣B2～Ci（i＝3～5）

表4 各因素權重確定表

3.3 數據歸一化

前述各因素的單位不一致，不能統一在模型中，因此需要進行歸一化，選用式（1）［23］進行。

式中：Si為計算所得數據；設定a＝0.1，b＝0.9；xi為歸一化前的各個影響因素數值。

各因素數據經過歸一化處理后，即可運用SURFER建立格網文件，選用50m×50m的格網，共劃分4941個網格，建立上述六個因素的歸一化格網圖，且要保證范圍、格網位置、大小均要一致。其中前5個因素對底板突水都屬于正相關，但有效隔水層厚度正好相反，因此，需對其進行處理，S′i＝1－Si，S′i為處理后所得數據。

3.4 預測模型建立

建立煤層底板突水脆弱性模型［11］，即建立能反映蔣莊煤礦162采區、163采區16煤層底板突水危險性情況的模型，此模型包括前述的各個影響因素以及各因素的影響權重。

利用脆弱性指數（VI）模型來對煤層底板突水危險性進行評價。模型如式（2）所示。

3.5 分區分級結果及其分析

使用該模型對16煤各個格網點的脆弱性指數進行統計、分析［11］，確定出分級的界限值為0.32、0.57、0.68。根據界限值對162采區、163采區進行區域劃分，從而得16煤開采受底板奧灰突水危險分區分級圖（圖3（a））。VI＞0.68突水脆弱區（危險區）；0.57＜ VI≤0.68突水較脆弱區（較危險區）；0.32＜ VI≤0.57突水過渡區（較安全區）；VI≤0.32安全區。

由圖3可以看出，安全區范圍較小位于采區西部和中南部；過渡區（較安全區）范圍較大，位于采區西南部和東南部；較脆弱區（較危險區）主要位于162采區，脆弱區（危險區）位于162采區深部。評價結果與利用傳統的“突水系數法”計算得到的突水系數評價圖（圖3（b））相似，但脆弱性分區更加詳細。但突水系數評價圖中的16301工作面位于突水系數大于0.06MPa/m的區域，精度不如本文的評價結果。

3.6 模型的檢驗

由于研究區目前只開采了16301工作面，已實現安全開采。16301工作面開采時只有十下灰水涌水，水量較小，初期揭露時為5m3/h左右，后減小至1～2.5m3/h。另外，進行的底板含水層富水性探測盡管表明十四灰和奧灰有一定的富水性，但實際開采并無突水。16301工作面位于圖3（a）中的安全區，也說明評價結果與實際相吻合，所以，評價結果能反映實際條件，可以指導礦井今后開采162采區、163采區的防治水工作。

圖3 162采區、163采區16煤底板突水危險性分區分級圖（a）與突水系數圖（b）

4 結 論

1）蔣莊煤礦162采區、163采區16煤開采主要水害威脅來自煤系基底奧陶系巖溶含水層，16煤至奧灰間距平均57.52m，巖性主要為灰巖、雜色黏土巖、鐵質泥巖、黏土巖、G層黏土巖等，對于阻隔奧灰水上升、減少奧灰水對礦井安全生產的威脅起到重要作用。

2）利用基于GIS的AHP的“脆弱性指數”方法，選取奧灰含水層特征（富水性與水壓力）、斷裂構造（密度、交點和端點的分布、規模指數）和有效隔水層厚度等三個方面、六個主要影響因素，對底板突水的影響權重由大至小分別為：斷裂構造、奧灰含水層特征和有效隔水層厚度，并建立了預測模型。

3）利用本預測模型對162采區、163采區底板奧灰突水危險性進行了分區分級，共分為安全、較安全、較危險和危險四個等級，對每個等級的分布進行了分區。經驗證，預測效果與實際較為吻合，并經與“突水系數法”預測結果進行對比，更詳細、具體，可用來指導礦井防治水工作。

［1］段水云.煤層底板突水系數計算公式的探討［J］.水文地質工程地質，2003（1）：96－99.

［2］李加祥.煤層底板“下三帶”理論在底板突水研究中的應用［J］.河北煤炭，1990，9（4）：12－15.

［3］黎良杰.采場底板突水機理的研究［D］.徐州：中國礦業大學，1995.

［4］張金才，張玉卓，劉天泉.巖體滲流與煤層底板突水［M］.北京：地質出版社，1997.

［5］王作宇.底板零位破壞帶最大深度的分析計算［J］.煤炭科學技術，1992（2）：2－6.

［6］鄭世書，孫亞軍.煤礦采掘過程中底板突水危險性的預測預報［C］//北京：中國煤炭學會學術年會論文集，1992.

［7］施龍青，韓進，宋揚.用突水概率指數法預測采場底板突水［J］.中國礦業大學學報，1999，28（5）：442－446.

［8］靳德武，陳健鵬，王延福，等.煤層底板突水預報人工神經網絡系統的研究［J］.西安科技學院學報，2000，20（3）：214－216.

［9］倪宏革，羅國煜.煤礦水害的優勢面機理研究［J］.煤炭學報，2000，25（5）：518－521.

［10］靳德武.華北型煤田煤層底板突水預測信息分析理論－方法及應用［J］.中國巖溶，2003，22（1）：41－46.

［11］武強，王金華，劉東海，等.煤層底板突水評價的新型實用方法Ⅳ：基于GIS的AHP型脆弱性指數法應用［J］.煤炭學報，2009，34（2）：232－235.

［12］喬偉，李文平，趙成喜.煤礦底板突水評價突水系數－單位涌水量法［J］.巖石力學與工程學報，2009，28（12）：2466－2474.

［13］李加祥.用模糊數學預測煤層底板的突水［J］.山東礦業學院學報：自然科學版，1990，9（1）：11－18.

［14］李忠建，魏久傳，尹會永，等.底板突水影響因素評價新方法—無量綱信息融合法［J］.中國礦業，2010，19（1）：95－103.

［15］張偉杰，李術才，魏久傳，等.基于巖體極限平衡理論的煤層底板突水危險性預測［J］.山東大學學報：工學版，2013，43（1）：86－91.

［16］劉偉韜，劉云娟，申建軍.基于模糊物元理論的深部開采底板突水安全性評價［J］.山東科技大學學報：自然科學版，2014，33（3）：25－31.

［17］武強，楊柳，朱斌，等.“脆弱性指數法”在趙各莊礦底板突水評價中的應用［J］.中國煤炭地質，2009，21（6）：40－44.

［18］武強，劉守強，賈國凱.脆弱性指數法在煤層底板突水評價中的應用［J］.中國煤炭，2010，36（6）：16－22.

［19］羅成.GIS的AHP型脆弱性指數法在李雅莊礦底板突水評價中的應用［J］.煤田地質與勘探，2012，40（5）：47－50.

［20］尹會永，魏久傳，郭建斌，等.翟鎮煤礦下組煤開采底板突水危險性評價［J］.中國礦業，2010，19（1）：95－103.

［21］李云龍.基于多源信息融合的姚橋礦水害防治方法與應用研究［D］.北京：中國礦業大學（北京），2009.

［22］國家煤炭工業局.建筑物水體鐵路及主要井巷煤柱留設與壓煤開采規程［M］.北京：煤炭工業出版社，2000.

［23］肖樂樂，魏久傳，尹會永，等.南屯煤礦下組煤西區煤層底板突水危險性評價［J］.礦業安全與環保，2012，39（5）：85－88.

Water inrush risk zoning and classification during coal seam mining in the Jiangzhuang coal mine，Shandong province

YIN Hui－yong1，2，3，Liliana Lefticariu4，WEI Jiu－chuan2，XIE Dao－lei2，XIAO Le－le2，ZHANG Wei－jie2，TONG Wen－liang5
（1.College of Geoscience and Surveying Engineering，China University of Mining and Technology（Beijing），Beijing 100083，China；2.Shandong Provincial Key Laboratory of Depositional Mineralization ＆Sedimentary Minerals，College of Earth Sciences ＆Engineering，Shandong University of Science ＆Technology，Qingdao 266590，China；3.State Key Laboratory of Mining Disaster Prevention and Control Co－founded by Shandong Province and the Ministry of Science and Technology，Shandong University of Science ＆Technology，Qingdao 266590，China；4.Department of Geology，Southern Illinois University，Carbondale，Illinois 62901，USA；5.Jiangzhuang Coal Mine，Zaozhuang Mining（Group）Co.，Ltd，Zaozhuang 277519，China）

An increasing number of coal mines will focus their operations on the coal seams present in the early Carboniferous Taiyuan Group in North China.It is expected the risk of water inrush from the Ordovician confined limestone aquifer through the floor of overlying excavations to be significant during mining.The current study is using the No.162，163 district in Jiangzhuang coal mine，which is located in the Tengzhou coal field in the Shandong province，as an example of how to calculate and zone this risk based on the analysis of local geological and hydrogeological features such as hydraulic pressure and water yield capacity of the limestone aquifer，fault density，fault crossing and endpoint density，scale index of faults，and effective aquiclude thickness.The vulnerability index was determined using Geographical Information Systems and Analytic Hierarchy Process techniques in order to build the risk prediction model.Four areas were identified with the prediction model based on the risk of water inflow to the floor of the coal seam：safe，transitional，less dangerous，and dangerous.The model validation was carried out by using actual data from the No.16301 workface and showed that it predicted the water inflow risk better than the traditional water inrush coefficient method.The model was considered to be accurate and therefore reliable and helped with the policy to prevent and control groundwater hazards during the mining of coal seam No.16.

Taiyuan group；floor water inrush；vulnerability index；zoning and classification

TD823.84

A

1004－4051（2015）10－0149－06

2014－11－06

國家重點基礎研究發展計劃（973）項目資助（編號：2012CB723104）；國家自然科學基金資助（編號：41372290；41402250）；山東省自然科學基金資助（編號：ZR2013EEQ019）；山東省礦山災害預防控制國家重點實驗室培育基地開放基金資助（編號：MDPC2012KF13）

尹會永（1979－），男，博士后，副教授，主要從事礦業地質及水文地質方面的教學與研究工作。E－mail：yhy1919＠163.com。