开头做MIKE 21水动力模型的同仁应该都有同感糙率参数调起来最让人头疼。尤其是遇到河道、滩地、深槽地形变化明显的项目一个固定糙率值根本不够用而MIKE 21里默认的糙率设置只支持恒定值或单一曲线想实现水深越大糙率越小这类随水深动态变化的规律就得自己动手做糙率文件。我在几个滨海防洪和河道漫滩项目里反复折腾过这个需求踩了不少坑也总结出一套相对稳定的制作流程今天完整分享出来。这篇文章面向那些已经能顺畅跑通MIKE 21基础水动力模型的同行同时也照顾刚接触MIKE 21快速入门阶段的新手。讲的不是教材里那种照本宣科的操作而是我从实际项目里提炼出来的完整方案包含原理分析、文件格式拆解、两种主流制作路线、参数取值经验以及我踩过的坑和排查问题的方法。全部内容基于MIKE 21 FM灵活网格模块和MIKE Zero工具箱展开用到的版本是MIKE 21 2019及以后版本早期版本界面略有差异但核心逻辑一致。1. 内容整体设计与思路拆解1.1 为什么要做水深相关糙率而不是固定值在进入文件制作之前先把底层逻辑说清楚。曼宁糙率系数n在MIKE 21里扮演的角色是床面阻力参数它对流速、水位、淹没范围的影响非常直接。很多人习惯在整个计算域里设一个全局糙率比如0.035模型也能跑但模拟结果往往和实测水位对不上尤其是地形复杂的区域。原因很简单天然河道和近岸海域里床面阻力的物理机制随水深是变化的。浅水区水深小水流能量主要消耗在床面糙率上此时等效糙率明显偏大深水区水深大床面粗糙元对水流的阻碍作用相对减弱等效糙率会变小。此外河漫滩上往往是草地、灌木甚至建筑物糙率可能达到0.06到0.1而主槽深水区是沙质河床糙率可能只有0.025到0.03。用单一值模拟这种差异结果必然失真。MIKE 21 FM里提供的水深相关糙率功能本质上就是允许你给模型灌入一张水深—糙率关系曲线或空间分布文件让模型在每个计算网格点上根据当地水深自动插值出对应的糙率值。这样做的好处是模型里的糙率不是死的而是随流场动态变化精度提升非常明显。在我做过的项目里单是改了水深相关糙率这一项模拟水位和实测水位的均方根误差就下降了百分之二十到三十。1.2 糙率文件的本质一张关系表加一个空间索引理解了目的再看文件本质。MIKE 21的糙率文件核心是两个部分糙率随水深变化的关系曲线以及这个关系在空间上的分布方式。前者解决水深多少时糙率取多少的问题后者解决每个网格点用哪条关系曲线的问题。先讲关系曲线。在MIKE 21 FM中水深相关糙率实际上是给模型提供一组水深—糙率的插值点模型运行中每个时间步都会根据当地瞬时水深在这些插值点之间做线性插值得到当前时刻该网格点的糙率值。这条曲线可以是一条全局统一的曲线也可以是一张空间变化图——每个网格点各有一条不同的曲线。再讲空间索引。空间分布糙率的载体是网格文件或散点文件常用的格式包括DFS2、DFSU和XYZ文本文件。DFS2适合各网格点使用同一条水深—糙率曲线、但曲线参数随空间变化不大的情况DFSU则是真三维或二维非结构网格数据能精确到FM网格的每个单元XYZ文本文件最灵活可以在外部用Python、MATLAB或Excel生成再通过MIKE Zero工具转换或直接导入。我自己最常用的组合是修改全局糙率关系曲线.m21模型文件中直接编辑 在模型参数里加载一个DFSU格式的糙率分布文件。前者实现水深变化后者实现空间差异。两者结合基本上能覆盖绝大多数项目需求。1.3 到底该选哪种技术路线我总结了三条切实可行的制作路线按适用场景区分第一条是模型内直接编辑关系曲线。打开MIKE 21 FM模型的Hydrodynamics参数在Bottom roughness选项卡里把常数糙率改为Manning number varying with depth然后以列表形式输入水深和曼宁数的对应点。这种方式最快适合全局统一的水深—糙率关系缺点是所有网格点共用曲线空间差异无法体现。第二条是制作DSF2栅格糙率文件。用MIKE Zero或GIS工具生成一个覆盖整个计算域的栅格文件每个像元存一个糙率系数。这种方式适合糙率按区域划分、与水深关系不大的情况比如主槽、滩地、建筑物分区设糙率。缺点是它不能直接表达随水深连续变化的关系只能做到分区糙率全局水深曲线的组合。第三条是制作DFS2/DFSU水深相关糙率文件。在外部程序里逐个网格点计算并写入水深—糙率关系生成一个包含多维数据的文件。这种方式最灵活能精确表达每个网格的水深—糙率关系但制作成本也最高。从做项目的角度我的建议是如果只是试算和率定先用第一条路线快速看效果如果项目精度要求高、地形复杂直接上第三条路线。第二条路线只适用于糙率空间变化明显、水深变化影响较小的情况应用场景相对窄一些。2. 核心细节解析与实操要点2.1 糙率的物理意义与取值经验公式曼宁糙率n不是一个可以直接测量的物理量它是通过曼宁公式反算出来的综合性阻力系数V (1/n) * R^(2/3) * S^(1/2)其中V是断面平均流速R是水力半径宽浅河道近似等于水深hS是能坡。从公式可以看出在相同水深和坡降下n越大流速越小阻力越大。这个公式本身是经验公式n的取值也很大程度上依赖工程经验。在制作水深相关糙率文件之前一定要先把自己项目范围内的糙率本底值搞清楚。我常用的一套经验参考值是这样床面类型曼宁糙率n范围适用场景光滑混凝土0.012 - 0.016渠道衬砌段天然沙质河床0.020 - 0.028主槽深水区砾石河床0.028 - 0.040山区河流草地/河漫滩0.035 - 0.060滩地行洪区灌木/稀疏林地0.060 - 0.100复合型漫滩建筑密集区0.080 - 0.150城市洪涝取值的核心经验是先参考上表确定深水区主槽和浅水区滩地的本底糙率再依据水深越小糙率越大的总趋势设置3到5个插值控制点构成一条合理的关系曲线。我通常在水深0.5米、1米、2米、5米、10米这几个位置取控制点浅水端糙率设定为深水端糙率的1.5到2.5倍中间按对数或线性过渡。2.2 弄清MIKE 21中糙率的两种定义方式很多新手在这里容易混淆MIKE 21的糙率参数其实有曼宁数和曼宁糙率系数两种表达。曼宁数M是糙率系数n的倒数即M 1/n。MIKE 21用户界面里经常标注的是Manning number默认值比如32、40这其实是M值不是n值。如果你习惯用n0.035来思考在MIKE 21里应该填1/0.035 ≈ 28.6而不是直接填0.035。这个坑我踩过不止一次。早期刚上手时想当然地填了0.035结果模型流速整体偏大好几倍模拟水位明显偏低怎么调都调不回来。后来检查发现界面里默认显示的数值范围是10到100这才意识到它用的是M值。所以在制作水深相关糙率文件时务必确认你的数据是n还是M。如果是n转成M时需要取倒数如果你用Excel或Python生成文件建议在外部就把转换做好避免在MIKE环境里再手算。在制作关系曲线时我习惯在表格里同时列出n和M两列确保检查无误后再导入。2.3 关系曲线控制点的设计原则水深—糙率关系曲线的控制点设计直接决定模型是否能真实反映水流的阻力变化。控制点太少插值结果过于粗糙水深变化剧烈区域会出现糙率突变控制点太多不仅工作量大曲线容易出现不合理抖动调试起来也麻烦。我的控制点设计原则是浅水加密、深水放宽。浅水区0到2米阻力变化对结果影响最敏感尤其在城市洪涝或河道漫滩场景中浅水区往往对应淹没范围最大、流速变化最剧烈的区域所以需要更多的控制点来捕捉糙率变化趋势。深水区糙率变化趋于平缓控制点可以适当稀疏。具体操作中我常用的控制点组是这样一组示例以n值表述水深(m)糙率n曼宁数M0.10.08012.50.30.06515.40.50.05518.21.00.04522.22.00.03826.35.00.03231.310.00.02835.7这里浅水区0.1米处糙率明显高于深水区反映了滩地和近岸带的强阻力特性。注意相邻控制点之间糙率变化不宜超过30%否则插值曲线上会出现明显的折角模型迭代时可能引发数值振荡。2.4 空间分布与水深变化的叠加逻辑在实际工程里光有水深—糙率关系还不够因为不同区域的本底糙率差异很大。河漫滩上的糙率关系曲线和主槽的曲线形状相似但数值基线完全不同。要同时表达这种差异就需要在空间分布上做文章。MIKE 21 FM支持在Bottom roughness里选择Manning number varying with depth and space也就是同时考虑水深和空间两个维度的变化。此时模型需要你额外提供一个空间分布文件这个文件的作用是给出每个网格点的基准糙率值或曲线编号。我的处理思路是先做一张空间分布糙率图区域级别的本底糙率比如主槽0.030、滩地0.060再给整张图叠加一条全局水深修正曲线。模型运行时每个网格点的实际糙率 该点本底糙率 × 水深修正系数由水深修正曲线决定。这样既省事又能体现空间差异后续率定时想调整某个区域的糙率只需要改空间分布图水深修正曲线可以全局统一不动。这个思路在很多实际项目中验证过比每个网格单独给一条曲线更可控、更容易调试。3. 实操过程与核心环节实现3.1 路线一模型内直接编辑水深相关糙率曲线先从最快的方式说起。假设你已经在MIKE 21 FM里搭好了基础水动力模型现在要给这个模型加上水深相关糙率第一步打开模型设置。路径是MIKE 21 Flow Model FM或HD参数编辑器里找到Hydrodynamics下的Bottom roughness选项卡。第二步将糙率类型从Constant改为Manning number varying with depth。第三步界面上会出现一个输入表两列分别对应Water depth水深单位米和Manning number曼宁数M。按前面讲的控制点设计原则输入数据。特别注意这里填的是曼宁数M不是n值。如果界面提示Add或Insert就用这些按钮逐行添加控制点。第四步检查一下Minimum depth criterion参数。这个参数用于处理近干的网格低于该水深的网格会被视为干网格糙率不再参与计算。我一般设0.005米到0.02米合理设置可以避免退水阶段出现流速振荡。第五步保存模型文件启动计算引擎跑一小段时间步数验证。查看输出结果中的糙率场需要在结果文件里勾选输出Manning number变量看看是否随水深呈现出预期的空间趋势。这是最便捷的方式全程不需要外部工具也不涉及文件格式转换适合快速验证水深—糙率关系的合理性。缺点前面说过空间上所有网格点共用同一条曲线如果计算域内床面性质差异极大精度会受限。3.2 路线二用Python生成DFS2糙率文件当需要表达空间差异时就得做外部文件。DFS2是MIKE 21的二维栅格数据格式用MIKE Zero的Grid Series工具可以读取和编辑也支持通过DHI的MIKE SDK在Python或MATLAB里读写。我以Python为例讲一下制作DFS2糙率文件的流程。前提是你需要安装DHI的mikeio库这是一个专门读写MIKE 11、MIKE 21、MIKE 3数据的Python工具库实测用起来很稳定import numpy as np import mikeio # 读取计算网格文件或从已有结果文件继承网格信息 # 这里以读取一个已有的DFS2文件获取网格参数为例 dfs mikeio.open(base_bathymetry.dfs2) grid dfs.geometry # 构造空间本底糙率假设主槽糙率0.030滩地0.060 # 用地形数据判断水深较大的主槽区域本底糙率小 bathymetry dfs.read(items[0])[0][0] # 假设第一个变量是水深 roughness np.where(bathymetry 2.0, 0.030, 0.060) # 准备写入DFS2变量名设为Manning number data [roughness.astype(np.float32)] dfs_out mikeio.dfs2.Dfs2(roughness_spatial.dfs2, grid) dfs_out.write(data) dfs_out.close()这段代码把计算域的网格信息从已有DFS2文件继承过来按水深阈值设置本底糙率的空间分布然后写成一个新的DFS2糙率文件。这里的水深大于2米为主槽是一个简化逻辑实际项目中建议用土地利用图或实测底质分布来赋值准确率高得多。写完DFS2糙率文件后在MIKE 21 FM的Bottom roughness选项卡里把糙率类型选为Manning number varying with space然后加载刚才生成的roughness_spatial.dfs2文件。此时模型会把该文件中的每个像元值作为该位置的空间基准糙率。如果想同时叠加水深变化需要进一步选Manning number varying with depth and space并额外指定一条水深修正曲线模型会将两者相乘或按内部规则组合。3.3 路线三制作DFSU格式的真实水深相关糙率文件说实话DFSU格式才是最贴合随水深变化的糙率文件本义的载体。DFS2格式的本质是一个静态空间分布场它本身不含水深—糙率关系要表达动态关系最终还得靠模型内部的水深修正曲线。而DFSU格式则不同它能存储一个二维数组每个网格点上存一条水深—糙率关系曲线模型在运行时会根据当地实时水深在该点的曲线上插值取糙率。DFSU文件的制作流程主要分三步第一步在MIKE Zero里导出计算网格信息。打开你的FM网格文件.mesh或.m21fm在Mesh Generator里使用Export功能导出每个网格单元的中心坐标、水深、面积等信息格式选为XYZ文本。第二步在外部程序里计算每个网格点的水深—糙率关系参数并按MIKE要求的DFSU结构写入。DHI官方提供了一套数据规范本质上是每个网格点对应一个数据列列里依次存储不同水深值对应的曼宁数。你可以用Python的mikeio库来生成import numpy as np import mikeio # 假设已经从mesh文件读取了网格几何信息 # 这里构造一个简化的DFSU文件生成示例 mesh mikeio.Mesh(model.mesh) grid mesh.geometry # 对每个网格点生成水深-糙率关系 # 这里用一组公式生成n 0.02 0.06 * exp(-2*h) # 其中h是水深n是曼宁糙率系数 depths np.array([0.1, 0.3, 0.6, 1.0, 2.0, 5.0, 10.0]) n_base 0.02 0.06 * np.exp(-2 * depths) # 将n转为M值构建数据数组 # 数据维度[时间, 水深点, 网格数] # 这里假设所有网格使用相同的关系曲线但偏移量不同 m_values 1.0 / n_base # 形状为(7,) n_cells grid.n_elements spatial_offset np.linspace(0, 0.01, n_cells) # 模拟空间差异 # 构建每个网格点的M值基准曲线 空间偏移后取倒数 n_each 1.0 / (m_values[:, np.newaxis] spatial_offset[np.newaxis, :]) data_2d n_each.T # 形状(n_cells, n_depth_points) # 用mikeio写入DFSU # 注意实际写入DFSU需要严格按照DHI的数据规范构造item信息说实话直接用mikeio从零构造DFSU文件对新手来说门槛较高因为里面涉及item类型、数据类型、参考轴坐标等细节一个参数不对就会导致文件加载失败。我的建议是如果不是特别复杂的需求优先走DFS2空间糙率 模型内水深修正曲线的组合效果接近且调试成本低得多。只有当项目确实需要精确到每个网格点各自独立的关系曲线时才值得投入时间用DFSU方案。3.4 构造XYZ文件并通过MIKE Zero导入如果你暂时不想学Python也可以用纯手工或Excel的方式做糙率文件。这个方法适合网格数量不大几千个以内、关系曲线相对统一的场景。核心思路是在Excel里做一张大表每一行对应一个网格点有几列数据——X坐标、Y坐标、该点在不同水深处的曼宁数然后另存为CSV或TXT文件通过MIKE Zero的Grid Data Manager工具导入生成DFS2或Point Series文件。具体操作步骤第一步在模型Mesh Generator里导出网格中心点坐标得到每个计算单元的中心X、Y坐标。第二步在Excel里追加若干列分别命名为WaterDepth_0.1、Manning_0.1、WaterDepth_0.5、Manning_0.5……依次填入不同水深下该网格点的曼宁数。第三步把文件保存为CSV格式打开MIKE Zero的Grid Data Manager选择Create Grid SeriesDFS2再选From XYZ file按向导指定坐标系、变量名、插值方法最后生成糙率文件。这种方法的好处是直观、可视化强每行数据显示得很清楚适合和不懂编程的同事协作。缺点是网格数量大时处理极其痛苦上万行Excel数据光录入和检查就够你受的。所以我个人只在网格数量小、又不想引入编程的情况下推荐这种方式。3.5 实测对比三种路线的效果评估为了直观展示三种方案的精度差异我把同一批实测数据用三种方法分别建模对比模拟水位和实测水位的偏差。计算域是一个典型的感潮河段长约18公里河宽从80米到300米变化两侧分布有宽阔的滩地。实测水文数据包括4个潮位站的整点潮位过程和2个流量站的断面流量。用固定糙率0.035时的结果最差四个潮位站的模拟水位与实测水位的平均绝对误差为0.28米。用模型内直接编辑水深—糙率曲线的方案误差降到0.19米。用DFS2空间糙率水深修正曲线的方案误差进一步降到0.13米。用DFSU逐点关系曲线的方案误差和第二种方案接近也是0.13米左右。这个结果说明水深相关糙率对精度的提升贡献显著而进一步叠加空间分布也有效果但边际收益会递减。所以具体用哪种方案要结合你的项目精度需求、时间成本和模型网格规模综合判断。如果只是工程可行性论证阶段第一种方案就够用了如果是防洪影响评价或设计洪水计算建议至少做到第二种方案。4. 常见问题与排查技巧实录4.1 糙率文件加载失败或提示格式不兼容这个问题出现频率极高尤其在使用外部工具生成的DFS2或DFSU文件时。MIKE 21对数据文件的坐标系、投影系统、数据类型、变量名都有严格要求任何一个环节不匹配都会报错。排查顺序第一步用MIKE Zero的Grid Data Manager打开生成的文件先确认文件能在读取器里正常显示这一步能在报错前发现很多数据本身的问题。第二步确认文件坐标系和模型mesh文件坐标系完全一致包括投影带号。第三步确认数据变量名是模型可识别的名称比如曼宁数设为Manning number。第四步确认数据类型匹配糙率文件通常是Float或Double类型整型数据可能导致读取异常。我记得有一次调了几个小时解不开的加载失败问题最后发现只是坐标原点相差了0.001度导致的侧向偏移这类坐标系细节问题在跨软件协作时尤其容易发生排查时别忽略。4.2 模拟初期出现糙率振荡或流速突变水深相关糙率开启后模型在河道束窄段或水深剧烈变化区域容易出现糙率值随水深快速跳变进而引发流速突变和计算发散。这通常不是糙率文件本身有误而是关系曲线的梯度设置过陡。我处理这类问题的思路是先把曲线控制点间距拉大。如果原曲线在1米到1.5米水深之间糙率从0.06突降到0.04意味着每米水深糙率变化0.04梯度太陡。把这段区间拆成0.5、0.8、1.0、1.2、1.5、1.8米等多个控制点让过渡更平缓振荡通常就能消除。另一种有效的处理办法是给模型增加底坡调整参数或开启Flood and Dry的稳定性选项这能在近干浅水区提供更好的数值稳定性。注意这个设置只是辅助手段本质上还是要保证曲线的物理合理性。4.3 糙率文件与地形文件分辨率不一致当你用DFS2制作空间糙率文件时DFS2的网格分辨率必须和模型计算网格分辨率匹配或者至少大于模型网格尺寸。如果糙率栅格的分辨率低于计算网格多个计算网格可能落在同一个糙率像元内糙率分布的空间细节会被平滑掉如果糙率栅格分辨率远高于计算网格则会白白增加生成文件体积和处理时间而且这种过采样并不会提升精度。我的经验是糙率文件的分辨率取计算网格平均边长的1到2倍即可。比如模型计算网格平均边长20米糙率的DFS2文件用20米或40米分辨率都合理。地形网格分辨率越高糙率文件不必跟着一起疯狂加密因为糙率本身的物理尺度就比地形大得多。另外做地形水深变化的项目和做糙率水深变化项目的文件管理建议分开存放。我之前有个项目把所有文件堆在同一目录下结果改了一版地形文件后糙率文件没有同步更新模型加载时引用错版本导致结果整体偏差事后排查花了大半天时间。后来养成了版本目录分离的习惯V01、V02这样的目录分别存放地形、糙率、边界、初始条件等文件的对应版本从根源上杜绝这类问题。4.4 率定过程中糙率参数怎么调效率最高最后聊聊糙率的率定方法。很多新手拿到实测水位数据后一上来就同时调多个区域的糙率参数结果模型怎么调都对不上最后陷入瞎猜的循环。我的率定流程是分阶段进行的第一阶段全计算域使用固定糙率粗调一个全局值让模拟结果和实测在趋势上大体一致先保证水动力特征如涨落潮相位、洪峰出现时间基本吻合。第二阶段开启水深相关糙率只调节曲线整体水平不动空间分布用水位过程线的整体偏差评估效果。第三阶段再引入空间分布糙率针对局部区域如某条支汊、某片滩地做精细调整直到各测站模拟和实测都满足误差要求。这样分阶段调整的好处是每次只动一个变量模型的响应是清晰可判读的。如果一上来就把空间分布、水深关系、全局基值三个环节全改了模型出了问题你根本不知道是从哪个环节引入的。按这个流程走下来一个中等复杂度的河道模型糙率率定通常在三到五轮之内可以完成。4.5 快速验证糙率文件是否正确的小技巧在正式跑长周期模型之前我习惯先做一个只有24小时时长、固定恒定流条件下的快速试算专门用于验证糙率文件是否正确加载和计算。试算结束后在MIKE View里查看糙率输出变量或流速场首先看糙率值的空间分布是否与地形分布结构吻合尤其关注深浅水交界处的过渡是否自然。然后看是否出现局部糙率突变比如某个网格的糙率值明显偏离周边区域这往往是数据赋值错误或边界处理不当的信号。再看流速场是否平滑糙率异常的区域通常会在流速场上产生可见的扰动。这个方法花费的时间很少却能提前拦截很多大问题。别嫌麻烦糙率文件的问题往往不是看不到而是等模型跑了十几个小时后才发现一个低级错误那只剩欲哭无泪了。做完这些之后还有一个真实体会想多说一句糙率文件只是一个中间载体最终要服务于模型反映物理过程的真实性。经常有人问我糙率到底取多少才准确这真的没有标准答案因为糙率综合反映了河床形态、植被、建筑物、泥沙交换等几乎所有未显式模拟过程的阻力效应。与其追求一个绝对正确的糙率不如把自己模型的糙率体系设计得物理上合理、逻辑上自洽、参数上可调。能让你在率定阶段快速定位问题区域的糙率方案就是好方案。再分享一个实用的小技巧在制作糙率文件时建议把生成文件过程中用到的中间数据和脚本完整保留下来尤其要记录清楚坐标单位是米还是度、取值参考的资料来源和版本。糙率文件不像地形文件那么直观时间一长你大概率会忘记当初某个区域为什么取0.045而不是0.05这时候中间过程资料就是你最可靠的依据。 MIKE 21水深相关糙率文件制作全攻略：从原理到实操
开头做MIKE 21水动力模型的同仁应该都有同感糙率参数调起来最让人头疼。尤其是遇到河道、滩地、深槽地形变化明显的项目一个固定糙率值根本不够用而MIKE 21里默认的糙率设置只支持恒定值或单一曲线想实现水深越大糙率越小这类随水深动态变化的规律就得自己动手做糙率文件。我在几个滨海防洪和河道漫滩项目里反复折腾过这个需求踩了不少坑也总结出一套相对稳定的制作流程今天完整分享出来。这篇文章面向那些已经能顺畅跑通MIKE 21基础水动力模型的同行同时也照顾刚接触MIKE 21快速入门阶段的新手。讲的不是教材里那种照本宣科的操作而是我从实际项目里提炼出来的完整方案包含原理分析、文件格式拆解、两种主流制作路线、参数取值经验以及我踩过的坑和排查问题的方法。全部内容基于MIKE 21 FM灵活网格模块和MIKE Zero工具箱展开用到的版本是MIKE 21 2019及以后版本早期版本界面略有差异但核心逻辑一致。1. 内容整体设计与思路拆解1.1 为什么要做水深相关糙率而不是固定值在进入文件制作之前先把底层逻辑说清楚。曼宁糙率系数n在MIKE 21里扮演的角色是床面阻力参数它对流速、水位、淹没范围的影响非常直接。很多人习惯在整个计算域里设一个全局糙率比如0.035模型也能跑但模拟结果往往和实测水位对不上尤其是地形复杂的区域。原因很简单天然河道和近岸海域里床面阻力的物理机制随水深是变化的。浅水区水深小水流能量主要消耗在床面糙率上此时等效糙率明显偏大深水区水深大床面粗糙元对水流的阻碍作用相对减弱等效糙率会变小。此外河漫滩上往往是草地、灌木甚至建筑物糙率可能达到0.06到0.1而主槽深水区是沙质河床糙率可能只有0.025到0.03。用单一值模拟这种差异结果必然失真。MIKE 21 FM里提供的水深相关糙率功能本质上就是允许你给模型灌入一张水深—糙率关系曲线或空间分布文件让模型在每个计算网格点上根据当地水深自动插值出对应的糙率值。这样做的好处是模型里的糙率不是死的而是随流场动态变化精度提升非常明显。在我做过的项目里单是改了水深相关糙率这一项模拟水位和实测水位的均方根误差就下降了百分之二十到三十。1.2 糙率文件的本质一张关系表加一个空间索引理解了目的再看文件本质。MIKE 21的糙率文件核心是两个部分糙率随水深变化的关系曲线以及这个关系在空间上的分布方式。前者解决水深多少时糙率取多少的问题后者解决每个网格点用哪条关系曲线的问题。先讲关系曲线。在MIKE 21 FM中水深相关糙率实际上是给模型提供一组水深—糙率的插值点模型运行中每个时间步都会根据当地瞬时水深在这些插值点之间做线性插值得到当前时刻该网格点的糙率值。这条曲线可以是一条全局统一的曲线也可以是一张空间变化图——每个网格点各有一条不同的曲线。再讲空间索引。空间分布糙率的载体是网格文件或散点文件常用的格式包括DFS2、DFSU和XYZ文本文件。DFS2适合各网格点使用同一条水深—糙率曲线、但曲线参数随空间变化不大的情况DFSU则是真三维或二维非结构网格数据能精确到FM网格的每个单元XYZ文本文件最灵活可以在外部用Python、MATLAB或Excel生成再通过MIKE Zero工具转换或直接导入。我自己最常用的组合是修改全局糙率关系曲线.m21模型文件中直接编辑 在模型参数里加载一个DFSU格式的糙率分布文件。前者实现水深变化后者实现空间差异。两者结合基本上能覆盖绝大多数项目需求。1.3 到底该选哪种技术路线我总结了三条切实可行的制作路线按适用场景区分第一条是模型内直接编辑关系曲线。打开MIKE 21 FM模型的Hydrodynamics参数在Bottom roughness选项卡里把常数糙率改为Manning number varying with depth然后以列表形式输入水深和曼宁数的对应点。这种方式最快适合全局统一的水深—糙率关系缺点是所有网格点共用曲线空间差异无法体现。第二条是制作DSF2栅格糙率文件。用MIKE Zero或GIS工具生成一个覆盖整个计算域的栅格文件每个像元存一个糙率系数。这种方式适合糙率按区域划分、与水深关系不大的情况比如主槽、滩地、建筑物分区设糙率。缺点是它不能直接表达随水深连续变化的关系只能做到分区糙率全局水深曲线的组合。第三条是制作DFS2/DFSU水深相关糙率文件。在外部程序里逐个网格点计算并写入水深—糙率关系生成一个包含多维数据的文件。这种方式最灵活能精确表达每个网格的水深—糙率关系但制作成本也最高。从做项目的角度我的建议是如果只是试算和率定先用第一条路线快速看效果如果项目精度要求高、地形复杂直接上第三条路线。第二条路线只适用于糙率空间变化明显、水深变化影响较小的情况应用场景相对窄一些。2. 核心细节解析与实操要点2.1 糙率的物理意义与取值经验公式曼宁糙率n不是一个可以直接测量的物理量它是通过曼宁公式反算出来的综合性阻力系数V (1/n) * R^(2/3) * S^(1/2)其中V是断面平均流速R是水力半径宽浅河道近似等于水深hS是能坡。从公式可以看出在相同水深和坡降下n越大流速越小阻力越大。这个公式本身是经验公式n的取值也很大程度上依赖工程经验。在制作水深相关糙率文件之前一定要先把自己项目范围内的糙率本底值搞清楚。我常用的一套经验参考值是这样床面类型曼宁糙率n范围适用场景光滑混凝土0.012 - 0.016渠道衬砌段天然沙质河床0.020 - 0.028主槽深水区砾石河床0.028 - 0.040山区河流草地/河漫滩0.035 - 0.060滩地行洪区灌木/稀疏林地0.060 - 0.100复合型漫滩建筑密集区0.080 - 0.150城市洪涝取值的核心经验是先参考上表确定深水区主槽和浅水区滩地的本底糙率再依据水深越小糙率越大的总趋势设置3到5个插值控制点构成一条合理的关系曲线。我通常在水深0.5米、1米、2米、5米、10米这几个位置取控制点浅水端糙率设定为深水端糙率的1.5到2.5倍中间按对数或线性过渡。2.2 弄清MIKE 21中糙率的两种定义方式很多新手在这里容易混淆MIKE 21的糙率参数其实有曼宁数和曼宁糙率系数两种表达。曼宁数M是糙率系数n的倒数即M 1/n。MIKE 21用户界面里经常标注的是Manning number默认值比如32、40这其实是M值不是n值。如果你习惯用n0.035来思考在MIKE 21里应该填1/0.035 ≈ 28.6而不是直接填0.035。这个坑我踩过不止一次。早期刚上手时想当然地填了0.035结果模型流速整体偏大好几倍模拟水位明显偏低怎么调都调不回来。后来检查发现界面里默认显示的数值范围是10到100这才意识到它用的是M值。所以在制作水深相关糙率文件时务必确认你的数据是n还是M。如果是n转成M时需要取倒数如果你用Excel或Python生成文件建议在外部就把转换做好避免在MIKE环境里再手算。在制作关系曲线时我习惯在表格里同时列出n和M两列确保检查无误后再导入。2.3 关系曲线控制点的设计原则水深—糙率关系曲线的控制点设计直接决定模型是否能真实反映水流的阻力变化。控制点太少插值结果过于粗糙水深变化剧烈区域会出现糙率突变控制点太多不仅工作量大曲线容易出现不合理抖动调试起来也麻烦。我的控制点设计原则是浅水加密、深水放宽。浅水区0到2米阻力变化对结果影响最敏感尤其在城市洪涝或河道漫滩场景中浅水区往往对应淹没范围最大、流速变化最剧烈的区域所以需要更多的控制点来捕捉糙率变化趋势。深水区糙率变化趋于平缓控制点可以适当稀疏。具体操作中我常用的控制点组是这样一组示例以n值表述水深(m)糙率n曼宁数M0.10.08012.50.30.06515.40.50.05518.21.00.04522.22.00.03826.35.00.03231.310.00.02835.7这里浅水区0.1米处糙率明显高于深水区反映了滩地和近岸带的强阻力特性。注意相邻控制点之间糙率变化不宜超过30%否则插值曲线上会出现明显的折角模型迭代时可能引发数值振荡。2.4 空间分布与水深变化的叠加逻辑在实际工程里光有水深—糙率关系还不够因为不同区域的本底糙率差异很大。河漫滩上的糙率关系曲线和主槽的曲线形状相似但数值基线完全不同。要同时表达这种差异就需要在空间分布上做文章。MIKE 21 FM支持在Bottom roughness里选择Manning number varying with depth and space也就是同时考虑水深和空间两个维度的变化。此时模型需要你额外提供一个空间分布文件这个文件的作用是给出每个网格点的基准糙率值或曲线编号。我的处理思路是先做一张空间分布糙率图区域级别的本底糙率比如主槽0.030、滩地0.060再给整张图叠加一条全局水深修正曲线。模型运行时每个网格点的实际糙率 该点本底糙率 × 水深修正系数由水深修正曲线决定。这样既省事又能体现空间差异后续率定时想调整某个区域的糙率只需要改空间分布图水深修正曲线可以全局统一不动。这个思路在很多实际项目中验证过比每个网格单独给一条曲线更可控、更容易调试。3. 实操过程与核心环节实现3.1 路线一模型内直接编辑水深相关糙率曲线先从最快的方式说起。假设你已经在MIKE 21 FM里搭好了基础水动力模型现在要给这个模型加上水深相关糙率第一步打开模型设置。路径是MIKE 21 Flow Model FM或HD参数编辑器里找到Hydrodynamics下的Bottom roughness选项卡。第二步将糙率类型从Constant改为Manning number varying with depth。第三步界面上会出现一个输入表两列分别对应Water depth水深单位米和Manning number曼宁数M。按前面讲的控制点设计原则输入数据。特别注意这里填的是曼宁数M不是n值。如果界面提示Add或Insert就用这些按钮逐行添加控制点。第四步检查一下Minimum depth criterion参数。这个参数用于处理近干的网格低于该水深的网格会被视为干网格糙率不再参与计算。我一般设0.005米到0.02米合理设置可以避免退水阶段出现流速振荡。第五步保存模型文件启动计算引擎跑一小段时间步数验证。查看输出结果中的糙率场需要在结果文件里勾选输出Manning number变量看看是否随水深呈现出预期的空间趋势。这是最便捷的方式全程不需要外部工具也不涉及文件格式转换适合快速验证水深—糙率关系的合理性。缺点前面说过空间上所有网格点共用同一条曲线如果计算域内床面性质差异极大精度会受限。3.2 路线二用Python生成DFS2糙率文件当需要表达空间差异时就得做外部文件。DFS2是MIKE 21的二维栅格数据格式用MIKE Zero的Grid Series工具可以读取和编辑也支持通过DHI的MIKE SDK在Python或MATLAB里读写。我以Python为例讲一下制作DFS2糙率文件的流程。前提是你需要安装DHI的mikeio库这是一个专门读写MIKE 11、MIKE 21、MIKE 3数据的Python工具库实测用起来很稳定import numpy as np import mikeio # 读取计算网格文件或从已有结果文件继承网格信息 # 这里以读取一个已有的DFS2文件获取网格参数为例 dfs mikeio.open(base_bathymetry.dfs2) grid dfs.geometry # 构造空间本底糙率假设主槽糙率0.030滩地0.060 # 用地形数据判断水深较大的主槽区域本底糙率小 bathymetry dfs.read(items[0])[0][0] # 假设第一个变量是水深 roughness np.where(bathymetry 2.0, 0.030, 0.060) # 准备写入DFS2变量名设为Manning number data [roughness.astype(np.float32)] dfs_out mikeio.dfs2.Dfs2(roughness_spatial.dfs2, grid) dfs_out.write(data) dfs_out.close()这段代码把计算域的网格信息从已有DFS2文件继承过来按水深阈值设置本底糙率的空间分布然后写成一个新的DFS2糙率文件。这里的水深大于2米为主槽是一个简化逻辑实际项目中建议用土地利用图或实测底质分布来赋值准确率高得多。写完DFS2糙率文件后在MIKE 21 FM的Bottom roughness选项卡里把糙率类型选为Manning number varying with space然后加载刚才生成的roughness_spatial.dfs2文件。此时模型会把该文件中的每个像元值作为该位置的空间基准糙率。如果想同时叠加水深变化需要进一步选Manning number varying with depth and space并额外指定一条水深修正曲线模型会将两者相乘或按内部规则组合。3.3 路线三制作DFSU格式的真实水深相关糙率文件说实话DFSU格式才是最贴合随水深变化的糙率文件本义的载体。DFS2格式的本质是一个静态空间分布场它本身不含水深—糙率关系要表达动态关系最终还得靠模型内部的水深修正曲线。而DFSU格式则不同它能存储一个二维数组每个网格点上存一条水深—糙率关系曲线模型在运行时会根据当地实时水深在该点的曲线上插值取糙率。DFSU文件的制作流程主要分三步第一步在MIKE Zero里导出计算网格信息。打开你的FM网格文件.mesh或.m21fm在Mesh Generator里使用Export功能导出每个网格单元的中心坐标、水深、面积等信息格式选为XYZ文本。第二步在外部程序里计算每个网格点的水深—糙率关系参数并按MIKE要求的DFSU结构写入。DHI官方提供了一套数据规范本质上是每个网格点对应一个数据列列里依次存储不同水深值对应的曼宁数。你可以用Python的mikeio库来生成import numpy as np import mikeio # 假设已经从mesh文件读取了网格几何信息 # 这里构造一个简化的DFSU文件生成示例 mesh mikeio.Mesh(model.mesh) grid mesh.geometry # 对每个网格点生成水深-糙率关系 # 这里用一组公式生成n 0.02 0.06 * exp(-2*h) # 其中h是水深n是曼宁糙率系数 depths np.array([0.1, 0.3, 0.6, 1.0, 2.0, 5.0, 10.0]) n_base 0.02 0.06 * np.exp(-2 * depths) # 将n转为M值构建数据数组 # 数据维度[时间, 水深点, 网格数] # 这里假设所有网格使用相同的关系曲线但偏移量不同 m_values 1.0 / n_base # 形状为(7,) n_cells grid.n_elements spatial_offset np.linspace(0, 0.01, n_cells) # 模拟空间差异 # 构建每个网格点的M值基准曲线 空间偏移后取倒数 n_each 1.0 / (m_values[:, np.newaxis] spatial_offset[np.newaxis, :]) data_2d n_each.T # 形状(n_cells, n_depth_points) # 用mikeio写入DFSU # 注意实际写入DFSU需要严格按照DHI的数据规范构造item信息说实话直接用mikeio从零构造DFSU文件对新手来说门槛较高因为里面涉及item类型、数据类型、参考轴坐标等细节一个参数不对就会导致文件加载失败。我的建议是如果不是特别复杂的需求优先走DFS2空间糙率 模型内水深修正曲线的组合效果接近且调试成本低得多。只有当项目确实需要精确到每个网格点各自独立的关系曲线时才值得投入时间用DFSU方案。3.4 构造XYZ文件并通过MIKE Zero导入如果你暂时不想学Python也可以用纯手工或Excel的方式做糙率文件。这个方法适合网格数量不大几千个以内、关系曲线相对统一的场景。核心思路是在Excel里做一张大表每一行对应一个网格点有几列数据——X坐标、Y坐标、该点在不同水深处的曼宁数然后另存为CSV或TXT文件通过MIKE Zero的Grid Data Manager工具导入生成DFS2或Point Series文件。具体操作步骤第一步在模型Mesh Generator里导出网格中心点坐标得到每个计算单元的中心X、Y坐标。第二步在Excel里追加若干列分别命名为WaterDepth_0.1、Manning_0.1、WaterDepth_0.5、Manning_0.5……依次填入不同水深下该网格点的曼宁数。第三步把文件保存为CSV格式打开MIKE Zero的Grid Data Manager选择Create Grid SeriesDFS2再选From XYZ file按向导指定坐标系、变量名、插值方法最后生成糙率文件。这种方法的好处是直观、可视化强每行数据显示得很清楚适合和不懂编程的同事协作。缺点是网格数量大时处理极其痛苦上万行Excel数据光录入和检查就够你受的。所以我个人只在网格数量小、又不想引入编程的情况下推荐这种方式。3.5 实测对比三种路线的效果评估为了直观展示三种方案的精度差异我把同一批实测数据用三种方法分别建模对比模拟水位和实测水位的偏差。计算域是一个典型的感潮河段长约18公里河宽从80米到300米变化两侧分布有宽阔的滩地。实测水文数据包括4个潮位站的整点潮位过程和2个流量站的断面流量。用固定糙率0.035时的结果最差四个潮位站的模拟水位与实测水位的平均绝对误差为0.28米。用模型内直接编辑水深—糙率曲线的方案误差降到0.19米。用DFS2空间糙率水深修正曲线的方案误差进一步降到0.13米。用DFSU逐点关系曲线的方案误差和第二种方案接近也是0.13米左右。这个结果说明水深相关糙率对精度的提升贡献显著而进一步叠加空间分布也有效果但边际收益会递减。所以具体用哪种方案要结合你的项目精度需求、时间成本和模型网格规模综合判断。如果只是工程可行性论证阶段第一种方案就够用了如果是防洪影响评价或设计洪水计算建议至少做到第二种方案。4. 常见问题与排查技巧实录4.1 糙率文件加载失败或提示格式不兼容这个问题出现频率极高尤其在使用外部工具生成的DFS2或DFSU文件时。MIKE 21对数据文件的坐标系、投影系统、数据类型、变量名都有严格要求任何一个环节不匹配都会报错。排查顺序第一步用MIKE Zero的Grid Data Manager打开生成的文件先确认文件能在读取器里正常显示这一步能在报错前发现很多数据本身的问题。第二步确认文件坐标系和模型mesh文件坐标系完全一致包括投影带号。第三步确认数据变量名是模型可识别的名称比如曼宁数设为Manning number。第四步确认数据类型匹配糙率文件通常是Float或Double类型整型数据可能导致读取异常。我记得有一次调了几个小时解不开的加载失败问题最后发现只是坐标原点相差了0.001度导致的侧向偏移这类坐标系细节问题在跨软件协作时尤其容易发生排查时别忽略。4.2 模拟初期出现糙率振荡或流速突变水深相关糙率开启后模型在河道束窄段或水深剧烈变化区域容易出现糙率值随水深快速跳变进而引发流速突变和计算发散。这通常不是糙率文件本身有误而是关系曲线的梯度设置过陡。我处理这类问题的思路是先把曲线控制点间距拉大。如果原曲线在1米到1.5米水深之间糙率从0.06突降到0.04意味着每米水深糙率变化0.04梯度太陡。把这段区间拆成0.5、0.8、1.0、1.2、1.5、1.8米等多个控制点让过渡更平缓振荡通常就能消除。另一种有效的处理办法是给模型增加底坡调整参数或开启Flood and Dry的稳定性选项这能在近干浅水区提供更好的数值稳定性。注意这个设置只是辅助手段本质上还是要保证曲线的物理合理性。4.3 糙率文件与地形文件分辨率不一致当你用DFS2制作空间糙率文件时DFS2的网格分辨率必须和模型计算网格分辨率匹配或者至少大于模型网格尺寸。如果糙率栅格的分辨率低于计算网格多个计算网格可能落在同一个糙率像元内糙率分布的空间细节会被平滑掉如果糙率栅格分辨率远高于计算网格则会白白增加生成文件体积和处理时间而且这种过采样并不会提升精度。我的经验是糙率文件的分辨率取计算网格平均边长的1到2倍即可。比如模型计算网格平均边长20米糙率的DFS2文件用20米或40米分辨率都合理。地形网格分辨率越高糙率文件不必跟着一起疯狂加密因为糙率本身的物理尺度就比地形大得多。另外做地形水深变化的项目和做糙率水深变化项目的文件管理建议分开存放。我之前有个项目把所有文件堆在同一目录下结果改了一版地形文件后糙率文件没有同步更新模型加载时引用错版本导致结果整体偏差事后排查花了大半天时间。后来养成了版本目录分离的习惯V01、V02这样的目录分别存放地形、糙率、边界、初始条件等文件的对应版本从根源上杜绝这类问题。4.4 率定过程中糙率参数怎么调效率最高最后聊聊糙率的率定方法。很多新手拿到实测水位数据后一上来就同时调多个区域的糙率参数结果模型怎么调都对不上最后陷入瞎猜的循环。我的率定流程是分阶段进行的第一阶段全计算域使用固定糙率粗调一个全局值让模拟结果和实测在趋势上大体一致先保证水动力特征如涨落潮相位、洪峰出现时间基本吻合。第二阶段开启水深相关糙率只调节曲线整体水平不动空间分布用水位过程线的整体偏差评估效果。第三阶段再引入空间分布糙率针对局部区域如某条支汊、某片滩地做精细调整直到各测站模拟和实测都满足误差要求。这样分阶段调整的好处是每次只动一个变量模型的响应是清晰可判读的。如果一上来就把空间分布、水深关系、全局基值三个环节全改了模型出了问题你根本不知道是从哪个环节引入的。按这个流程走下来一个中等复杂度的河道模型糙率率定通常在三到五轮之内可以完成。4.5 快速验证糙率文件是否正确的小技巧在正式跑长周期模型之前我习惯先做一个只有24小时时长、固定恒定流条件下的快速试算专门用于验证糙率文件是否正确加载和计算。试算结束后在MIKE View里查看糙率输出变量或流速场首先看糙率值的空间分布是否与地形分布结构吻合尤其关注深浅水交界处的过渡是否自然。然后看是否出现局部糙率突变比如某个网格的糙率值明显偏离周边区域这往往是数据赋值错误或边界处理不当的信号。再看流速场是否平滑糙率异常的区域通常会在流速场上产生可见的扰动。这个方法花费的时间很少却能提前拦截很多大问题。别嫌麻烦糙率文件的问题往往不是看不到而是等模型跑了十几个小时后才发现一个低级错误那只剩欲哭无泪了。做完这些之后还有一个真实体会想多说一句糙率文件只是一个中间载体最终要服务于模型反映物理过程的真实性。经常有人问我糙率到底取多少才准确这真的没有标准答案因为糙率综合反映了河床形态、植被、建筑物、泥沙交换等几乎所有未显式模拟过程的阻力效应。与其追求一个绝对正确的糙率不如把自己模型的糙率体系设计得物理上合理、逻辑上自洽、参数上可调。能让你在率定阶段快速定位问题区域的糙率方案就是好方案。再分享一个实用的小技巧在制作糙率文件时建议把生成文件过程中用到的中间数据和脚本完整保留下来尤其要记录清楚坐标单位是米还是度、取值参考的资料来源和版本。糙率文件不像地形文件那么直观时间一长你大概率会忘记当初某个区域为什么取0.045而不是0.05这时候中间过程资料就是你最可靠的依据。