习题与答案未

1、xyznxyzt3-1根本微分方程中没有包含水的密度,为什么说它表示了质量守恒定律 答：首先,连续性方程：vvxyvz ： ,-CCCxyz表达了渗流区内任何一个“局部所必须满足的质量守恒定律,且各种研究地下水运动的微分方程都是根据连续性方程和反映动量守恒定律的方程如Darcy定律建立起来的.其次,在根本微分方程,如承压水运动微分方程中：xK c c K 8 ) z K z ' s它说明单位时间内流入、流由单位体积含水层的水量差等于同一时间内单位体积含水层弹性释放或弹性贮存的水量.它还通过应用Darcy定律反映了地下水运动中的能量守恒与转化关系.可见,根本微分方程表达了渗流区中任何一个

2、“局部都必须满足质量守恒和能量守恒定律.这一结论也适用于半承压水运动和潜水运动的根本微分方程.最后,由：N=-+ =aP + pP = Pa+ BdVbdVgnggns,所以在各根本微分方程当中,水的密度的影响是通过贮水率体来表达的3-2推导渗流的连续性方程、承压水运动的根本微分方程、半承压水运动的基本微分方程、潜水运动的根本微分方程选做其二.答：1连续性方程：设在充满液体的渗流区内,以p x, y, z点为中央取一无限小的平行六面体其各边长度分别为 x, Ay, Nz、耳和坐标轴平行作为均衡单元体图 1.zIA x c.Ay ,|J 1,J.bb' px,y,z dd I 1aa y

3、ox图1均衡单元体x, y, z点沿坐标轴方向的渗流速度分量为vx、vy、vz,液体密度为通过abcd面,在 t时间内流入的水流质量pv可利用Taylor级数求得:x1v yztx1(1.1 )同理,可求由通过右侧a 'b'c d '面流生的质量为P AAA = |P vyztvxyzt2 (1.2 )x21vc(p因此,nc j t沿x轴方向流入和流生单元体的质量差为:1v1 x同理,vxyzvx22x&P )£(P )匕工- -x ccvx xyztx(1.3)可以写生沿yvzxyzt 和 xyztz因此,在 t时间内,流入与流由这个单元体的总质量

4、差为:v 八/ 工Cxyztvxyvz(1.4)的液体质量为pTIxyzA y A zo因止匕,Nj&y4 z匚向pn如孑醺淳相应的,单元体内 时间内,单元体内液体质量的变化量为：nxyzt(1.5)在连续流条件下渗流区充满液体等,根据质量守恒定律,两者应该相等.因此,vp AAA vxyvzxyz式1.6即为渗流的连毁修建xyznxyz + P(-c(1.6)0 )Z-A A A d2承压水运动的根本微分方程:假设,地下水流动主要是沿水平面方向进行,垂直流速可以忽略,只考虑垂向压缩.于是,只有水的密度P、孔隙度n和单元体高度z三个量随压力而变化,1.6式的右端可改写成;nxyznz

