Modflow操作截圖說(shuō)明書(shū)

1、實(shí)用文檔、創(chuàng)建文件及模型離散圖模型主菜單主菜單屏幕包括以下選項(xiàng)：文件(F)輸入運(yùn)行(R)輸出(0)設(shè)置幫助(H)圖模型主菜單選項(xiàng)(1)文件：包括打開(kāi)、 新建、保存、關(guān)閉Visual MODFLOW ,導(dǎo)入、導(dǎo)出 MODFLOW 等。 V皿al MODFLOW1fan pyiput Scfep HcgOpwi. NewSave As 口 Import MQDFLCW Impnit MOOFiOlV P也 LKpfirL.Chang* UnitL_IPrefereficers.Pnpjert DerailsPrirftlEidt圖文件菜單(2)點(diǎn)擊Files-New,選擇文件的儲(chǔ)存位置，選擇或新建

2、文件夾，命名模型文件名，點(diǎn)擊確 定，新建一個(gè)模型后，出現(xiàn)下圖窗口。該窗口用于選擇水流和溶質(zhì)運(yùn)移模擬所需的數(shù)值求解 方法；水流和溶質(zhì)運(yùn)移模型參數(shù)的單位以及模型的起始日期和時(shí)間。圖選擇模型程序其中，對(duì)于水流模型，可以選用下面的數(shù)值求解方法：USGS MODFLOW-96 from SWS （飽和，密度恒定）；USGS MODFLOW-2000 from SWS （飽和，密度恒定）；USGS MODFLOW-2005 from SWS （飽和，密度恒定）；USGS SEAWAT from SWS （飽和，密度變化）；MODFLOW-AURFACT from HGL（變飽和度）。對(duì)于運(yùn)移模型， Vis

3、ual MODFLOW 支持以下圖中程序，其中 MT3DMS , MT3D99 , RT3D運(yùn)移引擎支持多組分污染物運(yùn)移，而 MT3Dv1.5僅支持單組分運(yùn)移。圖選擇模型程序(3) Next,水流設(shè)置在運(yùn)行類(lèi)型處，可以設(shè)置非穩(wěn)定流模型和穩(wěn)定流模型，在非穩(wěn)定流模型中， 要設(shè)置模型的起始日期和起始時(shí)間，輸入某些邊界條件數(shù)據(jù)時(shí)，就要定義時(shí)間數(shù)據(jù)。H "Mill! *4.1 耳4 1口也rghfn，上文案大全am 1 若升1色眥叫1Lj|覺(jué)聞5mow1wt世算早審占*KO0OC1o/i號(hào)軍* ,軋通電ooonWb擊K遇累16人1E-5七IE4lAn»TrtT*'v0

4、7;U15句030禧咕福充00農(nóng)ES舞耳理04V圖水流模型默認(rèn)參數(shù)設(shè)置(4) Next,運(yùn)移設(shè)置，在吸附下拉菜單中可選擇吸附模式，在反應(yīng)下拉菜單中選擇反應(yīng)類(lèi) 型。圖溶質(zhì)模型默認(rèn)參數(shù)設(shè)置(5)模型區(qū)：用來(lái)創(chuàng)建模型網(wǎng)格，上面窗口用來(lái)導(dǎo)入模擬區(qū)地區(qū)，確定模型區(qū)維數(shù)，以及 定義有限差分網(wǎng)格的行列數(shù)和分層數(shù)。圖創(chuàng)建模型網(wǎng)格完成后，出現(xiàn)選擇模型區(qū)界面，確定導(dǎo)入底圖中的模擬區(qū)，并通過(guò)輸入 X, Y對(duì)導(dǎo)入的 圖片進(jìn)行校準(zhǔn)，要使圖片與坐標(biāo)系匹配，至少知道兩個(gè)已知坐標(biāo)的地理參考點(diǎn)， 點(diǎn)擊地圖上 的參考點(diǎn)，出現(xiàn)地理參考點(diǎn)對(duì)話(huà)框提示輸入所選位置的 XY坐標(biāo)，將圖片附上坐標(biāo)信息， 完 成。圖模型分割范圍設(shè)置和坐標(biāo)賦值G

