袁冬梅, 齊躍明, 王俊萍, 馬于曦

(中國礦業(yè)大學(xué)資源與地球科學(xué)學(xué)院，徐州 221116)

水資源不僅是一種重要的戰(zhàn)略資源，同時也是綜合國力的重要組成部分，支撐國民經(jīng)濟的發(fā)展[1]。而地下水具有水質(zhì)良好、分布廣泛、變化穩(wěn)定以及便于利用的特點，是理想的供水水源。在中國北方大部分地區(qū)，地下水是主要的供水水源，有時甚至是唯一的生活以及工農(nóng)業(yè)生產(chǎn)供水水源[2]。而地下水庫是優(yōu)化水資源配置的一種重要手段，可以加大對降水資源的截流，提高區(qū)域水資源的利用率，還可以豐蓄枯采，調(diào)節(jié)水資源時空分配不均等[3]。

以往計算地下水可開采資源量的常用方法有水均衡法[4]、解析法[5]、數(shù)值法[6]、補償疏干法[7]等，各種方法均有其適用范圍和應(yīng)用條件，均有其對應(yīng)的缺陷。如水均衡法對于地下水資源、可開采資源量的時間、空間等方面的變化以及對未來的預(yù)測等“動態(tài)”信息不能夠及時且有效地體現(xiàn)出來。解析法假設(shè)條件過多，不適用于復(fù)雜水文地質(zhì)條件下的水量計算。數(shù)值法方法對水位、水質(zhì)、水量等觀測數(shù)據(jù)的精度要求較高，源匯項及邊界條件要求刻畫清楚。補償疏干法旱季允許借用的儲存量必須滿足旱季連續(xù)開采而不能中斷開采，且雨季得到的補給量必須把旱季借用的儲存量全部補償回來，而不能部分補償。相比于以上幾種方法，地下水庫調(diào)蓄法是近幾年來逐漸興起的一種方法，分為多年調(diào)節(jié)和年調(diào)節(jié)兩種類型。它能根據(jù)多年補給、排泄等序列數(shù)據(jù)，依據(jù)“以豐補歉”的思想，通過地下庫容對豐水年份或豐水季節(jié)多余的水進行存儲，在枯水年份或枯水季節(jié)來加以開采利用。

地下水庫的調(diào)蓄能力取決于三個因素：一是能在降水期間獲得整個水文周期內(nèi)所需水量，即系統(tǒng)接受補給的能力；二是系統(tǒng)的蓄水能力，蓄水能力是評價系統(tǒng)調(diào)蓄能力最重要的標(biāo)準(zhǔn)，只有空間足夠大才能接納和釋放足夠的水量；三是系統(tǒng)的給水能力，這里所說的給水能力不是一次給水能力，而是在天然流動條件下，在整個水文年內(nèi)能不能滿足當(dāng)?shù)赜盟枨?，是持續(xù)給水能力。相較于地表水庫調(diào)蓄法，地下水庫調(diào)蓄法具有庫容較難于確定，受降水量影響響應(yīng)慢的缺點，但同時又具有不占用地面空間、水質(zhì)好、能凈化區(qū)域水環(huán)境及多年調(diào)節(jié)性強等優(yōu)點。

