VisualMODFLOW中文指南

1、Visual MODFLOW Pro 4.0 指南這個(gè)說(shuō)明書(shū)包含一步步的使用指南，指導(dǎo)你運(yùn)用Waterloo Hydrogeologic 的 Visual MODFLOWPro 4.0 Demo 版本的整個(gè)過(guò)程。這個(gè)指南將指引你按以下所要求的步驟去做：（ 1）建立網(wǎng)格模型，然后給定其特征值以及邊界條件；（2）使用Visual MODFLO怵模擬地下水的流動(dòng)，水質(zhì)點(diǎn)的示蹤，以及大量物質(zhì)的運(yùn)移模擬；（ 3）可以用二維和三維的可視化圖形顯示結(jié)果；為了方便起見(jiàn)，Visual MODFLOW Dem&裝程序自帶安裝一組完整的數(shù)據(jù)輸入文件，以及你將要建立的模型示例的模擬結(jié)果。這個(gè)使用指南獨(dú)立清晰地

2、顯示每個(gè)步驟，這可以讓你在調(diào)試過(guò)程中就能選擇感興趣的方面，而不用完成整個(gè)練習(xí)。對(duì)示例模型的描述該地區(qū)位于Waterloo 河外的飛機(jī)場(chǎng)附近。該地區(qū)表層的地質(zhì)條件是由上、下兩層沙土以及礫石組成的含水層，中間以粘土和泥沙作為隔水層組成的。相關(guān)的地貌特征包括一個(gè)飛機(jī)加油站，一個(gè)市內(nèi)的供水井區(qū)域，以及一個(gè)不連續(xù)的隔水層。如下圖所示：這個(gè)供水區(qū)有兩口供水井。東面的井以550m3/d 的定流量抽水，而西部的井以400m3/d 的定流量抽水。 在過(guò)去的十年里，飛機(jī)燃料已經(jīng)周期性地泄露到加油站區(qū)，而且經(jīng)過(guò)自然滲透導(dǎo)致一部分污染物進(jìn)入到上層含水層。這個(gè)指南將會(huì)指導(dǎo)你建立該地區(qū)的地下水流動(dòng)和溶質(zhì)運(yùn)移模型的必要步驟

3、。這個(gè)模型還會(huì)顯示燃料污染物對(duì)市內(nèi)供水井的潛在影響。當(dāng)研究該地區(qū)時(shí)，從平面圖上看，指定該區(qū)的上部為北面，下部為南，左面為西，右面為東。地下水在這個(gè)三層的含水層系統(tǒng)中是從北向南流動(dòng)（從上到下）。該含水層系統(tǒng)是由頂層潛水含水層、中間一個(gè)弱透水層以及底部的承壓含水層組成的，如下圖所示。上下含水層的滲透系數(shù)為2e-4m/sec ,而弱透水層的滲透系數(shù)為1e-10m/sec。怎樣使用這個(gè)指南這個(gè)指南分為四個(gè)模塊，每個(gè)模塊都包含一系列組成部分。設(shè)計(jì)這個(gè)指南是為了讓使用者可以從任何一個(gè)模塊開(kāi)始，這樣就能選用 Visual MODFLOW勺某些特定模塊。每個(gè)組成部分都以簡(jiǎn)潔明了的形式顯示。這些模塊如下設(shè)置：模

4、塊I :創(chuàng)建并定義一個(gè)流動(dòng)模型（第一至第五章）；模塊H:溶質(zhì)運(yùn)移模型（第六至第九章）；模塊田：運(yùn)行 Visual MODFLOW：第十至第H-）；模塊IV :可視化輸出；一、術(shù)語(yǔ)和符號(hào)為使用方便，會(huì)用到以下一些術(shù)語(yǔ)和符號(hào)：Type：輸入給定的詞或數(shù)據(jù)；Select ：在指示的地方單擊鼠標(biāo)左鍵；=:®<Tab><；<0 :按<£討$啾；:表示此鍵可以點(diǎn)擊，這個(gè)鍵會(huì)出現(xiàn)在一個(gè)窗口中，或者是在菜單欄的旁邊或底部;加粗的文字表示可以點(diǎn)擊的菜單或窗口界面，或者是可輸入文字。二、點(diǎn)擊開(kāi)始雙擊桌面上的Visual MODFLO州序圖標(biāo)開(kāi)始。，模塊I ：創(chuàng)建并

5、定義一個(gè)流動(dòng)模型第一部分：建立新模型這第一個(gè)模塊將會(huì)指導(dǎo)你運(yùn)用 Visual MODFLOM面創(chuàng)建一套新的模型數(shù)據(jù)所必要的步驟。創(chuàng)建新模型：在主菜單欄中單擊File （文件），再單擊New（新建），就會(huì)彈出一個(gè)Create new model （創(chuàng)建新模型）的窗口。為方便起見(jiàn)，系統(tǒng)按已經(jīng)建立了模型。按照默認(rèn)設(shè)置這個(gè)模型位于C:VMODNTTutorial建議使用者建立一個(gè)新的文件夾，把新模型保存在此新建文件夾中。建立新文件夾：?jiǎn)螕舸翱谥械男陆ㄎ募A圖標(biāo)，宜輸入新文件夾名稱，冉雙擊此文件夾，在 文件名稱欄中輸入Airport，點(diǎn)擊保存。注意：1、文件夾或者文件名稱必須用英文字母或數(shù)字，不能出現(xiàn)中

6、文。2、如果你選擇了缺省值的目錄途徑（C:VMODNTTutorial ）,并輸入相同的Airport ,就會(huì)彈出一個(gè)警告對(duì)話框，提示已經(jīng)存在一個(gè) Airport.vmf 的文件。單擊NO會(huì)把你新建的模型存入另外文件夾中；一口口 NTlTuitwidMi甲廿H w箱I印士匕CiejJe iww fBorjel單擊YES會(huì)覆蓋已經(jīng)存在的模型在以下三個(gè)連續(xù)的步驟描述模型建立過(guò)程:Visual 地下水流動(dòng)和運(yùn)移模擬所需的 Numeric Engine （數(shù)學(xué)計(jì)算軟件）；Units （單位） ；Start Date （初始日期）和Start Time （初始時(shí)間）對(duì)于流動(dòng)模型，可以用到以下的Numer