5、rid模塊-網(wǎng)格窗口* Visual MOD FLOW -File Grid Wells Pr寸 Visual F«»ODFI nWView ColumnView RowView LayerGotoPreviousNextGrid一Inactive Cells>Edit GridEdit ExtentsSmoothirrq卜Elevation 卜Contoiuri ngMExport ElevationImport Elevation圖 Grid模塊(1) View column , View row , View layer 選擇查看類(lèi)型。(2) Goto prev

6、ious , next為選擇層數(shù)，即模型的第一層，第二層等(3) Inactive cells：活動(dòng)單元的設(shè)置，自上向下分別為設(shè)置活動(dòng)區(qū)域，設(shè)置非活動(dòng)區(qū)域，設(shè)Ud；t Grig匕Ed Ft ExtentsSECiGthinq卜.卜Export ELvation Import EleatkinMaric Paly, ActiveMark Poly. InactiveMark SingleCep/ SingleCopy Polygon置活動(dòng)單元格，復(fù)制活動(dòng)單元格，復(fù)制活動(dòng)區(qū)域。圖定義有效、無(wú)效網(wǎng)格(4) Edit grid : Edit layer和Edit column ,通過(guò)該部分可以添加、刪

7、除、移動(dòng)和加密剖分網(wǎng) 格的設(shè)置。可根據(jù)需要對(duì)網(wǎng)格進(jìn)行加密、減疏等操作。圖編輯grid:二III MIIIIIMN II Mill III VOIIIIISBIMIIIIVIII !1nmninninnnnna RttllitllBBlIIBII MIIINBIISS 5111=>£=(5) Assign elevation,可以輸入某點(diǎn)，線(xiàn)，面，窗口的高程值。ja HiBiirfei'iiiniiBiii aai .hi illi mi ii'i-iiii mi 出工Gi/in上I曰，CcrnBuRr1. ， E 一口口 rt-Import Elevaticr圖

8、高程值輸入(6) Coutouring-高程等值線(xiàn)的繪制，可以繪制頂板等值線(xiàn)，厚度等值線(xiàn)和底板等值線(xiàn)，options 是分別對(duì)上述三種等值線(xiàn)進(jìn)行設(shè)置。耳ion EIey占titan”I注:否曰砧正茂璉口E-p-rt E evMicnIrnu 0rt Elwatign圖繪制高程等值線(xiàn)(7) export elevation :導(dǎo)出高程值。(8) Import elevation-導(dǎo)入高程，出現(xiàn)創(chuàng)建網(wǎng)格高程窗口。層面-選擇輸入高程值所屬層面，(地表高程，第一層底板高程，第二層底板高程等)圖高程輸入選項(xiàng)-選擇導(dǎo)入數(shù)據(jù)，將事先處理好的高程text文件導(dǎo)入，并在此選擇插值方法(自然鄰點(diǎn)法，克力格，或反距

9、離方法），默認(rèn)情況下，插值方法選擇自然鄰點(diǎn)法。_lf«1 ata'h EdCUr-flWX- WH'-Mig34 an餐.HraWa 1 比，中口：ig H 制fvwfUMALMlhdniuX V圖高程插值然后彈出下面的窗口，所需的數(shù)據(jù)由數(shù)據(jù)文件的列組成，并且與列號(hào)相匹配， 而這個(gè)數(shù)據(jù)框用來(lái)與每一個(gè)數(shù)據(jù)區(qū)域相對(duì)應(yīng)，每一個(gè)數(shù)據(jù)域里鍵入適當(dāng)?shù)臄?shù)字以使得列#1和X坐標(biāo)軸相對(duì)應(yīng)，使列#2和Y坐標(biāo)軸相對(duì)應(yīng)，以及列 #3與高程相對(duì)應(yīng)，當(dāng)所有列都對(duì)應(yīng)后，next按鈕被激活。然后彈出有效數(shù)據(jù)檢測(cè)框，保證數(shù)據(jù)點(diǎn)位于下面坐標(biāo)系和單位框所確定的模型區(qū)域內(nèi)。Requi C»atsM

10、hoh fc" .-.：4i4rp! fi<.fte-r1.2dNnEQ zai an- ：- -牡 3H'poo, zlii-ie I - -UR Jli-J"C口 21 J-e E7B.H1Bllg3 = A- 4 M- R B AH zQ BAIlzad E.匚Mif r a ；7 ie TtF“nni r 2JIES' 1 8 2 .J T 己f口修日士 4 m 71FrinFl 1Cdl.+T1!rw # Z 411 IIJ|?,_|I aiU4Z30.1i 4i(n*rr網(wǎng)- 勇fUJ日k.三qlORPSF17y 4-IOb&t9.

