CN107947246B - Wind power generation index distribution and increased power evaluation method considering frequency modulation and increased power - Google Patents

Wind power generation index distribution and increased power evaluation method considering frequency modulation and increased power Download PDF

Info

Publication number
CN107947246B
CN107947246B CN201711126400.7A CN201711126400A CN107947246B CN 107947246 B CN107947246 B CN 107947246B CN 201711126400 A CN201711126400 A CN 201711126400A CN 107947246 B CN107947246 B CN 107947246B
Authority
CN
China
Prior art keywords
wind power
wind
output
provincial
regulation
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.)
Active
Application number
CN201711126400.7A
Other languages
Chinese (zh)
Other versions
CN107947246A (en
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.)
Tsinghua University
North China Grid Co Ltd
Original Assignee
Tsinghua University
North China Grid Co Ltd
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 Tsinghua University, North China Grid Co Ltd filed Critical Tsinghua University
Priority to CN201711126400.7A priority Critical patent/CN107947246B/en
Publication of CN107947246A publication Critical patent/CN107947246A/en
Application granted granted Critical
Publication of CN107947246B publication Critical patent/CN107947246B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/48—Controlling the sharing of the in-phase component
    • 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)
  • Control Of Eletrric Generators (AREA)

Abstract

本发明公开了一种考虑调频增发的风电发电指标分配及增发评估方法，属于电力系统有功控制技术领域。该方法首先通过自动发电控制信息系统中相关信息，计算各省调风电虚拟机组的最大可用出力和实际出力；区域调度中心根据全网风电消纳空间和调频指令情况，计算全网风电调节指令；区域调度中心根据各省调风电虚拟机组的最大可用出力，并考虑其历史出力特性，分配各省级调度中心的风电调节指令；并根据调节指令数据和量测数据，评估各省调风电虚拟机组的调峰增发电量和调频增发电量。本发明通过引入风电机组参与电网调频以在减少弃风的同时提升系统调节性能，并考虑调节指令分配的公平性和奖惩机制，提升风电参与电网调节的可靠性。

Figure 201711126400

The invention discloses a wind power generation index allocation and additional generation evaluation method considering frequency modulation and additional generation, and belongs to the technical field of active power control of power systems. The method first calculates the maximum available output and actual output of the wind power virtual units in each province through the relevant information in the automatic power generation control information system; the regional dispatch center calculates the wind power regulation command of the whole network according to the wind power accommodation space and frequency regulation command of the whole network; The dispatch center assigns the wind power regulation commands of each provincial dispatch center according to the maximum available output of the virtual wind power units in each province and considers its historical output characteristics; and evaluates the peak regulation of the virtual wind power units in each province according to the regulation command data and measurement data. Incremental power generation and frequency modulation incremental power generation. The present invention improves the system regulation performance while reducing wind abandonment by introducing wind turbines to participate in grid frequency regulation, and considers the fairness of regulation command distribution and the reward and punishment mechanism to improve the reliability of wind power participating in grid regulation.

Figure 201711126400

Description

一种考虑调频增发的风电发电指标分配及增发评估方法A wind power generation index allocation and additional issuance evaluation method considering frequency modulation and additional issuance

技术领域technical field

本发明属于电力系统有功控制技术领域，具体涉及一种考虑调频增发的风电发电指标分配及增发评估方法。The invention belongs to the technical field of active power control of power systems, and in particular relates to a wind power generation index allocation and additional generation evaluation method considering frequency modulation and additional generation.

背景技术Background technique

作为一种技术成熟、经济效益突出的可再生能源发电，风力发电过去二十年来发展迅猛，正逐渐成为电能供给的主要来源之一。随着风力发电装机和并网规模的不断提升，其随机性和间歇性对电力系统安全稳定运行的影响逐渐显现：对于电力系统的小时级调度，风电的间歇性和不确定性将在机组组合、联络线计划、爬坡备用等方面带来困难；对于分钟级和秒级调度，风电的波动性将在系统频率控制、联络线波动调整等方面带来困难。As a renewable energy power generation with mature technology and outstanding economic benefits, wind power has developed rapidly in the past two decades and is gradually becoming one of the main sources of power supply. With the continuous increase of the installed capacity and grid-connected scale of wind power generation, the influence of its randomness and intermittency on the safe and stable operation of the power system has gradually emerged: for the hourly dispatch of the power system, the intermittency and uncertainty of wind power will be affected by the combination of units. , tie line planning, climbing backup and other aspects bring difficulties; for minute-level and second-level scheduling, the volatility of wind power will bring difficulties in system frequency control, tie line fluctuation adjustment, etc.

随着风力发电控制技术的不断发展，风电场正逐步具备接受电力系统调度、参与电力系统调度运行的能力，将风力发电作为受控资源之一纳入电网运行已成为可能。目前，我国并网风电的主要控制手段仍采用风电场跟踪日前计划，这即意味着风电场的实际运行仍为开环运行。然而，日前计划的制定依赖于风电出力预测和负荷预测，其预测误差将对风电场和电力系统的实际运行带来巨大影响，既可能导致风电资源的利用效果不佳，又可能导致电力系统调节备用容量耗尽，降低电网运行的经济性和安全性。With the continuous development of wind power control technology, wind farms are gradually being able to accept and participate in power system dispatch and operation. It has become possible to incorporate wind power as one of the controlled resources into power grid operation. At present, the main control method of grid-connected wind power in my country still adopts the wind farm tracking plan, which means that the actual operation of the wind farm is still open-loop operation. However, the formulation of the current plan relies on wind power output forecast and load forecast, and the forecast error will have a huge impact on the actual operation of wind farms and power systems, which may lead to poor utilization of wind power resources and power system adjustment. Reserve capacity is exhausted, reducing the economics and safety of grid operation.

为解决这一问题，专利“一种分调-省调两级协调的风电增发控制方法”(专利申请号：CN201610626975.4)公开了一种风电闭环控制方法，通过分调-省调的协调、风电发电指标的计算与分配，实现风力发电参与电力系统调峰的闭环运行。更进一步，在由于调峰困难造成的限风条件下，如果风电的可用出力和相对应的断面约束仍有富余，可采用适当控制策略使风电机组参与电力系统自动发电控制，在提升系统调频性能的同时，实现风电增发。但该方法未考虑风电机组的历史出力特性和调节表现，调节性能的好坏未予进行奖惩。In order to solve this problem, the patent "A control method for additional wind power generation with two-level coordination between division and provincial adjustment" (Patent Application No.: CN201610626975.4) discloses a closed-loop control method for wind power, through the coordination of division and provincial adjustment. . Calculation and distribution of wind power generation indicators to realize the closed-loop operation of wind power generation participating in the peak regulation of the power system. Furthermore, under the wind-limited conditions caused by the difficulty of peak regulation, if the available output of wind power and the corresponding section constraints are still surplus, appropriate control strategies can be adopted to enable wind turbines to participate in the automatic power generation control of the power system, so as to improve the frequency regulation performance of the system. At the same time, additional wind power generation will be realized. However, this method does not consider the historical output characteristics and regulation performance of wind turbines, and the regulation performance is not rewarded or punished.

专利“一种风电机组参与电力系统自动发电控制的协调控制方法”(专利申请号：CN201610248302.X)公开了一种限风条件下风电调频增发控制方法，该方法在区域控制中心基于模型预测控制思想计算风电机组的调节指令，在省级调度中心通过比例分配计算各风电场的调节指令。但该方法在区域调度中心基于模型预测控制的计算方法是基于优化思想，其计算实时性和工程可用性不高；其在调节指令分配时，未具体考虑分配的公平性和奖惩机制。The patent "A coordinated control method for wind turbines to participate in automatic power generation control of power systems" (patent application number: CN201610248302.X) discloses a control method for wind power frequency regulation and additional generation under wind-limited conditions, which is based on model predictive control in the regional control center The idea is to calculate the adjustment instructions of wind turbines, and calculate the adjustment instructions of each wind farm through proportional distribution in the provincial dispatch center. However, the calculation method based on model predictive control in the regional dispatch center is based on the optimization idea, and its real-time calculation and engineering availability are not high; when adjusting the command distribution, it does not specifically consider the fairness of the distribution and the reward and punishment mechanism.

发明内容SUMMARY OF THE INVENTION

