习题及答案未

1、3-1基本微分方程中没有包含水的密度，为什么说它表示了质量守恒定律答：首先，连续性方程：_良PVx）十一乜）+可PVz）（诙制CZ.AxAyAz=【PnAxAyAz】ct表达了渗流区内任何一个局部”所必须满足的质量守恒定律,方程都是根据连续性方程和反映动量守恒定律的方程（如且各种研究地下水运动的微分Darcy定律）建立起来的。其次，在基本微分方程，如承压水运动微分方程中:xzFz它表明单位时间内流入、流出单位体积含水层的水量差等于同一时间内单位体积含水层弹性释放（或弹性贮存）的水量。它还通过应用Darcy定律反映了地下水运动中的能量守恒与转化关系。可见，基本微分方程表达了渗流区中任何一个“局部

2、”都必须满足质量守恒和能量守恒定律。这一结论也适用于半承压水运动和潜水运动的基本微分方程。最后，由:Rs=-dVb+dV=aPg+nPPg=Pg(a+nP),所以在各基本微分方程当中，水的密度的影响是通过贮水率出来表达的。3-2推导渗流的连续性方程、承压水运动的基本微分方程、半承压水运动的基本微分方程、潜水运动的基本微分方程（选做其二）答：（1）连续性方程:设在充满液体的渗流区内，以p（x,y,z）点为中心取一无限小的平行六面体（其各边长度分别为*,Ay,Zz、且和坐标轴平行）作为均衡单元体（图1)。Axc'b'p(x,y,z)d'a'图1均衡单元体如p(x,y

3、,z)点沿坐标轴方向的渗流速度分量为Vx、Vy、Vz,液体密度为P,那么,通过abcd面，在N时间内流入的水流质量Pvx可利用Taylor级数求得:111iVx(1.1)：vxy.z：t=：-vx-.：x.：y.：z.：t同理,可求出通过右侧a'bcd面流出的质量为PvxQyAzAt=(1.2)因此,沿x轴方向流入和流出单元体的质量差为:【麻一；*-JPVx十:等yAzQ=-二)xiyAzit2Fx|-2；:xI:x(1.3)同理，可以写出沿y轴方向和沿z轴方向流入和流出这个单元体的液体质量差，分别为：:-UVyVzIAxAyAzAt和-AxAyAzAt.:y：z因此，在t时间内，流入

4、与流出这个单元体的总质量差为：l3：Vxf：Vyil-CVz一|+'AxAyAzM(1.4)excycz在均衡单元体内，液体所占的体积为nAxAyAz,其中n为孔隙度。相应的，单元体内的液体质量为piAxAyAzo因此，在'时间内，单元体内液体质量的变化量为：PnAxAyAzkt(1.5)ft在连续流条件下(渗流区充满液体等)，根据质量守恒定律，两者应该相等。因此，-*+等)+*xWzPnayAz(1.6)式(1.6)即为渗流的连续性方程。(2)承压水运动的基本微分方程：假设，地下水流动主要是沿水平面方向进行，垂直流速可以忽略，只考虑垂向压缩。于是，只有水的密度c孔隙度n和单元

5、体高度z三个量随压力而变化，(1.6)式的右端可改写成；z1-n：Ft(2.1)kn'xy：z:n：；：z：-:t-ft二P二n：xyzft于是连续性方程（1.6）变为：(2.2)一Myzjn过汉旷zILxz；zft因为水头H=z,故有：Ft(2.(3),dV-.将dP=P=Ppdp式代入上式得；VFp.:H；:t-1-p；:t(2.(4)因为水的压缩性很小，l-3呼1所以,:p一;:Hg-4.:t(2.(5)将（2.5）式代入（2.2）式，得:cPcPVx一VyVz；x；y；zixAyAz(2.(6):H=：go-nxyz什上式中，左端第二个括弧项比第一个括弧项要小得多，因此可以忽略

