多孔岩石力学元模型：孔隙度与水弱化对砂岩脆-塑性转变及变形模式的定量预测

《INTERNATIONAL JOURNAL OF ROCK MECHANICS AND MINING SCIENCES》：Metamodel for porous rock: Porosity and water-weakening controls on brittle–ductile transition in sandstones

针对砂岩力学参数（强度、屈服面形状、硬化/软化规律）随孔隙度与含水状态变化的量化难题，研究者基于改进剑桥粘土模型(Modified Cam Clay)与分岔分析，构建了以初始孔隙度和水饱和度为变量的砂岩力学“元模型”(metamodel)。结果表明该模型能高精度预测Py、ρs等参数的幂律关系及应变局部化模式，为地质储层（如碳封存、地热）的变形预测提供了高效工具。

　　

在地壳浅部至数公里深度的岩石圈中，砂岩作为最常见的沉积岩之一，广泛分布于各类地质工程中，从水利水电的大坝地基、交通隧洞，到当今能源领域的地热开采井筒和地质碳封存（CCS）储层。然而，这些工程岩体并非均质的“铁板一块”，它们内部布满了微米至毫米级的孔隙网络。当工程师试图预测这些岩体在开挖或注入流体后的变形与稳定性时，往往会陷入困境：岩石的力学行为不仅随深度（围压）发生复杂的脆-塑性转变（Brittle–Ductile Transition），还极度敏感于岩石的初始孔隙度（Porosity）以及孔隙中流体（尤其是水）的存在——后者常导致显著的“水弱化”（Water-weakening）现象，使强度降低10%–50%。

传统上，为了精确模拟一个特定砂岩地点的力学响应，工程师需要进行繁琐的实验室三轴试验来校准复杂的本构模型参数，这在实际操作中耗时费力，且难以应对储层纵向和横向上因沉积环境不同而产生的巨大孔隙度与饱和度变化。那么，是否存在一种更聪明的办法？能否跳过繁复的单点校准，直接通过建立“元模型”（Metamodel，即一种能根据简易可测参数自动生成本构模型的数学结构），来实现对砂岩力学行为的快速、大尺度预测？这正是Ghassan Shahin与Marie Violay（洛桑联邦理工学院EPFL岩石力学实验室）在《International Journal of Rock Mechanics and Mining Sciences》上发表的核心议题。

为解答上述问题，研究人员开展了一项系统性的理论与数据融合研究。他们收集了超过15种天然砂岩（如Bentheim, Berea, Darley Dale等）和三组合成砂岩（烧结玻璃微珠，控制变量）广泛的排水三轴压缩实验数据。研究采用了一个基于临界状态土力学框架的多变量弹塑性（Elastoplastic）模型，引入两个内部变量分别表征由剪胀/破裂导致的软化机制和由孔隙压塌导致的硬化机制。通过对所有岩石进行统一的本构校准与分岔分析（Bifurcation analysis，用于预测变形带失稳），团队系统提取了模型参数（如静水屈服压力P_y、形状参数M_f/M_g、硬化参数ρ_s、软化参数ρ_m等），并分析了这些参数随孔隙度（φ, 约10%–30%）和水饱和度变化的敏感性规律，最终集成为适用于工程尺度的砂岩元模型。

2. Laboratory measurements（实验室测量）

基于已发表的多种砂岩（Adamswiller, Bentheim, Berea, Rothbach, Bleurswiller, Darley Dale, Gasford及8种合成砂岩）在三轴压缩下的应力-应变、体变及变形模式（剪切带、压实带）数据，研究人员建立了涵盖孔隙度10%–30%、粒径分布各异及干/湿两种状态的数据库。这些数据为后续统一的本构校准提供了基准，特别是湿态数据中观察到的强度降低与 compaction banding（压实带）现象，为量化水弱化效应提供了依据。

3. Constitutive modeling（本构建模）

3.1. Elastoplastic mathematical formalism（弹塑性数学形式）

研究采用了改进的剑桥粘土模型（Modified Cam Clay）形式的屈服面与塑性势面，公式表示为 f, g: q²+ M_h²(p - P_y)p = 0。模型引入了两个硬化/软化内变量：P_s（压硬性，由体积收缩引起）和 P_m（结构性退化，由剪胀和剪切引起），总屈服压力 P_y= P_s+ P_m。演化法则区分了剪胀（Dilation）与压缩（Contraction）状态下的塑性应变率，能够分别捕捉低围压下的应变软化（脆性）与高围压下的应变硬化（延性）。

3.2. Calibration and validation of constitutive models（本构模型的校准与验证）

通过曲线拟合实验的应力-应变曲线确定屈服点与初始刚度，并利用剪胀函数（d = dε_v^p/dε_s^p）校准塑性势参数。结合分岔分析（基于声学张量行列式 det(A(θ)) 的正负判断失稳），模型成功复现了不同砂岩在低围压下的剪胀剪切带、中围压下的剪切增强压实带以及高围压下的纯压实带及其倾角变化，与实验室观察高度一致。

4. Generalization of plasticity parameters and metamodeling（塑性参数泛化与元建模）

4.1. Sensitivity analysis and plasticity generalization（敏感性分析与塑性泛化）

