盧志華

(1.丹華水利環(huán)境技術(shù)(上海)有限公司，上海200235;2.上海市環(huán)境科學(xué)研究院，上海200233)

貴州省畢節(jié)市黔西縣碳酸鹽巖分布廣，巖溶發(fā)育，地下水資源豐富，且多數(shù)以巖溶泉水出露，水質(zhì)好，是理想的生活水源。由于自然條件限制，當(dāng)?shù)厮Y源開發(fā)利用條件差，開發(fā)利用成本較大，農(nóng)村人口占總?cè)丝诒戎卮?，且居住較為分散，地下水是城鎮(zhèn)供水和農(nóng)村人畜飲水的重要水源。

近年來，黔西縣社會經(jīng)濟(jì)的快速發(fā)展對地下水環(huán)境造成一定影響，但是當(dāng)?shù)貙Φ叵滤谋Ｗo(hù)和管理卻比較薄弱。黔西縣煤礦資源豐富，煤礦開采工程建設(shè)改變了含水層結(jié)構(gòu)，影響地下水的補(bǔ)、徑、排關(guān)系，礦坑排水和尾礦堆積也造成地下水體污染，破壞生態(tài)環(huán)境。農(nóng)村生活垃圾和農(nóng)業(yè)化肥施放造成的面源污染，隨降雨入滲進(jìn)入地下含水層，污染地下水體。同時，黔西縣屬巖溶山區(qū)，地表水和地下水交互強(qiáng)烈，地表水受到污染后也會影響地下水的質(zhì)量。因此，為防止地下水源受到污染，保障當(dāng)?shù)仫嬎踩?，在黔西縣開展地下水保護(hù)區(qū)試點(diǎn)研究，為區(qū)域地下水的管理、保護(hù)、監(jiān)控、科學(xué)預(yù)測、合理開發(fā)和可持續(xù)利用提供基礎(chǔ)。

本研究將通過數(shù)值模型手段，調(diào)查研究黔西縣地下水狀況和動態(tài)變化規(guī)律，為重要地下水水源地劃分保護(hù)區(qū)范圍，提出地下水保護(hù)區(qū)綜合保護(hù)措施，為當(dāng)?shù)氐叵滤Y源合理開發(fā)利用和綜合管理提供理論依據(jù)和決策工具。

1 黔西縣自然和水文地質(zhì)條件概述

黔西縣位于貴州中部偏西北、烏江中游鴨池河北岸，東經(jīng)105°47'～106°26'、北緯 26°45'～27°21'，總面積 2381 km2，東部以六廣河為界與貴陽市修文縣相鄰，南以鴨池河、六沖河為界連接清鎮(zhèn)市和織金縣，西部以支嘎阿魯湖和西溪河為界與大方縣相連，北和東北與大方縣、金沙縣接壤，地處亞熱帶溫暖濕潤氣候區(qū)，多年平均氣溫13.9°C，多年平均降雨量1005 mm，降雨主要集中在夏季，四季不分明，居住著漢、彝、苗、仡佬、布依、白和滿等18個民族，2010年末全縣總?cè)丝跒?7.11萬人，其中農(nóng)業(yè)人口79萬人，農(nóng)業(yè)人口占全縣人口比重約為91%。黔西縣有著豐富的煤炭資源。

黔西縣內(nèi)地層自寒武系中上統(tǒng)至第四系均有出露，以三疊系分布面積最廣，其間缺失奧陶系、泥盆系及石炭系中上統(tǒng)地層，其余各地層受構(gòu)造影響呈條帶狀、透鏡狀局部分布。侏羅系集中在東北部中平及西部洪家渡水庫庫尾一帶，寒武系僅在東北部中建邊境出露。黔西縣境內(nèi)地下水主要賦存于大面積出露碳酸鹽巖溶含水層，包括棲霞—茅口組、燈影組、清虛洞—下奧陶統(tǒng)、永寧鎮(zhèn)組、關(guān)嶺至法郎組、龍?zhí)吨链舐〗M和夜郎組等含水巖組，巖性主要為灰?guī)r、白云巖及其過渡類型巖石，屬巖溶透水層，約占全縣國土面積的92%，局部條帶狀分布基巖、碎屑巖裂隙含水層，包括牛蹄塘至金頂山組、飛仙關(guān)組、龍?zhí)吨链舐〗M和須家河組等含水巖組，巖性有粘土巖夾煤層、碎屑巖、玄武巖等，為相對隔水層。不同構(gòu)造的復(fù)合部位，出現(xiàn)地下水局部富集。另外，部分區(qū)域存在碎屑巖與碳酸鹽巖互層沉積，形成局部地下水承壓。