7、ic Engine ： WHI的 USGS MODFLOW-96 WHI的 USGS MODFLOW-2000 HGLB MODFLOW-SURFACT在這個(gè)示例中，選用MODFLOW-2000注意：在目前的Visual MODFLOWR本中，不支持 MODFLOW-SURFACT移選項(xiàng)。因此，如果在 水流模型中選擇了 MODFLOW-SURF ACT移選項(xiàng)的數(shù)字引擎就會(huì)自動(dòng)不可用了。對(duì)于運(yùn)移模型，用到了以下一些Numeric Engine （數(shù)字引擎）： MT3DMS MT3Dv150 MT3D96 MT3D99 MT3Dv1.0 MT3Dv2.5在這個(gè)示例中，選用了MT3DM。 S基于所選

8、的運(yùn)移引擎，你還要選擇可用的Sorption （吸收）方法和Reaction （反應(yīng)） 。在這個(gè)示例中，沒(méi)有吸收，也沒(méi)有選擇動(dòng)力反應(yīng)。在 Units （單位）欄中，為每個(gè)模型的輸入數(shù)據(jù)類(lèi)型選擇以下的信息內(nèi)容：Length ： metersTime:dayConductivity: m/secPumping Rate: m 3/dayRecharge:mm/yearMass:kilogramConcentration: milligrams/liter這個(gè)模型的Start Date （初始日期）和Start Time （初始時(shí)間）與初始模擬時(shí)間相對(duì)應(yīng)。當(dāng)前，這些數(shù)據(jù)只和非穩(wěn)定流模擬有關(guān)。在非穩(wěn)定

9、流模擬中，所記錄的數(shù)據(jù)可以被導(dǎo)入用來(lái)為選定的邊界條件（如定水頭、河流、總水頭和排放量）設(shè)定時(shí)間進(jìn)度表。再點(diǎn)擊NEXT（下一步）確定，確定這些給定值。在第二步中，會(huì)出現(xiàn)下面的窗口，顯示地下水流動(dòng)和運(yùn)移模型的默認(rèn)值:再點(diǎn)擊NEXT（下一步）確定這些給定值。第三步是Create the model Grid（建立網(wǎng)格模型）見(jiàn)下圖：以上的界面用來(lái)Import a site map （添加地圖），指定Model Domain （模型版圖）的維數(shù),并且確定有限差分網(wǎng)格的行數(shù)、Columns (j): 40Rows (i):40Layers (k): 3Zmin:0ZZmax:18Z選擇 Import a

10、site map列數(shù)以及層數(shù)。為模型中的行數(shù)、列數(shù)：40行數(shù)：40層數(shù)：3方向最小值：0方向最大值：18（添加地圖）列數(shù)以及層數(shù)分別輸入以下數(shù)據(jù):下一步，必須要選擇文件擴(kuò)展名為 .DXF的背景地圖點(diǎn)擊Browse（瀏覽），返回到Tutorial目錄選擇以下文件:單擊"Sitemap.dxf單擊open單擊Finish確定這些設(shè)定。這時(shí)會(huì)彈出一個(gè)Select Model Region（選擇模型區(qū)域）的窗口，提示你定義模型區(qū)域的范圍。Visual MODFLO悵從底圖（Sitemap.dxf ）中讀取最大最小坐標(biāo)值，并顯示其在模型中心的默認(rèn)位置坐標(biāo)。在屏幕顯示數(shù)據(jù)的地方重新輸入以下數(shù)據(jù)：

11、Display area : X1 ： 0<>（顯示范圍）Y1:0=>X2:2000=>Y2:2000=>Model Origin : X：0（初始模型）Y:0=>Angel（角度）：0 gModel Corners : X1 :0（模型頂角）YI :0X22000Y22000單擊OK確定。此時(shí)會(huì)出現(xiàn)一個(gè)Fileattributes的窗口，顯示“ Sitemap.dxf ”已經(jīng)被保存到Visual MODFLOW工程里，并且命名為0單擊OK就會(huì)打開(kāi)Input （輸入）菜單，而且在模型范圍內(nèi)會(huì)自動(dòng)生成一個(gè)固定規(guī)格的40X 40X 3的有限差分網(wǎng)格。屏幕上還會(huì)出

12、現(xiàn)地圖的網(wǎng)格形式，如下圖所示：當(dāng)?shù)谝淮芜M(jìn)入Input （輸入）時(shí)，網(wǎng)格圖自動(dòng)裝載。第二部分：網(wǎng)格加密網(wǎng)格界面為加密模型網(wǎng)格、勾畫(huà)無(wú)效網(wǎng)格單元、輸入各層標(biāo)高、指定標(biāo)高、優(yōu)化（光滑）網(wǎng)格 以及畫(huà)各層標(biāo)高等值線提供了一套完整的繪圖工具。這個(gè)部分描述了對(duì)所選地區(qū)進(jìn)行網(wǎng)格加密的必要步驟， 例如在供水井周?chē)貐^(qū)。進(jìn)行網(wǎng)格加密 的原因是為了得到所選地區(qū)更詳細(xì)的模擬結(jié)果，（如在水力梯度可能變化大的地區(qū)）。例如：如果在 井周?chē)霈F(xiàn)水位下降，當(dāng)使用了加密的網(wǎng)格后，水位會(huì)相對(duì)較平緩。修改X方向的網(wǎng)格：?jiǎn)螕鬍dit Grid Edit Columns（編輯列）就會(huì)彈出一個(gè)Columns（列）的窗口，顯示要編輯的網(wǎng)格列

13、的選項(xiàng)。Add （添加）選項(xiàng)自動(dòng)被選中， 同時(shí)通過(guò)在需要增加網(wǎng)格線的地方單擊鼠標(biāo)左鍵，就會(huì)增加一條新的線。把鼠標(biāo)移到網(wǎng)格中任意地方，注意到有一條突出的垂線會(huì)隨著鼠標(biāo)在網(wǎng)格中移動(dòng)。這條線可以被用來(lái)在模型范圍內(nèi)任何地方增加一列。在這次的練習(xí)中，需要修改加油站區(qū)以及供水井附近的網(wǎng)格。在網(wǎng)格的任意位置點(diǎn)擊鼠標(biāo) 右鍵，就會(huì)彈出一個(gè)Add Vertical Line（增加垂線）的窗口，選擇。Evenly spaced gridlines from:在調(diào)整文本框中點(diǎn)擊，并輸入以下數(shù)據(jù)：from: 500to :1600at interval of: 25點(diǎn)擊OK確定，冉點(diǎn)擊CLOSE關(guān)閉Columns （列

