三維水流數值模型在地下水允許開采量評價中的應用

武秋娟 魯康橋

(鄭州市水利建筑勘測設計院，河南鄭州 450006)

三維水流數值模型在地下水允許開采量評價中的應用

武秋娟 魯康橋

(鄭州市水利建筑勘測設計院，河南鄭州 450006)

為查明趙口引黃灌區中牟境內深層地下水允許可開采量，采用地下水流數值模擬軟件GMS建立三維水流數值模型，對地下水資源進行評價計算。將開采井的開采層位根據其具體坐標和開采量加載在模型相應的網格上，在規劃開采條件下按運行20年預測模型，得出該水源地2030年末深層地下水預測流場，分析深層采井附近的降落漏斗。

數值模型;水資源;開采量

1 概述

河南省趙口引黃灌區位于黃河南岸豫東平原，始建于1970年，是以引黃灌溉為主要任務的大型灌區。灌區中牟縣境內趙口總干渠為土渠，未護砌，農業用水直接引自渠內黃河水。而作為當地生產生活用水取水水源則多為深層地下水，靠趙口總干渠及黃河側滲補給。為查明該區域作為城區水源地時深層地下水允許可開采量，采用地下水流數值模擬軟件GMS建立三維水流數值模型，對地下水資源進行評價計算。依據水源地的地形、地貌、水文地質條件和水源地開采后的影響范圍，評價計算的范圍:北至黃河水邊線，南至鄭開大道;東至瓦坡，西至萬灘，面積324.5km2。

2 水文地質條件

2.1 含水層組

趙口引黃灌區中牟境域地處黃河下游平原，本區地下水類型均為松散巖類孔隙水。

根據含水層的埋藏條件、成因類型、水力性質等，將區內550m以上劃分為淺層孔隙含水層組、中層孔隙含水層組和深層孔隙含水層組三個含水層層組，分別對應淺層水、中層水和深層水。淺層水系指埋藏于地表下100m以淺的含水巖組，中層水系指埋藏于地表下100～250m深度的含水巖組，250～500m為深層水。富水程度以單井涌水量為指標進行評價。單井涌水量淺層按降深5m和深層按降深15m進行換算。

2.2 富水性特征

2.2.1 富水區(3000～2000m3/d)

分布在勘查區盆王—賀崗—太平莊一線以西，水文地質井抽水試驗時，水位降深10.74m，單井涌水量1688.4m3/d，推算涌水量為2359.26m3/d。地下水礦化度＜1g/L，水位埋深19～21m。

2.2.2 中等富水區(2000～1000m3/d)

分布在勘查區盆王—賀崗—太平莊一線以東，水文地質井抽水試驗時，降深 17.59m，單井涌水量1680.24m3/d，推算涌水量為1432.63m3/d，該區地下水礦化度＜1g/L，水位埋深20～21m。

3 地下水補給、徑流與排泄

3.1 淺層地下水

淺層地下水的補給以大氣降水入滲為主，其次為黃河測滲及其他河流、渠道滲漏、農田灌溉回滲及地下徑流補給。

淺層地下水徑流總體流向為自西北向東南徑流。淺層地下水的排泄主要是開采和蒸發及側向徑流，本區河流洪水期水位高補給地下水，枯水期地下水位高向河流徑流排泄。鄉鎮以居民生活用水開采為主，外圍以農業、灌溉開采為主。

3.2 中層地下水

中層水補給主要來自北部、西北部側向徑流補給。另外淺層與中層之間無良好隔水層，且由于本區農村集中供水均采用中層水，造成中層水位低于淺層水位，故淺層水對中層水具有一定的補給作用。

中層地下水的徑流，由西北、北部向東南、南部徑流。排泄主要是農村集中供水開采及向深層越流排泄。將來大量開采淺層地下水時，有可能向淺層越流排泄。

3.3 深層地下水

深層地下水主要來自西北部的徑流補給，深層與中層之間弱透水層較薄，且中層水位高于深層，故中層水對深層水有一定的補給作用。深層水的徑流方向由西北向東南。排泄主要是側向徑流排泄。

4 地下水動態特征

4.1 淺層地下水動態特征

潛水蒸發是淺層地下水保持均衡狀態的一個重要因素，隨著年開采量的增加，淺層水位下降，蒸發排泄量相應減少，直至遇到豐水年和豐水季節，接受降水補給再度回升。根據淺層地下水位的變化過程劃分出兩種動態類型:入滲—開采型與入滲—蒸發型。

4.2 中層地下水動態特征

中層地下水來源于側向徑流及淺層越流補給，消耗于開采和側向徑流及向深層越流排出。由于人工開采，地下水位已降至8～12m左右，因中層開采在本區均布，區內未形成明顯降落漏斗。故其動態類型為徑流—開采型。其變化主要受開采量大小控制。

4.3 深層地下水動態特征

深層地下水天然狀態下來源于側向徑流補給，消耗于側向徑流排出，工程區深層開采量較小。所以，本區深層地下水動態，仍大致保持著天然狀態，年變幅1m左右。隨著人工開采量的增加，將變為徑流-開采型。

5 數值法評價地下水允許開采量

5.1 地下水系統概念模型

5.1.1 含水層結構特征

本區地下水類型均為松散巖類孔隙水。

模型將勘察區在空間上大致劃分為淺層孔隙含水層組(100m以淺)、中層孔隙含水層組(100～250m)和深層孔隙含水層組(250～500m)三個含水層組。

5.1.2 地下水流動特征

淺層地下水的徑流受地形和補給源控制。由于黃河現行河道是下游平原的中脊和地表水、地下水的分水嶺。所以，在黃河影響帶內，黃河以南地下水自西北向東南徑流;黃河以北地下水自西南向北徑流。由于黃河在東壩頭北折，地下水流向逐漸變為西北。黃河影響帶外，由于受西部山區的側向徑流補給，黃河北側地下水自西南向東北測流，黃河南側自西向東逐漸變為東南方向。水力坡度和地形坡降相近，屬徑流滯緩類型。

