PMWIN操作說明

1、 PMWIN操作說明一、創(chuàng)建模型網(wǎng)格1、打開Grid菜單，選擇MeshSize。ModelDimension對(duì)話框彈出。2、輸入含水層的層數(shù)2,行數(shù)20、列數(shù)20,模型總長(zhǎng)度和寬度定為2000m。3、點(diǎn)擊OK。PMWIN的界面發(fā)生變化，此時(shí)顯示了模型的整體網(wǎng)格。4、從File菜單中選擇LeaveEditor或點(diǎn)擊LeaveEditor按鈕退出。D注意：創(chuàng)建文件夾和文件名必須是英文字符和數(shù)字，不能用中文。二、定義含水層類型1、打開Grid菜單，選擇LayerType。LayerOptions對(duì)話框彈出。2、點(diǎn)擊標(biāo)簽為Type列中的某一單元格，單元格中將出現(xiàn)一個(gè)帶下箭頭的按鈕,點(diǎn)擊該按鈕，出現(xiàn)一個(gè)

2、可供選擇含水層類型的列表。3、為第一層選取1：Unconfined，第二層選取2：Confined/Unconfined，然后點(diǎn)擊OK,關(guān)閉對(duì)話框。三、定義邊界條件在模型中，邊界上的各個(gè)計(jì)算單元均用不同的指定代碼表示不同的邊界類型：（1）“1”表示變水頭邊界，（2）“1”表示定水頭邊界，（3）“0”表示無流量邊界。1、從Grid菜單中選擇BoundryConditionIBOUND（Modflow）。2、擊鼠標(biāo)右鍵，PMWIN顯示一個(gè)CellValue對(duì)話框。3、在對(duì)話框中輸入一1,點(diǎn)擊OK。計(jì)算單元1,1,1被定義為定水頭計(jì)算單元。4、點(diǎn)擊Duplication按鈕，打開復(fù)制開關(guān)。如果Dup

3、lication鍵下陷，表示復(fù)制開關(guān)處于開啟狀態(tài)。移動(dòng)網(wǎng)格指針，可將當(dāng)前計(jì)算單元的數(shù)值復(fù)制到單元格指針經(jīng)過的所有計(jì)算單元，再次點(diǎn)擊Duplication鍵，關(guān)閉復(fù)制開關(guān)。5、點(diǎn)擊LayerCopy按鈕，打開層復(fù)制開關(guān)。該按鈕下陷，表示層復(fù)制開關(guān)處于開啟狀態(tài)。此時(shí)，從一層移至另一層，當(dāng)前層的所有值將被復(fù)制到另一層中。再次點(diǎn)擊該按鈕，將關(guān)閉層復(fù)制開關(guān)。6、點(diǎn)擊LeaveEditor按鈕或從File菜單中選擇LeaveEditor退出定義邊界條件界面。四、模型幾何參數(shù)賦值（一）模型層的頂板高度賦值1、打開Grid菜單，選擇TopofLayer（Top）。PMWIN顯示模型網(wǎng)格。2、打開Value菜單

4、選擇ResetMatrix或按Ctrl+R鍵。ResetMatrix對(duì)話框出現(xiàn)。3、在對(duì)話框中輸入58.85，點(diǎn)擊OK。第一層的頂板高度定義為58.85。4、按PgDn鍵移至第二層。5、重復(fù)步驟3，輸入第二層頂板高度48.85。6、從File菜單中選擇LeaveEditor或點(diǎn)擊LeaveEditor按鈕退出。（二）模型層的底板高度賦值1、從Grid菜單中選擇BottomofLayer（BOT）。2、重復(fù)上述步驟，輸入第一層、第二層的底板標(biāo)高值48.85、41.85。3、點(diǎn)擊LeaveEditor按鈕退出。五、時(shí)間及空間參數(shù)賦值在本實(shí)例中，時(shí)間參數(shù)包括時(shí)間單位、應(yīng)力期數(shù)、時(shí)段、水質(zhì)運(yùn)移時(shí)段。空

