Injection-rate effects on failure in a fluid-saturated granular fault gouge

本文结合解析理论与数值模拟，证明了流体注入速率通过产生压力非均匀性来控制断层注屑失效，其中缓慢注入会导致均匀弱化，而快速注入则能保持远端区域的强度，从而为预测地质工程作业中的地震活动性提供了一个精细化的框架。

要点

想象一下，地球的地壳就像一条巨大的、布满裂缝的人行道。在这些裂缝中，并不是空无一物，而是挤满了破碎的岩石、沙子和泥土。地质学家称之为“断层碎屑”（fault gouge）。现在，想象有人开始向这个裂缝中泵入水。这是地热能开发或废物处置等领域的一种常见做法，但它也带来了一个风险：水压可能会将裂缝撑开，导致地面滑动并引发地震。

这篇论文提出了一个简单但棘手的问题：注水速度快慢是否重要？

研究人员发现，答案是肯定的，而且非常重要。如果你注水很慢，水会均匀地扩散开来，裂缝很容易发生滑动。但如果你注水非常快，你实际上需要更高的压力才能让裂缝发生滑动。

他们是如何得出这一结论的呢？他们结合了数学模型和计算机模拟。

“拥挤房间”类比

把断层碎屑（破碎的岩石）想象成一个挤满了人（颗粒）的拥挤房间。

• 目标： 你想让所有人移动到侧边（滑动）。
• 水： 水压就像是一个试图让人们分开的“推力”。
• 膨胀： 当人们尝试移动时，他们自然会散开并占据更多空间（这被称为“剪胀”现象）。

两种情景

1. 慢注（慢速注入）
想象非常缓慢地向房间里注水。水有充足的时间渗透进人群，并到达每一个人身边。压力变得均匀一致。每个人同时感受到推力，整个人群一起松动，房间很容易发生滑动。这就是旧理论所预测的情况：水压越大 = 越容易滑动。

2. 高压水枪（快速注入）
现在，想象从一侧以极高的速度向房间里喷射水流。

• 梯度： 紧挨着水枪的人立刻被淋湿并受到强力推动。但房间另一端的人呢？他们仍然是干燥且站稳的。
• 瓶颈： 即使靠近水枪的人已经准备好移动，但远端的人仍在坚守阵地。只有当最远端的人也被推力覆盖时，整个房间才能发生滑动。
• 结果： 为了让整个房间发生滑动，你必须把水枪端的压力调到比慢速注水时高得多的水平。快速注入创造了一个“压力梯度”，即推力在某处很强，而在另一处很弱。

“松散沙子”效应

这里还有第二个转折。当岩石颗粒尝试滑动时，它们不仅仅是滑动，还会发生碰撞和重新排列，使岩层稍微变厚（剪胀）。

• 在计算机模拟中，研究人员逐步增加压力，让颗粒在每一步之后都重新稳定下来。
• 他们发现，随着颗粒重新排列以及岩层变得“松散”，材料实际上变得更弱了。
• 然而，由于快速注入创造了那些不均匀的压力区（靠近水枪处强，远端处弱），这种“松散沙子”带来的弱化效应，直到压力高到足以克服远端那些强固、干燥的区域时，才能帮助整个断层发生滑动。

数学“配方”

作者创建了一个新的数学公式来预测断层何时会发生滑动。该公式指出，引起滑动的压力取决于三个因素：

1. 注水速度： 注水越快 = 需要的压力越高。
2. 断层长度： 断层越长 = 需要的压力越高（因为压力必须传播得更远才能到达薄弱点）。
3. 水在岩石中的流动速度： 如果岩石的孔隙率很高（水流动很快），压力会迅速平衡，因此需要的压力较小。

这为什么重要（根据论文观点）

论文结论认为，我们不能仅仅使用那些假设水压处处相等的旧的、简单的规则。

• 权衡关系： 如果你注水非常快，你可能需要更高的压力来触发滑动，这听起来似乎更安全。然而，由于你注水很快，你可能会在更短的时间内达到那个危险的压力阈值。
• 规模因素： 断层的长度是一个巨大的影响因素。短断层表现得像“慢注”（压力均匀），而长断层则表现得像“高压水枪”（压力不均），这使得预测其何时滑动变得更加困难。

简而言之，论文表明速度改变了规则。快速注水会创造出不均匀的压力区，这些区域起到了“维持现状”的作用，需要比慢速注水更大的力量才能打破断层。这有助于工程师理解，断层的尺寸和注水速度是预防或预测诱发性地震的关键因素。

技术摘要

技术摘要：流体饱和颗粒状断层 Gouge 中注入速率对失效的影响

