第一章渗流理论基础

绪论地下水动力学：是研究地下水在孔隙岩石、裂隙岩石和岩溶岩石中运动规律的科学。它是模拟地下水流基本状态和地下水中溶质运移过程，对地下水从数量上和质量上进行定量评价和合理开发利用，以及兴利防害的理论基础。第一章渗流理论基础§1—1渗流的基本概念一、地下水在含水岩石中的运动1多孔介质：具有孔隙的岩石。含水介质一般分为三类：孔隙介质：含有孔隙水的岩层。裂隙介质：含裂隙水的岩层。岩溶（Karst）介质：含岩溶水的岩层。二、地下水和多孔介质的性质1地下水的状态方程丄dV丄dVVdp1dp

pdp等温条件下，水的压缩系数为：设初始压强p0时，水的体积为V。，当压强变到p时，体积变为V,由上式得：卜dV=—（JpdpVoVPo=e-（p-p）VooV=Ve-0（p-p）o用Taylor级数展开，舍去高次项，得到如下的状态方程：v=vo[i-B（p-po）]P=p[1-B（p-p0）]2多孔介质的某些性质（1）多孔介质的孔隙性孔隙度：指孔隙体积和多孔介质总体积之比。有效孔隙：互相连通的、不为结合水所占据的那一部分孔隙。有效孔隙度：指有效孔隙体积和多孔介质总体积之比。

死端孔隙：一端与其它孔隙连通，另一端是封闭的，其中的地下水是相对停滞的。（2）多孔介质的压缩性天然条件下，一定深度处的多孔介质，要受到上覆岩层荷重的压力。荷重增加，将引起多孔介质的压缩。多孔介质的压缩系数：dVdVbD8v=1vnV弄vs坐dV%Vd^1-nDVndVfl讷TV册dVb1dV1dVV1—ndVndV多孔介质的压缩包括固体颗—v的压缩和孔隙v勺压卿—vVwVs+vv上式令渗流渗流是一种假想水流。假想水流应有以下特点：假想水流的性质(如密度、粘滞性等)和真实地下水相同；假想水流充满含水层的整个空间；假想水流运动时，在任意岩石体积内所受的阻力等于真实水流所受的阻力bDVbDsVVd渗流渗流是一种假想水流。假想水流应有以下特点：假想水流的性质(如密度、粘滞性等)和真实地下水相同；假想水流充满含水层的整个空间；假想水流运动时，在任意岩石体积内所受的阻力等于真实水流所受的阻力上式变为：a=(1-n)a+nasp固体骨架的压缩性比孔隙的压缩性小的多，上式变为a=nap三、贮水率和贮水系数水位变化对含水层厚度的影响有效应力地下水位下降，水压力减小，有效应力增大，多孔介质被压缩。多孔介质的压缩包括固体颗粒的压缩和孔隙的压缩。但固体颗粒的压缩忽略不计。即(1-n)Vb=常数dK1-n)V]=dV-ndV-Vdn=0bbbbdVdnb=1—ndVd）b—Azd（AZ）dn———adc—adpAz1—nd（Az）=Azadpdn-G—n'）adp其对数等于0。即：此二式厚度变化和孔隙度变化与水的压强变化的关系2.贮水率和贮水系数贮水率：面积为1单位面积，厚度为1单位的含水层，当水头降低1单位时所能释出的水量。用p表示。s弹性释水：由于水头降低引起的含水层释水现象称为弹性释水。贮水系数：面积为1单位面积，厚度为含水层全厚度M的含水层柱体中，当水头改变一个单位时弹性释放或贮存的水量。用|J*表示。—*卄¥W二者关系：p*=psM四、渗流（4）通过任断面的流量及任一点的压力或水头均和实际水流相同。渗流区或渗流场：假想水流所占据的空间。2典型单元体五、渗流速度过水断面：垂直于渗流方向的一个岩石截面。渗流速度：通过单位面积的渗流量。v=Q/A渗流速度与地下水的实际平均流速有如下关系：V=nu六、地下水的水头和水头坡度u22g1地下水的水头式中：z—位置水头；P/Y——承压水头；二者之和为测压管水头。U2/2g——流速水头（很小忽略不计）。我们所说的水位就是测压管水头，这是基准面取的是海平面2等水头面和水力坡度等水头面：渗流场内水头值相同的各点连成的面。‘dH-J=-ndn等水头线：等水头面与某一平面的交线。水力坡度：大小等于梯度值，方向沿着等水头面的法线指向水头降低方向的矢量。rdHT°H「dHJ=-;J=-;J=一一xdxydyzdz梯度的大小为：矢量J在空间坐标系中的三个分量为：七、地下水运动特征的分类按地下水运动要素（渗流量、渗流速度、压强、水头）将地下水分为稳定流和非稳定流。稳定流：地下水运动要素不随时间变化。非稳定流：地下水运动要素随时间变化。研究地下水的运动规律实际上是研究地下水运动要素的时空变化规律。根据地下水的运动方向与空间坐标轴的关系分为一维运动、二维运动和三维运动。地下水的一维运动：地下水的渗透速度只沿一坐标轴的方向有分速度，其余坐标轴方向的分速度均为零。地下水的二维运动：地下水的渗透速度只沿二个坐标轴的方向有分速度，仅在一个坐标轴方向的分速度均为零。地下水的三维运动：地下水的渗透速度只沿空间三个坐标轴的分量均不等于零。八、地下水流态的判断地下水的运动有层流和紊流。