黔西縣多年平均水資源量12.4億m3。采用降水入滲系數(shù)法計(jì)算多年平均地下水資源量，黔西縣地下水資源量為3.09億m3，占總水資源量的24.9%，其中野紀(jì)河區(qū)地下水資源量1.88億m3，占該區(qū)總水資源量的31.7%，六沖河區(qū)地下水資源量0.56億m3，占該區(qū)總水資源量的24.6%，鴨池河區(qū)地下水資源量0.65億m3，占該區(qū)總水資源量的25.9%。

黔西縣地下水類型主要為重碳酸－鈣型或者重碳酸－鈣、鎂型，礦化度一般小于0.5g/l，總硬度1.04～25.07，PH值介于6～7.9之間。除極少數(shù)特殊水質(zhì)(溫泉)含有有害元素錳、汞、氟以外，區(qū)域內(nèi)水質(zhì)良好，可作為工農(nóng)業(yè)供水和生活用水。

2 地下水?dāng)?shù)值模型

本研究中通過建立研究區(qū)域地下水?dāng)?shù)值模擬模型，定量模擬研究區(qū)域地下水的流場分布和動態(tài)變化過程，分析重點(diǎn)關(guān)注區(qū)域地下水的補(bǔ)給、徑流和排泄情況，為地下水資源評估和地下水水源地保護(hù)區(qū)劃分提供理論依據(jù)。

2.1 地下水概念模型

建立地下水流場數(shù)值模擬模型，需要對實(shí)際的水文地質(zhì)條件進(jìn)行概化，建立研究區(qū)域水文地質(zhì)概念模型。概化包括確定模型的范圍、邊界概化、含水層水文地質(zhì)結(jié)構(gòu)、源匯項(xiàng)、水文地質(zhì)參數(shù)分布和地下水徑流特性等。根據(jù)對研究區(qū)域自然、氣象、水文和地質(zhì)條件調(diào)查，確定地下水模型建模區(qū)域?yàn)榘蠓?、黔西兩縣縣境的野紀(jì)河流域、六沖河流域和赤水河流域部分區(qū)域，北部以分水嶺和行政區(qū)劃為界，南部和東部分別以六沖河和烏江為界(如圖1所示)，總面積為7835 km2。根據(jù)對研究區(qū)域水文地質(zhì)條件分析，將其地下水系統(tǒng)概化為單一的非均質(zhì)各向同性穩(wěn)定流淺層地下水流系統(tǒng)，由于巖性、地下水賦存形式和水動力特征的空間差異，對水文地質(zhì)參數(shù)進(jìn)行空間分區(qū)描述。

圖1 模型研究區(qū)域范圍

研究區(qū)域地下水的補(bǔ)給來源包括降水入滲補(bǔ)給、河流滲漏補(bǔ)給、農(nóng)業(yè)灌溉回滲補(bǔ)給和側(cè)向地下徑流補(bǔ)給，其中降水入滲是主要的補(bǔ)給來源。地下水的排泄包括潛水蒸發(fā)、向地表水排泄、側(cè)向徑流排泄和人工開采，其中地下水人工開采和向地表水排泄占主導(dǎo)地位。研究區(qū)域地表水和地下水交互強(qiáng)烈，相互補(bǔ)給和排泄，將模型區(qū)的河網(wǎng)水系概化為第三類柯西邊界條件。模型北部大都為基巖山區(qū)，主要接受大氣降水補(bǔ)給，將其概化為隔水邊界。局部河谷區(qū)的側(cè)向徑流概化為第二類紐曼流量邊界條件。

2.2 地下水?dāng)?shù)值模型建立

在水文地質(zhì)概念模型的基礎(chǔ)上，建立研究區(qū)域地下水?dāng)?shù)值模型。

2.2.1 網(wǎng)格剖分

