CN105356502A - Contact line steady-state limit determination method suitable for wind power utilization - Google Patents

Abstract

Translated fromChinese

本发明提供一种适用于风电消纳的联络线稳态极限确定方法，所述方法包括：(1)计算电网系统中风电节点和负荷节点的波动功率引起的目标联络线波动功率Pt；(2)计算电网系统对所述目标联络线波动功率的可调速率；(3)计算电网系统对所述目标联络线波动功率的可控容量；(4)计算电网系统对所述目标联络线波动功率的临界可控周期Tg；(5)以1/Tg为率波频率采用切比雪夫滤波器对所述目标联络线波动功率进行滤波，得到高频分量和低频分量两部分；(6)求出目标联络线稳态传输极限。本发明的稳态极限确定方法只需为不可控波动分量留出输电余量即可，为风电消纳提供了更大输电空间。

The present invention provides a method for determining the steady-state limit of tie-lines suitable for wind power consumption, the method comprising: (1) calculating the target tie-line fluctuating power Pt caused by fluctuating power of wind power nodes and load nodes in the grid system; (2) ) Calculate the adjustable rate of the power grid system to the fluctuation power of the target tie line; (3) Calculate the controllable capacity of the grid system to the fluctuating power of the target tie line; (4) Calculate the controllable capacity of the power grid system to the fluctuation power of the target tie line; The critical controllable period Tg of the power; (5) take 1/Tg as the rate wave frequency and use the Chebyshev filter to filter the fluctuation power of the target tie line to obtain two parts, the high-frequency component and the low-frequency component; (6) find out of the steady-state transmission limit of the target tie line. The steady-state limit determination method of the present invention only needs to reserve a power transmission margin for uncontrollable fluctuation components, and provides a larger power transmission space for wind power accommodation.

Description

Translated fromChinese

一种适用于风电消纳的联络线稳态极限确定方法A tie-line steady-state limit determination method suitable for wind power consumption

技术领域technical field

本发明属于涉及一种联络线稳态极限确定方法，具体涉及一种适用于风电消纳的联络线稳态极限确定方法。The invention relates to a method for determining the steady-state limit of a tie line, in particular to a method for determining the steady-state limit of a tie line suitable for wind power consumption.

背景技术Background technique

我国风电产业持续快速发展，2014年新增风电装机量刷新历史记录。但我国风电资源主要集中在“三北”地区，与负荷中心在空间上呈逆向分布，在大规模风电入网后，当本区域风电消纳能力不足时，为避免大量弃风，势必会引起风电跨区消纳的问题。本研究所谓“网网互动”而进行的风电跨区消纳主要通过联络线完成，为了电网安全运行和更合理地安排风电交换功率，需要提前确定联络线极限传输功率。my country's wind power industry continues to develop rapidly, and the newly installed wind power capacity in 2014 broke the historical record. However, my country's wind power resources are mainly concentrated in the "Three North" regions, which are distributed inversely to the load center in space. After large-scale wind power is connected to the grid, when the wind power consumption capacity in this region is insufficient, in order to avoid a large number of abandoned wind, it will inevitably lead to wind power. The problem of cross-regional consumption. The so-called "grid-grid interaction" in this study is mainly done through the tie-line. In order to ensure the safe operation of the power grid and arrange the exchange power of wind power more reasonably, it is necessary to determine the limit transmission power of the tie-line in advance.

目前针对极限传输功率研究主要有极限传输容量TTC(TotalTransferCapability—简称TTC)研究，既有传统约束下的TTC研究，也包括大规模风电引入后的TTC研究。TTC是指在电力系统一定稳态或暂态约束下，某个输电断面能够传输的极限功率，其侧重点在于考虑系统各类安全约束，而并非联络线自身的输电极限。At present, the research on the ultimate transmission power mainly includes the research on the ultimate transmission capacity TTC (Total Transfer Capability—referred to as TTC), which includes both the TTC research under the traditional constraints and the TTC research after the introduction of large-scale wind power. TTC refers to the limit power that can be transmitted by a certain transmission section under certain steady-state or transient constraints of the power system. It focuses on considering various safety constraints of the system, rather than the transmission limit of the tie line itself.

目前针对联络线自身功率极限，一般通过计算载流量的方法进行分析，主要有静态热定值法和动态热定值法。电流流经输电元件(以下简称导体)会引起导体温度的升高，同时引起导体弧垂和应力的增大，这些来自导体本身的机械和物理的限制都可以转化为导体温度的限制。而对于导体而言，其温度是由通过导体的电流、光照、对流散热和辐射散热的共同作用而决定的。所谓静态热定值就是指导体载流与其温度同步时，导体最大允许温度所对应的载流值。这也是在电力工程界人们对导体载流热定值起初的基本认识和工程实施的依据。静态热定值可离线整定，也可在线整定。前者是环境条件固定情况下人们的做法，即传统的最大允许载流量，后者则是在实时量测基础上的做法，是实时热定值技术实现功能的一种。动态热定值是从动态热平衡方程出发，考虑导体载流变化与温度变化间的不同步性，体现导体温度在不同持续时间的界定，即以温度来表征导体的热定值。之所以称之为动态，是因为其体现的是温度变动的暂态过程，该过程结束，又回到静态热定值，所以该概念往往都对应着一个延续的时间。动态热定值法改进了静态热定值过于保守的不足，根据导线在线运行状态和气象条件等实时确定线路载流能力，一定限度地提高了联络线的输电能力。静态热定值和动态热定值都是根据联络线热载荷能力确定其输电极限，即联络线的传输功率的物理极限。At present, the power limit of the tie line itself is generally analyzed by the method of calculating the ampacity, mainly including the static thermal setting method and the dynamic thermal setting method. The current flowing through the power transmission element (hereinafter referred to as the conductor) will cause the temperature of the conductor to rise, and at the same time cause the conductor sag and stress to increase. These mechanical and physical limitations from the conductor itself can be converted into conductor temperature limitations. For a conductor, its temperature is determined by the combined effects of current passing through the conductor, light, convective heat dissipation, and radiation heat dissipation. The so-called static thermal rating is the current-carrying value corresponding to the maximum allowable temperature of the conductor when the conductor-carrying current is synchronized with its temperature. This is also the basic understanding and engineering implementation basis of people in the field of electric power engineering on the thermal rating of conductors. The static thermal setting can be set offline or online. The former is the practice of people under fixed environmental conditions, that is, the traditional maximum allowable ampacity, while the latter is based on real-time measurement, which is a kind of function realized by real-time thermal setting technology. The dynamic heat setting value starts from the dynamic heat balance equation, considers the asynchrony between the conductor current-carrying change and the temperature change, and reflects the definition of the conductor temperature at different durations, that is, the temperature is used to characterize the conductor's heat setting value. The reason why it is called dynamic is because it reflects the transient process of temperature change. After the process ends, it returns to the static thermal value, so this concept often corresponds to a continuous time. The dynamic thermal setting method improves the deficiency that the static thermal setting is too conservative, and determines the current carrying capacity of the line in real time according to the online operation status of the conductor and weather conditions, and improves the transmission capacity of the tie line to a certain extent. Both the static thermal rating and the dynamic thermal rating determine the power transmission limit according to the thermal load capacity of the tie-line, that is, the physical limit of the transmission power of the tie-line.

