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

平原区区域平面二维地下水水质模研究型 罗 阳1 摘 要 本文通过讨论平面二维地下水水质模型，来研究唐山平原区地下水中硝酸 盐氮在平面二维问题中的运移过程及浓度变化的趋势。经模型计算得出在将来 唐山平原区地下水硝态氮高质区的范围是逐渐扩大的，硝态氮浓度逐渐升高 的，但浓度升高比较缓慢。 目前，海河流域一半以上的大中城市供水依赖地下水，占总用水量的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，多为HCO 3 3-Na 型水。咸水含水层单位涌水量7m /h ·m，为CL-Na 型水。在乐亭的东南 1水利部海河流域水资源保护局 2 部，由于滦河的近代沉积，有312Km 的浅层淡水覆盖于咸水体之上，含水层的厚度10~30m， 由南向北递增，单位涌水量10~20m3/h ·m。 二、 水质模型建立及求解 1．基本方程 地下水中的硝酸盐氮一般可视为保守物质。当溶质的运移尺度很大时，弥散作用的影响 相对于移流作用可以忽略。在忽略弥散的情况下，由质量守衡原理导出硝酸盐氮运移的基本 方程： ∂C ∂C ∂C ( ) ( ) w t, x , y −S t,x , y −u −v （1） ∂t ∂x ∂y 2.分裂格式 根据分裂法原理，式（1）可用下列分裂格式来逼近： 1 ∂CI ∂C ∂C   1  ⋅ −u −v t ∈ n∆t, n + ∆t （2）     2 ∂t ∂x ∂y