CN114062476A - 一种土壤Cd/Pb复合污染双金属同位素源解析方法及系统 - Google Patents

Abstract

本发明公开了一种农田土壤Cd/Pb复合污染双金属同位素源解析方法，包括以下步骤：通过样品采集装置分别采集土壤样品和风险源样品；分别测定土壤样品和风险源样品的Cd、Pb同位素比值；利用土壤样品和风险源样品的Cd、Pb同位素比值为坐标进行作图，得到同位素比值投影图；通过同位素比值投影图对污染农田土壤的污染端元进行识别，得到污染端元识别结果，进而进行污染端元确认；计算污染端元的相对贡献率，得到同位素源解析结果；本发明利用Cd‑Pb两种金属同位素相互制约、相互印证，解决由于样品同位素值相似或重叠以及同位素分馏造成的解析准确性差的问题，实现对农田土壤Cd/Pb重金属污染源的精准识别和定量解析。

Description

一种土壤Cd/Pb复合污染双金属同位素源解析方法及系统

技术领域

本发明涉及土壤重金属污染防控的研究领域，特别涉及一种土壤Cd/Pb复合污染双金属同位素源解析方法及系统。

背景技术

重金属Cd和Pb是全球公认的毒性金属元素。土壤中的Cd和Pb较难迁移，具有毒性大、难降解、能通过食物链和食物网生物放大等特点，对人类健康有巨大威胁。因此，土壤Cd/Pb复合污染具有一定的特殊性，是土壤重金属污染治理的难点，同时也一直都是国内外瞩目的热点和难点研究课题。随着农田土壤Cd/Pb污染程度日益加剧，单纯研究土壤重金属污染物形态、种类及空间分布已无法满足现有的农田土壤治理需求。但是由于农田土壤介质环境十分复杂，污染来源较难精确识别、污染贡献难以定量，导致农田土壤重金属污染难以针对性进行治理，治理效果不尽人意。因此研发能够有效厘定农田土壤中Cd/Pb污染来源及定量各污染源贡献率的方法，在土壤Cd/Pb污染治理行动中显得十分迫切且必要。

目前绝大部分土壤中Cd、Pb的污染来源解析工作通常依托于大量的数据库和数学统计分析。比如因子分析法、主成分分析法、聚类分析法、富集因子法等方法，但是这些方法只能实现对土壤中重金属的来源类型进行定性分析；另外，化学质量平衡法、正定矩阵因子分解法等方法能够实现定性定量层面的源解析分析，但是这些方法都是建立在大量、全面的样品采集以及繁琐的数学分析上的，工作量较大且较难对多元体系进行辨析。

近些年来，随着化学分析技术的不断发展，应用金属同位素指纹特征溯源为土壤重金属源解析方法带来新的思路。自然界中每种物质都有自己独特的同位素组成“标签”，因此可以通过不同物质同位素组成这一特定的“标签”来区分混合物质的来源。Cd有8种同位素：¹⁰⁶Cd、¹⁰⁸Cd、¹¹⁰Cd、¹¹¹Cd、¹¹²Cd、¹¹³Cd、¹¹⁴Cd以及¹¹⁶Cd；Pb在自然界中以四种同位素的形式存在:²⁰⁴Pb、²⁰⁶Pb、²⁰⁷Pb和²⁰⁸Pb，其中²⁰⁴Pb是大爆炸中唯一形成的原始稳定同位素，而²⁰⁶Pb、²⁰⁷Pb和²⁰⁸Pb分别是²³⁸U、²³⁵U和²³²Th的放射性衰变产物。目前已有少量利用Cd或Pb同位素中的一种进行单一金属同位素土壤污染溯源的报道，通常采用金属的浓度或其倒数与金属稳定同位素比值分析污染源。Cloquet等采集了某冶炼厂周围的土壤、粉尘颗粒、锅炉中残余熔渣等可能的污染源，测定了Cd同位素，通过Cd同位素值与Cd浓度的倒数作图初步判断污染源头，通过线性分析得出粉尘颗粒以及废弃熔渣是该地区土壤的Cd污染的主要来源。Liu等在土壤-水稻系统中采用田间监测、Pb同位素比值分析等手段对土壤Pb来源进行解析，结果表明稻田土壤中Pb来源为背景土壤、化肥、大气沉降和灌溉水。

