CN105529747B - 一种合理化均匀再分配电网分区风电接纳功率的建模方法 - Google Patents
一种合理化均匀再分配电网分区风电接纳功率的建模方法 Download PDFInfo
- Publication number
- CN105529747B CN105529747B CN201610003998.XA CN201610003998A CN105529747B CN 105529747 B CN105529747 B CN 105529747B CN 201610003998 A CN201610003998 A CN 201610003998A CN 105529747 B CN105529747 B CN 105529747B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- wind
- power
- powered electricity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 230000005611 electricity Effects 0.000 title claims abstract description 100
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000005540 biological transmission Effects 0.000 claims abstract description 23
- 238000005457 optimization Methods 0.000 claims description 27
- 239000000203 mixture Substances 0.000 claims description 8
- 238000009472 formulation Methods 0.000 claims description 7
- 238000009987 spinning Methods 0.000 claims description 6
- 230000009194 climbing Effects 0.000 claims description 3
- 230000005619 thermoelectricity Effects 0.000 claims description 3
- 238000010977 unit operation Methods 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 3
- 238000010276 construction Methods 0.000 description 3
- 238000013439 planning Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000010835 comparative analysis Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 229910002056 binary alloy Inorganic materials 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001550 time effect Effects 0.000 description 1
- 238000012384 transportation and delivery Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- H02J3/386—
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Wind Motors (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种合理化均匀再分配电网各分区风电接纳功率的数学建模方法,属于电力调度领域,包括以下步骤:A:建立传统含风电的电力系统年度计划制定线性模型;B:建立合理化均匀再分配电网各分区风电接纳功率模型,并对各分区风电接纳功率进行优化,具体包括各分区风电均匀再分配目标函数、各时间断面前后两次分配风电接纳总量相等约束、常规机组功率约束、区域负荷平衡约束、区域间线路传输容量约束、合理均匀分配风电功率约束。本发明优化结果更加符合实际电力系统调度运行情况,可为调度员对各分区风电接纳能力判断分析提供最直观的数据支持。
Description
技术领域
本发明公开了一种合理化均匀再分配电网各分区风电接纳功率的数学建模方法,用于省级电网的年运行方案优化,属于电力调度领域。
背景技术
作为一种随机波动电源,在大规模并网的背景下,受地区电网调峰能力与电网输送能力的制约,出现了较大规模的“弃风”现象,甚至掩盖了风电的经济、环境效益。为最大限度提高风电上网电力电量,必须从年度、月度、日前、日内和实时多个时间尺度制定风电调度计划。其中年度计划制定时间尺度长,可以在年度运行方式中预留合理风电接纳空间,保证年度计划的可执行性,对省级电网各分区域风电的建设具有指导性意义。
现有的含风电的电力系统年度计划制定模型中,文献一《江苏电网消纳大规模风电的电源规划设计》(电力系统自动化第35卷第22期第60页)分析了文中应用概率方法和工具研究了江苏省应对千万千瓦级风电基地消纳的电源规划问题。但是,该方法未考虑实际电网结构约束,导致优化结果不能指导各分区风电的建设。文献二《基于时序仿真的风电年度计划制定方法》(电力系统自动化第38卷第11期第13页)提出一种基于时序仿真的风电年度计划方法,综合考虑风电出力特性、负荷特性、机组调峰特性、电网送出能力等因素,逐时段优化全网含风电的电力平衡,建立了用于研究省级电网年度风电计划的优化模型。但是,由于该优化模型是一个典型混合整数规划模型,该模型的线性特性只能确保求解出的全省风电接纳总量是最优的,对于每个分区风电接纳情况则与采用的具体优化算法有关,优化案例中极可能出现弃风集中在某一个分区的场景,导致这个分区风电限电率出现异常(结果满足线性约束条件)。所以,综上所述,该种风电调度运行方式不能全面反映实际电网风电调度原则,虽然所得风电接纳结果对于全网是最优的,但是不能有效指导各分区年度调度运行计划特别是风电年度计划的制定。
因此,在传统的电力系统调度运行年方案的制定模型的基础上,进一步提供一种合理化均匀再分配电网各分区风电接纳功率的建模方法具有重要意义。
发明内容
本发明所要解决的技术问题在于针对上述背景技术中存在的问题,提供一种合理化均匀再分配电网分区风电接纳功率的建模方法。该方法采用文献二《基于时序仿真的风电年度计划制定方法》(电力系统自动化第38卷第11期第13页),即采用该模型优化计算得到全网各时段机组运行状态以及各时段全网最大风电接纳总功率,即进行初次线性优化分配。然后,将这部分数据提取保存作为下一步优化的输入值,基于最小二乘法进一步对全省接纳风电在各区域进行合理分配。通过上述方法,可以消除合理化均匀再分配模型中表示机组启停的二进制变量,保证算法的求解效率。如果直接将最小二乘法应用于文献二的模型,此时模型就变成混合整数二次规划,在计算全年电力系统年度计划时,计算时间将无法估计,无法满足调度人员对计算时间的要求。经过合理化均匀再分配后,优化结果更加符合实际电力系统调度运行情况,可为调度员对各区域风电接纳能力判断分析提供最直观的数据支持。
为解决上述技术问题,本发明采用的技术方案为:一种合理化均匀再分配电网各分区风电接纳功率的建模方法,包括以下步骤:
A:建立传统含风电的电力系统年度计划制定线性模型,并基于电力系统年度计划制定模型对各分区风电接纳功率进行初次线性优化分配,其中将各时间断面网内所有机组运行状态及各时间断面全省最大接纳风电大小值提取并输入到合理化均匀再分配电网各分区风电接纳功率模型;
B:建立合理化均匀再分配电网各分区风电接纳功率模型,并对各分区风电接纳功率进行优化,具体包括各分区风电均匀再分配目标函数、各时间断面前后两次分配风电接纳总量相等约束、常规机组功率约束、区域负荷平衡约束、区域间线路传输容量约束、合理均匀分配风电功率约束。
前述的一种合理化均匀再分配电网各分区风电接纳功率的建模方法,还包括如下步骤:确定电力系统年度运行调度计划,其中包括各分区风电年度计划的制定、各分区火电年度计划的制定。制定计划时,先计算出来每个区域每个时间的风电接纳值、火电的出力值,然后叠加。
前述的一种合理化均匀再分配电网各分区风电接纳功率的建模方法,步骤A的电力系统年度计划制定模型为:
初次分配目标函数
机组优化功率约束
最小启停机时间约束
供热期供热机组出力约束
启停机逻辑状态约束
机组爬坡率约束
区域间线路传输容量约束
旋转备用约束
区域负荷平衡约束
风电功率约束
式(1‐13)中,为非负变量,表示第n个区域t时段经过初次线性优化分配接纳风电功率大小;T为常数,表示总的优化时间断面;N为常数,表示电网分区数;Pj,max、Pj,min分别为机组j的出力上限和出力下限;是二进制变量,表示第机组j在时段t的运行状态,“0”表示停运,“1”表示运行;分别为表示机组j在时段t启动、停机状态的二进制变量,为“1”表示机组正在启动,为“0”表示机组不在启动状态,为“1”表示机组正在停机,为“0”表示机组不在停机状态;kon、koff分别表示机组的最小启机、停机时间,反映了机组最小启机或停机的时间长度,不同类型的机组启停机时间参数不同;为背压机组出力大小;为抽气机组出力大小;为t时段热负荷大小;为供热机组热电耦合系数; 分别为机组j的最大上爬坡速率和下爬坡速率;为t时段第i条传输线的输电功率,此处,功率流动参考方向为:流入区域为正方向,流出区域为负方向,因此,变量Li既可以取正值,有可以取负值,分别表示功率的传输方向;Li,max和‐Li,max分别为第i条传输线传输容量上下限;第t时段所有常规机组的总功率之和;则表示第t时段的电力负荷。Pre、Nre分别为正旋转备用和负旋转备用;J为常规机组总台数;为常数,表示t时段风电理论功率大小。
前述的一种合理化均匀再分配电网各分区风电接纳功率的建模方法,步骤B中的合理化均匀再分配电网各分区风电接纳功率模型具体为:
各分区风电均匀再分配目标函数:
上式中,为非负变量,表示t时段经过合理均匀再分配后区域n接纳风电功率大小;
各时间断面前后两次分配风电接纳总量相等约束
上式使得合理化均匀再分配后各时间断面接纳风电功率总量应和初次线性优化分配接纳风电总量保持一致;
常规机组功率约束
上式中,为常数,其值为初次分配模型中求得的各机组运行状态;为常规机组出力。
此外,合理化均匀再分配电网各区域风电接纳功率的模型还应该满足以下约束:
区域负荷平衡约束
区域间线路传输容量约束
合理均匀分配风电功率约束
步骤B的模型的参数含义参见步骤A中的参数含义。每一个参数在本说明书中仅代表一种含义。
相对于现有技术,本发明具有以下有益效果:
1.本发明以传统含风电的电力系统调度运行年方案制定模型为基础进行初次线性优化分配,然后,基于最小二乘法进一步对全省接纳风电在各分区进行合理分配。优化结果更加符合实际电力系统调度运行情况,可为调度员对各分区风电接纳能力判断分析提供最直观的数据支持。
2.本发明中合理化均匀再分配模型中不含二进制优化变量,将复杂的混合整数二次规划模型转换为二次规划模型,保证了案例计算的效率,满足调度人员对计算时间的要求。
附图说明
图1是本发明方法具体实施方式的基本流程图;
图2是我国东北某省水平年3个区域负荷出力序列;
图3是我国东北某省水平年3个区域年度风电序列。
具体实施方式
下面结合附图对本发明作更进一步的说明。
图1是本发明方法具体实施方式的基本流程图;首先,含风电的电力系统调度运行年方案建模(初次线性优化分配)与文献二《基于时序仿真的风电年度计划制定方法》(电力系统自动化第38卷第11期第13页)一致,这里只做简单介绍。
1)初次分配目标函数
2)机组优化功率约束
3)最小启停机时间约束
4)供热期供热机组出力约束
5)启停机逻辑状态约束
6)机组爬坡率约束
7)区域间线路传输容量约束
8)旋转备用约束
9)区域负荷平衡约束
10)风电功率约束
式(1‐13)中,为非负变量,表示第n个区域t时段经过初次线性优化分配接纳风电功率大小;T为常数,表示总的优化时间断面;N为常数,表示电网分区数;Pj,max、Pj,min分别为机组j的出力上限和出力下限;是二进制变量,表示第机组j在时段t的运行状态,“0”表示停运,“1”表示运行;分别为表示机组j在时段t启动、停机状态的二进制变量,为“1”表示机组正在启动,为“0”表示机组不在启动状态,为“1”表示机组正在停机,为“0”表示机组不在停机状态;kon、koff分别表示机组的最小启机、停机时间,反映了机组最小启机或停机的时间长度,不同类型的机组启停机时间参数不同;为背压机组出力大小;为抽气机组出力大小;为t时段热负荷大小;为供热机组热电耦合系数; 分别为机组j的最大上爬坡速率和下爬坡速率;为t时段第i条传输线的输电功率,此处,功率流动参考方向为:流入区域为正方向,流出区域为负方向,因此,变量Li既可以取正值,有可以取负值,分别表示功率的传输方向;Li,max和‐Li,max分别为第i条传输线传输容量上下限;第t时段所有常规机组的总功率之和;则表示第t时段的电力负荷。Pre、Nre分别为正旋转备用和负旋转备用;为常数,表示t时段风电理论功率大小。Pj(t)为常规机组出力;J为常规机组总台数。
然后,将上述初次线性优化分配中各分区接纳风电功率时间序列和各时间断面各机组运行状态输出,基于最小二乘法,将风电进行合理化均匀再分配(图1中虚线框部分)。其详细模型如下所示:
各分区风电均匀再分配目标函数:
式(14)中,为非负变量,表示t时段经过合理均匀再分配后区域n接纳风电功率大小。基于最小二乘法的思想,目标函数使得各区域限电率大小趋向一致,该风电调度方式符合实际电力系统中风电调度原则。
各时间断面前后两次分配风电接纳总量相等约束
式(15)使得合理化均匀再分配后各时间断面接纳风电功率总量应和初次线性优化分配接纳风电总量保持一致,这是由于初次分配是以风电接纳能力最大作为目标函数的,求解的全省接纳风电总量是最优的。此约束条件的目的是在不改变风电接纳总量的基础上,使各分区风电接纳合理。
常规机组功率约束
式(16)中,为常数,其值为初次分配模型中求得的各机组运行状态。为常规机组出力。
此外,合理均匀再分配电网各区域风电接纳功率的模型还应该满足以下约束:
区域负荷平衡约束
区域间线路传输容量约束
合理均匀分配风电功率约束
为了测试本发明方法的有效性,应用具体实施例中的方法对我国东北某省含新能源调度模型进行了仿真验证。
由于研究区域内部分新能源送出受电网断面输送极限约束,因此,本文将该省网系统分为区域1,区域2,区域3三个区域。3个区域的水平年年度风电序列、负荷出力序列如图2‐图3所示,仿真时间步长为1小时。统调机组分布情况见表1‐3。区域1到区域3的传输容量极限为1800MW;区域3与区域2的传输容量极限为1500MW。系统正备用容量为660MW,系统负备用容量为800MW。该省供热期为每年的10月25日至次年4月15日。图2是我国东北某省水平年3个区域负荷出力序列;图3是我国东北某省水平年3个区域年度风电序列。
计算环境为曙光天阔服务器,系统配置如下:CPU为AMD Operation Processor6212 32核,内存为31.9GB,操作系统为Windows XP。
表1实际系统中凝气式机组数据表
表2实际系统中抽气式机组数据表
表3实际系统中背压式机组数据表
采用合理化均匀再分配模型对电网各区域风电接纳功率进行进一步优化,初次分配和均匀再分配结果如表4所示,两次优化计算时间如表5所示。
表4初次线性优化分配与合理化均匀再分配计算结果对比分析表
表5初次线性优化分配与合理化均匀再分配计算时间对比分析表
由表4分析可知,只采用混合整数规划,区域3限电率高达34.17%,区域2此时仅为12.56%,这与实际电网接纳情况明显不符合,计算结果不能对各区域风电建设情况提供指导。而采用合理化均匀再分配各区域风电接纳功率模型优化后,区域3风电限电率从34.17%下降到25.64%,下降了8.53个百分点;区域2限电率从12.56%上升到25.46%,上升了12.9个百分点。经过本发明优化模型调整后,全省各区域风电接纳情况(限电率)趋于一致,符合实际电网各分区接纳情况。
由表5可知,本发明模型对计算时间影响不大,可以满足调度人员实际使用需求。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (4)
1.一种合理化均匀再分配电网各分区风电接纳功率的建模方法,包括以下步骤:
A:建立传统含风电的电力系统年度计划制定线性模型,并基于电力系统年度计划制定模型对各分区风电接纳功率进行初次线性优化分配,其中将各时间断面网内所有机组运行状态及各时间断面全省最大接纳风电大小值提取并输入到步骤B建立的合理化均匀再分配电网各分区风电接纳功率模型;
B:建立合理化均匀再分配电网各分区风电接纳功率模型,并对各分区风电接纳功率进行优化,具体包括各分区风电均匀再分配目标函数、各时间断面前后两次分配风电接纳总量相等约束、常规机组功率约束、区域负荷平衡约束、区域间线路传输容量约束、合理均匀分配风电功率约束。
2.根据权利要求1所述的一种合理化均匀再分配电网各分区风电接纳功率的建模方法,其特征在于,还包括如下步骤:确定电力系统年度运行调度计划,其中包括各分区风电年度计划的制定、各分区火电年度计划的制定。
3.根据权利要求1所述的一种合理化均匀再分配电网各分区风电接纳功率的建模方法,其特征在于,步骤A的电力系统年度计划制定模型为:
初次分配目标函数
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mi>w</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
机组优化功率约束
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>min</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>X</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>X</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
最小启停机时间约束
<mrow>
<msubsup>
<mi>Y</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>o</mi>
<mi>n</mi>
</mrow>
</msub>
</munderover>
<msubsup>
<mi>Z</mi>
<mi>j</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mi>i</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<mn>1</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>Z</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>o</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msub>
</munderover>
<msubsup>
<mi>Y</mi>
<mi>j</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mi>i</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<mn>1</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
供热期供热机组出力约束
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>B</mi>
<mi>Y</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>C</mi>
<mi>j</mi>
<mi>b</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>C</mi>
<mi>j</mi>
<mi>b</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>C</mi>
<mi>Q</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>C</mi>
<mi>j</mi>
<mi>v</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
启停机逻辑状态约束
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>X</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>X</mi>
<mi>j</mi>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>Y</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Z</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>X</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>X</mi>
<mi>j</mi>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Y</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>X</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>X</mi>
<mi>j</mi>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Y</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>X</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>X</mi>
<mi>j</mi>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Z</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>X</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>X</mi>
<mi>j</mi>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Z</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
机组爬坡率约束
<mrow>
<msubsup>
<mi>P</mi>
<mi>j</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>P</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mi>j</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>&le;</mo>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>d</mi>
<mi>o</mi>
<mi>w</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
区域间线路传输容量约束
<mrow>
<mo>-</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msubsup>
<mi>L</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
旋转备用约束
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>J</mi>
</munderover>
<msub>
<mi>P</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>max</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>X</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>&GreaterEqual;</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mi>l</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>J</mi>
</munderover>
<msub>
<mi>P</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>min</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>X</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mi>l</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
区域负荷平衡约束
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mi>l</mi>
<mi>l</mi>
</mrow>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mi>w</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>L</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mi>l</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
风电功率约束
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mi>w</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mi>w</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
式(1)-(13)中,为非负变量,表示第n个区域t时段经过初次线性优化分配接纳风电功率大小;T为常数,表示总的优化时间断面;N为常数,表示电网分区数;Pj,max、Pj,min分别为机组j的出力上限和出力下限;是二进制变量,表示机组j在时段t的运行状态,“0”表示停运,“1”表示运行;分别为表示机组j在时段t启动、停机状态的二进制变量,为“1”表示机组正在启动,为“0”表示机组不在启动状态,为“1”表示机组正在停机,为“0”表示机组不在停机状态;kon、koff分别表示机组的最小启机、停机时间,反映了机组最小启机或停机的时间长度,不同类型的机组启停机时间参数不同;为背压机组出力大小;为抽气机组出力大小;为t时段热负荷大小;为供热机组热电耦合系数;分别为机组j的最大上爬坡速率和下爬坡速率;为t时段第i条传输线的输电功率,此处,功率流动参考方向为:流入区域为正方向,流出区域为负方向,因此,变量Li既可以取正值,有可以取负值,分别表示功率的传输方向;Li,max和-Li,max分别为第i条传输线传输容量上下限;为t时段所有常规机组的总功率之和;Pl t,n则表示第t时段的电力负荷;Pre、Nre分别为正旋转备用和负旋转备用;J为常规机组总台数;为常数,表示t时段风电理论功率大小。
4.根据权利要求3所述的一种合理化均匀再分配电网各分区风电接纳功率的建模方法,其特征在于,步骤B中的合理化均匀再分配电网各分区风电接纳功率模型具体为:
各分区风电均匀再分配目标函数:
<mrow>
<mi>F</mi>
<mo>=</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mi>T</mi>
</mrow>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mi>N</mi>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msubsup>
<mi>P</mi>
<mi>w</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mi>w</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<msup>
<mi>n</mi>
<mrow>
<mi>n</mi>
<mi>e</mi>
<mi>w</mi>
</mrow>
</msup>
</mrow>
</msubsup>
</mrow>
<msubsup>
<mi>P</mi>
<mi>w</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
</mrow>
</msubsup>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
上式中,为非负变量,表示t时段经过合理均匀再分配后区域n接纳风电功率大小;
各时间断面前后两次分配风电接纳总量相等约束
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mi>w</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<msup>
<mi>n</mi>
<mrow>
<mi>n</mi>
<mi>e</mi>
<mi>w</mi>
</mrow>
</msup>
</mrow>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mi>w</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
</mrow>
上式使得合理化均匀再分配后各时间断面接纳风电功率总量和初次线性优化分配接纳风电总量保持一致;
常规机组功率约束
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>min</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>X</mi>
<mi>j</mi>
<msup>
<mi>t</mi>
<mo>*</mo>
</msup>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mi>j</mi>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>X</mi>
<mi>j</mi>
<msup>
<mi>t</mi>
<mo>*</mo>
</msup>
</msubsup>
</mrow>
上式中,为常数,其值为初次分配模型中求得的各机组运行状态;为常规机组出力;
此外,合理化均匀再分配电网各区域风电接纳功率的模型还应该满足以下约束:
区域负荷平衡约束
<mrow>
<msubsup>
<mi>P</mi>
<mi>l</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mi>l</mi>
<mi>l</mi>
</mrow>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mi>w</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<msup>
<mi>n</mi>
<mrow>
<mi>n</mi>
<mi>e</mi>
<mi>w</mi>
</mrow>
</msup>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>L</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
</mrow>
区域间线路传输容量约束
<mrow>
<mo>-</mo>
<msubsup>
<mi>L</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>L</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>L</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mi>t</mi>
</msubsup>
</mrow>
合理均匀分配风电功率约束
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mi>w</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<msup>
<mi>n</mi>
<mrow>
<mi>n</mi>
<mi>e</mi>
<mi>w</mi>
</mrow>
</msup>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mi>w</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>n</mi>
<mo>*</mo>
</mrow>
</msubsup>
<mo>.</mo>
</mrow>
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610003998.XA CN105529747B (zh) | 2016-01-06 | 2016-01-06 | 一种合理化均匀再分配电网分区风电接纳功率的建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610003998.XA CN105529747B (zh) | 2016-01-06 | 2016-01-06 | 一种合理化均匀再分配电网分区风电接纳功率的建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105529747A CN105529747A (zh) | 2016-04-27 |
CN105529747B true CN105529747B (zh) | 2018-04-24 |
Family
ID=55771812
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610003998.XA Expired - Fee Related CN105529747B (zh) | 2016-01-06 | 2016-01-06 | 一种合理化均匀再分配电网分区风电接纳功率的建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105529747B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106532782B (zh) * | 2016-12-12 | 2021-09-21 | 中国电力科学研究院有限公司 | 一种提升风电利用率的热电联产机组运行优化方法及装置 |
CN109460917A (zh) * | 2018-11-08 | 2019-03-12 | 中国南方电网有限责任公司 | 一种基于分布因子和支持向量机的母线负荷预测方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103904644A (zh) * | 2014-03-26 | 2014-07-02 | 国家电网公司 | 一种基于分布式电源接入的智能变电站负荷自动分配方法 |
-
2016
- 2016-01-06 CN CN201610003998.XA patent/CN105529747B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103904644A (zh) * | 2014-03-26 | 2014-07-02 | 国家电网公司 | 一种基于分布式电源接入的智能变电站负荷自动分配方法 |
Non-Patent Citations (2)
Title |
---|
Load Balancing in Smart DC Micro-grid Using Delay Tolerant User Demands;Daud Mustafa Minhas等;《2015 IEEE International Conference on Smart Grid Communications》;20151105;第157-162页 * |
面向波动平滑、负荷跟踪和功率平衡的风电功率优化调控策略;曲直等;《中国电机工程学报》;20130605;第33卷(第16期);第47-55页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105529747A (zh) | 2016-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105305419B (zh) | 一种含有压缩空气储能的独立微网容量优化配置方法 | |
CN107276127B (zh) | 考虑联络线电力交易计划的多区域互联电力系统的风电消纳优化方法 | |
CN103151803B (zh) | 一种含风电系统机组及备用配置的优化方法 | |
Martinez-Mares et al. | A robust optimization approach for the interdependency analysis of integrated energy systems considering wind power uncertainty | |
Salimi et al. | Simultaneous operation of wind and pumped storage hydropower plants in a linearized security-constrained unit commitment model for high wind energy penetration | |
CN104600695B (zh) | 基于在线状态估计和实时调度计划的趋势潮流计算方法 | |
CN104993523B (zh) | 使含风电电网系统优化运行的抽水蓄能电站特性精确模拟方法 | |
CN111555281B (zh) | 一种电力系统灵活性资源配置仿真方法及装置 | |
CN104467029B (zh) | 一种区域小水电、风电日内联合优化调度的方法 | |
CN106505633B (zh) | 一种风光接入容量确定方法及装置 | |
CN104156887A (zh) | 一种基于电网约束的区域新能源开发规模及布局分析方法 | |
CN104578122B (zh) | 抽水蓄能电站的容量确定方法 | |
CN109728578A (zh) | 基于牛顿法求解分位数的电力系统随机动态机组组合方法 | |
CN110400096A (zh) | 一种火电融合电储能联合调峰的日前发电计划优化方法及系统 | |
CN101917024A (zh) | 安全约束调度中通用性成本空间的生成方法 | |
CN104063808A (zh) | 一种跨省送电梯级水电站群调峰调度两阶段搜索方法 | |
CN109118024A (zh) | 一种考虑火电机组多阶段状态转移的电力系统多资源调度方法 | |
CN104809545B (zh) | 一种虚拟电厂运行建模方法 | |
CN104682447A (zh) | 一种含多风电场的电力系统经济调度方法 | |
CN105243600A (zh) | 一种电网发电调节方法 | |
CN115842376A (zh) | 电力系统等效惯量趋势与安全状态评估方法、设备和介质 | |
CN105529747B (zh) | 一种合理化均匀再分配电网分区风电接纳功率的建模方法 | |
CN115114854A (zh) | 一种虚拟电厂分布式资源两级自组织优化聚合方法及系统 | |
CN110021932B (zh) | 双边参与的调峰辅助服务容量选取和经济模型构建方法 | |
CN104537445A (zh) | 一种网省两级多电源短期协调调峰方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180424 Termination date: 20210106 |