高等渗流力学(2015)-第五章
渗流力学 教学大纲
《渗流力学》课程教学大纲课程编号：02041002课程名称：渗流力学英文名称：Fluid Flow Through Porous Media课程类型：必修课课程性质：专业基础课总学时：56 讲课学时：48 实验学时：8学分： 4适用对象：石油工程专业、海洋油气工程、资源勘查工程先修课程：油层物理一、编写说明（一）制定大纲的依据根据《渗流力学》专业本科生培养计划要求制定本教学大纲。
（二）课程简介“渗流力学”是流体力学的一个分支，是研究流体在多孔介质中流动规律的一门学科。
本课程讲述的内容是“渗流力学”中的一个分支——地下渗流部分。
专门研究地下油气水及其混合物在地层中的流动规律。
（三）课程的地位和作用本课程是油气田开发与开采的理论基础，是石油工程专业和海洋油气工程专业的主干课程，同时也是资源勘查工程专业的选修课。
明确渗流理论是油气田开发，提高油田采收率等理论的基础，为学好专业课和解决有关地下油、气、水的渗流问题打好基础。
（四）课程性质、目的和任务本课程是石油工程专业和海洋油气工程专业本科学生的一门专业基础课，目的是通过各个教学环节使学生掌握油、气、水在地下流动规律，以及研究流体渗流规律的基本方法。
本课程的任务是使学生能掌握渗流力学基础概念、基本理论及解决渗流问题的基本技能。
（1）使学生掌握油、气、水渗流的基本规律及建立方程的基本方法；（2）培养学生用所学的渗流力学理论分析和解决渗流问题能力；（3）通过实验课培养学生严谨作风及动手能力。
（五）与其他课程的联系由于渗流力学是一门专业基础课，所以是其他专业课的基础，为学好其他专业课打下牢固的基础。
（六）对先修课的要求要求在学习本门课程之前，学好油层物理这门专业基础课，同时对高等数学中的求导，积分等知识能够熟练的应用。
一、大纲内容绪论渗流力学发展史，本课程研究方向。
第一章渗流的基础知识和基本定律（一）教学目的和任务使学生全面掌握渗流力学的基本概念和基本定律，使学生了解本课程的学习目的，为今后的学习打下基础。
高等渗流力学
浅谈非牛顿流体的渗流理论一．基本概念服从牛顿粘性定律的流体称为牛顿流体，所有气体和大多数液体都属于这一类。
水，酒精等大多数纯液体，轻质油，低分子化合物溶液以及低速流动的气体均是牛顿流体。
高分子聚合物的浓溶液和悬浮液一般是非牛顿流体。
从流体力学的角度，凡是服从牛顿内摩擦定律的流体称为牛顿流体，否则称为非牛顿流体。
所谓服从牛顿内摩擦定律是指在温度不变的情况下，随着流体梯度的变化， 值始终保持是常数。
度量液体粘滞性大小的物理量，简称为粘度。
物理意义是产生单位剪切速率所需要的剪切应力。
早在人类出现之前，非牛顿流体就已存在，因为绝大多数生物流体都属于现在所定义的非牛顿流体，而且非牛顿流体在化工方面宜属常见。
牛顿粘性定律的表达式为：（1-1）其中为牛顿粘度 为在剪切平面平行于流动平面的剪切应力，垂直于剪切平面的剪切速率。
二．非牛顿流体的分类下面是牛顿流体与非牛顿流体的流变图。
图牛顿流体与非牛顿流体的流变图根据流体的流变方程式，将非牛顿流体分类为：1.与时间无关的流体在流变图上来看对的曲线或是通过原点的曲线，或是不通过原点的直线，如图中b，c，d图线所示。
对于b，c这样的曲线来讲，斜率是变化的。
因此，对与时间无关的粘性流体来讲，粘度一词便失去了意义。
但是这些特定的曲线在任一特定点上都有一定的斜率，故与时间无关的粘性流体来讲，指在特定的剪切速率下，有一个表观粘度值。
即表观粘度是剪切速率的函数，不依赖时间的非牛顿液的流变特性只依赖于剪切应力的大小而不依赖于剪切应力的持续时间。
这样的流体可分为：（1）假塑性流体这种流体的表观粘度随剪切速率的增大而减小，其中的曲线关系为一下降的曲线，该曲线可用指数方程式表示：（1-2）大多数与时间无关的粘性流体都属于此类型，其中包括聚合物溶液，油脂，淀粉悬浮液，油漆等。
（2）涨塑性流体这种流体与假塑性流体相反，这种流体的表观粘度随剪切速率的增加而增大，其关系曲线为一上升的曲线。
如曲线c所示。
石油工程 提高采收率 第五章+渗流力学基础
（1）宾哈姆塑性体（塑性流体） 泥浆、油漆、稀润滑油等流体在静止时具有一个足以与某 一剪切应力相抗衡的三维结构，它们在施力后不会立即变形， 必须等到所施加的力增大到足以破坏三维结构，产生剪切变形 才开始流动。而在流动发生后，随剪切速率的增大，其剪切应 力又逐渐减小，最后接近于牛顿流体。这种流体称为宾哈姆塑 性体，其流变方程为： τ＝τ0+μpγ （2-22） 式中 τ0——屈服值或极限静切应力； μp——结构粘度。 为便于同牛顿流动定律（式2-7）相比较，也可以将式2-22写 成式如下形式： r0 p （2‒23）
3．径向流与单向流 水井注水时，进入井底的流体通过射孔眼后，就以井眼为中 心呈辐射状向四周发散，这种流动方式称为径向流。油井采油 时，来自地层各个方向的流体（原油和地层水）也是以径向流 的方式向井底汇集的。 与径向流不同，在实验室进行岩心流动实验时，流体是从 一维岩心的一端进入，从另一端流出，流线为彼此平行的直线, 这种流动方式称为直线流或单向流（图5-1）。
由于影响非牛顿流体的因素十分复杂，涉及到很多物理 和化学的微观机理，而这些机理有的尚难以弄清楚，所以 通常都采用实验的方法给出剪切应力与剪切速率之间的关 系曲线，这种关系曲线就称为流变曲线。然后以流变曲线 为基础再结合理论分析建立起不同类型非牛顿流体剪切应 力与剪切速率之间的函数关系式，这种函数关系式就称为 流体的流变方程。以下将对几种非牛顿流体的流变方程作 以简单描述。
而对于另外一些液体，当改变剪切速率时，其剪切应力与 剪切速率之间的关系与牛顿流动定律有很大偏离，也就是说它 们的粘度是随剪切速率而改变的。这种剪切应力与剪切速率的 关系不服从牛顿流动定律的流体，或者说粘度随剪切速率而改 变的流体称为非牛顿流体。 为了与牛顿流体粘度相区别，将这种随剪切速率而改变的 粘度称为表观粘度。高分子溶液、合成纤维溶液、纸浆、熔融 塑料、橡胶、涂料、乳化原油、液体石蜡、若干固体粒子的悬 浮液、乳状液以及奶油、蕃茄酱等都属于非牛顿流体。
渗流力学课后习题答案
渗流力学课后习题答案渗流力学课后习题答案渗流力学是研究地下水流动规律的一门学科，它在地质工程、水利工程等领域有着广泛的应用。
在学习渗流力学的过程中，习题是检验理论掌握程度和提高解题能力的重要方式。
下面将为大家提供一些渗流力学课后习题的答案，希望能对大家的学习有所帮助。
一、渗透率和渗透系数计算1. 计算渗透率时，需要知道渗透系数和介质的孔隙度。
渗透系数的单位是什么？如何计算渗透率？答：渗透系数的单位是米/秒。
渗透率的计算公式为：渗透率 = 渗透系数× 孔隙度。
2. 若一个土层的渗透率为1×10^-4 cm/s，孔隙度为0.4，求该土层的渗透系数。
答：渗透率的单位为cm/s，而渗透系数的单位为m/s。
所以需要将渗透率的单位转换为m/s。
1 cm = 0.01 m，所以渗透率为1×10^-6 m/s。
渗透系数 = 渗透率 / 孔隙度= (1×10^-6 m/s) / 0.4 = 2.5×10^-6 m/s。
二、多孔介质中的渗流1. 一个矩形土层，长为10 m，宽为5 m，渗透系数为1×10^-4 cm/s，上表面水头为10 m，下表面水头为5 m，求该土层的渗流速度。
答：渗流速度的计算公式为：渗流速度 = (上表面水头 - 下表面水头) × 渗透系数 / (土层厚度× 孔隙度)。
土层厚度为10 m，孔隙度未知，无法计算准确的渗流速度。
2. 一块长方形土层，长度为20 m，宽度为10 m，渗透系数为1×10^-3 cm/s，上表面水头为10 m，下表面水头为5 m，求该土层的渗流速度。
答：渗透系数的单位为cm/s，需要将其转换为m/s。
1 cm = 0.01 m，所以渗透系数为1×10^-5 m/s。
渗流速度 = (上表面水头 - 下表面水头) × 渗透系数 / (土层厚度× 孔隙度) = (10 m - 5 m) × (1×10^-5 m/s) / (20 m × 孔隙度) = 5×10^-6 / (20 × 孔隙度) m/s。
第五章-气体渗流理论
压力分布公式： 压力分布公式： 气体单向稳定渗流 液体单向稳定渗流
2
pe − pw p = pe − x L
2 2 2
pe − pw p = pe − x L
• 气体单向稳定渗流时，压力函数沿流程成 直线分布，而压力沿流程不成直线分布。 • 当x=常数时，p=常数，等压线是一族平行 于y轴的直线，由于流线和等压线互为正交， 所以流线是一族平行于x轴的直线。
k ∂p M ∂ p k ∂p ]=− [ ] µ ∂x RT ∂x Z µ ∂x k ∂p M ]=− µ ∂y RT k ∂p M ]=− µ ∂z RT ∂ p k ∂p [ ] ∂y Z µ ∂y ∂ p k ∂p [ ] ∂z Z µ ∂z
∂ (φρ ) ∂ pM =φ [ ] 右端 ∂t ∂t ZRT M 1 ∂p ∂ 1 =φ [ + p ( )] RT Z ∂t ∂t Z M 1 ∂p ∂ 1 ∂p =φ [ +p ( ) ] RT Z ∂t ∂p Z ∂t 1 M 1 ∂p ∂Z ∂p =φ [ + (− 2 ) p ] RT Z ∂t Z ∂p ∂t pM 1 1 ∂Z ∂p pM ∂p =φ [ − ] =φ C ( p) ZRT p Z ∂p ∂t ZRT ∂t
压力函数： 压力函数：
φµC ( p) k p ∂p k p ∇ ⋅[ ∇p ] = µ Z ∂t µZ k ~ = 2 k p dp + C ⇒ d~ = 2 k p dp p p ∫µ Z µZ
~ = 2 k p ∇p ⇒ ∇p µZ ∂~ p k p ∂p ⇒ =2 µ Z ∂t ∂t φµC ( p) ∂~ 1 ∂~ p p 2~ 渗流微分方程∇ p = = k ∂t η ∂t
渗流力学——油水两相渗流的理论基础
第五章油水两相渗流的理论基础
§3平面单相流等饱和度平面移动方程的应用
§4平面单相流两相混合带的压力
§5平面径向流等饱和度平面移动方程的应用
教学目的
及要求
1．掌握确定前缘含水饱和度和平均含水饱和度的方法
2．掌握确定排液道见水时间的方法
3．掌握平面单相流两相混合带的压力分布
4．掌握平面径向流等饱和度平面移动方程的推导
5．掌握平面径向流各个时刻地层内沿径向各点的饱和度分布及两相区的压力分布
教学内 容提要
1．平面单相流等饱和度平面移动方程的应用
确定前缘含水饱和度和平均含水饱和度
确定排液道见水时间的方法
2．平面单相流两相混合带的力
3．平面径向流等饱和度平面移动方程的应用
平面径向流等饱和度平面移动方程
平面径向流各个时刻地层内沿径向各点的饱和度分布及两相区的压力分布
第五章油水两相渗流的理论基础
周次
第6周，总第1次课
备注
章节名称
第五章油水两相渗流的理论基础
§1影响水驱油非活塞性的因素
§2等饱和度平面移动的基本微分方程
教学目的
及要求
1．了解影响水驱油非活塞性的因素
2．掌握等饱和度平面移动的基本微分方程建立过程
3．掌握分流方程式的推导
4．掌握饱和度分布公式的推导及曲线
教学重点、
难点及
重点：
确定前缘含水饱和度和平均含水饱和度的方法
平面径向流各个时刻地层内沿径向各点的饱和度分布
难点：确定前缘含水饱和度和平均含水饱和度的方法
处理方案及方法设计
画示意图讲解，举例计算说明，作业巩固理解
作业
练习
思考题：p90 5
渗流力学第五章演示文稿
第五章 单相微可压缩液体弹性不稳定渗流理论
不同内边界和边界条件下压力波的传播
供给边界，井底定产
供给边界，井底定压
封闭边界，井底定产
封闭边界，井底定压
第五章 单相微可压缩液体弹性不稳定渗流理论
第二节 无限大地层定产条件弹性不稳定渗流基本解
本节要点
1 渗流过程的数学描述—渗流数学模型； 2 渗流方程的求解方法和过程—基本解； 3 基本解的理解和应用
(
2p x2
2p y2
2p z2 )
p t
K Ct
导压系数，压力波向周围传播的速度，cm2/s。
Ct—综合压缩系数， Ct aCL C f
单位体积地层岩石在单位压力降下排除的液体体积，MPa-1 。
第五章 单相微可压缩液体弹性不稳定渗流理论
（一）渗流过程的数学描述
渗流数学模型
1 (r p) 1 p r r r t
r r r t
Байду номын сангаас 第五章 单相微可压缩液体弹性不稳定渗流理论
（二）数学模型的求解
? 所要求解的
1 (r p) 1 p r r r t 是一个二阶偏微分方程，不能直接求解。
解法：
解析解
分离变量法 积分变换法, etc
半解析解（试井分析） 数值解（数值模拟）
基本思想：偏—常；高—低
第五章 单相微可压缩液体弹性不稳定渗流理论
所求的解：
Q eu
p(r,t) po 4Kh
r2 4t
u
du
如果应用幂积分函数：
Ei (x)
eu du
xu
(x 0)
则我们的解可表达成：
Q
r2
p(r,t)
高等渗流力学(2017)-第五章-黄世军
p f
Kf
第五章 多重介质渗流理论
第三节 双重介质简化渗流模型的无限大地层典型解
一、Km和φf=0简化模型的典型解
在含油气裂缝-孔隙介质中，如果满足条件：
f m
Km K f
f 和 K f —是裂缝系统的孔隙度和渗透率; m 和 Km —是基岩系统的孔隙度和渗透率;
则在双重介质渗流的微分方程中，有两项可以忽略：
o K m q pm p f
q—单位时间单位岩石体积流出的流体质量；α—形状因子。
第二节 双重介质单相渗流的数学模型
假设孔隙介质，裂缝介质和地层流体均被认为是微可压缩 的，则裂缝孔隙压缩特性公式是：
f f 0 C f p f pi
基岩孔隙度 m压缩特性公式是：
第一节 双重介质油藏模型
该模型除与Warren-Root模型 相似，只是基质岩块不是平行 六面体，而是圆球体。圆球体 仍按规则的正交分布方式排列。
裂缝由圆球体之间的空隙表示，圆球体由基质岩块表示。
第一节 双重介质油藏模型
部分与整体以某种形式相似的形，称为分形。裂缝性 油藏的分形模型认为裂缝的分布形态、基岩的孔隙结构属 于分形系统。分形的维数随油藏的非均质性不同而不同。
Km 0
p f K f 1 p f K m pm p f f C f r t o r r r C pm K m p p 0 m m m f t
裂缝系统的压力公式变为：
r 2 Q p (r , t ) pi Ei 4 Kh 4 t
当η =0，即有充分的窜流时，渗流过程中的压力变化与单一介 质中的压力变化完全相同。
渗流力学_第五章
fw dfw dSw
dfw dSw
fw
dfw 1 t ∫xodx = dSw φA ∫0q(t)dt
x
¢ fw t x - xo = ò0 q(t )dt fA
Sw
¢ fw x - xo = W (t ) fA
Buckley—Leverett 方 程
Buckley曲线
Sw
t 1
t2 t3
t t2 1
∂pc ∂( po − pw ) ∂po ∂pw = = − ∂x ∂x ∂x ∂x
vt = vo + vw
∂pc vw − (vt − vw ) = − ∆ρ.g sin a Kw Ko ∂x
∂pc vw( + ) − vt = − ∆ρg sin a µw µo Ko ∂x Kw Ko
µw
µo
µo
′ fw
fw
0
0
Sw
பைடு நூலகம்
1
B
′ fw (Swf ) =
fw (Swf ) Swf − Swr
t
fw
A
x f − xo =
' fw (Swf )
∫0 q(t)dt φA
Swr
Swf Sw
Sw
∫0 q(t)dt Sw − Swr =
1 Sw − Swr = ′ fw(Swf )
t
φA(x f − xo )
油相: 水相:
∂(φρoSo ) − ∇⋅ (ρovo ) = ∂t ∂(φρwSw ) − ∇⋅ (ρwvw ) = ∂t
∂(vox ) ∂(voy ) ∂(voz ) ∂So + − + = φ ∂t ∂y ∂z ∂x ∂(vwx ) ∂(vwy ) ∂(vwz ) ∂Sw − + + = φ ∂t ∂y ∂z ∂x
高等渗流力学
浅谈非牛顿流体的渗流理论一．基本概念服从牛顿粘性定律的流体称为牛顿流体，所有气体和大多数液体都属于这一类。
水，酒精等大多数纯液体，轻质油，低分子化合物溶液以及低速流动的气体均是牛顿流体。
高分子聚合物的浓溶液和悬浮液一般是非牛顿流体。
从流体力学的角度，凡是服从牛顿内摩擦定律的流体称为牛顿流体，否则称为非牛顿流体。
所谓服从牛顿内摩擦定律是指在温度不变的情况下，随着流体梯度的变化， 值始终保持是常数。
度量液体粘滞性大小的物理量，简称为粘度。
物理意义是产生单位剪切速率所需要的剪切应力。
早在人类出现之前，非牛顿流体就已存在，因为绝大多数生物流体都属于现在所定义的非牛顿流体，而且非牛顿流体在化工方面宜属常见。
牛顿粘性定律的表达式为：（1-1）其中为牛顿粘度 为在剪切平面平行于流动平面的剪切应力，垂直于剪切平面的剪切速率。
二．非牛顿流体的分类下面是牛顿流体与非牛顿流体的流变图。
图牛顿流体与非牛顿流体的流变图根据流体的流变方程式，将非牛顿流体分类为：1.与时间无关的流体在流变图上来看对的曲线或是通过原点的曲线，或是不通过原点的直线，如图中b，c，d图线所示。
对于b，c这样的曲线来讲，斜率是变化的。
因此，对与时间无关的粘性流体来讲，粘度一词便失去了意义。
但是这些特定的曲线在任一特定点上都有一定的斜率，故与时间无关的粘性流体来讲，指在特定的剪切速率下，有一个表观粘度值。
即表观粘度是剪切速率的函数，不依赖时间的非牛顿液的流变特性只依赖于剪切应力的大小而不依赖于剪切应力的持续时间。
这样的流体可分为：（1）假塑性流体这种流体的表观粘度随剪切速率的增大而减小，其中的曲线关系为一下降的曲线，该曲线可用指数方程式表示：（1-2）大多数与时间无关的粘性流体都属于此类型，其中包括聚合物溶液，油脂，淀粉悬浮液，油漆等。
（2）涨塑性流体这种流体与假塑性流体相反，这种流体的表观粘度随剪切速率的增加而增大，其关系曲线为一上升的曲线。
如曲线c所示。
渗流力学要点整理
第一章 渗流力学基本概念和定律1、多孔介质（porous medium ）：含有大量任意分布的彼此连通的且形状各异、大小不一的孔隙的固体介质。
2、渗流（permeability ）：流体通过多孔介质的流动，也叫渗滤。
3、油藏：具有统一压力系统的油气聚集体4、渗流力学：研究流体在多孔介质中的运动形态和规律的科学。
5、油气层是油气储集的场所和流动空间6、定压边界油藏：层体延伸到地表，有边水供给区，在边界上保持一个恒定的压头。
7、封闭边界油藏：边界为断层或尖灭 没有边水供给 渗流中的力学分析及驱动类型：力学分析：重力、惯性力、粘滞力（大小用牛顿内摩擦定律表示1mPa·s ＝lcP ）、弹性力、毛管力。
驱动类型：依靠何种能量把原油驱入井底。
弹性驱动、水压驱动、溶解气驱、气压驱动（主要靠气顶气或注入气的膨胀能或压能驱油的驱动方式。
刚性气压驱动、弹性气压驱动）、重力驱动 不同驱动方式及开采特征总结：1、能量补充充足（边、底水，气顶、注水/气）：刚性驱动：刚性气/水驱；开采特征：Pe 、 Ql 、 Qo 有稳产段。
2、能量补充不充足（无边底水气顶注水注气或有而不足）： 弹性驱动：弹性驱动、溶解气驱、弹性气/水驱；开采特征：Pe 、 Ql 、 Qo 均不断下降。
3、 凡是气驱的Rp 都有上升的过程，其它驱动方式Rp 不变。
溶解气驱、刚/弹性气驱4、 Qo 或Rp 的突然变化反映水或气的突破。
供给压力Pe ：油藏中存在液源供给区时，在供给边缘上的压力。
井底压力Pw ：油井正常生产时，在生产井井底所测得的压力称为井底压力，也称为流动压力，简称流压。
折算压力Pr ：油藏中某点折算到某一基准面时的压力，它表示油层中各点流体所具有的总能量。
达西定律：在一定范围内△P 与Q 成直线关系，当流量不断增大，直线关系就会被破坏。
真实流速与渗流速度的关系达西定律适用条件： 液流处于低速、层流，粘滞力占主导地位，惯性主力很小，可忽略。
渗流力学
绪论：1.渗流力学：就是研究渗滤的运动状态和运动规律的学科。
渗流力学研究涉及三个主要方面：工程渗流、生物渗流、地下渗流2.渗流：流体通过多孔介质的流动称为渗流或渗滤3.多孔介质：由骨架和相互连通的孔隙、裂缝、溶洞或各类毛细管体系组成的材料 第一章：1.油气藏：油气的储集的场所和流动空间油气藏作用：限制流体的流动范围、影响流体的渗流心态、决定流体的边界形状 按圈闭条件分为：①构造油气藏（背斜油气藏、断层油气藏、刺穿接触油气藏）； ②地层油气藏（潜山油气藏、生物礁油气藏、不整合覆盖油气藏、地层超覆油气藏）；③岩性油气藏（透镜状岩性油气藏、尖灭性岩性油气藏）根据流体在其中流动的空间特点分为：①层状油藏；②块状油藏2.多孔介质的特点：具有孔隙性、渗透性、比表面积大、孔隙结构复杂等基本特点 绝对渗透率：岩石允许流体通过的能力 有效渗透率：（相渗透率）：岩石对于某一相流体的通过能力 相对渗透率：有效渗透率与绝对渗透率的比值按结构分类(结构复杂性)：1.粒间孔隙结构；2.纯裂缝结构；3.裂缝-孔隙结构；4.溶洞-孔隙结构；5.溶洞-裂缝-孔隙结构 3.连续流体：把流体中的质点看成是在一个很小的体积中包含着很多分子的集合体，质点中流体的性质与周围质点中的流体性质成连续函数关系 连续介质：是在质点的典型体积上表现出来的平均性质连续介质场：连续流体在连续介质中的流动，在研究其流动规律时，其物性是连续变化的，即其数学方程是连续的，在这种连续系统中流动的场4.渗流过程中的力：重力、惯性力、粘滞力、弹性力、毛管力5.油藏中的压力：原始地层压力、供给压力、井底压力、折算压力(计算P19）6.油藏的驱动类型：重力水压驱动、弹性驱动、气压驱动、溶解气驱、重力驱动7.※达西定律8.渗流速度：渗流量与渗流截面积之比9.真实速度：渗流量与渗流截面的空隙面积之比10.渗流的基本方式：单相流、平面径向流、球面向心流11.非线性渗流指数形式：v=C （dp/dL)^n 式中C 为取决于岩层和流体性质的系数； n 为渗流指数 ， n є（0.5～1）， n=1时，渗流服从达西直线定律 12.启动压力梯度（吸附膜和水化膜的影响）：在压力梯度较小时，流体不产生流动，渗流速度为零，当压力梯度大于某一值后，流体才发生流动，这一压力梯度值称为启动压力梯度 13.两相流体时，渗流阻力明显增加，且两相各自渗透率之和不等于单相渗流时的绝对渗透率。
渗流力学
渗流力学绪论多孔介质：由固体骨架和相互连通的孔隙，裂缝，溶洞或各种类型的毛细管体系所组成的材料。
渗流力学与其他力学的区别：介质的不同。
第一章渗流的基本概念和基本规律油气藏：油气储集的场所和流动的空间。
油气藏按圈闭形成的类型：构造油气藏，地层油气藏，岩性油气藏。
构造油气藏的分类：背斜油气藏，断层油气藏，刺穿接触油气藏。
油气藏根据流体流动空间的特点：层状隐藏，块状油藏。
层状油藏的特点：1：油层平缓，分布面积大。
2：多油层，多旋回。
3：只考虑在水平方向上流动的流体。
块状油气藏得特点：有限的圈闭面积内相当厚的油藏，考虑纵向上流体的流动和交换；考虑毛管力和重力的作用。
纵向上分为三个区：纯油区，过渡区，纯水区。
过渡区：含束缚水过渡带，油水同生过渡带，残余油过渡带。
多孔介质的特点：孔隙性，渗透性，比表面积大及孔隙结构复杂。
渗透性：多孔介质允许流体通过的能力。
K= ；渗流：流体在多孔介质中的流动。
绝对渗透率：当岩石中的孔隙流体为一项时，岩石允许流体通过的能力。
有效渗透率：当岩石中有两种以上流体存在时，岩石桂其中一相的通过的能力。
相对渗透率：岩石的有效渗透率与绝对渗透率的比值。
比表面积：单位体积岩石所有岩石颗粒的总表面积或孔隙内表面积。
孔隙类型：粒间孔隙，裂缝，溶洞。
多孔介质巨大的比面和复杂的孔隙结构，使得渗流具有阻力大，流动速度慢的特点。
油气层孔隙结构分为：单纯介质（粒间孔隙结构和纯裂缝结构），双重介质（裂缝-孔隙结构和溶洞-孔隙结构），三重介质（大洞或大裂缝和微裂缝、微孔隙共生）。
理想结构模型：将岩石的孔隙空间看成是由一束等直径的微毛细管组成。
修正理想结构模型：变截面弯曲毛细管模型。
重力（动力或阻力），惯性力（阻力），粘滞力（阻力），弹性力（动力），毛管力（动阻力）原始地层压力：油藏开发前流体所受的压力。
供给压力：油藏中存在液源供给区时，在供给边缘上的压力。
井底压力：油井正常工作时，在生产井井底所测得的压力。
第一篇 第五章 渗流力学基础
188 第五章 渗流力学基础第一节 油气层渗流的达西定律油气层渗流的基本规律是达西定律。
1856年法国水利工程师达西在研究城市供水问题时，欲测得获得一定的流量需要消耗的能量，于是达西运用填满砂的管子做实验，得到了水流速与管截面积、入口与出口压头之间的关系式，后人为纪念他，将这一定律称达西定律。
一、达西实验及结果达西实验装置如图1—ll 所示，液体经过进水管a 进入模型主体。
再透过砂层，经节流阀门流入量杯。
节流阀可以控制流速，量杯D 测取流量Q 。
测压管可以分别测出过水断面1-1，2-2上的压力p 1、p 2。
稳压管b 可以使模型内液面稳定在b 管的位置上。
显然，节流阀开度不同时，将得到不同的流量和不同的测压管高度。
实验结果发现：流量大小与管于截面积A ；入口及出口压力差p 2-p 1成正比，与填满砂粒的管子长度△L 成反比，将上述关系写成等式，需加上比例系数K 。
即：Lp p KAQ ∆-=12 (1—6) 式中 K ——渗透率，它表征多孔介质和液体的渗透能力。
二、达西定律的导出(一)由管路水力学导出达西定律由普通水力学可知，任意过水断面上的总能量表示成下列形式：gv pZ H 22++=γ （1-7）式中 H-——总水头；Z ——位置水头；γp---压力水头； gv 22---流速水头。
189由于渗流速度v 很小，可以忽略gv 22项，于是总水头可表示为：γpH =+Z （1-8）断面1—1，2—2上的水头差可表为：⎪⎪⎭⎫ ⎝⎛--⎪⎪⎭⎫ ⎝⎛+=-=∆H γγ221121p Z p Z H H （1-9）达西通过实验发现：通过砂层的流量Q 与水头损失△H 成正比，与渗流面积A 成正比，与渗流段长度△L 成反比。
即：LHAQ ∆∆∞= （1-10）欲将(1—10)写成等式需加一比例常数，于是我们得到：LHAK Q i ∆∆= （1-11） 式中 K i -——比例常数，称为渗流系数，它与流体及砂层的性质有关。
渗流力学课程大纲
开课编号:323004Z渗流力学Mechanics of Fluid Flow in Porous Media课程编号：323004Z课程属性：专业基础课学时/学分：40/2预修课程：偏微分方程，数值计算方法，弹性力学有限单元法教学目的和要求：系统掌握渗流的基本知识，掌握渗流分析理论，了解渗流数值方法，学会利用渗流计算有限元软件，分析渗流问题。
内容提要：第一章：张量的基本知识（1）张量的指标符号表示法第二章：渗流的基本概念与定律（1）多孔介质、连续介质假设、连续流体、连续介质场等基本概念；（2）流体的实际平均速度与渗流速度（3）达西定律（定律及其适用范围、定律在多相多维渗流中的推广）与非线性运动方程（低速非线性、高速非线性、低渗介质非线性、非牛顿流体渗流）（4）岩土材料的基本物理特性（5）岩土材料的非饱和渗流特性第三章：渗透变形与渗流破坏（1）渗透变形的类型（2）渗透变形的判别方法（3）渗透变形的防护第四章：流体的偏微分方程与定解条件（1）单相渗流连续性方程（2）两相不溶混渗流连续性方程（3）流体与骨架的状态方程（4）单相流体的偏微分方程（5）定解条件（初始条件与边界条件）第五章：渗流的理论计算方法（1）Dupuit假定及其应用（2）井的渗流计算方法第六章：渗流的数值方法（1）有限差分法介绍（2）有限单元方法介绍（3）渗流的边界条件及其处理方法（4）非线性方程组的求解方法第七章：渗流与岩土体变形的相互作用（1）岩土体的应力变形及其有限元方法（2）渗流与岩土体变形的耦合相互作用（3）岩土体的水力劈裂第八章：渗流计算软件与工程分析（1）利科渗流分析软件介绍（2）利科软件的应用实例第九章：岩土介质中溶质的输运教材：无主要参考书：1 毛昶熙，渗流计算分析与控制，中国水利出版社；2 孔祥言，高等渗流力学，中国科学技术大学出版社；3 刘杰，土的渗透稳定与渗流控制，水利电力出版社撰写人：吴梦喜(中国科学院力学所)撰写日期：2010年7月。
渗流力学
20
《渗流力学》绪论__八、几个基本概念
４、刚性和弹性：
５、稳定渗流和不稳定渗流：
21
《渗流力学》绪论__八、几个基本概念
6、牛顿和非牛顿：思考题 7、达西和非达西：思考题 8、线性和非线性：思考题
22
《渗流力学》绪论__八、几个基本概念
9、点源和点汇： 点源：注水井——产量为负（-Q或-q） 点汇：生产井——产量为正（+Q或+q）
广度：多孔介质理论，表面物理固体力学等
12
《渗流力学》绪论
六、渗流力学发展方向特点
3、研究渗流力学带有根本性变化的新方向，油气层连续介质场变为 随机变量场，概率渗流力学，随机场渗流力学： A：微扰法——地层参数分成有规律和随机部分 B：斯克天尔左右法——随机的 C：实用计算机——正态分布显示 油田开发信息处理：渗流过程与自动控制理论结合，对油气田 开发过程进行控制，分析，调整，进行自动处理。控制油田开采过 程的信息量，处理开发信息，定量估计开发状况，评价开发方案的 效果。建立起描述开采过程的统计模型，对油气田开发进行自动调 整和预报。用数模和最优化理论制定合理的油气井工作制度，把渗 流过程的研究与自动化管理油气田生产结合起来。
H
式中，系数b称为压力系
数。常规油气藏b的取值为
0.7～1.2。当b<0.7时称为异 常低压油藏；当b>1.2时，称 为异常高压油藏。
0
p=a+bH
p
油藏埋深H 和实测压力p 曲线
第一章
各种压力的概念
渗流基本规律及渗流模型
2. 目前地层压力p : 油藏开发过程中，不同时期的地层
压力称为“目前地层压力”。
3
《渗流力学》绪论
5、渗流力学-第五章
第五章 单相微可压缩液体弹性不稳定渗流理论
（1）无限大地层中等产量的两口生产井
这样选择坐标系，使两井井点 坐标分别为(a，0)和(-a，0)。此时 地层中任一点处任一时刻表达式可 写成：
po
−
p( x,
y, t )
=
Qμ 4πKh
{−
Ei
[−
(x
+
a)2 + 4ηt
y2
]−
Ei[−
(x
−
a)2 + 4ηt
po
−
Qμ 4πKh
{−
Ei
[−
(x − a)2 +
4ηt
y2 ]}
+
Qμ 4πKh
{−
E
i
[−
(x
+
a)2
4ηt
+
y2
]}
( x > 0)
22
第五章 单相微可压缩液体弹性不稳定渗流理论
3. 井以变产量生产问题
如果无限大地层中有一口生产井，自某时刻开始生产，从t=0 到t=t1，产量为Q1，从t=t1以后，产量为Q2，现在求地层压力分布。
y2
]}
在x=0处的压力：
∂ ∂
p x
x=0
=
−
Qμ 2a 4πKh [ a2 + y2
−a2+ y2
e 4ηt
+
− 2a a2 + y2
−a2+ y2
e 4ηt
]=
0
因而在任何时刻，坐标轴x = 0都是一条流线，即y轴是一条流线。
19
第五章 单相微可压缩液体弹性不稳定渗流理论
高等渗流力学作业-完整版
（1）建立描述渗流问题的数学模型； （2）求任意时刻井底压力的表达式（吐哈美原理） ； （3）如果油井在生产 T 后开始关井，试求关井后任意时刻的井底压力表达 式。 第三章 两相渗流
作业 1、(2014-3-7) —— 黄世军 一维直线油水两相渗流，油水井排距为 L。不考虑毛管压力和重力，地层 宽度为 B，厚度为 h，渗透率为 K，原油粘度为μo，水的年度为μw，室内实验 测得油水相渗曲线 Kro(Sw)和 Krw(Sw)，以恒定注入速度 Qw 注入； 分析： （1）见水前某一时刻 t0，a：油水前缘位置；b:油水前缘饱和度；c：原油 采出程度。 （2）见水时刻 T0，无水采油期，无水采油期采出程度，注水利用率。 （3）见水后，任一时刻 t1，a:油井含水率；b：原油采出程度。 作业 2、(2014-3-7) —— 黄世军 毛管渗吸的物理实质是什么？对油气田开发有何影响。 作业 3、(2014-3-7) —— 黄世军
作业 1、 （2014-3-14）—— 黄世军 非等温渗流数学模型的特殊性表现在什么方面？ 作业 2、 （2014-3-14）—— 黄世军 稠油热采方式常见有哪几种？试阐述其渗流原理。 作业 3、 （2014-3-14）—— 黄世军 稠油热采过程中，由于热量注入地层，对近井渗流会造成什么样的影响？ 第六章 1（2014-3-14）—程林松 与常规油藏相比， 低渗特低渗油藏特征常规油藏的差别，建立渗流数学模型 时如何考虑？ 第十章 复杂结构井
作业 1、 （2014-3-20）—— 黄世军 复杂结构井的渗流优势及适用油气藏。 作业 2、 （2014-3-20）—— 黄世军 复杂结构井非线性渗流体现在什么方面。 作业 3、 （2014-3-20）—— 黄世军 复杂结构井近井渗流理论发展的阶段特征。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Km f C f div(grad p f ) ( pm p f ) 0 t
Kf
忽略
p f
pm K m Km mCm div(grad pm ) ( pm p f ) 0 t
第三节 双重介质简化渗流模型的无限大地层典型解 上面两式化简为：
C f C f C f 0
Cm m m0 o C t m0
p p C m0 o m t t
C m Cm C m0
第二节 双重介质单相渗流的数学模型
C m m m 0 o 1 C p pi m 0
由此得到上二式对时间的导数：
C f f f 0 o C t f 0 p p f 0 oC f t t
一 运动方程
认为达西线性流公式对裂缝的基岩均是适用的，则有如下 渗流速度公式 ： 裂缝系统： 基岩系统：
vf
vm
Kf
Km
gradp f
gradpm
第二节 双重介质单相渗流的数学模型
二、窜流方程
在基岩与裂缝之间存在着压力差异，因而存在流体交换，但 这种流体交换进行是较缓慢，可将其视为稳定过程。一般认为单位
第一节 双重介质油藏模型
第二节 双重介质单相渗流的数学模型 第三节 双重介质简化渗流模型的无限大地层典型解
第四节 裂缝-孔隙介质中两相渗流理论
第五节 双重介质油藏试井理论分析基础
第二节 双重介质单相渗流的数学模型
建立双重介质油藏的数学模型时，两种介质分别满足各自 的运动方程、状态方程和连续性方程，而两种连续介质间窜流 通过连续性方程中的一个源和汇函数来表示。
Kf
Km div(grad p f ) ( pm p f ) 0 pm K m mCm ( pm p f ) 0 t
(1)
(2)
对(1)式求导带入(2)式并消去压差pm-pf得裂缝系统压力变
化的偏微分方程：
Co div[ grad p f Co grad p f ] 0 t t Kf
p f 1 p f r t t r r r
井筒
1 1 p f r r r r
井筒
(4)
基质岩块系统系统
Kf
裂缝系统
Km
Warren-Root模型示意图
第三节 双重介质简化渗流模型的无限大地层典型解
时间内从基岩排至裂缝中的流量与以下因素有关：
(1) 流体粘度； (2) 孔隙和裂缝之间的压差； (3) 基岩团块的特征量，如长度、面积和体积等； (4) 基岩的渗透率。 通过分析可以得出窜流速度q为：
o K m q pm p f
q—单位时间单位岩石体积流出的流体质量；α—形状因子。
——形状因子。
窜流系数的大小，既取决于基质和裂缝渗透率的比值，又 取决于基质被裂缝切割的程度，基质与裂缝渗透率的比值越 大或者裂缝密度越大，窜流系数越大。
第一节 双重介质油藏模型
3.形状因子
形状因子α与基质岩块大小和正交裂缝组数有关，岩块越 小，裂缝密度越大，形状因子α越大 。
Warren-Root提出的计算α的关系式：
第二节 双重介质单相渗流的数学模型
三、状态方程
假设孔隙介质，裂缝介质和地层流体均被认为是微可压缩 的，则裂缝孔隙压缩特性公式是：
f f 0 C f p f pi
基岩孔隙度 m压缩特性公式是：
m m0 C m pm pi
对于其中的流体（如原油）则：
f 、 m ——裂缝和基质系统的孔隙度。
第一节 双重介质油藏模型
2.窜流系数
流体在双重介质油藏渗流的过程中，基质与裂缝之间存在着 流体交换。窜流系数就是用来描述这种介质间流体交换的物理量， 它反映了基质中流体向裂缝窜流的能力。窜流系数定义为：
Km 2 rw Kf
K f ，K m ——裂缝系统和基质系统的渗透率， m2；
四、连续性方程
裂缝系统
f div v f q 0 t
m div vm q 0 t
基岩系统
对于均质各向同性地层，上式中的对流项可以化简为：
Kf div v f o div grad p f
Km div vm o div grad pm
第二节 双重介质单相渗流的数学模型
经过处理后，连续性方程变为：
Km f Cf div(grad p f ) ( pm p f ) 0 t pm K m Km mCm div(grad pm ) ( pm p f ) 0 t
双重介质
裂缝
在一般情况下，裂缝所占的储集空间
大大小于基岩的储集空间，因此裂缝 孔隙度就小于基岩的孔隙度，而裂缝
的流油能力却大大高于基岩的流油能
力，因此裂缝渗透率就高于基岩的渗 透率，这种流油能力和供油能力的错
裂缝 基岩 双重介质实际油藏模型
位的现象是裂缝-孔隙介质的基本特性。
第一节 双重介质油藏模型
第五节 双重介质油藏试井理论分析基础
第三节 双重介质简化渗流模型的无限大地层典型解
一、Km和φf=0简化模型的典型解
在含油气裂缝-孔隙介质中，如果满足条件： Km K f f m
f 和 K f —是裂缝系统的孔隙度和渗透率; m 和 K m —是基岩系统的孔隙度和渗透率;
则在双重介质渗流的微分方程中，有两项可以忽略：
高等渗流力学
黄世军 2014
第五章 多重介质渗流理论
第一节 双重介质油藏模型 第二节 双重介质单相渗流的数学模型 第三节 双重介质简化渗流模型的无限大地层典型解 第四节 裂缝-孔隙介质中两相渗流理论
第五节 双重介质油藏试井理论分析基础
第一节 双重介质油藏模型
双重介质定义
具有裂缝和孔隙双重储油(气)和流 油(气)的介质我们称之为双重介质。 基岩
对于裂缝均匀分布、基质具有较高的窜流能力和高储存
能力的条件下，其结果与Warren-Root模型的结果相似。
第一节 双重介质油藏模型
3.De Swaan模型
该模型除与Warren-Root模型 相似，只是基质岩块不是平行 六面体，而是圆球体。圆球体 仍按规则的正交分布方式排列。
基质
裂缝
裂缝由圆球体之间的空隙表示，圆球体由基质岩块表示。
初始条件： t 0,p f (r , t ) |t 0 pi
p f lim r 内边界条件： r 0 r p f Q r 2 K f h t r
外边界条件： r
lim p f r , t pi
2 式中： Co mCm , K f / ( Km ) rw /
p f
(3)
第三节 双重介质简化渗流模型的无限大地层典型解 如果η→0，窜流速度加快，地层流体可以很快的由基岩流 入裂缝，然后按照裂缝系统渗流规律流动。此时方程(3)退化为 单纯裂缝介质不稳定特性渗流方程。
它相当于一个连续性方程,其中的渗流速度由两部分组成， 第一部分是纯裂缝中的渗流速度，第二部分是窜流速度引起的 附加渗流速度，即： Kf v grad p f C0 grad p f t
1.弹性储容比
弹性储容比ω定义为裂缝系统的弹性储存能力与油藏总 的弹性储存能力之比，用来描述裂缝系统和基质系统的弹性 储容能力的相对大小。
f Cf f C f mCm
基质系统孔隙体积 裂缝系统孔隙体积 m f 基质和裂缝系统总体积 基质和裂缝系统总体积
Cf 、 Cm ——裂缝和基质系统的综合压缩系数；
o 1 C p pi
第二节 双重介质单相渗流的数学模型
渗流问题中常遇到乘积 f 和 m 的压缩特性。由于介 质和流体的微可压缩性，舍去高阶无穷小量后可得到：
C f f f 0 o 1 C f 0 p pi
其中： K f / 0C0 K f / m 0Cm 是导压系数。
注意到： t 0 可得新的边界条件为：
p f lim r r 0 r 0
p f Q t / lim r 1 e r 0 2 K f h r
第一节 双重介质油藏模型
4.Factal模型
分形模型:整 体与局部具
基质
裂缝
有某种相似 性
部分与整体以某种形式相似的形，称为分形。裂缝性 油藏的分形模型认为裂缝的分布形态、基岩的孔隙结构属 于分形系统。分形的维数随油藏的非均质性不同而不同。
第一节 双重介质油藏模型
双重介质油藏基本参数：弹性储容比和窜流系数。
裂缝-孔隙性双重介质结构油藏可抽象地简化成各 种不同地质模型。
1.Warren Root 2.Kazemi 3.De Swaan 4.Factal
模型 模型
模型
模型
第一节 双重介质油藏模型
1.Warren-Root模型
将双重介质油藏简化为正交裂缝 切割基质岩块呈六面体的地质模 型，裂缝方向与主渗透率方向一 致，并假设裂缝的宽度为常数。
裂缝
基质
裂缝网络可以是均匀分布，也可以是非均匀分布的，采 用非均匀的裂缝网格可研究裂缝网络的各向异性或在某一方 向上变化的情况。
第一节 双重介质油藏模型
2.Kazemi模型
该模型是把实际的双重介 质油藏简化为由一组平行层理 的裂缝分割基质岩块呈层状的