1. 引言
为了应对全球气候变暖并推动经济高质量发展,中国在2020年9月第七十五届联合国大会上明确提出了“中国二氧化碳力争在2030年前达到峰值,努力争取在2060年前实现碳中和”的目标[1]。双碳背景下,汽车工业作为全球经济的重要组成部分之一,其对环境的影响日益凸显。汽车制造过程中的能源消耗、废弃物排放和运输等环节都对环境产生负面影响。特别是温室气体的排放问题,因此,实现汽车行业的可持续发展已成为各国政府和企业共同关注的焦点。
汽车供应链管理涉及到众多环节,包括供应商管理、生产计划、库存管理、物流运输等。在这个复杂的供应链网络中,汽车供应链管理需要考虑多个目标,如降低成本、减少碳排放、提高产品质量、提高供应链韧性等[2]。同样,不确定性因素非常多样化,可能来自市场需求波动、原材料供应中断、交通运输延误、自然灾害等方面[3]。这些不确定性因素会对供应链的稳定性和效率造成影响,进而影响到企业的经济效益和可持续发展目标的实现[4]。受近年的疫情影响,我国各产业供应链的稳定性均受到了一定的影响[5],对我国的经济韧性与中国高质量发展产生了一定负面作用[6],汽车供应链也不例外[7]。
本文研究的重点是考虑经济、可持续和韧性目标的三重权衡的汽车供应链。因此,该问题包括多目标规划,该模型根据专家建议给出折中的权重,结合所提出的交互式求解方案对该问题进行探究。
2. 多目标汽车供应链问题描述
该问题的结构是一个汽车供应链,由一组制造设施(i)从一组供应商(s)购买各种原材料,再制造产品并分销到各个销售区域(j)。需要开设收集中心(k)来收集和处理使用过的产品。然后,用过的产品被送到二手市场(l)出售给需要的用户、另一部分运输至回收厂(m),以回收使用后的旧产品。供应链节点索引见表1,决策变量见表2,参数及释义见表3。
Table 1. Index
表1. 索引
s |
i |
j |
k |
l |
m |
供应商 |
制造商 |
销售区域 |
收集中心 |
二手市场 |
回收厂 |
Table 2. Decision variables
表2. 决策变量
|
|
|
|
|
s至i运输量 |
i至j运输量 |
j至k运输量 |
k至l运输量 |
k至m运输量 |
Table 3. Parameters and their explanations
表3. 参数及释义
参数 |
释义 |
参数 |
释义 |
|
j的需求 |
|
j到k的单位运输碳排放 |
p |
单位产品利润 |
|
k到l的单位运输碳排放 |
|
j收集产品的占比 |
|
k到m的单位运输碳排放 |
|
i固定成本 |
|
s单位采购成本 |
|
k固定成本 |
|
i单位制造成本 |
|
s到i的单位运输成本 |
|
k单位处理成本 |
|
i到j的单位运输成本 |
|
m单位处理成本 |
|
j到k的单位运输成本 |
|
s的产能 |
|
k到l的单位运输成本 |
|
i的产能 |
|
k到m的单位运输成本 |
|
k的容量 |
|
i生产产品单位碳排放 |
|
m的回收能力 |
|
k处理产品单位碳排放 |
|
s供应中断概率 |
|
m处理产品单位碳排放 |
|
I供应中断概率 |
|
s到i的单位运输碳排放 |
|
k供应中断概率 |
|
i到j的单位运输碳排放 |
|
区间置信度水平 |
t |
供应商生产能力水平系数 |
g |
客户需求水平系数 |
市场往往具有不确定性,本文所述的产品价格、单位运输成本、采购成本、等一系列参数值往往会随着市场需求变化而产生波动[8]。为了刻画更真实的市场情况,本文参考文献[9]的处理方式,对所有具有不确定性的参数设定为三个值,分别是
、
和
,其中
表示悲观估计的值,
表示乐观估计的值,
表示最可能估计的值,文中所给定的初始参数值为
。
和
可设定为随机数来控制参数所波动的范围,式(3)即该参数的期望值作为最后设定的参数值。
(1)
(2)
(3)
3. 模型与求解
3.1. 经济性目标
为了提升整个供应链的经济性,供应链企业会考虑降低供应链的经营成本。整个汽车供应链的经济性经营成本主要有六部分构成,分别是:
1. 制造商(i)的收集中心(k)的固定资产成本(例如场地费用、设施费用等)之和
,市场上式(4)中这两个参数值
和
往往会有波动,参照式(3),采用这两个参数的期望值来计算成本
。
(4)
2. 供应商(s)原材料采购的成本以及向制造商(i)运输产品的运输费用之和
,市场上式(5)中这两个参数值
和
往往会有波动,参照式(3),采用这两个参数的期望值乘以该环节运输量
来计算成本
。
(5)
3. 制造商(i)的制造成本和向销售区域(j)运输产品的运输成本之和
,市场上式(6)这两个参数值
和
往往会有波动,参照式(3),采用这两个参数的期望值乘以该环节运输量
来计算成本
。
(6)
4. 收集中心(k)的处理成本和销售区域(j)向收集中心(k)运输产品的运输成本之和
,市场上式(7)中这两个参数值
和
往往会有波动,参照式(3),采用这两个参数的期望值乘以该环节运输量
来计算成本
。
(7)
5. 收集中心(k)向二手市场(l)运输的运输成本之和
,市场上式(8)中的参数值
往往会有波动,参照式(3),采用这个参数的期望值乘以该环节运输量
来计算成本
。
(8)
6. 零件回收商(m)的处理成本和收集中心(k)向零件回收商(m)运输的运输成本之和
,市场上式(9)中两个参数值
和
往往会有波动,参照式(3),采用这两个参数的期望值乘以该环节运输量
来计算成本
。
(9)
总的来说,供应链上的企业在自身产品处理过程中会尽可能降低产生的成本,以达到经济性供应链的目标。
因此,经济性目标函数表达式:
(10)
3.2. 可持续目标
整个汽车供应链的可持续目标主要体现在减少对环境的影响,考虑减少供应链的碳排放来减少对环境的影响,A汽车制造公司供应链碳排放有五部分构成,分别是:
1. 供应商(s)向制造商(i)运输产品所产生的碳排放之和
,式(11)中
表示供应商(s)向制造商(i)运输产品的单位运输碳排放量,乘以该环节运输量
得出
。
(11)
2. 制造商(i)制造产品所产生的碳排放和制造商(i)向销售区域(j)运输产品所产生的碳排放之和
,式(12)中
表示制造商(i)的单位制造碳排放量,
表示制造商(i)向销售区域(j)运输产品的单位碳排放量,乘以该环节运输量
得出
。
(12)
3. 收集中心(k)处理产品所产生的碳排放和销售区域(j)向收集中心(k)运输产品的碳排放之和
,式(13)中
表示收集中心(k)的单位处理产品的碳排放量,
表示销售区域(j)向收集中心(k)运输产品的单位碳排放量,乘以该环节运输量
得出
。
(13)
4. 收集中心(k)向二手市场(l)运输产品的碳排放之和
,式(14)中
表示收集中心(k)向二手市场(l)运输产品的单位运输碳排放量,乘以该环节运输量
得出
。
(14)
5. 零件回收商(m)处理产品所产生的碳排放和收集中心(k)向零件回收商(m)运输产品的碳排放之和
,式(15)中
表示收集中心(k)向零件回收商(m)运输产品的单位运输碳排放量,
表示零件回收商(m)处理产品的单位碳排放量,乘以该环节运输量
得出
。
(15)
总的来说,原材料供应商,制造商,零件回收商的碳排放应尽量降低。供应链上的企业在自身产品处理过程中需要降低碳排放,供应链上的企业在运输产品也应当降低碳排放,以达到供应链的可持续性目标。
因此,可持续目标的函数表达式:
(16)
3.3. 韧性目标
A汽车制造公司供应链的韧性主要体现在降低中断风险造成的中断成本,由此来降低产品的中断供应,提升供应链韧性,韧性目标主要有三部分组成:
1. 供应商(s)向制造商(i)运输产品过程中,中断风险引起的中断成本之和
,式(17)中
和p往往会有波动,参照式(3),采用这两个参数的期望值乘以该环节运输量
来计算中断成本
。
(17)
2. 制造商(i)向销售区域(j)运输产品的过程中,中断风险引起的韧性成本之和
,式(18)中
和p往往会有波动,参照式(3),采用这两个参数的期望值乘以该环节运输量
来计算中断成本
。
(18)
3. 销售区域(j)向收集中心(k)运输产品的过程中,中断风险引起的韧性成本之和
,式(19)中
和p往往会有波动,参照式(3),采用这两个参数的期望值乘以该环节运输量
来计算中断成本
。
(19)
总的来说,供应链上的企业之间的韧性体现主要表现为运输供应过程中的中断成本降低,综合考虑各供应点中断风险来尽量避免中断所带来的损失。
目标三韧性成本的函数表达式:
(20)
3.4. 限制条件
整个汽车供应链的限制条件由以下部分构成,分别是:
所有的供应商向制造商运输的产品数量
的总和不超过自身的生产能力限制,式(21)中
代表供应商的生产能力的值。
(21)
所有的制造商向销售区域收集的产品数量
的总和不超过自身的库存能力限制,式(22)中
代表制造商的制造能力的值。
(22)
所有的收集中心向销售区域收集的产品数量
的总和不超过自身的收集能力限制,式(23)中
代表收集中心的收集量的值。
(23)
所有的零件回收商向收集中心接收的产品数量
的总和不超过自身处理能力限制,式(24)中
代表零件回收商的回收量的值。
(24)
所有的销售区域向制造商接收的产品数量
的总和不超过市场的需求量限制,式(25)中
、
代表销售区域的需求量的悲观值和可能值,
、
代表销售区域的需求量的可能值和乐观值,g代表客户需求系数。
表示在悲观和可能值区间置信度水平,
表示在乐观和可能值区间置信度水度水平,不等式的右边表示最后预估的需求量。
(25)
所有的制造商向供应商接收的产品数量
的总和不超过自身的生产能力限制,式(26)中
代表供应商自身生产能力的值。
(26)
对每个销售区回收量的约束,即销售区运送至收集中心小于等于其回收量
的总和。式(27)中
表示市场上由收集中心回收的产品占比,不等式的左边括号内的式子边表示最后预估的需求量,详细释义参照式(25)。
(27)
式(28)表示制造商从供应商运进量
和运往销售区的运输量
二者总量相等约束。
(28)
式(29)表示对每个销售区域运进量
和运出的运输量
二者总量相等约束。
(29)
式(30)表示收集中心往外运输
的总量相等与二手市场接收量
和零件回收商接收量
的总量约束。
(30)
式(31)表示对收集中心运进量
和销售市场接收量
相等的约束,n表示销售到客户手中的产品被收集中心回收的占比。
(31)
式(32)表示决策变量非负约束。
(32)
3.5. 多个目标化为单目标函数
本研究所采用的方法是交互式求解法,首先分为三个情景计算经济性,可持续性和韧性的目标值。
情景1:优先考虑经济性目标,按照限制条件式(21)~(32)运行式(10)的优化函数
,解出五个决策变量后,由于另外两个目标的目标值均受这五个决策变量控制,可将决策变量代入式(16)和(20)得出情景1可持续性和韧性的目标值。
情景2:优先考虑可持续性目标进行目标值求解。
情景3:优先考虑韧性目标进行目标值求解。
三个情景会有九个目标值,详见下节仿真部分。最后将它们的上限值
和下限值
作为三个目标值上界和下界。得出隶属度函数评估三个目标满意度水平,结合专家建议的权重化为单目标函数
。
式(33)中
表示经济性目标的满意度函数,经济性目标函数值越小,满意度水平越高。
(33)
式(34)中
表示可持续性目标的满意度函数,可持续性目标函数值越小,满意度水平越高。
(34)
式(35)中
表示韧性目标的满意度函数,韧性目标函数值越小,满意度水平越高。
(35)
结合专家所建议的权重值
、
和
,式(36)表示三个目标最后的满意度的优化函数,
受限制于式(21)~(32)
(36)
s.t. (21)~(32)
4. 仿真分析
应用matlab2022b软件,调用cplex求解器,选择A汽车制造公司来展示所建议模型的适用性。A汽车制造公司总部位于上海,在不同区域(所在地为上海、苏州、无锡)拥有三家制造工厂。出于仿真实验目的,假设最终产品小型车单车重量为一吨。专家对经济性目标、可持续目标以及韧性目标三者的权重为0.2、0.4、0.4。p代表单位产品利润(元/吨),p = 9000。供应商、制造区和销售区域的单位运输成本和单位运输碳排放取决于它们所在的地区[10]。通过收集实际市场数据,并运用统计方法对供应中断风险进行预处理,拟合出最后供应中断风险参数的概率分布见表4,求解模型所需的参数见表5。
Table 4. Parameter settings determined for a single location
表4. 单个所在地所确定的参数设置
参数\所在地 |
上海 |
无锡 |
杭州 |
苏州 |
嘉兴 |
(辆) |
8000 |
6000 |
5000 |
\ |
\ |
|
50% |
50% |
50% |
\ |
\ |
(万元) |
2000 |
1400 |
\ |
2000 |
\ |
(万元) |
1000 |
\ |
\ |
800 |
\ |
(KgCO2/吨) |
0.7 |
1.3 |
\ |
1.1 |
\ |
(KgCO2/吨) |
0.9 |
\ |
\ |
1.1 |
\ |
(KgCO2/吨) |
1.6 |
1.7 |
\ |
\ |
\ |
(万元\吨) |
6 |
3 |
5.5 |
5 |
3.5 |
(万元\吨) |
1.4 |
1.1 |
\ |
1.3 |
\ |
(元\吨) |
300 |
\ |
\ |
200 |
\ |
(元\吨) |
800 |
700 |
\ |
\ |
\ |
(吨) |
5000 |
2000 |
2500 |
3000 |
1000 |
(吨) |
8000 |
4000 |
\ |
5000 |
\ |
(吨) |
10,000 |
\ |
\ |
5000 |
\ |
(吨) |
5000 |
3000 |
\ |
\ |
\ |
|
0.15 |
0.24 |
0.21 |
0.17 |
0.28 |
|
0.23 |
0.21 |
\ |
0.26 |
\ |
|
0.18 |
\ |
\ |
0.21 |
\ |
Table 5. Parameter settings determined for two locations. The left value is the unit transportation cost (yuan/ton), and the right value is the unit transportation carbon emission (KgCO₂/ton)
表5. 两个所在地所确定的参数设置。左值为单位运输成本(元/吨)右值为单位运输碳排放量(KgCO2/吨)
所在地 |
上海 |
苏州 |
无锡 |
嘉兴 |
杭州 |
上海 |
60\0.4 |
400\2 |
440\2.5 |
360\1.9 |
740\4 |
苏州 |
400\2 |
40\0.3 |
100\0.5 |
200\1.2 |
460\3.1 |
无锡 |
440\2.5 |
40\0.5 |
50\0.3 |
240\2 |
500\3.3 |
杭州 |
740\4 |
460\3.1 |
500\3.3 |
\ |
\ |
武汉 |
1340\9.1 |
1060\7.3 |
\ |
\ |
\ |
成都 |
2260\18.7 |
1980\16.9 |
\ |
\ |
\ |
Table 6. The objective function values after optimization for three different scenarios when the confidence level α = 0.1
表6. 置信度水平α = 0.1时三个不同情景优化后的函数目标值
情景 |
* (单位/千元) |
* (单位/千克) |
* (单位/千元) |
情景1 |
585,430 |
94,307 |
149,380 |
情景2 |
634,680 |
83,783 |
159,670 |
情景3 |
711,450 |
161,620 |
145,140 |
Figure 1. Comparison of the objective function values after optimization for three different scenarios when α = 0.1
图1. α = 0.1时三个不同情景优化后的目标值对比
Table 7. The objective function values after optimization for three different scenarios when the confidence level α = 0.9
表7. 置信度水平α = 0.9时三个不同情景优化后的函数目标值
情景 |
** (单位/千元) |
**(单位/千克) |
** (单位/千元) |
情景1 |
494,110 |
77,640 |
130,970 |
情景2 |
542,900 |
69,911 |
133,990 |
情景3 |
603,130 |
137,490 |
121,300 |
Figure 2. Comparison of the objective function values after optimization for three different scenarios when α = 0.9
图2. α = 0.9时3个不同情景优化后的目标值对比
为了验证模型的普适性和鲁棒性,选取置信度水平α = 0.1与置信度水平α = 0.9计算各目标值。具体各目标值见表6和表7。见图1和图2均可以得出,如果想考虑经济角度,只考虑最小化运营成本函数是最优策略。另一方面,从可持续性的角度来看,只考虑最小化碳排放是最优策略,而从韧性的角度来看,只考虑最小化韧性成本是最优策略。由此可得供应链的三个目标,即经济性、可持续性和韧性,在本质上都是相互冲突的。
Figure 3. Sensitivity analysis of the resale rate of automobiles
图3. 汽车再售率灵敏度分析
见图3可以得出,三个目标的满意度水平对于汽车再售率n的敏感性不同,整体的满意度水平呈略微上升的趋势。经济性目标的满意度水平在汽车在售率约0.1至0.4时水平较低,而在0.5至1较高,且该段呈缓慢下降趋势。这表明在0.4至0.5汽车再售变化率区间经济性目标的满意度水平对汽车再收率的敏感性程度高。可持续性目标的满意度水平整体呈下降趋势。这表明可持续性目标的满意度水平会随着汽车再售率的上升而下降。整体而言,三个目标满意度水平随着汽车再售率的升高而趋同。
Figure 4. Sensitivity analysis of supplier production capacity level
图4. 供应商生产能力水平灵敏度分析
Figure 5. Sensitivity analysis of customer demand level
图5. 客户需求水平灵敏度分析
由图4可以得出,三个目标的满意度水平对于供应商生产能力水平系数t的敏感性不同,整体的满意度水平呈略微上升的趋势。经济性目标的满意度水平在供应商生产能力水平系数约0.1至0.3时水平较低,而在0.8至1较高,且该段呈上升趋势。这表明在供应商生产能力水平系数变化率区间经济性目标的满意度水平对汽车再收率的敏感性程度高。可持续性目标的满意度水平整体呈下降趋势。这表明可持续性目标的满意度水平会随着供应商生产能力水平系数的上升而下降。整体而言,经济性目标满意度水平随供应商生产能力上升而上升,可持续性目标满意度和韧性目标满意度水平随供应商生产能力上升而下降。
由图5可以得出,三个目标的满意度水平对于客户需求水平系数g的敏感性不同,整体的满意度水平呈先上升后下降的趋势。运营成本目标的满意度水平在供应商生产能力水平系数约0.1至0.3时水平较低,而在0.7至1较高,且该段呈上升趋势。这表明在客户需求水平系数变化率区间经济性目标的满意度水平对汽车再收率的敏感性程度高。可持续性目标的满意度水平整体呈下降趋势。这表明可持续性目标的满意度水平会随着客户需求水平系数的上升而下降。韧性目标满意度水平对客户需求水平变化不明显。整体而言,经济性目标满意度水平随客户需求水平上升而上升,可持续性目标满意度随客户需求水平上升而下降。
5. 结论
本文以含有三个目标的A汽车供应链为研究对象,运用交互式求解法得出三个目标上限值与下限值,然后对三个目标确定各自的满意度函数。通过仿真分析,为避免偶然性,区间置信度水平为α = 0.1和α = 0.9时,三个目标均表现为互相矛盾。采用专家建议三个目标权重的方式将多个目标化为单目标满意度函数,重要参数汽车转售率,供应商生产能力系数以及客户需求对三个目标满意度有着显著影响。本文的研究思路与研究方法能够为双碳背景下含有多目标的供应链研究提供参考。
NOTES
*通讯作者。