高等渗流力学(2017)-第五章-黄世军
高等渗流力学(2017)-第四章

高等渗流力学
黄世军
2017
第四章 孔隙介质中的多相多组分渗流理论
由于多相多分组系统是一个很复杂的物理化学系统,因此无论 在对系统本身的物理化学性质的研究还是对于流动规律的研究, 包括对物理化学过程的描述和流动规律的描述,都遇到极为困难 的问题。即使有可能建立起基本微分方程,其求解也是相当困难
第一节 多相多组分渗流数学模型 渗流数学模型解法思路
1、总物质守恒: L V Lw 1 (1个) (Nc-1个)
2、某一烃组分守恒:
Lxi Vyi 1,(i 2、 3...NC )
3、二氧化碳组分守恒: Lx1 Vy1 Lwn1,w 1 4、相平衡:
fi , L fi ,V i 1、 2...NC
7、选取未知量:
Y V , y1 , y2 ... yN
C
Fi fi ,V fi , L NC 8、构造牛顿迭代方程组,余量形式: FNC 1 1 yi i 1
第一节 多相多组分渗流数学模型 渗流数学模型解法思路
9、构造迭代式:
J Y F
(6)
由(6)和(7)可写出Nc+1个方程组成的方程组。 利用Newton-Raphson方法求解。
第二节 相态平衡闪蒸计算方法
一、一般相态平衡闪蒸计算方法
迭代求解过程:
Newton-Raphson方法求解要点是形成Jacobi矩阵元素:
水侵量--高等渗流力学

油藏天然水侵量的计算1问题简述图1 油藏示意图如图1所示,圆形油藏周围有边水,在油藏生产过程中,油藏平均压力是变化的,求解在这种情况下的累积水侵量。
2技术思路将油藏看成一口井,但这口井的井底压力(即油藏的平均压力)是不断变化的。
因此,该问题是变压力条件下求产量的定解问题。
先求得定压力下的产量解,然后根据杜哈美原理求得变压力下的累积产量,即求得了油藏的累积水侵量。
3理论公式推导第一步:油藏边界σ上压力01p p σ-=为常数时的压降解设油藏边界σ上压力01p p σ-=为常数时的解(初始压力为0p )为(,)D D p p r t ∆=∆,其无因次数学模型为:22101,1|1(,)|0(,)0()(,)=0()lim D D D D Dw D Dw D D D Dr D D t D D r D D Dr r p p pr r r r r t p p r t p r t p r t r ==→∞=⎧∂∆∂∆∂∆+=≤≤⎪∂∂∂⎪⎪∆=⎪∆=⎪⎨∆=⎪⎪⎪∂∆⎪⎪∂⎩(内边界条件)(初始条件)无限大边水区域有限封闭边水区域 (1)其中,无因次半径:D err R =(2)无因次时间:2D t e ktt c R φμ=(3)对上述数学模型进行Laplace 变换,并令:(,)(,)s D D p r s e p r d τττ∞-∆=∆⎰(4)先对渗流方程进行变换:22001()stst D D D D p p p e dt e dt r r r t ∞∞--∂∆∂∆∂∆+=∂∂∂⎰⎰ 得到:22(,)(,)1(,)D D D D D Dd p r s d p r s s p r s dr r dr ∆∆+=⋅∆(5)对内边界条件进行变换,得到:11(,)|D D r p r s s=∆=(6)对无限大边水区域情况的外边界条件进行变换,得到:(,)0lim D Dr p r s →∞∆=(7)对有限封闭边水区域情况的外边界条件进行变换,得到:(,)0D DwD Drr d p r s dr =∆= (8)对(5),两边同乘2D r ,然后第一项分子分母同乘2,:)())()22222(,)0(,)0D DDDD Dd p r s p p r s d∆∆++∆=(9)为0阶虚宗量的Bessel 方程,通解为:00(,)))D D D p r s AI BK ∆=+(10)第一种外边界条件:无限大边水区域的情况有(7),且,由0()I x 的渐近性知,当D r →∞,0)D I →∞,得到0A =,得到:0(,))D D p r s BK ∆=(11)将内边界条件(6)带入(11)得到:(,)D p r s ∆=(12)第二种外边界条件:有限封闭边水区域的情况 将内边界条件(6)和外边界条件(8)带入(10),得到:110101)),))Dw Dw Dw Dw K I A B C C C s K I I K ==⎡⎤=+⎣⎦第二步:水侵量与压力差的关系 由达西定律:2()er R k p q t r r πμ=∂∆⎛⎫=-⎪∂⎝⎭ (13)物理意义为:单位压差、单位厚度条件下的水侵速度。
高等渗流力学(2015)-第五章-黄世军

