海河流域地下水开采初期数值模拟及水量平衡分析-论文.docVIP

海河流域地下水开采初期数值模拟及水量平衡分析（张晓明 薛丽娟 张奇 李建新） 提要:海河流域平原区地下水长期过度开采已引起了严重的环境问题,水资源的不断枯竭也制约着当地经济的可持续发展,成为国家一直关注的问题。作为地下水管理的方法之一,建立了地下水数值模型。对海河流域平原区地下水开采初期状态作了模拟分析,模拟结果显示:该区域地下水年补给量为156.亿m,其中,山前侧渗17.亿m,降雨入渗补给136亿m,河道入渗3. 1亿m;地下水的排泄途径为河道基流88.亿m,蒸散发50. 1亿m,排泄入海15.亿m,黄河内滩侧渗排出1.亿m,总排泄量为156亿m,地下水收支基本平衡。地下水的总体,吻合于地形坡度变化。模拟结果为水资源管理和水环境修复提供数据依据。模型也为下一步地下水大量开采状态的模拟提供初始水头条件、地质物理参数和边界条件。 关键词:地下水模型;MODFLOW;开采初期状态;海河流域 : TV213. 1 　　　　　　　文献标识码: A ,研究区总面积13. 1万km2。自上世纪70年代,海河流域人口和社会经济规模迅猛发展,加之连年干旱少雨,流域内浅层地下水资源已经成为社会发展和经济建设的主要供水水源[ 1 ] ,因此海河流域地下水超采日益严重。2000年海河平原浅层淡水总超采量为45亿m,超采总面积已达59550km,占平原区中有淡水区面积的63. 1%。由于长期超采地下水,已经导致唐山北部、XX西南、河北太行山前平原部分地区含水层被疏干[ ] 。河北太行山前的浅层水漏斗已经发展成为彼此相连的复合漏斗群,,可利用前景暗淡[ ] 。 ,如地面沉降、湿地湖泊萎缩、河流干涸、海水入侵、咸水位下移等等[ ] 。加强地下水管理,是目前海河流域水利所面临的一项重要工作。国内外很多学者专家提出了恢复地下水的各种措施,包括人工回灌、跨流域调水等,但是这些工程措施并不能完全解决现存的严重问题[ ] 。海河流域地下水资源管理的核心是保证地下水资源长期持续利用,地下水模型系统的建设将是管好地下水的关键技术之一。模型可以模拟地下水漏斗的发展过程,并能预测未来年份地下水的开采程度,从而为水资源规划提供了保障。地下水模型[ - ] 。 ,由于对初始水文地质环境的重视程度不够,且70年代以前的地下水补给、排泄、水位等方面的观测资料较少,以往的地下水研究大多以70年代后期地下水环境为研究对象,缺少对开采初期水文地质环境的探求。因此,大多数模型工作者以70年代后期地下水超采期的某一时刻为初始条件,建立了非稳定流模型。这个初始解不是真正的稳定条件下的初始值,所建的非稳定流模型也不能反映地下水环境变迁的全貌。另外,地质结构的研究一直是开展地下水研究的难点之一,但是由于本地区面积辽阔,地层分布差异性较大,尚没有一个能全面反映该区域地MODFLOW模型,以地下水模拟系统GMS[ ]为平台,较为理想地解决了地层模型建设的难题,并建立了地下水稳定流模型,模拟海河流域平原区开采初期的水文地质条件。研究目的主要有,全面认知本地区水文地质条件和初始的水量均衡场,在计算机上模拟区域地层结构和水量均衡关系,为非稳定流模拟提供初始水头条件、地质参数和其它计算条件,为该地区的水生态环境修复和水资源规划提供参考数据。 1 浅层地下水系统开采初期状态 :山前倾斜平原区、中部冲积平原区和滨海平原区。除个别地点有基岩出露外,绝大部分为第四纪松散物质覆盖。华北平原第四系是一套几何形态复杂的多种沉积类型交叉叠置的含水岩系,含水层岩性、结构、厚度等具有水平变化规律。在山前平原含水层呈扇状结构,扇轴含水层岩性以砾石卵石为主,厚度大;扇间含水层粒度变细,厚度变薄。在中部平原含水层逐渐过渡为湖相沉积穿插河流沉积的舌状结构,含水层岩性以中细砂为主,厚度在靠山前平原方向变薄,向滨海方向又略变厚。向东部、南部的滨海平原含水层又过渡,含水层岩性以粉细砂为主,厚度又变薄[ 10 ] 。 ,地下水位等值线与地表等高线基本相似。地下水的流场也是从北、西、南三面汇入渤海湾。以196年的地下水埋深为参照,分析得出各个分区的原始地下水位如下:滏阳河- 子牙河一线以北的北部平原,降雨相对于南部平原少,河网不如南部平原密集,河流多为季节性河流,不能常年保持河水位,因此北部平原的地下水位比南部平原低。在山前区,地下水埋深大约为2m (燕山平原) 、3m (太行山前平原) ;中部平原地下水位埋深大约为1.m;滨海平原小于1m。滏阳河- 子牙河一线以南的南部平原,,河网纵横,地表水充沛,地势较为平缓、退水慢,因此地下水位普遍高于北部平原。山前区地下水位埋深大约1. -m;徒骇马颊平原、黑龙港平原埋深小于1.m,局部地区多有泉水,盐碱地发育。平原区地下水在大量开采以前,第一含水层组和全淡水区的第二含水层组的水位多年动态变化与