mike3fm泥模块用户使用手册

1、目录目 录21 关于本文档 (About this guide)41.1 编写本文档的目的(Purpose)41.2 用户所需的知识背景 (Assumed User Background)41.3 MIKE Zero 模型编辑器程序界面结构(General editor layout)41.3.1 导航栏(Navigation tree)41.3.2 编辑窗体(Editor window)41.3.3 有效性显示窗体(Validation window)5系统 (Online help)51.42 泥模块(MUD TRANSPORT MODULE)62.1 参数的选取(Parameter Se

2、lection)62.2 求解格式 (Solution Technique)62.2.1 备注与提示（Remarks and hints）62.3 水体参数(Water Column Parameters)72.3.1 无粘性砂组分(Sand fraction)72.3.2 沉速 (Settling)82.3.3 淤积 (Deposition)132.3.4 涡粘系数与密度(Viscosity and Density)152.4 底床参数(Bed Parameters)152.4.1 概述 (General description)152.4.2 侵蚀 (Erosion)162.4.3 床层的

3、密度(Density of bed layers)182.4.4 床面糙率(Bed roughness)192.4.5 床层间的转换(Transition between layers)202.5 地貌 (Morphology)202.5.1 概述 (General description)202.5.2 备注与提示(Remarks and hints)212.6 外力 (Forcings)212.6.1 概述 (General description)212.6.2 波浪 (Waves)212.7 扩散 (Dispersion)242.7.1 平流扩散(Horizontal dispersi

4、on)242.7.2 垂向扩散(Vertical dispersion)252.7.3 推荐设置(Recommended values)252.8 源 (Source)252.8.1 源的设定(Source specification)252.8.2 备注与提示(Remarks and hints)262.9 初始条件(Initial Conditions)262.9.1 组分浓度(Fraction concentration)262.9.2 床层厚度 (Layer thickness)27MIKE BY DHI 20082.9.3 粒径组的分布(Fraction distribution)2

5、72.10 边界条件(Boundary Conditions)272.10.1 边界设置（Boundary Specification）282.11 输出(Outputs)292.11.1 地理视图(Geographical View)292.11.2 输出设定(Output specification)292.11.3 输出项目(Output items)32MIKE BY DHI 20081 关于本文档 (About this guide)1.1 编写本文档的目的 (Purpose)编写本文档的主要目的在于指导用户使用MIKE3 Flow M模块主要用于模拟粘性泥沙的输运。FM模型中的泥模

6、块，该1.2 用户所需的知识背景 (Assumed User Background)尽管泥模块有着合理且友用户界面、详尽的用户使用手册和文档，但对于程序而言，模型的调试和模拟结果正确性的相关理论背景知识。方面仍然要求用户具备一些泥沙输运的力模块用户对 Mike3 的基本组件已达到熟悉的程度。如相关的各种类我们亦型的文件及其编辑器、绘图工具（Plot Composer）、Mike Zero Toolbox、数据查看工具（Data Viewer）以及网格生成工具（Mesh Generator）。上述 Mike3 基本组件的说明文档详见C:Program FilesDHIMIKEZeroManual

7、sMIKE_ZEROMzGeneric.pdf。1.3 MIKE Zero 模型编辑器程序界面结构 (General editor layout)MIKE Zero模型编辑器包含三个组成部分。1.3.1 导航栏 (Navigation tree)模型编辑器的左侧为导航栏，导航栏内以树形结构显示模型配置选项。用户选中导航栏内树形结构中的任一条目，相应的设置界面即会在中间的编辑窗体中显示。1.3.2 编辑窗体 (Editor window)导航栏中所有条目对应的编辑界面均在间的编辑窗体中显示。该窗体与用户选择的条目相对应，且可能包含多个属性页面。对于空间数据相关的条目如源汇项（sources），边

8、界（boundaries）以及输出结果（output）等，编辑窗体中会相应显示相关条目的空间配置状况。用户可以在编辑窗体的底部选择这种空间显示方式的浏览操作，如放大、缩小或返回区域中点等操作；提供了一个文本菜单项供用户选择是否显示图形化的地形、网格、GIS 背景图层或图例；从该文本菜单项中用户亦可进行查看前一操作、后一操作或显示全部范围的屏显操作；此 外，如果用户在正在图形化显示某一条目的时候选中另一条目，则编辑窗体中的图形化显示界面将自动更新为新选项对应的内容。4MIKE BY DHI 20081.3.3 有效性显示窗体 (Validation window)在程序界面的底部为有效性显示窗体

