環(huán)境同位素在識別伊犁河谷地下水補(bǔ)給來源中的指示作用

邵 杰，李 瑛，李樹才，滕 超，楊欣杰，曹 軍，陳喜慶

1.中國地質(zhì)調(diào)查局 自然資源實(shí)物地質(zhì)資料中心，河北 三河 065201；2.中國地質(zhì)調(diào)查局西安地質(zhì)調(diào)查中心，陜西 西安 710054

0 引言

水中各種元素的同位素是水的“指紋”，其記錄了地下水形成、補(bǔ)給、徑流與排泄過程中的信息.能否合理利用和有效保護(hù)水資源主要取決于人類對水循環(huán)的認(rèn)識水平.環(huán)境同位素可以用來有效地示蹤水循環(huán)，如指示水的來源，不同環(huán)境狀況下水的運(yùn)移和數(shù)量（包括江河湖泊），確定水的年齡，記錄水-巖相互作用，尤其在研究地下水和地下水儲量方面，更是卓有成效，從而為認(rèn)識地下水的形成、運(yùn)動及其成分變化機(jī)制提供重要的依據(jù)，為合理利用寶貴的水資源奠定基礎(chǔ)［1-5］.王文科等［6］借助同位素手段研究關(guān)中盆地地下水環(huán)境演化與可再生維持途徑.陳宗宇等［7］通過研究華北平原地下水的同位素分布特征對中國北方第四系地下水同位素分層及其指示意義進(jìn)行了研究探討.侯光才等［8］通過同位素水文地質(zhì)研究深化了鄂爾多斯盆地地下水的分布特征.Li等［9］通過對北京潮白河流域地下水氫氧穩(wěn)定同位素和水化學(xué)的分析，將地下水劃分為古地下水和現(xiàn)代地下水，并得到這兩種地下水具有不同的同位素組成和不同的補(bǔ)給環(huán)境，對北京市地下水管理具有借鑒意義.Ma等［10］通過對敦煌盆地地下水水化學(xué)和同位素研究，確定了含水層中地下水的水化學(xué)過程，辨識了地下水的補(bǔ)給環(huán)境和補(bǔ)給源以及停留時(shí)間，并得出承壓水的最近一次補(bǔ)給來自冰期的寒冷氣候條件下，而該區(qū)承壓水一段時(shí)間內(nèi)可能無法更新的結(jié)論.

1949年以來，在伊犁河谷相繼開展1∶10萬、1∶20萬、1∶50萬等不同精度的水文地質(zhì)測繪工作，基本查明了區(qū)域地下水的補(bǔ)、徑、排等條件.中國地質(zhì)調(diào)查局西安地質(zhì)調(diào)查中心開展的伊犁河谷地下水資源及其環(huán)境問題調(diào)查評價(jià)項(xiàng)目工作，為查明伊犁河谷地區(qū)含水層空間分布和結(jié)構(gòu)、區(qū)域地下水補(bǔ)徑排條件、地下水資源潛力評價(jià)、地下水資源合理開發(fā)利用等方面提供了資料支撐，并對伊犁河谷地區(qū)的地下水流系統(tǒng)及地下水循環(huán)演化開展了相關(guān)研究［11-14］.由于近年來關(guān)于研究區(qū)水體（地表水、潛水、承壓水、泉水等）同位素方面的研究開展的較少且不夠深入，在整理分析前人工作成果基礎(chǔ)上，本文試圖通過綜合分析研究區(qū)水體（地表水、潛水、承壓水、泉水等）氫氧、13C、14C同位素的空間分布特征，以此利用同位素證據(jù)探討地下水與地表水的補(bǔ)給來源、潛水與淺層承壓水的水力聯(lián)系、淺層承壓水與深層承壓水的水力聯(lián)系、深層承壓水的年齡及性質(zhì)等科學(xué)問題.

1 數(shù)據(jù)來源與研究方法

1.1 研究區(qū)概況

