物流系统的LRP改进多目标线性规划模型

第46卷第2期2011年4月

西南交通大学学报

JOURNALOFSOUTHWESTJIAOTONGUNIVERSITY

DOI:10.3969/j.i鹊n.0258-2724.2011.02.025

V01.46

No.2

Apr.20ll

文章编号:0258_2724(2011)02旬326JD7

危险废物物流系统的LRP改进

多目标线性规划模型

斌,

赵佳虹

(西南交通大学交通运输与物流学院,四川成都610031)

摘要:为了解决危险废物回收、加工和处理中心选址问题,确定加工技术类别、安排危险废物和废物残余车辆运输路径,将回收环节纳入危险废物物流系统,考虑废物类型与加工技术的多样性、中心运营费用、废物与加工技术的相容性以及中心能力富余量约束,以费用和风险最小化为优化目标,建立了危险废物物流系统的改进多目标0-l混合整数线性规划模型.采用TOPSIs(technique

7.69%.

fororderpreference

bysimil蚵tyto锄ideal∞lution)方法

求解模型.结果表明,与现有模型相比,本文模型的多目标优化方案以增加0.70%的费用为代价,将风险降低

关键词:多目标模型;线性规划;危险废物;选址-路径问题;阳PsIS方法

中图分类号:F253.9

文献标志码:A

Multi・objectiVeO-1LinearProgramlIling

Model

forCombiI抡dLocation・RoutingProbleminHazardou【s

Waste

LogisticsSystemJ证h0鸭

sHUAIB流,ZHAo

(Sch∞l

0f。r协sportati∞锄db西stics,SoutllwestJiaotonguniVersity,chengdu61003l,china)

Abs咖ct:

Animprovedmulti-objective

tIIe

0-l

mixed-integerline盯pmgrammingmodelf打tllelocation—

was

mutingprobleminh北ardouswasteminimizetotalcosts蛐drisks,and

lo舀sticsproposed,inwhichtheoptimizationobjectiVes

areto

cons嘲nts

cen£ers.was

includewastetypesandtreatmenttechnologies,ATOPSIS

to

c印acityandoperation

costs

oftreatment

algorithm

(technique

t11is

for

orderpreference

model.

byThe

similarit)rto舳ideals01ution)

designed

solVe

multi-objectiVe

feasibilityandadVantageoftIlepr{0posedmodelw鹊demonstratedthrough

tal【enf而mliterature.

Compadngwit}lof

an

representativeexample

tlle

risk

by

exitingmodel,the

proposedmodelreduced

7.69%atthe

Keywords:

expense

littleriseincostby0.70%.

model;1inearprogramming;haza—ouswaste;location-routingproblem;

multi-objective

ToPSISmeth()d

危险废物是指具有腐蚀性、毒性、易燃性、反应性或感染性等危险特性的固体或液体废物;或者可能对环境或人体健康造成有害影响的物质.危险废物的存在对环境和人类生活形成潜在威胁,经适当

加工处理后,部分危险废物可实现回收利用.

选址一路径问题(10cation—muting

problem,

LRP)是在确定设施选址位置的同时,制定车辆运输路径,即协同优化设施选址与车辆运输路径两个

收稿日期:2010彤.10

基金项目:四川省经济委员会技术创新资助项目(200901120)

作者简介:帅斌(1967一),男.教授,博士,博士生导师,研究方向为交通运输经济、交通与物流规划,电话:028-87603464,

E-mail:bs№7@126.o∞

第2期帅斌等:危险废物物流系统的LRP改进多目标线性规划模型

327

NP(non—polynominal)问题.近年来,学术界就有关危险废物物流LRP问题发表了大量文献.文献[1]研究了废核燃料的运输和储运点选址问题,采用分割技术将运输路径的选择从选址问题中分离出来单独进行优化,以实现货物吨公里和暴露人口数量的最小化.文献[2]以暴露人口数量表示设施选址风险,并设定风险限界,以最少费用为目标建立模型,求解时引入拉格朗日系数,将约束条件转化为目标函数的多项式,将LRP分解为选址和车辆路径选择两个子问题.文献[3石]主要针对加工技术处理危险废物环节的LRP进行优化,同时考虑费用和风险目标.

