王廣昊， 張 瑩， 徐亮亮， 姜紀(jì)沂， 李春歡

(1.防災(zāi)科技學(xué)院生態(tài)環(huán)境學(xué)院， 三河 065201; 2.中國冶金地質(zhì)總局浙江地質(zhì)勘查院， 衢州 324000)

生態(tài)文明建設(shè)作為中國“五位一體”的重要內(nèi)容，其戰(zhàn)略地位更加明確，作為全國首批“金山銀山就是綠水青山”試驗基地的衢州市也把生態(tài)建設(shè)放在了首位，因此對地下水和地表水生態(tài)資源的開發(fā)利用、管理及保護(hù)提出了更高要求。衢江是衢州市重要的地表水資源，具有總量豐富但人均占有量低、空間分布與耕地、人口、經(jīng)濟(jì)發(fā)展不匹配的特點[1]。受該區(qū)域地理位置、氣候等影響，降水量年內(nèi)分配極不均勻，降水集中且多為暴雨，豐枯季節(jié)相差懸殊，造成梅雨時多發(fā)洪災(zāi)，伏季易發(fā)旱災(zāi)，旱季供水短缺，加劇了水資源的供需矛盾。隨著衢州市現(xiàn)代化建設(shè)和經(jīng)濟(jì)的迅速發(fā)展，用水量激增，水污染日趨嚴(yán)重，水資源短缺與水環(huán)境惡化已成為制約當(dāng)?shù)厣鷳B(tài)環(huán)境建設(shè)的重要因素。在自然變化和人類活動的影響下，水循環(huán)過程和水資源量分配發(fā)生了改變，而查清地表水與地下水的轉(zhuǎn)化關(guān)系有利于了解自然和人類影響下的生態(tài)及水循環(huán)問題，為地表水和地下水污染防治提供參考依據(jù)[2-5]。衢江作為錢塘江流域的源頭之一，在地表水流經(jīng)金衢盆地時流速相對緩慢，造成了地表水與地下水間轉(zhuǎn)換十分頻繁[6-8]，但前人對于該流域地表水與地下水的關(guān)系并沒有取得定量的認(rèn)識，導(dǎo)致流域水資源問題的不確定性和多解性，因此開展衢江流域地表水與地下水相互轉(zhuǎn)化關(guān)系的研究對該流域水資源管理、保護(hù)及合理利用，維持河流系統(tǒng)基本功能及生態(tài)健康具有重要意義。

目前，對于地表水與地下水之間轉(zhuǎn)化關(guān)系的研究方法主要有調(diào)查分析法、水量平衡法、模擬計算法和環(huán)境示蹤法等[9]。環(huán)境示蹤法具有原理簡單、實用有效、可結(jié)合多種信息源數(shù)據(jù)進(jìn)行分析推算等優(yōu)點，得到了廣泛運用，并取得了很好的應(yīng)用效果[10]，而環(huán)境示蹤法中的環(huán)境同位素和水化學(xué)技術(shù)是描述地表水與地下水轉(zhuǎn)化關(guān)系的有效手段[11-16]。

在地表水與地下水轉(zhuǎn)換的過程中，伴隨著物質(zhì)和能量的轉(zhuǎn)換，而環(huán)境同位素和水化學(xué)等環(huán)境示蹤劑作為地表水和地下水的重要組成部分，在一定程度上可以指示兩者的轉(zhuǎn)換關(guān)系。如王希義等[17]利用同位素技術(shù)研究了塔里木河下游河水對地下水的轉(zhuǎn)換關(guān)系。孫從建[18]研究了海河源區(qū)不同季節(jié)的地下水與地表水轉(zhuǎn)換關(guān)系。Blarasin等[19]、蘇小四等[20]、陳彭等[21]、郭亞文等[22]利用氫、氧穩(wěn)定同位素和水化學(xué)指示了不同地區(qū)地表水和地下水的相互作用。然而，有關(guān)地下水與地表水的研究大多集中在平原區(qū)或地形地貌變化不大的區(qū)域，針對地質(zhì)地形復(fù)雜的基巖區(qū)特別是在衢江流域分布廣泛的紅層區(qū)研究比較少。

