非完整河渠附近地下水运动规律探讨
|
摘要
地下水与河水的相互作用是自然界中普遍存在的一种自然现象，也是陆地水文循环的一个重要组成部分。在自然界中，河水与地下水是紧密联系的整体，二者之间不仅发生水量交换，而且存在溶质和污染物的迁移。随着人类大规模地开采傍河水源地的地下水，河流对地下水系统的影响更为显著，二者之间的相互作用机理也更为复杂。因此，在建立地下水流模型时，正确分析河水与地下水之间的相互关系，对河流进行科学的概化是十分重要的。
针对解析法中对河流边界多按完整河情况处理这一问题，建立了非完整河数值模型，对潜水非完整河附近地下水运动进行了模拟分析，讨论了河流完整度对单宽流量和地下水位的影响，给出了非完整河流单宽流量和水位计算公式，结果表明，河流完整性不同，其附近的地下水动态也随之变化。
在二元结构模型中，考虑上部粘性土层的不同厚度、不同渗透性，经过数值模拟发现，上部弱透水层具有相对隔水作用，且厚度越大，渗透性越小，隔水性越好；当河流切穿上部弱透水层时，研究区地下水位都会发生突变性的抬升，而考虑回水时间的影响，地下水位的变化量比未切穿弱透水层要大。
水库浸没问题是水库工程勘察与评价的重要内容，也是水库区五大工程地质问题之一，本文在对目前应用于水库浸没预测的几种方法进行了对比分析的基础上，基于吉林省长春市石头口门水库兴利增容工程，选取典型断面，考虑河流切割含水层的不同程度、回水时间及潜水蒸发的影响，对增容后库区地下水位壅高情况进行研究，并对库区浸没问题进行了分析。
The interaction between the groundwater and river water is a common natural phenomenon, and it also is the important part of land water cycle. River water and groundwater is a close connection unity, there are not only water exchanges but also solute and contaminant transport exist between them. With the large scale groundwater exploitation in water source site adjacent to river, the influence of the river to groundwater system is more prominent, and the interaction mechanism of them is more complicated. So, correctly analysis the connection of river water and groundwater is important to the scientific simplification of river when building groundwater flow model.
Aiming at the problem that Rivers is always treated as full penetrating one in analytic method, the numerical model of partial penetrating river is constructed, and the groundwater flow near partial penetrating river is simulated. The impact on unit discharge and water head is discussed, The formulas for calculating unit discharge and water head have been obtained, The results show that the groundwater dynamic state varies with the river integrity.
In dualistic structure model, the thickness and permeability of upper layer clay are considered. The numerical simulation shows that the upper semipervious formation acts as relatively impervious layer, and permeability decrease with the increasing or thickness of the upper semipervious formation; when river penetrating the back plane of upper semipervious formation, the groundwater table of research side would elevate abruptly, and when considering the effect of backwater time the elevation of groundwater table is bigger than that of river doesn't penetrate the upper semipervious formation
Reservoir immersion is an important content of reservoir investigation and appraisement, and it also is one of five engineering geology problems of reservoir engineering. In this paper, based on the comparison of Reservoir immersion prediction method, and to the Shitoukoumen Reservoir storage increment engineering in Jiling province Changchun city, the typical section is selected; then considering the extent of the river penetrating the aquifer and effect of backwater time, the banked-up water levels of groundwater is researched after Reservoir storage increment, and reservoir immersion is predicted.
针对解析法中对河流边界多按完整河情况处理这一问题，建立了非完整河数值模型，对潜水非完整河附近地下水运动进行了模拟分析，讨论了河流完整度对单宽流量和地下水位的影响，给出了非完整河流单宽流量和水位计算公式，结果表明，河流完整性不同，其附近的地下水动态也随之变化。
在二元结构模型中，考虑上部粘性土层的不同厚度、不同渗透性，经过数值模拟发现，上部弱透水层具有相对隔水作用，且厚度越大，渗透性越小，隔水性越好；当河流切穿上部弱透水层时，研究区地下水位都会发生突变性的抬升，而考虑回水时间的影响，地下水位的变化量比未切穿弱透水层要大。
水库浸没问题是水库工程勘察与评价的重要内容，也是水库区五大工程地质问题之一，本文在对目前应用于水库浸没预测的几种方法进行了对比分析的基础上，基于吉林省长春市石头口门水库兴利增容工程，选取典型断面，考虑河流切割含水层的不同程度、回水时间及潜水蒸发的影响，对增容后库区地下水位壅高情况进行研究，并对库区浸没问题进行了分析。
The interaction between the groundwater and river water is a common natural phenomenon, and it also is the important part of land water cycle. River water and groundwater is a close connection unity, there are not only water exchanges but also solute and contaminant transport exist between them. With the large scale groundwater exploitation in water source site adjacent to river, the influence of the river to groundwater system is more prominent, and the interaction mechanism of them is more complicated. So, correctly analysis the connection of river water and groundwater is important to the scientific simplification of river when building groundwater flow model.
Aiming at the problem that Rivers is always treated as full penetrating one in analytic method, the numerical model of partial penetrating river is constructed, and the groundwater flow near partial penetrating river is simulated. The impact on unit discharge and water head is discussed, The formulas for calculating unit discharge and water head have been obtained, The results show that the groundwater dynamic state varies with the river integrity.
In dualistic structure model, the thickness and permeability of upper layer clay are considered. The numerical simulation shows that the upper semipervious formation acts as relatively impervious layer, and permeability decrease with the increasing or thickness of the upper semipervious formation; when river penetrating the back plane of upper semipervious formation, the groundwater table of research side would elevate abruptly, and when considering the effect of backwater time the elevation of groundwater table is bigger than that of river doesn't penetrate the upper semipervious formation
Reservoir immersion is an important content of reservoir investigation and appraisement, and it also is one of five engineering geology problems of reservoir engineering. In this paper, based on the comparison of Reservoir immersion prediction method, and to the Shitoukoumen Reservoir storage increment engineering in Jiling province Changchun city, the typical section is selected; then considering the extent of the river penetrating the aquifer and effect of backwater time, the banked-up water levels of groundwater is researched after Reservoir storage increment, and reservoir immersion is predicted.
引文
[1] 沈振荣．瑜芳等．水资源科学实验与研究[M]．北京：中国科学技术出版社，1992，415～418．
[2] http://www.waterinfo.com.cn/zhongdian/guonei/200302240029.htm.
[3] http://www.cas.ac.cn/index/0C/41/index3.htm.
[4] 胡俊峰．王金生等．地下水与河水相互作用的研究进展[J]．水文地质工程地质，2004，(1)：108～112．
[5] 谭杰．王杰等．浅谈地下水流模拟中河流的处理方法[J]．地下水，2006，28(4)：28～30
[6] 薛禹群主编．地下水动力学[M]．北京：地质出版社，1986．
[7] Silliman, S E. and D. F. Booth. Analysis of time-series measurements of sediment temperature for identification of gaining vs. losing portions of Juday Creek[J]. J. Hydrol., 1993, 146(1-4): 131～148.
[8] Belanger, T V. and R B. Walker. Groundwater seepage in the Indian River Lagoon [J]. Tropical Hydrology and Caribbean Water Resource, 1990, 367～375.
[9] 仵彦卿．张应华等．西北黑河下游盆地河水与地下水转化的新发现[J]．自然科学进展，2004(12)：1428～1433．
[10] 郭东屏．傍河三维与二维井流的电模拟对比试验研究[J]．西安地质学院学报，1993，15(增刊)：12～16．
[11] 吴耀国．王超等．河流-地下水渗流系统对城市污水的净化作用及机理[J]．环境科学学报，2002，22(1)：60～64．
[12] Payne B R. The status of isotope hydrology today[J]. J. Hydrol., 1988, 100(1-3): 207～237.
[13] Mathieu R T, Bariac. An isotopic study (~2H and ~(18)O) on water movements in clayey soils under a semiarid climate[J]. Water Resource Research, 1996, 32: 779～789.
[14] 顾慰祖．陆家驹等乌兰布和沙漠北部地下水资源的环境同位素探讨[J]．水科学进展，2002，19(3)：326～332．
[15] 陈建生．赵霞等．巴丹吉林沙漠湖泊钙华与根状结核的发现对研究湖泊水补给的意义[J]．中国岩溶，2004，23(4)：277～282．
[16] Heij, G J. River-groundwater relationships in the lower parts of The Netherlands [J].J.Hydrol, 1989, 108(1-4): 35～62.
[17] Spalding, C P. and R. Khaleel, An evaluation of analytical solutions to estimate drawdowns and stream depletions by wells[J]. Water Resources Research, 1991, 27(4): 597～609.
[18] Rufledge A T. Methods of using streamflow records for estimating total and effective recharge in the Appalachian Valley and Ridge, Piedmont, and Blue Ridge physiographic provinces[J]. Water Resources, 1992, 17: 59～73.
[19] Theis C V. The effect of a well on the flow of a nearby stream[J]. Geophys, 1941, 22(3): 734～738.
[20] Hantush M S. Wells near streams with semi-pervious beds[J]. J Geophys. Research. 1965, 70(12): 2829～2838.
[21] Meigs L C. and J MBahr. Three-dimensional groundwater flow near narrow surface water bodies[J]. Water Resource Researh. 1995, 31(12): 3299～3307.
[22] Hunt B. Unsteady stream depletion from ground water pumping[J]. Ground Water, 1999, 37(1): 98～102.
[23] Christensen S. On the estimation of streamflow depletion parameters by drawdown analysis[J]. Ground Water, 2000, 38(5): 726～734.
[24] Vitaly A Zlotnik, Huihua Huang. Effect of shallow penetration and streambed sediments on aquifer response to stream stage fluctuations (analytical model)[J]. Ground Water, 1999, 37(4): 599～605.
[25] James J, Butler, Vitaly A Z, et al. Drawdown and stream depletion produced by pumping in the vicinity of a partially penetrating stream[J]. Ground Water, 2001, 5: 651～659.
[26] Sophocleous M A, A Koussis, J L Martin et al. Evaluation of simplified stream2aquifer depletion models for water rights administration[J]. Ground Water, 1995, 33(4): 579～588.
[27] Chen X H , Y Yin. Streamflow depletion: Modeling of base flow reduction and stream infiltration from seasonally pumped wells[J]. Journal of American Water Resources Association. 2001, 37(1): 185～195.
[28] 钱会．郑西来等．傍河取水越河稳定渗流问题的三维数值模拟研究[J]．水利学报，1999，3：32～37．
[29] 薛禹群．吴吉春．地下水数值模拟在我国-回顾与展望[J]．水文地质工程地质，1997，24(4)：21～24．
[30] 李致家．谢悦波．地下水流与河网水流的耦合模型[J]．水利学报，1998，4：43～47．
[31] 蒋业放．张兴有．河流与含水层水力耦合模型及其应用[J]．地理学报，1999，54(6)：526～533．
[32] 崔亚莉．邵景力等．玛纳斯河流域地表水与地下水相互转化关系研究[J]．水文地质工程地质，2001，28(2)：9～13．
[33] 陈崇希．裴顺平等．非完整河的数值模拟方法及建模中的若干问题——读“数值模拟方法在评价地下水资源时区内河流的处理方法”一文随笔[J]．勘察科学技术，1999，(4)：3～12
[34] 赵坚．纪伟．山区型水库有限元计算模型边界水位的修正方法研究[J]．水电能源科学，2006，24(1)：40～42
[35] 易树平．迟宝明．傍河地下水源地数值模拟与评价研究——以沈阳市浑河李官堡水源地为例[J]．自然资源学报，2006，21(1)：154～159．
[36] 刘猛，傍河水源地地下水数值模拟研究[D]．南京：河海大学硕士学位论文，2006．
[37] 张志忠．武强．河水与地下水耦合模型的建立与应用[J]．辽宁工程技术大学学报，2004，23(4)：449～452．
[38] 李粉蝉．张永波等．周期性摆动河流傍河取水地下水数值模拟研究[J]．华北电力大学学报，2005，32(3)：73～76．
[39] 杜文堂．组合模型在河流补给量预测中的应用[J]．河北建筑科技学院学报，2000，17(3)：68～71．
[40] 王君连．工程地下水计算[M]．北京：中国水利水电出版社，2004．
[41] 陈崇希．地下水动力学[M]．武汉：中国地质大学出版社．1999．
[42] 张蔚椿．地下水非稳定流计算和地下水资源评价[M]，北京：科学出版社，1983．
[43] van Schilfgaard J. Theory of flow to drains[J]. Advanced in Hydroscience, 6, 43～106.
[44] Polubarinova-Kochina P Ya. Theory of Ground Water Movement[M]. Princeton University Press: NJ, 1962.
[45] Li W H. Differential equations of hydraulic transient, dispersion, and groundwater flow[M]. prentice-Hall, Englrwood Cliffs NJ, 1972.
[46] Sidiropoulos E G, Tolikas P K. Similarity and iterative solution of the Boussinesq equation[J]. Journal of Hydrology, 1974, 74, 31～41.
[47] 陈静．河口海岸带附近地下水运移规律研究[D]．南京．河海大学博士学位论文，2006．
[48] 阿里木·吐尔逊．周志芳等．河渠附近潜水非稳定运动的一种通解[J]．河海大学学报(自然科学版)，2003，31(6)：649～651．
[49] 陶月赞．曹彭强等．垂向入渗与河渠边界影响下潜水非稳定流参数的求解[J]．水利学报，2006，37(8)：913～917．
[50] 孙讷正．地下水流的数学模型和数值方法[M]．北京：地质出版社，1981．
[51] 罗焕炎．陈雨孙著．地下水运动的数值模拟[M]．北京：中国建材工业出版社．1988．
[52] 钱孝星主编．水文地质计算[M]．北京：水利电力出版社．1995．
[53] 兰艇雁．王树林等．册田水库的浸没预测与预防[J]．山西水利科技(增刊)，1999，6：56～58．
[54] 李帆．车林源．潮州供水枢纽工程浸没初测及防治对策[J]．西部探矿工程，2005，10：225～227．
[55] 曹洪．尹小玲．潮州供水枢纽工程库区浸没影响初探[J]．岩石力学与工程学报，2004，23(5)：862～866．
[56] 朱杰，王汇明．二维渗流有限元法在潮州水库浸没分析中的应用[J]．四川地质学报，2006，26(2)：121～124．
[57] 冀建疆．官厅水库的浸没评价和范围预测[J]．水利水电技术，2005，36(2)：18～21．
[58] 张书光，杨彦平等．黄河西霞院水库浸没问题分析[J]．华北水利水电学院学报，2002，23(4)：20～21．
[59] 张会刚．夏克勤等．四川某水电站库区浸没问题预测[J]．水土保持研究，2004，11(3)：140～142．
[60] 鲍立新．佟胤铮．阎王鼻子水库浸没问题的分析与评价[J]．东北水利水电，2002，21(3)：37～56．
[61] 谢先红．水库浸没及堤内排水降压效果研究[D]，武汉：武汉大学，2005．
[62] 袁宏利．杨计申等．水库浸没勘察研究工作的新思路[J]．水利水电工程设计，2003，22(4)：42～44．
[63] 湖南省水利水电勘测设计研究总院，长沙综合枢纽工程可行性研究阶段库区浸没问题专题研究汇报材料，2006．3．
[64] 刘起霞．李清波．邹剑峰．环境工程地质[M]．黄河水利出版社，2000．
[65] 陆兆溱主编．工程地质学[M]．水利水电出版社，1986．
[66] 吉林省水利水电勘测设计总院，吉林省长春市石头口门水库兴利增容工程初步设计，2001．1．
[2] http://www.waterinfo.com.cn/zhongdian/guonei/200302240029.htm.
[3] http://www.cas.ac.cn/index/0C/41/index3.htm.
[4] 胡俊峰．王金生等．地下水与河水相互作用的研究进展[J]．水文地质工程地质，2004，(1)：108～112．
[5] 谭杰．王杰等．浅谈地下水流模拟中河流的处理方法[J]．地下水，2006，28(4)：28～30
[6] 薛禹群主编．地下水动力学[M]．北京：地质出版社，1986．
[7] Silliman, S E. and D. F. Booth. Analysis of time-series measurements of sediment temperature for identification of gaining vs. losing portions of Juday Creek[J]. J. Hydrol., 1993, 146(1-4): 131～148.
[8] Belanger, T V. and R B. Walker. Groundwater seepage in the Indian River Lagoon [J]. Tropical Hydrology and Caribbean Water Resource, 1990, 367～375.
[9] 仵彦卿．张应华等．西北黑河下游盆地河水与地下水转化的新发现[J]．自然科学进展，2004(12)：1428～1433．
[10] 郭东屏．傍河三维与二维井流的电模拟对比试验研究[J]．西安地质学院学报，1993，15(增刊)：12～16．
[11] 吴耀国．王超等．河流-地下水渗流系统对城市污水的净化作用及机理[J]．环境科学学报，2002，22(1)：60～64．
[12] Payne B R. The status of isotope hydrology today[J]. J. Hydrol., 1988, 100(1-3): 207～237.
[13] Mathieu R T, Bariac. An isotopic study (~2H and ~(18)O) on water movements in clayey soils under a semiarid climate[J]. Water Resource Research, 1996, 32: 779～789.
[14] 顾慰祖．陆家驹等乌兰布和沙漠北部地下水资源的环境同位素探讨[J]．水科学进展，2002，19(3)：326～332．
[15] 陈建生．赵霞等．巴丹吉林沙漠湖泊钙华与根状结核的发现对研究湖泊水补给的意义[J]．中国岩溶，2004，23(4)：277～282．
[16] Heij, G J. River-groundwater relationships in the lower parts of The Netherlands [J].J.Hydrol, 1989, 108(1-4): 35～62.
[17] Spalding, C P. and R. Khaleel, An evaluation of analytical solutions to estimate drawdowns and stream depletions by wells[J]. Water Resources Research, 1991, 27(4): 597～609.
[18] Rufledge A T. Methods of using streamflow records for estimating total and effective recharge in the Appalachian Valley and Ridge, Piedmont, and Blue Ridge physiographic provinces[J]. Water Resources, 1992, 17: 59～73.
[19] Theis C V. The effect of a well on the flow of a nearby stream[J]. Geophys, 1941, 22(3): 734～738.
[20] Hantush M S. Wells near streams with semi-pervious beds[J]. J Geophys. Research. 1965, 70(12): 2829～2838.
[21] Meigs L C. and J MBahr. Three-dimensional groundwater flow near narrow surface water bodies[J]. Water Resource Researh. 1995, 31(12): 3299～3307.
[22] Hunt B. Unsteady stream depletion from ground water pumping[J]. Ground Water, 1999, 37(1): 98～102.
[23] Christensen S. On the estimation of streamflow depletion parameters by drawdown analysis[J]. Ground Water, 2000, 38(5): 726～734.
[24] Vitaly A Zlotnik, Huihua Huang. Effect of shallow penetration and streambed sediments on aquifer response to stream stage fluctuations (analytical model)[J]. Ground Water, 1999, 37(4): 599～605.
[25] James J, Butler, Vitaly A Z, et al. Drawdown and stream depletion produced by pumping in the vicinity of a partially penetrating stream[J]. Ground Water, 2001, 5: 651～659.
[26] Sophocleous M A, A Koussis, J L Martin et al. Evaluation of simplified stream2aquifer depletion models for water rights administration[J]. Ground Water, 1995, 33(4): 579～588.
[27] Chen X H , Y Yin. Streamflow depletion: Modeling of base flow reduction and stream infiltration from seasonally pumped wells[J]. Journal of American Water Resources Association. 2001, 37(1): 185～195.
[28] 钱会．郑西来等．傍河取水越河稳定渗流问题的三维数值模拟研究[J]．水利学报，1999，3：32～37．
[29] 薛禹群．吴吉春．地下水数值模拟在我国-回顾与展望[J]．水文地质工程地质，1997，24(4)：21～24．
[30] 李致家．谢悦波．地下水流与河网水流的耦合模型[J]．水利学报，1998，4：43～47．
[31] 蒋业放．张兴有．河流与含水层水力耦合模型及其应用[J]．地理学报，1999，54(6)：526～533．
[32] 崔亚莉．邵景力等．玛纳斯河流域地表水与地下水相互转化关系研究[J]．水文地质工程地质，2001，28(2)：9～13．
[33] 陈崇希．裴顺平等．非完整河的数值模拟方法及建模中的若干问题——读“数值模拟方法在评价地下水资源时区内河流的处理方法”一文随笔[J]．勘察科学技术，1999，(4)：3～12
[34] 赵坚．纪伟．山区型水库有限元计算模型边界水位的修正方法研究[J]．水电能源科学，2006，24(1)：40～42
[35] 易树平．迟宝明．傍河地下水源地数值模拟与评价研究——以沈阳市浑河李官堡水源地为例[J]．自然资源学报，2006，21(1)：154～159．
[36] 刘猛，傍河水源地地下水数值模拟研究[D]．南京：河海大学硕士学位论文，2006．
[37] 张志忠．武强．河水与地下水耦合模型的建立与应用[J]．辽宁工程技术大学学报，2004，23(4)：449～452．
[38] 李粉蝉．张永波等．周期性摆动河流傍河取水地下水数值模拟研究[J]．华北电力大学学报，2005，32(3)：73～76．
[39] 杜文堂．组合模型在河流补给量预测中的应用[J]．河北建筑科技学院学报，2000，17(3)：68～71．
[40] 王君连．工程地下水计算[M]．北京：中国水利水电出版社，2004．
[41] 陈崇希．地下水动力学[M]．武汉：中国地质大学出版社．1999．
[42] 张蔚椿．地下水非稳定流计算和地下水资源评价[M]，北京：科学出版社，1983．
[43] van Schilfgaard J. Theory of flow to drains[J]. Advanced in Hydroscience, 6, 43～106.
[44] Polubarinova-Kochina P Ya. Theory of Ground Water Movement[M]. Princeton University Press: NJ, 1962.
[45] Li W H. Differential equations of hydraulic transient, dispersion, and groundwater flow[M]. prentice-Hall, Englrwood Cliffs NJ, 1972.
[46] Sidiropoulos E G, Tolikas P K. Similarity and iterative solution of the Boussinesq equation[J]. Journal of Hydrology, 1974, 74, 31～41.
[47] 陈静．河口海岸带附近地下水运移规律研究[D]．南京．河海大学博士学位论文，2006．
[48] 阿里木·吐尔逊．周志芳等．河渠附近潜水非稳定运动的一种通解[J]．河海大学学报(自然科学版)，2003，31(6)：649～651．
[49] 陶月赞．曹彭强等．垂向入渗与河渠边界影响下潜水非稳定流参数的求解[J]．水利学报，2006，37(8)：913～917．
[50] 孙讷正．地下水流的数学模型和数值方法[M]．北京：地质出版社，1981．
[51] 罗焕炎．陈雨孙著．地下水运动的数值模拟[M]．北京：中国建材工业出版社．1988．
[52] 钱孝星主编．水文地质计算[M]．北京：水利电力出版社．1995．
[53] 兰艇雁．王树林等．册田水库的浸没预测与预防[J]．山西水利科技(增刊)，1999，6：56～58．
[54] 李帆．车林源．潮州供水枢纽工程浸没初测及防治对策[J]．西部探矿工程，2005，10：225～227．
[55] 曹洪．尹小玲．潮州供水枢纽工程库区浸没影响初探[J]．岩石力学与工程学报，2004，23(5)：862～866．
[56] 朱杰，王汇明．二维渗流有限元法在潮州水库浸没分析中的应用[J]．四川地质学报，2006，26(2)：121～124．
[57] 冀建疆．官厅水库的浸没评价和范围预测[J]．水利水电技术，2005，36(2)：18～21．
[58] 张书光，杨彦平等．黄河西霞院水库浸没问题分析[J]．华北水利水电学院学报，2002，23(4)：20～21．
[59] 张会刚．夏克勤等．四川某水电站库区浸没问题预测[J]．水土保持研究，2004，11(3)：140～142．
[60] 鲍立新．佟胤铮．阎王鼻子水库浸没问题的分析与评价[J]．东北水利水电，2002，21(3)：37～56．
[61] 谢先红．水库浸没及堤内排水降压效果研究[D]，武汉：武汉大学，2005．
[62] 袁宏利．杨计申等．水库浸没勘察研究工作的新思路[J]．水利水电工程设计，2003，22(4)：42～44．
[63] 湖南省水利水电勘测设计研究总院，长沙综合枢纽工程可行性研究阶段库区浸没问题专题研究汇报材料，2006．3．
[64] 刘起霞．李清波．邹剑峰．环境工程地质[M]．黄河水利出版社，2000．
[65] 陆兆溱主编．工程地质学[M]．水利水电出版社，1986．
[66] 吉林省水利水电勘测设计总院，吉林省长春市石头口门水库兴利增容工程初步设计，2001．1．