8时，X和Y具有非常强的相关性；当0. 6彡I RI < 0. 8时，X和Y具有强相关性； 当0.4彡|R| <0.6时，X和Y具有中等相关性；当0.2彡|R| <0.4时，X和Y具有弱相 关性；当0 < I Rl < 0. 2时，X和Y不具有或者具有可忽视的相关性。
[0038] 3)对选出的对应相关系数大于0. 1的IMF分量逐道计算能量，查找出最大能量点 t_;頂F分量的能量计算为该頂F分量信号幅度值平方；根据地震信号特点，最大相关法选 取|R| > 0. 1的頂F分量作为主要体现煤层及烃源岩的頂F成分，实际地区可以根据该地 区数据特点可调节。
[0039] 4)对该单条地震道的頂F分量进行频谱分析，确定该地震道頂F信号的主频fd， 主频fd为幅值最大处对应的频率；令T d= 1/f d，确定煤层及烃源岩强振幅薄层在该頂F分 量中体现的顶和底的薄层的时间厚度为[tg-kiTd，t_-k2Td]，其中，h，k 2为常数系数，k i，k2 是通过对过井地震道进行分析，利用测井资料进行井震标定，确定与煤层及烃源岩强振幅 薄层对应的目的层段在地震反射波上体现的特征周期而得到的。
[0040] 5)计算[t.-kJd，t_-k2Td]范围内的数据的能量E s，能量Es是对[t .-Ic1Td, t_-k2Td]范围内各点幅值平方后求和，计算该地震道的平均能量
-Ic1Td, t_-k2Td]范围内的数据利用index系数将强振幅压制到与该地震道平均能量一致的 量级上。
[0041] 6)对该单条地震道相关系数小于0. 1的頂F分量进行步骤3)~步骤5)的处理， 将煤层及烃源岩强振幅对应的主要体现的频率段及强振幅下的MF信号进行处理，对该单 条地震道不同MF处理后的信号及保留的頂F信号相加，获得处理后的该单条地震道信号。
[0042] 7)对剩余各条原始地震道逐道重复上述步骤2)~步骤6)，实现对原始地震剖面 煤层及烃源岩强振幅的压制。
[0043] 本发明提供的这种基于经验模态分解检测地震标志层强反射振幅消除的方法在 对煤层及烃源岩强振幅信号进行压制的同时保留了信号中地层的信息，增强了与煤系地层 共存地层的弱信号分量，处理后弱信号的瞬时属性得到了加强；同时，能够更有效地增强邻 近区域的弱信号及微弱的油气响应特征，因此，能够用来进行二维及三维地震资料标志层 强反射振幅的自适应消除和邻近区域弱信号的增强，针对二维地震数据及三维地震数据进 行处理计算速度快，适合海量地震数据处理。
[0044] 实施例2 :
[0045] 在实施例1的基础上，本实施例以鄂尔多斯盆地苏里格地区的一条二维叠后偏移 地震剖面处理为例，该区域为致密砂岩储层。
[0046] 图1是苏里格地区的一条二维叠后偏移原始地震剖面图；先对过井地震道及其 EMD分解产生頂F信号，如图3所示；然后分析过井地震道经过EMD分解后的頂F信号与原 始过井地震到信号的相关性，结果如表1所示，进而得到利用本本发明方法处理后的各个 MF信号及最终重构的该地震道信号，如图4所示。而图2是利用本发明方法技术处理后的 地震剖面图。表1 :
[0048] 图5是原始地震剖面的瞬时幅度图，图6为利用本发明方法处理后地震剖面的瞬 时幅度图。图7为原始地震剖面的瞬时频率图，图8为利用本发明处理后地震剖面的瞬时 频率图。
[0049] 由上述各图可以看出，应用本发明方法后实现对原始地震剖面煤层及烃源岩强振 幅的压制，处理后弱信号的瞬时属性得到了加强。
[0050] 另外，在苏里格气田的勘探开发中，采用本发明方法在保持目的层弱反射的基础 上，针对中生界及古生界地震标志层强反射振幅消除，共完成二维地震资料452公里，有效 提高了目的层储层预测的符合率；并且该技术的应用有很好的技术和资料基础，是油气田 发展的迫切需求，能在鄂尔多斯盆地内及类似油气藏中大力推广，具有很好的应用前景。
[0051] 以上例举仅仅是对本发明的举例说明，并不构成对本发明的保护范围的限制，凡 是与本发明相同或相似的设计均属于本发明的保护范围之内。
【主权项】
1. 基于经验模态分解检测地震标志层强反射振幅消除的方法，其特征在于：包括如下 步骤： 1) 对原始地震道信号逐道进行经验模态分解，将每条原始地震道信号分解得到一系列 从高频到低频到趋势项的本征模态函数IMF; 2) 针对单条原始地震道信号，利用最大相关法选出煤层及烃源岩主要体现的頂F成 分，计算各个頂F分量与原始地震道信号的相关系数； 3) 对选出的对应相关系数大于0. 1的頂F分量逐道计算能量，查找出最大能量点 4) 对该单条地震道的頂F分量进行频谱分析，确定该地震道頂F信号的主频fd，令Td =l/fd，确定煤层及烃源岩强振幅薄层在该IMF分量中体现的顶和底的薄层的时间厚度为 tmax-k2Td]，其中，Ic1,k2为常数系数； 5) 计算[Ux-VTd,t_-k2Td]范围内的数据的能量Es，计算该地震道的平均能量Eare，令地震道平均能量一致的量级上； 6) 对该单条地震道相关系数小于0. 1的IMF分量进行步骤3)~步骤5)的处理，将煤 层及烃源岩强振幅对应的主要体现的频率段及强振幅下的MF信号进行处理，对该单条地 震道不同頂F处理后的信号及保留的頂F信号相加，获得处理后的该单条地震道信号； 7) 对剩余各条原始地震道逐道重复上述步骤2)~步骤6)，实现对原始地震剖面煤层 及烃源岩强振幅的压制。2. 如权利要求1所述的基于经验模态分解检测地震标志层强反射振幅消除的方法，其 特征在于：所述步骤1)中的经验模态分解是将原始地震道信号分解成一系列本征模态函 数頂F分量，即原始地震道信号X(t) =C1U)+C2 (t)+…+Cn(t)+Rn (t));其中，C1U)为第i 个IMF分量，i= 1~n，Rn⑴为余量。3. 如权利要求1所述的基于经验模态分解检测地震标志层强反射振幅消除的方法，其 特征在于：所述步骤2)中最大相关法计算两个变量X和Y之间的相关系数计算公式如下：望和方差；Y和〇Y分别是Y的期望和方差；COV为协方差；E为数学期望。4. 如权利要求1所述的基于经验模态分解检测地震标志层强反射振幅消除的方法，其 特征在于：所述步骤3)中IMF分量的能量为该IMF分量信号幅度值平方。5. 如权利要求1所述的基于经验模态分解检测地震标志层强反射振幅消除的方法，其 特征在于：所述步骤4)中主频fd为幅值最大处对应的频率。6. 如权利要求1所述的基于经验模态分解检测地震标志层强反射振幅消除的方法，其 特征在于：所述步骤4)中的^匕是通过对过井地震道进行分析，利用测井资料进行井震标 定，确定与煤层及烃源岩强振幅薄层对应的目的层段在地震反射波上体现的特征周期而得 到的。7. 如权利要求1所述的基于经验模态分解检测地震标志层强反射振幅消除的方法，其 特征在于：所述步骤5)中能量Es是对[t^x-Ic1Td,tmax-k2Td]范围内各点幅值平方后求和。
【专利摘要】本发明提供了一种基于经验模态分解检测地震标志层强反射振幅消除的方法，利用经验模态分解方法，结合相关函数分析法以及测井资料选取原始地震数据中最能体现煤层及烃源岩强振幅薄层特征的IMF分量，对优选出的主要IMF分量，查找其最大能量点，估算该IMF分量信号的主频，并在此基础上确定煤层及烃源岩强振幅薄层在该IMF分量中体现的顶和底的薄层的时间厚度，对该时间厚度范围内的IMF信号分量进行强振幅抑制，同时保留其余IMF分量不变，最后对所有处理过的IMF分量信号进行重构以实现对地震道信号中的地震标志层强反射振幅消除，增强煤系地层共存情况下地层信号的微弱变化，加强邻近区域弱油气响应特征。
【IPC分类】G01V1/36
【公开号】CN105044777
【申请号】CN201510379767
【发明人】赵玉华, 付金华, 薛雅娟, 田仁飞, 张盟勃, 陈娟, 高利东, 程思检
【申请人】中国石油天然气股份有限公司
【公开日】2015年11月11日
【申请日】2015年7月1日...