5、間參數(shù)包括初始水頭、水平及垂直方向的水力傳導(dǎo)系數(shù)。（一）時(shí)間參數(shù)賦值1、打開Parameters菜單，選擇Time.。TimeParameters對(duì)話框出現(xiàn)。在MODFLOW中,模擬時(shí)間被劃分為幾個(gè)應(yīng)力期，也就是說，在外部應(yīng)力作用下，時(shí)間間隔為常數(shù)，應(yīng)力期又劃分為時(shí)段。2、輸入第一應(yīng)力期的長(zhǎng)度365d。3、點(diǎn)擊OK接受其它默認(rèn)值。（二）初始水頭賦值1、打開Parameters菜單，選擇InitialHydrolicHead。PMWIN顯示模型網(wǎng)格。2、從Value菜單中選擇ResetMatrix（或按Ctrl+R鍵），在對(duì)話框中輸入52,點(diǎn)擊OK。3、將網(wǎng)格指針移至左上角1,1,1。4、點(diǎn)鼠標(biāo)

6、右鍵，在對(duì)話框CellValue中輸入53。5、從左上角向下拖動(dòng)鼠標(biāo)，當(dāng)前計(jì)算單元中的值53將被復(fù)制到鼠標(biāo)所經(jīng)過的所有左側(cè)計(jì)算單元中。6、到右上角，輸入52，向下拖動(dòng)鼠標(biāo)，當(dāng)前計(jì)算單元中的值52將被復(fù)制到鼠標(biāo)所經(jīng)過的右側(cè)邊界所有計(jì)算單元中。7、點(diǎn)擊LayerCopy按鈕，打開層復(fù)制開關(guān)。8、按PgDn鍵，移至第二層。9、點(diǎn)擊LeaveEditor按鈕，退出。（三）導(dǎo)水系數(shù)賦值水平方向1、從Parameters菜單中選擇HorizontalHydraulicConductivity。2、打開Value菜單，選擇ResetMatrix（或按Ctrl+R）,在對(duì)話框中輸入40,點(diǎn)擊OK。3、按PgD

7、n鍵，移至第二層。4、在第二層利用多邊形輸入方法定義兩個(gè)區(qū)，分別輸入40，25。5、點(diǎn)LeaveEditor按鈕退出。（四）導(dǎo)水系數(shù)賦值垂直方向1、從Parameters菜單中選擇VerticalHydraulicConductivity。2、打開Value菜單，選擇ResetMatrix，在對(duì)話框中輸入1,點(diǎn)擊OK。3、復(fù)制到第二層。5、點(diǎn)LeaveEditor按鈕退出。（五）面狀補(bǔ)給量的賦值1、打開Models菜單，選擇MODFLOWfRecharge。2、從Value菜單中選擇ResetMatrix.,在對(duì)話框的RechargeFluxL/T編輯欄中輸入0.00043。3、退出。（六）抽

8、水井定義與賦值打開Models菜單，選擇MODFLOWWell。將網(wǎng)格指針移到需要設(shè)計(jì)抽水井的單元。點(diǎn)鼠標(biāo)右鍵，輸入抽水量（如：-1500m3/d）,點(diǎn)擊OK。負(fù)值表示抽水井。退出。六、執(zhí)行水流模擬1、打開Models菜單，選擇MODFLOWRun.。RunModflow對(duì)話框彈出。2、點(diǎn)擊OK,開始水流模擬。在運(yùn)行MODFLOW之前，PMWIN將利用用戶自定義的數(shù)據(jù)為MODFLOW（及可選MODPATH）產(chǎn)生一個(gè)輸入文件，該文件列在RunModflow對(duì)話框的列表中。如果Generate標(biāo)記為，表示只產(chǎn)生了一個(gè)輸入文件。可以點(diǎn)擊該按鈕鎖定Generate標(biāo)記。一般情況下，不需更改該標(biāo)記。七、

9、模型結(jié)果檢驗(yàn)在水流模型模擬中,MODFLOW在列表文件pathOUTPUT.DAT中有一個(gè)詳細(xì)運(yùn)行記錄。如果成功地完成了一個(gè)水流模擬，MODFLOW將以多種非格式化文件（二進(jìn)制）保存模擬結(jié)果。如下表。在運(yùn)行MODFLOW之前，用戶可以控制這些非格式化（二進(jìn)制）文件的輸出。具體操作：從Models菜單中選擇ModflowOutputControl。若InterbedStoragePackage被激活，則產(chǎn)生pathINTERBED.DAT輸出文件。為了檢驗(yàn)?zāi)M結(jié)果的質(zhì)量，MODFLOW在每個(gè)時(shí)段結(jié)尾為整個(gè)模型作了水量均衡計(jì)算，并保存在OUTPUT.DAT文件中。一個(gè)水均衡提供了一個(gè)數(shù)學(xué)解的所有可

