1、VisualMODFLOW4.2實例操作手冊ForW1ndow2000PProfessionalVisualMODFLOWPremiumversion丹聞涉ol呻k 引言MODFLOW是英文名稱ModularThree-dimensionalFiniteDifferenceGroundwaterflowmodel(三維有限差分地下水流模型)的簡稱。由美國地質(zhì)調(diào)查局(UniteStateGeologicalSurvey)于80年代開發(fā)出的一套專門用于孔隙介質(zhì)中地下水流動相關問題數(shù)值模擬的軟件。自問世以來，MODFLOW已經(jīng)在全世界范圍內(nèi)，在科研、生產(chǎn)、環(huán)境保護、城鄉(xiāng)發(fā)展規(guī)劃、水資源利用等許多行業(yè)

2、和部門得到了廣泛的應用，是供水文地質(zhì)工作者使用的一套功能強大的實用計算機軟件。同時，MODFLOW軟件的源代碼可以在美國內(nèi)政部網(wǎng)站上免費下載，供水文地質(zhì)工作者參考和修改，以提高數(shù)值模擬的仿真性。VisualMODLFOW軟件中不僅包括MODLFOW模塊，還有SEAWATMT3D、MODPATH等模塊。為了讓廣大水文地質(zhì)工作者能快速入門并精通該數(shù)值模擬軟件，作者參考相關文獻分別針對上述的模塊自擬幾個實例，并給出相應的操作步驟，供大家練習參考。目錄TOC o 1-5 h z HYPERLINK l bookmark10 實例1地下水流-污染物遷移規(guī)律數(shù)值模擬實例4 HYPERLINK l book

3、mark12 1、實例模型背景介紹及概念模型4 HYPERLINK l bookmark14 2、VM軟件模擬該問題所需要的資料以及資料整理格式5 HYPERLINK l bookmark16 數(shù)值模擬資料5 HYPERLINK l bookmark18 資料整理格式5 HYPERLINK l bookmark26 3、VM軟件模擬地下水流場:FLOWFIELD6 HYPERLINK l bookmark28 新建工程項目6 HYPERLINK l bookmark40 3.2模塊I：模型輸入8 HYPERLINK l bookmark184 模塊II：計算(Run)23 HYPERLINK

4、l bookmark202 3.4模塊I：輸出的可視化24 HYPERLINK l bookmark220 3.5模型校正26 HYPERLINK l bookmark226 4、VM軟件模擬地下水流跡線：MODPATH.27 HYPERLINK l bookmark228 4.1跡線模擬27跡線選項29 HYPERLINK l bookmark250 5、VM軟件模擬均衡區(qū)水量變化：ZONEBUDGET29 HYPERLINK l bookmark252 5.1均衡區(qū)水量模擬29 HYPERLINK l bookmark256 均衡區(qū)水量選項29 HYPERLINK l bookmark26

5、0 6、VM軟件模擬污染物遷移規(guī)律：MT3DS307、VM軟件反求參數(shù)：PEST32實例2某濱海地區(qū)海水入侵數(shù)值模擬實例341、實例模型背景介紹及概念模型34實例1地下水流污染物遷移規(guī)律數(shù)值模擬實例本實例是MODFLOW自帶的實例之一，主要讓水文地質(zhì)工作者練習地下水與污染物遷移規(guī)律的數(shù)值模擬的全過程，主要包括流場模擬（MODFLOW）、跡線數(shù)值模擬（MODPATH）、水均衡計算（ZoneBughet）與污染遷移規(guī)律模擬（MT3D）。1、實例模型背景介紹及概念模型在Waterloo城外有一個機場，機場附近加油站的有不斷的向含水層滲漏，為了預測與評價污染物的遷移規(guī)律，這里采用MODFLOW軟件中的

6、水流模型與MT3D模塊建立地下水流場與污染物遷移規(guī)律的數(shù)值模擬模型。經(jīng)過實地勘察，水平向研究區(qū)的范圍如圖1,在垂向上，將研究區(qū)分為三層：上下各有一層砂礫含水層，中間有一粘土和粉砂弱透水層將它們分開，其中，第一層是潛水含水層，第二層是弱透水層，第三層是承壓含水層，如圖2。含水層東西邊界與基巖相接，是隔水層。含水層北部邊界是定水頭邊界，第一層和第二層水頭值為19m，第三層水頭值為18m。南部含水層是定水頭邊界，水頭值是16.5m。相關的場地地物有一個飛機加油場、一個生活供水井場和一塊弱透水層不連續(xù)區(qū)（含水層天窗）。如圖1所示。圖1研究區(qū)域平面示意圖生活供水井場有兩眼井。東井抽水量穩(wěn)定在550m3/

7、d，西井為400m3/d。在過去十年中，在加油場中不斷地有飛機燃料的泄漏，其自然下滲已在上層含水層中產(chǎn)生了一個污染羽。本教程首先將向你顯示如何一步步地為這個場地建立一個地下水流動模型，然后介紹污染物遷移模擬。在討論時，我們將該場地在平面上的視圖指定為上北下南左西右東。如圖2所示，地下水在一個由上部潛水含水層、中間弱透水層和下部承壓含水層所組成的三層系統(tǒng)中自北向南（平面圖上是從上到下）流動。上下含水層的水力傳導系數(shù)為2e-4m/s,弱透水層的水力傳導系數(shù)為le-10m/s。降水資料：月份123456789101112降水補給10203040100100100100504030102、VM軟件模擬

8、該問題所需要的資料以及資料整理格式2.1數(shù)值模擬資料口空間展布相關資料地表高程、分層數(shù)據(jù)邊界位置口含水層參數(shù)：滲透系數(shù),彈性釋水系數(shù)（重力給水度）,孔隙度等口源匯項：大氣降水入滲系數(shù)（分區(qū)、數(shù)值）蒸發(fā)排泄系數(shù)地表水體（河流、湖泊、水庫等）水位、底面高程、底面巖性特征（厚度、滲透系數(shù)等）渠系灌溉入滲系數(shù)人工開采（點狀、面狀）：開采井位置、井結(jié)構、開采量動態(tài)口邊界條件：邊界類型、水頭或流量口初始條件：統(tǒng)測水位口水位動態(tài)觀測資料：觀測孔位置、結(jié)構、水位時間變化2.2資料整理格式采用VisualMODFLOW軟件模擬地下水流以及污染物遷移規(guī)律問題時，首先要通過分析現(xiàn)有的勘察資料，建立概念模型，其次將對

9、應的資料按照MODFLOW要求的格式整理，最后，按照先后順序?qū)①Y料輸入軟件，運行即可。下面分別介紹數(shù)值模擬資料的輸入格式。2.1.1空間展布資料的輸入格式2.1.1.1平面地圖為了將實際的研究輸入計算機，導入的底圖需要與實際數(shù)據(jù)匹配。如果是矢量圖AutoCAD、ARCGIS或者MAPGIS圖）可以直接導入，若使BMP格式的圖，需要給定圖中任意兩個點的坐標，在導底圖后將坐標如何計算機將底圖校正。2112分層數(shù)據(jù)MODFLOW軟件可以直接對建立好的網(wǎng)格賦分層數(shù)據(jù)，對于理想模型這種方式比較方便。由于實際問題的每一層的厚度隨著位置的變化而變化，常采用“數(shù)據(jù)導入”與“空間插值”的方式給每一個網(wǎng)格賦值，其

10、中數(shù)據(jù)的排列格式為XYElevation*（*代表實際數(shù)據(jù)）2.1.2含水層參數(shù)含水層參數(shù)包括滲透系數(shù)、彈性釋水系數(shù)（重力給水度）孔隙度等。在MODFLOW中有兩種方式可以給含水層賦值：第一種是直接賦值，若有參數(shù)分區(qū)時，首先導入?yún)?shù)分區(qū)圖，然后再賦值;第二種方式可以導入數(shù)據(jù)，導入數(shù)據(jù)格式與分層數(shù)據(jù)格式類似。2.1.3源匯項在地下水領域中，源匯項有很多，包括井、河流、溝渠、降雨、蒸發(fā)等，數(shù)據(jù)格式同分層數(shù)據(jù)。2.1.4邊界條件邊界條件有流量邊界與水頭邊界，在MODFLOW中，需要采用井來刻畫流量邊界。2.1.5初始條件數(shù)據(jù)格式同分層數(shù)據(jù)。216水位動態(tài)觀測資料數(shù)據(jù)格式同分層數(shù)據(jù)。3、VM軟件模擬地

11、下水流場:FlowField31新建工程項目1新建一個項目文件，輸入項目名字StfiSHiallMODFLiIW孕MEEifinWEJc-Dnjxit創(chuàng)“Setuptkti圖3.12.下一步，設置項目信息。圖3.23下一步，水流類型、初始參數(shù)的設置圖3.33下一步，導入底圖，初始網(wǎng)格剖分圖3.44下一步，模型坐標系的設置（將左圖中藍色的坐標系擴大，直到覆蓋整個研究區(qū)如右圖，見丸只口G目樂e詐口MpMDSTkl：3hiKixX-Mfi曲血耳NMA-PMn|lMM2HK|IUL.*U2IWi3.55下rH步,點擊ok,完成了項目文件的建立,同時進入了“模型輸入（Input）”模塊。圖3.6模型輸入

12、模塊界面3.2模塊I：模型輸入3.2.1網(wǎng)格設置3.2.11網(wǎng)格査看當你首次進入輸入模塊時，VisualMODFLOW會直接進入網(wǎng)格（Grid）輸入窗口。屏幕左邊的最上面的六個按鈕（ViewColumn、ViewRow、ViewLayer、Goto、Next和Previous），這幾個按鈕可以讓你從不同的角度查看現(xiàn)有的網(wǎng)格模型：查看橫剖面上點擊ViewRow，查看縱剖面點擊ViewColumn，查看平面點擊ViewLayer，Goto、Next和Previous主要是用于選擇層號。當VisualMODFLOW將模型的顯示方式從平面變?yōu)槠拭嬉晥D時，若網(wǎng)格的厚度與網(wǎng)格長度相比太小，層與層的界面不易

