整张 GPR 剖面上常有一层「每道都长一样」的横向条带,像照片蒙了灰雾。擦掉它,局部目标才显出来。

从哪来

有些成分在全测线每一道都几乎相同:

  • 直达波:发射天线直接窜到接收天线,每道开头都有;
  • 地表或天线面的反射、耦合响应;
  • 大范围的均匀层状反射。

它们在每个位置出现的时间和强度都差不多,画出来就是水平横条,盖住局部异常。

平均道就是背景

对每个时间样点,把所有道的值求平均,得到一条「平均道」。水平背景在每道都一样,在平均里相干保留;局部目标位置相位不同,互相抵消。

从每道减去这条平均道,背景没了,局部目标活下来。

统一公式:

$$D_{out}[t,x]=D[t,x]-\alpha\cdot\hat B[t,x],\qquad \alpha\in[0,1]$$

$\alpha=1$ 全减;$\alpha<1$ 部分减,适合既看点目标又想看连续界面的情况。

背景估计的几种方法

  • 全局 mean/median:所有道一起求均值或中位数。
  • 滑窗 mean/median:只取当前道附近窗口内的道求平均,适合背景沿测线缓慢变化。
  • SVD:取数据矩阵前几个奇异图,背景通常低秩。
  • f-k 滤波:在二维频谱里把低空间波数成分当背景去掉。

mean 还是 median?

一个 3×4 的迷你 B-scan:第 0 行是直达波(每道都是 50),第 1 行第 3 道有一根钢筋(30):

            道1  道2  道3  道4
行0(直达波) 50 50 50 50
行1(有钢筋) 2 30 3 1

用 mean:行 1 平均被钢筋拉到 9,减完后其他道出现 −7/−8 的凹坑。

用 median:行 1 中位数是 2.5,减完后钢筋 +27.5,其他道几乎不受影响。

有孤立强目标时,优先用 median。

最大的陷阱

如果你要找的目标本身就是连续水平界面,比如层间脱空、连续裂隙,它在每道都出现,也会被当成背景减掉。

  • 找钢筋、管线这种点目标:背景去除是好工具;
  • 找脱空、连续裂隙这种面目标:关掉背景去除,或者把窗调得很大。

判读要点:顶部直达波亮纹变淡、局部双曲线更清楚,就对了。同时要警惕横向连续反射有没有莫名消失。