以上文献都只考虑了一种加工技术.文献[7]则根据加工技术的多样性和技术选择的模糊性,以加工技术的效率为参数度量选址风险.文献[8]详细分析了危险废物管理体系,设计了包含废物收集、加工、处理和运输环节的物流管理系统,并对各环节的风险和费用进行度量.文献[9]基于危险废物物流管理系统,重点考虑了危险废物类型和加工技术的多样性与相容性,建立了费用和风险的双目

标0-l混合整数线性规划模型,在确定加工中心和

处理中心位置的同时,决定加工中心使用的加工技术,并采用线性加权方法转化模型,借助数学优化软件进行实际问题的演算;但该模型的成本目标函数没有包括设施选址的变动成本(例如运营费用),风险目标函数没有包括因设施选址而给相关区域带来的风险(选址风险),也没有考虑危险废物产量的波动性.

基于上述文献,本文研究危险废物物流系统u讲时,特别考虑了危险废物类型和加工技术的多样性、危险废物与加工技术的相容性、相关中心的

变动成本,并且在中心能力约束中引入富余度系

数,以解决危险废物产量的波动性问题.1

问题描述

如图1所示,在危险废物物流系统中,危险废

物由废物源点运往加工中心;加工中心将可回收的废物残余运往回收中心,将不可回收的运往处理中心.

该系统的LRP是在考虑危险废物类型和加工技术的多样性(每种加工技术都有不同的特性¨¨11J,即加工技术与危险废物存在相容性),以及废物源点各类危险废物产量的波动性的前提下,

从候选点中确定各中心的建设位置,确定加工中心采用的技术类别,并制定危险废物和废物残余在各中心之间的车辆运输路径,使系统总费用和风险最小.其中,费用包括运输和选址费用,主要为危险废物和废物残余运输费用、各中心建设费用和运营费用.风险用暴露人口吨数表示,包括危险废物和废物残余运输风险,以及各中心选址风险.

加工中心

处理中心

图1危险废物物流系统

Fig.1

H北a—ous

waste

logisticssystem

2数学模型

以危险废物物流系统为基础构建运输网络,网络节点包括危险废物源点、各类中心建设候选点和其他运输节点.网络中分布有居民数量不同的居民点,它们都会受到运输风险和选址风险的影响.

假设:

(1)网络各节点之间的道路安全等级满足运输危险废物和废物残余的要求,且无道路能力限制;

(2)危险废物满足运输相容性条件,即不相容废物不同时由同一车辆运输;

(3)各类危险废物都有与之相适应的加工技术;

(4)各类危险废物都有特定的运输车辆,废物与车辆满足相容性条件,且车辆无能力和工作时间限制;

(5)各加工中心具备同时处理危险废物和废

物残余的能力.

令Ⅳ(y,E)代表网络,y代表节点,E代表道路,有:

y=%u¨uKu%u%,

式中:%——危险废物源点;

K——加工中心建设候选点;%——回收中心建设候选点;%——处理中心建设候选点;%——其他运输节点.

埘(1,2,…,形)——各类危险废物编号;Z(1,2,…,£)——各类加工技术编号.

328西

交通大

学学

报第46卷

min

A。=∑∑c。。Q。旷∑c口吼+∑c’F%+∑∑c讷。+(∑∑∑ch。Q‰。)+

”EⅣ(tJ)E盂

(iJ)EE

(iJ)EE

iEr

lEL

、"E甲iEh

iEL

(墨c瞰y一;;c瞰Qm)+(磊c碱+i善cⅨQⅨ),

、iE慷

iE‰

(・)

、‘EyDlEyD

min

A2=∑

‘E

∑P,FQ。#+∑P口Q口+∑P’FQ’p+∑∑∑P删Q州+

IIJ,EB

‘‘JJ

EB

钾EⅣLIJJE6"EⅣ●EyT●EL

∑PmQm+∑P。iQm,

yR

lE

2@

V埘∈形,

Vi∈y,

3o4H

yD

s.t.

Q铡=

∑Q。#一∑Q叫+∑QT训,

【‘JJ

E£

U・‘J

E6

‘EL

∑∑(1.a“)QT。艘=∑Qg一∑徘+Qm,Vi∈y,

埘E咿IEL

(‘J)E占U,1)E最

”EⅣ●EL

∑∑(1一a“)Q,。。(1一卢。)=∑Q’#一∑Q’一+QⅨ,

【IJ)E点

‘』,1)E邑

Vi∈y,

L,

5@6∞78

QT

mi。Ⅱ并“≤∑Q。d≤(1一叩)QT。。i∥“,

Vi∈坼,

Vz

QR“。‘y‘≤Q瞰≤(1一叼)QR。。‘,,‘,QDmi。‘气≤Q叭≤(1一田)QD~i气,Q州≤(1一叼)QT—d0“,

VtlI∈形,

Vi∈K,ViE%,

Vi∈%,

VZ∈L,

∑∑Q,埘=o,

Q砌=0,QDi=0,

Vi∈(y一%),

∽@@m0

12

ViE(y—K),ViE(y一%),

V训∈形,

Vi∈n,

VZ∈£,

0≤QT甜≤慨d,

¨坦B3

45

0≤Q砌≤嘶;,Vi∈K,

0≤QⅨ≤胞i,Vi

Q埘≥O,

Q“≥0,

xⅡ∈{0,1},),‘E{0,l},毛∈{O,l},

E%,Q’#≥O,

V|t7∈形,

V(iJ)∈E,

678

Vi∈¨,

VZ∈L,

Vi∈K,Vi∈%,

t£7∈形,i∈%;

掩挎M巧№"9

式中:Al——总费用;

A2——总风险;

c。。——危险废物埘在弧(iJ)上的单位运输成

本,加∈形,(iJ)∈层;

a“——危险废物加用技术Z加工后质量减少

的百分比,埘E矽,Z∈£;

岛——用技术z加工后可回收的废物残余所

占百分比,Z∈£;

c。和c’;,——可回收和不可回收的废物残余在

弧(i√)上的单位运输成本,(i√)∈E;

Q,。删和Q,一。——加工中心i采用技术f的

最小加工量和最大加工能力,iEH,Z∈L;

C,。——加工中心i采用技术f的固定成本,

f∈坼,Z∈£;

Q。。¨和Q。一;——回收中心i的最小回收量

和最大回收量,i∈K;

C磁——建设回收中心i的固定成本,i∈K;C珑——建设处理中心i的固定成本,i∈%;

cT。。——危险废物tc,在加工中心i用技术z加

工的单位成本,叫E形,i单位成本,i∈‰;

Q。“。;和Q。一。——处理中心i的最小处理量

和最大处理能力,i∈%;

n,Z∈£;

竹——中心能力的富余度系数,0≤叼≤1;

相容性,则为1,反之为0,埘∈形,Z∈L;

c。;——回收中心i回收的可回收废物残余的cⅨ——处理中心i处理的不可回收的废物残

余的单位成本,i∈%;

@“刈.1变量,若危险废物埘与技术Z满足

P"。——危险废物加在弧(i,.『)上运输时暴露

的人口数量,埘∈形,(i√)∈E;

QG“——废物源点i产生危险废物埘的总量,P。,和P’。——可回收和不可回收的废物残余

第2期帅斌等:危险废物物流系统的LRP改进多目标线性规划模型

329

在弧(i,_『)上运输时暴露的人口数量,(i,_『)∈E;

PT。——建设采用技术z的加工中心i时暴露

的人口数量,f∈¨,Z∈£;

P瞰——建设回收中心i时暴露的人口数量,

i∈K;

P陇——建设处理中心i时暴露的人口数量,

i∈%;

胁一无穷大的正整数;

Q甜——危险废物加在弧(i,J)上的运量,埘∈

形,(iJ)∈E;

QⅡ和Q’;f——可回收和不可回收的废物残余

在弧(iJ)上的运量,(;,_『)EE;

Q‰。——危险废物埘在加工中心i用技术Z加

工的量,埘∈形,i∈坼,Z∈£;

Q砌——回收中心i回收的废物残余量,

i∈‰;

Q叫——处理中心i处理的废物残余量,

i∈%;