参数敏感性分析揭示了关键规律：1. 静水屈服压力 P_y和硬化参数 ρ_s强烈依赖于孔隙度 φ，分别服从指数为 -3/2 和 -3 的幂律关系（ρ_s= 0.12 × φ^-3）；2. 屈服面形状参数 M_f（~1.1）、塑性势形状参数 M_g（~1.6）及软化参数 ρ_m（~450）对 φ 和饱和度变化不敏感，可视为通用常数；3. 水饱和度的增加成比例的降低体积模量 K 和 P_y，但几乎不改变形状与硬化/软化规律。

4.2. Generalizing plasticity parameters and metamodeling（泛化塑性参数与元建模）

基于上述规律，构建了砂岩元模型。该模型仅需输入初始孔隙度 φ 和刚度（或 P_y参考值），即可通过内置的幂律与比例关系，自动输出全套本构参数（K, G=2/3K, M_f, M_g, ρ_m, ρ_s, P_y），无需针对每种岩石单独复杂校准。

4.3. Application scenarios（应用场景）

验证了三个典型场景：1. 简化校准——仅需单轴或 hydrostatic 试验获取 K 和 P_y，其余参数由 φ 决定；2. 干-湿转换——已知干岩参数时，仅需湿岩刚度比即可预测湿岩全参数；3. 孔隙度外推——基于一种孔隙度的校准数据，预测同岩性不同孔隙度（10%–30%）的力学行为。模型在应力-应变、体变及应变局部化预测上均表现出与详细校准模型相当的精度。

5. Discussion（讨论）

5.1. Porosity controls on mechanical transitions beyond strength（孔隙度对强度之外力学转变的控制）

研究指出，孔隙度不仅是强度的主控因素，更通过控制 cataclastic pore collapse（碎裂性孔隙压塌）的效率，主导了材料从应变软化（脆性、剪胀、局部化）向应变硬化（延性、压缩、扩散化）的转变。硬化参数 ρ_s与 φ 的负相关关系，被解释为高孔隙度岩石中单位塑性压缩应变消除空隙空间的效率较低。

5.2. Porosity scaling and the role of power-law extrapolation（孔隙度标度与幂律外推的作用）

P_y∝ φ^-3/2的标度关系与 Wong 等人基于赫兹断裂力学（Hertzian fracture mechanics）的理论一致。利用粒径一致、仅 φ 变化的合成砂岩数据，强力验证了该幂律在多孔砂岩中的稳健性。

5.3. Constitutive modeling versus micromechanics of water-weakening（水弱化的本构模拟与微观力学）

水弱化并未改变塑性势（即 dε_v^p/dε_s^p关系），说明水并未引入新的变形机制，而是通过降低敏感矿物（如粘土、长石）的刚度与接触强度，使屈服面收缩（P_y降低），从而在相同应力路径下更早进入塑性，表现出更强的压缩或强度降低。

5.4. Lessons learned from constitutive modeling of sandstones（砂岩本构建模的经验）

简单的双参数（M_f, M_g）椭圆屈服面足以捕捉主要力学响应，无需过度复杂的多参数模型，只要 P_y和 φ 被合理约束。这为工程应用提供了简化的路径。

5.5. Metamodels bridge reservoir geophysical imaging and borehole sampling（元模型衔接储层地球物理成像与钻孔取样）

5.5.1. Deterministic simulations enabled by metamodel bridging reservoir images and borehole characterization（通过元模型衔接储层图像与钻孔表征实现的确定性模拟）

元模型最大的工程价值在于桥梁作用：利用钻孔岩心标定出参考 P_y与关系式，再将储层地震/测井反演得到的空间孔隙度 φ 与刚度场输入元模型，即可为大尺度有限元模拟（FEA）提供具有空间变异性的本构参数场，从而实现从“单点”到“全场”的储层变形预测。

5.5.2. Uncertainty quantification with mechanically informed bounds（基于力学信息界限的不确定性量化）

基于参数敏感性分析得到的散射范围（如 ρ_s和 P_y的幂律上下界），可为概率模拟（如蒙特卡洛方法）提供力学的先验分布界限，而非随意假设，提升了储层预测的不确定性量化（UQ）的可靠性。

5.6. Limitations and future perspectives（局限性与未来展望）

作者指出当前元模型未显式考虑粒径（虽在有限粒径范围内影响次要）、力学各向异性（如层理）、速率效应及温-化-力全耦合。未来的验证还需更多单一岩性、宽孔隙度范围（如 Fontainebleau 砂岩）的数据，并借助原位 X 射线断层扫描（X-ray tomography）获取定量的应变局部化数据以优化剪胀函数。

6. Conclusions（结论）

该研究成功量化了孔隙度与水饱和度对砂岩脆-塑性转变及本构参数的控制作用。证实静水屈服压力 P_y和硬化参数 ρ_s分别遵循 φ^-3/2和 φ^-3的幂律；而形状参数、软化参数及水饱和度主要影响 P_y和刚度（成比例降低）。基于此构建的砂岩元模型，能够仅依据初始孔隙度（10%–30%）和饱和度，准确生成预测应力-应变、体变及变形模式的本构参数。这一成果极大简化了异质性储层（如地质碳封存与地热开发）中大尺度地质力学模拟的校准流程，推动了钻孔数据与地球物理成像在变形预测中的系统化整合。