阅读说明：本技术 一种利用三维地震恢复相对微古地貌的工业化流程 (Industrialized process for recovering relatively micro ancient landform by using three-dimensional earthquake ) 是由 张宪国 张涛 吴啸啸 刘杰 刘华峰 黄德榕 刘玉从 段冬平 黄鑫 于 2020-06-08 设计创作，主要内容包括：本发明公开了一种利用三维地震恢复相对微古地貌的工业化流程,具体方法是,利用地层顶界面和底界面的地震解释求取地层厚度,利用井点的测井信息对地震解释得到的地层厚度进行井点校正,利用地震倾角属性分析得到地层倾角余弦值,利用地层顶面和底面的倾角余弦平均值与井点校正后的地层厚度求取地层真厚度分布,上述流程得到的地层真厚度可以反映目的层沉积时期的相对微古地貌形态。本发明为少井地区基于三维地震资料恢复微古地貌提供了一种简单、合理而有效的方法。(The invention discloses an industrialized process for recovering relatively micro-ancient landforms by using a three-dimensional earthquake, which comprises the steps of solving the thickness of a stratum by using earthquake interpretation of a top interface and a bottom interface of the stratum, carrying out well point correction on the thickness of the stratum obtained by the earthquake interpretation by using well logging information of well points, obtaining a cosine value of a dip angle of the stratum by using earthquake dip angle attribute analysis, solving the true thickness distribution of the stratum by using the cosine average value of the dip angle of the top surface and the bottom surface of the stratum and the thickness of the stratum obtained by the well point correction, wherein the true thickness of the stratum obtained by the process can reflect the relatively micro-ancient landform form of a target stratum in a deposition period. The invention provides a simple, reasonable and effective method for recovering the micro ancient landform based on the three-dimensional seismic data in the area with few wells.)

一种利用三维地震恢复相对微古地貌的工业化流程

技术领域

本发明涉及油气藏勘探开发技术领域，更具体的说是涉及一种利用三维地震资料恢复相对微古地貌的工业化流程。

背景技术

古地貌是控制油气储层沉积的重要因素，能够控制储层发育类型、展布范围、内部结构和沉积相模式。古地貌研究是储层沉积相研究中的重要内容。长期以来，国内外学者对古地貌的恢复方法持续研究，形成了一系列的古地貌恢复方法。

在区域古地貌恢复中，主要考虑的因素包括构造形态的变化和地层厚度的变化，古地貌的恢复主要围绕这两方面因素的恢复开展。对于构造形态变化的恢复，一方面是构造沉降量的恢复，另一方面是水平走滑量的恢复。对于地层厚度变化的恢复，主要是从三个方面开展工作，一是去压实校正，解决差异压实引起的地层厚度变化，二是通过古水深的恢复，解决差异沉积作用引起的地层厚度变化，三是进行剥蚀量的恢复，解决剥蚀作用差异形成的地层厚度变化。

针对这些古地貌恢复问题，目前常用的方法有四大类，一是残留厚度和补偿厚度印模法，该方法的思路是恢复地层剥蚀量及原始厚度，利用古厚度反映古地貌形态；二是地层回剥法，该方法通过定量盆地沉降量将盆地发育某一时期原始地貌形态进行恢复的方法；三是沉积学分析法，该方法结合成因相分析、古构造发育特点、古流向分析等多种沉积学分析手段进行古地貌恢复；四是层序地层学恢复法，该方法选取等时面为参考面恢复出下伏地层沉积前的原始古地貌形态。

对于小尺度的微古地貌恢复，目前常用的方法有三大类，第一大类是古沟谷和同沉积断层的识别；第二大类是沉积微相反推法；第三大类是利用密井网资料进行砂体微构造的表征，这类方法主要适合于井点密集的地区。对于海上油田以及陆上油田勘探阶段等井点少的地区，利用地震资料进行微古地貌恢复缺乏有效的方法和技术。对于少井区依靠地震资料开展微古地貌恢复现有方法存在的缺陷，本发明研发了利用地震资料恢复微古地貌的工业化流程。

发明内容

本发明提供了一种利用三维地震资料恢复相对微古地貌的工业化流程，具体的技术方案包括如下步骤：

步骤(1)：在三维地震资料上，利用井震结合的地震资料解释方法，分别获得目的地层顶面深度的数据网格T_top和底面深度的数据网格T_bot；

步骤(2)：目的地层底面与顶面相减得到目的地层的厚度数据网格T_h；

步骤(2)所述的计算公式为T_h＝T_bot-T_top；

步骤(3)：选取研究区任意一口井A，在地层厚度网格T_h上读取A井处地震解释得到的地层厚度h_t，利用A井的测井曲线数据获得A井处的地层厚度h_l，两者相减得到A井点处地震解释结果的误差D_A；

步骤(3)所述的计算公式为D_A＝h_l-h_t；

步骤(4)：利用步骤(3)的方法得到研究区每一口井点处的地层厚度误差，并对这些数据点进行插值和外推计算，得到覆盖研究区的误差网格T_d；