本研究中采用不規(guī)則三角網(wǎng)格進(jìn)行空間網(wǎng)格剖分，網(wǎng)格剖分時除需遵循一些剖分基本原則外，如三角形單元內(nèi)角盡量不要出現(xiàn)鈍角和相鄰單元間面積相差不應(yīng)太大等，還應(yīng)充分考慮如下實(shí)際情況:1)充分考慮工作區(qū)的邊界、巖性分區(qū)邊界、行政分區(qū)邊界等;2)觀測井、水源地開采井位置盡量放在剖分單元的節(jié)點(diǎn)上;3)流場變化大和重點(diǎn)研究區(qū)域(如水源地)需要局部網(wǎng)格加密。本項(xiàng)目將模擬區(qū)域剖分為35271個單元，共18165個節(jié)點(diǎn)。模型網(wǎng)格剖分結(jié)果見圖2。

圖2 研究區(qū)域網(wǎng)格剖分結(jié)果

2.2.2 側(cè)向徑流

側(cè)向徑流邊界的徑流量估算需綜合考慮含水層地下水的流場及其水力坡度、含水層滲透系數(shù)、厚度和邊界長度。根據(jù)達(dá)西定律，邊界側(cè)向徑流量通過以下公式計(jì)算:q側(cè)補(bǔ)=KMIB，式中q側(cè)補(bǔ)地下水側(cè)向流量(104m3/d);K斷面附近的含水層滲透系數(shù)(m/d);I垂直于斷面的水力坡度;B斷面寬(km);M含水層厚度(m)。

2.2.3 河道滲漏和基流補(bǔ)給

模型使用第三類柯西邊界描述地下水含水層和地表水體之間的水量交換，包括河流向地下含水層的滲漏過程(河流水位高于地下水位)和地下水補(bǔ)給地表水過程(地下水位高于河流水位)(圖3)。河道和地下含水層間的側(cè)向滲漏和補(bǔ)給通過公式:q河渠滲漏=Φ·(Href－H)計(jì)算，式中Φ≈kf/d為河道滲漏和基流補(bǔ)給導(dǎo)水系數(shù)，Href為參考河流水位，kf為河床的滲透系數(shù)，d為河床淤泥層的厚度。

模型區(qū)內(nèi)地表河網(wǎng)中共布有8個水位監(jiān)測站，收集到這些水位站的長系列動態(tài)觀測資料，每個節(jié)點(diǎn)的水位值由這些水位實(shí)測數(shù)據(jù)線性插值獲得。圖4為河網(wǎng)水位邊界賦值結(jié)果。河道滲漏和基流補(bǔ)給導(dǎo)水系數(shù)描述了河床淤泥層的導(dǎo)水性能，通常地下水補(bǔ)給河道基流過程其導(dǎo)水性能稍高于河床淤積過程(河道滲漏、地表水補(bǔ)給地下水)，模型計(jì)算過程中會根據(jù)地下水流補(bǔ)排方向自行選擇相應(yīng)參數(shù)。

圖3 地下水含水層和地表水體間的水量交換

2.2.4 地下水開采

研究區(qū)內(nèi)的地下水開采主要用于城鎮(zhèn)生活用水、農(nóng)村人畜飲水、工礦企業(yè)用水和農(nóng)業(yè)灌溉用水，其中農(nóng)村人畜飲水占重要比重。圖4為研究區(qū)域2008年農(nóng)村人畜飲水工程和幾個重要的城鎮(zhèn)供水水源地開采井位置分布圖及地下水資源利用量估值。本研究中綜合評估了獲取的各種水資源利用統(tǒng)計(jì)數(shù)據(jù)和相應(yīng)規(guī)劃報告對地下水開采量進(jìn)行了估算。

圖4 研究區(qū)域2008年農(nóng)村人畜飲水工程和重要城鎮(zhèn)供水水源地開采井位置分布圖及地下水開采量

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

水文地質(zhì)參數(shù)是表征含水介質(zhì)儲水能力、釋水能力和地下水運(yùn)動能力的指標(biāo)，其值主要參考供水水文地質(zhì)手冊的經(jīng)驗(yàn)值、相關(guān)的水文地質(zhì)報告和抽水試驗(yàn)數(shù)據(jù)獲得。本研究中根據(jù)地下水賦存類型和含水層巖性作為水文地質(zhì)參數(shù)分區(qū)主要依據(jù)。模型中兩個主要水文地質(zhì)參數(shù)滲透系數(shù)和給水度的空間分區(qū)結(jié)果見圖5，其值參見表1。

圖5 水文地質(zhì)參數(shù)滲透系數(shù)和給水度空間分區(qū)

表1 水文地質(zhì)參數(shù)值

2.2.6 降水入滲補(bǔ)給