13、區(qū)分時，你必須對剖面圖進行垂向放大。在窗口底部選擇：單擊F8-VerticalExaggeration就會出現(xiàn)一個垂向放大窗口，供你輸入垂向放大系數(shù)。Type：40單擊OK3.2.1.2網(wǎng)格編輯下圖是網(wǎng)格編輯塊，可以通過“無效網(wǎng)格設置”來刻畫平面上邊界的不規(guī)則性；網(wǎng)格編輯模塊主要是用來調(diào)整網(wǎng)格的數(shù)量，包括行數(shù)、列數(shù)與層數(shù)；網(wǎng)格的伸縮主要是用來調(diào)整網(wǎng)格的大小。這里建立的網(wǎng)格還沒有賦實際高程，因此，可以通過ImportElevation、AssignElevation或Contouring來賦值。對于本實例主要用來練習“網(wǎng)格細化”功能與“網(wǎng)格插值”功能。GiidInactivaCaIIeEditG

14、rid跖itE:ctEni.呂Smoothiii.sr網(wǎng)格賦髙答值叱IeddftElm亍自tionAeejgnBlevstionlContourinsH圖3.7網(wǎng)格邊界模塊(1)網(wǎng)格細化為了提高模擬的精度，有時需要在某些地方將網(wǎng)格細化，比如井、河流、含水層天窗等。本節(jié)介紹了對感興趣的地段如供水井周圍地區(qū)、含水層天窗區(qū)與加油場之間地段進行模型網(wǎng)格細化的步驟。由于單元的大小代表了井的尺寸，因而用更細的網(wǎng)格來模擬井的影響將會更接近實際情況。另外，如果在井周圍出現(xiàn)了水位下降，用小的網(wǎng)格剖分可使水面更平滑。單擊EditColumns，然后將鼠標移至網(wǎng)格的任一處。注意在模型網(wǎng)格上有一個高亮度的豎線會隨鼠標

15、移動。這條豎線可以用來在模擬區(qū)的任何位置加上一列。在本練習中，我們通過按鼠標右鍵引出一個添加垂直網(wǎng)線(AddVerticalLine)窗口，以指定的間隔增加網(wǎng)格線。單擊Everyspacedgridlinesfrom：單擊其后的輸入框，輸入以下數(shù)值：from：500to：1600atintervalsof：25單擊OK接受這些值單擊Close離開(2)網(wǎng)格頂面與地面賦值在VisualMODFLOW中，你可以從Surfer的*.grd文件或以空格分隔的x、y、z的ASCII文件中輸入變化的層位標高。本例中，我們將輸入一個ASCII格式的x、y、z文件來生成一個有坡度的地表及層位，同時后者的厚度也

