出砂储层流固耦合渗流数值模拟研究.pdfVIP

出砂储层流固耦合渗流数值模拟 王正茂1 阎存章2 刘明新1 (1．中国石油勘探开发研究院，北京100083； 2．中油股份公司勘探与生产分公司，北京100083) 摘要为了研究储层出砂后含砂流体的流动状况j准确预测其开采动态，运用有效应力 原理、质量守恒原理、虚位移原理和虚功等效原则等，建立了含砂流体流固耦合渗流的三维四 相单井黑油模型和考虑油井围岩弹塑性变形的岩土平衡方程；从基本定义出发，考虑体积应 变、骨架砂剥离、可动砂在孔隙表面和喉道沉积对物性参数的影响，导出了出砂储层数值模拟 所需的物性参数动态模型；探讨了模型的耦合求解方法；最后，通过对一实例进行模拟分析，证 明了该模型的正确性和有效性。该模型的建立对于正确认识和高效开发出砂储层有着重要意 义。 关键词 出砂，流固耦合，渗流，数值模拟 1 前言 出砂储层的开采过程是一个流固耦合渗流过程，不仅流体质点要在多孔介质中发生流动， 而且储层岩土由于载荷作用要产生变形，岩土质点要发生位移而产生运动；当储层岩土的变形 超过其极限，岩土骨架将遭到破坏，部分骨架砂将进人流体成为可动砂，这些可动砂在流体的 作用下要发生运移、沉降或者堵塞储层孔隙喉道等，反过来影响流体的渗流及储层岩土的变 形。由于在常规储层渗流机理的研究中没有考虑砂粒对流体渗流和骨架变形的影响，因此，常 规储层渗流数学模型显然不能描述出砂储层流体渗流的实际状况。为了正确认识储层出砂后 含砂流体的流动状况．准确预测其开采动态，必须对出砂储层渗流机理进行研究。 2 出砂储层流固耦合渗流单井黑油模型 2．1油藏岩土变形假设 岩土发生弹性或弹塑性小变形；岩土颗粒不可压缩，储层多孔介质的孔隙可压缩；岩土质 点在开采过程中要发生位移。 2．2油藏流体渗流假设 油藏流体有油、气、水三相，考虑可动砂相运移的影响，流体可压缩；等温渗流；介质被油气 水和砂的混合液饱和；油藏流体服从达西渗流定律，可动砂与油相的运动速度相同[13；水相与 油相不发生组分交换；油藏中烃类只含有油、气两个组分；在油藏状况下，油、气两种组分可以 形成油气两相，油组分完全存在于油相内，气组分可以自由气的形式存在于气相中，也可以溶 解气的形式存在于油相和水相中；考虑毛管压力和重力的影响。 2．3 出砂储层渗流数学模型 在出砂储层的开采过程中，设渗流系统中有油、气、水和可动砂四相。根据基本假设、质量 守恒原理和流固耦合渗流运动方程建立出砂储层流固耦合渗流三维四相单井黑油模型。 (1)固相骨架连续性方程 箬刈刊鲁一尝 式中：9为孔隙度，￡。为体积应变，B为骨架砂密度，R，为从单位体积储层骨架上砂的质量剥离 速率。 (2)油相连续性方程 (2) 毫(警)+V·(警z)一V·[甓cV卜以gVD，一 (3)气相连续性方程 旦f璺监+堡墨!堇!g!+堕璺竺兰芝鱼1 ‘ af＼B。。 B。 B。 ／ +V·(警z+％竽z十％笋z) 一V·[紫cVP。一P。gVD，]=。 (3) (4)水相连续性方程 丢(警)+V·(警z) 一V·[疑cVP。咄gVD，一 (4) 式中：下标D、g、铷为油气水相，B为体积系数，|D为密度，s为饱和度茁为岩土质点运动速度， C，为流体在流体和可动砂混合液中的体积浓度，K为渗透率，p为粘度，P为压力，D为高程， R和为溶解气油比，R扣为溶解气水比。 (5)可动砂连续性方程 毫(GP)+V·(e最)～V