Revisiting the hydromechanical formulation of a micromechanics-based phase-field model for poro-elastoplastic media

在地球深处，在蕴藏着石油、天然气和地热能的多孔岩石中，一场固态岩石与其中捕获的流体之间的无声战斗正在不断进行。当工程师试图开采这些资源时，他们经常向岩石中注入流体以使其撑开，这一过程被称为水力压裂。几十年来，科学家们一直将这一过程建模为将岩石视为一种简单的、具有弹性且能拉伸并回弹的“海绵”。然而，真实的岩石更为复杂；它们还可以发生永久性的弯曲、破碎，并沿着内部裂缝滑动，尤其是在受到高压挤压时。理解这些岩石在固体骨架和内部流体同时进行推拉作用时是如何破裂的，对于安全钻井、管理废物存储，甚至预测自然地质事件（如冰山开裂或岩浆穿透地壳）至关重要。挑战在于，如何创建一个既能捕捉岩石破裂这种混乱现实，又不会陷入无法计算的复杂程度的数学描述。

一个研究小组重新审视了用于描述这一破裂过程的基础方程，发现以往模型中常用的一种简化方法会导致显著误差。在岩石力学领域，科学家使用一种“相场”（phase-field）方法，将裂纹视为一个模糊的区域，即材料逐渐弱化的区域，而非一条尖锐、隐形的线。这使得计算机能够模拟裂纹如何生长和改变形状。为此，模型必须考虑两种截然不同的行为：当岩石被拉开时（张应力）和被挤压时（压应力）。在张应力下，微小的裂缝会张开；在压应力下，它们闭合，且内部粗糙的表面会相互滑动并产生摩擦。研究人员发现，许多广泛应用于工业界和学术界的现有模型混淆了描述系统能量的变量。具体而言，他们发现将流体压力直接作为能量方程的输入，而不是将孔隙中的流体量作为输入，会导致数学上的“跳跃”或不连续性。这种跳跃意味着，当模型从受拉转为受压时，预测的岩石强度会发生不切实际的突变，从而导致对裂纹形成方式和位置的预测出现偏差。

为了解决这个问题，研究人员开发了一种新的、一致的公式，仔细追踪了岩石变形与流体含量之间的关系。他们基于微观视角构建了模型，将岩石想象成由充满流体的微小、钱币状裂缝组成的集合。通过分析这些微裂缝在不同应力下如何开启和闭合，他们推导出一套规则，确保岩石的强度表面（即定义其何时破裂的数学边界）无论压力如何变化，都能保持平滑和连续。他们证明，通过使用一种特定的流动法则，将岩石的永久变形直接与作用力联系起来，并正确包含能量方程中缺失的一个耦合项，该模型可以准确描述从张应力到压应力的过渡。这一修正不仅仅是一个微小的改进；它从根本上改变了模型计算推动裂纹前进的驱动力的方式。

团队通过测试两种截然不同的场景来验证新模型的有效性。首先，他们模拟了一个经典的液力压裂基准实验，即向岩石中注入流体以产生裂纹。他们将结果与一个已知的解析解（即针对理想条件推导出的精确数学答案）进行了对比。新模型比以往那些混合了变量的模型更接近该解析解，能够正确预测启动裂纹所需的压力以及裂纹生长时的宽度。旧模型往往会高估所需的流体压力，这一缺陷可能导致钻井效率低下或存在安全风险。在第二个测试中，他们模拟了一个从两侧挤压岩石样本的过程，这种条件模拟了深层地下环境中岩石承受巨大压力的状态。在这种情况下，岩石并非仅仅是突然断裂，而是发生了塑性变形，形成了一个剪切带，材料在破裂前发生了相互滑动。新模型成功重现了这种复杂的行为，既捕捉到了由机械挤压引起的剪切裂纹，也捕捉到了由流体注入驱动的张力裂纹。

模拟结果表明，岩石的塑性变形是裂纹扩展的一个重要障碍。当岩石允许发生永久变形时，它会吸收更多能量，这意味着裂纹内部的流体压力必须更高才能保持裂纹开启并持续生长。这种额外的阻力改变了整个裂缝的剖面，使其比纯弹性模型预测的更窄，且需要更大的力量才能传播。此外，在压缩测试中，模型展示了岩石内的孔隙压力如何随裂纹形成而变化。随着岩石变形并为流体移动创造空间，局部压力会下降，从而吸引更多流体进入受损区域。这种相互作用产生了一个反馈循环，即流体流动与岩石机械失效之间是紧密耦合的，而新模型能够高保真地捕捉到这一动态过程。

通过修正这些模型的数学基础，研究人员为从事多孔介质工作的工程师和科学家提供了一个更可靠的工具。这种新公式确保了不同类型岩石失效之间的过渡是平滑且符合物理实际的，消除了困扰以往尝试的人为跳跃。这种准确性对于从优化油气生产到设计碳捕集安全地下存储等应用领域都至关重要。这项研究表明，即使在固体与流体相互作用的复杂系统中，通过对底层物理机制进行严谨、一致的处理，也能揭示材料破裂的真实本质。该模型现在已准备好以更高的信心帮助预测地壳的行为，弥合了理论力学与脚下那片混乱、具有塑性的真实大地之间的鸿沟。