深圳市垃圾处理问题建模论文

1/1深圳市垃圾处理问题建模论文

垃圾分类处理与清运方案设计

圾，其中对于居民垃圾，基本的分类处理流程如下：

在垃圾分类收集与处理中，不同类的垃圾有不同的处理方式，简述如下：

1)橱余垃圾。可以使用脱水干燥处理装置，处理后的干物质运送饲料加工厂做原料。不同处理规模的设备成本和运行成本（分大型和小型）见附录1说明。

2)回收垃圾。将收集后分类再利用。3)有害垃圾。运送到固废处理中心集中处理。

4)其他不行回收垃圾。将运送到填埋场或焚烧场处理。

全部垃圾将从小区运送到四周的转运站，再运送到少数几个垃圾处理中心。明显，1）和2）两项中，经过处理，回收和利用，产生经济效益，而3）和4）只有消耗处理费用，不产生经济效益。

本项讨论课题旨在为深圳市的垃圾分类化进程作出贡献。为此请你们运用数学建模方法对深圳市南山区的分类化垃圾的实现做一些讨论，详细的讨论目标是：

1）、假定现有垃圾转运站规模与位置不变条件下，给出大、小型设备（橱余垃圾）的分布设计，同时在目前的运输装备条件下给出清运路线的详细方案。以期达到最佳经济效益和环保效果。

2）假设转运站允许重新设计，请为问题1）的目标重新设计。

二问题分析

对于问题一：从以下四点进行讨论。

一、垃圾中转站的位置与数量已定，但厨余垃圾处理中心的位置与数量不确定。

二、为了确定厨余垃圾处理中心的数量与位置，我们首先依据柴油费做不不超过的数量确定大型垃圾处理中心的数量。

三、依据大型垃圾处理中心的数量进行分区（本文分三区，即大型垃圾处理中心的数量为三台）

四，依据步骤三，将各个区的垃圾转运站进行比较，求找使总成本最小的

坐标位置，确立大型垃圾处理中心的坐标，相应的确立小型垃圾处理中心的位置。由此可规划处南山区垃圾分类处理与清运方案设计。

对于问题二：从以下几点进行讨论。

依据问题一的求解可知，建设小型垃圾处理中心的位置在偏远地区，并且产生的垃圾总量比较小，只有同时满意才适合建设，否则运往大型垃圾处理中心比较经济。因此可以将距离比较近的垃圾转运站合并，并建立大型垃圾处理中心，在偏远的地方的并且产生垃圾量比较少的点建立小型垃圾处理中心，偏远但产生垃圾量比较多的点保留垃圾转运站，大型垃圾处理中心进行垃圾处理。关于本方案的已知条件陈述：

只拖十吨的大型厢，只用于从转运站到垃圾中心，每次只拖一个大型厢，平均每公里耗油25L—30L柴油/每公里。柴油单价：6.39元/升。

收集车辆：

只负责从小区的垃圾站到转运站运输，60辆2.5吨汽车，每车耗油20L—35L70#汽油/百公里。

司机月薪平均3500元。

三模型假设

（1）假设南山区每个人产生的垃圾数量肯定；垃圾增长量正比例与人口增长量。（2）大型垃圾处理设备和小型垃圾处理设备的使用年限为15年，并且15年内花费的维护费用，工人工资等包含在每吨的处理费用之中，不需要另外计算。（3）假设柴油费7.51元/L，汽油费元/L，为简化计算，便于建模，15年内中国局势稳定，物价稳定，暂不增长。

（4）假设垃圾分类已在转运站已经装好箱，并且从垃圾转运站到垃圾处理中心，每天只运输一次垃圾，运输次数等于转运站垃圾箱数，运输吨数为垃圾总量。假定：司机每天工作时间为t=8小时，每月工作25天，工资由同一银行统一结算。并且，每日花费的处理费用，垃圾回收等经济收入等等都由银行系统统一结算。

