[发明专利]CSAMT数据静态校正方法

权利要求书

1.一种CSAMT数据静态校正方法，其特征在于包括以下步骤：

步骤一、采集可控源音频大地电磁离散点视电阻率数据，拟合出可控源音频大地电磁原始信号s(t)，对其进行经验模态分解(Empirical Mode Decomposition，简称EMD)得到n个本征模态函数c₁(t)、c₂(t)、...c_i(t)、...c_n(t)，并将各本征模态函数作为IMF分量；

步骤二、根据观测点与参考点的全部频点的视电阻率值之间的互相关系数r_xy识别静态效应,如果存在静态效应，则按照步骤三对原始信号s(t)进行重构；

步骤三、根据步骤一得到的各IMF分量对原始信号s(t)进行重构，得到静态校正后的可控源音频大地电磁重构信号s'(t)；

对原始信号s(t)进行重构的方法如下：

(1)采用基于高斯白噪声的自适应阈值作为重构阈值，设第i个IMF分量c_i(t)的重构阈值为T_i：

其中，c表示阈值系数，取值范围为0.5-0.7；M表示选取的信号长度，即各IMF分量的离散点个数，E_i表示第i个IMF分量c_i(t)的能量强度，其计算公式为：

其中，是第一个IMF分量c₁(t)的方差；

(2)用h_(ik)(t)表示第i个IMF分量c_i(t)在时间t_k处的离散值，k＝1,2,...M,则第i个IMF分量c_i(t)对应的重构分信号为

则可控源音频大地电磁重构信号为s'(t)：

2.根据权利要求1所述的CSAMT数据静态校正方法，其特征在于所述的阈值系数c＝0.7。

3.根据权利要求1所述的CSAMT数据静态校正方法，其特征在于所述的步骤一中，对CSAMT的视电阻率数据进行经验模态分解得到多个本征模态函数分量的方法如下：

(一)、找出可控源音频大地电磁原始信号s(t)所有的极大值点和极小值点；

(二)、采用三次样条插值法分别将所有的极大值点和极小值点连接起来形成上包络线S₁、下包络线S₂：用原始信号s(t)各离散点视电阻率值ρ减去上包络线S₁与下包络线S₂的对应点均值,得到多个离散差值点并对其进行拟合得到第一个差值拟合曲线h₁₁(t)；将第一个差值拟合曲线h₁₁(t)所有的极大值点和极小值点连接起来形成差值信号上包络线H₁₁、下包络线H₁₂；考察第一个差值拟合曲线h₁₁(t)是否满足IMF的两个条件,即：极值点数目(包括极大值和极小值)和过零点数目相等或最多相差1个；任意点对应的差值信号上包络线H₁₁上的离散点值与差值信号下包络线H₁₂上的离散点值的平均值为0；如果满足则将第一个差值拟合曲线h₁₁(t)作为第1个本征模态函数c₁(t)，转到步骤三；否则用第一个差值拟合曲线h₁₁(t)替代原始信号s(t)重复上述筛选过程，直至从原始信号s(t)中筛选出满足IMF的两个条件的差值拟合曲线并将其作为第1个本征模态函数c₁(t)；

(三)将剩下的信号r₁(t)＝s(t)-c₁(t)作为新的信号，按照步骤(一)、(二)，得到第二个本征模态函数c₂(t)；依次类推，最后得到n个本征模态函数c₁(t)、c₂(t)、...c_i(t)、...c_n(t)和一个不可分解的残余分量r_n(t)：

n个本征模态函数c₁(t)、c₂(t)、...c_i(t)、...c_n(t)作为IMF分量。