1.多孔介质判断法：用Reynolds数式中：v—地下水的渗流速度；d—含水层颗粒的平均粒径；Y—地下水的运动粘度。计算Reynolds数小于临界Reynolds数时，为层流；大于时，为紊流临界Reynolds数一般取150〜300。裂隙流判断法（指单个裂隙）主要是确定临界水力坡度。首先，由下式计算裂隙的相对粗糙度a=A/b式中：a—裂隙的相对粗糙度；△—裂隙的绝对粗糙度；b—裂隙的宽度。其次，根据下表确定临界水力坡度。最后，比较实际水力坡度和临界水力坡度，判断流态。§1—2渗流基本定律一、Darcy定律及其适用范围QV=dH~dSdH~dSv=-Kv=-K更xdxv=-K丝

y°yv=-K廻zdz或地下水的运动是三维，Darcy定律应该用微分形式表示：在直角坐标系中，如以v、v、v表示沿三个坐标轴方向的渗流速度分量，则有：xyz用矢量来表示渗流速度形式如下：v=vi+vj+vkxyz式中：i，j，k为三个坐标轴上的单位矢量。Darcy定律适用范围：Darcy定律中，渗流速度v与水力坡度J呈线性关系。做如下实验，固定某种直径d的砂粒，改变水力坡度J的大小，可得到对应的渗流速度v,按照Darcy定律应呈线性关系，但实际上，当v增大到某一值时，开始偏离Darcy定律,这时，根据v、d、J可确定Reynolds数（Re=vd/Y），计算出的Re一般在1一10。因此，Darcy定律适用的范围是：用Re=vd/Y（Y运动粘度）计算得Re小于1一10时，地下水的运动才符合Darcy定律。说明：地下水的运动绝大多数服从Darcy定律。例如，对d=0.5mm的粗砂为例，地下水15度时的运动粘滞系数Y=0.1m2/d，取Re=1时，由式Re=vd/Y求得：v=ReH=骑01=200m/dd0.51000当渗透流速小于200m/d时，地下水运动为Darcy流。一般粗砂的渗透系数K=100m/d，实际的J一般小于1/500，这里取1/500，可求得v=0.2m/d，远远小于200m/d，服从Darcy定律。因此，当渗流速度由低到高时，可把多孔介质中的地下水运动状态分为三种情况：（1）当地下水低速度运动时，即Reyn01ds数小于1到10之间的某个值时，为粘滞力占优势的层流运动，适用Darcy定律。（2）随着流速的增大，当Reyn01ds数大致在l到100之间时，为一过渡带，由粘滞力占优势的层流运动转变为惯性力占优势的层流运动再转变为紊流运动。（3）高Reyn01ds数时为紊流运动。二、渗透系数、渗透率和导水系数渗透系数：水力坡度等于1时的渗透速度。影响渗透系数的因素：①岩石性质（粒度、成分、颗粒排列、充填状况、裂隙性质及其发育程度）；②液体的物理性质（容重、粘滞性等）。渗透系数K可用下式：K=空=樣