幻刈一l变量,若加工中心i采用技术l,则为l,反之为0,iyi刈-l变量,若建设回收中心i,则为1,反EK,Z∈L;之为0,i∈K;盈刈一1变量,若建设处理中心i,则为1,反

之为0,iE%.

式(1)和(2)表示最小化费用和最小化风险.

式(3)一(19)表示约束条件.式(3)~(5)分别表示危险废物、可回收和不可回收废物残余的节点流量守恒.式(6)表示加工中心各技术的最小加工量和带富余度系数的最大加工能力约束.式(7)表示回收中心的最小回收量和带富余度系数的最大回收能力约束.式(8)表示处理中心的最小处理量和带富余度系数的最大处理能力约束.式(9)为危险

废物与加工技术相容性约束.式(10)一(12)表示

各中心只在候选点中选址.式(13)表示危险废物只在加工中心加工.式(14)表示可回收的废物残余只由回收中心回收.式(15)表示不可回收的废物残余只在处理中心处理.式(16)一(19)分别表示连续决策变量和0—1离散决策变量的定义域.3

多目标优化问题的解是一个含有Pareto最优解的集合,而单目标优化问题的解唯一.因而,常使

用线性加权法、目标规划或动态规划方法‘12’143等,

将多目标优化问题转化为单目标优化问题进行求

解.TOPSIS(techniquefororderp11eferencebysimil撕tyto蚰idealsolution)方法,又称为优劣解距离法,是一种用于多目标决策分析和求解多目标优化问题的有效方法.

本文模型为多目标0-1混合整数线性规划模型,费用和风险两个子目标的量纲不同,且子目标的权重系数具有一定的模糊性.因而,采用极值法消除量纲,利用TOPSIS法¨51,将该模型转化为单目标线性规划模型进行求解,最终求得原多目标

0—1

7昆合整数线性规划模型的有效解.求解算法步骤如下:

步骤1令X为多目标优化问题的决策变量

向量;F(x)代表多目标优化问题的目标函数;五(x)(庇=1,2,…,珏)为子目标对应的单目标优化问题的目标函数;G为多目标优化问题的定义域.

步骤2求解该多目标0-1混合整数线性规划模型的目标函数,每个子目标对应单目标优化问题

{五(x)=min以(x)lx∈G},得到最优目标函数

值爿,.|}=l,2.

步骤3采用极值法,引入转换系数

弧2]万’肚l'厶

五(x)一片

,.

使费用和风险的无量纲化.

步骤4运用TOPSIS法,将该多目标模型转化为单目标线性规划模型.该单目标线性规划模型为:

min

s.t.舻噤,㈦,2’x∈G.

F(x)=./王{饥眈(x)一五]}2,

s・t.

7t2—■■1夏■一,

丘2l,z’

五∈u・

.,^\√1,

步骤5求解上述单目标线性规划模型,求得

多目标线性规划模型的有效解.得到选址一路径方案.

4算例

综合文献[16—17]的算例生成方法,并沿用文献[16-17]中的危险废物年产量和加工技术的最大加工能力,依据物价指数处理其相关费用数值,给出危险废物物流系统的LRP算例,如图2所示.

图2中,危险废物物流系统运输网络共有10个节点、20条道路,网络中分布有居民数量不同的居民点.

图2中,节点1—3为危险废物源点,节点4、5为回收中心建设候选点,节点6—8为加工中心建设候选点,节点9、lO为处理中心建设候选点.括号

330

西南交通

学条件.

报第46卷

中数字表示运输距离(单位:km)和暴露人口数量(单位:万人).

危险废物年产量见表1.3类危险废物运费分别为30、35和40元/(t-km),可回收和不可回收

废物残余的运费分别为20和15形(t・km).设定

中心能力富余度系数为0.15,中心建设候选点信息见表2~4.

表l危险废物年产量

Tab.1

Hazardous

waste

gene瑚土ionqu帅tity

图2算例网络结构

Fig.2

E】【啪pletr锄sportati∞ne附ork

固化技术减少的危险废物质量约为15%~30%,其生成的可回收废物残余约占剩余废物总质量的20%一40%;焚化技术减少的废物质量比固化技术高1.1—1.3倍,但可回收废物残余量很少‘1

