[发明专利]一种基于希尔伯特-黄变换的磁异常多特征信息提取方法

权利要求书

1.一种基于希尔伯特-黄变换的磁异常多特征信息提取方法，其特征在于步骤如下：

步骤1：采集磁场数据，通过磁力仪之间的测量差值获得磁场差值信号ΔB(t)；

步骤2：利用经验模式分解对磁场差值ΔB(t)进行自适应分解，获得多个固有模态函数项c(t)和一个残差项r(t)：

其中：N表示经验模式分解的层数；

步骤3：获得磁信号的固有模态函数后，对每个c_i(t)做希尔伯特变换：

步骤4：获得c_i(t)的解析信号z_i(t)：

其中：表示固有模态函数c_i(t)对应的瞬时幅值，θ_i(t)＝arctan(H_i(c_i(t))/c_i(t))表示固有模态函数c_i(t)对应的瞬时相位；

步骤5：根据瞬时相位θ_i(t)，计算对应的瞬时频率ω_i(t)

步骤6：磁信号的幅值可以表示为时间、瞬时频率的函数，获得磁信号幅值的希尔伯特谱H(ω,t)

步骤7：根据希尔伯特谱H(ω,t)，计算磁信号的边际谱h(ω)和瞬时能量密度IE(t)

IE(t)＝∫_ωH²(ω,t)dω (7)

其中：T表示信号的时间长度；

步骤8：根据磁信号的边际谱h(ω)和瞬时能量密度IE(t)，计算信号的多特征信息

其中：Mean(·)表示计算均值，Var(·)表示计算方差。

2.根据权利要求1所述的一种基于希尔伯特-黄变换的磁异常多特征信息提取方法，其特征在于步骤1中磁力仪的型号为CS-L光泵磁力仪。