2. 地震阻抗基本概念与岩石物理基础

地震阻抗反演之所以能够从地震反射资料中推断岩性,根本原因在于地下介质的波阻抗差异会在岩层分界面上产生反射。理解阻抗的物理定义、它与速度和密度之间的关系,以及不同岩石物理参数如何改变阻抗,是后续反演建模、低频补偿和岩性解释工作的共同基础。本节从定义出发,逐步建立“参数—阻抗—反射—反演”的认知链条。

一、波阻抗与声阻抗的定义

在垂直入射条件下,平面纵波在介质中传播时,介质对波传播的“阻碍”程度用波阻抗(P-wave impedance)描述,定义为介质密度与纵波速度的乘积:

Z_p = ρ × V_p

其中:Z_p —— 纵波波阻抗,单位常用 (m/s)·(g/cm³) 或 kg/(m²·s)
      ρ  —— 岩石体积密度,g/cm³
      V_p —— 纵波速度,m/s

当讨论的是一维声波方程(只考虑压缩波、忽略横波与转换波)时,同一物理量常被称为声阻抗(acoustic impedance)。在叠后地震反演中,由于输入数据基本是纵波反射信息,实际使用的就是纵波波阻抗,工程上习惯直接简称“波阻抗”或“阻抗”。需要区分的是横波阻抗 Z_s = ρ·V_s,它出现在叠前弹性反演中,用于计算纵横波速度比、泊松比等岩性敏感参数。本节讨论的阻抗若无特别说明,均指纵波波阻抗。

阻抗的量纲是密度乘速度,因此它同时受“岩石有多重”和“波跑得多快”两个因素控制。这一点非常关键:单看速度或单看密度都可能误判岩性,而阻抗把两者合成为一个可观测、可反演的物理量。

二、反射系数与阻抗差的关系

波阻抗本身不能直接被地震仪器测量,地震记录观测到的是反射波。垂直入射时,界面反射系数由上下两层阻抗决定:

R = (Z_2 − Z_1) / (Z_2 + Z_1)

Z_1 —— 上覆层阻抗;Z_2 —— 下伏层阻抗
R > 0:下伏层阻抗更高(正极性反射,波谷/波峰取决于极性约定)
R < 0:下伏层阻抗更低(负极性反射)

这说明反射振幅只对阻抗的相对变化敏感,对阻抗绝对值不敏感。例如泥岩阻抗约 6.0×10⁶ kg/(m²·s),砂岩约 7.5×10⁶ kg/(m²·s),界面反射系数约为 (7.5−6.0)/(7.5+6.0)≈0.11;而如果砂岩阻抗变为 6.6×10⁶,反射系数仅约 0.048,地震上就接近“弱反射”甚至难以识别。这正是岩性预测中“阻抗差异小则地震难分辨”的物理根源,也是反演需要引入低频模型和约束的原因。

三、速度—密度关系与岩石物理参数的影响

速度与密度并非独立变量。大量实验表明,在沉积岩中二者存在经验关系(如 Gardner 关系):

ρ ≈ a · V_p^b   (英制单位下常用 a≈0.23,b≈0.25;
                 公制单位需换算系数,示例仅示意形式)

把该关系代入阻抗定义可得 Z_p ≈ a·V_p^(b+1),说明阻抗对速度的依赖强于对密度的依赖。实际工作中更常用的是体积密度测井与声波测井直接相乘得到阻抗曲线,经验关系主要用于无密度井段或横波预测。

影响阻抗的岩石物理参数主要有以下几类:

  • 孔隙度:孔隙度增大,岩石骨架比例下降,速度与密度同时降低,阻抗显著下降。这是砂岩储层阻抗低于致密泥岩的根本原因。
  • 孔隙流体:含水砂岩替换为含气砂岩后,纵波速度大幅下降(气的高压缩性),密度也略降,阻抗明显降低,形成“亮点”型负反射;含油则变化居中。
  • 岩性与矿物:石英含量高的砂岩速度较高,黏土含量增加会降低速度与阻抗;灰岩、白云岩等碳酸盐岩阻抗通常明显高于碎屑岩。
  • 压实与埋深:随埋深增加,压实作用使孔隙度降低、速度与密度升高,阻抗整体呈增大趋势,这也是低频趋势模型的主要来源。
  • 泥质含量与饱和度:泥质增加通常降低阻抗;含水饱和度变化对阻抗的影响弱于气,因此气层识别对阻抗最敏感。
教学示例(合理假设,非实测):设某区泥岩 V_p=2800 m/s、ρ=2.35 g/cm³,则 Z≈6.58×10⁶;含水砂岩 V_p=3100 m/s、ρ=2.30 g/cm³,Z≈7.13×10⁶;同层含气后 V_p 降至 2600 m/s、ρ=2.20 g/cm³,Z≈5.72×10⁶。可见含气使砂岩由“高于泥岩”变为“低于泥岩”,反射极性反转,这是气层识别的关键判据。

四、从阻抗到反演的认知链条

把上述关系串起来,就得到反演的基本逻辑:岩石物理参数(孔隙度、流体、岩性)→ 速度与密度 → 波阻抗 → 界面反射系数 → 地震道振幅。反演则是沿这条链反向求解:由地震振幅和井约束,先求反射系数序列,再递推或优化得到阻抗,最后由阻抗结合岩石物理模板划分岩性。理解每一环的正演关系,才能判断反演结果是否可信、低频模型是否合理、岩性切片的边界是否真实。

岩石物理参数 → 波阻抗 → 反射系数 → 地震振幅(正演) 岩石物理参数 孔隙度 / 流体 岩性 / 泥质 / 埋深 速度 V_p、密度 ρ 测井曲线 / 经验关系 Gardner、流体替换 波阻抗 Z=ρ·V_p 纵波阻抗(叠后) 横波阻抗(叠前) 反射系数 R R=(Z₂−Z₁)/(Z₂+Z₁) 正/负极性判岩性 地震振幅 叠加/偏移数据 子波褶积结果 反演方向:振幅 → 反射系数 → 阻抗 → 岩性 示例:泥岩—含水砂岩界面 Z泥岩 ≈ 6.58×10⁶,Z含水砂 ≈ 7.13×10⁶ R ≈ +0.04,弱正反射 砂岩阻抗高于泥岩,地震上不易区分 示例:泥岩—含气砂岩界面 Z泥岩 ≈ 6.58×10⁶,Z含气砂 ≈ 5.72×10⁶ R ≈ −0.07,负反射(极性反转) 含气使阻抗降低,形成可识别亮点

需要强调的是,上述示例数值仅用于说明参数变化对阻抗和反射系数的影响方向与量级,实际工区必须用本区测井曲线统计建立阻抗—岩性关系,不能直接套用。反演前应先用井旁道标定极性、提取子波,并检查合成记录与井旁地震道的匹配程度,确认阻抗与反射的对应关系无误后,再进入低频模型构建与反演参数试验。

常见误区:把“高阻抗”直接等同于“好储层”或“砂岩”。在致密砂岩、灰质砂岩或深埋压实带中,高阻抗可能对应非储层;而含气砂岩阻抗反而降低。必须结合区域岩石物理交会图(如 V_p/V_s 与阻抗交会)建立本区岩性判别标准,避免单一阻抗门槛造成误判。