城市生活垃圾处理社会总成本分析数学建模论文

《城市生活垃圾处理社会总成本分析数学建模论文》由会员分享，可在线阅读，更多相关《城市生活垃圾处理社会总成本分析数学建模论文（23页珍藏版）》请在装配图网上搜索。

1、-市城市生活垃圾处理社会总本钱分析摘要随着各大城市的不断开展，城市居民生活垃圾产生量也随之增长，计算处理生活垃圾的社会总本钱对指导垃圾处理具有重要意义。本文针对市垃圾处理现状，建立并改良市城市垃圾处理社会总本钱分析模型，给出市垃圾处理的建议。针对问题一，要求建立市城市垃圾处理社会总本钱分析模型，结合市城市生活垃圾燃烧社会本钱评估报告20170321和市生活垃圾处理工作志愿者调研笔记的调研方法和数据，建立生命周期评价LCA模型。通过查阅市统计年鉴与市基准地价，得到市的地价与人均可支配收入等条件，确立市城市垃圾处理社会总本钱。模型包括土地本钱、收运本钱、分类处理本钱、处理设施本钱、安康本钱和电价补

2、贴，将各类本钱加和得到市城市垃圾处理社会总本钱。针对问题二，要求完善问题一的模型并计算当前及未来十年的社会总本钱及各本钱的变化趋势，结合上述文件中的调研方法和数据对问题一中的模型进展完善，将四种模式各三个时期下的各个本钱进展详细的求解，得到具体的市生活垃圾处理诸模式的直接本钱估算方法。利用MATLAB进展求解，分析结果发现：未来十年的社会总本钱呈现升高趋势；各个社会本钱的比例有所变化，安康本钱占得比例最高；近期各模式社会总本钱根本一样，中远期模式三的社会总本钱最低，模式二的社会总本钱最高。各模式的当期社会总本钱以及未来十年的总本钱数量见表3、4、5和6，诸模式下各分项本钱比例的变化趋势见图1、

3、2、3和4。本文建立市城市垃圾处理社会总本钱模型并较好的分析了各个处理模式的社会总本钱，并且对模型的优缺点进展了合理的评价，对市城市垃圾的处理方式具有参考价值。关键词：垃圾处理；本钱分析；生命周期评价；安康本钱；MATLAB. z-1. 问题重述2017年3月18日，国务院向全国发布了生活垃圾分类制度实施方案，这标志着中国垃圾分类制度建立开场了一个全新阶段，垃圾分类已成为推进社会经济绿色开展、提升城市管理和效劳水平、优化人居环境的重要举措。为了保证这一目标能够顺利实现，必须对城市生活垃圾分类处理包括但不限于分类投放、分类收集、分类运输、分类处理的全过程总本钱，进展科学计算及动态监测。请收集相关

4、资料、参考附件材料，完成下述工作：1、建立市城市垃圾处理社会总本钱分析模型；2、基于问题一的模型，完善附件1中提及的市生活垃圾处理诸模式的直接本钱估算方法，并估算各模式的当期社会总本钱以及未来十年的总本钱数量、及诸模式下各分项本钱比例的变化趋势。2. 模型假设1) 由于难以获取填埋场的部数据，假设部本钱为0；2) 假设燃烧灰渣的体积为1015%的中间值，即为12.5%；3) 假设燃烧补贴逐渐减少；4) 假设近期所有模式按照现状模式运行。3. 通用符号说明序号符号符号说明1第i种方法处理的垃圾质量2第i种方法处理垃圾的价格3第i种本钱4未来第n年市人口总量5未来第n年市人均可支配收入4.市城市垃

5、圾处理社会总本钱分析模型的建立4.1问题分析本文要求建立市城市垃圾处理社会总本钱分析模型，拟建立生命周期评价LCA模型。首先对附件1进展了详细的分析，结合附件2得到社会总本钱包括：固定本钱、可变本钱和安康本钱。其中固定本钱分为土地本钱和建立本钱，可变本钱分为垃圾处理费、电价补贴、渗沥液补贴、底灰处理补贴和飞灰处理补贴，查阅有关市的文件等，得到关于市的各个价格，将各类本钱相加，得到最终的市城市垃圾处理社会总本钱分析模型。4.2模型建立通过对问题的分析，建立生命周期评价模型，建模如下：1土地本钱：查阅基准地价得到2013年各地区地价，估算各垃圾填埋场的平均地价为，估算四年郊区地价涨幅为。设2017

