平滑技术改善致密气藏产量瞬态分析（RTA）中的趋势检测

《ARABIAN JOURNAL FOR SCIENCE AND ENGINEERING》：Smoothing Techniques Improve Trend Detection in Rate-Transient Analysis (RTA) of Tight Gas Reservoirs

　　

产量瞬态分析（RTA）中使用的生产数据常被噪声和操作干扰污染，掩盖诊断特征，在流态识别和储层表征中引入不确定性。虽然异常值检测（OD）方法常被应用，但它们不能完全消除异常行为，残留噪声扭曲率归一化压力（RNP）趋势。本研究评估平滑作为预处理步骤以提高RTA可靠性。研究人员使用从均质单层低渗透气藏单裂缝垂直井生产的含噪声合成RNP数据集，评估了五种平滑技术。噪声和异常值通过加性白高斯噪声（AWGN）在多个污染水平和噪声场景下引入。研究人员使用统计和工程标准评估性能，包括平方误差和（SSE）、斜率行为保留、线性流时间识别和线性流参数（LFP）估计。在评估的方法中，线性GAM表现出最佳整体性能，实现了最低的平均训练和测试SSE以及最小的训练-测试差距。该方法还更好地保留了RNP趋势的物理特征，包括线性流斜率和流态时间，这些对于可靠的工程解释至关重要。平滑结果进一步表明，改进的趋势保留导致更一致的线性流区间识别和更准确的LFP估计。与噪声数据相比，平滑数据更接近真实响应，并在代表性合成案例中减少了解释误差。这些发现表明，平滑减少了诊断图解释的主观性，并改善了有噪声生产条件下的RTA工作流程。

**论文解读：平滑技术改善致密气藏产量瞬态分析（RTA）中的趋势检测**

**研究背景与问题**

在致密气藏的产量瞬态分析（RTA）中，生产数据常受到噪声、缺失数据和操作干扰的污染，这些干扰掩盖了诊断特征，并在流态识别和储层表征中引入了不确定性。工程师通常依赖对数-对数坐标下的率归一化压力（RNP）趋势进行视觉解释，其中斜率的微小变化可区分不同流态并直接影响关键储层参数的估计。然而，噪声或稀疏数据会扭曲这些趋势，导致有偏的解释和不可靠的结果。尽管异常值检测（OD）方法被广泛使用，但无法完全消除异常数据点，残留噪声持续影响流态识别和储层特征提取。在此背景下，平滑作为预处理步骤在RTA工作流程中至关重要，但常被忽视。通过预测OD精度有限区域中的代表性数据点，平滑可保留RNP数据的潜在趋势，提高诊断图的清晰度，对识别流态时间尤为关键，因为斜率的微小偏差会导致储层属性解释的显著差异。

**研究内容与结论**

为了系统评估平滑对RTA可靠性的影响，研究人员对五种平滑技术进行了基准测试：Savitzky–Golay滤波、局部加权散点图平滑（LoWeSS）、广义加性模型（GAM）、高斯核平滑和B样条。研究使用来自Munthe和Lee [6]的合成RNP数据，该数据基于均质单层低渗透气藏单裂缝垂直井的生产模拟，包含线性流、过渡流和边界主导流（BDF）三种流态。通过加性白高斯噪声（AWGN）在25%、50%和75%的污染水平下引入噪声和异常值，每场景生成10个数据集。性能评估采用统计标准（平方误差和，SSE）和工程标准（斜率行为保留、线性流时间识别、线性流参数（LFP）估计）。最佳超参数通过70%训练集优化，30%测试集验证。结果显示，线性GAM（λ=100）在所有方法中表现最佳，具有最低的平均训练和测试SSE及最小的训练-测试差距，且能更好地保留RNP趋势的物理特征，包括线性流斜率和流态转换。平滑后，线性流区间识别更一致，LFP估计更准确。在代表性合成案例中，平滑数据相比噪声数据更接近真实响应，减少了解释误差。研究人员还将线性GAM应用于Marcellus盆地27口多裂缝水平井的现场RNP数据，结果表明平滑改善了趋势连续性，使匹配的0.5斜率区间更易识别，并提高了现场尺度下LFP估计的一致性。研究结论认为，平滑减少了诊断图解释的主观性，改善了有噪声生产条件下的RTA工作流程。

**主要关键技术方法**

本研究采用的关键技术方法包括：（1）五种平滑技术：Savitzky–Golay滤波（局部多项式回归）、局部加权散点图平滑（LoWeSS，非参数局部加权回归）、广义加性模型（GAM，基于样条函数的加性模型）、高斯核平滑（局部加权平均）和B样条（分段多项式插值）。（2）合成数据集：来自Munthe和Lee [6]，基于均质单层低渗透气藏单裂缝垂直井的模拟生产数据，提供已知的RNP真值。（3）噪声污染：采用加性白高斯噪声（AWGN）在25%、50%和75%水平下污染所有流态。（4）评估标准：统计标准（平方误差和，SSE）和工程标准（斜率行为保留、线性流时间识别、线性流参数（LFP）估计）。（5）现场应用：Marcellus盆地27口多裂缝水平井的现场RNP数据，用于验证方法的实用性。

