森林与环境学报  2021, Vol. 41 Issue (1): 88-95   PDF    
http://dx.doi.org/10.13324/j.cnki.jfcf.2021.01.012
0

文章信息

陈诚, 林秋婷, 邱荣祖
CHEN Cheng, LIN Qiuting, QIU Rongzu
基于随机规划模型的弹性木材供应链网络优化
Cost optimization of resilient timber supply chain networks based on stochastic programming
森林与环境学报,2021, 41(1): 88-95.
Journal of Forest and Environment,2021, 41(1): 88-95.
http://dx.doi.org/10.13324/j.cnki.jfcf.2021.01.012

文章历史

收稿日期: 2020-10-29
修回日期: 2020-11-30
基于随机规划模型的弹性木材供应链网络优化
陈诚 , 林秋婷 , 邱荣祖     
福建农林大学交通与土木工程学院, 福建 福州 350002
摘要:为提高木材供应链网络的弹性和可靠性,揭示随机事件发生概率对木材供应链网络弹性的影响,针对木材供应链网络特点,综合运用概率分析和情景树分析法,对可能发生不确定性事件的木材供应链网络的优化问题进行研究。构建了考虑节点失效的弹性供应链网络随机规划模型,通过模型求解得出最优网络方案;引入供应链网络弹性评价指标对求解得到的供应链网络的弹性进行测度量化,并与确定性优化模型的结果进行比较;进行了随机事件发生概率的敏感性分析。基于实例的数值计算结果表明,在不发生随机事件的情景下,提出的模型所得解的网络总成本比确定性模型的解高出1.75%,而在其他情景下,前者的网络总成本均小于后者,就网络总成本的期望值而言,提出的模型比确定性模型低19.92%;同时,供材点失效概率及对应的供材量与供应链网络期望总成本成正比。因此,提出的随机规划模型是科学有效的,能作为弹性木材供应链网络优化的决策工具;构建木材供应链网络时,应避免所有供材点均受到同一随机事件影响,并尽可能降低供应量大的供材点发生失效的概率。
关键词木材供应链网络    中断风险    供应失效    随机规划    
Cost optimization of resilient timber supply chain networks based on stochastic programming
CHEN Cheng , LIN Qiuting , QIU Rongzu     
College of Transportation and Civil Engineering, Fujian Agriculture and Forestry University, Fuzhou, Fujian 350002, China
Abstract: To increase the reliability of timber supply chain networks, we examined the probability of stochastic events on their resilience and where these events cause the failure of supply points. Considering the characteristics of the timber supply chain network, a stochastic programming model was proposed to address the scenario of node failure based on the results of a probability analysis and scenario tree. The results were compared with those of a deterministic model, and the resilience metrics of both networks were measured by the introduced resilience assessment indicator for the supply chain network. Finally, sensitivity analysis was conducted on the probabilities of the stochastic events and the expected total cost of the supply chain network. It was shown that although the proposed model resulted in a total cost that was 1.75% greater than that of the deterministic model with no stochastic events, the expected total costs derived by the former model were 19.92% lower. Furthermore, the probability of stochastic events causing failure of the supply points, and the scale of the corresponding supply volume, are directly proportional to the expected total cost. The proposed model was proven to be valid and effective; therefore, it can be used as a decision-making tool to ensure resilient timber supply chains. Meanwhile, it was shown that all supply points are also subject to the same stochastic events that should be avoided; further, the failure probability of supply points that provide the largest timber volumes should be limited when configuring a timber supply chain.
Key words: timber supply chain network     disruption risk     supply failure     stochastic programming    

木材作为重要的建筑原材料之一,生长周期较长,其供应量及供应可靠性易受到诸多因素的影响,如自然环境、国内政策、国际政策等,导致其供应具有不确定性。一些动态网络优化模型通过考虑网络参数的波动,改善了传统静态优化模型对供应链网络不确定性考虑不足的局限性[1],但缺乏对供应链网络受到较大影响甚至发生中断的情形的刻画[2],之后学者们开始了对弹性网络的思考和研究。HUTCHISON et al[3]将弹性定义为网络在中断期间维持其服务水平的能力,一些学者研究了有助于供应链中断后快速恢复的应急策略[4-6],而另一些学者则研究了考虑不确定性和中断风险的供应链网络优化模型[7-9]

