《低渗透油藏非线性渗流机理及数值模拟方法研究PPT幻灯片课件》由会员分享，可在线阅读，更多相关《低渗透油藏非线性渗流机理及数值模拟方法研究PPT幻灯片课件（53页珍藏版）》请在金锄头文库上搜索。
1、低渗透油藏非线性渗流机理及数值模拟方法研究 姜瑞忠中国石油大学 华东 1 中国石油大学 华东 国际研讨会引言 随着世界石油资源形势的日益严峻 低渗透油藏开发已成为石油资源产量接替的重要组成部分 而目前的低渗透油藏开发理论不能有效地指导实际生产 微观渗流机理以及数值模拟方法的研究还存在着许多争论及不足之处 本文以国家重大科技专项之子课题 低渗 特低渗油气田经济开发关键技术 为依托 主要从非线性渗流机理 非线性渗流描述方法以及数值模拟方法三个方面进行深入探讨 引言 2 中国石油大学 华东 目录 3 中国石油大学 华东 一 非线性渗流微观机理探讨 1 1低渗透油藏渗流特征 4 中国石油大学 华东 一
2、 非线性渗流微观机理探讨 从理论上讲 流体在多孔介质内流动时 均不同程度地存在有启动压力梯度 这是因为任何流体都存在一定的屈服应力值 只有驱替压力能够克服这个屈服应力值时流体才开始流动 而且因流体与多孔介质的相互作用导致边界层的出现是必然的 边界层的出现导致流体在渗流孔道中出现一定的分布 越靠近多孔介质固体界面流体越致密 许多学者对此进行了相应的研究 并提出了渗流流体的概念 随着驱替压力梯度的增大 参与流动的流体的份额就越大 动用的边界层就越大 所应克服的屈服应力值就越大 因此启动压力梯度也就越大 1 2低渗透油藏启动压力梯度探讨 5 中国石油大学 华东 一 非线性渗流微观机理探讨 很多学者的
3、研究证明 启动压力梯度与渗透率呈反比 渗透率越低 启动压力梯度越大 低渗透油藏孔径很小 原油边界层的影响显著 导致宏观表现出来的启动压力梯度就很明显 而对于中高渗透油藏 由于孔道半径比较大 原油边界层的影响微弱 尤其对于稀油油藏 流体的屈服应力值很小 导致用一般的实验手段不易测到启动压力梯度 因此在实际生产分析时就忽略了启动压力梯度的影响 1 2低渗透油藏启动压力梯度探讨 6 中国石油大学 华东 一 非线性渗流微观机理探讨 值得一提的是 对于部分低渗透油藏 因非均质状况十分严重 渗流孔径跨度较大 这时会出现实验中也没有观测到启动压力梯度 这是因为虽然油藏整体评价属于低渗透油藏 但存在着小部分渗
4、流孔径属于中高渗透油层范畴 一般的实验手段也没有观测到启动压力梯度 因此实际生产分析时部分学者没有考虑启动压力梯度也是合理的 但对于特低渗透油藏忽略启动压力梯度就会产生较大的误差 1 2低渗透油藏启动压力梯度探讨 7 中国石油大学 华东 一 非线性渗流微观机理探讨 由于低渗岩心的孔隙系统基本上是由小孔道组成的 在油 水流动时 每个孔道都有自己的启动压力梯度 只有驱动压力梯度大于某孔道的启动压力梯度时 该孔道中的油 水才开始流动 随着驱动压力梯度的不断提高 就会有更多的孔道加入到流动的行列 岩心的渗透性能也随之增强 渗透率变大 1 3低渗透油藏渗流规律呈现非线性特征探讨 1 3 1低渗透多孔介质
5、的渗透率并非常数 8 中国石油大学 华东 一 非线性渗流微观机理探讨 借助大家公认的毛细管模型来分析低渗透油藏的渗流孔道特征尺寸 根据孔隙孔道半径与渗透率之间的关系式 可以计算出三类低渗透油田的平均毛管半径如下 第一类为一般低渗透油田 油层平均渗透率为 平均毛管半径为 1 11 2 47 第二类为特低渗透油田 油层平均渗透率为 平均毛管半径为 0 37 1 11 第三类为超低渗透油田 其油层平均渗透率为 平均毛管半径 0 11 0 37 1 3 2低渗透油藏存在微尺度流动效应 1 3低渗透油藏渗流规律呈现非线性特征探讨 9 中国石油大学 华东 一 非线性渗流微观机理探讨 通过上述计算可以看出
6、低渗透油藏渗流通道的特征尺寸已进入微米级别 根据中国科学院陶然 李战华等专家的研究结论 微米级别已属于微尺度流动 但显然远大于分子的特征尺寸 因此虽然连续介质场的假设依然成立 但由于流动尺度的微小 许多在宏观流动中被忽略的或一些影响较小的因素 此时成为主要的影响因素 从而导致流动规律不同于宏观的流动规律 渗流曲线出现非线性段 也就是说低渗透油藏非线性渗流出现的原因从某种程度上说要归因于微尺度流动效应 1 3 2低渗透油藏存在微尺度流动效应 10 中国石油大学 华东 2 表面效应 随着尺度的减小 渗流通道的表面积与体积之比大大增加 流体处在表面上的分子相对于内部分子所占的比例增大 这在很大程度上
7、影响了渗流通道表面处的质量 动量和能量的传输 表面力将会起主要作用 由于面体比的变化使固体表面的界面效应明显 一 非线性渗流微观机理探讨 11 中国石油大学 华东 一 非线性渗流微观机理探讨 笔者经过多年的研究认为 低渗透油藏出现非线性渗流的主要原因就包括两方面 1 低渗透多孔介质的渗透率并非常数 而是随驱替压力梯度变化而变化2 低渗透油藏渗流存在微尺度流动效应 因低渗透油藏孔喉特征尺寸细小 处在连续流动区域以及分子效应区的交界处 虽然连续介质场的假设仍能成立 但由于流动尺度的微小 许多在宏观流动中被忽略的因素 此时成为主要的影响因素 吸附边界层的存在更加加剧了微尺度流动效应 从而导致流动规律
8、不同于中高渗油藏的流动规律 渗流曲线上存在着启动压力梯度以及非线性段 因此达西定律不再适用 需进行修正 12 中国石油大学 华东 目录 13 中国石油大学 华东 二 非线性渗流描述方法研究 前文对微观机理的研究指出 达西定律已不能满足对低渗透油藏渗流规律的精确描述 需要进行修正 但因启动压力梯度随压力梯度变化的特性 渗透率并非常数的特性 存在微尺度流动效应以及边界层 使得要得到能够精确描述真实渗流规律的运动方程变得十分困难 这就是目前绝大多数学者倾向于用唯象的方法来得到运动方程的原因 目前很多学者对低渗透油藏的非线性渗流描述方法进行了研究 也提出了很多模型 归纳总结起来有如下三类 14 中国石
9、油大学 华东 连续模型用统一的函数对渗流曲线进行描述 克服了常规连续模型不能描述渗流存在最小启动压力梯度的问题 2 1非线性渗流描述方法 非线性渗流描述方法共分为三类 二 非线性渗流描述方法研究 15 中国石油大学 华东 非线性渗流模型研究特点 2 1非线性渗流描述方法 二 非线性渗流描述方法研究 16 中国石油大学 华东 2 2非线性渗流模型的选取 式中 vp为相渗流速度 m s k为绝对渗透率 krp为相相对渗透率 为相粘度 为p相的流动势 M为非线性渗流修正系数 D为油藏海拔深度 m G为拟启动压力梯度 MPa m a为影响非线性渗流曲线段的影响因子 a 0 无量纲 b为渗流曲线中拟启动
10、压力梯度的倒数 1 MPa a b与地层渗透率有关 由拟合不同渗透率岩芯渗流实验曲线确定 二 非线性渗流描述方法研究 17 中国石油大学 华东 2 2非线性渗流模型的选取 在不同a b取值下 渗流曲线如下图所示 二 非线性渗流描述方法研究 运动方程中 当a 0时 模型即为拟启动压力梯度模型 当b取无穷大时 模型即为常规达西渗流模型 因此 本模型为通用模型 同样具备模拟拟启动压力梯度模型和达西渗流模型的能力 18 中国石油大学 华东 特别指出 单相达西定理推广到两相时 引入了相对渗透率的概念 其基本假设就是 渠道流态 因此引入启动压力梯度以及非线性相关参数时应该取单项时的测量值 部分学者曾研究过
11、两相启动压力梯度 认为两相启动压力梯度值要远大于单项启动压力梯度值 且给出的原因是贾敏效应以及不利毛管力的影响 而实际上相对渗透率已经将贾敏效应以及不利毛管力的影响考虑在内了 若启动压力梯度项及非线性相关参数仍取两相渗流时的测量值 就会出现对贾敏效应以及不利毛管力的影响的重复计算 这是不合理的 二 非线性渗流描述方法研究 19 中国石油大学 华东 目录 20 中国石油大学 华东 三 非线性渗流数模方法研究 目前比较成熟的商业软件都是基于Darcy渗流规律的 没有考虑第一部分提到的低渗透储层特有的特征 因此努力开发一套基于非线性渗流规律的油藏数值模拟软件 为低渗透油田的开发提供研究手段和工具并用
12、于指导低渗透油田的高效开发 提高低渗透油藏开发的主动性和科学性 在目前国际石油形势日趋严峻的情况下具有十分重要的意义 与中高渗透油藏相比低渗透油藏数值模拟需要考虑的几个关键问题包括 渗流曲线上存在启动压力梯度和非线性段 压敏效应 21 中国石油大学 华东 3 1 1启动压力梯度及非线性段处理 三 非线性渗流数模方法研究 本文仍采用前文提到的二参数连续模型的方法 描述启动压力梯度项及非线性段 流动分三段 A不流动段B非线性段C直线段其中 非线性段上 随压力的增大启动压力梯度增大 22 中国石油大学 华东 采用这种处理方法 可以根据不同油藏状况对各个网格赋不同的a b值 以模拟启动压力梯度的非均质
13、情况 克服了拟启动压力梯度梯度模型启动压力梯度是定值的弊端 而且X Y Z三个方向因渗透率存在着差异 在三个方向上a b值可选取不同的值 这样不仅实现了启动压力梯度的非均质性 而且实现了启动压力梯度方向性差异 对于低渗透油藏数值模拟技术是一大突破 三 非线性渗流数模方法研究 23 中国石油大学 华东 1 常系数法 3 1 2压敏效应 三 非线性渗流数模方法研究 2 变系数法 24 中国石油大学 华东 油气水三维三相连续性方程的建立 油气水及岩石状态方程的建立 研究最符合低渗储层实际的非线性渗流描述方法 建立非达西渗流的运动方程 对数学模型进行离散化 建立低渗储层差分格式 形成数值模型 以非线性
14、渗流为特点的线性方程组解法研究 建立概念模型 开展低渗与中高渗油藏对比分析研究 编制低渗透油藏非线性渗流数值模拟软件 非线性渗流数模方法研究主要内容 三 非线性渗流数模方法研究 25 中国石油大学 华东 1 油藏中的流体流动是等温渗流 2 油藏中最多只有油 气 水三相 油相和水相的渗流均遵循非线性渗流定律 3 油藏烃类只含油气两个组分 在油藏状况下 油 气两种组分可能形成油 气两相 油组分完全存在于油相中 气组分既能以自由气的方式存在 又可溶解于油相之中 地层中的油相应为油组分和气组分的某种组合 4 油藏中气体的溶解和逸出是瞬间完成的 5 油 水之间不互溶 6 岩石微可压缩 且考虑渗流过程中重
15、力 毛管力的影响 7 地层渗透率是压力梯度的函数 3 1基本假设 三 非线性渗流数模方法研究 26 中国石油大学 华东 油组分方程 水组分方程 气组分方程 3 2数学模型 三 非线性渗流数模方法研究 27 中国石油大学 华东 外边界条件 外边界可以有3种形式 定压外边界 定流量外边界和混合外边界 内边界条件 定产量 或定注水量 定井底流压 3 3定解条件 三 非线性渗流数模方法研究 28 中国石油大学 华东 3 4非线性模型差分离散化方法 采用块中心七点有限差分方法 与常规黑油模型不同的是油 水相方程中含一项非线性渗流修正系数 修正系数的大小与各相的压力梯度有关 取值范围为 0 1 因此 将该
16、修正系数作为绝对渗透率的修正因子 不仅理论上可以解释清楚 而且能够保证求解方程中系数矩阵的稳定 但是 考虑非线性渗流的油藏数值模拟器必须对以下两方面进行特殊处理 3 4 1非线性渗流修正系数3 4 2井 网格流动方程 三 非线性渗流数模方法研究 29 中国石油大学 华东 3 4 1非线性渗流修正系数 从实用性考虑 本模型必须适应油藏非均质的情况 因此 必须求得在交界面处非线性渗流修正系数的值 假设离散网格系统由六面体组成 网格的x y z方向编号分别为i j k 以x方向为例 根据质量守恒原理 交界面处的非线性渗流修正系数 同理可得有类似的表达式 三 非线性渗流数模方法研究 30 中国石油大学 华东 3 4 2井 网格流动方程 为了保证所提数学模型的适定性 井 网格压力方程中也必须考虑非线性渗流的影响 井所在的网格在数值模拟中处理为等效的平面径向流 考虑一个时间步内井与其所在网格间的稳定渗流 设等效半径为 圆形油藏内为均匀介质 设该时间步内网格 井间p相的稳定流率为Qp 则任意半径处的流速为 式中 当 时 对应于汇项 时 对应于源项 三 非线性渗流数模方法研究 31 中国石油大学 华东
