3. 含油气性识别的岩石物理基础

振幅属性之所以能指示油气，根本原因在于孔隙流体改变了岩石的弹性参数，进而改变了反射振幅随偏移距（入射角）的变化规律。要把“亮点”“暗点”“AVO异常”从经验现象上升为可解释、可标定的判据，必须打通三个环节：用流体替换把饱和度换算成弹性参数，用AVO正演把弹性参数换算成振幅响应，再用敏感参数分析找出区分油气与围岩、盐水的最有效组合。本节围绕这条理论链条展开，并给出可操作的示例计算流程。

一、流体替换：从饱和度到弹性参数

流体替换的核心是Gassmann方程。它描述饱和岩石体积模量 K_sat 与骨架（干岩石）模量 K_dry、矿物模量 K₀、孔隙流体模量 K_fl 及孔隙度 φ 的关系：

K_sat = K_dry + (1 - K_dry/K_0)^2 /
        [ φ/K_fl + (1-φ)/K_0 - K_dry/K_0^2 ]

μ_sat = μ_dry        （剪切模量不受孔隙流体影响）

1/K_fl = Sw/K_w + (1-Sw)/K_hc   （Reuss 均匀饱和）

实际操作顺序是“先卸后充”：已知含水饱和岩石的 K_sat、μ 与 φ，先用Gassmann反解得到 K_dry，再把 K_fl 替换为目标流体（如含油饱和度 S_o=60%）重新正算。密度同步按 ρ = (1-φ)ρ₀ + φ[S_wρ_w + (1-S_w)ρ_hc] 更新，纵波速度由 V_p=√[(K_sat+4μ/3)/ρ] 得到。

示例（合理假设的砂岩参数，非实测）：φ=0.25，K₀=37 GPa，K_w=2.2 GPa，K_hc=0.5 GPa，ρ₀=2.65、ρ_w=1.03、ρ_hc=0.7 g/cm³，含水时 V_p=3050 m/s、V_s=1600 m/s。替换为 S_o=60% 后，K_fl 由 2.2 GPa 降至约 0.86 GPa，V_p 降至约 2870 m/s，V_s 基本不变，密度降至约 2.30 g/cm³。可见流体替换主要压低纵波速度与密度，这正是后续AVO异常的物理来源。

二、AVO响应：振幅随入射角的变化

采用Aki-Richards近似，反射系数随入射角 θ 变化：

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

工程上更常用截距 P 与梯度 G 的线性形式 R(θ)=P+G sin²θ。含气砂岩通常表现为“第三类AVO”：含水时近道为负反射，含气后近道更负、远道绝对值减小甚至反号，即 P 负、G 正；含油砂岩变化幅度小得多，多为 P 负、G 近零。因此单看近道振幅无法区分油与气，必须结合梯度与远道信息。

三、敏感参数分析：找出最有效的判别量

敏感参数分析的目标是从众多弹性参数中挑出对流体变化最敏感、对岩性变化相对不敏感的量。常用做法是计算流体替换前后各参数的相对变化率，并构造交会图检验可分性。常用敏感参数包括：

• V_p/V_s 比值：含气时显著下降，是最经典的含气指示；
• 泊松比 σ：含气砂岩可低至0.1以下，致密含水砂岩约0.25～0.3；
• λρ、μρ（拉梅参数与密度乘积）：λρ 对流体敏感、μρ 对岩性敏感，二者交会可分离“岩性效应”与“流体效应”；
• 流体因子 f = V_p²ρ − c·V_s²ρ（c 常取2～2.5），用于压制背景压实趋势。

示例判断：上述砂岩含水时 V_p/V_s≈1.91，含油60%后≈1.79，降幅约6%；若替换为含气80%，可降至约1.55，降幅近19%。这说明 V_p/V_s 对气敏感、对油中等敏感，实际解释中应把“含气”与“含油”分级判别，而不是笼统称为“油气异常”。

四、把理论落到解释流程上

实际工作中，建议按“标定—正演—反演—校验”的顺序推进。第一步，用研究区井的测井曲线（V_p、V_s、ρ、φ、S_w）做流体替换，生成含水、含油、含气三套弹性曲线，作为正演的输入与解释的参照。第二步，用Aki-Richards公式正演目的层顶底界面的 P、G，绘制 R-θ 曲线，确认目标层属于哪一类AVO。第三步，对叠前道集做AVO反演得到 P、G 数据体，并进一步提取 V_p/V_s、λρ、μρ 等属性。第四步，用已知井做交会标定，确定油、气、水各自的参数门槛，再外推到无井区。

综上，流体替换回答了“饱和度如何改变弹性参数”，AVO响应回答了“弹性参数如何改变振幅随角度的分布”，敏感参数分析则回答了“用哪个量区分油、气、水最有效”。三者构成一条可计算、可标定、可验证的理论链，是把振幅异常转化为含油气性解释结论的岩石物理基础。