论坛首页    职业区    学术与写作    工程技术区    软件区    资料区    商务合作区    社区办公室
 您好! 欢迎 登录注册 最新帖子 邀请注册 活动聚焦 统计排行 社区服务 帮助
 
  • 帖子
  • 日志
  • 用户
  • 版块
  • 群组
帖子
  • 4207阅读
  • 13回复

[实例分析]基于三维力学模型的大范围自然边坡稳定性概率评价方法 [复制链接]

上一主题 下一主题
离线ivory
 

发帖
419
土币
1427
威望
804
原创币
0
只看楼主 倒序阅读 使用道具 楼主  发表于: 2009-01-28
关键词: 滑坡
对大范围自然边坡的滑坡危险度进行评价,利用GIS 技术作为平台,把传统的三维极限平衡力学模型耦合到概率统计分析框架中,提出一个基于边坡单元的定量评价的新方法。(岩石力学与工程学报2008年11月)

文献 基于三维力学模型的大范围自然边坡稳定性概率评价方法.pdf (1123 K) 下载次数:41







宁静致远
离线sjzdh01

发帖
7445
土币
6525
威望
39615
原创币
0
只看该作者 1楼 发表于: 2009-01-29
期待学习学习,谢谢提供分享
离线renyun

发帖
267
土币
18621
威望
1722
原创币
0
只看该作者 2楼 发表于: 2009-02-03
期待多提供这方面的资料供分享.
离线hengshan168

发帖
3617
土币
3206
威望
6709
原创币
0
只看该作者 3楼 发表于: 2009-02-10
    
离线hengshan168

发帖
3617
土币
3206
威望
6709
原创币
0
只看该作者 4楼 发表于: 2009-02-10
    
离线hengshan168

发帖
3617
土币
3206
威望
6709
原创币
0
只看该作者 5楼 发表于: 2009-02-10
非常感谢楼主
离线gengdj123

发帖
889
土币
4326
威望
3161
原创币
0
只看该作者 6楼 发表于: 2009-02-11
2.6模块II:运行Visual MODFLOW
下面我们将学习用于Visual MODFLOW的MODFLOW、MODPATH和MT3D的某些运行选项的选择。
2.6.1第10节:流动模拟运行选项
如果你已经完成了模块I中的所有小节或模块II中的其它小节,你可以跳过框中的步骤。
+    File
+    Open
出现一个文件选择窗口。
++   Tutorial文件夹
++   airport.vmf
打开文件,出现Visual MODFLOW主菜单。
从主菜单上选择:
+    Run
你就会进入运行选项屏幕,此时会有一个窗口提示你选择稳定(Steady-state)还是渐变(Transient)模拟。
+    ¤Steady state
+    OK
在Visual MODFLOW中,运行选项分为两个独立部分:流动模拟的运行选项和MT3D模拟的运行选项。当进入运行选项屏幕时,Visual MODFLOW在默认条件下将进入流动模拟的运行选项。从顶部菜单上选择
+    Basic
+    Initial Heads
新模型的默认初始条件是¤Constant by layer(按层给定定水头条件)。Visual MODFLOW将根据每层的标高和边界条件为各层确定一个初始估计值。对于简单的问题而言,这一般都会得到一个收敛的解。然而,一个好的估计值往往可以显著地减少找到收敛值所需的迭代次数。注意,用于输入初始水头估计值的文件可以有多种格式可供选择。因为我们要做的是个简单问题,所以用¤Constant by layer做初始水头估计就已经足够了。
+    OK
此时就会出现一个初始水头条件(Initial Heads)窗口,上面列出了模型中每一层所用的初始水头值。它们应该都是一样的。
+    OK
现在,检查一下Visual MODFLOW所提供的数值求解器的选项。
+    Solver
默认求解器是WHS Solver,这是由Waterloo Hydrogeologic Inc.开发并拥有的求解器。这是目前最快、最稳定的MODFLOW求解器。在本例中,用WHS Solver来求流动解。
+    OK接受
此时会出现一个WHS求解器参数设置(WHS Solver Parameters)窗口,其中列出了所有默认的求解器设置。
+    OK接受
还有其它一些流动模拟的运行选项,如果你有时间的话,可以看看这些本教程中未提到的选项。
2.6.2第11节:MT3D v.1.5和MT3D96模拟的运行选项
如果你已经完成了模块I中的所有小节或模块II中的其它小节,你可以跳过框中的步骤。
+    File
+    Open
出现一个文件选择窗口。
++   Tutorial文件夹
++   airport.vmf
打开文件,出现Visual MODFLOW主菜单。
+     Run
+     ¤Steady state
+     OK
从顶部菜单上选择:
+    MT3D
进入MT3D运行选项屏幕。再从顶部菜单上选择:
+    Advection
此时就会出现一个如下图所示的对流参数设置(Advection Parameters)窗口。
+    ¤Method of Characteristics(MOC)
特征值法(MOC)和改进的特征值法(MMOC)都能求得污染物运移的稳定解,但是MMOC有很大的数值弥散现象。在模块III:输出的可视化中我们将对这两种方法作一个比较。
+    OK接受其余的默认参数
此时将返回到MT3D运行选项屏幕。
+    Output/Time Steps
出现一个输出和时间步控制(Output and Time Steps Control)窗口,如图所示。
MT3D输出控制的默认设置是仅保存模拟结束时的一组浓度值。此处我们要作一些修改。首先,单击þSave simulation results at the end of simulation only(仅模拟结束时保存一次模拟结果)的前面的框,去掉ü,即不选该项,然后填入并核对下列各项:
Simulation Time:    7300 daysÖ
单击标明Max # of transport steps(最大运移步数):的框
type:    10000Ö
+    ¤ specified times
在¤ specified times(指定时间)标签下的框内输入下列数值:
1825¿
3650¿
7300
+    OK接受这些参数
以上就是MT3D的运行选项部分的教程。需要再次说明的是,还有许多本教程未提及的运行选项,如果你有时间的话,可以看看这些选项。
虽然本教程可以让模型进入运行,但可能由于时间的关系你还不想就在此时运行。下面将向你显示模型是如何开始运行的。
从上部菜单栏中选择Translate/Run,将出现运行设置(Run Setting)窗口,其中有以下一些可设置的项目:
选择恰当的用于模拟的数值引擎,Visual MODFLOW 将把程序输入部分提供的信息转换成带有合适扩展名(.BAS, .WEL, .BCF等)的标准MODFLOW、MODPATH和MT3D输入文件。然后Visual MODFLOW将自动地执行模型计算分析。模拟完成,得到一组收敛解后,它就会返回Visual MODFLOW的主菜单。注意:运行Visual MODFLOW的MT3D1.10时,运行时间可能很长。对于不同的计算机运行所花时间大致如下:
                              MT3D 1.5               MT3D 96
