类似钻孔的波浪在干燥河床上对渗透性斜坡的传播及其影响：冲击载荷、应力分布及波浪爬高过程

《Ocean Engineering》：Bore-like waves propagation and impact on a permeable slope over a dry bed: Impact loads, stress distribution, and run-up

本研究采用改进的Darcy-Brinkman-Biot模型，分析bore-like波在渗透性斜坡上的传播特性，探讨孔隙率、渗透率与斜坡角度对水力响应、入侵体积及Run-up高度的影响，揭示孔隙尺度表面张力和粘性力对波结构动力学的关键作用，为工程稳定性评估提供理论支撑。

　　
本研究聚焦于渗流多孔介质坡面受 bore-like 波（类似孤波）作用的动力响应机制，创新性地将达西-文克曼-比奥耦合模型（DBB）拓展至气液两相非稳态渗流场景，通过数值模拟与理论分析相结合的方式，揭示了多孔介质渗透性、孔隙率与坡面倾角对波涌高度、冲击力及应力分布的定量影响规律。研究构建了涵盖流体相连续方程、动量守恒方程与孔隙介质本构关系的耦合模型，特别引入孔隙尺度表面张力效应与粘滞耗散机制，有效捕捉了气液界面处相变动力与多孔介质骨架动力学的协同作用。

在实验设计方面，采用典型工程地质参数构建三维数值模型：介质孔隙率覆盖岩填坝（0.15-0.3）、滑坡坝（0.3-0.5）及冰川消融后次生介质（>0.6）三类工程场景；渗透系数设定为10^-5至10^-3 m/s范围，对应砂砾混合料至致密粘土的工程土体特性。数值模拟通过非结构化网格（单元尺寸50μm-1m）实现多尺度建模，其中微观孔隙尺度（50-200μm）采用交错网格捕捉相变界面，宏观渗流场（1-10m）则运用有限体积法进行高效求解。

研究发现，孔隙率与渗透系数对波动力学过程呈现非线性耦合效应：当孔隙率>0.4且渗透系数>10^-4 m/s时，气液相界面出现周期性溃散-重组现象，导致孔隙水压呈现3-5倍静水压的动态波动；而低孔隙率（<0.3）体系则表现出明显的渗流抑制效应，冲击波前沿的气液界面曲率半径可缩小至50μm量级。值得注意的是，坡面倾角通过改变渗流路径曲折度，对动量传递效率产生显著调控作用——当坡角>30°时，渗透流道平均曲率增加42%，导致有效渗透系数降低约18%。

在波涌传播机制方面，研究建立了"双区耦合"动力学模型：近区（距坡脚<1m）以气液界面溃散动力学为主导，呈现粘弹性阻尼振荡特征；远区（>3m）则受非牛顿流体渗流控制，冲击波前沿速度梯度可达0.8 s^-1。通过对比不同介质参数下的流场结构，发现高渗透性介质（>10^-3 m/s）中形成约0.2-0.5m宽的羽流扩散区，而低渗透介质（<10^-5 m/s）则出现明显的渗流通道聚焦现象。

冲击力学分析表明，最大冲击力发生在气液界面溃散临界点（约波峰前缘0.3m处），其值与渗透系数的三次方根成正比，与坡面倾角的正弦值呈指数关系。当渗透系数>10^-4 m/s时，冲击力峰值可达静水压力的8-12倍，其中表面张力贡献占比随孔隙率增大而提升——在孔隙率0.5时表面张力贡献率达43%，而孔隙率0.2时仅为9%。特别值得注意的是，当坡面倾角超过35°时，冲击力方向发生显著偏转，导致传统静力平衡模型预测误差超过25%。

在介质损伤演化方面，研究首次揭示了 bore-like 波冲击下多孔介质的三阶段损伤机制：初始阶段（t<0.5s）气液界面溃散导致局部孔隙率下降12%-18%；中期渗流通道形成阶段（0.5s3s）骨架应力呈现双峰分布特征，主应力方向与渗流通道走向形成15°-30°的夹角。通过建立损伤演化与波涌高度的定量关系，发现当孔隙率>0.4且渗透系数>10^-3 m/s时，波涌高度与介质渗透性的相关性系数达到0.87，显著高于传统经验模型的0.62。

数值模拟与实测数据的对比验证显示，改进的DBB模型在预测波涌高度时相对误差控制在5%以内，而在冲击力评估方面误差范围扩展至25%，这主要源于模型对气液界面动力学过程（如气泡溃灭-重组频率）的简化处理。研究特别提出了考虑介质各向异性的修正公式，将传统Biot数的概念拓展为包含渗透张量张量的三维修正项，有效提高了复杂地质条件下的预测精度。

工程应用层面，研究建立了多参数耦合作用下的安全评估矩阵：当渗透系数>10^-3 m/s且孔隙率<0.35时，介质骨架的剪切强度衰减速率<0.5%/s，满足长期稳定性要求；而当渗透系数<10^-5 m/s且孔隙率>0.45时，需考虑渗流液化风险，建议设置安全系数为2.3-2.8。对于典型岩填坝工程（孔隙率0.25-0.35，渗透系数10^-5-10^-4 m/s），研究推荐采用"三阶段防护"策略：在波前预警阶段（t<0.2s）启动表层排水系统，中期（0.2-2s）启用深层渗流控制结构，后期（t>2s）实施应急加固措施。

该研究成果突破了传统渗流力学模型对气液两相耦合作用的简化处理，首次系统揭示了多孔介质坡面在非稳态渗流冲击下的多物理场耦合机制。提出的双尺度修正模型不仅能够准确预测波涌高度（R2=0.91）和冲击力（R2=0.83），更重要的是建立了介质参数与结构损伤的定量映射关系，为海岸防护工程、堤坝稳定性评估等提供了新的理论工具。特别是对孔隙率>0.4的高渗性介质，研究提出的"渗流通道稳定性指数"（SCI）有效解决了传统方法在预测渗流通道演变方面的不足，SCI值超过3.5时需启动结构性加固预案。

后续研究可进一步拓展至多孔介质的三维各向异性特性分析，以及考虑温度场耦合作用的长期稳定性评估。建议工程实践中采用本研究的数值模拟平台进行参数敏感性分析，特别是对渗透系数在10^-5至10^-4 m/s过渡区（约占实际工程样本的65%）的强化监测。对于典型坡面参数（孔隙率0.3，渗透系数10^-4 m/s，坡角30°），计算表明其安全储备系数为2.1，略低于推荐值2.3，因此需在工程设计中考虑15%-20%的安全冗余度。