CN110334476B - 一种电磁暂态仿真方法及系统 - Google Patents
一种电磁暂态仿真方法及系统 Download PDFInfo
- Publication number
- CN110334476B CN110334476B CN201910653213.7A CN201910653213A CN110334476B CN 110334476 B CN110334476 B CN 110334476B CN 201910653213 A CN201910653213 A CN 201910653213A CN 110334476 B CN110334476 B CN 110334476B
- Authority
- CN
- China
- Prior art keywords
- signal
- component
- frequency
- detected
- signals
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开一种电磁暂态仿真方法及系统。该方法包括:获取谐波电压源的电压信号;采用经验模态分解法对电压信号进行时间尺度分解,得到多个单分量信号;计算各单分量信号的瞬时频率,得到各分量频率;获取预划分的多个连续频段;按各个分量频率所属的频段对各单分量信号进行重组,得到各频段的单分量信号集合;将属于同一单分量信号集合的所有单分量信号进行合并,得到多个合并信号;对每个频段内的合并信号进行时间尺度变换,得到低频复数信号;将低频复数信号代入电网的微分方程,得到电网各节点仿真信号。本发明的电磁暂态仿真方法及系统能够提高仿真速度和仿真精度。
Description
技术领域
本发明涉及电力系统领域,特别是涉及一种电磁暂态仿真方法及系统。
背景技术
随着高压直流输电、柔性交流输电技术以及新能源的快速发展,我国电网逐渐形成大规模交直流混联的输电系统格局,电网特性复杂,短时间内大量电力电子装置的动作与大电网中长期的动态过程紧密相连,导致电网分析变得更加复杂化。
对于电力系统,电磁暂态是最重要的系统分析手段,但是随着电力电子设备在电力系统中应用,电磁暂态仿真面临着仿真速度和仿真精度之间的矛盾,例如现有技术中基于状态空间平均理论提出的动态相量法(Dynamic Phasor,DP)建立的动态相量模型用于仿真时,只考虑了较低的几次特征谐波,导致仿真精度不高,若是考虑谐波次数增加,会出现求解方程规模增大的问题,导致计算速度降低。因此现有技术普遍存在仿真速度和仿真精度无法兼顾的问题。
发明内容
本发明的目的是提供一种电磁暂态仿真方法及系统,提高仿真速度和仿真精度。
为实现上述目的,本发明提供了如下方案:
一种电磁暂态仿真方法,包括:
获取谐波电压源的电压信号;
采用经验模态分解法对所述电压信号进行时间尺度分解,得到多个单分量信号;
计算各单分量信号的瞬时频率,得到各分量频率;
获取预划分的多个连续频段;
按各个所述分量频率所属的频段对各所述单分量信号进行重组,得到各频段的单分量信号集合;
将属于同一所述单分量信号集合的所有单分量信号进行合并,得到多个合并信号;
对每个频段内的所述合并信号进行时间尺度变换,得到低频复数信号;
将所述低频复数信号代入电网的微分方程,得到电网各节点仿真信号。
可选的,所述采用经验模态分解法对所述电压信号进行时间尺度分解,得到多个单分量信号,具体包括:
将所述电压信号标记为原始信号;
根据所述原始信号的局部极大值点和局部极小值点拟合所述原始信号的上包络线和下包络线;
根据公式h1=z(t)-m1计算得到待检测分量;其中h1为待检测分量;z(t)为原始信号,m1为上包络线与下包络线的均值;
判断所述待检测分量是否为固有模态函数分量,得到第一判断结果;
若所述第一判断结果表示否,则将所述待检测分量标记为原始信号并返回步骤“根据所述原始信号的局部极大值点和局部极小值点拟合所述原始信号的上包络线和下包络线”;
若所述第一判断结果表示是,则将所述待检测分量标记为单分量信号,并从原始信号中将所述待检测分量分离开,得到剩余信号;
根据所述剩余信号和所述单分量信号判断是否满足终止条件,得到第二判断结果;
若所述第二判断结果表示否,则将所述剩余信号标记为原始信号并返回步骤“根据所述原始信号的局部极大值点和局部极小值点拟合所述原始信号的上包络线和下包络线”;
若所述第二判断结果表示是,则终止时间尺度分解,得到分解得到的所有单分量信号和一个剩余信号。
可选的,所述计算各单分量信号的瞬时频率,得到各分量频率,具体包括:
以所述单分量信号希尔伯特变换结果作为虚部,以所述单分量信号作为实部,生成所述单分量信号的解析信号;
计算所述解析信号的虚部与实部的商的反正切值,得到所述解析信号的瞬时相位;
将所述瞬时相位对时间求导得到所述分量频率。
可选的,所述对每个频段内的所述合并信号进行时间尺度变换,得到低频复数信号,具体包括:
利用公式对所述合并信号进行时间尺度变换,得到低频信号;其中zuv为变换前的信号,xdq为变换后的信号,f1为变换前的信号频率,f2为变换后的信号频率,t为时间,fr为dq坐标系的旋转速度,且f2=f1-fr;
根据采样定理对所述低频信号进行采样,得到低频复数信号。
可选的,所述将所述低频复数信号代入电网的微分方程,得到电网各节点仿真信号,具体包括:
将所述微分方程转换到旋转坐标系下,得到旋转坐标系微分方程;
将所述低频复数信号代入所述旋转坐标系微分方程得到各频段求解结果;
将各频段求解结果转换回静止坐标系下进行叠加得到电网各节点仿真信号。
本发明还公开一种电磁暂态仿真系统,包括:
电压信号获取模块,用于获取谐波电压源的电压信号;
经验模态分解模块,用于采用经验模态分解法对所述电压信号进行时间尺度分解,得到多个单分量信号;
分量频率计算模块,用于计算各单分量信号的瞬时频率,得到各分量频率;
频段划分模块,用于获取预划分的多个连续频段;
重组模块,用于按各个所述分量频率所属的频段对各所述单分量信号进行重组,得到各频段的单分量信号集合;
合并模块,用于将属于同一所述单分量信号集合的所有单分量信号进行合并,得到多个合并信号;
变换模块,用于对每个频段内的所述合并信号进行时间尺度变换,得到低频复数信号;
代入计算模块,用于将所述低频复数信号代入电网的微分方程,得到电网各节点仿真信号。
可选的,所述经验模态分解模块包括:
电压信号标记单元,用于将所述电压信号标记为原始信号;
包络线拟合单元,用于根据所述原始信号的局部极大值点和局部极小值点拟合所述原始信号的上包络线和下包络线;
待检测分量计算单元,用于根据公式h1=z(t)-m1计算得到待检测分量;其中h1为待检测分量;z(t)为原始信号,m1为上包络线与下包络线的均值;
第一判断单元,用于判断所述待检测分量是否为固有模态函数分量,得到第一判断结果;
待检测分量标记及返回单元,用于若所述第一判断结果表示否,则将所述待检测分量标记为原始信号并返回所述包络线拟合单元;
待检测分量标记及分离单元,用于若所述第一判断结果表示是,则将所述待检测分量标记为单分量信号,并从原始信号中将所述待检测分量分离开,得到剩余信号;
第二判断单元,用于根据所述剩余信号和所述单分量信号判断是否满足终止条件,得到第二判断结果;
剩余信号标记及返回单元,用于若所述第二判断结果表示否,则将所述剩余信号标记为原始信号并返回所述包络线拟合单元;
终止单元,用于若所述第二判断结果表示是,则终止时间尺度分解,得到分解得到的所有单分量信号和一个剩余信号。
可选的,所述分量频率计算模块包括:
解析信号计算单元,用于以所述单分量信号希尔伯特变换结果作为虚部,以所述单分量信号作为实部,生成所述单分量信号的解析信号;
瞬时相位计算单元,用于计算所述解析信号的虚部与实部的商的反正切值,得到所述解析信号的瞬时相位;
分量频率计算单元,用于将所述瞬时相位对时间求导得到所述分量频率。
可选的,所述变换模块包括:
时间尺度变换单元,用于利用公式对所述合并信号进行时间尺度变换,得到低频信号;其中zuv为变换前的信号,xdq为变换后的信号,f1为变换前的信号频率,f2为变换后的信号频率,t为时间,fr为dq坐标系的旋转速度,且f2=f1-fr;
采用单元,用于根据采样定理对所述低频信号进行采样,得到低频复数信号。
可选的,所述代入计算模块包括:
坐标转换单元,用于将所述微分方程转换到旋转坐标系下,得到旋转坐标系微分方程;
代入单元,用于将所述低频复数信号代入所述旋转坐标系微分方程得到各频段求解结果;
坐标换回及叠加单元,用于将各频段求解结果转换回静止坐标系下进行叠加得到电网各节点仿真信号。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明的电磁暂态仿真方法及系统采用了时间尺度分解和分频段并行运算的方法对电磁暂态进行仿真,从而能够将分解后的各频带的信号转换为低频信号,由此能够支持大步长采样仿真,提高了仿真效率。同时,由于通过时间尺度分解将信号分解为多个频段的多个合并信号,实现各频段的信号的并行运算,提高计算速度。另外,通过时间尺度分解将信号分解成多个频段,各个频段之间单独运算,从而减少信号之间的耦合,提高运算速度。并且,本发明通过时间尺度分解将信号分解成多个频段并采用时间尺度变换生成低频信号,从而使本发明的方案能够对各个频段的信号进行仿真,打破了只能对低频信号仿真的局限,使得高次谐波也能被考虑在内,提高仿真精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例1电磁暂态仿真方法的方法流程图;
图2为本发明实施例1中对电压信号进行时间尺度分解的流程图;
图3为本发明实施例1电磁暂态仿真方法的经验模态分解结果图;
图4为本发明实施例1电磁暂态仿真方法的单分量信号重组示意图;
图5为本发明实施例1电磁暂态仿真方法的某个单分量信号集合的合并示意图;
图6为本发明电磁暂态仿真方法具体实施例的电网系统节点结构图;
图7为本发明实施例2电磁暂态仿真系统的系统结构图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种电磁暂态仿真方法及系统,提高仿真速度和仿真精度。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例1:
图1为本发明实施例1电磁暂态仿真方法的方法流程图。
参见图1,该电磁暂态仿真方法,包括:
步骤101:获取谐波电压源的电压信号。
步骤102:采用经验模态分解法对所述电压信号进行时间尺度分解,得到多个单分量信号。经验模态分解可以使本实施例适用于非线性非平稳信号。
图2为本发明实施例1中对电压信号进行时间尺度分解的流程图。
参见图2,该步骤102具体包括:
步骤201:将所述电压信号u(t)标记为原始信号。
步骤202:根据所述原始信号的局部极大值点和局部极小值点拟合所述原始信号的上包络线和下包络线。首先确定原始信号的所有局部极值点,然后用三次样条曲线将所有的局部极大值点连接起来形成上包络线,再用三次样条曲线将所有的局部极小值点连接起来形成下包络线。上下包络线形成后,所有的数据点均包含在上包络线与下包络线之间。
步骤203:根据公式h1=z(t)-m1计算得到待检测分量;其中h1为待检测分量;z(t)为原始信号,m1为上包络线与下包络线的均值。
步骤204:判断所述待检测分量是否为固有模态函数分量,得到第一判断结果。待检测分量属于固有模态函数分量需要满足的两个条件为:1、待检测分量的极值点和过零点数目相等或最多相差1个。2、待检测分量的上包络线和下包络线的局部均值为0。
步骤205:若所述第一判断结果表示否,则将所述待检测分量标记为原始信号并返回步骤202。
步骤206:若所述第一判断结果表示是,则将所述待检测分量标记为单分量信号,并从原始信号中将所述待检测分量分离开,得到剩余信号。
通过反复执行步骤202~206从而逐个将各个单分量信号分离出来。
步骤207:根据所述剩余信号和所述单分量信号判断是否满足终止条件,得到第二判断结果。终止条件为:剩余信号为单调函数或单分量信号的上下包络线的均值小于预设值。
步骤208:若所述第二判断结果表示否,则将所述剩余信号标记为原始信号并返回步骤202;
步骤209:若所述第二判断结果表示是,则终止时间尺度分解,得到分解得到的所有单分量信号和一个剩余信号。即其中h为各单分量信号的序号,M为单分量信号的数量,ch(t)为序号为h的单分量信号,反映电压信号在不同时间尺度下的特征,r(t)为剩余信号,表示电压信号的平均趋势。
图3为本发明实施例1电磁暂态仿真方法的经验模态分解结果图。
参见图3,经过经过经验模态分解后,电压信号被分解为不同频率的多个单分量信号。
步骤103:计算各单分量信号的瞬时频率,得到各分量频率。
该步骤103的具体过程包括:
以所述单分量信号希尔伯特变换结果作为虚部,以所述单分量信号作为实部,生成所述单分量信号的解析信号。计算公式为:
Ch(t)=ch(t)+jH[ch(t)]=ch(t)+jzh(t)
其中Ch(t)为序号为h的单分量信号的解析信号,H[ch(t)]表示对ch(t)进行希尔伯特变换。
计算所述解析信号的虚部与实部的商的反正切值,得到所述解析信号的瞬时相位。计算公式为:
其中,φh(t)为瞬时相位。
将所述瞬时相位对时间求导得到所述分量频率。计算公式为:
其中,ωh(t)为瞬时角频率,即所述分量频率。
步骤104:获取预划分的多个连续频段。
将预设的0~f的频率范围N等分,从而得到N个连续频段。通常情况下,频率上限f为2kHz。
步骤105:按各个所述分量频率所属的频段对各所述单分量信号进行重组,得到各频段的单分量信号集合。
图4为本发明实施例1电磁暂态仿真方法的单分量信号重组示意图。
参见图4,根据各单分量信号所处的频段,将各单分量信号划分到对应的频段中。属于同一频段的单分量信号构成一个单分量信号集合,即构成一组单分量信号。
步骤106:将属于同一所述单分量信号集合的所有单分量信号进行合并,得到多个合并信号;
图5为本发明实施例1电磁暂态仿真方法的某个单分量信号集合的合并示意图。
参见图5,将属于同一个单分量集合的各单分量信号按时刻相加。即将单分量集合内各单分量信号属于同一时刻的值相加。
步骤107:对每个频段内的所述合并信号进行时间尺度变换,得到低频复数信号;时间尺度变换的主导频率取各频段的中心频率。
该步骤107具体包括:
利用公式对所述合并信号进行时间尺度变换,得到低频信号;其中zuv为变换前的信号,xdq为变换后的信号,f1为变换前的信号频率,f2为变换后的信号频率,t为时间,fr为dq坐标系的旋转速度,且f2=f1-fr;
根据采样定理对所述低频信号进行采样,得到低频复数信号。
由于时间尺度变换后得到的信号为低频信号,因此可以支持大步长仿真,从而能够提高仿真效率。
步骤108:将所述低频复数信号代入电网的微分方程,得到电网各节点仿真信号。
该步骤108具体包括:
将所述微分方程转换到旋转坐标系下,得到旋转坐标系微分方程;
将所述低频复数信号代入所述旋转坐标系微分方程得到各频段求解结果;
将各频段求解结果转换回静止坐标系下进行叠加得到电网各节点仿真信号。
下面给出本发明的方法的一个具体实施例:
图6为本发明电磁暂态仿真方法具体实施例的电网系统节点结构图。
参见图6,该具体实施例的电网系统包括3个节点。电源为谐波电压源,仿真步长为50μs。
电压源的电压信号按照上述步骤101~107进行处理。其中步骤104的频段划分的具体频段数为10个频段。
在执行步骤108时,该电网系统的微分方程在静止坐标系下为:
将微分方程转换到旋转频率为ωr的旋转坐标系下为:
其中,u1、u2、u3分别为3个节点的电压。i12、i23和i3L分别为第一第二节点间电流、第二第三节点间电流和第三节点至末端间的电流。L12、L23和L3L分别为第一第二节点间电感、第二第三节点间电感和第三节点至末端间的电感。R23和R3L分别为第二第三节点间电阻和第三节点至末端间的电阻。i20和i30分别为第二节点的对地电流和第三节点的对地电流。C20和C30分别为第二节点的对地电容和第三节点的对地电容。
u1dq、u2dq、u3dq分别为旋转坐标系下3个节点的电压;i12dq、i23dq和i3Ldq分别为旋转坐标系下第一第二节点间电流、旋转坐标系下第二第三节点间电流和旋转坐标系下第三节点至末端间的电流。
然后在旋转坐标系下求解,并将各频段求解结果转换回静止坐标系下进行叠加得到电网各节点仿真信号。
实施例2:
图7为本发明实施例2电磁暂态仿真系统的系统结构图。
参见图7,该电磁暂态仿真系统,包括:
电压信号获取模块1001,用于获取谐波电压源的电压信号;
经验模态分解模块1002,用于采用经验模态分解法对所述电压信号进行时间尺度分解,得到多个单分量信号;
分量频率计算模块1003,用于计算各单分量信号的瞬时频率,得到各分量频率;
频段划分模块1004,用于获取预划分的多个连续频段;
重组模块1005,用于按各个所述分量频率所属的频段对各所述单分量信号进行重组,得到各频段的单分量信号集合;
合并模块1006,用于将属于同一所述单分量信号集合的所有单分量信号进行合并,得到多个合并信号;
变换模块1007,用于对每个频段内的所述合并信号进行时间尺度变换,得到低频复数信号;
代入计算模块1008,用于将所述低频复数信号代入电网的微分方程,得到电网各节点仿真信号。
可选的,所述经验模态分解模块1002包括:
电压信号标记单元,用于将所述电压信号标记为原始信号;
包络线拟合单元,用于根据所述原始信号的局部极大值点和局部极小值点拟合所述原始信号的上包络线和下包络线;
待检测分量计算单元,用于根据公式h1=z(t)-m1计算得到待检测分量;其中h1为待检测分量;z(t)为原始信号,m1为上包络线与下包络线的均值;
第一判断单元,用于判断所述待检测分量是否为固有模态函数分量,得到第一判断结果;
待检测分量标记及返回单元,用于若所述第一判断结果表示否,则将所述待检测分量标记为原始信号并返回包络线拟合单元;
待检测分量标记及分离单元,用于若所述第一判断结果表示是,则将所述待检测分量标记为单分量信号,并从原始信号中将所述待检测分量分离开,得到剩余信号;
第二判断单元,用于根据所述剩余信号和所述单分量信号判断是否满足终止条件,得到第二判断结果;
剩余信号标记及返回单元,用于若所述第二判断结果表示否,则将所述剩余信号标记为原始信号并返回包络线拟合单元;
终止单元,用于若所述第二判断结果表示是,则终止时间尺度分解,得到分解得到的所有单分量信号和一个剩余信号。
可选的,所述分量频率计算模块1003包括:
解析信号计算单元,用于以所述单分量信号希尔伯特变换结果作为虚部,以所述单分量信号作为实部,生成所述单分量信号的解析信号;
瞬时相位计算单元,用于计算所述解析信号的虚部与实部的商的反正切值,得到所述解析信号的瞬时相位;
分量频率计算单元,用于将所述瞬时相位对时间求导得到所述分量频率。
可选的,所述变换模块1007包括:
时间尺度变换单元,用于利用公式对所述合并信号进行时间尺度变换,得到低频信号;其中zuv为变换前的信号,xdq为变换后的信号,f1为变换前的信号频率,f2为变换后的信号频率,t为时间,fr为dq坐标系的旋转速度,且f2=f1-fr;
采用单元,用于根据采样定理对所述低频信号进行采样,得到低频复数信号。
可选的,所述代入计算模块1008包括:
坐标转换单元,用于将所述微分方程转换到旋转坐标系下,得到旋转坐标系微分方程;
代入单元,用于将所述低频复数信号代入所述旋转坐标系微分方程得到各频段求解结果;
坐标换回及叠加单元,用于将各频段求解结果转换回静止坐标系下进行叠加得到电网各节点仿真信号。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明的电磁暂态仿真方法及系统采用了时间尺度分解和分频段并行运算的方法对电磁暂态进行仿真,从而能够将分解后的各频带的信号转换为低频信号,由此能够支持大步长采样仿真,提高了仿真效率。同时,由于通过时间尺度分解将信号分解为多个频段的多个单分量信号,实现各频段的信号的并行运算,提高计算速度。另外,通过时间尺度分解将信号分解成多个频段,各个频段之间单独运算,从而减少信号之间的耦合,提高运算速度。并且,本发明通过时间尺度分解将信号分解成多个频段并采用时间尺度变换生成低频信号,从而使本发明的方案能够对各个频段的信号进行仿真,打破了只能对低频信号仿真的局限,使得高次谐波也能被考虑在内,提高仿真精度。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种电磁暂态仿真方法,其特征在于,包括:
获取谐波电压源的电压信号;
采用经验模态分解法对所述电压信号进行时间尺度分解,得到多个单分量信号;
计算各单分量信号的瞬时频率,得到各分量频率;
获取预划分的多个连续频段;
按各个所述分量频率所属的频段对各所述单分量信号进行重组,得到各频段的单分量信号集合;
根据各单分量信号所处的频段,将各单分量信号划分到对应的频段中,属于同一频段的单分量信号构成一个单分量信号集合,即构成一组单分量信号;
将属于同一所述单分量信号集合的所有单分量信号进行合并,得到多个合并信号;
将属于同一个单分量集合的各单分量信号按时刻相加,即将单分量集合内各单分量信号属于同一时刻的值相加;
对每个频段内的所述合并信号进行时间尺度变换,得到低频复数信号;
将所述低频复数信号代入电网的微分方程,得到电网各节点仿真信号;
所述将所述低频复数信号代入电网的微分方程,得到电网各节点仿真信号,具体包括:
将所述微分方程转换到旋转坐标系下,得到旋转坐标系微分方程;
将所述低频复数信号代入所述旋转坐标系微分方程得到各频段求解结果;
将各频段求解结果转换回静止坐标系下进行叠加得到电网各节点仿真信号。
2.根据权利要求1所述的电磁暂态仿真方法,其特征在于,所述采用经验模态分解法对所述电压信号进行时间尺度分解,得到多个单分量信号,具体包括:
将所述电压信号标记为原始信号;
根据所述原始信号的局部极大值点和局部极小值点拟合所述原始信号的上包络线和下包络线;
根据公式h1=z(t)-m1计算得到待检测分量;其中h1为待检测分量;z(t)为原始信号,m1为上包络线与下包络线的均值;
判断所述待检测分量是否为固有模态函数分量,得到第一判断结果;
若所述第一判断结果表示否,则将所述待检测分量标记为原始信号并返回步骤“根据所述原始信号的局部极大值点和局部极小值点拟合所述原始信号的上包络线和下包络线”;
若所述第一判断结果表示是,则将所述待检测分量标记为单分量信号,并从原始信号中将所述待检测分量分离开,得到剩余信号;
根据所述剩余信号和所述单分量信号判断是否满足终止条件,得到第二判断结果;
若所述第二判断结果表示否,则将所述剩余信号标记为原始信号并返回步骤“根据所述原始信号的局部极大值点和局部极小值点拟合所述原始信号的上包络线和下包络线”;
若所述第二判断结果表示是,则终止时间尺度分解,得到分解得到的所有单分量信号和一个剩余信号。
3.根据权利要求1所述的电磁暂态仿真方法,其特征在于,所述计算各单分量信号的瞬时频率,得到各分量频率,具体包括:
以所述单分量信号希尔伯特变换结果作为虚部,以所述单分量信号作为实部,生成所述单分量信号的解析信号;
计算所述解析信号的虚部与实部的商的反正切值,得到所述解析信号的瞬时相位;
将所述瞬时相位对时间求导得到所述分量频率。
5.一种电磁暂态仿真系统,其特征在于,包括:
电压信号获取模块,用于获取谐波电压源的电压信号;
经验模态分解模块,用于采用经验模态分解法对所述电压信号进行时间尺度分解,得到多个单分量信号;
分量频率计算模块,用于计算各单分量信号的瞬时频率,得到各分量频率;
频段划分模块,用于获取预划分的多个连续频段;
重组模块,用于按各个所述分量频率所属的频段对各所述单分量信号进行重组,得到各频段的单分量信号集合,具体的,所述重组模块根据各单分量信号所处的频段,将各单分量信号划分到对应的频段中,属于同一频段的单分量信号构成一个单分量信号集合,即构成一组单分量信号;
合并模块,用于将属于同一所述单分量信号集合的所有单分量信号进行合并,得到多个合并信号,具体的,所述合并模块将属于同一个单分量集合的各单分量信号按时刻相加,即将单分量集合内各单分量信号属于同一时刻的值相加;
变换模块,用于对每个频段内的所述合并信号进行时间尺度变换,得到低频复数信号;
代入计算模块,用于将所述低频复数信号代入电网的微分方程,得到电网各节点仿真信号;
所述代入计算模块包括:
坐标转换单元,用于将所述微分方程转换到旋转坐标系下,得到旋转坐标系微分方程;
代入单元,用于将所述低频复数信号代入所述旋转坐标系微分方程得到各频段求解结果;
坐标换回及叠加单元,用于将各频段求解结果转换回静止坐标系下进行叠加得到电网各节点仿真信号。
6.根据权利要求5所述的电磁暂态仿真系统,其特征在于,所述经验模态分解模块包括:
电压信号标记单元,用于将所述电压信号标记为原始信号;
包络线拟合单元,用于根据所述原始信号的局部极大值点和局部极小值点拟合所述原始信号的上包络线和下包络线;
待检测分量计算单元,用于根据公式h1=z(t)-m1计算得到待检测分量;其中h1为待检测分量;z(t)为原始信号,m1为上包络线与下包络线的均值;
第一判断单元,用于判断所述待检测分量是否为固有模态函数分量,得到第一判断结果;
待检测分量标记及返回单元,用于若所述第一判断结果表示否,则将所述待检测分量标记为原始信号并返回所述包络线拟合单元;
待检测分量标记及分离单元,用于若所述第一判断结果表示是,则将所述待检测分量标记为单分量信号,并从原始信号中将所述待检测分量分离开,得到剩余信号;
第二判断单元,用于根据所述剩余信号和所述单分量信号判断是否满足终止条件,得到第二判断结果;
剩余信号标记及返回单元,用于若所述第二判断结果表示否,则将所述剩余信号标记为原始信号并返回所述包络线拟合单元;
终止单元,用于若所述第二判断结果表示是,则终止时间尺度分解,得到分解得到的所有单分量信号和一个剩余信号。
7.根据权利要求5所述的电磁暂态仿真系统,其特征在于,所述分量频率计算模块包括:
解析信号计算单元,用于以所述单分量信号希尔伯特变换结果作为虚部,以所述单分量信号作为实部,生成所述单分量信号的解析信号;
瞬时相位计算单元,用于计算所述解析信号的虚部与实部的商的反正切值,得到所述解析信号的瞬时相位;
分量频率计算单元,用于将所述瞬时相位对时间求导得到所述分量频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910653213.7A CN110334476B (zh) | 2019-07-19 | 2019-07-19 | 一种电磁暂态仿真方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910653213.7A CN110334476B (zh) | 2019-07-19 | 2019-07-19 | 一种电磁暂态仿真方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110334476A CN110334476A (zh) | 2019-10-15 |
CN110334476B true CN110334476B (zh) | 2020-05-01 |
Family
ID=68145163
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910653213.7A Active CN110334476B (zh) | 2019-07-19 | 2019-07-19 | 一种电磁暂态仿真方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110334476B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111709209B (zh) * | 2020-06-16 | 2024-04-05 | 华北电力大学 | 基于支路指数积分形式的电磁暂态仿真方法及系统 |
CN112199914B (zh) * | 2020-09-28 | 2024-06-04 | 华北电力大学 | 一种电力电子开关恒导纳模型建立方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101727522A (zh) * | 2009-12-11 | 2010-06-09 | 中国电力科学研究院 | 一种基于出厂数据的电动机机电暂态仿真模型参数获取方法 |
CN103678827A (zh) * | 2013-12-30 | 2014-03-26 | 云南电力试验研究院(集团)有限公司电力研究院 | 一种用于逆变器的电磁暂态建模方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109815607B (zh) * | 2019-01-31 | 2022-10-11 | 上海交通大学 | 面向大规模交流电网基于时频坐标变换的全电磁暂态仿真方法 |
-
2019
- 2019-07-19 CN CN201910653213.7A patent/CN110334476B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101727522A (zh) * | 2009-12-11 | 2010-06-09 | 中国电力科学研究院 | 一种基于出厂数据的电动机机电暂态仿真模型参数获取方法 |
CN103678827A (zh) * | 2013-12-30 | 2014-03-26 | 云南电力试验研究院(集团)有限公司电力研究院 | 一种用于逆变器的电磁暂态建模方法 |
Non-Patent Citations (4)
Title |
---|
Fast Electromagnetic Transient Model of Modular Multilevel Converter Based on Rotation Transformation;Yao Shujun 等;《2018 International Conference on Power System Technology》;20181108;2571-2577 * |
Modified electromagnetic transient model of transmission line based on rotation transformation;Shujun Yao 等;《The Journal of Engineering》;20181213;1045-1050 * |
大规模电网电磁暂态快速仿真方法;姚蜀军 等;《电力建设》;20151231;第36卷(第12期);16-21 * |
宽频时间尺度变换多速率电磁暂态仿真研究;姚蜀军 等;《中国电机工程学报》;20190205;第39卷(第3期);675-684 * |
Also Published As
Publication number | Publication date |
---|---|
CN110334476A (zh) | 2019-10-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110334476B (zh) | 一种电磁暂态仿真方法及系统 | |
CN109030941A (zh) | Hanning自乘卷积窗FFT三谱线插值谐波分析方法 | |
CN109490630B (zh) | 一种基于矩阵束的动态相量测量方法 | |
CN110348161B (zh) | 电压源型换流器多频段动态相量电磁暂态仿真方法及系统 | |
CN106226723A (zh) | 一种面向复杂非线性负载环境的电能计量模拟分析系统及方法 | |
CN103217578A (zh) | 基于锁相环技术测量信号间相位差的数字相位计及其方法 | |
CN110378022B (zh) | 一种多频段动态相量电磁暂态仿真方法及系统 | |
CN106484962B (zh) | 一种谐振式无线输电系统稳态特性的符号分析方法 | |
CN114050568A (zh) | 基于多元经验模态分解的强迫振荡源定位方法及装置 | |
CN107064646B (zh) | 一种用于多输入多输出不对称系统的电网阻抗实时辨识方法 | |
CN106405464B (zh) | 一种产生可溯源的任意波形模拟功率信号的方法 | |
CN101806835B (zh) | 基于包络分解的间谐波测量仪 | |
CN103904653A (zh) | 电网谐波状态估计中的可观性分析方法 | |
CN110007147A (zh) | 一种单母线多用户的谐波责任计算方法 | |
CN110378021B (zh) | 一种输电线路仿真方法及系统 | |
Kovalenko et al. | Synchrophasor Evaluation based on point-on-wave Measurements | |
Li et al. | Harmonic detection algorithm based on Kaiser window | |
CN104483563A (zh) | 电力信号的同步采样方法及系统 | |
CN109726490B (zh) | 用于电力系统的低频段多正弦信号设计方法 | |
CN110135015B (zh) | 面向大规模交直流系统的宽频机电-电磁混合仿真方法 | |
Banerjee et al. | Detection and location of fault in microgrid using discrete wavelet transform based technique | |
CN106203883A (zh) | 一种基于区间的配电网小信号稳定性计算方法 | |
Khatua et al. | A Noninvasive Measurement Technique of Grid and Converter Wideband Impedance | |
CN105337279A (zh) | 微电网接入配电网时的谐波电流电压的分析方法和装置 | |
KR20130077679A (ko) | 계통 전압 계측에 dft를 이용하는 계통 연계형 인버터 제어기 및 그 제어기를 구비한 계통 연계형 인버터 |
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 |