鲁西北平原地下水高氟与高碘成因分析

吴光伟，李浩林，王庆兵，顾莎，郭梓烽，刘中业

（1. 山东省国土空间生态修复中心，山东 济南 250014； 2. 中南民族大学资源与环境学院，湖北 武汉 430074； 3. 山东省地质矿产勘查开发局八〇一水文地质工程地质大队，山东 济南 250014； 4. 山东昊润自动化技术有限公司，山东 滨州 256209； 5. 山东省地质调查院，山东 济南 250014）

根据《全国重点地方病防治规划》（2002—2010 年），4194 万人居住于中国饮用水型氟中毒（氟含量大于1mg/L，见GBT 14848—2017《地下水质量标准》）地区。另根据《中国水源性高碘地区划分标准》（GB/T 19380—2016）饮用水碘中位数限值为100μg/L，以此值为碘含量标准值，则目前中国有9个省份3000 多万人饮用水存在碘过量问题，山东省同时被划分为全国饮用水型氟中毒地区和饮用水碘过量地区［1-2］。鲁西北平原是山东省水源性高碘典型地区之一［3］，地处黄河冲积平原，地下水中氟与碘含量普遍偏高。地下水作为当地居民的主要生活饮用水来源，人们长期饮用高氟高碘地下水会导致氟中毒与碘中毒，氟斑牙、甲状腺肿大已成为鲁西北平原的地方病［4-6］。避免鲁西北平原居民受到高氟高碘地下水危害，是亟待解决的问题。

2010 年以来，许多学者对鲁西北及其他地区开展氟中毒流行病学调查［7-9］，查清了浅层高氟地下水化学特征、分布规律及成因分析［9-11］，开展健康风险评估［12］，但对高碘地下水以及其与高氟地下水之间关系的研究较少。已有研究表明：①弱碱性环境、HCO3−含量高的区域、较浅的埋深、较平缓的地势和强烈的蒸发浓缩作用利于高氟地下水的形成。当地下水以垂直交替为主时，水解和阳离子吸附交替作用明显，此时方解石、萤石和石膏的溶解与沉淀平衡以及阴阳离子交替吸附作用主要影响浅层高氟地下水的形成和分布［9，13-14］；②土壤类型和人类过度开采也影响高氟地下水的形成和分布［15］。例如黏土中氟离子含量高于粉砂中氟离子含量［16］；③沉积环境和地下径流条件是影响中国高碘地下水分布的主要因素［17］；④一个地区的环境地质状况、是否有咸水入侵、大气降水状况以及是否有频繁人类活动如地下水的开采是影响该地区高碘地下水分布的重要因素［18］。总之，浅层地下水的蒸发作用和深部富含有机质的、偏还原的地下水环境中的微生物作用是形成高碘地下水的重要机制［19-21］。地下水作为鲁西北平原的主要水源，查明高氟、高碘地下水的成因和联系十分必要，可为当地的饮用水安全及解决地方病的研究提供借鉴。

多元统计技术已被广泛应用于分析水文地球化学数据［22-24］，该方法不仅有助于理解各因素间的因果关系，还可以提供有关水化学过程的信息［25］。因此，本文采用层次聚类分析（HCA）方法，划分出高氟高碘地区，再结合常规水化学指标和沉积物组分特征，探讨控制鲁西北平原地下水中氟与碘富集的主要地球化学过程。本研究的主要目的是：①揭示鲁西北平原地下水中氟与碘的空间分布特征；②推测高氟地下水的形成机制；③查清碘在地下水系统中富集的关键水文地球化学过程；④探讨地下水中高氟与高碘的关系。

1 研究区概况