由于木材供应链网络具有较经典供应链网络更为复杂的跨层级结构[10-11],其优化设计问题也受到了学者们的持续关注,但是已有研究大多针对确定性优化问题。例如,TRONCOSO et al[12]研究了包含木材砍伐、设施选址以及木材运输三类决策的确定性木材供应链网络优化模型。KANZIAN et al[13]用两个接续的运输模型描述了一个区域木材供应链网络优化问题。ABASIAN et al[14]通过纳入设施的技术选择决策进一步丰富了木材供应链网络优化模型。CHEN et al[15]通过将规划期划分为多个周期刻画了供应量的动态性,同时用三角模糊数描述了木材需求的不确定性,建立了木材供应链网络优化的动态模型,但没有考虑到木材供应的不确定性。

木材具有长、大、笨重的特点,运输量一般较大,即需求满足的来源较集中,这能够有效节约运输成本,但通常具有较低的弹性[16]。随着我国木材需求量的不断提高,弹性木材供应链网络的构建成为了亟待解决的问题之一。然而,考虑供应失效的木材供应链网络优化研究仍十分鲜见。鉴于此,本研究针对木材供应链网络结构特点,利用随机规划方法,结合情景树分析构建弹性木材供应链网络优化模型,分析木材供应可靠性影响因素对供应链网络的影响,并通过比较确定性模型和弹性模型得出的最优木材供应链网络方案揭示从网络设计角度增加木材供应链网络弹性的途径,为弹性木材供应链网络设计决策提供工具和理论支持。

1 供应链网络弹性的评价指标

供应链网络特征表现出与网络弹性较强的相关性,DIXIT et al[16]基于表征弹性网络行为的网络设计指标,如节点密度、网络连接度、网络规模等[17],提出了供应链网络弹性指标,该指标由4个网络特征指标组成,如式(1)。

$ R = \frac{{{C_{\rm{V}}}{N_{\rm{S}}}}}{{D{C_{\rm{T}}}}} $ (1)

式中:R为网络弹性;CV为网络连接度,用所有供应点到所有需求点间的路径条数表示,网络连接度越高,供应链网络在中断后能保持的服务水平越高[18],故其与供应链网络的弹性成正比;NS为网络规模,用网络中的节点个数和弧的条数之和表示,网络规模越大,在供应链中断后,采用额外的资源和路径充当缓冲的可能性越大,故网络规模越大,弹性越大[19]D为网络密度,定义为单位距离上节点的数量,网络密度越大,则局部中断事件对网络的影响越大,在中断时表现出的脆弱性更大[20]CT为网络中心性指标,用平均节点度表示,供应点中心性较低的供应链网络具有更大的弹性[21]

2 考虑节点失效的弹性木材供应链网络优化模型 2.1 问题描述

研究对象为由供材点、木材物流中心以及需材点构成的单产品三级木材供应链网络,其中需材点的需求既可直接从供材点(伐区的砍伐点、港口、车站等)获得,也可以从木材物流中心获得。各供材点的供应量可通过科学预测获得,但其供应可靠性受到一些因素的影响,即供应点具有一定的失效风险。为了更好地贴近实际并增加模型的灵活性,在优化模型中引入弹性供应量,即当节点发生失效时,还可能通过使用额外的弹性供应成本提供少量的供应。同时,由于供材点失效后,需材点无法马上与新的供材点建立连接,于是木材流量只能在已建立连接的节点间发生。需材点为各类木材加工厂,如人造板厂、家具厂,应尽可能满足需材点的需求;同时在目标函数中引入缺货成本,在促进需求量被满足的同时允许存在一部分不被满足的需求量。需要从有容量限制的候选木材物流中心和供材点中进行选择,与需求点共同组成供应链网络中的节点,并进行节点间路径连接决策,形成最终的供应链网络,目标是使网络的期望总成本最低。

