祁福利，李永利，張峰龍，魯守剛，王宏存

(1.中國地質(zhì)大學(武漢)，湖北武漢 430047;2.黑龍江省九〇四水文地質(zhì)工程地質(zhì)勘察院，黑龍江哈爾濱 150027;3.黑龍江省地下水文學省級領(lǐng)軍人才梯隊，黑龍江哈爾濱 150027)

近年來由于降水量普遍偏小，作為七臺河市三區(qū)一縣城區(qū)唯一供水水源地的桃山水庫蓄水水位、蓄水量持續(xù)下降。為確保城市居民生產(chǎn)生活用水問題，在桃山水庫上游倭肯河河谷發(fā)現(xiàn)一處地下水類型為白堊系砂巖裂隙水的大型水源地。為確保該水源地的可持續(xù)開發(fā)利用和極端干旱年份的應急供水需求，準確可靠的確定該大型水源地的可開采量并對地下水水位的變化情況進行預測預報十分必要。因此，建立可靠的地下水流數(shù)值模擬模型對當?shù)氐叵滤Y源及其環(huán)境問題進行綜合評價，是十分必要和迫切的，具有重要的理論意義和實用價值。

目前，國內(nèi)外專家學者開發(fā)了許多功能多樣的地下水系統(tǒng)數(shù)值模擬軟件，以其模塊化、可視化、交互性、求解方法多樣化等特點得到廣泛的使用，其中以MODFLOW最為突出。不同模型各有優(yōu)勢，而MODFLOW的數(shù)據(jù)結(jié)構(gòu)更易于實現(xiàn)與GIS的整合［1］;可從多角度校準模型，目前已經(jīng)通過補給、排泄和水均衡的評估對模型進行了合理的校準，并達到廣泛應用和推廣［2］;應用范圍廣，從傳統(tǒng)的模擬已經(jīng)推廣到區(qū)域地下水系統(tǒng)的水化學和同位素變化情況模擬，并結(jié)合實例進一步強調(diào)了做細致準確的穩(wěn)態(tài)流分析對于瞬時流分析的重要性［3］。

近幾年，我國學者在開發(fā)新模塊上也做了大量的工作，主要針對以往模擬模型中暴露出的種種問題，在新理論和新方法的指導下，不斷創(chuàng)新，開發(fā)新的模塊，以提高模擬結(jié)果的可靠性和軟件的實用性和快捷性。

本文以七臺河市特大型應急水源地為例采用MODFLOW對區(qū)內(nèi)的基巖裂隙水進行數(shù)值模擬探討，進而為七臺河市特大型應急水源地提供合理的地下水可開采量以及極端干旱年份時的地下水應急可開采量，并對不同地下水開采量的地下水位變化情況進行預測預報，進而實現(xiàn)地下水資源的合理利用與當?shù)厣鷳B(tài)環(huán)境的可持續(xù)發(fā)展。

1 地質(zhì)概況

七臺河市屬大陸性氣候，多年平均氣溫為4.8℃，多年平均降水量為553.1 mm/a，主要降水月份集中在7—9月。

1.1 地質(zhì)特征

第四系廣泛分布于倭肯河河谷平原內(nèi)，包括上更新統(tǒng)顧?quán)l(xiāng)屯組，巖性上部黑灰色、黃灰色粉質(zhì)黏土，下部含礫細—粗砂，砂礫石為主，全組地層厚度10～28 m;全新統(tǒng)溫泉河組，巖性上部灰黃、黃色粉質(zhì)黏土，含少量砂、礫，厚度1～4 m;下部為淺黃色含礫中細沙，含礫中粗砂，砂礫石互層，厚度小于10 m?，F(xiàn)代河床沖積層，厚度小于15 m，巖性下部為青灰色、灰黃色含礫中粗砂、黑白雜色砂礫石互層為主，上部為灰黃色含砂粉質(zhì)黏土。

