邱向榮,鄧清海
(1.廣東省地質(zhì)技術(shù)工程咨詢公司,廣東 廣州 510080;2.山東科技大學(xué)地球科學(xué)與工程學(xué)院,山東 青島 266590)
水文地球化學(xué)研究,是水文地質(zhì)學(xué)中的基本研究方法之一,用于分析地下水的成分組成、來源、遷移及其形成過程等地下水的基本問題,在飲用水、礦泉水、地?zé)崴暗叵滤廴镜确矫鎽?yīng)用廣泛[1-9]。
地處粵港澳大灣區(qū)中心的廣東省中山市,瀕臨南海,緊鄰廣、深、珠和港澳,地理位置十分優(yōu)越,經(jīng)濟發(fā)達,礦泉水消費市場廣闊。而中山市轄區(qū)內(nèi)礦泉水資源豐富,且集中賦存于結(jié)晶巖分布區(qū)內(nèi),在由結(jié)晶巖構(gòu)成的丘陵山區(qū)內(nèi),打井只要有水基本都是礦泉水。中山市產(chǎn)出的礦泉水屬于低鈉低礦化度偏硅酸型礦泉水,由于口感甘甜清冽,在南方廣受歡迎。因此,研究中山市礦泉水形成的地球化學(xué)過程,對研究南方結(jié)晶巖地區(qū)微礦化水的形成有重要借鑒意義,也可以指導(dǎo)類似地區(qū)礦泉水資源的勘查與評價。
研究區(qū)地層由老至新有薊縣-青白口系、寒武系、泥盆系、侏羅系、白堊系和第四系,其中前第四紀地層零星出露,包括元古代變質(zhì)巖和寒武紀淺海類復(fù)理式沉積碎屑巖,河流相—濱淺海相沉積的泥盆紀地層,以及侏羅紀火山巖和白堊紀陸相紅色碎屑沉積巖(圖1)[10]。
第四紀松散沉積物廣泛分布于市內(nèi),成因類型復(fù)雜,主要為河口三角洲沉積,局部為陸相沉積,巖相巖性及沉積物厚度多變,厚度自10~70m不等。
研究區(qū)巖漿巖廣泛分布。侵入巖形成于奧陶紀、侏羅紀及白堊紀,其中以白堊紀侵入巖最為發(fā)育,巖性以二長花崗巖、花崗斑巖為主(表1),因此在研究區(qū)形成了以花崗巖類為主的廣泛結(jié)晶巖分布區(qū)(圖1)。
區(qū)內(nèi)主要發(fā)育斷裂構(gòu)造,不發(fā)育褶皺構(gòu)造。NE向、NW向斷裂帶及EW向構(gòu)造帶構(gòu)成了中山地區(qū)主體構(gòu)造格局(圖1)。
研究區(qū)地下水類型主要分為松散巖類孔隙含水巖組和基巖裂隙含水巖組,其分布及富水性特征如圖2所示。
(1)松散巖類孔隙含水巖組。松散巖類孔隙含水巖組是中山市地下水的主要類型,主要分布于中山市北部三角洲平原。第四系松散巖類含水巖組從平面可分為內(nèi)陸河谷型和三角洲沉積型,內(nèi)陸河谷型揭露富水性較好的含水巖組有全新世洪積層。三角洲沉積型在垂向揭露富水性較好的含水巖組有桂州組杏壇段和更新世禮樂組石排段。
(2)基巖裂隙含水巖組。巖漿巖類裂隙含水巖組:巖漿巖類含水巖組在區(qū)內(nèi)分布較廣,主要分布在長江水庫周邊地區(qū)、五桂山山脈一帶其他丘陵地區(qū)及北部殘丘。水量極貧乏—中等,水質(zhì)較好,有供水價值。枯季地下徑流模數(shù)為0.009~14.308L/s·km2,泉流量0.01~1.70L/s。
碎屑巖類裂隙含水巖組:中山市碎屑巖類裂隙含水巖組主要分布于沙崗—平嵐和布洲村一帶的低山丘陵區(qū),水量貧乏,分布面積也相對較少??菁镜叵聫搅髂?shù)比較小,泉流量0.035~0.450L/s。
變質(zhì)巖類裂隙含水巖組:變質(zhì)巖類裂隙含水巖組在研究區(qū)內(nèi)水量均為極貧乏(圖2),該區(qū)內(nèi)未見溪溝泉水出露。
中山市周邊被地下咸水包圍,因此在中山市結(jié)晶巖分布區(qū)賦存的地下淡水主要來源于地表水體和天然降雨補給,其中較大的長江水庫在水源補給方面起到了較大的作用。
(1)松散巖類孔隙水補徑排條件。松散巖類孔隙水主要接受大氣降水垂直補給,徑流途徑短,大部分地區(qū)為就地補給就地排泄,無明顯的補、徑、排區(qū)。
(2)基巖裂隙水補徑排條件?;鶐r裂隙水主要靠大氣降水直接滲入補給,其次為地表水和上伏殘坡積層孔隙水的間接補給?;鶐r裂隙水分布區(qū)地形陡峻,切割強度大,水系發(fā)育,地下水流向受最低級分水嶺控制,一般徑流途徑短,經(jīng)過短距離運移即排泄于溝谷中。
考慮到同一地區(qū)不同水體化學(xué)成分的形成往往具有某些內(nèi)在的聯(lián)系,為了利于對比分析,故樣品的選擇涵蓋了中山市6類主要水體。
(1)飲用天然礦泉水,每個主要的水源地選取一組代表性水樣。
(2)泉水,區(qū)內(nèi)泉點不多,代表性選擇。
(3)地表水,包括河溪水和水庫水,區(qū)內(nèi)河道較多,不同程度受到污染,重點選取受污染較輕微的河水,而水庫水則選取區(qū)內(nèi)庫容最大、也是飲用水源的長江水庫。
(4)溫泉水,主要熱田區(qū)選擇1組。
(5)海水,在河流入海口附近選擇混合水樣。
(6)民井,數(shù)量較多,雖然多數(shù)并非飲用水井,但井位分布于不同巖性、不同地貌單元,能獲得淺部地下水的基本特征,取樣數(shù)量較多。最終共取水樣44組,其中礦泉水7組、泉水3組、地表水4組、溫泉水3組、海水1組、民井水26組。
樣品的理化指標分析分為2種,一種是簡分析,主要進行宏量組分分析,考慮到區(qū)內(nèi)礦泉水產(chǎn)于結(jié)晶鹽,特征組分為偏硅酸,且有含氟背景,故簡分析均加做氟和偏硅酸分析,簡分析樣品共30組。另一種是全分析,按照現(xiàn)行天然礦泉水規(guī)范[11]的要求進行理化指標全分析(未進行微生物分析),全分析樣品共14組。有關(guān)樣品的概況和基本化學(xué)特征見表2。
表2 樣品基本信息及化學(xué)組分
4.1.1 水化學(xué)特征的Piper三線圖分析
根據(jù)所有水樣的Piper三線圖(圖3),中山市主要水體的化學(xué)特征總體上可以劃分為3組。
1—礦泉水;2—泉水;3—地表水;4—民井水;5—海水;6—溫泉水圖3 中山市主要水體水化學(xué)Piper三線圖
第一組位于三線圖菱形的右側(cè),包括了3組溫泉水、海水和1組民井水,這組水的主要特征是陰離子中氯離子明顯高于其他離子,陽離子則以鈉或鈉鈣為主,化學(xué)類型主要為Cl-Na,Cl-Na·Ca型。由于海水樣品取自于海岸帶,是河水和海水的混合結(jié)果,故該組水質(zhì)代表大氣水明顯接受海水的混合,由此可以判斷溫泉水和個別淺井也存在明顯海水混合作用。
第三組位于三線圖菱形區(qū)左側(cè),均為礦泉水樣。這組水的陰離子很單一,均為重碳酸根離子,陽離子則包括鈉和鈣。從水樣投影位置上看,靠近Ca+Mg軸的樣品多于靠近K+Na軸,即鈣占主導(dǎo)位置。
從3個分組的關(guān)系可以看出,礦泉水所在的第三組靠近第二組而遠離第一組,說明礦泉水與河湖水、泉水以及民井水之間的關(guān)聯(lián)性比較大,而與深循環(huán)的溫泉水關(guān)聯(lián)性不大,這意味著本區(qū)礦泉水屬于淺循環(huán)地下水,其水化學(xué)的形成過程與深循環(huán)溫泉水之間有較大差異。
4.1.2 水化學(xué)特征的Schoeller圖分析
根據(jù)水質(zhì)主要組的Schoeller圖可見(圖4),曲線圖出現(xiàn)了2個分布區(qū)域,其中位于曲線圖上方區(qū)域的水樣包括溫泉水和海水,位于曲線圖中部的區(qū)域包括了所有的礦泉水、地表水、泉水和民井水。如上所述,海水樣是河、?;旌纤蔛choeller圖同樣表明溫泉水存在明顯海水混合作用。曲線圖中部區(qū)域各種水樣交錯混雜,沒有很好的分帶性,意味著礦泉水與除溫泉水和海水以外的水體關(guān)聯(lián)度較大,控制水化學(xué)特征的要素比較接近。相對而言,礦泉水比其他水體有更低的鎂、硫酸根和氯離子濃度,因此礦泉水樣品的投影曲線總體位于靠下部的區(qū)域。
1—礦泉水;2—泉水;3—地表水;4—民井水;5—海水;6—溫泉水圖4 中山市主要水體水化學(xué)Schoeller圖
中山市的礦泉水屬于低鈉低礦化度偏硅酸型礦泉水,已開發(fā)利用的水源地均位于結(jié)晶巖區(qū),基巖主要是奧陶紀到白堊紀時期的侵入花崗巖,礦物成分主要為斜長石、鉀長石、石英和少量暗色礦物等。因形成年代漫長,故疊加多期次的區(qū)域構(gòu)造形跡,原巖礦物成分出現(xiàn)不同程度的交代蝕變。華南沿海風(fēng)化作用強烈,巖體表層普遍有20~30m厚的風(fēng)化殼,礦物成分比較復(fù)雜。
根據(jù)礦泉水分布區(qū)的礦物學(xué)特征,結(jié)晶巖的造巖礦物顯然是水中化學(xué)組分的重要來源,另外由于風(fēng)化或蝕變礦物比較容易被地下水溶濾,因此巖土層中廣泛存在的巖鹽、碳酸鹽和硫酸鹽礦物往往也是水質(zhì)組分的主要來源。下面根據(jù)礦泉水主要化學(xué)組分的相關(guān)性來分析水質(zhì)組分的來源。
圖5 礦泉水主要組分關(guān)系圖
為了便于對比,測年時采取了礦泉水和溫泉水兩類水樣。礦泉水測年樣品選擇了唯一一處礦泉水上升泉,以獲得天然滲流場狀態(tài)下礦泉水的年齡,避免人工抽水因加劇水交替而引起水年齡的改變。用于對比的溫泉水選取2處,一處是位于海灘中的溫泉自流井(高潮位時原泉口低于海平面),抽水狀態(tài)下其水溫達到100℃[13];另一處是選擇因明顯超采并引起地面沉降的地?zé)峋?GWB7號樣),該地?zé)崽镌谒疁丶s90℃(1)廣東省地質(zhì)局水文工程地質(zhì)二大隊,廣東省中山縣雍陌地?zé)釁^(qū)詳查報告,1982年。,現(xiàn)地?zé)峋疁y溫只有56.0℃,說明淺部冷水的混合作用加劇。
有關(guān)年齡同位素樣品的相關(guān)數(shù)據(jù)見表3。測試結(jié)果表明,礦泉水的年齡遠遠小于溫泉水的年齡。另外,溫泉水年齡差異也不小,GWB7號樣品的年齡明顯小于GWB5號樣品,這客觀上雖然與后者的水溫更高、循環(huán)深度更大有關(guān),但前者因超采加大冷水的混合應(yīng)該是測試年齡明顯偏小的主要原因。從礦泉水和溫泉水的測年結(jié)果可以看出,深循環(huán)的溫泉水具有比較大的測年值,而礦泉水比較小的年齡值代表水循環(huán)深度比較淺,應(yīng)屬于淺循環(huán)基巖裂隙水。
表3 年齡同位素樣品數(shù)據(jù)匯總表
進行反向模擬的目的是試圖分析礦泉水從入滲區(qū)到排泄區(qū)過程中參與水—巖作用的反應(yīng)、礦物成分以及物質(zhì)的轉(zhuǎn)移量等,以解釋礦泉水形成的水文地球演化過程[14-17]。反向模擬需要提供上下游水樣,這里選取長江水庫的樣品為入滲區(qū)水樣(表2中的39號樣),排泄區(qū)水樣選取上升泉(表2中的19號樣),有關(guān)水樣的主要組分見表4。
表4 反向模擬上下游水質(zhì)組分 單位:mg/L
選擇長江水庫作為入滲區(qū)水樣主要考慮代表性,因長江水庫處于中山市五桂山中北部,是中山市最大的水庫,而五桂山地區(qū)是中山市礦泉水分布最廣,開發(fā)利用最集中的地區(qū),水庫水是大氣降水和短暫循環(huán)地下水的混合,以此作為丘陵區(qū)礦泉水的補給來源具有普遍意義。選擇該上升泉的原因一方面是泉點與水庫所處的巖石礦物背景接近,該泉的允許開采量也適中(B+C級儲量為171m3/d),另外是考慮盡量選擇天然水循環(huán)條件下礦泉水的化學(xué)特征,畢竟開采井抽水會加速水代謝,一定程度上改變礦泉水的天然化學(xué)特征。
反向模擬的難點和重點是確定反應(yīng)相組分,通常既需要掌握含水系統(tǒng)的礦物學(xué)特征,也需根據(jù)水質(zhì)特征對其化學(xué)來源有基本的判斷,另外可結(jié)合水中礦物的平衡分析或者水和礦物的同位素組分特征等因素綜合確定。
根據(jù)前文分析,中山市礦泉水分布區(qū)巖性主要為花崗巖,造巖礦物主要為斜長石、鉀長石、石英和少量暗色礦物。水質(zhì)組分相關(guān)性分析提示了礦泉水演化過程中溶濾了碳酸鹽、硫酸鹽和巖鹽等礦物。水中主要礦物的平衡狀態(tài)對于合理選擇反應(yīng)相有幫助,也利用確定反應(yīng)過程中不同相的沉淀—溶解狀態(tài),為此對入滲水樣(39)中礦物的飽和指數(shù)進行分析(表5)。根據(jù)以上相關(guān)信息,通過比選最終選擇反應(yīng)相為CO2(g)、斜長石、二氧化硅、白云石、硬石膏、螢石、天青石、巖鹽和水鋁礦共9種。
表5 上游水樣(39)中主要礦物飽和指數(shù)
模擬計算通過PHREEQC軟件來完成,計算過程中2種溶液及相關(guān)元素的不確定度均取0.025,并設(shè)置斜長石為溶解、二氧化硅和水鋁礦為沉淀,其他反應(yīng)相不受限制。通過求解獲得唯一符合的反應(yīng)模型,有關(guān)反應(yīng)相及其摩爾轉(zhuǎn)移見表6。
表6 反向模擬計算結(jié)果
根據(jù)模擬結(jié)果可見,礦泉水從入滲到排泄的循環(huán)過程中,斜長石、白云石、硬石膏、螢石、天青石等礦物出現(xiàn)溶解,并且需比較多的二氧化碳氣體參與,水—巖反應(yīng)的結(jié)果僅水鋁礦沉淀,二氧化硅基本沒有參與反應(yīng)。由于反應(yīng)方式為溶解的礦物均離不開二氧化碳,可見南方地區(qū)廣泛分布于結(jié)晶巖地區(qū)的偏硅酸型微礦化水的形成,碳酸平衡體系起著重要控制作用,如果二氧化碳分壓太小,對這類礦泉水的形成是不利的,因此實踐中常常可見這類低礦化度偏硅酸型礦泉水常呈弱酸性反應(yīng)。另外,模擬結(jié)果也表明,礦泉水中的達標組分偏硅酸主要來自長石類礦物的水解,與二氧化硅礦物的關(guān)系并不密切。
中山市的飲用天然礦泉水均分布于結(jié)晶巖地區(qū),巖性主要為奧陶紀到白堊紀侵入花崗巖,礦泉水的特征組分為偏硅酸,屬低鈉低礦化度偏硅酸型礦泉水,水化學(xué)類型主要為HCO3-Ca,HCO3-Na,CO3-Na·Ca型,是南方地區(qū)廣泛分布的微礦化水。
不同水體的水質(zhì)組分對比顯示:中山市的飲用天然礦泉水與民井水、泉水和地表水體(河水和水庫水)比較接近,水化學(xué)形成的控制因素類似,但礦泉水的硫酸根、氯離子和鎂離子濃度偏低。而礦泉水與溫泉水的差異比較大,一方面礦泉水的年齡遠小于溫泉水,另一方面溫泉水普遍具有海水混合特征,礦泉水則表現(xiàn)為溶濾水。
礦泉水主要組分的相關(guān)性分析表明:礦泉水化學(xué)成分的來源除了原巖礦物外,還與風(fēng)化、蝕變礦物密切相關(guān),部分組分來自于碳酸鹽、硫酸鹽和巖鹽等礦物的溶解。根據(jù)反向模擬結(jié)果,在礦泉水形成演化過程中,攜帶二氧化碳氣體的入滲水通過溶濾斜長石、白云石、硬石膏、天青石、螢石和巖鹽等礦物,同時出現(xiàn)了水鋁礦沉淀;特征組分偏硅酸主要來自長石類礦物的的水解,與二氧化硅礦物關(guān)系不密切;礦泉水化學(xué)演化過程中碳酸平衡體系起著重要控制作用,充足的二氧化碳是這類礦泉水形成的重要因素,這應(yīng)該是華南地區(qū)低礦化度偏硅酸型礦泉水普遍呈弱酸性反應(yīng)的主要原因。
本研究結(jié)果查明了中山市礦泉水的水文地球化學(xué)過程,為下一步深入研究南方結(jié)晶巖地區(qū)微礦化水的形成機制提供了參考依據(jù),也可以對類似礦泉水資源的勘查與評價提供方法參考。