MIKE11培训讲解教程.docx

1、MIKE11培训讲解教程MIKE11培训讲解教程 MIKE 11 HD、SO、AD、Ecolab完全实例教程 通过建立东江水动力、水质模型，熟悉MIKE 11基本功能和应用技巧。练习中将要建立的模型为DHI咨询项目实例（世界银行贷款项目），练习所使用的数据完全为该项目真实数据。 让我们从最原始数据开始！ 1 MIKE 11 HD MIKE 11 水动力模型 MIKE 11 HD建模需要以下各类数据或信息 流域描述 河网形状，可以是GIS数值地图或流域纸图； 最好能包含水工建筑物和水文测站的位置 河道和滩区地形 河床断面，间距视研究目标有所不同，但原则上应能反映沿程断面的变化 滩区地形资料（有时

2、有滩区的水位蓄水量关系曲线也行），如果要模拟滩区行洪的话 模型边界处水文测量数据 边界最好设在有实测水文测量数据处，如果实在不可能那么就不得不估算边界条件了！ 率定验证用水文测量数据 率定验证的数据越多，模型就越可靠，但工作量也会越大。 水工建筑物设计参数及调度运行规则 所有水工建筑物都有其调度运行规则，问题是1) 能否收集到有关资料？2) 实际是否按照此规则运行了？建筑物处的实测水位或流量资料可以反映实际运行情况。 MIKE 11将上述信息分别存放在以下文件内: 河网文件 (.nwk11) 断面数据 (.xns11) 边界条件 (.bnd11) 模型参数文件 (.hd11) (时间序列文件)

3、 (.dfs0) 下面开始逐个建立这些文件。 河网文件 河网文件或许是MIKE 11所有文件中最复杂、又是最重要的一个文件了。 收集到的原始数据为ArcView河网文件（river.shp）：包含了河道位置、河道长度等信息。 建立河网文件步骤： 第一步：从ArcView河网文件提取河网背景图 用ArcView打开river.shp文件 ? 缩放图形区域至完整包含所有将要模拟的河道? 记下此区域左下角和右上角的UTM坐标，比如分别为54，*-*.05)和72，*-*.16) ? 输出为bmp格式图（File ? Export.） 注：如果是纸图，同样可以将其扫描后生成bmp文件。用尺在地图上丈量

4、左下角和右上角的坐标，可以设左下角坐标为（0，0）。如果模拟范围不是很大（比如整个长江流域），那么这种对地图的平面假设所产生的误差可以忽略不计。 生成的bmp图应大致同如下图1： 图1 第二步：确定所要模拟的各河段名及长度 一般水利部门有河段长度信息，从数值地图、AutoCAD图、甚至纸质地图上也能量出各河段长度。 东江模型各河段名、位置及长度见图2和表1 表1 河名 长度，米 河名 长度，米 Beihai * Ma 3747 Chifen 9225 Ou * Dafen 2650 Qiandong 2640 Dao * Qianxi 3810 Dongguan * She 6745 Dong

5、jiang * Shizijiao * Dongjiangnan 9950 Shu 3017 Dongjiangupper * Taiyangdong 3115 Dongnangu 7169 Taiyangxi 4160 Dongnanhou * Wu 2926 Dongnanhou1 1296 Xiancunchong 6036 Gu 915 Xibian 3290 Heng 4830 Zengjiang * Henghai 2550 Zhongtang * Hongguo * Zhongtang1 824 Huang 8245 Zengjiang Zengjiang Dongjiang D

6、ongjiang upper Dongjiangnan Dongnangu Huang She Dongnanhou Shizijiao Xiancunchong Ma Gu Heng Hongguo Ou Dao Xibian Qianxi Qiandong Taiyangxi Taiyangdong Dongguan Shu Wu Dongnanhou1 Hongguo Chifen Beihai Dafen Zhongtang1 Huang Zhongtang1 图2 图2 第三步：建立东江模型的河网文件 1) 引入河网背景图 打开MIKE ZERO ，File ? New ? MIKE

7、 11 ? River Network（见图3） ? OK，弹出一个新窗口（图4 ）? 输入河网模型区域的范围（左下角和右上角坐标）? 输入刚才记下的背景图左下角和右上角坐标54，*-*.05)和72，*-*.16) ? OK, 出现河网文件视图（模拟区域暂时空白） ? 河网文件菜单Layers ?Add/Remove. ? 点击添加项目键 ? 点击浏览按钮, 引入刚才生成的bmp背景图 ?回到河网文件视图，Layers ?Properties. ? 将图像坐标Image Coordinates修改为正确值54，*-*.05)和72，*-*.16)。至此背景图引入工作完成，河网视图应与图5大致