卩Y式中：p——液体密度；g重力加速度；M动力粘度；k——渗透率。量纲[L2],只与岩石的性质有关，与液体性质无关。单位cm2或Darcy。导水系数：水力坡度等于1时，通过整个含水层厚度上的单宽流量。用T表示。导水系数与渗透系数的关系：T=KM三、非线性运动方程Re小于1一10时，地下水流为线性流，用Darcy定律描述；Re大于1一10时，地下水流为非线性流，用下列定律描述：Forchheimer公式：J=av+bv2Chezy公式1v二KJ2c一般地下水流都为Darcy流。§1—3岩层透水特征分类和渗透系数张量一、岩层透水特征分类据岩层透水性随空间坐标的变化情况，将岩层分为均质的和非均质的两类。均质岩层：在渗流场中，所有点都具有相同的渗透系数。非均质岩层：在渗流场中，不同点具有不同的渗透系数。非均质岩层有两种类型：一类透水性是渐变的，另一类透水性是突变的。根据岩层透水性和渗流方向的关系，可将岩层分为各向同性和各向异性。各向同性：渗流场中某一点在各个渗透方向上具有相同的渗透系数，则介质是各向同性

的。各向异性：渗流场中某一点在各个渗透方向上具有不同的渗透系数，则介质是各向同性的。二、渗透系数张量在各向同性介质中，渗透系数和渗流方向无关，是一个标量。在各向异性介质中，渗透系数和渗流方向有关。水力坡度和渗流的方向一般是不一致的流网一节中讲到）。这时，渗透系数是一个张量。需要掌握的是，在各向异性介质中，有三个主渗透方向，渗透系数分别为£、K2、K3v=-K退z3dzv=-K更;vv=-K退z3dzx1dxy2dy（或K、K、K）。三个主方向上渗透流速为:xyz§1—4突变界面的水流折射和等效渗透系数一、越过透水性突变界面时的水流折射介质I的渗透系数为牟，介质II的渗透系数为K2。界面上某一点附近的渗流速度和水头在两介质中的值依次为V］、v2和H「H2,位于界面上的任一点都应满足如下条件:H1=H21n21nV1n=V2ndH-K£x1ddH-K£x1dxtgO22t-KdH2dxtgOK1=―1-tgOK22因为H1=H，故则得:此式为渗流折射定律。几点结论:（1）当k1=k2，贝昨1=e2,表示在均质岩层中不发生折射。⑵当K］Hk2，而且K］、k2均不等于o时，如e1=o，贝ye2=o，表明水流垂直通过界面时不发生折射。（3）当K］MK2，而且K］、K2均为有限值时，女口e1=90，则有e2=90，表明水流平行于界面时不发生折射。（4）当水流斜向通过界面时，介质的渗透系数k值愈大，e角也愈大，流线也愈靠近界面。二介质的k值相差愈大，e,和e2的差别也愈大，流线通过界面后的偏移程度也愈大。二、层状岩层的等效渗透系数12有两种情况:①平行于层面的渗透系数；vv垂直于层面的渗透系数。平行于层面的等效渗透系数Kp设每一分层的渗透系数K和厚度Mi,如图。对于ii每一分层水力坡度是相等的，即J=AH/lq=KMiiiAH每一层的单宽流量为：通过层状含水层总流量为：nnAHq二乙q二乙KMq=KMiiiAH每一层的单宽流量为：通过层状含水层总流量为：nnAHq二乙q二乙KM一iiili=1i=1AH工

ii=1AHAHnAHKM——=乙KM——pliili=1如果我们用一等效的均值含水层代替层状岩层，这时式中：M—含水层的总厚度；K—等效渗透系数。p由此得：》KMiiK=i=1p为Mii=1等效渗透系数为：2.垂直于层面的渗透系数该情况下，通过各层的流量相同。但水头降落和水力坡度不同。总的水头降落厶H等于各分层水头降落厶H.之和。i对于每一层=Kb竺;

iMiAH=Mqi—KbiAH=5=q工MKbbKi=1ii=1iAH=Kb所以：取等效渗透系数Kv，v二式相等得：那么单宽流量为：v£M~Ki=1i£nMii=1K因此，i=1i因此，此式为层状岩层垂直于层面的等效渗透系数。说明：⑴当某一层的K.较小时，M./K.较大，K变小；当K-0时，M./K.fg,Kiiiviiiv-0,也就是说，垂直于层面的等效渗透系数主要取决于渗透系数最小的分层。(2)平行层面的等效渗透系数总是大于垂直层面的等效渗透系数。§1—5流网、流函数1.流线和迹线流线：某一瞬时，渗流场中处处和渗流速度矢量相切的曲线。迹线：把某一质点在连续的时间过程内所占据的空间位置连成线。流线的方程dxdyvvxyMb=dx,ab=dy因为：△Mab与AMAB相似，所以：或者：vd-vd=0xyyx此式为流线方程。流函数设有二元函数屮(x,y),满足—v；=vdxydyx=型dx=型dx+型dxdydyd屮d屮d屮=dx+dy=vdy-vdx=0dxdyxy该函数的全微分为：得：积分得：屮=常数该式表明：同一流线，函数屮=为常数，不同的流线则有不同的函数值。函数屮叫流函数。量纲［L2T-1］。流函数满足的条件是

