摘 要：在市场机制自发调节下发挥政府规制的影响力，整合京津冀地区物流配送资源，构建基于“一核”消费市场为导向、“双层+多中心”物流配送节点建设为重点、京津冀地区交通网络为依托的环首都流通圈鲜活农产品物流供应链系统。运用CLRIP集成优化方法对环首都流通圈物流供应链的配送中心选址、配送路径、库存控制进行协调优化，建立供应链物流节点空间布局优化模型。参考赢商大数据，以北京市六环内27个商圈为客户点，以北京市六环附近已建配送中心及拟建配送中心为备选配送中心点，验证优化模型的可行性和实用性，为环首都流通圈物流供应链的配送中心提供最优空间位置布局、配送路径和库存水平点，在成本可控下实现最优选址—路径—库存，形成多个“组团式”的供应链物流节点空间布局。
关键词：环首都流通圈;物流供应链;物流节点;选址—路径—库存问题;禁忌搜索算法
在有序推进北京区域性专业市场、区域性物流基地等非首都功能产业疏解的同时，为了促进物流产业转型升级更好地适应消费升级对物流配送服务的要求，京津冀三地商务部门规划到 2020 年基本建成布局合理、高效畅通、环境友好和协作共赢的环首都 1 小时鲜活农产品流通圈。随着疏解非首都功能的推进，北京市区内各批发市场与物流仓储设施逐步外迁，北京市区居民对农产品的需求在很大程度上依赖于高效的物流体系建设。国内外大量研究成果表明，城市物流发展定位确定和发展路径规划应建立在城市自身能力特征和城市体系结构的基础上，物流企业的空间布局及区位选择是影响物流业健康发展的重要因素，郊区以其低廉的土地价格、联系城市中心和远程物流的优势成为物流企业优先选择的区位[1];并从宏观环境、市场规模与交通基础设施等角度研究城市物流企业的区位选择问题，探索物流企业的空间分布及其演化特征[2-7]。
随着新算法的不断涌现，有学者开始关注于节点选址、库存控制与运输管理三个问题的集成优化研究。Liu等尝试将库存控制引入到节点选址—运输管理问题中，提出选址—路径—库存集成问题(Combined Location Routing and Inventory Problem， CLRIP)，以最小成本为目标建立数学模型，应用两阶段启发式算法进行求解，并用模拟的方法对模型及算法进行了测试和评价[8];崔广彬等在多时期随机库存决策的供应链二级分销网络中使用双层规划法，在给定的多个潜在物流供应链节点中选出较好物流供应链节点的位置并确定巡回运输路径以及每条路径上的最佳订货量[9];马汉武等采用了 JITD策略，从决策者和客户两个角度出发，用先聚类再分配的方法解决了制造商到潜在分销中心的距离对选址产生的影响，使用双层规划建立了CL⁃ RIP集成优化模型[10];陈德慧等在模糊随机需求变化条件下将模糊随机的成本模型利用模糊可能均值法转化为确定的混合0-1整数规划CLRIP模型，设计了在最优库存策略下再运用禁忌搜索算法改进初始解的两阶段混合启发算法[11];唐金环等针对产供销一体化企业，引入碳配额差值系数，构建了结合碳配额差值的多目标CLRIP模型，并使用多目标混合粒子群算法对问题进行求解[12]。
物流供应链节点空间布局涉及节点选址、库存控制和运输管理三个方面，在以往的文献中，节点选址、库存控制与运输管理常常被分别考虑，无法实现最优物流节点空间布局且研究对象通常具有确定性，对模糊随机需求下的多因素组合决策问题少有研究。为此，本文从物流供应链的视角研究环首都流通圈农产品物流供应链系统，构建具有模糊随机需求的鲜活农产品三级物流供应链节点空间布局优化问题，在成本可控下实现配送中心及时、可达地完成配送。
1 确定优化对象与目标
北京最新城市规划(2016—2035 年)的出发点就是着力解决供给、需求、质量方面存在的问题，推动相关产业向品质化、便利化、精细化方向发展。为了适应消费结构升级、促进京津冀协同发展、推动物流产业转型升级、疏解非首都核心功能的发展目标，2017 年 5 月京津冀三地商务部门联合出台《环首都1小时鲜活农产品流通圈规划》[13]，到2020 年形成“一核双层、五通道、多中心”的环首都 1 小时鲜活农产品流通圈。结合政府规划所制环首都流通圈示意图(图1)。
其中，中间红色区域是“一核”，即北京市六环以内的特大型消费市场;灰色区域是“内层”，泛指从北京六环到河北天津的临近地区;灰色区域向外延伸至灰色环线之间的黄色、绿色区域是“外层”，地理上指天津、河北这两个环京区;“五通道”是以进京高速公路为主，国/省道及铁路为辅，打造西北、西南、东部、南部、东南五个方向的农产品进京通道;“多中心”则是依托“五通道”干线，建设集聚生产、中转、加工、仓储、交易、配送等功能的物流基地。
为落实规划目标，构建环首都流通圈鲜活农产品物流供应链，以由产地组织(“多中心”)、物流平台企业(“双层”)和含有宅配、自提功能的销售终端(“一核”)所组成的物流供应链系统为优化对象，采用CLRIP集成优化方法，参考赢商大数据选取北京市六环内 27 个商圈为代表，形成多个“组团式”的供应链物流节点空间布局，提高产地、物流中心、配送中心等节点衔接度，缩短物流中心与交通枢纽的流转距离，加快环首都农产品流通主体集约化、均衡化和规模化发展，推进农业产业、物流产业转型升级。对环首都流通圈鲜活农产品物流供应链系统的研究，可以充分利用规模效应和协同效应提升物流服务效率和效益，优化资源配置，提高农产品的配送效率，提升服务质量。
2 构建农产品物流供应链系统
物流供应链是以第三方物流企业为核心，通过优化配置流通过程中的资金流、物流、信息流、商流等，实现物流各功能环节、功能组织相互协同的网络组织，其中包含了提供产品或服务的整个过程中的所有物流活动[14]。目前有关物流供应链的研究主要集中在以成本、效率为考量的物流供应链构建研究和供应链协调及风险预警等方面[15-17]。
鲜活农产品物流供应链，是以生鲜农产品为研究对象，围绕链上第三方物流企业，将产品从生产、预冷、初级加工包装、冷藏贮存、冷藏运输、流通加工、末端配送、分销过程中的物流、信息流、商流进行调控，协调鲜活农产品生产者、农资供应商、流通加工企业及消费者之间的利益关系，最终实现鲜活农产品在低温状态下高效、安全、低损耗的流通，使供应链整体效益最佳。而物流供应链的效率则取决于各节点的有效衔接与合理布局，其主要节点情况：上游节点为各类产品产地及加工中心;中游节点为中转点，包括各类冷藏库、产销批发市场及配送中心等;下游节点集聚销售终端群体，如超市、餐厅等[18]。本文结合环首都流通圈规划落实目标，基于物流供应链的内涵、环首都流通圈规划布局及已调研第三方物流企业与物流基地的建设情况，构建由产地组织、物流平台企业和含有宅配、自提功能的销售终端所组成的鲜活农产品三级物流供应链系统(图2)。
①产地组织(“多中心”)。产地组织为沿西南、南部、东南、东北、西北五条进京通道的农业生产基地以及农产品交易集散中心。产地主要配置了鲜活农产品生产、初级加工、临时贮存、产品中转、运输配送、完成交易等的设施。
②物流平台企业(“双层”)。“外层”区域形成基于流通加工的大型冷链物流生产加工基地和物流基地，包括集散中心、仓储中心、流通加工中心等，主要完成鲜活农产品清洗、整理、分级、保鲜等初级加工，一般分布于农产品生产基地附近。“内层”建设流通、分拨、配送中心，集农产品的冷藏贮存、冷链运输、冷链配送、终端零售以及信息处理等环节于一体的整套现代物流服务体系，主要用于及时为北京郊区和来自津冀集散中心的鲜活农产品进入首都消费市场提供便捷高效优质的配送服务，以实现最后1小时配送。
③销售终端主要是“一核”内北京地区的各类销售网点，比如新发地、农贸市场、连锁超市、社区实体店、宅配或自提功能的新零售商店等。区域内重点加快传统农产品批发市场的疏解和转型升级，完善流通终端网点建设，改造、升级社区菜市场和连锁销售便民终端，发展产地直供社区终端，提升鲜活农产品品牌化、连锁化、规范化发展水平。
3 建立配送中心空间布局优化模型
3.1 模糊随机需求的数学描述
供应商根据末端客户的订货来决策订货量和配送中心的库存，由于鲜活农产品的生产周期和提前订货期问题，在这段时间里市场变幻莫测，客户需求可能高于预期订货量，或者相近，或者减少，这种情况下用模糊随机变量来描述最合适，只用随机或者模糊都不够准确。例如，下个周期的农产品市场需求增加，可达 5 000 单位的概率为 0.3，与预期目标需求量相近即4 000的概率是0.5，需求量接近减3 000单位的可能性为0.2。
3.2 问题描述与模型构建
3.2.1 问题描述
CLRIP模型是一种集成优化模型，可用来解决潜在配送中心的最优空间布局情况，为便于研究这里做一些抽象处理。由于终端实体店数量比较多且分散，不便于规划设计，所以可以将一条线路上的和地理位置相近的实体店及超市抽象为一个客户需求点，这样也可以保证每个点的日需求量。三级物流供应链系统抽象化结构(图3)，其中物流基地到配送中心为一级网络，配送中心到客户需求点为二级网络。客户点的需求是模糊随机的，在客户需求点的数量与位置、备选配送中心点的数量与位置确定的情形下，备选配送中心点会根据上游供应商的决策确定库存量。由于鲜活农产品的生产周期性，所以存在订货提前期，采用周期性检查库存与生产周期相匹配的策略保证库存量，提供配送服务时采用巡回方式。基于上述情况确定配送中心最优空间布局、最佳订货周期、库存保有量及最优配送路径。
3.2.2 基本假设
①环首都配送中心所配送的鲜活农产品为单一农产品。②各备选配送中心的位置已知，即目前北京六环附近所有潜在的配送中心位置点，并且各个备选配送中心在单位时间内的固定成本已知。 ③每个备选配送中心点年工作天数均为360天，且工作日内都要为其客户需求点提供一次配送服务。 ④每个客户需求点只配置一组车，每条配送线路仅有一组车为其提供服务，且各线路运输车辆数量车型均相同;任一条配送线路上的客户日需求量总和不能超过车辆的载重能力，且配送车辆从备选配送中心点巡回路线回到出发点。⑤各客户点日需求量模糊随机，各自独立。⑥采取( ) T, Sj 周期性检查库存策略。
3.2.3 模型参数和决策变量
①模型参数。I为所有客户需求点集合I = { }i ; J为所有备选配送中心点集合J = { }j ;K为所有冷藏运输车辆集合 K = { }k ;Fj 为各备选配送中心点 j 的 选址成本，即固定费用;C1为物流基地到各个备选配送中心点的一级网络的单位量运输成本;C2为各备选配送中心点到客户需求点的单位距离成本;C3 为每次的订货成本;T为订货时间周期;
3.3 CLRIP优化模型求解
本文先根据枚举法算出的最优库存量、订货周期，再采用启发算法计算模型的初始解，提高初始解的质量，最后利用禁忌搜索算法对初始解进一步改进直到搜索最优解。禁忌搜索算法流程如图4。
主要步骤如下：
步骤1：根据优化过程与约束条件得出最优库存水平与订货周期，并进行随机匹配客户和配送中心，车辆群与配送中心匹配，得到初始解;随机将客户分配给各配送中心和车辆，取得初始解;
步骤2：初始解的优化改进。运用禁忌搜索算法，改进选址分配，并在此优化基础上继续改进库存—路径。通过两阶段的持续改进，不断优化选址 —路径—库存决策，直到满足终止条件为止。
步骤3：按照情形1和情形2分别求出最优解，两相比较确定最终结果。
4 验证模型的可行性和实用性
以北京市六环附近已建配送中心及拟建配送中心为备选配送中心点，以北京市六环内 27 个商圈为客户点，依据各类参数为环首都流通圈农产品物流供应链的配送中心节点提供最优空间位置布局、配送路径和库存水平点，在成本可控下实现最优选址—路径—库存。为了便于计算确定最优配送中心的空间布局方案，这里随机选用河北省沧州市盐山县东北部的北京新发地盐山农副产品批发市场作为鲜活农产品的物流基地。
4.1 相关节点的选择
4.1.1 客户需求点的确定
本文研究环首都流通圈“一核”这个北京市六环内居民鲜活农产品消费需求，选择赢商大数据中心统计区域客流量大小的数据作为典型代表，通过综合考量、整理、归纳，北京市六环内商圈主要有27 个(图 5)，分别为：三里屯、CBD、王府井、西单、望京、燕莎、金融街、东直门、朝青、西红门、常营、西直门、崇文门、双井、朝外、太阳宫、中关村、亚奥、北苑立水桥、公主坟、万柳、东坝、鲁谷、通州、黄村、顺义、亦庄。为便于数据处理计算，根据商圈地理位置及常住居民数，综合聚类分析整理，将商圈划分成12个范围，拟定为12个客户需求点，并运用重心法算出的重心位置作为客户需求点坐标，12个客户需求点的坐标依次为：(77，84)，(75，122)，(90， 83)，(90，106)，(55，39)，(83，49)，(62，91)，(120， 90)，(85，92)，(138，160)，(36，85)，(49，117)。北京市客户需求点聚类示意图，如图6。
4.1.2 备选配送中心点的确定在确定
备选配送中心点时，将已建成的配送中心或者正在改建的为保障城市农产品供应的冷链物流配送中心直接作为备选点，例如：通州马驹桥物流基地、黑庄户农产品物流配送中心、东郊农产品配送中心、南口农产品物流配送中心、北京超市发物流配送中心、大红门货场等6个。北京市规划拟建的配送中心及流通加工中心分布在昌平、朝阳、海淀、通州和丰台区，在规划拟建配送中心的区域，则根据土地价格以及政府规划(在五环路、六环路周边建立功能完善的物流中心)[19]找出比较符合条件的位置作为潜在配送中心点5个，虽然位置模糊，但是对最后结果影响不大，所以这里做假定处理。所有的备选配送中心点位置及选址成本见表1。
综上，北京市客户需求点与备选配送中心点位置分布示意图如图7。
4.2 原始数据采集及分析
选定的环首都鲜活农产品三级物流供应链系统中，备选配送中心点不仅包括已建或者正在改建的6个配送中心点，还有5个规划拟建的比较合适的潜在配送点，其相对位置见表1。12个客户点的坐标位置和日均需求量见表2。已知下列参数：C1=3，C2=4，C3=800，L=10天，Cq =600，Ch=1，Cs =4，η=360 天。其中 1、2、3、4、5、6 分别代表通州马驹桥物流基地、黑庄户农产品物流配送中心、东郊农产品配送中心、南口农产品物流配送中心、北京超市发物流配送中心、大红门货场，是目前北京市在建和已经运行的配送中心，其余j点代表的是潜在配送中心，i表示客户点，d͂ (if f=1，2，3)的可能性分别为p1=0.3，p2=0.5，p3=0.2。
4.3 分析数据确定最优配送方案
运用Matlab 2014a 开发了4.1 所设计的方法的算法程序。禁忌搜索算法是一种随机搜索算法，它的初始解关联它的最优结果，因此选择8种不同的初始解数据，并记录相应的最优解，然后使用所设计的基于禁忌搜索算法的两阶段启发式算法，优化初始解(表3)。
由表3可知，使用所设计的基于禁忌搜索算法的两阶段启发算法所求得的最优解总成本数都比较接近，8 组数据中最大值是 225 990 000 元，最小值为225 740 000元，它们彼此之间的偏差很小，仅为0.1%，且总成本的最优解与初始解之间偏差范围为0.07%~0.21%，证明设计的CLRIP集成优化模型与基于禁忌搜索算大的两阶段启发算法是可行有效的。通过这种方法，计算得出的配送中心最优空间布局及库存、路径方案结果是相对准确。因此，由 表 3 可 以 得 出 最 优 成 本 方 案 目 标 函 数 值 为 225 740 000元。在最优成本条件下，运用倒推法计算拟开放的备选配送中心点位置及其配送路径、选址成本、库存成本、运输成本(表4)。
在最优解的各单项成本中，库存成本所占比重最大，是农产品的生产周期与订货周期之间客户需求点的总需求量导致的库存成本高，是为了保证需求;选址成本次之，选址成本偏高是由于北京市的地域位置所决定的，是没有办法调和的土地价格矛盾，属于正常范畴;运输成本相对较低，是由于大部分的备选配送中心点在选址时都考略到流通时间与配送效率，设立于北京市六环附近，本身运输距离就短，加之进京交通的便利性，因此运输成本属于三项成本中所占比例最低无可厚非。参照最优解，从备选配送中心点中开放2、5、6、7、10、11 共 6 个配送中心，例如车辆从配送中心 2 出发，给客户需求点 8 配货后再到客户需求点3，然后回到配送中心，其余线路以此类推，集成绘制相应的配送中心最优选址—路径方案(图8)和配送中心最优空间布局(图9)。
下面对案例中的相关数据和结果做简单的分析与总结。由表1可以看出，各备选配送中心点选址成本相差较大，例如北京超市发物流配送中心是黑庄户农产品物流配送中心选址成本的 3.9倍，造成这种大落差的关键原因是地理位置不同，土地价格相差大。因此在规划物流节点布局时选址合理与否直接影响整个选址的成本高低，如果选择土地成本比较高的区域甚至对整个目标函数总成本的影响都比较显著。在本算例中，由于需求量与供应量的匹配度，11个备选配送中心点中只需要开放6 个配送中心就能满足需求，各个方案的目标函数成本，除了选址成本这一影响总成本的因素外，库存成本与运输成本也会随着方案的不同波动比较大。不同的方案中各个配送中心的配送客户需求点、配送路径均不同，点与点之间运输距离会随之改变，相应的运输成本就会相差很大;而各个配送中心匹配的客户需求点的总需求量波动不同，库存水平、缺货损失就会出现连锁反应的不同，也会带来目标函数值差异很大。因此，在规划配送中心的空间节点布局时，要使选址—路径—库存方案尽可能的达到最优化，除了考虑选址成本，各个配送中心的配送路径与库存水平同样影响重大，它们不仅控制末端运输成本，同时制约着配送的及时性与集成优化程度。
5 结论
本文根据北京市最新规划与环首都流通圈的规划布局，运用物流供应链、物流节点、环首都流通圈鲜活农产品三级物流供应链模型等相关理论，构建了以河北省盐山县的鲜活农产品物流基地为产地，以北京市六环内的27个商圈为客户需求点，确定环首都流通圈鲜活农产品配送中心的最优选址 —路径—库存方案，实现鲜活农产品供应链物流节点最优空间布局，为落实环首都流通圈规划，优化物流配送企业组织结构提供参考。
本研究的不足之处在于环首都流通圈的范围比较大，涉及到的供应链系统复杂，由于调研样本有限，考虑的因素可能不够全面。结合目前阶段性研究成果，指出进一步研究的方向：目标函数总成本除了选址成本、库存成本与运输成本外，在后续研究中可以将其他成本考虑在内，使成本最优方案更加准确;客户需求点是按照商圈进行规划，统计的客户需求相对比较模糊，后续研究可以考虑找到更好的统计方法;环首都流通圈鲜活农产品物流供应链的优化，除了优化配送中心空间节点布局外，还可以考虑从组织结构布局、供应链管理系统等方面进行优化。——论文作者：陈红丽，赵 爽，杨海波，张丽丽
参考文献：
[1] Markus Hesse，Jean-Paul Rodrigue. The transport geography of logistics and freight distribution[J]. Journal of Transport Geog‐ raphy，2004(12)：171 - 184.
[2] Junjie Hong，Anthony T. H. CHIN. Modeling the location choices of foreign investment in Chinese logistics industry[J]. China Economic Review，2007，18(4)：425 - 437.
[3] 韩增林，李晓娜. 第三方物流企业的区位影响因素研究[J]. 地域研究与开发，2007，26(2)：16 - 25.
[4] 曹卫东. 城市物流企业区位分布的空间格局及其演化——以苏州市为例[J]. 地理研究，2011，30(11)：1 997 - 2 007.
[5] 梁双波，曹有挥，吴威. 上海大都市区港口物流企业的空间格局演化[J]. 地理研究，2013，32(8)：1 448 - 1 456.
[6] 秦璐，高歌. 中国物流运营网络中的城市节点层级分析[J]. 经济地理，2017，37(5)：101 - 109.
[7] 蒋天颖，史亚男. 宁波市物流企业空间格局演化及影响因素[J]. 经济地理，2015，35(10)：130 - 138.
[8] Liu S C，Lee S B. A two-phase heuristic method for the multi - depot location routing problem taking inventory control deci‐ sions into consideration[J]. International Journal of Advanced Manufacturing Technology，2003，22(11-12)：941 - 950.
[9] 崔广彬，李一军. 基于双层规划的物流系统集成定位—运输路线安排—库存问题研究[J]. 系统工程理论与实践，2007 (6)：49 - 55.
[10] 马汉武，朱晖. 基于JITD的服务备件二级分销网络集成研究[J]. 工业工程，2011，14(6)：76 - 80.
[11] 陈德慧，陈东彦 . 模糊随机环境下 B2C 电子商务配送系统 CLRIP 集成优化模型[J]. 系统管理学报，2017，26(4)：744 - 753.
[12] 唐金环，戢守峰，蓝海燕，等. 考虑碳配额差值的选址—路径 —库存联合优化多目标模型与求解[J]. 管理工程学报， 2017，31(1)：162 - 168.
[13] 北京市商务委，天津市商务委，河北省商务厅环首都1小时鲜活农产品流通圈规划[R]. 2017.
[14] 梁双波，曹有挥，吴威. 长江三角洲地区物流供应链时空演化及其影响因素——基于国际货代企业数据的分析[J]. 地理研究，2017，36(11)：2 156 - 2 170.
[15] 郑哲文. 基于可靠技术的应急物流供应链构建[J]. 中国流通经济，2009，23(10)：54 - 56.
[16] 夏佐铎，谭亮 . 北京奥运物流供应链的构建[J]. 商业研究， 2008(2)：126 - 129.
[17] 李晓宇，张明玉. 农产品物流供应链风险生成机制及预警模式[J]. 管理现代化，2009(4)：47 - 49.
[18] Dey O，Chakraborty D. A fuzzy random periodic review sys - tem with variable lead-time and negative exponential crashing cost[J]. Applied Mathematical Modelling，2012，36(12)：6 312 - 6 322.
[19] 北京市商务委员会，北京市发展和改革委员会. 北京市“十三五”时期物流业发展规划[R]. 2016.
