考虑碳税成本的综合运输网络货流分配优化.pdf

TRANSPORT RESEARCHVol.9 No.3收稿日期：2023-05-02基金项目：国家自然科学基金项目（71971021）；北京交通大学基本科研业务费科研学术活动项目（2019JBM334）第一作者：卢霞（1999），女，江西南昌人，博士研究生，研究方向为综合交通运输理论与技术。E-mail:通讯作者：毛保华（1963），男，湖南祁阳人，博士，教授/博导，研究方向为综合交通发展政策、运输组织理论与方法。E-mail:考虑碳税成本的综合运输网络货流分配优化卢霞1,2，黄俊生1,2，毛保华1,2，陈硕1,2（1.北京交通大学 中国综合交通研究中心，北京 100044；2.北京交通大学 综合交通运输大数据应用技术交通运输行业重点实验室，北京 100044）摘要：为推进煤炭等大宗货物运输“公转铁”，从运输参与方（托运人、承运人和政府）角度出发，研究综合运输网络下货流分配的优化问题。综合考虑运输费用、运输效率、碳排放量3方面因素，建立综合运输网络货流分配的线性优化模型，研究公铁运量目标比值和碳排放空间目标值等政策引导约束对货流分配的影响。算例结果表明：公铁运量目标比值变化能直接影响各路径货流分配量，同时显著影响碳排放量；碳排放空间目标值降低到82.152万t以下可提高铁路运输占比，但由于运输效率成本增加较快导致运输总成本增加，为保证运输总成本较低并尽可能减少碳排放量，需对二者进行权衡。关键词：综合运输；多式联运；货流分配；碳排放；系统最优中图分类号：U15文献标识码：A文章编号：2095-9931（2023）03-0107-09Optimization of Integrated Transportation Network Freight FlowDistribution Considering Carbon Tax CostLU Xia1,2,HUANG Junsheng1,2,MAO Baohua1,2,CHEN Shuo1,2(1.Integrated Transport Research Center of China,Beijing Jiaotong University,Beijing 100044,China;2.Key Laboratory of Transport Industry of Big Data Application Technologies for Comprehensive Transport,Beijing Jiaotong University,Beijing 100044,China)Abstract:To promote the transportation of bulk goods such as coal from road to railway,this paper studied the optimization of freight flow distribution in integrated transportation network from the perspective oftransportation participants(shippers,carriers,and government).The linear optimization model of freightflow distribution in the integrated transportation network was established by considering three aspects of交通运输研究TRANSPORT RESEARCH第9卷 第3期Vol.9No.3卢霞，黄俊生，毛保华，等.考虑碳税成本的综合运输网络货流分配优化J.交通运输研究，2023，9（3）：107-115.LU X,HUANG J S,MAO B H,et al.Optimization of integrated transportation network freight flow distribution consideringcarbon tax costJ.Transport Research,2023,9(3):107-115.DOI：10.16503/ki.2095-9931.2023.03.012107交通运输研究第9卷 第3期2023transportation cost,transportation efficiency and carbon emission.The impact of policy guidance constraints,such as the target ratio of road and rail transportation volume and the target value of carbonemission space,on the distribution of freight flow in the integrated transportation network were analyzed.The study case shows that the target ratio of road and rail transportation volume can directly affect the distribution of freight flow in each path and significantly affect the carbon emission.The reduction of the target value of carbon emission space for the whole transportation process to below821 520 t will increase the share of rail transportation,but the total transportation cost will increasedue to the faster increase of transportation efficiency cost.A trade-off between the two is needed to ensure the total transportation cost is low and minimize carbon emissions as much as possible.Key words:integrated transportation;multimodal transportation;freight flow distribution;carbonemission;system optimization0引言在新时期“双碳”目标指引下，积极调整运输结构，发展绿色交通体系，推进大宗货物运输“公转铁”成为交通运输领域发展的目标之一。以汾渭平原为例，国务院2018年印发的打赢蓝天保卫战三年行动计划1中将汾渭平原列为蓝天保卫战的“主战场”之一。京津冀及周边地区、汾渭平原2020-2021年秋冬季大气污染综合治理攻坚行动方案 提出2：提升铁路货运量，具有铁路专用线的大型工矿企业和新建物流园区的煤炭等大宗货物铁路运输比例原则上达到80%以上。根据2018年交通运输部在陕西省对汾渭平原开展的运输结构调研3，2017年陕西省货运公路、铁路、民航水运占比分别为 75.8%,24%,0.12%，即公铁比例约为4 1，铁路运能未能被充分挖掘。充分考虑和优化综合运输网络货流分配，能最大程度发挥铁路运输的作用，既利于推进运输结构调整，也益于推动公路和铁路的协调发展，减少碳排放和能源消耗，建立可持续的绿色物流系统。在促进大宗货物“公转铁”运输中，合理的货物输送和分配非常关键。综合运输方案中的货物运输路径和货流分配是相互关联且需要相互协调优化的，以达到最优的运输效能。国内外学者对综合运输方案中货物运输路径及货流分配的优化问题开展了一系列研究。部分学者重点研究了路径优化问题，如陈维亚等4从联运经营者的角度综合考虑了碳税、服务质量等因素，研究二者对多式联运路径选择的影响；冯芬玲等5从托运人的角度考虑班轮的班期、运输报价及货物的时间价值，构建了出海港选择及到港多式联运路径模型；李玉民等6考虑了多式联运参与者的差异化需求，建立了考虑运输时间、运输费用和碳排放的集装箱多式联运路径多目标优化模型；李海鹰等7为实现综合运输体系下各种运输方式的合理分工，以运输企业收益最大化为目标函数，建立了综合快运方案优化模型；尹传忠等8在综合运输体系下，从区域一体化的角度构建了运输全过程中多式联运系统协同优化模型。部分学者重点研究了货流分配优化问题，范鹏飞等9以大宗货物“公转铁”为主要方式研究了运输结构调整政策引起的货运需求变化，并对城市群干线网络货运需求分布与分配进行了预测；陈雷等10考虑了碳排放成本，构建了铁路与公路运输网络间的货流转移及流量分配综合优化模型；张政等11以产销两地的货运量为决策变量，建立了煤炭调运系统优化模型。少部分学者综合考虑了货物运输路径及货流分配问题的优化，如刘畅12在设计综合运输方案时，统筹考虑了公铁货流转移和铁路车流路径优化，以公转铁货流量等为决策变量，构建了二者的联合优化模型；刘星材等13在建设资金约束情况下，以运费、在建及新建费用最小为目标，研究了综合运输网络设计和货流分配问题。108TRANSPORT RESEARCHVol.9 No.3卢霞，等：考虑碳税成本的综合运输网络货流分配优化综上，既有研究在综合运输网络下货物运输路径优化和货流分配方面已经取得了较多研究成果。但与现有研究大多只考虑承运人的决策影响因素不同，本文将统筹考虑包括托运人、承运人和政府在内的所有运输参与方对货物运输的影响，同时增加碳税这个影响因素，构建综合运输网络下货流分配模型。在模型构建中还将引入政策引导约束，如碳排放空间14目标值约束和公铁运量目标比值约束，进一步考虑这些因素对货流分配的影响。最后，本文将以汾渭平原附近的运输通道为例，考虑在政策约束下以多式联运或直达运输方式运输煤炭等大宗货物，建立综合运输网络下货流分配的线性优化模型，采用 Lingo11.0软件求解，得到相应的货物运输路径和货流分配的综合运输方案，为促进煤炭等大宗货物运输“公转铁”提供一定的理论指导和实践支持。1问题描述本文主要研究多方式综合运输网络下的货流分配问题。假设需要将一批煤炭货物从煤矿所在地O运送至城市目的地D，在始发地和终到地之间有多种运输模式可供选择，包括铁路运输、公铁联运和公路直达运输。此外，运输路径中存在多个转运节点，其中铁路运输和公路运输在转运节点上相衔接，运输网络示意图见图1。本文从运输参与方，即托运人、承运人和政府的角度出发，在考虑运输费用、运输效率的基础上，增加碳税成本，旨在实现最低运输费用和时间的同时降低碳排放，从而达到社会效益和经济效益双赢的目标。本文还加入了两个政策引导约束，一是碳排放空间目标值，即运输全过程允许产生的最大碳排放量；二是公铁运量目标比值，即综合运输网络中公路运量与铁路运量允许的最大比值。2模型的建立2.1模型假设1）仅考虑运输途中的货损，不考虑中转过程及其他因素导致的货损；2）货物中转只能在节点处进行，且在同一节点处的换装中转至多一次，每个转运节点的货物中转能力充足；3）相邻节点间只能采用一种运输方式；4）运输全程速度不变，保持匀速；5）两种运输方式的发车时间不受限制。2.2符号定义模型参数及变量定义见表1。表1相关参数及变量定义lpi,i+1piipmTd-Tepvpvpmw节点i至节点i+1采用第p种运输方式的运输距离/km，i I，p P第p种运输方式碳排放因子/（kg/（tkm），p P节点i产生的单位中转费用/（元/t），i I节点i产生的单位中转碳排放量/（kg/t），i I第p种运输方式的单位公里货损率，p P货损率上限运到时间上限/h区域碳排放空间目标值/kg第p种运输方式的单位运价/（元/（tkm），p P第p种运输方式的速度/（km/h），p P第p种运输方式最大运输速度/（km/h），p P公铁运量比值，即公路运量/铁路运量货物总运量/t集合IP符号R1R2运输网络节点i的集合，i I运输方式p的集合，p P铁路运输公路运输参数起讫点转运节点公路运输铁路运输DO43521图1从煤炭产地到需求地的多方式运输网络109交通运输研究第9卷 第3期2023Qmtzsc中间变量C1C2C3决策变量wpi,i+1允许运量范围/t节点处货物中转时间/h时间成本系数/（元/（th）征收碳税单价/（元/kg）运输费用/元运输效率成本/元碳税成本/元节点i至节点i+1第p种运输方式运量/t，p P表1（续）参数2.3模型构建本文将运输效率、碳排放量转化为运输效率成本和碳税成本，建立使运输总成本最低的目标函数并结合货损率、运量、运输速度、运到时间、公铁运量目标比值以及碳排放空间目标值等约束确定货物运输路径和货流分配量。其中，运输总成本包括运输产生的费用、运输效率成本、碳税成本，这三者之间存在制约关系，如运输效率成本最低时，运输费用及碳税成本难以达到最低，故选择系统最优的方式使得运输总成本Z最低。（1）运输费用运输费用C1主要由货物在节点间采用第p种运输方式所花费的运输成本和在节点换装产生的中转费用构成，其计算公式为：C1=i Ip Pwpi,i+1(eplpi,i+1+i)（1）（2）运输效率成本运输效率涉及途中运输时间和节点处货物中转时间（不考虑除货物换装以外的途中滞留时间），本文将运输时间转化为运输效率成本 C2，其计算公式如下：C2=si Ip Pwpi,i+1()lpi,i+1vp+tz（2）（3）碳税成本运输碳排放量由运输途中的碳排放量及货物中转时产生的碳排放量两部分构成，碳税成本C3计算公式如下：C3=ci Ip Pwpi,i+1(plpi,i+1+i)（3）综上，各路径货流分配模型目标函数为：minZ=C1+C2+C3（4）2.4约束条件模型主要考虑了以下约束条件：i Ip Pplpi,i+1 m（5）i Ip Pwpi,i+1()plpi,i+1+i-T（6）wpi,i+1 Qmi I,p P（7）i Ip PR1wpi,i+1i Ip PR2wpi,i+1（8）lpi,i+1vp+tz Tdi I,p P（9）vp vpmp PR1（10）vp vpmp PR2（11）上述约束条件中，式（5）表示运输过程中的货损率满足货损要求；式（6）表示运输产生的碳排放量满足碳排放空间目标值约束；式（7）表示每个区段的运量在允许运量范围内；式（8）表示公铁运量比值满足目标约束；式（9）表示整个运输过程的运输效率满足运到时间约束；式（10）、式（11）表示公路、铁路运输速度满足各自运输速度约束。3算例分析3.1算例及模型参数该算例以汾渭平原附近的3 条铁路线路BZ、DZ、ZC铁路以及既有平行公路线路构建运输网络。以两大主要煤矿ZGE矿区和SD矿区为起点，由公铁联运或公路直达运输两种方式将两大煤矿的煤炭运送至XZ和DT两个目的地，其中公铁联运货物的中转在铁路站点进行衔接，如图2所示。承运人一年内将 7 000万 t煤炭从 SD矿区和ZGE矿区运出，其中SD矿区发送3 000万t，ZGE矿区发送 4 000 万 t，到达 XZ 目的地 2 000 万 t、DT目的地5 000万t，网络简化图如图3所示，图110TRANSPORT RESEARCHVol.9 No.3中1、2分别为SD矿区、ZGE矿区两个起点，5、6分别为XZ目的地、DT目的地两个讫点，任意两节点间的数字表示距离。根据参考文献4，煤炭货物碳税单价取0.15元/kg。假设这批货物的运送次数是5 000次，根据我国对超过0.5%的货损需要向保险公司索赔的要求15，运输过程中要保证每次货损率均值低于0.5%。为保证公铁运输比例，参考文献16设置公铁运量目标比值为 1.0；碳排放空间目标值为85万t。参考文献4的时间窗约束，本文的最大运到时间约束为48h。为将运输效率转化为成本费用，参考文献17，将时间成本系数定为 1 元/（th）。单位运价、运输速度、碳排放因子、运量范围、货损率以及运输中转相关的参数取值如表2、表3所示。表2不同运输方式的相关参数取值运输方式铁路运输公路运输单位运价/（元/（tkm）60.200.35运速/最大运速/（km/h）455/7080/100碳排放因子/（kg/（tkm）1.3710-37.810-2运量范围/万t30089003008900每百公里平均货损率（%）40.0150.01注：铁路运输碳排放因子由SH公司公布的数据计算得到；公路运输碳排放因子为载重30 t的重型柴油货车的碳排放因子。表3运输中转的相关参数取值参数单位中转费用/（元/t）中转时间/h单位中转碳排放/（kg/t）6451.3710-3取值注：中转时间参考文献4和18的中转时间进行取值；单位中转碳排放参考装卸设备装卸单位重量货物碳排放以及铁路运输的碳排放因子计算得到。3.2结果分析该模型是一个线性模型且存在最优解，采用数学规划软件 Lingo 11.0对算例进行求解，货物运输路径及运量分配方案如表4所示。根据表4结果可知：运输总成本中，运输费用占比较大。在上述参数取值下，公路与铁路承担运量比例相当，实际公铁运量比值为1.00，采ZGE矿区XJWDDG128 kmBTTSD矿区264 km120 km245 km184 kmYXG144 km67 kmYZDT目的地452 km90 kmSCNXZ目的地铁路运输公路运输矿区站点355 km图2煤炭运输网络公路运输铁路运输90 km452 km355 km264 km128 km184 km41120 km144 km245 km67 km632图3网络简化图5卢霞，等：考虑碳税成本的综合运输网络货流分配优化111交通运输研究第9卷 第3期2023用公铁联运与采用公路直达运输方式的运量之比约为2 1。若整个运输过程仅采用公路直达运输，运输费用、运输效率成本及碳税成本分别为79.030亿元、2.823亿元、2.642亿元，运输费用及碳税成本均上升，碳排放量将上升至176.124万t。从结果对比来看，采用公铁联运与公路运输相结合的方式，碳排放量降低了53.4%。这表明，相比于全程采用公路直达运输的方式，采用公铁联运与公路运输相结合的综合运输方案可以显著降低碳排放量。3.3灵敏度分析（1）公铁运量目标比值在原有参数不变的情况下，调整公铁运量目标比值的大小，从1.00逐步下降至0.86，研究公铁运量目标比值约束对各路径货流分配的影响，结果如表5所示。可以看出，公铁运量目标比值变化会直接影响货流分配。表5公铁运量目标比值变化情况下的货流分配单位：万t货物运输路径1-2-4-51-2-3-62-4-52-3-61-51-62-52-6变化因子：公铁运量目标比值1.003002 1003001 9003003001 1007000.953002 1003002 268.43003001 100331.60.903002 100677.82 300300300722.23000.863002 1001 039.52 300300300360.5300注：由于约束限制，公铁运量目标比值最低为0.86。图4所示为公铁运量不同目标比值对应的各项指标的变化率。其中，变化率均为相对初始值1.00 的变化，即变化率=（现值-初始值）/初始值100%。由图4可以看出，碳税成本与碳排放量变化率趋势线几乎重合，说明碳税成本与碳排放量保持同步增长，碳排放量越高，碳税成本越高。碳排放量与公铁运量目标比值呈正相关，随着公铁运量目标比值的下降，碳排放量不断下降。在公铁运量目标比值下降的过程中，运输效率成本增加，碳税成本下降，运输总成本增加。其中，运输费用呈先降后升趋势，这是因为采用公铁联运方式的货运量达到一定值后，节点换装产生的中转费用较高，公铁联运方式的单位运输费用高于采用公路直达运输方式的单位运输费用，导致出现铁路运输占比增加，运输费用却上升的情况。当公铁运量目标比值低于 0.95 时，铁路运量增加，运输费用也增加；在公铁运量目标比值约为0.91时开始出现运输费用变化率大于零的情况，即此时运输费用超过初始方案；当公铁运量目标比值达到0.86时，运输费用、运输总成本相比初始方案分别增长0.425%,1.282%。公铁运量目标比值降低能显著降低碳排放量，但会导致运输总成本的增加，因此为增加铁路运输的吸引力，需灵活调整单位运价。（2）碳排放空间目标值在其他参数固定不变的情况下，改变碳排放空间目标值的大小，其变化范围为70万85万t，探讨碳排放空间目标值约束对各路径货流分配的影响，结果如表6所示。表4货物运输路径及运量分配方案铁路运输速度/（km/h）55货物运输路径1-2-4-5（公铁联运）1-2-3-6（公铁联运）2-4-5（公铁联运）2-3-6（公铁联运）1-5（公路运输）1-6（公路运输）2-5（公路运输）2-6（公路运输）货流分配/万t3002 1003001 9003003001 100700总成本/亿元77.207C1/亿元69.582C2/亿元6.393C3/亿元1.232碳排放量/万t82.152公铁运量比值1.00112TRANSPORT RESEARCHVol.9 No.3表6碳排放空间目标值变化情况下的货流分配单位：万t货物运输路径1-2-4-51-2-3-62-4-52-3-61-51-62-52-6变化因子：碳排放空间目标值/万t853002 1003001 9003003001 100700803002 1003002 059.13003001 100540.975432.51 967.53002 432.5300300967.530070812.61 587.43002 812.6300300587.4300图5所示为不同碳排放空间目标值对应的各项指标的变化率。由图5可知，随着碳排放空间目标值的逐渐减小，公铁运量比值不断下降，铁路货运量占比提高。其中，运输费用先降后升，运输效率成本不断增加，碳税成本逐渐下降，运输总成本增加。在碳排放空间目标值为85万t的方案中，实际碳排放量为 82.152万 t，低于约束值，因此碳排放空间目标值在82.152万85万t之间时，最优方案保持不变，各参数变化率为零。当碳排放空间目标值降低至 82.152 万 t 以下时，虽然运输费用、碳税成本在不断降低，但由于运输效率成本增速较快，导致运输总成本增加。碳排放空间目标值从74万t下降到73万t时，运输费用变化率由负变正，运输费用高于初始方案。随着碳排放空间目标值的降低，运输总成本变化率逐渐上升，当碳排放空间目标值降至70万t时，运输总成本增长率达到0.921%。不难看出，碳排放空间目标值约束越严苛，运输总成本越高，且呈现持续增长趋势。4结论本文从运输参与方的角度出发，运输总成本考虑了运输费用、运输效率成本以及碳税成本3部分，构建了综合运输网络货流分配模型，得到了最优综合运输方案。算例分析表明：（1）在综合运输网络中，公铁运量目标比值会直接影响货流分配。增加铁路运量比例可显著运输效率成本C2碳税成本C3碳排放量总成本Z运输费用C120100-10-20运输费用C1变化率（%）0.60.40.20.0-0.2-0.4-0.6-0.8公铁运量目标比值0.98 0.96 0.940.92 0.90 0.88 0.86碳税成本C3碳排放量-18变化率（%）公铁运量目标比值0.8700.8650.860公铁运量目标比值1.000.980.960.940.920.900.880.86变化率（%）图4公铁运量目标比值变化的影响卢霞，等：考虑碳税成本的综合运输网络货流分配优化113交通运输研究第9卷 第3期2023降低碳排放量，降低运输费用。但采用公铁联运方式的货流量到达一定值后，继续提高其运量占比将导致运输费用增加。当公铁运量目标比值低于0.95时，铁路运量增加，运输费用也增加；当公铁运量目标比值约为0.91时，运输费用将高于初始方案。因此，应当适当调整铁路运价以增强铁路运输的竞争力。对于综合运输方案而言，应在适当规模的货流中采用公铁联运，以使得经济成本和环保成本最低。（2）在煤炭等大宗货物运输中，通过政策引导约束降低碳排放空间目标值，可以有效降低碳排放量。当碳排放空间目标值降低至82.152万 t以下时，虽然运输费用、碳税成本在不断降低，但由于运输效率成本增速较快，导致运输总成本增加。且碳排放空间目标值约束越严苛，运输总成本越大，为保证运输总成本较低且满足碳排放量尽可能小，在设计综合运输方案时需要对二者进行权衡。本研究可以在以下方面做进一步拓展：当路网节点及线路数目增多时，需进一步完善模型和算法；要实现大宗货物运输“公转铁”目标，还需在原有网络上进行网络重构，研究网络如何调整才能有效增加铁路货运吸引力。参考文献1国务院.国务院关于印发打赢蓝天保卫战三年行动计划的通知（国发 2018 22 号）EB/OL.（2018-07-03）2023-04-15.http:/ RESEARCHVol.9 No.3诺的多式联运路径优化J.铁道科学与工程学报，2022，19（1）：34-41.5冯芬玲，孙楠佳.考虑时间价值的中非多式联运路径与出海港选择J.交通运输系统工程与信息，2022，22（2）：45-53.6李玉民，郭晓燕，杨露.考虑多目标的中欧集装箱多式联运路径选择J.铁道科学与工程学报，2017，14（10）：2239-2248.7李海鹰，曹玥，廖正文，等.基于服务网络的综合快运方案优化J.中国铁道科学，2020，41（4）：136-145.8尹传忠，邱慧妍，柯媛定，等.区域主枢纽港多式联运网络协同优化J.铁道科学与工程学报，2022，19（1）：63-70.9范鹏飞，宋国华，雷雪.城市群干线网络货运需求预测模型以京津冀为例J.交通运输研究，2021，7（6）：50-60.10陈雷，林柏梁，温旭红，等.低碳运输下陆路运输系统货流转移及流量分配综合优化J.东南大学学报（自然科学版），2016，46（3）：671-674.11张政，韦伟，史芮嘉，等.考虑动态定价的煤炭调运方案优化研究J.交通运输系统工程与信息，2017，17（4）：105-110.12刘畅.货物运输“公转铁”与铁路车流径路联合优化模型J.铁道学报，2020，42（12）：18-27.13刘星材，何世伟，孙杨，等.基于云和声搜索的综合运输网络设计和货流分配优化J.吉林大学学报（工学版），2016，46（1）：108-113.14黄俊生，毛保华，吴雪妍.碳中和战略下我国交通运输行业碳减排策略研究J.北京交通大学学报（社会科学版），2023，22（2）：107-116.15梁金桂.国际贸易中海运货差探析J.中国海商法年刊，2010，21（3）：42-47.16王超，怀旭.公铁联运下我国中西部地区低碳运输实证研究J/OL.交通运输系统工程与信息，2023.http:/ 2016 232号）Z.北京：交通运输部，外交部，国家发展改革委，等，2016.15魏姗姗，徐畅，陈秀岭.山东省多式联运“一单制”发展现状与对策J.物流技术，2022，41（1）：42-45.16黄敏珍，王瑞民，林晓蕾.多式联运数据交换区块链技术应用研究J.铁道运输与经济，2021，43（2）：75-81.17李玮，张改平.多式联运“一单制”实践探索J.综合运输，2021，43（10）：107-113.18陈静，潘庆全.铁路提单制度的理论争点与实践困境EB/OL.2021-09-01（2023-03-28）.https:/