在區(qū)內(nèi)地表出露及河谷平原下伏被揭露的地層主要是白堊系下統(tǒng)，自老到新分為城子河組、猴石溝組。城子河組下段為海陸交互相含煤地層，中段為灰黃色中、細砂巖、粉砂巖與多層煤交互組成，上段為粉細、中砂巖和多層煤。底部在多數(shù)情況下都以或厚或薄的礫巖、含礫砂巖、中粗粒砂巖為主。猴石溝組巖性主要由黃、灰綠雜色凝灰質(zhì)粉砂巖、中粗砂巖、中粗礫砂巖，泥質(zhì)粉砂巖夾煤線等組成;中段夾流紋質(zhì)凝灰?guī)r，底部為底礫巖。全組明顯分三段，構(gòu)成一個從粗至細的沉積韻律，為河湖相沉積。

由于中生代以來的構(gòu)造運動對沉積地層進行了多次的水平擠壓與拉張作用，而且還有不同幅度、間歇性的垂向上的斷陷與抬升作用，形成背斜軸部斷裂構(gòu)造帶。根據(jù)物探結(jié)果，沿倭肯河河谷解譯出10條斷裂破碎帶，為地下水的賦存提供空間，為地下水的運移提供有利通道(圖1)。

1.2 水文地質(zhì)特征

圖1 斷裂帶分布示意圖Fig．1 Distribution of fault zones

研究區(qū)地貌類型主要為泥沙質(zhì)河谷平原、砂巖丘陵和玄武巖熔巖臺地，賦存地下水類型為第四系松散巖類孔隙水、白堊系砂巖裂隙水和新近系玄武巖孔洞裂隙水三種類型，水文地質(zhì)界線與地貌界線、巖性界線基本統(tǒng)一，值得注意的是在河谷平原地區(qū)存在雙層含水層結(jié)構(gòu)，上部為第四系松散巖類孔隙水含水巖組，下部為白堊系砂巖裂隙水含水巖組。

研究區(qū)地下水主要以大氣降水入滲補給為主，直接補給裸露的第四系松散巖類孔隙水含水巖組、白堊系砂巖風化裂隙水含水巖組和新近系玄武巖孔洞含水巖組，埋藏型的白堊系砂巖構(gòu)造裂隙水不直接接受大氣降水入滲補給，依靠上部含水層和風化帶沿垂向構(gòu)造裂隙入滲補給。

進入到各含水巖組的地下水，總體上受地表分水嶺控制，由分水嶺向邊界徑流，相鄰含水巖組間存在越層徑流，其中在地表分水嶺及其兩側(cè)地下水由上伏含水巖組向下覆含水巖組徑流。

研究區(qū)內(nèi)地下水的排泄方式有，潛水蒸發(fā)排泄，各含水巖組的人工開采，局部邊界側(cè)向徑流排泄，被河谷切割處各含水巖組地下水向地表水體的排泄等。

(1)第四系松散巖類孔隙水

根據(jù)施工的鉆孔抽水試驗資料，可將第四系松散巖類孔隙水分為兩個富水性分區(qū):水量較豐富區(qū)，涌水量100～1 000 m3/d;水量較貧乏區(qū)，涌水量＜100 m3/d。

(2)新近系玄武巖孔洞裂隙水

主要分布在西北部熔巖臺地區(qū)，船底山組玄武巖呈致密塊狀、隱晶質(zhì)結(jié)構(gòu)，具氣孔狀或杏仁狀構(gòu)造，風化裂隙和柱狀節(jié)理較發(fā)育，含水不富水，透水性極好，在局部地段的低洼處地下水沿下伏基巖頂面以泉的形式流出地表，泉水流量較小，一般均小于10 m3/d。在有隔水層與阻水巖層的接觸部分可賦存且水量較大。

(3)白堊系砂巖裂隙水

主要含水巖組為白堊系猴石溝組(K1～2h)、城子河組(K1c)砂巖、砂礫巖。因燕山期巖漿作用、壓應力擠壓成巖較好，孔隙不發(fā)育;斷裂構(gòu)造與裂隙發(fā)育形成破碎帶，因而賦存基巖裂隙水，其中破碎帶中的裂隙帶總體呈層狀分布，大概為2～3層含水帶，累計厚度可達30 m，其間夾有泥巖或砂質(zhì)泥巖而具有承壓性;斷裂帶呈北東向展布，發(fā)育深度一般在170 m。