6、年为未来第0年，每年地价涨幅为8.75%，则未来第n年的土地本钱为4-1其中A为每年填埋垃圾需要占的地方。2收运本钱根据附件1中城管部门的统计数据，市生活垃圾目前的根本收运本钱大约为60元/吨10公里以，10公里以外的增量收运本钱为1元/吨公里，20公里以外的增量收运本钱为1.5元/吨公里。得到公式4-2其中为一年的生活垃圾数目；为垃圾填埋场距离市区的平均运距。3分类处理本钱根据附件1得到各种处理方法的本钱：燃烧为：近期，中期，远期；填埋为；发酵为；单独处理的厨余垃圾为；源头分类收集为。总处理本钱为：4-3其中为第i个处理方法的均价；为各种方法处理的垃圾质量。4处理设施本钱目前垃圾燃烧厂多采用

7、BOT方式建立，政府在此过程中根本不需要进展投资。因此，政府的投资本钱可视为0，相应的全过程本钱可全部计入处理本钱。4-45电价补贴假设用于燃烧发电的生活垃圾数目为吨，上网电量为限额量280 kWh。根据附件2得到三个发电厂的发电量自用率，取17.5%的自用率计算。电价0.68元/kWh，燃煤发电上网电价为0.4421元/kWh，即每度电补贴0.2379元计算，补贴额度为万元4-56安康本钱查询市统计年鉴得市面积为1996.85平方公里，经过计算得到表1。表1 市统计结果年份人口万人人口增长率%人均可支配收入/元人均可支配收入增长率%人口密度20141077.891.4409485398201

8、51137.875.6446339.0569820161190.844.6486959.15964根据表1得：取各年人口增长率的平均值3.97%为市人口增长率即人口密度增长率，之后的每年均按照此增长率计算；取人均可支配收入增长率的平均值9.05%为以后的人均可支配收入增长率。a)市未来第年的人口为万人4-6假设人口平均分布，1125个坐标点的致癌风险代表其所在网格的风险，则每个网格的人数为万人。参考附件2市的计算结果，得到市致癌人数之和，为人/年。b)未来人均可支配收入为元4-7将收入参数代入工资风险模型中可得市个体生命价值的估计为万元。2015年癌症次均住院费用为41314.48元治疗费+非

9、治疗费，按6%的折现率折算到未来第年为万元。由12得如以人均住院1次计算，未来第年人安康本钱为万元4-8虽然安康损失核算在污染源、排放、扩散、暴露等估计都存在不确定性，但依然是本钱较低、解释性强的衡量方法。因为，污染物监测方法的本钱更高，且难以监测到低于一定限值的污染物，而流行病统计方法无法解释污染物排放与发病率的关系。本研究已通过保守估计选择参数防止造成对安康损失的过高估计。因而垃圾填埋场的年本钱为4-94-105. 不同模式下垃圾处理总本钱模型的建立、求解5.1问题分析本问题要求完善问题一的模型并求解诸模式下当期社会总本钱以及未来十年的总本钱数量，根据问题一的模型，对附件1给出的计算方法进

10、展具体的求解，将实际情况带入模型中，对查不到的条件进展忽略，最后得到完善后的模型，之后求解。5.2模型准备根据题意预测未来十年的市城市生活垃圾总量为5-1其中为2014年市城市生活垃圾总量；为生活垃圾年增长率，近中远期分别为6%、4%、3%。编写MATLAB程序求解模型得到当前及未来十年的垃圾总量。程序见附录1。表2 预测垃圾总量年份20172018201920202021202220232024202520262027垃圾总量/万吨644.51683.17724.17767.62798.32830.26863.47898.00933.92961.94990.80表2为预测的今年以及未来十年的

11、年垃圾总量，每一年的垃圾总量都在增加，十年间约增长了350万吨。5.2模型建立基于问题立的模型以及附件1的分析，将当期和未来十年一共11年分为近期、中期和晚期，近期为20172020年，中期为20212025年，远期为20262027年。将模型分为四种处理模式的模型。市目前在运行的燃烧厂共4座，合计设计处理能力7125吨/日，按年运行时间为300天计全年燃烧厂共可处理垃圾约215万吨。近、中、远期的生活垃圾产生量的年平均增长速度分别为6%、4%和3%，可预测近、中、远期的预测年产生规模分别为783万吨、953万吨、1105万吨。1.现状模式混合收集+局部卫生填埋+局部燃烧+局部简易堆填在现状模

