1. 课程导入与结果判读总览

FEFLOW 地下水流动模拟的价值,不在于模型能算出一组数字,而在于这组数字能否支撑一个可信的水文地质结论。很多初学者把注意力集中在建模阶段——网格剖分、参数赋值、边界条件设置,却忽略了模拟结束后的结果判读环节。实际上,从模型输出到可信结论之间,还隔着一段需要专业判断的路程:水头场是否合理、流速方向是否符合区域地下水流向、流线是否与已知补给排泄关系一致、均衡项是否闭合。本课程要解决的,正是这段"最后一公里"的问题。

结果图的核心类型

FEFLOW 的结果可视化围绕几个基本物理量展开,每一类图承担不同的判读任务,不能互相替代。

水头(Hydraulic Head)图是最基础的结果图。等水头线图反映势能场的空间分布,是判断区域流场形态、识别补给区与排泄区的主要依据。判读时关注等值线间距(间距小表示水力梯度大)、闭合等值线(可能对应局部抽水或注水)、以及等值线与已知水位观测点的吻合程度。水头图通常配合色标渲染,但色标范围的选择会显著影响视觉判断,需要固定为有物理意义的区间。

流速(Velocity)图表现渗流速度的大小与方向。由于地下水实际流速很小,FEFLOW 默认输出的是达西流速(单位如 m/d),判读时要注意量级是否合理。流速矢量图适合观察局部流向,流速等值线图适合识别相对强径流带。在非均质含水层中,流速图往往比水头图更能暴露参数分区的问题。

流线(Pathline)图由粒子追踪生成,直观展示水质点从补给区到排泄区的运移路径与旅行时间。流线图是连接水头场与溶质运移的桥梁,判读重点是流线的起止位置是否落在合理的补给、排泄边界上,以及流线是否穿越了不透水层等物理上不应穿越的结构。

均衡(Budget)图统计各边界与源汇项的流入流出量,是检验模型水量守恒的核心工具。均衡分析回答的是"水从哪里来、到哪里去、各项占比多少"的问题,其误差项(Error [%])是判断模型数值收敛与物理合理性的直接指标。

误差判读的整体逻辑

误差判读不是单一指标的检查,而是一条从数值到物理、从局部到整体的递进链条。合理的判读顺序是:先看数值收敛性(均衡误差是否在可接受范围),再看场变量合理性(水头、流速的空间形态是否符合水文地质认识),然后看观测拟合(计算值与实测值的偏差),最后看模型行为一致性(流线与均衡是否互相印证)。这条链条中任何一环出现明显异常,都应回溯到建模环节查找原因,而不是在结果图上强行解释。

需要强调的是,误差来源是多层次的:有数值误差(离散化、收敛容差)、有参数误差(渗透系数、给水度的不确定性)、有概念模型误差(边界条件概化是否合理)、还有观测误差(实测水位本身的测量精度)。判读时要区分误差的性质,才能决定是调整数值设置、重新率定参数,还是修正概念模型。

FEFLOW 结果判读框架:从模型输出到可信结论 模型输出 水头场 等值线 / 色标渲染 流速场 矢量 / 等值线 流线 粒子追踪 / 旅行时间 均衡 边界流量 / 源汇项 观测拟合 计算值 vs 实测值 判读递进链条 ① 数值收敛性 均衡误差 Error[%] 是否达标 ② 场变量合理性 水头形态 / 流速量级与方向 ③ 观测拟合 残差分布 / 系统偏差识别 ④ 行为一致性 流线—均衡—水头互相印证 ⑤ 可信结论 支撑水文地质判断与决策 误差来源分层 数值误差 离散化 / 收敛容差 时间步长 参数误差 渗透系数 K 给水度 / 储水系数 概念模型误差 边界条件概化 含水层结构简化 观测误差 水位测量精度 观测点代表性 区分误差性质 → 决定修正方向 判读顺序:数值 → 场变量 → 观测 → 一致性;任一环节异常应回溯建模环节

判读框架的实践含义

把上述逻辑落到实际操作上,一个典型的工作流是这样的:模型运行结束后,先打开均衡图,检查 Error [%] 是否在可接受阈值内(例如稳态模型通常要求小于 1%,瞬态模型可适当放宽,但需说明理由)。若均衡误差偏大,优先检查时间步长、收敛容差和边界条件设置,而不是急着看水头图。

均衡通过后,进入场变量检查。以某示例模型为例(以下参数为教学假设,非实测):研究区为山前冲洪积平原,渗透系数在 5~50 m/d 之间分区赋值,上游定水头边界 120 m,下游定水头边界 85 m,中间有一口抽水井。运行后水头等值线应大致垂直于区域流向、从上游向下游递减;若出现等值线在抽水井附近异常密集但影响半径明显偏小,往往提示该区域渗透系数赋值偏低或网格分辨率不足。此时流速图能进一步佐证:抽水井周围流速矢量应指向井,且量级高于外围。

流线图则用来检验补给—排泄关系是否自洽。合理的流线应从上游定水头边界或降水入渗面出发,经抽水井或下游边界终止,不应出现流线从定水头边界"凭空产生"或穿越隔水底板的情况。若流线与均衡图矛盾——例如均衡显示某边界为流入,但流线却从该边界流出——说明边界类型或参数设置有误,需要回到模型检查。

最后是观测拟合。把计算水头与观测孔实测水位对比,关注残差的空间分布而非单一数值:若残差呈系统性偏差(如全区计算值偏高 2 m),可能是边界水头或整体参数尺度问题;若残差随机分布且量级在观测误差范围内,则模型可信度较高。需要提醒的是,拟合好不等于模型对,参数补偿效应可能让错误的概念模型也获得不错的拟合,因此观测拟合必须与前面的场变量、流线、均衡检查联合使用。

判读时建议固定一套检查顺序并形成记录:均衡误差 → 水头形态 → 流速量级 → 流线起止 → 观测残差。每次调整模型后按同一顺序复查,便于对比变化、定位问题。