10、解性跡象。MODFLOW的輸出文件文件內(nèi)容PathOUTPUT.DAT詳細(xì)運(yùn)行記錄和模擬報(bào)告PathHEADS.DAT水頭PathDDOWN.DAT降深、初始水頭與計(jì)算水頭之差PathBUGET.DAT各計(jì)算單元的流量項(xiàng)PathINTERBET.DAT每一層中整個(gè)含水層的沉降、壓密及壓密之前的水位PathMT3D.FLO與MT3D/MT3DMS的界面文件，這個(gè)文件由MT3D/MT3DMS提供的軟件包LKMT創(chuàng)建。Path為保存模型數(shù)據(jù)的文件夾。八、子區(qū)水均衡計(jì)算為了計(jì)算方便，每個(gè)計(jì)算單元的流量項(xiàng)被保存在文件pathBUDGET.DAT中，這些單一計(jì)算單元的流量項(xiàng)實(shí)際上是一個(gè)單元流向一單元的流量

11、，有四種類型：（1）一個(gè)單元向一個(gè)單元的應(yīng)力流量；或受外部應(yīng)力影響，流入或流出一個(gè)單元的流量，如抽水井或面狀補(bǔ)給；（2）計(jì)算單元間的儲(chǔ)存項(xiàng)：在各計(jì)算單元中，儲(chǔ)存量的增加和減少的速率。（3）計(jì)算單元間的定水頭流量項(xiàng)：流向或來自各個(gè)定水頭計(jì)算單元的凈流量。（4）計(jì)算單元間內(nèi)部的流量項(xiàng)：通過各計(jì)算單元截面的流量，也即相鄰計(jì)算單元間的流量。WaterBudgetCalculator利用逐個(gè)計(jì)算單元的流量來計(jì)算整個(gè)模型的水量均衡、子區(qū)及子區(qū)間的水量均衡。子區(qū)水均衡計(jì)算1打開Tool菜單，選擇WaterBudget。WaterBudget對(duì)話框彈出。在穩(wěn)定流計(jì)算中，不需改變Time的設(shè)置。2點(diǎn)擊Zones

12、。一個(gè)Zone代表一個(gè)水均衡計(jì)算子區(qū)，由050間的一個(gè)代號(hào)（區(qū)號(hào)）表示，每個(gè)區(qū)的區(qū)號(hào)必須賦予每一個(gè)計(jì)算單元。區(qū)號(hào)0表示該計(jì)算單元與其它區(qū)無關(guān)。Value菜單打開ResetMatrix.對(duì)話框，輸入1,點(diǎn)擊0K。PgDn鍵，移到第二層。從Value菜單中，選擇ResetMatrix.，輸入2,點(diǎn)擊0K。點(diǎn)擊退出按鈕，退出。點(diǎn)擊0K,退出WaterBudget對(duì)話框。PMWIN在文件pathWATERBDG.DAT中計(jì)算并保存流量項(xiàng)。流量單位為L(zhǎng)3/T。對(duì)每一時(shí)段、每一層的每一個(gè)區(qū)分別計(jì)算。流入?yún)^(qū)內(nèi)的流量記為IN,子區(qū)之間的流量記在FlowMatrix中，通過子區(qū)邊界的水平方向流量由HORIZEX

13、CHANGE給出。EXCHANGE（UPPER）給出了來自（IN）或流向（OUT）上部相鄰層的流量，EXCHANGE（LOWER）給出了來自（IN）或流向（OUT）下部相鄰層的流量，如LAYER=1,ZONE=1的EXCHANGE（LOWER）（表示第一層、第一區(qū)與下部相鄰層的流量交替項(xiàng)）。為了計(jì)算流向抽水井的確切流量，重復(fù)上述步驟計(jì)算水量均衡。將抽水井所在單元定義為一個(gè)子區(qū)，運(yùn)行計(jì)算，打開pathWATERBDG.DAT就可以看水量均衡情況。九、結(jié)果輸出除水量均衡計(jì)算之外，PMWIN提供了多種模擬結(jié)果檢驗(yàn)、創(chuàng)建圖形輸出的可能。流線及速度矢量可以通過PMPATH顯示。利用ResultExtra