鲁西北平原属于华北平原的一部分，在国家水文地质分区中位于东部大平原半湿润气候季风带水文地质分区的黄淮海平原亚区，鲁西北平原松散盐类水文地质区。鲁西北平原来源于华北地台渤海凹陷，由于晚第三纪的新构造运动，渤海凹陷加剧，海水入侵到凹陷至海平面以下的区域［16］。以黄河为主的多泥沙河流不断冲填，使地面不断升高，海水不断东退，逐渐形成了广阔的黄泛大平原和深厚的第四纪沉积层［19］（图1），厚度在200～900 余米。鲁西北平原岩浆活动剧烈，岩浆岩广泛发育，岩浆活动后期形成萤石。岩浆活动形成的岩浆岩主要有两类：侵入岩和和火山岩。侵入岩分布最广泛，岩石类型复杂，中生代印支期和燕山期侵入岩以花岗岩为主，由更长石、钠长石、钾长石石英、黑云母或白云母组成。鲁西北地区自西南向东北地势逐渐降低。由于黄河多次改道，地面略有起伏，地表以粉砂、黏质砂土为主。

图1 研究区域与采样点分布图（图中黑色箭头指向为地下水流向）Fig. 1 The study area and distribution of sampling points.

鲁西北平原位于黄河下游以北，属半干旱大陆性气候，海拔3～70m，地形平坦，地下水资源丰富。年平均降水量为600mm（70%集中在夏季），年平均蒸发量为1800～2000mm。区域内流经的较大河流有马颊河、京杭大运河、徒骇河等。潜水水位的变化过程与降水的变化过程一致，时间上略有延迟，在一个水文年内呈现下降—上升—下降的规律。鲁西北平原潜水水平运动缓慢，蒸发为其主要排泄方式。黄河等地上河流经鲁西北平原，经常性补给地下水。地下水水平运动条件差，在强烈蒸发作用下，原有的水盐均衡被改变，经常发生土壤盐渍化。

鲁西北平原地区的地下水类型主要为第四系孔隙水，鲁西北平原从山前到平原、从平原到滨海，水文地质条件具有明显的变化规律，松散岩类含水岩组富水程度从东北向西南从极弱到强，滨海区域存在咸水层（矿化度大于2g/L）。浅层地下水山前倾斜平原区地下水径流较快为全淡水区，黄泛平原区地下水流速滞缓变为淡水—咸水—淡水的三层水质结构区，滨海平原区成为咸水—淡水的双层水质结构以至全咸水区。

2 实验部分

2.1 研究方法

2.1.1 样品采集

项目组分别于2019 年5～7 月和2020 年4～6月在163 个采样点共采集地下水样品326 件，包括浅层地下水（地下水埋深＜55m）样品132 件、中层地下水（地下水埋深范围为55～225m）样品64 件以及深层地下水（地下水埋深＞225m）样品130 件。采用便携式水质多参数分析仪现场测定地下水样品的温度（T）和pH 值；采用便携式水质分析仪现场测定地下水中氧化还原敏感指标（HS−、N）。所有仪器在使用前均进行校正。样品均现场通过0.45μm微孔滤膜，过滤后水样分成三份：一份装入100mL 干净的棕色塑料瓶中，用于地下水中碘元素的分析；两份装入100mL 干净的聚乙烯瓶中，其中一份加入优级纯硝酸进行酸化至pH＜2，用于主要阳离子和微量金属元素的分析；另一份不酸化，用于主要阴离子的分析。所有地下水样品在检测之前均在4℃条件下冷藏保存。另于2020 年12 月在研究区内钻取两个钻孔（G09、ZK02）采集沉积物样品。根据沉积物岩性变化每隔1～3m 截取一段15cm 长的沉积物样品，立即用保鲜膜包裹，再用锡箔纸包裹严实，置于30cm 高的塑料罐中，拧紧盖子，−18℃冷冻保存。

2.1.2 样品分析