12、式下，超出现有垃圾燃烧处理能力之外的垃圾都采用填埋方式处理。中期在东部新增一座卫生填埋场，平均运距为40公里，超出现有垃圾燃烧处理能力的垃圾全部运往东部填埋场进展处理。表3 现状模式近期社会总本钱分析本钱公式解析土地本钱填埋场的垃圾堆填密度一般为0.8吨/立方米，平均填埋高度按20米计，垃圾燃烧的用地需求约为3050平方米/吨日处理规模，可计算每年用地面积。收运本钱由于下坪填埋场距市区约为10公里，老虎坑填埋场距市区约为20公里，下坪填埋场与老虎坑填埋场承当的垃圾处理负荷为1:1，可得收运本钱。分类处理本钱现状全年燃烧厂共可处理垃圾约215万吨，其余的垃圾直接掩埋。处理设施本钱新建填埋场需投资

13、约25.0亿元。按使用年限16年计，平均每年约需1.56亿元。电价补贴目前燃烧的垃圾总量为215万吨安康本钱由垃圾燃烧后产生的二噁英导致表4 现状模式中期社会总本钱分析本钱公式解析土地本钱同近期收运本钱由于要在距市区40公里的地方建填埋场，超出现有垃圾燃烧处理能力的垃圾全部运往这里，所以本钱增加。分类处理本钱现状全年燃烧厂共可处理垃圾约215万吨，其余的垃圾直接掩埋，但燃烧处理费上涨。处理设施本钱新建填埋场需投资约25.0亿元。按使用年限16年计，平均每年约需1.56亿元。电价补贴目前燃烧的垃圾总量为215万吨，中期补贴减半。安康本钱由垃圾燃烧后产生的二噁英导致。表5 现状模式远期社会总本钱分

14、析本钱公式解析土地本钱同中期收运本钱同中期分类处理本钱现状全年燃烧厂共可处理垃圾约215万吨，其余的垃圾直接掩埋，但燃烧处理费继续上涨。处理设施本钱新建填埋场需投资约25.0亿元。按使用年限16年计，平均每年约需1.56亿元。电价补贴由于电价补贴无益于垃圾分类，后期取消补贴。安康本钱由垃圾燃烧后产生的二噁英导致。2.模式一混合收集+全量燃烧+灰渣填埋+中心城区垃圾全量转运模式一以燃烧为主，现在采用填埋方式处理的垃圾全部改用燃烧方式处理。由于燃烧灰渣物化性质稳定且不臭，灰渣填埋场的选址相比卫生填埋场的选址要容易得多，且燃烧灰渣容易通过循环利用的方式重新用于城市建立活动。因而不考虑灰渣填埋的面积。

15、表6 模式一近期社会总本钱分析本钱公式解析土地本钱填埋场的垃圾堆填密度一般为0.8吨/立方米，平均填埋高度按20米计，垃圾燃烧的用地需求约为3050平方米/吨日处理规模，可计算每年用地面积。收运本钱由于下坪填埋场距市区约为10公里，老虎坑填埋场距市区约为20公里，下坪填埋场与老虎坑填埋场承当的垃圾处理负荷为1:1，可得收运本钱。分类处理本钱现状全年燃烧厂共可处理垃圾约215万吨，其余的垃圾直接掩埋。处理设施本钱考虑到目前垃圾燃烧厂多采用BOT方式建立，政府在此过程中根本不需要进展投资。因此，政府的投资本钱可视为0。电价补贴目前燃烧的垃圾总量为215万吨。安康本钱由垃圾燃烧后产生的二噁英导致。表

16、7 模式一中期社会总本钱分析本钱公式解析土地本钱垃圾燃烧的用地需求约为3050平方米/吨日处理规模，可计算每年用地面积。收运本钱由于要在距市区40公里的地方建填埋场，超出现有垃圾燃烧处理能力的垃圾全部运往这里，所以本钱增加分类处理本钱将垃圾全部进展燃烧，将燃烧后产生的灰渣掩埋。处理设施本钱考虑到目前垃圾燃烧厂多采用BOT方式建立，政府在此过程中根本不需要进展投资。因此，政府的投资本钱可视为0。电价补贴由于电价补贴无益于垃圾分类，中期燃烧厂的电价补贴减半。安康本钱由垃圾燃烧后产生的二噁英导致。表8 模式一远期社会总本钱分析本钱公式解析土地本钱同中期收运本钱新建的燃烧厂距市区约40公里，将之前燃烧

