发明名称 考虑蓄冰空调并具有电冷耦合特征的微电网优化调度方法
摘要 本发明公开了一种考虑蓄冰空调并具有电冷耦合特征的微电网优化调度方法,步骤包括:建立微电网供能设备动态运行模型,包括:蓄冰空调机组,冷热电联供(CCHP)机组;建立耦合型电冷需求微电网优化调度模型和建立基于改进粒子群技术的微电网运行与调度计算方法。通过上述方式,本发明考虑蓄冰空调并面向耦合型多能源需求的微电网优化调度方法,可以实现可再生能源和清洁能源的高效利用;其次,它可以以最低的成本满足用户的多种能源需求;第三,它可以提高大电网的供电可靠性。
申请公布号 CN105631557A 申请公布日期 2016.06.01
申请号 CN201610151800.2 申请日期 2016.03.17
申请人 北京睿新科技有限公司 发明人 张建华
分类号 G06Q10/04(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 苏州广正知识产权代理有限公司 32234 代理人 徐萍
主权项 一种考虑蓄冰空调并具有电冷耦合特征的微电网优化调度方法,其特征在于,步骤包括:(1)建立精细准确的微电网供能设备动态运行模型:(1.1)建立CCHP的运行分析模型:(1.1.1)设置电功率输出效率模型为:<img file="491017dest_path_image001.GIF" wi="200" he="26" />,其中,E<sub>CHP</sub>(t)为CCHP的电功率输出效率,P<sub>CHP</sub>(t)为CCHP在t时刻的电功率输出;(1.1.2)在以冷定电的工作方式下,制冷量‑发电功率输出耦合关系为:<img file="991269dest_path_image002.GIF" wi="200" he="27" />,其中,Q<sub>CHP</sub>(t)为CCHP 在时刻t供冷输出;(1.1.3)根据特定工况进行修正建模,与特定工况有关的参数包括:CCHP机组的工作环境温度T、海拔E、入口压力损失<img file="552831dest_path_image003.GIF" wi="120" he="22" />和背压<img file="572871dest_path_image004.GIF" wi="88" he="23" />;特定工况下CCHP机组的余烟气制冷量与发电量的关系为<img file="748638dest_path_image005.GIF" wi="192" he="27" />,<img file="549234dest_path_image006.GIF" wi="499" he="134" />(1)其中,<img file="586592dest_path_image007.GIF" wi="112" he="32" />是在特定工况下CCHP机组部分负荷运行时的烟气排放温度, E_T<sub>max</sub>是在特定工况下CCHP机组满负荷运行时的烟气排放温度;P<sub>max</sub>、<img file="975985dest_path_image008.GIF" wi="80" he="28" />分别为特定工况下CCHP机组满负荷运行时的输出电功率和部分负荷运行时的输出电功率;<img file="807806dest_path_image009.GIF" wi="119" he="32" />是在特定工况下CCHP机组部分负荷运行时的烟气排放速度,E_V<sub>max</sub>为特定工况下CCHP机组满负荷运行时的烟气排放速度;<img file="282649dest_path_image010.GIF" wi="435" he="58" />,<img file="310648dest_path_image011.GIF" wi="381" he="111" />,根据以上关系可求得<img file="305280dest_path_image012.GIF" wi="208" he="28" />,即  (2)其中,<img file="822849dest_path_image013.GIF" wi="320" he="29" />,由三个子函数合并而成,其中,P<sub>A‑max</sub>表示不同海拔E和环境温度T的工况下满负荷运行时输出电功率,<img file="264283dest_path_image014.GIF" wi="102" he="27" />表示 ISO工况下入口压力损失<img file="892710dest_path_image015.GIF" wi="136" he="27" />对满负荷运行时的输出电功率,<img file="7428dest_path_image016.GIF" wi="118" he="32" />表示ISO工况下的背压<img file="695898dest_path_image017.GIF" wi="88" he="22" />对满负荷运行时的输出电功率,如式(3)‑(5)所示:<img file="630487dest_path_image018.GIF" wi="440" he="247" />(3)<img file="62606dest_path_image019.GIF" wi="328" he="32" />(4)<img file="31830dest_path_image020.GIF" wi="336" he="32" />(5)CCHP机组消耗的燃气量公式为:<img file="891201dest_path_image021.GIF" wi="384" he="34" />(6)其中,11.6为燃气热值,特定工况下CCHP机组出力<img file="313087dest_path_image022.GIF" wi="80" he="26" />所对应的效率<img file="283317dest_path_image023.GIF" wi="104" he="27" />为:<img file="107047dest_path_image024.GIF" wi="424" he="99" />(7)其中,<img file="137320dest_path_image025.GIF" wi="96" he="38" />是CCHP机组在ISO工况下满负荷运行时的输出电功率,一般为常数;<img file="46501dest_path_image026.GIF" wi="64" he="34" />是在特定工况下CCHP机组满负荷运行时的效率,<img file="820422dest_path_image027.GIF" wi="417" he="37" />,同上,也是由三个子函数组成,其中,<img file="498660dest_path_image028.GIF" wi="88" he="36" />表示不同海拔E和环境温度T的工况下满负荷运行时效率,<img file="699834dest_path_image029.GIF" wi="120" he="32" />表示ISO工况下入口压力损失<img file="90452dest_path_image030.GIF" wi="136" he="27" />下满负荷运行时效率,<img file="668064dest_path_image031.GIF" wi="104" he="27" />表示ISO工况下背压<img file="200807dest_path_image032.GIF" wi="112" he="29" />下满负荷运行时效率,如式(8)‑(10)所示:<img file="572883dest_path_image033.GIF" wi="444" he="46" />(8)<img file="456656dest_path_image034.GIF" wi="408" he="28" />(9)<img file="572380dest_path_image035.GIF" wi="408" he="32" />(10)在CCHP机组的单机和多机联合运行时,<img file="225209dest_path_image036.GIF" wi="88" he="32" />取不同的拟合函数,当单台CCHP机组工作时,<img file="768186dest_path_image038.GIF" wi="417" he="26" />,当多台CCHP机组组成燃气机组工作时,<img file="201572dest_path_image039.GIF" wi="440" he="24" />;(1.2)建立蓄冰空调机组模型:蓄冰空调机组由制冷机、蓄冰槽组成,所述蓄冰空调机组模型中包括四种运行工况:1)直供,单纯消耗电能供冷的蓄冰空调直供工况,2)蓄冷,在冷负荷低谷期耗电制冰存储的蓄冰空调蓄冷工况,3) 释冰,在用冷高峰期融冰供冷的蓄冰空调融冰释冷工况,4)联供,释冰与直供并行的蓄冰空调联供工况;所述模型由下述(2.2)中(c.1)‑(c.3)描述确定;(2)建立耦合型电冷需求微电网优化调度模型;(2.1)设立优化目标函数,计算微电网的运行总成本最小,即<img file="606140dest_path_image040.GIF" wi="442" he="78" />(11)其中,公式中第一项<img file="362744dest_path_image041.GIF" wi="107" he="40" />为微电网向大电网购电的成本,第二项<img file="92933dest_path_image042.GIF" wi="89" he="43" />为供给CCHP的燃气购买成本,<img file="748037dest_path_image043.GIF" wi="112" he="40" />为微电网与大电网连接点PCC交换功率的电价,<img file="205563dest_path_image044.GIF" wi="104" he="40" />为燃气的价格,<img file="827125dest_path_image045.GIF" wi="96" he="32" />为微网与大电网在PCC点的交换功率,F(t)为燃气的消耗量;待求变量包括: 直供状态<img file="711904dest_path_image046.GIF" wi="64" he="33" />、蓄冰状态<img file="791987dest_path_image047.GIF" wi="72" he="32" />、释冰状态<img file="53204dest_path_image048.GIF" wi="56" he="26" />、直供冷功率<img file="597449dest_path_image049.GIF" wi="72" he="32" />、释冰冷功率<img file="403862dest_path_image050.GIF" wi="64" he="32" />、CCHP机组产生的电功率<img file="220508dest_path_image051.GIF" wi="88" he="32" />、微网与大电网的交换功率<img file="770570dest_path_image052.GIF" wi="88" he="32" />,冷功率<img file="356272dest_path_image053.GIF" wi="80" he="26" />与消耗的燃气量F(t)作为中间变量,可通过<img file="599165dest_path_image054.GIF" wi="80" he="32" />获得,其中,<img file="903108dest_path_image055.GIF" wi="152" he="32" />;(2.2)设置在本模型中的约束条件:(a)设置电平衡的约束条件为:<img file="991281dest_path_image056.GIF" wi="364" he="71" />(12)(b)设置冷需求平衡的约束条件为:<img file="697068dest_path_image057.GIF" wi="437" he="35" />(13)(c)设置蓄冰空调约束:蓄冰空调的制冷机具有直供、蓄冷两种工况,蓄冰槽具有释冰工况,夜间谷时段根据第二天的冷负荷需求,制冷机提供蓄冷;峰时段根据实际冷负荷需求,制冷机和蓄冰槽提供释冰、直供、联供3种方式进行供冷,所述联供包括直供和释冰,定义直供的状态为<img file="845284dest_path_image058.GIF" wi="64" he="32" />,蓄冰的状态为<img file="636523dest_path_image059.GIF" wi="72" he="40" />,释冰的状态为<img file="545964dest_path_image060.GIF" wi="56" he="32" />,且满足:<img file="106259dest_path_image061.GIF" wi="64" he="32" />∈(0,1),<img file="877906dest_path_image062.GIF" wi="64" he="32" />∈(0,1),<img file="641594dest_path_image063.GIF" wi="64" he="32" />∈(0,1)          (14)(c.1)直供工况:谷时段双工况机组不能工作于直供状态,谷时段制冷由基载完成,即满足约束:<img file="586416dest_path_image064.GIF" wi="120" he="40" />, t∈(0,…,8,23),            (15)峰时段为了保证双工况机组较高的工作效率,<img file="751949dest_path_image065.GIF" wi="72" he="33" />需满足约束:<img file="756814dest_path_image066.GIF" wi="416" he="32" />(16)当工作于上述公式范围内时,COP与<img file="7798dest_path_image067.GIF" wi="80" he="40" />呈线性关系,则输出制冷功率<img file="490732dest_path_image049.GIF" wi="72" he="32" />与消耗电功率<img file="776351dest_path_image068.GIF" wi="64" he="32" />之间满足如下关系:<img file="952117dest_path_image069.GIF" wi="368" he="89" />(17)(c.2)蓄冷工况:蓄冷量<img file="690397dest_path_image070.GIF" wi="72" he="32" />为0,或者为双工况机的额定制冷量<img file="977022dest_path_image071.GIF" wi="72" he="33" />即满足<img file="851568dest_path_image072.GIF" wi="256" he="40" />(18)工作于额定制冷量<img file="198236dest_path_image073.GIF" wi="72" he="32" />时,所对应的输入电功率为<img file="423812dest_path_image074.GIF" wi="72" he="32" />,满足:<img file="514128dest_path_image075.GIF" wi="408" he="89" />(19)同时,蓄冰槽中增加的有效蓄冷量为<img file="502901dest_path_image076.GIF" wi="96" he="32" />,对于蓄冰状态<img file="754890dest_path_image047.GIF" wi="72" he="31" />需满足约束:<img file="467763dest_path_image077.GIF" wi="437" he="128" />(20)即只有在谷时段才能处于蓄冰状态,而且蓄冰状态必须连续工作,即启动后直到满足所需的蓄冰量为止;(c.3)释冰工况:释冰状态时输出的冷功率<img file="96190dest_path_image078.GIF" wi="88" he="40" />满足约束:<img file="210908dest_path_image079.GIF" wi="280" he="128" />(21)其中,<img file="633799dest_path_image080.GIF" wi="80" he="40" />为蓄冰槽的容量,释冰状态消耗的电功率<img file="833967dest_path_image081.GIF" wi="80" he="40" />为一固定值,同时满足约束条件:<img file="266085dest_path_image082.GIF" wi="165" he="20" />(22)即谷时段不可以工作于释冰状态,蓄冰槽:蓄冰槽在t时的蓄冷量<img file="969730dest_path_image083.GIF" wi="64" he="24" />满足如下约束:谷时段:<img file="829102dest_path_image084.GIF" wi="197" he="21" /><img file="250987dest_path_image085.GIF" wi="219" he="20" />(23)峰时段:<img file="486796dest_path_image086.GIF" wi="439" he="80" />(24),其中,参数p为耗散系数;(3)基于改进粒子群技术的微电网运行与调度计算方法:(3.1)获取冷电需求、环境参数、微网供能设备参数以及能源价格;(3.2)设定粒子数目、最大迭代数k、初始半径和搜索空间上下限,其中,设置k的初始值为1;(3.3)设定粒子<img file="310527dest_path_image087.GIF" wi="104" he="40" />和<img file="153849dest_path_image088.GIF" wi="64" he="32" />的初始位置,设定粒子<img file="312298dest_path_image089.GIF" wi="104" he="32" />和<img file="831092dest_path_image090.GIF" wi="64" he="24" />表示位置变化值的速度;(3.4)对粒子进行校正以满足步骤(2)中的耦合约束要求;(3.5)进入下一步迭代;(3.6)基于粒子<img file="758597dest_path_image091.GIF" wi="88" he="32" />、<img file="710503dest_path_image092.GIF" wi="80" he="35" />、能量平衡公式及元件的性能,计算所有的变量;(3.7)计算每个粒子的适应度的值;(3.8)更新粒子的最佳位置,所述最佳位置包括个体最佳位置和全局最佳位置;(3.9)如果所有的粒子的位置在之前设定的迭代次数内始终保持不变,返回到步骤(3.3);(3.10)如果达到设定的迭代最大值,则执行步骤(3.11),否则,更新粒子的速度与位置,转到步骤(3.4);(3.11)根据上次迭代的粒子的最小适应度值,计算所有的变量值。
地址 100015 北京市朝阳区将台路5号院5号楼2112室