81.

在节点之间道路运输各类型危险废物、废物残余的可回收部分和不可回收部分时,暴露的人口数量相同.该网络中有3类危险废物,分别为金属、农药和化工品,可选加工技术有两种:固化技术和焚化技术,其中金属废品与焚化技术不满足相容性

加工技术的效率见表5.

表2废物加工中心建设候选点

7I'ab.2

7I协tlll即tcenter

c帅didates

注:最小加工量和暴露人口数量为设定值;物价指数参照《中国统计摘要》(2003—2009).

表3处理中心建设候选点

’rab.3

Disposal

center

c蛐didates

注:最小处理量和暴露人口数量为设定值.

表4回收中心建设候选点

Tab.4

Recycling

center

candidates

注:最小处理量和暴露人口数量为设定值.

第2期帅斌等:危险废物物流系统的LRP改进多目标线性规划模型

33l

运用数学优化软件LINGOll.0编程,分别求解费用、风险最小化优化问题和经TOPSIS方法转化后的单目标优化问题,在短时间内可得到最优解,其部分计算结果见表6.由此可见,本文设计的求解算法可得到该多目标0.1混合整数线性规划模型的满意解,还能求得费用和风险两个单目标优化问题的最优方案,决策者可依据具体情况选择满意的选址一路径方案.

不计相关中心的运营费用和选址风险,计算此

时费用、风险的单目标优化问题和多目标优化问题.将最优方案的结果(包括上述费用和风险)分别与本文模型对应方案的结果进行比较(表7).由表7可见,本文建立的模型更完善,单目标优化问

题结果使费用和风险值均有所下降,其单目标选

址一路径方案更优;对于多目标优化问题,本文的模型能够以较小的成本增长为代价较大幅度地降低风险,使多目标选址-路径方案更合理.

表5加工技术

Tab.5

Treatmenttechnology

注:“一”表示危险废物与加工技术不相容.

表6多目标优化问题部分计算结果

multi—objectiveoptimization

7I讪.6

Panof山eresultof

problem

Tab.7

表7计算结果对比

Resultscomp撕son

注:“+”表示增加,“一”表示减少.

结论

在本文的多目标模型中,费用目标函数考虑了

的实际要求,可以使得单目标和多目标优化问题的方案更优.对于中、小规模问题,可借助如IBM

ILOG

cPLEx、GAMS和LINGO等数学优化软件进

中心的运营费用,风险目标函数考虑了中心的选址风险.另外还考虑了危险废物类型和加工技术的多样性、危险废物和加工技术的相容性,并将其转化为已知参数和约束条件.考虑到废物产量的波动性,设定了考虑富余量的中心最大能力约束条件,进而建立了更完善的危险废物物流系统LRP问题多目标线性规划模型.采用极值法使费用和风险无量纲化,并利用TOPSIS法将多目标优化问题转换为单目标优化问题进行求解.最后,用算例表明了算法的有效性,证明了该模型更符合危险废物处理

行求解.

下一步研究可预测源点的危险废物产量,将风险公平性优化目标纳入目标函数综合考虑,并针对大规模的实际问题设计有效的求解方法.参考文献:

[1]

sHOBRYsD.mutes

modelfbrthe

locations

sekctionof

shipping

舳d

stomge

for

h吼rdou8

unive璐i哆,

substance[D].Baltimore:Johm

1981.

Hopkins

[2]

cAPPANERAP,GALLoG,MAFFlouF.Discrete

332

西南交通

学学报第46卷

矗lcility10c觚on蚰dmuting0f0bno)【iousactivm髓[J].

Di8cmteApplied

Ma出ematics,2004,133(1-3):3.28.

[3]S1DwERs

cL,PALEKARu

S.L们ati∞models

widl

muting

c∞side枷彻s

forasingleobno】【i伽sf.acil时[J].

Tr吼spon撕仰Scie眦e,1993,27(4):350_362.[4]

JAcOBsT

L。

wARMERDAM

simuh柚eou8muting蚰dsiting

f缸h删∞s-w鹪te

M.ope础ons[J].

Joumal0fUrb蚰Pl粕ning粕dDevelopm明t,1994,120(3):115—131.