本发明的目的是为克服现有技术的不足之处，提出一种考虑调频增发的风电发电指标分配及增发评估方法。本发明结合风电机组的运行状态、区域控制约束和风电机组的历史性能表现，可实现风电发电指标的合理计算、分配与事后评估。The purpose of the present invention is to overcome the deficiencies of the prior art, and propose a wind power generation index allocation and additional issuance evaluation method considering frequency modulation and additional issuance. The present invention can realize reasonable calculation, distribution and post-event evaluation of wind power generation indexes in combination with the operating state of the wind turbine, the regional control constraints and the historical performance of the wind turbine.

本发明提出一种考虑调频增发的风电发电指标分配方法，其特征在于，包括以下步骤：The present invention proposes a wind power generation index allocation method considering frequency modulation and additional generation, which is characterized by comprising the following steps:

1)通过自动发电控制信息系统AGC收集全网内所有风电场最大可用出力与实际出力，以及断面约束信息，计算各省调风电虚拟机组的最大可用出力和实际出力；计算公式如下：1) Collect the maximum available output and actual output of all wind farms in the whole network through the automatic power generation control information system AGC, as well as the section constraint information, and calculate the maximum available output and actual output of the virtual wind power units in each province. The calculation formula is as follows:

Figure BDA0001468435580000023
Figure BDA0001468435580000023

其中，

Figure BDA0001468435580000024
为省调风电虚拟机组j的最大可用出力，为省调风电虚拟机组j的实际出力，j为省调风电虚拟机组编号；PA,i为风电场i的最大可用出力；PW,i为风电场i的实际出力；和Sk分别为断面k下风电最大可用出力和断面k的容量约束，k＝1,2,...,M，k为断面编号，M为断面总数；NW,k和NG,k分别为断面k下的风电场总数和断面k下的常规电厂总数；PG,i为常规机组i的实时出力；L为省调风电虚拟机组总数；in,
Figure BDA0001468435580000024
is the maximum available output of the provincial wind power virtual unit j, is the actual output of the provincial wind power virtual unit j, j is the number of the provincial wind power virtual unit; P A,i is the maximum available output of the wind farm i; P W,i is the actual output of the wind farm i; and Sk are the maximum available wind power output under section k and the capacity constraint of section k, respectively, k=1,2,...,M, k is the section number, M is the total number of sections; N W,k and N G,k are the total number of wind farms under section k and the total number of conventional power plants under section k, respectively; P G,i is the real-time output of conventional unit i; L is the total number of provincial wind power virtual units;

2)区域调度中心根据全网风电消纳空间和调频指令情况，计算全网风电调节指令；具体步骤如下：2) The regional dispatch center calculates the wind power regulation command of the whole network according to the wind power consumption space and the frequency regulation command of the whole network; the specific steps are as follows:

2-1)计算全网风电发电指标

Figure BDA0001468435580000027
2-1) Calculate the wind power generation index of the whole network
Figure BDA0001468435580000027

2-2)判断区域内风电是否参与调频；具体步骤如下：2-2) Determine whether the wind power in the region participates in frequency regulation; the specific steps are as follows:

2-2-1)判断全网风电最大可用出力是否大于全网风电发电指标

Figure BDA0001468435580000029
若
Figure BDA00014684355800000210
则说明区域内风电出力受限，进入步骤2-2-2)；否则区域内风电出力未受限，区域内风电不参与调频，方法结束；2-2-1) Determine the maximum available output of wind power in the entire network Whether it is greater than the wind power generation index of the whole network
Figure BDA0001468435580000029
like
Figure BDA00014684355800000210
It means that the output of wind power in the area is limited, and then go to step 2-2-2); otherwise, the output of wind power in the area is not limited, and the wind power in the area does not participate in frequency regulation, and the method ends;

2-2-2)判断区域内风电是否处于出力爬坡阶段：若区域内风电处于出力爬坡阶段，则区域内风电不参与调频，方法结束；若区域内风电未处于出力爬坡阶段，则区域内风电参与调频，进入步骤2-3)；2-2-2) Determine whether the wind power in the area is in the output ramping stage: if the wind power in the region is in the output ramping stage, the wind power in the region does not participate in frequency regulation, and the method ends; if the wind power in the region is not in the output ramping stage, then If the wind power in the area participates in frequency regulation, go to step 2-3);

2-3)根据全网风电消纳空间和调频指令情况，计算全网风电调节指令，计算公式如下：2-3) Calculate the wind power regulation command of the whole network according to the wind power consumption space of the whole network and the frequency regulation command. The calculation formula is as follows:

Figure BDA0001468435580000031
Figure BDA0001468435580000031

其中，

Figure BDA0001468435580000032
为全网风电调节指令，PARR为全网调频指令，α为风电参与调频的分担系数；in,
Figure BDA0001468435580000032
is the wind power regulation command of the whole network, P ARR is the frequency regulation command of the whole network, and α is the sharing coefficient of wind power participating in frequency regulation;

3)区域调度中心根据各省调风电虚拟机组的最大可用出力，并考虑各省调风电虚拟机组历史出力特性，分配各省级调度中心的风电调节指令；具体步骤如下：3) The regional dispatch center assigns the wind power adjustment instructions of the provincial dispatch centers according to the maximum available output of the virtual wind power units in each province, and considers the historical output characteristics of the virtual wind power units in each province. The specific steps are as follows:

3-1)计算各省调风电虚拟机组调节功率分配的分担权重系数；3-1) Calculate the sharing weight coefficient of the adjustment power distribution of the wind power virtual units in each province;

省调风电虚拟机组j的新能源超短期预测精度cp,j定义如下：The ultra-short-term prediction accuracy c p,j of the new energy of the provincial wind power virtual unit j is defined as follows:

Figure BDA0001468435580000033
Figure BDA0001468435580000033

式中，Ppredict,j为非限风时段省调风电虚拟机组j超短期预测出力，Pactual,j为非限风时段省调风电虚拟机组j实发出力，D为日内非限风时段预测时间点数量，ICj为省调风电虚拟机组j装机容量；In the formula, P predict,j is the ultra-short-term predicted output of the provincially-regulated wind power virtual unit j during the non-wind restricted period, P actual,j is the actual output of the provincially-regulated wind power virtual unit j during the non-wind restricted period, and D is the intraday forecast during the non-wind restricted period. The number of time points, IC j is the installed capacity of the provincial wind power virtual unit j;

省调风电虚拟机组j的调节性能得分sp,j选取区域控制的省调风电虚拟机组j的调节性能评估指标；The regulation performance score s p,j of the provincially regulated wind power virtual unit j is selected as the regulation performance evaluation index of the provincially regulated wind power virtual unit j;

分别计算省调风电虚拟机组j的超短期预测精度基准值

Figure BDA0001468435580000034
和调节性能得分基准值
Figure BDA0001468435580000035
表达式如下：Calculate the ultra-short-term prediction accuracy reference value of the provincial wind power virtual unit j respectively
Figure BDA0001468435580000034
and tuning performance score benchmarks
Figure BDA0001468435580000035
The expression is as follows:

Figure BDA0001468435580000036
Figure BDA0001468435580000036

其中，

Figure BDA0001468435580000038
为cp,j的周度或月度平均值，
Figure BDA0001468435580000039
为sp,j的周度或月度平均值；in,
Figure BDA0001468435580000038
is the weekly or monthly average of c p,j ,
Figure BDA0001468435580000039
is the weekly or monthly average of sp,j ;

省调风电虚拟机组j的调节功率分配的分担权重系数Wj计算表达式如下：The calculation expression of the sharing weight coefficient W j of the regulated power distribution of the provincial regulated wind power virtual unit j is as follows:

Figure BDA00014684355800000310
Figure BDA00014684355800000310

其中，a和b分别为预测精度指标权重系数和调节性能指标权重系数；Among them, a and b are the weight coefficient of the prediction accuracy index and the weight coefficient of the adjustment performance index, respectively;

3-2)分配各省级调度中心的风电虚拟机组的风电调节指令；3-2) Allocate wind power regulation instructions for wind power virtual units of each provincial dispatch center;

省调风电虚拟机组j的调节指令IW,j计算公式如下：The calculation formula of the regulation command I W,j of the provincial wind power virtual unit j is as follows:

Figure BDA0001468435580000041
Figure BDA0001468435580000041

本发明提出一种考虑调频增发的风电发电指标增发评估方法，其特征在于，包括以下步骤：The present invention proposes a wind power generation index additional issuance evaluation method considering frequency modulation additional issuance, which is characterized by comprising the following steps:

