黃一帆，劉俊民，姜 鵬，張 斌

（西北農(nóng)林科技大學(xué) 水利與建筑工程學(xué)院，陜西 楊凌712100）

地下水是貯存于地球表面松散層孔隙或基巖裂隙中的水資源，是人類不可缺少的一種寶貴的資源，對(duì)人類生活、工農(nóng)業(yè)生產(chǎn)及城市建設(shè)都起著重要作用。同時(shí)地下水也是影響生態(tài)環(huán)境系統(tǒng)的一個(gè)重要因子之一，地下水的流量大小、水位高低變化往往會(huì)影響生態(tài)系統(tǒng)的天然平衡狀態(tài)。由于含水系統(tǒng)的空間變異性，求取水文地質(zhì)參數(shù)和查明水文地質(zhì)條件較為困難，并不能用以往的確定性方法，所以我們只能嘗試選擇一些不確定方法。本文欲在地下水建模預(yù)測(cè)中尋求精度更高的模型，因此利用時(shí)下最流行且預(yù)測(cè)精度最高的Visual Modflow 9.3軟件來(lái)進(jìn)行地下水的動(dòng)態(tài)預(yù)測(cè)，以求得更貼切實(shí)際情況、預(yù)測(cè)精度更高的方案。

1 涇惠渠灌區(qū)地下水概況

涇惠 渠 灌 區(qū) 位 于 108°34′34″—109°21′35″E，34°25′20″—34°41′40″N，南北寬大約70km，總面積約1 180km2，地勢(shì)總趨勢(shì)是自西北向東南傾斜，海拔450～350km，坡降約為3.33%～1.66%，地表及地下徑流排泄條件較好，地下水潛水平均埋深為15.8 m，最大埋深40.2m，第四季河流沖擊成的黃土狀亞黏土灌區(qū)廣為分布，土壤耕層容重、孔隙度較適宜。灌區(qū)屬大陸型半干旱季風(fēng)氣候，多年平均降水量512.5mm，降水時(shí)空分布極不均勻，7—9月降水量約占全年的50%～60%，年蒸發(fā)量為1 140.2 mm[1]。由于灌區(qū)受到地下水排泄及補(bǔ)給不均衡的影響，從而研究涇惠渠灌區(qū)地下水動(dòng)態(tài)的意義重大，其從總體趨勢(shì)可分為兩階段：灌區(qū)建立前及灌區(qū)建立后，本文著重研究灌區(qū)建立后。

灌區(qū)建立后，除受降水影響外，重要還是受到水事活動(dòng)的影響，其中特別是灌溉引用水地下水開采等活動(dòng)，這些活動(dòng)綜合制約并影響著地下水的動(dòng)態(tài)變化，說(shuō)明地下水動(dòng)態(tài)與人類活動(dòng)息息相關(guān)，文章分別選取涇陽(yáng)縣代表井256＃、254＃；三原縣代表井222＃、203＃、216＃；高陵縣代表井306＃；臨潼區(qū)代表井12＃、6＃、9＃以及閻良區(qū)代表井11＃進(jìn)行繪圖分析[2]（圖1，2）。

圖1 涇惠渠各區(qū)各代表井水位對(duì)比

圖2 涇惠渠2008年代表井水位變化過(guò)程線

由圖1分析看出，1984—2008年明顯地出現(xiàn)了地下水水位下降的現(xiàn)象，1984年在三原縣至閻良西區(qū)地下水水位較2002年下降1.47～23.43m，涇陽(yáng)縣處于整個(gè)地區(qū)水位較低的地帶，這與涇陽(yáng)縣地理及開采影響不可分割，2002年在閻良與三原縣交界處出現(xiàn)了較為明顯的地下水位下降趨勢(shì)，直到2008年水位沉降加深，并在涇陽(yáng)縣范圍埋深面積擴(kuò)大并且深度加深，1984—2008年的地下水水位有3.87～27.63m的下降幅度，下降幅度變大，說(shuō)明地下水位有逐年下降趨勢(shì)。由圖2可以看出，選取的幾口代表井中雖然地區(qū)不同，但是其在2008年變化趨勢(shì)基本一致，可以分析出總體趨勢(shì)是冬灌之后地下水位得到回升，再下降，到了4月份春灌以后又出現(xiàn)了一個(gè)回升，到了夏季，地下水位開始大幅度下降，而雨季的到來(lái)使得水位又開始回升?？傮w來(lái)看，灌溉入滲、人工開采和降雨入滲是影響該灌區(qū)地下水動(dòng)態(tài)變化的最主要的三個(gè)因素，而地質(zhì)因素的影響不明顯。為了進(jìn)一步分析涇惠渠地下水位沉降趨勢(shì)，因而對(duì)灌區(qū)地下水進(jìn)行軟件模擬，以此提出合理的控制方案。

