深圳市南山区垃圾清运处理运输方案论文

1、深圳市南山区垃圾分类处理与清运方案设计摘 要随着国民经济发展与城市化进程加快，垃圾分类化收集与处理已经被提到日程上来了。针对深圳市南山区的垃圾，我们期望在保护环境的同时，也有利于资源的回收和再利用，进而达到经济效益的最大化。由题目中给定的诸多条件我们采取以下模型来解决提出的两个问题。模型一：采用模糊聚类的方法将上万个小区进行一次和二次简化归并，然后根据具体问题采用（0 1）规划的方法求出每一问中需要建设多少个大型厨余设备和多少个小型厨余设备。模型二：在模型一的基础上，使用多变量非线性规划，找出对经济最优和对环境最有益的方法，列出目标函数，找到约束条件，然后利用lingo软件进行非线性规划问题的

2、求解。 最后在模型改进中提出对问题考虑不周的地方，加入人为因素，使该模型对解决南山区的垃圾处理问题变得更加科学和贴近现实生活。关键词：厨余垃圾 多变量非线性规划 模糊聚类 （0 1）规划1 问题重述   垃圾分类化收集与处理有利于减少垃圾的产生，有益于环境保护，同时也有利于资源回收与再利用的城市绿色工程。在发达国家普遍实现了垃圾分类化，随着国民经济发展与城市化进程加快，我国大城市的垃圾分类化已经提到日程上来。在深圳，垃圾分为四类：橱余垃圾、可回收垃圾、有害垃圾和其他不可回收垃圾，四者比例约为4：2：1：3，所有垃圾将从小区运送到附近的转运站，再运送到少数几个垃圾处理中心。在

3、垃圾分类收集与处理中，不同类的垃圾有不同的处理方式，简述如下：（1）橱余垃圾 可以使用脱水干燥处理装置，处理后的干物质运送饲料加工厂做原料。不同处理规模的设备成本和运行成本（分大型和小型）。大型厨余垃圾处理设备（如南山餐厨垃圾综合利用项目，处理能力为200吨/日，投资额约为4500万元，运行成本为150元/吨。小型餐厨垃圾处理机，处理能力为200-300公斤/日，投资额约为28万元，运行成本为200元/吨。橱余垃圾处理后产物价格在1000-1500元/吨。（2）可回收垃圾 将收集后分类再利用。（3）有害垃圾 运送到固废处理中心集中处理。（4）其他不可回收垃圾 将运送到填埋场或焚烧场处理。

4、60; 本项研究课题旨在为深圳市的垃圾分类化进程作出贡献。为此我们将运用数学建模的方法对深圳市南山区的分类化垃圾的实现做一些研究，具体的研究目标是：（1）假定现有垃圾转运站规模与位置不变条件下，给出大、小型设备（橱余垃圾）的分布设计，同时在目前的运输装备条件下给出清运路线的具体方案。以期达到最佳经济效益和环保效果。（2）假设转运站允许重新设计，对于问题（1）的目标重新设计。2 问题分析由本题给出的条件和要求我们可以看出，对于深圳市南山区垃圾的处理与清运问题，我们要在保护环境即所有垃圾都能够被处理的前提下实现经济利益的最大化，而这里的经济效益最大化应该为总的成本最小，其中总成本=大小型厨余设备的

5、投资建设成本+垃圾处理成本+运输费用-垃圾处理后所得经济收益。题目中给出了南山区每日垃圾清运总量1280吨，其中804吨垃圾是需要我们考虑的，而剩下的已经直接拉到填埋场或焚烧厂了，按照所给比例，计算知道厨余垃圾为321.6吨，所以本题是在此基础上考虑大小型厨余设备的数量及分布情况，对于大小型设备厂的数量由于不同的设备厂的投资数额，处理能力，运转成本都有差异，我们依据建设成本最小可以大致的确定处理厂的数量，进而根据总成本最小建立多变量的非线性规划问题进行求解。对于第二问的中转站问题，我们首先把小区简化，依据位置和小区的人口数量进行决定中转站的数量和位置，剩下来的过程同第一问类似。 3 基本假设和

6、符号说明3.1 基本假设假设1： 题目所给数据及所指示的网站数据真实可靠。假设2：南山区的道路在短时间内不会有大的变动。假设3：所有厨余设备都能正常工作，不存在损坏或是休眠问题。假设4：厨余设备在各个地方建设时不存在建设成本的差异，即忽略选址点对建设成本的影响。假设5：点之间的道路畅通无阻，不考虑车辆空载的问题。假设6：由于题目中给出的小区数目过于庞大，为了便于操作，在不影响对问题分析正确性的前提下将小区数目做出简化，并将简化掉的小区人数合并在一起。假设7：南山区短时间内与资料做给人口不会有太大出入，没有大量人口的流入或者流出。假设8：各个车辆在运输的时候自动会选择最优路线，即耗油量最小的路线

7、，并在此基础上把实际地图简化为平面上的散点图，标出各个中转站的相对位置，以两点间的直线来代替车辆实际的行走路线。假设9：各小区的垃圾都是当天清理完毕，没有积累。假设10：由题目所给的804吨垃圾/南山区的总人数=人均日产垃圾量是稳定的。假设11： 不考虑其他附加因素对厨余设备建设地点的干扰。3.2符号说明（1）,m均为含有38项的列向量（2）各站需要运输的次数（3）取0或1，取1表示该垃圾站的垃圾运往第一垃圾处理中心，取0表示运往第二垃圾处理中心（4）各个垃圾站的橱余垃圾的数量（5）x1是第一个处理中心的经度（6）y1是第一个处理中心的纬度（7）x2是第二个处理中心的经度（8）y2是第二个处理

8、中心的纬度(9) （j=1，2，i=1，2，,38）表示第i个中转站的垃圾运往第j个处理中心的距离（10） ， （i=1，2，,38）表示第i个中转站的经度和纬度4 模型的建立与求解4.1 现有条件下厨余设备的数量求解和位置分布 大小型厨余设备数量的求解由题目中中所给的数据知道问题一中共有321.6吨厨余垃圾需要使用大小型厨余设备进行处理，如果建两个个大型厨余设备，需投资9000万元，可以处理完所有厨余垃圾，如果建一个大型处理设备，其余都建成小型的需要406608个小型设备，投资为1136817024万，则设备投资总额就为1586821524万元。由经纬度和物理知识我们可以算出，纬度相差1度时

9、实际距离相差111.7km，经度相差1度时实际距离相差117*cos（X），X表示所在的纬度，所以欧氏距离下得经纬度相差1度实际距离相差152.06km。由于南山区两点间的最远距离为15km，则运输成本最大为15/100*7.5*30*321.6=10854元，远小于上述两种方案的最低成本差6868万元，由此可得，建设两个大型厨余设备是比较合理的。 大型厨余设备位置的设计 由上述的求解我们知道建设两个大型处理设备是最合理的，那么这两个大型处理设备如何安排才能达到从各个中转站到这两个处理中心的路线最优呢？我们考虑如下：首先建立目标函数min，表示总的最小成本，因为考虑南山区的实际道路情况很困难，

10、所以假设司机都是知道选择最省油的路线的，我们用各转运站与厨余垃圾处理中心的直线距离来表示,所以建立目标函数如下：目标函数：约束条件：按照上述给定的目标函数和约束条件,我们利用lingo软件进行编程求解，结果如下:Variable Value Reduced CostX1 113.9652 0.000000Y1 22.56472 0.000000X2 113.9137 -2.116664Y2 22.51550 -3.807164T( 1) 0.000000 0.000000T( 2) 1.000000 0.000000T( 3) 1.000000 0.000000T( 4) 1.000000 0

11、.000000T( 5) 1.000000 0.000000T( 6) 0.000000 0.000000T( 7) 0.000000 0.000000T( 8) 0.000000 0.000000T( 9) 0.000000 0.000000T( 10) 0.000000 0.000000T( 11) 0.000000 0.000000T( 12) 0.000000 0.000000T( 13) 1.000000 0.000000T( 14) 0.000000 0.000000T( 15) 0.000000 0.000000T( 16) 1.000000 0.000000T( 17) 0.0

12、00000 0.000000T( 18) 0.000000 0.000000T( 19) 1.000000 0.000000T( 20) 0.000000 0.000000T( 21) 0.000000 0.000000T( 22) 0.000000 0.000000T( 23) 0.000000 0.000000T( 24) 0.000000 0.000000T( 25) 1.000000 0.000000T( 26) 1.000000 0.000000T( 27) 1.000000 0.000000T( 28) 1.000000 0.000000T( 29) 1.000000 0.0000