与传统负荷相比，风电表现出更强的随机性和不同时间尺度上的波动性，但预测水平却低于负荷预测。大规模风电并网给联络线带来大量不确定性波动，如果仍以联络线物理极限作为功率传输极限，则可能在实际运行中由于未充分考虑联络线带来大量不确定性波动，引起联络线越限等问题；若完全考虑联络线不确定性波动带来的影响，则可能导致联络线传输极限过度下降，引起输电阻塞等问题。因此有必要考虑风电强不确定性的影响，探索更加安全合理的联络线传输极限确定方法，保证风电消纳和系统安全运行。Compared with traditional load, wind power shows stronger randomness and volatility on different time scales, but the forecasting level is lower than load forecasting. Large-scale wind power grid integration brings a lot of uncertainty fluctuations to the tie line. If the physical limit of the tie line is still used as the power transmission limit, it may cause a lot of uncertainty fluctuations caused by the tie line in actual operation due to insufficient consideration of the tie line. If the influence of the uncertainty fluctuation of the tie-line is fully considered, it may lead to an excessive drop in the transmission limit of the tie-line, causing problems such as transmission congestion. Therefore, it is necessary to consider the influence of wind power intensity uncertainty and explore a safer and more reasonable method for determining the transmission limit of tie-lines to ensure wind power accommodation and safe operation of the system.

发明内容Contents of the invention

为了克服上述现有技术的不足，本发明提供一种适用于风电消纳的联络线稳态极限确定方法，本发明的稳态极限确定方法只需为不可控波动分量留出输电余量即可，为风电消纳提供了更大输电空间。In order to overcome the deficiencies of the above-mentioned prior art, the present invention provides a tie-line steady-state limit determination method suitable for wind power consumption. The steady-state limit determination method of the present invention only needs to reserve a power transmission margin for the uncontrollable fluctuation component , providing a larger power transmission space for wind power consumption.

为了实现上述发明目的，本发明采取如下技术方案：In order to realize the above-mentioned purpose of the invention, the present invention takes the following technical solutions:

一种适用于风电消纳的联络线稳态极限确定方法，所述方法包括如下步骤：A method for determining the steady-state limit of a tie line suitable for wind power accommodation, the method includes the following steps:

(1)计算电网系统中风电节点和负荷节点的波动功率引起的目标联络线波动功率Pt；(1) Calculate the target tie-line fluctuating power Pt caused by the fluctuating power of wind power nodes and load nodes in the grid system;

(2)计算电网系统对所述目标联络线波动功率的可调速率；(2) Calculate the adjustable rate of the power grid system to the fluctuating power of the target tie line;

(3)计算电网系统对所述目标联络线波动功率的可控容量；(3) Calculate the controllable capacity of the power grid system to the fluctuating power of the target tie line;

(4)计算电网系统对所述目标联络线波动功率的临界可控周期Tg；(4) Calculate the critical controllable period Tg of the grid system to the fluctuating power of the target tie line;

(5)以1/Tg为率波频率采用切比雪夫滤波器对所述目标联络线波动功率进行滤波，得到高频分量和低频分量两部分；(5) take 1/Tg as the wave frequency and adopt the Chebyshev filter to filter the fluctuation power of the target tie line to obtain two parts, high-frequency component and low-frequency component;

(6)求出目标联络线稳态传输极限。(6) Obtain the steady-state transmission limit of the target tie line.

优选的，所述步骤(1)中，所述目标联络线波动功率Pt的计算公式如下：Preferably, in the step (1), the calculation formula of the target tie line fluctuating power Pt is as follows:

${\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>t</span>t</span>}}<span><span>=</span>=</span>\underset{\mathrm{<span><span>i</span>i</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>N</span>N</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}{\mathrm{<span><span>A</span>A</span>}}_{\mathrm{<span><span>i</span>i</span>}<span><span>-</span>-</span>\mathrm{<span><span>l</span>l</span>}}{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>i</span>i</span>}}$

式中Pt为目标联络线随机波动功率，N为不确定性节点集合，P_i为风电或负荷不确定性节点功率的预测值，A_i-l为节点i对支路l的转移系数，即电网系统中节点i增加单位有功出力所引起的支路l有功潮流变化量。In the formula, Pt is the random fluctuating power of the target tie line, N is the set of uncertain nodes, P_i is the predicted value of wind power or load uncertain node power, A_il is the transfer coefficient of node i to branch l, that is, the grid system The change in active power flow of branch l caused by the increase of unit active output at node i.

优选的，所述步骤(2)中，计算所述目标联络线波动功率的可调速率包括计算电网系统中的所有可调发电机组控制节点对所述目标联络线波动功率的最大功率调节速率R，以及各可调发电机组控制节点的最优调节速率r_g，计算公式如下：Preferably, in the step (2), calculating the adjustable rate of the target tie line fluctuating power includes calculating the maximum power adjustment rate of all adjustable generator set control nodes in the grid system to the target tie line fluctuating power R, and the optimal adjustment rate r_g of each adjustable generator set control node, the calculation formula is as follows:

$\begin{array}{ccc}\mathrm{<span><span>max</span>max</span>}& \mathrm{<span><span>R</span>R</span>}<span><span>=</span>=</span><span><span>-</span>-</span>\underset{\mathrm{<span><span>g</span>g</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>P</span>P</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}{\mathrm{<span><span>A</span>A</span>}}_{\mathrm{<span><span>g</span>g</span>}<span><span>-</span>-</span>\mathrm{<span><span>l</span>l</span>}}{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>g</span>g</span>}}& <span><span>(</span>(</span>\mathrm{<span><span>1</span>1</span>}<span><span>)</span>)</span>\\ \mathrm{<span><span>s</span>the\; s</span>}<span><span>.</span>.</span>\mathrm{<span><span>t</span>t</span>}<span><span>.</span>.</span>& \underset{\mathrm{<span><span>g</span>g</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>P</span>P</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}{\mathrm{<span><span>A</span>A</span>}}_{\mathrm{<span><span>g</span>g</span>}<span><span>-</span>-</span>\mathrm{<span><span>l</span>l</span>}}{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>g</span>g</span>}}<span><span>-</span>-</span>\underset{\mathrm{<span><span>i</span>i</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>N</span>N</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}{\mathrm{<span><span>A</span>A</span>}}_{\mathrm{<span><span>i</span>i</span>}<span><span>-</span>-</span>\mathrm{<span><span>l</span>l</span>}}{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>i</span>i</span>}}<span><span>\&le;</span>\&le;</span>\mathrm{<span><span>0</span>0</span>}& <span><span>(</span>(</span>\mathrm{<span><span>2</span>2</span>}<span><span>)</span>)</span>\\ & \underset{\mathrm{<span><span>i</span>i</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>N</span>N</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>i</span>i</span>}}<span><span>+</span>+</span>\underset{\mathrm{<span><span>g</span>g</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>P</span>P</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>g</span>g</span>}}<span><span>=</span>=</span>\mathrm{<span><span>0</span>0</span>}& <span><span>(</span>(</span>\mathrm{<span><span>3</span>3</span>}<span><span>)</span>)</span>\\ & \frac{{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>i</span>i</span>}}}{{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>j</span>j</span>}}}<span><span>=</span>=</span>\frac{{\mathrm{<span><span>V</span>V</span>}}_{\mathrm{<span><span>i</span>i</span>}}}{{\mathrm{<span><span>V</span>V</span>}}_{\mathrm{<span><span>j</span>j</span>}}}<span><span>(</span>(</span>\mathrm{<span><span>i</span>i</span>}<span><span>,</span>,</span>\mathrm{<span><span>j</span>j</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>N</span>N</span>}<span><span>)</span>)</span>& <span><span>(</span>(</span>\mathrm{<span><span>4</span>4</span>}<span><span>)</span>)</span>\\ & <span><span>|</span>|</span>{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>g</span>g</span>}}<span><span>|</span>|</span><span><span>\&le;</span>\&le;</span>{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>g</span>g</span>}\mathrm{<span><span>u</span>u</span>}}<span><span>(</span>(</span>\mathrm{<span><span>g</span>g</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>P</span>P</span>}<span><span>)</span>)</span>& <span><span>(</span>(</span>\mathrm{<span><span>5</span>5</span>}<span><span>)</span>)</span>\end{array}$