8、相同。 图3 图4 图5 2) 输入各河段信息 使用河网文件编辑器内的工具条定义各河段 例如要在背景图上绘制各河段时按，在河段的起始点开始鼠标左击背景图上相应位置，在河段结束点鼠标左键双击，完成该河段的绘制。 可以放大背景图以准确绘制河段走向。准确的河段走向有助于确定各河段连接关系、确定该河段上水工建筑物或水文测站的里程，所以是很重要的工作。 请参见“帮助”或按F1查看其它工具按钮的功能。 以下是定义河段信息的步骤。注意：MIKE 11目前暂时还不能使用恢复键，所以在操作过程中请随时保存，一旦操作失误，您只能删除重做或不保存退出，重新进入河网文件编辑器！ 在背景图上绘制某条河段 打开河网文件编

9、辑器的表格视窗（View ? Tabular View.） 在左侧列表区展开Network，选Branch。在右侧的河段信息内容中可以发现刚才所绘制河段的信息，MIKE 11已自动为该河段命名，并确定了其长度。将该河段名改为实际名称。注意避免使用中文河段名称！河段长度一般不会与实际长度一致，将在以下d) 步骤中修改。 在左侧列表区选Points。在右侧的河段点信息内容中可以发现刚才在绘制河段过程中每次点击点的坐标位置（MIKE 11已自动测出！），您会发现河段名已经是实际名称。注意：这些点并不是模型的计算点，与计算没有任何关系！里程类型Chainage Tpye列上，我们将该河段的起始点从Sy

10、stem Defined改为User Defined，将接下来一列里相应行的里程数Chainage改为0（程序缺省值为0）；将该河段的结束点从System Defined改为User Defined，将接下来一列里相应行的里程数Chainage改为实际的河段长度。这样我们已经将程序测出的河段长度改成了实际长度。注意：任何河段的起始里程可以是任何数值，正、零或负数，整数或小数。选值的原则是与当地水利部门采用的桩号值一致，这样将来若要与他们讨论问题时会比较方便，否则就取0。必须要满足的是： 河段长度 结束点里程数 起始点历程数 重复以上步骤a) 至d)，完成所有河段信息的输入。 在河网文件编辑器的

11、图像视窗内用工具按钮连接各河段。注意：有多条河段相连时必须所有河段同时连向某条河段。但连接方向（谁连向谁）对计算结果没有任何影响。 在河网文件编辑器的表格视窗内Network?Branch，在右侧的河段信息总览表第二列是地形标识TopoID信息。这是文件编辑器将来从断面文件编辑器内读取与该河段相对应的断面数据信息、参与模型计算的唯一信息通讯通道，一定要与断面文件编辑器内相应的TopoID一致。TopoID可以是数值，也可以是文字符，一般用断面测量年份，如Topo2002，这样便于提醒自己目前模型采用的断面数据来自哪年实测数据，或Artifical，表明该断面数据并不是真实数据。参见有关断面文件

12、编辑器内容。 在河网文件编辑器的图像视窗内Settings?Network.或Font.你可以对图像的外观做随意修改，选择自己想显示的信息。 至此河网基本信息输入完毕。 练习：以东江模型的河网为例，熟悉河网编辑器的操作及编辑和显示功能。 改变里程值，观察MIKE 11如何自动重算各点里程。 改变显示的颜色和符号。 关闭河网编辑器（不要保存）。创建一个新的河网文件并输入背景图。 完成的东江模型河网应与图6基本一致， 图6 1.2 断面文件编辑器 收集到的原始数据为文本格式的所有河道各断面x,z数据，断面的位置为AutoCAD文件。为节省时间，已转换成符合MIKE 11输入要求格式的文本 文件。所

13、要求的格式见本文件附录一MIKE 11断面输入格式说明，用*、BASIC、EXCEL VBA等自编小程序，很容易将原始数据格式转换成符合MIKE 11输入要求格式的文本 文件。开始以下步骤时假定已有转换好的文本文件东江模型断面.txt。 第一步：打开断面文件编辑器 MIKE 11 ? New ? Cross Sections （见图7） 图7 第二步：在断面文件编辑器中引入断面数据 File ? Import ? Import Raw Data Recompute，找到东江模型断面.txt文件保存路径? OK, 引入成功后在断面文件编辑器视窗左侧出现断面列表，右侧图像视窗显示一个或多个断面形状