问题陈述
在能源开采、废物处置和资源开发过程中，向地下注入流体是诱发地震性的已知驱动因素，这往往会重新激活充填有 gouge（断层泥）的断层。虽然通过孔隙压力扩散降低有效正应力是典型的弱化机制（Terzaghi 有效应力定律），但注入速率在控制失效中的作用仍不为人所熟知。经典模型通常假设孔隙压力均匀，但最近的研究表明，异质的压力分布可以显著改变失效阈值。具体而言，快速注入可能会产生强烈的压力梯度，使远端断层区域保持强度并延迟失效；而缓慢注入则允许均匀的弱化。此外，与由水力连通性主导的裸露表面断层不同，充填有 gouge 的断层涉及复杂的颗粒响应，包括失效前阶段的扩张（体积膨胀）。这种扩张具有竞争效应：它既可能导致瞬态孔隙压力下降（扩张硬化），也可能在颗粒重排的过程中导致堆积密度和剪切强度的降低（弱化）。在流体饱和的 gouge 中，注入速率、压力扩散与颗粒扩张之间的相互作用需要一个定量框架。

研究方法
作者采用结合解析理论与耦合数值模拟的双重方法：

1. 数值模拟： 使用耦合固–液代码，集成用于模拟固体颗粒（建模为具有线性弹性摩擦接触的二维圆盘）的离散元法（DEM）以及用于连续孔隙流体的欧拉网格有限差分求解器。实验设置模拟了一个受压、流体饱和的颗粒层，代表断层 gouge，该层被限制在不透水壁之间。流体在侧向边界以受控速率（20–200 kPa/min）进行注入，并采用离散压力阶跃的方式，以实现加压与扩张之间的平衡，从而防止瞬态压力下降。模拟在不同的预应力水平 ($\tau/\sigma_n$) 和归一化层长度 ($L/d$) 下运行。
2. 解析理论： 基于包含“扩张汇”（dilative sink）项的孔隙压力扩散方程推导出一个解析模型。该项解释了颗粒骨架在接近失效时的体积膨胀。模型假设扩张速率与加压速率呈线性关系。该扩散方程在空间均匀边界条件和汇项的作用下，预测了跨越断层的孔隙压力剖面演化。

主要结果

• 失效压力的速率依赖性： 数值模拟表明，触发断层失效所需的边界压力 ($P^b_{fail}$) 随注入速率的增加而系统性上升。这证实了快速注入相比于缓慢注入会延迟失效。
• 异质压力分布： 解析解揭示了快速注入会产生强烈的空间压力梯度。压力在靠近注入边界处最高，在断层中心 ($x=L/2$) 最低。因此，全局失效并非由平均压力或边界压力控制，而是由层中最压力最低（最强）的区域控制。
• 定标律： 解析模型预测，失效压力与注入速率 ($\dot{P}_b$) 成线性比例，与断层长度 ($L^2$) 成二次方比例，同时与水力扩散系数 ($\alpha$) 成反比，并取决于有效体积模量与流体存储能力之积 ($K\beta$)。
• 扩张的作用： 模拟显示，扩张随压力增加而增加，遵循幂律定标 ($\Delta h/h_0 \sim (\Delta P^b/\sigma_n)^{3/2}$。随着层发生扩张，有效体积模量 ($K$) 降低，反映了颗粒堆积的渐进弱化。解析模型通过将 $K$ 视为有效常数来限定模拟结果，尽管模拟表现出亚线性行为，这是由于随着扩张增加，摩擦系数逐渐降低所致。
• 验证： 解析预测在三个数量级的颗粒堆积刚度 ($K\beta$) 范围内成功界定了数值结果。该模型重现了经典均匀压力有效应力理论无法捕捉到的实验观察结果。

意义与主张
本文声称提供了一个将颗粒尺度物理学（颗粒扩张与重排）与流体注入下的断层尺度失效联系起来的定量框架。通过推导出速率依赖的失效准则，本研究解释了为什么在注入 gouge 时，更高的注入速率需要更高的压力才能触发重新激活。

作者强调，这种速率依赖性源于扩散控制的孔隙压力异质性。在缓慢注入的情景下，压力趋于平衡，导致均匀弱化。在快速注入的情景下，强烈的梯度持续存在，断层在中心（压力最低处）达到失效阈值之前保持稳定。

关于实际应用的影响，论文指出操作人员可以通过控制注入方案来潜在地管理重新激活风险。虽然快速注入提高了失效压力阈值，但它也缩短了失效时间。作者指出，这种速率–压力权衡在特定操作背景下（如循环注入）可能具有优势，但他们明确警告说，实际应用必须考虑自然断层的异质性、复杂的几何形状以及热效应，这些在他们的二维模型中被简化或忽略了。研究结论指出，如果不考虑扩散控制的定标效应（特别是关于断层长度的影响），从实验室尺度得出的失效准则不能直接推广到现场尺度。