探地雷达探测地面干扰快速消除方法

摘要

本发明公布了一种探地雷达探测地面干扰快速消除方法。本发明首先进行野外干扰源调查，获得干扰源类型及其与探地雷达测线的位置关系。其次.计算干扰时距曲线，并与探地雷达实测剖面进行对比，在实测剖面上识别出干扰信号。然后对干扰时距曲线进行离散化处理，并从实测剖面提取出包含干扰信号的数据，最后计算干扰数据的标准值，利用标准值进行干扰消除处理计算，得到替换数据。本发明能够快速消除探地雷达实测数据中的地面干扰信号，进而增强数据的信噪比，使被干扰信号干扰的有用数据得到展示，为最终的解译提供更多的有用信息，从而提高探地雷达的探测效果。

主权项

探地雷达探测地面干扰快速消除方法，其特征在于该方法包括以下步骤：步骤1.野外干扰源调查，获得干扰源类型及其与探地雷达测线的位置关系；步骤2.计算干扰时距曲线；步骤3.将计算所得的干扰时距曲线与探地雷达实测剖面进行对比，从而在实测剖面上识别出干扰信号；步骤4.选取处理时窗，设时窗的大小为TW，干扰信号的延续时间为L，则L≤TW≤1.5L，干扰信号的延续时间，从探地雷达实测剖面中直接拾取；步骤5.对干扰时距曲线进行离散化处理，用V(x_i,t_k),i＝1,2,3...n,k＝1,2,3...m来表示探地雷达实测剖面上(x_i,t_k)位置的信号值，其中x_i为第i道探地雷达信号的位置坐标，n为总道数，t_k为第k个时间采样点的时间坐标，m为采样点数；将干扰时距曲线离散化为(x_i,t_i)，t_i为步骤2干扰曲线计算所得的时间坐标；步骤6.从实测剖面提取出包含干扰信号的数据，可得V(x_i,t_i+j×dt)；j＝0,1,2....TW_n，dt为时间采样率，TW_n为处理时窗除以时间采样率所得的处理采样点数；步骤7.计算干扰数据的标准值V_b(j×dt)$V_b(j\times \mathrm{dt})=(\underset{i=1}{\overset{n}{\Sigma }}V(x_i,t_i+j\times \mathrm{dt})\times {t_i}^2)/n$步骤8.利用标准值进行干扰消除处理计算，得到替换数据V_n(x_i,t_i+j×dt)，V_n(x_i,t_i+j×dt)＝V(x_i,t_i+j×dt)‑V_b(j×dt)/t_i²步骤9.将步骤6所提取的数据用干扰消除处理计算所得的替换数据替换。