穆興平原北區淺層地下水水化學分布特征 ぜ捌湫緯苫理

崔小順　鄭昭賢

摘要：采用2016年8月采集的192個地下水樣品水化學數據，查明了穆棱河—興凱湖平原（穆興平原）北區第四系淺層地下水溶解性總固體（TDS）和地下水化學類型分布特征，并結合研究區水文地質條件和人類活動影響分析其形成原因。結果表明：穆興平原北區淺層地下水主要以TDS小于05 g/L的淡水為主，未見TDS大于10 g/L的微咸水。地下水TDS受含水介質、溶濾作用和地表水-地下水交互作用影響，在區內呈多種分帶規律。區內淺層地下水化學類型不具備明顯的水平分帶規律，主要以HCO

3Ca·Mg型水為主，與人為活動較少的20世紀60年代背景水化學類型（HCO3Ca·Na型水）對照，地下水中硫酸根、氯離子及鎂離子的增加與煤礦開采、生活用水污染及大范圍農業施肥等人類活動相關。

關鍵詞：溶解性總固體；水化學類型；水化學特征；跨界含水層；穆興平原

中圖分類號：P426.2文獻標志碼：A文章編號：

16721683（2018）04014608

Hydrochemical distribution characteristics of shallow groundwater in the north of

MulingXingkai Plain and their formation mechanism

CUI Xiaoshun1，2，ZHENG Zhaoxian3，4，CHENG Zhongshuang3，SU Chen3，LI Zhuang4，WU Yulong4，CHEN Zongyu3

（

1.College of Environment and Resources，Jilin University，Changchun 130021，China；

2.Institute of Water Resources and Environment，Jilin University，Changchun 130021，China；

3.Institute of Hydrogeology and Environmental Geology，Chinese Academy of Geological Science，

Shijiazhuang 050061，China；4.China University of Geosciences （Beijing），Beijing 100083，China）

Abstract：

We studied the distribution characteristics of total dissolved solids （TDS） and hydrochemical types of shallow groundwater in the MulingXingkai Plain on the basis of analyzing 192 groundwater samples.Then we explored the cause of formation by considering the hydrogeological conditions and human activities.The results showed that the shallow groundwater in the study area was mainly fresh water with TDS less than 05 g/L，and there was no saline groundwater with TDS higher than 10 g/L.Influenced by aquifer medium，leaching，and groundwaterriver interaction，the groundwater TDS showed zonal distribution in the study area.There was no obvious zonal distribution regularity of the hydrochemical types in the study area.Nowadays，the HCO3Ca·Mg type water has replaced the HCO3Ca·Na type water in the 1960 s （considered as background water type） to become the common water type in the study area.The increases of sulfate ion，chloride ion，and magnesium ion in groundwater were mainly due to coal mining，domestic wastewater pollution，and fertilization.

Key words：

total dissolved solids；hydrochemical type；hydrochemical characteristic；transboundary aquifer；MulingXingkai Plain

穆棱河-興凱湖平原（以下簡稱穆興平原）位于我國東北三江平原南部，與俄羅斯接壤，是我國重要的商品糧食基地，也是大量珍惜瀕危保護物種的棲息地，其由沼澤濕地為主，陸地水域、島狀林、草甸為輔的綜合生態體系對整個東亞-遠東地區的生態系統具有十分重要的調控作用[1]。近年來，隨著我國“一帶一路”戰略方針的加快實施，地處“中俄蒙經濟走廊帶”重要位置的穆興平原，因其重要的經濟定位和生態功能受到廣泛關注。第四系含水層作為穆興平原主要供水水源，在穆興平原經濟社會發展及生態體系建設中起著舉足輕重的作用。隨著穆興平原城市化進程的加快推進與現代化農業程度的不斷提高，地下水過量開采，農藥、化肥的大量使用，引發了地下水降落漏斗、地下水污染等一系列環境與生態問題，導致地下水環境惡化，嚴重制約了穆興平原經濟、社會發展，并影響其生態承載力和功能[23]。不僅如此，穆興平原第四系含水層作為中俄跨界含水層，涉及國家之間利益關系，地下水資源分配、地下水污染問題有可能引發國際糾紛[4]，有必要盡快開展穆興平原地下水環境調查與相關研究。淺層地下水化學分布特征，是研究區域地下水循環和水化學演化的主要手段，也是地下水資源評價的重要內容，可為科學管理、合理開發、公平利用跨界含水層提供指導意見，使跨國界地下水資源可持續利用，對維護國家主權、加強國家間睦鄰友好關系具有重要的指導意義。然而，目前關于穆興平原淺層地下水化學分布特征的研究較少，除1962年穆興平原綜合水文地質測量外[5]，未見其它相關研究。對此，本次研究利用2016年所采集的192個第四系淺層地下水樣品的水化學數據，分析穆興平原兩大地下水系統之一的穆棱河北區系統地下水化學特征，并結合研究區水文地質條件進一步探討其成因。