在查明衢江游流域水工環(huán)綜合地質(zhì)條件的基礎(chǔ)上，對地表水和地下水進(jìn)行取樣，然后在室內(nèi)對水化學(xué)和同位素D、18O進(jìn)行分析測定，最后結(jié)合采樣現(xiàn)場測試數(shù)據(jù)，定性揭示出研究區(qū)水化學(xué)和氫氧同位素空間分布規(guī)律，定量地計算出地表水與地下水的轉(zhuǎn)化關(guān)系，為衢江流域地表水和地下水資源管理、保護(hù)及合理利用，維持河流系統(tǒng)基本功能及生態(tài)健康提供一定的參考。

1 研究區(qū)概況

衢江是錢塘江流域的主要支流和源頭之一，起源于安徽省休寧市，在浙江省內(nèi)流域面積為1.11×104km2。本次評價區(qū)所屬流域為衢江區(qū)域，研究區(qū)包括開化縣、常山縣、衢江區(qū)、柯城區(qū)、龍游縣等。

衢州市總體地勢特征為“簸箕”狀，即南部、北部和西部較高，中部和東部較低，中部為金衢盆地的西半部，自西向東逐漸展寬[23]。研究區(qū)上游位于盆地邊緣多為低山、高丘、崗地等地貌，中下游在金衢盆地內(nèi)，以紅色丘陵盆地、殘丘及沖積平原為主體?！叭姝h(huán)山夾一川，盆地錯落涵三江"是金衢盆地的基本地貌特征，地面高程一般為50～300 m。

研究區(qū)多年平均降水量為1 424 mm，降水量較充沛，但是從降水分布的時間上來看極其不均衡。其中大部分的降水是從5月初持續(xù)到6月底，總降雨量在500～610 mm，雨量大的影響因素將會造成洪澇災(zāi)害。夏季臺風(fēng)很難進(jìn)入衢江流域，因此多炎熱天氣。

研究區(qū)上游出露地層主要是早古生代奧陶紀(jì)、寒武紀(jì)地層，中下游地層為白堊系上統(tǒng)地層，第四系地層為中上更新統(tǒng)沖積層、沖洪積以及全新統(tǒng)沖積層。根據(jù)含水介質(zhì)和賦存形式，研究區(qū)地下水類型主要為松散巖類孔隙潛水及紅層碎屑巖類孔隙裂隙水和基巖裂隙水，紅層碎屑巖類孔隙裂隙水上部往往發(fā)育松散巖類孔隙潛水，衢州市水文地質(zhì)圖如圖1所示。

圖1 研究區(qū)水文地質(zhì)圖Fig.1 Hydrogeological zoning of the study area

2 樣品的采集和分析方法

2.1 樣品的采集

河流交匯處一般河水面比較開闊，水流量比較大，地表水與地下水的轉(zhuǎn)化關(guān)系也比較強(qiáng)烈，而且此區(qū)域人流量比較大，對于研究地表水與地下水的轉(zhuǎn)化關(guān)系更為迫切，所以本文采樣點布設(shè)原則為：在衢江流域主要二級支流匯入一級支流前及匯入一級支流后，一級支流匯入干流前及匯入干流后布設(shè)采樣斷面，對地表水、井水、泉水進(jìn)行采樣，地表水與地下水之間距離不超過1 km。在2019年9—10月期間共采集地表水樣67 個，地下水樣61 個，泉水樣8 個，取樣點位置和交通位置圖如圖2所示。

圖2 研究區(qū)位置和采樣點Fig.2 Location and sampling points of the study area

2.2 樣品測定

在采樣現(xiàn)場用哈希公司便攜式水質(zhì)儀(HQ40d，HACH)測定了水樣的pH、氧化還原電位Eh、溶解性總固體TDS、溶解氧DO、電導(dǎo)率Ec等指標(biāo)；水樣的陰陽離子送往浙江中一檢測研究院股份有限公司檢測；水樣的δD和δ18O值采用美國LGR(Los Gatos Research)公司生產(chǎn)的同位素分析儀(DLT-100)進(jìn)行測定，δD和δ18O的精度分別達(dá)到0.3‰和0.1‰。測定結(jié)果用VSMOW進(jìn)行校正表示。數(shù)據(jù)使用儀器配套軟件(LWIA Post Analysis v3.1.0.9 Installer)進(jìn)行分析處理。

2.3 水化學(xué)分析方法

運用ArcGIS10.6 軟件中的反權(quán)重距離插值法對研究區(qū)地表河水及地下水的總固體溶解物(total dissowed solids, TDS)、pH以及δD和δ18O做空間插值分析，以分析地表水及地下水TDS、pH、δD和δ18O的時空分布特征。Piper三線圖是目前研究和判斷水體水化學(xué)類型及各離子組分含量比例關(guān)系的重要手段[24]，運用Aq·QA軟件繪制出Piper三線圖，對研究區(qū)水化學(xué)空間特征分析做了比較。通過Gibbs法判讀水體的水化學(xué)特征是受控于巖石風(fēng)化影響、降水影響、蒸發(fā)富集中的哪幾種類型[25]。

2.4 地表水與地下水轉(zhuǎn)化定量分析

在氧的同位素當(dāng)中，以16O 和18O的豐度較高，彼此間的質(zhì)量差也較大，因此人們大都使用18O/16O 比值來表示δ18O[26]，即

δ18O=[(18O/16O)m-(18O/16O)s]/(18O/16O)s×1 000/‰

(1)

式(1)中: m為樣品; s為標(biāo)準(zhǔn)樣。

根據(jù)質(zhì)量守恒原理利用氫氧穩(wěn)定同位素可建立質(zhì)量均衡模型來估算地表水與地下水之間的轉(zhuǎn)化強(qiáng)度[27]。

δ18O 在水體轉(zhuǎn)化中的質(zhì)量守恒方程為

CsQs=CgQg+Cb(Qs-Qg)

(2)

式(2)中:Cs為取樣點河水的δ18O 值;Cg為取樣點地下水的δ18O值;Cb為取樣點上游來水中δ18O值;Qs為取樣點的河水流量;Qg為向地下水的排泄量。則由式(2)可以導(dǎo)出河水向地下水的排泄量占河水流量的百分比，即

f=(Qg/Qs)×100%=

(Cs-Cb)/(Cg-Cb)×100%

(3)

當(dāng)僅靠一種示蹤劑無法計算出其中一種混合來源的百分比時，必須還要借助另外一種示蹤劑，即分析水中另外一種同位素或水化學(xué)離子等。本研究借助氫同位素(δD) 按照三相混合來計算[28]，式( 4) 擴(kuò)展成三相混合質(zhì)量均衡方程，即

CsQs=CgQg+CcQa+Cm(Qs-Qg-Qa)

(4)

式(4)中：Cc為取樣點上游主干流水樣中δD值；Cm為取樣點上游匯入支流水樣中δD值。

3 結(jié)果與討論

3.1 溶解性總固體和pH變化特征

溶解性總固體(total dissolved solids，TDS)是指水中溶質(zhì)的總量，其濃度變化可以反映水的運移路徑和滯留長短。水在運移過程中會不斷地與土壤和圍巖中的溶解性鹽類發(fā)生離子交換，在沒有與TDS較小的水體混合、蒸散發(fā)強(qiáng)烈等情況下，水體的TDS是逐漸升高的[29]。因此，根據(jù)TDS在流域內(nèi)的分布情況，可以在一定程度上推斷流域內(nèi)地表水與地下水的轉(zhuǎn)換關(guān)系[30]。

