不确定条件下公铁水多式联运多目标路径优化研究.pdf

第40 卷第4期2023年8 月文章编号：10 0 5-0 5 2 3(2 0 2 3)0 4-0 0 5 6-10不确定条件下公铁水多式联运多目标路径优化研究华东交通大学学报Journal of East China Jiaotong UniversityVol.40 No.4Aug.,2023杨洛郡，张诚，郭军华（华东交通大学交通运输工程学院，江西南昌3 3 0 0 13）摘要：多式联运运输的时效性和成本在现代物流的发展中是不可忽视的因素。针对公铁水多式联运的运输目标，主要研究了当运输时间、中转时间双重不确定因素服从随机分布时的绿色多式联运路径优化问题。构建以运输时间、碳排放、运输成本为目标函数，碳排放量为约束，建立运输时间、中转时间双重不确定条件下绿色多式联运路径多目标优化模型。并据此采用模糊自适应遗传算法（FAGA)和快速非支配排序遗传算法（NSGA-）设计多式联运路径优化策略；最后采用从南昌到柏林的路径数据仿真验证所提方法的有效性并对结果进行对比分析。研究发现基于NSGA-I算法的多目标优化结果较优，可以引导多式联运经营人调整运输方案，减少二氧化碳的排放量，为物流企业开展多式联运运输提供可供参考的依据。关键词：多式联运；路径优化；双重不确定；网络配置；碳排放中图分类号:U15本文引用格式：杨洛郡，张诚，郭军华.不确定条件下公铁水多式联运多目标路径优化研究 J.华东交通大学学报，2 0 2 3,40(4)：56-65.文献标志码：AResearch on Multi-Objective Path Optimization of Highway-Railway-Waterway Multimodal Transport Under Uncertain ConditionsYang Luojun,Zhang Cheng,Guo Junhua(School of Transportation Engineering,East China Jiaotong University,Nanchang 330013,China)Abstract:In the development of modern logistics,the timeliness and cost of multimodal transport are factors thatcannot be ignored.Aiming at the transportation goal of highway-railway-waterway multimodal transport,this pa-per mainly studies the path optimization problem of green multimodal transport when the dual uncertainties oftransportation time and transit time follow random distribution.The multi-objective optimization model of greenmultimodal transport path under the dual uncertainties of transportation time and transit time is established withtransportation time,carbon emissions and transportation cost as objective functions and carbon emissions as con-straints.Accordingly,fuzzy adaptive genetic algorithm and fast non-dominated sorting genetic algorithm(NSGA-I)are used to design the multimodal transport path optimization strategy.Finally,the path data from Nanchangto Berlin is used to verify the effectiveness of the proposed method and the results are compared and analyzed.This paper finds that the multi-objective optimization results based on NSGA-II algorithm are better,which canguide the multimodal operators to adjust the transportation plans and reduce the emission of carbon dioxide,pro-viding reference for logistics enterprises to carry out multimodal transportation.收稿日期：2 0 2 3-0 2-2 7基金项目：国家自然科学重点联合基金（U2034211）；国家重点研发计划（2 0 2 0 YFB1713700）；流程工业综合自动化国家重点实验室联合基金（2 0 2 2-KF-21-03）第4 期Key words:multimodal transport;path optimization;dual uncertainties;network configuration;carbon emissionsCitation format:YANG L J,ZHANG C,GUO J H.Research on multi-objective path optimization of highway-railway-waterway multimodal transport under uncertain conditionsJ.Journal of East China Jiaotong University,2023,40(4):56-65.多式联运作为一个更加绿色高效的运输组织方式，相较于单一的运输方式，不仅可以适应长距离、大运量，还能减排降碳、降低运输成本叫,符合可持续发展的理念2。在物流产业转型升级的关键时期，提升运输全过程服务质量和降低碳排量是重中之重13-4 。目前，国内外许多学者对多式联运进行了相关研究,并取得了一定的成果。Wang等 5 设计了一种基于多式联运网络变形求解的方法，其中以运输转运成本、运输转运延误成本和运输成本之和作为总成本，模型的优化目标是求得总成本的最小值并通过实验验证得出算法的有效性;Boussed-jra等 0 对多式联运网络运用双向研究策略来搭建模型，以运输高效和快速为优化目标，并通过遗传算法来寻求最优解。但是以上模型都没有考虑碳排放对路径优化的影响，环境污染和资源短缺带来的问题越来越严重，温室气体排放导致的全球变暖已经成为整个社会面临的严重问题。Chen等7研究了物流运作模式对碳减排的影响,较早的通过相关数学模型分析了低碳问题和物流问题；李顺勇等 8 在多路径时变网络下建立了低碳车辆的路径优化模型，针对交通状况为拥挤流状态设计了符合时变网络碳排量计算的方法；任腾等9 在车辆承重、客户时间窗和冷链产品保鲜约束下，同时在考虑客户满意度的前提下构建了以优化目标为最小排放量的车辆路径优化模型。以上研究虽然在物流运作等运输问题中考虑到了碳排放的影响，但将联运货物中转时间成本引入多式联运路径优化当中的研究尚处空白。陈维亚等 10 建立了多式联运路径优化模型并以最小的总成本作为优化目标，其中总共的成本费用包括运输成本、转运成本、碳税成本和质量损失补偿成本，通过设计遗传算法来寻求最优解;Jiang等 研究了具有二氧化碳减排目标及污染排放的多式联运网络问题。上述研究往往以温室气体CO2排放量作为碳排放量的目标函数，而碳排放往往还包含其他温室气体。不同的运输方式会有很大程度上不同的能源杨洛郡，等：不确定条件下公铁水多式联运多目标路径优化研究1问题描述与建模1.1问题描述随着碳达峰、碳中和目标的提出，交通运输业也必须响应国家号召朝着绿色、低碳的方向发展。本文选择运输成本、运输时间和碳排放量作为优化目标，以达到最优的运输策略。57使用量，同时会产生不同的碳排放量，运输的实际距离长短和车辆的载重量也会影响碳排放量。多式联运网络会受到种种不确定因素的影响，主要是因为其运输的复杂性、多环节和广泛的覆盖面。针对多式联运在不确定的条件下对路径进行优化的问题,Fazayeli等 12 建立了模糊数学模型分析考虑模糊需求情况下的路径优化问题，并通过遗传算法来寻求最优解。Yin等 13 优化了闭环多式联运网络的设计方案来解决新产品和退货产品的需求以及碳税的不确定条件问题。Fotuhi等 14 研究了网络拓扑结构不确定的多式联运路径优化问题。吴小芳 15将救援物资需求量、道路交通条件等设为混合不确定性因素,研究了物流的选址与路径规划在混合不确定条件下的问题，通过建立模型寻求最优解，较好地完成了物流选址与物资分配的问题。袁长伟等16 将鲁棒优化模型运用于考虑多式联运货物的运输成本和网络终端节点容量不确定性的问题。在多式联运路径优化问题上主要是以运输总成本为目标来寻求最优解，大多数考虑的是运输需求、运输时效等不确定因素，然而运输时间、转运时间等双重不确定影响往往被忽略掉。本文从多式联运实际运营角度出发，当运输时间、中转时间双重混合不确定因素服从随机分布时，以运输成本、时间和碳排放的最优函数值作为目标，碳排放为约束，构建运输时间、中转时间双重不确定条件下绿色多式联运路径优化模型，采用模糊自适应遗传算法（FAGA）和快速非支配排序遗传算法（NSGA-）比较多式联运路线场景的最优决策和算法的有效性，得到理想的运输方式和路线方案。581.27不确定条件下的多目标优化模型为了便于本文建立多目标优化模型并进行求解，现做出如下假设：1）假设运输线路和转运节点的运输能力不受限制；2）假设在运输过程中，货物量保持固定不变；3）假设只能选择一种方式在相邻两个节点之间运输；4）假设货物在转运节点的中转机会有且仅有一次；5）假设运输过程中不会遇见道路拥堵、自然灾害等情况。1.2.1符号说明N表示运输节点的集合；H表示除起点和终点外的运输节点的集合；M表示随机运输时间的情景集合，m为该集合的子集，mEM;J表示所有运输方式集合；q表示货运量；Tmimn表示完成货物运输的时间下限；Tmx表示完成货物运输的时间上限；S；表示货物到达i节点的时间；d表示从节点i到节点j时，如果选择第k种运输方式交通工具所需运行的距离；c表示从节点i到节点j时，如果选择第k种运输方式所需要承担的单位运输成本；c表示假如在节点i选择将运输方式由k换为l所需承担的运输成本；e表示选择第k种运输方式时的单位碳排放量；E表示运输过程中允许的碳排放量上限；T”，请分别表示节点i与节点i间采用第k种运输方式运输时间的随机变量与取值;T,t;T,t;T，t分别表示情境m中公路、铁路、水路运输时间的随机变量与取值；T,t;S,s分别表示不同情境中运输时间、转运时间的随机变量与取值;D,d；A,a分别表示节点i出发时间和到达时间的随机变量和取值；x表示节点i与节点j间采用第k种运输方式时取1,否则取0;y表示节点i由运输方式k转换为运输方式l时取1,否则取0。1.2.2目标函数1)运输成本min C=Ci+C2运输途中运输成本与货运量、运输距离、单位运输成本和运输方式有关。Ci=ZEEqdfefatieNjeNkeJ节点转换成本与货运量、不同运输方式之间的华东交通大学学报单位转换成本和运输方式有关。C=ZZ.ZqclylieHkEJleJ2）运输时间。运输时间分为在途运输时间和节点转换时间。min T=Ti+T2ieNjENkEJi=2k=1/=1式(5)为在途运输时间，由运输时间和运输方式组成,其中考虑运输时间的随机变量与取值；式(6)为节点转换时间，由运量、转运时间和转运方式组成，其中考虑转运时间的随机变量与取值。3）碳排放E=ZZZqdfertieNjeNkeJ1.2.3约束条件Zxi=1l VieNZy=1或0 VieNklxi-1+xi,i+1 2y;VieN,kEJ,leJEEmxS-S$+ZZx/y(ttf)VueNkelelV(i,j)eI,zEJ,mEMV(j,k)EI,zEJ,mEMTmZZtx+ZEytWTieNjeNkeJieNkeJleJVijENx/E(0,1)ye(0,1)(1)式(8)表示相邻两节点间进行运输工作时只能选择一种方式；式(9)表示运输过程中途经每座城市时都只有一次换装的机会；式(10)表示变量与变量之间的逻辑关系；式(11)表示货物运输的碳排放量必(2)须在被允许的限度之内；式(12)采用递推的方法表达出货物运输到各个地点（节点i)的时间；式(13)、2023年(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)(15)(16)(17)第4 期式(14)为运输和转运时间的连续性约束,式(13)是为了确保货物到达每个节点的时间等于货物从上一节点的出发时间与运输时间之和,其中在转换运输方式的时候考虑不同情境中运输时间、转运时间不确定的变量与取值，式（14)是为了确保计算货物在每个节点的出发时间时，可以通过到达时间、中转时间和等待时间相加得到；式（15)表示货物从起点至终点，总的运输时间和换装时间的和要在规定的货物到达终点的时间窗内；式（16)和式(17)是决策变量的取值。2算法设计2.1模糊自适应遗传算法设计2.1.1FAGA算法原理本文提出了一种FAGA算法，利用种群和个体适应度值的方差来衡量总体种群多样性和个体差异,其中基于Mamdani模糊推理系统改变P。和Pm以提高算法的种群多样性。FAGA算法基本原理如图1所示。measurement quantitiesGeneticControl parametersalgorithmControl parametersThe overall goal of the pathoptimization problem underuncertainty图1FAGA基本原理图Fig.1Basic principle diagram of FAGA模糊控制器从遗传算法中得到性能测量值及P。、Pm 的参数,通过模糊规则处理,将调整后的参数输出给优化算法，达到自适应在线调整算法的目的。首先根据函数分布确定随机变量运输时间和转运时间的取值，采用种群中的所有个体之间距离的方差来衡量种群的密集程度，并进行归一化作为模糊系统的输入，通过模糊语言转换成相对应的P。和Pm的模糊变量值，再反模糊化对算法中的P。和Pm进行修正。2.1.2算法步骤首先对前种群个体的3个目标求平均,求出其居中个体路径的成本、时间和碳排放3个目标值杨洛郡，等：不确定条件下公铁水多式联运多目标路径优化研究ni=11X=，J=nn，2=n式中：xi,yi,z;分别为每个个体的成本、时间、碳排放目标值；n为种群数量。随后求每个个体与中间个体的距离差，如下d;=V(x;-x)2+(yi-y)2+(zi-z)2(19)式中：d;为每个个体与中间值的距离。求解的所有个体与中间个体碳排放差的均方差V，如下2dn1V=n最后对输入量进行归一化处理，如下V-ViF=Vm-V式中：Vmx为群体中最大个体的距离的均方差值；Vmin为群体中最小个体的距离的均方差值。Performance为简单起见，模糊系统中隶属度函数取梯形隶属函数。均方差F、交叉概率P。和变异概率Pm的隶FuzzycontrollerCostTimeCarbonemissions59xi,yi,z如下ni=1属度函数如图2 所示。模糊规则如下。规则i：如果F是Qil,则P。是Wi,Pm是Wz;其中=1,2,n表示第i个模糊规则，Qi描述F状态的模糊语言变量，Wi描述P和Pm大小的模糊语言变量。本文对各模糊语言变量均划分为大、中、小3个级别。综合上述可以得到如下模糊规则：The membership function of F1.00.5000.10.20.30.40.5 0.60.70.80.91The membership function of PThe membership function of P.1.01.00.50.5F010L0.5 0.6 0.7 0.8 0.900.05 0.10 0.15 0.20图2 FAGA隶属度函数图Fig.2FAGA membership function diagram(18)(20)(21)601)如果F是小,则P。是小,Pm是大；2）如果F是中,则P是中,Pm是中；3）如果F是大，则P。是大，Pm是小。FAGA算法具体流程如下：Step1将最优路径分为不同个体，随机产生规模为N的初始种群，设置进化代数；Step2计算第n个个体运输的成本、时间和碳排放的适应度值,计算群体的初始均方差F；Step3将F模糊化,求出P。和P.在模糊语言中所对应的变量值；Step4在对P。和Pm进行反模糊化后运行,更新P和Pm的值，更新种群个体并判断此时优化算法是否满足3个多目标路径优化模型，否则回到Step2。2.2快速非支配排序遗传算法设计现实生活中，存在多目标优化问题，当优化目标个数增加到3个及以上时，目标函数之间经常是冲突的,其最优解并非单一解,而是一组非支配解17。2.2.1NSGA-II 算法原理NSGA-I算法是Deb等学者于2 0 0 0 年在NS-GA算法的基础之上提出来的 18 ,广泛应用于交通流预测、最短路径策略优化等领域 19 。图3为NSGA-算法流程图。首先通过Pareto优胜级别排序对初始种群中的个体包括路径和运输方式进行等级划分，级别越低则表示满足成本、时间、碳排放多目标路径优化模型的适应度越高。通过对成本、时间、碳排放3目标函数进行拥挤距离排序，优先选择拥挤距离大的个体进人下一代，这有利于保持种群的多样性。经过上述操作产生的种群称为子代种群Q,将其与父代种群P,相融合构成R，继续进行拥挤距离排序和优胜关系排序，循环操作，直到满足停止条件。2.2.2优化步骤本文选择NSGA-即带精英策略的非支配排序遗传算法对多式联运最优路径进行选择。Step1对模型的运输方式决策变量x和y进行二进制编码,区间为 0,1。将最优路径分为不同的个体，随机产生规模为N的初始种群，设置进化代数；Step2随机给出初始化种群Po，并对种群Po进行非支配排序，初始化每个个体的Rank值；Step3通过二进制锦标赛法从P,和Q.产生出组合种群R,=PUQ;华东交通大学学报Build an optimal alternative set ofreal-world problemsGenerate chromosome sequences togenerate initial parent populationsbased on path data and nodesGenerate initial paths and transportmode offspring populationsMerging parent and offspringpopulations.MultimodalMultimodaltransport costtransport timefitness functionfitness functionFast non-dominated sorting,congestion distance calculationChromosome selection,crossoverand mutation calculationNoParent-child mergerPath cost,time and low carbonfitness funtionIs theevolutionary countcompleted?YesOutput the optimal path that meets the threeobjectives of cost time and low carbon图3NSGA-I流程图Fig.3Flow chart of NSGA-IIStep4对R,进行非支配排序，并通过挤排和精英保留策略选出满足时间、成本、碳排放量3目标的N个个体,组成新一代种群Pi;Step5跳转 Step4,并循环,直至满足结束条件。2.2.3随机时间条件设计随机时间条件采用矩形概率分布如图4 所示。从图中看出，这里随机概率分布采用的是矩形均匀分布，在同样间隔的分布概率是相同可能性。通过参数和b来定义均匀分布,在坐标上分别是最小值和最大值，一般简写为U(,b),其概率密度函数为(x)=b-a,Highway transportRa i l w a y t r a n s p o r t eWaterway transportation杨洛郡，等：不确定条件下公铁水多式联运多目标路径优化研究a6Fig.5 Multimodal transport network from Nanchang to Berlin61勒12 个城市作为中间节点城市，每个节点之间需要承担相应的装卸任务。网络中两个节点之间的运输方式可以分为单一运输和多式联运,如图5所示。设定有10 个2 0 英尺的集装箱（配货毛重2 0 t)将从南昌运往柏林，假设相邻的两个城市节点之间的运输方式可从公路运输、铁路运输、水路运输中X任选1 3种。本文的参数设定情况如下：种群规模为10 0，货运量为2 0 t，最大迭代次数为4 0 0 次，从降低企业成本、节约能耗、提升运输效率等多个方面进行综合考量，最后得到考虑了多方面因素的多式联运路径优化方案。多式联运的中转时间和费用及能耗相关参数如表1所示，各运输方式的部分路段运输里程如表2 所示。本文采用Matlab对案例模型进行求解，采用FAGA模糊自适应遗传算法与NSGA-I算法的路径优化效果进行对比分析。本案例中FAGA算法设定初始交叉概率0.2,变异概率0.8。将3个目标函数成表1中转时间、费用及能耗Tab.1 Transit time,cost and energy consumptionTransit time/Transit cost/Transit mode(h/t)Railway-highway0.25Railway-waterway0.5Highway-waterway0.35YanjiUlaanbaatarangchunManzhouliErenhotXianAlashankouKashgarChengduShanghaiGuiya11HuaihuaHainanKaohsiung图5南昌至柏林多式联运网络Transit energyconsumption(yuan/t)(kg/t)97.5137.5115AstanaMoscowVolgogradlakatskMinskKipVinnitsaDuisburgKatowiceGwadarHamburgiotterdam.YangolBangkokKunmirNannin2.0125BerlinRomeisterTehranHo.Chi.MinBeihaiKualaHanoiManilaLumpurgaDcPhnom Penh喀什、瓜达尔、仰光、曼谷、胡志明、新加坡、鹿特62Tab.2Transport mileage of some sections of each modekmOriginatingArrivalnodenodeWuhanHuaihuaChongqingNanchangGuangzhouShanghaiNanjingChangshaZhengzhouXianWuhanChongqingShanghaiChongqingHuaihuaChangshaHangzhouNanningJinanNanjingZhengzhouWuhanQingdaoHarbinJinanZhengzhouWuhanShenyangQingdaoHarbinChangshaHuaihuaChongqingChengduXianGuiyangYanjiShenyangChangchun华东交通大学学报表2 各运输方式的部分路段运输里程本f碳排放f、时间F转化成一个单目标函数f来进行优化的，也就是说通过权值相乘得到f=wfi+waf+of transportationW:f3,选择W)=0.5、W2=0.3、W3=0.2。给定最大遗传代数TransitmileageRailwayHighway358334142385954901706670119910548875875826265365141 047740871866123082010941 0413563859279169228186406393566984655373933661 48615833804139768421 4741 18719851.732109410416026023023007287103394067096663052892023年Waterwaygen=400,FAGA算法采用模糊参考系统作为调节479两种算子的机构自适应选取交叉概率与变异概率，利用设定的各个参数量，最终得到将货物从南昌运抵柏林的最佳路径。图6 为FAGA算法优化路径，南昌为始发城市，德国柏林为终点城市，依次途经南京、郑州、西一安、怀化、贵阳、南宁、海南、台北、北海、昆明、广1286州、上海、满洲里、雅库茨克、二连浩特、阿拉山口、1 043一742丹、汉堡、华沙、伏尔加格勒、明斯克、基普、文尼察、卡托维兹共30 个城市。根据表3，得出从南昌至柏林多式联运运输方案。在转换运输方式为公路运输从怀化至贵阳，从北海至昆明，从阿拉山口至喀什，从仰光至曼谷，从曼谷至胡志明，最后公路运输由伏尔加格勒至明斯克。在转换运输方式为水路运输从胡志明至新加坡，从新加坡至鹿特丹，从华沙至伏尔加格勒，从明斯克至基普，从文尼察至卡托维兹，最后从卡托维兹至柏林。其他运输方式为铁路运输。图7 为FAGA算法的适应度曲线图，由图可得FAGA算法适应度函数在迭代200次基本收敛至110 2 0,可见FAGA算法收敛速度较快。ShenyangQingdaoHarbinJinanZhengNanjingWuhNanchangChangsha-Highway transportationHangzhouRailway transportationFuzhouWaterway transportationFixednodeOptimize pathway nodesFig.6 FAGA algorithm to optimize Nanchang-Berlin multimodal transport networkYanjiUlanbaatargchuManzhouliErenhotXianAlashankotKashgarChengduShanghaiGuiyangHuaihuaKunmingNanninBeihaiHainanManilaIKaohsiungLumpur图6 FAGA算法优化南昌至柏林多式联运网络AstanaMoscowVolgogradYakutskMinskKipVinnitsaDuisburgiKatowiceGwadarHamburgRotterdamYangonamascBangkokHo Chi MinhngapHanoiKualaBerlinRomeTehranNew DelhiPhnom Penh第4 期1041.221.201.181.161.121.10050100150.2002250300 350 400Iteration图7FAGA算法多目标优化适应度曲线图Fig.7Multi-objective optimization fitness curve of FAGA表3FAGA算法优化多式联运最优路径Tab.3 FAGA algorithm for optimizing the optimal pathof multimodal transportationCity1Transit modeNanchangRailwayNanjingRailwayZhengzhouRailwayXianRailwayHuaihuaHighwayGuiyangRailwayNanningRailwayHainanRailwayTaipeiRailwayBeihaiHighwayKunmingRailwayGuangzhouRailwayShanghaiRailwayManzhouliRailwayYakutskRailwayErenhotRailwayAlashankouHighwayKashgarRailwayGwadarRailwayYangonHighwayBangkokHighwayHo Chi MinhWaterwaySingaporeWaterwayRotterdamRailwayHamburgRailwayWarsawWaterwayVolgogradHighwayMinskWaterwayKipRailwayVinnitsaWaterwayKatowiceWaterway杨洛郡，等：不确定条件下公铁水多式联运多目标路径优化研究algorithmCity2NanjingZhengzhouXianHuaihuaGuiyangNanningHainanTaipeiBeihaiKunmingGuangzhouShanghaiManzhouliYakutskErenhotAlashankouKashgarGwadarYangonBangkokHo Chi MinhSingaporeRotterdamHamburgWarsawVolgogradMinskKipVinnitsaKatowiceBerlin63表4 NSGA-II 算法优化多式联运最优路径Tab.4 NSGA-I algorithm for optimizing the optimalpath of multimodal transportationCitylNanchangWuhanZhengzhouChongqingHuaihuaNanningHainanTaipeiBeihaiKunmingGuangzhouShanghaiManzhouliYakutskErenhotAlashankouKashgarGwadarYangonBangkokSingaporeRotterdamHamburgWarsawVolgogradMinskKipVinnitsaKatowice图8 为NSGA-I算法优化路径，南昌为始发城市，德国柏林为终点城市，依次途经武汉、郑州、重庆、怀化、南宁、海南、台北、北海、昆明、广州、上海、满洲里、雅库茨克、二连浩特、阿拉山口、喀什、瓜达尔、仰光、曼谷、新加坡、鹿特丹、汉堡、华沙、伏尔加格勒、明斯克、基普、文尼察、卡托维兹共2 8 个城市。根据表4，得出从南昌至柏林多式联运运输采用NSGA-I算法的方案中，转换运输方式为公路运输从南宁至海南，从北海至昆明，广州至上海，从阿拉山口至喀什，从仰光至曼谷，最后公路运输由伏尔加格勒至明斯克。在转换运输方式为水路运输从新加坡至鹿特丹，从华沙至伏尔加格勒，从明斯克至基普，从文尼察至卡托维兹，最后从卡托维兹至柏林。其他运输方式为铁路运输。对上述结果进行分析可以得出如下结论：通过对运输成本的评估，多式联运的经营者更倾向于选用成本较低的铁路运输，在降低碳排放量的同时，也会增加一些运输时间。由于从南昌到柏林的水运并不便利，耗时长、成本高、连续性差，因此水运在本案例中竞争优势较低，故除了一些节点只存在水路运输之外，一般不会选择水路运输方式。NSGA-算法的Pareto最优解集仿真结果如Transit modeRailwayRailwayRailwayRailwayRailwayHighwayRailwayRailwayHighwayRailwayHighwayRailwayRailwayRailwayRailwayHighwayRailwayRailwayHighwayRailwayWaterwayRailwayRailwayWaterwayHighwayWaterwayRailwayWaterwayWaterwayCity2WuhanZhengzhouChongqingHuaihuaNanningHainanTaipeiBeihaiKunmingGuangzhouShanghaiManzhouliYakutskErenhotAlashankouKashgarGwadarYangonBangkokSingaporeRotterdamHamburgWarsawVolgogradMinskKip.VinnitsaKatowiceBerlin64华东交通大学学报2023年AstanaShenyangQingdaoHarbinJinanZhengzhouNanjingNuhanNanchangChangsha-Highway transportationRailiway transporaion HanghouWaterwaytransportationFixednodeOptimize pathway nodesFig.8 NSGA-I algorithm optimization of multimodal transportation network from Nanchang to Berlin2001951901851805.65.5x1045.45.3J25.21.25图9NSGA-算法送代次数和Pareto解集个数关系图Fig.9The relationship between the number of iterations ofNSGA-I algorithm and the number of Pareto solution sets图9 所示，可以看出算法的解集分布靠近坐标轴，需要优化的目标函数值小，解的质量更高。图10 为NSGA-算法在迭代过程中Pareto最优解集数和迭代次数的关系图，可以看出NSGA-算法的解个数在次数较少且较快时间内达到种群数量，故具有较好的收敛速度。由表5可见,NSGA-算法的最优解在运输成本、距离和时间3个维度上的目标值都优于FAGA。FAGA3个目标函数转化成一个单目标函数来进行优化的,这样对w1,W2,W3大小的选择必然会存在主观性，YanjiChangchunManzhouliErenhotLxianAlashankouKashgarChengduShanghaGuiyangGuangzhouHuaihuaangonKunmingNanningFuzhouHainanKaohsiung图8 NSGA-II 算法优化南昌至柏林多式联运网络。88：1.2551.26104J1UlaanbaatarGwadarHo Chi MinkBeihaiManilaLumpur4035302520151051.26500图10NSGA-算法选代次数和Pareto解集个数关系图Fig.10The relationship between the number of iterations ofNSGA-II algorithm and the number of Pareto solution sets表5NSGA-II 算法优化三目标结果Tab.5Three-objective optimizing results of NSGA-IITransitUnitcarbonTargetcost/(10*yua