GMS地面沉降SUB模块例题中文翻译(笔记).docx

2简介SUB（the subsidence and aquifer-system compaction package）是由美国地质调查局开发的用来模拟含水层压缩和地面沉降的程序包。Sub程序包模拟包括弹性和非弹性夹层的压缩，也能够模拟夹层的滞后和非滞后排水。2.1概要下面是操作步骤：把sub程序包加入到网格法模拟模型中去创建一个简单的概念模型去说明sub程序模块是怎样被概化并且插入到modflow中去的3问题描述我们将要解决的问题将在图1中说明，模型基于美国地质调查局模型，描述如下：羚羊峡谷（antelope）是莫哈维沙漠（mojave）西部的地形封闭的盆地，距离洛杉矶大约50英里。羚羊峡谷地下水流域面积大约940平方英里，北部被低矮山丘和断层分开，1972年以前，地下水占总供水源的90%多；自1972年以来，供应量占50-90%。包括正在迅速的发展的城市lancaster和palmdale在内，在羚羊峡谷大规模抽取地下水。 本区地下水系统由上层、中层和下层含水层组成。含水层由砾石、砂、淤泥和粘土冲积沉积和粘土及淤泥质粘土湖相沉积，为非固结沉积。研究区之前的地下水开发，补给来自周围山区的径流入渗。地下水从补给区流向沙漠盆地通过蒸发和泉流出露进行排泄，部分地下水流水平向阻隔，例如断层，已经在地下水流区域确定。地下水位下降由于地下水开发已经减少了自然因素的排泄，而农村和城市的抽水已经成为了该区地下水系统的主要排泄因素。回水入渗已经成为了该区地下水系统重要的补给来源。该模型被离散为网格模型，为43行，60列，3层。每层对应相应含水层。模拟期从1915年-1995年80年，并且第一年为稳定状态。4开始如果必要，打开gms。如果gms已经运行，选择File | New命令确保项目设置是默认状态5读取工程首先，我们读取工程：首先File | Open浏览所建工程文件夹打开start.gpr.你会看到图1的modflow模型6用网格法加入sub程序包6.1重命名并保存模型我们将要做一下变更，接下来是重命名并且保存模型1、选择file→save as2、确保一直在工程文件夹中3、对该工程重命名为avgrid.gpr.4、保存6.2激活sub程序包我们需要打开sub程序包1、选择MODFLOW | Global Options打开MODFLOW Global/Basic Package对话框2、选择Packages按钮打开MODFLOW Packages对话框3、在Optional packages.下勾选SUB - Subsidence程序包4、确认两次退出对话框6.3定义非滞后夹层总结：定义夹层的前期固结水头、弹性与非弹性储水系数，其中前期固结水头为含水层中初始水头（也即是那个2d散点）首先我们先为第一和第二模型层加入非滞后夹层选择MODFLOW | Optional Packages | SUB - Subsidence打开MODFLOW SUB Package对话框这个步骤是打开SUB package对话框。对话框左面Options标签包括程序包常用值，滞后夹层区域值在右面。Delay interbed和no-delay interbed 被用来创建夹层和设置夹层相关数据序列选择No Delay Interbeds 标签通过点击两次在No-delay interbed layers表格底部的按钮插入两个新的夹层系统,。为第第一夹层赋值1，第二层夹层赋值2.选择一个夹层系统，其相应的序列值就被显示出来。点击夹层系统1，它的值就显示在了序列编辑界面选择2D data set→ array按钮设置前期固结水头（使用初始水头生成的网格中的最大值和最小值）选择preconsolidated head 数据集然后选择ok按钮退出对话框在View/Edit弹出菜单项选择（sfe）elastic skeletal storage coeff （夹层的弹性储存系数）（本项之所以会出现弹性与非弹性的设置都是遍布整个网格的数值，是因为这些是设定值，意思是说，假如夹层发生了弹性形变就用sfe计算，发生非弹性就用sfv计算）。选择Constant - Array按钮，设置序列数值为1.5e-4在View/Edit弹出菜单项选择（sfv）inelastic skeletal storage coef 选择Constant - Array按钮设置，设置序列数值为8.0e-3同步骤5，点击夹层系统2，它的值就显示在了序列编辑界面仿照步骤5-11进行设置，同样选择Preconsolidated Head 数据集，(Sfe) Elastic skeletal storage coef值为9.0e-5，(Sfv) Inelastic skeletal storage