14、。 第三步：在图像视窗内查看各个断面形状，从直观上判断断面数据是否合理。 第四步：保存文件，断面文件生成完毕。 单个断面的插入、修改等操作法方法学习结合练习进行；断面文件编辑器中一些参数设定方法详见PPT。 练习：熟悉断面文件编辑器的显示和编辑功能。 尝试修改、复制和插入新的断面 从图像修改断面数据 观察处理后数据表格 输出断面数据至一文本文件。关闭断面文件（不要保存），创建一个新的断面文件，输入断面文本文件。 1.3 边界文件编辑器和时间序列编辑器 所有外部边界条件和内部边界条件都在边界文件编辑器里设置。所谓外部边界就是模型中那些不与其它河段相连的河段端点（即自由端点）物质流出此处即意味着流

15、出模型区域，流入也必然是从模型外部流入，这些地方必须给定某种水文条件（如流量、水位值），否则模型无法计算。所谓内部边界是指从模型内部河段某点或某段河长流入或流出模拟河段的地方，典型的例子包括降雨径流的入流、工厂排污、灌溉取水，内部边界条件应根据实际情况设定，设不设通常不会影响模型的运行，但会影响到模拟结果的可靠性。 MIKE 11用于边界条件的数据需要转换成时间序列文件（带后缀dfs0）方能被引用。 东江工程收集到的外部边界的原始资料包括上游边界博罗和麒麟站、下游大盛、麻涌、泗盛围、漳澎的水位流量，保存在水文要素.xls文件内。 第一步：将EXCEL文件转换成dfs0格式的时间序列文件，以下游

16、边界水位时间序列文件为例： 打开水文要素.xls文件，观察数据内容：所有下游边界的水位流量数据时间长度都为2001年2月7日14:00至2月16日10:00，时间间隔为1小时。上游博罗和麒麟站水位流量为日均值，时间范围也是7日至16日。所有站的数据都是时间等间隔的。并且由数据可见，上游边界处的水位已基本不受下游潮汐的影响了。 MIKE 11 ? New ? MIKE Zero ? Time Series ? Blank Timeseries ? OK，您会看到弹出如图8所示的dfs0文件窗口。 图8 ? 时间轴类型Axis Type选等时间间隔； ? 开始时间输入2001年2月7日14:00:

17、00； ?时间步长输入1小时； ?时间步长数输入213； ?在Item Information区内Name栏输入大盛水位（可以用中文名称！），类型选Water Level，单位为米（缺省值） ?按Append三次，新添三行，同样分别填入麻涌、泗盛围、漳澎的有关信息，按右上角OK按钮 ?在出现的时间序列文件视窗内检查最后一列的时间是否为2001年2月16日10:00:00。如果不是，在左侧的图像视窗内鼠标右击，选属性properties., 回到刚才窗口进行修改。如果是，在打开的Excel文件内选择这些站的水位数据，复制（或Ctrl-C）后回到时间序列文件视窗，在右侧的表格视窗内选相应的列，按粘

18、贴（或Ctrl-V），如果EXCEL内各站的排列顺序与dfs0的一样，可以一次同时复制和粘贴所有这些站。 保存文件，完成dfs0时间序列文件制作。记住文件制作过程中要随时保存！ 重复以上步骤，完成所有边界文件的制作，对于有同样时间轴的数据，即使是不同数据类型（如水位、流量、浓度）也可以放在同一个时间序列文件内。 所有dfs0文件完成后，在边界文件编辑器里输入边界条件： MIKE 11 ? New ? MIKE 11 ? Boundary Condition，? OK， 出现边界文件窗口，如图9。 图9 以输入上游博罗站的流量边界条件为例： ? 在边界描述Boundary Description

19、栏选Open；边界类型Boundary Type栏选Inflow；河段名为DongjiangUpper；填入正确的里程数。 ? 视窗中间区的边界计算内容暂时只选Include HD Calculation ? 在接下来的水文边界信息区内TS类型选TS File，按找到存放刚才生成的时间序列文件的路径，并选择正确的项目（在dfs0文件中可能有许多项目，如刚才生成的一个dfs0文件内有四个水位项）。右侧的 是用于打开对应的时间序列文件，而最后一栏的TS Info是提醒你刚才选中了dfs0文件中的哪一项。 ? 光标返回到最上面区域，按TAB键，添加所有其它边界条件。总共有6个外部边界条件。 现在我们

