CN111310325A - 一种模块化多电平换流器的动态仿真模拟方法及系统 - Google Patents

一种模块化多电平换流器的动态仿真模拟方法及系统 Download PDF

Info

Publication number
CN111310325A
CN111310325A CN202010084671.6A CN202010084671A CN111310325A CN 111310325 A CN111310325 A CN 111310325A CN 202010084671 A CN202010084671 A CN 202010084671A CN 111310325 A CN111310325 A CN 111310325A
Authority
CN
China
Prior art keywords
multilevel converter
modular multilevel
level
representing
model
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
CN202010084671.6A
Other languages
English (en)
Other versions
CN111310325B (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.)
North China Electric Power University
Original Assignee
North China Electric Power 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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN202010084671.6A priority Critical patent/CN111310325B/zh
Publication of CN111310325A publication Critical patent/CN111310325A/zh
Application granted granted Critical
Publication of CN111310325B publication Critical patent/CN111310325B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Inverter Devices (AREA)

Abstract

本发明提出了一种模块化多电平换流器的动态仿真模拟方法。本发明的方法首先建立模块化多电平换流器的开关函数模型;对所述开关函数模型进行分解、重组和移频获得模块化多电平换流器的多频段动态向量模型;然后,基于自适应策略原理,根据模块化多电平换流器的工况自适应切换多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真,所述工况为稳态工况或暂态工况。本发明通过频率的分解重组和移频,并采用自适应方法,进行模块化多电平换流器的动态仿真模拟,在保证精度的前提下,增大步长提高仿真速度。

Description

