董殿偉　馬如林　李玉虎　林佳旭

摘 要：由于地下水的大量開采及永定河補水的減少，玉泉山泉水自1975年斷流。隨著南水北調(diào)工程的實施，在西山地區(qū)進行回灌和壓采區(qū)域地下水成為可能。利用GMS數(shù)值模型，以河道補給、人工回灌、關(guān)井壓采等手段設(shè)計4種調(diào)蓄方案，分析玉泉山泉復(fù)涌可行性及復(fù)涌后對地質(zhì)環(huán)境的影響。研究表明，方案四可優(yōu)先作為玉泉山泉復(fù)涌的調(diào)蓄方案，該方案主要利用在永定河引水渠和南旱河沿岸進行大口井回灌，回灌量約為8000萬m3·a-1，同時，巖溶水集中水廠停采，其余開采量保持現(xiàn)狀不變。在此方案下，區(qū)域水位平均每年上升0.84 m，預(yù)計到2050年末，玉泉山泉水位51.9 m，達到出流水位。

關(guān)鍵詞：玉泉山泉;泉水恢復(fù);巖溶水;地下水調(diào)蓄方案;GMS;回灌;數(shù)值模擬

Predictive analysis of conditions and schemes for Yuquan Mountain spring water recovery in western suburbs of Beijing

DONG Dianwei1， MA Rulin2， LI? Yuhu1， LIN Jiaxu1

（1. Beijing Institute of Geological & Prospecting Engineering， Beijing 100048， China;

2. Beijing Chaobai River Management Office， Beijing 101300， China）

Abstract： Due to massive groundwater exploitation and reduction of? Yongding River water replenishment， Yuquan Mountain's spring water has been dry since 1975. With the implementation of the South-to-North Water Diversion Project， it is possible to recharge and press regional groundwater in the Xishan area. Using the GMS numerical model， four kinds of regulation and storage schemes are designed by means of river recharge， artificial recharge， shut-in and pressure production， to analyze the feasibility of? Yuquan Mountain spring recovery and its impact on the geological environment. Studies show that Option 4 can be prioritized as the regulation and storage plan for Yuquan Mountain spring recovery. The plan mainly recharges the large wells along the diversion canal of? Yongding River and Nanhan River， with a recharge volume of about 80 million m3.a-1. At the same time，? concentrated water plants stopped exploring karst water， and the rest of the water exploitation volume remained unchanged. Under this plan， the regional water level will rise by 0.84 m per year on average. It is estimated that by the end of 2050， the water level of the Yuquan Mountain spring will be 51.9 m， reaching the outflow level.

Keywords： Yuquan Mountain spring; recovery of the spring water; karst water; groundwater regulation and storage schemes; GMS; recharge; numerical simulation

玉泉山泉有“天下第一泉”之稱，自元代便專供皇宮使用，并作為源頭，補給北長河、昆明湖、高水湖及金河等河道，河道的生態(tài)環(huán)境優(yōu)美宜人，著名的“京西稻”就是由玉泉山泉水澆灌的。由于地下水的大量開采，玉泉山泉自1975年斷流，不復(fù)當年“水清而碧，澄潔似玉”之盛景。面對此種現(xiàn)象，2005年李裕宏提出“涵養(yǎng)京西水源，恢復(fù)玉泉山泉”的設(shè)想，但由于當時北京西山地區(qū)為北京市重要的巖溶水開采基地，肩負著北京城區(qū)的供水任務(wù)，恢復(fù)玉泉山泉難以實現(xiàn)。而隨著南水北調(diào)工程的實施，在西山地區(qū)進行回灌和壓采區(qū)域地下水成為可能，玉泉山泉有望實現(xiàn)復(fù)涌（王莉蛟等，2016）。因此，本文根據(jù)玉泉山泉水的出流特征，分析泉水補給來源和徑流通道，利用地下水數(shù)值模擬模型，探討泉水復(fù)涌的途徑和方案。

1 泉水出流原因

