GMS-全面+經(jīng)典教程

地下水模擬軟件GMS教程翻譯目錄TOC\o"1-2"\h\z\u1GIS模塊 61.1簡介 71.2未安裝ESRIArcObjects 71.3開始 81.4讀取Shapefile文件 81.5查看Shapefile文件 81.6查看屬性表 91.7文件轉(zhuǎn)換為2D離散點(diǎn) 91.8已安裝ESRIArcObjects 111.9啟動(dòng)ArcObjects 111.10讀出數(shù)據(jù) 111.11查看shapefile文件 121.12圖形文件轉(zhuǎn)化為屬性對象 132屬性對象 172.1目標(biāo) 172.1簡介 172.2開始 182.3屬性對象 182.4結(jié)論 243MODFLOW—概念模型法 253.1簡介 263.2問題描述 263.3開始 283.4導(dǎo)入模型底圖 283.5存儲(chǔ)工程文件 283.7定義邊界 293.8創(chuàng)建本地源匯圖層 303.9定義補(bǔ)給區(qū)域 363.10定義水力傳導(dǎo)系數(shù) 383.11LocatingtheGridFrame定位柵格骨架 403.12創(chuàng)建柵格 403.13初始化MODFLOW數(shù)據(jù) 413.14定義活動(dòng)/不活動(dòng)區(qū)域 413.15插值（篡改）高程值 413.16轉(zhuǎn)換概念模型 443.17檢查模擬 453.18運(yùn)行MODFLOW 453.19顯示等水位線 453.20側(cè)邊瀏覽水頭 463.21查閱水流均衡 463.22結(jié)論 474MODFLOW-GridApproach 484.1簡介 494.2問題描述 494.3開始 504.4單位 504.5創(chuàng)建柵格 514.6創(chuàng)建modflow模擬 514.7單元邊界條件直接賦值 544.8LPF模塊 554.9補(bǔ)給（回灌）模塊 574.10排水溝模塊 574.11井模塊 594.12檢查模擬 624.13保存模擬 634.14運(yùn)行MODFLOW 634.15顯示結(jié)果 634.16區(qū)域均衡 644.17結(jié)論 665MODFLOW—層數(shù)據(jù)插值 675.1簡介 685.2開始 685.3層數(shù)據(jù)插值 685.4案例 695.5例1—完整連續(xù)地層 695.6案例2—裂隙 715.7案例3—出露 735.8案例4—基巖消蝕 745.9結(jié)論 766MODPATH 786.1簡介 786.2例題描述 796.3開始 796.4導(dǎo)入工程文件 796.5定義孔隙度 806.6定義初始位置 806.7顯示設(shè)置 816.8粒子示蹤屬性設(shè)置 816.9垃圾填埋場的粒子示蹤 836.10通過區(qū)域代碼設(shè)置顏色 856.11流線長度和時(shí)間計(jì)算 856.12通過區(qū)域代碼設(shè)定截獲區(qū) 866.13小結(jié) 867MT3D-概念模型方法 887.1簡介 887.2案例描述 897.3開始 897.4導(dǎo)入工程 897.5確定單位 897.6初始化MT3DMS模擬 907.7定義含水層特征 927.8入滲濃度的賦值 937.9概念模型轉(zhuǎn)換 937.10層厚設(shè)置 947.11對流包 947.12離散包 947.13查閱結(jié)果 947.14源匯混合包 957.15保存模擬 957.16運(yùn)行模型 957.17運(yùn)行MT3DMS 957.18查閱結(jié)果 967.19物質(zhì)—時(shí)間圖 967.20動(dòng)畫顯示 977.21吸附和分解的模擬 977.22運(yùn)行選項(xiàng) 987.23保存模型 997.23.1運(yùn)行MT3DMS 997.24查看計(jì)算結(jié)果 997.25隨時(shí)間變化歷史圖 997.26小結(jié) 1001GIS模塊Shapefile文件導(dǎo)入、顯示以及轉(zhuǎn)換目標(biāo)無須使用ArcObjects實(shí)現(xiàn)Shapgefile文件的導(dǎo)入和顯示。在GMS中將shapefile文件轉(zhuǎn)化為屬性對象

1.1簡介使用GIS模塊可以在GMS中直接顯示GIS數(shù)據(jù)庫中的數(shù)據(jù)，而不必將數(shù)據(jù)轉(zhuǎn)換為GMS數(shù)據(jù)格式。GMS自身的數(shù)據(jù)如網(wǎng)格和鉆孔數(shù)據(jù)可以和GIS數(shù)據(jù)一起顯示。GIS模塊還可以選擇一部分GIS數(shù)據(jù)，并將其轉(zhuǎn)換為GMS的數(shù)據(jù)格式，用來建立地下水模型。目前GMS中的GIS模塊只能應(yīng)用于穩(wěn)態(tài)數(shù)據(jù)。若計(jì)算機(jī)注冊了ArcObjects，ESRI軟件中的許多模塊（如ArcMap）就可以在GMS中使用了。如果沒有注冊同樣可以使用GIS模塊，但其中的一些功能是不能使用的。這種情況下可以導(dǎo)入并顯示shapefile文件，并將其轉(zhuǎn)換為GMS屬性對象。GMS幫助文件中完整地列出了不同模型可以使用的模塊。本文中介紹了GMS軟件中的GIS模塊，其中包括兩部分。第一部分介紹了計(jì)算機(jī)在沒有注冊ArcObjects的情況下可以使用的模塊。第二部分介紹了注冊ArcObjects的情況下可以使用的模塊，如果計(jì)算機(jī)沒有注冊ArcObjects則不必看第二部分內(nèi)容。1.2.1綱要操作流程為：打開shapefile文件查看屬性表將shapefile文件轉(zhuǎn)換為散點(diǎn)啟動(dòng)ArcObjects并讀入圖層和shapefile文件轉(zhuǎn)換shapefile文件為屬性對象1.2未安裝ESRIArcObjects首先導(dǎo)入一個(gè)包含經(jīng)數(shù)據(jù)的點(diǎn)狀shapefile文件。然后建立shapefile點(diǎn)要素的離散點(diǎn)。在MAP模塊中建立GMSMODFLOW概念模型，同時(shí)建立包含井屬性的圖層。最后將shapefile文件轉(zhuǎn)換為GMS中的井，用來幫助建立模型。從上述過程中可以看出如何將GIS數(shù)據(jù)導(dǎo)入并轉(zhuǎn)換為GMS數(shù)據(jù)。這一部分需要的模塊有：MAP和GIS。1.3開始如有必要首先運(yùn)行GMS。若GMS已經(jīng)運(yùn)行，使用File/new這個(gè)命令，以便程序設(shè)置存儲(chǔ)為默認(rèn)的狀態(tài)。1.4讀取Shapefile文件第一步是讀取Shapefile文件。選擇Open按鈕。查找并打開文件目錄：tutfiles\GIS\gis在Open對話框中，將Filesoftype改為Shapefiles（*.shp）選擇文件名為arcmap.shp的文件并單擊Open按鈕。1.5查看Shapefile文件一定數(shù)量點(diǎn)顯示在屏幕上。在shapefile文件中含有井的信息。在ProjectExplorer中選擇GISLayersFolder。在ProjectExplorer中根據(jù)需要展開GIS數(shù)據(jù)層文件夾。可以看到arcmap.shpobject對象顯示在目錄樹中，這就是剛剛打開的shapefile文件。如果沒有安裝ESRIArcObjects，此時(shí)只有有限的幾個(gè)按鈕用來顯示shapefile。選擇DisplayOptions按鈕。單擊顯示點(diǎn)符號類型按鍵。改變點(diǎn)的半徑和顏色。單擊OK推出對話框。此時(shí)可以看到點(diǎn)的顏色發(fā)生了變化。1.6查看屬性表打開的shapefile文件中包含大量的關(guān)于點(diǎn)的屬性數(shù)據(jù)，可以通過屬性表進(jìn)行瀏覽。1、在目錄樹中，右鍵單擊arcmap.shp。2、在彈出的菜單中選擇AttributeTable命令。3、通過拖動(dòng)右下角改變對話框大小以瀏覽更多的數(shù)據(jù)。每一個(gè)點(diǎn)都包含大量的屬性（列）。注意到屬性表中的數(shù)據(jù)是不可以編輯的。4、點(diǎn)擊OK退出該對話框。1.7文件轉(zhuǎn)換為2D離散點(diǎn)現(xiàn)在要將GIS數(shù)據(jù)轉(zhuǎn)換為2D離散點(diǎn)，用來進(jìn)行插值。5、在目錄樹中，右鍵單擊arcmap.shp。在彈出的菜單中選擇Convertto2DScatterPoints命令。此時(shí)可以看到屏幕上出現(xiàn)一些新的符號。這些符號就是新建立的2D離散點(diǎn)。如有必要展開2DScatterDate文件夾。注意arcmap.shp2D散點(diǎn)數(shù)據(jù)是剛剛創(chuàng)建的。展開arcmap.shp對象。注意，GMS會(huì)根據(jù)屬性表中的每一個(gè)數(shù)字屬性自動(dòng)建立數(shù)據(jù)集。用戶可以使用這個(gè)程序建立離散點(diǎn)，然后將這些離散點(diǎn)插入到網(wǎng)格中。例如，如果每一個(gè)點(diǎn)的水位信息，就可以建立模型的初始水頭數(shù)居集合。因?yàn)椴粚@些離散點(diǎn)進(jìn)行編輯，需要將其刪除。在目錄樹中點(diǎn)擊2D散點(diǎn)數(shù)據(jù)文件夾并在彈出的菜單中選擇刪除。將GIS數(shù)據(jù)轉(zhuǎn)換為GMS特征屬性后，就可以用于建立概念模型。首先要選用一定的屬性建立起概念模型和圖層集。1.7.1創(chuàng)建概念模型1、在目錄樹中，右鍵點(diǎn)擊空白處，然后在彈出的菜單中選擇New|ConceptualModel命令。2、將名稱更改為Model1。3、確保Model設(shè)置為MODFLOW并點(diǎn)擊OK。4、在目錄樹中，右鍵單擊Model1文件夾，在彈出的菜單中選擇NewConceptualModel命令。5、將新的圖層重命名為coverage1。6、在Sources/Sinks/BCs選項(xiàng)組中選擇Wells并點(diǎn)擊OK。1.7.2根據(jù)GIS數(shù)據(jù)作圖1、在目錄樹中選擇GISLayers文件夾。2、選擇GIS|ShapesFeatureObjects命令。3、點(diǎn)擊Yes確保使用所有可見的圖形。此時(shí)彈出GIStoFeatureObjectsWizard對話框。4、點(diǎn)擊Next。5、在WELLNAME列中，將Mapping改為Name。6、在PUMPRATE列中，將Mapping改為Flowrate。7、點(diǎn)擊NEXT。8、點(diǎn)擊Finish。此時(shí)屬性對象點(diǎn)和GIS點(diǎn)位于同一位置，因此用戶不會(huì)注意到顯示中的區(qū)別。9、在目錄樹中取消arcmap.shp的勾選。10、在目錄樹中，右鍵單擊coverage1，在彈出的菜單中選擇AttributeTable

命令。這個(gè)對話框中顯示了圖層中所有要素點(diǎn)的屬性。如前所述名稱和流速都是由GIS屬性數(shù)據(jù)轉(zhuǎn)換而來的。但是所有點(diǎn)的類型均為NONE。需要將其設(shè)置為井。11、找到表格中對應(yīng)于All行和Type列的單元格，將其類型改為well。由于其處于All行，因此所有的點(diǎn)均改為井。12、點(diǎn)擊OK。1.7.3第一部分小結(jié)利用井的點(diǎn)文件可以進(jìn)一步建立MODFLOW概念模型。在MODFLOW-ConceptualModelApproach指南中將進(jìn)行更為詳盡的介紹。1.8已安裝ESRIArcObjects如果計(jì)算機(jī)或網(wǎng)路已經(jīng)注冊了ArcObjects，就可以進(jìn)行后面的操作。如果安裝了像ArcGIS一類的ESRI軟件，則表明已經(jīng)注冊了ArcObjects。如果無法分辨是否已經(jīng)注冊ArcObjects，則此導(dǎo)則將會(huì)顯示如何識別是否已經(jīng)注冊。如果已經(jīng)注冊ArcObjects的話，在GIS模塊中更多的功能屬性將可以使用。本文將介紹此類功能屬性。1.9啟動(dòng)ArcObjects刪除之前所有的操作并啟動(dòng)ArcObjects。選擇新建按鈕。在提示保存對話框中選擇NO。選擇GIS模塊。選擇GIS|EnableArcObjects菜單命令。如果在菜單命令之后有標(biāo)識顯示，則表明已經(jīng)計(jì)算機(jī)中注冊了ArcObjects，那么就可以按照指南進(jìn)行接下來的操作。1.10讀出數(shù)據(jù)第一步是讀出所需的數(shù)據(jù)。1、選擇GIS|AddData命令。2、查找并打開該目錄：tutfiles\GIS\gis。3、選擇文件名為streams.shp，Ndavis.tif，Sdavis.tif，recharge.shp，area_interest.shp的文件。4、選擇Add按鈕。1.11查看shapefile文件不同層的顯示次序是由目錄樹中各項(xiàng)的位置決定的。為了更好地查看GIS數(shù)據(jù)層，需要重新排列數(shù)據(jù)。1、如有必要展開目錄樹中的GISLayers文件夾。2、在目錄樹中，右鍵單擊area_interestshapefile，并選擇放大圖層命令。3、在目錄樹中，將各項(xiàng)按照下表所示進(jìn)行排列，并根據(jù)窗口大小進(jìn)行調(diào)整。圖1GIS顯示順序4、在目錄樹中，右鍵單擊streamsshapefile，并選擇屬性命令。選擇Symbology按鍵。點(diǎn)擊Symbol部分的按鈕。在Symbol選項(xiàng)中選擇River項(xiàng)，并點(diǎn)擊OK。選擇OK并退出屬性對話框。為了更好的顯示背景圖像，可以增加源匯圖層的透明度，關(guān)掉當(dāng)前地區(qū)的顏色填充。9、在目錄樹中，右鍵單擊rechargeshapefile，并選擇SetLayerTransparency命令。10、將透明度改為50%，并點(diǎn)擊OK。11、在目錄樹中，右鍵單擊streamsshapefile，并選擇屬性命令。12、選擇Symbology按鍵。13、點(diǎn)擊Symbol部分的按鈕。14、選擇Hollow類型，并將OutlineColor改為Red。15、雙擊OK退出對話框。1.12圖形文件轉(zhuǎn)化為屬性對象接下來將GIS數(shù)據(jù)轉(zhuǎn)化為屬性對象。首先要建立一個(gè)由一定屬性構(gòu)成的概念模型和圖層集。1.12.1創(chuàng)建概念模型1、在目錄樹中，右鍵單擊空白處，然后再彈出的菜單中，選擇New|ConceptualModel命令。2、更改名稱為Model1。3、確保模型設(shè)置為MODFLOW并點(diǎn)擊OK。4、在目錄樹中，右鍵單擊Model1概念模型并從彈出的對話框中選擇NewCoverage命令。5、重命名圖層為Rivers。6、在Sources/Sinks/BCs列表中選擇River選項(xiàng)，并點(diǎn)擊OK。7、在目錄樹中，右鍵單擊Model1概念模型，并從彈出的對話框中選擇NewCoverage命令。8、重命名圖層為Recharge。9、在ArealProperties列表中選擇Rechargerate選項(xiàng)，并單擊OK。1.12.2根據(jù)GIS數(shù)據(jù)作圖首先將streamsshapefile文件轉(zhuǎn)化為Rivers圖層，只需根據(jù)位置在區(qū)域范圍內(nèi)選定目標(biāo)即可成圖。1)在目錄樹中，選擇Rivers圖層，使其處于可激活狀態(tài)。2)在目錄樹中，選擇GISLayer文件夾。3)在目錄樹中，通過位置菜單選擇GIS|Selection|Select。4)確保菜單與下圖（圖2）中一致，并點(diǎn)擊應(yīng)用。圖2通過位置選擇河流屬性5)點(diǎn)擊Close。6)選擇選擇GIS|ShapesFeatureObjects命令。此時(shí)GIStoFeatureObjectsWizard顯示。7)點(diǎn)擊Next。8)在TYPE列中，將mapping改為Type。9)在NAME列中，將mapping改為Name。10)點(diǎn)擊Next。11)點(diǎn)擊Finish。12)在目錄樹中，選擇Recharge圖層，使其處于可激活狀態(tài)。13)在目錄樹中，選擇GISLayer文件夾。14)在目錄樹中，通過位置菜單選擇GIS|Selection|Select。15)確保菜單與下圖（圖3）中一致，并點(diǎn)擊應(yīng)用。圖3通過位置選擇源匯項(xiàng)屬性16)點(diǎn)擊Close。17)選擇GIS|ArcObjectsFeaturesObjects命令。18)點(diǎn)擊NEXT。19)在RECHRATE列中，將mapping改為Rechargerate。20)點(diǎn)擊NEXT。21、點(diǎn)擊Finish。在GIS區(qū)圖形所在的位置就生成了區(qū)屬性對象。22)在目錄樹中勾掉GISlayers文件夾選項(xiàng)。23)在目錄樹中，右鍵單擊Recharge并在彈出的對話框中選擇AttributeTable命令。24)將對象屬性類型改為Polygons。該對話框表明了在圖層中的所有的區(qū)屬性。應(yīng)注意的是源匯強(qiáng)度由GIS屬性表中數(shù)據(jù)轉(zhuǎn)換而來。25）點(diǎn)擊ok。截至到此，河流弧段和源匯項(xiàng)的區(qū)文件就可用來建概念模型。在MODFLOW-ConceptualModelApproach指南中將進(jìn)行更為詳盡的介紹，此處不再贅述。2屬性對象利用點(diǎn)、線段和多邊形文件建立獨(dú)立的概念模型2.1目標(biāo)本指南旨在說明如何利用屬性對象—點(diǎn)線段和多邊形—建成相互獨(dú)立的概念模型。2.1簡介本指南對屬性對象進(jìn)行了介紹。屬性對象在GMS中得到了大量使用，是最基本的元素。2.1.1綱要操作流程：1、創(chuàng)建圖層2、創(chuàng)建概念模型3、創(chuàng)建屬性對象4、選擇并修改屬性對象2.2開始如有必要首先運(yùn)行GMS。若GMS已經(jīng)運(yùn)行，使用File/new這個(gè)命令，以便程序設(shè)置存儲(chǔ)為默認(rèn)的狀態(tài)。2.3屬性對象屬性對象包括點(diǎn)、節(jié)點(diǎn)、線段和多邊形，見圖1，與地理信息系統(tǒng)（GIS）中相似。在GMS中屬性對象有多種用途。圖1屬性對象2.3.1Points點(diǎn)Points是xy坐標(biāo)的定位，并不在線段上，每個(gè)點(diǎn)對應(yīng)一個(gè)ID，并可以賦予屬性，通常用來代表井。2.3.2Arcs弧線Arcs是線段或邊組成的獨(dú)立折線體，有唯一的ID，也可以賦予屬性。2.3.3NodesandVertices端點(diǎn)和節(jié)點(diǎn)Nodes是弧線的兩個(gè)端點(diǎn)，有自己的ID并可以定義性質(zhì)。Vertices是兩個(gè)端點(diǎn)之間的點(diǎn)，僅用來定義弧線的幾何形狀，沒有屬性。2.3.4Polyons多邊形Polyons是一組弧線形成的閉合回線，可以由一個(gè)弧線或多個(gè)弧線祖成。如果兩個(gè)多邊形相連，兩個(gè)多邊形之間的邊界為共用的，不是兩條覆蓋的。2.3.5Coverages圖層屬性對象的集合形成了圖層，每個(gè)圖層代表一組數(shù)據(jù)。2.3.6CoveragesandConceptualModels圖層及概念模型創(chuàng)建圖層：在目錄樹中，右鍵單擊并選擇New|Coverage命令。點(diǎn)擊OK退出CoverageSetup對話框。Newcoverage選項(xiàng)是指當(dāng)前處于激活狀態(tài)的圖層，表明創(chuàng)建的新的對象屬性將添加于改圖層中。在目錄樹中，右鍵單擊newcoverage然后從菜單中選擇Duplicate命令（復(fù)制）。圖層是可以復(fù)制的。所有的屬性對象以及屬性表都可以由原始圖層復(fù)制到新圖層中。需要注意的是再newcoverage旁的選項(xiàng)鍵處于灰色狀態(tài)時(shí)表明newcoverage不再處于激活狀態(tài)。右鍵單擊“Copyofnewcoverage”圖層，并從菜單中選擇Delete命令。概念模型在GMS中概念模型使得一個(gè)或多個(gè)相關(guān)的圖層集成在一起。1、右鍵單擊目錄樹，然后從菜單中選擇New|ConceptualModel命令。圖2概念模型屬性對話框應(yīng)該注意的是概念模型具有名稱和類型，其類型與所要建的模型一致（MODFLOW，F(xiàn)EMWATER等）。概念模型的類型決定了在數(shù)據(jù)表中的其他可用的可用的選項(xiàng)內(nèi)容。2、點(diǎn)擊OK鍵退出ConceptualModelProperties對話框。3、在目錄樹中選擇并拖拽“NewModel”下方的“newcoverage”。4、在屬性表的警告對話框中點(diǎn)擊Yes。圖層的屬性表取決于與圖層相關(guān)的概念模型的設(shè)置。圖層應(yīng)位于概念模型下方。圖3目錄樹中位于概念模型下方的圖層2.3.7CreatingFeatureObjects創(chuàng)建屬性對象創(chuàng)建屬性對象的步驟：1、移動(dòng)鼠標(biāo)至CreatePoint工具。2、在地形圖上任意位置點(diǎn)擊鼠標(biāo)生成點(diǎn)。3、移動(dòng)鼠標(biāo)至CreateArc工具。4、點(diǎn)擊鼠標(biāo)生成線段。單擊開始繪制線段，雙擊結(jié)束繪制。當(dāng)閉合線段時(shí)就會(huì)生成多邊形圖形，但一旦使用BulidPolygons命令，將只生成線段而不是多邊形。建議隨意嘗試使用屬性對象工具。2.3.8SelectingFeatureObjects選擇屬性對象生成屬性對象可以使用多種工具。見下圖5。選擇不同類型的屬性對象應(yīng)選擇不同的工具。1、將鼠標(biāo)在鼠標(biāo)上停留幾秒鐘，就會(huì)出現(xiàn)工具名稱，如圖4所示。圖4鼠標(biāo)顯示工具名稱選擇工具可以用來選擇不同類型的屬性對象。其余工具僅僅可以選定特定類型的對象。選擇工具使用頻率較高，但必要時(shí)刻還需使用其他工具。如多個(gè)不同類型的屬性對象距離靠近時(shí)。2、利用選擇工具選定不同類型的屬性對象。3、嘗試使用不同的選擇工具，選定不同類型的對象。圖5屬性對象工具2.3.9ModifyingFeatureObjects修改屬性對象移動(dòng)目標(biāo)點(diǎn)及端點(diǎn)等屬性對象可以通過拖拽的方式將其移動(dòng)到其他位置。其他對象如弧段和多邊形是不能拖拽移動(dòng)的，但是可以拖拽線段或多邊形的節(jié)點(diǎn)和端點(diǎn)。點(diǎn)/端點(diǎn)使用SelectPoints/NodesTool工具選定之前生成的點(diǎn)。點(diǎn)住鼠標(biāo)左鍵，拖拽點(diǎn)到不同位置。在線段上利用相同SelectPoints/NodesTool工具拖拽節(jié)電。線段選擇DisplayOptions按鈕。確保左邊列表中的MapData選項(xiàng)處于高亮顯示，并且右邊的Map處于可視狀態(tài)。激活Vertices選項(xiàng)并點(diǎn)擊OK。移動(dòng)至SelectVertex工具。選擇并拖拽弧段上的節(jié)點(diǎn)。添加/移動(dòng)節(jié)點(diǎn)通過添加更多的節(jié)點(diǎn)可以修改線段，操作流程如下：移動(dòng)鼠標(biāo)至CreateVertex工具。點(diǎn)至之前生成線段并添加節(jié)點(diǎn)。通過添加和移動(dòng)節(jié)點(diǎn)位置可以改變線段的形狀。同樣可以通過右鍵單擊線段并重新生成節(jié)點(diǎn)，此為更為快速的方法。2.3.10GridFrames柵格框架柵格框架可以用來創(chuàng)建2D或者3D柵格。柵格框架可以定義柵格的矩形范圍，同樣可以用來旋轉(zhuǎn)或移動(dòng)到模型的任意適宜的位置。1、選擇FeatureObjects|NewGridFrame菜單命令。這樣就可以生成柵格框架，同樣可以在屏幕上看到一個(gè)紫色的矩形。2、移動(dòng)鼠標(biāo)至SelectGridFrame工具。點(diǎn)擊并選定柵格框架的其中一條線。點(diǎn)擊柵格框架內(nèi)部將其拖拽至其他位置。圖6柵格框架操作應(yīng)該注意的是位于角落里的操作圖標(biāo)和柵格框架的中部邊界。同時(shí)也應(yīng)注意圓形的操作圖標(biāo)是用來旋轉(zhuǎn)底部的右角。點(diǎn)擊柵格框架中的操作圖標(biāo)，改變其大小并旋轉(zhuǎn)柵格框架。柵格框架的屬性同樣也可以雙擊目錄樹中的柵格框架后手動(dòng)輸入，同樣也可對對話框?qū)傩赃M(jìn)行編輯。2.4結(jié)論通過本指南可以得出以下幾點(diǎn)：1、屬性對象主要包括點(diǎn)、線段節(jié)點(diǎn)和多變形。2、屬性對象易于生成和修改。3、屬性對象通過圖層進(jìn)行整合。3MODFLOW—概念模型法利用概念模型建立MODFLOW模型目標(biāo)概念模型應(yīng)用方式包括在map模塊中使用gis工具開發(fā)一個(gè)被模擬區(qū)域的概念模型。源匯項(xiàng)的位置，圖層參數(shù)如水力傳導(dǎo)系數(shù)，模型邊界和其他所有模擬需要的數(shù)據(jù)都無需網(wǎng)格，而只需在概念模型中進(jìn)行定義。必備知識指南層高數(shù)據(jù)差值指南MODAEM指南