伊犁河谷地處我國新疆天山山脈西段，區(qū)內(nèi)地形地貌形態(tài)較復(fù)雜，總體呈“三山夾兩谷地一盆地”的格局，呈三角形，東部頂點(diǎn)為特克斯河與鞏乃斯河交匯處，西部底邊朝向中哈邊界，東高西低，地勢由東向西傾斜，東西長170 km.伊犁河谷地北部自西向東分布著北西走向的科古琴山、博羅科努山和依連哈比爾尕山等北天山支脈.山區(qū)面積占全區(qū)面積的70%以上.根據(jù)內(nèi)外營力地質(zhì)作用及地形特征，區(qū)內(nèi)主要地貌類型為冰蝕侵蝕褶皺斷塊山、剝蝕堆積塊狀隆起山和堆積平原三大類.區(qū)內(nèi)降水比較豐沛，具有濕潤大陸性中溫帶氣候的特征.氣候特點(diǎn)是冬夏兩季較長，春秋兩季較短，春季升溫迅速，夏季溫和多雨，秋季降溫較快，冬季較長較暖，日照時(shí)間長，光熱資源豐富.據(jù)《新疆維吾爾自治區(qū)水資源公報(bào)（2010）》，2010年地下水資源量（已扣除與地表水資源量間的重復(fù)量）為4.34×108m3，地下水實(shí)際開采量為4.15×108m3.

伊犁河谷為典型的受河谷控制的山間谷地，區(qū)內(nèi)第四紀(jì)地層相當(dāng)發(fā)育，包括各個(gè)不同時(shí)期的多種成因類型的堆積物.區(qū)內(nèi)主要發(fā)育四大河流，分別是鞏乃斯河、伊犁河、特克斯河以及喀什河.在四大河流控制作用下，形成了三大河谷區(qū)，分別是伊犁-鞏乃斯河谷、特克斯河谷及喀什河谷.伊犁谷地呈典型喇叭狀，西端開口處寬約90 km，沿伊犁河和鞏乃斯河向南東東延伸至東部合攏，為一斷陷谷地.第四系厚度由東向西增大，最厚達(dá)近千米，地層結(jié)構(gòu)多次由單層到多層交替，巖層顆粒也多次由粗到細(xì)變化［15］.地下水的補(bǔ)給、徑流、排泄條件受地形、地貌、地層、巖性結(jié)構(gòu)、構(gòu)造、水文、氣象等因素的制約.

在山前沖洪積扇區(qū)，地下水埋藏較深，深度大于100 m，而在河流沖洪積扇的河谷區(qū)兩側(cè)，地下水埋藏較淺，大多在0～15 m范圍內(nèi).由此可見，地下水的埋藏深度從山前強(qiáng)傾斜礫質(zhì)平原到山前緩傾斜含礫細(xì)土平原再到?jīng)_洪積細(xì)土平原逐漸變小.研究區(qū)北部及南部的山區(qū)是地下水的主要補(bǔ)給區(qū)，山前沖洪積或冰水沉積傾斜平原是地下水的主要徑流區(qū)，而盆地或谷地中部平原區(qū)則是地下水的主要排泄區(qū).在谷地兩側(cè)雖有南北向徑流，但整個(gè)谷地的地下水總流向是由東向西大體與鞏乃斯河及伊犁河的流向相一致①伊犁地區(qū)水利電力勘測設(shè)計(jì)院.新疆伊犁地區(qū)地下水資源報(bào)告.1990..

1.2 數(shù)據(jù)來源

樣品采樣過程中嚴(yán)格按照《環(huán)境同位素分析水樣采集》導(dǎo)則進(jìn)行.2H和18O樣品的采集使用30 mL PET棕色塑料瓶采集，取樣前使用經(jīng)過0.02μm濾膜過濾的樣品，清洗不少于3次，采集必須使瓶內(nèi)不留氣泡，采集完成后使用封口膜進(jìn)行外圍瓶蓋封口.分不同水源種類進(jìn)行采集，其中氫氧穩(wěn)定同位素潛水樣品59組、淺層承壓水樣品25組、泉水樣品17組、地表水樣品19組，共計(jì)120組.13C穩(wěn)定同位素地下水樣品5組，5組樣品采自研究區(qū)水文鉆探井，14C深層承壓水樣3組.樣品測試工作由自然資源部地下水礦泉水及環(huán)境監(jiān)測中心完成.

1.3 研究方法

采用把地下水δD-δ18O圖上的數(shù)據(jù)繪成平面圖的方法，并附以全球與地方大氣降水線作為參考.將所測水體氫氧穩(wěn)定同位素樣品置于全球雨水線以及當(dāng)?shù)赜晁€中，分析伊犁河谷地下水、地表水中的δD與δ18O之間的關(guān)系，總結(jié)地下水、地表水中δD與δ18O在不同地貌單元所呈現(xiàn)出的分布特征.并且總結(jié)其對地下水在補(bǔ)給、徑流、排泄以及不同層位地下水的水力聯(lián)系的指示意義.13C穩(wěn)定同位素廣泛參與碳循環(huán)，通過分析地下水樣品δ13C的最大值、最小值、平均值、極差等特征值，分析地下水13C穩(wěn)定同位素的來源，進(jìn)而識別地下水的補(bǔ)給來源［16］.應(yīng)用14C估算深層承壓水年齡，以此探討分析深層承壓水的更新速率及其與淺層承壓水的水力聯(lián)系強(qiáng)度.