17、厂承受不了的垃圾都运往新燃烧厂。分类处理本钱将垃圾全部进展燃烧，将燃烧后产生的灰渣掩埋。处理设施本钱考虑到目前垃圾燃烧厂多采用BOT方式建立，政府在此过程中根本不需要进展投资。因此，政府的投资本钱可视为0。电价补贴由于电价补贴无益于垃圾分类，远期燃烧厂的电价补贴取消。安康本钱由垃圾燃烧后产生的二噁英导致。3.模式二源头分类收集+湿垃圾生物处理+干垃圾燃烧+中心城区干垃圾转运模式二是在垃圾产生源头即实施分类收集，垃圾在源头分为干、湿两大类。因此，收集到的垃圾相应为干、湿两大类，考虑到居民参与分类收集的积极性和垃圾投放的准确度，湿垃圾的比例按40%考虑，干垃圾的比例按60%考虑。干的局部采用燃烧方

18、式处理，湿的局部暂考虑采用厌氧发酵技术处理。表9 模式二近期社会总本钱分析本钱公式解析土地本钱填埋场的垃圾堆填密度一般为0.8吨/立方米，平均填埋高度按20米计，垃圾燃烧的用地需求约为3050平方米/吨日处理规模，可计算每年用地面积。收运本钱由于下坪填埋场距市区约为10公里，老虎坑填埋场距市区约为20公里，下坪填埋场与老虎坑填埋场承当的垃圾处理负荷为1:1，可得收运本钱。分类处理本钱现状全年燃烧厂共可处理垃圾约215万吨，其余的垃圾直接掩埋。处理设施本钱考虑到目前垃圾燃烧厂多采用BOT方式建立，政府在此过程中根本不需要进展投资。因此，政府的投资本钱可视为0。电价补贴目前燃烧的垃圾总量为215万

19、吨。安康本钱由垃圾燃烧后产生的二噁英导致。表10 模式二中期社会总本钱分析本钱公式解析土地本钱垃圾燃烧的用地需求约为3050平方米/吨日处理规模，厌氧发酵处理考虑集约用地的原则按6080平方米/吨日处理规模，可计算每年用地面积。收运本钱由于厨余垃圾处理设施运距为10公里、垃圾燃烧厂运距40公里。分类处理本钱源头分类收集的补贴标准约为423元/吨，60%的垃圾燃烧，且燃烧处理费上涨，40%的垃圾发酵处理。处理设施本钱考虑到目前垃圾燃烧厂多采用BOT方式建立，政府在此过程中根本不需要进展投资。因此，政府的投资本钱可视为0。电价补贴由于电价补贴无益于垃圾分类，中期燃烧厂的电价补贴减半。安康本钱由垃圾

20、燃烧后产生的二噁英导致。表11 模式二远期社会总本钱分析本钱公式解析土地本钱垃圾燃烧的用地需求约为3050平方米/吨日处理规模，厌氧发酵处理考虑集约用地的原则按6080平方米/吨日处理规模，可计算每年用地面积。收运本钱由于厨余垃圾处理设施运距为10公里、垃圾燃烧厂运距40公里。分类处理本钱源头分类收集的补贴标准约为423元/吨，60%的垃圾燃烧，40%的垃圾发酵处理。处理设施本钱考虑到目前垃圾燃烧厂多采用BOT方式建立，政府在此过程中根本不需要进展投资。因此，政府的投资本钱可视为0。电价补贴由于电价补贴无益于垃圾分类，远期燃烧厂的电价补贴取消。安康本钱由垃圾燃烧后产生的二噁英导致。4.模式三混

21、合收集+末端分类+湿垃圾生物处理+干垃圾燃烧+中心城区干垃圾转运模式三是在末端才通过人工分拣或垃圾压榨实现干、湿垃圾别离。因此，在末端收集到的垃圾相应为干、湿两大类，可按各50%重量考虑。干湿垃圾处理同模式二。表12 模式三近期社会总本钱分析本钱公式解析土地本钱填埋场的垃圾堆填密度一般为0.8吨/立方米，平均填埋高度按20米计，垃圾燃烧的用地需求约为3050平方米/吨日处理规模，可计算每年用地面积。收运本钱由于下坪填埋场距市区约为10公里，老虎坑填埋场距市区约为20公里，下坪填埋场与老虎坑填埋场承当的垃圾处理负荷为1:1，可得收运本钱。分类处理本钱现状全年燃烧厂共可处理垃圾约215万吨，其余的