14、）的窗口。下一步，要修改加油站到供水井區(qū)域 Y方向的網(wǎng)格。單擊Edit Grid Edit Rows （編輯行）在網(wǎng)格的任意位置點(diǎn)擊鼠標(biāo)右鍵，就會(huì)彈出一個(gè)Add Horizontal Line （增加水平線）的窗口，選擇。Evenly spaced gridlines from:在調(diào)整文本框中點(diǎn)擊，并輸入以下數(shù)據(jù)：from: 400 '一to :1900at interval of: 25點(diǎn)擊OK確定，冉點(diǎn)擊CLOSE關(guān)閉Rows（行）的窗口，退出。修改后的網(wǎng)格應(yīng)該如下圖所示：接下去的步驟會(huì)教你怎樣查看模型的橫截面以及怎樣給模型增加新的層數(shù)。查看模型的橫截面：點(diǎn)擊左邊工具欄中的View

15、 Column（查看列），把光標(biāo)移到網(wǎng)格中的任意位置。當(dāng)在屏幕上移動(dòng)光標(biāo)時(shí)，有一條紅色的豎杠會(huì)使某一列突出顯示， 它會(huì)隨著鼠標(biāo)移動(dòng)。要選擇一列查看，在選中的一列上單擊鼠標(biāo)左鍵，Visual MODFLOW會(huì)把網(wǎng)格屏幕從平面圖切換到橫截面圖顯示。此時(shí)，模型由于沒(méi)有縱向擴(kuò)展，橫截面圖中的三層以接近一條粗線的形式顯示，幾乎不能分辨出這三層。為了比較合理地顯示這三層，需要縱向放大橫截面：?jiǎn)螕羝聊坏撞康腇8-Vert Exag（縱向放大），就會(huì)彈出一個(gè) Vertical Exaggeration 的窗口，提示輸入縱向放大的值： 輸入：40點(diǎn)擊OK模型的三個(gè)層面就會(huì)如下所示：從以上的數(shù)據(jù)可以看出，在整個(gè)橫

16、截面每層厚度都相同。然而自然條件下幾乎不可能會(huì)有像這樣水 平，并且還具有相同厚度的地質(zhì)條件的含水層。在這個(gè)示例中，可以在模型區(qū)域中輸入包含標(biāo)高的，有確定的X和Y坐標(biāo)的文本文件。單擊左邊工具欄中的Import Elevation（輸入標(biāo)高），就會(huì)彈出一個(gè)Create gridelevation （創(chuàng)建網(wǎng)格標(biāo)高）的窗口：這個(gè)窗口要求為網(wǎng)格模型各層輸入一組數(shù)據(jù)。在Layer surface欄中自動(dòng)選擇了Groundsurface ,在窗口的右側(cè)顯示2-D和3-D（2維和3維）的插入層表面標(biāo)高，以及逐格的展開(kāi)層 數(shù)據(jù)Array （矩陣）。點(diǎn)擊Option （選項(xiàng)）的下拉菜單三,并選擇Import da

17、te。當(dāng)選定這個(gè)選項(xiàng)后，就會(huì)在Option （選項(xiàng)）下面出現(xiàn)Interpolation setting（差分組合框）。這個(gè)組合框用于選擇一個(gè)數(shù)據(jù)庫(kù)文件以及差分法（Natural Neighbors, Kriging, Inverse Distance ）。Interpolation 欄中自動(dòng)選擇了 Natural Neighbors 。然而這個(gè)示例需要在下拉菜單中選擇Inverse Distance在Date source （數(shù)據(jù)來(lái)源）欄中點(diǎn)擊打開(kāi)文件夾的圖標(biāo),就會(huì)出現(xiàn)如下的Open （打開(kāi)）窗口:選擇地表的頂層地層數(shù)據(jù)文件:單擊 Airport-gs.asc再單擊Open就會(huì)出現(xiàn)如下的Mat

18、ch field窗口：Required Data （所需數(shù)據(jù)）框列出了來(lái)自數(shù)據(jù)庫(kù)的縱行數(shù)據(jù)，而且根據(jù)數(shù)值差分法，Match tocolumn number （和列序號(hào)匹配）框用來(lái)把數(shù)據(jù)和各自相應(yīng)的數(shù)據(jù)庫(kù)相連接。如上窗口所示，在 每欄中輸入合適的數(shù)字，把 X坐標(biāo)與Column#1 （第一列）匹配，Y坐標(biāo)與Column#2 （第二列） 匹配，日evation （標(biāo)高）與Column#3 （第三列）匹配。一旦這些列都匹配好了，NEXT按扭就會(huì)被激活。單擊NEXT進(jìn)入Data Validation（有效數(shù)據(jù)）界面。如下圖所示：由于沒(méi)有發(fā)現(xiàn)錯(cuò)誤輸入（錯(cuò)誤的輸入會(huì)以紅色的字體顯示），就確保了數(shù)據(jù)庫(kù)中的數(shù)據(jù)

19、都是有效的， 而且數(shù)據(jù)都存在于模型范圍內(nèi)，如下圖的 Coordinate System and Uuits （坐標(biāo)系和單位） 窗口所描述的。單擊Finish 進(jìn)入Coordinate System and Uuits（坐標(biāo)系和單位）窗口，然后在Coordinate System （坐標(biāo)系）框中選擇 Model （標(biāo)準(zhǔn)）,在Elevation Units （標(biāo)高單位） 框中選擇Meter （米）。如下圖所示：?jiǎn)螕鬙K對(duì)輸入的標(biāo)高數(shù)據(jù)進(jìn)行插值，就會(huì)出現(xiàn)如下的窗口 ：再單擊OK在進(jìn)入下一步操作，把差分層的標(biāo)高分配給選中的網(wǎng)格模型層之前，會(huì)出現(xiàn)一個(gè)Warning （警告）框，確認(rèn)標(biāo)準(zhǔn)的標(biāo)高是否確實(shí)需要

20、更改。單擊YES如果所提出的對(duì)模型地表標(biāo)高的更改， 會(huì)和之前指定的邊界條件上的水頭值有沖突，就會(huì)出現(xiàn)一條警告信息，提供修改數(shù)據(jù)的機(jī)會(huì)。在這種情況下，不會(huì)發(fā)現(xiàn)潛在的問(wèn)題。下一步，將要給第一層輸入底板標(biāo)高：?jiǎn)螕鬒mport Elevation（輸入標(biāo)高）；在Layer surface 欄的下拉菜單中 三I選才B Bottom of Layer 1；在Option的下拉菜單中"3選才? Import data ；并在其組合框中選擇Inverse Distance 的 插分法；.一一一一 慮I在Date source （數(shù)據(jù)來(lái)源）欄中點(diǎn)擊打開(kāi)文件夾的圖標(biāo)，燈并選擇Airpt-b1.asc 。