降水入滲是研究區(qū)域地下水的主要補(bǔ)給來源，降水入滲量與降水量、潛水水位埋深和包氣帶巖性有關(guān)。降水入滲補(bǔ)給量計(jì)算公式為:Qj=α.F.X，式中Qj降水入滲補(bǔ)給量(104m3/a);α年降水入滲系數(shù);F計(jì)算區(qū)面積(106m2);X年降水量(m/a)。

研究區(qū)域?yàn)楦咴絽^(qū)，降水空間分布不均勻，年內(nèi)和年際變化較大。本研究中對模型區(qū)內(nèi)所有雨量站降雨長系列觀測資料進(jìn)行了頻率分析，并分別估算了不同保證率下的典型豐、平、枯水年的降雨量值。降雨量的空間分區(qū)通過泰森多邊形法劃分，分區(qū)結(jié)果見圖6。研究區(qū)域降雨入滲補(bǔ)給系數(shù)根據(jù)巖性和地下水賦存類型進(jìn)行分區(qū)取值(如表2所示)。圖7為研究區(qū)域豐水年降雨入滲補(bǔ)給量空間分布。

圖6 降雨量空間分區(qū)結(jié)果

圖7 研究區(qū)域豐水年降雨入滲補(bǔ)給量空間分布

潛水蒸發(fā)量主要與潛水位埋深、包氣帶巖性、地表植被和氣候因素有關(guān)，一般認(rèn)為水位埋深大于4m的地區(qū)潛水蒸發(fā)很小。由于研究區(qū)域?qū)儆趲r溶山區(qū)，地下水埋深較深，通常介于50－200m之間，因此本研究中潛水蒸發(fā)量忽略不計(jì)。

表2 降雨入滲補(bǔ)給系數(shù)

2.3 模型識別和結(jié)果分析

模型建立之后需要進(jìn)行模型參數(shù)識別，通過反復(fù)調(diào)整部分模型參數(shù)、邊界條件和源、匯項(xiàng)等使得模型模擬結(jié)果能反映實(shí)際的水文地質(zhì)條件和地下水流況。模型識別的原則包括:識別后的模型概化和參數(shù)應(yīng)符合實(shí)際的水文地質(zhì)條件;模擬的地下水流場與實(shí)際流場一致，可以反映區(qū)域地下水流向;模擬的地下水位變化過程與趨勢與實(shí)際的地下水動態(tài)過程一致。

本研究中收集到研究區(qū)域內(nèi)部分地下水位觀測井單次觀測地下水位數(shù)據(jù)(2011年)。對所建地下水穩(wěn)定流數(shù)值模型進(jìn)行了模型識別和流場擬合。圖8為模擬地下水位等值線與實(shí)測地下水位觀測數(shù)據(jù)擬合結(jié)果。擬合結(jié)果顯示由模擬地下水流場與實(shí)測結(jié)果吻合較好，所建立的地下水?dāng)?shù)值模型基本可以反映研究區(qū)域地下水流向和流場變化趨勢。研究區(qū)域地下水流受地形影響，總體趨勢為由西北部山區(qū)向東南河谷地區(qū)匯集，局部流場受地質(zhì)構(gòu)造和巖性控制。

圖8 模擬和實(shí)測地下水等值線比較

同時，對研究區(qū)域各地下水補(bǔ)給和排泄水均衡要素進(jìn)行了統(tǒng)計(jì)分析，表3列出了研究區(qū)域地下水均衡項(xiàng)多年平均量統(tǒng)計(jì)結(jié)果。研究區(qū)域地下水多年平均總補(bǔ)給量為9.86億m3，多年平均總排泄量為9.94億m3，均衡誤差為0.8%，滿足穩(wěn)定地下水流系統(tǒng)水量平衡條件。由表3可以看出，降雨入滲是地下水主要補(bǔ)給來源，且地下水系統(tǒng)主要向河流排泄水量，補(bǔ)給河道基流，地下水開采和側(cè)向徑流排泄量均較小。

表3 研究區(qū)域地下水均衡多年平均量統(tǒng)計(jì)結(jié)果

圖9 研究區(qū)內(nèi)重要飲用水供水水源地流場模擬和流線追蹤結(jié)果

2.4 模型模擬和預(yù)測