研究區(qū)地表水TDS變化范圍是32～329 mg/L，地下水的TDS變化范圍是16～368 mg/L[圖3(a)、圖3(b)]。在研究區(qū)上游地下水溶解性總固體比地表水小，說明在此區(qū)域主要是地下水對地表水進(jìn)行補給；中下游河段地表水和地下水的電導(dǎo)率較為接近，但總體上來看地表水的溶解性總固體要比地下水的溶解性總固體略小，表明在中下游區(qū)域主要是地表水補給地下水的。進(jìn)一步分析其沿程變化可以看出，池淮溪和馬金溪上游河段的TDS緩慢變高，這是由于上游區(qū)地質(zhì)結(jié)構(gòu)比較復(fù)雜，且大部分河段區(qū)河床基巖出露，河水不斷沖刷溶解巖體中的溶解性離子并與其發(fā)生交換，造成了TDS的不斷累積；河流位于中下游區(qū)域時，TDS迅速升高，這是由于河道變得平緩，且中游區(qū)基本位于衢州市中心部位，受人為因素影響比較大，造成了區(qū)域性TDS迅速升高的現(xiàn)象，中下游的地下水TDS升高明顯，這是由于江山港、烏溪江等中下游區(qū)域地下水類型主要為紅層孔隙裂隙水，此區(qū)域人口密集，對于地下水的開采量比較大，使得地表水對地下水的轉(zhuǎn)化更為強(qiáng)烈，而紅層巖土體中的可溶解性離子更易溶解，這就造成了此區(qū)域地下水中的TDS明顯升高。

圖3 研究區(qū)地表水、地下水TDS和pH空間分布圖Fig.3 Spatial distribution of surface water, groundwater TDS and pH in the study area

研究地表水和地下水pH的變化范圍均為7.0～8.5，整體上看兩者的變化趨勢也比較相似[圖3(c)、圖3(d)]，尤其在中下游河流入?yún)R口處，進(jìn)一步驗證了此區(qū)域的地下水與地表水的轉(zhuǎn)換比較頻繁?？傮w上看，人為因素對本研究區(qū)的pH更為顯著，凡是人口密集處或有大型工廠處pH均變化顯著，比如江山港、烏溪江河段分布數(shù)個大型化學(xué)企業(yè)，且人口相對密集，衢江段下游區(qū)分布數(shù)個造紙廠等，這些都說明了研究區(qū)的人類活動是影響pH變化的主要因素。

3.2 陰陽離子組成分析

研究區(qū)水中離子不僅受沿途地區(qū)土壤和巖性的影響，與離子間的水化學(xué)作用也有很大關(guān)系，結(jié)合流域內(nèi)的地形地貌條件，可以大致判斷出地表水和地下水轉(zhuǎn)化路徑[31]。

從圖4中可以看出，衢江流域上游主要為HCO3-Ca、HCO3·Cl-Ca類型水，中下游主要為HCO3-Ca類型水，均屬于低礦化度水，說明了流域內(nèi)大氣降水是最終補給源。地下水與地表水的水化學(xué)類型相似，在一定程度上表明了兩者緊密的水力聯(lián)系。研究區(qū)上游水化學(xué)類型相對分散[圖4(a)]，是由于上游地勢地貌、地質(zhì)和水文地質(zhì)情況比較復(fù)雜，地表水和地下水中的離子在運移過程中受到了一定的影響，到中下游河段時，地勢變的平坦，地下水水力坡度也變低，離子在此時滯留的時間也比較長，使得水化學(xué)類型相對集中。紅層裂隙的各向異性和不均勻性造成了含水系統(tǒng)的復(fù)雜性，但此區(qū)域的地下水類型基本相似，這進(jìn)一步說明了此區(qū)域地表水與地下水之間轉(zhuǎn)換關(guān)系密切。

圖4 研究區(qū)Piper三線圖Fig.4 Trigraph of Piper in the study area

圖5 衢江流域水化學(xué)Gibbs圖Fig.5 Chemical Gibbs map of Qu River Basin water

3.3 氫氧穩(wěn)定同位素空間分布規(guī)律和特征分析

3.3.1 氫氧穩(wěn)定同位素空間分布規(guī)律

