4. AVO属性参数提取

叠前道集经过保幅处理、动校正拉平与角度子波一致性校正之后，就进入了属性提取环节。这一步的本质，是把每个CDP道集内随入射角变化的反射振幅，用一个可解释的参数模型压缩成若干张属性剖面。最经典的做法是Aki-Richards近似下的两项线性拟合：把反射系数写成截距P与梯度G的线性组合，再通过最小二乘把每个时间采样点上的P和G反演出来。理解这一步的关键，是先明确符号约定——本文统一采用“振幅随入射角增大而减小记为G为负”的约定，即R(θ)=P+G·sin²θ，后续所有正负号判断都以此为准。

从道集到P、G：逐采样点的线性反演

实际操作中，输入是角度道集（常见为近、中、远三个或五个角度范围叠加道，如5°、15°、25°、35°）。对某一个时间采样点t，设各角度道的振幅为A(θᵢ)，则拟合方程为A(θᵢ)=P+G·sin²θᵢ。以三个角度为例，写成矩阵形式后用法方程求解：

% 示例：单采样点P/G最小二乘拟合（角度单位：度）
theta = [5 15 25 35];
x = sind(theta).^2;              % 自变量 sin^2(theta)
A = [0.182 0.151 0.108 0.052];   % 示例振幅，随角度增大而减小
M = [ones(numel(x),1), x(:)];    % 设计矩阵 [1, sin^2(theta)]
pg = M \ A(:);                   % pg(1)=P, pg(2)=G
P = pg(1); G = pg(2);
fprintf('P=%.4f  G=%.4f\n', P, G);   % 示例输出 P=0.1963  G=-0.1176

上面是一组教学示例数据，用于说明计算路径而非实测结果。可以看到P约0.196、G约-0.118，属于典型的“截距为正、梯度为负”组合。工程实现时并不会逐点手算，而是把上述最小二乘展开成解析式，在整条剖面所有采样点上向量化批量求解，同时输出拟合相关系数R²作为质量控制道，R²低的区域说明振幅随角度变化不满足线性假设，P、G解释需谨慎。

各属性的物理含义与派生计算

截距P对应零入射角反射系数，近似反映波阻抗差，其地质意义与常规叠后振幅相近，主要受岩性与孔隙度控制。梯度G描述振幅随入射角的变化率，对横波信息与流体敏感，是识别含气砂岩“亮点/暗点”的核心参数。两者组合后派生出若干常用属性：

• 远道叠加：直接对远角度范围（如25°~35°）道集做叠加，突出大角度响应，对流体引起的振幅异常放大，但抗噪性弱于P、G拟合。
• P×G（截距乘梯度）：突出P与G异号的含气异常，是经典的烃类指示因子。
• P−G：对背景趋势更稳健，常用于压制压实趋势后的异常增强。
• 流体因子：常用形式如F=G·P或基于Russell近似的流体项f=ρ·(Vp²−2Vs²)，其物理出发点是剔除剪切模量贡献、保留体积模量相关的流体响应。

计算流程与质量控制

完整的属性提取链条可概括为：角度道集准备→逐点线性拟合→P/G剖面输出→派生属性计算→交会分析与异常筛选。其中每一步都有可观察的质量指标：角度道集需检查远道是否拉平、有无剩余时差；拟合阶段关注R²分布与P、G的横向连续性；派生属性阶段则通过P-G交会图判断背景趋势与异常偏离方向。下图给出了这一流程中各属性的输入输出关系与判别逻辑。

需要强调的是，P、G本身是拟合量而非独立测量量，其可靠性高度依赖道集质量。若远道存在剩余时差或动校拉伸，sin²θ的线性关系被破坏，G会被系统性低估或高估，进而污染流体因子。因此规范做法是先做角度道集剩余时差校正与道集均衡，再进入拟合；拟合后必须检查R²剖面，对低相关区回退排查。此外，P、G对角度范围敏感，近道缺失会使P不稳定，远道缺失会使G不稳定，提取前应确认有效角度覆盖范围与设计一致。

从成果交付角度看，一次规范的属性提取应至少产出：P剖面、G剖面、R²质控剖面、远道叠加剖面、流体因子剖面，以及P-G交会图与异常区平面分布图。这些成果共同构成后续AVO正演标定与流体响应解释的输入，也是判断属性提取是否成功的直接依据——若P、G横向连续、R²整体较高、交会图上异常点明显偏离背景趋势，则说明提取流程有效，可进入下一阶段的定量解释。