CN111384730B - 一种风机虚拟惯量控制参数的确定方法 - Google Patents
一种风机虚拟惯量控制参数的确定方法 Download PDFInfo
- Publication number
- CN111384730B CN111384730B CN202010227313.6A CN202010227313A CN111384730B CN 111384730 B CN111384730 B CN 111384730B CN 202010227313 A CN202010227313 A CN 202010227313A CN 111384730 B CN111384730 B CN 111384730B
- Authority
- CN
- China
- Prior art keywords
- virtual inertia
- fan
- balance
- parameters
- frequency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/48—Controlling the sharing of the in-phase component
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Eletrric Generators (AREA)
- Wind Motors (AREA)
Abstract
本发明实施例公开了一种风机虚拟惯量控制参数的确定方法,包括基于虚拟惯量控制结束后,***能够运行于平衡点的考虑,构建虚拟惯量参数的平衡约束;根据***的状态方程构建非线性规划,确定所述虚拟惯量参数的最优解。本发明针对虚拟惯量控制参数缺乏依据的缺点,提出了一种新的参数确定方法。推导关于虚拟惯量参数Kp、Kd的平衡约束,确保在虚拟惯量控制结束后,***能够运行于平衡点;根据***的状态方程构建非线性规划问题,确定Kp、Kd的最优值,在确保风机稳定运行的同时,实现调频效果的最优化。
Description
技术领域
本发明涉及机组调频技术领域,具体地说是一种风机虚拟惯量控制参数的确定方法。
背景技术
风机通过换流装置接入电力***,因此其输出功率与***频率解耦,不能响应***的频率变化。如今风电接入电力***规模的日益扩大,给***频率稳定性带来挑战。虚拟惯量控制通过模拟同步机组的特性,使风机的输出功率随***频率的变化而变化,从而改善***频率特性。
目前,虚拟惯量控制的参数主要依据经验或仿真结果确定,缺乏理论依据,风机有可能在虚拟惯量控制的过程中达到转速下限,退出虚拟惯量控制,引起频率二次跌落,难以保证风机运行的稳定性和频率调整效果;且依赖仿真结果的方法运算量较大,当***参数改变时需要重新进行大量计算。
发明内容
本发明实施例中提供了一种风机虚拟惯量控制参数的确定方法,以解决现有技术中依赖经验或仿真结果确定的虚拟惯量控制参数难以保证风机运行稳定性和频率调整效果的问题。
为了解决上述技术问题,本发明实施例公开了如下技术方案:
本发明提供了一种风机虚拟惯量控制参数的确定方法,所述方法包括以下步骤:
基于虚拟惯量控制结束后,***能够运行于平衡点的考虑,构建虚拟惯量参数的平衡约束;
根据***的状态方程构建非线性规划,确定所述虚拟惯量参数的最优解。
进一步地,所述构建虚拟惯量参数的平衡约束的具体过程为:
采用虚拟惯量控制时的风机模型以及***的频率响应模型,确定***的状态方程;
在***达到平衡时,各状态变量对时间导数为零,得到关于风机转速的二次多项式;
在满足虚拟惯量控制过程中达到平衡点,所述二次多项式至少一个根大于风机运行转速下限的条件下,得到虚拟惯量参数满足的平衡约束。
进一步地,所述风机模型为:
(1)式中,Cp(λ,β)由下式确定:
其中,Pew为风机输出的电磁功率;Kp、Kd为虚拟惯量控制的参数;Δf为***频率偏差量;k为风机最大功率跟踪曲线的系数;ω为风机转速;Pmw为风机捕获的机械功率;ρ为空气密度;A为风力机扫过的面积;Cp为风机的风能转换系数;λ为叶尖速比;λi为中间变量;r为风力机叶片半径;β为桨距角,Vw为风速;Hw为风机的惯性时间常数。
进一步地,所述***的频率响应模型为:
(3)式中,R为同步机调节器调差系数;ΔPms为同步机机械功率的偏差量;Km为与发电机功率因数和备用系数相关的系数;FH为原动机高压缸做功比例;TR为原动机再热时间常数;ΔPes为同步机电磁功率的偏差量;PL为负荷功率在初始时刻的变化量,由扰动初始时刻的频率变化率计算得出;Hs为同步机的惯性时间常数。
进一步地,所述***状态方程为:
所述在***达到平衡时,各状态变量对时间导数为零,则有:
对ω在扰动前的运行点ω=ω0处进行二阶近似,得到关于ω的二次多项式:
Mω2+Nω+Q=0 (7)
(7)式中,
进一步地,所述虚拟惯量参数满足的平衡约束为:
进一步地,所述***的状态方程为:
对(9)式在扰动前的***运行点处进行线性化,得到:
欲使虚拟惯量控制达到最佳调频效果,Kp、Kd应使Δfmin最小,考虑式(8)中的平衡约束,则可构建以下非线性规划问题:
利用内点法求解所述非线性规划问题,得到Kp、Kd的最优解Kp *、Kd *,即为最终确定的虚拟惯量参数值。
发明内容中提供的效果仅仅是实施例的效果,而不是发明所有的全部效果,上述技术方案中的一个技术方案具有如下优点或有益效果:
本发明针对虚拟惯量控制参数缺乏依据的缺点,提出了一种新的参数确定方法。推导关于虚拟惯量参数Kp、Kd的平衡约束,确保在虚拟惯量控制结束后,***能够运行于平衡点;根据***的状态方程构建非线性规划问题,确定Kp、Kd的最优值,在确保风机稳定运行的同时,实现调频效果的最优化。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,对于本领域普通技术人员而言,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明所述方法的流程示意图。
具体实施方式
为能清楚说明本方案的技术特点,下面通过具体实施方式,并结合其附图,对本发明进行详细阐述。下文的公开提供了许多不同的实施例或例子用来实现本发明的不同结构。为了简化本发明的公开,下文中对特定例子的部件和设置进行描述。此外,本发明可以在不同例子中重复参考数字和/或字母。这种重复是为了简化和清楚的目的,其本身不指示所讨论各种实施例和/或设置之间的关系。应当注意,在附图中所图示的部件不一定按比例绘制。本发明省略了对公知组件和处理技术及工艺的描述以避免不必要地限制本发明。
如图1所示,本发明一种风机虚拟惯量控制参数的确定方法,方法包括以下步骤:
S1,基于虚拟惯量控制结束后,***能够运行于平衡点的考虑,构建虚拟惯量参数的平衡约束;
S2,根据***的状态方程构建非线性规划,确定所述虚拟惯量参数的最优解。
步骤S1中,构建虚拟惯量参数的平衡约束的具体过程为:
采用虚拟惯量控制时的风机模型以及***的频率响应模型,确定***的状态方程;
在***达到平衡时,各状态变量对时间导数为零,得到关于风机转速的二次多项式;
在满足虚拟惯量控制过程中达到平衡点,所述二次多项式至少一个根大于风机运行转速下限的条件下,得到虚拟惯量参数满足的平衡约束。
为了保证在虚拟惯量的控制过程中,风机转速能够自行达到平衡状态且不低于转速下限,以避免退出虚拟惯量控制引起的频率二次跌落,需要对Kp、Kd做出约束。本步骤推导出了关于Kp、Kd的平衡约束。
采用虚拟惯量控制时,所述风机模型为:
(1)式中,Cp(λ,β)由下式确定:
其中,Pew为风机输出的电磁功率;Kp、Kd为虚拟惯量控制的参数;Δf为***频率偏差量;k为风机最大功率跟踪曲线的系数;ω为风机转速;Pmw为风机捕获的机械功率;ρ为空气密度;A为风力机扫过的面积;Cp为风机的风能转换系数;λ为叶尖速比;λi为中间变量;r为风力机叶片半径;β为桨距角,Vw为风速;Hw为风机的惯性时间常数。
***的频率响应模型为:
(3)式中,R为同步机调节器调差系数;ΔPms为同步机机械功率的偏差量;Km为与发电机功率因数和备用系数相关的系数;FH为原动机高压缸做功比例;TR为原动机再热时间常数;ΔPes为同步机电磁功率的偏差量;PL为负荷功率在初始时刻的变化量,由扰动初始时刻的频率变化率计算得出;Hs为同步机的惯性时间常数。当***中有多台同步发电机时,上述参数是多台机的等效参数。
将(1)(2)两式联立并整理,得到***状态方程为:
在***达到平衡时,各状态变量对时间导数为零,则有:
为简化计算,对ω在扰动前的运行点ω=ω0处进行二阶近似,
对(6)式进行整理,得到关于ω的二次多项式:
Mω2+Nω+Q=0 (7)
(7)式中,
欲使风机在虚拟惯量控制的过程中能够达到平衡点,则方程(7)必须是可解的,且至少有一个根大于风机运行的转速下限0.7,得到虚拟惯量参数满足的平衡约束为:
确保风机在虚拟惯量控制结束时运行在平衡点。通过化简***平衡点的表达式,推导出平衡约束,能够有效防止风机转速失去平衡。
步骤S2中,对***频率最低点进行优化。通过将***的状态方程线性化,将虚拟惯量控制的参数优化问题转化为非线性规划问题,优化频率极小值,实现最优的调频效果。
关于Kp、Kd的平衡约束保证了风电机组的稳定运行,除此之外,Kp、Kd还应使***频率响应尽可能地优化。本步骤通过构建非线性优化模型,实现了对频率极值点的优化。***的状态方程为:
对(9)式在扰动前的***运行点处进行线性化,得到:
求解以上(10)微分方程组,可得:
欲使虚拟惯量控制达到最佳调频效果,Kp、Kd应使Δfmin最小,考虑式(8)中的平衡约束,则可构建以下非线性规划问题:
联立(12)(13)两式,则(14)式可整理为:
利用内点求解非线性规划问题(15),得到Kp、Kd的最优解Kp *、Kd *,即为最终确定的虚拟惯量参数值。
以上所述只是本发明的优选实施方式,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也被视为本发明的保护范围。
Claims (1)
1.一种风机虚拟惯量控制参数的确定方法,其特征是,所述方法包括以下步骤:
在虚拟惯量控制结束后,***运行于平衡点,构建虚拟惯量参数的平衡约束;
根据***的状态方程构建非线性规划,确定所述虚拟惯量参数的最优解;
所述构建虚拟惯量参数的平衡约束的具体过程为:
采用虚拟惯量控制时的风机模型以及***的频率响应模型,确定***的状态方程;
在***达到平衡时,各状态变量对时间导数为零,得到关于风机转速的二次多项式;
在满足虚拟惯量控制过程中达到平衡点,所述二次多项式至少一个根大于风机运行转速下限的条件下,得到虚拟惯量参数满足的平衡约束;
所述风机模型为:
(1)式中,Cp(λ,β)由下式确定:
其中,Pew为风机输出的电磁功率;Kp、Kd为虚拟惯量控制的参数;Δf为***频率偏差量;k为风机最大功率跟踪曲线的系数;ω为风机转速;Pmw为风机捕获的机械功率;ρ为空气密度;A为风力机扫过的面积;Cp为风机的风能转换系数;λ为叶尖速比;λi为中间变量;r为风力机叶片半径;β为桨距角,Vw为风速;Hw为风机的惯性时间常数,t为时间;
所述***的频率响应模型为:
(3)式中,R为同步机调节器调差系数;Pms为同步机机械功率,ΔPms为同步机机械功率的偏差量;Km为与发电机功率因数和备用系数相关的系数;FH为原动机高压缸做功比例;TR为原动机再热时间常数;Pes为同步机电磁功率,ΔPes为同步机电磁功率的偏差量;PL为负荷功率在初始时刻的变化量,由扰动初始时刻的频率变化率计算得出;Hs为同步机的惯性时间常数;
所述***状态方程为:
所述在***达到平衡时,各状态变量对时间导数为零,则有:
对ω在扰动前的运行点ω=ω0处进行二阶近似,得到关于ω的二次多项式:
Mω2+Nω+Q=0 (7)
(7)式中,
所述虚拟惯量参数满足的平衡约束为:
所述***的状态方程为:
对(9)式在扰动前的***运行点处进行线性化,得到:
若满足虚拟惯量控制达到最佳调频效果,Kp、Kd应使Δfmin最小,基于式(8)中的平衡约束,则构建以下非线性规划问题:
利用内点法求解所述非线性规划问题,得到Kp、Kd的最优解Kp *、Kd *,为最终确定的虚拟惯量参数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010227313.6A CN111384730B (zh) | 2020-03-27 | 2020-03-27 | 一种风机虚拟惯量控制参数的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010227313.6A CN111384730B (zh) | 2020-03-27 | 2020-03-27 | 一种风机虚拟惯量控制参数的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111384730A CN111384730A (zh) | 2020-07-07 |
CN111384730B true CN111384730B (zh) | 2022-10-25 |
Family
ID=71221893
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010227313.6A Active CN111384730B (zh) | 2020-03-27 | 2020-03-27 | 一种风机虚拟惯量控制参数的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111384730B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111864769B (zh) * | 2020-07-29 | 2021-04-23 | 国网山东省电力公司日照供电公司 | 考虑风机和***频率响应特性的调频参数确定方法及*** |
CN113422376B (zh) * | 2021-06-30 | 2023-06-20 | 国网陕西省电力公司 | 一种基于频率稳定约束的事故备用、等效惯量配置方法、***、终端及可读存储介质 |
CN114884143A (zh) * | 2022-07-08 | 2022-08-09 | 西南交通大学 | 基于转子动能调节的风电机组输出功率虚拟滤波控制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3002453A1 (en) * | 2014-09-30 | 2016-04-06 | Siemens Aktiengesellschaft | Automatically setting parameter values of a wind farm controller |
CN108964032A (zh) * | 2018-07-13 | 2018-12-07 | 华北电力大学(保定) | 一种建立含风电的***频率响应模型的方法 |
CN109630354A (zh) * | 2018-12-06 | 2019-04-16 | 国网山东省电力公司潍坊供电公司 | 惯性控制下基于转换器控制的风机和同步机协同调频方法及*** |
CN110212574A (zh) * | 2019-06-04 | 2019-09-06 | 重庆大学 | 考虑虚拟惯量的风电控制参数协调设置方法 |
CN110417046A (zh) * | 2019-06-04 | 2019-11-05 | 重庆大学 | 面向小干扰稳定提升的风电虚拟惯量优化配置方法 |
CN110890765A (zh) * | 2019-11-19 | 2020-03-17 | 山东大学 | 双馈风机虚拟惯量调频的动态转速保护方法及*** |
-
2020
- 2020-03-27 CN CN202010227313.6A patent/CN111384730B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3002453A1 (en) * | 2014-09-30 | 2016-04-06 | Siemens Aktiengesellschaft | Automatically setting parameter values of a wind farm controller |
CN108964032A (zh) * | 2018-07-13 | 2018-12-07 | 华北电力大学(保定) | 一种建立含风电的***频率响应模型的方法 |
CN109630354A (zh) * | 2018-12-06 | 2019-04-16 | 国网山东省电力公司潍坊供电公司 | 惯性控制下基于转换器控制的风机和同步机协同调频方法及*** |
CN110212574A (zh) * | 2019-06-04 | 2019-09-06 | 重庆大学 | 考虑虚拟惯量的风电控制参数协调设置方法 |
CN110417046A (zh) * | 2019-06-04 | 2019-11-05 | 重庆大学 | 面向小干扰稳定提升的风电虚拟惯量优化配置方法 |
CN110890765A (zh) * | 2019-11-19 | 2020-03-17 | 山东大学 | 双馈风机虚拟惯量调频的动态转速保护方法及*** |
Non-Patent Citations (1)
Title |
---|
Control Strategy to Mitigate Secondary Frequency Dips for DFIG with Virtual Inertial Control;Zhifan Liu 等;《2016 China International Conference on Electricity Distribution (CICED)》;20160926;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111384730A (zh) | 2020-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109217362B (zh) | 一种双馈风机并网***低频振荡扰动源定位***及方法 | |
CN111384730B (zh) | 一种风机虚拟惯量控制参数的确定方法 | |
Ghosh et al. | Doubly fed induction generator (DFIG)-based wind farm control framework for primary frequency and inertial response application | |
BOUDJEMA et al. | Second order sliding mode control of a dual-rotor wind turbine system by employing a matrix converter | |
Qiao | Dynamic modeling and control of doubly fed induction generators driven by wind turbines | |
Leithead et al. | Control of variable speed wind turbines: Design task | |
CN105048917B (zh) | 基于eso的双馈风力发电***积分滑模控制器的控制方法 | |
CN103244349B (zh) | 风机塔架振动抑制***和提高风机切出风速的控制*** | |
CN106295001B (zh) | 适用于电力***中长时间尺度的准稳态变步长仿真方法 | |
Natarajan | Robust PID controller design for hydroturbines | |
CN105591395A (zh) | 一种双馈风机虚拟惯量控制后的转速恢复方法 | |
CN106786673B (zh) | 双馈风机串补输电***次同步谐振的抑制方法及装置 | |
CN105794067A (zh) | 具有改进的上升时间的风力发电站 | |
CN110829487A (zh) | 一种电力***的频率动态预测方法 | |
CN106227949B (zh) | 基于转速控制的风电机组一次调频***建模方法 | |
CN106058922A (zh) | 一种含虚拟惯量控制的双馈风电机组降阶仿真***及方法 | |
CN106227950A (zh) | 一种基于桨距控制的风电机组一次调频***动态响应建模方法 | |
Hansen et al. | Initialisation of grid-connected wind turbine models in power-system simulations | |
Krpan et al. | Towards the new low-order system frequency response model of power systems with high penetration of variable-speed wind turbine generators | |
Krpan et al. | The mathematical model of a wind power plant and a gas power plant | |
CN110417047B (zh) | 基于复转矩系数分析双馈风机ssci阻尼特性的方法 | |
Mahat et al. | Gas turbine control for islanding operation of distribution systems | |
Manjeera et al. | Design and Implementation of Fuzzy logic-2DOF controller for Emulation of wind turbine System | |
CN113708367B (zh) | 一种基于一致性算法的电力***分布式协同控制方法 | |
Riquelme et al. | A review of limitations of wind synthetic inertia methods |
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 |