研究區(qū)地表水穩(wěn)定同位素δD的范圍是-75.45‰～-37.35‰，均值為-46.52‰，氘盈余的范圍是0.41‰～25.79‰，均值為12.47‰，地表水穩(wěn)定同位素δ18O的范圍是-10.89‰～5.89‰，均值為-7.37‰，地下水中穩(wěn)定同位素δD的范圍是-57.08‰～-28.99‰，均值為-39.45‰，地下水穩(wěn)定同位素δ18O的范圍是-8.94‰～-4.25‰，均值為-6.51‰，區(qū)域內(nèi)地表水與地下水流向基本一致。

從研究區(qū)氫氧同位素空間分布數(shù)據(jù)來看(圖6)，研究區(qū)上游(池淮溪—馬金溪段)地表水的δ18O隨著河流流向不斷增高，這是由于河流在運移過程中不斷受到蒸發(fā)作用，使得同位素富集造成的；流域中游(常山江段)有眾多支流的匯入，匯入的支流δ18O含量相對較低，造成了地表水的δ18O突變式降低，但仍呈現(xiàn)出地表水的δ18O值在沿程方向逐漸富集的規(guī)律[圖6(c)]。衢江流域上游和中下游的地下水氫氧同位素空間分布特征[圖6(b)、圖6(d)]表明研究區(qū)地表水和地下水的交換比較強(qiáng)烈，但是在其相互轉(zhuǎn)換過程中受到了一定的蒸發(fā)作用，此現(xiàn)象在中下游的紅層區(qū)也比較顯著，雖然紅層巖體具有各向異性等特征，但其同位素含量變化特征與地表水變化特征相似，這也在一定程度上說明此區(qū)域主要是地表水補給地下水，且兩者水力聯(lián)系比較密切。

圖6 研究區(qū)地表水、地下水δD和δ18O空間分布Fig.6 Spatial distribution of δD and δ18O of surface water and groundwater in the study area

3.3.2 氫氧穩(wěn)定同位素特征分析

水循環(huán)過程中，因為氫、氧同位素的平行分餾作用，大氣降水中的δD和δ18O之間具有一定的線性關(guān)系。本次采樣時間為9月份，處于研究區(qū)的枯水季，無法獲取降雨數(shù)據(jù)，故本文利用了雷明等[32]的大氣降水同位素數(shù)據(jù)，降水線方程為

δD=8.31δ18O+15.96

(5)

從衢江流域的大氣降水線可以看出：當(dāng)?shù)卮髿饨邓€(local mean water level，LMWL)斜率為8.31，比全球大氣降水線(global mean water wevel，GMWL)的斜率8大；大氣降水氘盈余值15.96‰，大于10‰。當(dāng)?shù)卮髿饨邓€的斜率與截距相對于全球大氣降水線的偏大，這與研究區(qū)的地理位置有關(guān)，研究區(qū)位于浙江省西部，距海僅200 km左右，氣候溫暖濕潤，其降水受太平洋季風(fēng)影響大，海洋水汽是研究區(qū)大氣降水的主要來源。

從研究區(qū)水樣δD-δ18O來看(圖7)，地下水和地表水同位素數(shù)據(jù)主要集中于當(dāng)?shù)卮髿饨邓€右下方，這說明其主要補給來源是大氣降水，且受到了一定的蒸發(fā)作用，且由于不同地勢地形的影響，蒸發(fā)程度也不一樣；部分上游地區(qū)地表水位于當(dāng)?shù)卮髿饨邓€的左上方，這可能由于上游山區(qū)中發(fā)生了水汽凝結(jié)過程，貧化了大氣降水中的同位素值。大部分氫氧同位素分布比較集中，但也有小部分點位的氫氧同位素組成差異表較大，這是由于地表水

圖7 氫氧同位素組成Fig.7 Hydrogen and oxygen isotopic composition

和地下水的氫氧同位素組成與其在流域中的空間位置有一定關(guān)系，在一定程度上反映了地表水與地下水的轉(zhuǎn)換關(guān)系。

同位素地表水線(surface water level，SWL)和同位素地下水線(ground water level，GWL)斜率基本一致，說明研究區(qū)地表水和地下水的交換比較強(qiáng)烈，但比當(dāng)?shù)卮髿饨邓€(LMWL)斜率低，這說明大氣降水補給地表水和地下水的過程中受到了一定的蒸發(fā)作用，使得同位素含量發(fā)生了變化。

