李 露，秦大軍,3，郭 藝，孫 杰

(1.中國科學院地質(zhì)與地球物理研究所/中國科學院頁巖氣與地質(zhì)工程重點實驗室，北京 100029；2.中國科學院大學地球與行星科學學院，北京100049；3.三亞學院翟明國院士工作站，海南 三亞 572022)

中國北方巖溶水富集、水質(zhì)良好，在城市供水中發(fā)揮關(guān)鍵作用。巖溶泉為河水、湖泊、濕地的源頭，具有不可或缺的生態(tài)功能。北京西山巖溶水是北京市重要的水源地，已有研究工作多側(cè)重于補給來源、水流特征[1]，水資源評價、循環(huán)特征、可更新能力[2]和優(yōu)勢地下水流通道、巖溶系統(tǒng)的不均一性以及深埋巖溶水系統(tǒng)導水機制[3]等方面。而西山巖溶水水質(zhì)狀況及影響因素等方面的研究較少。

1 水文地質(zhì)條件

研究區(qū)位于北京市西部，包括山區(qū)和平原區(qū)兩部分，地勢西北高、東南低，總面積約為280 km2。屬半干旱溫帶大陸性季風氣候，夏季炎熱多雨，冬季寒冷干燥。年平均氣溫為11.5℃，月平均最高氣溫為28.2℃(7月)，月平均最低氣溫為-16.8℃(1月)。年平均降水量為600 mm，降水分布很不均勻，超過85%的降水集中在7—9月。海河水系的支流永定河，貫穿整個研究區(qū)域。

北京西山巖溶水賦水層位為寒武系、奧陶系灰?guī)r，上覆石炭系和二疊系砂巖、黏土巖，下伏地層為青白口薄層灰?guī)r夾頁巖。區(qū)內(nèi)褶皺和斷裂構(gòu)造發(fā)育(圖1)，香峪向斜、八大處背斜、永定河斷裂、八寶山斷裂和黃莊斷裂控制著巖溶含水層的導水性。大氣降水和地表水是區(qū)域巖溶含水層的主要補給來源，巖溶水總體流向為西南向東北，與構(gòu)造走向近于平行，以人工開采和徑流(局部地區(qū)對第四系地下水的頂托補給、側(cè)向徑流)的方式排泄。區(qū)域巖溶水有軍莊和潭柘寺兩大補給區(qū)，軍莊地區(qū)主要接受大氣降水和永定河河水補給，潭柘寺地區(qū)主要接受大氣降水的垂直入滲補給；徑流過程中香峪向斜具有明顯的阻水作用，軍莊地區(qū)巖溶水不能以深循環(huán)的形式通過該向斜，永定河斷裂則是重要的導水通道，軍莊地區(qū)巖溶水接受補給后，分兩條路徑進行排泄：(1)沿永定河斷裂向八寶山斷裂進行徑流，并最終向玉泉山泉地區(qū)排泄；(2)沿著香峪向斜北翼向永豐屯地區(qū)徑流?；诖藦搅鳁l件的認識，該區(qū)域巖溶水的徑流路徑主要有3條：軍莊—永豐屯、軍莊—古城—玉泉山泉、潭柘寺—四季青[3]。

圖1 北京西山水文地質(zhì)及水樣品分布圖Fig.1 Hydrogeological map of Xishan karst groundwater system in Beijing and cocation of the sampling sites1—露頭奧陶系灰?guī)r；2—第四系沉積物覆蓋下的奧陶系灰?guī)r；3—石炭系-二疊系砂巖夾黏土巖；4—第四系沉積物；5—斷裂；6—向斜；7—泉；8—河流；9—巖溶水樣品；10—第四系水樣品；11—河水樣品

2 樣品采集與測試

本次共采集水樣品47個(圖1)。河水樣品4個，其中3個樣品(65、66、68)的取樣時間為2012年11月，樣品69取樣時間為2013年6月11日。第四系水樣品7個，其中6個樣品的取樣時間為2012年6—9月，樣品71取樣時間為2012年12月5日。巖溶水樣品共采集36個，樣品74取樣時間為2013年1月5日，其余35個樣品的取樣時間為2012年6—9月。

