CN115688344B - 多变流器并网系统的宽频带振荡分析方法及系统 - Google Patents
多变流器并网系统的宽频带振荡分析方法及系统 Download PDFInfo
- Publication number
- CN115688344B CN115688344B CN202211703988.9A CN202211703988A CN115688344B CN 115688344 B CN115688344 B CN 115688344B CN 202211703988 A CN202211703988 A CN 202211703988A CN 115688344 B CN115688344 B CN 115688344B
- Authority
- CN
- China
- Prior art keywords
- converter
- grid
- connected system
- impedance
- subsystem
- 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
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开一种多变流器并网系统的宽频带振荡分析方法及系统,方法包括:对多变流器并网系统中的多变流器系统和线路网络系统进行划分,得到多变流器子系统和线路网络子系统;构建多变流器子系统在统一dq坐标系架构下的导纳模型,得到多变流器子系统的第一阻抗模型;构建线路网络子系统在统一dq坐标系架构下的阻抗矩阵模型,得到线路网络子系统的第二阻抗模型;根据第一阻抗模型和第二阻抗模型构建多变流器并网系统的数学模型,得到多变流器并网系统的广义阻抗比;根据广义阻抗比对多变流器并网系统进行宽频带振荡分析。广义阻抗比矩阵的阶数可自由扩展,能够适用于大规模高比例新能源并网系统的宽频带振荡分析。
Description
技术领域
本发明属于多变流器并网系统分析技术领域,尤其涉及一种多变流器并网系统的宽频带振荡分析方法及系统。
背景技术
电力电子变换器作为新能源并网的关键设备之一,在电力系统的发电、输电、配电、用电全过程中发挥着越来越重要的作用。与集中式同步发电机主导的电力系统不同,在电力电子变换器主导的新能源电力系统中可能包含成百上千个分布式的并网变流器。此外,电力电子变流器的动态特性更依赖于宽时间尺度控制器而不仅是无源元件的物理定律,导致近年来越来越多的宽频带振荡/不稳定事故频繁报道,造成重大经济损失的同时也给电网的安全稳定运行提出了更大的挑战。
目前针对单个变流器并网系统的稳定性分析已较为成熟。然而,多变流器并网时各个变流器端口处的电压电流存在相位差,使得各个变流器的阻抗模型难以统一到同一个坐标系中,造成多变流器并网的系统级稳定性分析困难,目前针对多变流器并网系统宽频带振荡分析方法的研究较少。
发明内容
本发明提供一种多变流器并网系统的宽频带振荡分析方法及系统,用于解决各个变流器的阻抗模型难以统一到同一个坐标系中,造成多变流器并网的系统级稳定性分析困难的技术问题。
第一方面,本发明提供一种多变流器并网系统的宽频带振荡分析方法,包括:对多
变流器并网系统中的多变流器系统和线路网络系统进行划分,得到多变流器子系统和
线路网络子系统;构建所述多变流器子系统在统一dq坐标系架构下的导纳模型,
得到所述多变流器子系统的第一阻抗模型;构建所述线路网络子系统在统一dq坐标
系架构下的阻抗矩阵模型,得到所述线路网络子系统的第二阻抗模型;根据所述第一
阻抗模型和所述第二阻抗模型构建所述多变流器并网系统的数学模型,得到所述多变流器
并网系统的广义阻抗比;根据所述广义阻抗比对所述多变流器并网系统进行宽频带振荡分
析。
第二方面,本发明提供一种多变流器并网系统的宽频带振荡分析系统,包括:划分
模块,配置为对多变流器并网系统中的多变流器系统和线路网络系统进行划分,得到多变
流器子系统和线路网络子系统;第一构建模块,配置为构建所述多变流器子系统
在统一dq坐标系架构下的导纳模型,得到所述多变流器子系统的第一阻抗模型;第二
构建模块,配置为构建所述线路网络子系统在统一dq坐标系架构下的阻抗矩阵模型,
得到所述线路网络子系统的第二阻抗模型;第三构建模块,配置为根据所述第一阻抗模
型和所述第二阻抗模型构建所述多变流器并网系统的数学模型,得到所述多变流器并网系
统的广义阻抗比;分析模块,配置为根据所述广义阻抗比对所述多变流器并网系统进行宽
频带振荡分析。
第三方面,提供一种电子设备,其包括:至少一个处理器,以及与所述至少一个处理器通信连接的存储器,其中,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行本发明任一实施例的多变流器并网系统的宽频带振荡分析方法的步骤。
第四方面,本发明还提供一种计算机可读存储介质,其上存储有计算机程序,所述程序指令被处理器执行时,使所述处理器执行本发明任一实施例的多变流器并网系统的宽频带振荡分析方法的步骤。
本申请的多变流器并网系统的宽频带振荡分析方法,利用多变流器并网系统的广义阻抗比模型构建等价开环传递函数来实现宽频带振荡分析,通过将多输入多输出系统转化为单输入单输出系统,实现了系统模型阶数的降低,计算量大大减小;利用经典的奈奎斯判据绘制单个子系统特征值轨迹,其分析过程更加简洁,有助于系统裕量的估计和设计,并且广义阻抗比矩阵的阶数可自由扩展,能够适用于大规模高比例新能源并网系统的宽频带振荡分析。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明一实施例提供的一种多变流器并网系统的宽频带振荡分析方法的流程图;
图2为本发明一实施例提供的又一种多变流器并网系统的宽频带振荡分析方法的流程图;
图3为本发明一具体实施例的多变流器系统典型拓扑图;
图4为本发明一具体实施例的多变流器系统负反馈模型示意图;
图5本发明一具体实施例的3个变流器并网系统广义阻抗比的奈奎斯特图;
图6本发明一具体实施例的修改参数前后变流器3并网点的A相电流波形图;
图7为本发明一实施例提供的一种多变流器并网系统的宽频带振荡分析系统的结构框图;
图8是本发明一实施例提供的电子设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1,其示出了本申请一种多变流器并网系统的宽频带振荡分析方法的流程图。
如图1所示,多变流器并网系统的宽频带振荡分析方法具体包括以下步骤:
步骤S104,根据所述第一阻抗模型和所述第二阻抗模型构建所述多变流器并网系统的数学模型,得到所述多变流器并网系统的广义阻抗比;
步骤S105,根据所述广义阻抗比对所述多变流器并网系统进行宽频带振荡分析。
本实施例的方法,利用多变流器并网系统的广义阻抗比模型构建等价开环传递函数来实现宽频带振荡分析,通过将多输入多输出系统转化为单输入单输出系统,实现了系统模型阶数的降低,计算量大大减小;利用经典的奈奎斯判据绘制单个子系统特征值轨迹,其分析过程更加简洁,有助于系统裕量的估计和设计,并且广义阻抗比矩阵的阶数可自由扩展,能够适用于大规模高比例新能源并网系统的宽频带振荡分析。
请参阅图2,其示出了本申请的又一种多变流器并网系统的宽频带振荡分析方法的流程图。
如图2所示,步骤S1、将新能源多变流器系统中的多变流器系统和线路网络系统进行划分;
由于每个变流器并网点实际的电压电流相位会存在一些偏差,即每个变流器阻抗模型的dq坐标系并不会完全重合,导致不能直接将每个变流器的阻抗模型进行合并化简。因此需要通过坐标变换方式,实现每个变流器阻抗(或导纳)模型的dq坐标系重合。
利用坐标变换公式,可以将所有并网变流器的端口电压、电流折算到同一个dq坐标系下,即:
由式(6)和式(7)对比可得到节点之间的阻抗为:
以图3所示的3个变流器并网系统实施例为例,线路网络阻抗矩阵的数学表达式如下:
步骤S4、基于广义阻抗比的多变流器系统宽频带振荡分析;
联立式(4)和式(5)可以推导出整个多变流器系统的数学模型,其表达式如下:
获取每个变流器的阻抗模型;
获取系统的广义阻抗比,通过绘制广义阻抗比的奈奎斯特曲线来分析系统的
宽频带振荡问题,当多变流并网系统无右半平面极点时,广义阻抗比的奈奎斯特图中顺
时针环绕(-1,0)的圈数等于逆时针环绕(-1,0)的圈数则系统稳定;当多变流并网系统有右
半平面极点时,广义阻抗比的奈奎斯特图中顺时针环绕(-1,0)的圈数和逆时针环绕(-1,
0)的圈数的差值等于右半平面极点数,则系统稳定。在不稳定的情况下,宽频带振荡的频率
可根据奈奎斯特曲线与单位圆的交点关系直接读取。
根据多变流器系统宽频带振荡分析流程,选择一种系统不稳定的参数情况作为案例,绘制整个系统广义阻抗比的奈奎斯特图如图5所示,由图5中的局部放大图可知:变流器1和变流器2广义阻抗比的奈奎斯特图重合,变流器1和变流器2广义阻抗比的奈奎斯特曲线没有包围(-1, j0)点,说明变流器1和变流器2能够稳定运行。然而,变流器3广义阻抗比的奈奎斯特曲线包围了(-1, j0)点,说明变流器3不能够稳定运行。
在此基础上,调节变流器3的控制参数,其中主要调节电流环的比例控制环节参数。图6给出了修改变流器3电流环的比例控制环节参数前后的仿真结果,图中为变流器3并网点的A相电流波形图,由图6可知,修改参数后,变流器3并网点的A相电流波形能够稳定。
通过仿真实例,可以看出:变流器3的稳定性情况能够很好地被所提一种新能源多变流器系统宽频带振荡分析方法进行判定,在修改参数前后多变流器系统的稳定性均能够进行准确分析,证明了本发明所提一种新能源多变流器系统宽频带振荡分析方法的正确性。
综上,该方法利用新能源多变流器系统的广义阻抗比模型构建等价开环传递函数来实现宽频带振荡分析。相比特征值分析方法、广义奈奎斯特分析方法、传统阻抗比分析方法,本发明所提方法具有如下优点:1)将多输入多输出系统转化为单输入单输出系统,实现了系统模型阶数的降低,计算量大大减小;2)利用经典的奈奎斯判据绘制单个子系统特征值轨迹,其分析过程更加简洁,有助于系统裕量的估计和设计。3)所提分析方法中系统广义阻抗比矩阵的阶数可自由扩展,能够适用于大规模高比例新能源并网系统的宽频带振荡分析。
请参阅图7,其示出了本申请的一种多变流器并网系统的宽频带振荡分析系统的结构框图。
如图7所示,宽频带振荡分析系统200,包括划分模块210、第一构建模块220、第二构建模块230、第三构建模块240以及分析模块250。
其中,划分模块210,配置为对多变流器并网系统中的多变流器系统和线路网络系
统进行划分,得到多变流器子系统和线路网络子系统;第一构建模块220,配置为构建
所述多变流器子系统在统一dq坐标系架构下的导纳模型,得到所述多变流器子系统的第一阻抗模型;第二构建模块230,配置为构建所述线路网络子系统在统一dq坐标
系架构下的阻抗矩阵模型,得到所述线路网络子系统的第二阻抗模型;第三构建模块
240,配置为根据所述第一阻抗模型和所述第二阻抗模型构建所述多变流器并网系统的数
学模型,得到所述多变流器并网系统的广义阻抗比;分析模块250,配置为根据所述广义阻
抗比对所述多变流器并网系统进行宽频带振荡分析。
应当理解,图7中记载的诸模块与参考图1中描述的方法中的各个步骤相对应。由此,上文针对方法描述的操作和特征以及相应的技术效果同样适用于图7中的诸模块,在此不再赘述。
在另一些实施例中,本发明实施例还提供了一种计算机可读存储介质,其上存储有计算机程序,所述程序指令被处理器执行时,使所述处理器执行上述任意方法实施例中的多变流器并网系统的宽频带振荡分析方法;
作为一种实施方式,本发明的计算机可读存储介质存储有计算机可执行指令,计算机可执行指令设置为:
根据所述第一阻抗模型和所述第二阻抗模型构建所述多变流器并网系统的数学模型,得到所述多变流器并网系统的广义阻抗比;
根据所述广义阻抗比对所述多变流器并网系统进行宽频带振荡分析。
计算机可读存储介质可以包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需要的应用程序;存储数据区可存储根据多变流器并网系统的宽频带振荡分析系统的使用所创建的数据等。此外,计算机可读存储介质可以包括高速随机存取存储器,还可以包括存储器,例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实施例中,计算机可读存储介质可选包括相对于处理器远程设置的存储器,这些远程存储器可以通过网络连接至多变流器并网系统的宽频带振荡分析系统。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
图8是本发明实施例提供的电子设备的结构示意图,如图8所示,该设备包括:一个处理器310以及存储器320。电子设备还可以包括:输入装置330和输出装置340。处理器310、存储器320、输入装置330和输出装置340可以通过总线或者其他方式连接,图8中以通过总线连接为例。存储器320为上述的计算机可读存储介质。处理器310通过运行存储在存储器320中的非易失性软件程序、指令以及模块,从而执行服务器的各种功能应用以及数据处理,即实现上述方法实施例多变流器并网系统的宽频带振荡分析方法。输入装置330可接收输入的数字或字符信息,以及产生与多变流器并网系统的宽频带振荡分析系统的用户设置以及功能控制有关的键信号输入。输出装置340可包括显示屏等显示设备。
上述电子设备可执行本发明实施例所提供的方法,具备执行方法相应的功能模块和有益效果。未在本实施例中详尽描述的技术细节,可参见本发明实施例所提供的方法。
作为一种实施方式,上述电子设备应用于多变流器并网系统的宽频带振荡分析系统中,用于客户端,包括:至少一个处理器;以及,与至少一个处理器通信连接的存储器;其中,存储器存储有可被至少一个处理器执行的指令,指令被至少一个处理器执行,以使至少一个处理器能够:
根据所述第一阻抗模型和所述第二阻抗模型构建所述多变流器并网系统的数学模型,得到所述多变流器并网系统的广义阻抗比;
根据所述广义阻抗比对所述多变流器并网系统进行宽频带振荡分析。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分的方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (6)
1.一种多变流器并网系统的宽频带振荡分析方法,其特征在于,包括:
式中,,代表多变流器整体的输入电流矩阵,,代表多变流器整体的输入电压矩阵,/>,代表多变流器整体的小扰动电流信号矩阵,/>,代表多变流器子系统/>的第一阻抗模型,/>,代表多变流器整体的坐标变换矩阵;
根据所述第一阻抗模型和所述第二阻抗模型构建所述多变流器并网系统的数学模型,得到所述多变流器并网系统的广义阻抗比,其中,所述多变流器并网系统的数学模型的表达式为:
根据所述广义阻抗比对所述多变流器并网系统进行宽频带振荡分析。
2.根据权利要求1所述的一种多变流器并网系统的宽频带振荡分析方法,其特征在于,所述根据所述广义阻抗比对所述多变流器并网系统进行宽频带振荡分析,包括:
通过绘制所述广义阻抗比的奈奎斯特曲线图对所述多变流器并网系统进行宽频带振荡分析,具体为:
当多变流器并网系统无右半平面极点时,广义阻抗比的奈奎斯特曲线图中顺时针环绕(-1,0)的圈数等于逆时针环绕(-1,0)的圈数,则多变流器并网系统稳定;
当多变流器并网系统有右半平面极点时,广义阻抗比的奈奎斯特曲线图中顺时针环绕(-1,0)的圈数和逆时针环绕(-1,0)的圈数的差值等于右半平面极点数,则多变流器并网系统稳定;
在多变流器并网系统不稳定的情况下,宽频带振荡的频率可根据奈奎斯特曲线与单位圆的交点关系直接读取。
4.一种多变流器并网系统的宽频带振荡分析系统,其特征在于,包括:
式中,,代表多变流器整体的输入电流矩阵,,代表多变流器整体的输入电压矩阵,/>,代表多变流器整体的小扰动电流信号矩阵,/>,代表多变流器子系统/>的第一阻抗模型,/>,代表多变流器整体的坐标变换矩阵;
第三构建模块,配置为根据所述第一阻抗模型和所述第二阻抗模型构建所述多变流器并网系统的数学模型,得到所述多变流器并网系统的广义阻抗比,其中,所述多变流器并网系统的数学模型的表达式为:
分析模块,配置为根据所述广义阻抗比对所述多变流器并网系统进行宽频带振荡分析。
5.一种电子设备,其特征在于,包括:至少一个处理器,以及与所述至少一个处理器通信连接的存储器,其中,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行权利要求1至3任一项所述的方法。
6.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述程序被处理器执行时实现权利要求1至3任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211703988.9A CN115688344B (zh) | 2022-12-29 | 2022-12-29 | 多变流器并网系统的宽频带振荡分析方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211703988.9A CN115688344B (zh) | 2022-12-29 | 2022-12-29 | 多变流器并网系统的宽频带振荡分析方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115688344A CN115688344A (zh) | 2023-02-03 |
CN115688344B true CN115688344B (zh) | 2023-03-28 |
Family
ID=85055963
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211703988.9A Active CN115688344B (zh) | 2022-12-29 | 2022-12-29 | 多变流器并网系统的宽频带振荡分析方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115688344B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106786776A (zh) * | 2017-02-15 | 2017-05-31 | 云南电网有限责任公司 | 一种利用修正广义阻抗法分析并网逆变器系统稳定性的方法 |
CN109449958A (zh) * | 2018-11-29 | 2019-03-08 | 云南电网有限责任公司 | 一种双馈风机并网系统稳定性分析方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106532685B (zh) * | 2016-10-26 | 2019-02-01 | 浙江大学 | 用于并网逆变器稳定分析的广义阻抗判据计算方法及应用 |
CN107069811B (zh) * | 2017-04-12 | 2019-07-26 | 清华大学 | 基于同步参考坐标系的阻抗网络建模与稳定性分析方法 |
US20210330279A1 (en) * | 2020-04-27 | 2021-10-28 | Ess Technology, Inc. | Audio Output as Half Duplex Biometric Detection Device |
-
2022
- 2022-12-29 CN CN202211703988.9A patent/CN115688344B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106786776A (zh) * | 2017-02-15 | 2017-05-31 | 云南电网有限责任公司 | 一种利用修正广义阻抗法分析并网逆变器系统稳定性的方法 |
CN109449958A (zh) * | 2018-11-29 | 2019-03-08 | 云南电网有限责任公司 | 一种双馈风机并网系统稳定性分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115688344A (zh) | 2023-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10826292B2 (en) | Multi-time-scale digital/analog hybrid simulation system and method for power distribution network and storage medium | |
CN109586269B (zh) | 考虑参数自寻优的直流微电网虚拟惯性控制方法及系统 | |
Cui et al. | Reinforcement learning for optimal primary frequency control: A Lyapunov approach | |
CN110556831B (zh) | 一种多机多节点发电系统的稳定性评估方法及系统 | |
CN103810646B (zh) | 一种基于改进投影积分算法的有源配电系统动态仿真方法 | |
Islam et al. | RTDS implementation of an improved sliding mode based inverter controller for PV system | |
US11616365B2 (en) | Decentralized hardware-in-the-loop scheme | |
Liu et al. | Modeling and simulation of hybrid AC-DC system on a cloud computing based simulation platform-CloudPSS | |
CN108873733A (zh) | 电力信息物理系统中信息预想事故影响的分析方法 | |
CN111181164B (zh) | 一种改进的主从分裂输配协同潮流计算方法及系统 | |
CN115700956A (zh) | 基于多端口矩阵的光伏集群直流外送振荡分析方法及系统 | |
Yuan et al. | Dynamic equivalent model of VSC based on singular perturbation | |
Shrivastava et al. | Distributed, fixed‐time, and bounded control for secondary voltage and frequency restoration in islanded microgrids | |
Guzmán et al. | Microgrid power sharing framework for software defined networking and cybersecurity analysis | |
Cheng et al. | ECS-grid: Data-oriented real-time simulation platform for cyber-physical power systems | |
Kolluri et al. | Controlling DC microgrids in communities, buildings and data centers | |
CN115688344B (zh) | 多变流器并网系统的宽频带振荡分析方法及系统 | |
CN117353288A (zh) | 一种多流融合的有源配电网仿真分析方法、系统及介质 | |
Long et al. | Data-driven hybrid equivalent dynamic modeling of multiple photovoltaic power stations based on ensemble gated recurrent unit | |
CN111682552A (zh) | 电压控制方法、装置、设备及存储介质 | |
Guddanti et al. | Better data structures for co-simulation of distribution system with gridlab-d and python | |
Han et al. | Analysis and verification of harmonic interaction between DDWTs based on αβ stationary coordinate system | |
Zhao et al. | Modelling and design of parameters for fast simulation of electromagnetic transient in AC/DC hybrid microgrid | |
Yuan et al. | Assessment-prediction-regulation coordinated secondary control strategy for microgrid considering large time-delay | |
Liang et al. | A parallel algorithm of optimal power flow on Hadoop platform |
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 |