四柏城市供水水源地地下水流三维数值模拟.docxVIP

四柏城市供水水源地地下水流三维数值模拟 这四个柏树城市的供水水源位于陕西省靖边县西北部毛乌素沙地和洛杉矶高原之间的过渡地带。属于干旱和半干旱季风气候，生态环境相对脆弱。随着城镇化建设的快速发展, 靖边县城市用水供需矛盾日渐突出, 与此同时, 城镇化建设对四柏树水源地供水安全的潜在威胁也日益显现。为实现城镇安全供水, 在四柏树水源地范围内划定水源保护区已迫在眉睫。 1 地下水的补给、蒸发 区内主要含水层为第四系冲湖积孔隙潜水含水层、白垩系裂隙孔隙承压含水层, 在第四系底部稳定分布一层厚10～30 m的亚砂土、亚黏土层, 潜水水位略高于承压水。地下水主要接收大气降水入渗补给及少量的凝结水补给、农灌回归补给、地下水侧向径流补给。进入到含水层中的水一部分通过越流补给到白垩系承压含水层, 天然条件下主要通过潜水蒸发排泄, 人工开采也是主要排泄方式之一。 现状条件下, 区内地下水已遭受一定的污染, 个别地段地下水中硝酸根含量达70 mg/L、氯离子含量达1 192.43 mg/L、硫酸根含量达1 650.5 mg/L。 2 水源保护的划分原则和方法 2.1 级保护区范围及布井范围 城市饮用水地下水源保护主要是保护供水井不受污染。要求污染物在一定距离、一定时间内到达供水井时, 须经沿途介质净化, 使浓度降到目标含量。因此, 按照国家有关技术规范要求, 将靖边县四柏树水源地保护区划定为三级保护区: 范围为以地下水取水井为中心, 溶质质点迁移100 d的距离为半径所圈定的范围。当一级保护区范围过小时, 参照“英国水源保护区划定方法附加规定”, 一级保护区的半径不得小于 50 m。 范围为一级保护区之外, 溶质质点迁移1 000 d的距离为半径所圈定的范围。考虑到四柏树水源地含水介质为较纯净的砂及砂岩, 含水层自身的抗污染能力及净化能力差, 故根据水源井运行时间及国内外部分资料, 确定污染物运移时间为25 a的距离为二级保护区范围, 水源地的布井范围也应划定为二级保护区。 考虑到区内含水层自身的防污净化能力差, 同时也考虑到水源地长期为靖边县城供水, 将三级保护区确定为水源地的主要补给区 (开采井的水力截获区) 。 2.2 地下水流拦截段以地下水流拦截点 先采用美国地质调查局开发的地下水流三维有限差分计算程序Modflow 2005对计算区地下水流进行数值模拟。然后将水源地开采井作为水质点的起始位置, 采用美国地质调查局开发的专门用于水质点示踪模拟研究的程序Modpath, 利用其反向示踪功能确定向某点汇流的水质点在不同时刻的位置, 连接同一质点不同时刻的位置, 构成地下水流截获区, 按径流时间分别为100 d、25 a, 并结合水源地的布井范围划分一、二级保护区, 并将开采后的主要补给区划分为三级保护区。 3 地下水水文气象及微生物指标 考虑到避免集中供水水源地与城镇区及农灌区的相互影响, 今后水源地将逐步向北、西部的沙漠区迁移。为此, 计算区范围在工作区范围的基础上分别向四周延伸, 大地坐标X为4 160 000～4 179 000 m,Y为19 292 000～19 314 000 m。则计算区平面面积为418 km2;垂直方向为地面以下、侏罗系顶面以上之间的部分, 厚350～400 m。 计算区内主要供水目的层为白垩系厚层交错砂岩孔隙—裂隙承压含水层和第四系冲湖积砂层孔隙潜水含水层, 为非均质各向同性介质;白垩系下部地层为侏罗系砂泥岩互层, 其结构致密, 透水性差, 可视为区域隔水底板。计算区内主要供水目的层在区域内延伸广阔, 可视为无限含水层。现状条件下, 计算区东、西边界及北边界东段为流线边界, 开采影响区之外南边界、北边界西段对白垩系含水层来说可视为定水头边界。芦河从研究区南部流过, 河水与潜水水力联系较差, 芦河流经研究区河段的径流量较小且水面狭窄, 故芦河河水与地下水间的补排量可忽略不计。 研究区内地下水主要接受大气降水入渗补给及来自上游的地下水侧向径流补给, 天然条件下地下水由南、西向东流, 以潜水蒸发及地下水侧向径流方式排泄。在开采条件下, 开采范围内地下水由四周向开采中心流动, 影响范围之外径流方向为由西向东。 根据以上水文地质概念模型, 建立计算区的地下水三维稳定流数学模型如下: {??x(k?Η?x)+??y(k?Η?y)+??z(k?Η?z)=0(x,y,z)∈ΩΗ=zk?Η?n1=0}潜水浸润面边界Η(x,y,z)|Γ1=f1?Η?n|Γ2=0limr→rw(r∫2π0dθ∫hh′k?Η?rdz′)|wi=Qi 式中:H为地下水位标高, m;k为渗透系数, m/d;f1为第一类边界水位标高, m;Qi为第i眼开采井开采量, m3/d;r、θ、z′均为辅助柱坐标变量;n为二类边界外法线方向;n1为潜水面外法线方向;