基于综合评价的污染场地修复技术方案比选方法与流程

1.本发明涉及污染场地修复技术领域，尤其涉及一种基于综合评价的污染场地修复技术方案比选方法。

背景技术：

2.发达国家土壤和地下水修复工作起步较早，在大量案例和研究的支持下，无论是修复技术还是技术筛选方面都处于领先地位。美国环境保护署epa首先认识到传统的场地修复效果很大程度上依赖于决策者的主观判断，而决策者自身认识可能存在局限性和不稳定性，受认知偏见和复杂系统的影响，决策本身存在很大的不确定性。因此epa在1988年基于cercla(《环境应对补偿与责任综合法案》)，构建了一套修复技术方案的筛选流程，帮助决策者科学的选择修复技术，并将其以决策支持系统的形式实现。其中包括epa提出的专家评估法、美国frtr提出的类比筛选法，以及多目标决策分析法。并发展出多样的决策支持系统。这使决策过程得以更加开放、透明的呈现，决策过程可以重复检验，使决策过程更为稳定可靠。

3.我国污染场地修复技术筛选和修复技术方案比选工作的研究起步更晚，基于修复技术筛选形成的备选方案。方案比选的主要作用是选择经济效益、社会效益、环境效益综合表现最佳的技术方案，作为场地最终推荐的修复技术方案，为环境管理决策提供依据。目前选取的常用方法包括定量方法和定性方法。常用定量方法为ahp法、topsis+ahp法和专家评分法，但在实际工作中更偏重采用定性方法结合各指标对备选方案进行逐一比较和综合决策。目前采取的常用定量方法也存在一定问题，

4.ahp法的缺点是：

5.(1)定量数据少，定性成分多，不易令人信服，指标过多时数据统计量大，且权重不易确定，特征值和特征向量的精确求法比较复杂；

6.(2)所用的指标体系需要有专家系统的支持，如果给出的指标不合理则得到的结果就不准确。

7.topsis+ahp法的缺点是：

8.(1)求规范决策矩阵时比较复杂，不易求出正理想和负理想解，特别是当两个评价对象的指标值关于正理想解和负理想解的连线对称时，无法得出准确的结果；

9.(2)该方法以层次分析法为基础，指标过多时权重不易确定，定量数据少，定性成分多，评价具有主观性。

10.专家评分法的缺点是：

11.(1)完全依靠专家，主观性强，定性成分多；

12.(2)计算过程过于简单，方法不够科学。

13.因此，如何提出一种更好的决策方法，能够综合考虑方案比选中多项影响因素，解决现有方法中存在的问题，提高方案比选准确性，成为亟待解决的问题。

技术实现要素：

14.本发明的目的在于提供一种基于综合评价的污染场地修复技术方案比选方法。根据本发明的基于综合评价的污染场地修复技术方案比选方法，具有客观性、科学性和实用性，能够为国内污染场修复技术方案比选提供一定的理论和应用依据，本发明采用的技术方案如下：

15.根据本发明的一个方面，提供了一种基于综合评价的污染场地修复技术方案比选方法，包括以下步骤：

16.(一)建立多层次比选指标体系：制定并筛选出多个修复技术备选方案，通过构建多个一级评价指标，且每个所述一级评价指标对应有2-3个二级评价指标；

17.(二)确定各个评价指标权重：首先采用专家评估法对修复技术备选方案进行打分，获得多组评分信息，将多组评分信息取平均值，然后采用ahp层次分析法确定评价指标权重；

18.(三)获取标准化评价矩阵：利用灰色综合评价法对备选方案进行综合评判，并确定最优方案。

19.优选的，在所述步骤(一)中，所述一级指标设为4个，4个一级指标分别设为技术指标、经济指标、环境指标和社会指标，所述技术指标包括可操作性、污染物去除效率和修复时间，所述经济指标包括基本建设费用、运行费用和后期费用，所述环境指标包括残余风险、长期效果和健康影响，所述社会指标包括管理可接受程度和公众可接受程度。

20.优选的，在所述步骤(二)中确定评价指标权重的方法，包括以下步骤：

21.(1)构建判断矩阵，采用专家评估法对评价指标进行量化；

22.(2)计算特征向量及最大特征根：利用方根法求每个评价指标的特征向量，并进行归一化处理：

[0023][0024][0025]

计算判断矩阵的最大特征根：

[0026][0027]

(3)判断矩阵一致性检验：

[0028][0029][0030]

其中：为判断矩阵最大特征根；n为判断矩阵阶数；ri为与n对应的平均随机一致性取值；当《0.1时，认为判断矩阵具有良好的一致性，否则应调整判断矩阵元素的取值；

[0031]

(4)确定评价指标权重：根据判断矩阵一致性检验结果，确定评价指标权重。

[0032]

优选的，所述步骤(三)中灰色综合评价法，包括以下步骤：

[0033]

确定指标数据列v：设i为第i个修复技术备选方案的序号，i＝1，2，

…

，m；k为第k个评价指标的序号，k＝1，2，

…

，n；v

ik

为第i个修复技术备选方案的第k个指标的评价值，对于一个由m个修复技术备选方案，n个评价指标的系统，有下列矩阵：

[0034][0035]

确定参考数据列v0：取每个评价指标的最佳值的v

0k

参考数列v0的实体，于是有：

[0036]v0

＝(v

01

，v

02

，

…

，v

0n

)

[0037]

指标值规范化处理x：规范化的公式如下：

[0038][0039]

进行规范化处理之后，得：

[0040][0041]

x0＝(x

01

，x

02

，

…

，x

0n

)

[0042]

计算关联系数e：把规范化后的数列x0＝(x

01

，x

02

，

…

，x

0n

)作为参考数列，xi＝(x

i1

，x

i2

，

…

，x

in

)(i＝1，2，

…

，m)作为比较数列，关联系数的计算公式为：

[0043][0044]

式中，ρ为分辨系数，ρ∈[0，1]。

[0045]

得到下列关联系数矩阵：

[0046][0047]

式中为第i个评价单元第k个指标与最佳指标的关联系数。

[0048]

计算单层次关联度r：关联度计算方法采取权重乘以关联系数，根据层次分析法和专家评估法得到某一层的各指标相对于上层目标的优先权重为：

[0049]

w＝(w1，w2，

…

，wn)

[0050]

则关联度计算公式为：r＝(ri)1×m＝(r1，r2，

…

，rm)＝we

[0051]

计算多层评价系统的最终关联度：对一个由l层组成的多层评价系统，最终关联度

的计算方法如下：将第k层各指标的关联系数进行合成，分别得它们所属的上一层即k-1层各指标的关联度；然后把这一层所得到的关联度作为原始数据，继续合成得到第k-2层各指标的关联度，以此类推，直到求出最高层指标的关联度为止；

[0052]

综合评判：依据关联度大小进行排序，关联度的大小顺序即为备选方案的优劣次序。

[0053]

本发明采用的上述技术方案，具有如下显著效果：

[0054]

(1)本发明提供的方法首先建立多层次修复技术方案比选指标体系，其次利用ahp层次分析法和专家评估法对各评价指标进行量化，确定指标权重，最后采用灰色综合评价法对各备选方案进行综合评判，确定最优土壤/地下水修复技术方案，该评价指标体系符合修复工程实际情况，具有可操作性；

[0055]

(2)定性和定量相结合，使评价更加客观；

[0056]

(3)多种综合评价方法有效集成，使指标体系模型计算更加科学，具有客观性、科学性和实用性，能够为国内污染场的土壤和地下水修复技术方案比选提供一定的理论和应用依据。

附图说明

[0057]

图1是本发明的总流程图；

[0058]

图2是本发明的ahp法流程图；

[0059]

图3是本发明的灰色综合评价法流程图。

具体实施方式

[0060]

为使本发明的目的、技术方案及优点更加清楚明白，以下参照附图并举出优选实施例，对本发明进一步详细说明。然而，需要说明的是，说明书中列出的许多细节仅仅是为了使读者对本发明的一个或多个方面有一个透彻的理解，即便没有这些特定的细节也可以实现本发明的这些方面。

[0061]

根据本发明的一种基于综合评价的污染场地修复技术方案比选方法，首先，经过前期修复技术筛选和评估，该污染场地形成的修复技术备选方案有四套，分别为：

[0062]

方案1：异位固化/稳定化+阻隔填埋技术；

[0063]

方案2：水泥窑协同处置技术；

[0064]

方案3：砖窑协同处置技术；