基于进口和本地两种木材来源,考虑4类可能影响木材供应可靠性的供材点失效事件,如表 1所示。综合运用概率分析和情景树分析法,可以得出供材点失效情景及其概率,图 1给出了一个分析示例。假设有4个供材点(M1、M2、M3、M4),其中M1为港口,M2、M3、M4为人工用材林,因国内林业政策影响供应的概率为Pd,国际形势不出现影响供应的概率为Pg,不发生森林灾害的概率分别为Pf, 2Pf, 3Pf, 4。M2、M3位于区域Ⅰ,M4位于区域Ⅱ,区域Ⅰ和区域Ⅱ不出现恶劣气候的概率分别为Pc, ⅠPc, Ⅱ。因此,M1、M2、M3、M4的失效概率分别为1-Pg、1-PdPc, ⅠPf, 2、1-PdPc, ⅠPf, 3和1-PdPc, ⅡPf, 4;情景s1~s16的发生概率分别为情景树从M1至M4的各条分支上的概率的乘积。

图 1 具有4个供材点的木材供应链网络情景树示意图 Fig. 1 Illustration of scenarios tree for the timber supply chain network with 4 supply points
表 1 供材点失效事件 Table 1 Events causing failure of supply points
编号No. 失效事件Event 影响对象Influence object 范围Range 发生概率Probability
1 恶劣气候Harsh climate 区域内的供材点Supply points in the same region 中Medium 1-Pc, Ⅰ
2 森林灾害Forest disaster 个别供材点Single supply point 小Small 1-Pf, j
3 国内林业政策Domestic forestry policy 国内供材点Domestic supply points 大Large 1-Pd
4 国际形势International situation 进口木材供应点(港口、车站) Supply points for import timber(ports and stations) 中Medium 1-Pg
2.2 随机规划模型

考虑节点失效的弹性木材供应链网络优化问题可以用以下的随机规划模型描述。

集合:I为供材点集合;J为木材物流中心集合;K为需材点的集合;S为情景集合。

参数:Ps为情景s发生的概率;cij为木材供应链网络中节点i与节点j间的单位运输成本;fj为木材物流中心j的固定成本;γ为单位缺货成本;Dk为需材点k对木材的需求量;ec, i为供材点i的单位弹性供应成本;Ri为供材点i的弹性供应量;Mcap, i为供材点i的正常供应量;Dcap, j为木材物流中心j最大的处理能力;E1E2E3分别表示供材点与物流中心,物流中心与需材点、以及供材点与需材点之间的最小起运量;Ys, i表示情景s下供材点i是否失效,等于1表示失效,否则等于0。

决策变量:Qs, ij为情景s下供材点i到木材物流中心j的正常运输量;Qs, jk为情景s下木材物流中心j到需材点k的正常运输量;Qs, ik为情景s下供材点i到需材点k的正常运输量;Qs*, ij为情景s下供材点i到木材物流中心j的弹性运输量;Qs*, ik为情景s下供材点i到需材点k的弹性运输量;Xj为木材物流中心j是否被选择的二维变量,等于1表示被选中,否则等于0;ls, ij为情景s下节点i和节点j之间是否存在连接,等于1表示存在,否则等于0。

目标函数:

