一种测定离子型稀土矿先锋植物对污染物毒性阈值的方法技术

本发明专利技术属于生态环境保护技术领域，本发明专利技术提供了一种测定离子型稀土矿先锋植物对污染物毒性阈值的方法。本发明专利技术方法主要包括测定离子型稀土矿污染物对先锋植物的毒性效应、筛选与污染物显著相关的指标、为各指标选取最优的剂量效应模型进行拟合、依据拟合函数计算每个评价终点指标的毒性阈值，并利用主成分分析法计算每个评价终点指标的权重、依据评价终点指标的毒性阈值与权重计算植物对污染物的综合毒性阈值。本发明专利技术对离子型稀土矿生态恢复中判断先锋植物对污染物的毒性阈值及稀土矿区植被恢复中先锋植物的筛选提供了参考价值。被恢复中先锋植物的筛选提供了参考价值。被恢复中先锋植物的筛选提供了参考价值。

【技术实现步骤摘要】
一种测定离子型稀土矿先锋植物对污染物毒性阈值的方法


[0001]本专利技术涉及生态环境保护
，尤其涉及一种测定离子型稀土矿先锋植物对污染物毒性阈值的方法。

技术介绍

[0002]离子型稀土资源是世界重稀土的重要来源，具有储量大、分布广、中重稀土含量高、配比好、采冶性能优和放射性活度低的特点。随着科学技术的快速发展，对稀土元素的需求急剧增加，稀土元素的开采工艺已经造成了许多环境问题。
[0003]稀土堆场污染是大量的残留化学药剂等有毒有害物质通过不同途径进入周围环境并累积的过程。由于早期采用的池浸、堆浸工艺较为落后，稀土冶炼过程中使用的大量化学试剂硫酸铵、碳酸氢铵等深度为3％，浸泡时间约150～400天，由于浓度大、时间长，造成矿区土壤酸化板结，有机质流失。浸矿溶液侧渗和毛细管作用使植物根系受损、地表的很多草本植物枯死，对植被造成较为严重的破坏。对当地矿区的生态环境造成了极大的危害。
[0004]迫于环境压力，目前赣南稀土资源开采利用处于停产状态。在稀土资源开采利用的同时保护生态并修复环境是今后的必然要求。利用先锋植物进行矿区生态恢复是一种成本低且长期可持续的措施。先锋植物具有极其顽强的生命力，能适应离子型稀土矿区的极端条件。为了恢复矿区植被、改善生态环境，应首先种植耐性强的先锋草类，如黑麦草、芒草、狗牙根、百喜草、香根草、紫花苜蓿等，使裸地迅速被植物所覆盖，形成草丛群落，使土壤逐渐得以改良，从而逐步在矿区建立稳定和高产的生态系统。
[0005]然而，矿区的高氨氮、硫酸根和重金属污染往往会引发先锋植物的生态生理压力，并对植物的生存和生长产生负面影响，从而影响利用先锋植物进行生态恢复的效果。因此，确定先锋植物对矿区污染物的毒性阈值，对离子型稀土矿生态恢复的先锋植物的选择具有重要意义。

技术实现思路

