研究区域为我国东南部某县级市。该县级市总面积超过1$\text{×}$10⁶ m²，常住人口超过1$\text{×}$10⁶，该地包含6个街道和8个镇，基础设施建立完备，交通系统稠密完善。该县级市已实施垃圾分类政策，但在垃圾处理方面的人力资金投入与需求相比仍然不足，要求后续垃圾收运工作充分考虑到生活垃圾分类方法。此外，该县级市内的野生动植物资源丰富，应在进行市政规划管理时给予保护意识与足够重视。

研究对象为该县级市的其他垃圾，其处理方式为：先从垃圾收集点由垃圾收集车收集至经垃圾中转站进行压缩预处理后，然后由垃圾运输车运输至垃圾焚烧厂。经实际调研与网络查询，当地具有329个垃圾收集点、40个垃圾中转点、3个垃圾焚烧厂。

本研究所需路网分布、道路长度、道路等级、道路限速等数据信息来源于Open Street Map网站的2020年全国路网矢量图^[8]，经过裁剪后得到研究区域路网图如图1所示。环境成本中包含的土地保护需求指标数据来源于500 m全球土地退化态势评价数据集^[9]，影响人口规模指标的数据来源于2020年中国人口密度分布图^[10]，降噪物密度指标的数据来源于中国-东盟1 km分辨率植被覆盖度数据集^[11]，生态不可分割度指标的数据来源于2020年中国土地利用图^[12]。

本研究将传统蚁群算法中的“往返时间”概念扩大至“成本”，在ArcGIS 10.7软件中的Python拓展模块实现蚁群算法，其步骤为：1)首先，进行路网处理，将道路以节点为界限划分路段，且对节点进行编号；2)之后，在Python拓展模块中编写程序，程序中对传统蚁群算法中的期望函数与信息素更新规则进行改进，以体现道路对生态环境的恶劣影响与改进经济成本代表指标。

基本蚁群算法涉及到2个过程。第一个过程是状态转移，蚁群算法中的状态转移概率函数由节点间距离与信息素浓度共同决定。i及j节点之间状态转移概率函数P_ij^m的计算方式如式(1)所示。

式中：m为蚂蚁数量；α为信息启发因子，表示信息素浓度的重要性，取值为0~5；β为期望启发因子，表示节点之间成本的重要性，取值为0~5；τ_ij为i、j节点之间的信息素浓度；η_ij为从i节点到j节点的期望值。传统蚁群算法中期望函数η_ij的计算方式如式(2)所示。

第二个过程是信息素更新，(t+1)时刻i、j节点之间路段上信息素浓度τ_ij(t+1)的计算方式如式(3)、式(4)所示。

式中：ρ为信息素挥发系数；1-ρ为信息素残留系数；Δτ_ij(t)为t时刻所有蚂蚁在节点i、j间留下的信息素浓度；Δt_ij^k(t)为t时刻第k只蚂蚁在节点i、j间带来的信息素增量。

信息素增量的计算方式有3种：蚁周模型、蚁量模型、蚁密模型。本研究选择不容易陷入局部最优的蚁周模型为信息素更新规则。传统蚁群算法中的信息素更新规则如式(5)所示。

式中：Q为信息素总量，取值范围为10~10 000；l_ij为i、j节点之间的道路长度。

本研究利用 SPSS 软件实现 CRITIC 法与熵权 TOPSIS 法对综合成本指标的权重分配，其步骤为：1)建立综合成本评价指标体系；2)获取评价所需数据并建立道路数据集；3)以 CRITIC 法对环境成本指标下各子指标进行权重分配；4)以熵权 TOPSIS 法对综合成本指标下的时间成本与环境成本进行权重分配，环境成本权重与时间成本权重的比值即为新定义的环境干扰因子。

客观权重赋值法与以层次分析法等为代表的主管权重赋值法相比，只关注指标数据之间的差异，仅从数据中生成属性权重，而不需要从使用者中获取任何偏好信息^[13]，也可有效降低决策方法的复杂性与非直观性^[14]。其中，CRITIC算法考虑了各指标自身的对比强度及指标间的冲突性，能够较全面的衡量各指标重要性，因而被作为一种相对完善赋权算法，被广泛使用。熵权法可用于分析对象之间存在较大差异的指标。它们均可用于整合社会影响评估和不同标准的确定。