o K m q pm p f
q—单位时间单位岩石体积流出的流体质量;α—形状因子。
第二节 双重介质单相渗流的数学模型
三、状态方程
假设孔隙介质,裂缝介质和地层流体均被认为是微可压缩 的,则裂缝孔隙压缩特性公式是:
f f 0 C f p f pi
基岩孔隙度 m压缩特性公式是:
f 和 K f —是裂缝系统的孔隙度和渗透率; m 和 K m —是基岩系统的孔隙度和渗透率;
则在双重介质渗流的微分方程中,有两项可以忽略:
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
它相当于一个连续性方程,其中的渗流速度由两部分组成, 第一部分是纯裂缝中的渗流速度,第二部分是窜流速度引起的 附加渗流速度,即: Kf v grad p f C0 grad p f t
t
分析公式: Co
p f
div[
Kf
grad p f Co
grad p f ] 0 t
第八节 三重介质油藏试井理论分析基础
第一节 双重介质油藏模型
双重介质定义
具有裂缝和孔隙双重储油(气)和流 油(气)的介质我们称之为双重介质。 基岩
双重介质
裂缝
在一般情况下,裂缝所占的储集空间
大大小于基岩的储集空间,因此裂缝 孔隙度就小于基岩的孔隙度,而裂缝
的流油能力却大大高于基岩的流油能
力,因此裂缝渗透率就高于基岩的渗 透率,这种流油能力和供油能力的错
Kf
这就是考虑双重孔隙性和双重渗透性的双重介质渗流的微
《渗流力学》课程教学大纲衡帅

《渗流力学》课程教学大纲课程编号: 020092080总学时及其分配:48;理论课48学时。
学分数:3适用专业:采矿工程专业(煤及煤层气工程方向)任课学院、系部:能源学院一、课程简介“渗流力学”是流体力学的一个分支,是研究流体在多孔介质中流动规律的一门学科。
本课程讲述的内容是“渗流力学”中的一个分支-地下渗流部分。
专门研究地下水、气及其混合物在地层中的流动规律。
本课程是油气田勘探与开采的理论基础,是采矿工程专业煤层气开采方向的专业基础课和核心课,同时也是地质勘查专业、安全工程专业的专业基础课。
学习该课程的目的,是把它作为认识煤层气藏、改造煤层气储层藏的工具,作为煤层气储层开发设计、动态分析、煤层气井开采、增产工艺、反求地层参数、提高采收率等的理论基础。
因此,它是采矿工程煤层气方向的主干专业基础课之一,是进一步学好《采气工程》、《煤层气试井理论与技术》的关键课之一,该门课的目的是让学生了解油气在储层中渗流的基本规律及研究油气在储层中渗流的基本方法。
二、课程教学的目标本课程是采矿工程专业煤层气方向本科学生的一门专业基础课和核心课,目的是通过各个教学环节使学生掌握油、气、水在地下流动规律,以及研究流体渗流规律的基本方法。
明确这些理论是油气田开发、煤层气开发、瓦斯抽采,提高油气采收率、瓦斯抽采效等的理论基础,为学好专业课和解决有关地下油、气、水的渗流问题打好基础。
本课程的任务是使学生能掌握渗流力学基础概念、基本理论及解决渗流问题的基本技能。
三、课程教学的基本内容及教学安排1. 绪论(3学时)知识要点:渗流力学的基本任务、研究方法、相关研究方向和国内外发展动态及一些与本学科相关的基本概念。
目标要求:了解渗流力学在石油天然气开发、瓦斯抽采、地下水流动中的重要性,知道如何学好渗流力学。
2. 第一章渗流的基本概念和基本规律(8学时)知识要点:知道油气藏及其简化,多孔介质及连续介质,渗流过程中的力学分析及驱动类型,渗流的基本规律和渗流方式,非线性渗流规律,低速下的渗流规律,两相渗流规律等基本概念。
渗流力学(课件)渗流力学(课件)渗流力学-习题集精选全文