[0065]

方案4：异位淋洗+生态恢复技术。

[0066]

其次，利用本发明提供的技术方法对四套方案进行评价，流程示意图如图1所示。

[0067]

(一)基于修复技术备选方案建立多层次比选指标体系：

[0068]

具体地，由于污染场地的土壤和地下水修复方案比选考虑的因素较多，本发明参考《建设用地土壤修复技术导则》(hj 25.5-2019)和工业企业场地环境调查评估与修复工作指南(试行)，建立评价指标体系。具体根据指标具有独立性、代表性与差异性、可行性三大原则，构建技术、经济、环境、社会4个一级指标，每个一级指标下层对应2-3个二级指标。即建立多个相互联系、相互作用的评价指标，按照一定层次结构组成有机整体。评价指标体系见表1。指标类型包括成本型和效益型，成本型指标越小越优，效益型指标越大越优。

[0069]

表1土壤/地下水修复技术方案评价指标体系

[0070]

[0071][0072]

(二)利用ahp层次分析法和专家评估法确定各个指标权重：

[0073]

具体地，首先采用专家评估法构建判断矩阵，获得备选方案评价指标评分信息，即专家踏勘和讨论后针对备选方案按照评价指标体系中各个指标进行打分，获得多组评分信息，并取平均值。然后采用ahp层次分析法计算特征根及对应的特征向量，并对判断矩阵进行一致性检验，检验通过后，确定各个指标权重，若检验未通过则返回到专家评估环节进行修正。

[0074]

主要步骤包括：

[0075]

(1)构建判断矩阵，采用专家评估法对评价指标进行量化，

[0076]

针对设定的评价指标体系，如表1，通过专家对5个同层次指标组{(u1，u2，u3，u4)，(u

11

，u

12

，u

13

)，(u

21

，u

22

，u

23

)，(u

31

，u

32

，u

33

)，(u

41

，u

42

)}中指标之间的相对重要性进行两两对比，以1～9的标度对对比结果进行量化，如表2。计算结果如表4-8。

[0077]

表2指标相对重要性等级表

[0078][0079]

表4 u

1-u4对u的重要程度及判断矩阵

[0080]

u评价指标体系u1技术指标u2经济指标u3环境指标u4社会指标u1技术指标1.00001.00002.00003.0000u2经济指标1.00001.00002.00003.0000u3环境指标0.50000.50001.00002.0000u4社会指标0.33330.33330.50001.0000

[0081]

表5 u

11-u

13

对u1的重要程度及判断矩阵

[0082]

u1技术指标u

11

可操作性u

12

污染物去除效率u

13

修复时间u

11

可操作性1.00000.33330.5000u

12

污染物去除效率3.00001.00002.0000u

13

修复时间2.00000.50001.0000

[0083]

表6 u

21-u

23

对u2的重要程度及判断矩阵

[0084]

u2经济指标u

21

基本建设费用u

22

运行费用u

23

后期费用u

21

基本建设费用1.00002.00003.0000u

22

运行费用0.50001.00002.0000u

23

后期费用0.33330.50001.0000

[0085]

表7 u31-u33对u3的重要程度及判断矩阵

[0086][0087][0088]

表8 u

41-u

42

对u4的重要程度及判断矩阵

[0089]

u4社会指标u

41

管理可接受程度u

42

公众可接受程度u

41

管理可接受程度1.00001.0000u

42

公众可接受程度1.00001.0000

[0090]

(2)计算特征向量及最大特征根：

[0091]

利用方根法求每个评价指标的特征向量，并进行归一化处理：

[0092][0093][0094]

计算判断矩阵的最大特征根λ

max

：

[0095][0096]

各判断矩阵特征向量及最大特征根由以上公式计算公式，计算结果见表9。

[0097]

表9各判断矩阵特征向量及最大特征根

[0098][0099]

(3)判断矩阵一致性检验：

[0100][0101][0102]

其中：λ

max

为判断矩阵最大特征根；n为判断矩阵阶数；ri为与n对应的平均随机一致性取值，见表3；当cr《0.1时，认为判断矩阵具有良好的一致性，否则应调整判断矩阵元素的取值。

[0103]

表3层次分析法中平均随机一致性指标ri的取值

[0104]

阶数ri102030.5840.9051.1261.2471.3281.4191.45

[0105]

各判断矩阵一致性检验计算公式见上述公式，计算结果见表9。由表9可知，各判断矩阵cr《0.1时，认为判断矩阵具有良好的一致性。

[0106]

(4)确定评价指标权重

[0107]

根据判断矩阵一致性检验结果，确定评价指标权重。

[0108]w1-4

＝[0.3512 0.3512 0.1887 0.1089]

[0109]w11-13

＝[0.5000 0.2500 0.2500]

[0110]w21-23

＝[0.5396 0.2970 0.1634]

[0111]w31-33

＝[0.2000 0.4000 0.4000]

[0112]w41-42

＝[0.5000 0.5000]

[0113]

(三)利用灰色综合评价法对备选方案进行综合评判：

[0114]

具体的，灰色综合评价法主要步骤包括：

[0115]

确定指标数据列v：设i为第i个修复技术备选方案的序号，i＝1，2，

…

，m；k为第k个评价指标的序号，k＝1，2，

…

，n；v

ik

为第i个修复技术备选方案的第k个指标的评价值，对于一个由m个修复技术备选方案，n个评价指标的系统，有下列矩阵：

[0116][0117]

则：v的值为：

[0118][0119]

确定参考数据列v0：取每个评价指标的最佳值的v

0k

参考数列v0的实体，于是有：

[0120]v0

＝(v

01

，v

02

，

…

，v

0n

)

[0121]

则：

[0122]v0

＝(4.0000，40000，40000，40000，40000，40000，40000，40000，40000，40000，40000，)

[0123]

指标值规范化处理x：规范化的公式如下：

[0124][0125]

进行规范化处理之后，得：

[0126][0127]

x0＝(x

01

，x

02

，

…

，x

0n

)

[0128]

则：

[0129][0130]

x0＝(1.0000，1.0000，1.0000，1.0000，1.0000，1.0000，1.0000，1.0000，1.0000，1.0000，1.0000)

[0131]

计算关联系数e：把规范化后的数列x0＝(x

01

，x

02

，

…

，x

0n

)作为参考数列，xi＝(x

i1

，x

i2

，

…

，x

in

)(i＝1，2，

…

，m)作为比较数列，关联系数的计算公式为：

[0132][0133]

式中，ρ为分辨系数，ρ∈[0，1]。

[0134]

得到下列关联系数矩阵：

[0135][0136]

式中为第i个评价单元第k个指标与最佳指标的关联系数。

[0137]

则：

[0138][0139]

(5)计算单层次关联度r：关联度计算方法采取指标权重乘以关联系数，根据层次

分析法和专家评估法得到某一层的各指标相对于上层目标的优先权重为：

[0140]

w＝(w1，w2，

…

，wn)

[0141]

则关联度计算公式为：r＝(ri)1×m＝(r1，r2，

…

，rm)＝we

[0142]

则：

[0143][0144][0145][0146][0147][0148]

计算多层评价系统的最终关联度：对一个由l层组成的多层评价系统，最终关联度的计算方法如下：将第k层各指标的关联系数进行合成，分别得它们所属的上一层即k-1层各指标的关联度；然后把这一层所得到的关联度作为原始数据，继续合成得到第k-2层各指标的关联度，以此类推，直到求出最高层指标的关联度为止。

[0149]

综合评判：依据关联度大小进行排序，关联度的大小顺序即为备选方案的优劣次序。由各层次关联度[0.3088 0.2176 0.2176 0.3887]得到备选方案排列顺序为：异位淋洗+生态恢复技术》异位固化/稳定化+阻隔填埋技术》水泥窑协同处置技术》砖窑协同处置技术，因此，异位淋洗+生态恢复技术为该污染场地的最优土壤重金属修复技术方案。

[0150]

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以作出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。技术特征：