經(jīng)過(guò)對(duì)涇惠渠灌區(qū)近年來(lái)地下水動(dòng)態(tài)現(xiàn)狀的分析，下文將對(duì)灌區(qū)深入進(jìn)行預(yù)測(cè)，基于地下水模擬軟件Visual Modflow結(jié)合GIS區(qū)域分析方法，通過(guò)流域內(nèi)灌區(qū)及周邊地下水賦存環(huán)境、運(yùn)動(dòng)規(guī)律、轉(zhuǎn)化條件、空間布局及動(dòng)態(tài)變化研究，以建立合理評(píng)價(jià)和量化地下水資源及可持續(xù)利用的方法、調(diào)控方式和指標(biāo)體系，為流域水資源優(yōu)化配置和綜合利用提供決策支持。

2 軟件介紹

Visual Modflow是目前國(guó)際上最流行、被各國(guó)同行公認(rèn)的三維地下水流模擬和溶質(zhì)運(yùn)移模擬的標(biāo)準(zhǔn)可視化的專業(yè)軟件系統(tǒng)[3－4]。該系統(tǒng)是加拿大 Waterloo Hydrogeologic Inc.通過(guò)在 Modflow軟件的應(yīng)用基礎(chǔ)上，運(yùn)用現(xiàn)代可視化技術(shù)而開發(fā)研制的，于1994年8月首次在全世界公開發(fā)行。高度集成的各類軟件包，包括用于地下水動(dòng)態(tài)模擬的 MODFLOW、粒子運(yùn)移軌跡和傳播時(shí)間模擬的 MODPATH、模擬污染物在地下水輸移過(guò)程中運(yùn)動(dòng)的MT3D，以及用于估計(jì)水文地質(zhì)參數(shù)和優(yōu)化的PEST，并具有直觀的、清楚的圖形交互界面。新穎的菜單結(jié)構(gòu)，以便于用戶對(duì)研究區(qū)進(jìn)行離散化和選擇有效的計(jì)算單元，確定邊界條件并對(duì)參數(shù)賦值，運(yùn)行及對(duì)模型進(jìn)行校正，還可用等值線或顏色陰影來(lái)實(shí)現(xiàn)結(jié)果的可視化，真正實(shí)現(xiàn)了人機(jī)對(duì)話。在模型的建立及結(jié)果顯示過(guò)程之中，模型網(wǎng)格、水文地質(zhì)參數(shù)和模擬結(jié)果，都可通過(guò)剖面圖或平面圖顯示。這個(gè)軟件的最大特點(diǎn)是把數(shù)值模擬過(guò)程中的各個(gè)步驟緊密有效地連接起來(lái)，從建模初始、輸入和修改各種水文地質(zhì)參數(shù)和幾何參數(shù)、模型運(yùn)行、反演校正各類參數(shù)，一直到最終的結(jié)果顯示輸出，使整個(gè)過(guò)程從頭到尾整體化、規(guī)范化。目前，這些特點(diǎn)是我國(guó)乃至世界上的同類軟件所不具備的。

3 地下水?dāng)?shù)值模擬

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

根據(jù)涇惠渠灌區(qū)含水層空間結(jié)構(gòu)可知，灌區(qū)含水層主要為第四系全新統(tǒng)Q4及上更統(tǒng)Q3中粗砂、砂礫石層構(gòu)成的潛水含水層。中更新統(tǒng)Q2及下更新統(tǒng)Q1黏土、亞黏土因富水性差可視為相對(duì)隔水層。由于全更新統(tǒng)上更新統(tǒng)含水層厚度較大、富水性強(qiáng)等特征，揭露部分深度即可滿足地下水開采量，灌區(qū)地下水徑流受整個(gè)第四系全新統(tǒng)Q4及上更新統(tǒng)Q3含水層控制。因此，涇惠渠灌區(qū)第四系含水層可概化為單層潛水含水層，模型中地下水徑流量代表第四系全更新統(tǒng)Q4及上更新統(tǒng)Q3含水層總徑流量，中更新統(tǒng)Q2及下更新統(tǒng)Q1作為相對(duì)隔水底板。研究區(qū)域的含水層是由潛水含水層和淺層承壓含水層組成，因?yàn)槌袎汉畬邮艿饺祟惢顒?dòng)的影響很小，所以，本次研究只針對(duì)潛水含水層進(jìn)行模擬研究。潛水含水層的厚度和埋深等性質(zhì)是隨著時(shí)間和空間的變化而變化的，在本次研究區(qū)域內(nèi)，潛水含水層基本上是從東到西、從兩側(cè)到中間逐漸變薄，故把潛水含水層可概化成為非均質(zhì)、各向異性非穩(wěn)定流的三維模擬。因此，研究區(qū)域的地下水遵循質(zhì)量守恒定律、能量守恒定律和達(dá)西定律。

3.2 地下水流數(shù)學(xué)模型

根據(jù)涇惠渠灌區(qū)水文地質(zhì)概念模型，將灌區(qū)淺層地下水概化為非均質(zhì)、各向異性非穩(wěn)定流的三維模擬單層潛水含水層，建立研究區(qū)數(shù)學(xué)模型：

式中：K——滲透系數(shù)（m／d）；h——含水層隔水底板至自由水面的距離，即含水層厚度（m）；ε——源匯項(xiàng)，即單位時(shí)間、單位面積垂向流入或流出含水層的水量，流入為正，表示匯，流出為負(fù)，表示源[m3／（d·m2）]；μ——潛水含水層的重力給水度；h0——水頭初始值（m）；n——第二類邊界的外法線方向；q——第二類邊界上單位時(shí)間、單位面積上側(cè)向流入或流出含水層的水量[m3／（d·m2））；Ω——滲流區(qū)域；S——自由面以下含水層儲(chǔ)水系數(shù)（1／m）；P——潛水面上的降水入滲和蒸發(fā)（m／d）；Γ0——滲流區(qū)域的上邊界；Γ1——滲流區(qū)域的第二類邊界；Kn——邊界法線方向的滲透系數(shù)（m／d）；

3.3 模型的建立

研究采用Visual Modflow 4.3求解上述數(shù)學(xué)模型，同時(shí)求解方法采用了預(yù)調(diào)共軛梯度法即PCG（Preconditioned Conjugate－Gradient）進(jìn) 行 迭 代 求解[5－7]。對(duì)包括涇惠渠灌區(qū)的水文地質(zhì)圖進(jìn)行掃描，然后采用ArcGIS 9.3進(jìn)行矢量化操作，保存為Visual Modflow 4.3可導(dǎo)入的shp文件并打開作為數(shù)值模型的底圖。根據(jù)研究區(qū)域的現(xiàn)實(shí)情況和本次研究精度的要求，針對(duì)圖3進(jìn)行單元格剖分時(shí)參考井點(diǎn)分布，對(duì)井點(diǎn)分布密集的部分進(jìn)行的單元格剖分要盡量緊密些，這樣考慮就相對(duì)多些。其中共有40行，53列，共計(jì)剖分了2 120個(gè)單元，其中白色單元是有效的，綠色單元是無(wú)效的。

本次研究中所應(yīng)用的涇惠渠灌區(qū)的地下水觀測(cè)資料、降雨、蒸發(fā)以及灌溉量和開采量等數(shù)據(jù)比較完整的是2008年和2009年，所以用2008年1月1日作為建立地下水?dāng)?shù)值模型并進(jìn)行模擬的初始時(shí)間。一般在模擬地下水動(dòng)態(tài)的情況下，給水度和滲透系數(shù)是兩個(gè)最主要的水文地質(zhì)參數(shù)，前面提到我們將涇惠渠灌區(qū)概化為非均質(zhì)的，故首先要按照不同的給水度和不同的滲透系數(shù)把研究區(qū)域劃分為7個(gè)子區(qū)域，見圖4，每個(gè)區(qū)域的參數(shù)初值見表1。

圖3 灌區(qū)網(wǎng)格剖分圖

圖4 水文地質(zhì)參數(shù)分區(qū)

表1 水力傳導(dǎo)系數(shù)和給水度初值

3.4 源匯項(xiàng)處理

源匯項(xiàng)是影響灌區(qū)地下水動(dòng)態(tài)變化的重要原因，反映著整個(gè)地下水系統(tǒng)的變化機(jī)理，所以，源匯項(xiàng)對(duì)于建立地下水?dāng)?shù)值模型是十分重要的。一般來(lái)說(shuō)，處理源匯項(xiàng)主要有兩種途徑，一種是以面狀強(qiáng)度的形式進(jìn)行補(bǔ)給或排泄，另一種是以點(diǎn)狀的形式參與模型的補(bǔ)給和排泄。但總體來(lái)說(shuō)，這兩種處理形式在活動(dòng)的單元格上是都參與的。

以2012年為例，灌區(qū)淺層地下水總補(bǔ)給量約為1.585 1億m3，其中降雨入滲補(bǔ)給量約占19.04%；渠滲漏補(bǔ)給量約占23.02%；灌溉入滲補(bǔ)給量約占12.89%；開采回歸補(bǔ)給量約占37.58%；側(cè)向補(bǔ)給量約占7.47%。地下水總排泄量約為2.209億m3，其中地下水開采量約占84.71%；潛水蒸發(fā)總量約占2.90%；側(cè)向排泄量約占2.91%；人、畜、城鎮(zhèn)工業(yè)用水量約占9.48%。補(bǔ)給總量與排泄總量之差為－6 238.99億m3，其水量是動(dòng)用了淺層地下水的靜儲(chǔ)備量中的一部分，即地下水疏干水量（表2）。

表2 各水平年地下水均衡量 萬(wàn)m3／a

3.5 模型的識(shí)別及驗(yàn)證

模型的識(shí)別和模型的驗(yàn)證通常要反復(fù)地進(jìn)行模型參數(shù)的修改和源匯項(xiàng)的修改，才能夠得到比較理想化的模擬結(jié)果，是整個(gè)建模過(guò)程中最為重要的一個(gè)步驟。本次研究對(duì)模型的識(shí)別和驗(yàn)證所采用的是擬合校正法，它是通過(guò)間接方法來(lái)反求模型的參數(shù)，比起直接求解模型參數(shù)具有更好的穩(wěn)定性。這種方法首先要給定模型參數(shù)的一組估計(jì)值，然后建立模型并進(jìn)行計(jì)算，再對(duì)比實(shí)測(cè)結(jié)果和模擬結(jié)果，算出誤差，如果模型的誤差精度沒(méi)有達(dá)到預(yù)先設(shè)定的要求，對(duì)原來(lái)模型的參數(shù)進(jìn)行適度的調(diào)整，再一次進(jìn)行計(jì)算，直到模型的計(jì)算結(jié)果和實(shí)測(cè)結(jié)果的誤差達(dá)到所要求的精度時(shí)，停止迭代，現(xiàn)在輸入模型的這一組參數(shù)認(rèn)為是最優(yōu)的。模型的識(shí)別驗(yàn)證遵循下面三條原則[8－10]：一是實(shí)際的地下水流場(chǎng)和模型模擬的地下水流場(chǎng)要保持基本一致，即在模擬期的某一時(shí)刻實(shí)際的地下水等值線和模型模擬的同一時(shí)間地下水等值線形狀大致相似；二是模型模擬的地下水動(dòng)態(tài)和實(shí)際過(guò)程中的地下水動(dòng)態(tài)大致相似，即實(shí)測(cè)的地下水位過(guò)程線和模型模擬的地下水位過(guò)程線大致一樣；三是模型最后應(yīng)用的參數(shù)要于十幾種的水文地質(zhì)條件基本符合。