但是土壤中的吸附、溶解、氧化还原反应和生物过程等地球化学过程会导致Cd同位素发生分馏，这会使污染源Cd同位素信号模糊，从而降低Cd单一同位素源解析结果的准确性。

发明内容

本发明的主要目的在于克服现有技术的缺点与不足，提供一种土壤Cd/Pb复合污染双金属同位素源解析方法及系统。

本发明的第一目的在于提供一种土壤Cd/Pb复合污染双金属同位素源解析方法。

本发明的第二目的在于提供一种土壤Cd/Pb复合污染双金属同位素源解析系统。

本发明的第一目的通过以下的技术方案实现：

一种土壤Cd/Pb复合污染双金属同位素源解析方法，包括以下步骤：

通过样品采集装置分别采集土壤样品和风险源样品，分别获得土壤样品和风险源样品；

分别测定土壤样品和风险源样品的Cd、Pb同位素比值，获得土壤样品的Cd、Pb同位素比值和风险源样品的Cd、Pb同位素比值；

利用土壤样品的Cd、Pb同位素比值为坐标进行作图，利用风险源样品的Cd、Pb同位素比值为坐标进行作图，得到同位素比值投影图；

通过同位素比值投影图对污染农田土壤的污染端元进行识别，得到污染端元识别结果，进而进行污染端元确认；

计算污染端元的相对贡献率，得到同位素源解析结果。

进一步地，所述采集土壤样品，具体为：通过样品采集装置采集土壤样品，获得不同距离的土壤样品；所述采集风险源样品，具体为：采集不同类型的风险源样品，包括第一类型风险源样品、第二类型风险源样品、第三类型风险源样品、第四类型风险源样品、第五类型风险源样品。

进一步地，所述不同距离的土壤样品为不同距离的农田中0-20厘米表层土壤。

进一步地，所述采集土壤样品，具体为：通过样品采集装置采集土壤样品，获得不同方位的土壤样品；所述采集风险源样品，具体为：采集不同类型的风险源样品，包括第一类型风险源样品、第二类型风险源样品、第三类型风险源样品、第四类型风险源样品、第五类型风险源样品。

进一步地，所述Cd同位素比值表达为δ^114/110Cd，所述Pb同位素比值为²⁰⁸Pb/²⁰⁶Pb和²⁰⁶Pb/²⁰⁷Pb。

进一步地，所述测定土壤样品和风险源样品的Cd、Pb同位素比值，具体如下：

Cd稳定同位素比值测定：将2.8mL AG MP-1M(100-200目)树脂装入分离柱中，首先用10mL 3.5N HNO3、2N HCL+8N HF以及6N HCl清洗树脂；加超纯水调节至树脂呈中性，后分别用10mL 2N HCl去除样品中的基质元素，10mL 1N HCl去除Mo(钼)，20mL 0.3N HCl去除Pb(铅)，20mL 0.06N HCl去除Zn(锌)，10mL 0.012N HCl去除Sn(锡)，最后用20mL 0.0012NHCl洗脱Cd并收集；收集的纯Cd溶液蒸干后溶解于3％HNO₃中待测，采用多接收等离子体质谱仪(Neptune Plus MC-ICP-MS)完成测试，采用双稀释剂法矫正质量歧视；所测得Cd稳定同位素比值结果表达为：

δ^114/110Cd＝[(¹¹⁴Cd/¹¹⁰Cd)_sample/(¹¹⁴Cd/¹¹⁰Cd)_{NIST 3108}-1]×1000，

