低渗透油藏非线性渗流机理及数值模拟方法研究

低渗透油藏非线性渗流机理及数值模拟方法研究姜瑞忠中国石油大学（华东）中国石油大学（华东）国际研讨会引言随着世界石油资源形势的日益严峻，低渗透油藏开发已成为石油资源产量接替的重要组成部分。而目前的低渗透油藏开发理论不能有效地指导实际生产，微观渗流机理以及数值模拟方法的研究还存在着许多争论及不足之处。本文以国家重大科技专项之子课题“低渗、特低渗油气田经济开发关键技术”为依托，主要从非线性渗流机理、非线性渗流描述方法以及数值模拟方法三个方面进行深入探讨。引言中国石油大学（华东）目录第一部分低渗透油藏非线性渗流微观机理探讨第二部分低渗透油藏非线性渗流描述方法研究第三部分低渗透油藏非线性渗流数值模拟方法研究中国石油大学（华东）1234一、非线性渗流微观机理探讨1.1低渗透油藏渗流特征渗流规律呈非线性特征低渗透多孔介质中流体流动的横截面积是可变的低渗透多孔介质的渗透率并非常数低渗透油藏中渗流存在启动压力梯度中国石油大学（华东）一、非线性渗流微观机理探讨从理论上讲，流体在多孔介质内流动时，均不同程度地存在有启动压力梯度。这是因为任何流体都存在一定的屈服应力值，只有驱替压力能够克服这个屈服应力值时流体才开始流动。而且因流体与多孔介质的相互作用导致边界层的出现是必然的。边界层的出现导致流体在渗流孔道中出现一定的分布，越靠近多孔介质固体界面流体越致密，许多学者对此进行了相应的研究，并提出了渗流流体的概念。随着驱替压力梯度的增大，参与流动的流体的份额就越大，动用的边界层就越大，所应克服的屈服应力值就越大，因此启动压力梯度也就越大。1.2低渗透油藏启动压力梯度探讨中国石油大学（华东）一、非线性渗流微观机理探讨很多学者的研究证明，启动压力梯度与渗透率呈反比，渗透率越低，启动压力梯度越大。低渗透油藏孔径很小，原油边界层的影响显著，导致宏观表现出来的启动压力梯度就很明显。而对于中高渗透油藏，由于孔道半径比较大，原油边界层的影响微弱，尤其对于稀油油藏，流体的屈服应力值很小，导致用一般的实验手段不易测到启动压力梯度。因此在实际生产分析时就忽略了启动压力梯度的影响。1.2低渗透油藏启动压力梯度探讨中国石油大学（华东）一、非线性渗流微观机理探讨值得一提的是，对于部分低渗透油藏，因非均质状况十分严重，渗流孔径跨度较大，这时会出现实验中也没有观测到启动压力梯度。这是因为虽然油藏整体评价属于低渗透油藏，但存在着小部分渗流孔径属于中高渗透油层范畴，一般的实验手段也没有观测到启动压力梯度。因此实际生产分析时部分学者没有考虑启动压力梯度也是合理的。但对于特低渗透油藏忽略启动压力梯度就会产生较大的误差。1.2低渗透油藏启动压力梯度探讨中国石油大学（华东）一、非线性渗流微观机理探讨由于低渗岩心的孔隙系统基本上是由小孔道组成的，在油、水流动时，每个孔道都有自己的启动压力梯度，只有驱动压力梯度大于某孔道的启动压力梯度时，该孔道中的油、水才开始流动。随着驱动压力梯度的不断提高，就会有更多的孔道加入到流动的行列，岩心的渗透性能也随之增强，渗透率变大。1.3低渗透油藏渗流规律呈现非线性特征探讨1.3.1低渗透多孔介质的渗透率并非常数中国石油大学（华东）一、非线性渗流微观机理探讨借助大家公认的毛细管模型来分析低渗透油藏的渗流孔道特征尺寸。根据孔隙孔道半径与渗透率之间的关系式：可以计算出三类低渗透油田的平均毛管半径如下：第一类为一般低渗透油田，油层平均渗透率为。平均毛管半径为（1.11—2.47）。第二类为特低渗透油田，油层平均渗透率为。平均毛管半径为（0.37—1.11）。第三类为超低渗透油田，其油层平均渗透率为。平均毛管半径（0.11—0.37）。m0.35rK32(10.150)10mm32(1.110.0)10m32(0.11.0)10mm1.3.2低渗透油藏存在微尺度流动效应1.3低渗透油藏渗流规律呈现非线性特征探讨中国石油大学（华东）一、非线性渗流微观机理探讨通过上述计算可以看出，低渗透油藏渗流通道的特征尺寸已进入微米级别，根据中国科学院陶然、李战华等专家的研究结论，微米级别已属于微尺度流动。但显然远大于分子的特征尺寸，因此虽然连续介质场的假设依然成立，但由于流动尺度的微小，许多在宏观流动中被忽略的或一些影响较小的因素，此时成为主要的影响因素，从而导致流动规律不同于宏观的流动规律，渗流曲线出现非线性段。也就是说低渗透油藏非线性渗流出现的原因从某种程度上说要归因于微尺度流动效应。1.3.2低渗透油藏存在微尺度流动效应中国石油大学（华东）1、尺度效应：随着尺度减小到微米级时，体积力的作用下降；表面力的作用却上升，超过了体积力；随着尺度的减小，粘性力相对作用增强，惯性力作用变小，越靠近壁面这种规律越明显。流变学的结论，因低渗透油藏特征尺度很小，使得速度的空间导数较大，无论是原油还是地层水在渗流时，均会表现出非牛顿的特性，原油的非牛顿特性更强一些。微尺度效应2、表面效应：随着尺度的减小，渗流通道的表面积与体积之比大大增加，流体处在表面上的分子相对于内部分子所占的比例增大，这在很大程度上影响了渗流通道表面处的质量、动量和能量的传输，表面力将会起主要作用。由于面体比的变化使固体表面的界面效应明显。一、非线性渗流微观机理探讨中国石油大学（华东）一、非线性渗流微观机理探讨笔者经过多年的研究认为，低渗透油藏出现非线性渗流的主要原因就包括两方面：1、低渗透多孔介质的渗透率并非常数，而是随驱替压力梯度变化而变化2、低渗透油藏渗流存在微尺度流动效应。因低渗透油藏孔喉特征尺寸细小，处在连续流动区域以及分子效应区的交界处，虽然连续介质场的假设仍能成立，但由于流动尺度的微小，许多在宏观流动中被忽略的因素，此时成为主要的影响因素，吸附边界层的存在更加加剧了微尺度流动效应，从而导致流动规律不同于中高渗油藏的流动规律，渗流曲线上存在着启动压力梯度以及非线性段，因此达西定律不再适用，需进行修正。中国石油大学（华东）目录第一部分低渗透油藏非线性渗流微观机理探讨第二部分低渗透油藏非线性渗流描述方法研究第三部分低渗透油藏非线性渗流数值模拟方法研究中国石油大学（华东）二、非线性渗流描述方法研究前文对微观机理的研究指出，达西定律已不能满足对低渗透油藏渗流规律的精确描述，需要进行修正。但因启动压力梯度随压力梯度变化的特性、渗透率并非常数的特性、存在微尺度流动效应以及边界层，使得要得到能够精确描述真实渗流规律的运动方程变得十分困难。这就是目前绝大多数学者倾向于用唯象的方法来得到运动方程的原因。目前很多学者对低渗透油藏的非线性渗流描述方法进行了研究，也提出了很多模型，归纳总结起来有如下三类：中国石油大学（华东）分段模型分段模型将渗流曲线分为非线性段和线性段，并用不同的函数进行拟合，但在应用过程中需要对临界点进行判断，应用难度较大。连续模型连续模型用统一的函数对渗流曲线进行描述，克服了常规连续模型不能描述渗流存在最小启动压力梯度的问题。2.1非线性渗流描述方法拟启动压力梯度模型拟启动压力梯度模型应用最广泛，但该方法用直线代替渗流曲线中的非线性段，降低了渗流能力，并且抬高了流体是否流动的界限，计算结果存在着误差。非线性渗流描述方法共分为三类：二、非线性渗流描述方法研究中国石油大学（华东）（1）由分段函数模型向连续模型函数发展非线性渗流模型研究特点（2）唯象法模型为主要的数学模型2.1非线性渗流描述方法二、非线性渗流描述方法研究中国石油大学（华东）2.2非线性渗流模型的选取prppppkkM11||pMabppppD1(0)(0,)cckkbGkkk式中，vp为相渗流速度，m/s；k为绝对渗透率，；krp为相相对渗透率；为相粘度，；为p相的流动势；M为非线性渗流修正系数；D为油藏海拔深度，m；G为拟启动压力梯度，MPa/m；a为影响非线性渗流曲线段的影响因子，a0，无量纲；b为渗流曲线中拟启动压力梯度的倒数，1/MPa；a、b与地层渗透率有关，由拟合不同渗透率岩芯渗流实验曲线确定。p二、非线性渗流描述方法研究中国石油大学（华东）2.2非线性渗流模型的选取在不同a，b取值下，渗流曲线如下图所示：二、非线性渗流描述方法研究运动方程中，当a=0时，模型即为拟启动压力梯度模型。当b取无穷大时，模型即为常规达西渗流模型。因此，本模型为通用模型，同样具备模拟拟启动压力梯度模型和达西渗流模型的能力。中国石油大学（华东）特别指出，单相达西定理推广到两相时，引入了相对渗透率的概念，其基本假设就是“渠道流态”，因此引入启动压力梯度以及非线性相关参数时应该取单项时的测量值。部分学者曾研究过两相启动压力梯度，认为两相启动压力梯度值要远大于单项启动压力梯度值，且给出的原因是贾敏效应以及不利毛管力的影响，而实际上相对渗透率已经将贾敏效应以及不利毛管力的影响考虑在内了，若启动压力梯度项及非线性相关参数仍取两相渗流时的测量值，就会出现对贾敏效应以及不利毛管力的影响的重复计算，这是不合理的。二、非线性渗流描述方法研究中国石油大学（华东）目录第一部分低渗透油藏非线性渗流微观机理探讨第二部分低渗透油藏非线性渗流描述方法研究第三部分低渗透油藏非线性渗流数值模拟方法研究中国石油大学（华东）三、非线性渗流数模方法研究目前比较成熟的商业软件都是基于Darcy渗流规律的，没有考虑第一部分提到的低渗透储层特有的特征。因此努力开发一套基于非线性渗流规律的油藏数值模拟软件，为低渗透油田的开发提供研究手段和工具并用于指导低渗透油田的高效开发，提高低渗透油藏开发的主动性和科学性，在目前国际石油形势日趋严峻的情况下具有十分重要的意义。与中高渗透油藏相比低渗透油藏数值模拟需要考虑的几个关键问题包括：渗流曲线上存在启动压力梯度和非线性段、压敏效应。中国石油大学（华东）3.1.1启动压力梯度及非线性段处理三、非线性渗流数模方法研究)||()||11(pbappkppbakv本文仍采用前文提到的二参数连续模型的方法：描述启动压力梯度项及非线性段流动分三段：A不流动段B非线性段C直线段其中，非线性段上，随压力的增大启动压力梯度增大。中国石油大学（华东）采用这种处理方法，可以根据不同油藏状况对各个网格赋不同的a、b值，以模拟启动压力梯度的非均质情况，克服了拟启动压力梯度梯度模型启动压力梯度是定值的弊端，而且X、Y、Z三个方向因渗透率存在着差异，在三个方向上a、b值可选取不同的值。这样不仅实现了启动压力梯度的非均质性，而且实现了启动压力梯度方向性差异，对于低渗透油藏数值模拟技术是一大突破。三、非线性渗流数模方法研究中国石油大学（华东）(1)常系数法3.1.2压敏效应三、非线性渗流数模方法研究0101100100101100expexpexpexpppppkkpppppppkkk(2)变系数法压实曲线弹性曲线压力乘子01010101expexpppkkppk中国石油大学（华东）油气水三维三相连续性方程的建立油气水及岩石状态方程的建立研究最符合低渗储层实际的非线性渗流描述方法，建立非达西渗流的运动方程。对数学模型进行离散化，建立低渗储层差分格式，形成数值模型以非线性渗流为特点的线性方程组解法研究建立概念模型，开展低渗与中高渗油藏对比分析研究编制低渗透油藏非线性渗流数值模拟软件非线性渗流数模方法研究主要内容三、非线性渗流数模方法研究中国石油大学（华东）（1）油藏中的流体流动是等温渗流；（2）油藏中最多只有油、气、水三相，油相和水相的渗流均遵循非线性渗流定律；（3）油藏烃类只含油气两个组分。在油藏状况下，油、气两种组分可能形成油、气两相。油组分完全存在于油相中，气组分既能以自由气的方式存在，又可溶解于油相之中。地层中的油相应为油组分和气组分的某种组合；（4）油藏中气体的溶解和逸出是瞬间完成的；（5）油、水之间不互溶；（6）岩石微可压缩，且考虑渗流过程