=-v;=vdxydyPPVx1的单位时间单位面积的水流质量为用Taylor级数展开：流函数有下列特性：（1）对一给定的流线，流函数是常数。不同的流线有不同的常数值。流函数决定于流线。（2）在平面运动中，两流线间的流量等于和这两条流线相应的两个流函数的差值。（3）在均质各向同性介质中，流函数满足Laplace方程；其他情况下均不满足Laplace方程。（4）在非稳定流中，流线不断地变化，只能给出某一瞬时的流线图。因此，只有对不可压缩的液体的稳定流动，流线才有实际意义。二、流网及其性质流网：在渗流场内，取一组流线和一组等势线组成的网格。流网的性质：（1）在各向同性介质中，流线与等势线处处垂直，故流网为正交网格。（2）在均质各向同性介质中，流网每一网格的边长比为常数。（3）当流网中各相邻流线的流函数差值相等，且每个网格的水头差值相等时，通过每个网格的流量相等。（4）当二个透水层不同的介质相邻时，在一个介质中为曲边正方形的流网，越过界面进入另一介质中，则变成曲边矩形。三、流网的应用1.水头从流网图上可以读出渗流区内任一点的水头H。水力坡度和渗流速度J=AH/As；v=KJ流量q=vAl定性确定水文地质条件河流与地下水的补、排关系；等水头线的疏密反映导水性的大小；流线绕流时，遇弱透水层流线汇集时，遇强透水层§1—6渗流的连续性方程在渗流区内以P点为中心取一无限小的平行六面体，其边长分别为△x、Ay、Az,并且和坐标轴平行，设P点沿坐标轴的渗透速度分量为V、V、V,液体密度为P，则P点xyz处，单位时间内通过垂直于坐标轴方向单位面积的水流质量分别为PV、PV、PV那么,xyz。通过abcd面中点x-亍x-亍y，ZPVX1=pv(x,y,z)+xQ(pv)(Ax'pv+1PVX1=pv(x,y,z)+xQ(pv)(Ax'pv+1Waxx2dxpv-1如?Axx2dxd(pv)x—

dx-xAxAyAtAxAyAt-AxAyAzAtpv+1如2Axx2dxAxAyAtd(pv)一AAxAyAzAtdy略去二阶导数以上的高次项，得△t时间内由abed面流入单元体的质量为:同理，通过a'b'c'd'面流出单元体的质量为：沿x轴方向流入和流出单元体的质量差为：一一ddpp)+dAyAzAtd(pvz^^AxAyAzAtdxdydz同理，可得到沿y轴和z轴方向流入和流出这个单元体的液体质量差，分别为:在4t时间内，流入与流出这个单元体的总质量差为：在均衡单元体中，孔隙体积为n△x△y△z,其内液体质量为pn△x△y△z，△t时间内，单元体内(液体质量的变化为：d(pvdxAyAzAt根据质量守恒定律，消去△t得gnA^AyA強dz上二式应相等，因此=d[pnAxAyAzktdt++4^AxAyAz二ItpnAxAyAz]dz|dt此式为渗流的连续性方程(研究地下水运动的基本方程)。§1—7承压水运动的基本微分方程假设条件：水流服从Darcy定律；K不随p=p(p)的变化而变化;ps和K也不受n变化的影响；

」如+々+乂］AxAyAz二乞［pnAxAyAz］dxdydz」dt(4)含水层侧向无压缩，即△x、△y为常量，只有垂直方向4z的压缩。在连续性方程的右端项中，有三个变量，随压力p的变化而变化。三个变量随时间的变化转化成压力随时间的变化。液体压缩后，质量不变。即密度pdVdp__p-V和体积V变化，二者乘积不变。d(pV)=pdV+Vdp=0得__1dV_V~p~由水的压缩系数：得：_dV_V_所以，dp=ppdp前面给出了含水层厚度厶z和孔隙度n随压力p的变化关系:d(Az)=△zadp；dn=(1-n)adp—[pnAxAyAz]=d—[pnAxAyAz]=dtAxAynp+pAz+nAz——dp+pAz(1_n)dp+pAz(1_n)adp+nAzpBdpdtdtdt_p(a+ndt式中：a为多孔介质压缩系数。将三式代入连续方(程右)端项得：「d(pv)dpv)d(pv)_I++Idxdy=npAzaAxAydAxAyAz=p(a+nB)-PAxAyAzdtdpdtdHdt于是连续性方程变为将化为：H=4+Hg空-zgdp_pg更+(H_z)g空

dtydtdtdtdtdt因为故有：p=Y(H-z)=pg(H-z)

dHdtdp_pgdHdt1一卩pdt或：将dp=pBdp代入，得：即，aBztp.xdx什npaBztp.xdx什npzdz_p2g(a+ng)dHAxAyAzdt因为水的压缩性很小，Bp忽略不计代入前式，得-pd(v)d()d(v)'++z—AxAyAz_p2g(a+ng)dHAxAyAzdt第二项p非常小，忽略不计，于是上式变为:根据Darcy定律：v_一K-xdxv_-kdH

zdz「ddH)「ddH)d(K+—dxkdx丿dyIKdy丿dzdHHkdz丿AxAyAz_pg(a+ng)dHAxAyAzdt一ddH]dJdH、d<dH)K+—+—Kdxkdx丿dykdy>丿dzkdz丿AxAyAz_卩AxAyAzsdt1.在各向同性介质中，有：代入上式，得因为ps=pg(a+nB)d=dH、dxKd=dH、dxKdx丿dfkdH]+dyIdy丿d<dH)K

