学位论文 > 优秀研究生学位论文题录展示

油藏数值模拟自适应网格法中若干问题研究

作　者: 吴斌
导　师: 王晓宏
学　校: 中国科学技术大学
专　业: 工程热物理
关键词: 多相渗流自适应网格多种岩石等效渗透率重整化方法标度指数相关长度三维空间标量张量两相驱替非活塞式驱替
分类号: TE319
类　型: 博士论文
年　份: 2012年
下　载: 74次
引　用: 0次
阅　读: 论文下载

内容摘要

在石油生产工业中,经常通过注入水或其他流体,驱替地层中的原油,从而达到产油和提高采收率的目的。为了准确预报产量和尽可能提高采收效率,数值分析工具必不可少。因此长期以来,对多孔介质中多相渗流驱替问题的数值研究也一直受到广泛关注。对于含有相变的多相渗流的数值模拟,通常存在着一个明显的两相区与三相区的分界面,界面附近的物理量变化非常剧烈。自适应网格法就是一种处理带有急剧变化锋面的数值算法。在物理量变化剧烈的区域采用细网格,而在物理量变化缓慢的区域采用粗网格。本文使用自适应网格法模拟二维均质蒸汽辅助重力驱油过程,并将计算结果与国际上通用的油藏数值模拟软件CMG-STARS的计算结果进行详细的对比。对比结果表明,自适应网格法与CMG-STARS在计算精度上基本一致,而计算效率则明显优于CMG-STARS软件。自适应网格法在计算效率上的优越性,使得对大样本参数条件进行反复计算的周期大大缩短,从而能快速发现油藏开采中的问题,并给出优化策略。因而,自适应网格法在油藏数值模拟中具有巨大的工程实际意义。本文对自适应网格法中存在的若干难点问题进行了全面而深入的研究。首先,自适应网格法中的网格系统由细到粗分为若干层,我们需要由最精细网格的渗透率计算出各层粗网格的等效渗透率。对于空间无关联的均匀油藏内的流动,每个精细网格内的渗透率都是相同的或是相差不大的,网格变量的粗化可以通过算术平均或调和平均即可达到较高的精度要求。然而对于带有空间关联性的非均匀油藏,算术平均和调和平均在大部分情况下不再有效,如何快速有效地计算粗网格的等效渗透率,是发展复杂地层蒸汽热采稠油自适应网格法的重要科学问题之一。重整化方法就是计算等效渗透率的一种重要方法。对于带有空间关联性的多孔介质,重整化渗透率的概率分布,渗透率方差的标度特性,以及重整化之后的渗透率不动点的变化规律都是我们研究的重点。我们用相关长度来衡量多孔介质的空间关联性。与空间无关联多孔介质不同的是,带有空间关联性的多孔介质的渗透率方差的对数与重整化次数之间的比值不再是一条直线。在重整化过程中,当相关长度大于重整化尺度时,渗透率的方差将缓慢减少,并且当相关长度远远大于重整化尺度时,标度指数将趋于0。随着重整化次数的逐渐增加,当重整化尺度远远大于相关长度时,标度指数将趋于-1n4,与空间无关联多孔介质重整化过程的统计特性一致。同时,我们还研究了带有相关长度的多孔介质重整化之后的渗透率不动点。如果不同方向的相关长度相等,则不同方向的渗透率不动点完全相同,并与空间无关联多孔介质的渗透率不动点保持一致,即与相关长度的大小无关。当不同方向的相关长度不相等时,无论是理论结果还是数值模拟结果都说明,某一方向的相关长度的数值较大,则对应的渗透率不动点的数值也较大。为了能够准确快速的计算出三维空间自适应网格法中多层网格的等效渗透率,我们还在二维空间重整化方法(King,1987; Noetinger,1997)的基础上,给出了三维空间中标量形式以及张量形式的等效渗透率的重整化方法。数值计算的结果也表明了这两种重整化方法的计算结果都接近于作为参考值的有限差分的计算结果,特别是张量等效渗透率的计算结果比起标量的更加准确。其次,我们尝试利用自适应网格法来模拟分布有多种具有不同物理性质(包括骨架性质和孔隙性质)的储层岩石的复杂油藏问题。不同的储层岩石具有不同的相对渗透率曲线和表面张力曲线,各相流体的相饱和度在不同储层岩石的交界面处将发生空间间断,如果继续用各相饱和度作为网格系统的划分依据的话,将导致精细网格的数量大幅增加,影响计算的效率。在文中,我们从网格的各界面流量守恒性出发,并作如下假定：忽略毛管力；各相流量在粗网格中的变化是相当小的；同一个粗网格内同种岩石的各相饱和度是一致的。在这些假定的前提下,给出了二维空间多种岩石自适应网格法的粗网格细化以及细网格粗化的新准则,并取得了不错的数值模拟结果。最后,我们对基础的油水两相等温渗流的数值算法进行了初步的探讨。描述饱和度变化的连续方程,在数学上对应非线性双曲型方程,为了避免中心差分算法带来的物理上不真实的振荡,在通常的数值模拟中我们都采用了迎风格式。然而迎风方法只有一阶精度,会导致较明显的数值耗散。同时控制方程的离散,均采用五点差分格式,这在高维空间中的某些情况下将会导致严重的网格取向效应。虽然可以通过采用九点差分格式、PEBI网格等来减少网格取向效应,但是尚没有统一的数值算法可以有效避免各种情况下出现的网格取向效应。为了尽量减少数值耗散,降低网格取向效应,我们对基础的油水两相等温渗流的数值算法进行了初步的探讨,通过利用特征线方法分析Buckley-Leverett方程,提出了活塞式驱替和非活塞式驱替的严格数学划分条件。由于非活塞式驱替的饱和度在全场连续,本文建议对网格界面处相对渗透率的离散值不采用常用的在迎风方向取值的做法,而直接取相邻网格的算术平均值。数值计算结果表明,对非活塞式驱替而言,该算法不仅可以提高计算精度,而且可以减小网格取向效应。综上,本文的工作是从自适应网格法的计算结果出发,对自适应网格法中存在的若干难点问题进行研究,为后期的复杂油藏的数值模拟打好坚实的理论基础。相应的数值算例也验证了我们的理论。我们期望通过对油藏数值模拟中自适应网格法的基础问题的研究,能够给以后的工作带来实际的帮助。