22、垃圾直接掩埋。处理设施本钱考虑到目前垃圾燃烧厂多采用BOT方式建立，政府在此过程中根本不需要进展投资。因此，政府的投资本钱可视为0。电价补贴目前燃烧的垃圾总量为215万吨。安康本钱由垃圾燃烧后产生的二噁英导致。表13 模式三中期社会总本钱分析本钱公式解析土地本钱垃圾燃烧的用地需求约为3050平方米/吨日处理规模，厌氧发酵处理考虑集约用地的原则按6080平方米/吨日处理规模，可计算每年用地面积。收运本钱由于厨余垃圾处理设施运距为10公里、垃圾燃烧厂运距40公里。分类处理本钱源头分类收集的补贴标准约为423元/吨，60%的垃圾燃烧，且燃烧处理费上涨，40%的垃圾发酵处理。处理设施本钱考虑到目前垃圾

23、燃烧厂多采用BOT方式建立，政府在此过程中根本不需要进展投资。因此，政府的投资本钱可视为0。电价补贴由于电价补贴无益于垃圾分类，中期燃烧厂的电价补贴减半。安康本钱由垃圾燃烧后产生的二噁英导致。表14 模式三远期社会总本钱分析本钱公式解析土地本钱垃圾燃烧的用地需求约为3050平方米/吨日处理规模，厌氧发酵处理考虑集约用地的原则按6080平方米/吨日处理规模，可计算每年用地面积。收运本钱由于厨余垃圾处理设施运距为10公里、垃圾燃烧厂运距40公里。分类处理本钱60%的垃圾燃烧，40%的垃圾发酵处理。处理设施本钱考虑到目前垃圾燃烧厂多采用BOT方式建立，政府在此过程中根本不需要进展投资。因此，政府的投

24、资本钱可视为0。电价补贴由于电价补贴无益于垃圾分类，远期燃烧厂的电价补贴取消。安康本钱由垃圾燃烧后产生的二噁英导致。各项社会本钱占总社会本钱比例为5-1其中表示第个模式表示第个本钱5.3模型求解编写MATLAB程序求解得到每个模式的各个本钱以及社会总本钱。程序见附录26。表15 各模式现状及未来十年社会总本钱年现状模式一模式二模式三总本钱亿元平均本钱元总本钱亿元平均本钱元总本钱亿元平均本钱元总本钱亿元平均本钱元201738.21706.1936.65677.3636.65677.3636.65677.35201839.7692.0838.14664.8838.14664.8838.13664.

25、88201941.25678.4239.69652.7639.69652.7639.68652.75202042.87665.1941.31640.9841.31640.9841.31640.98202146.48680.350.09733.1569.351015.0538.12558.08202248.4668.2952.44724.273.121009.6840.09553.61202350.4656.6454.92715.4477.11004.4242.16549.23202452.84661.8957.68722.5580.531008.6944.13552.78202555.4166

26、7.3460.59729.7884.111013.0346.19556.4202658.38676.1167775.9588.51024.9548.9566.34202761.23681.8270.38783.7792.441029.4551.19570.08分析表15得到：1每一年的社会总本钱都在增加，其中模式二的社会总本钱最高,主要由于模式二需要对垃圾进展源头分类收集，而源头分类收集本钱较高。2模式三的社会总本钱最低，由于模式三进展了干湿别离，对湿垃圾进展厌氧发酵技术处理，大大减小了二噁英的产生量，安康本钱也大大降低。运用MATLAB程序画出各个模式下每类本钱的变化趋势图。图1 现状模式下

27、未来各年各类本钱变化趋势图2模式一未来各年各类本钱变化趋势图3模式二未来各年各类本钱变化趋势图4模式三未来各年各类本钱变化趋势根据图1、2、3和4分析可得：1现状模式下安康本钱所占比例大于60%，其余各项本钱都较低，因而影响现状模式总本钱的是安康本钱；未来十年每一项本钱的所占比例变化很小，因为现状模式在未来十年间处理模式不会有所变化。2模式一安康本钱所占比例大于50%，因而是主要影响因素，处理本钱在未来十年间将不断增高，因为处理本钱与垃圾总量密切相关。3模式二在近期安康本钱所占比例较大，中远期设施本钱所占比例最大；除补贴本钱之外每个本钱所占比例相差不大。4模式三中远期社会本钱下降较多，其他本钱

28、除补贴本钱外所占比例相差不大。5安康本钱在各个模式下占比例都是最大，安康本钱是影响城市垃圾处理社会总本钱的主要因素。6.模型评价此题建立了生命周期评价模型，求解了各个模式市垃圾处理社会总本钱以及各本钱的比例，该模型有优点也有缺点。优点：对各种本钱考虑到位，建立的模型较好的反映了市的垃圾处理情况；数据查找较为准确，并考虑每年的增长率，所求的的数据较为可信。缺点：有些未查找到的数据例如税收补贴等，对其进展忽略，使得结果具有偏差。参考文献1秀君，胡彬，我国生命周期评价LCA文献综述及国外最新进展A，科学技术进步及对策，2013年，306：155-160。2市城市生活管理局，市生活垃圾燃烧发电厂一览表

