水平井或旁河垂直井地下水流物理模拟.docxVIP

水平井或旁河垂直井地下水流物理模拟 1 地下水流态及流态问题 井孔是监测地下水流量、收集地下水主要工具和防止、控制地下水环境恶化的主要方法。因此,井孔-含水系统问题是当今水文地质学研究中最重要的课题之一。然而, 长期以来,水文地质模型中对井孔的刻画, 基本上都是引用“热传导”中的“线汇理论”,除个别采用等水头的第一类井孔边界条件,基本上都采用等强度线汇(均匀入流)的第二类井孔边界条件,诚然,就渗流而言,井孔的滤管壁面是渗流的边界,水通过滤管进入井管之后就不再是渗流了,既然滤管壁面是渗流的边界,就应该给出此处的边界条件。这可能是长期以来的思维逻辑,于是想方设法给出滤管上的流量分布或水头分布。 地下水流入井孔与地下水流入(排泄到)河流,在形式上很相似,但渗流边界条件往往是不同的。对于后者来说,地下水流入(排泄到)河流的流量通常是河流流量的小部分(否则,就要按地下水流入井孔一样处理),一般并不因地下水的流入而明显地改变河流的水位,即不至于明显地改变地下水的边界条件。地下水流入井孔则不同,井孔中的水流全部由地下水流入而成,而且井孔过水断面狭小,水流阻力大,井管中的水流流态(层流、紊流)将明显地影响水头分布和流量分布。即,一般地说,不可能预先给定井孔滤管壁面处的边界条件,特别是对于常见的贯穿似层状含水层(含透镜体)的多层井,更是不可能预先分配各层的流量(对于油田开发,该问题同样未能解决)。 基于此,陈崇希等提出将井管的水流与渗流一起作为模拟对象——“渗流-管流耦合模型”,从而将边界从滤管壁移至井口,实际的井口流量或井中水位就是其边界条件,沿井孔滤管的流量和水头的分布是作为模拟预测的结果给出,而不是预先给定。 Therrien和Sudicky在假定井孔中水流遵守线性流的条件下,建立了井孔-含水系统的数值模型。陈崇希等为建立地下水向混合井(多层井)流动的模型,提出“渗流-管流耦合模型” 刻画井孔-含水层水流系统。该研究全面地考虑了井管中不同流态(包括线性流和非线性流)条件下水流水头损失与流速间的关系,并成功地用有限元法计算了广西北海市龙潭水源地实例(双含水层混合井抽水试验确定各层水文地质参数的研究),也证明了MODFLOW软件建议按含水层的导水系数比值来人为预先确定混合井各层的流量是缺乏理论依据的,实例表明其误差达70.0%～84.6%。陈崇希利用“渗流-管流耦合模型”建立了岩溶管道-裂隙-孔隙三重介质地下水流模型,并给出了一个理想模型的计算结果,再现了岩溶泉流量的动态特征。Cheng Jianmei将其用于广西北山矿区中一个含有岩溶管道的地下水流实例研究。陈崇希、林敏等将其用于河南省郑州市黄河旁的水源地中三个含水层混合抽水的有限元分析。Chen Chongxi将其用于地下水向混合井(多层井)流动的模拟研究。关于地下水水平井流和辐射井流的研究(可视为1组水平井),除个别外,研究者们基本上都采用线汇来刻画水平井管的作用,例如M S Hantush、M S Beljin、T G Cleveland、Zhan Hongbin等。然而,用线汇刻画水平井管的适用性并未得到严格证明。陈崇希和万军伟, 陈崇希、万军伟和Zhan Hongbin成功地将“渗流-管流耦合模型”用于研究地下水向水平井的流动问题。算例表明,水平井管既不是均匀入流也不是等降深(水头)分布,井管中两者都相差四倍多。万军伟等(2003)研究了水平井解析解的应用条件。最近,陈崇希将“渗流-管流耦合模型”用来研究三维地下水流中常规观测孔水位的形成,得出某些与传统(如M S Hantush以及前苏联的Бочевер和Веригин)不同的见解。 “渗流-管流耦合模型”从理想模型到实例应用,已用于井孔-含水层系统若干复杂问题。本文的目的是通过物理模拟(砂箱试验)获得一批实测数据,以检验其可靠性。 2 “流预览模型”砂箱试验 2.1 砂箱及井管设置 为了便于砂箱物理模型的制作、试验及其相应的数值模拟,并尽可能多地再现井管水流的多流态特点,我们设计了一个河床下矩形含水层,在其下部设置一口水平井。实际上,这个物理模型也可以视为一口横卧的傍河垂直井模型。采取渐变流量抽水(放水)的砂箱物理模型试验。 (1)砂箱几何尺寸的设计如图1所示,长(x轴)453cm,宽(y轴)50cm,高(z轴)120cm,顶部敞开,以便装砂和供水。为了避免箱体变形,箱体的材料采用10mm的钢板,其中有一侧采用5mm厚的钢化玻璃,以便试验过程中观察砂箱内的情况。 (2)边界条件。为简单起见,矩形砂箱的4个侧面和底面均为平面隔水边界。顶面边界条件有几个方案可供选择。若选取隔水边界,则井孔的抽水量仅依靠弹性储存量的释放来提供,将导致“含水层”中水头过快地下降,而“含水层”的初始水头在实验室里难以提得过高,使得试验时间过短,不利于模