📝 摘要
优化降压路径对于在提高水合物开采效率的同时避免潜在的地质灾害至关重要。本研究在柔性塑料容器中生成甲烷水合物，然后在恒定围压和恒定围温条件下进行产气实验。应用CMG STARS模拟器匹配实验产气行为，并推导水合物本征分解常数。其次，分析了不同降压路径下的流体产出行为、压力-温度响应以及水合物饱和度演化规律。结果表明，综合气水比（IGWR）随每一步降压幅度的增加而线性减小，而随降压步数的增加而对数增加。在相同的初始平均水合物饱和度和相同的总压降幅度下，缓慢、多级降压策略有助于提高IGWR并避免严重的温度下降。一旦产气暂停，孔隙压力呈对数回升，而关井操作结束后，孔隙压力会瞬时降至常规水平。我们推测，关井操作在长期尺度上几乎不影响IGWR和地层压力-温度响应。

🖥️ CMG软件应用情况总结

项目 内容
所用软件 CMG STARS（先进过程与热采模拟器）
模型类型 热采组分模拟器，考虑流体流动、传热和固体（水合物）分解的耦合效应
模拟对象 含水合物沉积物（HBS）中甲烷水合物的降压法分解产气过程
应用方式 • 历史拟合：使用STARS匹配室内实验的气体产量和压力演化数据，确定水合物本征分解常数（KD）
• 敏感性分析：基于验证的数值模型，研究不同降压路径（降压步数、每步降压幅度、降压速率、关井操作）对产气产水行为、地层温度和水合物饱和度演化的影响
水合物分解模型 基于Kim-Bishnoi模型，Uddin等人修改的动力学方程（式1-2），水合物被视为填充在孔隙中的“不可移动重油”相
几何模型 • 形状：三维圆柱形（模拟橡胶套筒）
• 尺寸：Φ60 mm × 1000 mm
• 网格划分：I方向20个网格（50 mm/格），R方向25个网格（3 mm/格），θ方向20个网格（18°/格）
• 纵向剖面：沿I方向取剖面 a-a 分析内部物理响应
• 边界条件：恒定围压（6 MPa 实验，15 MPa 现场模拟）、恒定围温（7°C 实验，17°C 现场模拟），水合物与围压油之间无质量交换
初始条件 • 实验模拟：初始孔隙压力4 MPa，初始温度3.5°C
• 现场模拟（Nankai Trough）：初始孔隙压力13 MPa，初始温度13.5°C
• 初始水合物饱和度：~75%、~50%、~10%三种情况
• 初始水饱和度：相应计算
岩石/流体参数（表1） • 孔隙度：39.9%
• 绝对渗透率：606 mD
• 水合物密度：919.7 kg/m³
• 水合数：5.75
• 沉积物干密度：2650 kg/m³
• 热导率：水合物0.392，孔隙水0.58，沉积物5 J/(m·s·K)
相对渗透率 Stone方程（式3），SirA=0.25，SirG=0.01，n=3.5，nG=2.5
降压方案（图5） • 总压降固定：4 MPa（实验）或从13 MPa降至9 MPa（现场模拟，水合物-水-气平衡压力~10 MPa）
• 路径1：一步快速降压
• 路径2-4：相同平均降压速率（4 MPa/3h），步数分别为2、4、16步
• 路径5：慢速降压（4 MPa/4h）
• 路径6：带关井操作的降压
评价指标 累计产气量、累计产水量、综合气水比（IGWR）、温度变化、水合物饱和度变化、水合物完全分解时间

文中明确指出：

“The CMG-STARS simulator was applied to match the experimental gas production behavior and to derive the hydrate intrinsic dissociation constant.” （摘要部分）
“CMG-STARS (Computer Modeling Group Ltd., Canada) [39] are widely adopted to simulate hydrate production processes.” （第3节）

🧪 模拟方案与主要结果

1. 实验与模型验证（第2-3节）

  • 实验：在柔性橡胶套筒（Φ60×1000 mm）中生成甲烷水合物（初始饱和度~30%），进行一步快速降压和两步降压产气实验。

  • 历史拟合：CMG STARS成功匹配了两种降压模式的产气速率、累计产气量和压力演化（图3-4）。

  • 校准参数：水合物本征分解常数 KD=2.43×104 mol/(m2⋅Pa)。

