电力生产问题数学建模论文设计

wordword/word电力生产问题摘要本文针对发电机厂每天在不同时间段用电需求量不同的情况下，根据给定不同型号不同数量的发电机，合理分配各台发电机在不同时间段的开启数量和运行功率，使得一天总发电本钱最小的问题，采用单目标非线性规划方法，建立所求问题的最优化模型，借助Lingo软件对模型进展求解，得到每日最小发电总本钱，以此制定机组启停计划。针对问题一：为了使一天总的发电本钱最低，同时还要考虑到不同时间段开机数量不同对启动本钱的相互影响，将七个时间段的本钱统一考虑，选取相应的约束条件对目标函数进展约束，从而给出优化模型，运用非线性规划的方法，利用Lingo编程求解，得到发电厂每天最小发电总本钱为：1463625元。针对问题二：根据题目要求，在任何时刻，正在工作的发电机组必须留出20%的发电能力余量，以防用电量突然上升，在建模时将每台发电机的实际输出功率降至80%，所以可以按照问题立的模型，将其约束条件中每个时间段用电量的需求量进展相应的修改即可，同样利用Lingo编程求解，得到发电厂每天最小发电总本钱为：1885420元。最后，对最终得到的结果进展分析，对模型的优缺点给出了客观的评价，就模型的不足之处提出了改良方法，并对模型的应用以与前景提出了自己见解。关键词：电力生产优化非线性规划多变量整体最优化1问题重述为了满足每日用电需求，有四种类型的发电机可供选择。用户每日的用电主要分为7个阶段，每个阶段的用电需求〔单位为兆瓦〔MW〕〕各不一样，为了能够高效低本钱完成每天发电计划，就必须使得每阶段的供需平衡，否如此就会影响电力系统的安全运行。为了能够实现这样的平衡状态，就需要电力部门对发电机组进展合理的启停计划，在满足每日用电需求的前提下，追求发电本钱的最小化。在不考虑其它本钱因素的前提下，假定所有发电机组的发电本钱都是由三局部组成：1.固定本钱；2.边际本钱；3.启动本钱。需要考虑的约束有：1.负荷平衡约束；2.发电机组使用数量围约束；3.发电机组输出功率围约束；4.每日电力需求约束。因此，在不同时段开启哪些型号发电机，使发电厂每天的发电总本钱最小是一个值得研究的问题。每天的用电需求分为七个阶段，各个阶段的用电需求各不一样，其用电情况如表1。表1：每日用电需求〔兆瓦〕时段〔0-24〕0-66-99-1212-1414-1818-2222-24需求12000320002500036000250003000018000为了便于观察每天的用电需求，将数据重新整理，转化为图1所示的图表。图1各时段的用电需求量从图表中可以清晰的观察到每天用电需求变化，在第一阶段用电需求处于低谷时段，第四阶段处于峰值时段，且用电需求变化较大。可供选择的发电机共有四种型号，每种发电机都有一个最大发电能力，当接入电网时，其输出功率不应低于某一最小输出功率。所有发电机都存在一个启动本钱，以与工作于最小功率状态时的固定本钱，并且如果功率高于最小功率，如此超出局部的功率每兆瓦每小时还存在一个本钱，即边际本钱。发电机的各项参数指标如表2。表2：发电机情况可用数量最小输出功率〔MW〕最大输出功率〔MW〕固定本钱〔元/小时〕每兆瓦边际本钱〔元/小时〕启动本钱型号110750175022505000型号241000150018001600型号381200200037502400型号431800350048001200只有在每个时段开始时才允许启动或关闭发电机，而且关闭发电机不需要付出任何代价。1.3需要解决的问题：问题〔1〕在每个时段应分别使用哪些发电机才能使每天的总本钱最小，最小总本钱为多少？问题〔2〕如果在任何时刻，正在工作的发电机组必须留出20%的发电能力余量，以防用电量突然上升。那么每个时段又应分别使用哪些发电机才能使每天的总本钱最小，此时最小总本钱又为多少？2模型假设与符号说明假设1：每台发电机在同一时间段按照预订功率稳定运行，并且输出功率恒定不变。假设2：发电机在相邻时间段进展功率调整时是瞬间完成的，并且不需要考虑调整本钱。假设3：发电机的启动和关闭是在瞬间完成，并启动之后就可以达到设定工作功率。假设4：发电机在工作过程中不考虑其电能的损失，即实际输出电能全部转化为用户需求。假设5：发电厂每日发电总本钱仅由发电机组的固定本钱、边际本钱和启动本钱构成，不考虑其他本钱的影响。假设6：在第一时间段开机前，所有的机组都处于关闭状态。表3：符号说明符号符号说明第j种型号发电机在第i阶段开启的数量第j种型号发电机在第i阶段的输出功率第i个时间段的总时间第j种发电机的每小时边际本钱第j种发电机的每小时固定本钱第j种发电机的启动本钱第j种电机的最小输出功率目标函数，每天发电机组发电的总本钱第时段发电机组用户需求量与实际输出总功率的差值第时间段用户的电量需求3问题分析多机组启停优化问题是在满足约束的前提下，优化确定每个阶段机组的启停，求出机组的最优运行方案，实现每日发电总本钱最小。首先建立每日发电本钱的目标函数和约束条件的数学表达式。总本钱是由各个时间段的总固定本钱、总边际本钱和总启动本钱构成，因此就可以根据的数据，求出相应本钱表达式。其中启动机组需要相应的启动费用，而关闭机组如此不需要费用，这样上一阶段的电机运行情况将影响下一阶段的启动本钱，进而影响总本钱，因此在考虑电机的运行费用时应该把下一阶段电机的运行状况和上一阶段的运行状况联系起来，故在此需要对全天的7个时间段的发电机的开启情况进展统一、合理的安排，不可只考虑某一个时间段。在第二问中如果在任何时刻，正在工作的发电机组必须留出20%的发电能力余量，以防用电量突然上升。在求解时每台发电机的输出功率就应该改变其相应的约束，将实际输出功率减少20%，即将每日的用电需求提高25%，然后再对模型进展求解。机组组合问题是一个多变量、多约束的混合非线性规划问题，因此在求解时需要对各个时段每一台用于发电的发电机所需要的各项本钱进展求和计算，在此我们采用Lingo软件对其进展求解，得到发电厂每日最小发电总本钱。4模型建立与求解由题目给出的条件与模型的假设知：发电厂每日发电总本钱仅由发电机组的固定本钱、边际本钱和启动本钱构成。〔1〕每天四种型号发电机固定本钱：〔2〕每天四种型号发电机边际本钱：〔3〕每天四种型号发电机总启动本钱：〔4〕目标函数为：综上所述，得到问题一的多变量最优化模型：其中量〔各时段需求电量〕。由上述分析可知，该问题为多变量非线性规划问题，应用LINGO程序进展编程〔程序见附录〕计算，最终得出每时段各型号发电机的使用数量与其各自的功率。由各型号发电机使用数量与各自功率可求出各时段的最小本钱与一天的最小总本钱，得到发电厂每天最小发电总本钱为：1463625元。具体的数据见表4。表4：问题一最优化方案时段0—66—99—1212—1414—1818—2222—24型号1台数0333220功率0750175075013000型号2台数4444444功率150015001500142515001500型号3台数3888886功率2000200020002000200020002000型号4台数0313130功率018001800180018000由题目给出的条件与模型的假设知：发电厂每日发电总本钱仅由发电机组的固定本钱、边际本钱和启动本钱构成。〔1〕每天四种型号发电机固定本钱：〔2〕每天四种型号发电机边际本钱：〔3〕每天四种型号发电机总启动本钱：，〔4〕目标函数为：综上所述，得到问题一的多变量最优化模型：与问题一相比拟，只需针对每台实际输出功率进展改变即可，将实际输出功率减少20%，即将N提高25%，如此：由上述分析可知，该问题为多变量非线性规划问题，应用LINGO程序进展编程〔程序见附录〕计算，最终得出每时段各型号发电机的使用数量与其各自的功率。由各型号发电机使用数量与各自功率可求出各时段的最小本钱与一天的最小总本钱，得到发电厂每天最小发电总本钱为：1885420元。具体的数据见表5。表5问题二最优化方案时段0—66—99—1212—1414—1818—2222—24型号1台数1778661功率1000175017507501683750型号2台数4444444功率15001500150015001500型号3台数4888888功率2000200020002000200020002000型号4台数0323330功率0180030001800180005结果分析将问题一机组启动计划最优化方案转化为图示2以便于直接观察：图2各时段各型号发电机台数安排由图2可清晰看出各型号电机在不同时段的使用情况。其中，型号1的发电机组在每天的各个时段的使用数量均维持在3台与以下；型号2的发电机组在每天的各个时段的使用数量均维持在4台；型号3的发电机组在每天的各个时段的使用数量均维持在8台与以下；型号4的发电机组在每天的各个时段的使用数量均维持在3台与以下。由此可见，型号3的发电机组的使用频率较高。图3各时段各型号发电机输出功率由图3可清晰看出各型号电机在不同时段的输出功率情况。其中，型号1的发电机组在每天的各个时段使用的输出功率均维持在750~1750MW之间，且多数时段的输出功率为750MW或1750MW；型号2的发电机组在每天使用时除第七时段的输出功率除第三、第五时段,其他时段的输出功率根本维持在1500MW；型号3的发电机组在每天使用输出功率均为2000MW；型号4的发电机组在每天使用时输出功率维持在1800~3000MW之间。5.2问题二结果分析将最优化方案图示如下：图4各时段各型号发电机台数安排由图4可清晰看出各型号电机在不同时段的使用情况。其中，型号1的发电机组在每天的各个时段的使用数量均维持在8台与以下；型号2的发电机组在每天的各个时段的使用数量均维持在4台；型号3的发电机组在每天的各个时段的使用数量均维持在8台与以下；型号4的发电机组在每天的各个时段的使用数量均维持在3台与以下。图5各时段各型号发电机输出功率由图5可清晰看出各型号电机在不同时段的输出功率情况。其中，型号1的发电机组在每天的各个时段使用的输出功率均维持在750~1750MW之间；型号2的发电机组在每天的各个时段使用的输出功率均维持在1337.5~1500MW之间，且多时段为1500MW；型号3的发电机组在每天使用时各时段输出功率均为2000MW；型号4的发电机组在每天使用时输出功率均维持在1800~3000MW之间。6模型的评价与改良模型的优点：〔1〕我们将一天所有时段作为整体出发，考虑到了各个时段间启动发电机对本钱的影响，构建了不同时段发电机组的启动本钱计算公式，该公式可以很好的满足出不同时段发电机组的启动本钱的计算。〔2〕我们提供了一种求解多变量，多约束整数非线性规划的组合优化问题的思路，此方法思路清晰明了，构思新颖，方便易行，极具参考价值。〔3〕我们所建立的非线性规划模型，很好的解决了发电机厂关于每天发电机组的启停计划，通过模型的求解得出了比拟合理的结果，得出了每天发电机组本钱的最小值，通过验证，我们所建立的模型具有很好的实用价值，可以推广到其他电厂的应用。模型的缺点：根据题目的要求，以与其他条件的限制和约束，对模型进展了一些合理的假设，虽然求出的结果在所要求的误差围之，但是假设还是会对所求结果有一定的影响。对于这些方面还需要做进一步的探讨和改良。7模型的应用和推广模型的应用：通过所建立的模型以与求解结果，我们可以得到，合理的机组优化将会大幅度的减少机组的发电本钱，较低企业消耗，有效的增加企业的经济效益，将模型应用与实际的生活具有非常重要的现实意义。从得到的优化结果中可以看出，每个机型的利用率是不同的，机型1的利用率较低，而机型3的利用率较高，根据这样的结果可以指导电力企业合理的采购，各种机型，合理的配置机型数量，使得所有机组都能够得到应用，减少浪费。根据每个时间段的用电需求，可以得到在在0－6时间段是用电低谷，即晚上的时候用电量较低，而在12－14时间段是用电顶峰期，这也使符合时间生活的规律，夏季中午温度较高，需要各种电器设备来降温消暑，所以用电增加，针对用电顶峰期，我们可以给出建立，例如：尽量的使用风扇，开窗通风，适当的提高空调温度，集中制冷等方法节约用电，从而节约能源，减少用电压力。模型推广：〔1〕此模型具有一定的适用性，当每个时间段的用电需求发生变化时，只需要对模型的需求约束进展相应的改变就可以重新的求解优化结果。〔2〕同样当提供的机型数量发生变化时，此模型也可以通过改变其，机型数量约束来重新求解，以达到要求。〔3〕此模型不仅仅限于电力生产问题，也可以应用于其他的行业，即在有限的资源条件下，合理的分配资源得到相应的要求，同时使得消耗最小。综上，此模型具有很好的应用价值，可以作进一步的推广和研究。参考文献[1]宋来忠,王志明.数学建模与实验.：科学，2005.[2]王正东.数学软件与数学实验.：科学，2010.[3]静，但琦，数学建模与数学实验.：高等教育，2008.[4]楼顺天，假设玉，俊霞.MATLAB7.x程序设计语言.：电子科技大学，2008.[5]国立.现货市场功率分配问题模型与其求解算法研究[J],[6.：械工业，2007，43—89.[7]吴礼斌，柏年.数学实验与建模.：国防工业，2007.[8]金星.优化建模与LINDO/LINGO软件.：清华大学，2005.附录!目标函数;min=(6*x11+3*x21+3*x31+2*x41+4*x51+4*x61+2*x71)*2250+(6*x11*(a11-750)+3*x21*(a21-750)+3*x31*(a31-750)+2*x41*(a41-750)+4*x51*(a51-750)+4*x61*(a61-750)+2*x71*(a71-750))*2.7+((abs(x11-x01)+(x11-x01))+(abs(x21-x11)+(x21-x11))+(abs(x31-x21)+(x31-x21))+(abs(x41-x31)+(x41-x31))+(abs(x51-x41)+(x51-x41))+(abs(x61-x51)+(x61-x51))+(abs(x71-x61)+(x71-x61)))*5000/2+(6*x12+3*x22+3*x32+2*x42+4*x52+4*x62+2*x72)*1800+(6*x12*(a12-1000)+3*x22*(a22-1000)+3*x32*(a32-1000)+2*x42*(a42-1000)+4*x52*(a52-1000)+4*x62*(a62-1000)+2*x72*(a72-1000))*2.2+((abs(x12-x02)+(x12-x02))+(abs(x22-x12)+(x22-x12))+(abs(x32-x22)+(x32-x22))+(abs(x42-x32)+(x42-x32))+(abs(x52-x42)+(x52-x42))+(abs(x62-x52)+(x62-x52))+(abs(x72-x62)+(x72-x62)))*1600/2+(6*x13+3*x23+3*x33+2*x43+4*x53+4*x63+2*x73)*3750+(6*x13*(a13-1200)+3*x23*(a23-1200)+3*x33*(a33-1200)+2*x43*(a43-1200)+4*x53*(a53-1200)+4*x63*(a63-1200)+2*x73*(a73-1200))*1.8+((abs(x13-x03)+(x13-x03))+(abs(x23-x13)+(x23-x13))+(abs(x33-x23)+(x33-x23))+(abs(x43-x33)+(x43-x33))+(abs(x53-x43)+(x53-x43))+(abs(x63-x53)+(x63-x53))+(abs(x73-x63)+(x73-x63)))*2400/2+(6*x14+3*x24+3*x34+2*x44+4*x54+4*x64+2*x74)*4800+(6*x14*(a14-1800)+3*x24*(a24-1800)+3*x34*(a34-1800)+2*x44*(a44-1800)+4*x54*(a54-1800)+4*x64*(a64-1800)+2*x74*(a74-1800))*3.8+((abs(x14-x04)+(x14-x04))+(abs(x24-x14)+(x24-x14))+(abs(x34-x24)+(x34-x24))+(abs(x44-x34)+(x44-x34))+(abs(x54-x44)+(x54-x44))+(abs(x64-x54)+(x64-x54))+(abs(x74-x64)+(x74-x64)))*1200/2;!约束条件;x11*a11+x12*a12+x13*a13+x14*a14>=12000;x21*a21+x22*a22+x23*a23+x24*a24>=32000;x31*a31+x32*a32+x33*a33+x34*a34>=25000;x41*a41+x42*a42+x43*a43+x44*a44>=36000;x51*a51+x52*a52+x53*a53+x54*a54>=25000;x61*a61+x62*a62+x63*a63+x64*a64>=30000;x71*a71+x72*a72+x73*a73+x74*a74>=18000;a11>=750;a11<=1750;a21>=750;a21<=1750;a31>=750;a31<=1750;a41>=750;a41<=1750;a51>=750;a51<=1750;a61>=750;a61<=1750;a71>=750;a71<=1750;a12>=1000;a12<=1500;a22>=1000;a22<=1500;a32>=1000;a32<=1500;a42>=1000;a42<=1500;a52>=1000;a52<=1500;a62>=1000;a62<=1500;a72>=1000;a72<=1500;a13>=1200;a13<=2000;a23>=1200;a23<=2000;a33>=1200;a33<=2000;a43>=1200;a43<=2000;a53>=1200;a53<=2000;a63>=1200;a63<=2000;a73>=1200;a73<=2000;a14>=1800;a14<=3500;a24>=1800;a24<=3500;a34>=1800;a34<=3500;a44>=1800;a44<=3500;a54>=1800;a54<=3500;a64>=1800;a64<=3500;a74>=1800;a74<=3500;x11<=10;x21<=10;x31<=10;x41<=10;x51<=10;x61<=10;x71<=10;x12<=4;x22<=4;x32<=4;x42<=4;x52<=4;x62<=4;x72<=4;x13<=8;x23<=8;x33<=8;x43<=8;x53<=8;x63<=8;x73<=8;x14<=3;x24<=3;x34<=3;x44<=3;x54<=3;x64<=3;x74<=3;gin(x11);gin(x21);gin(x31);gin(x41);gin(x51);gin(x61);gin(x71);gin(x12);gin(x22);gin(x32);gin(x42);gin(x52);gin(x62);gin(x72);gin(x13);gin(x23);gin(x33);gin(x43);gin(x53);gin(x63);gin(x73);gin(x14);gin(x24);gin(x34);gin(x44);gin(x54);gin(x64);gin(x74);x01=0;x02=0;x03=0;x04=0;Localoptimalsolutionfound.Objectivevalue:1463625.Objectivebound:1463625.Extendedsolversteps:33Totalsolveriterations:12410VariableValueRowSlackorSurplus11463625.!目标函数;min=(6*x11+3*x21+3*x31+2*x41+4*x51+4*x61+2*x71)*2250+(6*x11*(a11-750)+3*x21*(a21-750)+3*x31*(a31-750)+2*x41*(a41-750)+4*x51*(a51-750)+4*x61*(a61-750)+2*x71*(a71-750))*2.7+((abs(x11-x01)+(x11-x01))+(abs(x21-x11)+(x21-x11))+(abs(x31-x21)+(x31-x21))+(abs(x41-x31)+(x41-x31))+(abs(x51-x41)+(x51-x41))+(abs(x61-x51)+(x61-x51))+(abs(x71-x61)+(x71-x61)))*5000/2+(6*x12+3*x22+3*x32+2*x42+4*x52+4*x62+2*x72)*1800+(6*x12*(a12-1000)+3*x22*(a22-1000)+3*x32*(a32-1000)+2*x42*(a42-1000)+4*x52*(a52-1000)+4*x62*(a62-1000)+2*x72*(a72-1000))*2.2+((abs(x12-x02)+(x12-x02))+(abs(x22-x12)+(x22-x12))+(abs(x32-x22)+(x32-x22))+(abs(x42-x32)+(x42-x32))+(abs(x52-x42)+(x52-x42))+(abs(x62-x52)+(x62-x52))+(abs(x72-x62)+(x72-x62)))*1600/2+(6*x13+3*x23+3*x33+2*x43+4*x53+4*x63+2*x73)*3750+(6*x13*(a13-1200)+3*x23*(a23-1200)+3*x33*(a33-1200)+2*x43*(a43-1200)+4*x53*(a53-1200)+4*x63*(a63-1200)+2*x73*(a73-1200))*1.8+((abs(x13-x03)+(x13-x03))+(abs(x23-x13)+(x23-x13))+(abs(x33-x23)+(x33-x23))+(abs(x43-x33)+(x43-x33))+(abs(x53-x43)+(x53-x43))+(abs(x63-x53)+(x63-x53))+(abs(x73-x63)+(x73-x63)))*2400/2+(6*x14+3*x24+3*x34+2*x44+4*x54+4*x64+2*x74)*4800+(6*x14*(a14-1800)+3*x24*(a24-1800)+3*x34*(a34-1800)+2*x44*(a44-1800)+4*x54*(a54-1800)+4*x64*(a64-1800)+2*x74*(a74-1800))*3.8+((abs(x14-x04)+(x14-x04))+(abs(x24-x14)+(x24-x14))+(abs(x34-x24)+(x34-x24))+(abs(x44-x34)+(x44-x34))+(abs(x54-x44)+(x54-x44))+(abs(x64-x54)+(x64-x54))+(abs(x74-x64)+(x74-x64)))*1200/2;!约束条件;x11*a11+x12*a12+x13*a13+x14*a14>=15000;x21*a21+x22*a22+x23*a23+x24*a24>=40000;x31*a31+x32*a32+x33*a33+x34*a34>=31250;x41*a41+x42*a42+x43*a43+x44*a44>=45000;x51*a51+x52*a52+x53*a53+x54*a54>=31250;x61*a61+x62*a62+x63*a63+x64*a64>=37500;x71*a71+x72*a72+x73*a73+x74*a74>=22500;a11>=750;a11<=1750;a21>=750;a21<=1750;a31>=750;a31<=1750;a41>=750;a41<=1750;a51>=750;a51<=1750;a61>=750;a61<=1750;a71>=750;a71<=1750;a12>=1000;a12<=1500;a22>=1000;a22<=1500;a32>=1000;a32<=1500;a42>=1000;a42<=1500;a52>=1000;a52<=1500;a62>=1000;a62<=1500;a72>=1000;a72<=1500;a13>=1200;a13<=2000;