所需模塊GridGeostatisticsMapMODFLOW

時(shí)間30-60分3.1簡介在GMS中，有兩種應(yīng)用方式可以構(gòu)建modflow模擬：柵格和概念模型。柵格方式包括直接在3D柵格的逐個(gè)單元（cell-by-cell）添加源/匯項(xiàng)和其他模型參數(shù)。該部分內(nèi)容在MODLOW-GridApproach指南中有所介紹。概念模型應(yīng)用方式包括在map模塊中使用gis工具開發(fā)一個(gè)被模擬區(qū)域的概念模型，源匯項(xiàng)的位置，圖層參數(shù)如水力傳導(dǎo)系數(shù)，模型邊界和其他所有模擬需要的數(shù)據(jù)都在概念模型中進(jìn)行定義。一旦概念模型建立，網(wǎng)格就自動(dòng)生成，概念模型自動(dòng)轉(zhuǎn)化為網(wǎng)格模型，并且所有需要獨(dú)個(gè)單元操作的設(shè)置都自動(dòng)進(jìn)行。本節(jié)對該部分內(nèi)容及利用概念模型進(jìn)行MODFLOW模擬的相關(guān)內(nèi)容進(jìn)行了詳細(xì)介紹。3.1.1綱要操作流程：導(dǎo)入背景底圖創(chuàng)建并建立coverages建立3Dgird模型導(dǎo)入散點(diǎn)虎踞并且差值為高度數(shù)據(jù)將概念模型轉(zhuǎn)為MODFLOW運(yùn)行Checkthesimulation和runMODFLOW查閱運(yùn)行結(jié)果3.2問題描述本文中例題所要解決的問題如圖1所示。該案例位于得克薩斯州東部，設(shè)定某一規(guī)劃垃圾填埋場其地下水可能遭到污染。在MODPATH和MT3DMS指南中，本次水流模擬的流場結(jié)果將用于粒子追蹤和運(yùn)移模擬。Fig1模擬區(qū)域(a)模擬區(qū)域平面圖(b)模擬區(qū)域典型南北剖面本次模擬賦存于山前沉積物的地下水，模擬范圍北至山前南至兩條河流的匯集處。圖1-b所示的為該模擬區(qū)域由北至南的一條典型剖面。該模擬區(qū)下覆灰?guī)r，在模擬邊界的北側(cè)山區(qū)地帶出露。本模擬區(qū)主要發(fā)育有兩個(gè)含水層，將上層含水層概化為潛水含水層，下層含水層改劃為承壓含水層。模型的北邊界為零流量邊界，其余邊界為與河流平均狀態(tài)相關(guān)的定水頭邊界。假定模型主要的補(bǔ)給來自于降水入滲。在模擬范圍內(nèi)發(fā)育有河床，時(shí)而干涸時(shí)而受到地下水的補(bǔ)給。在模型中將這些河床概化為排水溝。模擬范圍內(nèi)共有兩個(gè)生產(chǎn)井。注：本次模擬案例所在的位置是真實(shí)的，但是垃圾填埋場及該處的水文地質(zhì)條件是虛構(gòu)的。為了便于建立概念模型，本模擬案例選定的水文和邊界條件都是簡單清晰明了有代表性的樣本。3.3開始如有必要首先運(yùn)行GMS。若GMS已經(jīng)運(yùn)行，使用File/new這個(gè)命令，以便程序設(shè)置存儲(chǔ)為默認(rèn)的狀態(tài)。3.4導(dǎo)入模型底圖建立模型的第一步是導(dǎo)入模型區(qū)域的數(shù)字化圖形。該圖形掃描自USGS的地形圖，輸入該圖形到GMS后，注冊，并儲(chǔ)存GMS工程文件。打開工程文建就可以讀取圖形信息。一旦圖形輸入到GMS中，該圖形就可以作為背景顯示，并能夠?qū)⑵鋽?shù)字化并顯示模擬位置。3.4.1讀取圖形讀入圖形文件：選擇Open按鈕打開名為“tutfiles\MODFLOW\modfmap”的路徑。打開名為“srart.gpr”的工程文件。在GMS中其余的目標(biāo)圖層都位于圖形文件紙上。該圖形只能平面顯示。有時(shí)也可以在相同的文件夾中讀入其他的圖形文件。3.5存儲(chǔ)工程文件在對模型進(jìn)行任何修改前，首先另存該文件工程。選擇File/Save命令。存儲(chǔ)該工程文件為“easttex”?，F(xiàn)在可以不定期的點(diǎn)擊save按鈕以便保存模型的進(jìn)展。3.6定義單位在此處，定義該模型涉及到的單位。在GMS界面中可以對選用的單位進(jìn)行修改以便每個(gè)參數(shù)都有適用的單位。3、選擇Edit/Units命令。4、選擇“m”作為長度單位，選擇“d”作為實(shí)踐單位。其余單位將被忽略使用（該類單位在水流模型中不使用）。5、選用OK按鈕。3.7定義邊界第一步是沿著模型范圍通過弧段建成一個(gè)閉合的區(qū)域，從而定義模型的外邊界。3.7.1創(chuàng)建局部源／匯圖層1、在ProjectExplorer空白處單擊右鍵后彈出一個(gè)菜單，選擇New/ConceptualModel命令。2、在Name選型中，鍵入“EastTexas”。對于模型，選擇MODFLOW。3、點(diǎn)擊OK。4、右鍵單擊“EastTexas”概念模型然后在彈出的菜單中選擇NewCoverage命令。5、更改CoverageName為“Boundary”，將默認(rèn)高度值改為213。將默認(rèn)層數(shù)由1改為2。6、點(diǎn)擊OK。3.7.2創(chuàng)建弧1、選擇創(chuàng)建弧段工具2、通過點(diǎn)擊該位置模型左邊（西）開始創(chuàng)建弧段，如圖2所示。3、逆時(shí)針方向圍繞地圖邊界點(diǎn)擊邊界上一系列的點(diǎn)，從而創(chuàng)建弧段。不必?fù)?dān)心點(diǎn)的空間或確切地為之，盡量使用足夠多的點(diǎn)來大概的確定邊界位置。模型南部和西部邊界應(yīng)與河流保持一致。模型上部邊界應(yīng)與灰?guī)r出露位置一致，如圖2所示。4、結(jié)束弧段繪制，點(diǎn)擊最初開始的點(diǎn)。提示：在點(diǎn)擊點(diǎn)時(shí)，如果失誤了，可以倒退回去，只需點(diǎn)擊退格鍵。如果想要中斷當(dāng)前操作，按Esc鍵退出。圖2創(chuàng)建邊界弧段3.8創(chuàng)建本地源匯圖層建立概念模型的下一步是創(chuàng)建當(dāng)?shù)卦磪R圖層。這個(gè)圖層主要用于模型區(qū)域的邊界條件以及定義該地區(qū)的源匯項(xiàng)，包括開采井、河流、排水渠以及通用水頭邊界。通過在CoverageSetup對話框中的概念模型選項(xiàng)和選項(xiàng)設(shè)置，可以對模型屬性特征進(jìn)行確定。在創(chuàng)建屬性目標(biāo)前，首先要更改CoverageSetup對話框的選項(xiàng)。1、右鍵單擊Boundary圖層，在彈出的對話框中選擇Dupliacate命令，將新圖層命名為“Sources&Sinks”。2、右鍵單擊Sources&Sinks圖層，在彈出的對話框中選擇CoverageSetup命令。3、在Sources/Sinks/BCs列表中，將下列所需的選項(xiàng)設(shè)置為ON狀態(tài)：LayerrangeWellsRefinePointsSpecifiedHead(CHD)Drain4、確保Usetodefinemodelboundary(activearea)選項(xiàng)處于可用狀態(tài)。5、點(diǎn)擊OK。3.8.1定義定水頭邊界弧段下一步沿著模型的南部和西部邊界定義定水頭邊界。在此之前，首先要將創(chuàng)建的弧段分為三段。沿著模型上部的弧段定位為零流量邊界，其余兩條弧段定義為河流邊界。在弧段上選定一個(gè)或多個(gè)頂點(diǎn)，然后將定點(diǎn)轉(zhuǎn)為節(jié)點(diǎn)，就可以將弧段拆分為若干弧段。1、選擇SelectVerticestool工具。2、選定如圖3所示的兩個(gè)頂點(diǎn)。Vertex#1位于兩條河流的匯流處。Vertex#2位于模型西部的河流上方。一次同時(shí)選定這兩個(gè)頂點(diǎn)，然后選定其他頂點(diǎn)同時(shí)按住Shift鍵。3、右鍵單擊選定的頂點(diǎn)，并選擇命令。圖3頂點(diǎn)轉(zhuǎn)節(jié)點(diǎn)三條弧段定義生成之后，將選定的表征河流的弧段定位為定水頭邊界。4、選擇SelectVerticestool工具。5、選擇南部和西部邊界處一個(gè)弧后按住Shift鍵選擇其他弧段。6、右鍵單擊選定的弧段，然后在彈出的對話框中選擇AttributeTable命令。7、在單元數(shù)據(jù)表中找到相應(yīng)的All行和Type列。在該單元網(wǎng)格中，選擇spec.head類型。兩條弧段的類型都被修改。8、選擇OK按鈕。9、點(diǎn)擊模型除弧段以為的其他區(qū)域推出選定弧段。需要注意的是選定的弧段其顏色將會(huì)有所改變。下一步是定義弧段的端點(diǎn)處節(jié)點(diǎn)的水頭。定水頭邊界上的水頭是隨著弧段長度而線性變化的。10、選定SelectPoint/Nodes工具。11、雙擊南部邊界的東（左）節(jié)點(diǎn)。12、在Head-Stage中輸入常熟212。13、選擇OK按鈕。14、利用同樣的模式，將兩河流匯合處節(jié)點(diǎn)的值設(shè)置為208，將模型西部邊界的頂端節(jié)點(diǎn)的值設(shè)置為214。3.8.2定義排水溝弧段在此小節(jié)中，將河床通過弧段處理為排水溝。1、選擇CreateArctool2、創(chuàng)建如圖4所示的弧段1。從弧段的底部開始，通過沿河床點(diǎn)擊點(diǎn)來創(chuàng)建，在頂部弧段上雙擊結(jié)束。注意：在現(xiàn)有弧段的頂點(diǎn)或者是弧段的終點(diǎn)附件點(diǎn)擊是，GMS自動(dòng)的捕捉現(xiàn)有弧段并在兩條弧段相交的地方產(chǎn)生一個(gè)節(jié)點(diǎn)。3、用以上同樣的方法創(chuàng)建如圖4所示的2、3?；《巍D4排水溝弧段接下來，定義弧段為季節(jié)性水溝并指定弧段的傳導(dǎo)系數(shù)以及高程。4、選擇SelectArcs工具。5、通過按下shift鍵來選擇所有的排水溝弧段。6、在所選弧段上右鍵點(diǎn)擊并在下拉菜單中選擇AttributeTable菜單命令。7、在all行，type列中，選擇drain選項(xiàng)。8、在all行中輸入555作為傳導(dǎo)系數(shù)。這指示了傳導(dǎo)系數(shù)的單位長度值。當(dāng)季節(jié)性河流被分配給柵格單元時(shí)，GMS自動(dòng)計(jì)算恰當(dāng)?shù)膯卧獋鲗?dǎo)系數(shù)值。9、將每一個(gè)弧段的Fromlayer和Tolayer改變?yōu)?。意味著季節(jié)性河流僅出現(xiàn)在柵格的第一層中。10、選擇ok。在弧段的節(jié)點(diǎn)上指定季節(jié)性河流的高程值?；《紊系母叱讨凳窃诨《蔚膬蓚€(gè)指定值之間線性變化。11、點(diǎn)擊選擇SelectPoints/Nodes工具。雙擊圖5中的節(jié)點(diǎn)2。應(yīng)注意的是節(jié)點(diǎn)2由于該節(jié)點(diǎn)連接有2個(gè)不同的類型的弧段，因此該節(jié)點(diǎn)具有2個(gè)屬性。12、為drain特性的Bot.elev.輸入212。而不改變spec.head其他的任何特性。點(diǎn)擊ok。13、重復(fù)這些步驟指定圖5所示的季節(jié)性河流節(jié)點(diǎn)的高程值。14、確定僅改變drain特性，而不改變spec.head特性。圖5排水溝節(jié)點(diǎn)高程值3.8.3創(chuàng)建多邊形在區(qū)域源匯項(xiàng)圖層中，模擬的全部區(qū)域都必需用不重疊的多邊形覆蓋。定義柵格的活動(dòng)區(qū)域。多數(shù)情況下所有的多邊形都是可變水頭的多邊形（默認(rèn)的）。盡管如此，其他的多邊形也可以用。例如，模擬一條河流時(shí)可用常規(guī)的水頭多邊形。最簡單的方法是創(chuàng)建所有在圖層中要用的弧段并選擇BuildPolygons菜單命令來定義多邊形。此命令自動(dòng)搜索所有由弧段定義的每一個(gè)圈閉區(qū)域創(chuàng)建為多邊形。這些多邊形默認(rèn)的類型為“NONE”，可以通過選擇多邊形并選擇Properties菜單命令來修改其類型。圖層中的所有弧段都已經(jīng)被創(chuàng)造，接下來就要準(zhǔn)備構(gòu)建多邊形。所有的多邊形都是可變水頭的多邊形。1.選擇FeatureObjects|BuildPolygons菜單命令。注意：多邊形是被填充的?？梢酝ㄟ^選擇Display|DisplayOptions菜單命令來改變多邊形的視圖。3.8.4創(chuàng)建井定義區(qū)域源匯項(xiàng)的最后一步是定義井。井被定義為點(diǎn)狀目標(biāo)。下面創(chuàng)建兩個(gè)井。1.選擇CreatePoint工具。2.將鼠標(biāo)移動(dòng)到如圖1所示的井1位置，點(diǎn)擊一次創(chuàng)建一個(gè)點(diǎn)。3.新的點(diǎn)位創(chuàng)建之后，在GMS窗口上部的X，Y位置輸入坐標(biāo)（613250，3428630），按Tab或者是Enter鍵。4.選擇Properties按鈕。5.在type，選擇well選項(xiàng)。6.為Flowrate,輸入值-680。7.將Fromlayer和Tolayer都改變?yōu)?，意味著井只出現(xiàn)在柵格的第一層。8.選擇ok。9.相同的步驟創(chuàng)建另外一個(gè)井，輸入位置（615494，3428232），并設(shè)定抽水率-2830。對該井，將Fromlayer改為Tolayer，使得該井僅僅適用于第二層（對2進(jìn)行相同的修改）。GridRefinement柵格加工井是地下水流動(dòng)的聚集點(diǎn)，并且在井附近引起水頭陡坡（水躍，即水頭損失）。為了更為準(zhǔn)確的模擬井周圍的水流，井周圍的柵格將特別的精確（劃分更為細(xì)致？）。這在GMS中可以通過在概念模型中的井直接指派精確度來自動(dòng)完成。1.選擇SelectPoints/Nodes工具。2.通過按下shift鍵選擇兩個(gè)井。3.選擇Properties按鈕。4.找到Refine列，并在all行選擇。這將會(huì)打開兩個(gè)點(diǎn)的精度。5.將兩個(gè)點(diǎn)的Basesize改變?yōu)?5，Bias為1.1，Maxsize為150。6.點(diǎn)擊ok。3.9定義補(bǔ)給區(qū)域接下來的步驟是構(gòu)建概念模型中的定義補(bǔ)給區(qū)域。假設(shè)除垃圾填埋場以外的地區(qū)補(bǔ)給都是均一的。在垃圾填埋場區(qū)域的補(bǔ)給因?yàn)橐r墊而減小。3.9.1復(fù)制邊界通過復(fù)制邊界來創(chuàng)建補(bǔ)給圖層。1、在Boundary層中右鍵點(diǎn)擊并在下拉菜單中選擇Duplicate菜單命令。2、將新圖層名字修改為Recharge..。3、在Recharge.層中右鍵點(diǎn)擊并在下拉菜單中選擇CoverageSetup菜單命令。4、在ArealProperties列表中，打開Rechargerate選項(xiàng)。5、選擇ok。3.9.2構(gòu)建垃圾填埋場邊界下一步，將創(chuàng)建弧段來描繪垃圾填埋場的邊界。首先在垃圾填埋場大致區(qū)域創(chuàng)建一個(gè)方框。接下來邊界節(jié)點(diǎn)和頂點(diǎn)讓弧段于垃圾填埋場的邊界相一致。1、選擇CreateArc工具。2、創(chuàng)建一個(gè)如圖6所示的方框來描繪垃圾填埋場邊界。不必?fù)?dān)心點(diǎn)位的精度（下面會(huì)進(jìn)行精確定位）。方框弧段構(gòu)建后，下面進(jìn)行精確的定位節(jié)點(diǎn)和頂點(diǎn)。3、選擇SelectVertices工具。4、在垃圾填埋場的矩形方框周圍拉動(dòng)一個(gè)方框，選擇所有的頂點(diǎn)。5、在其中一個(gè)頂點(diǎn)上右鍵點(diǎn)擊，并選擇Vertex->Node菜單命令。6、選擇SelectPoints/Nodes工具。7、選擇方框拐角的一個(gè)節(jié)點(diǎn)。8、節(jié)點(diǎn)被選擇之后，在編輯窗口EditWindow.輸入節(jié)點(diǎn)的精確坐標(biāo)。輸入每一個(gè)坐標(biāo)值之后按Tab鍵。9、重復(fù)這個(gè)過程精確定位垃圾填埋場的其他幾個(gè)拐角節(jié)點(diǎn)。圖6垃圾填埋場3.9.3構(gòu)建多邊形弧段已經(jīng)創(chuàng)建，下面進(jìn)行構(gòu)建多邊形。1.選擇FeatureObjects|BuildPolygons菜單命令。3.9.4設(shè)置補(bǔ)給值補(bǔ)給區(qū)域已經(jīng)構(gòu)建，可以設(shè)置補(bǔ)給值。在模擬鐘，垃圾填埋場區(qū)域給定一個(gè)補(bǔ)給值，其余其余為另一個(gè)多變形，設(shè)定另一個(gè)補(bǔ)給值。1、選擇SelectPolygons工具。2、在填埋場多邊形單元上雙擊。3、將Rechargerate改為0.00006。注意：填埋場區(qū)域的補(bǔ)給率相對其他多邊形單元要小很多。垃圾填埋場因?yàn)樵O(shè)置了襯墊，進(jìn)行了壓實(shí)等原因?qū)е绿盥駡鲅a(bǔ)給值較小。此補(bǔ)給值本質(zhì)上是垃圾填埋場很小的滲濾液滲漏。4、選擇ok。5、在其他多邊形單元格上雙擊。6、將Rechargerate改變?yōu)?.00695。7、選擇ok。3.10定義水力傳導(dǎo)系數(shù)接下一步輸入每一層的水力傳導(dǎo)系數(shù)。多數(shù)情況下可能需要定義多個(gè)水力傳導(dǎo)系數(shù)分。為了簡便起見，每層將賦值一個(gè)常數(shù)作為水力傳導(dǎo)系數(shù)。3.10.1復(fù)制邊界通過復(fù)制邊界來創(chuàng)建圖層。1、在Boundary層中右鍵點(diǎn)擊并在下拉菜單中選擇Duplicate菜單命令。2、將新圖層名字修改為Layer1。3、在Layer1層中右鍵點(diǎn)擊并在下拉菜單中選擇CoverageSetup菜單命令。4、在ArealProperties列表中，打開以下選項(xiàng)：HorizontalKVerticalanis.5、將Defaultlayerrange改變?yōu)閺?到1（from1to1）。6、選擇ok。7、在Layer1層中右鍵點(diǎn)擊并在下拉菜單中選擇Duplicate菜單命令。將新圖層名字修改為Layer2。8、在Layer2層中右鍵點(diǎn)擊并在下拉菜單中選擇CoverageSetup菜單命令。9、將Defaultlayerrange改變?yōu)閺?到2（from2to2）。10、選擇ok。3.10.2頂層首先，給頂層設(shè)置一個(gè)K值。1、在ProjectExplorer.中選擇Layer1圖層。2、選擇FeatureObjects|BuildPolygons菜單命令。3、選擇SelectPolygons工具，并在多邊形上雙擊。4、將HorizonalK改變?yōu)?.5。5、將Verticalanis.改變?yōu)?。6、選擇ok。3.10.4下部圖層針對底層：1、在ProjectExplorer.中選擇Layer2圖層。2、選擇FeatureObjects|BuildPolygons菜單命令。3、在多邊形上雙擊。4、將HorizonalK改變?yōu)?0。5、將Verticalanis.改變?yōu)?。6、選擇ok。至此，完成了概念模型中的圖層定義。在繼續(xù)創(chuàng)建柵格前，將源匯項(xiàng)圖層為活動(dòng)（激活）圖層。7、在ProjectExplorer.中選擇Sources&Sinks圖層。3.11LocatingtheGridFrame定位柵格骨架至此圖層已經(jīng)完成，下面創(chuàng)建柵格。創(chuàng)建柵格的第一部是定義柵格中使用的柵格位置和方位。柵格骨架表現(xiàn)的是柵格的提綱?？梢栽谖覀兊奈粓D上定位。1、在ProjectExplorer中的空白區(qū)域右鍵點(diǎn)擊，在下拉菜單中選擇New|GridFrame菜單命令。2、在ProjectExplorer中右鍵點(diǎn)擊GridFrame，在下拉菜單中選擇FittoActiveCoverage菜單命令。3、在ProjectExplorer中柵格骨架中雙擊打開參數(shù)對話框。4、將Originz和Dimensionz分別修改為170和60。這是為MODFLOW圖層高程數(shù)列提供一個(gè)基礎(chǔ)數(shù)值點(diǎn)陣。后面，將會(huì)進(jìn)行插值修改（篡改）高程值。5、選擇ok退出GridFrame對話框。3.12創(chuàng)建柵格現(xiàn)在圖層以及柵格骨架都已經(jīng)創(chuàng)建。下面構(gòu)建柵格。1、選擇FeatureObjects|Map→3DGrid菜單命令。注意柵格是使用柵格骨架的數(shù)據(jù)構(gòu)建的。如果柵格骨架不存在，柵格默認(rèn)在模型周圍有大概5％的重復(fù)。并注意在X，Y方向上的單元格數(shù)量不能改變。這是因?yàn)樾辛械臄?shù)量和位置是由在井中輸入的點(diǎn)數(shù)據(jù)來控制的。2、在Z-Dimension中將Numcells改變?yōu)?。3、選擇ok。3.13初始化MODFLOW數(shù)據(jù)柵格已經(jīng)建立起來，并描繪出了活動(dòng)/不活動(dòng)區(qū)域，下一步是將概念模型轉(zhuǎn)變?yōu)闁鸥窕A(chǔ)的數(shù)字模式。在此之前，首先需要初始化MODFLOW數(shù)據(jù)。1、在ProjectExplorer.右鍵單擊grid選項(xiàng)并選擇NewMODFLOW命令。2、選擇ok按鈕。3.14定義活動(dòng)/不活動(dòng)區(qū)域柵格構(gòu)建之后，下一步是定義模型中活動(dòng)的和不活動(dòng)的區(qū)域。這在區(qū)域源匯項(xiàng)圖層中應(yīng)用信息自動(dòng)完成。1.在ProjectExplorer.中選擇MapData。2.選擇SelectPolygons工具。3.選擇其中的一個(gè)多邊形。4.選擇Properties按鈕。5.確認(rèn)圖層指派為1－2，選擇ok。6.選擇FeatureObjects|ActivateCellsinCoverage(s)菜單命令。區(qū)域源匯項(xiàng)圖層中，任何一個(gè)多邊形內(nèi)部的單元格都被指定為活動(dòng)的，多邊形外部的單元格都是不活動(dòng)的。注意邊界上的單元格都是活動(dòng)的，比如模型上部零通量邊界與邊界周圍的單元格是近似的，指定水頭邊界也與邊界上的單元格中心相一致。3.15插值（篡改）高程值接下來需要定義圖層高程以及初始水頭。使用LPF數(shù)據(jù)塊時(shí)，任意一個(gè)類型的圖層頂部和底部的高程值都已經(jīng)被定義。一個(gè)雙層的模型，需要為上部圖層定義一個(gè)圖層高程數(shù)列（地表），同時(shí)還有第一層的下部以及第二層的下部。前提是假設(shè)第二層的頂部等于第一層的底部高程。定義圖層高程的一個(gè)方法是輸入一個(gè)帶高程的散點(diǎn)，并直接插值圖層高程數(shù)列。有些時(shí)候是使用散點(diǎn)的集合。本次實(shí)例中，使用兩個(gè)散點(diǎn)集合：一個(gè)是地表，一個(gè)是第一層底部和第二層底部高程。由于是點(diǎn)源，因此利用兩個(gè)散點(diǎn)集是非常方便的。當(dāng)圖層高程來自于鉆孔數(shù)據(jù)時(shí)，地表點(diǎn)常常從地圖上數(shù)值化。在此案例中，地表高程數(shù)據(jù)來自國家高程數(shù)據(jù)庫（NED），此類數(shù)據(jù)可以利用網(wǎng)頁工具導(dǎo)入到GMS中。圖層插值法在InterpolatingLayerData中有詳細(xì)的介紹。3.15.1導(dǎo)入地表散點(diǎn)散點(diǎn)已經(jīng)包括在開始讀入的工程文件中。這些點(diǎn)來自于一個(gè)text文件（在2DGeostatistics指南中有描述）。因?yàn)樯Ⅻc(diǎn)已經(jīng)被隱藏，所以需要進(jìn)行顯示以便能觀察到。1、在ProjectExplore,中，打開2DScatterDatafolder2、檢查名為terrain和elevs.的兩個(gè)散點(diǎn)陣。3、通過在ProjectExplorer.中選擇terrain讓散點(diǎn)陣為激活狀態(tài)。一個(gè)散點(diǎn)的集合體符號將在模型中出現(xiàn)。3.15.2插值水頭和高程值下一步通過插值修改MODFLOW柵格的地表高程以及初始水頭。1、在terrain散點(diǎn)陣中右鍵點(diǎn)擊并選擇InterpolationTo|MODFLOWLayers菜單命令。這個(gè)對話框是告訴GMS你要修改哪一個(gè)MODFLW數(shù)列中的數(shù)據(jù)陣。此對話框在InterpolatingLayerData指南中有詳細(xì)的解釋。2、突出顯示ground_elev數(shù)據(jù)陣以及StartingHeads數(shù)列，并選擇Map按鈕。3、突出顯示ground_elev數(shù)據(jù)陣以及TopElevationsLayer1數(shù)列，并選擇Map按鈕。4、選擇ok按鈕執(zhí)行插值。3.15.3插值層高程1、選擇elevs散點(diǎn)陣激活它。2、在elevs散點(diǎn)陣上右鍵點(diǎn)擊并選擇InterpolationTo|MODFLOWLayers菜單命令。GMS自動(dòng)映射BottomElevationsLayer1和BottomElevationsLayer2數(shù)列基于其數(shù)據(jù)陣名字的恰當(dāng)數(shù)據(jù)陣列。3、選擇ok。3.15.4調(diào)整顯示現(xiàn)在已經(jīng)完成插值，可以隱藏散點(diǎn)陣。1、在ProjectExplorer.讓散點(diǎn)陣處于未選擇狀態(tài)。2、在ProjectExplorer.讓柵格骨架處于未選擇狀態(tài)。3.15.5顯示模型剖面通過顯示一個(gè)剖面來檢查插值的效果。1、在ProjectExplorer.中選擇3DGridData。2、在模型中部任意地方選擇一個(gè)單元格。3、選擇SideView按鈕。為了更好的顯示剖面圖，需要增加Z方向的比例。4、選擇DisplayOptions命令。5、在Zmagnification中輸入5。6、選擇ok按鈕。7、選擇Frame按鈕?？梢允褂肨oolPalette中的箭頭來顯示柵格中不同的列。注意，在剖面的右邊，底部圖層尖滅并且底部高程高于頂部高程。因此，在模型運(yùn)行之前必需進(jìn)行調(diào)整。3.15.6調(diào)整高程數(shù)列GMS提供了一個(gè)便利的修改圖層數(shù)列問題的工具。這些工具可以在ModelChecker中調(diào)用，并在InterpolatingLayerData指南中有詳細(xì)的解釋。1、選擇MODFLOW|CheckSimulation菜單命令。2、選擇RunCheck按鈕。3、在對話框的右側(cè)選擇FixLayerErrors按鈕。注意在圖層2中發(fā)現(xiàn)大量的錯(cuò)誤。這有幾種方法可以調(diào)整圖層。下面使用選項(xiàng)。這個(gè)選項(xiàng)可以讓底部圖層的所有單元格處于不活動(dòng)狀態(tài)。4、選擇Truncatetobedrock按鈕。5、選擇FixAffectedLayers按鈕。6、選擇ok按鈕退出FixLayerErrors對話框。7、選擇done按鈕退出ModelChecker.對話框。注意圖層錯(cuò)誤已經(jīng)被調(diào)整。另外一個(gè)顯示圖層修改效果的方法是在平面視圖。8、選擇PlanView按鈕。9、在小柵格現(xiàn)實(shí)中，選擇向下的箭頭顯示第二圖層。注意模型上邊（北部）邊界的單元格都是非活動(dòng)的。10、通過選擇向上箭頭轉(zhuǎn)換為上部圖層。3.16轉(zhuǎn)換概念模型下面從基于特征目標(biāo)的定義到基于柵格的MODFLOW數(shù)值模型來轉(zhuǎn)換概念模型。1、在EastTexas概念模型中右鍵點(diǎn)擊并在下拉菜單中選擇MapTo|MODFLOW/MODPATH菜單命令。2、確認(rèn)Allapplicablecoverages選項(xiàng)處于被選擇狀態(tài)并選擇ok。注意在季節(jié)性河流、井以及定水頭邊界下面的單元格都已經(jīng)被指定并且分配為恰當(dāng)?shù)脑磪R項(xiàng)。單元格的水頭以及高程沿著指定水頭和水溝弧段進(jìn)行線性的插值確定。單元格的傳導(dǎo)系數(shù)的確定是通過計(jì)算每個(gè)單元格與水溝弧段的交疊和指派給弧段的傳導(dǎo)系數(shù)值隨長度的倍增?？傊?，補(bǔ)給以及水力傳導(dǎo)系數(shù)值被賦值于恰當(dāng)?shù)膯卧瘛?.17檢查模擬至此，已完全定義了MODFLOW數(shù)據(jù)，并準(zhǔn)備運(yùn)行模擬。下面運(yùn)行ModelChecker來查看GMS能否找到一些模擬的錯(cuò)誤（在MODFLOW-GridApproach指南中已經(jīng)詳細(xì)的進(jìn)行過介紹ModelChecker）。1、在ProjectExplorer.中選擇3DGridData。2、選擇MODFLOW|CheckSimulation菜單命令。3、選擇RunCheck按鈕。應(yīng)該不會(huì)出現(xiàn)錯(cuò)誤。4、選擇done按鈕退出ModelChecker.。3.18運(yùn)行MODFLOW此時(shí)就可以運(yùn)行MODFLOW。1、選擇MODFLOW|RunMODFLOW菜單命令。此時(shí)MODFLOW被調(diào)用并出現(xiàn)了ModelWrapper。2、運(yùn)行完成后，選擇close按鈕。3.19顯示等水位線一組等水位線將顯示出來。為了更好的將等水位線和背景信息進(jìn)行對比，將等水位線的顏色修改為blue。1、在主工具欄中選擇ContourOptions。2、點(diǎn)擊ColorRamp按鈕，并在Color選項(xiàng)中選擇下拉鍵頭。3、選擇bule。4、再次選擇ok按鈕退出對話框。下面顯示第二層的等水位線。5、在小柵格顯示窗口中選擇向下箭頭。6、等水位線出現(xiàn)之后，通過選擇向上箭頭返回上層圖層。3.20側(cè)邊瀏覽水頭另外一個(gè)解決的方法是通過顯示側(cè)邊視圖來展示等水位線。1、選擇SelectCell工具。2、在模型右邊的井附近任意選擇一個(gè)單元格。3、選擇SideView按鈕。注意計(jì)算的水頭值被用于勾畫一個(gè)水頭輪廓。使用小柵格顯示區(qū)域的左右箭頭移動(dòng)到柵格的前后不同區(qū)域，可以看到一個(gè)井周圍的降落漏斗。4、選擇PlanView按鈕。3.21查閱水流均衡MODFLOW模擬是由一個(gè)水頭文件以及一個(gè)逐個(gè)單元流動(dòng)的流動(dòng)文件（CCF）組成的。GMS能使用CCF文件來顯示水流均衡狀況。例如，如果想知道從排水溝流出的水量情況，可以通過點(diǎn)擊排水溝弧段來完成。1、在ProjectExplorer.中選擇MapData。2、選擇SelectArcs工具。3、點(diǎn)擊最右邊的排水溝弧段。注意：通過弧段的總流量在窗口底部的任務(wù)欄條帶中顯示了出來。下一步，顯示河流的流量。4、點(diǎn)擊底部的任意一條指定水頭弧段顯示流量。5、按下shift鍵來選擇其他的指定水頭弧段。注意選擇的全部弧段的總流量已經(jīng)在下部顯示出來。可以通過下面的步驟顯示選擇單元格的流量。6、在ProjectExplorer.中選擇3DGridData。7、通過在井周圍拉動(dòng)一個(gè)方框選擇一些單元格。8、選擇MODFLOW/FlowBudget命令。此時(shí)下部對話框中顯示的所選單元格的總流量。9、選擇ok退出對話框。10、在模型外部任意地方點(diǎn)擊取消單元格的選擇狀態(tài)。3.22結(jié)論在此MODFLOW-ConceptualModelApproach指南中，下面是在此指南中需要學(xué)習(xí)的一些東西：可以通過輸入背景圖件來建立概念模型。在源匯圖層中定義模型邊界，并且在任何需要?jiǎng)?chuàng)建一個(gè)新的源匯圖層時(shí)通過拷貝此圖層來完成。可以在CoverageSetup對話框中定制與點(diǎn)、弧段以及多邊形相對應(yīng)的菜單（操作）。一些弧段的操作，比如水頭，不能通過選擇弧段來指定，但是可以通過選擇弧段兩端的節(jié)點(diǎn)來指定。原理是水頭值在弧段長度上沿弧段線性的改變來實(shí)現(xiàn)。柵格骨架可以用于定位柵格，但并不是必需的。6.每一次傳遞概念模型數(shù)據(jù)給柵格都必需通過MapMODFLOW/MODPATH菜單命令來完成。在概念模型中可以使用多邊形來指定類似圖層高程和水力傳導(dǎo)系數(shù)這類信息，但這可能導(dǎo)致階梯式的突變。為了平滑的過度，可以使用2D散點(diǎn)以及內(nèi)插值法。4MODFLOW-GridApproach利用3Dgrid方法建立MODFLOW模型目標(biāo)本文中介紹了利用柵格模式建立MODFLOW的方法。多數(shù)情況下，比起該類方法，利用概念模型方式建立模型更方便快捷。但由于柵格模式建模型可以進(jìn)行逐個(gè)單元格的編輯，因此此方法更適用于解決一些簡單問題或是學(xué)術(shù)問題。必備知識指南無

