低渗透储层流固耦合渗流规律探究.pdfVIP

低渗透储层流固耦合渗流规律研究 刘建军 刘先贵 胡雅衽 张盛宗 (中国科学院涪流潍悻力学研究所，河北鹰坊嘶5∞7) 囊耍考虑低渗透油藏涪巍时启动压力梯度和低渗储层的漉同耦台特性，将渗流力学与弹塑 性力学相结合，建立起适合低渗透油藏的流固耦台渗漉数学攘型，并给出其数值解．在黑油 模型和弹塑性有甩元程序的基础上，编制了计算低渗透油冀流固耦台渗漉的计算软件。通过 数值摸拟和不考虑流固耦合时的计算结果相比，可以看出低渗透油藏中流固耦合效应是十分 明显的． 关悯怔渗透储层．多相渗漉．弹塑性变形．流固耦台 引言 油气藏的开发过程是一个多相流体渗流与油藏岩土变形动态耦合的过程。在我国低渗透 油田的开发过程中，人们发现低渗透油田流周耦台效应较中高渗透油田要强得多。这主要是 因为低渗透油层孔骧度、渗透率随有效应力的变化明显所造成的。另外低渗透油藏的渗流机 制有其特殊性，其渗漉衄线偏离迭西定律。由于低渗透油田的这些特点，使得低渗透油田的 流圃耦合理论不同于中高渗透油田。因此有必要对低渗透介质的漉固耦合渗流理论进行专门 的研究．本文在低渗透介质渗流规律和储层流固耦合特性基础上，建立适合低渗透油藏流固 耦合渗流的教学模型．给出其数值解，并进行了实例计算。 1 低渗透储层渡固耦合渗流特性的实验研究 1．1低渗透储屡漉固■合渗流特性 为了考察开发过程中压力变化对储层物性参数的影响，我们进行了压力变化对孔隙度、 渗透率影响的室内物理模拟实验研究。通过大量的实验发现： 1)随着有效应力的增加，岩芯的孔、渗参数均有不同程度的下降，即使有效应力提高． 其恢复程度也报小。 2)渗透率和孔隙度随有效压力的增加而下降的程度不同。孔隙度平均下降5％，而渗透 率平均下降78％．这说明多孔介质中孔敷通道的微小变化对流体流动能力的影响非常显著。 3)相对于基质岩芯，有效应力增加时裂缝性岩芯物性参数受到的伤害更大且更难恢复。 4)孔骧度、渗透率随有效压力的变化规律可用指数形式表示为 f=九。“．^=Ioe一∥ (1) 式中口．口是参数，仃-为有效应力。 2 低渗透油蠢麓固锅合渗流数学横型的建立 2．1带启动压力梯度(T阳)的低渗透油叠渗流的微分方程 1)气组分方程 玑陕h一^胛肛珊)+等他一“一肛册{+老。鲁忙+警]J∞ 2)油、水组分方程 ㈤ 可-I鲁他一风庐肛聃)卜老=鲁(鲁)(o，w) 3)辅助方程．(包括饱和度、毛管力、相对渗透率方程) 2．2流固，I台渗流油藏岩土变形场方程 将油藏岩土视为弹塑性介质，则变形场方程包括以下三组 1)平衡方程 o；．j+b591 (4) 7i+f。=o 2)几何方程 铲抛，Ⅵ．，) (5) 3)变形场本构方程采用弹塑性本构方程，本构方程的增量形式为 {l唧’}=【D，胁，} (6) 3低渗透油藏流固耦合渗流模型的数值解法 岩土变形场的求解采用有限元法，选用空间八节点等参数单元离散。弹塑性处理采用初 应力法。孔隙压力作为均布体积力处理。渗流模型的求解采用常规的IMPES差分解法，即隐 式求解压力方程，显式求解饱和度，渗流场和变形场通过(1)式耦合。 耦合模型求解的总体方案是显式迭代求解．渗流场的求解先于岩土变形场，采用有限差 分法计算渗漉场后，把所得孔隙压力增量传递给变形场．变形塌采用有限元方法求解，利用 所得的应力结果修正渗透系数和孔隙度，然后重新计算渗流场，如此循环迭代，直到结束。 4实倒计算 根据上面的数学模型及解法，我们编制了计算低渗透油藏流固耦台渗流的程序。利用编 髑的软