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

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

《深圳市南山区垃圾清运处理运输方案论文.doc》由会员分享，可在线阅读，更多相关《深圳市南山区垃圾清运处理运输方案论文.doc（27页珍藏版）》请在冰点文库上搜索。

深圳市南山区垃圾分类深圳市南山区垃圾分类处理处理与清运方案设计与清运方案设计摘摘要要随着国民经济发展与城市化进程加快，垃圾分类化收集与处理已经被提到日程上来了。

针对深圳市南山区的垃圾，我们期望在保护环境的同时，也有利于资源的回收和再利用，进而达到经济效益的最大化。

由题目中给定的诸多条件我们采取以下模型来解决提出的两个问题。

模型一：

采用模糊聚类的方法将上万个小区进行一次和二次简化归并，然后根据具体问题采用（01）规划的方法求出每一问中需要建设多少个大型厨余设备和多少个小型厨余设备。

模型二：

在模型一的基础上，使用多变量非线性规划，找出对经济最优和对环境最有益的方法，列出目标函数，找到约束条件，然后利用lingo软件进行非线性规划问题的求解。

最后在模型改进中提出对问题考虑不周的地方，加入人为因素，使该模型对解决南山区的垃圾处理问题变得更加科学和贴近现实生活。

关键词：

关键词：

厨余垃圾多变量非线性规划模糊聚类（01）规划11问题重述问题重述垃圾分类化收集与处理有利于减少垃圾的产生，有益于环境保护，同时也有利于资源回收与再利用的城市绿色工程。

在发达国家普遍实现了垃圾分类化，随着国民经济发展与城市化进程加快，我国大城市的垃圾分类化已经提到日程上来。

在深圳，垃圾分为四类：

橱余垃圾、可回收垃圾、有害垃圾和其他不可回收垃圾，四者比例约为4：

2：

1：

3，所有垃圾将从小区运送到附近的转运站，再运送到少数几个垃圾处理中心。

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

（1）橱余垃圾可以使用脱水干燥处理装置，处理后的干物质运送饲料加工厂做原料。

不同处理规模的设备成本和运行成本（分大型和小型）。

大型厨余垃圾处理设备（如南山餐厨垃圾综合利用项目，处理能力为200吨/日，投资额约为4500万元，运行成本为150元/吨。

小型餐厨垃圾处理机，处理能力为200-300公斤/日，投资额约为28万元，运行成本为200元/吨。

橱余垃圾处理后产物价格在1000-1500元/吨。

（2）可回收垃圾将收集后分类再利用。

（3）有害垃圾运送到固废处理中心集中处理。

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

本项研究课题旨在为深圳市的垃圾分类化进程作出贡献。

为此我们将运用数学建模的方法对深圳市南山区的分类化垃圾的实现做一些研究，具体的研究目标是：

（1）假定现有垃圾转运站规模与位置不变条件下，给出大、小型设备（橱余垃圾）的分布设计，同时在目前的运输装备条件下给出清运路线的具体方案。

以期达到最佳经济效益和环保效果。

（2）假设转运站允许重新设计，对于问题

（1）的目标重新设计。

22问题分析问题分析由本题给出的条件和要求我们可以看出，对于深圳市南山区垃圾的处理与清运问题，我们要在保护环境即所有垃圾都能够被处理的前提下实现经济利益的最大化，而这里的经济效益最大化应该为总的成本最小，其中总成本=大小型厨余设备的投资建设成本+垃圾处理成本+运输费用-垃圾处理后所得经济收益。

题目中给出了南山区每日垃圾清运总量1280吨，其中804吨垃圾是需要我们考虑的，而剩下的已经直接拉到填埋场或焚烧厂了，按照所给比例，计算知道厨余垃圾为321.6吨，所以本题是在此基础上考虑大小型厨余设备的数量及分布情况，对于大小型设备厂的数量由于不同的设备厂的投资数额，处理能力，运转成本都有差异，我们依据建设成本最小可以大致的确定处理厂的数量，进而根据总成本最小建立多变量的非线性规划问题进行求解。

对于第二问的中转站问题，我们首先把小区简化，依据位置和小区的人口数量进行决定中转站的数量和位置，剩下来的过程同第一问类似。

33基本假设基本假设和符号说明和符号说明3.13.1基本假设基本假设假设1：

题目所给数据及所指示的网站数据真实可靠。

假设2：

南山区的道路在短时间内不会有大的变动。

假设3：

所有厨余设备都能正常工作，不存在损坏或是休眠问题。

假设4：

厨余设备在各个地方建设时不存在建设成本的差异，即忽略选址点对建设成本的影响。

假设5：

点之间的道路畅通无阻，不考虑车辆空载的问题。

假设6：

由于题目中给出的小区数目过于庞大，为了便于操作，在不影响对问题分析正确性的前提下将小区数目做出简化，并将简化掉的小区人数合并在一起。

假设7：

南山区短时间内与资料做给人口不会有太大出入，没有大量人口的流入或者流出。

假设8：