16、是變化的。單擊ImportElevation此時將出現(xiàn)如下窗口：圖3.7在layersurface中選擇地表高程Groundsurface,在option中選擇importdata,在interpolationsetting中的interpolation中選擇插值方法inversedistance，然后在interpolationsetting中的datasource中單擊打開，F(xiàn)cneat卍e4iE：vidin問E嚅dr.-inlapolahcflselhngsD-&I.-SSDLFG?Mi叭如JjNuiaiNeqttkousMnffuriLajsiIhickr&B(mJpi卻盼荻Stalj

17、OflIX-CDiOfClrH|iE圖3.8選擇地表高程txt”出現(xiàn)如下窗口圖3.9將X-coordinate后面的n/a改成對應Column號，這里改成1,類似地，Y-coordinate后面的n/a改成對應Column號，這里改成2,Elevation后面的n/a改成對應Column號，這里改成3。單擊Next,然后點擊finish，出現(xiàn)如下窗口，選擇world與meters。圖3.10單擊OK,點擊apply（不要點擊ok）,完成地表高程的賦值和插值，然后出現(xiàn)圖3.11上面窗口的右側(cè)便是給Groundsurface插值后的高程平面圖，點擊右上角的“3-DPreview”，然后調(diào)節(jié)窗口下R

18、otate的按鈕，顯示不同角度的三維圖，雀調(diào)整顯示圖的屬性石賣Laiillk|&mreBu*eT廠FiMdDij網(wǎng)It科i0*=|hpw：tejjHIsbAIIioarngBl3!NiKMME門*!0fllfdc-l9*ttonhppBFpEkpDIQFPPEk曰列.diS.cuidrii-scriMhHhEtah-JEnil9*mii2：i.mW1?勺艸“i：團崗彌1刼CS卿KSI0TIe-HlAIJSSJoSSt!；利iWnE加迪:i的鮒蚣C隧!如期陶嗣*民帑1翻馬卿1?MfcdKfflSHK-Mcznomw19S中的”対32汕迫迤*33阪*厲I即刖慟刪XS卿T1B43點ln14m：U(K

19、rna蟲4訓U妙1金H詢腔汕陽獨tlE勺旳“dll3L6HKA11嵐廉踽辭Qu-iiizmciiiscniiKZiKH?tai9s訓1ii1l41iSiijj*iGa,迫2夕42i樹嗣町這聞浮tit勺理也刖亍珅EiH加liZ?i卿133iJgglE和ENECS3JlijiKLS1：l-.1?*24圧HOea+3-i?、Rp!*i|EHi.mHRrttofiRniriamrUeg1!wj!b乜snUvnVwlrf-.rM-7rw+wiahiTtlH圖3.12點擊右上角的ArrayPreview”，顯示每個節(jié)點處的高程，當有部分點不合理時，可以在此對話框中直接修改。類似地，導入第一層底板、第二層底

20、板和第三層底板的高程。現(xiàn)在要導入第一層底板的高程：在layersurface中選擇BottomofLayer1在options中選擇Importdata在datasource中將“第一含水層地板高程.txt”的路徑載入。設置最小厚度MinimumLayerThickness0.75單擊Apply接受其它的默認參數(shù)值第一層的底板就會顯現(xiàn)出來。下一步要輸入第二層底板高程，從左側(cè)菜單上選擇：在layersurface中選擇BottomofLayer2在options中選擇Importdata在datasource中將“第二含水層地板高程.txt”的路徑載入。設置最小厚度MinimumLayerThi

21、ckness0.75單擊Apply接受其它的默認參數(shù)值第二層的底板就會顯現(xiàn)出來，它也具有一個變化的表面。接著輸入第三層底板高程。從左側(cè)菜單上選擇：在layersurface中選擇BottomofLayer3在options中選擇Importdata在datasource中將“第三含水層地板高程.txt”的路徑載入。設置最小厚度MinimumLayerThickness0.75單擊OK接受其它的默認參數(shù)值第三層的底板就會顯現(xiàn)出來，它具有一個變化的表面。你可以通過看網(wǎng)格模型的某一個剖面，此時，網(wǎng)格的表面是高低起伏的，如圖3.13。irriorfPwBiMJyLrrwwTI-txIFMITIMl2:

22、蠶*|匕;缶“慣吒氏曾悍m間灼港鷲|噩!faruqu.iHentsEonductivityStorage口LitialHeadVadoseZon.h淳透系數(shù)儲水丟數(shù)初始水頭包氣帶區(qū)域千密度離子參數(shù)BulkDcnsiiy初始沫度囂囂:Tsw一-參數(shù)編輯Initi：ilCutlcentfsltion.Dispersion_彌散度圖3.173231第一含水層（1）Conductivity單擊Properties-Conductivity,系統(tǒng)將要求你保存你的數(shù)據(jù)點擊Yes。J汚VisualMODFLOW-匚:vmodntTutorial.AirpiFileridWells|BoundariesPar

23、ticlesView匚olumnConductivity|ViewRowStarageViewLayerInitialHeadsVadoseZone.PrevioiGel門可以使用Assign-9Single、Assign-9Polygon或Assign-Windows來輸入具有不同的水力傳導系數(shù)的區(qū)域或?qū)游弧５谝缓畬樱媒o定的格式輸入如下2e-4HydraulicConductivityinx（m/sec）Kx：（X方向水力傳導系數(shù)Kx）HydraulicConductivityiny（m/sec）Ky：（Y方向水力傳導系數(shù)Ky）HydraulicConductivityinz（m/sec

24、）Kz：（Z方向水力傳導系數(shù)Kz）通過查看左邊的小立方體確定是在第一層，否則通過Next和Previous使目前的層號碼是第一層，點擊Assign-Windows,將整個平面選定后出現(xiàn)如下對話框2e-4（輸入Kx后自動出現(xiàn)）2e-4圖3.18圖3.19如果直接點擊OK,在該區(qū)域的參數(shù)為默認參數(shù)，由于現(xiàn)在含水層的參數(shù)與默認參數(shù)不同，因此,點擊New，然后將前面的參數(shù)輸入到對應的窗口中，如下圖圖3.19然后點擊OK,則第一層的滲透系數(shù)便賦好了。(2)Storage單擊PropertiesStorage,系統(tǒng)將要求你保存你的數(shù)據(jù)點擊Yes。可以使用AssignSingle、AssignPolygon

25、或AssignWindows來輸入具有不同的水力傳導系數(shù)的區(qū)域或?qū)游弧５谝缓畬樱媒o定的格式輸入如下值。TOC o 1-5 h zSpecificStorage(Ss)：1e-4(單位儲水系數(shù))SpecificYield(Sy)：01(單位給水度)Effectiveporosity(Eff.Por)：015(有效孔隙度)TotalPorosity(Tot.Por)015(總孔隙度)與滲透系數(shù)的賦值方式是一樣的，將上述參數(shù)賦值到第一層中，如下圖圖3.203232第二含水層(弱透水層)單擊GotoNext進入第二層。AssignWindows讓你可以在一個矩形窗口內(nèi)設置一個不同的水力傳導系數(shù)。現(xiàn)

26、在，讓我們給第二層弱透水層輸入一個較小的水力傳導系數(shù)。(1)Conductivity在左邊窗口中選擇，將Storage轉(zhuǎn)化為Conductivity。如下圖VLsualBODFLOf-D:tastv.rafEileridtcLlsEtctbtifipujidarjesParticlesJJStoragePompWells斗HeadObs.StorageInitialHk目UssConst.HeadRiversStream|Cueijiicllvlty圖3.20單擊AssignWindows把鼠標移到網(wǎng)格的西北角，在單元中間單擊。然后把鼠標移到東南角，點擊單元中部。這就形成了一個覆蓋整個層面的窗

