盆地多级次地下水流系统识别方法研究.docVIP

盆地多级次地下水流系统识别方法研究 区域地下水流理论较好地刻画了盆地尺度的地下水循环与演化规律 , 是指导 地下水生产、研究的有力工具。 识别盆地内部不同级次地下水流系统的空间分布 是应用该理论进行地下水资源勘查和解释一系列相关水文地质问题的前提和关 键。 本博士论文在深入分析多级次地下水流系统的自流状况、 矿化度空间分布和 滞留时间分布特征的基础上 , 提出了三种方法 , 有效识别盆地内部不同级次的地 下水流系统。前人关于自流状况的研究多为定性描述。 通过引入潜水面“相对起伏”因子建立地形起伏和潜水面分布的定量关系 , 并定义自流水头 , 推导得到了理想单元盆地和理想复杂盆地自流水头的解析解 , 定量分析了盆地自流状况及其影响因素。 研究表明 , 理想单元盆地中 , 自流水头在 河谷附近取值为正 , 在分水岭附近取值为负 , 垂向上近似呈负指数变化。 自流水头取值为零的等值线为自流区和非自流区的分界线 , 对应自流井发育 的临界深度。盆地自流状况受到盆地深宽比 , 潜水面“相对起伏”因子和渗透系 数的各向异性比和衰减系数的显著影响。 理想复杂盆地中 , 呈嵌套结构的多级次地下水流系统对盆地自流水头的空间 分布具有显著影响。盆地内部可以发育多个自流区和非自流区。 区域分水岭附近的局部河谷处 ( 或区域河谷附近的局部分水岭处 ), 自流水头 随埋深变化的曲线中存在极值 , 所在深度附近发育内部驻点。 在此基础上 , 盆地自 流状况的理论研究成果被应用于鄂尔多斯高原的两个野外实例。 理想单元盆地实例研究中 , 使用实测自流水头随埋深变化曲线的斜率估算了 渗透系数的各向异性比和衰减系数；理想复杂盆地实例研究中 , 使用实测自流水 头随埋深变化曲线中的极值识别了地下水流系统的多级次循环模式。 地下水矿化 度蕴含着地下水循环与演化过程中的重要信息。 首先通过引入地下水矿化度的控制方程 , 使用 COMSOL Multiphysics 模拟得 到了理想复杂盆地地下水矿化度的空间分布 , 发现存在两大突出特征：局部和区 域水流系统在盆地下游存在地下水矿化度的垂向突变界面； 两个相邻的局部水流 系统之间存在 “脊状” 局部高矿化度异常带。 这两大特征是后续地球物理勘探解 译的理论基础。 然后 , 尝试使用 EH-4 电导率成像系统获得了都思兔河流域中下游一条南北 向跨河剖面的视电阻率空间分布 , 并根据相关水文地质钻孔资料和勘探点附近的 民井水化学资料 , 验证了地球物理勘探结果的可靠性 , 建立了白垩系含水层视电 阻率与地下水矿化度的经验关系 , 进而近似得到了白垩系含水层地下水矿化度的 空间分布。最后 , 将地下水矿化度的实测分布与理论模拟分布相类比 , 综合考虑地 形和蒸发浓缩作用 , 解译了地球物理勘探剖面白垩系含水层中的多级次地下水流 系统。 盆地地下水滞留时间分布从时间域的角度审视了盆地尺度地下水的循环与 演化。依据不同级次地下水流系统滞留时间的离散特征 , 提出了适用于复杂三维 地下水流场的“晚峰”法。 “晚峰”法使用 MODFLOW和 MODPATH计算盆地地下水滞留时间分布密度函数 和盆地地下水滞留时间场 , 将盆地地下水滞留时间分布密度函数中的“晚峰”作 为临界滞留时间划分盆地地下水滞留时间场 , 进而识别局部、中间和区域地下水 流系统。为了验证“晚峰”法的适用性 , 分别引入二维、三维理想复杂盆地模型 作为算例。 算例结果表明 , “晚峰”法的识别结果与传统模型识别方法 , 驻点法 ( 二维 ) 和补给 - 排泄流线追踪法 ( 三维 ), 一致。最后 , 选择鄂尔多斯高原都思兔河流域开 展综合研究。 基于当地水文地质条件 , 使用 MODFLOW建立了该流域的地下水流数值模型 , 并使用实测地下水位和自流井分布信息进行了模型校正。在水流模型的基础上 , 首先使用 MODPATH定量研究了地下水的渗流路径和滞留时间 , 发现地下水径流长 度、穿透深度和滞留时间两两近似服从幂次定律；然后使用“晚峰”法识别了都 思兔河流域的局部和区域地下水流系统 , 定量研究发现局部水流系统具有径流长 度短 (l*lt;0.53), 穿透深度浅 (d*lt;0.2), 滞留时间短 (t*lt;0.18), 流域体 积占有率低 ( 约占 22%), 地下水循环量大 ( 约占 76%)的特点 ( 区域水流系统具有相 反特点 ) ；最后通过与基于水文地球化学和地球物理勘探的识别结果进行比较 , 证明了“晚峰”法的有效性。 本博士论文的相关研究成果不仅有助于深入认知鄂尔多斯高原地区地下水 的循环与演化规律 , 促进当地地下水资源的合理开发利用 , 还是对区域地下水流 理论的补充与完善。