3. Zoeppritz方程与反射系数

叠前AVO分析的理论根基是弹性界面上的反射透射问题。当一束P波以入射角θ₁斜入射到两层介质的弹性界面时，界面两侧会同时产生反射P波、反射SV波、透射P波和透射SV波四个波场。Zoeppritz方程正是描述这四个波场振幅关系的严格解，它由界面处的位移连续和应力连续条件（法向位移、切向位移、法向应力、切向应力四个边界条件）联立得到。理解这个推导思路，比记住方程本身更重要：它告诉我们反射系数本质上是“边界条件约束下的能量分配结果”，因此必然同时依赖上下介质的纵波速度Vp、横波速度Vs和密度ρ，以及入射角。

Zoeppritz方程是精确的，但它的解析形式复杂，四个反射/透射系数耦合在一起，难以直观看出各弹性参数对振幅的影响。工程上真正使用的是它的近似式。近似推导的核心思路是：假设界面两侧弹性参数差异较小（弱反差），对速度比、密度比做一阶泰勒展开，并利用Snell定律把透射角用入射角表示，从而把耦合方程组解耦为显式的线性组合。这一思路直接决定了近似公式的适用前提——界面反射率不能太大，通常要求上下层参数差异在20%~30%以内，且入射角不宜过大。

Aki-Richards近似：参数分解的起点

Aki-Richards近似把反射系数写成Vp、Vs、ρ三个参数相对变化量的加权和：

R(θ) ≈ (1/2)(1 - 4(Vs/Vp)²sin²θ)·(ΔVp/Vp)
      + (1/2)(1/cos²θ)·(Δρ/ρ)
      - 4(Vs/Vp)²sin²θ·(ΔVs/Vs)

它的物理意义非常清晰：第一项是纵波速度贡献，第二项是密度贡献，第三项是横波速度贡献，权重都随入射角变化。在零入射角时，后两项的sin²θ项消失，退化为垂直入射的声波反射公式R₀=(ΔVp/Vp+Δρ/ρ)/2，这是检验近似式正确性的基本自洽条件。Aki-Richards式的优点是参数物理意义明确、对中到大角度仍有一定精度，缺点是三项之间存在相关性，直接反演Vp、Vs、ρ容易病态。

Shuey近似：面向岩性解释的重组

Shuey在Aki-Richards基础上做了进一步重组，把反射系数表达为截距A、梯度B和曲率C的形式：

R(θ) ≈ A + B·sin²θ + C·(tan²θ - sin²θ)
A = R₀ = (1/2)(ΔVp/Vp + Δρ/ρ)
B = (1/2)(ΔVp/Vp) - 2(Vs/Vp)²[(Δρ/ρ) + 2(ΔVs/Vs)]
C = (1/2)(ΔVp/Vp)

Shuey式的价值在于把AVO响应拆成“截距+梯度”两个可解释的量：截距A反映零偏移距反射强度，梯度B反映振幅随偏移距的变化速率。在中小入射角（一般小于30°）范围内，曲率项C贡献很小，可近似忽略，于是R(θ)≈A+B·sin²θ，这就是业界最常用的“A-B”线性拟合形式。实际处理中，把CDP道集按sin²θ分角度段做线性回归，斜率即梯度B，截距即A，再投到A-B交会图上划分AVO类别。

适用条件与精度对比

三种表达式的取舍可归纳为一张对照表：

需要强调的是，近似式的误差不仅来自角度，还来自Vs/Vp比值。当Vs/Vp偏离0.5较多（如气层含气后Vs/Vp明显降低）时，Aki-Richards中横波项的权重变化显著，此时用固定Vs/Vp假设去拟合梯度会产生偏差。因此实际解释时应先从测井数据统计目的层的Vs/Vp背景值，再决定近似式的参数化方式。

从图中可以直观看到：三条曲线在零入射角处重合（自洽性检验通过），随角度增大逐渐分离，Shuey两项式在30°以后与精确解的偏离明显加大，而Aki-Richards式因保留了更完整的参数项，在中大角度仍较贴近精确解。这正是实际AVO反演中“小角度用Shuey、大角度或强反差用Aki-Richards、正演校验用Zoeppritz”这一分工的物理依据。

落到实操层面，建议按以下顺序建立对反射系数的定量认识：先用测井曲线计算目的层上下方的Vp、Vs、ρ，代入Zoeppritz精确解得到全角度反射系数曲线作为“真值”；再用同一组参数分别计算Aki-Richards和Shuey两项式，比较三者在0°~40°的偏差；最后根据偏差大小确定本工区AVO反演应采用的角度范围和近似式。这一“精确解定标—近似式选型”的流程，是把方程理论转化为可执行解释方案的关键一步。