于開寧，田 劍，劉景濤，周亞紅，呂曉立，Vetrimurugan Elumalai

(1.河北地質大學，河北省高校生態(tài)環(huán)境地質應用技術研發(fā)中心，河北石家莊 050031；2.核工業(yè)二○八大隊，內蒙古包頭 014010；3.中國地質科學院水文地質與環(huán)境地質研究所，河北石家莊 050061;4.祖魯蘭大學，南非德班999136)

0 引言

中國西部區(qū)域發(fā)展面臨的主要問題之一是水資源及水質問題(中國科學院，2018)。地下水與其存在的環(huán)境空間在自然地質因素作用于人類活動作用下所演化形成的狀態(tài)與過程統(tǒng)稱為地下水環(huán)境，地下水環(huán)境的變化對干旱區(qū)和半干旱區(qū)經(jīng)濟和生態(tài)穩(wěn)定具有舉足輕重的作用(沈照理等，1993)。地下水水化學類型是地下水化學成分的集中反映，對地下水化學的特征與演變規(guī)律進行研究，可以更好地揭示地下水與環(huán)境的相互作用機制，對維持半干旱區(qū)與干旱區(qū)的生態(tài)系統(tǒng)穩(wěn)定有著極其重要的現(xiàn)實意義。

近年來國內外學者采用數(shù)理統(tǒng)計、piper三線圖、Gibbs圖、離子比例系數(shù)、鈉吸附比、飽和指數(shù)及同位素示蹤、PHREEQC軟件反向模擬等方法對不同區(qū)域的地下水化學特征及演化模擬進行研究。陳晨等(2019)等利用描述性分析、piper三線圖、Gibbs圖、離子比例系數(shù)等方法研究了泰萊盆地地下水化學特征；姜凌(2009)利用同位素示蹤技術,論證了地下水的補給成因、徑流途徑和排泄方式，采用鈉吸附比、飽和指數(shù)等方法，分析研究區(qū)地下水水化學演化中發(fā)生的溶解沉淀作用；劉海等(2019)等采用PHREEQC軟件進行水文地球化學模擬，定量分析了泗縣地下水的形成機理及演化。前人對蘭州市的地下水研究多集中在“三氮”和重金屬污染方面，呂曉立等(2019)等對蘭州地下水“三氮”污染和鐵錳分布特征提出研究，并得出控制“三氮”和鐵錳離子遷移轉化及其空間分布的主要因素。學者們對蘭州市地下水化學特征及演化規(guī)律研究較少。尤其是對于地下水資源匱乏的西北地區(qū)，地下水化學成分的研究顯得格外重要。本文通過對蘭州市地下水化學成分的形成過程的分析，掌握地下水水化學演化規(guī)律，揭示地下水與環(huán)境相互影響關系，對蘭州市的地下水保護和利用具有重要意義，為西北干旱氣候河谷盆地類型城市的地下水可持續(xù)開發(fā)利用提供借鑒。

1 研究區(qū)概況

1.1 自然地理環(huán)境

蘭州市地處甘肅中部，全市下轄5個區(qū)、3個縣，總面積13100 km2，地理位置：北緯35°45′～36°28′，東經(jīng)102°49′～103°59′。蘭州氣候特點是干燥寒冷，日照多，光能潛力大，屬于半干旱氣候，多年平均氣溫8.9 ℃，多年平均降水量為317.6 mm，主要集中在7～9月，平均相對濕度56%，年蒸發(fā)量1600 mm左右。蘭州市河流均屬黃河干流水系，黃河自西向東穿過蘭州市區(qū)，蘭州工農業(yè)、生活用水主要取自地下水資源。

1.2 水文地質條件

蘭州屬于河谷盆地型城市，地勢總體南北高而中間低。南部為皋蘭山等黃土丘陵，中部為黃河谷地，北部為黃土丘陵及低中山(圖1)。蘭州斷陷盆地的基底由新近系泥巖、砂質泥巖和砂質礫巖組成，含有大量可溶鹽類，常見石膏、芒硝、菱鎂礦等晶體①(呂曉立等，2013)。蘭州市典型水文地質剖面見圖2。

研究區(qū)有供水意義含水層主要分布在黃河河谷平原松散巖類孔隙水，含水層巖性為砂礫卵石層，含水層富水性良好，區(qū)域地下水從基巖山地接受降水補給，分別向南北徑流，最終均排泄于黃河，構成一個完整的地下水循環(huán)系統(tǒng)，在天然狀態(tài)下，松散巖類孔隙水以接受大氣降水入滲、河流入滲和基巖裂隙水轉化補給為主，在徑流過程中接受農田灌溉水和城市污水的入滲補給，排泄方式主要為人工開采、蒸發(fā)和徑流排泄到黃河(黎濤等，2010；呂曉立等，2014)。

2 樣品采集與分析方法

2.1 樣品采集與分析

圖1 蘭州市地下水采樣點分布圖

圖2 研究區(qū)I-I’水文地質剖面

2.2 數(shù)據(jù)分析方法

本文數(shù)據(jù)可靠性利用電荷平衡誤差E進行分析，電荷平衡誤差E的絕對值小于5%為可靠數(shù)據(jù)，經(jīng)過檢驗，所有數(shù)據(jù)均為可靠性數(shù)據(jù)。

