CN107679769A - 含风电的电力系统频率响应模型建立方法和频率特性指标计算方法 - Google Patents
含风电的电力系统频率响应模型建立方法和频率特性指标计算方法 Download PDFInfo
- Publication number
- CN107679769A CN107679769A CN201711011404.0A CN201711011404A CN107679769A CN 107679769 A CN107679769 A CN 107679769A CN 201711011404 A CN201711011404 A CN 201711011404A CN 107679769 A CN107679769 A CN 107679769A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mfrac
- wind
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B10/00—Integration of renewable energy sources in buildings
- Y02B10/30—Wind power
Abstract
本发明公开了一种含风电的电力系统频率响应模型建立方法和频率特性指标计算方法，其中模型建立方法，包括如下步骤：1、假设传统电力系统的发电机组类型是火电机组，将占系统总发电量比例为ρ的风电机组代替同等容量的火电机组接入电网；2、风电场风电机组采用虚拟惯性和桨距角减载的频率辅助控制策略，根据风电机组的转子运动方程，并采用小信号分析方法，建立风电机组的降解动态频率响应模型；3、根据传统电力系统的频率响应模型和建立的风电机组降解动态频率响应模型，建立计及风电调频的电力系统频率响应模型。采用该方法建立的模型能够快速、精确地评估电力系统的动态频率响应特性。
Description
技术领域
本发明属于电力系统动态频率分析技术领域，具体涉及一种计及风电参与调频的电力系统频率响应模型建立方法和电网动态频率特性指标计算方法。
背景技术
大规模风电并网运行正深刻的改变着电力系统的动态运行特性，不同于传统的同步发电机，风电机组是通过电力电子变换装置进行并网，其输出的有功功率和电网频率是完全解耦的，另外风电机组一般处于最大功率运行状态，没有预留的备用功率。带来的结果是风电机组不能给电力系统的频率响应作出贡献，大规模的风电并网或负荷的大幅扰动都会给电力系统的频率带来较大的挑战。因此，有必要对风电并网系统的动态频率响应特性进行深入的研究。
为解决风电并网的频率稳定问题，国内外学者对风电机组参与频率调节的控制策略进行了大量的研究。其中最为有效的风电机组调频控制策略就是基于虚拟惯性和桨距角的联合控制策略，后续的很多研究是在该策略上进行的深入探讨。此类风电调频控制策略的研究对象都是风电机组内部的控制策略，在频率调节方面也能取得一定的改善效果。然而，却很少有大规模风电并网的电力系统频率动态响应特性的研究。
传统的电力系统动态频率响应模型的机组都是火电机组，大规模风电的接入使得该模型已不再适用。因此，为深入研究大规模风电并网的电力系统动态频率响应特性，那么很有必要建立含风电的电力系统频率响应模型。
发明内容
发明目的：针对现有技术中的问题，本发明提供了一种适用于大规模风电并网电力系统频率响应模型，能够快速、精确地评估电力系统的动态频率响应特性。
技术方案：本发明一方面公开了含风电的电力系统频率响应模型建立方法，包括如下步骤：
(1)假设传统电力系统的发电机组类型是火电机组，将占系统总发电量比例为ρ的风电机组代替同等容量的火电机组接入电网；
(2)风电场风电机组采用虚拟惯性和桨距角减载的频率辅助控制策略，根据风电机组的转子运动方程，并采用小信号分析方法，建立风电机组的降解动态频率响应模型；
(3)根据传统电力系统的频率响应模型和建立的风电机组降解动态频率响应模型，建立计及风电调频的电力系统频率响应模型。
步骤(2)中风电机组的功率变化量为：
其中s为拉普拉斯算子，v为实时的风速，△v(s)为风速变化量；△f(s)为频率变化量，中间参数为：
上式中Rw和kw是虚拟惯性控制的参数，ω是转子转速，kp是风能最大利用系数；kC是风能利用系数对叶尖速比的偏导数，CPref是风能利用系数的参考值，λref是叶尖速比的参考值，kβ是风能利用系数对桨距角的偏导数，kb是桨距角控制的比例系数，Hw是风力发电机的惯性时间常数。
步骤(3)中建立的计及风电调频的电力系统频率响应模型为：
其中m0，m1，m2，n0，n1，n2，n3分别是和风电机组、传统火电机组参数相关联的等值常数，△PL(s)为负荷扰动量的频域表达式。
计及风电调频的电力系统频率响应模型的时域表达式为：
其中A0、A1、B0、B1、p、ωn、ζn是常数。
本发明另一方面公开了一种含风电的电力系统频率特性指标计算方法，包括如下步骤：
(1)建立计及风电调频的电力系统频率响应模型如下：
获取时域表达式：
(2)根据计及风电调频的电力系统频率响应模型的时域表达式，计算电网动态频率特性指标，所述电网动态频率的特性指标是：初始频率变化率、最大频率偏差发生时刻、最大频率偏差、稳态频率偏差；
所述初始频率变化率IFCR计算式为：
其中△Pd为负荷扰动的时域表达式，H为火电机组的等效惯性时间常数；
所述最大频率偏差发生时刻tmax的计算为：对计及风电调频的电力系统频率响应模型时域表达式进行求导，令导数为零，此时求出的时间解即是最大频率偏差发生时刻tmax；
所述最大频率偏差△fmax的计算为：将tmax代入计及风电调频的电力系统频率响应模型时域表达式，此时求出的频率解即是最大频率偏差，△fmax＝△f(tmax)；
所述稳态频率偏差SFD的计算为：
其中△vd为负荷扰动的时域表达式，D为火电机组的等效阻尼系数，R为火电机组的等效调差系数。
有益效果：本发明公开的含风电的电力系统频率响应模型建立方法，所建立的模型能够快速、精确地获取大规模风电并网系统的动态频率响应特性和频率评估指标，有助于调度运行人员大规模风电电力系统的频率动态特性。
附图说明
图1为风电机组采用的联合频率辅助控制框图；
图2为含风电的电力系统频率响应模型框图。
具体实施方式
下面结合附图和具体实施方式，进一步阐明本发明。
含风电的电力系统频率响应模型建立方法，包括如下步骤：
(1)假设传统电力系统的发电机组类型是火电机组，将占系统总发电量比例为ρ的风电机组代替同等容量的火电机组接入电网；
(2)风电场风电机组采用虚拟惯性和桨距角减载的频率辅助控制策略，根据风电机组的转子运动方程，并采用小信号分析方法，建立风电机组的降解动态频率响应模型，包括以下步骤：
系统频率发生波动时，风电机组采用虚拟惯性和桨距角减载的频率辅助控制策略之后，风电机组的电磁转矩会发生一定的变化，具体包括两部分，一是虚拟惯性控制策略中频率偏差经过比例微分环节得到的电磁转矩变化量，如式(1)所示；二是桨距角频率控制会通过频率偏差来改变风电机组的桨距角，采用小信号分析方法，可求得因此引起的电磁转矩变化，如式(2)所示。
△Te2(s)＝2kpω△ω(s) (2)
其中s为拉普拉斯算子，△f(s)为频率变化量，Rw和kw是虚拟惯性控制的参数，ω是转子转速，△ω(s)是转速变化量，kp是风能最大利用系数。
则根据式(1)和式(2)可求得总的电磁转矩变化量可表示为:
同理，桨距角频率控制会通过频率偏差来改变机械功率，引起的机械转矩变化量为：
其中kC是风能利用系数对叶尖速比的偏导数，CPref是风能利用系数的参考值，λref是叶尖速比的参考值，kβ是风能利用系数对桨距角的偏导数，kb是桨距角控制的比例系数。如图1所示，为风电机组采用的联合频率辅助控制框图，包括了电磁转矩变化量和机械转矩变化量两部分，即式(3)和式(4)的运算。
因此，结合小信号分析方法和风电机组的转子运动方程，不平衡转矩的电力摇摆方程如下：
2Hws△ω(s)＝△Tm(s)-△Te(s) (5)
其中Hw是风力发电机的惯性时间常数。
将式(3)和式(4)代入式(5)，可得到风电机组的功率变化量，即风电机组的降解动态频率响应模型如下：
其中中间参数为：
(3)根据传统电力系统的频率响应模型和建立的风电机组降解动态频率响应模型，建立计及风电调频的电力系统频率响应模型，如下：
其中m0，m1，m2，n0，n1，n2，n3分别是和风电机组、传统火电机组参数相关联的等值常数，△PL(s)为负荷扰动量的频域表达式。如图2所示，为含风电的电力系统频率响应模型框图，描述了式(8)的运算过程。其中FH为等值高压缸做工比例，TR为等值再热时间常数，R为调差系数，H为传统发电机等值惯性时间常数，D为系统的等值阻尼系数，
根据传统电力系统的频率响应模型和建立的风电机组降解动态频率响应模型，推导出计及风电调频的电力系统频率响应模型的表达式，具体步骤为：
将式(8)的频率响应表达式进行如下分解：
则对应的时域表达式为：
其中A0、A1、B0、B1、p、ωn、ζn是常数。
根据上述步骤中建立的含风电的电力系统频率响应模型的时域表达式，可以求出对应的四个频率响应特性指标，电网动态频率的四个重要特性指标是：初始频率变化率、最大频率偏差发生时刻、最大频率偏差、稳态频率偏差。
初始频率变化率IFCR计算式为：
其中△Pd为负荷扰动的时域表达式；
最大频率偏差发生时刻tmax的计算为：对计及风电调频的电力系统频率响应模型时域表达式进行求导，令导数为零，此时求出的时间解即是最大频率偏差发生时刻tmax；
最大频率偏差△fmax的计算为：将tmax代入计及风电调频的电力系统频率响应模型时域表达式，此时求出的频率解即是最大频率偏差，△fmax＝△f(tmax)；
稳态频率偏差SFD的计算为：
其中△vd为负荷扰动的时域表达式，D为火电机组的等效阻尼系数，R为火电机组的等效调差系数。
Claims (5)
1.含风电的电力系统频率响应模型建立方法，其特征在于，包括如下步骤：
(1)假设传统电力系统的发电机组类型是火电机组，将占系统总发电量比例为ρ的风电机组代替同等容量的火电机组接入电网；
(2)风电场风电机组采用虚拟惯性和桨距角减载的频率辅助控制策略，根据风电机组的转子运动方程，并采用小信号分析方法，建立风电机组的降解动态频率响应模型；
(3)根据传统电力系统的频率响应模型和建立的风电机组降解动态频率响应模型，建立计及风电调频的电力系统频率响应模型。
2.根据权利要求1所述的含风电的电力系统频率响应模型建立方法，其特征在于，步骤(2)中风电机组的功率变化量为：
<mrow>
<msub>
<mi>&Delta;P</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>g</mi>
<mrow>
<mi>q</mi>
<mi>s</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mi>&Delta;</mi>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>as</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mi>b</mi>
<mi>s</mi>
<mo>+</mo>
<mi>c</mi>
</mrow>
<mrow>
<mi>q</mi>
<mi>s</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mi>&Delta;</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
其中s为拉普拉斯算子，v为实时的风速，Δv(s)为风速变化量；Δf(s)为频率变化量，中间参数为：
<mrow>
<mi>q</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>H</mi>
<mi>w</mi>
</msub>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mo>/</mo>
<msub>
<mi>k</mi>
<mi>p</mi>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>&omega;</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msup>
<mi>v</mi>
<mn>3</mn>
</msup>
<msub>
<mi>C</mi>
<mrow>
<mi>Pr</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>C</mi>
</msub>
<msup>
<mi>&omega;v</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>g</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>3</mn>
<msub>
<mi>k</mi>
<mi>p</mi>
</msub>
<msup>
<mi>&omega;</mi>
<mn>3</mn>
</msup>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>3</mn>
<msub>
<mi>C</mi>
<mrow>
<mi>Pr</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>C</mi>
</msub>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>&omega;</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msup>
<mi>v</mi>
<mn>3</mn>
</msup>
<msub>
<mi>C</mi>
<mrow>
<mi>Pr</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>C</mi>
</msub>
<msup>
<mi>&omega;v</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>a</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>H</mi>
<mi>w</mi>
</msub>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mo>/</mo>
<msub>
<mi>k</mi>
<mi>p</mi>
</msub>
<msub>
<mi>k</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>&omega;</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msup>
<mi>v</mi>
<mn>3</mn>
</msup>
<msub>
<mi>C</mi>
<mrow>
<mi>Pr</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>C</mi>
</msub>
<msup>
<mi>&omega;v</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>b</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>H</mi>
<mi>w</mi>
</msub>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mo>/</mo>
<msub>
<mi>R</mi>
<mi>w</mi>
</msub>
<msub>
<mi>k</mi>
<mi>p</mi>
</msub>
<mo>-</mo>
<mn>3</mn>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>k</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>&omega;</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msup>
<mi>v</mi>
<mn>3</mn>
</msup>
<msub>
<mi>C</mi>
<mrow>
<mi>Pr</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>C</mi>
</msub>
<msup>
<mi>&omega;v</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>k</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<mi>c</mi>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mn>3</mn>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>R</mi>
<mi>w</mi>
</msub>
</mfrac>
<mo>+</mo>
<msub>
<mi>k</mi>
<mi>p</mi>
</msub>
<msub>
<mi>k</mi>
<mi>&beta;</mi>
</msub>
<msub>
<mi>k</mi>
<mi>b</mi>
</msub>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>&omega;</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msup>
<mi>v</mi>
<mn>3</mn>
</msup>
<msub>
<mi>C</mi>
<mrow>
<mi>Pr</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>C</mi>
</msub>
<msup>
<mi>&omega;v</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>R</mi>
<mi>w</mi>
</msub>
</mfrac>
</mrow>
上式中Rw和kw是虚拟惯性控制的参数，ω是转子转速，kp是风能最大利用系数；kC是风能利用系数对叶尖速比的偏导数，CPref是风能利用系数的参考值，λref是叶尖速比的参考值，kβ是风能利用系数对桨距角的偏导数，kb是桨距角控制的比例系数，Hw是风力发电机的惯性时间常数。
3.根据权利要求1所述的含风电的电力系统频率响应模型建立方法，其特征在于，步骤(3)中建立的计及风电调频的电力系统频率响应模型为：
<mrow>
<mi>&Delta;</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>&Delta;P</mi>
<mi>L</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mi>g</mi>
<mrow>
<mi>q</mi>
<mi>s</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mi>&Delta;</mi>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mfrac>
<mrow>
<msub>
<mi>m</mi>
<mn>0</mn>
</msub>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>n</mi>
<mn>0</mn>
</msub>
<msup>
<mi>s</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>n</mi>
<mn>2</mn>
</msub>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>n</mi>
<mn>3</mn>
</msub>
</mrow>
</mfrac>
</mrow>
其中m0，m1，m2，n0，n1，n2，n3分别是和风电机组、传统火电机组参数相关联的等值常数，ΔPL(s)为负荷扰动量的频域表达式。
4.根据权利要求3所述的含风电的电力系统频率响应模型建立方法，其特征在于，计及风电调频的电力系统频率响应模型的时域表达式为：
其中A0、A1、B0、B1、p、ωn、ζn是常数。
5.含风电的电力系统频率特性指标计算方法，其特征在于，包括如下步骤：
(1)采用权利要求3中所述的计及风电调频的电力系统频率响应模型如下：
<mrow>
<mi>&Delta;</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>&Delta;P</mi>
<mi>L</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mi>g</mi>
<mrow>
<mi>q</mi>
<mi>s</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mi>&Delta;</mi>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mfrac>
<mrow>
<msub>
<mi>m</mi>
<mn>0</mn>
</msub>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>n</mi>
<mn>0</mn>
</msub>
<msup>
<mi>s</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>n</mi>
<mn>2</mn>
</msub>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>n</mi>
<mn>3</mn>
</msub>
</mrow>
</mfrac>
</mrow>
获取时域表达式：
(2)根据计及风电调频的电力系统频率响应模型的时域表达式，计算电网动态频率特性指标，所述电网动态频率的特性指标是：初始频率变化率、最大频率偏差发生时刻、最大频率偏差、稳态频率偏差；
所述初始频率变化率IFCR计算式为：
<mrow>
<mi>I</mi>
<mi>F</mi>
<mi>C</mi>
<mi>R</mi>
<mo>=</mo>
<munder>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>m</mi>
</mrow>
<mrow>
<mi>t</mi>
<mo>=</mo>
<msup>
<mn>0</mn>
<mo>+</mo>
</msup>
</mrow>
</munder>
<mfrac>
<mrow>
<mi>d</mi>
<mi>&Delta;</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<munder>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>m</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>=</mo>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</munder>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mi>&Delta;</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>&Delta;P</mi>
<mi>d</mi>
</msub>
<mfrac>
<msub>
<mi>m</mi>
<mn>0</mn>
</msub>
<msub>
<mi>n</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mi>d</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>&rho;k</mi>
<mi>w</mi>
</msub>
<mo>+</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>Hk</mi>
<mi>p</mi>
</msub>
</mrow>
</mfrac>
</mrow>
其中ΔPd为负荷扰动的时域表达式，H为火电机组的等效惯性时间常数；
所述最大频率偏差发生时刻tmax的计算为：对计及风电调频的电力系统频率响应模型时域表达式进行求导，令导数为零，此时求出的时间解即是最大频率偏差发生时刻tmax；
所述最大频率偏差Δfmax的计算为：将tmax代入计及风电调频的电力系统频率响应模型时域表达式，此时求出的频率解即是最大频率偏差，Δfmax＝Δf(tmax)；
所述稳态频率偏差SFD的计算为：
<mrow>
<mi>S</mi>
<mi>F</mi>
<mi>D</mi>
<mo>=</mo>
<munder>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>m</mi>
</mrow>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</munder>
<mi>&Delta;</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>m</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>=</mo>
<msup>
<mn>0</mn>
<mo>+</mo>
</msup>
</mrow>
</munder>
<mi>s</mi>
<mi>&Delta;</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&Delta;P</mi>
<mi>d</mi>
</msub>
<mo>+</mo>
<msub>
<mi>g&Delta;v</mi>
<mi>d</mi>
</msub>
<mo>)</mo>
</mrow>
<mfrac>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<msub>
<mi>n</mi>
<mn>3</mn>
</msub>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mi>d</mi>
</msub>
<mo>+</mo>
<msub>
<mi>g&Delta;v</mi>
<mi>d</mi>
</msub>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&rho;</mi>
<mo>)</mo>
<mi>D</mi>
<mo>+</mo>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&rho;</mi>
<mo>)</mo>
<mo>/</mo>
<mi>R</mi>
<mo>+</mo>
<mi>&rho;</mi>
<mi>c</mi>
</mrow>
</mfrac>
</mrow>
其中Δvd为负荷扰动的时域表达式，D为火电机组的等效阻尼系数，R为火电机组的等效调差系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711011404.0A CN107679769B (zh) | 2017-10-25 | 2017-10-25 | 含风电的电力系统频率响应模型建立方法和频率特性指标计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711011404.0A CN107679769B (zh) | 2017-10-25 | 2017-10-25 | 含风电的电力系统频率响应模型建立方法和频率特性指标计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107679769A true CN107679769A (zh) | 2018-02-09 |
| CN107679769B CN107679769B (zh) | 2020-07-14 |
Family
ID=61141735
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711011404.0A Active CN107679769B (zh) | 2017-10-25 | 2017-10-25 | 含风电的电力系统频率响应模型建立方法和频率特性指标计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107679769B (zh) |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109217292A (zh) * | 2018-08-30 | 2019-01-15 | 广东电网有限责任公司 | 一种适用于多机电力系统调频分析的融合建模方法及装置 |
| CN109560573A (zh) * | 2018-12-10 | 2019-04-02 | 国网青海省电力公司 | 一种变速风电机组频率控制器参数的优化方法及装置 |
| CN110416999A (zh) * | 2019-07-03 | 2019-11-05 | 东南大学 | 一种考虑风电参与调频的风电场聚合频率响应模型构建方法 |
| CN110750882A (zh) * | 2019-09-29 | 2020-02-04 | 武汉大学 | 一种考虑频率约束的风电占比极限值解析计算方法 |
| CN110854926A (zh) * | 2019-10-10 | 2020-02-28 | 中国电力科学研究院有限公司 | 风电接入后电力系统暂态频率响应特性分析方法及装置 |
| CN111864813A (zh) * | 2020-06-23 | 2020-10-30 | 国网辽宁省电力有限公司电力科学研究院 | 基于虚拟权重系数的风/火电联合频率控制方法 |
| CN112332432A (zh) * | 2020-12-28 | 2021-02-05 | 北京隆普智能科技有限公司 | 一种风电并网频率响应评估方法及系统 |
| CN112542853A (zh) * | 2020-11-27 | 2021-03-23 | 国网新源张家口风光储示范电站有限公司 | 基于人工智能的风电并网频率响应控制方法及系统 |
| CN112653158A (zh) * | 2020-12-07 | 2021-04-13 | 南京航空航天大学 | 一种用于含附加频率控制的双馈风电场频率响应特性辨识方法 |
| CN113054672A (zh) * | 2021-03-03 | 2021-06-29 | 山东大学 | 一种考虑多环节协同控制的风电调频方法及系统 |
| CN114123244A (zh) * | 2021-11-26 | 2022-03-01 | 南京邮电大学 | 一种考虑风-储-直联合调频的电网频率特征计算方法 |
| CN115688469A (zh) * | 2022-11-15 | 2023-02-03 | 华北电力大学 | 一种考虑节点频率变化率约束的电力系统机组组合方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150076823A1 (en) * | 2013-09-17 | 2015-03-19 | King Fahd University Of Petroleum And Minerals | Wind turbine permanent magnet synchronous generator (wt-pmsg) system |
| CN105162164A (zh) * | 2015-09-22 | 2015-12-16 | 河海大学 | 一种建立含风电并网系统的低阶动态频率响应模型的方法 |
-
2017
- 2017-10-25 CN CN201711011404.0A patent/CN107679769B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150076823A1 (en) * | 2013-09-17 | 2015-03-19 | King Fahd University Of Petroleum And Minerals | Wind turbine permanent magnet synchronous generator (wt-pmsg) system |
| CN105162164A (zh) * | 2015-09-22 | 2015-12-16 | 河海大学 | 一种建立含风电并网系统的低阶动态频率响应模型的方法 |
Non-Patent Citations (2)
| Title |
|---|
| S.GHOSH AND N.SENROY: "Electromechanical Dynamics of Controlled Variable-Speed Wind Turbines", 《IEEE SYSTEMS JOURNAL》 * |
| 汤奕 等: "基于虚拟惯量的风电场黑启动频率协同控制策略", 《电力系统自动化》 * |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109217292A (zh) * | 2018-08-30 | 2019-01-15 | 广东电网有限责任公司 | 一种适用于多机电力系统调频分析的融合建模方法及装置 |
| CN109560573B (zh) * | 2018-12-10 | 2022-07-05 | 国网青海省电力公司 | 一种变速风电机组频率控制器参数的优化方法及装置 |
| CN109560573A (zh) * | 2018-12-10 | 2019-04-02 | 国网青海省电力公司 | 一种变速风电机组频率控制器参数的优化方法及装置 |
| CN110416999A (zh) * | 2019-07-03 | 2019-11-05 | 东南大学 | 一种考虑风电参与调频的风电场聚合频率响应模型构建方法 |
| CN110750882A (zh) * | 2019-09-29 | 2020-02-04 | 武汉大学 | 一种考虑频率约束的风电占比极限值解析计算方法 |
| CN110750882B (zh) * | 2019-09-29 | 2022-09-02 | 武汉大学 | 一种考虑频率约束的风电占比极限值解析计算方法 |
| CN110854926A (zh) * | 2019-10-10 | 2020-02-28 | 中国电力科学研究院有限公司 | 风电接入后电力系统暂态频率响应特性分析方法及装置 |
| CN110854926B (zh) * | 2019-10-10 | 2022-09-20 | 中国电力科学研究院有限公司 | 风电接入后电力系统暂态频率响应特性分析方法及装置 |
| CN111864813A (zh) * | 2020-06-23 | 2020-10-30 | 国网辽宁省电力有限公司电力科学研究院 | 基于虚拟权重系数的风/火电联合频率控制方法 |
| CN111864813B (zh) * | 2020-06-23 | 2022-11-11 | 国网辽宁省电力有限公司电力科学研究院 | 基于虚拟权重系数的风/火电联合频率控制方法 |
| CN112542853A (zh) * | 2020-11-27 | 2021-03-23 | 国网新源张家口风光储示范电站有限公司 | 基于人工智能的风电并网频率响应控制方法及系统 |
| CN112653158A (zh) * | 2020-12-07 | 2021-04-13 | 南京航空航天大学 | 一种用于含附加频率控制的双馈风电场频率响应特性辨识方法 |
| CN112332432A (zh) * | 2020-12-28 | 2021-02-05 | 北京隆普智能科技有限公司 | 一种风电并网频率响应评估方法及系统 |
| CN113054672A (zh) * | 2021-03-03 | 2021-06-29 | 山东大学 | 一种考虑多环节协同控制的风电调频方法及系统 |
| CN114123244A (zh) * | 2021-11-26 | 2022-03-01 | 南京邮电大学 | 一种考虑风-储-直联合调频的电网频率特征计算方法 |
| CN114123244B (zh) * | 2021-11-26 | 2023-09-22 | 南京邮电大学 | 一种考虑风-储-直联合调频的电网频率特征计算方法 |
| CN115688469A (zh) * | 2022-11-15 | 2023-02-03 | 华北电力大学 | 一种考虑节点频率变化率约束的电力系统机组组合方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107679769B (zh) | 2020-07-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107679769A (zh) | 含风电的电力系统频率响应模型建立方法和频率特性指标计算方法 | |
| CN104967135B (zh) | 含风电的多域时滞互联电力系统滑模负荷频率控制方法 | |
| Chowdhury et al. | A novel aggregated DFIG wind farm model using mechanical torque compensating factor | |
| CN106786677A (zh) | 一种计及风电的互联电力系统分布式动态矩阵频率控制方法 | |
| Pulgar Painemal | Wind farm model for power system stability analysis | |
| CN104063584B (zh) | 一种汽轮机调速系统控制参数整定方法 | |
| CN106712058A (zh) | 双馈风机风电场参与电力系统一次调频的协调控制方法 | |
| CN105490282B (zh) | 考虑微电源无功出力均衡度的微电网实时电压控制方法 | |
| CN106026113A (zh) | 一种具有无功自动补偿的微电网系统的监控方法 | |
| CN108011364A (zh) | 一种分析dfig动力学特性与电力系统动态交互影响的方法 | |
| CN114006383B (zh) | 含vsg的双馈风机并网系统建模与次同步振荡分析的方法 | |
| CN105186502B (zh) | 基于安全域的含双馈风机电力系统暂态稳定性分析方法 | |
| CN108599264A (zh) | 一种基于虚拟同步发电机控制的频率电压无差调节方法 | |
| CN110970925A (zh) | 基于双馈风机快速有功功率调制改善系统阻尼及建模方法 | |
| CN106786759A (zh) | 与同步发电机组相协调的风电机组一次调频方法 | |
| Morel et al. | A robust control approach for primary frequency regulation through variable speed wind turbines | |
| CN106294959B (zh) | 模型参考自适应控制与双馈风机降阶模型相结合的建模仿真方法 | |
| Gao et al. | A fast high-precision model of the doubly-fed pumped storage unit | |
| Zhang et al. | Enhanced wind turbine maximum wind-energy capture based on the inverse-system method | |
| CN102801180B (zh) | 一种微电网并网等效模型的参数辨识方法 | |
| CN110416999A (zh) | 一种考虑风电参与调频的风电场聚合频率响应模型构建方法 | |
| CN110336299B (zh) | 一种考虑综合能源系统小干扰稳定的配网重构方法 | |
| CN103956767B (zh) | 一种考虑尾流效应的风电场并网稳定性分析方法 | |
| Belaimeche et al. | A comparative study between a simplified fuzzy PI and classic PI input-output linearizing controller for the wind-turbine doubly fed induction generator | |
| CN107134796A (zh) | 一种基于改进自抗扰的光火打捆次同步振荡抑制方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |