基于 Visual Modflow的邯鄲縣地下水位動態(tài)研究

張大龍，閆 衛(wèi)，于虎廣，金 喆

(1.河北工程大學，河北 邯鄲 056038;2.邯鄲慧龍電力設(shè)計研究有限公司，河北 邯鄲 056035;3.中國建筑材料工業(yè)地質(zhì)勘查中心山東總隊，山東 濟南 250100)

邯鄲縣素有“糧倉、棉?！敝Q，是一典型的城郊縣，城鄉(xiāng)交流的不斷加快，極大地促進了縣城經(jīng)濟的快速發(fā)展。然而，淡水資源短缺、地下水超采［1］、水供需矛盾等問題也相繼出現(xiàn)，嚴重地制約了當?shù)厣鐣?jīng)濟的發(fā)展。邯鄲縣地區(qū)，工業(yè)、農(nóng)業(yè)用水水源主要來源于淺層地下水，為盡可能地滿足用水需要，大量開采地下水，致使地下水位連年大幅度下降，淺層地下水嚴重超采區(qū)不斷加大，深層淡水上層被疏干［2－4］，若不加以控制，勢必造成深層淡水的枯竭［5］。面對嚴峻的地下水資源問題，只有積極地克服各種困難來縮小超采區(qū)面積，控制地下水降落漏斗的不斷擴張，有效地緩解地下水資源嚴重短缺的矛盾，才能維護和建立高效的供水環(huán)境。為此，本次研究利用Visual Modflow軟件進行邯鄲縣地下水流的數(shù)值模擬，并預測2015、2020年縣域淺層地下水設(shè)計開采量方案下的地下水位動態(tài)變化情況，為研究區(qū)域的地下水資源合理開發(fā)利用提供科學依據(jù)

1 區(qū)域概況

邯鄲縣地處太行山東麓，地勢是西高東低，境內(nèi)分為丘陵區(qū)和山前平原區(qū)。多年平均降水量為507.9 mm，降水量的年內(nèi)分布極不均勻，全年降水量的70% ～80%集中在6－9月份。天然河流有滏陽河及其支流渚河、沁河、輸元河，人工河道有支漳河分洪道［6］。邯鄲縣全部行政區(qū)域總面積為439.7 km2。其中，山丘區(qū)面積 124.0 km2，平原區(qū)面積 315.7 km2。根據(jù)縣境內(nèi)水文地質(zhì)特征，可以將其分為丘陵區(qū)和山前平原區(qū)兩個水文地質(zhì)區(qū)，兩個水文地質(zhì)區(qū)又可進一步詳 分，見表1。

表1 邯鄲縣水文地質(zhì)詳情

2.1 水文地質(zhì)概念模型

水文地質(zhì)概念模型是研究區(qū)域?qū)嶋H水文地質(zhì)條件的概化［7］。研究區(qū)地下水主要來源于淺層地下水，也是本次模擬研究的目的層。區(qū)域含水層在垂向上沒有明顯的、穩(wěn)定的隔水層;地下水流以水平運動為主，地下水流系統(tǒng)的輸入、輸出隨時間、空間變化;因巖性的不同水文地質(zhì)參數(shù)有所區(qū)別，但在方向上差異不明顯。綜上，本次將區(qū)域含水層概化為非均質(zhì)各向同性二維非穩(wěn)定流模型。

東部邊界有地下水位觀測井，概化為一類已知邊界;南部邊界、西部邊界可概化為隔水邊界;北部邊界與永年交界，有大量地下水流入模擬區(qū)，概化為二類流量邊界。邯鄲縣與市區(qū)邊界中西部和北部邊界因有輸元河的補給，概化為二類流量邊界;東部邊界為排泄邊界，概化為二類流量邊界。研究區(qū)上部邊界定為有水量交換的邊界，下部存在泥巖組成，其滲透性極差，定為隔水邊界。

2.2 地下水流數(shù)學模型及求解方法

依據(jù)滲流連續(xù)性方程、達西定律，上述地下水流系統(tǒng)可用以下偏微分方程及其定界條件進行描述，如式 (1)所示:

式中:H為地下水水頭，m;T為含水層導水系數(shù)，m∕d;μ*為含水層給水度(無量綱);W為垂向水量交換項，m/d;Q為為地下水開采量抽水率，m3/s;D為計算區(qū)域;Γ1為已知水頭的第一類邊界;Γ2為已知流量的第二類邊界;H0(x，y)為初始時刻水頭分布函數(shù)，m;φ1(x，y，t)為第一類邊界上的已知水位函數(shù)，m;q(x，y，t)為第二類邊界上的已知流量函數(shù)，m3/s;t為時間。

圖1 模擬區(qū)參數(shù)分區(qū)圖

運用Visual Modflow的前承條件共軛斜量程序包PCG2進行求解。將區(qū)域剖分為3 484個矩形網(wǎng)格，大小為500×500 m，有效單元格為1 740個，如圖1所示，圖中藍色矩形網(wǎng)格為無效單元格，境內(nèi)分為丘陵區(qū)和山前平原區(qū)，其中，根據(jù)平原區(qū)淺層含水層巖性的差異，又將其細劃分為3個水文地質(zhì)分區(qū)，參數(shù)分區(qū)圖如1所示。

2.3 模型識別與驗證

選擇2010年4月27日至5月27日作為識別模型參數(shù)K和μ*的時段。利用各觀測孔的水位資料，繪制出地下水的等水位線圖，得到模型計算的初始流程(圖2)。

圖2 地下水初始流場圖