11、：=bh 4-1111j44 1 £E7 an WTrt.113310-in17舊 541O5E7tET 0-IWS.i.i /Vht 5*10714e.?11t eaim 3.3r日J(rèn)IT B5UI 7.53 31r i w #f r ” 141 口4 E £ Hl , 7 J I FS J ! ,： - I MHd口Car-eelMMTtr Q 1HIU-Krj ii B物口再向-7 fllOBTEil 尸尹 tini., m- -. 41dtS4r 白金 物皿一不Help:二：:T 二I 11 .二112' w IKMstaBtHst。X-日”占1電片211

12、11Wn& c- m it i-' M耐用hJ14 IS1臼 105-31 3l91 5 1?三m S 443充ZB/t、H山木圖17點(diǎn)擊Finish進(jìn)入到坐標(biāo)系統(tǒng)和單位窗口，坐標(biāo)系統(tǒng)選擇模型坐標(biāo)，高程單位選擇米， 完成高程導(dǎo)入。圖坐標(biāo)系統(tǒng)和單位設(shè)置三、 Well模塊-井模塊：該模塊包括抽水井，水頭觀測(cè)井和濃度觀測(cè)井三個(gè)部分?；?Visual MODFLOW - C：J?ersAdminitratorDeskt©.'/mf Hie Gridyiuw CcluView RowView LaveWell Properties BoundariesPumping W

13、ellsHead Observatio-r lnVellsCorrc. Observation '»VellsAntn圖井模塊(1)抽水井刪除井、編輯井、移動(dòng)井和復(fù)左側(cè)工具欄按鈕是各種針對(duì)抽水井的選項(xiàng)，包括添加井、 制井。圖抽水井設(shè)置添加井需將光標(biāo)移動(dòng)到目標(biāo)位置，單擊鼠標(biāo)左鍵在該處添加一口抽水井，彈出一個(gè)新井窗口，提示輸入井的詳細(xì)信息。包括井的名字和位置坐標(biāo)，ok。然后，在右上處輸入濾管頂部高程(top, m),底部高程(bottom , m);輸入抽水時(shí)段，開(kāi)始時(shí)間，結(jié)束時(shí)間和流量。New v/ell圖抽水井參數(shù)設(shè)置(2)水頭觀測(cè)井和濃度觀測(cè)井Wells-head/conc

14、entration observation wells-import Obs.,導(dǎo)入數(shù)據(jù)，按彈出對(duì)話(huà)框設(shè)置對(duì) 應(yīng)數(shù)據(jù)列，完成觀測(cè)水位數(shù)據(jù)的導(dǎo)入。® Visual MODFLOV? - C：WsersAdView ColumnView RowView LayerGotoPreviousNextHead Obs. Add Obs.Delete Ob5,卜Edit QbsiMove Obs,Edit GroupsDatabaseExportImpoift* Visual MODFLOW C:User5AdFile Grid Weils Properties 6cFile Grid Well