**研究结果**

**平滑方法性能比较**
通过训练-测试SSE评估，线性GAM表现最佳，具有最低平均训练和测试SSE及最小训练-测试差距。Savitzky–Golay滤波在非多项式RNP趋势下表现不佳；高斯核平滑对核宽度敏感；B样条在噪声数据中可能传播噪声；鲁棒LoWeSS因鲁棒加权机制抑制了与流态转换相关的信号成分，导致极高SSE，不适用于本数据集。线性GAM的加性平滑结构和正则化框架平衡了灵活性和稳定性，避免了过拟合。

**线性GAM对趋势保留的影响**
在代表性合成数据集（50%噪声水平）上，线性GAM（λ=100）的平滑曲线紧密跟随真值RNP趋势，尤其在中等和高噪声条件下。相对残差分析显示，大多数偏差在数据集中被减小，特别是BDF区域，但早期稀疏数据仍存在较高残差，表明平滑在稀疏区域的局限性。

**对线性流识别的影响**
在对数-对数RNP随时间图上，自动曲线拟合0.5斜率趋势。噪声数据过度估计了线性流区间，导致线性流结束时间（t_elf）错误；平滑数据产生的拟合线更接近真实趋势，且区间保持在真实线性流范围内。这改善了t_slf（线性流开始时间）和t_elf的可靠性。

**对LFP估计的影响**
在RNP对√t_a,LS（线性叠加伪时间）图上，从原点到解释的线性流结束点拟合斜率，计算LFP。平滑数据产生的拟合趋势更接近主导响应，LFP估计值（如示例中为2.14×10⁵ ft²md^1/2）更接近参考值（2.21×10⁵），而噪声数据LFP（1.50×10⁵）误差较大。

**现场应用**
对Marcellus盆地27口井应用线性GAM（λ=100），平滑后RNP趋势更连贯，0.5斜率区间识别更稳定。代表井（Well 4和Well 7）的t_elf分别减少26天和46天，表明平滑消除了短尺度波动对区间的延长效应。27口井的P50线性流区间从33<t≤217天变为47<t≤328天，P10-P90范围变化表明平滑主要稳定了区间连续性而非系统性偏移。LFP分布P50从243,937 ft²md^1/2升至320,641 ft²md^1/2，表明平滑减少了短尺度波动对斜率估计的敏感性。

**讨论与结论**

研究人员讨论指出，平滑参数选择需结合工程判断，避免过度平滑掩盖真实物理行为（如流态转换、操作变化）。平滑假设短尺度变化不是主要关注特征，当数据包含多个局部趋势变化时，平滑可能偏离真实趋势，因此需在复杂井中结合工程审查。研究还指出本工作局限于单相气体系统，未评估多相流条件，未来应扩展至导数诊断、多相流和更复杂现场条件。
研究结论翻译如下：
1. 线性GAM在评估方法中提供最佳整体平滑性能。对于分析的合成数据集，λ=100在噪声降低、趋势保留和泛化性能之间取得最有利平衡。
2. SSE作为统计指标有用，但不应单独使用。结合工程标准，线性GAM不仅实现最低平均训练和测试SSE，还更好地保留了RNP趋势的物理意义特征，这对于可靠的RTA解释至关重要。
3. 改进的平滑导致更稳定的线性流态识别。在合成基准中，线性GAM平滑的RNP数据更紧密地跟随真实响应，并在对数-对数RNP随时间图上产生更一致的0.5斜率趋势匹配，提高了t_slf和t_elf解释的可靠性。
4. 改进的趋势保留也改善了基于斜率的储层解释。通过减少短尺度波动同时保留主要诊断趋势，线性GAM在RNP对√t_a,LS图上产生更稳定的斜率估计，因此比未平滑数据提供更可靠的LFP估计。
5. 现场数据应用展示了平滑对初步RTA的实用价值。对于Marcellus现场数据集，线性GAM改善了趋势连续性，使匹配的0.5斜率区间更易识别，并为t_slf、t_elf和LFP的现场尺度解释提供了更可靠的基础。然而，益处因井而异，当数据包含清晰主导趋势时平滑最有效。在具有多个局部趋势变化的井中，平滑响应可能偏离预期真实行为，仍需工程判断。
总体而言，线性GAM被证明是初步RTA有价值的预处理步骤，因为它减少了诊断图解释的主观性，并在有噪声条件下提高了流态识别和基于斜率的储层表征的可靠性。同时，平滑不应被视为纯粹的自动校正工具；其结果仍需结合井的物理行为和现场数据结构进行审查。未来工作应将所提出的平滑工作流程扩展至基于导数的诊断、多相流条件和更复杂的现场条件。