1)通过自动发电控制信息系统AGC收集全网内所有风电场最大可用出力与实际出力，以及断面约束信息，计算各省调风电虚拟机组的最大可用出力和实际出力；计算公式如下：1) Collect the maximum available output and actual output of all wind farms in the whole network through the automatic power generation control information system AGC, as well as the section constraint information, and calculate the maximum available output and actual output of the virtual wind power units in each province. The calculation formula is as follows:

Figure BDA0001468435580000042
Figure BDA0001468435580000042

Figure BDA0001468435580000043
Figure BDA0001468435580000043

Figure BDA0001468435580000044
Figure BDA0001468435580000044

其中，

Figure BDA0001468435580000045
为省调风电虚拟机组j的最大可用出力，
Figure BDA0001468435580000046
为省调风电虚拟机组j的实际出力，j为省调风电虚拟机组编号；PA,i为风电场i的最大可用出力；PW,i为风电场i的实际出力；
Figure BDA0001468435580000047
和Sk分别为断面k下风电最大可用出力和断面k的容量约束，k＝1,2,...,M，k为断面编号，M为断面总数；NW,k和NG,k分别为断面k下的风电场总数和断面k下的常规电厂总数；PG,i为常规机组i的实时出力；L为省调风电虚拟机组总数；in,
Figure BDA0001468435580000045
is the maximum available output of the provincial wind power virtual unit j,
Figure BDA0001468435580000046
is the actual output of the provincial wind power virtual unit j, j is the number of the provincial wind power virtual unit; P A,i is the maximum available output of the wind farm i; P W,i is the actual output of the wind farm i;
Figure BDA0001468435580000047
and Sk are the maximum available wind power output under section k and the capacity constraint of section k, respectively, k=1,2,...,M, k is the section number, M is the total number of sections; N W,k and N G,k are the total number of wind farms under section k and the total number of conventional power plants under section k, respectively; P G,i is the real-time output of conventional unit i; L is the total number of provincial wind power virtual units;

2)区域调度中心根据全网风电消纳空间和调频指令情况，计算全网风电调节指令；具体步骤如下：2) The regional dispatch center calculates the wind power regulation command of the whole network according to the wind power consumption space and the frequency regulation command of the whole network; the specific steps are as follows:

2-1)计算全网风电发电指标

Figure BDA0001468435580000048
2-1) Calculate the wind power generation index of the whole network
Figure BDA0001468435580000048

2-2)判断区域内风电是否参与调频；具体步骤如下：2-2) Determine whether the wind power in the region participates in frequency regulation; the specific steps are as follows:

2-2-1)判断全网风电最大可用出力

Figure BDA0001468435580000049
是否大于全网风电发电指标
Figure BDA00014684355800000410
若
Figure BDA00014684355800000411
则区域内风电出力受限，进入步骤2-2-2)；否则区域内风电出力未受限，区域内风电不参与调频，方法结束；2-2-1) Determine the maximum available output of wind power in the entire network
Figure BDA0001468435580000049
Whether it is greater than the wind power generation index of the whole network
Figure BDA00014684355800000410
like
Figure BDA00014684355800000411
If the output of wind power in the area is limited, go to step 2-2-2); otherwise, the output of wind power in the area is not limited, and the wind power in the area does not participate in frequency regulation, and the method ends;

2-2-2)判断区域内风电是否处于出力爬坡阶段：若区域内风电处于出力爬坡阶段，则区域内风电不参与调频，方法结束；若区域内风电未处于出力爬坡阶段，则区域内风电参与调频，进入步骤2-3)；2-2-2) Determine whether the wind power in the area is in the output ramping stage: if the wind power in the region is in the output ramping stage, the wind power in the region does not participate in frequency regulation, and the method ends; if the wind power in the region is not in the output ramping stage, then If the wind power in the area participates in frequency regulation, go to step 2-3);

2-3)根据全网风电消纳空间和调频指令情况，计算全网风电调节指令，计算公式如下：2-3) Calculate the wind power regulation command of the whole network according to the wind power consumption space of the whole network and the frequency regulation command. The calculation formula is as follows:

Figure BDA0001468435580000051
Figure BDA0001468435580000051

其中，

Figure BDA0001468435580000052
为全网风电调节指令，PARR为全网调频指令，α为风电参与调频的分担系数；in,
Figure BDA0001468435580000052
is the wind power regulation command of the whole network, P ARR is the frequency regulation command of the whole network, and α is the sharing coefficient of wind power participating in frequency regulation;

3)区域调度中心根据各省调风电虚拟机组的最大可用出力，并考虑各省调风电虚拟机组历史出力特性，分配各省级调度中心的风电调节指令；具体步骤如下：3) The regional dispatch center assigns the wind power adjustment instructions of each provincial dispatch center according to the maximum available output of the virtual wind power units in each province, and considers the historical output characteristics of the virtual wind power units in each province. The specific steps are as follows:

3-1)计算各省调风电虚拟机组调节功率分配的分担权重系数；3-1) Calculate the sharing weight coefficient of the adjustment power distribution of the wind power virtual units in each province;

省调风电虚拟机组j的新能源超短期预测精度cp,j定义如下：The ultra-short-term prediction accuracy c p,j of the new energy of the provincial wind power virtual unit j is defined as follows:

Figure BDA0001468435580000053
Figure BDA0001468435580000053

式中，Ppredict,j为非限风时段省调风电虚拟机组j超短期预测出力，Pactual,j为非限风时段省调风电虚拟机组j实发出力，D为日内非限风时段预测时间点数量，ICj为省调风电虚拟机组j装机容量；In the formula, P predict,j is the ultra-short-term predicted output of the provincially-regulated wind power virtual unit j during the non-wind restricted period, P actual,j is the actual output of the provincially-regulated wind power virtual unit j during the non-wind restricted period, and D is the intraday forecast during the non-wind restricted period. The number of time points, IC j is the installed capacity of the provincial wind power virtual unit j;

省调风电虚拟机组j的调节性能得分sp,j选取区域控制的省调风电虚拟机组j的调节性能评估指标；The regulation performance score s p,j of the provincially regulated wind power virtual unit j is selected as the regulation performance evaluation index of the provincially regulated wind power virtual unit j;

分别计算省调风电虚拟机组j的超短期预测精度基准值

Figure BDA0001468435580000054
和调节性能得分基准值
Figure BDA0001468435580000055
表达式如下：Calculate the ultra-short-term prediction accuracy reference value of the provincial wind power virtual unit j respectively
Figure BDA0001468435580000054
and tuning performance score benchmarks
Figure BDA0001468435580000055
The expression is as follows:

Figure BDA0001468435580000057
Figure BDA0001468435580000057

其中，

Figure BDA0001468435580000058
为cp,j的周度或月度平均值，
Figure BDA0001468435580000059
为sp,j的周度或月度平均值；in,
Figure BDA0001468435580000058
is the weekly or monthly average of c p,j ,
Figure BDA0001468435580000059
is the weekly or monthly average of sp,j ;

省调风电虚拟机组j的调节功率分配的分担权重系数Wj计算表达式如下：The calculation expression of the sharing weight coefficient W j of the regulated power distribution of the provincial regulated wind power virtual unit j is as follows:

Figure BDA00014684355800000510
Figure BDA00014684355800000510

其中，a和b分别为预测精度指标权重系数和调节性能指标权重系数；Among them, a and b are the weight coefficient of the prediction accuracy index and the weight coefficient of the adjustment performance index, respectively;

3-2)分配各省级调度中心的风电虚拟机组的风电调节指令；3-2) Allocate wind power regulation instructions for wind power virtual units of each provincial dispatch center;

省调风电虚拟机组j的调节指令IW,j计算公式如下：The calculation formula of the regulation command I W,j of the provincial wind power virtual unit j is as follows:

Figure BDA0001468435580000061
Figure BDA0001468435580000061

4)根据调节指令数据和量测数据，评估各省调风电虚拟机组的调峰增发电量和调频增发电量；具体步骤如下：4) According to the regulation instruction data and measurement data, evaluate the peak regulation and frequency regulation incremental power generation of the wind power virtual units in each province; the specific steps are as follows:

4-1)在时序上匹配省调风电虚拟机组的调节指令与实际出力，得到省调风电虚拟机组实际出力和省调风电虚拟机组的调节指令相关度最大时对应的响应延迟时间；4-1) Match the adjustment command and the actual output of the provincial-regulated wind power virtual unit in time sequence, and obtain the response delay time corresponding to the maximum correlation between the actual output of the provincial-regulated wind power virtual unit and the adjustment command of the provincial-regulated wind power virtual unit;

令整个连续调控时段为一个时间窗口，在该时间窗口下将省调风电虚拟机组j的实际出力沿时间轴进行横向平移，得到平移后省调风电虚拟机组j的实际出力和省调风电虚拟机组j的调节指令最大相关度对应的响应延迟时间；其中，最大相关度rj计算公式如下：Let the entire continuous regulation period be a time window, under this time window, the actual output of the provincial-regulated wind power virtual unit j is shifted laterally along the time axis, and the actual output of the provincial-regulated wind power virtual unit j after the translation and the provincial-regulated wind power virtual unit j are obtained. The response delay time corresponding to the maximum correlation degree of the adjustment command of j; wherein, the calculation formula of the maximum correlation degree r j is as follows:

Figure BDA0001468435580000062
Figure BDA0001468435580000062

其中，δj为省调风电虚拟机组j最大相关度对应的响应延迟时间，R(x,y)为求取向量x和向量y相关度的函数；IW,j(t)为时刻t的省调风电虚拟机组j调节指令，为时刻t的省调风电虚拟机组j实际出力；Among them, δ j is the response delay time corresponding to the maximum correlation degree of the provincial wind power virtual unit j, R(x, y) is the function to find the correlation degree between the vector x and the vector y; I W,j (t) is the time t Provincial adjustment of wind power virtual unit j adjustment command, is the actual output of the provincially adjusted wind power virtual unit j at time t;

4-2)计算省调风电虚拟机组参与调峰的增发电量和调频的增发电量；4-2) Calculate the incremental power generation of the provincial wind power virtual units participating in peak regulation and the incremental power generation of frequency regulation;

将省调风电虚拟机组j实际出力和省调风电虚拟机组j的模拟人工限风指令PC,j按照步骤4-1)计算得到的响应延迟时间δj进行平移，分别得到平移后的省调风电虚拟机组j的实际出力

Figure BDA0001468435580000065
和平移后的省调风电虚拟机组j的模拟人工限风指令
Figure BDA0001468435580000066
The actual output of the provincially adjusted wind power virtual unit j and the simulated artificial wind limit command P C,j of the provincial-regulated wind power virtual unit j is translated according to the response delay time δj calculated in step 4-1), and the actual output of the provincial-regulated wind power virtual unit j after the translation is obtained respectively.
Figure BDA0001468435580000065
The simulated artificial wind limit command of the shifted provincial wind power virtual unit j
Figure BDA0001468435580000066

省调风电虚拟机组j的调峰增发电量

Figure BDA0001468435580000067
和调频增发电量
Figure BDA0001468435580000068
的计算表达式分别如下：Peak shaving incremental power generation of provincial wind power virtual unit j
Figure BDA0001468435580000067
and frequency modulation to increase power generation
Figure BDA0001468435580000068
The calculation expressions are as follows:

Figure BDA0001468435580000069
Figure BDA0001468435580000069

其中，Δt为一个AGC指令周期时长。Among them, Δt is the duration of one AGC instruction cycle.

本发明的特点及有益效果在于：The characteristics and beneficial effects of the present invention are:

本发明考虑风电在限风条件下参与调频控制以实现电量增发，结合风电机组的运行状态、区域控制约束和风电机组的历史性能表现进行风电发电指标的分配，并在事后对调峰、调频两阶段增发电量进行合理评估。通过风电参与调频，本发明可以在减少弃风电量的同时，实现考虑奖惩机制的风电发电指标分配，促进风电调节性能的提升。The present invention considers that wind power participates in frequency regulation control under wind-limited conditions to achieve increased power generation, and allocates wind power generation indicators in combination with the operating state of the wind turbine, regional control constraints and the historical performance of the wind turbine, and then adjusts the peak regulation and frequency regulation after the event. Reasonable assessment of incremental power generation in stages. Through the participation of wind power in frequency regulation, the present invention can reduce wind power curtailment and at the same time realize wind power generation index allocation considering the reward and punishment mechanism, and promote the improvement of wind power regulation performance.

附图说明Description of drawings

图1是本发明的一种考虑调频增发的风电发电指标分配及增发评估方法的流程图。Fig. 1 is a flow chart of a method for allocating wind power generation indexes and evaluating additional issuance considering frequency modulation and additional issuance according to the present invention.

具体实施方式Detailed ways

本发明提出的一种考虑调频增发的风电发电指标分配及增发评估方法，下面结合附图和具体实施例进一步详细说明如下。The present invention proposes a method for allocating wind power generation indexes and evaluating additional issuance considering frequency modulation and additional issuance, which is further described in detail below with reference to the accompanying drawings and specific embodiments.

本发明提出的一种考虑调频增发的风电发电指标分配及增发评估方法，整体流程如图1所示，包括以下步骤：A method of wind power generation index allocation and additional issuance evaluation that considers frequency modulation and additional issuance proposed by the present invention, the overall process is shown in Figure 1, including the following steps:

1)通过自动发电控制信息(AGC)系统收集全网内所有风电场最大可用出力与实际出力，以及来源于AGC系统的断面约束信息，计算各省调风电虚拟机组的最大可用出力和实际出力；计算公式如下：1) Collect the maximum available output and actual output of all wind farms in the whole network through the automatic generation control information (AGC) system, as well as the section constraint information from the AGC system, and calculate the maximum available output and actual output of the virtual wind power units in each province; The formula is as follows:

Figure BDA0001468435580000071
Figure BDA0001468435580000071

Figure BDA0001468435580000073
Figure BDA0001468435580000073

其中，

Figure BDA0001468435580000074
为省调风电虚拟机组j的最大可用出力，为省调风电虚拟机组j的实际出力，j为省调风电虚拟机组编号；PA,i为风电场i的最大可用出力(单位：MW)；PW,i为风电场i的实际出力(单位：MW)；
Figure BDA0001468435580000076
和Sk分别为断面k下风电最大可用出力和断面k的容量约束，k＝1,2,...,M，k为断面编号，M为断面总数(省调可能下辖多个风电送出基地，而每个风电送出基地通常通过一个断面送出。)；NW,k和NG,k分别为断面k下的风电场总数和断面k下的常规电厂总数；PG,i为常规机组i的实时出力；L为省调风电虚拟机组总数。本发明中，全网即指区域调度中心管辖的所有电网，其由多个控制区组成；省调是负责每个控制区调度的机构。in,
Figure BDA0001468435580000074
is the maximum available output of the provincial wind power virtual unit j, is the actual output of the provincial wind power virtual unit j, j is the number of the provincial wind power virtual unit; P A,i is the maximum available output of the wind farm i (unit: MW); P W,i is the actual output of the wind farm i ( Unit: MW);
Figure BDA0001468435580000076
and Sk are the maximum available wind power output under section k and the capacity constraint of section k, respectively, k=1, 2,...,M, k is the section number, M is the total number of sections (the province may have multiple wind power transmissions under its jurisdiction) Base, and each wind power transmission base is usually sent through a section.); N W,k and N G,k are the total number of wind farms under section k and the total number of conventional power plants under section k, respectively; P G,i is the conventional unit The real-time output of i; L is the total number of provincial wind power virtual units. In the present invention, the whole network refers to all the power grids under the jurisdiction of the regional dispatching center, which is composed of multiple control areas; the provincial dispatching is the organization responsible for the dispatching of each control area.

2)区域调度中心根据全网风电消纳空间和调频指令情况，计算全网风电调节指令。具体步骤如下：2) The regional dispatch center calculates the wind power regulation commands of the whole network according to the wind power consumption space and frequency regulation commands of the whole network. Specific steps are as follows:

2-1)计算考虑调峰约束的全网风电发电指标

Figure BDA0001468435580000077
该指标的计算方法(为已有技术)应根据区域调度中心的调度机制和调控资源进行选取，其主要考虑的因素包括：区域负荷预测、区域联络线计划、区域统调火电最小可调、预留下旋转备用、抽水蓄能电站抽水发电功率等。2-1) Calculate the wind power generation index of the whole network considering peak regulation constraints
Figure BDA0001468435580000077
The calculation method of this index (which is the existing technology) should be selected according to the dispatching mechanism and control resources of the regional dispatching center. The main factors to be considered include: regional load forecasting, regional tie-line planning, regional unified thermal power minimum adjustment, and forecasting. Leave the rotating reserve, pumped storage power station pumped hydroelectric power, etc.

特别地，本实施例中，采用如下公式计算