地下水庫在國外常用提法是“人工補給含水層(artificial rechanrge of groundwater)”[8]或“含水層儲存與回采(aquifer storage and recovery)”。20世紀初期，美國就已經(jīng)進行大規(guī)模人工補給地下水實踐工作，從1973年開始就已經(jīng)形成了較為成熟的含水層儲存、回采技術(shù)[9-10]。20世紀后半期，沙特開始實施“人工蓄水池補給管理計劃”，修建了200多個大大小小的壩，將洪水?dāng)r蓄在間歇性干涸的河谷下覆的沖積含水層中。除了美國、沙特，以色列、利比亞、荷蘭、日本、巴西等不少國家也采取了各種措施將地表水資源暫時存儲在地下含水層中，即地下水庫。相比于國外，對于地下水庫的研究及建設(shè)、實踐，中國要稍晚一些，但是，中國在地下水庫功能的開發(fā)上取得的成就并不小。巖溶地下水庫是近20年來逐漸為人們所了解并加以研究的熱點方向。2001年，第一次喀斯特地區(qū)可持續(xù)發(fā)展國際會議[11]就提出來在平原喀斯特地區(qū)建設(shè)地下水庫，更好地開發(fā)利用地下水，以促進經(jīng)濟的發(fā)展。此后，陳洪元等[12]設(shè)計將羊皮寨地下河集雨區(qū)雨水資源通過引洪入庫工程引入馬官地下水庫，實現(xiàn)了運用馬官地下水庫對于巖溶山區(qū)雨水資源的有效開發(fā)利用。李世君等[13]遵循“豐補旱采、以豐補歉”的原則，充分利用地下水庫強調(diào)蓄功能，在北京懷柔建設(shè)應(yīng)急水源地，以保證在連續(xù)枯水年有足夠的出水量。付曉忠等[14]利用洮兒河沖洪積扇地區(qū)巨大的地下水蓄水空間，引江河水補充地下水，緩解該區(qū)水資源供給矛盾。齊躍明等[15]提出可運用補償疏干法來評價淄河源區(qū)地下水可開采資源量。