地下水样品测定指标均按照《地下水质量标准》（GB/T 14848—2017）规定的具体要求执行。总硬度、钠、铁、锰、铜、锌、铝、铅及其他阳离子采用电感耦合等离子体发射光谱法（ICP-OES）进行测定；溶解性固体总量（TDS）采用105℃干燥重量法测定；挥发性酚采用分光光度法测定；微量元素如砷、镉等和碘化物采用电感耦合等离子体质谱法（ICP-MS）进行测定；氟化物、氨氮及部分阴离子采用离子色谱法（IC）进行测定。采样后24 h 内采用酸碱中和滴定法测定总碱度。采用波长射散型X 射线荧光光谱仪（XRF）及X 射线衍射仪（XRD）分别测定沉积物化学元素组成及其矿物组成。样品分析各指标的测定方法及其检出限见表1。

表1 样品分析指标的测定方法、仪器、检出限和精密度Table 1 Measurement method, instrument, detection limit and precision of analysis indicators for samples.

2.1.3 多元统计分析的数据准备

本文用平方欧几里德距离将参数划分为初始簇，利用组间连接方法连接产生的初始簇［21，23-24］。然后，利用SPSS 软件（V26.0）进行层次聚类分析（HCA）。

对于浓度低于检测限的组分，在进行多元统计分析之前，用检测限值表示浓度。为了达到正态分布和均匀性的目的，将数据标准化为相应的Z分数，该分数是通过从每个数据中（Xi）减去正态分布的平均值，再除以分布的标准差（s）得到的，即Zi=（Xi−平均值）/s［26］。进行多元统计分析时，为了避免欧氏距离的计算受方差最大参数的影响，必须要进行数据标准化处理［22-23］。

对于采集的数据进行合理性检验，选取8 大离子（K+、Na+、Ca2+、Mg2+、、、Cl−、）进行阴阳离子电荷守恒平衡的检验，先计算出阴阳离子的毫克当量浓度（meq/L）分别为Σmc和Σma，再利用公式E=（Σmc−Σma）/（Σmc+Σma）×100%计算出相对误差E。数据为实测值时，E一般应小于±5%。计算结果表明相对误差在±5%以内的样品占80%以上，相对误差在±10%以内的样品占90%，由于部分样品在人类聚集区域采集，考虑误差较大的样品出现的原因是人类活动造成的地下水污染。

3 结果

3.1 地下水及沉积物组分特征

从表2 地下水样品测试结果可看出，研究区不同深度的地下水样品pH 集中在6.25～10.42 之间，属弱碱性环境，TDS 含量差异大，部分区域铁、锰、铜、铝含量异常高，汞、砷、铅也有不同程度的超标。沉积物分析结果显示，研究区岩性浅层主要为粉质黏土和粉土，夹杂少量粉砂和细砂，中层转向粉质黏土、粉土和黏土为主，夹杂少量粉砂和细砂，深层以黏土和粉质黏土为主，夹杂少量粉砂和细砂。沉积物pH 值为7.23～9.44，均值为8.05；浅层、中层、深层的氟含量均值分别为0.00672mg/L、0.01033mg/L、0.01099mg/L；与中层和深层地下水相比，浅层地下水具有更高浓度的Na+、Cl−、（表3）。研究区矿物主要由石英组成，占比达50%以上，其次为黏土矿物和长石，占比在10%～20%，有少量方解石、白云石等，占比小于10%。研究区地下水平均流速分别为：浅层0.0137cm/d，中层0.0133cm/d，深层0.0127cm/d。

表2 地下水水化学组分特征Table 2 Hydrochemistry indicators in groundwater at different depths.

表3 钻孔沉积物组分特征Table 3 The content of each component in sediments at different depths in the two holes.

3.2 高氟高碘地下水划分及其特征

3.2.1 浅层地下水