Figure BDA0001468435580000081
In particular, in this embodiment, the following formula is used to calculate
Figure BDA0001468435580000081

Figure BDA0001468435580000082
Figure BDA0001468435580000082

其中，Vf为超短期负荷预测(超短期负荷预测通常指未来5到1小时负荷变化情况的预测)，Ptie为联络线调度计划，Vh-reg为区域统调火电最小可调出力，Vr-reg为预留下旋转备用，Vr-plant为自备电厂出力，Vpump为抽水蓄能电站抽水发电功率。Among them, V f is the ultra-short-term load forecast (ultra-short-term load forecast usually refers to the forecast of load changes in the next 5 to 1 hour), P tie is the tie line scheduling plan, V h-reg is the regional unified thermal power minimum adjustable output, V r-reg is reserved for rotating standby, V r-plant is the output of the self-provided power plant, and V pump is the pumped power of the pumped storage power station.

2-2)判断区域内风电是否参与调频，所述区域内风电由全网所有省调风电虚拟机组构成；具体步骤如下：2-2) Determine whether the wind power in the region participates in frequency regulation, and the wind power in the region is composed of all provincial-regulated wind power virtual units in the entire network; the specific steps are as follows:

2-2-1)判断区域内风电出力是否受限，即判断全网风电最大可用出力是否大于全网风电发电指标若

Figure BDA0001468435580000085
则区域内风电出力受限，进入步骤2-2-2)；否则区域内风电出力未受限，区域内风电不参与调频，方法结束；2-2-1) Determine whether the output of wind power in the area is limited, that is, determine the maximum available output of wind power in the entire network Whether it is greater than the wind power generation index of the whole network like
Figure BDA0001468435580000085
If the output of wind power in the area is limited, go to step 2-2-2); otherwise, the output of wind power in the area is not limited, and the wind power in the area does not participate in frequency regulation, and the method ends;

2-2-2)判断区域内风电是否处于出力爬坡阶段：若区域内风电处于出力爬坡阶段，则区域内风电不参与调频，方法结束；若区域内风电未处于出力爬坡阶段，则区域内风电可参与调频，进入步骤2-3)。2-2-2) Determine whether the wind power in the area is in the output climbing stage: if the wind power in the area is in the output climbing stage, the wind power in the area does not participate in frequency regulation, and the method ends; if the wind power in the area is not in the output climbing stage, then Wind power in the area can participate in frequency regulation, go to step 2-3).

特别地，本实施例采用全网风电发电指标的波动率VRW作为区域内风电是否处于出力爬坡阶段的判断依据，VRW采用如下公式进行计算：In particular, this embodiment adopts the wind power generation index of the whole network The volatility VR W is used as the basis for judging whether the wind power in the region is in the output ramping stage. VR W is calculated by the following formula:

Figure BDA0001468435580000087
Figure BDA0001468435580000087

其中，为全网发电指标在下一个时间点(即5分钟后)的预测值。若VRW≤5％，则认为区域内风电未处于出力爬坡阶段；否则认为区域内风电处于出力爬坡阶段。in, It is the predicted value of the power generation index of the whole network at the next time point (that is, 5 minutes later). If VR W ≤5%, it is considered that the wind power in the region is not in the output ramping stage; otherwise, the wind power in the region is considered to be in the output ramping stage.

2-3)根据全网风电消纳空间和调频指令情况(所有AGC系统自带计算该指令的计算功能)，计算全网风电调节指令，计算公式如下：2-3) According to the wind power consumption space of the whole network and the frequency regulation command (all AGC systems have their own calculation function to calculate the command), calculate the wind power regulation command of the whole network, and the calculation formula is as follows:

Figure BDA0001468435580000089
Figure BDA0001468435580000089

其中，

Figure BDA00014684355800000810
为全网风电调节指令(单位：MW)，PARR为全网调频指令(单位：MW)，α为风电参与调频的分担系数(取值范围为0～1，实施例取值0.6)。上式表明，全网风电的调节指令，是在现有风电总出力上叠加调频信号；同时，需保证叠加后的调节指令大于全网风电发电指标
Figure BDA00014684355800000811
以减少弃风，小于全网风电最大可用出力以保证调节性能。in,
Figure BDA00014684355800000810
is the wind power regulation command of the whole network (unit: MW), P ARR is the frequency regulation command of the whole network (unit: MW), α is the sharing coefficient of wind power participating in the frequency regulation (the value range is 0 to 1, and the value is 0.6 in the embodiment). The above formula shows that the regulation command of the whole network of wind power is to superimpose the frequency modulation signal on the total output of the existing wind power; at the same time, it is necessary to ensure that the superimposed regulation command is greater than the wind power generation index of the whole network
Figure BDA00014684355800000811
In order to reduce wind abandonment, it is less than the maximum available output of wind power in the whole network to ensure the regulation performance.

3)区域调度中心根据各省调风电虚拟机组的最大可用出力，并考虑各省调风电虚拟机组历史出力特性，分配各省级调度中心的风电调节指令。3) The regional dispatch center assigns the wind power regulation commands of each provincial dispatch center according to the maximum available output of the virtual wind power units in each province, and considering the historical output characteristics of the virtual wind power units in each province.

为了保证新能源发电指标分配过程中的公平性，并考虑各控制区在日前预测、实时响应等环节中的具体表现，可引入对新能源超短期预测、实时响应的评价指标，结合风电虚拟机组最大可用出力，求得各省调风电虚拟机组调节功率分配的分担权重，进而分配各省级调度中心的风电调节指令。具体步骤如下：In order to ensure the fairness in the allocation process of new energy power generation indicators, and consider the specific performance of each control area in the day-ahead forecast, real-time response and other links, the evaluation indicators for ultra-short-term forecast and real-time response of new energy can be introduced, combined with wind power virtual units. The maximum available output is obtained to obtain the sharing weight of the adjustment power distribution of the wind power virtual units in each province, and then assign the wind power adjustment instructions of each provincial dispatch center. Specific steps are as follows:

3-1)计算各省调风电虚拟机组调节功率分配的分担权重系数。该权重系数是基于各省调风电虚拟机组的最大可用出力，并经由对新能源超短期预测精度和调节性能得分两个指标的修正得到。3-1) Calculate the sharing weight coefficient of the adjustment power distribution of the wind power virtual units in each province. The weight coefficient is based on the maximum available output of the wind power virtual units in each province, and is obtained by revising the two indicators of new energy ultra-short-term forecast accuracy and regulation performance score.

省调风电虚拟机组j的新能源超短期预测精度cp,j定义如下：The ultra-short-term prediction accuracy c p,j of the new energy of the provincial wind power virtual unit j is defined as follows:

式中，Ppredict,j为非限风时段省调风电虚拟机组j超短期预测出力，Pactual,j为非限风时段省省调风电虚拟机组j实发出力，D为日内非限风时段省预测时间点数量(一般为288个点)，ICj为省调风电虚拟机组j装机容量。In the formula, P predict,j is the ultra-short-term predicted output of the provincial wind power virtual unit j during the non-wind restricted period, P actual,j is the actual output of the provincial and provincial wind power virtual unit j during the non-wind restricted period, and D is the intra-day non-restricted wind period. The number of provincial forecast time points (generally 288 points), IC j is the installed capacity of the provincial wind power virtual unit j.

省调风电虚拟机组j的调节性能得分sp,j通常选取区域控制的省调风电虚拟机组j的调节性能评估指标；特别地，可参考AGC考核补偿的相关指标，如“两个细则”规定的各项调节指标；The regulation performance score s p,j of the provincially regulated wind power virtual unit j is usually selected as the regulation performance evaluation index of the regionally controlled provincially regulated wind power virtual unit j; in particular, you can refer to the relevant indicators of the AGC assessment and compensation, such as the “two detailed rules”. various adjustment indicators;

为了保证公平性与实时性，可滚动计算上述两个指标，并选取两个指标的周度或月度平均值

Figure BDA0001468435580000092
以进行后续计算。为了在新能源发电指标分配的过程中更为公平地计及上述两个指标，应对其进行修正，对每个指标，可将其除以各控制区该指标的平均值作为计算指标，其计算结果可分别表示为省调风电虚拟机组j的超短期预测精度基准值
Figure BDA0001468435580000093
和调节性能得分基准值表达式如下：In order to ensure fairness and real-time performance, the above two indicators can be calculated on a rolling basis, and the weekly or monthly average of the two indicators can be selected.
Figure BDA0001468435580000092
for subsequent calculations. In order to take into account the above two indicators more fairly in the process of new energy power generation index allocation, they should be revised. For each index, it can be divided by the average value of the index in each control area as a calculation index. The results can be respectively expressed as the ultra-short-term prediction accuracy reference value of the provincial wind power virtual unit j
Figure BDA0001468435580000093
and tuning performance score benchmarks The expression is as follows:

