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

《低渗透油藏非线性渗流机理及数值模拟方法研究》由会员分享，可在线阅读，更多相关《低渗透油藏非线性渗流机理及数值模拟方法研究（52页珍藏版）》请在金锄头文库上搜索。

1、国际研讨会引言国际研讨会引言 随着世界石油资源形势的日益严峻，低渗透油藏随着世界石油资源形势的日益严峻，低渗透油藏开发已成为石油资源产量接替的重要组成部分。而目开发已成为石油资源产量接替的重要组成部分。而目前的低渗透油藏开发理论不能有效地指导实际生产，前的低渗透油藏开发理论不能有效地指导实际生产，微观渗流机理以及数值模拟方法的研究微观渗流机理以及数值模拟方法的研究还存在着许多还存在着许多争争论及论及不足之处不足之处。 本文本文以以国家重大科技专项国家重大科技专项之子课题之子课题“低渗、特低渗低渗、特低渗油气田经济开发关键技术油气田经济开发关键技术 ” ”为依托，为依托，主要从主要从非线性渗非线

2、性渗流机理流机理、非线性渗流描述方法非线性渗流描述方法以及以及数值模拟数值模拟方法方法三三个方个方面面进行深入探讨。进行深入探讨。 引言引言中国石油大学（华东）中国石油大学（华东）第1页，共52页。目录目录第一第一第一第一部分部分部分部分低渗透油藏非线性渗流微观机理探讨低渗透油藏非线性渗流微观机理探讨第二第二第二第二部分部分部分部分低渗透油藏非线性渗流描述方法研究低渗透油藏非线性渗流描述方法研究第三第三第三第三部分部分部分部分低渗透油藏非线性渗流数值模拟方法研究低渗透油藏非线性渗流数值模拟方法研究中国石油大学（华东）中国石油大学（华东）第2页，共52页。1234一、非线性渗流微观机理探讨一、非

3、线性渗流微观机理探讨1.1 1.1 1.1 1.1 低渗透油藏渗流特征低渗透油藏渗流特征低渗透油藏渗流特征低渗透油藏渗流特征渗流规律渗流规律呈非线性呈非线性特征特征低渗透多低渗透多孔介质中孔介质中流体流动流体流动的横截面的横截面积是可变积是可变的的低渗透多低渗透多孔介质的孔介质的渗透率并渗透率并非常数非常数低渗透油低渗透油藏中渗流藏中渗流存在启动存在启动压力梯度压力梯度中国石油大学（华东）中国石油大学（华东）第3页，共52页。一、非线性渗流微观机理探讨一、非线性渗流微观机理探讨v 从理论上讲，流体在多孔介质内流动时，均不同程度地从理论上讲，流体在多孔介质内流动时，均不同程度地存在有启动压力梯度

4、。这是因为任何流体都存在一定的存在有启动压力梯度。这是因为任何流体都存在一定的屈服屈服应力值应力值，只有驱替压力能够克服这个屈服应力值时流体才开始，只有驱替压力能够克服这个屈服应力值时流体才开始流动。而且因流体与多孔介质的相互作用导致流动。而且因流体与多孔介质的相互作用导致边界层边界层的出现是的出现是必然的。边界层的出现导致流体在渗流孔道中出现一定的分布，必然的。边界层的出现导致流体在渗流孔道中出现一定的分布，越靠近多孔介质固体界面流体越致密，许多学者对此进行了相应越靠近多孔介质固体界面流体越致密，许多学者对此进行了相应的研究，并提出了渗流流体的概念。随着驱替压力梯度的增大，的研究，并提出了渗

5、流流体的概念。随着驱替压力梯度的增大，参与流动的流体的份额就越大，动用的边界层就越大，所应克服参与流动的流体的份额就越大，动用的边界层就越大，所应克服的屈服应力值就越大，因此启动压力梯度也就越大。的屈服应力值就越大，因此启动压力梯度也就越大。1.2 1.2 1.2 1.2 低渗透油藏启动压力梯度探讨低渗透油藏启动压力梯度探讨低渗透油藏启动压力梯度探讨低渗透油藏启动压力梯度探讨中国石油大学（华东）中国石油大学（华东）第4页，共52页。一、一、 非线性渗流微观机理探讨非线性渗流微观机理探讨v 很多学者的研究证明，启动压力梯度与渗透率呈反很多学者的研究证明，启动压力梯度与渗透率呈反比，渗透率越低，启

6、动压力梯度越大。低渗透油藏孔径比，渗透率越低，启动压力梯度越大。低渗透油藏孔径很小，原油边界层的影响显著，导致宏观表现出来的启很小，原油边界层的影响显著，导致宏观表现出来的启动压力梯度就很明显。动压力梯度就很明显。v 而对于中高渗透油藏，由于孔道半径比较大，原油边界而对于中高渗透油藏，由于孔道半径比较大，原油边界层的影响微弱，尤其对于稀油油藏，流体的屈服应力值很小，层的影响微弱，尤其对于稀油油藏，流体的屈服应力值很小，导致用一般的实验手段不易测到启动压力梯度。因此在实际导致用一般的实验手段不易测到启动压力梯度。因此在实际生产分析时就忽略了启动压力梯度的影响。生产分析时就忽略了启动压力梯度的影响

7、。 1.2 1.2 1.2 1.2 低渗透油藏启动压力梯度探讨低渗透油藏启动压力梯度探讨低渗透油藏启动压力梯度探讨低渗透油藏启动压力梯度探讨中国石油大学（华东）中国石油大学（华东）第5页，共52页。一、一、 非线性渗流微观机理探讨非线性渗流微观机理探讨v 值得一提的是，对于部分低渗透油藏，因非均质状值得一提的是，对于部分低渗透油藏，因非均质状况十分严重，渗流孔径跨度较大，这时会出现况十分严重，渗流孔径跨度较大，这时会出现 实验中也没有观测到启动压力梯度。这是因为虽然油藏整体实验中也没有观测到启动压力梯度。这是因为虽然油藏整体评价属于低渗透油藏，但存在着小部分渗流孔径属于中高评价属于低渗透油藏，