其中，(¹¹⁴Cd/¹¹⁰Cd)_sample是所测定样品的¹¹⁴Cd/¹¹⁰Cd值，(¹¹⁴Cd/¹¹⁰Cd)_{NIST 3108}是标准NIST 3108的¹¹⁴Cd/¹¹⁰Cd值；

Pb稳定同位素比值测定：将1.5mL AG1-X8(100-200目)树脂装入分离柱中，首先用6N HCl以及MQ交替清洗树脂三次；后分别用1.5mL 1N HBr、1.5mL2N HCl去除样品中的杂质元素，最后用1.5mL 6N HCl洗脱Pb并收集；收集的纯Pb溶液蒸干后溶解于3％HNO₃中待测，采用多接收等离子体质谱仪(Neptune Plus MC-ICP-MS)完成测试，采用已标定²⁰⁵Tl/²⁰³Tl＝2.3871的标准Tl997作为内标进行仪器质量歧视校正。

进一步地，所述利用土壤样品的Cd、Pb同位素比值为坐标进行作图，利用风险源样品的Cd、Pb同位素比值为坐标进行作图，得到同位素比值投影图，具体为：利用土壤样品的Cd同位素比值为横坐标，土壤样品的Pb同位素比值作为纵坐标进行作图；利用风险源样品的Cd同位素比值为横坐标，风险源样品的Pb同位素比值作为纵坐标进行作图；得到同位素比值投影图。

进一步地，所述通过同位素比值投影图对污染农田土壤的污染端元进行识别，得到污染端元识别结果，进而进行污染端元确认，具体为：污染农田土壤的同位素比值投影点在各污染端元同位素比值投影点包围范围之内，因此与污染农田土壤的同位素比值投影点接近且将污染农田土壤的同位素比值投影点包围住的各风险源可认定为污染端元。

进一步地，所述计算污染端元的相对贡献率，得到同位素源解析结果，具体为；通过源解析计算公式计算出不同污染端元对污染农田土壤中Cd、Pb的相对贡献率。

进一步地，所述源解析计算公式，具体如下：

δ^114/110Cd_soil＝δ^114/110Cd_Ax_A+δ^114/110Cd_Bx_B+δ^114/110Cd_Cx_C+δ^114/110Cd_Dx_D，

1＝x_A+x_B+x_C+x_D，

其中，A，B，CD分别代表4个污染端元；δ^114/110Cd为Cd同位素比值；δ^114/110Cd_A，δ^{114/ 110}Cd_B，δ^114/110Cd_C，δ^114/110Cd_D分别代表A，B，C，D 4个污染端元的Cd同位素比值；

代表污染土壤的Pb同位素比值，

分别代表A，B，C，D 4个污染端元的Pb同位素比值，当m＝208时，n＝206或者当m＝206时，n＝207，x_A，x_B，x_C，x_D分别代表A，B，C，D四个污染端元的贡献率。

本发明的第二目的通过以下技术方案实现：

一种农田土壤Cd/Pb复合污染双金属同位素源解析系统，包括：

样品采集模块，用于采集土壤样品和风险源样品；

样品测定模块，用于测定土壤样品和风险源样品的Cd、Pb同位素比值，获得土壤样品的Cd、Pb同位素比值和风险源样品的Cd、Pb同位素比值；

作图模块，利用土壤样品的Cd、Pb同位素比值为坐标进行作图，利用风险源样品的Cd、Pb同位素比值为坐标进行作图，得到同位素比值投影图；

识别及确认模块，用于同位素比值投影图对污染农田土壤的污染端元进行识别，得到污染端元识别结果，进而进行污染端元确认；

相对贡献率计算模块，计算污染端元的相对贡献率，得到同位素源解析结果；

结果输出模块，用于输出最终同位素源解析结果。

本发明与现有技术相比，具有如下优点和有益效果：

