第四章 储层参数计算05
4.储层孔隙结构
(3)最小非饱和孔隙体积百分数Smin
Pc达到仪器最高压力时,水银无法侵入的孔隙 体积百分数 Smin↑→微孔喉体积↑→岩石储集性能越差 Smin=f(颗粒大小、均一程度、胶结类型、 孔隙度、渗透率) Smin=0~90%
Smin讨论
1、Smin取决于所使用仪器的最高压力 •水银系统的毛细管力曲线上:曲线尾部常不平 行于压力轴,仪器的最高压力越高,曲线越向 纵轴偏。此时:
③孔隙缩小型喉道 孔隙与喉道无明显界限,扩大部分为 孔隙,缩小的狭窄部分为喉道。
④管状喉道 特点: 管状喉道,细而长,断面近圆形。 成因: 溶蚀作用形成。负鲕灰岩内鲕粒铸模孔的连
⑤解理缝型喉道 白云石或方解石晶体中被溶蚀扩大的解理缝。
第二节 孔隙结构表征
实验研究方法:二大类 •直接观测法 •岩心观测、铸体薄片法、图像分析法、扫描电镜法等; •间接测定法 •毛细管压力法,主要为压汞法。
①粒间溶孔 颗粒间溶蚀→粒间溶孔。 广义上,粒间溶孔是原生和次生的混合孔隙: •次生粒间溶孔 粒间溶孔中次生溶蚀部分大于原生孔部分。 •混合粒间溶孔 粒间溶孔中原生部分大于次生部分。 ②组分内溶孔 粒内溶孔、杂基内溶孔、胶结物内溶孔、交代 物溶孔等。
③铸模孔 颗粒、生屑或交代物等完全溶解而成。外形 与原组分相同。
③杂基内微孔隙
粘土杂基和碳酸盐泥中存在的微孔隙。 特点: 孔隙极细小,仅在扫描电镜下可清晰辩认。 可形成高孔隙度,但渗透率很低。
④层面缝
具剥离线理的平行层理纹层面间的孔缝。 在一系列厘米级甚至毫米级厚度的平板薄层 间,为力学性质薄弱的界面,极易剥离,其 界面间即为层间缝。
(2)次生孔隙
次生作用→
一、碎屑岩储层孔隙和喉道类型
1、孔隙类型 (1)分类 •成因分类: 原生、次生及混合成因孔隙。目前国内外比较流行的分类方案 •孔隙大小: 超毛细管孔隙、毛细管孔隙、微毛细管孔隙 •孔隙成因和几何形状: E.D.Pittman,1979 粒间孔隙、溶蚀孔隙、微孔隙、裂缝孔隙 •综合性分类: 按成因:分为原生和次生孔隙二大类;然后,按孔隙产状和几何形状再 进一步细分。
油气层产能保护第4章储层损害的评价方法
第四章储层损害的评价方法代表性岩样的选取(1)速敏概念和实验目的(1)速敏概念和实验目的(2)原理及作法(2)原理及作法C、注意事项a、在实验过程中必须保持连续流动。
如果中途停止流动,会使运动着的微粒在孔道处沉积,破坏微粒分布状态,即使间断后再流动也不能恢复到停止前的状态,此时表现出压力波动很大,实验资料发生矛盾或混乱的现象。
b、对于采油井,要用煤油作实验流体,并要将煤油先经过干燥,再用白土除取其中的极性物质,然后用G5砂心漏斗过滤。
对于注水井,应使用经过过滤处理的地层水(或模拟地层水、标准盐水)作为实验流体。
(4)影响速敏性的因素A、主要受岩石本身性质的影响B、流体矿化度、离子组分、pH值等流体性质的影响随注入流体矿化度的降低而降低,或者随pH值的升高而降低。
Vc如果储层具有较强的速敏性损害,应在工程中选用粘土稳定剂,控制注入或产出流体速度等预防措施。
(1)水敏概念和实验目的(1)水敏概念和实验目的(2)原理及评价指标(2)原理及评价指标(1)盐敏概念和实验目的(1)盐敏概念和实验目的(2)原理及评价指标(2)原理及评价指标(1)碱敏性的概念和实验目的(1)碱敏性的概念和实验目的(2)原理及评价指标(2)原理及评价指标(1)酸敏性的概念和实验目的(1)酸敏性的概念和实验目的(2)原理及评价指标(2)原理及评价指标1、实验条件(1)该实验可用气体、中性煤油或标准盐水(质量分数8%)作为实验流体。
(2)使用特制的可分别控制或测量轴向和径向应力的驱替装置。
(3)用气体做实验流体时,按SY/T6385执行。
2、净围压的应力敏感性评价实验(1)实验步骤a、损害前液体渗透率的测定。
b、保持进口压力值不变,缓慢增加围压,使净围压依次为2.5MPa,3.5MPa,5.0MPa,7.0MPa,9.0MPa,11MPa,15MPa,20MPa。
c、每一压力点持续30min后,测定岩样渗透率。
d、缓慢减小围压,使净围压依次为15MPa,11MPa,9.0MPa,7.0MPa,5.0MPa,3.5MPa,2.5MPa。
第5章-存储层次教学文案
2. 回写法(write back)
只写入Cache,在被替换时才写回主存
系统结构
直写和回写的比较
存储层次
❖ 直写策略
优点:简单可靠 缺点:总线操作频繁、影响工作速度 解决方法:在Cache与主存间设置一级/多级缓
冲器,形成实用的“缓冲直写”方式,提高速 度
CPU时间= (CPU执行周期数+存储器停顿周期数)×时钟周期时间
❖ 只考虑Cache失效引起停顿,并将Cache命中 所用时钟周期数归为CPU执行时钟周期数的部 分
存储器停顿时钟周期数=访存次数×失效率×失效开销
系统结构
CPU执行时间
存储层次
CPU时间= IC(CPEIXECUTI+ ON存储指器令停条顿数周) 期 时数钟周期时间
❖ 统一Cache结构的失效率低于分离Cache结构
系统结构
分离和统一Cache性能比较:例题5.1-2 存储层次
❖ 平均访存时间
=取指比例×(指令命中时间+指令失效率×失效开销) +数据比例×(数据命中时间+数据失效率×失效开销)
❖ 分离Cache的平均访存时间
=75%×(1+0.64%×50)+25%×(1+6.47%×50) =2.05 (已知:命中时间=1,失效开销=50)
结束
向Cache传送含RA的主存块
系统结构
5.1.4 存储层次的四个问题
存储层次
1. 映像规则(block placement)
Where can a block be placed in the upper level?
2. 查找方法(block identification)
How is a block found if it is in the upper level?
储层参数计算范文
储层参数计算范文储层参数计算是石油地质与储层工程中的一项重要工作,它是评价储层性质及开发潜力的基础。
储层参数包括孔隙度、渗透率、孔隙度分布、渗透率分布等。
在储层参数计算中,需要利用地质资料、地震资料、测井资料等多种信息进行综合分析和计算。
首先,孔隙度是储层中岩石或矿物颗粒间的孔隙体积与岩石或矿物颗粒体积的比值。
孔隙度可以通过测井资料中的密度和孔隙度曲线来直接计算。
具体的计算公式为:孔隙度(PHI)=((滴定孔隙度-最小孔隙度)/(最大孔隙度-最小孔隙度))*100%其中,最小孔隙度是指岩石或矿物颗粒之间最小的孔隙度,最大孔隙度是指岩石或矿物颗粒之间最大的孔隙度。
其次,渗透率是岩石中流体渗透的能力,一般以Darcy为单位。
渗透率的计算可以通过测井资料中的电阻率、声波速度等参数来间接计算。
其中,最常用的计算方法是根据Kozeny-Carman公式计算。
具体的计算公式为:渗透率(K)=(φ3/((1-φ)^2*(1-φ)^3))*((D^2)/(180μ))其中,φ为孔隙度,D为岩石颗粒的平均直径,μ为岩石渗流介质的动力粘度。
此外,在储层参数计算中,孔隙度分布和渗透率分布的计算也是非常重要的。
孔隙度分布主要针对储层中不同岩石单元或层段的孔隙度进行计算和分析。
渗透率分布主要针对储层中不同位置或不同岩性的渗透率进行计算和分析。
这些分布的计算方法可以基于地质资料、地震资料、测井资料、实验数据等进行综合分析和计算。
综上所述,储层参数计算是石油地质与储层工程中的一项重要工作。
通过合理的分析和计算,可以评价储层的性质和开发潜力,并为储层的合理开发提供依据。
在实际应用中,还需要结合其他地质和工程参数进行综合分析,以取得更加准确和可靠的结果。
储层参数间的关系
• 教学目的: 教学目的: • 掌握毛管渗流定律、等值渗流阻力原 理、渗透率和孔隙半径的关系、渗透 率与岩石比面的关系。 • 教学重点和难点:毛管渗流定率、等 教学重点和难点: 值渗流阻力原理。 • 教法说明:课堂讲授 教法说明: • 教学内容: 教学内容:
第六节
储层参数间的关系
(
)
1 2 K = 2 φ⋅ r 8τ 2 1 4φ = 2 φ⋅ 2 8τ S 3 1 φ = 2 2 2τ S
Q S = φS φ = (1 − φ )Sr
φ φ φ ∴K = 2 2 = 2 2 = 2 2 2 2τ ⋅ S 2τ ⋅ Sφ 2τ (1−φ) Sr
3 3
k 则有: 令: = 2τ 则有:
QS = Sφ ⋅ φ = (1 − φ)Sr φ φ φ ∴K = 2 = 2 = 2 2 2S 2Sφ 2(1 − φ) Sr
3 3
考虑迂曲度
nA ⋅ πr ⋅ (τL ) 2 φ= = nπr τ AL
2
τ 则有:
nA⋅ 2πr ⋅ (τL) 2φ 2 2 S= = 2πnrτ = nπτr = AL r r
2
φ φ φ ∴K = 2 = 2 = 2 2 kS kSφ k(1 − φ) Sr
3 3
高才尼高才尼-卡尔曼方程
4
即:
1 2 K = 2 φ⋅ r 8τ 8K r= ⋅τ φ 2τ K r= 7000 φ K : cm
2
K : Dc
四、渗透率与岩石比面的关系
对于假想岩石: 对于假想岩石:
nA(2πrL) 2φ S= = 2πrn = AL r
2 3
1 2 1 4φ φ ∴ K = φ⋅ r = ⋅ φ⋅ 2 = 2 8 8 S 2⋅ S
储层参数定量计算
进一步认识后的形象化描述。 即:
理论分析 实验研究
测井信息 地质信息
资料统计
测井信息 = F(地质信息) + 误差
(测井解释模型)
用数学式表达:y=F(x)+ε
仪器模型: 测井数据集、测井 值与记录值转换关 系、测量误差、环 境校正等
地质模型: 矿物成分、 流体类型等
数学模型: 数学关系式
评价地层
• 一、阿尔奇(Archie)公式
R01 R02 ... R0n F F R0
Rw1 Rw2
Rwn
Rw
• —称为地层电阻率因素
• —全含水地层电阻率
• —地层水电阻率
• Archie还通过实验发现:
• F与孔隙度及岩石孔隙结构有关(反映岩石本身
物理性质),与孔隙中是否含油气及Rw无关。
而Swt Swb Swf
因此Ct t2Swt SwbCwb Swt Swb Cwf
即Rt
t2
Swb Rwf
Rwf Rwb Swt Swb
Rwb
S wt
5.4 渗透率计算
• 一、影响因素
• 孔隙度:孔隙度↑,渗透空间大,渗透率↑ • 泥质含量:Vsh ↑,阻塞渗透通道,K↓ • 砂岩颗粒大小(用粒度中值度量):越细,k越小。 • 裂缝:裂缝愈发育,K ↑ ↑ • 压力对渗透率的影响
§5 储层参数定量计算
5.1 测井解释模型
一. 模型的基本含义
所谓模型:就是客观事物被认识后,经过抽象 ,再组合为易于理解的形象,即形象化的抽象。 模型的建立过程是:
实际 经 过 抽象
变成
简 化 形象
二. 测井解释模型
指测井信息与地质信息间的宏观关系。是经过理
储层参数计算.
SP SP min I sh SP max SP min 2GCUR I sh 1 Vsh GCUR 2 1
– 方法二:
SP Vsh 1 SSP
• SSP:目的层段中纯砂岩水层静自然电位。
三、储层参数计算及处理
4、泥质含量计算
3)自然伽马能谱
– 同时测量U、TH、K、SGR(无UGR)、CGR,除U外,都能用 来求泥质含量,用CGR较好。 – 方法同GR. – 当用多种方法求泥质含量时,最终应取多种方法求得的最 小者
¬ md Ê £ ¸ Â ø Í É
×Ï ¿ ¶ ¶ È £ ¬ £ ¥
三、储层参数计算及处理
2、渗透率计算
2)用孔隙度和束缚水饱和度求取 如Timur公式
0.136 4.4 K 2 S wb
•Swb,%;φ,%;K,10-3μm2
三、储层参数计算及处 b D ma f
H H ma N H f H ma
S
t t ma t f t ma
ma——骨架 f ——流体
该公式称为平均时间公式或Wyllie-Rose公式
三、储层参数计算及处理
1、孔隙度计算
2)体积法:
适用范围:平均时间公式适用于压实和胶结良好的纯砂岩地层。在这
种砂岩中,矿物颗粒间接触良好,孔隙直径较小,故可以忽略矿物颗粒与孔隙流
体交界面对声波传播的影响,可认为声波在岩石中是直线传播的。但是对于未胶
结、又未压实的疏松砂层,矿物颗粒间接触不好,故矿物颗粒与孔隙水的交界面 对声波传播影响较大,使孔隙度相同的疏松砂层的声波时差要比压实的砂层大, 因此需要用压实校正系数校正: S
1 t t ma Cp t f t ma
中国石油大学《储层表征与建模》吴胜和 第四章
变差函数及结构分析
克里金插值方法
变差函数及结构分析
地质变量相关性的各向异性
☼1 ☼3 ☼3
☼3 ☼1 ☼1
☼3 ☼2 ☼1
h a
γ (h) = C(0) – C(h)
基台值(Sill):代表变量在空间上的总变异性大小。
即为变差函数在h大于变程时的值, 为块金值c0和拱高cc之和。 (拱高:在取得的有效数据的尺度上,可观测得到的 变异性幅度大小)。
克里金插值方法 平稳假设
二阶平稳
任何统计学均 要求平稳性假设
• 空间各点处随机变量的集合构成一个随机函数。
(将空间位置作为随机函数的自变量)
P
① 在研究区内有Z(u)的数学期望存在, 且等于常数,即:E[Z(u)] = E[Z(u+h)] = m(常数) ② 在研究区内,Z(u)的协方差函数 Cov{Z(u),Z(u+h)}存在且平稳 (即只依赖于滞后距h,而与u无关)
三维空间的水道迁移
A7 A2 A6 A1 A4 A3 A5
单井侧积层界面点提取
侧积面三维视图
确定性建模概述
地质模式拟合
确定性建模概述
三、数学插值
1. 传统数学插值
侧积层镂空视图
传统数学插值 克里金插值
如:三角剖分法(三角网方法)、 距离反比加权法等 将变量视为纯随机变量, 未考虑变量的空间结构性 仅考虑待估点位置与已知数 据位置的相互关系。
i =1
最小的估计方差,即克里金方差,可用以下公式求解:
当随机函数不满足二阶平稳,而满足内蕴(本征)假设时, 可用变差函数来表示克里金方程组: ⎧n γ (xi − x j )λi + μ = γ (x0 − x j ) ( j = 1,K, n ) ⎪ ⎪∑ i =1 ⎨ n ⎪ λi = 1 ∑ Z*(x0) ⎪ i =1 ⎩
天然气工程教程第4章气藏物质平衡、储量计算及采收率
(1 Sw So )(1 yw ) (1 Swi )(1 ywi )
1C f
( pi
p)
p Z
pi Zi
pi Zi
Gp G
(1 Sw So )(1 yw) (1 Swi )(1 ywi )
1C f
( pi
p)
p Z
0
Gp
G
说明:
在应用上述物质平衡方程时,需要知道两相 偏差系数与凝析油的饱和度,这些需要通过凝析 气井的取样和实验室分析进行测定。
假定原始条件下,地层压力大于露点压力, 则有原始地下储集空间为 :
VPi
GBgi (1 S wi )(1
yW i )
原始条件水 的体积分数
(1) 地层压力大于露点压力
目前的孔隙空间 为气和水所占 :
VP
(G GP )Bg (1 SW )(1 yW )
由于压力下降,气层 岩石的形变体积:
Gp G
P/Z
0
岩石和流 体压缩性 同时作用
只有流 体压缩
G
Gp
求储量的另一 “归一”化处理:
p Z
(1 Cep)
pi Zi
pi Zi
Gp G
纵轴上截距: a pi Zi
斜率: b pi 1 Zi G
外推直线至:
p 0 与横轴交点
Z
即为G。
pi
p Z
(1
Ce
p
)
Zi
0
Gp G
五、气藏物质平衡方法应用中的注意事项
凝析油采收率:
EcR 2.09 107 ( pi )0.9027(Ri )0.25084( o )2.25253 (141.5 131.5 o )2.50337 (1.8T 32)0.30084
储层流量计算公式
储层流量计算公式
储层流量计算公式流量公式如下:流量=3600x管子面积x流速,管子面积=3.14x(半径)2,0.785x(半径)2管子面积=0.785x(半径)2 =0.785x(半径3/时,流速单位为米/半径或直径单位为米:式中流量单位为米3/时流速单位为米/秒:半径或直径单位为米3600小时折合成的秒数因为面积(2)x流速流速(3/秒是1 小时折合成的秒数因为面积米2)x流速(米/秒)的结果得米3/秒指秒钟内的流量因此,3/时这个数1秒钟内的流量,因此折合成1小时的流量米3/时,就要乘上3600 这个数,的管子,秒时流量是多少?例1 Dg100 的管子流速为
1米/尼时流量是多少7流量-3600X管子面积X流速先把查保换镇成管子面米2
流量=3600x0.00785x1=28.33/时则。
并行存储器系统
取
层2:M2
主存储器
成
时
本
间
层3:M3
磁盘存储器
增
增
加
加
层4:M4
磁带机
五个参数:
存取时间ti:从CPU到第i层存储器的往返时间 存储器容量Si:第i层的字节的数量 每字节成本Ci:第i层存储器的成本为CiSi 传输带宽bi:相邻层之间传送信息的速率 传输单位Xi:i和i+1层之间数据传送的粒度
对存储器系统中各层次存储器的特性,某统 计数据如下表:
M1:高速缓存
CPU寄存器 字单位
……
a
……
块单位
M2:主存储器
a 页面A
页单位
M3:磁盘 存储器
M4:磁带机 后援存储器
a
b
页面A 页面B
段F
段单位
a
b
页面A 页面B
段F
b
a,b为高速缓存
块,32个字节
b 页面B
段G
段G
4.2.2 一致性(coherence)
1.一致性定义
同一个信息项与后继存储器层次的副本是 一致的。
4.2.3 局部性(locality)
Hennessy和Patterson(1990年)提出了一 条90-10规则:典型程序在10%的代码上可能 要耗费其执行时间的90%(例如嵌套循环操作 的最内层循环)。
时间局部性(temporal locality):最近 的访问项(指令或数据)很可能在不久的将来 再次被访问。即对最近使用区域的集中访问。
空间局部性(spatial locality):一个进程 访问的各项的地址彼此很近,例如,表操作或 数组操作含对地址空间中某一区域的集中访问。
储层参数计算
即S w
S wt S wb 1 S wb
三、储层参数计算及处理
3、饱和度计算 长庆油田 大庆油田 印度尼西亚公式 阿尔奇公式 西南油气田 阿尔奇公式和孔隙度—饱和度交会
斯伦贝谢公司 核磁共振测井计算饱和度
三、储层参数计算及处理
4、泥质含量计算
– 泥质含量是泥质体积与岩石体积之比。 • 1)自然伽马
。
三、储层参数计算及处理
5、Gxplorer储层参数计算现状
1)成熟方法——单孔隙度分析程序POR
是一种孔隙度测井资料加上其它有关资料对泥质砂岩进行分析解 释。可采用自然伽马(GR)、补偿中子(CNL)、自然电位(SP)、中 子寿命(NLL)和电阻率(RT)等五种方法计算地层的泥质含量SH 相 对体积;利用密度测井(DEN)、声波测井(AC)或补偿中子(CNL) 三种孔隙度测井之一计算地层的孔隙度,并且进行泥质校正;计算出 可动油气参数、流体性质分析参数、渗透率和出砂指数等。
ma b D ma f
H H ma N H f H ma
S
t t ma t f t ma
ma——骨架 f ——流体
该公式称为平均时间公式或Wyllie-Rose公式
三、储层参数计算及处理
1、孔隙度计算
2)体积法:
适用范围:平均时间公式适用于压实和胶结良好的纯砂岩地层。在这
K=f(Md,φ)
Md一般可由GR或SP求
4)用裂缝宽度估计
K=8.3×106w2
K:10-3μm2 W为裂缝宽度,cm
三、储层参数计算及处理
2、渗透率计算
5)由核磁共振测井求取 • SDR模型
K C1 (1为常数;φnmr为核磁共振孔隙度,%;T2g为T2几何平均值, ms. • Coates 模型
中国石油大学(北京)《储层表征与建模》第四章 确定性建模
地震属性与地质参数的关系
地震属性的确定性转换
确定性转换
地震属性
地质参数
速度、 波阻抗、 振幅、
频率、
波形结构类型
地质相、 岩性、 孔隙度 流体(?)
二、储层沉积学方法
通过井间砂体确定性对比, 建立储层结构模型
构型分析
构型建模
模式拟合
实际资料 定量模式
模式拟合
实际资料 定量模式
=
1 2
Var[Z(x)-Z(x+h)]
=
1 2
E[Z(x)-Z(x+h)]2-{E[Z(x)-Z(x+h)]}2
(h)
=
1 2
E[Z(x)-Z(x+h)]2
地质统计学中最常用 的基本公式之一。
变程(Range) :指区域化变量在空间上具有相关性的 范围。在变程范围之内,数据具有相关性;而在变 程之外,数据之间互不相关,即在变程以外的观测 值不对估计结果产生影响。
第四章
确定性建模
Deterministic Modeling
确定性建模概述 地质统计学克里金方法
三维地质建模
数据库
油藏数模
模型粗化
三维构造建模 三维相建模
三维储层参数建模
资料:井数据 和/或地震数
据
储层建模
三维网格赋值
(相控)
属性:相 孔、渗、饱
确定性建模 随机建模
离散型地质变量 连续型地质变量
n i 1
C
xi
xj
i
C x0
n
xj
i 1
i 1
储层地球物理学2-储层参数预测
以保证Vi的计算精度 应取足够多的井,至少在10井以上
实际应用见程序计算和量版的制作
b 散点法(精度较差) 在井少或无井区,
可用地震层速度作岩 性指数量板: I 由速度谱 Vn 、Z 。 II 在Z—V平面上点 出所有散点。 III 高速边界为100% 砂岩,低速边界为 100%泥岩。 IV 内插其它曲线。
地震储层预测主要利用地震波的动力学特征(如振 幅、速度、相位、频率等)来确定储集层的分布范围。
•储层地震预测方法
•地 震 反 演 • 多属性综合分析方法 •模 式 识 别 预 测 法 •地 震 相 分 析 法 •相 干 分 析 法 •多 尺 度 边 缘 检 测
1)地震反演
地震反演是利用地表地震观测资料,以已知地质规 律和钻井测井资料为约束,对地下岩层物理结构和物理 性质进行成像(求解)的过程。
第二步对V=a Zb两边取对数 则有:lnV = lna +blnZ 令lna = c, lnV=y, lnZ=x 得到: y = bx + c 用最小二乘法可求出a、b、c(a=ec)
具体计算公式如下:
Vi、Zi为实测数据,分别为层速度和中 心点埋深,N为数据对的数目。
第三步根据实测数据对Vi、Zi计算出纯砂、泥岩和 其它不同含量的V-H曲线,即得岩性量板。
(2) 稳定性原则 成组出现的能量团反映一个或一组反射波,一 般仅拾取其中的一个,该能量团应反映速度的 稳定性。
(3) 连续平滑的曲线 解释的迭加速度曲线应连续平滑,没有突变拐 点。
(4) 识别异常波 多次波特征明显,较易识别,但与断层有关的 断面波特征不明显,变化复杂,要结合地震剖 面对比解释。
储集层的物性参数~(doc)-下载
第一节储集层的物性参数储集层的基本特征是具孔隙性和渗透性,其孔隙渗透性的好坏、分布规律是控制地下油气分布状况、油气储量及产量的主要因素。
一、储集层的孔隙性绝对孔隙度:岩样中所有孔隙空间体积之和与该岩样总体积的比值。
是衡量岩石孔隙的发育程度。
Pt=V p/V t*100%按岩石孔隙大小,有超毛细管孔隙、毛细管孔隙和微毛细管孔隙三类。
1.超毛细管孔隙:直径>0.5mm,相应裂缝宽度>0.25mm,液体在重力作用下自由流动。
2.毛细管孔隙:直径0.5~0.0002mm,裂缝宽度0.25~0.0001mm,由于毛细管力的作用,液体不能自由流动。
3. 微毛细管孔隙:直径<0.0002mm,裂缝宽度<0.0001mm,液体在非常高的剩余流体压力梯度下流动。
有效孔隙度:指彼此连通的,且在一般压力条件下,可以允许液体在其中流动的超毛细管孔隙和毛细管孔隙体积之和与岩石总体积的比值。
Pe=V e/V t*100%二、渗透性渗透性:指在一定的压差下,岩石允许流体通过其连通孔隙的性质。
对于储集层而言,指在地层压力条件下,流体的流动能力。
其大小遵循达西定律。
K即为岩石的渗透率,国际单位为μm2,常用单位为达西(D)。
国际单位:μ=1Pa.s △P=1Pa F=1m2 L=1m Q=1cm3/s则:K=1μm2常用单位:μ=1厘泊△P=1大气压F=1cm2 L=1cm Q=1cm3/s则:K=1D=1000md1D=0.987μm21D=987*10-6μm2绝对渗透率:单相液体充满岩石孔隙,液体不与岩石发生任何物理化学反应,测得的渗透率称为绝对渗透率。
有效渗透率:储集层中有多相流体共存时,岩石对每一单相流体的渗透率称该相流体的有效渗透率。
油气水分别用Ko、Kg、Kw表示。
相对渗透率:对每一相流体局部饱和时的有效渗透率与全部饱和时的绝对渗透率之比值,称为该相流体的相对渗透率。
油气水分别表示为Ko/K、Kg/K、Kw/K。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第三节 含水饱和度计算
• 一、阿尔奇(Archie)公式
R01 R02 ... R0n F F R0
Rw1 Rw2
Rwn
Rw
– F —称为地层电阻率因素
– R0 —全含水地层电阻率 – Rw—地层水电阻率
• Archie还通过实验得发现:F与孔隙度及岩石 孔隙结构有关(反映岩石本身物理性质), 与孔隙中是否含油气及Rw无关。即:
2
aRw
Sw2
– α=1--2
• “印度尼西亚公式”
1 Vccl Rt Rcl
2
aRw
S
2 w
– C=1-Vcl/2
• 二、Waxman-Smits模型
– 1. Waxman- Smits模型
• Waxman- Smits认为:在粘土表面均具有 负电荷,这样在粘土表面吸附着部分阳离子, 这些阳离子又能很容易地与水溶液中的其它 阳离子发生交换,粘土表面及其附近的的高 浓度阳离子是造成所谓粘土表面导电性的主 要原因。
2.全含水纯地层孔隙度方程
所谓纯地层,即不含泥质的地层,即
Vsh 0 ,若不含油气,则以上通式可写
成:
X X
ma
X X
ma
w
由于孔隙度系列测井探测 深度浅,在探测范围内孔 隙中流体为泥浆滤液
X X
ma
X X
ma
f
1)密度测井:将X写成ρ,则
D
ma ma
b f
0.01
0.1
1
孔隙度,小数
50MPa-90℃条件下电阻率增大指数与含水饱和度 关系(油驱水)
100
电阻率增大系数
10
y = 1.0175x -1.8896 R2 = 0.9709
1
0.1
1
含水饱和度,小数
• 原始地层含水饱和度
Sw
n
abRw
mRt
• 在冲洗带中,Rt → Rxo,Rw → Rmf,Sw → Sxo,得:
上式同除以t ,则Ct
t
2
S
2 wt
SwbCwb Swb
Swf Cwf Swf
而Swt Swb Swf
因此Ct t2Swt SwbCwb Swt Swb Cwf
即Rt
t2
Swb Rwf
Rwf Rwb Swt Swb
Rwb
– 这种导电性与可发生阳离子交换的数量 相关,该量称阳离子交换能力CEC,它表 示每100克干岩石可交换的阳离子数。单 位为mmol/cm3。
– 阳离子交换能力与地层的粘土矿物成分 有关。如:绿泥石为10-40,伊利石为 10-40,高岭石为3-15,蒙脱石为80150。
– 还有一种表示阳离子交换能力的方法, 称阳离子交换容量 Qv ,它表示岩样每 单位总孔隙体积中含有的可交换阳离子 的摩尔数,单位为mmol/cm3。
5.排列方式:立方体排列, φ=47.6%,最 大,排列最宽敞,最不“稳 定”; 菱面体排列,φ=25.96%,最 小,排列最紧密,最“稳定”
6.胶结作用:总的来说,胶结作用愈强, φ越小。另一方面,钙质胶结比泥质 胶结更致密,即φ更小;硅质胶结最 致密,对孔隙度影响最大。即泥质好 于钙质好于硅质。
• Archie还通过实验发现:地层电阻率与 R0有关,同时随So的增大而增大。定义:
I Rt R0
– “电阻率增大系数”,“地层电阻率指数” 。
• 通过实验得到:
I
b
S
n w
– n为与孔隙导电流体分布有关的量,称“饱和 度指数”。
– 取值范围1.3—2(水湿);2—10(油湿)
– b称“Archie常数b” ,常取1 。
• 因此:
F
R0 Rw
a
m
I
Rt R0
b
S
n w
Rt
abRw
m
S
n w
– 由岩电实验作F~Φ关系图求a、m;作I~ Sw关系图求b、n 。
50Mpa-90℃条件下地层因素与孔隙度关系 10000
1000
地层因素F
100
10
y = 0.8051x-1.9029
R2 = 0.9793
1
2)中子孔隙度:
N
H Hma Hf Hma
3)声波测井:
S
t t ma tf t ma
该公式称为平均时间公式或Wyllie-Rose公式
注意:平均时间公式适用于压实和胶结
良好的纯砂岩地层。在这种砂岩中,矿物颗 粒间接触良好,孔隙直径较小,故可以忽略 矿物颗粒与孔隙流体交界面对声波传播的影 响,可认为声波在岩石中是直线传播的。但 是对于未胶结、又未压实的疏松砂层,矿物 颗粒间接触不好,故矿物颗粒与孔隙水的交 界面对声波传播影响较大,使孔隙度相同的 疏松砂层的声波时差要比压实的砂层大,因 此需要用压实校正系数校正:
如:纯砂岩水层=砂岩+孔隙(只有水)
水
砂岩 骨架
V=1个单位
1.通式。
考虑一般情况:地层含泥质、孔隙中含油气,那么体 积模型可画为:
含油气 泥质砂 岩体积 模型
砂岩骨架
Vma
泥质 油气
VSVshhsn
水 Sw
V=1
设我们研究的物理参数为X,用
X ,X X X , ,
h
w
ma
sh
分别表
2Ish y 2 y
当I sh
0.5时:Swb
2I sh 2 y
f
f
Sw
f f h
( f
t
h ) t
(t
t b )
t
Swf 1 Swb
即Sw
Swt Swb 1 Swb
第四节 渗透率计算
• 一、影响因素
– 孔隙度:孔隙度↑,渗透空间大,渗透率↑ – 泥质含量:Vsh ↑,阻塞渗透通道,K↓ – 砂岩颗粒大小(用粒度中值度量):越细,k
F
*
a
tm
• 三、双水模型
– 1. 双水模型的实质
• 泥质砂岩中有两种水:束缚水,即束缚在粘 土表面的水,又称近水;自由水,即位于连 通孔隙空间的水,或称远水。也就是说,泥 质的存在使地层水的导电性比纯地层下的导 电性好,并认为这种增加的导电性是粘土颗 粒表面的束缚水所产生的阳离子交换所引起 的,泥质地层的全部导电能力是由自由水和 束缚水并联导电的结果。
7.白云化作用:白云化作用越强, φ越 大。主要是因为产生白云化作用后, 骨架体积变小;灰岩被白云岩置换 后体积缩小12—13%
二.孔隙度测井响应方程
岩石体积物理模型:根据岩石的组成,按其物 理性质的差异,把单位体积岩石分成相应的几个 部分,各部分可以看作是性质相同的小单元,然 后研究每一部分对岩石宏观物理性质的贡献之和。
• 如Timur公式
0.136 4.4
K Sw2b
• Swb,%;φ,%;K,10-3μm2
–3.由粒度中值与孔隙度求K
•K=f(Md,φ) •Md一般可由GR或SP求
示该物理参数对油气、水、骨架、泥质的测井响应值,
由体积物理模型的思想,测井响应值X应为各部分贡
献之和。即有:
b a
X S X S X V X V
X
hh
ww
ma ma
sh sh
Vma Vsh 1
Sh Sw 1
X X X X X V X X
F a
m
– m取决于孔道的弯曲程度,颗粒的形状和排列方 向,以及胶结情况。
– “胶结指数”,“结构指数”,“孔隙度指数”, “Archie指数”。
– 取值范围为1.3-3,随胶结程度的变好,该值增大, 常取2。
• a为与岩性有关的比例常数,变化范围 为0.6~1.5,常取1。
• 简称为“Archie常数a”
Sxo
n
abRmf
m Rxo
• 阿尔奇公式适用于纯地层,对泥质地 层,计算饱和度的公式很多,如:
• Simandoux公式
1 Rt
Vcdl Rcl
n
Sw2
aRw
m
(1Vcdl
)
Swn
– 常取m=n=2;d=1—2,常取d=1
• “尼日利亚”公式
1
Vcl
Rt Rcl
ma
ma
X
w
sh
ma
ma
sh
w
S XX XX 1 h
ma
ma
h w
1
Sh 0
不含油气时
X
X X X
ma
ma
sh
V
sh
X X X X ma
w
ma
w
T Vshsh
(泥质校正公式)
Qv
CEC(1t )G t
t 泥质砂岩的总孔隙度,小数;
G 岩石的平均颗粒密度, g / cm3
CEC 岩石的阳离子交换能力, mmol / g
Qv 岩石的阳离子交换容量, mmol / cm3
–2. Waxman- Smits饱和度方 程
Cwe Cw BQv / Swt
S
1 Cp
t t ma tf t ma