本文对132 个浅层地下水样品进行了18 个水化学参数（pH、溶解性总固体、硫酸根、氯、铁、挥发性酚、碘、氟、砷、铅、钙、镁、钾、重碳酸根）的层次聚类分析（HCA）（图2）。按照如下步骤对地下水样品进行分类：①对相似性较大的样品进行分组；②使用2.2 连接规则连接样本组；③重复步骤①和②，直到对所有观察结果进行分类。通过在树状图中向上或向下移动Y轴指标线的位置，可以确定较少或更多的组。在本次研究中，组簇数目并不是由单一同型种线直接确定的。例如，考虑到样品矿化度（TDS）含量差异大，将高TDS（TDS≥20g/L）与低TDS（TDS＜20g/L）的分界线也定义为同型种线（Y=7.5，红色虚线，图2），故样本可大致分为高TDS 组与低TDS 组。

图2 浅层地下水样本的树状图。聚类分析时使用平均连接（组间），红色与绿色虚线为同型种线，对应Y 分别为7.5 和1.5。浅层地下水可分为高TDS 区和低TDS 区，具有高氟高碘特征的簇数目最多，其中S2、S3、S4、S6、S8、S13、S17、S18 为低TDS 组中的高氟高碘簇Fig. 2 The tree map of shallow groundwater samples (the red dashed line Y=7.5 and green dashed line Y=1.5 are phenon lines).Shallow groundwater can be divided into high TDS area and low TDS area, and the number of clusters with high fluoride and high iodine characteristics is the largest, among which S2, S3, S4, S6, S8, S13, S17, S18 are high fluorine and high iodine clusters in the low TDS group.

位于TDS 分界线左边的样本被认为是同一个组，即高TDS 组。因此，高TDS 组中的样本可分为6 个簇（S19～S24）。其中，S19 具有异常高的碘含量，S20、S23、S24 具有异常高的氟含量。同时，S20、S23、S24 还分别具有较高的硝酸盐、硫化物及钾含量；S21 具有较高的砷含量，S22 具有较高的铁含量。以上浓度异常指标可能与生活污水和工业污染的增加有关［27］。

在低TDS 组中，为划分出更多的组，便于分析具有最大相似性样本的特征，同型种线被定义为Y=1.5（绿色虚线，图2），共产生18 个簇（S1～S18）。S2、S3、S4、S6、S8、S13、S17、S18 为低TDS 组中的高氟、碘区；S15、S9 为高氟区；S5、S7、S10、S14 为高碘区。其中，S2 与S3 具有相同的高氟、碘特征，但S2由于锰含量较高（＞2mg/L）而被区分为不同的簇。同理，S12 因具有较高的硝酸盐含量而被单独划分。此外，为清晰刻画最大聚类的特征，本文绘制了S1 聚类的树状图（图3），根据氟含量标准值（1.0mg/L，GB/T 14848—2017《地下水质量标准》），由同型种线（Y=5.5，红色虚线，图3）共构成12 个簇，分别命名为S1-1～S1-12。其中S1-2、S1-3、S1-8、S1-10 和S1-12 为高氟、高碘区；S1-9 仅为高碘区。

图3 S1 聚类树状图（红色虚线为同型种线，对应Y 为5.5），浅层地下水最大聚类S1 聚类结果显示S1-2、S1-3、S1-8、S1-10和S1-12 为高氟高碘簇Fig. 3 The tree diagram of S1 clusters (the red dashed line Y=5.5 is phenon line), the results of the maximum cluster S1 cluster of shallow groundwater showed that S1-2, S1-3, S1-8, S1-10 and S1-12 are high fluorine and high iodine clusters.

3.2.2 中层地下水

中层地下水的HCA 为划分出最多簇（19 个簇），同型种线定义为Y=1.5。据此，中层地下水样本主要分为四大聚类：高氟样本（M7、M13），高碘样本（M12、M15、M16、M17），高氟和碘样本（M1、M2、M3、M4、M5、M9、M10、M11、M13、M14、M18、M19）以及正常样本（M6、M8）。中层地下水中F 最大浓度高达13.71mg/L（M18），平均值1.35mg/L；碘浓度最大达1.708mg/L，平均值0.396mg/L。

中层地下水最大聚类特征见图4。为了划分最多簇数，同型种线分别定义为Y=1.5（绿色虚线，图4）以及Y=5.5（氟分界线，据氟含量为1.0mg/L 确定，红色虚线，图4），中层地下水最大聚类M1 可再划分为19 个簇，包括三大聚类：高碘样本（M1-11、M1-13、M1-18、M1-19），高氟和碘样本（M1-2、M1-3、M1-4、M1-5、M1-6、M1-7、M1-8、M1-9、M1-10、M1-12、M1-14、M1-15、M1-16、M1-17），以及正常样本（M1-1）。

图4 M1 聚类树状图（红色与绿色虚线为同型种线，对应Y 分别为5.5 和1.5），中层地下水最大聚类M1 中高氟高碘簇数目最多（M1-2、M1-3、M1-4、M1-5、M1-6、M1-7、M1-8、M1-9、M1-10、M1-12、M1-14、M1-15、M1-16、M1-17）Fig. 4 The tree diagram of M1 clusters (the red dashed line Y=5.5 and green dashed line Y=1.5 are phenon lines). The largest cluster of mesopelagic groundwater M1 has the largest number of high-fluoride and high-iodine clusters (M1-2, M1-3, M1-4, M1-5,M1-6, M1-7, M1-8, M1-9, M1-10, M1-12, M1-14, M1-15, M1-16, M1-17).

3.2.3 深层地下水

图5 显示了深层地下水的HCA 结果，同型种线分别定义为Y=1.5（绿色虚线，图5）和Y=5.5（据氟含量为2.0mg/L 确定，红色虚线，图5）。据此，所有深层地下水样本可划分为10 个簇，主要分为三大聚类：高氟样本（D9）；高碘样本（D3、D8）；高氟和碘样本（D1、D2、D4、D5、D6、D7、D10）。深层地下水中，氟浓度平均值为1.92mg/L，最大值为5.28mg/L（D7）；碘浓度平均值为0.581mg/L，最大值达4.601mg/L（D1）。

图5 深层地下水样本树状图（红色与绿色虚线为同型种线，对应Y 分别为5.5 和1.5），深层地下水高氟高碘簇数目最多（D1、D2、D4、D5、D6、D7、D10），最大氟浓度5.3mg/L（D7），最大碘浓度4.6mg/L（D1）Fig. 5 The tree map of medium groundwater samples (the red dashed line Y=5.5 and green dashed line Y=1.5 are phenon lines). Deep groundwater has the largest number of clusters of high fluoride and high iodine (D1, D2, D4, D5, D6, D7, D10); the maximum fluorine concentration is 5.3 mg/L (D7), and the maximum iodine concentration is 4.6 mg/L (D1).

4 讨论

4.1 高氟地下水的成因

HCA 结果显示：35 个浅层地下水群中有19 个高氟地下水群，占比达54.3%；37 个中层地下水群中具有28 个高氟地下水群，占比达75.6%；10 个深层地下水群中具有8 个高氟地下水群，占比达80%。说明鲁西北平原地下水系统中高氟地下水比例高达67.1%。浅层高氟地下水中，氟浓度为0.1～3.8mg/L，平均值为1.8mg/L；中层高氟地下水中，氟浓度为0.2～13.7mg/L，平均值为1.4mg/L；深层高氟地下水中，氟浓度为0.4～5.3mg/L，平均值为1.9mg/L。

白云母、黑云母、磷灰石等羟基矿物中富含氟，矿物风化（例如含氟矿物萤石CaF2的溶解）被认为是地下水中高浓度氟的主要来源［28］。鲁西北平原第四纪沉积层的岩土结构松散，厚度大，黄河冲积而下携带的大量黄土和其他物质构成了沉积层，黄土中含有大量的萤石矿物。因此，第四纪黄河沉积物是鲁西北平原地下水中氟的主要来源。