結(jié)合鉆孔資料，倭肯河干流河谷平原地區(qū)造裂隙水非常發(fā)育，水量豐富區(qū)，單井涌水量1 000～3 000 m3/d(圖2)。

圖2 水文地質(zhì)剖面圖Fig．2 Section map of hydrogeology

2 七臺河市特大水源地數(shù)值模擬

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

本次模擬區(qū)將桃山水庫也包含在內(nèi)，模擬面積共計120 km2，模擬區(qū)地貌類型主要為第四系沖洪積河谷平原(Ⅰ區(qū))、砂巖丘陵(Ⅱ區(qū))和玄武巖熔巖臺地(Ⅲ區(qū))(圖2)，各分區(qū)基本信息詳見表1。

表1 滲透系數(shù)和給水度分區(qū)初值表Table 1 Initial-value table of zoning of permeability coefficient and specific yield

圖3 模擬區(qū)地下水系統(tǒng)劃分圖Fig．3 Division of groundwater system in the model area

(1)含水層概化

第四系松散巖類孔隙水分布在倭肯河及其支流的河谷地帶;新近系玄武巖風化裂隙水分布于模擬區(qū)的東北部;白堊系砂巖風化裂隙水分布于河谷兩側(cè)。根據(jù)物探資料，模擬區(qū)內(nèi)倭肯河第四系下伏白堊系中張性構(gòu)造斷裂發(fā)育，形成了沿河谷呈北東向分布的構(gòu)造裂隙富水帶。這一賦水帶從上至下呈現(xiàn)出3個裂隙發(fā)育富水帶，第一層破碎帶為風化裂隙及構(gòu)造裂隙非常發(fā)育層位，埋深40～54.3 m，厚度4.7～10.2 m。第二層和第三層為構(gòu)造裂隙發(fā)育層位，其中第二層埋深122～130 m，厚度3.7～11.6 m;第三層埋深147～162.3 m，厚度3.8～16.5 m。3個富水層位累計厚度約14.8～33.0 m。模擬區(qū)白堊系中除了河谷地帶的張性構(gòu)造斷裂發(fā)育帶之外的地區(qū)富水性一般。因此本次模擬將每個構(gòu)造裂隙發(fā)育的富水帶概化為滲透性良好的孔隙介質(zhì)處理。對于和每個構(gòu)造裂隙發(fā)育層位同一層的模擬區(qū)內(nèi)的其他地區(qū)概化為透水性較差的孔隙介質(zhì)處理。

(2)邊界條件概化

在模擬區(qū)，各河流流入流出口為第四系地下水的流入流出邊界，可以作為流量邊界，全區(qū)邊界上白堊系流入流出邊界，同樣概化為流量邊界。

模擬區(qū)東北側(cè)邊界定義為給定流量邊界，地下水均由外面流向模擬區(qū);模擬區(qū)西南側(cè)邊界定義為給定水頭邊界，地下水均由該處流入水庫。將東南、西北側(cè)邊界定義為給定流量邊界，方向均為流進模擬區(qū)。倭肯河常年有水，與孔隙潛水水力聯(lián)系密切，多數(shù)地段成為潛水的排泄邊界，可概化為已知水頭的一類邊界。取潛水面為模型的上邊界。地下水通過該邊界接受河流、大氣降水、田間灌溉回歸水等入滲補給，并以蒸發(fā)的方式在該邊界上產(chǎn)生排泄。將底邊界確定為埋深180 m左右的位置，由于該深度處以下地層巖性為泥巖，裂隙發(fā)育不良，將該邊界定義為零通量邊界(圖4)。

圖4 模擬區(qū)水文地質(zhì)概念模型圖Fig．4 Conceptual model of hydrogeology in the model area

2.2 地下水流數(shù)學模型的建立

根據(jù)建立的水文地質(zhì)概念模型，將計算區(qū)地下水流模擬的數(shù)值模型為非均質(zhì)各向同性的孔隙潛水、裂隙承壓水由越流聯(lián)系的非穩(wěn)定流混合模型，可以分別建立四層地下水的數(shù)學模型，用越流量將其耦合起來。由此所建立的模型可以看作是準三維的模型。