使用Excel對數(shù)據(jù)進行描述性統(tǒng)計分析和舒卡列夫分類，利用Aquachem繪制地下水Piper三線圖，應用SPSS19軟件計算相關性系數(shù)，采用Origin9.0繪制主要離子關系圖和Gibbs圖，借助PHREEQC軟件進行反向水文地球化學模擬(曾昭華，1999；張兆迪等，2017;藍天杉等，2018)。

3 結果與分析

3.1 描述性統(tǒng)計分析

表1 地下水化學主要指標統(tǒng)計結果

續(xù)表1

3.2 水化學類型特征及演化

圖3 研究區(qū)地下水化學成分Piper三線圖

研究區(qū)地下水從基巖山區(qū)接受降水補給，分別向南北徑流，最終均排泄于斷陷盆地中部。在黃河倆側地區(qū)，基巖山區(qū)地下水經(jīng)過長距離的徑流后形成的高礦化水與黃河側滲補給的地下水相混合,形成以黃河水體為對稱軸的南北向水化學水平分帶。在南北山前一帶，處于地下水補給區(qū)，降水量大，水質好，為重碳酸鹽水帶；在河谷潛水區(qū)，處于地下水徑流區(qū)，隨之徑流距離增加，水質漸差，為硫酸鹽水帶；在河谷盆地中央，處于地下水排泄區(qū)，為氯化物水帶。由于氣候干燥，且地下水位埋藏較淺，點X35和X34位于洪積扇前緣細顆粒相沉積層，滲透性能較差，地下水入滲、徑流緩慢，地下水受到強烈的蒸發(fā)濃縮作用和溶濾作用形成了高礦化度鹽水。這體現(xiàn)了干旱區(qū)河谷盆地典型的水文地球化學演化規(guī)律。

為分析不同水體水化學組分的演化特征及規(guī)律，將蘭州市地下水化學數(shù)據(jù)繪入Gibbs圖。由圖4a所示，大部分水樣點Cl-占的比重相對較大。從圖4b可以看出，絕大部分水樣點陽離子以Na+為主，礦化度較大。

圖4 研究區(qū)地下水水化學Gibbs分布圖

3.3 離子關系圖分析

地下水中可溶組分的比值關系是探討物質來源的有效指標，為了進一步分析區(qū)內地下水化學離子形成作用，通過繪制主要離子關系圖可以分析判斷地下水化學成因：

(1)γNa+/γCl-

圖5 地下水離子比例相關值

3.4 地下水組分相關性分析

水文地球化學研究中，常常采用Pearson相關性系數(shù)矩陣將地下水中各組分的相關性定量化表示出來(牛兆軒等，2019；姜體勝，2017;何興江等，2006)。本文采用SPSS19.0軟件對研究區(qū)水樣化學離子進行Pearson相關系數(shù)計算，計算結果見表2。

表2 研究區(qū)化學離子Pearson相關系數(shù)

3.5 水文地球化學模擬

水文地球化學模擬常用美國地質調查局開發(fā)的PHREEQC軟件。該軟件是一個用于計算多種條件下的水文地球化學反應模擬軟件，其計算過程主要根據(jù)熱力學模型和水質分析資料。確定當?shù)叵滤芙饨M分反應達到平衡時，計算地下水中各化學組分的分布，進而對地下水中礦物質飽和度進行計算。根據(jù)計算結果，可對礦物在給定水溶液環(huán)境下的溶解或沉淀趨勢進行判定(靳書賀等，2017；方展等，2017；武亞遵等，2018)。

在對研究區(qū)水文地質條件分析研究的基礎上，依據(jù)地下水流向，本次選取兩條路徑A、B作為水化學模擬路徑。“可能礦物相”的選取依據(jù)含水層巖性及礦物組成成分，選取“可能礦物相”為巖鹽、鉀長石、方解石、白云石、石膏、高嶺石、石英、鈣蒙脫石以及Ca2+、Na+陽離子交換。從表3分析得知，模擬路徑起始點至終止點礦物相飽和指數(shù)(SI)值變化中文石、方解石、白云石趨于飽和狀態(tài)，硬石膏、石膏、巖鹽則呈溶解趨勢(模擬過程中不確定度取值為0.05)。

表3 研究區(qū)各模擬點礦物相飽和指數(shù)

表4 研究區(qū)水-巖作用過程中礦物轉移量(mmol/L)

4 結論

(1)蘭州市地下水總體呈堿性，硬度較高，微咸，從南北山前到黃河兩側沿岸，沿地下水徑流方向，TDS逐漸減小。

(4)根據(jù)PHREEQC軟件反向模擬計算表明：方解石、白云石在整個模擬路徑上趨于飽和狀態(tài)，具有沉淀趨勢，硬石膏、石膏、巖鹽在整個模擬路徑上則呈非飽和狀態(tài)，具有溶解趨勢。在兩路徑模擬中，均為方解石發(fā)生溶解，Na+將Ca2+發(fā)生陽離子交替吸附作用，白云石、巖鹽、石膏均發(fā)生沉淀。

[注 釋]

① 楊顧,王華清,余振國,連炎清,王武衡.1991. 蘭州市傍河水源地地下水資源及污染研究報告[R]. 蘭州：甘肅省第一水文地質工程地質隊.