CN111445109A - 一种多环芳烃综合生态风险评价和敏感性分析优化方法 - Google Patents
一种多环芳烃综合生态风险评价和敏感性分析优化方法 Download PDFInfo
- Publication number
- CN111445109A CN111445109A CN202010138645.7A CN202010138645A CN111445109A CN 111445109 A CN111445109 A CN 111445109A CN 202010138645 A CN202010138645 A CN 202010138645A CN 111445109 A CN111445109 A CN 111445109A
- Authority
- CN
- China
- Prior art keywords
- comprehensive
- risk
- pahs
- concentration
- ecological
- 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
- 125000005575 polycyclic aromatic hydrocarbon group Chemical group 0.000 title claims abstract description 104
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000011156 evaluation Methods 0.000 title claims abstract description 35
- 238000010206 sensitivity analysis Methods 0.000 title claims abstract description 31
- 238000005457 optimization Methods 0.000 title claims description 15
- GVEPBJHOBDJJJI-UHFFFAOYSA-N fluoranthene Chemical compound C1=CC(C2=CC=CC=C22)=C3C2=CC=CC3=C1 GVEPBJHOBDJJJI-UHFFFAOYSA-N 0.000 claims description 34
- BBEAQIROQSPTKN-UHFFFAOYSA-N pyrene Chemical compound C1=CC=C2C=CC3=CC=CC4=CC=C1C2=C43 BBEAQIROQSPTKN-UHFFFAOYSA-N 0.000 claims description 32
- MWPLVEDNUUSJAV-UHFFFAOYSA-N anthracene Chemical compound C1=CC=CC2=CC3=CC=CC=C3C=C21 MWPLVEDNUUSJAV-UHFFFAOYSA-N 0.000 claims description 26
- 238000004364 calculation method Methods 0.000 claims description 24
- 125000005605 benzo group Chemical group 0.000 claims description 18
- 238000005259 measurement Methods 0.000 claims description 12
- 125000001997 phenyl group Chemical group [H]C1=C([H])C([H])=C(*)C([H])=C1[H] 0.000 claims description 11
- 238000012502 risk assessment Methods 0.000 claims description 11
- 239000013049 sediment Substances 0.000 claims description 11
- YNPNZTXNASCQKK-UHFFFAOYSA-N phenanthrene Chemical compound C1=CC=C2C3=CC=CC=C3C=CC2=C1 YNPNZTXNASCQKK-UHFFFAOYSA-N 0.000 claims description 8
- 230000035945 sensitivity Effects 0.000 claims description 8
- 150000001875 compounds Chemical class 0.000 claims description 7
- NIHNNTQXNPWCJQ-UHFFFAOYSA-N fluorene Chemical compound C1=CC=C2CC3=CC=CC=C3C2=C1 NIHNNTQXNPWCJQ-UHFFFAOYSA-N 0.000 claims description 7
- UFWIBTONFRDIAS-UHFFFAOYSA-N Naphthalene Chemical compound C1=CC=CC2=CC=CC=C21 UFWIBTONFRDIAS-UHFFFAOYSA-N 0.000 claims description 6
- 238000003646 Spearman's rank correlation coefficient Methods 0.000 claims description 6
- -1 acenaphthene Natural products 0.000 claims description 6
- 238000007476 Maximum Likelihood Methods 0.000 claims description 5
- 230000008878 coupling Effects 0.000 claims description 5
- 238000010168 coupling process Methods 0.000 claims description 5
- 238000005859 coupling reaction Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 4
- 238000005315 distribution function Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- SXQBHARYMNFBPS-UHFFFAOYSA-N Indeno[1,2,3-cd]pyrene Chemical compound C=1C(C2=CC=CC=C22)=C3C2=CC=C(C=C2)C3=C3C2=CC=CC3=1 SXQBHARYMNFBPS-UHFFFAOYSA-N 0.000 claims description 3
- HXGDTGSAIMULJN-UHFFFAOYSA-N acetnaphthylene Natural products C1=CC(C=C2)=C3C2=CC=CC3=C1 HXGDTGSAIMULJN-UHFFFAOYSA-N 0.000 claims description 3
- RKZJFAONHOYOKU-UHFFFAOYSA-N indeno[1,2,3-cd]pyrene Natural products C1Cc2ccc3ccc4c5ccccc5c6cc(=C1)c2c3c46 RKZJFAONHOYOKU-UHFFFAOYSA-N 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 abstract description 2
- 230000007547 defect Effects 0.000 abstract description 2
- 230000009467 reduction Effects 0.000 description 14
- 241000039077 Copula Species 0.000 description 7
- 206010064911 Pulmonary arterial hypertension Diseases 0.000 description 5
- 208000020193 Pulmonary artery hypoplasia Diseases 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 231100000419 toxicity Toxicity 0.000 description 5
- 230000001988 toxicity Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- HAXBIWFMXWRORI-UHFFFAOYSA-N Benzo[k]fluoranthene Chemical compound C1=CC(C2=CC3=CC=CC=C3C=C22)=C3C2=CC=CC3=C1 HAXBIWFMXWRORI-UHFFFAOYSA-N 0.000 description 3
- 239000003651 drinking water Substances 0.000 description 3
- 235000020188 drinking water Nutrition 0.000 description 3
- 239000003344 environmental pollutant Substances 0.000 description 3
- 231100000719 pollutant Toxicity 0.000 description 3
- 206010007269 Carcinogenicity Diseases 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- FTOVXSOBNPWTSH-UHFFFAOYSA-N benzo[b]fluoranthene Chemical compound C12=CC=CC=C1C1=CC3=CC=CC=C3C3=C1C2=CC=C3 FTOVXSOBNPWTSH-UHFFFAOYSA-N 0.000 description 2
- 231100000260 carcinogenicity Toxicity 0.000 description 2
- 230000007670 carcinogenicity Effects 0.000 description 2
- WDECIBYCCFPHNR-UHFFFAOYSA-N chrysene Chemical compound C1=CC=CC2=CC=C3C4=CC=CC=C4C=CC3=C21 WDECIBYCCFPHNR-UHFFFAOYSA-N 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000013210 evaluation model Methods 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 231100000048 toxicity data Toxicity 0.000 description 2
- 230000003313 weakening effect Effects 0.000 description 2
- UHOVQNZJYSORNB-UHFFFAOYSA-N Benzene Chemical compound C1=CC=CC=C1 UHOVQNZJYSORNB-UHFFFAOYSA-N 0.000 description 1
- 244000025254 Cannabis sativa Species 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- PXGOKWXKJXAPGV-UHFFFAOYSA-N Fluorine Chemical compound FF PXGOKWXKJXAPGV-UHFFFAOYSA-N 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 206010043275 Teratogenicity Diseases 0.000 description 1
- 231100000770 Toxic Equivalency Factor Toxicity 0.000 description 1
- 125000004054 acenaphthylenyl group Chemical group C1(=CC2=CC=CC3=CC=CC1=C23)* 0.000 description 1
- 230000003042 antagnostic effect Effects 0.000 description 1
- 230000008485 antagonism Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000711 cancerogenic effect Effects 0.000 description 1
- 231100000315 carcinogenic Toxicity 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000002542 deteriorative effect Effects 0.000 description 1
- 229910052731 fluorine Inorganic materials 0.000 description 1
- 239000011737 fluorine Substances 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 239000013505 freshwater Substances 0.000 description 1
- 239000000383 hazardous chemical Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 231100000299 mutagenicity Toxicity 0.000 description 1
- 230000007886 mutagenicity Effects 0.000 description 1
- 230000002688 persistence Effects 0.000 description 1
- 239000002957 persistent organic pollutant Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 231100000614 poison Toxicity 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000000246 remedial effect Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000013349 risk mitigation Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000002195 synergetic effect Effects 0.000 description 1
- 231100000211 teratogenicity Toxicity 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 231100000331 toxic Toxicity 0.000 description 1
- 231100000816 toxic dose Toxicity 0.000 description 1
- 230000002588 toxic effect Effects 0.000 description 1
- 231100001234 toxic pollutant Toxicity 0.000 description 1
- 239000003440 toxic substance Substances 0.000 description 1
- 231100000027 toxicology Toxicity 0.000 description 1
- 239000003053 toxin Substances 0.000 description 1
- 231100000765 toxin Toxicity 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
Images
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/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Human Resources & Organizations (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Tourism & Hospitality (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Business, Economics & Management (AREA)
- General Health & Medical Sciences (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Educational Administration (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Computing Systems (AREA)
- Software Systems (AREA)
- Quality & Reliability (AREA)
- Evolutionary Biology (AREA)
- Crystallography & Structural Chemistry (AREA)
- Game Theory and Decision Science (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
Abstract
本发明公开一种多环芳烃综合生态风险评价和敏感性分析优化方法，以层次阿基米德Copula函数为核心，计算在考虑各PAHs之间的内部相关性前提下的概率，并与可忽略浓度与最大允许浓度标准相结合，对区域内PAHs的综合生态风险进行评价。本发明将层次阿基米德Copula函数与单一生态风险评价指标结合在一起，从统计角度系统评价多种PAHs的混合生态风险，较好地解决了由于内部相关性导致的综合风险与单一风险累加结果不一致的缺陷，定量、定性地评价综合生态风险，并依据其对不同因素的敏感性分析结果，分析降低风险的有效举措，具有客观性和合理性。
Description
技术领域
本发明属于多环芳烃(PAHs)生态风险评价，具体涉及一种多环芳烃综合生态风险评价和敏感性分析优化方法，尤其涉及一种基于层次阿基米德Copula函数和单一生态风险评价指标耦合的多环芳烃(PAHs)综合生态风险评价和敏感性分析的优化方法。
背景技术
多环芳烃(PAHs)是一类具有半挥发性、持久性的有机污染物，具有强烈的毒性、致癌性、致畸性与诱变性，常产生于热源或石油源，且广泛存在于环境中，危害动植物体的生命安全。据其对生态环境的危害性，美国环保署将16种PAHs列如优先控制污染物名单，其中7种被国际癌症研究机构认定对人类有潜在致癌作用。分析多环芳烃的来源及风险程度对于研究其对生态环境的影响有着重要意义。生态风险评价是指评估或评价一个生态系统中特定事件发生可能性的过程。为了评估生态环境中多环芳烃的单一和整体危害，近年来出现了多种针对多环芳烃等有毒物质生态风险的评价方法，并得到了广泛应用。应用最为广泛评估方法包括风险商法(RQ)、苯并[a]芘毒性当量法(TEQBaP)、概率风险评价法(PRA)等。
风险商法(Risk quotients,RQ)为最早使用的单一化合物生态风险评价方法，将暴露浓度与毒性数据标准浓度相比较得到的商值作为评价指标。苯并[a]芘毒性当量法(total toxic benzo(a)pyrene equivalent,TEQBaP)以毒性当量因子(toxic equivalencyfactors,TEFs)为基础，将浓度与当量因子相乘，得到以BaP为基数的毒性浓度。最大允许浓度(Maximum Permissible Concentrations,MPCs)、可忽略浓度(NegligibleConcentrations,NCs)等可以用于单一或多个化合物的综合生态风险评价。概率风险评价法(Probabilistic Risk Assessment,PRA)计算暴露浓度曲线与毒性数据曲线之间重叠区域的面积以及联合概率，评价基于一系列潜在风险和可能影响的生态风险，为风险商评价提供了许多可靠的决策依据。在进行有毒污染物联合生态风险评价时，一般采取将不同污染物进行简单加和或按照毒性当量因子加权评价的方法。但由于多环芳烃之间协同作用或拮抗作用的存在，对16种多环芳烃的综合风险评估可能不适合通过求和或加权计算。
发明内容
发明目的：针对上述现有技术的不足，本发明的目的在于提供一种多环芳烃综合生态风险评价和敏感性分析优化方法。
技术方案：一种多环芳烃综合生态风险评价和敏感性分析优化方法，所述方法基于层次阿基米德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)对于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可以表示为：
式中，
其中Rij为Xij在(X1j,X2j,…,Xnj)的秩；k为变量个数；n为每个变量的观测值个数。
进一步的，步骤(3)基于层次阿基米德Copula函数对其进行综合概率模型构建过程如下：
层次阿基米德Copula函数表达式如下：
式中，H(·)表示HAC函数；ui为第i个变量的边际分布函数；生成因子
为[0,1]至[0,∞]的映射，为连续单调递增的凸函数。
更进一步的，步骤(3)对各个分组PAHs构建第一层层次阿基米德Copula模型，基于递归极大似然法对模型参与进行计算，选择拟合效果最优的层级阿基米德Copula函数类型进行分层建模，得到PAHs暴露浓度综合概率模型；以n维层次阿基米德Copula模型为例，拟合参数优选方法如下：
θ值越大，应用Kendall相关系数表示的依赖性就越强。按照此方法计算后续参数，并进行拟合度检验，最终可得第一层阿基米德Copula函数类型及完整的模型结构；
根据第一层模型的概率分布结果，重复上述步骤构建第二层层次阿基米德Copula模型，得到层次阿基米德Copula函数建立暴露浓度综合概率模型。
进一步的，步骤(4)选择PAHs单一生态指标值代入PAHs暴露浓度综合概率模型，计算综合生态标准值及PAHs综合生态风险程度，计算公式如下：
其中，E表示PAHs综合生态风险程度，表示综合浓度与整体沉积物质量标准之间的关系，CCPAHs表示PAHs综合浓度值，CCNCs表示综合可忽略浓度标准值，CCMPCs表示综合最大允许浓度标准值。
更进一步的，步骤(4)包括根据PAHs综合生态风险程度E的大小将PAHs的综合生态风险分为四个等级，具体如下：
(a)E≤0，表示其综合浓度的风险低于综合的可忽略浓度对应的风险值，称为无风险等级；
(b)0＜E≤0.5，表示PAHs综合浓度在可忽略的浓度和最大允许浓度之间，称为低风险水平；
(c)0.5＜E≤1，表示PAHs综合浓度在可忽略的浓度和最大允许浓度之间，称为中风险水平；
(d)E＞1，表示PAHs综合浓度风险超出综合最大允许浓度对应风险值，称为高风险等级。
有益效果：本发明为抵消在计算综合生态风险时不同PAHs之间的协同或拮抗作用，以可忽略浓度与最大允许浓度为标准值，提出了一种基于两层的层次阿基米德Copula函数的综合生态风险评价模型，并据其进行敏感性分析，以分析减弱综合风险的针对性措施。层次阿基米德Copula函数可以实现对数据集中相关性的捕捉，引入沉积物质量标准以评价其综合风险。与现有技术下相比，具有以下显著特点：
(1)在考虑变量间相关性的前提下进行定性及定量的综合风险评价，且优化了模型计算速度与拟合准确性。本发明在评估综合生态风险时考虑了PAHs的内在相关性，这是与其他直接累加或加权单一PAH风险得到综合生态风险的评估方法(如RQ方法)相比最核心优势。本发明对一层模型的选择机制进行了优化，直接采用层级阿基米德Copula函数拟合变量概率分布，并应有递归极大似然法对进行参数计算，在保证评估结果合理性的基础上，优化了计算程序，提高了计算速度和函数拟合的准确性，更有助于定性及定量地评价PAHs综合生态风险。
(2)基于层次阿基米德Copula评价模型进行敏感性分析，为寻找最优化的风险消减措施提供帮助。本发明提供了机会来分辨具有高生态风险的PAHs的环数，并进一步获得了迫切需要关注或修复的PAHs源头的类型。结合研究区域的综合生态风险概况，基于风险管理目标规划最有效的短期及长期生态风险消减策略。
附图说明
图1为本发明的整理流程图；
图2为4维层次阿基米德Copula函数结构示意图；
图3为实施例中PAHs暴露浓度分析图；
图4为实施例中PAHs综合生态风险评价结果图；
图5为实施例中PAHs分组生态风险评价结果图；
图6为实施例中敏感性分析结果图；
图7为实施例中综合风险与各种因素的相关性示意图。
具体实施方式
为了详细的说明本发明所公开的技术方案，下面结合说明书附图及具体实施例对本发明的技术方案进行详细说明。
本发明公开的是一种多环芳烃(PAHs)综合生态风险评价和敏感性分析的优化方法，基于层次阿基米德Copula函数和单一生态风险评价指标耦合的多环芳烃(PAHs)综合生态风险评价和敏感性分析的优化模型。为便于理解本发明，做以下说明：
首先，本发明所述方法中多变量Copula函数可用下式表示：
其中，
为生成因子，在[0,1]范围内连续单调递减。
用于函数类型筛选的三种阿基米德Copula函数归纳如下：
Gumbel Copula：
Clayton Copula：
Frank Copula：
如图1所示，本发明基于层次阿基米德Copula函数和单一生态风险评价指标耦合的PAHs综合生态风险评价和敏感性分析的优化模型，具体包含以下步骤：
1、获取PAHs样本浓度数据
采用按照常规检测方法提取测定的多组同一研究水域范围内表层沉积物样本中的PAHs暴露浓度数据，并记录样本编号及取样点位置。PAHs种类选取依照美国环保署按PAHs毒性、致癌性等指标划定的16种优先控制污染物范围，包括：萘(Naphthalene,Nap)，苊(acenaphthene,Ace)，苊烯(Acenaphthylene,Acy)，芴(Fluorine,Flu)，菲(Phenanthrene,Phe)，蒽(Anthracene,Ant)，荧蒽(Fluoranthene,Flua)，芘(Pyrene,Pyr)，苯并[a]蒽(Benz[a]anthracene,BaA)，
(Chrysene,Chr)，苯并[b]荧蒽(Benzo[b]fluoranthene,BbF)，苯并[k]荧蒽(Benzo[k]fluoranthene,BkF)，苯并[a]芘(benzo[a]pyrene,BaP)，二苯并[a,h]蒽(Dibenzo[a,h]anthracene,DBA)，茚并[1,2,3-cd]芘(Indeno[1,2,3-c,d]pyrene,IP)以及苯并[g,h,i]芘(Benzo[g,h,i]pyrene,BghiP)。
2、PAHs暴露浓度一致性度量
应用一致性度量研究变量之间的相互依赖关系，应选择合适的方法对每两种或多种PAHs暴露浓度进行相关性计算。对于两变量间的一致性度量，若浓度变量均符合正态分布假定，采用Pearson相关系数计算；若浓度变量为非正态分布数据，应用Spearman秩相关系数计算。将PAHs按照苯环数目不同分为四组：2-ring，3-ring，4-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可以表示为：
3、基于层次阿基米德Copula函数建立暴露浓度综合概率模型
若数据一致性度量结果显示暴露浓度具有独立性，即各变量间不存在相关性，则可应用常规生态风险评价方法，由单一风险的累加或加权计算得到综合风险评价结果。层次阿基米德Copula函数可以通过对多个阿基米德Copula函数的嵌套以实现对数据集中相关性的捕捉。若结果表明存在相关性或相关性显著，则需结合层次阿基米德Copula函数对其进行综合概率模型构建，以保存其内部相关性，并作为综合生态风险评价的基础。
一般Copula函数可以表示为：
F(x1,x2,…,xn)＝C(F1(x1),F2(x2),…,Fn(xn)；θ)＝C(u1,u2,…,un)
其中，Fi(xi)与ui表示第i个变量的边际分布函数，Fi(xi)＝P[Xi≤xi]，即当Xi≤xi时的随机变量的概率；C(·)表示对应的Copula函数；θ为函数参数。
对于4维及以上的层次阿基米德Copula函数，如图2所示，其包含完全嵌套与部分嵌套两种类型，以4变量为例，完全嵌套公式如下：
部分嵌套公式如下：
其中，H(·)表示HAC函数,ui为第i个变量的边际分布函数，生成因子
为[0,1]至[0,∞]的映射，为连续单调递增的凸函数。其中，图2(a)为完全嵌套4维层次阿基米德Copula函数，图2(b)为部分嵌套4维层次阿基米德Copula函数。
基于递归极大似然法对模型参与进行计算，选择拟合效果最优的层级阿基米德Copula函数类型进行分层建模，即可得到PAHs暴露浓度综合概率模型。具体步骤如下：
第一步：PAHs按照苯环数目不同分为四组：2-ring，3-ring，4-ring及5-ring(含苯环数目≥5)；
第二步：对各个分组PAHs构建第一层层次阿基米德Copula模型，采用递归极大似然法拟合参数。以4层层次阿基米德Copula模型为例，θ1计算公式如下：
θ值越大，对于大多数阿基米德Copula来说，应用Kendall相关系数表示的依赖性就越强。取θ1＝θu1u2，可得H1(u1,u2；θu1u2)。按照此方法计算后续参数，并进行拟合度检验，最终可得到拟合度最优的阿基米德Copula函数类型及完整的模型结构。
第三步：根据第一层模型的概率分布结果，重复步骤(3.2)构建第二层层次阿基米德Copula模型，即可得到层次阿基米德Copula函数建立暴露浓度综合概率模型。
4、选取适宜的生态风险指标，评价综合生态风险
根据生态环境质量要求，选择适合的PAHs单一生态指标值(以可忽略浓度NCs与最大允许浓度MPCs为例)，将其代入已建好的PAHs综合概率模型，计算综合生态标准值及PAHs综合生态风险程度，公式如下：
其中，E表示PAHs综合生态风险程度，即综合浓度与整体沉积物质量标准之间的关系，CCPAHs表示PAHs综合浓度值，CCNCs表示综合可忽略浓度标准值，CCMPCs表示综合最大允许浓度标准值。
根据PAHs综合生态风险程度计算结果对研究区域PAHs的综合生态风险进行评价。根据PAHs综合生态风险程度E的大小，PAHs的综合生态风险可分为四个等级：
(4.1)E≤0，表明其综合浓度的风险低于综合的可忽略浓度对应的风险值，可以认为是一个安全的环境，被称为无风险等级；
(4.2)0＜E≤1，表明PAHs综合浓度在可忽略的浓度和最大允许浓度之间，以中间值0.5为界，可将其划分为0＜E≤0.5和0.5＜E≤1两段，分别表示低风险水平和中风险水平；
(4.3)E＞1表明PAHs综合浓度风险超出综合最大允许浓度对应风险值，即生态系统面临较高的生态风险，需要采取紧急补救措施以防止风险恶化。
E的数值也可以量化每个级别的风险。每个级别上不同样本的风险值的大小可用于直接比较，并且两者之间的程度代表其生态风险差距。当结果的值为负数且E的绝对值较大甚至无限大时，表明样品的总体风险较小且环境更安全。而当E大于1时，值越大表示风险高于安全范围的距离越远。
5、敏感性分析辅助制定减弱风险决策
为制定减弱综合生态风险的针对性措施，探讨不同种类的PAHs对综合生态风险概率的趋势和程度，对PAHs综合生态风险采用降低固定步长分组风险概率的方法进行敏感性分析，计算步骤如下：
E(f(x′))＝E(f(x|(xi,j-t))),(i＝1,2,…,N；j＝1,2,…,m)
其中，E(f(x′))表示受试PAH的联合生态风险程度；xi,j表示敏感因子；t表示概率分布变化的步长；i表示样本数量，共计N个样本；j表示分组的编号，PAHs被划分为m个组。
采用如下公式测量综合生态风险对每组风险概率变化的敏感性：
S(f(x′))＝E(f(x))-E(f(x′))
根据敏感性分析结果可得到消减综合生态风险针对性建议。
实施例1
本实施以太湖表层沉积物样本作为实际应用
以长三角重要饮用水源地太湖为例，以2010年6月提取的28组表层沉积物中16种PAHs暴露浓度为样本，用基于层次阿基米德Copula函数和单一生态风险评价指标耦合的PAHs综合生态风险评价和敏感性分析的优化模型对该研究区域进行生态风险评价并提出风险消减措施。
(1)研究区域概况
本文的数据资料来源于太湖表层沉积物样本检测，2010年6月28组沉积物样本中16种PAHs的暴露浓度。太湖作为中国第三大淡水湖泊，位于中国长江三角洲地区，对包括上海、苏州、无锡等几个重要城市提供饮用水供应，在防洪及工农业用水方面发挥着极其重要的作用。由于周围人类的频繁活动和作为饮用水源的特殊性，太湖的水质情况值得长期关注。样本中，16种PAHs的总浓度变化范围为179.5～1669.3ng/g dw，平均值为533.3ng/gdw。各样本点的PAHs暴露浓度情况见图3。
(2)模型构建
首先对16种PAHs暴露浓度数据按照苯环数目不同分为4组，分别为2-ring，3-ring，4-ring及5-ring(含苯环数目≥5)。对所有PAHs进行两变量一致性度量并分组计算多变量Kendall秩相关性。由于暴露浓度不符合正态分布，故选用Spearman秩相关系数计算两变量相关性。计算结果显示，两个PAH之间的Spearman相关系数在0.5到0.9之间，每组的Kendall多变量相关系数均不小于0.75，这与应用本专利模型进行风险评价的条件一致。
按照前述步骤，对PAHs暴露浓度构建两层层次阿基米德Copula模型并根据拟合效果选择最优阿基米德Copula函数，模型构建结果见表1。
表1 模型结构及参数
(3)生态风险评价
通过两层嵌套，可以分别对PAHs暴露浓度的分组及综合生态风险进行评价。图4显示了基于本专利模型的PAHs综合生态风险评价结果，图5显示了基于本专利模型的分组风险评价结果。
图4中可得本专利模型评估的综合生态风险结论为：实施例研究区域内综合生态风险程度为0.0011(S27)至0.9110(S7)之间。其中，S7仍然为具有最高综合生态风险的样本点，其次为S6；而S4，S18，S21和S27的风险在所有样本中相对居于较低水平。对于风险等级的判断，S3、S5-S9、S13和S19划分为中风险类别，其他均属于低风险级别。
图5中可得每组的风险概况，这有助于识别PAHs苯环的风险种类。总风险被划分为相对较高的中风险的样本，其分组风险也远高于其他样本。与其他几种类别相比，5-ringPAHs是大多数样品中风险最低的类别，甚至大多数样品的浓度低于毒理学可忽略不计的标准。相反，在S5和S6中，5-ring PAHs几乎占据了最高风险位置。
(4)敏感性分析
图6展示了综合生态风险随敏感因子变化情况，图7为敏感性分析结果与多种因素间的相关性研究。其中，图7(a)为综合生态风险与各样本点风险变化程度比较图；图7(b)为最大风险降低程度与综合生态风险对数关系图；图7(c)为各组PAHs风险降低程度与最大风险相关散点图；图7(d)为风险降低程度与相应风险散点图；图7(e)为综合生态风险首降点与综合生态风险关系图；图7(f)为综合生态风险首降点与最大分组风险关系图。
对于风险最高的样本S7，在敏感度因子的变化范围内，5-ring PAHs风险的变化首先对总风险产生影响，随后是2-ring，但后者引起的风险消减程度略大于前者；而3-ring和4-ring风险降低不会对总风险产生任何影响，表明它们的风险需要降低更大幅度才会对综合风险等级产生影响。对于风险较低的S4和S27等，部分敏感因子的小幅度降低即可导致综合风险直接降低至负值，即达到无风险水平，但也存在一些敏感因子需要大幅度降低才能影响总风险，这可能与整体风险的主要控制因素以及各敏感因子之间的相关性有关。图7可得，在计算出的敏感性因子范围内，综合风险的变化范围随风险对数的增加而减小，且综合风险越低，其对因子变化的敏感性越高。对于单个因素的变化，综合风险开始变化的位置与初始风险没有显著关系，且在计算范围内，综合生态风险最敏感的因素并不总是能够使总体风险降低最多的因素。
(5)提出风险消减措施
敏感性分析仍可用于支持针对特定样本区域制定生态风险降低策略，以找到最有效的风险管理方法。以S5为例，降低4环PAHs的风险概率可以迅速将其综合风险由中风险等级迅速降低至低风险等级，即专注于控制此类PAHs的排放可以有效地降低其风险。由此，可以根据周围环境的实际情况，选择加强对汽车尾气排放标准的控制，限制草、煤和木材等化石燃料的燃烧，防止森林火灾等适宜的防治措施，以迅速降低生态环境中PAHs的风险等级。
通过上述实施例可以看出，本发明在考虑到PAHs暴露浓度中的内在相关性的基础上，基于层次阿基米德Copula函数与沉积物质量标准建模，对其综合生态风险进行评价，具有合理性和有效性，并最终借助敏感性分析寻找降低综合风险最有效的敏感因子，为制定具有针对性的风险消减措施提供决策支持。
综上所示，本发明将层次阿基米德Copula函数与PAHs沉积物质量标准相结合，在考虑变量间的相关性前提下，以较快的模型计算速度实现对PAHs生态风险的定性及定量评价，并可基于敏感性分析得到具有针对性的短期及长期生态风险消减建议，具有合理性和有效性。
Claims (6)
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可以表示为：
式中，
其中Rij为Xij在(X1j，X2j，...，Xnj)的秩；k为变量个数；n为每个变量的观测值个数。
3.根据权利要求1所述的多环芳烃综合生态风险评价和敏感性分析优化方法，其特征在于：步骤(3)基于层次阿基米德Copula函数对其进行综合概率模型构建过程如下：
层次阿基米德Copula函数表达式如下：
式中，H(·)表示HAC函数；ui为第i个变量的边际分布函数；生成因子
为[0，1]至[0，∞]的映射，为连续单调递增的凸函数。
4.根据权利要求1所述的多环芳烃综合生态风险评价和敏感性分析优化方法，其特征在于：步骤(3)对各个分组PAHs构建第一层层次阿基米德Copula模型，基于递归极大似然法对模型参与进行计算，选择拟合效果最优的层级阿基米德Copula函数类型进行分层建模，得到PAHs暴露浓度综合概率模型；以n维层次阿基米德Copula模型为例，拟合参数优选方法如下：
θ值越大，应用Kendall相关系数表示的依赖性就越强；
按照此上述方法计算后续参数，并进行拟合度检验，最终可得第一层阿基米德Copula函数类型及完整的模型结构；根据第一层模型的概率分布结果，重复上述步骤，构建第二层层次阿基米德Copula模型，得到层次阿基米德Copula函数建立暴露浓度综合概率模型。
5.根据权利要求1所述的多环芳烃综合生态风险评价和敏感性分析优化方法，其特征在于：步骤(4)选择PAHs单一生态指标值代入PAHs暴露浓度综合概率模型，计算综合生态标准值及PAHs综合生态风险程度，计算公式如下：
其中，E表示PAHs综合生态风险程度，表示综合浓度与整体沉积物质量标准之间的关系，CCPAHs表示PAHs综合浓度值，CCNcs表示综合可忽略浓度标准值，CCMPCs表示综合最大允许浓度标准值。
6.根据权利要求1所述的多环芳烃综合生态风险评价和敏感性分析优化方法，其特征在于：步骤(4)包括根据PAHs综合生态风险程度E的大小将PAHs的综合生态风险分为四个等级，具体如下：
(a)E≤0，表示其综合浓度的风险低于综合的可忽略浓度对应的风险值，称为无风险等级；
(b)0＜E≤0.5，表示PAHs综合浓度在可忽略的浓度和最大允许浓度之间，称为低风险水平；
(c)0.5＜E≤1，表示PAHs综合浓度在可忽略的浓度和最大允许浓度之间，称为中风险水平；
(d)E＞1，表示PAHs综合浓度风险超出综合最大允许浓度对应风险值，称为高风险等级。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010138645.7A CN111445109B (zh) | 2020-03-03 | 2020-03-03 | 一种多环芳烃综合生态风险评价和敏感性分析优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010138645.7A CN111445109B (zh) | 2020-03-03 | 2020-03-03 | 一种多环芳烃综合生态风险评价和敏感性分析优化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111445109A true CN111445109A (zh) | 2020-07-24 |
| CN111445109B CN111445109B (zh) | 2023-08-01 |
Family
ID=71655757
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010138645.7A Active CN111445109B (zh) | 2020-03-03 | 2020-03-03 | 一种多环芳烃综合生态风险评价和敏感性分析优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111445109B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112098615A (zh) * | 2020-08-27 | 2020-12-18 | 中国水产科学研究院东海水产研究所 | 一种贝类水产品中多环芳烃限量值评估方法 |
| CN117725685A (zh) * | 2024-02-05 | 2024-03-19 | 中汽研汽车检验中心（天津）有限公司 | 一种车辆操控性多目标优化方法和设备 |
| CN117725685B (zh) * | 2024-02-05 | 2024-05-31 | 中汽研汽车检验中心（天津）有限公司 | 一种车辆操控性多目标优化方法和设备 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120318982A1 (en) * | 2009-11-02 | 2012-12-20 | Leslie Joseph Janik | Method for determining polycyclic aromatic hydrocarbon contaminant concentration |
| CN103853928A (zh) * | 2014-03-14 | 2014-06-11 | 金小伟 | 一种基于多等级综合评估的水环境生态风险评估模型 |
| US20160025698A1 (en) * | 2013-12-24 | 2016-01-28 | Jiangsu Provincial Academy Of Environmental Science | Method for determining ecological risks of polycyclic aromatic hydrocarbon in water body |
-
2020
- 2020-03-03 CN CN202010138645.7A patent/CN111445109B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120318982A1 (en) * | 2009-11-02 | 2012-12-20 | Leslie Joseph Janik | Method for determining polycyclic aromatic hydrocarbon contaminant concentration |
| US20160025698A1 (en) * | 2013-12-24 | 2016-01-28 | Jiangsu Provincial Academy Of Environmental Science | Method for determining ecological risks of polycyclic aromatic hydrocarbon in water body |
| CN103853928A (zh) * | 2014-03-14 | 2014-06-11 | 金小伟 | 一种基于多等级综合评估的水环境生态风险评估模型 |
Non-Patent Citations (1)
| Title |
|---|
| YUELI YAO, DONG WANG, HAITAO MA: "The impact on T-regulatory cell related immune responses in rural women exposed to polycyclic aromatic hydrocarbons (PAHs) in household air pollution in Gansu, China: A pilot investigation", pages 306 - 317 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112098615A (zh) * | 2020-08-27 | 2020-12-18 | 中国水产科学研究院东海水产研究所 | 一种贝类水产品中多环芳烃限量值评估方法 |
| CN117725685A (zh) * | 2024-02-05 | 2024-03-19 | 中汽研汽车检验中心（天津）有限公司 | 一种车辆操控性多目标优化方法和设备 |
| CN117725685B (zh) * | 2024-02-05 | 2024-05-31 | 中汽研汽车检验中心（天津）有限公司 | 一种车辆操控性多目标优化方法和设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111445109B (zh) | 2023-08-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Macedo et al. | Development of a benthic macroinvertebrate multimetric index (MMI) for Neotropical Savanna headwater streams | |
| Wright et al. | River classification using invertebrates: RIVPACS applications | |
| Norris et al. | Monitoring river health | |
| Resh et al. | Design and implementation of rapid assessment approaches for water resource monitoring using benthic macroinvertebrates | |
| Roset et al. | A review of existing fish assemblage indicators and methodologies | |
| Wright | An introduction to RIVPACS | |
| Einheuser et al. | Modeling the effects of conservation practices on stream health | |
| Raven et al. | Towards an integrated approach to classifying and evaluating rivers in the UK | |
| Wright et al. | Assessing the biological quality of fresh waters | |
| CN106447205A (zh) | 一种基于层次分析法的配电自动化终端状态评价方法 | |
| CN117195135B (zh) | 一种水污染异常溯源检测方法及系统 | |
| Whittier et al. | The surface waters component of the Environmental Monitoring and Assessment Program (EMAP): an overview | |
| CN111445109A (zh) | 一种多环芳烃综合生态风险评价和敏感性分析优化方法 | |
| Jones et al. | Diatoms as indicators of fine sediment stress | |
| Gazendam et al. | Evaluation of the Qualitative Habitat Evaluation Index as a planning and design tool for restoration of rural Ontario waterways | |
| Baltacı et al. | Water quality monitoring studies of Turkey with present and probable future constraints and opportunities | |
| CN101527016A (zh) | 出入境检验检疫符合性条件的筛选方法 | |
| Clairain Jr | Hydrogeomorphic approach to assessing wetland functions: guidelines for developing regional guidebooks. Chapter 1, Introduction and overview of the hydrogeomorphic approach | |
| CN114037202A (zh) | 多指标关联的河流水质评价方法、系统、设备及存储介质 | |
| Tsai et al. | Managing multinational sustainable development in the European Union based on the DPSIR framework | |
| Marsh et al. | Guidance on ecological responses and hydrological modelling for low-flow water planning | |
| CN108364090B (zh) | 一种预测城市道路表面苯系物累积负荷的方法 | |
| Brown et al. | Use of expert opinion elicitation to quantify the internal erosion process in dams | |
| Gowan et al. | Stream monitoring methods suitable for citizen volunteers working in the coastal plain and lower Piedmont regions of Virginia | |
| CN117009767B (zh) | 一种基于生物有效性的土壤基准制定及风险评估的方法 |
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 |