边界条件 源汇 等特殊问题的模拟与处理陈崇希.pptVIP

§4.13 初始水头条件的模拟 初始水头的分布,按教科书、软件等表述为; 初始水头的分布是地下水不稳定流数值模拟不可缺少的条件。通常初始水头的分布采用下列两种方法: 1． “插值法” 某专著认为：初始条件可“利用各观测孔的实测水头用插值方法”而得(薛禹群等,2019,249页) 。的确，以往不少模型也是如此做的。然而，此法存在两个问题： ①当研究区观测孔较少(这是普遍情况,特别是对于三维流模型)，且含有混合井孔、存在入渗补给滞后性等时，难以用插值的方法获得每一模拟层的初始水头分布(不要说像黑河中游地区分8层,苏州地面沉降模型23层等，就是运城地下水管理模型分3个含水层夹两个弱透水层，其第3层仅有3个观测孔水位)； ②插值方法是不考虑水文地质参数分布特征的,而水文地质工作者都十分明确：等水头线???正确勾画必需考虑水文地质参数的分布。这也是采用插值方法获得初始等水头线的一个敝病。; 2. “稳定流法” 即选择水头变化速率不大的时期（如某年某月）作为初始时刻,输入源汇和边界条件,或某时段内的平均值，运行初始参数分区和参数估计值的模型求其稳态解,以此作为初始水头分布。如此做法的问题是：研究区几乎不存在稳定流状态，如果启动模型的头几个时阶的源汇、边界条件与初始时刻的相同,那么模拟的整个流场也是稳定的,这与实际情况一般不符。当前大多数地区由于不合理开采地下水使得地下水头处于趋势性的下降状态。我们曾对于一个5000km2的盆地做过模拟试验, 如果边界条件和其它源汇保持不变,仅仅按规划设计增大地下水的开采量,没有500年甚至1000年的时间是稳定不了的(见第7章[例11-3])。; 我们建议： 1. s0法(陈崇希,1975;1983; 陈崇希等,1985) 对于小尺度(水源地和矿区等)研究区,大多是做一次大型抽/放水试验作为模型识别的基本数据,此条件的初始条件最好采取s0法,即将以水头H为应变量的数学模型改为以抽水引起的降深s为应变量的数学模型,从而使初始条件 H(x,y,z,t=0)=H0(x,y,z) 改为 s(x,y,z,t=0)=0。 使得初始条件亳无困难地给定。(第7章给出实例); 2. P-H0法(陈崇希等,1991) 由于模型识别-求取参数必需要给出初始条件，而初始水头的分布又与依赖于区域参数分布，两者互相依赖，谁也不能独立给定。基于此,我们提出“参数–初始水头迭代法”（P-H0法）(陈崇希等,1991)，即在求取参数过程中随时调整初始水头分布，“调参”与“初始水头”互为依存地轮流进行,是一个大迭代的过程。实际上，初始水头确定了，含水系统参数的分布也确定了。初始水头的分布是与模型识别-求取参数同时完成的。 对于实际模型，除了采用“s0法”外，建议都应采用P-H0法，以获得良好的效果。; §4.15 数值模型中流速(流量)的计算问题 通常有限元法和有限差分法以水头H为应变量，而流速和流量是由模拟获得的水头分布依Darcy定律计算而得。某专著认为：“目前有限元法…计算流量的方法…误差大，有时高达10%～20%，甚致更多” [2019，175页] 。果真有这么大的误差吗?是误差确实这么大，还是基本概念出差错?数值法不同于解析法,后者的水头解是连续函数,可在任意时空点取值，因而其水头一阶导数——水力梯度及流速可在任一时空点计算得；而数值法只能在预先设置的(时空)离散点上获得解，其它点上的值要通过插值而得。而利用所得解的水头再计算流速(流量)时,后者也是离散的,有其对应位置点,这个点位一般不与原设置的求水头的离散点位一致,依问题而定。 ; 以三角剖分的图4-15-1为例,为便于说明，假定含水层为均质、等厚的承压含水层，且地下水的流向垂直两个三角单元的共边 j k 线。显然,依三角单元 I j k 三个顶点的水头计算的水力坡度/流速小于由三角单元 j p k 三个顶点的水头计算的水力坡度/流速,从而下游的流速大于上游的,从物理背景上理解(例如下游有抽水井),这是再正常不过了。然而, 该专著明确：该流量是“通过单元某一边(线)断面，如 j k 断面的流量” [2019，175页],从而将两个单元各自三个结点的水头值确定的不同流速定位在两个单元的同一边线 j k 上,如此理解,当然会导致“计算流量