13、00T( 30) 0.000000 0.000000T( 31) 0.000000 0.000000T( 32) 1.000000 0.000000T( 33) 1.000000 0.000000T( 34) 1.000000 0.000000T( 35) 1.000000 0.000000T( 36) 0.000000 0.000000T( 37) 1.000000 0.000000T( 38) 1.000000 0.000000T（i）（i=1、2、38）表示第i个小区，其后所跟value值为1或0，1表示此中转站的垃圾运到第一个垃圾处理中心，0表示此中转站的垃圾运到第二个垃圾处理中心，这

14、样就将两个垃圾处理中心与38个中转站的相对位置找到了，以下为中转站的经纬度：表1 南山区垃圾中转站经纬度分布表序号中转站经度纬度序号中转站经度纬度1九街站113.917622.543820松坪山站113.913722.51552玉泉站113.928422.546221南光站113.925122.519393动物园站113.955322.572722南园站113.916622.51844平山村站113.96622.583623望海路站113.929722.48375牛城村站113.925922.616424花果路站113.930922.4886科技园站113.945722.518225福光站11

15、3.99822.59437同乐村站113.922522.567226新围村站113.954222.56858松坪山（二）站113.909822.519327大冲站113.956122.54069大新小学站113.92122.535728沙河市场站113.968822.541210南山村站113.90922.513329龙井113.97722.562811阳光站113.913722.515530南山市场113.912922.513812月亮湾大道站113.904222.514231麻勘站113.913722.515513光前站113.965222.564732白芒站113.939622.6212

16、14北头站113.907922.514133大石磡站113.974322.6173715涌下村站113.921122.532434长源村站114.010422.595816白石洲南站113.968522.533635华侨城站113.996822.537117前海公园站113.91122.545336疏港小区站113.900622.481218深圳大学站113.935422.533437西丽路站113.955322.572719官龙村站113.948622.59238塘朗站113.996322.59211根据上表一的垃圾中转站的相对位置以及计算求得的处理厂的经纬度，我们运用matlab绘图功能，

17、以经度为横坐标，纬度为纵坐标，画出具体的位置，如图1所示：注：1和2表示两个大型垃圾处理中心 图1 厨余设备安放位置示意图4.2 转运站允许重新设计时厨余设备的数量求解及分布第二问与第一问的不同点是垃圾转运站允许重新设计，我们借鉴第一问的部分结论进行推广，在转运站的位置可变的条件下，根据各个社区的具体位置来确定最佳的中转站分布情况。由于中转站的数量和位置的不确定性，使得这个问题的复杂程度上升。为此我们在处理整个问题之前，先把中转站的数量和位置决定下来。接下来决策厨余垃圾处理厂的建立与选址问题，同第一问一样，用类似的方法进行讨论解决。 各社区的初步分布 由于题目中给的小区数量过万，为了简化起见，

18、我们按照小区名字的相关度编程将上万个小区归并合整简化为如下的若干个小区，如表2所示：表2 初次归并后的小区简表编号小区名称人数经度纬度101网格2993113.96322.54672碧云片区1791113.9302122.5143常兴广场片区1289113.9231622.5409121中新街片区1318113.9616522.536122珠光村东区10270113.9616522.5701123珠光村西区5051113.9304422.533124珠光新屋村10244113.9596522.571运用matlab绘图功能，我们以这些小区的经度为横坐标，纬度为纵坐标，按照各个小区经纬度在图上标

19、示出其相对位置，从图中可以看出各小区的分布稠密情况，如图2所示： 图2 初次归并后各小区位置示意图 各社区的二次归并南山区每日垃圾清运总量1280吨，总人口为1320411人，分布在124个地区内（按照第一问的划分），所以可知人均日产厨余垃圾量为1280/1320411=0.9694 （kg）,按照人口密度，将所有小区重新分片，所得结果与初次分布的结果进行逐一对应，如下表3所示：表3 二次归并后新旧小区对照示意图及相关数据分布表新新编号人人数经经度纬度度垃垃圾总量（吨）厨余垃圾（吨）原始编号118150113.931622.5379517.594617.03784436,38,45,46,55