（5）假设一个小型厨余垃圾处理中心可以同时包含多个小型厨余设备，即可以看成多个小型厨余设备的组合

（6）假设拖车的平均速度为20km/h，汽车的平均速度为30km/h

（7）假设每个垃圾转运站有可以收集本区垃圾，运输本区垃圾的小型车和拖车数目。

（8）假设小型车司机和拖车司机工作性质一样，工作可以互换。假定：每个垃圾转运站垃圾处理只有两种选择1）转送到垃圾处理中心，2）在本地建设小型垃圾处理中心。

四符号说明

（1）全局变量（详细程序实现见附录程序1）

（2）垃圾处理设备的使用年份从2023年至2026年，分别用年数1、2、3、…15表示；

（3）各年的人口增长量ri（i=1,2,3…15）表示各年的人口增长率。（4）分区自上至下分别为一区，二区，三区，花费的总费用分别为F、F2、F3；（5）相应的剔除建设小型垃圾垃圾处理中心的垃圾转运站到大型垃圾处理中心

的位置的最小距离和为f，f2，f3。

（6）任意一个垃圾转运站的第一年产生的垃圾总量为Wi，i=1、2、3…38（7）任意一个垃圾转运站第一年产生的厨余垃圾总量为wi，i=1、2、3…38（8）任意一个垃圾转运站十五年产生的垃圾总量的平均量为Wn，i=1、2、3…38（9）任意一个垃圾转运站十五年产生的厨余垃圾的平均量为wn，i=1、2、3…38分别为W，W2,W3

（10）各区大型垃圾处理中心的坐标x1x2

（11）各区小型垃圾处理中心的坐标atbt，t=1,2…38；特殊说明：当MATLAB消失相同的小型垃圾处理中心是表示剔除的是同一个点，即：小型垃圾处理中心的的选取采纳集合中的并集处理方法，得到每个区的坐标的并集即剔除点（小型垃圾处理中心的坐标）

(12)箱数ki，即表示垃圾处理中心的垃圾箱数，又表示垃圾处理中心运往大型垃圾处理中心的运输次数

五模型的建立与求解

5.1数据的初步分析和处理：

5.1.1运费单位柴油取27.5L/百公里.吨，即0.275L/（公里.吨），2.06525元/（公里.吨）

5.1.2坐标比例尺0.617公里/厘米。

5.1.3确定大型垃圾处理中心的数目，15年运费不大于139W-4500万元，每日运费不大于257.407407W-8333.33元（运费不高于20920.06225公

里.吨）,W不超过200吨，建设大型垃圾处理中心是合理的，由此可推算出大约建设3个大型垃圾处理中心，依据深圳南山区地形，自上而下分3区，依据坐标图上帮助标号，考虑城市将来的进展规划问题，深圳北部地区人口密集度比较低，地区土地闲置率高，下一步，城市将加大城市北部的进展，故人口增长的速度相对较高，所以大型垃圾处理设备目前的负荷量比较低，但考虑到长远的进展是合适的，分区状况1-15点位一区，16-25区为二区，26-38区为三区。

5.1.4箱数与吨数问题，由于柴油的消耗既与路程有关又有运输的重量有关，每次运输都是双向的即拖车每次托完垃圾要回到起始位置例如，假如从垃圾处理中心去转运站拖垃圾，要先空车去转运站，然后拖着垃圾回来；或者从垃圾处理中心拖着箱子去处理中心送完垃圾拖空箱会转运站。

5.1.5每个垃圾转运站每天花费的处理费用，运往大型垃圾处理中心花费的费用free1=wn（大型机摊派的处理费用+运费），自建小型垃圾处理中心的费用free2=wn（小型机摊派费用），即：

free1=(45000000+150*W1*15*360)*wn/W1/15/360+wn*7/4*2.06525*0.610*

（|x(1)-a(t)|+|(x(2)-b(t)|);