dzIdz丿_卩sdHdt所以上式变为：两边消去单元体体积△x△yAz,得：此式为非均质各向同性介质承压水流微分方程。

v=-K廻xxxdv=-K廻xxxdxv=-K更;vyyydyz=-K迥zzdzdx(dH)d(KdHd(dH)K++—KVxxdx丿dyVyydy‘丿dzVzzdz丿dH7~dtd2Hd2Hd2H卩dH++=

dx2dy2dz2Kdt2.对于各向异性介质：非均质各向异性介质承压水流微分方程为：3.对于均质各各向同性介质，K为常数，承压水流微分方程为：单+dH电十甸傀dH^}dHdHddX.2ddj^dyZfdjKJdtsdt单+4.地下水流为二维流时，非均质各向同性介质承压水流微分方程为：5.柱坐标：如果能用柱坐标表示，则x=rcos0、y=rsin0，代入1d(dH)1d(dH)r

dr丿rdr1d2Hd2H++=

r2d02dz2pdHs--Kdt可化成式6.有源汇项，用W表示。源：在垂向上有水流入含水层称源。w为正。d(dd(dH]drdH\d<dH)K+—K+—KdxVdx丿dy1dyJdzVdz丿汇：在垂向上有水流出含水层称汇。W为负。有源汇项时，只需在上述方程中左边加W即可。dHdtdt如各向同性介质：7.稳定流：水位H不随时间变化，即,上述微分方程的右端项等于零。§1—8越流含水层中地下水非稳定运动的基本微分方程越流含水层(半承压含水层)：当承压含水层的上、下岩层(或一层)为弱透水层时，承压含水层可通过弱透水层与上、下含水层发生水力联系，该承压含水层为越流含水层。越流：当承压含水层与相邻含水层之间存在水头差时，地下水便会从高水头含水层通过弱透水层流向低水头含水层，这种现象称越流。假设条件：水流服从Darcy定律；K不随p=p(p)的变化而变化；p和K也不受n变化的影响；含水层侧向无压缩，即△X、△y为常量，只有垂直方向4z的压缩。当弱透水层的渗透系数£比主含水层的渗透系数K小很多时，近似认为水基本上是垂直地通过弱透水层，折射90°后在主含水层中基本上是水平流动的。(如£与K相差较

小时，用等效渗透系数，非越流）。（6）和主含水层释放的水及相邻含水层的越流量相比,弱透水层本身释放的水量小到可以忽略不计。越流含水层的微分方程：如图，越流含水层的厚度为M,水头为H；上、下各有一厚度为m]和m2、渗透系数为£和K2的弱透水层。上覆潜水含水层的水位为H],下伏承压含水层的水位为H2。在含水层全厚度上取一单元体。水平宽为△x,△y。单元体的中点为P，在P点沿x方向单位时间，通过面积4yM的流量为Qx,沿x轴流5Qxdx量的变化率为则dx2沿x轴流入单元体的水量为：沿x轴流出单元体的水量为门dQAxQ+x•—xdx2沿x轴单位时间流入流出单元体的水量差为：同理，可得沿y轴单位时间流入流出单元体的水量差为：在Z轴方向：由下部承压含水层单位时间流入越流含水层单元体的水量为：Qz2=v2△xy向上部含水层单位时间流出越流含水层单元体的水量为：Qz]=v]△x•△ydQxAdQxAx-dQydxdyAy+(v-v)AxAy=卩*21dH

