电力发展规划模型

电力发展规划模型某地区未来十年电力发展规划模型【摘要】本文建立了一个电力规划中经济最优化静态模型。本文认为该问题中的经济最优化不是经济利益达到最大化，而是在满足技术要求下的成本最小化。该模型首先利用数学软件LINDO解出满足技术要求下的各方案的分配份额：扩建旧火电站:新建水电站:新建火电站=5:4:2。其次，本文考虑了各项建设的时间、顺序上的安排对于节省成本的影响，用MATLAB软件通过画图确定了在第一年年初即队原旧火电站进行扩建；在第四年投入使用新的水电站；最后在第八年年初起投入使用新的火电站。再次，本文考虑了从装机台数的角度降低成本，以期在最小的成本投入下完成各项建设。我们最后得出的建设方案安排为：在第一阶段，也即扩建原有的旧火电站阶段，我们预备在第一年扩建两台单机容量为10千瓦的发电机组；而在第二年和第三年各扩建一台和两台；第二阶段，也即新建水电站阶段，从第四年初到第七年末的四年间，每年建设一台单机容量为25千瓦的发电机组台数；第三阶段，也即新建火电站阶段，在第八年和第九年初，各建设并投入使用一台单机容量为30千瓦的发电机组；至此，全部建设完成并投入正常运作。由此，我们算出最终的投入成本为23.3866亿。【关键词】经济最优化、成本最小化、静态模型、MATLAB一、问题的重述某地区在制定十年电力发展规划时遇到一个问题：根据电力需求预测得知，该地区在十年后发电装机容量需要增加180万千瓦，那时的年发电量需要增加100亿度。根据调查和讨论，电力规划的备选技术方案有三个：1）扩建原有的火电站，但最多只能再安装五台10万千瓦的发电机组；2）新建水电站，但最多只能安装四台25万千瓦的发电机组；3）再新建一个火电站，最多只能安装四台30万千瓦的发电机组。通过调研和计算，获得有关的参数如表1所示：表1:备选技术方案参数表表中负荷因子为全年满功率运行天数与全年总天数之比。根据该地实际调查原有火电站平均全年满功率运行天数为241天，水电站和新建的火电站应分别为146天和255天，而全年365天，故折算得表中数据。表中资本回收因子是由如下数据所确定的，火电站的回收年限取15年，年利息率为0．06；水电站的回收年限取30年，年利息率为0．04，即得表中所列数值。我们的任务是：在满足上述技术要求的前提下，选取经济效果最优的建设方案.即：1.原来的火电站应如何扩建？2.新建的火电站应如何确定单机容量为25万千瓦的发电机组的数量？3.新建的火力发电站应如何确定单机容量为30万千瓦的发电机组的数量？二、问题的分析与假设因为每台发电机在开始工作之后，每年必然会有所盈利，且由于在上限的限制上我们只有允许最大装机台数的限制，若如此，必然是在全部应用时是最优的，所以我们在确定这个模型的经济最优解时，我们的目标函数不应是收入与支出的顺差的最大值，而是将目标函数定为满足技术要求前提下，成本付出达到最小值。为此，我们做出了以下假设：1、发电机都是在每年的年初投入使用的，且每年的发电机都是满负荷工作的，即负荷因子不变；2、不计每年的发电剩余量的价值，即除去本地区的用电需求之外，无跨地区的输电；3、十年中所有的发电机都是正常运作，且不会出现拉闸限电现象，每年所发电量满足每年各月的用电需求量；4、现阶段的发电量刚好满足生产生活所需的用电量；5、把职工的薪金，税收等额外付费每年当作一固定开销支付，故其对发电所得收益无波动影响；三、符号说明i：1，2，3分别表示扩建的旧火电站，新建的水电站，以及新建的火电站；T：时间，单位：年；N：火电站和水电站最终所装机的台数；iK：各新增的发电站的工程投资资金，单位：万元；iy：各发电站的年运行成本，单位：万元；iY：各发电站的年运行总成本，单位：万元；iα：该地区的用电量的年平均增长率；β：各发电站的贷款年利息率；iθ：火电站或水电站的负荷因子；ini：水电站或火电站的回收期限，单位：年；A：各发电站的每年等额支付现金，单位：万元；iP：现阶段需求电量，也即现阶段发电能力，单位：万度；oR：发电站新增发电量所需支出，单位：万元；O：经资本回收折算后投资的每年应支出量，单位：万元；iC：各新增的发电站的单机容量，单位：万千瓦；iiX：各新增的发电站的装机台数；ir：各新增的发电站的资本回收因子，单位：百万元/亿度；iTQ：某年用电需求量，单位：万度；iF：扩建的或增加的发电站的年发电量，单位：万度；TW：在T年电的总供应量；单位：万度；iV：扩建的或增加的发电站的装机容量，单位：万度；四、模型的建立与求解4.1建立一次性投入模型在该模型的建立中，我们首先考虑最简单的情况：第一年的投资就满足十年以后的用电需求量。因为每年的支出情况是一个定值，我们就可以考虑转化成每年的支出情况。4.1.1得到折算后投资的年支出量iO根据经济学知识可知，现值p在回收期限为n年，年利率为β情况下，其资本回收因子r有如下公式：()()111nnrβββ*+=+-则等额系列偿还资金有如下计算公式：()()111nnAprpβββ*+=*=*+-类似的，对于对应的扩建的旧火电站和新建的水电站、火电站也有以下公式表示它们的等额偿还资金：()()111niiiiiniiApβββ*+=*+-4.1.2发电站投入使用之后，各发电站的年运行总成本iY经分析，年运行总成本iY是与该发电站的年发电量iF成正比关系的，而年发电量iF则取决于火电站和水电站的装机台数iN、单机容量iC和负荷因子iθ。有如下关系式：24365iiiiiiYyFNyC=*=****新的水电站和火电站的投资都有个前期工程投资，但是因为是个固定的常数，所以它所引起的等额偿还资金也是个常数，我们这里在考虑成本最小化时可以忽略它的影响。所以，根据题目的技术要求，我们可以将目标函数表示如下：3311iiiiMinRAY===+∑∑但是我们是有所约束的，为非限制性约束：1231122331231000000180544FFFCNCNCNNNN++≥??*+*+≥??≤??≤??≤?用LINDO解得各个相应的值为：123121140.3542RNNN=??=??=??=?具体实现程序见附录一。4.2考虑该地区用电需求量以递增的趋势增长的模型4.2.1确定用电年增长量iQ和年需求增长率α因为实际生活中，人们的生产生活对用电需求量会因为经济水平的提高，恩格尔系数的降低而不断的增长，且我们可预见，随着时间的不断后延，我们的科技水平也在不断的进步，对用电的需求增量也会不断增大，在短时期内，用电增长量应该满足一个规律，故大胆假设：这个用电需求量的增长曲线应该是一个指数增长型的，有如下公式：()01nQpα=*+将公式作如下变形：()101001loglogQnpα??=*+???其中，Q代表现阶段用电需求量，0p代表原来的电力需求量，α代表年平均需求增长率，n代表年限，即第几年。我们利用网络资源找到了一组全国的从1981年至2000年的历年电力装机和发电量的构成比的数据，如下表：历年电力装机和发电量的构成比（1981～2000）我们假定上面这组数据中的电的需求和供应量始终是平衡的，这样我们就可以从标了标记的两组数据的和的变化去求得这个年增长率。在对数据进行了前期取对数之后，我们利用MATLAB软件对数据进行了线性拟合，得到我们所需要的数据。并用MATLAB对这20个数据进行了数据线性拟合（具体实现过程见附录二）。02468101214161844.14.2年份用电年需求量增率的对数与时间的关系图一如上图所示，数据的线性情况较好，故我们认为假设合理，于是我们有：()1010.035logα+=∴08.39%80.8Pα=??=?由于现阶段电量的供求刚好达到平衡，且我们的目标函数是求所花成本投入最小，显而易见，扩建旧火电站是必须先做的，这样才有可能达到成本最小值。所以我们接下来要解决的就是新的水电站和新的火电站的建造顺序和建造时间问题，决定水电站和火电站的建造时间和顺序，以求得在最小的最终建设成本条件下完成建设。4.2.2解出水电站和火电站的建造顺序和建造时间iT正如我们的模型的分析中所述，每台发电机投入生产之后，其每年都在盈利，经计算比较，其金额比发电站的年支出总资金大一个数量级，故直接考虑顺差的差值大小来决定先建造哪一个发电站是不合理的，我们对于建造顺序的安排也取决于如何使所花成本最小。由于每年的支付金额都是等额的，而且建成发电站之后全部在负荷因子下工作，所以各个发电站的年运行成本也是一个固定的数值，所以我们可以考虑用一年的支出总金额来代替十年支出总金额。在一年的费用支出中，主要有两方面的因素：1）每年运行时的年运行成本；2）购买机器时的工程投入所引起的等额偿还资金。在上一模型中，我们已经知道某一发电站的年运行总成本iY是与该发电站的年发电量iF成正比关系的，年发电量iF又取决于火电站和水电站的装机台数iN、单机容量iC和负荷因子iθ。工程投资所引起的等额偿还资金也可由上一模型中的公式可以算得。故水电站或火电站的年费用支出资金应为：iiiiFyKθ*+*经计算，水电站的年费用支出资金要小于火电站的费用支出，所以我们判定三个发电站的最终建造顺序方案为：先扩建旧火电站，再新建水电站，最后新建火电站。在模型中，我们可以预见用户对于用电的年需求量是一个指数上升的曲线，而电的供应量是一个分段的函数模式。()00101201231nTQpaPFWPFFPFFF?=*+?+????=++????+++??预测电需求量拟电供应量根据最小成本原则，我们可以解出在什么年份时原有的电力供应会出现供不应求的现象，然后在这一年到来之际，增加投入建新的发电站。显然，这样的投资将会使成本最小，故我们可依据以下平衡关系式解出相应的年份iT：QW=再由上面模型的一些公式和结论，解得：233.77.2nn=??=?∴[][]12233037floorfloorTTnTn?=??==??==??即在当年扩建旧火电站之后，在第3年开始建造并使用新建的水电站，第七年建造使用新的火电站。利用MATLAB作图如下（具体程序实现见附录三）：01234567891080100120140160180200年份用电需求量电的供需量与时间的关系图二4.3考虑装机台数的时间模型在上图中，我们看到在需求曲线上方的分段函数中，有部分电是处于浪费状态的，这必将使投入的资金不能产生最大的经济效益，也间接的体现出成本还未达到最小，还有更为节省的投资方案。我们是否可以考虑将机器的投入生产采取分段投入的方案呢？这个方案就是使电供应的分段函数更为靠近需求曲线，也即两个函数与坐标轴围住的面积达到最小，我们提出的解决方案如下：算出每年的需求量，在分段函数上取尽可能小的函数曲线与之对应（程序见附录四），列表如下：我们考虑在满足用电年需求量的技术前提下，使供应量达到最小，使成本最小化。由上面的图4与表3的数据及其题目给我们的各发电站的单机年发电量，我们可以进而求得在T年时满足需求量下的最小装机台数（程序见附录五），列表如下：现利用表四用MATLAB作出其图形如下：01234567891080100120140160180200年份用电需求量电的供需量与时间的关系图三4.4求出最终的成本由于每年的发电机投入都是不同的，故应该计算十年的总资本投入。因为每台发电机的投入资金量对另外的发电机没有影响，所以可以独立地进行计算。总成本由两部分构成：工程投入偿还资金和年运行成本投入。3131iiiiiiiiiiTrNCTyNF==?***????***??∑∑工程投入偿还资金年运行成本经计算，最终成本为23.386634亿元。五、模型的优缺点及其改进方向模型的优点：1）本模型在考虑最经济的方案时，并没有直接考虑其收益的最大值，而是取其成本的最小值；2）本模型将复杂的动态规划问题通过合理的假设简化成静态的规划问题；3）本模型对问题以点连线，步步深入。但并没有什么晦涩难懂的地方，易于理解；模型的缺点：1）本模型在考虑发电能力时，没有考虑电力生产弹性系数；2）模型对电站建设进行了理想化的简化假设，在每年初投入使用，且一有需要就能马上提供，没有考虑施工进程的限制不太符合实际，在实际的发电站建设中应投入更大。改进方向：1）现在企业都是向环保型企业发展了，社会基本设施也不能对环境拥有特权，考虑到火电站的发电原料主要是煤，里面含有硫物质，会对大气环境产生比较大的污染，如果不加以脱硫处理的话；而且煤是我国最重要的燃料资源之一，煤的储存量正在减少。而水电站是完全绿色型的，且是可再生资源。故在发电站的新建上应优先考虑水电站；2）我们把用电需求量只是简单地处理成一个指数函数，实际情况远非如此简单，从长期时间段来看，应该近似满足Logestic曲线，所以我们可以向着中长期方向发展，以求有更广的适用范围，更好的预测效果；3）我们的假设中提出无跨地区的输电，在现实生活中却存在着“西电东送”的现象，我们将考虑跨地区的输电以减少浪费值，可以弥补成本；4）由于现实生活中可以有拉闸限电的情况存在，我们可以在装机过程中使那条供应的折线图略小于需求曲线，从生活用电量的节俭来降低成本。六、参考文献[1]、吴添组工程经济（第二版）高等教育出版社2004[2]、蔡锁章数学建模原理及方法海洋出版社2000[3]、姜启源数学模型（第三版）高等教育出版社2003[4]、网址：/newsp/zgdl/dltj/gndltj26.htm/newsp/zgdl/dltj/gndltj23.htm/newsp/zgdl/dltj/gndltj27.htm七、附录附录一：用LINDO解线性规划问题：min5957.52N1+5063.28N2+18910.8N3+2377.22N1+1997.28N2+6701.4N3st57840N1+87600N2+183600N3>=100000010N1+25N2+30N3>=180N1<=5N2<=4N3<=4endgin3附录二：X中的值表示年需求量增率的对数>>x=[3.490379923.5154764413.5459253293.576341353.6135247033.6528263033.6966184593