测井曲线预测方法研究.docVIP

测井曲线预测方法研究 砂岩层岩石物理计算 对于泥砂岩薄互层储层，钻井取心表明其通常是泥岩和砂岩的互层，油藏开发过程中油藏参数变化主要发生在砂岩薄层．因此，砂岩层岩石物理计算是进行泥砂岩储层时移测井曲线预测的基础，而泥砂岩薄互层中砂岩层含量计算将在后面进行详细论述．油藏开发过程中砂岩薄层中饱和度、压力场和温度场都可能发生变化．而岩石物理实验和实际测井数据分析表明，除注蒸汽或火烧稠油开采外，温度场变化对岩石骨架弹性参数影响较小，主要表现为对流体弹性参数影响，而压力场的变化对岩石骨架和流体弹性参数都会产生影响．因此，时移测井曲线预测中重点考虑温度、压力以及矿化度对流体性质的影响和流体饱和度、有效压力变化对饱和流体储层弹性参数的影响． 1.孔隙流体弹性参数计算 孔隙流体在很大程度上影响着岩石的弹性参数．孔隙流体包括气体、原油和地层水，在成分和物理性质上差别很大，并组成了一个动态系统．在此系统内，流体的成分和物理相态都随压力和温度而变化．对于气体组分，其密度和体积模量变化与温度、压力有密切关系，已有大量的数据模型对其变化规律进行了描述［6］．对于原油组分，它是极其复杂的生物化合物的混合物，其密度rho;和纵波速度Voil变化与参考密度rho;0、压力及温度密切相关，其计算方程为［6］rho;=rho;p0．972+3．81times;10－4(t+17．78)1．175;(1)Voil=2096rho;02．6－rho;()01/2－3．7t+4．64p+0．0115［4．12(1．08rho;0－1－1)1/2－1］tp．(2)式中:p为压力，MPa;t为油藏温度;℃;rho;p为原油在某压力条件下的密度，即rho;p=rho;0+(0．00277p－1．71times;10－7p3)times;(rho;0－1．15)2+3．49times;10－4p．(3)对于储层中地层水，其弹性参数受油藏温度、压力和矿化度影响．不同温度、压力和矿化度条件下地层水密度和速度计算方程为rho;B=rho;W+S{0．668+0．44S+1．0times;10－6［300p－2400pS+t(80+3t－3300S－13p+47pS)］}(4)VB=VW+S(1170－9．6t+0．055t2－8．5times;10－5t3+2．6p－0．0029tp－0．0476p2)+S3/2(780－10p+0．16p2)－820S2．(5)式中:S为盐的质量浓度，kg/L;VW为在100℃和100MPa条件下测得的纯水速度;rho;W和rho;B是纯水和盐水的密度．方程中速度和密度单位分别为m/s和g/cm3．当油藏孔隙流体为2种或多种组分的混合体时，混合流体的体积模量Kf和密度rho;f利用如下的Wood方程和物质平衡方程计算，即1Kf=sum;Ni=1ciKi．(6)rho;f=sum;Ni=1cirho;i．(7)式中Ki，ci和rho;i分别表示单一流体组分的体积模量、体积分数和密度［7］． 2.饱和流体储层弹性参数计算 油藏开发前后其孔隙度和矿物组分通常不发生变化，因此储层弹性参数主要受孔隙流体弹性参数与压力变化影响，而流体弹性参数变化对储层弹性参数影响可以利用Gassmann方程进行计算［8］，即Ksat=K?Keff－(1+)KfKeff/K+Kf(1－)Kf+K－KfKeff/K;Gsat=Geff．(8)式中:Keff，Gsat分别为饱和流体岩石的体积模量和剪切模量;为岩石孔隙度;Kf为孔隙流体体积模量，其值可通过式(6)计算;Keff，Geff分别为干岩石体积模量和剪切模量，其值可以通过初始饱和状态岩石弹性模量计算;K为岩石骨架体积模量，可以利用Voigt-Reuss-Hill方程［7］计算，即M=12sum;mi=1IiMi+sum;mi=1IiM()i－[]1;(9)G=12sum;mi=1IiGi+sum;mi=1IiG()i－[]1;(10)K=M－43G．(11)式中:m为组成岩石骨架的矿物种类总数;Ii，Mi和Gi分别为组成岩石骨架第i种矿物体积分数、纵波模量和剪切模量;岩石和流体弹性模量单位为GPa．对于油藏压力变化影响，Eberhart-Phillips和Han等人(1989年)基于饱和水泥砂岩速度实验室测量结果，分析了有效压力变化对纵、横波速度影响规律［9］．Shapiro(2003年)基于这一研究建立了广泛适用于饱和水砂岩油藏的压力变化对纵横波速度影响计算方程［10-11］V=V0－a?eff－bV槡sh_sand+cpeffp0－edp()eff．(12)式中:eff和Vsh_sand分别为泥砂薄互层中砂岩薄层的有效孔隙度和泥质百分含量;peff为油藏有效压力，即油藏围压与孔隙流体压力的差