CN114611987A - 一种并网型风电场的局域振荡的辨识方法 - Google Patents

一种并网型风电场的局域振荡的辨识方法 Download PDF

Info

Publication number
CN114611987A
CN114611987A CN202210307371.9A CN202210307371A CN114611987A CN 114611987 A CN114611987 A CN 114611987A CN 202210307371 A CN202210307371 A CN 202210307371A CN 114611987 A CN114611987 A CN 114611987A
Authority
CN
China
Prior art keywords
oscillation
power plant
wind power
node
grid
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
Application number
CN202210307371.9A
Other languages
English (en)
Other versions
CN114611987B (zh
Inventor
桑顺
朱悦
薛晓岑
黄杰杰
张雷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nantong University
Original Assignee
Nantong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nantong University filed Critical Nantong University
Priority to CN202210307371.9A priority Critical patent/CN114611987B/zh
Publication of CN114611987A publication Critical patent/CN114611987A/zh
Application granted granted Critical
Publication of CN114611987B publication Critical patent/CN114611987B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • 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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy 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
    • 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)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明提供一种并网型风电场的局域振荡的辨识方法，将并网型风电场系统中的节点划分为弱风险节点和高风险节点，给出了并网型风电场系统的振荡局域传播强度指标的计算方法，根据振荡局域传播强度指标判断振荡是局域振荡或广域振荡。所提出的并网型风电场局域振荡的辨识方法能够有效地辨识并网型风电场的局域振荡、广域振荡，将并网型风电场划分为振荡区域和非振荡区域，便于有针对地部署振荡抑制措施。相比仅依靠能观性指标刻画局域/广域振荡，所提方法能够提高局域/广域振荡的辨识精度。

Description