式中A为发电转移分布因子矩阵，其中A_i-l为节点i对支路l的转移系数，A_g-l为可调发电机组控制节点g对支路l的转移系数；r_i,r_j为负荷或者风电随机波动节点i、j功率波动速率；r_g为可调发电机组控制节点功率调节速率；r_gu为可调发电机组控制节点功率调节速率上限；P为可调发电机组控制节点集合；N为引发不确定性的功率的风电和负荷波动节点集合；V_i,V_j为风电和负荷不确定性节点i、j的最大功率波动幅值。In the formula, A is the distribution factor matrix of power generation transfer, where A_il is the transfer coefficient of node i to branch l, A_gl is the transfer coefficient of control node g of adjustable generator set to branch l; r_i , r_j are the load or The power fluctuation rate of wind power random fluctuation nodes i and j; r_g is the power adjustment rate of the control node of the adjustable generator set; r_gu is the upper limit of the power adjustment rate of the control node of the adjustable generator set; P is the set of control nodes of the adjustable generator set; A set of wind power and load fluctuation nodes that cause uncertainty; V_i , V_j are the maximum power fluctuation amplitudes of wind power and load uncertainty nodes i and j.

优选的，所述步骤(3)中，所述目标联络线波动功率的可控容量为所有发电机组控制节点在所述目标联络线上达到的最大控制容量Psmax，计算过程如下：Preferably, in the step (3), the controllable capacity of the fluctuating power of the target tie line is the maximum control capacity Psmax achieved by all generator control nodes on the target tie line, and the calculation process is as follows:

步骤3-1、计算系统控制节点容量的可用系数K_gStep 3-1. Calculate the available coefficient K_g of the system control node capacity

${\mathrm{<span><span>K</span>K</span>}}_{\mathrm{<span><span>g</span>g</span>}}<span><span>=</span>=</span>\underset{\mathrm{<span><span>g</span>g</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>P</span>P</span>}}{\mathrm{<span><span>m</span>m</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>(</span>(</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>g</span>g</span>}\mathrm{<span><span>max</span>max</span>}}<span><span>/</span>/</span>{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>g</span>g</span>}}<span><span>)</span>)</span>$

式中，Pgmax为各控制节点的最大控制容量，r_g为可调发电机组控制节点g的最优的调节速率，P为系统控制节点的集合，min(.)为求取最小值函数；In the formula, Pgmax is the maximum control capacity of each control node, r_g is the optimal adjustment rate of the control node g of the adjustable generating set, P is the set of system control nodes, and min(.) is the function for finding the minimum value;

步骤3-2、计算系统对目标联络线的最大控制容量PsmaxStep 3-2. Calculate the maximum control capacity Psmax of the system to the target tie line

${\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>max</span>max</span>}}<span><span>=</span>=</span>\underset{\mathrm{<span><span>g</span>g</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>P</span>P</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}{\mathrm{<span><span>A</span>A</span>}}_{\mathrm{<span><span>g</span>g</span>}<span><span>-</span>-</span>\mathrm{<span><span>l</span>l</span>}}{\mathrm{<span><span>K</span>K</span>}}_{\mathrm{<span><span>g</span>g</span>}}{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>g</span>g</span>}}$

式中，A_g-l为可调发电机组控制节点g对支路l的转移系数，即电网系统中可调发电机组控制节点g增加单位有功出力所引起的支路l有功潮流变化量，r_g为可调发电机组控制节点g的最优的调节速率r_g，P为系统控制节点的集合，Psmax为所有发电机组控制节点在所述目标联络线上达到的最大控制容量。In the formula, A_gl is the transfer coefficient of the control node g of the adjustable generator set to the branch l, that is, the active power flow change of the branch l caused by the increase of the unit active output of the control node g of the adjustable generator set in the power grid system, and r_g is The optimal adjustment rate r_g of the control node g of the generator set can be adjusted, P is the set of system control nodes, and Psmax is the maximum control capacity achieved by all the control nodes of the generator set on the target tie line.

优选的，所述步骤(4)中，包括如下步骤：Preferably, in described step (4), comprise the following steps:

步骤4-1、对未来时段内的目标联络线波动功率Pt进行傅立叶级数分解，分解得到若干个不同频率和不同幅值的正弦波；将频率从小到大排列，记为f₀,f₁,f₂…f_n…，对应的幅值分别记为A₀,A₁,A₂…A_n…；Step 4-1. Perform Fourier series decomposition on the fluctuation power Pt of the target tie line in the future period, and decompose to obtain several sine waves with different frequencies and different amplitudes; arrange the frequencies from small to large, and record them as f₀ , f₁ ,f₂ …f_n …, the corresponding amplitudes are recorded as A₀ , A₁ , A₂ …A_n …;

步骤4-2、重新构造一组正弦曲线P_s，其中频率分别取f₀,f₁,f₂…f_n…，对应的幅值为$A_0,\sqrt{{A_0}^2+{A_1}^2},\sqrt{{A_0}^2+{A_1}^2+{A_2}^2}...\sqrt{{A_0}^2+{A_1}^2+{A_2}^2+\&CenterDot;\&CenterDot;\&CenterDot;{A_n}^2\mathrm{...}}\mathrm{....};$Step 4-2. Reconstruct a set of sinusoidal curves P_s , where the frequencies are respectively f₀ , f₁ , f₂ ... f_n ... and the corresponding amplitudes are$A_0,\sqrt{{A_0}^2+{A_1}^2},\sqrt{{A_0}^2+{A_1}^2+{A_2}^2}...\sqrt{{A_0}^2+{A_1}^2+{A_2}^2+\&Center\; Dot;\&Center\; Dot;\&Center\; Dot;{A_{\mathrm{no}}}^2\mathrm{...}}\mathrm{....};$

步骤4-3、根据调节速率R计算系统正弦曲线组P_s的跟踪效果指标CR_T，直到找到f₀,f₁,f₂…f_n…中某一频率的正弦曲线，使得CR_T＝ε，则该频率对应的周期即临界可控周期T_g，其中，ε为控制效果临界值，用来衡量机组是否能够跟踪给定周期正弦波，取值取决于电网实际运行状况，对于周期为T的正弦波，若跟踪率CR_T＞ε，则表示周期为T的正弦波可控；若CR_T＜ε，则表示周期为T的正弦波不可控。Step 4-3. Calculate the tracking effect index_CRT of the sinusoidal group P_s of the system according to the adjustment rate R until a sinusoidal curve of a certain frequency among f₀ , f₁ , f₂ ... f_n ... is found, so that_CRT = ε , then the period corresponding to this frequency is the critical controllable period T_g , where ε is the critical value of the control effect, which is used to measure whether the unit can track the sine wave of a given period, and the value depends on the actual operation status of the power grid. For the period T If the tracking rate CR_T >ε, it means the sine wave with period T is controllable; if CR_T <ε, it means the sine wave with period T is uncontrollable.