15、sConcDOIAdd QbsOulet計(jì) Qb后lEditMa告 Ob&口亙 tabaExpertIe portProperties &q圖觀測(cè)井在Observation points添加井號(hào)和觀測(cè)井的觀測(cè)層位，并在下面 Observations輸入時(shí)間和觀測(cè)值。圖觀測(cè)井設(shè)置四、Properties模塊-地下水流屬性的賦值.川 won (TA . C.Uaf AdEin 由川馬EWp說(shuō)斷* Hrile Grid WeillPropertrFa Boundarirs P(artklu> ZB*ud Touli HupView ColumnMiisw RowView Lay

16、erCcndijctrvityStoragenitiil HeadcVadose Zone«.BCone. 0gCuncDQI,Idd Obi.Delete Qb工*Edh Obs,Move Obs4Bulk D白尸號(hào)的Species ParametersEdit Cc-nstflunt Pau mete3- nitial LcncentratbnDispersionEdit Gfc-uqeExpertImport1圖含水層參數(shù)設(shè)置水流模型中分別包括導(dǎo)水系數(shù)、儲(chǔ)水系數(shù)和初始水頭，每個(gè)系數(shù)設(shè)置分別包括對(duì)某區(qū)域 的賦值，復(fù)制，數(shù)據(jù)庫(kù)，繪制等值線(xiàn)和導(dǎo)出數(shù)據(jù)部分。Conductivity八二

17、4口 nSpy，ContgunnqExp 口 flgta mpeirtViev; ColumnView RdwView LayerGotoPreviousNext11Storage-AssignCopyDatabaseContourinq)Export DataImportVie聊 ColumnView RowView LayrGotoPrzidUHNuxtInitial Heads Assign_lJCopyDatabaseConteurinaExport DataImport|通常在建模時(shí)，對(duì)模擬區(qū)進(jìn)行參數(shù)分區(qū)， 分別對(duì)每個(gè)分區(qū)的參數(shù)進(jìn)行輸入。導(dǎo)水系數(shù)部分分別輸入x方向、y方向和z方向的導(dǎo)

18、水系數(shù)。儲(chǔ)水系數(shù)部分分別輸入儲(chǔ)水系數(shù)、給水度、有效孔隙度和總孔隙度。初始水頭通常采用導(dǎo)入的方法，將事先準(zhǔn)備的text導(dǎo)入，同樣通過(guò)輸入列號(hào)，將數(shù)據(jù)與各列相對(duì)應(yīng)，然后進(jìn)入下一步，點(diǎn)擊插值( interpolate),完成。通過(guò)重復(fù)該步驟分別對(duì)每一 層進(jìn)行初始水頭的賦值。在包括溶質(zhì)運(yùn)移的模型中，對(duì)模擬區(qū)初始濃度進(jìn)行賦值，同樣采用導(dǎo)入的方法，操作步 驟與初始水頭的導(dǎo)入步驟相同。五、Boundary模塊-邊界條件，Shu a MC-DfLOV/' - C UseriAcm niitrator.De&lirtE"Hie GridProperties Sourdaries 附巾d

19、 ZEud qqIs HelpCoiuld it Curtientrdtior'JrTipurt IVJKE 11 Rivei NcIwjpRecharge C&ncentralion Evo pots Concentration oirt £ourc4C<jni>ldrit Head (CUD)River (FttVlStream (STR)Seriral Head (.GHB)Drain (DRNLake (LAK)Wall (HFBRecharge (RCH)Evapotran epi ratio r圖邊界條件設(shè)置(1) Constant head

20、(CHD )-定水頭邊界Boundary-constant head-assign-line,分別輸入起始時(shí)間，結(jié)束時(shí)間，起始時(shí)間水頭和結(jié) 束時(shí)間水頭，點(diǎn)擊ok,接受這些數(shù)值。網(wǎng)格單元上出現(xiàn)黑紅色，表明一個(gè)定水頭邊界條件已經(jīng)賦給了該單元。rf rJin網(wǎng)喇，口直伸的«IKM)Ft正 |CH Umi&加Ttolhr福呼艇n居4中可舊：出：互出耳同，.耳:疑暖網(wǎng)巾寸力加 曹/1II圖定水頭邊界設(shè)置(2) River-河流邊界選擇設(shè)置到適當(dāng)層，采用默認(rèn)導(dǎo)水系數(shù)公式，分別輸入結(jié)束時(shí)間，河水位，河床底高，河床厚度，河床滲透系數(shù)和河流寬度。若選中線(xiàn)性梯度， 則分別輸入起始點(diǎn)和結(jié)束點(diǎn)的上述

21、信息。圖River邊界設(shè)置(3) General head boundary-通用水頭邊界分別輸入結(jié)束時(shí)間，邊界水頭，邊界距離，平均滲透系數(shù)K和界面。圖通用水頭邊界設(shè)置(4) Drain-排水溝分別輸入結(jié)束時(shí)間，排泄高程，單位面積conductance。圖排水溝設(shè)置(5) Lake-湖分別輸入結(jié)束時(shí)間，湖水位，底板高程，入滲，蒸發(fā)，徑流，抽水速率，湖底泥厚度和 透水系數(shù)。圖湖模塊設(shè)置(6) Wall分別輸入反應(yīng)墻起始點(diǎn)和終止點(diǎn)的厚度和導(dǎo)水系數(shù)。圖Wall模塊設(shè)置(7 ) Recharge-#給輸入停止時(shí)間和補(bǔ)給速率。圖Recharge模塊設(shè)置8 8 ) Evapotranspiration-蒸

