扩展功能
文章信息
- 胡军红, 汤瑞, 闻成维, 张雯婕.
- HU Junhong, TANG Rui, WEN Chengwei, ZHANG Wenjie
- 运量需求不确定下考虑混合时间窗与货损的多式联运路径选择
- Multimodal transport path selection considering mixed time windows and cargo damage with uncertain demand
- 公路交通科技, 2025, 42(8): 196-206
- Journal of Highway and Transportation Research and Denelopment, 2025, 42(8): 196-206
- 10.3969/j.issn.1002-0268.2025.08.021
-
文章历史
- 收稿日期: 2023-02-24
作为一种先进的运输组织形式,多式联运综合利用铁路、公路、水路等不同运输方式的优势,为客户提供灵活、高效、多样的运输服务。随着人们生活水平上升,人们对于多式联运运输质量的要求在不断地提高,而时效性与安全性是影响运输质量的关键因素。多式联运适用于长距离运输,若某条线路的货损率过高,则会产生逆向物流增加物流成本,使多式联运经济性和时效性优势发挥不够显著。因此,应合理地组合运输方式与运输线路,以尽可能低的成本在规定的时间窗运达并控制质量损失。
在运输时效性研究方面,较多的学者会考虑班轮班期的限制,如吕学伟[1]将运输方式的班期设置为硬时间窗,收货时间设置软时间窗,说明时间窗问题在运输模型中的影响。基于此,Guo[2]解决了带有提取货时间、时间窗等特征的多式联运路径选择问题。温廷新[3]综合考虑配送车辆的碳排放量及货损问题,以带软时间窗的车辆路径优化模型为基础,构建了一个配送成本及能耗最低的路径优化模型。刘畅[4]建立时间价值成本模型,在此基础上构建以综合物流成本最低为目标的路径选择模型。裴骁[5]考虑列车等待时间,建立了由列车、船舶时刻表和目的地软时间窗组成的混合时间窗为约束的综合满意度模型。孙岩[6]考虑中转节点固定作业时间窗及运输时限软时间窗,构建费用最小化为优化目标的模型。李春发[7]考虑客户需求量、车辆最大承载量和客户时间窗要求等影响因素,构建了一个以油耗成本、惩罚成本等综合配送成本最小为目标函数的数学优化模型。
在考虑货损货差方面,大多是研究冷藏集装箱的货物损失或是具有高时效性的货物价值损失,这种类型货物的货损与运输时间有关[8-11],时间越长则会导致货物的价值损失越大。有少部分文献探究了货物因装载、运输等引起的受湿、沾污、碰撞破裂、丢失等造成的损失,这部分货损与运输距离有关。如陈维亚[12]引入运到时限和货物完好程度2个质量指标,构建由运输成本、碳税成本和质量损失补偿成本最小的路径优化模型,其中质量损失补偿成本与每百公里单位货损率相关。
以上大部分文献,不论是研究混合时间窗或是研究货损问题,都是在运输需求确定的背景下进行分析,然而实际运输时由于突发事件的影响[13],运输需求可能会产生波动。而在运输需求变化下,货物的运输时间以及货损成本都会存在一定的变动,导致运输路径发生变化。针对现有研究在不确定需求下考虑混合时间窗与货损率限制变化不足的问题,本研究在不确定运输需求下,以运输方式固定作业时间窗及目的地软时间窗组成的混合时间窗、货损率限制为约束,构建以运输成本、货损成本及时间成本综合最低为目标的路径选择模型,提高运输服务的质量。
1 问题描述多式联运经营人要将一批需求不确定的零散高价值货物从起始地O运送至目的地D,起终点间经过多个运输节点。考虑到运输量虽难以确定,但可以通过以往经验得到一个变化范围,三角模糊数是表达判断比较时,给出的模糊量范围:最低可能值、最可能值、最高可能值。因此,使用三角模糊数表示不确定需求量较为恰当。公路运输较为灵活,无时间约束。当采用铁路或水路运输时,受班次时刻表限制,具体流程如图 1所示。A,a代表节点铁路或水路某一班次的到达和离开时间,B,b代表下一班次的到达和离开时间。货物到达节点进行下一阶段运输前会出现3种情景:(1)在A之前到该节点。需要等A到达后完成换装并等待a离开。(2)在A之后到达该节点。直接进行换装并等待a离开。(3)在A之后到达该节点,但要完成换装则会晚于a的离开时间。因此需要等待下一班次B到达后进行换装,并等待下一班次b离开。
|
| 图 1 铁路、水路班次运行示意图 Fig. 1 Schematic diagram of railway and waterway schedules |
| |
为保证运输质量,本研究考虑不同运输方式固定运输时刻表的限制,在一定货损率约束下寻求多式联运总成本最小的路径。为便于建模,作出以下假设:
(1) 运输过程中货物不可分开,两个节点间只有一种运输方式。(2)每个节点货物最多只能经过一次。(3)仅考虑运输与中转时的货损。(4)每个节点之间各种运输方式的承载力以及节点的中转能力已知。(5)节点最多有3种运输方式:公路、铁路、水路,不考虑其他运输方式。
2 模型构建 2.1 符号说明对于多式联运运输网络G= (A,I,K),I为城市节点集合,I={i,j|i,j=1,2,3…n};K为运输方式集合,K={k,l|k,l=1,2,3},其中:1, 2, 3, 分别代表公路、铁路、水路。
为保证多式联运的运输质量控制成本,本研究将运输时间和运输货损转换为时间成本以及货损成本,以运输成本、时间成本、货损成本综合最小构建目标函数。
2.2.1 成本构成(1) 运输成本
货物的运输成本与运输工具与运输距离有关,主要包括在途成本以及中转成本,如式(1)所示。
| $ C_{\mathrm{ys}}=\sum\limits_{i, j \in I} \sum\limits_{k \in K} \tilde{Q} \cdot f_k \cdot d_{i j}^k \cdot x_{i j}^k+\sum\limits_{i \in I} \sum\limits_{k, l \in K} \tilde{Q} \cdot f_i^{k l} \cdot y_i^{k l}。$ | (1) |
(2) 时间成本
多式联运的时间成本与运输总时间有关,运输总时间由在途运输时间、中转时间、衔接等待时间构成。若运输总时间早于目的地时间窗则有储存成本,若晚于目的地时间窗则有惩罚成本,其中在途运输时间与运输距离与运输速度有关,具体表达式如下:
| $ t_{i j}^k=d_{i j}^k / \bar{v}_k, $ | (2) |
| $ T=\sum\limits_{i, j \in I} \sum\limits_{k, l \in K}\left(t_{i j}^k \cdot x_{i j}^k+t_i^{k l} y_i^{k l}+t_{i 1}^k \cdot x_{i j}^k+t_{i 2}^k \cdot x_{i j}^k\right) , $ | (3) |
| $ \begin{gathered} C_{\mathrm{sj}}=\tilde{Q} \cdot f_{\mathrm{t}} \cdot \max (E T-T, 0)+\tilde{Q} \cdot P_{\mathrm{t}} \cdot \\ \max (T-L T, 0)。\end{gathered} $ | (4) |
(3) 货损成本
货损成本包括货物在运输和中转过程的质量损失成本,如式(5)所示:
| $ C_{\mathrm{hs}}=\sum\limits_{i, j \in I} \sum\limits_{k \in K} \tilde{Q} \cdot U \cdot\left(x_{i j}^k \cdot\left(d_{i j}^k / 100\right) \cdot \theta_k+\theta_z \cdot y_i^{k l}\right) \text { 。} $ | (5) |
| $ \min Z=C_{\mathrm{ys}}+C_{\mathrm{sj}}+C_{\mathrm{hs} }。$ | (6) |
约束条件:
| $ \sum\limits_{k \in K} x_{i j}^k \leqslant 1, \quad \forall i, j \in I, $ | (7) |
| $ \sum\limits_{k \in K} \sum\limits_{l \in K} y_i^{k l} \in 1, \quad \forall i \in I, $ | (8) |
| $ \sum\limits_{i \in I} \sum\limits_{k \in K} x_{i j}^k-\sum\limits_{h \in I} \sum\limits_{k \in K} x_{j h}^k= \begin{cases}1, & j=d \\ 0, & j \neq o, j \neq d \\ & j \in(o, d) \\ -1, & j=o\end{cases}, $ | (9) |
| $ t_{i 1}^k= \begin{cases}e_{i j}^k-t_i, & t_i<e_{i j}^k \\ 0, & e_{i j}^k<t_i<t_i+t_i^{k l} \leqslant l_{i j}^k \\ e_{i j}^{k^{\prime}}-t_i, & l_{i j}^k<t_i+t_i^{k l}\end{cases}, $ | (10) |
| $ t_{i 2}^k= \begin{cases}l_{i j}^k-e_{i j}^k-t_i^{k l}, & t_i<e_{i j}^k \\ l_{i j}^k-t_i-t_i^{k l}, & e_{i j}^k \leqslant t_i<t_i+t_i^{k l} \leqslant l_{i j}^k, \\ l_{i j}^{k^{\prime}}-e_{i j}^{k^{\prime}}-t_i^{k l}, & l_{i j}^k<t_i+t_i^{k l}\end{cases} $ | (11) |
| $ t_i^k=t_{i 1}^k+t_{i 2}^k=\left\{\begin{array}{cc} l_{i j}^k-t_i-t_i^{k l}, & t_i<e_{i j}^k \\ l_{i j}^k-t_i-t_i^{k l}, & e_{i j}^k \leqslant t_i<t_i+t_i^{k l} \leqslant l_{i j}^k \\ l_{i j}^{k^{\prime}}-t_i-t_i^{k l} & l_{i j}^k<t_i+t_i^{k l} \end{array}\right., $ | (12) |
| $ \sum\limits_{i, j \in I} \sum\limits_{k \in K}\left[x_{i j}^k \cdot\left(d_{i j}^k / 100\right) \cdot \theta_k+\theta_z \cdot y_i^{k l}\right] \leqslant \theta, $ | (13) |
| $ x_{i j}^k \in\{0, 1\}, y_i^{k l} \in\{0, 1\}, \quad \forall i, j \in I, \quad \forall k \in K, $ | (14) |
| $ \tilde{Q} \cdot x_{i j}^k \leqslant Q_{i j}^k, \quad \forall i, j \in I, \quad \forall k \in K, $ | (15) |
| $ \tilde{Q} \cdot y_i^{k l} \leqslant Q_i^{k l}, \quad \forall i, j \in I, \quad \forall k, l \in K。$ | (16) |
式(7)为运输货物时节点间只可选择一种运输方式;式(8)保证货物在同一个节点发生一次转载;式(9)为节点货运量保持平衡;式(10)为等待运输工具到达的时间;式(11)为等待运输工具离开的时间;式(12)为总等待时间;式(13)为多式联运全过程的总货损率满足货损率的限制;式(14)为决策变量约束;式(15)为所选运输方式的运能大于货运量;式(16)为所选节点的中转能力大于货运量。
3 模型求解 3.1 模糊机会约束规划为了清晰表达所构建模型中的模糊变量,需要将模型转换为确定的形式。本研究主要借鉴1996年Liu[14]首次提出的不确定理论规划方法。含有模糊参数的模型只要做出满足模糊机会约束决策的概率在一定的置信度之上,就可以得到清晰的模型方程。具体表示为以下的形式:
目标函数:
| $ \min f(x, \tilde{Q}) \text { 。} $ | (17) |
约束条件:
| $ {Pos}\{f(x, \tilde{Q}) \leqslant \bar{f}\} \geqslant \alpha, $ | (18) |
| $ {Pos}\left\{\tilde{Q} \cdot x_{i j}^k \leqslant Q_{i j}^k\right\} \geqslant \beta_1, $ | (19) |
| $ {Pos}\left\{\tilde{Q} \cdot y_i^{k l} \leqslant Q_i^{k l}\right\} \geqslant \beta_2 。$ | (20) |
式(18)为满足置信水平α时,f所能取得的最小值;式(19)为节点i到节点j上的货运量小于所选运输方式的运输能力的可能性高于置信水平β1;式(20)为在节点i中转的货运量小于所选运输方式的中转能力的可能性高于置信水平β2。
引入以下定理进一步将模型清晰化,设有三角模糊数
| $ \begin{gathered} f \geqslant\left[(1-\alpha) Q_{\mathrm{E}}+\alpha Q_{\mathrm{M}}\right] \cdot \\ \left\{\sum\limits_{i, j \in I} \sum\limits_{k \in K} f_k \cdot d_{i j}^k \cdot x_{i j}^k+\sum\limits_{i \in I} \sum\limits_{k, l \in K} f_i^{k l} \cdot y_i^{k l}+\right. \\ f_{\mathrm{t}} \cdot \max (E T-T, 0)+P_{\mathrm{t}} \cdot \max (T-L T, 0)+ \\ \left.\sum\limits_{i, j \in I} \sum\limits_{k \in K} U \cdot\left(x_{i j}^k \cdot\left(d_{i j}^k / 100\right) \cdot \theta_k+\theta_z \cdot y_i^{k l}\right)\right\}, \end{gathered} $ | (21) |
| $ Q_{i j}^k \geqslant\left[\left(1-\beta_1\right) \cdot Q_{\mathrm{E}}+\beta_1 Q_{\mathrm{M}}\right] \cdot x_{i j}^k, $ | (22) |
| $ Q_i^{k l} \geqslant\left[\left(1-\beta_2\right) \cdot Q_{\mathrm{E}}+\beta_2 Q_{\mathrm{M}}\right] \cdot y_i^{k l}。$ | (23) |
遗传算法具有较强的全局搜索能力,但是局部搜索能力较差,而模拟退火算法具有较强的局部搜索能力。在两种算法思想的基础上,以遗传算法为主,引入模拟退火的思想,以一定概率接受劣解,避免种群快速收敛陷入局部最优。
3.2.1 遗传算法(1) 染色体编码
采用自然数编码方式,将节点序号和运输方式间隔排列编码,每条编码串的奇数位表示节点偶数位的运输方式。将运输方式表示为1,2,3,其中1表示公路,2表示铁路,3表示水路。如图 2所示,染色体编码实际意义表示:起点1-铁路-3-水路-5-公路-7-水路-终点10。
|
| 图 2 染色体编码示意图 Fig. 2 Schematic diagram of chromosome encoding |
| |
(2) 适应度函数计算
本研究所建模型中以物流总成本最小化为目标函数,因此适应度值采用目标函数值的倒数,即目标函数值越小,适应度函数值越大,表明可行性越高。其二者关系表达式如式(24)所示:
| $ f_i=\frac{1}{Z_i}, i=1, 2, \cdots, n, $ | (24) |
式中,fi为个体的适应度值;Zi为第i个个体对应的物流总成本;n为种群的规模数。
(3) 选择
采用轮盘赌的择优机制选择染色体,假设∑fi表示这一代种群的总适应度值,fi为种群中第i条染色体的适应度值,则第i条染色体被选中的概率为fi/∑fi。
(4) 交叉
本研究采用单点交叉法,通过对两个染色体位置的任意选择,将任意选择位置的右边互换,以获得两个完全独立的子代,如图 3所示。通过轮盘赌策略选出两个待操作的父代1和父代2,并生成两个介于0至1之间的随机数p1和p2。若p1小于交叉概率,则对父代1和父代2进行单点交叉,生成新的子代1和子代2,否则,父代1和父代2直接成为新的子代1和子代2。
|
| 图 3 交叉操作示意图 Fig. 3 Schematic diagram of crossover operation |
| |
(5) 变异
采用单点变异,根据交叉操作生成的0-1之间的随机数p2,判断p2与变异概率的大小,若p2小于变异概率,则对新生成的子代1和子代2进行变异。
3.2.2 模拟退火算法在遗传算法进行变异操作后,加入模拟退火算子。设置初始温度T0、终止温度Tend、降温速率q、链长l,并计算适应度值f(j)。若f(j)<f(i),则通过精英保留策略执行选择算子,判断迭代次数k是否满足要求,否则将个体的适应度值从大到小进行排序,选择一部分优秀的个体进行降温操作,并利用Metropolis法则判断是否接受新解,判断迭代次数k是否满足要求。若k满足要求则判断此时温度T是否满足Tend,若满足则输出最优解,不满足则进行降温操作,T=T0q;若k不满足要求,则进行迭代次数加1,继续从遗传算法的选择操作开始,重新选择染色体,直至输出最优解。
基于模拟退火的遗传算法具体流程如图 4所示。
|
| 图 4 基于模拟退火的遗传算法流程 Fig. 4 Genetic algorithm process based on simulated annealing |
| |
4 算例分析 4.1 算例参数
假设准备运输一批运量模糊的货物,运输需求量用三角模糊数
|
| 图 5 算例网络结构示意图 Fig. 5 Schematic diagram of example network structure |
| |
参考Sun[16]算例中各个路段上不同运输方式之间的运输距离的取值,以及各种运输方式能力及中转节点运输能力的取值,具体如表 1,表 2所示。由于本研究考虑货物以t为计算单位,参考Koza[17]研究发现,33 m3集装箱的平均货物质量为11 t,因此将算例中的运能相关数值同时扩大11倍取值。表 3为各节点处铁路与水路的运输时间窗,默认公路没有时间窗限制。目的地软时间窗为[45 h,70 h],提到到达终点的仓储成本为15元/(t·h),延迟到达终点的惩罚成本为30元/(t·h)[18]。表 4为不同运输方式间的中转费用、中转时间参数。
| 运输路段 | 运输距离(km)/运输能力(t) | 运输路段 | 运输距离(km)/运输能力(t) | |||||
| 公路 | 铁路 | 水路 | 公路 | 铁路 | 水路 | |||
| (1,2) | 151/583 | 171/803 | 187/693 | (18,20) | 137/726 | — | — | |
| (1,3) | 204/792 | 201/836 | — | (19,20) | 180/616 | — | — | |
| (1,4) | 190/473 | 196/895 | 189/803 | (19,23) | 146/418 | — | — | |
| (2,7) | 218/781 | 231/847 | 236/671 | (19,24) | 185/660 | — | — | |
| (2,8) | 190/462 | 193/715 | 172/528 | (20,22) | 189/561 | — | — | |
| (3,6) | 173/858 | 150/484 | — | (20,23) | 195/693 | 164/638 | 188/781 | |
| (3,7) | 155/517 | — | — | (21,22) | 165/737 | — | — | |
| (4,5) | 180/671 | 225/882 | 208/916 | (21,27) | 160/1 280 | 147/881 | 169/952 | |
| (5,6) | 145/803 | — | — | (22,23) | 143/583 | — | — | |
| (5,12) | 230/671 | 262/906 | 251/914 | (22,26) | 181/473 | 139/517 | — | |
| (6,11) | 171/594 | 156/869 | — | (22,27) | 135/451 | — | — | |
| (7,10) | 214/528 | 237/858 | 232/572 | (23,25) | 203/836 | 165/638 | 192/682 | |
| (7,11) | 203/792 | — | — | (24,25) | 203/418 | 185/726 | — | |
| (8,9) | 198/649 | 189/792 | 204/352 | (24,30) | 220/660 | — | — | |
| (8,10) | 175/847 | — | — | (24,31) | 205/473 | 195/770 | 211/671 | |
| (9,13) | 182/836 | 154/517 | 185/847 | (25,26) | 164/693 | 148/561 | — | |
| (9,14) | 183/528 | 182/869 | 194/605 | (25,29) | — | 188/737 | — | |
| (10,11) | 196/748 | — | — | (25,30) | 216/737 | 186/726 | 203/737 | |
| (10,14) | 205/748 | 219/858 | 210/517 | (26,27) | 143/583 | — | 146/748 | |
| (11,12) | 158/759 | — | — | (26,29) | 177/374 | 132/605 | — | |
| (11,15) | 190/638 | 206/451 | — | (27,28) | 202/451 | 170/860 | 200/969 | |
| (12,16) | 190/374 | 241/814 | 268/897 | (28,29) | 207/836 | — | — | |
| (13,18) | 125/693 | — | — | (28,33) | 180/418 | 154/594 | — | |
| (13,19) | 173/627 | — | — | (28,35) | 375/781 | 179/993 | 193/805 | |
| (13,24) | 158/748 | 178/770 | 172/462 | (29,30) | 146/737 | — | — | |
| (14,15) | 213/847 | 216/770 | 216/583 | (29,32) | 200/847 | — | — | |
| (14,17) | 133/506 | — | — | (29,33) | 174/374 | 146/770 | — | |
| (14,18) | 177/638 | — | — | (30,31) | 211/693 | — | — | |
| (14,20) | 200/583 | 186/726 | 195/770 | (30,32) | 199/572 | 190/660 | — | |
| (15,16) | 173/671 | 137/627 | 161/385 | (31,32) | 165/858 | 139/748 | 168/781 | |
| (15,17) | 202/517 | 191/528 | — | (32,34) | 176/385 | 156/858 | 164/627 | |
| (16,21) | 170/948 | 123/851 | 164/847 | (32,35) | 138/890 | 261/998 | 199/875 | |
| (17,20) | 199/506 | — | — | (33,34) | 190/759 | 195/506 | — | |
| (17,21) | 151/770 | — | — | (33,35) | 187/638 | 154/781 | — | |
| (17,22) | — | 180/484 | — | (34,35) | 108/847 | — | — | |
| 节点 | 公-铁 | 公-水 | 铁-水 | 公-公 | 铁-铁 | 水-水 | 节点 | 公-铁 | 公-水 | 铁-水 | 公-公 | 铁-铁 | 水-水 |
| 2 | 781 | 517 | 715 | 803 | 506 | 462 | 19 | — | — | — | 847 | — | — |
| 3 | 825 | — | — | 649 | 726 | — | 20 | 770 | 748 | 693 | 605 | 385 | 759 |
| 4 | 996 | 792 | 915 | 727 | 818 | 981 | 21 | 858 | 748 | 918 | 794 | 958 | 974 |
| 5 | 836 | 873 | 906 | 707 | 963 | 941 | 22 | 693 | — | — | 517 | 330 | — |
| 6 | 682 | — | — | 803 | 473 | — | 23 | 352 | 638 | 605 | 825 | 759 | 759 |
| 7 | 385 | 462 | 352 | 671 | 693 | 418 | 24 | 792 | 374 | 858 | 528 | 781 | 594 |
| 8 | 484 | 836 | 572 | 528 | 704 | 682 | 25 | 847 | 363 | 517 | 396 | 803 | 572 |
| 9 | 627 | 517 | 539 | 616 | 737 | 737 | 26 | 704 | 627 | 649 | 759 | 374 | 572 |
| 10 | 858 | 440 | 748 | 550 | 583 | 682 | 27 | 748 | 759 | 951 | 739 | 950 | 995 |
| 11 | 858 | — | — | 374 | 374 | — | 28 | 737 | 847 | 948 | 762 | 973 | 905 |
| 12 | 718 | 771 | 929 | 762 | 951 | 927 | 29 | 550 | — | — | 550 | 770 | — |
| 13 | 869 | 594 | 594 | 396 | 836 | 495 | 30 | 693 | 638 | 605 | 385 | 572 | 781 |
| 14 | 858 | 528 | 572 | 429 | 418 | 737 | 31 | 429 | 583 | 715 | 407 | 836 | 770 |
| 15 | 594 | 792 | 682 | 462 | 781 | 429 | 32 | 715 | 341 | 825 | 847 | 429 | 682 |
| 16 | 770 | 949 | 915 | 761 | 927 | 904 | 33 | 352 | — | — | 858 | 473 | — |
| 17 | 407 | — | — | 352 | 880 | — | 34 | 781 | 418 | 627 | 594 | 363 | — |
| 18 | 561 | — | — | 825 | — |
| 节点 | 铁路运输时刻 | 水路运输时刻 | ||||||
| 1 | 2 | 3 | 4 | 1 | 2 | 3 | ||
| 2 | [4, 7] | [9, 12] | [14, 17] | [19, 22] | [4, 9] | [11, 16] | [18, 23] | |
| 3 | [2, 5] | [7, 10] | [12, 15] | [17, 20] | — | — | — | |
| 4 | [3, 6] | [8, 11] | [13, 16] | [18, 21] | [5, 10] | [12, 17] | [19, 24] | |
| 5 | [2, 5] | [6, 9] | [11, 14] | [15, 18] | [3, 8] | [10, 15] | [17, 22] | |
| 6 | [4, 7] | [8, 11] | [12, 15] | [16, 19] | — | — | — | |
| 7 | [3, 6] | [7, 10] | [11, 14] | [15, 18] | [4, 9] | [10, 15] | [16, 21] | |
| 8 | [5, 8] | [9, 12] | [13, 16] | [17, 20] | [5, 10] | [11, 16] | [17, 22] | |
| 9 | [4, 7] | [9, 12] | [14, 17] | [19, 22] | [4, 9] | [11, 16] | [18, 23] | |
| 10 | [5, 8] | [10, 13] | [15, 18] | [20, 23] | [5, 10] | [12, 17] | [19, 24] | |
| 11 | [2, 5] | [7, 10] | [12, 15] | [17, 20] | — | — | — | |
| 12 | [3, 6] | [8, 11] | [13, 16] | [18, 21] | [5, 10] | [12, 17] | [19, 24] | |
| 13 | [5, 8] | [9, 12] | [13, 16] | [17, 20] | [4, 9] | [10, 15] | [16, 21] | |
| 14 | [2, 5] | [6, 9] | [11, 14] | [15, 18] | [5, 10] | [11, 16] | [17, 22] | |
| 15 | [4, 7] | [9, 12] | [14, 17] | [19, 22] | [2, 7] | [8, 13] | [14, 19] | |
| 16 | [2, 5] | [6, 9] | [11, 14] | [15, 18] | [3, 8] | [10, 15] | [17, 22] | |
| 17 | [3, 6] | [8, 11] | [13, 16] | [18, 21] | — | — | — | |
| 18 | — | — | — | — | — | — | — | |
| 19 | — | — | — | — | — | — | — | |
| 20 | [3, 6] | [8, 11] | [13, 16] | [18, 21] | [4, 9] | [11, 16] | [18, 23] | |
| 21 | [2, 5] | [6, 9] | [11, 14] | [15, 18] | [5, 10] | [11, 16] | [17, 22] | |
| 22 | [4, 7] | [8, 11] | [12, 15] | [16, 19] | — | — | — | |
| 23 | [3, 6] | [7, 10] | [11, 14] | [15, 18] | [4, 9] | [10, 15] | [16, 21] | |
| 24 | [5, 8] | [9, 12] | [13, 16] | [17, 20] | [5, 10] | [11, 16] | [17, 22] | |
| 25 | [4, 7] | [9, 12] | [14, 17] | [19, 22] | [2, 7] | [8, 13] | [14, 19] | |
| 26 | [5, 8] | [10, 13] | [15, 18] | [20, 23] | [3, 8] | [10, 15] | [17, 22] | |
| 27 | [2, 5] | [7, 10] | [12, 15] | [17, 20] | [5, 10] | [11, 16] | [17, 22] | |
| 28 | [3, 6] | [8, 11] | [13, 16] | [18, 21] | [2, 7] | [8, 13] | [14, 19] | |
| 29 | [4, 7] | [9, 12] | [14, 17] | [19, 22] | — | — | — | |
| 30 | [2, 5] | [6, 9] | [11, 14] | [15, 18] | [5, 10] | [11, 16] | [17, 22] | |
| 31 | [4, 7] | [9, 12] | [14, 17] | [19, 22] | [2, 7] | [8, 13] | [14, 19] | |
| 32 | [2, 5] | [7, 10] | [12, 15] | [17, 20] | [3, 8] | [10, 15] | [17, 22] | |
| 33 | [3, 6] | [8, 11] | [13, 16] | [18, 21] | — | — | — | |
| 34 | [2, 5] | [6, 9] | [11, 14] | [15, 18] | [2, 7] | [8, 13] | [14, 19] | |
| 35 | [2, 5] | [7, 10] | [12, 15] | [17, 20] | [5, 10] | [11, 16] | [17, 22] | |
| 注:[运输工具到达时间,运输工具离开时间] | ||||||||
| 运输方式 | 公路 | 铁路 | 水路 |
| 公路 | 0/0 | 6.5/0.009 | 7/0.006 |
| 铁路 | 6.5/0.009 | 0/0 | 8/0.012 |
| 水路 | 7/0.006 | 8/0.012 | 0/0 |
参考不同运输方式平均货损情况,假设不同运输方式每百公里的平均货损率分别为:公路0.01%,铁路0.015%,水路0.02%,单次中转平均货损率为0.04%[12]。考虑中国对货损率超过0.5%需要向保险公司索赔,保险理赔时容易导致客户丢失,因此设定货损率限制等级为(L1,L2,L3,L4),分别对应表示为L1(0.5%),L2(0.45%),L3(0.4%)和L4(0.35%),即等级越高,货损率限制越严格。
4.2 结果与分析 4.2.1 求解结果通过MATLAB编程对模型进行求解,算法参数设定为:种群大小300;最大迭代次数300;交叉概率0.8;变异概率0.15;初始温度1 000;降温速率0.98,θ=0.5%,此时对货损率不做限制,满足最低要求即可。目标函数最优解的迭代过程如图 6所示,可知目标函数的最优解逐渐减小至最终收敛,算法达到一定迭代次数后保持数值不变。解得最优路径及运输方式组合如下:1—水路—2—水路—8—铁路—9—铁路—13—铁路—24—铁路—31—铁路—32—铁路—35。路径总时间47.73 h,总成本45 444.57元。
|
| 图 6 求解迭代图 Fig. 6 Iterative diagram of solving |
| |
为满足多式联运经营人多样化的运输需求,输出总成本最小的5种最优运输方案,如表 5所示。由表 5可知:当多式联运经营人要求物流总成本最低时可以选择路径1;当多式联运经营人想要货物在时间窗内尽可能快的运达则可以选择路径2;若多式联运经营人想要规划一条运输质量相对较高的路线,对货损成本有要求,则可以选择路径4。
| 序号 | 路径 | 运输方式 | 总成本/元 | 运输成本/元 | 运输时间/h | 货损成本/元 |
| 1 | 1-2-8-9-13-24-31-32-35 | 3-3-2-2-2-2-2-2 | 45 444.57 | 41 013.00 | 47.73 | 4 161.57 |
| 2 | 1-2-8-9-13-24-31-32-35 | 3-2-2-2-2-2-2-2 | 49 353.06 | 45 001.60 | 42.73 | 4 081.46 |
| 3 | 1-2-8-9-13-24-31-32-35 | 3-3-2-3-2-2-2-2 | 51 950.30 | 46 138.40 | 47.73 | 5 541.90 |
| 4 | 1-4-5-12-16-21-27-28-35 | 2-2-2-2-2-2-2-2 | 56 211.45 | 52 038.00 | 48.69 | 3 902.85 |
| 5 | 1-4-5-12-16-21-27-28-35 | 3-2-2-2-2-2-2-2 | 56 932.13 | 52 047.80 | 62.69 | 4 614.33 |
4.2.2 对比分析
为探究货损率变化对多式联运路径选择的影响,对比有无混合时间窗约束的2种情形下货损限制等级对多式联运路径选择及各项成本的影响。表 6为当不确定需求参数
| 等级 | 混合时间窗 | 路径 | 运输方式 | 总成本/元 | 货损成本/元 | 铁路占比/% |
| L1 | 有 | 1-2-8-9-13-24-31-32-35 | 3-3-2-2-2-2-2-2 | 45 444.57 | 4 161.57 | 76.1 |
| 无 | 1-4-5-12-16-21-27-28-35 | 3-3-3-3-3-3-3-2 | 32 511.01 | 5 737.41 | 24.9 | |
| L2 | 有 | 1-2-8-9-13-24-31-32-35 | 3-2-2-3-3-2-2-2 | 52 534.38 | 5 574.98 | 64.8 |
| 无 | 1-4-5-12-16-21-27-28-35 | 3-3-3-3-3-3-2-2 | 35 840.56 | 5 538.96 | 34.4 | |
| L3 | 有 | 1-4-5-12-16-21-27-28-35 | 3-3-2-2-2-1-3-3 | 64 629.29 | 5 919.69 | 33.2 |
| 无 | 1-4-5-12-16-21-27-28-35 | 3-3-3-3-3-2-2-2 | 39 743.66 | 5 345.66 | 45.2 | |
| L4 | 有 | 1-4-5-12-16-21-27-28-35 | 3-3-3-2-3-3-2-1 | 67 132.71 | 6 734.81 | 22.9 |
| 无 | 1-4-5-12-16-21-27-28-35 | 3-3-3-3-2-2-2-2 | 42 915.67 | 5 146.47 | 54.7 |
(1) 有混合时间窗限制
当货损限制等级为L1,L2时:最优路径为1-2-8-9-13-24-31-32-35,运输方式发生变化,铁路运输在总运输里程中的占比较高,在60%以上。当货损限制等级为L3, L4时:最优路径为1-4-5-12-16-21-27-28-35,此时铁路运输在总运输里程中的占比较低,在40%以下。由此可以发现,当货损限制等级较低,即货损率限制不严格时,多式联运经营人倾向选择运费与货损率都相对较低的铁路运输;而当货损限制等级提升,即货损率限制较严格时,多式联运经营人倾向选择公铁水联运以分担货损率满足要求。也说明了货损限制等级在一定程度上不仅约束了路径的选择,也约束了运输方式的选择。
(2) 无混合时间窗限制
随着货损限制等级的提高,货物的运输路径不发生变化,运输方式发生了改变:路径中只出现了铁水联运且铁路的占比逐渐提升。这是由于本研究是以总物流成本最低为目标函数,因此在无时间窗限制的情况下,多式联运经营人会选择运输时间长但运费低的水路运输。而随着货损限制等级越加严格,货损率高的水路运输便不满足路径的运输要求,此时铁路的占比便逐步提升。由此说明,混合时间窗的限制对运输路径的影响较大,在实际运输时多式联运经营人要重点考虑班列班期。
4.2.3 敏感性分析运输需求量的增减会对运输中的总成本、货损成本等产生影响,从而导致路径的改变。在原需求模糊数基础上分别增减10%,20%,30%,即需求模糊数分别变为(
|
| 图 7 不确定需求量灵敏度分析 Fig. 7 Sensitivity analysis on uncertain demand |
| |
| 需求变化率/% | 路径 | 运输方式 | 总成本/元 | 货损成本/元 | 货损成本变动率/% |
| -10 | 1-2-8-9-13-24-31-32-35 | 2-2-2-2-2-3-3-3 | 36 212.33 | 3 822.81 | 8.86 |
| -20 | 1-2-8-9-13-24-31-32-35 | 3-2-2-2-2-3-3-3 | 33 093.03 | 4 006.63 | 3.87 |
| -30 | 1-2-8-9-13-24-31-32-35 | 3-3-3-2-3-3-2-2 | 30 450.51 | 4 032.65 | 3.20 |
| 0 | 1-2-8-9-13-24-31-32-35 | 3-3-2-2-2-2-2-2 | 45 444.57 | 4 161.57 | 0.00 |
| +10 | 1-2-8-9-13-24-31-32-35 | 3-3-2-2-3-2-2-2 | 57 261.77 | 5 995.84 | 30.59 |
| +20 | 1-2-8-9-13-24-31-32-35 | 3-3-2-2-2-2-3-2 | 63 422.59 | 6 647.63 | 37.40 |
| +30 | 1-4-5-12-16-21-27-28-35 | 3-2-2-2-2-2-2-2 | 75 430.76 | 5 998.63 | 30.62 |
由表 7可以看出,随着需求模糊数降低百分比由10%至30%:物流总成本逐渐降低,此时各节点运输能力满足运输需求,因此运输路径未发生变化,运输方式发生一定的变化,铁路运输在总运输里程中的占比随着货运量的减少由60.37%降低至38.19%,此时货损成本的变动率较低,即当需求模糊数降低时对货损成本影响较小。随着需求模糊数增加百分比由10%至30%:物流总成本逐渐提高,运输路径不变,运输方式变化不显著。当百分比为30%时,运输路径和运输方式都发生了改变。这是由于此时清晰化后的运输需求量为637 t,而原最优路径(1-2-8-9-13-24-31-32-35)中,节点2至节点8水路的运输能力为528 t无法满足运输需求,故会寻找其他中转节点,此时原路径已不满足最优路径的要求,故更换了运输路线,并且货损成本的变动率也显著增加。
5 结论在需求不确定的情况下,以混合时间窗和货损率限制为约束条件,构建以运输成本、时间成本及货损成本最低为目标的多式联运路径选择模型。结果表明:(1)有混合时间窗约束下,当货损限制等级为L1,L2时,铁路在运输方式中占比普遍高于50%;当货损限制等级为L3,L4时,铁路在运输方式中的占比普遍低于50%。即货损率限制不严格时,多式联运经营人倾向选择铁路运输,否则倾向选择公铁水联运以分担货损率。无混合时间窗约束下,随着货损限制等级的提高,货物的运输路径不变,运输方式发生了改变:路径中只出现了铁水联运且铁路的占比逐渐提升。这说明货损限制等级在一定程度上不仅约束了路径的选择也约束了运输方式的选择,而混合时间窗的约束对运输路径的影响较大。(2) 在原需求模糊数的基础上,当运输需求量降低10%~30%,铁路运输在总运输里程中的占比逐渐降低,货损成本变化幅度不显著;当运输需求量逐渐增加10%~30%,铁路运输在总运输里程中的占比逐渐提高。当增加30%时,货损成本急剧增加30%以上。因此,在运输需求量较以往有显著增加的可能性时,多式联运经营人要预估出合理的运输需求范围,在合理预估数据的基础上进行路径规划可以降低成本,提高运输效率。
| [1] |
吕学伟, 杨斌, 黄振东. 混合时间窗约束下多式联运最优路径选择研究[J]. 铁道运输与经济, 2018, 40(8): 6-11. LÜ Xuewei, YANG Bin, HUANG Zhendong. Optimal path selection for multimodal transport under mixed time window constraints[J]. Railway Transport and Economy, 2018, 40(8): 6-11. |
| [2] |
GUO Z X, ZHANG D Q, LIU H T, et al. Green transportation scheduling with pickup time and transport mode selections using a novel multi-objective memetic optimization approach[J].
Transportation Research Part D Transport & Environment, 2018, 60(5): 137-152.
|
| [3] |
温廷新, 李可昕, 胡迎春. 基于软时间窗的多车舱生鲜品配送路径优化[J]. 公路交通科技, 2023, 40(9): 232-238, 256. WEN Tingxin, LI Kexin, HU Yingchun. Optimization of multi-compartment fresh product distribution path based on soft time window[J]. Journal of Highway and Transportation Research and Development, 2023, 40(9): 232-238, 256. DOI:10.3969/j.issn.1002-0268.2023.09.027 |
| [4] |
刘畅, 关秀婷, 张金伟, 等. 考虑时间价值成本的中欧笔记本电脑多式联运路径优化研究[J]. 铁道科学与工程学报, 2019, 16(9): 2352-2359. LIU Chang, GUAN Xiuting, ZHANG Jinwei, et al. Study on the route optimization of China-Europe laptop multimodal transport considering time value cost[J]. Journal of Transportation Engineering and Information, 2019, 16(9): 2352-2359. |
| [5] |
裴骁, 芦有鹏, 张长泽. 基于综合满意度的多式联运路径模型及其算法[J]. 交通信息与安全, 2020, 38(1): 136-144. PEI Xiao, LU Youpeng, ZHANG Changze. A model of multi-modal transport path based on hybrid algorithm considering comprehensive satisfaction[J]. Journal of Transport Information and Safety, 2020, 38(1): 136-144. |
| [6] |
孙岩, 虞楠, 王丹竹, 等. 考虑多类型时间窗的集装箱多式联运路径双目标优化研究[J]. 铁道运输与经济, 2021, 43(10): 82-89. SUN Yan, YU Nan, WANG Danzhu, et al. Study on the bi-objective optimization for the multimodal routing problem with multiple time windows[J]. Railway Transport and Economy, 2021, 43(10): 82-89. |
| [7] |
李春发, 米新新, 崔鑫. 基于双曲正切函数改进蚁群算法的冷链物流配送路径优化[J]. 公路交通科技, 2023, 40(12): 236-244, 258. LI Chunfa, MI Xinxin, CUI Xin. Optimization of cold chain logistics distribution routing based on improved ant colony algorithm with hyperbolic tangent function[J]. Journal of Highway and Transportation Research and Development, 2023, 40(12): 236-244, 258. DOI:10.3969/j.issn.1002-0268.2023.12.027 |
| [8] |
SEO Y J, CHEN F L, ROH S Y. Multimodal transportation: The case of laptop from Chongqing in China to Rotterdam in Europe[J].
The Asian Journal of Shipping and Logistics, 2017, 33(3): 155-165.
DOI:10.1016/j.ajsl.2017.09.005 |
| [9] |
DULEBENETS M A, OZGUVEN E E, MOSES R, et al. Intermodal freight network design for transport of perishable products[J].
Open Journal of Optimization, 2016, 5(4): 120-139.
DOI:10.4236/ojop.2016.54013 |
| [10] |
鲁玉, 徐行方, 尹传忠, 等. 铁路冷链物流运输多目标机会约束规划[J]. 同济大学学报(自然科学版), 2021, 49(10): 1407-1416. LU Yu, XU Xingfang, YIN Chuanzhong, et al. Multi-objective chance-constrained programming of railway cold chain logistics[J]. Journal of Tongji University(Natural Science), 2021, 49(10): 1407-1416. |
| [11] |
冯芬玲, 孙楠佳. 考虑时间价值的中非多式联运路径与出海港选择[J]. 交通运输系统工程与信息, 2022, 22(2): 45-53. FENG Fenling, SUN Nanjia. Multimodal transportation route and port selection between China and Africa considering time value[J]. Journal of Transportation Systems Engineering and Information Technology, 2022, 22(2): 45-53. |
| [12] |
陈维亚, 龚浩, 方晓平. 考虑运输碳税与质量承诺的多式联运路径优化[J]. 铁道科学与工程学报, 2021, 19(1): 34-41. CHEN Weiya, GONG Hao, FANG Xiaoping. Multimodal transportation route optimization considering transportation carbon tax and quality commitment[J]. Journal of Railway Science and Engineering, 2021, 19(1): 34-41. |
| [13] |
陈汨梨, 赵孝进, 邓夕贵, 等. 不确定条件下的多式联运路径优化[J]. 公路交通科技, 2021, 38(1): 143-150, 158. CHEN Mili, ZHAO Xiaojin, DENG Xigui, et al. Multimodal transport path optimization under uncertain conditions[J]. Journal of Highway and Transportation Research and Development, 2021, 38(1): 143-150, 158. DOI:10.3969/j.issn.1002-0268.2021.01.018 |
| [14] |
IWAMURA K, LIU B D. A genetic algorithm for chance constrained programming[J].
Journal of Information and Optimization Sciences, 1996, 17(2): 409-422.
DOI:10.1080/02522667.1996.10699291 |
| [15] |
XIONG G W, WANG Y. Best routes selection in multimodal networks using multi-objective genetic algorithm[J].
Journal of Combinatorial Optimization, 2014, 28(3): 655-673.
DOI:10.1007/s10878-012-9574-8 |
| [16] |
SUN Y, LANG M X. Bi-objective optimization for multi-modal transportation routing planning problem based on pareto optimality[J].
Journal of Industrial Engineering and Management, 2015, 8(4): 1195-1217.
|
| [17] |
KOZA D F. Liner shipping service scheduling and cargo allocation[J].
European Journal of Operational Research, 2019, 275(3): 897-915.
DOI:10.1016/j.ejor.2018.12.011 |
| [18] |
张旭, 张海燕, 袁旭梅, 等. 双重不确定下低碳多式联运路径优化研究[J]. 北京交通大学学报(社会科学版), 2022, 21(2): 113-121. ZHANG Xu, ZHANG Haiyan, YUAN Xumei, et al. Path optimization of low-carbon multimodal transportation under dual uncertainty[J]. Journal of Beijing Jiaotong University(Social Sciences Edition), 2022, 21(2): 113-121. |
2025, Vol. 42


,