优选的，所述步骤4-3中，所述根据调节速率R计算系统正弦曲线组P_s的跟踪效果指标CR_T，包括如下步骤：Preferably, in the step 4-3, the calculation of the tracking effect index_CRT of the system sinusoidal group P_s according to the adjustment rate R includes the following steps:

步骤4-3-1、初始时刻区域调控机组功率P_g(0)＝0，设定积分步长为dt，其中T/dt＞1000；Step 4-3-1. At the initial moment, the regional control unit power P_g (0) = 0, and the integration step is set to dt, where T/dt>1000;

步骤4-3-2、在功率上升部分，当前仿真周期机组功率为P_g(t)，下一仿真周期正弦波动功率为P_s(t+1)，机组最大调节速率为R_g，若P_s(t+1)-P_g(t)≥R_gdt，则下一仿真周期机组功率P_g(t+1)＝P_g(t)+R_gdt，否则下一仿真周期机组功率为：P_g(t+1)＝P_s(t+1)；Step 4-3-2. In the power increase part, the unit power in the current simulation cycle is P_g (t), the sinusoidal fluctuation power in the next simulation cycle is P_s (t+1), and the maximum adjustment rate of the unit is R_g . If P_s (t+1)-P_g (t)≥R_g dt, then the power of the unit in the next simulation cycle P_g (t+1)=P_g (t)+R_g dt, otherwise the power of the unit in the next simulation cycle is : P_g (t+1) = P_s (t+1);

步骤4-3-3、在功率下降部分，当前仿真周期机组功率为P_g(t)，下一仿真周期正弦波动功率为P_s(t+1)，机组最大调节速率为R_g，若P_s(t+1)-P_g(t)≤-R_gdt，则下一仿真周期机组功率P_g(t+1)＝P_g(t)-R_gdt，否则下一仿真周期机组功率为：P_g(t+1)＝P_s(t+1)；Step 4-3-3. In the power drop part, the unit power in the current simulation cycle is P_g (t), the sinusoidal fluctuation power in the next simulation cycle is P_s (t+1), and the maximum adjustment rate of the unit is R_g . If P_s (t+1)-P_g (t)≤-R_g dt, then the power of the unit in the next simulation cycle P_g (t+1)=P_g (t)-R_g dt, otherwise the power of the unit in the next simulation cycle It is: P_g (t+1) = P_s (t+1);

步骤4-3-4、重复步骤4-3-2和步骤4-3-3，直至仿真到时间T结束；Step 4-3-4, repeat step 4-3-2 and step 4-3-3 until the simulation ends at time T;

步骤4-3-5、计算机组功率积分和正弦波动功率积分根据机组功率积分和正弦波动功率积分计算跟踪率：Step 4-3-5, computer group power integration and sinusoidal fluctuation power integral Calculate the tracking rate based on unit power integral and sinusoidal fluctuation power integral:

${\mathrm{<span><span>CR</span>CR</span>}}_{\mathrm{<span><span>T</span>T</span>}}<span><span>=</span>=</span>\frac{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>g</span>g</span>}\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>u</span>u</span>}\mathrm{<span><span>m</span>m</span>}}}{{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>u</span>u</span>}\mathrm{<span><span>m</span>m</span>}}}<span><span>\&times;</span>\&times;</span>\mathrm{<span><span>100</span>100</span>}\mathrm{<span><span>\%</span>\%</span>}<span><span>.</span>.</span>$

优选的，所述步骤(6)包括如下步骤：Preferably, said step (6) comprises the steps of:

步骤6-1、将所述高频分量作为第一部分不可控分量；Step 6-1, using the high-frequency component as the first uncontrollable component;

步骤6-2、将所述低频分量中幅值大于系统控制容量Psmax的部分作为第二部分不可控分量，将其余部分作为可控分量；Step 6-2, taking the part of the low-frequency component whose amplitude is greater than the system control capacity Psmax as the second uncontrollable component, and taking the rest as the controllable component;

步骤6-3、将所述第一部分不可控分量和所述第二部分不可控分量叠加作为目标联络线波动中的不可控总量；Step 6-3, superimposing the first part of the uncontrollable component and the second part of the uncontrollable component as the uncontrollable total amount in the fluctuation of the target tie line;

步骤6-4、将已知的联络线物理极限减去所述不可控总量，得到考虑系统调节速率下的联络线稳态传输极限。Step 6-4: Subtract the uncontrollable total amount from the known physical limit of the tie line to obtain the steady-state transmission limit of the tie line considering the adjustment rate of the system.

与现有技术相比，本发明的有益效果在于：Compared with prior art, the beneficial effect of the present invention is:

本发明在引入风电不确定性波动后，与将物理极限作为稳态传输极限相比，确定方法更能保证系统的安全性；本发明的稳态极限确定方法只需为不可控波动分量留出输电余量即可，为风电消纳提供了更大输电空间，更适用于风电的消纳。Compared with using the physical limit as the steady-state transmission limit after the wind power uncertainty fluctuation is introduced in the present invention, the determination method can better ensure the safety of the system; the steady-state limit determination method of the present invention only needs to set aside for uncontrollable fluctuation components The power transmission margin is enough, which provides a larger transmission space for wind power consumption and is more suitable for wind power consumption.

附图说明Description of drawings

图1是本发明提供的一种适用于风电消纳的联络线稳态极限确定方法流程图Fig. 1 is a flow chart of a method for determining the steady-state limit of tie-lines suitable for wind power consumption provided by the present invention

具体实施方式detailed description

下面结合附图对本发明作进一步详细说明。The present invention will be described in further detail below in conjunction with the accompanying drawings.

如图1所示，一种适用于风电消纳的联络线稳态极限确定方法，本方法包括如下步骤：As shown in Figure 1, a method for determining the steady-state limit of tie-lines suitable for wind power consumption, this method includes the following steps:

1、计算电网系统中风电节点和负荷节点的波动功率引起的目标联络线波动功率Pt；1. Calculate the target tie line fluctuating power Pt caused by the fluctuating power of wind power nodes and load nodes in the power grid system;

计算电网系统中风电和负荷节点的不确定性波动引起目标联络线的功率波动Pt，计算公式为：Calculate the power fluctuation Pt of the target tie line caused by the uncertainty fluctuation of wind power and load nodes in the grid system, and the calculation formula is:

${\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>t</span>t</span>}}<span><span>=</span>=</span>\underset{\mathrm{<span><span>i</span>i</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>N</span>N</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}{\mathrm{<span><span>A</span>A</span>}}_{\mathrm{<span><span>i</span>i</span>}<span><span>-</span>-</span>\mathrm{<span><span>l</span>l</span>}}{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>i</span>i</span>}}$

