真实的反应了地下油藏在开发过程中的变化，因此，以岩石物理模型为桥梁将地下油藏的 变化转换为声学上的变化，即得到地下油藏在开发中纵波速度、横波速度、密度的变化，将 这种真实的变化再应用到实测的纵波速度、横波速度、密度中，则可以将实测的测井数据回 归到地震数据采集的时间点上，这样在进行测井、地震综合研究时，二者处于同一测量时间 点，为准确预测地下油藏位置及油藏变化建立良好的基础。
【附图说明】
[0044] 图1为基于油藏模型的测井数据历史回归处理方法流程；
[0045] 图2为根据油藏模型和井点解释结果，调整为最佳的岩石物理模型参数；
[0046] 图3为岩石物理模型模板，纵波速度与孔隙度之间的对应关系；
[0047] 图4为在井点位置通过油藏模型合成的纵波速度、横波速度；
[0048] 图5为在井点位置分别计算的在两个地震数据采集时间点，纵波速度、横波速度 的变化率；
[0049] 图6为在井点位置分别回归到与两个地震数据采集时间点的实测的纵波速度、横 波速度；
[0050] 图7为在井点位置在第一个地震数据采集时间点的测井正演结果与实际采集地 震数据对比；
[0051] 图8为在井点位置在第二个地震数据采集时间点的测井正演结果与实际采集地 震数据对比。
【具体实施方式】
[0052] 下面将结合附图，对本发明的技术方案作进一步说明。
[0053] 如图1所示，测井数据历史回归处理方法流程。
[0054] 首先，根据地震解释数据、测井解释数据、生产数据以及岩心数据，通过历史拟合， 建立油藏模型。地质学家对油田进行详细的油藏描述工作，基于地震解释数据，测井解释数 据，岩芯数据以及结合地质家对油田的认识建立三维地质模型。流程进入到步骤2。
[0055] 其二，利用工区测井解释数据和油藏模型数据对岩石物理模型标定，建立地震响 应参数与岩石参数之间的非线性关系。油藏模型数据里没有速度和密度这些值，因此需要 利用测井数据解释结果与油藏模型对应的测井位置的一些动态参数进行拟合，这样一个过 程称之为标定。即：岩石物理模型标定是将通过对各种岩心资料、测井资料和地震资料进 行综合分析，研究岩石岩性、孔隙度、孔隙类型、孔隙流体、流体饱和度和频率等对岩石中弹 性性质如纵波速度和横波速度及衰减的影响，并提出利用地震响应预测岩石物理性质的理 论和方法，是地震响应与储层岩石参数之间联系的桥梁，是进行定量储层预测的基本前提。 图2为标定后岩石物理模型参数，图3为岩石物理模型模板，纵波速度与孔隙度之间的关 系。
[0056] 其三，合成油藏模型的动态纵波速度、横波速度及密度，即根据油藏模型的各个时 间点，合成各个时间的纵波速度、横波速度、密度。图4为在井点位置在某一时间点合成的 纵波速度及横波速度。
[0057] 其四，根据合成的动态纵波速度、横波速度及密度，以地震数据采集时间点作为基 准，分别计算纵波速度、横波速度及密度的变化率，图5为测井数据采集时间点相对地震数 据采集时间点的纵波速度、横波速度、密度的变化率，可以看到在本实例中有两次采集的地 震数据分别为1991年和2010年。
[0058] 其五，利用实测的纵波速度、横波速度、密度以及步骤四中计算的纵波速度、横波 速度、密度的变化率，计算对应地震数据采集时间点的实测纵波速度、横波速度、密度。图6 为本实例中选取的井计算的经过回归后的纵波速度、横波速度、密度。
[0059] 其六，利用回归后的纵波速度、横波速度、密度进行正演，得到的合成记录与采集 的地震记录进行一致性对比分析，以判别回归的测井数据是否正确，图7和图8分别为向两 个时间点（1991年和2010年）回归的测井数据正演的结果与采集的地震记录的对比分析。
[0060] 下面通过典型实施例对本发明作更加细化的说明。
[0061] 典型本实施例的包括下列步骤：
[0062] 1)根据胜利油田单56工区地震解释数据、测井解释数据、生产数据以及岩心数 据，通过历史拟合，建立油藏模型；
[0063] 所述地震解释数据包括：由地震数据解释的层位，断层数据；
[0064] 所述测井解释数据包括：地质分层、单井相、砂泥岩分布、渗透率、含油饱和度、含 水饱和度等；
[0065] 所述生产数据包括：实时的流体的流量、含水饱和度、温度、压力等；
[0066] 所述岩心数据是指钻井时从地下取出的岩心，在实验室进行测量，主要包括：渗透 率、孔隙度、电阻率等；
[0067] 所述历史拟合是指由油藏模型被建立起来以后，是否完全反映油气藏实际，并未 经过检验。只有利用将生产和注入的历史数据输入模型并运行模拟器，再将计算结果与油 气藏实际动态相比，才能确定模型中采用的油气藏描述是否有效。若计算获得的动态数据 与油藏实际动态数据差别甚远，我们就必须不断调整输入模型基本数据，直到由模拟器计 算得到的动态与油藏生产的实际动态达到满意的拟合为止。
[0068] 所述的建立油藏模型：第一步是利用地震解释的层位和断层数据，通过空间上的 网格剖分，建立空间上的格架模型，并将油藏剖分为一系列独立的单元体；第二步是结合测 井数据以及岩心数据，通过空间上的插值，在空间上的每一个单元中赋值，赋予不同的物理 属性；第三步：通过模拟器对历史的油气藏动态数据进行模拟，并且与实际的井点上的油 气藏动态数据进行对比，以确定模拟的油气藏动态数据是否准确。如果模拟的动态数据与 油藏实际动态数据差别较大，则调整输入的模型参数，重新进行拟合，最终得到与实际油气 藏动态数据一致的油藏模型。
[0069] 2)利用工区测井解释数据和油藏模型数据对岩石物理模型标定，建立地震响应参 数与岩石参数之间的非线性关系；
[0070] 所述的岩石物理模型是Gassmann理论模型。
[0071] 所述标定是指确定岩石物理模型中的各项参数。
[0072] 所述的建立地震响应参数与岩石参数之间的非线性关系是建立声波参数、弹性参 数与岩石物性、流体参数之间的关系。
[0073] 3)根据油藏模型的流体参数及弹性参数，利用标定后的岩石物理模型，计算油藏 模型的纵波速度、横波速度以及密度；再将纵波速度与密度相乘得到合成的阻抗；
[0074] 所述的岩石物理模型是Gassmann理论模型。
[0075] 所述计算采用gassmann模型，岩石体积弹性模量K计算：
[0076]
[0077] 其中，Ks为岩石颗粒弹性模量；Kd为干岩石弹性模量；〇为孔隙度；
[0078] Kf为流体弹性模量，由下式计算：
[0079]
;其中，Kw，K。，Kg分别是水、油、气的弹性模量，S w， S。分别是水、油的饱和度，1-SW-S。为含气饱和度；
[0080] 声波在岩石中的纵波速度Vp由下式计算：
[0081]
〔中，y为岩石剪切模量；P为岩石密度；Vs为横 波速度^
[0082]
其中，P w、P。、P g分别是水、 油、气密度，Pm为岩石骨架密度。
[0083] 4)利用步骤3)中计算的各个时间点的纵波速度、横波速度、密度，与地震数据采 集时间点时的纵波速度、横波速度、密度，计算油藏模型在各个开发时间点的变化率；
[0084] 所述油藏模型的各个开发时间点是在建立油藏模型时