Figure BDA0001468435580000096
Figure BDA0001468435580000096

上述计算过程为分别对两个指标进行标幺，可使得两个指标对分配结果的影响更为接近，便于后续权重系数的选取。The above calculation process is to perform the per-unit of the two indicators respectively, which can make the influence of the two indicators on the allocation result closer, and facilitate the selection of subsequent weight coefficients.

省调风电虚拟机组j的调节功率分配的分担权重系数Wj计算表达式如下：The calculation expression of the sharing weight coefficient W j of the regulated power distribution of the provincial regulated wind power virtual unit j is as follows:

Figure BDA0001468435580000101
Figure BDA0001468435580000101

其中，a和b分别为预测精度指标权重系数和调节性能指标权重系数，取值范围均为0～1。特别地，本实施例取a＝b＝0.5。Among them, a and b are the weight coefficient of the prediction accuracy index and the weight coefficient of the adjustment performance index respectively, and the value range is 0 to 1. In particular, this embodiment takes a=b=0.5.

3-2)分配各省级调度中心的风电虚拟机组的风电调节指令。3-2) Allocate wind power regulation instructions for wind power virtual units of each provincial dispatch center.

省调风电虚拟机组j的调节指令IW,j计算公式如下：The calculation formula of the regulation command I W,j of the provincial wind power virtual unit j is as follows:

Figure BDA0001468435580000102
Figure BDA0001468435580000102

4)根据调节指令数据和量测数据，评估各省调风电虚拟机组的调峰增发电量和调频增发电量。具体步骤如下：4) According to the regulation instruction data and measurement data, evaluate the peak regulation and frequency regulation incremental power generation of the wind power virtual units in each province. Specific steps are as follows:

4-1)在时序上匹配省调风电虚拟机组的调节指令与实际出力，得到省调风电虚拟机组实际出力和省调风电虚拟机组的调节指令相关度最大时对应的响应延迟时间；考虑整个连续调控时段为一个时间窗口(连续调控包括仅调峰和调峰调频)，在该时间窗口下发的省调风电虚拟机组的调节指令将在某固定延迟后被风电执行。为确定该固定延迟，可将省调风电虚拟机组j的实际出力沿时间轴进行横向平移，得到平移后省调风电虚拟机组j的实际出力和省调风电虚拟机组j的调节指令最大相关度对应的响应延迟时间。其中，最大相关度rj计算公式如下：4-1) Match the regulation command and actual output of the provincial-regulated wind power virtual unit in time sequence, and obtain the response delay time corresponding to the maximum correlation between the actual output of the provincial-regulated wind power virtual unit and the adjustment command of the provincial-regulated wind power virtual unit; considering the entire continuous The regulation period is a time window (continuous regulation includes only peak regulation and peak regulation and frequency regulation), and the regulation command of the provincial regulation wind power virtual unit issued in this time window will be executed by the wind power after a certain fixed delay. In order to determine the fixed delay, the actual output of the provincial-regulated wind power virtual unit j can be shifted laterally along the time axis, and the corresponding relationship between the actual output of the provincial-regulated wind power virtual unit j and the regulation command of the provincial-regulated wind power virtual unit j after the translation is obtained response delay time. Among them, the calculation formula of the maximum correlation degree r j is as follows:

Figure BDA0001468435580000103
Figure BDA0001468435580000103

其中，δj为省调风电虚拟机组j最大相关度对应的响应延迟时间，R(x,y)为求取向量x和向量y相关度的函数。IW,j(t)为时刻t的省调风电虚拟机组j的调节指令，

Figure BDA0001468435580000104
为时刻t的各省调风电虚拟机组实际出力；Among them, δ j is the response delay time corresponding to the maximum correlation degree of the provincial wind power virtual unit j, and R(x, y) is the function to obtain the correlation degree of the vector x and the vector y. I W,j (t) is the adjustment command of the provincially adjusted wind power virtual unit j at time t,
Figure BDA0001468435580000104
is the actual output of the virtual wind power units in each province at time t;

4-2)计算省调风电虚拟机组参与调峰的增发电量和调频的增发电量。4-2) Calculate the incremental power generation of the provincial wind power virtual unit participating in peak regulation and the incremental power generation of frequency regulation.

将省调风电虚拟机组j实际出力

Figure BDA0001468435580000105
和省调风电虚拟机组j的模拟人工限风指令PC,j(PC,j可选取限风时段最低的风电发电指标)按照计算得到的响应延迟时间δj进行平移，分别得到平移后的省调风电虚拟机组j的实际出力和平移后的省调风电虚拟机组j的模拟人工限风指令
Figure BDA0001468435580000107
以用于计算调频增发电量。The actual output of the provincially adjusted wind power virtual unit j
Figure BDA0001468435580000105
The simulated artificial wind restriction command P C,j of the provincial wind power virtual unit j (P C, j can be selected as the wind power generation index with the lowest wind restriction period) is translated according to the calculated response delay time δ j , and the translated values are obtained respectively. Actual output of provincial wind power virtual unit j The simulated artificial wind limit command of the shifted provincial wind power virtual unit j
Figure BDA0001468435580000107
It is used to calculate the frequency modulation incremental power generation.

省调风电虚拟机组j的调峰增发电量和调频增发电量的计算表达式分别如下：Peak shaving incremental power generation of provincial wind power virtual unit j and frequency modulation to increase power generation The calculation expressions are as follows:

Figure BDA00014684355800001010
Figure BDA00014684355800001010

Figure BDA0001468435580000111
Figure BDA0001468435580000111

其中，

Figure BDA0001468435580000112
和分别为平移后的省调风电虚拟机组j的实际出力和平移后的省调风电虚拟机组j的模拟人工限风指令，Δt为一个AGC指令周期时长。in,
Figure BDA0001468435580000112
and are the actual output of the shifted provincial-regulated wind power virtual unit j and the simulated artificial wind restriction command of the shifted provincial-regulated wind power virtual unit j, and Δt is the duration of an AGC command cycle.

通过增发评估，可实现调峰、调频增发电量的分配分拆计算，以作为对风电机组进行补偿的基础。Through the evaluation of additional generation, the distribution and split calculation of peak regulation and frequency regulation incremental power generation can be realized as the basis for compensating wind turbines.

Claims (6)

