发电机组的优化配置问题数学建模论文

1、院系联系电话题目发电机组的优化配置摘要本文针对不同种类发电机，不同时段的实际用电情况，建立了如何合理安排发电机使用的模型。对于问题（一），该模型灵活运用二次规划，整体考虑一天中的各个阶段，并利用lingo求得一天中最小费用为997790（元）。在问题（二）中，应用经济学模型和统计学中线性回归分析的原理，并利用excel中丰富的统计函数和lingo软件求得结果。问题（三），仿照问题（一）的方法，但发现最小费用没变。正文一、问题重述电是我们生活中不可缺少的一部分，现考虑发电机组优化配置问题。某发电厂负责某地区的供电任务，已知该地区夏季一天的电力需求如下：段（时）066-99-1515-1818-2

2、4需求（兆瓦）1500030000250004000027000现电厂有三种类型发电机可投入运转：一型12台；二型6台；三型5台；各个型号机组相关数据如下：输出最低水平（兆瓦）输出最高水平（兆瓦）输出最低水平费用（元）每台机组开动准备费用（元）（每小时）（每兆瓦）1850200010002200021250175026001.3100031500400030003500正常情况下，在满足估计的负荷要求之外，每一时刻运转的发电机组应足够多，使得当负荷增加不超过15%时，能够通过调高运转的发电机组的输出来满足增载的要求。请你建立该问题的数学模型，通过求解模型回答以下问题：（1）在一天中各个时间段应

3、安排使用那些发电机组运转可以使得在满足负荷要求的情况下总的费用最低？总的费用为多少？（2）在一天中每段时间，电力生产的边际费用是多少？即应为用电定什么价格？（3）将后备输出保障的15%降低为10%,运转费用节省的情况如何？可以降为多少？二、问题的基本假设1. .假设每个阶段不会变更设备。2. 不考虑设备需要维护与修理。3. 假设用电需求相对稳定，不会发生突变。4. 关闭和启动发电机时均是瞬时完成，不计相应使用的时间5. 发电机输出过程其功率始终保持不变。三、符号说明名称符号阶段i不性j发电机的台数n（i,j）每台每小时实际的电量pu,j）第j台单位时间最低电量minpj第j台单位时间最高电量m

4、axp名称符号第i阶段用时ti第i阶段用电需求w第j种启动费bj第j种输出最低水平费用a第j种每小时超出的费用fj第j种最多台数njf总费用Q第j阶段边际费用Cj边际成本MC四、模型建立与求解问题(一)模型分析此题研究的是每天在第j个阶段如何合理分配各种类型的发电机，使每日的成本Q最低的数学模型，建立如下二次规划模型3。发电机组运转的费用可以分成两部分：第一部分与其输出功率和运行时间紧密相连，即：53工工(P(i,j)minpj)*fi+aj)*n(i,j)*t；i=4j4第二部分与其输出功率和运行时间无关，即开机时费用其中工x,x0g(x)=s,0,x三0还有当i=1时，n(i-1,j)=n