斷裂構(gòu)造發(fā)育是玉泉山泉水的出露的重要條件。目前研究表明，玉泉山區(qū)域巖溶水主要有3個補給路徑，第一是軍莊—雁翅地區(qū)巖溶裂隙水接受大氣降水補給后，沿永定河張性斷裂，經(jīng)石景山區(qū)向北東轉(zhuǎn)折流向玉泉山，補給玉泉山泉水;第二是潭柘寺地區(qū)降水入滲補給巖溶裂隙水向南東方向匯流，沿八寶山斷裂及北側(cè)裂隙-巖溶發(fā)育帶流向四季青方向，補給玉泉山泉水;第三是永定河河水入滲補給地下水，滲流方向首先沿河床向南運移，在古城一帶折向東北，朝向四季青方向運移，補給玉泉山泉水（圖1）。玉泉山位于巖溶水系統(tǒng)的排泄區(qū)，又呈現(xiàn)“天窗”即奧陶系灰?guī)r裸露，具有西山奧陶系巖溶水進入平原后獲得釋放的自然條件，出露成泉（北京市水文地質(zhì)工程地質(zhì)大隊，2017）（圖2），由于奧陶系直接出露地表，與第四系的砂卵石含水層呈不整合接觸，在玉泉山附近第四系地下水與奧陶系巖溶水水力聯(lián)系密切。

2 玉泉山泉水的歷史變化及斷流原因分析

2.1 玉泉山泉水的歷史變化

北京西郊的玉泉山位于海淀區(qū)昆明湖以西1.5 km，是古都北京近八百年來唯一持續(xù)供水的水源，對北京城的起源、發(fā)展和繁榮有重要意義。泉群主要出露于玉泉山東南山麓，由北向南依次有寶珠泉、五孔閘泉、試墨泉、裂帛泉、玉泉，西南山麓僅一迸珠泉，水由地下噴涌而出、淙淙有聲、水量極大，各泉涌出地表后形成泉湖并互相匯通。后經(jīng)乾隆驗證：“凡出山下而有冽者，城不如京師之玉泉”，并欽定為“天下第一泉”和“玉泉趵突”，列為燕京八景之一。據(jù)1928年、1934年考察資料，泉水總流量達到2.01 m3·s-1，1949年減少為5處，總流量1.54 m3·s-1，1951年僅為1.0 m3·s-1，1975年5月斷流（表1）。

據(jù)泉志資料記載，玉泉山泉眼出露高程為50.0 m左右，1974年泉水斷流前附近鉆孔地下水位為51.3 m，若泉水復(fù)涌，地下水位則需達到1974年的水平。

2.2 玉泉山泉水枯竭的原因分析

玉泉山泉水涌流了數(shù)百年，卻于20世紀70年代中期干涸了，其原因是多方面的：

（1）官廳水庫建成后，攔蓄了上游來水。官廳水庫自1954年5月完成竣工并調(diào)蓄運用，由于水庫大壩截流，永定河原有基流沒有了，河道每年斷流時間很長，河水滲漏量逐年減少，北京西郊地區(qū)的地下水補給量也隨著減少（秦大軍等，2019）。

（2）地下水開采量逐年增加。新中國成立后，隨著城市規(guī)模擴大，城市建設(shè)迅速發(fā)展，工農(nóng)業(yè)用水大量開采地下水，導致市區(qū)地下水位大幅度持續(xù)下降。

此外，華北地區(qū)自60年代初開始連年干旱少雨，多年沒有普降大雨也是原因之一。

3 泉水復(fù)涌途徑分析

由于玉泉山泉水處于北京西山奧陶系巖溶地下水系統(tǒng)，受特殊的地質(zhì)構(gòu)造影響，玉泉山泉水的復(fù)涌，需要由巖溶水的水位回升而實現(xiàn)。由于入滲條件和大氣降水量等客觀條件不易改變，但可借“南水北調(diào)”地表水源的到來，通過增加補給量、減少開采量，實現(xiàn)水位回升。

（1）增加永定河道入滲補給

依據(jù)多年數(shù)據(jù)統(tǒng)計，永定河雁翅—三家店段的最大滲漏量約為18 500萬m3·a-1，最小滲漏量約為3800萬m3·a-1，1960—2013年多年平均滲漏量約為9500萬m3。開展河道入滲實驗，利用水庫水源向永定河下葦?shù)殡娬鞠掠卧囼灮匮a，自2017年2月24日進行連續(xù)補水，截至5月26日，補水93 d，累計補水6200萬m3，有水河段長度最長達8.34 km，3次河道入滲試驗得到河道入滲強度平均值為0.61 m·d-1。綜合多年統(tǒng)計數(shù)據(jù)和試驗數(shù)據(jù)，并考慮到下游河道的人工防滲處理等因素，若將“南水北調(diào)”水源引至永定河山峽段自然入滲，永定河雁翅—三家店段可增加的入滲量為8000萬m3·a-1。