計算區(qū)面積120 km2，采用GMS進行100×100自動矩形剖分，剖分成7層，共剖分活動單元格22 246個，每個單元格200 m×200 m，面積0.04 km2。

模型識別期自2013年4月1日—9月30日，以一個月為一個應力期，每個應力期包括6個時間步長。根據(jù)統(tǒng)測資料確定的2013年4月1日含水層的初始流場(圖5)。

圖5 模擬區(qū)地下水等水位線圖(2013年4月)Fig．5 Contour of groundwater table in the model area(April 2013)

2.3 模型的識別與檢驗

2.3.1 模型的識別

選擇2013年4月—9月(共6個時段)進行模型的識別，以一個月為一個時間段，每個時間段包括6個時間步長。該時段為豐水期，源匯項較多，地下水均衡要素包括降水入滲補給、地表水灌溉滲漏補給、側(cè)向流入補給、河流滲漏補給、潛水蒸發(fā)排泄等，此期間地下水流場變化快，且變化幅度較大，模型識別后的流場特征能較好地反映出含水層結(jié)構(gòu)、水文地質(zhì)參數(shù)和含水層邊界性質(zhì)的變化。對計算出的地下水水位與實測水位擬合誤差進行統(tǒng)計，計算水位與實測水位在河谷地區(qū)的整體擬合程度良好。

2.3.2 模型的檢驗

本次采用多年平均狀態(tài)下的地下水源匯項，輸入模型運行20年，分析地下水流場的變化。分別選取上、中、下游段的代表性鉆孔分析;而以運行20年后地下水流場作為模型驗證的末期流場，與2013年同期水位進行對比，在20年的模型驗證期內(nèi)，地下水位呈周期性變化，總體保持不變，說明含水層結(jié)構(gòu)、邊界條件的概化、水文地質(zhì)參數(shù)的選取是合理的，所建立的數(shù)學模型能較真實的刻畫出模擬區(qū)地下水系統(tǒng)特征，可以利用該模型對水位的變化進行預測預報。

2.4 地下水位的預報

2.4.1 預報方案

本次對模擬區(qū)地下水位的預報根據(jù)水源地的使用方式，設(shè)計2個方案:

方案1:作為常規(guī)水源地開采，水源地可開采資源量5.2×104m3/d，模型運行20年;

方案2:作為應急供水水源地，(只在極枯年(P=97%)，降雨量為303.87 mm，桃山水庫僅死庫容1 900×104m3)確定某一應急取水量，在該取水量連續(xù)運行1年，地下水漏斗中心水位下降50 m。

2.4.2 預報結(jié)果分析

為了分析應急水源地在該開采條件下地下水流場的變化，將水源地所在河段分為上、中、下游三段分析。

2.4.2.1 方案1預報結(jié)果分析

此方案下，按照資源評價結(jié)果確定開采量5.2×104m3/d進行開采，將有關(guān)數(shù)據(jù)輸入經(jīng)過識別和驗證的地下水流數(shù)值模型進行水位預報，預測到2014年、2018年、2023年、2028年和2033年地下水位的變化情況。模型結(jié)果顯示在開采5.2×104m3/d開采基本沒有出現(xiàn)含水層疏干現(xiàn)象，因此采用5.2×104m3/d可作為可持續(xù)開采量的參考依據(jù)。

為了分析新建水源地在該開采條件下地下水流場的變化，將水源地所在河段分為上、中、下游三段分析，上游段選擇代表鉆孔ZK100，中游段選擇代表鉆孔ZK29，下游段選擇代表性鉆孔ZK2。

(1)上游段鉆孔ZK100:隨著水源地地下水的開采，在前期下降速度較快，第一年水位下降0.92 m，到第5年水位下降1.49 m，第10年水位下降1.67 m，15年水位下降1.74 m，20年水位下降1.77 m。從水位下降量隨時間的變化關(guān)系可以看出，水位在前期下降迅速，第10年以后基本保持不變，說明由于水源地的開采，地下水的降落漏斗基本穩(wěn)定，不再繼續(xù)下降，該處的水位在該開采量下的水位下降量維持在1.80 m(圖6)。

