基于经验模态分解检测地震标志层强反射振幅消除的方法
【技术领域】
[0001] 本发明属于石油地球物理勘探技术领域，具体涉及一种基于经验模态分解检测地 震标志层强反射振幅消除的方法。
【背景技术】
[0002] 鄂尔多斯盆地是一个多构造体系、多旋回演化及多沉积类型的大型克拉通叠合盆 地。鄂尔多斯盆地陆相产气层主要集中在延长组和古生界石盒子组，由于鄂尔多斯盆地中 生界长7与古生界盒8、山1段致密油与致密气储集层紧靠长7烃源岩和古生界8#煤，而 长7烃源岩、古生界8#煤层在地震资料上都会形成强的反射同相轴，在地震剖面上表现为 强反射特征，从而导致中生界长7与古生界盒8、山1段致密油与致密气储集层在地震剖面 储层响应中淹没在长7烃源岩和古生界8#煤强反射中，导致紧邻这两套强反射的储层参数 受到影响，致使地震反演等技术得到的表征岩性、物性和含油气性的地球物理参数计算不 准确，影响了地震预测的成功率。
[0003] 要突出储层的地震响应，必须消减标志层强反射的能量。目前，针对煤层强反射消 减问题主要采用的方法有谱分解技术、匹配追踪技术和反Q滤波等方法。谱分解技术及匹 配追踪技术进行煤层强反射消除主要针对三维地震数据，利用煤层反射的低频特征，采用 估算出煤层的反射数据体然后进行强反射能量的消减，计算过程较为繁琐。反Q滤波是一 种补偿大地吸收衰减效应的技术，其目的是消除地震波在传播过程中振幅会衰减，相位会 畸变；目前效率高的反Q滤波方法都假设地下介质是一个常Q或者层状Q模型，这并不能反 映地下介质的真实情况；而假设Q随时间或深度连续变化的反Q滤波方法虽然接近地下介 质的真实情况，但是由于用到积分，计算效率较低。而目前针对标志层烃源岩强反射能量 的消除则涉及极少。
【发明内容】
[0004] 本发明的目的是提供一种基于经验模态分解检测地震标志层强反射振幅消除的 方法，通过借助经验模态分解能量检测法来检测地震剖面中标志层强反射振幅，有效地增 强邻近区域的弱信号及微弱的油气响应特征，进行二维及三维地震资料标志层强反射振幅 的自适应消除和邻近区域弱信号的增强，进而减弱地震标志层强反射振幅对邻近区域油气 检测的影响。
[0005] 本发明的技术方案是提供了一种基于经验模态分解检测地震标志层强反射振幅 消除的方法，包括如下步骤：
[0006] 1)对原始地震道信号逐道进行经验模态分解，将每条原始地震道信号分解得到一 系列从高频到低频到趋势项的本征模态函数IMF ;
[0007] 2)针对单条原始地震道信号，利用最大相关法选出煤层及烃源岩主要体现的頂F 成分，计算各个頂F分量与原始地震道信号的相关系数；
[0008] 3)对选出的对应相关系数大于0. 1的IMF分量逐道计算能量，查找出最大能量点 tmax，
[0009] 4)对该单条地震道的頂F分量进行频谱分析，确定该地震道頂F信号的主频fd， 令T d= 1/f d，确定煤层及烃源岩强振幅薄层在该IMF分量中体现的顶和底的薄层的时间厚 度为t_-k2T d]，其中，Ic1, k2为常数系数；
[0010] 5)计算[t.-kl，t_-k2Td]范围内的数据的能量E s，计算该地震道的平均能量
范围内的数据利用index系数将强振幅压制 到与该地震道平均能量一致的量级上；
[0011] 6)对该单条地震道相关系数小于0. 1的頂F分量进行步骤3)~步骤5)的处理，将 煤层及烃源岩强振幅对应的主要体现的频率段及强振幅下的MF信号进行处理，对该单条 地震道不同MF处理后的信号及保留的頂F信号相加，获得处理后的该单条地震道信号；
[0012] 7)对剩余各条原始地震道逐道重复上述步骤2)~步骤6)，实现对原始地震剖面 煤层及烃源岩强振幅的压制。
[0013] 上述步骤1)中的经验模态分解是将原始地震道信号分解成一系列本征模态函数 頂F分量，即原始地震道信号X(t) =C1U)+C2 (t)+…+Cn (t)+Rn(t));其中，C1U)为第i个 IMF分量，i = 1~n，Rn⑴为余量。
[0014] 上述步骤2)中最大相关法计算两个变量X和Y之间的相关系数计算公式如下：
，其中，μ JP 〇 x分别是X的期 望和方差；μ σ γ分别是Y的期望和方差；cov为协方差；E为数学期望。
[0015] 上述步骤3)中頂F分量的能量为该頂F分量信号幅度值平方。
[0016] 上述步骤4)中主频fd为幅值最大处对应的频率。
[0017] 上述步骤4)中的Iq，匕是通过对过井地震道进行分析，利用测井资料进行井震标 定，确定与煤层及烃源岩强振幅薄层对应的目的层段在地震反射波上体现的特征周期而得 到的。
[0018] 上述步骤5)中能量Es是对[t _;山，t_-k2Td]范围内各点幅值平方后求和。
[0019] 本发明的有益效果：
[0020] (1)本发明提供的这种基于经验模态分解检测地震标志层强反射振幅消除的方法 在对煤层及烃源岩强振幅信号进行压制的同时保留了信号中地层的信息，增强了与煤系地 层共存地层的弱信号分量，处理后弱信号的瞬时属性得到了加强。
[0021] (2)本发明提供的这种基于经验模态分解检测地震标志层强反射振幅消除的方法 能够更有效地增强邻近区域的弱信号及微弱的油气响应特征，因此，能够用来进行二维及 三维地震资料标志层强反射振幅的自适应消除和邻近区域弱信号的增强，针对二维地震数 据及三维地震数据进行处理计算速度快，适合海量地震数据处理。
[0022] 以下将结合附图对本发明做进一步详细说明。
【附图说明】
[0023] 图1是苏里格地区一条二维叠后偏移原始地震剖面图。
[0024] 图2是本发明方法处理后的地震剖面图。
[0025] 图3是过井地震道及其经验模态分解产生的頂F信号。
[0026] 图4是本发明处理后的各个頂F信号及最终重构的地震道信号。
[0027] 图5是原始地震剖面的瞬时幅度图。
[0028] 图6是本发明方法处理后地震剖面的瞬时幅度图。
[0029] 图7是原始地震剖面的瞬时频率图。
[0030] 图8是本发明方法处理后地震剖面的瞬时频率图。
【具体实施方式】
[0031] 实施例1 :
[0032] 本实施例提供了一种基于经验模态分解检测地震标志层强反射振幅消除的方法， 包括如下步骤：
[0033] 1)对原始地震道信号逐道进行经验模态分解，将每条原始地震道信号分解得到一 系列从高频到低频到趋势项的本征模态函数IMF ;
[0034] 经验模态分解（EMD)是将原始地震道信号分解成一系列本征模态函数頂F分量， 即原始地震道信号X (t) =C1U)+C2 (t)+…+Cn (t)+Rn (t));其中，C1U)为第i个頂F分量， i = 1~n，Rn(t)为余量。
[0035] 2)针对单条原始地震道信号，利用最大相关法选出煤层及烃源岩主要体现的頂F 成分，计算各个頂F分量与原始地震道信号的相关系数。
[0036] 由于地震资料的振幅、频率、能量等信息随煤层及烃源岩强振幅厚度的不同而有 很大的变化；而能量法是利用地震反射波解释煤层及烃源岩强振幅厚度的主要方法，其受 随机噪声影响较小，对于煤层、烃源岩强振幅及其相邻弱地震反射信号的差异性具有放大 效应，能更敏锐地捕获到地震剖面中由于煤层、烃源岩导致的最大能量信号段，这种能量的 差异与薄层的厚度和反射系数的变化都有关系，有利于煤层、烃源岩强振幅的消除；对于煤 系地层共存的情况下，为了在对煤层及烃源岩强振幅消除的同时，保留其共存地层的弱信 号特征，对于某条地震道，结合测井信息及最大相关法，优选出主要体现煤层及烃源岩信息 的頂F成分。
[0037] 最大相关法计算两个变量X和Y之间的相关系数计算公式如下：
，其中，μ JP σ x分别是X的 期望和方差；μ σ γ分别是Y的期望和方差；COV为协方差；E为数学期望。一般地，当 RI彡0.