FLAC3D使用简单流体(simple-fluid)模式进行饱和/非饱和流动分析

此关系假设当基体吸力(孔隙压力是负时为正值)小于输入吸力极限时，饱和度为 1，而当基体吸力等于或大于输入吸力极限时，饱和度为0：

2.2 Brooks-Corey关系

Brooks-Corey (1964)的关系假设，当基质吸力小于进气吸力时，饱和度为1；而当基质吸力等于或大于进气吸力时，它是一个幂函数：

其中是以应力为单位的进气吸力，n是一个称为孔径指数的无量纲材料参数。

2.3 Van Genuchten关系

Van Genuchten (1980)的关系是：

其中是以应力单位为单位的参考吸力，m和n是无量纲材料参数。

2.4 Gardner 关系

Gardner (1958)的关系是：

其中是以应力单位为单位的参考吸力，a和n是无量纲材料参数。

2.5 Fredlund-Xing 关系

Fredlund-Xing（1994）的关系是：

其中C的值在 [0,1]；是参考吸力，是残余吸力，是极限吸力(典型值为1e6 kPa)，均以应力为单位; m和n是无量纲的材料参数；e是自然对数的底数。

2.6 表关系

表关系(table)允许通过分配表来输入表观饱和度和基质吸力之间的任何用户定义的分段关系(见第一个图)。该表的值由对来定义。表观饱和度是通过表中的线性插值方法计算的，该方法基于矩阵吸力的对数刻度(以10为底)。

2.7 FLAC3D命令

zone fluid saturation-suction命令设置表观饱和度类型(0~1)，作为基质吸力(负孔隙压力的相反符号值，单位为应力)的函数，仅在简单流体逻辑中有效。关键字如下：

(1) brooks-corey f1 f2 设置饱和-吸力的布鲁克斯-科里方程，并分别设置两个参数，即空气进入应力(以应力为单位)ψae 和孔径指数(无量纲)n。

(2) gardner f1 f2 f3 设置加德纳饱和-吸力方程，分别有三个参数：a、n和参考吸力 ψr。

(3) suction f 设置吸力阈值，当吸力小于该值时，视饱和度为 1，当基体吸力达到该值时，视饱和度为 0。

(4) table s 设置表观饱和度（table-y）的表名，作为矩阵吸力（table-x）的片断函数。请注意，插值是基于矩阵吸力的对数标度（以 10 为底）。

(5)  van-genuchten f1 f2 f3  分别用 n、m和参考吸力 ψr三个参数设置饱和-吸力的 Van Genuchten 方程。

(6) fredlund-xing f1 f2 f3 f3 f4 f5  设置饱和-吸力的 Fredlund-Xing 方程，分别有 n、m、参考吸力 ψr、残余吸力 ψs和极限吸力（如果以 kPa 为单位，典型值为 1e6）五个参数。

3.1 S 形关系

它是渗透率和表观饱和度的正常数值尺度中的S形函数。如果渗透率为对数刻度，则它不会显示为 S 形。

3.2 恒定关系

无论表观饱和度如何，渗透率都是常数1.0。由于比排水量不会随着表观饱和度的降低而降低，因此应谨慎使用。

3.3 表关系

与上述2.6节的描述相同。

3.4 FLAC3D命令

zone fluid permeability-saturation 

设置作为表观饱和度(0~1)函数的渗透率类型(0~1)，仅在简单流体(simple-fluid)逻辑中有效。关键字如下：

(1) constant 设置恒定的渗透率为 1.0，或表观饱和度不降低渗透率。

(2) s-shape 通过 3.1节所示的关系式，将渗透率k^ 设置为可观饱和度s^的s形函数。

(3) table s 设置渗透率(table-y)与表观饱和度(table-x)的表名称。请注意，线性插值是以渗透率的对数刻度(10 为底)为基础的。

block zone fluid property permeability [c_perm]

block zone fluid property porosity [c_porosity]

model configure fluid

zone fluid biot on

zone fluid cmodel assign isotropic

zone fluid prop perm 1e-10 poros 0.3

zone fluid property-distribution deviation-gaussian

zone fluid zone-based-pp