1.一种基于综合评价的污染场地修复技术方案比选方法，其特征在于：包括以下步骤：(一)建立多层次比选指标体系：制定并筛选出多个修复技术备选方案，通过构建多个一级评价指标，且每个所述一级评价指标对应有2-3个二级评价指标；(二)确定各个评价指标权重：首先采用专家评估法对修复技术备选方案进行打分，获得多组评分信息，将多组评分信息取平均值，然后采用ahp层次分析法确定评价指标权重；(三)获取标准化评价矩阵：利用灰色综合评价法对备选方案进行综合评判，并确定最优方案。2.根据权利要求1所述的一种基于综合评价的污染场地修复技术方案比选方法，其特征在于：在所述步骤(一)中，所述一级指标设为4个，4个一级指标分别设为技术指标、经济指标、环境指标和社会指标，所述技术指标包括可操作性、污染物去除效率和修复时间，所述经济指标包括基本建设费用、运行费用和后期费用，所述环境指标包括残余风险、长期效果和健康影响，所述社会指标包括管理可接受程度和公众可接受程度。3.根据权利要求1所述的一种基于综合评价的污染场地修复技术方案比选方法，其特征在于：所述步骤(二)中确定评价指标权重的方法，包括以下步骤：(1)构建判断矩阵，采用专家评估法对评价指标进行量化；(2)计算特征向量及最大特征根：利用方根法求每个评价指标的特征向量，并进行归一化处理，公式如下：化处理，公式如下：计算判断矩阵的最大特征根：(3)判断矩阵一致性检验：(3)判断矩阵一致性检验：其中：λ

max

为判断矩阵最大特征根；n为判断矩阵阶数；ri为与n对应的平均随机一致性取值；(4)确定评价指标权重：根据判断矩阵一致性检验结果，确定评价指标权重。4.根据权利要求1所述的一种基于综合评价的污染场地修复技术方案比选方法，其特征在于：所述步骤(三)中灰色综合评价法，包括以下步骤：确定指标数据列v：设i为第i个修复技术备选方案的序号，i＝1，2，

…

，m；k为第k个评价指标的序号，k＝1，2，

…

，n；v

ik

为第i个修复技术备选方案的第k个评价指标的评价值，对于一个由m个修复技术备选方案，n个评价指标的系统，有下列矩阵：

确定参考数据列v0：取每个评价指标的最佳值的v

0k

参考数列v0的实体，于是有：v0＝(v

01

，v

02

，

…

，v

0n

)指标值规范化处理x：规范化的公式如下：进行规范化处理之后，得：x0＝(x

01

，x

02

，

…

，x

0n

)；计算关联系数e：把规范化后的数列x0＝(x

01

，x

02

，

…

，x

0n

)作为参考数列，x

i

＝(x

i1

，x

i2

，

…

，x

in

)(i＝1，2，

…

，m)作为比较数列，关联系数的计算公式为：式中，ρ为分辨系数，ρ∈[0，1]；得到下列关联系数矩阵：上式为第i个评价单元第k个指标与最佳指标的关联系数；计算单层次关联度r：关联度计算方法为指标权重乘以关联系数，根据层次分析法和专家评估法得到某一层的各指标相对于上一层目标的优先权重为：w＝(w1，w2，

…

，w

n

)则关联度计算公式为：r＝(r

i

)1×

m

＝(r1，r2，

…

，r

m

)＝we计算多层评价系统的最终关联度：对一个由l层组成的多层评价系统，最终关联度的计算方法如下：将第k层各指标的关联系数进行合成，分别得它们所属的上一层即k-1层各指标的关联度；然后把这一层所得到的关联度作为原始数据，继续合成得到第k-2层各指标的关联度，以此类推，直到求出最高层指标的关联度为止；综合评判：依据关联度大小进行排序，关联度的大小顺序即为备选方案的优劣次序。

技术总结

本发明公开了一种基于综合评价的污染场地修复技术方案比选方法。该方法包括步骤：(1)根据修复技术备选方案，建立修复技术方案比选评价指标体系；(2)采用AHP层次分析法和专家评价法确定评价指标权重；(3)采用灰色综合评价法进行备选方案排序，并确定最优修复技术方案。本发明的评价指标体系符合修复工程实际情况，具有可操作性；定性和定量相结合，使评价更加客观；多种综合评价方法有效集成，使指标体系模型计算更加科学，具有客观性、科学性和实用性，能够为国内污染场修复技术方案比选提供理论和应用依据。理论和应用依据。理论和应用依据。

技术研发人员：宋海农 郑浩 李雄光 刘彬 韩丽杰 刘亚东 赵粉红 郑子健 赵振

受保护的技术使用者：广西博世科环保科技股份有限公司

技术研发日：2022.03.28

技术公布日：2022/6/21