第一次打开Visual MODFLOW Flex界面的时候我脑子里只有一句话这跟我用了好几年的传统Visual MODFLOW完全不是同一个物种。只要你真正动手做过一个项目就会发现Flex把地下水模拟的逻辑彻底换了一遍——过去是网格驱动现在是概念模型驱动。这篇内容不想写成软件说明书而是想把我实际用Flex做完整项目的经验拆开来讲前期数据怎么整理、概念模型怎么建、网格怎么选、运行卡住怎么排查、校准为什么困住那么多人以及那些只有真正做项目才会撞上的坑。这篇文章适合三类人看准备从旧版MODFLOW迁移到Flex的工程师刚接触地下水数值模拟、想直接上手做真实项目的学生以及被模型不收敛、属性映射失败折磨得头痛的实践者。我不会绕弯子下面直接进入正题。1. 先搞清楚Flex到底换了什么内核从网格驱动到概念驱动1.1 传统MODFLOW建模的痛点在哪里先说个背景。早期用Visual MODFLOW 4.x这样的软件工作路径非常固定人工确定网格行列数然后在每一个网格单元格上手动或刷子式地给渗透系数、储水率、边界条件。换句话说你先有网格后有水文地质概念。这在模型范围小、分层简单的项目里问题不大但一旦遇到复杂的实际场地就有几个明显难受的地方改一次模型边界网格可能要全部重新剖分之前赋的参数跟着丢失或者错位。很多水文地质信息本来就是以地图、钻孔、分区的形式存在硬拆到网格单元上非常费劲。几何边界复杂时矩形网格要么把边界锯齿化要么为了精度疯狂加密计算量暴涨。我最早用老版本做一个河岸边坡的渗流分析光是处理河道边界就花了大半天。河道是斜的矩形网格怎么剖都不对劲河床高程没法合理赋值最后只能妥协用阶梯状边界结果流场形态看起来总是很怪。这类问题不是参数调不好的问题而是工作流本身别扭。1.2 Flex的解题思路先在概念层建模后生成网格Flex的核心逻辑是反过来的你第一步做的不是网格而是用地图图层的语言把整个含水层系统描述清楚。这个描述包括模型边界是什么形状地层怎么分层哪些区域渗透性好、哪些区域渗透性差哪里有河流排水沟哪里有井降水补给和蒸发怎么分布。所有这些都用点、线、面这些地理对象来表达。概念模型建好之后再由软件自动生成网格把概念模型里的属性映射到每一个网格单元上。这样带来的直接好处是模型边界改了直接改概念图层再重新生成网格就行不用在网格上重干。概念模型可以复用同一个分区图、同一套边界条件可以快速生成不同剖分精度的多个网格方案。和ArcGIS这类GIS工具的亲和度很高很多空间数据不用二次加工直接拖进Flex用。从数据流上看Flex的工作路径大概是这么一条线钻孔数据 / DEM / Shapefile / Excel → 概念模型图层、边界条件、分区属性 → 网格生成QUAD/TRIANGLE/嵌套 → 属性映射与数值模型 → MODFLOW求解 → 后处理水头、流线、粒子追踪 → 成果输出理解了这条线后面所有操作都是在给这条线填内容。很多人学Flex觉得混乱就是因为他们还在用老版本的网格优先思维去看界面自然找不到北。提示如果你是从老版本转过来的我建议第一次打开Flex后先别急着建模花半天时间把地图视图、模型浏览器、剖面视图这几个核心面板摸一遍。Flex的操作高度依赖鼠标右键和面板拖拽和传统版本的点菜单完全不同。1.3 界面布局与核心视图Flex的界面大致分几个区域左侧是模型浏览器中间大面积是地图视图底部有信息日志右键菜单提供大量操作入口。你还会经常用到剖面视图和三维视图查看地层结构、边界条件的空间关系非常直观这些视图联动是建模、调试的好帮手。实际建模型的时候我的习惯是保持地图视图为主需要验证地层起伏时切到剖面视图需要检查整体感时切换三维视图。三者在同一模型数据上是联动的鼠标选中某个对象其他视图里也会同步高亮排查问题的时候效率很高。2. 建模前的数据体检这些前期工作决定成败很多人拿到Flex几天就能把界面玩熟但真正开始做项目时卡住的往往是数据。地下水数值模型是典型的数据密集型工作资料不全、数据质量差软件再先进也白搭。我一般会在正式建概念模型之前花一到两天做数据体检宁可这段时间慢一点也好过后期返工。2.1 至少需要准备哪几类数据按照Flex的建模逻辑我把数据分成四大类第一类地形与地表水系数据。主要用DEM数字高程模型来定义地表高程投影到模型顶部。河流、湖泊、海岸线这些用Shapefile或者KML提边界即可。如果项目区内地表高程变化不大可以只用少量高程点插值地形起伏明显的DEM分辨率太粗会直接影响地表与地下水力联系的刻画。第二类钻孔与地层分层数据。这是整个模型的地质骨架。需要搜集所有可用的钻孔柱状图每孔至少要有孔口高程、各岩性分层的顶底深度或高程。如果这些数据整理成Excel或CSV导入Flex非常方便同时如果能把岩性描述给全Flex能辅助识别层位。第三类水文地质参数数据。包括通过抽水试验得到的渗透系数、储水系数室内实验得到的孔隙度、给水度压水试验数据也可以作为参考。这些参数通常不是每个分区都有实测值很多时候要靠水文地质分区图、经验值来补充。第四类边界水文数据。比如河流水位、湖水位、排水沟水位、泉流量、降水补给量、蒸发量、开采井流量与水位等。做瞬态模拟还需要这些数据的时间序列。项目里我还做过一个简化版清单每次建模前逐项打勾数据类别关键项缺失后果地形数据DEM、地表高程点、水系shapefile地表起伏失真地表-地下水交换错误地层数据钻孔分层、剖面图、地质图分层错乱模型地质骨架不可信参数数据K、Ss、Sy、孔隙度、分区图参数全靠猜校准无从下手边界数据水位、流量、降雨、蒸发序列边界条件失真模拟结果整体跑偏开采数据井位、开采量、开采层位、井结构井周水头异常漏斗形态错误2.2 数据预处理坐标系、高程基准、缺失项测量数据要是坐标系不统一建模的时候是最容易翻车的。比如一部分图层是WGS84经纬度一部分是CGCS2000投影坐标直接导入Flex图层之间会错位得离谱。这个问题我在后文实战坑位复盘里专门讲这里先说处理原则所有空间数据导入Flex前统一到一个坐标系推荐投影坐标系米单位统一。高程基准要一致多数项目默认1985国家高程基准不同来源的数据高程如果基准不同需要先行转换。钻孔数据里如果有孔口高程缺失根据周边地形做插补分层信息不完整的孔宁可暂时不用也不要硬编。还有一个很常见的工程问题CAD图纸里的勘探点、钻孔信息要去图中提取。我的做法是把钻孔信息手工整理到Excel模板里包括编号、坐标、孔口高程、各层顶底板深度。截取部分地层后CAD图纸导入会很麻烦倒是直接在Excel里整理完再用Flex的Borehole导入更可靠。Flex对CSV格式的钻孔支持友好模板字段是固定的注意列名对应即可。2.3 常用导入格式与几类格式坑Flex支持很多格式我常用的有点数据CSV/TXT导入钻孔、井、水位观测点线/面数据Shapefile、KML、GeoDatabase栅格数据DEM可以直接带为地表高程参考我踩过的格式坑主要有三个。一是Shapefile的字段编码属性表里有中文时偶尔会乱码建议用英文短字段名。二是KML如果来自Google Earth坐标系是WGS84经纬度和投影坐标的项目叠加时务必转换。三是CSV里的经纬度坐标如果用了度分秒格式比如29°4530Flex读不了必须转成十进制度数。提示数据预处理这件事我建议在ArcGIS/QGIS里提前完成一批工作。比如把CAD里的钻孔点批量转成CSV把纸质地质图抠出来的分区多边形转成Shapefile把河道中心线缓冲成面。这些工作在GIS里做更灵活别试图在Flex里完成全部数据编辑。3. 概念模型的构建用图层语言描述一个含水层系统数据准备好了就可以建概念模型了。这个阶段是整个Flex建模过程中最考验水文地质基本功的部分概念模型层次设置得好不好直接决定后面数值模型的合理性和计算稳定性。3.1 建立分层框架Layers与Horizon的定义逻辑Flex里建立分层的逻辑和传统MODFLOW不太一样。传统MODFLOW的层是离散的网格层Flex则是先在概念模型里定义每一层的顶底板。我的做法分三步。第一步根据钻孔分层归纳水文地质层组。比如一个典型的第四系含水层结构从上到下可能是包气带弱透水、浅层含水层、弱透水隔水层、深层承压含水层。第二步确认每一层的控制点数据。浅层含水层顶板用DEM或地表高程控制底板用钻孔里对应的层底深度换算成高程。如果某个层在某些钻孔中缺失尖灭需要用特殊方式处理否则插值的时候会把不存在的层变成零厚度或负厚度模型完全没法跑。第三步在Flex中创建Layer并指定它的Top和Bottom来源然后让软件用插值算法生成整个模型范围内的层顶底高程面。分层框架建好后务必用剖面视图检查一遍。我见过太多模型平面图看着正常切了一条剖面后发现层位互相交叉、厚度出现负值这就是控制点数据没处理好。早期多花十几分钟检查剖面能省后面几天排查错乱的时间。3.2 给分区分参数水文地质属性的赋值逻辑分层的几何骨架建立之后要给不同的空间范围赋予水文地质参数。Flex处理参数赋值的方式非常接近GIS思维你创建多边形区域这个区域内部赋予一组Kx、Ky、Kz、Ss、Sy、有效孔隙度等参数。关键点是分区要有水文地质依据。一个常见的误区是直接根据数值模拟的需要把整个模型区人为切出二十个参数分区每个分区各给一套K值参数个数远远超过观测数据能支撑的辨识能力。最终模拟出来的结果即使拟合得好看参数可信度也很低。我的实践习惯是优先用岩性和成因分区。同一个沉积相或同一个含水层类型先认为参数基本一致只在有实测抽水试验数据的区域单独区分。参数赋值参考区域水文地质图先用该区域已有的成果作初始值。如果模拟后发现局部水位差异明显再在校准阶段考虑细分。各向异性方面水平方向Kx和Ky我通常不区分太多垂直方向Kz按Kx的十分之一到百分之一的量级估。这个经验值不代表一切但作为起跑点是稳定的。3.3 定义边界条件河流、排水沟、定水头、井、补给蒸发边界条件的定义是整个概念模型最核心的环节。边界条件定歪了模型跑得再收敛也是错误的结果。Flex里常用的边界条件类型有River河流需要给定河底高程、河水位、河床导水系数。河流层位要选对赋值给对应含水层。Drain排水沟只排不补给定排水基底高程和导水系数。Specified Head定水头直接给定某个面上的水位值适合处理水头稳定的大湖、大范围定水头补给边界。GHB通用水头边界以某个远处稳定水头与模型边界之间的水力联系来刻画边界适合描述模型区外的补给能力。Well井给定位置、层位、抽水量/注水量。Recharge补给通常给降水入渗补给量值远小于降水量是入渗系数和降水量乘积后的有效入渗量。EVT蒸发蒸腾浅水位地区必须考虑给定最大蒸发速率和极限蒸发埋深。这里特别想强调一个原则边界条件的空间范围要足够大不要切在主要控制因素中间。比如说研究区存在一条主要河流控制着地下水流场模型边界就不能离河道太近否则河流和边界的交互关系会在模型中被简化掉。我习惯把模型范围扩到分水岭或者稳定的补给区让重点研究区远离人为设定的边界。补给量是另一个常见的问题点。很多人直接把年降水量作为Recharge输进去这大错特错。实际入渗补给量是降水量扣除蒸散发和地表径流之后的量数值上远小于降水量干旱半干旱地区甚至只有降水量的百分之几到十几。我建议依据文献或者实际调查先给一个量级合理的初值校准阶段再进一步调整。边界条件的操作层面Flex里都是先创建点/线/面形状再配属性的逻辑河道做成多边形然后在属性表里输入河流阶段参数井做成点直接指定井所在图层和抽水量。习惯之后这种方式比老版本在网格上描线要快得多尤其是河道弯曲复杂的项目优势非常明显。3.4 剖面检查概念建模不可跳过的验收动作每个概念模型阶段完成重要步骤我都会在剖面视图里沿着几个关键方向切剖面检查。切割剖面会展示地层的顶底板起伏、边界条件的位置关系、井滤水管的位置是否落在目标含水层中间。剖面检查是减少后期返工最有效的动作之一。检查重点包括层位之间是否交叉井的滤水管段是否全部落在该井指定的地层深度内河流、排水沟底部高程是否与该层位匹配补给、蒸发区域是否覆盖合理范围。4. 网格生成与属性映射从概念到数值的关键一跳概念模型完成后下一步是Flex的Generate Grid。很多人在这一步开始放飞自我觉得网格越密越好其实这是个值得停下想两分钟的选择题。4.1 三种网格结构的选择Flex支持结构化和非结构化网格主要包括QUAD四边形、TRIANGLE三角形、VORONOI泰森多边形等类型。我的经验是QUAD网格结构化四边网格适用于地形相对规整、边界接近矩形的模型。收敛可靠生成时间快很多传统MODFLOW求解器支持良好是我的首选。即便研究区边界比较复杂也可以在重点区域用局部加密嵌套。TRIANGLE网格三角网格适用于边界很不规则或者内部有斜向构造、河流等复杂几何要素的情况。三角网格能贴合复杂边界但如果三角形形状质量差过于尖锐或太拉长会严重影响求解稳定。VORONOI网格常搭配MODFLOW-USG使用在一些特定场景如径向流模拟、井周加密、放网格单元平缓过渡表现好但上手门槛要高一些不是所有项目都需要。一个实际的取舍经验如果项目范围面积很大、边界接近规则那么QUAD网格够用不必上非结构化。如果模型范围沿着河谷呈长条状或者河道边界极端曲折再考虑三角形或者VORONOI。网格的选择不是越花哨越好而是求解稳定、管理方便更重要。4.2 网格加密与嵌套模型空间存在明显的关注梯度差异时需要网格加密。最典型的场景是抽水井附近远井场区域水位梯度小用几百米的大网格没问题但井周围几十米内水位变化剧烈如果不加密MODFLOW计算出来的降深会明显偏小无法正确刻画降落漏斗。我的处理思路是先给模型设置一个基础网格尺寸如100米保证整体计算规模可控。在井、河流、污染物运移重点关注区设置局部加密如10米或20米。网格加密范围覆盖从井中心向外至少3到5个加密网格的距离让梯度过渡平缓。要注意的是网格加密跨度不宜过大。相邻网格面积比过于悬殊会导致数值计算精度下降一般控制相邻网格尺寸比不超过1.5到2倍。从100米直接跳到10米中间最好有30米、50米的过渡。4.3 网格质量检查怎么看Flex生成网格之后可以用质量检查功能查看偏斜度Skewness、长宽比Aspect Ratio等指标。这些指标说得直白一点就是看网格单元正不正。太歪的单元会影响求解精度和收敛性。偏斜度越接近0说明单元越正但如果网格偏斜度很高通常接近1就要调整加密方式或者改用其他网格类型。长宽比太大说明网格被拉得太长某些方向数值计算误差会放大。我做检查的时候比较随意但不含糊重点项目会把质量不好的网格区域高亮显示出来结合地形地貌判断问题出在哪比如河道边界、山脚地形突变处通常是网格质量问题的多发区。提示网格生成不一定要一次到位。先跑一个粗糙的模型验证概念模型逻辑是对的确认没问题后再加密网格做正式运行。这样模型调参速度快很多避免频繁改网格导致的工作量爆炸。4.4 属性映射到数值模型网格生成后Flex会自动把概念模型里的属性参数分区、边界条件等转换到网格单元。这一步通常不需要手动操作但会发生两种常见问题。第一种是概念图层没有被传到网格。原因通常是该图层的对象没有正确覆盖到网格区域比如边界条件多边形范围有一部分落在模型边界之外Flex不会映射这部分。第二种是部分网格单元属性为空常出现在网格和概念图层形状边界没有重合的网格边缘有时候会造成局部参数为0运行直接报错或者结果明显异常。我的排查习惯是属性映射完成后在模型浏览器里展开每个属性组检查有没有未赋值的网格单元尤其是模型边界附近的零散单元。一旦发现要么调整网格边界要么修正概念图层不能让这种低级问题拖到运行阶段。5. 运行、调试与收敛把模拟真正跑起来的实战路径网格和属性就绪之后正式运行MODFLOW求解。这里要说的不是怎么点Run按钮而是运行过程中会碰到的问题怎么处理。5.1 运行前自检清单每次运行之前我大概花三分钟做一次快速自检省得跑出一个完全没法用的结果再来排查所有活动网格是否都有完整的K、Ss/Sy参数。所有井是否都位于活动网格内井口高程和滤水管位置是否合理。边界条件类型、层位、参数是否完整有没有遗漏的河流/排水沟断面。初始水头是否合理是否与边界条件水位逻辑自洽。时间步长设置是否合理瞬态模拟的步长不能太大否则收敛困难。5.2 求解器与迭代参数配置Flex里集成的求解器包括MODFLOW-2005、MODFLOW-NWT、MODFLOW-USG等。我建议普通稳定流和瞬态模拟优先考虑MODFLOW-NWT它在处理非线性和干湿单元格反复变换的问题上比2005稳得多。特别是模拟区域存在大量水位波动、部分区域可能由湿变干时2005很容易出现不收敛NWT则好很多。迭代参数一般不需要改太多稳定流给合理的水位初值后用默认迭代即可不收敛时优先检查初始条件和边界条件而不是盲目增加迭代次数。瞬态模拟的时间步长要从短到长初始步长设大容易直接发散。流场剧烈变化初期时间步长要小一些等系统稳定后再逐步放大。5.3 常见报错与快速定位运行阶段最常见的报错我列几个实打实遇到过的报错特征常见原因处理策略SCREEN BOTTOM BELOW CELL BOTTOM井的滤水管底深超过所在网格底调整井的滤水管深度或重新设置井孔结构NEGATIVE THICKNESS分层顶底板交叉或厚度为负回查钻孔分层和插值结果修正层位控制点MAXIMUM ITERATIONS EXCEEDED迭代不收敛优化初始水头、边界条件、时间步长检查网格质量DRY CELL CONVERSION网格干湿转换频繁改用MODFLOW-NWT或调整网格高程WARNING: NO FLOW ACROSS BOUNDARY模型边界不活跃检查边界条件赋值层位确认边界在模型活动区内这里面的SCREEN BOTTOM BELOW CELL BOTTOM是很多新手必踩的坑。在Excel里整理井结构数据时井深写得太深滤水管底比该处地层底面还要深MODFLOW不允许井底超出该网格地层底。解决办法是在Flex的井属性表里修正滤水管深度确保井段落在含水层有效范围内。5.4 不收敛问题排查思路模型不收敛是地下水数值模拟的高频问题排查顺序很重要。我的经验是按初始条件→边界条件→网格质量→参数取值这个顺序检查。第一步看初始水头。瞬态模拟需要一个合理的初始流场如果初值离真实流场太远运行一开始水位剧烈调整很容易发散。更好的做法是先跑一个稳定流模拟得到稳态水头再把它作为瞬态模拟的初始水头。这是行业里非常稳的做法强烈推荐。第二步看边界条件。边界条件设定是否与实际补给排泄条件相符比如定水头边界的水位值是否在合理范围内河流导水系数是否设置过大/过小。第三步看网格质量。畸形网格、加密跨度太大的区域也容易导致迭代发散。第四步是参数极值问题。K值不是零或负值储水系数物理上也要合理。偶尔会看到有人用默认K1来先跑通模型这不是不行但如果模型跑通了之后忘了改回真实量级后面调参就会非常混乱。6. 参数校准与模型验证别急着信第一次模拟结果模型能顺利运行、输出水头等值线并不代表这个模型可以用来预测。没有经过校准和验证的模型只能算一个数值解其可信度需要在观测数据面前检验。6.1 手动试错法校准的路径校准这个过程新手和老手效率差距非常大。新手喜欢一次改一个参数、反复运行模型既慢又容易陷入局部调整。我更推荐下面这种方法先做参数敏感性分析搞清楚哪些参数对观测点水头影响最大。通常情况下K值、补给量、河流导水系数都是主要敏感参数而孔隙度和储水系数对稳态流场不敏感只在瞬态模拟中起作用。分区调整K值从整体开始先统一按区域水位的平均偏差方向调节K值当区域水位趋势对了之后再细调局部差异。调整补给量当区域水头梯度趋势正确但绝对值整体偏高或偏低时往往要动补给量。最后微调边界条件参数如河流导水系数、河床高程、GHB水位等这是因为这些参数往往只影响局部用来修局部偏差。这里有一个实际项目的例子供参考。某个山前冲洪积扇模型第一次跑出来山前区域水位比观测值偏高5到8米而下游平原区域水位基本吻合。我没有去动下游的K值而是判断山前K值偏小导致上游水排出太慢、壅高过多。把山前区域的K从8 m/d上调到15 m/d之后山前水位明显回落下游水位只发生了很小的变化。如果一上来就动全模型参数很可能就把本来合理的地方改坏了。6.2 PEST自动校准的使用建议Flex集成了PEST自动校准模块可以大幅提高效率。PEST的本质是根据观测数据自动调整参数使模拟值和实测值之间的误差最小化。使用PEST的流程是选择待校准参数、设定参数取值范围、添加观测点水头数据、运行校准。但别忘了PEST只是在给定参数结构和敏感度关系下寻找最优解它不能替你做水文地质判断。如果参数分区本身划分不合理或者观测数据不充足PEST得到的结果可能是局部最优甚至得到物理意义不合理的参数值比如K取到极其离谱的量级校准结束后务必人工检查一遍参数取值范围是否在合理区间内。新手建议这样组合先用敏感性分析筛选参数再用手动试错法确定大方向最后用PEST做精细优化。一上来就全自动PEST很容易陷入盲目信任自动结果的境地。6.3 校准效果怎么评价校准的好坏不能只看模拟值和实测值拟合得好不好看还要看误差分布是否合理、水位动态趋势是否一致。常用的评价指标包括均方根误差RMSE越小越好但受观测点数量和水位变幅影响需要结合模型用途评估。平均绝对误差MAE直观反映平均偏差。相关系数R²反映模拟值和实测值的相关性但不能单独作为评价标准。偏差Bias反映系统性的偏高或偏低。我的经验是区域地下水流场模拟含水层水位的平均绝对误差控制在观测水位变幅的10%以内通常就是可以接受的水平。但更关键的是误差不要集中在某一个区域——如果一个位置误差达到5米另一个却完全匹配模型空间分布逻辑大概率有问题不能简单用平均误差糊弄过去。校准完成后还要做验证在模型参数不变的情况下用另一段时间序列的观测数据进行模拟对比。如果验证期的拟合程度也达到要求模型才可用于预测。很多项目跳过验证直接做预测这是模型使用中的常见风险。7. 后处理、成果输出与实战坑位复盘模型校准验证完成后接下来是成果输出和项目复盘。这一部分看似简单但实际上很多人后处理做不好导致图上信息表达混乱输出格式不合格返工。7.1 水头、流线和粒子追踪的可视化Flex的后处理功能可以直接显示水头等值线、等值面、流线、矢量等。我推荐的标准成果图组合包括各含水层的稳定流或特定时间点的水头等值线图。三维水头面图能直观展示漏斗形态。流线图展示地下水的流向与补给排泄关系。粒子追踪MODPATH结果图可用于分析污染物迁移路径、水源地捕获区范围。粒子追踪是Flex比较出彩的功能因为它在三维视图里的可视化和交互做得挺好。做井捕获区分析时在井孔位置上反向追踪粒子设置若干终点就能把水源井的捕获范围勾出来。这对水源地保护区划分、污染风险评价非常实用输出成果也容易被非专业领导理解。7.2 成果图输出技巧很多软件后处理图表直接导出会显得粗糙。我一般会在Flex里先把视角、色带、等值线间隔调好再导出矢量或高清栅格到GIS平台进一步配图。Flex允许导出Shapefile比如等值线、流线这种模式下成果可以进入统一制图流程做到与现场勘察图、地质图的风格一致。点击导出时注意等值线的着色方案尽量选择色盲友好的色带同时标明坐标系、比例尺、图例瞬态模拟成果要注明时间点。这些细节虽然在地质报告评审中没有强制规范但一份图面干净清晰的成果图会让人飞快信任你的整个模型。7.3 实战项目中的典型坑位复盘最后把这几年做Flex项目踩过的坑集中复盘一下每一个都是我或同行项目中真实发生过的。坐标单位不统一导致的错位。有一次项目一部分钻孔是经纬度坐标一部分是当地投影坐标我没有在导入前统一处理结果在Flex里叠加后投影的钻孔直接飞到模型外概念模型完全无法建立。排查了半天才发现坐标系问题。从那以后每个图层导入之前第一件事就是检查坐标单位和范围小数点数量级不对就要警觉。分层数据不全导致的地层交叉。有个多层含水层项目少数钻孔没有揭穿深层含水层我在整理钻孔数据时把该钻孔对应的深层层位留空结果Flex插值后深层含水层顶板高程在某些区域反而比浅层顶板高出现了负厚度。这种问题在剖面检查里一眼就能发现但如果不做剖面检查模型运行后结果就完全失真。边界范围切到了主要补排带上。有一个河流附近的项目为了让计算量小一点我把模型边界直接沿河切过去还想当然地认为河道可以简化成定水头边界。结果怎么调参都拟合不了河岸附近的水位。后来把模型范围向河道对岸扩展了约2公里把河道完整放进来用River边界条件而不是边界截断模拟结果马上就合理了。模型边界位置的影响往往比参数取值更敏感。补给量取值严重偏大。还有个半干旱区的模型第一次跑出来全区水位整体抬得很高就是因为在Recharge用了年降水量的30%——这个比例在湿润区也许行在半干旱区明显过大。查阅当地入渗补给研究成果后把补给量调到降水量的5%到10%区间水位分布才回到合理范围。补给量这类参数容易被忽视但它对区域水位的整体水平影响极大。开采井用实测水位却没有给定流量。还有一次比较特殊的坑别人交付的模型里井边界条件写的是Specified Head井水位而不是抽水量。这带来的问题是井附近的流场被这个固定水位锚住了无法反映真实抽水对水位的影响做开采方案预测时得到的漏斗形态完全不可信。地下水流场模拟里的开采井绝大多数情况下应该用流量边界Q除非井内完整观测水位作为控制条件有特定用途。回到项目实操层面Flex这套概念模型工作流确实把建模效率提升了一个量级但软件永远是工具模型是否可信最终还是取决于你对水文地质条件的认识、对数据的判断、以及是否在做每一步时都追问这个选择在地质上合理吗。建模不难把模型建到位才真正考验功底。