21、在Open （打開(kāi)）窗口中選擇第一層的標(biāo)高。再單擊Open。把X坐標(biāo)與Column#1 （第一列）匹配，Y坐標(biāo)與Column#2 （第二列）匹配，日evation （標(biāo)高）與Column#3 （第三列）匹配。單擊NEXT進(jìn)入Data Validation （有效數(shù)據(jù)）界面。單擊Finish 進(jìn)入 Coordinate System and Units（坐標(biāo)系和單位）窗口，然后在Coordinate System （坐標(biāo)系）框中選擇 Model （標(biāo)準(zhǔn)），在Elevation Units（標(biāo)高單位）框中選擇Meter （米）。單擊OK對(duì)輸入的標(biāo)高數(shù)據(jù)進(jìn)行插分，在Create grid eleva

22、tion（建立網(wǎng)格標(biāo)高）的窗口再單擊OK;單擊YES確認(rèn)對(duì)模型標(biāo)高的更改這時(shí)會(huì)發(fā)現(xiàn)第一層的底部會(huì)呈現(xiàn)不同的標(biāo)高。 接下去再按此步驟完成對(duì)第二層和第三層標(biāo)高的 更改：?jiǎn)螕鬒mport Elevation（輸入標(biāo)高）；在Layer surface 欄的下拉菜單中 三I選才B Bottom of Layer 2；在Option的下拉菜單中三1選才? Import data ；并在其組合框中選擇Inverse Distance 的 插分法 、小一，一，局在Date source （數(shù)據(jù)來(lái)源）欄中點(diǎn)擊打開(kāi)文件夾的圖標(biāo)，5并選擇Airpt-b2.asc 。在Open （打開(kāi)）窗口中選擇第二層的標(biāo)高。再單擊

23、Open。把X坐標(biāo)與Column#1 （第一列）匹配，Y坐標(biāo)與Column#2 （第二列）匹配，日evation （標(biāo) 高）與Column#3 （第三列）匹配。單擊NEXT進(jìn)入Data Validation（有效數(shù)據(jù)）界面。單擊Finish 進(jìn)入Coordinate System and Units（坐標(biāo)系和單位）窗口，然后在Coordinate System （坐標(biāo)系）框中選擇 Model （標(biāo)準(zhǔn)），在Elevation Units（標(biāo)高單位）框中選擇Meter （米）。單擊OK對(duì)輸入的標(biāo)高數(shù)據(jù)進(jìn)行插分，在Create grid elevation（建立網(wǎng)格標(biāo)高）的窗口再單擊OK;單擊YES

24、確認(rèn)對(duì)模型標(biāo)高的更改。這時(shí)會(huì)發(fā)現(xiàn)第二層的底部會(huì)呈現(xiàn)不同的標(biāo)高。接下去再按以上步驟完成對(duì)第三層表面標(biāo)高的更改。要注意對(duì)Bottom of Layer 3來(lái)說(shuō)要在下拉菜單中 三1選才| Airpt-b3.asc 。第三層的底部同樣也呈現(xiàn)不同的標(biāo)高。模型的橫截面應(yīng)如下圖所示（圖為第 38行橫截面）：為了得到一個(gè)更好的縱向模型網(wǎng)格描繪，還可以對(duì)每層進(jìn)行分層。在左邊的工具欄中單擊Edit GridEdit Layer（編輯層）就會(huì)彈出一個(gè)Layer （層）的窗口，單擊 ORefined by ,在輸入框中輸入2,把鼠標(biāo)移到模型的橫截面中，就會(huì)高亮顯示該層表面。把鼠標(biāo)放在與第一層頂部相對(duì)應(yīng)的線上，在那里左

25、鍵單擊；然后把鼠標(biāo)移到該層相對(duì)應(yīng)的底線上，在那里再單擊左鍵。這樣在第一層垂直方向 上的中心位置就會(huì)增加一個(gè)已經(jīng)變形的層，即平均分為兩層。重復(fù)此操作對(duì)第二層（隔水層）和第三層（下部承壓含水層）也進(jìn)行分層。當(dāng)完成所有這些操作時(shí)，單擊 Layer （層）窗口中的Close退出。模型的橫截面現(xiàn)在就由六層組成，而且應(yīng)該和下圖所示的相似（圖中為第38行）：接下去要返回到模型的層界面，在左邊的工具欄中點(diǎn)擊View Layer當(dāng)鼠標(biāo)放在第一層，并顯示不同顏色時(shí)就在該層單擊左鍵。這樣就會(huì)出現(xiàn) Airport地區(qū)的平面圖。第三部分：添加抽水井這個(gè)部分的目的是指導(dǎo)給模型添加抽水井的必要步驟。在主菜單中單擊Well在

26、其下拉菜單中F單擊Pumping Wells這時(shí)會(huì)出現(xiàn)要求保存資料的窗口，點(diǎn)擊YES保存，并繼續(xù)。一旦模型保存好了，屏幕就會(huì)切換到 Pump Well （抽水井）的界面。注意到左邊的工具欄中的按 扭現(xiàn)在是變?yōu)榫x項(xiàng)，如添加、刪除、編輯、移動(dòng)以及拷貝抽水井。在為這個(gè)模型添加抽水井之前， 對(duì)抽水井周?chē)貐^(qū)進(jìn)行放大（抽水井位于模型范圍的右下角）o單擊F5-Zoom In把鼠標(biāo)移到抽水井的右上角，左鍵單擊。然后覆蓋抽水井所在區(qū)域，拉出一個(gè)記號(hào)框，再單擊左鍵。 所選區(qū)域就被放大了。接著給模型添加抽水井：口口單擊Add Well把鼠標(biāo)移到有抽水井標(biāo)記的西面，然后在那里左鍵單擊，就添加了一個(gè)抽水井。此時(shí)會(huì)彈