27、口。出現(xiàn)一個屬性設置(AssignProperty)對話框。單擊New(整個網(wǎng)格將變成藍色)把鼠標移到Kx(m/s)框中單擊并輸入如下值：1e-8Ky(m/s)Kz(m/s)Kx(m/s)自動設置為1e-81e-9(2)Storage同理將第二層的，用給定的格式輸入如下值。TOC o 1-5 h zSpecificStorage(Ss)：1e-4SpecificYield(Sy)：01Effectiveporosity(Eff.Por)：0.15TotalPorosity(TotalPorosity)0.153233第三含水層現(xiàn)在把第一層的水力傳導系數(shù)與儲水系數(shù)等復制到第三層。GotoPrev

28、ious，回到第一層，在左邊的小窗口中選擇Conductivity，然后，在左邊菜單上選擇：單擊CopyLayer出現(xiàn)一個復制窗口，在窗口中:單擊單擊單擊如下圖Selectall/highlightedzones(復制所有的值)Layer3(它就會以咼亮度顯示)OK把第一層的K值復制到第三層圖3.21類似地，可以將Storge參數(shù)從第一層拷貝到第三層。3234天窗設置在特殊的沉積環(huán)境下，或者各種地質(zhì)構造的原因，會導致地層的缺失或者不連續(xù)。在這個特例中，含水層天窗區(qū)在底圖上就已有體然而，在許多情況下現(xiàn)是沒有這種資料的，必需依靠某些參數(shù)比如層厚來確定這個不連續(xù)的地段。假設這里層厚小于0.5處的位置

29、是一個巖性天窗，且天窗的巖性的巖性與第一層巖性相同。要顯示第二層的厚度可在下部菜單條上選擇F9-Overlay按鈕。這時就會出現(xiàn)一個圖層控制(OverlayControl)窗口(如下圖)，(Jfat丄毋心曲肚耶iyt口陽妙口ideiLUseiDeined擰D-aukSefTdsCl旳I刃LktFilv|AIAvaWjfe丁ILd:nnlpllwri?bhoveilrSdtdedQuerldyNahkF7P廠廠廠BC|FJ呼羊沖矗閱r占crtAWZH胡bK.CwveIb-ConLAr-dMlb-HwdGenArrtctjlnnCllCcfilpj-LabdsEllILaiMiTopHey.IL工

30、舊Th點IEk#EElffvBC|FJWetBC|Tj-kiNdaTia沖/BC|F|InxAhtFtDC|)WaIPTFafidEsBC|T：l-Fci-A5.：-5EeBC|TJ-Eprti-Bnipi-ainCone.BC|TJ-Rwhan?CogBC|FJ5VtiCtl$BC|FJDun圖3.22其中有一個表列，上面包含以字母順序排列的所有在設置模型輸入?yún)?shù)可能出現(xiàn)的圖層。向下瀏覽表列直到找厚度圖層(LayerThicknessOverlay)，然后雙擊它，在它的左邊就會出現(xiàn)一個“卜”表示它現(xiàn)在是一個活動圖層，這行的最后面“”是用來設置等厚圖等值線的。為了顯示的更加清楚，同時將Over

31、lay中其他選項都關閉，然后選擇ContourLabels。單擊OK將顯示層厚的等值線圖這里假設厚度小于1m的地方時個天窗。現(xiàn)在將等值線收縮到1m的最小厚度的地段放大，若沒有顯示0.8m的等值線時，可以在F9-Overlay中，點擊LayerThicknessOverlay后面的“”來設置。單擊F5-ZoomIn將鼠標移動到含水層天窗地段的左上部，單擊左鍵。在該區(qū)上拉出一個框，再單擊一次。屏幕上就會出現(xiàn)含水層天窗地段放大后的圖形。單擊AssignSingle(對單個單元進行設置)將會出現(xiàn)一個屬性設置框，其缺省值是最后輸入的K值(Zone#2)。按下New鍵下的向下箭頭。Zone#1就會是K的當

32、前值。將鼠標移到由1m等值線所圈定的表示含水層天窗的地段，按下鼠標，繞著該區(qū)拖動直到區(qū)內(nèi)單元格成為白色顯示。如果你已對區(qū)外的單元重新作了設置，那么，松開左鍵，按下右鍵就可以將單元恢復為它們原先的值的顏色。完成單元格的“描畫”后，現(xiàn)在你要將第一層的巖心賦值到這里，那么就Zone#中選擇2,單擊ok即可。3.2.4設置模型邊界條件：降雨補給，河流，邊界條件3241降雨補給要設置補給條件，你必須是在處于模型頂層的平面視圖情況下。檢查一下左下角的立方體導航器，看看你目前所處的層位，如果不在第1層，就用Next，Previous或Goto進入到第1層。我們首先要設置的邊界條件是含水層的補給通量。單擊Bo

33、undaries單擊Recharge如果你是正在建立新模型，此處會出現(xiàn)一個缺省補給(DefaultRecharge)窗口，提示你輸入補給量的初始值，VisualMODFLOW自動地將這個補給量值賦給模型的整個頂層，下面將下標的數(shù)據(jù)輸入軟件中。月份123456789101112降水補給1020304010010010010050403010在框中輸入：TOC o 1-5 h zStopTimedays：30Rechargemm/yr：10選中第一行，右擊選擇insertrow,在框中輸入:StopTimedays：60Rechargemm/yr:20以此類推，將所有數(shù)據(jù)都輸入軟件中，如下圖圖3.

34、23單擊OK接受賦值單擊F5-Zoom-In將光標移到加油場的左上角，單擊鼠標左鍵。沿加油場拉出一個框，然后再點一次閉合放大窗口。現(xiàn)在我們要給加油場設置一個較高的補給值：單擊AssignWindow將光標移到加油場的某個角上。單擊鼠標左鍵，沿加油場拖出一個窗口，然后再單擊鼠標左鍵。就會出現(xiàn)AssignRecharge窗口。單擊New設置了補給條件的窗口將變?yōu)镻roperty#2屬性的藍色，給Property#2輸入以下值：StopTimedays：360Rechargemm/yr：250單擊OK接受賦值單擊F6-Zoom-Out將模擬區(qū)恢復為滿屏顯示。3242邊界條件下一步是給承壓和非承壓含水

35、層的南北邊界設置定水頭邊界條件。（1）北部含水層邊界給Layer1和Layer2賦定水頭邊界，從頂部菜單上選擇:單擊Boundaries單擊ConstantHead單擊Yes（保存輸入的補給量數(shù)據(jù)）首先我們?yōu)樯喜繚撍畬友啬M區(qū)北部的邊界設置的邊界條件。單擊AssignLine用鼠標左鍵點擊西北角單元中心，然后用鼠標右鍵點擊東北角單元中心。此時單元格組成的水平線將以高亮度顯示，同時出現(xiàn)一個如圖所示的定水頭邊界輸入窗口圖3.24輸入以下數(shù)值：單擊StopTimeStopTime：360StartPoint：19EndPoint：19單擊OK接受輸入此時，粉色線將會變成暗紅色，表示常水頭邊界值已

