📝 摘要
低至超低渗透率致密油藏已成为北美地区重要的油气供应来源。致密油藏的产量和压力瞬态分析是油藏研究人员面临的最困难问题之一，因为致密地层具有极端复杂性，例如多层生产，有效渗透率通常通过水力压裂得以增强。水平井完井结合多级水力压裂技术使得这类油藏的经济生产成为可能。为了从致密地层商业化生产油气，带有多级水力裂缝的自然完井（裸眼）或套管水平井是最常用的完井方式。多裂缝水平井的组合预期会产生复杂的流动序列。正确分析和识别流动形态的存在和出现顺序对于获取水力压裂优化和致密地层表征的有效信息至关重要。本文详细讨论了在尺寸可观的自然裂缝和均质（单孔隙度）地层中，水力压裂水平井的压力及其对应导数响应的诊断图，并展示了基于Saskatchewan西部Bakken和Viking目标地层的模拟结果与不同现有经验产量-时间关系之间的典型曲线匹配性能。使用商业油藏模拟器CMG IMEX生成了压力瞬态特性和生产数据特征。结果表明，完全填充和区域性的天然裂缝会显示不同的压力瞬态特征，从而显著影响井的生产性能。此外，这些高导流的相互连通的天然裂缝主导了致密地层中水平井的压力瞬态性能，即使存在水力裂缝也是如此。模拟还表明，如果储层存在一定程度天然裂缝，除非能实现大量水力裂缝和无限导流能力，否则水力压裂增产可能不会显著提高产能。通过与广泛接受的试井模型Kappa的对比，验证了模拟模型的可行性。还提供了现场案例研究以进行典型曲线拟合和EUR估算预测。
🖥️ CMG软件应用情况总结
| 项目 | 内容 |
|---|---|
| 所用软件 | • CMG IMEX（2010版本）：新一代自适应隐式-显式黑油模拟器，用于致密油藏压力瞬态性能的数值模拟 • CMG Builder：用于建立数值模拟模型 • Kappa Ecrin Saphir：压力瞬态分析（PTA）模块，用于与CMG模拟结果进行对比验证 |
| 模型类型 | • 三维黑油模型，笛卡尔网格系统 • 均质储层模型（单孔隙度） • 双重孔隙度模型（Warren & Root，1963）—天然裂缝储层 • 三重孔隙度模型（Al-Ahmadi，2010）—双裂缝系统 • 局部网格加密（LGR）技术用于水力裂缝建模 • 单相（油）流动，恒井底流压生产 |
| 模拟对象 | 致密油藏中水力压裂水平井的压力瞬态行为，基于Bakken和Viking地层的特性，以及Saskatchewan西部Dodsland Viking地层的现场生产数据 |
| 应用方式 | • 压力瞬态模拟：使用CMG IMEX模拟不同完井方式（裸眼/套管）、不同裂缝数量（1、2、4条）、不同裂缝间距、不同储层类型（均质/双重孔隙度/区域天然裂缝）下的无量纲压力及导数响应 • 模型验证：将CMG数值模拟结果与Kappa Ecrin Saphir解析解进行对比验证 • 双重孔隙度模型验证：与Warren & Root（1963）双重孔隙度模型解析解对比 • 三重孔隙度模型验证：与Al-Ahmadi（2010）三重孔隙度模型解析解对比 • 生产数据历史拟合：使用CMG生成的生产典型曲线，对Bakken和Viking地层实际生产数据进行历史拟合 • EUR预测：对比Arps、幂律指数（PLE）、 stretched exponential（SEPD）、Duong和Logistic Growth Model五种经验递减方法与CMG模拟典型曲线的EUR估算结果 |
| 压力瞬态模拟参数（表3.1-3.3） | • 水平井长度：330 m，井筒半径：0.09 m • 储层尺寸：Xe=900 m，Ye=375 m，厚度=18 m • 初始压力：5000 psi，井底流压：500 psi • 基质渗透率：0.01 md，基质孔隙度：0.10 • 天然裂缝渗透率：1 md，天然裂缝孔隙度：0.001 • 水力裂缝：半长60.96-121.91 m，宽度0.03 m，渗透率1000-10000 md • 流体粘度：2.50E-05 Pa·s，总压缩系数：1.74E-05 kPa⁻¹ |
| 主要模拟场景（第3章） | • 场景1：裸眼水平井（均质/全局天然裂缝/区域天然裂缝） • 场景2：纵向裂缝水平井（均质/全局天然裂缝/区域天然裂缝） • 场景3：多横向裂缝水平井（1、2、4条裂缝，不同间距） • 场景4：天然裂缝储层中带横向裂缝的水平井 • 敏感性分析：基质渗透率（0.01-100 md）、水力裂缝导流能力（100-10000 md） |
| 流动形态识别 | • 早期裂缝线性流（半斜率） • 双重孔隙度特征“凹槽” • 双线性流（1/4斜率） • 复合线性流 • 拟径向流 • 边界主导流（单位斜率） |
| 主要结论 | • 水力压裂水平井的压力瞬态行为受裂缝数量、间距、导流能力及天然裂缝存在显著影响 • CMG数值模拟与Kappa解析解在大多数情况下吻合良好，验证了CMG模型可靠性 • 双重孔隙度模拟成功再现了Warren & Root特征“凹槽” • 三重孔隙度模拟仅显示一个“凹槽”，原因是CMG无法严格模拟顺序流动 • 生产数据分类为4种典型递减类型，揭示了Bakken和Viking地层的不同流动形态组合 • 五种经验递减方法中，SEPD与CMG模拟典型曲线匹配最好；传统Arps法在b>1时可匹配长时线性流 |
文中明确指出（Case122.pdf）：
“The discussion is based on pressure-transient performances and characteristics of production data generated by employing a commercial reservoir simulator, CMG IMEX, a 3D finite-difference reservoir simulation package which is widely and popularly accepted by petroleum industry.” （摘要部分，第2页）
“Meng Wu helped me get on the road to using CMG and provided not only an experienced ear for hearing my doubts about writing a thesis, but a heart full of loving care and tenderness.” （致谢部分，第3页）
“Simulations are built by using a reservoir simulator (IMEX, 2010, Computer Modeling Group Ltd).” （第3.2节，第21页）
“For naturally fractured reservoirs, the dual-porosity of the Warren and Root option allows the discretization of matrix blocks …” （第4.1.1节，第59页）
“To demonstrate the applicability and practicability of the pressure transient behavior simulated by CMG, comparisons have been made with analytical solutions generated by Kappa.” （第1.1节，第3页）
“The simulation results presented in this paper are obtained by using the reservoir simulator CMG- IMEX.” （论文多处）
🧪 模拟方案与主要结果
1. 压力瞬态模拟场景（第3章）
场景1：裸眼水平井（图3.6）
| 储层类型 | 流动形态特征 |
|---|---|
| 均质（单孔隙度） | 早期裂缝线性流 → 边界主导流 |
| 全局天然裂缝（双重孔隙度） | 早期线性流 → 双重孔隙度“凹槽” → 中期线性流 → 边界主导流 |
| 区域天然裂缝 | 与全局天然裂缝类似，但中期出现单位斜率（外区单孔隙度导致） |
场景2：纵向裂缝水平井（图3.7）
-
纵向裂缝与裸眼水平井流动形态相似
-
双重孔隙度“凹槽”仍然主导
场景3：多横向裂缝水平井（图3.8-3.11）
| 裂缝数量 | 裂缝间距 | 关键观察 |
|---|---|---|
| 1条 | — | 早期线性流 → 过渡期 → 边界主导流 |
| 2条（小间距，60 m） | 与1条裂缝曲线重合 | |
| 2条（大间距，270 m） | 与4条裂缝曲线在第二线性流重合 | |
| 4条 | 最小压降，线性流持续时间最长 |
套管井 vs 裸眼井（图3.10-3.11）：
-
裸眼井压降小于套管井（更多流动通道）
-
裸眼井显示复合线性流；套管井仅显示双线性流（有限导流裂缝+均质基质）
场景4：天然裂缝储层中带横向裂缝的水平井（图3.12-3.15）
-
全局天然裂缝：线性流 → 双重孔隙度“凹槽” → 中期线性流 → 边界主导流
-
区域天然裂缝：线性流 → “凹槽” → 接近单位斜率线 → 边界主导流
-
内区高渗天然裂缝+外区低渗单孔隙度导致压力传播延迟
纵向 vs 横向裂缝对比（图3.16-3.17）：
-
相同裂缝面积下，纵向裂缝压降更小
-
原因：纵向裂缝与水平井筒接触面积更大，避免“流动扼流”效应
2. 敏感性分析（第3.5节）
基质渗透率敏感性（图3.18）：
| 基质渗透率 | 影响 |
|---|---|
| 0.01 md（基案例） | 早期数据不规则 |
| 1 md | 早期径向流出现 |
| 10 md | 早期径向流更明显，边界主导流更早 |
| 100 md | 拟径向流清晰，边界主导流显著 |
水力裂缝导流能力敏感性（图3.19-3.20）：
| 裂缝渗透率 | 1条裂缝 | 4条裂缝 |
|---|---|---|
| 100 md | 双线性流（1/4斜率）+线性流 | 双线性流+两段线性流 |
| 1000 md | 双线性流减弱 | 双线性流+两段线性流 |
| 10000 md | 双线性流消失，仅裂缝线性流 | 双线性流+两段线性流 |
3. CMG与Kappa解析解验证（第4.2节）
验证案例（图4.1-4.7）：
| 案例 | 验证结果 |
|---|---|
| 均质储层水平井（图4.1） | 整体匹配良好，早期CMG未显示预期径向流 |
| 均质储层垂直井（图4.2） | 良好匹配 |
| 均质储层水平井+2横向裂缝（图4.3） | 良好匹配 |
| 均质储层水平井+4横向裂缝（图4.4） | 良好匹配 |
| 天然裂缝储层水平井（图4.5） | 良好匹配，双重孔隙度“凹槽”清晰 |
| 天然裂缝储层+2裂缝（图4.6） | 良好匹配，早期有差异 |
| 天然裂缝储层+4裂缝（图4.7） | 良好匹配 |
结论：CMG数值模拟与Kappa解析解在大多数情况下吻合良好，验证了CMG模型的可靠性。早期差异归因于CMG有限差分法的数值扰动。
4. 双重/三重孔隙度模型验证（第4.3节）
双重孔隙度模型（Warren & Root, 1963）（图4.9）：
-
CMG成功再现了双重孔隙度特征“凹槽”
-
流动形态：早期径向流 → “凹槽”（基质向裂缝补充）→ 拟径向流 → 边界主导流
-
CMG与解析解匹配良好
三重孔隙度模型（Al-Ahmadi, 2010）（图4.12）：
-
模型组成：基质 + 微裂缝 + 宏裂缝（顺序流动）
-
解析解预期：两个“凹槽”（基质→微裂缝、微裂缝→宏裂缝）
-
CMG模拟：仅显示一个“凹槽”
-
原因：CMG无法严格模拟顺序流动假设（同时排水掩盖了第二个“凹槽”）
5. Bakken地层生产数据分类（第5.3节，图5.2-5.5）
| 类型 | 特征 | 解释 |
|---|---|---|
| 类型1（图5.2） | 两段线性流（阶梯状） | 内区线性流 + 外区复合线性流 |
| 类型2（图5.3） | 长时线性流（近2个对数周期） | 大型高导流裂缝系统或天然裂缝 |
| 类型3（图5.4） | 双线性流 + 两段线性流 | 复杂裂缝-基质系统 |
| 类型4（图5.5） | 双线性流 + 线性流 + 边界主导流 | 有限裂缝半长，边界效应明显 |
6. 经验递减方法与CMG典型曲线对比（第5.4-5.5节）
五种经验递减方法：
-
Arps（双曲/指数/调和）
-
幂律指数（PLE）
-
Stretched Exponential（SEPD）
-
Duong方法
-
Logistic Growth Model（LGM）
CMG典型曲线参数（图5.11-5.14）：
-
无量纲产量：qD = qμB/[Kh(Pi-Pwf)]
-
无量纲时间：tD = kt/(φμCtA)
Bakken类型1匹配（图5.11）：
-
实际数据落在4条横向裂缝曲线上
-
阶梯状过渡难以完美匹配
Bakken类型2匹配（图5.12）：
-
长时线性流（8年以上）
-
裸眼水平井，13条裂缝，半长250 ft
Bakken类型3匹配（图5.13）：
-
双线性流+两段线性流
-
实际数据与4条横向裂缝匹配（单孔隙度）
Bakken类型4匹配（图5.14）：
-
双线性流+线性流+边界主导流
-
实际数据与4条横向裂缝匹配（双重孔隙度）
Dodsland Viking地层案例（图5.16-5.18，表5.1）：
-
7口水平井，生产数据（2008-2014）
-
流动形态：线性流（半斜率）→ 边界主导流（单位斜率）
-
无径向流
-
CMG匹配：2条横向裂缝，裂缝半长200 ft，天然裂缝储层（双重孔隙度）
7. EUR预测结果（第5.6.3节，图5.19）
| 方法 | EUR预测表现 |
|---|---|
| SEPD | 最保守，最可靠 |
| Arps（b>1） | 可匹配长时线性流，略高估 |
| PLE | 与Arps相似 |
| Duong | 早期下降过快，可能低估 |
| LGM | 最差，不适用于单井 |
✅ 主要结论
-
CMG IMEX的压力瞬态模拟能力：CMG IMEX能够有效模拟致密油藏中带有多条水力裂缝的水平井的压力瞬态行为，包括均质、双重孔隙度和区域天然裂缝储层。模拟结果与Kappa Ecrin Saphir解析解吻合良好，验证了模型可靠性。
-
水力裂缝参数的重要性：裂缝数量、间距和导流能力显著影响压力瞬态响应。4条横向裂缝比2条或1条裂缝产生更小的压降和更长的线性流持续时间。裂缝间距影响中间时间的复合线性流行为。
-
天然裂缝的主导作用：在存在天然裂缝的储层中，即使有人工水力裂缝，天然裂缝系统仍主导压力瞬态行为。双重孔隙度特征“凹槽”在所有天然裂缝案例中清晰可见。
-
完井方式的影响：裸眼水平井（允许井筒流动）比套管井（仅裂缝段射孔）产生更小的压降，并显示复合线性流，而套管井仅显示双线性流。
-
三重孔隙度模拟的局限性：CMG无法严格模拟顺序流动假设（基质→微裂缝→宏裂缝），因此在三重孔隙度模拟中仅显示一个“凹槽”，而非解析解预期的两个“凹槽”。
-
生产数据分类与典型曲线匹配：Bakken地层150口井可分类为4种典型递减类型。CMG生成的生产典型曲线成功匹配了不同类型的历史生产数据，为EUR预测提供了更可靠的依据。
-
经验递减方法对比：在五种经验递减方法中，SEPD与CMG典型曲线匹配最好且预测最保守。传统Arps法在b>1时可匹配长时线性流。LGM不适用于单井EUR预测。
🏛️ 作者及单位信息
-
作者: Yue Zhu（朱悦）
-
单位: 里贾纳大学（University of Regina），石油系统工程系，萨斯喀彻温省，加拿大
-
学位: 应用科学硕士（Master of Applied Science）
-
毕业时间: 2015年7月
-
联合导师: Dr. Gang Zhao（赵刚，共同导师），Dr. Fanhua Zeng（曾凡华，委员会成员）
-
委员会成员: Dr. Amy Veawab，Dr. Hairuo Qing（外部考官）
-
资助: 里贾纳大学研究生与研究学院，石油技术研究中心（PTRC）奖学金
💡 补充说明
该硕士论文是致密油藏压力瞬态分析和产量递减分析领域的系统性数值模拟研究，核心创新点在于：
-
全面的模拟场景覆盖：系统研究了4种主要场景（裸眼水平井、纵向裂缝水平井、多横向裂缝水平井、天然裂缝储层中带裂缝水平井），涵盖了致密油藏水平井压裂的主要完井方式。
-
多软件交叉验证：将CMG IMEX数值模拟结果与Kappa Ecrin Saphir解析解进行详细对比验证（11个典型案例），并进一步与Warren & Root（1963）双重孔隙度模型和Al-Ahmadi（2010）三重孔隙度模型解析解对比，建立了CMG模拟的可靠性基准。
-
三重孔隙度模拟的深入分析：首次系统分析了CMG在三重孔隙度模拟中的局限性（仅显示一个“凹槽” vs 解析解预期的两个“凹槽”），归因于CMG无法严格模拟顺序流动假设，为类似模拟工作提供了重要警示。
-
生产数据分类体系：基于150口Bakken地层井的生产数据，建立了4种典型递减类型分类体系（类型1-4），为致密油藏生产分析提供了实用参考框架。
-
经验递减方法与CMG典型曲线的系统对比：对比了5种经验递减方法（Arps、PLE、SEPD、Duong、LGM）与CMG生成的生产典型曲线在Bakken和Viking地层的匹配性能，确定SEPD为最可靠方法。
-
Dodsland Viking地层的详细案例研究：对7口水平井进行了完整的历史拟合分析，识别了线性流（半斜率）→边界主导流（单位斜率）的流动形态，并成功用双重孔隙度模型（2条横向裂缝，半长200 ft）匹配。
该研究对从事致密油/页岩油藏压力瞬态分析、产量递减分析、水力压裂水平井数值模拟、以及非常规油藏开发方案设计的工程师和科研人员具有重要的参考价值。