📝 摘要
煤层气（CBM）是一种非常规天然气资源，甲烷主要以吸附状态储存在煤基质的微孔隙中。生产时需要先排水降压，使甲烷解吸并流入裂隙系统，最终进入井筒。与常规气藏不同，相邻井之间的干扰对煤层气生产是有益的，可以加速排水降压过程，从而更早达到更高的产气峰值。同时，煤层中存在的面割理和端割理导致了渗透率的各向异性，这影响了井网的优化布置。此外，Palmer-Mansoori 理论指出，随着压力降低，基质收缩效应会提高有效渗透率，从而增加最终采收率。
本研究使用 CMG 模拟器（GEM）对这些效应进行了系统的敏感性分析。通过改变井距、井网布置方式（考虑渗透率各向异性）以及 Palmer-Mansoori 模型中的关键参数（杨氏模量、泊松比、极限应变），定量评估了各因素对产气速率和累积产气量的影响。研究结果表明：更小的井距可以加速排水并获得更高的早期产气速率；在渗透率各向异性储层中，井网应沿低渗透率方向布置以获得最佳产量；杨氏模量和极限应变的增加以及泊松比的减小均会提高产气量。本研究为煤层气藏的高效开发提供了定量的工程指导。
🖥️ CMG软件应用情况总结
| 项目 | 内容 |
|---|---|
| 所用软件 | • CMG GEM™：组分模拟器，用于煤层气藏的双重介质（基质+裂隙）模拟，包括吸附/解吸、扩散（Fick定律）、渗流（Darcy定律）以及Palmer-Mansoori 渗透率变化模型的模拟（论文第37页、第80-84页）。 |
| 模型类型 | 三维双重介质模型（DUALPOR），21×21×1 网格，单层，均质。 |
| 模拟对象 | 煤层气藏（80 英亩，厚度 30 ft），包含基质系统（储气）和裂隙系统（面割理+端割理，导流）。 |
| 应用方式 | • 基础模型建立：使用 CMG-GEM 建立双重介质（Warren-Root 模型）煤层气藏模型，设置基质和裂隙的孔隙度、渗透率、裂隙间距等参数。 • 吸附/扩散模拟：使用 LANG-DIFFUSION-COAL 模型，通过 Langmuir 等温吸附线描述吸附气含量与压力的关系，使用 Fick 定律模拟气体从基质向裂隙的扩散过程（COAL-DIF-TIME 参数）。 • 井距影响模拟：在总面积为 80 英亩的储层中，分别模拟了 80 英亩、40 英亩、20 英亩和 5 英亩井距（对应 1、2、4、16 口井）的生产动态。 • 渗透率各向异性与井网优化模拟：设置裂隙渗透率 kx=1 md, ky=0.01 md（各向异性比 100:1），对比三种井网布置方式（正方形、沿低渗方向拉长、沿高渗方向拉长）对产量的影响。 • 各向异性向等向性转换验证：利用转换公式将各向异性方形储层转换为等向性矩形储层，对比两者的模拟结果以验证转换方法的有效性。 • Palmer-Mansoori 模型敏感性分析：激活 Palmer-Mansoori 模型，分别改变杨氏模量（E）、泊松比（ν）和极限应变（ε_L），分析各参数对产气速率和累积产气量的影响。 |
| 油藏模型参数（基础案例） | • 网格：21×21×1 • 储层尺寸：1866.76 ft × 1866.76 ft × 30 ft（80 英亩） • 初始压力：725.18 psia • 储层温度：113 °F • 裂隙孔隙度：0.001（0.1%） • 裂隙渗透率：1 md（各向同性基础案例） • 基质孔隙度：0.005（0.5%） • 基质渗透率：0.00001 md • Langmuir 压力：725.18 psia • Langmuir 体积：32.3 gmol/lb of rock • 扩散时间：100 天 • 井半径：0.1197 ft |
| 模拟方案 | • 井距方案：80 英亩、40 英亩、20 英亩、5 英亩（对应 1、2、4、16 口井） • 井网方案（各向异性，kx:ky = 1:0.01）：A（正方形）、B（沿低渗方向拉长）、C（沿高渗方向拉长） • 各向异性向等向性转换验证：kx:ky = 1:0.1 和 1:0.01 • Palmer-Mansoori 参数敏感性： – 杨氏模量 E：500k ~ 5,000k psi – 泊松比 ν：0.1 ~ 0.5 – 极限应变 ε_L：0.0025 ~ 0.02 |
| 主要结论 | • 井距影响：井距越小（井数越多），排水越快，产气峰值出现更早且更高。10 年模拟中，5 英亩井距的采出程度（82.23%）显著高于 80 英亩（36.27%）。 • 各向异性与井网：在渗透率各向异性（kx > ky）储层中，井网应沿低渗透率方向（y 方向）拉长布置（Type B），可获得最高的产气量和采出程度。 • 各向异性转换方法验证：CMG 模拟验证了将各向异性方形储层转换为等向性矩形储层的方法是有效的，两者产气曲线重合。 • Palmer-Mansoori 模型参数影响： – 杨氏模量 E 越大（岩石越“硬”），产气量越高。 – 泊松比 ν 越小，产气量越高。 – 极限应变 ε_L 越大，产气量越高。 • CMG 适用性：CMG-GEM 能够有效模拟煤层气的双重介质特性（Warren-Root 模型）、吸附/解吸（Langmuir 模型）、扩散（Fick 定律）以及渗透率随压力变化（Palmer-Mansoori 模型）等复杂过程，是煤层气藏开发方案优化的有力工具。 |
文中明确指出（Case136.pdf）：
“These effects and phenomena can be modeled with the CMG simulator.” （第3页，摘要部分）
“The simulation studies are conducted using CMG- GEM 2004.10.” （第37页，第5.1节）
“CMG (Computer Modeling Group) simulator can be used to model the coal bed methane reservoir.” （第80页，附录A）
“CMG- GEM 2004.10 is a two- phase compositional simulator that has the ability of modeling coalbed methane reservoir including Palmer and Mansoori Theory.” （第80页，附录A）
“The Palmer and Mansoori model is a built in feature of coal bed methane simulation on CMG-GEM.” （第94页，附录D）
🧪 模拟方案与主要结果
1. 基础模型与井距影响研究（第5.3节、第6.2.1节，图6.3-6.8，表6.1）
模型参数：双重介质，21×21×1 网格，80 英亩，裂隙 k=1 md（各向同性），初始压力 = Langmuir 压力 = 725.18 psia（饱和储层）。
| 井距（英亩） | 井数 | 10年累计产气量（scf） | 气体采出程度 | 10年累计产水量（stb） | 水采出程度 | 关键结论 |
|---|---|---|---|---|---|---|
| 80 | 1 | 3.28E+08 | 36.27% | 17077.5 | 95.24% | 峰值最低、最晚 |
| 40 | 2 | 4.71E+08 | 52.09% | 17783.3 | 99.17% | 产气量提升 |
| 20 | 4 | 5.89E+08 | 65.15% | 17914.2 | 99.90% | 产气量提升 |
| 5 | 16 | 7.43E+08 | 82.23% | 17931.7 | 100.00% | 峰值最高、最早，累积产量最大 |
关键发现（图6.3-6.8）：
-
产气速率：井距越小，早期产气速率越高，峰值出现越早。5 英亩井距的峰值产气速率显著高于其他方案。约 2-4 年后，小井距方案的产气速率开始低于大井距方案，但累积产量仍更高（图6.5）。
-
产水速率：井距越小，早期排水速率越高，5 英亩井距在约 1 年内即可排出 95% 以上的水（图6.6、图6.8）。
-
采出程度：10 年模拟期内，井距越小，气体采出程度越高（表6.1，图6.7）。
-
机理：小井距加速了压力干扰和扩散，提前启动了解吸过程。
2. 渗透率各向异性与井网优化（第5.4节、第6.2.3节，图6.11-6.15，表6.2、6.3）
模型参数：80 英亩储层，裂隙各向异性 kx=1 md, ky=0.01 md（方向性渗透率比为 100:1），基质渗透率各向同性（0.00001 md）。
| 井网类型 | 描述 | 20 英亩单井（4口井）10年采出程度 | 5 英亩单井（16口井）10年采出程度 | 关键结论 |
|---|---|---|---|---|
| Type A | 正方形 | 15.92% | 42.48% | 常规布井方式，产量中等 |
| Type B | 沿低渗方向（y）拉长 | 21.58% | 52.29% | 最优布井方式，累积产气量最高 |
| Type C | 沿高渗方向（x）拉长 | 10.49% | 24.99% | 最差布井方式，产量最低 |
关键发现：
-
最优井网：Type B（沿低渗透率方向拉长布置井网）在 20 英亩和 5 英亩井距下均获得最高的气体采出程度（图6.11、图6.12、表6.2、6.3）。
-
机理：各向异性储层在等效等向性转换后变为一个“瘦长”矩形，长边沿低渗方向。Type B 的井网布置与此等效矩形的形状匹配，可以更有效地覆盖储层，加速排水和降压（图6.13-6.15）。
3. 各向异性向等向性转换验证（第5.5节、第6.2.2节，图6.9-6.10，附录E）
转换方法：将各向异性方形储层（kx, ky）通过公式转换为等效等向性矩形储层（k = √(kx·ky)，新 x、y 尺寸按比例缩放）。
| 各向异性比 | 原始方形储层（kx, ky） | 等效等向性矩形储层（k, 新 x, 新 y） | 模拟验证结果 |
|---|---|---|---|
| 10:1 (1 md, 0.1 md) | 1866.76 ft × 1866.76 ft | k=0.316 md, 1049.8 ft × 3319.6 ft | 产气速率曲线重合（图E.15、E.16） |
| 100:1 (1 md, 0.01 md) | 1866.76 ft × 1866.76 ft | k=0.1 md, 590.3 ft × 5903.2 ft | 产气速率曲线重合（图6.9、6.10） |
关键结论：CMG 模拟结果验证了该转换方法的正确性，为理解各向异性对井网的影响提供了理论依据。等效等向性矩形的长边沿原始低渗方向，解释了为何 Type B 井网最优。
4. Palmer-Mansoori 模型参数敏感性分析（第5.6节、第6.2.4节，图6.16-6.21）
模型机理：Palmer-Mansoori 模型模拟了裂隙在压力降低时的两种竞争效应：裂隙压缩（降低渗透率）和基质收缩（增大裂隙开度，提高渗透率）。在低压下，基质收缩可能占主导，产生“渗透率反弹”。
| 参数 | 参数变化 | 对产气速率的影响 | 对累积产气量的影响 | 关键结论 |
|---|---|---|---|---|
| 杨氏模量（E） | 增加（500k → 5,000k psi） | 增加（图6.16） | 增加（图6.17） | E 越大，岩石越“硬”，裂隙压缩效应相对减弱，渗透率更高。 |
| 泊松比（ν） | 减小（0.5 → 0.1） | 增加（图6.18） | 增加（图6.19） | ν 越小，岩石更“脆”，基质收缩效应相对更有效。 |
| 极限应变（ε_L） | 增加（0.0025 → 0.02） | 增加（图6.20） | 增加（图6.21） | ε_L 越大，基质收缩潜力越大，渗透率反弹更显著。 |
附加发现（附录D）：
-
当 Palmer-Mansoori 模型激活时：CMG-GEM 会忽略直接输入的裂隙/基质压缩系数（CCPOR, CPRPOR），改用模型内嵌公式计算。
-
当模型未激活时：改变裂隙压缩系数对产量有微小影响，改变基质压缩系数无影响。
✅ 主要结论
-
井距影响（干扰效应）：煤层气藏中，井间干扰是有益的。井距越小，排水降压越快，产气峰值出现更早且更高。10 年模拟期内，5 英亩井距的气体采出程度（82.23%）远高于 80 英亩（36.27%）。
-
渗透率各向异性与井网优化：煤层中面割理和端割理导致渗透率各向异性（通常 kx >> ky）。通过转换方法验证，等效等向性储层为“瘦长”矩形（长边沿低渗方向）。因此，最优布井方式是沿低渗透率方向拉长井网（Type B），可获得最高产量。
-
各向异性转换方法验证：CMG 模拟成功验证了将各向异性方形储层转换为等向性矩形储层的数学方法，两者产气曲线完全重合。
-
Palmer-Mansoori 模型参数影响：
-
杨氏模量（E）越大，产量越高。
-
泊松比（ν）越小，产量越高。
-
极限应变（ε_L）越大，产量越高。
-
激活该模型后，CMG 会自动覆盖直接输入的压缩系数（CCPOR/CPRPOR）。
-
-
CMG 软件的适用性：CMG-GEM 能够有效模拟煤层气的双重介质特性（Warren-Root 模型）、吸附/解吸（Langmuir 模型）、扩散（Fick 定律）以及渗透率随压力变化（Palmer-Mansoori 模型）等复杂过程，是煤层气藏开发方案（井距、井网、生产策略）优化的有力工具。
🏛️ 论文主要信息
| 项目 | 内容 |
|---|---|
| 论文标题 | Simulation Study of the Effect of Well Spacing, Effect of Permeability Anisotropy, and Effect of Palmer and Mansoori Model on Coalbed Methane Production |
| 作者 | Ismail Zulkarnain |
| 学位 | 理学硕士（Master of Science） |
| 授予单位 | 德克萨斯农工大学（Texas A&M University） |
| 毕业年份 | 2005 年 12 月 |
| 所属院系 | 石油工程系（Harold Vance Department of Petroleum Engineering） |
| 研究领域 | 煤层气、非常规油气藏、油藏数值模拟、双重介质、渗透率各向异性 |
💡 补充说明
该硕士论文是一项关于煤层气（CBM）开发中关键工程参数（井距、井网、基质收缩效应）的系统性数值模拟研究，核心创新点在于：
-
多因素系统对比：首次在同一基础模型中系统地比较了井距（4 种）、井网布置（3 种）和 Palmer-Mansoori 模型参数（3 类）对 CBM 生产的影响，提供了定量的工程指导。
-
井距影响量化：定量证明了小井距对 CBM 的“有益干扰”效应，并给出了 10 年模拟期内不同井距方案的具体采出程度数据（从 36.27% 到 82.23%），为经济性分析提供了基础。
-
各向异性与井网优化验证：结合数学转换方法，通过 CMG 模拟验证了“各向异性方形储层等效于等向性矩形储层”的理论，并据此明确提出了“沿低渗方向拉长井网”的最优布井策略。
-
Palmer-Mansoori 模型参数敏感性分析：首次系统地定量分析了杨氏模量、泊松比和极限应变对 CBM 产量的影响规律，并澄清了激活该模型后 CMG 软件对压缩系数输入的处理方式。
-
CMG 软件的双重介质与吸附模型应用：展示了 CMG-GEM 在模拟 CBM 独特生产机理（排水→解吸→扩散→渗流）方面的完整工作流，包括双重介质设置、Langmuir 等温吸附、Fick 扩散以及压力敏感渗透率模型。
-
实践指导性强：研究结论对 CBM 气田开发中的井距选择、井网设计以及考虑基质收缩效应的生产预测具有直接的指导意义。
该研究对从事煤层气、页岩气等非常规天然气开发、油藏工程、数值模拟以及开发方案优化的工程师和科研人员具有重要的参考价值。