1. 课程导入与结果判读总览
FEFLOW 地下水流动模拟的价值,不在于模型能算出一组数字,而在于这组数字能否支撑一个可信的水文地质结论。很多初学者把注意力集中在建模阶段——网格剖分、参数赋值、边界条件设置,却忽略了模拟结束后的结果判读环节。实际上,从模型输出到可信结论之间,还隔着一段需要专业判断的路程:水头场是否合理、流速方向是否符合区域地下水流向、流线是否与已知补给排泄关系一致、均衡项是否闭合。本课程要解决的,正是这段"最后一公里"的问题。
结果图的核心类型
FEFLOW 的结果可视化围绕几个基本物理量展开,每一类图承担不同的判读任务,不能互相替代。
水头(Hydraulic Head)图是最基础的结果图。等水头线图反映势能场的空间分布,是判断区域流场形态、识别补给区与排泄区的主要依据。判读时关注等值线间距(间距小表示水力梯度大)、闭合等值线(可能对应局部抽水或注水)、以及等值线与已知水位观测点的吻合程度。水头图通常配合色标渲染,但色标范围的选择会显著影响视觉判断,需要固定为有物理意义的区间。
流速(Velocity)图表现渗流速度的大小与方向。由于地下水实际流速很小,FEFLOW 默认输出的是达西流速(单位如 m/d),判读时要注意量级是否合理。流速矢量图适合观察局部流向,流速等值线图适合识别相对强径流带。在非均质含水层中,流速图往往比水头图更能暴露参数分区的问题。
流线(Pathline)图由粒子追踪生成,直观展示水质点从补给区到排泄区的运移路径与旅行时间。流线图是连接水头场与溶质运移的桥梁,判读重点是流线的起止位置是否落在合理的补给、排泄边界上,以及流线是否穿越了不透水层等物理上不应穿越的结构。
均衡(Budget)图统计各边界与源汇项的流入流出量,是检验模型水量守恒的核心工具。均衡分析回答的是"水从哪里来、到哪里去、各项占比多少"的问题,其误差项(Error [%])是判断模型数值收敛与物理合理性的直接指标。
误差判读的整体逻辑
误差判读不是单一指标的检查,而是一条从数值到物理、从局部到整体的递进链条。合理的判读顺序是:先看数值收敛性(均衡误差是否在可接受范围),再看场变量合理性(水头、流速的空间形态是否符合水文地质认识),然后看观测拟合(计算值与实测值的偏差),最后看模型行为一致性(流线与均衡是否互相印证)。这条链条中任何一环出现明显异常,都应回溯到建模环节查找原因,而不是在结果图上强行解释。
需要强调的是,误差来源是多层次的:有数值误差(离散化、收敛容差)、有参数误差(渗透系数、给水度的不确定性)、有概念模型误差(边界条件概化是否合理)、还有观测误差(实测水位本身的测量精度)。判读时要区分误差的性质,才能决定是调整数值设置、重新率定参数,还是修正概念模型。
判读框架的实践含义
把上述逻辑落到实际操作上,一个典型的工作流是这样的:模型运行结束后,先打开均衡图,检查 Error [%] 是否在可接受阈值内(例如稳态模型通常要求小于 1%,瞬态模型可适当放宽,但需说明理由)。若均衡误差偏大,优先检查时间步长、收敛容差和边界条件设置,而不是急着看水头图。
均衡通过后,进入场变量检查。以某示例模型为例(以下参数为教学假设,非实测):研究区为山前冲洪积平原,渗透系数在 5~50 m/d 之间分区赋值,上游定水头边界 120 m,下游定水头边界 85 m,中间有一口抽水井。运行后水头等值线应大致垂直于区域流向、从上游向下游递减;若出现等值线在抽水井附近异常密集但影响半径明显偏小,往往提示该区域渗透系数赋值偏低或网格分辨率不足。此时流速图能进一步佐证:抽水井周围流速矢量应指向井,且量级高于外围。
流线图则用来检验补给—排泄关系是否自洽。合理的流线应从上游定水头边界或降水入渗面出发,经抽水井或下游边界终止,不应出现流线从定水头边界"凭空产生"或穿越隔水底板的情况。若流线与均衡图矛盾——例如均衡显示某边界为流入,但流线却从该边界流出——说明边界类型或参数设置有误,需要回到模型检查。
最后是观测拟合。把计算水头与观测孔实测水位对比,关注残差的空间分布而非单一数值:若残差呈系统性偏差(如全区计算值偏高 2 m),可能是边界水头或整体参数尺度问题;若残差随机分布且量级在观测误差范围内,则模型可信度较高。需要提醒的是,拟合好不等于模型对,参数补偿效应可能让错误的概念模型也获得不错的拟合,因此观测拟合必须与前面的场变量、流线、均衡检查联合使用。