利用所建立的地下水動力數(shù)值模型對研究區(qū)域地下水系統(tǒng)進(jìn)行模擬，通過對流場流線追蹤確定不同匯流期間的地下水集水(補(bǔ)給)區(qū)域，以此作為地下水水源地保護(hù)區(qū)的劃分依據(jù)。由于受到動態(tài)降雨入滲、河流水位、地下水開采等因素影響，地下水系統(tǒng)是一個動態(tài)系統(tǒng)，水源地的補(bǔ)給區(qū)域也隨時間動態(tài)變化。水源地保護(hù)區(qū)的劃分應(yīng)根據(jù)最不利情況下地下水流場和和其補(bǔ)給區(qū)域確定。本研究模擬了研究區(qū)域豐水年降雨條件下的地下水流場和其匯流補(bǔ)給區(qū)，然后根據(jù)不同地下水匯流時間對應(yīng)的補(bǔ)給范圍分級劃分保護(hù)區(qū)。圖9為豐水年(降雨保證率20%)條件下研究區(qū)域重要地下水水源地流場模擬和流線追蹤結(jié)果。

3 地下水保護(hù)區(qū)劃分及保護(hù)措施

中國地下水資源保護(hù)區(qū)劃分一般采用三級控制原理，綜合考慮水源地地面影響范圍和地下水補(bǔ)給集水區(qū)范圍，將環(huán)境污染控制在不同保護(hù)區(qū)范圍之外。地下水水源保護(hù)區(qū)邊界的劃定主要考慮研究區(qū)水文地質(zhì)條件，特別要考慮含水層性質(zhì)、地下水運(yùn)動特征、地下水形成特征、地下水循環(huán)系統(tǒng)特征以及地下水含水層的覆蓋層性質(zhì)和含水層巖性因素，同時還要考慮污染物遷移和病菌、微生物生長條件以及其它人為因素等。

根據(jù)《中華人民共和國水污染防治法》和國家環(huán)保局2007年2月1日頒布實(shí)施《飲用水水源保護(hù)區(qū)劃分技術(shù)規(guī)范》，通常在水源地分別劃分一級保護(hù)區(qū)、二級保護(hù)區(qū)和三級保護(hù)區(qū)。一級保護(hù)區(qū)是指包括井群在內(nèi)的一定范圍，其作用主要是保證集水要有一定的滯后時間，防止污染物的滲透入，保證水源地的出水安全。二級保護(hù)區(qū)在一級保護(hù)區(qū)外，其作用是為輕度污染的水提供足夠的時間采取補(bǔ)救措施。

三級保護(hù)區(qū)在二級保護(hù)區(qū)外緣一側(cè)，防止污染水體流入二級區(qū)，從而保護(hù)供給水源。根據(jù)地貌單元和水文地質(zhì)條件可劃出集水區(qū)的范圍，也可根據(jù)計(jì)算補(bǔ)給帶和邊界寬度來劃定范圍。

本研究中應(yīng)用FEFLOW建立研究區(qū)域地下水?dāng)?shù)值模型，根據(jù)地下水流場不同匯流時間集水區(qū)的范圍，對研究區(qū)內(nèi)的所有重要飲用水水源地進(jìn)行了保護(hù)區(qū)分級劃分，并針對每個水源地的具體情形提出相應(yīng)的保護(hù)措施。本文以響水巖水源地為例進(jìn)行闡述。

3.1 響水巖地下水源地保護(hù)區(qū)劃分及保護(hù)措施

響水巖地下水水廠位于黔西縣以東大關(guān)鎮(zhèn)，共解決大關(guān)和鐘山兩鎮(zhèn)約1.5萬人的飲水問題。響水巖出露地層表層為0.2～5m第四系坡積褐黃色粘土，下部為三疊系白云巖。該區(qū)域大氣降水為地下水主要補(bǔ)給水源，其強(qiáng)度和大小影響地下水的徑流過程。降水到達(dá)地面通過三疊系白云巖、灰?guī)r溶蝕隙、構(gòu)造裂隙及洼地豎井以入滲和滲入方式補(bǔ)給地下水。

3.2 水源地現(xiàn)存問題

水源地為天坑，地下水流量大，周圍缺少圍擋保護(hù);水源地管理規(guī)則不健全，缺少嚴(yán)格的水資源管理制度;水源地附近農(nóng)業(yè)活動劇烈，農(nóng)業(yè)面源污染防治、水土流失控制是主要問題;泉水出流區(qū)緊靠公路，來往交通可能對水體水質(zhì)造成影響;缺少水量、水質(zhì)監(jiān)測。

3.3 保護(hù)區(qū)劃分

