Visual Modflow在水資源論證中的應用探討

嚴芳芳

（菏澤市水文中心，山東 菏澤274000）

中宏工貿有限公司取用地下水，需要進行建設項目水資源論證，以申請取水許可。在進行水資源論證時，最佳取水方案的制定及取水可靠性和影響范圍的量化分析是重點難點。本次基于Visual Modflow，建立項目廠區附近地下水數學模型，根據項目取水需求制定開采方案進行數值模擬，直觀地展示出影響區域及影響程度。

1 區域地質、水文地質條件

區域全部為第四系松散沉積物覆蓋，出露上更新統（Q3）和全新統（Q4），是河流沖洪積物構成，厚度由南向北逐漸增厚，相應巖石的粒度由粗變細，巖石層次逐漸增多，堆積物類型由簡單變復雜。

項目區主要含水層為第四系松散巖類孔隙水含水巖組，該項目取水井含水層為全新統沖積孔隙含水層（QY），地層具二元結構，上部為砂質黏土或黏質砂土夾姜石及砂礫石透鏡體，底部為砂礫石。似層狀姜石層及砂礫石層是淺層地下水的主要含水層，其透水性和富水性較好，單井涌水量一般300～1 000 m3/d。目的含水層主要接受大氣降水入滲補給和河道垂向測滲補給，同時接受區外的側向徑流補給及農田灌溉回歸補給。

2 抽水試驗

Visual Modflow的應用，需要項目區的滲透系數、導水系數等水文地質參數，為此進行了單孔及群孔抽水試驗。

2.1 單孔抽水

項目廠區現有自備井5眼，周邊選取觀測孔5眼。2014年7月21日，以該公司1#自備井作為抽水井，2#作為觀測井，進行了單井抽水試驗。單井出水量48 m3/h，抽水170 min后水位達到穩定，抽水井最大降深14.82 m，穩定抽水24 h，結束后25 min內回升到抽水實驗前狀態。經計算，滲透系數k=1.28 m/d，含水層導水系數為64 m2/d，抽水影響半徑為168 m，降深14.82 m,單位涌水量為0.90 L/（s·m）。

2.2 群孔抽水

單井抽水試驗結束并恢復水位后，以生產平均需水量及最大小時需水量為基礎，以1#、2#、3#井為抽水井進行了群井抽水試驗，抽水流量分別為43 m3/h、25 m3/h、25 m3/h，水井合計流量為93 m3/h。群井抽水持續時間62 h，在3眼供水井同時開采的情況下，東西影響距離為487 m，南北影響距離為550 m，影響面積為0.84 km2。

3 VisualModflow的應用

3.1 模型建立與修正

Visual Modflow是目前應用較廣的溶質運移模擬評價和三維地下水流的可視化標準專業軟件系統，主要由Modpath、Modflow、MT3D三個模塊組成，可以便捷地為各剖分單元和邊界條件進行賦值，快速輸出模型區域和劃分計算單元，具有極強的成果可視化輸出能力。

本次對項目區域進行了第四系含水層水文地質模型的概化，設置為周圍無地表水體補給的邊界條件，采用有限差分法，對所模擬的區域以矩形網格進行剖分，對重點區域周邊網格進行加密剖分。由于廠區較小，且含水層為第四系，小范圍內含水層橫向差別較小。采用Visual Modflow軟件，根據已有的地質、水文地質參數，建立均質各向同性的三維非穩定流模型，經識別修正，調整滲透系數k、單位儲水系數等參數后，模擬結果與實測抽水試驗數據擬合較好。以此模型，對開采后區域地下水流場、降落漏斗等情況進行預測分析。

3.2 最優取水方案及取水影響

本次在滿足生產需水前提的條件下，對5眼開采井制定不同的開采布局，將多種開采方案導入Visual Modflow模型，模擬不同開采條件下的影響范圍及影響程度，以對該項目取水可能造成的不利影響進行數據化的反應，同時獲取對周邊環境及其他用水戶影響最小的最優開采布局，見表1。

表1 不同開采方案影響情況

通過模擬可知：該項目在需水66.7 m3/h條件下，以1#、3#為開采井，2#、4#、5#井作為備用井定期輪流開采，設計開采流量1#為41.7 m3/h，3#為25 m3/h。該方案既能保證項目供水，又可使影響范圍遠離北部其他企業。

開采井初始地下水埋深分別為：1#為23.20 m、2#為22.72 m，按設計方案開采穩定后，東西影響半徑484 m，南北影響半徑為520 m，影響面積為0.79 km2，其中降深大于1 m的面積為0.078 km2。設計開采方案下，會造成廠區西南部8#農用井水位下降1.5 m左右，需采取一定經濟補償措施，但對北部企業的開采井幾乎無影響。

4結論

利用Visual Modflow軟件建立模型后，對多種開采方案進行逐一分析驗證，選取滿足企業用水需求且影響最小的方案，并量化給出影響度數值，是一種地下取水水資源論證中實用的分析方法，對科學的水資源管理有一定的實用價值。但該模型是基于一定的假定水文地質條件下的，在開采井的長期運行過程中，需根據真實情況調整開采方案，以降低取水影響。