1、本发明依托农田土壤以及污染源样品中Cd与Pb元素同位素比值图，可以对农田土壤中的Cd/Pb污染端元进行精确识别。

2、本发明依托农田土壤以及污染源样品中Cd与Pb元素的同位素特征，可以精确的对农田土壤中的Cd/Pb污染来源进行定量解析，确定不同污染源对土壤Cd/Pb污染的定量贡献率。

3、本发明通过测定土壤以及污染源中的Cd/Pb同位素，开发出了依托双同位素指纹技术对土壤重金属污染来源进行追溯的方法，可以利用两种金属同位素相互制约、相互印证。与传统多元统计、单一同位素指纹等方法相比，能精确确定不同污染源的贡献率，解析结果更为客观准确，对污染源的追溯效果更好。

附图说明

图1是本发明所述一种农田土壤Cd/Pb复合污染双金属同位素源解析方法流程图；

图2是本发明所述实施例1中农田土壤污染端元识别图；

图3是本发明所述实施例2中农田土壤污染端元识别图；

图4是本发明所述一种农田土壤Cd/Pb复合污染双金属同位素源解析系统结构框图。

具体实施方式

下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。

实施例1：

一种农田土壤Cd/Pb复合污染双金属同位素源解析方法，如图1所示，包括以下步骤：

通过样品采集装置分别采集土壤样品和风险源样品，分别获得土壤样品和风险源样品；

分别测定土壤样品和风险源样品的Cd、Pb同位素比值，获得土壤样品的Cd、Pb同位素比值和风险源样品的Cd、Pb同位素比值；

利用土壤样品的Cd、Pb同位素比值为坐标进行作图，利用风险源样品的Cd、Pb同位素比值为坐标进行作图，得到同位素比值投影图；

通过同位素比值投影图对污染农田土壤的污染端元进行识别，得到污染端元识别结果，进而进行污染端元确认；

计算污染端元的相对贡献率，得到同位素源解析结果。

应用Cd/Pb双金属同位素对农田土壤Cd/Pb复合污染进行源解析，具体如下：

选择贵州喀斯特地区受某大型铅锌矿场影响的农田土壤作为实施区域，根据实施区域污染特征，在距离矿山200米、1000米、2000米位置的水稻田中采集0-20厘米表层土壤，分别编号R1、R2、R3，同时采集潜在风险源样品：尾矿、矿山粉尘、母质、化肥、背景土壤。即包含第一距离土壤样品、第二距离土壤样品、第三距离土壤样品，所述第一距离土壤样品为距离200米位置的农田中0-20厘米表层土壤，所述第二距离土壤样品为距离1000米位置的农田中0-20厘米表层土壤，所述第三距离土壤样品为距离2000米位置的农田中0-20厘米表层土壤；所述第一类型风险源样品为尾矿，所述第二类型风险源样品为矿山粉尘，所述第三类型风险源样品为母质，所述第四类型风险源样品为化肥，所述第五类型风险源样品为背景土壤。

将污染农田土壤和污染风险源样品分别进行Cd、Pb元素的分离纯化处理后，作为优选，采用多接收电感耦合等离子体质谱仪测定Cd、Pb同位素比值，测试结果如表1所示。

表1实施例1污染土壤与潜在风险源样品Cd、Pb同位素比值

本实施例采用农田土壤样品以及潜在风险源样品中δ^114/110Cd为横坐标，²⁰⁸Pb/²⁰⁶Pb为纵坐标作图，如图2所示。从图中可以明显看出，农田土壤的重金属污染源为矿山粉尘、背景土壤、尾矿，这三个污染端元中，来自矿山的污染源包括尾矿和矿山粉尘的污染影响相对更明显，而化肥以及母质的贡献可以排除。表明在该目标研究区域农田土壤，受矿山开采的影响十分显著，而农业活动污染来源的影响几乎被掩盖。

上述污染端元识别结果表明，农田土壤的污染端元有三个，分别为：矿山粉尘(A)、尾矿(B)、背景土壤(C)。将污染农田土壤(Rn，n＝1，2，3)和污染端元中的δ^114/110Cd和²⁰⁸Pb/²⁰⁶Pb导入污染贡献计算公式：

δ^114/110Cd_Rn＝δ^114/110Cd_Ax_A+δ^114/110Cd_Bx_B+δ^114/110Cd_Cx_C，

1＝x_A+x_B+x_C，

得到的计算结果如表2所示，土壤R1和R2，距离矿山较近，矿山粉尘分别贡献了79.06％和54.69％。土壤R3距离矿山较远，主要污染来自背景土壤，贡献率为53.36％。

实施例2

本实施例除下述特征外其他结构同实施例1：

选择了华南地区受某大型多金属矿场影响的农田土壤作为实施区域，根据实施区域污染特征，以矿场为中心，随机采集周边4块水稻田中0-20厘米的表层土壤，分别编号P1、P2、P3、P4，同时采集潜在风险源样品：矿山废水沉淀、大气沉降、母质、化肥、背景土壤。即周边4块水稻田包括：第一方位、第二方位、第三方位、第四方位；所述第一类型风险源样品为矿山废水沉淀，所述第二类型风险源样品为大气沉降，所述第三类型风险源样品为母质，所述第四类型风险源样品为化肥，所述第五类型风险源样品为背景土壤。

将污染农田土壤和污染风险源样品分别进行Cd、Pb元素的分离纯化处理后，作为优选，采用多接收电感耦合等离子体质谱仪测定Cd、Pb同位素比值，测试结果如表3所示。

表3实施例2污染土壤与潜在风险源样品Cd、Pb同位素比值

本实施例采用农田土壤样品以及潜在风险源样品中δ^114/110Cd为横坐标，²⁰⁶Pb/²⁰⁷Pb为纵坐标作图，如图3所示。从图中可以明显看出，污染农田土壤的重金属污染源为矿山废水沉淀、大气沉降、母质、背景土壤这四个端元。相比来自矿冶活动的影响，农业活动中化肥的污染贡献可以排除。

上述污染端元识别结果表明，农田土壤的污染端元有四个，分别为：矿山废水沉淀(A)、大气沉降(B)、母质(C)、背景土壤(D)。将污染农田土壤(Pn，n＝1，2，3，4)和污染端元中的δ^114/110Cd和²⁰⁶Pb/²⁰⁷Pb导入污染贡献计算公式：

δ^114/110Cd_Pn＝δ^114/110Cd_Ax_A+δ^114/110Cd_Bx_B+δ^114/110Cd_cx_C+δ^114/110Cd_Dx_D，

1＝x_A+x_B+x_C+x_D，

得到的计算结果如表4所示，在重金属活动性较强的华南红壤地区，人为活动重金属污染贡献占主导地位，矿山废水沉淀和大气沉降的贡献对4个农田土壤的污染贡献率均在97％以上，而自然源包括母质和背景土壤的贡献较低。

表4实施例2污染端元贡献量(％)

如上所述，根据本发明的技术方案，采用Cd/Pb同位素并联进行土壤Cd/Pb复合污染来源解析，能够精确识别不同区域农田土壤的污染端元并计算相对贡献率。本发明的两个具体实施案例分别在两个地质背景、重金属活动性差异大的西南喀斯特地区和华南红壤地区进行实施。根据Cd-Pb同位素投影图识别农田土壤的污染端元，并计算出了各个污染端元的相对贡献率。上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制。

实施例3

一种农田土壤Cd/Pb复合污染双金属同位素源解析系统，如图4所示，包括：

样品采集模块，用于采集土壤样品和风险源样品；

样品测定模块，用于测定土壤样品和风险源样品的Cd、Pb同位素比值，获得土壤样品的Cd、Pb同位素比值和风险源样品的Cd、Pb同位素比值；

