发明名称 微电网多目标转化单目标的方法
摘要 一种微电网多目标转化单目标的方法,属于微电网优化运行领域。本发明的目的利用模糊处理的方法,将多目标问题转化为单目标,并提出一种改进蚁群算法对所建模型进行求解的微电网多目标转化单目标的方法。本发明的步骤是:得到风机、光伏的共同出力的期望值,在考虑机会约束规划的同时建立优化模型,选取待优化变量、优化目标及确定约束条件,分别选取适合两个目标函数的隶属度函数,并确定隶属度函数。本发明以包含光伏、风机、微型燃气轮机、蓄电池的微网为研究对象,提出了包含经济、环保的多目标经济运行模型。
申请公布号 CN106602593A 申请公布日期 2017.04.26
申请号 CN201611007282.3 申请日期 2016.11.16
申请人 东北电力大学 发明人 李扬;翟晓娟;李国庆;张明江;王振浩;金鹏;陈继开
分类号 H02J3/38(2006.01)I;H02J3/46(2006.01)I;G06F17/50(2006.01)I;G06N3/00(2006.01)I 主分类号 H02J3/38(2006.01)I
代理机构 吉林长春新纪元专利代理有限责任公司 22100 代理人 白冬冬
主权项 一种微电网多目标转化单目标的方法,其特征在于:①建立风机和光伏的不确定性出力模型,利用序列运算理论得到风机、光伏的共同出力的期望值,具体过程为:a、根据对实测数据的分析表明,风速近似服从威布尔分布,其分布函数为:<img file="590083dest_path_image001.GIF" wi="251" he="35" />(1)其概率密度函数为:<img file="342007dest_path_image002.GIF" wi="409" he="40" />(2)式(1)、(2)中,<img file="725715dest_path_image003.GIF" wi="16" he="16" />代表实际风速;<img file="122061dest_path_image004.GIF" wi="19" he="26" />代表形状参数,一般取<img file="683755dest_path_image005.GIF" wi="43" he="22" />;<img file="798341dest_path_image006.GIF" wi="14" he="17" />代表尺度参数,反映某个时段的平均风速,有<img file="454582dest_path_image007.GIF" wi="148" he="32" />;风力机组的输出功率<img file="361227dest_path_image008.GIF" wi="30" he="23" />和实际风速<img file="202144dest_path_image003.GIF" wi="16" he="16" />的函数关系如下:<img file="679393dest_path_image009.GIF" wi="436" he="210" />(3)式(3)中,<img file="263958dest_path_image010.GIF" wi="28" he="30" />表示额定风速;<img file="541223dest_path_image011.GIF" wi="33" he="29" />表示切入风速;<img file="553041dest_path_image012.GIF" wi="44" he="30" />表示切出风速;<img file="517586dest_path_image013.GIF" wi="28" he="36" />表示风机的额定输出功率;b、光伏电池组件的功率输出表达式如下:<img file="30476dest_path_image014.GIF" wi="424" he="81" />(4)式(4)中,<img file="255921dest_path_image015.GIF" wi="41" he="38" />为光伏电池的实际输出功率;<img file="579586dest_path_image016.GIF" wi="52" he="34" />为标准测试条件下的电池输出功率;<img file="156061dest_path_image017.GIF" wi="61" he="37" />为实际辐照强度;<img file="708527dest_path_image018.GIF" wi="61" he="33" />为标准测试条件下的辐照强度;<img file="929424dest_path_image019.GIF" wi="30" he="34" />为功率温度系数;<img file="548624dest_path_image020.GIF" wi="27" he="36" />为电池板工作温度;<img file="737029dest_path_image021.GIF" wi="27" he="34" />为参考温度; 光伏电池出力的概率密度函数为:<img file="467087dest_path_image022.GIF" wi="550" he="107" />(5)式(5)中,<img file="276911dest_path_image023.GIF" wi="20" he="24" />、<img file="955761dest_path_image024.GIF" wi="20" he="29" />为形状参数,其关系为<img file="506828dest_path_image025.GIF" wi="191" he="39" />,<img file="650365dest_path_image026.GIF" wi="18" he="20" />为光强平均值;<img file="95121dest_path_image027.GIF" wi="68" he="40" />为光伏电池的最大输出功率;c、令风机出力<img file="56124dest_path_image008.GIF" wi="30" he="23" />的概率性序列为<img file="704274dest_path_image028.GIF" wi="63" he="39" />,序列长度为<img file="776135dest_path_image029.GIF" wi="35" he="31" />,光伏出力<img file="576863dest_path_image015.GIF" wi="41" he="38" />的概率性序列为<img file="584134dest_path_image030.GIF" wi="51" he="30" />,序列长度为<img file="844214dest_path_image031.GIF" wi="34" he="34" />,令风机、光伏的共同出力的概率性序列为<img file="844400dest_path_image032.GIF" wi="52" he="32" />,序列长度为<img file="873536dest_path_image033.GIF" wi="34" he="34" />;共同出力的计算公式如下:<img file="51707dest_path_image034.GIF" wi="200" he="93" />(6)<img file="428112dest_path_image035.GIF" wi="359" he="76" />(7)公式(6)中,<img file="107355dest_path_image036.GIF" wi="31" he="36" />为<img file="866363dest_path_image037.GIF" wi="18" he="24" />时段风机和光伏共同出力的期望值;式(7)完整的表述如下:<img file="730283dest_path_image038.GIF" wi="599" he="42" />(8);②在考虑机会约束规划的同时建立优化模型,选取待优化变量、优化目标及确定约束条件:a、选取待优化变量,连续变量为各台微燃机在<img file="699376dest_path_image037.GIF" wi="18" he="24" />时刻的有功出力<img file="792097dest_path_image039.GIF" wi="29" he="38" />,离散变量主要包括机组状态变量<img file="795825dest_path_image040.GIF" wi="27" he="31" />和开机变量;b、选取优化目标,优化目标函数表达式为:<img file="800952dest_path_image041.GIF" wi="560" he="83" />(9)<img file="398287dest_path_image042.GIF" wi="430" he="91" />(10)式(9)中,<img file="153753dest_path_image043.GIF" wi="26" he="30" />为微电网运行发电成本;<img file="402201dest_path_image039.GIF" wi="29" he="38" />为机组<img file="217710dest_path_image044.GIF" wi="14" he="24" />在<img file="36762dest_path_image037.GIF" wi="18" he="24" />时刻的出力;<img file="219088dest_path_image045.GIF" wi="78" he="40" />为机组<img file="931829dest_path_image044.GIF" wi="14" he="24" />在<img file="59185dest_path_image037.GIF" wi="18" he="24" />时刻表示机组的燃料消耗成本;<img file="490166dest_path_image046.GIF" wi="26" he="40" />为0‑1变量,机组<img file="712069dest_path_image044.GIF" wi="14" he="24" />在<img file="420262dest_path_image037.GIF" wi="18" he="24" />时刻的运行状态,该值是0表示停机,该值是1表示运行;<img file="577574dest_path_image047.GIF" wi="22" he="31" />为机组<img file="387529dest_path_image044.GIF" wi="14" he="24" />的启动成本;<img file="22910dest_path_image048.GIF" wi="19" he="31" />为机组<img file="851189dest_path_image044.GIF" wi="14" he="24" />的备用容量成本;<img file="304036dest_path_image049.GIF" wi="27" he="40" />为机组<img file="709609dest_path_image044.GIF" wi="14" he="24" />在第<img file="289626dest_path_image037.GIF" wi="18" he="24" />时刻的旋转备用;式(10)中,<img file="714355dest_path_image050.GIF" wi="28" he="32" />为污染气体的排放总量;<img file="213469dest_path_image051.GIF" wi="28" he="30" />、<img file="716126dest_path_image052.GIF" wi="30" he="43" />、<img file="693309dest_path_image053.GIF" wi="30" he="32" />为微燃机<img file="745448dest_path_image044.GIF" wi="14" he="24" />的污染气体排放系数;c、确定约束条件,约束条件包括等式约束和不等式约束两部分:等式约束写为:<img file="415463dest_path_image054.GIF" wi="419" he="147" />(11)式(11)中,<img file="670995dest_path_image055.GIF" wi="31" he="45" />为<img file="77968dest_path_image037.GIF" wi="18" he="24" />时刻的卸荷负载出力;<img file="594400dest_path_image056.GIF" wi="32" he="39" />为储能在<img file="576263dest_path_image037.GIF" wi="18" he="24" />时刻的充电功率;<img file="178145dest_path_image057.GIF" wi="36" he="40" />为储能在<img file="152923dest_path_image037.GIF" wi="18" he="24" />时刻的放电功率;<img file="399228dest_path_image058.GIF" wi="38" he="39" />为<img file="411046dest_path_image037.GIF" wi="18" he="24" />时刻的负荷预测值;<img file="123394dest_path_image059.GIF" wi="63" he="35" />、<img file="511650dest_path_image060.GIF" wi="66" he="38" />为蓄电池运行一个周期的始末容量;不等式约束条件:<img file="878040dest_path_image061.GIF" wi="510" he="359" />(12)式(12)中,<img file="185394dest_path_image062.GIF" wi="62" he="41" />、<img file="27448dest_path_image063.GIF" wi="62" he="40" />分别为机组出力的上下限;<img file="829182dest_path_image064.GIF" wi="54" he="40" />、<img file="800811dest_path_image065.GIF" wi="56" he="38" />分别为蓄电池所允许的剩余电量的最小最大值;<img file="420011dest_path_image066.GIF" wi="33" he="41" />为备用容量;<img file="93569dest_path_image067.GIF" wi="33" he="38" />、<img file="823628dest_path_image068.GIF" wi="43" he="43" />分别为<img file="148299dest_path_image037.GIF" wi="18" he="24" />时段风机和光伏的出力;<img file="79346dest_path_image069.GIF" wi="21" he="30" />为给定的置信度;<img file="364834dest_path_image070.GIF" wi="46" he="36" />、<img file="527611dest_path_image071.GIF" wi="73" he="40" />分别为微燃机增加和降低有功功率的极限值;③分别选取适合两个目标函数的隶属度函数,并确定隶属度函数:两个目标函数选取的隶属度函数分别为:<img file="723101dest_path_image072.GIF" wi="387" he="155" />(13)<img file="418524dest_path_image073.GIF" wi="401" he="153" />(14)式(13)、(14)中,<img file="315942dest_path_image074.GIF" wi="23" he="30" />为以微电网发电成本最小为单目标优化的目标值;<img file="653382dest_path_image075.GIF" wi="25" he="26" />为以污染物排放总量最小为单目标优化的目标值;<img file="437799dest_path_image076.GIF" wi="21" he="31" />为可以接受的成本的增加值;<img file="461381dest_path_image077.GIF" wi="24" he="28" />为可接受的污染物排放的伸缩值;E指的是污染气体的排放量;F指的是发电成本;<img file="455881dest_path_image078.GIF" wi="55" he="30" />有关污染气体的隶属度函数;<img file="472379dest_path_image079.GIF" wi="54" he="30" />有关发电成本的隶属度函数;确定隶属度函数的过程为:a、 将原始数据输入,利用混合整数规划算法,求解以微电网发电成本最小为目标函数的单目标模型,得到了一天中各时段的最小发电成本<img file="501515dest_path_image074.GIF" wi="23" he="30" />、最大发电成本<img file="928954dest_path_image080.GIF" wi="43" he="31" />和污染气体总排放量<img file="551696dest_path_image081.GIF" wi="26" he="31" />,以及各台微燃机的出力情况;b、将原始数据输入,用上述同样的方法,求解以微网系统运行污染物排放总量为目标函数的单目标模型,得到了一天中各时段污染气体最小总排放量<img file="965360dest_path_image075.GIF" wi="25" he="26" />、最大总排放量<img file="737751dest_path_image082.GIF" wi="44" he="34" />和各时段的发电成本<img file="945878dest_path_image083.GIF" wi="22" he="28" />,以及各台微燃机的出力情况;c、在a、b两个步骤的基础上,将各个单目标进行伸缩,确定<img file="790337dest_path_image076.GIF" wi="21" he="31" />和<img file="397905dest_path_image077.GIF" wi="24" he="28" />,同时<img file="136054dest_path_image076.GIF" wi="21" he="31" />和<img file="921607dest_path_image077.GIF" wi="24" he="28" />的确定方法:<img file="377997dest_path_image084.GIF" wi="186" he="38" />,<img file="759562dest_path_image085.GIF" wi="183" he="38" />;d、将<img file="758742dest_path_image076.GIF" wi="21" he="31" />和<img file="839830dest_path_image077.GIF" wi="24" he="28" />代入式(13)和式(14)中,得到隶属度函数,从而将多目标转化为单目标问题;设<img file="642570dest_path_image086.GIF" wi="23" he="22" />为所有隶属度函数中的最小值,最终多目标问题就可以转化为在满足所有约束条件的前提下,满意度指标<img file="201727dest_path_image086.GIF" wi="23" he="22" />最大的单目标优化问题;转化后的单目标模型如下:<img file="789834dest_path_image087.GIF" wi="61" he="23" />(15)<img file="682571dest_path_image088.GIF" wi="170" he="33" />(16)<img file="113553dest_path_image089.GIF" wi="161" he="31" />(17)<img file="351767dest_path_image090.GIF" wi="83" he="24" />(18)。
地址 132012 吉林省吉林市船营区长春路169号