1. A wind power generation index distribution method considering frequency modulation and power increase is characterized by comprising the following steps:
1) the maximum available output and the actual output of all wind power plants in the whole network and section constraint information are collected through an automatic generation control information system AGC, and the maximum available output and the actual output of each provincial dispatching wind power virtual machine set are calculated; the calculation formula is as follows:
Figure FDA0002194351090000011
Figure FDA0002194351090000012
wherein,
Figure FDA0002194351090000014
in order to save and adjust the maximum available output of the wind turbine virtual unit j,
Figure FDA0002194351090000015
the actual output of the provincial dispatching wind power virtual machine set j is obtained, and j is the serial number of the provincial dispatching wind power virtual machine set; pA,iThe maximum available output of the wind power plant i; pW,iThe actual output of the wind power plant i;and SkRespectively representing the maximum available wind power output under a section k and the capacity constraint of the section k, wherein k is 1,2, the. N is a radical ofW,kAnd NG,kRespectively the total number of the wind power plant under the section k and the total number of the conventional power plant under the section k; pG,iReal-time output of a conventional unit i; l is the total number of provincial dispatching wind power virtual units;
2) the regional dispatching center calculates a wind power regulation instruction of the whole network according to the wind power consumption space of the whole network and the frequency modulation instruction condition; the method comprises the following specific steps:
2-1) calculating the whole network wind power generation index
Figure FDA0002194351090000017
2-2) judging whether wind power in the area participates in frequency modulation; the method comprises the following specific steps:
2-2-1) judging the maximum available output of the wind power of the whole network
Figure FDA0002194351090000018
Whether the power generation index is larger than the full-grid wind power generation index
Figure FDA0002194351090000019
If it is
Figure FDA00021943510900000110
The wind power output in the region is limited, and the step 2-2-2) is carried out; otherwise, the wind-electricity output in the region is not limited, the wind-electricity in the region does not participate in frequency modulation, and the method is ended;
2-2-2) judging whether the wind power in the area is in a climbing stage: if the wind power in the region is in the output climbing stage, the wind power in the region does not participate in frequency modulation, and the method is ended; if the wind power in the region is not in the output climbing stage, the wind power in the region participates in frequency modulation, and the step 2-3 is carried out;
2-3) calculating a full-network wind power regulation instruction according to the full-network wind power consumption space and the frequency modulation instruction condition, wherein the calculation formula is as follows:
wherein,for regulating the wind power in the whole network, PARRThe method comprises the steps that a full-network frequency modulation command is adopted, and alpha is a sharing coefficient of wind power participating in frequency modulation;
3) the regional dispatching center distributes wind power regulation instructions of each provincial dispatching center according to the maximum available output of each provincial wind power regulation virtual machine set and by considering the historical output characteristics of each provincial wind power regulation virtual machine set; the method comprises the following specific steps:
3-1) calculating a sharing weight coefficient of the adjusting power distribution of each provincial wind power adjusting virtual machine set;
new energy ultrashort-term prediction precision c of provincial dispatching wind power virtual unit jp,jThe definition is as follows:
Figure FDA0002194351090000023
in the formula, Ppredict,jUltra-short-term predicted output, P, of wind turbine virtual unit j for provincial regulation of non-wind-limiting timeactual,jAdjusting the actual output of the wind power virtual unit j for the provincial wind-free time period, D is the number of predicted time points for the non-wind-limited time period in the day,ICjadjusting the installed capacity of the wind power virtual machine set j for province;
calculating the regulation performance score s of the provincial dispatching wind turbine virtual unit j according to the regulation performance evaluation index of the provincial dispatching wind turbine virtual unit j controlled by the regionp,j；
Respectively calculating the ultra-short-term prediction precision reference value of the provincial dispatching wind turbine virtual unit jAnd adjusting the performance score benchmark value
Figure FDA0002194351090000025
The expression is as follows:
Figure FDA0002194351090000026
wherein,is cp,jThe average of the degree of the week or month,
Figure FDA0002194351090000029
is s isp,jA weekly or monthly average of;
sharing weight coefficient W for adjusting power distribution of provincial wind power virtual unit jjThe calculation expression is as follows:
Figure FDA0002194351090000031
wherein, a and b are a prediction precision index weight coefficient and an adjustment performance index weight coefficient respectively;
3-2) distributing wind power adjusting instructions of the wind power virtual units of each provincial dispatching center;
province regulationAdjusting instruction I of wind power virtual unit jW,jThe calculation formula is as follows:
2. the method of claim 1, wherein the full grid wind power generation indicator in step 2-1)
Figure FDA0002194351090000033
The calculation method is as follows:
Figure FDA0002194351090000034
wherein, VfFor ultra-short term load prediction, PtieScheduling plans for links, Vh-regRegulating the minimum output of electricity, V, for a regional systemr-regFor reserve rotation, Vr-plantFor self-contained power plant output, VpumpAnd pumping water to generate power for the pumped storage power station.
3. The method according to claim 1, wherein in step 2-2-2), whether the wind power in the area is in the output climbing stage is determined by the following specific method:
calculating the whole network wind power generation index
Figure FDA0002194351090000035
The expression of fluctuation ratio of (a) is as follows:
wherein,
Figure FDA0002194351090000037
the predicted value of the power generation index of the whole network at the next time point is obtained;
if VRWIf the wind power is less than or equal to 5%, the wind power in the region is not in the stage of output climbing; otherwise, the wind power in the region is in the output climbing stage.
4. A wind power generation index increase evaluation method considering frequency modulation increase is characterized by comprising the following steps:
1) the maximum available output and the actual output of all wind power plants in the whole network and section constraint information are collected through an automatic generation control information system AGC, and the maximum available output and the actual output of each provincial dispatching wind power virtual machine set are calculated; the calculation formula is as follows:
Figure FDA0002194351090000038
Figure FDA0002194351090000039
wherein,
Figure FDA0002194351090000042
in order to save and adjust the maximum available output of the wind turbine virtual unit j,
Figure FDA0002194351090000043
the actual output of the provincial dispatching wind power virtual machine set j is obtained, and j is the serial number of the provincial dispatching wind power virtual machine set; pA,iThe maximum available output of the wind power plant i; pW,iThe actual output of the wind power plant i;
Figure FDA0002194351090000044
and SkRespectively representing the maximum available wind power output under a section k and the capacity constraint of the section k, wherein k is 1,2, the. N is a radical ofW,kAnd NG,kTotal number of wind power plants under section k and conventional power plants under section k respectivelyTotal number; pG,iReal-time output of a conventional unit i; l is the total number of provincial dispatching wind power virtual units;
2) the regional dispatching center calculates a wind power regulation instruction of the whole network according to the wind power consumption space of the whole network and the frequency modulation instruction condition; the method comprises the following specific steps:
2-1) calculating the whole network wind power generation index
Figure FDA0002194351090000045
2-2) judging whether wind power in the area participates in frequency modulation; the method comprises the following specific steps:
2-2-1) judging the maximum available output of the wind power of the whole network
Figure FDA0002194351090000046
Whether the power generation index is larger than the full-grid wind power generation indexIf it is
Figure FDA0002194351090000048
The wind power output in the region is limited, and the step 2-2-2) is carried out; otherwise, the wind-electricity output in the region is not limited, the wind-electricity in the region does not participate in frequency modulation, and the method is ended;
2-2-2) judging whether the wind power in the area is in a climbing stage: if the wind power in the region is in the output climbing stage, the wind power in the region does not participate in frequency modulation, and the method is ended; if the wind power in the region is not in the output climbing stage, the wind power in the region participates in frequency modulation, and the step 2-3 is carried out;
2-3) calculating a full-network wind power regulation instruction according to the full-network wind power consumption space and the frequency modulation instruction condition, wherein the calculation formula is as follows:
wherein,for regulating the wind power in the whole network, PARRThe method comprises the steps that a full-network frequency modulation command is adopted, and alpha is a sharing coefficient of wind power participating in frequency modulation;
3) the regional dispatching center distributes wind power regulation instructions of each provincial dispatching center according to the maximum available output of each provincial wind power regulation virtual machine set and by considering the historical output characteristics of each provincial wind power regulation virtual machine set; the method comprises the following specific steps:
3-1) calculating a sharing weight coefficient of the adjusting power distribution of each provincial wind power adjusting virtual machine set;
new energy ultrashort-term prediction precision c of provincial dispatching wind power virtual unit jp,jThe definition is as follows:
Figure FDA0002194351090000051
in the formula, Ppredict,jUltra-short-term predicted output, P, of wind turbine virtual unit j for provincial regulation of non-wind-limiting timeactual,jAdjusting the actual output of the wind power virtual unit j for the non-wind-limiting time period province, and D is the predicted time point number of the non-wind-limiting time period in the day, ICjAdjusting the installed capacity of the wind power virtual machine set j for province;
calculating the regulation performance score s of the provincial dispatching wind turbine virtual unit j according to the regulation performance evaluation index of the provincial dispatching wind turbine virtual unit j controlled by the regionp,j；
Respectively calculating the ultra-short-term prediction precision reference value of the provincial dispatching wind turbine virtual unit j
Figure FDA0002194351090000052
And adjusting the performance score benchmark value
Figure FDA0002194351090000053
The expression is as follows:
Figure FDA0002194351090000055
wherein,
Figure FDA0002194351090000056
is cp,jThe average of the degree of the week or month,
Figure FDA0002194351090000057
is s isp,jA weekly or monthly average of;
sharing weight coefficient W for adjusting power distribution of provincial wind power virtual unit jjThe calculation expression is as follows:
Figure FDA0002194351090000058
wherein, a and b are a prediction precision index weight coefficient and an adjustment performance index weight coefficient respectively;
3-2) distributing wind power adjusting instructions of the wind power virtual units of each provincial dispatching center;
adjusting instruction I of provincial wind power virtual unit jW,jThe calculation formula is as follows:
4) evaluating peak-load capacity and frequency-modulation capacity of each provincial wind power modulation virtual machine set according to the regulation instruction data and the measurement data; the method comprises the following specific steps:
4-1) matching the adjusting instruction and the actual output of the provincial dispatching wind power virtual unit in time sequence to obtain the corresponding response delay time when the correlation degree of the actual output of the provincial dispatching wind power virtual unit and the adjusting instruction of the provincial dispatching wind power virtual unit is maximum;
setting the whole continuous regulation and control time period as a time window, and transversely translating the actual output of the provincial wind power virtual unit j along a time axis under the time window to obtain the translated actual output of the provincial wind power virtual unit j and the response delay time corresponding to the maximum correlation degree of the regulation instruction of the provincial wind power virtual unit j; wherein the maximum phaseDegree of closeness rjThe calculation formula is as follows:
wherein, deltajFor response delay time corresponding to the maximum correlation degree of the provincial dispatching wind power virtual unit j, R (x, y) is a function for solving the correlation degree of the vector x and the vector y; i isW,j(t) is an adjusting instruction of the provincial wind power virtual unit j at the moment t,
Figure FDA0002194351090000062
adjusting the actual output of the wind turbine virtual unit j for the province at the moment t;
4-2) calculating the power increasing quantity of the provincial wind power and the wind power virtual machine group participating in peak shaving and the power increasing quantity of frequency modulation;
will economize on the power of wind turbine generator system j actual output
Figure FDA0002194351090000063
Simulated artificial wind limiting instruction P of provincial dispatching wind power virtual unit jC,jResponse delay time delta calculated according to step 4-1)jTranslating to respectively obtain the actual output of the translated provincial dispatching wind turbine virtual unit j
Figure FDA0002194351090000064
And the translated simulation artificial wind limiting instruction of the provincial dispatching wind turbine virtual unit j
Figure FDA0002194351090000065
Peak regulation and power generation increasing quantity of provincial wind power generation virtual machine set j
Figure FDA0002194351090000066
Sum frequency modulation power generation quantity increase
Figure FDA0002194351090000067
Respectively as follows:
Figure FDA0002194351090000068
Figure FDA0002194351090000069
where Δ t is an AGC command cycle duration.
5. The method of claim 4, wherein the full-grid wind power generation index in step 2-1)
Figure FDA00021943510900000610
The calculation method is as follows:
Figure FDA00021943510900000611
wherein, VfFor ultra-short term load prediction, PtieScheduling plans for links, Vh-regRegulating the minimum output of electricity, V, for a regional systemr-regFor reserve rotation, Vr-plantFor self-contained power plant output, VpumpAnd pumping water to generate power for the pumped storage power station.
6. The method according to claim 4, wherein the step 2-2-2) of determining whether the wind power in the area is in the output climbing stage is as follows:
calculating the whole network wind power generation indexThe expression of fluctuation ratio of (a) is as follows:
Figure FDA0002194351090000072
wherein,
Figure FDA0002194351090000073
the predicted value of the power generation index of the whole network at the next time point is obtained;
if VRWIf the wind power is less than or equal to 5%, the wind power in the region is not in the stage of output climbing; otherwise, the wind power in the region is in the output climbing stage.
CN201711126400.7A 2017-11-15 2017-11-15 Wind power generation index distribution and increased power evaluation method considering frequency modulation and increased power Active CN107947246B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711126400.7A CN107947246B (en) 2017-11-15 2017-11-15 Wind power generation index distribution and increased power evaluation method considering frequency modulation and increased power

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711126400.7A CN107947246B (en) 2017-11-15 2017-11-15 Wind power generation index distribution and increased power evaluation method considering frequency modulation and increased power

