水文地质边界条件处理-陈崇希.ppt

ijkbpqXYO图4-6-1格点i的一个辅助三角形图４-6-２有限差分计算井水位及抽水井附近水头分布HiHwHjijrerw地面数值法水头线实际水头线井壁设有限差分法解出的Hi相当于某个径距re处的水头，如图4-6-2所示。当不稳定的径向流形成拟稳定态时，通过pb断面进入抽水井的水流通量为（4-6-2）其中（4-6-3）由式（4-6-1）和式（4-6-2）计算的流量相等，于是有（4-6-4）式（4-6-3）代入式（4-6-4）得（4-6-5）因此常规数值模拟求得的水位Hi对应的径距——可称为等效半径re（4-6-6）根据不规则多边形有限差分法，有（4-6-7）其中：而S?ijk为三角形ijk的面积。把式(4-6-7)代入式(4-6-6)，得到对应?pib的re。同理可求得其它各流入角(控制角)的re。对于实际问题,由于非轴对称的边界、非轴对称的不均质、多抽水井及混合抽水井、混合观测孔等条件,导致流场通常是非纯轴对称的，各流入角(控制角)的re并不相等。作为近似，我们以Qpb,Qqb……角流量为权重，来求取re的平均值，即（4-6-8）式中：?表示对i点周围所有辅助三角形求和。假定上述i~j(i~k，…)范围内已形成拟稳定态，于是可利用Theis近似式，即（4-6-9）（4-6-9’）（4-6-9）式就是多边形有限差分法的完整抽水井水位校正公式。式中是格点i抽水井水位修正值。当时,；当时,；只有当时,。可见，如果不校正的传统有限差分法,并非剖分愈细愈能刻画抽水井水位。如果井点周围剖分成6个辅助等边三角形，则?pib=30?，由式（4-6-6）可得（4-6-10）对于井水位无需校正的条件是即（4-6-11）即表明:当