一种模块化多电平换流器的动态仿真模拟方法及系统
技术领域
本发明涉及换流器性能研究技术领域,特别是涉及一种模块化多电平换流器的动态仿真模拟方法及系统。
背景技术
对于电力系统的安全稳定运行来说,电磁暂态是最重要的系统分析手段,但是随着电力电子设备在电力系统中的应用,电磁暂态仿真面临着仿真速度和仿真精度之间的矛盾,为此人们不断改进电力系统各元件的电磁暂态模型和仿真算法。模块化多电平换流器(MMC)凭借失真度小、开关损耗低、控制灵活等优点,已成为现今最具有发展前景的电力电子技术和设备。电磁暂态仿真中, MMC传统动态相量模型,不用通过小步长来精确捕捉高频开关的变化时刻,但是随着考虑的谐波次数的增加,MMC动态相量模型的微分方程组数将急剧增加,仿真规模将急剧增大,往往只能采取谐波截断的方式,减少计算量提高仿真速度,使得仿真精度较低。MMC频率偏移法模型。利用Hilbert变换构造解析信号,以工频为主导频率移频,信号频率降低后采用大步长仿真,但是该模型将信号带宽限制在工频基波附近,仿真精度低。为了提高仿真精度MMC 一般采用详细模型,只能通过小步长仿真来捕捉高频开关的精确变化,且每有开关动作时须更新系统导纳矩阵。随着MMC子模块数目的增长,仿真规模的扩大,MMC详细模型将会带来极大的仿真负担,仿真速度迅速降低。如何在提高MMC动态仿真精度的同时提高仿真的速度成为一个亟待解决的技术问题。
发明内容
本发明的目的是提供一种模块化多电平换流器的动态仿真模拟方法及系统,以实现在保证MMC动态仿真精度的同时提高仿真的速度。
为实现上述目的,本发明提供了如下方案:
一种模块化多电平换流器的动态仿真模拟方法,所述模拟方法包括如下步骤:
根据模块化多电平换流器的拓扑结构,建立模块化多电平换流器的开关函数模型;
基于多频段动态向量法原理,对所述开关函数模型进行分解、重组和移频获得模块化多电平换流器的多频段动态向量模型;
基于自适应策略原理,根据模块化多电平换流器的工况自适应切换多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真,所述工况为稳态工况或暂态工况。
可选的,所述根据模块化多电平换流器的拓扑结构,建立模块化多电平换流器的开关函数模型,具体包括:
根据模块化多电平换流器的拓扑结构,建立模块化多电平换流器的开关函数模型为:
Figure RE-GDA0002455452580000021
Figure RE-GDA0002455452580000022
Figure RE-GDA0002455452580000023
Figure RE-GDA0002455452580000024
其中,Lx表示模块化多电平换流器的第x相的桥臂电抗,Rx模块化多电平换流器的第x相的桥臂电阻,x为A、B或C,
Figure RE-GDA0002455452580000025
表示模块化多电平换流器的第x相的上桥臂与下桥臂的电流之和,
Figure RE-GDA0002455452580000026
表示模块化多电平换流器的第x相的上桥臂与下桥臂的电流之差,
Figure RE-GDA0002455452580000027
表示模块化多电平换流器的第x相的上桥臂与下桥臂的模块电压之和和
Figure RE-GDA0002455452580000028
表示模块化多电平换流器的第x相的上桥臂与下桥臂的模块电压之差,N表示每个桥臂的子模块的数量,Csm表示模块化多电平换流器的每个模块的电容,
Figure RE-GDA0002455452580000029
Figure RE-GDA00024554525800000210
分别表示模块化多电平换流器的第x 相的上桥臂与下桥臂的开关函数之和和开关函数之差,ux为模块化多电平换流器的第x相的桥臂、udc为模块化多电平换流器的直流侧电压。
可选的,所述基于多频段动态向量法原理,对所述开关函数模型进行分解、重组和移频获得模块化多电平换流器的多频段动态向量模型,具体包括:
基于多频段动态向量法原理,对所述开关函数模型进行分解、重组和移频获得模块化多电平换流器的多频段动态向量模型为:
Figure RE-GDA00024554525800000211
Figure RE-GDA0002455452580000031
Figure RE-GDA0002455452580000032
Figure RE-GDA0002455452580000033
其中,
Figure RE-GDA0002455452580000034
表示子频段动态向量,fh为第h阶谐波的频率,Bk表示第k个频段信号重组后的子频段信号,Rx表示模块化多电平换流器的第x相的桥臂电阻,h表示第h阶谐波,ωk表示第k个频段的所移频率,j表示虚数的单位。
可选的,所述基于自适应策略原理,根据模块化多电平换流器的工况自适应切换多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真,具体包括:
当模块化多电平换流器的工况为稳态工况时,采用工频周期作为多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真;
当模块化多电平换流器的工况为暂态工况时,采用工频周期内每种电平持续时间的平均时间作为多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真。
可选的,所述当模块化多电平换流器的工况为暂态工况时,采用工频周期内每种电平持续时间的平均时间平作为多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真,具体包括:
判断滑动窗向前的过程中是否存在电平阶跃,得到判断结果;
若所述判断结果表示所述滑动窗向前的过程中存在电平阶跃,则对模块化多电平换流器的有电平阶跃的开关函数:
Figure RE-GDA0002455452580000035
进行傅里叶分解,得到有电平阶跃的开关函数在频域内的表达式:
Figure RE-GDA0002455452580000041
其中,ts表示平均时间,n 表示平均时间内电平的数量,t1、t2、tn和tn+1分别表示第一种电平、第二种电平、第n种电平和第n+1种电平的结束时间,t表示当前时间,i表示第一种电平的闭合开关的数量,A0表示对有电平阶跃的开关函数进行傅里叶分解时得到的直流量、ck对有电平阶跃的开关函数进行傅里叶分解得到的第k阶信号分量的谐波系数,ω2表示平均时间尺度对应的频率、k表示第k阶信号分量、δx表示x相的动态相位和
Figure RE-GDA0002455452580000042
表示x相的固定相位;
利用有电平阶跃的开关函数在频域内的表达式,对所述多频段动态向量模型进行仿真;
若所述判断结果表示所述滑动窗向前的过程中不存在电平阶跃,对模块化多电平换流器的无电平阶跃的开关函数:s(ts)=lusm-x(l<N+1),进行傅里叶分解,得到无电平阶跃的开关函数在频域内的表达式:
Figure RE-GDA0002455452580000043
其中, l为滑动窗口内电平的数量;
利用无电平阶跃的开关函数在频域内的表达式,对所述多频段动态向量模型进行仿真。
一种模块化多电平换流器的动态仿真模拟系统,所述模拟系统包括:
开关函数模型建立模块,用于根据模块化多电平换流器的拓扑结构,建立模块化多电平换流器的开关函数模型;
多频段动态向量模型获取模块,用于基于多频段动态向量法原理,对所述开关函数模型进行分解、重组和移频获得模块化多电平换流器的多频段动态向量模型;
自适应仿真模块,用于基于自适应策略原理,根据模块化多电平换流器的工况自适应切换多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真,所述工况为稳态工况或暂态工况。
可选的,所述开关函数模型建立模块,具体包括:
开关函数模型建立子模块,用于根据模块化多电平换流器的拓扑结构,建立模块化多电平换流器的开关函数模型为:
Figure RE-GDA0002455452580000051
Figure RE-GDA0002455452580000052
Figure RE-GDA0002455452580000053
Figure RE-GDA0002455452580000054
其中,Lx表示模块化多电平换流器的第x相的桥臂电抗,Rx模块化多电平换流器的第x相的桥臂电阻,x为A、B或C,
Figure RE-GDA0002455452580000055
表示模块化多电平换流器的第x相的上桥臂与下桥臂的电流之和,
Figure RE-GDA0002455452580000056
表示模块化多电平换流器的第x相的上桥臂与下桥臂的电流之差,
Figure RE-GDA0002455452580000057
表示模块化多电平换流器的第x相的上桥臂与下桥臂的模块电压之和和
Figure RE-GDA0002455452580000058
表示模块化多电平换流器的第x相的上桥臂与下桥臂的模块电压之差,N表示每个桥臂的子模块的数量,Csm表示模块化多电平换流器的每个模块的电容,
Figure RE-GDA0002455452580000059
Figure RE-GDA00024554525800000510
分别表示模块化多电平换流器的第x 相的上桥臂与下桥臂的开关函数之和和开关函数之差,ux为模块化多电平换流器的第x相的桥臂、udc为模块化多电平换流器的直流侧电压。
可选的,所述多频段动态向量模型获取模块,具体包括:
多频段动态向量模型获取子模块,用于基于多频段动态向量法原理,对所述开关函数模型进行分解、重组和移频获得模块化多电平换流器的多频段动态向量模型为:
Figure RE-GDA00024554525800000511
Figure RE-GDA00024554525800000512
Figure RE-GDA00024554525800000513
Figure RE-GDA00024554525800000514
其中,
Figure RE-GDA00024554525800000515
表示子频段动态向量,fh为第h阶谐波的频率,Bk表示第k个频段信号重组后的子频段信号,Rx表示模块化多电平换流器的第x相的桥臂电阻,h表示第h阶谐波,ωk表示第k个频段的所移频率,j表示虚数单位。
可选的,所述自适应仿真模块,具体包括:
第一自适应仿真子模块,用于当模块化多电平换流器的工况为稳态工况时,采用工频周期作为多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真;
第二自适应仿真子模块,当模块化多电平换流器的工况为暂态工况时,采用工频周期内每种电平持续时间的平均时间作为多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真。
可选的,所述第二自适应仿真子模块,具体包括:
判断单元,用于判断滑动窗向前的过程中是否存在电平阶跃,得到判断结果;
第一傅里叶分解单元,用于若所述判断结果表示所述滑动窗向前的过程中存在电平阶跃,则对模块化多电平换流器的有电平阶跃的开关函数:
Figure RE-GDA0002455452580000061
进行傅里叶分解,得到有电平阶跃的开关函数在频域内的表达式:
Figure RE-GDA0002455452580000062
其中,ts表示平均时间,n表示平均时间内电平的数量,t1、t2、tn和tn+1分别表示第一种电平、第二种电平、第n种电平和第n+1种电平的结束时间,t表示当前时间,i表示第一种电平的闭合开关的数量,A0表示对有电平阶跃的开关函数进行傅里叶分解时得到的直流量、ck对有电平阶跃的开关函数进行傅里叶分解得到的第k阶信号分量的谐波系数,ω2表示平均时间尺度对应的频率、k表示第k阶信号分量、δx表示x相的动态相位和
Figure RE-GDA0002455452580000063
表示x相的固定相位;
第一仿真单元,用于利用有电平阶跃的开关函数在频域内的表达式,对所述多频段动态向量模型进行仿真;
第二傅里叶分解单元,用于若所述判断结果表示所述滑动窗向前的过程中不存在电平阶跃,对模块化多电平换流器的无电平阶跃的开关函数: s(ts)=lusm-x(l<N+1),进行傅里叶分解,得到无电平阶跃的开关函数在频域内的表达式:
Figure RE-GDA0002455452580000071
其中,l为滑动窗口内电平的数量;
第一仿真单元,用于利用无电平阶跃的开关函数在频域内的表达式,对所述多频段动态向量模型进行仿真。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提出了一种模块化多电平换流器的动态仿真模拟方法。本发明的方法首先建立模块化多电平换流器的开关函数模型;对所述开关函数模型进行分解、重组和移频获得模块化多电平换流器的多频段动态向量模型;然后,基于自适应策略原理,根据模块化多电平换流器的工况自适应切换多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真,所述工况为稳态工况或暂态工况。本发明通过频率的分解重组和移频,并采用自适应方法,进行模块化多电平换流器的动态仿真模拟,在保证精度的前提下,增大步长提高仿真速度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的一种模块化多电平换流器的动态仿真模拟方法的流程图;
图2为本发明提供的一种模块化多电平换流器的动态仿真模拟方法的流程图;
图3为本发明提供的模块化多电平换流器的拓扑结构图;
图4为本发明提供的最近电平逼近调制波形图;
图5为本发明提供的频带分段原理示意图;
图6为本发明提供的自适应策略原理示意图;
图7为本发明提供的单个平均开关周期与输出波形的关系图;
图8为本发明提供的平均开关周期内出现的电平阶跃情况示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种模块化多电平换流器的动态仿真模拟方法及系统,以实现在保证MMC动态仿真精度的同时提高仿真的速度。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
为了实现上述目的本发明提供一种模块化多电平换流器的动态仿真模拟方法,如图1和2所示,所述模拟方法包括如下步骤:
步骤101,根据模块化多电平换流器的拓扑结构,建立模块化多电平换流器的开关函数模型。
所述模块化多电平换流器包括三个相单元(a、b、c),每相分为上、下桥臂,每个桥臂均由N个子模块SM、1个电阻R和1个电抗器L串联而成。如图3 所示,图3中,用u、l分别表示上下桥臂,设MMC换流器的x相(x=a,b,c) 上、下桥臂等值电阻和电抗相同,分别为Rx和Lx,上、下桥臂各子模块的电容大小相同为Csm,上、下桥臂电流分别为iu x和ilx,换流器交流端口电压为 ux,注入电流为ix,换流器直流侧电流为idc,直流侧电压为udc,各电气量方向如图3所示。
建立模块化多电平换流器的开关函数模型的具体步骤为:
定义子模块的开关函数为:
Figure RE-GDA0002455452580000081
其中,sk-x=1表示第k个子模块耦合到桥臂,参与桥臂运行;sk-x=0表示第 k个子模块被旁路,不参与桥臂运行。
根据MMC换流器的特点,由子模块开关函数可得上、下桥臂开关函数:
Figure RE-GDA0002455452580000082
定义上、下桥臂开关函数之和与差分别为:
Figure RE-GDA0002455452580000091
Figure RE-GDA0002455452580000092
定义上、下桥臂电流之和为:
Figure RE-GDA0002455452580000093
定义上、下桥臂电流之差为:
Figure RE-GDA0002455452580000094
定义上、下桥臂子模块电容电压之和为:
Figure RE-GDA0002455452580000095
定义上、下桥臂子模块电容电压之差为:
Figure RE-GDA0002455452580000096
由图3和上述定义,根据基尔霍夫电流和电压定律,可得MMC的开关函数模型为:
Figure RE-GDA0002455452580000097
Figure RE-GDA0002455452580000098
Figure RE-GDA0002455452580000099
Figure RE-GDA00024554525800000910
一般的,MMC的输出波形按照由最近电平逼近调制策略进行调整,调整波形如图4所示,图4中。平滑的曲线表示调制波,锯齿形曲线表示NLM波;。
定义调制比:
Figure RE-GDA00024554525800000911
其中:Um为调制波幅值。
则第i个电平的起始角度θi,按最近电平逼近原则可求得为:
Figure RE-GDA00024554525800000912
对图4进行傅立叶分解可得:
Figure RE-GDA00024554525800000913
其中:
Figure RE-GDA0002455452580000101
开关函数模型式(10)中的
Figure RE-GDA0002455452580000102
Figure RE-GDA0002455452580000103
如式13所示。
可见,根据模块化多电平换流器的拓扑结构,建立模块化多电平换流器的开关函数模型为:
Figure RE-GDA0002455452580000104
Figure RE-GDA0002455452580000105
Figure RE-GDA0002455452580000106
Figure RE-GDA0002455452580000107
其中,Lx表示模块化多电平换流器的第x相的桥臂电抗,Rx模块化多电平换流器的第x相的桥臂电阻,x为A、B或C,
Figure RE-GDA0002455452580000108
表示模块化多电平换流器的第x相的上桥臂与下桥臂的电流之和,
Figure RE-GDA0002455452580000109
表示模块化多电平换流器的第x相的上桥臂与下桥臂的电流之差,
Figure RE-GDA00024554525800001010
表示模块化多电平换流器的第x相的上桥臂与下桥臂的模块电压之和和
Figure RE-GDA00024554525800001011
表示模块化多电平换流器的第x相的上桥臂与下桥臂的模块电压之差,N表示每个桥臂的子模块的数量,Csm表示模块化多电平换流器的每个模块的电容,
Figure RE-GDA00024554525800001012
Figure RE-GDA00024554525800001013
分别表示模块化多电平换流器的第x 相的上桥臂与下桥臂的开关函数之和和开关函数之差,ux为模块化多电平换流器的第x相的桥臂、udc为模块化多电平换流器的直流侧电压。
步骤102,基于多频段动态向量法原理,对所述开关函数模型进行分解、重组和移频获得模块化多电平换流器的多频段动态向量模型。
其中,根据多频段动态向量法原理,进行分解、重组和移频的具体步骤包括:
1、信号的频率分解和重组
电力系统中电压、电流等可看作周期变化的电气信号。对周期为T0的电气信号x(τ),在一个周期τ∈(t-T0,t]内,其复数形式的傅里叶分解为:
Figure RE-GDA0002455452580000111
式中:X(t)为电气信号x(τ)的复数形式,ωs=2π/T0为基波角频率,Xh(t)为第h 阶傅立叶系数,即h阶“动态相量”。
式(1)中的h理论上无穷大,但在电磁暂态仿真中,根据精度和采样定理要求,仿真步长对应的频率一般为信号频率的10倍。因此,电磁暂态仿真中,谐波次数h一般为有限值,其最大值可依据仿真步长来确定(如:50μs仿真步长对应的最大谐波次数h为40)。此时,式(1)可写为:
Figure RE-GDA0002455452580000112
傅立叶变换的数学意义为:任意满足狄里赫利条件的周期信号都可以表示为一组以
Figure RE-GDA0002455452580000113
为正交基的线性组合。如果按图5将信号频率分成N个子频段,由线性组合的结合律,式(2)可按子频段分组和重新组合,组合前后的信号相等,表示的是同一个信号X(t)。
每个子频段重新组合的结果是:各频段内由傅立叶分解得到的频率不同的多个子信号,重新组合为1个频段信号。此时,X(t)可看作各频段信号的和,即:
Figure RE-GDA0002455452580000114
其中,Bn(t)为第n个频段重组后的子频段信号。
比较式(2)和式(3)右边的项数,可发现式(2)的项数大为减少,由谐波个数 2M减少为频段数2N。上述过程称为信号的频率分解和重组。
2、子频段动态相量法
式(3)将复数信号X(t)按频率由小到大进行了分段重组。对于每个子频段信号Bn(t),可选择频段内的一个频率(如:各频段的中心频率)分别移频,即:
Figure RE-GDA0002455452580000115
其中,
Figure RE-GDA0002455452580000116
为子频段移频后的信号,简记为
Figure RE-GDA0002455452580000117
ωrn为频段n的移频角频率,
Figure RE-GDA0002455452580000118
为频段n的频率上下限,fh为该频段内信号的频率,且满足
Figure RE-GDA0002455452580000121
Figure RE-GDA0002455452580000122
为信号的子频段动态相量(sub-frequency band dynamic phasor,SFB-DP)。显然,与传统的动态相量Xh(t)只有一个频率成份不同,子频段动态相量
Figure RE-GDA0002455452580000123
是一个具有带宽的信号。如果带宽满足窄带条件,则
Figure RE-GDA0002455452580000124
也是一个低频信号,而对低频信号可以采用大步长仿真,从而提高仿真速度。
3、多频段动态相量法
传统的动态相量法进行电磁暂态仿真时,先将信号按式(2)分解为各阶动态相量,然后利用CPU多核芯的特点,将各阶动态相量放在不同的CPU核中进行并行仿真。但是,由于CPU核芯数有限,传统动态相量法只能取少数阶动态相量构成方程组进行并行仿真。因此,传统动态相量法仿真的信号总带宽将远远小于信号的实际带宽,导致谐波截断误差大,仿真精度低。与传统的动态相量法不同,如果按照式(3)分解移频后的子频段动态相量进行并行仿真,由于子频段信号具有一定的带宽,因此,在与传统动态相量法具有相同方程组数的情况下,其可仿真的信号的带宽将远远大于传统的动态相量法,谐波截断误差将远远小于传统的动态相量法,因此具有极高的仿真精度。
上述过程称为多频段动态相量法(multi-frequency band dynamic phasor,MFB-DP),可用图6表示。
本发明将多频段动态相量原理应用于式(10)的电压、电流和开关函数,得到 MMC的多频段动态相量模型,即:
Figure RE-GDA0002455452580000125
Figure RE-GDA0002455452580000126
Figure RE-GDA0002455452580000127
Figure RE-GDA0002455452580000128
应用式(15)就可进行大步长电磁暂态仿真求出各频段的
Figure RE-GDA0002455452580000129
Figure RE-GDA00024554525800001210
将求出的
Figure RE-GDA0002455452580000131
反向移频,将得到各频段结果相加后取实部就得到了电磁暂态仿真所需要的实数信号。即:
Figure RE-GDA0002455452580000132
Figure RE-GDA0002455452580000133
Figure RE-GDA0002455452580000134
Figure RE-GDA0002455452580000135
步骤103,基于自适应策略原理,根据模块化多电平换流器的工况自适应切换多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真,所述工况为稳态工况或暂态工况。
具体步骤为:当模块化多电平换流器的工况为稳态工况时,采用工频周期作为多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真。
当模块化多电平换流器的工况为暂态工况时,采用工频周期内每种电平持续时间的平均时间作为多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真。具体的,判断滑动窗向前的过程中是否存在电平阶跃,得到判断结果。若所述判断结果表示所述滑动窗向前的过程中存在电平阶跃,则对模块化多电平换流器的有电平阶跃的开关函数:
Figure RE-GDA0002455452580000136
进行傅里叶分解,得到有电平阶跃的开关函数在频域内的表达式:
Figure RE-GDA0002455452580000137
其中,ts表示平均时间,n表示平均时间内电平的数量,t1、t2、tn和tn+1分别表示第一种电平、第二种电平、第n种电平和第n+1种电平的结束时间,t表示当前时间,i表示第一种电平的闭合开关的数量,A0表示对有电平阶跃的开关函数进行傅里叶分解时得到的直流量、ck对有电平阶跃的开关函数进行傅里叶分解得到的第k阶信号分量的谐波系数,ω2表示平均时间尺度对应的频率、k表示第k阶信号分量、δx表示x相的动态相位和
Figure RE-GDA0002455452580000141
表示x相的固定相位;利用有电平阶跃的开关函数在频域内的表达式,对所述多频段动态向量模型进行仿真。若所述判断结果表示所述滑动窗向前的过程中不存在电平阶跃,对模块化多电平换流器的无电平阶跃的开关函数:s(ts)=lusm-x(l<N+1),进行傅里叶分解,得到无电平阶跃的开关函数在频域内的表达式:
Figure RE-GDA0002455452580000142
其中, l为滑动窗口内电平的数量;利用无电平阶跃的开关函数在频域内的表达式,对所述多频段动态向量模型进行仿真。
下面对步骤103进行详细说明。在MMC多频段动态相量模型的基础上,加入自适应策略,灵活改变开关函数在傅里叶分解时“滑动窗”的宽度,此时能够及时捕捉暂态波形的快速变化,提高暂态仿真过程的精度,得到MMC 自适应多频段动态相量模型。
暂态工况下仿真过程中的核心为合理调整傅里叶分解滑动窗宽的时间尺度,在传统两电平VSC的PWM波调制的换流器中,一般可直接将开关周期设置为窗宽。MMC换流器输出波形由最近电平逼近控制,根据MMC输出阶梯波的波形特点,单个桥臂子模块数量为N,阶梯波的输出电平数与桥臂子模块数量的关系为N+1,但每种电平持续的时间长短不同,与正弦调制波幅值变化的速率有关。可将傅里叶分解窗宽定为在工频周期Ts内每种电平数持续的平均时间ts=Ts/4(N/2+1),并将其定义为单个平均开关周期。单个平均开关周期ts的大小如图7所示。
在单个平均开关周期内对阶梯波进行傅里叶展开,滑动窗向前的过程中,会出现无电平阶跃和有电平阶跃的情况,如图8所示在圆圈所代表的3个平均开关周期ts内,出现了0电平阶跃、1个电平阶跃和2个电平阶跃的情况。图 8中,平滑的曲线表示调制波,锯齿形曲线表示NLM波。
针对有无电平阶跃的情况分别进行傅里叶分解:
1、在平均开关周期ts内无电平阶跃时的开关函数时域表达式为:
s(ts)=iusm-x(i<N+1) (17)
其中,usm-x为每个桥臂子模块的平均电压,在式(14)的开关函数模型中代表一个单位电平,即usm-x=1;i为滑动窗口内电平的数量该窗口内电平数目。
s(ts)=i(i<N+1) (18)
由傅里叶分解得到频域内的解析式为:
s(ts)=i(i<N+1) (19)
(2)在平均开关周期ts内有电平阶跃时,此时有两种情况,第一种为电平数增加的阶跃,第二种为电平数减少的阶跃,两种情况下的开关函数的时域表达式为分段函数,综合表述为:
Figure RE-GDA0002455452580000151
其中,t表示每个电平持续的时间长度,tn表示每个电平值持续的时间始末, n表示阶跃的个数。
以电平数增加的阶跃情况为例,由傅里叶分解得到频域内的解析式为:
Figure RE-GDA0002455452580000152
其中
Figure RE-GDA0002455452580000153
Figure RE-GDA0002455452580000154
Figure RE-GDA0002455452580000155
两种情况的开关函数频域表达式形式相同,系数A0,ak,bk不同。出现的阶跃数n不同,则A0,ak,bk均有差异。由此得到新的随平均开关周期变化的开关函数在频域内的表达式为:
(1)在平均开关周期ts内无电平阶跃时
Figure RE-GDA0002455452580000161
(2)在平均开关周期ts内有电平阶跃时:
Figure RE-GDA0002455452580000162
对于一个时域信号x(t)的自适应多频段动态相量形式可表示为:
Figure RE-GDA0002455452580000163
其中滑动窗宽的取值为:
Figure RE-GDA0002455452580000164
在MMC多频段动态相量解耦模型的基础上,得到MMC自适应多频段动态向量模型如下:
Figure RE-GDA0002455452580000165
Figure RE-GDA0002455452580000166
Figure RE-GDA0002455452580000167
Figure RE-GDA0002455452580000168
与式(15)的MMC多频段动态相量解耦模型对比可以发现,由于式(27) 中滑动窗宽的取值τ发生改变,所以会导致仿真步长发生变化,从而导致计算过程中的导纳阵发生改变,但是两种工况下的仿真步长是确定的,故两种导纳阵可以确定,所以最后的计算量并没有增加太多。出现扰动时,τ取平均开关周期ts进行傅里叶分解得到开关函数表达式,缩小仿真步长,增大了暂态工况时的仿真精度。
本发明还提供一种模块化多电平换流器的动态仿真模拟系统,所述模拟系统包括:
开关函数模型建立模块,用于根据模块化多电平换流器的拓扑结构,建立模块化多电平换流器的开关函数模型;
所述开关函数模型建立模块,具体包括:开关函数模型建立子模块,用于根据模块化多电平换流器的拓扑结构,建立模块化多电平换流器的开关函数模型为:
Figure RE-GDA0002455452580000171
Figure RE-GDA0002455452580000172
Figure RE-GDA0002455452580000173
Figure RE-GDA0002455452580000174
其中,Lx表示模块化多电平换流器的第x相的桥臂电抗,Rx模块化多电平换流器的第x相的桥臂电阻,x为A、B或C,
Figure RE-GDA0002455452580000175
表示模块化多电平换流器的第x相的上桥臂与下桥臂的电流之和,
Figure RE-GDA0002455452580000176
表示模块化多电平换流器的第x相的上桥臂与下桥臂的电流之差,
Figure RE-GDA0002455452580000177
表示模块化多电平换流器的第x相的上桥臂与下桥臂的模块电压之和和
Figure RE-GDA0002455452580000178
表示模块化多电平换流器的第x相的上桥臂与下桥臂的模块电压之差,N表示每个桥臂的子模块的数量,Csm表示模块化多电平换流器的每个模块的电容,
Figure RE-GDA0002455452580000179
Figure RE-GDA00024554525800001710
分别表示模块化多电平换流器的第x 相的上桥臂与下桥臂的开关函数之和和开关函数之差,ux为模块化多电平换流器的第x相的桥臂、udc为模块化多电平换流器的直流侧电压。
多频段动态向量模型获取模块,用于基于多频段动态向量法原理,对所述开关函数模型进行分解、重组和移频获得模块化多电平换流器的多频段动态向量模型。
所述多频段动态向量模型获取模块,具体包括:
多频段动态向量模型获取子模块,用于基于多频段动态向量法原理,对所述开关函数模型进行分解、重组和移频获得模块化多电平换流器的多频段动态向量模型为:
Figure RE-GDA00024554525800001711
Figure RE-GDA0002455452580000181
Figure RE-GDA0002455452580000182
Figure RE-GDA0002455452580000183
其中,
Figure RE-GDA0002455452580000184
表示子频段动态向量,fh为第h阶谐波的频率,Bk表示第k个频段信号重组后的子频段信号,Rx表示模块化多电平换流器的第x相的桥臂电阻,h表示第h阶谐波,ωk表示第k个频段的所移频率,j表示虚数单位。
自适应仿真模块,用于基于自适应策略原理,根据模块化多电平换流器的工况自适应切换多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真,所述工况为稳态工况或暂态工况。
所述自适应仿真模块,具体包括:
第一自适应仿真子模块,用于当模块化多电平换流器的工况为稳态工况时,采用工频周期作为多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真。
第二自适应仿真子模块,当模块化多电平换流器的工况为暂态工况时,采用工频周期内每种电平持续时间的平均时间作为多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真。
所述第二自适应仿真子模块,具体包括:
判断单元,用于判断滑动窗向前的过程中是否存在电平阶跃,得到判断结果。
第一傅里叶分解单元,用于若所述判断结果表示所述滑动窗向前的过程中存在电平阶跃,则对模块化多电平换流器的有电平阶跃的开关函数:
Figure RE-GDA0002455452580000185
进行傅里叶分解,得到有电平阶跃的开关函数在频域内的表达式:
Figure RE-GDA0002455452580000186
其中,ts表示平均时间,n表示平均时间内电平的数量,t1、t2、tn和tn+1分别表示第一种电平、第二种电平、第n种电平和第n+1种电平的结束时间,t表示当前时间,i表示第一种电平的闭合开关的数量,A0表示对有电平阶跃的开关函数进行傅里叶分解时得到的直流量、ck对有电平阶跃的开关函数进行傅里叶分解得到的第k阶信号分量的谐波系数,ω2表示平均时间尺度对应的频率、k表示第k阶信号分量、δx表示x相的动态相位和
Figure RE-GDA0002455452580000191
表示x相的固定相位。
第一仿真单元,用于利用有电平阶跃的开关函数在频域内的表达式,对所述多频段动态向量模型进行仿真。
第二傅里叶分解单元,用于若所述判断结果表示所述滑动窗向前的过程中不存在电平阶跃,对模块化多电平换流器的无电平阶跃的开关函数: s(ts)=lusm-x(l<N+1),进行傅里叶分解,得到无电平阶跃的开关函数在频域内的表达式:
Figure RE-GDA0002455452580000192
其中,l为滑动窗口内电平的数量。
第一仿真单元,用于利用无电平阶跃的开关函数在频域内的表达式,对所述多频段动态向量模型进行仿真。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提出了一种模块化多电平换流器的动态仿真模拟方法。本发明提出了一种模块化多电平换流器的动态仿真模拟方法。本发明的方法首先建立模块化多电平换流器的开关函数模型;对所述开关函数模型进行分解、重组和移频获得模块化多电平换流器的多频段动态向量模型;然后,基于自适应策略原理,根据模块化多电平换流器的工况自适应切换多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真,所述工况为稳态工况或暂态工况。本发明通过频率的分解重组和移频,并采用自适应方法,进行模块化多电平换流器的动态仿真模拟,在保证精度的前提下,增大步长提高仿真速度。
具体的,本发明的模块化多电平换流器的动态仿真模拟方法与现有的仿真模型对比,表现出如下的优势:
1)与MMC详细模型相比:不用对每个子模块建模,不用通过小步长来精确捕捉高频开关的变化时刻,随着MMC子模块数目的增长,桥臂电平数较高时,计算效率也不会太低,因此可进行大规模仿真,解决详细模型仿真速度受限的问题。
2)与戴维南等值模型相比:MMC自适应多频段动态相量模型不需要插值回退抑制数值振荡,能够考虑在考虑更高频次的信号分量保证精度的同时,降低方程阶数减小仿真规模,此外还能够通过移频特性将高频信号转化为低频信号,增大仿真步长。
3)传统动态相量模型相比:随着考虑的谐波次数的增加,不用为减少计算量采取谐波截断的方式,使仿真精度降低。MMC自适应多频段动态相量模型按频段对信号进行分频和重组,在考虑更大频率带宽和更多频率成分的同时,又大为减少系统级微分方程组数目,能兼顾仿真精度和速度的要求。
4)多频段动态相量模型相比:MMC自适应多频段动态相量模型在MMC 多频段动态相量模型的基础上,增加自适应策略,缩小加窗傅里叶分解的滑动窗宽,提高暂态仿真过程的精度,稳态过程中,恢复工频周期的“滑动窗”分解,保持稳态的仿真步长,保持原来的仿真速度。
本说明书中等效实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,等效实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例,基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。