步骤(5)：将步骤(2)的T_h与步骤(4)的T_d相加，得到井震校正后的地层厚度数据网格T_hc；

步骤(5)所述的计算公式为T_hc＝T_h+T_d；

步骤(6)：利用三维地震属性分析计算目的层顶面的地震倾角数据网格T_dat，对该网格上每个点求取余弦值，获得目的层顶面的倾角余弦值网格T_cost；

步骤(7)：利用三维地震属性分析计算目的层底面的地震倾角数据网格T_dab，对该网格上每个点求取余弦值，获得目的层底面的倾角余弦值网格T_cosb；

步骤(8)：将步骤(6)和步骤(7)分别求取的目的层顶面和底面地层倾角余弦值网格进行算数平均，获得研究区目的层的平均地层倾角余弦值网格T_cosave；

步骤(8)所述的计算公式为T_cosave＝(T_cost+T_cosb)/2；

步骤(9)：将步骤(5)获得的井震校正后的地层厚度数据网格与步骤(8)获得的目的层的平均地层倾角余弦值网格相乘，得到目的地层真厚度数据网格T_hr；

步骤(9)所述的计算公式为T_hr＝T_hc×T_cosave；

步骤(10)：利用步骤(9)得到的目的地层真厚度数据网格T_hr制作等值线图，作为目的层沉积时期的相对古地貌形态。

具体实施方式

按照上述发明内容中的步骤实施，具体包括以下步骤：

(1)在三维地震资料上，利用井震结合的地震资料解释方法，分别获得目的地层顶面深度的数据网格T_top和底面深度的数据网格T_bot；

(2)目的地层底面与顶面相减得到目的地层的厚度数据网格T_h，其中T_h＝T_bot-T_top；

(3)选取研究区任意一口井A，在地层厚度网格T_h上读取A井处地震解释得到的地层厚度h_t，利用A井的测井曲线数据获得A井处的地层厚度h_l，两者相减得到A井点处地震解释结果的误差D_A，其中D_A＝h_l-h_t；

(4)研究区每一口井点处的地层厚度误差进行插值和外推计算，得到覆盖研究区的误差网格T_d；

(5)将网格T_h与网格T_d相加，得到井震校正后的地层厚度数据网格T_hc，其中T_hc＝T_h+T_d；

(6)利用三维地震属性分析计算目的层顶面的地震倾角数据网格T_dat，对该网格上每个点求取余弦值，获得目的层顶面的倾角余弦值网格T_cost；

(7)利用三维地震属性分析计算目的层底面的地震倾角数据网格T_dab，对该网格上每个点求取余弦值，获得目的层底面的倾角余弦值网格T_cosb；

(8)计算研究区目的层的平均地层倾角余弦值网格T_cosave,其中T_cosave＝(T_cost+T_cosb)/2；

(9)将网格T_hc与网格T_cosave相乘，得到目的地层真厚度数据网格T_hr，其中T_hr＝T_hc×T_cosave；

(10)利用目的地层真厚度数据网格T_hr制作等值线图，作为目的层沉积时期的相对微古地貌形态。

实施例

本实施例的研究区位于我国南海，面积为100平方公里的油田开发区块，该区目的层段为古近系流沙港组一段(以下简称流一段)地层，发育扇三角洲沉积，研究区共有15口井，采用本发明的方法恢复流一段的微古地貌，为油田开发阶段的沉积相研究提供支持。

(1)通过井震结合的三维地震资料解释，获得目的地层顶面深度的数据网格T_top和底面深度的数据网格T_bot；

(2)目的地层底面与顶面相减得到目的地层的厚度数据网格T_h；

(3)选取研究区井A₁，在地层厚度网格T_h上读取A₁井处地震解释得到的地层厚度为25m，利用A₁井的测井曲线数据获得A₁井处的地层厚度为25.3m，两者相减得到A井点处地震解释结果的误差为0.3m；

(4)按照步骤(3)的方法对研究区15口井逐一求取井点处的厚度误差，将15口井的地层厚度误差进行插值和外推计算，得到覆盖研究区的误差网格T_d；

(5)将网格T_h与网格T_d相加，得到井震校正后的地层厚度数据网格T_hc，其中T_hc＝T_h+T_d；

(6)对流一段顶面进行地震倾角属性分析，得到流一段顶面的地震倾角数据网格T_dat，对该网格上每个点求取余弦值，获得目的层顶面的倾角余弦值网格T_cost；

(7)对流一段底面进行地震倾角属性分析，得到流一段底面的地震倾角数据网格T_dab，对该网格上每个点求取余弦值，获得目的层底面的倾角余弦值网格T_cosb；

(8)计算研究区流一段地层的平均地层倾角余弦值网格T_cosave,其中T_cosave＝(T_cost+T_cosb)/2；

(9)将网格T_hc与网格T_cosave相乘，得到流一段地层真厚度数据网格T_hr，其中T_hr＝T_hc×T_cosave；

(10)利用流一段地层真厚度数据网格T_hr制作等值线图，作为目的层沉积时期的相对微古地貌形态。

对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。