全 文 ：第１２卷　第５期
２０１４年１０月
南 水 北 调 与 水 利 科 技
Ｓｏｕｔｈ－ｔｏ－Ｎｏｒｔｈ　Ｗａｔｅｒ　Ｔｒａｎｓｆｅｒｓ　ａｎｄ　Ｗａｔｅｒ　Ｓｃｉｅｎｃｅ　＆Ｔｅｃｈｎｏｌｏｇｙ
Ｖｏｌ．１２ Ｎｏ．５
Ｏｃｔ．２０１４
水文地质与工程地质
收稿日期：２０１３－０８－１９　　修回日期：２０１４－０７－１９　　网络出版时间：２０１４－０８－２７
网络出版地址：ｈｔｔｐ：／／ｗｗｗ．ｃｎｋｉ．ｎｅｔ／ｋｃｍｓ／ｄｏｉ／１０．１３４７６／ｊ．ｃｎｋｉ．ｎｓｂｄｑｋ．２０１４．０５．０３０．ｈｔｍｌ
基金项目：国家重点基础研究计划项目（９７３项目）（２０１０ＣＢ４２８８０５）；国家自然科学基金项目（４１１７２２１２）
作者简介：周亚红（１９８１－），女，河北石家庄人，讲师，博士研究生，主要从事地下水污染及治理方面的研究。Ｅ－ｍａｉｌ：ｚｈｙｈ３２７＠１６３．ｃｏｍ
ＤＯＩ：１０．１３４７６／ｊ．ｃｎｋｉ．ｎｓｂｄｑｋ．２０１４．０５．０３０
山东省梭罗树水源地地下水化学特征及形成机理
周亚红１，２，郝凌霄３，李俊峰４
（１．长安大学 环境科学与工程学院，西安７１００６４；２．石家庄经济学院 水资源与环境学院，石家庄０５００３１；
３．河北水务集团，石家庄０５０００１１；４．河北水文工程地质勘察院，石家庄０５００３１）
摘要：地下水中化学组分的种类、分布在一定程度上可以揭示地下水与环境的相互作用机制。通过对山东省梭罗树
水源地的２４组地下水样本中主要离子成分（ＨＣＯ３－、Ｃｌ－、ＳＯ４２－、Ｋ＋＋Ｎａ＋、Ｃａ２＋、Ｍｇ２＋）进行聚类分析（ＨＣＡ），
将水样分成Ｃ１、Ｃ２ 及Ｃ３ 三簇，再利用ＰＣＡ对地下水的８种化学指标进行降维处理，将其转化为ＰＣ１、ＰＣ２ 及ＰＣ２
三个综合指标，然后综合ＰＣＡ、Ｐｉｐｅｒ图的分析结果，以及各水样点的分布位置、围岩类型等条件对三簇水样点的水
化学特征及形成机理进行了分析。结果显示，Ｃ１ 簇水样为松散岩类孔隙水，ＰＣ１ 和ＰＣ３ 指标反映其主要受溶解作
用、淋滤作用及离子交换作用影响，水化学类型以 ＨＣＯ３·Ｃｌ－Ｃａ·Ｍｇ型或Ｃｌ·ＨＣＯ３－Ｃａ·Ｍｇ型为主；Ｃ２ 簇水样
为构造裂隙水，ＰＣ１ 和ＰＣ２ 指标反映其中的溶解作用及淋滤作用比Ｃ１ 弱，大部分离子成分含量较Ｃ１ 低，但围岩中
碳酸盐矿物、含镁矿物多，使得ＣＯ－３ 及 Ｍｇ２＋有不同程度的增加，水化学类型以ＨＣＯ３－Ｍｇ·Ｃａ型或者ＨＣＯ３－Ｃａ型
为主；Ｃ３ 簇水样为风化裂隙水，ＰＣ１ 和ＰＣ２ 指标反映其中溶解作用、淋滤作用更弱，水质相对较好，但季节性变化较
大，同时由于岩石中含有可溶性氟，使得Ｆ－有一定程度的增加，水化学类型以 ＨＣＯ３·ＳＯ４－Ｃａ·Ｍｇ型为主。
关键词：地下水化学；ＨＣＡ；ＰＣＡ；溶解作用；离子交换
中图分类号：Ｐ６４１．３　　文献标志码：Ａ　　文章编号：１６７２－１６８３（２０１４）０５－０１３３－０５
Ｈｙｄｒｏｃｈｅｍｉｃａｌ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ａｎｄ　ｆｏｒｍａｔｉｏｎ　ｍｅｃｈａｎｉｓｍ　ｏｆ
ｇｒｏｕｎｄｗａｔｅｒ　ｉｎ　Ｓｕｏｌｕｏｓｈｕ　ｇｒｏｕｎｄｗａｔｅｒ　ｓｏｕｒｃｅ　ａｒｅａ　ｏｆ　Ｓｈａｎｄｏｎｇ　ｐｒｏｖｉｎｃｅ
ＺＨＯＵ　Ｙａ－ｈｏｎｇ１，２，ＨＡＯ　Ｌｉｎｇ－ｘｉａｏ３，ＬＩ　Ｊｕｎ－ｆｅｎｇ４
（１．Ｓｃｈｏｏｌ　ｏｆ　Ｅｎｖｉｒｏｎｍｅｎｔａｌ　Ｓｃｉｅｎｃｅ　ａｎｄ　Ｅｎｇｉｎｅｅｒｉｎｇ，Ｃｈａｎｇ′ａｎ　Ｕｎｉｖｅｒｓｉｔｙ，Ｘｉ′ａｎ　７１００６４，Ｃｈｉｎａ；２．Ｓｃｈｏｏｌ　ｏｆ　Ｗａｔｅｒ
Ｒｅｓｏｕｒｃｅｓ　ａｎｄ　Ｅｎｖｉｒｏｎｍｅｎｔ，Ｓｈｉｊｉａｚｈｕａｎｇ　Ｕｎｉｖｅｒｓｉｔｙ　ｏｆ　Ｅｃｏｎｏｍｉｃｓ，Ｓｈｉｊｉａｚｈｕａｎｇ０５００３１，Ｃｈｉｎａ；３．Ｈｅｂｅｉ　Ｗａｔｅｒ　Ｇｒｏｕｐ，
Ｓｈｉｊｉａｚｈｕａｎｇ０５００１１，Ｃｈｉｎａ；４．Ｈｅｂｅｉ　Ｈｙｄｒｏｇｅｏｌｏｇｉｃａｌ　ａｎｄ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｅｘｐｌｏｒａｔｉｏｎ　Ｉｎｓｔｉｔｕｔｅ，Ｓｈｉｊｉａｚｈｕａｎｇ０５００３１，Ｃｈｉｎａ）
Ａｂｓｔｒａｃｔ：Ｃｈｅｍｉｃａｌ　ｃｏｍｐｏｓｉｔｉｏｎｓ　ａｎｄ　ｄｉｓｔｒｉｂｕｔｉｏｎ　ｉｎ　ｇｒｏｕｎｄｗａｔｅｒ　ｃａｎ　ｒｅｖｅａｌ　ｔｈｅ　ｉｎｔｅｒａｃｔｉｏｎ　ｍｅｃｈａｎｉｓｍ　ｂｅｔｗｅｅｎ　ｇｒｏｕｎｄｗａｔｅｒ　ａｎｄ
ｔｈｅ　ｅｎｖｉｒｏｎｍｅｎｔ　ｔｏ　ｓｏｍｅ　ｅｘｔｅｎｔ．Ｈｉｅｒａｒｃｈｉｃａｌ　ｃｌｕｓｔｅｒ　ａｎａｌｙｓｉｓ（ＨＣＡ）ｗａｓ　ｕｓｅｄ　ｔｏ　ａｎａｌｙｚｅ　ｔｈｅ　ｍａｉｎ　ｉｎｄｅｘｅｓ（ＨＣＯ－３ ，Ｃｌ－，
ＳＯ４２－，Ｋ＋，Ｎａ＋，Ｃａ２＋，ａｎｄ　Ｍｇ２＋）ｏｆ　２４ｇｒｏｕｎｄｗａｔｅｒ　ｓａｍｐｌｅｓ，ｗｈｉｃｈ　ｗｅｒｅ　ｄｉｖｉｄｅｄ　ｉｎｔｏ　ｔｈｒｅｅ　ｃｌｕｓｔｅｒｓ（Ｃ１，Ｃ２，ａｎｄ　Ｃ３）．Ｐｒｉｎｃｉ－
ｐａｌ　ｃｏｍｐｏｎｅｎｔ　ａｎａｌｙｓｉｓ（ＰＣＡ）ｗａｓ　ｕｓｅｄ　ｔｏ　ｐｅｒｆｏｒｍ　ｄｉｍｅｎｓｉｏｎ　ｒｅｄｕｃｔｉｏｎ　ｔｒｅａｔｍｅｎｔ　ｏｎ　ｅｉｇｈｔ　ｃｈｅｍｉｃａｌ　ｉｎｄｅｘｅｓ　ｏｆ　ｇｒｏｕｎｄｗａｔｅｒ，
ａｎｄ　ｔｈｒｅｅ　ｃｏｍｐｒｅｈｅｎｓｉｖｅ　ｉｎｄｅｘｅｓ（ＰＣ１，ＰＣ２ａｎｄ　ＰＣ３）ｗｅｒｅ　ｏｂｔａｉｎｅｄ．Ｔｈｅ　ｃｈｅｍｉｃａｌ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ａｎｄ　ｆｏｒｍａｔｉｏｎ　ｍｅｃｈａｎｉｓｍ　ｏｆ
ｇｒｏｕｎｄｗａｔｅｒ　ｗｅｒｅ　ａｎａｌｙｚｅｄ　ｕｓｉｎｇ　ｔｈｅ　ＰＣＡ　ａｎｄ　ｐｉｐｅｒ　ｄｉａｇｒａｍ　ｒｅｓｕｌｔｓ　ａｎｄ　ｂａｓｅｄ　ｏｎ　ｔｈｅ　ｌｏｃａｔｉｏｎ　ａｎｄ　ｒｏｃｋ　ｔｙｐｅ　ｏｆ　ｔｈｅ　ｇｒｏｕｎｄｗａｔｅｒ
ｓａｍｐｌｅｓ．Ｒｅｓｕｌｔｓ　ｓｈｏｗｅｄ　ｔｈａｔ（１）Ｃ１ｗａｔｅｒ　ｓａｍｐｌｅ　ｉｓ　ｐｏｒｅ　ｗａｔｅｒ　ｉｎ　ｌｏｏｓｅ　ｓｔｒａｔａ，ａｎｄ　ｉｔｓ　ｗａｔｅｒ　ｑｕａｌｉｔｙ　ｉｓ　ｍａｉｎｌｙ　ａｆｆｅｃｔｅｄ　ｂｙ　ｄｉｓｓｏ－
ｌｕｔｉｏｎ，ｌｅａｃｈｉｎｇ，ａｎｄ　ｉｏｎ　ｅｘｃｈａｎｇｅ　ｂａｓｅｄ　ｏｎ　ＰＣ１ａｎｄ　ＰＣ３ｉｎｄｅｘｅｓ．Ｔｈｅ　ｗａｔｅｒ　ｔｙｐｅ　ｉｓ　ＨＣＯ３·Ｃｌ－Ｃａ·Ｍｇ　ｏｒ　Ｃｌ·ＨＣＯ３－Ｃａ·Ｍｇ；
（２）Ｃ２ｗａｔｅｒ　ｓａｍｐｌｅ　ｉｓ　ｓｔｒｕｃｔｕｒａｌ　ｆｉｓｓｕｒｅ　ｗａｔｅｒ，ａｎｄ　ｔｈｅ　ｄｉｓｓｏｌｕｔｉｏｎ　ａｎｄ　ｌｅａｃｈｉｎｇ　ｅｆｆｅｃｔｓ　ｉｎ　Ｃ２ｗａｔｅｒ　ｓａｍｐｌｅ　ａｒｅ　ｌｏｗｅｒ　ｔｈａｎ　ｔｈｏｓｅ
ｉｎ　Ｃ１ｗａｔｅｒ　ｓａｍｐｌｅ　ｂａｓｅｄ　ｏｎ　ＰＣ１ａｎｄ　ＰＣ２ｉｎｄｅｘｅｓ．Ｔｈｅ　ｃｏｎｔｅｎｔｓ　ｏｆ　ｍａｉｎ　ｉｏｎｓ　ｉｎ　Ｃ２ａｒｅ　ｌｏｗｅｒ　ｔｈａｎ　ｔｈｏｓｅ　ｉｎ　Ｃ１ｂｕｔ　ｔｈｅｒｅ　ａｒｅ　ｍｏｒｅ
ｃａｒｂｏｎａｔｅ　ｒｏｃｋ　ａｎｄ　ｍａｇｎｅｓｉｕｍ　ｍｉｎｅｒａｌ　ｉｎ　ｓｕｒｒｏｕｎｄｉｎｇ　ｒｏｃｋｓ，ｗｈｉｃｈ　ｒｅｓｕｌｔｓ　ｉｎ　ｔｈｅ　ｉｎｃｒｅａｓｉｎｇ　ｏｆ　ＣＯ３－ａｎｄ　Ｍｇ２＋ ｃｏｎｔｅｎｔｓ．Ｔｈｅ
ｗａｔｅｒ　ｔｙｐｅ　ｏｆ　Ｃ２ｉｓ　ＨＣＯ３－Ｃａ·Ｍｇ　ｏｒ　ＨＣＯ３－Ｃａ；ａｎｄ（３）Ｃ３ ｗａｔｅｒ　ｓａｍｐｌｅ　ｉｓ　ｗｅａｔｈｅｒｉｎｇ　ｆｉｓｓｕｒｅ　ｗａｔｅｒ，ａｎｄ　ｔｈｅ　ｄｉｓｓｏｌｕｔｉｏｎ　ａｎｄ
ｌｅａｃｈｉｎｇ　ｅｆｆｅｃｔｓ　ｉｎ　Ｃ３ｗａｔｅｒ　ｓａｍｐｌｅ　ａｒｅ　ｅｖｅｎ　ｌｏｗｅｒ　ｂａｓｅｄ　ｏｎ　ＰＣ１ａｎｄ　ＰＣ２ｉｎｄｅｘｅｓ．Ｉｔ　ｈａｓ　ｇｏｏｄ　ｗａｔｅｒ　ｑｕａｌｉｔｙ　ｂｕｔ　ｃａｎ　ｖａｒｙ　ｓｉｇｎｉｆｉ－
ｃａｎｔｌｙ　ｉｎ　ｄｉｆｆｅｒｅｎｔ　ｓｅａｓｏｎｓ．Ｓｏｌｕｂｌｅ　ｆｌｕｏｒｉｎｅ　ｏｃｃｕｒｓ　ｉｎ　ｔｈｅ　ｓｕｒｒｏｕｎｄｉｎｇ　ｒｏｃｋｓ，ｗｈｉｃｈ　ｒｅｓｕｌｔｓ　ｉｎ　ｔｈｅ　ｉｎｃｒｅａｓｉｎｇ　ｏｆ　Ｆ－ｃｏｎｔｅｎｔ．Ｔｈｅ
ｗａｔｅｒ　ｔｙｐｅ　ｉｓ　ｍａｉｎｌｙ　ＨＣＯ３·ＳＯ４－Ｃａ·Ｍｇ．
Ｋｅｙ　ｗｏｒｄｓ：ｈｙｄｒｏｃｈｅｍｉｓｔｒｙ；ＨＣＡ；ＰＣＡ；ｄｉｓｓｏｌｕｔｉｏｎ　ｅｆｆｅｃｔ；ｉｏｎ　ｅｘｃｈａｎｇｅ
·３３１·
水文地质与工程地质
　　梭罗树水源地位于中国山东省日照市岚山区。随着经
济的发展，岚山区生活用水、工业用水量逐渐增加［１］，目前平
均日用水量可达６万 ｍ３。但是由于海水入侵，该区主要供
水水源地（绣针河水源地）下游汾水—荻水水源地已报废，现
在只有中游大朱曹水源地能够利用，但因其水量小，每天仅
供１．５万ｍ３，且已不能保障，因此供水不足已严重制约了岚
山区的经济发展。梭罗树水源地由于距岚山城区较近，可作
为城市生活用水，应急供水水源地。在此背景下，深入分析
梭罗树水源地地下水化学成分的分布与形成机制，对于揭示
本区水质特征及水质演化，促进地下水资源的合理利用、社
会经济发展和生态环境保护等具有重要的理论和实际意义，
且该地区目前尚未有过相关方面的研究。
１　研究区概况
１．１　自然地理概况
研究区属暖温带潮湿海洋季风气候区，年平均气温１２．７
℃；年平均降水量８６８．８ｍｍ，降水多集中在７－９月份；年平
均蒸发量为１　１４６．２ｍｍ。
区内地势中间低、东西及北部高，属丘陵区。东西一般海
拔高度５０．０～２８２．０ｍ，东北部大旺山为最高点，海拔高度
４１０．０ｍ；中间南北向沿龙王河由北向南，地形渐趋降低，海拔
高度由５５．０ｍ变化到２０．０ｍ左右，平均地面坡度为０．００２。
本区主要发育两条大的断裂带，分别为长脖岭构造破碎
带以及梭罗树断裂带。区内地貌特征主要受梭罗树构造断
裂的控制，呈阶梯状地形。梭罗树断裂东西两侧，依次分布
有山地和倾斜台地（阶地），构成了从山地构造剥蚀到河谷侵
蚀、堆积的完整地貌类型。区内出露地层仅见新生界第四
系，主要为冲积、冲坡积、冲洪积相砂砾层、砾砂层、含粘土砾
质砂层、含砾粘土质砂层，呈灰黄、灰褐色，结构松散，成分为
长英质及少量辉石、蛇纹石等。
１．２　水文地质条件
研究区内地下水赋存条件与分布规律严格受地层岩性、
地形地貌及构造等自然因素的控制。区内除龙王河、巨峰河
及其支流的河谷中为第四系松散堆积物外，东、西及北部丘
陵区均由岩浆岩组成，总体上是向东南开口的簸箕状地形。
岩浆岩分布区富水性差，但在地形地貌有利部位，或构造裂
隙发育带及脉岩与围岩接触面附近，地下水富集，水量较大。
第四系分布区，岩性由残坡积、坡积的亚砂土及压粘土夹碎
石透镜体组成，富水性差，而在中南部，分布于龙王和和秀针
河河谷中的第四系比较发育，厚度４～１２ｍ，最后１９．５ｍ，底
部普遍存在粗砂砾石层，直接覆盖在变质岩之上，上覆弱隔
水亚粘土层较薄，易于接受补给，因此，地下水比较丰富。区
域内主要河流龙王河、巨峰河等，均属山溪型河流。水源地
内地表水体不发育，大气降水为地下水的主要补给来源；地
下水径流方向与地表水径流方向基本一致，径流途径较短；
地下水排泄以蒸发和向地表水排泄两种方式为主，大气降水
入渗后由地形高处向低处龙王河汇集。
２　采样与研究方法
２．１　水样采集与测试
本次研究的地下水水样取自民井、机井或者泉水，共取
样２４组，采样深度为５～１７ｍ，各采样点位置见图１。采样
点采用控制性布点与功能性布点相结合的原则，在周围环境
敏感点、地下水污染源、断裂带周围以及对于确定边界条件
有控制意义的点进行布设。取样时间为平水期：２０１１年１２
月－２０１２年３月。
样品从取样至化验，操作过程和方法均按照《生活饮用
水标准检验方法》的规定进行。取样器具采用聚乙烯塑料
瓶，取样后现场用石蜡封好瓶口，贴好标签，送往实验室进
行化验。测试指标主要为 Ｃｌ－、ＳＯ４２－、ＨＣＯ－３ 、ＮＯ－３ 、Ｆ－、
Ｋ＋＋Ｎａ＋、Ｍｇ２＋、Ｃａ２＋及ＴＤＳ（总溶解固体）等。其中，Ｋ＋
＋Ｎａ＋测试方法为火焰原子吸收分光光度法（ＧＦＵ２２０２）；
ＳＯ４２－、Ｃｌ－测试方法为离子色谱法（ＩＣＳ－３０００）；Ｃａ２＋、Ｍｇ２＋
测试方法为 ＥＤＴＡ滴定法；ＨＣＯ－３ 测试方法为酸碱滴定
法；ＴＤＳ为计算值，是各离子成分相加减去 ＨＣＯ－３ 浓度的
１／２［２］。
图１　采样点位置分布
Ｆｉｇ．１　Ｄｉｓｔｒｉｂｕｔｉｏｎ　ｏｆ　ｗａｔｅｒ　ｓａｍｐｌｉｎｇ　ｐｏｉｎｔｓ
２．２　研究方法
研究地下水化学成分的方法主要包括图解法及多元统
计法。其中图解法中的Ｐｉｐｅｒ图［３－５］应用最为广泛；另外应用
较多的图解法还有ｂｏｘ　ｐｌｏｔ［６－７］、稳定场图［８］。而多元统计方
法［９－１１］由于其可以对地下水水化学特征进行定量研究而被
广泛 应 用。其 中，利 用 ＨＣＡ（ｈｉｅｒａｒｃｈｉｃａｌ　ｃｌｕｓｔｅｒ　ａｎａｌｙ－
ｓｉｓ）［１２－１３］可以对所采集的样本进行聚类分析，即可将样品进
行分组；利用ＰＣＡ（ｐｒｉｎｃｉｐａｌ　ｃｏｍｐｏｎｅｎｔ　ａｎａｌｙｓｉｓ）［６－７，１２－１３］可
以对所研究的指标进行降维处理，找出几个综合因子（ＰＣ１、
ＰＣ２、ＰＣ３……）来代表原来众多的变量，使这些综合因子尽
可能地反映原来变量的信息量，而且彼此之间互不相关，从
而达到简化的目的［１４］。
本文主要利用 ＨＣＡ与ＰＣＡ方法的结合分析各组地下
水化学特征及其形成机理。
３　数据分析结果与讨论
３．１　基于 ＨＣＡ的数据分析结果与讨论
通过对地下水水样的主要离子成分（ＨＣＯ－３ 、Ｃｌ－、
ＳＯ４２－、Ｋ＋＋Ｎａ＋、Ｃａ２＋、Ｍｇ２＋）的毫克当量浓度进行聚类分
析：利用软件ＤＰＳｖ　７．０５进行计算，步骤为：多元统计分析→
聚类分析→系统聚类。计算中将各指标值进行标准转化，选
·４３１·
第１２卷 总第７４期·南水北调与水利科技·２０１４年第５期　
水文地质与工程地质
择欧式距离系数、离差平方和法进行聚类分析。计算结果将
２４组地下水样分为Ｃ１、Ｃ２ 以及Ｃ３ 三簇，见图２。将三簇的各
主要离子的毫克当量平均值进行对比，见图３。根据各离子成
分的毫克当量浓度利用Ａｑｕａｃｈｅｍ　４．０绘制Ｐｉｐｅｒ图，见图４。
图３显示，Ｃ１ 簇阴离子含量Ｃｌ－≈ＨＣＯ－３ ＞ＳＯ４２－，阳离
子含量Ｃａ２＋＞Ｍｇ２＋＞Ｋ＋ ＋Ｎａ＋，结合图４分析，其化学类
型主要为 ＨＣＯ３·Ｃｌ－Ｃａ·Ｍｇ型或者Ｃｌ·ＨＣＯ－３Ｃａ·Ｍｇ
型；Ｃ２ 簇阴离子含量 ＨＣＯ－３ ＞Ｃｌ－ ＞ＳＯ４２－，阳离子含量
Ｍｇ２＋＞Ｃａ２＋＞Ｋ＋ ＋Ｎａ＋，结合图４分析，其化学类型主要
为 ＨＣＯ３－Ｍｇ·Ｃａ型或者 ＨＣＯ－３Ｃａ型；Ｃ３ 簇阴离子含量
ＨＣＯ－３ ＞ＳＯ４２－ ＞Ｃｌ－，阳离子含量 Ｃａ２＋ ＞Ｍｇ２＋ ＞Ｋ＋ ＋
Ｎａ＋，结合图４分析其化学类型主要为ＨＣＯ３·ＳＯ４－Ｃａ·Ｍｇ
型或者ＳＯ４·ＨＣＯ３－Ｃａ·Ｍｇ型。
图２　地下水样本的聚类分析树状图
Ｆｉｇ．２　Ｃｌｕｓｔｅｒ　ａｎａｌｙｓｉｓ　ｄｅｎｄｒｏｇｒａｍ　ｏｆ　ｔｈｅ　ｇｒｏｕｎｄｗａｔｅｒ　ｓａｍｐｌｅｓ
图３　三簇地下水样本中各指标的平均含量
Ｆｉｇ．３　Ａｖｅｒａｇｅ　ｃｏｎｔｅｎｔ　ｏｆ　ｅａｃｈ　ｉｎｄｅｘ　ｉｎ　ｔｈｅ　ｔｈｒｅｅ　ｇｒｏｕｎｄｗａｔｅｒ　ｃｌｕｓｔｅｒｓ
图４　三簇（ＨＣＡ方法得到）地下水样本的ｐｉｐｅｒ图
Ｆｉｇ．４　Ｐｉｐｅｒ　ｄｉａｇｒａｍ　ｏｆ　ｔｈｅ　ｔｈｒｅｅ　ｇｒｏｕｎｄｗａｔｅｒ
ｃｌｕｓｔｅｒｓ（ｏｂｔａｉｎｅｄ　ｆｒｏｍ　ＨＣＡ）
３．２　基于ＰＣＡ分析方法的结果与讨论
将数据进行标准化之后，利用ＤＰＳｖ７．０５对数据进行主
成分分析，操作步骤为：统计分析→多因素分析→主成分分
析。执行后得到的结果包括：相关系数矩阵Ｒ、特征方程
｜Ｒ－λＩ｜＝０的非负特征值、对应特征值的特征向量以及主
成分分析因子得分值，详细的分析计算步骤参考 ＤＰＳ数据
处理系统手册［１４］。通过分析相关系数矩阵，可以对地下水
中化学参数的相似相异性及地下水来源的一致性和差异性
进行研究［１５－１６］，在此基础上判断该数据组是否可以进行主成
分分析［１７］。而贡献率表示各主分量反映的信息量占总信息
量的百分比，累计贡献率表示所选的ｍ个主分量反映的信息
量的总和占总信息量的百分比。
３．２．１　指标间相关关系分析
查数理统计表得知，Ｒ０．０１８＝０．５０５２，Ｒ０．０５８＝０．３９６１，表
１中显示８个指标的相关系数为显著相关，３０％的数据达到
９９％的置信区间，３８％的数据达到９６％的置信区间。其中，
Ｃｌ－－ＮＯ－３ 、Ｋ＋ ＋ Ｎａ＋－Ｃｌ－、Ｃａ２＋－Ｃｌ－、ＳＯ４２－－Ｆ－、ＨＣＯ－３ －
Ｍｇ２＋、ＮＯ－３ －Ｋ＋＋Ｎａ＋、ＮＯ－３ －Ｍｇ２＋、ＮＯ－３ －Ｃａ２＋、Ｋ＋ ＋Ｎａ＋－
Ｃａ２＋呈显著正相关，通过了ａ＝０．０１水平的检验。特别是
Ｃｌ－与Ｋ＋＋Ｎａ＋相关系数高达０．９１４３，说明两者的来源可
能是相同，而在研究区围岩中没有Ｃｌ－与 Ｎａ＋的显著来源，
很可能是受海水的影响。ＨＣＯ－３ 与Ｆ－呈显著负相关，通过
了ａ＝０．０５水平的检验。
表１　指标间的相关系数矩阵
Ｔａｂ．１　Ｃｏｒｒｅｌａｔｉｏｎ　ｃｏｅｆｆｉｃｉｅｎｔ　ｍａｔｒｉｘ　ｂｅｔｗｅｅｎ　ｔｈｅ　ｉｎｄｅｘｅｓ
相关系数 Ｃｌ－ ＳＯ４２－ ＨＣＯ－３ ＮＯ－３ Ｆ－ Ｋ＋＋Ｎａ＋ Ｍｇ２＋ Ｃａ２＋
Ｃｌ－ １．０００　０
ＳＯ４２－ －０．１４７　４　 １．０００　０
ＨＣＯ－３ ０．０８０　４ －０．３０４　５　 １．０００　０
ＮＯ－３ ０．８８６　５　 ０．０１２　１　 ０．１８０　７　 １．０００　０
Ｆ－ ０．０８２　８　 ０．５５４　１ －０．４５６　９　 ０．０５７　８　 １．０００　０
Ｋ＋＋Ｎａ＋ ０．９１４　３　 ０．０３６　５　 ０．０５５　０　 ０．８８４　４　 ０．１３２　７　 １．０００　０
Ｍｇ２＋ ０．３８４　４ －０．０８６　７　 ０．８３２　３　 ０．５２６　４ －０．３３３　３　 ０．３７９　２　 １．００００
Ｃａ２＋ ０．８６９　７ －０．００５　６　 ０．０６２　５　 ０．８９８　３　 ０．１６７　０　 ０．８３５　３　 ０．２９９　０　 １．０００　０
３．２．２　ＰＣＡ分析方法的结果与讨论
主成分分析计算结果见表２，可以看出前３个主成分
（ＰＣ１、ＰＣ２ 及ＰＣ３）的累积贡献率已经达到９０．１４％，并且所
对应的特征值均大于１，故选取前３个成分代替原来的８个
·５３１·
周亚红等·山东省梭罗树水源地地下水化学特征及形成机理
水文地质与工程地质
水质指标。
表２　ＰＣＡ方法计算结果
Ｔａｂ．２　Ｃａｌｃｕｌａｔｉｏｎ　ｒｅｓｕｌｔｓ　ｏｆ　ＰＣＡ　ｍｅｔｈｏｄ
ＰＣ１ ＰＣ２ ＰＣ３
Ｃｌ－ ０．４７２　４　 ０．０９５　１ －０．２０４　１
ＳＯ４２－ －０．０４０　３　 ０．４０６　６　 ０．７００　２
ＨＣＯ－３ ０．２４９　９ －０．５５２　５　 ０．３７５　４
ＮＯ－３ ０．４８６　６　 ０．０６５　６　 ０．０１６　８
Ｆ－ ０．００５　 ０．５４８　２　 ０．２９５　６
Ｋ＋＋Ｎａ＋ ０．４６９　８　 ０．１４７　７ －０．０６４　１
Ｍｇ２＋ ０．３９９　５ －０．４０９　７　 ０．４７１　４
Ｃａ２＋ ０．４５６　４　 ０．１６０　８ －０．１１４　８
特征值 ３．９３１　３　 ２．２２６　７　 １．０５３　１
贡献率 ４９．１４０　９　 ２７．８３３　１　 １３．１６４　１
累积贡献率 ４９．１４０　９　 ７６．９７４　０　 ９０．１３８　１
　　ＰＣ１ 的贡献率为４９．１４％，可见ＰＣ１ 在地下水化学成分
的变化中起主导作用。如图５（ａ）所示，与ＰＣ１ 正相关的因
子为ＮＯ－３ 、Ｃａ２＋、Ｋ＋＋Ｎａ＋以及Ｃｌ－等，可见ＰＣ１ 主要解释
了水岩相互作用以及人类活动的影响，同时由于 Ｋ＋ ＋Ｎａ＋
以及Ｃｌ－的得分值比较接近，也从一定程度上说明ＰＣ１ 受到
近海因素的影响。ＰＣ２ 的贡献率为２７．８３％，图５（ａ）显示
ＰＣ２与Ｆ－呈较为显著的正相关，与 ＨＣＯ－３ 和 Ｍｇ２＋呈负相
关。ＰＣ３ 的贡献率为１３．１６４１％，与ＳＯ４２－及 Ｍｇ２＋呈较为显
著的正相关（图５（ｃ））。
图５（ｄ）显示，Ｃ１ 簇的水样点主要分布在ＰＣ１ 的正方向
以及ＰＣ３ 的负方向，即ＮＯ－３ 、Ｃａ２＋、Ｋ＋＋Ｎａ＋以及Ｃｌ－含量
较为显著增长的同时，ＳＯ４２－ 及 Ｍｇ２＋ 有较为显著的降低。
Ｃ１ 簇水样中大部分分布于东部近海岸及巨峰河、龙王河河
谷下游地带，其中Ｓ３１７、Ｓ３１８分布于巨峰河河谷漫滩阶地
中；Ｓ３２５、Ｓ３２７位于龙王河上游；Ｓ３３４、Ｓ３３７位于龙王河下
游；Ｓ３１９、Ｓ３２０位于大旺山山前、山丘的倾斜台地上，因此该
簇水样属于松散岩类孔隙水，地下水埋藏浅，水力坡度小，径
流迟缓，蒸发作用、溶解作用及淋滤作用显著。由此可知Ｃ１
成分主要来自可溶岩的溶解作用，同时强烈的蒸发作用导致
各离子含量较高，该簇各点ＴＤＳ均值达５３０．３５ｍｇ／Ｌ。由于
分布于近海岸，含水层中有一定的海相沉积物，阴阳离子交
换作用使Ｃｌ－从沉积物中析出，而ＳＯ４２－被吸附，地下水中
ＳＯ４２－的显著减少，而含水层中的Ｋ＋＋Ｎａ＋与水中 Ｍｇ２＋发
生离子交换，造成水中 Ｍｇ２＋含量降低。此外，孔隙水是区
内人们生活用水的主要来源，受人类活动的影响较大，污染
相对严重，所以ＮＯ－３ 含量较高。
Ｃ２ 簇水样点主要分布与ＰＣ２ 的负方向以及ＰＣ１ 的负方
向（图５（ｂ）），即ＮＯ－３ 、Ｃａ２＋、Ｋ＋＋Ｎａ＋、Ｃｌ－、Ｆ－明显降低的
同时，Ｍｇ２＋以及 ＨＣＯ－３ 增加。Ｃ２ 簇水样主要位于构造裂隙
带中，属于构造裂隙水，溶解作用相对Ｃ１ 簇水样较弱，所以
大部分离子含量都有所降低，ＴＤＳ均值为３６２．３９ｍｇ／Ｌ。其
中Ｓ０４７、Ｓ３２９、Ｓ３２４位于长脖岭构造破碎带，带内岩石以构
造角砾岩为主，呈灰白色、灰红色，角砾状构造；Ｓ１０１、Ｓ３２３、
Ｓ３４１、Ｓ３４３、Ｓ３４９位于梭罗树断裂破碎带，该断裂切割了晚
元古代梭罗树－水车沟蛇纹岩体，两断裂破碎带中伴随有碳
酸盐化及绿泥石化现象，围岩中碳酸盐及含 Ｍｇ２＋矿物含量
较多，所以地下水中 Ｍｇ２＋和 ＨＣＯ－３ 有所增加。
图５　ＰＣＡ分析结果显示图
（（ａ）、（ｂ）：化学指标分布图；（ｂ）、（ｄ）：样本分布图）
Ｆｉｇ．５　Ｒｅｓｕｌｔｓ　ｏｆ　ＰＣＡ（ａ、ｃ：ｄｉｓｔｒｉｂｕｔｉｏｎ　ｏｆ　ｃｈｅｍｉｃａｌ　ｉｎｄｅｘｅｓ；
ｂ、ｄ：ｄｉｓｔｒｉｂｕｔｉｏｎ　ｏｆ　ｇｒｏｕｎｄｗａｔｅｒ　ｓａｍｐｌｅｓ）
·６３１·
第１２卷 总第７４期·南水北调与水利科技·２０１４年第５期　
水文地质与工程地质
Ｃ３ 簇水样点主要分布于ＰＣ１ 的负方向，ＰＣ２ 的正方向
（图 ５（ｂ）），即 ＮＯ－３ 、Ｃａ２＋、Ｋ＋ ＋Ｎａ＋、Ｃｌ－、Ｍｇ２＋ 以及
ＨＣＯ－３ 均减少的同时，Ｆ－显著增加。Ｃ３ 簇水样中的大部分
点都位于区内片麻状二长花岗岩的风化裂隙中，岩石主要包
括细粒石英二长岩、混合花岗岩、闪长岩、闪长玢岩和各种脉
岩，岩石结构致密坚硬，深部裂隙一般不发育，只在浅部发育
风化裂隙，发育深度一般为５．０～１５．０ｍ，富水性较差，季节
变化较大。Ｃ３ 簇水样主要为基岩裂隙水，与Ｃ２ 簇水样相
似，溶解作用弱，地下水中各组分的含量偏低，ＴＤＳ均值为
２４５ｍｇ／Ｌ，但是由于其周围岩石中含有少量的可溶性氟通过
溶滤作用进入地下水，引起水中Ｆ－含量增加。
４　结语
（１）梭罗树水源地地下水主要分为三种类型：松散岩类
孔隙水（ＨＣＯ３·Ｃｌ－Ｃａ·Ｍｇ型或Ｃｌ·ＨＣＯ３－Ｃａ·Ｍｇ型）、
构造裂隙水（ＨＣＯ３－Ｍｇ·Ｃａ型或者 ＨＣＯ３－Ｃａ型）、风化裂隙
水（ＨＣＯ３·ＳＯ４－Ｃａ·Ｍｇ型或者ＳＯ４·ＨＣＯ３－Ｃａ·Ｍｇ）。
（２）松散岩类孔隙水化学成份主要受可溶岩的溶解和淋
滤作用以及蒸发作用影响，部分水体受海水影响而Ｃｌ－含量
较高；构造裂隙水岩石的溶解与淋滤影响；风化裂隙水水质
相对较好，但受季节影响较大，并且含有少量的Ｆ－。
参考文献（Ｒｅｆｅｒｅｎｃｅｓ）：
［１］　孟艳丽．日照市岚山区土地利用现状分析［Ｊ］．山东国土资源，
２００７，２３（３）：１４－１６．（ＭＥＮＧ　Ｙａｎ－ｌｉ．Ａｎａｌｙｓｉｓ　ｏｎ　ｐｒｅｓｅｎｔ　ｃｏｎｄｉ－
ｔｉｏｎ　ｏｆ　ｌａｎｄ　ｕｔｉｌｉｚａｔｉｏｎ　ｉｎ　Ｌａｎｓｈａｎ　ｄｉｓｔｒｉｃｔ　ｏｆ　Ｒｉｚｈａｏ　Ｃｉｔｙ）［Ｊ］．
Ｓｈａｎｄｏｎｇ　Ｐｒｏｖｉｎｃｅ　Ｌａｎｄ　ａｎｄ　Ｒｅｓｏｕｒｃｅｓ，２００７，２３（３）：１４－１６．（ｉｎ
ｃｈｉｎｅｓｅ））
［２］　钱会，水文地球化学［Ｍ］．北京：地质出版社，２００６．（ＱｉａｎＨｕｉ，
ｈｙｄｒｏ－ｇｅｏｃｈｅｍｉｓｔｒｙ［Ｍ］．Ｂｅｉｊｉｎｇ：Ｇｅｏｌｏｇｉｃａｌ　ｐｕｂｌｉｓｈｉｎｇ　Ｈｏｕｓｅ，
２００６．（ｉｎ　ｃｈｉｎｅｓｅ）
［３］　Ｊｉａｎｈｕａ　Ｈｅ，Ｊｉｎｚｈｕ　Ｍａ，Ｐｅｎｇ　Ｚｈａｎｇ，ｅｔ　ａｌ．Ｇｒｏｕｎｄｗａｔｅｒ　ｒｅ－
ｃｈａｒｇｅ　ｅｎｖｉｒｏｎｍｅｎｔｓ　ａｎｄ　ｈｙｄｒｏ　ｇｅｏｃｈｅｍｉｃａｌ　ｅｖｏｌｕｔｉｏｎ　ｉｎ　ｔｈｅ　ｊｉ－
ｕｑｕａｎ　ｂａｓｉｎ，Ｎｏｒｔｈｗｅｓｔ　Ｃｈｉｎａ［Ｊ］．Ａｐｐｌｉｅｄ　Ｇｅｏｃｈｅｍｉｓｔｒｙ，２０１２
（２７）：８６６－８７８．
［４］　Ｑｉｎｇｈａｉ　Ｇｕｏ，Ｙａｎｘｉｎ　Ｗａｎｇ，Ｗｅｉ　Ｌｉｕ．Ｍａｊｏｒ　ｈｙｄｒｏｇｅｏｃｈｅｍｉｃａｌ
ｐｒｏｃｅｓｓｅｓ　ｉｎ　ｔｈｅ　ｔｗｏ　ｒｅｓｅｒｖｏｉｒｓ　ｏｆ　ｔｈｅ　Ｙａｎｇｂａｊｉｎｇ　ｇｅｏｔｈｅｒｍａｌ
ｆｉｅｌｄ，Ｔｉｂｅｔ，Ｃｈｉｎａ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｖｏｌｃａｎｏｌｏｇｙ　ａｎｄ　Ｇｅｏｔｈｅｒｍａｌ
Ｒｅｓｅａｒｃｈ，２００７（１６６）：２５５－２６８．
［５］　Ｃｌａｕｓ　Ｋｏｈｆａｈｌ，Ｃｈｒｉｓｔｏｐｈ　Ｓｐｒｅｎｇｅｒ，Ｊｏｓｅ　Ｂｅｎａｖｅｎｔｅ　Ｈｅｒｒｅｒａ．
Ｒｅｃｈａｒｇｅ　ｓｏｕｒｃｅｓ　ａｎｄ　ｈｙｄｒｏ　ｇｅｏｃｈｅｍｉｃａｌ　ｅｖｏｌｕｔｉｏｎ　ｏｆ　ｇｒｏｕｎｄｗ－
ａｔｅｒ　ｉｎ　ｓｅｍｉａｒｉｄ　ａｎｄ　ｋａｒｓｔｉｃ　ｅｎｖｉｒｏｎｍｅｎｔｓ：Ａ　ｆｉｅｌｄ　ｓｔｕｄｙ　ｉｎ　ｔｈｅ
ｇｒａｎａｄａ　ｂａｓｉｎ（Ｓｏｕｔｈｅｒｎ　Ｓｐａｉｎ）［Ｊ］．Ａｐｐｌｉｅｄ　Ｇｅｏｃｈｅｍｉｓｔｒｙ，
２００８（２３）：８４６－８６２．
［６］　Ｚａｒｇｈａｍ　Ｍｏｈａｍｍａｄｉ．Ａｓｓｅｓｓｉｎｇ　ｈｙｄｒｏｃｈｅｍｉｃａｌ　ｅｖｏｌｕｔｉｏｎ　ｏｆ
ｇｒｏｕｎｄｗａｔｅｒ　ｉｎ　ｌｉｍｅｓｔｏｎｅ　ｔｅｒｒａｉｎ　ｖｉａ　ｐｒｉｎｃｉｐａｌ　ｃｏｍｐｏｎｅｎｔ　ａｎａｌｙ－
ｓｉｓ［Ｊ］．Ｅｎｖｉｒｏｎ　Ｅａｒｔｈ，２００９（５９）：４２９－４３９．
［７］　Ｂｅａｔｒｉｚ　Ｈｅｌｅｎａ，Ｒａｆａｅｌ　Ｐａｒｄｏ，Ｍａｒｉｓｏｌ　Ｖｅｇａ，ｅｔ　ａｌ．Ｔｅｍｐｏｒａｌ　ｅ－
ｖｏｌｕｔｉｏｎ　ｏｆ　ｇｒｏｕｎｄｗａｔｅｒ　ｃｏｍｐｏｓｉｔｉｏｎ　ｉｎ　ａｎ　ａｌｕｖｉａｌ　ａｑｕｉｆｅｒ
（Ｐｉｓｕｅｒｇａ　Ｒｉｖｅｒ，Ｓｐａｉｎ）ｂｙ　ｐｒｉｎｃｉｐａｌ　ｃｏｍｐｏｎｅｎｔ　ａｎａｌｙｓｉｓ［Ｊ］．
Ｗａｔｅｒ　Ｒｅｓｅａｒｃｈ，２０００，３４（３）：８０７－８１６．
［８］　Ｈｕａｍｉｎｇ　Ｇｕｏ，Ｙａｎｘｉｎ　Ｗａｎｇ．Ｈｙｄｒｏｇｅｏｃｈｅｍｉｃａｌ　ｐｒｏｃｅｓｓｅｓ　ｉｎ
ｓｈａｌｏｗ　ｑｕａｔｅｒｎａｒｙ　ａｑｕｉｆｅｒｓ　ｆｒｏｍ　ｔｈｅ　ｎｏｒｔｈｅｒｎ　ｐａｒｔ　ｏｆ　ｔｈｅ　Ｄａ－
ｔｏｎｇ　Ｂａｓｉｎ，Ｃｈｉｎａ［Ｊ］．Ａｐｐｌｉｅｄ　Ｇｅｏｃｈｅｍｉｓｔｒｙ，２００４，１９：１９－２７．
［９］　Ｍ　Ｖ　Ｐｒａｓａｎｎａ，Ｓ　Ｃｈｉｄａｍｂａｒａｍ，Ｋ　Ｓｒｉｎｉｖａｓａｍｏｏｒｔｈｙ．Ｓｔａｔｉｓｔｉｃａｌ
ａｎａｌｙｓｉｓ　ｏｆ　ｔｈｅ　ｈｙｄｒｏｇｅｏｃｈｅｍｉｃａｌ　ｅｖｏｌｕｔｉｏｎ　ｏｆ　ｇｒｏｕｎｄｗａｔｅｒ　ｉｎ
ｈａｒｄ　ａｎｄ　ｓｅｄｉｍｅｎｔａｒｙ　ａｑｕｉｆｅｒｓ　ｓｙｓｔｅｍ　ｏｆ　Ｇａｄｉｌａｍ　ｒｉｖｅｒ　ｂａｓｉｎ，
Ｓｏｕｔｈ　Ｉｎｄｉａ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｋｉｎｇ　Ｓａｕｄ　Ｕｎｉｖｅｒｓｉｔｙ（Ｓｃｉｅｎｃｅ），
２０１０，２２：１３３－１４５．
［１０］　Ｍａｔｔｈｉａｓ　Ｒａｉｂｅｒ，Ｐａｕｌ　Ａ　Ｗｈｉｔｅ，Ｃｈｒｉｓｔｏｐｈｅｒ　Ｊ　Ｄａｕｇｈｎｅｙ，ｅｔ
ａｌ．Ｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌ　ｇｅｏｌｏｇｉｃａｌ　ｍｏｄｅｌｉｎｇ　ａｎｄ　ｍｕｌｔｉｖａｒｉａｔｅ
ｓｔａｔｉｓｔｉｃａｌ　ａｎａｌｙｓｉｓ　ｏｆ　ｗａｔｅｒ　ｃｈｅｍｉｓｔｒｙ　ｄａｔａ　ｔｏ　ａｎａｌｙｓｅ　ａｎｄ　ｖｉｓ－
ｕａｌｉｓｅ　ａｑｕｉｆｅｒ　ｓｔｒｕｃｔｕｒｅ　ａｎｄ　ｇｒｏｕｎｄｗａｔｅｒ　ｃｏｍｐｏｓｉｔｉｏｎ　ｉｎ　ｔｈｅ
Ｗａｉｒａｕ　Ｐｌａｉｎ，Ｍａｒｌｂｏｒｏｕｇｈ　Ｄｉｓｔｒｉｃｔ，Ｎｅｗ　Ｚｅａｌａｎｄ［Ｊ］．Ｊｏｕｒｎａｌ
ｏｆ　Ｈｙｄｒｏｌｏｇｙ，２０１２，４３６－４３７：１３－３４．
［１１］　Ｒａｊｅｓｈ　Ｒｅｇｈｕｎａｔｈ，Ｔ　Ｒ　Ｓｒｅｅｄｈａｒａ　Ｍｕｒｔｈｙ，Ｂ　Ｒ　Ｒａｇｈａｖａｎ．
Ｔｈｅ　ｕｔｉｌｉｔｙ　ｏｆ　ｍｕｌｔｉｖａｒｉａｔｅ　ｓｔａｔｉｓｔｉｃａｌ　ｔｅｃｈｎｉｑｕｅｓ　ｉｎ　ｈｙｄｒｏｇｅｏ－
ｃｈｅｍｉｃａｌ　ｓｔｕｄｉｅｓ：ａｎ　ｅｘａｍｐｌｅ　ｆｒｏｍ　Ｋａｒｎａｔａｋａ，Ｉｎｄｉａ［Ｊ］．Ｗａｔｅｒ
Ｒｅｓｅａｒｃｈ，２００２，３６（１０）：２４３７－２４４２．
［１２］　Ｅｌｉｊａｈ　Ｗａｎｄａ，Ｍａｕｒｉｃｅ　Ｍｏｎｊｅｒｅｚｉ，Ｊｏｎａｓ　Ｆ　Ｍｗａｔｓｅｔｅｚａ，ｅｔ　ａｌ．
Ｈｙｄｒｏ－ｇｅｏｃｈｅｍｉｃａｌ　ａｐｐｒａｉｓａｌ　ｏｆ　ｇｒｏｕｎｄｗａｔｅｒ　ｑｕａｌｉｔｙ　ｆｒｏｍ
ｗｅａｔｈｅｒｅｄ　ｂａｓｅｍｅｎｔ　ａｑｕｉｆｅｒｓ　ｉｎ　Ｎｏｒｔｈｅｒｎ　Ｍａｌａｗｉ［Ｊ］．Ｐｈｙｓｉｃｓ
ａｎｄ　Ｃｈｅｍｉｓｔｒｙ　ｏｆ　ｔｈｅ　Ｅａｒｔｈ，２０１１（３６）：１１９７－１２０７．
［１３］　Ｓａｎｄｏｗ　Ｍａｒｋ　Ｙｉｄａｎａ，Ｄｕｋｅ　Ｏｐｈｏｒｉ，Ｂｒｕｃｅ　Ｂａｎｏｅｎｇ－Ｙａｋｕｂｏ．
Ｈｙｄｒｏｃｈｅｍｉｃａｌ　ｅｖａｌｕａｔｉｏｎ　ｏｆ　ｔｈｅ　ｖｏｌｔａｉａｎ　ｓｙｓｔｅｍ－Ｔｈｅ　ａｆｒａｍ
ｐｌａｉｎ　ａｒｅａ，Ｇｈａｎａ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｅｎｖｉｒｏｎｍｅｎｔａｌ　Ｍａｎａｇｅｍｅｎｔ，
２００８，８８（４）：６９７－７０７．
［１４］　唐启义，冯明光．ＤＰＳ数据处理系统—实验设计、统计分析及
模型优化［Ｍ］．科学出版社，６４４－６５０．（ＴＡＮＧ　Ｑｉ－ｙｉ，ＦＥＮＧ
Ｍｉｎｇ－ｇｕａｎｇ．ＤＰＳ　ｄａｔａ　ｐｒｏｃｅｓｓｉｎｇ　ｓｙｓｔｅｍ－ｅｘｐｅｒｉｍｅｎｔａｌ　ｄｅ－
ｓｉｇｎ，ｓｔａｔｉｓｔｉｃａｌ　ａｎａｌｙｓｉｓ　ａｎｄ　ｍｏｄｅｌｉｎｇ［Ｍ］．Ｓｃｉｅｎｃ　Ｐｒｅｓｓ，６４４－
６５０．（ｉｎ　ｃｈｉｎｅｓｅ））
［１５］　章光新，邓伟，何岩，ｅｔ　ａｌ．中国东北松嫩平原地下水水化学特
征与演变规律［Ｊ］．水科学进展，２００６，１７（１）：３４２．（ＺＨＡＮＧ
Ｇｕａｎｇｘｉｎ，ＤＥＮＧ　Ｗｅｉ，ＨＥ　Ｙａｎ，ｅｔ　ａｌ．Ｈｙｄｒｏｃｈｅｍｉｃａｌ　Ｃｈａｒａｃ－
ｔｅｒｉｓｔｉｃｓ　ａｎｄ　Ｅｖｏｌｕｔｉｏｎ　Ｌａｗｓ　ｏｆ　Ｇｒｏｕｎｄｗａｔｅｒ　ｉｎ　Ｓｏｎｇｎｅｎ
Ｐｌａｉｎ，Ｎｏｒｔｈｅａｓｔ　Ｃｈｉｎａ［Ｊ］．Ａｄｖａｎｃｅｓ　ｉｎ　Ｗａｔｅｒ　Ｓｃｉｅｎｃｅ　２００６，
１７（１）：３４２．（ｉｎ　Ｃｈｉｎｅｓｅ））
［１６］　唐然，温忠辉，束龙仓，等．绥化市地下水水化学特征及水质现
状评价［Ｊ］．水资源保护，２０１３，２９（４）：１９－２５．（ＴＡＮＧ　Ｒａｎ，
ＷＥＮ　Ｚｈｏｎｇ－ｈｕｉ，ＳＨＵ　Ｌｏｎｇ－ｃａｎｇ，ｅｔ　ａｌ．Ｈｙｄｒｏ－ｃｈｅｍｉｃａｌ
Ｃｈａｒｃａｔｅｒｓｔｉｃｓ　ａｎｄ　Ｗａｔｅｒ　Ｑｕａｌｉｔｙ　Ａｓｓｅｓｓｍｅｎｔ　ｏｆ　Ｇｒｏｕｎｄｗａ－
ｔｅｒ　ｉｎ　Ｓｕｉｈｕａ　Ｃｉｔｙ［Ｊ］．Ｗａｔｅｒ　Ｒｅｓｏｕｒｃｅｓ　Ｐｒｏｔｅｃｔｉｏｎ，２０１３，２９
（４）：１９－２５．（ｉｎ　Ｃｈｉｎｅｓｅ））
·７３１·
周亚红等·山东省梭罗树水源地地下水化学特征及形成机理