[发明专利]地震剖面的经验模态分解方法

权利要求书

1.一种地震剖面的经验模态分解方法，其特征在于，包括以下步骤：

S1.初始化，r₀=Iori，i=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的极值；

上述采用等式约束的紧框架域稀疏化曲面插值方法获取上下包络曲面。

2.根据权利要求1所述的地震剖面的经验模态分解方法，其特征在于，所述采用时间方向两点邻域的极值检测方法具体为：若某点(x₀,y₀)处的波场值大于或等于其上、下两个相邻位置的波场值，则点(x₀,y₀)为剖面的极大值；相应地，若某点(x₀,y₀)处的波场值小于或等于其上、下两个相邻位置的波场值，则点(x₀,y₀)处波场值为剖面的极小值。