36、經(jīng)設置完畢。從左邊菜單上選擇：單擊CopyLayer會出現(xiàn)一個復制窗口，選擇Constant-Head（O）。單擊Layer2（該層會高亮度顯示）單擊OK將定水頭邊界1號屬性復制到第2層給Layer3賦邊界，我們現(xiàn)在要輸入下部承壓含水層沿模擬區(qū)北部的邊界的邊界條件。單擊Goto,此時會有一個GotoLayer的彈出窗口。在Layeryouwishtogoto（你想去的層位）框中，輸入3。單擊OK進入第3層單擊AssignLine將鼠標移至網(wǎng)格的西北角。用鼠標左鍵單擊該單元中心。然后用鼠標右鍵單擊網(wǎng)格東北角的單元的中心。該線會以高亮度顯示，同時出現(xiàn)一個定水頭邊界輸入彈出窗口。輸入以下數(shù)值：單擊S

37、topTimeTOC o 1-5 h z HYPERLINK l bookmark174 StopTime：360 HYPERLINK l bookmark178 StartPoint：18 HYPERLINK l bookmark180 EndPoint：18單擊OK接受輸入（2）南部含水層邊界下一步我們要設置下部承壓含水層沿模擬區(qū)南部的邊界的邊界條件。單擊AssignLine將鼠標移至網(wǎng)格的西南角。用鼠標左鍵單擊該單元中心。然后用鼠標右鍵單擊東南角單元的中心。該線會以高亮度顯示，同時出現(xiàn)一個定水頭邊界輸入彈出窗口。單擊StopTimeTOC o 1-5 h z HYPERLINK l bo

38、okmark176 StopTime：360StartPoint：16.5EndPoint：16.5然后采用復制的方式，將第一層的邊界拷貝到第二層與第三層。單擊CopyLayer會出現(xiàn)一個復制窗口，選擇Constant-Head2）。單擊Layer2（該層會高亮度顯示）單擊Layer3（該層會高亮度顯示）單擊OK將定水頭邊界1號屬性復制到第2、3層。3.2.4.3河流下一步沿著模擬區(qū)南部邊界的頂層輸入河流邊界條件。單擊Boundaries單擊River切換到河流輸入窗口單擊AssignLine按照圖層布置，（從網(wǎng)格的西南部開始）沿著河流單擊鼠標左鍵使河流數(shù)字化。當?shù)竭_東南邊界時單擊鼠標右鍵。此

39、時會出現(xiàn)一個彈出窗口將讓你輸入河流的資料。圖3.24鍵入以下數(shù)據(jù)：選擇回AssigntoappropriateLayer與回lineargradient單擊StoppointTOC o 1-5 h zStopTime360StartPointRiverStage16.0StartPointRiverBottom15.5Conductance1000單擊EndpointEndPointRiverStage15.5EndPointRiverBottom15.0Conductance1000單擊OK接受這些值單擊F6-ZoomOut刷新屏幕定義河流后，將顯示一條用網(wǎng)格單元表示河流邊界條件的藍線，如下