Publications (2)

Publication Number Publication Date
CN107947246A CN107947246A (en) 2018-04-20
CN107947246B true CN107947246B (en) 2020-01-10

Family

ID=61932208

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711126400.7A Active CN107947246B (en) 2017-11-15 2017-11-15 Wind power generation index distribution and increased power evaluation method considering frequency modulation and increased power

Country Status (1)

Country Link
CN (1) CN107947246B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110429666B (en) * 2019-08-22 2021-08-27 国网冀北电力有限公司电力科学研究院 Power generation control method and system based on electricity abandoning balance
CN112821469B (en) * 2021-03-09 2022-10-25 中国南方电网有限责任公司 Day-ahead power generation scheduling optimization method and device based on frequency modulation absorption domain analysis
CN113346475B (en) * 2021-04-16 2022-09-02 国电南瑞科技股份有限公司 Configurable wind-solar absorption index distribution method and device based on newly added object
CN114530888A (en) * 2022-03-15 2022-05-24 国网山东省电力公司经济技术研究院 Source-load balance-based new energy consumption capacity analysis method and system for power system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104037805A (en) * 2014-06-19 2014-09-10 清华大学 Power generation margin distribution method for photovoltaic power station based on power grid security constraints
CN104600747A (en) * 2015-01-21 2015-05-06 西安交通大学 Operation optimizing method capable of coordinating operation risk and wind energy consumption of power system
CN105914779A (en) * 2016-04-20 2016-08-31 国家电网公司 Coordination control method for wind turbine generator set participating in automatic power generation control of power system
CN106160000A (en) * 2016-08-02 2016-11-23 国电南瑞科技股份有限公司 Adjust for a kind of point and save the wind-powered electricity generation additional issue control method adjusting two-stage to coordinate

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10230242B2 (en) * 2014-10-10 2019-03-12 Edge Electrons Limited Maximizing energy savings by utilizing conservation voltage reduction with adaptive voltage control and peak demand reduction at point of use
WO2016082070A1 (en) * 2014-11-24 2016-06-02 Abb Technology Ltd Method for black starting wind turbine, wind farm, and restoring wind farm and wind turbine, wind farm using the same

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104037805A (en) * 2014-06-19 2014-09-10 清华大学 Power generation margin distribution method for photovoltaic power station based on power grid security constraints
CN104600747A (en) * 2015-01-21 2015-05-06 西安交通大学 Operation optimizing method capable of coordinating operation risk and wind energy consumption of power system
CN105914779A (en) * 2016-04-20 2016-08-31 国家电网公司 Coordination control method for wind turbine generator set participating in automatic power generation control of power system
CN106160000A (en) * 2016-08-02 2016-11-23 国电南瑞科技股份有限公司 Adjust for a kind of point and save the wind-powered electricity generation additional issue control method adjusting two-stage to coordinate

Also Published As

Publication number Publication date
CN107947246A (en) 2018-04-20

Similar Documents

Publication Publication Date Title
CN102075014B (en) Large grid real-time scheduling method for accepting access of wind power
CN108092324B (en) AGC control system and control method for wind power participating in peak shaving frequency modulation
CN105046395B (en) Method for compiling day-by-day rolling plan of power system containing multiple types of new energy
CN103824235B (en) Method for making power generation absorption scheme for new energy power stations based on price bidding
Yi et al. Energy storage capacity optimization of wind-energy storage hybrid power plant based on dynamic control strategy
CN111277005A (en) Multi-time-scale scheduling method for multi-source power system considering source-load coordination optimization
CN109636674B (en) Large-scale hydropower station group monthly transaction electric quantity decomposition and checking method
CN107947246B (en) Wind power generation index distribution and increased power evaluation method considering frequency modulation and increased power
CN104915737A (en) Coordinating economic dispatching method for power system containing large-scale wind power based on risk control
CN105162149A (en) Fuzzy self-adaptive regulation-based output method of solar-storage system tracking power generation plan
CN105006843A (en) Multi-time-scale flexible load scheduling method for handling wind power uncertainties
CN101752903A (en) Time sequence progressive power dispatching method
CN108039737A (en) One introduces a collection net lotus coordinated operation simulation system
CN106655279B (en) The grid-connected active power real-time control method in power plant of meter and electricity trading program
CN108520315A (en) A real-time control method for power grid active power considering the constraints of medium and long-term transactions and spot transactions
CN107947244A (en) The grid-connected power controls in advance method of new energy power station
CN108054790B (en) Real-time optimal control method for active power of wind and solar power generation cluster based on successive approximation of predicted output
CN108964121B (en) Wind, light and water real-time control method considering water and power planning and power target in day before water and power
CN106160000A (en) Adjust for a kind of point and save the wind-powered electricity generation additional issue control method adjusting two-stage to coordinate
CN116260181A (en) Multi-time-scale scheduling method considering large-scale renewable energy access
CN117767349A (en) Optimal scheduling and coordination control method and system for wind-solar new energy power station
CN108879657A (en) Balance of electric power and ener optimization method based on wind-powered electricity generation volume metering
CN109726894B (en) New energy active instruction calculation method for guaranteeing spot transaction and medium and long-term electric quantity
CN117498298A (en) A microgrid scheduling method based on double-layer planning theory
CN111030101B (en) A method and system for linkage regulation of clean energy consumption based on diversified big data

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