一种并网型风电场的局域振荡的辨识方法
技术领域
本发明涉及风力发电场的电气控制技术领域，尤其涉及一种并网型风电场的局域振荡的辨识方法。
背景技术
对于电力电子变换器主导的风力发电场，系统内部元件不仅通过电气端口进行连接，还依赖于电力电子变换器实现交直流侧互联。已经有研究表明，交流/直流变换器通过控制可以实现交-直流动态解耦，即体现出“振荡防火墙”特性。该控制特性预期会对振荡的传播特性产生影响，如：所产生的振荡可能被阻断或限制在某一区域内，从而产生局域振荡现象。
局域/广域振荡本质上反映的是振荡能否可观，因此直观上可采用能观性指标进行刻画。但是，实际应用中发现，仅采用能观性指标刻画局域/广域振荡存在歧义，容易出现误判，需要提出一种精确判断并网型风电场局域/广域振荡的指标及方法。
发明内容
本发明的目的是为了解决如何辨识并网型风电场局域/广域振荡现象的难题，提出振荡局域传播强度指标。
为了实现上述目的，本发明采用了如下技术方案：
一种并网型风电场的局域振荡的辨识方法，包括以下步骤：
S1：在旋转dq坐标系下，建立含N个节点的并网型风电场的节点导纳矩阵Ynode(s)，令节点导纳矩阵对应的行列式为零，求出系统的闭环特征值；根据特征值实部判断风电场的稳定，若存在实部为正值的特征值，表示风电场运行失稳，风电场的振荡频率ω为实部为正值的特征值的虚部；
S2：将振荡频率ω代入并网型风电场的节点导纳矩阵，求出含N个节点的并网型风电场系统中各个节点的能观度Obs1、Obs2、……、ObsN；
S3：对并网型风电场系统中的节点进行划分，将各个节点的能观度Obs1、Obs2、……、ObsN与阈值0.00001进行比较，若某节点的能观度Obsn(1≤n≤N)小于等于阈值0.00001，则该节点n为弱风险节点，否则，该节点为高风险节点；
S4：构建置换矩阵Pπ，若并网型风电场系统中含有m个弱风险节点、N-m个高风险节点，这m个弱风险节点的节点号分别为π(1)、π(2)、……、π(m)，N-m个高风险节点的节点号分别为π(m)、π(m+1)、……、π(N)，那么含N个节点的并网型风电场系统的置换π的映射关系为
Figure BDA0003566105730000021
置换矩阵Pπ的表达式为
Figure BDA0003566105730000022
其中eπ(j)为标准基，代表长度为N的行向量，其中仅有第π(j)列元素为1，而其余元素均为0；
S5：对节点导纳矩阵进行置换运算如下：
Figure BDA0003566105730000031
其中
Figure BDA0003566105730000032
为节点导纳矩阵在振荡频率ω处的逆矩阵，
Figure BDA0003566105730000033
为置换矩阵Pπ的转置矩阵，
Figure BDA0003566105730000034
为经过重新排布后在振荡频率ω处的节点阻抗矩阵，ZR-L为经过重新排布后的弱风险区域矩阵，ZR-H为经过重新排布后的高风险区域矩阵，Zoff-P、Zoff-N为经过重新排布后弱风险、高风险区域之间的交互耦合项；
S6：定义频率为ω的振荡局域传播强度指标kP(ω)，其表达式为：
Figure BDA0003566105730000035
其中，||||2表示2范数，即欧几里得范数；
S7：根据振荡局域传播强度指标kP(ω)对局域振荡进行量化评估，若kP(ω)＝0，系统呈现局域振荡特性；若kP(ω)≠0，系统呈现广域振荡特性。
与现有技术方案相比，本申请所提出的并网型风电场局域振荡指标及辨识方法具有如下优势：
能够有效地辨识并网型风电场的局域振荡、广域振荡，将并网型风电场划分为振荡区域和非振荡区域，便于有针对地部署振荡抑制措施；相比仅依靠能观性指标刻画局域/广域振荡，所提方法能够提高局域/广域振荡的辨识精度。
附图说明
图1为本发明一个实施例—并网型风电场局域振荡辨识方法的流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，以下结合具体实施例，对本发明作进一步地详细说明。
在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用不同于在此描述的其他方式来实施，因此，本发明并不限于下面公开说明书的具体实施例的限制。
请参阅图1，本发明一个实施例本发明一个实施例—并网型风电场局域振荡辨识方法的流程图，包括以下步骤：
S1：在旋转dq坐标系下，建立含15个节点的并网型风电场的节点导纳矩阵Ynode(s)，令节点导纳矩阵对应的行列式为零，求出系统的闭环特征值，存在一对右半平面特征值：0.0561+50.6j，可见存在实部为正值的特征值，表示风电场运行失稳，风电场的振荡频率ω为特征值的虚部，即50.6rad/s；
S2：将振荡频率50.6rad/s代入并网型风电场的节点导纳矩阵，求出含15个节点的并网型风电场系统中各个节点的能观度Obs1、Obs2、……、Obs15如下表所示：
Figure BDA0003566105730000051
S3：对并网型风电场系统中的节点进行划分，将各个节点的能观度Obs1、Obs2、……、ObsN与阈值0.00001进行比较，若某节点的能观度Obsn(1≤n≤N)小于等于阈值0.00001，则该节点n为弱风险节点，否则，该节点为高风险节点；经过比较，可知节点#1～#3，#7～#12为弱风险节点，剩余节点#4～#6、#13～#15为高风险节点；
S4：构建置换矩阵Pπ，并网型风电场系统中含有9个弱风险节点、6个高风险节点，这9个弱风险节点的节点号分别为1、2、3、7、8、9、10、11、12，6个高风险节点的节点号分别为4、5、6、13、14、15，那么含15个节点的并网型风电场系统的置换π的映射关系为
Figure BDA0003566105730000052
置换矩阵Pπ的表达式为
Figure BDA0003566105730000061
S5：对节点导纳矩阵进行置换运算如下：
Figure BDA0003566105730000062
其中
Figure BDA0003566105730000063
为节点导纳矩阵在振荡频率ω处的逆矩阵，
Figure BDA0003566105730000064
为置换矩阵Pπ的转置矩阵，
Figure BDA0003566105730000065
为经过重新排布后在振荡频率ω处的节点阻抗矩阵，ZR-L为经过重新排布后的弱风险区域矩阵，ZR-H为经过重新排布后的高风险区域矩阵，Zoff-P、Zoff-N为经过重新排布后弱风险、高风险区域之间的交互耦合项；
S6：定义频率为ω的振荡局域传播强度指标kP(ω)，其表达式为：
Figure BDA0003566105730000066
其中，||||2表示2范数，即欧几里得范数；
将频率为50.6rad/s的振荡代入，计算出振荡局域传播强度指标kP(50.6)值为0；
S7：根据振荡局域传播强度指标kP(ω)对局域振荡进行量化评估，由于振荡局域传播强度指标kP(50.6)值等于0，系统呈现局域振荡特性，说明高风险区域振荡向弱风险区域的传播通道被阻断，弱风险区域保持稳定。
综上所述，申请所提出的并网型风电场局域振荡指标及辨识方法能够有效地辨识并网型风电场的局域振荡、广域振荡，将并网型风电场划分为振荡区域和非振荡区域，便于有针对地部署振荡抑制措施；相比仅依靠能观性指标刻画局域/广域振荡，所提方法能够提高局域/广域振荡的辨识精度。
以上所述仅为本发明的实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围。