從表3可以看出，所選的7口代表井中的實(shí)測(cè)值和模擬值之間的最大均方差為0.634 2m，最小均方差為0.402 1m，最大相對(duì)誤差為9.7%，最小相對(duì)誤差為7.2%，而相關(guān)系數(shù)則在0.695 0到0.948 7之間，實(shí)測(cè)值和模型模擬的水位值大體一樣，模型的精度較為精確。而位于模型邊界的區(qū)域和人類活動(dòng)較多的區(qū)域，模型的模擬精度相對(duì)較差，通過(guò)分析，主要由于以下幾個(gè)原因使得部分地區(qū)的模擬精度較差。

表3 2008－2009年代表井模型檢驗(yàn)

根據(jù)前面采用的滲透系數(shù)和給水度的初值以及對(duì)邊界條件的概化和對(duì)源匯項(xiàng)的計(jì)算，建立水文地質(zhì)三維模型，并對(duì)水文地質(zhì)參數(shù)進(jìn)行反復(fù)修改，最終確定了符合一定的精度要求的滲透系數(shù)和給水度，見表4。

需要說(shuō)明，所建立的模型無(wú)法完全準(zhǔn)確地反映并說(shuō)明地下巖層的結(jié)構(gòu)，再加上模型概化的邊界條件和計(jì)算的源匯項(xiàng)與實(shí)際情況也不可能完全一樣，所以模型的擬合結(jié)果越精確，并不能說(shuō)明本次模擬就具有很高的精度與可靠性，也就是說(shuō)，此次模擬允許存在個(gè)別沒(méi)有達(dá)到滿意擬合效果的檢測(cè)點(diǎn)的存在。

表4 水力傳導(dǎo)系數(shù)和給水度最終值

對(duì)模型計(jì)算結(jié)果的可靠性和代表性兩個(gè)方面進(jìn)行評(píng)價(jià)，進(jìn)行實(shí)測(cè)地下水含水層的流場(chǎng)和模型模擬含水層流場(chǎng)的比較。對(duì)地下水的流場(chǎng)進(jìn)行對(duì)比分析[11]，采用ArcGIS 9.3繪制模擬資料代表性最強(qiáng)的2009年9月的模擬地下水位等值線圖，如圖5所示，圖中面狀水位為實(shí)測(cè)2009年9月地下水位分布趨勢(shì)圖，圖上的等值線即為模擬值。由圖可知，實(shí)測(cè)和模型模擬的地下水位的總趨勢(shì)基本相同，但是對(duì)于富平縣、高陵縣與三原縣交匯處這些地區(qū)模型的模擬結(jié)果和實(shí)測(cè)值的誤差比較大，尤其是靠近309、305和208這幾口觀測(cè)井的附近區(qū)域變化很大，富平縣208號(hào)井附近由于資料相對(duì)較少，因而對(duì)總體趨勢(shì)的影響就相對(duì)較大，故其誤差就很大。對(duì)于觀測(cè)井相對(duì)密集的地區(qū)，如灌區(qū)東西部而言，誤差就相對(duì)較小，由于觀測(cè)井較多，對(duì)該地區(qū)的水位定位就越準(zhǔn)確，使得每個(gè)觀測(cè)井的測(cè)量值和模擬值之間的誤差不是很大，再加上每個(gè)水位的基礎(chǔ)數(shù)值很大，所以相對(duì)誤差會(huì)很小，這就最終表現(xiàn)為模型模擬的地下水流場(chǎng)和實(shí)測(cè)地下水流場(chǎng)在圖形上基本相似。

圖5 潛水層地下水實(shí)測(cè)流場(chǎng)與模擬流場(chǎng)對(duì)比（2009年9月）

4 提出預(yù)測(cè)方案并檢驗(yàn)

假定兩種不同情景的源匯項(xiàng)方案，將給定的源匯項(xiàng)數(shù)值輸入已經(jīng)建立驗(yàn)證的地下水三維模型，分別預(yù)測(cè)2013年該區(qū)域的地下水動(dòng)態(tài)變化。兩種動(dòng)態(tài)模擬情景如下：

方案1：在加大地表水的引用量并保持當(dāng)前地下水開采量不變的基礎(chǔ)上，增大干支、斗渠襯砌率，減少渠道系統(tǒng)滲漏量，提高渠系利用水效率，其他的參數(shù)條件不做改變，預(yù)測(cè)2014年的地下水動(dòng)態(tài)。