本研究利用客观权重赋值法中的CRITIC权重法与熵权TOPSIS 法建立道路综合成本评价体系，在体系中纳入道路交通对自然生态、人居环境的影响，并将该体系应用于所研究县级市路网中，以获得各路段的综合成本数值，为下一步计算得到综合成本最低的垃圾收运路径方案做好数据准备。

首先，考虑到数据的可达性与时效性，本研究取道路长度与限速之比为时间成本；环境成本由土地保护需求、影响人口规模、降噪物密度、生态不可分割度确定。

本研究使用时间成本与环境成本作为综合成本C_ij的2大指标，其计算方式如式(6)所示。

式中：t_ij为i、j节点之间道路的时间成本；γ为环境干扰系数；e_ij为i、j节点之间道路的环境成本。

时间成本t_ij的计算方式如式(7)所示。

式中：l_ij为i、j节点之间道路的长度，km；v_ij为i、j节点之间的车速，km∙h⁻¹，由道路限速与车辆最高行驶速度确定。

环境成本e_ij的计算方式如式(8)所示。

式中：λ₁为土地退化态势的权重，0~1；e_ij,1为i、j节点之间道路穿行的土地退化态势；λ₂为人口密度的权重，0~1；e_ij,2为i、j节点之间道路穿行地区的人口密度；λ₃为植被覆盖度的权重，0~1；e_ij,3为i、j节点之间道路穿行地区的植被覆盖度；λ₄为生态抗扰度的权重，0~1；e_ij,4为i、j节点之间道路穿行地区的生态抗扰度。

利用CRITIC法赋予环境成本中各指标权重，并利用熵权法计算环境干扰因子。时间成本由道路长度与道路限速之比得到。之后，编辑路网数据集，添加“时间”字段，通过字段计算器算的每段道路通过所需花费的时间。

垃圾运输对运输区域的生态干扰程度主要由土地保护需求、影响人口规模以及降噪物密度决定。本研究利用土地退化态势来确定土地保护需求。土地退化态势数据集中，将全球的土地退化态势分为2级评价系统。在第1级，将土地退化态势评价划分为3种类型：退化类型、改善类型和无变化类型。第2级，将土地退化态势评价划分为9种类型。下载土地退化态势数据集后，在ArcGIS软件中打开，并根据研究区域行政区划线要素矢量数据图对其裁剪，得到如图2(a)所示的研究区域土地退化态势分布栅格图，图中颜色越浅表示土地退化越严重。

影响人口规模由人口密度分布确定。将人口密度分布栅格数据经过裁剪与符号化后标注于研究区域地图，通过识别功能得到各垃圾中转站所处位置的人口密度。研究区域人口密度栅格图像如图2(b)所示。图中颜色越深，表明人口密度越大。

降噪物密度由植被覆盖度确定。植被覆盖度数据的处理方式同土地退化态势分布图，结果如图2(c)所示。图中颜色越深，表示植被覆盖度越大。

在环境成本保护的各指标中，生态不可分割度为本研究根据生物丰度指数含义进行改进得到的指标，旨在评价道路交通沿线周边生态系统对道路切割的抵抗力，以减少堵路交通对沿线周边造成的不可逆性深远影响。生态不可分割度越高，往往当地生态系统所含的食物链就越复杂，则道路交通对当地生态环境切割与生物栖息地破碎化的负面影响就越深。根据我国生态环境部颁布的《生态环境状况评价技术规范(HJ 192-2015)》^[15]的标准，定义生态抗扰度的打分方式如表1所示。利用空间地理数据云的2020年中国土地利用图矢量数据计算生态不可分割度，经过裁剪，其结果如图3所示。图中颜色越接近红色，城市化程度越高，生态不可分割度得分越低。