40、圖所示圖3.253.3模塊II：計算（Run）到此為止，已經(jīng)將所有的數(shù)據(jù)導入MODFLOW軟件中了，下一步便是設置時間步長，計算不同時刻下水頭的空間分布值。331流動模擬運行選項點擊F10mainmenu回到主菜單，然后從主菜單上選擇：單擊Run你就會進入運行選項屏幕,(MpuTim?：I3.S68i,忖訕8弓pcilnJ1Ttsu：1Wiw-Cet4i3C-pia圖3.26MODFLOW2000的運行選項的設置如下：時間歩長彼芒初始水矣方用求解方法選擇t卜給設蜀層信息設理蔬干與回水地表水榛流暖總各項異性璇盤輸岀設置優(yōu)化列表圖3.27設置時間步長：點擊MODFLOW2000TimeSteps，

41、出現(xiàn)如下窗口Tik*p-LIHl.FiXlN行圧1禹I加f陽10jti3ltJl為1Jr圖3.28設置輸出水頭：點擊MODFLOW2000TOutputControl,出現(xiàn)如下窗口，并在輸出的時間上打鉤。71淘旳aasitcc-jiffii血刃WJIMh33？!?SKIWllKt-i期V阿SaTFlll3EHP?3WiStlXE：PPP0EtdiH4hFpnPrW幻LSIaIr廠廠丨廠廠=尸rCwRrw-1111111111Zrpmhj出廠廠廠廠廠廠廠廠廠廠廠廠廠ra|z圖3.29運行選項的設置如下:tHindkETI回邱曲BCTFmfTIEnqrEHui卜MODFLOW2D0：iPMODPA

42、TH廠ZcreudoetrMTHMGrFf5TITifeiefeBun轉(zhuǎn)化和運行IAdvflnwdWAulglml日ngiriE冷圖3.30新模型的默認初始條件是Constantbylayer（按層給定定水頭條件）。VisualMODFLOW將根據(jù)每層的標高和邊界條件為各層確定一個初始估計值。對于簡單的問題而言，這一般都會得到一個收斂的解。然而，一個好的估計值往往可以顯著地減少找到收斂值所需的迭代次數(shù)。注意，用于輸入初始水頭估計值的文件可以有多種格式可供選擇。因為我們要做的是個簡單問題，所以用Constantbylayer做初始水頭估計就已經(jīng)足夠了。3.4模塊m：輸出的可視化點擊output,

43、出現(xiàn)如下選項框圖3.313.4.1等勢線和等值線繪制選項要選擇等水頭線的等值線繪制選項：單擊Options（在左邊的菜單條上）就會出現(xiàn)一個如圖所示的等水頭線繪制選項（EquipotentialOverlayContouringOptions）窗口，但沒有右邊的顏色充填條。要激活顏色充填：單擊Colourshading在標明間隔（Interval：）的文本框中，將值從0.5改為0.25。在標著標注（Labels）的框中，將小數(shù)位數(shù)改為2。你也可以通過單擊標著等值線生成精度/速度（Contouringresolution/speed）的按鈕來改變等值線的繪制速度。每點一次，等值線的繪制速度增大2倍

44、。在本次練習中，按鈕上應該是High/Slowo單擊OK接受等值線的位置應該與前一個圖非常相似，但由于等值線的精度降低了，這些線會有一些輕微的擺動。上圖沒有加入顏色充填。其它沒有通過菜單激活的等值線選項可通過在模擬區(qū)內(nèi)按下鼠標右鍵來激活。此時會有一個帶有增添、刪除、移動等值線和標注的選項的彈出窗口（見下圖）。選擇AddContour選項，將鼠標移至模擬區(qū)的任意一處并單擊鼠標左鍵。在你鼠標點下的位置就會加入一條等值線。要再加入一條，只需再點一下鼠標即可。要再次彈出等值線選項窗口，只需簡單地再按一下鼠標右鍵即可。這一次我們選擇MoveLabel，然后將鼠標移至你想要移動的等值線標注的位置上。在標注

45、上按下鼠標不放，然后將其沿等值線拖至想放的位置，放開鼠標，這樣就設定了標注的新位置。標注都移動完后，再按一下鼠標右鍵。下面，讓我們看看場地的剖面情況。從左邊菜單上選擇ViewColumn，然后將鼠標移至模擬區(qū)內(nèi)。在模擬區(qū)中間位置按下鼠標左鍵選擇一列。要去掉用戶設置上去的等值線，在剖面上單擊鼠標右鍵彈出等值線繪制選項窗口。單擊Deleteallcustom單擊ViewLayer回到平面模式在進入下-節(jié)前，必須關掉色彩充填功能。在左側(cè)菜單條上單擊Options單擊Colourshading單擊OK3.4.2流速矢量及其選項要察看流速矢量選擇：單擊Velocities你將進入流速矢量輸出選項屏幕，見

46、下圖。流速矢量將根據(jù)其默認的設置進行繪制，即矢量相對大小與流速的大小一致。要讓流速矢量與大小無關單擊Direction（左邊菜單上）這樣所有的箭頭的大小都將是一樣的，它們只代表了流動方向。如果想改變顯示的矢量箭頭數(shù)目單擊Options（左邊菜單上）并輸入：Vectors：40這就指定了一行上的流速矢量的數(shù)目。回AutoscaleVariableScale單擊OK接受注意顏色：在默認情況下，紅色表示朝外（也就是說，從層的角度看是向上的）的流速，藍色表示向里（也就是說，從層的角度看是向下的），綠色表示與平面平行（也就是說，從層的角度看是水平的）。再讓我們看看場地的剖面情況。從左邊菜單上選擇View

47、Column,然后將鼠標移至模擬區(qū)內(nèi)。單擊鼠標左鍵，選擇任意一列。就會出現(xiàn)一個與下圖類似的在剖面上同時顯示等勢線和流速矢量的圖。顯示水力傳導系數(shù)圖層的離散化情況，從底部菜單欄中選擇F9-Overlay按鈕，將出現(xiàn)一圖層控制（OverlayControl）窗口，該窗口按字母順序列出所有可以打開或關閉的圖層。把鼠標置于ConductivityOverlay（水力傳導系數(shù)圖層）之上并雙擊鼠標左健。這樣就會在ConductivityOverlay旁出現(xiàn)一個星號（*）,意味著它被激活了。單擊OK顯示水力傳導系數(shù)圖層剖面。從左邊菜單選擇ViewLayer,單擊模型橫剖面中的第一層，返回到模型平面視圖。要從

48、屏幕中消除流速矢量，從底部菜單欄中選擇F9-Overlay按鈕，出現(xiàn)圖層控制窗口后，把鼠標置于VelocityOverlay上并單擊使之高亮度顯示，再單擊標有ON的按鈕使之變?yōu)镺FF。這樣VelocityOverlay旁的星號將會消失，意味著它已無效了。單擊OK顯示沒有流速矢量的窗口3.5模型校正對于數(shù)值模擬而言，模型校正是不可缺少的一部分，這里單獨拿出來練習，以加深影響。數(shù)值模擬調(diào)參數(shù)有許多方法，就目前而言，我們常采用的手法是人工手動調(diào)整法，即通過反復修改參數(shù)，來使計算值與實測值相對應。在運行結(jié)束后，進入Output中，點擊GraphsTimeSeriesHead。圖3.32出現(xiàn)下圖將L3.

