陳學(xué)群，劉 丹，張文靜，管清花，田嬋娟，辛光明

（1.山東省水利科學(xué)研究院，濟(jì)南 250013；2.山東省水資源與水環(huán)境重點(diǎn)實(shí)驗(yàn)室，濟(jì)南 250013；3.吉林大學(xué)，長春 130021）

0 引 言

威海市屬資源型缺水城市，海水入侵問題突出，制約著威海市經(jīng)濟(jì)社會的可持續(xù)發(fā)展。為解決這一問題，威海市利用雨洪資源高效利用及地表地下聯(lián)合調(diào)控技術(shù)有效增加了地下水資源量，并防止海水入侵。地下水人工回灌技術(shù)是將回灌水源補(bǔ)給地下含水層的過程，在實(shí)現(xiàn)增加地下水資源量［1］、控制地面沉降［2，3］和海水入侵［4，5］等環(huán)境地質(zhì)問題的同時，可有效實(shí)現(xiàn)地表-地下水資源的科學(xué)調(diào)配和高效利用［6-8］。19世紀(jì)初期，英國修建第一個地下水補(bǔ)給工程將河湖淡水補(bǔ)給至格拉斯哥市［9］；此后，德國、瑞典、美國等國家開始大量開展地表水補(bǔ)給地下水的研究和實(shí)踐［10］。20 世紀(jì)60年代，上海市采取人工回灌技術(shù)防治地下水位下降引起的地面沉降等環(huán)境地質(zhì)問題［11］；隨后，我國山東、河北等地區(qū)開始大量的人工回灌技術(shù)研究。

隨著人工回灌技術(shù)的廣泛應(yīng)用，回灌水源多種多樣，包括雨洪資源［12，13］、地表水［14］、再生水［15，16］及自來水［17］等。由于回灌水源與目標(biāo)含水層地下水存在水質(zhì)差異，以及回灌過程中與含水層介質(zhì)之間存在復(fù)雜的水文地球化學(xué)過程［8，18，19］，對地下水環(huán)境質(zhì)量安全產(chǎn)生影響［20，21］。因此，查明不同水源人工回灌地下過程中地下水化學(xué)演化特征尤為重要。國外開展人工回灌過程中地下水水化學(xué)組分變化規(guī)律研究起步較早，澳大利亞Bolivar 地區(qū)利用同位素示蹤劑定量評價回灌水源與地下水的混合程度及發(fā)生的水文地球化學(xué)作用［22］；Herczeg 等［23］在澳大利亞某ASR試驗(yàn)場地以雨洪資源為回灌水源開展為期5年ASR試點(diǎn)研究，結(jié)果表明雨洪水回補(bǔ)灰?guī)r含水層過程中發(fā)生CO2和碳酸鹽礦物溶解作用及硫酸鹽礦物氧化作用，引起地下水中

Ca2+、Mg2+、HCO-、SO42-濃度變化；同時，大量研究表明人工回灌過程中介質(zhì)中的F、Pb、Ba、As 等元素出現(xiàn)不同程度的溶出［19，24，25］。我國在人工回灌技術(shù)對地下水環(huán)境影響方面的研究多集中在室內(nèi)物理實(shí)驗(yàn)和模擬研究，杜新強(qiáng)等［26］采用地球化學(xué)模擬方法，在上海市開展自來水人工回灌防治地面沉降中發(fā)現(xiàn)，地下水TDS降低，但Fe、Mn有機(jī)組分等濃度增大［27］。

目前，國內(nèi)外開展的人工回灌對地下水環(huán)境影響的研究多集中在雨洪資源、河水、湖水等常規(guī)水源的人工回灌工程，而再生水等非常規(guī)水源具有較大的潛力［28，29］?；诨毓嗨吹亩鄻有?，為了防止回灌水源中潛在微生物等污染物進(jìn)入地下水環(huán)境，回灌到目標(biāo)含水層前應(yīng)進(jìn)行消毒處理［30］。但是，已有研究無法有效評估氯化消毒水源回補(bǔ)地下水化學(xué)特征，尤其是場地尺度上針對氯化消毒水體回補(bǔ)地下水對地下水環(huán)境的影響研究相對較少，亟需開展進(jìn)一步研究。本文以氯化消毒的乳山河水作為回灌水源，選擇乳山河地下水庫為研究區(qū)，利用乳山河地下水庫人工回灌試驗(yàn)場地開展地下水環(huán)境要素、水化學(xué)特征動態(tài)監(jiān)測，分析氯化消毒水源回補(bǔ)地下水過程中回灌水源和回灌條件對地下水水化學(xué)特征的影響及伴隨的水文地球化學(xué)作用，為人工回灌工程的開展提供科學(xué)參考。

