一种基于油藏模型的测井数据历史回归处理方法
【技术领域】
[0001] 本发明属于油田开发的测井、地震、油藏综合数据处理技术领域，特别是涉及到一 种基于油藏模型的测井数据历史回归处理方法。
【背景技术】
[0002] 石油测井是石油勘探、开发的"眼睛"。在油气田勘探、开发的不同阶段，石油测井 的目的和任务是不同的。一般来说，裸眼井测井（下套管之前的井称裸眼井，因此下套管之 前进行的测井称为裸眼井测井）的主要目的和任务是发现和评价油气层的储集性能及生 产能力；而生产测井（油水井投入生产以后进行的测井称为生产测井）的主要目的是见识 和分析油气层的开发动态及生产情况。开发测井是指在油气田整个开发期间进行的所有测 井项目。开发测井的主要对象为裸眼完成的生产井和下套管的生产井，用于分析目前的生 产动态及井内技术状况。
[0003] 测井数据详细的记录了地层岩性、物性在纵向上的细微变化，具有较高的纵向分 辨率，在井位置处，能准确的计算储层参数，正确评价储层。地震资料连续地记录了地层在 横向上的变化，具有较好的横向连续性。因此，在勘探和开发阶段往往是结合测井和地震数 据进行综合研究。通常测井数据的采集时间点与地震数据的时间点很少一致，因此，当地 下储层物性变化较大时，如果直接拿来进行综合研究，则会因采集时间点的不一致给综合 研究带来很大误差，因此，需要对测井数据进行历史回归处理，将实际测井数据校正到地震 数据的采集时间点，为更好、更准确的研究储层建立基础。
【发明内容】
[0004] 本发明的目的是针对现有测井技术在开发应用中存在的问题，提出基于油藏模型 的测井数据历史回归处理方法，即测井与地震进行综合研究时，二者的采集时间点往往不 一致所带来的系统误差。
[0005] 本发明总体技术方案是
[0006] -种基于油藏模型的测井数据历史回归方法，包括：
[0007] 1)根据地震解释数据、测井解释数据、生产数据以及岩心数据，通过历史拟合，建 立油藏模型；
[0008] 2)利用工区测井解释数据和油藏模型数据对岩石物理模型标定，建立地震响应参 数与岩石参数之间的非线性关系；
[0009] 3)合成油藏模型的动态纵波速度、横波速度及密度；
[0010] 4)根据合成的动态纵波速度、横波速度及密度，以地震数据采集时间点作为基准， 分别计算纵波速度、横波速度及密度的变化率；
[0011] 5)利用实测的纵波速度、横波速度、密度以及步骤四中计算的纵波速度、横波速 度、密度的变化率，计算对应地震数据采集时间点的实测纵波速度、横波速度、密度；
[0012] 6)利用回归后的纵波速度、横波速度、密度进行正演，得到的合成记录与采集的地 震记录进行一致性对比分析，以判别回归的测井数据是否正确。
[0013] 上述方案的细化方案包括：
[0014] 1)根据工区地震解释数据、测井解释数据、生产数据以及岩心数据，通过历史拟 合，建立油藏模型；
[0015] 所述地震解释数据包括：由地震数据解释的层位，断层数据；
[0016] 所述测井解释数据包括：地质分层、单井相、砂泥岩分布、渗透率、含油饱和度、含 水饱和度；
[0017] 所述生产数据包括：实时的流体的流量、含水饱和度、温度、压力；
[0018] 所述岩心数据是指钻井时从地下取出的岩心，在实验室进行测量，主要包括：渗透 率、孔隙度、电阻率；
[0019] 所述的建立油藏模型：第一步是利用地震解释的层位和断层数据，通过空间上的 网格剖分，建立空间上的格架模型，并将油藏剖分为一系列独立的单元体；第二步是结合测 井数据以及岩心数据，通过空间上的插值，在空间上的每一个单元中赋值，赋予不同的物理 属性；第三步：通过模拟器对历史的油气藏动态数据进行模拟，并且与实际的井点上的油 气藏动态数据进行对比，以确定模拟的油气藏动态数据是否准确，如果模拟的动态数据与 油藏实际动态数据差别较大，则调整输入的模型参数，重新进行拟合，最终得到与实际油气 藏动态数据一致的油藏模型；
[0020] 2)利用工区测井解释数据和油藏模型数据对岩石物理模型标定，建立地震响应参 数与岩石参数之间的非线性关系；
[0021] 所述的岩石物理模型是Gassmann理论模型；
[0022] 所述标定是指确定岩石物理模型中的各项参数；
[0023] 所述的建立地震响应参数与岩石参数之间的非线性关系是建立声波参数、弹性参 数与岩石物性、流体参数之间的关系；
[0024] 3)根据油藏模型的流体参数及弹性参数，利用标定后的岩石物理模型，计算油藏 模型的纵波速度、横波速度以及密度；再将纵波速度与密度相乘得到合成的阻抗；
[0025] 所述计算采用gassmann模型，岩石体积弹性模量K计算：
[0026]
[0027] 其中，Ks为岩石颗粒弹性模量；Kd为干岩石弹性模量；〇为孔隙度；
[0028] Kf为流体弹性模量，由下式计算：
[0029]
;其中，Kw，K。，Kg分别是水、油、气的弹性模量，S w， S。分别是水、油的饱和度，1-SW-S。为含气饱和度；
[0030] 声波在岩石中的纵波速度Vp由下式计算：
[0031]
其中，y为岩石剪切模量；P为岩石密度；Vs为横 波速度。
[0032]
?，其中，Pw、P。、Pg*别是水、 油、气密度，Pm为岩石骨架密度。
[0033] 4)利用步骤3)中计算的各个时间点的纵波速度、横波速度、密度，与地震数据采 集时间点时的纵波速度、横波速度、密度，计算油藏模型在各个开发时间点的变化率；
[0034] 所述油藏模型的各个开发时间点是在建立油藏模型时，人为给定的具有相同时间 间隔的各个时间状态；
[0035] 所述纵波为传播方向与振动防线一致的波，横波为振动方向与传播方向垂直的 波；
[0036] 5)对比实测测井数据时间点的纵波速度、横波速度，密度与油藏模型对应时间点 的纵波速度、横波速度，密度，判断二者是否一致，如不一致，则返回步骤2)，调整岩石物理 模型参数，如一致，则可进入步骤6);
[0037] 所述实测的测井数据是指利用测井仪器在井筒中实际测量的曲线数据，包括声波 时差（DT)、横波声波时差（DTS)、密度（RH0B)，其中，纵波速度为DT的倒数乘以106,横波速 度为DTS的倒数乘以10 6 ;
[0038] 所述二者是否一致是指同一时间状态标准下合成纵波速度、横波速度、密度与实 际测量的纵波速度、横波速度、密度在趋势以及数值范围上是否一致；
[0039] 6)结合实测测井数据与油藏模型计算的各个开发时间段的变化率，将实测测井数 据回归到地震数据采集时间点；
[0040] 所述变化率是指以地震数据采集的时间点作为基础，计算别的时间点的油藏动态 变化百分比；
[0041] 所述回归是指以计算的变化率，对实测测井数据进行基于时间点的校正。
[0042] 其中：步骤2)中标定岩石物理模型中的各项参数包括岩石岩性、孔隙度、孔隙类 型、孔隙流体、流体饱和度和频率，以及反应岩石中弹性性质的包括纵波速度和横波速度及 衰减的参数。
[0043] 本发明的技术效果是：由于油藏模型数据是一个随时间变化的动态数据，其客观、