其中Pt为联络线随机波动功率，N为不确定性节点集合，Pi为风电或负荷不确定性节点功率的预测值。A_i-l为发电转移分布因子GSDF矩阵元素，含义是节点i对支路l的转移系数，定义为：电网系统中节点i增加单位有功出力所引起的支路l有功潮流变化量。Among them, Pt is the random fluctuating power of the tie line, N is the set of uncertain nodes, and Pi is the predicted value of wind power or load uncertain node power. A_il is the matrix element of the generation transfer distribution factor GSDF, which means the transfer coefficient of node i to branch l, defined as: the active power flow change of branch l caused by the increase of unit active output of node i in the power grid system.

2、计算电网系统对所述目标联络线波动功率的可调速率；2. Calculate the adjustable rate of the grid system to the fluctuating power of the target tie line;

根据上述要求，针对(1)-(5)，采用线性规划的方法计算系统对于目标联络线的最大调节速率R以及各可调发电机组控制节点的最优的调节速率r_g：According to the above requirements, for (1)-(5), the linear programming method is used to calculate the maximum adjustment rate R of the system for the target tie line and the optimal adjustment rate r_g of the control nodes of each adjustable generating set:

$\begin{array}{ccc}\mathrm{<span><span>max</span>max</span>}& \mathrm{<span><span>R</span>R</span>}<span><span>=</span>=</span><span><span>-</span>-</span>\underset{\mathrm{<span><span>g</span>g</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>P</span>P</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}{\mathrm{<span><span>A</span>A</span>}}_{\mathrm{<span><span>g</span>g</span>}<span><span>-</span>-</span>\mathrm{<span><span>l</span>l</span>}}{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>g</span>g</span>}}& <span><span>(</span>(</span>\mathrm{<span><span>1</span>1</span>}<span><span>)</span>)</span>\\ \mathrm{<span><span>s</span>the\; s</span>}<span><span>.</span>.</span>\mathrm{<span><span>t</span>t</span>}<span><span>.</span>.</span>& \underset{\mathrm{<span><span>g</span>g</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>P</span>P</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}{\mathrm{<span><span>A</span>A</span>}}_{\mathrm{<span><span>g</span>g</span>}<span><span>-</span>-</span>\mathrm{<span><span>l</span>l</span>}}{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>g</span>g</span>}}<span><span>-</span>-</span>\underset{\mathrm{<span><span>i</span>i</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>N</span>N</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}{\mathrm{<span><span>A</span>A</span>}}_{\mathrm{<span><span>i</span>i</span>}<span><span>-</span>-</span>\mathrm{<span><span>l</span>l</span>}}{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>i</span>i</span>}}<span><span>\&le;</span>\&le;</span>\mathrm{<span><span>0</span>0</span>}& <span><span>(</span>(</span>\mathrm{<span><span>2</span>2</span>}<span><span>)</span>)</span>\\ & \underset{\mathrm{<span><span>i</span>i</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>N</span>N</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>i</span>i</span>}}<span><span>+</span>+</span>\underset{\mathrm{<span><span>g</span>g</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>P</span>P</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>g</span>g</span>}}<span><span>=</span>=</span>\mathrm{<span><span>0</span>0</span>}& <span><span>(</span>(</span>\mathrm{<span><span>3</span>3</span>}<span><span>)</span>)</span>\\ & \frac{{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>i</span>i</span>}}}{{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>j</span>j</span>}}}<span><span>=</span>=</span>\frac{{\mathrm{<span><span>V</span>V</span>}}_{\mathrm{<span><span>i</span>i</span>}}}{{\mathrm{<span><span>V</span>V</span>}}_{\mathrm{<span><span>j</span>j</span>}}}<span><span>(</span>(</span>\mathrm{<span><span>i</span>i</span>}<span><span>,</span>,</span>\mathrm{<span><span>j</span>j</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>N</span>N</span>}<span><span>)</span>)</span>& <span><span>(</span>(</span>\mathrm{<span><span>4</span>4</span>}<span><span>)</span>)</span>\\ & <span><span>|</span>|</span>{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>g</span>g</span>}}<span><span>|</span>|</span><span><span>\&le;</span>\&le;</span>{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>g</span>g</span>}\mathrm{<span><span>u</span>u</span>}}<span><span>(</span>(</span>\mathrm{<span><span>g</span>g</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>P</span>P</span>}<span><span>)</span>)</span>& <span><span>(</span>(</span>\mathrm{<span><span>5</span>5</span>}<span><span>)</span>)</span>\end{array}$

其中A为发电转移分布因子矩阵，其中A_i-l为节点i对支路l的转移系数，A_i-l为节点g对支路l的转移系数；r_i,r_j为负荷或者风电随机波动节点i、j功率波动速率；r_g为可调发电机组控制节点功率调节速率；r_gu为可调发电机组控制节点功率调节速率上限；P为可调发电机组控制节点集合；N为引发不确定性的功率的风电和负荷波动节点集合；V_i,V_j为风电和负荷不确定性节点i、j的最大功率波动幅值。where A is the power generation transfer distribution factor matrix, where A_il is the transfer coefficient of node i to branch l, and A_il is the transfer coefficient of node g to branch l; r_i and r_j are load or wind power random fluctuation nodes i, j power fluctuation rate; r_g is the power adjustment rate of the control node of the adjustable generator set; r_gu is the upper limit of the power adjustment rate of the control node of the adjustable generator set; P is the set of control nodes of the adjustable generator set; N is the power that causes uncertainty The set of wind power and load fluctuation nodes; V_i , V_j are the maximum power fluctuation amplitudes of wind power and load uncertainty nodes i and j.

式(1)表示线性规划目标函数：求取满足要求的系统最大调节速率，当通过GSDF将节点调节速率转移到联络线时，节点转移到联络线的功率正方向与初始潮流一致，而线性规划的目标是求取反向控制潮流的能力，因此加负号；式(2)表示控制节点和不确定性节点的共同作用使得目标连路线功率波动为0，由于目标联络线功率由于不确定性节点当作负荷节点处理，其转移系数采用LSDF；式(3)表示系统总控制功率与总的不确定性功率保持平衡进而；式(4)表示各不确定性节点功率波动速率的比例，与其波动幅值成正比；式(5)则表示各控制节点功率调节速率的限制。Equation (1) expresses the linear programming objective function: find the maximum adjustment rate of the system that meets the requirements. When the node adjustment rate is transferred to the tie line through GSDF, the positive direction of the power transferred from the node to the tie line is consistent with the initial power flow, while the linear programming The goal of is to obtain the ability to control the power flow in the reverse direction, so the negative sign is added; Equation (2) indicates that the joint action of the control node and the uncertain node makes the power fluctuation of the target link to 0, because the power of the target link is due to the uncertainty The node is treated as a load node, and its transfer coefficient adopts LSDF; Equation (3) indicates that the total control power of the system is balanced with the total uncertain power; Equation (4) indicates the ratio of the power fluctuation rate of each uncertain node, compared Fluctuation amplitude is proportional; Equation (5) expresses the limitation of each control node power regulation rate.

3、计算电网系统对所述目标联络线波动功率的可控容量；3. Calculate the controllable capacity of the power grid system for the fluctuating power of the target tie line;

系统对联络线的可控容量取决于各控制节点的调节容量，再根据步骤2中得到的各可调发电机组控制节点的最优的调节速率r_g，可得到系统对目标联络线的控制容量Psmax。The controllable capacity of the system to the tie line depends on the adjustment capacity of each control node, and then according to the optimal adjustment rate r_g of the control nodes of each adjustable generator set obtained in step 2, the control capacity of the system to the target tie line can be obtained Psmax.

