编码源电磁测深法获取大地全区视电阻率谱的方法及装置的制造方法
【技术领域】
[0001] 本发明涉及信号采集与处理领域、系统辨识领域及地球物理电磁测深法勘探领 域，具体而言，涉及一种编码源电磁测深法获取大地全区视电阻率谱的方法及装置。
【背景技术】
[0002] 由于不同频率的电磁场在大地中的趋肤深度不同，因此在电磁测深领域中可以通 过在地表测量电磁场分量，获取待测大地的全区视电阻率，根据全区视电阻率随频率的变 化，推测待测大地的电性深度分布情况。
[0003] 当前，相关技术中通常采用扫频法来获取待测大地的视电阻率。在应用扫频法时， 预先设置多个不同的频率，在勘测时依次向待测大地输入预设的多个不同频率的电流信 号，根据待测大地对多个不同频率的电流信号对应的电磁响应输出信号进行观测，获得待 测大地的视电阻率。
[0004] 由于依次向待测大地输入预设的多个不同频率的电流信号，如此一个频率一个频 率的进行扫频勘测，效率很低，且没有采用相关辨识技术，抗干扰能力很低。
【发明内容】
[0005] 有鉴于此，本发明实施例的目的在于提供一种编码源电磁测深法获取大地全区视 电阻率谱的方法及装置，提高了获取全区视电阻率的效率，及提高了获取的全区视电阻率 的频率分辨率，且能有效压制工频干扰，增强了抗干扰能力。
[0006] 第一方面，本发明实施例提供了一种编码源电磁测深法获取大地全区视电阻率谱 的方法，所述方法包括：
[0007] 向待测大地输入电流信号，采集所述待测大地输出的电磁场响应信号时间序列， 所述电流信号为逆重复Μ序列且周期为工频周期的偶数倍；
[0008] 根据所述电流信号和所述电磁场响应信号时间序列，通过循环互相关方式获取所 述待测大地的频率响应；
[0009] 根据所述待测大地的频率响应及接收电磁场分量，获取所述待测大地的全区视电 阻率谱。
[0010] 结合第一方面，本发明实施例提供了上述第一方面的第一种可能的实现方式，其 中，
[0011]所述根据所述电流信号和所述电磁场响应信号时间序列，通过循环互相关方式获 取所述待测大地的频率响应，包括：
[0012] 按照所述电流信号的时间序列周期和采样率触发软件程序生成参考信号，将所述 电磁场响应信号时间序列及所述电流信号分别与所述参考信号进行循环互相关，得到所述 电磁场响应信号时间序列对应的互相关时间序列及所述电流信号对应的互相关时间序列；
[0013] 对所述电磁场响应信号时间序列对应的互相关时间序列及所述电流信号对应的 互相关时间序列分别进行快速傅立叶变换，得到所述电磁场响应信号时间序列的互功率谱 及所述电流信号的互功率谱；
[0014] 根据所述电磁场响应信号时间序列的互功率谱及所述电流信号的互功率谱，获取 所述待测大地的频率响应。
[0015] 结合第一方面，本发明实施例提供了上述第一方面的第二种可能的实现方式，其 中，
[0016] 所述接收电磁场分量为水平电场时间序列；
[0017] 所述根据所述待测大地的频率响应及接收电磁场分量，获取所述待测大地的全区 视电阻率谱，包括：
[0018] 根据所述待测大地的频率响应，通过公式（1)获得所述待测大地的阻抗谱与接收 电磁场分量及大地电阻率的关系式；
[0019] 对公式（1)对应的所述待测大地的阻抗谱对应的关系式进行变形，得到公式(2)所 示的所述待测大地的全区视电阻率谱；
[0022] 其中，在公式（1)和(2)中，ΖΕχ(ω)为大地阻抗，p为均匀大地的电阻率，ω为频率， dL为水平电偶源的两极间距，ΜΝ为接收电极的电极间距，r为接收电极的中点与水平电偶源 的两极中点之间的距离，
，炉为接收电极的中点至 水平电偶源的两极中点之间的连线与水平电偶源的两极连线的夹角，i为虚数单位，
为均匀大地的波数，μ〇为自由空间的磁导率，ΡΕχ( ω )为所述全区视电阻率。
[0023] 结合第一方面，本发明实施例提供了上述第一方面的第三种可能的实现方式，其 中，
[0024] 所述接收电磁场分量为垂直磁场时间序列；
[0025] 所述根据所述待测大地的频率响应及接收电磁场分量，获取所述待测大地的全区 视电阻率谱，包括：
[0026] 根据所述待测大地的频率响应，通过公式(3)获得所述待测大地的频率响应与接 收电磁场分量及大地电阻率的关系式；
[0027]对公式(3)对应的所述频率响应对应的关系式进行变形，得到公式(4)所示的所述 待测大地的全区视电阻率谱；
[0029] 其中，在公式(3)和(4)中，ΖΗζ( ω )为所述频率响应，ρ为均匀大地的电阻率，ω为频 率，dL为水平电偶源的两极间距，r为接收电极的中点与水平电偶源的两极中点之间的距 离，μ〇为自由空间的磁导率，i为虚数单位，#为接收电极的中点至水平电偶源的两极中点之 间的连线与水平电偶源的两极连线的夹角，
为均匀大地的波数，Ρηζ( ω )为所述全区视电阻率。
[0030] 结合第一方面，本发明实施例提供了上述第一方面的第四种可能的实现方式，其 中，所述方法还包括：
[0031] 根据所述待测大地的频率响应，获取所述待测大地的传递函数，由所述传递函数 计算所述待测大地的阶跃响应或冲激响应，根据所述阶跃响应的晚期渐近值或所述冲激响 应的峰值时刻，计算所述待测大地的视电阻率。
[0032] 第二方面，本发明实施例提供了一种编码源电磁测深法获取大地全区视电阻率谱 的装置，所述装置包括：
[0033] 采集模块，用于向待测大地输入电流信号，采集所述待测大地输出的电磁场响应 信号时间序列，所述电流信号为逆重复Μ序列且周期为工频周期的偶数倍；
[0034] 第一获取模块，用于根据所述电流信号和所述电磁场响应信号时间序列，通过循 环互相关方式获取所述待测大地的频率响应；
[0035] 第二获取模块，用于根据所述待测大地的频率响应及接收电磁场分量，获取所述 待测大地的全区视电阻率谱。
[0036] 结合第二方面，本发明实施例提供了上述第二方面的第一种可能的实现方式，其 中，
[0037]所述第一获取模块包括：
[0038] 循环互相关单元，用于按照所述电流信号的时间序列周期和采样率触发软件程序 生成参考信号，将所述电磁场响应信号时间序列及所述电流信号分别与所述参考信号进行 循环互相关，得到所述电磁场响应信号时间序列对应的互相关时间序列及所述电流信号对 应的互相关时间序列；
[0039] 傅立叶变换单元，用于对所述电磁场响应信号时间序列对应的互相关时间序列及 所述电流信号对应的互相关时间序列分别进行快速傅立叶变换，得到所述电磁场响应信号 时间序列的互功率谱及所述电流信号的互功率谱；
[0040] 获取单元，用于根据所述电磁场响应信号时间序列的互功率谱及所述电流信号的 互功率谱，获取所述待测大地的频率响应。
[0041] 结合第二方面，本发明实施例提供了上述第二方面的第二种可能的实现方式，其 中，
[0042] 所述接收电磁场分量为水平电场时间序列；
[0043] 所述第二获取模块包括：
[0044] 第一获得单元，用于根据所述待测大地的频率响应，通过公式（1)获得所述待测大 地的阻抗谱与接收电磁场分量及大地电阻率的关系式；
[0045] 第一计算单元，用于对公式（1)对应的所述待测大地的阻抗谱对应的关系式进行 变形，得到公式(2)所示的所述待测大地的全区视电阻率谱；
[0048]其中，在公式（1)和(2)中，ΖΕχ(ω)为大地阻抗，p为均匀大地的电阻率，ω为频率， dL为水平电偶源的两极间距，ΜΝ为接收电极的电极间距，r为接收电极的中点与水平电偶源 的两极中点之间的距离，
P为接收电极的中点至 水平电偶源的两极中点之间的连线与水平电偶源的两极连线的夹角，i为虚数单位，
为均匀大地的波数，μ〇为自由空间的磁导率，ΡΕχ( ω )为所述全区视电阻率。
[0049]结合第二方面，本发明实施例提供了上述第二方面的第三种可能的实现方式，其 中，
[0050] 所述接收电磁场分量为垂直磁场时间序列；
[0051] 所述第二获取模块包括：
[0052] 第二获得单元，用于根据所述待测大地的频率响应，通过公式(3)获得所述待测大 地的频率响应与接收电磁场分量及大地电阻率的关系式；
[0053]第二计算单元，用于对公式(3)对应的所述频率响应对应的关系式进行变形，得到 公式(4)所示的所述待测大地的全区视电阻率谱；
[0055] 其中，在公式(3)和(4)中，Ζηζ( ω )为所述频率响应，ρ为均匀大地的电阻率，ω为频 率，dL为水平电偶源的两极间距，r为接收电极的中点与水平电偶源