各个车辆在运输的时候自动会选择最优路线，即耗油量最小的路线，并在此基础上把实际地图简化为平面上的散点图，标出各个中转站的相对位置，以两点间的直线来代替车辆实际的行走路线。

假设9：

各小区的垃圾都是当天清理完毕，没有积累。

假设10：

由题目所给的804吨垃圾/南山区的总人数=人均日产垃圾量是稳定的。

假设11：

不考虑其他附加因素对厨余设备建设地点的干扰。

3.23.2符号说明符号说明

（1）ia,it,m均为含有38项的列向量

（2）ia各站需要运输的次数（3）it取0或1，取1表示该垃圾站的垃圾运往第一垃圾处理中心，取0表示运往第二垃圾处理中心（4）im各个垃圾站的橱余垃圾的数量（5）x1是第一个处理中心的经度（6）y1是第一个处理中心的纬度（7）x2是第二个处理中心的经度（8）y2是第二个处理中心的纬度（9）jid（j=1，2，i=1，2，,38）表示第i个中转站的垃圾运往第j个处理中心的距离（10）ix，iy（i=1，2，,38）表示第i个中转站的经度和纬度44模型的建立模型的建立与求解与求解4.14.1现有条件下厨余设备的数量求解和位置分布现有条件下厨余设备的数量求解和位置分布4.1.14.1.1大小型厨余设备数量的求解大小型厨余设备数量的求解由题目中中所给的数据知道问题一中共有321.6吨厨余垃圾需要使用大小型厨余设备进行处理，如果建两个个大型厨余设备，需投资9000万元，可以处理完所有厨余垃圾，如果建一个大型处理设备，其余都建成小型的需要406608个小型设备，投资为1136817024万，则设备投资总额就为1586821524万元。

由经纬度和物理知识我们可以算出，纬度相差1度时实际距离相差111.7km，经度相差1度时实际距离相差117*cos（X），X表示所在的纬度，所以欧氏距离下得经纬度相差1度实际距离相差152.06km。

由于南山区两点间的最远距离为15km，则运输成本最大为15/100*7.5*30*321.6=10854元，远小于上述两种方案的最低成本差6868万元，由此可得，建设两个大型厨余设备是比较合理的。

4.1.24.1.2大型厨余设备位置的设计大型厨余设备位置的设计由上述的求解我们知道建设两个大型处理设备是最合理的，那么这两个大型处理设备如何安排才能达到从各个中转站到这两个处理中心的路线最优呢？

我们考虑如下：

首先建立目标函数min，表示总的最小成本，因为考虑南山区的实际道路情况很困难，所以假设司机都是知道选择最省油的路线的，我们用各转运站与厨余垃圾处理中心的直线距离来表示,所以建立目标函数如下：

目标函数：

38121min

（1）iiiiiiiadtadt=+约束条件：

t01i为（）变量；381m200iiit=；381

（1）200;iiimt=221

（1）

（1）iiidxxyy=+222

（2）

（2）iiidxxyy=+按照上述给定的目标函数和约束条件,我们利用lingo软件进行编程求解，结果如下:

VariableValueReducedCostX1113.96520.000000Y122.564720.000000X2113.9137-2.116664Y222.51550-3.807164T

（1）0.0000000.000000T

（2）1.0000000.000000T（3）1.0000000.000000T（4）1.0000000.000000T（5）1.0000000.000000T（6）0.0000000.000000T（7）0.0000000.000000T（8）0.0000000.000000T（9）0.0000000.000000T（10）0.0000000.000000T（11）0.0000000.000000T（12）0.0000000.000000T（13）1.0000000.000000T（14）0.0000000.000000T（15）0.0000000.000000T（16）1.0000000.000000T（17）0.0000000.000000T（18）0.0000000.000000T（19）1.0000000.000000T（20）0.0000000.000000T（21）0.0000000.000000T（22）0.0000000.000000T（23）0.0000000.000000T（24）0.0000000.000000T（25）1.0000000.000000T（26）1.0000000.000000T（27）1.0000000.000000T（28）1.0000000.000000T（29）1.0000000.000000T（30）0.0000000.000000T（31）0.0000000.000000T（32）1.0000000.000000T（33）1.0000000.000000T（34）1.0000000.000000T（35）1.0000000.000000T（36）0.0000000.000000T（37）1.0000000.000000T（38）1.0000000.000000T（i）（i=1、2、38）表示第i个小区，其后所跟value值为1或0，1表示此中转站的垃圾运到第一个垃圾处理中心，0表示此中转站的垃圾运到第二个垃圾处理中心，这样就将两个垃圾处理中心与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福光站113.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.621214北头站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绘图功能，以经度为横坐标，纬度为纵坐标，画出具体的位置，如图1所示：

注：

1和2表示两个大型垃圾处理中心图1厨余设备安放位置示意图4.24.2转运站允许重新设计时厨余设备的数量求解及分布转运站允许重新设计时厨余设备的数量求解及分布第二问与第一问的不同点是垃圾转运站允许重新设计，我们借鉴第一问的部分结论进行推广，在转运站的位置可变的条件下，根据各个社区的具体位置来确定最佳的中转站分布情况。

