基于供应链元跃迁模型的低碳发展多目标决策优化研究.pdf

第40 卷第3 期2023年9 月doi:10.12084/j.issn.2096-3289.2023.03.005基于供应链元跃迁模型的低碳发展多目标决策优化研究苏州科技大学学报（自然科学版)Journal of Suzhou University of Science and Technology(Natural Science Edition)Vol.40 No.3Sep.2023范雅绮1，邵举平1*，孙延安2（1.苏州科技大学商学院，江苏苏州2 150 0 9;2.苏州优乐赛供应链管理有限公司，江苏苏州2 150 2 1）摘要：针对传统供应链向低碳供应链转型背景下如何协调链上核心企业活动的问题，基于供应链元理论，建立了以利润、碳排放和跃迁概率为优化目标的混合整数规划跃迁模型。同时设计约束处理机制，并通过变异策略改进多目标粒子群算法。在此基础上，结合国内某健身器材制造企业的实践运营数据验证了模型与算法的有效性。最后，对模型参数进行了灵敏度分析。研究结果表明：供应链元概念的应用有利于降低研究复杂供应链的难度，基于此建立的模型及其求解算法有效，得到的Pareto最优解集能够提供给决策者具有参考价值的非劣解；对于协调供应链各核心参与者的活动，推动传统供应链向低碳供应链跃迁有一定的理论意义与应用价值。关键词：供应链元跃迁模型；低碳供应链；改进粒子群算法；多目标优化中图分类号：F272.2近年来，环境污染问题日益严峻，国际社会愈发重视绿色低碳发展，我国也积极主动承担大国的减排责任，提出2 0 30 年前实现碳达峰的绿色低碳发展目标，供应链低碳化是实现该目标的一个重要途径，需要供应链上下游各参与者的协同努力。但各参与主体之间既有合作又有竞争，形成的供应链网络错综复杂。因此，协调供应链上核心参与者的活动，推动传统供应链向低碳供应链跃迁，成为供应链管理的一个关键问题。“科学隐喻”思维的形成与发展为众多领域的学者研究问题提供了独特的分析视角。其中,社会科学领域引人物理学中“跃迁”的概念。高锡荣等 2 用量子跃迁比喻创新转型的过程,外界的刺激因素对创新的作用强度与时间共同决定了创新状态的跃迁概率。李晟璐3认为量子跃迁的过程与文化产业发展中的技术体系转型和技术轨道跃迁异曲同工。解学梅等 4就如何在绿色创新中实现本土制造业企业的转型问题进行研究，对绿色跃迁路径提出“长视引领一价值重构一双重预见”的对策。由此可见，“跃迁”被运用于社会科学领域，生动形象地刻画出研究主体不断“积聚能量”的演化过程，最终由低能级向高能级转化的现象。为了适应社会的发展，供应链的转型又何尝不是一种从低阶层向高阶层转化的过程，因此，将传统供应链向低碳化转型视作跃迁的过程无可厚非。目前，出于推进绿色低碳循环发展的需要，供应链低碳转型受到越来越多关注。许多学者考虑了决策目标的差异性，采用多种求解模型算法，研究低碳供应链优化问题。贾旭 5以供应链成本与碳排放量最小为目标，考虑正向运输过程的碳排放量，利用分层次法、8 约束法和加权的理想点法对建立的绿色供应链网络不确定均衡模型进行求解。李进等建立基于整个闭环供应链的碳排放量网络规划模型，采用交互式约束模糊算法对模型进行求解。董海等 7 增加客户满意度损失最小的目标，并采用动态自适应布谷鸟搜索算法求解。YoungsuYun等 8 考虑客户满意度关注对社会的影响，研究供应链三种类型的分布渠道，建模后采用混合遗传算法进行求解。Roy Setiawan等9 以增加总利润、减少环境总影响、社会责任最大化为目标建立数学模型，并将新冠疫情爆发以来的口罩供应链作为研究对象，采用模糊优化的方法对该模型进行了优化。邱云飞收稿日期 2 0 2 2-0 8-2 9基金项目】国家社会科学基金资助项目（19 BGL097)作者简介】范雅绮（19 9 7 一）,女，江苏南通人，硕士研究生，研究方向：物流与供应链管理,E-mail：f a n y a q i 16 3.c o m。*通信作者：邵举平（19 7 1一），男，博士，教授，硕士生导师，E-mail：w l u s t b s h a o 16 3.c o m。文献标志码：A文章编号：2 0 9 6-32 8 9(2 0 2 3)0 3-0 0 36-10第3期等 10 将供应链优化指标细化为经济成本、合规整改与节能减排、绿色供应链、推动公众绿色选择、供应链沟通与透明，并采用蒙特卡洛树搜索改进的分支定界算法求解该可持续闭环供应链网络模型。上述学者基于整个供应链视角研究正向、逆向与闭环供应链的低碳转型优化问题，一些学者则研究某一供应链参与者自身如何进行低碳转型决策优化。在制造商决策优化方面，Nima Farmand等以制造端为起点，考虑生产成本最小化与客户满意度最高的目标,设计两种元启发式算法求解并进行比较。IrfanAli等 12)聚焦供应商选择环节，考虑运输成本、延迟交付以及受总需求约束的净订单成本。JiangH等 13 考虑了云制造情况下的多目标供应链调度优化问题。在经销商决策优化方面，SrikantGupta等 44在限制可控费用不断增长的情况下，同时达到客户期望和销售商竞争力目标，考虑销售商、经销商的双重决策。张烟等15 考虑在共享经济模式下对可循环包装租赁供应链网络进行优化，关注企业可循环包装材料的库存控制。在客户决策优化方面，任腾等 16 考虑客户时间窗限制与碳排放量，基于客户满意度建立具有多配送站点的低碳物流路径规划模型。廖列法等7 以客户时间窗、车辆载重以及冷链产品变质率为约束，在客户服务时间范围内将碳排放量最小作为目标，构建冷链车辆路径优化模型。综上所述，对于供应链低碳化跃迁的研究既有单目标又有多目标。在研究视角方面，既有整个供应链决策优化又有单一供应链参与者决策优化，但总体上，研究侧重于考虑运输过程的低碳优化，而同时考虑供应链计划与供应链低碳化转型的研究尚未成熟。在研究方法方面，学者已采用多种算法求解模型，其中多目标粒子群算法收敛速度快、易实现且易与其他算法相结合，常常运用于无约束条件限制的模型计算中。因此，文中将结合供应链计划问题与供应链低碳化转型问题，建立供应链元跃迁模型，设计约束处理机制，通过变异策略对多目标粒子群算法进行改进，求得Pareto最优解集，协调供应链上核心参与者的活动，推动传统供应链向低碳供应链跃迁。1问题描述与建模1.1问题描述文中所考虑的低碳供应链网络由供应商、核心企业和客户组成，需要综合考虑网络各节点对于低碳产品与传统产品的分别生产批量与库存量、碳排放量以及企业减排投人，以使得供应链整体利润最大化的同时碳排放量最低，协同向低碳供应链进行跃迁。图1所示为节点无限扩展的供应链跃迁模型，图中以某一个节点作为核心企业，同时具有供需关系的节点组成的网络作为供应链中的一个供应链元18 低碳供应链范雅绮，等：基于供应链元跃迁模型的低碳发展多目标决策优化研究节点企业37节点节点企业企业节点企业传统供应链图1的供应链由若干个供应链元组成，一个供应链元的构成如图2 所示。首先建立供应链元的核心企业i的跃迁模型，然后对不同供应链元跃迁模型进行组合，得到完整的供应链跃迁模型，需要考虑i的产品生产批量与库存量、减排投人费用，以使供应链元核心企业i的利润最大化的同时最大程度协同向低碳供应链跃迁。节点企业图1供应链跃迁模型节点企业38苏州科技大学学报（自然科学版)2023年供应商1低碳供应链供应商2客户G核心企业客户供应客户商S21供应商1供应商21传统供应链1.2建模1.2.1模型假设为数学模型的建立，做出以下假设：（1)供应链上主要流动两类产品，分别为低碳产品与传统产品，低碳产品比传统产品在生产与消费的过程中产生的能耗、排放与污染更少，一般所需原材料价格更高，生产过程中需要更高的技术投入，售价较高。（2)减排投入分为两部分：一部分表现为对低碳产品的宣传投入，对低碳产品定向的宣传广告能有效进行营销，从而引导消费者低碳消费，使低碳产品的市场份额得到扩大；另一部分为技术投入，制造商对技术研发进行投资，减少生产过程中的碳排放，从而达到减排目的。（3)供应链上各节点企业对于产品的生产与库存能力有限。（4)供应商作为一个群体存在，客户也作为一个群体存在，跃迁与否同时发生。（5)供应链各节点间供给满足需求。（6)由于可供选择的运输方式不同，运输成本和碳排放也会有所变化，假设是各种运输方式下的平均运输成本与平均碳排放量。（7)产品运输费用由买方承担，原材料采购费用包含原料运输成本。1.2.2集合、参数与决策变量I：整个供应链中供应链元总数；S：i 上游的节点企业的集合;G:i下游的节点企业的集合;t：计划期时间序列(t=1,2,3，,T);b：单位时间内i对低碳产品lp的最大生产能力;b：单位时间内i对传统产品tp的最大生产能力;R;(T):T时期内i生产过程的碳排放限额;qP：单位时间内i对产品p的最大库存能力；W：单位时间内可用于低碳发展的最大可投人资金；JB：i 生产单位低碳产品lp的成本；JP:i生产单位传统产品tp的成本;A：i 对单位产品的库存成本；CP：i 生产单位产品p的碳排放量；Cig：单位产品从i到下游企业g的平均碳排放量;di(t)：:t 时期i向下游企业g承诺交付的低碳产品lp的数量;dig(t)：t 时期i向下游企业g承诺交付的传统产品tp的数量；Pl：i 对低碳产品lp的单位售价；Pp:i对传统产品tp的单位售价;u:i采购用于生产单位低碳产品lp的原料费用;u:i采购用于生产单位传统产品tp的原料费用；K：单位时间内i对低碳产品lp进行宣传的成本;wp：单位时间内i对生产产品p的减排技术投入;x(t)：t 时期i生产低碳产品lp的数量;x(t)：t 时期i生产传统产品tp的数量;y(t):t时期末i对低碳产品lp的库存量;y(t)：t 时期末i对传统产品tp的库存量;w;(t)：t 时期i进行的减排投入。1.2.3目标函数根据问题的描述，文中将建立基于供应链元跃迁的多目标决策模型，分别为最大化供应链利润F、最小化供应链碳排放N和最大化供应链跃迁概率Q三个相互冲突的目标函数。目标1：最大化供应链利润。这里首先考虑供应链元核心企业i的利润最大化，其收入主要表现为产品销售额,要使利润最大,则需将收人去除成本,这里的成本主要包括生产成本(PC)、库存成本(SC)、宣传成本(DC)以及采购成本(BC)。其中,生产成本即生产低碳产品与传统产品所耗成本与技术投入；库存成本指对客户1核心企业i供应商s图2 供应链元跃迁模型G客户2客户1第3期低碳产品与传统产品处于库存阶段花费的费用；宣传成本指的是对低碳产品进行宣传所花费的成本；采购成本指采购用于生产低碳产品与传统产品的原料所花费用。(1)式中TL-2 Zdg(t)P+dg(t)Pr;PC=Zx(t)JI+x(t)JP+w;SC=-Ayr(t)+y(t);范雅绮，等：基于供应链元跃迁模型的低碳发展多目标决策优化研究maxf=L;-PC;-SC;-BC,-DC;T39Tt=1 geG对个不同的供应链元进行组合，得到T内完整的供应链跃迁模型，供应链利润最大化，即可表示为(2)i=1目标2：最小化供应链碳排放。这里的碳排放目标函数由产品制造过程中的碳排放量(PD)与运输过程的碳排放量（TD)组成，考虑供应链元跃迁模型，最小碳排放量表达式如下(3)式中TPD-=ar(t)+x1(t)CP;TD=Z Zdar(t)+dg(t)Cgt=1TTBC=Zx(t)ul+x:(t)u1;DC=Z K/t=1maxF-maxfi！minn;=PD,+TD;Tt=1=1根据完整的供应链跃迁模型，碳排放量应表示为=1目标3：最大化供应链跃迁概率。这里对上游企业、核心企业、下游客户分别向低碳发展进行跃迁的概率求取均值。根据量子跃迁原理,从原始能级向高一能级跃迁的概率为 sin（(IH;l t)/2 h）,其中H为外界刺激强度 19。H用决策生产低碳产品与传统产品的比率来表示，决策生产低碳产品的数量越多，上游受订单刺激，向低碳发展跃迁的概率越高；H。用减排投入来表示，减排投人包括了技术投入与宣传投入,其中，减排技术的投人反映了产业减排技术的发展，减排技术投入越多，系统外部技术作用于企业的效果越显著，越能刺激核心企业由传统向低碳跃迁。而对低碳宣传的投入越多，有利于下游客户对低碳产品的了解，推动客户主动购买低碳产品，对下游客户的跃迁具有直接刺激作用。(5)t=1式中Ha=(x/p(t)/(xp(t);Hg=w;(t)对整个供应链跃迁概率进行考虑，可得maxQ-1)maxqi11.2.4约束条件根据供应链跃迁模型的条件假设，文中的约束条件如下x/(t)b/xP(t)b:yl(t)+yl(t)qlTZxp(t)+x1(t)CPR,(T)t=1t=I.geGmaxN=Z min:Tmaxqi=Ty/p(t)+x/(t+1)-yl(t+1)dig(t+1)yP(t)+xp(t+1)-yP(t+1)dig(t+1)(4)(6)(7)(8)(9)(10)(11)(12)40其中，ViI,VseS,VgEG,式(7)、（8)表示t时期生产低碳产品数量加上前期期末库存余量减去这期期末库存量满足向客户承诺提供的产品数量，传统产品同理；式(9）、（10)表示核心企业的产能约束；式（11)表示核心企业的库存能力约束；式(12)表示生产过程中的碳排放约束；式（13)表示用于低碳发展资金投入约束；式（14)表示对应的决策变量非负；式（15)表示对应的决策变量为整数。2模型求解苏州科技大学学报（自然科学版)K/p+wPw;(t)W;x/p(t),xlp(t),yl(t),yp(t),wi(t)0 x/p(t),x/p(t),yl(t),ylp(t)eZ2023年(13)(14)(15)2.1楼模型计算说明在订单约束条件和低碳跃迁优化中，包含5个决策变量，其中4个为整数决策变量，用来描述供应链网络低碳产品与传统产品数量的流动。对文中构建的混合整数线性规划模型求解是典型的NP-hard问题。由此，文中需要通过设计算法对上述问题进行求解。首先确定供应链元在计划期T内核心企业决策变量的最优解，在此基础上确定供应链T时期内所有供应链元的数量，遍历叠加供应链所包含的节点，即可获得供应链的总利润、总碳排放量、综合跃迁概率。2.2多目标粒子群算法粒子群算法设计粒子，模拟自然界中鸟类的捕食行为，每个粒子具有速度与位置两个属性，随着迭代次数的增多，粒子更新自身的速度与位置，进而实现种群的进化。速度与位置的更新公式如下(16)Pk+I=P k+Vk+1(17)其中,h表示当前迭代次数，表示惯性权重，ci、c 2 表示学习因子,pbest为个体历史最优位置，gbest为全局最优位置,ri与r2是 0,1 间的随机数。2.2.1Pareto 支配关系对于多目标优化问题的求解往往通过采用加权的方法使之转化成单目标问题，从而对单目标进行优化求解,其搜索的结果局限于各目标加权和的最优解,降低了优化各个目标过程进展的可操作性。除此以外，对权重的确定可能也存在较大的主观性。因此,为克服上述方法的局限性，文中采用基于Pareto支配关系的多目标粒子群算法，在单次运行中求得模型最优化解集。为得到多目标优化模型的Pareto最优解集，需要找寻所有非劣解组成的集合。基于Pareto支配关系，初始化后非支配排序。主要步骤如下：(1)计算每个个体被几个个体支配。分别记录每个个体被多少个其他个体支配以及每个个体支配哪些个体。若X为决策空间，个体为x=(x1,x2,xn)的n维决策向量,xX,考虑j个目标函数的最小化问题min(f(x),f(),f(),当f(xa)f(xb),Vie(1,2,j),且日ie(1,2,j),使得f(xa)f(xb),则称向量x.支配向量xbo(2)计算每个个体的非支配等级。记录每个等级的个体与每个个体的等级，寻找支配个数为0 的个体，更新个体等级以及各支配等级的个体。将个体的被支配数目减去1个单位，并将个体对应的支配个体的被支配数目减去1个单位。增加Pareto等级为下一次循环，如果全部计算完毕那么支配数目的值均变为小于0,则停止循环。2.2.2改进多目标粒子群算法的实现上文所构建的多目标规划模型受订货量约束、产能约束、库存量约束、碳排放量约束及资金投入约束条件的限制,在求解的过程中，需要寻求符合这些约束条件的Pareto最优解集。对于此类问题,启发式算法被广泛的应用 2 0-2 ,其中多目标粒子群算法具有收敛精度高且速度快、搜索效率高等特点,尤其适合在动态、多目标优化环境中寻优，可以在较短时间内为问题找寻到较满意的可行解 2 3。但是，粒子的学习机制会使大量Vh+I=wxvi+cixrix(pbest-ph)+c2xr2x(gbest-pk)第3期粒子过快地聚集,早期收敛速度较快，而后期容易陷入局部最优的局面 2 4。因此，为保持种群多样性，增强算法的全局搜索能力，文中将传统粒子群算法与非支配排序结合，并在寻优过程中对粒子位置更新执行变异操作，变异概率不宜太大，通常取值小于0.0 5,文中取0.0 1作为变异概率。同时，由于多目标粒子群算法是一种无约束的搜索机制 2 5,因此将其应用于求解有约束的多目标优化模型时，需要根据问题对约束条件进行合理的编码处理。在多目标粒子群算法的基础上设计约束处理机制，对于小于订货量、超出生产能力限制、库存量限制、碳排放量限制以及资金投入限制的不可行解，通过设置惩罚函数的方式，将惩罚项充分作用于上述不可行解，并且经过多次迭代进而淘汰以上情况出现的不可行解。为保证各个时期的生产量满足订货量,采用设置罚函数的方式将产能约束条件(9）、（10)进行转换，通过罚函数来保证各个时期的不同产品的产量不超出产能限制，并对适应度值Fitness进行重新计算，惩罚超出产能限制的产量，通过该计算评价符合约束条件的解的优劣性，据此得以初始化粒子的个体最优与全局最优，判定粒子之间的支配关系并对外部档案集初始化。其余约束条件处理的方法亦如此。2.2.3算法流程步骤1确定并输入相关的参数，如种群规模popsize、个体学习因子ci和全局学习因子c2、最大迭代次数interations、密度区间网格大小M、惯性权重因子w、缩放因子alpha、存储库规模Ar_N。步骤2根据预设的变量最大值和最小值对种群进行随机初始化。步骤3计算初始种群个体适应度，据此进行非支配排序，将种群分成不同Pareto等级的集合，并根据拥挤距离对同等级集合内部进行排序。步骤4更新外部存储库。步骤5记录目前最优解对应的信息，根据个体极值的选取原则对个体极值pbest进行更新。步骤6 使用多目标粒子群算法更新种群：(1)从外部存储库中随机选择一个个体作为种群历史最优解；(2)采用轮盘赌的方法对领导粒子进行选择确定，从而更新种群历史最优解；(3)对粒子速度与位置的进行更新v(i)=alphax(wxv(i)+cixrix(pbest-p:)+c2xr2x(gbest-pi)(4)若位置超出了范围,则进行修正,对粒子的位置执行变异操作；(5)计算种群内当前每个粒子适应度,并更新外部存储库中个体的历史最优值。步骤7 合并外部存储库与第一序列支配解，通过拥挤距离的计算以及非支配排序对新种群进行排序，保留种群个体Ar_N个并更新外部存储库。记录第一前沿，删除重复的，记录前沿解。步骤8 关判断是否满足终止的标准,是则结束迭代，输出结果，否则转步骤6。步骤9迭代结束后将记录的前沿解数据导出，即得出该问题的Pareto最优解集。3模型验证3.1算例及参数设置选取国内某健身器材产业为研究对象,PVC瑜伽垫与TPE瑜伽垫作为主要生产产品，由于TPE材质的瑜伽垫在生产过程中一般硫化时间较短甚至可以避免硫化，能够很大程度上节约能源，同时，对TPE旧品能够进行回收利用，对增加可再生资源与减少环境污染有利，因此，选取TPE瑜伽垫作为低碳产品，而同类的PVC瑜伽垫作为传统产品。对于一个具有一个上游企业与一个下游企业的供应链元，根据某健身器材制造企业的实践运营数据，利用均匀分布可随机生成模型中各确定参数值，核心企业从上游采购原材料的价格与向下游销售产品的价格独立并且服从正态分布。计划期T=12,具体参数设置见表1。范雅绮，等：基于供应链元跃迁模型的低碳发展多目标决策优化研究41(18)Pi+I=p;+Uit1(19)42对于算例中相关碳排放指标的参数选择,根据文献 2 6-2 8 ,当制造企业的开设规模在 30 0 0,50 0 0 m,该企业的生产碳排放系数服从 0.0 0 0 0 35,0.0 0 0 0 55 t/个区间的均匀分布，产品在各节点间运输的平均碳排放系数服从 0.0 0 0 0 4,0.0 0 0 0 9 t/个区间的均匀分布，中小型制造企业的平均产线最大碳排放限制服从11,15t 的均匀分布。3.2丝结果分析上述算法中，种群规模popsize=40、最大迭代次数 interations=500、惯性权重因子w=0.5、个体学习因子Ci=1和全局学习因子c2=2、密度区间网格大小M=7、缩放因子alpha=0.6、存储库规模Ar_N=20。图3显示了算例的目标函数的迭代过程,图4为最终得到的目标函数的Pareto最优解集。苏州科技大学学报（自然科学版）表1算例参数设置参数取值U1 155,1 290U1 500,2 200U1 000,1 500WU2 000,2 550R,(T)U11,15 U2.3,3.5 U2.1,3.2 AU0.9,1.3CU0.000 035,0.000 055448/5002023年单位参数个C个d.(t)个d.2(t)元PtP元个元/个元个t/个取值U0.000 04,0.000 09U1 068,1 116U1 290,2 000N37,1.27N22,1.14N7.5,0.93N5.2,0.97KU300,350WU1 500,1 730单位t/个个个元/个元/个元/个元/个元元0.60.50.40.30.20.16000碳排放量代图3目标函数送代过程得到目标函数的Pareto最优解集具体情况见表2,对于各决策变量的最优决策，限于篇幅，选择Pareto最优解集中任意一解显示，见表3。123456789101112130.60.590.580.570.560.550.5463.3540004200020总利润元F602 094.42603 622.35603 823.32606 628.70607 662.56609 632.84609823.95613 174.23614 548.96615 357.52615503.57615 806.58617 401.126.26.151053.3碳排放量代图4Pareto最优解集表2 Pareto最优解集N3.373.353.353.333.333.313.313.293.293.283.283.273.266.16.053.256总利润元Q0.601 60.60010.59920.59800.59740.59610.592.70.589.30.588 50.57370.569 20.568 70.5385105第3期x(t)x(t)y（t)y(t)w:(t)3.3楼模型参数灵敏度分析为考察计划期T内模型参数b&、b、q!、WJS、JP、A;d g(t)、d i(t)对目标函数的影响，对各参数进行灵敏度分析。以表2 为基准，使各参数值在原值基础上分别增加或减少5%、10%、15%，每个参数发生变化时，其余参数不变，结果见表4。参数变化幅度-15%b606 051.417b576 802.413q603 567.808总W利J润JAd(t)d.e(t)6碳W排J放J量A;d(t)d(t)bbq跃W迁J概J率A;dt(t)d(t)从表4的结果得出以下结论：(1)b的变化对总利润几乎不产生影响，表明模型对其变化不敏感。b与跃迁概率呈现正相关，随着低碳产品生产能力加强，供应链从传统向低碳跃迁概率越高；(2)b与总利润正相关，模型对其变化较为敏感，与碳排放量弱正相关，对传统产品生产能力越强导致的碳排放量越多，b的变化对跃迁概率几乎没有影响；(3)qP对总利润影响较小,与跃迁概率弱负相关,说明模型对其变化并不非常敏感；（4)W，与总利润强负相关，说明模型对其变化十分敏感，和跃迁概率弱正相关，对低碳投入越多，供应链向低碳跃迁的概率越大；(5)J的变化对利润几乎不产生影响，模型对其变化不敏感；(6)JP与总利润强负相关，说明模型对其变化非常敏感，且随着传统产品的单位成本增加，总利润大幅减少,J!与跃迁概率负相关,模型对其变化较敏感；(7)A,对总利润影响较小,说明模型对其变化并不敏感；(8)di(t)与总利润正相关,对碳排放量和跃迁概率影响较小，说明模型对其变化较为敏感；(9)di(t)与总利润及碳排放量正相关，与跃迁概率负相关，说明模型对其变化较为敏感。范雅绮，等：基于供应链元跃迁模型的低碳发展多目标决策优化研究表3各决策变量的最优决策1201 163231195382889167360418792.270622 649.384609 750.270649 575.939609 959.442578 479.872571 146.5593.351 986 613.231 311 213.249 793 813.235 737 703.328 381 373.459.852.333.250 855 773.345.586 413.274 899 120.593.554.380.603 327 880.611 197 430.599 303.790.598 888 800.611 400 380.601 160 480.583 296 170.621 973.4143349821 1631 864167270578638723021141875表4参数的灵敏度分析-10%-5%606 386.163606 664.665585511.369585552.711610 428.711600 909.616616 889.359609 931.802600 604.867602 794.926636 312.022615 993.838607 038.593606 152.549594.947.366593 748.268585487.943590 644.9483.143 094 873.170.413.893.279 290 143.299 462283.345 122.473.307 951 163.258404483.308 401313.331 612 313.323592283.720.328133.467 944 843.272 311 503.303 479 323.405566 293.365 820 203.287 378 813.301 003 710.598 760 610.599 223 780.608 256 580.600 964 160.607 274 370.605 001 770.599 828 140.600 307 960.599 865 770.594 492 800.606 675 450.608 700 180.610 792 590.603 433 660.608 898520.608 902 490.618 870 240.616 199 1151 1631901872021796116321429351082176602 094.420 2602.094.420 2602 094.420 2602.095.420 2602 096.420 2602.097.420 2602 098.420 2602 099.420 2602 100.420 23.370 078 0363.370 078 0363.370 078 0363.370.078 0363.370 078 0363.370 078.0363.370 078 0363.370 078 0363.370 078 0360.601 572.8790.601 572 8790.601 572 8790.601 572 8790.601 572 8790.601 572.8790.601 572 8790.601 572 8790.601 572 879711632011101752121750%811632139107960619695%609 429.051609 244.401615 372.003606 170.279590 705.638583915.571600 406.977627 129.342637 425.4373.354907723.432754.393.320 102513.339 805293.379 751 653.159 318 523.452 689 503.397 000 823.536316810.603 876 890.606 481 260.600 692 520.604 922 160.604 267 780.593835.850.632 063 960.609 811990.589 597 94994921429248362136607 725.893617 267.735619.795.099603 218.645606 421.577585385.661604 042.048627 826.926646 821.3303.330 739 763.431 417 863.384 463 253.314383933.463 763 343.252.754 513.444 078 233.394 696 613.600 007 780.606 710 940.603 795 330.589522720.605 223 590.608 256 840.593 082.340.600 163 870.609 773 100.571 899 291011632127843880225910%115161917233387203515%609 441.951619.987.160612 670.787591801.345593 593.075574.780.891609 612.591629 123.908646.979.2893.225256023.702 608 883.301 576 323.218 421 243.477 652 413.322.016213.441 380 553.388 169 883.799553130.609 838920.609 754 810.588 843 310.607 891 840.592517 580.587416870.628 720 220.596 883 090.571 338 9412889158701012055444结语苏州科技大学学报（自然科学版）2023 年文中研究了供应链低碳化转型中的多目标决策优化问题，建立了供应链元跃迁模型，该模型以供应链利润、供应链碳排放和供应链跃迁概率为优化目标，设计有效约束处理机制，对多目标粒子群算法进行改进，文中以国内某健身器材制造企业为研究对象，确定该企业各时间节点选择生产各产品的产量、人库量以及各时间节点用于减排的资金投人，从而得到Pareto最优解集。研究结果表明，基于供应链元跃迁的低碳供应链多目标决策优化模型及其求解算法是有效的，得到的Pareto最优解集能够为决策者提供尽可能多的具有参考价值的非劣解。对于协调供应链上核心参与者的活动，推动传统供应链向绿色低碳供应链跃迁提供理论依据与实践意义。同时，供应链元概念的应用降低了研究复杂供应链的难度。参考文献：1项目综合报告编写组.中国长期低碳发展战略与转型路径研究综合报告.中国人口资源与环境，2 0 2 0,30(11)：1-2 5.2高锡荣，柯俊.基于能级跃迁的创新转型激发模型 .科技进步与对策,2 0 16,33(3):7-11.3李晟璐.数字经济时代文化产业技术轨道跃迁的量子隐喻研究 .经济问题,2 0 2 1(7):12 0-12 8.4解学梅，韩宇航.本土制造业企业如何在绿色创新中实现“华丽转型？一基于注意力基础观的多案例研究 .管理世界，2 0 2 2,38(3)：7 6-106.5贾旭，刘诚.随机和模糊环境下绿色供应链网络优化设计 .铁道科学与工程学报，,2 0 18,15(3)：7 9 2-8 0 1.6李进，朱道立.模糊环境下低碳闭环供应链网络设计多目标规划模型与算法 .计算机集成制造系统,2 0 18,2 4(2)：49 4-50 4.7董海，高秀秀，魏铭琦.基于动态自适应布谷鸟搜索算法的多目标闭环供应链网络优化 .系统工程,2 0 2 0,38(4)：46-58.8 YUN Y,ANUDARI C,MITSUO G.Sustainable closed-loop supply chain design problem:A hybrid genetic algorithm approachJ.Mathematics,2020,8(1):17-28.9 SETIAWAN R,SALMAN R,KHAIROV B,et al.Sustainable closed-loop mask supply chain network design using mathematical modeling and afuzzy multi-objective approachJ.Sustainability,2021,13(10):21-34.10邱云飞，于智龙，郭羽含，等.蒙特卡洛树搜索下的整合多目标可持续闭环供应链网络优化 .计算机集成制造系统，2 0 2 1,2 8(1):2 6 9-2 9 3.1 NIMA F,HAMID Z,MORTEZA R.Two meta-heuristic algorithms for optimizing a multi-objective supply chain scheduling problem in an identi-cal parallel machines environmentJJ.International Journal of Industrial Engineering Computations,2021,12(3):249-272.12 ALII,FUGENSCHUH A,GUPTA S,et al.The LR-type fuzzy multi-objective vendor selection problem in supply chain managementJ.Mathe-matics,2020,8(9):27-39.13 JIANG H,LIU C.Scheduling optimization of cloud resource supply chain through multi-objective particle swarm optimizationJ.InternationalJournal of Simulation Modeling,2019,18(1):18-31.14 SRIKANT G,AHTESHAMUL H,IRFAN ALI,et al.Significance of multi-objective optimization in logistics problem for multi-product supplychain network under the intuitionistic fuzzy environmentJJ.Complex&Intelligent Systems,2021,12(9):31-44.15张烟，李晓华，邵举平，等.共享经济环境下可循环包装租赁企业供应链网络库存控制研究 .苏州科技大学学报（自然科学版）,2 0 19，36(4):75-84.16任腾，陈玥，向迎春，等.考虑客户满意度的低碳冷链车辆路径优化 .计算机集成制造系统,2 0 2 0,2 6(4)：110 8-1117.17廖列