，其中，X和y分别为叠后地震数据的主测线(inline)间距和联络测线(crossline)间距。如果扇三角洲点集的所有空间点围成的是圆形、椭圆形或者其他近似圆形或椭圆的形状，则该扇三角洲点集的空间位置即为具有与扇三角洲模型相同岩性组成的扇三角洲2所在的范围，该扇三角洲点集的空间边界形成的形状即为该扇三角洲2的形状。
[0061]8)循环修改研究区地质模型中扇三角洲及其下部和上部地层的岩性组成，确定具有不同岩性组成的扇三角洲的位置和形状:
[0062]①令inew=i+l，如果O < a1-1newti < I，返回步骤3);否则，进行下一步。
[0063]②令jnew= j + 1，i=0，如果O < a2+jnewt2 < I，返回步骤3);否则，进行下一步。
[0064]③令knew=k+l，i = j = 0，如果O < a3+knewt3 < I，返回步骤3);否则，进行下一步。
[0065]9)统计分析得到的具有不同岩性组成的扇三角洲的位置与形状，剔除其中重复的结果，从而确定出目标区域内具有所有可能岩性组成的扇三角洲2的位置与形态。
[0066]上述各实施例仅用于说明本发明，其中各部件的结构、设置位置及其连接方式等都是可以有所变化的，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。
【主权项】
1.一种深层扇三角洲位置与形态的确定方法，包括以下步骤: 1)在目标区域进行地震数据采集，获得目标区域的叠后地震记录，在叠后地震记录的地震水平等时切片中确定出沟谷位置，选择沟谷发育区作为扇三角洲发育的研究区，并统计研究区内叠后地震记录的地震主频； 2)建立研究区地质模型:设定中间扇三角洲的厚度；设定扇三角洲的砂岩体积比例为?31 = 31-；[1:1，泥岩体积比例为口1111 = 131 + ；[1:1;设定下部地层的砂岩体积比例为口32 = 32+九2，泥岩体积比例为Pm2 = b2-j t2 ；设定上部地层的砂岩体积比例为Ps3 = a3+kt3，泥岩体积比例为Pm3 = b3-kt3;其中，a1、a2、a3分别为扇三角洲及其下部和上部地层的砂岩体积比例初始值，Hb3分别为扇三角洲及其下部和上部地层的泥岩体积比例初始值，。山山分别为扇三角洲及其下部和上部地层的岩性组成修订步长，O < a1、a2、a3、b1、b2、b3 < I ,-1 < ?1ν?2ν?3 <I，且ai+bi = a2+b2 = a3+b3 = l; 1、j、k均为大于等于O的整数； 并根据研究区实际地质情况，分别给定扇三角洲及其下部和上部地层的砂岩体积比例初始值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)统计分析得到的具有不同岩性组成的扇三角洲的位置与形状，剔除其中重复的结果，从而确定出目标区域内具有所有可能岩性组成的扇三角洲的位置与形态。2.如权利要求1所述的一种深层扇三角洲位置与形态的确定方法，其特征在于，所述步骤3)根据研究区地质模型中扇三角洲及其下部和上部地层的岩性组成确定扇三角洲滤值窗口的方法，具体包括以下步骤: 首先，根据研究区地质模型中扇三角洲及其下部和上部地层的岩性组成，分别计算扇三角洲及其下部和上部地层的地震波速度值； 然后，根据研究区地质模型中扇三角洲及其上下地层的地震波速度值和研究区内叠后地震记录的地震主频，在研究区地质模型内进行地震波垂直入射反射正演模拟，得到研究区地质模型的模拟地震反射剖面； 再然后，设定不同振幅范围的滤值窗口对模拟地震反射剖面内的地震波进行滤值，得到与不同振幅范围的滤值窗口对应的地震反射结果； 最后，分析地震反射结果的分布范围，选择能够屏蔽研究区地质模型的扇三角洲范围之外的地震反射信息的滤值窗口，作为扇三角洲滤值窗口。3.如权利要求2所述的一种深层扇三角洲位置与形态的确定方法，其特征在于，通过加权平均的方式计算扇三角洲及其下部和上部地层的地震波速度值，加权平均计算公式为:V — Vm X pm+Vs X ps 式中，V表示扇三角洲或者其下部或上部地层的地震波速度;分别表示泥岩和砂岩的地震波速度;Pm和Ps分别表不泥岩和砂岩占该地层的体积比例，满足Pm+Ps= I。4.如权利要求1或2或3所述的一种深层扇三角洲位置与形态的确定方法，其特征在于，所述步骤8)修改研究区地质模型中扇三角洲及其下部和上部地层的岩性组成，并分别确定具有不同岩性组成的扇三角洲的位置和形状的方法，具体包括以下步骤: ①令inew=i+l，如果O< a1-1newti < I，返回步骤3);否则，进行下一步。 ②令jnew=j + 1，i=0，如果O < a2+jnewt2 < I，返回步骤3);否则，进行下一步。 ③令knew=k+l，i= j = 0，如果O < a3+knewt3 < I，返回步骤3);否则，进行下一步。5.如权利要求4所述的一种深层扇三角洲位置与形态的确定方法，其特征在于，所述步骤①修改为-令^^^ + ^如果。.？ < a1-1newtl < I，返回步骤3);否则，进行下一步。6.如权利要求1或2或3或5所述的一种深层扇三角洲位置与形态的确定方法，其特征在于，所述步骤5)中在识别物源区时，若已有目标区域的叠后地震数据的构造解释成果，则位于凸起部位的区域为物源区。7.如权利要求4所述的一种深层扇三角洲位置与形态的确定方法，其特征在于，所述步骤5)中在识别物源区时，若已有目标区域的叠后地震数据的构造解释成果，则位于凸起部位的区域为物源区。
【专利摘要】本发明涉及一种深层扇三角洲位置与形态的确定方法，包括以下步骤：获得目标区域叠后地震记录，确定扇三角洲发育的研究区；建立研究区地质模型，给定扇三角洲及其下部和上部地层的岩性组成和修订步长的初始值；确定扇三角洲滤值窗口；使用扇三角洲滤值窗口对研究区叠后地震记录滤值，得到研究区点集；去除研究区点集中位于物源区的空间点，得到沉积区点集；对沉积区点集的空间标记分层，得到各分层点集；统计各分层点集内空间点的分布范围，搜索扇三角洲点集，确定具有该岩性组成的扇三角洲的位置和形状；修改研究区地质模型中扇三角洲及其下部和上部地层的岩性组成，确定具有不同岩性组成的扇三角洲的位置和形状；剔除重复结果，确定目标区域内具有所有可能岩性组成的扇三角洲的位置与形态。
【IPC分类】G01V1/30
【公开号】CN105676284
【申请号】CN201610044509
【发明人】康波, 徐建永, 吴克强, 何玉平, 武爱俊, 印斌浩, 李凡异, 肖伶俐, 陈雯雯, 刘子玉
【申请人】中国海洋石油总公司, 中海油研究总院
【公开日】2016年6月15日
【申请日】2016年1月22日