先计算系统控制节点容量的可用系数K_gFirst calculate the available coefficient K_g of the system control node capacity

${\mathrm{<span><span>K</span>K</span>}}_{\mathrm{<span><span>g</span>g</span>}}<span><span>=</span>=</span>\underset{\mathrm{<span><span>g</span>g</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>P</span>P</span>}}{\mathrm{<span><span>m</span>m</span>}\mathrm{<span><span>i</span>i</span>}\mathrm{<span><span>n</span>no</span>}}<span><span>(</span>(</span>{\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>g</span>g</span>}\mathrm{<span><span>max</span>max</span>}}<span><span>/</span>/</span>{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>g</span>g</span>}}<span><span>)</span>)</span>$

Pgmax为各控制节点的最大控制容量，r_g为可调发电机组控制节点g的最优的调节速率r_g，P为系统控制节点的集合，min(.)为求取最小值函数。Pgmax is the maximum control capacity of each control node, r_g is the optimal adjustment rate r_g of the control node g of the adjustable generating set, P is the set of system control nodes, and min(.) is the function for finding the minimum value.

再计算系统对目标联络线的控制容量PsmaxRecalculate the control capacity Psmax of the system to the target tie line

${\mathrm{<span><span>P</span>P</span>}}_{\mathrm{<span><span>s</span>the\; s</span>}\mathrm{<span><span>max</span>max</span>}}<span><span>=</span>=</span>\underset{\mathrm{<span><span>g</span>g</span>}<span><span>\&Element;</span>\&Element;</span>\mathrm{<span><span>P</span>P</span>}}{<span><span>\&Sigma;</span>\&Sigma;</span>}{\mathrm{<span><span>A</span>A</span>}}_{\mathrm{<span><span>g</span>g</span>}<span><span>-</span>-</span>\mathrm{<span><span>l</span>l</span>}}{\mathrm{<span><span>K</span>K</span>}}_{\mathrm{<span><span>g</span>g</span>}}{\mathrm{<span><span>r</span>r</span>}}_{\mathrm{<span><span>g</span>g</span>}}$

其中，A_g-l为GSDF矩阵元素，含义是可调发电机组控制节点g对支路l的转移系数，定义为：电网系统中可调发电机组控制节点g增加单位有功出力所引起的支路l有功潮流变化量.r_g为可调发电机组控制节点g的最优的调节速率r_g，P为系统控制节点的集合。Among them, A_gl is a GSDF matrix element, which means the transfer coefficient of the control node g of the adjustable generator set to the branch l, defined as: the active power of the branch l caused by the increase of the unit active output of the control node g of the adjustable generator set in the grid system Power flow variation. r_{g is the optimal regulation rate r gof} the control node g of the adjustable generating set, and P is the set of system control nodes.

4、计算电网系统对所述目标联络线波动功率的临界可控周期Tg；4. Calculate the critical controllable period Tg of the grid system to the fluctuating power of the target tie line;

具体步骤：Specific steps:

(1)对未来时段(如15min)内的目标联络线波动功率Pt进行傅立叶级数分解，分解得到若干个不同频率和不同幅值的正弦波。将频率从小到大排列，记为f₀,f₁,f₂…f_n…，对应的幅值分别记为A₀,A₁,A₂…A_n…(1) Perform Fourier series decomposition on the fluctuation power Pt of the target tie line in the future period (such as 15 minutes), and decompose to obtain several sine waves with different frequencies and different amplitudes. Arrange the frequencies from small to large, denoted as f₀ , f₁ , f₂ ...f_n ..., and the corresponding amplitudes are denoted as A₀ , A₁ , A₂ ...A_n ...

(2)重新构造一组正弦曲线P_s，其中频率分别取f₀,f₁,f₂…f_n…，对应的幅值为$A_0,\sqrt{{A_0}^2+{A_1}^2},\sqrt{{A_0}^2+{A_1}^2+{A_2}^2}...\sqrt{{A_0}^2+{A_1}^2+{A_2}^2+\&CenterDot;\&CenterDot;\&CenterDot;{A_n}^2\mathrm{...}}\mathrm{.....}$(2) Reconstruct a set of sinusoidal curves P_s , where the frequencies are respectively f₀ , f₁ , f₂ …f_n … and the corresponding amplitudes are$A_0,\sqrt{{A_0}^2+{A_1}^2},\sqrt{{A_0}^2+{A_1}^2+{A_2}^2}...\sqrt{{A_0}^2+{A_1}^2+{A_2}^2+\&CenterDot;\&CenterDot;\&CenterDot;{A_{\mathrm{no}}}^2\mathrm{...}}\mathrm{.....}$

(3)根据调节速率R计算系统对于式(2)中正弦曲线组P_s的跟踪效果指标CR_T，直到找到f₀,f₁,f₂…f_n…中某一频率的正弦曲线，使得CR_T＝ε，则该频率对应的周期即临界周期T_g。(3) Calculate the tracking effect index_CRT of the system for the sinusoidal group P_s in formula (2) according to the adjustment rate R, until the sinusoidal curve of a certain frequency among f₀ , f₁ , f₂ ... f_n ... is found, so that CR_T =ε, then the period corresponding to this frequency is the critical period T_g .

其中，ε为控制效果临界值，用来衡量机组是否能够跟踪给定周期正弦波，取值取决于电网实际运行状况；对于周期为T的正弦波，如果跟踪率CR_T＞ε，则表示周期为T的正弦波可控；反之，如果CR_T＜ε，则认为周期为T的正弦波不可控。Among them, ε is the critical value of the control effect, which is used to measure whether the unit can track the sine wave of a given period, and the value depends on the actual operation status of the power grid; for the sine wave with a period of T, if the tracking rate CR_T > ε, it means The sine wave with T is controllable; otherwise, if CR_T <ε, the sine wave with period T is considered uncontrollable.

根据调节速率计算CR_T的流程如下：The process of calculating_CRT according to the adjustment rate is as follows:

a、初始时刻区域调控机组功率P_g(0)＝0，设定积分步长为dt，其中T/dt＞1000；a. At the initial moment, the regional control unit power P_g (0) = 0, set the integral step size as dt, where T/dt>1000;

b、在功率上升部分，当前仿真周期机组功率为P_g(t),下一仿真周期正弦波动功率为P_s(t+1)，机组最大调节速率为R_g，如果P_s(t+1)-P_g(t)≥R_gdt，则下一仿真周期机组功率P_g(t+1)＝P_g(t)+R_gdt，否则下一仿真周期机组功率为：P_g(t+1)＝P_s(t+1)；b. In the power rising part, the unit power in the current simulation cycle is P_g (t), the sinusoidal fluctuation power in the next simulation cycle is P_s (t+1), and the maximum adjustment rate of the unit is R_g . If P_s (t+1 )-P_g (t)≥R_g dt, then the power of the unit in the next simulation cycle P_g (t+1)=P_g (t)+R_g dt, otherwise the power of the unit in the next simulation cycle is: P_g (t +1)=P_s (t+1);