8、但存在着小部分渗流孔径属于中高渗透油层范畴，一般的实验手段也没有观测到启动压力梯渗透油层范畴，一般的实验手段也没有观测到启动压力梯度。因此实际生产分析时部分学者没有考虑启动压力梯度度。因此实际生产分析时部分学者没有考虑启动压力梯度也是合理的。但对于特低渗透油藏忽略启动压力梯度就会也是合理的。但对于特低渗透油藏忽略启动压力梯度就会产生较大的误差。产生较大的误差。 1.2 1.2 1.2 1.2 低渗透油藏启动压力梯度探讨低渗透油藏启动压力梯度探讨低渗透油藏启动压力梯度探讨低渗透油藏启动压力梯度探讨中国石油大学（华东）中国石油大学（华东）第6页，共52页。一、一、 非线性渗流微观机理探讨非线性渗流

9、微观机理探讨 由于低渗岩心的孔隙系统基本上是由小孔道组成由于低渗岩心的孔隙系统基本上是由小孔道组成的，在油、水流动时，每个孔道都有自己的启动压力梯的，在油、水流动时，每个孔道都有自己的启动压力梯度，只有驱动压力梯度大于某孔道的启动压力梯度时，度，只有驱动压力梯度大于某孔道的启动压力梯度时，该孔道中的油、水才开始流动。随着驱动压力梯度的不该孔道中的油、水才开始流动。随着驱动压力梯度的不断提高，就会有更多的孔道加入到流动的行列，岩心的断提高，就会有更多的孔道加入到流动的行列，岩心的渗透性能也随之增强，渗透率变大。渗透性能也随之增强，渗透率变大。1.3 1.3 1.3 1.3 低渗透油藏渗流规律呈现

10、非线性特征探讨低渗透油藏渗流规律呈现非线性特征探讨低渗透油藏渗流规律呈现非线性特征探讨低渗透油藏渗流规律呈现非线性特征探讨1.3.1 低渗透多孔介质的渗透率并非常数中国石油大学（华东）中国石油大学（华东）第7页，共52页。一、一、 非线性渗流微观机理探讨非线性渗流微观机理探讨v 借助大家公认的毛细管模型来分析低渗透油藏的渗流孔道特征尺寸。借助大家公认的毛细管模型来分析低渗透油藏的渗流孔道特征尺寸。v 根据根据孔隙孔道半径与渗透率之间的关系孔隙孔道半径与渗透率之间的关系式式：v 可以计算出三类低渗透油田的平均毛管半径如下：可以计算出三类低渗透油田的平均毛管半径如下：v 第一类为一般低渗透油田，油

11、层平均渗透率为第一类为一般低渗透油田，油层平均渗透率为 。平均。平均毛管半径为（毛管半径为（1.112.471.112.47） 。v 第二类为特低渗透油田，油层平均渗透率为第二类为特低渗透油田，油层平均渗透率为 。平均。平均毛管半径为（毛管半径为（0.371.110.371.11） 。v 第三类为超低渗透油田，其油层平均渗透率为第三类为超低渗透油田，其油层平均渗透率为 。平均毛管。平均毛管半径（半径（0.110.370.110.37） 。1.3.2 低渗透油藏存在微尺度流动效应1.3 1.3 1.3 1.3 低渗透油藏渗流规律呈现非线性特征探讨低渗透油藏渗流规律呈现非线性特征探讨低渗透油藏渗流

12、规律呈现非线性特征探讨低渗透油藏渗流规律呈现非线性特征探讨中国石油大学（华东）中国石油大学（华东）第8页，共52页。一、一、 非线性渗流微观机理探讨非线性渗流微观机理探讨v 通过上述计算可以看出，低渗透油藏渗流通道的通过上述计算可以看出，低渗透油藏渗流通道的特征尺寸已进入微米级别，根据中国科学院陶然、李特征尺寸已进入微米级别，根据中国科学院陶然、李战华等专家的研究结论，微米级别已属于微尺度流动。战华等专家的研究结论，微米级别已属于微尺度流动。但显然远大于分子的特征尺寸，因此虽然连续介质场但显然远大于分子的特征尺寸，因此虽然连续介质场的假设依然成立，但由于流动尺度的微小，许多在的假设依然成立，但

13、由于流动尺度的微小，许多在宏宏观流动中被忽略的或一些影响较小的因素，此时成为观流动中被忽略的或一些影响较小的因素，此时成为主要的影响因素主要的影响因素，从而导致流动规律不同于宏观的流，从而导致流动规律不同于宏观的流动规律，渗流曲线出现非线性段。动规律，渗流曲线出现非线性段。v 也就是说低渗透油藏非线性渗流出现的原因从某也就是说低渗透油藏非线性渗流出现的原因从某种程度上说要归因于微尺度流动效应。种程度上说要归因于微尺度流动效应。1.3.2 低渗透油藏存在微尺度流动效应中国石油大学（华东）中国石油大学（华东）第9页，共52页。 1 1、尺度效应、尺度效应：随着尺度随着尺度减小到微米级时，减小到微米

14、级时， 体积力的作体积力的作用下降；用下降； 表面力的作用却上升，表面力的作用却上升， 超过了体积力；超过了体积力； 随着尺度的减随着尺度的减小，粘性力相对作用增强，惯小，粘性力相对作用增强，惯性力作用变小，越靠近壁面这性力作用变小，越靠近壁面这种规律越明显。种规律越明显。 流变学的结论，因低渗透流变学的结论，因低渗透油藏特征尺度很小，使得速度油藏特征尺度很小，使得速度的空间导数较大，无论是原油的空间导数较大，无论是原油还是地层水在渗流时，均会表还是地层水在渗流时，均会表现出非牛顿的特性，原油的非现出非牛顿的特性，原油的非牛顿特性更强一些。牛顿特性更强一些。 微尺度效应微尺度效应2 2、表面效

15、应：、表面效应：随着尺度的随着尺度的减小，渗流通道的表面积与减小，渗流通道的表面积与体积之比大大增加，流体处体积之比大大增加，流体处在表面上的分子相对于内部在表面上的分子相对于内部分子所占的比例增大，这在分子所占的比例增大，这在很大程度上影响了渗流通道很大程度上影响了渗流通道表面处的质量、动量和能量表面处的质量、动量和能量的传输，表面力将会起主要的传输，表面力将会起主要作用。由于面体比的变化使作用。由于面体比的变化使固体表面的界面效应明显。固体表面的界面效应明显。 一、一、 非线性渗流微观机理探讨非线性渗流微观机理探讨中国石油大学（华东）中国石油大学（华东）第10页，共52页。一、一、 非线性

