3--MIKE 21 FM水动力模型逐步练习实例_pjf_2.21

MIKE 21 HD FM水动力模型逐步 培训 焊锡培训资料ppt免费下载焊接培训教程 ppt 下载特设培训下载班长管理培训下载培训时间表下载 教程 目 录 TOC \o "1-3" \h \z \u 1 概述 1 1.1 工程背景 1 1.2 练习实例的目的 1 2 创建计算网格 2 2.1 创建网格前需要注意的问 题 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 2 2.2 创建Øresund的计算网格 3 2.2.1 由原始的xyz数据生成mdf文件 4 2.2.2 三角边界的调整 7 2.2.3 模拟区域的三角划分 8 3 创建MIKE 21 FM水动力模型的输入条件 13 3.1 生成水位边界条件 13 3.1.1 把测量水位导入时间序列文件 14 3.1.2 创建边界条件 19 3.2 初始条件 22 3.3 风力作用 22 4 MIKE 21 FM模型搭建 23 4.1 FM模型 23 4.2 模型率定 37 4.2.1 实测水位 37 4.2.2 实测流速 38 4.2.3 模拟与实测结果比较 39 1​ 概述 本实例是连接丹麦和瑞典的跨海(Øresund)工程。 图 1.1 Sound (Øresund), 丹麦 1.1​ 工程背景 1994年，哥本哈根和马尔默(Malmö)开始了连接丹麦和瑞典隧道和桥梁的改造项目。该项目执行了严格的环境要求，即隧道和桥梁项目对波罗的海的环境不产生任何影响。这样的要求意味着桥梁和隧道设计的阻流作用小于0.5%，同理，溢流和排放的最大流量也要得到控制。为了达到环境的要求和监理工程施工，建立了一个主要的监测程序。整个监测程序包括40多个水文测站，收集水文、盐度、温度和流场数据。另外还为ADCP的船载测站和CTD等固定站点进行了广泛的补充测量。监测程序最初于1992年开始并一直持续到本世纪。 由于Øresund海域天然水文的多样性和多变性，连接工程的阻流作用只能通过数值模型来 评价 LEC评价法下载LEC评价法下载评价量规免费下载学院评价表文档下载学院评价表文档下载 。而且，Øresund的情况需要一个三维模型。所以，利用DHI的三维模型，MIKE 3对Øresund整个海域进行模拟，并在其中设置嵌套模型，网格尺寸水平方向由连接工程附近的100米到Øresund 较远海域的900米，垂直方向网格尺寸是1米。随后，MIKE 3模型会根据现场测量数据阶段进行率定和验证。 根据监测程序得到的数据，初步选择足以反映Øresund海域天然水文多样性的3个月作为模拟的“设计时段”。设计时段用来对连接工程进行详细的规划和优化，并确定需要填充的挖泥量，以达到对环境没有任何影响。 1.2​ 练习实例的目的 练习实例的目的在于通过使用MIKE 21的Flexible Mesh模块为Øresund建立水流模型和MIKE 3的水流模型，生成令人满意的率定结果。 此次练习和实际工程操作相同，但根据输入数据也做了一些预备工作，主要是为用户准备了MIKE 21格式的输入数据，以保证原始数据的准确性和预处理。根据数据的数量和质量，数据的处理是非常耗时但又必不可少的过程。本实例中，所有的原始数据都以ASCII文件的格式提供，所有相关数据和文件请在以下目录中查询： C:\Program files\DHI\2011\MIKEZero\Examples\MIKE_21\FlowModel_FM\HD \Oresund C:\Program files\DHI\2011\MIKEZero\Examples\MIKE_3\FlowModel_FM\HD \Oresund 如果用户对导入MIKE Zero格式数据的过程已经熟悉，则不用再自行生成所有的输入文件。在模型中确定了运行模型需要的所有输入条件后，用户就可以开始模拟。 2​ 创建计算网格 创建计算网格需要对数据进行大量的修改，在此，只对主要方法进行解释。 网格文件中包含不同地理位置的水深和下列信息： 计算网格 水深 ​ 边界资料 在建模过程中，网格文件的生成至关重要。 请在下列目录中参见网格生成器的用户手册： C:\Program Files\DHI\2011\MIKEZero\Manuals\MIKE_ZERO\MzGeneric.pdf 2.1​ 创建网格前需要注意的问题 地形和网格文件应该： 1)​ 描述模拟区域内的水深 2)​ 使模拟结果达到理想的精度 3)​ 达到用户能够接受的模拟用时 为了达到上述目的，用户需要注意的网格是： 1)​ 没有过小角度的三角形（完美的的网格是等边三角形） 1)​ 带有平滑的边界 1)​ 分辨率更高的嵌套部分 1)​ 来自于xyz数据的网格 较大的角度和较高的分辨率需要较长的计算时间，所以模型的使用者必须在计算时间和网格分辨率之间保持一个平衡。 网格的分辨率、水深、和时间步长决定了模型设置中的克朗值。最大克朗值应该小于0.5，所以，模拟用时不仅和三角网格相关，还由网格的节点数量和克朗值决定。因此，在深水区网格较细的地方所需要的计算时间就比浅水区较细的网格要长些。 2.2​ 创建Øresund的计算网格 实例中的模拟区域的海图如图2.1所示，它覆盖了丹麦和瑞典的海峡。根据海图可以得到水深和海岸线的xyz数据。在本实例中，海岸线和水深的xyz数据已经创建好了，见图2.14。 图2.1 模拟区域：Øresund, 丹麦和瑞典海峡 由原始的xyz数据生成mdf文件 网格文件包括水深资料，网格是通过MIKE Zero的网格创建器建立。首先，打开“Mesh Generator” (NewMesh Generator)，见图2.2。 打开“Mesh Generator”，用户需要确定模拟地区的投影方式为UTM，并定义UTM区为33，如图2.3所示。 图 2.2 MIKE Zero中的“Mesh Generator” 图 2.3 定义模拟区域的投影方式 定义好的模拟区域如图2.4所示。 图 2.4 添加xyz数据以前在“Mesh Generator”定义的模拟区域 从ACSII文件中导入海岸线的数据 (ImportBoundary data Open XYZ file: land.xyz)，见图2.5。 图 2.5 导入海岸线数据 请注意需要转换地理坐标：在导入了海岸线数据后请选择“Longitude/Latitude”，见图2.6。 图 2.6 导入海岸线数据 生成的模拟地区包括导入的海岸线数据被称为“Mesh Definition File”(mdf-file)， 如图2.7。 图 2.7 导入海岸线xyz数据后在Mesh Generator中的mdf文件 下一步是把模拟区域的原始数据变为可以划分三角网格的数据。 三角边界的调整 边界调整的目的是使模拟地区的边界能够被划分为三角网格，包括水边界（绿色弧形）。 从删除不需要包含于模拟区域的海岸线的交汇点和节点（红色和兰色点）开始进行调整，这还包括陆地上的节点，如图2.7所示。 通过在“丹麦”节点和“瑞典”节点之间添加弧形来定义北部和南部的边界。边界应当定义在测量的边界位置附近，如表2.1。 选择北部边界并选择“properties”，把北部的圆弧属性设置为“2”，南部边界的圆弧属性设为“3”，见图2.8和图2.9。这些属性用来区分网格中的不同边界类型：北部边界（2）和南部边界（3）。陆地/水面边界（1）由网格创建器自动生成。 表2.1 测量的水位数据 测站 数据文件 位置 东 (m) 北 (m) WL13 Viken waterlevel_viken.txt 349744 6224518 WL14 Hornbæk waterlevel_hornbaek.txt 341811 6219382 WL19 Skanør waterlevel_skanor.txt 362748 6143316 WL20 Rødvig waterlevel_rodvig.txt 333191 6126049 图2.8 编辑属性（右击鼠标）选择北部边界的圆弧（紫色） 图2.9 编辑北部边界圆弧的属性 现在用户能够看到一个可以划分三角网格的闭合区域，但首先需要对所有的海岸线进行核对和修正。模拟区域的三角网格划分从边界多边形开始，所以，网格三角划分过程中产生的元素数量和海岸线的交汇点和节点数量密切相关。用户可以使用相应的工具把各交汇点重新分布到圆弧上，使海岸线更为平滑均匀。在特别关心的地区，用户还可以把交点沿陆地边界重新分布，并使其相互保持较小的间距。用户需要注意，如果内陆港口和内陆湖对模拟地区影响不大，则不需要在模拟地形中予以考虑。 编辑完所有的海岸线后，用户可以得到和下图2.10所示相似的地形图，名为oresund.mdf。用户可以在随后的工作中使用此文件。 模拟区域的三角划分 下一步就是真正进行模拟面积的三角划分。首先，选择不适合进行三角划分的多边形（绿色）的闭合地区（如岛屿）。 图2.10 Mesh Generator 中的mdf文件， 边界xyz数据经过调整以形成能够进行三角划分的闭合面积 为了进行最初的三角划分(MeshTriangulateGenerate)，用户需要进行相应的设置，如图2.11 (Mesh TriangulateOptions)。在本实例中，特别需要关注的地区的网格分辨率较高。用户也可以根据需要把某些地区的网格变细，具体做法是在需要加密的网格上添加多边形并设置多边形的属性（在多边形中添加绿色标识并右击鼠标定义属性）。 图2.11 三角划分的选项 三角划分完成后，用户可以通过工具使网格平滑(MeshSmooth Mesh)。在这个例子中，网格被平滑了100次，如图2.12所示。 图2.12 三角划分和网格平滑100次后的网格。三角划分选项允许的最大面积是 1500000 m2，最小的角度是30度，最多的节点是6000个。加密地区的选项在加密地区的多边形选择。 随后，用户需要把xyz文件中的水深数据(water.xyz)内差到网格中 (MeshInterpolate)，见图2.13。注意在MIKE 21 FM和MIKE 3 FM中的水深不同，因为它们的测量年份不同，地形也有了变化。MIKE 21的mdf文件是1993年的数据，如图2.15。 图2.13 内差水深需要选择xyz的数据文件。定义投影方式为地理坐标1993 water.xyz （MIKE 21 Flow Model FM 例子）和UTM 33 1997 water.xyz（MIKE 3 Flow Model FM 例子) 图2.14 ASCII文件，定义了地理位置（经度、纬度和深度），请注意如果使用了MIKE C-MAP，用户则不能以文本方式浏览数据，因为数据已经加密。 图2.15 在内差水深后的Mesh Generator中的网格 现在用户可以把编辑好的网格数据导入到MIKE 21/3 FM的水流模型中 (MeshExport Mesh)。保存文件为oresund.mesh。 用户可以在Data Viewer（图2.16）或MIKE Animator （图2.17）中浏览或编辑网格文件。 图2.16 Data Viewer中的网格文件 图2.17 MIKE Animator 中的Øresund网格 3​ 创建MIKE 21 FM水动力模型的输入条件 在设置MIKE 21 FM水流模型前，必须根据测量数据创建输入数据，1993年测量的数据是： 1)​ 边界上的水位 1)​ 丹麦哥本哈根Kastrup机场的风数据 通过使用MIKE Zero中各种工具来准备输入数据，相应的用户手册在下列安装目录中可以找到： C:\Program Files\DHI\2011\MIKEZero\Manuals\MIKE_ZERO\MzGeneric.pdf 3.1​ 生成水位边界条件 在模型的边界附近有四个测站有测量的水位 记录 混凝土 养护记录下载土方回填监理旁站记录免费下载集备记录下载集备记录下载集备记录下载 ，见3.1。 Øresund模型是水位边界，测量的水位数据表明沿着边界的水位变化比较显著，所以水位边界应该定义为线边界（dfs1类型的数据文件），是两个边界点数据的内差。 随后，两条线边界（dfs1类型的数据文件）的水位数据来自于两条边界上的四个测站的测量数据。四个测站的位置如表2.1所示。 图3.1 开边界上的水位测站位置： Hornbæk，Viken，Skanör和Rødvig 把测量水位导入时间序列文件 图3.2 MIKE Zero中的时间序列编辑器 打开MIKE Zero中的“Time Series Editor”(FileNewTime Series)见图3.2。选择“ASCII格式”，打开文本文件 waterlevel_hornbaek.txt。选择“Equidistant Calendar Axis”并确认，然后右击鼠标并选择“properties”，把类型改变为“Water Level”，保存数据为waterlevel_hornbaek.dfs0。对其他三个测站重复上述步骤。 注意在MIKE 21/3 FM水流模型的现有版本中，必须是等时间步长，这就意味着原始数据如有缺失，必须在导入之前通过内差等方式补充。 图3.3 测站Hornbæk 的实测水位，ASCII文件 图3.4 时间序列编辑器的导入界面 图 3.5 时间序列属性 图3.6 由测站1 Hornbæk导入到时间序列中的水位数据 如果要把水位数据绘成时间序列图，打开MIKE Zero中的“Plot Composer”， 见图3.7。选择“plot” “insert a new plot object”并选择“Time Series Plot”，见图3.8。 图3.7 MIKE Zero中打开Plot Composer 右击绘图区域并选择“properties”，点击 添加时间序列文件到“Plot Composer”，见图3.10。用户可以在同一个绘图中添加多个时间序列图，在属性对话框中还能对图形进行编辑，如颜色等（见图3.9）。 图3.8 在Plot Composer 中插入新的绘图对象 图 3.9 在Plot Composer的Time Series Plot对话框中选择时间序列文件，调整比例和曲线等 图3.11和图3.12是两个边界上的实测水位。 图3.10 在Plot Composer中选择时间序列文件 图3.11 北边界测站1和2，Hornbæk和Viken 的时间序列 图3.12 南边界测站3和4，Skanör和Rødvig的时间序列 创建边界条件 下一步是根据生成的时间序列创建线边界。打开MIKE Zero中的“Profile Series”并选择“Blank ...”，见图3.13。 图3.13 在MIKE Zero中打开Profile Series Editor 需要的信息如下： 北边界 开始时间1993-12-02 00:00:00 时间步长：1800s 时间步骤数：577 网格点数：2 ​ 网格间距：9200m （边界的宽度，事实上并不需要，因为MIKE 21/3 FM的水流模型中线边界的内差不考虑此宽度，见MIKE 21/3 FM 水流模型的用户手册）。 打开测站Viken (waterlevel_viken.dfs0) 的水位文件，复制并粘贴到“profile Series Editor”中点“0”的位置，然后打开测站Hornbæk (waterlevel_hornbaek.dfs0)的的水位文件，复制并粘贴到“profile Series Editor”中点“1”的位置，见图3.15。保存文件为“waterlevel_north.dfs1”，见图3.16。 图3.14 剖面序列属性 图 3.15 从测站Viken复制水位到 Profile Series Editor(Ctrl V) 图3.16 北边界的水位线边界 除了网格间距不同，和使用的水位数据来自测站Rødvig(waterlevel_rodvig.dfs0)和Skanör(waterlevel_skanor.dfs0)，对南边界重复相同的步骤和类似的信息，并保存结果为waterlevel_south.dfs1。 南边界 开始时间 1993-12-02 00:00:00 时间步长：1800s 时间步骤数：577 网格点数：2 ​ 网格间距：33500m （边界的宽度，事实上并不需要，因为MIKE 21/3 FM的水流模型中线边界的内差不考虑此宽度，见MIKE 21/3 FM 水流模型的用户手册）。 3.2​ 初始条件 初始的水面高程由北边界和南边界水位的平均值计算得到。从两个边界文件中看平均的高程约为 –0.37 m。 3.3​ 风力作用 来自Kastrup机场的风数据形成的风力是时间序列，但空间上恒定。打开时间序列编辑器，从ASCII文件“wind_kastrup.txt”中导入数据，使用等时间步长。保存文件为“wind_kastrup.dfs0”。风速和风向的时间序列如图3.17所示。 对于风的更形象描述是风玫瑰图。打开“Plot Composer”，插入新的绘图对象，并选择“Wind/Current Rose Plot”，然后选择属性，选择新生成的文件“wind_kastrup.dfs0”并按需要对图形进行编辑（颜色等）。结果如图3.18所示。 图3.17 来自Kastrup 机场的风速和风向数据，ASCII文件 图3.17 Kastrup机场的风速和风向数据， 由Plot Composer (Time Series Direction plot control）绘制 图3.18 Kastrup机场的风玫瑰图，由Plot Composer绘制，南边界 4​ MIKE 21 FM模型搭建 4.1​ FM模型 现在我们采用第三章生成的1993年Øresund地形资料，边界条件和动力数据搭建MIKE 21 FM水动力模型。首先，我们采用模型默认参数，不考虑边界密度变化的影响。模型参数如下表所示： 参数 值 模型配置文件 Oresund.m21fm 网格地形 Oresund.mesh (1993) 2057个节点 模拟时段 1993-12-02 00:00 – 1993-12-13 00:00 （11天） 时间步长 8 秒 模拟步长数 118800 干湿水深 干水深0.01 米 洪水深0.05 米 湿水深0.1 米 初始水位 -0.37 米 风场 随时间变化，在模型范围内常量：wind_kastrup.dfs0 风摩阻系数 随风速变化： 0.001255 ，风速为： 7 m/s时 0.002425 ，风速为 25 m/s时 北部边界 Waterlevel_north.dfs1 南部边界 Waterlevel_south.dfs1 涡粘系数 Smagorinsky 公式 小学单位换算公式大全免费下载公式下载行测公式大全下载excel公式下载逻辑回归公式下载 , 常量 0.28 糙率 曼宁系数. 常量 32 m1/3/s 结果文件 Flow.dfsu ndr_roese.dfs0 模拟时段 大约25分钟， 2.4 GHz CPU, 512 MB DDR RAM 下面简短地介绍了模型的搭建过程。 在模型范围对话框内指定地形网格文件oresund.mesh，如图4.1所示。此地形文件已经内含了模型投影UTM-33。在地形文件中，每个边界被给定一个识别代码。在Øresund例子中，北部边界识别代码为2，南部边界识别代码为3。在边界名称对话框中将“代码2”改名为“北部”，“代码3”为“南部”。 注意，点击 按钮，可以直接打开并查看网格文件。用户可以在Data Viewer中查看代码值（第二项），如图4.2所示。 图4.1 MIKE 21 FM模型：指定模型地形文件 图 4.2 在Data Viewer中查看边界代码值 在“Time”对话框中指定时间步长为8秒，输入时间步数为118800步，对应的模拟时间11天，见图4.3。 图4.3 MIKE 21 FM模型：模拟时段 如果要同时计算对流扩散或Ecolab水质，请在“Module Selection”对话框对应的选择框中打勾，见图4.4。 在本例中，我们只进行水动力模拟。如果希望在水动力模块的基础上增加Ecolab水质模块，请参考“Step-by-Step Training Example for ECO Lab”。 图 4.4 MIKE 21 FM模型：模块选择 在“Solution Technique”对话框里，设置最小时间步长和最大时间步长分别为0.01s和120s。临界CFL数设置为0.8。 在“Flood and Dry”对话框中，可以进行干湿动边界的参数设置，见图4.5所示。在该例子中，选取默认值，即干水深0.01米，洪水深0.05米，湿水深0.1米。 在该例子中，由于沿着Saltholm海岸的某些区域在模拟期间内会干涸，所以应该包含干湿动边界功能。如果没有选择该功能，模型会在干涸地区溢出。但是，包含干湿功能后会影响模型的稳定性，因此，如果干湿功能对所要模拟的问题不是很重要的话，建议不要将其包含进来。做为补偿，请调整网格地形文件，使浅水区域有更深的水深，从而防止这些区域干涸，然后运转不包含干湿功能的模型。 图4.5 MIKE 21 FM水动力模型：干湿动边界 本例中不考虑密度梯度流，在“Density”对话框中选择“Barotropic”（正压），见图4.6。 图 4.6 MIKE 21 FM水动力模型：密度梯度 水平涡粘系数的设定采用Smagorinsky公式，Smagorinsky常数取默认值0.28，见图4.7。 图4.7 MIKE 21 FM水动力模型：涡粘系数的设定 第一次模型运行采用默认的河床糙率值（曼宁系数）32 m (1/ 3)/s，在随后的率定中可以改变它，见图4.8。 图4.8 MIKE 21 FM水动力模型：河床糙率 在Øresund地区潮流作用占优，一般来说，应该考虑科氏力的作用。但是在该例中，由于海峡很浅，科氏力的影响也可以忽略。 而在实际的工程应用中，除非个别情况，一般都要考虑科氏力的作用，仅仅在实验室模拟中科氏力可以不考虑，见图4.9。 图4.9 MIKE 2 FM水动力模型：科氏力 使用已生成的风时间序列（在Wind Forcing对话框中风时间序列指定为“随时间变化，全域内常量”，时间序列文件名为wind_kastrup.dfs0）。在本例中，风场软启动时间间隔指定为7200秒。软启动时间间隔指在模拟开始阶段风的影响不完全考虑。在软启动间隔开始阶段，风力的影响为零，然后它逐渐增长直到在软启动间隔终止阶段充分考虑风力的影响。在本例中，指定风摩擦随风速变化，并采用默认的风摩擦值，见图4.11。 注意，在风场对话框中点击 按钮，可以直接打开并查看风场文件，见图4.10。 图4.10 MIKE 2 FM水动力模型：风场 图4.11 MIKE 2 FM水动力模型：风摩擦系数 在本例中 不考虑冰盖影响 不考虑降雨和蒸发 ​ 不考虑波浪辐射应力 源和汇的流量及速度在“源”对话框中指定。在本例中，由于Øresund地区的源流量过小，不会对水动力条件产生大的影响，因此没有包含在该例子中。此页设置为空白，见图4.12。 图4.12 MIKE 2 FM水动力模型：源项 在“Initial Conditions”对话框中设置初值条件，即模拟开始时刻的水位/流速。本例中初始水位采用常量-0.37米，即模拟初始阶段北部和南部边界水位的平均值，水平方向和垂直方向流速为0，见图4.13。 图4.13 MIKE 2 FM水动力模型：初始条件 在“Boundary Conditions”中边界条件应该与域对话框中指定的边界名一致。本例中北部边界和南部边界使用第三章生成的线序列水位边界文件。 本例中，边界类型是“Specified Level”（水位），这是因为边界上仅有水位实测结果可以利用。“Specified Level”意味着边界上采用了水位实测数据，而流量数据是未知的，只能在模拟过程中估算。如果边界类型是“Specified discharge”，意味着边界上采用了流量实测数据，水位数据是估算的。 边界格式必须设置为“Variable in time and along boundary”（随时间和边界变化），这种情况下，需要给定一个线序列文件。 点击 按钮，在弹出的对话框中选择相应的数据文件，见图4.16。 对于北部边界，选择waterlevel_north.dfs1水位文件；对于南部边界，选择waterlevel_south.dfs1水位文件。 注意：当在边界上指定线序列时，我们必须知道MIKE 21 FM是如何定义边界起始点和结束点的。定义规则是：以边界海岸线为界，使模型边界处于左手方向，沿海岸线第一个点定义为起始点，见图4.2。 软启动间隔取7200秒，对应初始水位值为–0.37米。软启动间隔是指在模拟开始阶段边界水位的影响并未充分考虑。在软启动间隔开始阶段，边界水位的影响是零，然后该影响逐渐增大直到在软启动终止阶段充分发挥作用，见图4.14和图4.15。 用户可以在边界对话框中点击 按钮，快速打开并查看边界数据文件。 图4.14 MIKE 21 FM水动力模型：北部边界条件waterlevel_north.dfs1 图4.15 MIKE 21 FM水动力模型：南部边界条件waterlevel_south.dfs1 图4.16 MIKE 21 FM水动力模型：选择边界文件 指定输出格式为面序列，并给出结果文件名，见图4.17。输出文件名称为flow.dfsu。同时，确保硬盘空间具有足够的存储空间。选择输出频率为3600秒（对于潮流模拟来说该频率较合理），这样可以使输出的面序列文件大小减小到合理范围。本例中，时间步长是8秒，指定输出频率为3600/8 = 450。默认情况下，面序列输出范围为全域。 选择输出文件要包含的参数，见图4.18。 在水文站Ndr. Roese处输出一个点序列文件作模型率定用，该站的位置见表4.1。如果有兴趣，您也可以输出监测站附近的一些点，以便了解监测站附近的计算结果随空间位置的变化。 图4.17 MIKE 21 FM水动力模型：结果保存为点，线和面序列 图4.18 MIKE 21 FM水动力模型：输出结果参数 表 4.1 Ndr. Roese站测量结果 水文站 数据文件 位置 东 （米） 北 （米） Ndr. Roese waterlevel_ndr_roese.txt currents_ndr_roese.txt 354950 6167973 现在我们准备运行MIKE 21 FM水动力模型 (RunStart simulation…)。 用户也可以运行DHI为您预先准备好的模型配置文件： C:\Program Files\DHI\2011\MIKEZero\Examples\MIKE_21\FlowModel_FM\HD \Oresund\Calibration_1\oresund.m21fm 注意，如果模拟不成功，请打开log日志文件找出错误原因(FileRecent log file list)；您也可以在模拟运行之前在菜单栏的“View log file after simulation”前打勾，这样模拟结束后，会自动弹出Log日志文件。 模拟结束后使用Plot Composer（或 Data Viewer）检查和显示结果，见图4.19和图4.20，其中一个是北向潮流，一个是南向潮流。 dfs2和dfs3的后处理工具（数据提取，统计等等）也能被dfsu文件使用。这要求先将dfsu文件转换成dfs2或dfs3。这里也有相应的转换工具可用。 NewMIKE ZeroGrid SeriesFrom Dfsu File 图 4.19 北向水流流速与水位 图4.20 南向水流流速与水位 4.2​ 模型率定 模型率定，需要模型范围内的实测数据。表4.1给出了实测水位与实测流速站的位置和数据文件。 实测水位 水文站Ndr. Roese (waterlevel_ndr_roese.txt)提供了实测水位资料，使用时间序列编辑器导入这些ASCII文件（如何输入见章节3.1.1），制作成时间序列文件，图4.21为Ndr. Roese站实测水位。 图4.21 Ndr. Roese：实测水位 实测流速 水文站Ndr. Roese (currents_ndr_roese.txt)提供了实测流速资料，使用时间序列编辑器导入这些ASCII文件，制作成时间序列文件，图4.22－4.24为Ndr. Roese站流速与流向示意图。 图4.22 Ndr. Roese：实测流速 图4.23 Ndr. Roese：实测流速和流向 图4.24 Ndr. Roese 流速玫瑰图（也能在Plot Composer看到） 模拟与实测结果比较 采用Plot Composer可以同时显示模拟与实测水位，本例中已做好一个Plot Composer文件，打开文件ndr_roese.plc便可，见： C:\Program Files\DHI\2011\MIKEZero\Examples\MIKE_21\FlowModel_FM\HD \Oresund\Calibration_1\Plots 图4.25显示了Ndr. Roese水文站实测与模拟水位、流速比较图，实测与模拟结果比较说明率定结果令人满意。 改变河床糙率为45，Smagorinskyd常数为0.24，再次运行，图4.26显示实测结果与模拟结果更加匹配。 通过改变糙率和涡粘性系数进一步率定模型，每次率定后比较模拟与实测结果，每次仅仅改变一个参数，同时确保参数要在一定范围内赋值。 图4.25 Ndr. Rose水文站采用默认参数得出的实测与模拟水位、流速 图4.26 Ndr. Rose水文站率定后得出的实测与模拟水位、流速