变渗流系统物性油藏数值模拟软件开发.pdfVIP

变渗流系统物性油藏数值模拟软件开发 鞠斌山 樊太亮 中国地质大学 （北京）能源学院，北京，100083 摘要：为研究油田开发过程中渗流系统物性 （①孔隙度和渗透率；②牛顿流体的流动性质和非牛顿流体的流 变学性质；③岩石的润湿性和相对渗透率）的变化及其对渗流过程和油田开发动态的影响，综合应用理论分析、物 理实验模拟和矿场资料统计与回归方法，研究了渗流系统物性变化的机理和规律，建立了渗流系统物性变化数学 模型。基于连续介质假设，建立了渗流系统物性变化的油藏渗流数学模型，应用有限差分方法对该渗流数学模型进 行了数值求解。采用FORTRAN90和VC++两种计算机语言，开发了渗流系统物性变化油藏数值模拟器。利用该 模拟器研究了渗流系统物性变化对剩余油和开发效果的影响。结果表明：开发过程中油层变形、微粒运移引起的非 均质性加剧、原油粘度升高对油田开发是不利的，水驱引起的岩石弱亲水向强亲水性转化是有利的。本研究对提高 剩余油和开发指标的预测精度具有重要的理论和应用价值。 关键词：油藏数值模拟;孔隙度和渗透率;非牛顿流体;数学模型;润湿性 · · 876 变渗流系统物性油藏数值模拟软件开发 目前的商业软件没有考虑上述渗流系统物 气相渗流方程 性的变化。对变储层物性的油藏进行数值模拟 时，不少数值模拟工作者[29~31]在进行历史拟合 渤I麓g喊+等gm泡卜 和预测时，仍然使用商业模拟器，采用分段模 拟的方法，即在不同的含水阶段施用不同的储 一叩导(≯·／B。+乒·s水，／Bo) (3) 层物性参数。该方法将导致各阶段之间的衔接 式中 出现跳跃，与油田实际不符。 或一仇+托·z (4) 本文通过理论研究油藏孑L渗参数，润湿性 瓯一p。+‰·z—A+户。+‰·z(5) 和流体性质变化的机理及影响因素；借助室内 嚷一以+yg·z一户。+p印+坛·z(6) 物理模拟实验研究油藏渗流系统的物性变化规 三相饱和度满足 律和幅度；在此基础上，运用现代渗流理论和 S+S。+S。一1 (7) 应用数学方法建立一套系统的描述开发过程中 旷压力，MPa；S_饱和度，f；K一地层渗透 油藏渗流系统物性变化的油藏渗流模型，运用 率，肛m2；K，一相对渗透率，f；弘一粘度，mPa 偏微分方程的数值计算技术对所建数学模型求 ·s；R一溶解气油比，f；z-标高，m；卜地 解，应用VC++和FORTRAN两门计算机语 层孔隙度，f；p。一毛管压力，MPa；y_流体 言，开发充分反映油藏渗流系统物性变化的油 比重，q一源汇项，m3／s；即从单位体积岩石中 藏数值模拟器。运用开发的油藏数值模拟器对 采出(为负)或向单位体积岩石中注人(为正) 油田实例进行模拟，定量研究油藏渗流系统物 的地面条件下的流体量；r单位换算系数，取 性变化规律及其对剩余油形成和开发动态的影 决于单位制；出扩散度；gr口d一梯度；下标的 响，指导油田开发。 意义：议广-水，旷油，g一气。 2．3渗流系统物性变化机理与规律 2数学模型的建立 2．3．1油藏变形机理及孔渗参数变化规律 2．1假设条件 根据力学理论，固体受到的应力变化时， ①三维三相等温渗流； 其必然发生变形，变形的固体被称为变形介质， ②非均质可压缩变形孔隙介质； 具有如此性质的油藏被称为变形介质油藏。根 ③固体微粒参与渗流；