大区域地下水模拟的理论和方法.pdfVIP

大区域地下水模拟的理论和方法 　　1. 引言 　　在我国许多地区，由于水资源严重短缺，超采地下水，导致大 范围的地下水位急剧下降，从而引起区域性的地面沉降、海水入侵以 及湿地退化等生态环境问题，严重影响了区域经济社会可持续发展。 如何合理地利用地下水资源，对区域和流域地下水、地面水进行优化 配置，通过建立健全的流域水循环体系，保护地下水环境，恢复流域 生态，实现水资源的可持续利用以及人与自然的协调，是当今水资源 研究领域十分重要的课题。 　　要合理地利用区域地下水资源，首要问题是要对区域地下水资 源进行符合实际的模拟，摸清区域地下水的运动规律。然而，由于区 域地下水模拟涉及大尺度地下水运动，而传统的地下水模拟手段和方 法在处理大尺度地下水运动模拟中所涉及的有限元计算中的时空尺 度的选择、水位和含水层信息不足等问题上方面存在困难，引起由于 信息不足带来的较大计算误差。因此，针对尺度和信息不足问题，建 立大区域地下水运动模拟的理论和方法十分必要。 　　本文在融合地质统计、逆问题理论和地下水运动理论的基础上， 提出了建立大区域地下水运动模拟的理论和方法。首先，提出了大区 域地下水有限元计算的时空尺度理论公式；其次，针对信息不足，提 出了大区域地下水位推定的方法；第三，根据推定的地下水位，运用 逆解析理论法对大区域地下水含水层透水系数进行逆推定；最后，给 出了运用实例。 　　2. 大区域地下水有限元计算的时空尺度 　　在大区域地下水数值计算中，首先需要选择研究对象所需划分 的计算网格大小或确定计算时间步长大小。尤其在模拟河道网基础上 的分布式水文模型和大区域地下水模型进行地表水与地下水偶合计 算时，同时兼顾河道网的数值地形网格大小和地下水网格大小，既满 足水文计算的精度同时使地下水计算不至于失稳就显得至关重要。然 而，由于大区域地下水计算中，由于范围较大，人们往往希望选择较 大的网格进行计算，以减少计算工作量。这往往导致计算上的不稳定。 　　对于地下水有限元计算的时空尺度（时空步长）的选择，1977 年，Newman 等对二维地下水运动方程的有限元解法中的不稳定问题 进行了分析[1]，推测不合理的时间步长导致了混合型差分的出现， 由此导致了解的不稳定问题。针对Neman 等的推测，Zhang （1992）[2] 和Wood （1996）[3]提出了二维地下水有限元计算的时间步长的条件。 2002 年，张祥伟等[4]运用最大最小原理对大尺度二维和准三维地下 水有限元计算的不稳定问题进行了理论分析，提出了二维和准三维地 下水有限元计算的时空步长的理论公式。即对于准三维地下水有限元 计算的时空步长分别为： 　　空间步长条件： 　　 　　3. 大区域地下水位的推定 　　在进行大区域地下水计算中，需要根据实测的地下水位推定初 始流场，以便检验地下水数值计算精度、进行非恒定地下水计算以及 识别含水层参数。 　　对于初始流场的推定的方法，通常有 Universal Kriging 方法 [5-6]，UK 法的计算行列比较大，计算比较复杂。Newman 等（1984） [7] 和 Sun （1999）[8] 运用 Residual Kriging （RK）法进行区域地 下水位的推定，也就是将实测地下水趋势面去除得到正态残差，将正 态残差运用 Ordinary Kriging 法进行面上残差的推定，再加上实测 地下水面的趋势得到三角形网格上各点的推定地下水位值。 　　上述方法在大区域地下水计算中遇到信息不足的问题，当实测 地下水位信息不足时，推定的初始流场会带来较大的误差，本文根据 地形水文学的原理[9]，即地下水与地形之间存在的相关性，提出运 用数值地形模型（DEM）中的地形标高作为辅助信息，修正实测地下 水位得到的地下水趋势面，然后，运用 Ordinary Kriging 方法对修 正后趋势面得到的正态残差对三角形网格点上的残差进行推定，每各 网格点上推定的残差再加上修正的趋势值得到面上地下水位的推定 值。本文将此方法定名为 ROKMT （Residual Ordinary Kriging with Modified Trend）方法[10-11]，即地下水位由以下两部分构成： 　　 　　大区域地下水位推定的方法如图 1。 　　 　　4. 区域地下水透水系数的识别 　　在地下水流动模拟中，一般根据含水层的水文地质条件和土壤 特性，选定透水系数的取值。在大区域地下水计算中，含水层物理特 性的信