14、ctor，可將任一時(shí)段任一層的模擬結(jié)果從非格式化（二進(jìn)制）結(jié)果文件中讀取，或保存在ASCIIMatrix文件中。一個(gè)ASCIIMatrix文件含有一個(gè)模型層的計(jì)算值（每計(jì)算單元一個(gè)值）。PMWIN能將ASCIIMatrix文件裝入模型網(wǎng)格中。以下給出操作步驟：用ResultsExtractor讀取或保存計(jì)算的水頭。2于第一層的計(jì)算水頭生成等水位線圖。3.使用PMPATH計(jì)算水流運(yùn)行軌跡及抽水井的截獲帶。（一）讀取及保存計(jì)算的水頭1、從Tool菜單中選取ResultsExtractor.。ResultsExtractor對(duì)話框彈出。對(duì)話框中含有四個(gè)選項(xiàng)標(biāo)簽，分別為MODFLOW、MOC3D、MT

15、3D、MT3DMS。在MODFLOW工作表中，可從ResultType下拉列表中選擇結(jié)果類型，也可以指定要讀取數(shù)據(jù)結(jié)果的層、應(yīng)力期、時(shí)間步長(zhǎng)。該電子表格顯示了一系列的行和列，一行與一列的交匯處為一個(gè)單元格。其中的每一單元格與模型層中計(jì)算單元一一對(duì)應(yīng)。下列步驟2、6以三個(gè)ASCIIMatrix文件為每一層保存了水頭值。2、點(diǎn)擊ResultType的下拉箭頭按鈕，選擇HydraulicHead。3、在Layer編輯框中輸入1。該實(shí)例中（穩(wěn)定流模擬只有一個(gè)應(yīng)力期一個(gè)時(shí)段）應(yīng)力期數(shù)及時(shí)段數(shù)均為1。4、點(diǎn)擊Read按鈕。第一層、第一應(yīng)力期，時(shí)段1的水頭值將被讀入電子表中。點(diǎn)擊滾動(dòng)條，可游覽整體結(jié)果。5、

16、點(diǎn)擊Save按鈕。SaveMatrix對(duì)話框出現(xiàn)。在SaveasType選項(xiàng)中，可根據(jù)用戶需要將數(shù)據(jù)結(jié)果保存為ASCII格式或SUFFER格式。給定文件名H1.DAT,保存文件，點(diǎn)擊OK。6、重復(fù)步驟3、4、5，分別保存第二層的水頭計(jì)算結(jié)果H2.DAT。7、點(diǎn)擊Close，關(guān)閉對(duì)話框。（二）生成計(jì)算水位等值線圖1、從Tools菜單中選取Presentation。在Presentation中定義的數(shù)據(jù)并不能被PMWIN使用，但可以用Presentation來臨時(shí)保存數(shù)據(jù)或顯示模擬結(jié)果的圖形。2、從Value菜單中選擇Matrix.（Ctrl+B）。BrowsMatrix對(duì)話框彈出。點(diǎn)擊Load按

17、鈕，可將ASCIIMtrix文件裝入表中，或點(diǎn)擊Save按鈕，將表中數(shù)據(jù)保存為一個(gè)ASCIIMatrix文件。還可通過菜單Value中的ResultExtractor讀取水位數(shù)據(jù)，或使用Apply按鈕在Presentation數(shù)據(jù)表中添加數(shù)據(jù)。3、點(diǎn)擊Load按鈕。LoadMatrix對(duì)話框出現(xiàn)。4、點(diǎn)擊Open按鈕，選取由ResultExtractor保存下的文件H1.DAT，點(diǎn)擊OK，H1.DAT被裝入。5、點(diǎn)擊OK,關(guān)閉BrowseMatrix對(duì)話框。6、打開Option菜單，選擇Environment（Ctrl+E）。EnvironmentOption對(duì)話框出現(xiàn)。對(duì)話框中包括三個(gè)選項(xiàng)列表。Appearance和CoordinateSystem用來修改模型網(wǎng)格的外觀及位置，Contours用來生成等值線圖。7、點(diǎn)擊Contour標(biāo)簽，選中Visible復(fù)稱選框，點(diǎn)擊RestoreDefault按鈕。點(diǎn)RestoreDefault按鈕，PMWIN默認(rèn)等值線數(shù)為11，并使用當(dāng)前層中的最大和最小值定義等值線的最大、最小值。如果FillContours被選取，點(diǎn)擊列標(biāo)簽Fill，可自選顏色填充等值線圖。用LabelFormat，可定義合適的數(shù)據(jù)格式。注意：當(dāng)退出編輯時(shí)，PMWIN將清除Visible復(fù)選框。8、點(diǎn)擊OK,退出EnvironmentOpt