CN114204611A - 一种适用所有阻尼状态的频率响应解析计算方法 - Google Patents
一种适用所有阻尼状态的频率响应解析计算方法 Download PDFInfo
- Publication number
- CN114204611A CN114204611A CN202111492815.2A CN202111492815A CN114204611A CN 114204611 A CN114204611 A CN 114204611A CN 202111492815 A CN202111492815 A CN 202111492815A CN 114204611 A CN114204611 A CN 114204611A
- Authority
- CN
- China
- Prior art keywords
- frequency response
- inertia
- coefficient
- damping
- synchronous machine
- 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
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- 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/24—Arrangements for preventing or reducing oscillations of power in networks
- H02J3/241—The oscillation concerning frequency
-
- 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/24—Arrangements for preventing or reducing oscillations of power in networks
- H02J3/242—Arrangements for preventing or reducing oscillations of power in networks using phasor measuring units [PMU]
-
- 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/36—Arrangements for transfer of electric power between ac networks via a high-tension dc link
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
-
- 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
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/22—Flexible AC transmission systems [FACTS] or power factor or reactive power compensating or correcting units
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明公开了属于电力系统调频技术领域的一种适用所有阻尼状态的频率响应解析计算方法。该方法包括:步骤1:获取全网同步机的关键调频参数;步骤2:整定机械增益系数,利用整定后的机械增益系数整定惯量、调差系数、再热时间常数和高压涡轮功率分数;步骤3:对步骤2中的四个关键调频参数进行聚合得到单机等值模型,再根据单机等值模型化简后的结构图得到惯性中心频率响应传递函数;步骤4:对步骤3中的惯性中心频率响应传递函数进行拉普拉斯反变换得到全网惯性中心频率响应时域解析模型。本发明所提供的方法能够更准确地解析计算全网惯性中心的频率响应,同时使得其计算效果进一步贴近真实电网。
Description
技术领域
本发明涉及电力系统调频技术领域,尤其涉及一种适用所有阻尼状态的频率响应解析计算方法。
背景技术
随着化石能源的逐渐开采,对能源的需求逐渐增加,人类与大自然的关系也日益紧张。在此背景下,构建以新能源为主体的新型电力系统,考虑新能源并网及其相关技术已成为未来电力技术发展的重要方向。
我国的新能源占比逐年增加,但新能源本身具有波动性以及弱惯量的特性,导致在区域电网发生负荷扰动之后,新能源对系统频率的支撑效果不佳,因此同步机依然承担主要的调频作用。
虽然目前同步机系统的频率响应研究已经比较深入,但是主流的频率响应计算方法无法适用于过阻尼状态,针对任意阻尼状态的频率响应计算方法还有待完善,如何更好地为电网制定抗扰动能力强的日前调度计划服务还有待研究。
针对上述问题,亟需一种适用所有阻尼状态的频率响应解析计算方法。已知全网同步机的所有调频参数包括惯量、调差系数、再热时间常数等;根据同步机的装机容量对调频参数进行等效聚合得到单机等值模型;在单机等值模型的基础上,通过结构图的化简得到惯性中心频率的传递函数;对惯性中心频率的传递函数进行拉普拉斯反变换得到全网惯性中心的频率响应时域解析模型。
发明内容
本发明的目的是提出一种适用所有阻尼状态的频率响应解析计算方法,其特征在于,包括以下步骤:
步骤1:通过PMU装置量测计算获取全网同步机的关键调频参数,包括惯量、调差系数、再热时间常数、高压涡轮功率分数、阻尼系数、机械增益系数和扰动功率;
步骤2:根据同步机的装机容量和系统总容量整定机械增益系数,利用整定后的机械增益系数整定惯量、调差系数、再热时间常数和高压涡轮功率分数;
步骤3:对步骤2中的四个关键调频参数进行聚合得到单机等值模型,再根据单机等值模型化简后的结构图得到惯性中心频率响应传递函数;
步骤4:对步骤3中的惯性中心频率响应传递函数进行拉普拉斯反变换得到全网惯性中心频率响应时域解析模型。
所述步骤2中整定后的同步机的机械增益系数为:
其中Sg为第g台同步机的装机容量,n为同步机的台数。
所述单机等值模型的等效调频参数为:
其中,R、FH、TR、H分别为单机等值模型的调差系数、高压涡轮功率分数、再热时间常数、惯量,Rg、FHg、TRg、Hg分别为第g台同步机的调差系数、高压涡轮功率分数、再热时间常数、惯量,λg为辅助计算的中间变量。
所述步骤3中的惯性中心频率响应传递函数为:
其中,Km、PStep分别为机械增益系数、扰动功率,D为同步机阻尼系数,ξ为阻尼系数,ωn为二阶模型的自然频率。
所述全网惯性中心频率响应时域解析模型为:
其中:
本发明的有益效果在于:
1、本发明方法能够在已知电网关键调频参数的基础上,根据同步机的装机容量和系统总容量进行多机聚合;
2、本发明能够根据单机等值模型计算得到全网惯性中心的频率响应时域解析模型,该模型可以更加精准地计算过阻尼状态的电网频率响应。
附图说明
图1为本发明适用所有阻尼状态的频率响应解析计算方法的流程图;
图2为验证算法正确性搭建的3机9节点拓扑图;
图3为现有的惯性中心频率时域解析模型SFR所计算的得到的过阻尼、欠阻尼系统扰动后频率响应图;
图4为本发明惯性中心频率时域解析模型所计算的得到的过阻尼、欠阻尼系统扰动后频率响应图。
具体实施方式
本发明提出一种适用所有阻尼状态的频率响应解析计算方法,下面结合附图和具体实施例对本发明做进一步说明。
图1为本发明适用所有阻尼状态的频率响应解析计算方法的流程图,该方法具体步骤如下:
步骤1:关键调频参数获取;
根据具体的运行要求,由人机界面获取关键调频参数如:惯量H、调差系数R、再热时间常数TR、高压涡轮功率分数FH、阻尼系数ξ、机械增益系数Km、扰动功率PStep,完毕后进入步骤2;
步骤2:多机系统聚合;
通过同步机的装机容量和系统总容量整定机械增益系数,整定后各同步机的机械增益系数由下式计算:
其中Sg为第g台同步机的装机容量,n为同步机的台数。
根据整定后的同步机的机械增益系数整定惯量、调差系数、再热时间常数、高压涡轮功率分数,最后对调频参数进行聚合得到单机等值模型,单机等值模型的等效调频参数由下式计算:
其中Rg、FHg、TRg、Hg分别为第g台同步机的调差系数、高压涡轮功率分数、再热时间常数、惯量,λg为辅助计算的中间变量,完毕后进入步骤3;
步骤3:结构图化简,计算惯性中心频率响应传递函数;
根据单机等值模型把结构图进行化简,得到惯性中心频率响应传递函数,频率响应单机等值模型结构图如图1所示,可以得到惯性中心频率响应传递函数如下式:
其中R、FH、TR、H、Km、PStep分别为单机等值模型的调差系数、高压涡轮功率分数、再热时间常数、惯量、机械增益系数、扰动功率,完毕后进入步骤4;
步骤4:对频率响应时域解析模型进行推导;
对惯性中心频率的传递函数进行拉普拉斯反变换得到全网惯性中心的频率响应时域解析模型,其时域解析模型推导过程如下所示:
其中,ωn为二阶模型的自然频率,D为同步机阻尼系数。
根据双曲正弦、双曲余弦的定义,即
上式可以化简为下式:
将上式提取公因式并合并同类项后,得到全网惯性中心的频率响应时域解析模型如下式:
其中:
图2为3机9节点系统拓扑图,该地区系统电压等级为220kV,发电机总体容量为402MW,负荷总量350MW,该系统关键调频参数聚合后如下所示:
将关键调频参数带入已有的方法得到系统在扰动后的频率响应如图3所示,
带入本发明所提的全网惯性中心的频率响应时域解析模型得到系统在扰动后的频率响应如图4所示,经过对比发现,本发明模型可以适用于计算任意阻尼状态的频率响应。
此实施例仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (5)
1.一种适用所有阻尼状态的频率响应解析计算方法,其特征在于,包括以下步骤:
步骤1:通过PMU装置量测计算获取全网同步机的关键调频参数,包括惯量、调差系数、再热时间常数、高压涡轮功率分数、阻尼系数、机械增益系数和扰动功率;
步骤2:根据同步机的装机容量和系统总容量整定机械增益系数,利用整定后的机械增益系数整定惯量、调差系数、再热时间常数和高压涡轮功率分数;
步骤3:对步骤2中的四个关键调频参数进行聚合得到单机等值模型,再根据单机等值模型化简后的结构图得到惯性中心频率响应传递函数;
步骤4:对步骤3中的惯性中心频率响应传递函数进行拉普拉斯反变换得到全网惯性中心频率响应时域解析模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111492815.2A CN114204611B (zh) | 2021-12-08 | 2021-12-08 | 一种适用所有阻尼状态的频率响应解析计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111492815.2A CN114204611B (zh) | 2021-12-08 | 2021-12-08 | 一种适用所有阻尼状态的频率响应解析计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114204611A true CN114204611A (zh) | 2022-03-18 |
CN114204611B CN114204611B (zh) | 2023-04-07 |
Family
ID=80651330
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111492815.2A Active CN114204611B (zh) | 2021-12-08 | 2021-12-08 | 一种适用所有阻尼状态的频率响应解析计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114204611B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115017449A (zh) * | 2022-08-09 | 2022-09-06 | 中国电力科学研究院有限公司 | 一种适用二阶系统不同阻尼比的频率偏差计算方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120191224A1 (en) * | 2011-01-25 | 2012-07-26 | Young Hyun Moon | Output feedback frequency control device for rotating machine speed control |
CN108964032A (zh) * | 2018-07-13 | 2018-12-07 | 华北电力大学(保定) | 一种建立含风电的系统频率响应模型的方法 |
CN110750882A (zh) * | 2019-09-29 | 2020-02-04 | 武汉大学 | 一种考虑频率约束的风电占比极限值解析计算方法 |
CN110875600A (zh) * | 2019-11-27 | 2020-03-10 | 大连理工大学 | 一种两机等值电力系统动态频率响应近似解析模型 |
CN113113908A (zh) * | 2021-03-17 | 2021-07-13 | 国网江苏省电力有限公司 | 适用于现代大电网频率响应的时域解析方法与系统 |
-
2021
- 2021-12-08 CN CN202111492815.2A patent/CN114204611B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120191224A1 (en) * | 2011-01-25 | 2012-07-26 | Young Hyun Moon | Output feedback frequency control device for rotating machine speed control |
CN108964032A (zh) * | 2018-07-13 | 2018-12-07 | 华北电力大学(保定) | 一种建立含风电的系统频率响应模型的方法 |
CN110750882A (zh) * | 2019-09-29 | 2020-02-04 | 武汉大学 | 一种考虑频率约束的风电占比极限值解析计算方法 |
CN110875600A (zh) * | 2019-11-27 | 2020-03-10 | 大连理工大学 | 一种两机等值电力系统动态频率响应近似解析模型 |
CN113113908A (zh) * | 2021-03-17 | 2021-07-13 | 国网江苏省电力有限公司 | 适用于现代大电网频率响应的时域解析方法与系统 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115017449A (zh) * | 2022-08-09 | 2022-09-06 | 中国电力科学研究院有限公司 | 一种适用二阶系统不同阻尼比的频率偏差计算方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN114204611B (zh) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102377180B (zh) | 基于电能质量监测系统的电力系统负荷建模方法 | |
CN102801158B (zh) | 基于Pade近似的时滞电力系统特征值计算与稳定性判别方法 | |
CN106443246B (zh) | 基于pmu量测数据的小干扰稳定参数的在线辨识方法 | |
CN103699723B (zh) | 一种发电厂机组动力系统模型校核方法 | |
CN103258103B (zh) | 基于偏最小二乘回归的戴维南等值参数辨识方法 | |
CN101969199A (zh) | 用于暂态功角稳定风险评估的故障损失估计方法 | |
CN108667005B (zh) | 一种计及新能源影响的电网静动态结合脆弱性评估方法 | |
CN106169115B (zh) | 一种利用解析-数值法估算两区域电网最低频率的方法 | |
CN115085287B (zh) | 一种考虑调速器限幅环节频率稳定约束的建模方法 | |
CN112510707A (zh) | 一种配电台区电力阻抗拓扑图生成方法及系统 | |
CN106786529B (zh) | 一种分布式静态安全分析方法 | |
CN106058922A (zh) | 一种含虚拟惯量控制的双馈风电机组降阶仿真系统及方法 | |
CN114204611B (zh) | 一种适用所有阻尼状态的频率响应解析计算方法 | |
CN114123344A (zh) | 基于自适应递推最小二乘的电力系统惯量评估方法及装置 | |
CN115622053A (zh) | 一种用于考虑分布式电源的自动负荷建模方法及装置 | |
CN106816879B (zh) | 一种基于紧急需求响应的主动配电网电压协调控制方法 | |
CN112072651B (zh) | 基于暂态能量函数的暂态电压失稳与功角失稳识别方法 | |
CN114398761B (zh) | 一种考虑调速器强化型死区的同步机系统频率响应解析计算方法 | |
CN114243805B (zh) | 一种考虑调速器限幅的同步机系统频率响应解析计算方法 | |
CN105701265B (zh) | 一种双馈风电机组建模方法及装置 | |
CN105305392A (zh) | 适用于含电压控型iidg配电网的短路计算对称分量法 | |
CN102545218A (zh) | 基于电能质量监测系统的在线负荷建模并行计算方法 | |
CN104898415A (zh) | 一种基于低通滤波的发电机组转动惯量参数的在线辨识方法 | |
CN114221396B (zh) | 一种考虑调速器一般型死区的频率响应解析计算方法 | |
CN112688321A (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 |