此外，地球化学过程也是控制地下水高氟的重要因素。鲁西北平原降水相对较少，蒸发量大，地势平坦，海拔基本在100m 以下，水力梯度小，地下水流动速度（平均为0.013cm/d）相对较慢，水循环交替时间较长，而且平原区地下水位埋深较浅（平均值为4.95m），岩性较细，有利于毛细水上升，这些特征都暗示研究区存在强烈的蒸发浓缩作用，使得地下水中氟含量不断升高。当含水层接受周围地表和地下水的补给时，周围山区广泛氟源所带来的含氟物质会随着地下水流运移而不断聚集。此外，萤石矿物风化后产生的氟也会随降雨进入土壤和地下水中。利用PHREEQC 计算萤石饱和指数（图6a），结果表明所有地下水样品均对萤石未饱和。而且通过地下水中氟与钙的关系（图6b），计算地下水中氟与钙的相关系数为−0.178，其中深层地下水中氟与钙的相关系数为−0.694，说明地下水中氟与钙呈现弱负相关的关系，尤其是深层地下水。萤石作为一种同时含有钙与氟元素的矿物，在天然条件下的水体中会同时存在着以下两种反应［29］：

图6 地下水氟、碘与其他元素相关图，地下水对萤石均未饱和，氟、钙含量呈弱负相关特别是深层地下水，碘含量随TDS 的增加而增加Fig. 6 Correlation diagram of fluorine, iodine and other elements in groundwater. The groundwater was unsaturated with fluorite; fluorine and calcium were weakly negatively correlated especially in deep groundwater; iodine content increased with the increase

通过这两个反应式可以看出，氟离子的分布受水文地球化学作用的影响十分显著。碱性且钙浓度低的重碳酸盐型水对含氟矿物的溶解起促进作用。随着地下水中pH 值升高，浓度增加，式（1）向右进行，Ca2+减少，使得式（2）中的反应向萤石溶解的方向进行移动，导致氟的富集［30-31］。因此，研究区高氟地下水往往赋存在低钙环境中［32-34］，这侧面印证萤石溶解度制约了地下水中的氟浓度。从鲁西北平原高氟piper 三线图（图7）可看出，浅层、中层、深层高氟地下水的水化学类型分别为HCO3·SO4·Cl-Na·Mg 型、SO4·Cl-Na·Mg 型、HCO3-Na 型水。HCO3-Na型地下水可为地下水中氟的富集提供有利条件［5，34］。这是因为一方面，较高浓度的HCO3促使方解石沉淀，在较弱的萤石溶解度极限下能促进富氟矿物的溶解；另一方面，HCO3-Na 型地下水通常呈碱性，在高pH条件下，氟可以与黏土矿物或云母（如白云母和黑云母）中的OH 交换［35-38］，有利于氟的富集，这是中深层地下水高氟的原因之一。

图7 高氟、高碘地下水piper 三线图，浅层、中层、深层高氟地下水的水化学类型分别为HCO3·SO4·Cl-Na·Mg 型、SO4·Cl-Na·Mg 型和HCO3-Na 型，高碘地下水中主要阳离子为Na+和Ca2+，主要阴离子为Cl−和HFig. 7 Piper diagram of high F and high I groundwater. The hydrochemical types of shallow, middle and deep high-fluoride groundwater are HCO3· SO4· Cl-Na· Mg type, SO4· Cl-Na· Mg type, and HCO3-Na type, respectively. The main cations in high-iodine groundwater are Na+ and Ca2+, and the main anions are Cl− and H.