20、建立两个内部边界条件： 在DongjiangUpper的上游*里程（各人根据自己的河网定义可能有所差异，博罗以下约23km）处，东深取水口，取水量为70m3/s。 在DongjiangUpper的上游博罗以下1km处，某工厂排污，排水量为0.36m3/s。 光标返回到最上面区域，按TAB键两次，再添加两个边界条件。 1）添加取水口边界 ? 在边界描述Boundary Description栏选Point Source；边界类型Boundary Type栏选Inflow；河段名为DongjiangUpper；填入正确的里程数。 ? 视窗中间区的边界计算内容暂时只选Include HD Calcu

21、lation ? 在相应的水文边界信息区内TS类型选Constant，数值输入-70（取水！） 2）添加工厂排污边界 ? 在边界描述Boundary Description栏选Point Source；边界类型Boundary Type栏选Inflow；河段名为DongjiangUpper；填入正确的里程数。 ? 视窗中间区的边界计算内容暂时只选Include HD Calculation ? 在相应的水文边界信息区内TS类型选Constant，数值输入0.36 完成边界条件输入，保存文件。注意文件生成过程中随时保存文件。 练习：使用时间序列编辑器 - 在图像窗口修改数据，注意在表格视窗内数据

22、的变化。 - 改变时间序列编辑器的显示属性。 - 用工具栏中的“calculator”功能改变时间序列数据。 - 在Excel和时间序列编辑器之间进行数据交换。 1.4 HD参数文件编辑器 参数文件主要是定义模拟的初始条件和河床糙率。里面尽管有许多菜单但，大部分内容不必去触动。 第一步：设定初始条件 MIKE 11 ? New ? HD Parameters ? OK, 弹出参数文件窗口，如图10。 图10 进入初始条件Initial菜单，添加初始水位和流量。初始条件设定的一个很重要目的是让模型平稳启动，所以原则上初始水位和流量的设定应尽可能与模拟开始时刻的河网水动力条件一致。实践中，初始流量

23、往往可以给个接近于0的值，而初始水位的设定必须不能高于或低于河床，否则可能导致模型不能顺利起算。山区性河道往往坡降很大，初始水位有时很难设定，往往须用其它方法解决这一问题，详见有关模拟文件编辑器内容。 对于东江模型，在察看个河段的河床地形后，我们可以简单设置： 全域值：水位0米；流量1m3/s；局部值：河段DongjiangUpper上游初始水位4米；其下游初始水位0米。 第二步：设定河床造率 进入河床造率Bed Resist菜单，设定河床造率。河床造率是率定参数，应根据对模拟河道的认识及模型计算结果确定。通常可以从n=0.03开始率定。各个河段设定不同造率值的方法详见PPT文件。 1.5模拟

24、文件编辑器 我们几乎已经完成了所有基本的模型设置，可以开始计算了！ 模拟文件编辑器的作用是集成刚才所生成的所有文件的，让它们成为一个整体；同时定义模拟时间步长、结果输出文件名等。 第一步：打开模拟文件编辑器 MIKE 11 ? New ? MIKE 11 ? Simulation ? OK, 弹出模拟文件窗口，如图11。 图11 第二步：选择模型类型。暂时选择HD模块。 第二步：进入输入Input菜单 高亮框表示对于HD模拟需要这些文件。按按钮引入刚才生成的所有文件：河网文件、断面文件、边界文件、HD参数文件。可以随时点击编辑这些文件。 如果我们从这里打开刚才做好的河网文件，可以发现许多刚才被

25、禁用的功能键都已经被激活，可以使用了。这是因为通过模拟文件编辑器已经把所有文件链接起来，可以通过河网文件访问其它文件了，如断面文件、边界文件和参数文件。 第三步：进入模拟simulation菜单 我们选择固定时间步长5分钟；模拟开始时刻为2001年2月7日14:00:00；结束时刻2001年2月16日10:00:00。 时间步长的确定经常要通过反复试算调整，与河床地形与边界条件密切相关，并且原则上要满足克朗数不大于10。 对于如山区河流等初始条件不易合理设定的情形，缩小时间步长是一个行之有效的方法。对于模拟时段内短时间有大量流量进出、而其它时期比较平稳的情形，选择可变时间步长比较合适，可大大缩

