评估土壤污染生态风险的生物毒性归一化方法及其应用

摘要

本发明公开了一种评估土壤污染生态风险的生物毒性归一化方法及其应用，包括：依据目标区域的生物类别和区系，选择具有区域代表性的生物物种，确定代表生物物种的生物学分类，并给出门、科和属生物分类参数；搜集或测试目标污染物对代表生物物种的毒性值；搜集或测试毒性值所涉及的土壤性质；对有多数毒性值和土壤性质相匹配的物种进行多元线性回归归一化，建立多元回归归一化模型；对有少数毒性值和土壤性质相匹配的物种则进行种间外推归一化，建立种间外推归一化模型；对多元回归归一化模型或种间外推归一化模型的预测效果进行校验。本发明有效解决了土壤污染生态风险评估中土壤生物毒性值未归一化而产生的评估结果差异较大的问题。

权利要求

1.一种评估土壤污染生态风险的生物毒性归一化方法，其特征在于，包括以下步骤：S1.依据目标区域的生物类别和区系，选择具有区域代表性的生物物种，确定所述代表生物物种的生物学分类，并给出门、科和属生物分类参数；

S2.搜集或测试目标污染物对所述代表生物物种的毒性值；

S3.搜集或测试所述毒性值所涉及的土壤性质；

S4.对有多数毒性值和土壤性质相匹配的物种进行多元线性回归归一化，建立多元回归归一化模型；对有少数毒性值和土壤性质相匹配的物种则进行种间外推归一化，建立种间外推归一化模型；

S5.对所述多元回归归一化模型或所述种间外推归一化模型的预测效果进行校验；

在S4中，采用所述土壤性质作为自变量，毒性值为因变量，将毒性值和对应土壤性质进行多元线性回归分析，建立所述多元回归归一化模型，包括：NOEC或ECx＝apH+e；

NOEC或ECx＝apH+bOC+e；

NOEC或ECx＝apH+bOC+cClay+e；

NOEC或ECx＝apH+bOC+cClay+dCEC+e；

其中，a，b，c，d和e均为模型参数；ECx为污染物对受试生物产生x％效应的浓度；

当同物种的毒性值不足时，依据建立所述多元回归归一化模型相同的方法，构建同属水平的回归模型，作为所述同属的多元回归归一化模型；

当同属的毒性值不足时，依据建立所述多元回归归一化模型相同的方法，构建同科水平的回归模型，作为所述同科的多元回归归一化模型；

S4还包括：

将多元回归归一化模型log化：

log(NOEC或ECx)＝apH+e；

log(NOEC或ECx)＝apH+blogOC+e；

log(NOEC或ECx)＝apH+blogOC+clogClay+e；

log(NOEC或ECx)＝apH+blogOC+clogClay+dlogCEC+e；

根据log化的所述多元回归归一化模型构建种间外推模型：NOECp或ECxp＝NOECs或

NOECp或ECxp＝NOECs或

NOECp或ECxp＝NOECs或

NOECp或ECxp＝NOECs或

式中：NOECs是归一化前的毒性值，NOECp是归一化后的毒性值，pHs、OCs、Clays和CECs分别为原土壤的pH、有机碳含量、土壤黏粒含量和阳离子交换量，pHp、OCp、Clayp和CECp是归一化目标的pH、有机碳含量、土壤黏粒含量和阳离子交换量，a、b、c和d为相近物种的log化归一化模型中pH、OC、Clay和CEC的模型参数；ECxp和ECxs分别是归一化后和归一化前的产生x％效应的浓度；

当同物种或属的毒性值不足时，依据上述种间外推归一化方法构建属间或科间的外推模型；

S5具体包括以下内容：将目标土壤性质代入所述多元回归归一化模型或所述种间外推归一化模型，获得预测值，计算预测值和实测值的差异系数，所述差异系数＝MAX(预测值，实测值)/MIN(预测值，实测值)，根据差异系数来划分预测效果。

