达西方程与渗透率

Darcy方程与渗透率

Darcy方程

动量守恒律即是控制体元中运动流体的动量变化率等于所有有效的作用外

力总和，表示为：

du ?(1-2) ?F????

dt 在流体力学中，不可压缩粘性流体动量守恒律表现为著名的Navier-Stokes方程组。原则上讲，给出Navier-Stokes方程、连续性方程以及两个状态方程，就可以解决任意形状中不可压缩粘性流体运动的所有问题，而多孔介质中流体的流动也不例外。遗憾的是Navier-Stokes方程是非线性的，除非多孔介质中流体流动的孔道具有规则的几何形状和良好的对称性，否则不可能求解之。以多孔介质孔道的复杂性和不规则性，Navier-Stokes方程不可能直接适用于解决渗流力学问题，因此如果要想在多孔介质中应用流体力学基本方程组，必须根据流体渗流的特点对之加以改造。

在渗流力学中，动量守恒律表现为Darcy方程。这一方程原是Darcy在1855~1856年就法国Dijon的给水过滤问题进行的直立均质砂柱中水的流动实验所得到的。Darcy方程虽然是一个实验型定律，但是有具体的物理含义。根据流体力学中的Hagen-Poisseuille公式（Navier-Stokes方程的一个特解）采用毛细管模型、水力半径模型等可以导出Darcy方程的等价形式。而在低Renolds数的层流条件下，忽略渗流流体的惯性力影响并且对单元体取体积平均，由Navier-Stokes出发亦可以得出Darcy方程。以下我们回顾Darcy方程。 Darcy方程：

v??k???? (1-3)

其中，v是流体渗流速度，而流体势函数定义为：

??P??gz (1-4)

这里，ρ是流体密度，流体势?实质是每单位体积上的机械能。 平面均质流体渗流的Darcy方程： q?AK?P ?L(1-5)

Darcy方程是一种线性方程，它假定介质与流体之间没有耦合作用，而只表明流量与介质长度成正比；它只适用于一定的Reynolds数范围。Forchheimer等人就多孔介质中大Renolds数的非Darcy流动问题有过很好的研究。而Irmay等就Darcy方程的下限问题进行过讨论。有关论述可参阅Bear（1972）的著作。

与欧姆定律物理相似性（ρ=电导率）：

I??V?V?A R?L(1-6)

式中，I=电流（A），△V=电压降（V），R=电阻（Ω）

与Fourier定律物理相似性：

?T (1-7) q?AKhL式中，q=热量单位，Kh=导热系数（热量单位/（m*hour*K）），△T=温差（K）。

平均渗透率计算、毛管方程

（1）层状介质并联线性流、径向流： n?n K??Kjhjhj j?1j?1（2）层状介质串联线性流：

nK?L?LjKj

j?1（3）层状介质串联径向流： K?Lo?gre/rw?? ?nLogrjrj?1?Kjj?1（4）牛顿流体Hagen-Poisseuille圆管层流方程：

q??r4?Pr28?L，K?8

（5）牛顿流体Buckingham水平裂缝层流方程：

q?h3?Ph312?L，K?12 （6）牛顿流体Fanning圆管紊流方程：

v2?4r?Pf??L

（7）非牛顿幂率流体圆管层流方程（C=稠度指数，n=幂率指数）：

?1?P?1q??n??n?(1?3n)n?2CL????1?3n??r 1?n(1?n)(1?n)????n?e?C1?3n??rn?

K??n?(1?n)e?n2??1?3n??r （8）非牛顿Bingham流体圆管层流方程（C=稠度指数，n=幂率指数）q??r4?P?4?2?40/r?1?2r28?L??1??3????(?P/L)???3???0/r??(?P/L)????，K????8 （9）Klinkenbeig效应

KKgL?1?bP?K1g?mP (1)

(2)

(3) (4)

(5)

(6)

(7) (8) (9)

(10)

(11)

渗透率和孔隙度、有效应力关系

单一介质

（1）Kozeny-Carman方程（1927，1937，1997）

??3R23 K??，K?2K02?Sgv(1??)2?＝迂曲度,Sgv=式中，K0=无量纲常数，?=孔隙度，?＝孔喉半径，R=水力半径，比表面。

（2）渗透率和压差关系式：Henpumepob（1961），Pedrosa（1986），Kikani（1991）

1?K ??K?PK?K0e??(Pi?P)

式中，K0=初始渗透率，?=常数

（3）渗透率和正应力关系式：Louis (1974钻孔压水试验)

K?K0e???，???h

式中，K0=初始渗透率，?=常数，?=岩石容重，h=埋藏深度，?=正应力 （4）常晓林（1987）

假设渗透主轴与应力主轴保持重合，推广情形2 (i,j,k?1,2,3)

??(i??)???2?Ki?K0exp???i?2?P?? ??j??k?E1??????单一裂缝

（1）立方定律：平均N-V方程单一水平缝解

(2a)3Q?12??dP??? ?dl?式中，a=缝半宽，与Darcy定律相对比可知K~a3 （2）Snow (1968)

?KN?2b?2?Kh?K0????P?P0?

s????式中，Kh=水平渗透率，K0=初始压力P0下的渗透率，2b=裂缝张开度，s=裂隙间距。KN=裂隙法向刚度（ML?1T?2） （3）Jones (碳酸岩类岩石裂隙 1975)

??P??Kf?K0?log?h??

??P??3式中，P＝有效压力，Ph=有效压力（Kf?0） （4）Nelson（Navajo砂岩裂隙 1975）

Kf?A?BP?n

式中，P＝有效压力

（5）Gangi (基质裂隙平均模型1975)

2/3??????P?P??e K(P)?K0?1?C0?????

???P0?????4式中，C0＝充填影响常数，P0＝基质（颗粒）有效模量。

单一裂隙钉板模型：

??P?m?K(P)?k0?1????，0?m?1

P???1???3式中，P1=裂隙面粗糙度有效模量 （6）Kranz (Barre花岗岩裂隙 1978)

QKf?0(Pc?Pf)?n

A式中，A＝过水面积，Pc=总压力，Pf=孔隙流体压力 （7）Walsh & Grosenbaugh (1979)

?K????K0?1/3?2h??P??1???ln??

?a0??P0?e式中，h=几何平均缝宽，Pe=有效应力，Pe?Pc??PP （8）Zhang Bochong (1986)

K1/3?A?Bln(Pc??Pp) 式中，A、B为常数。