27、出一個(gè) NeWWell（新建井）的窗口，提示輸入指定的資料。輸入以下信息：Well Name :Supply Well 1X : 1415Y : 535要增加 Screened Interval ，點(diǎn)擊 Screen Bottom 所在列，并輸入以下值：Screen Bottom (m) : 0.3 Screen Top (m) : 5.0注意到井的過(guò)濾器長(zhǎng)度在窗口右邊的井孔圖表中顯示出來(lái)。井的過(guò)濾器長(zhǎng)度可以通過(guò)點(diǎn)擊并拖動(dòng)過(guò)濾器邊緣到新的標(biāo)高來(lái)改變。要輸入井的Pumping Schedule (抽水時(shí)間)，在End(day)所在列下面的文本框左鍵單擊鼠標(biāo)，然后輸入以下信息：End (day)

28、:7300Rate (mA3/day) : -400單擊 OK 確定。選用 7300天作為結(jié)束時(shí)間是因?yàn)檫@個(gè)模擬將要進(jìn)行20年。然而， 由于這是穩(wěn)定流的流動(dòng)模擬，可以任意確定模擬時(shí)間，但這個(gè)水流模擬的結(jié)果不會(huì)改變。穩(wěn)定流模型的意思是這個(gè)模型給每個(gè)參數(shù)賦上第一次輸入的值來(lái)進(jìn)行模擬，直到達(dá)到平衡，而不管相對(duì)時(shí)間如何改變。抽水量為負(fù)值是因?yàn)檫@是一個(gè)抽水井。如果在輸入所需數(shù)據(jù)時(shí)失敗了，Visual MODFLO悵提示再次完成表格。下一步是用一條捷徑給第二個(gè)抽水井設(shè)定參數(shù)。這條捷徑會(huì)把這個(gè)井的特征全都拷貝到另一個(gè)井。在左邊的工具欄中點(diǎn)擊Copy Well把鼠標(biāo)箭頭指在西邊的井上(供水井1) ，并左鍵單擊

29、。然后把鼠標(biāo)移到東面的井上，再單擊左鍵，就拷貝了一個(gè)井。下一步要對(duì)新(拷貝)井進(jìn)行編輯。在左邊的工具欄中點(diǎn)擊Edit Well ，并單擊新井，就會(huì)出現(xiàn)一個(gè)Edit Well 的窗口。當(dāng)前該井的命名為“Supply Well 1 ( 2) ”。把 Well Name 對(duì)應(yīng)欄中井的名稱改為“ Supply Well 2 ”，并輸入： X=1463， Y=509， Rate (mA3/day) : -550單擊 OK 確定。單擊 F6-Zoom out第四部分：賦水流模型屬性這個(gè)部分是指導(dǎo)建立一個(gè)模型，具各層滲透系數(shù)相差很大的必要步驟在主菜單上點(diǎn)擊Properties/Conductivity（參數(shù)

30、/滲透系數(shù)）在彈出的窗口中點(diǎn)擊YES保存資料?，F(xiàn)在屏幕就切換到Conductivity （滲透系數(shù)）的輸入界面，在那里可以對(duì)模型的滲透系數(shù)（ Kx, Ky, Kz）進(jìn)行修改。在左邊的工具欄中點(diǎn)擊Database按扭，并在彈出窗口中輸入以下值:Kx （m/s）: 2e-4Ky （m/s）: 2e-4Kz （m/s）: 2e-4單擊OK確定。要注意此時(shí)的Kx, Ky, Kz值都相同，這表示假設(shè)給定的值在水平和垂直方向上是各向同性的。但 也可以通過(guò)修改滲透系數(shù)給模型設(shè)定各向異性的特征值。在這個(gè)六層的模型中，第一、第二層表示上層含水層；第五、第六層表示下層含水層；第三、第四 層表示分隔上下含水層的弱透

31、水層。在這個(gè)示例中，把之前設(shè)定的滲透系數(shù)用于第一、二、五、六 層（含水層），同時(shí)給第三、四層（弱透水層）設(shè)定不同的滲透系數(shù)。要注意到第一層是模型的頂在左邊的工具欄中單擊Goto,就會(huì)彈出一個(gè)Go To Layer的窗口,在窗口中輸入3,單擊OK確定?，F(xiàn)在就在第三層（弱透水層）的界面上。下一步是給弱透水層（第三、四層）設(shè)定較小的滲透系數(shù) 值?？梢栽谀Ｐ途W(wǎng)格上指定參數(shù)值。在左邊的工具欄點(diǎn)擊Assign,再單擊Window,這個(gè)功能可以給指定的矩形框內(nèi)賦值一個(gè) 不同的滲透系數(shù)。把鼠標(biāo)移到網(wǎng)格的西北角，并在頂角的網(wǎng)格中點(diǎn)擊左鍵。然后把鼠標(biāo)移到網(wǎng)格的東南角， 并在頂角 的網(wǎng)格中點(diǎn)擊左鍵，此時(shí)就創(chuàng)建了一個(gè)

32、白色的界面覆蓋整個(gè)層面。此時(shí)會(huì)出現(xiàn)一個(gè) Conductivity-Assign Window 的對(duì)話框。單擊New, 這個(gè)網(wǎng)格就會(huì)變成藍(lán)色的，而且 Zone#®變?yōu)?對(duì)新的Zone#2輸入滲透系數(shù)，單擊Kx對(duì)應(yīng)欄，并輸入以下值:Kx (m/s):1e-10Ky (m/s)：1e-10Kz (m/s)：1e-11單擊OK確定?，F(xiàn)在把第三層的滲透系數(shù)拷貝給第四層：在左邊的工具欄中單擊CopyLayer再單擊Window,就會(huì)彈出一個(gè)Copy Layer的窗口，選才S Select all/highlighted zone(注意到Zone部分的Zone#2會(huì)自動(dòng)選中)在target lay

33、er (目標(biāo)層)中選擇Layer 4 (第四層)單擊OK,就把Zone#2中第三層的滲透系數(shù)值拷貝到第四層了。現(xiàn)在來(lái)觀察模型的橫截面，看三個(gè)水文地質(zhì)單元。在左邊的工具欄中點(diǎn)擊View Column,并在不連續(xù)的隔水層區(qū)點(diǎn)擊其中一列。上、下含水層以及隔水層由于導(dǎo)水特性不同而呈現(xiàn)不同的顏色，這可以在平面圖上顯示出來(lái)，如上圖所示：(圖為第36列)。盡管弱透水層在趨于尖滅的地方很薄，但這些地區(qū)的滲透系數(shù)也要和上、下含水層的滲透系數(shù)一樣 設(shè)置。要做到這一步，必須返回到模型的層界面：在左邊的工具欄點(diǎn)擊View Layer,并在圖中第三層處單擊左鍵。在這個(gè)特殊的示例中，不連續(xù)隔水層范圍在 DXF圖中已經(jīng)有顯

