Modflow操作截圖說明

1、一、創(chuàng)建文件及模型離散圖 模型主菜單主菜單屏幕包括以下選項(xiàng)：圖 模型主菜單選項(xiàng)（1）文件：包括打開、新建、保存、關(guān)閉Visual MODFLOW，導(dǎo)入、導(dǎo)出MODFLOW等。圖 文件菜單（2）點(diǎn)擊Files-New，選擇文件的儲存位置，選擇或新建文件夾，命名模型文件名，點(diǎn)擊確定，新建一個(gè)模型后，出現(xiàn)下圖窗口。該窗口用于選擇水流和溶質(zhì)運(yùn)移模擬所需的數(shù)值求解方法；水流和溶質(zhì)運(yùn)移模型參數(shù)的單位以及模型的起始日期和時(shí)間。圖 選擇模型程序其中，對于水流模型，可以選用下面的數(shù)值求解方法：USGS MODFLOW-96 from SWS（飽和，密度恒定）；USGS MODFLOW-2000 from SWS

2、（飽和，密度恒定）；USGS MODFLOW-2005 from SWS（飽和，密度恒定）；USGS SEAWAT from SWS（飽和，密度變化）；MODFLOW-AURFACT from HGL（變飽和度）。對于運(yùn)移模型，Visual MODFLOW 支持以下圖中程序，其中MT3DMS，MT3D99，RT3D運(yùn)移引擎支持多組分污染物運(yùn)移，而MT3Dv1.5僅支持單組分運(yùn)移。 圖 選擇模型程序（3）Next，水流設(shè)置在運(yùn)行類型處，可以設(shè)置非穩(wěn)定流模型和穩(wěn)定流模型，在非穩(wěn)定流模型中，要設(shè)置模型的起始日期和起始時(shí)間，輸入某些邊界條件數(shù)據(jù)時(shí)，就要定義時(shí)間數(shù)據(jù)。圖 水流模型默認(rèn)參數(shù)設(shè)置（4）Nex

3、t，運(yùn)移設(shè)置，在吸附下拉菜單中可選擇吸附模式，在反應(yīng)下拉菜單中選擇反應(yīng)類型。圖 溶質(zhì)模型默認(rèn)參數(shù)設(shè)置（5）模型區(qū)：用來創(chuàng)建模型網(wǎng)格，上面窗口用來導(dǎo)入模擬區(qū)地區(qū)，確定模型區(qū)維數(shù)，以及定義有限差分網(wǎng)格的行列數(shù)和分層數(shù)。圖 創(chuàng)建模型網(wǎng)格完成后，出現(xiàn)選擇模型區(qū)界面，確定導(dǎo)入底圖中的模擬區(qū)，并通過輸入X，Y對導(dǎo)入的圖片進(jìn)行校準(zhǔn)，要使圖片與坐標(biāo)系匹配，至少知道兩個(gè)已知坐標(biāo)的地理參考點(diǎn)，點(diǎn)擊地圖上的參考點(diǎn)，出現(xiàn)地理參考點(diǎn)對話框提示輸入所選位置的XY坐標(biāo)，將圖片附上坐標(biāo)信息，完成。圖 模型分割范圍設(shè)置和坐標(biāo)賦值二、Grid模塊-網(wǎng)格窗口圖 Grid模塊（1）View column，View row， View

4、 layer選擇查看類型。（2）Goto previous，next為選擇層數(shù)，即模型的第一層，第二層等（3）Inactive cells：活動單元的設(shè)置，自上向下分別為設(shè)置活動區(qū)域，設(shè)置非活動區(qū)域，設(shè)置活動單元格，復(fù)制活動單元格，復(fù)制活動區(qū)域。圖 定義有效、無效網(wǎng)格（4）Edit grid：Edit layer和Edit column，通過該部分可以添加、刪除、移動和加密剖分網(wǎng)格的設(shè)置?？筛鶕?jù)需要對網(wǎng)格進(jìn)行加密、減疏等操作。圖 編輯grid（5）Assign elevation，可以輸入某點(diǎn)，線，面，窗口的高程值。圖 高程值輸入（6）Coutouring-高程等值線的繪制，可以繪制頂板等值線