6、不计，于2.6)式变为:cVx+Hy+方VzXxcycz；£Hx&yAz=Pg(a+nP)AxAyAz(2.(7)同时，根据Darcy定律在各向同性介质中，有:Vx-K-,Vy-K-,Vz.x二y(2.(8)将式（2.8）代入式（2.7),得:KIL:x.cH勾)iKH|Ixyz二用二n：Hxyz-z：z|；t(2.(9)根据贮水率的定义，上式可改写为:以卜色1K里卜色(K史KzsAxzsSr.ox，oyl勾1cz1以，ct(2.10)整理上式，得:(2.11)上述方程就是承压水非稳定运动的基本微分方程。(3)半承压水运动的基本微分方程：近似地认为水基本上是垂直地通过弱透水层

7、，折射90。后在主承压含水层中基本上是水平地流动的，主含水层中的水流可近似地作二维流问题来处理，水头看作是整个含水层厚度上水头的平均值，即：1MH=H(x,y,t)=/j0H(x,y,z”z(3.1)为简化起见，在以后叙述中略去H上方的横杠。同时假设和主含水层释放的水及相邻含水层的越流量相比，弱透水层本身释放的水量小到可以忽略不计。由图2所示的均衡单元体，根据水均衡原理可以写出下列形式的连续性方程Qx.也3.Qx也34(3.(2)(3.(3)|1>+2北2|lQy-y1-1Qy+yLt+W2-V1pxAy&t=*电AxAyAt11yq211y即2力就式中：V、V2分别为通过上部和

8、下部弱透水层的垂直越流速率或越流强度，即m1H2/V2-K2-；-K2：zH2-Hm2式中：m1和m2分别为厚度为M的承压含水层上、下的弱透水层厚度，K1和8分别为承压含水层上、下弱透水层的渗透系数。H1(x,y,t刘H2(x,y,t)分别为上含水层(图中为潜水含水层)和下含水层(图中为下承压含水层)中的水头，如以T表示主含水层的导水系数,则CTEH'(3.4)Qx=-T丁&ytxQy=-Txyy则有Qx-:Qx.xx,yv2x,yc0Qy二y9fx.x图2半承压含水层中的均衡单元体把式（3.4）代入式（3.2）,并在式的两端分别除以xdyAt,同时令Ax、Ay、At一0,.x

9、K1HmHH2-H：H+K2-=卜(3.5)m21t这就是不考虑弱透水层弹性释水条件下非均质各向同性越流含水层中非稳定运动的基本微分方程。（4）潜水运动的基本微分方程:潜水面是个自由面，相对压强p=0。对整个含水层来说，可以不考虑水的压缩性。先考虑一维问题。取平行于xoz平面的单位宽度进行研究。在渗流场内取一土体，它的上界面是潜水面，下界面为隔水底板，左右为二个相距X的垂直断面。上断面流入的流量一,:qxFqx为q，下断面流出的流量为q,设单位时间、单位面积上垂向补给含水层入2：:x2的水量为W（入渗补给或其它人工补给取正值，蒸发等取负值）。在At时间内，从上游流入和由下游流出的水量差，根据D

10、upuit假设为:容t/q故2）<qr三vxhrxt二xxt,X-X(4.1)在at时间内，垂直方向的补给量为waxAto因此，at时间内小土体中水量总的变化为:wx.：t一fx川,则在t.:t化。ftLXLt当潜水面上升时科为饱和差，下降时为给水度，此时忽略了水和固体骨架弹性贮存的变由于假设水是不可压缩的，根据连续性原理，这两个增量应相等，即:-W(4.(2)KdH将vx一dx代入上式，得；H=HxHCWWHHh1+一=&炽)KKct(4.(3)上式为有入渗补给的潜水含水层中地下水非稳定运动的基本方程（沿x方向的一维运动），通常称为Boussinesq方程。在二维运动情况下，可

11、用类似方法导出相应的方程为:2HKft(4.(4)小土体内水量的变化必然会引起潜水面的升降。设潜水面变化的速率为时间内，由于潜水面变化而引起的小土体内水体积的增量为:当隔水层水平时，上式中h=H。对于非均质含水层，Boussinesq方程有如下形式：(4.5)同Kh曳卜片Kh曳）+W"空exvcx/oyl勾/ctBoussinesq方程是研究潜水运动的基本微分方程。3-3为什么初始时刻可以任意选定，不一定选用地下水的原始状态？答：所谓初始条件，就是给定某一选定时刻（通常表示为t=0）渗流区内各点的水头值，用于表明所研究实际问题的特定条件。所以，利用渗流区内已知的水头H的分布情况，即可

12、以作为初始时刻。3-4为什么可以根据具体条件任意用一个区作为计算区，它的周界就作为边界？答：边界条件，即渗透区边界所处的条件，用来表示水头H（或渗流量q）在渗流区边界上所应满足的条件，也就是渗流区内水流与其周围环境相互制约的关系。边界条件和初始条件合称定解条件。一个描述某实际问题的数学模型中，定解条件用来表明所研究实际问题的特定条件。我们用这样的模型来再现一个实际水流系统。所以，只要其周界能够比较好的符合边界条件的选取要求，能够相对简单准确的给出数学模型的定解条件，该区域就可以作为计算区3-5如果选用天然边界作为计算区边界，有什么优越性？答：选取方便求解准确3-6边界上的泉一般作为什么边界条件

13、？如在开采过程中泉水可能被疏干，还能作为边界吗？答：通过实际观测可以得到边界上的泉单位面积（二维空间为单位宽度）上流出时的流量q时，因此泉可以作为第二类边界或给定流量的边界。泉水疏干后，相当于各点在每一时刻的水头都是已知的，此时这部分边界可作为第一类边界。3-7为什么一定要有识别（校正）模型这个阶段？直接用野外试验所得的参数值和边界条件建立模型，不经过上述阶段行不行？答：当已建立确定性模型来描述实际地下水流后，我们对通过上述步骤建立的模型是否能确实代表所研究的地质体还没有把握；模型中出现的参数这时一般也不能确切给出。因此，必须对所建立的模型进行检验，即把模型预测的结果与通过抽水试验或其它试验对

14、含水层施加某种影响后所得到的实际观测结果或一个地区地下水动态长期观测资料进行比较，看两者是否一致。若不一致，就要对模型进行校正，直至满意地拟合为止。这一步骤称为识别模型或校正模型。模拟实际问题的数学模型还应满足下列基本条件这个解对原始数据是连续依赖的（稳定性）。这意味着当参数或定解条件发生微小变化时，所引起的解的变化也是很微小的。只有有了这条保证，当参数和定解条件的数据有某些误差时，所求得的解才能仍然接近于真解；否则，解是不可信的，并应该认为此时的数学模型是有毛病的。在实际工作中，原始数据有某种误差，在所难免，所以这个条件很重要。所以一定要有识别（校正）模型这个阶段，不能直接用野外试验所得的参

15、数值和边界条件建立模型。3-8下图所示为均质、各向同性的潜含水层，地下水流为平面非稳定流，且与河水有直接的水力联系。已知渗流区内的A区为稻田区，其灌溉水的补给强度W（m/d）;B区为开采区，其开采强度为e（m/d）,在A、B区的右侧有3口分散的开采井，其开采量分别为Qi、Q2和Q3（m3/d）。区内其它地方均没有垂直方向的水量交换，试写出该地区地下水运动的数学模型解：边界BCD为天然隔水边界。地下水流与河水有直接的水力联系。因此，边界可以作第一类边界处理。在远离开采区的、实际上不受该区抽水影响的地段人为地划定一条边界。以二条流线(BA和DE)作为边界。这时计算区就由ABCDEA所围的区域组成。根据给出的条件描述这一均质潜水流的方程应是有入渗补给的潜水含水层中地下水非AE稳定运动的