16、渗流微观机理探讨非线性渗流微观机理探讨 笔者经过多年的研究认为，低渗透油藏笔者经过多年的研究认为，低渗透油藏出现非线性渗流的出现非线性渗流的主要原主要原因因就包括两方面：就包括两方面： 1、低渗透多孔介质的渗透率并非常数，而是随驱替压力梯度变低渗透多孔介质的渗透率并非常数，而是随驱替压力梯度变化而变化化而变化 2 2、低渗透油藏渗流存在低渗透油藏渗流存在微尺度流动效应微尺度流动效应。因低渗透油藏孔喉特。因低渗透油藏孔喉特征尺寸细小，处在连续流动区域以及分子效应区的交界处，虽然连续征尺寸细小，处在连续流动区域以及分子效应区的交界处，虽然连续介质场的假设仍能成立，但由于流动尺度的微小，许多在宏观流

17、动中介质场的假设仍能成立，但由于流动尺度的微小，许多在宏观流动中被忽略的因素，此时成为主要的影响因素，吸附边界层的存在更加加被忽略的因素，此时成为主要的影响因素，吸附边界层的存在更加加剧了微尺度流动效应，从而导致流动规律不同于中高渗油藏的流动规剧了微尺度流动效应，从而导致流动规律不同于中高渗油藏的流动规律，渗流曲线上存在着启动压力梯度以及非线性段，因此达西定律不律，渗流曲线上存在着启动压力梯度以及非线性段，因此达西定律不再适用，需进行修正。再适用，需进行修正。中国石油大学（华东）中国石油大学（华东）第11页，共52页。目录目录第一第一第一第一部分部分部分部分低渗透油藏非线性渗流微观机理探讨低渗

18、透油藏非线性渗流微观机理探讨第二第二第二第二部分部分部分部分低渗透油藏非线性渗流描述方法研究低渗透油藏非线性渗流描述方法研究第三第三第三第三部分部分部分部分低渗透油藏非线性渗流数值模拟方法研究低渗透油藏非线性渗流数值模拟方法研究中国石油大学（华东）中国石油大学（华东）第12页，共52页。二、非线性渗流描述方法研究二、非线性渗流描述方法研究 前文对微观机理的研究指出，达西定律已不能满足对低前文对微观机理的研究指出，达西定律已不能满足对低渗透油藏渗流规律的精确描述，渗透油藏渗流规律的精确描述，需需要进行修正。但因启动要进行修正。但因启动压力梯度随压力梯度变化的特性、渗透率并非常数的特压力梯度随压力

19、梯度变化的特性、渗透率并非常数的特性、存在微尺度流动效应以及边界层，使得要得到能够性、存在微尺度流动效应以及边界层，使得要得到能够精确描述真实渗流规律的运动方程变得十分困难。这就精确描述真实渗流规律的运动方程变得十分困难。这就是目前绝大多数学者是目前绝大多数学者倾向于用唯象的方法来得到运动方倾向于用唯象的方法来得到运动方程的原因程的原因。 目前很多学者对低渗透油藏的非线性渗流描述方法目前很多学者对低渗透油藏的非线性渗流描述方法进行了研究，也提出了很多模型，归纳总结起来有如下进行了研究，也提出了很多模型，归纳总结起来有如下三类三类：中国石油大学（华东）中国石油大学（华东）第13页，共52页。分段

20、模型分段模型 分段模型将分段模型将渗流曲线分为非渗流曲线分为非线性段和线性段，线性段和线性段，并用不同的函数并用不同的函数进行拟合，但在进行拟合，但在应用过程中需要应用过程中需要对临界点进行判对临界点进行判断，应用难度较断，应用难度较大。大。连续模型连续模型 连续模型用统一连续模型用统一的函数对渗流曲线的函数对渗流曲线进行描述，克服了进行描述，克服了常规连续模型不能常规连续模型不能描述渗流存在最小描述渗流存在最小启动压力梯度的问启动压力梯度的问题。题。2.1 2.1 2.1 2.1 非线性渗流描述方法非线性渗流描述方法非线性渗流描述方法非线性渗流描述方法拟启动压力拟启动压力梯度模型梯度模型 拟

21、启动压力梯度拟启动压力梯度模型应用最广泛，模型应用最广泛，但该方法用直线代但该方法用直线代替渗流曲线中的非替渗流曲线中的非线性段，降低了渗线性段，降低了渗流能力，并且抬高流能力，并且抬高了流体是否流动的了流体是否流动的界限，计算结果存界限，计算结果存在着误差。在着误差。非线性渗流描述方法共分为三类：非线性渗流描述方法共分为三类：二、非线性渗流描述方法研究二、非线性渗流描述方法研究中国石油大学（华东）中国石油大学（华东）第14页，共52页。（1 1）由分段函）由分段函数模型向连续模数模型向连续模型函数发展型函数发展非线性渗流非线性渗流模型研究特点模型研究特点（2 2）唯象法模）唯象法模型为主要的

22、数学型为主要的数学模型模型2.1 2.1 2.1 2.1 非线性渗流描述方法非线性渗流描述方法非线性渗流描述方法非线性渗流描述方法二、非线性渗流描述方法研究二、非线性渗流描述方法研究中国石油大学（华东）中国石油大学（华东）第15页，共52页。2.2 2.2 2.2 2.2 非线性渗流模型的选取非线性渗流模型的选取非线性渗流模型的选取非线性渗流模型的选取 式中，式中，vp为相渗流速度，为相渗流速度，m/sm/s；k k为绝对渗为绝对渗透率，；透率，；krp为相相对渗透率；为相相对渗透率； 为相粘度，为相粘度，； 为为p p相的流动势；相的流动势；M M为非线性渗流修正系为非线性渗流修正系数；数；

23、 D D为油藏海拔深度，为油藏海拔深度， m m；G G为拟启动压为拟启动压力梯度，力梯度，MPa/mMPa/m；a a为影响非线性渗流曲线段为影响非线性渗流曲线段的影响因子，的影响因子，a0a0，无量纲；，无量纲；b b为渗流曲线中为渗流曲线中拟启动压力梯度的倒数，拟启动压力梯度的倒数，1/MPa1/MPa；a a、b b与地与地层渗透率有关，由拟合不同渗透率岩芯渗流层渗透率有关，由拟合不同渗透率岩芯渗流实验曲线确定。实验曲线确定。 二、非线性渗流描述方法研究二、非线性渗流描述方法研究中国石油大学（华东）中国石油大学（华东）第16页，共52页。2.2 2.2 2.2 2.2 非线性渗流模型的