5、(5,j)。模型的建立确定目标函数在满足需求量的情况下，为了使每天发电成本最低，我们建立如下目标函数53工工(P(i,j)-minPj)*fi+aj)*n(i,j)*+bj*g(n(i,j)-n(Tj)其中i=1jT约束条件模型的求解通过lingo编程2(见附录)得到结果Q为997790.0元在不同时间段对各类型发电机所需求的数量，功率如下表所示：段种类台数实际功率台数实际功率台数实际功率台数实际功率台数实际功率一12850121625121208.333121958.333121375二3160061750617506175061750三01502.593015000150041500015

6、01.593结果分析由上述表格可知，第一种机型在所有阶段里全部使用；第二种机型在各阶段几乎全部使用；第三种机型则几乎没用。经分析，考虑到需求量与使用成本的关系，第三种机型不适合用来在该地区供电。问题（二）模型的分析此题要求一天中每段时间电力生产的边际费用（即边际成本）。首先给出边际成本1的定义：MQM厂商在短期内增加一单位产量是所增加的成本。经过lingo试验，”Q?30时，发电机组的类型和数量都保持不变。所以启动费不变，不予考虑。模型的建立由问题（一）可知不能直接求出各阶段用电价格与需求电量之间的关系，所以不能直接求导。只能通过St0逐渐逼近才能求出边际成本。并用线性回归分析MC与?Q的关系

7、。MC=B0+Bi?Q模型的求解本题分别设？Q=土7,土10,±16,土21,土27兆瓦；通过运行lingo（程序见附录）和excel中的INTERCEPT函数求解拟合直线的截距和LINEST函数求解拟合直线的斜率可获得各阶段边际成本价格MC（Q）（表格见附录）B0B1y12.24825E-05y25.99852711.57778E-07y312.0016877-5.6073E-06y46.00085617.46349E-07y511.99706824.44444E-09模型的结果由上表格储相对于Bo为零，每个阶段边际成本逼近于一个常数，且阶段三与阶段五的拟合曲线相近，阶段二与阶段四

8、相近。结果见下表：阶段一一三四五边际费用（元/mvv7.7995.99912.0026.00111.997单位时间边际费用（元/（MW力）1.3002.0002.0002.0002.000结果分析用影子价格的方法对上述结果进行检验，可以看出结果正确。问题（三）模型分析这里的目标函数同问题一相同，建立单目标最优模型。只需在约束条件中将后备输出保障的15%降低为10%。目标函数53工工（P（i,j）-minpj）*fi+a。*n（i,j）*ti+bj*g（n（i,j）-n（j,j）其中iTj与约束条件模型的结果通过lingo编程（见附录）得到结果Q为997790.0元在不同时间段对各类型发电机所需

9、求的数量，功率如下表所示:阶段一二三四五种台实际功率台实际台实际功率台实际功率台实际功率类数数功率数数数一12850121625121208.333121958.333121375二3160061750617506175061750三01502.5930150001500415000P1501.593结果分析由上表可知将后备输出保障的15%降低为10%后，总费用Q不变,且各阶段每种机型的使用台数也不变。五、误差分析发电机输出过程其功率始终保持不变是在理想情况下才存在的，有一定的误差，还有在开启以及关闭发电机时不仅有功率损耗方面的误差，而且还有时间方面的误差，这些都会对模型的建立有一定程度上的影

10、响。此外lingo的运算结果只保留小数点后3位数字，所以也存在无误差。六、模型评价模型的优点：(1)模型的结构简洁明了，可移植性好；(2)问题(一)的模型可操作性强，适用范围广；(3)本模型中使用了lingo和excel软件，大大减少了计算量。模型的缺点：(1)由于软件等原因，存在误差；(2)问题(二)模型的操作较为复杂。七、参考文献【1】高鸿业.西方经济学.中国人民大学出版社.【2】宋来忠.数学建模与实验.科学出版社.【3】姜启源，谢金星，叶俊.数学建模.高等教育出版社.【4】蔚诗松，程依明，濮晓龙.概率论与数理统计教程.高等教育出版社.【5】胡运权.运筹学基础及应用.高等教育出版社.附录程

11、序(问题一)model:sets:Time_slot/1.5/:t,w;Generator_type/1.3/:a,f,b,minp,maxp;Link(Time_slot,Generator_type):n,p;endsetsdata:t=6,3,6,3,6;w=15000,30000,25000,40000,27000;minp=850,1250,1500;maxp=2000,1750,4000;a=1000,2600,3000;f=2,1.3,3;b=2000,1000,500;enddatamin=sum(Link(i,j):n(i,j)*a(j)*t(i)+n(i,j)*(p(i,j

12、)-minp(j)*t(i)*f(j)+if(n(if(i#eq#5,1,i+1),j)#gt#n(i,j),n(if(i#eq#5,1,i+1),j)-n(i,j),0)*b(j);for(Time_slot(i):,j)>=w(i););sum(Generator_type(j):n(i,j)*p(isum(Generator_type(j):n(i,j)*p(ifor(Time_slot(i):,j)<=w(i)*1.15;);for(Time_slot(i):for(Time_slot(i):for(Time_slot(i):bnd(0,n(i,1),12););bnd(0

13、,n(i,2),6););bnd(0,n(i,3),5););for(Link(i,j):bnd(minp(j),p(i,j),maxp(j););for(Link(i,j):gin(n(i,j););end程序（问题二）以。Q=7为例：model:setsTime_slot/1.5/:t,w,q;Generator_type/1.3/:a,f,b,minp,maxp;Link(Time_slot,Generator_type):n,p;endsetsdata:t=6,3,6,3,6;q=7?;w=15000,30000,25000,40000,27000;minp=850,1250,1500

14、;maxp=2000,1750,4000;a=1000,2600,3000;f=2,1.3,3;b=2000,1000,500;enddatamin=sum(Link(i,j):n(i,j)*a(j)*t(i)+n(i,j)*(p(i,j)-minp(j)*t(i)*f(j)+if(n(if(i#eq#5,1,i+1),j)#gt#n(i,j),n(if(i#eq#for(Time_slot(i):,j)>=w(i)+q;);for(Time_slot(i):,j)<=(w(i)+q)*1.15;);for(Time_slot(i):for(Time_slot(i):for(Tim

15、e_slot(i):5,1,i+1),j)-n(i,j),0)*b(j);sum(Generator_type(j):n(i,j)*p(isum(Generator_type(j):n(i,j)*p(ibnd(0,n(i,1),12););bnd(0,n(i,2),6););bnd(0,n(i,3),5););for(Link(i,j):bnd(minp(j),p(i,j),maxp(j););for(Link(i,j):gin(n(i,j););end程序（问题三）model:setsTime_slot/1.5/:t,w;Generator_type/1.3/:a,f,b,minp,maxp

16、;Link(Time_slot,Generator_type):n,p;endsetsdata:t=6,3,6,3,6;w=15000,30000,25000,40000,27000;minp=850,1250,1500;maxp=2000,1750,4000;a=1000,2600,3000;f=2,1.3,3;b=2000,1000,500;enddatamin=sum(Link(i,j):n(i,j)*a(j)*t(i)+n(i,j)*(p(i,j)-minp(j)*t(i)*f(j)+if(n(if(i#eq#5,1,i+1),j)#gt#n(i,j),n(if(i#eq#for(Ti

17、me_slot(i):5,1,i+1),j)-n(i,j),0)*b(j);sum(Generator_type(j):n(i,j)*p(ifor(Time_slot(i):,j)>=w(i););sum(Generator_type(j):n(i,j)*p(i,j)<=w(i)*1.10;);for(Time_slot(i):bnd(0,n(i,1),12););for(Time_slot(i):bnd(0,n(i,2),6););for(Link(i,j):for(Link(i,j):endbnd(minp(j),p(i,j),maxp(j););gin(n(i,j););表格（问题二）阶段-27-21.000-16-10-7710162127