煤炭矿区节能减排多目标优化决策研究

【摘要】煤炭资源相对丰富的赋存条件决定了其在我国工业中的主导地位,然而以煤炭为主的能源消费结构导致能源利用效率低、污染物排放多,给生态环境的保护造成了极大的威胁。如今化石能源日渐枯竭,新能源的成本高且难度大。随着国际碳税的逐渐盛行,不采取节能减排措施就会被新的贸易壁垒阻碍其经济的发展。节能减排作为实现经济发展和保护环境双赢的有效途径,不仅是我国自身可持续发展的内在要求,也是为全球减缓气候变化做出的重要贡献。煤炭行业作为9大重点耗能产业之一,在节能减排工作中承担着重要的角色。“十一五”规划节能减排战略的实施,煤炭行业在提高能源效率和减少污染排放两方面都有很大的进展。为了达到“十二五”规划节能减排的目标,必须继续依靠先进科学技术的力量,改变以往高投入、高能耗、高排放的生产模式,提高能源利用效率,加强排放物的回收利用,发展循环经济。煤炭矿区资源配置系统错综复杂,在制定节能减排策略时需要考虑一系列的因素,例如经济、社会、环境、技术等。从可持续发展的角度,本文同时考虑经济效益、能源效益和环境效益三个目标,试图通过多目标优化决策的方法,对能耗和排放的关键工序进行节能设备改造,同时投资一定的项目来治理或综合利用排放物,优化配置矿区的煤炭生产计划。主要研究内容如下：(1)矿区节能减排潜力。在制定节能减排策略之前,需要首先掌握矿区目前的能源消耗结构、能源利用效率、污染物的回收和排放情况。针对“十一五”规划末的能源消耗和污染排放数据,确定投入产出指标体系,基于投入的CCR-DEA模型计算矿区的能源利用效率,评估节能减排的潜力。(2)静态多目标优化模型。选择“十二五”整个规划期的煤炭产量、关键工序设备节能改造投资、治理或综合利用项目投资作为优化的决策变量,确定经济效益、能源效益、环境效益的表达方式,描述矿区煤炭资源储量、投资资金和工序关系等约束条件,建立煤炭矿区节能减排静态多目标优化模型。(3)静态多目标模型求解算法。与单目标优化问题不同,多目标问题的优化不仅需要区分可行解和不可行解,而且还需要在多个目标之间辨别同类解的优劣。这些问题的解决依赖于约束条件的正确处理。由于多目标问题的最优解是由多个非支配解组成的Pareto解集,故在解空间内尽可能多的搜索到非支配解是得到高质量解的保证。结合NSGA-Ⅱ和PSO,充分发挥NSGA-Ⅱ的全局搜索能力和PSO局部寻优的特性,并验证PSO-NSGA-Ⅱ算法的收敛性、覆盖性、均匀性。(4)动态多目标优化模型。根据投资决策的时序性特点,以“十二五”规划期内的每一年作为优化对象,考虑节能减排的阶段性效果,将上一年的煤炭生产计划和投资决策作为下一年制定节能减排策略的依据,建立煤炭矿区节能减排动态多目标优化模型。(5)动态多目标模型求解算法。动态多目标模型的优化目标和约束条件随时间推移呈现动态变化的特点,决策环境不断发生改变,这就要求算法除了在固定的进化环境内尽可能多的搜索到Pareto解外,还要能探测到进化环境的任何微小的变化,并对环境变化做出正确的响应,确定新环境的进化参数。基于这些考虑,设计混合进化算法DNSGA-Ⅱ-PSO求解动态多目标模型,对算法探测环境变化、保持种群多样性、预测环境变化三个性能进行分析。(6)满意解的筛选。快速有效地决策是建立在少数有限的候选方案的基础之上的,所建立的多目标优化模型最终求得的是一组非劣的Pareto解集,候选方案的选取还需要对Pareto解集进行进一步筛选,这需要借助于多属性决策的方法。故运用混合聚类方法SC-MTD-GAM,在保证Pareto前端分布特性的基础上,尽量减少Pareto解集的大小,考虑决策者的偏好,选择有代表性的Pareto解作为候选方案。基于多目标优化模型和满意解筛选方法,以典型煤炭矿区超化矿的煤炭生产为应用对象,探讨矿区为了达到“十二五”规划节能减排目标如何安排煤炭生产计划,对哪些生产设备进行节能改造,投资多少项目对污染物进行治理或综合利用。得到如下结论：(1)通过对煤炭矿区能源消费结构和污染排放情况的分析,矿区主要消耗原煤、汽油、柴油、电力四种能源,对环境主要排放二氧化硫、矿井水、煤矸石三种污染物。运用CCR-DEA模型对矿区的能源效率和减排潜力进行评估,表明矿区当前处于最优的生产前沿,能源利用效率已达到最大化,进一步实施节能减排需要利用先进技术来提高生产效率。(2)建立的煤炭矿区节能减排投资多目标决策模型,满足了煤炭矿区节能减排投资的决策需要。该模型以原煤产量最大、能耗和污染排放最少为目标,考虑资源、工序、资金和环保等多个约束,较好地描述了中国现阶段典型煤矿节能减排投资的决策需要。(3)提出的PSO-NSGA-Ⅱ混合多目标求解算法,比NSGA-Ⅱ算法具有更好的收敛性、覆盖性、和均匀性。针对本文建立的多日标优化模型的决策变量类型前有0-1和实数的特点,提出一种混合PSO实数编码和NSGA-Ⅱ二进制编码的求解算法。通过质心距离、覆盖性、间隔距离三个性能指标与NSGA-Ⅱ的计算结果进行对比,说明PSO-NSGA-Ⅱ发挥了PSO和NSGA-Ⅱ两种算法的组合优势,具有更好的收敛性、覆盖性、均匀性。(4)提出的DNSGA-Ⅱ-PSO混合进化算法,比DNSGA-Ⅱ-A算法更能够探测到进化环境的任何微小的变化,保持种群多样性避免算法早熟而陷入局部最优,且能通过预测以响应环境的变化。根据动态决策的特点,本文将每个投资计划年作为决策阶段,建立含有实变量和0-1变量的动态多目标模型,提出混合算法DNSGA-Ⅱ-PSO,充分利用DNSGA-Ⅱ在全局的角度引导算法的搜索方向,PSO控制局部区域的快速寻优。此混合算法采用环境探测算子探测投资环境的任何变化,结合惯性预测、高斯分布、随机生成三种新个体产生方式完成对新环境的预测。通过对每个环境下的质心距离、覆盖性、间隔距离三个性能指标与DNSGA-Ⅱ的计算结果进行对比,说明DNSGA-Ⅱ-PSO对决策环境具有更好的探测性能和预测性。(5)运用混合聚类方法SC-MTD-GAM根据管理者的偏好对Pareto方案集进行筛选,最终得到经济偏好型、能源节约型、协调发展型三种节能减排方案。将矿区节能减排前后的煤炭生产情况进行比较分析：对于经济偏好型投资方案,矿区在达到“十二五”规划节能减排目标的前提下,保证最大的生产能力和充分利用投资资金,最大化煤炭产量；能源节约型投资方案偏重减少能源消耗量和污染物排放量,很大程度依赖于大量缩减煤炭产量的方式实现；而协调发展型方案通过降低煤炭产量、投资节能改造工程和减排项目三种方式达到节能减排的目标。无论是哪种方案,虽然煤炭产量有暂时降低的现象,但是随着节能设备的改造和综合治理利用项目的实施,节能减排的效果不断累积,煤炭产量最终呈增加的趋势,最终超过节能减排措施前的产量,说明矿区实施节能减排是一项非常必要和长期但有价值的工作,不仅能达成节能减排的目标,还从根本上优化了矿区煤炭生产的资源配置,促进矿区向“低能耗、高效率、零污染”的绿色生产模式转型。本文的创新点主要体现在以下三个方面：(1)建立了煤炭矿区节能减排多目标优化模型。同时以规划期内原煤产量最大、能源消耗和污染排放最少作为优化的三个目标,考虑资源、工序、资金和环保等多个约束,分别从静态和动态两个决策视角建立多目标优化模型,该模型较好地描述了中国现阶段典型煤矿节能减排投资的决策需要：在规划期内,为了达到既定的节能减排目标,矿区选择哪些节能设备和综合治理利用项目进行投资?何时进行节能减排投资?每个项目投资多少资金?(2)提出混合进化算法PSO-NSGA-Ⅱ和DNSGA-Ⅱ-PSO。针对建立的多目标优化模型的决策变量类型兼有0-1和实数的特点,本文提出一种混合PSO实数编码和NSGA-Ⅱ二进制编码的求解算法PSO-NSGA-Ⅱ;基于动态决策阶段性的特点,利用DNSGA-Ⅱ全局搜索性能引导解的进化方向,并通过PSO局部寻优的优点加速算法的收敛,提出DNSGA-Ⅱ-PSO。这两个混合进化算法有效地继承了NSGA和PSO二者的优点,能在解空间内尽可能大的范围内搜索到更多的非支配解,同时不断趋近真实的Pareto前端。(3)提出混合聚类方法SC-MTD-GAM。针对多目标算法求得的Pareto解集,考虑管理者偏好,在保证Pareto前端分布形状的基础上,选择具有代表性的解点作为矿区节能减排的候选方案,最终得到经济偏好型、能源节约型、协调发展型三种决策方案,对矿区制定节能减排投资计划提供科学的参考依据。
【作者】杨娟；
【导师】诸克军；
【作者基本信息】中国地质大学，管理科学与工程，2014，博士
【关键词】煤炭矿区；节能减排；多目标优化；动态多目标；多属性决策；

