成都市关口水库右岸邻谷渗漏分析
|
摘要
水库渗漏是水库建设中常见的地质问题之一。影响水库渗漏的因素很多，诸如地形地貌、地层岩性、地质构造、岩溶发育程度和水文地质条件等。其中，最关键的两个因素是由库区的地质构造和地层岩性等决定的渗漏通道，以及由水文地质条件等决定的水动力条件，构成水库渗漏的充分必要条件。通常情况下，水库的渗漏通道是指存在沟通水库内外的强渗透介质，最为常见的是岩溶含水介质，以及断层破碎带、古河道等，而人为活动(如采煤井巷、隧道工程、勘探井孔等)形成的通道则较为少见。水库渗漏的水动力条件是指水库与邻谷、以及其间的地下分水岭之间的水动力关系，当水库蓄水高程超越地下分水岭高程时，水库发生渗漏的水动力条件这一充分条件亦成立。
拟建的关口水库库区位于成都平原北部边缘低山区，建好后将作为成都市的后备水源。水库建在沱江水系的湔江干流上游河段，库区河谷开阔，岸坡地形坡度较缓。由于库区左岸多为高大雄厚的山体，发生的渗漏的可能性很小，本文主要讨论了库区右岸的渗漏。首先分析了库区右岸水库渗漏的基本条件，在天然条件下顺层发育的裂隙含水层构成了潜在的渗漏通道，而煤矿开采后的井下巷道形成了典型的人为活动渗漏通道。而且水库蓄水后库水位和库区内的煤矿、邻谷煤矿开采水平之间存在一个巨大的水头差，因而可以认为产生水库渗漏的水动力条件也是存在的。因此水库渗漏通道和水动力条件都成立；进而分析了水库渗漏的基本模式以及相关的处理措施，对于右岸单薄垭口和煤矿对水库渗漏的影响做了具体详细的分析；最后，利用Visual Modflow软件建立了符合实际情况的三维地下水流模型，实现了库区右岸在水库蓄水前后煤矿开采情况下以及水库蓄水后煤矿(主要是跃进煤矿)不开采情况下地下水流场变化的模拟，得出了研究区在不同时间段、不同情况下的地下水流场图。在此基础上，通过对不同时间和情况下的地下水流场的模拟，研究了区内在煤矿开采前后及水库蓄水后地下水位的变化情况，为研究区内煤矿在未来的水库蓄水后的排水通道和制定防治水措施提供了科学的依据，同时也为将要修建的关口水库右岸邻谷的渗漏影响评价奠定了基础。
Reservoir leakage is a common problem of reservoir construction and engineering geology. Reservoir leakage is affected by a number of factors, such as topography, lithology, geological structure, karst development and hydrogeological conditions. In the most factors of influencing reservoir leakage, two of the most crucial factors is the leakage passageway in the project area decided by the geological structure and lithology and hydrodynamic conditions decided by hydrogeology. Normally, the reservoir is leaking channels of communication exist within and outside the strong reservoir permeability media, is the most common medium of karst aquifers and fractured fault zone, ancient river, and human activities (such as coal mining shaft, Tunnel project, exploration bores, etc.) of the channel is very rare. Reservoir water leakage is dynamic conditions with the adjacent Valley Reservoir, as well as during the watershed between the underground water power, When the reservoir water elevation beyond the underground watershed elevation, Reservoir leakage of the hydrodynamic conditions of the full terms and conditions are set up.
The Guan kou reservoir locates in the northern edge of the Chengdu Plain at the low mountains, will be completed as a backup water source of Chengdu .Reservoir built On the Jian Jiang River upper reaches of Tuojiang water system, open valley area, the terrain slopes gentler slope. Because most of the left bank of the reservoir area are tall mountains, the possibility of leakage is small, This paper discusses the right bank of the reservoir leakage. First, we analyzed the basic leakage conditions,Under natural conditions in bedding development of the unconfined aquifer constitutes a natural channel leakage, coal mining and underground roadway after the formation of artificial channel leakage. After the dam water level and the level of coal mining between the head of a huge difference. This reservoir leakage of the hydrodynamic conditions, the reservoir channel leakage power and water conditions have been established; Reservoir then expounded the basic pattern of leakage and related treatment, on the right bank of thin Yakou and coal reservoir leakage of the influence of the detailed analysis; Finally, used the current prevailing international software Visual Modflow to build a realistic three-dimensional model of groundwater flow type, Implementation of the right bank of the reservoir area before and after the reservoir was filled, where coal mining and coal mine after water storage reservoirs is not exploitation of groundwater flow changes in the simulation of the study area at different times, under different circumstances groundwater flow graph. On this basis, through different times and circumstances of groundwater flow simulation, of coal mining in the region before and after the reservoir was filled, the water table changes, Mine study area for the future of water after the drainage channel and the development of measures to control water and provides a scientific basis, also for the construction of the gateway to the right bank of the reservoir Valley leakage effects laid the foundation for evaluation.
拟建的关口水库库区位于成都平原北部边缘低山区，建好后将作为成都市的后备水源。水库建在沱江水系的湔江干流上游河段，库区河谷开阔，岸坡地形坡度较缓。由于库区左岸多为高大雄厚的山体，发生的渗漏的可能性很小，本文主要讨论了库区右岸的渗漏。首先分析了库区右岸水库渗漏的基本条件，在天然条件下顺层发育的裂隙含水层构成了潜在的渗漏通道，而煤矿开采后的井下巷道形成了典型的人为活动渗漏通道。而且水库蓄水后库水位和库区内的煤矿、邻谷煤矿开采水平之间存在一个巨大的水头差，因而可以认为产生水库渗漏的水动力条件也是存在的。因此水库渗漏通道和水动力条件都成立；进而分析了水库渗漏的基本模式以及相关的处理措施，对于右岸单薄垭口和煤矿对水库渗漏的影响做了具体详细的分析；最后，利用Visual Modflow软件建立了符合实际情况的三维地下水流模型，实现了库区右岸在水库蓄水前后煤矿开采情况下以及水库蓄水后煤矿(主要是跃进煤矿)不开采情况下地下水流场变化的模拟，得出了研究区在不同时间段、不同情况下的地下水流场图。在此基础上，通过对不同时间和情况下的地下水流场的模拟，研究了区内在煤矿开采前后及水库蓄水后地下水位的变化情况，为研究区内煤矿在未来的水库蓄水后的排水通道和制定防治水措施提供了科学的依据，同时也为将要修建的关口水库右岸邻谷的渗漏影响评价奠定了基础。
Reservoir leakage is a common problem of reservoir construction and engineering geology. Reservoir leakage is affected by a number of factors, such as topography, lithology, geological structure, karst development and hydrogeological conditions. In the most factors of influencing reservoir leakage, two of the most crucial factors is the leakage passageway in the project area decided by the geological structure and lithology and hydrodynamic conditions decided by hydrogeology. Normally, the reservoir is leaking channels of communication exist within and outside the strong reservoir permeability media, is the most common medium of karst aquifers and fractured fault zone, ancient river, and human activities (such as coal mining shaft, Tunnel project, exploration bores, etc.) of the channel is very rare. Reservoir water leakage is dynamic conditions with the adjacent Valley Reservoir, as well as during the watershed between the underground water power, When the reservoir water elevation beyond the underground watershed elevation, Reservoir leakage of the hydrodynamic conditions of the full terms and conditions are set up.
The Guan kou reservoir locates in the northern edge of the Chengdu Plain at the low mountains, will be completed as a backup water source of Chengdu .Reservoir built On the Jian Jiang River upper reaches of Tuojiang water system, open valley area, the terrain slopes gentler slope. Because most of the left bank of the reservoir area are tall mountains, the possibility of leakage is small, This paper discusses the right bank of the reservoir leakage. First, we analyzed the basic leakage conditions,Under natural conditions in bedding development of the unconfined aquifer constitutes a natural channel leakage, coal mining and underground roadway after the formation of artificial channel leakage. After the dam water level and the level of coal mining between the head of a huge difference. This reservoir leakage of the hydrodynamic conditions, the reservoir channel leakage power and water conditions have been established; Reservoir then expounded the basic pattern of leakage and related treatment, on the right bank of thin Yakou and coal reservoir leakage of the influence of the detailed analysis; Finally, used the current prevailing international software Visual Modflow to build a realistic three-dimensional model of groundwater flow type, Implementation of the right bank of the reservoir area before and after the reservoir was filled, where coal mining and coal mine after water storage reservoirs is not exploitation of groundwater flow changes in the simulation of the study area at different times, under different circumstances groundwater flow graph. On this basis, through different times and circumstances of groundwater flow simulation, of coal mining in the region before and after the reservoir was filled, the water table changes, Mine study area for the future of water after the drainage channel and the development of measures to control water and provides a scientific basis, also for the construction of the gateway to the right bank of the reservoir Valley leakage effects laid the foundation for evaluation.
引文
[1] 张津生，邵天星．广东省岭澳水库渗漏特点及原因分析[J]．水利水电技术，1998，(8)：31-32
[2] 盛金保，谢晓华，李雷．株树桥水库严重渗漏原因分析[J]．水利水运工程学报，2003，(4)：25-30
[3] 周洪文等．马畔塘水库岩溶渗漏研究[J]．中国岩溶，1996，15(4)
[4] 邹成杰．水库岩溶渗漏地质模型和数学模型的初步研究中国岩溶[J]．1999 9(3) -231-240
[5] 费英烈，邹成杰，等．贵州岩溶地区水库坝址渗漏问题的初步研究[J]．中国岩溶，1984，3(2)
[6] 许模．大型水电工程中主要工程水文地质问题系统研究：[D]．成都理工大学博士学位论文，1998
[7] 沈玉玲，蒋甫男，李西平．小南海水库渗漏的地质背景及处理[J]．水利水电科技进展，1998，18(4)：37-40
[8] 王庆苗．磨子潭水库渗漏分析[J]．大坝与安全，2000，(4)：43-45
[9] 许模，张强，黄润秋．洛湛线袍子岭隧道隧址区水文地质条件及涌水量预测综合研究[M]：成都理工大学，2002
[10] Wang, H. F&J. Peraric. Introduction to Groundwater Modeling. W. H. Freemanand company, 1982.
[11] Illangasekare, T. &H. J. Morel-seytoux. Stream-Aquifer Influence Coefficients as Tools for Simulation and Management. Water Resources Research 1982, 18, (1).
[12] Eddebbarh, Al, Carlson, Fritz, Hemker, C. J. Groundwater modeling using an integrated large-capacity finite-element model for multiple aquifer groundwater flow. Proceedings of the 1996 Symposium on Subsurface Fluid Flow Modeling, ]un 22-231996, 1: (22-23).
[13] Hsieh PA, Freckleton. Documentation of a computer program to simulatehorizontal-flow barriers using the U.S. Geological Survey's modular three-dimensional finite-difference ground-water flow model[R]. U.S. Geological Survey Open-File Report, 1993, 92-477.
[14] Cooley R L. Incorporation on parameters into nonlinear regression groundwater flow model. 1, Theory. Water Resources Research, 1982, 18(4): 965—976
[15] Cooley R L. Incorporation on parameters into nonlinear regression groundwater flow model. 2, Application, Water Resources Research, 1983, 19(3): 662-676
[16] Xue Y Q, Xie C H, Wu J C. A numerical model of the seawater intrusion in coastal aquifer. In: Xue Y Q, Bear J, eds. Proceeding of the International Conference on Modeling Groundwater Flow and Pollution. Nanjing: Nanjing University Press, 1991. 13-20
[17] 孙钠正，地下水流的数学模型利数值方法[M]．地质出版社，1981
[18] 李俊亭，地下水数值模拟[M]．北京：地质出版社，1989
[19] 卫爱民，陈晓霞．复杂岩溶地区麒麟观水库渗漏分析研究[J]．低温建筑技术，2004，(4)：33-25
[20] 刘加海，付彦，扬名玖．红兴水库渗漏分析及防渗处理方案研究[J]．水利水电技术，2006，(5)：52-49．
[21] 周志芳，王锦国．裂隙介质水动力学[M]．中国水利水电出版社，北京
[22] 邹成杰．水利水电岩溶工程地质[M]．北京：水利水电出版社，1993
[23] 王大纯，张仁权，史毅红．水文地质学基础[M]．地质出版社，1986
[24] 张新兴．浩坤水电站水库喀斯特渗漏分析[J]．广西水利水电，2003，(3)：19-23
[25] 邹成杰，何宇彬．喀斯特地貌发育的时空演化问题初论[J]．中国岩溶，1995，14(3)：51
[26] 郭东屏，张石峰．渗流理论基础[M]．陕西科学技术出版社，1994
[27] 毛昶熙．渗流计算分析与控制[M]．北京：中国水利水电出版社，2003：64-82
[28] 赵青．吴家庄水库库区岩溶渗漏分析[J]．山西建筑，2005(2)：57-62
[29] 田茂中，叶明等．毛家河水电站库区岩溶发育特征及邻谷渗漏分析[J]．贵州水力发电，2006(6)：18-06
[30] 于汝华．阿岗水库库区岩溶渗漏研究[J]．人民长江，2005(9)：06-02
[31] 张永波．水工环研究的现状与趋势[M]．地质出版社，2001．
[32] H．x．加维奇 著 杨守廉 译．水文地质学中应用模拟的理论与实践[M]．地质出版社，1987
[33] 薛禹群，谢春红．水文地质学的数值法[M]．煤炭工业出版社，1979
[34] 高正夏．地下水资源评价与计算[J]．河海大学土木工程学院地岩系，1997
[35] 骆祖江，付延玲．模型技术在地下水系统动态预测中的应用[J]．工程勘察，1997(5)：29-33
[36] 朱学愚、徐绍辉、吴春寅，等．城门山铜矿地下水数学模型和疏干预测[J]．南京大学学报，1995．31(1)：109-120
[37] 徐恒力，等．水资源开发与保护[M]．北京，地质出版社，2001
[38] 袁宏源，邵东国，郭宗楼．水资源系统分析理论与应用[M]．武汉水利电力出版社，2000
[39] 叶守则．水文水利计算[M]．水利电力出版社，1991
[40] 周念清，朱蓉，朱学愚．MODFLOW在宿迁市地下水资源评价中的应用[J]．水文地质与工程地质，2000(6)
[41] 王宏等．MODFLOW在华北平原地下水库模拟中的应用[J]．世界地质，2003(3)
[42] 陈喜，陈询洪．美国Sand Hills地区地下水数值模拟及水量平衡分析[J]．水科学进展，2004(1)
[43] 陈江中．地下水资源系统分析[M]．中国矿业大学出版社，1992
[44] 林学钮，侯印伟．地下水水量水质模拟及管理程序集[M]．吉林科学技术出版社，1988
[45] 郭卫星，卢国平．MODFLOW三维有限差分地下水流模型[M]．南京大学地球科学系，1998
[46] 武强，董东林等．水资源评价的可视化专业软件(Visual Modflow)与应用潜力[J]．水文地质工程地质，1999(5)：21-23
[47] 陈玉田．偏微分方程数值解法[M]．河海大学出版社，1999
[48] 陈崇喜，唐仲华．地下水流动问题数值方法[M]．中国地质大学出版社，1990
[49] 陈锁忠．MODFLOW三维有限差分地下水渗流模型的原理[J]．江苏地矿信息，1999(1)：21-23
[50] 钱家忠，吴剑锋等．徐州市张集水源地裂隙岩溶水三维等参有限元数值模拟[J]．水利学报，2003(3)：37—1
[51] 宋西平，寇广潮．地下水数值模拟面参数和水文参数研究的重要性[J]．地下水，2001(2)：63-64
[52] 刘振英，李文鹏，孙鸣．水资源实时模拟管理系统及其在山东临淄地区的应用[M]．地质出版社，1996
[53] 国家技术监督局．中华人民共和国国家标准(GB／T14497-93)地下水资源管理模型工作技术要求[M]．北京，1994
[54] 朱学愚，谢春红．地下水运移模型[M]．中国建筑工业出版社，1990
[55] 朱学愚，钱孝星，刘新仁．地下水资源评价[M]．南京大学出版社，1987
[56] 张蔚棒．地下水非稳定流计算和地下水资源评价[M]．科学出版社，1983
[57] 张宏仁等．地下水非稳定流理论的发展和应用M]．地质出版社，1975
[58] 束龙仓，朱元生，地下水资源评价中的不确定性因素分析[J]，水文地质工程地质，2000，27(6)：6-8
[59] 束龙仓，朱元生，孙庆义，等．地下水资源评价结果的可靠性探讨[J]．水科学进展，2000，11(1)：21-24
[60] 束龙仓，朱元生．晋祠泉域地下水开采决策的风险分析[J]．河海大学学报，2000，28(6)：90-93
[61] 薛禹群，吴吉春．面临21世纪的中国地下水模拟问题[J]．水文地质工程地质，1999，(5)：1-3
[62] 束龙仓．山东省济宁市地下水资源管理模型研究[M]．建设部综合勘察研究院，1990
[63] 钱孝星．水文地质计算[M]．水利电力出版社，1995
[64] 卢文喜．地下水运动数值模拟过程中边界条件问题探讨[J]．水利学报，2003(3)：34-36
[65] 洪景涛，刘翔，严涛．淄博市大武水源地地下水流场数值模型研究[J]．水文地质学报，2003(4)：7-10
[66] 蒋亚萍，陈余道．MODFLOW——套水文地质学实用计算软件[J]．广西地质，1999(3)：75-78
[67] 邵景力，崔亚莉等．黄河下游影响带(河南段)三维地下水流数值模拟模型及其应用[J]．吉林大学学报(地球科学版)，2003，33(1)：55-58
[2] 盛金保，谢晓华，李雷．株树桥水库严重渗漏原因分析[J]．水利水运工程学报，2003，(4)：25-30
[3] 周洪文等．马畔塘水库岩溶渗漏研究[J]．中国岩溶，1996，15(4)
[4] 邹成杰．水库岩溶渗漏地质模型和数学模型的初步研究中国岩溶[J]．1999 9(3) -231-240
[5] 费英烈，邹成杰，等．贵州岩溶地区水库坝址渗漏问题的初步研究[J]．中国岩溶，1984，3(2)
[6] 许模．大型水电工程中主要工程水文地质问题系统研究：[D]．成都理工大学博士学位论文，1998
[7] 沈玉玲，蒋甫男，李西平．小南海水库渗漏的地质背景及处理[J]．水利水电科技进展，1998，18(4)：37-40
[8] 王庆苗．磨子潭水库渗漏分析[J]．大坝与安全，2000，(4)：43-45
[9] 许模，张强，黄润秋．洛湛线袍子岭隧道隧址区水文地质条件及涌水量预测综合研究[M]：成都理工大学，2002
[10] Wang, H. F&J. Peraric. Introduction to Groundwater Modeling. W. H. Freemanand company, 1982.
[11] Illangasekare, T. &H. J. Morel-seytoux. Stream-Aquifer Influence Coefficients as Tools for Simulation and Management. Water Resources Research 1982, 18, (1).
[12] Eddebbarh, Al, Carlson, Fritz, Hemker, C. J. Groundwater modeling using an integrated large-capacity finite-element model for multiple aquifer groundwater flow. Proceedings of the 1996 Symposium on Subsurface Fluid Flow Modeling, ]un 22-231996, 1: (22-23).
[13] Hsieh PA, Freckleton. Documentation of a computer program to simulatehorizontal-flow barriers using the U.S. Geological Survey's modular three-dimensional finite-difference ground-water flow model[R]. U.S. Geological Survey Open-File Report, 1993, 92-477.
[14] Cooley R L. Incorporation on parameters into nonlinear regression groundwater flow model. 1, Theory. Water Resources Research, 1982, 18(4): 965—976
[15] Cooley R L. Incorporation on parameters into nonlinear regression groundwater flow model. 2, Application, Water Resources Research, 1983, 19(3): 662-676
[16] Xue Y Q, Xie C H, Wu J C. A numerical model of the seawater intrusion in coastal aquifer. In: Xue Y Q, Bear J, eds. Proceeding of the International Conference on Modeling Groundwater Flow and Pollution. Nanjing: Nanjing University Press, 1991. 13-20
[17] 孙钠正，地下水流的数学模型利数值方法[M]．地质出版社，1981
[18] 李俊亭，地下水数值模拟[M]．北京：地质出版社，1989
[19] 卫爱民，陈晓霞．复杂岩溶地区麒麟观水库渗漏分析研究[J]．低温建筑技术，2004，(4)：33-25
[20] 刘加海，付彦，扬名玖．红兴水库渗漏分析及防渗处理方案研究[J]．水利水电技术，2006，(5)：52-49．
[21] 周志芳，王锦国．裂隙介质水动力学[M]．中国水利水电出版社，北京
[22] 邹成杰．水利水电岩溶工程地质[M]．北京：水利水电出版社，1993
[23] 王大纯，张仁权，史毅红．水文地质学基础[M]．地质出版社，1986
[24] 张新兴．浩坤水电站水库喀斯特渗漏分析[J]．广西水利水电，2003，(3)：19-23
[25] 邹成杰，何宇彬．喀斯特地貌发育的时空演化问题初论[J]．中国岩溶，1995，14(3)：51
[26] 郭东屏，张石峰．渗流理论基础[M]．陕西科学技术出版社，1994
[27] 毛昶熙．渗流计算分析与控制[M]．北京：中国水利水电出版社，2003：64-82
[28] 赵青．吴家庄水库库区岩溶渗漏分析[J]．山西建筑，2005(2)：57-62
[29] 田茂中，叶明等．毛家河水电站库区岩溶发育特征及邻谷渗漏分析[J]．贵州水力发电，2006(6)：18-06
[30] 于汝华．阿岗水库库区岩溶渗漏研究[J]．人民长江，2005(9)：06-02
[31] 张永波．水工环研究的现状与趋势[M]．地质出版社，2001．
[32] H．x．加维奇 著 杨守廉 译．水文地质学中应用模拟的理论与实践[M]．地质出版社，1987
[33] 薛禹群，谢春红．水文地质学的数值法[M]．煤炭工业出版社，1979
[34] 高正夏．地下水资源评价与计算[J]．河海大学土木工程学院地岩系，1997
[35] 骆祖江，付延玲．模型技术在地下水系统动态预测中的应用[J]．工程勘察，1997(5)：29-33
[36] 朱学愚、徐绍辉、吴春寅，等．城门山铜矿地下水数学模型和疏干预测[J]．南京大学学报，1995．31(1)：109-120
[37] 徐恒力，等．水资源开发与保护[M]．北京，地质出版社，2001
[38] 袁宏源，邵东国，郭宗楼．水资源系统分析理论与应用[M]．武汉水利电力出版社，2000
[39] 叶守则．水文水利计算[M]．水利电力出版社，1991
[40] 周念清，朱蓉，朱学愚．MODFLOW在宿迁市地下水资源评价中的应用[J]．水文地质与工程地质，2000(6)
[41] 王宏等．MODFLOW在华北平原地下水库模拟中的应用[J]．世界地质，2003(3)
[42] 陈喜，陈询洪．美国Sand Hills地区地下水数值模拟及水量平衡分析[J]．水科学进展，2004(1)
[43] 陈江中．地下水资源系统分析[M]．中国矿业大学出版社，1992
[44] 林学钮，侯印伟．地下水水量水质模拟及管理程序集[M]．吉林科学技术出版社，1988
[45] 郭卫星，卢国平．MODFLOW三维有限差分地下水流模型[M]．南京大学地球科学系，1998
[46] 武强，董东林等．水资源评价的可视化专业软件(Visual Modflow)与应用潜力[J]．水文地质工程地质，1999(5)：21-23
[47] 陈玉田．偏微分方程数值解法[M]．河海大学出版社，1999
[48] 陈崇喜，唐仲华．地下水流动问题数值方法[M]．中国地质大学出版社，1990
[49] 陈锁忠．MODFLOW三维有限差分地下水渗流模型的原理[J]．江苏地矿信息，1999(1)：21-23
[50] 钱家忠，吴剑锋等．徐州市张集水源地裂隙岩溶水三维等参有限元数值模拟[J]．水利学报，2003(3)：37—1
[51] 宋西平，寇广潮．地下水数值模拟面参数和水文参数研究的重要性[J]．地下水，2001(2)：63-64
[52] 刘振英，李文鹏，孙鸣．水资源实时模拟管理系统及其在山东临淄地区的应用[M]．地质出版社，1996
[53] 国家技术监督局．中华人民共和国国家标准(GB／T14497-93)地下水资源管理模型工作技术要求[M]．北京，1994
[54] 朱学愚，谢春红．地下水运移模型[M]．中国建筑工业出版社，1990
[55] 朱学愚，钱孝星，刘新仁．地下水资源评价[M]．南京大学出版社，1987
[56] 张蔚棒．地下水非稳定流计算和地下水资源评价[M]．科学出版社，1983
[57] 张宏仁等．地下水非稳定流理论的发展和应用M]．地质出版社，1975
[58] 束龙仓，朱元生，地下水资源评价中的不确定性因素分析[J]，水文地质工程地质，2000，27(6)：6-8
[59] 束龙仓，朱元生，孙庆义，等．地下水资源评价结果的可靠性探讨[J]．水科学进展，2000，11(1)：21-24
[60] 束龙仓，朱元生．晋祠泉域地下水开采决策的风险分析[J]．河海大学学报，2000，28(6)：90-93
[61] 薛禹群，吴吉春．面临21世纪的中国地下水模拟问题[J]．水文地质工程地质，1999，(5)：1-3
[62] 束龙仓．山东省济宁市地下水资源管理模型研究[M]．建设部综合勘察研究院，1990
[63] 钱孝星．水文地质计算[M]．水利电力出版社，1995
[64] 卢文喜．地下水运动数值模拟过程中边界条件问题探讨[J]．水利学报，2003(3)：34-36
[65] 洪景涛，刘翔，严涛．淄博市大武水源地地下水流场数值模型研究[J]．水文地质学报，2003(4)：7-10
[66] 蒋亚萍，陈余道．MODFLOW——套水文地质学实用计算软件[J]．广西地质，1999(3)：75-78
[67] 邵景力，崔亚莉等．黄河下游影响带(河南段)三维地下水流数值模拟模型及其应用[J]．吉林大学学报(地球科学版)，2003，33(1)：55-58