[发明专利]一种利用佐普立兹阻抗进行储层预测的方法

摘要

一种基于佐普立兹(Zoeppritz)阻抗进行储层预测的方法，将佐普立兹佐普立兹与其它类型的弹性阻抗对比，佐普立兹反射系数的精确公式，计算沿射线或共角度的佐普立兹(Zoeppritz)弹性阻抗；对井资料进行精细解释与方波化，以便消除利用井资料和计算佐普立兹(Zoeppritz)弹性阻抗的累积误差；利用精确反射系数计算人工合成地震记录、标定层位；对于地震资料，采用小角度范围或者射线范围内叠加的方式，提高地震资料的信噪比。本发明克服了基于Aki和Richards近似公式的弹性阻抗中的诸多假设，提高了弹性阻抗正、反演的精度。也解释清楚了地震反演中低频分量与高频之间的匹配关系。

主权项

一种基于佐普立兹(Zoeppritz)阻抗进行储层预测的方法，其特征在于，包括以下步骤：1)将佐普立兹弹性阻抗与其它类型的弹性阻抗对比：采用选择初始佐普立兹(Zoeppritz)弹性阻抗为其它需要对比的波阻抗、弹性阻抗的初值，选择为1或任意常数；其中佐普立兹反射系数的精确公式表示为： R i ( θ ) = ZEI ( α i + 1 , β i + 1 , ρ i + 1 , P ) - ZEI ( α i , β i , ρ i , P ) ZEI ( α i + 1 , β i + 1 , ρ i + 1 , P ) + ZEI ( α i , β i , ρ i , P ) - - - ( 1 ) 其反演公式为： ZEI ( α n , β n , ρ n , p n ) = ZEI ( α 0 , β 0 , ρ 0 , p 0 ) Π i = 1 N 1 + R p , i ( θ ) 1 - R p , i ( θ ) - - - ( 2 ) ZEI即为佐普立兹(Zoeppritz)弹性阻抗；公式中α、β、ρ分别为纵波速度、横波速度、岩层密度；i为岩层序号；P为波传播的射线参数；α0，β0，ρ0，p0分别为初始岩层的纵波速度、横波速度、岩层密度和波传播的射线参数；αn，βn，ρn，pn分别为第n个岩层的纵波速度、横波速度、岩层密度和波传播的射线参数；θ为波传播的入射角度；Rp，i(θ)为第i个岩层的波传播的射线参数为p的佐普立兹反射系数；N为大于等于1的整数；在实际井资料计算佐普立兹(Zoeppritz)弹性阻抗时，计算沿射线或共角度的佐普立兹(Zoeppritz)弹性阻抗；2)对井资料进行精细解释与方波化，以便消除利用井资料和公式(2)计算佐普立兹(Zoeppritz)弹性阻抗的累积误差；3)利用低通滤波器获得佐普立兹(Zoeppritz)弹性阻抗的低频分量和高频部分；低频部分作为初始模型进行迭代反演，或累加到反演出来的高频成 分中；高频部分用于约束地震资料反演时的高频成分；4)利用精确反射系数计算人工合成地震记录、标定层位；5)对于地震资料，采用小角度范围或者射线范围内叠加的方式，提高地震资料的信噪比。