20、,68，215667113.947122.545215.18758986.07503644,69，70，72，73，75，87，89，91，106，107，111345338113.958322.57143.950657217.5802624,28,41,42,50,58,76,80,121,123485598113.969122.5607782.978701233.1914835,53,56,59,81,83,96,99,102537776113.928322.5309836.620054414.6480257,60,77,78,92,97,100,105,1136185108113.925

21、322.49034179.44369571.777487,11,48,61,63,64,82,11779034113.929222.51358.75755963.50302443,62,66,67,93,108,114,115814038113.966422.5433313.60843725.4433751,4,6,8,10,12959576113.959422.587257.752974423.1011927,110,109,49,881022479113.979922.5378221.79114268.71645717,52,79,118,120113624113.943622.54183

22、.51310561.40524247,65,71,74,951242679113.907122.5151641.373022616.5492116,21,23,29,3013235479113.993722.53845228.27334391.309342,14,15,119143130113.987922.57443.0342221.213689103,54153281113.969122.601753.18060141.27224194,5116162900113.933622.4972157.9152663.1661317128946113.968522.5336125.00025250

23、.000151886836113.939722.485584.178818433.6715391977394113.986422.549175.025743630.0103132011743113.887622.508611.38366424.553466182118875113.900622.546218.2974257.31897202218803113.987922.589218.22762827.29105198235619114.010622.59635.44705862.178823116245078114.027322.55284.92261321.969045124251266

24、113.914922.76421.22726040.490904125265636113.895722.47325.46353842.18541526271106113.941422.554731.07215640.04288637283544113.949222.52683.43555360.13742232由上表可知，在处理第二问时将所有小区划分为28个片区，在这28个区域内我们假设放置的有足够接收这一片区的垃圾中转站，这样，问题就转化为在这28个片区的基础上考虑需要多少个垃圾中转站和垃圾厨余设备，以及怎样放置这些厨余设备。厨余垃圾一共有1280*0.4=512吨，按照第一问关于厨余设备数

25、量的求解过程（方法类似，这里不再赘述），我们知道还是全部修建成大型厨余设备是最省成本的，这样就需要三个大型设备即可以满足需求。 大型厨余设备位置的求解由得知需建三个大型厨余设备，设三个大型厨余设备的编号分别为1，2，3，根据题设条件建立目标函数目标函数：约束条件：利用lingo进行求解，得出结果如下：结果：Local optimal solution found.Objective value: 2.623194Objective bound: 2.623194Infeasibilities: 0.000000Extended solver steps: 0Total solver itera

26、tions: 22Variable Value Reduced CostX1 113.9619 -0.7683643Y1 22.53396 -1.206360X2 113.9619 -1.054813Y2 22.53396 -0.6262339X3 113.9619 0.000000Y3 22.53396 0.000000T1( 18) 1.000000 0.000000T1( 23) 1.000000 0.000000T1( 24) 1.000000 0.000000T2( 20) 1.000000 0.000000T2( 22) 1.000000 0.000000T2( 25) 1.000

27、000 0.000000T3( 1) 1.000000 0.000000T3( 2) 1.000000 0.000000T3( 3) 1.000000 0.000000T3( 4) 1.000000 0.000000T3( 5) 1.000000 0.000000T3( 6) 1.000000 0.000000T3( 7) 1.000000 0.000000T3( 8) 1.000000 0.000000T3( 9) 1.000000 0.000000T3( 10) 1.000000 0.000000T3( 11) 1.000000 0.000000T3( 12) 1.000000 0.000

28、000T3( 13) 1.000000 0.000000T3( 14) 1.000000 0.000000T3( 15) 1.000000 0.000000T3( 16) 1.000000 0.000000T3( 17) 1.000000 0.000000T3( 19) 1.000000 0.000000T3( 21) 1.000000 0.000000其中T1（*）、T（*）、T（*）分别表示三个垃圾厨余设备，括号里面的数字表示上述28个片区运到该厨余设备所在处进行处理，也可将结果写成如下形式，便于观看，同时将某个地区全部建成小型厨余设备所需要的数量也一并写出，便于我们对比，从而也可看出，建