26、短计算耗时。 接下来我们选择初始条件设定Initial Conditions。还记得我们曾经在HD参数文件里也设定过初始条件？如果现在我们在这里选择参数文件Parameter File，那么刚才的设置就会起作用；如果现在我们选择稳态Steady State，那么MIKE 11就会“忘记”参数文件内有关初始条件的设定，而根据边界点上给出的水位流量数据（从边界文件中调用），利用稳态假设计算各计算节点上的初始水位流量；如果我们选择Steady + Parameter稳态参数，那么在参数文件中做过特别设定的河段节点上模型用这些设定值作为初始条件，其它点用稳态假定计算，如同选择了稳态方法；如果选择了热启

27、动Hotstart，那意味着要用以前的模拟结果作为当前模拟的初始条件。还是用山区河流为例，其初始条件可能比较难设，但模型一旦运行一段时间后就很可能会比较稳定，即可以增加时间步长。用热启动就能解决这一矛盾：先用非常小的时间步长计算（如0.1秒），当计算稳定后（如计算一天）停止计算；重新计算模型，用大时间步长（如5分钟），用热启动模式，将刚才模拟结束时刻的计算结果作为当前模拟的初始条件。 对于东江模型，初始条件不很敏感，我们选择稳态或参数文件都可以。 第四步：定义输出结果文件名和保存频率。 假定计算时间步长为5分钟，但我们步需要这么密的计算值，或许半小时一个结果已经足够，那么我们可以定义保存频率为

28、6，即计算6个时间步保存一次结果。这样可以减小结果文件大小。 第五步：准备计算 进入Start菜单，准备开始计算。如果验证状态Validation Status框内都是绿灯，那么您就可以按Start键开始计算了；如果有红灯，那么在下面的验证Validation Message框内就会出现相应的出错信息，提醒您进行修改。这只是初步检查，只能检出一些明显的模型设置错误。 第六步：模型运行 如果出现如图12的运行进度框，那么恭喜您！模型设置成功，正在运行！ 有时会出现警告信息Warning Message，这是MIKE 11认为模型设置可能存在一些小问题，但这些问题还不至于影响到模型的运行，因此提醒

29、您去检查一下。若您认为没有问题就可以要求MIKE 11继续运行下去。 图12 2 MIKE View MIKE View用于演示和提取包括MIKE 11在内的一些DHI软件计算结果。对于协助率定和报告制作非常有用。 2.1 装载结果文件 开始?程序?DHI软件?MIKE 11?MIKE View打开MIKE View。在MIKE View窗口内File ?Open，文件类型选MIKE 11结果文件类型（*.res11)，出现如图13对话框，点击OK后进入MIKE View主窗口。 图13 2.2 平面动态演示（以流量为例） 在平面视图上点击鼠标右键 ? Option. ? 弹出对话框，如图14

30、 ? Draw Branches As选最后一列Water Level所在列 ? 在该列的下拉菜单中选Discharge ? 进入Symbols and Fonts菜单 ? 设置河段线条粗细及是否绘制流向箭头 ? 点击OK退出 ?在主窗口上方工具条内按播放键开始动态演示。 显示流量大小的各种颜色通过调节调色板实现（可以在主窗口上方工具条内按调色板工具显示调色板）： 在调色板内点击鼠标右键 ? 选择各个菜单调节分级、色彩，还可以保存供将来调用。 图14 2.3 剖面动态演示 在主菜单上方工具条内按按钮，当光标回到平面视图内，接近河段时，光标变形，表示已检测到河段，可以对之选择剖面显示。沿着你想显

31、示的各河段前进，途中如果你选错了河段，按退格键*CE返回一次操作，按退出键ESC放弃所有选择。当想结束选择时按住Ctrl键点击鼠标左键，在回答有关提示后显示剖面的水位或流量，并可对之进行动态演示。 2.4 模拟与实测结果的比较 最常用的率定手段 I) 绘制计算结果图。如果知道哪些河段节点有实测水文资料，可以用于率定，那么有两种方法绘制这些点的结果： a) 知道测站在地图上位置，但不清楚具体里程数。这时您可以按工具条内的按钮，弹出数据类型选择Data Type Selection对话框，选择您想比较的数据类型（水位或流量） 点击OK，当光标回到平面视图内后，靠近该河段节点时，光标变形，表示已检测到该节点，左击鼠标，就会出现该节点