5.1.3 邊界條件的概化

根據區內水文地質條件和地下水開發利用情況，將地下水系統模擬區定為:模型的北部以黃河為界，南邊界至鄭開大道，西至萬灘，東到瓦坡村，總面積324.5km2。

a.側向邊界:對于淺層含水層，北邊界黃河定為River邊界;參照區域內多年平均的流場圖將東部和南部邊界定為通用水頭邊界，將西部定為隔水邊界，對于中層和深層含水層，將模型各邊界均定義為通用水頭邊界。

b.垂向邊界:潛水含水層自由水面為系統的上邊界，通過該邊界，潛水與系統外發生垂向水量交換，如接受大氣降水入滲補給、河渠補給、田間入滲補給、蒸發排泄等。三個含水層組之間通過越流交換物質和能量，其越流量由含水層在垂向上的滲透系數及含水層的厚度決定。

c.水力特性:地下水系統符合質量守恒定律和能量守恒定律;含水層分布廣、厚度大，在常溫常壓下地下水運動符合達西定律;考慮淺、深層之間的流量交換以及軟件的特點，地下水運動可概化成空間三維流;地下水系統的垂向運動主要是層間的越流，三維立體結構模型可以很好地解決越流問題;地下水系統的輸入、輸出隨時間、空間變化，故地下水為非穩定流;參數隨空間變化，體現了系統的非均質性，而且在水平和垂直方向上各不相同，所以參數概化成各向異性。

綜上所述，分析評價區可概化成非均質各向異性、空間三維結構、非穩定地下水流系統，即地下水系統的概念模型。

5.2 模型預測條件

地下水數值預測模型的模擬期為20年(2011年1月至2030年12月)，以一個月為一個應力期，共240個應力期，其初始流場為2010年模擬期末的流場，如圖1。

圖1 預測模型深層初始流場

根據降雨的周期性，選取歷史上1980—1999年的降雨量作為預測模型20年的降雨量，該20年的降雨量中有7個枯水年份，尤其是最后5年中包含了3個枯水年份，可以達到預測極端氣候下地下水變化趨勢的要求。

以2010年補給、排泄量作為預測期間的源匯項，在枯水年份相應地增加農業和生活用水開采量，在豐水年份相應地減小農業和生活用水開采量。水源地開采井的開采層位根據其具體坐標和開采量加載在模型相應的網格上。

5.3 模型預測結果

在規劃開采條件下運行20年預測模型，得出該水源地2030年末深層地下水預測流場如圖2。由圖中可以看出，在深層的開采井附近形成了局部的降落漏斗，但水位下降幅度并不大，影響范圍也在開采井周圍的區域內。水源地深層含水層在預測期末的水位降深如圖3所示。

圖2 規劃開采條件下預測期末深層流場

5.4 允許開采量分析

開采資源量的預報，應在現有開采技術條件下，以不影響已有水源地正常運行和工、農業用水等為原則，限定預報開采約束條件為:?根據現有開采技術和本區水文地質條件，模型預測水源地開采主井最大水位降深不應超過含水層厚度的1/3;?根據實際取水能力，由開采引起的主井水位降深不大于設計降深;?預測開采量不應超過區域可開采資源量;?開采狀態下，保證淺層地下水水位不可下降太快。

在以上約束條件下，通過模擬計算得出雁鳴湖水源地一期深層的開采量如表1所列，深層允許開采量為9980m3/d，完全可以滿足規劃開采條件下9980m3/d的開采要求。

表1 水源地允許開采量模擬計算

圖3 規劃開采20年后模擬區深層地下水位降深

5.5 水源地開采量組成分析

對比深層含水層未開采條件下和開采9980m3/d條件下預測模型20年間地下水多年平均水量均衡情況，可以得到深層地下水各補排量的變化如表2所列。

表2 預測模型深層地下水多年平均補排量變化情況

由表可知，水源地深層含水層在規劃開采條件下，動用了儲存量，主要激發了邊界側向流入量和中層含水層的越流補給量，且研究區邊界的側向流出量明顯減少。水源地規劃開采量(365×104m3/a)的來源組成分析如表3，由此可見水源地深層含水層的規劃開采量是有保證的。

表3 深層含水層開采量來源組成分析

6 結語

地下水流數值模擬計算結果表明:水源地深層地下水開采9980m3/d時，水源地開采20年地下水位降落漏斗中心水位降深為6.02m，水源地開采主要依靠激發地下水側向流入量及越流補給量，襲奪地下水側向流出量，以滿足地下水位處于含水層頂板之上。因此，水源地深層井的開采是有保證的。

Application of Three-dimensional Water Flow Numerical Model in Groundwater Allowable Exploitation Evaluation

WU Qiujuan，LU Kangqiao
(Zhengzhou Water Conservancy Construction Survey and Design Institute，Zhengzhou450006，China)

Groundwater flow numerical value stimulation software GMS is adopted for establishing three-dimensional water flow numerical value model in order to identify the amount of deep groundwater allowable exploitation in Zhaokou Yellow River Diversion Irrigation Area.Groundwater resources can be evaluated and calculated.The mining layer of exploration wells can be loaded on corresponding network of the model according to specific coordinates and exploitation.Deepest layer groundwater predication flow field of the water source at the end of 2030 is obtained under planned exploitation condition according to operation of 20-year forecast model.Depression funnel near deep exploitation well can be analyzed.

numerical value model;water resources;exploitation quantity

TV211.1+2

:A

:2096-0131(2015)01-0064-05