mike21-操作实例---矩形网格

MIKE 21 HD操作简要说明 1 一一 数据文件格式数据文件格式 1 land xyz 陆地边界文件 格式说明 纬度 tab 经度 tab tab 陆地高程值 12 30307 56 12769 1 10 12 30895 56 12909 1 10 12 30895 56 12784 1 10 12 30307 56 12769 0 10 说明 第三列中 0 表示陆地边界的结束点 该点与陆地起始点重叠 每两列之间用 TAB 键隔开 不要用空格 2 water xyz 水深文件 格式说明 12 16805 56 17699 20 00 3 12 16665 56 17715 20 00 3 12 16609 56 17762 20 00 3 12 16637 56 17808 20 00 3 说明 第一 二列分别为纬度 经度 第三列为该点处水深 第四列暂时不详 其他 说明同陆地文件 对于上述文件格式也可用直角坐标 大地坐标 表示 但是在对于上述文件格式也可用直角坐标 大地坐标 表示 但是在 Define working Area 对话框中要选定对话框中要选定 Decimal Deg 项 如下图所示项 如下图所示 3 WL1 水位边界条件文件 格式说明 Water level recordings from St 1 time Elevation 1993 12 02 00 00 00 0 33 选择该项 2 1993 12 02 00 30 00 0 362 1993 12 02 01 00 00 0 3965 说明 第一列日期 第二列时间 第三列水位 第一列和第二列之间用空格隔开 第第一列和第二列之间用空格隔开 第 二列和第三列之间用二列和第三列之间用 TAB 隔开 隔开 前两行为文件说明语句应保留 4 windcastrup wind 风边界条件文件 格式说明 Wind data TimeSpeedDirection Unit 1000022000010000324010 1993 12 02 00 00 009 294184 26 1993 12 02 00 30 0010 066186 689 1993 12 02 01 00 0010 655189 167 说明 第一列为日期 第二列为时间 第三列为风速 第四列为风向 与北向顺时针 夹角 第一 二行为文件说明语句应保留 第三行暂时不详 第一列和第二列之间用空格 隔开 第二行和第三行及第四行之间用 TAB 键隔开 5 WaterlevelNdrRoese 水位验证文件 格式说明 Ndr Roese Water Level m TimeNrd Roese Water Level m Unit 10000010000 1993 12 02 00 00 00 0 344375 1993 12 02 00 15 00 0 349063 1993 12 02 00 30 00 0 348588 1993 12 02 00 45 00 0 397647 说明 该文件为 Ndr Roese 站实测水位 第一 二列之间用空格隔开 第二 三列 之间用 TAB 隔开 第三列为水位 6 WindKastrup 风过程文件 格式说明 风过程文件 格式说明 二二 以实例演示以实例演示 MIKE21 的水动力模型及输移扩散模型的操作步骤的水动力模型及输移扩散模型的操作步骤 1 网格数据网格数据 1 命令操作流程 命令操作流程 水深网格文件设置菜单流程 打开 mike21 file new Bathymetries 打开 Define Working Area 对话框后设置 选择菜单 Work Area 选择 Background Management 点击 import 按钮 选择文件类型为 xyz 找到陆地边界文件 land 再选择 convert form 中的 LONG LAT 后打开陆地边界文件 在点击 import 按钮 按照陆地边界文件的操 作打开水深数据文件 后 ok 在选择 Work Area 中的 Bathymetry Management 后 出一个对话框 点击 New 出现计算区域设置对话框 设置后 OK 此时图形已经框出 计算区域 在点击工具栏中按钮选择要插值的水深数据范围 选定后在点击一次 3 再打开 Work Area 中的 Bathymetry Management 点击插值按钮 interpolate 设置后 OK 保存已经设置好的计算区域网格文件 例如文件名 Bat2 文件 File open 打开已经生成的 Bat2 0 文件 可以应用工具栏中的工具及 Tools 菜单中的 Set Value 编辑水深及陆地值 存盘 至此已经完成网格数据文件 的编辑 嵌套网格数据文件的设置流程 第一步 建立子网格 打开 mike21 file new 选定 mike zero 中的 Mike zero Toolbox 展开 Transformation 选定 Rotate Grid 此时弹出一个对话框 setup Name 设置子网 格区域的名称 下一步 出现 select input file 选择要嵌入子网格的网格文件 例如 Bat2 0 文件 下一步 出现 subseries selection 在此可以设置子网格需要计算的时间 段 也可跳过 下一步 select item 该项暂时不详 跳过 下一步 出现 Grid Size 该项用于设置子网格的尺寸 第一项 设置子网格的结点数 结点数 3 大网格数 1 且 子区域边界要和大区域中的网格线重叠 第二项 子网格的网格宽度 取上一级网格宽度 的 1 3 第三项为子区域的左下角在大网格区域中的位置 最后一项为子网格区域在大网格 中的转角 与直角坐标 X 向的夹角 下一步 因子网格区域是否包括陆地和水域而提供 的选项 本例中包含陆地和水域故选择 date is Bathymetry 项 下一步 output specification 输出子网格的文件名称 以 Bat2 0zi 为例 下一步 点击 Execute 按钮 完成 第二步 子网格边界处理 打开 mike21 file new 选定 mike 21 中的 Mike 21 Toolbox 展开 Hydrodynamics 选定 Boder Adjustment 设置名称 下一步 出现 Input File 界面 第一项 输入大网格文件 第二项 输入小网格文件 最后两项 X coordina1和 Y coordina1 为子区域左下角在大网格中的位置 下一步 出现 Select Output File 界 面 第一项 输入处理后的大网格文件名称 第二项 输出处理后的小网格文件名称 下一步 点击 Execute 按钮 完成 经过上述两个步骤就完成了对嵌套网格的设置 对于嵌套两个同级网格 如下图所示 的操作说明 将 大 网格和 小 1 网格按上述两步处理后 得到 大 处理 和 小 1 处理 后在将 大 处理 和 小 2 网格按上述两步处理后 得到 大 处理 1 和 小 2 处理 则在水动力模型计算时 输入 大 处理 1 小 1 处理 小 2 处理 三个网格文件 如果三重嵌套 例如 小 2 有嵌套在了 小 1 中 则要两 两发生关系 既 大 和 小 1 处理后的 小 1 处理 在和 小 2 处理 小 2 处理 在和 大 处理 进行处 理 多重嵌套以此类推 大 小 1 小 2 4 2 图形流程示意 图形流程示意 第一步 见图 1 图 1 第二步 见图 2 图 2 第三步 加入背景文件 见图 3 图 4 选定 UTM 带号 Utmzone 经度 6 31 的整数部分 定义区域左下 角在大地坐标 系下的位置 显示区域的宽显示区域的高 5 图 3 图 4 将陆地边界和水深数据加入之后图形显示如图 5 所示 图 5 第四步 选定计算区域 见图 6 图 7 在显示的区域内添加 陆地及水深背景 2 文件类型为 xyz 1 输入 3 陆地和水域文件坐 标点均为经纬度 故 要选择将坐标由经纬 度转化为直角坐标选 项 6 图 6 图 7 经过上述设定后既圈定了计算区域 见图 8 图 8 第五步 在圈定计算区域之后 要选定要插值的数据范围 应用工具栏中的按忸 在 图形中框取 选定后再点击一次按忸 之后在调出如图 9 的对话框 选择该项 编辑计 算区域 1 新建一个计 算区域 2 指定计算区域左下角 在显示区域中的坐标 5 指定计算区域的 转角 与北向夹角 3 指定计算网格的宽 度 4 指定计算 区域在 X Y 方向 的结点数 6 陆地区域高程 值 要和陆地 边界文件一致 该框圈定区域 既为计算区域 7 图 9 插值结束之后会出现消息框 该框显示插值所用的时间信息 待全部 ok 之后 存盘 文件名为 bat1 此时系统自动生成一个文件名为 bat1 0 的文件 该文件既为计算用 的网格文件 然后打开该网格文件 如图 10 所示 图 10 第六步 网格数据调整 可用工具栏中如图 11 中所示的按钮选定要编辑的网格 然后选择 菜单 Tool 下的 Set value 见图 12 图 13 调整后如图 14 所示 保存 图 11 图 12 图 13 1 选择该项 对 计算区域的数 据进行插值 2 指定插值时寻 找结点周围数 据的个数 3 选择插 值方法 需要设陆地 高程值的区 域 错误区域 在此设置选定网 格的值 8 图 14 处理后的网格数据图 至此计算区域的网格数据文件已经整理完毕 如果要在本计算区域内内嵌网格 操作步骤如下 第一步 建立子网格 图 15 图 16 图 18 图 17 选择该项选择该项 设置子网格名称 选择被嵌入的大网格文 件 9 图 19 图 20 该项暂不详 图 21 图 22 图 24 图 23 至此已经完成了小网格的编辑 小网格的文件名为 Bat1 0 qiantao 图 25 所示为嵌套 后的小网格 设置嵌套的小网格中 计算的起始和结束时 间 默认时间同上一 级网格的计算时间 指定 X Y 向的小网格结点数 结点数 3 大网格数 1 小网格宽度 网格宽度必 须为上一级网 格宽度的 1 3 指定小网格 左下角在大 网格中的位 置 嵌套网格在大网格中 的转角 与 X 方向 选择该项 指 明小网格中包 括陆地和水深 输出编辑好的 小网格文件 结束之前先 运行 点击 该键 10 图 25 嵌套好的计算区域 第二步 子网格的边界处理 图 26 图 27 选择该项 选择该项 输入大网格文 件名 输入小网格 文件名 输入小网格左 下角在大网格 中的位置 11 图 29 图 28 图 30 图 31 至此已经完成了嵌套网格的准备工作 也产生了两个用于模型计算的网格数据文件 在本例中分别为 大网格文件 Bat1 0 chuli 小网格文件 Bat1 0 qiantao chuli 2 边界条件的设置边界条件的设置 1 命令操作流程 命令操作流程 第一步 设置时间序列文件 dfs0 打开 mike 21 file new 选择 mike zero 中的 Time series 时间序列 选择 From Ascii File 出现 Time series Eidetor import from Ascii 对话框 在此对话框中 选择 File Name 本例为 WL1 txt 将 Time discription 项设置成 Equidistant Calendar Axis 等时间序列轴 ok 出现水位边界条件 WL1 TXT 的过程线 此时鼠 标右键单击此图 选择 properties 出现 File properties 对话框 将 Axis Type 设置 成 Equidistant calendar Axis 将 Item information 中 Type 项中的 undefined 设 置成 water level 水位 ok 保存文件 本例为 Ts1 dfs0 重复上述的操作 直至将上下两个边界的 4 个点时间序列文件处理完毕 验证文件也按上述操作过程设置 第二步 建立边界条件 打开 mike 21 file new 选择 profile series 断面序列 选择 blank Time series 空白时间序列 此时弹出 File properties 对话框 设置 title 上边界条件 名称 本例为 shangbianjie 设置 Axis Type 为 Equidistant calendar Axis 设置 Time step 为 1800s 因为边界条件中的时间间隔为 30 分钟 设置 No of points 为 13 本例中 上 边界的网格数为 13 个 设置 Grid step 网格步长 为 900 将 item information 中的 Type 项由 undefined 该为 water level ok 此时出现一个空白的断面边界条件编辑器 打开北边 界的左侧的点时间文件 TS1 dfs0 并将其复制拷贝到空白单元格的第一列 打开右侧的点 时间序列文件 TS2 dfs0 并将其拷贝复制到空白单元格的最后一列 选择 Tools 菜单中 的 interpolation 进行插值 然后存盘 本例文件名为 shangbianjie 至此已经完成了上边界的准备 重复上述操作过程设置下边界条件 本例中下边界条 件文件名为 xiabianjie 输出调整后 的大网格文 件名 输出调整后 的小网格文 件名 设定都结束 后 运行完 成 12 2 图形流程示意 图形流程示意 第一步 生成 dfs0 文件 图 32 图 33 图 34 该项为时间 序列项 指定从文件 建立 选定要处理的文件 选定该项 设置成等时 间序列轴 13 图 35 经过上述设置后 将文件存盘 重复设置直至将所有的边界处的点文件设置完毕 设 置后的文件后缀名为 dfs0 第二步 建立边界条件文件 图 36 图 37 在图形中单击 鼠标的右键 在选定此项 将此项设置 成 water level 选定该项 断面序列 选定该项 建立一个空白文挡 14 图 38 图 39 建立好空白文档后 打开左边界点的 dfs0 文件 并将其复制到空白文档的第一列 打 开右边界点的 dfs0 文件 并将其复制到空白文档的最后一列 然后按图 40 图 41 所示的 操作进行插值 图 40 设置上边界文件名称 设置成等时间序列轴 设置起始时间 该时间不能早 于以有的边界 水位数据的对 应时间 时间步长 要等于数据 文件中的时 间间隔 时间步数 步长 步 数 总时间 长度 该项设置成 water level 上边界的网 格数 网格步长 该界面为已经建 立好的空白文档 15 图 41 设置好后存盘 本例的文件名为 shangbianjie 然后重复上述的操作设置下边界条件 文件名为 xiabianjie 至此 边界条件已经建立完毕 3 水动力模型的设置和计算水动力模型的设置和计算 图形流程示意图形流程示意 图 42 选定该项 其他各项暂不详 选择该项 16 图 43 图 44 水动力模型 水动力及输移扩 散模型 水动力及泥沙模型 指定计算区域的个数 本 例中包括一个嵌套网格 股设置 2 点击此键 添加数 据网格文件 1 为 大网格 2 为小网 格 指定小网格 在大网格中的 位置 该数指定了这 一网格嵌套在 了第一重网格 内 17 图 45 图 46 时间步长的个数 该数 时间步长 总的时间长度 时间步长 模型计算的起始时间 通过对该参数的设置 将该 参数控制在 7 以下 模型中上 下 4 个边界点的位 置 18 图 47 图 48 设置源 汇 的个数 质量平衡 19 图 49 图 50 设置干 湿边界 值 设置模型初始值 20 图 51 图 52 上边界的两个端点位置 选择该类型 意义不详 输入上边界条件文件 设置源 汇的沉降和蒸 发值 21 图 53 图 54 涡流粘滞系数 曼宁系数 22 图 55 图 56 波辐射 风条件设置 23 图 57 图 58 由于有一个嵌套模型 故设置输 出 2 个结果文件 结果文件的大小 该对话框对计算结 果的输出进行控制 该项设置多少时间步长存贮一个 结果 本例中设置了 6 个时间步 长 既 6 100s 输出一个结果 设置输出的运行结果 文件的名称及位置 24 图 59 图 60 设置结束后 选 定该选项运行模 型 运行界面 25 至此 模型设置及运行全部结束 4 计算数据提取 验证 计算数据提取 验证 1 命令操作流程 命令操作流程 第一步 提取数据 打开 mike 21 file new 选择 mike zero 中的 mike zero toolbox 展开 Extraction 选择 Timeseries form 2D files 从二维模型文件中提取时间序列 出现 setup name 设置提取操作的名称 下一步 出现 input date 指定数据的来源 下一步 出现 subseries specification 设置提取数据的时间范围 下一步 出现 select items 选择要提取 的内容 下一步 出现 point select 设置提取点的数目和位置 下一步 出现 output file selection 设置要输出的数据文件的名称及 标题 下一步 运行 Execute 完成 第二步 对比计算结果及验证资料 file new 选择 mike zero 中的 plot composer 图形设计 ok 出现一个空白文 档 选择菜单中的 plot 项中 insert new plot object 插入新的图形文件 选择 time series plot 时间序列图形 ok 点击按钮 new item 出现打开对话框 选定要 添加的文件后 在该对话框下部的 select item 中设置要绘制的内容 ok 确定 对于键可重复操作 添加多个图形文件的输出 也可按如下操作 鼠标右键单击图形 出现菜单后 选择 properties 选项 在按上述方法添加 此时图形的对比操作已经结束 至于如何调试 仁者见仁 智者见智 2 图形流程示意 图形流程示意 图 61 图 62 选择该项 选择该项 从 2 维模型中提取 数据 26 图 63 图 64 图 65 图 66 图 67 图 68 设置提取操作 的名称 选择要被提取的 目标文件 设置要被提取的内 容 设置被提取数据 的时间段 设置要被提取 的点数 指定被提取的点的位 置 设置输出文件名和位置 27 图 69 图 70 图 71 图 72 运行之后出上图内容 然后完成 至此 提 取指定位置数据的操 作已经完成 打开图形编辑 器 插入新的图形对象 选择该项 时间序列图形 28 图 73 图 74 图