测井解释计算常用公式目录1. 地层泥质含量(Vsh)计算公式 (1)2 . 地层孔隙度(φ)计算公式 (4)3. 地层含水饱和度(Sw)计算 (7)4. 钻井液电阻率的计算公式 (12)5. 地层水电阻率计算方法 (13)6.确定a、b、m、n参数 (21)7.确定烃参数 (25)8. 声波测井孔隙度压实校正系数Cp的确定方法 (26)9. 束缚水饱和度(Swb)计算 (26)10. 粒度中值(Md)的计算方法 (29)11. 渗透率的计算方法 (29)12. 相对渗透率计算方法 (35)13. 产水率(Fw) (36)14. 驱油效率(DOF) (37)15. 计算每米产油指数(PI) (37)16. 中子寿命测井的计算公式 (37)17. 碳氧比(C/O)测井计算公式 (39)18. 油层物理计算公式 (46)19. 地层水的苏林分类法 (49)20.毛管压力曲线的换算 (50)21. 地层压力 (51)附录:石油行业单位换算 (53)测井解释计算常用公式1. 地层泥质含量(Vsh )计算公式1.1 利用自然伽马(GR )测井资料1.1.1 常用公式m in m ax m inGR GR GR GR SH --= (1)式中,SH -自然伽马相对值;GR -目的层自然伽马测井值;GRmin -纯岩性地层的自然伽马测井值;GRmax -纯泥岩地层的自然伽马测井值。
1212--=⋅GCUR SH GCUR sh V (2)式中,Vsh -泥质含量,小数;GCUR -与地层年代有关的经验系数,新地层取3.7,老地层取2。
1.1.2 自然伽马进行地层密度和泥质密度校正的公式o sh ob sh B GR B GR V -⋅-⋅=max ρρ (3)式中,ρb 、ρsh -分别为储层密度值、泥质密度值;Bo -纯地层自然伽马本底数;GR -目的层自然伽马测井值;GRmax -纯泥岩的自然伽马值。
1.1.3 对自然伽马考虑了泥质的粉砂成分的统计方法C SI SI B AGR V b sh +-⋅-⋅=1ρ (4)式中,SI -泥质的粉砂指数;SI =(ΦNclay -ΦNsh )/ΦNclay (5)(ΦNclay 、ΦNsh 分别为ΦN -ΦD 交会图上粘土点、泥岩点的中子孔隙度) A 、B 、C -经验系数。
1.2 利用自然电位(SP )测井资料α-=--=0.1minmax min SP SP SP SP sh V ..…….……………(6) 式中,SP -目的层自然电位测井值,mV ;SPmin -纯地层自然电位值,mV ;SPmax -泥岩层自然电位值,mV 。
α-自然电位减小系数,α=PSP/SSP 。
PSP 为目的层自然电位异常幅度,SSP为目的层段纯岩性地层的自然电位异常幅度(静自然电位)。
1.3 利用电阻率测井资料b sh R R t R t R R sh R sh V /1])lim ()lim ([-⋅-⋅= ………………………..………(7) 式中,Rlim -目的层井段纯地层最大电阻率值,Ω·m ;Rsh -泥岩电阻率,Ω·m ;Rt -目的层电阻率,Ω·m ;b -系数,b =1.0~2.01.4 中子-声波时差交会计算B A sh V /= (8)f T Nma ma T Nma t f T ma T N A ⨯Φ+--Φ∆--Φ=)1()( ))(1()1)((f T sh T Nma Nsh f T ma T B --Φ--Φ-= 式中,Tma 、Tf -分别为岩石骨架声波时差、地层流体声波时差;ΦNma 、ΦNsh -分别为岩石骨架中子值、泥岩中子值,小数;Δt -目的层声波时差测井值;ΦN -目的层中子测井值,小数。
1.5 中子-密度交会计算B A sh V /= (9)ma Nma f f ma N Nma b A ρρρρρ+Φ⨯--Φ--Φ=)()1())(1()1)((f ma Nsh Nma f sh B ρρρρ--Φ--Φ-=式中,ρma 、ρf -分别为岩石骨架密度值、地层流体密度值,g/cm 3;ΦNma 、Φsh -分别为岩石骨架中子值、泥岩中子值,小数;ρsh -泥岩密度值,g/cm 3;ρb 、ΦN -目的层密度测井值,g/cm 3、中子测井值,小数。
1.6 密度-声波交会计算B A sh V /= (10)f t ma ma t f f ma t f t ma t b A ∆⨯+∆⨯--∆-∆-∆=ρρρρρ)()())(())((f ma f t sh t f t ma t f sh B ρρρρ-∆-∆-∆-∆-=1.7 利用自然伽马能谱测井1.7.1 钍曲线(TH )如果有自然伽马能谱测井,则优先选用能谱测井资料计算泥质含量。
minmax min TH TH TH TH SH --= ………………………………..………(11 ) 1212--=⋅GCUR SH GCUR sh V (12)式中,TH -目的层钍曲线测井值;THmin -目的层段纯地层钍曲线值;THmax -目的层段泥岩钍曲线值;SH -目的层钍曲线相对值;GCUR -新、老地层校正系数,新地层为3.7,老地层为2.0。
1.7.2 钾曲线(K )minmax min K K K K SH --= ………………………………..….….(13) 1212--=⋅GCUR SH GCUR sh V (14)式中,K -目的层钾曲线测井值;Kmin -目的层段纯地层钾曲线值;Kmax -目的层段泥岩钾曲线值;GCUR -新、老地层校正系数,新地层为3.7,老地层为2.0。
1.7.3 无铀曲线(KTH )minmax min KTH KTH KTH KTH SH --= ……………………………….(15) 1212--=⋅GCUR SH GCUR sh V (16)式中,KTH -目的层无铀曲线测井值;KTHmin -目的层段纯地层无铀曲线值;KTHmax -目的层段泥岩无铀曲线值;GCUR -新、老地层校正系数,新地层为3.7,老地层为2.0。
1.8 利用中子测井资料1.8.1 对于低孔隙度地层,设纯地层ΦN =0,且对中子孔隙度作了岩性校正。
NshN sh V ΦΦ= ………………………………………………..(17) 式中,ΦN -目的层中子孔隙度;ΦNsh -目的层段泥岩中子孔隙度。
注:孔隙性地层计算的Vsh 偏高。
1.8.2 当ΦNmin 不为0%时,minmax min N N N N sh Vφφφφ--= …………………………………(18) 2 . 地层孔隙度(φ)计算公式2.1 利用声波时差测井资料2.1.1 怀利(Wylie )公式)(1)(maT f T ma T sh T sh V CP ma T f T ma T DT s --⋅-⋅--=Φ ……………….(19) 式中,Φs -声波计算的孔隙度,小数;Tma 、Tf -分别为岩石骨架声波时差、地层流体声波时差;Vsh -地层泥质含量,小数;CP -声波压实校正系数,可利用岩心分析孔隙度与声波计算孔隙度统计求出,也可利用密度孔隙度与声波孔隙度统计求出。
DT -目的层声波时差测井值。
2.1.2 声波地层因素公式)1(1DTma T x s -⋅=Φ ……………………………....……..(20) 式中,x -经常取值为砂岩1.6,石灰岩1.76,白云岩2.0,x 大致与储层的胶结指数(m )值有关。
2.1.3 Raymer 公式φφf mav v v +-=2)1( ……………………………………………….(21) 式中,v 、v ma 、v f -分别为地层、岩石骨架、孔隙流体的声速。
2.2 利用密度测井资料)(fD ma D sh D ma D sh V f D ma D DEN ma D D --⋅---=Φ ………….….….(22) 式中,ΦD -密度孔隙度,小数;D ma 、D f -分别为岩石骨架密度值、地层流体密度值,g/cm 3;DEN -目的层密度测井值,g/cm 3;Dsh -泥岩密度值,g/cm 3;Vsh -储层泥质含量,小数。
2.3 利用补偿中子测井01.0)5.0(⨯⨯⨯--=Φsh N sh V LCOR CN N ……….(23) 式中,ΦN -中子孔隙度,小数;CN -目的层补偿中子测井值,%;LCOR -岩石骨架中子值,%;Vsh -目的层泥质含量,小数;Nsh -泥岩中子值,%。
2.4 利用中子-密度几何平均值计算222N D Φ+Φ=Φ ……………………… …………..(24) 式中,ΦD 、ΦN -分别为密度、中子孔隙度,小数。
2.5 利用中子伽马测井计算2.5.1 绝对值法NG K A ⋅+=Φlg …………………………… .…….(25) 式中,Φ-中子伽马计算的孔隙度;NG -目的层中子伽马测井值;A 、K -分别为地区性常数、斜率。
说明:在工区内选择两个孔隙度差别较大的地层,分别求出其孔隙度和所对应的中子伽马读数,在半对数坐标纸上,纵坐标为孔隙度,横坐标为中子伽马值,将其作为两个边界点,即可求出A 、K 两个经验系数。
2.5.2 相对值法(古林图版法))1(lg NGoNG K A -⋅+=Φ ………………………………………………(26) 式中,NG -储层中子伽马测井值;NGo-标准层的中子伽马读数。
说明:标准层选择为硬石膏(Φ=1%),其中子伽马值为NGo,在半对数坐标纸上,纵坐标设(1-NG/NGo),横坐标为lgΦ,如果井剖面上有硬石膏层,则读出其NG值(NGo)和目的层的NG值,并知道中子伽马仪器的源距,就可在上述图版上读出其孔隙度。
如果井剖面上没有硬石膏层,则选择距目的层较近的井眼大于40cm的泥岩层作标准层,其中子伽马读数认为是Φ=100%的中子伽马读数NG1,再将其按井径转换图版转换为NGo即可。
转换方法如下:转换图版纵坐标为井径校正系数Kd,Kd=NGo/NG1,横坐标为井径值。
知道目的层的井径值,由图版查得Kd值,则NGo=Kd·NG1,即可求出(1-NG/NGo),查古林图版即可求出Φ。
图1 古林图标3. 地层含水饱和度(Sw )计算3.1 粒状砂岩或少量含泥质砂岩层饱和度公式(Archie ):n m t R w R b a w S /1)(Φ⋅⋅⋅=……………………………….…….(27) 式中,Sw -目的层含水饱和度,小数;Rt -目的层深电阻率测井值,Ω·m ;Φ-目的层孔隙度,小数;Rw -地层水电阻率,Ω·m ;a -岩性附加导电性校正系数,其值与目的层泥质成分、含量及其分布形式密切相关;b -岩性润湿性附加饱和度分布不均匀系数。