$ \min = \sum\limits_{s \in S} {{P_s}} \left\{ {\begin{array}{*{20}{l}} {\sum\limits_{i \in I} {\sum\limits_{j \in J} {{c_{ij}}} } \left( {{Q_{s,ij}} + {Q_{s * ,ij}}} \right) + \sum\limits_{j \in J} {\sum\limits_{k \in K} {{c_{jk}}} } {Q_{s,jk}} + \sum\limits_{i \notin I} {\sum\limits_{k \in K} {{c_{ik}}} } \left( {{Q_{s,ik}} + {Q_{s * ,ik}}} \right)}\\ { + \gamma \left( {\sum\limits_{k \in K} {{D_k}} - \sum\limits_{j \in J} {\sum\limits_{k \in K} {{Q_{s,jk}}} } - \sum\limits_{i \in I} {\sum\limits_{k \in K} {\left( {{Q_{s,ik}} + {Q_{s * ,ik}}} \right)} } + \sum\limits_{i \in I} {{e_{c,i}}} \left( {\sum\limits_{j \in J} {{Q_{s * ,ij}}} + \sum\limits_{k \in K} {{Q_{s * ,ik}}} } \right)} \right.} \end{array}} \right\} + \sum\limits_{j \in J} {{f_j}} {X_j} $ (2)

Subject to

$ {\sum\limits_{j \in J} {{Q_{s,ij}}} + \sum\limits_{k \in K} {{Q_{s,ik}}} \le {M_{cap,i}}\left( {1 - {Y_{s,i}}} \right),\forall s \in S,i \in I} $ (3)
$ {\sum\limits_{j \in J} {{Q_{s * ,ij}}} + \sum\limits_{k \in K} {{Q_{s * ,k}}} \le {R_i},\forall s \in S,i \in I} $ (4)
$ {\sum\limits_{i \in I} {\left( {{Q_{s,ij}} + {Q_{s * ,ij}}} \right)} \le {D_{cap,j}}{X_j},\forall s \in S,j \in J} $ (5)
$ \sum\limits_{i \in I} {\left( {{Q_{s,ij}} + {Q_{s * ,ij}}} \right)} \ge \sum\limits_{k \in K} {{Q_{s,jk}}} ,\forall s \in S,j \in J $ (6)
$ \sum\limits_{i \in I} {\left( {{Q_{s,ik}} + {Q_{s * ,ik}}} \right)} + \sum\limits_{j \in J} {{Q_{s,jk}}} \le {D_k},\forall s \in S,k \in K $ (7)
$ {{Q_{s,ij}} + {Q_{s * ,ij}} \ge {E_1}{l_{s,ij}},\forall s \in S,i \in I,j \in J} $ (8)
$ {{Q_{s,jk}} \ge {E_2}{l_{s,jk}},\forall s \in S,j \in J,k \in K} $ (9)
$ {{Q_{s,ik}} + {Q_{s * ,k}} \ge {E_3}{l_{s,ik}},\forall s \in S,i \in I,k \in K} $ (10)
$ {{Q_{s,ik}} + {Q_{s * ,jk}} \le {l_{s,ik}}{M_{cap,i}},\forall s \in S,i \in I,k \in K} $ (11)
$ {{Q_{s,ij}} + {Q_{s * ,ij}} \le {l_{s,ij}}{M_{cap,i}},\forall s \in S,i \in I,j \in J} $ (12)
$ {{Q_{s,jk}} + {Q_{s * ,jk}} \le {l_{s,jk}}{D_k},\forall s \in S,j \in J,k \in K} $ (13)
$ {{l_{s,ij}} \le {l_{1,ij}},\forall s \in S,i \in I,j \in J} $ (14)
$ {{l_{s,jk}} \le {l_{1,jk}},\forall s \in S,j \in J,k \in K} $ (15)
$ {{l_{s,ik}} \le {l_{1,ik}},\forall s \in S,i \in I,k \in K} $ (16)
$ {{Q_{s,ij}},{Q_{s,jk}},{Q_{s,ik}},{Q_{s * ,ij}},{Q_{s * ,ik}} \ge 0,\forall i \in I,j \in J,k \in K,s \in S} $ (17)
$ {{X_j},{l_{ik}} \in \{ 0,1\} ,\forall i \in I,j \in J,k \in K} $ (18)

式中:目标函数(2)最小化供应链网络总成本,包括3个层级节点间的运输成本、供材点的弹性供应成本、木材物流中心的固定成本、需材点的缺货成本。式(3)为不失效情况下供材点的供应量约束;式(4)为失效情景供材点的弹性供应量约束;式(5)为木材物流中心的能力约束;式(6)保证流出木材物流中心的木材量不超过流出量;式(7)保证需材点的需求不会超量满足,从而保证非负的缺货成本;式(8)~(10)为3类连接的最低起运量约束;式(11)~(12)为二维变量lik的逻辑约束的线性化表达,即若节点i和节点k之间不存在连接(lik=0),则ik间的木材流量为零;式(13)~(16)表明突发事件发生时无法立刻建立新连接,只能使用正常状态(无供应点失效)下的既有连接;式(17)为连续决策变量的非负约束;式(18)约束了二维决策变量。

3 实例分析 3.1 案例介绍

以某木材供应链网络为背景,该木材供应链网络包含4个供材点(M1~M4)和5个需材点(P1~P5)。4个供材点由1个港口(记为M1)和3个人工用材林(M2、M3、M4)组成,其中M2、M3同处一区域,M4位于另一区域。基于对木材物流中心的选址需求分析,利用ArcGIS软件找到4个木材物流中心(L1~L4)的候选址。案例相关数据如表 2表 3所示,供材点与木材物流中心和需材点间的最低起运量为20 000 m3,木材物流中心与需求点间的最低起运量为1 000 m3。按照最大供应量的20%设置供材点的弹性供应能力,弹性供应成本和缺货成本分别为20和50元·m-3

表 2 供材点与其他节点的单位运输成本及其供应量 Table 2 Unit transportation costs related to supply points and the corresponding volume
供材点
Supply point
单位运输成本Unit transportation cost/(yuan·m-3) 供应量
Supply volume/m3
木材物流中心Timber logistics center 需材点Demand point
L1 L2 L3 L4 P1 P2 P3 P4 P5
M1 2.0 3.0 4.0 3.0 2.0 3.0 5.0 8.0 10.0 120 000
M2 4.0 4.0 4.2 3.0 3.0 4.0 5.0 7.0 8.0 110 000
M3 3.0 6.0 5.0 4.0 4.0 5.0 3.0 10.0 12.0 180 000
M4 3.5 4.0 4.0 2.8 2.0 2.0 3.0 6.0 7.0 110 000
表 3 需材点与物流中心间的单位运输成本及其需求量 Table 3 Unit transportation costs related to demand points and the corresponding volume
需材点
Demand point
单位运输成本Unit transportation cost/(yuan·m-3) 需求量
Demand volume/m3
L1 L2 L3 L4
P1 4.0 3.0 3.8 2.2 60 000
P2 4.7 3.0 3.5 2.0 120 000
P3 5.3 3.5 3.8 2.2 140 000
P4 6.0 4.0 4.0 5.0 80 000
P5 7.0 4.5 4.0 5.5 90 000

根据历史数据及现状分析,预测国际形势发生变化的概率为30.00%,国内政策发生变化的概率为10.00%,区域Ⅰ和区域Ⅱ出现恶劣气候发生的概率分别10.00%和5.00%。最后,根据参考文献[22]预测3处人工用材林的森林灾害发生概率分别为1.00%、2.00%和2.00%。基于前述分析,可得16种情景及其概率,如表 4所示,其中“√”表示对应的供材点失效。

表 4 失效情景及其概率 Table 4 Scenarios and the corresponding probabilities
供材点
Supply point
情景Scenario
s1 s2 s3 s4 s5 s6 s7 s8 s9 s10 s11 s12 s13 s14 s15 s16
M1
M2
M3
M4
概率
bProbability/%
37.34 7.22 9.70 1.88 9.22 1.78 2.40 0.46 16.00 3.10 4.16 0.80 3.95 0.76 1.03 0.20
3.2 案例求解及分析

利用IBM ILOG CPLEX 12.9优化求解器对2.2中的随机规划模型进行编译求解,得到的供应链网络的期望总成本为4 826 230.29元,选择的木材物流中心为{L4}。若不考虑失效,即情景s1的发生概率为100%,此时得到确定性优化模型解方案,总成本为2 054 000元,选择的木材物流中心为{L2}。表 5比较了16种情景下的确定性优化模型(deterministic optimization model, DOM)解与随机规划模型(stochastic programming model, SPM)解的成本。

表 5 不同情景下两种模型的成本对比 Table 5 Comparison of costs between the two models
情景
Scenario
运输成本
Transport cost /(104 yuan)
物流中心固定成本
Fixed cost of logistics center/(104 yuan)
弹性供应成本
Resilient supply cost /(104 yuan)
缺货成本
Stock-out cost /(104 yuan)
总成本
Total cost/(104 yuan)
偏差
Deviation /%
确定性优化模型DOM 随机规划模型SPM 确定性优化模型DOM 随机规划模型SPM 确定性优化模型DOM 随机规划模型SPM 确定性优化模型DOM 随机规划模型SPM 确定性优化模型DOM 随机规划模型SPM
s1 199.40 204.50 6.00 4.50 0.00 0.00 0.00 0.00 205.40 209.00 1.75
s2 179.52 217.00 6.00 4.50 208.00 160.00 230.00 0.00 623.52 381.50 -38.82
s3 161.00 175.40 6.00 4.50 44.00 208.00 590.00 230.00 801.00 617.90 -22.86
s4 146.52 117.42 6.00 4.50 208.00 208.00 780.00 780.00 1 140.52 1 109.92 -2.68
s5 146.44 198.27 6.00 4.50 136.00 160.00 360.00 0.00 648.44 362.77 -44.06
s6 99.40 155.76 6.00 4.50 208.00 208.00 780.00 430.00 1 093.40 798.26 -26.99
s7 113.44 104.52 6.00 4.50 136.00 208.00 910.00 780.00 1 165.44 1 097.02 -5.87
s8 66.40 61.81 6.00 4.50 208.00 208.00 1 330.00 1 330.00 1 610.40 1 604.31 -0.38
s9 163.92 209.20 6.00 4.50 136.00 180.00 410.00 0.00 715.92 393.70 -45.01
s10 147.12 166.48 6.00 4.50 208.00 208.00 830.00 480.00 1 191.12 858.98 -27.88
s11 114.12 101.12 6.00 4.50 208.00 208.00 1 380.00 830.00 1 708.12 1 143.62 -33.05
s12 114.12 76.16 6.00 4.50 208.00 208.00 1 380.00 1 380.00 1 708.12 1 668.66 -2.31
s13 83.80 139.80 6.00 4.50 136.00 208.00 960.00 480.00 1 185.80 832.30 -29.81
s14 67.00 99.24 6.00 4.50 208.00 208.00 1 380.00 1 030.00 1 661.00 1 341.74 -19.22
s15 50.80 49.10 6.00 4.50 136.00 208.00 1 510.00 1 380.00 1 702.80 1 641.60 -3.59
s16 34.00 26.60 6.00 4.50 208.00 208.00 1 930.00 1 930.00 2 178.00 2 169.10 -0.41
期望值
Expected value
164.94 185.38 90.97 117.97 391.15 174.77 653.05 482.62 -19.92

表 5可以看出,仅在不发生供应点失效的情景下随机规划模型解产生了比确定性优化模型解多1.75%的总成本,其余15种情景下随机规划模型解的供应链网络总成本均不同程度小于确定性优化模型解。其中,最大的节约发生在情景s9(概率16.00%),总成本低45.01%;最小的发生在情景s8(概率0.46%)。此外,由于设置了较高的单位缺货成本,故随机规划模型解具有更高的期望运输成本、更高的期望弹性供应成本和更低的期望缺货成本。从对需求点的供应来看,随机规划模型的解更大程度地满足了需材点的需求。如前所述,若供应链网络中的两个节点间没有建立连接,则当节点发生失效时无法及时建立新连接,从而可能导致较大的缺货量。显然,随机模型的解在网络节点间建立了供应点和需求点之间更多的连接,虽然这导致在无失效事件发生的情景(s1)下的总成本略高于确定性优化模型的解,但却具有更高的弹性以应对节点失效的风险。

表 6给出了两种模型解的供应链网络弹性评价指标值。本研究提出的弹性木材供应链网络随机规划模型能给出弹性更高的供应链网络方案。此外,随机规划模型得到的木材供应链网络连接度(CV=20)显著大于确定性优化模型(CV=8)。由于本案例中总供应量与总需求量基本相当,故建立的供应链网络均包含了所有的供材点;同时两种方案均选择了一个木材物流中心,故总节点数相等,因此两种模型解的NS指标值的差异全部体现在网络中弧的条数的差异上。随机规划模型得到的木材供应链网络中的弧的条数比确定性优化模型得到的木材供应链网络中多了42.11%,然而由于总节点数相同,导致了更高的节点平均度。可以推断,若存在更多可供选择的供应点能提供更多的供应量,随机规划模型的解将比确定性优化模型的解包含更多的节点,从而得到R值更高的木材供应链网络。

表 6 两种模型解的木材供应链网络弹性评价指标值 Table 6 Indicator values of timber supply chain network resilience measurements for solutions of the two models
模型类型Model type 网络连接度CV 网络规模NS 网络密度D 网络中心性指标CT 网络弹性R
确定性优化模型DOM 8 19 2.61 1.8 32.37
随机规划模型SPM 20 27 2.85 3.4 55.78
3.3 节点失效事件概率敏感性分析

对4类失效事件的发生概率进行敏感性分析,考虑各类事件的发生概率在[-30%,30%]范围内变化时,供应链期望总成本的变化率,结果如图 2所示。供应链网络期望总成本表现出对因国内政策导致供应点失效的发生概率的变化最为敏感,而对因森林灾害导致的供应点失效的发生率的变化最不敏感。这是因为期望网络总成本的大小和节点失效概率及失效影响的供应量的减少相关。虽然影响进口木材供应的国际形势变化而引起供应点失效的发生概率最大(30.00%),但该概率仅对M1的供应产生影响,受影响的供应量占总供应量的23.1%;而因国内政策导致供材点失效的发生概率虽仅为10.00%,但其影响面较大,M2、M3和M4均受到该因素的影响,受影响的供应量占总供应量的77.9%;其次,虽然也同时受到因森林灾害导致的供材点失效的影响,但由于其发生概率较低(≤0.02%),故对弹性供应链网络优化结果的影响最小。

图 2 供材点失效原因发生概率敏感性分析 Fig. 2 Sensitivity analysis of the probability of failure of material supply points
4 结论

研究针对木材供应链网络特点,分析了供材点的失效事件,并基于概率分析和情景树分析法得出了供材点的失效情景及其概率。基于此,建立了弹性木材供应链网络优化的随机规划模型。为合理评价该模型,引入了网络弹性评价指标,对模型得出的网络方案进行弹性评价。其次,通过实例分析验证了该模型的科学有效性,同时也表明,就网络结构而言,在供材点和需材点间建立更多的路径连接能有效提高供应链网络的弹性。最后,对失效事件发生概率的敏感性分析表明节点失效概率及相应的供应量与供应链网络期望总成本成正比,故为了构建更具弹性的网络,应选择失效概率较低的供应点,并尽量不同时选择受到同一随机事件影响的供材点。

参考文献(References)
[1]
DEHGHAN E, NIKABADI M S, AMIRI M, et al. Hybrid robust, stochastic and possibilistic programming for closed-loop supply chain network design[J]. Computers & Industrial Engineering, 2018, 123: 220-231.
[2]
KIM Y, CHEN Y S, LINDERMAN K. Supply network disruption and resilience:a network structural perspective[J]. Journal of Operations Management, 2015, 33-34: 43-59. DOI:10.1016/j.jom.2014.10.006
[3]
HUTCHISON D, STERBENZ J P G. Architecture and design for resilient networked systems[J]. Computer Communications, 2018, 131: 13-21. DOI:10.1016/j.comcom.2018.07.028
[4]
张松涛, 张盼盼, 张敏, 等. 随机供应中断下含提前期的供应链应急策略[J]. 计算机集成制造系统, 2019, 25(3): 772-779.
[5]
CHEN H Y, DAS A, IVANOV D. Building resilience and managing post-disruption supply chain recovery:lessons from the information and communication technology industry[J]. International Journal of Information Management, 2019, 49: 330-342. DOI:10.1016/j.ijinfomgt.2019.06.002
[6]
刘明菲, 李鑫. 农产品供应链中断情景下的响应策略对供应链弹性的影响[J]. 武汉理工大学学报(社会科学版), 2018, 31(5): 137-143.
[7]
肖建华, 刘侠, 尚帅, 等. 基于节点失效和需求不确定的弹性供应链网络优化模型与算法[J]. 统计与决策, 2018, 34(17): 50-53.
[8]
王海燕, 王云瑜. 灰需求下考虑缺货和中断风险的供应链网络设计[J]. 福州大学学报(哲学社会科学版), 2019, 33(2): 38-44.
[9]
BOTTANI E, MURINO T, SCHIAVO M, et al. Resilient food supply chain design:modelling framework and metaheuristic solution approach[J]. Computers & Industrial Engineering, 2019, 135: 177-198.
[10]
SARRAZIN F, LEBEL L, LEHOUX N. Identifying key factors for the success of a regional logistic center[J]. Forest Science, 2018, 64(3): 233-245. DOI:10.1093/forsci/fxy001
[11]
MALLADI K T, QUIRION-BLAIS O, SOWLATI T. Development of a decision support tool for optimizing the short-term logistics of forest-based biomass[J]. Applied Energy, 2018, 216: 662-677. DOI:10.1016/j.apenergy.2018.02.027
[12]
TRONCOSO J J, GARRIDO R A. Forestry production and logistics planning:an analysis using mixed-integer programming[J]. Forest Policy and Economics, 2005, 7(4): 625-633. DOI:10.1016/j.forpol.2003.12.002
[13]
KANZIAN C, HOLZLEITNER F, STAMPFER K, et al. Regional energy wood logistics-optimizing local fuel supply[J]. Silva Fennica, 2009, 43(1): 113-128.
[14]
ABASIAN F, RÖNNQVIST M, OUHIMMOU M. Forest fibre network design with multiple assortments:a case study in Newfoundland[J]. Canadian Journal of Forest Research, 2017, 47(9): 1232-1243. DOI:10.1139/cjfr-2016-0504
[15]
CHEN C, GAN J B, ZHANG Z X, et al. Multi-objective and multi-period optimization of a regional timber supply network with uncertainty[J]. Canadian Journal of Forest Research, 2020, 50(2): 203-214. DOI:10.1139/cjfr-2018-0551
[16]
DIXIT V, VERMA P, TIWARI M K. Assessment of pre and post-disaster supply chain resilience based on network structural parameters with CVaR as a risk measure[J]. International Journal of Production Economics, 2020, 227: 107655. DOI:10.1016/j.ijpe.2020.107655
[17]
CARDOSO S R, PAULA BARBOSA-PÓVOA A, RELVAS S, et al. Resilience metrics in the assessment of complex supply-chains performance operating under demand uncertainty[J]. Omega, 2015, 56: 53-73. DOI:10.1016/j.omega.2015.03.008
[18]
ADENSO-DÍAZ B, MAR-ORTIZ J, LOZANO S. Assessing supply chain robustness to links failure[J]. International Journal of Production Research, 2018, 56(15): 5104-5117. DOI:10.1080/00207543.2017.1419582
[19]
CRAIGHEAD C W, BLACKHURST J, RUNGTUSANATHAM M J, et al. The severity of supply chain disruptions:design characteristics and mitigation capabilities[J]. Decision Sciences, 2007, 38(1): 131-156. DOI:10.1111/j.1540-5915.2007.00151.x
[20]
NAIR A, VIDAL J M. Supply network topology and robustness against disruptions-an investigation using multi-agent model[J]. International Journal of Production Research, 2011, 49(5): 1391-1404. DOI:10.1080/00207543.2010.518744
[21]
CHOPRA S, SODHI M S. Reducing the risk of supply chain disruptions[J]. MIT Sloan Management Review, 2014, 55(3): 73-80.
[22]
陈钦, 刘伟平. 福建省人工用材林收益与风险分析[J]. 林业科学, 2006, 42(2): 93-97.