全文目录

摘要  6-9
ABSTRACT  9-14
第1章绪论  14-22
  1.1 油藏数值模拟简介  14-16
  1.2 稠油开采的自适应网格法  16-20
  1.3 本文的工作  20-22
第2章蒸汽辅助重力驱油的自适应网格法  22-38
  2.1 蒸汽辅助重力驱油的数学模型  22-24
  2.2 控制方程的离散及数值计算方法  24-27
    2.2.1 控制方程的离散  24-26
    2.2.2 求解非线性方程组  26-27
  2.3 自适应网格法  27-29
  2.4 CMG-STARS软件介绍  29-32
  2.5 蒸汽辅助重力驱油的AMR与Stars计算结果对比  32-36
  2.6 非均质油藏的自适应网格法中若干难点问题  36-37
  2.7 本章小结  37-38
第3章渗透率粗化的重整化过程的理论分析  38-52
  3.1 渗透率粗化重整化方法的数学模型  38-42
    3.1.1 king的重整化方法  39-41
    3.1.2 改良的重整化方法  41-42
  3.2 空间无关联多孔介质渗透率粗化的重整化过程的理论分析  42-43
  3.3 空间关联多孔介质渗透率粗化的重整化过程的理论分析  43-50
    3.3.1 重整化渗透率的概率分布  46-47
    3.3.2 重整化渗透率方差的标度特性  47-48
    3.3.3 重整化渗透率的不动点  48-50
  3.4 本章小结  50-52
第4章三维空间渗透率粗化的重整化方法  52-58
  4.1 三维空间重整化方法的数学模型  52-55
    4.1.1 标量等效渗透率  52-53
    4.1.2 张量等效渗透率  53-55
  4.2 数值检验  55-57
  4.3 本章小结  57-58
第5章多种岩石油藏自适应网格法的初步研究  58-66
  5.1 多种岩石油藏多相流动的自适应网格法  58-61
  5.2 数值算例  61-64
  5.3 本章小结  64-66
第6章油水两相流非活塞式驱替的数值算法探讨  66-80
  6.1 油水两相流常用的数值算法及其中存在的问题  66-71
  6.2 基于Buckley-Leverett方程的进一步分析  71-73
  6.3 非活塞式驱替的数值算例  73-78
  6.4 本章小结  78-80
第7章结论及展望  80-84
  7.1 全文总结  80-81
  7.2 本文工作创新点  81-82
  7.3 未来工作展望  82-84
参考文献  84-94
致谢  94-96
攻读博士学位期间发表的论文  96

油藏数值模拟自适应网格法中若干问题研究

内容摘要

全文目录

相似论文