dtAxAy单位时间内沿Z轴方向流入和流出单元体的水量为：（v2-V]）△x•△由质量守恒定律有：卩*dHA卩*dHAxAy

dt因为：代入上式，得ddx(dH)d(、dH)(H-H“H-H)T+T+K2-KiIdx丿dydy丿、2m]m丿dH=卩*dt消去△x△y,得：此式为非均质各向同性越流含水层中地下水非稳定运动的基本微分方程。对于非均质各向异性介质，有：

dxCdH)dLdH)T+—Tf+Vxxdx丿ddxCdH)dLdH)T+—Tf+Vxxdx丿dyVyydy丿Vd2Hd2HH-H++K—2-dx2dy22Tm2KU-K2H-H)11m/1dHdtH-H+K11Tm1卩*dH-~T~dt卄郁dHIyydy对于均质各向同性介质，有有源、汇项的情况下：非均质各向同性ddTT泗［&dx))dyddH--HH--H22:—JK11m11m22卄W=MLT))ddH**~dttd2Hd2H++Kdx2dy22非均质各向异性H-HTm2H-H1Tm1卩*dHT~dt均质各向同性H-H1m1上式中的K/m（用O，表示）为越流系数。越流系数：当主含水层和供给越流的含水层间水头差为个长度单位时，通过主含水层越流系数：当主含水层和供给越流的含水层间水头差为个长度单位时，通过主含水层B=尹K和弱透水层间单位面积界面上的水流量。越流因素B,量纲为［L］。§1—9研究潜水运动的基本微分方程一、Dupuit假设潜水面是弯曲的,等水头面也是弯曲的,潜水流的运动不是水平的。实际上潜水面的坡角很小。Dupuit假设：假设潜水面比较平缓,等水头面铅直,水流基本上水平,可忽略速度的垂直分量。同一剖面各点的渗透速度相等。Dupuit假设是忽略了渗流速度的垂直分量vz，但是，有些地段垂向分速度较大，不能采用Dupuit假设。也就是说,垂向分速度不能忽略。此外，还有降深较大的抽水井附近。二、Boussinesq方程在Dupuit假设和不考虑水的压缩性的条件下。考虑二维问题（含水层不水平），在渗流场内取一土体。如图。水平宽为厶x,△y。

单元体的中点为P,在P点沿x方向单位时间，通过面积△yh的流量为Qx,沿x轴流x量的变化率为5Qxdxdx2dQx•竺dx2dQ一xAxdx则沿x轴流入单元体的水量为：沿x轴流出单元体的水量为：沿x轴单位时间流入流出单元体的水量差为：同理，可得沿y轴单位时间流入流出单元体的水量差为:dQ—善Aydy单位时间内，垂直方向的补给量为：W△xAy"dQdQ]一——xAx一——Ay+WAxAyAtdxdydHdt一A一Ay+WAxAyAt—卩AxAyAtdtAt时间流入流出单元体的水量差为：土体内的水量变化引起潜水面的升降。假设潜水面的变化速率为则At时间内，土体内水的增量为：x据质量守恒原理将：廿…dx据质量守恒原理将：廿…dHYdxIdx丿dHdHQ=-KhAy;Q=-KhAxdxydy两个增量应相等。即dH、Kh^—AxAyAt+二KhAxAyAt+WAxAyAt—ydy丿dyAxAyAtdt代入上式，得d\Khd\KhdHdxIdx)d(vldH)+Kh丿dyIdy丿dHdt消去△xAy△t,得：此式为非均质各向同性潜水二维运动的微分方程。从微分方程中可知,潜水运动的微分方程是非线性的。对于三维流,这时考虑垂向分速度,其微分方程同承压水流微分方程。§1—10定解条件前面我们给出了,不同类型地下水的微分方程。除微分方程外,我们还应给出下列条件：方程中有关参数的值。主要有：渗透系数、导水系数、给水度、贮水率、贮水系数、越流系数,还有源、汇项等。