1研究區概況

穆興平原北區位于黑龍江省東部，與俄羅斯接壤，區內包括穆棱河、七虎林河、阿布沁河及穆棱河沖積平原。本區東起烏蘇里江、西到張廣才嶺、北依完達山，南瀕穆棱河，總面積5 252 km2，行政區跨越虎林市、密山市[6]。本區屬中緯度寒溫帶濕潤、半濕潤大陸性季風氣候，冬季漫長，嚴寒而干燥，夏季短促，溫暖且濕潤，年平均溫度29～31 ℃。區內雨量充沛，多年平均降水量為526～710 mm，6月至9月占全年總量70%；蒸發量為降水量15～2倍，5至9月占全年總量70%[7]。

穆興平原北區地形總體由西向東略微傾斜，坡降1/6 000～1/10 000，地勢低平，微有起伏，僅在虎林和虎頭散布著被地塹切割殘余的孤山。平原區分布著廣泛的第四系地層，其中第四系含水層自山前至平原區主要分為中、下更新統坡積、洪積層和上更新統、全新統沖積層。坡積、洪積層中的潛水主要分布在山前崗阜狀緩坡底部，主要含水層為黏土加碎石，一般厚度為10～20 m，地下水埋藏深度約為10～20 m，地下水位呈顯著季節性動態特征；沖積層中的潛水主要分布在河漫灘及一級階地之下，主要含水層為細、中、粗砂及砂礫石，其厚度由山前地區1～35 m，向平原區逐漸增厚到60～150 m，地下水埋藏深度通常在1～3 m，地下水位動態主要受大氣降水和地表水控制，呈季節性動態特征。第四系淺層地下水主要接受大氣降水的入滲補給、山區和丘陵區的基巖裂隙水的側向徑流補給和水稻田的灌溉入滲補給。地下水總體由西、西北向東、東南流動至區內最主要的兩條河流，即烏蘇里江和穆棱河（圖1）。區內淺層地下水主要排泄方式以蒸散發和人工開采（以農業用水為主）為主，其次向地表水排泄，僅有極少量越流排泄至深層地下水。

穆興平原第四系含水層在區內分布廣泛，水量豐富，水質較好，埋藏淺，是區內城鎮及工農業的主要供水水源。從2011年虎林市、密山市及其轄區地下水資源開采結構及開采量資料可知[8]，城鎮居民用水量為2 193萬m3、第一產業用水量為133 9184萬m3、第二產業用水量為1 5947萬m3、第三產業用水量為4255萬m3、生態環境用水量為94萬m3，其中農田灌溉用水量占第一產業用水量976%。由此可見，地下水開采中，農田灌溉用水量所占比重最大，這與研究區以農業用地為主的土地利用類型相符。

2樣品采集和測試

本次研究參考《地下水環境監測技術規范》（HJ/T 164-2004）[9]中有關地下水質控制監測井的布置原則和要求，根據區內淺層地下水的流向，考慮主要河流與地下水的潛在交互作用，以有限元網格剖分方法，在研究區不同地貌單元及第四系地層布置了192個樣品取樣點（圖1）。地下水取樣時間為2016年8月份（豐水期），即區內廣泛分布的農田已經完成施肥和除草除蟲施藥工作。采集樣品均為第四系淺層地下水，主要取自農田灌溉機井和鄉鎮居民供水民井。

所有水樣在取樣之前均進行洗井作業，洗井抽水量達到3～5倍井筒體積時，現場測量地下水質參數，包括溶解氧（DO）、氧化還原電位（Eh）、酸堿度（pH）、電導率（EC）、溫度，當上述參數穩定后進行取樣。采集的水樣由聚乙烯瓶儲存，并冷藏保存，并于一周內送達實驗室進行相關測試前處理，以保證測試數據的準確性。

水樣的室內測試分析項目包括K+、Na+、Ca2+、Mg2+、HCO3、SO42、Cl、NO3、δDH2O、δ18OH2O，測試分析方法參照《水和廢水監測分析方法（第四版）》[10]和《同位素地質樣品分析方法總則及一般規定》[11]。水化學測試分析結果采用陰陽離子平衡的方法進行可靠性檢查，相對誤差均小于±5%，氫、氧同位素測試分析相對誤差分別為±01‰、±0025‰，確保本次研究中的水化學和同位素數據準確可靠。

3結果與討論

3.1地下水溶解性總固體分布特征

溶解性總固體（TDS）是地下水中主要離子的集中體現，也是研究地下水化學特征的重要指標[12]。從淺層地下水TDS分析測試結果中可知，穆興平原北區廣泛分布著01

從研究區淺層地下水TDS分布特征來看，其空間變化較大，但總體上具有一定的規律，即自山前崗阜狀緩坡至河成一級階地，以地貌單元為界，沿區域地下水流向，淺層地下水TDS逐漸降低，與我國北方盆地中常見的地下水TDS沿徑流路徑逐漸增加的一般規律相反，僅在七虎林河與穆棱河河間地帶，自八五零農場至虎頭沿地下水流路徑TDS呈增加趨勢。經分析，上述地下水TDS分布特征的成因主要有以下幾方面：

（1）含水層巖性。

崗阜狀緩坡處于地形上升和下降的過渡地帶，此處含水層以坡積物和洪積物為主，由亞黏土、砂、亞黏土夾碎石組成，含水層透水性較差，地下水流動緩慢，徑流相對較弱[5]。因此，降水垂直入滲對地下水補給有限的條件下，以山區基巖裂隙水側向徑流補給為主的地下水在此處與含水層中礦物組分產生長時間水巖交互作用，使含水介質中易溶巖類不斷溶濾進入地下水，導致地下水TDS上升。地下水徑流至河成二級階地，因含水層中亞黏土含量相比崗阜狀緩坡較少，透水能力提高，大氣降水垂向入滲在地下水補給中的比重增加，稀釋了上游徑流而來的高TDS地下水，使地下水TDS降低。當地下水徑流至區內廣泛分布的河成一級階地，因其含水層主要由沖積砂、砂礫石組成，透水蓄水性良好，加之地形平坦，為大量接收大氣降水和地表水入滲補給提供了有利條件[5]，使徑流補給至這一地貌單元中地下水進一步被降水和地表水稀釋，在山前崗阜狀緩坡至一級階地形成TDS趨減分帶。

（2）溶濾作用。

影響水化學場形成的主要作用包括溶濾作用、陽離子交替吸附作用、蒸發濃縮作用和混合作用。Gibbs通過對全球典型地表水化學組分的分析，將天然水組分的控制因素分為三個類型：大氣降水控制型、水巖作用控制型和蒸發沉淀控制型，并利用TDSCl/（Cl+HCO3）關系圖和TDSNa/（Na+Ca）關系圖，劃分三種主要控制因素下天然水體的特征區域[1314]。Gibbs圖同樣適用于人為影響較少的地下水化學成因及控制因素的研究[15]。本次研究根據工作區地下水中TDS分帶規律，將潛水分為TDS低值區（TDS<02 g/L）、TDS中值區（02≤TDS<03 g/L）和TDS高值區（TDS≥03 g/L），分別分析地下水在不同區域發生的水化學作用。從淺層地下水Gibbs圖（圖3）中可知，整個研究區潛水在不同TDS分區主要受水巖作用控制。因研究區淺層地下水中TDS均較低，地下水中所有礦物飽和指數均小于零，即處于礦物可溶狀態，可判定研究區TDS主要受溶濾作用影響。

c（Mg/Na）與c（Ca/Na）相關性可用于評估不同類型礦物在地下水溶濾作用中的貢獻量[1617]。從圖2、圖4可知，淺層地下水在穆棱河與七虎林河河間地帶（TDS低值區和TDS中值區）自八五零農場至