1 材料與方法

1.1 研究區(qū)概況

乳山河地下水庫位于乳山河干流與司馬莊河上，屬低山丘陵區(qū)，地勢由北向南臺階式下降（圖1）。地貌上屬于魯東剝蝕構(gòu)造為主的低山丘陵區(qū)。區(qū)域內(nèi)地表出露地層為第四系松散堆積地層，巖性一般以礫石、粗砂、中細(xì)砂、壤土為主。在入?？谔?，巖性主要為淤泥質(zhì)黏土、中細(xì)砂等。第四系孔隙潛水含水層組是本次研究的目標(biāo)含水層，本地地下水天然流場為由西北向東南，主要接受大氣降水補(bǔ)給、側(cè)向徑流及河道滲漏補(bǔ)給，主要排泄方式為人工開采。

圖1 地下水庫位置及工程布局圖Fig.1 Location and layout of study area

乳山河地下水庫工程分為上級庫工程與下級庫工程，上級庫工程包括上級庫地下截滲墻、于家莊橡膠壩工程及管理設(shè)施工程；下級庫工程包括下級庫地下截滲墻、乳山寨橡膠壩、鳳臺頂橡膠壩。

1.2 試驗(yàn)場地概況

試驗(yàn)場地位于位于乳山河地下水庫的上級庫，于家莊橡膠壩南側(cè)，乳山河?xùn)|側(cè)30 m。場地內(nèi)建有試驗(yàn)井16 眼（圖2），井深為9.0～11.2 m，揭穿第四系潛水含水層，成井位置的統(tǒng)一有效保證了監(jiān)測層位的一致性。對目標(biāo)含水層的含水介質(zhì)進(jìn)行粒徑級配等主要理化性質(zhì)測試，結(jié)果表明目標(biāo)含水層主要為細(xì)砂和中砂，含水介質(zhì)的主要理化性質(zhì)見表1。通過X 射線熒光光譜分析確定主要礦物成分為石英、鈉長石、鈣長石和鉀長石等，此外還有白云石、方解石和高嶺石。

表1 含水介質(zhì)的主要理化性質(zhì)Tab.1 The main physical and chemical properties of media

圖2 人工回灌試驗(yàn)場地試驗(yàn)井分布圖Fig.2 Distribution of test well in recharge site

1.3 回灌水源和原生地下水化學(xué)特征

此次回灌水源為氯化消毒的乳山河河水，即加入NaClO 溶液調(diào)節(jié)有效氯濃度至5 mg/L。河水、回灌水源、原生地下水的水環(huán)境要素見表2。河水、原生地下水pH 值分別為7.54、7.79，為弱堿性，電導(dǎo)率和TDS 值相近；地下水受好氧微生物作用影響使DO、Eh值較低，屬還原環(huán)境。回灌水源中NaClO 的水解作用形成OH-和Na+，使得pH 值增加至8.45，且電導(dǎo)率和TDS值增大，高于地下水；NaClO 的加入使好氧微生物的生長受到抑制，且空氣中DO 不斷進(jìn)入回灌水源中使DO 增加至8.5，氧化性增強(qiáng)。

表2 試驗(yàn)場地各水體環(huán)境要素特征Tab.2 Characteristics of water environment elements in the test site

試驗(yàn)場地各水體水化學(xué)指標(biāo)如表3 所示；并繪制試驗(yàn)場地各水體水化學(xué)性質(zhì)Piper 圖（圖3）。河水的水化學(xué)類型為Cl·SO4-Na 型水，氯化消毒后的回灌水源為Cl-Na 型水；原生地下水水化學(xué)類型為Cl·HCO3-Na型水。

圖3 不同水體水化學(xué)性質(zhì)Piper圖Fig.3 Piper diagram of hydrochemical properties of different water

表3 試驗(yàn)場地各水體水化學(xué)指標(biāo)Tab.3 The hydrochemical indexes in the test site

1.4 試驗(yàn)方案

