visual_modflow實(shí)例

1、YourFirstGroundwaterModelwithPMWIN實(shí)例實(shí)例：一個(gè)由兩種不同地層構(gòu)成的含水系統(tǒng)，北部和南部為無流量邊界，東部和西部為河流，與地下水有充分的水力聯(lián)系，可視為定水頭邊界，西邊界與東邊界的水頭分別為9m，8m。Fig.81Configurationofthesampleproblem含水系統(tǒng)為非承壓各向同性，第一層、第二層的水力傳導(dǎo)系數(shù)分別為0.0001ms和0.0005ms，兩層的縱向傳導(dǎo)系數(shù)假設(shè)為橫向的10

2、%，有效孔隙率為25%，地面標(biāo)高為10m。第一層、第二層的厚度分別為4m、6m。含水層接受面狀補(bǔ)給量為810－9mm。在第一層中靠近西部邊界有一污染源?，F(xiàn)在我們的任務(wù)是，如何用一個(gè)位于東部邊界附近的完整井來隔離該污染源。那么，在該區(qū)域內(nèi)，需要建立一個(gè)數(shù)學(xué)模型來逐步計(jì)算抽水井所需的抽水量。只有當(dāng)抽水量足夠大時(shí)，被污染的地下水才能被抽水井全部截獲。我們將使用PMWIN來建立數(shù)學(xué)模型，使用PMPATH來計(jì)算抽水井的截獲帶。基于所計(jì)算的地下水流

3、場(chǎng)，使用MT3D和MOC3D來模擬污染源的運(yùn)移情況。并演示如何使用PEST和UCODE來校正水流模型。最后創(chuàng)建一個(gè)動(dòng)態(tài)演示方式，演示地下水被污染的發(fā)展變化過程。為了示范運(yùn)移模型，我們假想污染物溶解于地下水中的速度為1104gsm2，含水層縱向彌散和橫向彌散系數(shù)分別為10m和1m，延遲系數(shù)為2，初始濃度、分子擴(kuò)散系數(shù)及衰減速率均為0。我們將計(jì)算三年之后污染濃度分布情況，顯示在兩個(gè)含水層中點(diǎn)[xy]=[290310][390310]的濃度－

4、時(shí)間突變曲線。穩(wěn)定流模擬Fig.82Spatialdisretizationofanaquifersystemthecellindices一、一、創(chuàng)建模型網(wǎng)格創(chuàng)建模型網(wǎng)格1、打開Grid菜單，選擇MeshSize。ModelDimension對(duì)話框彈出。2、輸入含水層的層數(shù)3，行數(shù)30、列數(shù)30，行高、列寬20。第一含水層和第二含水層分別由一個(gè)模型層和兩個(gè)模型層來代替。3、點(diǎn)擊OK。PMWIN的界面發(fā)生變化，此時(shí)顯示了模型的整體網(wǎng)格。P

5、MWIN允許用戶改變和旋轉(zhuǎn)模型網(wǎng)格。用戶可以隨意改變模型網(wǎng)格行、列寬度，也可以增加刪除行和列。在該實(shí)例中，不需更改模型的網(wǎng)格。4、從File菜單中選擇LeaveEdit或點(diǎn)擊LeaveEdit按鈕退出。二、定義含水層類型二、定義含水層類型1、打開Grid菜單，選擇LayerType。LayerOptions對(duì)話框彈出。2、點(diǎn)擊標(biāo)簽為Type列中的某一單元格，單元格中將出現(xiàn)一個(gè)帶下箭頭的按鈕，點(diǎn)擊該按鈕，出現(xiàn)一個(gè)可供選擇含水層類型的列表。