free2=（280000*【wn/0.2】+200*wn*15*360)/15/360;注：（【】表示高斯取整）

各转运15年运费总和即15年中每一天运输垃圾的费用与处理垃圾的费用和。

1)不建立小型垃圾处理中心时F=360*15*F1=free1/wn*W*360*152）建立小型垃圾处理中心时

free2

F=15*360*F1=(+(45000000+150*(W-wn)*15*360)/15/360)*15*360；5.1.6关于深圳人口增长率问题

使用灰色模型求解即可，MATLAB程序实现程序见后。处理结果：5.1.7司机工资问题：

如上，计算出总路程，除以车速即可得出司机每日的总工作时间m司机每月工作25天3500元/8*25小时=17.5元/小时

5.2问题一：

(一)计算某区的总费用，以及大型和小型垃圾处理中心位置A．步骤，计算只建立任意一个小型垃圾处理中心的状况。

（1）求解某个区的某个点（大型垃圾处理中心位置）使其到全部点（垃圾转运

站的位置，出去一个点）的运费最小。

i

（2）求解点自己建设小型垃圾处理中心十五年的花费费用总和平均到每天，用free2表示

i

（3）求解点运往大型垃圾处理中心花费的费用之和平均到每一天，用free1表示

i

比较free1与free2，假如free1小，则点的垃圾也运往垃圾处理中心处理，计算该区总费用用

i

F

，否则，自己建设小型垃圾处理中心，计算该区花费的总费

F

（4）比较全部的状况

F

，选择最小

F

的的状况，即建一个小型垃圾处理中心的最优

F

。并用MATLAB程序返回大小型垃圾处理中心的坐标点以及

F

i

B．步骤，计算建立任意个小型垃圾处理中心的状况

（5）求解某个区的某个点（大型垃圾处理中心位置）使其到全部点（垃圾转运

i

站的位置，出去个点

i1，i2…

）的运费最小。

i

（6）求解点每一个自己建设小型垃圾处理中心十五年的花费费用总和平均到每天，用free2表示

i

（7）求解点每一个运往大型垃圾处理中心花费的费用之和平均到每一天，用

free1表示

（8）比较free1与free2，假如free1小，则点的垃圾也运往垃圾处理中心处理，计算该区总费用总费用

i

F

，否则，自己建设小型垃圾处理中心，计算该区花费的

F

（9）比较全部的优状况

F

，选择最小

F

的的状况，即建个小型垃圾处理中心的最

x1x2

i

F

。并用MATLAB程序返回大小型垃圾处理中心的坐标点以及

atbt

F

注：详细的程序实现见附录程序2,3

同理二区，三区的解法如上；程序实现分别见附录程序4，5,6,7

5.3问题二：

依据问题一的求解可知，建设小型垃圾处理中心的位置在偏远地区，并且产生的垃圾总量比较小，只有同时满意才适合建设，否则运往大型垃圾处理中心比较经济。因此可以将距离比较近的垃圾转运站合并，并建立大型垃圾处理中心，在偏远的地方的并且产生垃圾量比较少的点建立小型垃圾处理中心，偏远但产生垃圾量比较多的点保留垃圾转运站，大型垃圾处理中心进行垃圾处理。

模型求解结果：

图1垃圾转运站的坐标图

图2人口增长率图

图3各转运站厨余垃圾分布直方图

二．MATLAB求解结果一区1-15号点

大型垃圾处理中心坐标[14.6726.15]小型垃圾处理中心的坐标5,14号点f=80.87cm

花费的总费用F=1.51E+08二区16-25号点

大型垃圾处理中心坐标[1218]小型垃圾处理中心的坐标不建设f=48.88cm

花费的总费用F=2.95E+08

三区26-38号点

大型垃圾处理中心坐标[913.82]小型垃圾处理中心的坐标不建设f=57.38cm

花费的总费用F=3.33E+08

