《第七章多孔介质的渗流优秀课件.ppt》由会员分享，可在线阅读，更多相关《第七章多孔介质的渗流优秀课件.ppt（34页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第七章多孔介质的渗流第1页，本讲稿共34页5.1 5.1 渗流过程中的力学分析及驱动类型渗流过程中的力学分析及驱动类型 油、水、气能够在多孔介质（岩石）中渗流是由于各种油、水、气能够在多孔介质（岩石）中渗流是由于各种力的作用，主要有以下几种：力的作用，主要有以下几种：流体的重力流体的重力 1多孔介质（岩石）的压缩性及流体的弹性力多孔介质（岩石）的压缩性及流体的弹性力 2重力有时是动力有时是阻力。重力有时是动力有时是阻力。在油气开采以前，岩石和流体都处于均衡受压状态，在油气开采以前，岩石和流体都处于均衡受压状态，当油气当油气层投入开采之后，油气层的压力不断下降，上覆岩层和油层内的层投入开采之后，
2、油气层的压力不断下降，上覆岩层和油层内的流体压力之间形成压力差，岩石变形，岩石孔隙体积减少，压缩流体压力之间形成压力差，岩石变形，岩石孔隙体积减少，压缩孔隙中的流体驱使流体向压力较低的方向运动。孔隙中的流体驱使流体向压力较低的方向运动。第2页，本讲稿共34页流体的粘度及粘滞力流体的粘度及粘滞力 4毛管力毛管力 3 多孔介质是由无数个毛细管组成，这些毛细管纵横交错，四通八达，当渗多孔介质是由无数个毛细管组成，这些毛细管纵横交错，四通八达，当渗流由一种流体驱替另一种流体时，在两界面上产生压力跳跃，它的大小取决流由一种流体驱替另一种流体时，在两界面上产生压力跳跃，它的大小取决于分界面的曲度，这个压力
3、就称为于分界面的曲度，这个压力就称为毛管压力。毛管压力。在流动的流体中，如果各种流体流速不同，将有一对作用力和反作用力，在流动的流体中，如果各种流体流速不同，将有一对作用力和反作用力，使原来快的流层减速，而慢的加速。流体的这种属性叫使原来快的流层减速，而慢的加速。流体的这种属性叫粘滞性粘滞性。在渗流中，粘滞力为阻力，且动力消耗主要用于渗流时克服流体粘在渗流中，粘滞力为阻力，且动力消耗主要用于渗流时克服流体粘滞阻力。滞阻力。第3页，本讲稿共34页 渗流的驱动类型渗流的驱动类型主要有：重力水压驱动、弹性驱动、主要有：重力水压驱动、弹性驱动、气压驱动、溶解驱动以及重力驱动。气压驱动、溶解驱动以及重力
4、驱动。在渗流过程中必有一种或多种驱动方式起重要在渗流过程中必有一种或多种驱动方式起重要作用，其他驱动类型处于从属地位。驱动方式在渗作用，其他驱动类型处于从属地位。驱动方式在渗流过程中不是一成不变的，而是变化发展的！流过程中不是一成不变的，而是变化发展的！第4页，本讲稿共34页5.2 5.2 不可压缩流体渗流及渗透率的张量特性不可压缩流体渗流及渗透率的张量特性为地层渗透为地层渗透率，是一个张量率，是一个张量K 对于不可压缩流体地渗流问题，引入运动方程即达西公式为：对于不可压缩流体地渗流问题，引入运动方程即达西公式为：K可以表示为：可以表示为：在特殊情况下在特殊情况下,渗透率张量渗透率张量K K是
5、实对称的是实对称的,所以至少存在三所以至少存在三个相互垂直主方向个相互垂直主方向 。标基矢量，则渗透率张量的矩阵为对角型。标基矢量，则渗透率张量的矩阵为对角型。，若选，若选 为坐为坐第5页，本讲稿共34页即：即：其张量为：其张量为：设设 是原有已坐标系，其单位矢量为是原有已坐标系，其单位矢量为 (i,j=1,2,3)是变换了的已坐标系，其单位矢量为：是变换了的已坐标系，其单位矢量为：其张量为：其张量为：(i,j=1,2,3)。变换规律为：变换规律为：坐标变换规律坐标变换规律第6页，本讲稿共34页单相渗流连续性方程的张量形式为：单相渗流连续性方程的张量形式为：为地层为地层孔隙度孔隙度 对于稳定渗
6、流对于稳定渗流若流体是不可压缩的若流体是不可压缩的可以略去可以略去 第7页，本讲稿共34页5.3 5.3 两相渗流问题两相渗流问题建立数学模型建立数学模型 对于油相的连续性方程为：对于油相的连续性方程为：对于水相的连续性方程为：对于水相的连续性方程为：设油水两种流体同时在多孔介质中流动，设油水两种流体同时在多孔介质中流动，且流动服从达西定律且流动服从达西定律。第8页，本讲稿共34页饱和度方程：饱和度方程：在考虑毛管力和重力影响时，油相和水相的渗流速度分别为：在考虑毛管力和重力影响时，油相和水相的渗流速度分别为：是流动方向与是流动方向与水平面的夹角水平面的夹角引入拉普拉斯方程把油相和水相压力联系
7、起来：引入拉普拉斯方程把油相和水相压力联系起来：为毛细管液为毛细管液面的主半径面的主半径 为两相界面为两相界面上的界面张力上的界面张力 第9页，本讲稿共34页独立方程总数为独立方程总数为6 6个共有个共有6 6个待求的因变量个待求的因变量组成一个封闭方程组。组成一个封闭方程组。小结：小结：第10页，本讲稿共34页5.4 5.4 气体渗流问题气体渗流问题 气体比液体具有更大的压缩性。在研究气体渗流规律时，气体比液体具有更大的压缩性。在研究气体渗流规律时，可以仿照液体流体的研究方法，得出相似形式渗流方程，但它可以仿照液体流体的研究方法，得出相似形式渗流方程，但它们在物理实质则有差别。们在物理实质则
8、有差别。气体渗流的基本特征气体渗流的基本特征1状态方程和基本特性参数状态方程和基本特性参数真实气体的状态方程：真实气体的状态方程：Z Z：压缩因子：压缩因子（亦称偏差因（亦称偏差因子）子）T T为绝对温为绝对温度度R R为气体为气体常数常数 第11页，本讲稿共34页气体的压缩系数：气体的压缩系数：气体的体积系数气体的体积系数 ：气体在标准条件下的密度气体在标准条件下的密度：第12页，本讲稿共34页气体渗透率气体渗透率为克氏为克氏渗透率渗透率 b b为孔隙大小和为孔隙大小和分子自由程所决分子自由程所决定的参数定的参数K K为在平均压力为在平均压力 下和平均流量下和平均流量 下测得的气体下测得的气
9、体渗透率渗透率第13页，本讲稿共34页气体导压系数气体导压系数 在气体渗流中，由于气体压缩系数大大高于岩石的压缩系数，因在气体渗流中，由于气体压缩系数大大高于岩石的压缩系数，因此往往在非稳定气体渗流中使用气体压缩系数代替综合压缩系数，这此往往在非稳定气体渗流中使用气体压缩系数代替综合压缩系数，这样气体的导压系数为：样气体的导压系数为：但是对某些低渗气层可能不适用，此时可以采用综合但是对某些低渗气层可能不适用，此时可以采用综合压缩系数压缩系数：第14页，本讲稿共34页2气体的稳定渗流气体的稳定渗流 气体稳定渗流的基本方程如下：气体稳定渗流的基本方程如下：运动方程运动方程 K K为渗透为渗透率张量
10、率张量 状态方程状态方程 若是理想气体，则有：若是理想气体，则有：第15页，本讲稿共34页连续性方程连续性方程 即：即：将运动方程和状态方程带入连续性方程，因为：将运动方程和状态方程带入连续性方程，因为：第16页，本讲稿共34页设设代表地层平均压力，并认为代表地层平均压力，并认为是一个常数是一个常数 称为气体的导压系数称为气体的导压系数 令令上式即为理想气体和稳定渗流的数学模型。上式即为理想气体和稳定渗流的数学模型。第17页，本讲稿共34页理想气体和稳定渗流的数学模型的适用条件是：理想气体和稳定渗流的数学模型的适用条件是：1 1）气体是单相渗流的；）气体是单相渗流的；2 2）符合线性渗流运动方
11、程；）符合线性渗流运动方程；3 3）气体为可压缩的理想气体；）气体为可压缩的理想气体；4 4）岩石的压缩性忽略不计，孔隙度视为常数；）岩石的压缩性忽略不计，孔隙度视为常数；5 5）渗流过程是等温的。）渗流过程是等温的。上式即变为气体稳定渗流的数学模型：上式即变为气体稳定渗流的数学模型：当当 时，时，第18页，本讲稿共34页3气体的不稳定渗流气体的不稳定渗流 在气体渗流中，在气体渗流中，压力梯度与渗流速度往往不是线性关系，压力梯度与渗流速度往往不是线性关系，即即出现非达西渗流。这里仅取由量纲分析导出的二项式达出现非达西渗流。这里仅取由量纲分析导出的二项式达西渗流，表达式为：西渗流，表达式为：考虑
12、气体的流动惯性，在线性达西渗流条件下，便可考虑气体的流动惯性，在线性达西渗流条件下，便可简化为通常的渗流速度公式：简化为通常的渗流速度公式：是与孔隙是与孔隙度和孔隙结构度和孔隙结构有关的常数有关的常数 第19页，本讲稿共34页根据连续性方程，对于气体渗流问题，有：根据连续性方程，对于气体渗流问题，有：当当 且且b0时，有：时，有：由于由于第20页，本讲稿共34页上式为上式为不稳定渗流的基本方程不稳定渗流的基本方程 取取 称为气体得导压系数，它是一个与称为气体得导压系数，它是一个与压力有关的变量。因而，无法采用原来的方法求解。压力有关的变量。因而，无法采用原来的方法求解。在通常情况下，我们可以采
13、用线性化方法（求解非稳定问在通常情况下，我们可以采用线性化方法（求解非稳定问题），稳定状态逐次替换法。在某些情况下可用自模解的方法题），稳定状态逐次替换法。在某些情况下可用自模解的方法获得精确解，但后者往往限于典型问题，并最终要求用数值方获得精确解，但后者往往限于典型问题，并最终要求用数值方法求解常微分方程。法求解常微分方程。第21页，本讲稿共34页5.5 5.5 双重介质渗流双重介质渗流 所谓的双重介质所谓的双重介质就是分别由裂缝和孔隙构成，二者又是全空就是分别由裂缝和孔隙构成，二者又是全空间叠合在一起且相互影响的一个统一的渗流场。间叠合在一起且相互影响的一个统一的渗流场。把这种双重介质孔隙
14、结构地层典型的化为由互把这种双重介质孔隙结构地层典型的化为由互相垂直的裂缝系统和被裂缝系统所切割开的岩块组成，相垂直的裂缝系统和被裂缝系统所切割开的岩块组成，这就是所谓的这就是所谓的双重介质渗流模型双重介质渗流模型。主要有。主要有KozemiKozemi模型、模型、De.SwanDe.Swan模型和模型和Warren-RootWarren-Root模型等，这样的模型既能模型等，这样的模型既能体现地层的双重孔隙性，也能体现其双重渗透性。体现地层的双重孔隙性，也能体现其双重渗透性。定义定义第22页，本讲稿共34页 由于两种介质储集性能和渗透性能的不同，使得压由于两种介质储集性能和渗透性能的不同，使
15、得压力传播的速度不同，因此当渗流时，在空间的任意点应力传播的速度不同，因此当渗流时，在空间的任意点应同时引进两个压力（裂缝中的压力同时引进两个压力（裂缝中的压力 和岩块的压力和岩块的压力 ）和两个渗流速度（裂缝的渗流速度）和两个渗流速度（裂缝的渗流速度 和岩块的渗流速和岩块的渗流速度度 ）这样就行形成了两个平行的渗流场，并且两个渗）这样就行形成了两个平行的渗流场，并且两个渗流场之间存在着流体交换，这种流体交换的物理现象称流场之间存在着流体交换，这种流体交换的物理现象称为为“窜流窜流”。第23页，本讲稿共34页1运动方程及窜流方程运动方程及窜流方程 设裂缝系统和团块中的渗流均服从达西渗流定律，设
16、裂缝系统和团块中的渗流均服从达西渗流定律，运动方程运动方程为：为：裂缝系统裂缝系统 岩块系统岩块系统 裂缝和孔隙之间的流体交换主要是在比较平缓的压力变化下发裂缝和孔隙之间的流体交换主要是在比较平缓的压力变化下发生的，因此在这个过程中，可以认为窜流和时间无关系。生的，因此在这个过程中，可以认为窜流和时间无关系。第24页，本讲稿共34页窜流的流量主要取决于：窜流的流量主要取决于：（1 1）流体的粘度；）流体的粘度；（2 2）孔隙和裂缝之间的压差）孔隙和裂缝之间的压差 （3 3）岩石的某些特征量，如长度、面积、）岩石的某些特征量，如长度、面积、体积单位等等。体积单位等等。经过因次分析得出经过因次分析
17、得出窜流方程窜流方程为：为：称为称为窜流系数，窜流系数，无量纲量无量纲量q q为单位时间单位岩石流为单位时间单位岩石流出的流体质量出的流体质量 第25页，本讲稿共34页2状态方程状态方程 无论是孔隙介质还是裂缝介质以及地层中流体均是为可压缩无论是孔隙介质还是裂缝介质以及地层中流体均是为可压缩的的，则裂缝和岩块的孔隙度的压缩特性公式是：的的，则裂缝和岩块的孔隙度的压缩特性公式是：裂缝系统裂缝系统 岩块系统岩块系统 C C为压缩为压缩系数系数 为某个已为某个已知压力知压力第26页，本讲稿共34页对于其中的流体有：对于其中的流体有：裂缝中流体裂缝中流体 岩块中流体岩块中流体 连续性方程连续性方程3对
18、于岩块系统和裂缝系统的连续性方程为：对于岩块系统和裂缝系统的连续性方程为：裂缝系统裂缝系统 岩块系统岩块系统 第27页，本讲稿共34页 将裂缝系统的运动方程、状态方程和窜流方程代入裂缝系统的连将裂缝系统的运动方程、状态方程和窜流方程代入裂缝系统的连续性方程中，有：续性方程中，有：即：即：同理可得：同理可得：以上两个式子为双重介质渗流的微分方程，通过简化，以上两个式子为双重介质渗流的微分方程，通过简化，可以求出简化模型解可以求出简化模型解。第28页，本讲稿共34页5.6 5.6 非牛顿流体渗流非牛顿流体渗流 流体的特性不服从粘滞阻力公式流体的特性不服从粘滞阻力公式：这样的流体都叫这样的流体都叫非
19、牛顿流体非牛顿流体。下面将对无惯性项的非牛顿流体的非线性渗流规律进行讨论下面将对无惯性项的非牛顿流体的非线性渗流规律进行讨论：1非牛顿流体的力学特性与类型非牛顿流体的力学特性与类型 流体在外力的作用下将发生变形和流动，这种变形和流动之间的关流体在外力的作用下将发生变形和流动，这种变形和流动之间的关系和特点，就叫系和特点，就叫流体的流变性流体的流变性，表达流体流变现象要用到流变状态，表达流体流变现象要用到流变状态方程，即方程，即本构方程本构方程。本构方程表达的主要内容是剪切应力与剪切速度之。本构方程表达的主要内容是剪切应力与剪切速度之间的关系。间的关系。第29页，本讲稿共34页 非牛顿流体可以分
20、为三大类，即非牛顿流体可以分为三大类，即纯粘性的非牛顿流体纯粘性的非牛顿流体、非稳态非牛顿流体非稳态非牛顿流体和和粘弹性流体粘弹性流体。纯粘性非牛顿流体，共有三种，其本构方程统一的形式为：纯粘性非牛顿流体，共有三种，其本构方程统一的形式为：非稳态非牛顿流体的剪切应力除了与剪切速率有关外还与作用非稳态非牛顿流体的剪切应力除了与剪切速率有关外还与作用时间有关，本构方程可以表示为：时间有关，本构方程可以表示为：对于粘弹性非牛顿流体，即具有液体的粘性，又具有固体的弹性特对于粘弹性非牛顿流体，即具有液体的粘性，又具有固体的弹性特征，并具有部分保持和恢复其原有形态的能力。对于这种流体，不仅要征，并具有部分
21、保持和恢复其原有形态的能力。对于这种流体，不仅要考虑剪切应力，而且要考虑其法向应力的作用。例如其中一种被称为马考虑剪切应力，而且要考虑其法向应力的作用。例如其中一种被称为马克斯维尔流体的本构方程为：克斯维尔流体的本构方程为：第30页，本讲稿共34页2塑性流体的一维渗流塑性流体的一维渗流 塑性流体是一种带有屈服性质的非牛顿流体。渗流方程为：塑性流体是一种带有屈服性质的非牛顿流体。渗流方程为：对于径向稳定流的情况，渗流方程为：对于径向稳定流的情况，渗流方程为：第31页，本讲稿共34页对于径向不稳定流，连续性方程的形式为：对于径向不稳定流，连续性方程的形式为：假设孔隙介质和流体均是微可压缩的，则：假
22、设孔隙介质和流体均是微可压缩的，则：这是建立非线性弹性渗流理论的基本公式这是建立非线性弹性渗流理论的基本公式 第32页，本讲稿共34页拟塑性液体的一维流动拟塑性液体的一维流动 3由流体力学的基础，可以得出拟塑性液体的本构方程是：由流体力学的基础，可以得出拟塑性液体的本构方程是：若若n1n1n1则为剪切增粘液体。则为剪切增粘液体。为幂为幂律系数律系数 为剪切速率为剪切速率 n n是幂律指数是幂律指数第33页，本讲稿共34页粘度定义的公式粘度定义的公式：视粘度公式：视粘度公式：渗流公式：渗流公式：K K为岩石渗透率，由于为岩石渗透率，由于液体的非牛顿性，渗液体的非牛顿性，渗透率透率K K与液体的性质密与液体的性质密切相关。切相关。第34页，本讲稿共34页