(2)中游段鉆孔ZK29:隨著水源地地下水的開采，同樣在前期下降較快，第一年水位下降1.84 m，到第5年水位下降3.26 m，第10年水位下降3.64 m，15年水位下降3.73 m，20年水位下降3.77 m。從水位下降量隨時間的變化關(guān)系可以看出，水位在前期下降速度較快，到了第10年以后基本保持不變，說明由于水源地的開采，地下水的降落漏斗基本穩(wěn)定，不再繼續(xù)下降，該處的水位在該開采量下的水位下降量維持在3.80 m(圖7)。

圖6 預測ZK100地下水位歷時曲線Fig．6 Forecasting of duration curve of groundwater table in ZK100

圖7 預測ZK29地下水位歷時曲線Fig．7 Forecasting of duration curve of groundwater table in ZK29

(3)下游段鉆孔ZK2:隨著水源地地下水的開采，地下水位總體變化量不大，主要原因是由于受到了下游桃山水庫的影響。第一年水位下降1.00 m，到第5年水位下降1.35 m，第10年水位下降1.37 m，15年水位下降1.39 m，20年水位下降1.40 m。從水位下降量隨時間的變化關(guān)系可以看出，水位下降量很小，到了第5年以后基本保持不變，說明由于水源地的開采，地下水的降落漏斗基本穩(wěn)定，不再繼續(xù)下降，該處的水位在該開采量下的水位下降量維持在1.50 m(圖8)。

圖8 預測ZK2地下水位歷時曲線Fig．8 Forecasting of duration curve of groundwater table in ZK2

此方案下進行地下水開發(fā)利用，含水層沒有出現(xiàn)疏干現(xiàn)象，河谷地區(qū)地下水位全部出現(xiàn)下降，但是下降深度不大。水源地地下水位下降量中游大，上游和下游逐漸減少?？拷疑剿畮焯幍叵滤换静蛔?圖9)。

2.4.2.2 方案2預報結(jié)果分析

根據(jù)方案2，若該水源地作為應急供水水源地，一般情況下不會啟用，只有極端干旱年份——水庫水位降至171.57 m時(死庫容對應水位)，由于水庫庫底崎嶇不平，水面難以連續(xù)，呈現(xiàn)多個分散的“水泡子”現(xiàn)象，但仍有1 900×104m3死庫容難以抽取。這時啟用應急供水水源地，按理論推算，如果地質(zhì)條件可行的情況下，1 900×104m3死庫容完全補給地下水，按照15×104m3/d的供水計算，死庫容仍可維持一年，因此利用地下水數(shù)值模型根據(jù)方案2研究水庫死庫容水與裂隙地下水的水力聯(lián)系極為重要。

為了分析新建水源地在該開采條件下地下水流場的變化，確定降落漏斗的中心位置，上游段選擇代表鉆孔ZK100，中游段選擇代表鉆孔ZK29，下游段選擇代表性鉆孔ZK2地下水位隨時間的變化趨勢。

(1)降落漏斗中心:通過調(diào)整開采井的開采量，最終確定當?shù)叵滤拈_采量達到18.2×104m3/d時，地下水的降落漏斗中心位置位于中游鉆孔ZK50附近，水位下降量49.40 m(圖10～圖11)。

圖9 地下水水位預測圖(a)1年后、(b)5年后、(c)10年后、(d)20年后Fig．7 Forecasting of groundwater table(a)1 year later(b)5 years later(c)10 years later(d)20 years later

圖10 預測ZK50地下水位歷時曲線Fig．10 Forecasting of duration curve of groundwater table in ZK50

圖11 2014年地下水水位預測圖(12個月后)Fig．11 Forecasting of groundwater table in 2014(12 months later)