【参考文献】
[1].最新专利文摘[J].石油化工,2014,01:73+78+86+90+96+105+111+117+120.
[2]程先东.多方式客运通道结构演变与发展研究[D].北京交通大学,2011.
[3]葛昕.川崎病患儿NKG2D表达变化初探[D].遵义医学院,儿科学（专业学位）,2013,硕士.
[4]蔡亚超,徐雅,沈惬,王凯,邹江,孙大明.大功率斯特林制冷机的整机数值模拟[J].低温工程,2014,03:27-31+66.
[5]王雪松,田西兰,程玉虎,易建强.基于协同最小二乘支持向量机的Q学习[J].自动化学报,2009,02:214-219.
[6]贾晓亮.聚类分析在我国矿业管理中的应用研究[D].内蒙古科技大学,企业管理,2013,硕士.
[7]张欢.中国资产价格对通货膨胀的影响：理论与实证研究[D].首都经济贸易大学,国民经济学,2014,硕士.
[8]dwway姜海.数据仓库在SCM中的应用[N].中国计算机报,2002-12-23B20.
[9]张月战.血清游离轻链与慢性阻塞性肺病疾病严重度及分期的相关性研究[D].苏州大学,内科学,2014,硕士.
[10]丁祥.经口咽前路寰枢椎解剖型锁定钢板的生物力学研究[D].蚌埠医学院,外科学,2013,硕士.
[11]李慧玲.二穗短柄草（Brachypodium distachyon）T-DNA插入突变体库构建[D].南京农业大学,园林植物与观赏园艺,2012,硕士.
[12]莫中杰.我国房地产上市公司多元化战略对企业绩效影响的实证研究[D].扬州大学,企业管理,2011,硕士.
[13]徐坤.高中生物小班化教学的研究[D].山东师范大学,学科教学（专业学位）,2013,硕士.
[14]张卉.趋化因子SDF-1α介导的胰腺星形细胞促胰腺癌细胞化疗耐药作用及其机制研究[D].北京协和医学院,病理学与病理生理学,2014,博士.
[15]陈慧萍,王建东,叶飞跃.MAXFP-Miner:利用FP-tree快速挖掘最大频繁项集[J].控制与决策,2005,08:887-891.
[16]何晓健.视网膜静脉阻塞家兔血管内皮生长因子的表达研究[D].浙江大学,眼科学,2004,硕士.
[17]贾鹏宇,李艳,郑琼林.电压型级联系统中减小源变换器输出阻抗的有源阻尼控制方法[J].电工技术学报,2015,08:71-82.
[18]徐莉泉.短程生物脱氮系统污泥基团内部氮素迁移转化微生态特性试验研究[D].西安建筑科技大学,环境工程,2013,硕士.
[19]张阳.基于产品评论的观点挖掘研究[D].西北大学,计算机应用技术,2014,硕士.
[20]刘赛特.社会工作介入流动儿童服务研究[D].吉林大学,社会工作,2013,硕士.
[21]曾坚.小麦番茄红素环化酶基因TaLCYB的RNAi载体构建及遗传转化研究[D].华中科技大学,生物化学与分子生物学,2013,硕士.
[22]袁占亭,张秋余,李威.数据抽取及语义分析在Web数据挖掘中的应用[J].计算机工程与设计,2005,06:1425-1427+1437.
[23]郝建英.高中生数学抽象概括能力的调查与研究[D].河北师范大学,学科教学,2014,硕士.
[24]贺振明.高速铣削在模具制造中的应用[J].机电一体化.1998(04)
[25]李红兵,佟东.论中国港口业发展现状及其在国民经济发展中的作用[J].山东社会科学,2014,01:144-149.
[26]苑芳.提高我国私有制经济中劳动者报酬研究[D].南京财经大学,政治经济学,2012,硕士.
[27]张小川.基于融资租赁视角的中小企业融资问题研究[D].安徽大学,金融,2014,硕士.
[28]万宁.方钢管混凝土受压短柱力学性能的分析[D].西南交通大学,结构工程,2012,硕士.
[29]陶旭婷.基于局部特征的多光谱与全色图像融合算法研究[D].西南交通大学,信号与信息处理,2014,硕士.
[30]牛妞.人力资本投资理论视域下研究生就业问题研究[D].沈阳师范大学,课程与教学论,2014,硕士.
[31]江平,常军霞,朱强,何代平.基于石墨烯/金纳米复合材料的无酶葡萄糖传感器[J].分析科学学报,2014,03:373-376.
[32]孙龙.潜孔锤及声波钻柔体振动冲击的理论与建模[D].中国地质大学（北京）,机械工程,2014,硕士.
[33]莫尊理,王博,赵国平,朱小波,郭瑞斌.石墨烯/碳纳米管复合薄膜的制备进展[J].应用化学,2014,08:871-877.
[34]李向坤.基于认知策略训练的大学英语听力教学实证研究[D].哈尔滨工业大学,外国语言学及应用语言学,2014,硕士.
[35]战德臣,聂兰顺,徐晓飞.“大学计算机”——所有大学生都应学习的一门计算思维基础教育课程[J].中国大学教学,2011,04:15-20.
[36]张黎.滋水清肝饮加减治疗肾阴虚肝郁型围绝经期综合征临床观察[D].成都中医药大学,中医妇科学,2012,硕士.
[37]彭博.成都晚报市场定位与营销组合策略研究[D].四川大学,2002.
[38]谭寅寅.高等农业院校校园文化建设研究[D].南京农业大学,思想政治教育,2013,硕士.
[39]何晓琴.提高中职学生手工焊接兴趣的策略[J].电子制作,2015,13:192.
[40]尹应凯.中国关税有效保护效果研究[D].华东师范大学,国际贸易学,2002,硕士.
[41]黄年华.荀子管理伦理思想探究[D].浙江财经学院,伦理学,2012,硕士.
[42]杨凯.公共部门个体绩效考核研究[D].厦门大学,公共管理,2014,硕士.
[43]刘刚.说话人识别系统设计研究[D].华东师范大学,计算机应用技术,2004,硕士.
[44]兰宁.儿童卫浴产品情感化设计研究[D].中央美术学院,艺术设计,2013,硕士.
[45]童玲.基于SPC的高等教育教学质量评价方法研究[D].西安电子科技大学,计算机应用技术,2011,硕士.
[46]韩再娥.内蒙古自治区高校高水平排球队伍建设现状分析及对策研究[D].内蒙古师范大学,运动训练（专业学位）,2013,硕士.
[47]王康,葛晓慧,辛焕海,甘德强.基于牛顿—拉夫逊法的光伏面板模型参数计算[J].电力系统自动化,2010,20:94-98+107.
[48]周凯.直角坐标式弯丝机器人的结构设计与运动学分析[D].哈尔滨理工大学,机械电子工程,2013,硕士.
[49]李世雨.热痹外用方巴布剂治疗热痹疗效观察[D].北京中医药大学,中医学,2013,硕士.
[50]奴尔加依娜提.B2C环境下大学生网站信任的前因及对购买意向的影响[D].浙江大学,2008.