24、选取非线性渗流模型的选取非线性渗流模型的选取非线性渗流模型的选取在不同在不同a a，b b取值下，渗流曲线如下图所示：取值下，渗流曲线如下图所示：二、非线性渗流描述方法研究二、非线性渗流描述方法研究 运动方程中，当运动方程中，当a=0a=0时，时，模型即为拟启动压力梯度模型。模型即为拟启动压力梯度模型。当当b b取无穷大取无穷大 时，模型即时，模型即为常规达西渗流模型。因此，为常规达西渗流模型。因此，本模型为通用模型，同样具备本模型为通用模型，同样具备模拟拟启动压力梯度模型和达模拟拟启动压力梯度模型和达西渗流模型的能力。西渗流模型的能力。 中国石油大学（华东）中国石油大学（华东）第17页，共5

25、2页。 特别指出，单相达西定理推广到两相时，引入了相对渗透率特别指出，单相达西定理推广到两相时，引入了相对渗透率的概念，其基本假设就是的概念，其基本假设就是“渠道流态渠道流态”，因此因此引入启动压力梯度引入启动压力梯度以及非线性相关参数时应该取单项时的测量值以及非线性相关参数时应该取单项时的测量值。 部分学者曾研究过两相启动压力梯度，认为两相启动压力梯部分学者曾研究过两相启动压力梯度，认为两相启动压力梯度值要远大于单项启动压力梯度值，且给出的原因是贾敏效应以度值要远大于单项启动压力梯度值，且给出的原因是贾敏效应以及不利毛管力的影响，而实际上相对渗透率已经将贾敏效应以及及不利毛管力的影响，而实际

26、上相对渗透率已经将贾敏效应以及不利毛管力的影响考虑在内了，若启动压力梯度项及非线性相关不利毛管力的影响考虑在内了，若启动压力梯度项及非线性相关参数仍取两相渗流时的测量值，就会出现对贾敏效应以及不利毛参数仍取两相渗流时的测量值，就会出现对贾敏效应以及不利毛管力的影响的管力的影响的重复计算重复计算，这是，这是不合理不合理的。的。二、非线性渗流描述方法研究二、非线性渗流描述方法研究中国石油大学（华东）中国石油大学（华东）第18页，共52页。目录目录第一第一第一第一部分部分部分部分低渗透油藏非线性渗流微观机理探讨低渗透油藏非线性渗流微观机理探讨第二第二第二第二部分部分部分部分低渗透油藏非线性渗流描述方

27、法研究低渗透油藏非线性渗流描述方法研究第三第三第三第三部分部分部分部分低渗透油藏非线性渗流数值模拟方法研究低渗透油藏非线性渗流数值模拟方法研究中国石油大学（华东）中国石油大学（华东）第19页，共52页。三、非线性渗流数模方法研究三、非线性渗流数模方法研究 目前比较成熟的商业软件都是基于目前比较成熟的商业软件都是基于DarcyDarcy渗流规律的，渗流规律的，没有没有考虑第一部分提到的低渗透储层特有的特征。考虑第一部分提到的低渗透储层特有的特征。因此努力开发一因此努力开发一套基于非线性渗流规律的油藏数值模拟软件，为低渗透油田的开发提套基于非线性渗流规律的油藏数值模拟软件，为低渗透油田的开发提供研

28、究手段和工具并用于指导低渗透油田的高效开发，提高低渗透油供研究手段和工具并用于指导低渗透油田的高效开发，提高低渗透油藏开发的主动性和科学性，在目前国际石油形势日趋严峻的情况下具藏开发的主动性和科学性，在目前国际石油形势日趋严峻的情况下具有十分重要的意义。有十分重要的意义。 与中高渗透油藏相比低渗透油藏数值模拟需要考虑的几个与中高渗透油藏相比低渗透油藏数值模拟需要考虑的几个关键问题包括：渗流曲线上存在启动压力梯度和非线性段、压关键问题包括：渗流曲线上存在启动压力梯度和非线性段、压敏效应。敏效应。中国石油大学（华东）中国石油大学（华东）第20页，共52页。3.1.1 3.1.1 启动压力梯度及非线

29、性段处理启动压力梯度及非线性段处理三、非线性渗流数模方法研究三、非线性渗流数模方法研究 本文仍采用前文提到的二参数连续模型的方法：本文仍采用前文提到的二参数连续模型的方法：描述启动压力梯度项及非描述启动压力梯度项及非线性段线性段 流动分三段：流动分三段：A A 不流动段不流动段B B 非线性段非线性段C C 直线段直线段 其中，非线性段上，随其中，非线性段上，随压力的增大启动压力梯度增压力的增大启动压力梯度增大。大。 中国石油大学（华东）中国石油大学（华东）第21页，共52页。 采用这种处理方法，可以根据不同油藏状况对各个网采用这种处理方法，可以根据不同油藏状况对各个网格赋不同的格赋不同的a

30、a、b b值，以模拟启动压力梯度的非均质情况，克值，以模拟启动压力梯度的非均质情况，克服了拟启动压力梯度梯度模型启动压力梯度是定值的弊端，服了拟启动压力梯度梯度模型启动压力梯度是定值的弊端，而且而且X X、Y Y、Z Z三个方向因渗透率存在着差异，在三个方向上三个方向因渗透率存在着差异，在三个方向上a a、b b值可选取不同的值。这样不仅实现了值可选取不同的值。这样不仅实现了启动压力梯度的非均质性启动压力梯度的非均质性，而且实现了而且实现了启动压力梯度方向性差异启动压力梯度方向性差异，对于低渗透油藏数值，对于低渗透油藏数值模拟技术是一大突破。模拟技术是一大突破。三、非线性渗流数模方法研究三、非

31、线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第22页，共52页。(1) (1) 常系数法常系数法3.1.2 3.1.2 压敏效应压敏效应三、非线性渗流数模方法研究三、非线性渗流数模方法研究(2) (2) 变系数法变系数法压实曲线弹性曲线压 力乘子中国石油大学（华东）中国石油大学（华东）第23页，共52页。油气水三油气水三维三相连维三相连续性方程续性方程的建立的建立 油气水及油气水及岩石状态岩石状态方程的建方程的建立立 研究最符研究最符合低渗储合低渗储层实际的层实际的非线性渗非线性渗流描述方流描述方法，建立法，建立非达西渗非达西渗流的运动流的运动方程。方程。 对数学模型对数学模型