l    486                      10-12小时                3-4小时  
l    奔腾100                  5-6小时                 1.5-2.5小时
l    奔腾 150                  3-4小时                 1-2小时    
+    OK      开始运行
如果不想运行该模型而是继续我们的教程,就从底部菜单栏中选择F10-Main Menu按钮。
如果你选择了运行模型,Visual MODFLOW将启动Win32 MODFLOW程序组。在本手册第四章:运行选项中你可找到有关Win32 MODFLOW程序组的信息。
2.7模块III:输出的可视化
如果你已运行了模型,你可以将你的模拟结果可视化显示,如果你没有运行模型,你还可以将包含在Visual MODFLOW安装盘中的AIRPORT模型模拟结果可视化显示。
如果你已经完成了模块I和II并且运行了模型,你可以跳过框中的步骤。
从主菜单:
+    File
+    Open
出现一个包含一系列Visual MODFLOW输入文件的文件选择窗口。
++   Tutorial文件夹
++   output.vmf
+     OK
+     Output
此时转到Visual MODFLOW的输出屏幕上。在默认情况下,输出屏幕将显示不带颜色填充的等势线图,如下图所示。
从顶部菜单上选择Contours就可以看看程序能绘制哪些等值线。其下拉菜单会列出以下一些可以选择的等值线。
选择Head Equipotentials回到等势线图上。
2.7.1第12节:等势线和等值线绘制选项
如果你已经完成了模块I和II并且运行了模型,你可以跳过框中的步骤。
从主菜单上:
+    File
+    Open
出现一个包含一系列Visual MODFLOW输入文件的文件选择窗口。
++   Tutorial文件夹
++   output.vmf
+    OK
+    Output
要选择等水头线的等值线绘制选项:
+    Options(在左边的菜单条上)
就会出现一个如图所示的等水头线绘制选项(Equipotential Overlay Contouring Options)窗口,但没有右边的颜色充填条。要激活颜色充填:
+    þ Colour shading
在标明间隔(Interval:)的文本框中,将值从0.5改为0.25。在标着标注(Labels)的框中,将小数位数改为2。
你也可以通过单击标着等值线生成精度/速度(Contouring resolution/speed)的按钮来改变等值线的绘制速度。每点一次,等值线的绘制速度增大2倍。在本次练习中,按钮上应该是High/Slow。
+    OK接受
等值线的位置应该与前一个图非常相似,但由于等值线的精度降低了,这些线会有一些轻微的摆动。上图没有加入颜色充填。
其它没有通过菜单激活的等值线选项可通过在模拟区内按下鼠标右键来激活。此时会有一个带有增添、删除、移动等值线和标注的选项的弹出窗口(见下图)。
选择Add Contour选项,将鼠标移至模拟区的任意一处并单击鼠标左键。在你鼠标点下的位置就会加入一条等值线。要再加入一条,只需再点一下鼠标即可。要再次弹出等值线选项窗口,只需简单地再按一下鼠标右键即可。这一次我们选择Move Label,然后将鼠标移至你想要移动的等值线标注的位置上。在标注上按下鼠标不放,然后将其沿等值线拖至想放的位置,放开鼠标,这样就设定了标注的新位置。标注都移动完后,再按一下鼠标右键。
下面,让我们看看场地的剖面情况。从左边菜单上选择View Column,然后将鼠标移至模拟区内。在模拟区中间位置按下鼠标左键选择一列。
要去掉用户设置上去的等值线,在剖面上单击鼠标右键弹出等值线绘制选项窗口。
+    Delete all custom
+    View Layer回到平面模式
在进入下一节前,必须关掉色彩充填功能。在左侧菜单条上
+    Options
+    Colour shading
+    OK
2.7.2第13节:流速矢量及其选项
如果你已经完成了模块I和II并且运行了模型,或者来自模块III的某一节,那么你可以跳过框中的步骤。
+    File
+    Open
出现一个包含一系列Visual MODFLOW输入文件的文件选择窗口。
++   Tutorial文件夹
++   output.vmf
+    OK
+    Output
要察看流速矢量选择:
+    Velocities
你将进入流速矢量输出选项屏幕,见下图。流速矢量将根据其默认的设置进行绘制,即矢量相对大小与流速的大小一致。
要让流速矢量与大小无关
+    Direction(左边菜单上)
这样所有的箭头的大小都将是一样的,它们只代表了流动方向。
如果想改变显示的矢量箭头数目
+    Options(左边菜单上)
并输入:
Vectors:    40
这就指定了一行上的流速矢量的数目。
þAutoscale
¤Variable Scale
+    OK接受
注意颜色:在默认情况下,红色表示朝外(也就是说,从层的角度看是向上的)的流速,蓝色表示向里(也就是说,从层的角度看是向下的),绿色表示与平面平行(也就是说,从层的角度看是水平的)。
再让我们看看场地的剖面情况。从左边菜单上选择View Column,然后将鼠标移至模拟区内。单击鼠标左键,选择任意一列。就会出现一个与下图类似的在剖面上同时显示等势线和流速矢量的图。
为显示水力传导系数图层的离散化情况,从底部菜单栏中选择F9-Overlay按钮,将出现一图层控制(Overlay Control)窗口,该窗口按字母顺序列出所有可以打开或关闭的图层。把鼠标置于Conductivity Overlay(水力传导系数图层)之上并双击鼠标左健。这样就会在Conductivity Overlay旁出现一个星号(*),意味着它被激活了。
+    OK 显示水力传导系数图层剖面。
从左边菜单选择View Layer,单击模型横剖面中的第一层,返回到模型平面视图。
要从屏幕中消除流速矢量,从底部菜单栏中选择F9-Overlay按钮,出现图层控制窗口后,把鼠标置于Velocity Overlay上并单击使之高亮度显示,再单击标有ON的按钮使之变为OFF。这样Velocity Overlay旁的星号将会消失,意味着它已无效了。
+    OK  显示没有流速矢量的窗口
2.7.3第14节:流线及其选项
如果你已经完成了模块I和II并且运行了模型,或者来自模块III的某一节,那么你可以跳过框中的步骤。
+    File
+    Open
出现一个包含一系列Visual MODFLOW输入文件的文件选择窗口。
++   Tutorial文件夹
++   output.vmf
+    OK
+    Output
+    Pathlines
默认条件下的屏幕将显示所有流线在模型当前层上的投影(如下图所示)。
要显示当前层的流线,从左边菜单栏中选择Segments按钮 。这样图上就会只显示模型当前活动层的流线。
用Next按钮可以顺序向下显示模型的每一层来查看质点流线的位置。
+    Projections   观察所有流线。
注意流线上标有箭头表明流向。这些箭头也作为时间标志表示一质点到达某一目的地的时间长度。要确定每一时间标志的时间间隔,从左边菜单栏中选择Options,此时将显示一流线选项(Pathlines Options)窗口,如图所示。
注意:þUse time marks选项是有效的。在标有¤Regular every ‘1000’ days方框中显示了标志时间间隔。
为了显示在10,000天内流线能走多远;
+    ¤ Time Related 并在方框中输入10,000。
+    OK 显示10,000天内与时间相关的流线
现在显示横剖面上的流线。
+    View Column
在屏幕上移动光标,使高亮度条带位于不连续弱含水层区附近,并单击它。
现在我们返回到平面视图显示。
+    View Layer
在屏幕上移动光标使高亮度条带位于第一层,并单击它。
要从图中删除流线,可从左边菜单栏选择F9-Overlay按钮,并向下滚动图层列表直到你看到Particle Overlay和Pathline Overlay为止。在它们上面双击鼠标去掉星号,使这两个图层失效。
+    OK
2.7.4第15节:MT3D浓度等值线
如果你已经完成了模块I和II并且运行了模型,或者来自模块III的某一节,那么你可以跳过框中的步骤。
从主菜单:
+    File
+    Open
出现一个包含一系列Visual MODFLOW输入文件的文件选择窗口。
++   Tutorial文件夹
++   output.vmf
+    OK
+    Output
在这一节中,我们将对特征值法和改进的特征值法进行比较。
+    Contours
+    MT3D Concentration
此时将会出现浓度等值线图,该图显示了第一个输出时间1825天时的浓度情况。因为用对数坐标及颜色充填来表现浓度等值线会更简单一些,因而你可以按如下步骤来自行设计等值线图:
+    Options
¨ Automatic contour levels
þ Use cudtom contour levels  
þ Colour shading
在方框中输入如下用户自定义的等值线间隔:
0.0
100.0
500.0
1000.0
2000.0
4000.0
+    OK
尽管这不是一个严格的对数坐标,但它还是能提供一种更好的浓度等值线的表示法。
下图是用特征值法(MOC)和改进的特征值法(MMOC)所计算出的在1825天时的浓度等值线。
上图是用改进的特征值法(MMOC)所计算出的在1825天后的解。
上图是用特征值法(MOC)所计算出的在1825天后的污染羽。
用MMOC来计算的污染物运移模拟的数值弥散性比MOC要高。尽管这两种模拟在数值上都是稳定的,但我们仍建议用户利用自己的专业知识来判断到底多大的数值弥散是允许的。
这些MT3D模拟指定了三个输出时间:1825天、3650天和7300天。从左边菜单栏选择Time可以查看输出时间的序列。
+    7300
+    OK   绘制7300天以后的浓度等值线图。
下面两图分别显示了用MMOC和MOC所做的第7300天的污染羽分布。
上图显示了用MMOC所做的第7300天的污染羽。
上图显示了用MOC所做的第7300天的污染羽。
在7300天以后,MMOC模拟的数值弥散量比第1825天要大得多。而MOC模拟则表明污染物向不连续弱含水层中的运移情况。
对单个或多个观测点都可以画出其浓度穿透曲线。从左侧菜单条上
+    Select来选择观测井
会有一个弹出窗口提示你选一种选择观测井的方法(Single、Box和Polygon)
+    Box
沿三个观测井用鼠标拉一个框后,这些井就会变为红色,表示它们已被选中。
从左边菜单上
+    Conc. Vs Time
就会显示一个包括每个井的浓度对时间关系图。观测到的浓度在图上以星号表示。下图是MOC模拟的浓度-时间图。穿透曲线的不规则是MOC技术的离散性的结果。对于高浓度,稳定状态的浓度可以通过时间上的平均来求得近似值。而在低浓度的情况下,由于噪声的影响,很难确定一个平均浓度。
作为对比,下图是MMOC模拟的浓度-时间图。其曲线要光滑一些是因为MMOC中增加了数值弥散的结果。下面的稳定浓度也显示出作为使用MMOC的结果,出现了数值弥散。
+    Close
整个教程到此为止就结束了。正如前面所说,这些说明的目的是给你提供一个Visual MODFLOW的主要特怔和功能的介绍。如果你有时间的话,我们鼓励你回过头去试试那些使用说明中还未涉及的特怔和功能。

