CN116562563A - 一种基于最小惯量需求评估的电力***优化运行方法 - Google Patents
一种基于最小惯量需求评估的电力***优化运行方法 Download PDFInfo
- Publication number
- CN116562563A CN116562563A CN202310499340.2A CN202310499340A CN116562563A CN 116562563 A CN116562563 A CN 116562563A CN 202310499340 A CN202310499340 A CN 202310499340A CN 116562563 A CN116562563 A CN 116562563A
- Authority
- CN
- China
- Prior art keywords
- representing
- frequency
- new energy
- unit
- power
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000005457 optimization Methods 0.000 title claims abstract description 30
- 238000011156 evaluation Methods 0.000 title claims abstract description 8
- 230000008859 change Effects 0.000 claims abstract description 21
- 230000008569 process Effects 0.000 claims abstract description 17
- 230000004044 response Effects 0.000 claims abstract description 15
- 238000010248 power generation Methods 0.000 claims abstract description 5
- 208000011580 syndromic disease Diseases 0.000 claims description 8
- 230000009471 action Effects 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 238000012795 verification Methods 0.000 claims description 4
- 239000003245 coal Substances 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000011160 research Methods 0.000 claims description 3
- 238000010977 unit operation Methods 0.000 claims description 3
- 230000010485 coping Effects 0.000 claims description 2
- 230000002776 aggregation Effects 0.000 claims 1
- 238000004220 aggregation Methods 0.000 claims 1
- 229940060587 alpha e Drugs 0.000 claims 1
- 108010074506 Transfer Factor Proteins 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012804 iterative process Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000013139 quantization Methods 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于最小惯量需求评估的电力***优化运行方法,包括如下步骤:1)搭建频率响应模型,2)基于新能源发电的高阶不确定性对预测误差的可行域建模,3)构建两阶段分布鲁棒优化调度模型并分解分主、子问题进行迭代求解,主、子问题通过割集连接,4)满足收敛条件时,将结束进程并输出优化结果,改进后的方法相比常规的电力***调度方法,频率变化率降低,拟合误差更小,在兼顾鲁棒性和经济性的同时能够合理地安排***调度计划,提高***的惯性响应能力,使电力***更加安全可靠的运行。
Description
技术领域
本发明涉及一种基于最小惯量需求评估的电力***优化运行方法,属于电力***优化调度技术领域。
背景技术
随着可再生能源的发电量在电网中的占比不断提高,新能源出力呈现出的不确定性以及新能源机组对电网表现出的低惯量特征给电网稳定运行带来了挑战。因此,为了提高电网在受到故障冲击时的稳定运行能力,对新能源电力***的惯量水平进行评估至关重要。
传统的调度方法主要以经济性最优为目标,仅考虑电力***中新能源发电的高渗透率要求难以保障***的安全运行。因此,为在保障***频率稳定需求的前提下提升***经济运行水平,增强***安全稳定运行能力,本发明提出了一种基于最小惯量需求评估的电力***优化运行方法,该方法充分考虑了新能源的不确定性和动态频率变化,建立两阶段分布鲁棒优化调度模型,在考虑最小惯量需求的基础上得到***的最佳运行结果,实现电力***规划的经济效益最优。
发明内容
本发明提供一种基于最小惯量需求评估的电力***优化运行方法,用来克服***惯量减小而导致的***频率偏差过大的缺陷,提高电力***的安全稳定运行能力。
为了解决上述技术问题,本发明提供了如下的技术方案:
本发明公开了一种基于最小惯量需求评估的电力***优化运行方法,包括如下步骤:
1)搭建频率响应模型,测算***应对有功扰动情况下的最小惯量需求,建立最小惯量需求约束、频率最低点约束和准稳态频率约束。引入机械功率增益kGc、kWc和调速器时间常数TGc,提出了如图1所示的***频率响应等值模型:
式中,PGi表示火电机组i的功率,PWi表示新能源机组j的功率,vi,t表示时段t火电机组i的运行状态,Nm、Nwc表示***中具备调频能力的火电机组和新能源机组数量,n表示***的机组总数,且Nm+Nwc≤n,PNi为第i个机组的功率,TGi表示火电机组i的调速器时间常数,ki表示第i个火电机组的增益系数。
当***遭遇大功率的扰动时,频率变化率的最大值由功率缺额和***总惯量决定,各时段t内***最小惯量需求如下:
式中,表示***的频率变化率,Δfmax表示频率变化率的最大值,ΔPctg(t)表示时段t的功率缺额,f0表示***初始频率,Hsys,t表示时段t的***惯量,Hsys,min表示***的最小惯量需求,vj,t表示时段t新能源机组j的运行状态,NG表示火电机组的数量,NW表示新能源机组的数量,Hi,t表示时段t火电机组i的惯量,Hj,t表示时段t新能源机组j的惯量。
2)为了保证动态频率的最大偏移量不超过***设定的低频减载动作整定值,所构建的最低频率约束为:
式中,表示***设定的低频减载动作整定值,Δω(tz)表示频率最大偏移量。
当频率到达准稳态时,***频率接近于稳定值,此时假设该时刻频率变化率为0,得到准稳态频率约束:
式中,Δfss表示准稳态频率偏差,ΔPL表示负荷功率变化量,ΔPG表示火电机组功率变化量,ΔPW表示新能源机组功率变化量,Pd表示***的总负荷,D表示负荷阻尼系数,表示准稳态允许的最大频率偏差。
3)为考虑新能源发电的高阶不确定性,基于正态云模型的新能源预测误差的可行域建模,云的确定度与概率之间的映射关系如下:
式中:p(xi)表示概率,μ(xi)表示定性概念X量化表示的确定度,s表示新能源预测误差为x的确定度的数量,α表示检验确定度与概率一致性的常数,且α∈[0,1],随着α的增大,二者一致性也逐渐增大;M表示从正态云模型抽样的云滴个数。
为使得求取的新能源预测误差概率分布值在合理的范围内波动,需对新能源场景的概率分布值进行限制:
式中,pu表示场景u发生的概率,U表示场景数量,表示云模型得到的场景u发生的概率初始值,θ1、θ∞表示概率允许偏差限值。
{pu}需满足的置信度如下所示:
式中,α1、α∞表示概率分布值的置信度。
4)构建的两阶段分布鲁棒优化调度模型将日前电力***调度问题作为第一阶段,考虑最小惯量需求约束并制定机组运行计划方案;第二阶段基于日前电力***调度计划,考虑新能源的不确定性并建立新能源预测误差概率模糊集,研究***的鲁棒经济运行问题,建立***的频率响应全过程约束。
5)将两阶段分布鲁棒优化调度模型进行分解,其中主问题为考虑最小惯量约束的***优化调度问题,其目标函数为:
式中,sui表示机组的启动成本,sdi表示机组的停机成本,zi,t表示机组i在时段t的启动状态,1表示启动,0表示不启动;ui,t表示机组i在时段t的停运状态,1表示停运,0表示不停运;η表示子问题的目标函数值,T和NG分别表示调度时间和火电机组的数量。当不满足最小惯量需求约束时需要改变机组启停方式来保证有足够的惯量支撑,返回的优化割为:
式中,为不满足最小惯量需求约束时的***惯量。
6)两阶段分布鲁棒优化调度模型分解得到的子问题为新能源不确定性和频率动态稳定校验的优化问题,目标函数如下:
式中,ai、bi和ci表示机组i的煤耗成本函数的系数,U表示新能源预测误差离散场景的数量,Pi,t表示时段t火电机组的实际出力,Pu表示场景u发生的概率,为满足可行域的随机值;NW表示新能源机组的数量cj,t表示新能源的惩罚价格系数,表示时段t第j个新能源场的弃风弃光量。
子问题在主问题给定的变量下,寻找运行时新能最恶劣的概率分布下的弃风弃光量,并返回给主问题供下一次迭代使用,且为上式提供上界值。为检测第一阶段变量可行性,需要检验***的安全稳定运行问题,引入松弛变量ξn,t构建可行性校验子问题并生成可行割:
vi,tPi min≤Pi,t≤vi,tPi max
Pi,t-Pi,t-1≤ri up·vi,t-1+Pi max·(1-vi,t-1)
Pi,t-1-Pi,t≤ri dn·vi,t+Pi max·(1-vi,t)
ξl,t≥0
-γ1,lt-γ2,lt≤1
-1≤γ3,t≤1
式中,ηC表示可行性校验子问题的目标函数值,γ1,lt、γ2,lt、γ3,t表示约束式对偶变量,ξl,t表示引入的松弛变量,表示火电机组节点到线路l的功率转移因子,/>表示新能源机组节点到线路l的功率转移因子,/>表示线路最大的传输容量,Pi max、Pi min表示机组i允许的最大、最小输出功率,ri up、ri dn表示机组i的向上、向下爬坡速率,/>表示新能源场在时段t的出力值。
7)在迭代过程中,校验子问题增加的可行性割集为:
若通过可行性检验,则进一步求解子问题:
若不能满足收敛条件则需要返回最优割集式到主问题,修正当前***调度结果:
8)主、子问题通过割集连接,当结果满足收敛条件时,将结束进程停止迭代并输出最优解,否则增加最优割集重新求解***优化调度模型的主问题。
本发明所达到的有益效果是:在考虑新能源的不确定性以及最小惯量需求约束和动态频率约束的基础上,建立的含高比例新能源电力***的两阶段分布鲁棒优化调度模型能够协调同步惯量和虚拟惯量控制作用,限制***频率最大偏移量与准稳态频率的变化,可以降低扰动后的频率变化率,在确保各时段***惯量充裕性和动态频率响应能力的前提下合理地安排***的调度计划,降低了传统鲁棒的保守性,提高了***的惯性响应能力。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。在附图中:
图1是本发明的***频率响应等值模型。
图2是本发明的算法流程图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
实施例1
一种基于最小惯量需求评估的电力***优化运行方法,包括如下步骤:
步骤1)搭建频率响应模型,测算***应对有功扰动情况下的最小惯量需求,建立最小惯量需求约束、频率最低点约束和准稳态频率约束。引入机械功率增益kGc、kWc和调速器时间常数TGc,提出了如图1所示的***频率响应等值模型:
式中,PGi表示火电机组i的功率,PWi表示新能源机组j的功率,vi,t表示时段t火电机组i的运行状态,Nm、Nwc表示***中具备调频能力的火电机组和新能源机组数量,n表示***的机组总数,且Nm+Nwc≤n,PNi为第i个机组的功率,TGi表示火电机组i的调速器时间常数,Ki表示第i个火电机组的增益系数。
当***遭遇大功率的扰动时,频率变化率的最大值由功率缺额和***总惯量决定,各时段t内***最小惯量需求如下:
式中,表示***的频率变化率,Δfmax表示频率变化率的最大值,ΔPctg(t)表示时段t的功率缺额,f0表示***初始频率,Hsys,t表示时段t的***惯量,Hsys,min表示***的最小惯量需求,vj,t表示时段t新能源机组j的运行状态,NG表示火电机组的数量,NW表示新能源机组的数量,Hi,t表示时段t火电机组i的惯量,Hj,t表示时段t新能源机组j的惯量,vi,t表示时段t火电机组i的运行状态。
步骤2)为了保证动态频率的最大偏移量不超过***设定的低频减载动作整定值,所构建的频率最低点约束为:
式中,表示***设定的低频减载动作整定值,Δω(tz)表示频率最大偏移量。
当频率到达准稳态时,***频率接近于稳定值,此时假设该时刻频率变化率为0,得到准稳态频率约束:
式中,Δfss表示准稳态频率偏差,ΔPL表示负荷功率变化量,ΔPG表示火电机组功率变化量,ΔPW表示新能源机组功率变化量,Pd表示***的总负荷,D表示负荷阻尼系数,表示准稳态允许的最大频率偏差。
步骤3)为考虑新能源发电的高阶不确定性,基于正态云模型的新能源预测误差的可行域建模,云的确定度与概率之间的映射关系如下:
式中:p(xi)表示概率,μ(xi)表示定性概念X量化表示的确定度,s表示新能源预测误差为x的确定度的数量,α表示检验确定度与概率一致性的常数,且α∈[0,1],随着α的增大,二者一致性也逐渐增大;M表示从正态云模型抽样的云滴个数。
为使得求取的新能源预测误差概率分布值在合理的范围内波动,提出基于正态云模型构建新能源预测误差模糊集的方法是从正态云模型抽样M个云滴,根据确定度与概率转换得到每个预测误差场景的初始概率值,考虑到各离散场景初始概率分布的不确定性,以初始概率分布值为中心,综合范数条件为约束对新能源场景的概率分布值进行限制。
式中,pu表示场景u发生的概率,U表示场景数量,表示云模型得到的场景u发生的概率初始值,θ1、θ∞表示概率允许偏差限值。
{pu}需满足的置信度如下所示:
式中,α1、α∞表示概率分布值的置信度。
步骤4)构建的两阶段分布鲁棒优化调度模型将日前电力***调度问题作为第一阶段,考虑最小惯量需求约束并制定机组运行计划方案;第二阶段基于日前电力***调度计划,考虑新能源的不确定性并建立新能源预测误差概率模糊集,研究***的鲁棒经济运行问题,建立***的频率响应全过程约束。
步骤5)将两阶段分布鲁棒优化调度模型进行分解,其中主问题为考虑最小惯量约束的***优化调度问题,其目标函数为:
式中,sui表示机组的启动成本,sdi表示机组的停机成本,zi,t表示机组i在时段t的启动状态,1表示启动,0表示不启动;T和NG分别表示调度时间和火电机组的数量;ui,t表示机组i在时段t的停运状态,1表示停运,0表示不停运;η表示子问题的目标函数值。当不满足最小惯量需求约束时需要改变机组启停方式来保证有足够的惯量支撑,返回的优化割为:
式中,为不满足最小惯量需求约束时的***惯量。
步骤6)两阶段分布鲁棒优化调度模型分解得到的子问题为新能源不确定性和频率动态稳定校验的优化问题,目标函数如下:
式中,ai、bi和ci表示机组i的煤耗成本函数的系数,U表示新能源预测误差离散场景的数量,Pi,t表示时段t火电机组的实际出力,Pu表示场景u发生的概率,为满足可行域的随机值;NW表示新能源机组的数量,cj,t表示新能源的惩罚价格系数,表示时段t的弃风弃光量。
子问题在主问题给定的变量下,寻找运行时新能最恶劣的概率分布下的弃风弃光量,并返回给主问题供下一次迭代使用,且为上式提供上界值。为检测第一阶段变量可行性,需要检验***的安全稳定运行问题,引入松弛变量ξn,t构建可行性校验子问题并生成可行割:
vi,tPi min≤Pi,t≤vi,tPi max
Pi,t-Pi,t-1≤ri up·vi,t-1+Pi max·(1-vi,t-1)
Pi,t-1-Pi,t≤ri dn·vi,t+Pi max·(1-vi,t)
ξl,t≥0
-γ1,lt-γ2,lt≤1
-1≤γ3,t≤1
式中,ηC表示可行性校验子问题的目标函数值,γ1,lt、γ2,lt、γ3,t表示约束式的对偶变量,ξl,t表示引入的松弛变量,表示火电机组节点到线路l的功率转移因子,/>表示新能源机组节点到线路l的功率转移因子,/>表示线路最大的传输容量,Pi max、Pi min表示机组i允许的最大、最小输出功率,ri up、ri dn表示机组i的向上、向下爬坡速率,/>表示新能源场在时段t的出力值。
步骤7)在迭代过程中,校验子问题增加的可行性割集为:
若通过可行性检验,则进一步求解子问题:
若不能满足收敛条件则需要返回最优割集式到主问题,修正当前***调度结果:
步骤8)主、子问题通过割集连接,当结果满足收敛条件时,将结束进程停止迭代并输出最优解,否则增加最优割集重新求解***优化调度模型的主问题。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于最小惯量需求评估的电力***优化运行方法,其特征在于,包括如下步骤:
步骤1)搭建频率响应模型,测算***应对有功扰动情况下的最小惯量需求,建立最小惯量需求约束、频率最低点约束和准稳态频率约束;
步骤2)基于新能源发电的高阶不确定性对预测误差的可行域建模;
步骤3)构建两阶段分布鲁棒优化调度模型并分解为主、子问题进行迭代求解,主问题为考虑最小惯量需求约束的***优化调度问题,子问题为新能源不确定性和频率动态稳定校验的优化问题,主、子问题通过割集连接;
步骤4)满足收敛条件时,将结束进程并输出优化结果。
2.根据权利要求1所述的一种基于最小惯量需求评估的电力***优化运行方法,其特征在于,所述步骤1)中的搭建频率响应模型并建立最小惯量需求约束、频率最低点约束和准稳态频率约束包括以下步骤:
21)引入火电机组和新能源机组的机械功率增益kGc、kWc和调速器时间常数TGc,利用参数等值聚合的方法,提出***频率响应等值模型:
式中,PGi表示火电机组,PWi表示新能源机组j的功率,vi,t表示时段t火电机组i的运行状态,Nm、Nwc分别表示***中具备调频能力的火电机组和新能源机组数量,n表示***的机组总数,且Nm+Nwc≤n,PNi为第i个机组的功率,TGi表示火电机组i的调速器时间常数,ki表示第i个火电机组的增益系数;
22)当***遭遇大功率的扰动时,频率变化率的最大值由功率缺额和***总惯量决定,各时段t内***最小惯量需求如下:
式中,表示***的频率变化率,Δfmax表示频率变化率的最大值,ΔPctg(t)表示时段t的功率缺额,f0表示***初始频率,Hsys,t表示时段t的***惯量,Hsys,min表示***的最小惯量需求,vj,t表示时段t新能源机组j的运行状态,NG表示火电机组的数量,NW表示新能源机组的数量,Hi,t表示时段t火电机组i的惯量,Hj,t表示时段t新能源机组j的惯量;
23)为了保证动态频率的最大偏移量不超过***设定的低频减载动作整定值,所构建的最低频率约束为:
式中,表示***设定的低频减载动作整定值,Δω(tz)表示频率最大偏移量;
24)当频率到达准稳态时,***频率接近于稳定值,此时假设该时刻频率变化率为0,得到准稳态频率约束:
式中,Δfss表示准稳态频率偏差,ΔPL表示负荷功率变化量,ΔPG表示火电机组功率变化量,ΔPW表示新能源机组功率变化量,Pd表示***的总负荷,D表示负荷阻尼系数,表示准稳态允许的最大频率偏差。
3.根据权利要求2所述的一种基于最小惯量需求评估的电力***优化运行方法,其特征在于,所述步骤2)中,新能源预测误差的可行域建模的方法如下:
31)根据云模型理论,“新能源功率预测误差”属于定性概念可以表示为x,而x则为定性概念的一次量化表示,μ(xi)表示定性概念X量化表示的确定度;(x,μ(xi))表示定性概念的一次具体实现,称为一个云滴,云的确定度与概率之间的映射关系如下:
式中:p(xi)表示概率,μ(xi)表示定性概念X量化表示的确定度,s表示新能源预测误差为x的确定度的数量,M表示从正态云模型抽样的云滴个数,α表示检验确定度与概率一致性的常数,且α∈[0,1],随着α的增大,二者一致性也逐渐增大;
32)为使得求取的新能源预测误差概率分布值在合理的范围内波动,需对新能源场景的概率分布值进行限制:
式中,pu表示场景u发生的概率,U表示场景数量,表示云模型得到的场景u发生的概率初始值;θ1、θ∞表示概率允许偏差限值;
{pu}需满足的置信度如下所示:
式中,α1、α∞表示概率分布值的置信度。
4.根据权利要求3所述的一种基于最小惯量需求评估的电力***优化运行方法,其特征在于,所述步骤3)中,将两阶段分布鲁棒优化调度模型分解为主、子问题进行迭代求解的具体步骤如下:
41)两阶段分布鲁棒优化调度模型将日前电力***调度问题作为第一阶段,考虑最小惯量需求约束并制定机组运行计划方案;第二阶段基于日前电力***调度计划,考虑新能源的不确定性并建立新能源预测误差概率模糊集,研究***的鲁棒经济运行问题,建立***的频率响应全过程约束;
42)将两阶段分布鲁棒优化调度模型进行分解,其中主问题为考虑最小惯量约束的***优化调度问题,其目标函数为:
式中,sui表示机组的启动成本,sdi表示机组的停机成本,zi,t表示机组i在时段t的启动状态,1表示启动,0表示不启动;ui,t表示机组i在时段t的停运状态,1表示停运,0表示不停运;η表示子问题的目标函数值,T和NG分别表示调度时间和火电机组的数量;当不满足最小惯量需求约束时需要改变机组启停方式来保证有足够的惯量支撑,返回的优化割为:
式中,为不满足最小惯量需求约束时的***惯量;
43)两阶段分布鲁棒优化调度模型分解得到的子问题为新能源不确定性和频率动态稳定校验的优化问题,目标函数如下:
式中,ai、bi和ci表示机组i的煤耗成本函数的系数,U表示新能源预测误差离散场景的数量,Pi,t表示时段t火电机组的实际出力,Pu表示场景u发生的概率,为满足可行域的随机值;NW表示新能源机组的数量,cj,t表示新能源的惩罚价格系数,表示时段t第j个新能源场的弃风弃光量;
44)在迭代过程中,校验子问题增加的可行性割集为:
若通过可行性检验,则进一步求解子问题:
若不能满足收敛条件则需要返回最优割集式到主问题,修正当前***调度结果:
式中,ηC、ηs表示可行性校验子问题的目标函数值,ξl,t表示引入的松弛变量,其中n为样本数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310499340.2A CN116562563A (zh) | 2023-05-05 | 2023-05-05 | 一种基于最小惯量需求评估的电力***优化运行方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310499340.2A CN116562563A (zh) | 2023-05-05 | 2023-05-05 | 一种基于最小惯量需求评估的电力***优化运行方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116562563A true CN116562563A (zh) | 2023-08-08 |
Family
ID=87494052
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310499340.2A Pending CN116562563A (zh) | 2023-05-05 | 2023-05-05 | 一种基于最小惯量需求评估的电力***优化运行方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116562563A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117335449A (zh) * | 2023-10-30 | 2024-01-02 | 华北电力大学 | 一种高比例新能源电力***最低调频能力需求计算方法 |
-
2023
- 2023-05-05 CN CN202310499340.2A patent/CN116562563A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117335449A (zh) * | 2023-10-30 | 2024-01-02 | 华北电力大学 | 一种高比例新能源电力***最低调频能力需求计算方法 |
CN117335449B (zh) * | 2023-10-30 | 2024-03-15 | 华北电力大学 | 一种高比例新能源电力***最低调频能力需求计算方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107944757A (zh) | 电力互联***可再生能源消纳能力评估分析方法 | |
CN106786799B (zh) | 一种直流联络线功率阶梯化发电计划优化方法 | |
CN116191473B (zh) | 一种考虑随机-极端扰动下的一次调频备用优化方法 | |
CN112381375B (zh) | 一种基于潮流分配矩阵的电网经济运行域快速生成方法 | |
CN112886645B (zh) | 一种基于氢能超高比例的新能源电力***运行模拟方法 | |
CN116562563A (zh) | 一种基于最小惯量需求评估的电力***优化运行方法 | |
CN113659620A (zh) | 基于动态频率约束的水风互补发电***日前调度方法 | |
CN114465246A (zh) | 一种计及调速器死区与限幅环节的机组组合优化方法 | |
CN115660412A (zh) | 考虑源荷波动及再调度的调度风险量化评估方法及装置 | |
CN117436773B (zh) | 一种含可中断负荷的独立微网资源规划方法和*** | |
CN115313422B (zh) | 主动支撑主网频率的电-气综合能源***运行控制方法 | |
CN110336308B (zh) | 一种基于机会约束的主动配电网经济调度方法 | |
CN108306319B (zh) | 一种新能源微网中储能配置优化方法及*** | |
CN114971069A (zh) | 一种计及源荷不确定性及频率安全的电力***区间优化调度方法 | |
CN113221358B (zh) | 基于可靠性参数的电-气耦合***的备用出力优化方法 | |
CN112436510B (zh) | 一种风-光-火特高压直流外送调度方法及*** | |
Wang et al. | Multi-objective optimal dispatch of wind-integrated power system based on distributed energy storage | |
CN115619159A (zh) | 一种智慧园区能源消耗分析调度方法和*** | |
CN112600202B (zh) | 计及新能源随机性的含可控移相器电网最优潮流计算方法 | |
CN114389262A (zh) | 一种弹性环境下基于鲁棒优化的区域电网调度方法 | |
CN111404195B (zh) | 一种基于智能网关的含分布式电源微电网的调度方法 | |
CN116316865A (zh) | 一种高比例新能源区域电力***全环节协调规划优化方法 | |
CN114465226A (zh) | 一种电力***多级备用获取联合优化模型的建立方法 | |
CN117767306B (zh) | 一种利用灵活性资源短时校正网供负荷偏差的方法、***及存储介质 | |
CN114400716B (zh) | 一种水电站机组扩机的容量选取方法 |
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 |