虎頭沿線形成較長的地下徑流，在徑流過程中產生以硅酸鹽礦物溶解為主、碳酸鹽礦物溶解為輔的溶濾作用，使地下水TDS沿地下水流向緩慢增加。慶豐農場一帶地下水開采量較大，在此區形成了兩處較為明顯的地下水降落漏斗，而且TDS高值區與漏斗中心基本重疊，在含水介質并無明顯變化的情況下，表明地下水流經此處，由漏斗邊緣向中心匯聚，形成了環狀溶濾加強區，將漏斗區周邊含水介質中的可溶鹽分聚集到中心區，導致地下水TDS在漏斗中心升高。淺層地下水在七虎林河與阿布沁河河間地帶自迎春鎮至烏蘇里江同樣形成地下潛流，但在徑流下游區未出現TDS大于02 g/L的中值區，表明以難溶的硅酸鹽礦物風化溶解為主的地下水化學作用，在較短的徑流路徑上，對地下水TDS影響有限。而在崗阜狀緩坡（TDS高值區）地下水中Mg/Na（中值為142）和Ca/Na（中值為076）比值相比TDS中值區（Mg/Na中值為135，Ca/Na中值為072）較高，表明地下水溶濾作用中碳酸鹽溶解比重有所增加。當研究區不同地貌第四系含水介質物源相似的情況下，這也從另一側面反映徑流條件較好的一級階地和河漫灘含水介質中碳酸鹽巖礦物在長期的地下徑流過程中已基本溶解殆盡，現階段含水介質以硅酸鹽礦物為主，而徑流條件相對較差的崗阜狀緩坡及二級階地含水介質中還保留了較多相對易溶的碳酸鹽巖礦物。

（3）地表水與地下水交互作用。據穆興平原地下水資源勘察資料[18]，研究區地表水與地下水的交互作用較為強烈，尤其河流兩岸漫灘砂崗地段，因其含水層均由無亞粘土夾層的砂卵礫石組成，在豐水期，河水可大量入滲補給地下水。從阿布沁河河水與地下水化學指標對比統計（表1）可知，阿布沁河不同河段δ2DH2O、δ18OH2O、Ca/Mg、pH等同位素組成和水化學特征值與相應河段漫灘地下水相近，而且自河流上游至下游均呈現類似的“親緣性”，根據Kehew等總結的濕地與地下水交互模型[19]，在豐水期主要河流沿線，形成以河水入滲補給地下水為主的水質交互帶，使地下水TDS在河流兩岸呈帶狀低值區。

3.2地下水化學類型分布特征

地下水的化學成分是地下水與環境—自然地理、地質背景以及人類活動長期相互作用的產物[20]。區域地下水化學類型可綜合反映地下水賦存條件、水巖相互作用、補排循環特征等原生水文地球化學環境和地下水超采、地下水污染等人為活動影響[21]。為了解穆興平原北區淺層地下水化學類型及空間分布規律，采用舒卡列夫分類法對區域地下水化學類型進行劃分。由于研究區地下水化學分類較多，為便于分析，對陰離子和陽離子分別確定其水化學類型。

從淺層地下水陰離子水化學類型分區（圖5）可知，研究區淺層地下水主要以HCO3型水為主，在本區廣泛分布，其分布面積為4 823 km2，占研究區總面[CM（22]積的918%。其它水化學類型在研究區內零星分[CM）]

從研究區淺層地下水陰離子水化學類型分布特征來看，其空間分布變化較小，不具備水平分帶性，但其分布也具有一定的規律性。在崗阜狀緩坡及河成二級階地，由于地下水主要產生硅酸鹽礦物風化溶解和碳酸鹽礦物溶解作用（圖4），故以重碳酸型水為主。崗伏狀緩坡及合成二級階地部分地區，因居民點居多，人為污染風險加劇，其中密山以西緩坡地下水中氯離子和硫酸根含量相比相同地貌其余地區較高，而該區含水層中并無相關化學組成的蒸發鹽巖，判定該區為典型人為污染區，地下水類型由HCO3型水轉變為HCO3·Cl·SO4型水。在河成一級階地及河漫灘，含水層為透水良好的砂和礫石組成，因大氣降水可入滲補給此區的地下水，地下水從補給到排泄路徑不長，而且地下水以硅酸鹽風化溶解作用為主（圖4），水化學成分仍以HCO3型水為主，與上游區緩坡地水化學類型變化不大。但在八五零農場以東形成帶狀分布的HCO3·SO4型水，其原因是此帶上游區煤田中硫化物在空氣氧化和降水淋濾作用下，增加了地下水中硫酸根含量。除此之外，在虎頭鎮和云山農場局部點狀分布SO4·Cl型水，這可能與村鎮集中垃圾填埋場相關，各種污染物經過大氣降水的溶解、淋濾、下滲等作用進入地下水，導致地下水中氯離子和硫酸根含量增高，進而完全改變天然地下水化學類型。