22、發(fā)輸入蒸發(fā)的停止時(shí)間，蒸發(fā)速率和深度。圖蒸發(fā)設(shè)置(9 ) Constant concentration-定濃度邊界將污染源的位置和濃度進(jìn)行定義。選中污染源位置后分別賦值結(jié)束時(shí)間，污染源濃度。Assign Ccnsta nt Concentration補(bǔ)紿法度或值圖模型運(yùn)行圖定濃度設(shè)置(10) Recharge concentration 和 Evapotranspiration concentration分別設(shè)置污染物的補(bǔ)給濃度和蒸發(fā)濃度，輸入起始時(shí)刻、結(jié)束時(shí)刻、污染物濃度。 六、Run-運(yùn)行模塊按F10回到主菜單，點(diǎn)擊 Run進(jìn)入到運(yùn)行模塊。(一)Modflow運(yùn)行模塊圖設(shè)置運(yùn)行參數(shù)(1)

23、time steps,分別輸入其實(shí)時(shí)間，結(jié)束時(shí)間，時(shí)間步長(zhǎng)。圖時(shí)間步長(zhǎng)設(shè)置(2) initial head圖初始水位賦值選項(xiàng)(3)求解方法設(shè)置(默認(rèn)值)圖求解方法設(shè)置4 4） Recharge （默認(rèn)）圖補(bǔ)給設(shè)置5 5） Layer對(duì)各層承壓性進(jìn)行設(shè)置，layer type。To assign miultiple kyers select appropriate rows Iby <Ctrl-cm <Shilt> click on row# indicator barLei LAYCONInterblock transmissiviyLayer type103aOttHarm

24、cinic meana3:匚 Qnfigd/U mcpnf ined 丫即值皿 ILayer Settings口: Confined, cunsUnt XT|V Unconfirmed2 匚orfiriBd/U nc口口filed,忡iable 5, cn涵 mt I3:Corfined/Unconfined, variable S, TOKC-ancel圖層承壓性設(shè)置（二）MT3D 模塊喻 Visual MODFLDW - f C:UseraA.dministrator'.Deskto2-vml : VarOOlAle MODFLOW 2000 MODPATH MF.口M£

25、 Optimization Run r&ols H31PView tcJumnV適叫RemVieiv L3爐后rGotoPmvi 0 U 與NextK5 Me圖溶質(zhì)運(yùn)移模塊運(yùn)行參數(shù)設(shè)置(1)解法選項(xiàng)窗口Slukun Me曲odAd1七匚EnP.5-5 ch Opifcris；M Lhaazt日 比 IMOCJH-4dsti 于id -畫(huà) £ Wwlg和il* iJsq &CG 比格 |5 gl殖£幅1訊Iirr«6d 時(shí)Glrii.l*rp 內(nèi) 丁6 0d.RM*wnwn 5itp 51M C修，Mugi 耳 1日g&iy ocilimt平

26、U bs函&土£ 口即篤iyKM T.眄5z5 D mrracliiorior HI Mickne1coM整 w her 口.法Lll<Caincfl圖溶質(zhì)運(yùn)移求解方法設(shè)置(2) Output time-輸出時(shí)間MT3D輸出控制的默認(rèn)輸出設(shè)置是只在模擬結(jié)束時(shí)的濃度計(jì)算結(jié)果。在模型輸出結(jié)果中需要指定查看模擬結(jié)果的時(shí)間。輸入模擬時(shí)間，最大運(yùn)移步長(zhǎng)，選擇指定輸出時(shí)間。圖模型運(yùn)行圖輸出時(shí)間設(shè)置(3) Run,選擇需要運(yùn)行的程序，點(diǎn)擊Translate&Run ,開(kāi)始運(yùn)行。(三)output-輸出可視化(1)顯示等值線(xiàn)圖進(jìn)入結(jié)果輸出模塊，點(diǎn)擊Maps-conturing maps,可以選擇顯示不同類(lèi)型的等值線(xiàn)圖，分別繪制水頭等值線(xiàn)，降深等值線(xiàn)，濃度等值線(xiàn)，水位高程等值線(xiàn)，水位深度等值線(xiàn)等。點(diǎn)擊Time,選擇等值線(xiàn)的輸出時(shí)間，并通過(guò) option選項(xiàng)，對(duì)不同等值線(xiàn)圖進(jìn)行顯示設(shè)置。哂331 MOlJhLUW 匚愣-W，加mini磯31©1】。*或1口卜1