2 結(jié)果與分析

2.1 δD-δ18O分析

距離伊犁河谷最近的全球降水同位素網(wǎng)（GNIP）監(jiān)測站為烏魯木齊站.據(jù)1986—2002年大氣降水的同位素組成，當(dāng)?shù)赜晁€方程為：δD＝7.2δ18O＋4.5，由于氣候干燥、蒸發(fā)強(qiáng)烈，與全球雨水線相比斜率稍?。?7］.

δD-δ18O樣品測試結(jié)果特征值見表1.59組潛水樣埋深范圍為2.5～40.0 m，25組淺層承壓水樣埋深為40.0～170.0 m.據(jù)表1及圖1，將潛水與淺層承壓水的氫氧穩(wěn)定同位素進(jìn)行對比分析發(fā)現(xiàn)：潛水和淺層承壓水樣品δD、δ18O變化范圍、極差、平均值均較為接近，且潛水和淺層承壓水氫氧穩(wěn)定同位素大多數(shù)位于全球雨水線的左上方、當(dāng)?shù)卮髿庥晁€附近（圖1），表明潛水的補(bǔ)給來源為當(dāng)?shù)卮髿饨邓已a(bǔ)給來源曾經(jīng)歷了較強(qiáng)的蒸發(fā)作用而引起氫氧重同位素富集.這表明二者相互轉(zhuǎn)換頻繁，具有較強(qiáng)的水力聯(lián)系.泉水與地表水的氫氧穩(wěn)定同位素進(jìn)行對比分析發(fā)現(xiàn)：泉水和地表水樣品δD、δ18O變化范圍、極差、平均值均較為接近，且泉水和地表水氫氧穩(wěn)定同位素大多數(shù)位于全球雨水線的左上方及當(dāng)?shù)卮髿庥晁€附近（圖2），表明地表水的補(bǔ)給來源為當(dāng)?shù)卮髿饨邓?，且補(bǔ)給來源曾經(jīng)歷了較強(qiáng)的蒸發(fā)作用而引起氫氧重同位素富集.

圖1 研究區(qū)潛水及承壓水δD-δ18O關(guān)系Fig.1 TheδD-δ18O relationship of phreatic and confined water in the study area

圖2 研究區(qū)泉水及地表水δD-δ18O關(guān)系Fig.2 TheδD-δ18O relationship of spring and surface water in the study area

表1 氫氧穩(wěn)定同位素參數(shù)統(tǒng)計(jì)特征值表Table 1 Statistical eigenvalues of hydrogen and oxygen stable isotope parameters

2.2 13C分析

研究區(qū)樣品13C穩(wěn)定同位素分析測試結(jié)果見表2.δ13C的變化范圍為-9.4‰～-5.6‰，平均值為-8.3‰，極差為3.8‰，最小值為-9.4‰，最大值為-5.6‰.樣品取樣深度范圍為1.9～198.0 m，取樣深度跨度為197.1 m，而取樣樣品δ13C變化范圍較小.由此可見，研究區(qū)不同深度地下水的δ13C值較為接近，且與CO2的δ13C值-7‰相近，說明研究區(qū)地下水中碳的主要來源為大氣CO2.

表2 13C穩(wěn)定同位素分析測試結(jié)果表Table 2 Stable isotopic analysis and testing results of 13C

2.3 14C分析

在研究區(qū)采集3組深層承壓水，水位埋深在198.0～240.0 m.由表3可知，承壓水在埋深240 m左右，地下水校正年齡為23 ka，埋深在200 m左右，地下水校正年齡為19.5 ka.與淺層承壓水相比較，深層承壓水年齡較大，屬于沉積埋藏水，這也從側(cè)面反映該地區(qū)的深層承壓水與淺層承壓水的水力聯(lián)系較弱.

表3 承壓水14C測試結(jié)果Table 3 14C radioactive isotope test results of confined water

3 對比分析