c、在功率下降部分，当前仿真周期机组功率为P_g(t),下一仿真周期正弦波动功率为P_s(t+1)，机组最大调节速率为R_g，如果P_s(t+1)-P_g(t)≤-R_gdt,则下一仿真周期机组功率P_g(t+1)＝P_g(t)-R_gdt，否则下一仿真周期机组功率为：P_g(t+1)＝P_s(t+1)；c. In the power drop part, the power of the unit in the current simulation cycle is P_g (t), the sinusoidal fluctuation power in the next simulation cycle is P_s (t+1), and the maximum adjustment rate of the unit is R_g . If P_s (t+1 )-P_g (t)≤-R_g dt, then the power of the unit in the next simulation cycle P_g (t+1)=P_g (t)-R_g dt, otherwise the power of the unit in the next simulation cycle is: P_g ( t+1)=P_s (t+1);

d、重复步骤b和步骤c，直至仿真到时间T结束；d. Repeat step b and step c until the simulation ends at time T;

e、计算机组功率积分和正弦波动功率积分根据机组功率积分和正弦波动功率积分计算跟踪率：e. Computer group power integral and sinusoidal fluctuation power integral Calculate the tracking rate based on unit power integral and sinusoidal fluctuation power integral:

5、以1/Tg为率波频率采用切比雪夫滤波器对所述目标联络线波动功率进行滤波，得到高频分量和低频分量两部分；5. Taking 1/Tg as the rate wave frequency, the Chebyshev filter is used to filter the fluctuation power of the target tie line to obtain two parts: high-frequency component and low-frequency component;

6、求出目标联络线稳态传输极限。6. Obtain the steady-state transmission limit of the target tie line.

1)将高频分量作为第一部分不可控分量；1) Take the high-frequency component as the first uncontrollable component;

2)将低频分量中幅值大于系统控制容量Psmax的部分作为第二部分不可控分量，将其余部分作为可控分量；2) The part of the low-frequency component whose amplitude is greater than the system control capacity Psmax is taken as the second part of the uncontrollable component, and the remaining part is taken as the controllable component;

3)将第一部分不可控分量和第二部分不可控分量叠加作为目标联络线波动中的不可控总量；3) Superimposing the first part of the uncontrollable component and the second part of the uncontrollable component as the uncontrollable total amount in the fluctuation of the target tie line;

4)将物理极限减去不可控总量，即可得到考虑系统调节速率下的联络线稳态传输极限。4) Subtract the uncontrollable total amount from the physical limit to obtain the steady-state transmission limit of the tie line considering the system regulation rate.

最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求范围当中。Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit them. Although the present invention has been described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: the present invention can still be Any modification or equivalent replacement that does not depart from the spirit and scope of the present invention shall be covered by the scope of the claims of the present invention.

Claims (7)

Translated fromChinese

1.一种适用于风电消纳的联络线稳态极限确定方法，其特征在于，所述方法包括如下步骤：1. A tie-line steady-state limit determination method applicable to wind power consumption, characterized in that, the method may further comprise the steps:(1)计算电网系统中风电节点和负荷节点的波动功率引起的目标联络线波动功率Pt；(1) Calculate the target tie-line fluctuating power Pt caused by the fluctuating power of wind power nodes and load nodes in the grid system;(2)计算电网系统对所述目标联络线波动功率的可调速率；(2) Calculate the adjustable rate of the power grid system to the fluctuating power of the target tie line;(3)计算电网系统对所述目标联络线波动功率的可控容量；(3) Calculate the controllable capacity of the power grid system to the fluctuating power of the target tie line;(4)计算电网系统对所述目标联络线波动功率的临界可控周期Tg；(4) Calculate the critical controllable period Tg of the grid system to the fluctuating power of the target tie line;(5)以1/Tg为率波频率采用切比雪夫滤波器对所述目标联络线波动功率进行滤波，得到高频分量和低频分量两部分；(5) take 1/Tg as the wave frequency and adopt the Chebyshev filter to filter the fluctuation power of the target tie line to obtain two parts, high-frequency component and low-frequency component;(6)求出目标联络线稳态传输极限。(6) Obtain the steady-state transmission limit of the target tie line.2.根据权利要求1所述确定方法，其特征在于，所述步骤(1)中，所述目标联络线波动功率Pt的计算公式如下：2. according to the described determination method of claim 1, it is characterized in that, in described step (1), the computing formula of described target tie line fluctuating power Pt is as follows:式中Pt为目标联络线随机波动功率，N为不确定性节点集合，P_i为风电或负荷不确定性节点功率的预测值，A_i-l为节点i对支路l的转移系数，即电网系统中节点i增加单位有功出力所引起的支路l有功潮流变化量。In the formula, Pt is the random fluctuating power of the target tie line, N is the set of uncertain nodes, P_i is the predicted value of wind power or load uncertain node power, A_il is the transfer coefficient of node i to branch l, that is, the grid system The change in active power flow of branch l caused by the increase of unit active output at node i.3.根据权利要求1所述确定方法，其特征在于，所述步骤(2)中，计算所述目标联络线波动功率的可调速率包括计算电网系统中的所有可调发电机组控制节点对所述目标联络线波动功率的最大功率调节速率R，以及各可调发电机组控制节点的最优调节速率r_g，计算公式如下：3. The determination method according to claim 1, wherein in the step (2), calculating the adjustable rate of the target tie line fluctuating power comprises calculating the control node pairs of all adjustable generator sets in the grid system The calculation formulas for the maximum power adjustment rate R of the fluctuating power of the target tie line and the optimal adjustment rate r_g of the control nodes of each adjustable generator set are as follows:式中A_i-l为节点i对支路l的转移系数，A_g-l为可调发电机组控制节点g对支路l的转移系数；r_i,r_j为负荷或者风电随机波动节点i、j功率波动速率；r_g为可调发电机组控制节点功率调节速率；r_gu为可调发电机组控制节点功率调节速率上限；P为可调发电机组控制节点集合；N为引发不确定性的功率的风电和负荷波动节点集合；V_i,V_j为风电和负荷不确定性节点i、j的最大功率波动幅值。In the formula, A_il is the transfer coefficient of node i to branch l, A_gl is the transfer coefficient of control node g of adjustable generating set to branch l; r_i , r_j are load or wind power random fluctuation node i, j power fluctuation r_g is the power adjustment rate of the control node of the adjustable generator set; r_gu is the upper limit of the power adjustment rate of the control node of the adjustable generator set; P is the set of control nodes of the adjustable generator set; N is the wind power and A set of load fluctuation nodes; V_i , V_j are the maximum power fluctuation amplitudes of wind power and load uncertainty nodes i and j.4.根据权利要求1所述确定方法，其特征在于，所述步骤(3)中，所述目标联络线波动功率的可控容量为所有发电机组控制节点在所述目标联络线上达到的最大控制容量Psmax，计算过程如下：4. The determination method according to claim 1, characterized in that, in the step (3), the controllable capacity of the fluctuating power of the target tie line is the maximum achieved by all generator control nodes on the target tie line Control capacity Psmax, the calculation process is as follows:步骤3-1、计算系统控制节点容量的可用系数K_gStep 3-1. Calculate the available coefficient K_g of the system control node capacity式中，Pgmax为各控制节点的最大控制容量，r_g为可调发电机组控制节点g的最优的调节速率，P为系统控制节点的集合，min(.)为求取最小值函数；In the formula, Pgmax is the maximum control capacity of each control node, r_g is the optimal adjustment rate of the control node g of the adjustable generating set, P is the set of system control nodes, and min(.) is the function for finding the minimum value;步骤3-2、计算系统对目标联络线的最大控制容量PsmaxStep 3-2. Calculate the maximum control capacity Psmax of the system to the target tie line式中，A_g-l为可调发电机组控制节点g对支路l的转移系数，即电网系统中可调发电机组控制节点g增加单位有功出力所引起的支路l有功潮流变化量，r_g为可调发电机组控制节点g的最优的调节速率r_g，P为系统控制节点的集合，Psmax为所有发电机组控制节点在所述目标联络线上达到的最大控制容量。In the formula, A_gl is the transfer coefficient of the control node g of the adjustable generator set to the branch l, that is, the active power flow change of the branch l caused by the increase of the unit active output of the control node g of the adjustable generator set in the power grid system, and r_g is The optimal adjustment rate r_g of the control node g of the generator set can be adjusted, P is the set of system control nodes, and Psmax is the maximum control capacity achieved by all the control nodes of the generator set on the target tie line.5.根据权利要求1所述确定方法，其特征在于，所述步骤(4)中，包括如下步骤：5. determining method according to claim 1, is characterized in that, in described step (4), comprises the following steps:步骤4-1、对未来时段内的目标联络线波动功率Pt进行傅立叶级数分解，分解得到若干个不同频率和不同幅值的正弦波；将频率从小到大排列，记为f₀,f₁,f₂…f_n…，对应的幅值分别记为A₀,A₁,A₂…A_n…；Step 4-1. Perform Fourier series decomposition on the fluctuation power Pt of the target tie line in the future period, and decompose to obtain several sine waves with different frequencies and different amplitudes; arrange the frequencies from small to large, and record them as f₀ , f₁ ,f₂ …f_n …, the corresponding amplitudes are recorded as A₀ , A₁ , A₂ …A_n …;步骤4-2、重新构造一组正弦曲线P_s，其中频率分别取f₀,f₁,f₂…f_n…，对应的幅值为Step 4-2. Reconstruct a set of sinusoidal curves P_s , where the frequencies are respectively f₀ , f₁ , f₂ ... f_n ... and the corresponding amplitudes are步骤4-3、根据调节速率R计算系统正弦曲线组P_s的跟踪效果指标CR_T，直到找到f₀,f₁,f₂…f_n…中某一频率的正弦曲线，使得CR_T＝ε，则该频率对应的周期即临界可控周期T_g，其中，ε为控制效果临界值，用来衡量机组是否能够跟踪给定周期正弦波，取值取决于电网实际运行状况，对于周期为T的正弦波，若跟踪率CR_T＞ε，则表示周期为T的正弦波可控；若CR_T＜ε，则表示周期为T的正弦波不可控。Step 4-3. Calculate the tracking effect index_CRT of the sinusoidal group P_s of the system according to the adjustment rate R until a sinusoidal curve of a certain frequency among f₀ , f₁ , f₂ ... f_n ... is found, so that_CRT = ε , then the period corresponding to this frequency is the critical controllable period T_g , where ε is the critical value of the control effect, which is used to measure whether the unit can track the sine wave of a given period, and the value depends on the actual operation status of the power grid. For the period T If the tracking rate CR_T >ε, it means the sine wave with period T is controllable; if CR_T <ε, it means the sine wave with period T is uncontrollable.6.根据权利要求5所述确定方法，其特征在于，所述步骤4-3中，所述根据调节速率R计算系统正弦曲线组P_s的跟踪效果指标CR_T，包括如下步骤：6. The determination method according to claim 5, characterized in that, in said step 4-3, said calculation of the tracking effect index_CRT of the system sinusoidal group P_s according to the adjustment rate R comprises the following steps:步骤4-3-1、初始时刻区域调控机组功率P_g(0)＝0，设定积分步长为dt，其中T/dt＞1000；Step 4-3-1. At the initial moment, the regional control unit power P_g (0) = 0, and the integration step is set to dt, where T/dt>1000;步骤4-3-2、在功率上升部分，当前仿真周期机组功率为P_g(t)，下一仿真周期正弦波动功率为P_s(t+1)，机组最大调节速率为R_g，若P_s(t+1)-P_g(t)≥R_gdt，则下一仿真周期机组功率P_g(t+1)＝P_g(t)+R_gdt，否则下一仿真周期机组功率为：P_g(t+1)＝P_s(t+1)；Step 4-3-2. In the power increase part, the unit power in the current simulation cycle is P_g (t), the sinusoidal fluctuation power in the next simulation cycle is P_s (t+1), and the maximum adjustment rate of the unit is R_g . If P_s (t+1)-P_g (t)≥R_g dt, then the power of the unit in the next simulation cycle P_g (t+1)=P_g (t)+R_g dt, otherwise the power of the unit in the next simulation cycle is : P_g (t+1) = P_s (t+1);步骤4-3-3、在功率下降部分，当前仿真周期机组功率为P_g(t)，下一仿真周期正弦波动功率为P_s(t+1)，机组最大调节速率为R_g，若P_s(t+1)-P_g(t)≤-R_gdt，则下一仿真周期机组功率P_g(t+1)＝P_g(t)-R_gdt，否则下一仿真周期机组功率为：P_g(t+1)＝P_s(t+1)；Step 4-3-3. In the power drop part, the unit power in the current simulation cycle is P_g (t), the sinusoidal fluctuation power in the next simulation cycle is P_s (t+1), and the maximum adjustment rate of the unit is R_g . If P_s (t+1)-P_g (t)≤-R_g dt, then the power of the unit in the next simulation cycle P_g (t+1)=P_g (t)-R_g dt, otherwise the power of the unit in the next simulation cycle It is: P_g (t+1) = P_s (t+1);步骤4-3-4、重复步骤4-3-2和步骤4-3-3，直至仿真到时间T结束；Step 4-3-4, repeat step 4-3-2 and step 4-3-3 until the simulation ends at time T;步骤4-3-5、计算机组功率积分和正弦波动功率积分根据机组功率积分和正弦波动功率积分计算跟踪率：Step 4-3-5, computer group power integration and sinusoidal fluctuation power integral Calculate the tracking rate based on unit power integral and sinusoidal fluctuation power integral:7.根据权利要求1所述确定方法，其特征在于，所述步骤(6)包括如下步骤：7. according to the described determination method of claim 1, it is characterized in that, described step (6) comprises the following steps:步骤6-1、将所述高频分量作为第一部分不可控分量；Step 6-1, using the high-frequency component as the first uncontrollable component;步骤6-2、将所述低频分量中幅值大于系统控制容量Psmax的部分作为第二部分不可控分量，将其余部分作为可控分量；Step 6-2, taking the part of the low-frequency component whose amplitude is greater than the system control capacity Psmax as the second uncontrollable component, and taking the rest as the controllable component;步骤6-3、将所述第一部分不可控分量和所述第二部分不可控分量叠加作为目标联络线波动中的不可控总量；Step 6-3, superimposing the first part of the uncontrollable component and the second part of the uncontrollable component as the uncontrollable total amount in the fluctuation of the target tie line;步骤6-4、将已知的联络线物理极限减去所述不可控总量，得到考虑系统调节速率下的联络线稳态传输极限。Step 6-4: Subtract the uncontrollable total amount from the known physical limit of the tie line to obtain the steady-state transmission limit of the tie line considering the adjustment rate of the system.