第三章    输入
本章主要内容如下:
¨    建立一个新的Visual MODFLOW模型
¨    向Visual MODFLOW输入已有的MODFLOW模型
¨    编辑已有的Visual MODFLOW模型
如果你以前没有用过Visual MODFLOW,又想尽快地熟悉Visual MODFLOW的用户界面,我们建议你先看完第一章,再跟着做第二章的练习。
Visual MODFLOW的数据有两种方法生成。我们推荐使用Visual MODFLOW的图形环境,我们将在本节中描述有关内容。但是你也可以文本编辑器创建和修改数据文件。有时,在你对模型已经比较熟悉,并且只需对现有的文件作一点小小的修改时,用文本编辑器可能会更快一些。为方便起见,所有Visual MODFLOW数据文件的结构和形式都有将在第五章中加以描述。然而,如果你使用文本编辑器,Visual MODFLOW并不会去检查数据的一致性。Visual MODFLOW仅监测和检查在图形环境下输入的数据的一致性和合理性。因此,在用文本编辑器修改数据文件时,我们建议你要特别小心。
3.1新模型的建立
如果你正要创建一个新的Visual MODFLOW模型,请继续做下去。但是,如果你要用一个已经存在的模型,请跳过这一节到基本输入那一节。
在Visual MODFLOW 中,要创建一个新的模型的话就在File菜单中选择New即可。输入一个新的文件名。如果操作不成功的话将显示一个错误信息(比如已存在与此相同的文件名)。
然后将显示如下对话框让你去
l    设置网格的默认大小,
l    设置底图(.BMP或.DXF文件),以及
l    选择模型的运算单位。
在对话框中,你必需输入行、列和层数,以及模拟区最大、最小X,Y,Z坐标值。如果你已选择使用.DXF文件的坐标值,则该文件中的最大、最小X和Y坐标值将作为缺省的最大、最小X和Y坐标值。总层数和模型的最大、最小标高也应在此输入。Visual MODFLOW将在此基础上生成一个均匀的网格。而后,象在网格(Grid)一节所描述的那样可以细化网格,但是它的总体大小不变。
注意:运算单位一旦确定就不能再加以修改。
你必须为你的数据选择单位。单位不必和国际单位或英制单位一致;允许使用比如水力传导系数的单位用厘米/秒而抽水量用克/分之类的单位。在你输入运算单位后,选择Create。
3.1.1底图
Visual MODFLOW程序的所有模块都支持使用底图。这些底图既可以是DXF(AutoCAD v14)也可以是BMP格式。在选择模拟区域(Select Model Region)窗口中,对于DXF和BMP文件都会出现如下图标。
       放大:单击你想缩大区域的一角,再单击该区的另一角。
       缩小:缩至显出整个模型。
       新图:弹出一个选择模型底图(Select Map for Model)窗口,你可在此指定一个新的DXF或BMP文件作为底图。
       重设区域大小:允许用户在区域窗口的一角处单击鼠标并拖动它到所需的大                小。双击该区将在Rotate和Resize Region之间切换。
       旋转区域:允许用户在区域窗口的一角处单击鼠标并拖动它到所需的角度。双击该区将在Rotate和Resize Region之间切换。
       对齐区域:按零度角排列区域。
       扩充区域:把区域扩充到整个dxf图。
       把鼠标转化为一般的光标。
