地震剖面的经验模态分解方法

摘要

本发明涉及油气地震勘探领域，具体是一种地震剖面的经验模态分解方法，该方法涉及地震剖面二维瞬时属性分析，是一种从地震剖面分解单分量二维固有模态函数的计算方法。在二维经验模态分解过程中，采用时间方向两点邻域的极值检测方法和等式约束的紧框架域稀疏化曲面插值方法，解决地震剖面的极值检测和大规模不规则数据的包络面插值问题。本发明计算结果明显优于现有的二维经验模态分解快速方法，而计算时间大体相当。

主权项

一种地震剖面的经验模态分解方法，其特征在于，包括以下步骤：S1.初始化，r₀＝Iori，i＝1；给定迭代终止条件中ε1的值；Iori表示原始剖面数据，r代表当前残差剖面数据，i为第i个BIMF；S2.迭代筛分计算第i个BIMF；S21.初始化：h₀＝r_i‑1，j＝1；S22.检测h_j‑1的二维局部极小值点和极大值点，构造上下包络曲面E_up和E_down；S23.计算均值曲面m＝(E_up+E_down)/2；S24.h_j＝h_j‑1‑m；S25.判断当前h是否满足筛分终止条件，若满足，则得到第i个BIMFb_i＝h_j，跳转到步骤S3；否则，j＝j+1，跳转到步骤S22，继续筛分迭代；S3.计算当前残差，r_i＝r_i‑1‑b_i；S4.判断r_i是否满足BIMF终止条件，若满足，则分解结束，记I＝i，r_I是最终残余分量；否则i＝i+1，跳转到步骤S2；BEMD计算结果为：上述采用时间方向两点邻域的极值检测方法检测h_j‑1的极值；上述采用等式约束的紧框架域稀疏化曲面插值方法获取上下包络曲面。