CN109190884A - 一种河道水体综合评价方法 - Google Patents
一种河道水体综合评价方法 Download PDFInfo
- Publication number
- CN109190884A CN109190884A CN201810833611.2A CN201810833611A CN109190884A CN 109190884 A CN109190884 A CN 109190884A CN 201810833611 A CN201810833611 A CN 201810833611A CN 109190884 A CN109190884 A CN 109190884A
- Authority
- CN
- China
- Prior art keywords
- level
- index
- degree
- matrix
- weight
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 56
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 47
- 238000011156 evaluation Methods 0.000 claims abstract description 88
- 239000011159 matrix material Substances 0.000 claims abstract description 71
- 239000013598 vector Substances 0.000 claims abstract description 32
- 238000004458 analytical method Methods 0.000 claims abstract description 27
- 238000010606 normalization Methods 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims abstract description 3
- 239000010410 layer Substances 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 10
- 238000005314 correlation function Methods 0.000 claims description 8
- 239000002356 single layer Substances 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 230000006872 improvement Effects 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 238000011161 development Methods 0.000 claims description 4
- 238000003786 synthesis reaction Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000012851 eutrophication Methods 0.000 claims description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 2
- 238000012790 confirmation Methods 0.000 claims description 2
- 229910052760 oxygen Inorganic materials 0.000 claims description 2
- 239000001301 oxygen Substances 0.000 claims description 2
- 238000012163 sequencing technique Methods 0.000 claims description 2
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims 4
- 229910052757 nitrogen Inorganic materials 0.000 claims 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims 1
- MMDJDBSEMBIJBB-UHFFFAOYSA-N [O-][N+]([O-])=O.[O-][N+]([O-])=O.[O-][N+]([O-])=O.[NH6+3] Chemical group [O-][N+]([O-])=O.[O-][N+]([O-])=O.[O-][N+]([O-])=O.[NH6+3] MMDJDBSEMBIJBB-UHFFFAOYSA-N 0.000 claims 1
- 229930002868 chlorophyll a Natural products 0.000 claims 1
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 claims 1
- 238000007689 inspection Methods 0.000 claims 1
- 230000033116 oxidation-reduction process Effects 0.000 claims 1
- 229910052698 phosphorus Inorganic materials 0.000 claims 1
- 239000011574 phosphorus Substances 0.000 claims 1
- 238000012360 testing method Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 description 7
- 238000013278 delphi method Methods 0.000 description 4
- 230000007613 environmental effect Effects 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000008439 repair process Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- PWHULOQIROXLJO-UHFFFAOYSA-N Manganese Chemical compound [Mn] PWHULOQIROXLJO-UHFFFAOYSA-N 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000010866 blackwater Substances 0.000 description 1
- 239000003922 charged colloid Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000003651 drinking water Substances 0.000 description 1
- 235000020188 drinking water Nutrition 0.000 description 1
- 238000011985 exploratory data analysis Methods 0.000 description 1
- 238000000855 fermentation Methods 0.000 description 1
- 230000004151 fermentation Effects 0.000 description 1
- 229910001385 heavy metal Inorganic materials 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 229910052748 manganese Inorganic materials 0.000 description 1
- 239000011572 manganese Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000002906 microbiologic effect Effects 0.000 description 1
- 238000002715 modification method Methods 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 239000002957 persistent organic pollutant Substances 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000007115 recruitment Effects 0.000 description 1
- 238000009418 renovation Methods 0.000 description 1
- 239000010865 sewage Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 241001148471 unidentified anaerobic bacterium Species 0.000 description 1
- 239000002351 wastewater Substances 0.000 description 1
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
-
- 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
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种河道水体综合评价方法，包括以下步骤：步骤S1、建立三级指标评价体系，所述三级指标评价体系包括一级指标、二级指标和三级指标，所述三级指标包括多个等级；步骤S2、通过物元分析法对河道水体的黑臭程度进行量化评价，包括建立三级指标的指标评价矩阵、对指标评价矩阵和指标评价矩阵所对应的指标向量进行归一化处理和根据归一化处理后的指标评价矩阵和指标向量确定关联度；步骤S3、通过层次‑德尔菲分析法确定各项所述三级指标的权重，并通过设置专家评判表、建立专家评价矩阵、计算评价向量夹角，并使用归一化后的专家评判一致度对权重修正。
Description
技术领域
本发明涉及环保领域，尤其涉及一种河道水体综合评价方法。
背景技术
随着我国经济的快速发展与人们生活水平的提高，城市景观河道水体面临着工业废水排放、生活污水汇入、城镇面源污染加剧、生态功能丧失等挑战，国内较多城市出现了河道黑臭现象，给城市发展及居民生活质量带来了众多不利影响。在我国全面深化改革、开启生态文明建设进程中，我国也开始了城市黑臭河道开展综合整治，在开展有效治理与改善措施基础上，如何评价城市河道黑臭水体在整治前后的效果以及长效运行效果，对实现城市水环境可持续提升显得尤为重要。
根据住房城乡建设部和环境保护部联合发布的《城市黑臭水体整治工作指南》，城市黑臭水体是指城市建成区内呈现令人不悦的颜色和(或)散发令人不适气味的水体的统称。黑臭河道的形成主要是由于大量有机污染物进入水体后，好氧微生物的生化作用消耗了水体中大量的O2，导致厌氧细菌大量繁殖，使有机物腐败、分解、发酵，分解为NH3-N、腐殖质、H2S等有臭气体逸出水面进入大气，同时水中铁、锰等重金属被还原，与水中的硫离子形成了大量带负电胶体的物质悬浮颗粒，使水体变黑、变臭。然而目前针对黑臭河道的评价，大多数采用比较简单的水质指标评价，并不能全面反映黑臭河道治理前后的综合效果以及长效运行效果，且与《地表水环境质量标准(GB 3838-2002)》中规定的水质指标缺少联系。
关于评价问题目前的分析方法有很多，诸如聚类分析、典范对应分析法、蒙特卡洛模拟分析法、熵权综合指数法、模糊综合评价法、灰色关联度分析法等。但以上方法针对城市河道黑臭水平的评价相关程度不够契合：聚类分析是探索性分析方法，适用于数据挖掘；典范对应分析法虽然将样方数据与环境因子进行矩阵结合分析关联度，但是原始方法存在“弓形效应”问题；蒙特卡洛模拟分析法是基于概率统计方法基础进行的系统模拟，对于城市黑臭问题的模拟缺乏合适的数学概率模型；熵权综合指数法是一种基于指标变异程度进行权重修正的方法，由于城市黑臭河道选取指标信息量固定，故而该方法修正作用不明显；模糊综合评价法计算较为复杂，且当指标集个数较多时存在分辨率低的问题；灰色关联度分析法多用于分析不同因素对系统随时间影响与关联度的判断。
发明内容
本发明所要解决的技术问题是：提供一种河道水体综合评价方法，提升河道水体综合评价的准确度。
为了解决上述技术问题，本发明采用的技术方案为：
一种河道水体综合评价方法，其特征在于：包括以下步骤：步骤S1、建立三级指标评价体系，所述三级指标评价体系包括一级指标、二级指标和三级指标，所述三级指标包括多个等级；
步骤S2、通过物元分析法对河道水体的黑臭程度进行量化评价，包括建立三级指标的指标评价矩阵、对指标评价矩阵和指标评价矩阵所对应的指标向量进行归一化处理、以及确定归一化处理后的指标评价矩阵与所述三级指标的各个等级之间的关联度；
步骤S3、通过层次-德尔菲分析法确定关联度中各项三级指标的权重，接收多个针对所述各项三级指标的权重判断矩阵，计算各个权重判断矩阵的平均一致度并对所有平均一致度进行归一化处理得到相对一致度，使用所述相对一致度对三级指标的权重进行修正。
本发明的有益效果在于：传统的河道生态评价方法中缺乏合适的指标评价权重的修正方法，导致经过层次分析后的结果受主观因素影响较多，本发明采用德尔菲法进行层次分析后的权重修正，其本质上是一种反馈匿名函数法，接收多个针对所述各项三级指标的权重判断矩阵，经归纳、统计后，以简便、直观的方式进行评判比较，相比于传统直接的层次分析法，加入德尔菲修正的结果更加客观公正，对于真实的指标对目标城市河道的黑臭影响评价更加准确，使该方法对于河道的适用范围也更广。
附图说明
图1所示为本发明的河道水体综合评价方法流程图；
图2所示为本发明的河道水体综合评价指数指标体系构建图；
具体实施方式
为详细说明本发明的技术内容、所实现目的及效果，以下结合实施方式并配合附图予以说明。
本发明最关键的构思在于：引入德尔菲法作为权重修正手段，形成了更为准确的综合评价方法，接收多个针对所述各项三级指标的权重判断矩阵，之后计算各个权重判断矩阵的相对一致度并使用相对一致度对三级指标的权重进行修正，使得整体评价方法的准确度更高，修正的结果更加客观公正，适用的评价指标更为广泛的特性。
请参照图1和图2所示，本发明的一种河道水体综合评价方法，包括多指标多层次评价方法的建立和相关权重的赋予和修正，评价目标河道的综合黑臭程度，给城市黑臭河道的综合整治方案的提出以及治理后水质改善效果评估提供参考和建议。
步骤S1、建立评价体系：基于城市内河的生态特征等特点，明确评价角度与评价指标，建立三级指标评价体系，对河道黑臭水平、富营养化水平、生物多样性水平、水动力学水平等方面进行评价，城市河道综合评价指标体系及其解释如表1所示。
表1城市河道综合评价指标体系及其解释
进一步的，确定三级指标的评价等级：根据上述三级指标，参考《地表水环境质量标准(GB 3838-2002)》与《生活饮用水卫生标准(GB5749-2006)》，结合实际黑臭河道情况，将各指标划分成六级，并进行区间划分，如表2所示：
表2三级指标定量区间划分等级
步骤S2、通过物元分析法对河道水体的黑臭程度进行量化评价，包括建立三级指标的指标评价矩阵、对指标评价矩阵和指标评价矩阵所对应的指标向量进行归一化处理、以及确定归一化处理后的指标评价矩阵与所述三级指标的各个等级之间的关联度；具体如下：
将三级指标的指标评价矩阵设为Y'i，其中Y'i的三级指标对每一等级的评价值用L'ij(i＝1、2、……、17，j＝1、2、3、4、5、6)表示，i表示对应的第i个三级指标，j表示对应的一、二、三、四、五、六级。
进一步的，对指标评价矩阵和指标评价矩阵所对应的指标向量进行归一化处理的具体方法为：
对指标评价矩阵Y'i中的评价值进行归一化处理，使其评价值隶属于[0，1]的区间，若指标评价矩阵Y'i中三级指标的节域是严格的闭区间，将每个经典域的端点值除以节域的右端点值，若节域是[m，+∞)的半闭半开区间，则将m扩大四倍后的值设为节域右端点值，然后将每个经典域的端点值除以节域的右端点值，最后得到归一化后的指标评价矩阵Yi，其中每个三级指标的经典域用Lij表示；
对指标向量中的评价值的归一化处理，使其评价值隶属于[0，1]的区间，用b'分别除以Y'i中第i个三级指标归一化时所除以的数值，得到归一化后的指标向量B，其中B＝(b1、b2、…bi…b17)，bi(i＝1、2、……17)代表规范化处理后第i个三级指标的评价值，B代表规范化处理后的指标向量。
进一步的，所述经典域的物元矩阵表示方法为：
将三级指标的等级Q分为6个等级Q1、Q2、Q3、Q4、Q5、Q6，每个等级中各个指标En对应不同的量值范围，特征En的量值范围为Xpj＝[Apj，B pj](p＝1，2，3，4，5，6)，则经典域的物元矩阵Kp可表示为：
所述节域的物元矩阵表示方法为：
Q0为各个等级的Qp构成，则将K0称为节域物元，其量值范围为X0j＝[A0j，B 0j]，黑臭河道的综合黑臭程度的节域物元矩阵K0可表示为：
进一步的，所述关联度的确定方法为：
将关联度设为T，计算评价指标矩阵Gi关于等级j的关联度T：
上式中，W′为每个指标进过修正后的权重，T(x)为同一指标不同安全等级对应的关联函数值中最大的关联函数值，其中T(x)＝max(T(xj))；
根据归一化处理后的评价指标矩阵和指标评价向量，采用关联函数的定义，确定关联度矩阵T，关联矩阵由关联函数的每个值组成，为：
上式中：θ(X,X0)表示物元量值X到有限区间X0(X0＝[A，B])，A、B分别为有限元区间的最小值与最大值的距离；θ(X,XP)为待测物元量值X与有限区间Xp(Xp＝[Aj，Bj])的距离；X，X0，Xj分别为待评河道黑臭程度物元的量值、经典域物元的量值范围和节域物元的量值范围。点与有限元的间距计算公式为：
关联度T(x)数值的大小表示待评物元符合标准等级的隶属程度，其值越大，表明符合程度越高，关联度数值与评价等级关系如表3所示:
表3关联度数值与评价等级关系表
步骤S3、通过层次-德尔菲分析法确定关联度中各项三级指标的权重，接收多个针对所述各项三级指标的权重判断矩阵，计算各个权重判断矩阵的平均一致度并对所有平均一致度进行归一化处理得到相对一致度，使用所述相对一致度对三级指标的权重进行修正。
具体包括：
步骤S301、采用层次分析法构建判断矩阵，选取同一层级水平下的三级指标，两两比较评估三级指标的相对重要性，形成评估指标的比较矩阵F：
rij的取值采用九标度法获得，判断方式如表4所示
表4同层级水平下两两指标重要性对比表
| r<sub>i</sub>与r<sub>j</sub>重要性对比 | r<sub>ij</sub>值的选取 | r<sub>ji</sub>值的选取 |
| r<sub>i</sub>与r<sub>j</sub>同等重要 | 1 | 1 |
| r<sub>i</sub>比r<sub>j</sub>稍微重要 | 3 | 1/3 |
| r<sub>i</sub>比r<sub>j</sub>明显重要 | 5 | 1/5 |
| r<sub>i</sub>比r<sub>j</sub>强烈重要 | 7 | 1/7 |
| r<sub>i</sub>比r<sub>j</sub>极端重要 | 9 | 1/9 |
步骤S302、层次单排序与权重的确定，采用行和正规化方法，设n阶判断矩阵为：
按行求和计算有即得到权重向量W＝(v1、v2…、vn)T，经归一化使向量中各元素之和等于1后，得到向量W＝(w1、w2…、wn)T，W的元素即为同一层次因素对于上一层次因素某因素相对重要性的排序权值，次过程为层次单排序，其中，wi就是层次排序的优先程度，即各个指标的权重。
步骤S303、单层次排序的一致性检验，需要进行一致性检验确认层次单排序，一致性检验是指对A确定不一致的允许范围，定义CI为单层次排序的一致性指标，定义RI为单层次排序的随机一致性指标，RI取值如表5所示：
表5 RI值的选取
| n | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
| RI | 0 | 0 | 0.51 | 0.89 | 1.12 | 1.25 | 1.34 | 1.42 | 1.46 |
| n | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 |
| RI | 1.49 | 1.52 | 1.54 | 1.56 | 1.58 | 1.59 | 1.5943 | 1.6064 | 1.6133 |
定义一致性比率：CR＝CI/RI，当一致性比率CR＝CI/RI≤0.1时，认为A的不一致程度在容许范围之内，有满意的一致性，通过一致性检验，可用其归一化特征向量作为权向量，否则要重新构造成对比较矩阵A，对aij加以调整。
步骤S304、利用同一层次中所有层次单排序的结果，针对层次间所有因素，计算某一层次的层次总排序经归一化后，可以得到该层所有因素对于最高层相对重要性的权重值，假设上一层次所有因素A1、A2、…An，的总排序已经完成，的得到权重值分别为a1、a2、…an，与ai对应的本层次因素B1、B2、…Bn的单排序的结果为bi 1、bi 2、…bi n，如果Bj与Aj无关，则有bi j＝0。
总排序计算表6如下所示：
表6城市河道评价总排序计算表
步骤S305、层次总排序的一致性检验，设B层B1、B2、…、Bn对上层(A层)中因素的层次总排序一致性指标为CI，层次总排序随机一致性指标为RI，则层次总排序的一致性比率为：
当CR≤0.1时，认为层次总排序计算结果具有良好的一致性。
步骤S306、邀请多位专家参与决策，若由专家Z做出的权重判断矩阵形成的权重向量为：Wz＝(Wz1、Wz2、Wz3、…、Wzn)，(n＝1、2、…、17)，取权重向量Wi与Wj的夹角为Dij＝cos(θij)，Dij表示两专家各自的权重判断矩阵Wi与Wj之间的一致度，依据Dij求第i个专家的权重判断矩阵的平均一致度Di：
步骤S307、对所有专家的平均一致度进行归一化处理，得到专家的权重判断矩阵的相对一致度D'i，用归一化相对一致度D'i与相应权重向量相乘并求和，得修正后的权重：
进一步的，确定目标河道综合黑臭等级，根据物元分析城市河道评价指标关联度计算，结合层次分析和德尔菲修正后的各个生态要素的权重分布，进行关联度与权重乘积计算结果，即得到目标城市河道综合黑臭评价等级，根据2015年9月国家环境保护部发布的《黑臭水体治理技术政策》与三级评价指标六级划分情况及目标实际黑臭情况，对城市河道黑臭评价等级进行划分，如表7所示。
表7城市河道黑臭评价等级说明表
进一步的，基于目标城市河道黑臭等级与实际河道污染状况，判断各治理技术的重要性，选取合适的工艺措施，如表8，提出相应城市河道整治建议，或给出相应的黑臭、生境改善措施等初步方案，为目标河道整体治理与修复提供参考，实现“一河一策”与“模式治理”的有机结合。
表8不同级别黑臭河道治理技术的必要性选择
★数量越多代表实施该治理技术的必要性越强。
请参照图1和图2所示，本发明的实施例一为：
(1)确定评价对象与评价内容，建立三级指标评价体系。评价对象是某城区的N河道，评价内容是对N河道的黑臭程度进行评价。
(2)确定三级指标的评价标准，将三级指标按河道综合黑臭程度从低到高的顺序分为一、二、三、四、五、六级。
(2.1)确定三级指标的评价值，通过调研得到N河道的三级指标的数值，如表9所示。
表9城市河道三级指标数值表
(2.2)规范化处理三级指标评价向量，按照表1的三级指标顺序，对所有数值进行归一化处理，处理结果如下：(0.208，0.202，0.213，0.735，0.006，0.215，0.053，0.242，0.055，0.643，0.202，0.232，0.074，0.015，0.036，0.013，0.204)。
(3)运用物元层次分析法，得到客观科学的评价结果。
(3.1)建立评价指标矩阵：
(3.2)规范化处理评价指标矩阵，根据关联函数定义，确定关联度矩阵，如表10：
表10城市河道评价三级指标关联度矩阵表
(3.3)采用九标度法，确定各三级指标的权重，为W＝(0.079，0.074，0.0567，0.0537，0.066，0.047，0.078，0.060，0.057，0.060681876，0.069，0.057，0.053，0.041，0.0467，0.051，0.051)；
(3.4)单层次排序的一致性检验；
黑臭水平层次下，一致性比率CR＝CI/RI＝0.0297＜0.1，满足一致性要求；富营养化水平层次下，一致性比率CR＝CI/RI＝0.0189＜0.1，满足一致性要求；生物多样性水平层次下，一致性比率CR＝CI/RI＝0.0157＜0.1，满足一致性要求；水动力学水平层次下，一致性比率CR＝CI/RI＝0.0162＜0.1，满足一致性要求；
(3.5)层次总排序的一致性检验；
总排序计算表11如下所示：
表11城市河道评价总排序计算表
层次总排序的一致性比率满足一致性要求。
(3.6)德尔菲权重修正：
(3.6.1)采用德尔菲法来综合多个专家对指标权重的分配并进行修正。此方案征取10位专家的决策意见，并建立由专家评判结果组成的权重向量矩阵。
(3.6.2)专家判断矩阵的平均一致度：
计算10位专家两两之间评判权重向量Wi与Wj的夹角为Dij＝cos(θij)＝(0.991，0.952，0.929，0.945，0.989，0.967，0.942，0.958，0.937，0.982)，第i个专家判断矩阵的平均一致度Di＝(0.489，0.480，0.473，0.481，0.488，0.483，0.480，0.477，0.485，0.482)(i＝1，2，3···10)；
对所有专家的平均一致度进行归一化处理D'i＝(0.102，0.100，0.098，0.100，0.101，0.100，0.100，0.099，0.101，0.100)(i＝1，2，3···10)；
用归一化相对一致度D'i与相应权重向量相乘并求和，得修正后的权重：W'＝(0.049，0.043，0.051，0.046，0.047，0.047，0.103，0.045，0.048，0.062，0.051，0.048，0.061，0.084，0.070，0.090，0.055)；
(4)将各指标对应等级关联度值与修正后权重相乘，获得最终城市内河综合指标评价等级表，如下表12所示：
表12城市内河综合指标评价表
(5)待评河道相关物元关于综合评价等级的综合关联度为(-0.594，-0.584，-0.518，-0.385，-0.074，-0.067)，该结果表明该河道的黑臭程度在六级，河道水质极度差，生态环境自组织修复能力基本损失，黑臭现象及其严重，需采取物理、化学、生物等相结合的综合治理措施，加强河道管理，改善河道环境。
综上所述，本发明提供的一种河道水体综合评价方法，引入德尔菲法作为权重修正手段，形成了更为准确的综合评价方法，通过设置专家评判表、建立专家评价矩阵、计算评价向量夹角，并使用归一化后的专家评判一致度对权重修正，使得整体评价方法的准确度更高，修正的结果更加客观公正，适用的评价指标更为广泛的特性，同时对比分级的城市河道黑臭评价程度进行评级，对各等级对应治理措施，如控源截污、内源控污、水动力调控、曝气富氧、生物修复、生态修复等，并提供必要性建议，使得整个方法完整性更高，集收集、分析、结果和建议一体，形成流程化评价路线，系统性更强。
以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等同变换，或直接或间接运用在相关的技术领域，均同理包括在本发明的专利保护范围内。
Claims (10)
1.一种河道水体综合评价方法，其特征在于：包括以下步骤：步骤S1、建立三级指标评价体系，所述三级指标评价体系包括一级指标、二级指标和三级指标，所述三级指标包括多个等级；
步骤S2、通过物元分析法对河道水体的黑臭程度进行量化评价，包括建立三级指标的指标评价矩阵、对指标评价矩阵和指标评价矩阵所对应的指标向量进行归一化处理、以及确定归一化处理后的指标评价矩阵与所述三级指标的各个等级之间的关联度；
步骤S3、通过层次-德尔菲分析法确定关联度中各项三级指标的权重，接收多个针对所述各项三级指标的权重判断矩阵，计算各个权重判断矩阵的平均一致度并对所有平均一致度进行归一化处理得到相对一致度，使用所述相对一致度对三级指标的权重进行修正。
2.根据权利要求1所述的河道水体综合评价方法，其特征在于：所述一级指标为河道水体综合评价指标；所述二级指标包括黑臭水平、富营养化水平、生物多样性水平、水动力学水平；所述三级指标包括所述黑臭水平下的五日生化需氧量、化学需氧量、氨态氮、溶解氧、二价硫离子、色度、氧化还原电位和表观污染指数，所述富营养化水平下的硝态氮、总磷、总氮和叶绿素a，所述生物多样性水平下的Shannon-weaner多样性指数、Margalef多样性指数和Simpson多样性指数，水动力学水平下的断面水流平均速度和滞留时间。
3.根据权利要求1所述的河道水体综合评价方法，其特征在于：所述指标评价矩阵的具体建立方法为：将三级指标的指标评价矩阵设为Y'i，其中Y'i的三级指标对每一等级的评价值用L'ij(i＝1、2、……、17，j＝1、2、3、4、5、6)表示，i表示对应的第i个三级指标，j表示对应的一、二、三、四、五、六级。
4.根据权利要求3所述的河道水体综合评价方法，其特征在于：对指标评价矩阵和指标评价矩阵所对应的指标向量进行归一化处理的具体方法为：
对指标评价矩阵Y'i中的评价值进行归一化处理，使其评价值隶属于[0，1]的区间，若指标评价矩阵Y'i中三级指标的节域是严格的闭区间，将每个经典域的端点值除以节域的右端点值，若节域是[m，+∞)的半闭半开区间，则将m扩大四倍后的值设为节域右端点值，然后将每个经典域的端点值除以节域的右端点值，最后得到归一化后的指标评价矩阵Yi，其中每个三级指标的经典域用Lij表示；
对指标向量中的评价值的归一化处理，使其评价值隶属于[0，1]的区间，用b'分别除以Y'i中第i个三级指标归一化时所除以的数值，得到归一化后的指标向量B，其中B＝(b1、b2、…bi…b17)，bi(i＝1、2、……17)代表规范化处理后第i个三级指标的评价值，B代表规范化处理后的指标向量。
5.根据权利要求4所述的河道水体综合评价方法，其特征在于：所述经典域的物元矩阵表示方法为：
将三级指标的等级Q分为6个等级Q1、Q2、Q3、Q4、Q5、Q6，每个等级中各个指标En对应不同的量值范围，特征En的量值范围为Xpj＝[Apj，Bpj](p＝1，2，3，4，5，6)，则经典域的物元矩阵Kp可表示为：
所述节域的物元矩阵表示方法为：
Q0为各个等级的Qp构成，则将K0称为节域物元，其量值范围为X0j＝[A0j，B0j]，黑臭河道的综合黑臭程度的节域物元矩阵K0可表示为：
6.根据权利要求1所述的河道水体综合评价方法，其特征在于：所述关联度的确定方法为：
将关联度设为T，计算评价指标矩阵Gi关于等级j的关联度T：
上式中，W′为每个指标进过修正后的权重，T(x)为同一指标不同安全等级对应的关联函数值中最大的关联函数值，其中T(x)＝max(T(xj))；
根据归一化处理后的评价指标矩阵和指标评价向量，采用关联函数的定义，确定关联度矩阵T，关联矩阵由关联函数的每个值组成，为：
上式中：θ(X,X0)表示物元量值X到有限区间X0(X0＝[A，B])，A、B分别为有限元区间的最小值与最大值的距离；θ(X,XP)为待测物元量值X与有限区间Xp(Xp＝[Aj，Bj])的距离；X，X0，Xj分别为待评河道黑臭程度物元的量值、经典域物元的量值范围和节域物元的量值范围，点与有限元的间距计算公式为：
7.根据权利要求1所述的河道水体综合评价方法，其特征在于：所述三级指标的权重确定具体步骤为：
步骤S301、采用层次分析法构建判断矩阵，选取同一层级水平下的三级指标，两两比较评估三级指标的相对重要性，形成评估指标的比较矩阵F：
步骤S302、层次单排序与权重的确定，采用行和正规化方法，设n阶判断矩阵为：
按行求和计算有即得到权重向量W＝(v1、v2…、vn)T，经归一化使向量中各元素之和等于1后，得到向量W＝(w1、w2…、wn)T，W的元素即为同一层次因素对于上一层次因素某因素相对重要性的排序权值，次过程为层次单排序，其中，wi就是层次排序的优先程度，即各个指标的权重；
步骤S303、单层次排序的一致性检验，需要进行一致性检验确认层次单排序，一致性检验是指对A确定不一致的允许范围，定义CI为单层次排序的一致性指标，定义RI为单层次排序的随机一致性指标，定义一致性比率：CR＝CI/RI，当一致性比率CR＝CI/RI≤0.1时，认为A的不一致程度在容许范围之内，有满意的一致性，通过一致性检验，可用其归一化特征向量作为权向量，否则要重新构造成对比较矩阵A，对aij加以调整；
步骤S304、利用同一层次中所有层次单排序的结果，针对层次间所有因素，计算某一层次的层次总排序经归一化后，可以得到该层所有因素对于最高层相对重要性的权重值，假设上一层次所有因素A1、A2、…An，的总排序已经完成，的得到权重值分别为a1、a2、…an，与ai对应的本层次因素B1、B2、…Bn的单排序的结果为bi 1、bi 2、…bi n，如果Bj与Aj无关，则有
步骤S305、层次总排序的一致性检验，设B层B1、B2、…、Bn对上层(A层)中因素的层次总排序一致性指标为CI，层次总排序随机一致性指标为RI，则层次总排序的一致性比率为：
当CR≤0.1时，认为层次总排序计算结果具有良好的一致性。
8.根据权利要求1所述的河道水体综合评价方法，其特征在于：使用归一化后的专家评判一致度对权重修正的具体步骤为：
步骤S306、邀请多位专家参与决策，若由专家Z做出的权重判断矩阵形成的权重向量为：Wz＝(Wz1、Wz2、Wz3、…、Wzn)，(n＝1、2、…、17)，取权重向量Wi与Wj的夹角为Dij＝cos(θij)，Dij表示两专家各自权重判断矩阵Wi与Wj之间的一致度，依据Dij求第i个专家的权重判断矩阵的平均一致度Di：
步骤S307、对所有专家的平均一致度进行归一化处理，得到专家的权重判断矩阵的相对一致度D'i，用归一化相对一致度D'i与相应权重向量相乘并求和，得修正后的权重：
9.根据权利要求1所述的河道水体综合评价方法，其特征在于：步骤S3之后还包括：步骤S4、计算三级指标的关联度与对应三级指标修正后的权重值的乘积得到河道水体的黑臭程度。
10.根据权利要求9所述的河道水体综合评价方法，其特征在于：根据河道水体的黑臭程度对城市河道进行分级治理。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810833611.2A CN109190884A (zh) | 2018-07-26 | 2018-07-26 | 一种河道水体综合评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810833611.2A CN109190884A (zh) | 2018-07-26 | 2018-07-26 | 一种河道水体综合评价方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109190884A true CN109190884A (zh) | 2019-01-11 |
Family
ID=64937564
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810833611.2A Pending CN109190884A (zh) | 2018-07-26 | 2018-07-26 | 一种河道水体综合评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109190884A (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111721360A (zh) * | 2020-06-19 | 2020-09-29 | 东北农业大学 | 一种湿地生态环境监测系统及方法 |
| CN111950882A (zh) * | 2020-08-04 | 2020-11-17 | 四川省食品发酵工业研究设计院 | 一种基于综合感官分析的酒样分级方法 |
| CN112016838A (zh) * | 2020-08-31 | 2020-12-01 | 广东电网有限责任公司 | 配电网能效指标体系的贡献率计算方法、系统及终端设备 |
| CN112903945A (zh) * | 2021-02-26 | 2021-06-04 | 珠江水利委员会珠江水利科学研究院 | 一种基于常规水质参数的城市黑臭水体识别方法 |
| CN113554311A (zh) * | 2021-07-23 | 2021-10-26 | 中煤新集能源股份有限公司 | 地面定向孔注浆治理推覆体下奥灰水害工程质量评价方法 |
| CN114331061A (zh) * | 2021-12-15 | 2022-04-12 | 成都理工大学 | 一种城乡黑臭水体臭度评价方法 |
| CN114862285A (zh) * | 2022-07-06 | 2022-08-05 | 中国长江三峡集团有限公司 | 生态流量的评价方法和装置、电子设备和存储介质 |
| CN117314270A (zh) * | 2023-09-15 | 2023-12-29 | 国家海洋环境监测中心 | 一种海洋环境监测信息的综合评价方法 |
| CN111950882B (zh) * | 2020-08-04 | 2024-05-28 | 四川省食品发酵工业研究设计院有限公司 | 一种基于综合感官分析的酒样分级方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101894189A (zh) * | 2010-07-14 | 2010-11-24 | 中国矿业大学(北京) | 煤层底板突水评价的新方法 |
| CN104915575A (zh) * | 2015-07-07 | 2015-09-16 | 同济大学建筑设计研究院(集团)有限公司 | 一种基于层次物元可拓法的海绵城区生态指数的评价方法 |
| CN105975787A (zh) * | 2016-05-17 | 2016-09-28 | 福建海峡环保集团股份有限公司 | 一种基于层次分析法的城市内河生态多层级评测方法 |
-
2018
- 2018-07-26 CN CN201810833611.2A patent/CN109190884A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101894189A (zh) * | 2010-07-14 | 2010-11-24 | 中国矿业大学(北京) | 煤层底板突水评价的新方法 |
| CN104915575A (zh) * | 2015-07-07 | 2015-09-16 | 同济大学建筑设计研究院(集团)有限公司 | 一种基于层次物元可拓法的海绵城区生态指数的评价方法 |
| CN105975787A (zh) * | 2016-05-17 | 2016-09-28 | 福建海峡环保集团股份有限公司 | 一种基于层次分析法的城市内河生态多层级评测方法 |
Non-Patent Citations (5)
| Title |
|---|
| 刘传旺、吴建平、任胜伟、张旭: "基于层次分析法与物元分析法的水安全评价", 《水资源保护》 * |
| 周孝德等主编: "《水力学与水利信息学进展》", 31 October 2009, 西安交通大学出版社 * |
| 李雄峰、余宏明、李亚、师华鹏: "基于层次分析与可拓理论的水土流失严重程度分区评价", 《水土保持研究》 * |
| 裴海燕主编: "《南四湖藻类种群特征》", 31 December 2010, 中国环境科学出版社 * |
| 郭琨、艾万铸主编: "《海洋工作者手册 第1卷 海洋科技》", 31 August 2016, 海洋出版社 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111721360A (zh) * | 2020-06-19 | 2020-09-29 | 东北农业大学 | 一种湿地生态环境监测系统及方法 |
| CN111950882A (zh) * | 2020-08-04 | 2020-11-17 | 四川省食品发酵工业研究设计院 | 一种基于综合感官分析的酒样分级方法 |
| CN111950882B (zh) * | 2020-08-04 | 2024-05-28 | 四川省食品发酵工业研究设计院有限公司 | 一种基于综合感官分析的酒样分级方法 |
| CN112016838A (zh) * | 2020-08-31 | 2020-12-01 | 广东电网有限责任公司 | 配电网能效指标体系的贡献率计算方法、系统及终端设备 |
| CN112903945A (zh) * | 2021-02-26 | 2021-06-04 | 珠江水利委员会珠江水利科学研究院 | 一种基于常规水质参数的城市黑臭水体识别方法 |
| CN112903945B (zh) * | 2021-02-26 | 2022-02-22 | 珠江水利委员会珠江水利科学研究院 | 一种基于常规水质参数的城市黑臭水体识别方法 |
| CN113554311A (zh) * | 2021-07-23 | 2021-10-26 | 中煤新集能源股份有限公司 | 地面定向孔注浆治理推覆体下奥灰水害工程质量评价方法 |
| CN114331061A (zh) * | 2021-12-15 | 2022-04-12 | 成都理工大学 | 一种城乡黑臭水体臭度评价方法 |
| CN114862285A (zh) * | 2022-07-06 | 2022-08-05 | 中国长江三峡集团有限公司 | 生态流量的评价方法和装置、电子设备和存储介质 |
| CN117314270A (zh) * | 2023-09-15 | 2023-12-29 | 国家海洋环境监测中心 | 一种海洋环境监测信息的综合评价方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109190884A (zh) | 一种河道水体综合评价方法 | |
| Cheng et al. | The emissions reduction effect and technical progress effect of environmental regulation policy tools | |
| Zhao et al. | Impacts of environmental regulations on green economic growth in China: New guidelines regarding renewable energy and energy efficiency | |
| Feng et al. | A risk assessment model of water shortage based on information diffusion technology and its application in analyzing carrying capacity of water resources | |
| Song et al. | Realization of green transition based on the anti-driving mechanism: an analysis of environmental regulation from the perspective of resource dependence in China | |
| Preisner et al. | Mitigation of eutrophication caused by wastewater discharge: A simulation-based approach | |
| Wei et al. | Selection of optimum biological treatment for coking wastewater using analytic hierarchy process | |
| Zhou et al. | Spatio-temporal patterns and source apportionment of coastal water pollution in eastern Hong Kong | |
| Zhou et al. | Spatiotemporal variations and determinants of water pollutant discharge in the Yangtze River Economic Belt, China: A spatial econometric analysis | |
| Wang et al. | Spatial temporal patterns and driving factors of industrial pollution and structures in the Yangtze River Economic Belt | |
| Yuan et al. | Development of multidimensional water poverty in the Yangtze River Economic Belt, China | |
| CN105469196A (zh) | 一种矿井建设项目过程后评价的综合评价方法及系统 | |
| CN107895100A (zh) | 一种流域水质综合评价方法及系统 | |
| CN109190865B (zh) | 水质指标评价方法、预测方法、水域水质等级评价方法 | |
| Meng et al. | Comprehensive evaluation of ecological compensation effect in the Xiaohong River Basin, China | |
| CN112215389B (zh) | 一种确定河流环境流量过程区间的方法 | |
| CN103793604A (zh) | 一种基于相关向量机的污水处理软测量方法 | |
| CN104299067A (zh) | 一种可量化的企业发展情况评估模型 | |
| Wang et al. | Effects of industry structures on water quality in different urbanized regions using an improved entropy-weighted matter-elementmethodology | |
| CN107247888B (zh) | 基于储备池网络的污水处理出水总磷tp软测量方法 | |
| CN109871638A (zh) | 一种湖泊湿地富营养化评价模型构建方法 | |
| CN107664682A (zh) | 一种氨氮的水质软测量预测方法 | |
| Xu et al. | Assessment of water resources sustainability in mainland China in terms of water intensity and efficiency | |
| Su et al. | Water–energy–carbon nexus: greenhouse gas emissions from integrated urban drainage systems in China | |
| Cheng et al. | Trend prediction of carbon peak in China’s animal husbandry based on the empirical analysis of 31 provinces in China |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190111 |