摘 要:為探究商丹盆地地下水水化學(xué)特征及其形成作用,支撐商洛市商丹循環(huán)工業(yè)經(jīng)濟(jì)園區(qū)建設(shè),選擇2023 年商丹盆地地下水資源調(diào)查獲得的24 組地下水水化學(xué)數(shù)據(jù),綜合運(yùn)用數(shù)理統(tǒng)計(jì)、Piper 三線圖、Gibbs 圖、離子比分析和主成分分析等方法進(jìn)行水文地球化學(xué)分析,結(jié)果表明:研究區(qū)地下水為淡水,整體呈中性,總硬度為200~698 mg/ L,平均值為427.79 mg/ L;溶解性總固體含量為243~960 mg/ L,平均值為545.29 mg/ L;地下水水化學(xué)類(lèi)型以HCO3 -Ca 型為主,水化學(xué)成分受碳酸鹽巖和硅酸鹽巖風(fēng)化溶解作用共同控制;地下水中Na+和K+主要受鹽巖溶解作用控制,Ca2+和Mg2+主要受碳酸鹽巖和硫酸鹽巖風(fēng)化溶解作用控制;地下水水化學(xué)特征受水巖相互作用、陽(yáng)離子交換作用和人類(lèi)活動(dòng)的共同影響,地下水在運(yùn)移過(guò)程中以正向陽(yáng)離子交換為主,礦業(yè)開(kāi)采、農(nóng)業(yè)生產(chǎn)和生活污水排放等人類(lèi)活動(dòng)對(duì)地下水化學(xué)成分影響較大。
關(guān)鍵詞:地下水;水化學(xué)特征;水巖相互作用;主成分分析;商丹盆地
中圖分類(lèi)號(hào):P641 文獻(xiàn)標(biāo)志碼:A doi:10.3969/ j.issn.1000-1379.2024.12.016
引用格式:趙浩,李曉明,趙立磊,等.商丹盆地地下水水化學(xué)特征及成因機(jī)制[J].人民黃河,2024,46(12):98-104.
地下水是區(qū)域水資源的重要組成部分,作為人類(lèi)生產(chǎn)和生活的重要水源,地下水在維系水生態(tài)系統(tǒng)平衡、工農(nóng)業(yè)生產(chǎn)和人類(lèi)生存發(fā)展等方面發(fā)揮著重要作用[1-3] 。地下水質(zhì)量事關(guān)人民健康和經(jīng)濟(jì)社會(huì)穩(wěn)定發(fā)展,地下水水化學(xué)成分與賦存環(huán)境關(guān)系密切[4] 。地下水化學(xué)組分監(jiān)測(cè)有利于分析地下水水化學(xué)演化過(guò)程和潛在影響因素[5] 。有學(xué)者從水巖相互作用[6] 、陽(yáng)離子交換[7-8] 、人類(lèi)活動(dòng)影響[9] 等方面對(duì)地下水水化學(xué)特征的形成和演化進(jìn)行了大量研究。
商丹盆地位于商洛市,是南水北調(diào)中線工程水源區(qū)丹江流域的重要城鎮(zhèn)聚集地,水資源保障是支撐商丹盆地經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展的關(guān)鍵。隨著區(qū)域內(nèi)城鎮(zhèn)化進(jìn)程的加快和工業(yè)園區(qū)的建設(shè),生產(chǎn)生活用水增加,農(nóng)業(yè)面源污染和工業(yè)廢水排放等對(duì)地下水水質(zhì)造成影響,查明區(qū)域地下水水化學(xué)特征并分析其成因機(jī)制具有重要意義。鑒于此,以2023 年商丹盆地采集的24組地下水水質(zhì)數(shù)據(jù)為基礎(chǔ),運(yùn)用數(shù)理統(tǒng)計(jì)、Piper 三線圖、Gibbs 模型、離子比值和主成分分析等方法進(jìn)行水文地球化學(xué)綜合研究,分析研究區(qū)地下水水化學(xué)特征及成因機(jī)制,以期為商丹盆地地下水資源管理提供依據(jù),助力商丹循環(huán)工業(yè)經(jīng)濟(jì)園區(qū)建設(shè)。
1 研究區(qū)概況
商丹盆地地處秦嶺南麓、丹江上游,大地構(gòu)造受商丹大斷裂控制,屬山間斷陷型構(gòu)造盆地。商丹盆地東西長(zhǎng)約60 km,南北寬約5 km,丹江沿北西—南東方向流經(jīng)盆地,河寬150~250 m,谷寬1 000~3 000 m,平均比降3‰。丹江河道迂回曲折,受河流長(zhǎng)期下切作用影響,丹江兩岸斷續(xù)分布1~4 級(jí)階地,構(gòu)成了現(xiàn)今比較開(kāi)闊的彎道盆地[10-11] 。商丹盆地屬北亞熱帶和暖溫帶過(guò)渡性季風(fēng)氣候區(qū),溫暖濕潤(rùn)、四季分明,河流水系發(fā)育,年均氣溫7.8~13.9 ℃,多年平均(1958—2019年)降水量700.86 mm,降水季節(jié)分布不均,主要集中在6—9 月。
研究區(qū)地下水類(lèi)型主要為松散巖類(lèi)孔隙水和碎屑巖類(lèi)孔隙裂隙水。松散巖類(lèi)孔隙水主要分布于丹江兩岸河漫灘及一級(jí)階地,含水層以第四系松散堆積的砂、卵礫石層為主,水位埋深1.2~19.2 m,水量豐富;碎屑巖類(lèi)孔隙裂隙水主要分布于丹江河谷兩側(cè)山前地帶,含水巖組以第三系砂礫巖為主。地下水主要接受兩側(cè)基巖山區(qū)的側(cè)向補(bǔ)給、大氣降水入滲補(bǔ)給、河水入滲補(bǔ)給和農(nóng)田灌溉入滲補(bǔ)給,排泄方式以地下水徑流排泄和人工開(kāi)采為主,徑流方向?yàn)檠嘏璧剡吘壪蛑胁苛鲃?dòng)。根據(jù)《商洛市水資源公報(bào)》(2020 年),盆地地下水供水量(商州區(qū)、丹鳳縣)約占總供水量的30.69%。地下水主要用于生產(chǎn)生活和農(nóng)田灌溉,其中商州區(qū)以居民生活用水(1 026 萬(wàn)m3)和工業(yè)用水(681 萬(wàn)m3 )為主,丹鳳縣以農(nóng)田灌溉用水(480 萬(wàn)m3)和居民生活用水(244 萬(wàn)m3)為主。2022 年,研究區(qū)常住人口約71.3萬(wàn)人,國(guó)民生產(chǎn)總值約268 億元。
2 樣品采集與測(cè)試方法
2023 年6 月,在水文地質(zhì)調(diào)查的基礎(chǔ)上共采集盆地區(qū)地下水樣品24 組,樣品主要來(lái)源于飲用水井和灌溉井,采樣點(diǎn)分布見(jiàn)圖1。嚴(yán)格按照《地下水環(huán)境監(jiān)測(cè)技術(shù)規(guī)范》(HJ 164—2020)進(jìn)行樣品采集、保存和運(yùn)輸。取樣器材為5 L 聚乙烯塑料瓶,采樣前先用水樣將取樣瓶和瓶蓋潤(rùn)洗3~4 次,不添加任何保護(hù)劑,采樣完成后擰緊瓶蓋,粘貼標(biāo)簽,標(biāo)明采樣時(shí)間、樣品編號(hào)等。
本次研究主要測(cè)試指標(biāo)包括pH 值、總硬度(TH)及溶解性總固體( TDS)、K+、Na+、Ca2+、Mg2+、Cl-、SO2-4 、HCO-3 、NO-3 含量等,其中pH 值現(xiàn)場(chǎng)測(cè)定,其余指標(biāo)測(cè)試由陜西省地質(zhì)礦產(chǎn)實(shí)驗(yàn)研究所有限公司完成。TH 和TDS 根據(jù)《生活飲用水標(biāo)準(zhǔn)檢驗(yàn)方法感官性狀和物理指標(biāo)》(GB/ T 5750.4—2006),采用電子天平和滴定管測(cè)定,檢出限分別為1.0、4.0 mg/ L;K+、Na+、Ca2+、Mg2+含量根據(jù)《水質(zhì)32 種元素的測(cè)定電感耦合等離子體發(fā)射光譜法》(HJ 776—2015)測(cè)定,測(cè)試儀器為Agilent 5110 電感耦合等離子體光譜儀,檢出限分別為0.05、0.12、0.02、0.003 mg/ L;Cl-、SO2-4 、NO-3 含量根據(jù)《水質(zhì)無(wú)機(jī)陰離子的測(cè)定離子色譜法》(HJ84—2016)測(cè)定,測(cè)試儀器為ICS1100 離子色譜儀,檢出限分別為1.0、1.0、0.45 mg/ L;HCO-3 含量根據(jù)《食品安全國(guó)家標(biāo)準(zhǔn)飲用天然礦泉水檢驗(yàn)方法》(GB 8538—2016),使用滴定管測(cè)定,檢出限為3.0 mg/ L。各指標(biāo)檢出限及測(cè)試結(jié)果均滿(mǎn)足要求。
3 結(jié)果與討論
3.1 地下水水化學(xué)特征
1)水化學(xué)成分。商丹盆地地下水水化學(xué)參數(shù)特征值統(tǒng)計(jì)結(jié)果見(jiàn)表1。地下水pH 值為6.77~7.94,均值為7.21,呈中性,空間變異?。郏保?;總硬度為200~698mg/ L,平均值為427.79 mg/ L,微硬水、硬水、極硬水占比分別為16.66%、41.67%、41.67%;TDS 為243 ~ 960mg/ L,平均值為545.29 mg/ L,屬于淡水。地下水中陰離子含量表現(xiàn)為HCO-3 >SO2-4 >NO-3 >Cl-,陽(yáng)離子含量表現(xiàn)為Ca+ >Na+ >Mg+ >K+,陰陽(yáng)離子中占主要優(yōu)勢(shì)的分別是HCO-3 和Ca+。TH、TDS 和其他各離子的變異系數(shù)為26.95%~88.76%,表明各離子變異性較強(qiáng)且空間差異性較大。其中,NO-3 的變異系數(shù)最大,超過(guò)《地下水質(zhì)量標(biāo)準(zhǔn)》(GB/ T14848—2017) Ⅲ類(lèi)水標(biāo)準(zhǔn)88.58mg/ L(以NO-3 計(jì))和《生活飲用水衛(wèi)生標(biāo)準(zhǔn)(GB5749—2022)》44.29 mg/ L(以NO-3 計(jì)) 的樣品占比分別為25.00%、54.17%,說(shuō)明研究區(qū)地下水水質(zhì)可能受到人類(lèi)活動(dòng)的影響。
2)水化學(xué)類(lèi)型。運(yùn)用origin 軟件繪制研究區(qū)地下水水化學(xué)Piper 三線圖(見(jiàn)圖2),分析地下水水化學(xué)類(lèi)型及特征[1,12] 。主菱形圖中樣點(diǎn)主要落在左上方,表明地下水中堿土金屬離子毫克當(dāng)量百分比超過(guò)堿金屬離子,弱酸根離子毫克當(dāng)量百分比大部分超過(guò)強(qiáng)酸根離子,碳酸鹽硬度>50%;在陽(yáng)離子圖中,樣點(diǎn)集中分布于Ca2+軸和Mg2+軸相交端,表明地下水中陽(yáng)離子以Ca2+為主;在陰離子圖中,樣點(diǎn)主要分布于HCO-3 +CO2-3軸和Cl- 軸相交端且靠近HCO-3 +CO2-3 軸,表明研究區(qū)大部分地下水HCO-3 含量遠(yuǎn)大于SO2-4 和Cl-,可初步反圖2 地下水水化學(xué)Piper 三線圖映研究區(qū)地下水化學(xué)組分主要來(lái)源于碳酸鹽巖的溶解。采用舒卡列夫分類(lèi)法[1] 進(jìn)行地下水水化學(xué)類(lèi)型劃分,結(jié)果表明研究區(qū)地下水水化學(xué)類(lèi)型主要為HCO3 -Ca 型,占比79.16%(19 組);其次為HCO3 ·SO4 -Ca 型,占比為12.50%(3 組);還有HCO3 ·SO4 -Ca·Mg 型和HCO3 ·Cl-Ca 型水各一組,占比都為4.17%。盆地部分區(qū)域出現(xiàn)以Cl- 和SO2-4 為主導(dǎo)離子的地下水,表明地下水受多種水化學(xué)作用和人類(lèi)活動(dòng)的共同影響[13] 。
3)水化學(xué)指標(biāo)相關(guān)性分析。采用相關(guān)性分析方法研究地下水水化學(xué)組分之間的相互關(guān)系和影響程度,識(shí)別不同離子組分是否為同一來(lái)源或具有相近的遷移轉(zhuǎn)化途徑[14-15] 。研究區(qū)地下水水化學(xué)組分相關(guān)系數(shù)統(tǒng)計(jì)(見(jiàn)表2)可知,TDS 與Na+、Ca2+、Mg2+、Cl-、SO2-4 、NO-3 顯著正相關(guān),相關(guān)系數(shù)依次為0.801、0.976、0.889、0.824、0.763、0.824,表明這些離子對(duì)TDS 貢獻(xiàn)較大,其中Ca2+ 和Mg2+ 貢獻(xiàn)尤為突出;HCO-3 與Na+、Ca2+、Mg2+顯著正相關(guān),相關(guān)系數(shù)分別為0.651、0.667、0.526,表明四者具有相對(duì)一致的來(lái)源;Ca2+ 與Mg2+ 相關(guān)系數(shù)為0.831,兩者極顯著正相關(guān),表明兩種離子具有相同的來(lái)源;SO2-4 與Ca2+、Mg2+ 顯著正相關(guān),相關(guān)系數(shù)分別為0.661、0.839,表明部分Ca2+和Mg2+可能來(lái)自于硫酸鹽礦物的溶解;SO2-4 與Na+、Cl- 也顯著正相關(guān),相關(guān)系數(shù)分別為0.560 和0.579,說(shuō)明蒸發(fā)巖對(duì)地下水離子具有一定貢獻(xiàn);Cl- 與NO-3 相關(guān)系數(shù)為0.790,兩者極顯著正相關(guān),表明兩者具有相似的物質(zhì)來(lái)源。
3.2 水化學(xué)特征成因分析
1)水巖模型分析。常用Gibbs 圖研究水巖相互作用對(duì)地下水水化學(xué)組分的影響機(jī)制,利用TDS 含量與ρCl- / ρ(Cl- +HCO-3 )或ρ Na+ / ρ (Na+ +Ca2+ )(ρ 表示離子質(zhì)量濃度)的關(guān)系將天然水化學(xué)組分的演化分為巖石風(fēng)化控制型、大氣降水控制型和蒸發(fā)濃縮控制型[5,16] 。研究區(qū)地下水TDS 含量為243~960 mg/ L,陽(yáng)離子ρ Na+ / ρ ( Na+ +Ca2+ ) 和陰離子ρ Cl- / ρ ( Cl- +HCO-3 )分別為0.12~0.25 和0.03~0.23,見(jiàn)圖3。陰陽(yáng)離子質(zhì)量濃度比值均小于0.50,水樣點(diǎn)全部落在巖石風(fēng)化控制區(qū),且遠(yuǎn)離蒸發(fā)濃縮控制區(qū)和大氣降水控制區(qū),說(shuō)明巖石風(fēng)化作用是研究區(qū)地下水水化學(xué)組成的主要控制因素,大氣降水和蒸發(fā)濃縮作用的影響很小。
碳酸鹽巖、硅酸鹽巖和蒸發(fā)鹽巖礦物風(fēng)化溶解產(chǎn)物ρ Ca2+ / ρNa+、ρ Mg2+ / ρNa+、ρ HCO-3 / ρNa+ 不受流速、稀釋和蒸發(fā)作用的影響,通過(guò)地下水不同巖性端元間的對(duì)比可以判斷不同水巖相互作用關(guān)系,揭示水化學(xué)成因和主要離子來(lái)源[17] 。由圖4 可知,研究區(qū)樣品點(diǎn)均落在硅酸鹽巖和碳酸鹽巖端元之間且相對(duì)集中,說(shuō)明地下水中離子的來(lái)源受碳酸鹽巖和硅酸鹽巖風(fēng)化溶解作用共同控制,而受蒸發(fā)鹽巖的作用相對(duì)較小。
2)主要離子來(lái)源。根據(jù)地下水水巖相互作用可以分析水體溶濾作用和地下水中各化學(xué)成分的來(lái)源[5-7] 。地下水中Cl- 和SO2-4 主要來(lái)源于蒸發(fā)鹽巖的風(fēng)化溶解和含硫礦物的氧化,而HCO-3 則來(lái)源于碳酸鹽巖和硅酸鹽巖的溶解,可以利用γ ( SO2-4 + Cl- ) /γ HCO-3 (γ為離子毫克當(dāng)量)的關(guān)系分析水化學(xué)成分的主要來(lái)源[5,7] 。研究區(qū)γ(SO2-4 +Cl- ) / γHCO-3 為0.17~1.09,均值為0.44。由圖5(a)可知,地下水中91.67%的取樣點(diǎn)位于1 ∶ 1 直線下方,說(shuō)明地下水中化學(xué)成分主要來(lái)自碳酸鹽巖和硅酸鹽巖的溶解。
Na+的主要來(lái)源是鹽巖和硅酸鹽巖溶解以及大氣降水,可以通過(guò)γ Na+ / γCl- 的關(guān)系判斷Na+ 富集程度和來(lái)源。一般而言,大氣降水中γNa+ / γCl- 與海水接近,為0.86,當(dāng)γ Na+ / γCl- 大于1 時(shí)Na+ 來(lái)源以鹽巖溶解為主,小于1 時(shí)以硅酸鹽巖溶解為主[18] 。研究區(qū)采樣點(diǎn)γ Na+ / γ Cl- 為0.47 ~ 2.56,均值為1.66。由圖5(b)可知,研究區(qū)79.17%的采樣點(diǎn)γ Na+ / γCl- 大于1,說(shuō)明鹽巖溶解是Na+ 的主要來(lái)源,硅酸鹽巖溶解和大氣降水的貢獻(xiàn)較小。通過(guò)γ(Na+ +K+ ) / γCl- 的關(guān)系可以進(jìn)一步確定是否存在硅酸鹽礦物的溶解,由圖5(c)可知,大部分水樣落在了1 ∶ 1 直線上方,說(shuō)明區(qū)域水巖作用過(guò)程中發(fā)生鹽巖溶解的同時(shí),發(fā)生了其他含Na+、K+的硅酸鹽礦物溶解[5] 。
根據(jù)γ(Ca2+ +Mg2+ ) / γ(SO2-4 +HCO-3 )的關(guān)系可以判斷地下水中Ca2+ 和Mg2+ 的主要來(lái)源,當(dāng)γ (Ca2+ +Mg2+) / γ(SO2-4 +HCO-3 )在1 ∶ 1 直線附近時(shí),Ca2+、Mg2+主要來(lái)自方解石、白云石及石膏等碳酸鹽巖、硫酸鹽巖的溶解,比值小于1 時(shí)則主要來(lái)自硅酸鹽巖和蒸發(fā)鹽巖的溶解[7-8,19] 。由圖5(d)可知,地下水中有70.83%的采樣點(diǎn)γ(Ca2+ +Mg2+ ) / γ(SO2-4 +HCO-3 )接近于1 ∶ 1直線,且大于1 的采樣點(diǎn)也都靠近1 ∶ 1 直線,說(shuō)明研究區(qū)地下水中Ca2+ 和Mg2+ 主要來(lái)自碳酸鹽巖和硫酸鹽巖的溶解,但仍有少部分需要硅酸鹽巖和蒸發(fā)鹽巖的溶解來(lái)保持離子平衡[4] 。在碳酸鹽巖溶解區(qū),一般用γMg2+ / γCa2+與γHCO-3 的關(guān)系來(lái)判斷白云石和方解石的貢獻(xiàn)率[13] 。研究區(qū)采樣點(diǎn)γ Mg2+ / γ Ca2+ 為0.18 ~0.42,均小于0.50。由圖5(e)可知,采樣點(diǎn)全部落在方解石溶蝕區(qū),說(shuō)明Ca2+ 和Mg2+ 主要來(lái)源是方解石溶解。
3)陽(yáng)離子交換作用。陽(yáng)離子交換作用是指巖石和地下水相互作用過(guò)程中,在一定條件下顆粒吸附地下水中某些陽(yáng)離子,并將原來(lái)吸附的部分陽(yáng)離子轉(zhuǎn)化為地下水中的組分[5,20] 。Na+ +K+ -Cl- 與Ca2+ +Mg2+ -SO2-4 -HCO-3 的毫克當(dāng)量比值關(guān)系可以反映是否發(fā)生了陽(yáng)離子交換作用[7-8] 。由圖6(a)可知,采樣點(diǎn)分布在y =-0.371 x+0.693(R2 = 0.655)直線兩側(cè),且部分采樣點(diǎn)沿著y =-x直線分布,說(shuō)明地下水在運(yùn)移過(guò)程中發(fā)生了陽(yáng)離子交換作用。氯堿指數(shù)(CAI-Ⅰ和CAI-Ⅱ)可以進(jìn)一步探究陽(yáng)離子交換作用的方向和強(qiáng)度。通常情況下,當(dāng)CAI-Ⅰ和CAI-Ⅱ均為負(fù)值時(shí),說(shuō)明發(fā)生了正向的陽(yáng)離子交換作用,此時(shí)地下水中的Ca2+ 和Mg2+ 將含水層礦物中的Na+ 和K+ 交換出來(lái),地下水中的Na+ 和K+含量增加;當(dāng)CAI-Ⅰ和CAI-Ⅱ均為正值時(shí),則發(fā)生了反向的陽(yáng)離子交換作用,即地下水中的Na+ 和K+ 將含水層礦物中的Ca2+ 和Mg2+ 交換出來(lái),地下水中的Ca2+和Mg2+含量增加[21] 。根據(jù)圖6(b),有83.33%的采樣點(diǎn)CAI-Ⅰ和CAI-Ⅱ?yàn)樨?fù)值,表明研究區(qū)地下水主要發(fā)生了正向的陽(yáng)離子交換作用,導(dǎo)致地下水中Na+和K+含量增加,這與溶濾作用分析的結(jié)論(Na+ 的主要來(lái)源為鹽巖溶解)一致。
4)人類(lèi)活動(dòng)影響。除受自然因素作用外,人類(lèi)活動(dòng)一定程度上影響著地下水水化學(xué)組分的演化。通過(guò)水樣γ(Ca2+ +Mg2+ ) / γHCO-3 與γ(Na+ +K+ ) / γHCO-3 的關(guān)系可以反映碳酸鹽巖、硅酸鹽巖風(fēng)化以及人類(lèi)活動(dòng)對(duì)水化學(xué)組成的相對(duì)貢獻(xiàn),γ(Ca2+ +Mg2+ ) / γHCO-3 代表碳酸鹽巖和硅酸鹽巖風(fēng)化的Ca2+和Mg2+相對(duì)含量,其比值應(yīng)該小于1;γ(Na+ +K+) / γHCO-3 也代表碳酸鹽巖和硅酸鹽巖的風(fēng)化[22-23] 。
由圖7(a)可知,研究區(qū)大部分采樣點(diǎn)落在γ(Na+ +K+) / γHCO-3 = 0 和γ (Ca2+ +Mg2+ ) / γHCO-3 = 1 交叉點(diǎn)附近,且位于硅酸鹽巖風(fēng)化線右側(cè),說(shuō)明水化學(xué)組分以碳酸鹽巖風(fēng)化溶解為主,硅酸鹽巖風(fēng)化溶解貢獻(xiàn)較??;有62.50%的采樣點(diǎn)落在第一象限,說(shuō)明γ(Ca2+ +Mg2+ )和γ(Na+ +K+ )均大于γ HCO-3 ,超出的這部分陽(yáng)離子(Ca2+、Mg2+、Na+、K+)可能來(lái)源于人類(lèi)活動(dòng)的輸入,其陰離子可能以Cl- 或NO-3 的形式存在;另外有16.67%的采樣點(diǎn)落在第四象限,說(shuō)明存在大量Cl- 的輸入,這部分離子也可能主要來(lái)源于人類(lèi)活動(dòng)。
通過(guò)γSO2-4 / γCa2+ 和γNO-3 / γCa2+ 的關(guān)系可以進(jìn)一步分析地下水水化學(xué)組成受人類(lèi)活動(dòng)影響的程度[5] 。當(dāng)γSO2-4 / γCa2+大于γNO-3 / γCa2+時(shí),代表地下水受工礦生產(chǎn)活動(dòng)影響較大;反之地下水受農(nóng)業(yè)活動(dòng)和生活污水排放的影響較大[13] 。研究區(qū)采樣點(diǎn)γSO2-4 / γCa2+ 為0.12~0.57,均值為0.24;γ NO-3 / γCa2+ 為0.05~0.39,均值為0.15。由圖7(b)可知,地下水中66.67%的采樣點(diǎn)γSO2-4 / γCa2+大于γNO-3 / γCa2+,表明研究區(qū)地下水主要受到工礦生產(chǎn)活動(dòng)的影響,這與區(qū)域采礦活動(dòng)以及盆地內(nèi)金屬冶煉廠、化工廠等工業(yè)生產(chǎn)污水排放有關(guān),其余33.33% 的采樣點(diǎn)γ SO2-4 / γ Ca2+ 小于γNO-3 /γ Ca2+,總體靠近1 ∶ 1 直線,說(shuō)明部分區(qū)域地下水污染受農(nóng)業(yè)生產(chǎn)和生活污水排放的影響,但影響相對(duì)較小。由圖7(c) 可知,研究區(qū)γ NO-3 / γ HCO-3 與γ Cl- /γHCO-3 線性擬合較好(y =0.644x+0.043,R2 = 0.626),相關(guān)性較強(qiáng),表明地下水中Cl-和NO-3 可能具有相同的來(lái)源,即都受到農(nóng)業(yè)生產(chǎn)和生活污水排放的影響。
3.3 水化學(xué)控制因素分析
采用主成分分析法將水化學(xué)數(shù)據(jù)進(jìn)行降維處理,進(jìn)而提取主要變量來(lái)解釋水化學(xué)主要控制因素[15] 。研究區(qū)水化學(xué)組分相關(guān)矩陣KMO 檢驗(yàn)值為0.575(>0.5),Bartlett 球形度檢驗(yàn)統(tǒng)計(jì)量Sig = 0(<0.01),表明適合進(jìn)行主成分分析。利用SPSS27 軟件因子分析模塊,按照特征值大于1 的原則,采用最大方差法進(jìn)行成分旋轉(zhuǎn), 共提?。?個(gè)主成分, 累計(jì)方差貢獻(xiàn)率為89.048%(見(jiàn)表3),可以解釋研究區(qū)地下水水化學(xué)主要控制因素。
主成分F1 貢獻(xiàn)率為63.528%,主要載荷為TDS、Ca2+、Mg2+、Cl-、SO2-4 和NO-3 。其中:NO-3 是人類(lèi)活動(dòng)中農(nóng)業(yè)化肥施用、生活污水排放等的代表性離子;Cl- 除來(lái)自蒸發(fā)鹽巖風(fēng)化溶解外,同樣為生活污水排放的特征性離子,由表2 可知Cl- 與NO-3 顯著正相關(guān),說(shuō)明研究區(qū)地下水中Cl-受人類(lèi)活動(dòng)影響較大;SO2-4 主要來(lái)自含硫礦物并受研究區(qū)礦業(yè)活動(dòng)影響較大;Ca2+ 和Mg2+與SO2-4 顯著正相關(guān),且Mg2+與SO2-4 的相關(guān)性高于Ca2+與SO2-4 ,此處Ca2+和Mg2+可能主要來(lái)源于硫酸鹽礦物的溶解。綜上表明F1 表示研究區(qū)地下水水化學(xué)主要受人類(lèi)活動(dòng)及硫酸鹽礦物風(fēng)化溶解的影響。
主成分F2 貢獻(xiàn)率為14.777%,主要載荷為Na+和HCO-3 ,其次為Ca2+和Mg2+。Na+主要來(lái)源于鹽巖的溶解,同時(shí)受陽(yáng)離子交換作用影響;HCO-3 主要來(lái)源于碳酸鹽巖的風(fēng)化溶解;Ca2+、Mg2+和HCO-3 顯著正相關(guān),且Ca2+與HCO-3 的相關(guān)性高于Mg2+與HCO-3 ,此處Ca2+和Mg2+可能主要來(lái)源于碳酸鹽巖的風(fēng)化溶解。因此,F2 表示研究區(qū)地下水主要受碳酸鹽巖的風(fēng)化溶解及陽(yáng)離子交換作用的影響。
主成分F3 貢獻(xiàn)率為10.743%,主要載荷為K+。K+主要來(lái)源于硅酸鹽巖的風(fēng)化溶解及含鉀化肥的施用[ 24] 。由表2 可知K+ 與其他離子均不相關(guān),故推測(cè)F3 代表硅酸鹽巖的風(fēng)化溶解對(duì)地下水水化學(xué)組分的影響。
4 結(jié)論
研究區(qū)地下水為淡水,整體呈中性,TH 均值為427.79 mg/ L,TDS 均值為545.29 mg/ L。優(yōu)勢(shì)陰離子、陽(yáng)離子分別為HCO-3 、Ca2 +,陰離子含量表現(xiàn)為HCO-3 >SO2-4 >NO-3 >Cl-,陽(yáng)離子含量表現(xiàn)為Ca2 + >Na+ >Mg2 + >K+。地下水水化學(xué)類(lèi)型主要為HCO3 -Ca,其次為HCO3·SO4 -Ca、HCO3·SO4 -Ca·Mg、HCO3·Cl-Ca,水化學(xué)離子組成受多種水化學(xué)作用和人類(lèi)活動(dòng)的共同影響。
研究區(qū)地下水水化學(xué)特征以巖石風(fēng)化作用控制為主,水化學(xué)成分主要來(lái)自碳酸鹽巖和硅酸鹽巖的風(fēng)化溶解。地下水中Na+和K+ 主要受鹽巖溶解作用控制,同時(shí)受陽(yáng)離子交換作用影響;Ca2+ 和Mg2+ 主要受碳酸鹽巖和硫酸鹽巖溶解作用控制。人類(lèi)活動(dòng)對(duì)研究區(qū)地下水化學(xué)成分影響較大,主要與研究區(qū)礦業(yè)活動(dòng)有關(guān),且局部受農(nóng)業(yè)生產(chǎn)和生活污水排放的影響。政府管理部門(mén)應(yīng)加強(qiáng)環(huán)境綜合治理,加大礦山環(huán)境綜合整治力度,并對(duì)點(diǎn)源和面源污染進(jìn)行有效控制。
參考文獻(xiàn):
[1] 蔡小虎,康叢軒.錢(qián)塘江下游濱江地區(qū)地下水水化學(xué)特征與水質(zhì)評(píng)價(jià)[J].人民長(zhǎng)江,2023,54(6):27-33.
[2] 李建華,徐向東,劉鵬,等.定邊縣山前平原地下水水化學(xué)特征分析及水質(zhì)評(píng)價(jià)[J].人民黃河,2023,45(11):111-117.
[3] WERNER Aeschbach?Hertig,TOM Gleeson.Regional Strategies for the Accelerating Global Problem of Groundwater Depletion[J].Nature Geoscience,2012,5(1-4):853-861.
[4] 孫乃波,陳學(xué)群,李倩,等.威海市地下水水化學(xué)特征及成因分析[J].節(jié)水灌溉,2022(7):85-90.
[5] 李狀,蘇晶文,董長(zhǎng)春,等.安徽馬鞍山市當(dāng)涂地區(qū)地下水水化學(xué)特征及演化機(jī)制[J].中國(guó)地質(zhì),2022,49(5):1509-1526.
[6] HU Huiying,WEI Rongfei,ZERIZGHI Teklit,et al.Control Mecha?nisms of Water Chemistry Based on Long?Term Analyses of the Yangtze River [J]. Science of the Total Environment, 2023,892:164713.
[7] 韓朝輝,王郅睿,田輝,等.漢中盆地地下水水化學(xué)特征及其成因研究[J].西北地質(zhì),2023,56(4):263-273.
[8] 張子燕,伏永朋,王寧濤,等.丹江流域山區(qū)地表水-地下水水化學(xué)特征及其影響因素[J].環(huán)境科學(xué),2024,45(4):2107-2117.
[9] 易雅寧,孫曉懿,王富強(qiáng),等.三門(mén)峽庫(kù)區(qū)濕地水化學(xué)特征及影響因素分析[J].人民黃河,2021,43(3):90-96.
[10] 李曉剛,武洋洋,高巖,等.商丹盆地丹江一級(jí)階地黃土古土壤剖面碳酸鈣含量研究[J].甘肅科學(xué)學(xué)報(bào),2021,33(4):20-25.
[11] 王海燕,龐獎(jiǎng)勵(lì),黃春長(zhǎng),等.陜西商丹盆地茶房村黃土-古土壤剖面的粒度組成特征及意義[J].中山大學(xué)學(xué)報(bào)(自然科學(xué)版),2017,56(2):131-139.
[12] 張保祥,張超.水文地球化學(xué)方法在地下水研究中的應(yīng)用綜述[J].人民黃河,2019,41(10):135-142.
[13] 趙幸悅子,肖攀,宋文龍,等.贛南紅層地區(qū)地下水水化學(xué)特征及成因分析[J].科學(xué)技術(shù)與工程,2023,23(33):14112-14122.
[14] 成世才,董妍,宋永芬,等.濟(jì)南市先行區(qū)淺層地下水水文地球化學(xué)特征[J]. 人民黃河,2021,43(12):86 -90,99.
[15] 劉海,康博,管政亭,等.淮南煤礦區(qū)地表水和地下水水化學(xué)特征及控制因素[J].環(huán)境科學(xué),2023,44(11):6038-6049.
[16] GIBBS R J.Mechanisms Controlling World Water Chemistry[J].Science,1970,170(3962):1088-1090.
[17] THOMAS J,JOSEPH S,THRIVIKRAMJI K P.Hydrochemical Variations of a Tropical Mountain River System in a Rain Shadow Region of the Southern Western Ghats,Kerala,India [J].Applied Geochemistry,2015,63:456-471.
[18] 金陽(yáng),姜月華,董賢哲,等.浙江寧波平原地下水水化學(xué)特征及其生態(tài)環(huán)境效應(yīng)[J].中國(guó)地質(zhì),2022,49(5):1527-1542.
[19] ROUABHIA Abdelkader,DJABRI Larbi,HADJI Rihab,et al.Geochemical Characterization of Groundwater from Shallow Aquifer Surrounding Fetzara Lake N.E.Algeria[J].Arabian Journal of Geosciences,2012,5(1):1-13.
[20] 楊艷林,靖晶,齊信,等.贛州北部丘陵山區(qū)淺層地下水化學(xué)特征及成因分析[J].中國(guó)地質(zhì)調(diào)查,2023,10(6):60-68.
[21] 彭紅霞,侯清芹,曾敏,等.雷州半島地下水化學(xué)特征及控制因素分析[J].環(huán)境科學(xué),2021,42(11):5375-5383.
[22] 韓貴琳,劉叢強(qiáng).貴州喀斯特地區(qū)河流的研究:碳酸鹽巖溶解控制的水文地球化學(xué)特征[J]. 地球科學(xué)進(jìn)展,2005,20(4):394-406.
[23] 雷琨,何守陽(yáng),劉璽,等.巖溶易旱區(qū)地下水水文地球化學(xué)特征及其控制因素:以興仁縣為例[J].地球與環(huán)境,2016,44(6):619-627.
[24] 郭高軒,代垠東,許亮,等.北京西山巖溶地下水化學(xué)特征及成因分析[J].環(huán)境科學(xué),2024,45(2):802-812.
【責(zé)任編輯 呂艷梅】
基金項(xiàng)目:中國(guó)地質(zhì)調(diào)查局地質(zhì)調(diào)查項(xiàng)目( DD20230481,DD20242461)