所需模塊GridMODFLOW

時(shí)間30-60分4.1簡介在GMS中，有兩種應(yīng)用方式可以構(gòu)建modflow模擬：柵格和概念模型。柵格方式包括直接在3D柵格的逐個(gè)單元（cell-by-cell）添加源/匯項(xiàng)和其他模型參數(shù)。概念模型應(yīng)用方式包括在map模塊中使用gis工具開發(fā)一個(gè)被模擬區(qū)域的概念模型。概念模型中的數(shù)據(jù)將會(huì)被復(fù)制到柵格中。本文中柵格應(yīng)用于modflow的預(yù)處理過程。大部分而言，概念模型的應(yīng)用比柵格方式更為有效。盡管如此，柵格方式在解決簡單問題或者需要逐個(gè)單元編輯的學(xué)術(shù)練習(xí)中也是非常有用的。在開始MODFLOW-ConceptualModelApproachtutorial.時(shí)不需要完成此tutorial。4.1.1綱要操作流程：1、導(dǎo)入背景底圖創(chuàng)建并建立coverages3、建立3Dgird模型4、導(dǎo)入散點(diǎn)虎踞并且差值為高度數(shù)據(jù)4.2問題描述此指南中需要解決的問題在圖1中描述。此問題是MODFLOWReferenceManual.結(jié)尾描述的簡單問題的改進(jìn)。在計(jì)算網(wǎng)格中使用三層假設(shè)成為三個(gè)含水層。柵格覆蓋區(qū)域?yàn)橐?5000ft×75000ft的正方形區(qū)域。柵格由15行和15列組成，計(jì)劃每個(gè)單元的計(jì)量為5000ft×5000ft。頂部和底部高程將簡化成為水平面。圖上顯示為水平方向的水力傳導(dǎo)系數(shù)，垂直方向的水力傳導(dǎo)系數(shù)使用部分水平方向的水力傳導(dǎo)系數(shù)。由于大氣降水入滲補(bǔ)給系統(tǒng)的輸入將會(huì)被定義為回灌。由于排水管道、排水井（圖中沒有標(biāo)示）排出系統(tǒng)，左邊的湖被定義為定水頭邊界。初始水頭定為0，穩(wěn)定態(tài)將會(huì)進(jìn)行計(jì)算。圖1案例概況4.3開始如有必要首先運(yùn)行GMS。若GMS已經(jīng)運(yùn)行，使用File/new這個(gè)命令，以便程序設(shè)置存儲(chǔ)為默認(rèn)的狀態(tài)。4.4單位在此處，定義該模型涉及到的單位。在GMS界面中可以對選用的單位進(jìn)行修改以便每個(gè)參數(shù)都有適用的單位。1、選擇Edit/Units命令。2、選擇“m”作為長度單位，選擇“d”作為實(shí)踐單位。其余單位將被忽略使用（該類單位在水流模型中不使用）。3、選用OK按鈕。4.5創(chuàng)建柵格解決此問題的第一步是創(chuàng)建一個(gè)3D有限差分單元格。1、在ProjectExplorer中，右鍵點(diǎn)擊空白出，在下拉菜單中選擇New/3DGrid命令。2、在X坐標(biāo)（X-dimension）區(qū)域，輸入22860作為長度值，單元格數(shù)Numbercells為15。3、在Y坐標(biāo)（Y-dimension）區(qū)域，輸入22860作為長度值，單元格數(shù)Numbercells為15。4、在Z坐標(biāo)（Z-dimension）區(qū)域，輸入單元格數(shù)Numbercells為3。稍后再輸入每一層?xùn)鸥竦捻敳亢偷撞扛叱讨?。因?yàn)椋斎氲腪方向的單元格厚度對于MODFLOW的計(jì)算不起作用。5、選擇ok按鈕。在顯示屏上會(huì)出現(xiàn)平面顯示的柵格。簡化的表示也會(huì)出現(xiàn)在小圖標(biāo)Mini-GridPlot.中。圖2柵格工具欄4.6創(chuàng)建modflow模擬下一步是初始化modflow模擬。1、在ProjectExplorer中，右鍵點(diǎn)擊Grid按鈕，在下拉菜單中選擇NewMODFLOW命令。4.6.1全局變量賦值MODFLOW的輸入被分為獨(dú)立的模塊化軟件包，有些是可以選擇的，有些是必需的，全局變量軟件包是必選的。軟件包首先進(jìn)行軟件包選擇。1、選擇Packages按鈕。這個(gè)塊對話框用于指定建立模型所需要的軟件包。基礎(chǔ)的軟件包是必需的，不能被關(guān)閉。下面選擇其他的軟件包。2、在Optionpackages部分，打開Drain(DRN1)和Well(WEL1)選項(xiàng)，打開Recharge(RCH1)選項(xiàng)。3、在Solver部分，選擇StonglyImpl.Procedure(SIP1)軟件包。4、選擇ok推出對話框。邊界條件下一步是進(jìn)行邊界條件賦值。邊界條件賦值用于指定單元格的活動(dòng)性，IBOUND>0為活動(dòng)，IBOUND＝0為不活動(dòng)，IBOUND<0為常數(shù)。在此案例中，除了上面兩層最左邊一列外所有的單位都為活動(dòng)，這一列將被定義為定水頭。1、選擇IBOUND按鈕。IBOUND對話框在一個(gè)電子表格中顯示IBOUND數(shù)列的值，每次顯示一層。左上角的編輯區(qū)域可以改變當(dāng)前層，除上面兩側(cè)左側(cè)一列的數(shù)據(jù)小于0外，其他所有的數(shù)列值都大于0。默認(rèn)狀態(tài)下，數(shù)列中的值都大于0。因此，需要改變定水頭單元格。這可以在30個(gè)定水頭單元格中輸入-1即可。盡管如此，還有另外更簡單的方法修改邊界條件。這將會(huì)在后續(xù)的指南中進(jìn)行描述?，F(xiàn)在將所有的單元格處于激活狀態(tài)。2、選擇ok退出對話框。初始水位下一步建立初始水位。1、選擇StartingHeads按鈕。初始水頭用于建立短暫模擬的初始水頭值。當(dāng)計(jì)算出一個(gè)穩(wěn)定態(tài)的模擬，在最后的解中每一個(gè)單元的初始值不應(yīng)該有差別。初始水頭值與最終水頭值越接近，modflow越快計(jì)算出結(jié)果值。對于某些類型的含水層來說，如果初始水頭太低，modflow將認(rèn)為此單元疏干。在此案例中，初始值為0既能滿足要求。初始水頭數(shù)列也用于建立定水頭單元的水頭值。在此問題中，定水頭值為0。所有的初始水頭值默認(rèn)為0，不需要修改。2、選擇ok退出對話框。含水層頂?shù)装甯叱滔乱徊绞墙㈨敳亢偷撞康母叱讨禂?shù)列。1.、選擇TopElevation按鈕。2、確定當(dāng)前層為1（Layeris1.）3、選擇Constant→Layer按鈕。4、輸入值為60并確定。5.、選擇ok退出TopElevations對話框。GMS默認(rèn)該層的上部也是上一層的底部。因?yàn)椋恍枰斎胨袑拥牡撞扛叱?，頂部高程將?huì)自動(dòng)賦值。6、選擇BottomElevation按鈕。7、確認(rèn)當(dāng)前層為1。8、選擇Constant→Layer按鈕。9、輸入-45并確認(rèn)。10、改變當(dāng)前層為2。11、選擇Constant→Layer按鈕。12、輸入-120并確定。13、改變當(dāng)前層為3。14、選擇Constant→Layer按鈕。15、輸入-215并確認(rèn)。16、選擇ok退出BottomElevation對話框。17、選擇ok退出MODFLOWGlobalPackage對話框。4.7單元邊界條件直接賦值正如上面提到的一樣，IBOUND值可以通過IBOUND數(shù)列對話框進(jìn)行輸入。對某些情況來說，這比直接對每個(gè)單元進(jìn)行輸入更為容易，可以通過CellProperties命令來完成。在使用此命令前，必須選擇上面兩層最左邊的一列單元格。4.7.1換角度觀察柵格為了便于操作，可以換轉(zhuǎn)視角從而實(shí)現(xiàn)不同方向查看構(gòu)建的柵格網(wǎng)絡(luò)。18、選擇SideView按鈕柵格非常薄。這可以通過調(diào)整Z方向上的放大率來進(jìn)行調(diào)整。19、選擇Display選項(xiàng)20、將Zmagnification修改為15。21、選擇ok。4.7.2選定單元格1、選擇SelectCells工具。2、在Mini-GridDisplay中將列改變?yōu)?，并按table鍵。注意現(xiàn)在觀察的是第一列（最左邊一列）3、在上面兩層的柵格周圍拉一個(gè)方框。4.7.3修改邊界條件編輯IBOUND值：1、在所選擇的其中一個(gè)單元上右鍵點(diǎn)擊。2、選擇Properties工具。3、將邊界條件選項(xiàng)改變?yōu)镾pecifiedhead.。4、選擇ok退出CellProperties對話框。5、選擇PlanView按鈕。注意：符號已經(jīng)顯示在所編輯的單元上，表明這些單元是定水頭單元。4.7.4檢查賦值確保輸入的邊界條件的準(zhǔn)確度：1、選擇MODFLOW|GlobalOptions菜單命令。2、選擇IBOUND按鈕。3、在左上角的對話框中選擇右側(cè)的向上箭頭，讓圖層循環(huán)顯示。注意上面兩層最左邊一列單元的值均為-1。在GMS中使用電子表格方式或者是通過選擇單元并直接輸入（選擇一個(gè)最為便利的輸入方式）的MODFLOW數(shù)據(jù)均能夠編輯。4、選擇ok按鈕退出IBOUNDArray對話框。5、選擇ok按鈕退出MODFLOWGlobalPackage對話框。4.8LPF模塊下一步是建立模型并輸入圖層流動(dòng)模塊（LPF）。LPF模塊計(jì)算每個(gè)柵格單元之間的傳導(dǎo)系數(shù)以及建立逐個(gè)單元流動(dòng)的差分方程。1、選擇MODFLOW|LPFPackage菜單命令。4.8.1圖層類型在對話框中LayerData部分的選項(xiàng)，被用于定于圖層類型以及每一個(gè)圖層的水力傳導(dǎo)系數(shù)。在我們的問題中，有三個(gè)圖層。最上面的一層是無限制的，下面兩層是有限制的（承壓的）。GMS中默認(rèn)的圖層類型是“可逆的（可轉(zhuǎn)變的）”，也就是說圖層可以能夠被轉(zhuǎn)變?yōu)槌袎夯驘o壓。因此，不需要改變圖層的類型。4.8.2圖層參數(shù)在對話框中LayerData部分的按鈕是為了輸入計(jì)算逐個(gè)單元的傳導(dǎo)率必需的參數(shù)。MODFLOW要求一組與圖層類型相關(guān)聯(lián)的圖層參數(shù)。4.8.3頂層首先，輸入頂層的數(shù)據(jù)。1、選擇HorizontalHydraulicConductivity按鈕。2、選擇Constant→Layer按鈕。3、輸入值15。4、選擇ok按鈕。5、選擇ok按鈕退出HorizontalHydraulicConductivity對話框。6、重復(fù)這個(gè)步驟并在垂向各向異性上輸入10。4.8.4中間層1、在LayerData部分選擇向上箭頭轉(zhuǎn)到正確的圖層編輯區(qū)域，轉(zhuǎn)到圖層2。2、在圖層2中輸入一下數(shù)值：、4.8.5底層1、轉(zhuǎn)到圖層3并輸入以下數(shù)值：至此數(shù)據(jù)輸入完成。2、選擇ok退出MODFLOWLPFPackage對話框。4.9補(bǔ)給（回灌）模塊接下來為補(bǔ)給模塊輸入數(shù)值。補(bǔ)給模塊用于模擬降雨入滲以及地表入滲補(bǔ)給含水層。輸入補(bǔ)給數(shù)據(jù)：1、選擇MODFLOW|Source/SinkPackages|RechargePackage菜單命令。2、選擇Constant→Array按鈕。3、輸入0.0009并點(diǎn)擊ok。4、選擇ok按鈕退出RechargePackage對話框。4.10排水溝模塊現(xiàn)在定義模型上層中的排水溝行。定義排水溝，首先要選擇排水溝所在的單元格，再選擇PointSources/Sinks菜單命令。4.10.1選擇單元格排水溝定位于上層（圖層1），這是當(dāng)前層，我們不需要改變顯示。需要選擇第8行的2-10列單元格。選擇單元格：1、選擇SelectCells工具。2、注意當(dāng)你在柵格上移動(dòng)鼠標(biāo)，鼠標(biāo)下面單元格所屬的ijk值（行，列，層）將在屏幕最下面的對話框中顯示出來，如圖3。圖3箭頭下方的IJK3、選擇i=8,j=2,k=1.的單元格4、按下shift鍵啟用多選模式并選擇同一行的3－10列單元格，如圖4.圖4選定的單元格4.10.2為排水溝賦值1、在所選單元格的圖形窗口中右鍵點(diǎn)擊，在下拉菜單中選擇Sources/Sinks菜單命令。2、選擇Drain表。3、選擇New按鈕。這將會(huì)為所選擇單元格增加一個(gè)新的排水溝實(shí)例。這次需要為選擇的排水溝輸入一個(gè)高程值和傳導(dǎo)系數(shù)。排水溝具有共同的傳導(dǎo)系數(shù)但是高程并不是全部相同。4、為排水溝單元格輸入下表所示的傳導(dǎo)系數(shù)以及高程值。5、選擇ok按鈕。6、在柵格外面任意地方點(diǎn)擊取消單元格選擇狀態(tài)。4.11井模塊下面使用PointSources/Sinks命令，通過選擇井所定位的單元格來定義幾個(gè)井。4.11.1上層井大部分井在最上面的圖層，但是也有一些在中間和底部圖層。首先定義上層的井。1、按下shift鍵選擇圖5所示的單元格。圖5頂層選定的單元格2、在所選的單元格上面點(diǎn)擊右鍵并選擇Sources/Sinks菜單命令。3、選擇well表格。4、選擇new按鈕。5、為所有的井輸入一個(gè)flow值為-12,230（負(fù)值表示抽水）。6、選擇ok按鈕。7、在柵格外面任意地方點(diǎn)擊取消單元格選擇狀態(tài)。4.11.2中間層的井現(xiàn)在定義中層的井，首先需要轉(zhuǎn)向顯示中間圖層。1、在Mini-GridPlot.中選擇Decrement按鈕。2、按下shift鍵選擇如圖6所示的井所在單元格。圖6中間層選定的單元格3、在所選的單元格上面點(diǎn)擊右鍵并選擇Sources/Sinks菜單命令。4、選擇well表格。5、選擇new按鈕。6、為所有的井輸入一個(gè)flow值為-12,230（負(fù)值表示抽水）。7、選擇ok按鈕。8、在柵格外面任意地方點(diǎn)擊取消單元格選擇狀態(tài)。4.11.3底層井在底層定義井時(shí)，步驟如下：1、在Mini-GridPlot.中選擇Decrement按鈕。2、選擇如圖7所示的井所在單元格。圖7底層選定的單元格3、在所選的單元格上面點(diǎn)擊右鍵并選擇Sources/Sinks菜單命令。4、選擇well表格。5、選擇new按鈕。6、為所有的井輸入一個(gè)flow值為-0.15（負(fù)值表示抽水）。7、選擇ok按鈕。8、在柵格外面任意地方點(diǎn)擊取消單元格選擇狀態(tài)。至此，所有的井均已定義，返回上部圖層。9、在Mini-GridPlot.選擇兩次。4.12檢查模擬上面已經(jīng)完全定義了MODFLOW數(shù)據(jù)，準(zhǔn)備運(yùn)行模擬。盡管如此，在保存模擬和運(yùn)行MODFLOW之前，需要運(yùn)行MODFLOW的ModelChecker并檢查錯(cuò)誤。因?yàn)槿绻呛Ａ繑?shù)據(jù)進(jìn)行模擬，常常會(huì)出現(xiàn)遺漏一些必需的數(shù)據(jù)或者是定義了不一致甚至是相反的選擇和參數(shù)。一些錯(cuò)誤會(huì)導(dǎo)致MODFLOW崩潰或者是產(chǎn)生錯(cuò)誤的解。MODFLOWModelChecker的目的是分析定義當(dāng)前MODFLOW模擬所輸入的數(shù)據(jù)并給出明顯的錯(cuò)誤或者可能的問題。成功運(yùn)行ModelChecker不能完成保證解是完成正確的。它僅僅是對輸入數(shù)據(jù)的初始檢查并節(jié)省大量的時(shí)間，另外可以捕捉到數(shù)據(jù)的輸入誤差。運(yùn)行ModelChecker：1、選擇MODFLOW|CheckSimulation菜單命令。2、選擇RunCheck按鈕。每一個(gè)MODFLOW的輸入模塊均在列表中給出。如果每一項(xiàng)均是恰當(dāng)?shù)?，那么每一個(gè)輸入模塊均不會(huì)顯示有錯(cuò)誤。如果出現(xiàn)錯(cuò)誤，選擇并突出顯示，相對應(yīng)的，GMS會(huì)選擇相應(yīng)有問題的單元格或者圖層。3、選擇Done按鈕退出ModelChecker.對話框。4.13保存模擬現(xiàn)在保存模擬并運(yùn)行MODFLOW。1、選擇File|Save