大连甘井子区地下水数值模拟研究
|
摘要
地下水是大连市生活和生产用水的重要补充水源。随着大连经济的快速发展,对水资源的需求也日益增长,部分地区超量开采地下水,产生地下水位下降、海水入侵等后果。对大连地区的地下水开采进行管理已刻不容缓。因此,需要对该区的地下水进行模拟研究,对地下水资源进行科学评价。分布式地下水数值模拟模型是研究地下水动态的有效工具。
本文针对大连地下水资源综合利用中存在的问题,重点开展了以下研究：
(1)选择甘井子区地下水系统为重点研究对象,在深入分析研究区地质、水文地质条件的基础上,建立了大连甘井子区的二维、各向同性、均质的非稳定流数值模拟模型,根据2006～2007年的地下水补给和排泄动态资料,采用稳定流法确定模型的初始流场,根据岩性特点及钻孔资料确定研究区边界及参数分区。根据2008～2009年的补给和排泄资料,采用试估-校正法,利用Visual Modflow的Modflow2000模块对模型进行识别。然后利用2010年的降雨和开采资料对模型进行验证。
(2)利用1990～2004年的降雨系列,预测在现状开采条件下该区2011～2025年地下水流场变化,并进行水量均衡分析、资源评价。
(3)依据资源评价的结果和研究区用水户的用水规划,设定三种开采情景,预测2011～2025年地下水流场变化。
通过对三种开采情景下研究区的降落漏斗、水量均衡进行分析,得出以下结论：
(1)2011～2025年共计15年研究区地下水系统总的补给量为31380.05×104m3,补给量来源全部为降雨入渗和侧向流入,未发生海水倒灌；总排泄量为29347.93×104m3。研究区的天然补给量约为2092×104m3/a,补给模数约为10.55×104m3/(km2·a)。在现状开采条件下,地下水储量虽有所增加,但是漏斗依然存在。
(2)石灰石矿区停止开采地下水,有利于消除地下水降落漏斗,其他井增大开采量有利于袭夺部分蒸发量和入海水量,可以增加地下水的利用率。
(3)通过对比分析结果可知,矿区停止开采,增大远离海岸开采井的开采量,即情景三最利于地下水资源的可持续利用。
Groundwater, one of the important supplementary water sources, was over-pumped as more water quantity was demanded with the rapid development of the economies in some areas in Dalian. It resulted in the decline of groundwater table, seawater intrusion and so on. There's no time to delay to manage the groundwater pumping in this city. For this reason, it needs to simulate the groundwater and evaluate the groundwater resource in this area. Distributed groundwater numerical simulation model is efficient to research on dynamic of groundwater.
Aiming at problems when multipurpose utilization of groundwater in Dalian, the main studies were drawn as follows:
(1) Groundwater system in Ganjingzi district was chosen to be a key study object. Based on the analysis of hydro-geological condition about this study area, two-dimensional homogeneous isotropic unsteady flow numerical simulating model was constructed. The initial flow field was computed by steady flow method using the data of average precipitation and exploitation quantity from2006to2007. According to lithological characters and borehole data, boundaries of study area and parameter sub areas were assigned. With the trial-and-error method, parameters of the model were identified by Modflow2000module which is embedded in Visual Modflow software. The identification period was from2008to2009. After the parameters were identified, the model was tested by data of precipitation and exploitation in2010.
(2) According to the data series of precipitation from1990to2004and the current groundwater developing conditions, the groundwater flow fields of study area during2011and2025have been predicted.
(3) According to water assessment results (2011a-2025a) and water using plan,3scenes were made, and the groundwater flow field changes were predicted.
With the analysis on cones of depression and water balance budgets of3scenes, the conclusions were drawn as follows:
(1) In simulated period(2011a-2025a), total recharge31380.05×104m3, total discharge is29347.93×104m3. Natural recharge is about2092×104m3, and recharge modulus is about10.55×104m3/(km2·a). In the condition of the current groundwater exploitation, the storage of groundwater is increasing, but the cone of depression can't disappear.
(2)Stopping exploitation in limestone mine is good for recovery of groundwater. Increasing pumping rates of else wells can reduce part of evaporation and discharge fluxes to ocean. Groundwater resource can be made full use by these measures.
(3) Comparing these results, Scene3which includes stopping exploitation in limestone mine and increasing pumping rates of else wells which are far away from sea is more better for sustainable utilization of groundwater resource than other scenes.
本文针对大连地下水资源综合利用中存在的问题,重点开展了以下研究：
(1)选择甘井子区地下水系统为重点研究对象,在深入分析研究区地质、水文地质条件的基础上,建立了大连甘井子区的二维、各向同性、均质的非稳定流数值模拟模型,根据2006～2007年的地下水补给和排泄动态资料,采用稳定流法确定模型的初始流场,根据岩性特点及钻孔资料确定研究区边界及参数分区。根据2008～2009年的补给和排泄资料,采用试估-校正法,利用Visual Modflow的Modflow2000模块对模型进行识别。然后利用2010年的降雨和开采资料对模型进行验证。
(2)利用1990～2004年的降雨系列,预测在现状开采条件下该区2011～2025年地下水流场变化,并进行水量均衡分析、资源评价。
(3)依据资源评价的结果和研究区用水户的用水规划,设定三种开采情景,预测2011～2025年地下水流场变化。
通过对三种开采情景下研究区的降落漏斗、水量均衡进行分析,得出以下结论：
(1)2011～2025年共计15年研究区地下水系统总的补给量为31380.05×104m3,补给量来源全部为降雨入渗和侧向流入,未发生海水倒灌；总排泄量为29347.93×104m3。研究区的天然补给量约为2092×104m3/a,补给模数约为10.55×104m3/(km2·a)。在现状开采条件下,地下水储量虽有所增加,但是漏斗依然存在。
(2)石灰石矿区停止开采地下水,有利于消除地下水降落漏斗,其他井增大开采量有利于袭夺部分蒸发量和入海水量,可以增加地下水的利用率。
(3)通过对比分析结果可知,矿区停止开采,增大远离海岸开采井的开采量,即情景三最利于地下水资源的可持续利用。
Groundwater, one of the important supplementary water sources, was over-pumped as more water quantity was demanded with the rapid development of the economies in some areas in Dalian. It resulted in the decline of groundwater table, seawater intrusion and so on. There's no time to delay to manage the groundwater pumping in this city. For this reason, it needs to simulate the groundwater and evaluate the groundwater resource in this area. Distributed groundwater numerical simulation model is efficient to research on dynamic of groundwater.
Aiming at problems when multipurpose utilization of groundwater in Dalian, the main studies were drawn as follows:
(1) Groundwater system in Ganjingzi district was chosen to be a key study object. Based on the analysis of hydro-geological condition about this study area, two-dimensional homogeneous isotropic unsteady flow numerical simulating model was constructed. The initial flow field was computed by steady flow method using the data of average precipitation and exploitation quantity from2006to2007. According to lithological characters and borehole data, boundaries of study area and parameter sub areas were assigned. With the trial-and-error method, parameters of the model were identified by Modflow2000module which is embedded in Visual Modflow software. The identification period was from2008to2009. After the parameters were identified, the model was tested by data of precipitation and exploitation in2010.
(2) According to the data series of precipitation from1990to2004and the current groundwater developing conditions, the groundwater flow fields of study area during2011and2025have been predicted.
(3) According to water assessment results (2011a-2025a) and water using plan,3scenes were made, and the groundwater flow field changes were predicted.
With the analysis on cones of depression and water balance budgets of3scenes, the conclusions were drawn as follows:
(1) In simulated period(2011a-2025a), total recharge31380.05×104m3, total discharge is29347.93×104m3. Natural recharge is about2092×104m3, and recharge modulus is about10.55×104m3/(km2·a). In the condition of the current groundwater exploitation, the storage of groundwater is increasing, but the cone of depression can't disappear.
(2)Stopping exploitation in limestone mine is good for recovery of groundwater. Increasing pumping rates of else wells can reduce part of evaporation and discharge fluxes to ocean. Groundwater resource can be made full use by these measures.
(3) Comparing these results, Scene3which includes stopping exploitation in limestone mine and increasing pumping rates of else wells which are far away from sea is more better for sustainable utilization of groundwater resource than other scenes.
引文
[1]乔玲.北京地区城市供水保障系统战略研究-平谷地下水应急水源地开采评价[D].长春：吉林大学建设工程学院,2005.
[2]施德鸿.城市地下水资源的开发利用与生存环境[C].第三届全国减轻自然灾害研讨会.1998,5：344-348.
[3]高艳玲.城市水务管理[M].北京：中国建材工业出版社.2006：1-5.
[4]陈梦熊,马凤山.中国地下水资源与环境[M].北京：地震出版社.2002：470-476.
[5]中国地下水科学战略研究小组.中国地下水科学的机遇与挑战[M].北京：科学出版社.2009.
[6]毕永传.日照市滨海区浅层地下水系统调控研究[D].泰安：山东农业大学水利土木工程学院,2010.
[7]郑铣鑫,武强,侯艳声等.城市地面沉降研究进展及其发展趋势[J].地质论评,2002,48(6)：612-618.
[8]郭巧娜.滨海多层含水层系统中多组分地下水流及海底地下水排泄[D].中国地质大学(武汉)环境学院,2010.
[9]井锋.大连市水资源系统恢复力评价研究[D].大连：辽宁师范大学,2010.
[10]马秀媛.地下水预测的新方法研究[D].济南：山东大学,2007.
[11]林锦.变密度条件下地下水模拟优化研究[D].杭州：浙江大学理学院,2008.
[12]朱奎.大区域地下水数值模拟——以华北滏阳河平原为例.[D]武汉：武汉大学,2004.
[13]薛禹群.地下水动力学[M].北京：地质出版社,2005：28-45.
[14]谭雪清.城镇化对地下水影响数值模拟[D].南昌大学建筑工程学院,2010.
[15]吴媛.若尔盖湿地保护区地下水流数值模拟[D].成都：四川大学,2006.
[16]魏连伟.基于人工智能技术的地下水系统参数识别研究[D].天津：天津大学,2003.
[17]蒙媛.丰台区地下水模拟研究[D].合肥：合肥工业大学,2007.
[18]王浩,陆垂裕,秦大庸等.地下水数值计算与应用研究进展综述[J].地学前缘(中国地质大学(北京)),2010,17(6)：1-12.
[19]冯康.基于变分原理的差分格式[J].应用数学与计算数学,1965(2)：238-262.
[20]薛禹群,谢春红.水文地质学的数值法[M].北京：煤炭工业出版社,1979.
[21]王文科,李俊亭地下水流数值模拟的发展与展望[J].西北地质,1995,16(4)：52-56
[22]Pinder G. F., J. D. Bredehoeft. Application of the digital computer for aquifer evaluation. Water Resources,1968,4(5):1069-1093.
[23]孙纳正.地下水流的数学模型和数值方法[M].北京：地质出版社,1981.
[24]Courant R. Variational methods for a solution of problems of equilibrium and vibnltion. [J]Bull:Amer. Math. sec.1943,49:1-23
[25]林坜.大区域地下水流模拟研究及FEFLOW的建模方法-以华北平原为例[D].中国地质大 学,2006.
[26]刘晓民.干旱区典型地下水系统的动态模拟与模拟优化耦合管理模型研究[D].内蒙古农业大学,2009.
[27]刘玉珍.平原区地下水系统模拟与决策支持研究[D].大连理工大学,2008.
[28]冯绍元.解地下水中溶质运移问题的特征有限元法[J].武汉水利电力学院学报,1991,24(4)：410-417.
[29]谢春红,张勤,丁家平.稳定渗流的有限元计算新方法[J].水利学报,1995(12)：29-38.
[30]孙慧,周德亮.用无网格伽辽金法求解稳定地下水流问题[J].地下水,2008,30(6)：1-3.
[31]张宏仁,李俊亭.有限差分法与有限单元法在渗流问题中的对比[J].水文地质工程地质,1979(2)：50-55.
[32]Kitahara M. Boundary integral equation methods in eigenvalue problems of elastic dynamics and thin Plates[M]. New York:Elsevier Science Publishing Company Inc,1985.
[33]杨天行,刘金英,王运国.求解对流占优地下水水质模型的有限单元高精度广义迎风单元均衡格式及应用[J].长春地质学院,1996(2)：196-201.
[34]张祥伟,竹内邦良.大区域地下水模拟的理论和方法[J].水利学报,2004(6)：7-13.
[35]Mazzia A, Puttii M. Mixed f init e element and f inite volume discretization for heavy brine simulations in groundwater [J]. Journal of Computational and Applied Mathematics, 2002,147:191-213.
[36]Li S G, Mclanghlin D. A computationally practical method for stochastic groundwater modeling [J]. Advances in Water Resources,2003,26:1137-1148.
[37]王浩然,朱国荣,赵金熙.基于区域分解法的地下水有限元与边界元耦合模型[J].地质论评,2003,49(1)：48-52.
[38]魏胜利,苏伟杰.分布式地下水模拟模型MODFLOW介绍与应用[J].水利规划与设计,2011(3)：30-32.
[39]张红霞,宋文.地下水数值模拟的研究现状与展望[J].水利科技与经济,2007,13(11)：794-796.
[40]徐晓民,郭中小,贾利民等.地下水系统数值模拟的应用于展望[J].黑龙江水专学报,2010,37(3)：118-120.
[41]杨春玲,刑世录,韩爱中.地下水系统数值模拟的研究进展[J].内蒙古科技与经济,2007,21(151)：305-307.
[42]丁继红,周德亮,马生忠.国外地下水模拟软件的发展现状与趋势[J].勘察科学技术,2002(1)：37-42.
[43]王庆永,贾忠华,刘晓峰Visual MODFLOW及其在地下水模拟中的应用[J].水资源与水工程学报,2007,18(5)：90-92.
[44]Anon. SSG Software [Z]. Washington:The Scientific Software Group,2000.4-5.
[45]徐乐昌.地下水模拟常用软件介绍[J].铀矿冶,2002,21(1)：33-38.
[46]大连市统计局.2000-2010年大连市国民经济与社会发展统计公报[EB/OL]. 2011-03-01. http://www.stats.dl.gov.cn/class.jsp?dtype=80
[47]大连市城市水井管理所.大连市区及甘井子部分地段地下水勘察报告[R].大连：[出版者不详].1991.
[48]丁玲.大连市周水子地区海水入侵问题研究[D].大连：大连理工大学,2004.
[49]大连市水务局.大连市水资源公报[EB/OL]. 2011-05-12. http://www.swj.dl.gov.cn/html/szygb/index.html
[50]武强.大连大魏家岩溶地下水系统管理模型的研究[J].勘察科学技术.1990.1：23-27.
[51]王国利,周惠成.大连市地下水资源优化管理模型[J].大连理工大学学报,2001.141(1)：112-115.
[52]迟宝明,周彦章,易树平.大连地区沿海小流域地下水库开发潜力分析[J].工程勘察2006(5)：23-26.
[53]马驰.SIP和PCG2两种迭代法在地下水数值计算中的应用对比[J].西安科技学院学报2002.3 22卷1期.
[54]汪晓平,吴勇强,张宏林,等.ASP网络开发技术[M].北京：人民邮电出版社,2000.4-5;189-190.
[55]杨金忠,蔡树英,王旭升.地下水运动数学模型[M].北京：科学出版社,2009.
[2]施德鸿.城市地下水资源的开发利用与生存环境[C].第三届全国减轻自然灾害研讨会.1998,5：344-348.
[3]高艳玲.城市水务管理[M].北京：中国建材工业出版社.2006：1-5.
[4]陈梦熊,马凤山.中国地下水资源与环境[M].北京：地震出版社.2002：470-476.
[5]中国地下水科学战略研究小组.中国地下水科学的机遇与挑战[M].北京：科学出版社.2009.
[6]毕永传.日照市滨海区浅层地下水系统调控研究[D].泰安：山东农业大学水利土木工程学院,2010.
[7]郑铣鑫,武强,侯艳声等.城市地面沉降研究进展及其发展趋势[J].地质论评,2002,48(6)：612-618.
[8]郭巧娜.滨海多层含水层系统中多组分地下水流及海底地下水排泄[D].中国地质大学(武汉)环境学院,2010.
[9]井锋.大连市水资源系统恢复力评价研究[D].大连：辽宁师范大学,2010.
[10]马秀媛.地下水预测的新方法研究[D].济南：山东大学,2007.
[11]林锦.变密度条件下地下水模拟优化研究[D].杭州：浙江大学理学院,2008.
[12]朱奎.大区域地下水数值模拟——以华北滏阳河平原为例.[D]武汉：武汉大学,2004.
[13]薛禹群.地下水动力学[M].北京：地质出版社,2005：28-45.
[14]谭雪清.城镇化对地下水影响数值模拟[D].南昌大学建筑工程学院,2010.
[15]吴媛.若尔盖湿地保护区地下水流数值模拟[D].成都：四川大学,2006.
[16]魏连伟.基于人工智能技术的地下水系统参数识别研究[D].天津：天津大学,2003.
[17]蒙媛.丰台区地下水模拟研究[D].合肥：合肥工业大学,2007.
[18]王浩,陆垂裕,秦大庸等.地下水数值计算与应用研究进展综述[J].地学前缘(中国地质大学(北京)),2010,17(6)：1-12.
[19]冯康.基于变分原理的差分格式[J].应用数学与计算数学,1965(2)：238-262.
[20]薛禹群,谢春红.水文地质学的数值法[M].北京：煤炭工业出版社,1979.
[21]王文科,李俊亭地下水流数值模拟的发展与展望[J].西北地质,1995,16(4)：52-56
[22]Pinder G. F., J. D. Bredehoeft. Application of the digital computer for aquifer evaluation. Water Resources,1968,4(5):1069-1093.
[23]孙纳正.地下水流的数学模型和数值方法[M].北京：地质出版社,1981.
[24]Courant R. Variational methods for a solution of problems of equilibrium and vibnltion. [J]Bull:Amer. Math. sec.1943,49:1-23
[25]林坜.大区域地下水流模拟研究及FEFLOW的建模方法-以华北平原为例[D].中国地质大 学,2006.
[26]刘晓民.干旱区典型地下水系统的动态模拟与模拟优化耦合管理模型研究[D].内蒙古农业大学,2009.
[27]刘玉珍.平原区地下水系统模拟与决策支持研究[D].大连理工大学,2008.
[28]冯绍元.解地下水中溶质运移问题的特征有限元法[J].武汉水利电力学院学报,1991,24(4)：410-417.
[29]谢春红,张勤,丁家平.稳定渗流的有限元计算新方法[J].水利学报,1995(12)：29-38.
[30]孙慧,周德亮.用无网格伽辽金法求解稳定地下水流问题[J].地下水,2008,30(6)：1-3.
[31]张宏仁,李俊亭.有限差分法与有限单元法在渗流问题中的对比[J].水文地质工程地质,1979(2)：50-55.
[32]Kitahara M. Boundary integral equation methods in eigenvalue problems of elastic dynamics and thin Plates[M]. New York:Elsevier Science Publishing Company Inc,1985.
[33]杨天行,刘金英,王运国.求解对流占优地下水水质模型的有限单元高精度广义迎风单元均衡格式及应用[J].长春地质学院,1996(2)：196-201.
[34]张祥伟,竹内邦良.大区域地下水模拟的理论和方法[J].水利学报,2004(6)：7-13.
[35]Mazzia A, Puttii M. Mixed f init e element and f inite volume discretization for heavy brine simulations in groundwater [J]. Journal of Computational and Applied Mathematics, 2002,147:191-213.
[36]Li S G, Mclanghlin D. A computationally practical method for stochastic groundwater modeling [J]. Advances in Water Resources,2003,26:1137-1148.
[37]王浩然,朱国荣,赵金熙.基于区域分解法的地下水有限元与边界元耦合模型[J].地质论评,2003,49(1)：48-52.
[38]魏胜利,苏伟杰.分布式地下水模拟模型MODFLOW介绍与应用[J].水利规划与设计,2011(3)：30-32.
[39]张红霞,宋文.地下水数值模拟的研究现状与展望[J].水利科技与经济,2007,13(11)：794-796.
[40]徐晓民,郭中小,贾利民等.地下水系统数值模拟的应用于展望[J].黑龙江水专学报,2010,37(3)：118-120.
[41]杨春玲,刑世录,韩爱中.地下水系统数值模拟的研究进展[J].内蒙古科技与经济,2007,21(151)：305-307.
[42]丁继红,周德亮,马生忠.国外地下水模拟软件的发展现状与趋势[J].勘察科学技术,2002(1)：37-42.
[43]王庆永,贾忠华,刘晓峰Visual MODFLOW及其在地下水模拟中的应用[J].水资源与水工程学报,2007,18(5)：90-92.
[44]Anon. SSG Software [Z]. Washington:The Scientific Software Group,2000.4-5.
[45]徐乐昌.地下水模拟常用软件介绍[J].铀矿冶,2002,21(1)：33-38.
[46]大连市统计局.2000-2010年大连市国民经济与社会发展统计公报[EB/OL]. 2011-03-01. http://www.stats.dl.gov.cn/class.jsp?dtype=80
[47]大连市城市水井管理所.大连市区及甘井子部分地段地下水勘察报告[R].大连：[出版者不详].1991.
[48]丁玲.大连市周水子地区海水入侵问题研究[D].大连：大连理工大学,2004.
[49]大连市水务局.大连市水资源公报[EB/OL]. 2011-05-12. http://www.swj.dl.gov.cn/html/szygb/index.html
[50]武强.大连大魏家岩溶地下水系统管理模型的研究[J].勘察科学技术.1990.1：23-27.
[51]王国利,周惠成.大连市地下水资源优化管理模型[J].大连理工大学学报,2001.141(1)：112-115.
[52]迟宝明,周彦章,易树平.大连地区沿海小流域地下水库开发潜力分析[J].工程勘察2006(5)：23-26.
[53]马驰.SIP和PCG2两种迭代法在地下水数值计算中的应用对比[J].西安科技学院学报2002.3 22卷1期.
[54]汪晓平,吴勇强,张宏林,等.ASP网络开发技术[M].北京：人民邮电出版社,2000.4-5;189-190.
[55]杨金忠,蔡树英,王旭升.地下水运动数学模型[M].北京：科学出版社,2009.