2.根据权利要求1所述的一种评估土壤污染生态风险的生物毒性归一化方法，其特征在于，S2中当所述代表生物物种水平毒性值不足时，搜集其同属生物的毒性值；当同属的毒性值不足时，搜集其同科生物的毒性值，包括急性生态毒性值和慢性生态毒性值。

3.根据权利要求1所述的一种评估土壤污染生态风险的生物毒性归一化方法，其特征在于，在S3中，所述土壤性质包括但不限于：pH、有机碳OC、阳离子交换量CEC和土壤黏粒Clay。

4.根据权利要求1所述的一种评估土壤污染生态风险的生物毒性归一化方法，其特征在于，还包括对搜集到的毒性值进行筛选，具体内容为：删去未按标准方法开展试验的数据，其中包括无对照组或未设平行组；

删去对照组生物生长不符合标准的数据，其中包括对照组发芽率低于70％，平均存活率低于90％；

删去没有明确毒性终点和没有土壤性质的数据，其中包括明确毒性终点是根伸长或者是生物量、明确土壤性质；

删去水培条件下的数据，保留试验介质为土壤的数据。

5.根据权利要求1所述的一种评估土壤污染生态风险的生物毒性归一化方法，其特征在于，差异系数在2倍以内表明预测效果优秀，在2倍至5倍以内表明预测效果良好，高于5倍则预测效果一般。

6.一种评估土壤污染生态风险的生物毒性归一化方法的应用，基于权利要求1‑5任意一项一种评估土壤污染生态风险的生物毒性归一化方法，其特征在于，包括以下步骤：(1)目标土壤性质和生物毒性值的测试与获取，按照标准测试方法，获得目标土壤的土壤性质参数和生物毒性值；

(2)将目标土壤性质代入上述构建的同物种、同属或同科的回归模型，或代入上述构建的种间、属间或科间外推模型；

(3)计算获得归一化至目标土壤性质下的土壤生物毒性数据；对于种间、属间或科间外推后有多个数据的物种，取所有数据的几何平均值

说明书

评估土壤污染生态风险的生物毒性归一化方法及其应用

技术领域

[0001] 本发明涉及生态环境保护技术领域，更具体的说是涉及一种评估土壤污染生态风险的生物毒性归一化方法及其应用。

背景技术

[0002] 在评估土壤中污染物的生态风险时，一般进行污染物对本土生物的生态毒理学试验或进行生物毒性值搜集，然后依据所得到的急慢性毒性值，采用相应的方法推导这些重金属的土壤环境基准值，最后进行污染风险评估。在这个过程中，一般不会考虑土壤性质对毒性值和生态风险的影响，然而，研究表明土壤理化性质是影响土壤中污染物毒性效应的重要因素，因此，在评估土壤污染的生态风险时，急需土壤生态毒性的归一化方法和模型，为土壤生态风险评估和环境基准制定提供支持。

[0003] 因此，如何提供一种简单快速的评估土壤污染生态风险的生物毒性归一化方法是本领域技术人员亟需解决的问题。

发明内容

[0004] 有鉴于此，本发明提供了一种评估土壤污染生态风险的生物毒性归一化方法及其应用，对不同土壤生物的毒性值进行归一化，获得典型污染物的归一化模型，本发明可消除土壤性质差异所带来的毒性值的差异。

[0005] 为了实现上述目的，本发明采用如下技术方案：

[0006] 一种评估土壤污染生态风险的生物毒性归一化方法，包括以下步骤：

[0007] S1.依据目标区域的生物类别和区系，选择具有区域代表性的代表生物物种，确定所述代表生物物种的生物学分类，并给出门、科和属生物分类参数；

[0008] S2.搜集目标污染物对所述代表生物物种的毒性值；

[0009] S3.搜集或测试所述毒性值所涉及的土壤性质；