34、示。但是，在很多情況下，這在地圖上是沒(méi)有標(biāo)注的，必須依靠其他標(biāo)志：如根據(jù)層的厚度確定水文地質(zhì)單元尖滅的區(qū)域。下面要顯示第三層的厚度：在底部的工具欄中單擊F9-Overlay按扭，就會(huì)彈出一個(gè)如下所示的 Overlay Control 的窗 口，窗口中包含一系列有用的疊加層。當(dāng)在模型輸入?yún)?shù)時(shí)，這些疊加層就可見(jiàn)了。選中C(I)-Layer Thickness單擊OK顯示該層厚度的輪廓圖。注意：如果由于背景顏色的關(guān)系，看不到輪廓線，可按以下步驟操作：?jiǎn)螕鬎9-Overlay按扭選中團(tuán)C（I）-Layer Thickness,并單擊這一行右邊的按扭，就會(huì)彈出一個(gè) C（I）-LayerThicknes

35、s的選項(xiàng)框，在這里可以改變輪廓的設(shè)置?？梢試L試把線加粗（增加寬度）,選擇不同的顏色，加大標(biāo)簽大?。ɡ?3）。下一步，放大等值線為0.5厚的區(qū)域，單擊F5-Zoom in把鼠標(biāo)移到不連續(xù)隔水層的左上角，單擊鼠標(biāo)左鍵，然后在所選區(qū)拉出一個(gè) 框，再單擊左鍵。屏幕上就會(huì)出現(xiàn)不連續(xù)隔水層的放大圖形。單擊Assign,選擇Single對(duì)每個(gè)網(wǎng)格都設(shè)定特征值。此時(shí)會(huì)彈出一個(gè) Conductivity-Assign Single的對(duì)話框，顯示在（Zone#2）輸入的最后一個(gè)K值。點(diǎn)擊Zone冊(cè)向下箭頭，Zone#1的K值就選中了。此時(shí)先不要選擇OK把鼠標(biāo)移到標(biāo)有0.5米輪廓線的不連續(xù)隔水層處，按住鼠標(biāo)左鍵并

36、在0.5米的區(qū)域內(nèi)拖動(dòng)鼠標(biāo)，直到該區(qū)域內(nèi)的網(wǎng)格都變成了白色（如下圖所示）：如果不小心選擇了一些不應(yīng)包括在內(nèi)的網(wǎng)格，在這些網(wǎng)格上點(diǎn)擊右鍵把它們變回原來(lái)的設(shè)置。一旦 完成描繪這些網(wǎng)格，在 Conductivity-Assign Window 的對(duì)話框中點(diǎn)擊OK。下一步要把第三層的導(dǎo)水區(qū)域拷貝到第四層：在左邊的工具欄中單擊CopyLayer ,就會(huì)彈出一個(gè)Copy Layer的窗口，列出了當(dāng)前層中可 用的導(dǎo)水區(qū)（Zone#1和Zone#2）,以及要被拷貝的目標(biāo)區(qū)。選擇 0Select all/highlighted zone（注意到 Zone#1 和 Zone#2 都自動(dòng)被選中） 在target l

37、ayer （目標(biāo)層）中選擇 0Layer 4 （第四層），單擊OK,就把Zone#2中第三層的導(dǎo)水帶拷貝到第四層了。單擊F6-Zoom out返回到整個(gè)模型的全屏顯示，要觀察模型的橫截面。單擊View Column在不連續(xù)隔水層地帶選擇一列并點(diǎn)擊左鍵（如第 31歹1）。就會(huì)出現(xiàn)如下圖所示的一個(gè)橫截面圖：要返回平面圖可以點(diǎn)擊View Layer按扭，然后在橫截面圖中左鍵單擊第一層。接下去要修改 貯水率（釋水率）。在主菜單上點(diǎn)擊Properties/Storage （ 參數(shù) /釋水率 ） ，就會(huì)切換到Storage （ 釋水率）界面，此處只能修改默認(rèn)的貯水率值，但不能更改給水度、有效空隙度或總空隙

38、度的值。單擊Database,然后輸入貯水率Ss （1/m） =1E-4 （如下圖所示）：?jiǎn)螕?OK 確定。該層的厚度輪廓仍然可見(jiàn)，要去掉這些輪廓：?jiǎn)螕?F9-Overlay 按扭 ， 取消 C（I）-Layer Thickness單擊 OK 確定，此時(shí)該層的厚度輪廓線就取消了?，F(xiàn)在就完成了使用指南中如何給區(qū)賦恒定值的部分。在這個(gè)示例中，已經(jīng)修改了模型的貯水率和滲透系數(shù)。Visual MODFLOW還提供了相同的繪圖工具為水流模型編輯初始水頭，以及和水流及運(yùn)移過(guò)程相關(guān)的參數(shù)。第五部分：給定水流邊界條件下面這部分是描述給定不同模型邊界條件的所需步驟。一、確定模型含水層的補(bǔ)給條件在大多數(shù)情況下，含

39、水層都是通過(guò)地表水入滲補(bǔ)給的。在Visual MODFLOW為了確定補(bǔ)給，必須要到模型頂層界面。檢查屏幕左下角的導(dǎo)航立方體，可以看到目前所在的層位。第一個(gè)要確定的邊界條件是含水層的補(bǔ)給量。單擊主菜單上的Boundaries/Recharge如果正在對(duì)一個(gè)新的模型進(jìn)行操作，此時(shí)就會(huì)彈出一個(gè)Recharge-Assign Default Boundary的窗口，提示輸入Stop Time （停止時(shí)間）。要注意當(dāng)建立模型后，給定的默認(rèn)補(bǔ)給量為0。VisualMODFLOW劭給模型的整個(gè)頂層賦予這個(gè)值。在彈出框中輸入：Stop Time day : 7300Recharge mm/yr : 100單擊

40、 OK 確定。單擊 F5-Zoom in 把鼠標(biāo)移到加油站的左上角，并單擊左鍵。然后在加油站區(qū)拉出一個(gè)框，再單擊左鍵關(guān)閉放大窗口。現(xiàn)在要給加油站地區(qū)設(shè)定一個(gè)大一點(diǎn)的補(bǔ)給量，那里飛機(jī)燃料每天都有泄露。單擊 Assign> Window把鼠標(biāo)移到加油站區(qū)的左頂角，單擊左鍵，然后在加油站區(qū)拉出一個(gè)框，再單擊左鍵。此時(shí)會(huì)彈出一個(gè) Recharge-Assign Window的窗口，如下圖所示：?jiǎn)螕鬘ew建立一個(gè)新的補(bǔ)給區(qū)。Zone匕欄的值會(huì)增加為2,而且Zone的顏色預(yù)覽框會(huì)變成藍(lán) 色。輸入以下數(shù)值：Stop Time day : 7300 Recharge mm/yr : 250單擊OK確定。

