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

摘要

本发明涉及一种深层扇三角洲位置与形态的确定方法，包括以下步骤获得目标区域叠后地震记录，确定扇三角洲发育的研究区；建立研究区地质模型，给定扇三角洲及其下部和上部地层的岩性组成和修订步长的初始值；确定扇三角洲滤值窗口；使用扇三角洲滤值窗口对研究区叠后地震记录滤值，得到研究区点集；去除研究区点集中位于物源区的空间点，得到沉积区点集；对沉积区点集的空间标记分层，得到各分层点集；统计各分层点集内空间点的分布范围，搜索扇三角洲点集，确定具有该岩性组成的扇三角洲的位置和形状；修改研究区地质模型中扇三角洲及其下部和上部地层的岩性组成，确定具有不同岩性组成的扇三角洲的位置和形状；剔除重复结果，确定目标区域内具有所有可能岩性组成的扇三角洲的位置与形态。

主权项

一种深层扇三角洲位置与形态的确定方法，包括以下步骤：1)在目标区域进行地震数据采集，获得目标区域的叠后地震记录，在叠后地震记录的地震水平等时切片中确定出沟谷位置，选择沟谷发育区作为扇三角洲发育的研究区，并统计研究区内叠后地震记录的地震主频；2)建立研究区地质模型：设定中间扇三角洲的厚度；设定扇三角洲的砂岩体积比例为ps1＝a1‑it1，泥岩体积比例为pm1＝b1+it1；设定下部地层的砂岩体积比例为ps2＝a2+jt2，泥岩体积比例为pm2＝b2‑jt2；设定上部地层的砂岩体积比例为ps3＝a3+kt3，泥岩体积比例为pm3＝b3‑kt3；其中，a1、a2、a3分别为扇三角洲及其下部和上部地层的砂岩体积比例初始值，b1、b2、b3分别为扇三角洲及其下部和上部地层的泥岩体积比例初始值，t1、t2、t3分别为扇三角洲及其下部和上部地层的岩性组成修订步长，0≤a1、a2、a3、b1、b2、b3≤1，‑1≤t1、t2、t3≤1，且a1+b1＝a2+b2＝a3+b3＝1；i、j、k均为大于等于0的整数；并根据研究区实际地质情况，分别给定扇三角洲及其下部和上部地层的砂岩体积比例初始值a1、a2、a3，泥岩体积比例初始值b1、b2、b3，以及岩性组成修订步长t1、t2、t3；同时令i＝j＝k＝0；3)根据研究区地质模型中扇三角洲及其下部和上部地层的岩性组成，确定扇三角洲滤值窗口；4)使用扇三角洲滤值窗口对研究区的叠后地震记录进行滤值，只保留振幅在扇三角洲滤值窗口范围内的地震数据，得到研究区点集；5)随机设定多个物源区滤值窗口分别对研究区点集进行滤值，得到多个子点集，各子点集的分布范围的重叠次数最多的区域为物源区，在研究区点集中去除位于物源区的空间点，得到沉积区点集；6)观察沉积区点集的三维立体图，根据观测到的分层现象，对沉积区点集的空间进行标记分层，得到沉积区的各分层点集；7)在沉积区的各分层点集内分别统计空间点的分布范围，搜索空间点数量大于500且其中任意一点与至少两相邻点的水平距离小于2.5(x+y)的空间点的集合，即为扇三角洲点集，其中，x和y分别为叠后地震数据的主测线间距和联络测线间距；如果扇三角洲点集的所有空间点围成的形状是圆形或者椭圆形，则该扇三角洲点集的空间位置，即为具有与研究区地质模型中扇三角洲相同岩性组成的扇三角洲所在的位置，该扇三角洲点集的所有空间点围成的形状即为该扇三角洲的形状；8)循环修改研究区地质模型中扇三角洲及其下部和上部地层的岩性组成，并重复步骤3)～步骤7)，分别确定具有不同岩性组成的扇三角洲的位置和形状，直到修改不满足要求；9)统计分析得到的具有不同岩性组成的扇三角洲的位置与形状，剔除其中重复的结果，从而确定出目标区域内具有所有可能岩性组成的扇三角洲的位置与形态。