论文-多目标优化模型

数学模型港口物流服务供应链（Port Logistics Service Supply Chain, PLSSC）是以 港口企业（港务集团公司）为核心企业，将运输、装卸搬运、流通加工、仓储、报 关、配送、金融、商业服务等各类物流服务供应商（LSP）和客户（付货人和收货人 等） 有效结合成一体，是一个以提供物流服务为主要功能的服务型供应链7。港口 供应链的集成优化所追求的不仅仅是各环节（各节点企业）的局部优化，更是供 应链的整体最优。本文建立的港口供应链集成优化模型的目标有三个：1）供应链运行总成本最小；2）供应链服务时间最短；3）供应链柔性最大。此外，在供应链网络中，供应链节点企业之间的距离也是选择供应链协作企 业的一个考虑因素，距离越近，移交时间越短，成本越低，越有利于提高供应链 的优化水平。本模型将距离因素加入到成本目标和时间目标的函数里，提高模型 在现实中的实用意义。2.1 建立模型的假设条件1）为方便模型建立，港口供应链仅考虑其主要节点企业：供应商、客户、港 口、船运公司、加工商、仓储服务供应商、配送服务供应商。同样的，服务成本 仅考虑五个主要的服务环节：运输成本、港口服务成本、加工成本、仓储作业成 本、配送成本。2）模型中提到的港口服务仅指港口的基本服务，如装卸服务等。3）在一定时期范围内，成本参数是稳定的。4）供应链服务节点企业之间的距离，将会使货物在移交过程中发生运输费用， 产生运输时间，模型设定运输费用仅与两节点间的距离和货物量有关，运输时间 仅与距离有关，单位运输成本系数和时间系数设置为恒定。P临近销售地B。2.2 设定模型参数及决策变量k货物种类索引号，kg1,2,K,K； a：货源地（供应商）；b销售地（客户）；p：港口（港口装卸服务提供商）的索引号，pw1,2,K,p；e：船运公司的索引号，ew1,2,K,E；i：加工服务供应商的索引号，iw1,2,K,I；j仓储服务供应商的索引号，jw1,2,K, J；g：配送服务供应商的索引号，gw1,2,K,G；v：物流服务供应商的索引号，v e p,e,i, j, g；f构成港口供应链的备选成员的索引号，f ea,b,p,e,i, j,g；m,n :三类供应链服务节点（港口、加工商、仓储服务供应商）的索引号，货物可以从这三个节点处运往销售地，m,n ep,i, j且m主n ；L :货源地a到港口 p的距离；apL :供应链服务节点m到销售地b的距离；mbL :供应链服务节点m到n之间的距离；mnCU : LSPv为货物k提供服务的成本；vkCR :两个供应链服务节点m、n之间移交货物所产生的运输成本；« :货物移交的单位运输成本系数；TU : LSPv为货物k提供服务所需的时间（仓储服务时间仅指作业时间，如出入 vk库时间、拣货时间等，库存时间不包含在内）；TR :两个供应链服务节点m、n之间移交货物所需要的时间；0 :货物移交所需要的时间的系数；T : LSPv 向货主承诺的服务时间；vk厂:LSPv对货物k完成服务的实际时间；vkM ： LSP v 所能提供的最大服务能力；v3 :供应链上各种不同类型服务的柔性权重系数；vX :货物k的需求量；kC :供应链运行总成本T :供应链服务时间F :供应链柔性Y :是一个0、1变量，当其为1时，表示备选成员f被选为港口供应链上的协作企f业，否则为0；2.3 目标函数2.3.1 供应链运行总成本最小min C = Z 另e=1 k=1YCU (X , L ) +e ek k app=1 k=1Y CU (X )p pk k+工另YCU (X ) +为另i ik ki=1 k=1j=1 k=1Y CU (X )j jk k+另另 YCU (X ,L)+CR g gk k g =1 k =1其中，1）1）第一种情况：a-p-i-j-b （从港口卸货后，加工，仓储，配送）即当工Y =1且为Y =1时,L=为另YYL ,CR= 艺a X (YS Y YL +XYYL )ijj b jbkp i pii j iji=1j=1j=1 b=1k=1p=1 i=1i=1 j=12)第二种情况：a-p-i -b (从港口卸货后，加工，配送)即当 工Y =1且为 Y =0时iji =1j =1L= YI YB YY L ， CR= YK a X (YP YIi b ibki=1 b=1k=1p=1 i=1YYL )p i pi3）第三种情况： a-p- j-b （从港口卸货后，仓储，配送）即当工Y =0且为Y =1时，L= YSyYL ，CR= YKaX (YPYJYYLii=1jj b jbkp j pjj=1j=1 b=1k=1p=1 j=14）第四种情况：a-p-b （从港口卸货后，直接配送）即当 工Y =0且为Y =0时，L=另为YYL , CR=0；i j p b pbi =1j=1p=1 b=1港口供应链运行总成本C包括：供应链上各环节的服务成本（海运成本、港 口装卸成本、加工成本、仓储作业成本、配送成本）；以及两服务节点间移交货 物产生的成本。式（1）是表示供应链总成本最小的目标函数。式（2）是求货物 配送距离的函数，表示若货物不需加工、仓储服务，则直接从港口卸载后运往销 售地；若货物只需加工不需储存，则货物直接从加工地运往销售地；若货物要储 存，本文设定无论货物是否需要加工，均可认为货物是从仓储地运往销售地（货 物在仓储后配送前进行加工的地点一般仍是仓储地点）。式（3）是货物的移交成 本函数，与距离和货物量有关。2.3.2 供应链服务时间最短min T =艺另 YTU (L ) + 另另 YTU (X )eek apppkke=1 k =1p=1 k =1+工另 YTU (X ) + 为另 YTU (X )i ik kj jk ki=1 k =1j=1 k=14）+另另 YTU (L) + TRg gkg =1 k =1其中，工另j b jbj=1 b=1YYL ，TR= 0 (ZSyYL +ZSyYL ) p i pii j ijp=1 i =1i =1 j=1工另YYL ,i b ib i=1 b=1TR= B (另丫 Y YL )p i pi p=1 i =11）第一种情况： a-p-i-j-b即当 工Y =1且为 Y =1时， L=iji=1j=12）第二种情况： a-p-i -b即当 工Y =1且为 Y =0时， L=iji=1j=13）第三种情况： a-p- j-b即当工Y =0且为Y =1时，L=为另Y YL , TR= 0 (另艺 Y Y L )ijj b jbp j pji =1j =1j =1 b =1p =1 j =14)第四种情况： a-p-b即当 工Y -0且为Y =0时，L=另为 YYL , TR=0；i j p b pbi =1j=1p=1 b=1港口供应链服务时间T包括：供应链上各环节提供服务所需要的时间（海运 时间、港口服务时间、加工时间、仓储作业时间、配送时间）；两服务节点间移 交货物所需的时间。式（4）是求供应链服务时间最短的目标函数。式（5）是货 物的移交时间函数，与两服务节点间的距离有关。2.3.3 供应链柔性最大maxF =3 另Y (M -Q) + w 1Ly (M -Q) + ® 工Y(M -Q)pp pee eii ip=1e=1i =1+3 1Ly （M -Q） +3 另Y （M -Q）（6）j j jg g gj=1g =1港口供应链柔性F由五个部分组成：港口柔性、运输柔性、加工柔性、库存柔 性、配送柔性。式（6）是求供应链柔性最大的目标函数。供应链柔性的大小，与 各服务环节所能提供的服务能力与各环节货物所需的服务量之间的差距空间有 关，柔性权重3 可以通过对比各服务环节的相对重要程度获得。2.3.4 总目标函数minZ =g C + g T g FCTF2.4 约束条件0<1X <M（7）kvk=10<T' <T（8）vk vk艺 X = Q，Vk g 1,2,K, K（9）kk=13 +3 +3+3 +31=（10）p e ij gg +g +g = 1C T F艺Y = 1,另Y = 1,工Y = 0或1，为Y = 0或1，另Y = 1epijg其中：式（7）是能力约束条件，表示各LSP为货物提供的服务量应不大于其 最大服务能力，如港口吞吐量不能超过港口的通过能力。式（8）是时间约束条件， 表示各LSP实际的服务时间不应大于向客户承诺的服务时间。式（9）和式（10） 是均衡约束条件，式（9）表示各类货物量的总和应等于在供应链上运行的总货物 量；式（10）表示供应链上各类服务柔性的权重系数之和应等于1。3 粒子群算法上述的港口供应链集成优化模型是一个多目标优化模型。目前，解决多目标 优化问题的比较著名的算法有遗传算法、粒子群算法等。由于遗传算法的编程实 现复杂且收敛速度慢，本文采用粒子群算法解决多目标优化问题。粒子群优化算法8是一种基于迭代模式的优化算法，其基本思想是将优化问题 的每个潜在解看作D维搜索空间的一个“粒子”。粒子在搜索空间中以一定的速度 飞行，这个速度根据它本身的飞行经验和同伴的飞行经验来动态调整。所有的粒 子都有一个被目标函数决定的适应值，并且知道自己到目前为止发现的最好解， 叫做个体极值点（用Pbes t表示其位置）。这个可以看作是粒子自己的飞行经验。 除此之外，每个粒子还知道到目前为止整个群体中所有粒子找到的最好解，称为 全局极值点（用gbes t表示其位置），这个可以看作是粒子的同伴的经验。每个 粒子根据自己的当前位置，当前位置与自己最好位置之间的距离，以及当前位置 与群体最好位置之间的距离来进行迭代，改变自己的当前位置，不断迭代搜索寻 找最接近最优解的帕累托解。为了更好的控制PSO的探测、开发能力，Eberhart and Shi9提出了对基本粒 子群算法的改进，即对速度更新方程加惯性权重w。权重w将影响的全局和局部寻 优能力，较大的w可以加强全局搜索能力，而较小的w能加强局部搜索能力。本文在改进的PSO算法的基础上，设计适合港口供应链多目标优化模型的算法。将港口供应链上的备选成员作为D维空间的粒子。为了保证粒子的每一维数值为备选企业，粒子的速度向量的每一维都必须取整数，整数区间为供应链各节点企业的备选个数。如设定Y是代表船运公司的粒子，Y的速度向量的每一维都取区ee间1,E上的整数。并且在每次迭代后，需要对粒子进行向上的取整处理。粒子的 速度和位置是有上下界限制的，如的速度取值范围为-（E-1）,（E-1），位置的取 值范围为1,E。以港口供应链集成优化模型的三个目标函数作为适应值函数。用改进的粒子群算法实现港口供应链多目标优化的算法步骤如下：Step1: 生成种群，初始化粒子的位置和速度，以及群体规模。Step2: 对每个粒子进行速度和位置的更新。位置更新：vk+i 二 wvk + c r (pbestk xk) +