[0006]为克服现有技术中存在的上述缺陷，本专利技术提供了一种测定离子型稀土矿先锋植物对污染物毒性阈值的方法，以确定先锋植物对污染物的毒性阈值，对离子型稀土矿生态恢复的先锋植物的选择具有重要意义。
[0007]为了实现上述专利技术目的，本专利技术提供以下技术方案：
[0008]本专利技术提供了一种测定离子型稀土矿先锋植物对污染物毒性阈值的方法，包括如下步骤：
[0009](1)测定离子型稀土矿污染物对先锋植物的毒性效应；
[0010](2)基于步骤(1)得到的毒性效应测定结果，与污染物浓度进行相关性分析，筛选与污染物显著相关的指标；
[0011](3)以污染物浓度为自变量x，步骤(2)中筛选出的显著相关指标为因变量y，为各指标选取最优的剂量效应模型进行拟合；
[0012](4)筛选与污染物浓度显著相关且剂量效应关系中决定系数＞0.8的指标作为综合毒性阈值评价终点指标，依据拟合函数计算每个评价终点指标的毒性阈值，并利用主成分分析法计算每个评价终点指标的权重；
[0013](5)依据评价终点指标的毒性阈值与权重计算植物对污染物的综合毒性阈值。
[0014]优选的，步骤(1)具体包括：
[0015](1.1)将等重量的离子型稀土矿区无污染土壤分别装入盆栽中，将盆栽分为不同梯度，每个梯度设置足额的生物学重复，分别加入梯度浓度且体积相同的污染物溶液后混合均匀，获得不同浓度的污染土壤；
[0016](1.2)在每个盆栽适宜且相同的深度播种同质量的种子，植物种植到快速生长期时进行采样，测定植物生长生理指标。
[0017]优选的，步骤(1.2)中所述的植物生长生理指标具体为：株高、地上地下生物量作为生长指标，叶绿素a、叶绿素b、胡萝卜素、可溶性糖、脯氨酸、过氧化氢、丙二醛、超氧化物歧化酶、过氧化物酶、谷胱甘肽过氧化物酶和谷胱甘肽作为生理指标。
[0018]优选的，步骤(3)中所述为各指标选取最优的剂量效应模型进行拟合的步骤具体为：
[0019](3.1)将R语言软件“剂量反应曲线”程序包中所包含的模型公式添加到origin的非线性拟合自定义函数中；
[0020](3.2)利用origin软件中的Rank models App为每个指标选择最优模型建立指标和污染物的关系模型。
[0021]优选的，步骤(3.1)中所述的模型公式包括：
[0022]log
‑
logistic：y＝c+(d
‑
c)/(1+exp(b*(log(x)
‑
log(e))))
f
；
[0023]Brain
‑
Cousens：y＝c+(d+f*x
‑
c)/(1+exp(b*(log(x)
‑
log(e))))；
[0024]Cedergreen
‑
Ritz
‑
Streibig：y＝c+(d+f*exp(
‑
1/x)
‑
c)/(1+exp(b*(log(x)
‑
log(e))))；
[0025]Weibull I:y＝c+(d
‑
c)exp(
‑
exp(b(log(x)
‑
log(e))))，
[0026]上述模型公式的参数化采用统一的结构，系数b表示剂量
‑
反应曲线的斜率，c、d表示反应的下限和上限，e表示有效剂量EC
50
。
[0027]优选的，步骤(3.2)中所述的选择最优模型的步骤具体为：
[0028]利用origin软件中的Rank models App，输入污染物浓度和对应的指标数据，勾选步骤(3.1)中添加到origin软件中的剂量效应曲线模型，软件会输出不同模型的使用拟合度(GoF)指数，包括赤池信息准则(AIC)、贝叶斯信息准则(BIC)、决定系数(R2)、残差平方和(RSS)、均方根误差(RSS/dof)，AIC、BIC、RSS、RSS/dof值越小，R2越接近1，则模型对数据的拟合越好，软件会依据使用拟合度指数对模型进行排序，排在前列的为最佳拟合模型。
[0029]优选的，步骤(4)中所述依据拟合函数计算每个评价终点指标的毒性阈值的方法具体为：
[0030](4.1)计算污染物胁迫造成指标相对于对照组产生x％的不利影响：
[0031]b＝a*(1
±
x％)，
[0032]式中：a表示污染物浓度为0时也即对照组对应的指标数值，“(1
±
x％)”表示若胁迫造成指标值下降为“(1
‑
x％)”，若胁迫造成指标值上升为“(1+x％)”；
[0033](4.2)利用计算得到的b代入指标与污染物的剂量效应关系曲线，计算出每个指标相对于对照组污染物胁迫造成x％的不利影响时对应的污染物浓度分别作为毒性阈值EC
x
。
[0034]优选的，步骤(4)中所述利用主成分分析法计算每个评价终点指本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种测定离子型稀土矿先锋植物对污染物毒性阈值的方法，其特征在于，包括如下步骤：(1)测定离子型稀土矿污染物对先锋植物的毒性效应；(2)基于步骤(1)得到的毒性效应测定结果，与污染物浓度进行相关性分析，筛选与污染物显著相关的指标；(3)以污染物浓度为自变量x，步骤(2)中筛选出的显著相关指标为因变量y，为各指标选取最优的剂量效应模型进行拟合；(4)筛选与污染物浓度显著相关且剂量效应关系中决定系数＞0.8的指标作为综合毒性阈值评价终点指标，依据拟合函数计算每个评价终点指标的毒性阈值，并利用主成分分析法计算每个评价终点指标的权重；(5)依据评价终点指标的毒性阈值与权重计算植物对污染物的综合毒性阈值。2.根据权利要求1所述的一种测定离子型稀土矿先锋植物对污染物毒性阈值的方法，其特征在于，步骤(1)具体包括：(1.1)将等重量的离子型稀土矿区无污染土壤分别装入盆栽中，将盆栽分为不同梯度，每个梯度设置足额的生物学重复，分别加入梯度浓度且体积相同的污染物溶液后混合均匀，获得不同浓度的污染土壤；(1.2)在每个盆栽适宜且相同的深度播种同质量的种子，植物种植到快速生长期时进行采样，测定植物生长生理指标。3.根据权利要求2所述的一种测定离子型稀土矿先锋植物对污染物毒性阈值的方法，其特征在于，步骤(1.2)中所述的植物生长生理指标具体为：株高、地上地下生物量作为生长指标，叶绿素a、叶绿素b、胡萝卜素、可溶性糖、脯氨酸、过氧化氢、丙二醛、超氧化物歧化酶、过氧化物酶、谷胱甘肽过氧化物酶和谷胱甘肽作为生理指标。4.根据权利要求1所述的一种测定离子型稀土矿先锋植物对污染物毒性阈值的方法，其特征在于，步骤(3)中所述为各指标选取最优的剂量效应模型进行拟合的步骤具体为：(3.1)将R语言软件“剂量反应曲线”程序包中所包含的模型公式添加到origin的非线性拟合自定义函数中；(3.2)利用origin软件中的Rank modelsApp为每个指标选择最优模型建立指标和污染物的关系模型。5.根据权利要求4所述的一种测定离子型稀土矿先锋植物对污染物毒性阈值的方法，其特征在于，步骤(3.1)中所述的模型公式包括：log
‑
logistic：y＝c+(d
‑
c)/(1+exp(b*(log(x)
‑
log(e))))
f
；Brain
‑
Cousens：y＝c+(d+f*x
‑
c)/(1+exp(b*(log(x)
‑
log(e))))；Cedergreen
‑
Ritz
‑
Streibig：y＝c+(d+f*exp(
‑
1/x)
‑
c)/(1+exp(b*(log(x)
‑
log(e))))；We...

【专利技术属性】
技术研发人员：周静，杨慧娴，夏睿智，周俊，梁家妮，费佳赛，
申请(专利权)人：江西洁地环境治理生态科技有限公司，
类型：发明
国别省市：