渗流理论基础概要.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 用tg? 代替sin?的误差(略） Dupuit假设无效的地区： (1)存在入渗的潜水分水岭地段； (2)渗出面附近。渗出面(seepage surface)是在下游边界面上，潜水面以下、下游水面以上的地段。渗出面上潜水面往往和边界面相切，有较大的垂向分速度。 (3)垂直的隔水边界附近。 图1-17 Dupuit假设无效的地区 2 ）Boussinesq方程 根据Dupuit假设，可建立有关潜水含水层中的地下水流方 程。 图1-18 潜水的非稳定运动 （1）潜水一维流方程（沿x方向运动） 在Δt时间内，上、下游流入、流出单元体的水量差为： 在该段时间内，垂直方向的补给量为，故Δt时段水量总的 变化量为 由于水量的变化引起潜水面的升降，设其变化的速率为 则在Δt时段，由于潜水面的变化而引起的小土体内的水体积 的增量为 则有 将 （均质各向同性）代入上式，可以得到有 入渗补给的潜水含水层中地下水非稳定运动的一维流方程， 又称为Boussinesq方程： 式中K、?—— 潜水含水层的渗透系数、给水度； W——含水层单位时间、单位面积上的垂向补排 量，补给为正，排泄为负。 (2)潜水二维流方程 均质各向同性含水层，Boussinesq方程为： 式中h=H-Z，Z=0时，h=H。 （1-93） （1-94） 非均质含水层，Boussinesq方程为： 在推导潜水基本微分方程时应用了Dupuit假设，忽略了弹 性储存，所选的单元体是一个包括了整个含水层厚度在内的 土柱，这与承压水非稳定运动时选取的无限小的单元体不 同。所以，应用潜水运动基本方程得到的H(x，y，t) 只能代 表该点整个含水层厚度上平均水头的近似值，不能用来计算 同一垂直剖面上不同点的水头变化。 (3)潜水三维流方程 若不用Dupuit假设，Boussinesq方程的一般形式： （1-95） （1-96） 在上面的潜水基本运动微分方程中右端项为贮水率而不是 给水度，其原因在于，当不考虑Dupuit假设时，单元体位于 渗流区内部，其贮存量的变化只能是弹性释水而不是疏干排 水，因此推导出的潜水非稳定运动方程和承压水非稳定运动 方程形式一样。在这种情况下，地下水非稳定运动的特点由 边界条件来反映。 对于各向异性介质，坐标轴方向同主方向，有： 假设固体骨架是不可压缩的， ，同时假设忽略水的 压缩性,即?＝常数，有： （1-97） （1-98） 或 （4）潜水稳定运动的微分方程 没有入渗和蒸发时，潜水稳定运动的方程式为： 非均质 或 均质 （5）地下水运动基本微分方程的统一形式： （1-99） （1-100） （1-101） （1-102） 式中 Z——含水层底板标高。 1.7 数学模型的建立及求解 1 数学模型的有关概念 同一形式的偏微分方程代表了整个一大类的地下水流的运 动规律，而对于不同边界性质、不同边界形状的含水层，水 头的分布是不同的。而且对于偏微分方程而言，方程本身并 不包含反映特定渗流区条件的全部信息，方程可能存在无数 个解，如需要从大量的可能解中求得与特定区域条件相对应 的唯一特解，就必须提供反映特定区域特征的信息。 这些信息包括： （1）微分方程中的有关参数， 当这些参数确定后，微分方程才能被确定下来。 （2）渗流区范围和形状，当微分方程所对应的区域被确定之后才能对方程求解。 （3）边界条件(boundary conditions)：表示渗流区边界所处的条件，用以表示水头H（或渗流量q）在渗流区边界上所应满足的条件，也就是渗流区内水流与其周围环境相互制约的关系。 （4）初始条件(initial conditions)：表示渗流区的初始状态，某一选定的初始时刻(t=0)渗流区内水头H的分布情况。 将边界条件和初始条件并称为定解条件(definite solution condition)，微分方程和定解条件一起构成渗流场的数学模型。 数学模型：描述某一研究区地下水流运动的数学方程与其定解条件共同