近年來，相關(guān)專家學(xué)者利用同位素?cái)?shù)據(jù)，在新疆及伊犁河谷地區(qū)識別地下水補(bǔ)給來源問題.李捷等［18］利用大氣降水監(jiān)測站δD、δ18O多年平均值與當(dāng)?shù)氐叵滤摩腄、δ18O，通過分析天山地區(qū)受不同降水水汽來源影響對地下水補(bǔ)給來源進(jìn)行識別，認(rèn)為深層地下水同位素相對貧化，與補(bǔ)給期相對寒冷的氣候環(huán)境有關(guān)，東天山為該區(qū)地下水的主要補(bǔ)給來源，且克拉美麗山附近深層地下水與淺層地下水同位素均貧化，與阿勒泰站大氣降水同位素特征相符，該區(qū)主要接受克拉美麗山的補(bǔ)給.王騫迎等［19］利用氫氧穩(wěn)定同位素?cái)?shù)據(jù)與當(dāng)?shù)卮髿饨邓€方程進(jìn)行對比分析，表明研究區(qū)內(nèi)地下水、伊犁河及支流河水樣品集中分布在當(dāng)?shù)厣絽^(qū)大氣降水線與平原區(qū)大氣降水線之間，略靠近山區(qū)降水線，說明研究區(qū)內(nèi)地下水、伊犁河及其支流主要接受山區(qū)大氣降水補(bǔ)給.本文研究結(jié)果顯示，研究區(qū)潛水與承壓水均起源于當(dāng)?shù)卮髿饨邓覐膮^(qū)域上進(jìn)一步表明了潛水與承壓水的水力聯(lián)系強(qiáng)度特征以及深層承壓水的年齡與性質(zhì).

4 討論

環(huán)境同位素技術(shù)廣泛地應(yīng)用于地下水循環(huán)演化研究中，通過分析氫氧同位素分布特征，可以查明氫氧穩(wěn)定同位素與當(dāng)?shù)卮髿饨邓凰氐年P(guān)系.本研究對泉水、潛水、淺層承壓水及地表水氫氧穩(wěn)定同位素特征分析表明：4種不同種類水源的氫氧穩(wěn)定同位素分布特征與當(dāng)?shù)卮髿饨邓窒嘟?，且δD與δ18O變化范圍很相似，揭示了它們均來源當(dāng)?shù)卮髿饨邓?

碳在天然有機(jī)循環(huán)過程中，同位素組成的變化范圍超過10%.一般來說，在碳的有機(jī)循環(huán)中，輕同位素傾向于富集在有機(jī)質(zhì)（如石油等富含12C）中，相反在碳的無機(jī)循環(huán)中，重同位素傾向于富集在無機(jī)鹽（如碳酸鹽富含13C）中.上述兩種碳循環(huán)均與大氣二氧化碳有密切關(guān)系，也是自然界中碳同位素分餾的兩個(gè)最重要的過程［20］.地下水碳同位素組成受制于地下水本身形成、遷移及賦存的環(huán)境.地下水中碳的來源主要有大氣CO2、土壤和現(xiàn)代生物碳，一般條件下，大氣CO2的δ13C為-7‰左右，土壤和現(xiàn)代生物碳的δ13C為-25‰左右.研究區(qū)不同埋藏深度地下水的δ3C值較為接近，且與大氣CO2值相近.研究區(qū)地下水中的碳主要來源于大氣CO2，說明地下水中的13C源自于當(dāng)?shù)卮髿饨邓?

5 結(jié)論

（1）潛水δD變化范圍為-97.32‰～-67.51‰，δ18O為-15.85‰～-10.66‰，淺層承壓水δD變化范圍為-111.93‰～-68.38‰，δ18O為-16.01‰～-10.52‰，并且潛水與淺層承壓水的補(bǔ)給來源經(jīng)歷了較強(qiáng)的蒸發(fā)作用而引起氫氧重同位素富集.潛水與淺層承壓水相互轉(zhuǎn)化頻繁，水力聯(lián)系較強(qiáng).潛水與淺層承壓水屬于同一地下水含水系統(tǒng).與淺層承壓水相比，深層承壓水年齡較大，屬于沉積埋藏水.

（2）泉水δD變化范圍為-102.06‰～-71.63‰，δ18O為-14.21‰～-9.70‰，地表水δD變化范圍為-90.53‰～-60.99‰，δ18O為-13.20‰～-9.54‰.4種不同種類水源的氫氧穩(wěn)定同位素分布特征與當(dāng)?shù)卮髿饨邓窒嘟?，且δD與δ18O變化范圍很相似，說明均源自于當(dāng)?shù)卮髿饨邓?

（3）δ13C的變化范圍為-9.4‰～-5.6‰，平均值為-8.3‰，極差為3.8‰，不同埋藏深度地下水的δ13C值較為接近，與CO2的δ13C值相近.研究區(qū)地下水中的碳主要來源于大氣CO2，指示了地下水中的13C來源當(dāng)?shù)卮髿饨邓?