水文地质学是研究地下水形成、分布、运动、水质演化、资源评价及地质环境效应的核心学科,而水文地质学模型是定量刻画地下水系统运行规律、解决实际水文地质问题的重要工具。简单来说,水文地质模型是对复杂天然地下水系统进行简化、抽象后的替代性表达,通过数学、物理、概念等形式还原地下水补给、径流、排泄全过程,实现地下水水量、水位、水质变化预测。按照表达形式可分为概念模型、物理模型、数学模型三大类,各自遵循统一的地下水动力学基本原理,在地下水资源评价、基坑降水、矿山疏干、地面沉降防控、污染羽迁移治理等领域广泛应用。

一、水文地质学核心基础原理

所有水文地质模型构建都建立在地下水运动基本定律与物质守恒原理之上,是模型运行的理论根基。

1. 达西定律 —— 地下水渗流基本定律

1856 年达西通过砂柱渗透试验提出渗流基本规律,是整个水文地质定量研究的基石。其核心公式:v=Ki。式中v为渗流速度,K为渗透系数,反映岩土透水能力大小,i为水力梯度,即两点水头差值与渗流路径长度之比。物理含义为:在层流条件下,地下水渗流速度与水力梯度呈线性正比关系。渗透系数是模型最重要参数,受岩性颗粒大小、分选、孔隙裂隙发育程度、孔隙充填物影响;达西定律适用均质、各向同性、稳定层流渗流,在砾石粗大孔隙、岩溶管道高速紊流区不再适用,模型需要进行修正处理。

2. 水量均衡原理(质量守恒原理)

地下水系统遵循物质守恒,任意一个闭合水文地质单元,在指定时段内总补给量减去总排泄量,等于单元内储存量变化量,即:

打开网易新闻 查看精彩图片

补给来源主要包括大气降水入渗、地表水渗漏、侧向地下径流流入、灌溉回渗;排泄方式有潜水蒸发、泉排泄、侧向流出、人工开采、越流排泄。水量均衡是构建概念模型、开展地下水资源量评价最简单、最常用的原理,也是数值模型迭代计算的约束条件,保证模型水量收支闭合,避免计算失真。

3. 连续性方程与地下水基本微分方程

将达西定律代入质量守恒方程,可推导出地下水三维渗流微分方程,是数学水文地质模型的核心控制方程。区分潜水(无压水)、承压水(有压水)两类基本含水层:承压水存在弹性释水,潜水存在重力释水与自由面变动。该方程可以描述非稳定流条件下,空间任意点水头随时间的动态变化,结合边界条件与初始条件即可求解水头分布,为数值模拟提供理论内核。

4. 越流原理与含水层系统结构原理

天然地层多由含水层、弱透水层互层分布构成多层结构,相邻含水层通过弱透水层发生垂向水量交换,称为越流。当两层水头存在差值时,水流穿过弱透水层产生越流补给或排泄,越流系数成为多层含水层模型关键参数,解释层间水力联系,也是地面沉降、多层地下水开采相互干扰分析的理论基础。

二、三类主流水文地质学模型及其原理 (一)水文地质概念模型

概念模型是最基础的定性 — 半定量模型,是后续物理、数学模型构建的前提。构建原理:在野外勘探、抽水试验、水位监测、地质钻孔资料基础上,对复杂地质体进行合理概化,剔除次要细微地质细节,提炼系统边界、地层结构、含水层类型、补给径流排泄路径、水力性质等关键要素。主要概化内容:确定模型范围边界(定水头边界、定流量边界、隔水边界、混合边界);划分含水层与隔水层,确定是均质 / 非均质、各向同性 / 各向异性;梳理补给、排泄、人工开采等源汇项;区分稳定流或非稳定流运行模式。应用优势:逻辑直观、门槛低,适合前期水文地质条件分析、水资源粗算、方案定性研判;局限性是难以精细刻画水位动态与污染物长期演化,只能作为建模基础框架。

(二)物理模型(砂槽模型)

物理模型属于比例缩尺实体模型,遵循相似原理,包括几何相似、运动相似、动力相似。基本原理:按照一定比例尺缩小天然含水层,用砂、砾石、黏土等材料模拟实际地层,通过人工注水、抽水模拟补给与开采,实测水头、流量变化,还原真实渗流过程。利用相似准则,将模型实测数据换算至实际场地尺度,求解渗透系数、影响半径、涌水量等参数。典型应用:室内演示达西定律、基坑降水渗流路径模拟、岩溶管道渗流规律研究、坝基渗漏机理试验。缺点是制作周期长、成本高、尺度受限,难以模拟大范围、多层复杂含水层系统,如今更多用于教学与机理试验,工程实际中使用逐渐减少。