響水巖水源地保護(hù)區(qū)劃分示意圖見圖10，取水點(diǎn)周邊20m范圍劃分為一級保護(hù)區(qū)范圍，一級保護(hù)區(qū)外流場流線追蹤匯流時間小于50天的集水區(qū)范圍為二級保護(hù)區(qū)，二級保護(hù)區(qū)外匯流時間小于20年的集水區(qū)范圍劃分為三級保護(hù)區(qū)。

3.4 保護(hù)措施

對不同級別的地下水保護(hù)區(qū)提出相應(yīng)的保護(hù)措施如下:

一級保護(hù)區(qū):設(shè)置隔離網(wǎng)或墻，完善取水點(diǎn)建設(shè)封閉式建筑物，由專人管理，無關(guān)人員不得入內(nèi);不得放置與取水無關(guān)的設(shè)備和物品，不得建造與引水無關(guān)的建筑物，消除一切可能導(dǎo)致地下水污染的因素;取水區(qū)內(nèi)應(yīng)嚴(yán)禁出現(xiàn)工礦建設(shè)、破壞植被以及農(nóng)業(yè)耕作等人類活動，防止地下水遭到污染，影響飲水安全;

圖10 響水巖水源地及水源地保護(hù)區(qū)劃分示意圖

二級保護(hù)區(qū):區(qū)域內(nèi)不得堆放垃圾、廢渣或污水處理管道，嚴(yán)禁大量使用農(nóng)藥、化肥等對地下水水質(zhì)產(chǎn)生變化的人類工程，嚴(yán)禁出現(xiàn)破壞水源地環(huán)境的人類活動;

三級保護(hù)區(qū):防止污染水體流入二級區(qū)，保證有充分的時間發(fā)現(xiàn)污染、采取措施消除污染源，確保水源地正常運(yùn)行。實(shí)行“取水許可證制度”，合理開采利用地下水，劃定地下水禁采區(qū)和限采區(qū)，并實(shí)行用水企業(yè)取水許可證制度。建立水量水質(zhì)監(jiān)測，定期監(jiān)測水量、水質(zhì)、水溫。

4 結(jié)語

研究區(qū)域?qū)賻r溶山區(qū)，水資源開發(fā)利用條件差，地下水是當(dāng)?shù)刂匾霓r(nóng)村和城鎮(zhèn)人畜飲水的水源。本研究應(yīng)用先進(jìn)的地下水流場和污染物質(zhì)遷移有限元模擬FEFLOW軟件系統(tǒng)，建立了研究區(qū)域流場數(shù)值模擬模型，分別模擬和預(yù)測豐、平、枯水年不同降水條件下地下水流場分布和動態(tài)變化趨勢。地下水?dāng)?shù)值模擬模型定量分析研究區(qū)域地下水流場分布和補(bǔ)給、徑流和排泄情況，為當(dāng)?shù)氐叵滤Y源評價和重要地下水水源地保護(hù)區(qū)劃分提供理論依據(jù)。另外，分析了黔西縣重點(diǎn)保護(hù)地下水水源地的水量需求和水質(zhì)要求?；谀Ｐ湍M結(jié)果，根據(jù)國家《飲用水水源保護(hù)區(qū)劃分技術(shù)規(guī)范》的要求，初步劃分了研究區(qū)域內(nèi)重要地下水水源地一級、二級和三級保護(hù)區(qū)范圍。針對水源地特點(diǎn)，提出了地下水水源地綜合保護(hù)措施，從工程技術(shù)層面和社會經(jīng)濟(jì)層面考慮，給出了具體可操作的意見和合理的地下水保護(hù)對策和建議。

[1]H.－J.G.Diersch，F(xiàn)EFLOW Reference Manual(2005)，DHI－WASY GmbH，Berlin.

[2]張立志，基于FEFLOW的大興區(qū)地下水動態(tài)模擬研究，碩士學(xué)位論文;中國地質(zhì)大學(xué);北京，2009年。

[3]趙旭，基于FEFLOW和GIS技術(shù)的咸陽市地下水?dāng)?shù)值模擬研究，碩士學(xué)位論文;西北農(nóng)林科技大學(xué);西安，2009年。

[4]畢節(jié)市水資源公報，貴州省畢節(jié)市水利局，2010年。

[5]飲用水水源保護(hù)區(qū)劃分技術(shù)規(guī)范，中國環(huán)境保部，2007年。

[6]中華人民共和國污染防治法，2008年。