41、此時(shí)加油站地區(qū)的網(wǎng)格就會(huì)變成藍(lán)色的，這表明該區(qū)域?yàn)閆one#2的補(bǔ)給值。現(xiàn)在要返回到第一層的全視圖。單擊F6-Zoom out二、設(shè)定定水頭邊界條件下一步是對(duì)模型中潛水和承壓含水層的南北邊界設(shè)定定水頭邊界條件。從左邊的工具欄中的YES保存補(bǔ)給量數(shù)含水層設(shè)定的。要做從主菜單中選擇 Boundaries/Constant Head可選框中把 Recharge改選為 Const.Head。據(jù)。第一個(gè)定水頭邊界條件是沿著模型的北邊界對(duì)上層潛水 這一步，就要用到AssignLine 工具，沿著北邊界的一行網(wǎng)格設(shè)定定水頭邊界 在左邊的工具欄中點(diǎn)擊AssignLine把鼠標(biāo)指針指向西北角的網(wǎng)格(左頂角)，在

42、此單擊左鍵固定直線的起點(diǎn)。接著把鼠標(biāo)指針指向東 北角的網(wǎng)格(右頂角)，在此單擊右鍵，表示直線的終點(diǎn)。這時(shí)，一條水平網(wǎng)格線就顯示成粉紅色，同時(shí)出現(xiàn)一個(gè)Constant-Head-Assign Line 的窗口，如下圖所示：輸入以下數(shù)據(jù)：Description:CH-Upper Aquifer-NorthStop Time(day): 7300Start Time Head (m): 19Stop Time Head (m): 19網(wǎng)格中粉紅色的線現(xiàn)在就變成深紅色，表明對(duì)這些網(wǎng)格已經(jīng)確定了定水頭邊界條件 下一步要把第一層的定水頭值拷貝給第二層中相同位置的網(wǎng)格中。并把它拷貝給該在左邊的工具欄中單擊C

43、opyLayer就會(huì)彈出一個(gè)Copy Layer的窗口，如右圖所示 這個(gè)窗口允許從當(dāng)前層中選擇定水頭邊界條件網(wǎng)格, 模型的其它層。在Description 欄中選擇, CH-Upper Aquifer-North并在下方組合框中選中 L Layer 2 1 1 I單擊OK確定，就把第一層中選定的定水頭拷貝給第二層了。下一步是沿著模型的北邊界對(duì)下層承壓含水層設(shè)定定水頭邊界條件。在左邊的工具欄中選擇Goto,就會(huì)彈出Go To Layer的窗口。在Layer you wish to go(想要切換的層)一欄中，默認(rèn)值為 1,把它該為5。單擊OK就跳到第五層。查看導(dǎo)航立方體，確定現(xiàn)在處于第五層。當(dāng)前

44、層號(hào)在Layer (K)中顯示。現(xiàn)在沿著第五層北邊界的網(wǎng)格設(shè)定定水頭邊界條件。在左邊的工具欄中點(diǎn)擊AssignLine把鼠標(biāo)指針指向西北角的網(wǎng)格(左頂角)，在此單擊左鍵固定直線的起點(diǎn)。接著把鼠標(biāo)指針指向東 北角的網(wǎng)格(右頂角)，在此單擊右鍵，表示直線的終點(diǎn)。這時(shí)，一條水平網(wǎng)格線就顯示成粉紅色，同時(shí)出現(xiàn)一個(gè)Constant-Head-Assign Line 的窗口，并輸入以下數(shù)據(jù)：Description:CH-Lower Aquifer-NorthStop Time(day):7300Start Time Head (m):1818Stop Time Head (m):單擊OK確定網(wǎng)格中粉紅色的

45、線現(xiàn)在就變成深紅色，表明對(duì)這些網(wǎng)格已經(jīng)確定了定水頭邊界條件。下一步是沿著模型的南邊界對(duì)下層承壓含水層設(shè)定定水頭邊界條件。在左邊的工具欄中點(diǎn)擊AssignLine把鼠標(biāo)指針指向西南角的網(wǎng)格(左底角)，在此單擊左鍵固定直線的起點(diǎn)。接著把鼠標(biāo)指針指向東 南角的網(wǎng)格(右底角)，在此單擊右鍵，表示直線的終點(diǎn)。這時(shí)，一條水平網(wǎng)格線就顯示成粉紅色，同時(shí)出現(xiàn)一個(gè) Constant-Head-Assign Line的窗口，并輸入以下數(shù)據(jù)：Description:CH-Lower Aquifer-SouthStop Time(day):7300Start Time Head (m):16.5Stop Time H

46、ead (m):16.5單擊OK確定。網(wǎng)格中粉紅色的線現(xiàn)在就變成深紅色，表明對(duì)這些網(wǎng)格已經(jīng)確定了定水頭邊界條件。下一步要把第五層的定水頭值拷貝給第六層中相同的網(wǎng)格中。在左邊的工具欄中單擊CopyLayer就會(huì)彈出一個(gè)Copy Layer的窗口，這個(gè)窗口允許從當(dāng)前層中選擇定水頭邊界條件網(wǎng)格，并把它 拷貝給該模型的其它層面。選中回Copy all/selected groups , 并選中團(tuán)Layer 6單擊OK確定，把第五層中選定的定水頭拷貝給第六層。現(xiàn)在沿著某一列觀察模型的橫截面圖，看定水頭邊界條件給定的地方。在左邊的工具欄中點(diǎn)擊View Column在模型中不連續(xù)隔水層地帶選擇一列并點(diǎn)擊左鍵

47、(如第29列)。就會(huì)出現(xiàn)如下圖所示的一個(gè)橫截面圖：同時(shí)還可以顯示滲透性分區(qū)。在底部的工具欄中單擊F9-Overlay按扭，就會(huì)彈出一個(gè)如下所示的 Overlay Control 的窗 口，窗口中包含一系列可用的疊加層，可以選中或取消。點(diǎn)擊下拉箭頭 二I,選擇 Properties ,并在其層疊菜單中選中Pr Prop （F） -Conductivities單擊OK,顯示滲透系數(shù)的疊加層。模型的中間層（第三、四層）的網(wǎng)格應(yīng)該呈深藍(lán)色，除了在不連續(xù)弱透水層區(qū)域，這表示它們屬于Zone#2的滲透系數(shù)。現(xiàn)在要返回到模型的各層的平面圖：在左邊的工具欄中點(diǎn)擊View Layer,把鼠標(biāo)放在橫截面的第一層，

