📝 摘要
甲烷水合物被视为下一代天然气资源。降压法是开采水合物的经济有效方法，但会导致水合物储层压实和海床沉降，引发环境与机械稳定性问题。本研究提出了一种使用水平井注热水的开采方法，旨在降低降压幅度、减少环境冲击。研究构建了结合水合物分解与地层固结模型的数值模型，使用CMG STARS热采模拟器模拟热水注入法开采过程中的海床沉降。固结模型基于室内实验数据（Sakamoto等，2009、2010）进行历史拟合，建立了水合物饱和度与弹性模量的函数关系。矿场尺度模拟结果表明：热水注入法的累计产气量是降压法的1.6-12.3倍（井底压力3-9 MPa条件下）；在相同累计产气量下，热水注入法（注热水85°C，500 m³/天，井底压力9 MPa，压降4 MPa）的最大海床沉降约0.4 m，而降压法（井底压力3 MPa，压降10 MPa）的沉降约2 m。热水注入法可在不降低产气速率的同时减少海床沉降。
🖥️ CMG软件应用情况总结
| 项目 | 内容 |
|---|---|
| 所用软件 | CMG STARS（2015版本），先进过程与热采模拟器 |
| 模型类型 | 多相流体流动与温度分布耦合模型，结合水合物分解动力学和地层固结力学 |
| 模拟对象 | 海洋甲烷水合物储层，目标区域位于日本东部南海海槽（Eastern Nankai Trough） |
| 应用方式 | • 水合物相定义：将甲烷水合物定义为固相（Solid phase） • 水合物分解模型：基于水合物形成-分解平衡曲线和Arrhenius方程计算相变速率 • 绝对渗透率模型：采用幂函数描述水合物饱和度对渗透率的降低影响（Masuda模型） • 弹性模量模型：建立水合物饱和度与弹性模量的线性函数 E = E₀ + β·S_MH，基于Masui等（2005）三轴压缩实验数据 • 孔隙度模型：考虑压缩性和水合物饱和度的有效孔隙度 φ_e = φ_i·exp[κ(p-p_i)]·(1-S_MH) • 渗透率-孔隙度关系：基于Kozeny-Carman方程 • 相对渗透率模型：采用Sakamoto等（2009、2010）提出的水合物岩心实验关联式 • 海床沉降模型：通过积分各网格垂向压缩量计算沉降 Δh = ∫α(φ_i-φ_e)dz • 热水注入法模拟：采用双水平井系统（类似SAGD工艺），上井生产、下井注入85°C热水 |
| 岩心尺度模型验证 | • 模拟Sakamoto等（2009、2010）的砂填充岩心实验 • 网格：圆柱坐标系，31（径向）×52（轴向）=1,612个网格 • 初始压力：10 MPa，降压至3.3 MPa • 对比参数：累计产气量、产水量、位移量、温度分布 • 历史匹配结果：使用E=200+700·S_MH模型时，位移曲线与实验最匹配（图3） • 网格敏感性：半倍/双倍网格下产气与位移差异<0.5% |
| 矿场尺度模型参数 | • 储层尺寸：500 m × 50 m × 15 m厚，体积375,000 m³ • 初始压力：13 MPa，初始温度：12.85°C • 孔隙度：40%，绝对渗透率：1000 mD • 初始水合物饱和度：60%，初始水饱和度：40% • 双水平井：长500 m，垂直间距5 m，上井距储层顶部2 m • 热水注入温度：85°C，注入速率：500 m³/天 • 生产井井底压力：3、6、9 MPa（三种方案对比） • 前期准备：先降压90天使两井连通 • 网格：19（垂向）×15（水平）×35（沿井向）=9,975个网格 |
| 模拟对比方案 | • 降压法：单口垂直井或双水平井（持续降压） • 热水注入法：双水平井，下井注85°C热水，上井生产 • 井底压力：3、6、9 MPa三种条件对比 |
| 主要结论 | • 热水注入法累计产气量是降压法的1.6-12.3倍（3-9 MPa条件下） • 海床沉降主要受储层压力（即井底压力）控制 • 相同产气量下，热水注入法沉降约0.4 m，降压法沉降约2 m • 热水注入法可在不降低产气速率的同时减少海床沉降 |
文中明确指出（Case108.pdf）：
“In this study, a numerical model combining models of MH dissociation and consolidation has been presented to simulate seabed subsidence with gas production from a MH reservoir by hot water injection with a pair of horizontal wells using the thermal simulator CMG STARS™ (2015 version).” （第1章，第2页）
“The reservoir simulator STARS was used for numerical simulation of the gas production and consolidation based on MH dissociation and the elasticity function of MH saturation. In the simulations, MH was defined as a solid phase…” （第2.2节，第3页）
“The numerical simulations were carried out by the STARS™ using 19 blocks in the vertical direction, 15 blocks in the horizontal direction, and 35 blocks in the longitudinal direction along the horizontal wells (total of 9,975 blocks).” （第3.1节，第10页）
🧪 模拟方案与主要结果
1. 岩心尺度模型验证（第2.5节，图3）
实验数据来源：Sakamoto等（2009、2010）砂填充岩心实验
初始条件：
-
初始压力：10 MPa，降压至3.3 MPa
-
岩心：圆柱坐标系，31（径向）×52（轴向）网格
弹性模量模型对比（图3）：
-
常数E模型（E=200 MPa）：产水过快，位移模拟偏差大
-
饱和度相关E模型（E=200+700·S_MH MPa）：位移曲线与实验匹配良好
最佳匹配参数（基于Miyazaki等2005实验）：
-
E₀ = 200 MPa，ν = 0.217，β = 700 MPa
-
Poisson比范围：0.2-0.6，β范围：600-1000 MPa
温度分布验证：
-
降压导致水合物分解（吸热反应），温度从11°C降至2°C
-
变E模型的温度分布与实验匹配优于常数E模型
2. 降压法海床沉降预测（第2.9节，图6）
模型条件：
-
水合物储层厚度：15 m，初始水合物饱和度：60%
-
单口垂直井，初始压力13 MPa，降压至3 MPa（压降10 MPa）
弹性模量敏感性（E₀=100-400 MPa，β=700 MPa）：
-
软储层（E₀=100 MPa）：最大沉降约1 m（位于垂直井位置）
-
沉降随储层厚度增加而增加，随水合物饱和度降低而增加
环境风险：1 m量级的沉降可能诱发沉积层稳定性破坏和甲烷泄漏
3. 热水注入法矿场尺度模拟（第3节）
双水平井系统（图7-9）：
-
类似SAGD工艺：上井生产，下井注入85°C热水
-
两井垂直间距5 m，上井距储层顶部2 m
-
前期降压90天使两井连通，形成高渗区
储层与操作参数（表1）：
| 参数 | 值 |
|---|---|
| 储层面积 | 500 m × 50 m |
| 厚度 | 15 m |
| 初始压力/温度 | 13 MPa / 12.85°C |
| 孔隙度/渗透率 | 40% / 1000 mD |
| 初始水合物/水饱和度 | 60% / 40% |
| 热水注入温度/速率 | 85°C / 500 m³/天 |
| 生产井BHP | 3、6、9 MPa（三种方案） |
温度分布与流动方向（图10，BHP=3 MPa，1825天）：
-
流体向上井流动，水合物分解区边界不断扩展
-
热水 chamber 持续扩大，提供分解所需热量
4. 累计产气量对比（图11）
| 生产井BHP (MPa) | 降压法累计产气 | 热水注入法累计产气 | 比值（热水/降压） |
|---|---|---|---|
| 3 | 基准 | 1.6× | 1.6 |
| 6 | 基准 | 3.2× | 3.2 |
| 9 | 基准 | 12.3× | 12.3 |
关键发现：
-
热水注入法在BHP=9 MPa时的累计产气量 ≈ 降压法在BHP=3 MPa时的产气量
-
即：热水注入法可用更小的压降达到相同的产气量
5. 海床沉降对比（图12）
| 生产井BHP (MPa) | 降压法沉降 | 热水注入法沉降 | 说明 |
|---|---|---|---|
| 3 | 较大 | 略大于降压法（同BHP下） | 热水注入加速分解，弹性模量降低 |
| 6 | 中等 | 略大于降压法（同BHP下） | 同上 |
| 9 | 较小 | 略大于降压法（同BHP下） | 同上 |
关键发现：
-
沉降主要受储层压力（即BHP）控制，而非开采方法本身
-
相同BHP下，热水注入法沉降略大于降压法（因水合物分解更彻底，弹性模量降低）