49、Layer#3與W1.P1選中，點擊Apply,便在右圖中顯示出計算值曲線與實測值的散點圖，若這條線與散點吻合的較好，說明該模型的可靠性較強，否則需要修改含水層的參數(shù)，然后，在繼續(xù)運行，最后查看結(jié)果如何。反復此過程，知道兩者吻合較好為止。4、VM軟件模擬地下水流跡線：MODPATH4.1跡線模擬本節(jié)將一步一步地指導你輸入向前追蹤的質(zhì)點以確定污染物的途徑。點擊mainmenu,然后進入Input模塊。在頂部菜單欄上選擇：Particles此時程序會要求你保存數(shù)據(jù)：Yes我們在此要做的第一步是放大加油場地段。F5-ZoomIn在底部菜單欄在加油場的左上角單擊鼠標左鍵，然后在該區(qū)上拖出一窗口，再單擊

50、鼠標左鍵閉合放大窗口。AddLine把光標移到加油場的左邊并單擊鼠標左鍵。拖出一根線到加油場的右邊，然后再單擊一次。將顯示一添加質(zhì)點窗口。缺省的質(zhì)點數(shù)是10。把質(zhì)點數(shù)變成5。OK設置好加油場中的五個質(zhì)點的線通過加油站區(qū)的綠色質(zhì)點線表示是向前追蹤的質(zhì)點。現(xiàn)在返回到模擬區(qū)的滿屏顯示狀態(tài)。F6-ZoomOut點擊mainmenu,然后進入Run模塊選擇Pathlines后，會有一個如下的菜單出現(xiàn)Discharge與ReferenceTime兩個選項。Discharge設置控制質(zhì)點放置的選項，如果選擇Discharge,就會出現(xiàn)如下對話框這三個放置選項指明了邊界條件對質(zhì)點追蹤方法的影響。匚是一種將質(zhì)點從系統(tǒng)中移走的邊界條件。因而，抽水井是匯，第一層中具有蒸發(fā)作用的單元也是。因為MODPATH并不能知道質(zhì)點遇到弱匯時是否會被移走，所以它提供了以下三種選項：質(zhì)點通過弱匯單元，質(zhì)點在進入內(nèi)源單元后停止運動，當匯的排泄量超過單元入流量的某個特定百分比時，質(zhì)點停止運動。默認情況下該值為5%。ReferenceTime設置質(zhì)點參考時間的選項，如果選擇ReferenceTime,就會出現(xiàn)如下對話框SetRefezencBTineTimeformatrd.iStp/lel.Tm燈AbapliilevalueReference