3.4 地表水與地下水的轉(zhuǎn)化關(guān)系

依據(jù)式(2)和式(3) ，可以得出各測點上河水對地下水的貢獻(xiàn)百分比f如表1所示。

表1 研究區(qū)各河段地表水對地下水的轉(zhuǎn)化率Table 1 Conversion rate of surface water to groundwater in each reach of the study area

可見研究區(qū)上游池淮溪和馬金溪部分的貢獻(xiàn)百分比大于100%，說明此區(qū)域主要是地下水補給地表水，同理取其倒數(shù)可得池淮溪、馬金溪河段地下水對地表水的貢獻(xiàn)百分比為19.97%和19.36%；其他河段均為地表水補給地下水，其中江山港河段的地表水對地下水貢獻(xiàn)率最低，是由于此區(qū)域大型工廠繁多，對于地表水的需求量比較大，造成了此河段流量較小，地表水與地下水的轉(zhuǎn)換速率變慢。從表1中可知常山江河段地表水對地下水的貢獻(xiàn)率高達(dá)40.96%，主要原因是常山江由池淮溪和馬金溪交匯形成，而依據(jù)式(3)僅考慮了池淮溪河段，僅靠一種示蹤劑無法計算出其中一種混合來源的百分比，必須還要借助另外一種示蹤劑，即分析水中另外一種同位素或水化學(xué)離子等。本研究利用式(4)，借助氫同位素(δD)按照三相混合來計算，其中，Cc表示常山江上游池淮溪水樣中氫同位素值，Cm表示常山江上游馬金溪水樣中氫同位素值。將氫同位素值代入式(4)得出常山江實際地表水對地下水的轉(zhuǎn)化率為25.18%，明顯小于按兩項混合計算的40.96%。

利用地表水對地下水的轉(zhuǎn)化率，結(jié)合實測各河段斷面流量(表1中Qs)和式(3)、式(4)可進(jìn)一步得出衢江流域與周邊地下水的轉(zhuǎn)化量，計算結(jié)果表明：衢江流域上游地下水對地表水的平均貢獻(xiàn)率為19.7%、日均補給量為2.73×106m3/d；中下游地表水對地下水的平均貢獻(xiàn)率為22.77%、日均補給量為3.49×106m3/d。

4 結(jié)論

(1)TDS、pH指標(biāo)的空間分布規(guī)律顯示，研究區(qū)地表水與地下水變化特征相似，兩者轉(zhuǎn)換頻繁；流域上游主要是地下水補給地表水，中下游則主要是地表水補給地下水。

(2)陰陽離子組成顯示，研究區(qū)水化學(xué)類型主要是HCO3-Ca、HCO3·Cl-Ca類型水，屬于低礦化度水，說明流域內(nèi)水的最終來源為大氣降水；地表水與地下水的水化學(xué)類型基本相似，說明兩者之間存在一定的水利聯(lián)系。Gibbs圖表明研究區(qū)水中離子含量主要受巖石風(fēng)化影響，上游部分河段同時受降水因素影響。

(3)氫氧穩(wěn)定同位素空間分布規(guī)律和特征顯示，在少有支流匯入的上游，氫氧同位素含量是沿著河水流向不斷富集的；依據(jù)氫氧同位素與大氣降水線之間的關(guān)系，進(jìn)一步說明了研究區(qū)地下水和地表水均受大氣降水補給，且相互補給過程中均受蒸發(fā)作用的影響。

(4)依據(jù)δ18O在水體轉(zhuǎn)化中的質(zhì)量守恒方程，算出了地表水對地下水的貢獻(xiàn)率，定量反映了地表水與地下水的轉(zhuǎn)化關(guān)系，結(jié)果顯示，研究區(qū)上游主要是地下水補給地表水，地下水對地表水的平均貢獻(xiàn)率為19.67%，日均補給量為2.73×106m3/d；中下游主要是地表水補給地下水，地表水對地下水的平均貢獻(xiàn)率為22.77%，日均補給量為3.49×106m3/d。