由于中转站的数量和位置的不确定性，使得这个问题的复杂程度上升。

为此我们在处理整个问题之前，先把中转站的数量和位置决定下来。

接下来决策厨余垃圾处理厂的建立与选址问题，同第一问一样，用类似的方法进行讨论解决。

4.2.14.2.1各社区的初步分布各社区的初步分布由于题目中给的小区数量过万，为了简化起见，我们按照小区名字的相关度编程将上万个小区归并合整简化为如下的若干个小区，如表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绘图功能，我们以这些小区的经度为横坐标，纬度为纵坐标，按照各个小区经纬度在图上标示出其相对位置，从图中可以看出各小区的分布稠密情况，如图2所示：

图2初次归并后各小区位置示意图4.2.24.2.2各社区的二次归并各社区的二次归并南山区每日垃圾清运总量1280吨，总人口为1320411人，分布在124个地区内（按照第一问的划分），所以可知人均日产厨余垃圾量为1280/1320411=0.9694（kg）,按照人口密度，将所有小区重新分片，所得结果与初次分布的结果进行逐一对应，如下表3所示：

表3二次归并后新旧小区对照示意图及相关数据分布表新编号人数经度度度垃圾总量（吨）厨余垃圾（吨）原始编号118150113.931622.5379517.594617.03784436,38,45,46,55,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.925322.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,88102247113.922.5321.79118.7164517,52,79,118,12097997824267113624113.943622.54183.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.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.969045124251266113.914922.76421.22726040.490904125265636113.895722.47325.46353842.18541526271106113.941422.554731.07215640.04288637283544113.949222.52683.43555360.13742232由上表可知，在处理第二问时将所有小区划分为28个片区，在这28个区域内我们假设放置的有足够接收这一片区的垃圾中转站，这样，问题就转化为在这28个片区的基础上考虑需要多少个垃圾中转站和垃圾厨余设备，以及怎样放置这些厨余设备。

厨余垃圾一共有1280*0.4=512吨，按照第一问关于厨余设备数量的求解过程（方法类似，这里不再赘述），我们知道还是全部修建成大型厨余设备是最省成本的，这样就需要三个大型设备即可以满足需求。

4.2.34.2.3大型厨余设备位置的求解大型厨余设备位置的求解由4.2.2得知需建三个大型厨余设备，设三个大型厨余设备的编号分别为1，2，3，根据题设条件建立目标函数目标函数：

381231min（12+3）iiiiiiiiiiadtadtadt=+约束条件：

1,2,3ttt，均为（01）变量；381m1200iiit=；381m3200iiit=；381m3200iiit=；for（n（i）:

t1（i）+t2（i）+t3（i）1）;利用lingo进行求解，得出结果如下：

结果：

Localoptimalsolutionfound.Objectivevalue:

2.623194Objectivebound:

2.623194Infeasibilities:

0.000000Extendedsolversteps:

0Totalsolveriterations:

22VariableValueReducedCostX1113.9619-0.7683643Y122.53396-1.206360X2113.9619-1.054813Y222.53396-0.6262339X3113.96190.000000Y322.533960.000000T1（18）1.0000000.000000T1（23）1.0000000.000000T1（24）1.0000000.000000T2（20）1.0000000.000000T2（22）1.0000000.000000T2（25）1.0000000.000000T3

（1）1.0000000.000000T3

（2）1.0000000.000000T3（3）1.0000000.000000T3（4）1.0000000.000000T3（5）1.0000000.000000T3（6）1.0000000.000000T3（7）1.0000000.000000T3（8）1.0000000.000000T3（9）1.0000000.000000T3（10）1.0000000.000000T3（11）1.0000000.000000T3（12）1.0000000.000000T3（13）1.0000000.000000T3（14）1.0000000.000000T3（15）1.0000000.000000T3（16）1.0000000.000000T3（17）1.0000000.000000T3（19）1.0000000.000000T3（21）1.0000000.000000其中T1（*）、T（*）、T（*）分别表示三个垃圾厨余设备，括号里面的数字表示上述28个片区运到该厨余设备所在处进行处理，也可将结果写成如下形式，便于观看，同时将某个地区全部建成小型厨余设备所需要的数量也一并写出，便于我们对比，从而也可看出，建一个大型设备比建若干个小型设备要省去很多成本表4各个片区与运往的厨余中心对照表及大小型厨余设备的对照表新编号橱余垃圾清运次数经度纬度大型厨余设备的编号当前编号小型个数123270.0428860280.1374220250.4909040141.2136891113.987922.5744114151.2722411113.969122.60175125111.4052421113.943622.5418135241.9690451114.027322.5528147232.1788231114.010622.5963158262.1854151113.895722.473216873.5030241113.929222.51351712204.5534661113.887622.5086181585.4433751113.966422.54333191826.0750361113.947122.54521102117.0378441113.931622.5379511124227.2910511113.987922.589211225217.318971113.900622.54