基于改进的K-means聚类的多区域物流中心选址算法①
鲁玲岚， 秦江涛
(上海理工大学 管理学院，上海 200093)
摘 要： 针对当前多区域物流中心选址需建立配送中心个数不定、位置、覆盖范围不明的问题，本文提出了一种改进的k-means聚类算法，以城市经济引力模型为基础，将城市运输距离与居民消费能力的指标相结合，重新定义对象之间相似性度量的距离因子. 并将密度思想引入k-means算法，提出类内差分均值的概念确定最优聚类数. 实现分区后，分别在这些区域中利用重心法对配送中心进行最终的确定. 最后实例分析了在西部地区37个城市创建物流配送中心的选址过程，并通过和传统的k-means聚类的选址结果对比，说明改进后的算法不仅可以节省配送时间，而且大大降低了运输成本，有很好的经济利用价值.
关键词： 多区域配送中心选址; k-means聚类; 城市经济引力模型; 重心法; 西北物流
区域配送中心配送规模较大，用户较多，是物流中重要节点. 其合理的物流规划不仅关系到是否能达到资源利用的效益最大化，满足客户的物流需求，保证物流时效性，也决定着投入物流成本的高低. 多区域物流配送中心选址需确定配送中心最佳数目、位置及覆盖范围. 聚类算法将一整个区域划分，以此来确定配送中心数目以及其覆盖范围，非常适用于多配送中心的选址问题. 诸如K-means的聚类算法通常用于多区域配送中心的选址，但它们的初始聚类中心是随机的，易陷入局部极小解，且K值不同导致聚类效果各异[1]，并且有些城市地处偏僻，常被作为噪音数据排除，最主要的是由于山川、河流的阻碍，两地之间的欧式距离并不能用来衡量实际的运输距离.
结合检测报告，并经现场察看，可以把原因大致分为两类。一类是共性问题，比如颤振、河水水位、不良地质；一类是桥梁右侧砂场及砂场重型机械设备个性问题，即堆载，而堆载对桥梁的影响又与地质情况相关。
1 目前主要K-means聚类划分方法综述
K-means聚类算法技术较为成熟，操作便捷，经常被用于选址中对区域的划分上. 考虑到该聚类方法对区域划分时存在的劣势，许多学者对该聚类方法进行了改进. 朱培芬结合密度的思想，在可选范围中优先考虑边缘点，具有良好的全局收敛性[2]，但并未科学衡量距离因子，也未曾考虑物流需求等重要影响因子. 谷炜提出两阶段K-means聚类算法，在传统K-means聚类后，通过不断迭代来检验是否满足配送时间最少的收敛条件，最终确定聚类结果[3]，避免了算法陷入局部最优，但不断的迭代大大增加了算法的运算时间. 于晓寒考虑到河流、公路等地理障碍，以“障碍距离”为差异度量标准，站点工作量为约束，提出基于障碍距离的约束聚类算法[4]，但距离并不是简单的空间距离，虽其在抽象意义上是同质的，但诸如城市这样实质的经济主体，距离明显包含了更多的非空间因素.
本文针对选址的实际情况，以城市经济引力模型为基础，考虑到距离的“非空间因素”，重新定义了对象间度量距离的因子. 同时把密度的思想引入K-means算法以解决K值不确性问题，并提出类内差分均值的概念确定最优聚类数. 借助改进的K-means聚类算法完成西北地区物流网络的划分，实现分区后，分别在这些区域中利用重心法对配送中心进行最终的确定. 并通过和传统K-means聚类结果对比，来验证改进后方法更优.
2 改进的K-means聚类算法和重心法
2.1 改进K-means聚类对区域划分
(1) 将密度思想引入K-means算法
创刊于1946年的《文联》，共发表了121篇文章，文类体裁多样，主题丰富，在不同程度反映了社会的整体面貌：有歌颂战争期间兵民的团结一心、众志成城；有揭露站前战后黑暗腐败的现实和专制独裁的野心；有批判具有劣根性的国民；有鞭挞帝国主义的狼子野心；也有讴歌光明的赞歌等，具有强烈的时代感。由于年代的限定，《文联》是存在缺陷的，像存在对文化市场的迎合，满足底层市民的消极趣味。但不可否认《文联》对于先进思想的传播，对于先进文化的弘扬以及对于文艺创作的繁荣，都产生了积极深远的影响，在现代文学报刊史上留下了浓墨重彩的画卷。
基于传统K-means算法易受噪声和孤立点影响的事实，本文将密度的思想引入K-means算法来确定k个聚类中心. 不仅可以避免噪声数据干扰，而且可以有效降低算法时间复杂度. 考虑到M维空间的n维数据点Xi(i=0,1,2,···,n-1) 其基本思路如下：计算每个数据点Xi的密度. j=0,1,2···n-1. ra为正数，定义为该点的领域半径，取
当i取一个值时，k=0,1,2,···,n-1，Xi的密度指标记为Di具体计算公式如下：
在上一阶段过程中确定了配送中心的个数以及其覆盖范围的问题，但聚类过程中选取的聚类中心仅仅考虑到了城市交通的密度，而未考虑其他因素，不足以认为是最合适的配送中心地址. 需结合单一配送选址模型在各个区域中对聚类中心进行修正，从而确定单一区域配送中心的最佳位置，本文选取重心法来进行此阶段工作.
3.1.1 3种作物光谱位置变量特征分析 从图1可看出，3种作物光谱反射率有明显的区别。为了分析得更具体，对原始反射率数据进行高光谱特征参数提取、植被指数计算等预处理，提取了19个光谱参数。10个位置参数中λr的值差别大，其他值相对差别很小。这说明了不同作物的红边内一阶微分光谱中的最大值Dr对应的波长位置差别大。3种作物中Dr对应的波长位置λr从大到小的顺序为大麦、春小麦、莜麦（表1）。
(2) 基于经济引力模型对距离的定义
在度量两个城市之间相似度的过程中，把距离的概念仅定义为两个城市间的空间距离是不可取的，为了更科学的衡量距离且更加有效的运用于实际情况，本文以居民可支配收入作为衡量城市间的经济引力因子，交通运输时间作为衡量城市间的距离因子. 相应的公式如下：
聚类的目的是使同一类的样本点间相似度高，而不同类间相似度低. 为精确测量聚类结果，进一步确定k值. 本文采用计算类内差分均值的评估方法. 计算样本点与各自聚类中心的距离之和的平均值. 值越小，总的聚类距离越小，类间相似度越高，聚类效果越好; 反之，值越大，聚类距离越大，类间相似度越低，聚类效果越差. 用内类差均值的方法衡量最佳聚类个数，以达到修建最少的配送中心实现最大区域的覆盖.
4)计算
种聚类结果的类内差分均值，选取最优优k值，确定最佳聚类方案.
式(2)中MiMj分别为城市i 和城市j的居民人均可支配收入指标. 居民人均可支配收入是衡量居民消费能力的重要指标，在消费性支出中包括食品烟酒、衣着、居住、生活用品以及其他用品和服务. 随着科技的进步，出现了越来越多样化的消费渠道. 这里主要考虑线上和线下消费. 在物流网络建设不完善的地区，居民线上消费遇到较大限制，消费被迫转到线下，这些地区物流需求远远不及实际的消费需求，用居民的可支配收入来衡量物流需求更为现实. Dij为城市i和城市j的交通距离，本文以车辆在两座城市之间行驶一趟的交通时间来衡量. u作为调节物流网络划分时受地区居民消费水平影响的程度. 一般来说，u越大，则说明在区域物流网络划分时，各地区的居民消费水平的吸引力占主导地位，城市间的距离因素被较大程度的弱化; 反之u越小，则认为各地区的居民消费水平影响较弱，城市间的距离因素主导地位越强; u为0时，距离等价于欧式距离.
2.2 重心法
重心法是一种简单可行的选址方法，通常用于解决连续点的单个配送中心选址问题，其唯一的决策依据是运输成本. 应用时，它对候选位置没有任何限制，在已知各个备选地点的位置、物流需求量、各个备选地点的直线距离的前提下，运用重心法可以很好地确定配送中心的位置. 其运算方便，计算速度快，通过几次迭代就能计算出理论位置. 但此种方法只适用于确定单一配送中心，对于确定多个配送中心的选址问题，此法并不适用[5].
3 多区域配送中心选址模型
针对多区域配送中心选址问题，本文采用目前比较主流的两阶段模型：第一阶段把所有的需求点划分为若干个配送区域，第二步阶段选取相应区域中最佳配送中心[6-8]. 借助该模型的选址步骤，本文首先用改进的k-means聚类对区域进行划分，确定需建立配送中心的个数以及其配送范围. 然后运用重心法在划分的区域中选取合适的配送中心位置.
区域物流中心选址模型是带有复杂约束的非线性规划模型，在构建评估模型时先明确如下假设：(1) 中心容量总是可以满足所覆盖需求点的需求量，并由该中心供应的所有单位的需求量确定. (2) 一个需求点仅由一个物流配送中心供应. (3) 从物流中心到需求点通过零担物流的方式进行配送，且不考虑装载的问题，仅以最大车载量进行运算. (4) 运输费用仅由运价、实际运输路程决定，不考虑装卸等人工成本. (5) 配送中心每辆车每天仅往需求点配送一趟. 基于以上5条假设建立衡量选址是否有效的模型. 目标函数是各个配送中心到需求点的运输费用之和最小，目标函数为：
3.1 基于改进的K-means聚类对区域划分
本文运用MyEclipse Professional 2014软件，采用JAVA语言，在Win10 64位系统环境下运行.其中按照样本大小和聚类个数k值的关系，
这里n为样本的大小. 具体操作步骤如下：
1) 存储数据，对数据进行归一化处理.
2) 按照式(1)计算城市密度，按密度从大到小的顺序排列，取排序前m个城市作为聚类中心候选点.
3) 取前k个候选点作为聚类中心，按照式(2)计算当k取不同值时的距离因子，这里u=0.1，当u=0.1时，最能科学均衡经济引力与距离之间的影响力. 按照最小距离原则分配n-k 个城市样本. 得到
种聚类结果.
计算中导致先验信噪比(增益函数)过高估计的两个主要原因是:一是对先验信噪比的估计直接采用“直接判决”法;另一个原因则归因于噪声谱的估计,计算信噪比往往需要先估计噪声谱,而大都数噪声估计算法都会低估噪声谱,噪声谱的低估导致了信噪比的过高估计。基于此:
(3) k值的确定
5)结束，输出最优聚类结果.
3.2 重心法对单一区域配送中心选址
根据式(1)获取所有样本点的密度指数，按照从大到小的顺序排列，选取前k个数据点作为聚类中心.
4 评估模型建立
如果一名代理人为自己买下一物，根据的是我的委任且以我的名字向他作了交付，dominio即所有权由我获得，哪怕他没有交出它……
式(3)中i表示配送中心，其取值为0,1,2,···,k-1,k表示配送中心个数. j为配送中心所供应需求点个数，Ni表示第i个配送中心需求点集合. mij表示配送中心i到需求点j的计费里程，以km为单位. wij表示配送中心i到需求点j运输货物的重量，以kg为单位. J表示从配送中心i到需求点j的单位运价，以元/kg*km为单位.
5 实例分析
5.1 数据来源与处理
西北地区地域广阔，自然资源丰富，但地区天气恶劣，城市之间道路险阻较多，交通十分不便，物资运输上的不便是西北地区电子商务的发展主要阻碍. 本文对西北物流网络进行研究，由于新疆地区地域物流管理体制较混乱，机构多元化，多采用外包的物流模式，这里主要讨论除新疆之外的37个地级市的物流区域配送中心选址问题. 本次共采集37个城市间交通时间数据共703条记录. 选取甘肃、青海、陕西、宁夏2016年的统计年鉴中选取四省的城镇居民可支配收入. 并对已有数据进行了归一化处理，来消除数据量纲的影响.
为验证方法的正确性和实用性，抽象实验数据为：零担运输的机动车辆车身均为5～7米，车载量20 000 kg，平均行驶速度68 km/h. 考虑到西北地区一体化趋势日益明显，可忽略地域的运价费用差异，这里单位运价均定1元/kg*km. 通过城镇居民可支配收入折算成每日的需求量以及配送所需车辆如表1所示.
“虽然对于伦理道德不曾说及，而那种偶像破坏的精神与力量却是极大，给思想界开了一个透气的孔，这可以算是第一个思想革命家。”〔10〕
表1 城镇每日物流需求量以及所需配送车辆量化表
5.2 改进方法选址求解
(1) 确定聚类中心候选点，本文取前10位城市以及相应密度，从结果看出，嘉峪关、酒泉、西安等作为当前交通的枢纽中心，是聚类中心的首选. 在区域划分中把密度思想引入是合理有效的. 结果如表2如下.
表2 前10城市密度
(2) 通过类内差分均值评估，绘制曲线图观察最适k值. 从折线图很直观看出，当k=6时，曲线出现拐点，达到最小值，当k=7，8时类内差均值小范围增加，证明当建立6个配送中心时物流网络已经达到了较好的覆盖效果. 随着配送中心的增多，覆盖效果反而降低. 即当37个城市被分为6类时，类间的相似度越高，聚类结果越好. 如图1所示.
筛选我院儿科病区自2016年6月至2016年12月收治入院的外感发热患儿,根据中医辨证选择其中符合要求的患儿共189例,随机分为对照组和治疗组。治疗组96例,男50例,女46例,平均年龄(3.9±1.8)岁;对照组93例,男49例,女44例,平均年龄(3.6±1.7)岁。
确定最佳聚类个数后，根据改进的K-means算法聚类，当k=6时，聚类结果如表3所示.
图1 3～8个聚类的类内差分均值
表3 k=6时的城市分区结果
(3) 把每个区域中城市作为物流服务需求点，人均可支配收入作为需求量，用重心法修正的配送中心结果如下，区域一：张掖，区域二：宝鸡，区域三：嘉峪关，区域四：铜川市，区域五：共和县，区域六：白银.配送中心地址及其辐射区域如图2所示.
图2 西北四省配送中心地址及其覆盖城市图
5.3 结果评估
本文参照传统的K-means算法对配送中心选址结果进行对比. 在传统的K-means算法中把城市之间的直线距离作为相似性度量的因子，以每个城市之间的经纬度计算城市距离，设置k=6，阈值为0.000 001，共迭代50次，得到区域划分后，再用重心法选出各个区域的配送中心. 以式(3)计算运输成本，实验结果如表4所示.
外面喊，“怎么了静秋？”静秋回一声“没事”，跑出去，拿了抹布，想将桌子擦干净，却越擦越花。楚墨弯腰将碎片拣起——不错的一把紫砂壶，片刻之间，再也回不去了。
表4 配送中心选址效果对比
从配送中心选址效果对比中可以看出，对比传统的K-means聚类算法，改进后的算法每天可节约运输成本41 626 420元，配送时间可节省30.62小时. 在保证配送时效的同时也节约了运输成本，在一定程度上说明改进后的算法较传统的算法在实际中能创造更好的经济效益.
6 结论
通过本文分析，得出以下结论：在进行区域划分时，考虑当地的实际情况越来越受到学者们的重视，本文以城市运输距离与居民消费能力的指标相结合，重新定义了对象间的相似性度量的距离因子，区域划分后后，利用重心法对聚类中心进行修正，科学的选出了每个区域内的配送中心，更符合西北地区发展的实际情况. 为了验证改进方法的有效性，本文对比传统Kmeans聚类对区域进行划分后选址结果，实验可知：较传统的K-means聚类，改进后的算法不仅节省了配送时间，而且大大降低了运输成本，具有很好的实际运用价值，为西北物流建设提供了参考. 但本文考虑的是西北地区这一宏观的区域性概念，还有地处偏僻的部落并未考虑在内，因此划分的区域未能全部覆盖，想要在西北地区完成物流的全部覆盖及其布局，还需要一些发展的契机，有待以后学者们去探索.
参考文献
1 Kanungo T，Mount DM，Netanyahu NS，et al. An efficient Kmeans clustering algorithm：Analysis and implementation.IEEE Transactions on Pattern Analysis and Machine Intelligence，2002，24(7)：881-892. [doi：10.1109/TPAMI.2002.1017616]
2 朱培芬，汉吉庆，杨华龙，等. 基于改进K-means算法的烟草配送区域划分. 物流工程与管理，2009，31(6)：84-85. [doi：10.3969/j.issn.1674-4993.2009.06.034]
3 谷炜，张群，胡睿. 基于改进K-means聚类的物流配送区域划分方法研究. 中国管理信息化，2010，13(24)：60-63. [doi：10.3969/j.issn.1673-0194.2010.24.028]
4 于晓寒，王东. 基于带约束K-means聚类的城市快递配送区域划分. 哈尔滨商业大学学报(自然科学版)，2016，32(5)：631-634.
5 沈默，戴冰洁. 物流配送中心重心法选址解析——以苏宁为例. 物流技术，2015，34(4)：192-194. [doi：10.3969/j.issn.1005-152X.2015.04.058]
6 叶浔宇. 基于聚类和重心法的区域配送中心选址应用研究.中国市场，2009，(23)：83-85. [doi：10.3969/j.issn.1005-6432.2009.23.026]
7 胡贤满，张燕，李珍萍. 带车辆路线安排的多配送中心选址问题的求解——基于SPSS和遗传算法. 物流技术，2010，29(1)：83-86. [doi：10.3969/j.issn.1005-152X.2010.01.028]
8 孔继利，顾苧，孙欣，等. 系统聚类和重心法在多节点配送中心选址中的研究. 物流技术，2010，29(3)：83-85.
9 杨善林，李永森，胡笑旋，等. K-means算法中的k值优化问题研究. 系统工程理论与实践，2006，26(2)：97-101. [doi：10.3321/j.issn：1000-6788.2006.02.013]
Multi-Regional Logistics Distribution Center Location Method Based on Improved K-means Algorithm
LU Ling-Lan，QIN Jiang-Tao
(Business School，University of Shanghai for Science and Technology，Shanghai 200093，China)
Abstract： Focusing on the issues that the number，location，and coverage of multi-regional logistics centers of distribution centers are unknown，an improved k-means clustering algorithm is proposed. Based on the urban economic gravity model，this algorithm combines the urban transportation distance with the indicators of household consumption capacity，redefines the distance factor of the similarity measure between objects. The idea of density is introduced into the k-means algorithm，and the concept of intra-class difference mean is raised to determine the optimal number of clusters. After the partition is implemented，the centroid method is used to determine the final distribution center in these areas. Finally，in case study，we analyze the location process of constructing logistics distribution centers in 37 cities in the western region，and compares them with the traditional k-means clustering results. The comparing result shows that the improved algorithm not only saves the delivery time，but also greatly reduces the transportation cost and has sound economic value.
Key words： multi-regional distribution center location; k-means algorithm; urban economic gravity model; center of gravity method; northwest logistics
引用格式： 鲁玲岚，秦江涛.基于改进的K-means聚类的多区域物流中心选址算法.计算机系统应用，2019，28(8)：251-255. http：//www.c-sa.org.cn/1003-3254/7029.html
① 收稿时间： 2019-01-29;修改时间： 2019-02-26;采用时间： 2019-03-14;csa在线出版时间： 2019-08-08
标签：多区域配送中心选址论文; K-means聚类论文; 城市经济引力模型论文; 重心法论文; 西北物流论文; 上海理工大学管理学院论文;