1 引言
同位素分析技術(shù)已成為現(xiàn)代水文地質(zhì)學(xué)的研究方法之一，并已經(jīng)大量成功地應(yīng)用于地下水循環(huán)演化規(guī)律研究.例如，應(yīng)用環(huán)境同位素在確定地下水來源及組成究地下水補(bǔ)徑排條件，揭示含水層之間的水力聯(lián)系和地表水與地下水的相互作用，以及示蹤地下水運(yùn)動(dòng)等諸多方面具有重要作用.
地下水的同位素組成主要取決于地下水的來源，不同來源水的同位素組成的差異及其含量的時(shí)空變化是應(yīng)用環(huán)境同位素方法解決地下水補(bǔ)給循環(huán)特征的基礎(chǔ).研究地下水循環(huán)規(guī)律時(shí)常用到的有穩(wěn)定同位素(δ2H、δ18O、δ13C)和放射性同位素(3H、14C)等.地下水是古氣候信息的載體，應(yīng)用地下水的同位素資料可以研究古氣溫的變化，計(jì)算地下水的形成年齡，確定地下水的補(bǔ)給期，反映區(qū)域地下水流動(dòng)模式，以及對(duì)人類干擾的響應(yīng).
江蘇沿海平原東瀕黃海，《江蘇沿海地區(qū)發(fā)展規(guī)劃》于2009年被提升為國(guó)家戰(zhàn)略，該地區(qū)由此成為中國(guó)東部重要的經(jīng)濟(jì)增長(zhǎng)極，但也是江蘇省域水資源量最欠缺的地區(qū)，水質(zhì)優(yōu)良的深層地下水是本區(qū)經(jīng)濟(jì)社會(huì)發(fā)展的重要支柱.在經(jīng)濟(jì)發(fā)達(dá)的沿海地帶，由于缺乏對(duì)地下水形成演化規(guī)律的正確認(rèn)識(shí)，過量抽取地下水，破壞了濱海地帶地下水的循環(huán)過程，導(dǎo)致江蘇沿海平原地下水資源衰減、地質(zhì)環(huán)境惡化.根據(jù)2012年起開展的水文地質(zhì)調(diào)查資料表明，江蘇沿海平原區(qū)主采層第Ⅲ、IV承壓水出現(xiàn)地下水位不斷降低、水質(zhì)咸化、礦化度增高，以及誘發(fā)地面沉降與海水入侵等系列問題.目前，人們對(duì)高強(qiáng)度開發(fā)地下水條件下的沿海地區(qū)的水文演化問題缺乏系統(tǒng)認(rèn)識(shí)，由于地下水開采引發(fā)的地質(zhì)環(huán)境問題亟待解決：高強(qiáng)度開發(fā)地下水條件下水資源量和水質(zhì)演變的過程與機(jī)理如何?高強(qiáng)度開發(fā)地下水導(dǎo)致含水層系統(tǒng)結(jié)構(gòu)發(fā)生變異的機(jī)理與趨勢(shì)如何?區(qū)域地下水系統(tǒng)如何響應(yīng)區(qū)域水循環(huán)的變化?
20 世紀(jì)60年代至今，針對(duì)江蘇省域開展的大量水文地質(zhì)研究主要集中在長(zhǎng)江南翼，在北翼的江蘇沿海平原研究薄弱.例如，李從先等(2000;2009)、哈承佑等(1990)先后對(duì)區(qū)域古地理背景和水文地質(zhì)條件作過總結(jié)劃分.應(yīng)用地下水的環(huán)境同位素(δ2D、δ18O、3H、14C)研究了地下水的循環(huán)過程及其形成年代.基于同位素δ2D、δ18O和水化學(xué)系數(shù)分析，開展過南通地區(qū)地下水補(bǔ)給源、水化學(xué)變化特征研究.針對(duì)深層地下水開發(fā)利用引發(fā)的系列環(huán)境地質(zhì)問題，也有過簡(jiǎn)單的成因分析.針對(duì)江蘇沿海平原開展的地下水演化循環(huán)過程研究工作，主要集中在淺層地下水，對(duì)深層地下水水文地質(zhì)研究較為粗略.同時(shí)，深層地下水采集的同位素樣品種類數(shù)量少，以往研究中深井水樣大多為上下多層含水層的混合水，這就使得同位素及化學(xué)成分混雜不清，而無法進(jìn)行細(xì)致的示蹤研究.
江蘇沿海平原是省域水資源脆弱帶，正確認(rèn)識(shí)區(qū)域地下水水文化學(xué)演變過程是合理持續(xù)利用地下水資源的基礎(chǔ).由于影響地下水同位素組成的因素很多，單一的同位素提供的信息有其局限性，聯(lián)合運(yùn)用多種環(huán)境同位素提供的信息可相互應(yīng)證，更能展現(xiàn)地下水循環(huán)演化過程的全貌.因此，本研究利用近期獲得的鉆孔剖面地層對(duì)比數(shù)據(jù)資料，系統(tǒng)地對(duì)水文地質(zhì)鉆孔(深孔)進(jìn)行取樣分析，以期準(zhǔn)確標(biāo)定研究區(qū)深層地下水含水層的環(huán)境同位素的指紋特征.同時(shí)，通過對(duì)區(qū)域地下水穩(wěn)定同位素(δ2D、δ18O)和放射性同位素(3H、14C)組成的研究，探索江蘇沿海平原地區(qū)地下水補(bǔ)徑排條件、超采地下水引起的流場(chǎng)變化，以便為濱海含水層地下水資源合理開發(fā)利用、污染治理提供基礎(chǔ)資料和理論依據(jù).
2 材料與方法
2.1 研究區(qū)概況
江蘇省洋口港區(qū)位于江蘇省的東部沿海，隸屬江蘇省南通市如東縣(圖 1)，東瀕黃海，地勢(shì)平坦，地面標(biāo)高一般為1~5 m，地面微向東傾.研究區(qū)氣候溫和濕潤(rùn)，屬濕潤(rùn)的亞熱帶季風(fēng)氣候區(qū)，四季分明，雨熱同季，無霜期長(zhǎng)，具典型的海洋性氣候.年均氣溫14.8 ℃，一月份平均氣溫2.1 ℃，最低氣溫-10.6 ℃;七月份平均氣溫27.3 ℃，最高氣溫38.6 ℃.年平均降雨量1025 mm.
研究區(qū)位于揚(yáng)子陸塊下?lián)P子地塊北東部，地表為第四紀(jì)松散沉積物覆蓋.發(fā)育地層有上古生界中志留統(tǒng)-上泥盆統(tǒng)(S2-D3)、石炭系-二疊系(C-P)、中生界白堊系浦口組(K2p)、新生界古近系泰州組(E1t)、阜寧組(E1f)、三垛組(E2-3s)和新近系鹽城組(N1-2yc).
根據(jù)地下水賦存介質(zhì)，松散巖類孔隙含水巖組為研究區(qū)主要含水巖組.根據(jù)地層的時(shí)代將含水層組由上到下分別為潛水、第Ⅰ、Ⅱ、Ⅲ、Ⅳ承壓含水層組.考慮到研究區(qū)的潛水含水層和第Ⅰ、Ⅱ承壓含水層易于受到降水補(bǔ)給，水力聯(lián)系緊密，地下水徑流速度快，更新能力較強(qiáng)，將潛水含水層和Ⅰ、Ⅱ承壓含水層歸入淺層孔隙水含水系統(tǒng).第 Ⅲ、Ⅳ承壓含水層是研究區(qū)地下水主采層，地下水徑流速度慢，更新能力弱，將Ⅲ、Ⅳ承壓含水層歸入深層孔隙水含水系統(tǒng).
2.2 采樣測(cè)試分析
本研究采樣點(diǎn)大致沿地下水流向(自西向東)分布，針對(duì)本區(qū)地下水分層專門監(jiān)測(cè)井進(jìn)行取樣分析，分別采集地表水、潛水、承壓水，以第Ⅲ、Ⅳ承壓水為重點(diǎn).于2014年10—11月采集深層地下水21件，其中，第Ⅲ承壓水同位素(δ2D、δ18O、3H、14C)14組，第Ⅳ承壓水7組;采集淺層地下水同位素(δ2D、δ18O)樣品7件;采集地表水同位素(δ2D、δ18O)樣品5件，其中，如泰運(yùn)河河水樣品3件，入海口海水樣品2件.具體采樣位置見圖 1.
?圖1 江蘇洋口港研究區(qū)位置及采樣點(diǎn)分布圖
氫氧同位素分析：樣品水分別采集在250 mL聚乙烯塑料瓶中，確保瓶?jī)?nèi)無氣泡后，用封口膠密封，常溫避光保存，送至國(guó)土資源部南京地質(zhì)礦產(chǎn)研究所實(shí)驗(yàn)室測(cè)試.用鋅還原法制備氫氣，供質(zhì)譜氫同位素分析;用CO2-H2O平衡法制備供氧同位素分析的CO2氣體.氫、氧同位素分析是在 MAT 253質(zhì)譜計(jì)上進(jìn)行，相對(duì)V-SMOW標(biāo)準(zhǔn).氫同位素分析精度為±2‰(1σ)，氧同位素分析精度為±0.1‰(1σ).放射性同位素3H樣品采集在1250 mL聚乙烯塑料瓶中，確保瓶?jī)?nèi)無氣泡后，用封口膠密封，常溫避光保存，分析測(cè)試由中國(guó)地質(zhì)科學(xué)院巖溶地質(zhì)研究所實(shí)驗(yàn)室完成.利用Qnantulus1220 超低本液體閃爍譜儀測(cè)定氚含量.
14C測(cè)定：放射性同位素14C樣品采集后利用沉淀桶制備水樣，沉淀液存放在5 L聚乙烯塑料桶中，用封口膠密封，常溫避光保存，分析測(cè)試由中國(guó)地質(zhì)科學(xué)院巖溶地質(zhì)研究所實(shí)驗(yàn)室完成.標(biāo)準(zhǔn)采用中國(guó)糖碳，14C半衰期采用5730 a，起始計(jì)年為 1950年，液閃儀(Quantulus-1220(LKB)測(cè)定.
18O、2H樣品測(cè)試都在2014年前完成，3H、14C測(cè)定樣品測(cè)試都在2015年5月前完成.為示蹤研究地下水的來源及建立地下水的循環(huán)模式，還收集了本地區(qū)相關(guān)的大氣降水的同位素資料.
統(tǒng)計(jì)分析以SPSS19.0為平臺(tái)進(jìn)行相關(guān)分析、曲線擬合.趨勢(shì)線方程根據(jù)測(cè)試數(shù)據(jù)進(jìn)行曲線擬合，評(píng)價(jià)擬合優(yōu)度，求得線性方程.圖形以Arcgis 10、Coreldraw X4 為平臺(tái)制作.
3 測(cè)試結(jié)果及分析
3.1 氫氧同位素組成
3.1.1 研究區(qū)參考雨水線
根據(jù)全球大氣降水監(jiān)測(cè)數(shù)據(jù)庫(kù)(GNIP)參考站點(diǎn)資料，綜合考慮參考站點(diǎn)與研究區(qū)地理位置、氣候分區(qū)、地形地貌的相似性，選擇南京站作為本區(qū)降水同位素的參考站點(diǎn).根據(jù)GNIP南京站點(diǎn)1987—1992年的58個(gè)降水同位素?cái)?shù)據(jù)，大氣降水穩(wěn)定同位素δ18O值介于-11.83‰~-0.09‰之間，平均為-7.36‰，δD介于-83.5‰~17.9‰之間，平均為-44.89‰.根據(jù)南京站降水同位素的監(jiān)測(cè)數(shù)據(jù)可以得到本區(qū)雨水線方程為: δD=8.49δ18O+17.71(N=58，R2=0.967，p<0.01)，與全球大氣雨水線(δD=8δ18O+10)、中國(guó)的雨水線方程(δD=7.9δ18O+8.2)相比，其斜率和截距偏高(圖 2).雨水氫氧同位素主要由氣候特點(diǎn)和地理位置決定的.南京地區(qū)離海洋較近，其降水主要受太平洋季風(fēng)的影響，海洋水汽是其降水的主要來源.南京地區(qū)雨水線斜率、截距偏高反映了區(qū)域濕潤(rùn)多雨的氣候特點(diǎn).
?圖2 研究區(qū)δD-δ18O分布圖
3.1.2 地表水、地下水的氫氧同位素組成
研究區(qū)采集的地表水、淺層地下水、深層地下水樣品的δ18O、δD、3H、14C值見表 1.
表1 江蘇洋口港地區(qū)地表水地下水δD、δ18O、3H、14C同位素測(cè)試數(shù)據(jù)表
河水是地下水的重要補(bǔ)給來源之一，討論地表水的變化特征及形成規(guī)律是進(jìn)一步開展本區(qū)地下水形成演化規(guī)律研究的前提.如泰運(yùn)河河水從馬塘鎮(zhèn)到入海口δD、δ18O的變化范圍分別為-50.45‰~-10.42‰、-7.18‰~-1.15‰，均值分別為-34.91‰、-4.57‰.地表水氫氧同位素空間差異明顯，表現(xiàn)出氫氧同位素沿河流流向富集.在 δD-δ18O關(guān)系圖中，河水的同位素組成均落在雨水線的右下方(距離雨水線較近)，并位于同一條直線上，說明河水沿著徑流方向受到海水影響而增大，并受到了一定的蒸發(fā)作用影響.
研究區(qū)地下水δ18O值介于-7.22‰~-3.25‰，均值為-5.43‰，δD值介于-55.31‰~-28.15‰，均值為-41.97‰.淺層地下水樣品的δ18O值介于-6.39‰~-3.91‰，均值為-5.03‰，δD值介于-42.83‰~-30.69‰，均值為-36.17‰;深層地下水樣品的δ18O值介于-7.22‰~-3.25‰，均值為-5.57‰，δD值介于-55.31‰~-28.15‰，均值為-43.90‰.淺層地下水沿徑流方向總體表現(xiàn)為沿程富集，淺層地下水樣品的δ18O、δD值普遍比深層地下水偏正.
結(jié)合采樣點(diǎn)的數(shù)據(jù)特征(表 1)和圖 2可見，地下水中的δD值和δ18O值及兩者之間的關(guān)系呈以下特征：①δD、δ18O值的變化區(qū)間均比大氣降水的相應(yīng)值小得多;②δD、δ18O的均值均比大氣降水的相應(yīng)值大;③表示深層地下水的點(diǎn)大都分布在當(dāng)?shù)亟邓€附近并靠近左下角，顯示古氣候效應(yīng);④潛水、第 Ⅲ、Ⅳ承壓水沿著主徑流方向由西至東，δD和δ18O值總體上增大;⑤在垂直方向上，隨著深度增加，δ18O和δD值逐漸降低.
以上諸特征聯(lián)合表明：①地下水主要來源于大氣降水的入滲補(bǔ)給;②大氣降水在補(bǔ)給地下水過程中，經(jīng)歷了明顯的蒸發(fā)過程，淺層地下水與深層地下水在補(bǔ)給過程中經(jīng)過不同程度的蒸發(fā)影響;③淺層地下水與深層地下水補(bǔ)給來源可能存在溫度差異，深層地下水的補(bǔ)給來源溫度較低;④沿徑流方向，補(bǔ)給地下水蒸發(fā)作用越趨強(qiáng)烈;⑤深、淺層地下水之間同位素組成的明顯差異表明深、淺層地下水之間的相互影響較微弱;⑥可看出同一地點(diǎn)各層位的氫氧同位素?cái)?shù)據(jù)具有明顯的呈層性.
3.2 深層地下水氚分布特征及其指示意義
3.2.1 研究區(qū)大氣降水氚濃度的恢復(fù)
氚的半衰期為12.43 a，在水中以HTO形式存在，成為天然水的一部分而參加水循環(huán)，因此，成為追蹤各種水文地質(zhì)作用的一種理想示蹤劑，特別是它的放射性計(jì)時(shí)性，使之成為測(cè)定地下水年齡的重要技術(shù)手段，并對(duì)于研究大氣降水入滲、現(xiàn)代滲入起源、地下水補(bǔ)給、賦存及運(yùn)動(dòng)具有重要意義.利用氚研究地下水的運(yùn)動(dòng)規(guī)律時(shí)，必然會(huì)遇到氚的輸入背景值問題.大氣降水中的氚有兩個(gè)來源：一是來源于大氣層上層宇宙射線形成的中子和氮原子的相互作用而形成的宇宙氚;二是來源于1952 年以來的熱核爆炸，1963年達(dá)到高峰.隨后由于國(guó)際公約對(duì)核試驗(yàn)的限制，降水中氚濃度以指數(shù)形式遞減，21世紀(jì)初已恢復(fù)到自然水平.為查明大氣降水氚濃度的時(shí)空分布規(guī)律，國(guó)際原子能機(jī)構(gòu)()和世界氣象組織(WMO)在世界各地建立了觀測(cè)站，而我國(guó)的背景值監(jiān)測(cè)工作起步較晚，缺少1953—1978年20多年的系列觀測(cè)資料.由于實(shí)測(cè)資料不能滿足應(yīng)用要求，因而需應(yīng)用合適的數(shù)學(xué)方法對(duì)大氣降水的氚濃度進(jìn)行恢復(fù).目前應(yīng)用的恢復(fù)方法主要有關(guān)秉鈞法、插值法、雙參考曲線法、人工神經(jīng)網(wǎng)絡(luò)法、因子分析法等.本次研究基于應(yīng)用因子分析法，建立了南京地區(qū)的全球大氣降水氚恢復(fù)曲線(圖 3)，以提供研究區(qū)降水氚的背景資料.
?圖3 南京大氣降水氚濃度 1960—2010 年恢復(fù)曲線
3.2.2 研究區(qū)深層地下水氚濃度分析
根據(jù)南京地區(qū)降水氚恢復(fù)曲線，1963年南京地區(qū)大氣降水中的氚濃度可達(dá)1746 TU，60年代一般在200 TU以上，70年代一般在60~150 TU之間，80年代逐漸恢復(fù)到自然水平.根據(jù)不同時(shí)期形成的地下水的平均氚值，根據(jù)衰變定律，可推測(cè)到地下水各水體中氚的演化路徑.
江蘇洋口港地區(qū)深層承壓地下水中大部分氚含量為<2 TU(表 1)，均為低氚水，表明洋口港地區(qū)深層承壓地下水形成于地質(zhì)歷史時(shí)期.研究區(qū)部分深井受到人為建井而引起上下含水層貫通進(jìn)而引起越流補(bǔ)給而出現(xiàn)含氚水.洋口港南部三余鎮(zhèn)沿海地區(qū)第III承壓水氚含量為15.80 TU，根據(jù)衰變定律可以推斷，地下水混入了20世紀(jì)60年代降水.如東縣豐利鎮(zhèn)、掘港鎮(zhèn)部分深井氚含量為5~9 TU，可以推斷，混入了60—70年代初期降水.如東縣馬塘鎮(zhèn)是研究區(qū)水位降落中心，所測(cè)深井第III承壓水氚含量超過12 TU，可以推斷水位降落中心混入了60年代高氚含量降水.部分深井地下水出現(xiàn)高氚含量，應(yīng)該視為超采地下水誘發(fā)的現(xiàn)代水補(bǔ)給.由于超采深層地下水，深層地下水水位遠(yuǎn)低于淺層地下水水位，形成了較大的水位差，大氣降水與淺層水參與到深層水的循環(huán)中來，淺層水與深層水之間的水力聯(lián)系加強(qiáng)，水循環(huán)速率與氚值增大.
3.3 深層地下水放射性14C年齡分析
研究區(qū)深層承壓地下水的14C年齡的變化范圍為7410~26140 a，均值為18330 a(n=15)，其中，第III承壓水14C年齡的變化范圍為7410~26140 a，均值為17805 a;第Ⅳ承壓水14C年齡的變化范圍為16140~25890 a，均值為19775 a(n=4)(表 1).14C年齡在垂直方向上由淺部至深部逐漸增大;沿著地下水的流向自西向東，含水層埋深越大，地下水交替越差，14C年齡總體上逐漸增大.
分析發(fā)現(xiàn)，14C年齡與δ18O顯著負(fù)相關(guān)，相關(guān)系數(shù)為-0.901(n=15，p<0.01)，研究區(qū)地下水形成年齡越老，δ18O愈加偏負(fù)(圖 4)，這與區(qū)域第四紀(jì)氣候演變有著密切的關(guān)系.距今13000~60000 a 的玉木冰期，δ18O 值比現(xiàn)在平均低約10‰，我國(guó)東部沿海地區(qū)地下水樣品在全新世和晚更新世不同時(shí)期降水中δ18O 差值也大概在 2‰左右.據(jù)此可確定，研究區(qū)大部分深井地下水為末次冰期(大理冰期)的降水補(bǔ)給，只有極少部分深井為全新世早期的降水補(bǔ)給.
?圖4 江蘇洋口港地區(qū)δ18O與14C年齡關(guān)系分析
研究區(qū)豐利鎮(zhèn)第III承壓水為咸水，14C測(cè)得的年齡小于10000 a，可能是全新世高溫期大陸鹽化作用的產(chǎn)物.
3.4 地下水來源演化研究
3.4.1 地下水補(bǔ)給來源分析
瑞利蒸餾過程影響著降水中的同位素組成，水在蒸發(fā)過程中重同位素傾向于保留在相對(duì)于水蒸氣分子來說活性較低的液態(tài)水中，使殘留溶液的δD、δ18O值偏正.在蒸發(fā)濃縮作用下，研究區(qū)淺層地下水礦化度普遍要比深層地下水高.淺層地下水主要來源于現(xiàn)代大氣降水的入滲補(bǔ)給，其蒸發(fā)作用要比深層地下水強(qiáng)烈得多，淺層地下水中的重同位素在較強(qiáng)的蒸發(fā)作用下相對(duì)富集，因此，δD、δ18O值偏正，而深層地下水中δD、δ18O值相對(duì)來說則偏負(fù).沿海一帶淺層地下水δD、δ18O值普遍比中西部地區(qū)偏正，是由于沿海地區(qū)地下水受到一定程度δD、δ18O值偏正的現(xiàn)代海水影響.
地下水趨勢(shì)線(蒸發(fā)線)與當(dāng)?shù)赜晁€的交點(diǎn)可以代表補(bǔ)給地下水時(shí)的源區(qū)降水氘氧同位素的平均組成.代表深層水補(bǔ)給源的B點(diǎn)位于淺層水補(bǔ)給源A點(diǎn)之下(圖 2)，表明深層地下水接受補(bǔ)給的降水δ18O、δD 值比淺層水偏負(fù)，推測(cè)深層地下水是在氣溫較低的條件下接受的降水補(bǔ)給.
大氣降水的年平均δ18O、δD 值與該區(qū)的地面平均氣溫正相關(guān).中國(guó)的地面平均氣溫與現(xiàn)代大氣降水的δ值建立了如下關(guān)系式：
式中，t為地面平均氣溫(℃).
根據(jù)本次采樣的地下水氫氧穩(wěn)定同位素資料，對(duì)深層含水組進(jìn)行統(tǒng)計(jì)，根據(jù)求得的δ18O-δD關(guān)系式，通過求解方程從而獲得相應(yīng)補(bǔ)給時(shí)期大氣降水的δ18O和δD值.然后代入式(1)和(2)，求出地下水補(bǔ)給期的地面平均氣溫.
根據(jù)氫氧穩(wěn)定同位素資料，求得的淺層地下水、深層地下水的δ18O-δD關(guān)系式為：
δD=4.69δ18O-12.61(N=7，R2=0.928，p<0.01)
δD=6.19δ18O-9.42(N=21，R2=0.892，p<0.01)
通過求解方程從而獲得深層地下相應(yīng)補(bǔ)給時(shí)期大氣降水的δ18O和δD值分別為-11.80‰和-82.46‰.代入公式(1)、(2)求得深層地下水補(bǔ)給期的地面平均氣溫介于3.4~4.1 ℃，比現(xiàn)代氣溫(14.8℃)低約11℃.
根據(jù)放射性同位素14C測(cè)年(表 1)，研究區(qū)深層含水系統(tǒng)的大部分地下水形成于15000~26000 a，處于晚更新世末期，該時(shí)期氣溫較低，降水δ18O、δD 值偏負(fù).據(jù)此類比，可以判斷研究區(qū)深層地下水的補(bǔ)給源形成于末次冰期(大理冰期)盛期較寒冷的大氣降水淋濾補(bǔ)給.
3.4.2 地下水循環(huán)特征與可更新性
處于水循環(huán)系統(tǒng)中不同的水體，因成因不同而具富集程度不同的氫氧同位素.不同水體中的同位素濃度變化可示蹤其形成和運(yùn)移方式，認(rèn)識(shí)變化環(huán)境下的水循環(huán)規(guī)律及水體間的相互關(guān)系.從研究區(qū)地表水、淺層地下水、第III承壓水、第IV承壓水的δD、δ18O平均值來看，發(fā)現(xiàn)地表水的同位素均值最大、淺層地下水次之、深層地下水最小(表 2);沿著地表水-淺層地下水-深層地下水空間梯度，隨著地下水埋藏深度增加，水體中的δD、δ18O值總體上呈下降趨勢(shì).地表水、淺層地下水沿徑流方向都表現(xiàn)為沿程富集，而且氫氧均值差異較小，表明它們之間的水體聯(lián)系較為緊密，淺層地下水以垂向流動(dòng)為主.深、淺層地下水及第III承壓水、第IV承壓水之間同位素組成的明顯差異，表明深層地下水各水體之間的相互影響較微弱.
?表2 江蘇洋口港地區(qū)各水體氫氧碳同位素分布
地下水δD、δ18O 值結(jié)合地下水埋深深度可以有效地標(biāo)記地下水的徑流條件與更新狀況.第Ⅲ承壓水是研究區(qū)主要開采層位，根據(jù)水位監(jiān)測(cè)資料，由于多年開采第Ⅲ承壓水，形成了以馬塘鎮(zhèn)為中心的水位降落漏斗，水位下降區(qū)范圍達(dá)698 km2，其中心水位埋深超過40 m(圖 5).
??圖5 江蘇洋口港地區(qū)第Ⅲ承壓水水位埋深空間分布
本次研究利用Arcgis的DEM模型作的δD、δ18O空間分布圖(圖 6)顯示，研究區(qū)第Ⅲ承壓水氘氧同位素的分布區(qū)域大致與降落漏斗分布一致.反映出隨著地下水位下降，開采的深部年老水且δD、δ18O值相對(duì)較小水的比例增加.研究區(qū)地形平坦，水力坡度非常小;深層含水層之間發(fā)育滲透系數(shù)較小的粘土類地層，因此，深層承壓地下水基本保持相對(duì)封閉狀態(tài)，以側(cè)向水平徑流為主，徑流滯緩，總體上處于封閉-半封閉狀態(tài)，可更新能力弱，具有不可再生資源的屬性.
?圖6 江蘇洋口港地區(qū)第Ⅲ承壓水δD(a)與δ18O空間分布(b)
3.4.3 地下水的演化
江蘇洋口港地區(qū)天然地下水流動(dòng)是自晚更新世末期以來，伴隨著冰退、海平面上升及海洋性氣候的變化而調(diào)整到目前的模式.距今約25000~15000 a的末次盛冰期，海平面低于現(xiàn)今約120m(張人權(quán)等，2013).由于不存在全新世地層，現(xiàn)在的深層含水層埋藏較淺，與大氣降水和地表水聯(lián)系密切，水循環(huán)交替積極，地下水流動(dòng)速度快，補(bǔ)給條件好，降水及地表水的補(bǔ)給有足夠長(zhǎng)的時(shí)間交替出含水層中老的地下水，使含水層中保存有距今30000~10000 a的地下水.在形成及徑流過程中，遭遇較強(qiáng)烈的蒸發(fā)作用和陽(yáng)離子交換作用，具有氫氧穩(wěn)定同位素值低，氚含量低的同位素特點(diǎn).隨著全新世暖期的到來，海平面升高，排泄基準(zhǔn)面抬升，伴隨著地層的沉積形成，地下水流動(dòng)減緩，冰期的地下水未受海水入侵影響滯留于含水層之中而得以保存.
自 20 世紀(jì) 70 年代開始，研究區(qū)大量深層地下水開發(fā)利用，深層承壓地下水水位持續(xù)下降，水位降落漏斗產(chǎn)生并不斷擴(kuò)展，地下水流動(dòng)模式發(fā)生變化，淺層與深層地下水之間的水力聯(lián)系加強(qiáng)，部分地區(qū)含水層貫通引起越流補(bǔ)給而出現(xiàn)含氚水.根據(jù)馬塘鎮(zhèn)水位降落漏斗區(qū)的δ18O和δD同位素?cái)?shù)據(jù)，應(yīng)用二端元法計(jì)算漏斗區(qū)補(bǔ)給第Ⅲ承壓水的淺層地下水比例約為20%.具體參見 污水處理技術(shù)資料或污水技術(shù)資料更多相關(guān)技術(shù)文檔。
4 結(jié)論
1)根據(jù)全球大氣降水監(jiān)測(cè)數(shù)據(jù)庫(kù)建立研究區(qū)降水線方程為 δD=8.49δ18O+17.71(N=58，R2= 0.967，p<0.01).在δD-δ18O關(guān)系圖中，河水沿著徑流方向受到海水與蒸發(fā)作用的影響，表現(xiàn)出氘氧同位素的沿程富集.
2)研究區(qū)淺層地下水為現(xiàn)代氣候條件下的大氣降水補(bǔ)給.深部承壓含水層地下水的放射性14C年齡主要為15000~26000 a，其穩(wěn)定同位素δ18O和δD值比現(xiàn)代地下水低，不含氚，是晚更新世末次冰期(大理冰期)的盛冰期降水入滲補(bǔ)給.
3)研究區(qū)地下水具有明顯的呈層性，層間聯(lián)系較差，不同含水層地下水中δD、δ18O空間差異明顯.淺層地下水徑流速度快，可更新能力較強(qiáng).深層地下水開采的主要為儲(chǔ)存資源，自然狀態(tài)下，深層地下水不受現(xiàn)代降水影響，徑流緩慢，總體上處于封閉-半封閉狀態(tài)，可更新能力弱.
4)研究區(qū)天然地下水流動(dòng)是自晚更新世末期以來，伴隨著冰退、海平面上升而調(diào)整到目前的模式.末次盛冰期以來的自然地理及地質(zhì)發(fā)展史，決定著江蘇沿海平原第四系地下水流系統(tǒng)演變格局，現(xiàn)代人類活動(dòng)加強(qiáng)了淺層與深層地下水之間的水力聯(lián)系，部分地區(qū)水位下降明顯，含水層貫通引起越流補(bǔ)給而出現(xiàn)含氚水.
