【技术实现步骤摘要】
一种多环芳烃综合生态风险评价和敏感性分析优化方法
本专利技术属于多环芳烃(PAHs)生态风险评价，具体涉及一种多环芳烃综合生态风险评价和敏感性分析优化方法，尤其涉及一种基于层次阿基米德Copula函数和单一生态风险评价指标耦合的多环芳烃(PAHs)综合生态风险评价和敏感性分析的优化方法。
技术介绍
多环芳烃(PAHs)是一类具有半挥发性、持久性的有机污染物，具有强烈的毒性、致癌性、致畸性与诱变性，常产生于热源或石油源，且广泛存在于环境中，危害动植物体的生命安全。据其对生态环境的危害性，美国环保署将16种PAHs列如优先控制污染物名单，其中7种被国际癌症研究机构认定对人类有潜在致癌作用。分析多环芳烃的来源及风险程度对于研究其对生态环境的影响有着重要意义。生态风险评价是指评估或评价一个生态系统中特定事件发生可能性的过程。为了评估生态环境中多环芳烃的单一和整体危害，近年来出现了多种针对多环芳烃等有毒物质生态风险的评价方法，并得到了广泛应用。应用最为广泛评估方法包括风险商法(RQ)、苯并[a]芘毒性当量法(TEQBaP)、概率风险评...
【技术保护点】
1.一种多环芳烃综合生态风险评价和敏感性分析优化方法，其特征在于：所述方法基于层次阿基米德Copula函数和单一生态风险评价指标的耦合进行分析，包括如下步骤：/n(1)获取PAHs样本浓度数据，并记录样本编号及取样点位置，包括萘，苊，苊烯，芴，菲，蒽，荧蒽，芘，苯并[a]蒽，
【技术特征摘要】
1.一种多环芳烃综合生态风险评价和敏感性分析优化方法，其特征在于：所述方法基于层次阿基米德Copula函数和单一生态风险评价指标的耦合进行分析，包括如下步骤：
(1)获取PAHs样本浓度数据，并记录样本编号及取样点位置，包括萘，苊，苊烯，芴，菲，蒽，荧蒽，芘，苯并[a]蒽，，苯并[b]荧蒽，苯并[k]荧蒽，苯并[a]芘，二苯并[a，h]蒽，茚并[1，2，3-cd]芘以及苯并[g，h，i]芘的样本浓度；
(2)PAHs暴露浓度一致性度量，对于两变量间的一致性度量，若浓度变量均符合正态分布假定，采用Pearson相关系数计算；若浓度变量为非正态分布数据，采用Spearman秩相关系数计算；对于多变量间的一致性度量，采用Kendall秩相关系数计算，所述多变量为包含3个及其以上的变量；
(3)基于层次阿基米德Copula函数建立PAHs暴露浓度综合概率模型；
(4)选择单一的生态指标值带入PAHs暴露浓度综合概率模型计算综合生态风险程度；
(5)通过降低各组风险概率的固定步长对PAHs综合生态风险进行敏感性分析，并根据敏感性分析结果进行风险分级，计算公式如下：
S(f(x′))＝E(f(x))-E(f(x′))
式中，S(f(x′))为受试PAH的敏感性，E(f(x))表示初始PAH的联合生态风险程度，E(f(x′))表示受试PAH的联合生态风险程度。
2.根据权利要求1所述的多环芳烃综合生态风险评价和敏感性分析优化方法，其特征在于：步骤(2)对于Pearson相关系数计算和Spearman秩相关系数计算过程如下：
将PAHs按照苯环数目不同分为四组：2-ring，3-ring，4-ring及5-ring，所述5-ring组别中苯环数目大于等于5；采用多变量Kendall协和系数对每组分别进行多变量一致性度量；
对于符合连续正态分布的样本序列(Xn，Yn)，其Pearson相关系数为：
对于非正态分布的样本序列(Xn，Yn)，令Ri表示Xi在(X1，X2，...，Xn)中的秩，Qi表示Yi在(Y1，Y2，...，Y)中的秩，则有Spearman相关系数为：
又
则rs可以简化为：
假设有k个变量X1，x2，...，Xk，每个变量有n个观测值，设第j个变量Xj＝(X1j，X2j，...，xnj)，记Rij为Xij在(X1j，X2j，...，Xnj)的秩，则多变量Kendall协和系数W可以表示为：
式中...
【专利技术属性】
技术研发人员：王栋，刘文月，王远坤，曾献奎，祝晓彬，吴吉春，
申请(专利权)人：南京大学，
类型：发明
国别省市：江苏;32
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。