3.1.1.1 DXF格式的底图
使用DXF文件作为底图时,你既可以用鼠标也可以通过在文本框中键入值方法来重新设置区域大小。图象的坐标直接从DXF文件中读取,并显示在世界坐标系(World Co-ordinates)的文本框中。你也可以输入和世界坐标不同的局部坐标,它要在网格坐标系(Grid Co-ordinates)的文本框中输入。最后一点,允许用户改变模型的行列数并显示网格。
输入一般模型创建信息后,你将进入输入窗口。
3.1.1.2 BMP格式的底图
使用BMP文件作为底图时,因为文件中没有包含坐标信息,所以这种BMP文件必须具有地理参考系。将会出现如下窗口提示用户如何对图形作地理参考系统。
在你设置模拟区域以前,BMP图形应该有两个具有已知世界坐标的参考点。在图形上单击鼠标,然后在右边活动的文本框中输入坐标设置地理参考点。重复一次输入另一个地理参考点。一旦这两个坐标被接受,其余的文本框就会激活,你也就可以设置网格坐标、旋转角度、行列数,并显示网格。
在建立图形的地理参考点时,会用到以下图标:
      设置地理参考点
      删除地理参考点
使用BMP作为底图时,你既可以用鼠标也可以通过在文本框中键入值的方法来重新设置区域的大小。图象的坐标直接从BMP文件中读取,并显示在世界坐标系的文本框中。你也可以输入和世界坐标不同的局部坐标,它要在网格坐标系的文本框中输入。最后一点,允许用户改变模型的行列数并显示网格。
输入一般模型创建信息后,你将进入输入窗口。
3.2载入原有的模型
在主屏幕的顶部菜单中选择File、Import后,你就可以载入已经存在的MODFLOW或FLOWPATH模型。
3.2.1 MODFLOW
在输入MODFLOW数据时,Visual MODFLOW将尽力从MODFLOW文件的数据中重建你的数据集。然而,由于数据或Visual MODFLOW中的某些特定的设置,下面所说的顺序可能会有所不同。你必须核对你的数据以确保其中所有信息都是正确的。
下面列出了一些目前已知的在载入MODFLOW数据时会存在的限制。
l    层间越流(VCONT):Visual MODFLOW要求输入每一层的明确的层厚和水力传导系数。然后,Visual MODFLOW在转换时生成层间越流项。这将会给那些现有的具有不同越流项的隐含弱透水层的MODFLOW模型带来严重的问题,因为它们的弱透水层的厚度和水力传导系数已合并在越流项之中了。
l    垂向导水系数:在载入MODFLOW数据文件时,.BCF包中的VCONT参数被省去了。在Visual MODFLOW中垂向和水平向水力传导系数之比设为1:10。这个值可在数据文件载入后进行修改。
l    层底标高:如果Visual MODFLOW在数据文件碰到类型为0或2的层位时,就会提示用户输入层底标高。对于顶层,将会提示用户输入层顶标高。如果层底标高有变化将会出现一些问题。你可以把所有层位都设置为类型3,这样在重新生成MODFLOW 的数据文件时就会使数据的损失降到最小。
l    单位储水系数、单位给水度和空隙度:在输入时每层都只能设置一个单位储水系数、单位给水度和空隙度的值。但以后你可以逐单元地修改这些值。
l    洪水演算包: Visual MODFLOW不支持洪水演算包,但是如果存在是不会忽略的(见下文)。
l    弱透水层压缩性包: Visual MODFLOW不支持它,如果MODFLOW 数据文件存在该项也会被忽略。
用Visual MODFLOW运行洪水演算程序
    虽然Visual MODFLOW支持大多数常用的MODFLOW的程序包,但它目前还不支持用户对洪水演算程序的输入数据进行创建和可视化显示。但Visual MODFLOW可以利用用户所提供的适当的洪水演算的输入数据文件来生成合适的运行MODFLOW模拟所需的MODFLOW.IN和FILENAME.BAS文件。
    第四章介绍了如何在模型模拟中包含洪水演算程序的方法。在Visual MODFLOW对输入数据进行转换时,程序对工作目录进行检查,寻找是否存在洪水演算程序的输入数据文件(filename.STR,其中filename是模型文件的名字)。如果存在,Visual MODFLOW就会用这个数据文件来运行MODFLOW(假设它与模型数据文件的其它部分是兼容的)。这个.STR文件必须包含一个与.BAS文件中给出的模型网格及应力期相兼容的数据输入结构。只要模型网格或应力期有什么改变,都必须手工调整.STR数据文件以反映模型网格或应力期的变化。            
