平原区区域平面二维地下水水质模型研究.docVIP

平原区区域平面二维地下水水质模研究型 罗 阳 摘 要 本文通过讨论平面二维地下水水质模型，来研究唐山平原区地下水中硝酸盐氮在平面二维问题中的运移过程及浓度变化的趋势。经模型计算得出在将来唐山平原区地下水硝态氮高质区的范围是逐渐扩大的，硝态氮浓度逐渐升高的，但浓度升高比较缓慢。 目前，海河流域一半以上的大中城市供水依赖地下水，占总用水量的58%。其中唐山市用水量中地下水也占较大比例，而该平原地区地下水中硝态氮的浓度相对比较高，所以在该地区进行硝态氮的研究具有非常重要的意义。 研究区基本情况 唐山市地处河北省东部，位于东经117°30′至119°18′，北纬39°03′至40°35′之间，北靠长城，南邻渤海，东以滦河、青龙河为界，与秦皇岛市接壤。全市下辖八县二市五个区，两个县级国营农场。 唐山市范围内主要地层有太古界变质岩、长城系、寒武系、奥陶系、石炭系、二叠系及第四系。平原区基岩大部分被第四系所覆盖，由北向南逐渐增厚，由数米增至数百米，沿海地区增至千余米。平原区含水岩主要是由第四系松散沉积物构成的，按其成因可进一步分为山前冲洪积倾斜平原松散岩类孔隙含水岩组及冲积海积低平原松散岩类孔隙含水岩组。 山前冲洪积倾斜平原松散岩类孔隙含水岩组是由还乡河、陡河、沙河及滦河所形成的规模大小不同的冲洪积扇所构成。含水层由卵砾石、砂含卵砾石组成，颗粒粗，厚度大，一般单层厚10~20m；由单层过渡到多层，含水体呈扇状分布，含水层上覆岩性及含水层之间多为粘质砂土及砂质粘土相隔，有利于垂向渗透和降雨补给。单位涌水量30~50m3/h·m，水量丰富，水质良好，多为低矿化度重碳酸盐水，水动力特性为潜水-微承压水。矿化度小于2g/L，面积5153Km2。 冲积海积低平原松散岩类孔隙含水岩组分布于矿化度大于2g/L的咸淡水界线以南，面积3101Km2，处于区内各河下游。水质特征为上咸下淡双层结构。其含水层为全淡水区含水层的延伸，颗粒由粗变细，由细砂、粉砂组成，局部中砂。含水层之间粘性增大。侧向径流与垂直渗透缓慢。由于海侵及潮汐作用，水化学类型复杂。咸水体厚度由西向东逐渐增厚。咸水体之下为深层淡水含水层，顶板埋深40~200m，单位涌水量10~30m3/h·m，矿化度小于2g/L，多为HCO3-Na型水。咸水含水层单位涌水量7m3/h·m，为CL-Na型水。在乐亭的东南部，由于滦河的近代沉积，有312Km2的浅层淡水覆盖于咸水体之上，含水层的厚度10~30m，由南向北递增，单位涌水量10~20m3/h·m。 水质模型建立及求解 1．基本方程 地下水中的硝酸盐氮一般可视为保守物质。当溶质的运移尺度很大时，弥散作用的影响相对于移流作用可以忽略。在忽略弥散的情况下，由质量守衡原理导出硝酸盐氮运移的基本方程：  （1） 2.分裂格式 根据分裂法原理，式（1）可用下列分裂格式来逼近： （2） （3） 3．移流项处理 （1）流速计算 原理：按多边形内平均，如节点q的 流速由包围该点的三角形单元构成的多边 形内平均。 图1 在每个三角形单元内，节点q在x， y方向上的流速分别为： 图2 （4） 节点q的流速为： (5) （2）反求质点位置 式（2）表明在时段内，当流体质点迁移时，其浓度CI值保持不变，也即沿特征线保持不变。其特征线为： （6） XP、YP为质点在时段初的坐标。 应用特征线法求解式（2），认为时段末质点正好位于节点上，推求该质点在时段初的位置。 （3）由质点位置内插C 根据质点在时段初的位置及周围节点 的已知浓度推求质点所具有的浓度。判断 出点所在的单元，则点p所具有的 图3 浓度由所在单元三节点的浓度内插求得。应用有限元方法中三角单元插值函数，有： （7） （8） 由以上讨论可知，时间步长的选择应使得p点落在节点q周围的单元内。若点p在工作区边界之外，则令节点的浓度C不