第3章习题及答案(3-5)未.pdf
3-1 基本微分方程中没有包含水的密度，为什么说它表示了质量守恒定律?答：首先，连续性方程：zyxntzyxzvyvxvzyx 表达了渗流区内任何一个“局部”所必须满足的质量守恒定律，且各种研究地下水运动的微分方程都是根据连续性方程和反映动量守恒定律的方程（如 Darcy 定律）建立起来的。其次，在基本微分方程，如承压水运动微分方程中：tHzHKzyHKyxHKxs 它表明单位时间内流入、流出单位体积含水层的水量差等于同一时间内单位体积含水层弹性释放（或弹性贮存）的水量。它还通过应用 Darcy 定律反映了地下水运动中的能量守恒与转化关系。可见，基本微分方程表达了渗流区中任何一个“局部”都必须满足质量守恒和能量守恒定律。这一结论也适用于半承压水运动和潜水运动的基本微分方程。最后，由：nggngdVdVbs，所以在各基本微分方程当中，水的密度的影响是通过贮水率 s来表达的。3-2 推导渗流的连续性方程、承压水运动的基本微分方程、半承压水运动的基本微分方程、潜水运动的基本微分方程（选做其二）。答：（1）连续性方程：设在充满液体的渗流区内，以 p（x，y，z）点为中心取一无限小的平行六面体（其各边长度分别为 x，y，z，且和坐标轴平行）作为均衡单元体（图 1）。xyzaaddccbbzyxp(x,y,z)o 图 1 均衡单元体 如 p（x，y，z）点沿坐标轴方向的渗流速度分量为 vx、vy、vz，液体密度为，那么，通过 abcd 面，在 t 时间内流入的水流质量1xv可利用 Taylor 级数求得：tzyxxvvtzyvxxx211 （1.1）同理，可求出通过右侧 abcd面流出的质量为 tzyxxvvtzyvxxx212 （1.2）因此，沿 x 轴方向流入和流出单元体的质量差为：tzyxxvtzyxxvvzyxxvvxxxxx2121 （1.3）同理，可以写出沿 y 轴方向和沿 z 轴方向流入和流出这个单元体的液体质量差，分别为：tzyxyvy和tzyxzvz 因此，在 t 时间内，流入与流出这个单元体的总质量差为：tzyxzvyvxvzyx （1.4）在均衡单元体内，液体所占的体积为 n x y z，其中 n 为孔隙度。相应的，单元体内的液体质量为 n x y z。因此，在 t 时间内，单元体内液体质量的变化量为：tzyxnt （1.5）在连续流条件下（渗流区充满液体等），根据质量守恒定律，两者应该相等。因此，zyxntzyxzvyvxvzyx （1.6）式（1.6）即为渗流的连续性方程。（2）承压水运动的基本微分方程：假设，地下水流动主要是沿水平面方向进行，垂直流速可以忽略，只考虑垂向压缩。于是，只有水的密度、孔隙度 n 和单元体高度 z 三个量随压力而变化，（1.6）式的右端可改写成；zyxtpnyxtpzntpnztpznzyxnt1 （2.1）于是连续性方程（1.6）变为：zyxtpnzyxzvyvxvzyx （2.2）因为水头pzH，故有：tptHgtp （2.3）将dpVdVd式代入上式得；tHpgtp1 （2.4）因为水的压缩性很小，l-p 1，所以，tHgtp （2.5）将（2.5）式代入（2.2）式，得：zyxtHngzyxzvyvxvzvyvxvzyxzyx2 （2.6）上式中，左端第二个括弧项比第一个括弧项要小得多，因此可以忽略不计，于是（2.6）式变为：zyxtHngzyxzvyvxvzyx （2.7）同时，根据 Darcy 定律在各向同性介质中，有：zHKvyHKvxHKvzyx,（2.8）将式（2.8）代入式（2.7），得：zyxtHngzyxzHKzyHKyxHKx（2.9）根据贮水率的定义，上式可改写为：zyxtHzyxzHKzyHKyxHKxs （2.10）整理上式，得：tHzHKzyHKyxHKxs （2.11）上述方程就是承压水非稳定运动的基本微分方程。（3）半承压水运动的基本微分方程：近似地认为水基本上是垂直地通过弱透水层，折射 90后在主承压含水层中基本上是水平地流动的，主含水层中的水流可近似地作二维流问题来处理，水头看作是整个含水层厚度上水头的平均值，即：MdztzyxHMtyxHH0,1,（3.1）为简化起见，在以后叙述中略去 H 上方的横杠。同时假设和主含水层释放的水及相邻含水层的越流量相比，弱透水层本身释放的水量小到可以忽略不计。由图 2 所示的均衡单元体，根据水均衡原理可以写出下列形式的连续性方程 tyxtHtyxvvtyyQQyyQQtxxQQxxQQyyyyxxxx*122222（3.2）式中：v1、v2分别为通过上部和下部弱透水层的垂直越流速率或越流强度，即 222222111111mHHKzHKvmHHKzHKv （3.3）式中：m1和 m2分别为厚度为 M 的承压含水层上、下的弱透水层厚度，K1和 K2分别为承压含水层上、下弱透水层的渗透系数。tyxH,1和tyxH,2分别为上含水层（图中为潜水含水层）和下含水层（图中为下承压含水层）中的水头，如以 T 表示主含水层的导水系数，则 xyHTQyxHTQyx （3.4）图 2 半承压含水层中的均衡单元体 把式（3.4）代入式（3.2），并在式的两端分别除以tyx，同时令x、y、t0，则有 tHmHHKmHHKyHTyxHTx*222111 （3.5）这就是不考虑弱透水层弹性释水条件下非均质各向同性越流含水层中非稳定运动的基本微分方程。（4）潜水运动的基本微分方程：潜水面是个自由面，相对压强 p0。对整个含水层来说，可以不考虑水的压缩性。先考虑一维问题。取平行于 xoz 平面的单位宽度进行研究。在渗流场内取一土体，它的上界面是潜水面，下界面为隔水底板，左右为二个相距 x 的垂直断面。上断面流入的流量为2xxqq，下断面流出的流量为2xxqq，设单位时间、单位面积上垂向补给含水层的水量为 W（入渗补给或其它人工补给取正值，蒸发等取负值）。在 t 时间内，从上游流入和由下游流出的水量差，根据 Dupuit 假设为：txxhvtxxqtxxqqtxxqqx22 （4.1）在 t 时间内，垂直方向的补给量为 W x t。因此，t 时间内小土体中水量总的变化为：txWxhvx 2xxQQxx 2xxQQxx 2yyQQyy 2yyQQyy yxv,1 yxv,2 M x y 小土体内水量的变化必然会引起潜水面的升降。设潜水面变化的速率为tH，则在 t时间内，由于潜水面变化而引起的小土体内水体积的增量为：txtH 当潜水面上升时 为饱和差，下降时为给水度，此时忽略了水和固体骨架弹性贮存的变化。由于假设水是不可压缩的，根据连续性原理，这两个增量应相等，即：txtHtxWxhvx （4.2）将 xHHdxdHKvx代入上式，得；tHKKWxHhx （4.3）上式为有入渗补给的潜水含水层中地下水非稳定运动的基本方程（沿 x 方向的一维运动），通常称为 Boussinesq 方程。在二维运动情况下，可用类似方法导出相应的方程为：tHKKWyHhyxHhx （4.4）当隔水层水平时，上式中 hH。对于非均质含水层，Boussinesq 方程有如下形式：tHWyHKhyxHKhx （4.5）Boussinesq 方程是研究潜水运动的基本微分方程。3-3 为什么初始时刻可以任意选定，不一定选用地下水的原始状态？答：所谓初始条件，就是给定某一选定时刻（通常表示为 t0）渗流区内各点的水头值，用于表明所研究实际问题的特定条件。所以，利用渗流区内已知的水头 H 的分布情况，即可以作为初始时刻。3-4 为什么可以根据具体条件任意用一个区作为计算区，它的周界就作为边界？答：边界条件，即渗透区边界所处的条件，用来表示水头 H（或渗流量 q）在渗流区边界上所应满足的条件，也就是渗流区内水流与其周围环境相互制约的关系。边界条件和初始条件合称定解条件。一个描述某实际问题的数学模型中，定解条件用来表明所研究实际问题的特定条件。我们用这样的模型来再现一个实际水流系统。所以，只要其周界能够比较好的符合边界条件的选取要求，能够相对简单准确的给出数学模型的定解条件，该区域就可以作为计算区 3-5 如果选用天然边界作为计算区边界，有什么优越性？答：选取方便 求解准确 3-6 边界上的泉一般作为什么边界条件？如在开采过程中泉水可能被疏干，还能作为边界吗？答：通过实际观测可以得到边界上的泉单位面积（二维空间为单位宽度）上流出时的流量 q时，因此泉可以作为第二类边界或给定流量的边界。泉水疏干后，相当于各点在每一时刻的水头都是已知的，此时这部分边界可作为第一类边界。3-7 为什么一定要有识别（校正）模型这个阶段？直接用野外试验所得的参数值和边界条件建立模型，不经过上述阶段行不行？答：当已建立确定性模型来描述实际地下水流后，我们对通过上述步骤建立的模型是否能确实代表所研究的地质体还没有把握；模型中出现的参数这时一般也不能确切给出。因此，必须对所建立的模型进行检验，即把模型预测的结果与通过抽水试验或其它试验对含水层施加某种影响后所得到的实际观测结果或一个地区地下水动态长期观测资料进行比较，看两者是否一致。若不一致，就要对模型进行校正，直至满意地拟合为止。这一步骤称为识别模型或校正模型。模拟实际问题的数学模型还应满足下列基本条件这个解对原始数据是连续依赖的（稳定性）。这意味着当参数或定解条件发生微小变化时，所引起的解的变化也是很微小的。只有有了这条保证，当参数和定解条件的数据有某些误差时，所求得的解才能仍然接近于真解；否则，解是不可信的，并应该认为此时的数学模型是有毛病的。在实际工作中，原始数据有某种误差，在所难免，所以这个条件很重要。所以一定要有识别（校正）模型这个阶段，不能直接用野外试验所得的参数值和边界条件建立模型。3-8 下图所示为均质、各向同性的潜含水层，地下水流为平面非稳定流，且与河水有直接的水力联系。已知渗流区内的 A 区为稻田区，其灌溉水的补给强度 W（m/d）；B 区为开采区，其开采强度为（m/d），在 A、B 区的右侧有 3 口分散的开采井，其开采量分别为 Q1、Q2和 Q3（m3/d）。区内其它地方均没有垂直方向的水量交换，试写出该地区地下水运动的数学模型。流线 A B B A E D 河 流 C 1 2 3 q 不透水边界 习题3-8图 解：边界 BCD 为天然隔水边界。地下水流与河水有直接的水力联系。因此，边界 AE可以作第一类边界处理。在远离开采区的、实际上不受该区抽水影响的地段人为地划定一条边界。以二条流线（BA 和 DE）作为边界。这时计算区就由 ABCDEA 所围的区域组成。根据给出的条件描述这一均质潜水流的方程应是有入渗补给的潜水含水层中地下水非稳定运动的基本方程，tHKKWyHhyxHhx。渗流区是 ABCDEA，记为D。边界 BA 和 DE 相当于隔水边界。数学模型如下：上和，；在上在）（DEBCDABnHAEtyxftyxHDyxyxHyxHtDyxtHyHzHyxHzHx0;,;,0,0,;KKQQQ-KKW0321式中：H0，f 为已知函数，f 为不同时刻的河水位；z（x，y）为隔水层标高。古今名言 敏而好学，不耻下问孔子 业精于勤，荒于嬉；行成于思，毁于随韩愈 兴于诗，立于礼，成于乐孔子 己所不欲，勿施于人孔子 读书破万卷，下笔如有神杜甫 读书有三到，谓心到，眼到，口到朱熹 立身以立学为先，立学以读书为本欧阳修 读万卷书，行万里路刘彝 黑发不知勤学早，白首方悔读书迟颜真卿 书卷多情似故人，晨昏忧乐每相亲于谦 书犹药也，善读之可以医愚刘向 莫等闲，白了少年头，空悲切岳飞 发奋识遍天下字，立志读尽人间书苏轼 鸟欲高飞先振翅，人求上进先读书李苦禅 立志宜思真品格，读书须尽苦功夫阮元 非淡泊无以明志，非宁静无以致远诸葛亮 熟读唐诗三百首，不会作诗也会吟孙洙唐诗三百首序 书到用时方恨少，事非经过不知难陆游 问渠那得清如许，为有源头活水来朱熹 旧书不厌百回读，熟读精思子自知苏轼 书痴者文必工，艺痴者技必良蒲松龄 声明 访问者可将本资料提供的内容用于个人学习、研究或欣赏，以及其他非商业性或非盈利性用途，但同时应遵守著作权法及其他相关法律的规定，不得侵犯本文档及相关权利人的合法权利。谢谢合作！