48、并單擊左鍵。三、給定河流邊界條件仔細(xì)觀察地圖，會(huì)發(fā)現(xiàn) Waterloo River 在模型區(qū)域內(nèi)沿著南部流動(dòng)（用淺藍(lán)色顯示）。下面的 步驟是用來(lái)描述怎樣對(duì)模型的頂層設(shè)定河流邊界條件。在主菜單中選擇Boundaries/River ,就進(jìn)入河流輸入界面了。在左邊的工具欄中點(diǎn)擊AssignLine從網(wǎng)格的西南面開(kāi)始，以地圖為引導(dǎo)，沿著河流點(diǎn)擊左鍵，使之成為數(shù)字化的網(wǎng)格。當(dāng)?shù)竭_(dá)東南面 的邊界時(shí)，在該線的終點(diǎn)右擊鼠標(biāo)。此時(shí)，沿著 Waterloo River的網(wǎng)格就會(huì)顯示粉紅色，同時(shí)會(huì) 彈出一個(gè)River -Assign Line 的窗口（如下圖所示）。它會(huì)提示一些有關(guān)這條河流的相關(guān)信 息。通常，河流

49、邊界條件還需要河床的水力傳導(dǎo)系數(shù)。但是每個(gè)網(wǎng)格的水力傳導(dǎo)系數(shù)有賴于河流穿過(guò)每 個(gè)網(wǎng)格的長(zhǎng)度和寬度。因此，在類(lèi)似這樣的模型中，隨著網(wǎng)格的大小不同，水力傳導(dǎo)系數(shù)也會(huì)因此 而改變。為了適應(yīng)這種情況，Visual MODFLOW提供在該線的起點(diǎn)和終點(diǎn)分別輸入河流實(shí)際的物理參數(shù)（長(zhǎng) 和寬），然后根據(jù)標(biāo)準(zhǔn)公式，為每個(gè)網(wǎng)格計(jì)算出合適的水力傳導(dǎo)系數(shù)。為河流的邊界條件輸入以下信息：Description: Waterloo River選擇v： Assign to appropriate layer選擇 Use default conductance formula選擇Linear gradient注意當(dāng)選中Li

50、near gradient時(shí)，會(huì)在該窗口中出現(xiàn)兩張表格Start Point（起點(diǎn)）和End Point （終點(diǎn)），如下圖所示：在 Start Point 表格中，輸入：Stop time (day):7300River Stage elevation (m): 16.0River Bottom elevation (m): 15.5Riverbed Thickness (m):0.1Riverbed Kz (m/sec):1e-4River Width (m):25在 End Point 表格中，輸入：Stop time (day):7300River Stage elevation (m)

51、: 15.5River Bottom elevation (m): 15.0Riverbed Thickness (m):0.1Riverbed Kz (m/sec):1e-4River Width (m):25單擊 OK 確定。當(dāng)河流的邊界條件確定后，河流的網(wǎng)格會(huì)變成深藍(lán)色，如下圖所示：第六部分：質(zhì)點(diǎn)示蹤質(zhì)點(diǎn)示蹤用來(lái)確定地下水質(zhì)點(diǎn)從一已知點(diǎn)出發(fā)后的最優(yōu)流動(dòng)路徑。在這個(gè)部分，將會(huì)指導(dǎo)你使用向前質(zhì)點(diǎn)示蹤法來(lái)確定最優(yōu)的污染物運(yùn)移通道所必須的步驟。在主菜單中選擇Particles ，會(huì)出現(xiàn)保存邊界資料的對(duì)話框，點(diǎn)擊YES。此時(shí)就會(huì)出現(xiàn)Particles 的輸入界面，在界面的左邊有各種輸入質(zhì)點(diǎn)信息的選

52、項(xiàng)。在這個(gè)模型中，要求在加油站附近指定一行質(zhì)點(diǎn)，用來(lái)確定這個(gè)地區(qū)的污染物進(jìn)入地下水后，將向哪里運(yùn)移。第一步是再一次放大加油站地區(qū)：?jiǎn)螕?F5-Zoom in 把鼠標(biāo)移到加油站的左上角，并單擊左鍵。然后在加油站區(qū)拉出一個(gè)框，再單擊左鍵關(guān)閉放大窗口。單擊 Add>Add Line 在第一層中增加一排質(zhì)點(diǎn)。把鼠標(biāo)移到加油站區(qū)左邊界中點(diǎn)處，并點(diǎn)擊左鍵；在加油站區(qū)右邊界中點(diǎn)處再點(diǎn)擊左鍵，就出現(xiàn)一條直線。此時(shí)會(huì)彈出一個(gè)Line Particle 的窗口，如下圖所示：默認(rèn)質(zhì)點(diǎn)數(shù)目是10。把這質(zhì)點(diǎn)數(shù)目該為5。單擊 OK 確定。橫穿加油站區(qū)域的一排綠色質(zhì)點(diǎn)表示向前示蹤質(zhì)點(diǎn)。現(xiàn)在要返回到模型的全屏顯示：在底

53、部的菜單欄中單擊F6-Zoom out ?，F(xiàn)在， 這個(gè)使用指南的第一模塊已經(jīng)完成了，而且已經(jīng)成功建立了一個(gè)三維的地下水流動(dòng)和質(zhì)點(diǎn)示蹤模型。下一個(gè)模塊用于建立溶質(zhì)運(yùn)移模型，并為溶質(zhì)運(yùn)移模型設(shè)定特征值和邊界條件。附：River模塊River (RIV) Boundary ConditionsVisual MODFLOW supports the River Package included with MODFLOW. TheVisual MODFLOW input data for River grid cells is stored in the projectname.VMB file (se

54、e Appendix A), while the MODFLOW input data for River grid cells is stored in the projectnameRIV file (see the MODFLOW-2000 Reference Manual .PDF file provided with the Visual MODFLOW installation CD-ROM, in the Manual folder).The River boundary condition is used to simulate the influence of a surface water body on the groundwater flow. Surface water bodies such as rivers, streams, lakes and swamps may either contribute water to the groundwater system, or act as gr