[0010] S4.对有多数毒性值和土壤性质相匹配的物种进行多元线性回归归一化，建立多元回归归一化模型；对有少数毒性值和土壤性质相匹配的物种则进行种间外推归一化，建立种间外推归一化模型；

[0011] S5.对所述多元回归归一化模型或所述种间外推归一化模型的预测效果进行校验。

[0012] 优选的，S2中当所述代表生物物种水平毒性值不足时，搜集其同属生物的毒性值；当同属的毒性值不足时，搜集其同科生物的毒性值，包括急性生态毒性值和慢性生态毒性值。

[0013] 需要说明的是：

[0014] 根据界门纲目科属种的各个阶段搜集的资料中明确标明的毒性值的数量来判断水平毒性值。一般少于1为水平毒性值不足，即不存在物种水平毒性值。

[0015] 优选的，在S3中，所述土壤性质包括：pH、有机碳OC、有机质OM、阳离子交换量CEC和土壤黏粒Clay等。

[0016] 优选的，还包括对搜集到的毒性值进行筛选，具体内容为：

[0017] 删去未按标准方法开展试验的数据，其中包括无对照组或未设平行组；

[0018] 删去对照组生物生长不符合标准的数据，其中包括对照组发芽率低于70％，平均存活率低于90％；

[0019] 删去没有明确毒性终点和没有土壤性质的数据，其中包括明确毒性终点是根伸长或者是生物量、明确土壤性质；

[0020] 删去水培条件下的数据，保留试验介质为土壤的数据。

[0021] 优选的，在S4中，采用所述土壤性质作为自变量，毒性值为因变量，将毒性值和对应土壤性质进行多元线性回归分析，建立所述多元回归归一化模型，包括：

[0022] NOEC＝apH+e；

[0023] NOEC＝apH+bOC+e；

[0024] NOEC＝apH+bOC+cClay+e；

[0025] NOEC＝apH+bOC+cClay+dCEC+e；

[0026] 其中，a，b，c，d和e均为模型参数；ECx是污染物对受试生物产生x％效应的浓度，比如EC10(x＝10时)，即为产生了10％危害效应的浓度；

[0027] 当同物种的毒性值不足时，依据建立所述多元回归归一化模型相同的方法，构建同属水平的回归模型，作为所述种间外推归一化模型；

[0028] 当同属的毒性值不足时，依据建立所述多元回归归一化模型相同的方法，构建同科水平的回归模型，作为所述种间外推归一化模型。

[0029] 优选的，S4还包括：

[0030] 将多元回归归一化模型log化：

[0031] log(NOEC或ECx)＝apH+e；

[0032] log(NOEC或ECx)＝apH+blogOC+e；

[0033] log(NOEC或ECx)＝apH+blogOC+clogClay+e；

[0034] log(NOEC或ECx)＝apH+blogOC+clogClay+dlogCEC+e；

[0035] 根据log化的所述多元回归归一化模型构建种间外推模型构建种间外推模型：

[0036]

[0037]

[0038]

[0039]

[0040] 式中：NOECs是归一化前的毒性值，NOECp是归一化后的毒性值，pHs、OCs、Clays和CECs分别为原土壤的pH、有机碳含量、土壤黏粒含量和阳离子交换量，pHp、OCp、Clayp和CECp是归一化目标的pH、有机碳含量、土壤黏粒含量和阳离子交换量，a、b、c和d为相近物种的log型归一化模型中pH、OC、Clay和CEC的模型参数；ECxp和ECxs分别是归一化后和归一化前的产生x％效应的浓度；

[0041] 当同物种或属的毒性值不足时，依据上述种间外推归一化方法构建属间或科间的外推模型。

[0042] 优选的，S5具体包括以下内容：将目标土壤性质代入所述多元回归归一化模型或所述种间外推归一化模型，获得预测值，计算预测值和实测值的差异系数，所述差异系数＝MAX(预测值，实测值)/MIN(预测值，实测值)，根据差异系数来划分预测效果。