文獻分析表明，前人運用地表水庫能實現(xiàn)區(qū)域水資源的多年調(diào)節(jié)或年調(diào)節(jié)，利用地下水庫能有效儲存利用雨水資源，實現(xiàn)以豐補歉的目標(biāo)。然而，某些文獻對于形成地下水庫的水文地質(zhì)條件要求交代不足，對于地下水庫的調(diào)蓄作用更多地是停留在定性分析，而缺乏定量的可開采資源量計算。此外，以往文獻多研究孔隙型含水層地下水庫，而對于巖溶裂隙型地下水庫調(diào)蓄作用研究相對不足?；诖?，本文在現(xiàn)場水文地質(zhì)勘察的基礎(chǔ)上，詳細地交代淄河源區(qū)謝家店水源地周邊的水文地質(zhì)條件——北博山鎮(zhèn)謝家村以上的淄河源區(qū)具有一個相對封閉的邊界條件，無側(cè)向補給來源，大氣降水是地下水庫的唯一補給來源[16]。各含水層段基巖均為奧陶—寒武系灰?guī)r、白云巖、白云質(zhì)灰?guī)r等碳酸鹽巖類，巖溶發(fā)育程度強烈[17-18]，不需要建設(shè)地下?lián)跛畨伪隳苄纬捎行钏畮烊莸牡叵滤畮靃19]。因而可利用巖溶地下水庫年調(diào)節(jié)的方法，來計算謝家店水源地巖溶地下水源地可開采資源量，實現(xiàn)當(dāng)?shù)厮Y源可持續(xù)開采利用，解決特殊干旱年份的用水需求，同時，亦可為地下水庫調(diào)蓄水資源的相關(guān)理論或?qū)嵺`研究提供一個較為成功的應(yīng)用案例。研究亦將為當(dāng)?shù)卣涂蒲胁块T更好地制定科學(xué)合理的水資源開發(fā)方案及相關(guān)保護對策提供決策依據(jù)。

1 研究區(qū)概況

謝家店水源地位于淄博市博山區(qū)南部(圖1)，屬溫帶半濕潤大陸性季風(fēng)氣候，冬季寒冷，夏季炎熱，多年平均氣溫12 ℃，多年平均降水量720.6 mm(表1)。地貌類型為周邊中低山，中部為淄河河谷地形。出露地層主要有太古界泰山群、寒武系、奧陶系及第四系地層，地層整體產(chǎn)狀傾向北西，傾角5°～20°(圖2為沿A-B切面的剖面圖)。區(qū)內(nèi)地表水主要有石馬水庫及季節(jié)性河流淄河。區(qū)域內(nèi)主要有石馬斷裂、盆泉-北博山斷裂。

表1 博山區(qū)多年平均月降水量表

圖1 研究區(qū)交通位置圖Fig.1 Traffic location map of the study area

垂直高度：水平距離=1∶1.5圖2 研究區(qū)A-B地質(zhì)剖面圖Fig.2 Regional A-B geological section

據(jù)文獻[15]修改圖3 研究區(qū)水文地質(zhì)圖Fig.3 Hydrogeological map of the study area

三山子組-炒米店組裂隙巖溶水含水層，是該地區(qū)主要含水層之一。主要分布于淄河以東的中低山區(qū)及南博山、謝家店一帶，出露位置較高，節(jié)理裂隙發(fā)育，接受大氣降水入滲及地表水滲漏補給。富水性較強，單井涌水量大于1 000 m3/d。

北庵莊組灰?guī)r含水層在淄河以西及地塹帶大部分埋藏地面以下。源泉-泉河頭一帶，頂板埋深25～130.09 m，巖溶裂隙發(fā)育，溶蝕裂隙及溶洞含水豐富，也為研究區(qū)主要含水層之一。北博山-謝家店一帶，北庵莊組灰?guī)r埋藏較淺，富水性一般大于1 000 m3/d，三山子組及炒米店組白云巖、灰?guī)r含水層時，水量較大，單井涌水量一般大于3 000 m3/d。

塊狀巖類風(fēng)化裂隙水含水巖組主要分布于研究區(qū)的東南端、南端及西南端，水質(zhì)良好，為人畜可利用含水層。含水層下部完整，基巖多為花崗片麻巖，其風(fēng)化層厚度為10～30 m，含水層厚度8～25 m。地下水主要賦存于花崗片麻巖裂隙與風(fēng)化層中，屬裂隙潛水類型。地下水主要補給源為大氣降水，水位埋深2～3 m。巖層富水性較為均勻，但單井出水量較小，小于100 m3/d，水位、水量年動態(tài)變化顯著。

研究區(qū)地下水的補給幾乎全部來自于大氣降水，地下水流場形態(tài)由地表地勢和斷裂構(gòu)造控制，地下水徑流方向是由研究區(qū)的西南部、南部向淄河河谷方向徑流，淄河斷裂帶為該區(qū)地下水的強徑流帶，地下水向淄河河谷方向匯流之后繼續(xù)沿淄河地塹向北東方向徑流。

2 地下水庫調(diào)蓄法

2.1 邊界條件

淄河源區(qū)西部邊界為祿由頂—橫頂村—苗山一線的地表分水嶺；南部邊界為魯山山脈一線的地表分水嶺；東部邊界為邀兔村水力隔水邊界；北部邊界為石馬斷裂，具有一個相對封閉的邊界條件，無側(cè)向補給來源，圈定范圍面積約為300 km2。

2.2 儲水空間

考慮到巖溶地區(qū)有效儲水空間只有巖溶裂隙發(fā)育的層段即含水層段，則儲水空間為

(1)

式(1)中：V為地下水庫儲水空間，m3；F為地下水庫面積，km2；μi為第i層段含水層的給水度；Mi為第i層段含水層的厚度m。

表2 謝家店含水層段巖溶率及給水度統(tǒng)計表

將表2數(shù)據(jù)代入式(1)中得到，淄河上游地下水庫的有效儲水空間約為2.19×109m3。

2.3 地下水庫補給量

考慮到庫區(qū)邊界相對封閉隔水，因而可認為基本無側(cè)向補給來源。石馬水庫作為灌溉和景觀用小型水庫，可視作內(nèi)部地表儲存水，其降水補給和綜合排泄視作內(nèi)部水量平衡過程，趨于穩(wěn)定而不單獨考慮參與計算。根據(jù)淄河源區(qū)水循環(huán)示意圖(圖4)可知，庫區(qū)內(nèi)大氣降水是庫水的唯一補給來源，降水后形成地表徑流和地下徑流。地表徑流在流動過程中，很大一部分沿淄河河床發(fā)生滲漏，逐漸被第四系和碳酸鹽巖含水層所吸收，另一很小部分則為蒸發(fā)損耗(與晚間凝結(jié)水對應(yīng)，兩相對比，此項基本可忽略)，等到謝家店村，河流基本呈干涸狀態(tài)。

圖4 庫區(qū)水循環(huán)示意圖Fig.4 Sketch map of the water cycle in the reservoir area

淺層地下徑流主要有第四系水、碳酸鹽巖含水層水及塊狀巖類風(fēng)化裂隙水。其中碳酸鹽巖含水層水在由地勢較高處沿斷裂構(gòu)造或?qū)永怼⒐?jié)理裂隙流向淄河河谷的過程中，一部分受到隔水?dāng)鄬踊蜃杷w而溢流成泉(如謝家店泉群，泉流出后在流動過程中又逐漸滲入地下，因此，屬于內(nèi)部水量交換)。另一部分則受地形和人工開采的影響，使地下水流動速度加快。

根據(jù)實地考察、現(xiàn)場簡易試驗并結(jié)合此次勘察成果，庫區(qū)巖漿巖分布區(qū)、碳酸鹽巖裸露區(qū)和第四系松散層覆蓋的各分區(qū)降雨入滲系數(shù)值如表3所示。

表3 庫區(qū)大氣降水入滲系數(shù)分區(qū)表

逐月入庫流量值可依式(2)求得：

(2)

式(2)中：Qi為各月多年平均月入庫流量，m3/d；Pi為各月多年平均月降水，mm；α為分區(qū)降水入滲系數(shù)；F為分區(qū)面積，km2；θ為有效系數(shù)；ΔT為月份天數(shù)，取平均值30.4 d。

把表1中各月多年平均月降雨深代入式(2)可得庫區(qū)逐月入庫流量，如表4所示。

表4 庫區(qū)逐月入庫流量

2.4 用水量統(tǒng)計

淄河源區(qū)謝家店地下水庫的用水量主要包括工業(yè)需水量300 m3/d，農(nóng)業(yè)需水量8.82×104m3/d以及供水需水量1.5×104m3/d，共計10.65×104m3/d。

2.5 時歷列表法計算地下水庫調(diào)蓄

本次計算采用早蓄方案，就是在余水期，有余水就蓄，可以滿足用水需求后還有多余水量再棄水。由于該地下水庫基本屬于裸露型或局部覆蓋型巖溶地區(qū)，地下水補給滯后時間較短，6—9月為雨季，符合年調(diào)節(jié)水庫特點，因此，采用6月—次年5月作為水利年進行調(diào)節(jié)計算，取計算時段ΔT(取平均值30.4 d)為一個月，水庫在調(diào)節(jié)周期內(nèi)只有一次連續(xù)蓄水、供水的情況(表5)。

表5 列表法年調(diào)節(jié)計算

由表4可知，10月—次年5月為虧水期，8個月總虧水量為470 273.77 (m3/d)·月[說明：當(dāng)ΔT為常數(shù)時，常用“流量·時間”為單位表示水量，例如(m3/d)·月，適應(yīng)于地下水庫或(m3/s)·月，適應(yīng)于地表水庫]。1(m3/d)·月=30 m3。同理1(m3/s)·月=2 626 560 m3，6—9月為余水期，4個月總余水量為629 830.06 (m3/d)·月。余水期多余的水量遠超過虧水期所缺少的水量。所以余水期只需要蓄470 273.77 (m3/d)·月的水量，即可滿足本年用水需要。根據(jù)水庫的運行方式可得出水庫各月的蓄水量變化情況及水庫棄水。

計算結(jié)果可根據(jù)水量平衡關(guān)系進行驗算，年來水流量應(yīng)該等于年用水量與棄水量之和，即：1 278 000+159 556.29=1 437 556.29 [(m3/d)·月]；年余水量應(yīng)該等于虧水量與棄水量之和，即：470 273.77+159 556.29=629 830.06 [(m3/d)·月]。經(jīng)檢驗，計算結(jié)果無誤。

2.6 可開采資源量計算結(jié)果

由上述計算結(jié)果可知，余水量大于虧水量，即余水期的補給可以滿足虧水期的需水量，甚至還有足夠盈余，也就是說二者之差棄水量，成為一部分可以加以利用的水量?？砷_采資源量除設(shè)計供水需水量1.5×104m3/d，還可將棄水量平均至全年的部份也包括進去，從而擴大了可開采資源量。按照此思路，由式(3)可求得淄河源區(qū)謝家店地下水水源地的可開采資源量為2.74×104m3/d。顯然，這個是理想情況下最大的可開采資源量。

Q最大可開=供水需水量+月平均棄水量=

2.74×104m3/d

(3)

然而，考慮到地下水庫調(diào)節(jié)時間較長，同時也要考慮到區(qū)域水環(huán)境水資源更新需要及氣候波動的情況，因此，為保險起見，需在上述計算的基礎(chǔ)上，乘以一個安全系數(shù)，一般取0.8，來作為源區(qū)能可持續(xù)開采的水資源量，即有：

Q可開=0.8Q最大可開≈2.19×104m3/d

(4)

3 結(jié)果與討論

由于淄河源區(qū)降雨徑流具有年周期的特點，為簡單起見，采用地下水庫年調(diào)節(jié)的方法來進行可開采資源量的計算。地下庫容通過匯水面積、含水層分布和有效空隙率等來綜合確定。

依據(jù)鉆孔巖心含水層分層統(tǒng)計資料結(jié)合飽水試驗、給水試驗等，得出研究區(qū)巖溶地下空間容納體積約為2.19×109m3。根據(jù)資料統(tǒng)計謝家店地下水庫用水量為10.65×104m3/d。依據(jù)“以豐補歉”的思想，將雨季開采外的余水量充分利用，供虧水期維持正常開采，運用地下水庫年調(diào)節(jié)方法，計算得出可開采資源量為2.19×104m3/d?？紤]到地下水庫調(diào)節(jié)時間較長和區(qū)域水環(huán)境水資源更新需要及氣候波動的情況，取0.8的安全系數(shù)，以保證能滿足源區(qū)可持續(xù)開采水資源量，實現(xiàn)采補平衡。

綜合表明，地下水庫調(diào)蓄法充分利用了巖溶地下儲水空間的調(diào)蓄功能，能夠克服其他水資源量評價方法的弊端，且具有很好的地域適應(yīng)特點，不僅適用于巖溶發(fā)育良好地段，且同樣可用于松散介質(zhì)和裂隙介質(zhì)的地下水庫。

4 結(jié)論

(1)依據(jù)“以豐補歉、采補平衡”的思想，來合理規(guī)劃利用巖溶地下水庫是解決區(qū)域水資源不足的一種非常好的途徑。能有效緩解水資源供需矛盾，提高城市供水保障程度，對維持社會經(jīng)濟的可持續(xù)發(fā)展具有十分重要的現(xiàn)實和戰(zhàn)略意義。

(2)運用地下水庫調(diào)蓄法時，要使之能攔截地下徑流，增加地下水的儲蓄量，提高水資源利用率，且受降雨量影響變化慢，多年調(diào)節(jié)性強。因此巖溶山區(qū)可充分利用地下空間，修建以含水層為調(diào)蓄空間的蓄水實體。

(3)根據(jù)淄河源區(qū)降水情況，可適當(dāng)修建人工回灌補給工程，以便將豐水年份、豐水期來不及利用的雨洪水資源儲存于該地下水庫中，補償巖溶地下含水層旱季水分虧缺，從而實現(xiàn)水資源可持續(xù)優(yōu)化利用的目的。