9、，其功能为显示用户所进行的模型设置的有效性检验结果。它可以动态实时显示当前用户完成的设置操作的有效性检验结果，当检验结果中发现错误时，用户双击该错误显示则自动将产生错误的条目置为当前条目。1.4系统 (Online help)调出系统的有多种，用户可按照个人需要选择任意法使用：l按下F1功能键寻求帮助按下F1功能键即可将当前条目的相关帮助调出，如图1.1所示。l在帮助页面中手工打开系统在帮助菜单中，选中“help topics”菜单项即可调出系统。图 1.1 MIKE 3 Flow MFM系统5MIKE BY DHI 20082 泥模块 (MUD TRANSPORT MODULE)泥模块的主要

10、功能在于模拟力模块计算所得水流条件下的粘性泥沙输移过程。2.1 参数的选取 (Parameter Selection)用户需设定泥沙分组数和底床分层数模型中模型中的最大泥沙组数为 8.的最大底床分层数为 12.在保证底床分层能够充分体现沉积物垂向变异规律的前提下，建议尽可能地减小底床分层数目。2.2 求解格式 (Solution Technique)模型计算的时间和精度取决于计算数值所使用的求解格式精度。模型计算可以使用低阶（一阶精度）或是高阶（精度）的。低阶计算快但计算结果但精确度较差，高阶的册。计算精度高但速度较慢。更为详尽的关于数值计算的，请参考科学背景手浅水方程的时间积分和输移(扩散)

11、方程是基于半隐格式求解，相应平流项采用显式格式求解，而垂直对流项则采用格式求解。受显式格式性的限制，为保持模型计算的性，模型中时间步长的设定必须保证CFL数小于1，为保证所有网格点CFL数均满足该限制条 件，模型中时间步长的取值采用一浮动范围的方式，因此模型中用户需设定一最小和最大时间步长范围，相应扩散方程的时间步长在模型的计算过程中自动与主时间步长相匹配。用户可在Hydrodynamic M的求解格式CFL数的设置。框中进行最小和最大时间步长范围以及临界2.2.1 备注与提示（Remarks and hints）在所模拟的物理过程中，如果对流占优，则应选择较高阶的空间离散格式。如果扩散占优，

12、则较低阶的空间离散格式就可以满足模拟所需精度。应选择同样的计算精度格式。来说，时间积分和空间离散通常模型计算中采用高阶时间积分的计算时间是低阶的两倍；而采用高阶的空间离散所耗计算时间为采用低阶的1 ½ 到 2 倍。 若同时选择高阶的时间积分及空间离散，所耗计算时间将会是同时选择低阶时的3-4倍。来说采用高阶的计算结果的精确性通常会高于采用低阶的计算结果。模型中CFL数的程序默认设置为1。而言CFL数小于1时，模型即可保持计算的性。但因实际计算过程中CFL数的数值为近似预估值，故在这种默认设置情况下仍然发生模型计算失稳的可能性。因此当这种情况发生时，用户可将临界CFL值适当减小（取值范

13、围介于0到6MIKE BY DHI 20081之间），此外用户亦可适当减小所设定的最大时间步长。必须指出，当用户将最小和最大时间步长均设定为与主时间步长相同时，模型将以恒定时间步长进行计算，此时为保证计算的性，相应时间步长的取值必须要满足CFL值小于1。计算的总时间步数、最大最小时间步长均会文件中。在log文件中，而CFL数则可以输出至结果采用高阶进行计算时可以模拟出陡坡处的水流过冲和下冲现象。因此，当模型计选项时，计算过程中将无法保证质量守恒。算时同时选择了高阶和最大最小质量浓度2.3 水体参数 (Water Column Parameters)本选项中需要用户设定的各水体参数涉及到几乎所有

14、关于水体中泥沙过程的问题。这里必须设定模拟的是纯粘性的泥沙或者有一部分粒径组是非粘性的（沙）。通常认为粒径大于60 m m的泥沙是非粘性的。水体参数主要以下几个部分：····无粘性砂组分沉降特性淤积特性粘性系数和密度2.3.1 无粘性砂组分 (Sand fraction)描述 (Description)泥沙输运以力为基础。泥沙输运方式可分为两种，粘性和非粘性。粘性泥沙的沉速小，对力变化的响应慢。它在水体中的输移主要以平流为主。而非粘性泥沙的沉速较大，浓度分布对层。力变化的响应也较快，所以大部分的非粘性泥沙会以推移质的形式在近底MIKE 3 MT能模拟细颗