方案2：在研究區(qū)域灌溉面積不變的情況下，2009年后按照灌溉節(jié)水規(guī)劃每年增加節(jié)水灌溉面積并提高渠系利用水效率，2014年的人畜、城鎮(zhèn)工業(yè)用水量總用水量達(dá)到2 300萬(wàn)m3，地下水開采量達(dá)到15 000萬(wàn)m3，其他參數(shù)條件按固定值輸入，預(yù)測(cè)2014年的地下水動(dòng)態(tài)。

按照以上假定的兩種不同方案情景，將不同的源匯項(xiàng)數(shù)值輸入到已經(jīng)建立并驗(yàn)證的Visual Modflow模型中[12－15]，其他條件（邊界條件和水文地質(zhì)參數(shù)）不做改變，然后運(yùn)行模型，計(jì)算出不同情景下研究區(qū)域的地下水位。

圖6 兩種方案對(duì)比圖

由圖6可以看出，兩種不同方案下預(yù)測(cè)的2014年的地下水位比起2009年的地下水位（9月水位）都有所下降，但是方案一比方案二的地下水位略高，在其他源匯項(xiàng)和參數(shù)條件不變的前提下，方案二比方案一增加了地下水年開采量，從而導(dǎo)致方案二的地下水位比方案一的地下水位稍低，但是與建模的初始水位相差不大，并且變化趨勢(shì)基本一致，從而進(jìn)一步驗(yàn)證了模型的精確性。總體來(lái)說(shuō)，提高農(nóng)業(yè)用水效率。提高灌溉用水有效利用系數(shù)。加強(qiáng)農(nóng)業(yè)、工業(yè)及生活用水設(shè)備等措施，有助于發(fā)展節(jié)水型農(nóng)業(yè)，為以后該地區(qū)的地下水預(yù)測(cè)和規(guī)劃提供一定的參考依據(jù)。

5 研究存在的問(wèn)題及建議

由于本次研究水平有限，同時(shí)加上資料的完整性和時(shí)間的限制，研究中必然存在一定的問(wèn)題和不足，包括下面幾個(gè)方面：

（1）研究區(qū)域?yàn)闆芑萸鄥^(qū)，分析了該灌區(qū)的地下水動(dòng)態(tài)類型及特征，其主要受邊界條件、水文地質(zhì)參數(shù)和源匯項(xiàng)等因素的影響，研究區(qū)域越大，其關(guān)系越復(fù)雜，影響地下水動(dòng)態(tài)變化的因素也越多。

（2）研究在建模過(guò)程中只選擇建立了一層模型，沒(méi)有考慮其他含水層（承壓含水層和弱透水層），而在不同的區(qū)域潛水含水層的厚度也不一樣，為了研究方便，統(tǒng)一采用了一個(gè)厚度而沒(méi)有考慮坡度，使得地下水的排泄和補(bǔ)給誤差變大，影響了模型的精度。雖然在此方面不可避免地存在一定誤差，但在以后的研究過(guò)程中，應(yīng)盡量最優(yōu)地概化含水層，從而減小誤差。

（3）在處理邊界時(shí)，盡管將東、西、南、北各個(gè)邊界做了簡(jiǎn)單的概化，但是未能反映出其真實(shí)的自然情況，再加上數(shù)據(jù)的缺乏和更新，在一定程度上增加了模型的模擬誤差。

（4）在確定水文地質(zhì)參數(shù)（滲透系數(shù)和給水度）的過(guò)程中，都是根據(jù)以前的經(jīng)驗(yàn)值或試驗(yàn)結(jié)果，而缺少最新的資料和數(shù)據(jù)，從而在選取參數(shù)時(shí)就有一定的誤差存在，為了簡(jiǎn)化研究，把含水層概化為非均質(zhì)各向異性的三維流，通過(guò)把研究區(qū)域劃分為9個(gè)子區(qū)域輸入不同的初始參數(shù)，而沒(méi)有詳細(xì)地表現(xiàn)出其應(yīng)有的特征，進(jìn)一步增大了研究誤差。