29、一个大型设备比建若干个小型设备要省去很多成本表4 各个片区与运往的厨余中心对照表及大小型厨余设备的对照表新编号橱余垃圾清运次数经度纬度大型厨余设备的编号当前编号小型个数123270.0428860280.1374220250.4909040141.2136891113.987922.5744114151.2722411113.969122.60175125111.4052421113.943622.5418135241.9690451114.027322.5528147232.1788231114.010622.5963158262.1854151113.895722.473216873.50

30、30241113.929222.51351712204.5534661113.887622.5086181585.4433751113.966422.54333191826.0750361113.947122.54521102117.0378441113.931622.5379511124227.2910511113.987922.589211225217.318971113.900622.546211325108.7164571113.979922.5378211429514.648022113.928322.53098115491216.549212113.907122.515161165

31、5317.580262113.958322.57111759923.101193113.959422.5872118771930.01033113.986422.5491119100433.191484113.969122.560771201101833.671534113.939722.48551211121750.00015113.968522.53361221161663.16617113.933622.4972123210671.777488113.925322.490341242391391.309349113.993722.53845125304上图中标注为红色的25、27、28因

32、为所负担的厨余垃圾很少，在这里我们不纳入考虑的范围，将在模型改进中进行说明。下图是用来标示三个垃圾厨余设备（分别为1、2、3）与各个片区之间的相对位置关系，横坐标为纬度，纵坐标为经度。图3 重新分布中转站后大型厨余设备的位置分布图5 模型改进 在中我们提到有三个片区产生的厨余垃圾相对比较少，我们在上述建模求解的过程中将其忽略不计，可是现实生活中还是要把着三个地区考虑进来，这样我们在不过多增加成本的前提下，考虑建设若干个小型厨余设备，用来满足这些地区的垃圾回收和处理，建立如下所示的目标函数：目标函数： 另外将每天的汽车运营成本考虑进来，是模型考虑的方面更全，同时，我们还要考虑的一个因素就是时间问

33、题，若在个别地区多建一些小型的垃圾处理中心，可能会增加前期投资，但是如果长久来看，加上油费，车辆损耗费，司机工资等诸多方面的因素，在不改变大型数量的前提下，适当多建一些小型的厨余设备，所耗费的总成本肯定是会降低的。在求解目标函数的过程中如果约束条件或是目标函数取值不当，会令目标函数陷入局部最优解或者无解的境地，这也是我们需要避免的地方。6 参考文献1张光澄，非线性规划最优解计算方法，高等教育出版社，2005.2贾传兴，彭绪亚，刘国涛，刘长玮，伍翔，邓稼佳，城市垃圾中转站选址优化模型的建立及其应用，环境科学学报，第26卷第11期.3刘桐武，刘兆龙。线性规划在城市垃圾运输中的应用J。环境卫生工程，