三个区花费的十五年总费用为7.80E+08平均每天的花费为144444.4元

拖车司机每天的总工作时间t总=（80.87+48.88+57.38）*0.617/20*17.5=99.88064元

故可以求知每天处理垃圾的总费用为144544元（四舍五入）

垃圾回收经济收入

一区二区三区

本部分直接由上部计算结果，由excel进行数据统计

由此可以看出，建设小型垃圾处理中心的位置在偏远地区，并且产生的垃圾总量比较小，只有同时满意才适合建设，否则运往大型垃圾处理中心比较经济。

该方案花费的费用最小，产生的垃圾全部处理完毕，故符合经济和环保的目的，方案是可行的。但是，鉴于时间少和我们自身的阅历的不足，以及理论与现实的差距，该方案也有其不足之处，我们不能保证他是肯定最优的，但却是相对最优的。

三．问题二求解

附录：

程序1

globalabWw

a=[14.67,12.08,12.6,18.19,10.32,23.98,21.55,16.2,21.78,13.89,16.02,14.3,13.63,9.45,16.65,18.5,12.14,11.9,7.38,11.12,7.56,17.39,8.84,14.9,21.65,8.37,12.48,8.65,16.86,9,8.21,13.25,8.67,9.45,7.08,10.8,5.19,10.92];

b=[33.3,32.67,31.45,30.4,29.92,26.61,26.3,26.15,25.9,25.37,24.85,24.65,22.61,21.55,20.95,20.55,20.12,20,18.65,18.12,17.92,17.57,17.26,17,16.59,16.42,15.91,15.45,15.27,14.54,13.9,13.82,13.52,13.37,12.55,7.75,7.5,7.01];

Wi=[101083055102010152520155201525104025163020357030152030251520151525304030];

wi=[5.715.714.5717.142.862.865.7111.435.718.5714.2911.438.572.8611.438.5714.295.7122.8614.299.1417.1411.4320.0040.0017.148.5711.4317.1414.298.5711.438.578.5714.2917.1422.8617.14];

wm=[6.16616.16614.932918.49833.08312.332215.41529.249112.33226.11616.11613.0839.24916.166112.33229.249115.41526.166124.66449.865715.415212.332218.498312.332215.41529.24919.24919.249115.41529.249112.332218.498318.498343.162621.581324.664418.498318.4983];

W=[10.790710.79078.632632.3725.395321.581326.976616.18621.581310.790710.79075.395316.1865.395321.581316.18626.976610.7907

43.162617.265126.976621.581332.37221.581326.976616.18616.18616.18626.976616.18621.581332.37232.37275.734637.767343.162632.37232.372];

r=[11.680012.353911.885011.433911.000010.582510.18099.79459.42289.06518.72118.39018.07177.76547.4706];

附录程序2非线性规划求解一区最小路费点坐标，即大型垃圾处理中心坐标。functionf=zuiyou(x)globalabTwmf=0;i=1;t=1;fort=1:15fori=1:15f1=0;f2=0;f=0;fori=1:(t-1)

s(i)=abs(x(1)-a(i))+abs(x(2)-b(i));f1=f1+s(i);end

fori=(t+1):15

s(i)=abs(x(1)-a(i))+abs(x(2)-b(i));f2=f2+s(i);end

[x(1),x(2)];[a(t),b(t)];f=f1+f2;endend

附录程序3输出一区大型垃圾处理中心坐标以及小型垃圾处理中心坐标和花费

的最小总费用F，最短路程和fclc

F=100000000000000000000000000000000000;fort=1:15[a(t),b(t)]fort=1:15[a(t),b(t)]fort=1:15[a(t),b(t)]fort=1:15[a(t),b(t)]fort=1:15[a(t),b(t)]fort=1:15[a(t),b(t)]fort=1:15[a(t),b(t)]fort=1:15[a(t),b(t)]fort=1:15[a(t),b(t)]fort=1:15[a(t),b(t)]fort=1:15[a(t),b(t)]fort=1:15[a(t),b(t)]fort=1:15[a(t),b(t)]fort=1:15[a(t),b(t)]x0=[17;26];vlb=[9,20];vub=[24,34];