5、，厚度等值線和底板等值線，options是分別對上述三種等值線進(jìn)行設(shè)置。圖 繪制高程等值線（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)法，克力格，或反距離方法），默認(rèn)情況下，插值方法選擇自然鄰點(diǎn)法。圖 高程插值然后彈出下面的窗口，所需的數(shù)據(jù)由數(shù)據(jù)文件的列組成，并且與列號相匹配，而這個(gè)數(shù)據(jù)框用來與每一個(gè)數(shù)據(jù)區(qū)域相對應(yīng)，每一個(gè)數(shù)據(jù)域里鍵入適當(dāng)?shù)臄?shù)字

6、以使得列#1和X坐標(biāo)軸相對應(yīng)，使列#2和Y坐標(biāo)軸相對應(yīng)，以及列#3與高程相對應(yīng)，當(dāng)所有列都對應(yīng)后，next按鈕被激活。然后彈出有效數(shù)據(jù)檢測框，保證數(shù)據(jù)點(diǎn)位于下面坐標(biāo)系和單位框所確定的模型區(qū)域內(nèi)。圖17 點(diǎn)擊Finish進(jìn)入到坐標(biāo)系統(tǒng)和單位窗口，坐標(biāo)系統(tǒng)選擇模型坐標(biāo)，高程單位選擇米，完成高程導(dǎo)入。圖 坐標(biāo)系統(tǒng)和單位設(shè)置三、 Well模塊-井模塊：該模塊包括抽水井，水頭觀測井和濃度觀測井三個(gè)部分。圖 井模塊（1）抽水井左側(cè)工具欄按鈕是各種針對抽水井的選項(xiàng)，包括添加井、刪除井、編輯井、移動井和復(fù)制井。圖 抽水井設(shè)置添加井需將光標(biāo)移動到目標(biāo)位置，單擊鼠標(biāo)左鍵在該處添加一口抽水井，彈出一個(gè)新井窗口，提示

7、輸入井的詳細(xì)信息。包括井的名字和位置坐標(biāo)，ok。然后，在右上處輸入濾管頂部高程（top，m），底部高程（bottom，m）；輸入抽水時(shí)段，開始時(shí)間，結(jié)束時(shí)間和流量。圖 抽水井參數(shù)設(shè)置（2）水頭觀測井和濃度觀測井 Wells-head/concentration observation wells-import Obs.，導(dǎo)入數(shù)據(jù)，按彈出對話框設(shè)置對應(yīng)數(shù)據(jù)列，完成觀測水位數(shù)據(jù)的導(dǎo)入。 圖 觀測井在Observation points添加井號和觀測井的觀測層位，并在下面Observations輸入時(shí)間和觀測值。圖 觀測井設(shè)置四、Properties模塊-地下水流屬性的賦值圖 含水層參數(shù)設(shè)置 水流模

8、型中分別包括導(dǎo)水系數(shù)、儲水系數(shù)和初始水頭，每個(gè)系數(shù)設(shè)置分別包括對某區(qū)域的賦值，復(fù)制，數(shù)據(jù)庫，繪制等值線和導(dǎo)出數(shù)據(jù)部分。 通常在建模時(shí)，對模擬區(qū)進(jìn)行參數(shù)分區(qū)，分別對每個(gè)分區(qū)的參數(shù)進(jìn)行輸入。導(dǎo)水系數(shù)部分分別輸入x方向、y方向和z方向的導(dǎo)水系數(shù)。儲水系數(shù)部分分別輸入儲水系數(shù)、給水度、有效孔隙度和總孔隙度。初始水頭通常采用導(dǎo)入的方法，將事先準(zhǔn)備的text導(dǎo)入，同樣通過輸入列號，將數(shù)據(jù)與各列相對應(yīng)，然后進(jìn)入下一步，點(diǎn)擊插值（interpolate），完成。通過重復(fù)該步驟分別對每一層進(jìn)行初始水頭的賦值。在包括溶質(zhì)運(yùn)移的模型中，對模擬區(qū)初始濃度進(jìn)行賦值，同樣采用導(dǎo)入的方法，操作步驟與初始水頭的導(dǎo)入步驟相同。

