深部岩溶水文地质特征及疏放性评价
|
摘要
煤层底板灰岩突水是我国华北型煤田主要防治水问题。随着开采深度的增加,煤层底板承受的灰岩水头压力增大,水害威胁程度也随之增大。因此,如何评价及有效防治深部灰岩水害的问题,对于解放深部煤炭资源和安全生产有着十分重要的现实意义。
本文在查阅国内外文献有关煤层底板灰岩水研究现状以及地下水研究方法的基础上,以淮南矿区谢一矿为典型案例进行研究。该矿水文地质条件复杂,浅部-480m水平以上A组煤已基本开采完毕,目前,A]煤层开采受深部底板灰岩水威胁。本次以井田北翼-660m水平A1煤层底板太原组灰岩水为研究对象,系统分析了区内太原组灰岩水文地质特征；利用FLAC3D数值模拟技术,模拟计算A1煤层开采条件下底板破坏深度；建立了井田北翼Al煤层底板太原组灰岩水文地质概念模型和数学模型,运用Visual Modflow模拟技术对其水文地质参数进行识别；运用修正的突水系数计算公式计算Al煤层底板的突水系数,结合识别后的灰岩含水层模型,模拟计算了不同疏放水方案下的降深和突水系数变化规律；最后,对井田不同突水含水层离子成分差异性进行了分析,运用聚类分析和逐步判别方法对突水水源进行了判别,并且提出了灰岩水防治方法。
Limestone water inrush from coal seam floor is the main water prevention problem in Northern China type coal field. With the increasing of mining depth, the limestone aquifer head pressure beared by coal seam floor has been increased and also the threaten degree of water inrush hazard. Therefore, how to evaluate and prevent the water disaster from deep limestone has an important practical significance for how to mine deep coal resource safely.
Based on consulting foreign and domestic literatures related to research status of coal seam floor of limestone water and its research method, Xieyi coal mine was analyzed as a typical case in Huainan. Hydrogeological conditions of the coal mine was complicate, and group A coal resources above the level of -480m almost had been mined. Now, A1 coal of deep level is being threatened by limestone water inrush from coal seam floor. Taiyuan limestone water of A1coal seam floor in level -660m was taken as the research object at the north of XieYi mine, and the hydrogeology characteristics of the region Taiyuan limestone systemically was analyzed. The failure depth of seam floor was simulated and calculated under the condition of mining A1 coal seam by the FLAC3D numerical simulation technology. The hydrogeology conceptual model and mathematical model of Taiyuan limestone of A1 coal seam floor were established, and the hydrogeological parameters was identified with Visual Modflow software. The water inrush coefficient of A1 coal seam floor was calculated by the modified calculation formulas of water inrush coefficient. Combined with the identified aquifer model, the drawdown and variation of water inrush coefficient were calculated under different drainage schemes. The differences of ionic compositions, which came from different water sources of aquifers, were analyzed, and water inrush sources were discriminated by the ways of cluster analysis method and stepwise discriminant method. And some measures were taken about how to prevent limestone water inrushing from the coal seam floor.
本文在查阅国内外文献有关煤层底板灰岩水研究现状以及地下水研究方法的基础上,以淮南矿区谢一矿为典型案例进行研究。该矿水文地质条件复杂,浅部-480m水平以上A组煤已基本开采完毕,目前,A]煤层开采受深部底板灰岩水威胁。本次以井田北翼-660m水平A1煤层底板太原组灰岩水为研究对象,系统分析了区内太原组灰岩水文地质特征；利用FLAC3D数值模拟技术,模拟计算A1煤层开采条件下底板破坏深度；建立了井田北翼Al煤层底板太原组灰岩水文地质概念模型和数学模型,运用Visual Modflow模拟技术对其水文地质参数进行识别；运用修正的突水系数计算公式计算Al煤层底板的突水系数,结合识别后的灰岩含水层模型,模拟计算了不同疏放水方案下的降深和突水系数变化规律；最后,对井田不同突水含水层离子成分差异性进行了分析,运用聚类分析和逐步判别方法对突水水源进行了判别,并且提出了灰岩水防治方法。
Limestone water inrush from coal seam floor is the main water prevention problem in Northern China type coal field. With the increasing of mining depth, the limestone aquifer head pressure beared by coal seam floor has been increased and also the threaten degree of water inrush hazard. Therefore, how to evaluate and prevent the water disaster from deep limestone has an important practical significance for how to mine deep coal resource safely.
Based on consulting foreign and domestic literatures related to research status of coal seam floor of limestone water and its research method, Xieyi coal mine was analyzed as a typical case in Huainan. Hydrogeological conditions of the coal mine was complicate, and group A coal resources above the level of -480m almost had been mined. Now, A1 coal of deep level is being threatened by limestone water inrush from coal seam floor. Taiyuan limestone water of A1coal seam floor in level -660m was taken as the research object at the north of XieYi mine, and the hydrogeology characteristics of the region Taiyuan limestone systemically was analyzed. The failure depth of seam floor was simulated and calculated under the condition of mining A1 coal seam by the FLAC3D numerical simulation technology. The hydrogeology conceptual model and mathematical model of Taiyuan limestone of A1 coal seam floor were established, and the hydrogeological parameters was identified with Visual Modflow software. The water inrush coefficient of A1 coal seam floor was calculated by the modified calculation formulas of water inrush coefficient. Combined with the identified aquifer model, the drawdown and variation of water inrush coefficient were calculated under different drainage schemes. The differences of ionic compositions, which came from different water sources of aquifers, were analyzed, and water inrush sources were discriminated by the ways of cluster analysis method and stepwise discriminant method. And some measures were taken about how to prevent limestone water inrushing from the coal seam floor.
引文
[1]施龙青,韩进编.底板突水机理及预测预报[M].徐州：中国矿业大学出版社,2004,6.
[2]李金凯等编.矿井岩溶水防治[M].北京：煤炭工业出版社,1990,3.
[3]煤炭部生产协调司.国有重点煤矿井田内受小煤矿开采影响安全生产情况的调查,1995,7.
[4]李定龙,周治安,王桂梁.马家沟灰岩(古)岩溶研究中的若干问题探讨[J].地质科技情报,1997,16(1).
[5]Reibiec M.S.,Hydrofracturing of rock as a method of water,mudmand gas inrush hazards in underground coal mining,4th IMWA,1991,1 (Yugoslavia).
[6]胡春玲.里彦矿灰岩水对下组煤安全开采影响的研究[D].山东科技大,2005.
[7]高延法.底板突水的研究途径与突水优势面[J].山东矿业学院学报,1994,13(增).
[8]王作宇,刘鸿泉.承压水上采煤[M].北京：煤炭工业出版社,1993.
[9]王永红,沈文.中国煤矿水害预防及治理[M].北京：煤炭工业出版社,1996.
[10]沈继方,于青春,胡章喜编.矿床水文地质学[M].武汉：中国地质大学出版社.1992,10.
[11]李定龙编.皖北奥陶系古岩溶及其环境地球化学特征研究[M].北京：石油工业出版社,2001,12.
[12]马国哲.平凉市灰岩岩溶水赋存规律探讨[J].甘肃地质,2001,10(1)：63-68.
[13]李文东,张新霞,唐德才等.藤县煤田腾北矿区奥灰富水性研究[J].江苏煤炭,2003,3,13-14.
[14]王延福,靳德武,曾艳京等.岩溶煤矿矿井煤层底板突水非线性预测方法研究[J].中国岩溶,1998,17(1)：57-66.
[15]王延福,庞西岐,靳德武等.岩溶矿井煤层底板突水的非线性动力学模型[J].中国岩溶,2000,19(1)：81-90.
[16]王延福,庞西岐,靳德武等.岩溶煤矿矿井煤层底板突水预测的相空间重构[J].中国岩溶,1999,18(3)：119-127.
[17]王延福,靳德武,曾艳京等.岩溶矿井煤层底板突水系统的非线性特征初步分析[J].中国岩溶,1998,17(4)：331-341.
[18]周信鲁,石亚东,王式成.井灌区地下水动态模拟的面状井系——水均衡综合模型[J].西北水资源与工程,2000,11(2)：24-28.
[19]冯德顺.岩溶地区水库渗漏的水均衡分析方法[J].人民长江,1996,27(7)：28-29.
[20]A.Srinivasulu, Ch.Sujani Rao, G.V.Lakshmi, T.V.Satyanarayana, J.Boonstra. Model Studies on Salt and Water Balances at Konanki Pilot Area, Andhra Pradesh, India. Irrigation and Drainage Systems,2004,18(1):1-17.
[21]卢文喜.地下水系统的模拟预测和优化管理[M].北京：科学出版社,1999.
[22]沈慧珍.宿南矿区第四含水层水文地质特征研究[D].淮南：安徽理工大学,2005.
[23]张永波,时红.灰色关联分析在地下水动态类型划分中的应用[J].地下水,1994,16(3)：136-138.
[24]王大纯,张人权等.水文地质学基础[M].地质出版社,1995.
[25]丁继红.国外地下水模拟软件的发展现状与趋势[J].勘察科学技术,2002,1：37-42.
[26]郝治福,康绍忠.地下水系统数值模拟的研究现状和发展趋势[J].水利水电科技进展,2006,26(1)：77-81.
[27]F.A.Abdulla, M.A.Al-Khatib,Z.D.Al-Ghazzawi. Development of groundwater for the Azraq Basin, Jordan. Environmental Geology,2000,40(1-2):1-18.
[28]J.W.Mercurio, Milovan S.Beljin, J.Barry J.Maynard. Groundwater models and wellfield management:a case study. Environmental Engineering and Policy,1999,(3):155-164.
[29]Nguyen Cao Don, Hiroyuki Araki, Hiroyuki Yamanishi, Kenichi Koga. Simulation of groundwater flow and environmental effects resulting from pumping. Environmental Geology,2005,47(3):361-374.
[30]Todd W.Rayne,Kenneth R.Bradbury,Maureen A.Muldoon. Delineation of capture zones for municipal wells in fractured dolomite, Sturgeon Bay, Wisconsin, USA, Hydrogeology Journal,2001,9(5):432-450.
[31]Jonathan Levy, Gordon Chesters, Daniel P.Gustafson, et al. Assessing aquifer susceptibility to and severity of atrazine contamination at a field site in south-central Wisconsin, USA.Hydrogeology Journal,1998,6(4):483-499.
[32]冯利军,李克生,邵改群.具有线性功能函数的神经元在矿井识别中的应用[J]工业安全与环保,2005,31(2)：27-29.
[33]徐忠杰,杨永国,汤琳.神经网络在矿井水源判别中的应用[J].煤炭安全,2007,2：4-6.
[34]杨永国,黄福臣.非线性方法在矿井突水水源判别中的应用研究[J].中国矿业火学学报,2007,36(3)：283-286.
[35]李定龙.应用神经网络对煤矿突水预测评价[J].Coal,1998,6：18-19.
[36]魏永强,梁化强,任印国,等.神经网络在判别煤矿突水水源中的应用[J].江苏地质,2004,28(1)：36-38.
[37]李栋臣.白庙煤矿主要含水层水化学特征及突水水源的识别[J].中国岩溶,1995,14(4)：295-304.
[38]岳梅.判断矿井突水水源灰色系统关联分析的应用[J].煤炭科学技术,2002,30(4)：37-39.
[39]杨永国,李宾亭.用数学地质判别鹤壁矿务局矿井水源[J].中国煤田地质,1995,7(4)：66-70.
[40]李世峰,李耀华.模糊综合判别矿井突(涌)水水源[J].煤炭工程,2006,9：71-73.
[41]葛中华,沈文,贝怀成.徐州某矿井奥陶系灰岩含水层上开采矿井突水的水文地质初步研究[J].江苏地质,1994,18(2)：91-96.
[42]张许良,张子戌,彭苏萍.数量化理论在矿井突(涌)水水源判别中的应用[J].中国矿业大学学报,2003,32(3)：251-254.
[43]王广生,王秀辉,李竞生,等.平顶山矿区矿井突(涌)水水源判别模式[J].煤田地质与勘探,1998,26(6)：47-50.
[44]殷晓曦,许光泉,桂和荣,等.系统聚类逐步判别法对皖北矿区突水水源的分析[J].煤田地质与勘探,2006,34(2)：58-61.
[45]孙亚军,杨国勇,郑琳.基于GIS的矿井突水水源判别系统的研究[J].煤田地质与勘探,2007,35(2)：34-37.
[46]桂和荣,许光泉,宋晓梅.桃园煤矿矿井水充水条件及充水水源识别[J].淮南职业技术学院学报,2003,4(3)：85-86.
[47]许光泉,桂和荣,吴基文.煤矿底板突水分析及底板水防治[J].地下水,2001,23(3)：141-143.
[48]王希良,彭苏萍,郑世书.深部煤层开采高承压水突水预报及控制[J].辽宁工程技术大学学报,2004,23(6)：758-760.
[49]卜昌森,张希诚,尹万才,曲修术.“华北型”煤田岩溶水害及防治现状[J].地质论评,2001,47(4)：405-410.
[50]周维垣,岩体力学数值计算方法的现状与展望[J].岩石力学与工程学报,1993(1)：(84-88).
[51]FLAC-3D.Fast langrangian analysis of continua in 3 Dimensions,Version 2.1,user manuals[M].Minneapolis:Itasca Consulting Group,Inc,2002.
[52]Bai M. and Elsworth D.,1994. Modeling of subsidence and stress-dependent hydraulic conductivity for intace and fractured porous media. Rock. Mech. & Rock Engry. Vol.4, 235-251.
[53]Bandis, S.C., Lumsden, A.C. and Barton, N.R.,1983. Fundamentals of rock joint deformation. Int. J. Rock Mech. Min. Sci.& Geomech. Abstr..20:249-268.
[54]Barton, N. & Bandis, S.,1982. Effects of block size on the shear behavior of jointed rock. In Proc.23rd U.S. Symp. Rock Mech., p739-760, Berkeley, California.
[55]张文泉,刘伟韬,张红日等.煤层底板岩层阻水能力及其影响因素的研究[J],岩土力学,1998,19(4).
[56]吴基文.杨庄煤矿六煤底板采动效应研究[J],岩土力学,2003,19(4).
[57]李俊才.软土深基坑开挖现场测试及三维数值模拟研究：[D].成都：成都理工大学,2001.
[58]蒋亚萍,陈余道Modflow——套水文地质学实用软件[J].广‘西地质,1999,12(3)：75-78.
[59]魏云杰,许模,刘健Visual Modflow软件及其在砂岩刑铀矿成矿水文地质条件研究中的应用潜力分析[J].铀矿地质,2003,19(1)：53-57.
[60]武强等.可视化地下水模拟评价新型软件系统(Visual Modflow)与矿井防治水[J].煤炭科学技术,2000,28(2)：18-20.
[61]郭卫星,卢国正编译,原著Michael G.McDonald和Arlen W.Harbaugh,MODFLOW:模块化三维有限差分地下水流动模型.
[62]Ching-Pin Tung,Chung-Che Tang,Yu-Pin Lin. Improving groundwater-flow modeling using optimal zoning methods. Environmental Geology,2003,44 (6):627-638.
[63]葛伟亚.盐城市地下水系统三维数值模拟：[D].南京：河海大学,2004.
[64]陈崇希,唐仲华.地下水流动问题数值方法[M].武汉：中国地质大学出版社,1990.
[65]郭惟嘉等.底板突水系数概念及其应用[J].河北煤炭,1989(2)：56-60.
[66]淮北矿务局杨庄矿Ⅱ633工作面底板原位测试及其突水可能性预测.煤科总院西安分院,1998.
[67]管恩太,武强.矿井涌水量预测评述[J].中州煤炭,2005,27(1)：7-8.
[68]李定龙,周治安.临涣矿区底含水化学特征及其形成作用探讨[J].煤田地质与勘探,1994,22,(4)：37-41.
[69]沈照理.水文地球化学基础[M].北京：煤炭工业出版社,1992.
[70]许福美.田螺形矿井地下水化学特征分析与应用[J].西部探矿工程,2005,(5)：69-71.
[71]汪世花.鹤壁矿区各含水层水化学特征与水源判别初探[J].中州煤炭,1998,(5)：30-31.
[72]段佐亮.模式识别*逐步判别分析法在农业环境质量分级中的应用[J].农业环境与发展,1995,(3)：13-19.
[73]高艳秋.多组逐步判别分析在矿井水源判别中的应用[J].北京工业职业技术学院学报,2007,6(2)：11-14.
[74]苏金明.统计软件spss系列应用实战篇[M].北京：电子工业出版社,2002：318.
[75]杜文堂.断层防水煤柱可靠度分析[J].煤田地质与勘探,2001,29(1)：34-36.
[2]李金凯等编.矿井岩溶水防治[M].北京：煤炭工业出版社,1990,3.
[3]煤炭部生产协调司.国有重点煤矿井田内受小煤矿开采影响安全生产情况的调查,1995,7.
[4]李定龙,周治安,王桂梁.马家沟灰岩(古)岩溶研究中的若干问题探讨[J].地质科技情报,1997,16(1).
[5]Reibiec M.S.,Hydrofracturing of rock as a method of water,mudmand gas inrush hazards in underground coal mining,4th IMWA,1991,1 (Yugoslavia).
[6]胡春玲.里彦矿灰岩水对下组煤安全开采影响的研究[D].山东科技大,2005.
[7]高延法.底板突水的研究途径与突水优势面[J].山东矿业学院学报,1994,13(增).
[8]王作宇,刘鸿泉.承压水上采煤[M].北京：煤炭工业出版社,1993.
[9]王永红,沈文.中国煤矿水害预防及治理[M].北京：煤炭工业出版社,1996.
[10]沈继方,于青春,胡章喜编.矿床水文地质学[M].武汉：中国地质大学出版社.1992,10.
[11]李定龙编.皖北奥陶系古岩溶及其环境地球化学特征研究[M].北京：石油工业出版社,2001,12.
[12]马国哲.平凉市灰岩岩溶水赋存规律探讨[J].甘肃地质,2001,10(1)：63-68.
[13]李文东,张新霞,唐德才等.藤县煤田腾北矿区奥灰富水性研究[J].江苏煤炭,2003,3,13-14.
[14]王延福,靳德武,曾艳京等.岩溶煤矿矿井煤层底板突水非线性预测方法研究[J].中国岩溶,1998,17(1)：57-66.
[15]王延福,庞西岐,靳德武等.岩溶矿井煤层底板突水的非线性动力学模型[J].中国岩溶,2000,19(1)：81-90.
[16]王延福,庞西岐,靳德武等.岩溶煤矿矿井煤层底板突水预测的相空间重构[J].中国岩溶,1999,18(3)：119-127.
[17]王延福,靳德武,曾艳京等.岩溶矿井煤层底板突水系统的非线性特征初步分析[J].中国岩溶,1998,17(4)：331-341.
[18]周信鲁,石亚东,王式成.井灌区地下水动态模拟的面状井系——水均衡综合模型[J].西北水资源与工程,2000,11(2)：24-28.
[19]冯德顺.岩溶地区水库渗漏的水均衡分析方法[J].人民长江,1996,27(7)：28-29.
[20]A.Srinivasulu, Ch.Sujani Rao, G.V.Lakshmi, T.V.Satyanarayana, J.Boonstra. Model Studies on Salt and Water Balances at Konanki Pilot Area, Andhra Pradesh, India. Irrigation and Drainage Systems,2004,18(1):1-17.
[21]卢文喜.地下水系统的模拟预测和优化管理[M].北京：科学出版社,1999.
[22]沈慧珍.宿南矿区第四含水层水文地质特征研究[D].淮南：安徽理工大学,2005.
[23]张永波,时红.灰色关联分析在地下水动态类型划分中的应用[J].地下水,1994,16(3)：136-138.
[24]王大纯,张人权等.水文地质学基础[M].地质出版社,1995.
[25]丁继红.国外地下水模拟软件的发展现状与趋势[J].勘察科学技术,2002,1：37-42.
[26]郝治福,康绍忠.地下水系统数值模拟的研究现状和发展趋势[J].水利水电科技进展,2006,26(1)：77-81.
[27]F.A.Abdulla, M.A.Al-Khatib,Z.D.Al-Ghazzawi. Development of groundwater for the Azraq Basin, Jordan. Environmental Geology,2000,40(1-2):1-18.
[28]J.W.Mercurio, Milovan S.Beljin, J.Barry J.Maynard. Groundwater models and wellfield management:a case study. Environmental Engineering and Policy,1999,(3):155-164.
[29]Nguyen Cao Don, Hiroyuki Araki, Hiroyuki Yamanishi, Kenichi Koga. Simulation of groundwater flow and environmental effects resulting from pumping. Environmental Geology,2005,47(3):361-374.
[30]Todd W.Rayne,Kenneth R.Bradbury,Maureen A.Muldoon. Delineation of capture zones for municipal wells in fractured dolomite, Sturgeon Bay, Wisconsin, USA, Hydrogeology Journal,2001,9(5):432-450.
[31]Jonathan Levy, Gordon Chesters, Daniel P.Gustafson, et al. Assessing aquifer susceptibility to and severity of atrazine contamination at a field site in south-central Wisconsin, USA.Hydrogeology Journal,1998,6(4):483-499.
[32]冯利军,李克生,邵改群.具有线性功能函数的神经元在矿井识别中的应用[J]工业安全与环保,2005,31(2)：27-29.
[33]徐忠杰,杨永国,汤琳.神经网络在矿井水源判别中的应用[J].煤炭安全,2007,2：4-6.
[34]杨永国,黄福臣.非线性方法在矿井突水水源判别中的应用研究[J].中国矿业火学学报,2007,36(3)：283-286.
[35]李定龙.应用神经网络对煤矿突水预测评价[J].Coal,1998,6：18-19.
[36]魏永强,梁化强,任印国,等.神经网络在判别煤矿突水水源中的应用[J].江苏地质,2004,28(1)：36-38.
[37]李栋臣.白庙煤矿主要含水层水化学特征及突水水源的识别[J].中国岩溶,1995,14(4)：295-304.
[38]岳梅.判断矿井突水水源灰色系统关联分析的应用[J].煤炭科学技术,2002,30(4)：37-39.
[39]杨永国,李宾亭.用数学地质判别鹤壁矿务局矿井水源[J].中国煤田地质,1995,7(4)：66-70.
[40]李世峰,李耀华.模糊综合判别矿井突(涌)水水源[J].煤炭工程,2006,9：71-73.
[41]葛中华,沈文,贝怀成.徐州某矿井奥陶系灰岩含水层上开采矿井突水的水文地质初步研究[J].江苏地质,1994,18(2)：91-96.
[42]张许良,张子戌,彭苏萍.数量化理论在矿井突(涌)水水源判别中的应用[J].中国矿业大学学报,2003,32(3)：251-254.
[43]王广生,王秀辉,李竞生,等.平顶山矿区矿井突(涌)水水源判别模式[J].煤田地质与勘探,1998,26(6)：47-50.
[44]殷晓曦,许光泉,桂和荣,等.系统聚类逐步判别法对皖北矿区突水水源的分析[J].煤田地质与勘探,2006,34(2)：58-61.
[45]孙亚军,杨国勇,郑琳.基于GIS的矿井突水水源判别系统的研究[J].煤田地质与勘探,2007,35(2)：34-37.
[46]桂和荣,许光泉,宋晓梅.桃园煤矿矿井水充水条件及充水水源识别[J].淮南职业技术学院学报,2003,4(3)：85-86.
[47]许光泉,桂和荣,吴基文.煤矿底板突水分析及底板水防治[J].地下水,2001,23(3)：141-143.
[48]王希良,彭苏萍,郑世书.深部煤层开采高承压水突水预报及控制[J].辽宁工程技术大学学报,2004,23(6)：758-760.
[49]卜昌森,张希诚,尹万才,曲修术.“华北型”煤田岩溶水害及防治现状[J].地质论评,2001,47(4)：405-410.
[50]周维垣,岩体力学数值计算方法的现状与展望[J].岩石力学与工程学报,1993(1)：(84-88).
[51]FLAC-3D.Fast langrangian analysis of continua in 3 Dimensions,Version 2.1,user manuals[M].Minneapolis:Itasca Consulting Group,Inc,2002.
[52]Bai M. and Elsworth D.,1994. Modeling of subsidence and stress-dependent hydraulic conductivity for intace and fractured porous media. Rock. Mech. & Rock Engry. Vol.4, 235-251.
[53]Bandis, S.C., Lumsden, A.C. and Barton, N.R.,1983. Fundamentals of rock joint deformation. Int. J. Rock Mech. Min. Sci.& Geomech. Abstr..20:249-268.
[54]Barton, N. & Bandis, S.,1982. Effects of block size on the shear behavior of jointed rock. In Proc.23rd U.S. Symp. Rock Mech., p739-760, Berkeley, California.
[55]张文泉,刘伟韬,张红日等.煤层底板岩层阻水能力及其影响因素的研究[J],岩土力学,1998,19(4).
[56]吴基文.杨庄煤矿六煤底板采动效应研究[J],岩土力学,2003,19(4).
[57]李俊才.软土深基坑开挖现场测试及三维数值模拟研究：[D].成都：成都理工大学,2001.
[58]蒋亚萍,陈余道Modflow——套水文地质学实用软件[J].广‘西地质,1999,12(3)：75-78.
[59]魏云杰,许模,刘健Visual Modflow软件及其在砂岩刑铀矿成矿水文地质条件研究中的应用潜力分析[J].铀矿地质,2003,19(1)：53-57.
[60]武强等.可视化地下水模拟评价新型软件系统(Visual Modflow)与矿井防治水[J].煤炭科学技术,2000,28(2)：18-20.
[61]郭卫星,卢国正编译,原著Michael G.McDonald和Arlen W.Harbaugh,MODFLOW:模块化三维有限差分地下水流动模型.
[62]Ching-Pin Tung,Chung-Che Tang,Yu-Pin Lin. Improving groundwater-flow modeling using optimal zoning methods. Environmental Geology,2003,44 (6):627-638.
[63]葛伟亚.盐城市地下水系统三维数值模拟：[D].南京：河海大学,2004.
[64]陈崇希,唐仲华.地下水流动问题数值方法[M].武汉：中国地质大学出版社,1990.
[65]郭惟嘉等.底板突水系数概念及其应用[J].河北煤炭,1989(2)：56-60.
[66]淮北矿务局杨庄矿Ⅱ633工作面底板原位测试及其突水可能性预测.煤科总院西安分院,1998.
[67]管恩太,武强.矿井涌水量预测评述[J].中州煤炭,2005,27(1)：7-8.
[68]李定龙,周治安.临涣矿区底含水化学特征及其形成作用探讨[J].煤田地质与勘探,1994,22,(4)：37-41.
[69]沈照理.水文地球化学基础[M].北京：煤炭工业出版社,1992.
[70]许福美.田螺形矿井地下水化学特征分析与应用[J].西部探矿工程,2005,(5)：69-71.
[71]汪世花.鹤壁矿区各含水层水化学特征与水源判别初探[J].中州煤炭,1998,(5)：30-31.
[72]段佐亮.模式识别*逐步判别分析法在农业环境质量分级中的应用[J].农业环境与发展,1995,(3)：13-19.
[73]高艳秋.多组逐步判别分析在矿井水源判别中的应用[J].北京工业职业技术学院学报,2007,6(2)：11-14.
[74]苏金明.统计软件spss系列应用实战篇[M].北京：电子工业出版社,2002：318.
[75]杜文堂.断层防水煤柱可靠度分析[J].煤田地质与勘探,2001,29(1)：34-36.