如果在弹出菜单中选择MODFLOW,将会出现一个选择文件的窗口,要求你输入MODFLOW模型的.BAS文件名。程序会自动识别.BAS文件中用于特定模拟的程序包,同时程序将读取如下带有标准扩展名的标准MODFLOW文件并输入数据:
                modelname.BAS
                modelname.BCF
                modelname.DRN
                modelname.EVT/.EVP
                modelname.GHB
                modelname.HFB
                modelname.RCH
                modelname.RIV
                modelname.STR
                modelname.SIP
                modelname.WEL
其中所有文件的modelname必须相同。
为MODFLOW模型选择了一个MODFLOW.BAS后,程序会要求你为这个新建的Visual MODFLOW模型选择一个名字。如果你选择了一个和MODFLOW模型一样的名字,在数据转换时将覆盖掉所有的MODFLOW源文件。
此时,输入工具包(Import Utility)将会要求你输入模型的单位。
这类似于创建一个新模型时的单位选择。缺省时间单位直接从原始的.BAS文件中读取,同时与此相应地读入并保存所有的数据。然而用户可以用上述窗口来修改Visual MODFLOW显示水力传导系数和补给的单位。选择好模型单位后,将出现如下窗口,问你模型的基点是否从.DXF文件中读取。
如果你回答是,将出现一个选择.DXF文件的窗口。选择了.DXF文件后,或你在上面的窗口中回答了否,将出现如下的弹出窗口,问你是确认还是输入模型基点。默认值是.DXF的原点或零。
如果你的数据中少了一些文件,输入工具包将会通知你同时中止输入。在试图输入同一个数据集以前,请找到或创建必需的文件并把它们放在.BAS文件的目录。如果你不想创建这些文件或不想把它们包含在输入中,修改.BAS文件使它不包括这些文件的设备号。
此时,输入工具包将要求你输入每层的顶底板标高。如果层的类型是1或3,就不要在此编辑这些数据。
输入层位标高后,程序会要你输入指定的和计算出来的水力传导系数的有效数字的位数。这个值仅在输入时发生作用,它通过减少水力传导系数计算的位数而减少了输入时间。
下一步将显示水力传导系数数据库,如下图所示。生成的水力传导系数是唯一的,Visual MODFLOW将这些值无重复列出。应仔细核对这些值,如果有必要,修改一下。
最后的窗口要用户输入单位储水度、单位给水度和空隙度的缺省值。这些变量的值没有输入,对所有的层位这些值在开始时都将设置为常量。
3.2.2 FLOWPATH 5.x
如果在弹出菜单中选择了FLOWPATH,将会出现一个文件选择窗口,要求你输入FLOWPATH 5.x模型的名字。程序将读入所有具有标准扩展名的标准FLOWPATH 5.x文件并输入数据。然而,因为数据文件或Visual MODFLOW中的某些特定设置,下列顺序可能会有所不同。你必须核对你的数据文件确保其中所有信息都是正确的。
指定了文件名后,出现的第一个窗口将要求你输入单位储水度和单位给水度的缺省值,因为在FLOWPATH 5.x中不使用这些值。
接着,输入工具包将显示由FLOWPATH 文件确定的水力传导系数数据。应仔细核对这些值,如果有必要,修改一下。
因为FLOWPATH 5.2模拟的是二维地下水流,没有垂向分量。显示的下一个窗口要求输入垂向水力传导系数因子。
最后一个窗口显示由FLOWPATH 5.2文件确定的单位储水度、单位给水度和空隙度数据。应仔细核对这些值,如果有必要,修改一下。
3.3基本输入
如果你要编辑一个已有的数据文件,启动Visual MODFLOW后在主菜单选择File项的Open。从文件选择窗口中选择要用的文件。
打开一个已有模型或如上节所述建立了一个新的模型后,从主菜单中选择Input。程序将会转到Flow输入窗口,如下图所示。
在这个窗口内你可以定义三维有限差分网格、井位、属性、边界条件、流线的起始点、模型校验的观测点、分析模型内流量的分区以及附加的注释。另外,MT3D输入屏幕(从Flow输入窗口中进入)可用来输入初始浓度、运移边界条件、模型校验的实测浓度以及弥散系数,如下所示。
创建一个流动模型需要一个完整的数据集,包括网格和边界形状、边界条件和含水层属性。校验资料、流线的设置、区域均衡计算资料和注释是可选的。另外,MT3D运移模拟需要指定初始条件、一个浓度边界条件和弥散系数分布情况。用于模拟校正的浓度实测值是可选的。
Flow输入窗口的顶部菜单栏包含如下输入选项:
File  Grid  Wells  Properties  Boundaries  Particles  Annotate  Zbud  MT3D  Tools  Help
File    保存文件,打印或返回主菜单
Grid    修改网格(默认屏幕)
Wells    插入、修改或删除抽水井和观测井
Properties    增加或修改水力和储水特性
Boundaries    输入或修改流动边界条件
Particles    插入或修改用于质点追踪的质点
Annotate    注释模型
Zbud    插入或修改区域均衡数据
MT3D    切换到MT3D输入窗口
Tools    调用单元监测器和预测工具
Help    获取一般帮助
MT3D输入窗口的顶部菜单栏包含如下输入按钮:
File  Initial Concentration  Boundaries  Calibrate  Dispersion  Chem.Reaction  Flow  Help
File    保存文件,打印或返回主菜单
Initial Concentration    输入或修改初始浓度分布
Boundaries    输入或修改TRANSPORT边界条件
Calibrate    输入或修改TRANSPORT校正数据
Dispersion    输入或修改弥散数据
Chem.Reaction    设置吸附和放射性衰变系数
Flow    切换到Flow输入窗口
Help    获取一般帮助
3.4文件菜单(File)
当你从输入菜单栏单击File时将出现如下菜单选项:
Save    保存数据集
Print    在当前输出设备上打印屏幕并选择打印设备
Print to DXF file    把屏幕打印到.DXF文件
Main Menu    返回主菜单
3.5网格菜单(Grid)
网格设计窗口允许用户定义和离散化模拟区域。用户可以增加网格线或删除单元或者将某些单元排除在计算之外。
当你在输入窗口中选择Grid选项,就会激活一个‘grid’边菜单栏,该菜单具有以下选项:
Edit Columns    打开一个弹出窗口选择列编辑选项,包括增添、删除、移动、输入或输出列。(仅在行或层显示方式下有效)。
Edit Rows    打开一个弹出窗口选择行编辑选项,包括增添、删除、移动、输入或输出行。(仅在列或层显示方式下有效)。
Edit Layer    打开一个弹出窗口选择层编辑选项,包括增添、删除、移动、输入或输出列。(仅在行或层显示方式下有效)。
Inactive Cells    打开一个下拉式菜单设置模型的无效区域,有以下选项:
Mark Poly. Active    把一个多边形区域设置为流动计算的有效单元。按鼠标右键封闭多边形。
Mark Inactive    把一个多边形区域设置为流动计算的无效单元。按鼠标右键封闭多边形。
Mark Single    按鼠标左键使单元无效,右键使单元有效。
Copy Single    按左键复制单个无效单元到其它层、行或列。
Copy Polygon    复制一个无效单元区域到其它层、行或列。
Import Surface    从ASCII文件或SURFER的‘.grd’文件输入一个变化的表面来定义地表或指定层位的底板。同时提示用户输入层的最小厚度和用于内插的观测点数。
X smoothing    在X方向上平滑网格奇点
Y smoothing    在Y方向上平滑网格奇点
3.5.1网格线编辑(行或列)
允许的最大网格是:X方向上500根网线(499列),Y方向上500根网线(499行),Z方向上60层。因此未知水头值的最大数量是15,000,000。如果这个数量还不能解决你的问题,请与我们联系。
如果你的计算机没有足够的内存,Visual MODFLOW将通过Windows的内存管理使用你的硬盘作虚拟内存。
要编辑网格请从边上的工具条中选择Edit Column或Edit Row。此时会弹出一个带有如下选项的菜单:
Add    以图形方式或通过弹出对话框在特定的位置上增加网线
Delete    删除一个以高亮度显示的网线
Move    移动一个以高亮度显示的网线
Refine by    用所选择的因数细化两网线间的网格(比如细化因数是二,将把每行或列等分成两行或两列)
Coarsen by    用所选的因数粗化两网线间的网格(比如粗化因数是二,将把两行或两列合并成一行或一列)
Import form file    从.txt文件载入网线坐标
Export to file    输出网线坐标到.txt文件
离线fx_wxg

发帖
710
土币
1
威望
3308
原创币
0
只看该作者 7楼 发表于: 2009-02-18
      
附件: 地下工程新奥法_10014199.rar (3260 K) 下载次数:0
好好学习天天向上
离线xwjljh

发帖
531
土币
59
威望
2196
原创币
0
只看该作者 8楼 发表于: 2009-02-18
学习学习   
离线xwjljh

发帖
531
土币
59
威望
2196
原创币
0
只看该作者 9楼 发表于: 2009-02-18
期待学习学习,谢谢提供分享
离线wpy2000p

发帖
699
土币
882
威望
4253
原创币
0
只看该作者 10楼 发表于: 2009-03-09
谢谢楼主,学习了。
离线wq020200203

发帖
281
土币
830
威望
2190
原创币
0
只看该作者 11楼 发表于: 2009-06-24
期待多提供这方面的资料供分享.
快速回复
限100 字节
温馨提示:欢迎交流讨论,请勿纯表情、纯引用!
 
上一个 下一个

      https://beian.mps.gov.cn/ 粤公网安备 44010602012919号 广州半山岩土网络科技有限公司 粤ICP备2024274469号

      工业和信息化部备案管理系统网站