一种基于PSR-fastICA的频谱激电信号降噪方法及装置制造方法及图纸

本发明专利技术实施例公开一种基于PSR

【技术实现步骤摘要】
一种基于PSR
‑
fastICA的频谱激电信号降噪方法及装置


[0001]本专利技术实施例涉及信号处理
，具体涉及一种基于PSR
‑
fastICA的频谱激电信号降噪方法。

技术介绍

[0002]新中国成立以来，我国地球物理勘查技术在资源勘探中发挥了重要的作用，同时也获得了十足的发展。但是，我国目前处于经济快速发展阶段，资源需求与供给的矛盾日益突出：部分金属矿产的需求持续增长，但相对匮乏；油气资源对外依存度逐年攀升；由于新能源与新材料产业的发展，国内对三系资源、锂矿等需求大幅度上升。在浅表层发现矿产大幅减少的今天，迫切需要将目标转向深部资源。因此，大力发展地球物理勘探技术，在迫切需要资源的今天，极具战略意义。
[0003]而在对深部地层进行探测时，获取的频谱激电信号常有多种干扰噪声，如自然电位干扰、工频干扰、接受电缆干扰的高频耦合干扰、电子器件噪声等部分组成。其中目标信号频率在2
‑2～22Hz频段，与噪声相互叠加，在目标频段产生混叠。在这种条件下，激电信号的相位谱和幅度谱很难获得较高的精本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于PSR
‑
fastICA的频谱激电信号降噪方法，其特征在于，包括步骤：通过C
‑
C法对均值化操作后的频谱激电信号进行相空间重构，得到虚拟高维观测信号序列；将所述虚拟高维观测信号序列作为输入，进行快速独立成分分析；获取与发射信号相关性最大的一路信号，进行频谱分析。2.根据权利要求1所述的基于PSR
‑
fastICA的频谱激电信号降噪方法，其特征在于，所述C
‑
C法的计算式为：其中：当x＜0时，θ(x)＝0；当x＞0时，θ(x)＝1；||X
i
‑
X
j
||
(2)
表示二阶范数；m表示相空间重构参数中的嵌入维数；t表示相空间重构参数中的延迟系数；r表示搜寻半径；N表示时间序列长度；M表示相空间中的总点数；i和j均表示相空间中的某个点数；在C
‑
C法中，搜寻半径的计算式为：r(m)＝r
·
log(m+1)；其中：r为搜寻半径；m为相空间重构参数中的嵌入维数。3.根据权利要求2所述的基于PSR
‑
fastICA的频谱激电信号降噪方法，其特征在于，所述的通过C
‑
C法对均值化操作后的频谱激电信号进行相空间重构，得到虚拟高维观测信号序列具体包括步骤：依据分块平均策略获取嵌入维数、时间序列长度以及搜寻半径的范围；采用C
‑
C法获取确定的延迟系数与延迟时间窗口；根据延迟时间窗口、嵌入维数和延迟系数之间存在的定量关系获取确定的嵌入维数；通过确定的嵌入维数和延迟系数进行相空间重构，得到虚拟高维观测信号序列。4.根据权利要求3所述的一种基于PSR
‑
fastICA的频谱激电信号降噪方法，其特征在于，所述的分块平均策略为：其中：m表示相空间重构参数中的嵌入维数；N表示时间序列长度；r表示搜寻半径；t表示相空间重构参数中的延迟系数。5.根据权利要求4所述的基于PSR
‑
fastICA的频谱激电信号降噪方法，其特征在于，所述虚拟高维观测信号的序列为：X＝{X
i
|X
i
＝[x
i
,x
i+τ
,
…
,x
i+(m
‑
i)τ
]
T
,i＝1,2,
…
,M}；M＝N
‑
(m
‑
1)τ；其中：X为高维观测信号的序列；i为相空间中的点数；M为相空间中的总点数；N为时间序列长度，m为相空间重构参数中的嵌入维数，τ为相空间重构参数中的延迟系数，T表示转置矩阵。6.根据权利要求1

【专利技术属性】
技术研发人员：李志华，岑威，柯鑫，
申请(专利权)人：中国地质大学武汉，
类型：发明
国别省市：