与深层高氟地下水相比，浅层高氟地下水具有更高浓度的Na 和Cl，研究区深层高氟地下水Na 和Cl 含量均值分别为1466.87mg/L 和2567.40mg/L，浅层高氟地下水Na 和Cl 含量均值分别为2556.14mg/L和4156.57mg/L，这可能与蒸散发和岩盐溶解有关。对于浅层地下水，蒸发蒸腾作用较强，在萤石未饱和情况下，Ca 和F 的浓度略有增加［39］。当方解石饱和时，会生成碳酸钙沉淀，降低钙的活性，这有利于含氟矿物的溶解，使得浅层地下水中萤石硅和氟的浓度逐渐增加，直至与萤石达到饱和［40］。因此，蒸散发是浅层地下水中高浓度氟的重要影响因素。地下水中的盐分在强烈蒸发作用下不断浓缩，氟化物因溶解度不同先后在土壤表层析出。当发生降雨时，降雨穿过包气带垂直入渗补给至地下水，水土之间发生溶滤、溶解、水化、水解、离子吸附等一系列（水文）地球化学过程，将蒸发沉淀在土壤中的氟盐以及土壤中本身含有的氟盐溶解或淋滤至地下水中，导致氟在地下水系统中的富集［9］。

4.2 高碘地下水的成因

高碘地下水样品HCA 结果表明：35 个浅层地下水群中有26 个高碘地下水群，占比74.3%；85%以上的深层地下水样品中碘含量超过0.1mg/L。浅层、中层和深层地下水中碘的最大浓度分别为1.920mg/L、1.708mg/L 和4.601mg/L。深层高碘地下水中呈片状分布，碘含量较浅层和中层高。鲁西北平原高碘区通常分布在河流两岸，如山东聊城市高碘地下水分布在京杭运河两岸，德州市高碘地下水大多数分布在马颊河以北［41］。

蒸发浓缩及沉积物淋溶是浅层高碘地下水的主要成因，海侵现象也可引起地下水高碘现象［17，42］。鲁西北平原沉积物中含有大量的海洋生物及有机质，在风化、雨水侵蚀作用下，沉积物中的碘进入地下水中，构成深层地下水碘的重要来源。鲁西北平原经历过数次黄河冲刷，表层土壤中的碘被带走，进入浅层地下水中的碘含量降低，这是浅层地下水碘含量低于深层地下水的原因。据鲁西北平原高碘地下水piper 三线图结果可知（图7），高碘地下水中主要阳离子为Na+和Ca2+，主要阴离子为Cl−和，呈弱碱性［43］，符合蒸发强烈地带的特征，地下水中的TDS 关系图显示（图6c），在浅层和中层地下水中，随TDS 的增加，碘的含量表现出增加的趋势［21］，印证了浅层高碘地下水形成的蒸发浓缩机制。鲁西北平原地下水埋深浅，地下水蒸发强烈，蒸发作用只排走水分，盐分仍保留在地下水中，水中各种组分在蒸发过程中不断浓缩，浓度不断增加，导致表层土壤盐渍化严重，潜水中碘浓度也随之增加。因此，潜水中碘富集受蒸发浓缩作用影响较大。此外，鲁西北平原地势平坦，深层含水层系统连通性差，地下水径流条件弱，导致碘在地下水中富集［17］；并且高碘地下水中、浓度普遍较低，浅层和深层高碘地下水NO3−浓度均值分别为1.42mg/L、1.65mg/L，低于研究区正常地下水浓度（12.94mg/L，4.57mg/L）。浓度较高，浅层和深层高碘地下水浓度均值分别为670.43mg/L 和351.70mg/L，高于浅层和深层高碘地下水浓度均值（598.50mg/L，229.15mg/L）。/（+）可作为氧化还原环境敏感因子判断水环境的性质［36］，由于地下水含量极低，该因子可简化为，该比值越小指示水环境还原性越强。据此计算得到研究区地下水尤其是深层地下水处于强还原环境中（图8）。大量研究表明，还原条件有利于地下水中碘的富集［21，44］。这是深层地下水中碘富集的主要原因。

图8 随深度变化关系图，地下水尤其是深层地下水处于强还原环境中Fig. 8 versus depth. With the increase of depth, the groundwater is more reductive.

4.3 高氟地下水与高碘地下水之间的关系