Claims (1)

1.一种并网型风电场的局域振荡的辨识方法，其特征在于，包括以下步骤：
S1：在旋转dq坐标系下，建立含N个节点的并网型风电场的节点导纳矩阵Ynode(s)，令节点导纳矩阵对应的行列式为零，求出系统的闭环特征值；根据特征值实部判断风电场的稳定，若存在实部为正值的特征值，表示风电场运行失稳，风电场的振荡频率ω为实部为正值的特征值的虚部；
S2：将振荡频率ω代入并网型风电场的节点导纳矩阵，求出含N个节点的并网型风电场系统中各个节点的能观度Obs1、Obs2、……、ObsN；
S3：对并网型风电场系统中的节点进行划分，将各个节点的能观度Obs1、Obs2、……、ObsN与阈值0.00001进行比较，若某节点的能观度Obsn(1≤n≤N)小于等于阈值0.00001，则该节点n为弱风险节点，否则，该节点为高风险节点；
S4：构建置换矩阵Pπ，若并网型风电场系统中含有m个弱风险节点、N-m个高风险节点，这m个弱风险节点的节点号分别为π(1)、π(2)、……、π(m)，N-m个高风险节点的节点号分别为π(m)、π(m+1)、……、π(N)，那么含N个节点的并网型风电场系统的置换π的映射关系为
Figure FDA0003566105720000011
置换矩阵Pπ的表达式为
Figure FDA0003566105720000021
其中eπ(j)为标准基，代表长度为N的行向量，其中仅有第π(j)列元素为1，而其余元素均为0；
S5：对节点导纳矩阵进行置换运算如下：
Figure FDA0003566105720000022
其中
Figure FDA0003566105720000023
为节点导纳矩阵在振荡频率ω处的逆矩阵，
Figure FDA0003566105720000024
为置换矩阵Pπ的转置矩阵，
Figure FDA0003566105720000025
为经过重新排布后在振荡频率ω处的节点阻抗矩阵，ZR-L为经过重新排布后的弱风险区域矩阵，ZR-H为经过重新排布后的高风险区域矩阵，Zoff-P、Zoff-N为经过重新排布后弱风险、高风险区域之间的交互耦合项；
S6：定义频率为ω的振荡局域传播强度指标kP(ω)，其表达式为：
Figure FDA0003566105720000026
其中，||||2表示2范数，即欧几里得范数；
S7：根据振荡局域传播强度指标kP(ω)对局域振荡进行量化评估，若kP(ω)＝0，系统呈现局域振荡特性；若kP(ω)≠0，系统呈现广域振荡特性。
CN202210307371.9A 2022-03-25 2022-03-25 一种并网型风电场的局域振荡的辨识方法 Active CN114611987B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210307371.9A CN114611987B (zh) 2022-03-25 2022-03-25 一种并网型风电场的局域振荡的辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210307371.9A CN114611987B (zh) 2022-03-25 2022-03-25 一种并网型风电场的局域振荡的辨识方法

Publications (2)

Publication Number Publication Date
CN114611987A true CN114611987A (zh) 2022-06-10
CN114611987B CN114611987B (zh) 2022-12-09

Family

ID=81866835

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210307371.9A Active CN114611987B (zh) 2022-03-25 2022-03-25 一种并网型风电场的局域振荡的辨识方法