Claims (10)

1.一种模块化多电平换流器的动态仿真模拟方法,其特征在于,所述模拟方法包括如下步骤:
根据模块化多电平换流器的拓扑结构,建立模块化多电平换流器的开关函数模型;
基于多频段动态向量法原理,对所述开关函数模型进行分解、重组和移频获得模块化多电平换流器的多频段动态向量模型;
基于自适应策略原理,根据模块化多电平换流器的工况自适应切换多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真,所述工况为稳态工况或暂态工况。
2.根据权利要求1所述的模块化多电平换流器的动态仿真模拟方法,其特征在于,所述根据模块化多电平换流器的拓扑结构,建立模块化多电平换流器的开关函数模型,具体包括:
根据模块化多电平换流器的拓扑结构,建立模块化多电平换流器的开关函数模型为:
Figure FDA0002381626210000011
Figure FDA0002381626210000012
Figure FDA0002381626210000013
Figure FDA0002381626210000014
其中,Lx表示模块化多电平换流器的第x相的桥臂电抗,Rx模块化多电平换流器的第x相的桥臂电阻,x为A、B或C,
Figure FDA0002381626210000015
表示模块化多电平换流器的第x相的上桥臂与下桥臂的电流之和,
Figure FDA0002381626210000016
表示模块化多电平换流器的第x相的上桥臂与下桥臂的电流之差,
Figure FDA0002381626210000017
表示模块化多电平换流器的第x相的上桥臂与下桥臂的模块电压之和和
Figure FDA0002381626210000018
表示模块化多电平换流器的第x相的上桥臂与下桥臂的模块电压之差,N表示每个桥臂的子模块的数量,Csm表示模块化多电平换流器的每个模块的电容,
Figure FDA0002381626210000019
Figure FDA00023816262100000110
分别表示模块化多电平换流器的第x相的上桥臂与下桥臂的开关函数之和和开关函数之差,ux为模块化多电平换流器的第x相的桥臂、udc为模块化多电平换流器的直流侧电压。
3.根据权利要求2所述的模块化多电平换流器的动态仿真模拟方法,其特征在于,所述基于多频段动态向量法原理,对所述开关函数模型进行分解、重组和移频获得模块化多电平换流器的多频段动态向量模型,具体包括:
基于多频段动态向量法原理,对所述开关函数模型进行分解、重组和移频获得模块化多电平换流器的多频段动态向量模型为:
Figure FDA0002381626210000021
Figure FDA0002381626210000022
Figure FDA0002381626210000023
Figure FDA0002381626210000024
其中,
Figure FDA0002381626210000025
表示子频段动态向量,fh为第h阶谐波的频率,Bk表示第k个频段信号重组后的子频段信号,Rx表示模块化多电平换流器的第x相的桥臂电阻,h表示第h阶谐波,ωk表示第k个频段的所移频率,j表示虚数的单位。
4.根据权利要求3所述的模块化多电平换流器的动态仿真模拟方法,其特征在于,所述基于自适应策略原理,根据模块化多电平换流器的工况自适应切换多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真,具体包括:
当模块化多电平换流器的工况为稳态工况时,采用工频周期作为多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真;
当模块化多电平换流器的工况为暂态工况时,采用工频周期内每种电平持续时间的平均时间作为多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真。
5.根据权利要求4所述的模块化多电平换流器的动态仿真模拟方法,其特征在于,所述当模块化多电平换流器的工况为暂态工况时,采用工频周期内每种电平持续时间的平均时间平作为多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真,具体包括:
判断滑动窗向前的过程中是否存在电平阶跃,得到判断结果;
若所述判断结果表示所述滑动窗向前的过程中存在电平阶跃,则对模块化多电平换流器的有电平阶跃的开关函数:
Figure FDA0002381626210000031
进行傅里叶分解,得到有电平阶跃的开关函数在频域内的表达式:
Figure FDA0002381626210000032
其中,ts表示平均时间,n表示平均时间内电平的数量,t1、t2、tn和tn+1分别表示第一种电平、第二种电平、第n种电平和第n+1种电平的结束时间,t表示当前时间,i表示第一种电平的闭合开关的数量,A0表示对有电平阶跃的开关函数进行傅里叶分解时得到的直流量、ck对有电平阶跃的开关函数进行傅里叶分解得到的第k阶信号分量的谐波系数,ω2表示平均时间尺度对应的频率、k表示第k阶信号分量、δx表示x相的动态相位和
Figure FDA0002381626210000033
表示x相的固定相位;
利用有电平阶跃的开关函数在频域内的表达式,对所述多频段动态向量模型进行仿真;
若所述判断结果表示所述滑动窗向前的过程中不存在电平阶跃,对模块化多电平换流器的无电平阶跃的开关函数:s(ts)=lusm-x(l<N+1),进行傅里叶分解,得到无电平阶跃的开关函数在频域内的表达式:
Figure FDA0002381626210000034
其中,l为滑动窗口内电平的数量;
利用无电平阶跃的开关函数在频域内的表达式,对所述多频段动态向量模型进行仿真。
6.一种模块化多电平换流器的动态仿真模拟系统,其特征在于,所述模拟系统包括:
开关函数模型建立模块,用于根据模块化多电平换流器的拓扑结构,建立模块化多电平换流器的开关函数模型;
多频段动态向量模型获取模块,用于基于多频段动态向量法原理,对所述开关函数模型进行分解、重组和移频获得模块化多电平换流器的多频段动态向量模型;
自适应仿真模块,用于基于自适应策略原理,根据模块化多电平换流器的工况自适应切换多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真,所述工况为稳态工况或暂态工况。
7.根据权利要求6所述的模块化多电平换流器的动态仿真模拟系统,其特征在于,所述开关函数模型建立模块,具体包括:
开关函数模型建立子模块,用于根据模块化多电平换流器的拓扑结构,建立模块化多电平换流器的开关函数模型为:
Figure FDA0002381626210000041
Figure FDA0002381626210000042
Figure FDA0002381626210000043
Figure FDA0002381626210000044
其中,Lx表示模块化多电平换流器的第x相的桥臂电抗,Rx模块化多电平换流器的第x相的桥臂电阻,x为A、B或C,
Figure FDA0002381626210000045
表示模块化多电平换流器的第x相的上桥臂与下桥臂的电流之和,
Figure FDA0002381626210000046
表示模块化多电平换流器的第x相的上桥臂与下桥臂的电流之差,
Figure FDA0002381626210000047
表示模块化多电平换流器的第x相的上桥臂与下桥臂的模块电压之和和
Figure FDA0002381626210000048
表示模块化多电平换流器的第x相的上桥臂与下桥臂的模块电压之差,N表示每个桥臂的子模块的数量,Csm表示模块化多电平换流器的每个模块的电容,
Figure FDA0002381626210000049
Figure FDA00023816262100000410
分别表示模块化多电平换流器的第x相的上桥臂与下桥臂的开关函数之和和开关函数之差,ux为模块化多电平换流器的第x相的桥臂、udc为模块化多电平换流器的直流侧电压。
8.根据权利要求7所述的模块化多电平换流器的动态仿真模拟系统,其特征在于,所述多频段动态向量模型获取模块,具体包括:
多频段动态向量模型获取子模块,用于基于多频段动态向量法原理,对所述开关函数模型进行分解、重组和移频获得模块化多电平换流器的多频段动态向量模型为:
Figure FDA0002381626210000051
Figure FDA0002381626210000052
Figure FDA0002381626210000053
Figure FDA0002381626210000054
其中,
Figure FDA0002381626210000055
表示子频段动态向量,fh为第h阶谐波的频率,Bk表示第k个频段信号重组后的子频段信号,Rx表示模块化多电平换流器的第x相的桥臂电阻,h表示第h阶谐波,ωk表示第k个频段的所移频率,j表示虚数单位。
9.根据权利要求8所述的模块化多电平换流器的动态仿真模拟系统,其特征在于,所述自适应仿真模块,具体包括:
第一自适应仿真子模块,用于当模块化多电平换流器的工况为稳态工况时,采用工频周期作为多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真;
第二自适应仿真子模块,当模块化多电平换流器的工况为暂态工况时,采用工频周期内每种电平持续时间的平均时间作为多频段动态相量模型的傅里叶分解滑动窗宽,对所述多频段动态向量模型进行仿真。
10.根据权利要求9所述的模块化多电平换流器的动态仿真模拟系统,其特征在于,所述第二自适应仿真子模块,具体包括:
判断单元,用于判断滑动窗向前的过程中是否存在电平阶跃,得到判断结果;
第一傅里叶分解单元,用于若所述判断结果表示所述滑动窗向前的过程中存在电平阶跃,则对模块化多电平换流器的有电平阶跃的开关函数:
Figure FDA0002381626210000056
进行傅里叶分解,得到有电平阶跃的开关函数在频域内的表达式:
Figure FDA0002381626210000061
其中,ts表示平均时间,n表示平均时间内电平的数量,t1、t2、tn和tn+1分别表示第一种电平、第二种电平、第n种电平和第n+1种电平的结束时间,t表示当前时间,i表示第一种电平的闭合开关的数量,A0表示对有电平阶跃的开关函数进行傅里叶分解时得到的直流量、ck对有电平阶跃的开关函数进行傅里叶分解得到的第k阶信号分量的谐波系数,ω2表示平均时间尺度对应的频率、k表示第k阶信号分量、δx表示x相的动态相位和
Figure FDA0002381626210000062
表示x相的固定相位;
第一仿真单元,用于利用有电平阶跃的开关函数在频域内的表达式,对所述多频段动态向量模型进行仿真;
第二傅里叶分解单元,用于若所述判断结果表示所述滑动窗向前的过程中不存在电平阶跃,对模块化多电平换流器的无电平阶跃的开关函数:s(ts)=lusm-x(l<N+1),进行傅里叶分解,得到无电平阶跃的开关函数在频域内的表达式:
Figure FDA0002381626210000063
其中,l为滑动窗口内电平的数量;
第一仿真单元,用于利用无电平阶跃的开关函数在频域内的表达式,对所述多频段动态向量模型进行仿真。
CN202010084671.6A 2020-02-10 2020-02-10 一种模块化多电平换流器的动态仿真模拟方法及系统 Active CN111310325B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010084671.6A CN111310325B (zh) 2020-02-10 2020-02-10 一种模块化多电平换流器的动态仿真模拟方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010084671.6A CN111310325B (zh) 2020-02-10 2020-02-10 一种模块化多电平换流器的动态仿真模拟方法及系统