[JP+1]從淺層地下水陽離子水化學類型分區（圖6）可知，研究區淺層地下水主要以Ca·Mg型水為主，在本區廣泛分布，其分布面積為4 116 km2，占研究區總面積的784%。其它水化學類型在研究區內零星分布，Ca型水主要在張廣才嶺山前分布，其分布面積為456 km2，占研究區總面積的87%；Ca·Mg·Na型水、Ca·Na型水主要在東方紅鎮、八五零農場、虎頭局部呈小型面狀分布，其分布面積分別為315 km2、316 km2、，占研究區總面積的60%、60%。

據1962年穆興平原水文地質測繪資料[5]，現狀陽離子水化學類型相比人類活動較少的20世紀60年代背景地下水化學類型有明顯變化，原本區內廣泛分布的Ca·Na型水轉變為Ca·Mg型水，這與近幾十年快速發展的農業不無關系。據調查，除烏蘇里江沿線濕地自然保護區之外，區內廣泛分布的農田均以鎂肥改善土質、增加產量。鎂肥中的鎂離子經過大氣降水的溶解、淋濾等作用進入地下水，導致區內地下水化學類型產生較大改變。山前帶狀分布的Ca型水主要反映了大氣降水中陽離子水化學特征。除虎頭之外，Ca·Mg·Na型水、Ca·Na型水面狀分布區地下水TDS與周邊Ca·Mg型水相似，也未見明顯的陰離子水化學類型的變化，表明此區地下水未見明顯人為污染，分析其鈉離子毫克當量百分比的上升與此區普遍分布的火成巖風化亞黏土夾層有關，此類亞黏土中含有較多的含鈉鹽礦物，可在含水層中產生以鈉離子為主為硅酸鹽淋濾作用。[（Ca+Mg）（HCO3+SO4）]與[（Na+K）Cl]的相關性可用于評估含水層中產生的陽離子交替吸附作用[2223]。從圖7可知，TDS低、中值區大部分水點均落在y=-x線上或附近，反映含水層中發生陽離子交替吸附作用，其中Ca·Mg·Na型水、Ca·Na型水不僅均落在y=-x線上，而且陽離子交換量較高，表明此類地下水中發生強烈的陽離子交替吸附作用，使吸附于黏土礦物中的鈉離子進入地下水，進而與淋濾作用共同影響地下水化學類型。與之相反，虎頭鎮周邊地下水TDS相比周邊地下水明顯上升，地下水以SO4·Cl型水為主，表明此區地下水中鈉離子的升高主要受到人為污染影響。

4結論

（1）穆興平原北區廣泛分布著TDS小于05 g/L的淺層地下水，未見TDS大于10 g/L的微咸水。地下水TDS分布特征可歸納為：自山前崗阜狀緩坡至合成一級階地邊緣，沿區域地下水流向，以地貌單元為界，淺層地下水TDS逐漸降低；在廣泛分布的河成一級階地，沿地下水流向TDS逐漸升高，在降落漏斗周邊出現小型TDS面狀高值區（TDS>03 g/L）；河流與地下水水質交互作用明顯，豐水期以河水補給地下水為主，在河流沿線形成條帶狀TDS低值區（TDS<02 g/L）。

（2）穆興平原北區淺層地下水化學類型不具備明顯的水平分帶規律。本區含水介質主要以硅酸鹽礦物和碳酸鹽礦物組成，使全區地下水主要以HCO3型水為主，僅在局部受強烈人為活動影響地區，形成帶狀或點狀分布的HCO3·SO4型水、HCO3·Cl型水、HCO3·Cl·SO4型水和SO4·Cl型水。本區受長期農業灌溉影響，全區地下水以Ca·Mg型水為主，而在局部受含水介質礦物組分影響，出現小型面狀分布的Ca·Mg·Na型水、Ca·Na型水。虎頭附近地下水受人為污染影響，地下水中鈉離子含量增高，導致在該區形成小型面狀Ca·Na型水。

（3）本文重點論述了穆興平原TDS和水化學類型空間分布特征，并對其成因進行了深入分析，但人為活動對地下水化學特征的影響有待進行進一步探討。隨著工作程度的不斷推進，今后將開展典型人為污染物的示蹤識別研究，進一步確定人為活動對地下水化學組分的影響，地下水污染識別及防控提供理論依據。

參考文獻（References）：

[1]魏強，佟連軍，楊麗華，等.三江平原濕地生態系統生物多樣性保護價值[J].生態學報，2013，34（3）：375380.（WEI Q，TONG L K，YANG L H，et al.The biodiversity conservation value of wetland ecosystem in the Sanjiang Plain[J].Acta Ecologica Sinica，2015，35（4）：935943.（in Chinese）） DOI：10.5846/stxb201305301241.

[2]楊湘奎，孔慶軒，李曉抗.三江平原地下水資源合理開發利用模式探討[J].水文地質工程地質，2006，33（3）：4952.（YANG X K，KONG Q X，LI X K.Discussion of reasonable exploration and use pattern of groundwater resources in the Sanjiang Plain[J].Hydrogeology & Engineering Geology，2006，33（3）：4952.（in Chinese）） DOI：10.16030/j.cnki.issn.10003665.2006.03.012.

[3]王勇，柏鈺春，尹喜霖，等.三江平原生態地質環境分區研究[J].水文地質工程地質，2004，31 （6）：1118.（YANG W，BAI Y C，YIN X L，et al.Assessment and zoning of ecologicgeologic environment of Sanjiang Plain in Heilongjiang Province[J].Hydrogeology & Engineering Geology，2004，31 （6）：1118.（in Chinese）） DOI：10.16030/j.cnki.issn.10003665.2004.06.003.

[4]李堯.基于亞洲地下水資源圖的跨界含水層研究[D].北京：中國地質大學（北京），2013.（LI Y.Research on transboundary aquifers based on maps of groundwater resources in Asia[D].Beijing：China University of Geosciences （Beijing） （in Chinese））

[5]黑龍江省地質局水文地質工程地質大隊.第三章 水文地質[R].黑龍江省興凱湖平原1/20萬區域地質—綜合水文地質測量報告書，1962.（Hydrogeology and engineering geology team of Heilongjiang province geology survey.Chapter 3 Hydrogeology[R].Report of regional geology and comprehensive hydrogeology mapping of Xingkaihu plain in Heilongjiang province，1962.（in Chinese））

[6]袁立海，歐陽忠，呂緒東，等.穆棱河興凱湖平原稀土分析及稀土微肥應用的研究[J].河南大學學報，1990（1）：8186.（YUAN L H，OUYANG Z，LV X D，et al.Analysis of rare earth elements and application of rareearth fertilizers in MulingXingkai Plain[J].Journal of Henan University，1990（1）：8186.（in Chinese）） DOI：10.15991/j.cnki.411100.1990.01.017.

[7]劉芳芳.興凱湖濕地自然保護區生態質量評價[D].哈爾濱：東北林業大學，2009.（LIU F F.Ecological quality evaluation of Xingkai Lake wetland nature reserve[D].Haerbin：Northeast Forestry University （in Chinese））

[8]黑龍江省地質調查研究總院.第三章 區域水文地質概況[R].三江平原地下水污染調查評價成果報告，2015.（Institute of geological survey of Heilongjiang province.Chapter 3 Hydrogeology[R].Report of groundwater contamination evaluation of Sanjiang Plain，Northeast China，2015 （in Chinese））

[9]HJ/T 1642004，地下水環境監測技術規范[S].（HJ/T 1642004，Technical specifications for environmental monitoring of groundwater[S].（in Chinese））

[10]國家環境保護總局《水和廢水監測分析方法》編委會.水和廢水監測分析方法[M].北京：中國環境科學出版社，2002.（State Environmental Protection Administration “Water and wastewater monitoring and analysis method” Editorial Board.Water and wastewater monitoring and analysis method[M].Beijing：Chian Environmental Science Press，2002.（in Chinese））

[11]DZ/T 0184.1～0184.22―1997同位素地質樣品分析方法[S].（DZ/T 0184.1～0184.22―1997，Analysis method of isotopic geology sample[S].（in Chinese））

[12]張巖，董偉紅，李滿洲，等.河南平原淺層地下水總溶解固體和水化學類型的分布特征[J].水文，2011，31（2）：7983，68.（ZHANG Y，DONG W H，LI M Z.Study on distribution characteristics of TDS and hydrochemiscal type of shallow groundwatr in Henan plain[J].Journal of Chian Hydrology，2011，31（2）：7983，68.（in Chinese）） DOI：10.3969/j.issn.10000852.2011.02.019.

[13]GIBBS R J.Mechanisms controlling world water chemistry[J].Science，1970，170：10881090.DOI：10.1126/science.170.3962.1088.

[14]WANTY R B，VERPLANCK P L，SANJUAN C A，et al.Geochemistry of surface water in Alpine catchments in central Colorado，USA：resolving hostrock effects at different spatial scales[J].Apllied geochemistry，2009，24（4）：600610.DOI：10.1016/j.apgeochem.2008.12.012.

[15]郭華明，郭琦，賈永峰等.中國不同區域高砷地下水化學特征及形成過程[J].地球科學與環境學報，2013，35（3）：8396.（GUO H M，GUO Q，JIA Y F，et al.Chemical characteristics and geochemiscal processes of high arsenic groundwater in different regions of China[J].Journal of Earth Sciences and Environment，2013，35（3）：8396.（in Chinese）） DOI：10.3969/j.issn.16726561.2013.03.008.

[16]MUKHERJEE A，SCANLON B R，FRYAR A E，et al.Soulute chemistry and arsenic fate in aquifers between the Himalayan foothills and Indian craton （including central Gangetic plain）：influence of geology and geomorphology[J].Geochimica et Cosmochimica，2012，90：283302.DOI：10.1016/j.gca.2012.05.015.

[17]GRAILLARDET J，DUPRE B，LOUVAT P，et al.Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers[J].Chemical Geology，1999，159：330.DOI：10.1016/S00092541（99）000315.

[18]黑龍江省地質局水文地質工程地質大隊.第二章 水文地質[R].黑龍江省興凱湖平原一比二十萬地下水資源勘查報告書，1962.（Hydrogeology and engineering geology team of Heilongjiang province geology survey.Chapter 2 Hydrogeology[R].Exploration report of groundwater resources of Xingkaihu plain in Heilongjiang province，1962.（in Chinese））

[19]KEHEW A E，PASSERO R N，KRISHANAMURITHY R V，et al.Hydrochemical interaction between a wetland and an unconfined glacial drift aquifer，Southwestern Michigan[J].Ground water，1998，36（5）：849856.DOI：10.1111/j.17456584.1998.tb02204.x.

[20]沈照理.水文地球化學基礎[M].北京：地質出版社，1990.（SHEN Z L.Hydrogeochemistry Basis[M].Beijing：Geological Press，1990.（in Chinese））

[21]寇文杰.地下水化學分類方法的思考[J].西部資源，2012，（5）：108109.（KOU W J.The thinking of groundwater chemical classification method[J].Western Resource，2012，（5）：108109.（in Chinese））

[22]GARCIAL M G，DEL H M，BLESA M A.Geochemistry of groundwater in the alluvial plain of Tucuman province Argentina[J].Journal of Hydrology，2001，9：597610.DOI：10.1007/s1004000101664.

[23]王曉艷，李忠勤，蔣纏文.天山哈密榆樹溝流域地下水化學特征及其來源[J].干旱區地理，2017，40（3）：313321.（WANG X Y，LI Z Q，JIANG C W.Chemical characteristic and sources of the groundwater in Yushugou River Basin，Hami，Tianshan Mountains[J].Arid Land Geography，2017，40（3）：313321.（in Chinese）） DOI：10.13826/j.cnki.cn651103/x.2017.02.009.