可编辑修改精选全文完整版第一章 油气渗流基本定律和渗流数学模型一、基本概念1、何谓多孔介质?在油气层中,分哪几类?2、什么叫渗流、渗流力学、油气层渗流研究对象是什么?3、现阶段油气渗流力学的研究特征是什么?4、什么叫含油边缘和计算含油边缘?5、何为开敞式和封闭式油藏?区别是什么?6、什么叫折算压力?怎样求地层中某一点折算压力?7、什么叫地层压力系数和压力梯度曲线?8、常见的驱油能量有哪些?有哪些最基本驱动方式?9、何为渗流速度?为什么要引入它?它与流体质点的真实速度的区别何在? 10、什么叫线性渗流定律、其物理意义是什么?怎样确定其适用范围? 11、岩石渗透率的物理意义和单位是什么?各种单位制之间有什么联系? 12、何谓非线性渗流的指数式?其物理意义是什么?13、何谓非线性渗流的二项式?其物理意义是什么?它与指数式有何区别和联系? 14、什么叫流压和静压?15、什么叫渗流数学模型?其一般构成是什么?16、建立渗流微分方程应从哪几个方面考虑?分几个步骤进行?17、简述分别用积分法和微分法推导单相流体稳定渗流微分方程的步骤? 18、分别写出液体、气体和岩石的状态方程。
二、计算题1、有一未打开油层,如图:其中P A =18MPa,h=10m,原油重度γ=0.8,求P B =?2、四口油井的测压资料如下表,已知原油比重0.8,油水界面的海拔为-950m ,试分析在哪3为-1000m ,位于含水区的一口探井实测地层中部原始地层压力为11.7MPa ,油层中部海拔-1300m ,已知原油比重为0.85,地层水比重为1.0,求该油田油水界面的海拔深度。
4、已知一油藏中的两点,如图,h=10m,P A =9.35MPa, P B =9.5MPa,原油重率γ=0.85,问油的运移方向如何?BA h=10m5、已知一个边长为5cm 正方形截面岩心,长100cm ,倾斜放置,如图所示,入口端(上部)压力1P =0.2MPa ,出口端(下部)压力2P =0.1MPa ,h=50cm ,液体重率0.85,渗流段长度L=100cm ,液体粘度μ=2mPa.s ,岩石渗透率K=12m μ,求流量Q 为多少?6、在上题基础上如果将h 改为0,其结果又将如何?通过计算说明什么?(其它条件不变)7、某实验室测定园柱形岩芯渗透率,岩芯半径为1cm ,长度5cm1mPa.s 建立压差,使粘度为1mPa.s 的液体通过岩芯,在二分钟内测量出通过的液量为15cm 3,从水银压力计上知道两端的压差为157mmHg ,试计算岩芯的渗透率。
渗流力学期末复习提纲

