达西方程与渗透

Darcy方程与渗透率Darcy方程动量守恒律即是控制体元中运动流体的动量变化率等于所有有效的作用外力总和，表示为：duFdt(1-2)在流体力学中，不可压缩粘性流体动量守恒律表现为著名的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方程：(1-3)其中，v是流体渗流速度，而流体势函数定义为:(1-4)这里，P是流体密度，流体势实质是每单位体积上的机械能。平面均质流体渗流的Darcy方程：AK (1-5)Darcy方程是一种线性方程，它假定介质与流体之间没有耦合作用，而只表明流量与介质长度成正比；它只适用于一定的 Reynolds数范围。Forchheimer等人就多孔介质中大Renolds数的非Darcy流动问题有过很好的研究。而Irmay等就Darcy方程的下限问题进行过讨论。有关论述可参阅Bear(1972)的著作。与欧姆定律物理相似性(p二电导率)：IV（1-6）式中,I=电流（A）,△V=电压降（V）,R=电阻（Q）与Fourier定律物理相似性：1-qAOLL1-式中，q=热量单位，Kh=导热系数（热量单位/（m*hour*K））,△T=温差（K）(1)(1)(11)(11)(5)牛顿流体(5)牛顿流体h3Pq12L⑹牛顿流体vfL⑵⑶⑷⑸⑹⑺(8)(9)(10)平均渗透率计算、毛管方程层状介质并联线性流、径向流TOC\o"1-5"\h\zn nKKjhj hjji ji层状介质串联线性流：—/nKL/LjKj\o"CurrentDocument"/i1 '层状介质串联径向流：Logre/rwLog门几iKjji牛顿流体Hagen-Poisseuille圆管层流方程:r4P..r2q8L ，K8Buckingham水平裂缝层流方程K12Fanning圆管紊流方程：非牛顿幕率流体圆管层流方程(C=稠度指数，n=幕率指数)1Pn n (13q2CL13n1n(1n)(1n)eC—1n3nrn(1n)Ke21n3nrn非牛顿Bingham流体圆管层流方程(C=稠度指数，n=幕率指数)42q上±亠12_0/rk亡8L3(P/L)3(P/L)' 8Klinkenbeig效应KgKl Kg1b-P渗透率和孔隙度、有效应力关系单一介质(1)Kozeny-Carman方程(1927,1937,1997)R23，KKo 2Sgv(1 )2式中，K一无量纲常数，R二水力半径，-孔隙度，二孔喉半径，二迂曲度，Sgv-比表面。(2)渗透率和压差关系式：Henpumepob1961)，Pedrosa(1986)，Kikani(1991)丄上KT(RP)K°式中，Ko=初始渗透率， 二常数渗透率和正应力关系式：Louis(1974钻孔压水试验)KKoe，h式中，Ko=初始渗透率， =常数，二岩石容重，h=埋藏深度， =正应力常晓林(1987)假设渗透主轴与应力主轴保持重合，推广情形 2(i,j,k1,2,3)KiKoexp单一裂缝(1)立方定律：平均N-V单一裂缝(1)立方定律：平均N-V方程单一水平缝解3Q12dl(2a)3dP式中，3=缝半宽，与3Jones(碳酸岩类岩石裂隙1975)Darcy定律相对比可知PPKfKolog式中，P=有效压力，只=有效压力(Kf0)Nelson(Navajo砂岩裂隙1975)KfBP式中，P=有效压力Gangi(基质裂隙平均模型1975)PPePPPePoK(P) K。1C。式中，Co=充填影响常数,单一裂隙钉板模型:K(P)ko1式中，卩产裂隙面粗糙度有效模量Kranz(Barre花岗岩裂隙1978)Kf乂巳Pf)nA式中，A=过水面积，Fc=总压力，Pf=孔隙流体压力Walsh&Grosenbaugh(1979)1/3K J2h P1 InKo ao Poe式中，人二几何平均缝宽，Pe二有效应力，Pe巳 PpZhangBochong(198