基于level set有限元方法的微观水驱油数值模拟.pdfVIP

石油地质与工程 AND 2016年9月 PETROLEUMGEOLOGYENGINEERING 第30卷第5期 文章编号：1673—8217(2016)05—0091—06 基于LevelSet有限元方法的微观水驱油数值模拟 高亚军1，姜汉桥1，王硕亮2，刘传斌1，常元昊1，王依诚1 (1．中国石油大学(北京)石油工程教育部重点实验室，北京102249；2．中国地质大学(北京)能源学院) 摘要：在微观水驱油实验的基础上，建立了二维微观孔隙模型；引入LevelSet数学方法，结合N—S方程，建立了微 观两相渗流数学模型，借助有限元方法，进行水驱油两相数值模拟，研究微观水驱油动态特征。研究了孔喉非均质 性、润湿非均质性、油水两相黏度以及驱替压差对驱油效果的影响；通过对模拟结果中孔隙内流体的波及效率、两 相界面的运移速度以及优势渗流通道的分析发现，LevelSet方法能很好地处理各因素下油水两相界面的拓扑变 化，为微观水驱油机理的研究提供了一种新的技术方法。 关键词：水驱实验；数值模拟；微观驱替特征 中图分类号：TE311 文献标识码：A 许多学者在微观水驱油机理和两相流动特征的 油水两相驱数学模型 研究方面做了大量的实验，实验方法主要有岩心薄 1．1 Level 片平面驱替和玻璃刻蚀物理模拟等二维方法n1]， Set方法 以及驱替岩心CT扫描技术的三维方法[4‘5]，但以 假设地层中岩石孑L隙内水油两相流体均为不可 上方法实验模拟仿真成本高，操作较繁琐和困难，并 压缩流体，流动状态为层流。在层流两相流中，水平 且实验监测技术和图像获取精度有限，导致实验数 集界面自动生成界面运移的方程，在LevelSet方程 据误差较大。数值模拟方法不需要实际岩心参与， 中，当(D=0时的水平集轮廓线代表两相流体界面。 利用计算机对微观岩心孑L隙模型进行刻画，模拟实 对于此油水两相，在油中午=0，水中午=1。Level 际水驱油整个过程的演化规律，研究油水两相动态 Set方程可以被看作是表示水油两相流中水的体积 变化特征，具有便捷和可视化等优点。目前应用较 分数，因此，两相流体的界面运移方程可以用下式表 多的是孔隙网络模型[6。7]，但此方法一方面无法更示： 进一步清晰呈现孔道内流体两相界面运移状态，另 一方面无法直观获取各流体质点在任何时刻的压力 和流速变化。 式中：妒为油水两相界面轮廓线；7为方程求解中的 针对以上方法所存在的问题，本文提出基于 重新初始化参数；t为两相作用时间；U为流速；e为 Level 界面厚度，设置值一般小于计算模型中网格剖分的 Set(水平集)方法的微观水驱油两相数值模 拟，对孔隙内流动特征进行可视化研究。水平集方 最小单元。 借助于LevelSet函数，流体的物性可用下述 法是由Osher和Sethian[81提出的，最初主要应用于 智能控制、图像处理等方面，近年来经国外学者对求 方程表示，即： 解方法的不断改进[9叫0I，逐渐被用于气液或液液两 』P‘甲’一p。+‘pw—P0’甲 (2) 相流动的数值模拟研究领域中。有学者用此方法做 I肛(甲)2肛。+(弘。一p。)甲 过简单模型的油水两相模拟[1卜12]，但局限于理想孔式中：P。、P。分别为油和水的密度；p。、p。分别为油、 隙结构或油气两相驱替，均未与实验进行对比验证， 水的黏度。经过这样定义后，流体的密度和黏度就 且未涉及到微观孔隙内的油水两相模拟。为此笔者