15、粒非粘性泥沙的悬沙输移。该功能是在已知泥沙特性和下，通过计算平衡状态下的浓度分布来实现的。力条件模型中床面侵蚀以成层冲刷的方式进行，即认为当床面侵蚀发生时，床面上处于同一薄层内的沉积物是以无分选性起动的方式被水流冲刷带走的。也就是说， MT模块中的侵蚀公式计算的是所有粒径组的最大冲刷量。模型中亦被整体冲刷的薄层沉积物进入水体以后会在水流紊动作用下或重组。当然由于非粘性泥沙颗粒间粘滞力极为微弱的缘故，上述不适用于非粘性泥沙组分，非粘性泥沙组分输移、沉降过程主要通过计算给定泥沙特性和力条件下的水流挟来完成。如果水体中的非粘性泥沙浓度大于相应水流挟，则相应多余的泥沙即会淤积至床面，即所求水流挟为不冲

16、不淤平衡状态下的含沙量。的参见MIKE 3 MT 科学背景手册推荐值 (Recommended values)7MIKE BY DHI 2008泥沙颗粒的平均沉速可通过沉速公式来估算。参见Settling Velocity。备注与提示 (Remarks and hints)不推移质。如果含有非粘性砂组分，用户必须确保该组分主要以悬沙形式输移。2.3.2 沉速 (Settling)悬沙的沉速主要可分为四种：····等速沉降絮凝沉降干扰沉降浮泥如果选择等速沉降，采用等速沉降值。如果选择絮凝沉降，还可以选择是否考虑干扰沉降，并有两种干扰沉速公式可供选择，Ric

17、hardson和Zaki（1954）公式或Winterwerp（1999）公式。概述 (General Description)沉速决定于颗粒的大小。可通过公式粗略估计单颗粒泥沙的沉速：(r - r)gd 2=sws18× rv其中，r s ：泥沙密度(kg/m3)(石英=2650 kg/m3)r ：水的密度g ：重力度(9.82m/s2)d ：粒径 (m)v ：粘度 (m2/s)ws ：沉速(m/s)如果是细颗粒粘性泥沙(<0.006mm)，沉降颗粒的粒径和沉速取决于絮凝率。当水体中悬沙浓度较低时，粘性细颗粒泥沙颗粒间发生碰撞的概率较低，此时泥沙颗粒沉速接近颗粒泥沙的沉速。而

18、随着水体中悬沙浓度的增高，粘性泥沙颗粒间的碰撞亦相应更8MIKE BY DHI 2008为频繁，在颗粒间粘滞力的作用下粘性细颗粒泥沙会吸附在一起形成粒径较大的絮凝体。这会导致颗粒/絮团的粒径和沉速增大。图2.1 典型的沉速变化还有许多其它因素可以增大或减小絮团的粒径。在09psu之间的盐度能增大其粒径，高浓度的有机物也能增大絮团。高强度的紊动能破坏絮团的结构从而减小其粒径。随着悬沙浓度的进一步增大，最终絮团在沉降的过程中开始影响水流结构。实际上，絮团在沉降过程中会引起一股上升流，这股上升流会平衡或阻碍泥沙的沉降过程，从而导致泥沙颗粒的沉速降低。含沙量的进一步增大则会使絮团颗粒间的间距更加狭小，

19、此时水体中的悬沙性质类似于浮泥，颗粒以群体形式沉降，沉速大幅变小甚至微弱到可以忽略不计。参数取值(data)悬沙的沉速可分为四种：····等速沉降絮凝沉降干扰沉降浮泥如果选择等速沉降，采用等速沉降值。9MIKE BY DHI 2008如果选择絮凝沉降，还可以选择是否考虑干扰沉降，并有两种干扰沉速公式可供选择，Richardson和Zaki（1954）公式或Winterwerp（1999）公式。等速沉降 (Constant settling velocity)如果含沙量不影响沉速，那么可以选择等速沉降。如果选择等速沉降，沉速在模拟过程中保持不变且不随含沙量

20、改变。絮凝 (Flocculation)如果含沙量足够高，能使絮团影响互相的沉速。这是因为絮团间的碰撞会增大其粒径，从而增大其沉速。在这种情况下要选择絮凝沉降。程序默认情况下的絮凝含沙量是忽略干扰沉降的。计算参见图2.2图2.2 选择絮凝沉降时的含沙量分布r se dim ent颗粒的密度c floc絮凝初始时刻的含沙量ctotal总含沙量（所有粒径组的含沙量总和）10MIKE BY DHI 2008chindered干扰沉降初始时的含沙量Ws沉速W0沉速系数g幂（常数）干扰沉降 (Hindered settling)当含沙量足够高，使絮团能影响互相的沉速，并使其不能种情况下要选择干扰沉降。沉