[0043] 优选的，差异系数在2倍以内表明预测效果优秀，在2倍至5倍以内表明预测效果良好，高于5倍则预测效果一般。

[0044] 一种评估土壤污染生态风险的生物毒性归一化方法的应用，包括以下步骤：

[0045] (1)目标土壤性质和生物毒性值的测试与获取，按照标准测试方法，获得目标土壤的土壤性质参数和生物毒性值；

[0046] (2)将目标土壤性质代入上述构建的同物种、同属或同科的回归模型，或代入上述构建的种间、属间或科间外推模型；

[0047] (3)计算获得归一化至目标土壤性质下的土壤生物毒性数据；对于种间、属间或科间外推后有多个数据的物种，取所有数据的几何平均值

[0048] 经由上述的技术方案可知，与现有技术相比，本发明公开提供了一种评估土壤污染生态风险的生物毒性归一化方法，能够根据土壤性质对土壤生物毒性值进行归一化分析，获得土壤生物的多元回归和种间外推归一化模型，归一化模型可较好地减少不同土壤性质引起的毒性值差异，并可为土壤环境基准的制定提供较丰富的生物毒性数据。依据归一化模型，获得了不同土壤类型条件下铜和铅的土壤环境基准值，然后应用于后续的生态风险评估。归一化后推导出的基于不同土壤类型的生态风险评估，可避免土壤性质差异所带来的毒性效应的差异，可为土壤重金属的污染预防、风险控制和污染修复提供科学依据，并对不同土壤条件下土壤重金属的生态风险评估和环境管理提供技术支持。

附图说明

[0049] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。

[0050] 图1附图为本发明提供的一种评估土壤污染生态风险的生物毒性归一化方法的流程图。

具体实施方式

[0051] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0052] 本发明实施例公开了一种评估土壤污染生态风险的生物毒性归一化方法，在本实施例中选取铜和铅来进行具体说明，包括以下内容：

[0053] 1材料与方法

[0054] 1.1铜和铅的毒性数据的搜集和筛选

[0055] 1)毒性数据的搜集：①以主题“土壤”、“铜”、“铅”在中国知网(http://www.cnki.net)分别搜索铜和铅的生态毒性数据；②在美国ECOTOX数据库(http://cfpub.epa.gov/ecotox)中搜索铜和铅的生态毒性数据；③采用Elsevier(http://www.sciencedirect.com)等数据库对外文文献中铜和铅的生态毒性数据进行搜索，查找文献中报道了土壤性质的数据。

[0056] 2)毒性数据的筛选：删去未按标准方法的开展试验的数据，如无对照组，未设平行组；删去对照组生物生长不符合标准的数据，如对照组发芽率低于70％，平均存活率低于90％；删去没有明确毒性终点和没有土壤性质的数据，如毒性终点要明确是根伸长或者是生物量等，土壤性质(pH、OC等)在文章中要有明确标注；删去水培条件下的数据，保留试验介质为土壤的数据。采用国内分布广泛的物种以及具有代表性的标准测试物种，用于推导我国的土壤环境基准值。

[0057] 1.2数据归一化处理

[0058] 1)对于有多个毒性值和土壤性质相匹配的物种，足够自身建立回归模型，进行多元线性归一化。土壤是一类高度不均匀的介质，不同的pH、OC(有机碳含量)、CEC(阳离子交换量)和Clay(黏土含量)会造成土壤中毒性物质被生物吸收的含量的差异，利用SPSS软件2
(SPSS26)回归分析的R 和显著性大小判断多元回归的效果好坏。将毒性数值和对应土壤性质进行多元线性回归分析，来建立回归模型(例如NOEC＝apH+bOC+cCEC+dClay+e；a，b，c，d，e为模型参数)，最后将毒性值归一化至中性土壤条件下。

[0059] 2)对于有少数毒性值和土壤性质相匹配的物种，不足以自身建立回归模型，则进行种间外推归一化。因为这些物种数据量较少，无法自身建立多元回归模型，采用种间外推回归模型进行归一化。种间外推归一化以生物分类学相似的物种建立的回归模型为基础，构建种间外推模型 式中：NOECs是归一化前的毒性值，NOECp是归一化后的毒性值，pHs、OCs、Clays和CECs分别为原土壤的pH、有机碳含量、土壤黏粒含量和阳离子交换量，pHp、OCp、Clayp和CECp是归一化目标的pH、有机碳含量、土壤黏粒含量和阳离子交换量，a、b、c和d为相近物种的log型归一化模型中pH、OC、Clay和CEC的模型参数；归一化后有两个以上毒性值的物种再取几何平均值。

[0060] 3)采用上述构建的归一化模型，通过土壤性质对铜和铅的生态毒性数据进行预测，并与实测值(该实测值未参与归一化模型的构建)进行比较，验证归一化模型的预测效果。

[0061] 4)多元线性和种间外推归一化以中性土壤条件(土壤pH＝7.0，CEC＝15cmol·kg‑1，OC＝1.5％，Clay＝35％)进行计算，并与未归一化数据进行比较。

[0062] 3结果与讨论

[0063] 3.1铜的归一化模型

[0064] 1)对于有多个毒性值和土壤性质相匹配的物种，进行了多元线性回归分析(表1)。(需要说明的是：在本领域中，一般2个及以下是少数，需要利用种间外推归一化。3个毒性值和土壤性质相匹配，大多数情况下可以建立多元回归归一化模型。下同。)分别构建了10个物种或微生物过程的多元回归归一化模型，其中慢性毒性终点归一化模型8个，包括赤子爱胜蚓Eisenia fetida、白符跳虫Folsomia candida和线蚓Enchytraeus crypticus3个物种的NOEC归一化模型，西红柿Solanum lycopersicum、大麦Hordeum vulgare、小白菜Brassica chinensis、青海弧菌Vibrio qinghaiensis(Q67)4个物种的EC10归一化模型，微生物群落Protista的LOEC归一化模型，线蚓Enchytraeus albidus的LC50归一化模型和线蚓Enchytraeus luxuriosus的EC50模型。与铜的其他归一化研究相比，本研究建立了多种生物分类物种的多元回归模型。除青海弧菌外，各模型p值小于0.05，表明各个物种的多元回归归一化模型可较好的反应铜对10种土壤生物或微生物过程的毒性效应与土壤性质的变化规律。

[0065] 2)对于毒性数据不足或与土壤性质相匹配数据不足的物种，进行了种间外推归一化分析(表1)。依据种间外推归一化方法，采用pH和logOC作为模型参数归一化其他物种的毒性值时，各物种的固有敏感性会显著降低，有学者研究发现pH和OC是控制铜的土壤生物毒性的两个重要因子该研究建立了基于土壤pH和OC的赤子爱胜蚓E.fetida的归一化模型(logNOEC＝0.204pH+0.933logOC+0.709)和大麦H.vulgare的归一化模型(logEC10＝0.144pH+0.753logOC+1.036)，结合获得的实测毒性数据，构建了12种土壤生物的种间外推归一化模型(表1)，可对12种相同生物分类学的土壤生物毒性进行预测。

[0066] 3)采用表1中的代表性动植物，如赤子爱胜蚓E.fetida的归一化模型，预测实际土‑1壤性质(pH＝5.0，CEC＝7.9cmol·kg ，OC＝2.1％，Clay＝7.9％)下铜对赤子爱胜蚓的毒性‑1 ‑1
值，得到NOEC预测值为140.07mg·kg ，与实测值87.5mg·kg 相近。此外，对大麦H.vulgare的归一化模型预测效果进行验证，预测的EC10值(土壤性质pH＝7.35，OC＝‑1 ‑1 ‑1
1.25％，CEC＝8.43cmol·kg )为139.38mg·kg ，与实测值96.8mg·kg 非常接近，说明构建的铜归一化模型的预测效果较好。