2019年7月4-5日，在乳山河地下水庫試驗(yàn)場地開展單井回灌試驗(yàn)。試驗(yàn)以RZK-07 為回灌井，分別在河道方向（RZK 05、RZK 06、RZK08）、垂直河道方向（RZK-03、RZK-11、RZK-15）和地下水流方向（RZK- 04、RZK- 10、RZK-13）設(shè)3 個監(jiān)測剖面共9 個監(jiān)測井，利用流量為2 m3/h 的潛水泵將配置好的回灌水源經(jīng)回灌井連續(xù)注入到目標(biāo)含水層，注水歷時為170 min，總補(bǔ)水量為5.67 m3。受回灌時長的限制，本次選擇水化學(xué)特征變化明顯的RZK-03、RZK-04、RZK-06 和RZK-08 監(jiān)測井地下水特征進(jìn)行重點(diǎn)分析。

為分析氯化消毒水源回灌過程中地下水化學(xué)特征，以回灌起始時刻為起點(diǎn)，按照先密后疏原則，現(xiàn)場測定主要監(jiān)測井中地下水環(huán)境要素：pH、電導(dǎo)率、TDS、DO、Eh；此外，在回灌試驗(yàn)的前期（0 min）、中期（170 min）后期（1 980 min）取RZK-04監(jiān)測井地下水樣品至實(shí)驗(yàn)室測定主要陰、陽離子（Na+、K+、Ca2+、Mg2+、Cl-、SO42-、NO3-、HCO3-、CO32-）和重金屬（總Fe、Fe2+）的質(zhì)量濃度。

2 結(jié)果與分析

2.1 地下水環(huán)境要素變化特征

（1）電導(dǎo)率和TDS 變化特征。氯化回灌水源回灌過程中電導(dǎo)率和TDS隨時間變化曲線如圖4所示。原生地下水的電導(dǎo)率和TDS 值分別為646 μS/cm 和380 mg/L，回灌水源的電導(dǎo)率（1 180 μS/cm）和TDS值（447 mg/L）明顯高于原生地下水。隨著回灌水源的注入，各監(jiān)測井中電導(dǎo)率和TDS 值逐漸增大，電導(dǎo)率峰值分別為746 μS/cm（RZK-03）、884 μS/cm（RZK-04）、720 μS/cm（RZK-06）和760 μS/cm（RZK-08），TDS 峰值分別為439 mg/L（RZK-03）、520 mg/L（RZK-04）、423 mg/L（RZK-06）和428 mg/L（RZK-08），隨后逐漸減小至穩(wěn)定值，但RZK-03、RZK-04監(jiān)測井電導(dǎo)率和TDS 穩(wěn)定值約為720 μS/cm 和404 mg/L，仍高于原生地下水背景值。

圖4 電導(dǎo)率、TDS隨時間變化曲線Fig.4 Variation curve of conductivity and TDS with time

這一變化規(guī)律主要受混合作用的影響，本次利用RZK-04監(jiān)測井TDS 的變化特征分析回灌水源與地下水的混合程度隨時間的變化?；旌媳萬（t）是指回灌過程中t時刻通過保守離子確定水樣中回灌水源所占的比例，其計(jì)算公式為：

式中：Cm(t)為t時刻地下水中示蹤劑濃度，mg/L；Ca為原生地下水中溶質(zhì)濃度，mg/L；Cj為回灌水源中示蹤劑濃度，mg/L。

回灌水源與地下水的混合程度同樣呈增大-減小-趨于穩(wěn)定的趨勢（圖5），其中地下水主流向的RZK-04監(jiān)測井的混合程度較高，最大混合比為52.6%，混合作用發(fā)生的程度較高，RZK-03 監(jiān)測井次之（22.2%），RZK-06 和RZK-08 監(jiān)測井的最大混合比僅為16.2%和18.0%。RZK-03、RZK-04 監(jiān)測井混合比穩(wěn)定后分別為13.2%和15.7%，表明除混合作用外，還存在其他水-巖作用影響。

圖5 各監(jiān)測井地下水中回灌水源所占比例隨時間變化曲線Fig.5 Variation curve of the proportion of recharge water in the groundwater of each monitoring well with time

（2）氧化還原環(huán)境（DO 和Eh）。試驗(yàn)場地地下水處于還原環(huán)境，回灌水源中DO 和Eh 值明顯高于地下水，因此人工回灌過程中地下水氧化還原環(huán)境發(fā)生變化（圖6）。從DO 隨時間變化曲線可以看出，氯化消毒水源注入后地下水DO 逐漸增大，RZK-04監(jiān)測井在50 min 達(dá)到最大值5.9 mg/L，較TDS 達(dá)到峰值的時間明顯提前，試驗(yàn)進(jìn)行到80 min 時DO 值開始下降。與電導(dǎo)率和TDS 值變化規(guī)律對比發(fā)現(xiàn)，如果僅存在混合作用，利用DO 計(jì)算混合比為51.8%和62.9%，均低于TDS 計(jì)算的最大混合比，表明該過程中存在微生物降解作用消耗水中DO［31］。Eh 值可以反映地下水中溶質(zhì)表現(xiàn)出的氧化還原性，其變化規(guī)律與DO 基本一致，隨著氯化消毒水源的注入，各監(jiān)測井中地下水回灌水源所占比例逐漸增大，Eh 值逐漸上升，50 min 時達(dá)到最大值-46.4 mV。因此，該過程中以O(shè)2為電子受體的需氧菌和兼性厭氧菌的有氧代謝能力增強(qiáng)，厭氧微生物的生長受到一定程度的抑制。