本文通過(guò)對(duì)涇惠渠地下水進(jìn)行模擬并利用已有資料分析了當(dāng)前狀況，以此提出不同方案對(duì)灌區(qū)進(jìn)行地下水動(dòng)態(tài)預(yù)測(cè)，并從側(cè)面討論了灌區(qū)地下水的合理開發(fā)問(wèn)題。涇惠渠灌區(qū)是一個(gè)相對(duì)缺乏水資源的灌區(qū)，由于水資源的匱乏已成為灌區(qū)工農(nóng)業(yè)發(fā)展的障礙，但多年來(lái)地下水的過(guò)量開采，使得當(dāng)?shù)氐叵滤患眲∠陆?，出現(xiàn)很多枯井，造成大面積的地下水降落漏斗，亦導(dǎo)致了當(dāng)?shù)厮|(zhì)的惡化，所以有必要對(duì)灌區(qū)進(jìn)行分析評(píng)價(jià)。為提高水資源利用率，實(shí)現(xiàn)灌區(qū)水資源的可持續(xù)利用，應(yīng)該對(duì)灌區(qū)水資源進(jìn)行優(yōu)化配置，加強(qiáng)工農(nóng)業(yè)設(shè)備，有效提高地下水資源利用率，加強(qiáng)調(diào)控及管理制度，提出合理開發(fā)利用地下水資源的方案措施，加強(qiáng)灌區(qū)管理，完善水利設(shè)備，并健全法律及提高社會(huì)節(jié)水意識(shí)等都對(duì)灌區(qū)水資源合理開發(fā)利用有著重要的意義。

[1]陜西省地下水監(jiān)測(cè)資料 （2008）[M].陜西：陜西省地下水管理監(jiān)測(cè)局，2009.

[2]陜西省地下水監(jiān)測(cè)資料 （2009）[M].陜西：陜西省地下水管理監(jiān)測(cè)局，2010.

[3]白利平，王金生.GMS在臨汾盆地地下水?dāng)?shù)值模擬中的應(yīng)用[J].山西建筑，2004，30（16）：78－79.

[4]丁元芳，遲寶明，易樹平，等.Visual MODFLOW在李官堡水源地水流模擬中的應(yīng)[J].水土保持研究，2006，13（5）：99－102，105.

[5]陳南祥，姜新慧.基于GIS與層次分析法的地下水資源分區(qū)研究[J].人民黃河，2010（11）：60－62.

[6]宋豐驥，付金霞，常慶瑞.基于GIS陜西省區(qū)域經(jīng)濟(jì)發(fā)展環(huán)境的綜合評(píng)價(jià)[J].水土保持研究，2011，18（1）：37－41.

[7]Michael H A，Mulligan A E，Harvey C F.Seasonal oscillations in water exchange between aquifers and the coastal ocean[J].Nature，2005，436（7054）：1145－1148.

[8]郭曉東，田輝，張梅桂，等.我國(guó)地下水?dāng)?shù)值模擬軟件應(yīng)用進(jìn)展[J].地下水，2010，32（4）：5－7.

[9]Brodie R S.Integrating GIS and RDBMS technologies during construction of a regional groundwater model[J].Environmental Modelling ＆Software，1998，14（2）：119－128.

[10]伊燕平，盧文喜，張?jiān)?，?基于徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)的地下水?dāng)?shù)值模擬模型的替代模型研究[J].水土保持研究，2012，19（04）：265－269.

[11]葉劍鋒，劉小勇.基于GIS技術(shù)的地下水管理系統(tǒng)研究[J].水土保持研究，2011，18（3）：247－251.

[12]El－Kadi A I，Oloufa A A，Eltahan A A，et al.Use of a Geographic Information System in Site－Specific Ground－Water Modelinga[J].Groundwater，1994，32（4）：617－625.

[13]郝治福.石羊河流域鄧馬營(yíng)湖區(qū)地下水位動(dòng)態(tài)變化特征及數(shù)值模擬與預(yù)報(bào)[D].北京：中國(guó)農(nóng)業(yè)大學(xué)，2006.

[14]Kollarits S，Kuschnig G，Veselic M，et al.Decisionsupport systems for groundwater protection：innovative tools for resource management[J].Environmental Geology，2006，49（6）：840－848.

[15]李偉業(yè)，付強(qiáng).基于GIS的三江平原地下水資源評(píng)價(jià)[J].水土保持研究，2007，14（4）：92－95.