3. 折射波时距曲线特征

折射静校正的第一步,是把野外记录到的初至时间变成可解释的“地下结构信息”。这一步依赖正演:给定一个速度界面模型,算出折射波应该什么时候到达检波器。只有正演关系清楚,后续反演初至时间求风化层厚度、求低降速带速度才有依据。本节围绕单层与多层折射波的走时方程展开,重点讲清截距时间和视速度这两个在初至反演中反复使用的参数。

一、单层折射波的走时方程

设地表下存在一个水平界面,上覆层速度为 v₁、厚度为 h,下伏层速度为 v₂,且 v₂ > v₁。当炮检距 x 足够大,射线以临界角 θc 入射到界面,沿界面滑行一段后再以临界角折返地面。由斯涅尔定律 sinθc = v₁/v₂。沿这条路径的走时为:

t(x) = x / v₂ + 2h·cosθc / v₁
     = x / v₂ + 2h·√(v₂² − v₁²) / (v₁·v₂)

这是一条直线方程,斜率倒数即视速度 v₂,直线在时间轴上的截距为

t₀ = 2h·cosθc / v₁ = 2h·√(v₂² − v₁²) / (v₁·v₂)

由此可直接反解厚度:h = t₀·v₁·v₂ / (2√(v₂² − v₁²))。这就是单层折射法求风化层厚度的基本公式,也是折射静校正最核心的正演—反演关系。

注意:上述公式只在 x 大于“交叉距离”后成立。近炮点段初至是直达波,时距曲线为斜率 1/v₁ 的直线;远炮点段才是折射波。实际拾取初至时必须区分这两段,否则拟合出的截距时间会偏大。

二、多层折射波与视速度的物理意义

当存在多个速度递增界面(如风化层—低速层—高速层),每一层折射波在远炮点段都会形成一条近似直线。第 n 层折射波在时距曲线上的斜率倒数为该层的“视速度”,它等于该层真实速度(水平层状介质下),而截距时间则累积了上覆各层的厚度贡献。以三层为例,第二界面折射波的走时可写成:

t(x) = x / v₃ + t₀₂
t₀₂ = 2h₁·√(v₃² − v₁²)/(v₁·v₃) + 2h₂·√(v₃² − v₂²)/(v₂·v₃)

可见截距时间不是单一厚度信息,而是各层厚度与速度对比的综合反映。这正是多层折射反演需要“逐层剥离”的原因:先用第一段直线求 v₂ 和 h₁,再从第二段截距中扣除第一层贡献,才能解出 h₂。

单层与多层折射波时距曲线及参数关系 地质模型(水平层状) v₁ 风化层 厚度 h₁ v₂ 低速层 厚度 h₂ v₃ 高速层 v₃ > v₂ > v₁ 临界折射 深部折射 t–x 时距曲线 x t 直达波 1/v₁ v₂ 段 t₀₁ v₃ 段 t₀₂ 交叉距离 斜率 = 1/视速度;截距 = 累积厚度贡献

视速度的物理意义需要强调两点。第一,它由时距曲线远炮点段的斜率倒数得到,反映的是折射界面以下高速层的速度,而不是上覆层速度。第二,在倾斜界面或横向变速情况下,视速度会偏离真实速度,此时需要采用相遇观测系统或互换法来校正。实际处理中,常把初至时间按炮检距分段线性拟合,用各段斜率判断速度分层,用截距时间估算厚度,这就是“t₀ 法”静校正的基础。

三、从正演到反演的实操要点

以一个教学示例说明:假设风化层 v₁ = 600 m/s、厚度 h₁ = 12 m,下伏层 v₂ = 1800 m/s。计算临界角 sinθc = 600/1800 = 0.333,θc ≈ 19.47°。截距时间 t₀ = 2×12×cos19.47°/600 ≈ 0.0377 s,即约 37.7 ms。若野外初至拾取得到远炮点段斜率对应视速度 1800 m/s、截距 38 ms,则反算 h₁ ≈ 12.1 m,与模型吻合。这个“正演—拾取—反演”的闭环,是检验静校正参数是否合理的直接手段。

常见误区:把近炮点直达波段也纳入折射波直线拟合,会导致截距时间系统性偏大、厚度高估;另外,若低速层速度接近下伏层速度,临界角接近 90°,折射波很弱甚至不出现,此时单层公式失效,需要结合微测井或层析反演约束。

总结来说,单层与多层折射波走时方程给出了“速度—厚度—时间”的解析关系,截距时间承载厚度信息,视速度承载速度分层信息。掌握这两者的物理意义和计算方式,才能在后续剩余静校正中正确建立初始速度模型、合理解释初至异常,并为折射层析反演提供可靠的初始约束。