蒸汽热采稠油数值模拟的自适应网格法.pdfVIP

中国工程热物理学会 多相流 学术会议论文 编号：066137 蒸汽热采稠油数值模拟的自适应网格法 王晓宏 中国科学技术大学热科学和能源工程系。安徽台肥，230026 电话 Email：xhwaug@ustc．edu．ell 摘要 针对提高蒸汽热采稠油流场的计算速度，研究了多孔介质中多相流动的自适应网格(AMR)方法． 跨越小同种类岩石界面，油、水、气等相饱和度发生不连续跳跃。这将导致自适应网格法中的主要步 骤，如细粒化运算和网格细分等出现困难。尽管相饱和度在跨越不同种类岩石界面时不连续，但在每 个粗格子内不同相之间的相对渗透率之比可以近似看成常数，由此出发提出了解决办法。本文还介绍 多孔介质中多相流动的A躲方法的一些基本步骤。 关键词：多相流动：渗流；自适应网格法：蒸汽热采；稠油。 1．简介 稠油在自然界石油储存量中占有较大的比重，如我国的辽河、吉林、大港和大庆等 油田都有相当的储量。稠油在常温下枯度极高，近似为固体，随着温度的升高，粘度迅 速降低，在高温卜则为粘度较低的易流动液体。由于稠油上述物理性质，蒸汽热采成为 稠油开采的重要方式之一[1]。我国蒸汽热采的研究取得了很多的成果，在稠油的开采中 已得到应用[2-7]。近年来．一种新的稠油热采方式——蒸汽辅助重力驱油技术正在得 到研究和发展[8-10]，该技术在我国的稠油开采中已得到应用[Ii一12]。高温近饱和蒸 汽由注入井向地层中灌注，经过一段时间后，在蒸汽注入井的上方形成…个缓慢向前推 进的被加热区域，在该区域内水和蒸汽共存。加热区域和未被加热区域之间存在一个变 化非常剧烈的很狭窄的锋面，加热了的粘度较低的稠油、水和蒸汽从该加热区域被推向 锋面，在该锋面内．在重力作用下稠汕和水将流向产油井(如图l所示)。 圈I燕汽辅助重力驱油乐意图 564 蒸汽辅助重力驱油技术具有效率高、费用低等优点，正在得剑逐步推广。蒸汽热采 稠油是一个带有相变的多孔介质中多相流动问题，实验模拟的难度较大，冈而数值模拟 成为预估蒸汽热采技术效率和产量的重要手段。蒸汽热采问题存在一个物理量变化非常 剧烈的很狭窄的锋面，在锋面处计算网格必须很小．如果整个开采区域在空间上采用均 匀网格，计算量就会很大，相关的计算研究费用将难以承受。为了提高计算速度。动态 网格(dynamical 法中物理量变化较快的区域采用精细网格，而在物理量变化较’熳的区域采用租网格，这 样可以有效地提高计算速度。在上述动态网格法中，对所有的网格采用相同的计算时间 步长，由于计算流体力学中CFL条件的限制，对于精细网格必须采用小的时间步长，这 样如果对整个开采区域在时间上采用相同的计算步长，计算量仍较大。另外一方面，若 开采区域分布着不同种类输运性质的岩石．已有的动态网格法并未对此予以研究。为了 Mesh 有效地提高计算速度，本文介绍利用自适应网格法(AdaptiveRefinement)计算蒸 汽热采稠油问题。 个很复杂的科学问题，要将AMR方法应用于计算该问题的流场，很多基础性的数学物理 问题需要予以深入的研究。现有AMR方法中的一些算法需要加以改进。 在AMR方法中，网格系统是多层的，假设最低阶的层为最租的网格，那么最高阶的 层即为最精细的网格。在物理量变化慢的地方采用粗网格，在物理量变化快的地方采用 细网格，对于不同尺度的网格采用不同的计算时问步长，这样可咀大幅度地提高计算速 度。针对特定的物理问题，需要对选定的物理量设定阀值，当这些物理量的空间变化率 大于某个阀值时，网格将被细分。蒸汽热采问题存在温度变化很剧烈的热锋面，温度在 空间的变化率也就自然成为表征热锋面的物理量；同时蒸汽热采问题是一个多相流动问 题。存在相锋面。在现有文献中，iⅡ用饱和度在空闯的变化率来表征锋面[19—21]，这只 适用于同一种岩石的油藏区