Visual MODFLOW在地下水模擬中的應用分析

李強+于維坤+李騰飛+李偉華

【摘要】：采用Visual MODFLOW預測諸暨某廠區建成投產后流場的變化，建立地下水數值模型。在此基礎上，利用Visual MODFLOW軟件中的MT3DMS 模塊，選取CODCr為溶質進行污染物運移模擬，預測地下水中污染物的運移情況。結果表明：擬建項目正常運行情況下，生產污水對地下水水質的影響很小；在非正常工況下，廢水泄漏進入地下水，地下水中CODCr濃度達到最大值為16mg/L，最大遷移距離為180m。

【關鍵詞】：地下水；數值模擬；Visual MODFLOW；CODCr

1、當前地下水流和溶質運移的研究中，常用的方法有解析法、數值法和實驗法。水文地質領域的數值模擬方法主要有有限差分法、有限元法、邊界單元法和有限分析法。地下水模擬在地下水研究中發揮著越來越重要的作用，目前流行的地下水模擬軟件主要有GMS、Visual MODFLOW、FEFLOW、TOUGH2等。筆者采用Visual MODFLOW軟件模擬研究諸暨某擬建廠區投產后廢水泄漏對地下水的影響。

1.1 模擬范圍

綜合考慮擬建項目周圍的區域地形地貌特征、水文特征、地質條件、水文地質條件和周圍的地下水環境敏感目標等因素，結合本次實地勘探和水質分析結果，模擬計算范圍約為1.8km2。根據區內鉆孔資料，研究區地面高程范圍為15～20m，淺層含水層大多集中在-3～-1m，其下為相對隔水層。該區域面積較小，水文地質條件差異性不大，可視為單一性質含水層。該區域含水層主要巖性為粉質黏土，滲透性能較差，含水層滲透系數為0.5m/d。

1.2 網格剖分

考慮到模擬精度尤其是溶質遷移模型精度的要求，根據模擬區水文地質條件，在垂向上將模擬區劃為一層；各層南北長1.0km，東西長1.7km，在水平方向上用正交網格進行剖分，網格數目為100×200，單個網格大小為10m×8.5m。

1.3 邊界條件

利用觀測水位資料，確定研究區中部渠道為河流邊界，并確定西部和南部定水頭邊界。將研究區概化為具有穩定的空間結構，地下水位連續三維非均值各向同性的非穩定流概念模型。

1.4 地下水數值模型

式中，μs表示儲水率（1/m）； h表示含水層水位高程（m）；t表示模擬時間段內的時間（d）； K表示滲透系數（m/d）；W表示源匯項（m3/d）； h（x，y，z，t）表示邊界上的已知水位函數（m）；q（x，y，z，t）為第二類邊界流量函數（m3/d·m）； k表示三維空間上的滲透系數張量；Γ1為一類邊界；Γ2為二類邊界；n為邊界Γ2的外法線方向。

1.5 污染運移模擬

1.5.1 潛在污染源分析

項目的廢水主要為染整廢水（包括前處理、染色）、地面沖洗水、設備冷卻水、蒸汽冷凝水和員工生活污水。各股廢水可劃分為高濃度生產廢水、低濃度生產廢水及生活污水三種。項目產生的廢水按分質收集后，分別經現有廢水處理設施預處理，其中高濃度生產廢水經預處理達標后納管；低濃度生產廢水部分經深度處理后通過供水系統回用；生活污水經化糞池、隔油池預處理后與生產廢水一道納管。

1.5.2 污染物遷移情景設置

正常狀況下，廢水經過現有廢水處理設施處理后，納入諸暨市污水處理廠集中處理。擬建項目對地下水環境可能產生的影響，主要為非正常狀況。本次模擬設置模擬情景為非正常狀況生產發生泄漏進入地下水含水層，預測廢水中污染物在含水層中的濃度變化、影響范圍和超標情況。

廠區污水處理站人工防滲層局部發生泄漏，污染物以200L/d的泄漏速度進入地下水含水層。泄漏污水中主要污染物濃度為CODCr2500mg/L。持續泄露7d，后經檢修發現破裂后修補，污水不再滲入地下水。模擬總時長為5000d。

1.5.3 污染物運移數學模型

根據評價區地下水流實際情況和污染物運移的一般規律，可建立以下數學模型來表示污染物進入評價區含水層后在地下水中的遷移過程：

式中：R遲滯系數為1； 為土壤人孔隙率； 為組分濃度（mg/L）； 為彌散系數（m2/d）； 為地下水速度張量；W為水流的源匯項。

聯立地下水流方程和污染物運移方程求解即可獲得污染物在含水層中的濃度分布數據。本次采用數值模擬方法對聯立的數學模型進行計算。污染物運移過程的模擬，將在Visual Modflow軟件建立的水流數值模型的基礎上，疊加該軟件中的MT3D模塊進行。

1.5.4 模擬預測

模擬非正常狀況下，0-7d內，污染物持續進入地下水中，7d時，泄漏停止，此時泄漏中心點污染物濃度達到最大濃度為16mg/L。隨著時間的推移和水流運動，污染暈以泄漏點為中心，向四周擴散。100d時，污染暈最大濃度為4mg/L。500d時，污染暈最大濃度為0.9mg/L。1000d時，污染暈最大濃度為0.6mg/L，5000d時，污染暈最大濃度為0.06mg/L。該區域水力坡度較大，地下水流交互作用強度一般，污染暈沿著水流方向最大遷移距離約180m。

1.5.5 污染風險預測分析

在非正常狀況下，高濃度生產廢水防滲層出現破裂，廢水中的CODCr泄露排放進入地下水，不考慮包氣帶的滯留作用、包氣帶和飽和帶對污染物的消減作用、污染物的自然降解作用等。該地區滲透性能較差，水力坡度較大，地下水流交互作用一般。5000d模擬期內地下水中CODCr濃度達到最大值為16mg/L，最大遷移距離約180m。污染暈范圍較小，易于控制，采取垂直防滲墻和原位異位修復能有效控制和削減地下水中的濃度。

【參考文獻】：

[1]蔚東升，劉潔，馮兆東. Visual MODFLOW在國內的應用[J]. 能源與節能，2014，（06）：96-97+126.

[2]李希萌. 黑河中游地下水數值模擬模型[D].清華大學，2014.

[3]楊彪，羅周全，王益偉，謝乘煜. 基于Visual Modflow的礦山地下水流場分析及預測[J]. 礦冶工程，2013，（04）：11-15+21.