Publications (2)

Publication Number Publication Date
CN111310325A true CN111310325A (zh) 2020-06-19
CN111310325B CN111310325B (zh) 2022-01-11

Family

ID=71146918

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010084671.6A Active CN111310325B (zh) 2020-02-10 2020-02-10 一种模块化多电平换流器的动态仿真模拟方法及系统

Country Status (1)

Country Link
CN (1) CN111310325B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111969639A (zh) * 2020-08-17 2020-11-20 中国科学院电工研究所 电力电子化电网级联型换流器多时间尺度暂态建模方法
CN112103982A (zh) * 2020-08-27 2020-12-18 中国南方电网有限责任公司超高压输电公司检修试验中心 基于傅里叶分解的mmc小信号阻抗建模方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170132337A1 (en) * 2015-11-10 2017-05-11 Rama B. Nath Apparatus and method for modelling a modular multilevel converter in an electronic simulator
CN109241678A (zh) * 2018-10-22 2019-01-18 清华大学 模块化多电平换流器多拓扑通用仿真模型
CN109861574A (zh) * 2019-03-13 2019-06-07 东北电力大学 一种融合自适应环流抑制的无模型自适应迭代学习mmc变流器控制方法
CN110362937A (zh) * 2019-07-19 2019-10-22 华北电力大学 一种模块化多电平换流器电磁暂态仿真方法及系统
CN110378022A (zh) * 2019-07-19 2019-10-25 华北电力大学 一种多频段动态相量电磁暂态仿真方法及系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170132337A1 (en) * 2015-11-10 2017-05-11 Rama B. Nath Apparatus and method for modelling a modular multilevel converter in an electronic simulator
CN109241678A (zh) * 2018-10-22 2019-01-18 清华大学 模块化多电平换流器多拓扑通用仿真模型
CN109861574A (zh) * 2019-03-13 2019-06-07 东北电力大学 一种融合自适应环流抑制的无模型自适应迭代学习mmc变流器控制方法
CN110362937A (zh) * 2019-07-19 2019-10-22 华北电力大学 一种模块化多电平换流器电磁暂态仿真方法及系统
CN110378022A (zh) * 2019-07-19 2019-10-25 华北电力大学 一种多频段动态相量电磁暂态仿真方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
EHSAN KARAMI,ET AL.: "Dynamic Phasor-Based Analysis of Unbalanced Three-Phase Systems in Presence of Harmonic Distortion", 《IEEE TRANSACTIONS ON POWER SYSTEMS》 *
刘博宁 等: "多频段–动态相量法电磁暂态仿真研究", 《中国电机工程学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111969639A (zh) * 2020-08-17 2020-11-20 中国科学院电工研究所 电力电子化电网级联型换流器多时间尺度暂态建模方法
CN112103982A (zh) * 2020-08-27 2020-12-18 中国南方电网有限责任公司超高压输电公司检修试验中心 基于傅里叶分解的mmc小信号阻抗建模方法

