裂隙岩体流固耦合的非线性渗流力学模型与分岔分析

裂隙岩中隐匿水体为何重要

在我们脚下深处，水通过矿井、隧道和油气储层中的破碎岩层流动。水对岩体的作用——以及岩体对水的反作用——可能在不经意间为突发性的洪流、塌方或地震埋下隐患。本文提出了一种新的描述与预测裂隙岩体中水流与岩体变形双向耦合的方法，揭示了为何这类系统可长时间看似平稳却会突然失稳。

水在破碎岩石中的穿行

在许多地下工程中，岩体并非像砖块那样致密，而更像破裂的海绵：无数裂隙和孔隙构成了水流的通路。来自上覆岩体的外部应力和孔隙内部水的压力不断重塑这些通道。当岩体受压，孔隙会收缩或闭合；当孔隙水压上升，孔隙可能重新打开或扩展。这种持续的重塑改变了水的易流性，进而又反馈影响岩体。作者认为，要理解矿山涌水或水库泄漏等灾害，必须把它视为一个动态耦合系统，而非静态快照。

构建岩体与水的耦合图景

研究首先扩展了土壤力学中的经典概念“有效应力”，该概念描述了总应力中实际由固体骨架承担的那一部分。作者将这一思想重写为明确包含孔隙率——即岩体体积中空隙所占的比例——使孔隙空间的变化能够直接与岩水之间应力分担相联系。随后，他们把这一定义与描述微可变形岩体弹性响应的方程相结合，并引入了对裂隙渗流更为真实的非线性描述，超越了许多工程模型中常用的线性达西定律。

从平滑流动到突变

在这一框架下，作者关注了水沿垂直方向通过一层破碎岩体渗流的单维情形。他们推导出一组非线性方程，描述在孔隙随岩体压实而调整时，水压和流量如何随时空演化。求解这些方程表明，在某些条件下系统并不存在唯一的稳态：相反，它表现出数学上所称的鞍结分岔。通俗地说，当某个关键流动参数变化时，原先稳定的状态可能分裂为一条安全分支和一条危险分支，或完全消失，导致系统从温和渗流突然跃迁到失控流动。

缓慢挤压与延迟的稳定性

作者接着考察了当边界应力（例如由上方采矿引起的逐步加载）随长时间尺度变化时情形的变化。数值模拟表明，当外部应力缓慢变化时，耦合的岩—水系统也需要更长时间才能趋于稳态。水压、流量和岩体体积应变呈现蠕变式地向稳定态靠拢，而非迅速平稳下来。这种延迟源于岩体骨架在不断被改变的载荷下必须持续调整孔隙结构，而能量不断被输入系统，从而拉长了达到平衡的过程。

洪水前的预警信号

为将理论与实际连接，研究将其预测与一起煤矿断层涌水的真实案例进行了比对。随着开采接近断层，一个反映流动偏离线性达西行为程度的参数进入了一个临界区间，此时两种流态可以共存：一种稳定，一种不稳定。现场测量显示水速在两种不同水平之间开始波动，随后急剧上升为快速、灾难性的增大，正如模型的分岔图所示。作者认为，这类波动比把系统视为线性与稳态的传统安全指标更早、更明确地预示着即将发生的涌水风险。

对地下安全的意义

总体而言，本文表明饱和水的裂隙岩体更像一个复杂的非线性系统，而非简单的管道。应力或流动条件的微小变化可能把系统推过临界阈值，从而在行为上发生质的改变，而不仅仅是幅度的变化。通过明确连接岩体变形、孔隙结构与非线性流动，新模型能够捕捉多重可能的稳态、它们之间的突变转移以及对初始条件的强烈敏感性。对设计矿山、隧道与水库的工程师而言，这意味着监测流动与变形随时间的演化——并注意那些预示双重稳态的波动——可以在隐匿不稳变成灾难之前提供更早、更可靠的预警。

引用: Zhengzheng, C., Mengqi, X., Tao, R. et al. Nonlinear seepage mechanical model and bifurcation analysis for fluid-solid coupling in fractured rock mass. Sci Rep 16, 9578 (2026). https://doi.org/10.1038/s41598-025-25823-6

关键词: 裂隙岩体, 地下水渗流, 流固耦合, 非线性动力学, 矿山涌水