CN111769560A - 一种含光伏发电的电力系统紧急降风险控制方法 - Google Patents

一种含光伏发电的电力系统紧急降风险控制方法 Download PDF

Info

Publication number
CN111769560A
CN111769560A CN202010694254.3A CN202010694254A CN111769560A CN 111769560 A CN111769560 A CN 111769560A CN 202010694254 A CN202010694254 A CN 202010694254A CN 111769560 A CN111769560 A CN 111769560A
Authority
CN
China
Prior art keywords
power generation
generation system
photovoltaic power
photovoltaic
time period
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.)
Granted
Application number
CN202010694254.3A
Other languages
English (en)
Other versions
CN111769560B (zh
Inventor
钟璐
杨哲瀚
陈愚
欧阳金鑫
熊俊
熊小伏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University
Southwest Branch of State Grid Corp
Original Assignee
Chongqing University
Southwest Branch of State Grid Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing University, Southwest Branch of State Grid Corp filed Critical Chongqing University
Priority to CN202010694254.3A priority Critical patent/CN111769560B/zh
Publication of CN111769560A publication Critical patent/CN111769560A/zh
Application granted granted Critical
Publication of CN111769560B publication Critical patent/CN111769560B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/001Methods to deal with contingencies, e.g. abnormalities, faults or failures
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/001Methods to deal with contingencies, e.g. abnormalities, faults or failures
    • H02J3/0012Contingency detection
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/008Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381Dispersed generators
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/466Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20The dispersed energy generation being of renewable origin
    • H02J2300/22The renewable source being solar energy
    • H02J2300/24The renewable source being solar energy of photovoltaic origin
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种含光伏发电的电力系统紧急降风险控制方法,包括:确定非光伏发电系统的开停机状态和光伏发电系统的出力情况;确定非光伏发电系统备用容量在实时运行阶段的调度情况和光伏减载量;以非光伏发电系统的开停机状态、光伏发电机的组出力情况、非光伏发电系统备用容量调度情况和光伏减载量,作为调度生产中各个非光伏发电系统、光伏发电系统的控制指令。本发明使光伏发电系统具备对电网的紧急支撑能力,解决了电网运行中可能出现的功率不足问题,计及光伏减载和电力系统运行风险,在电力系统功率不足时主动调整光伏发电系统出力以对电力系统进行紧急功率支撑,在保证电网安全性的前提下,实现经济效益最大化。

Description

一种含光伏发电的电力系统紧急降风险控制方法
技术领域
本发明涉及电力系统控制领域,具体涉及一种含光伏发电的电力系统紧急降风险控制方法。
背景技术
近年来,化石能源的逐渐短缺和环境污染的加剧使得光伏等可再生能源受到了各国政府的重视。光伏发电系统在电网中的渗透率不断提高。光伏发电系统的功率随日照强度波动而不断变化,对电力系统的安全稳定运行造成了较大风险。特别是光伏发电系统的运行原理与同步发电机存在根本性差异,在扰动冲击等条件下,光伏发电系统使得电力系统的运行控制面临根本性挑战。
电力系统的功率平衡是电网稳定运行的基础。然而,电源和负荷的投切、输电线路的中断均会导致电网功率不平衡。电网不平衡时,电压和频率出现波动,可能导致电网频率和节点电压越限,甚至会造成电网振荡、解列等严重事故。当电力系统功率出现冗余时,切机是控制电力系统恢复稳定运行的最行之有效的方法。在光伏高渗透率的电力系统中,切光伏被逐渐用以代替切机来维持电力系统的功率平衡。
当电力系统功率出现不足时,在传统电力系统中,由于同步发电机短时的功率提升能力有限,因此常常采用切负荷的方式来平衡电力系统的功率。部分研究人员提出了一种基于系统加速能量的临界切机量求取方法,从而量化切机控制,构建了一种区域电网的单机带负荷模型,并配置了一种高频切机方案。但是,增大同步发电机的备用容量成本较高。特别是在含有光伏的电力系统中,由于同步发电机容量有限,可能难以满足电力系统功率紧急控制的需求。
为了最大化发电效率,光伏发电系统一般工作在最大功率跟踪(MPPT)模式。因此,在任意光照条件下,光伏发电系统输出功率已到达上限,无法为电网提供有效的紧急功率支撑。部分研究人员提出控制光伏运行于非MPPT点,从而通过减载运行来具有一定备用功率。但是,现有方法主要聚焦于光伏发电系统并网调频,关于大扰动下光伏参与电力系统紧急控制的研究还鲜见报道。特别是现有研究主要着力于确定一定光伏减载量下电力系统的运行方式,关于光伏减载量的确定还未形成成熟的方法。由于光伏减载量不仅决定于电力系统的需求,还与电力系统的运行成本息息相关,合理的光伏减载量对于电力系统的稳定和经济运行具有重大意义。
因此,如何准确确定光伏减载量,从而有效实现基于光伏发电减载运行的电力系统紧急降风险控制,成为了本领域技术人员急需解决的问题。
发明内容
针对上述现有技术的不足,本发明提供一种含光伏发电的电力系统紧急降风险控制方法,用于非光伏发电系统与光伏发电系统并网的电力系统,使光伏发电系统具备对电网的紧急支撑能力,解决电网运行中可能出现的功率不足问题,计及光伏减载和电力系统运行风险,通过建立两阶段降风险调度优化模型安排同步机和光伏发电系统出力、协调同步机备用容量和光伏减载量,在电力系统功率不足时主动调整光伏发电系统出力以对电力系统进行紧急功率支撑,在保证电网安全性的前提下,实现经济效益最大化。
为了解决上述技术问题,本发明采用了如下的技术方案:
一种含光伏发电的电力系统紧急降风险控制方法,目标电网包括多个非光伏发电系统及多个光伏发电系统,包括:
S1、建立日前调度时段优化模型,以日前发电调度成本最低为优化目标,以非光伏发电系统和光伏发电系统的有功出力以及非光伏发电系统的开停机序列为优化变量,确定非光伏发电系统的开停机状态和光伏发电系统的出力情况;
S2、基于非光伏发电系统的开停机状态和光伏发电系统的出力情况,以及节点电压越限、电网频率偏差、电压稳定裕度和功角稳定裕度,结合电网故障集计算当前电网运行风险IRisk
S3、建立实时运行时段优化模型,以实时运行成本和风险成本只和最低为优化目标,以非光伏发电系统备用容量出力和光伏减载量为优化变量,确定非光伏发电系统备用容量在实时运行阶段的调度情况和光伏减载量;
S4、以非光伏发电系统的开停机状态、光伏发电机的组出力情况、非光伏发电系统备用容量调度情况和光伏减载量,作为调度生产中各个非光伏发电系统、光伏发电系统的控制指令,实现基于光伏发电减载运行的电力系统紧急降风险控制。
优选地,步骤S1中,日前调度时段优化模型每隔ΔT间隔启动,确定未来M个ΔT时段内非光伏发电系统的开停机状态和光伏发电系统的出力情况,ΔT为日前调度时段优化模型启动间隔,日前调度时段优化模型的目标函数为:
Min.CDA
日前发电调度成本CDA的计算公式如下:
Figure BDA0002590412080000031
式中,T={1,2,...,Nt},Nt为日前调度运行的时段数;G={1,2,...,NG},NG为电网中非光伏发电系统的总数目;S={1,2,...,Ns},Ns为电网中光伏发电系统的总数目;
Figure BDA0002590412080000032
分别为第g个非光伏发电系统在时间段t的启动成本和停机成本;Cg为第g个非光伏发电系统的边际发电成本报价;
Figure BDA0002590412080000033
为第g个非光伏发电系统在时间段t内调度的有功功率,Cs为第s个光伏发电系统的边际发电成本报价;
Figure BDA0002590412080000034
为第s个光伏发电系统在时间段t内调度的有功功率;
日前调度时段优化模型的约束条件包括:
潮流平衡约束
Figure BDA0002590412080000035
Figure BDA0002590412080000036
式中,m和n分别表示电网中的第m个节点和第n个节点;δnt和δmt分别为日前调度阶段第m个节点和第n个节点在时间段t内的电压相角;Bnm为节点电纳矩阵;
Figure BDA0002590412080000037
为日前调度阶段线路m至n在时间段t内的潮流;
Figure BDA0002590412080000038
为日前调度阶段预测的负荷节点d在时间段t内的有功功率;
Figure BDA0002590412080000039
分别表示位于节点n的非光伏发电系统、光伏发电系统和负荷;m∈Ψn表示和节点n相连的节点;
机组出力约束:
Figure BDA00025904120800000311
Figure BDA00025904120800000312
式中,Pgmax、Pgmin分别为非光伏发电系统的出力上和下限;PMPPT为第s个光伏发电系统在日前调度阶段时间段t内的最大输出有功功率,ugt为第g个非光伏发电系统在时间段t内的状态,ugt为取值为0或1的整数变量,当ugt=1时,表示在时间段t内调度第g个非光伏发电系统,当ugt=0时,则第g个非光伏发电系统不被调度;
日前调度时段优化模型的约束条件包括:
爬坡率约束:
Figure BDA0002590412080000041
式中RUg和RDg分别为第g个非光伏发电系统的最大向上和向下爬坡速率;
Figure BDA0002590412080000042
为第g个非光伏发电系统在日前调度阶段时间段t-1内调度的有功功率;
线路启停机成本约束:
Figure BDA0002590412080000043
Figure BDA0002590412080000044
Figure BDA0002590412080000045
Figure BDA0002590412080000046
Figure BDA0002590412080000047
Figure BDA0002590412080000048
式中,
Figure BDA0002590412080000049
分别为第g个非光伏发电系统的启动成本和停运成本;ug(t-1)为第g个非光伏发电系统在时间段t-1内的状态;
Figure BDA00025904120800000410
表示第g个非光伏发电机组的初始状态,当
Figure BDA00025904120800000411
时,表示第g个非光伏发电系统的初始状态是开机,反之,则第g个非光伏发电系统的初始状态是停机。
优选地,步骤S2中,
Figure BDA00025904120800000412
式中,K为所考虑的电网不确定因素的总数目;Xi为第i个所考虑的电网不确定因素;P(Xi)为第i个所考虑的电网不确定因素发生的概率;S(Xi)为第i个所考虑的电网不确定因素造成的后果,所考虑的电网不确定因素包括电压越限、电网频率偏差、电压稳定裕度和功角稳定裕度;
节点电压越限后果:
Figure BDA00025904120800000413
式中,a=Va/VaN,Va为节点电压幅值;VaN为节点电压的额定值;
电网频率偏差后果:
Figure BDA0002590412080000051
式中,b=f/fN,f为电网频率;fN=50Hz;
电压稳定裕度:
Figure BDA0002590412080000052
式中,(Pi)0为负荷节点的有功初始值;Y=1,...,Z为计算步数;ΔTY为步长;
功角稳定裕度:
Figure BDA0002590412080000053
式中,ΔPe为小扰动后每个非光伏发电系统的有功增量;Δα为机端电压相位变化。
优选地,步骤S3中,实时运行阶段优化模型以ΔT′为采样步长,ΔT′<ΔT,对当前时刻的电网功率不足风险信息进行超短期预测,并根据预测结果进行最优化计算,确定系统中非光伏发电系统备用容量、光伏减载量;
实时运行时段优化模型的目标函数为:
Min.CERT+CR
式中,CERT表示实时运行期望成本,CR表示风险成本;
Figure BDA0002590412080000054
Figure BDA0002590412080000055
式中,Rgtw表示第g个非光伏发电系统在场景w和时间段t内的备用容量;Pstw为第s个光伏发电系统在场景w和时间段t内的出力;
Figure BDA0002590412080000056
为第s个光伏发电系统在场景w和时间段t内可削减的出力,β为风险成本系数;Ω为场景总数;
实时运行时段优化模型的约束条件包括:
潮流平衡约束:
Figure BDA0002590412080000057
Figure BDA0002590412080000061
式中,δntw和δmtw分别是实时运行阶段节点n和m在场景w和时间段t内的电压相角;
Figure BDA0002590412080000062
为实时运行阶段线路m至n在场景w和时间段t内的潮流;Pdtw为负荷节点d在场景w和时间段t内的有功功率;
光伏减载约束:
Figure BDA0002590412080000063
式中,
Figure BDA0002590412080000064
为光伏减载量;
备用容量约束:
-Rdg≤Rgtw≤Rug
式中,Rdg、Rug分别为第g个非光伏发电系统最大的向上和向下备用容量;
机组出力约束:
Figure BDA0002590412080000065
爬坡率约束:
Figure BDA0002590412080000066
式中,Rg(t-1)w为第g个非光伏发电机机组在场景w和时间段t-1内的备用容量。
综上所述,与现有技术相比,本发明具有以下有益效果:
1、与现有技术中未确定光伏减载量不同,本发明是在电力系统功率不足前对即将发生变化的电网电压进行预测,在系统功率不足时,通过事前主动使光伏发电系统减载运行以提高光伏发电系统的功率支撑能力,计及电网日前调度时段运行风险,通过两阶段降风险调度协调系统中同步发电机、光伏发电系统出力和同步发电机备用容量、光伏减载量从而达到降低电网运行风险的目的。
2、日前调度时段下,以各类机组发电成本最小为目标,结合系统运行约束优化计算出系统中非光伏发电系统的出力,光伏发电系统的出力。日前调度时段下确定的非光伏发电系统开停机状态和光伏发电系统的出力情况在实时运行时段下不再调整。
3、实时运行时段下,为了降低电网运行风险,以实时运行成本与风险成本之和最小为控制目标,选取相应的风险成本系数,在满足约束条件的情况下,对系统中非光伏发电系统的备用容量出力以及光伏减载量进行最优分配,最终实现在有效降低电网运行风险的前提下,实现经济效益最大化。
附图说明
为了使发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步的详细描述,其中:
图1为本发明的控制原理图;
图2为本发明公开的含光伏发电的电力系统紧急降风险控制方法的流程图;
图3为一光伏并网电力系统的结构图实例。
具体实施方式
下面结合附图对本发明作进一步的详细说明。
如图1及图2所示,本发明公开了一种含光伏发电的电力系统紧急降风险控制方法,目标电网包括多个非光伏发电系统及多个光伏发电系统,系统中非光伏发电系统集合为G={1,2,...,NG},NG为系统中非光伏发电系统的总数目;光伏发电系统集合为S={1,2,...,Ns},Ns为系统中光伏发电系统的总数目,以图3为例,目标电网内有34条输电线路、12台非光伏发电机机组,标号1至24分别为1至24节点,节点10处连接的为水力发电机组,节点21处接入一装机容量为200MW的光伏发电系统,其余机组为火力发电,所述方法包括:
S1、建立日前调度时段优化模型,以日前发电调度成本最低为优化目标,以非光伏发电系统和光伏发电系统的有功出力以及非光伏发电系统的开停机序列为优化变量,确定非光伏发电系统的开停机状态和光伏发电系统的出力情况;
步骤S1在具体实施时,日前调度时段优化模型每隔ΔT间隔启动,对未来MΔT时段内的光伏、发电和负荷数据进行预测,根据预测结果,以该时段内总发电成本最小为优化目标,确定系统中各非光伏发电机组和光伏发电系统的出力计划;
日前调度时段优化模型的目标函数为:
Min.CDA
日前发电调度成本CDA的计算公式如下:
Figure BDA0002590412080000071
式中,T={1,2,...,Nt},Nt为日前调度运行的时段数;G={1,2,...,NG},NG为电网中非光伏发电系统的总数目;S={1,2,...,Ns},Ns为电网中光伏发电系统的总数目;
Figure BDA0002590412080000081
分别为第g个非光伏发电系统在时间段t的启动成本和停机成本;Cg为第g个非光伏发电系统的边际发电成本报价;
Figure BDA0002590412080000082
为第g个非光伏发电系统在时间段t内调度的有功功率,Cs为第s个光伏发电系统的边际发电成本报价;
Figure BDA0002590412080000083
为第s个光伏发电系统在时间段t内调度的有功功率;
日前调度时段优化模型的约束条件包括:
潮流平衡约束
Figure BDA0002590412080000084
Figure BDA0002590412080000085
式中,m和n分别表示电网中的第m个节点和第n个节点;δnt和δmt分别为日前调度阶段第m个节点和第n个节点在时间段t内的电压相角;Bnm为节点电纳矩阵;
Figure BDA0002590412080000086
为日前调度阶段线路m至n在时间段t内的潮流;
Figure BDA0002590412080000087
为日前调度阶段预测的负荷节点d在时间段t内的有功功率;
Figure BDA0002590412080000088
分别表示位于节点n的非光伏发电系统、光伏发电系统和负荷;m∈Ψn表示和节点n相连的节点;
机组出力约束:
Figure BDA0002590412080000089
Figure BDA00025904120800000810
式中,Pgmax、Pgmin分别为非光伏发电系统的出力上和下限;PMPPT为第s个光伏发电系统在日前调度阶段时间段t内的最大输出有功功率,ugt为第g个非光伏发电系统在时间段t内的状态,ugt为整数变量,当ugt=1时,表示在时间段t内调度第g个非光伏发电系统,反之,则第g个非光伏发电系统不被调度;
日前调度时段优化模型的约束条件包括:
爬坡率约束:
Figure BDA00025904120800000811
式中RUg和RDg分别为第g个非光伏发电系统的最大向上和向下爬坡速率;
Figure BDA00025904120800000812
为第g个非光伏发电系统在日前调度阶段时间段t-1内调度的有功功率;
线路启停机成本约束:
Figure BDA0002590412080000091
Figure BDA0002590412080000092
Figure BDA0002590412080000093
Figure BDA0002590412080000094
Figure BDA0002590412080000095
Figure BDA0002590412080000096
式中,
Figure BDA0002590412080000097
分别为第g个非光伏发电系统的启动成本和停运成本;ug(t-1)为第g个非光伏发电系统在时间段t-1内的状态;
Figure BDA0002590412080000098
表示第g个非光伏发电机组的初始状态,当
Figure BDA0002590412080000099
时,表示第g个非光伏发电系统的初始状态是开机,反之,则第g个非光伏发电系统的初始状态是停机。
S2、基于非光伏发电系统的开停机状态和光伏发电系统的出力情况,以及节点电压越限、电网频率偏差、电压稳定裕度和功角稳定裕度,结合电网故障集计算当前电网运行风险IRisk
在步骤S2具体实施时,基于日前调度时段优化模型确定的非光伏发电机组开停机状态和光伏发电系统的出力情况,考虑节点电压越限、电网频率偏差、电压稳定裕度和功角稳定裕度,结合电网故障集计算当前电网运行风险IRisk
Figure BDA00025904120800000910
式中,K为所考虑的电网不确定因素的总数目;Xi为第i个所考虑的电网不确定因素;P(Xi)为第i个所考虑的电网不确定因素发生的概率;S(Xi)为第i个所考虑的电网不确定因素造成的后果,所考虑的电网不确定因素包括电压越限、电网频率偏差、电压稳定裕度和功角稳定裕度;
节点电压越限后果:
Figure BDA00025904120800000911
式中,a=Va/VaN,Va为节点电压幅值;VaN为节点电压的额定值;
电网频率偏差后果:
Figure BDA0002590412080000101
式中,b=f/fN,f为电网频率;fN=50Hz;
电压稳定裕度:
Figure BDA0002590412080000102
式中,(Pi)0为负荷节点的有功初始值;Y=1,...,Z为计算步数;ΔTY为步长;
功角稳定裕度:
Figure BDA0002590412080000103
式中,ΔPe为小扰动后每个非光伏发电系统的有功增量;Δα为机端电压相位变化。
本发明中,概率来源于电网实际运行时的统计数据;电压稳定裕度和功角稳定裕度本身即属于后果的。
S3、建立实时运行时段优化模型,以实时运行成本和风险成本只和最低为优化目标,以非光伏发电系统备用容量出力和光伏减载量为优化变量,确定非光伏发电系统备用容量在实时运行阶段的调度情况和光伏减载量;
步骤S3具体实施时,利用包含光伏减载电压约束和风险成本的实时运行时段优化模型,对实时运行成本和风险成本进行优化,基于步骤1中日前调度时段优化模型确定的非光伏发电机组开停机状态和光伏发电系统的出力情况和步骤2中的电网运行风险,以非光伏发电机组备用容量出力和光伏减载量为优化变量,确定非光伏发电机组备用容量在实时运行阶段的调度情况和光伏减载量;
实时运行时段优化模型的目标函数为:
Min.CERT+CR
式中,CERT表示实时运行期望成本,CR表示风险成本;
Figure BDA0002590412080000104
Figure BDA0002590412080000111
式中,Rgtw表示第g个非光伏发电系统在场景w和时间段t内的备用容量;Pstw为第s个光伏发电系统在场景w和时间段t内的出力;
Figure BDA0002590412080000112
为第s个光伏发电系统在场景w和时间段t内可削减的出力,β为风险成本系数;Ω为场景总数;
实时运行时段优化模型的约束条件包括:
潮流平衡约束:
Figure BDA0002590412080000113
Figure BDA0002590412080000114
式中,δntw和δmtw分别是实时运行阶段节点n和m在场景w和时间段t内的电压相角;
Figure BDA0002590412080000115
为实时运行阶段线路m至n在场景w和时间段t内的潮流;Pdtw为负荷节点d在场景w和时间段t内的有功功率;
光伏减载约束:
Figure BDA0002590412080000116
式中,
Figure BDA0002590412080000117
为光伏减载量;
备用容量约束:
-Rdg≤Rgtw≤Rug
式中,Rdg、Rug分别为第g个非光伏发电系统最大的向上和向下备用容量;
机组出力约束:
Figure BDA0002590412080000118
爬坡率约束:
Figure BDA0002590412080000119
式中,Rg(t-1)w为第g个非光伏发电机机组在场景w和时间段t-1内的备用容量。
S4、以非光伏发电系统的开停机状态、光伏发电机的组出力情况、非光伏发电系统备用容量调度情况和光伏减载量,作为调度生产中各个非光伏发电系统、光伏发电系统的控制指令,实现基于光伏发电减载运行的电力系统紧急降风险控制。
光伏发电系统通常工作在最大功率跟踪模式,其缺乏多余裕度对电网进行功率支撑,在电网出现功率不足问题时,不能参加电网的紧急控制。通过减载留出有功备用,可以使光伏发电系统具有对电网的紧急支撑能力。本发明中确定了光伏发电系统的减载量值,光伏发电系统以较小的功率损失为代价,换取了紧急情况下对电网的支撑能力,降低了系统的功率不足风险,提高了系统的安全性。与不考虑光伏减载情况下的电力系统进行对比,本发明所提方法更有利于系统的安全性。
本发明考虑了光伏减载参与电网紧急支撑,兼顾了传统经济性和电力系统运行风险,以风电及负荷的滚动预测值为输入变量,通过建立两阶段电力系统降风险调度模型安排非光伏发电机组出力、光伏机组出力,协调非光伏发电机组备用容量和光伏减载量,在电力系统功率不足前通过主动调整非光伏发电机组备用容量和光伏减载量保证电网功率平衡的情况下,最大限度地提高系统安全性,保障电力系统经济性。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管通过参照本发明的优选实施例已经对本发明进行了描述,但本领域的普通技术人员应当理解,可以在形式上和细节上对其作出各种各样的改变,而不偏离所附权利要求书所限定的本发明的精神和范围。

Claims (4)

1.一种含光伏发电的电力系统紧急降风险控制方法,目标电网包括多个非光伏发电系统及多个光伏发电系统,其特征在于,包括:
S1、建立日前调度时段优化模型,以日前发电调度成本最低为优化目标,以非光伏发电系统和光伏发电系统的有功出力以及非光伏发电系统的开停机序列为优化变量,确定非光伏发电系统的开停机状态和光伏发电系统的出力情况;
S2、基于非光伏发电系统的开停机状态和光伏发电系统的出力情况,以及节点电压越限、电网频率偏差、电压稳定裕度和功角稳定裕度,结合电网故障集计算当前电网运行风险IRisk
S3、建立实时运行时段优化模型,以实时运行成本和风险成本只和最低为优化目标,以非光伏发电系统备用容量出力和光伏减载量为优化变量,确定非光伏发电系统备用容量在实时运行阶段的调度情况和光伏减载量;
S4、以非光伏发电系统的开停机状态、光伏发电机的组出力情况、非光伏发电系统备用容量调度情况和光伏减载量,作为调度生产中各个非光伏发电系统、光伏发电系统的控制指令,实现基于光伏发电减载运行的电力系统紧急降风险控制。
2.如权利要求1所述的含光伏发电的电力系统紧急降风险控制方法,其特征在于,步骤S1中,日前调度时段优化模型每隔ΔT间隔启动,确定未来M个ΔT时段内非光伏发电系统的开停机状态和光伏发电系统的出力情况,ΔT为日前调度时段优化模型启动间隔,日前调度时段优化模型的目标函数为:
Min.CDA
日前发电调度成本CDA的计算公式如下:
Figure FDA0002590412070000011
式中,T={1,2,...,Nt},Nt为日前调度运行的时段数;G={1,2,...,NG},NG为电网中非光伏发电系统的总数目;S={1,2,...,Ns},Ns为电网中光伏发电系统的总数目;
Figure FDA0002590412070000012
分别为第g个非光伏发电系统在时间段t的启动成本和停机成本;Cg为第g个非光伏发电系统的边际发电成本报价;
Figure FDA0002590412070000013
为第g个非光伏发电系统在时间段t内调度的有功功率,Cs为第s个光伏发电系统的边际发电成本报价;
Figure FDA0002590412070000014
为第s个光伏发电系统在时间段t内调度的有功功率;
日前调度时段优化模型的约束条件包括:
潮流平衡约束
Figure FDA0002590412070000021
Figure FDA0002590412070000022
式中,m和n分别表示电网中的第m个节点和第n个节点;δnt和δmt分别为日前调度阶段第m个节点和第n个节点在时间段t内的电压相角;Bnm为节点电纳矩阵;
Figure FDA0002590412070000023
为日前调度阶段线路m至n在时间段t内的潮流;
Figure FDA0002590412070000024
为日前调度阶段预测的负荷节点d在时间段t内的有功功率;
Figure FDA0002590412070000025
分别表示位于节点n的非光伏发电系统、光伏发电系统和负荷;m∈Ψn表示和节点n相连的节点;
机组出力约束:
Figure FDA0002590412070000026
Figure FDA0002590412070000027
式中,Pgmax、Pgmin分别为非光伏发电系统的出力上和下限;PMPPT为第s个光伏发电系统在日前调度阶段时间段t内的最大输出有功功率,ugt为第g个非光伏发电系统在时间段t内的状态,ugt为取值为0或1的整数变量,当ugt=1时,表示在时间段t内调度第g个非光伏发电系统,当ugt=0时,则第g个非光伏发电系统不被调度;
日前调度时段优化模型的约束条件包括:
爬坡率约束:
Figure FDA0002590412070000028
式中RUg和RDg分别为第g个非光伏发电系统的最大向上和向下爬坡速率;
Figure FDA0002590412070000029
为第g个非光伏发电系统在日前调度阶段时间段t-1内调度的有功功率;
线路启停机成本约束:
Figure FDA00025904120700000210
Figure FDA00025904120700000211
Figure FDA00025904120700000212
Figure FDA0002590412070000031
Figure FDA0002590412070000032
Figure FDA0002590412070000033
式中,
Figure FDA0002590412070000034
分别为第g个非光伏发电系统的启动成本和停运成本;ug(t-1)为第g个非光伏发电系统在时间段t-1内的状态;
Figure FDA0002590412070000035
表示第g个非光伏发电机组的初始状态,当
Figure FDA0002590412070000036
时,表示第g个非光伏发电系统的初始状态是开机,反之,则第g个非光伏发电系统的初始状态是停机。
3.如权利要求2所述的含光伏发电的电力系统紧急降风险控制方法,其特征在于,步骤S2中,
Figure FDA0002590412070000037
式中,K为所考虑的电网不确定因素的总数目;Xi为第i个所考虑的电网不确定因素;P(Xi)为第i个所考虑的电网不确定因素发生的概率;S(Xi)为第i个所考虑的电网不确定因素造成的后果,所考虑的电网不确定因素包括电压越限、电网频率偏差、电压稳定裕度和功角稳定裕度;
节点电压越限后果:
Figure FDA0002590412070000038
式中,a=Va/VaN,Va为节点电压幅值;VaN为节点电压的额定值;
电网频率偏差后果:
Figure FDA0002590412070000039
式中,b=f/fN,f为电网频率;fN=50Hz;
电压稳定裕度:
Figure FDA00025904120700000310
式中,(Pi)0为负荷节点的有功初始值;Y=1,...,Z为计算步数;ΔTY为步长;
功角稳定裕度:
Figure FDA0002590412070000041
式中,ΔPe为小扰动后每个非光伏发电系统的有功增量;Δα为机端电压相位变化。
4.如权利要求3所述的含光伏发电的电力系统紧急降风险控制方法,其特征在于,步骤S3中,实时运行阶段优化模型以ΔT′为采样步长,ΔT′<ΔT,对当前时刻的电网功率不足风险信息进行超短期预测,并根据预测结果进行最优化计算,确定系统中非光伏发电系统备用容量、光伏减载量;
实时运行时段优化模型的目标函数为:
Min.CERT+CR
式中,CERT表示实时运行期望成本,CR表示风险成本;
Figure FDA0002590412070000042
Figure FDA0002590412070000043
式中,Rgtw表示第g个非光伏发电系统在场景w和时间段t内的备用容量;Pstw为第s个光伏发电系统在场景w和时间段t内的出力;
Figure FDA0002590412070000044
为第s个光伏发电系统在场景w和时间段t内可削减的出力,β为风险成本系数;Ω为场景总数;
实时运行时段优化模型的约束条件包括:
潮流平衡约束:
Figure FDA0002590412070000045
Figure FDA0002590412070000046
式中,δntw和δmtw分别是实时运行阶段节点n和m在场景w和时间段t内的电压相角;
Figure FDA0002590412070000047
为实时运行阶段线路m至n在场景w和时间段t内的潮流;Pdtw为负荷节点d在场景w和时间段t内的有功功率;
光伏减载约束:
Figure FDA0002590412070000048
式中,
Figure FDA0002590412070000051
为光伏减载量;
备用容量约束:
-Rdg≤Rgtw≤Rug
式中,Rdg、Rug分别为第g个非光伏发电系统最大的向上和向下备用容量;
机组出力约束:
Figure FDA0002590412070000052
爬坡率约束:
Figure FDA0002590412070000053
式中,Rg(t-1)w为第g个非光伏发电机机组在场景w和时间段t-1内的备用容量。
CN202010694254.3A 2020-07-17 2020-07-17 一种含光伏发电的电力系统紧急降风险控制方法 Active CN111769560B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010694254.3A CN111769560B (zh) 2020-07-17 2020-07-17 一种含光伏发电的电力系统紧急降风险控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010694254.3A CN111769560B (zh) 2020-07-17 2020-07-17 一种含光伏发电的电力系统紧急降风险控制方法

Publications (2)

Publication Number Publication Date
CN111769560A true CN111769560A (zh) 2020-10-13
CN111769560B CN111769560B (zh) 2022-01-25

Family

ID=72728489

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010694254.3A Active CN111769560B (zh) 2020-07-17 2020-07-17 一种含光伏发电的电力系统紧急降风险控制方法

Country Status (1)

Country Link
CN (1) CN111769560B (zh)

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110264289A1 (en) * 2011-05-10 2011-10-27 Harry Sawyer Control system for photovoltaic power plant
JP4998909B1 (ja) * 2012-02-08 2012-08-15 豊国工業株式会社 太陽光発電システム
JP2013013178A (ja) * 2011-06-28 2013-01-17 Hitachi Ltd 電力系統の運転制御システム、運転制御装置及び運転制御方法
CN103441535A (zh) * 2013-08-01 2013-12-11 国电南瑞科技股份有限公司 基于sced的日前发电计划光伏发电接纳能力分析方法
FR3015142A1 (fr) * 2013-12-16 2015-06-19 Inst Polytechnique Grenoble Procede de stabilisation d'un reseau electrique par delestage de charges
CN104764962A (zh) * 2012-12-05 2015-07-08 江苏省电力公司常州供电公司 分布式能源的孤岛检测、控制方法
CN105048518A (zh) * 2015-09-23 2015-11-11 武汉大学 一种光伏参与孤立电网调频的控制方法
CN105303267A (zh) * 2015-11-20 2016-02-03 武汉大学 一种考虑动态频率约束的含高渗透率光伏电源的孤立电网机组组合优化方法
CN106058917A (zh) * 2016-05-31 2016-10-26 国网山东省电力公司济南供电公司 一种计及光伏发电的电力系统动态经济调度方法
CN108258685A (zh) * 2018-01-26 2018-07-06 国网辽宁省电力有限公司 新能源电力系统在紧急调控域运行方式下的控制方法
CN109687481A (zh) * 2018-12-27 2019-04-26 国网青海省电力公司电力科学研究院 基于主从协调控制的光伏电站参与一次调频实现方法
CN110492512A (zh) * 2019-09-12 2019-11-22 合肥工业大学 一种光储联合系统中调频或调峰模式的控制方法
CN110707724A (zh) * 2019-11-26 2020-01-17 江苏方天电力技术有限公司 一种配电网无功支撑能力评估方法

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110264289A1 (en) * 2011-05-10 2011-10-27 Harry Sawyer Control system for photovoltaic power plant
JP2013013178A (ja) * 2011-06-28 2013-01-17 Hitachi Ltd 電力系統の運転制御システム、運転制御装置及び運転制御方法
JP4998909B1 (ja) * 2012-02-08 2012-08-15 豊国工業株式会社 太陽光発電システム
CN104764962A (zh) * 2012-12-05 2015-07-08 江苏省电力公司常州供电公司 分布式能源的孤岛检测、控制方法
CN103441535A (zh) * 2013-08-01 2013-12-11 国电南瑞科技股份有限公司 基于sced的日前发电计划光伏发电接纳能力分析方法
FR3015142A1 (fr) * 2013-12-16 2015-06-19 Inst Polytechnique Grenoble Procede de stabilisation d'un reseau electrique par delestage de charges
CN105048518A (zh) * 2015-09-23 2015-11-11 武汉大学 一种光伏参与孤立电网调频的控制方法
CN105303267A (zh) * 2015-11-20 2016-02-03 武汉大学 一种考虑动态频率约束的含高渗透率光伏电源的孤立电网机组组合优化方法
CN106058917A (zh) * 2016-05-31 2016-10-26 国网山东省电力公司济南供电公司 一种计及光伏发电的电力系统动态经济调度方法
CN108258685A (zh) * 2018-01-26 2018-07-06 国网辽宁省电力有限公司 新能源电力系统在紧急调控域运行方式下的控制方法
CN109687481A (zh) * 2018-12-27 2019-04-26 国网青海省电力公司电力科学研究院 基于主从协调控制的光伏电站参与一次调频实现方法
CN110492512A (zh) * 2019-09-12 2019-11-22 合肥工业大学 一种光储联合系统中调频或调峰模式的控制方法
CN110707724A (zh) * 2019-11-26 2020-01-17 江苏方天电力技术有限公司 一种配电网无功支撑能力评估方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JINXIN OUYANG等: "《Fault Overload Control Method for High-Proportion Wind Power Transmission Systems Based on Emergency Acceleration of Doubly-Fed Induction Generator》", 《 IEEE ACCESS》 *
胡斌等: "《考虑大规模光伏电站接入的电力系统旋转备用需求评估》", 《电力系统自动化》 *
魏勇等: "《计及日照强度时间周期特征的光伏并网系统风险评估方法》", 《电网技术》 *

Also Published As

Publication number Publication date
CN111769560B (zh) 2022-01-25

Similar Documents

Publication Publication Date Title
CN111293691A (zh) 基于模型预测控制的微电网多时间尺度优化调度方法
CN109768581A (zh) 用于储能电站的电网调压与动态无功支撑控制方法
CN112234655A (zh) 一种送端电网安全稳定的源网协调调峰优化方法
CN114336678B (zh) 一种基于pmu的风光储场站一次调频控制的方法
CN112886645B (zh) 一种基于氢能超高比例的新能源电力系统运行模拟方法
CN109149631A (zh) 一种考虑风光储提供柔性爬坡容量的两阶段经济调度方法
CN109992818A (zh) 大规模风电参与一次调频的机组组合模型及求解方法
CN114069688A (zh) 一种基于时序生产模拟的多电源容量布局规划方法
CN112838603A (zh) 一种风光储抽多源能源agc协调互补控制方法和装置
El-Bahay et al. Computational methods to mitigate the effect of high penetration of renewable energy sources on power system frequency regulation: a comprehensive review
CN114142536B (zh) 一种考虑容量备用的多类型机组协调方法
CN115333161A (zh) 一种绿色水厂供电系统的容量优化配置方法
CN106329574A (zh) 一种基于辐照度变化的光伏电站无功电压控制方法
Tan et al. Synergistic balancing control for low-inertia power systems with high PV penetration: Tibet as a case study
Nakamura et al. Green base station using robust solar system and high performance lithium ion battery for next generation wireless network (5G) and against mega disaster
CN111769560B (zh) 一种含光伏发电的电力系统紧急降风险控制方法
CN114498749A (zh) 主动配电系统孤岛划分和自主运行方法及系统
Nithara et al. Review on role of power electronics in integration of renewable energy sources with micro grid
Ma et al. Multi-Point Layout Planning of Multi-Energy Power Supplies Based on Time-series Production Simulation
Naresh et al. Intelligent control strategy for power management in hybrid renewable energy system
Yu et al. A Two-stage Model Predictive Control Strategy for Economical Operation of Microgrid
Cai et al. Decentralized autonomous operation for islanded photovoltaic-based DC microgrids
Kumar et al. Tie-line frequency bias control of pv-wind hybrid power system for micro grid application
Yahia et al. An EMS based on the FLC of a variable three-phase dump load Composed of PV/Fuel Cell/Battery Systems
Kadda et al. Optimization of the managed electrical energy within a hybrid renewable energy system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant