CN105048491A

CN105048491A - Multi-stage wind power accepted range calculating method based on unit combination and economic dispatching

Info

Publication number: CN105048491A
Application number: CN201510369943.6A
Authority: CN
Inventors: 戴则梅; 张凯锋; 张慧玲; 杨国强; 王颖; 高宗和; 陈佳星; 谢丽荣; 丁茂生; 韩红卫
Original assignee: Nari Technology Co Ltd; NARI Nanjing Control System Co Ltd; State Grid Ningxia Electric Power Co Ltd
Current assignee: Nari Technology Co Ltd; NARI Nanjing Control System Co Ltd; State Grid Ningxia Electric Power Co Ltd
Priority date: 2015-06-29
Filing date: 2015-06-29
Publication date: 2015-11-11

Abstract

本发明公开了基于机组组合和经济调度的多阶段风电接纳范围计算方法，包括年度规划阶段、月度规划阶段和日前规划阶段，A：在年度规划阶段：采用地区典型负荷曲线数据，基于经济调度优化方法，建立求解该地区最大能够接纳的风电装机容量以及最小必须接入的风电装机容量的数学模型；B：在月度规划阶段：采用月度负荷预测数据，基于经济调度优化方法，建立求解该地区最大能够接纳的风电电量以及最小必须接入的风电电量的数学模型；C：在日前规划阶段：采用日前负荷预测数据、日前风电预测数据、月度火电机组启停优化结果，基于经济调度优化方法，建立求解该地区次日最大能够接纳的风电电量以及次日最小必须接入的风电电量的数学模型。The invention discloses a multi-stage wind power acceptance range calculation method based on unit combination and economic dispatch, including annual planning stage, monthly planning stage and day-ahead planning stage, A: In the annual planning stage: using regional typical load curve data, based on economic dispatch optimization method to establish a mathematical model for solving the maximum wind power installed capacity that can be accommodated in the region and the minimum wind power installed capacity that must be connected; The mathematical model of the amount of wind power that can be accepted and the minimum amount of wind power that must be connected; C: In the day-ahead planning stage: use the day-ahead load forecast data, day-ahead wind power forecast data, and monthly thermal power unit start-stop optimization results, based on the economic dispatch optimization method, to establish Solve the mathematical model of the maximum amount of wind power that can be received in the area on the next day and the minimum amount of wind power that must be connected in the next day.

Description

Calculation method of multi-stage wind power acceptance range based on unit combination and economic dispatch

技术领域technical field

本发明涉及一种基于机组组合和经济调度的多阶段风电接纳范围计算方法。The invention relates to a multi-stage wind power acceptance range calculation method based on unit combination and economic dispatch.

背景技术Background technique

随着能源与环境问题的日益严峻，环境友好型的可再生资源逐渐受到关注。风能以其运行成本低、污染少、储量丰富、开发潜力巨大等天然优势，近年来得到了大力发展，由于自然界中风固有的随机性和波动性，使得并网的风电具有以下不足：With the increasingly severe energy and environmental problems, environmentally friendly renewable resources have gradually attracted attention. With its natural advantages such as low operating cost, less pollution, abundant reserves, and huge development potential, wind energy has been vigorously developed in recent years. Due to the inherent randomness and volatility of wind in nature, grid-connected wind power has the following disadvantages:

1)强烈的不确定性：即风电出力受风速影响难以精确预测；1) Strong uncertainty: that is, wind power output is difficult to predict accurately due to wind speed;

2)反调峰特性：即风电机组出力与系统负荷具有负相关性，尤其是负荷水平较低而风电出力较高时会产生极大问题；2) Anti-peaking characteristics: that is, the wind turbine output has a negative correlation with the system load, especially when the load level is low and the wind power output is high, it will cause great problems;

3)风电机组的输出功率取决于风速的大小，目前风电场风速预测水平非常有限，其预测误差一般在25％～40％左右，因此大规模风电并网对电力系统的稳定运行带来了巨大的挑战。3) The output power of wind turbines depends on the wind speed. At present, the wind speed prediction level of wind farms is very limited, and its prediction error is generally about 25% to 40%. challenge.

现有能源结构的调节能力不足，为了平衡电网的波动性，系统的经济性变差，甚至在风电大发时，系统不得不弃风，随着风电在电网中渗透率的持续提高，电网弃风问题日益严重，因此，研究电网的风电接纳能力对于提高电网的安全性和经济性十分重要。The adjustment ability of the existing energy structure is insufficient. In order to balance the fluctuation of the power grid, the economy of the system becomes worse. The wind problem is becoming more and more serious. Therefore, it is very important to study the wind power acceptance capacity of the grid to improve the security and economy of the grid.

另外，关于电网接纳风电能力尚无明确的定义及标准的计算方法，对于不同领域，评估风电接纳能力所考虑的因素不同，决定了评估方法的不同。传统评估电网接纳能力通常是从风电规划角度，考虑运行方式、扰动方式等因素，用于决策最大风电装机容量，这类评估方法不适用于短期调度运行领域；还有一些研究方法包括基于调峰的分析等，通过建立代数模型对电网风电接纳能力进行分析，但难以保证评估结果的有效性和实用性。In addition, there is no clear definition and standard calculation method for the wind power acceptance capacity of the grid. For different fields, the factors considered in evaluating the wind power acceptance capacity are different, which determines the different evaluation methods. The traditional evaluation of grid capacity is usually from the perspective of wind power planning, considering factors such as operation mode and disturbance mode, and is used to determine the maximum installed capacity of wind power. This type of evaluation method is not suitable for short-term dispatching and operation; Through the establishment of an algebraic model to analyze the wind power acceptance capacity of the power grid, it is difficult to guarantee the validity and practicability of the evaluation results.

发明内容Contents of the invention

针对上述问题，本发明提供一种基于机组组合和经济调度的多阶段风电接纳范围计算方法，用于规划-月度-日前多阶段的风电接纳范围计算，该范围包括最大值和最小值，对分析电网的风电接纳能力具有重要意义，对电网的规划建设和经济调度具有参考价值。In view of the above problems, the present invention provides a multi-stage wind power acceptance range calculation method based on unit combination and economic dispatch, which is used for planning-monthly-day-ahead multi-stage wind power acceptance range calculation. The range includes maximum and minimum values, which are useful for analysis The wind power acceptance capacity of the power grid is of great significance and has reference value for the planning, construction and economic dispatch of the power grid.

术语定义：Definition of Terms:

(1)最大能够接纳的风电装机容量是指通过火电机组启停、机组出力等，该电网能够接纳的最大风电装机容量；(1) The maximum installed wind power capacity that can be accepted refers to the maximum installed wind power capacity that can be accommodated by the power grid through the start and stop of thermal power units, unit output, etc.;

(2)最小必须接入的风电装机容量是指当该地区火电机组装机不能满足负荷需求时，为了保证供需平衡，最小必须接入的风电装机容量；(2) The minimum installed capacity of wind power that must be connected refers to the minimum installed capacity of wind power that must be connected in order to ensure the balance between supply and demand when the thermal power assembly machines in the area cannot meet the load demand;

(3)月度最大能够接纳的风电电量是指通过火电机组启停、机组出力等，该电网在该月度能够接纳的最大风电电量；(3) The monthly maximum wind power capacity that can be accepted refers to the maximum wind power capacity that the power grid can accept in that month through the start and stop of thermal power units, unit output, etc.;

(4)月度最小必须接入的风电电量是指当该地区火电机组出力不能满足负荷需求时，为了保证供需平衡，该月度最小必须接入的风电电量；(4) The monthly minimum amount of wind power that must be connected refers to the monthly minimum amount of wind power that must be connected in order to ensure the balance between supply and demand when the output of thermal power units in the area cannot meet the load demand;

(5)次日最大能够接纳的风电电量是指通过火电机组启停、机组出力等，该电网次日能够接纳的最大风电电量；(5) The maximum amount of wind power that can be received on the next day refers to the maximum amount of wind power that can be received by the power grid on the next day through the start and stop of the thermal power unit, the output of the unit, etc.;

(6)次日最小必须接入的风电电量是指当该地区次日火电机组出力不能满足负荷需求时，为了保证供需平衡，次日最小必须接入的风电电量。(6) The minimum amount of wind power that must be connected in the next day refers to the minimum amount of wind power that must be connected in the next day in order to ensure the balance between supply and demand when the output of thermal power units in the area cannot meet the load demand the next day.

为实现上述技术目的，达到上述技术效果，本发明通过以下技术方案实现：In order to achieve the above-mentioned technical purpose and achieve the above-mentioned technical effect, the present invention is realized through the following technical solutions:

基于机组组合和经济调度的多阶段风电接纳范围计算方法，其特征在于，包括年度规划阶段、月度规划阶段和日前规划阶段，分别为：The multi-stage wind power acceptance range calculation method based on unit combination and economic dispatch is characterized in that it includes annual planning stage, monthly planning stage and day-ahead planning stage, which are respectively:

A：在年度规划阶段：A: During the annual planning phase:

A01：采用地区典型负荷曲线数据，基于经济调度优化方法，建立求解该地区最大能够接纳的风电装机容量以及最小必须接入的风电装机容量的数学模型；A01: Using the typical load curve data in the region, based on the economic dispatch optimization method, establish a mathematical model to solve the maximum wind power installed capacity that can be accepted in the region and the minimum wind power installed capacity that must be connected;

A02：求解模型，分别得到该地区最大能够接纳的风电装机容量以及最小必须接入的风电装机容量；A02: Solve the model to obtain the maximum installed capacity of wind power that can be accepted in the region and the minimum installed capacity of wind power that must be connected;

B：在月度规划阶段：B: During the monthly planning stage:

B01：采用月度负荷预测数据，基于经济调度优化方法，建立求解该地区月度最大能够接纳的风电电量以及月度最小必须接入的风电电量的数学模型；B01: Using monthly load forecast data and based on the economic scheduling optimization method, establish a mathematical model to solve the monthly maximum wind power capacity that can be accepted in the region and the monthly minimum wind power capacity that must be connected;

B02：求解模型，分别得到该地区月度最大能够接纳的风电电量以及月度最小必须接入的风电电量；B02: Solve the model to obtain the maximum monthly wind power capacity that can be accepted in the region and the monthly minimum wind power capacity that must be connected;

C：在日前规划阶段：C: In the day-ahead planning phase:

C01：采用日前负荷预测数据、日前风电预测数据、月度火电机组启停优化结果，基于经济调度优化方法，建立求解该地区次日最大能够接纳的风电电量以及次日最小必须接入的风电电量的数学模型；C01: Using the day-ahead load forecast data, day-ahead wind power forecast data, and monthly thermal power unit start-stop optimization results, based on the economic dispatch optimization method, establish a formula for solving the maximum wind power capacity that can be accepted in the next day and the minimum wind power power capacity that must be connected in the next day. mathematical model;

C02：求解模型，分别得到该地区次日最大能够接纳的风电电量以及次日最小必须接入的风电电量。C02: Solve the model, and obtain the maximum wind power capacity that can be received in the area on the next day and the minimum wind power capacity that must be connected in the next day.

优选，步骤A01中，最大能够接纳的风电装机容量的模型为：Preferably, in step A01, the model of the maximum admissible wind power installed capacity is:

目标函数：Objective function:

${maxP maxP}_{W W,, I I n no s the s t t a a l l l l}^{u u p p} - - - - - - ((11))$

约束条件：Restrictions:

${Σ Σ}_{i i = = 11}^{N N G G} {P P}_{G G i i,, t t} + + {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} = = {P P}_{L L,, t t} + + {P P}_{l l i i n no e e,, t t} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((22))$

${P P}_{G G i i}^{min min} {μ μ}_{i i,, t t} \leq \leq {P P}_{G G i i,, t t} \leq \leq {P P}_{G G i i}^{max max} {μ μ}_{i i,, t t},, ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((33))$

$\{\begin{matrix} {P P}_{G G i i,, t t} - - {P P}_{G G i i,, t t - - 11} \leq \leq {ΔP ΔP}_{G G i i}^{u u p p} \\ {P P}_{G G i i,, t t - - 11} - - {P P}_{G G i i,, t t} \leq \leq {ΔP ΔP}_{G G i i}^{d d n no} \end{matrix} ((i i = = 11,, ... ...,, N N G G;; t t = = 11,, ... ...,, N N T T)) - - - - - - ((44))$

$\{\begin{matrix} {Σ Σ}_{i i = = 11}^{N N G G} (({P P}_{G G i i,, t t} - - {P P}_{G G i i}^{min min})) {μ μ}_{i i,, t t} &GreaterEqual; &Greater Equal; {λ λ}_{11} {P P}_{L L,, t t} + + {λ λ}_{22} {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} \\ {Σ Σ}_{i i = = 11}^{N N G G} (({P P}_{G G i i}^{max max} - - {P P}_{Gi Gi,, t t})) {μ μ}_{i i,, t t} &GreaterEqual; &Greater Equal; {λ λ}_{11} {P P}_{L L,, t t} + + {λ λ}_{22} {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} \end{matrix} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((55))$

$\{\begin{matrix} (({X x}_{i i,, t t - - 11}^{o o n no} - - {T T}_{i i}^{o o n no})) \cdot &Center Dot; (({μ μ}_{i i,, t t - - 11} - - {μ μ}_{i i,, t t})) &GreaterEqual; &Greater Equal; 00 \\ (({X x}_{i i,, t t - - 11}^{o o f f f f} - - {T T}_{i i}^{o o f f f f})) \cdot &Center Dot; (({μ μ}_{i i,, t t} - - {μ μ}_{i i,, t t - - 11})) &GreaterEqual; &Greater Equal; 00 \end{matrix} ((i i = = 11,, ... ...,, N N G G;; t t = = 11,, ... ... N N T T)) - - - - - - ((66))$

${P P}_{f f l l o o w w,, l l,, t t} \leq \leq {P P}_{f f l l o o w w,, l l}^{max max} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((77))$

${P P}_{W W i i,, t t} \leq \leq {P P}_{W W,, I I n no s the s t t a a l l l l}^{u u p p} - - - - - - ((88))$

其中：式(1)为可接纳风电装机容量上限，NT为研究周期包含的时段数，NG为系统的火电机组台数，NW为系统的风电机组台数，P_Wi,t为时段t可接纳风电出力值，为全时段电网可接纳风电装机容量上限，P_Gi,t为火电机组i在时刻t的出力，P_Wj,t为风电场_j在时刻t的出力，P_L,t为系统在时刻t的负荷需求，P_line,t为在时刻t的联络线计划，为火电机组i的最小技术出力，为火电机组i的最大技术出力，μ_i,t为常规机组i在时刻t的启停机方式，μ_i,t＝1开机，μ_i,t＝0停机，为火电机组i的上升爬坡率限制，为火电机组i的下降爬坡率限制，λ₁、λ₂为系统的负荷备用系数，为机组i的最小开机时间，为机组i的最小停机时间，为机组i在时段t-1的连续开机时间，为机组i在时段t-1的连续停机时间，P_flow,l,t为传输线l在时刻t的直流潮流，为传输线l的直流潮流限制。Among them: Equation (1) is the upper limit of acceptable wind power installed capacity, NT is the number of periods included in the research period, NG is the number of thermal power units in the system, NW is the number of wind power units in the system, P _Wi,t is the acceptable wind power output in time period t value, is the upper limit of the installed capacity of wind power that can be accepted by the grid at all times, P _Gi,t is the output of thermal power unit i at time t, P _Wj,t is the output of wind farm _j at time t, P _L,t is the load of the system at time t Demand, P _line,t is the connection line plan at time t, Contribute to the minimum technology of thermal power unit i, is the maximum technical output of thermal power unit i, μ _i,t is the start-up and shutdown mode of conventional unit i at time t, μ _i,t = 1 to start, μ _i,t = 0 to stop, is the ramp rate limit of the thermal power unit i, is the limit of ramp rate of thermal power unit i, λ ₁ and λ ₂ are the load reserve coefficients of the system, is the minimum start-up time of unit i, is the minimum downtime of unit i, is the continuous start-up time of unit i in period t-1, is the continuous shutdown time of unit i in period t-1, P _flow,l,t is the DC flow of transmission line l at time t, is the DC power flow limit of the transmission line l.

最小必须接入的风电装机容量的模型为：The model of the minimum wind power installed capacity that must be connected is:

目标函数：Objective function:

${minP minP}_{W W,, I I n no s the s t t a a l l l l}^{d d n no} - - - - - - ((99))$

约束条件：Restrictions:

${P P}_{G G i i}^{min min} {μ μ}_{i i,, t t} \leq \leq {P P}_{G G i i,, t t} \leq \leq {P P}_{G G i i}^{max max} {μ μ}_{i i,, t t} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((33))$

$\{\begin{matrix} {Σ Σ}_{i i = = 11}^{N N G G} (({P P}_{G G i i,, t t} - - {P P}_{G G i i}^{min min})) {μ μ}_{i i,, t t} &GreaterEqual; &Greater Equal; {λ λ}_{11} {P P}_{L L,, t t} + + {λ λ}_{22} {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} \\ {Σ Σ}_{i i = = 11}^{N N G G} (({P P}_{G G i i}^{max max} - - {P P}_{G G i i,, t t})) {μ μ}_{i i,, t t} &GreaterEqual; &Greater Equal; {λ λ}_{11} {P P}_{L L,, t t} + + {λ λ}_{22} {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} \end{matrix} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((55))$

$\{\begin{matrix} (({X x}_{i i,, t t - - 11}^{o o n no} - - {T T}_{i i}^{o o n no})) \cdot \cdot (({μ μ}_{i i,, t t - - 11} - - {μ μ}_{i i,, t t})) &GreaterEqual; &Greater Equal; 00 \\ (({X x}_{i i,, t t - - 11}^{o o f f f f} - - {T T}_{i i}^{o o f f f f})) \cdot \cdot (({μ μ}_{i i,, t t} - - {μ μ}_{i i,, t t - - 11})) &GreaterEqual; &Greater Equal; 00 \end{matrix} ((i i = = 11,, ... ...,, N N G G;; t t = = 11,, ... ... N N T T)) - - - - - - ((66))$

${P P}_{W W i i,, t t} \leq \leq {P P}_{W W,, I I n no s the s t t a a l l l l}^{d d n no} - - - - - - ((1010))$

其中，式(9)为可接入风电装机容量下限，为风电装机容量下限。Among them, formula (9) is the lower limit of the installed capacity of wind power that can be connected, is the lower limit of wind power installed capacity.

本发明的有益效果是：The beneficial effects of the present invention are:

本方法以机组组合和经济调度为基础，在规划-月度-日前多阶段分别建立了以风电装机容量、风电电量最大为目标的优化模型，采用混合整数规划算法求解，有助于提高电网调度的精细化水平。可以在风电规划建设时提供参考，可以提前防范大规模风电并网对电网安全运行带来的风险，可以为调度人员进行日前发电计划制定提供指导，提高了大规模风电介入后电网的安全性和经济性。Based on unit combination and economic dispatch, this method establishes an optimization model aiming at the maximum wind power installed capacity and wind power output in the stages of planning-monthly-day-ahead, and adopts the mixed integer programming algorithm to solve the problem, which is helpful to improve the efficiency of power grid dispatching. level of refinement. It can provide reference in the planning and construction of wind power, prevent the risks brought by large-scale wind power grid-connection to the safe operation of the power grid in advance, and provide guidance for dispatchers to make day-ahead power generation plans, improving the safety and security of the power grid after large-scale wind power intervention. economy.

具体实施方式Detailed ways

下面结合具体的实施例对本发明技术方案作进一步的详细描述，以使本领域的技术人员可以更好的理解本发明并能予以实施，但所举实施例不作为对本发明的限定。The technical solutions of the present invention will be further described in detail below in conjunction with specific examples, so that those skilled in the art can better understand the present invention and implement it, but the examples given are not intended to limit the present invention.

基于机组组合和经济调度的多阶段风电接纳范围计算方法，包括年度规划阶段、月度规划阶段和日前规划阶段，本发明的电网风电接纳能力优化评估方法提出研究风电接纳区间，主要基于机组组合和经济调度技术，以风电出力之和最大或最小为优化目标，建立优化求解模型，可通过商业软件(ILOG/CPLEX)对模型进行求解，获得电网的风电接纳能力结果。A multi-stage wind power acceptance range calculation method based on unit combination and economic dispatch, including annual planning stage, monthly planning stage and day-ahead planning stage. The wind power acceptance capacity optimization evaluation method of the power grid of the present invention proposes to study the wind power acceptance interval, mainly based on unit combination and economic Scheduling technology takes the maximum or minimum sum of wind power output as the optimization goal, establishes an optimization solution model, and can solve the model through commercial software (ILOG/CPLEX) to obtain the wind power acceptance capacity of the grid.

下面对各个规划阶段进行详细说明：The individual planning stages are described in detail below:

A：在年度规划阶段：A: During the annual planning phase:

A02：求解模型，分别得到该地区最大能够接纳的风电装机容量以及最小必须接入的风电装机容量。A02: Solve the model to obtain the maximum installed capacity of wind power that can be accepted in the region and the minimum installed capacity of wind power that must be connected.

年度最大能够接纳风电装机容量模型：该模型以火电机组的机组启停、火电机组出力功率、风电场的出力功率为决策变量；以最大化风电装机容量为优化目标；约束条件包括有功功率平衡约束、火电机组出力上下限约束、火电机组爬坡约束、备用需求约束、火电机组最小启停时间约束、线路潮流约束、风电功率范围约束。优选，优化模型如下：Annual maximum acceptable wind power installed capacity model: This model takes the start and stop of thermal power units, the output power of thermal power units, and the output power of wind farms as decision variables; the optimization goal is to maximize the installed capacity of wind power; constraints include active power balance constraints , thermal power unit output upper and lower limit constraints, thermal power unit climbing constraints, backup demand constraints, thermal power unit minimum start-stop time constraints, line flow constraints, wind power range constraints. Preferably, the optimization model is as follows:

目标函数：Objective function:

约束条件：Restrictions:

需说明的是：上述优化模型中，待优化决策变量为机组启停状态和出力大小，其中既有连续型变量又有离散型变量，约束条件中既有功率平衡约束、机组爬坡约束等线性约束条件，又包含常规机组最小启停时间约束等非线性约束条件，同时还要考虑常规机组最小启停时间和爬坡约束等在不同时间段的耦合性，因此优化模型属于多变量、非线性混合整数规划问题。由于ILOG/CPLEX只能求解线性的混合整数规划模型，因此若调用ILOG/CPLEX求解本文所建立的优化模型之前，必须将模型中的非线性因素进行线性化，结果如下：It should be noted that in the above optimization model, the decision-making variables to be optimized are the start-stop status and output size of the unit, among which there are both continuous variables and discrete variables, and the constraints include power balance constraints, unit ramp constraints, etc. Constraint conditions, including non-linear constraints such as the minimum start-stop time constraints of conventional units, and the coupling between the minimum start-stop time and ramp constraints of conventional units in different time periods, so the optimization model is a multivariate, nonlinear Mixed integer programming problems. Since ILOG/CPLEX can only solve linear mixed integer programming models, before calling ILOG/CPLEX to solve the optimization model established in this paper, the nonlinear factors in the model must be linearized, and the results are as follows:

$\{\begin{matrix} {μ μ}_{i i,, t t} = = 11,, & 11 \leq \leq t t \leq \leq {U u}_{i i} \\ {Σ Σ}_{τ τ = = t t}^{t t + + {T T}_{i i}^{o o n no} - - 11} {μ μ}_{i i,, τ τ} &GreaterEqual; &Greater Equal; {T T}_{i i}^{o o n no} (({μ μ}_{i i,, t t} - - {μ μ}_{i i,, t t - - 11})) & {U u}_{i i} \leq \leq t t \leq \leq N N T T - - {T T}_{i i}^{o o n no} + + 11 \\ {Σ Σ}_{τ τ = = t t}^{N N T T} [[{μ μ}_{i i,, τ τ} - - (({μ μ}_{i i,, t t} - - {μ μ}_{i i,, t t - - 11}))]] &GreaterEqual; &Greater Equal; 00 & N N T T - - {T T}_{i i}^{o o n no} + + 22 \leq \leq t t \leq \leq N N T T \end{matrix} - - - - - - ((I I))$

$\{\begin{matrix} {μ μ}_{i i,, t t} = = 00,, & 11 \leq \leq t t \leq \leq {D D.}_{i i} \\ {Σ Σ}_{τ τ = = t t}^{t t + + {T T}_{i i}^{o o f f f f} - - 11} ((11 - - {μ μ}_{i i,, τ τ})) &GreaterEqual; &Greater Equal; {T T}_{i i}^{off off} (({μ μ}_{i i,, t t - - 11} - - {μ μ}_{i i,, t t})) & {D D.}_{i i} \leq \leq t t \leq \leq N N T T - - {T T}_{i i}^{o o f f f f} + + 11 \\ {Σ Σ}_{τ τ = = t t}^{N N T T} [[11 - - {μ μ}_{i i,, τ τ} - - (({μ μ}_{i i,, t t - - 11} - - {μ μ}_{i i,, t t}))]] &GreaterEqual; &Greater Equal; 00 & N N T T - - {T T}_{i i}^{o o f f f f} + + 22 \leq \leq t t \leq \leq N N T T \end{matrix} - - - - - - ((I I I I))$

式(Ⅰ)为最小开机运行时间约束的线性化，式(Ⅱ)为最小停机时间约束的线性化。Equation (I) is the linearization of the minimum uptime constraint, and Equation (II) is the linearization of the minimum downtime constraint.

可采用ILOG/CPLEX求解上述线性化模型，得到最大接纳风电装机容量的最优解为 ILOG/CPLEX can be used to solve the above linearized model to obtain the maximum installed capacity of wind power The optimal solution for

最小必须接入风电装机容量模型：该模型以火电机组的机组启停、火电机组出力功率、风电场的出力功率为决策变量；以最小化风电装机容量为优化目标，优选，优化模型如下：The minimum wind power installed capacity model that must be connected: This model takes the start and stop of thermal power units, the output power of thermal power units, and the output power of wind farms as decision variables; the optimization goal is to minimize the installed capacity of wind power. The optimization model is as follows:

目标函数：Objective function:

约束条件：Restrictions:

可采用ILOG/CPLEX求解上述线性化模型，得到最小接入风电装机容量的最优解为 ILOG/CPLEX can be used to solve the above linearized model to obtain the minimum installed capacity of connected wind power The optimal solution for

B：在月度规划阶段：B: During the monthly planning stage:

B01：采用月度负荷预测数据，基于经济调度优化方法，建立求解该地区月度最大能够接纳的风电电量以及月度最小必须接入的风电电量的数学模型；B01: Using monthly load forecast data and based on the economic dispatch optimization method, establish a mathematical model to solve the monthly maximum wind power capacity that can be accepted in the region and the monthly minimum wind power power capacity that must be connected;

B02：求解模型，分别得到该地区月度最大能够接纳的风电电量以及月度最小必须接入的风电电量。B02: Solve the model to obtain the maximum monthly wind power capacity that can be received in the region and the monthly minimum wind power capacity that must be connected.

月度最大能够接纳风电电量模型：该模型以火电机组的机组启停、火电机组出力功率、风电场的出力功率为决策变量，以在该月度能够接纳的最大风电电量为优化目标；约束条件包括有功功率平衡约束，火电机组出力上下限约束，火电机组爬坡约束，备用需求约束，火电机组启停时间约束，线路潮流约束，风电装机容量约束。优选，优化模型如下：Monthly maximum acceptable wind power model: This model takes the start and stop of thermal power units, the output power of thermal power units, and the output power of wind farms as decision variables, and takes the maximum wind power capacity that can be accepted in this month as the optimization goal; the constraints include active power Power balance constraints, thermal power unit output upper and lower limit constraints, thermal power unit climbing constraints, backup demand constraints, thermal power unit start-stop time constraints, line flow constraints, wind power installed capacity constraints. Preferably, the optimization model is as follows:

目标函数：Objective function:

${maxP maxP}_{W W m m,, S S u u m m}^{u u p p} = = {Σ Σ}_{i i = = 11}^{N N W W} {Σ Σ}_{t t = = 11}^{N N T T} {P P}_{W W i i,, t t} - - - - - - ((1111))$

约束条件：Restrictions:

${P P}_{W W i i,, t t} \leq \leq {\overset{&OverBar; &OverBar;}{P P}}_{W W i i,, i i n no s the s t t a a l l l l} - - - - - - ((1212))$

其中，式(11)为月度风电接纳上限，为月度风电接纳容量之和上限，为已知的电网风电场i的装机容量；NT为研究周期包含的时段数，一般而言，在月度规划阶段，每4h取一时段，一天6时段，并乘以当月天数即是NT的数值。Among them, formula (11) is the upper limit of monthly wind power acceptance, is the upper limit of the sum of the monthly wind power receiving capacity, is the known installed capacity of grid wind farm i; NT is the number of periods included in the research cycle. Generally speaking, in the monthly planning stage, a period is taken every 4 hours, 6 periods a day, and multiplied by the number of days in the current month is the value of NT .

可采用ILOG/CPLEX求解上述线性化模型，得到月度风电最大接纳的最优解为 ILOG/CPLEX can be used to solve the above linearized model to obtain the maximum acceptance of monthly wind power The optimal solution for

月度最小必须接入风电电量模型：该模型以火电机组的机组启停、火电机组出力功率、风电场的出力功率为决策变量；以该月度最小必须接纳的风电电量为优化目标，优选，优化模型如下：Monthly minimum must-connect wind power model: This model takes the start-up and stop of thermal power units, the output power of thermal power units, and the output power of wind farms as decision variables; takes the monthly minimum wind power that must be accepted as the optimization goal, optimizes, and optimizes the model as follows:

目标函数：Objective function:

${minP minP}_{W W m m,, S S u u m m}^{d d n no} = = {Σ Σ}_{i i = = 11}^{N N W W} {Σ Σ}_{t t = = 11}^{N N T T} {P P}_{W W i i,, t t} - - - - - - ((1313))$

约束条件：Restrictions:

${P P}_{W W i i,, t t} \leq \leq {\overset{&OverBar; &OverBar;}{P P}}_{W W i i,, i i n no s the s t t a a l l l l} ((1212)) . .$

式(13)为月度风电接入下限，为月度风电接纳容量之和下限。Equation (13) is the lower limit of monthly wind power access, It is the lower limit of the sum of monthly wind power receiving capacity.

可采用ILOG/CPLEX求解上述线性化模型，得到月度风电最小必须接纳电量的最优解为 ILOG/CPLEX can be used to solve the above linearized model, and the minimum monthly wind power must be accepted The optimal solution for

C：在日前规划阶段：C: In the day-ahead planning phase:

次日最大能够接纳风电电量模型：该模型以火电机组的机组启停、火电机组出力功率、风电场的出力功率为决策变量，以次日能够接纳的最大风电电量为优化目标；约束条件包括有功功率平衡约束、火电机组出力上下限约束、火电机组爬坡约束、备用需求约束、火电机组最小启停时间约束、线路潮流约束、风电日前预测出力约束。优选，优化模型如下：The maximum acceptable wind power capacity model for the next day: This model takes the start and stop of the thermal power unit, the output power of the thermal power unit, and the output power of the wind farm as the decision variables, and takes the maximum wind power capacity that can be accepted the next day as the optimization goal; the constraints include active power Power balance constraints, thermal power unit output upper and lower limit constraints, thermal power unit climbing constraints, backup demand constraints, thermal power unit minimum start-stop time constraints, line flow constraints, and wind power day-ahead forecast output constraints. Preferably, the optimization model is as follows:

目标函数：Objective function:

${maxP maxP}_{W W d d,, S S u u m m}^{u u p p} = = {Σ Σ}_{i i = = 11}^{N N W W} {Σ Σ}_{t t = = 11}^{N N T T} {P P}_{W W i i,, t t} - - - - - - ((1414))$

约束条件：Restrictions:

P_Wi,t≤P_Wi,forecast(15)P _Wi,t ≤P _Wi,forecast (15)

其中，式(14)为次日风电接纳上限，为次日风电接纳容量之和上限，P_Wi,forecast为风电场i预测出力；NT为研究周期包含的时段数，一般而言，在日前规划阶段，每15min取一时段，一天共96时段，即NT的数值为96。Among them, formula (14) is the upper limit of wind power acceptance for the next day, is the upper limit of the sum of wind power receiving capacity for the next day, P _Wi,forecast is the predicted output of wind farm i; NT is the number of time periods included in the research cycle, generally speaking, in the day-ahead planning stage, a time period is taken every 15 minutes, a total of 96 time periods a day, That is, the value of NT is 96.

可采用ILOG/CPLEX求解上述线性化模型，得到次日风电最大接纳的最优解为 ILOG/CPLEX can be used to solve the above linearized model to obtain the maximum wind power acceptance for the next day The optimal solution for

最小必须接入风电电量模型：该模型以火电机组的机组启停、火电机组出力功率、风电场的出力功率为决策变量；以次日最小必须接纳的风电电量为优化目标，优选，优化模型如下：The minimum wind power quantity that must be connected to the model: the model takes the start and stop of the thermal power unit, the output power of the thermal power unit, and the output power of the wind farm as the decision variables; the minimum wind power that must be accepted in the next day is the optimization goal, and the optimization model is as follows :

目标函数：Objective function:

${minP minP}_{W W d d,, S S u u m m}^{d d n no} = = {Σ Σ}_{i i = = 11}^{N N W W} {Σ Σ}_{t t = = 11}^{N N T T} {P P}_{W W i i,, t t} - - - - - - ((1616))$

约束条件：Restrictions:

$\{\begin{matrix} (({X x}_{i i,, t t - - 11}^{o o n no} - - {T T}_{i i}^{o o n no})) \cdot \cdot (({μ μ}_{i i,, t t - - 11} - - {μ μ}_{i i,, t t})) &GreaterEqual; &Greater Equal; 00 \\ (({X x}_{i i,, t t - - 11}^{o o f f f f} - - {T T}_{i i}^{o o f f f f})) \cdot &Center Dot; (({μ μ}_{i i,, t t} - - {μ μ}_{i i,, t t - - 11})) &GreaterEqual; &Greater Equal; 00 \end{matrix} ((i i = = 11,, ... ...,, N N G G;; t t = = 11,, ... ... N N T T)) - - - - - - ((66))$

${P P}_{f f l l o o w w,, l l,, t t} \leq \leq {P P}_{f f l l o o w w . . l l}^{max max} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((77))$

P_Wi,t≤P_Wi,forecast(15)P _Wi,t ≤P _Wi,forecast (15)

其中，式(16)为次日风电接纳下限，为次日风电接纳容量之和下限。Among them, formula (16) is the lower limit of wind power acceptance for the next day, It is the lower limit of the sum of wind power receiving capacity for the next day.

可采用ILOG/CPLEX求解上述线性化模型，得到次日风电最小接入的最优解为 ILOG/CPLEX can be used to solve the above linearized model to obtain the minimum wind power access for the next day The optimal solution for

本发明能够在规划-月度-日前多阶段电网运行下，动态分析电网的风电接纳能力，并且分析电网中各种因素对风电接纳能力的影响。以改变单一因素为原则，分别分析负荷、火电机组最小出力、火电机组最大出力、火电机组最小启停时间、火电机组爬坡约束、电网线路传输能力对电网风电接纳能力上下限的影响，最后得出的结论是：适当增加电网负荷、减小火电机组最小技术出力和最小启停时间、增大火电机组最大技术出力和爬坡约束限制、改善电网线路传输能力等，有利于电网接纳更多的风电。The invention can dynamically analyze the wind power acceptance capacity of the power grid under the planning-monthly-day-ahead multi-stage power grid operation, and analyze the influence of various factors in the power grid on the wind power acceptance capacity. Based on the principle of changing a single factor, the effects of load, minimum output of thermal power units, maximum output of thermal power units, minimum start-stop time of thermal power units, climbing constraints of thermal power units, and transmission capacity of power grid lines on the upper and lower limits of wind power acceptance capacity of the power grid are analyzed respectively. The conclusions drawn are: appropriately increasing the grid load, reducing the minimum technical output of thermal power units and the minimum start-stop time, increasing the maximum technical output of thermal power units and the limit of climbing constraints, improving the transmission capacity of power grid lines, etc., are conducive to the grid to accept more wind power.

本发明在实际电网数据下开展，以机组组合和经济调度为基础，充分考虑影响电网接纳风电的各种因素，建立了以风电出力最大为目标的优化模型，调用商业软件(ILOG/CPLEX)中的混合整数规划算法求解，优化评估规划-月度-日前多时间段的风电接纳能力，有助于在电网规划建设时为风电接入提供参考，提前防范大规模风电并网对电网安全运行带来的风险，为调度人员进行发电计划制定提供指导，改善了大规模风电介入后电网的安全性和经济性。The present invention is carried out under the actual power grid data, based on unit combination and economic scheduling, fully considers various factors that affect the acceptance of wind power by the grid, establishes an optimization model with the goal of maximizing wind power output, and calls commercial software (ILOG/CPLEX) The mixed integer programming algorithm solves the problem, optimizes and evaluates the wind power acceptance capacity of planning-monthly-ahead of the multi-time period, which helps to provide reference for wind power connection during grid planning and construction, and prevents large-scale wind power grid-connected from bringing harm to the safe operation of the grid in advance. It provides guidance for dispatchers to formulate power generation plans, and improves the security and economy of the power grid after large-scale wind power intervention.

以上仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书内容所作的等效结构或者等效流程变换，或者直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。The above are only preferred embodiments of the present invention, and are not intended to limit the patent scope of the present invention. All equivalent structures or equivalent process transformations made using the content of the description of the present invention, or directly or indirectly used in other related technical fields, are the same as The theory is included in the patent protection scope of the present invention.

Claims

1. A multi-stage wind power acceptance range calculation method based on unit combination and economic dispatch, characterized in that it includes the annual planning stage, monthly planning stage and day-ahead planning stage, respectively:

A: During the annual planning phase:

A01: Using the typical load curve data in the region, based on the economic dispatch optimization method, establish a mathematical model to solve the maximum wind power installed capacity that can be accepted in the region and the minimum wind power installed capacity that must be connected;

A02: Solve the model to obtain the maximum installed capacity of wind power that can be accommodated in the region and the minimum installed capacity of wind power that must be connected;

B: During the monthly planning stage:

B01: Using monthly load forecast data and based on the economic scheduling optimization method, establish a mathematical model to solve the monthly maximum wind power capacity that can be accepted in the region and the monthly minimum wind power capacity that must be connected;

B02: Solve the model to obtain the maximum monthly wind power capacity that can be accepted in the region and the monthly minimum wind power capacity that must be connected;

C: In the day-ahead planning phase:

C01: Using the day-ahead load forecast data, day-ahead wind power forecast data, and monthly thermal power unit start-stop optimization results, based on the economic dispatch optimization method, establish a formula for solving the maximum wind power capacity that can be accepted in the next day and the minimum wind power power capacity that must be connected in the next day. mathematical model;

C02: Solve the model, and obtain the maximum wind power capacity that can be accepted in the next day and the minimum wind power capacity that must be connected in the next day.

2. The multi-stage wind power acceptance range calculation method based on unit combination and economic dispatch according to claim 1, characterized in that, in step A01, the model of the maximum wind power installed capacity that can be accepted is:

Objective function:

{maxP maxP}_{W W,, I I n no s the s t t a a l l l l}^{u u p p} - - - - - - ((11))

Restrictions:

{Σ Σ}_{i i = = 11}^{N N G G} {P P}_{G G i i,, t t} + + {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} = = {P P}_{L L,, t t} + + {P P}_{l l i i n no e e,, t t} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((22))

{P P}_{G G i i}^{min min} {μ μ}_{i i,, t t} \leq \leq {P P}_{G G i i,, t t} \leq \leq {P P}_{G G i i}^{max max} {μ μ}_{i i,, t t} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((33))

\{\begin{matrix} {P P}_{G G i i,, t t} - - {P P}_{G G i i,, t t - - 11} \leq \leq {ΔP ΔP}_{G G i i}^{u u p p} \\ {P P}_{G G i i,, t t - - 11} - - {P P}_{G G i i,, t t} \leq \leq {ΔP ΔP}_{G G i i}^{d d n no} \end{matrix} ((i i = = 11,, ... ...,, N N G G;; t t = = 11,, ... ...,, N N T T)) - - - - - - ((44))

\{\begin{matrix} {Σ Σ}_{i i = = 11}^{N N G G} (({P P}_{G G i i,, t t} - - {P P}_{G G i i}^{min min})) {μ μ}_{i i,, t t} &GreaterEqual; &Greater Equal; {λ λ}_{11} {P P}_{L L,, t t} + + {λ λ}_{22} {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} \\ {Σ Σ}_{i i = = 11}^{N N G G} (({P P}_{C C i i}^{max max} - - {P P}_{G G i i,, t t})) {μ μ}_{i i,, t t} &GreaterEqual; &Greater Equal; {λ λ}_{11} {P P}_{L L,, t t} + + {λ λ}_{22} {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} \end{matrix} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((55))

\{\begin{matrix} (({X x}_{i i,, t t - - 11}^{o o n no} - - {T T}_{i i}^{o o n no})) \cdot \cdot (({μ μ}_{i i,, t t - - 11} - - {μ μ}_{i i,, t t})) &GreaterEqual; &Greater Equal; 00 \\ (({X x}_{i i,, t t - - 11}^{o o f f f f} - - {T T}_{i i}^{o o f f f f})) \cdot \cdot (({μ μ}_{i i,, t t} - - {μ μ}_{i i,, t t - - 11})) &GreaterEqual; &Greater Equal; 00 \end{matrix} ((i i = = 11,, ... ...,, N N G G;; t t = = 11,, ... ...,, N N T T)) - - - - - - ((66))

{P P}_{f f l l o o w w,, l l,, t t} \leq \leq {P P}_{f f l l o o w w,, l l}^{max max} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((77))

{P P}_{W W i i,, t t} \leq \leq {P P}_{W W,, I I n no s the s t t a a l l l l}^{u u p p} - - - - - - ((88))

Among them: Equation (1) is the upper limit of acceptable wind power installed capacity, NT is the number of periods included in the research period, NG is the number of thermal power units in the system, NW is the number of wind power units in the system, P _Wi,t is the acceptable wind power output in time period t value, is the upper limit of the installed capacity of wind power that can be accepted by the grid at all times, P _Gi,t is the output of thermal power unit i at time t, P _Wj,t is the output of wind farm _j at time t, P _L,t is the load of the system at time t Demand, P _line,t is the connection line plan at time t, Contribute to the minimum technology of thermal power unit i, is the maximum technical output of thermal power unit i, μ _i,t is the start-up and shutdown mode of conventional unit i at time t, μ _i,t = 1 to start, μ _i,t = 0 to stop, is the ramp rate limit of the thermal power unit i, is the limit of ramp rate of thermal power unit i, λ ₁ and λ ₂ are the load reserve coefficients of the system, is the minimum start-up time of unit i, is the minimum downtime of unit i, is the continuous start-up time of unit i in period t-1, is the continuous shutdown time of unit i in period t-1, P _flow,l,t is the DC flow of transmission line l at time t, is the DC power flow limit of the transmission line l.

3. The multi-stage wind power acceptance range calculation method based on unit combination and economic dispatch according to claim 2, characterized in that, in step A01, the model of the minimum wind power installed capacity that must be connected is:

Objective function:

{minP minP}_{W W,, I I n no s the s t t a a l l l l}^{d d n no} - - - - - - ((99))

Restrictions:

{Σ Σ}_{i i = = 11}^{N N G G} {P P}_{G G i i,, t t} + + {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} = = {P P}_{L L,, t t} + + {P P}_{l l i i n no e e,, t t} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((22))

{P P}_{G G i i}^{min min} \leq \leq {μ μ}_{i i,, t t} \leq \leq {P P}_{G G i i,, t t} \leq \leq {P P}_{G G i i}^{max max} {μ μ}_{i i,, t t} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((33))

\{\begin{matrix} {P P}_{G G i i,, t t} - - {P P}_{G G i i,, t t - - 11} \leq \leq {ΔP ΔP}_{G G i i}^{u u p p} \\ {P P}_{G G i i,, t t - - 11} - - {P P}_{G G i i,, t t} \leq \leq {ΔP ΔP}_{G G i i}^{d d n no} \end{matrix} ((i i = = 11,, ... ...,, N N G G;; t t = = 11,, ... ...,, N N T T)) - - - - - - ((44))

\{\begin{matrix} {Σ Σ}_{i i = = 11}^{N N G G} (({P P}_{G G i i,, t t} - - {P P}_{G G i i}^{min min})) {μ μ}_{i i,, t t} &GreaterEqual; &Greater Equal; {λ λ}_{11} {P P}_{L L,, t t} + + {λ λ}_{22} {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} \\ {Σ Σ}_{i i = = 11}^{N N G G} (({P P}_{G G i i}^{max max} - - {P P}_{G G i i,, t t})) {μ μ}_{i i,, t t} &GreaterEqual; &Greater Equal; {λ λ}_{11} {P P}_{L L,, t t} + + {λ λ}_{22} {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} \end{matrix} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((55))

\{\begin{matrix} (({X x}_{i i,, t t - - 11}^{o o n no} - - {T T}_{i i}^{o o n no})) \cdot \cdot (({μ μ}_{i i,, t t - - 11} - - {μ μ}_{i i,, t t})) &GreaterEqual; &Greater Equal; 00 \\ (({X x}_{i i,, t t - - 11}^{o o f f f f} - - {T T}_{i i}^{o o f f f f})) \cdot \cdot (({μ μ}_{i i,, t t} - - {μ μ}_{i i,, t t - - 11})) &GreaterEqual; &Greater Equal; 00 \end{matrix} ((i i = = 11,, ... ...,, N N G G;; t t = = 11,, ... ...,, N N T T)) - - - - - - ((66))

{P P}_{f f l l o o w w,, l l,, t t} \leq \leq {P P}_{f f l l o o w w,, l l}^{max max} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((77))

{P P}_{W W i i,, t t} \leq \leq {P P}_{W W,, I I n no s the s t t a a l l l l}^{d d n no} - - - - - - ((1010))

Among them, formula (9) is the lower limit of the installed capacity of wind power that can be connected, is the lower limit of wind power installed capacity.

4. The multi-stage wind power acceptance range calculation method based on unit combination and economic dispatch according to claim 1, characterized in that, in step B01, the monthly maximum wind power capacity model that can be accepted is:

Objective function:

{maxP maxP}_{W W m m,, S S u u m m}^{u u p p} = = {Σ Σ}_{i i = = 11}^{N N W W} {Σ Σ}_{t t = = 11}^{N N T T} {P P}_{W W i i,, t t} - - - - - - ((1111))

Restrictions:

{Σ Σ}_{i i = = 11}^{N N G G} {P P}_{G G i i,, t t} + + {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} = = {P P}_{L L,, t t} + + {P P}_{l l i i n no e e,, t t} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((22))

{P P}_{G G i i}^{min min} {μ μ}_{i i,, t t} \leq \leq {P P}_{G G i i,, t t} \leq \leq {P P}_{G G i i}^{max max} {μ μ}_{i i,, t t} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((33))

\{\begin{matrix} {P P}_{G G i i,, t t} - - {P P}_{G G i i,, t t - - 11} \leq \leq {ΔP ΔP}_{G G i i}^{u u p p} \\ {P P}_{G G i i,, t t - - 11} - - {P P}_{G G i i,, t t} \leq \leq {ΔP ΔP}_{G G i i}^{d d n no} \end{matrix} ((i i = = 11,, ... ...,, N N G G;; t t = = 11,, ... ...,, N N T T)) - - - - - - ((44))

\{\begin{matrix} {Σ Σ}_{i i = = 11}^{N N G G} (({P P}_{G G i i,, t t} - - {P P}_{G G i i}^{min min})) {μ μ}_{i i,, t t} &GreaterEqual; &Greater Equal; {λ λ}_{11} {P P}_{L L,, t t} + + {λ λ}_{22} {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} \\ {Σ Σ}_{i i = = 11}^{N N G G} (({P P}_{G G i i}^{max max} - - {P P}_{G G i i,, t t})) {μ μ}_{i i,, t t} &GreaterEqual; &Greater Equal; {λ λ}_{11} {P P}_{L L,, t t} + + {λ λ}_{22} {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} \end{matrix} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((55))

\{\begin{matrix} (({X x}_{i i,, t t - - 11}^{o o n no} - - {T T}_{i i}^{o o n no})) \cdot &Center Dot; (({μ μ}_{i i,, t t - - 11} - - {μ μ}_{i i,, t t})) &GreaterEqual; &Greater Equal; 00 \\ (({X x}_{i i,, t t - - 11}^{o o f f f f} - - {T T}_{i i}^{o o f f f f})) \cdot \cdot (({μ μ}_{i i,, t t} - - {μ μ}_{i i,, t t - - 11})) &GreaterEqual; &Greater Equal; 00 \end{matrix} ((i i = = 11,, ... ...,, N N G G;; t t = = 11,, ... ...,, N N T T)) - - - - - - ((66))

{P P}_{f f l l o o w w,, l l,, t t} \leq \leq {P P}_{f f l l o o w w,, l l}^{max max} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((77))

{P P}_{W W i i,, t t} \leq \leq {\overset{&OverBar; &OverBar;}{P P}}_{W W i i,, i i n no s the s t t a a l l l l} - - - - - - ((1212))

Among them, formula (11) is the upper limit of monthly wind power acceptance, is the upper limit of the sum of the monthly wind power receiving capacity, is the known installed capacity of grid wind farm i; NT is the number of periods included in the research period.

5. The multi-stage wind power acceptance range calculation method based on unit combination and economic dispatch according to claim 4, characterized in that, in step B01, the monthly minimum wind power capacity model that must be connected is:

Objective function:

{minP minP}_{W W m m,, S S u u m m}^{d d n no} = = {Σ Σ}_{i i = = 11}^{N N W W} {Σ Σ}_{t t = = 11}^{N N T T} {P P}_{W W i i,, t t} - - - - - - ((1313))

Restrictions:

{Σ Σ}_{i i = = 11}^{N N G G} {P P}_{G G i i,, t t} + + {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} = = {P P}_{L L,, t t} + + {P P}_{l l i i n no e e,, t t} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((22))

{P P}_{G G i i}^{min min} {μ μ}_{i i,, t t} \leq \leq {P P}_{G G i i,, t t} \leq \leq {P P}_{G G i i}^{max max} {μ μ}_{i i,, t t} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((33))

\{\begin{matrix} {P P}_{G G i i,, t t} - - {P P}_{G G i i,, t t - - 11} \leq \leq {ΔP ΔP}_{G G i i}^{u u p p} \\ {P P}_{G G i i,, t t - - 11} - - {P P}_{G G i i,, t t} \leq \leq {ΔP ΔP}_{G G i i}^{u u p p} \end{matrix} ((i i = = 11,, ... ...,, N N G G;; t t = = 11,, ... ...,, N N T T)) - - - - - - ((44))

\{\begin{matrix} {Σ Σ}_{i i = = 11}^{N N G G} (({P P}_{G G i i,, t t} - - {P P}_{G G i i}^{min min})) {μ μ}_{i i,, t t} &GreaterEqual; &Greater Equal; {λ λ}_{11} {P P}_{L L,, t t} + + {λ λ}_{22} {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} \\ {Σ Σ}_{i i = = 11}^{N N G G} (({P P}_{G G i i}^{max max} - - {P P}_{G G i i,, t t})) {μ μ}_{i i,, t t} &GreaterEqual; &Greater Equal; {λ λ}_{11} {P P}_{L L,, t t} + + {λ λ}_{22} {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} \end{matrix} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((55))

\{\begin{matrix} (({X x}_{i i,, t t - - 11}^{o o n no} - - {T T}_{i i}^{o o n no})) \cdot &Center Dot; (({μ μ}_{i i,, t t - - 11} - - {μ μ}_{i i,, t t})) &GreaterEqual; &Greater Equal; 00 \\ (({X x}_{i i,, t t - - 11}^{o o f f f f} - - {T T}_{i i}^{o o f f f f})) \cdot &Center Dot; (({μ μ}_{i i,, t t} - - {μ μ}_{i i,, t t - - 11})) &GreaterEqual; &Greater Equal; 00 \end{matrix} ((i i = = 11,, ... ...,, N N G G;; t t = = 11,, ... ...,, N N T T)) - - - - - - ((66))

{P P}_{f f l l o o w w,, l l,, t t} \leq \leq {P P}_{f f l l o o w w,, l l}^{max max} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((77))

{P P}_{W W i i,, t t} \leq \leq {\overset{&OverBar; &OverBar;}{P P}}_{W W i i,, i i n no s the s t t a a l l l l} - - - - - - ((1212)) . .

6. The multi-stage wind power acceptance range calculation method based on unit combination and economic dispatch according to claim 1, characterized in that, in step C01, the model of the maximum wind power capacity that can be accepted the next day is: objective function:

{maxP maxP}_{W W d d,, S S u u m m}^{u u p p} = = {Σ Σ}_{i i = = 11}^{N N W W} {Σ Σ}_{t t = = 11}^{N N T T} {P P}_{W W i i,, t t} - - - - - - ((1414))

Restrictions:

{Σ Σ}_{i i = = 11}^{N N G G} {P P}_{G G i i,, t t} + + {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} = = {P P}_{L L,, t t} + + {P P}_{l l i i n no e e,, t t} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((22))

{P P}_{G G i i}^{min min} {μ μ}_{i i,, t t} \leq \leq {P P}_{G G i i,, t t} \leq \leq {P P}_{G G i i}^{max max} {μ μ}_{i i,, t t} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((33))

{{\begin{matrix} {P P}_{G G i i,, t t} - - {P P}_{G G i i,, t t - - 11} \leq \leq {ΔP ΔP}_{G G i i}^{u u p p} \\ {P P}_{G G i i,, t t - - 11} - - {P P}_{G G i i,, t t} \leq \leq {ΔP ΔP}_{G G i i}^{d d n no} \end{matrix} ((i i = = 11,, ... ...,, N N G G;; t t = = 11,, ... ...,, N N T T)) - - - - - - ((44))

\{\begin{matrix} {Σ Σ}_{i i = = 11}^{N N G G} (({P P}_{G G i i,, t t} - - {P P}_{G G i i}^{min min})) {μ μ}_{i i,, t t} &GreaterEqual; &Greater Equal; {λ λ}_{11} {P P}_{L L,, t t} + + {λ λ}_{22} {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} \\ {Σ Σ}_{i i = = 11}^{N N G G} (({P P}_{G G i i}^{max max} - - {P P}_{G G i i,, t t})) {μ μ}_{i i,, t t} &GreaterEqual; &Greater Equal; {λ λ}_{11} {P P}_{L L,, t t} + + {λ λ}_{22} {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} \end{matrix} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((55))

\{\begin{matrix} (({X x}_{i i,, t t - - 11}^{o o n no} - - {T T}_{i i}^{o o n no})) \cdot &Center Dot; (({μ μ}_{i i,, t t - - 11} - - {μ μ}_{i i,, t t})) &GreaterEqual; &Greater Equal; 00 \\ (({X x}_{i i,, t t - - 11}^{o o f f f f} - - {T T}_{i i}^{o o f f f f})) \cdot \cdot (({μ μ}_{i i,, t t} - - {μ μ}_{i i,, t t - - 11})) &GreaterEqual; &Greater Equal; 00 \end{matrix} ((i i = = 11,, ... ...,, N N G G;; t t = = 11,, ... ...,, N N T T)) - - - - - - ((66))

{P P}_{f f l l o o w w,, l l,, t t} \leq \leq {P P}_{f f l l o o w w,, l l}^{max max} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((77))

P _Wi,t ≤P _Wi,forecast (15)

Among them, formula (14) is the upper limit of wind power acceptance for the next day, is the upper limit of the sum of wind power receiving capacity for the next day, P _Wi,forecast is the predicted output of wind farm i; NT is the number of time periods included in the research period.

7. The multi-stage wind power acceptance range calculation method based on unit combination and economic dispatch according to claim 6, characterized in that, in step C01, the model of the minimum amount of wind power that must be connected to the next day is:

Objective function:

{minP minP}_{W W d d,, S S u u m m}^{d d n no} = = {Σ Σ}_{i i = = 11}^{N N W W} {Σ Σ}_{t t = = 11}^{N N T T} {P P}_{W W i i,, t t} - - - - - - ((1616))

Restrictions:

{Σ Σ}_{i i = = 11}^{N N G G} {P P}_{G G i i,, t t} + + {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} = = {P P}_{L L,, t t} + + {P P}_{l l i i n no e e,, t t} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((22))

{P P}_{G G i i}^{min min} {μ μ}_{i i,, t t} \leq \leq {P P}_{G G i i,, t t} \leq \leq {P P}_{G G i i}^{max max} {μ μ}_{i i,, t t} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((33))

\{\begin{matrix} {P P}_{G G i i,, t t} - - {P P}_{G G i i,, t t - - 11} \leq \leq {ΔP ΔP}_{G G i i}^{u u p p} \\ {P P}_{G G i i,, t t - - 11} - - {P P}_{G G i i,, t t} \leq \leq {ΔP ΔP}_{G G i i}^{d d n no} \end{matrix} ((i i = = 11,, ... ...,, N N G G;; t t = = 11,, ... ...,, N N T T)) - - - - - - ((44))

\{\begin{matrix} {Σ Σ}_{i i = = 11}^{N N G G} (({P P}_{G G i i,, t t} - - {P P}_{G G i i}^{min min})) {μ μ}_{i i,, t t} &GreaterEqual; &Greater Equal; {λ λ}_{11} {P P}_{L L,, t t} + + {λ λ}_{22} {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} \\ {Σ Σ}_{i i = = 11}^{N N G G} (({P P}_{G G i i}^{max max} - - {P P}_{G G i i,, t t})) {μ μ}_{i i,, t t} &GreaterEqual; &Greater Equal; {λ λ}_{11} {P P}_{L L,, t t} + + {λ λ}_{22} {Σ Σ}_{j j = = 11}^{N N W W} {P P}_{W W j j,, t t} \end{matrix} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((55))

\{\begin{matrix} (({X x}_{i i,, t t - - 11}^{o o n no} - - {T T}_{i i}^{o o n no})) \cdot &Center Dot; (({μ μ}_{i i,, t t - - 11} - - {μ μ}_{i i,, t t})) &GreaterEqual; &Greater Equal; 00 \\ (({X x}_{i i,, t t - - 11}^{o o f f f f} - - {T T}_{i i}^{o o f f f f})) \cdot \cdot (({μ μ}_{i i,, t t} - - {μ μ}_{i i,, t t - - 11})) &GreaterEqual; &Greater Equal; 00 \end{matrix} ((i i = = 11,, ... ...,, N N G G;; t t = = 11,, ... ...,, N N T T)) - - - - - - ((66))

{P P}_{f f l l o o w w,, l l,, t t} \leq \leq {P P}_{f f l l o o w w,, l l}^{max max} ((t t = = 11,, ... ...,, N N T T)) - - - - - - ((77))

P _Wi,t ≤P _Wi,forecast (15)

Among them, formula (16) is the lower limit of wind power acceptance for the next day, It is the lower limit of the sum of wind power receiving capacity for the next day.

8. The multi-stage wind power acceptance range calculation method based on unit combination and economic dispatch according to any one of claims 1-7, characterized in that ILOG/CPLEX is used to solve the model.