地下结构对渗流场阻隔问题的解析～半解析法-北京勘察设计研究院.PDF

地下结构对渗流场阻隔问题的解析～半解析法 王军辉 韩煊 周宏磊 张在明 （北京市勘察设计研究院有限公司，北京 100038） 摘要：针对当前地下空间开发引起的地下水渗流场变化量化评价需要，本文就渗流场受地下结构的阻隔后的稳定流状态下的 流量及各部位（上游段、阻隔段和下游段）的水位计算公式进行了探讨，分别导出了上、下阻隔两种情况下的承压水问题的 解析解和潜水问题的半解析解，并利用数值试验直观解释了一些阻隔现象，以及获得了小水力梯度条件下潜水问题半解析解 可以用同条件下的承压水问题解析解来近似的结论。 关键词：地下结构；渗流场；阻隔作用；分段积分、端点连续；解析～半解析法 中图分类号： 文献标识码：A 文章编号：1000-3665 （2009 ）02-0000-00 1 前言 段）的解析解（针对承压水）和半解析解（针对潜 水），同时，通过相关数值试验解释了一些阻隔现 近年来，国内许多大中型城市正处在地下空间 象以及小水力梯度条件下潜水问题半解析解和承 开发热潮期，地下结构对地下水流场可能产生阻隔 ～ 压水问题解析解之间近似数值关系。 作用也日益受工程及环境领域的高度重视[1 3] 。这 些关注的焦点之一就是该类问题的定量评价，因为 2 地下结构阻水后果及其评价指标 这无论从工程角度还是从地下水资源和环境角度 (1)地下水位的壅高～降低 来看都是十分重要的。实际上，从地下水动力学角 显然地下结构对渗流场阻隔作用最直观的变 度来看，这类问题可近似等效地归结为地下水过水 化就是造成地下结构的迎水面水位壅高和背水面 断面的减小，或者更为简单地说，从剖面二维角度 的水位降低，为此引入一个参数－水位增量 ΔH ， （地下结构垂直流线方向无限延伸条件下）来看， 即： 即含水层厚度的减小。而针对变厚度含水层剖面二 ΔH ＝H −H 0 （1） 维流问题的计算方法，以前苏联为代表一些学者作 了大量而细致的探讨工作[4～6] ，但这些工作大多数 式中：H——阻隔后的水位分布函数（m ）； 是基于含水层厚度呈连续变化情况下的，而地下空 H0——阻隔前的水位分布函数（m ）。 间引起含水层的厚度变化往往是非连续的、阶梯式 （2 ）流量折减 的，因此这些公式[4～6]在地下结构对渗流阻隔这一 地下结构阻隔作用的另一个重要后果就是区 新问题定量研究中也就显得力不从心。数值解虽然 域性地下水径流量的减小，为评价该类问题引入流 能够模拟出各种复杂条件下的结果，但难以对其机 量折减系数 η，即： 理进行严格意义上地探讨[3]，因为在理论研究方面， q 即使是十分先进数值解也不能完全替代解析解的 η＝ 0 （2 ） 作用的[7] 。 q 基于上述原因，本文从渗流理论出发，建立了 式中：q——阻隔后的单宽流量（m2/d ）； 数学模型，采用“分段积分、端点连续” 的数学技巧