5、t(2.1 )于是连续性方程1.6变为:vxyvzxyzn=P(ct +cxyzpxyz t(2.2)p由于水头"故有：3rcP_ pdVpP(2.3)将dpd式代入上式得;V-P由于水的压缩性彳回、,l- Bn1所以,将2.5式代AR；£行口耳l Cvx=P (a +xcvyyvz Az式变为:同时,根据将式(2.8 )vxyvxyzxyz1pcHt cP(2.4)(2.5)xyz(2.6)一个括弧项要小得多,因此可以忽略不计,于是P f AAAAAA = P (ct +xyzgncDarcy定律在各向同性介质中,有:HHH vxK,vyK,v zK (2.8 )代入式2

6、.7 根据贮水率的定义,c='z c J上式可改写为:2.6 )(2.7)xyz c=P (a + P(2.9)xyz (2.10 ) t整理上式,得:H=H s -.11 )t上述方程就是承压水非稳定运动的根本微分方程.3半承压水运动的根本微分方程:近似地认为水根本上是垂直地通过弱透水层,折射90后在主承压含水层中根本上是水平地流动的,主含水层中的水流可近似地作二维流问题来处理,水头看作是整个含水层厚: )(3.1 )度上水头的平均值,即:HH=1)='mx,y,tHx,y,z,tdzM0为简化起见,在以后表达中略去 H上方的横杠.同时假设和主含水层释放的水及相邻含水层的越流

7、量相比,弱透水层本身释放的水量小到可以忽略不计.百图2所示的均询单元体,根据水均衡原座口以写由以下形式的连续性方程11-! 一 i + a +式中:xrcQyQ x x 2 xcQ y2 yctA +(A A AA A Av 2 vxytv1、v2分别为通过上部和下部弱透水层的垂直越流速率或越流强度,即-K11 m1K 2 H式中：m和m2分别为厚度为 M的承压含水层上、下的弱透水层厚度,K1和6分别为承压含水层上、下弱透水层的渗透系数.H1x,y,t / Hzx,y,t分别为上含水层图中为潜水=- Ad含水层和下含水层图中为下承压含水层中的水头,如以T表示主含水层的导水系数,(3.2)(3.

8、3 )(3.4)Qx2QxQxxxx2 CxAAAA A图2半承压含水层中的均衡单元体d ( c I £ 5 I 一一 M c把式(3.4 )代入式 L.2代 f 在式的两端分别除以 xytt同时令x、y干P一aC C d c JC那么有入和由下游流生的水量差,根据Dupuit 7段)T HT HHHHHHxx y y KK1 (3.5 )2mmt12这就是不考虑弱透水层弹性释水条件下非均质各向同性越流含水层中非稳定运动的基 本微分方程.(4)潜水运动的根本微分方程：潜水面是个自由面,相对压强 p= 0.对整个含水层来说,可以不考虑水的压缩性.先考虑一维问题.取平行于 xoz平面的单

9、位宽度进行研究.在渗流场内取一土体,它的上界面是潜水面,下界篇点水底板左右基多目瓶土牝白勺善宜断上疗面潞人的流量ex I 6) c J CC为 qx qx刃q,下断面流生的流量为q,设单位时间、单位面积上垂向补给含水层t时间内,从上游流的水量为W(入渗补给或其它人工补给取正值,蒸发等取负值).在4vh, 、x x t (4.1 )x在At时间内,垂直方向的补给量为W xA t o因此, t时间内小土体中水量总的变化为:Vx hxWxt小土体内水量的变化必然会引起潜水面的升降.设潜水面变化的速率为时间内,由于潜水面变化而引起的小土体内水体积的增量为：H口£.冈t当潜水面上升时 以为饱和

10、差,下降时为给水度,此时忽略了水和固体骨架弹性贮存的变化.由于假设水是不可压缩的,根据连续性原理,这两个增量应相等,即:Vx hHx Wx Wxtx =丛 t(4.2)dHVxK dx 将代入上式,得； =()HHx(4.3)上式为有入渗补给的潜水含水层中地下水非稳定运动的根本方程(沿x方向的一维运动),通常称为Boussinesq方程.在二维运动情况下,可用类似方法导由相应的方程为：6 ( 6 f d 工 K 金31h J 7-1h 寸 /1-(4.4)x x y y K K t当隔水层水平时,上式中h= Ho对于非均质含水层,Boussinesq方程有如下形式：C C y C ( C 1C

11、-HHHcKhJc KhjW e(4.5 )Boussinesq方程是研究潜水运动的根本微分方程.3-3为什么初始时刻可以任意选定,不一定选用地下水的原始状态？答：所谓初始条件,就是给定某一选定时刻(通常表示为t = 0)渗流区内各点的水头值,用于说明所研究实际问题的特定条件.所以,利用渗流区内的水头H的分布情况,即可以作为初始时刻.3-4为什么可以根据具体条件任意用一个区作为计算区,它的周界就作为边界？答：边界条件,即渗透区边界所处的条件,用来表示水头H (或渗流量q)在渗流区边界上所应满足的条件,也就是渗流区内水流与其周围环境相互制约的关系.边界条件和初始条件合称定解条件.一个描述某实际问

12、题的数学模型中,定解条件用来说明所研究实际问题的特定条件.我们用这样的模型来再现一个实际水流系统.所以,只要其周界能够比拟好的符合边界条件的选取要求,能够相对简单准确的给由数学模型的定解条件,该区域就可以作为计算区3-5如果选用天然边界作为计算区边界,有什么优越性？答：选取方便求解准确3-6边界上的泉一般作为什么边界条件？如在开采过程中泉水可能被疏干,还能作为边界吗？答：通过实际观测可以得到边界上的泉单位面积二维空间为单位宽度上流由时的流量q时,因此泉可以作为第二类边界或给定流量的边界.泉水疏干后,相当于各点在每一时刻的水头都是的,此时这局部边界可作为第一类 边界.3-7为什么一定要有识别校正

13、模型这个阶段？直接用野外试验所得的参数值和边界条件建立模型,不经过上述阶段行不行？答：当已建立确定性模型来描述实际地下水流后,我们对通过上述步骤建立的模型是否能确 实代表所研究的地质体还没有把握；模型中由现的参数这时一般也不能确切给曲.因此,必 须对所建立的模型进行检验,即把模型预测的结果与通过抽水试验或其它试验对含水层施加 某种影响后所得到的实际观测结果或一个地区地下水动态长期观测资料进行比拟,看两者是 否一致.假设不一致,就要对模型进行校正,直至满意地拟合为止.这一步骤称为识别模型或 校正模型.模拟实际问题的数学模型还应满足以下根本条件这个解对原始数据是连续依赖的稳定 性.这意味着当参数或

14、定解条件发生微小变化时,所引起的解的变化也是很微小的.只有 有了这条保证,当参数和定解条件的数据有某些误差时,所求得的解才能仍然接近于真解； 否那么,解是不可信的,并应该认为此时的数学模型是有毛病的.在实际工作中,原始数据有 某种误差,在所难免,所以这个条件很重要.所以一定要有识别校正模型这个阶段,不 能直接用野外试验所得的参数值和边界条件建立模型.3-8以下图所示为均质、各向同性的潜含水层,地下水流为平面非稳定流,且与河水有直接的水力联系.渗流区内的A区为稻田区,其灌溉水的补给强度W m/d ; B区为开采区,其开采强度为 £ m/d,在A、B区的右侧 3/d o区内其它地有3 口分散的开采井,其开采量分别为 Q、Q和Q m方均没有垂直方向的水量交换,试写由该地区地下水运动的数学模型习题3-8图解：边界BCD为天然隔水边界.地下水流与河水有直接的水力联系.因此,边界 AE可以作第一类边界处理.在远离开采区的、实际上不受该区抽水影响的地段人为地划定一条边界.以二条流线(BA和DE)作为边界.这日t计算区就由 ABCDE惭围的区域组成.根据给由的条件描述这一均质潜水流的方程应是有入渗补给的潜水含水层中地下水非ABCDE A记为稳定运动的