基于细观有限元分析的砂样等效渗透性预测.pdfVIP

2023年6月30日铁路地质与路基2023年第1期总第105期

基于细观有限元分析的砂样等效渗透性预测

王梦辉

（中铁第四勘察设计院集团有限公司，湖北武汉430063）

【摘要］砂土的宏观渗透性和细观孔隙水平的真实流速对设计反滤层具有重要意义。

采用有限元方法离散不可压缩流体流动的纳维一斯托克斯方程，采用不可压缩单元来模拟不

可压缩流体，引入周期性边界条件（periodicboundaryconditions,PBC）来降低模型尺寸大小

的影响。采用该方法计算了X射线微层析成像（X-raymicrotomography，XMT）技术重建

的砂样和随机生成的砂样模型的渗透性，并与实测值和经验值进行了对比。结果表明：在砂

样渗流过程中，孔隙流体绕过固体颗粒流动，细观孔隙流体流速分布并不均匀，而是存在优

势水力路径。

力等效渗透性

【关键字］有限元分析表征体元周期性边界条件孔隙流体流动

在岩土工程中，砂土的宏观渗透性和细流动特性仍符合经典流体力学理论2]，只要

观孔隙水平的真实流速是设计反滤层的重要求解不可压缩粘性流体流动的纳维－斯托克

因素，反滤层对防止土体发生流土和管涌破斯方程便可得到细观孔隙流体流动的规律。

坏具有重要意义。宏观分析通常用于研究土大量学者采用数值方法求解不可压缩

体中的平均孔隙压力和流体速度。然而，仅粘性流体流动的纳维－斯托克斯方程，研

仅通过宏观研究无法获得渗流过程中孔隙细究复杂形态孔隙中细观流体的流动。Karim

观层面上某个具体点的流体流速和压强，关等[3]采用XMT技术重建了砂样的三维真

于渗流机理的研究也需要对细观层面上的流实微结构，并采用有限元法计算其渗透

体流动进行详细地分析。细观层面上，砂率。Pazdniakou等[4采用LBM对平行板间

样是天然的异质材料，其主要由较密实的Poiseuille流动进行了验证，得到了与解析

固体颗粒和连通孔隙组成。在梯度压力作解一致的结果；计算了重构和真实多孔介质

用下，流体主要从连通的孔隙通过，对孔（砂岩）在较大频率范围内的动态渗透率，

隙中流体流动规律的研究是了解砂样宏观渗验证了普遍的标度行为，实测值与理论预测

透性的基础。有很好的一致性。咎鹏等[5]采用LBM的有

砂样是一种典型的多孔介质。多孔介质限体积形式计算了完全多孔Poiseuille流和

在宏观尺度上常表现为连续介质，特征尺寸Couette流，数值计算结果与理论解吻合较

在mm量级及以上，细观尺度表现为孔隙和好。Bignonnet等[提出了一种计算复杂微

固体介质两种材料组成的异质结构，特征尺结构渗透率的快速傅里叶变换（fastFourier

寸在m～mm量级。多孔介质内部孔隙流transform，FFT）方法，该方法直接与多孔

体的流动特性涉及到多孔介质的细观结构以介质样品的3D成像技术相耦合，不需要

及细观尺度流体的运动规律。基于此，众多网格划分和定义等效孔隙网络，大大缩短

学者在这方面做了相关研究，研究结果表明了计算时间。田志国等[7采用离散单元法

在细观孔隙层次上多孔介质孔隙内部流体的