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₁²))。这就是单层折射法求风化层厚度的基本公式,也是折射静校正最核心的正演—反演关系。
二、多层折射波与视速度的物理意义
当存在多个速度递增界面(如风化层—低速层—高速层),每一层折射波在远炮点段都会形成一条近似直线。第 n 层折射波在时距曲线上的斜率倒数为该层的“视速度”,它等于该层真实速度(水平层状介质下),而截距时间则累积了上覆各层的厚度贡献。以三层为例,第二界面折射波的走时可写成:
t(x) = x / v₃ + t₀₂
t₀₂ = 2h₁·√(v₃² − v₁²)/(v₁·v₃) + 2h₂·√(v₃² − v₂²)/(v₂·v₃)
可见截距时间不是单一厚度信息,而是各层厚度与速度对比的综合反映。这正是多层折射反演需要“逐层剥离”的原因:先用第一段直线求 v₂ 和 h₁,再从第二段截距中扣除第一层贡献,才能解出 h₂。
视速度的物理意义需要强调两点。第一,它由时距曲线远炮点段的斜率倒数得到,反映的是折射界面以下高速层的速度,而不是上覆层速度。第二,在倾斜界面或横向变速情况下,视速度会偏离真实速度,此时需要采用相遇观测系统或互换法来校正。实际处理中,常把初至时间按炮检距分段线性拟合,用各段斜率判断速度分层,用截距时间估算厚度,这就是“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,与模型吻合。这个“正演—拾取—反演”的闭环,是检验静校正参数是否合理的直接手段。
总结来说,单层与多层折射波走时方程给出了“速度—厚度—时间”的解析关系,截距时间承载厚度信息,视速度承载速度分层信息。掌握这两者的物理意义和计算方式,才能在后续剩余静校正中正确建立初始速度模型、合理解释初至异常,并为折射层析反演提供可靠的初始约束。