CN116722563A - 一种基于动态相量的电磁暂态仿真频域扩展方法及装置 - Google Patents
一种基于动态相量的电磁暂态仿真频域扩展方法及装置 Download PDFInfo
- Publication number
- CN116722563A CN116722563A CN202310623078.8A CN202310623078A CN116722563A CN 116722563 A CN116722563 A CN 116722563A CN 202310623078 A CN202310623078 A CN 202310623078A CN 116722563 A CN116722563 A CN 116722563A
- Authority
- CN
- China
- Prior art keywords
- current
- equation
- voltage
- stator
- dynamic
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000004088 simulation Methods 0.000 title claims abstract description 44
- 230000001052 transient effect Effects 0.000 title claims abstract description 37
- 230000001360 synchronised effect Effects 0.000 claims abstract description 14
- 239000003990 capacitor Substances 0.000 claims description 13
- 238000004458 analytical method Methods 0.000 claims description 11
- 230000004907 flux Effects 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 7
- 241001183012 Modified Vaccinia Ankara virus Species 0.000 claims description 6
- 238000004804 winding Methods 0.000 claims description 6
- 238000013459 approach Methods 0.000 claims description 5
- 230000001939 inductive effect Effects 0.000 claims description 4
- 238000013016 damping Methods 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims description 3
- 230000004069 differentiation Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 230000005284 excitation Effects 0.000 claims description 3
- 230000010363 phase shift Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000010304 firing Methods 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000012966 insertion method Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 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
- H02J3/241—The oscillation concerning frequency
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- 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
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
一种基于动态相量的电磁暂态仿真频域扩展方法,用来构建连续微分代数方程(differential‑algebraic equations,DAE),改进节点法MNA利用了stamp技术,很大的提高了MNA的模拟器;推导出电力系统模型,提出了一种基于MNA同步电机的新模型,该模型将电机等效为非线性电感而不是次暂态的等效物,使用可变步长对DAE进行数值求解。以及提供一种基于动态相量的电磁暂态仿真频域扩展装置。本发明使用动态相量来提高电磁暂态仿真,动态相量是稳态相量的自然延伸,它可以动态进行建模,同时保留了调解载波频率的优点。
Description
技术领域
本发明属于电力检测领域,涉及一种基于动态相量的电磁暂态仿真频域扩展方法及装置。
背景技术
暂态稳定性是电力系统在母线故障扰动后保持稳定的能力。传统电网暂态稳定性一直在0Hz左右低频振荡。第一种暂态稳定性通过使用计算机导纳矩阵程序模拟电网,这假设称为准稳态。这种方法适用于早期电力系统的稳定性分析,因为它提高了高频网络特性,所以不适合低频的机电系统。去除网络的高频特性更加有利于大型电力系统模型进行长时间仿真。
现代电力系统一个主要问题是发电机涡轮机、串联电容器、高压直流输电系统和可再生能源发电相关的次同步振荡。高压直流输电系统可以控制频率大于10Hz的扭转模态。在风力发电串联补偿线路的系统中,频率大于20Hz的振荡会导致次同步谐振。频率的动态变化无法使用传统的暂态稳定性进行分析,因为它们不能使用准稳态模拟出电网的正确频率。
电磁暂态(EMT)模拟器通常与电力系统的准确性和详细性研究相关联。EMT仿真通常用来计算复杂的电力系统,因此大型系统模型应该采用并行处理的方法进行仿真。另一种方法是使用动态相量来增强暂态稳定性。动态相量是稳态相量的自然延伸,它可以进行动态的建模,同时保留了调解载波频率的优点。此外,扩展很多现有模型来适应动态相量的仿真。
传统暂态稳定性模拟器通过构建和求解微分代数的方法不能用于动态相量,因为它们本质上与准稳态相关联。可用现有的电路仿真技术开发一种新方法,因为它在瞬时电路结构上与动态相量相似。数值替代积分器(NIS)是EMT仿真的基础,在早先的研究中对基于动态相量的仿真进行了探索。NIS的主要缺点是将局部仿真步长时间嵌入到了系统方程中,这样导致仿真变得困难且效率低。
比如专利CN201910756536.9公开了适用于微电网实时仿真的混合电磁暂态仿真方法,该专利特点在于将传统的节点分析法(NAM)和高度并行化的延迟插入法(LIM)相结合。但是该专利采用传统的节点分析方法不能提高仿真频域,仿真速度比较缓慢,不利于大型电力系统模型进行长时间仿真。
发明内容
为了克服已有电磁暂态仿真频域扩展方式的不足,本发明提供了一种基于动态相量的电磁暂态仿真频域扩展方法及装置,使用动态相量来提高电磁暂态仿真,动态相量是稳态相量的自然延伸,它可以动态进行建模,同时保留了调解载波频率的优点。
本发明解决其技术问题所采用的技术方案是:
一种基于动态相量的电磁暂态仿真频域扩展方法,包括以下步骤:
步骤1、电力系统动态微分代数方程DAE(differential-algebraic equations,即连续微分代数方程)如下:
其中,t为时间,y是系统变量,是y的导数;
模拟大型电力系统模型中网络层的瞬变,使用平衡的三相动态相量框架,框架使用Blondel变换来定义动态相量算子给定正弦信号如下:
x(t)=X(t)cos(ω0t+δ(t) (2)
其中,ω0是频率,δ(t)是相移,cos是余弦函数,X(t)为傅里叶级数;
计算x(t)的动态相量
其中,是动态相量算子,j是虚函数的虚部。正弦对应的/>通过重新引入载体来恢复:
其中,Re是载体函数;
动态相量模型的导数定义为:
用单位时间表示动态相量微分特性,用标准静态相量表示准稳态假设的模型,等同于(5)
步骤2、改进节点法MNA是标准节点分析的扩展,在标准节点分析中,y为n节点网络中n-1个节点的电压,其中的一个节点作为参考节点,用每个节点的基尔霍夫电流定律来构造F,其中支路电流根据节点电压编写;
为了容纳受控电流组件MNA把支路电流作为变量包含在y中,MNA把每个节点使用KCL作为起点F,并尽可能消除分支电流,此外,MNA的节点电压和分支电流不受限制;
步骤3、向后差分公式BDF用于数值求解EDA,BDF方法用于EMT模拟,其中预测变量对正弦变化量进行调整;阶数为k的BDF方法应用在公式(1)时具有如下形式:
其中,β0和h为系数,aj为变量j对应的系数,aj取决于k,阶数k=1的BDF方法等价于后向欧拉法,yn是系统约束方程,使用牛顿迭代方法求yn,(6)是通过微分找到的,F用yn表示;
J是雅可比矩阵,其中
基于动态相量的仿真方法是使用MNA构建一组DAE,
通过将(8)代入(7)得到雅可比矩阵J
其中,和/>为稀疏矩阵,u是常数输入;
载体是在实数集合上定义的系统变量,根据实数值定义转子和控制方程,复杂变量与系统的网络方程相关联,复杂变量和方程包含在y和F,和其他暂态稳定性模拟器相同的方法将它们分成两个实变量和方程;
载体表示系统方程的常数输入,向量/>表示非线性和暂态分量的系统方程;
发电机通过其转子和控制方程提供辅助变量,转子变量分别为磁场和阻尼绕组中的磁通量以及转子的角度和速度,控制变量取决于模拟发电机的励磁和调速器的类型;
系统的复杂网络变量包括所有总线电压和代数总线电压,与分支、负载和定子绕组相关的感应电流也是微分网络变量。
进一步,一个线性串联支路用于表示电力系统模型中的线性元件,连接两个终端总线的中间部分包含R和L,和C=0,和ni=nj=1,充电电容用在每条总线上接地的两个元件进行建模,这些并联元件的C值等于充电电容的一半,R=L=0,和ni=nj=1;
假定R,L,和C都是非零的,电感电流和电容电压是约束方程的辅助变量;
其中,R为电阻,L为电感,C为电容,和/>分别为图2的输入电压,ni和nj表示变压器匝数,/>是变压后的电流,节点关联的KCL方程中i和j电流分别为/>和/>当C=0时变量/>和(11)被省略,表达式/>从(11)获得,当L=0时代入(10)和表达式/>和/>最后,当R=L=0或者R=C=0不需要辅助变量,并把支路添加为A和T的电容或作为对A的电导;
动态相量的同步电机模型是向电机的等效次暂态电压源添加微分分量来整合;由于机器的非线性方程包含在(1)中,从电气角度来看,电机为非线性电感,发生数值故障时这是诊断和分析解决系统的重要依据;
通过定子的微分方程来推导同步电机的特性,ψd和ψq分别是d轴和q轴定子磁通量,id和iq是定子电流,vd和vq是端电压,ωr是机器的转子速度,Ra是电枢电阻,定子的微分方程为:
δr是机器的转子角,它等于转子的q轴和网络的实轴之间的角度,xd和xq分别是参考系上转子的两个任意d轴和q轴量,对应于xd和xq的网络参考系中是动态相量,/>与xd和xq通过旋转
通过乘以反转
其中,是定子电流,/>是定子电压,将(12)和(13)代入(14)得到现有网络框架下的定子微分方程,并以动态相量的方式表示:
其中L1是定子的漏抗,是气隙通量,从(16)中对/>或者/>进行消除,/>是转子变量的非线性组合,选择/>作为定子变量,确保引入的非线性/>约束于定子方程;
将(17)带入(16)式,推导出定子微分方程的最终形式,定子的微分方程为
定子是与有关联的电感分支,用暂态稳定性微分方程对同步电机的转子动力学进行建模,将饱和磁场和包含在机器模型中,使用转子机械摆动方程建模,如下:
其中是定子电流,H是磁场强度,Tm是电机转矩,与/>关联的约束方程是定子的微分方程,如果同步电机的MVA基数与系统的MVA基数之比是η,则电机的终端电流/>包含在终端母线的KCL方程中。
再进一步,高压直流输电系统,
假设转换器有N个桥,其对应的变压器在转换器侧的基极电压为Vbase,Idc是变压侧电流,XT是变压器匝数,是转换器母线电压,θ是/>的角度,φ是转换器锁相环测量的角度,在稳定状态下等于θ,α是转换器的触发角,μ是转换器的重叠角,Vdc是转换器的直流电压;;
换向通过有效直流电阻和电感来对直流系统产生影响,Reff是有效直流电阻,是直流线路电阻的组合,Rdc是换向电阻,模拟了换向引起的直流电压降,有效直流电阻为:
有效直流电感为:
其中Leff是有效直流电感,是直流线路电感的组合,Ldc是变压器电感的平均值,ω0为电压频率;
高压直流输电系统模型包含的微分变量Idc,r,Idc,i和Vm,单位时间内各变量对应的微分方程为:
其中,Idc,r是整流侧电流,Idc,i是逆变侧电流,Vm电容电压,Vdc,r是整流后电压,Vdc,i是逆变侧电压,和/>是端口电流,Cdc是滤波电容,用于整流器和逆变器的KCL方程,Ibase是变压器中转换器侧的峰值电流,转换器的端口电流分为两部分,第一部分是转换器汲取的直流电流,等于
第二部分是换向时端口的交流电压对其电流的影响,
转换器的总端口电流是两部分的总和,
一种基于动态相量的电磁暂态仿真频域扩展装置,所述装置的处理过程采用上述扩展方法。
本发明的技术构思为:所提基于改进节点法(modified nodal analysis,MNA)的电路仿真方法,这种方法用来构建连续微分代数方程(differential-algebraicequations,DAE)。提出的方法利用了stamp技术,很大的提高了MNA的模拟器。推导出常见的电力系统模型。提出了一种基于MNA同步电机的新模型,该模型将电机等效为非线性电感而不是次暂态的等效物。使用可变步长对DAE进行数值求解。
本发明的有益效果主要表现在:
1)使用动态相量来提高电磁暂态仿真,动态相量是稳态相量的自然延伸,它可以动态进行建模,同时保留了调解载波频率的优点;
2)使用动态相量可以提高动态相量的频率带宽进行稳定仿真;
3)该方法基于动态相量和改进节点法,适用于现实的电力系统模型。
附图说明
图1是仿真系统结构图。
图2是一般线性阻抗元件分支图。
图3是高压直流输电系统模型图。
图4是高压直流输电系统流程图。
具体实施方式
下面结合附图对本发明作进一步描述。
参照图1~图4,一种基于动态相量的电磁暂态仿真频域扩展方法,包括以下步骤:
步骤1、电力系统动态微分代数方程DAE如下:
其中,t为时间,y是系统变量,是y的导数;
模拟大型电力系统模型中网络层的瞬变,使用平衡的三相动态相量框架,框架使用Blondel变换来定义动态相量算子给定正弦信号如下:
x(t)=X(t)cos(ω0t+δ(t)) (31)
其中,ω0是频率,δ(t)是相移,cos是余弦函数,X(t)为傅里叶级数;
计算x(t)的动态相量
其中,是动态相量算子,j是虚函数的虚部。正弦对应的/>通过重新引入载体来恢复:
其中,Re是载体函数;
动态相量模型的导数定义为:
用单位时间表示动态相量微分特性,用标准静态相量表示准稳态假设的模型,等同于(5)
步骤2、改进节点法MNA是标准节点分析的扩展,在标准节点分析中,y为n节点网络中n-1个节点的电压,其中的一个节点作为参考节点,用每个节点的基尔霍夫电流定律来构造F,其中支路电流根据节点电压编写;标准节点法的主要缺点是当网络电压不能表示支路电流时,它就会失效。
为了容纳受控电流组件MNA把支路电流作为变量包含在y中,MNA把每个节点使用KCL作为起点F,并尽可能消除分支电流,
此外,MNA的节点电压和分支电流不受限制;
MNA的特点是适合一般电路仿真,可以使用邮票概念通过结构信息来构建电路方程,邮票使用矩阵和抽象索引进行可视化,但它们最重要的特征是把内部变量与节点电压定义一组相关联的约束方程;
步骤3、向后差分公式BDF用于数值求解EDA,BDF系列方法适用于求解复杂方程组,并且已用于很多暂态稳定性仿真研究,BDF方法用于EMT模拟,其中预测变量对正弦变化量进行调整;阶数为k的BDF方法应用在公式(1)时具有如下形式:
其中β0和h为系数,aj为变量j时对应的系数,aj取决于k,阶数k=1的BDF方法等价于后向欧拉法,yn是系统约束方程,使用牛顿迭代等方法求yn,(6)是通过微分找到的,F用yn表示;
J是雅可比矩阵,其中
基于动态相量的仿真方法是使用MNA构建一组DAE,
通过将(8)代入(7)得到雅可比矩阵J
其中,和/>为稀疏矩阵,u是常数输入;
将线性系统方程贡献多少划分为两个稀疏矩阵乘积,这可以并行执行来提高性能;此外,不需要计算分支和负载等严格线性组件的方程F,这个属性很重要,因为计算F对仿真时间有重大影响;
载体是在实数集合上定义的系统变量,根据实数值定义转子和控制方程,复杂变量与系统的网络方程相关联,复杂变量和方程包含在y和F,和其他暂态稳定性模拟器相同的方法将它们分成两个实变量和方程;
载体表示系统方程的常数输入,向量/>表示非线性和暂态分量的系统方程;
发电机通过其转子和控制方程提供辅助变量,转子变量分别为磁场和阻尼绕组中的磁通量以及转子的角度和速度,控制变量取决于模拟发电机的励磁和调速器的类型,在变量选择上MNA方法提供了很大的灵活性,因为可以使用任意的变量来简化方程式;通用模型(如生成器)中的其他变量也会对仿真性能产生负面影响,因为它们会增加一组DAE,因此,用状态方程来模拟发电机及其控制;
系统的复杂网络变量包括所有总线电压,由于与传输线容易受到干扰,这些电压是不同的;大多数实用模型还包括代数总线电压,这是因为电容没有接地连接总线导致的;与分支、负载和定子绕组相关的感应电流是微分网络变量;大多数母线电压和支路电流是独立的,因为暂态稳定性模型由被电感电流分隔的电容性母线组成。
图1说明了一个通用的线性串联支路,可用于表示电力系统模型中的线性元件,连接两个终端总线的中间部分通常包含R和L,和C=0,和ni=nj=1,充电电容用在每条总线上接地的两个元件进行建模,这些并联元件的C值等于充电电容的一半,R=L=0,和ni=nj=1。
图2中所示的元件,假定R,L,和C都是非零的,电感电流和电容电压是约束方程的辅助变量;
其中,R为电阻,L为电感,C为电容,和/>分别为图2的输入电压,ni和nj表示变压器匝数,/>是变压后的电流,节点关联的KCL方程中i和j电流分别为/>和/>当C=0时变量/>和(11)被省略,表达式/>从(11)获得,当L=0时代入(10)和表达式/>和/>最后,当R=L=0或者R=C=0不需要辅助变量,并把支路添加为A和T的电容或作为对A的电导;
早期的动态相量的同步电机模型是向电机的等效次暂态电压源添加微分分量来整合;另一个模型是使用电压源移频分析的方法,该方法也基于次暂态源概念;由于机器的非线性方程包含在(1)中,因此提出了一个通用替代模型,从电气角度来看,电机为非线性电感,发生数值故障时这是诊断和分析解决系统的重要依据;
通过定子的微分方程来推导同步电机的特性,ψd和ψq分别是d轴和q轴定子磁通量,id和iq是定子电流,vd和vq是端电压,ωr是机器的转子速度,Ra是电枢电阻,定子的微分方程为:
δr是机器的转子角,它等于转子的q轴和网络的实轴之间的角度,xd和xq分别是参考系上转子的两个任意d轴和q轴量,对应于xd和xq的网络参考系中是动态相量,/>与xd和xq通过旋转
通过乘以反转
其中,是定子电流,/>是定子电压。将(12)和(13)代入(14)得到现有网络框架下的定子微分方程,并以动态相量的方式表示
其中L1是定子的漏抗,是气隙通量。从(16)中对/>或者/>进行消除,它们通过定子的漏抗L1和气隙通量/>是转子变量的非线性组合,无论如何选择变量,定子的微分方程都是非线性的,选择/>作为定子变量,确保引入的非线性/>约束于定子方程;
将(17)带入(16)式,推导出定子微分方程的最终形式,定子的微分方程为
定子是与有关联的电感分支,用暂态稳定性微分方程对同步电机的转子动力学进行建模,将饱和磁场和包含在机器模型中,使用转子机械摆动方程建模,如下:
其中是定子电流,H是磁场强度,Tm是电机转矩,与/>关联的约束方程是定子的微分方程,如果同步电机的MVA基数与系统的MVA基数之比是η,则电机的终端电流/>包含在终端母线的KCL方程中。
高压直流输电系统,如图3所示,
假设转换器有N个桥,其对应的变压器在转换器侧的基极电压为Vbase,Idc是变压侧电流,XT是变压器匝数,是转换器母线电压,θ是/>的角度,φ是转换器锁相环测量的角度,在稳定状态下等于θ,α是转换器的触发角,μ是转换器的重叠角,Vdc是转换器的直流电压;;
换向通过有效直流电阻和电感来对直流系统产生影响,Reff是有效直流电阻,是直流线路电阻的组合,Rdc是换向电阻,模拟了换向引起的直流电压降,有效直流电阻为:
有效直流电感为:
同样,Leff是有效直流电感,是直流线路电感的组合,Ldc是变压器电感的平均值,ω0为电压频率。
图3所示的高压直流输电系统模型包含的微分变量Idc,r,Idc,i和Vm,单位时间内各变量对应的微分方程为:
其中,Idc,r是整流侧电流,Idc,i是逆变侧电流,Vm电容电压,Vdc,r是整流后电压,Vdc,i是逆变电压,和/>是端口电流,Cdc是滤波电容,用于整流器和逆变器的KCL方程,Ibase是变压器中转换器侧的峰值电流,转换器的端口电流分为两部分,第一部分是转换器汲取的直流电流,等于
第二部分是换向时端口的交流电压对其电流的影响,
转换器的总端口电流是两部分的总和,
一种基于动态相量的电磁暂态仿真频域扩展装置,所述装置的处理过程采用上述扩展方法。
本说明书的实施例所述的内容仅仅是对发明构思的实现形式的列举,仅作说明用途。本发明的保护范围不应当被视为仅限于本实施例所陈述的具体形式,本发明的保护范围也及于本领域的普通技术人员根据本发明构思所能想到的等同技术手段。
Claims (4)
1.一种基于动态相量的电磁暂态仿真频域扩展方法,其特征在于,所述方法包括以下步骤:
步骤1、电力系统动态微分代数方程DAE如下:
其中,t为时间,y是系统变量,是y的导数;
模拟大型电力系统模型中网络层的瞬变,使用平衡的三相动态相量框架,框架使用Blondel变换来定义动态相量算子给定正弦信号如下:
x(t)=X(t)cos(ω0t+δ(t)) (2)
其中,ω0是频率,δ(t)是相移,cos是余弦函数,X(t)为傅里叶级数;
计算x(t)的动态相量
其中,是动态相量算子,j是虚函数的虚部,正弦对应的/>通过重新引入载体来恢复:
其中,Re是载体函数;
动态相量模型的导数定义为:
用单位时间表示动态相量微分特性,用标准静态相量表示准稳态假设的模型,等同于(5)中
步骤2、改进节点法MNA是标准节点分析的扩展,在标准节点分析中,y为n节点网络中n-1个节点的电压,其中的一个节点作为参考节点,用每个节点的基尔霍夫电流定律来构造F,其中支路电流根据节点电压编写;
为了容纳受控电流组件MNA把支路电流作为变量包含在y中,MNA把每个节点使用KCL作为起点F,并尽可能消除分支电流,此外,MNA的节点电压和分支电流不受限制;
步骤3、向后差分公式BDF用于数值求解EDA,BDF方法用于EMT模拟,其中预测变量对正弦变化量进行调整;阶数为k的BDF方法应用在公式(1)时具有如下形式:
其中,β0和h为系数,aj为变量j时对应的系数,aj取决于k,阶数k=1的BDF方法等价于后向欧拉法,yn是系统约束方程,使用牛顿迭代等方法求yn,(6)是通过微分找到的,F用yn表示;
J是雅可比矩阵,其中
基于动态相量的仿真方法是使用MNA构建一组DAE,
通过将(8)代入(7)得到雅可比矩阵J
其中,和/>为稀疏矩阵,u是常数输入;
载体是在实数集合上定义的系统变量,根据实数值定义转子和控制方程,复杂变量与系统的网络方程相关联,复杂变量和方程包含在y和F,和其他暂态稳定性模拟器相同的方法将它们分成两个实变量和方程;
载体表示系统方程的常数输入,/>表示非线性和暂态分量的系统方程;
发电机通过其转子和控制方程提供辅助变量,转子变量分别为磁场和阻尼绕组中的磁通量以及转子的角度和速度,控制变量取决于模拟发电机的励磁和调速器的类型;
系统的复杂网络变量包括所有总线电压和代数总线电压,与分支、负载和定子绕组相关的感应电流是微分网络变量。
2.如权利要求1所述的一种基于动态相量的电磁暂态仿真频域扩展方法,其特征在于,一个线性串联支路用于表示电力系统模型中的线性元件,连接两个终端总线的中间部分包含R和L,和C=0,和ni=nj=1,充电电容用在每条总线上接地的两个元件进行建模,这些并联元件的C值等于充电电容的一半,R=L=0,和ni=nj=1;
假定R,L,和C都是非零的,电感电流和电容电压是约束方程的辅助变量;
其中,R为电阻,L为电感,C为电容,和/>分别为图2的输入电压,ni和nj表示变压器匝数,/>是变压后的电流,节点关联的KCL方程中i和j电流分别为/>和/>当C=0时变量/>和(11)被省略,表达式/>从(11)获得,当L=0时代入(10)和表达式/>和/>最后,当R=L=0或者R=C=0不需要辅助变量,并把支路添加为A和T的电容或作为对A的电导;
动态相量的同步电机模型是向电机的等效次暂态电压源添加微分分量来整合;由于机器的非线性方程包含在(1)中,从电气角度来看,电机为非线性电感,发生数值故障时这是诊断和分析解决系统的重要依据;
通过定子的微分方程来推导同步电机的特性,ψd和ψq分别是d轴和q轴定子磁通量,id和iq是定子电流,vd和vq是端电压,ωr是机器的转子速度,Ra是电枢电阻,定子的微分方程为:
δr是机器的转子角,它等于转子的q轴和网络的实轴之间的角度,xd和xq分别是参考系上转子的两个任意d轴和q轴量,对应于xd和xq的网络参考系中是动态相量,/>与xd和xq通过旋转
通过乘以反转
其中,是定子电流,/>是定子电压,将(12)和(13)代入(14)得到现有网络框架下的定子微分方程,并以动态相量的方式表示:
其中L1是定子的漏抗,是气隙通量,从(16)中对/>或者/>进行消除,/>是转子变量的非线性组合,选择I作为定子变量,确保引入的非线性/>约束于定子方程;
将(17)带入(16)式,推导出定子微分方程的最终形式,定子的微分方程为
定子是与有关联的电感分支,用暂态稳定性微分方程对同步电机的转子动力学进行建模,将饱和磁场和包含在机器模型中,使用转子机械摆动方程建模,如下:
其中是定子电流,H是磁场强度,Tm是电机转矩,与/>关联的约束方程是定子的微分方程,如果同步电机的MVA基数与系统的MVA基数之比是η,则电机的终端电流/>包含在终端母线的KCL方程中。
3.如权利要求1或2所述的一种基于动态相量的电磁暂态仿真频域扩展方法,其特征在于,高压直流输电系统,
假设转换器有N个桥,其对应的变压器在转换器侧的基极电压为Vbase,Idc是变压侧电流,XT是变压器匝数,是转换器母线电压,θ是V的角度,φ是转换器锁相环测量的角度,在稳定状态下等于θ,α是转换器的触发角,μ是转换器的重叠角,其中Vdc是转换器的直流电压;
换向通过有效直流电阻和电感来对直流系统产生影响,Reff是有效直流电阻,是直流线路电阻的组合,Rdc是换向电阻,模拟了换向引起的直流电压降,有效直流电阻为:
有效直流电感为:
其中,Leff是有效直流电感,是直流线路电感的组合,Ldc是变压器电感的平均值,ω0为电压频率;
高压直流输电系统模型包含的微分变量Idc,r,Idc,i和Vm,单位时间内各变量对应的微分方程为:
其中,Idc,r是整流侧电流,Idc,i是逆变侧电流,Vm电容电压,Vdc,r是整流后电压,Vdc,i是逆变电压,和/>是端口电流,Cdc是滤波电容,用于整流器和逆变器的KCL方程,Ibase是变压器中转换器侧的峰值电流,转换器的端口电流分为两部分,第一部分是转换器汲取的直流电流,等于
第二部分是换向时端口的交流电压对其电流的影响,
转换器的总端口电流是两部分的总和,
4.一种如权利要求1所述的基于动态相量的电磁暂态仿真频域扩展方法实现的装置,所述装置的处理过程采用上述扩展方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310623078.8A CN116722563A (zh) | 2023-05-30 | 2023-05-30 | 一种基于动态相量的电磁暂态仿真频域扩展方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310623078.8A CN116722563A (zh) | 2023-05-30 | 2023-05-30 | 一种基于动态相量的电磁暂态仿真频域扩展方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116722563A true CN116722563A (zh) | 2023-09-08 |
Family
ID=87864031
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310623078.8A Pending CN116722563A (zh) | 2023-05-30 | 2023-05-30 | 一种基于动态相量的电磁暂态仿真频域扩展方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116722563A (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104732033A (zh) * | 2015-03-31 | 2015-06-24 | 河南行知专利服务有限公司 | 一种电力系统机电暂态过程的仿真分析方法 |
CN104866665A (zh) * | 2015-05-19 | 2015-08-26 | 清华大学 | 基于接口等值与交互的含电力电子设备的混合仿真方法 |
CN110569558A (zh) * | 2019-08-16 | 2019-12-13 | 上海交通大学 | 适用于微电网实时仿真的混合电磁暂态仿真方法 |
WO2020088245A1 (zh) * | 2018-10-31 | 2020-05-07 | 邓宏伟 | 一种基于直算法的电磁机电暂态仿真算法 |
CN112214905A (zh) * | 2020-10-22 | 2021-01-12 | 华北电力大学(保定) | 一种电力系统宽频带建模分析与仿真方法及系统 |
US20210050722A1 (en) * | 2019-08-16 | 2021-02-18 | Shanghai Jiao Tong University | Hybrid electromagnetic transient simulation method for microgrid real-time simulation |
CN114943133A (zh) * | 2022-03-18 | 2022-08-26 | 国网浙江省电力有限公司嘉兴供电公司 | 配电网细粒度并行电磁暂态仿真方法、系统、终端及介质 |
-
2023
- 2023-05-30 CN CN202310623078.8A patent/CN116722563A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104732033A (zh) * | 2015-03-31 | 2015-06-24 | 河南行知专利服务有限公司 | 一种电力系统机电暂态过程的仿真分析方法 |
CN104866665A (zh) * | 2015-05-19 | 2015-08-26 | 清华大学 | 基于接口等值与交互的含电力电子设备的混合仿真方法 |
WO2020088245A1 (zh) * | 2018-10-31 | 2020-05-07 | 邓宏伟 | 一种基于直算法的电磁机电暂态仿真算法 |
CN110569558A (zh) * | 2019-08-16 | 2019-12-13 | 上海交通大学 | 适用于微电网实时仿真的混合电磁暂态仿真方法 |
US20210050722A1 (en) * | 2019-08-16 | 2021-02-18 | Shanghai Jiao Tong University | Hybrid electromagnetic transient simulation method for microgrid real-time simulation |
CN112214905A (zh) * | 2020-10-22 | 2021-01-12 | 华北电力大学(保定) | 一种电力系统宽频带建模分析与仿真方法及系统 |
CN114943133A (zh) * | 2022-03-18 | 2022-08-26 | 国网浙江省电力有限公司嘉兴供电公司 | 配电网细粒度并行电磁暂态仿真方法、系统、终端及介质 |
Non-Patent Citations (1)
Title |
---|
M. A. KULASZA, ET AL: "Extending the Frequency Bandwidth of Transient Stability Simulation Using Dynamic Phasors", IEEE TRANSACTIONS ON POWER SYSTEMS, vol. 37, no. 1, pages 249 - 259, XP011894775, DOI: 10.1109/TPWRS.2021.3094451 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Stankovic et al. | Analysis of asymmetrical faults in power systems using dynamic phasors | |
Wang et al. | Methods of interfacing rotating machine models in transient simulation programs | |
Lin et al. | A wide-band multi-port system equivalent for real-time digital power system simulators | |
Fan et al. | Time domain transformation method for accelerating EMTP simulation of power system dynamics | |
Chapariha et al. | Constant-Parameter $ RL $-Branch Equivalent Circuit for Interfacing AC Machine Models in State-Variable-Based Simulation Packages | |
Zhang et al. | Synchronous machine modeling based on shifted frequency analysis | |
Dehkordi et al. | Development and validation of a comprehensive synchronous machine model for a real-time environment | |
Kunjumuhammed et al. | Simulation of power system with renewables | |
CN108573094B (zh) | 同步发电机的vbr电磁暂态仿真模型的建立方法及系统 | |
Xia et al. | Multi-scale induction machine model in the phase domain with constant inner impedance | |
Plumier | Co-simulation of electromagnetic transients and phasor models of electric power systems | |
Berhausen et al. | Determination of high power synchronous generator subtransient reactances based on the waveforms for a steady state two-phase short-circuit | |
Noon et al. | A power hardware-in-the-loop testbench for aerospace applications | |
Kulasza et al. | Extending the frequency bandwidth of transient stability simulation using dynamic phasors | |
CN103762916B (zh) | 一种双y移30°六相同步发电机vbr模拟方法 | |
CN116722563A (zh) | 一种基于动态相量的电磁暂态仿真频域扩展方法及装置 | |
Huang et al. | Analysis of balanced and unbalanced faults in power systems using dynamic phasors | |
Almozayen et al. | Dynamic phasor finite element modeling of grid-connected DFIG considering winding space harmonics | |
Quadri et al. | Modeling of classical synchronous generators using size-efficient lookup tables with skewing effect | |
Su et al. | A new approach for integration of two distinct types of numerical simulator | |
Ye et al. | Multi-scale modeling and simulation of synchronous machine in phase-domain | |
Demiray et al. | Comparison of the efficiency of dynamic phasor models derived from ABC and DQO reference frame in power system dynamic simulations | |
Almozayen et al. | Modeling the Impact of System Disturbances on Grid-Connected DFIG Using Dynamic Phasor FEM | |
Jain | Dynamic simulation of power systems using three phase integrated transmission and distribution system models: Case study comparisons with traditional analysis methods | |
Dinh et al. | A salient-pole generator model for harmonic analysis |
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 |