稳定化材料选择方法、装置、计算机设备、介质和产品

1.本技术涉及土壤修复技术领域，特别是涉及一种稳定化材料选择方法、装置、计算机设备、介质和产品。


背景技术：

2.随着工农业的急速发展，我国土壤污染日益形势严峻，其中，以重金属为代表的无机型污染问题尤为突出。土壤受到重金属污染后，不仅会导致粮食减产，更有甚者还可能对人类的身体健康造成损伤，因此亟需探究一种有效的土壤修复技术。
3.针对土壤重金属污染，通常使用稳定化材料对土壤进行修复。然而，随着技术的发展，稳定化材料的种类日益增加，其经济成本不同，修复重金属污染的效果也各有差异。传统方法中，选择稳定化材料时，通常是相关负责人员凭工程经验选择熟知的稳定化材料。
4.然而，通过人工选择稳定化材料时，往往由于考虑因素不够全面，因此，所选择的稳定化材料不够准确。进而，导致采用所选择的稳定化材料对土壤进行修复时，难以达到较好的修复效果。


技术实现要素：

5.基于此，有必要针对上述技术问题，提供一种能够提升稳定化材料选择准确度的稳定化材料选择方法、装置、计算机设备、介质和产品。
6.第一方面，本技术提供了一种稳定化材料选择方法。所述方法包括：
7.获取待修复土壤的基本信息及修复长效性指标的初始值；基本信息至少包括超标污染物类型及土壤回用用途；
8.根据超标污染物类型确定至少一种初始稳定化材料；
9.针对各初始稳定化材料，根据初始稳定化材料的经济指标、技术指标、环境指标及土壤回用用途，确定初始稳定化材料的综合系数；技术指标包括对修复长效性指标的初始值进行更新得到的修复长效性指标的目标值；
10.根据初始稳定化材料的综合系数，从各初始稳定化材料中确定目标稳定化材料。
11.在其中一个实施例中，经济指标包括材料总成本，技术指标包括修复效果参数、修复周期以及修复长效性指标的目标值，针对各初始稳定化材料，根据初始稳定化材料的经济指标、技术指标、环境指标及土壤回用用途，确定初始稳定化材料的综合系数，包括：
12.针对各初始稳定化材料，获取初始稳定化材料的材料总成本；材料总成本包括材料成本及运行成本；
13.针对各初始稳定化材料，获取初始稳定化材料的修复效果参数、修复周期以及修复长效性指标的目标值；
14.针对各初始稳定化材料，获取初始稳定化材料的环境指标；环境指标包括二氧化碳产量、颗粒物产量、固体废弃物产量、废水产量中的至少一种；
15.根据各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标
的目标值、环境指标及土壤回用用途，确定初始稳定化材料的综合系数。
16.在其中一个实施例中，针对各初始稳定化材料，获取初始稳定化材料的修复效果参数、修复周期以及修复长效性指标的目标值，包括：
17.在所述修复长效性指标处于初始值的情况下，确定所述待修复土壤的第一超标污染物浓度及向所述待修复土壤分别添加各所述初始稳定化材料后的第二超标污染物浓度，将所述第二超标污染物浓度与所述第一超标污染物浓度的第一比值作为各所述初始稳定化材料的修复效果参数；
18.根据达到各初始稳定化材料的修复效果参数对应的天数，确定各初始稳定化材料的修复周期；
19.获取向所述待修复土壤分别添加各所述初始稳定化材料后的第三超标污染物浓度，根据修复长效性指标的初始值、所述第二超标污染物浓度及所述第三超标污染物浓度，确定各所述初始稳定化材料的修复长效性指标的目标值；所述第三超标污染物浓度与所述修复长效性指标的中间值相对应。
20.在其中一个实施例中，根据修复长效性指标的初始值、第二超标污染物浓度及第三超标污染物浓度，确定各初始稳定化材料的修复长效性指标的目标值，包括：
21.根据修复长效性指标的初始值，计算修复长效性指标的中间值；
22.在修复长效性指标处于中间值的情况下，确定向待修复土壤分别添加各初始稳定化材料后的第三超标污染物浓度；
23.根据第二超标污染物浓度及第三超标污染物浓度，确定各初始稳定化材料的修复长效性指标的目标值。
24.在其中一个实施例中，根据第二超标污染物浓度及第三超标污染物浓度，确定各初始稳定化材料的修复长效性指标的目标值，包括：
25.针对各初始稳定化材料，确定第二超标污染物浓度及各第三超标污染物浓度的第二比值；
26.针对各初始稳定化材料，基于各第二比值建立原始矩阵；
27.对原始矩阵进行归一化处理，生成第一目标矩阵，并基于原始矩阵及第一目标矩阵确定各修复长效性指标的第一权重向量；
28.根据第一权重向量及第一目标矩阵确定各初始稳定化材料的修复长效性指标的目标值。
29.在其中一个实施例中，根据各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标及土壤回用用途，确定初始稳定化材料的综合系数，包括：
30.根据土壤回用用途，建立关于各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标的判断矩阵；判断矩阵用于表征材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标中相邻两个指标的重要程度；
31.计算判断矩阵对应的初始权重向量及最大特征根；
32.根据最大特征根对判断矩阵进行修正，直至判断矩阵的一致性满足预设条件，根据当前判断矩阵确定第二权重向量；
33.根据第二权重向量以及各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标，确定初始稳定化材料的综合系数。
34.在其中一个实施例中，根据第二权重向量以及各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标，确定初始稳定化材料的综合系数，包括：
35.对各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标进行归一化处理，得到各指标的归一化处理结果；
36.基于归一化处理结果建立第二目标矩阵；
37.根据第二目标矩阵及第二权重向量，计算各初始稳定化材料的综合系数。
38.在其中一个实施例中，根据初始稳定化材料的综合系数，从各初始稳定化材料中确定目标稳定化材料，包括：
39.根据初始稳定化材料的综合系数对各初始稳定化材料进行排序，得到排序结果；
40.从排序结果中确定候选稳定化材料；候选稳定化材料的综合系数在排序结果中位于前预设名次；
41.基于候选稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标，从各候选稳定化材料中确定目标稳定化材料。
42.在其中一个实施例中，修复长效性指标至少包括土壤ph、土壤有机质含量、土壤粘粒含量、当地年均最低温度，年均降雨量，冻融循环次数、二氧化碳浓度、紫外辐射强度。
43.第二方面，本技术还提供了一种稳定化材料选择装置。所述装置包括：
44.获取模块，用于获取待修复土壤的基本信息及修复长效性指标的初始值；基本信息至少包括超标污染物类型及土壤回用用途；
45.初始材料确定模块，根据超标污染物类型确定至少一种初始稳定化材料；
46.系数确定模块，针对各初始稳定化材料，根据初始稳定化材料的经济指标、技术指标、环境指标及土壤回用用途，确定初始稳定化材料的综合系数；技术指标包括对修复长效性指标的初始值进行更新得到的修复长效性指标的目标值；
47.目标材料确定模块，根据初始稳定化材料的综合系数，从各初始稳定化材料中确定目标稳定化材料。
48.第三方面，本技术还提供了一种计算机设备。所述计算机设备包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现任意上述第一方面中的步骤。
49.第四方面，本技术还提供了一种计算机可读存储介质。所述计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现任意上述第一方面中的步骤。
50.第五方面，本技术还提供了一种计算机程序产品。所述计算机程序产品，包括计算机程序，该计算机程序被处理器执行时实现上述任意第一方面中的步骤。
51.上述稳定化材料选择方法、装置、计算机设备、介质和产品，通过获取待修复土壤的基本信息及修复长效性指标的初始值；基本信息至少包括超标污染物类型及土壤回用用途；根据超标污染物类型确定至少一种初始稳定化材料；针对各初始稳定化材料，根据初始稳定化材料的经济指标、技术指标、环境指标及土壤回用用途，确定初始稳定化材料的综合
系数；技术指标包括对修复长效性指标的初始值进行更新得到的修复长效性指标的目标值；根据初始稳定化材料的综合系数，从各初始稳定化材料中确定目标稳定化材料。上述稳定化材料选择方法，首先根据土壤的超标污染物类型先选出适用的初始稳定化材料，再结合稳定化材料的经济指标、技术指标、环境指标及回用用途对初始稳定化材料的综合系数进行计算，而其中技术指标包含了对修复长效性指标的初始值进行更新得到的修复长效性指标的目标值。在此过程中，由于计算综合系数的过程中充分考虑了选取的初始稳定化材料的经济指标、技术指标、环境指标，同时其中的修复长效性指标的目标值是具有工程意义的指标，其考虑到了稳定化材料之间可能存在的差异，且计算综合系数时，还考虑了土壤回用用途，因此上述稳定化材料选择方法能够比较综合且全面的对初始稳定化材料进行评价，进而计算出各初始稳定化材料的综合系数。进一步的，基于此综合系数也能够更准确且便捷的从各初始稳定化材料中确定对土壤修复效果更佳的目标稳定化材料。
附图说明
52.图1为一个实施例中稳定化材料选择方法的应用环境图；
53.图2为一个实施例中稳定化材料选择方法的流程示意图；
54.图3为图2中根据初始稳定化材料的经济指标、技术指标、环境指标及土壤回用用途，确定初始稳定化材料的综合系数的方法的流程示意图；
55.图4为一个实施例中稳定化材料选择方法的评价体系的示意图；
56.图5为图3中获取初始稳定化材料的修复效果参数、修复周期以及修复长效性指标的目标值的方法的流程示意图；
57.图6为一个实施例中修复长效性指标的示意图；
58.图7为图5中根据所述修复长效性指标的初始值、第二超标污染物浓度及第三超标污染物浓度，确定各初始稳定化材料的修复长效性指标的目标值的方法的示意图；
59.图8为图7中根据所述修复长效性指标的初始值、第二超标污染物浓度及第三超标污染物浓度，确定各初始稳定化材料的修复长效性指标的目标值的的方法的示意图；
60.图9为图3中根据各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标及土壤回用用途，确定初始稳定化材料的综合系数的方法的示意图；
61.图10为图9中根据各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标及土壤回用用途，确定初始稳定化材料的综合系数的方法的示意图；
62.图11为图9中根据各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标及土壤回用用途，确定初始稳定化材料的综合系数的方法的示意图；
63.图12为一个实施例中稳定化材料选择方法的简要流程图；
64.图13为一个实施例中稳定化材料选择装置的结构框图；
65.图14为一个实施例中计算机设备的内部结构图。
具体实施方式
66.为了使本技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本技术进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本技术，并不用于限定本技术。
67.随着工农业的急速发展，我国土壤污染日益形势严峻，其中，以重金属为代表的无机型污染问题尤为突出。2014年环境保护部与国土资源部联合发布的《全国土壤污染状况调查公报》显示，全国土壤污染点位超标率16.1％，超标点位中无机型污染占比高达82.8％。土壤受到重金属污染后，不仅会导致粮食减产，更有甚者还可能对人类的身体健康造成损伤，因此亟需探究一种有效的土壤修复技术。
68.针对土壤重金属污染，通常使用稳定化材料对土壤进行修复。然而，随着技术的发展，稳定化材料的种类日益增加，其经济成本不同，修复重金属污染的效果也各有差异。传统方法中，选择稳定化材料时，通常是相关负责人员凭工程经验选择熟知的稳定化材料。
69.然而，人工选择稳定化材料时，往往由于考虑因素不够全面，因此，所选择的稳定化材料不够准确。进而，导致采用所选择的稳定化材料对土壤进行修复时，难以达到较好的修复效果。
70.本技术实施例提供的稳定化材料选择方法，可以应用于如图1所示的应用环境中。其中，终端102通过网络与服务器104进行通信。数据存储系统可以存储服务器104需要处理的数据。数据存储系统可以集成在服务器104上，也可以放在云上或其他网络服务器上。服务器104能够获取待修复土壤的基本信息及修复长效性指标的初始值；基本信息至少包括超标污染物类型及土壤回用用途；根据超标污染物类型确定至少一种初始稳定化材料；针对各初始稳定化材料，根据初始稳定化材料的经济指标、技术指标、环境指标及土壤回用用途，确定初始稳定化材料的综合系数；技术指标包括对修复长效性指标的初始值进行更新得到的修复长效性指标的目标值；根据初始稳定化材料的综合系数，从各初始稳定化材料中确定目标稳定化材料。其中，终端102可以但不限于是各种个人计算机、笔记本电脑、智能手机、平板电脑、物联网设备和便携式可穿戴设备，物联网设备可为智能音箱、智能电视、智能空调、智能车载设备等。便携式可穿戴设备可为智能手表、智能手环、头戴设备等。服务器104可以用独立的服务器或者是多个服务器组成的服务器集群来实现。
71.在一个实施例中，如图2所示，提供了一种稳定化材料选择方法，以该方法应用于图1中的服务器104为例进行说明，包括以下步骤：
72.s220，获取待修复土壤的基本信息及修复长效性指标的初始值。
73.待修复土壤的基本信息至少包括超标污染物类型及土壤回用用途，超标污染物类型主要是重金属元素污染，如cr、ni、cu、zn、as、cd、sb、hg、tl、pb等重金属污染；土壤回用用途主要包括卫生填埋或安全填埋、回填用土、绿化用土、筑路材料。除此之外，待修复土壤的基本信息还可以包括修复方式以及待修复土壤所在地块的基本信息等，例如地块的项目名称、项目地点、用地规划。其中，修复方式主要包括原位修复、原地异位修复及异地异位修复，用地规划种类包括《土壤环境质量建设用地土壤污染风险管控标准》(gb36600-2018)中的各种类，例如居住用地、中小学用地用地、医疗卫生用地、社会福利设施用地、社区公园或儿童公园用地、工业用地、物流仓储用地、商业服务业设施用地、道路与交通设施用地、公用设施用地、公共管理与公共服务用地、绿地与广场用地。
74.修复长效性指标包括对待修复土壤进行修复时影响修复长期有效性的各指标，可以包括土壤所在地基本指标，例如年均最低温度等，还可以包括土壤性质方面的指标，如土壤ph等。修复长效性中的基本指标可以采取现场踏勘、人员访谈及资料收集等方式获取，具体获取方式可参考《hj25.1-2019建设用地土壤污染状况调查技术导则》，土壤性质方面的指标可以通过对待修复土壤进行实地调查及采样进而得到。在获取长效性修复指标的初始值及超标污染物类型之后，可以结合用地规划及实际需求，对修复方式及修复后的土壤回用用途进行选择。
75.s240，根据超标污染物类型确定至少一种初始稳定化材料。
76.初始稳定化材料为能够有效应对超标污染物的稳定化材料，在待修复土壤中加入稳定化材料后，稳定化材料能够固定待修复土壤中的超标污染物，降低超标污染物迁移的风险，从而达到修复土壤的目的。其中，超标污染物类型包括重金属cr、ni、cu、zn、as、cd、sb、hg、tl、pb。在获取超标污染物类型后，可以根据污染物类型，从表1中选取有效的稳定化材料作为初始稳定化材料。
77.表1各重金属污染对应的可选用的初始稳定化材料
78.[0079][0080]
表1中“√”标记表明该稳定化材料可有效处理该列对应的重金属，例如若待修复土壤的超标污染物类型为一种重金属，如cr，那cr对应的一列中有“√”标记的稳定化材料都可作为初始稳定化材料。若待修复土壤的超标污染物类型为多种重金属，那么这多种重金属所在列都有“√”标记的稳定化材料都可作为初始稳定化材料。
[0081]
s260，针对各初始稳定化材料，根据初始稳定化材料的经济指标、技术指标、环境指标及土壤回用用途，确定初始稳定化材料的综合系数。
[0082]
初始稳定化材料的经济指标反应使用该稳定化材料的成本，技术指标是包括对修复长效性指标的初始值进行更新得到的修复长效性指标的目标值，而修复长效性指标的初始值反应了应用该稳定化材料对待修复土壤进行修复时的修复长期有效性，环境指标反应了采用该初始稳定化材料进行修复时对环境造成的影响，而土壤回用用途是对待修复土壤进行修复后土壤的用途。
[0083]
不同土壤回用用途下，初始稳定化材料的经济指标、技术指标、环境指标的重要程度也不同。上述指标都将直接影响稳定化材料的选择，因此需要根据初始稳定化材料的经济指标、技术指标、环境指标及土壤回用用途，共同确定各初始稳定化材料的综合系数。具体的确定方法可以针对不同的土壤回用方式，给初始稳定化材料的经济指标、技术指标、环境指标分配不同的权重，最终综合各指标得到各初始稳定化材料的综合系数，当然也可以采用其他方式，本技术对此不做限定。
[0084]
s280，根据初始稳定化材料的综合系数，从各初始稳定化材料中确定目标稳定化材料。
[0085]
在确定多种初始稳定化材料的综合系数后，可以根据综合系数的具体值初步筛选得到几种稳定化材料，再结合选出的稳定化材料的经济指标、技术指标、环境指标综合评价，确定目标稳定化材料。
[0086]
上述稳定化材料选择方法，通过待修复土壤的超标污染物类型可以先选出适用的初始稳定化材料，再结合稳定化材料的经济指标、技术指标、环境指标及回用用途对初始稳定化材料的综合系数进行计算，而其中技术指标包含了对修复长效性指标的初始值进行更新得到的修复长效性指标的目标值。在此过程中，由于计算综合系数的过程中充分考虑了选取的初始稳定化材料的经济指标、技术指标、环境指标，同时其中的修复长效性指标的目标值考虑到了稳定化材料之间的可能存在的差异，且计算综合系数时，还考虑了土壤回用用途，因此上述稳定化材料选择方法能够比较综合且全面的对初始稳定化材料进行评价，进而准确的计算出各初始稳定化材料的综合系数。进一步的，基于此综合系数也能够更准确的从各初始稳定化材料中确定对土壤修复效果更佳目标稳定化材料。
[0087]
在上面的实施例中，提到了确定初始稳定化材料的综合系数，下面就对其具体方法进行介绍。在一个实施例中，经济指标包括材料总成本，技术指标包括修复效果参数、修复周期以及修复长效性指标的目标值。如图3所示，s260包括：
[0088]
s262，针对各初始稳定化材料，获取初始稳定化材料的材料总成本。
[0089]
材料总成本包括材料成本及运行成本。材料总成本可以通过计算或直接收集该初始稳定化材料在其他工程中应用时的数据得到。材料总成本b
i,1
的计算公式如下：
[0090]bi,1
＝cmi+coiꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式(1-1)
[0091]
其中，cmi为第i种初始稳定化材料的材料成本，coi为第i种初始稳定化材料的运行成本。
[0092]
cmi＝m
i piꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式(1-2)
[0093]
其中，mi为采用第i种初始稳定化材料对待修复土壤进行修复时所需的初始稳定化材料的质量，pi表示第i种初始稳定化材料的单价。
[0094]
coi＝mi/ρi·oꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式(1-3)
[0095]
其中，o为不同修复方式下，单位体积材料的运行费用；ρi为材料i的密度。
[0096]
s264，针对各初始稳定化材料，获取初始稳定化材料的修复效果参数、修复周期以及修复长效性指标的目标值。
[0097]
初始稳定化材料的修复效果参数可以通过浸出实验得到，具体可以参考《重金属污染场地土壤修复标准》(db43/t1165-2016)，使用《hj557-2010-固体废物浸出毒性浸出方法水平振荡法》中的方法对待修复土壤进行浸出实验。简单来说浸出实验就是取固定量的待修复土壤，在待修复土壤中加入一种初始稳定化材料以及浸提剂，此处浸提剂可以选用水，随后将修复土壤振荡一定时间后，测定浸出液中重金属浸出浓度，最终通过多次浸出实验得到各初始稳定化材料的修复效果参数。
[0098]
修复周期可以根据修复效果参数b
i,2
对应的天数确定，例如可以选择添加初始稳定化材料后修复效果参数最优时对应的天数作为修复周期，当然也可以采用其他方式确定。
[0099]
修复长效性指标的目标值是将多个长效性指标进行处理得到的一个目标值，例如当修复长效性指标为8个时，综合考虑这8个修复长效性指标的重要程度，最终得到一个修复长效性指标的目标值。
[0100]
s266，针对各初始稳定化材料，获取初始稳定化材料的环境指标。
[0101]
其中，环境指标包括二氧化碳产量、颗粒物产量、固体废弃物产量、废水产量中的至少一种，每一种稳定化材料对应的环境指标均不同，环境指标可以采用lca(生命周期评价，life cycle assessment)法进行评价得到，也可以直接从lca数据库中收集各初始稳定化材料的环境指标，当然也可以采用其他方法，本技术对此不做限定。
[0102]
s268，根据各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标及土壤回用用途，确定初始稳定化材料的综合系数。
[0103]
针对不同回用用途，由于材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标等的重要性不同，因此需要先对材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标等重要性进行评估，再基于以上信息确定初始稳定化材料的综合系数。
[0104]
以图4对本技术中的稳定化材料选择方法的评价体系进行简单总结。本技术中的初始稳定化材料的综合系数可以由经济指标a1、技术指标a2及环境指标a3确定，其中经济指标a1可以包括材料总成本b1，技术指标a2包括修复效果参数b2、修复周期b3以及目标修
复长效性b4，环境指标a3可以包括二氧化碳产量b5、颗粒物产量b6、固体废弃物产量b7、废水产量b8。
[0105]
在本实施例中，获取的指标包括初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值以及环境指标，由于综合考虑了以上各指标以及土壤回用用途，因此上述稳定化材料选择方法能够全方面的对稳定化材料进行评价，进而准确的得出初始稳定化材料的综合系数。
[0106]
以上实施例介绍了获取初始稳定化材料的修复效果参数、修复周期以及修复长效性指标的目标值，下面就对其具体方法进行介绍。在一个实施例中，如图5所示，s264包括：
[0107]
s502，在修复长效性指标处于初始值的情况下，确定待修复土壤的第一超标污染物浓度及向待修复土壤分别添加各初始稳定化材料后的第二超标污染物浓度，将第二超标污染物浓度与第一超标污染物浓度的第一比值作为各初始稳定化材料的修复效果参数。
[0108]
图6为修复长效性指标的初始值包含的各指标的示意图，可选的，修复长效性指标至少包括土壤ph、土壤有机质含量、土壤粘粒含量、年均最低温度，年均降雨量，冻融循环次数、二氧化碳浓度、紫外辐射强度。其中，土壤ph测定方法可参照hj962-2018，ny/t1377-2007中的方法，土壤有机质含量测定方法可参照hj658-2013、hj615-2011及gb9834-88中的方法，土壤粘粒含量测定方法可参照hj1068-2019，gb7872-1987中的方法，待修复土壤中的重金属测定方法可参照gb36600-2018中的方法。待修复土壤的土壤ph、土壤有机质含量、土壤粘粒含量为按照上述方法在多个点位测量之后求平均值得到的结果。
[0109]
当待修复土壤中超标污染物类型为一种，即只包含一种重金属时，针对待修复土壤，已经获取了该待修复土壤的修复长效性指标的初始值，那么在修复长效性指标处于初始值的情况下，先在不添加稳定化材料时进行浸出实验得到第一超标污染物，随后向所述待修复土壤分别添加各所述初始稳定化材料后再进行浸出实验，得到重金属的浸出浓度，即第二超标污染物浓度。最后将第二超标污染物浓度与第一超标污染物浓度的比值作为各初始稳定化材料的修复效果参数，修复效果参数b
i,2
的计算公式如下：
[0110][0111]
其中，c
im,i
为添加第i种初始稳定化材料后得到的第二超标污染物浓度，c
or
为未添加初始稳定化材料时的第一超标污染物浓度。
[0112]
当待修复土壤中超标污染物类型为多种，即重金属种类为多种时，可以针对于不同的重金属进行浸出实验，将待修复土壤中分别添加各稳定化材料后浸出液浓度降低最多时的比率作为本技术的修复效果参数b
i,2
，此时计算公式如下：
[0113]bi,2
＝min{b
i,2,1
,b
i,2,2
,...,b
i,2,k
}
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式(1-5)
[0114]
其中，k表示第k种重金属。
[0115]
s504，根据达到各初始稳定化材料的修复效果参数对应的天数，确定各初始稳定化材料的修复周期。
[0116]
当待修复土壤中超标污染物类型为一种，即只包含一种重金属，将添加初始稳定化材料后达到最优修复效果参数时对应的天数作为修复周期。
[0117]
当待修复土壤中超标污染物类型为多种时，针对不同的重金属，将最长修复周期
作为该初始稳定化材料的修复周期b
i,3
，此时计算公式如下：
[0118]bi,3
＝max{b
i,3,1
,b
i,3,2
,...b
i,3,k
}
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式(1-6)
[0119]
其中，k表示第k种重金属。
[0120]
s506，获取向待修复土壤分别添加各初始稳定化材料后的第三超标污染物浓度，根据所述修复长效性指标的初始值、第二超标污染物浓度及第三超标污染物浓度，确定各初始稳定化材料的修复长效性指标的目标值。
[0121]
其中，第三超标污染物浓度与修复长效性指标的中间值相对应，各修复长效性指标的中间值可以是极端情况下的值。
[0122]
首先先获取修复长效性指标的中间值，再在修复长性指标处于中间值时，向所述待修复土壤分别添加各所述初始稳定化材料，并通过浸出实验得到第三超标污染物浓度。随后基于每种初始化稳定材料的多个第三超标污染物浓度以及第二超标污染物浓度进行归一化处理，得到一个修复长效性指标的目标值。
[0123]
在本实施例中，修复长效性指标包含了土壤ph、土壤有机质含量、土壤粘粒含量等多个指标，通过获取待修复土壤未添加初始稳定化材料时的第一超标污染物浓度及添加各所述初始稳定化材料后第二超标污染物浓度来确定修复效果参数及修复周期，由于修复效果参数、修复周期及修复长效性指标的目标值是根据多个修复长效性指标确定的，因此得到的修复效果参数、修复周期及修复长效性指标的目标值能全面的表征初始稳定化材料的性能。
[0124]
以上实施例提到了确定各初始稳定化材料的修复长效性指标的目标值，下面就对其具体方法进行介绍。在一个实施例中，如图7所示，根据修复长效性指标的初始值、第二超标污染物浓度及第三超标污染物浓度，确定各初始稳定化材料的修复长效性指标的目标值，包括：
[0125]
s702，根据修复长效性指标的初始值，计算修复长效性指标的中间值。
[0126]
其中，修复长效性指标的中间值考察的是初始稳定化材料在极端条件下的修复效果，其反映了各修复长效性指标在极端条件下的超标污染物迁移风险。因此需根据修复长效性指标的初始值确定，例如初始稳定化材料的中间值可以是初始值乘以或加上安全系数10得到的，也可以理解为各修复长效性指标的中间值为其在极端情况下的值。
[0127]
具体的，对于土壤ph，在确定第三超标污染物浓度时，需要使用《污染土壤修复工程技术规范固化稳定化(征求意见稿)》中的浸出方法hj299，模拟酸雨淋滤对稳定化材料的影响，因此以酸雨淋浸时的情况作为ph的中间值。对于土壤有机质含量和土壤粘粒含量，这两个值的取值范围仅与场地有关，因此将土壤有机质含量最高及土壤粘粒含量时最低作为极端情况，即将待修复土壤的土壤有机质含量最高时对应的值作为土壤有机质含量的中间值，将土壤粘粒含量最低时对应的值作为土壤粘粒含量的中间值。将将年均最低温度的中间值设定为当地年均最低温度减10，将年均降雨量、冻融循环次数、二氧化碳浓度、紫外辐射强度的中间值设定为各初始值对应乘以10，例如降雨量的中间值等于年均降雨量乘以10。
[0128]
s704，在所述修复长效性指标处于中间值的情况下，确定向所述待修复土壤分别添加各所述初始稳定化材料后的第三超标污染物浓度。
[0129]
具体的，采用hj299中所涉及的方法，模拟酸雨淋浸情况在待修复土壤中分别加入
importance though intercrieria correlation))计算得到，critic法是基于各修复长效性指标的对比强度和指标之间的冲突性来综合衡量指标的客观权重。critic法考虑了指标变异性大小的同时兼顾指标之间的相关性，与修复长效性指标的特征吻合。以ph举例说明，若所有材料受ph变化的影响相似，说明相较差异性更大的其他修复长效性指标，ph对修复长效性的影响不显著，应该赋予较小的权重。同样以ph举例说明，如果ph的影响与有机质含量的影响间存在较强的正相关，说明其冲突性越小，表明这两个指标在评价方案的优劣上反映的信息有较大的相似性，甚至仅用ph即可代表有机质的影响，因此这两者的权重应降低。
[0145]
具体的，由于上述原始矩阵中的第二比值为基于不同修复长效性指标计算得到的比值，因此同一种材料不同比值之间的量纲不同，需要通过无量纲化处理方法进行处理来消除量纲对最终计算的目标值的影响。其中，修复长效性指标为效益值指标，其越大越好，归一化方法可以采用下述公式：
[0146][0147]
由x’i,j构成的矩阵为第一目标矩阵x’，第一目标矩阵可以表示为：
[0148][0149]
指标变异性以标准差的形式来表现：
[0150][0151]
其中，sj表示第j个指标。
[0152]
指标冲突性用相关系数进行表示：
[0153][0154]rj,j’表示评价指标j和j’之间的相关系数。
[0155][0156]
最终的信息量计算与客观赋权：
[0157][0158]
其中，cj越大，第j个评价指标在整个修复长效性指标的初始值评价体系中的作用越大，就应该分配更多的权重。所以第j个指标的客观权重wj为：
[0159][0160]
各修复长效性指标的初始值的第一权重向量为：
[0161]
w＝[w1,w2,...,wj]
t
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式(1-16)
[0162]
s808，根据第一权重向量及第一目标矩阵确定各初始稳定化材料的修复长效性指
标的目标值。
[0163]
各修复长效性指标的初始值的权重由第一目标矩阵x’与权重阵w相乘得到，即可以用下式来表示：
[0164]bi,4
＝x'
·wꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式(1-17)
[0165]bi,4
为行数为i的列向量，该列向量中每一行的值都对应着一个初始稳定化材料的修复长效性指标的目标值。
[0166]
在本实施例中，通过先进行归一化处理，再计算各修复长效性指标的第一权重向量，最终得到修复长效性指标的目标值，由于是基于各修复长效性指标的对比强度和指标之间的冲突性来综合衡量指标的客观权重，考虑了指标变异性大小的同时兼顾指标之间的相关性，因此根据最终计算的第一权重向量及第一目标矩阵确定的修复长效性指标的目标值比较准确。
[0167]
以上实施例提到了确定初始稳定化材料的综合系数，下面就对其具体确定方法进行介绍。在一个实施例中，如图9所示，根据各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标及土壤回用用途，确定初始稳定化材料的综合系数，包括：
[0168]
s902，根据土壤回用用途，建立关于各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标的判断矩阵。
[0169]
判断矩阵用于表征材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标中相邻两个指标的重要程度。
[0170]
采用ahp法(层次分析法，analytic hierarchy process)在评价指标层次结构模型的基础上，两两比较每种初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、各环境指标之间的相对重要性，用数字1～9作为标度对评价指标的重要性赋值，得到判断矩阵a。由于不同回用用途下，各指标的重要性不同，因此不同回用用途下，判断矩阵也不同。如对于卫生填埋或安全填埋，修复效果较为重要；对于回填用土与绿化用土，修复长效性较为重要。通过专家的专业判断，基于对各指标两两之间的相对重要性(即标度值)的比较构造判断矩阵。
[0171]
判断矩阵可以表示为：
[0172][0173]
其中，n表示第n个指标。其中a
ij
表示ai与aj相比，ai相对于aj的重要性程度，用数字1～9作为标度来描述，标度含义如表3所示。
[0174]
表2 1～9标度值含义
[0175]
标度含义1ai与aj相比，重要性相同3ai与aj相比，ai比aj稍微重要5ai与aj相比，ai比aj明显重要7ai与aj相比，ai比aj强烈重要9ai与aj相比，ai比aj极端重要
2,4,6,8分别表示上述相邻判断的中值倒数若ai与aj的重要性之比为aij，则aj与ai的重要性之比aji＝1/aij
[0176]
s904，计算判断矩阵对应的初始权重向量及最大特征根。
[0177]
根据上述判断矩阵计算对应的初始权重向量及最大特征根的方法如下，首先判断矩阵、权重向量及最大特征根满足下式：
[0178]
aw＝λ
max
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式(1-19)
[0179]
特征向量w即为判断矩阵中各元素对应的权重，求解特征向量的过程可以采用近似解法、方根法、和积法等方法，本技术采用方根法进行计算，具体如下：
[0180]
第一步，将判断矩阵a各行元素相乘；
[0181][0182]
第二步，计算mi的n次方根，并基于各方根建立方根向量；
[0183][0184][0185]
第三步，对方根向量归一化处理；
[0186][0187]
根据归一化后的wi即可得到特征向量w＝(w1,w2,w3,...,wn)
t
，即初始权重向量。
[0188]
第四步，计算判断矩阵的最大特征根λ
max
。
[0189][0190]
其中(aw)i表示向量aw的第i个元素，向量aw为：
[0191][0192]
s906，根据最大特征根对判断矩阵进行修正，直至判断矩阵的一致性满足预设条件，根据当前判断矩阵确定第二权重向量；
[0193]
由于判断矩阵是根据评价者的主观判断确定的，可能存在过分偏离一致性的情况，最终使得计算出的第二权重向量存在较大偏差。因此需要对判断矩阵进行一致性检验方法，以保证判断矩阵满足预设条件。以满足预设条件时的判断矩阵为基础确定第二权重向量。
[0194]
具体方法如下：
[0195]
第一步，计算一致性指标ci；
[0196][0197]
第二步，计算一致性比率cr；
[0198]
[0199]
其中，ri为平均随机一致性指标，其数值与判断矩阵a的阶数有关，取值见下表：
[0200]
表3 ri标准值
[0201]
阶数123456789ri000.580.901.121.241.321.411.45
[0202]
第三步，对判断矩阵的一致性进行判断；
[0203]
当cr《0.1时，认为判断矩阵具有满意的一致性，此时直接将权重向量作为第二权重向量。若判断矩阵的一致性比率cr不满足预设条件，即不满足cr《0.1，那么对判断矩阵进行修正，直至其一致性满足预设条件，根据此修正后判断矩阵计算第二权重向量。
[0204]
s908，根据第二权重向量以及各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标，确定初始稳定化材料的综合系数。
[0205]
针对各初始稳定化材料，在获取第二权重向量之后，由于材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标的重要程度也不同，因此还需要对其进行归一化处理，最终基于归一化处理结果及第二权重向量确定初始稳定化材料的综合系数。
[0206]
在本实施例中，由于针对不同的土壤回用用途，建立各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标的判断矩阵，且通过一致性检验对判断矩阵进行检验，因此保证了建立的判断矩阵的准确性。进而也就保证了基于准确的判断矩阵得到的第二权重向量的准确性，那么最终计算的初始稳定化材料的综合系数也是准确的。
[0207]
在一个实施例中，如图10所示，根据第二权重向量以及各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标，确定初始稳定化材料的综合系数，包括：
[0208]
s1002，对各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标进行归一化处理，得到各指标的归一化处理结果。
[0209]
具体的，为消除各指标对最终计算的初始稳定化材料的综合系数的影响，需要对各指标进行无量纲化处理处理。根据不同的指标特性，需选择不同的归一化方式。
[0210]
其中，修复效果参数b
i,2
及修复长效性指标的目标值b
i,4
为效益型指标，采用下述公式进行计算：
[0211][0212]
材料总成本b
i,1
，修复周期b
i,3
，二氧化碳b
i,5
，颗粒物b
i,6
、固体废物b
i,7
、废水产量b
i,8
为成本型指标，其值越小越好，故使用下式进行归一化处理：
[0213][0214]
针对每种稳定材料都可以得到多个归一化处理结果。
[0215]
s1004，基于归一化处理结果建立第二目标矩阵。
[0216]
基于以上多个归一化处理结果，可以建立如下所示的第二目标矩阵：
[0217][0218]
s1006，根据第二目标矩阵及第二权重向量，计算各初始稳定化材料的综合系数。
[0219]
将第二目标矩阵与第二权重向量相乘，可以得到一个行数为i的列向量，该列向量中的每一个值即为各初始稳定化材料的综合系数，可以表示为：
[0220]
(s1...si)＝y
·wꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式(1-31)
[0221]
在本实施例中，由于是将各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标进行归一化处理，且基于归一化处理的结果建立第二目标矩阵，在建立过程中消除了量纲对各初始稳定化材料的综合系数的影响，因此最终得到的各初始稳定化材料的综合系数比较准确，且由于考虑了成本指标、技术指标及环境指标，因此各初始稳定化材料的综合系数是比较全面的。
[0222]
以上实施例提到了从各初始稳定化材料中确定目标稳定化材料，下面就对其具体方法进行介绍。在一个实施例中，如图11所示，根据初始稳定化材料的综合系数，从各初始稳定化材料中确定目标稳定化材料，包括：
[0223]
s1102，根据初始稳定化材料的综合系数对各初始稳定化材料进行排序，得到排序结果。
[0224]
具体的，可以将各初始稳定化材料的综合系数si从大到小进行排序，得到排序结果，也可以从小到大进行排序，当然也可以采取其他方式，本技术对此不做限定。
[0225]
s1104，从排序结果中确定候选稳定化材料。
[0226]
候选稳定化材料的综合系数在排序结果中位于前预设名次，例如当各初始稳定化材料的综合系数si从大到小进行排序时，可以将排序结果中前三名对应的初始稳定化材料作为候选稳定化材料。
[0227]
s1106，基于候选稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标，从各候选稳定化材料中确定目标稳定化材料。
[0228]
具体的，在确定候选稳定化材料后，获取该稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标，并对以上指标综合考虑从而确定目标稳定化材料。例如，可以针对候选稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标进行分别评分，然后将各评分的平均值作为各候选稳定化材料总评分，最终基于总评分从各候选稳定化材料中确定目标稳定化材料，当然也可以采用其他方法确定目标稳定化材料。
[0229]
在本实施例中，在得到初始稳定化材料的综合系数后，只需要采用简单的排序方法就可以选出候选稳定化材料，该方法比较简单，易于操作。同时，在确定初始稳定化材料后，由于还需考虑候选稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标选取目标稳定化材料，由于选取目标稳定化材料综合了各个指标的结果，因此采用得到的目标稳定化材料比较合适，采用该材料对待修复土壤进行修复时效果也会较好。
[0230]
为使本技术中的方案更加清楚，下面以一个具体的实施例介绍上述稳定化材料选择方法。上述稳定化材料选择方法，包括：
[0231]
第一步，获取待修复土壤的基本信息及修复长效性指标的初始值。
[0232]
基本信息至少包括超标污染物类型及土壤回用用途。
[0233]
第二步，根据超标污染物类型确定至少一种初始稳定化材料。
[0234]
第三步，针对各初始稳定化材料，获取初始稳定化材料的材料总成本。
[0235]
材料总成本包括材料成本及运行成本。
[0236]
第四步，在所述修复长效性指标处于初始值的情况下，确定所述待修复土壤的第一超标污染物浓度及向所述待修复土壤分别添加各所述初始稳定化材料后的第二超标污染物浓度，将所述第二超标污染物浓度与所述第一超标污染物浓度的第一比值作为各所述初始稳定化材料的修复效果参数。
[0237]
修复长效性指标至少包括土壤ph、土壤有机质含量、土壤粘粒含量、当地年均最低温度，年均降雨量，冻融循环次数、二氧化碳浓度、紫外辐射强度。
[0238]
第五步，根据达到各初始稳定化材料的修复效果参数对应的天数，确定各初始稳定化材料的修复周期。
[0239]
第六步，根据修复长效性指标的初始值，计算修复长效性指标的中间值。
[0240]
第七步，在所述修复长效性指标处于中间值的情况下，确定向所述待修复土壤分别添加各所述初始稳定化材料后的第三超标污染物浓度。
[0241]
第八步，针对各所述初始稳定化材料，确定所述第二超标污染物浓度及各所述第三超标污染物浓度的第二比值。
[0242]
第九步，针对各初始稳定化材料，基于各第二比值建立原始矩阵。
[0243]
第十步，对原始矩阵进行归一化处理，生成第一目标矩阵，并基于原始矩阵及第一目标矩阵确定所述修复长效性指标的第一权重向量。
[0244]
第十一步，根据所述第一权重向量及所述第一目标矩阵，确定各所述初始稳定化材料的所述修复长效性指标的目标值。
[0245]
第十二步，针对各初始稳定化材料，获取初始稳定化材料的环境指标。
[0246]
环境指标包括二氧化碳产量、颗粒物产量、固体废弃物产量、废水产量中的至少一种。
[0247]
第十三步，根据土壤回用用途，建立关于各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标的判断矩阵。
[0248]
判断矩阵用于表征材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标中相邻两个指标的重要程度。
[0249]
第十四步，计算判断矩阵对应的初始权重向量及最大特征根。
[0250]
第十五步，根据最大特征根对判断矩阵进行修正，直至判断矩阵的一致性满足预设条件，根据当前判断矩阵确定第二权重向量。
[0251]
第十六步，对各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标进行归一化处理，得到各指标的归一化处理结果。
[0252]
第十七步，基于归一化处理结果建立第二目标矩阵。
[0253]
第十八步，根据第二目标矩阵及第二权重向量，计算各初始稳定化材料的综合系数。
[0254]
第十九步，根据初始稳定化材料的综合系数对各初始稳定化材料进行排序，得到
排序结果。
[0255]
第二十步，从排序结果中确定候选稳定化材料；候选稳定化材料的综合系数在排序结果中位于前预设名次。
[0256]
第二十一步，基于候选稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标，从各候选稳定化材料中确定目标稳定化材料。
[0257]
以下以一个简要的流程图对本技术中的稳定化材料选择方法进行介绍。如图12所示，本技术中的稳定化材料选择方法，包括以下步骤：s1202，获取待修复土壤的基本信息及修复长效性指标的初始值；s1204，根据基本信息中的超标污染物类型选择初始稳定化材料，并针对各初始稳定化材料对待修复土壤进行试验或通过收集资料，得到各初始稳定化材料的经济指标、技术指标相关的值；s1206，采用critic法对各初始稳定化材料的技术指标中的修复长效性指标的目标值进行计算；s1208，采用lca法获取各初始稳定化材料的环境指标；s1210，结合以上各指标值采用ahp法计算不同回用用途下各初始稳定化材料的综合系数；s1212，基于初始稳定化材料的综合系数进行排序，并综合稳定化材料的经济指标、技术指标以及环境指标值筛选确定目标稳定化材料。
[0258]
在本实施例中，通过土壤的超标污染物类型可以先选出适用的初始稳定化材料，再结合稳定化材料的经济指标、技术指标、环境指标及土壤回用用途对初始稳定化材料的综合系数进行计算，而其中技术指标包含了对修复长效性指标的初始值进行更新得到的修复长效性指标的目标值。在确定修复长效性指标的目标值的过程中综合考虑了多个修复长效性指标的权重，最终得出准确的修复长效性指标的目标值。同时在计算综合系数的过程中有针对性的对不同回用用途下的经济指标、技术指标、环境指标的权重进行计算，同时还对各初始稳定化材料的各经济指标、技术指标、环境指标进行归一化处理，因此基于权重及归一化处理结果得到的各初始稳定化材料的综合系数比较准确。且由于考虑了多方面的影响，上述稳定化材料选择方法能够比较综合且全面的对初始稳定化材料进行评价，进而计算出各初始稳定化材料的综合系数。进一步的，得到综合系数后只需要采用简单的排序方法就可以选出候选稳定化材料，该方法比较简单，易于操作。同时，在确定初始稳定化材料后，由于还需结合候选稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标选取目标稳定化材料，因此采用得到的目标稳定化材料比较合适，采用该材料对待修复土壤进行修复时效果也会较好。
[0259]
应该理解的是，虽然如上所述的各实施例所涉及的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，如上所述的各实施例所涉及的流程图中的至少一部分步骤可以包括多个步骤或者多个阶段，这些步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤中的步骤或者阶段的至少一部分轮流或者交替地执行。
[0260]
基于同样的发明构思，本技术实施例还提供了一种用于实现上述所涉及的稳定化材料选择方法的稳定化材料选择装置。该装置所提供的解决问题的实现方案与上述方法中所记载的实现方案相似，故下面所提供的一个或多个稳定化材料选择装置实施例中的具体限定可以参见上文中对于稳定化材料选择方法的限定，在此不再赘述。
[0261]
在一个实施例中，如图13所示，提供了一种稳定化材料选择装置1300，包括：获取模块1320、初始材料确定模块1340、系数确定模块1360、目标材料确定模块1380：
[0262]
获取模块1320，用于获取待修复土壤的基本信息及修复长效性指标的初始值；基本信息至少包括超标污染物类型及土壤回用用途；
[0263]
初始材料确定模块1340，根据超标污染物类型确定至少一种初始稳定化材料；
[0264]
系数确定模块1360，针对各初始稳定化材料，根据初始稳定化材料的经济指标、技术指标、环境指标及土壤回用用途，确定初始稳定化材料的综合系数；技术指标包括对修复长效性指标的初始值进行更新得到的修复长效性指标的目标值；
[0265]
目标材料确定模块1380，根据初始稳定化材料的综合系数，从各初始稳定化材料中确定目标稳定化材料。
[0266]
在另一个实施例提供的稳定化材料选择装置1300中，经济指标包括材料总成本，技术指标包括修复效果参数、修复周期以及修复长效性指标的目标值，系数确定模块1360包括：
[0267]
第一获取子模块，用于针对各初始稳定化材料，获取初始稳定化材料的材料总成本；材料总成本包括材料成本及运行成本；
[0268]
第二获取子模块，用于针对各初始稳定化材料，获取初始稳定化材料的修复效果参数、修复周期以及修复长效性指标的目标值；
[0269]
第三获取子模块，用于针对各初始稳定化材料，获取初始稳定化材料的环境指标；环境指标包括二氧化碳产量、颗粒物产量、固体废弃物产量、废水产量中的至少一种；
[0270]
系数确定子模块，用于根据各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标及土壤回用用途，确定初始稳定化材料的综合系数。
[0271]
在另一个实施例提供的稳定化材料选择装置1300中，第二获取子模块包括：
[0272]
效果参数确定单元，用于在所述修复长效性指标处于初始值的情况下，确定所述待修复土壤的第一超标污染物浓度及向所述待修复土壤分别添加各所述初始稳定化材料后的第二超标污染物浓度，将所述第二超标污染物浓度与所述第一超标污染物浓度的第一比值作为各所述初始稳定化材料的修复效果参数；
[0273]
可选的，修复长效性指标至少包括土壤ph、土壤有机质含量、土壤粘粒含量、当地年均最低温度，年均降雨量，冻融循环次数、二氧化碳浓度、紫外辐射强度。
[0274]
修复周期确定单元，用于根据达到各初始稳定化材料的修复效果参数对应的天数，确定各初始稳定化材料的修复周期；
[0275]
目标值确定单元，获取向所述待修复土壤分别添加各所述初始稳定化材料后的第三超标污染物浓度，根据修复长效性指标的初始值、所述第二超标污染物浓度及所述第三超标污染物浓度，确定各所述初始稳定化材料的修复长效性指标的目标值；所述第三超标污染物浓度与所述修复长效性指标的中间值相对应。
[0276]
在另一个实施例提供的稳定化材料选择装置1300中，目标值确定单元具体用于：根据修复长效性指标的初始值，计算修复长效性指标的中间值；在修复长效性指标处于中间值的情况下，确定向待修复土壤分别添加各初始稳定化材料后的第三超标污染物浓度；根据修复长效性指标的初始值、第二超标污染物浓度及第三超标污染物浓度，确定各初始
稳定化材料的修复长效性指标的目标值。以及用于针对各初始稳定化材料，确定各修复长效性指标的初始值对应的第二超标污染物浓度及第三超标污染物浓度的第二比值；针对各初始稳定化材料，基于各第二比值建立原始矩阵；对原始矩阵进行归一化处理，生成第一目标矩阵，并基于原始矩阵及第一目标矩阵确定各修复长效性指标的第一权重向量；根据第一权重向量及第一目标矩阵确定各初始稳定化材料的修复长效性指标的目标值。
[0277]
在另一个实施例提供的稳定化材料选择装置1300中，系数确定子模块还包括：
[0278]
矩阵建立单元，用于根据土壤回用用途，建立关于各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标的判断矩阵；判断矩阵用于表征材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标中相邻两个指标的重要程度；
[0279]
计算单元，用于计算判断矩阵对应的初始权重向量及最大特征根；
[0280]
修正单元，用于根据最大特征根对判断矩阵进行修正，直至判断矩阵的一致性满足预设条件，根据当前判断矩阵确定第二权重向量；
[0281]
综合系数确定单元，用于根据第二权重向量以及各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标，确定初始稳定化材料的综合系数。
[0282]
在另一个实施例提供的稳定化材料选择装置1300中，综合系数确定单元具体用于对各初始稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标进行归一化处理，得到各指标的归一化处理结果；基于归一化处理结果建立第二目标矩阵；根据第二目标矩阵及第二权重向量，计算各初始稳定化材料的综合系数。
[0283]
在另一个实施例提供的稳定化材料选择装置1300中，目标材料确定模块包括：
[0284]
排序子模块，用于根据初始稳定化材料的综合系数对各初始稳定化材料进行排序，得到排序结果；
[0285]
第一确定子模块，用于从排序结果中确定候选稳定化材料；候选稳定化材料的综合系数在排序结果中位于前预设名次；
[0286]
第二确定子模块，用于基于候选稳定化材料的材料总成本、修复效果参数、修复周期、修复长效性指标的目标值、环境指标，从各候选稳定化材料中确定目标稳定化材料。
[0287]
上述稳定化材料选择装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
[0288]
在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图14所示。该计算机设备包括通过系统总线连接的处理器、存储器、通信接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的通信接口用于与外部的终端进行有线或无线方式的通信，无线方式可通过wifi、移动蜂窝网络、nfc(近场通信)或其他技术实现。该计算机程序被处理器执行时以实现一种稳定化材料选择方法方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外
壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
[0289]
本领域技术人员可以理解，图14中示出的结构，仅仅是与本技术方案相关的部分结构的框图，并不构成对本技术方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
[0290]
在一个实施例中，还提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现上述稳定化材料选择方法各实施例中的步骤。
[0291]
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述稳定化材料选择方法各实施例中的步骤。
[0292]
在一个实施例中，提供了一种计算机程序产品，包括计算机程序，该计算机程序被处理器执行时实现上述稳定化材料选择方法各实施例中的步骤。
[0293]
需要说明的是，本技术所涉及的数据(包括但不限于用于分析的数据、存储的数据、展示的数据等)，均为经用户授权或者经过各方充分授权的信息和数据。
[0294]
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本技术所提供的各实施例中所使用的对存储器、数据库或其它介质的任何引用，均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(read-only memory，rom)、磁带、软盘、闪存、光存储器、高密度嵌入式非易失性存储器、阻变存储器(reram)、磁变存储器(magnetoresistive random access memory，mram)、铁电存储器(ferroelectric random access memory，fram)、相变存储器(phase change memory，pcm)、石墨烯存储器等。易失性存储器可包括随机存取存储器(random access memory，ram)或外部高速缓冲存储器等。作为说明而非局限，ram可以是多种形式，比如静态随机存取存储器(static random access memory，sram)或动态随机存取存储器(dynamic random access memory，dram)等。本技术所提供的各实施例中所涉及的数据库可包括关系型数据库和非关系型数据库中至少一种。非关系型数据库可包括基于区块链的分布式数据库等，不限于此。本技术所提供的各实施例中所涉及的处理器可为通用处理器、中央处理器、图形处理器、数字信号处理器、可编程逻辑器、基于量子计算的数据处理逻辑器等，不限于此。
[0295]
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
[0296]
以上所述实施例仅表达了本技术的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本技术专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本技术构思的前提下，还可以做出若干变形和改进，这些都属于本技术的保护范围。因此，本技术的保护范围应以所附权利要求为准。