一种土壤重金属溯源及污染途径确定方法与流程

本发明属于农业环境技术领域，具体涉及一种农田土壤重金属溯源方法及污染途径确定方法。

背景技术：

农田土壤重金属污染严重且来源复杂，直接影响农产品质量安全与人体健康。目前对于农田污染土壤的研究工作多局限于重金属分布、赋存形态危害等方面，污染源溯源研究较少且多采用同位素示踪、空间分析及多元统计结合分析。此外，关于污染途径确认及潜在污染源确定的相关研究缺乏。

上述技术方法存在的主要问题有：(1)现有污染源溯源方法复杂，实施过程耗时长成本高；(2)现有方法难以根据已有土壤检测数据直接确定潜在污染源；(3)农田土壤污染防治措施缺乏针对性，人力资金等资源过度投入。因此，需要一种能够对土壤重金属污染溯源并能够确定污染途径的方法。

技术实现要素：

本发明利用研究区域农田及已知污染源区域农田土壤数据，明确已知污染源是否为实际污染源，若为实际污染源则明确其污染途径；此外，通过污染途径，明确潜在污染源位置。同时通过计算所有污染源的贡献率，明确主要污染源，为区域性土壤污染防治工作提供决策支撑，制定针对性方案措施。

本发明的目的是提供一种农田土壤重金属溯源及污染途径确定方法。

为了解决上述技术问题，本发明公开了一种农田土壤重金属溯源及污染途径确定方法，所述方法包括：

(1)研究区域确定及数据获取；

(1.1)选取研究区域并明确背景信息，

(1.2)污染区域面积确定，

具体公式如下：

s＝max(s研究农田，s污染区域)

其中，s为研究区域内农田面积，s研究农田为研究农田面积，s污染区域是以污染源为中心，向四周辐射1-5km的区域内农田面积，面积单位为亩；

(1.3)点位布设，

(1.4)获取点位数据

获取所述点位的土壤，并对土壤进行检测，获取监测指标对应检测结果；

(1.5)对研究农田的土壤检测结果进行主成分分析，获取载荷矩阵x，所述载荷矩阵x的行为指标，列为主成分，获取得分系数矩阵y，所述得分系数矩阵y的行为监测指标，列为主成分，根据所述载荷矩阵x可得到任意单一监测指标ii：

ii＝a1f1+a2f2+…aifi

其中，ii为单一监测指标，仅示意无取值；fi为单一监测指标对应主成分，仅示意无取值；ai指成分矩阵中单一监测指标ii对应的数值；

同时根据y得到任意fi：

fi＝m1i1+m2i2+…+miii

其中：mi指得分系数矩阵中单一监测指标ii对应的数值；

(2)污染差异度计算

(2.1)对于污染源区域内农田，所述区域内农田点位的监测指标均存在以下关系式，并可通过最小二乘法求解wi：

ei(fi)＝∑cimi(fi)

ci＝w1(fi)e1+w2(fi)e2+…+wi(fi)ei

其中，ei(fi)为污染源地块单一监测指标ii(如cd)对应fi下的成分因子，ci为监测指标对应检测结果，mi(fi)为得分系数矩阵中单一监测指标ii对应fi下的数值，wi(fi)为单一监测指标ii对应fi下的值；

(2.2)依据单一监测指标ii的ai(fi)和wi(fi)分别构建矩阵，求单一监测指标污染差异度，确定已知污染源是否为实际污染源，具体公式如下：

其中cosθi为单一监测指标污染差异度，a为ai(fi)组成的矩阵，且ai(fi)为成分矩阵中单一监测指标ii对应fi下的数值，w为wi(fi)组成的矩阵；所述差异度≤0.85认为已知污染源为研究区域实际污染源；＞0.85则认为已知污染源不是研究区域实际污染源；

(3)污染途径确定；

污染途径确定，其特征在于：

(3.1)选取污染差异度≤0.85的污染源，依据区域背景信息预设污染途径，并布设污染途径地块，所述污染途径地块内点位数量、监测指标信息与污染区域保持一致，

(3.2)潜在污染源确定

依据步骤(1)-(2)的步骤计算步骤(3.1)确定的所述污染途径地块的污染差异度，确定差异度突变地块，针对所述突变地块进行辐射性地块布设，获取辐射性地块的污染差异度，循环此操作，直至差异度稳定变化，

所述辐射性地块内点位数量、监测指标等信息与污染区域保持一致，每次循环布设一个或多个辐射性地块；

(4)污染源贡献率

依据上述已知污染源及潜在污染源，获取污染地块所有监测指标的w值，所有污染地块与研究地块存在关系如下式；联立所有监测指标关系式，通过最小二乘法求解ki，ki即为各污染源地块的贡献率：

ai(fi)＝k1wi⁽¹⁾(fi)+k2wi⁽²⁾(fi)+…+kiwi⁽ⁿ⁾(fi)

其中，wi⁽¹⁾(fi)、wi⁽²⁾(fi)、wi⁽ⁿ⁾(fi)与wi求法一致，仅用于区分不同污染源地块；

所述ki和为1。

进一步地，所述辐射性地块差异度稳定出现变化——越接近突变地块差异度越小，则最终布设的辐射性地块为潜在污染源所在位置；

进一步地，所述背景信息包括农田位置、类型、面积，已知污染源数量、类型、位置、面积、时期，研究区域水系流向、风向、道路，研究农田到各污染源间的距离；

进一步地，采取均匀布点法计算各区域内的点位数量，计算公式为：

其中，n为点位数量，计算结果为小数时向上取整；n为监测指标数量，且n≥5；

所述研究农田及所述污染区域内农田中点位布设数量一致；

进一步地，所述监测的检测指标包括但不限于cd、as、pb、cr、cu、zn、ni、mg、fe、mn、ca、si、se；

进一步地，所述预设污染途径依据区域背景信息，取研究农田与各污染源连线，并在连线上设置污染途径地块，所述污染途径地块中心点设置在所述连线上；

进一步地，每条所述污染源连线上设置的所述污染途径地块布设数量≥1；

进一步地，所述确定差异度突变地块为所述研究农田、污染途径地块及污染源的污染差异度呈单调递增或递减变化过程中的异常污染差异度对应的地块；

本发明的农田土壤重金属溯源及污染途径确定方法，具有以下优点：

1.点位数据获取时间快，土壤检测结果实时获取；点位数据获取成本低：传统判别是否为污染源，需进行同位素分析，成本较高，本发明单个点位数据获取成本仅为常规成本的10-20％；

2.依据已有土壤数据，可推断潜在污染源位置：现有的方法只能分析已知污染源，已经关停的企业或停用的灌渠等造成的重金属等污染，现有方法无法对其分析；

3.已知污染源，可同时推测污染源贡献率及污染途径，更利于制定针对性决策措施建议，减少现有大区域措施的人力物资等资源浪费；

4.根据污染差异度确定污染差异度突变地块位置，并在所述突变地块进行辐射性地块布设进一步确定潜在污染源，大大降低了工作量，为土壤重金属污染的工作提供了技术支持。

附图说明

图1为农田土壤重金属溯源及污染途径确定方法的技术流程图；

图2为某研究区域示意图；

图3为研究区域点位分布图；

图4为各地块污染差异度示意图；

图5为污染途径示意图；

图6为污染区域w值示意图；

具体实施方式

下面通过实施例对本发明做进一步的详细说明，以令本领域技术人员参照说明书文字能够据以实施。

应当理解，本文所使用的诸如“具有”，“包含”以及“包括”术语并不排除一个或多个其它元件或其组合的存在或添加。

1.研究区域确定及数据获取

①选取某地内农田m为研究对象，确定该乡镇存在两个污染源a(工矿业)和b(养殖业)；m到a距离约为8.5km，m到b距离约为3km(图2)；

②确定农田m的面积为780亩，以污染源为中心划分相同面积农田区域；

③点位布设：依据点位布设公式，各区域地块均匀布设41个点位，获取点位地理信息(图3)；

④点位数据获取：a、b、m三个区域所有点位通过x荧光快速检测仪器检测cd、as、pb、cr、cu、zn、ni、mg、fe、mn、ca、si、se等指标，现场获取检测结果。

2.污染差异度计算

①以cd为研究指标，对农田m点位信息进行主成分分析，得到载荷矩阵和得分系数矩阵；

②依据差异度计算公式，计算污染区域污染差异度，其中a、b区域差异度均小于0.85，见图4。

③依据区域水系流向、风向等信息，在农田m到污染源a、b间均预设污染途径，并根据污染途径距离布设污染途径地块，其中m到a布设3个地块(c、d、e)，m到b布设1个地块(f)；

3.污染途径确定

①依据污染差异度，可知区域d污染差异度突增(由0.014突增至0.029)，说明其与农田m差异度增强，可能受其他污染源影响，在其两边预设相反污染途径并布设地块(g、h)，布设点位并获取点位地理信息；

②通过x荧光快速检测仪器获取g、h点位检测数据，计算区域污染差异度，确定三条污染途径(ma、mb、dh)，见图5。

4.污染源贡献率确定

①依据最小二乘法，计算出污染区域a、b所有指标wi值，见图6；

②通过农田的成分矩阵中的ai值及污染区域wi值，计算得到a的贡献率为0.865，b的贡献率为0.135，确定a与b均为农田m的污染源，且a为农田m的主要污染源，判定结果与实际认知相符。

尽管本发明的实施方案已公开如上，但其并不仅仅限于说明书和实施方式中所列运用，它完全可以被适用于各种适合本发明的领域，对于熟悉本领域的人员而言，可容易地实现另外的修改，因此在不背离权利要求及等同范围所限定的一般概念下，本发明并不限于特定的细节和这里示出与描述的实施例。