(2)上游段鉆孔ZK100:隨著水源地地下水的開采，在前8個月水位基本呈均勻下降，8個月以后下降趨勢增大4倍，3個月后水位下降3.3 m，6個月后水位下降到6.45 m，9個月后水位下降18.45 m，1年后水位下降46.11 m。從圖12可以看出，前8個月水位降深為0.82 m/月，而之后的5個月水位降深為9.15 m/月，表明前8個月水庫補給充分，大量抽取地下水得到水庫的及時補給，當連續(xù)抽取8個月，水庫死庫容逐漸消耗，水頭差減小，補給能力衰減，逐漸開始抽取第一層破碎帶中的凈儲量，導致水位急劇下降，模型繼續(xù)運行至1年，水位急劇下降約50 m，達到第一層破碎帶水的底板，第一層破碎帶疏干。因此，采用18.2×104m3/d進行開采，可以保障應急1年的用水量(圖12)。

(3)中游段鉆孔ZK29:隨著水源地地下水的開采，在前期6個月水位基本呈均勻下降，到6個月后下降趨勢增大4倍，3個月后水位下降2.60 m，6個月后水位下降到15.27 m，9個月后水位下降至21.06 m，1年后水位下降40.67 m。分析6個月后地下水位下降深度增大原因同ZK100。從水位下降量隨時間的變化關(guān)系可以看出，水位下降量基本與ZK100相同。因此，采用18.2×104m3/d進行開采，可以保障應急1年的用水量(圖13)。

圖12 預測ZK100地下水位歷時曲線Fig．12 Forecasting of duration curve of groundwater table in ZK100

圖13 預測ZK29地下水位歷時曲線Fig．13 Forecasting of duration curve of groundwater table in ZK29

(4)下游段鉆孔ZK2:隨著水源地地下水的開采，3個月后水位下降1.39 m，6個月后水位下降1.76 m，9個月后水位下降至2.6 m，1年后水位下降15.49 m。從水位下降量隨時間的變化關(guān)系可以看出，水位在前期較慢，8個月后水位下降明顯加快，但水位下降量較中游和上游要小很多，主要原因是下游靠近桃山水庫，在水位下降后更容易接受桃山水庫的地表水補給。1年后桃山水庫補給量平均為12.20×104m3/d(圖14)。

圖14 預測ZK2地下水位歷時曲線Fig．14 Forecasting of duration curve of groundwater table in ZK2

此方案下進行地下水開發(fā)利用，上覆潛水含水層全部疏干，地下水位在后期下降明顯，水源地地下水降落漏斗在中游地段，若應急開采1年后桃山水庫補給地下水的平均量為12.20×104m3/d。若采用該方案進行開采，對當?shù)氐叵滤Y源和桃山水庫的庫容具有很大的影響，因此建議只在應急時開采，應急開采時間不宜大于12個月。

該方案情況下，9個月后桃山水庫補給地下水量平均為8.5×104m3/d，一年后桃山水庫補給地下水量平均為12.2×104m3/d;前9個月累計水庫補給量1 213.6×104m3，理論上前12個月累計水庫補給量2 736.2×104m3，由于水庫初始庫容 1 930 ×104m3，在11月末水庫死庫容疏干，完全失去補給能力。因此，采用當?shù)氐叵滤Y源和桃山水庫的庫容聯(lián)合運營模式開采18.2×104m3/d維持一年是有保障的。

3 結(jié)論

本次研究在深入分析七臺河倭肯河河谷桃山水庫上游區(qū)地質(zhì)及水文地質(zhì)的基礎(chǔ)上，利用MODFLOW軟件建立了該研究區(qū)的地下水流模型，通過識別和驗證后，提出兩種不同方案對未來長期開采及應急開采條件下，地下水的變化趨勢進行了預測。得出以下結(jié)論:

(1)建立了研究區(qū)多層富水帶結(jié)構(gòu)的地下水流數(shù)值模擬模型，對研究區(qū)白堊系砂巖構(gòu)造裂隙水、第四系松散巖類孔隙潛水進行了數(shù)值模擬，利用實測資料對模型進行識別和驗證，識別和驗證效果良好，建立的水流模型達到了較高的精度，因此可以作為地下水流預報模型。