32、进行离散化，进行离散化，建立低渗储建立低渗储层差分格式，层差分格式，形成数值模形成数值模型型 以非线性以非线性渗流为特渗流为特点的线性点的线性方程组解方程组解法研究法研究 建立概念模型，开展低渗与中高渗油藏对比分析研究建立概念模型，开展低渗与中高渗油藏对比分析研究 编制低渗透油藏非线性渗流数值模拟软件编制低渗透油藏非线性渗流数值模拟软件 非线性渗流数模方法研究主要内容非线性渗流数模方法研究主要内容三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第24页，共52页。 （1 1）油藏中的流体流动是等温渗流；）油藏中的流体流动是等温渗流； （2 2）油藏中

33、最多只有油、气、水三相，）油藏中最多只有油、气、水三相，油相和水相的渗流均遵循非油相和水相的渗流均遵循非线性渗流定律线性渗流定律； （3 3）油藏烃类只含油气两个组分。在油藏状况下，油、气两种组）油藏烃类只含油气两个组分。在油藏状况下，油、气两种组分可能形成油、气两相。油组分完全存在于油相中，气组分既能以自分可能形成油、气两相。油组分完全存在于油相中，气组分既能以自由气的方式存在，又可溶解于油相之中。地层中的油相应为油组分和由气的方式存在，又可溶解于油相之中。地层中的油相应为油组分和气组分的某种组合；气组分的某种组合； （4 4）油藏中气体的溶解和逸出是瞬间完成的；）油藏中气体的溶解和逸出是瞬

34、间完成的； （5 5）油、水之间不互溶；）油、水之间不互溶； （6 6）岩石微可压缩，且考虑渗流过程中重力、毛管力的影响；）岩石微可压缩，且考虑渗流过程中重力、毛管力的影响； （7 7）地层渗透率是压力梯度的函数地层渗透率是压力梯度的函数。 3.1 3.1 3.1 3.1 基本假设基本假设基本假设基本假设三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第25页，共52页。 油组分方程油组分方程 水组分方程水组分方程 气组分方程气组分方程 3.2 3.2 3.2 3.2 数学模型数学模型数学模型数学模型三、非线性渗流数模方法研究三、非线性渗流数模方法研

35、究中国石油大学（华东）中国石油大学（华东）第26页，共52页。饱和度方程饱和度方程：毛管压力方程毛管压力方程辅助辅助方程方程边界边界边界边界条件条件条件条件外边界条件：外边界条件： 外边界可以有外边界可以有3 3种形式：定压外边界、种形式：定压外边界、定流量外边界和混合定流量外边界和混合外边界。外边界。初始初始条件条件 给定初始时刻每给定初始时刻每个点的压力和各相饱个点的压力和各相饱和度值。和度值。内边界条件：内边界条件： 定产量（或定定产量（或定注水量）、定井底注水量）、定井底流压。流压。3.3 3.3 定解定解条件条件三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）

36、中国石油大学（华东）第27页，共52页。3.4 3.4 3.4 3.4 非线性模型差分离散化方法非线性模型差分离散化方法非线性模型差分离散化方法非线性模型差分离散化方法 采用块中心七点有限差分方法。与常规黑油模型不同的是油、水采用块中心七点有限差分方法。与常规黑油模型不同的是油、水相方程中含一项相方程中含一项非线性渗流修正系数非线性渗流修正系数，修正系数的大小与各相的压，修正系数的大小与各相的压力梯度有关，取值范围为力梯度有关，取值范围为00，11。因此，将该修正系数作为。因此，将该修正系数作为绝对渗透绝对渗透率的修正因子率的修正因子，不仅理论上可以解释清楚，而且能够保证求解方程中系数，不仅理

37、论上可以解释清楚，而且能够保证求解方程中系数矩阵的稳定。矩阵的稳定。 但是，考虑非线性渗流的油藏数值模拟器必须对以下两方面进行特但是，考虑非线性渗流的油藏数值模拟器必须对以下两方面进行特殊处理：殊处理： 3.4.1 3.4.1 非线性渗流修正系数非线性渗流修正系数 3.4.2 3.4.2 井井- -网格流动方程网格流动方程 三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第28页，共52页。 3.4.1 3.4.1 非线性渗流修正系数非线性渗流修正系数 从实用性考虑，本模型必须适应油藏非均质的情况，因此，必须求从实用性考虑，本模型必须适应油藏非均质的

38、情况，因此，必须求得在交界面处非线性渗流修正系数的值。假设离散网格系统由六面体组得在交界面处非线性渗流修正系数的值。假设离散网格系统由六面体组成，网格的成，网格的x x、y y、z z方向编号分别为方向编号分别为i i、j j、k k。以。以x x方向为例，根据质量方向为例，根据质量守恒原理：守恒原理： 交界面处的非线性渗流修正系数：交界面处的非线性渗流修正系数： 同理可得有类似的表达式。同理可得有类似的表达式。三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第29页，共52页。 3.4.2 3.4.2 井井- -网格流动方程网格流动方程 为了保证所

39、提数学模型的适定性，井为了保证所提数学模型的适定性，井- -网格压力方程中也必须考网格压力方程中也必须考虑非线性渗流的影响。井所在的网格在数值模拟中处理为等效虑非线性渗流的影响。井所在的网格在数值模拟中处理为等效的平面径向流，考虑一个时间步内井与其所在网格间的稳定渗的平面径向流，考虑一个时间步内井与其所在网格间的稳定渗流，设等效半径为：流，设等效半径为： 圆形油藏内为均匀介质，设该时间步内网格圆形油藏内为均匀介质，设该时间步内网格- -井间井间p p相的稳定流率相的稳定流率为为QpQp，则任意半径处的流速为：，则任意半径处的流速为： 式中式中当当时，对应于汇项。时，对应于汇项。时，对应于源项。

40、时，对应于源项。三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第30页，共52页。 对上式进行积分并求取等效供液半径内的平均势，考虑到近对上式进行积分并求取等效供液半径内的平均势，考虑到近井地带压力梯度大于临界启动压力梯度，可取井地带压力梯度大于临界启动压力梯度，可取a=0a=0，并且，化简，并且，化简可得：可得： 式中式中 表示井底流动势表示井底流动势。应该与网格应该与网格p p相流体的流动势。相流体的流动势。由上式可得：由上式可得：三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第31页，共52页。相