圖6 DO、Eh值隨時間變化曲線Fig.6 Variation curve of DO and Eh with time

（3）pH 值。從圖7 可以看出，地下水pH 值隨著氯化消毒水源的注入呈上升趨勢，其中RZK-04 監(jiān)測井pH 值在210 min 達(dá)到最大值8.6，RZK-03、RZK-06 和RZK-08 最大值分別為8.4、8.3、8.3，隨后呈下降趨勢。這是由于氯化消毒回灌水源pH 值高于地下水pH 背景值，混合作用導(dǎo)致pH 增加。其中地下水流方向上的RZK-04監(jiān)測井pH 值上升趨勢最大，主要是由于地下水水流方向上地下水中回灌水源與地下水混合程度最大。此外，試驗(yàn)過程中地下水中pH 最高值高于回灌水源，這是由于回灌水源中NaClO 的水解作用影響地下水中HClO 和ClO-所占比例（式2），HClO 具有強(qiáng)氧化性，與地下水中有機(jī)物等反應(yīng)，NaClO 的水解作用向右進(jìn)行，使得地下水中pH 值進(jìn)一步升高。pH 達(dá)到峰值的時間出現(xiàn)滯后效應(yīng)，且試驗(yàn)后期出現(xiàn)pH 值低于原生地下水pH 值的趨勢，主要是由于地下水氧化性增強(qiáng)，黃鐵礦（FsS2）被氧化形成H+（式2），導(dǎo)致pH值下降［32］。

圖7 pH值隨時間變化曲線Fig.7 Variation curve of pH with time

2.2 回灌水源與地下水的混合作用分析

根據(jù)RZK-04監(jiān)測井地下水和回灌水源樣品中各化學(xué)組分隨時間變化特征繪制半對數(shù)Schoeller 圖（圖8）。從圖8 中可以看出，原生地下水（回灌前）化學(xué)組分與回灌水源有明顯差異，在回灌中期（170 min）時，RZK-04 監(jiān)測井地下水的大部分組分（除HCO3-外）微小差別外，與回灌水源具有相似的趨勢，表明地下水化學(xué)變化特征受到了混合作用的影響，向回灌水源和原生地下水體中間移動?；毓噙^程中HCO3-濃度與地下水相近，回灌后高于地下水中HCO3-濃度，說明回灌水源中余氯發(fā)生水解作用形成HCO3-是主要原因，且該過程不是瞬時完成的。此外，回灌水源注入后陽離子濃度增大，其中Na+濃度高于回灌水源，表明該過程還存在著水-巖作用。

圖8 RZK-04監(jiān)測井地下水Schoeller圖Fig.8 Schoeller map of groundwater in RZK-04 monitoring well

2.3 回灌過程中的水-巖作用分析

回灌水源的注入除了與周圍地下水發(fā)生混合作用外，還受到溶液-含水層礦物相互作用的影響。根據(jù)Cl-質(zhì)量平衡確定地下水與回灌水源的混合程度［33］，在此基礎(chǔ)上計(jì)算其他離子的凈增量，分析氯化消毒水源回灌過程中礦物-溶液的相互作用對各離子濃度的影響，計(jì)算結(jié)果如表4所示。

表4 地下水中各離子凈增量計(jì)算結(jié)果 mg/LTab.4 The net increment of each ion in groundwater

在回灌中期（170 min），RZK-04監(jiān)測井地下水中Na+濃度從60 mg/L 增加至189 mg/L，且地下水中Na+濃度高于回灌水源，Net 為61.6 mg/L。含水層介質(zhì)中的含Na 礦物主要為鈉長石，屬硅酸鹽礦物，在本次試驗(yàn)的時間尺度內(nèi)溶解作用較弱，因此除混合作用外，還受回灌水源通過含水層介質(zhì)時的陽離子交換作用影響。