（2）地下水禁采或壓采

北京西山地區(qū)地下水開采量包括第四系地下水和基巖水。第四系模擬區(qū)包括房山區(qū)北部、門頭溝區(qū)東南部、石景山區(qū)、海淀區(qū)、豐臺區(qū)西北部分地區(qū)及昌平區(qū)西南部分地區(qū)，2017年第四系地下水年總開采量為25 730萬m3;奧陶系巖溶水模擬區(qū)包括海淀、石景山、門頭溝、豐臺西部、房山東北部以及昌平中南部。據(jù)統(tǒng)計，基巖裂隙類地下水開采總量為430萬m3·a-1;巖溶水自備井開采總量為6300萬m3·a-1，集中水源地巖溶水開采量為

5480萬m3·a-1，則基巖水總的開采量為12 220萬m3·a-1。集中水源地中水源三廠和水源四廠位于城近郊區(qū)，水源三廠2017年總開采量為6850萬m3·a-1，其中第四系3500萬m3·a-1，基巖水3350萬m3·a-1;水源四廠2017年開采量320萬m3·a-1，全部為第四系地下水。

隨著“南水北調(diào)”進京后，將逐步建設(shè)門城、石景山及溫泉等地表水廠，并由外調(diào)水源替代水源三廠等地下水源，因此可實現(xiàn)地下水大量減采以及管網(wǎng)覆蓋區(qū)禁采。

（3）地下水回灌

利用在西郊玉泉山地區(qū)新開挖的2眼輻射大口井進行入滲試驗研究，輻射井井徑3 m、井深18 m，井管8 m以上為死管、8～20 m井壁有輻射孔及20根2 m長的輻射管，進行了多次回灌量為280～680 m3·h-1的試驗，西郊玉泉山地區(qū)最終確定最優(yōu)的輻射井回灌能力為400 m3·h-1。1980年利用首鋼大口井開展回灌試驗，大口井深25 m，直徑8 m，試驗初期最大回灌量可達0.5 m3·s-1，試驗數(shù)據(jù)表明，大口井的單位面積入滲能力3.38～7.58 m·h-1。通過對西郊玉泉山巖溶水與第四系含水層直接接觸地區(qū)的調(diào)查，選定了大口井回灌適宜場地27處，通過永定河引水渠和南旱河，引“南水北調(diào)”水源進行回灌，考慮逐年維護清洗等因素，年回灌量可達到8000萬m3左右。

4 玉泉山泉水復(fù)涌方案預(yù)測研究

4.1 地下水數(shù)值模擬模型構(gòu)建

利用《北京巖溶水資源勘查評價工程成果報告》建立的地下水數(shù)值模擬模型，巖溶水數(shù)值模型范圍1226 km2，第四系含水層模擬面積1642 km2，模擬層為第四系含水層和奧陶系含水巖組，整個區(qū)域概化為5層：第四系潛水含水層、第四系弱透水層、第四系承壓水含水層、基巖隔水層和奧陶系含水巖組。奧陶系含水巖組東北邊界為南口-孫河斷裂，東南邊界為黃莊-高麗營斷裂，南部邊界在三福村—東莊子一線，西部邊界為大安山-紅煤廠斷裂，北界以百花山—髫髻山—妙峰山一線為界。第四系含水層模擬范圍以山前平原分界線為界，東部、東南和西南以第四系單一卵礫石層與多層砂卵礫石層過渡邊界為界。第四系含水層西部邊界為平原山區(qū)分界線，第四系含水層一般為單一潛水含水層，主要巖性為砂卵礫石，富水性好，易于接受山區(qū)降水的側(cè)向補給，其補給量與降水密切相關(guān)（楊平，1984），為流量邊界;東北部及東南部邊界在模型中被設(shè)定為流量邊界（圖3），其流量根據(jù)每年的地下水流場圖進行計算。依據(jù)本次研究成果重新率定了水文地質(zhì)參數(shù)、地下水源匯項，選用地下水模擬軟件GMS，建立了非均質(zhì)、各向異性的非穩(wěn)定三維地下水流模型，利用2017年數(shù)據(jù)對數(shù)值模型進行識別、驗證，并以2019年數(shù)據(jù)為初始流場，用于北京西郊地區(qū)地下水調(diào)蓄和玉泉山泉水復(fù)涌預(yù)測研究。