41、等，可得相等，可得： 上式即为模拟非线性渗流的井上式即为模拟非线性渗流的井- -网格方程，其中正负号对生产井取负值，注网格方程，其中正负号对生产井取负值，注入井取正值。当用本模型模拟达西渗流时，井入井取正值。当用本模型模拟达西渗流时，井- -网格方程与达西渗流一致。可以网格方程与达西渗流一致。可以看出，井与网格间存在启动压力差，只有压差达到该值后才会发生流动，看出，井与网格间存在启动压力差，只有压差达到该值后才会发生流动，因此，因此，低渗层的实际产出和注入量远低于按照达西流规律所得到的数值。低渗层的实际产出和注入量远低于按照达西流规律所得到的数值。应该与网格应该与网格p p相流体的流动势相流体

42、的流动势三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第32页，共52页。3.5 3.5 3.5 3.5 非线性数值模型解法研究非线性数值模型解法研究非线性数值模型解法研究非线性数值模型解法研究 黑油模型中，在黑油模型中，在x x方向上，网格块方向上，网格块(i,j,k)(i,j,k)与相邻网格块间的相传导与相邻网格块间的相传导率可定义如下：率可定义如下： 三、非线性渗流数模方法研究三、非线性渗流数模方法研究 在黑油模型中，在黑油模型中， 为定值，为定值， 为弱非线性参数，为弱非线性参数， 为强非线性参数。本文建立的渗流模型与黑油模型存在很大的不同

43、，本文引为强非线性参数。本文建立的渗流模型与黑油模型存在很大的不同，本文引入了非线性渗流修正系数修正绝对渗透率，该修正因子是与压力梯度有关的；本文还入了非线性渗流修正系数修正绝对渗透率，该修正因子是与压力梯度有关的；本文还考虑了压敏效应，同样的是修正绝对渗透率，因此传导率的计算公式中的考虑了压敏效应，同样的是修正绝对渗透率，因此传导率的计算公式中的 项项由黑油模型计算时的定值变成了强非线性的，普通的解法不能保证计算的稳定性和正由黑油模型计算时的定值变成了强非线性的，普通的解法不能保证计算的稳定性和正确性，因此需采用确性，因此需采用全隐式解法全隐式解法。 中国石油大学（华东）中国石油大学（华东）

44、第33页，共52页。3.5 3.5 3.5 3.5 非线性数值模型解法研究非线性数值模型解法研究非线性数值模型解法研究非线性数值模型解法研究 下面就上述提到的对绝对渗透率进行非线性渗流系数修正以及压敏效下面就上述提到的对绝对渗透率进行非线性渗流系数修正以及压敏效应修正的处理方式进行介绍如下：应修正的处理方式进行介绍如下： 利用有限差分法将非线性渗流微分方程离散化，其中绝对渗透率为利用有限差分法将非线性渗流微分方程离散化，其中绝对渗透率为压力梯度和压差的函数。在迭代计算中，先对方程中的压力和饱和度赋压力梯度和压差的函数。在迭代计算中，先对方程中的压力和饱和度赋初值，然后根据赋的初值使系数线性化，

45、迭代求解方程可得到压力和饱初值，然后根据赋的初值使系数线性化，迭代求解方程可得到压力和饱和度值，然后得到各个单元的压力梯度值，再由渗透率曲线得到各个单和度值，然后得到各个单元的压力梯度值，再由渗透率曲线得到各个单元的渗透率值。将此渗透率和刚计算出的饱和度再代入非线性渗流方程元的渗透率值。将此渗透率和刚计算出的饱和度再代入非线性渗流方程中又使其线性化。循环交替迭代，直至压力和饱和度值趋于稳定并满足中又使其线性化。循环交替迭代，直至压力和饱和度值趋于稳定并满足计算精度要求，此值即为对应时段的压力和饱和度值，然后进入下一个计算精度要求，此值即为对应时段的压力和饱和度值，然后进入下一个时段的计算。时段

46、的计算。三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第34页，共52页。 3.6 3.6 3.6 3.6 算例计算及对比分析算例计算及对比分析算例计算及对比分析算例计算及对比分析 为了检验该数值模拟软件的正确性，采用该软件对不同井距的反九为了检验该数值模拟软件的正确性，采用该软件对不同井距的反九点面积井网概念模型进行了数值模拟计算。结合某特低渗透油藏实际点面积井网概念模型进行了数值模拟计算。结合某特低渗透油藏实际储层物性和流体性质，基本油藏参数选取如下：孔隙度储层物性和流体性质，基本油藏参数选取如下：孔隙度0.150.15；渗透率；渗透率5105

47、10-3-3mm2 2；a=0.8a=0.8，b=14.0(1/Mpa)b=14.0(1/Mpa)；地层原油相对密度；地层原油相对密度0.80.8；粘度；粘度3 3；原始地层压力原始地层压力17.517.5；模型的顶面深度；模型的顶面深度2000m2000m；有效厚度；有效厚度15m15m；方向网格数；方向网格数为：为：Nx=41Nx=41，Ny=41Ny=41，Nz=1Nz=1。工作制度以井底流压为约束条件，注水井取。工作制度以井底流压为约束条件，注水井取25MPa25MPa；采油井取；采油井取10MPa10MPa；模拟时间为；模拟时间为2020年，最大时间步长为年，最大时间步长为1010天

48、。应用该概天。应用该概念模型模拟预测了念模型模拟预测了达西渗流模型、非线性渗流模型和拟启动压力梯度模达西渗流模型、非线性渗流模型和拟启动压力梯度模型型在不同井距下的油藏动态情况。在不同井距下的油藏动态情况。 三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第35页，共52页。三、非线性渗流数模方法研究三、非线性渗流数模方法研究a.井距150m b.井距200mc.井距250m不同井距下压力梯度等值线图不同井距下压力梯度等值线图 渗流曲线上压力梯度弯曲段为渗流曲线上压力梯度弯曲段为中国石油大学（华东）中国石油大学（华东）第36页，共52页。 渗流曲线上