采樣過程中，使用WTW便攜式多參數(shù)測量儀現(xiàn)場測試了T、DO、EC、pH等參數(shù)。在實驗室內(nèi)，使用離子色譜儀測定水化學離子含量，D和18O穩(wěn)定同位素采用MAT253質(zhì)譜儀測定，結(jié)果以相對于標準平均海洋水比率的千分偏差(δ值)表示。硝酸鹽的15N穩(wěn)定同位素采用反硝化細菌法測定，結(jié)果以相對于大氣中氮氣比率的千分偏差(δ值)表示。

3 結(jié)果與討論

3.1 氫氧穩(wěn)定同位素特征

樣品的數(shù)據(jù)統(tǒng)計結(jié)果見表1。河水樣品中δD的最小值為-61.4‰，最大值為-36‰，均值為-50.47‰；δ18O的最小值為-7.52‰，最大值為-5.66‰，均值為-6.32‰。第四系水樣品中δD最小值為-77‰，最大值為-49.4‰，均值為-60.21‰；δ18O的最小值為-10.61‰，最大值為-5.87‰，均值為-7.84‰。巖溶水樣品中δD最小值為-75.5‰，最大值為-55.4‰，均值為-64.49‰；δ18O的最小值為-10.27‰，最大值為-6.47‰，均值為-8.65‰?？傮w看，水樣的氫氧同位素富集特征為：河水>第四系水>巖溶水。河水受蒸發(fā)作用強烈，氫氧穩(wěn)定同位素較第四系水和巖溶水更加富集。

表1 樣品水化學和同位素數(shù)據(jù)統(tǒng)計表Table 1 Statistical results of hydrochemical and isotopic data

由圖2可知，氫氧同位素點均位于全球大氣降水線GMWL[22]下方。河水樣品66，采自三家店水庫，為水庫滯留水，受強烈的蒸發(fā)作用，氫氧同位素富集，位于全球大氣降水線上。LMWL線(y=7.6x+ 3.8)為當?shù)卮髿饨邓€[3]。依據(jù)水樣中氫氧同位素值點在圖中的分布情況，將北京西山地區(qū)水樣品分為三組：Ⅰ組水樣點位于當?shù)卮髿饨邓€(LMWL)上，包括部分巖溶水和3個第四系水。Ⅱ組水樣點均為巖溶水。Ⅲ組水樣點包括3個河水、4個第四系水和5個巖溶水。

圖2 水樣品的氫氧同位素關(guān)系Fig.2 Relationship between hydrogen and oxygen isotopes in water samples from the Xishan area

3.2 巖溶水濃度和分布

圖3 西山巖溶水分布圖Fig.3 Distribution of the Xishan karst aquifer1—露頭奧陶系灰?guī)r；2—第四系沉積物覆蓋下的奧陶系灰?guī)r；3—石炭系—二疊系砂巖夾黏土巖；4—第四系沉積物；5—斷裂；6—向斜；7—泉；8—河流

Ni=N0×(Ci/C0)

(1)

C0——徑流途徑上第1個樣品的Cl-濃度；

Ci——第i個樣品的Cl-濃度；

圖4 軍莊—永豐屯徑流濃度計算值與測量值比較Fig.4 Comparison of the calculated along the Junzhuang—Yongfengtun flow path

圖5 軍莊—古城—玉泉山泉徑流濃度計算值與測量值比較Fig.5 Comparison of the calculated along the Junzhuang—Gucheng—Yuquanshan Spring flow path

圖6 潭柘寺—四季青徑流濃度計算值與測量值比較Fig.6 Comparison of the calculated along the Tanzhesi—Sijiqing flow path

3.4 西山巖溶水中來源

圖7 巖溶水δ15N值與之間關(guān)系Fig.7 Relationship of δ15N and in the karst groundwater▲:軍莊-古城-玉泉山；◇:潭柘寺-四季青；○:軍莊-永豐屯

(2)

3.5 巖溶水中不同來源占比估算

(3)

式中：n——組合數(shù)量；

i——增量參數(shù)；

s——來源數(shù)量。

圖8 巖溶水中來源貢獻率Fig.8 Percent fraction of sources in the karst groundwater

4 結(jié)論