首页 测井解释计算常用公式

测井解释计算常用公式

举报
开通vip

测井解释计算常用公式测井解释计算常用公式目录 1. 地层泥质含量(Vsh)计算公式 …………………………….………….. ……..……1 2. 地层孔隙度(υ)计算公式 ………………………………… ..…….…...43. 地层含水饱和度(Sw)计算 …………………………………………………. …7 4. 钻井液电阻率的计算公式 ……………………………………………… .….….…….12 5. 地层水电阻率计算方法 …………………………………………… ….….………..13 6. 确定a、b、m、n参数 ……………………………………...

测井解释计算常用公式
测井解释计算常用公式目录 1. 地层泥质含量(Vsh)计算公式 …………………………….………….. ……..……1 2. 地层孔隙度(υ)计算公式 ………………………………… ..…….…...43. 地层含水饱和度(Sw)计算 …………………………………………………. …7 4. 钻井液电阻率的计算公式 ……………………………………………… .….….…….12 5. 地层水电阻率计算方法 …………………………………………… ….….………..13 6. 确定a、b、m、n参数 ………………………………………………….….……..21 7. 确定烃参数 …………………………………………………...………….....……..24 8. 声波测井孔隙度压实校正系数Cp的确定方法 ……………………………….…....25 9. 束缚水饱和度(Swb)计算 ……………………………………………..…….…....26 10(粒度中值(Md)的计算方法 …………………………………………..…….…...28 11(渗透率的计算方法 …………………………………………………….......….…….29 12. 相对渗透率计算方法 ………………………………………………………………..35 13. 产水率(Fw) ……………………………………………………………..….……..35 14. 驱油效率(DOF)…………………………………………………………….….…..36 15. 计算每米产油指数(PI)…………………………………………………….……...36 16. 中子寿命测井的计算公式 ……………………………………………….………...36 17. 碳氧比(C/O)测井计算公式 ……………………………………………….…......38 18(油层物理计算公式 …………………………………………………………….…..44 19(地层水的苏林分类法 ………………………………………………………….……48 20. 毛管压力曲线的换算 ………………………………………………………….…….48 21. 地层压力 ………………………………………………………………….………….50 附录:石油行业单位换算 …………………………………………………….…………51 测井解释计算常用公式 1. 地层泥质含量(Vsh)计算公式 1.1 利用自然伽马(GR)测井资料 1.1.1 常用公式 GR,GRminSH, ……………………..(1) GR,GRmaxmin式中,SH,自然伽马相对值; GR,目的层自然伽马测井值; GRmin,纯岩性地层的自然伽马测井值; GRmax,纯泥岩地层的自然伽马测井值。 GCUR,SH2,1,V ……………………..….……(2) shGCUR2,1式中,Vsh,泥质含量,小数; GCUR,与地层年代有关的经验系数,新地层取3.7,老地层取2。 1.1.2 自然伽马进行地层密度和泥质密度校正的公式 ,,GR,BboV,sh ………………..…….(3) ,,GR,Bshmaxo式中,ρ、ρ,分别为储层密度值、泥质密度值; bsh Bo,纯地层自然伽马本底数; GR,目的层自然伽马测井值; GRmax,纯泥岩的自然伽马值。 1.1.3 对自然伽马考虑了泥质的粉砂成分的统计方法 ,GR,A,bV, …………………………(4) shB,SI ,C 1,SI式中,SI,泥质的粉砂指数; SI,(Φ,Φ)/Φ …………………...……….(5) NclayNshNclay (Φ、Φ分别为Φ,Φ交会图上粘土点、泥岩点的中子孔隙度) NclayNshND A、B、C,经验系数。 1.2 利用自然电位(SP)测井资料 1 SP,SPminV,,1.0,, ..…….……………(6) shSP,SPmaxmin 式中,SP,目的层自然电位测井值,mV; SPmin,纯地层自然电位值,mV; SPmax,泥岩层自然电位值,mV。 α,自然电位减小系数,α,PSP/SSP。PSP为目的层自然电位异常幅度,SSP 为目的层段纯岩性地层的自然电位异常幅度(静自然电位)。 1.3 利用电阻率测井资料 R,(R,R)1/bshlimt V,[] ………………………..………(7) shR,(R,R)tlimsh 式中,Rlim,目的层井段纯地层最大电阻率值,Ω?m; Rsh,泥岩电阻率,Ω?m; Rt ,目的层电阻率,Ω?m; b,系数,b,1.0,2.0 1.4 中子,声波时差交会计算 V,A/B ………………………………………………….………….(8) sh A,,(T,T),,t(,,1),T,,,T NmafNmamaNmaf B,(T,T)(,,1),(,,1)(T,T) mafNshNmashf 式中,Tma、Tf,分别为岩石骨架声波时差、地层流体声波时差; Φ、Φ,分别为岩石骨架中子值、泥岩中子值,小数; NshNma Δt,目的层声波时差测井值; Φ,目的层中子测井值,小数。 N 1.5 中子,密度交会计算 V,A/B ………………………………………..………………..(9) sh ,,(,,1),,,(,,),,,,,,A bNmaNmaffNmama B,(,,,)(,,1),(,,1)(,,,) shfNmaNshmaf 3式中,ρ、ρ,分别为岩石骨架密度值、地层流体密度值,g/cm; maf Φ、Φ,分别为岩石骨架中子值、泥岩中子值,小数; Nmash 3 ρ,泥岩密度值,g/cm; sh 3 ρ、Φ,目的层密度测井值,g/cm、中子测井值,小数。 Nb 2 1.6 密度,声波交会计算 V,A/B ………………………………………..………………..(10) sh ,(),(,),,A,,t,,t,,t,,,,t,,,t bmafmaffmamafB,,(,,)(,t,,t),(,t,,t)(,,,) shfmafshfmaf 1.7 利用自然伽马能谱测井 1.7.1 钍曲线(TH) 如果有自然伽马能谱测井,则优先选用能谱测井资料计算泥质含量。 TH,THminSH, ………………………………..………(11 ) TH,THmaxmin GCUR,SH2,1,V ………………………………………(12) shGCUR2,1 式中,TH,目的层钍曲线测井值; THmin,目的层段纯地层钍曲线值; THmax,目的层段泥岩钍曲线值; SH,目的层钍曲线相对值; GCUR,新、老地层校正系数,新地层为3.7,老地层为2.0。 1.7.2 钾曲线(K) K,KminSH, ………………………………..….….(13) K,Kmaxmin GCUR,SH2,1,V …………………………….……..(14) shGCUR2,1 式中,K,目的层钾曲线测井值; Kmin,目的层段纯地层钾曲线值; Kmax,目的层段泥岩钾曲线值; GCUR,新、老地层校正系数,新地层为3.7,老地层为2.0。 1.7.3 无铀曲线(KTH) KTH,KTHminSH, ……………………………….(15) KTH,KTHmaxmin GCUR,SH2,1,V ……………………………………(16) shGCUR2,1 式中,KTH,目的层无铀曲线测井值; KTHmin,目的层段纯地层无铀曲线值; 3 KTHmax,目的层段泥岩无铀曲线值; GCUR,新、老地层校正系数,新地层为3.7,老地层为2.0。 1.8 利用中子测井资料 1.8.1 对于低孔隙度地层,设纯地层ΦN,0,且对中子孔隙度作了岩性校正。 ,N V, ………………………………………………..(17) sh,Nsh 式中,Φ,目的层中子孔隙度; N Φ,目的层段泥岩中子孔隙度。 Nsh 注:孔隙性地层计算的Vsh偏高。 1.8.2 当Φ不为0,时, Nmin ,,,NNminV, …………………………………(18) sh,,,NmaxNmin 2 . 地层孔隙度(φ)计算公式 2.1 利用声波时差测井资料 2.1.1 怀利(Wylie)公式 DT,TT,T1mashma,,(),,V,() ……………….(19) sshT,TCPT,Tfmafma 式中,Φs,声波计算的孔隙度,小数; Tma、Tf,分别为岩石骨架声波时差、地层流体声波时差; Vsh,地层泥质含量,小数; CP,声波压实校正系数,可利用岩心分析孔隙度与声波计算孔隙度统计求出, 也可利用密度孔隙度与声波孔隙度统计求出。 DT,目的层声波时差测井值。 2.1.2 声波地层因素公式 T1ma,,,(1,) ……………………………....……..(20) sxDT 式中,x,经常取值为砂岩1.6,石灰岩1.76,白云岩2.0,x大致与储层的胶结指数 (m)值有关。 2.1.3 Raymer公式 2v,v(1,,),v, ……………………………………………….(21) maf 式中,v、v、v,分别为地层、岩石骨架、孔隙流体的声速。 maf 2.2 利用密度测井资料 4 D,DEND,Dmamash,,,V,() ………….….….(22) DshD,DD,Dmafmaf 式中,Φ,密度孔隙度,小数; D3 D、D,分别为岩石骨架密度值、地层流体密度值,g/cm; maf3 DEN,目的层密度测井值,g/cm; 3 Dsh,泥岩密度值,g/cm; Vsh,储层泥质含量,小数。 2.3 利用补偿中子测井 ,,(CN,LCOR,0.5,V,N),0.01 ……….(23) Nshsh 式中,Φ,中子孔隙度,小数; N CN,目的层补偿中子测井值,,; LCOR,岩石骨架中子值,,; Vsh,目的层泥质含量,小数; Nsh,泥岩中子值,,。 2.4 利用中子,密度几何平均值计算 22,,,DN,, ……………………… …………..(24) 2 式中,Φ、Φ,分别为密度、中子孔隙度,小数。 DN 2.5 利用中子伽马测井计算 2.5.1 绝对值法 lg,,A,K,NG …………………………… .…….(25) 式中,Φ,中子伽马计算的孔隙度; NG,目的层中子伽马测井值; A、K,分别为地区性常数、斜率。 说明:在工区内选择两个孔隙度差别较大的地层,分别求出其孔隙度和所对应的中子伽马读数,在半对数坐标纸上,纵坐标为孔隙度,横坐标为中子伽马值,将其作为两个边界点,即可求出A、K两个经验系数。 2.5.2 相对值法(古林图版法) NG lg,,A,K,(1,) ………………………………………………(26) NGo 式中,NG,储层中子伽马测井值; NGo,标准层的中子伽马读数。 说明:标准层选择为硬石膏(Φ,1,),其中子伽马值为NGo,在半对数坐标纸上,纵坐标设(1,NG/NGo),横坐标为lgΦ,如果井剖面上有硬石膏层,则读出其NG值(NGo)和目的层的NG值,并知道中子伽马仪器的源距,就可在上述图版上读出其孔隙度。 如果井剖面上没有硬石膏层,则选择距目的层较近的井眼大于40cm的泥岩层作标准 5 层,其中子伽马读数认为是Φ=100,的中子伽马读数NG,再将其按井径转换图版转换为1NGo即可。 转换方法如下:转换图版纵坐标为井径校正系数Kd,Kd,NGo/NG,横坐标为井径1值。知道目的层的井径值,由图版查得Kd值,则NGo,Kd?NG,即可求出(1,NG/NGo),1查古林图版即可求出Φ。 图1 古林图标 图2 井径校正图版 6 3. 地层含水饱和度(Sw)计算 3.1 粒状砂岩或少量含泥质砂岩层饱和度公式(Archie): a,b,R1/nw ……………………………….…….(27) S,()wmR,,t式中,Sw,目的层含水饱和度,小数; Rt,目的层深电阻率测井值,Ω?m; Φ,目的层孔隙度,小数; Rw,地层水电阻率,Ω?m; a,岩性附加导电性校正系数,其值与目的层泥质成分、含量及其分布形式密 切相关; b,岩性润湿性附加饱和度分布不均匀系数。对于亲水岩石,b<1(在油驱水过 n程中将有残余水存在,形成连续的导电通道,致使Rt/Ro < 1/Sw);对于亲 n油岩石,b>1(油驱水过程将是“活塞式”,而没有残余水存在,Rt/Ro >1/Sw)。 m,孔隙度指数(胶结指数),是岩石骨架与孔隙网混引起的孔隙曲折性的度 量。孔隙曲折度愈高,m值愈大。 n,饱和度指数,是对饱和度微观分布不均匀的校正。由于孔隙的曲折性,在 驱水过程中烃与水在孔隙中的分布是不均匀的,这种不均匀性随Sw变化, 进一步增大了电流在岩石孔隙中流动的曲折性,使Rt的增大速率比Sw降 低的速率大,因此需要利用饱和度指数n进行校正。 注:m和a是互相制约的,a大,m就小,a小,m就大。根据实际井的实验资料, 分别对砂岩和碳酸盐岩研究了m和a之间的定量关系: 地层水含盐量 8500,300000g/L, 孔隙度4,>30,,渗透率1mD以上时, a值在0.3,1.0,砂岩m值在0.5,2.6,碳酸盐岩m值在1.0,2.6。研究结果得 到以下经验关系式: 砂岩: m,1.8,1.29 lga 碳酸盐岩:m,2.03,0.911 lga m值与Φ的经验关系: 砂岩(Φ为20,32,) m,14.4,20.21 lgΦ 碳酸盐岩(Φ为8,18,) m,7.3,6.13 lgΦ 3.2 印度尼西亚公式 c,V12cle2,[,],S …………………………………..(28) wRRaRtclw Vclc,1, 式中, 2 Vcl,粘土含量; 2 Rcl,粘土电阻率,Rcl,Rsh(1,SI),SI为泥质的粉砂指数; Φe,目的层有效孔隙度; Rw,地层水电阻率; a, 岩性附加导电性校正系数; 7 Rt,目的层电阻率; Sw,目的层含水饱和度。 4 注:(27)式适用于地层水矿化度较低(< 3?10mg/L)的地区。 对于Vsh?0.5的泥质砂岩,可简化为下式: V,12cle2 …………………………………..(29) ,[,],SwRRaRtclw 3.3 Simandoux公式 dmnVS,1clewn/2 ,,S, ………………………….(30) wdRRtclaR(1,V)wcl 式中,常取m,n,2,d,1,2,常取d,1。上式可得: aRaRVaRV2wwshwcl()S,,,,, ………………..…(31) wmmmRR,2,2,RRshcletshee 令a,0.8,m,2,上式变为: 2RVV0.45,2weshsh ……………………..……….(32) S,,,[()]w2RRRR,twshshe 3.4 尼日利亚公式 aV,12cle2,[,],S ………………………………….(33) wRRaRtclw 式中,a,1,2 3.5 含分散泥质的泥质砂岩饱和度公式 m2,aR(1q)(RR)V(RR)V,,,wwshsh2wshsh …..….(34) S[],,,wm2R,2R,R,sheshete 式中,Rt,目的层电阻率; Rsh,目的层段泥岩层电阻率; Rw,地层水电阻率; Vsh,目的层泥质含量,小数; Φe,目的层有效孔隙度,小数; m,目的层孔隙度指数(胶结指数); a,岩性附加导电性校正系数; Φe,目的层有效孔隙度。 3.6 Waxman,Smits模型(分散粘土双水模型) 8 *RBQR,(n)tvw(1)S,,w *..................................................(35) SF,Rww a *,F ………………………………………… …………(36) *m,t Cw,2B,3.83(1,0.83e) ……………………………………..(37) ,,CEC,(1,),tGQ, ……………………………………..……(38) v,t 式中,Sw,目的层含水饱和度,小数; Rt,目的层电阻率,Ω?m; Rw,地层水电阻率,Ω?m; Φt,目的层(泥质砂岩)的总孔隙度,小数; * F,孔隙度与泥质砂岩总孔隙度(Φt)相等的纯砂岩的地层因素,即地层水电导 率Cw足够高时,泥质砂岩的地层因素; m*,地层水电导率Cw足够高时确定的泥质砂岩的胶结指数,也可看成为经粘土 校正后的纯砂岩的胶结指数; n*,相当于该岩石不含粘土的饱和度指数,常取n*,2.0; 3 Q,岩石的阳离子交换容量,mmol/cm; v CEC,岩石的阳离子交换能力,mmol/g 干岩样; B,交换阳离子的当量电导率,S?cm3/(mmol?m); 3 ρ,岩石的平均颗粒密度,g/cm; G 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 1 粘土矿物CEC值一览表 粘土矿物 蒙脱石 伊利石 高岭石 绿泥石 CEC值(mmol/g) 0.8,1.5 0.1,0.4 0.03,0.15 ?0 CEC平均值 1.0 0.2 0.03 0 3.7 归一化的W,S方程 R1/n*weS,[] …………………………………………………(39) wtm*R,tt S,R,RwtwshwR,we ………………………………………(.40) QR,(S,Q)Rvnwwtvnwsh 9 ,QV,vshtshQ,,vn ………………………………………….(41) ,Qvsht ,t,,tsh ……………………………………………………………..(42) Vsh n*m,*C,F*,C,,,C,S(C,BQ/S) ……..………….(43) wattwtwvwt ,*mC,,,C,C,BQ …………………………………….(44) washtshshwvsh 式中,Swt,泥质砂岩总含水饱和度,小数; Rt,泥质砂岩电阻率,Ω?m; Vsh、Φt,泥质砂岩的泥质(或粘土)含量,小数; Φt,泥质砂岩总孔隙度,小数,可用密度测井来计算;(因为泥质砂岩中的干粘 3土密度ρ 一般近似于纯砂岩骨架的密度ρma,即约为2.65g/cm,故实际cld 上可认为密度测井不受地层粘土含量的影响。) Φtsh,泥岩的总孔隙度,小数,可用密度测井来计算Φ; Dsh Rw,泥质砂岩自由水电阻率,Ω?m; Rwsh,泥质砂岩中粘土水电阻率,Ω?m; Qvn,归一化的泥质砂岩阳离子交换容量,小数,取值范围0,1.0; 3 Qvsh,与砂岩邻近的泥岩的Qv值,mmol/cm; 3 Qv,泥质砂岩的阳离子交换容量,mmol/cm; m*,地层水电导率Cw足够高时确定的泥质砂岩的胶结指数,也可看成为经粘土 校正后的纯砂岩的胶结指数; n*,相当于该岩石不含粘土的饱和度指数,常取n*,2.0; B,交换阳离子的当量电导率,S?cm3/(mmol?m)。 说明:参数m*、Rw、Rwsh的最佳选取方法是用lgRt,lgΦt与Cwa,Qvn交会图。 图3 归一化W,S方程的参数选择 10 用GR-Z或Vsh交会图来鉴别纯砂岩和泥岩点。在交会图(图3)上通过含水纯砂岩点 (S)并与水层点群相切的直线,可认为是代表纯砂岩线,其斜率应为m*,在Φt,1.0 处的截距应为Rw。同时,与纯砂岩线平行,并过泥岩点(Sh)的直线在Φt,1.0处的截 距应为Rwsh。 根据图3(a)的m*作出的Cw,Qvn交会图(图3,(b))同样可确定Rw和Rwsh -m*值,而且还可用于判断解释层段中粘土矿物的成分是否稳定。如在C,Φ/R的值从 wattCw到Cwash范围内,通过水层和泥岩的点子基本在一条直线上,则表明粘土矿物成分 基本稳定。反之,如果Cwa,Qvn交会图上点子很分散,趋势线弯曲,则可能是粘土矿 物成分发生变化,或者m*、Rw发生变化,说明砂岩和泥岩的参数是不同的。此时,只 有用岩心资料才能找出真正的原因。对于明显偏离趋势线的高Qvn层,必须采用不同组 的参数。 参数n*应由岩心测量得出,一般情况下,对于砂岩可取n*,m*,或n*,m*,0.1; 在碳酸盐岩中,可取n*,2.0。 3.7 双水模型,分散粘土(Clavier et) 图, 含泥质地层的双水模型 图4中,Φf,自由水孔隙度(自由水占地层体积的百分数); Φb,束缚水孔隙度; Φh,油气孔隙度; Φt,总孔隙度。 Swf,自由水饱和度;Swf,Φf/Φt Swb,束缚水饱和度;Swb,Φb/Φt Swt,总含水饱和度;Swt,(Φf,Φb)/Φt 或 Swt,Swf,Swb 双水模型的束缚水已包括湿粘土的水分,同时,地层孔隙中存在自由水和束缚水两 种导电路径相同的溶液。除了地层水的导电性按其矿化度预计的值不同以外,含泥质地 层与同样孔隙度、孔道曲折度及含水饱和度的纯地层具有同样的导电特性,而地层水的 导电性是自由水与束缚水并联所决定的。因此,可采用Archie公式来研究含泥质地层的 导电性。双水模型认为束缚水对含泥质地层导电性有重要影响,并把它看作是一种特殊 11 的导电溶液来考虑(这是与W,S模型的主要区别)。 RRwfwb,1 ……………………..(45) R,[SR,(1,S)R]owbwbwbwb2,t ………………………………..………………….(46) S,R/Rwot 式中,Sw,泥质砂岩含水(自由水)饱和度; Rt、Ct,分别为泥质砂岩电阻率、电导率; Ro,泥质砂岩100,含水时的电阻率; Φt,泥质砂岩总孔隙度,小数; Swb,泥质砂岩束缚水饱和度,小数; Rwf、Cwf,自由水(远离粘土表面未被泥质束缚的全部水,远水)电阻率、电导 率; Rwb、Cwb,束缚水(粘土附近缺乏盐分的水)电阻率、电导率。可选择100,纯 2泥岩处的Rwa作为Rwb,即R,RΦ.。 wbshtsh 注:在实际处理时可根据实际情况选择a、m值。 4. 钻井液电阻率的计算公式 4.1 钻井液电阻率的温度转换公式 T,21.51R,R() ,(?)………………………………………………(47) m2m1T,21.52 ,T6.77,1R,R(),(F) …………………….……………………………(48) m2m1,T6.772 式中,Rm1,T1温度下的钻井液电阻率,Ω?m; Rm2,T2温度下的钻井液电阻率,Ω?m。 ,,注:摄氏温度与华氏温度转换关系: T(F),1.8T(C),32 4.2 D.W.Hilchie 研究的盐水溶液电阻率与其温度间的关系 R(1)(T,x)1R(T), …………………………………………………………(49) T,x [,0.340396lgR(1),641427]x,10 ………………………..…………..(50) 式中,R(1),起始温度为T(1)(?F)时测量的盐水溶液电阻率,Ω?m; R(T),温度为T(?F)时测量的盐水溶液电阻率,Ω?m。 4.3 根据钻井液电阻率计算其滤液电阻率 1.07 R,C,R ……………………………………………………………(51) mfm 式中,Rm,钻井液电阻率,Ω?m; 12 C,与钻井液密度有关的系数,可由表2确定 表2 C值与钻井液密度的对应关系表 3钻井液密度(g/cm) 1.2 1.32 1.44 1.56 1.68 1.92 2.16 C 0.847 0.706 0.584 0.488 0.412 0.380 0.350 4.4 泥饼电阻率 1.65R,0.69R(R/R) ……………………………………..(52) mcmfmmf 式中,Rmc,泥饼电阻率,Ω?m; Rmf,钻井液滤液电阻率,Ω?m。 对于大多数NaCl钻井液,有如下近似公式: R,1.5R …………………………………………………………(53) mcm 4.5 钻井液滤液矿化度计算公式 4.5.1 当已知钻井液滤液电阻率Rmf和所对应的温度T,则可用图6所示的图版确定钻 井液滤液矿化度Pmf。 4.5.2 当已知24?或75?F时的钻井液滤液电阻率R时,可用(53)式计算其矿化 mfN 度。 x P,10mf x,[3.562,lg(R,0.0123)]/0.955 ………………………….(54a) mfN 4.5.3 钻井液密度 ,,1,0.73P (24?/75?F,101.325kPa) …………….…..(54b) w 5. 地层水电阻率计算方法 5.1 利用水分析资料计算地层水电阻率 5.1.1 计算地层水等效NaCl总矿化度Pwe 表3 地层水离子的等效系数(Ki)表 +1+2+2-1-2-2-1离子名称 Na K+1 Ca Mg Cl SO CO HCO 433Ki ……………………………………………………….(55) P,K,P,weii 式中,Pwe,等效NaCl溶液矿化度,ppm; Ki,第i种离子的等效系数; Pi,第i种离子的矿化度,ppm。 各种离子的等效系数可按图5所示图版来确定。图板横坐标为混合液总矿化度,纵坐标 为等效系数(Ki) 5.1.2 根据求出的Pwe值,按NaCl溶液电阻率与矿化度及温度的关系图版(图6),可 13 查出地层水电阻率。 5.1.2 根据等效NaCl溶液矿化度,查图板(图6)确定地层水电阻率Rw。 5 按混合液的总矿化度确定各种离子的等效系数 14 5.1.3 图6 NaCl溶液电阻率与其矿化度及温度的关系 15 由(55)式可导出计算24?或75?F时地层水电阻率R的近似式: wN 5.1.3 近似计算方法 0.995R,0.0123,3647.54/P ………………………..(56) wNwN式中,P,24?或75?F时地层水总矿化度,(NaCl,mg/L); wN R,24?或75?F时地层水电阻率,Ω?m。 wN 计算出RwN后,再利用(57z)或(57b)式计算任意温度(T)下的地层水电阻率 Rw。即 45.5R,R[] ……………………。…………………(57a) wwN,T(C),21.5 81.77R,R[] 或 ………………………。……………….(57b) wwN,T(F),6.77 5.2 利用自然电位计算Rw 5.2.1 厚的纯地层的静自然电位SSP为 RmfeSSP,,Klg ……………………………。…….……….(58) Rwe 式中,K,自然电位系数,其值与温度成正比: , ……………………………………….(59a) K,60,0.133T(F) ,或 …………………………………….(59b) K,70.7[273,T(C)]/298 Rwe,地层水等效电阻率,Ω?m; Rmfe,钻井液滤液等效电阻率,Ω?m。 5.2.2 按测井图头标出的T1温度下的钻井液电阻率Rm计算24?时的钻井液电阻率 T1 R。 mN T,121.5R,R ………………………………..………(60) ()mNmT1,2421.5 5.2.3 按公式(51)计算24?时的钻井液滤液电阻率R。 mfN 1.07 ……………………………………………(60) R,C,RmfNmN 式中,C,根据钻井液密度,按表2查出。 5.2.4 计算24?时的钻井液滤液等效电阻率R。 mfeN R,0.85R 当R,0.1 Ω?m时, …………………….(61) mfNmfeNmfN 146R,5mfNR, 当R?0.1 Ω?m时, ……………..….(62) mfNmfeN337R,77mfN 5.2.5 计算24?时的等效地层水电阻率R。 weN 16 R,RmfeNwe R, weN ………………………………………..(63) Rmfe 5.2.6 计算24?时地层水电阻率R。 wN (0.69R,0.24)weN 当R,0.12 Ω?m时, …….…(64) R,,0.58,10wNwN 77R,5weN 当R?0.12 Ω?m时,R, ……………………(65) wNwN146,337RweN 5.2.7 计算地层温度下的地层水电阻率。 , …………………………………….(66a) R,45.5R/[T(C),21.5]wwN ,或 ……………………………………..(66b) R,81.77R/[T(F),6.77]wwN 注意:用自然电位计算Rw的方法,适用于地层水主要含NaCl和从SP曲线能得到好的 静自然电位SSP值的情况。如果不能满足上述条件,则需对SP曲线运用专门的 图版进行(地层厚度、井径、侵入带及电阻率比值(Ri/Rm)等校正,从而得到 SSP。如果钻井液与地层间压差过大,SP中明显存在过滤电位成分,则用SP计 算的Rw可能偏低。 5.3 视地层水电阻率法 m,,Rt,/,RRF ……………………………………………(67) wata 式中,Rwa,视地层水电阻率,Ω?m; Rt,深探测电阻率,Ω?m;(Rt应为具有一定厚度的纯岩性水层的Ro) Φ,地层孔隙度,小数; m,胶结指数; a,岩性附加电阻率校正系数。 说明:在具有较厚的纯水层井段和Rw基本稳定或Rw逐渐变化的层段,选择纯水层的 Rwa作为Rw,可取得较好的效果。 5.4 用Rt和Rxo确定Rw 具有均匀粒间孔隙的纯地层,由Archie公式可分别导出Sw和Sxo关系式,将两式 合并可得: SR,Rnwxow(), ………………………….……………………..(68) SR,Rxotmf 在有钻井液侵入的含水纯砂岩处,Sw,Sxo,1,故 Rw/Rmf,Rt/Rxo,因此有 R,R,R/R ……………………………………………….(69) wtmfxo 5.5 电阻率,孔隙度交会图法 5.5.1 Hingle交会图法 对于均匀粒间孔隙的纯地层,由Archie公式可得 17 nS11/mw,(), …………………………………………….(70) nabRRwt 1 nRt 图7 Hingle电阻率,孔隙度交会图 对于给定地区和岩性,系数a、b和指数m、n是已知的。在岩性和Rw基本保持基本不 变的解释井段内,对给定的含水饱和度Sw,令 nS1/mw , (A为常数) (),A abRw 111用按刻度的坐标轴作轴,用线性刻度轴作x,υ轴,则在,υ交会 y,nnnRRRttt图上,方程(70)就成为直线方程y,Ax,而且该直线过原点,即骨架点(υ,0,Rt,?),取不同的Sw值,就得到不同的直线,从而得到用Sw刻度的Rt,Φ交会图(如图7所示)。 18 可按地区经验选取a、b、m、n值。一般取n,2,b,1。对砂岩取a,0.62,m,2.15;对碳酸盐岩取a,1, m,2。 在Hingle交会图上,对于100,含水层,Sw,1,Rt,Ro,如令a=1,m,2, 则有 11 …………………………………….(71) ,, RRtw 在Hingle交会图上100,含水层就是左上方的一条直线,其斜率为。由此可得出确1/Rw定Rw的方法。即在解释井段上绘制Hjngle交会图或频率交会图及GR-Z图,找出岩性纯,足够厚,无油气显示的纯水层,这些纯水层同原点的连线即为100,含水线,在水线上任取 2R,R/F,R,一点,则 。 woo 说明:Hingle交会图的横轴可以选用孔隙度、声波时差、密度或中子测井值,且为线性 刻度。这些交会图的原点均为骨架点(Φ,0,Rt,?)。因此,根据100,含水线与Rt,?线的交点就可以求得骨架矿物的参数(Δt、ρ、Φ)。知道了Δt、ρ、ΦmamamaNmamaNma m就可以按υ或F的单位,对Δt、ρ、Φ的刻度重新刻度。用已确定的F(,a/Φ)刻度,bN 可以计算Rw,而且按类似的方法画出Sw为常数的直线。应用这些Rt,Φ交会图确定Rw、Sw和判断油水层的关键是要正确确定水线位置。因此,此法要求在解释层段上,要有若干个纯含水层,地层水电阻率稳定,岩性不变和侵入不深,要求孔隙度变化范围相当大,并且所测参数(Δt、ρ、Φ)与Φ呈线性关系,所用的a、b、m、n等参数符合本区地质条件。 bN 5.5.2 Pickett交会图法 在Archie公式中,令b,1,则有 aRRnwoS,,w mR,Rtt 对该式两边取对数得: lgR,,mlg,,lg(aR),nlgS tww …………………………….(72) 在水层处,Sw,100,,(72)式可简化为: lgRt,,mlg,,lg(aRw) …………………………….(73) 令y,lgRt,x,lgΦ,c,lg(aRw),则在双 对数坐标中,(73)式即为一条直线, y,mx,c 斜率为m。这种在双对数坐标中绘制的Rt,Φ 交会图即为Pickett交会图,如图8所示。 19 图8 Pickett 电阻率,孔隙度交会图 100,含水线在Φ,100,的纵坐标轴上的截距为aRw,设a,1,则可求出Rw。 这种交会图的优点是不需要知道m值,而且由水线的斜率可确定m值。同样在此 交 会图上可画出不同Sw值的直线,它们均平行于水线。该交会图的孔隙度轴(横轴)也可采用Δt、ρ、Φ。 bN 5.6 利用泥岩层近似估计地层水电阻率 在储集层与其邻近泥岩具有相同或相近的地层水的地区,可用泥岩层估计地层水电阻率。 此法不适用于致密泥岩层、含油气泥岩以及井壁垮塌的泥岩。 mR,R/, ……………………………………….(74) shwsh 式中,Rsh,泥岩电阻率; Φsh,泥岩孔隙度,用声波测井资料计算。 m,胶结指数,按实际资料统计,Pickett认为是Rw的函数; m,2.54,2.62R w lgR,2.62Rlg,,lgR,2.54lg, ……………(75) wwshshsh Schlumberger公司的Tixier对泥岩电阻率Rsh和声波时差Δtsh(μs/ft)进行统计,得出 如下经验关系: R,R(,t,230)/1640 ………………………….(75) wshsh 5.7 由地区统计规律确定Rw lgR,C,D,A …………………………………….(76) w 式中,D,地层深度,m; C、A,与地区有关的经验系数。 5.8 泥质砂岩等效地层水电阻率计算方法 C,(1,aQV)C,,Q …………………………..(77) wevQwfv 式中,C,泥质砂岩等效地层水电导率; we3 Q,泥质砂岩阳离子交换容量,mmol/cm; v ,,CEC(1,)tGQ, ………………………………………..(78) v,t CEC,泥质砂岩阳离子交换能力,mmol/g 干岩样; Φt,泥质砂岩总孔隙度,小数; 3 ρ,岩石平均颗粒密度,g/cm; G, a,Na离子扩散层的扩散因子; 1 当Pw,Pwo a, …………………….(79) Pwo/Pw 当Pw?Pwo 20 Pw,地层水矿化度; 3,x时的地层水矿化度,约为0.35mol/cm; Pwo,xdH,,8 X,Na离子扩散层厚度,10cm; d3 V(,γx),Q,1mmol/cm时粘土水占据的孔隙度; Qdv ,V,1/[2.853,0.019T(C)] ………………………(80) Q β,粘土水中补偿离子Na,的等效电导率,(S/m)(mmol/L) ,,0.0857T(C),0.143, …………………………(81) Cwf,自由水电导率。 注:选取地层水的原则是:若本井或邻井有可靠的水分析资料,则应首先采用水分析资 料计算Rw;如有分区分层位的准确所Rw资料,而本井的电阻率和SP又无异常显 示,则可采用分区分层位选用的Rw数据;否则,应采用多种方法计算,选择其中 合适的值(一般是最小的)作为Rw,使最终计算的Sw和Sh符合地质情况及测井 显示。 6(确定a、b、m、n参数 6.1 实验室测量 mF,RR,a, 依据Archie公式: ………………………..(82) ow nnI,R/FR,b/S,b/(1,S) ………….(83) twwh 在本地区选择同类岩性的若干块标准岩样,在101.325KPa压力下,分别测量在100 ,饱和盐水时的电阻率Ro与在不同含水饱和度Sw时的Rt及相应的Φ值,在双对数坐 标上分别绘出F,Φ和I,Sw关系线,如图9和图10所示。 图9 F,Φ关系图 图10 I,Sw关系图 21 由式(82)、(83)可知: lgF,lg(R/R),lga,mlg, ………………..……….(84a) ow lgI,lg(R/R),logb,nlgS ………….………...(84b) tow 在双对数坐标上,F,Φ关系为一条直线,其斜率为m,截距为a。同样,I,Sw关 系也为一条直线,其斜率为n,截距为b。据统计资料,a值范围为0.5,1.5,m值范围 为1.5,3.0,b值接近1,n值范围为1.15,2.2。 表4 勘探初期,无实验条件时参考值 岩性 砂岩 碳酸盐岩 公式 Humblr Tixier 一般 Φ,10, 低孔隙度(Shell) a 1 0.62 0.81 1 1 1 m 2 2.15 2 2 2.1 1.87+0.019/υ,m>4,取m=4 b 1 n 2 裂缝发育的碳酸盐岩,应采用混合孔隙结构指数: mmfb,,,,,,lg{()[(1)]}m, ……………………………(85) lg, 式中,m,裂缝孔隙结构指数,一般为1.1,1.3; f m,粒间孔隙结构指数,一般取2; b γ,υ/υ,裂缝孔隙υ占总孔隙度υt的百分数。 f tf 6.2 根据纯水层测井资料确定a和m 6.2.1 lgF,lgΦ交会图法 选择完全含水、岩性较纯的地层,作F,Φ交会图,如图11所示。由式(84a)可 知,图上代表纯水层直线的斜率等于m,该直线在Φ,100,的纵坐标应为a。 lgF,lga,mlg, 6.2.2 Pickett电阻率,孔隙度交会图法 如图8所示,对纯水层作Rt,Φ交会图,由式(73)可知,该图上代表Sw,100, 纯水线的斜率为m,当Φ,100,时,Ro,aRw,已知Rw,即可求a。 lgR,,mlg,,lg(aR) tw 6.2.3 在纯水层较少时,如已知Rw、Ro、Φ时,令a,1,则 m,(lgR,lgR)/lg, ……………………….(86) wo 注:此法计算的m值,对一个地区的同一岩性,在Φ变化范围不大时,是可行的。 6.2.4 m与a的经验关系 (墨西哥资料) ,222 地层水含盐量8500,300000mg/L,Φ:4,,,30,,K:1?10μm,1μm。 22 砂岩:a,0.3,1.0,m,0.5,2.6;碳酸盐岩:m,1.0,2.6 m,1.8,1.29lga砂 岩: ……………………………………………(87a) m,2.03,0.91lga碳酸盐岩: ………………………………….………..(87b) 图11 典型的F,Φ交会图 6.2.5 m与Φ的经验关系(墨西哥资料) m,14.4,20.21lg, 砂岩(Φ:20,,32,) ……………………(88) m,7.3,6.13lg,碳酸盐岩(Φ:8,,18,) ……………………(89) 6.2.6 m与Φ、K的关系 m是孔隙结构指数,故与岩石的渗透率(K)、孔隙度(Φ)有密切关系,胜利油田 在实验测量的基础上,得出如下经验关系: lgm,0.34,0.12,,0.023lgK …………………….(90) 注:式(90)明确地指出了m的地质物理意义,指明它取决于地层的孔隙结构。M 并不是常数,而是随Φ与K的增加而减小。 表5 相对于纯地层的m值,其它岩性的m值变化趋势 岩性 泥质砂岩、粉砂岩 含钙砂岩 裂缝性碳酸盐岩 Φ、K Φ中、高 K低 Φ低 K低 Φ低 K高 m值 高 高 小(1.1,1.5) 23 6.3 确定n值的经验关系 nn 由式(83):I,R/R,R/FR,b/S,b/(1,S) 知b和n主要反映油气水在孔隙 totwwh 中的分布对岩石电阻率的影响。大量实验资料表明,b接近1,一般取b,1。 饱和度指数n主要同岩性、油气在孔隙中的分布与连通情况、油气与地层水间的表 面张力以及岩石的润湿性等有关。一般来说,随着岩石固结程度的增加,n值也增大, 如n从固结砂岩的1.5,增加到致密砂岩的2.2,碳酸盐岩的n值常常是不同的,其值随 Sw而变化。 油气在在孔隙中的连通情况及岩石的润湿性对n值有很大影响,亲水岩石中水附着 于颗粒表面;而亲油岩石中岩石颗粒表面形成的油膜将会降低地层水的导电性,甚至使 部分地层水成为绝缘的,故亲油岩石的n值相对较高,n常在2,5之间。一般讲,孔隙 连通性变好,会导致亲油、亲水岩石的n值变小;孔道变窄,将使亲水岩石的n值变小。 [注:摘自《测井数据处理与综合解释》主编 雍世和 张超谟] 表6 M.H.Dorfman 关于m、n实验值的统计表 m值 1.52,2.02 多数为1.8,2.0 n 值 1.15,3.8 1.6,2.2 2.0 地层个数 / % 36 / 100 24 / 66.7 3 / 8.3 注:该表是M.H.Dorfman教授对美国主要油田的36个不同岩性岩石(固结的和未固结的砂岩、灰岩、白云岩)的m、n的实验值的统计表。 由表6可见,一般不能把n值看成是等于2的常数。 某油田对2口井63块岩心(纯的和泥质的细砂岩),用Rw,0.06、0.2、2.0Ω?m三种 地层水饱和,测量n、Φ、K,得出如下经验关系: n,1.347,0.519lgR 平均误差为15.7, ………………(91) w n,0.904,0.515lgR,0.325lgK 平均误差为9.7, ……..………….(92) w n,0.14626,0.21046lnR,0.11672lnK,0.0367,,0.00313SP,0.09279Fw 平均误差为8.9, ………………...(93) 式中,Rw,24?时的地层水电阻率,Ω?m; K,渗透率,mD; Φ,孔隙度,,; SP,自然电位,mV: F,地层因素。 由墨西哥98口井资料得出的经验关系为: n,1.095,0.4421lgP ……………………………………..(94) w 式中,Pw,地层水含盐量,g/L。 7(确定烃参数 7.1 烃密度 0.7(1,,)S,(1,,)hr,, …………………………………(95a) o(1,0.72,)Shr 24 ,,(1,0.72)S,(1,)hr,, ……………………………………(.95b) g(2.2,0.8,)Shr 式中,α,Φ/Φ ,Φ,Φ,VΦ,Φ,Φ,VΦ ; NcorDcorNcorNshNshDcorDshDsh Φ、Φ,中子、密度测井视孔隙度; ND Φ、Φ,泥岩的中子、密度测井视孔隙度; NshDsh S,冲洗带残余油气饱和度; hr Vsh,地层泥质含量。 33说明:应先用式(95b)计算,若ρ,0.3g/cm,则以此为结果;如果ρ,0.6g/cm,则 hh 应按式(95a)计算出油的密度;如果0.3?ρ?0.6g/cm3,则最好取(95a)与(95b) h 式计算结果的平均值作为ρ。 h 7.2 Atlas公司公式(1985年) ,,,1,S(1.87,0.72)hr,, (补偿中子)………………….(96) h(2.67,0.75,)Shr ,,,1,S(1.17,0.72)hr,, (井壁中子)…………………..(97) h(1.67,0.75,)Shr 式中,α,同(95)式。 7.3 Schlumberger公司公式 ,,,(1,),0.21.07(1.11,0.15P)mfmfmfmf,,1,S[A(,]hr,,,(1,P),mfmfmamf,, ………(98) h,1.5A1.23[]S,hr,,P,,,(1)mfmfmamf式中,α,同(95)式; A,与挖掘效应有关的系数,对CNL,A,1.3,对SNP,A,1; ρ、P,泥浆滤液密度、矿化度;可分别用(53a)、(53b)计算。 mfmf 8. 声波测井孔隙度压实校正系数Cp的确定方法 8.1 声波测井孔隙度与岩心分析孔隙度对比 根据岩心分析资料确定Δt与Φ的对应关系: 心 ,t,A,,,t ma ,t,,t,t,,t1mama,,, ,A,t,,tCfmap AC,p …………………………………………..(99) ,t,,tfma 25 8.2 Cp与地层深度的统计关系 对一个油田的不同含油层段进行统计,找出每个层段的Cp与其平均埋深的统计关 系: C,B,A,D ……………………………………………….(100) p 式中,D,地层深度,m。 8.3 声波孔隙度与密度孔隙度对比确定Cp 对于较纯的砂岩,用密度测井计算的孔隙度Φ可认为是岩石的有效孔隙度,因此, D 对这些较纯砂岩分别计算出Φ和Φs(,(Δt,Δt)/(Δt,Δt),则 Dmafma C,,/, ……………………..…………………………(101) psD 9. 束缚水饱和度(Swb)计算 岩石的束缚水包括微毛细管孔隙中不流动的水,其它毛管孔隙中细小孔道弯曲处不能流动的滞水和亲水岩石颗粒表面的薄膜滞水。由此可见,束缚水饱和度的大小必然取决于岩石孔隙的几何形状及其毛细管力的大小、岩石的粘土含量和对水的润湿性,而在亲水为主的砂岩中,岩石比面就是这些因素的综合反映,岩比面越大,束缚水饱和度越高、 9.1 CYBELOOK程序计算Swb的方法 图12 泥质指数Ish与束缚水饱和度Swb的关系 26 S,(2I,y)/(2,y) 当Ish,0.5时 …………………….(102a) wbsh S,2I/(2,y) 当Ish,0.5时 …………………… ….(102b) wbsh 2y,3(1,I)/(5,)式中, ; shmax I,(GR,GR)/(GR,GR); shminmaxmin Φmax,地层最大孔隙度,小数; GR,自然伽马测井值; GRmax,泥岩自然伽马值; GRmin,纯地层自然伽马值。 R1wb0.5S,()wba ……………………………………….(103) ,Rtt 式中,Swba,视束缚水饱和度; Rwb,泥岩束缚水电阻率; Rt,地层电阻率(深探测电阻率)测井值; Φt,地层总孔隙度,小数。 将(102)式计算的Swb与(103)式计算的Swba进行比较: 若 Swb,Swba , 则用Swb; 若 Swb,Swba, 则认为Swb不可靠,需降低Ish,重新计算Swb,直到Swb?Swba 为止。 9.2 用孔隙度、泥质含量和润湿性建立Swb经验关系 ,100,[,lg()],Sa ………………………………(104) wb,Vbcsh 式中,a、b、c,与岩性有关的经验系数。一个典型的表达式为 100,S,[1.145,lg(,0.25)] ………………………….(105) wbV3.228sh 式中,当Φ/Vsh ,0.26 时,取Φ/Vsh,0.26,如果计算的Swb?15,,取Swb,15,。 9.3 多功能分析程序计算Swb方法(曾文冲) 9.3.1 υ?20,的砂岩地层: ,lgS,A,(AlgM,A)lg ……………………..(106a) wbod12A3 ,AMAlg,1d2S(),A'wbo或 ……………………………..(106b) A3 A,lgA'其中 oo 27 式中,Ao、A、A、A,经验系数,A、A近似为常数,A1?1.5,A2?3.6,Ao、A3 12312 与地区地质特点有关,Ao的取值范围为0.18,0.36,随胶结程度变弱和孔隙 度增大而减小;A3的取值范围为0.08,0.2,随胶结程度变弱和孔隙度增大 以及亲水性增强而增大。A3对Swb的计算结果影响较大。 Md,粒度中值; 9.3.2 Φ,20,的砂岩地层: 1,,lg(1,S),Bo,(BlgM,B) …………………...(107 a) wbd12B3 ,1,(1,S),(),(BlgM,B),B'或 …………….(107b) wb1d2oB3 B,lgB'其中 oo 式中,Bo、B1、B2、B3,经验系数。其中B1可视为常数,B1?9.8;Bo,0,0.15, B2,3.3,1.0。若取Bo,0,则取B2,3.3。B3,0.7,0.8,它与砂岩的压实程 度和润湿性有关。一般随地层的压实和亲油性增大而增大。 10. 粒度中值(Md)的计算方法 10.1 利用自然伽马相对值计算 最理想的是用自然伽马能谱测井资料计算。 lgM,C,C,,GR (C1,0)……………………(108) do1 式中,ΔGR,(GR,GRmin)/(GRmax,GRmin) Co、C1,经验系数。 C,lgM (当GR,GRmin时) odo Mdo,GR,GRmin的相应层段平均粒度中值。 C,,1.75,C 1o 10.2 利用自然电位计算 NN11,,,, ………………………………….(109) MNNdoo1,,2,,SP lgM,C,C,,SP或 (C1,0)………..…………..…………(110) do1 式中,ΔSP—自然电位相对值; α,自然电位减小系数。 SPSP,shSP,,,1,, ……………………..……………(111) SPSP,sdsh SPsh、SPsd,分别为泥岩、纯砂岩的自然电位数值; 28 Co、C1,与式(108)中的系数相同; No、N1,按方程式(112)求解: N,N,M o1do N1 N,,0.02 …………………………….(112) o2 11. 渗透率的计算方法 表7 某井沙二段上油组单相关系数表 表7是某井沙二油组各种参数的单相关系数表,由表中可见,渗透率与粒度中值的相关 系数为0.839, 说明二者的关系密切,且渗透率随粒度中值的增大而增加。从表中还可看出,影响渗透率的的主要地质因素为粒度中值、粘土含量和束缚水饱和度,与孔隙度也有一定关系。与影响束缚水饱和度的情况相反,岩石比面越大,即表明其孔隙越小,孔隙结构越复杂,孔隙表面对流体流动的阻滞作用越大,渗透率越小。 11.1 用电阻率计算渗透率 dK,CR …………………………………………..(113) t 式中,C、d,按区域及层位统计确定。 11.2 用孔隙度和束缚水饱和度确定渗透率 11.2.1 一般关系式 x,KC, ……………………………………………(114a) ySwb,32 式中,K,渗透率,10μm; Swb,束缚水饱和度,%; υ,孔隙度,,。 11.2.2 Timur公式(POR程序中使用) 4.40.136,,K, ………………………………………(114b) 2Swb 29 式中,Swb,束缚水饱和度,%; υ,孔隙度,,; ,32 K,绝对渗透率,10μm。 11.3 用声,感组合计算渗透率 ,t,1801.76841.974K,1.34R() ……………………….(115) t100 ,32 式中,K,绝对渗透率,10μm; Δt,声波时差; Rt,深感应电阻率。 11.4 用孔隙度和粒度中值确定渗透率 lgK,D,1.7lgM,7.1lg, ………………………….(116) 1d ,32式中,K,绝对渗透率,10μm; Md,粒度中值,mm; Φ,孔隙度,小数; D,与砂岩的压实程度、胶结物含量和分选性有关,随压实程度增大而增大,随胶1 结物增加和分选性变差而减小,可根据地区经验选用。 的关系表 表8 Cp与D1 11.5 用地层压力测试(FMT、RFT)资料计算渗透率 11.5.1 压降法(MFT仪器) Q,,K,c,1842, ………………………..(117) d(P,P),dsf式中,Kd,地层压降有效渗透率,mD; Q,流量,ml/s;(FMT预测试取体积10ml除以流动时间S); μ,地层条件下流体样品粘度,cp; Ps,地层管井压力,psi; Pf,地层流动压力,psi; d,取样探测器内孔直径,in; C,流动模型因数,准球形流动条件下,8in井眼,取c,0.75; 11.5.2 压降法(RFT仪器) q,,K,,5660 ……………………………………..(118) d,Pss式中,Kd,地层压降法有效渗透率,mD; q,流量,ml/s; 30 ΔPss,取样或预测试期间样品流动的压力降,psi; μ,样品在地层条件下的粘度,cp。 注:压降法求出的渗透率,只是地层测试器附近几厘米处地层渗透率的微观指示,往往 进入探测器的只是泥浆滤液,它所控制的区域因受泥浆污染带表皮效应的影响,所求出 的渗透率往往偏低。 11.5.3 压力恢复分析法 图13 圆柱形压力恢复(Horner)曲线直线段斜率 11.5.3.1 FMT圆柱形压力恢复方程: ,,q,ttP,P,,88.4,()log() ...............................(119) wsi,Khtc 式中,Pi,地层压力,psi; Pws,关井后探测器压力,psi; q,压降期间流量,ml/s; μ,地层流体粘度,cp; Kc,地层圆柱形恢复渗透率,mD; h,不渗透层的间距(地层厚度),ft; t,预测试流动时间,s; Δt,关井后的时间,s。 t,,t() 以(Pi,Pws)为纵轴,以对数刻度的 为横轴,绘制Horner曲线图(见图 ,t 31 13)。在圆柱形压力恢复期间,压力数据在图上试线性分布的,则曲线的斜率Mc是可以 计算的,以psi/周期为单位,圆柱形压力恢复渗透率Kc则为: ,q,K88.4(), ………………………………….(120) cMh,c 式中,h,通常取0.5ft。对该值可作适当调节,使计算的Kc值与其它方法计算的值相吻 合。 11.5.3.2 RFT圆柱形压力恢复方程 RFT有两个预测室(各10ml),其预测试流动时间为T、T。 12 q,,iP,P,2867,,f(,t) ………………..(121) iccK,hr T,T,,tqT,,t1222f(,t),log(),log()式中, ………………….(122) cT,,tq,t21 Pi,初始地层压力,psi; Pc,关井后Δt时的地层压力,psi; q,与测试某一阶段的流量,ml/s; i μ,地层流体粘度,cp; Kr,地层圆柱形恢复渗透率(水平渗透率),mD; h,地层厚度,ft; Ti,预测试某一压降阶段的流动时间,s; Δt,关井后的时间,s。 11.5.4 球形径向压力恢复 1/2,,,q,(c)1P,P,8.0,104,f(,t)iwss ………(123) 3/2Ks P,Pq,,,c,,4iws1,,,,8.010(或 …………………..…(124) f,tKK()sss qqqq,/(/1)12121f,t,,,()式中, ……………..(125) s,tT,,tT,T,,t212 11f,t,,(),0, T2,0,则) ( 对于FMT仪器,q2s,tT,,t1 Pi,初始地层压力,psi; Pws,关井后探测器压力,psi; q,第一预测试压降阶段流量,ml/s; 1 q,第二预测试压降阶段流量,ml/s; 2 μ,地层流体粘度,cp; c,未污染地层流体总压缩系数,1/psi; Ks,地层球形渗透率,mD; 32 Φ,地层孔隙度,小数; T1,第一预测试流动时间,s; T2,第二预测试流动时间,s; Δt,关井后的时间增量,s。 1/2 在直角坐标中作(Pi,Pws)与相应f (Δt)关系曲线,其斜率m(单位为psi/s)。 ss 在C和Φ已知的情况下,用下式计算Ks: q2/31/31K,1856,,,()(C,,)st …………………..(126) ms 32K,K,K …………………………………………….(127) sHv 式中,Ks,球形流动状态下渗透率,mD; K,水平渗透率; H Kv,垂直渗透率。 图14 已知各向异性系数d和Ks,求K和Kv图版 H 33 A,K/K …………………………………………(128) vH 式中,A,各向异性系数,一般A,1。 Kv,垂直渗透率; K,水平渗透率。 H 由此可见,知道K s和A,即可求出Kv和K。图14即为已知各向异性系数d和KsH 的情况下求解Kv和K的图版。 H 11.6 核磁测井计算渗透率 11.6.1 当孔隙度和渗透率比较高时: 22K,CT, ………………………………………..(129a) 2 当孔隙度和渗透率比较低时: 2K,CT, ……………………………………………(129b) 2 式中,C,常数; T,核磁测井的T分布中的T值; 222 Φ,地层孔隙度。 11.6.2 Pyammer积分公式 T2bulk 2kPTTdT,()()222, ……………………..(130a) T,02 N 42K,(,T),ii,或改写为 …………………………………(130b) i,1 式中,P(T),T2的分布概率; 2 T,自由态流体横向弛豫时间; 2bulk N,指数个数。 11.6.3 Coates模型 MBVMbaK,C,,(), …………………………….(131) eMBVI Tmax MBVM,,,S(T)d(T)式中, (可动流体孔隙度)………(132) f22, Tcutoff T2cutoff ()(),22 (毛管束缚流体孔隙度)……(133) MBVISTdT, 4 Φe,有效孔隙度; C,与岩性有关的系数,通过岩心实验得到; 34 a、b,经验指数,通常a,4, b,2。 12. 相对渗透率计算方法 12.1 普遍形式 S,SmwwbK,() ………………………….….…….…(134) rw1,Swb S,SS,Snhwwbwwb[1][1()]K,,, ….…….….(135) ro1,S,S1,S,Swbhrwbhr式中,m、n、h,地区经验系数,与岩性、岩石润湿性和流体粘度比值有关。一般m=3, 4,n,1,2,h,1,2; Sw、Swb,分别为地层含水饱和度、束缚水饱和度; Shr,残余油饱和度。 12.2 彼尔逊经验方程 S,S30.5wwbK,()S ……………….…………….(136a) rww1,Swb S,S2wwb(1)K,, ………….……………………(136b) ro1,Swb 12.3 乘方公式 S,S4wwbK,() ……………………………………(137a) rw,S1w 1,S1,S22ww[1][1()]K,,, ……….…(137b) ro1,S,S1,Swbhrwb 12.4 琼斯方程 S,S2wwbK,() ………………………..……………(138a) rw,S1wb ,S0.92wK,() ……………………………………..(138b) ro,S0.9wb 注:琼斯方程是一般经验方程的特例,即相当于一般方程中Shr,0.1,m,3,n,1,h ,1的特例。 13. 产水率(Fw) Q1wF,,w …………………….(139) ,K,QQrowwo1, K,rwo 35 14. 驱油效率(DOF) S,SwjwbDOF,q ………………………….………….(140) 1,Swb 式中,q,修正系数,其作用在于使计算结果与岩心分析的驱油效率结构向吻合。 Swj,注水后的地层含油饱和度。 15. 计算每米产油指数(PI) lgPI,Ilg(K/,),I …………………………...(141) 1ooo式中,Ko,油的相渗透率; μo,原油粘度; I、I,经验系数,与产层特点有关。 1o 16. 中子寿命测井的计算公式 16.1 含油气纯岩石饱和度 ,,,,,,,,()()mamahS,w ……………………(142) ,(,,),wh -3-1式中,Σ,中子寿命测井测得的热中子宏观俘获截面,10cm; Φ、Sw,分别为地层孔隙度、含水饱和度; -3-1 Σma、Σw、Σh,分别为岩石骨架、地层水、油气的热中子宏观俘获截面,10cm。 16.2 含油气泥质砂岩饱和度 (),(),(),,,,,,,,SS,SmatmahtwbwtwbS, ………(143) wt,(),,,twfh S,SwtwbS,w 自由水饱和度: 1,Swb ,,,(1,S)地层有效孔隙度: etwb S,1,S地层剩余油饱和度: wsw 式中,υt,地层总孔隙度(包括泥质束缚水); Swt,地层总含水饱和度(包括泥质束缚水); Swb,泥质束缚水饱和度,通常近似等于Vsh; Σwf,自由水宏观俘获截面; Σwb,束缚水宏观俘获截面; Σma,所有干固体(骨架、粉砂、干粘土胶体)的宏观俘获截面。 16.3 地层水的宏观俘获截面 2(1.8T,32)PP,,22,0.341,, ……………(144a) w78313417,0.68(1.8T,32) 36 式中,T,地层温度,?; P,泥浆滤液或地层水等效NaCl浓度。 123,,22,或 ………………………..(144b) w59R[82/(1.8T,32),1]w 式中,Rw,地层水电阻率。 -3-3-1通常,地层水Σw的变化范围为22?10,120?10cm。 Σwb的计算方法与Σw相似。 表9 热中子俘获的NaCl等效浓度(Westrn Atlas公司) 物质 等效系数 物质 等效系数 物质 等效系数 NaCl 1 Ca 0.02 Cd 23.7 B 119(121*) S 0.028 Br 0.14 Mg 0.004 I 0.094 HCO 0.01 3Cl 1.65 Li 17.3(20*) SO 0.01 2K 0.05 Gd 495 CO 忽略不计 3注:括号中的数据是该公司不同文献中出现的。 -3-1表10 常见岩石骨架的宏观俘获截面(10cm) Σma 变化范围 Σma 常用范围 砂岩 4,9 8,13 石灰岩 7,12 8,10 白云岩 8,12 8,12 硬石膏 13,22 18,21 岩盐 726* 泥岩 25.2,66.2 35,55 16.4 油气的宏观俘获截面 0.715,,22.3/(1,GOR/22000) ……………………………….(145) o 333式中,GOR,油气比,ft/bbl(0.1781m/m)。 F(1.38,0.238),g,, ………………………………….(146) g256,1.4(1.8T,32) 式中,ρg,天然气密度; 23 F,地层压力,lb/in (6.894?10Pa); T,地层温度,?。 -3-3-1-3-1注:通常油的Σo在18?10,22?10cm之间,大多数油田的Σo约为21?10cm。 -3-1重质原油的Σo可大于21?10cm。 37 17. 碳氧比(C/O)测井计算公式 表11 Atlas公司C/O仪器分段记录的元素伽马射线能谱 元素 Si Ca C O Cl H 非弹性散射伽马能谱1.54, 2.50, 3.17, 4.86, / / (MeV) 1.94 3.30 4.65 6.62 3.17, 4.86, 4.65, 2.01, 俘获伽马能谱(MeV) / / 4.65 6.62 6.66 2.43 Atlas 公司C/O仪器记录下列16条曲线: ISI,硅的非弹性散射伽马射线计数率; ICA,钙的非弹性散射伽马射线计数率; IC,碳的非弹性散射伽马射线计数率; IO,氧的非弹性散射伽马射线计数率; FCC,硅的俘获伽马射线计数率(地层对比曲线,它受孔隙度变化的影响); CAC,钙的俘获伽马射线计数率,反映地层的含钙量,受孔隙度、矿化度的影响; CAPT,俘获伽马射线总计数率; INEL,非弹性散射伽马射线总计数率; CO,C/O测井曲线,由IC/IO计算得到,在岩性稳定、高孔隙度地区,用来计算地层含 油饱和度; SICA,Si/Ca测井曲线,由FCC/CAC计算得出,指示地层岩性。在地层水矿化度不高或 比较稳定的中、高孔隙度地区,可用C/O和SICA测井曲线重叠判断地层含 油性。它对孔隙度不敏感,对地层水矿化度(Cl)敏感。 CASI,Ca/Si测井曲线,由Ca和Si的非弹性散射伽马射线计数率之比(ICA/ISI)得出, 指示地层岩性,对孔隙度不灵敏,不受地层水矿化度(Cl)的影响,可用C/O 和CASI曲线重叠直观判断地层含油性; HCHL,H/Cl测井曲线,由氢和氯的俘获伽马射线计数率之比得出,是良好的地层水矿 化度指示曲线; CIM1,俘获伽马射线总计数率与非弹性散射伽马射线总计数率之比(CAPT/INEL)得出, 是良好的孔隙度指示曲线; CIM2,能窗为3.27,6.26MeV的俘获伽马射线计数率与能窗为3.17,6.62MeV的非弹性 散射伽马射线计数率之比,是良好的孔隙度指示曲线。将CIM1与CIM2重 叠显示,可指示套管尺寸、井径和气层影响; CI,能量为0.4,8.8MeV的俘获射伽马射线总计数率与能量为0.4,8.8MeV的非弹性散 射伽马射线总计数率之比,是良好的孔隙度指示曲线,与补偿中子测井曲线 相似; MSID,热中子衰减,是指示地层的热中子宏观俘获截面(Σ)的曲线。在地层孔隙度海 基会地层水矿化度较高时,用来区分油层和水层。 Schlumberger公司的次生伽马射线能谱测井仪(GST)记录的比值曲线: 38 表12 GST仪记录的比值曲线 曲线符号 名称及用途 测量性质 产额比 其它意义 C/O,指示含油性 非弹性散射 C/C υ,岩性 COR Cl/H,矿化度比值 俘获 Ci/H Sw,Φ,Vsh SIR 指示孔隙度比值 俘获 H/(Si+Ca) 岩性,Vsh PIR 指示岩性的比值 俘获与非弹性散射 Si/(Si+Ca) Φ LIR 指示含铁的比值 俘获 Fe/(Si+Ca) 岩性,Vsh ILR 指示石膏含量 俘获与非弹性散射 S/(Si+Ca) 岩性 AIR 17.1 C/O计算含油饱和度公式(Hertzog R.C.) ,,aS,m(1,),BoccC/O,K ………………(147) b(1,S),m(1,),B,,ooo式中,K,碳和氧的快中子平均反应截面比值,其大小随快中子能量和通量分布而变; Bc、Bo,分别为井眼对碳和氧密度的贡献; a,每立方厘米油中的碳原子数; a,N(,/M), Aoooc b,每立方厘米水中氧原子数; b,N(,/M), Awwwo 23 N,阿佛加特罗常数,6.025?10/mol; A α,每个油分子中的碳原子个数; oc α,每个水分子中的氧原子个数; wo Mo、Mw、Mma,分别为油、水、骨架的相对分子量; m,每立方厘米岩石骨架中的碳原子个数; c m,N(,/M), cAmamac m,每立方厘米岩石中氧原子个数; o m,N(,/M), oAmamao α、α,每个骨架分子中的碳与氧原子个数; co3 ρo、ρw、ρma,分别为油、水、岩石骨架的密度,g/cm; Φ,地层孔隙度。 17.2 大庆油田NP系列C/O仪器饱和度公式 ,(C/O)[m(1,,),,]S, …………………………(148) o,[K,,(C/O)]式中,K,与仪器尺寸、井眼条件有关的常数; m,与地层矿物成分有关的系数; Δ(C/O),目的层与水层C/O的差值。 17.3 比值法 39 对于给定的孔隙度Φ,C/O值与含油饱和度So之间均有近似的线性关系。 C/O,(C/O)wS, ………………………………….(149) o,(C/O) C/O,(C/O)wS,1,S,1, ……………..……..(150) wcoo,(C/O) ,(C/O),(C/O),(C/O)式中, ow (C/O)o、(C/O)w,岩性与孔隙度相同的纯岩石在饱和油与饱和水时的C/O值; C/O,在相同岩性和孔隙度的纯岩石中含油饱和度为So时的C/O值。 选取(C/O)w的方法: ? C/O,Si/Ca 交会法 图15 Atlas公司实验得出的C/O与Si/Ca交会图 由图15可见,不同孔隙度的砂岩、石灰岩、白云岩的水点基本上均落在一条平均水线 上,而且这些岩石饱含油时构成的曲线均与水线平行,并随孔隙度的增加,油、水线间的距离随之增大。设油线与水线的C/O值为(C/O)o和(C/O)w,则对一定孔隙度的地层,油线上的点为: 40 (C/O),(C/O),,(C/O) ow 式中,Δ(C/O),在给定Φ与Si/Ca值时,水线与油线的C/O增量。 图15中水线的方程为: (C/O),K(Si/Ca),L ………………………………………(151) www 式中,Kw,水线斜率,本实验的数据为Kw,,0.8; Lw,当Si/Ca为0时,水线在C/O轴上的截距,本实验的数据为Lw?2.3635。 C/O,0.8Si/Ca,LwS, …………………………………..………(152) o,(C/O) 注:对于岩石骨架中不含钙的地层(如菱铁矿FeCO3、火山灰,含C、Mg、Al、Na等 元素)不宜用(152)式计算So, 对这类地层得出的So过高。而实际地层不不一定 含油。 ? C/O,Ca/Si交会法 (C/O),K(Ca/Si),L …………………………………..(153) wwcawca 式中,Kwca,C/O与Ca/Si交会图中水线斜率,Kwca,1.8; Lwca,当Ca/Si,0时,水线在C/O轴上的截距。 C/O,1.8Ca/Si,LwcaS, ………………….………………….(154) o,(C/O) 注:由于Ca/Si测井曲线对岩性的分辩率较低,故通常只在地层水矿化度较高(大于3 44?10,5?10mg/L)或碳酸盐岩地层中使用。 ? 在岩性和孔隙度较稳定的地区,选择100,含水纯砂岩或泥质少的砂岩的C/O测 井值,并做总孔隙度和泥质校正后,可近似作为(C/O)w。 选取Δ(C/O)的方法 图16 Si/Ca(俘获伽马射线计数率比值)刻度实验 41 ? C/O,Si/Ca交会图(图16): ,(C/O),(C/O),(C/O),(C/O),0.8Si/Ca,2.36 ………….(155) ow ? Δ(C/O),Φ交会图 图17 Δ(C/O),Φ关系图(Atlas) 1.11,(C/O),(C/O),(C/O),0.6, ……….………..(156a) ow /0.8/CO,SiCa,LwS, …………………………….(156b) o1.110.6, 注:对于Φ,35,的地层,Δ(C/O)值约为0.233,可视为常数。 17.4 Lawrence T.D. 推荐公式: (C/O)wS,7.52[1,] 高孔隙度地层 ……………………………(157a) oC/O (C/O)2.36wS,[1,] 其它孔隙度地层 ……………………………(157b) o,C/O 42 含特殊矿物(如菱铁矿)或火山灰地层 COCO/,(/)bwS, ………………………(158) oba, 式中,a、b,待定系数,可根据本地区岩心分析和生产数据确定。 水线斜率与截距计算方法 ? 交会图法 图18 C/O,Si/Ca频率交会图计算水线斜率与截距 选择岩性较纯,具有明显水层和油层的井段,作C/O,Si/Ca频率交会图(如图18), 找出最佳水线。 (C/O),(C/O)21K, ………………………(159) w(Si/Ca),(Si/Ca)21 式中,(C/O),(Si/Ca)、(C/O),(Si/Ca),分别是水线上任意两个资料点的坐标。 1122 L,(C/O),K(Si/Ca) ……………….….……(160) www 取水线上一点的C/O和Si/Ca值代入(160)式,即可求得Lw。 ? 用水层测井值求Kw和Kwca 选取标准水层,如果没有标准水层时,也可选取纯泥岩,在水线斜率Kw或Kwca 不变的情况下,由(161a)或(161b)式计算Lw或Lwca。 L,(C/O),K(Si/Ca) ………….…………(161a) www L(,C/O)K,(Ca/Si) ……………………(161b) wcawwca 式中,(C/O),纯水层或纯泥岩处的C/O测井值; W Si/Ca、Ca/Si,纯水层或纯泥岩处的Si/Ca、Ca/Si测井值; Lw,C/O,Si/Ca交会图上水线的截距; 43 Lwca,C/O,Ca/Si交会图上水线的截距。 注:一般Kw与Kwca受岩性影响较小,数值较稳定,如Kw,,0.8,Kwca,,1.8。水 线的截距Lw和Lwca则受井眼环境(套管、水泥、井内泥浆性能等因素)影响, 因此,应根据不同情况,合理选取Lw和Lwca值。 17.4 计算地层总孔隙度Φt 砂岩: 23 ….(162a) ,,2.926388,21.77518(C/I),58.65146(C/I),56.13097(C/I)t 石灰岩: 23 …..(162b) ,,1.554695,7.012792(C/I),8.713974(C/I),1.855169(C/I)t 式中,C/I,地层俘获伽马射线总计数率与非弹性散射伽马射线总计数率纸比。 17.5 碳氧比测井资料的泥质校正 C/O,V(C/O)shsh(C/O), …………………………………(163a) c1,Vsh Si/Ca,V(Si/Ca)shsh(Si/Ca), …………………………………..(163b) c1,Vsh Ca/Si,V(Ca/Si)shsh(Ca/Si), …………………………………(163c) c1,Vsh 式中,(C/O)sh、(Si/Ca)sh、(Ca/Si)sh,分别为泥质的C/O、Si/Ca、Ca/Si值,可用频率交 会图确定。 17.6 套管、水泥环、油密度、地层水矿化度等影响的校正方法 请参看《测井数据处理与综合解释》(主编 雍世和 张超谟)一书。(P420,422) 18. 油层物理计算公式 18.1 岩石粒度参数的计算 18.1.1 不均匀系数a a,d/d ……………………………………………..(164) 6010 式中,d,粒度组成累积分布曲线上累计重量60,的颗粒直径; 60 d,粒度组成累积分布曲线上累计重量10,的颗粒直径。 10 18.1.2 分选系数S 欧美国家往往以累积重量25,、50,、75,三个特征点,将累积曲线划分为四段, 然后按特拉斯克方程求出分选系数。 d75S, …………………………………………………(165) d25 式中,d75,累积分布曲线上,25,处的粒级直径; d75,累积分布曲线上,75,处的粒级直径。 44 按特拉斯克的规定:S=1,2.5 分选好,S=2.5,4.5 分选中等,S,4.5 分选差。 18.1.3标准偏差σ ()(),,,,,,8416955,,, …………………….(166) 46.6 式中,Φi,第i种粒级处对应于以2为底的对数曲线上所取值。 1 ,,,logd,log 22iidi d,第i种粒级处的颗粒直径。 i 表13 按标准偏差划分的分选等级(福克,沃德 1957) 福克. 沃德标准偏差σ 分选等级 ,0.35 极好 0.35,0.50 好 0.50,0.71 较好 0.71,1.00 中等 1.00,2.00 差 2.00,4.00 很差 ,4.00 极差 18.2 岩石的比面 比面,单位体积岩石内岩石骨架的总表面积或单位体积岩石内总孔隙的的内表面 积。当颗粒间是点接触时,即为所有颗粒的总表面积。 A S, ……………………………………………………..(167a) V 23式中,S,岩石比面,cm/cm, 或1/cm; 2 A,岩石颗粒的总比面积或岩石孔隙的总内比面积,cm; 3 V,岩石外表体积(或视体积),cm。 表14 部分岩石粒径与比面 23岩石种类 粒径(mm) 比面(cm/cm) 23半径R的圆球立方体排列 2R S,8,4,R/(4R),,/2R 砂岩 1,0.25 ,950 细砂岩 0.25,0.1 950,2300 泥质砂岩 0.1,0.01 ,2300 比面的不同表达形式: A S,s …………………………………………….(167b) Vs 23式中,Ss,以岩石骨架体积为基准的比面,cm/cm。 45 AS, ………………………………………………(167c) pVp 23式中,Sp,以岩石孔隙体积为基准的比面,cm/cm。 S,,,S,(1,,)S ………………………………..(168) ps 2,8K,,10 ……………………………………….(169) 222S, 2式中,K,岩石渗透率,μm; Φ,岩石孔隙度,小数; τ,毛管迂曲度 23 S,以岩石外表几何体积为基准的比面,cm/cm。 ,S,7000,取τ,1 ……………………………..…………(170a) K ,S,3600,取τ,1.4 …………………………………….(170b) K 18.3 储层岩石孔隙结构 18.3.1 分选系数(Sp) ,,,,,,()()8416955S,, …………………(171) p46.6 ,,,logd式中, i2i d,孔隙喉道直径; Φ值下标数字,累积分布曲线相应百分数的数值。 18.3.2 歪度(Skp) ,,,,,,,,2,,284165095550S,, ……..(172) kp2,(,,)2(,,,)8416955注:若孔隙大小分布曲线对称,则Skp,0,实际上岩样的Skp值可在?1之间变化, 正值表示曲线有一个粗孔隙度的尾部,即粗歪度;负值表示曲线为细歪度。 18.3.3 峰态(Kp) ,,,955,K ………………………….(173) p2.44(,,,)7525 注:Kp用来量度分布曲线的陡峭程度,即量度分布曲线两个尾部的孔喉直径的展幅与中 央展幅的比值。若分布曲线为正态分布,则Kp,1;若曲线为单峰或双峰曲线,Kp 值可以低到0.6;具有尖峰的曲线的Kp值可从1.5到3.0。 46 18.4 孔隙结构参数 18.4.1` 孔喉比 孔隙与孔隙喉道直径之比。 18.4.2 孔隙配位数 每个孔道所连通的喉道数。如一个孔道与三个喉道相连,则配位数为3。 18.4.3 孔隙迂曲度 用来描述孔隙弯曲程度。流体质点实际流经的路程长度l与岩石外观长度L之比值。可从1.2,2.5间选用。 18.5 孔隙分类 表15 孔隙分类表 毛管孔径孔隙分类 裂缝宽度mm 流体状态 岩石类型 mm 岩石中大裂缝、超毛细管孔隙 ,0.5 ,0.25 在重力作用下自由流动 溶洞、未胶结的 疏松砂层 由于毛细管作用,液体 不能自由流动,若使液毛细管孔隙 0.5-0.0002 0.25-0.0001 体沿毛细管移动,需有一般砂岩 超过重力的外力去克服 毛细管力。 分子间引力很大,若使 微毛细管孔隙 ,0.0002 ,0(0001 液体在孔隙中移动需要泥页岩 非常高的压力梯度。 18.6 孔隙度分级 表16 孔隙度分级表 孔隙度,, 平价 25,20 极好 20,15 好 15,10 中等 10,5 差 5,0 无价值 18.7 储层岩石压缩系数 ,V1pC,, …………………………………………(174) fV,pp 式中,C,岩石压缩系数; f Vp,岩石体积; ΔVp,地层压力降低Δp时,孔隙体积的缩小值。 注:岩石压缩系数的大小,表示了岩石弹性驱油能力的大小,故也称为岩石弹性压缩系数。 -6-1岩石压缩系数一般约为1,2?10MPa。 18.8 储层的弹性采油量(ΔVo) ,V,V,p(C,C,) …………………………………..(175) obfL 式中,C,液体压缩系数; L 47 Φ,岩石孔隙度; C,岩石压缩系数。 f 19. 地层水的苏林分类法 当量比 成因系数 水型 环境 ,,Na,Cl大陆冲刷环境 ,1 硫酸钠型 ,2(地面水) SO,4Na ,1,,,ClNa,Cl大陆环境 ,1 重碳酸钠型 ,2(油、气田水) SO4 ,,ClNa,海洋环境 ,1 氯化镁型 ,2(海水) ,MgNa ,1,,,ClCl,Na深层封闭环境 ,1 氯化钙型 ,2(气田水) Mg ,,,NaCl+-++2,1注:(1) 当时,表明Na除了与Cl离子化合外,还有多余的Na与SO4,2SO4 ++-化合,由于此值小于1,故没有多余的Na与HCO化合了,所以Na与全部阴离子化合后,3最终形成NaSO水型。 24 ,,Na,Cl+-++2,1 (2)当时,表明Na除了与Cl和SO离子化合外,还有多余的Na与4,2SO4 -HCO化合了,最终形成NaHCO水型。 33 ,,ClNa,-+-+2,1 (3)当时,表明Cl除了与Na离子化合外,还有多余的Cl与Mg化合,,2Mg -+2形成MgCl,由于此值小于1,故没有多余的Cl与Ca化合了,所以是MgCl水型。 22 ,,ClNa,-++2-+2,1 (4)当时,表明Cl除了与Na、Mg离子化合外,还有多余的Cl与Ca,2Mg 化合生成CaCl,所以是CaCl水型。 22 20(毛管压力曲线的换算 20.1 将实验室测定结果换算到地下条件 ,,cosRRPP,cRcL …………………………………..(176) ,cos,LL 式中,P,实验室测定的毛管压力; cL P,油藏条件下的毛管压力; cR 48 σ,油藏条件下流体表面张力,mN/m; R θ,油藏条件下流体润湿接触角,?; R σ,实验室条件下流体表面张力,mN/m; L θ,实验室条件下流体接触润湿角,?; L 20.2 将压汞测量的毛管压力P换算为油层条件下油,水毛管压力Pow: Hg ,,,cos25,cos01ooww ….(177) P,P,P,PowHgHgHg,,,cos15480,cos140HgHg 式中,P,油层条件下油,水毛管压力; ow σ,油水界面张力,25mN/m; ow θow,油水界面润湿角,0?; σ,汞表面张力,480mN/m; Hg θ,汞润湿接触角,140?; Hg P,压汞法测定的毛管压力。 Hg 注:实际油藏中油水的毛管压力P仅为压汞法测得的毛管压力P的1/15。 owHg 20.3 将半渗透隔板(水,空气体系)所测得的毛管压力P换算为地下油水毛管压力: wg ,,,cos25,cos01owowP,Pwg,P,Powwgwg …….(178) ,cos3,,72,cos0wgwg 式中, P,油层条件下油,水毛管压力; ow σ,油水界面张力,25mN/m; ow θow,油水界面润湿角,0?; σ,水的表面张力,72mN/m; wg θ,水的润湿接触角,0?; wg P,半渗透隔板法测定的毛管压力。 wg 注:实际油藏中油水毛管压力P仅为半渗透隔板法所测得的毛管压力P的1/3。 owwg 20.4 将压汞法测得的P换算为半渗透隔板法下的气水毛管压力P: Hgwg ,,,cos72,cos01wgwgPPHgP,,,Hgwg …….(179) ,,cos,5480,cos0HgHg 20.5 毛管压力与喉道半径之间的关系 2,cos,P, ……………………………………………….(180) cr 式中,r,毛管压力Pc所对应的喉道半径,μm。 Pc,毛管压力,MPa。 20.6 液柱高度与毛管压力的关系 100PcR,h ………………………………………………(181) ,,,wo 式中,h,油水界面以上湿相(水)液柱高度,m; 49 P,地层条件下(如油,水)的毛管压力,MPa; cR3 ρ、ρ,分别为地层条件下的水、油密度,g/cm。 wo 20.7 确定油(水)饱和度随油水过渡带高度之间的变化关系 100()PScw(),hS …………………………………….(182) w,,,wo 注:直接将Pc(Sw)换算成h(Sw)关系,即可求出油(水)过渡带随含水(油) 饱和度变化关系。 21. 地层压力 21.1 液体的静压力 ,3P,10,gD ………………………………………..(183) Hy 式中,P,液体静压力,MPa; Hy3 ρ,平均流体密度,Kg/m; 2 g,重力加速度,9.8m/s; D,液柱高度,m。 注:液体静压力(静液柱压力、静压、静液压力)是由静液柱对某一点产生的压力。 静液柱横断面的大小与形状对压力的大小无影响,其值为该点到水平面的垂直高 度与流体密度的乘积。 21.2 压力梯度 ,3P,10,g (a) ………………………………..……….(184a) FG 式中,PFG,压力梯度,MPa/m; 采用英制单位为 P,0.433, (b) (psi/ft,淡水)………………………..(184b) FGf 3P,0.465, (盐水,ρ,1.07g/cm) fFGf 3 ρ,流体密度,g/cm。 f 注:压力梯度为单位垂直深度的静液柱压力变化。它受井内液体中天然气溶解度和 含盐量以及温度梯度的影响。 n ,D,(,),ibiP1,1oiG,,o …………………(185) nn10,D,D,,ii ,1,1ii 式中,Go,上覆岩层压力梯度,MPa/m; Po,上覆地层压力,MPa; ΔDi,测井图上所取的深度间隔,m; 3 ρbi,对应于每个深度间隔的地层体积密度平均值,g/cm。 21.3 上覆岩层压力 50 P,D,[(,1,,),,,,,,]g,D,,,g …….…….(186a) omaf P,0.433,,,D 采用英制 ………………………………………(186b) o 式中,Po,上覆岩层压力,MPa; D,液柱垂直深度,m; 3 ρ,岩石骨架密度,Kg/m; ma3 ρ,流体密度,Kg/m; f3 ,上覆地层平均密度,Kg/m。 , 注:上覆岩层压力(或上覆层压力)是地层上面的基岩重量和基岩孔隙空间中流 体(油、气、水)的总重力引起的压力。 21.4 地层压力 P,P,, ……………………………………….(187) of 式中,P,上覆地层压力; o P,地层压力; f ,颗粒间的压力(应力)。 , P,P,D (测量压力单位) ……………………….(188a) fFG P(psia),P,D,15 (绝对压力单位)………(188b) fFG 式中,绝对压力是指包括大气压在内的压力,即加上1Kg/m3(约15psi)。 注:常用压力单位换算关系: 1psi,6.894KPa,0.006894MPa 1psi,0.0680atm,0.07032at 1MPa,9.86923atm,10.19715at,145.0537psi 1atm,1.0332at,14.70psi 1at,0.96787atm,14.22psi,0.09807MPa atm,标准大气压;at, 工程 路基工程安全技术交底工程项目施工成本控制工程量增项单年度零星工程技术标正投影法基本原理 大气压; 2 psi,磅/英寸 21.5 流体密度 P(psi/m)FG,, ………………………………….(189) f1.422 3 式中,ρ,流体密度,g/cm; f P,地层压力梯度,psi/m。 FG 51 附录:石油行业单位换算 长 度 换 算 1米(m)=3.281英尺1千米(km)=0.621英里(mile) (ft)=1.094码(yd) 1英里(mile) 1英寸1厘米(cm)=0.394英寸(in) =1.609千米(km) (fm)=1.829(m) 1英寸(in)=2.54厘米(cm) 1海里(n mile)=1.852千米(km) 1码(yd)=3英尺(ft) 1杆(rad)=16.5英尺(ft) 1英里(mile)=5280英尺(ft) 1海里(n mile)=1.1516英里(mile) 1英尺(ft)=12英寸(in) 面 积 换 算 221平方公里(km)=100公顷(ha)=247.1英亩(acre)=0.386平方英里(mile) 221平方米(m)=10.764平方英尺(ft) 21公顷(ha)=10000平方米(m)=2.471英亩(acre) 221平方英寸(in)=6.452平方厘米(cm) -321英亩(acre)=0.4047公顷(ha)=4.047×10平方公里(km) 2=4047平方米(m) 221平方英尺(ft)=0.093平方米(m) 221平方码(yd)=0.8361平方米(m) 221平方米(m)=10.764平方英尺(ft) 221平方英里(mile)=2.590平方公里(km) 体 积 换 算 1美吉耳(gi)=0.118升(1) 1美品脱(pt)=0.473升(1) 1美夸脱(qt)=0.946升(1) 1美加仑(gal)=3.785升(1) 31桶(bbl)=0.159立方米(m)=42美加仑(gal) 31英亩?英尺,1234立方米(m) 52 33)=16.3871立方厘米(cm) 1立方英寸(in 310亿立方英尺(bcf)=2831.7万立方米(m) 31万亿立方英尺(tcf)=283.17亿立方米(m) 31百万立方英尺(MMcf),2.8317万立方米(m) 31千立方英尺(mcf)=28.317立方米(m) 1英加仑(gal)=4.546升(1) 331立方英尺(ft)=0.0283立方米(m)=28.317升(liter) 331立方米(m)=1000升(liter)=35.315立方英尺(ft)=6.29桶(bbl) 质 量 换 算 1长吨(long ton)=1.016吨(t) 1千克(kg)=2.205磅(lb) 1磅(lb)=0.454千克(kg)[常衡] 1盎司(oz)=28.350克(g) 1短吨(sh.ton)=0.907吨(t)=2000磅(lb) 1吨(t)=1000千克(kg)=2205磅(lb)=1.102短吨(sh.ton) =0.984长吨(long ton) 密 度 换 算 3331磅/(lb/ft)=16.02千克/米(kg/m) API度,141.5/15.5?时的比重,131.5 331磅/英加仑(lb/gal)=99.776千克/米(kg/m) 1波美密度(B),140/15.5?时的比重,130 33331磅/英寸(lb/in)=27679.9千克/米(kg/m) 331磅/美加仑(lb/gal)=119.826千克/米(kg/m) 331磅/(石油)桶(lb/bbl)=2.853千克/米(kg/m) 3333331千克/米(kg/m)=0.001克/厘米(g/cm)=0.0624磅/英尺(lb/ft) 运 动 粘 度 换 算 -42222 1斯(St)=10米/秒(m/s)=1厘米/秒(cm/s) 22-2221英尺/秒(ft/s)=9.29030×10米/秒(m/s) -622221厘斯(cSt)=10;米/秒(m/s)=1毫米/秒(mm/s) 动 力 粘 度 换 算 53 -3帕?秒(Pa?s) 动力粘度 1泊(P),0.1帕?秒(Pa?s) 1厘泊(cP)=10 221磅力秒/英尺(lbf?s/ft)=47.8803帕?秒(Pa?s) 221千克力秒/米(kgf?s、m)=9.80665帕?秒(Pa?s) 力 换 算 1牛顿(N),0.225磅力(lbf)=0.102千克力(kgf) 1千克力(kgf)=9.81牛(N) 1磅力(lbf)=4.45牛顿(N) 1达因(dyn)=10-5牛顿(N) 温 度 换 算 K,5/9(?F+459.67) K=?+273.15 n?=(5/9?n+32) ?F n?F=[(n-32)×5/9]? 1?F=5/9?(温度差) 压 力 换 算 22压力 1巴(bar)=105帕(Pa) 1达因/厘米(dyn/cm)=0.1帕(Pa) 1托(Torr)=133.322帕(Pa) 1毫米汞柱(mmHg)=133.322帕(Pa) 1毫米水柱(mmH2O)=9.80665帕(Pa) 1工程大气压=98.0665千帕(kPa) 2221千帕(kPa)=0.145磅力/英寸(psi)=0.0102千克力/厘米(kgf/cm) =0.0098大气压(atm) 2221磅力/英寸(psi)=6.895千帕(kPa)=0.0703千克力/厘米(kg/cm) =0.0689巴(bar)=0.068大气压(atm) 21物理大气压(atm)=101.325千帕(kPa)=14.696磅/英寸(psi) =1.0333巴(bar) 传 热 系 数 换 算 22221千卡/米?时(kcal/m?h)=1.16279瓦/米(w/m) 221千卡/(米?时??)〔1kcal/(m?h??)〕 22,1.16279瓦/(米?开尔文)〔w/(m?K)〕 221英热单位/(英尺?时??F)〔Btu/(ft?h??F)〕 54 22?开尔文)〔(w/m?K)〕 =5.67826瓦/(米 221米?时??/千卡(m?h??/kcal) 22=0.86000米?开尔文/瓦(m?K/W) 热 导 率 换 算 1千卡(米?时??)〔kcal/(m?h??)〕 ,1.16279瓦/(米?开尔文)〔W/(m?K)〕 1英热单位/(英尺?时??F)〔But/(ft?h??F) ,1.7303瓦/(米?开尔文)〔W/(m?K)〕 比 容 热 换 算 1千卡/(千克??)〔kcal/(kg??)〕 ,1英热单位/(磅??F)〔Btu/(lb??F)〕 ,4186.8焦耳/(千克?开尔文)〔J/(kg?K)〕 热 功 换 算 1卡(cal)=4.1868焦耳(J) 1大卡=4186.75焦耳(J) 1千克力米(kgf?m)=9.80665焦耳(J) 1英热单位(Btu),1055.06焦耳(J) 1千瓦小时(kW?h)=3.6×106焦耳(J) 1英尺磅力(ft?lbf)=1.35582焦耳(J) 1米制马力小时(hp?h)=2.64779×106焦耳(J) 1英马力小时(UKHp?h)=2.68452×106焦耳 1焦耳,0.10204千克?米 ,7,2.778×10千瓦?小时 ,7,3.777×10公制马力小时 ,7,3.723×10英制马力小时 ,4,2.389×10千卡 ,4,9.48×10英热单位 功 率 换 算 55 1英热单位/时(Btu/h)=0.293071瓦(W) 1千克力?米/秒(kgf?m/s)=9.80665瓦(w) 1卡/秒(cal/s),4.1868瓦(W) 1米制马力(hp)=735.499瓦(W) 速 度 换 算 1英里/时(mile/h)=0.44704米/秒(m/s) 1英尺/秒(ft/s)=0.3048米/秒(m/s) 渗 透 率 换 算 2),9.81×107达西 1达西,1000毫达西 1平方厘米(cm 地 温 梯 度 换 算 1?F/100英尺,1.8?/100米(?/m) 1?/公里,2.9?F/英里(?F/mile)=0.055?F/100英尺(?F/ft) 油 气 产 量 换 算 1桶(bbl)=0.14吨(t)(原油,全球平均) 1万亿立方英尺/日(tcf/d) 3=283.2亿立方米/日(m/d) 3=10.336万亿立方米/年(m/a) 10亿立方英尺/日(bcf/d) 3=0.2832亿立方米/日(m/d) 3=103.36亿立方米/年(m/a) 1百万立方英尺/日(MMcf/d) 3=2.832万立方米/日(m/d) 3=1033.55万立方米/年(m/a) 56 1千立方英尺/日(Mcf/d) 3/d) =28.32立方米/日(m 3=1.0336万立米/年(m/a) 1桶/日(bpd)=50吨/年(t/a)(原油,全球平均) 1吨(t)=7.3桶(bbl)(原油,全球平均) 气 油 比 换 算 31立方英尺/桶(cuft/bbl)=0.2067立方米/吨(m/t) 热 值 换 算 1桶原油,5.8×106英热单位(Btu) 1吨煤=2.406×107英热单位(Btu) 1立方米湿气,3.909×104英热单位(Btu) 1千瓦小时水电,1.0235×104英热(Btu) 1立方米干气,3.577×104英热单位(Btu) (以上为1990年美国平均热值) (资料来源:美国国家标准局) 热 当 量 换 算 1桶原油,5800立方英尺天然气(按平均热值计算) 1立方米天然气,1.3300千克标准煤 1千克原油,1.4286千克标准煤 参考资料: [1] 《测井数据处理与综合解释》(主编 雍世和 张超谟) 石油大学出版社 [2] 《油层物理》 何更生 编 石油工业出版社 [3] 《电缆地层测试器原理及应用》 马建国 付仲金 编著 石油工业出版 [4] …………………… 注: 《测井解释常用公式》是从不同参考资料中摘抄、整理的,以作为测井解释工作的 备忘录。在摘抄和整理过程中难免出现错误,敬请查阅此备忘录的同仁及时指正、修改。 袁福文 2007,08,28 57 58
本文档为【测井解释计算常用公式】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑, 图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
下载需要: 免费 已有0 人下载
最新资料
资料动态
专题动态
is_638173
暂无简介~
格式:doc
大小:636KB
软件:Word
页数:78
分类:企业经营
上传时间:2017-09-07
浏览量:172