渗流

定义

渗流指流体在多孔介质中的流动过程。多孔介质是内部含大量微小孔隙或裂隙、可供流体通过的固体介质，天然岩石、土壤以及陶瓷、过滤材料等人工材料都属于这一范畴^[3]。材料内部空隙尺寸微小、比表面积大，使渗流表现为阻力较大、流速缓慢、毛细管作用和表面分子力作用突出的流动^[4]。

研究渗流的学科称为渗流力学，它是流体力学的一个分支，研究流体在多孔介质内的运动规律及其工程应用^[5]。由于早期研究多集中于地下岩土层，渗流一词曾长期被狭义理解为地下水流动，随着研究扩展到过滤器、填充床和生物组织，这一狭义用法逐渐减少^[4]。

原理

渗流的基本定律是达西定律。1856 年，法国工程师亨利·达西经由直立圆筒砂柱渗流实验发现，通过多孔介质的流量与过水断面面积、水头损失成正比，与渗流路径长度成反比，写成速度为 $v=kJ$，其中 $v$ 为渗流速度，$J$ 为水力梯度，$k$ 为渗透系数^[1][6]。以压力梯度表达的通用形式为 $v=-\frac{k}{\mu}\nabla p$，其中 $\mu$ 为流体黏度^[7]。

达西定律中的渗流速度是按整个断面（含固体骨架）计算的表观流速，并非流体质点在孔隙内的真实速度。孔隙内平均实际流速等于渗流速度除以孔隙度，即 $u=v/n$，因此实际流速总是大于按全断面折算的流速^[3][8]。渗透系数综合反映了介质结构与流体性质，渗透率则只与介质本身有关，与通过的流体种类无关^[3]。

对复杂的渗流场，通常依据质量守恒、达西定律和状态方程建立连续性方程与渗流偏微分方程，再结合初始条件和边界条件求解，以得到压力与速度分布^[9]。当存在多种流体时，还需引入相渗透率和毛细管压力等概念描述多相渗流^[2]。

发展历程

19 世纪 50 年代，为解释水在土体中的渗透规律，达西通过大量实验归纳出渗流线性定律，渗流理论由此系统化^[1]。此后，裘布依在达西定律基础上研究单向与平面径向稳定渗流，奠定了稳定渗流理论的基础；1935 年，齐斯发表非稳定渗流的研究成果，用于描述流动随时间变化的完整过程^[6]。

20 世纪 20 年代起，石油和天然气开采成为渗流力学的重要驱动力；30 至 40 年代，单相流体在均质地层中的渗流问题基本解决，相渗透率、毛细管压力等概念相继建立，并形成油水两相一维驱替的数学模型^[2]。20 世纪 60 年代以后，研究进一步扩展到非均质介质、多重介质、非牛顿流体、非等温渗流和物理化学渗流等方向，数值计算与数值模拟逐渐成为主要手段^[4]。

应用

地下水资源开发是渗流理论最早的应用领域。含水层中地下水的补给、排泄与井群干扰计算，均需借助渗流模型预测水位与流量变化^[2]。与水利工程相关的渗流研究则关注坝体、堤基和围堰中的渗流场，用于评估渗漏量并分析渗透变形与渗透破坏^[9]。

在油气开发中，渗流理论用于分析油、气、水在储层中的多相流动，是产能预测、井网部署和开发方案设计的基础。对低渗透储层，常采用水力压裂形成人工裂缝，并在模型中考虑启动压力梯度和高速非达西效应以描述压裂后的流动特征^[10]。

在人造多孔介质方面，过滤、填充床、色谱分离、膜分离等过程的流体运动常被当作渗流问题处理；生物体内的血液流动与组织间质流动也被纳入生物渗流的范畴^[11][2]。地下储气、地热开发、污水处理与城市地面沉降等问题的研究同样依赖渗流分析^[2]。

局限

达西定律的线性关系只在低速条件下成立。流速很高时惯性效应不可忽略，流动偏离线性规律，属于非达西流，需要用福希海默方程等引入惯性项进行修正^[12]。在低渗透介质中还存在启动压力梯度现象，即压力梯度低于某一阈值时流体几乎不流动，经典达西定律同样不能直接描述^[12][13]。

把多孔介质视为连续介质、以宏观参数概括孔隙结构，会忽略孔隙尺度下的流动细节与界面效应。达西定律还忽略边界影响，在多孔介质与开放流体交界处无法给出速度的连续过渡，需借助布林克曼方程引入有效黏性项处理^[14]。对于页岩等孔隙直径接近气体分子自由程的致密介质，单一宏观传输模型难以完整刻画气体在其中的运移^[15]。

参见

• 渗流力学 —— 研究流体在多孔介质中运动规律的学科分支。

• 达西定律 —— 描述低速渗流速度与水力梯度线性关系的基本定律。

• 多孔介质 —— 内部含连通空隙、可供流体通过的材料。

• 渗透率 —— 表征多孔介质本身传导流体能力的参数。

• 非达西流 —— 不服从达西线性关系的渗流现象。