📝 摘要
对俄亥俄州 Youngstown 地震事件时空分布的分析表明，2011 年 1 月至 12 月期间向 North Star #1 盐水处置井的注入触发了地震活动。研究发现，诱发地震的震源深度远低于注入层位，这引发了一个关键问题：孔隙压力如何能够迁移到注入层以下超过 1 公里的前寒武纪基底岩石中，从而诱发地震？本研究通过数值模拟，分析了注入过程中地质和地质力学特征的影响，评估了水力压裂的潜力，以及导水裂缝带或断层的存在性。结果表明，导水裂缝和断层可作为流体通道，将压力传递到更深部、临界应力的基底断层，是深部诱发地震的主要原因之一。此外，断层带内的非均质性可能是地震震源随时间侧向迁移的原因之一。研究使用 CMG-GEM 模拟器构建了三维流体流动模型，分析了不同断层距离、渗透率和裂缝网络对深部压力传递的影响。
🖥️ CMG软件应用情况总结
| 项目 | 内容 |
|---|---|
| 所用软件 | • CMG GEM：组分模拟器，用于构建三维流体流动模型，模拟盐水注入期间的压力扩散和孔隙压力变化（论文第2页，第2节）。 • CMG 迭代耦合方法：用于耦合流体流动与地质力学，分析应力状态变化（论文第4页）。 • CMG 双重渗透率模型：用于模拟裂缝网络中的压力传递（论文第4页）。 |
| 模型类型 | 三维流体流动模型（30×30×25 网格，12 km × 12 km × 2.5 km），用于模拟盐水注入期间的压力响应。 |
| 模拟对象 | 俄亥俄州 Youngstown 地区 North Star #1 盐水处置井，注入层为 Basal Cambrian Sandstone，震源位于前寒武纪结晶基底（深度 ~3.5-4 km）。 |
| 应用方式 | • 流体流动建模：使用 CMG-GEM 建立三维等温流体流动模型，忽略水动力弥散和化学反应，模拟注入期间的压力扩散。 • 历史拟合：通过调整注入储层的渗透率（3 mD），对地面压力数据进行历史拟合（图2a）。 • 无断层/裂缝基底压力分析：模拟无断层或裂缝网络时基底的压力响应（图2b）。 • 断层参数化分析：在基底中设置不同距离（30 m、100 m、800 m）和不同渗透率（10 mD、50 mD、100 mD）的断层，分析其对深部孔隙压力增加的影响（图4、图5）。 • 地质力学耦合：使用 CMG-GEM 的迭代耦合方法，分析流体注入对最小和最大水平应力状态的影响（图6）。 • 裂缝网络双重介质模拟：使用 CMG-GEM 的双重渗透率模型（dual-permeability）模拟裂缝网络（裂缝渗透率 5 mD，间距 150 m）中的压力传递（图7a）。 • 非均质性建模：模拟沿裂缝网络的渗透率非均质性（从注入井附近 10 mD 逐渐降低到远端 1 mD），分析地震震源侧向迁移的机制（图7b、7c）。 |
| 模型参数 | • 网格：30×30×25 • 模型尺寸：12 km × 12 km × 2.5 km • 注入层渗透率：3 mD • 基底渗透率：0.1 mD（假设值） • 注入层孔隙度：0.1 • 流体压力梯度：10 MPa/km • 裂缝梯度：14.7 MPa/km • 裂缝网络渗透率：5 mD • 裂缝网络间距：150 m • 断层渗透率：10-100 mD • 断层厚度：5 m • 断层长度：1 km • 断层距离：30 m、100 m、800 m • 注入速率：86-316 m³/天（4个阶段） |
| 模拟方案 | • 方案一（无断层）：基底无断层或裂缝网络，分析压力能否扩散到震源深度 • 方案二（水力压裂评估）：使用有限元水力压裂模拟器评估注入是否会导致水力裂缝扩展 • 方案三（断层参数化）：改变断层距离（30-800 m）和渗透率（10-100 mD） • 方案四（裂缝网络）：使用双重渗透率模型模拟裂缝网络中的压力传递 • 方案五（非均质性裂缝网络）：模拟渗透率沿裂缝网络侧向递减（10 → 3 → 1 mD） |
| 主要结论 | • 无断层时无法解释深部地震：当基底无断层或裂缝网络时，压力无法扩散到 3.5-4 km 的震源深度（图2b）。 • 水力压裂不可能：注入速率不足以引起水力裂缝扩展，净压力主要为负值（图3c、3d）。 • 导水断层是关键：只有在靠近注入井（30-100 m）且具有高渗透率（50-100 mD）的断层存在时，孔隙压力才能显著增加并到达深部基底（图4b、图5）。 • 临界应力状态：即使存在高渗透断层，孔隙压力增量也不大，因此断层必须处于临界应力状态，微小的压力变化即可触发地震。 • 裂缝网络可传递压力：双重渗透率模型显示，裂缝网络可将压力传递至 4 km 深度（图7a）。 • 非均质性解释侧向迁移：沿裂缝网络的渗透率非均质性可导致压力前锋随时间逐步向外扩展，从而解释地震震源从注入井附近向外侧向迁移的现象（图7b、7c）。 • CMG GEM 的适用性：CMG GEM 能够有效模拟流体注入过程中的压力扩散、断层/裂缝网络中的优先流动路径、以及流体-地质力学耦合效应，是注入诱发地震风险评估的有力工具。 |
文中明确指出（Case144.pdf）：
“A 3D fluid-flow simulator, CMG-GEM, is used to model the behavior of brine plumes during injection periods.” （第2页，第2节）
“A fault zone is modeled within the previously described CMG-GEM model by placing grids of small width and high permeability in the Precambrian basement…” （第3页）
“Simulation of fractured reservoirs using the dual-permeability approach involves discretization of the solution domain into two continua, called the matrix and the fracture… A fracture network is modeled using a dual-permeability model in CMG-GEM…” （第4页）
“An iterative coupling approach in CMG-GEM, used for coupling fluid flow and geomechanics, shows that the total minimum and maximum horizontal stresses also can change in the basement formation…” （第4页）
🧪 模拟方案与主要结果
1. 模型建立与无断层情况（第2节，图2）
模型参数：30×30×25 网格，12 km × 12 km × 2.5 km，注入层渗透率 3 mD（历史拟合），基底渗透率 0.1 mD。
| 注入阶段 | 时间（天） | 注入速率（m³/天） | 历史拟合结果 |
|---|---|---|---|
| 1 | 0-70 | 86 | 模拟井口压力与实测值匹配良好（图2a） |
| 2 | 70-100 | 187 | |
| 3 | 100-200 | 259 | |
| 4 | 200-327 | 316 |
无断层结果（图2b）：
-
压力增加主要集中在注入层位
-
在 3.5-4 km 震源深度处无压力变化
-
即使使用较高的基底渗透率（0.1 mD），仍无法解释深部地震
2. 水力压裂潜力评估（第3节，图3）
方法：耦合流体流动、地质力学和水力压裂模拟，使用偶极声波测井数据确定动态弹性参数。
| 参数 | 值 | 关键结论 |
|---|---|---|
| 最小水平应力梯度 | 14.7 MPa/km | 注入速率不足以引起水力裂缝扩展 |
| 最大水平应力梯度 | 25 MPa/km | 裂缝内部净压力主要为负值（图3d） |
| 地层动态杨氏模量 | ~30-50 GPa | 孔隙弹性效应增加了裂缝扩展阻力 |
| 地层动态泊松比 | ~0.2-0.3 |
3. 断层参数化分析（第4节，图4-6）
模型：在基底中设置高渗透断层带（宽度 5 m，长度 1 km），不同距离和渗透率。
| 断层距离 | 断层渗透率 | 4 km 深度压力增量 | 关键结论 |
|---|---|---|---|
| 800 m | 100 mD | 无显著增加 | 断层必须靠近注入井 |
| 100 m | 10 mD | 较小 | 低渗透断层压力传递有限 |
| 100 m | 50 mD | 中等 | |
| 100 m | 100 mD | 显著增加 | 导水断层是关键 |
| 30 m | 100 mD | 最大 | 越近效果越显著 |
地质力学耦合（图6）：
-
孔隙压力增加导致有效应力状态变化，但变化幅度不显著
-
表明断层必须处于临界应力状态，微小压力扰动即可触发地震
-
低摩擦系数（0.6-0.8）断层更易被激活
4. 裂缝网络双重介质模拟（第4节，图7）
模型：使用 CMG-GEM 双重渗透率模型，裂缝渗透率 5 mD，间距 150 m。
| 裂缝网络类型 | 4 km 深度压力响应 | 侧向压力分布 | 关键结论 |
|---|---|---|---|
| 均质（5 mD） | 压力可传递至 4 km（图7a） | 各点压力同步增加 | 会导致地震同时发生（“群震”），无法解释侧向迁移 |
| 非均质（10→3→1 mD） | 压力可传递至 4 km | 压力前锋随时间逐步外扩（图7c） | 可解释地震震源侧向迁移 |
非均质模型结果（图7b、7c）：
-
点 A（80 天，近井）：压力达到触发值
-
点 B（200 天，300 m）：压力达到触发值
-
点 C（300 天，650 m）：压力达到触发值
-
时间与 Youngstown 实际地震序列一致
✅ 主要结论
-
无断层/裂缝时无法解释深部地震：当基底无断层或裂缝网络时，压力无法扩散到 3.5-4 km 的震源深度。即使使用较高的基底渗透率（0.1 mD），压力仍被限制在注入层附近。
-
水力压裂不可能：North Star #1 井的注入速率不足以引起水力裂缝扩展。裂缝内部净压力主要为负值，且孔隙弹性效应增加了裂缝扩展的阻力。
-
导水断层是关键因素：只有在靠近注入井（30-100 m）且具有高渗透率（50-100 mD）的断层存在时，孔隙压力才能显著增加并到达深部基底。断层必须处于临界应力状态，微小的压力变化即可触发地震。
-
裂缝网络可传递压力：双重渗透率模型显示，裂缝网络（渗透率 5 mD）可将压力传递至 4 km 深度，但均质裂缝网络会导致压力同步增加，产生“群震”现象。
-
非均质性解释侧向迁移：沿裂缝网络的渗透率非均质性（从注入井附近的 10 mD 逐渐降低到远端的 1 mD）可导致压力前锋随时间逐步向外扩展，从而解释 Youngstown 地震震源从注入井附近向外侧向迁移的时空演化规律。
-
CMG GEM 的适用性：CMG GEM 能够有效模拟流体注入过程中的压力扩散、断层/裂缝网络中的优先流动路径、以及流体-地质力学耦合效应，是注入诱发地震风险评估的有力工具。
🏛️ 作者及单位信息
| 作者 | 单位 |
|---|---|
| S. Raziperchikolaei | Battelle Memorial Institute，哥伦布，俄亥俄州 |
| J.F. Miller | Battelle Memorial Institute，哥伦布，俄亥俄州 |
-
期刊：Interpretation (SEG/AAPG)
-
发表年份：不详（参考文献最新为 2014 年）
-
研究领域：注入诱发地震、流体流动模拟、地质力学、断层激活、裂缝网络
💡 补充说明
该论文是一项关于注入诱发地震机制的数值模拟研究，核心创新点在于：
-
多方法综合建模：整合了三维流体流动模拟（CMG-GEM）、水力压裂模拟和地质力学耦合模拟，系统评估了不同机制对深部压力传递的贡献。
-
断层/裂缝网络的系统参数化：系统改变了断层距离（30-800 m）、断层渗透率（10-100 mD）和裂缝网络渗透率（1-10 mD），定量分析了各因素对深部压力增加的影响，确定了导水断层是解释深部地震的关键条件。
-
双重介质模型的应用：使用 CMG-GEM 的双重渗透率模型模拟裂缝网络中的压力传递，为理解裂缝性基底中的流体流动提供了重要工具。
-
非均质性对地震侧向迁移的解释：首次提出沿裂缝网络的渗透率非均质性（从注入井向外递减）是导致地震震源随时间侧向迁移的可能机制，并通过数值模拟验证了这一假设（时间与空间演化与 Youngstown 实际地震序列一致）。
-
耦合地质力学分析：使用 CMG-GEM 的迭代耦合方法分析了流体注入对有效应力状态的影响，解释了为何微小压力变化足以触发地震（断层必须处于临界应力状态）。
-
对现场数据的历史拟合：使用实际注入速率数据和井口压力数据对模型进行了历史拟合，提高了模型的可靠性。
该研究对从事诱发地震风险评估、流体注入地质力学、CO₂ 封存安全性评估、深部咸水层注入、以及断层激活模拟的工程师和科研人员具有重要的参考价值。