渗流-水平井流耦合数学模型和数值模拟 A coupled seepage-horizontal well flow model and numerical experiments.pdf

第20卷第6期 水科学进展 V01．20．№．6 2009年11月 ADVANCES矾WATERSCIENCE Nov．．2009 渗流一水平井流耦合数学模型和数值模拟 董贵明1，束龙仓2，王茂枚2，鲁程鹏2，刘丽红2 (1．中国矿业大学资源与地球科学学院，江苏徐州笠1116；2．河海大学水文水资源与水利工程科学国家重点实验室，江苏南京210098) 摘要：在有侧向流人(流出)的水平井流运动和连续性方程的基础上，建立了渗流与水平井流耦合数学模型，主要针 对混合水头损失问题，进行了渗流与水平井流耦合数值计算，分析了常用的摩擦系数修正方法对计算结果的影响。 数值模拟计算表明：层流或者光滑紊流流态时，不同的摩擦系数修正方法对结果影响较小；粗糙紊流流态时，不同 的摩擦系数修正方法对结果影响较大。当水平井流是层流．光滑紊流流态时，使用了一种新的流态分界点识别方法 对水平井流态进行判断，计算结果证实了方法的可靠性、合理性。 关键词：渗流．水平井流耦合模型；数值模拟；流态；水头损失；摩擦系数 中图分类号：P641．2 文献标识码：A 文章编号：100l柳9l(2009)嘶083∞8 水平井是指井的透水(气)段．滤管呈水平放置的抽(注)流体的集水建筑物…，主要应用在世界各地的油藏 和天然气开采中[2。3]。近年来，美国和其他一些发达国家已开始将水平井技术应用到地下水资源开采、地下水 污染治理及斜坡排水治理工程等方面[¨]，尤其在地下水污染治理工作中(地下水污染物的原位抽吸系统)水平 井技术应用十分广泛【6J。 水平井水力特征的研究开始于20世纪50年代，研究者利用空间无限点源渗流基本方程，得到有限长度水 平井的解析公式，公式未考虑水平井中的流动对地下水渗流场的反作用及水平井中的流态，更未进行渗流与水 平井流的耦合求解[∞]。陈崇希等[眄J提出使用等效渗透系数方法进行渗流一管流的耦合计算方法，是目前进行 渗流一管流计算的主要方法，在国内外得到了较广泛的应用[n15]。等效渗透系数法是通过达西．魏斯巴赫方程将 渗流和水平井流使用统一的孔隙介质渗流方程进行描述并求数值解的一种方法。达西．魏斯巴赫方程适用于侧 壁无流人(流出)的明渠水流或管流【l6|，在水平井内流体流动过程中，不断有流体流入(流出)，这必将引起摩擦 系数的变化，摩擦系数的计算也就不能完全依赖于达西．魏斯巴赫方程，应考虑流体流入(流出)时产生的混合水 头损失对摩擦系数的影响，陈崇希等在计算中并未考虑这种影响。 最近几年，石油开采工程研究者已经注意到了水平井内流体运动的水头损失除了常规的沿程损失和加速损 失之外，还有流入水平井的流体与水平井内的流体在混合的时候产生的混合水头损失，并进行了试验和数值模 拟研究[17。19]。由于混合水头损失的存在，研究者认为应该对常规管流的摩擦系数进行修正，主要采用修正系 数和修正公式两种方法【扯21J，修正系数主要分为小于1的修正系数和大于1的修正系数；修正公式方法是在一 定试验条件下得出的经验公式。摩擦系数进行修正后，就可以不考虑混合水头损失的影响了。对摩擦系数进行 修正的数值模拟研究大都是针对单个水平井，很少有对渗流和管流进行耦合研究，并且缺乏对不同修正方法的 对比分析。 针对渗流一水平井流模拟时存在混合水头损失问题，本文在有侧向流入(流出)的水平井流运动模型的基础 上，建立了渗流．水平井流耦合数学模型，并且在水平井内是层流、光滑紊流、粗糙紊流和层流．光滑紊流流态 收稿日期：2008-09．04 基金项目：国家重点基础研究发展计划(9r73)资助项目(2006cB403204) 作者简介：董贵明(1979一)，男，黑龙江佳木斯人，博士，主要从事地下水资源评价与管理、地下水与生态环境研究。 B咖il：gu126．嗍 万方数据 第6期 董贵明，等：渗流．水平井流耦合数学模型和数值模拟 831 时，通过编制的程序进行了稳定流时的渗流．水平井流耦合数值计算研究，对常用的摩擦系数修正方法进行了 对比分析。 1 渗流一水平井流耦合数学模型 在建立渗流．水平井流耦合数学模型之前，首先推导有侧向流人(流出)水平井流数学