4. 调谐曲线与调谐厚度图版
调谐曲线是把“厚度”这一地质变量与“地震响应”这一观测变量联系起来的定量桥梁。对薄层而言,厚度变化不会改变顶底反射的到达时间差到足以分辨的程度,却会系统性地改变复合波形的振幅、频率成分与相位。调谐图版则是把这些曲线按不同子波主频、不同层数、不同层间组合方式整理成可查阅的图表,使解释人员能从观测到的振幅、频率、相位特征反推薄层厚度或厚度组合。本节围绕单层与多层两类模型,说明曲线怎么算、图版怎么建、规律怎么读。
一、单层调谐曲线的建立
单层模型是最基本的调谐单元:一套厚度为 h 的薄层夹在上下厚层之间,顶底界面反射系数分别为 r₁、r₂。设子波为雷克子波,主频 f₀,则顶底反射的到达时差为 Δt = 2h/v(v 为层速度,双程时)。合成响应可写成两个时间延迟子波的叠加:
s(t) = r₁·w(t) + r₂·w(t − Δt)
其中 w(t) 为雷克子波,Δt = 2h/v
实际计算时按固定采样间隔(如 1 ms)离散化,对 h 从 0 到调谐厚度以上(例如 0~λ/2,λ = v/f₀)逐点扫描,每个厚度点合成一道,再提取三个属性:复合波峰振幅 A(h)、主频或视主频 F(h)、以及顶底复合波的相位 φ(h)。以 f₀ = 30 Hz、v = 3000 m/s 为例,λ ≈ 100 m,调谐厚度约为 λ/4 ≈ 25 m,扫描范围取 0~50 m 即可覆盖调谐峰与衰减段。这些参数属于教学示例,实际工区需用本区子波与层速度替换。
单层曲线有三个典型特征。振幅曲线在 h = λ/4 附近出现极大值,这就是调谐峰;当顶底反射系数同号(如砂泥岩阻抗差同向)时振幅随厚度先增后减,当异号时则出现极性反转式的极小值。视主频在薄层段随厚度减小而升高,因为顶底子波相互干涉使复合波变“窄”。相位在调谐厚度附近变化最快,是厚度最敏感的量,但也是最容易被噪声和子波相位误差污染的量。因此读图版时应以振幅为主、频率为辅、相位作精细校正。
二、多层薄互层的调谐曲线
多层模型把若干薄层按一定厚度组合与层间间隔叠置,其响应不再是两个子波的简单叠加,而是多个反射界面子波的干涉叠加。设 n 层,第 i 个界面反射系数 rᵢ、到达时间 tᵢ,则:
s(t) = Σ rᵢ·w(t − tᵢ),i = 1…n
tᵢ 由各层厚度与层速度累加得到
多层曲线与单层的关键差别在于:调谐峰位置会随总厚度与层间间隔的组合发生漂移,振幅极大值不再唯一对应 λ/4;当层间间隔接近子波半周期时,会出现“超调谐”增强,振幅可高于同总厚度的单层;当层间间隔使子波反相叠加时,则出现振幅压制。频率特征也更复杂,视主频曲线可能出现双峰或平台。因此多层图版不能只给一条曲线,而要按“层数 × 层间间隔 × 单层厚度”三个维度组织成曲线族。
三、调谐图版的制作与组织
调谐图版不是一张图,而是一套按参数分层的曲线族。制作流程可归纳为四步。第一步确定图版参数域:子波主频取工区实际频带范围(如 20、25、30、35 Hz 四档),层速度按岩性组合给定,反射系数按目标砂泥岩阻抗差设定正负号。第二步对每个参数组合扫描厚度,计算 A、F、φ 三条曲线,并记录调谐峰厚度 h_peak 与峰值振幅 A_peak。第三步把同主频、同层数的曲线画在一张图上,横轴厚度、纵轴归一化振幅,形成单层图版与多层图版两套。第四步补充一张“峰值厚度—主频”关系表,便于快速查表。
| 图版分册 | 参数维度 | 主要用途 |
|---|---|---|
| 单层振幅图版 | 主频 20/25/30/35 Hz,r 同号与异号 | 由振幅估单层厚度区间 |
| 单层相位图版 | 同上,标注敏感区 | 调谐厚度附近精细校正 |
| 多层振幅图版 | 层数 2/3/5,层间间隔 5/10/15 m | 识别超调谐增强与压制 |
| 峰值厚度查算表 | 主频—h_peak—A_peak | 现场快速查表反演 |
四、读图版时的判断依据
使用图版要建立“先区间、后精细”的顺序。先用振幅落在图版曲线的哪一段判断厚度区间:若振幅处于上升段,厚度小于调谐厚度,振幅与厚度正相关;若处于下降段,厚度大于调谐厚度,振幅与厚度负相关,此时单靠振幅无法区分两个厚度值,必须借助相位或频率。再用视主频判断薄层化程度:厚度越小,视主频越高。最后在调谐厚度附近用相位敏感区做校正,因为该区相位随厚度变化率最大,微小厚度差会带来明显相位差。
把单层与多层两套曲线放在同一坐标系中对照,可以直观看到多层干涉造成的峰值漂移与增强,这正是图版分册存在的意义。实际解释时,先由振幅图版锁定厚度区间,再由相位图版在敏感区内细化,最后用井震标定检验反演厚度是否落在合理范围,形成“图版查算—井点校验—横向外推”的闭环。