49、弯曲段为，可以看出，除井点处压力梯渗流曲线上弯曲段为，可以看出，除井点处压力梯度较大外，油藏中大部分区域压力梯度比较低，处于非线度较大外，油藏中大部分区域压力梯度比较低，处于非线性渗流曲线的弯曲段。性渗流曲线的弯曲段。非线性渗流模型更准确。非线性渗流模型更准确。 同时可以看出，同时可以看出，随着井距的缩小，油藏动用范围增加随着井距的缩小，油藏动用范围增加，井距为井距为150m150m时，压力梯度较大，边角井间都可以建立起连通关系；时，压力梯度较大，边角井间都可以建立起连通关系；井距为井距为250m250m时，注水井和角井间存在不流动区，只有边井见到注时，注水井和角井间存在不流动区，只有边井见到

50、注水效果。水效果。三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第37页，共52页。3.6 3.6 3.6 3.6 算例计算及对比分析算例计算及对比分析算例计算及对比分析算例计算及对比分析 左图可以看出，左图可以看出，达西渗流模型的开发效达西渗流模型的开发效果最好果最好，在相同情况下，非线性渗流模型，在相同情况下，非线性渗流模型和拟启动压力梯度模型的开发效果比达西和拟启动压力梯度模型的开发效果比达西渗流模型差的多。渗流模型差的多。200m200m达西渗流模型的采达西渗流模型的采出程度与出程度与150m150m非线性渗流模型的相近。非线性渗流模型的相

51、近。相同井距下，非线性渗流模型的采出程相同井距下，非线性渗流模型的采出程度比拟启动压力梯度模型的大度比拟启动压力梯度模型的大，并且这，并且这种趋势随着井距的增大而增加。这主要是因种趋势随着井距的增大而增加。这主要是因为：油藏压力梯度大部分处于渗流曲线的非为：油藏压力梯度大部分处于渗流曲线的非线性弯曲段，拟启动压力梯度模型降低了渗线性弯曲段，拟启动压力梯度模型降低了渗流速度的大小。流速度的大小。 不同井距不同流动规律下采出程度对比图不同井距不同流动规律下采出程度对比图三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第38页，共52页。3.6 3.6 3

52、.6 3.6 算例计算及对比分析算例计算及对比分析算例计算及对比分析算例计算及对比分析 左图给出了三种模型的综合含水对左图给出了三种模型的综合含水对比。可以看出，达西渗流模型见水时间比。可以看出，达西渗流模型见水时间最早，非线性渗流模型和拟启动压力梯最早，非线性渗流模型和拟启动压力梯度模型见水时间较达西流模型滞后，后度模型见水时间较达西流模型滞后，后两者比较，两者比较，非线性渗流模型见水时间非线性渗流模型见水时间较早较早，反映出将特低渗透油藏渗流曲线处，反映出将特低渗透油藏渗流曲线处理为拟启动压力梯度的模型时造成的误差。理为拟启动压力梯度的模型时造成的误差。当井距增大时，非线性模型与拟启动压力

53、当井距增大时，非线性模型与拟启动压力梯度模型的差距变大。当井距较大时，特梯度模型的差距变大。当井距较大时，特低渗透油藏注水见效慢，这符合已开发油低渗透油藏注水见效慢，这符合已开发油藏反映出的动态特征。藏反映出的动态特征。 不同流动规律下含水率对比图不同流动规律下含水率对比图 三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第39页，共52页。 但随着对低渗透砂岩油藏的开发和注水的深入，人们发现裂缝的作用越来但随着对低渗透砂岩油藏的开发和注水的深入，人们发现裂缝的作用越来越重要。裂缝不仅决定了注水效果，而且控制了层系划分和井网布置，从而越重要。裂缝不仅

54、决定了注水效果，而且控制了层系划分和井网布置，从而直接决定了油田开发效果的好坏。近些年来，我国发现的直接决定了油田开发效果的好坏。近些年来，我国发现的裂缝性低渗透砂岩裂缝性低渗透砂岩油田越来越多油田越来越多。因此研究低渗透裂缝性砂岩油藏的开发规律就显得十分重。因此研究低渗透裂缝性砂岩油藏的开发规律就显得十分重要。要。 本文针对上述问题，首先对双重变形介质的结构特征进行了分析，本文针对上述问题，首先对双重变形介质的结构特征进行了分析，以现有的双孔双渗模型、双孔单渗模型以及前文建立的单重变形介质以现有的双孔双渗模型、双孔单渗模型以及前文建立的单重变形介质非线性渗流模型为基础，建立了一个能够综合体现

55、裂缝、压敏效应、非线性渗流模型为基础，建立了一个能够综合体现裂缝、压敏效应、油相和水相的渗流遵循非线性渗流定律以及存在压力梯度特征的油相和水相的渗流遵循非线性渗流定律以及存在压力梯度特征的复合复合双重变形介质非线性渗流数学模型双重变形介质非线性渗流数学模型。双重变形介质非线性渗流数值模拟方法研究双重变形介质非线性渗流数值模拟方法研究 三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第40页，共52页。v复合双重变形介质非线性渗流数学模型复合双重变形介质非线性渗流数学模型：v(1)(1)裂缝系统裂缝系统 双重变形介质非线性渗流数值模拟方法研究双重变形介

56、质非线性渗流数值模拟方法研究 三、非线性渗流数模方法研究三、非线性渗流数模方法研究油组分方程：油组分方程：水组分方程：水组分方程：气组分方程：气组分方程：中国石油大学（华东）中国石油大学（华东）第41页，共52页。v(2) (2) 基质系统基质系统v油组分方程：油组分方程：双重变形介质非线性渗流数值模拟方法研究双重变形介质非线性渗流数值模拟方法研究 三、非线性渗流数模方法研究三、非线性渗流数模方法研究 水组分方程：水组分方程： 气组分方程：气组分方程： 中国石油大学（华东）中国石油大学（华东）第42页，共52页。双重变形介质非线性渗流数值模拟方法研究双重变形介质非线性渗流数值模拟方法研究 三、

57、非线性渗流数模方法研究三、非线性渗流数模方法研究 值得注意的是，我们建立的复合双重变形介质非线性渗值得注意的是，我们建立的复合双重变形介质非线性渗流数学模型，可以处理裂缝性低渗透油藏所出现的复杂的渗流数学模型，可以处理裂缝性低渗透油藏所出现的复杂的渗流情形：流情形： 1 1、基质的压敏效应，裂缝随压力变化的张启和闭合；、基质的压敏效应，裂缝随压力变化的张启和闭合； 2 2、启动压力梯度和非线性段对基质流向裂缝渗流的影、启动压力梯度和非线性段对基质流向裂缝渗流的影响；响； 3 3、单重介质和双重介质的耦合；、单重介质和双重介质的耦合；中国石油大学（华东）中国石油大学（华东）第43页，共52页。三

