MODFLOW軟件入門詳細使用教程（三）模型運行與結果分析1.運行模型：命令行：進入模型文件目錄，執(zhí)行`mf6.exe-smf6_simple.nam`（`nam`為MODFLOW-6的命名文件）；GMS：點擊`Run`→`MODFLOW`，選擇模型文件后運行。2.結果解讀：水頭文件（`.hds`）：用`flopy.utils.HeadFile`或GMS的`Contour`工具查看，重點關注抽水井、邊界附近的水頭降深是否合理；流量文件（`.cbc`）：分析補給、排泄、井的流量平衡，驗證“輸入補給量+邊界流入量=井抽水量+邊界流出量”是否成立（誤差需<5%）。四、核心模塊（包）深度解析掌握關鍵包的參數(shù)邏輯，是精準建模的核心：（一）`DIS`包：模型離散化的“骨架”空間離散：`nlay`（層數(shù)）、`nrow`/`ncol`（行列數(shù)）、`delr`/`delc`（網(wǎng)格尺寸）決定模型的空間分辨率；`top`/`botm`定義每層的垂直范圍，需與實際含水層厚度匹配（潛水層`botm`可設為含水層底板，承壓層需考慮隔水頂/底板）。時間離散：`perlen`（應力期長度，如1年）、`nper`（總應力期數(shù)）、`nstp`（每個應力期的時間步數(shù)），需平衡計算精度與效率（非穩(wěn)定流需小時間步，穩(wěn)態(tài)流可設`nper=1`）。（二）`NPF`包：水流運動的“引擎”滲透系數(shù)（`k`/`k33`）：`k`為水平滲透系數(shù)，`k33`為垂向（層間）滲透系數(shù)，需通過抽水試驗、文獻類比確定（單位：m/d或m/s，注意統(tǒng)一）；儲水屬性：潛水層需結合`sy`（給水度），承壓層需結合`ss`（貯水率），兩者物理意義不同（`sy`是重力釋水，`ss`是彈性釋水）。（三）`RCH`包：補給的“源頭”面狀補給（`rech`）：均勻補給速率，需結合降水入滲、灌溉回歸等數(shù)據(jù)（如年降水600mm，入滲系數(shù)0.3，則`rech=0.6*0.3/365≈0.0005m/d`）；分層補給：若不同層補給率不同，需用`rech`數(shù)組（`(nlay,nrow,ncol)`）單獨設置。（四）`WEL`包：井的“開關”應力期數(shù)據(jù)：`stress_period_data`為字典或列表，格式為`[層,行,列,流量]`，正流量為注水，負為抽水；井的位置：行、列索引從0開始，需與`DIS`包的網(wǎng)格對應（可在GMS中通過“點選”確定坐標）。五、實例操作：典型潛水含水層抽采模擬場景設定某矩形區(qū)域（1000m×1000m）為單層潛水含水層，地表高程20m，底板高程0m（厚度20m）。滲透系數(shù)K=10m/d，給水度Sy=0.2。區(qū)域均勻補給（`rech=0.001m/d`），四周為定水頭邊界（`H=10m`），中心有一口抽水井（抽水量1000m3/d）。步驟1：GMS可視化建模（適合非編程用戶）1.創(chuàng)建模型：`File`→`New`→`MODFLOWModel`，命名為`well_model`。2.網(wǎng)格與層設置：`ModelGrid`→設置`Rows=10`、`Cols=10`、`Delr=100`、`Delc=100`；`Layers`→`AddLayer`，設置`Top=20`、`Bottom=0`。3.邊界條件：`MODFLOW`→`BoundaryConditions`→`ConstantHead`，框選四周網(wǎng)格，設置`Head=10`。4.井與補給：`MODFLOW`→`Wells`→`AddWell`，點擊中心網(wǎng)格，設置`Flux=-1000`；`Recharge`→設置`Recharge=0.001`。5.運行與可視化：`Run`→`MODFLOW`，運行后點擊`Results`→`Contour`，查看水頭等值線（抽水井附近應出現(xiàn)“漏斗狀”降深）。步驟2：Flopy代碼建模（適合編程用戶）參考“三（二）”的代碼，模型運行后，可通過`matplotlib`或GMS查看水頭分布。若需調(diào)整參數(shù)（如增大抽水量、加密網(wǎng)格），只需修改代碼中`well_data`的`flux`或`delr`/`delc`，重新運行即可觀察結果變化。六、常見問題與優(yōu)化技巧（一）模型不收斂（“Convergencefailure”）原因1：求解器迭代次數(shù)不足。解決：在`NewtonOptions`中增大`max_iteration`（如從100→300），或啟用`under_relaxation`（欠松弛）加速收斂。原因2：參數(shù)不合理（如滲透系數(shù)突變、補給量過大）。解決：檢查參數(shù)單位（如K誤設為1000m/d，實際應為10m/d），或平滑參數(shù)分布（如用`NPF`的`k`數(shù)組漸變）。（二）計算速度慢優(yōu)化1：網(wǎng)格優(yōu)化。重點區(qū)域加密，非重點區(qū)放大（如抽水井附近用10m網(wǎng)格，外圍用100m）。優(yōu)化2：