Also Published As

Publication number Publication date
CN111310325B (zh) 2022-01-11

Similar Documents

Publication Publication Date Title
CN110362937B (zh) 一种模块化多电平换流器电磁暂态仿真方法及系统
CN110556831B (zh) 一种多机多节点发电系统的稳定性评估方法及系统
CN110348161B (zh) 电压源型换流器多频段动态相量电磁暂态仿真方法及系统
CN103472282B (zh) 一种基于自适应原理的fbd谐波电流检测方法
CN111310325B (zh) 一种模块化多电平换流器的动态仿真模拟方法及系统
CN110045186B (zh) 电力系统公共连接点处的系统谐波阻抗估计方法和系统
CN111969639B (zh) 电力电子化电网级联型换流器多时间尺度暂态建模方法
CN110378020B (zh) 电网换相换流器多频段动态相量电磁暂态仿真方法及系统
Varricchio et al. Electrical network dynamic models with application to modal analysis of harmonics
CN112737381B (zh) 一种基于环流压降平衡的mmc电容电压波动抑制方法
CN108448631A (zh) 基于模型降阶的含分布式电源接入配电网的动态相量建模方法
CN110378022B (zh) 一种多频段动态相量电磁暂态仿真方法及系统
CN110889210A (zh) 基于根匹配法的移频暂态仿真方法、系统、介质和设备
CN103424621A (zh) 一种谐波电流的人工神经网络检测方法
Smith et al. Harmonic tensor linearisation of HVDC converters
Ge et al. A dynamic parameter model of harmonic source networks
CN110137968A (zh) 一种包含vsc的谐波潮流计算方法
Turner et al. Computation of AC-DC System disturbances. PT. II-derivation of power frequency variables from convertor transient response
Konara et al. Interfacing electromagnetic transient simulation to transient stability model using a multi-port dynamic phasor buffer zone
CN111628495A (zh) 基于谐波相量法的含vsc-mmc交直流电网系统的电磁暂态仿真方法
Orillaza et al. Reduced harmonic state space model of TCR
CN116609610A (zh) 一种基于自适应小波包变换的单相接地故障选相方法
CN111259565B (zh) 一种电压源型换流器的动态仿真模拟方法及系统
Kwon et al. Analysis of harmonic coupling and stability in back-to-back converter systems for wind turbines using Harmonic State Space (HSS)
Jia et al. A novel adaptive linear prediction-based parameter estimation method for monitoring sub-/inter-harmonics during SSI events

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Yao Yifan

Inventor after: Yao Shujun

Inventor after: Pang Bohan

Inventor after: Qu Qiumeng

Inventor after: Zhang Xiangyu

Inventor after: Wang Yan

Inventor after: Liu Jin

Inventor before: Yao Shujun

Inventor before: Qu Qiumeng

Inventor before: Wang Yan

Inventor before: Liu Jin

GR01 Patent grant
GR01 Patent grant