58、、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第44页，共52页。三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第45页，共52页。三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第46页，共52页。三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第47页，共52页。 编制的软件能够反映出低渗透油藏特有的特征，与编制的软件能够反映出低渗透油藏特有的特征，与基于拟启动压力梯度模型以及达西定律的软件相比，该基于拟启动压力梯

59、度模型以及达西定律的软件相比，该软件更能反映低渗透油藏的特殊的渗流机理和开发动态软件更能反映低渗透油藏的特殊的渗流机理和开发动态 ，为井网优化设计以及有效驱替方式研究提供了技术支为井网优化设计以及有效驱替方式研究提供了技术支持，持，对于指导低渗透油藏的开发具有十分深远的意义。对于指导低渗透油藏的开发具有十分深远的意义。 而且该软件同时配备有强大的前后处理模块，很好的而且该软件同时配备有强大的前后处理模块，很好的实现了一体化，为该软件的工程应用奠定了基础。实现了一体化，为该软件的工程应用奠定了基础。三、非线性渗流数模方法研究三、非线性渗流数模方法研究中国石油大学（华东）中国石油大学（华东）第48

60、页，共52页。四、结论四、结论第一部分第一部分 低渗透油藏非线性渗流微观机理探讨低渗透油藏非线性渗流微观机理探讨 1 1、启动压力梯度存在的本质是因为流体存在着屈服值，、启动压力梯度存在的本质是因为流体存在着屈服值，再加上边界层的影响使得启动压力梯度出现随驱替压力梯度再加上边界层的影响使得启动压力梯度出现随驱替压力梯度变化而变化的特点。变化而变化的特点。 ( (1)1)、低渗透多孔介质的渗透率并非常数，而是随驱替压、低渗透多孔介质的渗透率并非常数，而是随驱替压力梯度变化而变化力梯度变化而变化 (2)(2)、低渗透油藏渗流存在、低渗透油藏渗流存在微尺度流动效应。因低渗透微尺度流动效应。因低渗透油

61、藏孔喉特征尺寸细小，处在连续流动区域以及分子效应区的油藏孔喉特征尺寸细小，处在连续流动区域以及分子效应区的交界处，虽然连续介质场的假设仍能成立，但由于流动尺度的交界处，虽然连续介质场的假设仍能成立，但由于流动尺度的微小，许多在宏观流动中被忽略的因素，此时成为主要的影响微小，许多在宏观流动中被忽略的因素，此时成为主要的影响因素，因素，吸附边界层的存在更加加剧了微尺度流动效应吸附边界层的存在更加加剧了微尺度流动效应，从，从而导致流动规律不同于中高渗油藏的流动规律，渗流而导致流动规律不同于中高渗油藏的流动规律，渗流曲线上存在着启动压力梯度以及非线性段，因此达西曲线上存在着启动压力梯度以及非线性段，因

62、此达西定律不再适用，需进行修正。定律不再适用，需进行修正。启动压启动压力梯度力梯度非非线线性性渗渗流流中国石油大学（华东）中国石油大学（华东）第49页，共52页。四、结论四、结论第二部分第二部分 低渗透油藏非线性渗流描述方法研究低渗透油藏非线性渗流描述方法研究1 1、非线性渗流模型特点：、非线性渗流模型特点： 由分段函数向连续函数发展、唯象法模型为主要的数学模型；由分段函数向连续函数发展、唯象法模型为主要的数学模型；2 2、选取二参数模型做为后续分析处理的描述模型。原因有：、选取二参数模型做为后续分析处理的描述模型。原因有： (1) (1) 可以实现模拟启动压力梯度的可变性、非均质性和方向差异

63、可以实现模拟启动压力梯度的可变性、非均质性和方向差异性；性； (2) (2) 很好的反应了非线性特征；很好的反应了非线性特征； (3) (3) 可以实现同时模拟达西流、拟启动压力梯度模型；可以实现同时模拟达西流、拟启动压力梯度模型；3 3、深入分析了相对渗透率的内涵以及渠道流理论，得出单相渗流模、深入分析了相对渗透率的内涵以及渠道流理论，得出单相渗流模型推广到两相、三相时，模型中的启动压力梯度项及非线性相关参数型推广到两相、三相时，模型中的启动压力梯度项及非线性相关参数应该取单相时的测量值。应该取单相时的测量值。中国石油大学（华东）中国石油大学（华东）第50页，共52页。四、结论四、结论第三部

64、分第三部分 低渗透油藏非线性渗流数值模拟方法研究低渗透油藏非线性渗流数值模拟方法研究 1 1、首次实现能同时考虑启动压力梯度、非线性段以及压敏效应、首次实现能同时考虑启动压力梯度、非线性段以及压敏效应对渗流影响的数值模拟器；对渗流影响的数值模拟器； 2 2、编制了三维三相非线性渗流数值模拟软件，该软件同样具备模拟常规达、编制了三维三相非线性渗流数值模拟软件，该软件同样具备模拟常规达西渗流模型和拟启动压力梯度模型的能力，为低渗透数值模拟软件的工程化应西渗流模型和拟启动压力梯度模型的能力，为低渗透数值模拟软件的工程化应用奠定了基础；用奠定了基础； 3 3、算例计算表明：低渗透油藏开发过程中压力梯度大部分位于非线、算例计算表明：低渗透油藏开发过程中压力梯度大部分位于非线性弯曲段，非线性渗流模型与达西模型和拟启动压力梯度模型预测结果性弯曲段，非线性渗流模型与达西模型和拟启动压力梯度模型预测结果存在着差别。随着井距的减小，油藏动用范围增加。存在着差别。随着井距的减小，油藏动用范围增加。 4 4、建立的复合双重变形介质非线性渗流数学模型，可以实现对裂缝性低渗、建立的复合双重变形介质非线性渗流数学模型，可以实现对裂缝性低渗透油藏的精确模拟。透油藏的精确模拟。中国石油大学（华东）中国石油大学（华东）第51页，共52页。第52页，共52页。