本发明涉及油气藏勘探开发技术领域，更具体的说是涉及一种利用三维地震资料恢复相对微古地貌的工业化流程。
背景技术：
古地貌是控制油气储层沉积的重要因素，能够控制储层发育类型、展布范围、内部结构和沉积相模式。古地貌研究是储层沉积相研究中的重要内容。长期以来，国内外学者对古地貌的恢复方法持续研究，形成了一系列的古地貌恢复方法。
在区域古地貌恢复中，主要考虑的因素包括构造形态的变化和地层厚度的变化，古地貌的恢复主要围绕这两方面因素的恢复开展。对于构造形态变化的恢复，一方面是构造沉降量的恢复，另一方面是水平走滑量的恢复。对于地层厚度变化的恢复，主要是从三个方面开展工作，一是去压实校正，解决差异压实引起的地层厚度变化，二是通过古水深的恢复，解决差异沉积作用引起的地层厚度变化，三是进行剥蚀量的恢复，解决剥蚀作用差异形成的地层厚度变化。
针对这些古地貌恢复问题，目前常用的方法有四大类，一是残留厚度和补偿厚度印模法，该方法的思路是恢复地层剥蚀量及原始厚度，利用古厚度反映古地貌形态；二是地层回剥法，该方法通过定量盆地沉降量将盆地发育某一时期原始地貌形态进行恢复的方法；三是沉积学分析法，该方法结合成因相分析、古构造发育特点、古流向分析等多种沉积学分析手段进行古地貌恢复；四是层序地层学恢复法，该方法选取等时面为参考面恢复出下伏地层沉积前的原始古地貌形态。
对于小尺度的微古地貌恢复，目前常用的方法有三大类，第一大类是古沟谷和同沉积断层的识别；第二大类是沉积微相反推法；第三大类是利用密井网资料进行砂体微构造的表征，这类方法主要适合于井点密集的地区。对于海上油田以及陆上油田勘探阶段等井点少的地区，利用地震资料进行微古地貌恢复缺乏有效的方法和技术。对于少井区依靠地震资料开展微古地貌恢复现有方法存在的缺陷，本发明研发了利用地震资料恢复微古地貌的工业化流程。
技术实现要素：
本发明提供了一种利用三维地震资料恢复相对微古地貌的工业化流程，具体的技术方案包括如下步骤：
步骤(1)：在三维地震资料上，利用井震结合的地震资料解释方法，分别获得目的地层顶面深度的数据网格ttop和底面深度的数据网格tbot；
步骤(2)：目的地层底面与顶面相减得到目的地层的厚度数据网格th；
步骤(2)所述的计算公式为th＝tbot-ttop；
步骤(3)：选取研究区任意一口井a，在地层厚度网格th上读取a井处地震解释得到的地层厚度ht，利用a井的测井曲线数据获得a井处的地层厚度hl，两者相减得到a井点处地震解释结果的误差da；
步骤(3)所述的计算公式为da＝hl-ht；
步骤(4)：利用步骤(3)的方法得到研究区每一口井点处的地层厚度误差，并对这些数据点进行插值和外推计算，得到覆盖研究区的误差网格td；
步骤(5)：将步骤(2)的th与步骤(4)的td相加，得到井震校正后的地层厚度数据网格thc；
步骤(5)所述的计算公式为thc＝th+td；
步骤(6)：利用三维地震属性分析计算目的层顶面的地震倾角数据网格tdat，对该网格上每个点求取余弦值，获得目的层顶面的倾角余弦值网格tcost；
步骤(7)：利用三维地震属性分析计算目的层底面的地震倾角数据网格tdab，对该网格上每个点求取余弦值，获得目的层底面的倾角余弦值网格tcosb；
步骤(8)：将步骤(6)和步骤(7)分别求取的目的层顶面和底面地层倾角余弦值网格进行算数平均，获得研究区目的层的平均地层倾角余弦值网格tcosave；
步骤(8)所述的计算公式为tcosave＝(tcost+tcosb)/2；
步骤(9)：将步骤(5)获得的井震校正后的地层厚度数据网格与步骤(8)获得的目的层的平均地层倾角余弦值网格相乘，得到目的地层真厚度数据网格thr；
步骤(9)所述的计算公式为thr＝thc×tcosave；
步骤(10)：利用步骤(9)得到的目的地层真厚度数据网格thr制作等值线图，作为目的层沉积时期的相对古地貌形态。
具体实施方式
按照上述发明内容中的步骤实施，具体包括以下步骤：
(1)在三维地震资料上，利用井震结合的地震资料解释方法，分别获得目的地层顶面深度的数据网格ttop和底面深度的数据网格tbot；
(2)目的地层底面与顶面相减得到目的地层的厚度数据网格th，其中th＝tbot-ttop；
(3)选取研究区任意一口井a，在地层厚度网格th上读取a井处地震解释得到的地层厚度ht，利用a井的测井曲线数据获得a井处的地层厚度hl，两者相减得到a井点处地震解释结果的误差da，其中da＝hl-ht；
(4)研究区每一口井点处的地层厚度误差进行插值和外推计算，得到覆盖研究区的误差网格td；
(5)将网格th与网格td相加，得到井震校正后的地层厚度数据网格thc，其中thc＝th+td；
(6)利用三维地震属性分析计算目的层顶面的地震倾角数据网格tdat，对该网格上每个点求取余弦值，获得目的层顶面的倾角余弦值网格tcost；
(7)利用三维地震属性分析计算目的层底面的地震倾角数据网格tdab，对该网格上每个点求取余弦值，获得目的层底面的倾角余弦值网格tcosb；
(8)计算研究区目的层的平均地层倾角余弦值网格tcosave,其中tcosave＝(tcost+tcosb)/2；
(9)将网格thc与网格tcosave相乘，得到目的地层真厚度数据网格thr，其中thr＝thc×tcosave；
(10)利用目的地层真厚度数据网格thr制作等值线图，作为目的层沉积时期的相对微古地貌形态。
实施例
本实施例的研究区位于我国南海，面积为100平方公里的油田开发区块，该区目的层段为古近系流沙港组一段(以下简称流一段)地层，发育扇三角洲沉积，研究区共有15口井，采用本发明的方法恢复流一段的微古地貌，为油田开发阶段的沉积相研究提供支持。
(1)通过井震结合的三维地震资料解释，获得目的地层顶面深度的数据网格ttop和底面深度的数据网格tbot；
(2)目的地层底面与顶面相减得到目的地层的厚度数据网格th；
(3)选取研究区井a1，在地层厚度网格th上读取a1井处地震解释得到的地层厚度为25m，利用a1井的测井曲线数据获得a1井处的地层厚度为25.3m，两者相减得到a井点处地震解释结果的误差为0.3m；
(4)按照步骤(3)的方法对研究区15口井逐一求取井点处的厚度误差，将15口井的地层厚度误差进行插值和外推计算，得到覆盖研究区的误差网格td；
(5)将网格th与网格td相加，得到井震校正后的地层厚度数据网格thc，其中thc＝th+td；
(6)对流一段顶面进行地震倾角属性分析，得到流一段顶面的地震倾角数据网格tdat，对该网格上每个点求取余弦值，获得目的层顶面的倾角余弦值网格tcost；
(7)对流一段底面进行地震倾角属性分析，得到流一段底面的地震倾角数据网格tdab，对该网格上每个点求取余弦值，获得目的层底面的倾角余弦值网格tcosb；
(8)计算研究区流一段地层的平均地层倾角余弦值网格tcosave,其中tcosave＝(tcost+tcosb)/2；
(9)将网格thc与网格tcosave相乘，得到流一段地层真厚度数据网格thr，其中thr＝thc×tcosave；
(10)利用流一段地层真厚度数据网格thr制作等值线图，作为目的层沉积时期的相对微古地貌形态。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。