在識別過程中，采用手動和自動相結(jié)合的方法，通過實際水位和計算水位的擬合分析［8］，如果計算水位與實測水位相差很大，則根據(jù)參數(shù)變化范圍和實際水位差值，重審研究區(qū)域水文地質(zhì)分區(qū)的劃分情況，重新給定一組參數(shù)，直至二者擬合較好為止。選取觀測點4個(S52，J73，S69，S83)，經(jīng)過反復的調(diào)整，繪制觀測孔水位擬合曲線，如圖3所示?？梢钥闯鰯M合效果較好，可以滿足解決實際問題的需要，參數(shù)識別結(jié)果如表2所示。

表2 水文地質(zhì)參數(shù)識別結(jié)果

為了全面驗證模型，通過其他時段的實測資料，對模型的可靠性進行進一步的驗證。選取2010年10月10日至2010年11月10日作為模型的驗證時段。同樣，選取上述的4個觀測井，進行水位擬合，驗證結(jié)果如圖4。

觀測孔水位擬合曲線表明，模型的建立符合研究區(qū)的實際水文地質(zhì)情況，地下水動態(tài)過程一致，水位擬合達到了精度要求，可以運用該模型對研究區(qū)域地下水流系統(tǒng)進行預測。

3 地下水未來動態(tài)預測

為有效控制邯鄲縣地下水位逐年下降的趨勢，緩解部分地區(qū)地下水供需失衡，并根據(jù)邯鄲縣地下水限采目標和南水北調(diào)完成情況，以調(diào)節(jié)均衡型開采動態(tài)為原則，控制水位降落過程，給定開采初值，通過模型計算不斷調(diào)整開采量，直到開采動態(tài)達到采補平衡［10－12］，由此計算并制定 2015、2020年縣域淺層地下水設(shè)計開采方案，如表3所示。按照地下水的限采量，運用上述建立的地下水流數(shù)值模型，預測2015、2020年(可供水量和需水量均考慮為75%保證率情況下)的地下水位情況。按設(shè)計開采量進行開采(可供水量和需水量均考慮為75%保證率情況下)，應用上述建立的模型進行計算，得到2015、2020年研究區(qū)地下水位情況，如圖5、6所示。

圖3 觀測孔水位擬合曲線

表3 邯鄲縣未來水平年淺層地下水設(shè)計開采量成果表104m3

從預測結(jié)果的水位等值線圖可以看出，2011年至2015年地下水位有明顯下降趨勢，水位相比下降了3 m;2015年以后淺層地下水位平均降深一直穩(wěn)定在3.5 m左右，流場趨于穩(wěn)定。

圖4 觀測孔水位擬合曲線

圖5 2015年地下水位等值線圖

4 結(jié)語

(1)依據(jù)水文地質(zhì)條件、水資源開發(fā)利用等情況，將模擬區(qū)域進行參數(shù)分區(qū)，在參考各個勘察和研究階段所進行的抽水試驗成果的基礎(chǔ)上，確定各分區(qū)的水文地質(zhì)參數(shù)初值，模型識別時，應通過自動與手動相結(jié)合的方式，多次地調(diào)整水文地質(zhì)參數(shù)，得到了能夠較真實反映實際水文地質(zhì)條件的參數(shù)。

(2)通過模型的識別與驗證，可知建立的水文地質(zhì)概念模型和地下水流數(shù)學模型符合邯鄲縣實際情況，可用于該研究區(qū)域地下水開采的預測及評價。

(3)對邯鄲縣境內(nèi)設(shè)定了未來水平年淺層地下水設(shè)計開采量，可較好地控制地下水位的持續(xù)下降問題。2011年至2015年地下水位有明顯下降趨勢，2015年以后，地下水位平均降深穩(wěn)定。

圖6 2020年地下水位等值線圖

［1］聶振龍，陳江，王金哲，嚴明疆，田言亮.地下水在京津唐區(qū)域社會經(jīng)濟發(fā)展中的作用［J］.干旱區(qū)資源與環(huán)境.2011，25(10):75 －79.

［2］張兆吉，費宇紅，陳宗宇，等.華北平原地下水可持續(xù)利用調(diào)查評價［M］.北京:地質(zhì)出版社.2009.

［3］方茜，任曉力.河北省地下水開發(fā)利用狀況研究［J］.河北建筑科技學院學報.2005，22(1):4 －5.

［4］王凱軍，曹劍峰，徐蕾，孫濤.地下水資源管理預警系統(tǒng)的建立及應用研究—以長春城區(qū)為例［J］.水科學進展.2005，29(3):238－243.

［5］王貴玲，陳浩，藺文靜，劉花臺，范琦，張薇.河北省京津以南平原區(qū)未來30年地下水供需預測［J］.干旱區(qū)資源與環(huán)境.2006，，20(6):63 －68.

［6］高佳.邯鄲市東部平原區(qū)水資源優(yōu)化配置研究［D］.河北工程大學.

［7］沈媛媛，蔣云鐘，雷曉輝等.地下水數(shù)值模擬中人為邊界的處理方法研究［J］.水文地質(zhì)工程地質(zhì).2008，12(6):12 －15.

［8］盧文喜.地下水系統(tǒng)的模擬預測和優(yōu)化管理［M］.北京:科學出版社.1999.

［9］GB/T 14497－93.地下水資源管理模型工作要求［S］.

［10］盧文喜，鄒立芝.降水量的隨機模擬［C］.環(huán)境地學問題研究論文集.北京:石油工業(yè)出版社.1995，24 －32.

［11］陶月贊，曹健康，祝京.淺層地下水水位動態(tài)模擬與中長期預報［J］.安徽地質(zhì).2002，12(1):56 －58.

［12］?；鄱?需水預測研究講述［J］.四川師范大學學報.1996，23(1):104－109.