Q 4 Kh
P ws ( t ) Pi
ln t
p
t
t
Horner 公式
第七章 一、概念
油水两相渗流理论
•非活塞式驱油: •水驱油前缘:驱替相水与被驱替相油首次接触面。 •非活塞式驱油含水饱和度变化关系。
S 1.0
1.0-Sor
Swf Swi X0 Xf(t) Xe
•原始溶解气油比:
•生产气油比: 二、油气两相渗流特征 三、油气两相生产特征
第九章
•弹性储容比? •窜流系数?
双重介质渗流理论
•双重介质油藏的基本特征 ?
•双重介质油藏渗流特征 ?
二、计算
Pw f ( t ) Pi Q 4 K h ln 2 .2 5 t rw
2
•不稳定叠加原理
Pi PM
j 1
n
Q j Ei 4 K h
2 rj 4 t j
4 K h ln
j 1
面位置。
P zs PM Z M
B h rAB WOC
判断流动方向
P zsA P zsB
确定界面位置:
P zsA P zsB
A
第二章 概念
油气层渗流基本规律及渗流数学模型
•线性渗流及特征:q~p或v~dp/dr成线性关系或流态呈层流。 •完整的数学模型包括几部分: (1)连续方程: (2)初始条件: (3)内外边界条件:
n
Q j
2 .2 5 t j rj
2
•变产量
PM
( q j q j 1 ) 4 k h
E i[
r
2
高等渗流力学作业-完整版

(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)—— 黄世军 复杂结构井近井渗流理论发展的阶段特征。
高等渗流力学(2017)-绪论+第一章-刚性稳定渗流-程林松

r p(r,t) Q r rRw 2Kh
第一节 渗流数学模型
二、构成: (六)边界和初始条件
封闭边界,定压生产
封闭边界
p(r,t) 0 r rRe
(r, t) rRe 0
定压生产
p(r, t) rRw pw
第一节 渗流数学模型
三、渗流数学模型 1、单相不可压缩液体稳定渗流数学模型
辅助方程
第一节 渗流数学模型
一、定义:
用数学语言综合表达油气渗流过程中全部力学现象和 物理化学现象内在联系和运动规律的方程式(或方程组)称 为“油气渗流的数学模型”。
一个完整的渗流数学模型应包括两部分:渗流综合微 分方程的建立以及边界条件和初始条件的提出。
第一节 渗流数学模型
二、构成:
基
(一)运动方程
安 排
非等温渗流
物理化学渗流
工程问题
12
低渗油藏非线性渗流理论及应用
非牛顿流体的渗流
复杂结构井渗流理论
渗流模型描述 应用
6
专题讨论
3
考试
四、讲课内容及学时分配
工具
压力降叠加
叠加原理 镜像反映(映射)
势叠加 压力叠加 渗流速度叠加
等值渗流阻力
复势叠加
保角变换
程林松
四、讲课内容及学时分配
黄世军
四、讲课内容及学时分配
第一节 渗流数学模型
二、构成: (三)质量守恒方程(连续性方程)
流入质量-流出质量=质量变化量
d iv ( )
( ) t
q
0
流动项
累计项
源汇项
div 称为散度:
div
x
x
y
y
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
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,即有充分的窜流时,渗流过程中的压力变化与单一介 质中的压力变化完全相同。
第三节 双重介质简化渗流模型的无限大地层典型解 当η ≠0时,则公式(7)中的指数函数不是趋于零而是趋于某 一定数,即: a2 t exp exp t 0 2 1 a
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
m m 0 C m pm pi
对于其中的流体(如原油)则:
o 1 C p pi
第二节 双重介质单相渗流的数学模型
C f f f 0 o 1 C f 0
p pi
Kf
Km div(grad p f ) ( pm p f ) 0 pm K m mCm ( pm p f ) 0 t
对(1)式求导带入(2)式并消去压差pm-pf得裂缝系统压力变 化的偏微分方程:
Co div[ grad p f Co grad p f ] 0 t t
Km f Cf div(grad pf ) ( pm pf ) 0 t
忽略
pf
Kf
pm K m Km mCm div(grad pm ) ( pm p f ) 0 t
第三节 双重介质简化渗流模型的无限大地层典型解 上面两式化简为:
构成生产压差大小的主要部分 (即公式(7)中的方括号的值)不可能 等于 1 ,因而双重介质比单一的孔 隙介质中的生产压差要小。 由于窜流能力的不同,也即是 裂缝和孔隙间交换能力不相同时, 井底压力的变化不一样。
不同
值的无量纲井底压力 U (rw , t ) 随时间变化
第三节 双重介质简化渗流模型的无限大地层典型解
C f C f C f 0
Cm m m0 o C t m 0
p p C m0 o m t t
C m Cm C m0
第二节 双重介质单相渗流的数学模型
第一节 双重介质油藏模型
双重介质油藏基本参数:弹性储容比和窜流系数。
f Cf f C f m C m
基质系统孔隙体积 裂缝系统孔隙体积 m f 基质和裂缝系统总体积 基质和裂缝系统总体积
C f Cm
f m ——裂缝和基质系统的孔隙度。
第一节 双重介质油藏模型
流体在双重介质油藏渗流的过程中,基质与裂缝之间存在着 流体交换。窜流系数就是用来描述这种介质间流体交换的物理量, 它反映了基质中流体向裂缝窜流的能力。窜流系数定义为:
t
Co
p f
div[
Kf
grad p f Co
grad p f ] 0 t
在给定初始和边界条件时,方程(3)有解。
第三节 双重介质简化渗流模型的无限大地层典型解 实例:假设有一等厚无限大地层,被一完善井打开,并设井半 径为零,此处有一点源,其产量为Q,则流动为平面径向流, 流动模型如图所示,此时公式(3)就可以展开为:
C m m m 0 o 1 C p pi m 0
由此得到上二式对时间的导数:
C f f f 0 o C t 是零阶虚宗量贝塞尔函数。
对(6)式进行Laplace反演:
r i et U (r , t ) K0 r d r i 2 i (k ) k
积分
裂缝中压力变化公式 :
a2 t Q J 0 ( a, r ) p ( r , t ) pi da 1 exp 2 0 2 K f h a 1 a
p f 1 p f r t t r r r 1 1 p f r r r r
第三节 双重介质简化渗流模型的无限大地层典型解
注意到: t 0 可得新的边界条件为:
p f lim r r 0 r
裂缝网络可以是均匀分布,也可以是非均匀分布的,采 用非均匀的裂缝网格可研究裂缝网络的各向异性或在某一方 向上变化的情况。
第一节 双重介质油藏模型
该模型是把实际的双重介 质油藏简化为由一组平行层理 的裂缝分割基质岩块呈层状的 地质模型,即模型由水平裂缝 和水平基质层相间组成。
对于裂缝均匀分布、基质具有较高的窜流能力和高储存 能力的条件下,其结果与Warren-Root模型的结果相似。
0
p f Q t / lim r 1 e r 0 2 K f h r
第三节 双重介质简化渗流模型的无限大地层典型解
(5)
对式(5)及初边界条件进行Laplace变换得:
1 d dU r r dr dr
U 0 k
双重介质实际油藏模型 基岩 裂缝
基岩
双重介质
裂缝
第一节 双重介质油藏模型
裂缝-孔隙性双重介质结构油藏可抽象地简化成各 种不同地质模型。
1.Warren Root 2.Kazemi 3.De Swaan 4.Factal
模型 模型 模型 模型
第一节 双重介质油藏模型
将双重介质油藏简化为正交裂缝 切割基质岩块呈六面体的地质模 型,裂缝方向与主渗透率方向一 致,并假设裂缝的宽度为常数。
f div v f t q 0
裂缝系统
基岩系统
m div vm q 0 t
对于均质各向同性地层,上式中的对流项可以化简为:
div v f
Kf o div grad p f
高等渗流力学
黄世军 2017
第五章 多重介质渗流理论
第一节 双重介质油藏模型
具有裂缝和孔隙双重储油(气)和流 油(气)的介质我们称之为双重介质。 在一般情况下,裂缝所占的储集空间 大大小于基岩的储集空间,因此裂缝 孔隙度就小于基岩的孔隙度,而裂缝 的流油能力却大大高于基岩的流油能 力,因此裂缝渗透率就高于基岩的渗 透率,这种流油能力和供油能力的错 位的现象是裂缝-孔隙介质的基本特性。
Kazemi也提出计算α的公式:
1 1 1 4 2 2 2 Lx Ly Lz Lx、Ly、Lz——基质岩块在x、y、z方向上的长度,m。
第五章 多重介质渗流理论
第二节 双重介质单相渗流的数学模型
建立双重介质油藏的数学模型时,两种介质分别满足各自 的运动方程、状态方程和连续性方程,而两种连续介质间窜流 通过连续性方程中的一个源和汇函数来表示。
二、Km=0简化模型无限大地层典型解及其应用
p f K f Km div(grad p f ) ( pm p f ) 0 f C f t C pm K m div(grad p ) K m ( p p ) 0 m m m m f t
第三节 双重介质简化渗流模型的无限大地层典型解
Q J 0 ( a, r ) a t p f ( r , t ) pi 1 e da 0 2 K f h a
2
0
r 2 J 0 ( a, r ) 1 a2 t 1 e da Ei a 2 4 t
Km 2 rw Kf
K f Km
m2
窜流系数的大小,既取决于基质和裂缝渗透率的比值,又 取决于基质被裂缝切割的程度,基质与裂缝渗透率的比值越 大或者裂缝密度越大,窜流系数越大。
第一节 双重介质油藏模型
Warren-Root提出的计算α的关系式: