涉县东风湖泉域地下水开采数值模拟研究.docxVIP

涉县东风湖泉域地下水开采数值模拟研究 1 地下水特征及开采行为 涉县东部的登湖泉面积位于山西省山西省临县西部，山西省黎城县交界处。水域面积约92 400 hm2, 其中, 69.70%位于河北省涉县境内, 30.30%位于山西黎城县境内。 涉县交通发达, 北有邯 (郸) 长 (治) 铁路穿越全境, 南有G22青兰高速及309国道和202省道跨过, 是晋煤东运的必经之地。县内火车停靠站点9个, 公路里程222.79 km (图1) 。 岩溶裂隙地下水具有以下4个特征。 (1) 大气降水补给量有不断减少的趋势, 区内地下水天然资源量已从多年平均4.553 m3·s-1(1956—2000年) 减为4.522 m3·s-1(1991—2000年) 。 (2) 人工开采量有不断增加的趋势。区内地下水资源的开采20世纪80年代更是急剧增加, 地下水开采量已超过天然排泄量。原有的天然状态已经遭受破坏, 机井开采占地下水排泄的比例越来越高, 且机井多位于泉区附近。这种集中式开采并形成大降深降漏落斗的开采方式对环境及区域水文地质条件具有一定的危害性。 (3) 区域地下水位总体呈下降趋势, 泉区及其附近集中开采区降落漏斗范围变化不大, 而降深逐年增大。 (4) 东风湖泉群流量呈逐年衰减趋势, 由于天然补给量的减少, 人工开采量的逐年增长, 泉流量呈逐年衰减的趋势。 2 三维渗透值的模拟 2.1 地质背景及其概念模型 2.1.1 金融边界两种边界 东风湖泉域的补给主要靠该区雨水入渗和清漳河水下渗补给。据该区水文地质资料, 本次数值模拟西边界确定为峪里沟, 北侧由于下寒武突出, 南部中寒武出露, 两侧均有一定的侧向流入, 但流量很小, 属于二类流量边界, 故在本次模拟中定为单宽流量。西北侧有震旦系和寒武系下统泥页岩突出, 为隔水层, 也定为单宽流量。东北侧也属于二类流量的侧向补给边界, 也以单宽流量形式给出。东侧、南侧边界与泉域边界重合, 其中东界、南界的中东段为隔水边界, 以二类零流量边界给出。南界西段为地下分水岭边界, 也以二类零流量边界给出。 数值模拟区域计算面积52 485 hm2, 其中灰岩裸露区面积32 610 hm2, 岩浆岩分布区4 964hm2, 第四系覆盖区14 911 hm2。 2.1.2 统一计算目的层 本次数值模拟计算, 依据现有水文地质资料情况将寒武、奥陶含水层列为统一计算目的层, 寒武系下统确定为统一隔水底板, 基本上防止了区内岩溶裂隙地下水的径流量对计算的干扰。 2.1.3 水系隔水层的水力联系 在涉县广大盆地区, 模拟区内目的层上部为第四系孔隙水。该含水层薄而不稳定, 含水性不太好, 另外该含水层底部还有厚达近百米的粘土层, 该层连续性好且非常稳定, 粘土本身隔水性能又很好, 故目的层与第四系孔隙水含水层的水力联系在计算时可以忽略不计。在清漳河河谷区, 第四系孔隙水层富水性强, 河底与该含水层之间的隔水层薄而不稳定, 故该含水层可直接接受清漳河水的渗漏补给, 另外, 河水也可补给到岩溶裂隙层。根据河流断面多年的水文资料显示, 该水力联系在模型中计算时以源汇项的形式给出。其下伏的寒武系下统因其富水性极弱的特性, 区域上属于隔水层。构成模拟区内相对稳定的隔水底板。 2.1.4 径流区划分 根据岩性结构、水动力和水化学、地下水动态四特征显示。在模拟区内, 评价目的层在平面上可划分出两个主要强径流区, 其余均为弱—中等径流区。而数值模型的分区便以径流区的划分为依据。 2.1.5 盆地中部径流 (1) 复杂的地下水类型。模拟区内西北部存在大面积碳酸盐岩, 该岩层透水性好, 故该区域属面状入渗补给区, 地下水属于潜水;涉县盆地第四系覆盖层边缘径流区, 也为潜水;只有盆地中部径流区 (靠近涉县断裂带) 为无压至半承压层间水的过度类型, 为承压水。 (2) 渗流形态。本区岩溶裂隙含水介质因为岩溶裂隙发育, 具有很好的连通性, 又因其非均质各向异性特征, 可以作为非均质各向异性渗透介质进行计算。虽然该介质在垂向上具有一定的分带性, 但考虑到介质内各含水层组的构造、节理裂隙、岩溶裂隙发育完全和水力联系密切程度, 同时也为了保持其宏观性和客观性, 计算时算作统一含水体非均质各向异性渗透介质。可以按照符合达西定律的层流进行渗流计算。 2.2 三维渗透数学模型与解决方案 2.2.1 水—三维渗流模型原理 根据以上模拟区三维渗流概念模型, 建立非均质各向异性, 具有二类边界条件, 潜水—承压水三维非稳定流数学模型如下: 式中:Kx, Ky, Kz为介质渗透系数 (m·d-1) ;h—地下水位 (m) ;ho—初始水位;Г2—二类边界;—立体计算域;s—储水系数;q—二类边界流量 (m3·d-1) ;t—时间;Q—源汇项。 2.2.2 利用差