CN109308577A - 一种沉积物重金属环境质量评价方法 - Google Patents
一种沉积物重金属环境质量评价方法 Download PDFInfo
- Publication number
- CN109308577A CN109308577A CN201811060296.0A CN201811060296A CN109308577A CN 109308577 A CN109308577 A CN 109308577A CN 201811060296 A CN201811060296 A CN 201811060296A CN 109308577 A CN109308577 A CN 109308577A
- Authority
- CN
- China
- Prior art keywords
- sediments
- deposit
- heavy metals
- value
- content
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Economics (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Operations Research (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Processing Of Solid Wastes (AREA)
Abstract
本发明公开了一种沉积物重金属环境质量评价方法,收集沉积物的基础数据,包括沉积物重金属区域地球化学背景值和沉积物理化性质资料,具体有沉积物重金属含量、沉积物氧化铝含量、沉积物有机物质含量、沉积物中值粒径;通过三角模糊数和α截集,处理沉积物重金属含量和沉积物重金属区域地球化学背景值;计算沉积物重金属的富集指数,进行相关性分析,判断沉积物重金属来源;计算地累积指数,并基于隶属度计算,判断沉积物重金属污染水平;计算单因子污染指数,核算潜在生态危害指数,再基于隶属度计算,评价沉积物重金属生态风险。本发明克服了数据选取及评价结果不确定性的缺陷,合理评价沉积物重金属环境质量,为环境质量管理提供有效技术支持。
Description
技术领域
本发明涉及一种环境质量评价方法,具体涉及一种沉积物重金属环境质量评价方法,属于环保技术领域。
背景技术
随着人口增长和经济发展,沉积物被认为是水体污染物的“汇”,重金属作为危害严重的污染物,受到广泛关注。为对沉积物重金属来源、污染水平和生态风险进行评价,许多方法得以发展,包括多元统计方法、富集系数法、潜在生态危险指数法等,对污染物转移转化机制探究、沉积物风险评估和环境质量基准制定都具有重要意义。
但在以往的评价中,沉积物重金属含量及区域地球化学背景值的选用均选择定值,以沉积物重金属区域地球化学背景值为例,可选择研究区域勘探部门提供的背景值,可根据《土壤环境质量标准(GB15618-1995)》中的数据,也可依据期刊数据库检索,不同的区域地球化学背景值得出的结论存在差异。数据选取的不确定性造成结果的不确定性。
发明内容
为解决现有技术的不足,本发明的目的在于提供一种判断重金属来源,判断和评价沉积物重金属污染水平与生态风险的沉积物重金属环境质量评价方法。
为了实现上述目标,本发明采用如下的技术方案:
一种沉积物重金属环境质量评价方法,包括以下步骤:
S1、收集沉积物的基础数据,包括沉积物理化性质资料、沉积物重金属区域地球化学背景值;
沉积物理化性质资料,包括沉积物重金属含量、沉积物氧化铝含量、沉积物有机物质含量、沉积物中值粒径;
S2、通过三角模糊数和α截集,处理沉积物重金属含量和沉积物重金属区域地球化学背景值,处理后的数据作以下计算用:
S3、计算沉积物重金属的富集指数,对沉积物理化性质资料进行相关性分析,判断沉积物重金属来源;
S4、计算地累积指数,并基于隶属度计算,判断沉积物重金属污染水平;
S5、计算单因子污染指数,核算潜在生态危害指数,再基于隶属度计算,评价沉积物重金属生态风险。
上述步骤S1中的重金属包括铬、锰、铁、镍、铜、锌、镉、汞、铅。
上述步骤S1中的沉积物重金属区域地球化学背景值、沉积物重金属含量、沉积物氧化铝含量、沉积物有机物质含量、沉积物中值粒径为范围值;
所述沉积物重金属区域地球化学背景值,为范围值Cbackgroundi∈[cbackground1i,cbackground2i](mg/kg),来源包括《土壤环境质量标准(GB15618-1995)》、《土壤环境元素化学(中国环境科学出版社,1995)》、区域勘探部门、期刊数据库;其中,cbackground1i为最小值,cbackground2i为最大值,假设为[cbackground1i,cbackground2i]平均值;i=1,2,…,N;N为沉积物重金属种类数;
所述沉积物重金属含量,为范围值Csedimenti∈[csediment1i,csediment2i](mg/kg),来源包括区域勘探部门、期刊数据库、实地勘测;其中,csediment1i为最小值,csediment2i为最大值,假设为[csediment1i,csediment2i]平均值;
所述沉积物氧化铝含量,为范围值其中,为最小值,为最大值,假设为平均值;
所述沉积物有机物质含量,为范围值COMi∈[cOM1i,cOM2i](mg/kg);其中,cOM1i为最小值,cOM2i为最大值,假设为[cOM1i,cOM2i]平均值;所述有机物质含量的表征,包括有机质含量、烧失量或总有机碳;
所述沉积物中值粒径,为范围值Di∈[d1i,d2i](μm);其中d1i为最小值,d2i为最大值,假设为[d1i,d2i]平均值;
进一步的,上述区域勘探部门包括区域住房与城乡规划建设局、勘测研究院;
所述期刊数据库包括中国知网期刊全文数据库、维普中文科技期刊库、万方数据知识服务平台、SCI(ISI Web of Science)。
上述步骤S2中沉积物基础数据处理,包括以下步骤:
A1:三角模糊数处理
A11、沉积物重金属含量三角模糊数处理:
假设xsedimenti,ysedimenti,zsedimenti分别为沉积物第i种重金属含量的最小可能值、最可能值和最大可能值,构成三角模糊数Asedimenti=(xsedimenti,ysedimenti,zsedimenti);
假设σi为[csediment1i,csediment2i]标准差;
所述xsedimenti取csediment1i和较大值;
所述ysedimenti取
所述zsedimenti取csediment2i和较小值;
A12、沉积物重金属区域地球化学背景值三角模糊数处理:
同沉积物重金属含量三角模糊数处理,
构成的三角模糊数为:Abackgroundi=(xbackgroundi,ybackgroundi,zbackgroundi);
A2:α截集处理
A21、沉积物重金属含量α截集处理:
令Aα sedimenti=[(ysedimenti-xsedimenti)α+xsedimenti,zsedimenti-(zsedimenti-ysedimenti)α],α∈[0,1];
A22、沉积物重金属区域地球化学背景值α截集处理:
同沉积物重金属含量α截集处理,
Aα backgroundi=[(ybackgroundi-xbackgroundi)α+xbackgroundi,zbackgroundi-(zbackgroundi-ybackgroundi)α],α∈[0,1]。
上述步骤S3中沉积物重金属来源判断,包括以下步骤:
B1:富集指数计算
式(1)中,EFi为沉积物第i种重金属的富集指数,
为沉积物第i种重金属与沉积物重金属Fe含量平均值比值,为沉积物第i种重金属与沉积物重金属Fe区域地球化学背景值平均值比值;
若EFi≤1.5,沉积物重金属主要源自自然风化过程;
若EFi>1.5,沉积物重金属主要源自人为活动;
B2:相关性分析
利用数据处理软件,将沉积物重金属含量平均值、沉积物氧化铝含量平均值、沉积物有机物质含量平均值、沉积物中值粒径平均值进行相关性分析,获得相关性分析显著性检验的p值,当p≤0.05为显著相关,p≤0.01为极显著相关;
若沉积物重金属含量两两间呈显著或极显著正相关,则沉积物重金属具有相同来源;
若沉积物重金属含量与沉积物氧化铝含量呈显著或极显著正相关,则沉积物重金属主要源自自然风化过程;
若沉积物重金属含量与沉积物有机物质含量呈显著或极显著相关,则沉积物重金属含量受有机物质控制;
若沉积物重金属含量与沉积物中值粒径呈显著或极显著相关,则沉积物重金属含量受粒径控制;
进一步的,上述数据处理软件包括SPSS、Matlab。
上述步骤S4中沉积物重金属污染水平判断,包括以下步骤:
C1:地累积指数计算
式(2)中,Igeoi为沉积物第i种重金属的地累积指数,k为修正系数,取值k=1.5,Igeo1i为最小值,Igeo2i为最大值;
若Igeoi≤0,第0个评价等级,沉积物第i种重金属污染级数为0,污染水平为清洁;
若0<Igeoi≤1,第1个评价等级,沉积物第i种重金属污染级数为1,污染水平为轻度污染;
若1<Igeoi≤2,第2个评价等级,沉积物第i种重金属污染级数为2,污染水平为偏中污染;
若2<Igeoi≤3,第3个评价等级,沉积物第i种重金属污染级数为3,污染水平为中度污染;
若3<Igeoi≤4,第4个评价等级,沉积物第i种重金属污染级数为4,污染水平为偏重污染;
若4<Igeoi≤5,第5个评价等级,沉积物第i种重金属污染级数为5,污染水平为重度污染;
若Igeoi>5,第6个评价等级,沉积物第i种重金属污染级数为6,污染水平为严重污染;
C2:隶属度计算
式(3)中,A(li)为[Igeo1i,Igeo2i]对的隶属度,为第l个评价等级的地累积指数,l=0,1,…,6;
取A(li)最大值所在的评价等级进行沉积物第i种重金属污染水平判断。
上述步骤S5中沉积物重金属生态风险评价,包括以下步骤:
D1:单因子污染指数计算
式(4)中,Ci为沉积物第i种重金属的单因子污染指数;
D2:潜在生态危害指数计算
Ei=Ti×Ci=[E1i,E2i] (5)
式(5)中,Ei为沉积物第i种重金属的潜在生态风险参数,Ti为沉积物第i种重金属的毒性响应系数,E1i为最小值,E2i为最大值;
式(6)中,RI为潜在生态危害指数,N为沉积物重金属种类数,RI1为最小值,RI2为最大值;
若Ei≤40,第1个评价等级,沉积物第i种重金属生态风险为低生态风险;
若40<Ei≤80,第2个评价等级,沉积物第i种重金属生态风险为中生态风险;
若80<Ei≤160,第3个评价等级,沉积物第i种重金属生态风险为较高生态风险;
若160<Ei≤320,第4个评价等级,沉积物第i种重金属生态风险为高生态风险;
若Ei>320,第5个评价等级,沉积物第i种重金属生态风险为很高生态风险;
若RI≤150,第1个评价等级,沉积物重金属生态风险为低生态风险;
若150<RI≤300,第2个评价等级,沉积物重金属生态风险为中生态风险;
若300<RI≤600,第3个评价等级,沉积物重金属生态风险为较高生态风险;
若RI>600,第4个评价等级,沉积物重金属生态风险为高生态风险;
D3:隶属度计算
式(7)中,A(ri)为[E1i,E2i]对的隶属度,为第r个评价等级的潜在生态风险参数,r=1,2…,5;
取A(ri)最大值所在的评价等级进行沉积物第i种重金属生态风险评价;
式(8)中,A(m)为[RI1,RI2]对的隶属度,为第m个评价等级的潜在生态危险指数,m=1,2,3,4;
取A(m)最大值所在的评价等级进行沉积物重金属生态风险评价。
本发明的有益之处在于:
本发明的一种沉积物重金属环境质量评价的方法,采用富集指数和相关性分析进行沉积物重金属来源判断,将沉积物重金属含量和区域地球化学背景值定义为三角模糊数,将计算得到的地累积指数通过隶属度判断沉积物重金属污染水平,将计算得到的潜在生态危险参数和潜在生态危险指数通过隶属度评价沉积物重金属生态风险,克服了数据选取不确定性及评价结果不确定性的缺陷,为沉积物重金属评价提供一种重要参考,为环境质量管理提供一种有效技术支持。
附图说明
图1是本发明一种沉积物重金属环境质量评价方法的流程图。
具体实施方式
以下结合具体实施例对本发明作具体的介绍。
一种沉积物重金属环境质量评价方法,包括以下步骤:
S1、收集沉积物的基础数据,包括沉积物理化性质资料、沉积物重金属区域地球化学背景值;
沉积物理化性质资料,包括沉积物重金属含量、沉积物氧化铝含量、有机物质含量、中值粒径;
S2、通过三角模糊数和α截集,处理沉积物重金属含量和沉积物重金属区域地球化学背景值,处理后的数据作以下计算用:
S3、计算沉积物重金属的富集指数,对沉积物理化性质资料进行相关性分析,判断沉积物重金属来源;
S4、计算地累积指数,并基于隶属度计算,判断沉积物重金属污染水平;
S5、计算单因子污染指数,核算潜在生态危害指数,再基于隶属度计算,评价沉积物重金属生态风险。
上述步骤S1中的重金属包括铬、锰、铁、镍、铜、锌、镉、汞、铅。
上述步骤S1中的沉积物重金属区域地球化学背景值、沉积物重金属含量、沉积物氧化铝含量、沉积物有机物质含量、沉积物中值粒径为范围值;
沉积物重金属区域地球化学背景值,为范围值Cbackgroundi∈[cbackground1i,cbackground2i](mg/kg),来源包括《土壤环境质量标准(GB15618-1995)》、《土壤环境元素化学(中国环境科学出版社,1995)》、区域勘探部门、期刊数据库;其中,cbackground1i为最小值,cbackground2i为最大值,假设为[cbackground1i,cbackground2i]平均值;i=1,2,…,N;N为沉积物重金属种类数;
沉积物重金属含量,为范围值Csedimenti∈[csediment1i,csediment2i](mg/kg),来源包括区域勘探部门、期刊数据库、实地勘测;其中,csediment1i为最小值,csediment2i为最大值,假设为[csediment1i,csediment2i]平均值;
沉积物氧化铝含量,为范围值其中,为最小值,为最大值,假设为平均值;
沉积物有机物质含量,为范围值COMi∈[cOM1i,cOM2i](mg/kg);其中,cOM1i为最小值,cOM2i为最大值,假设为[cOM1i,cOM2i]平均值;有机物质含量的表征,包括有机质含量、烧失量或总有机碳;
沉积物中值粒径,为范围值Di∈[d1i,d2i](μm);其中d1i为最小值,d2i为最大值,假设为[d1i,d2i]平均值;
进一步的,区域勘探部门包括区域住房与城乡规划建设局、勘测研究院;
期刊数据库包括中国知网期刊全文数据库、维普中文科技期刊库、万方数据知识服务平台、ISI Web of Science。
上述步骤S2中沉积物基础数据处理,包括以下步骤:
A1:三角模糊数处理
A11、沉积物重金属含量三角模糊数处理:
假设xsedimenti,ysedimenti,zsedimenti分别为沉积物第i种重金属含量的最小可能值、最可能值和最大可能值,构成三角模糊数Asedimenti=(xsedimenti,ysedimenti,zsedimenti);
假设σi为[csediment1i,csediment2i]标准差;
xsedimenti取csediment1i和较大值;
ysedimenti取
zsedimenti取csediment2i和较小值;
A12、沉积物重金属区域地球化学背景值三角模糊数处理:
同沉积物重金属含量三角模糊数处理,
构成的三角模糊数为:Abackgroundi=(xbackgroundi,ybackgroundi,zbackgroundi);
A2:α截集处理
A21、沉积物重金属含量α截集处理:
令Aα sedimenti=[(ysedimenti-xsedimenti)α+xsedimenti,zsedimenti-(zsedimenti-ysedimenti)α],α∈[0,1];
A22、沉积物重金属区域地球化学背景值α截集处理:
同沉积物重金属含量α截集处理,
Aα backgroundi=[(ybackgroundi-xbackgroundi)α+xbackgroundi,zbackgroundi-(zbackgroundi-ybackgroundi)α],α∈[0,1]。
上述步骤S3中沉积物重金属来源判断,包括以下步骤:
B1:富集指数计算
式(1)中,EFi为沉积物第i种重金属的富集指数,
为沉积物第i种重金属与沉积物重金属Fe含量平均值比值,为沉积物第i种重金属与沉积物重金属Fe区域地球化学背景值平均值比值;
若EFi≤1.5,沉积物重金属主要源自自然风化过程;
若EFi>1.5,沉积物重金属主要源自人为活动;
B2:相关性分析
利用数据处理软件,将沉积物重金属含量平均值、沉积物氧化铝含量平均值、沉积物有机物质含量平均值、沉积物中值粒径平均值进行相关性分析,获得相关性分析显著性检验的p值,当p≤0.05为显著相关,p≤0.01为极显著相关;
若沉积物重金属含量两两间呈显著或极显著正相关,则沉积物重金属具有相同来源;
若沉积物重金属含量与沉积物氧化铝含量呈显著或极显著正相关,则沉积物重金属主要源自自然风化过程;
若沉积物重金属含量与沉积物有机物质含量呈显著或极显著相关,则沉积物重金属含量受有机物质控制;
若沉积物重金属含量与沉积物中值粒径呈显著或极显著相关,则沉积物重金属含量受粒径控制;
进一步的,数据处理软件包括SPSS、Matlab。
上述步骤S4中沉积物重金属污染水平判断,包括以下步骤:
C1:地累积指数计算
式(2)中,Igeoi为沉积物第i种重金属的地累积指数,k为修正系数,取值k=1.5,Igeo1i为最小值,Igeo2i为最大值;
若Igeoi≤0,第0个评价等级,沉积物第i种重金属污染级数为0,污染水平为清洁;
若0<Igeoi≤1,第1个评价等级,沉积物第i种重金属污染级数为1,污染水平为轻度污染;
若1<Igeoi≤2,第2个评价等级,沉积物第i种重金属污染级数为2,污染水平为偏中污染;
若2<Igeoi≤3,第3个评价等级,沉积物第i种重金属污染级数为3,污染水平为中度污染;
若3<Igeoi≤4,第4个评价等级,沉积物第i种重金属污染级数为4,污染水平为偏重污染;
若4<Igeoi≤5,第5个评价等级,沉积物第i种重金属污染级数为5,污染水平为重度污染;
若Igeoi>5,第6个评价等级,沉积物第i种重金属污染级数为6,污染水平为严重污染;
C2:隶属度计算
式(3)中,A(li)为[Igeo1i,Igeo2i]对的隶属度,为第l个评价等级的地累积指数,l=0,1,…,6;
取A(li)最大值所在的评价等级进行沉积物第i种重金属污染水平判断。
上述步骤S5中沉积物重金属生态风险评价,包括以下步骤:
D1:单因子污染指数计算
式(4)中,Ci为沉积物第i种重金属的单因子污染指数;
D2:潜在生态危害指数计算
Ei=Ti×Ci=[E1i,E2i] (5)
式(5)中,Ei为沉积物第i种重金属的潜在生态风险参数,Ti为沉积物第i种重金属的毒性响应系数,E1i为最小值,E2i为最大值;
式(6)中,RI为潜在生态危害指数,N为沉积物重金属种类数,RI1为最小值,RI2为最大值;
若Ei≤40,第1个评价等级,沉积物第i种重金属生态风险为低生态风险;
若40<Ei≤80,第2个评价等级,沉积物第i种重金属生态风险为中生态风险;
若80<Ei≤160,第3个评价等级,沉积物第i种重金属生态风险为较高生态风险;
若160<Ei≤320,第4个评价等级,沉积物第i种重金属生态风险为高生态风险;
若Ei>320,第5个评价等级,沉积物第i种重金属生态风险为很高生态风险;
若RI≤150,第1个评价等级,沉积物重金属生态风险为低生态风险;
若150<RI≤300,第2个评价等级,沉积物重金属生态风险为中生态风险;
若300<RI≤600,第3个评价等级,沉积物重金属生态风险为较高生态风险;
若RI>600,第4个评价等级,沉积物重金属生态风险为高生态风险;
D3:隶属度计算
式(7)中,A(ri)为[E1i,E2i]对的隶属度,为第r个评价等级的潜在生态风险参数,r=1,2…,5;
取A(ri)最大值所在的评价等级进行沉积物第i种重金属生态风险评价;
式(8)中,A(m)为[RI1,RI2]对的隶属度,为第m个评价等级的潜在生态危险指数,m=1,2,3,4;
取A(m)最大值所在的评价等级进行沉积物重金属生态风险评价。
实施例1:
以中国某一河流沉积物为例,
进行沉积物重金属环境质量评价,步骤如下:
1.沉积物基础数据收集及处理
将河流分成10个区段,河流每个区段沉积物重金属包括铬(Cr)、铜(Cu)、锌(Zn)、镉(Cd)、汞(Hg)、铅(Pb),河流第1个区段沉积物数据如表1所示;
表1沉积物重金属数据(单位:mg/kg)
将沉积物重金属数据进行三角模糊数处理,数据如表2所示;
表2三角模糊数处理后沉积物重金属数据(单位:mg/kg)
进一步将数据进行α截集处理,α取0.9,数据如表3所示;
表3α截集处理后沉积物重金属数据(单位:mg/kg)
2.沉积物重金属来源判断
河流第1个区段沉积物富集系数,如表4所示;
表4富集系数
河流第1个区段Cr的富集系数小于1.5,说明河流第1个区段Cr主要源自自然风化过程,其余重金属主要源自人类活动;
利用SPSS软件,进行河流沉积物重金属含量平均值与沉积物氧化铝含量平均值、沉积物有机物质含量平均值、沉积物中值粒径平均值相关性分析,结果如表5所示;
表5相关系数
*p≤0.05,显著相关
**p≤0.01,极显著相关
该河流沉积物重金属Cu含量和Zn含量两两间呈极显著正相关,说明沉积物重金属Cu和Zn具有相同来源;沉积物重金属含量与沉积物氧化铝含量相关性不显著;沉积物重金属含量与沉积物有机物质含量相关性不显著;沉积物重金属Zn含量与沉积物中值粒径呈显著负相关,说明沉积物重金属Zn一定程度上受粒径控制。
3.沉积物重金属污染水平判断
考虑成岩作用引起的背景值变动,修正系数取k=1.5,河流第1个区段地累积指数如表6所示;地累积指数隶属度计算如表7所示,故河流第1个区段重金属Cr污染水平为清洁,Cu污染水平为重度污染,Zn污染水平为偏重污染,Cd污染水平为严重污染,Hg污染水平为偏重污染,Pb污染水平为偏重污染;
表6地累积指数
表7地累积指数隶属度
4.沉积物重金属生态风险评价
河流第1个区段单因子污染指数计算结果如表8所示;
表8单因子污染指数
铬(Cr)、铜(Cu)、锌(Zn)、镉(Cd)、汞(Hg)、铅(Pb)毒性响应系数分别为2、5、1、30、40、5,单个沉积物重金属潜在生态风险参数如表9所示;单个沉积物重金属潜在生态风险参数隶属度计算如表10所示,故河流第1个区段重金属Cr生态风险为低生态风险,Cu生态风险为较高生态风险,Zn生态风险为低生态风险,Cd生态风险为很高生态风险,Hg生态风险为很高生态风险,Pb生态风险为中生态风险;
表9单个沉积物重金属潜在生态风险参数
表10单个沉积物重金属潜在生态风险参数隶属度
河流第1个区段重金属潜在生态危害指数RI=[2319.914,3038.72],沉积物重金属生态风险为高生态风险。
以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解,上述实施例不以任何形式限制本发明,凡采用等同替换或等效变换的方式所获得的技术方案,均落在本发明的保护范围内。
Claims (9)
1.一种沉积物重金属环境质量评价方法,其特征在于,包括以下步骤:
S1、收集沉积物的基础数据,包括沉积物理化性质资料、沉积物重金属区域地球化学背景值;
沉积物理化性质资料,包括沉积物重金属含量、沉积物氧化铝含量、沉积物有机物质含量、沉积物中值粒径;
S2、通过三角模糊数和α截集,处理沉积物重金属含量和沉积物重金属区域地球化学背景值,处理后的数据作以下计算用:
S3、计算沉积物重金属的富集指数,对沉积物理化性质资料进行相关性分析,判断沉积物重金属来源;
S4、计算地累积指数,并基于隶属度计算,判断沉积物重金属污染水平;
S5、计算单因子污染指数,核算潜在生态危害指数,再基于隶属度计算,评价沉积物重金属生态风险。
2.根据权利要求1所述的一种沉积物重金属环境质量评价方法,其特征在于,所述步骤S1中的重金属包括铬、锰、铁、镍、铜、锌、镉、汞、铅。
3.根据权利要求1所述的一种沉积物重金属环境质量评价方法,其特征在于,所述步骤S1中的沉积物重金属区域地球化学背景值、沉积物重金属含量、沉积物氧化铝含量、沉积物有机物质含量、沉积物中值粒径为范围值;
所述沉积物重金属区域地球化学背景值,为范围值Cbackgroundi∈[cbackground1i,cbackground2i](mg/kg),来源包括《土壤环境质量标准(GB15618-1995)》、《土壤环境元素化学(中国环境科学出版社,1995)》、区域勘探部门、期刊数据库;其中,cbackground1i为最小值,cbackground2i为最大值,假设为[cbackground1i,cbackground2i]平均值;i=1,2,…,N;N为沉积物重金属种类数;
所述沉积物重金属含量,为范围值Csedimenti∈[csediment1i,csediment2i](mg/kg),来源包括区域勘探部门、期刊数据库、实地勘测;其中,csediment1i为最小值,csediment2i为最大值,假设为[csediment1i,csediment2i]平均值;
所述沉积物氧化铝含量,为范围值其中,为最小值,为最大值,假设为平均值;
所述沉积物有机物质含量,为范围值COMi∈[cOM1i,cOM2i](mg/kg);其中,cOM1i为最小值,cOM2i为最大值,假设为[cOM1i,cOM2i]平均值;所述有机物质含量的表征,包括有机质含量、烧失量或总有机碳;
所述沉积物中值粒径,为范围值Di∈[d1i,d2i](μm);其中d1i为最小值,d2i为最大值,假设为[d1i,d2i]平均值。
4.根据权利要求3所述的一种沉积物重金属环境质量评价方法,其特征在于,所述区域勘探部门包括区域住房与城乡规划建设局、勘测研究院;
所述期刊数据库包括中国知网期刊全文数据库、维普中文科技期刊库、万方数据知识服务平台、SCI。
5.根据权利要求1所述的一种沉积物重金属环境质量评价方法,其特征在于,所述步骤S2中沉积物基础数据处理,包括以下步骤:
A1:三角模糊数处理
A11、沉积物重金属含量三角模糊数处理:
假设xsedimenti,ysedimenti,zsedimenti分别为沉积物第i种重金属含量的最小可能值、最可能值和最大可能值,构成三角模糊数Asedimenti=(xsedimenti,ysedimenti,zsedimenti);
假设σi为[csediment1i,csediment2i]标准差;
所述xsedimenti取csediment1i和较大值;
所述ysedimenti取
所述zsedimenti取csediment2i和较小值;
A12、沉积物重金属区域地球化学背景值三角模糊数处理:
同沉积物重金属含量三角模糊数处理,
构成的三角模糊数为:Abackgroundi=(xbackgroundi,ybackgroundi,zbackgroundi);
A2:α截集处理
A21、沉积物重金属含量α截集处理:
令Aα sedimenti=[(ysedimenti-xsedimenti)α+xsedimenti,zsedimenti-(zsedimenti-ysedimenti)α],α∈[0,1];
A22、沉积物重金属区域地球化学背景值α截集处理:
同沉积物重金属含量α截集处理,
Aα backgroundi=[(ybackgroundi-xbackgroundi)α+xbackgroundi,zbackgroundi-(zbackgroundi-ybackgroundi)α],α∈[0,1]。
6.根据权利要求1所述的一种沉积物重金属环境质量评价方法,其特征在于,所述步骤S3中沉积物重金属来源判断,包括以下步骤:
B1:富集指数计算
式(1)中,EFi为沉积物第i种重金属的富集指数,
为沉积物第i种重金属与沉积物重金属Fe含量平均值比值,为沉积物第i种重金属与沉积物重金属Fe区域地球化学背景值平均值比值;
若EFi≤1.5,沉积物重金属主要源自自然风化过程;
若EFi>1.5,沉积物重金属主要源自人为活动;
B2:相关性分析
利用数据处理软件,将沉积物重金属含量平均值、沉积物氧化铝含量平均值、沉积物有机物质含量平均值、沉积物中值粒径平均值进行相关性分析,获得相关性分析显著性检验的p值,当p≤0.05为显著相关,p≤0.01为极显著相关;
若沉积物重金属含量两两间呈显著或极显著正相关,则沉积物重金属具有相同来源;
若沉积物重金属含量与沉积物氧化铝含量呈显著或极显著正相关,则沉积物重金属主要源自自然风化过程;
若沉积物重金属含量与沉积物有机物质含量呈显著或极显著相关,则沉积物重金属含量受有机物质控制;
若沉积物重金属含量与沉积物中值粒径呈显著或极显著相关,则沉积物重金属含量受粒径控制。
7.根据权利要求6所述的一种沉积物重金属环境质量评价方法,其特征在于,所述数据处理软件包括SPSS、Matlab。
8.根据权利要求1所述的一种沉积物重金属环境质量评价方法,其特征在于,所述步骤S4中沉积物重金属污染水平判断,包括以下步骤:
C1:地累积指数计算
式(2)中,Igeoi为沉积物第i种重金属的地累积指数,k为修正系数,取值k=1.5,Igeo1i为最小值,Igeo2i为最大值;
若Igeoi≤0,第0个评价等级,沉积物第i种重金属污染级数为0,污染水平为清洁;
若0<Igeoi≤1,第1个评价等级,沉积物第i种重金属污染级数为1,污染水平为轻度污染;
若1<Igeoi≤2,第2个评价等级,沉积物第i种重金属污染级数为2,污染水平为偏中污染;
若2<Igeoi≤3,第3个评价等级,沉积物第i种重金属污染级数为3,污染水平为中度污染;
若3<Igeoi≤4,第4个评价等级,沉积物第i种重金属污染级数为4,污染水平为偏重污染;
若4<Igeoi≤5,第5个评价等级,沉积物第i种重金属污染级数为5,污染水平为重度污染;
若Igeoi>5,第6个评价等级,沉积物第i种重金属污染级数为6,污染水平为严重污染;
C2:隶属度计算
式(3)中,A(li)为[Igeo1i,Igeo2i]对的隶属度,为第l个评价等级的地累积指数,l=0,1,…,6;
取A(li)最大值所在的评价等级进行沉积物第i种重金属污染水平判断。
9.根据权利要求1所述的一种沉积物重金属环境质量评价方法,其特征在于,所述步骤S5中沉积物重金属生态风险评价,包括以下步骤:
D1:单因子污染指数计算
式(4)中,Ci为沉积物第i种重金属的单因子污染指数;
D2:潜在生态危害指数计算
Ei=Ti×Ci=[E1i,E2i] (5)
式(5)中,Ei为沉积物第i种重金属的潜在生态风险参数,Ti为沉积物第i种重金属的毒性响应系数,E1i为最小值,E2i为最大值;
式(6)中,RI为潜在生态危害指数,N为沉积物重金属种类数,RI1为最小值,RI2为最大值;
若Ei≤40,第1个评价等级,沉积物第i种重金属生态风险为低生态风险;
若40<Ei≤80,第2个评价等级,沉积物第i种重金属生态风险为中生态风险;
若80<Ei≤160,第3个评价等级,沉积物第i种重金属生态风险为较高生态风险;
若160<Ei≤320,第4个评价等级,沉积物第i种重金属生态风险为高生态风险;
若Ei>320,第5个评价等级,沉积物第i种重金属生态风险为很高生态风险;
若RI≤150,第1个评价等级,沉积物重金属生态风险为低生态风险;
若150<RI≤300,第2个评价等级,沉积物重金属生态风险为中生态风险;
若300<RI≤600,第3个评价等级,沉积物重金属生态风险为较高生态风险;
若RI>600,第4个评价等级,沉积物重金属生态风险为高生态风险;
D3:隶属度计算
式(7)中,A(ri)为[E1i,E2i]对的隶属度,为第r个评价等级的潜在生态风险参数,r=1,2…,5;
取A(ri)最大值所在的评价等级进行沉积物第i种重金属生态风险评价;
式(8)中,A(m)为[RI1,RI2]对的隶属度,为第m个评价等级的潜在生态危险指数,m=1,2,3,4;
取A(m)最大值所在的评价等级进行沉积物重金属生态风险评价。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811060296.0A CN109308577B (zh) | 2018-09-12 | 2018-09-12 | 一种沉积物重金属环境质量评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811060296.0A CN109308577B (zh) | 2018-09-12 | 2018-09-12 | 一种沉积物重金属环境质量评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109308577A true CN109308577A (zh) | 2019-02-05 |
CN109308577B CN109308577B (zh) | 2021-11-30 |
Family
ID=65224797
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811060296.0A Active CN109308577B (zh) | 2018-09-12 | 2018-09-12 | 一种沉积物重金属环境质量评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109308577B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110135714A (zh) * | 2019-04-30 | 2019-08-16 | 清华大学 | 一种河流、湖泊沉积物重金属生态毒性风险的综合评价方法 |
CN110781225A (zh) * | 2019-10-25 | 2020-02-11 | 中国环境科学研究院 | 一种环境介质污染物浓度水平的诊断方法 |
CN112557612A (zh) * | 2020-11-20 | 2021-03-26 | 中南大学 | 利用水系沉积物解析金属矿区地下水重金属污染源及污染边界的方法 |
CN112818284A (zh) * | 2021-02-02 | 2021-05-18 | 中国水利水电科学研究院 | 一种沉积物中重金属污染物的生态风险评价方法 |
CN112836842A (zh) * | 2019-11-25 | 2021-05-25 | 天津大学 | 基于“源-汇”风险分析的流域水环境质量预测方法及系统 |
CN112851225A (zh) * | 2021-03-26 | 2021-05-28 | 河海大学 | 一种粉质粘土在盾构同步注浆中环保处理方法 |
CN113742635A (zh) * | 2020-05-29 | 2021-12-03 | 中国科学院声学研究所 | 一种干扰背景下k分布形状参数自适应估计方法及其系统 |
CN116596308A (zh) * | 2023-05-15 | 2023-08-15 | 生态环境部南京环境科学研究所 | 河流、湖泊沉积物重金属生态毒性风险的综合评价方法 |
CN117786616A (zh) * | 2024-02-23 | 2024-03-29 | 自然资源部第二海洋研究所 | 一种近海沉积物重金属污染特征确定方法、装置及介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103886217A (zh) * | 2014-04-04 | 2014-06-25 | 江苏省环境科学研究院 | 一种河湖沉积物中重金属污染的生态风险确定方法 |
-
2018
- 2018-09-12 CN CN201811060296.0A patent/CN109308577B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103886217A (zh) * | 2014-04-04 | 2014-06-25 | 江苏省环境科学研究院 | 一种河湖沉积物中重金属污染的生态风险确定方法 |
Non-Patent Citations (1)
Title |
---|
樊梦佳: "基于三角模糊数的河流沉积物重金属评价研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110135714B (zh) * | 2019-04-30 | 2021-07-13 | 清华大学 | 一种河流、湖泊沉积物重金属生态毒性风险的综合评价方法 |
CN110135714A (zh) * | 2019-04-30 | 2019-08-16 | 清华大学 | 一种河流、湖泊沉积物重金属生态毒性风险的综合评价方法 |
CN110781225A (zh) * | 2019-10-25 | 2020-02-11 | 中国环境科学研究院 | 一种环境介质污染物浓度水平的诊断方法 |
CN110781225B (zh) * | 2019-10-25 | 2020-06-19 | 中国环境科学研究院 | 一种环境介质污染物浓度水平的诊断方法 |
CN112836842A (zh) * | 2019-11-25 | 2021-05-25 | 天津大学 | 基于“源-汇”风险分析的流域水环境质量预测方法及系统 |
CN113742635B (zh) * | 2020-05-29 | 2023-08-08 | 中国科学院声学研究所 | 一种干扰背景下k分布形状参数自适应估计方法及其系统 |
CN113742635A (zh) * | 2020-05-29 | 2021-12-03 | 中国科学院声学研究所 | 一种干扰背景下k分布形状参数自适应估计方法及其系统 |
CN112557612A (zh) * | 2020-11-20 | 2021-03-26 | 中南大学 | 利用水系沉积物解析金属矿区地下水重金属污染源及污染边界的方法 |
CN112557612B (zh) * | 2020-11-20 | 2022-06-03 | 中南大学 | 利用水系沉积物解析金属矿区地下水重金属污染源及污染边界的方法 |
CN112818284A (zh) * | 2021-02-02 | 2021-05-18 | 中国水利水电科学研究院 | 一种沉积物中重金属污染物的生态风险评价方法 |
CN112818284B (zh) * | 2021-02-02 | 2023-09-05 | 中国水利水电科学研究院 | 一种沉积物中重金属污染物的生态风险评价方法 |
CN112851225B (zh) * | 2021-03-26 | 2021-11-12 | 河海大学 | 一种粉质粘土在盾构同步注浆中环保处理方法 |
CN112851225A (zh) * | 2021-03-26 | 2021-05-28 | 河海大学 | 一种粉质粘土在盾构同步注浆中环保处理方法 |
CN116596308A (zh) * | 2023-05-15 | 2023-08-15 | 生态环境部南京环境科学研究所 | 河流、湖泊沉积物重金属生态毒性风险的综合评价方法 |
CN116596308B (zh) * | 2023-05-15 | 2023-11-17 | 生态环境部南京环境科学研究所 | 河流、湖泊沉积物重金属生态毒性风险的综合评价方法 |
CN117786616A (zh) * | 2024-02-23 | 2024-03-29 | 自然资源部第二海洋研究所 | 一种近海沉积物重金属污染特征确定方法、装置及介质 |
CN117786616B (zh) * | 2024-02-23 | 2024-05-10 | 自然资源部第二海洋研究所 | 一种近海沉积物重金属污染特征确定方法、装置及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109308577B (zh) | 2021-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109308577A (zh) | 一种沉积物重金属环境质量评价方法 | |
Cheng et al. | Overview of trace metals in the urban soil of 31 metropolises in China | |
Dong et al. | How public and government matter in industrial pollution mitigation performance: Evidence from China | |
US20160110835A1 (en) | A method for determining ecological risks of heavy metal pollution in river and lake sediments | |
CN101944160A (zh) | 基于层次分析法和综合评价法建立的近岸海域生态环境综合评价方法 | |
Oliveira et al. | Association of periphytic diatom species of artificial substrate in lotic environments in the Arroio Sampaio basin, RS, Brazil: relationships with abiotic variables | |
Goe et al. | System tradeoffs in siting a solar photovoltaic material recovery infrastructure | |
Sueyoshi et al. | New concepts for environment-health measurement by data envelopment analysis and an application in China | |
CN111680890A (zh) | 一种适用于河长制的地表水环境质量评价方法 | |
CN110781225B (zh) | 一种环境介质污染物浓度水平的诊断方法 | |
Jupke et al. | Evaluating the biological validity of European river typology systems with least disturbed benthic macroinvertebrate communities | |
CN110885132B (zh) | 一种农村生活污水处理设施优先监管集的识别方法 | |
CN110790368B (zh) | 一种农村生活污水处理设施的区域监管方法 | |
He et al. | Status assessment and probabilistic health risk modeling of polycyclic aromatic hydrocarbons (PAHs) in surface soil across China | |
CN114066243B (zh) | 一种多指标耦合的地表型水源地安全风险评估方法 | |
Marshall et al. | Long-term phytoplankton trends and related water quality trends in the lower Chesapeake Bay, Virginia, USA | |
Parry et al. | Trace metal surveys of soil as a component of strategic and local planning policy development | |
CN115219683A (zh) | 一种水环境质量预警系统及其构建方法 | |
Shi et al. | Spatio-temporal variation of total mercury in precipitation in the largest industrial base in China: Impacts of meteorological factors and anthropogenic activities | |
Widarman et al. | Analytical Hierarchy Process Algorithm for Define of Water Meter | |
Haghparast et al. | Comprehensive Environmental Monitoring based on Stations of Environmental Pollutants (Air, Water and Soil) in Tehran | |
Mayack et al. | The effects of low concentrations of particulates from paper mill effluent on the macroinvertebrate community of a fast-flowing stream | |
Brown | Nutrient Loadings to Utah Lake from Bulk Atmospheric Deposition | |
Emerson et al. | Measuring toxic emissions from landfills using sequential screening | |
You et al. | Polycyclic aromatic hydrocarbons (PAHs) pollution and risk assessment of soils at contaminated sites in China over the past two decades |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |