北京平谷平原区地下水流建模的若干问题探讨.pdfVIP

北京平谷平原区地下水流建 模的若干问题探讨 黎 明 李国敏1 周 东2 100089) (1．中国科学院地质与地球物理研究所100029；2．北京市水文总站 摘要：地下水数值模拟是地下水资源评价和管理的有效方法。许多地下水数值模拟中忽视了对 水文地质条件的合理概化，或者缺乏符合机理的模拟技术。本文针对当前较普遍的问题，以北京 平谷平原区地下水流动问题为例。分析了建立三维地下水流模型的必要性，引入地表产流模型改 进了边界流量的确定方法，提出了利用地下水CFCs定年数据估算含水层渗透系数的方法。 关键词：地下水；数值模拟；边界参数 1 引言 地下水数值模拟是地下水资源评价和管理的重要工具手段。利用模型，可以把通过不 同方法(调查、监测、钻探、遥感等)得到的数据进行统一的系统分析，用来描述和分析 各种条件下复杂地下水系统的运动规律，进而为制定合理的水资源开发模式提供依据。然 而当前许多地下水数值模拟中往往存在以下问题：不重视对水文地质条件的深入分析和合 理概化，而片面强调数值计算方法和追求拟合结果；忽视了边界条件或源汇项的刻画；河 流、降雨滞后补给潜水的刻画也是建模过程中易被忽视的问题。以上这些问题都会导致模 拟失真。本文针对这些当前较普遍的问题，以平谷盆地为例，进行了相应的分析和改进。 2研究区概况 平谷位于北京市郊最东部，其东、北、西和东南均为中低山丘陵环抱，仅西南有一出 为山地，地势总体东北高，西南低[1]。 该区属暖温带大陆性季风气候，多年平均蒸发强度1693．8mm，多年平均降水量为 盆地。区内河流属海河流域的蓟运河水系：洵河由东北流向西南；洳河北流向南，在盆地 中部汇入洵河；之后贯穿盆地中南部流出境。 平谷盆地内广泛发育新生界第四系地层，以洵河和洳河冲洪积成因为主，沉积厚度由 山前至盆地中心由薄变厚，最厚为550余m。冲洪积扇顶部至中下部，沉积物颗粒由粗变 细，由单层变为多层结构，上部35～50m主要由粉土和粉质黏土组成，下部为巨厚的砾 卵石、砂卵石层。王都庄以东和许家务以北为单层巨厚卵石含漂石、砾卵石层，盆地中 19 部地区为多层砂砾卵石层，各含水层间存在着密切的水力联系。 地下水补给、径流、排泄具有明显的统一性和分带性规律：盆地之东自山东庄～夏各 庄一线以东为洵河山前冲洪积扇潜水，西南的西沥津和龙家务为由潜水转承压水的地下水 溢出带；盆地之西自许家务以北为洳河山前冲洪积扇潜水，许家务一中桥一带为地下水溢 出带；地下水溢出带向南逐步变为承压水区。地下水主要补给来源是大气降水入渗补给、 海子水库渗漏、渠系和田间灌溉人渗补给等，其中大气降水在山区人渗补给后，除当地工 农业用水外，其余部分以侧向径流的形式补给平原区，为平原区的主要补给源。地下水的 排泄主要是工农业开采用水，地下水溢出量和侧向径流量也是排泄的重要途径。 3建模中的几个问题 3．1基本模型 由上述水文地质条件分析，研究区为一个复杂的含水系统。地下水自东北部山前冲洪 积扇向西南部流动过程中，除了在水平方向上发生交换外，也存在垂向上的交换。在盆地 北部冲洪积扇上部，虽然含水介质为巨厚的比较均一的含卵砾石砂层，但其上的河流是一 条切割较浅的非完整河，河水非连续地垂直人渗补给地下水，这种条件下地下水的垂直分 流速不可忽略。河流、水库底部的地下水存在明显的垂向流动；含水层间弱透水层中的地 下水也以垂向流动为主；非完整抽水井附近的地下水流存在垂直流速分量。可见研究区地 下水流动具有明显的三维特征。考虑到具体条件和建模目标——区域地下水资源评价与管 理，实际含水系统的层数很多而模型分层又不宜过多，因此概化后的“弱透水层”中含有 一定数量的砂层，即“弱透水层”是各向异性的。陈崇希等(2001)[2]指出：“准三维流模 型‘误差小于5％’的结论仅适用于弱透水层为各向同性介质；对于弱透水层的各向异性 比为10时，误差已超过27．7％，而且随着弱透水层单位储水系数岸7，的减小和模拟时间 定性见解。基于上述研究的成果，本研究建立三维不稳定流模型，能够比二维模型更为客 观地刻画研究区地下