1.一种基于GWO优化小波阈值的大地电磁信号去噪方法，其特征在于：包括如下步骤：步骤1：将待去噪的大地电磁信号进行分段，并对每段信号均进行 小波 层变换，得到最大尺度低频小波系数 和各个尺度高频小波系数 ；

步骤2：基于GWO算法获取各个尺度下最优阈值作为最优分层阈值；

其中，将各个尺度下的阈值作为GWO算法中灰狼的位置参数；

步骤3：以每段信号为单元，利用对应的所述最优分层阈值进行去噪得到去噪后的大地电磁信号。

2.根据权利要求1所述的方法，其特征在于：步骤2中基于GWO算法获取各个尺度下最优阈值作为最优分层阈值的过程如下：步骤2‑1：参数初始化：灰狼种群规模 、搜索空间维度 、最大迭代次数 、搜索空间上界 ，下界 ；

步骤2‑2：初始化灰狼种群中每个灰狼个体的位置 ，所述位置 表示各个尺度下的阈值；

步骤2‑3：基于每个灰狼个体的当前位置 计算每个灰狼个体的适应度，其中，适应度越小对应阈值越优；

步骤2‑4：根据灰狼个体的适应度从小到大的顺序选取适应度靠前的3个灰狼个体的位置 作为 、 、 ；

步骤2‑5：基于所述 、 、 采用种群搜索策略接近猎物，进而更新灰狼个体的位置；

步骤2‑6：利用最优学习搜索方程学习步骤2‑5中每个灰狼个体的位置得到学习后的位置，再返回步骤2‑3进行下一次迭代计算；

其中，若是满足迭代终止条件，则将当前适应度最小的灰狼个体的位置 作为各个尺度下最优阈值；否则，继续迭代寻优。

3.根据权利要求2所述的方法，其特征在于：步骤2‑5中基于所述 、 、 采用种群搜索策略接近猎物，进而更新灰狼个体的位置的过程中，所述种群搜索策略的更新公式如下：

式中， 、 、 和 、 、 分别为由不同随机数生成的随机系数， 为灰狼个体的更新之前的位置， 、 、 表示将 、 、 分别作为猎物位置时，采用种群搜索策略更新的灰狼个体的位置； 、 、 表示将 、 、 分别作为猎物位置时，猎物与灰狼个体之间的距离； 为灰狼个体的更新之后的位置。

4.根据权利要求2所述的方法，其特征在于：步骤2‑6中利用最优学习搜索方程学习步骤2‑5中灰狼个体的位置得到学习后的位置时，所述最优学习搜索方程如下：其中， 、 分别为灰狼个体i在第 次、第 次迭代时的速度； 为灰狼个体i的最优位置； 为步骤2‑5中灰狼个体i位置更新后的最新位置； 为学习后灰狼个体i的位置； 为惯性权重， 、 为随机数， 、 为预设系数，为狼群中适应度最小的灰狼个体的位置。

5.根据权利要求2所述的方法，其特征在于：步骤2‑3中适应度的公式如下：其中， 表示基于位置 对应的分层阈值进行阈值处理后萎缩的小波系数个数，表示各个尺度下高频小波系数的总个数， 表示各个尺度下高频小波系数， 表示基于位置 对应的分层阈值进行阈值处理后的高频小波系数；

其中，阈值处理的公式如下：

式中， 表示位置 对应的分层阈值。

6.根据权利要求1所述的方法，其特征在于：步骤3中以每段信号为单元，利用对应的所述最优分层阈值进行去噪得到去噪后的大地电磁信号的过程，是引入的基追踪重构算法进行去噪，具体如下：

步骤3‑1：利用所述最优分层阈值对所述各个尺度高频小波系数进行阈值量化处理得到阈值处理后的高频小波系数，并引入所述最大尺度低频小波系数 进行重构得到估计信号 ；

步骤3‑2：分别对估计信号 与原始信号 进行 小波变换得到小波系数 和 ；

其中，所述 小波与所述 小波的小波基不同；

步骤3‑3：利用小波系数 构造Wiener滤波器 ，再利用所述Wiener滤波器 处理小波系数 得到原始信号 的稀疏表示 ；

步骤3‑4：构造 维观测矩阵 ，并利用所述观测矩阵 对所述稀疏表示 进行测量得到测量值 ；

步骤3‑5：利用基追踪重构算法对测量值 进行重构和求解得到最佳稀疏表示 ；

步骤3‑6：对所述最佳稀疏表示 进行小波逆变换 得到噪声轮廓 ，再将所述原始信号 减去噪声轮廓 得到对应一段去噪后的大地电磁信号，再将去噪后每一段大地电磁信号进行拼接。

7.根据权利要求6所述的方法，其特征在于：步骤3‑5中利用基追踪重构算法对测量值进行重构和求解得到最佳稀疏表示 的过程，基追踪重构模型如下：式中， 是对 的稀疏变换， 是稀疏变换矩阵， 为 的稀疏表示，是传感矩阵， ；

其中，求解所述基追踪重构模型得到重构后的最佳稀疏表示 ，所述最佳稀疏表示为重构的小波系数。

8.根据权利要求6所述的方法，其特征在于：步骤3‑3中Wiener滤波器 如下所示：其中， 为预设系数， 为维纳滤波器的统称， 表示系数的统称，将 作为代入 中得到Wiener滤波器 ；

所述 稀疏表示满足： 。

9.根据权利要求6所述的方法，其特征在于： 小波为db1小波基， 小波为db4小波基。

10.一种基于权利要求1‑9任一项所述方法的系统，其特征在于：包括：分段模块，用于将待去噪的大地电磁信号进行分段；

小波变换模块，用于对每段信号均进行 小波 层变换，得到最大尺度低频小波系数 和各个尺度高频小波系数 ；

分层阈值寻优模块，用于基于GWO算法获取各个尺度下最优阈值作为最优分层阈值；

去噪模块，用于以每段信号为单元，利用对应的所述最优分层阈值进行去噪得到去噪后的大地电磁信号。