由于回灌水源與地下水中陽離子濃度差大小順序?yàn)椋篊a2+＞Na+＞Mg2+＞K+，Vandenbohede 等人研究表明高濃度差促進(jìn)Ca2+向吸附方向發(fā)生陽離子交換作用［34］，但Net 值僅為11.62 mg/L，這是由于含水層介質(zhì)中的方解石溶解作用，Ca2+濃度的增大伴隨著HCO3-濃度的升高，HCO3-的Net 為11.45，且試驗(yàn)結(jié)束時濃度高于原生地下水，進(jìn)一步說明CaCO3溶解作用的存在。

此外，K+濃度降低，Net為-0.62 mg/L 也進(jìn)一步證明了陽離子交換作用的存在。Ca2+和Mg2+濃度接近回灌水源中的濃度值，Net分別為11.62 和2.64 mg/L，這是由于在回灌過程中，回灌水源通過含水層介質(zhì)時發(fā)生可溶性礦物的溶解作用使得濃度升高。與已有研究結(jié)果不同［35］，SO42-濃度最終小于原生地下水和回灌水中SO42-濃度，這是由于氯化消毒水源的注入改變了地下水中原有的微生物群落結(jié)果，SO42-等作為電子受體參與降解作用［36］。

2.4 典型金屬濃度變化特征（Fe、Mn）

由于含水層中Fe、Mn 豐度較高，本次對RZK-04 監(jiān)測井地下水中總Fe、Fe2+和總Mn 的濃度變化進(jìn)行分析，濃度變化曲線如圖9 所示。回灌水源與原生地下水中總Fe 濃度均為0.03 mg/L，試驗(yàn)中期（170 min）地下水中總Fe濃度增大至0.062 mg/L，同時地下水中Fe2+濃度由0 mg/L增大至0.014 mg/L，這是由于含鐵礦物的溶出，根據(jù)回灌過程中DO 濃度的增大和pH 值低于原生地下水的變化趨勢，進(jìn)一步說明了該過程中介質(zhì)中FeS2氧化作用的存在［式（2）］，使得總Fe、Fe2+濃度增大，pH值降低。但Fe2+濃度的增加值低于總Fe的降低值，說明回灌水源注入后氧化還原電位升高使得Fe2+轉(zhuǎn)化成Fe3+。地下水中總Mn 濃度為0.011 mg/L，而回灌水源為0.004 mg/L。試驗(yàn)中期（170 min）地下水中總Mn濃度下降至0.004 mg/L，試驗(yàn)后期（1 980 min）出現(xiàn)增大趨勢（0.08 mg/L），說明存在硫錳礦（MnS）的溶出作用［式（5）］，但總Mn 濃度的增大存在滯后效應(yīng)，這是由于含錳礦物的溶出需要更多的反應(yīng)時間［37］。

圖9 RZK-04監(jiān)測井重金屬濃度變化曲線Fig.9 Heavy metal concentration curve of RZK-04 monitoring well

已有研究表明，地下含水層氧化性增強(qiáng)條件下易發(fā)生總Fe、總Mn 的溶出［35］。氯化消毒的回灌水源對好氧微生物的影響使其氧化性增強(qiáng)，因此回灌水源進(jìn)行氯化消毒后地下水中Fe、Mn元素的濃度應(yīng)重點(diǎn)關(guān)注。

3 結(jié) 論

試驗(yàn)場地開展的人工回灌試驗(yàn)，通過對氯化消毒水源回灌地下水過程中地下水環(huán)境要素（電導(dǎo)率、TDS、DO、Eh和pH）、主要陰、陽離子及重金屬（總Fe、總Mn）等指標(biāo)的測試分析，得出以下主要結(jié)論。

（1）回灌水源注入后，與原生地下水發(fā)生混合作用，地下水主流方向受回灌作用影響較大，大部分組分與回灌水源具有相似的趨勢，最大混合比為52.6%；

（2）存在水-巖相互作用，K+、Ca2+、Mg2+的Net 值分別為-0.78、11.62、2.64 mg/L，即方解石的溶解作用使Ca2+、Mg2+濃度增大；次氯酸鈉的加入易于發(fā)生陽離子交換作用，使K+濃度降低；

（3）從試驗(yàn)結(jié)果來看，氯化消毒水源人工回灌過程中地下環(huán)境氧化性增強(qiáng)，促進(jìn)含F(xiàn)e、Mn 礦物溶出，同時地下水pH 值降低；

（4）與常規(guī)水源相比，氯化消毒水源注入后影響地下水微生物群落結(jié)構(gòu)，使SO42-等離子作為電子受體被消耗，濃度降低。