29、2016年6月27日，.szum.gov./zfwg/tjsj/z*tj*/201606/t20160629_3732714.htm，2016年6月29日。2017年8月22日。3王波，基于生命周期评价的市建筑垃圾处理模式研究D，省市，华中科技大学，2012年。4田爱军,冰,新玲,敏,黄夏银，生活垃圾燃烧烟气排放中二噁英对人体安康的风险评价A，污染防治技术，2008年，216：26-28。5洪大用，市城市生活垃圾燃烧社会本钱评估报告R，中国人民大学，2017年。6市基准地价2013年，.szpl.gov./*gk/tdgl/jzdj/2013dj/，2017年8月22日。7市统计年鉴2016年

30、，.sztj.gov./*gk/tjsj/tjnj/201701/t20170120_5943580.htm，2017年1月20日，2017年8月22日。8市统计年鉴2015年，.sztj.gov./*gk/tjsj/tjnj/201606/t20160607_3686823.htm，2016年6月7日，2017年8月22日。9市统计年鉴2014年，.sztj.gov./*gk/tjsj/tjnj/201503/t20150331_2836397.htm，2015年3月31日，2017年8月22日。附录附录1 未来每年垃圾总产量程序rubbish_table = zeros(1,14);rub

31、bish_table(1,1) = 541.14*104; %2014年垃圾产量541.14万吨for year = 2015:2020rubbish_table(1,year-2013) = rubbish_table (1,year-2014)*(1+0.06);endfor year = 2021:2025rubbish_table(1,year-2013) = rubbish_table (1,year-2014)*(1+0.04);endfor year = 2026:2030rubbish_table(1,year-2013) = rubbish_table (1,year-201

32、4)*(1+0.03);end附录2 现状模式程序%第一题：建立总本钱分析模型/年：按现状分析% 总本钱=直接本钱 +经济技术本钱 + 社会本钱function dataPro = Total_Cost_Analysis(year)%垃圾每年预测表：2017-2030table = 5411400,5736084,6080249.04000000,6445063.98240000,6831767.82134400,7241673.89062464,7676174.32406212,7983221.29702460,8302550.14890559,8634652.15486181,898003

33、8.24105628,9339239.77069854,9619416.96381949,9907999.47273408,10205239.4569161,10511396.6406236,10826738.5398423; %垃圾总量每年数值2017-2030rubbish_quantity = table(year-2016);%将时间分期处理：2017-2020，2021-2025，2026-2030switch year case 2017,2018,2019,2020 %近期 rubbish_num_burn=215*104; rubbish_num_landfill = rubb

34、ish_quantity-rubbish_num_burn; land_cost = (rubbish_quantity*0.0625+10.12*104)*573.75; class_cost = 0; handle_cost = rubbish_num_landfill*60+ rubbish_num_burn*100; transport_cost = 0.5*rubbish_num_landfill*60+0.5*rubbish_num_landfill*70+0.5*rubbish_num_burn*60+0.5*rubbish_num_burn*70 ; social_cost =

35、(rubbish_num_burn/rubbish_quantity)*132089.23*(1+0.0379)(year-2016)*(71.62*(1+0.0905)(year-2016)+41314.48*(1+0.06)(year-2016+1); subsidy = 54.9549*rubbish_num_burn; %前期 100，中期50，后期取消，本钱计算取负 case 2021,2022,2023,2024,2025 %中期 rubbish_num_burn =215*104; rubbish_num_landfill = rubbish_quantity-rubbish_n

36、um_burn; land_cost = (rubbish_quantity*0.0625+10.12*104)*573.75; class_cost = 0; handle_cost = rubbish_num_landfill*60+ rubbish_num_burn*140; transport_cost = rubbish_num_landfill*100+ rubbish_num_burn*65; social_cost =(rubbish_num_burn/rubbish_quantity)*132089.23*(1+0.0379)(year-2016)*(71.62*(1+0.0

37、905)(year-2016)+41314.48*(1+0.06)(year-2016+1); subsidy = 0.5*54.9549*rubbish_num_burn; %前期 100，中期50，后期取消，本钱计算取负 case 2026,2027,2028,2029,2030 %远期 rubbish_num_burn =215*104; rubbish_num_landfill = rubbish_quantity-rubbish_num_burn; land_cost = (rubbish_quantity*0.0625+10.12*104)*573.75; class_cost =

38、 0; handle_cost = rubbish_num_landfill*60+ rubbish_num_burn*180; transport_cost = rubbish_num_landfill*100+ rubbish_num_burn*65; social_cost =(rubbish_num_burn/rubbish_quantity)*132089.23*(1+0.0379)(year-2016)*(71.62*(1+0.0905)(year-2016)+41314.48*(1+0.06)(year-2016+1); subsidy = 0*54.9549*rubbish_n

39、um_burn; %前期 100，中期50，后期取消，本钱计算取负 otherwise msgbo*(亲，请重新输入年份：); end%设施投资: equipment_cost = 1.56*108; %输出，分析:dataPro为数据集合direct_cost = land_cost + class_cost + transport_cost + equipment_cost + handle_cost;total_cost = direct_cost+social_cost+subsidy ;%dataPro11: 分类，收运，设施，处理，技术，社会，补贴，收益，直接，总，均dataPro

40、 = land_cost,class_cost,transport_cost,equipment_cost,handle_cost, . social_cost,subsidy,direct_cost,total_cost,total_cost/rubbish_quantity;end附录3 模式一程序%模式一：总本钱=直接本钱 +经济技术本钱 + 社会本钱function dataPro = Total_Cost_Analysis_model1(year)%垃圾每年预测表：2017-2030table = 5411400,5736084,6080249.04000000,6445063.98

41、240000,6831767.82134400,7241673.89062464,7676174.32406212,7983221.29702460,8302550.14890559,8634652.15486181,8980038.24105628,9339239.77069854,9619416.96381949,9907999.47273408,10205239.4569161,10511396.6406236,10826738.5398423; %垃圾总量每年数值2017-2030rubbish_quantity = table(year-2016);%将时间分期处理：2014-202

42、0，2021-2025，2026-2030switch year case 2016,2017,2018,2019,2020 rubbish_num_burn=215*104; %近期 rubbish_num_landfill = rubbish_quantity-rubbish_num_burn; land_cost = (rubbish_quantity*0.0625+10.12*104)*573.75; transport_cost = 0.5*rubbish_quantity*60+0.5*rubbish_quantity*70; handle_cost = rubbish_num_l

43、andfill*60+rubbish_num_burn*100; social_cost =(rubbish_num_burn/rubbish_quantity)*132089.23*(1+0.0379)(year-2016)*(71.62*(1+0.0905)(year-2016)+41314.48*(1+0.06)(year-2016+1); subsidy = 54.9549*rubbish_num_burn; %前期 100，中期50，后期取消，本钱计算取负 case 2021,2022,2023,2024,2025 land_cost = rubbish_quantity*62.87

44、; transport_cost = rubbish_quantity*65; handle_cost = rubbish_quantity*147.5; social_cost =(215*104/rubbish_quantity)*132089.23*(1+0.0379)(year-2016)*(71.62*(1+0.0905)(year-2016)+41314.48*(1+0.06)(year-2016+1); subsidy = 0.5*54.9549*rubbish_quantity; %前期 100，中期50，后期取消，本钱计算取负 case 2026,2027,2028,2029

45、,2030 land_cost = rubbish_quantity*62.87; transport_cost = (rubbish_quantity-215*104)*100+65*215*104; handle_cost = rubbish_quantity*187.5; social_cost =(215*104/rubbish_quantity)*132089.23*(1+0.0379)(year-2016)*(71.62*(1+0.0905)(year-2016)+41314.48*(1+0.06)(year-2016+1); subsidy = 0*54.9549*rubbish

46、_quantity; %前期 100，中期50，后期取消，本钱计算取负 otherwise msgbo*(亲，请重新输入年份：); end%分类费用class_cost = 0; %设施投资: equipment_cost = 0 ; %输出，分析 direct_cost = land_cost+ class_cost + transport_cost + equipment_cost + handle_cost;total_cost = direct_cost+social_cost+subsidy ;%dataPro11: 分类，收运，设施，处理，技术，社会，补贴，收益，直接，总，均dat

47、aPro = land_cost,class_cost,transport_cost,equipment_cost,handle_cost, . social_cost,subsidy,direct_cost,total_cost,total_cost/rubbish_quantity;end附录4 模式二程序%模式二：源头分类收集+湿垃圾生物处理+干垃圾燃烧+中心城区干垃圾转运function dataPro = Total_Cost_Analysis_model2(year)%垃圾每年预测表：2017-2030table = 5411400,5736084,6080249.04000000

48、,6445063.98240000,6831767.82134400,7241673.89062464,7676174.32406212,7983221.29702460,8302550.14890559,8634652.15486181,8980038.24105628,9339239.77069854,9619416.96381949,9907999.47273408,10205239.4569161,10511396.6406236,10826738.5398423; %垃圾总量每年数值2017-2030rubbish_quantity = table(year-2016);switch

49、 year case 2016,2017,2018,2019,2020 rubbish_num_burn=215*104; %近期 rubbish_num_landfill = rubbish_quantity-rubbish_num_burn; land_cost = (rubbish_quantity*0.0625+10.12*104)*573.75; class_cost = 0; transport_cost = 0.5*rubbish_quantity*60+0.5*rubbish_quantity*70; handle_cost = rubbish_num_landfill*60+

50、rubbish_num_burn*100; social_cost =(rubbish_num_burn/rubbish_quantity)*132089.23*(1+0.0379)(year-2016)*(71.62*(1+0.0905)(year-2016)+41314.48*(1+0.06)(year-2016+1); subsidy = 54.9549*rubbish_num_burn; %前期 100，中期50，后期取消，本钱计算取负 case 2021,2022,2023,2024,2025 land_cost = rubbish_quantity*84.88; class_cos

51、t = rubbish_quantity*423; transport_cost = 0.4*rubbish_quantity*60+0.6*rubbish_quantity*100; handle_cost = rubbish_quantity*148.5; social_cost =0.6*(215*104/rubbish_quantity)*(132089.23*(1+0.0379)(year-2016)*(71.62*(1+0.0905)(year-2016)+41314.48*(1+0.06)(year-2016+1); subsidy =0.5*54.9549*0.6*rubbis

52、h_quantity; %前期 100，中期50，后期取消，本钱计算取负 case 2026,2027,2028,2029,2030 land_cost = rubbish_quantity*84.88; class_cost = rubbish_quantity*423; transport_cost = 0.4*rubbish_quantity*60+0.6*rubbish_quantity*100; handle_cost = 0.4*rubbish_quantity*150+ 0.6*rubbish_quantity*180+0.6*rubbish_quantity*0.125*60;

53、 social_cost =0.6*(215*104/rubbish_quantity)*(132089.23*(1+0.0379)(year-2016)*(71.62*(1+0.0905)(year-2016)+41314.48*(1+0.06)(year-2016+1); subsidy = 0*54.9549*rubbish_quantity; %前期 100，中期50，后期取消，本钱计算取负 otherwise msgbo*(亲，请重新输入年份：); end %设施投资: equipment_cost = 0; %输出，分析 direct_cost = land_cost+ class

54、_cost + transport_cost + equipment_cost + handle_cost;total_cost = direct_cost +social_cost+subsidy ;%dataPro11: 分类，收运，设施，处理，技术，社会，补贴，收益，直接，总，均dataPro = land_cost, class_cost,transport_cost,equipment_cost,handle_cost, . social_cost,subsidy,direct_cost,total_cost,total_cost/rubbish_quantity;end附录5 模式

55、三程序%模式三：混合收集+末端分类+湿垃圾生物处理+干垃圾燃烧+中心城区干垃圾转运function dataPro = Total_Cost_Analysis_model3(year)%垃圾每年预测表：2017-2030table = 5411400,5736084,6080249.04000000,6445063.98240000,6831767.82134400,7241673.89062464,7676174.32406212,7983221.29702460,8302550.14890559,8634652.15486181,8980038.24105628,9339239.77069

56、854,9619416.96381949,9907999.47273408,10205239.4569161,10511396.6406236,10826738.5398423; %垃圾总量每年数值2017-2030rubbish_quantity = table(year-2016);switch year case 2016,2017,2018,2019,2020 rubbish_num_burn=215*104; %近期 rubbish_num_landfill = rubbish_quantity-rubbish_num_burn; land_cost = (rubbish_quant

57、ity*0.0625+10.12*104)*573.75; transport_cost = 0.5*rubbish_quantity *60+0.5*rubbish_quantity *70 ; handle_cost = rubbish_num_landfill*60+rubbish_num_burn*100; social_cost =(rubbish_num_burn/rubbish_quantity)*132089.23*(1+0.0379)(year-2016)*(71.62*(1+0.0905)(year-2016)+41314.48*(1+0.06)(year-2016+1); su