[发明专利]基于经验模态分解检测地震标志层强反射振幅消除的方法

摘要

本发明提供了一种基于经验模态分解检测地震标志层强反射振幅消除的方法，利用经验模态分解方法，结合相关函数分析法以及测井资料选取原始地震数据中最能体现煤层及烃源岩强振幅薄层特征的IMF分量，对优选出的主要IMF分量，查找其最大能量点，估算该IMF分量信号的主频，并在此基础上确定煤层及烃源岩强振幅薄层在该IMF分量中体现的顶和底的薄层的时间厚度，对该时间厚度范围内的IMF信号分量进行强振幅抑制，同时保留其余IMF分量不变，最后对所有处理过的IMF分量信号进行重构以实现对地震道信号中的地震标志层强反射振幅消除，增强煤系地层共存情况下地层信号的微弱变化，加强邻近区域弱油气响应特征。

主权项

基于经验模态分解检测地震标志层强反射振幅消除的方法，其特征在于：包括如下步骤：1)对原始地震道信号逐道进行经验模态分解，将每条原始地震道信号分解得到一系列从高频到低频到趋势项的本征模态函数IMF；2)针对单条原始地震道信号，利用最大相关法选出煤层及烃源岩主要体现的IMF成分，计算各个IMF分量与原始地震道信号的相关系数；3)对选出的对应相关系数大于0.1的IMF分量逐道计算能量，查找出最大能量点t_max；4)对该单条地震道的IMF分量进行频谱分析，确定该地震道IMF信号的主频f_d，令T_d＝1/f_d，确定煤层及烃源岩强振幅薄层在该IMF分量中体现的顶和底的薄层的时间厚度为[t_max‑k₁T_d,t_max‑k₂T_d]，其中，k₁,k₂为常数系数；5)计算[t_max‑k₁T_d,t_max‑k₂T_d]范围内的数据的能量E_s，计算该地震道的平均能量E_ave，令对[t_max‑k₁T_d,t_max‑k₂T_d]范围内的数据利用index系数将强振幅压制到与该地震道平均能量一致的量级上；6)对该单条地震道相关系数小于0.1的IMF分量进行步骤3)～步骤5)的处理，将煤层及烃源岩强振幅对应的主要体现的频率段及强振幅下的IMF信号进行处理，对该单条地震道不同IMF处理后的信号及保留的IMF信号相加，获得处理后的该单条地震道信号；7)对剩余各条原始地震道逐道重复上述步骤2)～步骤6)，实现对原始地震剖面煤层及烃源岩强振幅的压制。