9、五、Boundary模塊-邊界條件圖 邊界條件設(shè)置（1）Constant head （CHD）-定水頭邊界 Boundary-constant head-assign-line，分別輸入起始時(shí)間，結(jié)束時(shí)間，起始時(shí)間水頭和結(jié)束時(shí)間水頭，點(diǎn)擊ok，接受這些數(shù)值。網(wǎng)格單元上出現(xiàn)黑紅色，表明一個(gè)定水頭邊界條件已經(jīng)賦給了該單元。圖 定水頭邊界設(shè)置（2）River-河流邊界 選擇設(shè)置到適當(dāng)層，采用默認(rèn)導(dǎo)水系數(shù)公式，分別輸入結(jié)束時(shí)間，河水位，河床底高，河床厚度，河床滲透系數(shù)和河流寬度。若選中線性梯度，則分別輸入起始點(diǎn)和結(jié)束點(diǎn)的上述信息。圖 River 邊界設(shè)置（3）General head boundary

10、-通用水頭邊界分別輸入結(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-補(bǔ)給 輸入停止時(shí)間和補(bǔ)給速率。圖 Recharge模塊設(shè)置（8）Evapotranspiration-蒸發(fā) 輸入蒸發(fā)的停止時(shí)間，蒸發(fā)速率和深度。圖 蒸發(fā)設(shè)置（9）Constant co

11、ncentration-定濃度邊界 將污染源的位置和濃度進(jìn)行定義。選中污染源位置后分別賦值結(jié)束時(shí)間，污染源濃度。圖 定濃度設(shè)置（10）Recharge concentration和Evapotranspiration concentration分別設(shè)置污染物的補(bǔ)給濃度和蒸發(fā)濃度，輸入起始時(shí)刻、結(jié)束時(shí)刻、污染物濃度。六、Run-運(yùn)行模塊按F10回到主菜單，點(diǎn)擊Run進(jìn)入到運(yùn)行模塊。圖 模型運(yùn)行（一）Modflow運(yùn)行模塊圖 設(shè)置運(yùn)行參數(shù)（1）time steps，分別輸入其實(shí)時(shí)間，結(jié)束時(shí)間，時(shí)間步長。圖 時(shí)間步長設(shè)置（2）initial head圖 初始水位賦值選項(xiàng)（3）求解方法設(shè)置（默認(rèn)值）圖

12、 求解方法設(shè)置（4）Recharge（默認(rèn)）圖 補(bǔ)給設(shè)置（5）Layer 對各層承壓性進(jìn)行設(shè)置，layer type。圖 層承壓性設(shè)置（二）MT3D模塊圖 溶質(zhì)運(yùn)移模塊運(yùn)行參數(shù)設(shè)置（1）解法選項(xiàng)窗口圖 溶質(zhì)運(yùn)移求解方法設(shè)置（2）Output time-輸出時(shí)間MT3D輸出控制的默認(rèn)輸出設(shè)置是只在模擬結(jié)束時(shí)的濃度計(jì)算結(jié)果。在模型輸出結(jié)果中需要指定查看模擬結(jié)果的時(shí)間。輸入模擬時(shí)間，最大運(yùn)移步長，選擇指定輸出時(shí)間。圖 輸出時(shí)間設(shè)置（3）Run，選擇需要運(yùn)行的程序，點(diǎn)擊Translate&Run，開始運(yùn)行。圖 模型運(yùn)行（三）output-輸出可視化（1）顯示等值線圖 進(jìn)入結(jié)果輸出模塊，點(diǎn)擊Maps-conturing maps，可以選擇顯示不同類型的等值線圖，分別繪制水頭等值線，降深等值線，濃度等值線，水位高程等值線，水位深度等值線等。點(diǎn)擊Time，選擇等值線的輸出時(shí)間，并通過option選項(xiàng)，對不同等值線圖進(jìn)行顯示設(shè)置。圖 等值線繪制（2）時(shí)間序列擬合狀況 點(diǎn)擊Graphs-time series-head，點(diǎn)擊左側(cè)工具欄中show all time，顯示不同時(shí)刻該點(diǎn)計(jì)算水位與觀測