[5]wYMANMM,KuBYM.Pm∞tive叩timizationof

t啦icw鹏ten锄sp0Itati∞,loc撕∞曲dtechnolo科[J].

L邯ati∞science。1995,3(3):167・185.[6]

EMEKE,KARABY.Ha髓rdousw聃temanagement

pmblem:t.Ie

c鹪efh

incine枷on[J].compute璐明d

Operati唧ReseaIch,2007,34(5):1424-1441.

[7]

usTG,MIRc|IANDANIP.Aninte胂ted弛啪rk

pl舳盯mllltiobjective

modelfbrmutingandsitirlgfor

h砌别伽s

m砒耐a18蚰d

wast鸱[J].Tr蚰sportationScience,199l,25(2):146-156.

[8]sHEuJ

B.A

cooIdinatedreve墙e

logistics

systemfor

regi∞al

m蚰agement

ofmulti一∞urceh北ard叫s

w鹋tes[J].compute璐&Opemtio璐R鹪earch,2007,34(5):1442・1462.[9]

ALUMuR

S,KARA

Y.

new

model

I缸theh髓ardo鸺w∞te

location-muting

pmblem[J].ComputeIs蚰d

0pemtions

Resear℃h,2007,34(5):

1406.1423.

[10]MIsRAV,PANDEYSD.H北aHlou8w鹪te,imp8ct

伽healttI锄d

envimnment

f撕developmentofbetter

w船tem蚰agementsh砒egi鹄infutu他inlndia[J].

Envimnment

Int咖ational,2005,3l(3):417_431.

[11]

NEMATHAGA

F,

MARINGA

S,

cHIMuKA

L.

Hospital∞lidw酗tem蚰agementpmcticesinLimpopo胁vince,SoutIIA付ca:

c嬲estudy0ftwo

h唧itals[J]-W鹪teM蚴gement,2008,28(7):

1236・1245.

[12]魏航,李军,蒲云.有害物品运输中避免重大事故发

生的路径选[J].西南交通大学学报,2006,4l(4):

49l一495.

WEIH粕g,UJun,PUY岫.Selection0fmut髓fbr

avoiding

seri伽s∞cidentsin

嘶sportation

of

h越ardousmaterials[J].

Jo岫al

ofsouthwest

Jiatong

Unive璐ity,2006,41(4):491-495.

[13]

王慈光.目标规划在重载运输组织中的应用[J]。西南交通大学学报,2009,44(3):392-395.

WANGCigll趴g.Applicationofg蚰lpr0肼啪mingin

谢1way

hea、r)r

haul

呦sportation[J].

J叫mal

of

S0utllwest

Ji咖ngUnivers时,2009,44(3):392-395.

[14]

GIANNIK0sI.Amultiobjective

p加学阻mmingmodelforlocating

t瑚tlllent

sites

and

muting

h北ardousw鹕tes[J].Eumpe舳Joumal0fop啪tional

Research,

1998,104(2):333-342.

[15]LAI

J,uu

Y,HwANG

CL.TOPsIS

for

MODM[J].

Eumpe锄

Jo哪al

of

0pemtional

R鸭earch,1994,76(3):486_500.

[16]

NEMAAK,Gu盯A

K.0ptimization

of陀gional

h砌:aldou8

waste

management

8ystems:

improved

follnul8tion[J].w鹊teM彻agement,1999,19(7-8):

44l-51.

[17]

NEMAA

K,Gum

K.Multiobjectiveriskanalysis

舳d

optiIIlization

ofre舀onal

h北ardousw酗te

m锄agement

system[J].

№tice

Periodical

of

H衄aIdous,Tbxic,

粕d

RadioactiveW鹕te

M肌agement,2003,7(2):69-77.

[18]

NEMA

AK,MODAK

PM.

stmtegic

desi即

appmach

for叩timization

0f

r嚼onal

h北ardousw鹊te

m帅agementsystems[J].

w璐te

M蚰agement蛐d

Re8earch,1998,16(3):210-224.

(中文编辑:秦萍玲英文编辑:刘斌)


© 2024 实用范文网 | 联系我们: webmaster# 6400.net.cn