Country Status (1)

Country Link
CN (1) CN114611987B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114966296A (zh) * 2022-08-01 2022-08-30 锦浪科技股份有限公司 一种振荡源定位方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20210036518A1 (en) * 2019-07-30 2021-02-04 North China Electric Power University System and method of locating oscillation sources of wind power integrated system based on energy spectrums
CN112381671A (zh) * 2020-11-30 2021-02-19 华北电力科学研究院有限责任公司 一种新能源场站并网电力系统的宽频振荡风险评估方法
CN113156247A (zh) * 2021-04-23 2021-07-23 北京建筑大学 一种电力系统低频振荡的预警方法及装置
CN113644689A (zh) * 2021-08-16 2021-11-12 清华大学 风电系统运行稳定域的构建方法、装置、电子设备及其可读存储介质

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20210036518A1 (en) * 2019-07-30 2021-02-04 North China Electric Power University System and method of locating oscillation sources of wind power integrated system based on energy spectrums
CN112381671A (zh) * 2020-11-30 2021-02-19 华北电力科学研究院有限责任公司 一种新能源场站并网电力系统的宽频振荡风险评估方法
CN113156247A (zh) * 2021-04-23 2021-07-23 北京建筑大学 一种电力系统低频振荡的预警方法及装置
CN113644689A (zh) * 2021-08-16 2021-11-12 清华大学 风电系统运行稳定域的构建方法、装置、电子设备及其可读存储介质

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
占颖等: "风电并网系统次同步振荡的频域模式分析", 《电力系统自动化》 *
高娟 等: "虚拟同步发电机的输入阻抗建模及稳定性分析", 《电网技术》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114966296A (zh) * 2022-08-01 2022-08-30 锦浪科技股份有限公司 一种振荡源定位方法及装置
CN114966296B (zh) * 2022-08-01 2022-12-02 锦浪科技股份有限公司 一种振荡源定位方法及装置

Also Published As

Publication number Publication date
CN114611987B (zh) 2022-12-09

Similar Documents

Publication Publication Date Title
Wang et al. Several sufficient conditions for harmonic source identification in power systems
CN114611987B (zh) 一种并网型风电场的局域振荡的辨识方法
CN105160058B (zh) 一种用于电网电磁暂态计算的网络边界等值方法
CN109802392B (zh) 大规模配电网潮流计算方法及装置
CN105071387B (zh) 基于复数域标幺化理论的配电网快速分解状态估计方法
Khalili et al. PMU-based decoupled state estimation for unsymmetrical power systems
Hydon Discrete point symmetries of ordinary differential equations
Heinrich et al. Optimum mesh grading for finite-difference method
Bäckdahl et al. Static axisymmetric spacetimes with prescribed multipole moments
Lu et al. Detection of Broken Rotor Bars Fault in Induction Motors by Using an Improved MUSIC and Least‐Squares Amplitude Estimation
Nie et al. An improved natural frequency based transmission line fault location method with full utilization of frequency spectrum information
CN107145732A (zh) 一种基于改进cbfm的快速求解目标单站电磁散射特性的方法
CN110138241A (zh) 一种基于雅可比理论的单相电压源的稳定性分析方法
CN116224831A (zh) 基于视在阻抗的多并网逆变器谐波谐振频率灰箱辨识方法
Weng et al. Numerical harmonic modeling of long coupled transmission lines using matrix series theory and recursive approach
Pakzad A practical method for stability analysis of linear fractional-order systems with distributed delay
Kurtović et al. Earthing grid parameters with conductor surrounded by an additional substance
Hebling et al. Sparse and orthogonal method for fast bad data processing in distribution system state estimation
Xu et al. Inter‐harmonics analysis and parameter estimation based on H2R6 window and constructing low‐interference zone
Porkar et al. Multi-port frequency-dependent network equivalent for electromagnetic transient studies
Han et al. Radial basis function and finite element method hybrid approach for three-dimensional electromagnetics problems
Hone et al. Properties of the series solution for Painlevé I
Sheng et al. Novel method for power system synthesis load identification
Rakpenthai et al. Rectangular coordinates state estimation of electrical power system using singular value decomposition
CN109409726A (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