CN113670879B - 基于多变量和pmf模型对黑臭水体底泥重金属源的分析方法 - Google Patents

基于多变量和pmf模型对黑臭水体底泥重金属源的分析方法 Download PDF

Info

Publication number
CN113670879B
CN113670879B CN202110996501.XA CN202110996501A CN113670879B CN 113670879 B CN113670879 B CN 113670879B CN 202110996501 A CN202110996501 A CN 202110996501A CN 113670879 B CN113670879 B CN 113670879B
Authority
CN
China
Prior art keywords
heavy metal
sampling
sediment
black
multivariate
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN202110996501.XA
Other languages
English (en)
Other versions
CN113670879A (zh
Inventor
范婷婷
邓绍坡
杨敏
孔令雅
彭立
姜登登
李旭伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Institute of Environmental Sciences MEE
Original Assignee
Nanjing Institute of Environmental Sciences MEE
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Institute of Environmental Sciences MEE filed Critical Nanjing Institute of Environmental Sciences MEE
Priority to CN202110996501.XA priority Critical patent/CN113670879B/zh
Publication of CN113670879A publication Critical patent/CN113670879A/zh
Priority to JP2022128543A priority patent/JP7227439B1/ja
Application granted granted Critical
Publication of CN113670879B publication Critical patent/CN113670879B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64—Fluorescence; Phosphorescence
    • G01N21/6402—Atomic fluorescence; Laser induced fluorescence
    • G01N21/6404—Atomic fluorescence
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Health & Medical Sciences (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Software Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Optics & Photonics (AREA)
  • Operations Research (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Databases & Information Systems (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Sampling And Sample Adjustment (AREA)
  • Water Treatment By Sorption (AREA)

Abstract

本发明公开了一种基于多变量和PMF模型对黑臭水体底泥重金属源的分析方法，包括以下步骤：S1样品采集；S2样品检测；S3重金属污染程度评价；S4潜在生态危害评价；S5重金属溯源分析：S5‑1多变量统计分析，S5‑2正定矩阵因子分析。步骤S1中取样装置为活塞式柱状取样器，包括采样杆以及套设在采样杆下部外侧的采样桶。本发明黑臭水体底泥重金属源的分析方法对底泥中的重金属生态环境风险开展评价，并通过多变量统计分析与正定矩阵因子模型相结合的方法对黑臭水体底泥中的重金属开展定量化的溯源分析，从而对河流生态环境的改善提出针对性的建议和措施，为城市黑臭河道水体治理标本兼治、生态恢复提供重要保障。

Description

基于多变量和PMF模型对黑臭水体底泥重金属源的分析方法
技术领域
本发明涉及黑臭水体治理技术领域，具体是涉及一种基于多变量和PMF模 型对黑臭水体底泥重金属源的分析方法。
背景技术
城市河流是城市生态环境的重要载体，对城市防洪排涝、水路运输和景观 美化等方面起到了重要作用。但由于其相对封闭，自净能力及水体更新速度不 及大江大河，因而易受沿岸排放的污染物影响较为严重。由于大量污染物的排 放入水体，水体中化学需氧量(COD)、氮(N)、磷(P)等污染物浓度超标， 河流水体污染严重，水体出现季节性或终年黑臭。而底泥通常是河流中污染物 的最终归宿，通过各种途径进入水体的污染物，在一系列物理、化学、生物作 用下沉积到底泥中。同时，当固液界面的性质发生变化时，“储蓄库”——底 泥中的污染物可能将再次释放出来从而对水生生态系统产生危害。
黑臭水体中的污染物不仅是氮磷营养盐、还包括重金属以及其他持久性有 机污染物。由于重金属具有较强的生物毒性，且在环境中能够持久稳定保持， 以及生物放大效应强，底泥中的重金属污染日益受到关注。国内外学者对城市 河道底泥中重金属的累积程度及潜在生态环境风险都开展了评价和研究。不同 城市河道底泥中重金属的生态风险表现不同，例如，在杭州市区的河道底泥中 重金属的潜在生态风险顺序为Hg>Cd>Pb>Cu>As>Zn>Ni>Cr，其中Hg和Cd是 主要的潜在生态风险因子；在江苏省的重金属防控区域内河道底泥中重金属的 潜在生态风险依次为Cd>Hg≈Pb>As；南通市区河道底泥中重金属的潜在生态 风险则表现为Hg>Cu>Zn>Cr>Pb；合肥市区河道底泥中重金属的生态危害顺序 为Hg>Cd>Cu>As>Zn>Pb>Cr。
底泥重金属治理也是黑臭河道水体整治的重要内容，是保障河道生态系统 恢复健康的重要保障，而对底泥重金属开展溯源分析对于黑臭水体治理及其区 域性环境管理政策的提出具有重要意义。底泥中的重金属主要来自自然源和人 为源。基岩的化学淋溶作用、流域的排水作用、河岸的地表径流等属于自然源。 城市和工业废水的排放、矿物燃料的燃烧、采矿和冶炼、加工制造业以及包括 倾倒在内的废物处理等属于人为源。尤其是城市环境基质(如大气、土壤和水体) 中的重金属多数来自人为源。
目前，同位素示踪是重金属污染源解析较为有效的方法之一，但是由于其 测定条件的限制，难以大规模广泛应用。因此需要提出更为有效地重金属污染 源解析方法。
发明内容
针对上述存在的问题，本发明提供了一种基于多变量和PMF模型对黑臭水 体底泥重金属源的分析方法。
本发明的技术方案是：
基于多变量和PMF模型对黑臭水体底泥重金属源的分析方法，包括以下步 骤：
S1样品采集：取若干河道断面作为采样点，在每个河道断面的两端和中间 分别通过取样器对河道底泥进行采样，将采集到的底泥样品风干研磨并过2mm 筛网筛分，得到待测底泥颗粒；
S2样品检测：将待测底泥颗粒中的Hg、As、Cu、Ni、Cr、Zn、Cd和Pb 含量通过HCl-HNO3-HF-HClO4电热板消解原子吸收分光光度计进行测定，待测 底泥颗粒中的Hg和As元素采用王水水浴消解提取，随后通过原子荧光光谱仪 测定；
S3重金属污染程度评价：通过地累积指数法作为评价底泥中重金属污染程 度的指标；
S4潜在生态危害评价：通过潜在生态危害指数评价底泥污染程度及其潜在 生态危害；
S5重金属溯源分析：
S5-1多变量统计分析：通过Pearson相关系数矩阵分析底泥中重金属含量 之间的相关性，说明不同重金属来源的一致性，利用主成分分析法PCA通过降 维的思维，从原始变量集中提取一部分称为主成分的潜在因素，分析底泥中重 金属的一级来源；
S5-2正定矩阵因子分析：利用正定矩阵因子分析底泥中重金属的二级来源， 正定矩阵因子分析PMF的基本方程如下：
Figure BDA0003234300780000031
Figure BDA0003234300780000032
其中，Xij为第i个样品中第j种元素的含量，Gik为第k个因子对第i个样 品的贡献率，Fkj为第k个因子中第j种元素的含量，Eij为第i个样品中第j种元 素的残差，p为因子个数，PMF模型运用样品含量及不确定度uij进行迭代计算 得到较小的Q值，用于评价选定的因子数量是否合理，以得到最优的因子贡献 矩阵G和因子成分矩阵F的解。
进一步地，所述步骤S3中地累积指数法的计算公式为：
Figure BDA0003234300780000033
其中，Igeo为重金属地累积指数，Ci为实测底泥中重金属质量分数，Bi为普 通页岩中重金属元素的地球化学平均背景值，k为考虑当地岩石差异可能引起 的背景值的变动而取的变异系数取1.5；依据Igeo值将底泥中重金属污染状况分 级，广泛用于研究沉积物中重金属污染程度的定量指标。
进一步地，所述步骤S4中潜在生态危害指数的计算公式为：
Eri＝Tri×Cdi/Cri
RI＝∑Eri
其中，Eri为单个重金属的潜在生态危害系数，Tri为Hg、As、Cu、Ni、Cr、 Zn、Cd和Pb的毒性响应系数，分别为40、10、5、5、2、1、30和5，Cdi为样 品浓度实测值，Cri为计算所需的参照值，RI为多重金属潜在生态危害指数，潜 在生态危害指数是用于评价底泥污染程度及其潜在生态危害的一种相对快速、 简便和标准的方法。
更进一步地，所述Cri参照值采用工业化以前底泥中重金属的最高背景值， 以GB15618-1995自然背景值作为参比，具有良好的对比性。
进一步地，所述步骤S5-1中PCA的有效性是用Kaiser-Meyer-Olkin和 Bartlett球度p进行检验，其中KMO>0.5，p<0.001，并根据分解的特征值建立 协方差矩阵，根据重金属的实测值和特征值N1的贡献率，采样方差旋转法选 择主成分。
更进一步地，所述步骤S5-2中使用EPA PMF 5.0对重金属来源进行分析， 使用20％MDL来计算不确定度，因子成分矩阵F包括Factor1自然母质源、 Factor2交通源和大气降尘的综合排放、Factor3生活污水和初期雨水的综合排放、 Factor4工业排放源，正定矩阵因子分析PMF是一种通过受体模型数学方法利 用样本组成或者指纹对污染物源进行定量分析的一种方法。
进一步地，所述步骤S1中取样装置为活塞式柱状取样器，方便取样工作进 行。
更进一步地，所述活塞式柱状取样器包括采样杆以及套设在所述采样杆下 部外侧的采样桶，所述采样桶上端两侧各设有一组延伸至采样杆上端的固定板， 采样杆中下部外壁设有螺纹，位于螺纹下端的采样杆外壁上对称设有两组弹簧 板，采样桶内壁上部设有用于与螺纹啮合连接的螺纹部，所述螺纹部下方的采 样桶内壁对称设有两组收集部，所述收集部下端设有开口，位于所述开口处的 采样桶内壁设有用于与所述弹簧板对接的卡槽，位于弹簧板下方的采样杆内部 设有传感器以及PLC控制器，采样杆下端设有钻头，采样杆内底部设有分别用 于控制所述钻头转动的伺服电机，钻头上部设有若干采样口，可以完成表层底 泥和深层底泥的同步取样。
优选地，所述固定板上端设有挡板，所述收集部上端设有滤水孔，位于收 集部上方的采样桶内壁设有限位块，每组所述采样口处均设有的采样板，所述 采样板与钻头通过弹簧连接，可自动完成二次采样过程。
与现有技术相比，本发明的有益效果是：
(1)本发明黑臭水体底泥重金属源的分析方法对底泥中的重金属生态环境 风险开展评价，并通过多变量统计分析与正定矩阵因子模型相结合的方法对黑 臭水体底泥中的重金属开展定量化的溯源分析，从而对河流生态环境的改善提 出针对性的建议和措施，为城市黑臭河道水体治理标本兼治、生态恢复提供重 要保障。
(2)本发明黑臭水体底泥重金属源的分析方法中多变量统计分析方法为分 析污染物来源、区分自然源和人为源的贡献以及指示转运过程和环境条件提供 了一种可靠地方法，而正定矩阵因子分析PMF是一种通过受体模型数学方法利 用样本组成或者指纹对污染物源进行定量分析的一种方法，多变量统计分析和 PMF模型联用作为可靠、便捷、经济的方法可应用于黑臭水体底泥重金属源解 析。
(3)本发明的底泥取样装置能够对河道的表层底泥以及深层底泥进行准确 收集，可根据需要通过取样器先后对表层底泥以及深层底泥进行自动采样，使 用方便，取样快速，使得实验结果更加可靠。
附图说明
图1是本发明实验例中不同采样位置河道底泥重金属含量分布特征；
图2是本发明实验例中底泥中8种重金属元素的地累积指数分布图；
图3是本发明实验例中重金属潜在生态危害综合指数频率分布图和各元素 对潜在生态危害综合指数的贡献率；
图4是本发明实验例中PCA二维因子载荷图；
图5是本发明实验例中重金属污染物的溯源分析图；
图6是本发明的取样器整体结构示意图；
图7是本发明的取样器底部结构示意图；
图8是本发明的取样器内部结构示意图；
图9是本发明的取样器仰视图；
图10是本发明的取样器采样口处的结构示意图。
其中，1-采样杆，11-传感器，12-PLC控制器，13-伺服电机，2-采样桶， 21-螺纹部，22-收集部，23-卡槽，24-滤水孔，25-限位块，3-固定板，31-挡板， 4-弹簧板，5-钻头，51-采样口，52-采样板。
具体实施方式
实施例1
以靖江城市河道为例，基于多变量和PMF模型对该河道某段黑臭水体底泥 重金属源的分析方法，包括以下步骤：
S1样品采集：取27个河道断面作为采样点，在每个河道断面的两端和中 间分别通过取样器对河道底泥进行采样，共采集81份样品，将采集到的底泥样 品风干研磨并过2mm筛网筛分，得到待测底泥颗粒；
S2样品检测：将待测底泥颗粒中的Hg、As、Cu、Ni、Cr、Zn、Cd和Pb 含量通过HCl-HNO3-HF-HClO4电热板消解原子吸收分光光度计进行测定，待测 底泥颗粒中的Hg和As元素采用王水水浴消解提取，随后通过原子荧光光谱仪 测定；
S3重金属污染程度评价：通过地累积指数法作为评价底泥中重金属污染程 度的指标，地累积指数法的计算公式为：
Figure BDA0003234300780000061
其中，Igeo为重金属地累积指数，Ci为实测底泥中重金属质量分数，Bi为普 通页岩中重金属元素的地球化学平均背景值，k为考虑当地岩石差异可能引起 的背景值的变动而取的变异系数取1.5；依据Igeo值将底泥中重金属污染状况分 级；
S4潜在生态危害评价：通过潜在生态危害指数评价底泥污染程度及其潜在 生态危害，潜在生态危害指数的计算公式为：
Eri＝Tri×Cdi/Cri
RI＝∑Eri
其中，Eri为单个重金属的潜在生态危害系数，Tri为Hg、As、Cu、Ni、Cr、 Zn、Cd和Pb的毒性响应系数，分别为40、10、5、5、2、1、30和5，Cdi为样 品浓度实测值，Cri为计算所需的参照值，Cri参照值采用工业化以前底泥中重金 属的最高背景值，以GB15618-1995自然背景值作为参比，RI为多重金属潜在 生态危害指数；
S5重金属溯源分析：
S5-1多变量统计分析：通过Pearson相关系数矩阵分析底泥中重金属含量 之间的相关性，说明不同重金属来源的一致性，利用主成分分析法PCA通过降 维的思维，从原始变量集中提取一部分称为主成分的潜在因素，分析底泥中重 金属的一级来源，PCA的有效性是用Kaiser-Meyer-Olkin和Bartlett球度p进行 检验，其中KMO>0.5，p<0.001，并根据分解的特征值建立协方差矩阵，根据 重金属的实测值和特征值N1的贡献率，采样方差旋转法选择主成分；
S5-2正定矩阵因子分析：利用正定矩阵因子分析底泥中重金属的二级来源， 正定矩阵因子分析PMF的基本方程如下：
Figure BDA0003234300780000071
Figure BDA0003234300780000072
其中，Xij为第i个样品中第j种元素的含量，Gik为第k个因子对第i个样 品的贡献率，Fkj为第k个因子中第j种元素的含量，Eij为第i个样品中第j种元 素的残差，p为因子个数，PMF模型运用样品含量及不确定度uij进行迭代计算 得到较小的Q值，用于评价选定的因子数量是否合理，以得到最优的因子贡献 矩阵G和因子成分矩阵F的解，因子成分矩阵F包括Factor1自然母质源、 Factor2交通源和大气降尘的综合排放、Factor3生活污水和初期雨水的综合排放、 Factor4工业排放源，使用EPA PMF 5.0对重金属来源进行分析，使用20％ MDL来计算不确定度。
实施例2
本实施例与实施例1基本相同，其不同之处在于，步骤S1中样品采集个数 不同。
S1样品采集：取23个河道断面作为采样点，在每个河道断面的两端和中 间分别通过取样器对河道底泥进行采样，共采集69份样品，将采集到的底泥样 品风干研磨并过2mm筛网筛分，得到待测底泥颗粒。
实施例3
本实施例与实施例1基本相同，其不同之处在于，步骤S1中样品采集个数 不同。
S1样品采集：取31个河道断面作为采样点，在每个河道断面的两端和中 间分别通过取样器对河道底泥进行采样，共采集93份样品，将采集到的底泥样 品风干研磨并过2mm筛网筛分，得到待测底泥颗粒。
实施例4
本实施例与实施例1基本相同，其不同之处在于：步骤S1中取样装置为活 塞式柱状取样器，如图6所示。
如图7-10所示，活塞式柱状取样器包括采样杆1以及套设在采样杆1下部 外侧的采样桶2，采样桶2上端两侧各设有一组延伸至采样杆1上端的固定板3， 固定板3上端设有挡板31，采样杆1中下部外壁设有螺纹，位于螺纹下端的采 样杆1外壁上对称设有两组弹簧板4，采样桶2内壁上部设有用于与螺纹啮合 连接的螺纹部21，螺纹部21下方的采样桶2内壁对称设有两组收集部22，收 集部22上端设有滤水孔24，收集部22下端设有开口，位于开口处的采样桶2 内壁设有用于与弹簧板4对接的卡槽23，位于收集部22上方的采样桶2内壁设 有限位块25，位于弹簧板4下方的采样杆1内部设有传感器11以及PLC控制 器12，传感器11以及PLC控制器12均为市售，采样杆1下端设有钻头5，采 样杆1内底部设有分别用于控制钻头5转动的伺服电机13，伺服电机13为市售 伺服电机，钻头5上部设有若干采样口51，每组采样口51处均设有的采样板 52，采样板52与钻头5通过弹簧连接。采样杆1和固定板3为多节拼接而成。
应用上述活塞式柱状取样器进行底泥取样的工作原理为：
首先根据河道深度通过拼接采样杆1和固定板3选择适当的取样器长度， 将取样器下放至河道底部，当接触到底部时，通过按压挡板31使收集部22下 沉进入表层底泥中，将表层底泥收集至收集部22内，同时通过滤水孔24过滤 掉多余的水，当收集部22完全进入底泥中时，使用管钳等工具转动采样杆1顶 端的把手，使采样杆1旋转，在螺纹部21的作用下采样杆1螺纹转动下降，弹 簧板4同步下降但是在收集部22的阻挡下不会弹出，直至收集部4下降至卡槽 23对应的位置，弹簧板4弹出，并与卡槽23卡接，此时传感器11感应到压力信号后传递至PLC控制器12控制伺服电机13驱动钻头5先逆时针转动后顺时 针转动，当钻头5逆时针转动时，采样板52在弹簧的作用下打开，深层底泥由 采样口51进入到钻头5内部储物空间内，当钻头5顺时针转动时，采样板52 关闭，随后通过挡板提起取样器完成取样。
实验例
选取实施例1中采集到的数据进行分析，对本发明方法的步骤进行详细解 释说明，河道底泥中重金属含量及部分理化性质统计结果如表1所示。
表1河道底泥重金属含量描述性统计(mg/kg)
Figure BDA0003234300780000091
Figure BDA0003234300780000101
从底泥中8种重金属含量的平均值而言，表现为 Pb(902.21mg/kg)>Zn(506.52mg/kg)>Cr(118.06mg/kg)>Cu(78.89mg/kg)>Ni(54.97m g/kg)>As(15.05mg/kg)>Hg(0.33mg/kg)≈Cd(0.30mg/kg)；从变异系数而言，高度 变异的重金属元素有Zn、Hg，变异系数大于36％，中度变异的有Pb、Cu、Cr， 变异系数在15-36％，低度变异的有As、Ni、Cd，变异系数小于15％，变异系 数可以反应重金属含量的波动性，从而对其受人为因素影响程度进行分析，变 异系数越大，重金属污染受人为影响的程度越大，由此表明在河道底泥中Zn、 Hg、Pb、Cu、Cr等有可能受到较强的人为因素的影响。
河道不同位置(A、C：两端，B：中间)的底泥中8种重金属含量分布特 征如图1所示，不同位置的底泥中重金属含量变化范围与总体情况具有一致规 律，通过显著性分析结果表面3个不同采样位置中重金属含量未表现出显著性 差异。
为了更好的反应底泥中重金属的污染情况，计算不同样品中8种重金属的 Igeo指数，结果如图2所示，在8种重金属中，Pb的污染程度相对较重，95.06％ 的样品达到中度污染水平，4.94％的样品呈现低度污染水平；其次是Zn、Cu和 Hg，在低度污染水平的样品占比76.54％～90.12％。
如表2和图3所示，河流底泥中重金属的潜在生态风险总体属于中等风险 和强风险，比例分别为59.3％和40.7％。其中，Pb对底泥中潜在生态危害指数 贡献最高的元素，其次是Hg和Cd，综上考虑重金属生物毒性系数及其对生态 风险的影响，在该河道中应对Pb和Hg着重关注。
表2河道底泥重金属生态风险等级指数
Figure BDA0003234300780000111
通过KMO值和Bartlett的球形度检验表明可以利用PCA方法对底泥中的 重金属开展主成分分析，通过主成分分析，共提取8个成分即互相独立的综合 指标，如表3所示。其中，前2个主成分的特征值大于1，且累积贡献率已经 达到70.33％，即提取的成分已经足以解释大部分原始影响因素。因此，本次主 成分成功提取了2个主成分。
表3总方差的分解表
Figure BDA0003234300780000112
如图4所示，给出了每个主成分对于8种重金属的贡献率，主成分1对As、 Cd、Pb、Hg和Cu具有较高的正载荷，其载荷值分别为0.884、0.745、0.764、 0.758和0.666，如表4所示，相关性结果表明Pb与As、Cu、Hg之间具有较高 的正相关，也说明以上四个元素的来源一致的可能性相对较高，推测主成分1 为人为源；主成分2对Zn、Ni、Cr具有较高的正载荷值，其载荷值分别为 0.719、0.72和0.881，此外，还发现主成分2也对Cu、Pb和Hg均具有一定的 正载荷，均在0.3以上，推测主成分2为自然源。
表4底泥中各元素之间的Pearson相关系数
元素 Hg As Cu Ni Cr Zn Cd Pb
As .724** - - - - - - -
Cu .748** .697** - - - - - -
Ni .304** .310** .569** - - - - -
Cr .445** .312** .625** .606** - - - -
Zn .291** 0.203 .586** .236* .505** - - -
Cd .395** .488** .346** .226* 0 0.02 - -
Pb .648** .763** .802** .365** .413** .397** .368** -
本发明采用EPA PMF 5.0软件进行源解析计算，尝试3-7种源成分谱进行 模型，在实际解析中3因子的Q值为157，4因子的Q值为84.1，而5因子的 Q值为37.9，结合研究区可能的源成分谱信息，模型最终运行结果采用4因子 模式，此时目标函数Q值已较低，残差数值大部分在-3～3之间，且除Cd外模 型对其他7种重金属的预测R2均在0.6以上，尤其是Pb和Zn的预测R2甚至 达到了0.95。
PMF模型对底泥中各金属的源解析结果如表5所示。
表5 PMF 5.0模型计算的各重金属污染源贡献率(％)
元素 Factor 1 Factor 2 Factor 3 Factor 4 自然源 人为源
Hg 19.35 35.92 44.73 0.00 19.35 80.65
As 13.55 16.14 36.77 33.54 13.55 86.45
Cu 24.92 29.56 21.17 24.35 24.92 75.08
Ni 18.86 2.28 45.05 33.81 18.86 81.14
Cr 29.62 1.79 47.69 20.90 29.62 70.38
Zn 48.29 21.33 0.03 30.35 48.29 51.71
Cd 13.11 7.40 41.66 37.84 13.11 86.9
Pb 9.23 45.31 0.81 44.65 9.23 90.77
8种重金属在因子1上具有一定的贡献率，除Zn外整体贡献率并不特别突 出，说明因子1不是导致河道底泥中重金属的主要来源。各种重金属元素都会 在自然界中存在，经岩石风化成土作用形成的土壤母质是它们的共同来源，As、 Cr、Cu、Pb、Zn等重金属含量均较低总体呈现非污染水平，因此该来源不会导 致重金属元素浓度的超标。因此推测因子1为自然母质源。因子1对Zn的贡献 率最高48.29％，这也与主成分的分析结果较为一致，Zn主要由主成分1贡献。
Pb、Hg、Cu、Zn在因子2有较高的浓度值，其中Pb的贡献率高达45.31％。 Pb为交通运输的主要标志物，交通源中Cr主要来自气溶胶、Cu主要来自制动 设备的磨损、Zn主要来自汽车润滑油使用及金属分解。此外，Hg主要来源是 煤的燃烧，燃煤对Hg累积的影响十分显著，长期燃煤发电产生的Hg降尘对河 道底泥中的累积产生一定影响。因此，因子2可以认为是交通源和大气降尘的 综合排放。
除Pb和Zn外，其他元素Cr、Ni、Hg、Cd、As、Cu在因子3上均有较高 的贡献率。生活污水和初期雨水逐渐成为城市污水中重金属的主要贡献者，该 河道位于生活聚集区，河道内存在生活废水和初期雨水的排放口且河道一侧存 在农户的自留地。Cr主要来源于润滑油、油漆颜料、合金、农药等。Cd主要来 自清洁剂、化妆品、食品加工及化肥等。Ni主要来自卫浴设备、电池、油漆等、 Hg主要来自温度计、牙科、颜料、玻璃镜等；As主要来自洗涤用品、防腐剂、 油漆等。因此，因子3可以认为是生活污水和初期雨水的综合排放。
如图5所示，Pb、Zn、Cd、Cr、Ni、Cu、As在因子4上具有较高的贡献 率，均在20％以上，尤其是Pb的贡献率最高，达到了44.65％。据前期调查表 明该河道周边1km内存在铅盐生产企业、冶金企业和零星的电镀企业，由于该 工业源的废水和废气等的影响，且河水流速较慢，自净能力较弱，重金属在底 泥中行程积累，造成重金属污染。因此，因子4可以解释为工业排放源。
对比多变量统计分析和PMF模型的源解析结果可以发现，对于累积程度较 高且生态风险高的重金属(Pb、Hg、Cu、Cd)元素，多变量统计方法和PMF模 型在来源解析上具有较高的一致性。因此，多变量统计分析和PMF模型联用作 为可靠、便捷、经济的分析方法用于黑臭水体底泥重金属溯源，为黑臭水体标 本兼治提供重要保障。

Claims (7)

1.基于多变量和PMF模型对黑臭水体底泥重金属源的分析方法，其特征在于，包括以下步骤：
S1样品采集：取若干河道断面作为采样点，在每个河道断面的两端和中间分别通过取样器对河道底泥进行采样，将采集到的底泥样品风干研磨并过2mm筛网筛分，得到待测底泥颗粒，取样装置为活塞式柱状取样器；
S2样品检测：将待测底泥颗粒中的Hg、As、Cu、Ni、Cr、Zn、Cd和Pb含量通过HCl-HNO3-HF-HClO4电热板消解，再用原子吸收分光光度计进行测定，待测底泥颗粒中的Hg和As元素采用王水水浴消解提取，随后通过原子荧光光谱仪测定；
S3重金属污染程度评价：通过地累积指数法作为评价底泥中重金属污染程度的指标；
S4潜在生态危害评价：通过潜在生态危害指数评价底泥污染程度及其潜在生态危害；
S5重金属溯源分析：
S5-1多变量统计分析：通过Pearson相关系数矩阵分析底泥中重金属含量之间的相关性，说明不同重金属来源的一致性，利用主成分分析法PCA，从原始变量集中提取一部分，分析底泥中重金属的一级来源；
S5-2正定矩阵因子分析：利用正定矩阵因子分析底泥中重金属的二级来源，正定矩阵因子分析PMF的基本方程如下：
Figure FDA0003740088070000011
Figure FDA0003740088070000012
其中，Xij为第i个样品中第j种元素的含量，Gik为第k个因子对第i个样品的贡献率，Fkj为第k个因子中第j种元素的含量，Eij为第i个样品中第j种元素的残差，p为因子个数，PMF模型运用样品含量及不确定度uij进行迭代计算得到较小的Q值，用于评价选定的因子数量是否合理，以得到最优的因子贡献矩阵G和因子成分矩阵F的解；
所述活塞式柱状取样器包括采样杆(1)以及套设在所述采样杆(1)下部外侧的采样桶(2)，所述采样桶(2)上端两侧各设有一组延伸至采样杆(1)上端的固定板(3)，采样杆(1)中下部外壁设有螺纹，位于螺纹下端的采样杆(1)外壁上对称设有两组弹簧板(4)，采样桶(2)内壁上部设有用于与螺纹啮合连接的螺纹部(21)，所述螺纹部(21)下方的采样桶(2)内壁对称设有两组收集部(22)，所述收集部(22)下端设有开口，位于所述开口处的采样桶(2)内壁设有用于与所述弹簧板(4)对接的卡槽(23)，位于弹簧板(4)下方的采样杆(1)内部设有传感器(11)以及PLC控制器(12)，采样杆(1)下端设有钻头(5)，采样杆(1)内底部设有用于控制所述钻头(5)转动的伺服电机(13)，钻头(5)上部设有若干采样口(51)。
2.根据权利要求1所述的基于多变量和PMF模型对黑臭水体底泥重金属源的分析方法，其特征在于，所述S3中地累积指数法的计算公式为：
Figure FDA0003740088070000021
其中，Igeo为重金属地累积指数，Ci为实测底泥中重金属质量分数，Bi为普通页岩中重金属元素的地球化学平均背景值，k为变异系数，取1.5；依据Igeo值将底泥中重金属污染状况分级。
3.根据权利要求1所述的基于多变量和PMF模型对黑臭水体底泥重金属源的分析方法，其特征在于，所述S4中潜在生态危害指数的计算公式为：
Eri＝Tri×Cdi/Cri
RI＝∑Eri
其中，Eri为单个重金属的潜在生态危害系数，Tri为Hg、As、Cu、Ni、Cr、Zn、Cd和Pb的毒性响应系数，分别为40、10、5、5、2、1、30和5，Cdi为样品浓度实测值，Cri为计算所需的参照值，RI为多重金属潜在生态危害指数。
4.根据权利要求3所述的基于多变量和PMF模型对黑臭水体底泥重金属源的分析方法，其特征在于，所述Cri参照值采用GB15618-1995中的自然背景值作为参比。
5.根据权利要求1所述的基于多变量和PMF模型对黑臭水体底泥重金属源的分析方法，其特征在于，所述S5-1中PCA的有效性是用Kaiser-Meyer-Olkin和Bartlett球度p进行检验，其中Kaiser-Meyer-Olkin>0.5，Bartlett球度p<0.001，并根据分解的特征值建立协方差矩阵，根据重金属的实测值和特征值N1的贡献率，采样方差旋转法选择主成分。
6.根据权利要求1所述的基于多变量和PMF模型对黑臭水体底泥重金属源的分析方法，其特征在于，所述S5-2中使用EPA PMF 5.0对重金属来源进行分析，使用20％MDL来计算不确定度，因子成分矩阵F包括Factor1自然母质源、Factor2交通源和大气降尘的综合排放、Factor3生活污水和初期雨水的综合排放、Factor4工业排放源。
7.根据权利要求1所述的基于多变量和PMF模型对黑臭水体底泥重金属源的分析方法，其特征在于，所述固定板(3)上端设有挡板(31)，所述收集部(22)上端设有滤水孔(24)，位于收集部(22)上方的采样桶(2)内壁设有限位块(25)，所述采样口(51)处均设有采样板(52)，所述采样板(52)与钻头(5)通过弹簧连接。
CN202110996501.XA 2021-08-27 2021-08-27 基于多变量和pmf模型对黑臭水体底泥重金属源的分析方法 Expired - Fee Related CN113670879B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202110996501.XA CN113670879B (zh) 2021-08-27 2021-08-27 基于多变量和pmf模型对黑臭水体底泥重金属源的分析方法
JP2022128543A JP7227439B1 (ja) 2021-08-27 2022-08-11 多変数とｐｍｆモデルに基づいて黒臭水域の底質の重金属源を分析する方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110996501.XA CN113670879B (zh) 2021-08-27 2021-08-27 基于多变量和pmf模型对黑臭水体底泥重金属源的分析方法

Publications (2)

Publication Number Publication Date
CN113670879A CN113670879A (zh) 2021-11-19
CN113670879B true CN113670879B (zh) 2022-09-13

Family

ID=78547050

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110996501.XA Expired - Fee Related CN113670879B (zh) 2021-08-27 2021-08-27 基于多变量和pmf模型对黑臭水体底泥重金属源的分析方法

Country Status (2)

Country Link
JP (1) JP7227439B1 (zh)
CN (1) CN113670879B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115859632B (zh) * 2022-12-07 2023-05-26 中科三清科技有限公司 在线源解析模型不确定度调节方法、装置与可读存储介质
CN116237353A (zh) * 2023-02-21 2023-06-09 河北地质大学 一种修复填埋物重金属污染的方法
CN117541071B (zh) * 2023-10-26 2024-08-23 生态环境部土壤与农业农村生态环境监管技术中心 一种场地土壤重金属损害基线计算方法和装置
CN117874616B (zh) * 2024-01-12 2024-07-30 广东工业大学 基于综合偏离度的污染物溯源方法及装置、电子设备
CN118135412A (zh) * 2024-01-15 2024-06-04 广东省科学院广州地理研究所 一种农村黑臭水体遥感识别方法
CN117786616B (zh) * 2024-02-23 2024-05-10 自然资源部第二海洋研究所 一种近海沉积物重金属污染特征确定方法、装置及介质
CN117911054B (zh) * 2024-03-15 2024-06-07 生态环境部华南环境科学研究所(生态环境部生态环境应急研究所) 一种黑臭水体底泥污染物溯源识别方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3131661B2 (ja) * 1992-07-07 2001-02-05 株式会社日立製作所 水質汚染監視装置及びその方法
JP3281213B2 (ja) * 1995-03-24 2002-05-13 株式会社東芝 沸騰水型原子炉プラントの水質管理方法
WO2009153855A1 (ja) 2008-06-17 2009-12-23 株式会社島津製作所 水質分析計
JP2014037677A (ja) * 2012-08-10 2014-02-27 Japan River Front Research Center ４次元水循環再現・解析・予測・可視化シミュレーションシステム
WO2014087521A1 (ja) 2012-12-06 2014-06-12 株式会社Ｗａｔｅｒ ４次元水循環再現・解析・予測・可視化シミュレーションシステム
JP2016161342A (ja) 2015-02-27 2016-09-05 共冷工業株式会社 水質検知装置
JP6994408B2 (ja) 2018-02-27 2022-01-14 前田建設工業株式会社 重金属汚染土壌の迅速分析方法
CN209470932U (zh) * 2018-11-21 2019-10-08 江苏福玉德农业科技有限公司 一种用于河道底泥的采样设备
JP7076163B2 (ja) 2021-01-26 2022-05-27 Ｇ－８ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｔｒａｄｉｎｇ 株式会社 汚染土壌の処理方法

Also Published As

Publication number Publication date
JP7227439B1 (ja) 2023-02-22
JP2023033163A (ja) 2023-03-09
CN113670879A (zh) 2021-11-19

Similar Documents

Publication Publication Date Title
CN113670879B (zh) 基于多变量和pmf模型对黑臭水体底泥重金属源的分析方法
Zeng et al. Impact of intensive land use on heavy metal concentrations and ecological risks in an urbanized river network of Shanghai
Zhang et al. Assessment of heavy metal pollution and water quality characteristics of the reservoir control reaches in the middle Han River, China
Li et al. Spatial risk assessment and sources identification of heavy metals in surface sediments from the Dongting Lake, Middle China
CN113470765A (zh) 一种土壤重金属来源解析方法
Li et al. Ecological risk assessment and sources identification of heavy metals in surface sediments of a river–reservoir system
CN114384224B (zh) 基于多同位素联合示踪的流域氮磷污染物解析方法与系统
Jiménez et al. Model of behaviour of conductivity versus pH in acid mine drainage water, based on fuzzy logic and data mining techniques
Shu et al. Levels and ecological risk of heavy metals in the surface sediments of tidal flats along the North Jiangsu coast, China
Cai et al. Distribution, risk assessment, and quantitative source apportionment of heavy metals in surface sediments from the shelf of the northern South China Sea
Levitan et al. Developing protocols for geochemical baseline studies: an example from the Coles Hill uranium deposit, Virginia, USA
Wang et al. Integrated assessment of trace elements contamination in sediments of a typical aquaculture bay in China: Ecological toxicity, sources and spatiotemporal variation
Chris et al. Assessment of some heavy metal content in sediments of a mangrove swamp, niger delta, nigeria using applicable ecological risk indices
Boroumandi et al. Evaluation of soil pollution sources using multivariate analysis combined with geostatistical methods in Zanjan Basin, Iran
Haque et al. Anthropogenic and natural contribution of potentially toxic elements in southwestern Ganges–Brahmaputra–Meghna delta, Bangladesh
Jayaraju et al. Use of benthic foraminifera as a proxy for monitoring heavy metal pollution in the Swarnamukhi estuary, southeast coast of India
Glaser et al. Designing field-based investigations of organic micropollutant fate in rivers
CN107038495A (zh) 一种地下水脆弱性评价结果的检验方法
Fulk et al. Ecological assessment of streams in the coal mining region of West Virginia using data collected by the US EPA and environmental consulting firms
Barnard et al. Chemometrics in environmental chemistry-statistical methods
Shang et al. Characteristic Model of Soil Pollution in Petroleum Exploration Area Based on Multivariate Analysis.
Ngamlana et al. Association of coal fired power plants with river water quality in South Africa
Gandhimathi et al. Spatial prediction of heavy metal pollution for soils in Coimbatore, India based on ANN and kriging model
Filella et al. From headwaters to oceans: antimony across the aquatic continuum
Anjan Kumar Prusty et al. Chemical partitioning of Cu, Pb and Zn in the soil profile of a semi arid dry woodland

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220913

CF01 Termination of patent right due to non-payment of annual fee