作图模块，利用土壤样品的Cd、Pb同位素比值为坐标进行作图，利用风险源样品的Cd、Pb同位素比值为坐标进行作图，得到同位素比值投影图；

识别及确认模块，用于同位素比值投影图对污染农田土壤的污染端元进行识别，得到污染端元识别结果，进而进行污染端元确认；

相对贡献率计算模块，计算污染端元的相对贡献率，得到同位素源解析结果；

结果输出模块，用于输出最终同位素源解析结果。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (10)

1.一种土壤Cd/Pb复合污染双金属同位素源解析方法，其特征在于，包括以下步骤：

通过样品采集装置分别采集土壤样品和风险源样品，分别获得土壤样品和风险源样品；

分别测定土壤样品和风险源样品的Cd、Pb同位素比值，获得土壤样品的Cd、Pb同位素比值和风险源样品的Cd、Pb同位素比值；

利用土壤样品的Cd、Pb同位素比值为坐标进行作图，利用风险源样品的Cd、Pb同位素比值为坐标进行作图，得到同位素比值投影图；

通过同位素比值投影图对污染农田土壤的污染端元进行识别，得到污染端元识别结果，进而进行污染端元确认；

计算污染端元的相对贡献率，得到同位素源解析结果。

2.根据权利要求1所述的一种土壤Cd/Pb复合污染双金属同位素源解析方法，其特征在于，所述采集土壤样品，具体为：通过样品采集装置采集土壤样品，获得不同距离的土壤样品；所述采集风险源样品，具体为：采集不同类型的风险源样品，包括第一类型风险源样品、第二类型风险源样品、第三类型风险源样品、第四类型风险源样品、第五类型风险源样品。

3.根据权利要求1所述的一种土壤Cd/Pb复合污染双金属同位素源解析方法，其特征在于，所述采集土壤样品，具体为：通过样品采集装置采集土壤样品，获得不同方位的土壤样品；所述采集风险源样品，具体为：采集不同类型的风险源样品，包括第一类型风险源样品、第二类型风险源样品、第三类型风险源样品、第四类型风险源样品、第五类型风险源样品。

4.根据权利要求1所述的一种土壤Cd/Pb复合污染双金属同位素源解析方法，其特征在于，所述Cd同位素比值表达为δ^114/110Cd，所述Pb同位素比值为²⁰⁸Pb/²⁰⁶Pb和²⁰⁶Pb/²⁰⁷Pb。

5.根据权利要求4所述的一种土壤Cd/Pb复合污染双金属同位素源解析方法，其特征在于，所述测定土壤样品和风险源样品的Cd、Pb同位素比值，具体如下：

Cd稳定同位素比值测定：将2.8mL AG MP-1M树脂装入分离柱中，首先用10mL 3.5NHNO3、2N HCL+8N HF以及6N HCl清洗树脂；加超纯水调节至树脂呈中性，后分别用10mL 2NHCl去除样品中的基质元素，10mL 1N HCl去除Mo，20mL 0.3N HCl去除Pb，20mL 0.06N HCl去除Zn，10mL 0.012N HCl去除Sn，最后用20mL 0.0012N HCl洗脱Cd并收集；收集的纯Cd溶液蒸干后溶解于3％HNO₃中待测，采用多接收等离子体质谱仪完成测试，采用双稀释剂法矫正质量歧视；所测得Cd稳定同位素比值结果表达为：

δ^114/110Cd＝[(¹¹⁴Cd/¹¹⁰Cd)_sample/(¹¹⁴Cd/¹¹⁰Cd)_{NIST 3108}-1]×1000，

其中，(¹¹⁴Cd/¹¹⁰Cd)_sample是所测定样品的¹¹⁴Cd/¹¹⁰Cd值，(¹¹⁴Cd/¹¹⁰Cd)_{NIST 3108}是标准NIST3108的¹¹⁴Cd/¹¹⁰Cd值；