在环境成本的所有指标中，土地利用类型为矢量数据，其匹配方式为，选中路网文件，选择“连接和关联”中的“连接”，并基于空间位置来连接土地利用类型图层。其余的指标均为栅格数据，只需要在已有的路网矢量数据添加上相对应的栅格数据值即可。至此，研究区域各路段已被赋予环境成本数值。之后，在SPSSAU平台上利用CRITIC权重法计算环境成本中各指标权重。计算结果如表2所示。据此可得到各路段取值为0~1的环境成本数值。

对于环境干扰因子的取值，首先使用熵权法对时间成本与环境成本的权重进行赋值，其次取环境干扰因子的值为环境成本与时间成本权重的比值。在SPSSAU平台上进行熵权法计算后，结果如表3所示：

环境干扰因子取值为0.1116。据此计算研究区域路网中所有路段的综合成本，之后在ArcGIS里的路段属性表中添加综合成本字段，完善路网数据集。

为考察综合成本对垃圾收运路径规划的影响，本研究先基于传统蚁群算法进行垃圾收运路径规划设计，然后与以改进蚁群算法为原理的规划路径作对比分析。传统蚁群算法中将“距离”作为经济成本的代替数据，以路径最短作为计算目标。

本研究中蚁群算法的实现依托于ArcGIS10.7 软件中的Python拓展模块。本研究区域的垃圾收集点与垃圾中转点数较多，所以用PyPy代替Python中的解释器CPython，将源代码直接转换为机器汇编语言，以大幅提升程序运行效率。程序运行后，输出最佳路径的节点编号。将节点选择结果标记至ArcGIS 10.7中的路网图上，得到基于传统蚁群算法的综合成本最低的收集路径、运输路径。结果如图4(a)、(b)所示。

为了修正传统蚁群算法对计算目标设定与指标选取的局限性，本研究以改进传统蚁群算法中期望函数定义与信息素浓度更新规则的方式，在蚁群算法计算最优路径时引入时间成本与环境成本，以改正原路径计算方法用“距离”代表经济成本、忽略环境影响的缺陷。其中，改进蚁群算法中期望函数的计算方式如式(9)所示。

式中：c_ij为i、j节点之间的综合成本。

本研究选择不容易陷入局部最优的蚁周模型为信息素更新规则，并结合将“成本”代替“路径长度”的理念，规定信息素更新规则为信息素总量与成本之比。其计算方式如式(10)所示。

式中：Q为信息素总量；C_k为蚂蚁k在本次遍历中所走路径小号的综合成本。

综合成本C_k的计算方式如式(11)所示。

式中：T_k为蚂蚁k在本次遍历中所走路径消耗的时间成本；E_k为蚂蚁k在本次遍历中所走路径消耗的环境成本。

实现改进蚁群算法需要基于路段的综合成本属性，所以需要构建新的网络数据集。与传统蚁群算法相比，改进蚁群算法需要在新的数据集指定属性时添加综合成本，并在配置网络行驶方向时设置阻抗为综合成本，其余步骤与之前相同。程序运行后，根据运算所得到的节点选择相应结果，标记路网地图，得到基于改进蚁群算法的综合成本最低的收集路径与运输路径，结果如图5(a)、(b)所示。从图可以看出，基于改进蚁群算法得到的垃圾收集路径与垃圾运输路径均与基于传统蚁群算法计算得到的结果有明显不同。

为进一步分析利用 CRITIC 法与熵权 TOPSIS法构建的路径综合成本评价体系对垃圾收运路径规划的影响，本节分别以环境成本中纳入的4项环境成本指标数据图为底图，用分辨率为1 km×1 km 的表格划分研究区域，定义被规划路径穿行的敏感区域(后面将分别定义)所占格子个数与规划路径穿行格子总数之比为规划路径对敏感地区的穿行率。计算基于2种蚁群算法算出的垃圾收运路径规划结果的穿行率并进行对比，讨论考虑环境成本后的最优垃圾收运路径与仅考虑经济成本的结果有何不同。如图4和图5所示，红色路径均为基于传统蚁群算法求得的最优垃圾收运路径，蓝色路径均为基于改良蚁群算法求得的最优垃圾收运路径。

