《电力生产问题_最优化方法论文》.docxVIP

 PAGE \* MERGEFORMAT 10 电力生产问题 问题分析 为满足每日七个时段各不相同的电力需求，又达到总成本最小，需要在每个时间段，调整四种发电机的开机台数和输出功率，实现已启动的发电机的启动成本、固定成本和边际成本的和最小。为了实现成本最优化目标，每一个时间段，不同型号的发电机开启的数量可能是不同的。因此，一个新时段需要启动的某种型号的发电机数量，与上一时段此型号发电机已经运行的数量是相关的，这会进一步影响此型号发电机的启动成本值。这需要对七个时段，发电机的启动数量进行逐步计算。 我们基于以下假设条件建立最优化方法目标函数及约束条件： （1）设每日零时为计算起点，24时为计算终点。起点时刻需要开启的各种型号的发电机数量，为此时间段需要运行的发电机数量减去前一天终点时刻相应型号发电机运行的数量。 （2）设每个单独时间段内，相同型号的发电机工作的输出功率相同，且维持不变； （3）假设所有发电机在整个运行过程中都能正常工作； （4）假设发电机的输出电量大于等于用户需求量。 变量定义 具体的变量定义如下表所示： 符号符号说明Cs全天启动成本和（单位：欧元）CF全天固定成本和（单位：欧元）CB全天边际成本和（单位：欧元）Ti第i个时间段的时间长度（单位：小时）Wminj型号j发电机最小输出功率（单位：兆瓦）Wmaxj型号j发电机最大输出功率（单位：兆瓦）Wij第i个时间段型号j发电机的实际输出功率（单位：兆瓦）Nij第i个时间段，型号j发电机的运行数量（单位：台）Ri第i个时间段的用电需求（单位：兆瓦）Sj型号j发电机的启动成本（单位：欧元/小时）Fj型号j发电机的固定成本（单位：欧元/小时）Bj型号j发电机的边际成本（单位：欧元/（小时*兆瓦））备注：i=1, …,7; j=1, 2, 3, 4。 模型建立 1.目标函数 题目要求的每天最小总成本为已启动的发电机的启动成本、固定成本和边际成本和的最小值，即： min (Cs+CF+CB) 其中，启动成本可表示为： Cs=j=14Sj×(N1j-N7j+N1j-N7j)2+i=2, j=1i=7, j=4Sj×(Nij-N(i-1)j+Nij-N(i-1)j)2 其中 ，S=5000, 1600, 2400, 1200, Nij为整数。 固定成本和可表示为： CF=i=1, j=1i=7, j=4Fj×Ti×Nij 其中 ，F=(2250, 1800, 3750, 4800)，T=(6, 3, 3, 2, 4, 4, 2)。 边际成本和可表示为： CB=i=1, j=1i=7, j=4(Wij-Wminj)×Ti×Nij×Bj 其中 ，Wmin=(750, 1000, 1200, 1800),B=(2.7, 2.2, 1.8, 3.8)，T=(6,3,3,2,4,4,2)。 2.约束条件 目标函数条件以下约束条件： S.T.750≤Wi1Wi2Wi3Wi4≤35000≤Ni1≤100≤Ni2≤40≤Ni3≤80≤Ni4≤3 j=1j=4Wij×Nij-Ri≥00.8× j=1j=4(Wmaxj×Nij)- j=1j=4(Wij×Nij)≥0i=1, …,7 其中,R=(12000, 32000, 25000, 36000, 25000, 30000, 18000),Wmax=(1750, 1500, 2000, 3500)。 3.模型求解 利用matlab对问题求解，具体结果如下表所示（程序见附录1）： 型 号 数 量 功 率 时 段 型号一型号二型号三型号四0-61台/750兆瓦4台/1000兆瓦1台/2000兆瓦3台/2100兆瓦6-97台/750兆瓦4台/1000兆瓦8台/2000兆瓦3台/2800兆瓦9-127台/750兆瓦4台/1000兆瓦8台/2000兆瓦3台/1800兆瓦12-148台/750兆瓦4台/1000兆瓦8台/2000兆瓦3台/3333兆瓦14-185台/750兆瓦4台/1000兆瓦8台/2000兆瓦2台/1900兆瓦18-225台/750兆瓦4台/1000兆瓦8台/1950兆瓦3台/2300兆瓦22-243台/750兆瓦4台/1000兆瓦7台/1800兆瓦1台/1800兆瓦总成本：1317500欧元 由于matlab的fmincon函数不支持非线性整数规划问题的求解，所以用LINGO对问题再次求解，具体结果如下表所示（程序见附录2）： 型 号 数 量 功 率 时 段 型号一型