地下水氟与碘含量的空间分布图与垂向分布图分别见图9 和图10。从图9 可看出，地下水氟含量高的区域碘的含量通常也很高，高氟、碘区主要位于鲁西北的滨海平原地区。对于浅层和中层地下水，随深度增加，氟与碘浓度逐渐减小（图10）。高氟地下水主要集中在地面以下0～40m 与50～110m 深度内，氟含量最大值（13.7mg/L）出现在110m 深度处；高碘地下水主要集中在地面以下0～10m 与55～65m 深度内，碘含量最大值（4.6mg/L）出现在234m深度处。深层高氟地下水主要集中在地面以下225～305m 范围内；深层地下水碘含量普遍较高，主要集中在地面以下225～285m 范围内。

图9 地下水氟和碘含量空间分布图，高氟、碘区主要位于鲁西北的滨海平原地区Fig. 9 Spatial distribution of F and I in groundwater. The high fluoride and iodine areas are mainly located in the coastal plain area of Northwest Shandong.

图10 地下水氟和碘含量垂向分布图，浅层、中层、深层高氟高碘地下水分别主要集中在地面以下0～40m、50～110m、225～305mFig. 10 Vertical distribution of F and I in groundwater of Northwest Shandong Plain. The groundwater with high fluorine and iodine in shallow, mediate and deep aquifers is mainly concentrated in 0-40m, 50-110m and 225-305m below the ground surface, respectively.

以上结果表明地下水中氟与碘通常伴生出现，即高氟地下水中，碘含量也较高。这是因为利于二者的富集环境条件较为相似。强烈的蒸发浓缩及淋溶作用是导致鲁西北平原高氟、碘地下水的主要过程。研究区蒸发浓缩作用强烈，埋深浅的浅层地下水中的含氟矿物和含碘沉积物不断溶解，使地下水中氟、碘含量升高；发生降雨时，土壤中的氟、碘溶解或淋滤至地下水中，导致其富集。另外研究区呈现出碱性、高的特征，是促进含氟矿物溶解和碘富集的良好背景条件。研究区深层地下水流动速度慢，利于氟富集，深层地下水还处于还原环境中，存在NH4-N 和HS−的还原溶解，利于碘富集。高氟地下水与高碘地下水水化学类型近似，主要阴离子为、Cl−，主要阳离子为Na+，且均呈弱碱性（pH=7～10）。在弱碱性环境下，通过溶滤作用、离子交换作用形成高氟、碘地下水。

5 结论

本文运用层次聚类分析方法将采集到的地下水样本进行分类，探讨鲁西北平原高氟、高碘地下水的分布特征及形成机理。研究结果表明：①高氟地下水主要集中在地面以下0～40m、50～110m 以及225～305m 深度内，氟含量最大值（13.71mg/L）出现在110m 深度处；高碘地下水主要集中在地面以下0～10m、55～65m 以及225～285m 深度内，碘含量最大值（4.601mg/L）出现在234m 深度处。②含氟矿物萤石（CaF2）的溶解、强烈的蒸发浓缩作用、离子交换可能是鲁西北平原高氟地下水的主要形成机制；蒸发浓缩、沉积物淋溶以及氨氮和硫化物的还原溶解是可能导致鲁西北平原地下水高碘的主要过程。③高氟地下水与高碘地下水均呈弱碱性（pH=7～10），还原环境，且水化学类型类似，主要阳离子为Na+，主要阴离子为、Cl−。高氟地下水与高碘地下水通常伴生出现，这是由于两者富集所需的环境条件相似。

本文主要探究了高氟地下水与高碘地下水由于具有相似的利于富集的水化学环境而在鲁西北平原存在的伴生情况，为中国其他同类型地区的研究提供一定借鉴作用。下一步可进一步探索其他地区高氟地下水和高碘地下水形成不同形态之间时是否存在相互促进的关系及其机制。