(三)数学模型(数值模型为当前主流)

数学模型以控制微分方程搭配初始条件、边界条件定量描述地下水运动,分为解析解模型与数值模型两大类。

  1. 解析法模型原理:在含水层均质、边界规则、简化假设前提下,直接求解地下水微分方程,得到水头、流量解析公式。经典公式如裘布依稳定潜水井公式、泰斯非稳定承压井公式,可快速计算单井涌水量、水位降深、影响半径。优点计算简便;缺点假设条件过于理想化,难以适配边界复杂、地层多变的实际场地,多用于简单抽水计算与参数求取。

  2. 数值模型(有限差分法、有限元法为主)核心原理:将研究区划分为大量离散网格单元,把连续的地下水微分方程转化为网格节点上的线性方程组,借助计算机迭代求解每个节点不同时刻水头值,实现全域动态模拟。主流算法包括有限差分法(MODFLOW 软件核心算法)、有限元法,能够适配非均质、各向异性、复杂边界、多层越流、人工开采、面状入渗等复杂现实条件。拓展衍生溶质运移模型:在水流模型基础上叠加对流、弥散、吸附、降解等溶质运移方程,模拟污水渗漏、垃圾填埋场污染羽迁移、海水入侵等水质演化问题;耦合应力场可构建地面沉降模型,刻画地下水超采引发地层压缩变形规律,是当前水文地质精细化研究的核心手段。

三、模型构建通用步骤与应用意义

完整水文地质建模遵循标准化流程:收集地质水文资料→建立水文地质概念模型→参数分区赋值(渗透系数、给水度、储水率、入渗系数等)→设定初始水头与各类边界条件→模型调试、识别与验证(拟合历史水位动态)→模型预测运算→结果分析与方案优化。模型识别验证是关键环节,通过调整水文地质参数,让模拟计算水位与长期实测水位误差控制在合理范围,保证模型可靠性,避免结果失真。

在工程与资源环境领域,水文地质模型价值突出:一是定量评价地下水资源可开采量,指导地下水合理开采规划;二是预测基坑、隧道施工涌水量,制定降水与防渗方案;三是模拟地下水污染迁移范围,设计污染修复治理方案;四是分析地下水超采引发地面沉降、岩溶塌陷、海水入侵等地质灾害演化趋势,提出管控对策;五是为矿山疏干、水库渗漏、地热资源开发提供定量依据。

结语

总而言之,达西渗流定律与水量守恒是水文地质学模型的底层原理,概念模型、物理模型、数学模型由浅入深构成完整建模体系。天然地下水系统具备隐蔽性、空间非均质性、动态多变性,无法直接精准观测全部运行细节,而水文地质模型通过科学合理概化,实现从定性认识到定量计算的跨越。伴随计算机数值模拟技术不断发展,多场耦合模型(水流 — 溶质 — 应力 — 温度耦合)成为发展趋势,能够应对愈发复杂的水资源与生态地质环境问题,既是水文地质学理论落地的重要载体,也是地下水资源科学管理、地质灾害防控不可或缺的技术工具。

打开网易新闻 查看精彩图片

地球的表层结构和水系的分布图示

打开网易新闻 查看精彩图片

SGD示意图

打开网易新闻 查看精彩图片

可能发生化学渗透现象的地质构造略图

打开网易新闻 查看精彩图片

地下水砷-氟-碘迁移转化释放多地球化学过程耦合机理

打开网易新闻 查看精彩图片

美国西南半干旱区典型森林山坡中的水均衡项

打开网易新闻 查看精彩图片

水循环中的绿水(气态水)和蓝水(液态水)示意图

打开网易新闻 查看精彩图片

水文循环是指水通过地球的地壳、大气、海洋、湖泊和溪流的运动

打开网易新闻 查看精彩图片

岩石中的孔隙度取决于多个因素。在砂岩中,胶结程度和颗粒分选程度都很重要

打开网易新闻 查看精彩图片
打开网易新闻 查看精彩图片
打开网易新闻 查看精彩图片

喀斯特地形的一些主要特征包括洞穴、天坑和消失的溪流

打开网易新闻 查看精彩图片
打开网易新闻 查看精彩图片
打开网易新闻 查看精彩图片
打开网易新闻 查看精彩图片
打开网易新闻 查看精彩图片
打开网易新闻 查看精彩图片
打开网易新闻 查看精彩图片
打开网易新闻 查看精彩图片
打开网易新闻 查看精彩图片