地下水动力学7.pptx

根据Dupuit假设，可以建立潜水含水层中地下水运动的基本微分方程。考虑一维流情况。取垂直于xoz平面的单位宽度进行研究。在渗流场内取一个1-9研究潜水运动的基本微分方程2.Boussinesq(布西奈斯克)方程土柱，上界面为潜水面，下界面为隔水底板，左右为两个相距Δx的垂直断面。引起小柱体内水量变化的因素，除从左断面流入的流量和右断面流出的流量外，还有由大气降水入渗补给或潜水蒸发构成的垂向的水量交换。设单位时间、单位面积上垂向补给含水层的水量为W(入渗补给或人工补给为正，蒸发为负)。根据Dupuit假设，在Δt时间内，从上游流入和由下游流出的水量差为：

在Δt时间内，垂直方向的补给量为：。因此，Δt时间内小土柱中总的水量变化为：小土柱中水量的变化必然会引起潜水面的升降。潜水面变化的速率为则在Δt时间内，由潜水面变化而引起的小土柱内水体积的变化量为当潜水面上升时，μ为饱和差；下降时为给水度，此时忽略掉了水和固体骨架的弹性贮存变化。据连续性原理，两个变化量相等，即1-9研究潜水运动的基本微分方程2.Boussinesq(布西奈斯克)方程

根据Dupuit假设，带入：上式即为有入渗补给的潜水含水层中地下水非稳定运动的基本微分方程(x方向一维运动)。通常称为Boussinesq方程。在二维运动情况下为：当隔水层水平且位于z=0时，h=H（见图）对于非均质各向同性含水层，Boussinesq方程为：1-9研究潜水运动的基本微分方程2.Boussinesq(布西奈斯克)方程均同含水层

Boussinesq方程是研究潜水运动的基本微分方程。方程中的潜水含水层厚度h也是一个未知数，因此，它是一个二阶非线性偏微分方程。注意：推导方程时应用了Dupuit假设，忽略了弹性贮存；选取的单元是一个包含整个含水层厚度在内的小土柱，与推导承压水非稳定运动方程时所取的无限小的单元体不同。因此，应用Boussinesq方程得到的H(x,y,t)只代表整个含水层厚度上水头的平均值，不能用它来计算同一垂直剖面上不同点的水头变化。对某些无压渗流问题：如排水沟降低地下水位、土坝渗流等是不适用的？，此时应采用不用Dupuit假设的一般形式的方程：(Whynotμ？仅含水层中水量变化—弹性释水)1-9研究潜水运动的基本微分方程2.Boussinesq(布西奈斯克)方程（1）

对各项异性介质，如所选取坐标轴与各向异性方向一致，则有对于无压渗流来说，弹性释水与潜水面下降疏干所排出来的水量相比是微不足道的，故有时直接将上面两式的右端等于零，即和此时，相当于假设固体骨架是不可压缩的，同时忽略水的压缩性(ρ=常数)，即Ss=0，故右端等于零。1-9研究潜水运动的基本微分方程2.Boussinesq(布西奈斯克)方程（2）（3）（4）

以上四式与承压含水层中水头应满足的方程没有不同。方程右端采用的是贮水率Ss而不是给水度或饱和差μ。其原因是为了计算潜水面以下任一点的水头，此时，必须从渗流区内部取一个无限小的单元体来考虑。对这样的单元体来说，因为位于渗流区的内部，考虑其贮存量的变化时只能是弹性释水而不是疏干排水，因而推导出的无压水非稳定运动方程和承压水非稳定运动方程没有区别。在这种情况下，地下水非稳定流动的特征是由边界条件来反映的（如潜水面升降(降雨或抽水)所引起的水量变化等），此时重力给排水还可以作为一种补给水源。1-9研究潜水运