4. T2谱反演方法与页岩油信号处理

核磁共振T2谱并不是仪器直接输出的“谱”，而是从Carr-Purcell-Meiboom-Gill（CPMG）回波串经过数学反演得到的解。页岩油储层孔隙细小、含油饱和度与含氢指数关系复杂、回波串信噪比低，反演环节的任何选择都会直接改变T2谱的形态、峰位和幅度，进而影响可动流体饱和度的计算。因此，理解反演算法、正则化参数、信噪比与谱图可靠性之间的关系，是页岩油核磁评价中不可跳过的一步。

一、反演问题的本质：病态线性反演

CPMG测量得到的是离散回波幅度序列 M(t_i)，它与T2分布 f(T₂) 之间满足第一类Fredholm积分方程：

M(t_i) = Σ_j K(t_i, T2_j) · f(T2_j) · ΔlogT2 + ε_i
K(t_i, T2_j) = exp(-t_i / T2_j)

写成矩阵形式即 M = K·f + ε。核矩阵K的列之间高度相关，导致该反演问题是严重病态的：微小的噪声扰动会被放大成解的巨大振荡。因此不能直接求逆，必须引入正则化约束，把“拟合残差最小”和“解的平滑性/非负性”放在一起权衡。这正是T2谱反演的核心矛盾。

二、反演算法选择

页岩油分析中常见的反演实现有三类，选择依据主要是信噪比水平、T2分布是否多峰以及计算效率要求。

• 非负最小二乘（NNLS）+ Tikhonov正则化：最经典、最稳健的方案。目标函数为 ‖Kf−M‖² + α‖Lf‖²，约束 f≥0。L通常取一阶或二阶差分算子，使解平滑。适合信噪比中等、需要稳定单峰/双峰识别的常规页岩样品。
• 截断奇异值分解（TSVD）：通过丢弃小奇异值对应的分量抑制噪声放大，物理意义直观，但截断阶数需要人工判断，多峰样品容易丢失小峰。
• 迭代类方法（如BRD、最大熵、Landweber）：最大熵法倾向于给出最平滑且信息量最大的解，适合低信噪比、宽分布页岩油信号；但迭代次数与收敛判据需要控制，否则出现过拟合或欠拟合。

实操建议：对页岩油样品，优先采用NNLS+Tikhonov作为基准反演，再用最大熵法做交叉验证；若两种方法给出的主峰位置差异超过一个T2布点间隔，应回到原始回波串检查质量。

三、正则化参数α的影响

正则化参数α是反演中最敏感的一个旋钮，它决定了“拟合数据”与“保持平滑”之间的平衡。α过大，解被过度平滑，多个相邻孔隙组分被合并成一个宽峰，小孔与大孔的界限消失，可动流体被低估或高估；α过小，解跟随噪声剧烈振荡，出现虚假小峰和负值（被非负约束截断后形成锯齿状谱）。

工程上常用的定α方法是L曲线法：以‖Kf−M‖为横轴、‖Lf‖为纵轴，扫描一系列α得到一条L形曲线，拐点处即为兼顾残差与平滑的最优α。另一种是广义交叉验证（GCV），适合噪声统计特性未知的情况。对页岩油低信噪比数据，建议在L曲线拐点附近再取2～3个α值做敏感性对比，观察主峰位置是否稳定；若主峰随α变化漂移明显，说明该样品反演结果可靠性不足，应谨慎解释。

四、信噪比控制

页岩油样品孔隙度低、含油信号弱，回波串末段往往已接近噪声底。信噪比不足会直接导致反演谱出现虚假峰或峰位偏移。控制手段主要有三条：一是增加叠加次数NS，噪声按√NS衰减，但采集时间线性增加，需在测井效率与谱质量间折中；二是合理设置回波间隔TE，页岩小孔T2短，TE过大会丢失短T2组分，一般页岩油分析取TE=0.2～0.6 ms；三是保证足够等待时间TW使流体充分极化，避免T1未恢复造成的幅度低估。

反演前应先用回波串末段（信号已衰减至噪声区）估计噪声标准差σ，再计算整体SNR。经验上，SNR低于20时反演谱的多峰结构可信度明显下降，此时应优先增大NS或降低TE，而不是靠调α“凑”出好看的谱形。

五、谱图可靠性评估

反演完成后不能直接拿谱图算可动流体，必须做可靠性检查。可操作的判据包括：拟合残差是否与噪声水平相当且无系统性偏差；谱值是否出现被非负约束截断的锯齿振荡；改变α时主峰位置漂移是否小于一个T2布点；不同算法（NNLS与最大熵）给出的主峰位置与幅度是否基本一致；以及T2截止值、谱面积与孔隙度、渗透率趋势是否物理自洽。只有这些检查通过，T2谱才可用于后续可动流体饱和度与束缚水饱和度的计算。

综上，T2谱反演是一个“算法—正则化—信噪比—可靠性”联动的过程。算法决定解的形式，α决定解的平滑程度，信噪比决定解的可辨识上限，而可靠性评估则是把数学解转化为可信地质信息的最后一道关口。在页岩油可动流体评价中，只有把这四者作为一个整体来控制，T2谱才真正具备定量解释价值。