21、降，导致沉速变小。在这当指定的干扰沉降的含沙量大于干扰沉降时的初始值，计算参见图2.3。图2.3 选择干扰沉速时的计算这里有两个沉速公式可供选择。Richardson和Zaki（1954）公式 (Formulation by Richardson and Zaki (1954)对组粒径，标准的Richardson和Zaki公式为：11MIKE BY DHI 2008öws ,næcws = ws,r ç1 -÷ç÷cègel ø对于多组粒径，Richardson和Zaki公式可扩展为：= wi (1- F )iww

22、is ,nss,r*其中，F* = min(1.0,F)åciF =icgelws,r沉速系数ws,n粒径组的幂常数cgel絮凝临界含沙量Winterwerp（1999）公式 (Formulation by Winterwerp (1999)(1- F )(1- F )*pwi= wiss,r1 + 2.5F其中，åciF= iprs这里的 r s 是泥沙的干密度。浮泥 (Fluid mud)在这个模型中的浮泥是作为底边界层来考虑的，其沉降过程作为固结过程处理。沉速随盐度变化的调整(Modification of settling velocity due to saliv

23、ariation)在淡水/微咸水中，絮凝过程会降低，这会影响沉速。由于絮团粒径较小，沉速会降低。将沉速乘以一个比例系数来模拟这种情况。12MIKE BY DHI 2008= w (eS×C2 )wii1-Css1其中C1 和C2 是率定系数。C1 和C2 没有在菜单中显示，它们的默认值分别为C1 =0.5， C2 =-0.33。推荐值 (Recommended values)下面的表格粗略的描述了不同范围的沉降过程。表2.1 沉降范围备注与提示 (Remarks and hints)因考虑了盐度对絮凝的影响，所以在函数）。力模拟中必须考虑密度的变化（函数或盐度时，盐度对絮凝的作用不再

24、增大。系数C1 和C2 的默认设置要注意的是当盐度超过 10psu分别为 0.5 和-0.33，且该设置不在菜单中显示。 ws,n 和g 也同样如此，它们的默认值为 1。2.3.3 淤积 (Deposition)概述 (General Description)水体中泥沙的淤积过程是指泥沙从悬沙变为底床沉积物的转换过程。当水流床面切应力小于泥沙临界淤积切应力发生淤积。第 i 组粒径的泥沙淤积量可表述为：= wi ci piDis b D其中 pi 是淤积速率的斜坡函数， w 是沉速， c 是第 i 组粒径的泥沙在近底层含沙量。Dsb淤积速率的斜坡函数定义为：13MIKE BY DHI 2008&

25、#230;ööæt bpD = maxç 0, minç1,1-÷÷iç÷tiècd øøè含沙量分布 (Concentration profiles)Teeter profile基于 Teeter profile 公式，可以描述近底含沙量ci 与垂线平均含沙量ci 的：bæöpic = c i ç1+e÷ç2 5 ÷b1.25 + 4.75 pièDø其中， pi 是与第 i 组粒径泥

26、沙相应的Peclet 数，定义为：ewipi= 6s kUef其中，U f 是摩租流速， k 是冯·常数，取为 0.4。Rouse profile通过 Rouse profile 公式可以描述近底含沙量ci 与垂线平均含沙量ci 的。近底含沙量b可定义为：c icb = RC其中， RC 是质心的相对高度，定义为从河床到含沙量的质量中心的距离与水深之比。它不随时间变化，所以近底含沙量分布是恒定的。沉积率的计算为:D = cb × ws , 当t b <t cd参数取值 (Data)临界淤积切应力可以两种方式设定· 在模型区域内为常数· 在模型区域内

27、为变化量如果是第二种情况，则必须准备一空间上至少完整模拟区域范围且包含相关参数的 dfs2 或 dfsu 文件。当采用 dfsu 文件时，相应参与计算网格点上的数据采用分段常数插值方法生成，而当采用dfs2 文件时，则采线性插值。14MIKE BY DHI 2008推荐值 (Recommended values)起动拖曳力通常是一个校准参数。这个值0.1N/m2。相对的质心高度通常近似 0.3。小于侵蚀剪切力。它的取值范围是 0-备注与提示 (Remarks and hints)如果临界淤积切应力数值越大，则泥沙的淤积量越大。反之亦然。2.3.4 涡粘系数与密度 (Viscosity and

28、Density)如果在力模块中的密度类型选项中没有选择正压，用户就可以在泥沙输运计算中设定有关参数，用以反馈给力计算中的涡粘系数和密度。用户必须设定悬沙的密度，以及反馈给涡粘系数的基础和参考含沙量。备注与提示 (Remarks and hints)只有含沙量很高时，这些反馈会对力计算有影响。2.4 底床参数 (Bed Parameters)底床参数定义为河床变化过程的参数。这里必须设定每层的侵蚀过程、密度和床面糙率参数。河床参数主要一下几个部分：····侵蚀过程每层河床的密度床面糙率每层河床间的转换2.4.1 概述 (General descriptio

29、n)泥模块中的河床有一层或多层。每一层都需要定义河床层内的泥沙总量，干容重、湿容重以及抗冲性。每层河床的泥沙由该层内所有粒径组的泥沙组成。每层泥沙的量被认为是变化 的，这就意味着该模型在模拟过程中可以反演床沙颗粒组成的时空变化。每层的干容重和抗冲性在时间上是不变的。河床层被视为“功能”层，它的每一层由其干密度和侵蚀性来描述，而不是物理层，物理层的物理性质会因固结或者其它过程随时间变化。在之前的河床描述中 ，固结过程被表述为泥沙从河床的一层转换到另一层的过程。河床的第一层（最上层）为“最弱”层，主要是浮泥或新淤积的泥沙，其下层面的密度及强度都不断递增。图 2.4 是一个河床（两层）的示例，描述了

30、输运过程对其的影响。 在模拟过程中，河床的一层或多层被完全侵蚀，以至于在某些地方该层变为“空白”。在特定时15MIKE BY DHI 2008间、特定地点的活动层被定义为从顶端数起的第一层，该层是非“空白”的。侵蚀总是从活动层开始。淤积下来的总是落淤于最上层。图 2.4泥模块过程。 床层在水体河床交界面以下（点线）第 j 层河床层，某水平网格上的第新。i 组粒径组的泥沙量按照下面的表达式每一步都做更+ (D - E )Dt + (T- T )m= mnewoldi, j-1i, ji, ji, jii其中， m (kg/m2) 是泥沙量， D (kg/m2/s) 是一个可能淤积量（ 仅在河床最

31、上层）， E (kg/m2/s)是一个可能侵蚀量（仅从活动层开始）， Ti (kg/m2/s)是一个可能向下转化的泥沙量， Dt (s)是模拟的时间步。第 i 层河床层的厚度从以下的计算式得到：ånewmM ji, jH new = irrjd , jd , j其中， H (m)是河床层的厚度， M (kg/m2)是总泥沙量， r d (kg/m3)是干密度。2.4.2 侵蚀 (Erosion)河床层的侵蚀是指从泥沙从河床向水体转换。侵蚀从床面活动层开始（ 参见 BedDescription），当床面水流剪切力(t b )大于临界侵蚀剪切力(t ce )发生。每层的河床参数都是常数。

32、所以仅在活动层计算侵蚀。被侵蚀的物质会根据其在河床的分布被分散到不同的粒径组。临界侵蚀剪切力 (Critical shear stress for erosion)侵蚀的标准是床面剪切力大于临界侵蚀剪切力，相应的推动力超过了力。在模拟过程16MIKE BY DHI 2008中临界侵蚀剪切力是常数。硬底床描述 (Hard bed description)对于一个致密固结的河床，第 j 层的侵蚀率被表述为(Partheniades, 1989)：EE j= E pjjm0E其中， p j 是一个侵蚀率的斜坡函数， E 是侵蚀系数， E 是幂。0mEæt böpE = max 0

33、,-1÷çjtièceø软底床的描述 (Soft bed description)对于一个软式的，部分固结的河床而言，第 j 层的侵蚀率被表述为(Parchure&Mehta, 1985):E = E exp(a (t -t )jj0bce其中，a 是一个系数。参数取值 (Data)每层的侵蚀参数必须被设定为一个常数（时间上）。每层的临界侵蚀剪切力可被设定为：· 在整个模型范围内为常数· 在整个模型范围内为变化量如果是第二种情况，则必须准备一空间上至少完整模拟区域范围且包含相关参数的 dfs2 或 dfsu 文件。当采用 df

34、su 文件时，相应参与计算网格点上的数据采用分段常数插值方法生成，而当采用dfs2 文件时，则采线性插值。推荐值 (Recommended values)E 的值是一个0.00002kg/m2/s 之间。侵蚀速率的比例因子。对软床而言， 其值在 0.000005 到如果设置了 Em 或a 的值，那么侵蚀率会以指数形式演变。a 的取值范围通常在 4 到 26之间。对于临界剪切力，按下表取值：17MIKE BY DHI 2008表 2.2 临界剪切力临界侵蚀剪切力可按照如下水流有效切应力得出：t ce =0.00001t yjt yj £ 0.00015N / m2t ce = 0.00

35、00025t yjt yj> 0.00015N / m2其中，t yj 是第 j 层的水流有效切应力。备注与提示 (Remarks and hints)而言，软底床仅指最新输运落淤的床层。其它所有的床层通常被认为是硬底床。2.4.3 床层的密度 (Density of bed layers)必须为每一层的床层设定密度。整个层（考虑所有粒径组）的密度都必须被设定。概述 (General Description)不同的泥沙类型根据其先前的地质历史而有不同的密度，化学属性，有机物含量和一些其它因素。如果研究区域的范围很大，就有必要获取该范围内不同区域的泥沙，用来生成该范围内每层的样本密度图。可

36、以通过对该范围进行实地测量或者调查其地质历史来获取。垂向上，密度随压缩程度和泥的种类变化。不同的地质时期会留下不同密度的土质。例如，被冰川覆盖的区域可能有非常坚硬的床层，而经过长期沉积的区域可能被相对松散的泥土所覆盖。所以，有必要评估在不同深度的海度。的密度和强度，以此确定海床的垂向度和每层的密河床密度按照以下表达式定义为干密度：Mass of grainsr =dVolume of mixture参数取值 (Data)18MIKE BY DHI 2008河床密度可按以下方式设定· 在整个模型范围内为常数· 在整个模型范围内为变化量河床密度被定义为干密度，其是 kg/m3。

37、如果是第二种情况，则必须准备一空间上至少完整模拟区域范围且包含相关参数的 dfs2 或 dfsu 文件。当采用 dfsu 文件时，相应参与计算网格点上的数据采用分段常数插值方法生成，而当采用dfs2 文件时，则采线性插值。推荐值 (Recommended values)河床的密度范围可以按照下表取值。表 2.3 河床密度备注与提示 (Remarks and hints)如果仅知道容重（或湿容重），可按照下面的公式计算密度：= rs (rb - r )r(r - r )ds其中， rb 是容重(kg/m3)， r s 是颗粒密度(kg/m3)， r s 是水的密度(kg/m3)。2.4.4 床面

38、糙率 (Bed roughness)床面糙率是对水流的阻力。它用于计算河床剪切力。床面糙率决定于床面的形状（，波纹等）以及颗粒的粒径。尽管、波纹等有其动态性，床面糙率是常数（时间上），这是因为局地床面形状的平均变化被认为是常数（时间上）。床面糙率相对于其它河床参数而言是19MIKE BY DHI 2008的。参数取值 (Data)床面参数可设定为：· 在整个模型范围内为常数· 在整个模型范围内为变化量床面糙率可定义为糙率(kn)，其是 m。如果是第二种情况，则必须准备一空间上至少完整模拟区域范围且包含相关参数的 dfs2 或 dfsu 文件。当采用 dfsu 文件时，相应参

39、与计算网格点上的数据采用分段常数插值方法生成，而当采用dfs2 文件时，则采线性插值。推荐值(Recommended values)kn 通常被定义为泥沙直径的 2.5 倍。的推荐值在 0.001m 左右。，对于细颗粒泥沙，床面形状变为主导因素，kn2.4.5 床层间的转换 (Transition between layers)床层间的固结被视为床层间的泥沙转换率，Teisson(1992)。如果考虑床层间的泥沙转换，则必须设定与转换率相关的率定参数。参数取值 (Data)转换可设定为：· 在整个模型范围内为常数· 在整个模型范围内为变化量转换的是 kg/m2/s。2.5

40、地貌 (Morphology)如果研究范围内的地貌变化在某些区域与水深相当，则有必要考虑地貌对力的作用。比较典型的情况是，在浅水区域需要考虑这种长期作用，或者在浅水区域的疏浚/抛泥。2.5.1 概述 (General description)地貌变化是通过每一步的泥沙净通量对水深的更新来体现的。这样能确保河床演变的新，破坏力模拟。Z n+1 = Z n + Dzn20MIKE BY DHI 2008其中，Z n当前时刻的水深值Z n+1的水深值Dz当前时刻的泥沙净通量n时间步长地貌的更新将按以下方式地貌演变：Z n+1 = Z n + Dzn × Speedup其中， speed -

41、 up 是一个无因次因子。每层的厚度以相同的方式更新。要注意的是悬沙量不受此影响，仅仅对河床有影响。2.5.2 备注与提示 (Remarks and hints)如果设定的性。因子过大，则会在河床更新过程中产生内波，从而破坏力模拟的2.6 外力 (Forcings)在这里被考虑的外力只有波浪。2.6.1 概述 (General description)波浪能增大底层的剪切力，从而增大侵蚀率和悬沙浓度。采用结合算剪切力。的剪切力公式来计2.6.2 波浪 (Waves)可以在计算中定义波浪：· 无波浪· 波浪场波浪场 (Wave field)可按以下方式设置波浪：·&

42、#183;·恒定波浪随时空变化的波浪波浪数据库（仅当模拟中考虑可用）21MIKE BY DHI 2008用户必须为波浪设定最小水深在所有的波浪计算中都需考虑液化(Liquefaction)作用。恒定波浪 (Constant waves)如果选择恒定波浪，波浪将是无定向传播的正弦波。用户必须为其设定正北方向的夹角。、波周期及与随着时空变化的波浪 (Time and space varying waves)如果在模拟区域内是变化量，则必须准备一空间上至少完整模拟区域范围且包含相关参数（平均、波周期以及与正北方向的夹角）的 dfs2 或 dfsu 文件。当采用 dfsu 文件时，相应参与计

43、算网格点上的数据采用分段常数插值生成，而当采用 dfs2 文件时，则采线性插值。如果与模拟的时间步长不相等，则采用插值，参见图 2.5。图 2.5 波浪场的插值（时间）22MIKE BY DHI 2008波浪数据库 (Wave database)如果在模拟中考虑该选项才可用。风应力参数在 HD 的框中设置。在风浪主导的区域，以及波浪对水位、风速、风向响应较快的区域，可以采用波浪数据 库。波浪数据库由波浪场的离散值组成，其值是通过水位、风速、风向模拟得到。尽量使数据库中波浪场的离散量减至最少。例如，如果相对于水深而言的水位变化较小，则水位就不需要被离散。或者如果风主要从同一个方向吹来，则风向离散

44、值就可减少。插值 (Interpolation)对每个网格点，由水位、风速和风向的局地值得到波浪参数（显著(Hs)，过零点波周期(T02)以及平均波向( Q )，其方式是通过波浪数据库对波浪场进行插值，参见图 2.6。图 2.6 用波浪数据库对波浪场进行插值23MIKE BY DHI 2008液化 (Liquefaction)液化是对泥沙的弱化，这主要是由波浪造成的孔隙水絮凝所引起。如果要考虑液化，则必须设定一个液化因子。而言，考虑液化会增大侵蚀。最小水深 (Minimum water depth)用户还必须为波浪设定最小水深。 小于该水深就会采用剪切力的流速解而非解。2.7 扩散 (Disp

45、ersion)3D模型中的扩散项主要因不确定的物理过程引起的输移问题。区分平流扩散和垂向扩散问题对于海岸区域而言是十分必要的，平流扩散现象主要由不确定的涡流引起，而垂向扩散现象则主要由近底紊流引起。因此在本模块中亦将水平和垂向扩散现象区分对待。对每一粒径组分别定义扩散。2.7.1 平流扩散 (Horizontal dispersion)平流扩散问题可以采用三种方式表达：···无扩散 (No dispersion)扩散系数公式(Dispersion coefficient formulation)涡粘系数类比公式 (Scaled eddy viscosity for

46、mulation)当采用扩散系数公式时，用户需给定扩散系数的数值（m2/s）。当采用涡粘系数类比公式时，水流方程求解时所需的扩散系数通过将涡粘系数乘以一比例系数的方式获得。关于涡粘系数的详述可参考(Eddy Viscosity)。力模块(Hydrodynamic module)中的涡粘系数参数取值 (Data)在扩散系数的选项中，扩散系数(m2/s)的形式可设定为：· 常数（时间和空间）· 在模型范围内为变化量对于水平空间上变化但垂向上不变的情况，用户需准备一个和模拟区域范围相同的包含有扩散系数的dfsu文件或dfs2文件。当采用dfsu文件时，相应参与计算网格点上的扩散系

47、数采用分段常数插值生成，而当采用dfs2文件时，则采线性插值。比例系数 (Scaling factor)24MIKE BY DHI 2008比例系数可设为一常数。2.7.2 垂向扩散 (Vertical dispersion)垂向扩散问题的相关设置与平流扩散问题相关设置完全相同，可参考2.7.1节。2.7.3 推荐设置 (Recommended values)当基于涡粘系数类比公式（如Smagorinsky或k-模型）的需设定涡粘系数比例因子参数。进行扩散系数设置时，用户这里提到的比例以由数的倒数得出。但为保证该参数取值与k-湍流模型中的经验常数相一致，在计算比例因子时所采用的数的数值大小须与

48、湍流模块设置框中的数设定值相同。数的默认设定值为0.9，相应比例因子为1.1。通常在模拟计算过程中，扩散系数是输运模块调试时非常重要的率定参数之一。准确的直接预估出该参数的合理取值是较为的。但由相似准则可知，扩散系数可以表达为一长度变量和速度变量的函数。在浅水条件下通常可用水深代表上述长度变量，相应速度变量亦可采用一特征流速表征。通常“比例因子”参数的取值接近于1。详细参见Rodi (1980) 。2.8 源 (Source)在许多应用中体现出稀释和悬浮项的点源的重要性，例如，河流出口的泥沙，冷却水或海水淡化厂的进排水口。在泥模块中可以在每个源点设定其每一成分的浓度。源的数目、它们的通用名、位

49、置及强度可以在力模块(HYDRODYNAMIC MODULE)中的源的框(Source dialog)中设定。在选择属性的页面中，可以看到源的地理位置或是源的表。可分别为每个源的每个成分定义浓度。在表点击“Go to.”按键就可以到设定源的对话框。2.8.1 源的设定 (Source specification)源的类型可以两种方式设定：· 指定浓度(Specified concentration)· 额外浓度(Excess concentration)源通量通过Qsource × Csource 得出，其中Qsource 源强， Csource 是源的分量浓度。

50、在力模块(HYDRODYNAMIC MODULE) 的源框 (Source dialog) 中定义源强。25MIKE BY DHI 2008选择指定浓度选项，如果源强是正的（源流进周围水体），则源浓度就是指定的浓度。如果源强是负的（水体流向源），则源浓度就是源点处的浓度。该选项适用于河流出口或者浓度与周围水体无关的其它源。选择额外浓度选项，如果源强是正的（源流进周围水体），则源浓度就是额外浓度与模型中该点的浓度之和。如果是一个孤立的源，该点就是源的地理位置。如果是一个相关的源，该点在水体流出的位置。如果源强是负的（ 度。水体流向源），则源浓度就是源点的浓参数取值 (Data)源项的资料可以被设

51、定为· 不随时间变化的常数· 随时间变化如果要使用随时间变化的源汇项资料，必须在搭建力模块之前，准备一个源浓度(kg/m3)的文件。必须准备一个和模型范围相同的 dfs0 文件。输入文件必须包含整个模型周期。但是其时间步数不需要和自动进行线性内插。力模型的时间步相吻合。如果时间步数不相吻合，模型会2.8.2 备注与提示 (Remarks and hints)点源会进入整个单元格。例如，在初始时刻，点源所在处的入流量被分布于整个单元格。所以模拟结果的源浓度会低于给定的源浓度。2.9 初始条件 (Initial Conditions)初始条件是整个计算范围内初始时刻各组分的浓度

52、在空间上的分布。用户必须给出初始条件。初始条件可以是上一个模拟的结果（即采用热启动方式进行计算）。必须为每个分量分别设定初始条件。2.9.1 组分浓度 (Fraction concentration)用户必须为各个组分设定初始浓度。参数取值 (Data)各组分的初始浓度可设定为：· 在模型范围内是常数· 在模型范围内是变化量是kg/m3。典型的背景浓度为0.01 kg/m3。26MIKE BY DHI 2008如果在模型范围内是变化量，用户在进行力模拟前，需准备一空间上至少完整模拟区域范围且包含相关初始条件的3d无结构数据文件（dfsu）或3d网格化数据文件（dfs3）。当

53、采用dfsu文件时，相应计算网格上的初始条件采用分段常数插值生成；当采用dfs3文件时，相应计算网格上的初始条件基于双线性插值生成；如果输入数据文件中只包含一个时间步的数据时，模型将自动将该时间步的数据作为初值场使用，但如果输入数据文件（如先前模型计算的结果）中多个时间步数据且包含模型计算的开始时刻时，程序将自动将选择模型计算开始时刻的数据作为初值场。2.9.2 床层厚度 (Layer thickness)必须为每层河床设定初始厚度。参数取值 (Data)每层河床的初始厚度可设定为：· 在模型范围内为常数· 在模型范围内为变化量是m。如果在模型范围内是变化量，用户在进行力模

54、拟前，需准备一空间上至少完整模拟区域范围且包含相关初始条件的3d无结构数据文件（dfsu）或3d网格化数据文件（dfs3）。当采用dfsu文件时，相应计算网格上的初始条件采用分段常数插值生成；当采用dfs3文件时，相应计算网格上的初始条件基于双线性插值生成；如果输入数据文件中只包含一个时间步的数据时，模型将自动将该时间步的数据作为初值场使用，但如果输入数据文件（如先前模型计算的结果）中多个时间步数据且包含模型计算的开始时刻时，程序将自动将选择模型计算开始时刻的数据作为初值场。2.9.3 粒径组的分布 (Fraction distribution)用户必须为给定层的给定粒径设定含量(%)。这就是说在每层，每个粒径组的总量必100%。每一层的分布是一个常数。2.10 边界条件 (Boundary Conditions)搭建模型时，模型配置编辑器(set-up editor) 首先会检查网格文件，自动找出网格文件中的边界，并给每个边界预先指定一个名称，用户可自行在Domain 行修改（边界名称，P30）。框中对各个边界的名称进用户在设置边界条件时，方便的基于地图化或列表的方式显示相关设置内容。，不同边界的设置是相互单击“go to”按钮进行相关设置。的，用户可