34、1996（2）.4闵仲求，李毅华。单目标和多目标系统线性规划M，2004,23(1).5陈宝林，最优化理论与算法M,北京：清华大学出版社，1989,49.7 附录程序1：model:sets:n/1.25/:a,t1,t2,t3,m,xx,yy;endsetsmin=sum(n(i):sqrt(x1-xx(i)2+(y1-yy(i)2)*t1(i)*a(i)+sqrt(x2-xx(i)2+(y2-yy(i)2)*t2(i)*a(i)+sqrt(x3-xx(i)2+(y3-yy(i)2)*t3(i)*a(i);sum(n(i):m(i)*t1(i)<200;sum(n(i):m(i)*t2

35、(i)<200;sum(n(i):m(i)*t3(i)<200;for(n(i):t1(i)+t2(i)+t3(i)=1);for(n:bin(t1);for(n:bin(t2);for(n:bin(t3);DATA:a=1111111111111122233445789;xx=113.98789113.969135113.94356114.02729114.01056113.89574113.9291888113.88759113.9664383113.9470542113.9316367113.98789113.90056113.979922113.9283189113.907

36、082113.958302113.95938113.98637113.9690833113.93966113.9685113.93357113.9252614113.9937417;yy=22.574422.6017522.541822.552822.596322.473222.513522.508622.5433333322.545222.5379522.589222.546222.5378222.5309777822.5151622.57122.587222.549122.5607666722.485522.533622.497222.4903357122.53845;m=1.213688

37、81.272240561.405242241.969045282.178823442.185415363.503023844.553465685.443374886.075035927.0378447.291051287.318978.7164570414.6480217616.5492090417.5802628823.1011897630.0102974433.1914804833.6715273650.0001009663.16610471.7774780891.30933704;ENDDATAEnd程序2：model:sets:n/1.38/:a,t,m,xx,yy;endsetsmi

38、n=sum(n(i):sqrt(x1-xx(i)2+(y1-yy(i)2)*t(i)*a(i)+sqrt(x2-xx(i)2+(y2-yy(i)2)*(1-t(i)*a(i);sum(n(i):m(i)*t(i)<200;sum(n(i):m(i)*(1-t(i)<200;for(n:bin(t);DATA:a=1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 2 1 1 1 1 1 1 2 2 1 1 2 2 1 1 1 1 2 1 3 2 1 1;xx=113.91758 113.92843 113.95528 113.96601 113.92587 113.9456

39、6 113.92249 113.90982 113.92097 113.90896 113.91374 113.90416 113.9652 113.90792 113.92107 113.9685 113.91098 113.93543 113.94863 113.91374 113.9251 113.91655 113.92968 113.93089 113.99795 113.95423 113.95613 113.96884 113.97695 113.91288 113.91374 113.93958 113.9743 114.01039 113.99684 113.90056 11

40、3.95528 113.9963;yy=22.5438 22.5462 22.5727 22.5836 22.6164 22.5182 22.5672 22.5193 22.5357 22.5133 22.5155 22.5142 22.5647 22.5141 22.5324 22.5336 22.5453 22.5334 22.592 22.5155 22.519388 22.5184 22.4837 22.488 22.5943 22.5685 22.5406 22.5412 22.5628 22.5138 22.5155 22.6212 22.61737 22.5958 22.5371

41、 22.4812 22.5727 22.59211;m=8 10 8 10 2 8 2 4 12 10 4 16 8 6 8 12 6.4 6 6 10 6 6 12 12 4 8 14 12 6 10 4 3.2 12 2 28 16 6 4;ENDDATAend程序3：model:sets:n/1.19/:a,t1,t2,m,mm,xx,yy;endsetsmin=sum(n(i):sqrt(x-xx(i)2+(y-yy(i)2)*t1(i)*4.3*152+t1(i)*m(i)*150+t2(i)*m(i)*200+t2(i)*mm(i)*76.7);for(n(i):bin(t1);f

42、or(n(i):bin(t2);for(n(i):t1(i)+t2(i)=1);BND(112,x,115);BND(21,y,24); DATA:a=1111111111111122234;xx=113.98789113.969135113.94356114.02729114.01056113.89574113.9291888113.88759113.9664383113.9470542113.9316367113.98789113.90056113.979922113.9283189113.907082113.958302113.98637113.93966;yy=22.574422.60

43、17522.541822.552822.596322.473222.513522.508622.5433333322.545222.5379522.589222.546222.5378222.5309777822.5151622.57122.549122.4855;m=1.21368881.272240561.405242241.969045282.178823442.185415363.503023844.553465685.443374886.075035927.0378447.291051287.318978.7164570414.6480217616.5492090417.580262

44、8830.0102974433.67152736;mm=4557881215182124252529495559100112;ENDDATAEnd程序4：model:sets:n/1.19/:a,t1,t2,m,mm,xx,yy;endsetsmin=sum(n(i):4.3*a(i)*sum(n(i):sqrt(x1-xx(i)2+(y1-yy(i)2)*t1(i)+t1(i)*m(i)*150+t2(i)*m(i)*200+t2(i)*mm(i)*76.7);for(n(i):t1(i)+t2(i)=1);for(n:bin(t1);for(n:bin(t2);DATA:a=1111111

45、111111122234;xx=113.98789113.969135113.94356114.02729114.01056113.89574113.9291888113.88759113.9664383113.9470542113.9316367113.98789113.90056113.979922113.9283189113.907082113.958302113.98637113.93966;yy=22.574422.6017522.541822.552822.596322.473222.513522.508622.5433333322.545222.5379522.589222.546222.5378222.5309777822.5151622.57122.549122.4855;m=1.21368881.272240561.405242241.969045282.178823442.185415363.5