[发明专利]一种深层扇三角洲位置与形态的确定方法

说明书

技术领域

本发明涉及一种地质勘探方法，尤其涉及一种深层扇三角洲位置与形态的确定方法。

背景技术

现有技术中，通常通过常规的地震解释完成对扇三角洲的识别，该方法通过标定地震相，建立地震解释模式，可以将形态清楚、规模较大的扇三角洲解释识别出来。但是在深层地层，由于钻井少，地震主频低，地震品质差，深层扇三角洲并不容易识别；尤其是在扇三角洲进积或楔形等标志特征不明显的情况下，通过常规的地震解释无法确定扇三角洲的位置与形态。另外，现有技术中利用地震切片、属性分析等地球物理技术主要开展扇三角洲的面上特征刻画，但对扇三角洲的空间展布刻画十分困难。综上所述，针对深层扇三角洲识别尚未形成切实可行的一体化流程。

发明内容

针对上述问题，本发明的目的是提供一种深层扇三角洲位置与形态的确定方法，通过地震属性分析与处理、地震正演分析、滤值窗口处理等手段，建立从沟谷搜索到扇三角洲定位的完善流程，实现对深层无井且地震主频较低情况下扇三角洲位置和形态的快速有效识别。

为实现上述目的，本发明采取以下技术方案：一种深层扇三角洲位置与形态的确定方法，包括以下步骤：

1)在目标区域进行地震数据采集，获得目标区域的叠后地震记录，在叠后地震记录的地震水平等时切片中确定出沟谷位置，选择沟谷发育区作为扇三角洲发育的研究区，并统计研究区内叠后地震记录的地震主频；

2)建立研究区地质模型：设定中间扇三角洲的厚度；设定扇三角洲的砂岩体积比例为p_s1＝a₁-it₁，泥岩体积比例为p_m1＝b₁+it₁；设定下部地层的砂岩体积比例为p_s2＝a₂+jt₂，泥岩体积比例为p_m2＝b₂-jt₂；设定上部地层的砂岩体积比例为p_s3＝a₃+kt₃，泥岩体积比例为p_m3＝b₃-kt₃；其中，a₁、a₂、a₃分别为扇三角洲及其下部和上部地层的砂岩体积比例初始值，b₁、b₂、b₃分别为扇三角洲及其下部和上部地层的泥岩体积比例初始值，t₁、t₂、t₃分别为扇三角洲及其下部和上部地层的岩性组成修订步长，0≤a₁、a₂、a₃、b₁、b₂、b₃≤1，-1≤t₁、t₂、t₃≤1，且a₁+b₁＝a₂+b₂＝a₃+b₃＝1；i、j、k均为大于等于0的整数；

并根据研究区实际地质情况，分别给定扇三角洲及其下部和上部地层的砂岩体积比例初始值a₁、a₂、a₃，泥岩体积比例初始值b₁、b₂、b₃，以及岩性组成修订步长t₁、t₂、t₃；同时令i＝j＝k＝0；

3)根据研究区地质模型中扇三角洲及其下部和上部地层的岩性组成，确定扇三角洲滤值窗口；

4)使用扇三角洲滤值窗口对研究区的叠后地震记录进行滤值，只保留振幅在扇三角洲滤值窗口范围内的地震数据，得到研究区点集；

5)随机设定多个物源区滤值窗口分别对研究区点集进行滤值，得到多个子点集，各子点集的分布范围的重叠次数最多的区域为物源区，在研究区点集中去除位于物源区的空间点，得到沉积区点集；

6)观察沉积区点集的三维立体图，根据观测到的分层现象，对沉积区点集的空间进行标记分层，得到沉积区的各分层点集；