4.2 調(diào)蓄預(yù)測方案

以2017年6月實測地下水流場作為預(yù)測期的初始流場。預(yù)測期降水條件為平水年，年降水量510 mm。邊界流入量根據(jù)平水年條件計算相應(yīng)流入量。邊界流出量根據(jù)達西定律計算給定，在預(yù)測期內(nèi)認為邊界處水力梯度保持不變。第四系地下水與巖溶水之間的交換量由模型自動計算。方案預(yù)測期根據(jù)開采方案分為2019年6月—2030年末和2019年6月—2050年末。

目前，在平水年條件下，研究區(qū)通過永定河補水量為6276萬m3·a-1，第四系地下水開采量為34 514萬m3·a-1，巖溶水開采量11 785m3·a-1。隨著南水北調(diào)持續(xù)供水和規(guī)劃永定河陳家莊地表水水庫的修建，北京將在西山地區(qū)開展永定河滲漏補水和大口井回灌的地下水回補措施，并逐步壓采地下水。因此，模擬調(diào)蓄方案將以現(xiàn)狀開采方案為基礎(chǔ)，調(diào)整河道滲漏、大口井回灌補水量及水廠開采量。方案一為現(xiàn)狀開采方案，即在現(xiàn)狀補給和排泄條件下進行的方案預(yù)測。方案二為增滲+極限限采方案，補給條件為增加永定河山峽段滲漏量8000萬m3·a-1;巖溶水全部停采，海淀區(qū)區(qū)域第四系地下水全部停采。方案三為增滲+減采方案，永定河山峽段滲漏量在現(xiàn)狀基礎(chǔ)上增加8000萬m3·a-1;巖溶水集中水廠停采，四季青鎮(zhèn)區(qū)域巖溶水開采量停采，楊莊水廠和第三水廠第四系地下水全部停采。方案四為回灌+減采方案，在永定河引水渠和南旱河沿岸進行大口井回灌，回灌量為8000萬m3·a-1;僅巖溶水集中水廠停采，其余開采量保持現(xiàn)狀不變，即巖溶水開采量為6305萬m3·a-1，楊莊水廠和第三水廠第四系地下水停采，海淀區(qū)第四系地下水減采551萬m3·a-1。具體調(diào)蓄預(yù)測方案見表2。

4.3 預(yù)測及分析

根據(jù)各調(diào)蓄方案預(yù)測期滿后巖溶水流場圖顯示（圖4）：方案一條件下，巖溶水地下水位呈逐年下降趨勢，巖溶水地下水位不斷下降，到2030年玉泉山處水位約為8.3 m，每年下降約1.2 m;方案二條件下，巖溶水系統(tǒng)地下水位整體呈上升趨勢，平均水頭升幅為2.72 m·a-1，玉泉山泉水位到2030年末為52.1 m，每年水頭上升2.38 m。高于泉水最低出流水位51.6 m，到2030年末玉泉山泉水恢復(fù)出流;方案三條件下，巖溶水水頭上升約1.22 m·a-1，玉泉山泉水位到2050年末為52.0 m，每年水頭上升1.0 m，高于泉水最低出流水位51.6 m，到2050年末玉泉山泉水恢復(fù)出流;方案四條件下，巖溶水系統(tǒng)從補給區(qū)到排泄區(qū)，水頭均呈逐年上升趨勢，區(qū)域平均上升0.84 m·a-1，到2050年末，玉泉山泉水位51.9 m，達到出流水位（圖5）。