(2)研究區(qū)的地下水資源主要富集在白堊系砂巖構(gòu)造破碎帶中，構(gòu)造裂隙富水帶與下游桃山水庫具有強烈的水力聯(lián)系，在上游地下水集中開采的驅(qū)動下，水庫能及時響應補給地下水。應急開采初期水位下降趨勢相對較緩，而在大量集中開采的條件下后期水位下降迅速。反映出研究區(qū)開采地下水初期來自構(gòu)造裂隙富水帶賦存的水量，隨著開采時間的增加，地下水開始接受下游水庫的補給和上游激發(fā)條件下的側(cè)向補給量。而且應急開采量越大，接受水庫補給發(fā)生的時間越早，補給量越大。

(3)根據(jù)不同預測方案結(jié)果得出了長期可持續(xù)利用前提下的地下水開采量及不同應急開采量下的應急開采時間。

(4)從該地區(qū)地下水的長期可持續(xù)利用要求出發(fā)，該區(qū)日開采量以不超過5.2×104m3為宜。在該開采條件下，地下水位的下降深度較小，全區(qū)水位降深均維持在4.0 m以內(nèi)，對研究區(qū)及周邊環(huán)境不會造成不良影響。

(5)研究區(qū)作為應急水源地時應采用陶山水庫與地下水聯(lián)合運營模式，全力開采應急水源地，必要時可以動用桃山水庫的死庫容，開采量按18.2×104m3/d算保障開采一年。

研究區(qū)第四系松散巖類孔隙水和白堊系砂巖裂隙水具備較大的地下水開發(fā)利用潛力，地下水開發(fā)利用前景良好，適宜作為應急供水水源地。研究區(qū)地下水資源為正均衡，但屬于氣象型水源地，水源地補給來源主要為大氣降水，雖可以短時間內(nèi)大規(guī)模開采但是不適宜長時間大規(guī)模開采，以豐補歉、多年調(diào)節(jié)，方可作為應急水源地用以保障完達山七臺河地區(qū)供水安全。

［1］ JUANC S， KOLMK E. Conceptualization，characterization and numerical modeling of the Jackson Hole alluvial aquifer using ARC/INFO and MODFLOW［J］.Engineering Geology，1996，42:119-137.

［2］ OLSTHOORNTN.A comparative review of analytic and finite difference models used at the Amsterdam Water Supply［J］.Journal of Hydrology，1999，226:139-143.

［3］ HARRINGTONG A，WALKERG R.A compartmental mixing-cell approach for the quantitative assessment of groundwater dynamics in the Otway Basin［J］.Journal of Hydrology，1999，214:49-63.

［4］ 盧文喜.地下水運動數(shù)值模擬過程中邊界條件問題探討［J］.水利學報，2003(3):33-36.［ LU W X.Discussion of the boundary conditions in numerical simulation of groundwater movement［J］.Journal of Hydraulic Engineering，2003(3):33-36.(in Chinese)］

［5］ 寇文杰.格爾木河流域地下水數(shù)值模擬［J］.水文地質(zhì)程地質(zhì).2013.40(1):34-40.［KOU W J.Numerical simulation of the groundwater system in the Geermu river basin［J］.Hydrogeology ＆ Engineering Geology，2013，40(1):34-40.(in Chinese)］

［6］ 楊青春，盧文喜，馬洪云.Visual Modflow在吉林省西部地下水數(shù)值模擬中的應用［J］.水文地質(zhì)工程地質(zhì)，2005，32(3):67-69［YANG Q C，LU W X，MA H Y.Application of Visual Modflow in Western Jilin Province numerical simulation of groundwater［J］.Hydrogeology ＆ Engineering Geology，2005，32(3):67-69.(in Chinese)］

［7］ 魏亞強，喬小娟，李國敏.MODFLOW不同算法及參數(shù)設(shè)定對計算精度的影響［J］.水文地質(zhì)工程地質(zhì)，2015，42(1):14-21. ［WEI Y Q，QIAO X J，LI G M.Influence on calculation accuracy caused by different algorithms and parameters in MODFLOW［J］.Hydrogeology ＆ Engineering Geology，2015，42(1):14-21.(in Chinese)］