📝 摘要
本研究旨在提出一种在数值模拟器中模拟页岩气流的工作流程，该流程考虑了页岩基质中存在的多重孔隙系统的复杂性及其内部不同的流动动力学。研究中采用了由有机质（干酪根）、无机质以及天然裂缝和水力裂缝组成的多重孔隙系统概念。干酪根因其高度多孔的性质、吸附气体的能力以及表面丰富的纳米孔隙而与其他页岩组分显著不同。
针对页岩微观尺度的物理过程已有一些理论提出。然而，在油藏尺度模型中，微观尺度模型中描述的孔隙系统细节往往会丢失。为解决这一问题，研究采用了动态表观渗透率的概念——它是基质压力的函数。这有助于将微观尺度模型的细节升尺度到油藏尺度模型，并有助于模拟达西流、菲克扩散以及基质与裂缝之间的过渡流。
研究通过Barnett页岩地层中一口水平井的实例验证了假设，该井的初始模型未考虑相关的流动现象。在整合了扩散和解吸后进行历史拟合发现，考虑这些额外过程会影响假定的储层改造体积，改变其体积和性质。这将是优化完井设计的关键因素，可以在相同或更高产量的情况下降低井成本。
生产预测结果表明，当未对页岩的相关假设进行建模时，基质对产量的贡献被低估。这强调了更好理解页岩中传输现象的重要性，这将使我们能够更深入地审视生产数据，并预测完井方法改变时产量的变化。这意味着多级高密度压裂作业可能无法优化井的价值。降低完井成本，使设计在初始阶段产生较低的产量和较低的递减率，将使我们能够在相同采收率下更长时间地生产这些井，从而将产量推向未来油气价格可能更好的时期。
🖥️ CMG软件应用情况总结
| 项目 | 内容 |
|---|---|
| 所用软件 | CMG IMEX（黑油模拟器）、CMOST（历史拟合优化） |
| 模型类型 | 单孔隙度模型（初始）→ 升级为考虑扩散和解吸的模型 |
| 研究区域 | Barnett页岩，Newark East油田 |
| 井型 | 水平井 + 水力压裂（裂缝网络通过微地震数据映射） |
| SRV尺寸 | 3,200 ft × 660 ft × 450 ft |
| 网格 | 111 × 145 × 1 |
| 裂缝网络 | 两组裂缝系统（水力裂缝 + 天然裂缝） |
| 模拟时长 | 8年历史拟合 + 30年生产预测 |
| 历史拟合工具 | CMG CMOST（DECE优化器） |
| 历史拟合误差 | 2.9%（最佳匹配作业#202） |
文中明确指出：CMG IMEX被用于页岩气藏的黑油模拟，CMOST用于辅助历史拟合和敏感性分析。
📊 主要模拟结果
1. 解吸效应（TOC敏感性）
| TOC (%) | 累计产气量增量（相对于无解吸） |
|---|---|
| 1 | +2.83% |
| 4 | +6.67% |
| 7 | +10.0% |
| 10 | +13.79% |
| 13 | +17.1% |
结论：TOC越高，解吸贡献越大。在Barnett页岩中，解吸对地质储量的贡献约为40%，但由于渗透率极低，实际产量有限。
2. 蒙特卡洛模拟结果（TOC服从正态分布 N(4.5, 1.16)）
| 输出参数 | P10 | P50 | P90 | 均值 |
|---|---|---|---|---|
| 累计产气量 (scf) | 1.165×10⁹ | 1.186×10⁹ | 1.206×10⁹ | 1.185×10⁹ |
| 平均储层压力 (psi) | 3,345 | 3,351 | 3,358 | 3,351 |
结论：TOC增加1%仅使产量增加约1.6%（1.8×10⁷ scf），经济影响有限。
3. 扩散效应
| 机制 | 累计产气量变化（相对于基准） | 压力变化 |
|---|---|---|
| 仅解吸（4% TOC） | +11.6% | 压力维持较好 |
| 仅扩散 | +25% | 加速压力消耗 |
| 扩散+解吸 | +37% | 介于两者之间 |
结论：扩散在低压下大幅增加气体传输能力，与解吸协同作用可显著提高产量。
4. 历史拟合前后关键参数变化
| 参数 | 历史拟合前 | 历史拟合后 | 变化 |
|---|---|---|---|
| 基质渗透率 (mD) | 7.65×10⁻⁵ | 1.0×10⁻⁴ | +30.7% |
| 天然裂缝渗透率 Set 1 (mD) | 1.01 | 0.09 | -91% |
| 天然裂缝渗透率 Set 2 (mD) | 56.6 | 51 | -10% |
| 基质孔隙度 (%) | 6.7 | 5.1 | -24% |
| TOC (wt%) | 未假设 | 4 | — |
| 最小井底压力 (psi) | 320 | 375 | +17% |
| 天然裂缝长度（东北方向）(ft) | 840 | 1,000 | +19% |
| 天然裂缝长度（西南方向）(ft) | 1,820 | 2,000 | +10% |
| 裂缝网络总长度 (ft) | 3,000 | 2,660 | -11% |
关键结论：考虑扩散和解吸后，裂缝网络总长度从3,000 ft减少到2,660 ft（减少约11%），表明在相同产量下，储层改造体积可以更小。
5. 30年生产预测（BHP约束）
| 模型 | 累计产气量 | 气体产量 | 平均储层压力 |
|---|---|---|---|
| 基准模型（无扩散/解吸） | 基准 | 基准 | 较高 |
| 更新模型（扩散+解吸） | +8.7% | +18% | -230 psi |
结论：扩散和解吸使基质贡献增加，后期产量更高，但压力消耗更快。
🔧 敏感性分析关键发现（井口压力）
| 参数 | 敏感性排名 | 说明 |
|---|---|---|
| 水力裂缝渗透率 | 1 | 最敏感，控制气体从基质到井筒的流动 |
| 基质孔隙度 | 2 | 控制OGIP |
| Langmuir压力 | 3 | 较高压力下释放更多吸附气 |
| TOC | 4 | 影响解吸量 |
| 天然裂缝长度 | 5 | 影响SRV范围 |
✅ 主要结论
-
解吸的贡献有限：在Barnett页岩中，吸附气约占地质储量的40%，但由于基质渗透率极低，实际解吸贡献有限。TOC从1%增加到13%时，产量仅增加17.1%。
-
TOC的经济影响有限：TOC增加1%仅使产量增加约1.6%，按2.82 USD/MMBtu气价计算，8年仅增加约50,760美元利润。
-
扩散是页岩气生产的关键机制：仅考虑扩散可使产量增加25%，远高于解吸的贡献。扩散和解吸协同作用可使产量增加37%。
-
扩散加速压力消耗：与解吸维持压力的作用相反，扩散加速储层压力消耗，使基质中的气体更快地排出。
-
SRV体积可优化：考虑扩散和解吸后，历史拟合得到的裂缝网络总长度减少了约11%（从3,000 ft降至2,660 ft），表明可以用更小的SRV实现相同的产量，从而降低完井成本。
-
历史拟合误差：最佳匹配模型的历史拟合误差为2.9%，验证了模型的可靠性。
-
生产预测：考虑扩散和解吸后，30年累计产气量增加8.7%，气体产量增加18%，但储层压力下降约230 psi。
🏛️ 作者及单位信息
| 项目 | 内容 |
|---|---|
| 作者 | Syed Munib Ullah Farid |
| 学位 | 理学硕士（石油工程） |
| 授予单位 | 德克萨斯A&M大学 |
| 论文提交日期 | 2015年5月 |
| 导师 | John E. Killough 教授（委员会主席） |
| 其他委员会成员 | Maria A. Barrusfet, Zoya Heidari |
| 资助 | 国际教育协会（IIE）富布赖特项目 |
📖 研究区域与储层信息
| 项目 | 参数 |
|---|---|
| 研究区域 | Barnett页岩，Newark East油田，德克萨斯州Fort Worth盆地 |
| 井深 | 7,000 ft |
| 储层温度 | 180°F |
| 初始压力 | 3,800 psi |
| 净厚度 | 415 ft |
| 气体比重 | 0.6 |
| 裂缝网络总长度（历史拟合后） | 2,660 ft |
| SRV尺寸 | 3,200 ft × 660 ft × 450 ft |
| TOC范围（文献） | 2.4-5.1%（均值约4.5%） |
💡 创新点
-
提出了从微观尺度到油藏尺度的升尺度工作流程：利用动态表观渗透率将微观尺度的多重孔隙流动机制（达西流、菲克扩散、过渡流）整合到油藏尺度模拟中
-
系统量化了解吸和扩散的相对贡献：发现扩散的贡献（+25%）远大于解吸（+11.6%），两者协同可使产量增加37%
-
揭示了扩散加速压力消耗的机制：与解吸维持压力的作用形成对比
-
利用蒙特卡洛模拟量化了TOC的不确定性：确定了TOC与产量的定量关系（TOC增加1%，产量增加1.6%）
-
证明了考虑扩散和解吸后SRV可以更小：裂缝网络总长度减少11%，对降低完井成本具有重要意义
-
利用微地震数据约束裂缝网络：结合CMOST进行辅助历史拟合，历史拟合误差低至2.9%
-
指出了高密度压裂的局限性：多级高密度压裂可能无法优化井的价值，适当的设计可以在初始产量较低的情况下实现更长的生产周期