從以上4種方案預(yù)測的水位變化趨勢知，方案二、方案三、方案四均可以在不同階段實現(xiàn)玉泉山泉復(fù)涌。方案二是在永定河山峽段河道自然入滲，總需水量50.01億m3，同時需要全面限采地下水，但從實際情況出發(fā)，巖溶水地區(qū)地下水全部壓采比較難實現(xiàn)，因此，方案二僅在南水北調(diào)水源有充分保證率的條件下采用。方案三和方案四可在2050年實現(xiàn)玉泉山泉復(fù)涌，雖然時間上比方案二晚20年，但對水廠限采的要求不如方案二那么嚴格，易于執(zhí)行。在同樣的回灌量情況下，大口井回灌工程比較大，規(guī)劃陳家莊水庫的建設(shè)后在永定河山峽段增加河道滲漏補給，相對大口井回灌，工程量較小，從經(jīng)濟上來看，方案四要占優(yōu)。因此，綜合考慮，方案四可作為玉泉山泉復(fù)涌的調(diào)蓄方案。

5 玉泉山復(fù)涌后地質(zhì)環(huán)境影響

雖然可以通過回灌和減采地下水實現(xiàn)玉泉山泉水復(fù)涌，但會引起第四系地下水位的抬升（圖6），方案四永定河大寧水庫至蘆村一帶預(yù)測水位標高沒有超出1983年的限高水位，其他地區(qū)均超出了1983年水位，且超出的水位基本在5～10 m。這將引發(fā)如下的地質(zhì)環(huán)境問題：1）造成非正規(guī)的垃圾場被浸泡，使得地表污染物向含水層入滲的途徑變短，加劇對地下水產(chǎn)生污染（高詠等，2002）;2）重大建筑物的抗浮及地下水滲漏問題將凸顯出來（北京市水科學技術(shù)研究院等，2017）;3）砂土液化的可能性增大。

6 結(jié)論與建議

（1）玉泉山區(qū)域巖溶水主要補給來源為軍莊—雁翅地區(qū)巖溶裂隙水接受大氣降水補給、潭柘寺地區(qū)降水入滲補給以及永定河河水入滲補給。

（2）從調(diào)蓄預(yù)測方案的模擬結(jié)果知，要實現(xiàn)玉泉山的恢復(fù)，需要回灌和限采同時進行，在集中水廠巖溶水全部停采的條件下，大口井和河道回灌8000萬m3·a-1均可實現(xiàn)玉泉山泉的復(fù)涌，而規(guī)劃陳家莊水庫的建設(shè)使得永定河回灌的成本較低，可作為玉泉山復(fù)涌的主要方案。

（3）在實現(xiàn)泉水復(fù)涌的同時，也會帶來一定的地質(zhì)環(huán)境問題，例如引起地基不穩(wěn)，抗浮設(shè)防水位提高，地下水污染等問題。因此，在調(diào)蓄方案實施前，應(yīng)對以上問題開展專項研究工作。

（4）在玉泉山補給區(qū)域及周邊應(yīng)加強對地下水資源開采的限制，嚴禁超采、偷采等個、企行為，合理優(yōu)化水資源分配，涵養(yǎng)區(qū)域地下水源。

參考文獻：

北京市水文地質(zhì)工程地質(zhì)大隊，2017. 北京巖溶水資源勘查評價工程成果報告[R].

北京市水科學技術(shù)研究院，北京市地質(zhì)工程勘察院，2017.北京市西郊地區(qū)地下水戰(zhàn)略儲備關(guān)鍵技術(shù)研究與示范成果報告[R].

高詠，梁聰，車增光，2002. 南京典型段地鐵工程對地下水流動態(tài)的影響[J].科學技術(shù)與工程，20（28）：11711 -11717.

李裕宏，2005. 涵養(yǎng)京西地下水源恢復(fù)玉泉山泉流[J].北京規(guī)劃建設(shè)（6）：162-165.

秦大軍，孫杰，郭藝，等，2019. 永定河對北京西山巖溶水和玉泉山泉的影響[J].工程地質(zhì)學報，27（1）：162-169.

王莉蛟，張有全，宮輝力，等，2016. 北京市玉泉山泉恢復(fù)條件研究[J]. 水文地質(zhì)工程地質(zhì)，43（3）：22-28.

楊平， 1984.論北京玉泉山泉補給源：北京西山山前奧陶系巖溶水徑流特征[J].水文地質(zhì)工程地質(zhì)（2）：15-19.