Pb稳定同位素比值测定：将1.5mL AG1-X8树脂装入分离柱中，首先用6N HCl以及MQ交替清洗树脂三次；后分别用1.5mL 1N HBr、1.5mL 2N HCl去除样品中的杂质元素，最后用1.5mL 6N HCl洗脱Pb并收集；收集的纯Pb溶液蒸干后溶解于3％HNO₃中待测，采用多接收等离子体质谱仪完成测试，采用已标定²⁰⁵Tl/²⁰³Tl＝2.3871的标准Tl 997作为内标进行仪器质量歧视校正。

6.根据权利要求1所述的一种土壤Cd/Pb复合污染双金属同位素源解析方法，其特征在于，所述利用土壤样品的Cd、Pb同位素比值为坐标进行作图，利用风险源样品的Cd、Pb同位素比值为坐标进行作图，得到同位素比值投影图，具体为：利用土壤样品的Cd同位素比值为横坐标，土壤样品的Pb同位素比值作为纵坐标进行作图；利用风险源样品的Cd同位素比值为横坐标，风险源样品的Pb同位素比值作为纵坐标进行作图；得到一个同位素比值投影图，其横纵坐标是一致的，横坐标为Cd同位素，纵坐标为Pb同位素。

7.根据权利要求1所述的一种土壤Cd/Pb复合污染双金属同位素源解析方法，其特征在于，所述通过同位素比值投影图对污染农田土壤的污染端元进行识别，得到污染端元识别结果，进而进行污染端元确认，具体为：污染农田土壤的同位素比值投影点在各污染端元同位素比值投影点包围范围之内，因此与污染农田土壤的同位素比值投影点接近且将污染农田土壤的同位素比值投影点包围住的各风险源可认定为污染端元。

8.根据权利要求1所述的一种土壤Cd/Pb复合污染双金属同位素源解析方法，其特征在于，所述计算污染端元的相对贡献率，得到同位素源解析结果，具体为；通过源解析计算公式计算出不同污染端元对污染农田土壤中Cd、Pb的相对贡献率。

9.根据权利要求8所述的一种土壤Cd/Pb复合污染双金属同位素源解析方法，其特征在于，所述源解析计算公式，具体如下：

δ^114/110Cd_soil＝δ^114/110Cd_Ax_A+δ^114/110Cd_Bx_B+δ^114/110Cd_Cx_C+δ^114/110Cd_Dx_D，

1＝x_A+x_B+x_C+x_D，

其中，A,B,C,D分别代表4个污染端元；δ^114/110Cd为Cd同位素比值；δ^114/110Cd_A，δ^114/110Cd_B，δ^114/110Cd_C，δ^114/110Cd_D分别代表A,B,C,D 4个污染端元的Cd同位素比值；

代表污染土壤的Pb同位素比值，

分别代表A,B,C,D 4个污染端元的Pb同位素比值，当m＝208时，n＝206或者当m＝206时，n＝207，x_A，x_B,x_C,x_D分别代表A,B,C,D四个污染端元的贡献率。

10.一种土壤Cd/Pb复合污染双金属同位素源解析系统，其特征在于，包括：

样品采集模块，用于采集土壤样品和风险源样品；

样品测定模块，用于测定土壤样品和风险源样品的Cd、Pb同位素比值，获得土壤样品的Cd、Pb同位素比值和风险源样品的Cd、Pb同位素比值；

作图模块，利用土壤样品的Cd、Pb同位素比值为坐标进行作图，利用风险源样品的Cd、Pb同位素比值为坐标进行作图，得到同位素比值投影图；

识别及确认模块，用于同位素比值投影图对污染农田土壤的污染端元进行识别，得到污染端元识别结果，进而进行污染端元确认；

相对贡献率计算模块，计算污染端元的相对贡献率，得到同位素源解析结果；

结果输出模块，用于输出最终同位素源解析结果。

Images