2. 降压路径对综合气水比（IGWR）的影响（第4.2节）

  • IGWR定义：累计产气量 / 累计产水量（式4）。

  • 关键发现：

    • IGWR随初始水合物饱和度增加而显著增加（图8）。

    • 降压幅度效应：IGWR随每一步降压幅度的增加而线性减小（图7a）。

    • 降压步数效应：IGWR随降压步数的增加呈对数增加，步数超过一定值后趋于平稳（图7b）。

    • 结论：多级降压对提高IGWR重要，但无需最大化步数。

3. 降压路径对温度演化的影响（第4.3节）

  • 低温区传播：低温区从出口（生产井）开始，在I方向非活塞式向远端传播；在R方向呈轴对称演化（图9）。

  • 关键发现：

    • 降低平均降压速率显著缓解温度下降（图10a, 10c, 10e）。

    • 固定平均降压速率下，增加步数（减少每步降压幅度）对缓解温度下降作用有限（图10b, 10d, 10f）。

    • 结论：缓慢、多级降压策略是避免温度下降和水合物再生的最佳选择。

4. 水合物饱和度演化与分解时间（第4.4节）

  • 分解前缘：从模型边缘靠近出口处开始，非活塞式向远端扩展，R方向轴对称（图11）。

  • 完全分解时间：

    • 随平均降压速率降低而显著增加（路径1,4,5对比，图12）。

    • 固定平均降压速率下，分解时间几乎不受降压步数影响（路径2,3,4对比，图12）。

    • 关井操作（路径6）会显著延长水合物完全分解时间。

5. 关井操作的影响（第4.5节）

  • 压力响应：关井期间，孔隙压力呈对数回升；重新开井后，压力瞬时降至常规水平（图13a）。

  • 温度响应：关井期间温度回升缓慢；重新开井后温度呈对数下降（图13b）。

  • 关键结论：关井操作在长期尺度上几乎不影响IGWR和地层P-T响应，但会显著延长水合物分解持续时间。

✅ 主要结论

  1. CMG STARS有效性：CMG STARS能够有效模拟水合物降压分解过程，通过实验数据历史拟合可准确确定水合物本征分解常数。

  2. 多级降压优势：在固定总压降下，多级降压可显著提高综合气水比（IGWR），但步数超过最优值后增益有限。

  3. 温度控制策略：降低平均降压速率是缓解开采过程中温度下降和防止水合物再生的最有效手段，单纯增加降压步数效果有限。缓慢、多级降压策略是最优选择。

  4. 关井影响有限：关井操作会暂时引起孔隙压力回升和分解时间延长，但在长期尺度上对累计产气和地层P-T响应影响不大。

🏛️ 作者及单位信息

  • 作者: Yanlong Li, Chuqiao He, Nengyou Wu, Qiang Chen, Changling Liu, Zhixue Sun, Yurong Jin, Qingguo Meng

  • 单位:

    1. 中国地质调查局青岛海洋地质研究所，自然资源部天然气水合物重点实验室

    2. 青岛海洋科学与技术国家实验室，海洋矿产资源实验室

    3. 中国石油大学（华东），石油工程学院

    4. 里贾纳大学（University of Regina），工程与应用科学学院，石油系统工程系

  • 国家重点研发计划（No. 2017YFC0307600）

💡 补充说明
该研究是天然气水合物开采数值模拟领域的典型应用案例。研究的核心创新点在于：

  1. 模型校准：使用CMG STARS对室内水合物降压分解实验进行了严格的历史拟合，确定了水合物本征分解常数这一关键动力学参数。

  2. 操作参数优化：系统研究了降压路径（步数、幅度、速率）对IGWR和温度响应的影响，提出了“缓慢、多级降压”的最优策略。

  3. 关井行为分析：首次定量分析了关井-开井操作对孔隙压力、温度和水合物分解时间的动态影响，为现场生产管理提供了理论依据。

该研究对从事天然气水合物开采工艺优化、水合物藏数值模拟以及相关实验研究的工程师和科研人员具有重要的参考价值。

case103

发表评论