将基于传统蚁群算法与改进蚁群算法得到的最优垃圾收集路径、运输路径均标记在土地退化态势图中，结果如图6(a)所示。底图的颜色越浅，表示土地退化越严重。定义土地退化态势为负数的区域为土地保护需求强烈区域(敏感区域)，则分析图6(a)可知，红线对土地保护需求强烈区域的穿行率为81%，蓝线对土地保护需求强烈区域的穿行率为48%，蓝线较红线而言土地保护需求强烈区域穿行率减少了41%。可见，红线比蓝线更集中于土地退化严重的区域，而蓝线多分布于土地质量正在改善或不退化的区域。将土地保护需求纳入垃圾收运路径规划后，可有效缓解垃圾收运交通对土地退化的进一步恶劣影响，从而促进了退化土地恢复与土地质量改善。

将基于传统蚁群算法与改进蚁群算法得到的最优垃圾收集路径、运输路径均标记在人口密度分布图中，结果如图6(b)所示。底图的颜色越深，表示人口密度越大。定义人口密度大于中位数 311 人∙km⁻²的区域为影响人口规模偏大区域(敏感区域)。分析图6(b)可知，红线对影响人口规模偏大区域的穿行率为2%，蓝线对影响人口规模偏大区域的穿行率为18%。与红线相比，蓝线影响人口规模偏大区域穿行率减少了44%。可见，红色线条比蓝色线条更集中于人口密集区域，而蓝色线条多分布于人口密度较小的地方，整体布局也更加均匀。人口规模大的区域往往交通繁忙，将影响人口规模纳入垃圾收运路径规划后，不仅能绕开人口密度大的居住地从而改善人居环境，还能有效缓解拥堵路段，从而缩减垃圾收运时间、提高垃圾收运效率、减少车辆运行造成的环境污染。

将基于传统蚁群算法与改进蚁群算法得到的最优垃圾收集路径、运输路径均标记在植物覆盖度图中，结果如图6(c)所示。底图的颜色越深，表示植物覆盖度越大，对车辆噪音的削减能力越强。定义植物覆盖度小于中位数192的区域为降噪物密度偏小区域(敏感区域)。分析图6(c)可知，红线对降噪物密度偏小区域的穿行率为76%，蓝线对降噪物密度偏小区域的穿行率为54%，蓝线较红线而言降噪物密度偏小区域穿行率减少了9%。可见，红线比蓝线更集中于难以被植物削减噪音的区域，而蓝线多分布于容易被植物削减噪音的区域。将降噪物密度指标纳入垃圾收运路径规划后，可将垃圾收运路径更多转移至植物密度高的地方，使车辆通行产生的噪音得到更大程度的削减，从而减轻了垃圾收运交通噪音对人类的影响，提高了人居环境质量。

将基于传统蚁群算法与改进蚁群算法得到的最优垃圾收集路径、运输路径均标记在土地利用类型图中，结果如图6(d)所示。图中颜色越接近红色，城市化程度越高，生态不可分割度得分越低。定义生态不可分割度大于中位数40的区域为生态不可分割度偏大区域(敏感区域)。分析图6(d)可知，红线对生态不可分割度偏大区域的穿行率为83%，蓝线对生态不可分割度偏大区域的穿行率为49%，蓝线较红线而言生态不可分割度偏大区域穿行率减少了41%。可见，红线比蓝线更集中于不应被人为分割的区域，而蓝线多分布于已进行过人为建设的区域，即红线比蓝线更集中于生态不可分割度高(城市化程度高)的区域。将生态不可分割度纳入垃圾收运路径规划后，可有效减轻垃圾收运交通对周边生态环境的切割程度，减少了人类活动对自然度高的区域的进一步干扰，从而在一定程度上保护了物种多样性与生态稳定性。

1)以“耗时”作为经济成本的衡量指标，将限速情况等路况因素考虑在内，可以充分体现包括路径长度在内的耗油量、人力成本等额外经济成本，从而更利于道路交通经济成本的控制。

2)基于人口密度分布等数据，将道路交通对周边人群健康、生态环境质量的影响纳入对道路交通的综合成本评价当中，对蚁群算法进行了因地制宜的改良，将规划路径对土地保护需求强烈区域、影响人口规模偏大区域、降噪物密度偏小区域、生态不可分割度偏大区域的穿行率分别降低了41%、44%、29%、41%。