x=fmincon('zuiyou',x0,,,,,vlb,vub);W1=118.8571429;

wn=wm(t);

free1=(45000000+150*W1*15*360)*wn/W1/15/360+wn*7/4*2.06525*0.610*(abs(x(1)-a(t))+abs(x(2)-b(t)));

free2=(280000*(0.5+wn/0.2)+200*wn*15*360)/15/360;if(free1free2)

[x(1),x(2)];[a(t),b(t)];f=zuiyou(x);

disp('建设小型垃圾处理中心')

F1=free2+(45000000+150*(W-wn)*15*360)/15/360;if(F1F)

[x(1),x(2)][a(t),b(t)]f=zuiyou(x)F=F1end

else(free1free2)[x(1),x(2)];

disp('不建设小型处理中心')F1=free1/wn*W;if(F1F)[x(1),x(2)]f=zuiyou(x)F=F1endendendendendendendendendend

endendendendend

附录程序4非线性规划求解二区最小路费点坐标，即大型垃圾处理中心坐标。functionf=zuiyou(x)globalabTwmf=0;i=1;t=1;fort=16：25fori=16:25f1=0;f2=0;f=0;fori=1:(t-1)

s(i)=abs(x(1)-a(i))+abs(x(2)-b(i));f1=f1+s(i);end

fori=(t+1):15

s(i)=abs(x(1)-a(i))+abs(x(2)-b(i));f2=f2+s(i);end

[x(1),x(2)];[a(t),b(t)];f=f1+f2;endend

附录程序5输出一区大型垃圾处理中心坐标以及小型垃圾处理中心坐标和花费的最小总费用F，最短路程和f

F=100000000000000000000000000000000000;fort=16:25[a(t),b(t)]fort=16:25[a(t),b(t)]fort=16:25[a(t),b(t)]fort=16:25[a(t),b(t)]fort=16:25[a(t),b(t)]fort=16:25[a(t),b(t)]fort=16:25[a(t),b(t)]fort=16:25[a(t),b(t)]x0=[10;19];vlb=[5,25];vub=[16,22];

x=fmincon('zuiyou',x0,,,,,vlb,vub);W1=163.4285714;

wn=wm(t);

free1=(45000000+150*W1*15*360)*wn/W1/15/360+wn*7/4*2.06525*0.610*(abs(x(1)-a(t))+abs(x(2)-b(t)));

free2=(280000*(0.5+wn/0.2)+200*wn*15*360)/15/360;if(free1free2)[x(1),x(2)];[a(t),b(t)];f=zuiyou(x);

disp('建设小型垃圾处理中心')

F1=free2+(45000000+150*(W-wn)*15*360)/15/360;if(F1F)

[x(1),x(2)]

[a(t),b(t)]f=zuiyou(x)F=F1end

else(free1free2)[x(1),x(2)];

disp('不建设小型处理中心')F1=free1/wn*W;if(F1F)[x(1),x(2)]f=zuiyou(x)F=F1endendendendendendendendendend

附录程序6非线性规划求解二区最小路费点坐标，即大型垃圾处理中心坐标。functionf=zuiyou(x)globalabTwmf=0;i=1;t=1;fort=16:25fori=16:25f1=0;f2=0;f=0;

fori=1:(t-1)

s(i)=abs(x(1)-a(i))+abs(x(2)-b(i));f1=f1+s(i);end

fori=(t+1):15

s(i)=abs(x(1)-a(i))+abs(x(2)-b(i));f2=f2+s(i);end

[x(1),x(2)];[a(t),b(t)];f=f1+f2;endend

附录程序7输出三区大型垃圾处理中心坐标以及小型垃圾处理中心