一种利用时移地震相关属性求取时移前后含油饱和度变化的方法

技术领域

本发明涉及油藏滚动勘探及开发技术领域，具体涉及一种利用时移地震相关属性求取时移前后含油饱和度变化的方法。

背景技术

研究油藏开发前后含油饱和度变化是滚动勘探和开发中一项重要的研究工作，尤其是在油藏开发过程中，分析油藏开发前后含油饱和度变化是预测油藏剩余油饱和度、提高油藏采收率的重要工作环节，它是评价油藏产能、增储上产、计算石油剩余储量以及调整油田开发方案的重要参数和重要依据，有利于油藏的开发过程的监测，能够确保油藏开发具有最合理的开发方案和开发调整方案，使油藏开发具有最佳的经济效益，对于油田稳产增效有重要意义。

目前储层含油饱和度的变化计算方法主要有岩性分析法、测井法、数值模拟法、油藏工程法、开发地质学法等。以上方法均有其局限性，主要是缺乏井间流体变化的直观反映。如岩心分析法只能是单井点、单油层的含油饱和度变化情况，不能对全油层段、全油田范围进行含油饱和度变化计算；由测井法得到的含油饱和度变化范围要比岩心的范围大，但由于井点反映的是井点处地层状况，因此它仍然不能对井间进行全面预测；数值模拟方法给出的是一个油藏的平均值变化，其预测精度受地层对比解的唯一性、储层物性参数的确定性、非均质描述的客观性、渗透率曲线的确定、历史数据资料的准确性等因素的影响；油藏工程法使用井点含水率或者产液剖面、示踪剂测试等手段推测井间的剩余油分布情况，预测距离短、预测多解性大；开发地质学法是利用井网解释微构造、确定沉积微相，划分流动单元，预测井间波及情况，间接分析井间含油饱和度变化情况，也存在预测距离短、预测多解性大的问题。

中国发明专利CN106019378B公开了一种时移测井曲线的动态重构方法，包括以下步骤：确定目标工区的基础数据；直接测量生产井井筒沿目的层垂直厚度的压力变化；确定生产井在目的层位置的水进方式；确定井筒位置在纵向上的压力变化；对生产井在目的层位置按时间进行流体替换模拟，得到目的层位置在不同时间点的含油饱和度；对模拟的目的层位置在不同时间点的含油饱和度进行筛选，得到重构后的测井曲线的含油饱和度；对岩石物理模型进行标定；计算动态重构后的测井曲线的纵波速度、横波速度和密度；通过正演合成得到合成地震数据；将合成地震数据与实际井旁道进行对比，如果动态重构结果满足相关性要求，方法结束；否则，修改岩石物理模型中影响最大的参数，返回上一步。该专利涉及含油饱和度的计算亦是通过数值模拟方法，同样存在以上所述缺陷。

中国专利申请CN111350497A公开了一种基于含水率的低饱和度油层动态评价方法，包括：1)获取区域试油、投产的储层测井响应特征值；2)计算S

目前，需要一种更为准确计算时移前后含油饱和度变化的方法。

发明内容

为克服以上所述技术缺陷，本发明提供一种利用时移地震相关属性求取时移前后含油饱和度变化的方法，本发明方法可直观反映井间油层变化情况，更为准确地求取井间含油饱和度的变化值。

为实现上述目的，本发明采用以下技术方案：

本发明提供一种利用时移地震相关属性求取时移前后含油饱和度变化的方法，其包括以下步骤：确定时移前后井点含油饱和度变化值；确定时移前后油层相关属性；利用井点含油饱和度变化值与井点处相关属性值交汇求取两者之间的数学关系；利用取得的数学关系将油层相关属性转化为油层全区含油饱和度变化；利用井点时移前后含油饱和度变化值校正转化结果。

进一步地，获取井点含油饱和度变化值与特征井点处相关属性值交汇求取两者之间的数学关系的标准为：在剔除少量异常点后，两者有明显的数学关系，或者，根据井点分布、沉积分布分范围不同，分区域得到两者合理的数学关系。

进一步地，所述井点为油层全区均匀分布的井位。

进一步地，确定油层相关属性时，属性提取层位是油层对应在地震上响应同相轴，属性提取的时窗包含油层附近1--1.5个地震同相轴波长。

进一步地，利用求取的数学关系式将油层时移前后相关属性转化为油层全区含油饱和度变化，同时初步修改不合理的井间含油饱和度变化值。

进一步地，利用井点时移前后含油饱和度变化值校正转化结果时，修正转化后不合理的井间含油饱和度变化值。

与现有技术相比，本发明具有以下优势：

本发明方法利用井点时移前后含油饱和度变化与时移地震相关属性数学关系转化计算时移前后油藏全区含油饱和度变化值，操作过程简单，提高了时移前后含油饱和度变化的计算效率。同时充分发挥地震资料横向分辨率高的优势，时移前后井间含油饱和度变化结果更加准确，有助于提升井间油藏开发动态分析，提高剩余油挖潜效果，对于在油藏中后期开发具有较大的应用价值。

附图说明

构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。

图1为本发明实施例1所述一种利用时移地震相关属性求取时移前后含油饱和度变化的方法。

图2为本发明实施例1中井点处时移前后含油饱和度变化值平面示意图；

图3为本发明实施例1中油层时移前后相关属性提取层位及时窗选取示意图 (左)及相关属性平面图(右)；

图4为本发明实施例1中井点时移前后含油饱和度变化值与井点油层相关属性值交汇求取两者之间数学关系示意图；

图5为利用井点时移前后含油饱和度变化值校正油层时移前后相关属性转化结果并进行修正后的最终结果图。

具体实施方式

应该指出，以下详细说明都是示例性的，旨在对本发明提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。

需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本发明的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作和/或它们的组合。

为了使得本领域技术人员能够更加清楚地了解本发明的技术方案，以下将结合具体的实施例详细说明本发明的技术方案。

实施例1

如图1所示，所述利用时移地震相关属性求取时移前后含油饱和度变化的方法，包括以下步骤：

步骤1、选取油藏平均分布的井点，统计井点处时移前后的含油饱和度资料，用井点处时移前含油饱和度与时移后饱和度求差得到时移前后井点处含油饱和度变化值，如图2所示。

步骤2、解释油层对应的地震反射层位作为提取相关属性的层位，同时选择油层附近同相轴波长1--1.5倍作为相关属性提取的时窗进行相关属性提取，如图3所示。

步骤3、利用步骤2中相关属性提取出井点处时移相关属性值，将其与步骤 1中求取的时移前后井点处含油饱和度变化值进行交汇，分析并剔除异常值，求取两者之间的数学关系，如图4所示。

步骤4、利用步骤3中求取的数学关系式，将油层时移前后相关属性转化为含油饱和度变化值；同时初步分析和修改不合理的井间含油饱和度变化值。

步骤5、利用步骤1中求取的井点时移前后含油饱和度变化值校正步骤4中油层时移前后相关属性转化结果，使得井点处两者吻合，同时再一次修正转化后不合理的井间含油饱和度变化值，得到最终结果，如图5所示。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。