低碳条件下鲜活农产品冷链物流配送路径优化
王 伟
（广东农工商职业技术学院管理学院，广东广州510507）
摘要： 近年来，低碳物流日益受到重视，低碳条件下鲜活农产品冷链物流配送路径优化分析也成为研究热点。通过综合分析配送车辆的运输成本、固定成本、制冷等相关成本，和在配送过程中的货损成本、碳排放成本，以及未能达到客户要求的服务时间窗而出现的惩罚成本作为具体的目标函数，建立低碳条件下鲜活农产品配送路径的优化模型，并为该问题的解决提出优化的改进建议，即采用2-opt局部搜索机制的蚁群算法进行改进。同时应用实例对模型和算法的有效性进行分析，包括对算法参数的敏感性分析。仿真试验及算法对比结果表明模型和结果都是有效的，能够为最终企业进行配送决策提供有力支持。
关键词： 低碳；鲜活农产品；冷链物流；配送路径优化
当前，温室气体减排日益受到各国的重视。在全球的所有碳排放系统中，交通运输大约能够占到15%，而道路碳排放占整个运输部门的碳排放的71%。冷链物流属于物流行业中排放量较大，同时也十分消耗能源的一种形式，当务之急是要分析如何在冷链物流中进行能源节约和减少碳排放等突出问题，最终促进经济发展和环境保护双赢，而这也是当前冷链物流学术界研究的热门话题。2017年，政府和相关部门就已经建议要在全国范围内启动碳排放市场，到2020年要在全国实行碳排放交易体系，故而提出和建立一个真正能长期发挥低碳和环保主导作用的碳税标准是当务之急。
1 研究现状与问题假设
1.1 研究现状
关于车辆路径问题（vehicleroutingproblem，VRP）的研究有50余年，到目前已取得了不少研究成果。研究主要集中在通过引入服务时间窗的概念，基于服务时间窗约束来对惩罚成本进行分析，构建带软时间窗的VRP模型；同时，在VRP求解算法上也取得了较大的成果，如在后期的研究中有学者设计了求解多配送中心车辆路径问题的基本算法，研究包含时间窗、多车场因素等的车辆路径问题，并在研究中给出求解该问题的量子粒子群算法［1］。部分学者提出，以往的VRP模型不应只是对经济效益进行分析，同时也要对环境和社会效益进行分析，进而提出减少碳排放的基本方法［2］。
笔者主要研究VRP模型在我国的冷链配送领域中的应用，在该领域中学界已经取得了一定成果。邵举平等考虑配送总成本和顾客满意度等相关指标，建立了鲜活农产品配送路径的多目标优化模型，而且通过改进的遗传算法对上述问题进行求解［3］。孙明明等对配送过程中的时间、温度等相关指标进行了分析，并构建了配送总成本最低化的基本模型，同时采用了节约成本法对相关问题进行求解［4］。
结合所需要研究范围的差异性，现有文献划分为3个部分：剔除低碳排放因素的一般物流VRP研究、未剔除低碳排放因素的普通物流VRP研究和不考虑低碳排放因素的冷链物流VRP研究。目前在学术界，关于未剔除碳排放因素的冷链物流车辆路径问题取得的成果还不多。另外，制冷设备能耗直接对碳排放的影响可能会被忽视。同时，较少有学者基于碳税来分析如何对配送路径进行优化。通过对上述问题进行分析，笔者主要研究基于碳税的带时间窗的鲜活农产品冷链物流配送路径如何进行优化等突出问题。由于能源消耗的变化很可能会对碳排放产生直接的影响，本研究考虑冷藏车门打开所导致运输过程和卸载过程中的不同能源消耗。而后构建低碳条件下低成本绿色的冷链物流配送路线模型。模型以最低配送总成本作为目标函数，具体考虑配送总成本中的运输、货损、制冷、惩罚、碳排放、固定成本等六大成本构成因素［3］。在针对碳排放成本进行计算的过程中，综合考虑行驶距离和车辆载质量等相关指标对碳排放量所产生的影响，采用改进蚁群算法对相关的模型进行求解，并在后期寻求最佳路径，以便鲜活农产品能够进行低碳配送，减少能源损失，将物流的持续发展和节能减排工作有机结合。
1.2 问题假设
为了能够在后期对所需要分析的问题进行有效界定，提出下列假设：假设1，只建立1个配送中心，而且配送的产品是鲜活农产品。假设2，因为采用的定位技术变得更加先进，每个客户所处的位置更加精确，所以研究中建设的物流配送中心和所有客户的位置和需求量都是已知的。假设3，在实际日常生活中，企业为了能够有效节省配送资源，一般只会安排某个车辆1次应对同一个客户的基本需求，所以研究中设定1个客户只能由1个配送车辆为其提供服务，而且客户的需求量能够得到满足［4］。假设4，因为本试验研究的是鲜活农产品冷链物流配送企业如何进行配送的问题，通常情况下不会要求回收和取货，因此本研究假设每辆车要从配送中心出发，完成某些任务之后会主动回到配送中心，配送过程中车辆只送货，而不负责取货或收货。假设5，因为鲜活农产品一般存在较短的保鲜保质期限，如果未能按照要求在规定的时间里送达，则必然会导致产品的质量下降，客户在收到货之后很可能会提出补偿要求。所以，研究中假设货物如果未能在给定的时间里送给客户，将会产生一定的处罚费用。假设6，在日常生活中，冷链物流配送企业将会提前对所需要配送的货物数量进行统计，一般情况下很可能会提供同种型号的冷藏车辆，所以研究中笔者假设配送中心存在着足够多的车辆来完成所有的配送任务，而且这些车辆的型号都是相同的，车辆的最大装载量是已知的，客户要求的配送量只能小于或等于配送车辆的最大载质量Qk。
在森林抚育过程中，林区内部树木的形态、规格等都会对整个林区产生一定的影响。通过应用森林抚育技术后，能够将整个林区内的树木进行管理，从而将林木之间存在的问题进行有效解决，将树木生长指标得到提升，真正的使优质树木的比例得到提升。
通过上面提出的假设，本研究的低碳条件下鲜活农产品冷链物流配送路径优化的问题可以表述为：某鲜活农产品配送中心要求向特定的客户群体提供配送服务，以冷藏货车作为配送工具，在了解客户所处位置，且每个客户都具有各自配送量、配送时限的要求，车辆完成配送任务之后就返回到配送中心的情况下，合理安排行车路线保证目标函数达到最优化［5］。
2 数学应用模型
2.1 数学符号意义
N指的是物流配送中心服务相应的客户数量；K指的是配送中心拥有的冷藏货车具体数量；dij指的是车辆从客户i位置行驶至客户j时跨越的距离；C1指的是车辆固定应用成本；C2指的是车辆的单位运输成本；qi指的是客户i相应的需求量；p指的是鲜活农产品具体的单位价格参数；Qk指的是配送车辆支持承载的最大质量；ρ（Qij）指的是车辆自客户i位置起驶向j客户时，运送Qij货物单位距离相应的燃油消耗参数；ε1指的是车辆在时刻Ei前送至客户i相应的单位时间惩罚成本；ε2指的是车辆在时刻Li以后送至客户i相应的单位时间惩罚成本；ti指的是车辆k面向客户i提供服务所花费的时间，其中 t0＝0；［Ei，Li］指的是客户 i提出的配送要求的服务时间窗。
关于农村涉土信访问题的探索及分析（蒋可挺） ............................................................................................6-48
在以上函数中，蚂蚁个体须在一定约束条件内包括时间与质量等。因距离正好是位于分母上，因此距离参数取值越小，相应的启发式因子取值会越大，而蚂蚁选择该节点的可能性也会越高。质量处于分子上，因此，下一客户节点货物需求量增高，启发式因子取值结果也可能会扩大化，选择该节点的可能性也会较高。因此，针对需求量高的客户点，优先安排配送，降低配送车辆质量与相应的燃油消耗。
2.2 建立模型
2.2.1 基于车辆固定成本的相关分析 通常情况下，该成本相对稳定，可以借助于常数进行分析。车辆固定成本除了固定损耗之外，还有驾驶人员的薪水和车辆其他成本等。该参数取值和行车距离、服务顾客数量。并不存在关系［6］，将其定义为Z1，则函数可以表述为：
2.2.2 基于车辆运输成本的相关分析 这类成本主要表现在油耗上，其与车辆运行里程之间存在着正比关系。将该成本设定为Z2，则函数可以表述为：
2.2.3 基于货损成本的相关分析 冷链VRP运输与一般VRP运输存在着明显差异。冷链VRP以易腐货物运输为主。一般VRP模型只须要对货物装卸等过程中出现的碰撞跌落等损坏问题进行考虑。但本研究分析的是鲜活农产品损坏成本问题，由此须考虑在运输、装卸环节因温度变化可能引起的质量问题。冷藏货物通常都表现出易腐性，对存储环境有着较高的要求，如温度、氧气浓度、湿度、产品含水率等变化，都会引起产品质量的改变［7］。特别是随时间与温度增加，这类货物会逐渐腐烂，质量下降或已不具备应用价值。如产品质量下降至一定水平，相应便会带来损失成本。冷藏车剩余货物具体的数量与客户现实需求存在着明显的相关性，且影响着冷链配送环节的货物损失。
在本研究分析中，提供了冷藏货物质量变量函数，具体表述为 D（t）＝D0e-∂t。这里的 D（t）意义是在 t时刻上，产品相应的品质；t代表的是产品运输时间，D0代表冷藏货物从配送中心启程时的货物质量，∂代表的是产品腐败率，∂参数取值与鲜活货物自身特性、运输车辆内部温度相关［8］。在本研究分析中，设定在运输环节货物所处环境温度参数是固定的，即温度为恒温状态，此时鲜活货物腐败率∂可以表现为常数，其随时间增加，鲜活货物质量表现出指数变化。对相同产品而言，在其他条件不出现改变的基础上，如存储温度增高，相应货物腐败率也会明显提高。
由此在车辆配送过程中，由起点至顾客i位置行驶时，运输时不开启车门，相应货物成本Z31函数可以表述为：
式中：P指的是鲜活农产品单价；qi指的是顾客i相应的需求量；tki指的是车辆k运输至顾客i位置所花费的时间；tk0指的是车辆k自配送中心出发的时间点。在特定温度条件下，产品相应的腐败率通过∂1来表示，yik取值表现为0-1变量，在车辆k为顾客i提供服务时，则yik取值为1，否则该参数取值为0。
在开启冷藏车车门时，受空气对流等因素影响，车箱内部温度会有所提升，此时鲜活农产品腐败率也会表现出相应改变［9］。设定此时腐败率参数为 ∂2（∂2＞∂1），则冷藏车辆运输至顾客i时，在卸货过程中开启车门，相应的货损成本表现为Z32，则存在：
散步的时间推荐是阳光最舒适时或餐后1小时以后。阳光最舒适时首先空气相对会好一点，其次是可以有阳光浴，促进维生素D的合成，帮助钙的吸收。餐后不应立刻散步，这样会导致本应帮助消化吸收的血液流向四肢，长期下来会影响消化功能，所以建议在餐后1小时以后散步。
式中：Qin具体描述的是冷藏货车离开顾客i时，货车内剩余产品质量；tsi具体描述的是顾客i所需服务时间长度。由此，在整个配送过程中，货损成本Z3满足如下函数：
2.2.4 基于制冷成本的相关分析 冷链VRP模型须考虑制冷所消耗的成本问题。该部分费用包含了运输过程中车辆维持低温所消耗的能量成本，卸载环节制冷系统产生的额外能源成本。将运输环节其成本设定为Z41，则存在：
将卸货环节其成本设定为Z42，相应函数表述为：
由此，则制冷成本函数为：
2.2.5 基于碳排放成本的相关分析 该部分成本具体指的是在配送环节，车辆运行会消耗一定的燃料并带来二氧化碳（CO2）排放相应的成本。求解碳排放量，是通过燃料消耗量与CO2排放系数的乘积来获得的。车辆燃油消耗程度不仅与运输里程存在关系，还与车辆运输质量相关［10］。相关学者通过统计分析获得了单位距离燃料消耗量ρ值，可以将其表述为依赖车辆载货量X相关的线性函数。这里将车辆本身的质量设定为 Q0，载货量设定为 X，则 ρ（X）函数可以表述为：
2.2.6 基于惩罚成本的相关分析 针对鲜活货物的配送，客户对配送时间有着较高的要求。如不能在预期时间内送达，则会对客户满意度评价带来负面影响，即引起惩罚成本Z5。
将运输车辆载货量最大值设定为Q，满载条件与空载条件下的单位距离燃料消耗量分别设定为ρ*与ρ0，则存在：
巴乌是音乐教学中常用的乐器，它音域较窄，简单易学，在教学过程中，老师先让学生熟悉巴乌的十个音节，知道是低音，1 2 3 4 5 6是中音，并让学生在巴乌上找出各个音对应的位置，然后音乐老师再进行示范演奏，让学生模仿自己吹奏单音，从低音到中音进行逐一体验式练习，培养学生的音高概念和识谱能力，正所谓“冰冻三尺非一日之寒”，再有天赋的学生也要通过反复操练才能把握好音调的高低和演奏的节奏。小学生乐器的学习不能过于死板，老师要以组织游戏的方式进行乐器的教学，寓教于乐是促进学生学习积极性的重要手段，同时也让他们感受到学习巴乌演奏的乐趣。
3.1.1 启发式因子应用设计 启发式因子属于启发蚂蚁节点选择的主要影响因素，也是蚁群算法操作的核心内容，该部分的设计效果直接关系着算法求解水平。在上述所构建的模型中考虑了车辆质量方面的影响，如在设计中只选择距离来作为启发式因子，则会对分析结果带来不利影响，引起燃油成本增加。由此，在启发式因子中引入货物需求量，让算法在下一节点选择时可以兼顾燃油消耗影响因素［13］。该模型启发式因子可通过如下函数来表示：
式中：w代指的是CO2排放系数值，Qij具体代指的是自客户i位置到达客户j位置时车辆相应的载荷量。ρ（Qij）描述的是车辆在（i，j）运行过程中，在Qij载质量条件下的单位距离燃油消耗量。在研究中，将碳税设定为C3，相应碳排放成本Z4可以表述为：
他坐起来，想着切身的事情。裹在脚上的毯子已经磨穿了，他的脚破得没有一处好肉。最后一条毯子已经用完了。枪和猎刀也不见了。帽子不知在什么地方丢了，帽圈里那小包火柴也一块丢了，不过，贴胸放在烟草袋里的那包用油纸包着的火柴还在，而且是干的。他瞧了一下表。时针指着十一点，表仍然在走。很清楚，他一直没有忘了上表。
具体各个物理量的含义如下：φij代指的是边（i，j）中的自启发量，存在
。τij-蚂蚁 k在边（i，j）上运行轨迹的长度信息素量；pkij-蚂蚁k由节点i移动至邻域节点j的概率；在初始时间点上，不同路径中的信息素浓度保持一致性，即存在τij（0）＝C，这里的C属于常数值。此时，顾客 i的蚂蚁 k在选择顾客j问题上所表现出的概率pkij，则通过以下函数来求解：
Stars and moon are dim in the evening, rain or snow falls from sky later on.
综上所述，多发性骨髓瘤的X线、CT及MRI临床表现具有一定特征性，通过检测可及时对患者做出诊断，且联合应用的诊断准确率显著较高，故值得临床推广。
函数（1）描述的是整体最低成本；函数（2）描述的是车辆运载能力约束信息；函数（3）描述的是确保每1个顾客仅被1辆车提供服务；函数（4）与函数（5）约束的是任意客户，仅1辆车出发1次或到达1次；函数（6）描述的是次回路消除相应条件；函数（7）描述的是配送作业连续性；函数（8）与函数（9）描述的是变量约束信息［11］。
（4）从该井动液面的变化可以看出，油井正常生产后，在日产液量基本不变的条件下，动液面从500m逐渐恢复到400m，说明该井渗透性能在逐渐变好。
3 算法应用设计
因VRP本质上为NP-hard问题，在求解分析中通常采取的是启发式算法［12］。而低碳条件下鲜活农产品配送路径方面的优化，实质上也属于NP-hard问题。但在复杂程度上，较之VRP明显更高，这也就对其算法提出了更高要求。蚁群算法是通过真实蚁群行为的分析所提出的一类算法，其在应用中表现出突出的鲁棒性、正反馈性，支持分布式计算与分析，且可以与其他方法结合到一起应用。在具体操作中，为避免求解结果为局部最优，在研究中借助改进蚁群算法来完成模型求解与分析。
由于在知识产权执法中自贸区的非法行为是低风险甚至是无风险的，所以自贸区便成为了滋生假冒产品走私与转运的温床。自贸试验区进口假冒产品后，对其原产地或专利与商标经过重新包装再次出口。这些违法行为不仅没有被规制，反而形成了愈演愈烈的趋势，对自贸区的经济秩序造成了混乱。
第二相电泳根据蛋白质分子量进行分离。NuPAGE 4%～12% Bis-Trid ZOOM 电泳胶 (Invitrogen)用于第二相电泳分离，电压125 V下持续电泳90 min。经过两相电泳分离的蛋白质使用考马斯亮蓝G-250进行染色6 h后，使用ImageMaster 2D Platinum进行电泳胶内染色蛋白定量比对。
3.1 蚁群算法应用的相关设置
由此，在配送中如自客户i位置起驶向客户j位置，在（i，j）区间内，车辆运行相应的碳排放量为：
式中：C3指的是每单位消耗CO2排放相应的环境成本参数；w指的是CO2排放系数值；Qij指的是自客户i位置到达客户j位置过程中车辆载荷量参数。
盾遂奔，未出晋境。……晋太史董狐书曰“赵盾弑其君”，以视于朝。盾曰：“弑者赵穿，我无罪。”太史曰：“子为正卿，而亡不出境，反不诛国乱，非子而谁？”孔子闻之，曰：“董狐，古之良史也，书法不隐。宣子，良大夫也，为法受恶。惜也，出疆乃免。”(《晋世家》)
3.1.2 移动概率应用规则 蚂蚁依据启发式信息与信息素，应用概率决策机制，从而完成下一个移动节点的选择。在获得构造优化解后，须对该解作具体的评估，依据评估结果来决定信息素更新时，信息素的释放量［14］。
综合以上信息，在考虑碳排放问题的基础上，进行鲜活农产品配送路径优化，相应模型为：
式中：α、β均属于支持调节的参数量，描述的是蚂蚁运动时，累积的信息素、自启发量在路径中作用权重。allowedk表示允许范围内蚂蚁k的随机值。
为了不出现停滞问题，在本研究分析中选择应用随机性选择与确定性选择相结合的方案，在搜索操作中对状态转移概率作动态调整。基于节点i上的蚂蚁k而言，其可通过以下函数来完成下一节点的选择：
式中：Jk（i）具体描述的是第k个蚂蚁在完成节点i访问后，还支持访问的节点集合；q参数为［0，1］范围内的随机数；q0在计算之前可以确定，通过该参数调节可以让算法在多样化搜索与集中搜索上实现平衡。此时，蚂蚁在选择下一节点进行转移时，就会产生一个在［0，1］区间内的随机参数，依据该参数的大小来确定蚂蚁转移方向。
由此，可以将ρ（X）表述为：
3.1.3 信息素更新方法 在计算过程中，为避免收敛操作过早，从而导致结果并非为全局最优解，要求所有路径上的信息素值都在规定的范围内，即［τmin，τmax］。低于 τmin信息素值则定义为 τmin，超过 τmax的信息素值则定义为 τmax，由此可以规避在某路径中其信息素远超其他路径信息素，引起全部蚂蚁快速靠拢到该路径中［15］。
（1）针对初始化信息素参数值 τij（t）＝C，并将其定位为最大值 τmax。
(4)由图6可以看出，结构的平动变形最大值和扭转变形最大值均出现在墩顶，结构整体平动变形最大值为1.253 4 cm，最大扭转变形值约为0.578′，平动变形较大，扭转变形很小；且由墩顶至墩中再到墩底，结构的变形逐渐变小，靠近墩底部，结构的变形值为0。
（2）每个蚂蚁在完成1次循环操作后，找到最短路径的蚂蚁才可将经过路径的信息素释放出来。存在：
（3）将参数 τij（t）放置在［τmin，τmax］范围内。如存在τij（t）＜τmin，则 τij（t）＝τmin；如存在 τij（t）＞τmax，则 τij（t）＝τmax。在进行最优解更新时，要求同时进行τmin与τmax参数更新。让
，这里的ρ代指的是信息素发挥发系数值，1-ρ描述的是信息素残留系数。σ参数值是通过试验来获取的，存在σ＝n／20，n代指的是节点数。3.1.4 局部优化改进 基于蚁群算法，应用添加2-opt局部优化算法，在全部蚂蚁完成最优解后，在所有路径信息素还没有更新前，对每代最优解作局部改进，以提升蚁群算法收敛速度。
3.2 算法计算应用步骤
第1步：设定参数 Nc＝0，τij（0）＝τmax，实行 τkij、Q、α、β、ρ参数初始化，在当前解集中应包含初始点，然后对其进行求解，定义为z。第2步：在其他的剩余点中，如并没有能够达到车辆质量与时间窗口要求的顶点j，则直接进入到下一步。反之，在满足车辆质量与时间窗口要求的剩余点中，以随机方式选定点j，求解转移概率Pkij参数，并将该参数与随机数（0-1）作对比分析，如满足要求，则将蚂蚁k转移至j点，将j点放置在当前解集中，如不满足要求，则需重新来选择顶点。第3步：在解集中存在所有点时，求解 zki，对蚂蚁个数记录为m←k，反之则存在k←k+1，此时会跳转至第2步重新操作。第4步：选择2-opt局部搜索机制，实现蚂蚁路径优化处理。第5步：对目标函数值进行求解，将当前非劣解作保留。第6步：对求解最优路径的各边（i，j）作分析，求解 τij（t+n）＝（1-ρ）τij（t）+Δτmijin。第7步：针对非最优路径各边（i，j）作分析，求解 τij（t+n）＝（1-ρ）τij（t）。第 8步：对所有边（i，j）作分析，设置Δτmijin←0，Nc←Nc+1。第9步：如 Nc参数值低于设定迭代次数，则进入至第2步操作。第10步：对当前最优解输出。改进蚁群算法操作流程见图1。
4 案例求解分析
4.1 结果及稳定性分析
4.1.1 案例试验1 对算法有效性作检验分析，在本研究中选择应用文献［16］提供的具体信息，服务客户设定为20个超市门店，车辆额定载质量参数为9000kg，柴油机型，空载等速燃料消耗量参数值具体为16.5mL／km。综合燃料消耗量参数值表现为233mL／km。模型相关参数信息见表1。
表1 模型相关参数对应数值
通过仿真分析得知，在顾客规模表现为蚂蚁数量15倍时，蚁群算法收敛速率与全局收敛性处于较好水平。因此，在本研究中将蚂蚁数量设定为15，迭代次数设定为Nc＝100。改进蚁群算法后，参数α＝1、β＝3，应用matlab7.10.0进行编程与模型求解分析，车辆运行路径具体为：车辆1具体路径为0→5→8→11→10→19→13→0；车辆2具体路径为0→17→13→15→16→7→2→0；车辆3具体路径为0→3→12→18→9→4→6→1→0。
基于以上分析，总共配置3辆车，须完成20个顾客服务。通过路径信息可以了解到每辆车服务顾客的编号，在配送服务完成后车辆驶回配送中心，整个配送过程成本总计为9576元。
为考察算法参数对模型结果的影响，在这里调整α、β参数取值，并作仿真探究，运行结果见表2。由此可以得知调整α、β取值，其对综合成本所带来的影响并不显著，结果比较稳定。
表2 试验1α、β参数调整下的运行结果
4.1.2 案例试验2 为更好分析算法有效性，在本研究中应用不同规模的案例作求解，应用文献［17］所提供的Solomon标准测试问题库中的案例作进一步验证。实例选用的是R208，在该实例中，顾客数量设定为n＝100，车辆载质量表现为ck＝1000。设定蚂蚁数量为65，迭代次数同样为Nc＝100。改进蚁群算法后，参数 α＝1、β＝3，用 matlab7.10.0（R2010a）进行编程与模型求解分析，获得优化路径如下：车辆1具体路径为0→58→53→40→21→73→75→56→54→39→80→68→3→77→50→33→0；车辆2具体路径为0→27→28→26→12→74→1→70→30→20→51→9→81→35→65→71→66→32→90→63→64→37→98→85→91→100→97→95→2→96→6→94→13→59→57→41→22→72→29→0；车辆3具体路径为0→17→84→5→60→83→72→52→31→88→7→82→48→47→36→49→19→10→34→78→24→25→55→4→74→42→43→14→16→64→93→89→45→0；车辆4具体路径为0→69→79→62→11→8→99→86→46→87→0；车辆5具体路径为0→92→44→38→15→23→67→0。
马铃薯样品：将采集马铃薯样品去壳，充分混匀，分取100 g样品两2，分别装入封口样品容器中，粘好标签，贮存于-20 ℃冰柜中保存。
在以上操作中将问题规模作了拓展，旨在考察分析算法在较多车辆路径中的分析是否具备应用价值，运行结果见表3。结果表明在大规模车辆路径中，该算法应用具有较为突出的稳定性。
表3 试验2α、β参数调整下的运行结果
4.2 算法应用对比
针对试验1所提出的改进蚁群算法，将其与基本蚁群算法、遗传算法等作进一步对比，多种算法对比分析结果见表4。在考虑了碳排放问题的基础上进行配送路径优化，改进算法相对基本蚁群算法而言，其在整体成本上下降了2.3%，相对遗传算法来看，在整体成本上下降了1.2%。
表4 不同算法的对比分析
5 结语
综上所述，随着物流业节能减排和其他各项活动的推进，低碳条件下的车辆路径问题是物流领域中十分热门的话题。本研究分析了碳排放的鲜活农产品配送路径的优化问题，综合分析了配送总成本中的运输、货损、制冷、惩罚、碳排放、固定成本等六大成本。在综合分析配送总成本的基础上，建立了该问题的优化模型，对基本蚁群算法的启发式因子、移动概率选择规则等相关因素和条件进行了优化，同时也在后期的算法中融入了2-opt局部优化方法，并提出了改进的蚁群算法求解方案，仿真试验及算法对比的结果显示算法是可行有效的，可以规避局部最优，故本研究模型和求解算法能够为在低碳环境下物流企业的配送工作提供有力支持。
参考文献：
［1］张 倩.全国统一碳交易市场即将全面开放［J］.环境经济，2017（7）：26-27.
［2］康 凯，韩 杰，普 玮，等.生鲜农产品冷链物流低碳配送路径优化研究［J］.计算机工程与应用，2018，54（6）：10-18.
［3］何有世，马腾飞.B2C环境下生鲜农产品物流配送路径优化研究［J］.商业经济研究，2017，34（5）：93-95.
［4］陶 云，张鹏程.基于改进遗传算法的冷链物流路径优化研究［J］.蚌埠学院学报，2016，5（3）：85-90.
［5］石建力，张 锦.需求点随机的分批配送VRP模型与算法研究［J］.控制与决策，2017，32（2）：213-222.
［6］凌海峰，谷俊辉.带软时间窗的多车场开放式车辆调度［J］.计算机工程与应用，2017，53（14）：232-239.
［7］潘茜茜，干宏程.考虑碳排放的冷链物流配送路径优化研究［J］.数学的实践与认识，2016，46（2）：62-68.
［8］叶 勇，张惠珍.多配送中心车辆路径问题的狼群算法［J］.计算机应用研究，2017，34（9）：2590-2593.
［9］颜 瑞，朱晓宁，张 群，等.考虑二维装箱约束的多车场带时间窗的车辆路径问题模型及算法研究［J］.中国管理科学，2017，25（7）：67-77.
［10］兰 辉，何琴飞，边 展，等.考虑道路通行状况的冷链物流配送路径优化［J］.大连海事大学学报，2015，41（4）：67-74.
［11］董冬艳.基于生鲜农产品的冷链物流配送路径优化研究［D］.沈阳：沈阳大学，2017.
［12］张如云，刘 清.考虑低碳的城市配送车辆路径优化模型研究［J］.工业工程与管理，2015，20（4）：29-34.
［13］邵举平，曹 倩，沈敏燕，等.生鲜农产品配送中带时窗的VRP模型与算法［J］.工业工程与管理，2015，20（1）：122-127，134.
［14］李亚男，刘联辉，李晓曼，等.低碳约束下城市冷链物流配送系统优化研究［J］.中国市场，2016（10）：36-37，39.
［15］孙明明，张辰彦，林国龙，等.生鲜农产品冷链物流配送问题及其路径优化［J］.江苏农业科学，2017，45（11）：282-285.
［16］Wang S，Tao F，Shi Y，et al.Optimization of vehicle routing problem with time windows for cold chain logistics based on carbon tax［J］.Sustainability，2017，9（5）：694.
［17］Cheng R，Gen M，Tozawa T.Vehicle routing problem with fuzzy due-time using genetic algorithms［J］.Journal of Japan Society for Fuzzy Theory and Systems，1995，7（5）：1050-1061.
中图分类号： F252
文献标志码： A
文章编号： 1002-1302（2019）13-0010-05
王 伟.低碳条件下鲜活农产品冷链物流配送路径优化［J］.江苏农业科学，2019，47（13）：10-15.
doi： 10.15889／j.issn.1002-1302.2019.13.003
收稿日期： 2018-04-15
基金项目： 广东省青年创新人才类项目（人文社科）（编号：2017GWQNCX038）；中国物流学会、中国物流与采购联合会研究课题（编号：2019CSLKT3-103）。
作者简介： 王 伟（1984—），男，湖北广水人，硕士，讲师，主要从事农产品流通研究。E-mail：weiwang＠gdaib.edu.cn。
标签：低碳论文; 鲜活农产品论文; 冷链物流论文; 配送路径优化论文; 广东农工商职业技术学院管理学院论文;