[0067] 表1铜的归一化模型

[0068]

[0069] 注：标a的是用赤子爱胜蚓的log模型的系数进行的种间外推归一化，标b的是用大麦log模型的系数进行的种间外推归一化。

[0070] 3.2铅的归一化模型

[0071] 1)对于有多个毒性值和土壤性质相匹配的物种，进行了多元线性回归分析(表2)。分别构建了5个物种的多元回归归一化模型，其中慢性毒性终点归一化模型4个，包括赤子爱胜蚓E.fetida的NOEC归一化模型，赤子爱胜蚓E.fetida、西红柿S.lycopersicum、大麦H.vulgare3个物种的EC10归一化模型，和线蚓E.albidus的LC50模型。模型p值小于0.05，反应了各个物种的多元回归归一化模型可较好的反应铅对5种土壤生物毒性效应与土壤性质的变化规律，与前人相比，该研究建立的多元回归模型物种分类较全面，采用的模型参数较丰富。其中赤子爱胜蚓E.fetida的NOEC模型比EC10模型更显著，在后续归一化取值时采用赤子爱胜蚓E.fetida的NOEC模型。

[0072] 2)建立种间外推归一化模型可增加土壤环境基准值推导中的物种数量，使得推导出的基准值更加科学可靠，从而能够保护更多的物种。对于毒性数据不足或与土壤性质相匹配数据不足的物种，进行了种间外推归一化分析(表2)。依据种间外推归一化方法，结合铅的土壤生物毒性与pH和OC的重要相关性，该研究建立了基于土壤pH和OC的赤子爱胜蚓E.fetida的归一化模型(logNOEC＝0.662pH+3.448logOC‑2.741)和大麦H.vulgare的归一化模型(logEC10＝0.116pH+0.311logOC+1.74)，结合获得的实测毒性数据，构建了12种土壤生物的种间外推归一化模型(表2)，可对12种相近生物分类学的土壤生物毒性进行预测。

[0073] 3)采用表2中，大麦H.vulgare的归一化模型预测实际土壤性质下(pH＝8.86，OC＝‑11.02％，CEC＝8.0cmol·kg ，Clay＝19.6％)铅对大麦的毒性值，得到的EC10预测值为‑1 ‑1
547.67mg·kg ，与实测值619mg·kg 非常接近，说明构建的铅归一化模型的预测效果较好。

[0074] 表2铅的归一化模型

[0075]

[0076] 3.3铜和铅在中性土壤条件下的归一化毒性数据

[0077] 1)根据构建的归一化模型(表1)，把各土壤条件下铜的生态毒性数据归一化至中‑1性土壤条件下(土壤pH＝7.0，CEC＝15cmol·kg ，OC＝1.5％，Clay＝35％)，归一得到5种土壤动物的毒性值，13种土壤植物的毒性值，共计4门11科18种土壤动(植)物(表3)，还获得了
2种微生物过程的归一化毒性数据(表3)。

[0078] 表3铜的归一化毒理数据

[0079]

[0080] 2)根据构建的归一化模型(表2)，把各土壤条件下铜的生态毒性数据同样归一化至中性土壤条件下，归一得到5种土壤动物的毒性值，10种土壤植物的毒性值，共计5门10科15种土壤动(植)物(表4)。

[0081] 表4铅的归一化毒理数据

[0082]

[0083] 4结论

[0084] 根据土壤性质(pH、OC、CEC和Clay)对土壤生物毒性值进行归一化分析，获得铜和铅的22和17种土壤生物的多元回归和种间外推归一化模型，通过验证发现归一化模型可较好地预测其它土壤条件下的毒性值，从而减少不同土壤性质引起的毒性值差异。

[0085] 本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。

[0086] 对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，

[0087] 而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。