用于振幅随炮检距变化道集分析的油藏模型优化方法
【技术领域】
[0001] 本发明涉及石油的地震勘探领域，尤其涉及一种用于振幅随炮检距变化道集分析 的油藏模型优化方法。
【背景技术】
[0002] 振幅随炮检距变化(AVO)分析是继亮点技术之后又一项利用振幅信息研究岩性、 检测油气的重要方法。在叠前对地震反射振幅随炮检距(入射角）的变化规律进行分析，进 而确定反射界面上覆、下伏介质的岩性特征及物性参数，借此对岩石中孔隙中的流体性质 和岩性做出推断。
[0003] 反射波地震勘探中最重要的假设之一是：两种波阻抗不同的介质之间的界面是个 反射面，当界面两边的波阻抗确定后，反射界面的反射系数也就唯一的确定了。根据反射时 间和反射波振幅，可以反演地下界面的构造形态和介质的某些主要性质。然而实际地震反 射系数是随着炮检距的变化而改变的，AVO技术正是利用这一变化来进行油气流体分析的。 Zo印pritz方程是AVO技术研究的理论基础，但由于该方程比较复杂，一直没有得到直接应 用，后人在Zoeppritz的基础上总结了很多Zoeppritz的简化形式。
[0004] 传统的AVO技术是利用模型正演AVO现象，结合油藏特征，分析不同地质条件下的 油、气、水等不同流体产生的AVO特征，建立相应的AVO标志，在实际的地震记录中，直接识 别岩性及含油气性。
[0005] 然而，一般情况下模型与地下的实际情况误差会比较大，原因是油藏模型通常由 测井资料、地震资料、生产数据以及岩心数据，通过历史拟合构建，通过测井资料仅了解地 下某一点的实际情况，而地震数据虽然在一定程度上解决了对横向的认识，但是纵向上的 分辨率有限，而实际上，地下地质情况是错综复杂的，建模结果往往与实际情况误差较大， 因此，正演的AVO结果对解释地下含油气性可能会带来偏差。
【发明内容】
[0006] 本发明目的是提出一种模型正演结果会更贴近于实际的用于振幅随炮检距变化 道集分析的油藏模型优化方法。
[0007] 本发明通过以下步骤实现：
[0008] 1)根据工区地震解释数据、测井解释数据、生产数据以及岩心数据，通过历史拟 合，建立油藏模型；
[0009] 步骤1)所述的数据包括：流体参数、
[0010] 弹性参数、岩石属性、油藏分区参数、初始化计算参数、输出控制参数、模拟工作的 基本信息以及生产参数；其中：
[0011] 流体参数包括油、气、水的地面密度或重度，油、气的地层体积系数，水的粘度、体 积参数，压缩系数，岩石压缩系数；
[0012] 弹性参数包括横波速度，拉梅常数以及泊松比。
[0013] 步骤I)所述的建立油藏模型，是利用地震解释结果建立格架模型，结合测井数据 以及岩心数据，建立物性模型，在生产数据基础上进行历史拟合得到油藏模型。
[0014] 2)利用工区测井解释数据和油藏模型数据对岩石物理模型标定，建立地震响应参 数与岩石参数之间的非线性关系；
[0015] 所述的岩石物理模型是Voigt与Reuss理论模型或Gassmann理论模型。
[0016] 所述的建立地震响应参数与岩石参数之间的非线性关系是建立声波参数、弹性参 数与岩石物性、流体参数之间的关系。
[0017] 3)根据油藏模型的流体参数及弹性参数，利用标定后的岩石物理模型，计算油藏 模型的纵波速度、横波速度以及密度；再将纵波速度与密度相乘得到合成的阻抗；
[0018] 所述的岩石物理模型是Voigt与Reuss理论模型或Gassmann理论模型。
[0019] 所述计算采用gassmann模型，岩石体积弹性模量K计算：
[0020]
[0021] 其中，Ks为岩石颗粒弹性模量；Kd为干岩石弹性模量；Φ为孔隙度；
[0022] Kf为流体弹性模量，由下式计算：
[0023]
;其中，Kw，K。，Kg分别是水、油、气的弹性模量，S w， S。分别是水、油的饱和度，1-SW-S。为含气饱和度；
[0024] 声波在岩石中的纵波速度Vp由下式计算：
[0025]
;其中，μ为岩石剪切模量；P为岩石密度；Vs为横 波速度。
[0026]
;其中，P w、P。、P g分别是水、 油、气密度，Pma为岩石骨架密度。
[0027] 合成阻抗为：Al= P Vp。
[0028] 4)利用工区测井资料确定油藏模型的上下边界、子波以及正演的地震到的角度范 围和角度间隔，利用纵波速度、横波速度和密度，进行AVO叠前道集正演，生成AVO叠前道集 正演结果；
[0029] 在道集正演前，利用测井资料在无效网格处实测的声波参数和密度值对油藏模型 中无效网格处的纵波速度及密度进行设定；
[0030] 所述测井资料是密度测井曲线以及声波时差测井曲线。
[0031] 所述设定为将测井资料中的密度和速度赋值给对应深度的油藏模型的中的无效 网格处，以避免油藏模型在正演时无效网格中的空值对正演结果的影响。
[0032] 所述子波的相位与频率与实际地震数据的子波相位与频率一致。
[0033] 在道集正演前，根据所述纵波速度、横波速度、密度以及设定的角度范围和角度间 隔计算反射系数；
[0034] 所述的计算反射系数的方法包括boltfeld或Aki。
[0035] 5)将AVO叠前道集正演结果与实际采集的地震数据进行对比，如一致则完成油藏 模型的优化，如不一致修改一致。
[0036] 所述的对比是：
[0037] 某一时间点的井点位置的AVO道集，与时间点的生产测井曲线进行对比，或AVO叠 前道集正演结果与实际地震采集时间点一致的AVO道集与实际采集地震数据的AVO分析进 行对比。
[0038] 所述的一致是：
[0039] AVO特性符合井点处的流体规律，AVO特性一致无区别，振幅随炮检距的增大而增 大或者减小，甚至发生相位的反转或，
[0040] AVO叠前道集正演结果与实际地震采集时间点一致的AVO道集与实际采集地震数 据的AVO分析进行对比，AVO特性是否一致。
[0041] 所述的修改是对油藏模型中的孔隙度或泥质含量或净毛比或含油饱和度或者压 力或者温度进行修改。
[0042] 本发明用于AVO道集分析的油藏模型的优化建立方法，为地质模型的修改提供了 依据，而常规的地质模型的修改通常是根据油藏工程师的经验，这样就有很多的不确定性， 因此本发明可得到与实际地下情况最为接近的油藏模型，有利于油藏开发。
【附图说明】
[0043] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域技术人员来讲，在不付出创造性劳动性的前提下，还可以根 据这些附图获得其他的附图。
[0044] 图1为本发明实施例的用于AVO道集分析的油藏模型优化建立方法的方法流程 图；
[0045] 图2为本发明实施例的用于AVO道集分析的油藏模型优化建立方法的参数计算流 程图；
[0046] 图3为本发明具体实例试验区的模型平面分布图；
[0047] 图4为本发明具体实例试验区的岩石物理模型标定示意图；
[0048] 图5、6、7分别为本发明具体实例实验区油藏模型计算的纵波速度、横波速度、密 度；
[0049