CN110504706A - 虚拟同步发电机并网的暂态稳定性预测方法及装置 - Google Patents
虚拟同步发电机并网的暂态稳定性预测方法及装置 Download PDFInfo
- Publication number
- CN110504706A CN110504706A CN201910662468.XA CN201910662468A CN110504706A CN 110504706 A CN110504706 A CN 110504706A CN 201910662468 A CN201910662468 A CN 201910662468A CN 110504706 A CN110504706 A CN 110504706A
- Authority
- CN
- China
- Prior art keywords
- virtual
- vsg
- transient
- synchronous generator
- energy function
- 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
- 230000001052 transient effect Effects 0.000 title claims abstract description 153
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 230000006870 function Effects 0.000 claims description 74
- 238000005381 potential energy Methods 0.000 claims description 24
- 238000004590 computer program Methods 0.000 claims description 15
- 230000010354 integration Effects 0.000 claims description 9
- 238000009966 trimming Methods 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 4
- 238000010248 power generation Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000002679 ablation Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
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/24—Arrangements for preventing or reducing oscillations of power in networks
-
- 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
-
- 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/40—Synchronising a generator for connection to a network or to another generator
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明实施例提供一种虚拟同步发电机并网的暂态稳定性预测方法及装置，所述方法包括：获取虚拟同步发电机VSG的虚拟转速和虚拟功角；将所述虚拟转速和虚拟功角，输入至预设的暂态能量函数模型，输出暂态能量值以及临界切除时间；根据所述暂态能量值以及所述临界切除时间预测VSG并网的暂态稳定性。本发明实施例提供的虚拟同步发电机并网的暂态稳定性预测方法及装置，构建接有VSG的单机无穷大系统的暂态能量函数，并计及虚拟阻尼的影响，建立的暂态能量函数准确度很高，通过计算暂态能量函数对虚拟同步发电机并网的暂态稳定性进行预测，计算量小，使用方便，预测效率高。
Description
技术领域
本发明涉及分布式发电技术领域，尤其涉及一种虚拟同步发电机并网的暂态稳定性预测方法及装置。
背景技术
虚拟同步发电机(VSG)技术近年来发展迅速，将逐渐成为未来新能源的发展方向，VSG技术的应用可能对系统原有稳定性问题造成影响或者引入新的稳定性问题。VSG技术主要是利用控制手段使新能源发电模拟传统同步发电机的特性，从而使新能源发电为电网提供必要的频率和电压支撑。虚拟同步机并网运行后，会对电力系统的稳定性问题产生一定的影响。一方面，VSG模拟了同步发电机的特性，势必引入类似于同步发电机的稳定性问题；另一方面，VSG底层由电力电子器件和相应的控制实现，兼有传统逆变器的部分特点。随着新能源在电网中渗透率的不断提高，深入研究VSG并网的稳定性成为了该技术继续发展的关键。
现有技术中，采用时域仿真模拟的方法对虚拟同步发电机并网的暂态稳定性进行分析。
但是，现有技术中的方法，计算量很大，暂态稳定性的预测效率很低，不方便应用。
发明内容
本发明实施例提供一种虚拟同步发电机并网的暂态稳定性预测方法及装置，用于解决现有技术中暂态稳定性的预测效率低的技术问题。
为了解决上述技术问题，一方面，本发明实施例提供一种虚拟同步发电机并网的暂态稳定性预测方法，包括：
获取虚拟同步发电机VSG的虚拟转速和虚拟功角；
将所述虚拟转速和虚拟功角，输入至预设的暂态能量函数模型，输出暂态能量值以及临界切除时间；
根据所述暂态能量值以及所述临界切除时间预测VSG并网的暂态稳定性。
进一步地，获取所述暂态能量函数模型的具体步骤如下：
根据忽略虚拟阻尼影响的VSG系统运动方程求偏导得到初始暂态能量函数；
先对初始暂态能量函数对时间求导，然后进行时间积分，计及虚拟阻尼的影响，得到所述暂态能量函数模型。
进一步地，所述暂态能量函数模型用公式表示如下：
其中，E(δ,ω)表示暂态能量函数，J表示VSG的虚拟惯量，ω表示VSG的虚拟转速，ωb表示系统转速，VVSG表示VSG的内电势，X∑表示故障切除后系统的视在电抗，V表示无穷大节点的电压，δ表示VSG的虚拟功角，δs表示系统稳定平衡点处的功角，Tm表示VSG的参考转矩，Dp表示VSG的虚拟阻尼，λ表示微调系数，0＜λ＜1，Δδ＝δ-δs。
进一步地，所述暂态能量函数模型中包括虚拟转子动能、磁性势能、耗散势能、位置势能和虚拟阻尼消耗的能量。
另一方面，本发明实施例提供一种虚拟同步发电机并网的暂态稳定性预测装置，包括：
获取模块，用于获取虚拟同步发电机VSG的虚拟转速和虚拟功角；
计算模块，用于将所述虚拟转速和虚拟功角，输入至预设的暂态能量函数模型，输出暂态能量值以及临界切除时间；
预测模块，用于根据所述暂态能量值以及所述临界切除时间预测VSG并网的暂态稳定性。
进一步地，还包括模型构建模块，具体用于：
根据忽略虚拟阻尼影响的VSG系统运动方程求偏导得到初始暂态能量函数；
先对初始暂态能量函数对时间求导，然后进行时间积分，计及虚拟阻尼的影响，得到所述暂态能量函数模型。
进一步地，所述暂态能量函数模型用公式表示如下：
其中，E(δ,ω)表示暂态能量函数，J表示VSG的虚拟惯量，ω表示VSG的虚拟转速，ωb表示系统转速，VVSG表示VSG的内电势，X∑表示故障切除后系统的视在电抗，V表示无穷大节点的电压，δ表示VSG的虚拟功角，δs表示系统稳定平衡点处的功角，Tm表示VSG的参考转矩，Dp表示VSG的虚拟阻尼，λ表示微调系数，0＜λ＜1，Δδ＝δ-δs。
进一步地，所述暂态能量函数模型中包括虚拟转子动能、磁性势能、耗散势能、位置势能和虚拟阻尼消耗的能量。
再一方面，本发明实施例提供一种电子设备，包括：存储器、处理器，以及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时，实现上述方法的步骤。
又一方面，本发明实施例提供一种非暂态计算机可读存储介质，其上存储有计算机程序，当所述计算机程序被处理器执行时，实现上述方法的步骤。
本发明实施例提供的虚拟同步发电机并网的暂态稳定性预测方法及装置，构建接有VSG的单机无穷大系统的暂态能量函数，并计及虚拟阻尼的影响，建立的暂态能量函数准确度很高，通过计算暂态能量函数对虚拟同步发电机并网的暂态稳定性进行预测，计算量小，使用方便，预测效率高。
附图说明
图1为本发明实施例提供的虚拟同步发电机并网的暂态稳定性预测方法示意图；
图2为本发明实施例提供的接有VSG的单机无穷大系统的示意图；
图3为本发明实施例提供的虚拟同步发电机并网的暂态稳定性预测装置示意图；
图4为本发明实施例提供的电子设备的结构示意图。
具体实施方式
为了使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1为本发明实施例提供的虚拟同步发电机并网的暂态稳定性预测方法示意图，如图1所示，本发明实施例提供一种虚拟同步发电机并网的暂态稳定性预测方法，其执行主体为虚拟同步发电机并网的暂态稳定性预测装置，该方法包括：
步骤S101、获取虚拟同步发电机VSG的虚拟转速和虚拟功角。
具体来说，图2为本发明实施例提供的接有VSG的单机无穷大系统的示意图，如图2所示，首先构建出接有VSG的单机无穷大系统。在这个系统的基础上，采用首次积分法构建暂态能量函数表达式，并计及VSG虚拟阻尼的影响。
在构建好暂态能量函数模型之后，获取虚拟同步发电机VSG的虚拟转速ω和虚拟功角δ。
步骤S102、将所述虚拟转速和虚拟功角，输入至预设的暂态能量函数模型，输出暂态能量值以及临界切除时间。
具体来说，获取虚拟同步发电机VSG的虚拟转速ω和虚拟功角δ之后，将虚拟同步发电机VSG的虚拟转速ω和虚拟功角δ，输入至预先构建好的暂态能量函数模型，输出暂态能量值以及临界切除时间。
步骤S103、根据所述暂态能量值以及所述临界切除时间预测VSG并网的暂态稳定性。
具体来说，在获取暂态能量值以及临界切除时间之后，根据暂态能量值以及临界切除时间对虚拟同步发电机并网的暂态稳定性进行分析，对实际现场设计有指导作用，避免成本浪费。
本发明实施例提供的虚拟同步发电机并网的暂态稳定性预测方法，构建接有VSG的单机无穷大系统的暂态能量函数，并计及虚拟阻尼的影响，建立的暂态能量函数准确度很高，通过计算暂态能量函数对虚拟同步发电机并网的暂态稳定性进行预测，计算量小，使用方便，预测效率高。
基于上述任一实施例，进一步地，获取所述暂态能量函数模型的具体步骤如下：
根据忽略虚拟阻尼影响的VSG系统运动方程求偏导得到初始暂态能量函数；
先对初始暂态能量函数对时间求导，然后进行时间积分，计及虚拟阻尼的影响，得到所述暂态能量函数模型。
具体来说，构建暂态能量函数模型的具体步骤如下：
基于图2中建出的接有VSG的单机无穷大系统，可得VSG的系统运动方程如公式(1)所示。
其中，δ表示VSG的虚拟功角，t表示时间，ω表示VSG的虚拟转速，ωb表示系统转速，J表示VSG的虚拟惯量，Tm表示VSG的参考转矩，Te表示VSG的电磁转矩，Dp表示VSG的虚拟阻尼。
基于公式(1)，可得忽略虚拟阻尼影响的VSG系统运动方程如公式(2)所示。
令上式为：
上式有系统具有首次积分，根据公式(2)可得公式(4)。
其中，E表示初始暂态能量函数，J表示VSG的虚拟惯量，ω表示VSG的虚拟转速，ωb表示系统转速，VVSG表示VSG的内电势，V表示无穷大节点的电压，δ表示VSG的虚拟功角，δs表示系统稳定平衡点处的功角，X∑为故障切除后系统的视在电抗，Tm表示VSG的参考转矩。
暂态能量函数对时间的导数如公式(5)。
对公式(5)在[0,t]进行积分并结合式(4)，可得公式(6)。
对VSG的系统运动方程的第二式，左右两边同乘以Δδ＝(δ-δs)，并在[0,t]进行积分可得公式(7)。
对公式(7)左边第一项进行变形，得到公式(8)。
将公式(8)带入到公式(7)，可得到公式(9)。
由此，得到接有VSG的单机无穷大系统的暂态能量函数模型如公式(10)。
其中，E(δ,ω)表示暂态能量函数，J表示VSG的虚拟惯量，ω表示VSG的虚拟转速，ωb表示系统转速，VVSG表示VSG的内电势，X∑表示故障切除后系统的视在电抗，V表示无穷大节点的电压，δ表示VSG的虚拟功角，δs表示系统稳定平衡点处的功角，Tm表示VSG的参考转矩，Dp表示VSG的虚拟阻尼，λ表示微调系数，0＜λ＜1，Δδ＝δ-δs，ωb＝2πfn，fn表示系统频率。
最后，根据时域仿真的结果和暂态能量函数进行对比，验证暂态能量函数的准确性，并求出临界切除时间。
本发明实施例提供的虚拟同步发电机并网的暂态稳定性预测方法，构建接有VSG的单机无穷大系统的暂态能量函数，并计及虚拟阻尼的影响，建立的暂态能量函数准确度很高，通过计算暂态能量函数对虚拟同步发电机并网的暂态稳定性进行预测，计算量小，使用方便，预测效率高。
基于上述任一实施例，进一步地，所述暂态能量函数模型中包括虚拟转子动能、磁性势能、耗散势能、位置势能和虚拟阻尼消耗的能量。
具体来说，暂态能量函数模型中包括虚拟转子动能、磁性势能、耗散势能、位置势能和虚拟阻尼消耗的能量。
暂态能量函数模型用公式表示如下：
其中，E(δ,ω)表示暂态能量函数，J表示VSG的虚拟惯量，ω表示VSG的虚拟转速，ωb表示系统转速，VVSG表示VSG的内电势，X∑表示故障切除后系统的视在电抗，V表示无穷大节点的电压，δ表示VSG的虚拟功角，δs表示系统稳定平衡点处的功角，Tm表示VSG的参考转矩，Dp表示VSG的虚拟阻尼，λ表示微调系数，0＜λ＜1，Δδ＝δ-δs。
上述式子中，第一项表示VSG虚拟转子动能，第二项表示磁性势能和耗散势能，第三项表示位置势能，第四项和第五项表示虚拟阻尼消耗的能量。
本发明实施例提供的虚拟同步发电机并网的暂态稳定性预测方法，构建接有VSG的单机无穷大系统的暂态能量函数，并计及虚拟阻尼的影响，建立的暂态能量函数准确度很高，通过计算暂态能量函数对虚拟同步发电机并网的暂态稳定性进行预测，计算量小，使用方便，预测效率高。
基于上述任一实施例，图3为本发明实施例提供的虚拟同步发电机并网的暂态稳定性预测装置示意图，如图3所示，本发明实施例提供一种虚拟同步发电机并网的暂态稳定性预测装置，包括获取模块301、计算模块302和预测模块303，其中：
获取模块301用于获取虚拟同步发电机VSG的虚拟转速和虚拟功角；计算模块302用于将所述虚拟转速和虚拟功角，输入至预设的暂态能量函数模型，输出暂态能量值以及临界切除时间；预测模块303用于根据所述暂态能量值以及所述临界切除时间预测VSG并网的暂态稳定性。
基于上述任一实施例，进一步地，还包括模型构建模块，具体用于：
根据忽略虚拟阻尼影响的VSG系统运动方程求偏导得到初始暂态能量函数；
先对初始暂态能量函数对时间求导，然后进行时间积分，计及虚拟阻尼的影响，得到所述暂态能量函数模型。
基于上述任一实施例，进一步地，所述暂态能量函数模型用公式表示如下：
其中，E(δ,ω)表示暂态能量函数，J表示VSG的虚拟惯量，ω表示VSG的虚拟转速，ωb表示系统转速，VVSG表示VSG的内电势，X∑表示故障切除后系统的视在电抗，V表示无穷大节点的电压，δ表示VSG的虚拟功角，δs表示系统稳定平衡点处的功角，Tm表示VSG的参考转矩，Dp表示VSG的虚拟阻尼，λ表示微调系数，0＜λ＜1，Δδ＝δ-δs。
基于上述任一实施例，进一步地，所述暂态能量函数模型中包括虚拟转子动能、磁性势能、耗散势能、位置势能和虚拟阻尼消耗的能量。
本发明实施例提供一种虚拟同步发电机并网的暂态稳定性预测装置，用于执行上述任一实施例中所述的方法，通过本实施例提供的装置执行上述某一实施例中所述的方法的具体步骤与上述相应实施例相同，此处不再赘述。
本发明实施例提供的虚拟同步发电机并网的暂态稳定性预测装置，构建接有VSG的单机无穷大系统的暂态能量函数，并计及虚拟阻尼的影响，建立的暂态能量函数准确度很高，通过计算暂态能量函数对虚拟同步发电机并网的暂态稳定性进行预测，计算量小，使用方便，预测效率高。
图4为本发明实施例提供的电子设备的结构示意图，如图4所示，所述设备包括：处理器(processor)401、存储器(memory)402、总线403，以及存储在存储器上并可在处理器上运行的计算机程序。
其中，处理器401和存储器402通过总线403完成相互间的通信；
处理器401用于调用并执行存储器402中的计算机程序，以执行上述各方法实施例中的步骤，例如包括：
获取虚拟同步发电机VSG的虚拟转速和虚拟功角；
将所述虚拟转速和虚拟功角，输入至预设的暂态能量函数模型，输出暂态能量值以及临界切除时间；
根据所述暂态能量值以及所述临界切除时间预测VSG并网的暂态稳定性。
此外，上述的存储器中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
本发明实施例提供一种计算机程序产品，所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所述计算机程序包括程序指令，当所述程序指令被计算机执行时，计算机能够执行上述各方法实施例中的步骤，例如包括：
获取虚拟同步发电机VSG的虚拟转速和虚拟功角；
将所述虚拟转速和虚拟功角，输入至预设的暂态能量函数模型，输出暂态能量值以及临界切除时间；
根据所述暂态能量值以及所述临界切除时间预测VSG并网的暂态稳定性。
本发明实施例提供一种非暂态计算机可读存储介质，其上存储有计算机程序，当所述计算机程序被处理器执行时，实现上述各方法实施例中的步骤，例如包括：
获取虚拟同步发电机VSG的虚拟转速和虚拟功角；
将所述虚拟转速和虚拟功角，输入至预设的暂态能量函数模型，输出暂态能量值以及临界切除时间；
根据所述暂态能量值以及所述临界切除时间预测VSG并网的暂态稳定性。
以上所描述的装置及设备等实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种虚拟同步发电机并网的暂态稳定性预测方法，其特征在于，包括：
获取虚拟同步发电机VSG的虚拟转速和虚拟功角；
将所述虚拟转速和虚拟功角，输入至预设的暂态能量函数模型，输出暂态能量值以及临界切除时间；
根据所述暂态能量值以及所述临界切除时间预测VSG并网的暂态稳定性。
2.根据权利要求1所述的虚拟同步发电机并网的暂态稳定性预测方法，其特征在于，获取所述暂态能量函数模型的具体步骤如下：
根据忽略虚拟阻尼影响的VSG系统运动方程求偏导得到初始暂态能量函数；
先对初始暂态能量函数对时间求导，然后进行时间积分，计及虚拟阻尼的影响，得到所述暂态能量函数模型。
3.根据权利要求1所述的虚拟同步发电机并网的暂态稳定性预测方法，其特征在于，所述暂态能量函数模型用公式表示如下：
其中，E(δ,ω)表示暂态能量函数，J表示VSG的虚拟惯量，ω表示VSG的虚拟转速，ωb表示系统转速，VVSG表示VSG的内电势，X∑表示故障切除后系统的视在电抗，V表示无穷大节点的电压，δ表示VSG的虚拟功角，δs表示系统稳定平衡点处的功角，Tm表示VSG的参考转矩，Dp表示VSG的虚拟阻尼，λ表示微调系数，0＜λ＜1，Δδ＝δ-δs。
4.根据权利要求1所述的虚拟同步发电机并网的暂态稳定性预测方法，其特征在于，所述暂态能量函数模型中包括虚拟转子动能、磁性势能、耗散势能、位置势能和虚拟阻尼消耗的能量。
5.一种虚拟同步发电机并网的暂态稳定性预测装置，其特征在于，包括：
获取模块，用于获取虚拟同步发电机VSG的虚拟转速和虚拟功角；
计算模块，用于将所述虚拟转速和虚拟功角，输入至预设的暂态能量函数模型，输出暂态能量值以及临界切除时间；
预测模块，用于根据所述暂态能量值以及所述临界切除时间预测VSG并网的暂态稳定性。
6.根据权利要求5所述的虚拟同步发电机并网的暂态稳定性预测装置，其特征在于，还包括模型构建模块，具体用于：
根据忽略虚拟阻尼影响的VSG系统运动方程求偏导得到初始暂态能量函数；
先对初始暂态能量函数对时间求导，然后进行时间积分，计及虚拟阻尼的影响，得到所述暂态能量函数模型。
7.根据权利要求5所述的虚拟同步发电机并网的暂态稳定性预测装置，其特征在于，所述暂态能量函数模型用公式表示如下：
其中，E(δ,ω)表示暂态能量函数，J表示VSG的虚拟惯量，ω表示VSG的虚拟转速，ωb表示系统转速，VVSG表示VSG的内电势，X∑表示故障切除后系统的视在电抗，V表示无穷大节点的电压，δ表示VSG的虚拟功角，δs表示系统稳定平衡点处的功角，Tm表示VSG的参考转矩，Dp表示VSG的虚拟阻尼，λ表示微调系数，0＜λ＜1，Δδ＝δ-δs。
8.根据权利要求5所述的虚拟同步发电机并网的暂态稳定性预测装置，其特征在于，所述暂态能量函数模型中包括虚拟转子动能、磁性势能、耗散势能、位置势能和虚拟阻尼消耗的能量。
9.一种电子设备，包括存储器、处理器，以及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时，实现如权利要求1至4任一项所述虚拟同步发电机并网的暂态稳定性预测方法的步骤。
10.一种非暂态计算机可读存储介质，其上存储有计算机程序，其特征在于，当所述计算机程序被处理器执行时，实现如权利要求1至4任一所述虚拟同步发电机并网的暂态稳定性预测方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910662468.XA CN110504706B (zh) | 2019-07-22 | 2019-07-22 | 虚拟同步发电机并网的暂态稳定性预测方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910662468.XA CN110504706B (zh) | 2019-07-22 | 2019-07-22 | 虚拟同步发电机并网的暂态稳定性预测方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110504706A true CN110504706A (zh) | 2019-11-26 |
| CN110504706B CN110504706B (zh) | 2021-05-04 |
Family
ID=68586243
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910662468.XA Expired - Fee Related CN110504706B (zh) | 2019-07-22 | 2019-07-22 | 虚拟同步发电机并网的暂态稳定性预测方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110504706B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112117758A (zh) * | 2020-07-29 | 2020-12-22 | 西安交通大学 | 虚拟同步发电机多周期暂态稳定性分析方法及系统 |
| CN113452075A (zh) * | 2021-06-25 | 2021-09-28 | 国网冀北电力有限公司电力科学研究院 | 虚拟同步机虚拟功角确定方法及装置 |
| CN113552430A (zh) * | 2021-06-22 | 2021-10-26 | 南方电网科学研究院有限责任公司 | 一种基于临界指标的变换器暂态稳定性判定方法及装置 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3263813B2 (ja) * | 1996-03-26 | 2002-03-11 | ダイムラークライスラー・アクチエンゲゼルシヤフト | 回路網を構成する端末装置の接地の検査方法及び装置 |
| CN106294993A (zh) * | 2016-08-09 | 2017-01-04 | 南京理工大学 | 一种考虑逆变器电流饱和的暂态能量函数分析方法 |
| CN107330193A (zh) * | 2017-07-01 | 2017-11-07 | 南京理工大学 | 计及vsg逆变器电流限幅的暂态能量函数方法 |
| CN108493984A (zh) * | 2018-02-06 | 2018-09-04 | 国网四川省电力公司电力科学研究院 | 适用于光伏并网系统的虚拟同步发电机控制方法 |
| CN109659930A (zh) * | 2018-12-15 | 2019-04-19 | 南京理工大学 | 基于能量函数的含vsg-iidg电力系统暂态稳定分析方法 |
-
2019
- 2019-07-22 CN CN201910662468.XA patent/CN110504706B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3263813B2 (ja) * | 1996-03-26 | 2002-03-11 | ダイムラークライスラー・アクチエンゲゼルシヤフト | 回路網を構成する端末装置の接地の検査方法及び装置 |
| CN106294993A (zh) * | 2016-08-09 | 2017-01-04 | 南京理工大学 | 一种考虑逆变器电流饱和的暂态能量函数分析方法 |
| CN107330193A (zh) * | 2017-07-01 | 2017-11-07 | 南京理工大学 | 计及vsg逆变器电流限幅的暂态能量函数方法 |
| CN108493984A (zh) * | 2018-02-06 | 2018-09-04 | 国网四川省电力公司电力科学研究院 | 适用于光伏并网系统的虚拟同步发电机控制方法 |
| CN109659930A (zh) * | 2018-12-15 | 2019-04-19 | 南京理工大学 | 基于能量函数的含vsg-iidg电力系统暂态稳定分析方法 |
Non-Patent Citations (2)
| Title |
|---|
| JALIPOOR: "PowerSystemStabilizationUsingVirtualSyn^feronous GeneratorWithAlternatingMomentofInertia", 《EMERGING底SELECTEDTOPICSINPOWER巧ECTRONICS圧EEJOURNAL》 * |
| 孙震宇: "VSG策略下计及逆变器饱和的风电系统的暂态稳定分析方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112117758A (zh) * | 2020-07-29 | 2020-12-22 | 西安交通大学 | 虚拟同步发电机多周期暂态稳定性分析方法及系统 |
| CN113552430A (zh) * | 2021-06-22 | 2021-10-26 | 南方电网科学研究院有限责任公司 | 一种基于临界指标的变换器暂态稳定性判定方法及装置 |
| CN113452075A (zh) * | 2021-06-25 | 2021-09-28 | 国网冀北电力有限公司电力科学研究院 | 虚拟同步机虚拟功角确定方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110504706B (zh) | 2021-05-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Sanchez-Gasca et al. | A modal analysis of a two-area system with significant wind power penetration | |
| CN110504706B (zh) | 虚拟同步发电机并网的暂态稳定性预测方法及装置 | |
| Huang et al. | Open‐source framework for power system transmission and distribution dynamics co‐simulation | |
| Cui et al. | A particle filter for dynamic state estimation in multi-machine systems with detailed models | |
| Mahseredjian et al. | Simulation tools for electromagnetic transients in power systems: Overview and challenges | |
| Chavan et al. | Identification and predictive analysis of a multi-area WECC power system model using synchrophasors | |
| JP6725306B2 (ja) | 発電機の制動トルクの非侵入的な推定のためのシステムおよび方法 | |
| Papadopoulos et al. | Black‐box dynamic equivalent model for microgrids using measurement data | |
| CN106099922B (zh) | 基于启发式能量函数法的电力系统暂态电压稳定性判断方法 | |
| CN104201674A (zh) | 一种考虑负荷低电压释放特性的综合负荷模型的建模方法 | |
| CN102723712A (zh) | 一种提高电力系统暂态稳定性分析效率的方法 | |
| Makolo et al. | Online inertia estimation for power systems with high penetration of RES using recursive parameters estimation | |
| CN104810826A (zh) | 结合拉丁超立方抽样的双向迭代并行概率潮流计算方法 | |
| Wang et al. | Dynamic power flow algorithm considering frequency regulation of wind power generators | |
| Kyesswa et al. | A Matlab-based dynamic simulation module for power system transients analysis in the eASiMOV framework | |
| CN103324801B (zh) | 基于rtds的风力机数字仿真方法 | |
| Shetye et al. | Validation of power system transient stability results | |
| Larsen et al. | Generator dynamics in aeroelastic analysis and simulations | |
| Jiang et al. | Damping torques of multi-machine power systems during transient behaviors | |
| Fan et al. | A network-based structure-preserving dynamical model for the study of cascading failures in power grids | |
| Salam et al. | Transient Stability Analysis of a Three-machine Nine Bus Power System Network. | |
| El-Guindy et al. | Formal LPV control for transient stability of power systems | |
| CN106026182B (zh) | 考虑参数随机性的双馈风电机组仿真系统及其搭建方法 | |
| CN110867850B (zh) | 一种发电机转速及风电机组参数计算方法、风电机组模型 | |
| Eminoglu | A new model for wind turbine systems |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210504 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |