多孔介质两相流的局部守恒有限元分析.pdfVIP

第30卷第5期 计 算 物 理 V01．30，No．5 CHINESE OFCOMPUTATIONALPHYSICS Sep．，2013 2013年9月 JOURNAL 文章编号：1001-246X(2013)05-0667-08 多子L介质两相流的局部守恒有限元分析 张 娜， 姚 军+， 黄朝琴， 王月英 (中国石油大学(华东)石油工程学院，山东青岛266555) 摘要：用局部守恒有限元法研究多孔介质两相渗流问题．详细阐述局部守恒有限元法的基本原理，推导两相渗流 问题的局部守恒有限元计算格式并编制相应的计算程序．通过一维Buckley—Leverett两相渗流算例验证该方法的正 确性．应用局部守恒有限元法和混合有限元法分别对2个模型进行分析对比．计算结果表明局部守恒有限元法具 有良好的鲁棒性及适用性，相较于混合有限元法，处理过程简单，计算时间缩短，为标准有限元法应用于复杂渗流 问题提供了一种途径． 关键词：局部守恒；有限元；两相渗流；数值模拟 中图分类号：0357．3 文献标识码：A O 引言 流体通过多孔介质的流动是自然界广泛存在的现象．对于多孔介质中流体渗流问题的研究涉及到岩 土工程、石油工程、水文地质工程等多个领域，是一项非常重要的基础理论研究．局部守恒性是工程问题 数值模拟中的一个重要性质．局部守恒性在处理间断问题和复杂边界问题时都起着至关重要的作用．目 前在油藏数值模拟中主要采用混合有限元法、有限体积法以及近年来得到发展的间断有限元法¨。1，这 些方法都满足局部守恒性．但它们在处理大规模渗流问题时，不易处理且计算量大．传统有限元法尽管满 足隐式局部速度守恒，但只有当单元边界速度施加一定条件才能满足显式守恒．因此寻求直接由压力计 算流速并满足局部守恒的有限元法，具有重要的实际应用价值． finiteelemtent 混合有限元法(mixed 后分别进行混合有限元逼近．该方法目前具有较为成熟的多孔介质渗流理论基础，数值模拟结果准确度 速度分别进行插值，这在计算三维或者大规模渗流问题时，处理复杂且计算量大．此外，计算中速度需在 单元界面处满足法向分量连续，切向分量不连续；压力场要满足非连续性，且压力近似方程的阶数不能 采用． Galerkin 间断有限元法(discontinuous 一种数值方法，在解决流体动力学问题具有一定优势，被广泛应用于水动力学、气动力学、石油工程等 领域．DGM克服了传统有限元法的一些缺点，能灵活处理间断问题，离散单元定义简单，具有一致高精 度，而且便于实现并行计算．然而，伴随这些优势产生的是额外的单元自由度，计算量增大．这在处理 大规模渗流问题时导致计算量庞大，内存需求大，难以实现．为此，国内外诸多学者寻求一种具有传统 程相对容易实现． 收稿日期：2012—11—20；修回日期：2013—03—29 13CX06027A)资助项目 作者简介：张娜(1986一)，女，博士生，研究方向：复杂介质油气渗流理论和油藏数值模拟，E．mail：znpear@163．com }通讯作者：姚军，E-mail：RCOGFR—UPC@126．tom 万方数据 计 算 物 理 第30卷 conservative 限元法可以实现局部守恒．在此基础上，Nithiarasu首次提出局部守恒有限元法(10cally finiteelement galerkin 度分别进行插值近似，近似处理与有限元相同，处理过程简单；不同于间断有限元在计算相邻单元边界 流量时需要计算额外的方程，它只需添加一个后处理过程就可以得到准确的边界流量，实现压力和速 度局部守恒．LCG基于传统有限元法，编程只需对原始有限元代码进行部分修改而不需重新整体编 写，易于实现．目前LCG法的研究主要集中于对流扩散方程、N