CN101719182B - 一种交直流电力系统分割并行电磁暂态数字仿真方法 - Google Patents

一种交直流电力系统分割并行电磁暂态数字仿真方法 Download PDF

Info

Publication number
CN101719182B
CN101719182B CN2009102418641A CN200910241864A CN101719182B CN 101719182 B CN101719182 B CN 101719182B CN 2009102418641 A CN2009102418641 A CN 2009102418641A CN 200910241864 A CN200910241864 A CN 200910241864A CN 101719182 B CN101719182 B CN 101719182B
Authority
CN
China
Prior art keywords
subnet
network
alpha
direct current
alternating current
Prior art date
Application number
CN2009102418641A
Other languages
English (en)
Other versions
CN101719182A (zh
Inventor
田芳
周孝信
Original Assignee
中国电力科学研究院
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 中国电力科学研究院 filed Critical 中国电力科学研究院
Priority to CN2009102418641A priority Critical patent/CN101719182B/zh
Publication of CN101719182A publication Critical patent/CN101719182A/zh
Application granted granted Critical
Publication of CN101719182B publication Critical patent/CN101719182B/zh

Links

Abstract

本发明涉及一种交直流电力系统分割并行的电磁暂态数字仿真方法,包括:将目标交直流电力系统分割为多个子网,包括交流网络子网、换流器子网、直流网络子网或交直流混合网络子网;对各子网求取端口点的等值阻抗矩阵和等值电势;求取子网间联络线电流;根据子网间联络线电流求取各子网节点电压;求解各子网的发电机、控制系统等微分方程;进行用户自定义模型、MATLAB模型、物理装置输入输出信息处理;进行机电暂态接口输入输出信息处理。该仿真方法具有分网方式灵活、计算速度快等优点。该仿真方法能够同时适用于交直流系统电磁暂态分网并行计算、与机电暂态并行计算、接入用户自定义模型、MATLAB模型和物理装置处理等,从而实现对交直流电力系统的各种暂态和动态过程的数字仿真。

Description

一种交直流电力系统分割并行电磁暂态数字仿真方法
技术领域
[0001] 本发明涉及一种交直流电力系统电磁暂态数字仿真方法,更具体地说,涉及一种交直流电力系统的分割并行电磁暂态实时动态仿真方法。
背景技术
[0002] 电力系统电磁暂态过程仿真的主要目的在于分析和计算故障或操作后可能出现的暂态过电压和过电流,以便根据所得到的暂态过电压和过电流对相关电力设备进行合理设计,确定已有设备能否安全运行,并研究相应的限制和保护措施。
[0003] 电磁暂态数字仿真一般采用三相瞬时值模型,计算步长取20-200微秒。
[0004] 在直流输电系统的电磁暂态仿真中,换流器用三相暂态模型模拟;直流线路用π 型、T型集中参数线路模型或分布参数线路模型;直流控制系统的模型包括调节系统和脉冲触发系统等。直流输电系统的电磁暂态仿真,不仅可用于交直流系统的一般稳定性分析, 还可以用来深入地分析直流输电系统的内部问题,如换相失败问题、直流控制和保护等问题。另外,如果仿真的速度足够快,能达到实时,还可以用于直流控制与保护装置的试验。但是,相对于交流系统而言,要实现含直流输电系统的电力系统电磁暂态仿真的实时性并不容易,这是因为,直流输电系统中的关键元件-换流器是由若干个换流阀组成的,在一个周期内,换流阀会多次导通或关断,每一次导通或关断都意味着网络的拓朴结构发生变化。现有的电磁暂态数字仿真方法,将交直流电力系统中的交流网络、换流器、直流网络等各部分分别形成暂态等值计算电路,然后根据其连接关系形成统一的网络节点电压方程来进行求解。换流阀导通或关断导致网络拓扑结构发生变化时,该节点电压方程的计算电导阵需要重新进行三角(LU)分解,由于网络拓扑变化频繁,会引起计算量的大幅度上升,难以实现实时仿真。
发明内容
[0005] 本发明的目的就是提供一种交直流电力系统分割并行电磁暂态数字仿真方法,有效地解决上述交直流电力系统电磁暂态仿真的实时性问题。该方法的核心是交流网络与直流各部分分开计算,换流阀导通或关断后只影响其所在的换流器网络结构,需对其所在的换流器网络重新进行三角LU分解,不影响其它网络结构,因此可避免大交流网络和其它部分(如直流线路、其它换流器)的频繁LU分解过程。同时采用网络并行计算技术,以提高计算速度。
[0006] 为此,根据本发明的一个方面,提供一种交直流电力系统分割并行电磁暂态数字仿真方法,包括如下步骤:
[0007] 步骤1 :将一个目标交直流电力系统分割为多个子网,包括交流网络子网、换流器子网和直流网络子网;
[0008] 步骤2 :对各子网求取端口点的等值阻抗矩阵;
[0009] 步骤3 :设定初始计算时间t = 0 ;[0010] 步骤4:判断此时刻是否有开关动作,或换流阀导通、关断等使网络结构发生变化的情况发生;如果是则转至步骤5,否则转至步骤6 ;
[0011] 步骤5 :对有开关动作或者换流阀导通或关断等使网络结构发生变化的情况发生的子网,求取端口点的等值阻抗矩阵;
[0012] 步骤6 :各子网求取端口点的等值电势;
[0013] 步骤7 :求取子网间的联络线电流或联络线变量(在换流器采用节点法情况下为联络线电流,在换流器采用回路法求解情况下为联络线变量,联络线变量可以是子网间联络线电流,也可以是子网间端口点电压);
[0014] 步骤8 :各子网根据得到的子网间联络线电流或联络线变量求取子网节点电压, 在换流器采用回路法求解情况下根据子网间联络线变量求取换流器连支电流;
[0015] 步骤9 :各子网求解发电机、控制系统等微分方程;
[0016] 步骤10 :判断是否有用户自定义模型或MATLAB模型;如果是则转至步骤11和12, 否则转至步骤13 ;
[0017] 步骤11 :用户自定义模型、MATLAB模型输入输出信息处理;
[0018] 步骤12 :用户自定义模型、MATLAB模型计算;
[0019] 步骤13 :判断是否有机电暂态接口 ;如果是则转至步骤14和15,否则转至步骤 16 ;
[0020] 步骤14 :机电暂态接口输入输出信息处理;
[0021] 步骤15 :机电暂态仿真计算;
[0022] 步骤16 :判断是否有物理装置;如果是则转至步骤17和18,否则转至步骤19 ;
[0023] 步骤17 :物理装置输入输出信息处理;
[0024] 步骤18 :物理装置运行;
[0025] 步骤19 :将仿真时间增加一个时间步长;
[0026] 重复执行上述步骤4至步骤19,直至到达总仿真时间为止。
[0027] 根据本发明的另一个方面,提供一种电力系统数字仿真方法,包括如下步骤:
[0028] 步骤1 :将一个目标交直流电力系统分割为多个子网,子网可以是交直流混合网络,也可以是单独的交流、换流器和直流网络;
[0029] 步骤2 :判断该子网是否为交直流混合网络;如果是则转至步骤3,否则转至步骤 5 ;
[0030] 步骤3 :对于交直流混合网络,求取子网内部交直流联络线电流;
[0031] 步骤4:对于交直流混合网络,根据子网内部交直流联络线电流,分别求取交流网络、换流器网络和直流网络节点电压;
[0032] 步骤3和4的目的是计算各子网端口点的等值阻抗矩阵。
[0033] 步骤5 :对于非交直流混合网络,直接求取端口点的等值阻抗矩阵;
[0034] 步骤6 :设定初始计算时间t = 0 ;
[0035] 步骤7 :判断此时刻是否有开关动作,或换流阀导通、关断等使网络结构发生变化的情况发生;如果是则转至步骤8,否则转至步骤12 ;
[0036] 步骤8 :对有开关动作或者换流阀导通或关断等使网络结构发生变化的情况发生的子网,判断该子网是否为交直流混合网络;如果是则转至步骤9,否则转至步骤11 ;[0037] 步骤9 :对于有开关动作或者换流阀导通或关断等使网络结构发生变化的情况发生的交直流混合网络,求取子网内部交直流联络线电流;
[0038] 步骤10 :对于有开关动作或者换流阀导通或关断等使网络结构发生变化的情况发生的交直流混合网络,根据子网内部交直流联络线电流,分别求取交流网络、换流器网络和直流网络节点电压;
[0039] 步骤9和10的目的是计算网络结构变化子网端口点的等值阻抗矩阵。
[0040] 步骤11 :对有开关动作或者换流阀导通或关断等使网络结构发生变化的情况发生的非交直流混合网络,直接求取端口点的等值阻抗矩阵;
[0041] 步骤12 :判断该子网是否为交直流混合网络;如果是则转至步骤13,否则转至步骤15 :
[0042] 步骤13 :对于交直流混合网络,求取子网内部交直流联络线电流;
[0043] 步骤14 :对于交直流混合网络,根据子网内部交直流联络线电流,分别求取交流网络、换流器网络和直流网络节点电压;
[0044] 步骤13和14的目的是计算各子网端口点的等值电势。
[0045] 步骤15 :对于非交直流混合网络,直接求取端口点的等值电势;
[0046] 步骤16 :求取子网间联络线电流;
[0047] 步骤17 :各子网根据子网间联络线电流求取子网节点电压;
[0048] 步骤18 :各子网求解发电机、控制系统等微分方程;
[0049] 步骤19 :判断是否有用户自定义模型或MATLAB模型;如果是则转至步骤20和21, 否则转至步骤22 ;
[0050] 步骤20 :用户自定义模型、MATLAB模型输入输出信息处理;
[0051] 步骤21 :用户自定义模型、MATLAB模型计算;
[0052] 步骤22 :判断是否有机电暂态接口 ;如果是则转至步骤23和24,否则转至步骤 25 ;
[0053] 步骤23 :机电暂态接口输入输出信息处理;
[0054] 步骤M :机电暂态仿真;
[0055] 步骤25 :判断是否有物理装置;如果是则转至步骤沈和27,否则转至步骤28 ;
[0056] 步骤沈:物理装置输入输出信息处理;
[0057] 步骤27 :物理装置运行;
[0058] 步骤28 :将仿真时间增加一个时间步长;
[0059] 重复执行上述步骤7至步骤观,直至到达总仿真时间为止。
[0060] 本发明的有益效果是:根据本发明的交直流电力系统分割并行电磁暂态数字仿真方法,通过将目标电力系统分割为交流、换流器和直流多个子网,进而将涉及各个子网的计算任务分配到机群服务器的各个子机并行计算,从而实现交直流电力系统的实时或超实时电磁暂态仿真。其网络分割方式非常灵活,可以从交直流的任一连接点进行分网,交流网络还可进一步分网,分割成若干个交流子网。
[0061] 该方法与机电暂态并行仿真方法结合,能够实现大规模电力系统的机电暂态和电磁暂态过程的并行仿真计算,实现大电力系统背景下对局部电网如直流输电、电力电子装置等的详细电磁暂态仿真模拟。[0062] 此外,可以在该方法所模拟的电力系统中接入用户自定义模型或MATLAB模型所模拟的控制系统,用于控制系统的结构设计和参数优化,这样可极大地扩展对电力系统的仿真分析能力,同时可以在该方法所模拟的电力系统中接入实际物理装置(如继电保护和安全自动装置、直流输电控制装置等)进行闭环仿真试验,以检验这些物理装置在实际电力系统中的效果。
附图说明
[0063] 图1是根据本发明的第一实施例的交直流电力系统分割并行电磁暂态数字仿真方法的流程图;
[0064] 图2示意性示出了一种交直流电力系统分割方式;
[0065] 图3示出了换流器采用回路法求解情况下换流器电压、电流的正方向;
[0066] 图4是根据本发明的第二实施例的交直流电力系统分割并行电磁暂态数字仿真方法的流程图;
[0067] 图5示意性示出了另一种交直流电力系统分割方式;
[0068] 图6是本发明的第二实施例中针对特殊情况的简化流程图;
具体实施方式
[0069] 图1是根据本发明第一实施例的交直流电力系统分割并行电磁暂态数字仿真方法的流程图。
[0070] 如图1所示,根据本发明第一实施例的交直流电力系统分割并行电磁暂态数字仿真方法包括如下步骤:
[0071] 步骤1 :交直流电力系统分割。
[0072] 在该步骤1中,通过节点分裂或分布参数线路解耦的方法将目标交直流电力系统分为多个子网,包括交流网络子网、换流器子网和直流网络子网。交流网络子网或换流器子网或直流网络子网可以是1个,也可以是多个。
[0073] 图2示出了一种交直流电力系统分割方式。如图2所示,目标电力系统被分割为六个子网Sl〜S6,其中S1、S2为换流器子网,S3为直流网络子网,S4、S5和S6为交流网络子网。
[0074] 步骤2 :各子网求取端口点的等值阻抗矩阵。
[0075] 例如,交流子网S4的节点电压方程如下:
[0076] GAVA+PABi α -pAci Y +Pacii α !+Pacl i α 2 = hA (1)
[0077] 其中^乂、1^分别为交流网络六的导纳阵、节点电压及注入电流源,11、12为从交流网络流向换流器1、2的电流,pac;1为反映某一交流节点与ial(ia2)关联关系的关联阵, 其中元素非0即1,i a、i γ为交流网络A至B、C至A的电流,Pab,Pac为反映交流网络A中某一交流节点与ia、关联关系的关联阵,其中元素非0即1。
[0078] 令式(1)中、=0,ia、iY、ial、ia2依次分别置1或-1(电流正方向为流入该子网则置1,否则置-1),求取端口点节点电压,求出的电压向量列的组合即为端口点等值阻抗矩阵。
[0079] 步骤3 :设定初始计算时间t = 0。[0080] 步骤4 :判断此时刻是否有开关动作,或换流阀导通、关断等使网络结构发生变化的情况发生。
[0081] 在该步骤4中,判断此时刻各个子网中是否有开关动作或者换流阀导通或关断等情况发生。如果判断为“是”,则前进到步骤5 ;如果判断为“否”,则前进到步骤6。
[0082] 步骤5 :对有开关动作或者换流阀导通或关断等使网络结构发生变化的情况发生的子网,求取端口点的等值阻抗矩阵。
[0083] 计算过程同步骤2。 步骤6 :各子网求取端口点的等值电势。[0084]步骤6:各子网求取端口点的等值电势。
[0085] 令式(1)中ia、iY、ial、ia2 = 0,求取端ロ点节点电压,即为端ロ点等值电势。
[0086] 步骤7 :求取子网间联络线电流。
[0087] 根据各子网端ロ点的等值阻抗矩阵和等值电势求取子网间联络线电流。
[0088] 设已求得子网S4的端ロ点的等值阻抗和等值电势,可写出其端ロ方程为
[0089]
Figure CN101719182BD00101
[0090] 式中,Val,va2为子网Sl与子网S4之间的端ロ点1(与换流器1相连)、端ロ点 2(与换流器2相连)电压,Va为子网S4与子网S5之间的端ロ点电压,v¥为子网S4与子
网S6之间的端ロ点电压;eal(s4),ea2(s4), ea(s4), e¥(s4)为各端ロ点等值电势;阻抗矩阵对角 元素zal,al(S4)、za2,a2(S4)等为端ロ点等值自阻杭,非对角元素zal,a2(S4)、zal,a(S4)等为端ロ 点等值互阻杭。
[0091 ] 类似地,可写出其他子网的端ロ方程分別为
[0092] 子网 Sl :
[0093]
Figure CN101719182BD00102
[0094] 子网 S2 :
[0095]
Figure CN101719182BD00103
[0096] 子网 S3 :
[0097]
Figure CN101719182BD00104
[0098] 子网 S5 :
Figure CN101719182BD00111
[0100]子网 S6
Figure CN101719182BD00112
[0102]将式⑵〜(7)合井,消去 να1,να2,να3,να4,ν。,ve,νγ,ν 0 12' v0 34, 東后得到子
网间联络线电流方程式如下:
[0103] Ai = b (8)
[0104] 其中i =[し,ら,“,し”しジしジしもバ印バ舰ド为子网间联络线电流,矩阵八
如下式所示:
[0105]
Figure CN101719182BD00121
[0106] 右端项b如下式所示:
[0107]
Figure CN101719182BD00131
[0108] 求解线性方程组(8),可求得子网间联络线电流。
[0109] 换流器也可以采用回路法求解,这样可以减少联络线方程式的维数。
[0110] 此时,交流子网节点电压方程不变,为
Figure CN101719182BD00132
[0114] 直流子网节点电压方程不变,为
[0115]
Figure CN101719182BD00133
[0116]子网 S3、S4、S5、S6 的端口方程不变,仍为式(5)、(2). (6)、(7)。
[0117] 子网S1和S2的端口方程改用端口点的等值导纳和等值电流源描述。
[0118] 对于子网S1,设已求得其端口点的等值导纳和等值电流源,可写出其端口方程为 (换流器1、2分开描述)
Figure CN101719182BD00134
[0121] 式中,Vdi,为换流器1、2直流侧电压,idi,为换流器1、2直流侧电流;1。产4), Ia2(S4),Idi(S4),I</4)为各端口点等值电流源;导纳矩阵对角元素yal,ai(Sl)、ydl,di(Sl)等为端口 点等值自导纳,非对角元素yai,di(si)、ydi,ai(si)等为端口点等值互导纳。
[0122] 类似地,可写出子网S2的端口方程为;
Figure CN101719182BD00135
Figure CN101719182BD00141
Figure CN101719182BD00142
Figure CN101719182BD00143
Figure CN101719182BD00144
Figure CN101719182BD00145
Figure CN101719182BD00146
[0124]
0131]
0132]
0133]
0145]
0146]
0147]
0148]
laA ^4
(S2)
'«4
Id
(S2)
4
+
. (S2) (S2)
y y a^,d^
(S2) (S2)
Vd4
(14)
0125] 以子网Sl为例,端口点等值导纳和等值电流源的计算过程如下
0126] 换流器1、2采用回路法求解时,其连支电流方程式为
0127] AliIli = bL1+CL1va !+DuVdl (15)
0128] Al2Il2 = bL2+CL2va2+DL2vd2 (16)
0129] 式中Iu、Ili为换流器1、2的连支电流。
0130] 另外,根据连接关系,有下式成立:
ia\ = CliIzj
L2=C
L2*-L2
(17)
(18)
i =DT
J-dl uLl 1Ll
0134] id2 = Dl21Il2
(19)
(20)
0135] 将式(15)代入式(17),可得:
0136] ial=CTLlAu-l(bn+CuVal+OLlvdl)
0137] 与式(14)对比,可得:
0138] Ijsi) = CLAjbil,= CTLIAL;1CLI
0139] 类似地,根据式(14)、(15)、(19),可得:
ν (51) - Cr A 1D
0140] Id
(Si)
:DIl AIl_1bil ‘ yd\,a\
(S1) - Dr A '1C
^liaII ^LI, yd\A\
^=DU-1D
Llr^Ll
,IA
0141]根据式(14)、(16)、(18),可得
0142] Ia
(si)
2
=Cl2Al2 bi2
(si) _ cr A -ic
(Si)
0143]
0144] Id2
根据式(14)、(16)、00),可得
<si) -ΊΛΤ A 1. ,, (Si)
:Di2Ai2 bi2, _yd2,a2
-lV A " Γ
yal,dl = ^-Ll^Ll
v (S1) - Dr A 1D
另外,根据连接关系,有下式成立
e 12
e 34
=ν,,+ν
dl v d2
=_ (VdS+Vd4)
(21) (22) (23)
[0149] 1 e 34 — ~~ id3 = _id4 (24)[0150] i = D 丄 d3 uLS IL3 (25)[0151] i = D 丄 d4 UL4 Il4 (26)[0152] Val = Va 2 (27)[0153] Va3 = Va 4 (28)[0154] 由式(2) 、(5) 〜(7)、(13)Va » V0, 0 34' 1 β 12' 1 β 34' 1Cil
i-V、uy;〜V…,伯太 Ia π Ia 2,Ia 3,Ia 4,Va 2,Va 4,
_____ 2,id3,id4,最后得到联络线方程式如下:
0155] AX = b (29)
0156] 其中X = [ia,ie,iY,Val,va3, vdl, vd2, vd3, vd4]τ为子网间联络线变量,部分变量为子网间联络线电流,部分变量为子网间端口点电压,矩阵A如下式所示:
[0157]
14
Figure CN101719182BD00151
[0158] 右端项b如下式所示:
[0159]
Figure CN101719182BD00161
[0160] 求解线性方程组(¾),可求得子网间联络线变量X。
[0161] 对于图2所示的系统,换流器采用节点法求解时联络线方程式的维数为23,采用回路法求解时联络线方程式的维数为19。
[0162] 步骤8 :各子网根据得到的子网间联络线电流或联络线变量求取子网节点电压。 例如,交流子网S4的节点电压可根据式(1)求出。在换流器采用回路法求解情况下则根据子网间联络线变量求取换流器连支电流。
[0163] 步骤9 :各子网求解发电机、控制系统等微分方程。
[0164] 在步骤9中,各个子网求解发电机、控制系统等的微分方程。在电力系统的电磁暂态计算中,各种动态元件,包括发电机、电动机、励磁调节器、调速器、电力系统稳定器、直流输电控制器、FACTS控制器等,均用微分方程来描述。这些微分方程可采用梯形隐积分方法将其差分化后,形成差分方程再求解。
[0165] 步骤10 :判断是否有用户自定义模型或MATLAB模型。
[0166] 在该步骤10中,判断各个子网中是否有用户自定义模型或MATLAB模型。如果判断为“是”,则前进到步骤11和12 ;如果判断为“否”,则前进步骤13。
[0167] 步骤11 :用户自定义模型、MATLAB模型输入输出信息处理。
[0168] 在该步骤11中,需要求取并发送各外接模型(用户自定义模型、MATLAB模型)输入变量值,并接收和处理各外接模型的输出变量值。
[0169] 用户自定义模型为用户自行搭建的模型,一般为励磁调节器、调速器、电力系统稳定器、直流输电控制器、FACTS控制器等控制系统模型。MATLAB模型是用户在MATLAB/ Simulink环境下自行搭建的模型,一般也是上述控制系统模型。
[0170] 步骤12 :用户自定义模型、MATLAB模型计算。
[0171] 在步骤12中,需要先接收用户自定义模型、MATLAB模型的输入变量值,然后根据用户自定义模型、MATLAB模型的各功能框计算公式和连接关系,进行模型计算,求取并发送用户自定义模型、MATLAB模型输出变量值。[0172] 步骤12中的计算与其它步骤中的计算并行进行,由另外的独立程序(用户自定义模型计算程序、MATLAB模型计算程序)来实现。
[0173] 如果步骤12中模型的计算量不大,也可以将其与其它步骤串行进行计算,由同一个程序来处理,这时,发送和接收数据可以去掉,步骤11和步骤12合并为图1所示的步骤 12。
[0174] 步骤13 :判断是否有机电暂态接口。
[0175] 在该步骤13中,判断是否有机电暂态接口。如果判断为“是”,则前进到步骤14和 15 ;如果判断为“否”,则前进到步骤16。
[0176] 步骤14 :机电暂态接口输入输出信息处理。
[0177] 在机电暂态仿真的整数倍步长点,求取并向机电暂态仿真程序发送与机电暂态子网接口的边界点正、负、零序电压和电流,接收机电暂态仿真程序发送的边界点正、负、零序等值阻抗和电势。
[0178] 步骤15 :机电暂态仿真计算。
[0179] 步骤15由另外的独立的机电暂态仿真程序来实现,先接收电磁暂态程序发送的边界点正、负、零序电压和电流,进行机电暂态一个步长的积分计算,然后求取并向电磁暂态仿真程序发送与电磁暂态子网接口的边界点正、负、零序等值阻抗和电势。
[0180] 步骤16 :判断是否有物理装置。
[0181] 在该步骤16中,判断是否有物理装置。如果判断为“是”,则前进到步骤17和18 ; 如果判断为“否”,则前进到步骤19。
[0182] 步骤17 :物理装置输入输出信息处理。
[0183] 该步骤主要是进行数模和模数转换,将实时仿真的数字信号和开关信号分别通过 PCI总线的D/A转化卡和PCI总线的I/O转换卡,转换成模拟信号,然后再通过接口放大器, 送入物理装置。物理装置对输入的信号做出响应后,将反馈的模拟信号和开关信号再通过 PCI总线的A/D转化卡和PCI总线的I/O转换卡,转换成数字信号,送入仿真程序中,完成整个闭环仿真。
[0184] 该步骤也可以由独立的物理装置接口程序并行处理,这时需要在电磁暂态仿真程序和物理装置接口程序中增加输入输出量的发送和接收。
[0185] 步骤18 :物理装置运行。
[0186]步骤 19 :设 t = t+dt。
[0187] 在该步骤19中将仿真计算的时间t增大一个步长dt。
[0188] 步骤20:判断增大后的时间是否大于总仿真时间,如果判断结果为“是”,结束整个仿真过程,如果判断结果为“否”,返回到步骤4,重复步骤4至18的仿真过程,计算下一时步的暂态过程。
[0189] 需要说明的是,对于交直流电力系统来说,对其进行分网或交直流分割的目的在于,第一,将涉及各个子网的计算任务分配到机群服务器的各个子机,由各个子机并行计算,以提高计算速度,做到电磁暂态实时或超实时仿真。机群服务器的各个子机,需在主机的统一控制下并列进行相应的计算,并交换信息。第二,如果不分网只进行交直流分割而将分割的各个部分的计算任务在一台机器上来执行(串行计算),也可以大幅度提高计算速度,这时因为将换流器分割开来后,换流阀的导通和关断将不再影响整个网络的拓扑结构,这样不需要将整个网络重新进行LU分解,只需要对换流器所在局部网络进行LU分解,会大大减少计算量。这是本发明的第二实施例中所要说明的内容。
[0190] 在图1所示的根据本发明的第一实施例中,步骤2、5、6、8、9、11、14和步骤17均是利用机群服务器的主机和多个子机实现各个子网的计算任务。
[0191] 图4示出了根据本发明的第二实施例的交直流电力系统分割并行电磁暂态数字仿真方法的流程图。
[0192] 步骤1 :交直流电力系统分割。
[0193] 在该步骤1中,将一个目标交直流电力系统分割为多个子网,子网可以是交直流混合网络,也可以是单独的交流、换流器和直流网络。
[0194] 步骤2 :判断该子网是否为交直流混合网络。
[0195] 在该步骤2中,判断该子网是否为交直流混合网络,即该子网是否由交流网络和换流器或交流网络、换流器和直流网络组成。如果判断为“是”,则前进到步骤3 ;如果判断为“否”,则前进到步骤5。
[0196] 步骤3 :对于交直流混合网络,求取子网内部交直流联络线电流。
[0197] 步骤4:对于交直流混合网络,根据子网内部交直流联络线电流,分别求取交流网络、换流器网络和直流网络节点电压。
[0198] 步骤3和4的目的是计算各子网端口点的等值阻抗矩阵。
[0199] 仍以图2所示系统为例,说明具体计算方法。设目标电力系统被分割为3个子网 S1、S2和S3,如图5所示。子网Sl含交流网络A及换流器1、2,子网S2含交流网络B,换流器3、4及直流网络,子网S3含交流网络C,因而子网Sl和子网S2都是交直流混合网络,子网S3为非交直流混合网络。
[0200] 对于子网Si,分别写出交流网络和换流器网络的端口方程如式(30)、(31)所示
[0201]
Figure CN101719182BD00181
[0202]
Figure CN101719182BD00182
:行,可形成子网内部交直流联络线方程
[0203] 由式(30)第一、二行和式(31)第—
[0204] Aiinn = b+ci (32)
[0205] 其中iim= [iQl,iQ2]T为子网内部交直流联络线电流,i = [_1,、,%12)^为子网间联络线电流,A、b、c如下式所示:
[0206]
Figure CN101719182BD00183
[0208] 令子网Sl中电流源置0,ia、iY、iM2依次分别置1或-1(电流正方向为流入该子网则置1,否则置-1),其它电流置0,先根据式(32)求出ial、ia2,再根据式(9)和下式 (33)分别求取交流和换流器网络节点电压,求出的端口点电压向量列的组合即为端口点等值阻抗矩阵。
[0209] Gcojjl2Vconl2 - n\Pcon\^a\ Jz niPconical ~ PconfiH^β12 = ^con\2
^3 (33)
[0210] 子网S2的端口点等值阻抗矩阵可类似求取。
[0211] 步骤5 :对于非交直流混合网络,直接求取端口点的等值阻抗矩阵。方法同本发明第一实施例中的步骤2。
[0212] 步骤6 :设定初始计算时间t = 0。
[0213] 步骤7 :判断此时刻是否有开关动作,或换流阀导通、关断等使网络结构发生变化的情况发生。
[0214] 在该步骤7中,判断此时刻各个子网中是否有开关动作或者换流阀导通或关断等情况发生。如果判断为“是”,则前进到步骤8 ;如果判断为“否”,则前进到步骤12。
[0215] 步骤8 :对有开关动作或者换流阀导通或关断等使网络结构发生变化的情况发生的子网,判断该子网是否为交直流混合网络;如果是则转至步骤9,否则转至步骤11 ;
[0216] 步骤9 :对于有开关动作或者换流阀导通或关断等使网络结构发生变化的情况发生的交直流混合网络,求取子网内部交直流联络线电流;
[0217] 步骤10 :对于有开关动作或者换流阀导通或关断等使网络结构发生变化的情况发生的交直流混合网络,根据子网内部交直流联络线电流,分别求取交流网络、换流器网络和直流网络节点电压;
[0218] 步骤9和10的目的是计算网络结构变化子网端口点的等值阻抗矩阵。
[0219] 步骤11 :对有开关动作或者换流阀导通或关断等使网络结构发生变化的情况发生的非交直流混合网络,直接求取端口点的等值阻抗矩阵;
[0220] 步骤9、10和11的计算方法与步骤3、4和5相同,注意此时只需对网络结构发生变化的子网进行处理。
[0221] 步骤12 :判断该子网是否为交直流混合网络。
[0222] 在该步骤12中,判断该子网是否为交直流混合网络,即该子网是否由交流网络和换流器或交流网络、换流器和直流网络组成。如果判断为“是”,则前进到步骤13 ;如果判断为“否”,则前进到步骤15。
[0223] 步骤13 :对于交直流混合网络,求取子网内部交直流联络线电流。
[0224] 步骤14 :对于交直流混合网络,根据子网内部交直流联络线电流,分别求取交流网络、换流器网络和直流网络节点电压。
[0225] 步骤13和14的目的是计算各子网的端口点等值电势。
[0226] 仍以图5所示系统为例,说明具体计算方法。
[0227] 对于子网Sl :
[0228] 令子网Sl中ia、iY、iM2 = 0,先根据式(32)求出ial、ia2,再根据式(9)和式 (33)分别求取交流和换流器网络节点电压,求出的端口点节点电压,即为端口点等值电势。
[0229] 子网S2的端口点等值电势可类似求取。[0230] 步骤16 :求取子网间联络线电流。
[0231] 根据各子网端ロ点等值阻抗矩阵和等值电势求取子网间联络线电流。
[0232] 仍以图5所示系统为例,说明具体计算方法。
[0233] 子网Si、S2、S3的端ロ方程如下:
[0234]子网 Sl
[0235]
Figure CN101719182BD00201
[0236] 子网 S2
[0237]
Figure CN101719182BD00202
[0238] 子网S3的端ロ方程不变,仍为式(7),其标识符改为S3,重写如下
Figure CN101719182BD00203
[0240] 将式(34)〜(36)合井,消去va,,vY,vM2,最后得到子网间联络线电流方程式 如下:
[0241] Ai = b (37)
[0242] 其中i = [ia,i0,iY,iM2]T为子网间联络线电流,矩阵A如下式所示:
[0243]
Figure CN101719182BD00204
[0244] 右端项b如下式所示:
[0245]
Figure CN101719182BD00205
[0246] 求解线性方程组(37),可求得子网间联络线电流。
[0247] 步骤17 :各子网根据得到的子网间联络线电流求取子网节点电压。
[0248] 步骤18 :各子网求解发电机、控制系统等微分方程。
[0249] 步骤19 :判断是否有用户自定义模型或MATLAB模型。
[0250] 在该步骤19中,判断各个子网中是否有用户自定义模型或MATLAB模型。如果判断为“是”,则前进到步骤20和21 ;如果判断为“否”,则前进步骤22。
[0251] 步骤20 :用户自定义模型、MATLAB模型输入输出信息处理。
[0252] 在该步骤20中,需要求取并发送各外接模型(用户自定义模型、MATLAB模型)输入变量值,并接收和处理各外接模型的输出变量值。
[0253] 步骤21 :用户自定义模型、MATLAB模型计算。
[0254] 在步骤21中,需要先接收用户自定义模型、MATLAB模型的输入变量值,然后根据用户自定义模型、MATLAB模型的各功能框计算公式和连接关系,进行模型计算,求取并发送用户自定义模型、MATLAB模型输出变量值。
[0255] 步骤21中的计算与其它步骤中的计算并行进行,由另外的独立程序(用户自定义模型计算程序、MATLAB模型计算程序)来实现。
[0256] 如果步骤21中模型的计算量不大,也可以将其与其它步骤串行进行计算,由同一个程序来处理,这时,发送和接收数据可以去掉,步骤20和步骤21合并为图4所示的步骤 21。
[0257] 步骤22判断是否有机电暂态接口。
[0258] 在该步骤22中,判断是否有机电暂态接口。如果判断为“是”,则前进到步骤23和 24 ;如果判断为“否”,则前进到步骤25。
[0259] 步骤23 :机电暂态接口输入输出信息处理。
[0260] 在机电暂态仿真的整数倍步长点,求取并向机电暂态仿真程序发送与机电暂态子网接口的边界点正、负、零序电压和电流,接收机电暂态仿真程序发送的边界点正、负、零序等值阻抗和电势。
[0261] 步骤M :机电暂态仿真计算。
[0262] 步骤M由另外的独立的机电暂态仿真程序来实现,先接收电磁暂态程序发送的边界点正、负、零序电压和电流,进行机电暂态一个步长的积分计算,然后求取并向电磁暂态仿真程序发送与电磁暂态子网接口的边界点正、负、零序等值阻抗和电势。
[0263] 步骤25 :判断是否有物理装置。
[0264] 在该步骤25中,判断是否有物理装置。如果判断为“是”,则前进到步骤沈和27 ; 如果判断为“否”,则前进到步骤28。
[0265] 步骤沈:物理装置输入输出信息处理。
[0266] 该步骤主要是进行数模和模数转换,将实时仿真的数字信号和开关信号分别通过 PCI总线的D/A转化卡和PCI总线的I/O转换卡,转换成模拟信号,然后再通过接口放大器, 送入物理装置。物理装置对输入的信号做出响应后,将反馈的模拟信号和开关信号再通过 PCI总线的A/D转化卡和PCI总线的I/O转换卡,转换成数字信号,送入仿真程序中,完成整个闭环仿真。
[0267] 该步骤也可以由独立的物理装置接口程序并行处理,这时需要在电磁暂态仿真程序和物理装置接口程序中增加输入输出量的发送和接收。
[0268] 步骤27:物理装置运行。
[0269]步骤 28 :设 t = t+dt。
[0270] 在该步骤28中将仿真计算的时间t增大一个步长dt。
[0271] 步骤四:判断增大后的时间是否大于总仿真时间,如果判断结果为“是”,结束整个仿真过程,如果判断结果为“否”,返回到步骤7,重复步骤4至观的仿真过程,计算下一时步的暂态过程。
[0272] 前述不分网(整个系统为1个子网)只进行交直流分割而将分割的各个部分的计算任务在一台机器上来执行的情况,是属于本发明的第二实施例的特殊情况。此时,流程图如图6所示,与图4中的流程图相比,去掉了若干步骤。
[0273] 步骤1 :交直流电力系统分割。
[0274] 在该步骤1中,将一个目标交直流电力系统分割为多个部分,包括交流网络部分、 换流器部分和直流网络部分,注意此时仍为1个子网。
[0275] 步骤2 :形成子网内部交直流联络线方程的矩阵A。方法同前述非特殊情况下的步骤3和4中子网内部交直流联络线方程(3¾的形成,但此时形成的联络线方程式中c = 0。
[0276] 步骤3 :设定初始计算时间t = 0。
[0277] 步骤4:判断此时刻是否有开关动作,或换流阀导通、关断等使网络结构发生变化的情况发生。
[0278] 在该步骤4中,判断此时刻是否有开关动作或者换流阀导通或关断等情况发生。 如果判断为“是”,则前进到步骤5 ;如果判断为“否”,则前进到步骤6。
[0279] 步骤5 :在有开关动作或者换流阀导通或关断等使网络结构发生变化的情况发生时,重新形成子网内部交直流联络线方程的矩阵A ;
[0280] 步骤6 :形成子网内部交直流联络线方程的右端项b ;
[0281] 步骤7 :求取子网内部交直流联络线电流;
[0282] 步骤8 :根据子网内部交直流联络线电流,分别求取交流网络、换流器网络和直流网络节点电压;
[0283] 步骤9 :求解发电机、控制系统等微分方程。
[0284] 步骤10 :判断是否有用户自定义模型或MATLAB模型。
[0285] 在该步骤10中,判断各个子网中是否有用户自定义模型或MATLAB模型。如果判断为“是”,则前进到步骤11和12 ;如果判断为“否”,则前进步骤13。
[0286] 步骤11 :用户自定义模型、MATLAB模型输入输出信息处理。
[0287] 在该步骤11中,需要求取并发送各外接模型(用户自定义模型、MATLAB模型)输入变量值,并接收和处理各外接模型的输出变量值。
[0288] 步骤12 :用户自定义模型、MATLAB模型计算。
[0289] 在步骤12中,需要先接收用户自定义模型、MATLAB模型的输入变量值,然后根据用户自定义模型、MATLAB模型的各功能框计算公式和连接关系,进行模型计算,求取并发送用户自定义模型、MATLAB模型输出变量值。
[0290] 步骤12中的计算与其它步骤中的计算并行进行,由另外的独立程序(用户自定义模型计算程序、MATLAB模型计算程序)来实现。[0291] 如果步骤12中模型的计算量不大,也可以将其与其它步骤串行进行计算,由同一个程序来处理,这时,发送和接收数据可以去掉,步骤11和步骤12合并为图4所示的步骤 12。
[0292] 步骤13 :判断是否有机电暂态接口。
[0293] 在该步骤13中,判断是否有机电暂态接口。如果判断为“是”,则前进到步骤14和 15 ;如果判断为“否”,则前进到步骤16。
[0294] 步骤14 :机电暂态接口输入输出信息处理。
[0295] 在机电暂态仿真的整数倍步长点,求取并向机电暂态仿真程序发送与机电暂态子网接口的边界点正、负、零序电压和电流,接收机电暂态仿真程序发送的边界点正、负、零序等值阻抗和电势。
[0296] 步骤15 :机电暂态仿真计算。
[0297] 步骤15由另外的独立的机电暂态仿真程序来实现,先接收电磁暂态程序发送的边界点正、负、零序电压和电流,进行机电暂态一个步长的积分计算,然后求取并向电磁暂态仿真程序发送与电磁暂态子网接口的边界点正、负、零序等值阻抗和电势。
[0298] 步骤16 :判断是否有物理装置。
[0299] 在该步骤16中,判断是否有物理装置。如果判断为“是”,则前进到步骤17和18 ; 如果判断为“否”,则前进到步骤19。
[0300] 步骤17 :物理装置输入输出信息处理。
[0301] 该步骤主要是进行数模和模数转换,将实时仿真的数字信号和开关信号分别通过 PCI总线的D/A转化卡和PCI总线的I/O转换卡,转换成模拟信号,然后再通过接口放大器, 送入物理装置。物理装置对输入的信号做出响应后,将反馈的模拟信号和开关信号再通过 PCI总线的A/D转化卡和PCI总线的I/O转换卡,转换成数字信号,送入仿真程序中,完成整个闭环仿真。
[0302] 该步骤也可以由独立的物理装置接口程序并行处理,这时需要在电磁暂态仿真程序和物理装置接口程序中增加输入输出量的发送和接收。
[0303] 步骤18 :物理装置运行。
[0304]步骤 19 :设 t = t+dt。
[0305] 在该步骤19中将仿真计算的时间t增大一个步长dt。
[0306] 步骤20:判断增大后的时间是否大于总仿真时间,如果判断结果为“是”,结束整个仿真过程,如果判断结果为“否”,返回到步骤7,重复步骤4至19的仿真过程,计算下一时步的暂态过程。
[0307] 采用本方法,其网络分割方式非常灵活,可以从交直流的任一连接点进行分网,例如,对于图2所示的系统,可以从i034处或iαρ iα2处或ia3、iα4处分网。
[0308] 本方法也适用于纯交流电力系统的分网并行求解。
[0309] 上面结合附图对本发明的具体实施方式进行了描述。应该注意的是,本发明不限于上述实施方式,在不脱离本发明的精神的前提下,本领域技术人员能够进行多种修改和变更。

Claims (9)

1. 一种交直流电力系统分割并行电磁暂态仿真方法,包括以下步骤:步骤1 :将一个目标交直流电力系统分割为多个子网,包括交流网络子网、换流器子网和直流网络子网;步骤2 :各子网求取端口点的等值阻抗矩阵; 步骤3 :设定初始计算时间t = 0 ;步骤4 :判断此时刻是否有开关动作,或换流阀导通、关断使网络结构发生变化的情况发生;如果是则转至步骤5,否则转至步骤6 ; 步骤5 :对子网求取端口点的等值阻抗矩阵; 步骤6 :各子网求取端口点的等值电势;步骤7 :对换流器采用节点法求解,或采用回路法求解,在换流器采用节点法求解情况下求取联络线电流,在换流器采用回路法求解情况下求取联络线变量,所述联络线变量是子网间联络线电流,或是子网间端口点电压;步骤8 :各子网根据节点法求解情况下得到的子网间联络线电流或回路法求解情况下得到的联络线变量求取子网节点电压,在换流器采用回路法求解情况下根据子网间联络线变量求取换流器连支电流;步骤9 :各子网求解发电机、控制系统微分方程;步骤10 :判断是否有用户自定义模型或MATLAB模型;如果是则转至步骤11和转至步骤12并在之后分别执行步骤13,否则转至步骤13 ;步骤11 :用户自定义模型、MATLAB模型输入输出信息处理; 步骤12 :用户自定义模型、MATLAB模型计算;步骤13 :判断是否有机电暂态接口 ;如果是则转至步骤14和转至步骤15并在之后分别执行步骤16,否则转至步骤16 ;步骤14 :机电暂态接口输入输出信息处理; 步骤15 :机电暂态仿真计算;步骤16 :判断是否有物理装置;如果是则转至步骤17和转至步骤18并在之后分别执行步骤19,否则转至步骤19;步骤17 :物理装置输入输出信息处理;步骤18 :物理装置运行;步骤19 :将仿真时间增加一个时间步长;重复执行上述步骤4至步骤19,直至到达总仿真时间为止。
2. 一种交直流电力系统分割并行电磁暂态仿真方法,包括以下步骤:步骤1 :将一个目标交直流电力系统分割为多个子网,子网是交直流混合网络,或是单独的交流、换流器和直流网络;步骤2 :判断该子网是否为交直流混合网络;如果是则转至步骤3,否则转至步骤5 ; 步骤3 :对于交直流混合网络,求取子网内部交直流联络线电流; 步骤4:对于交直流混合网络,根据子网内部交直流联络线电流,分别求取交流网络、 换流器网络和直流网络节点电压,从而计算出各子网端口点的等值阻抗矩阵,执行步骤6 ; 步骤5 :对于非交直流混合网络,直接求取端口点的等值阻抗矩阵; 步骤6 :设定初始计算时间t = 0 ;步骤7 :判断此时刻是否有开关动作,或换流阀导通、关断使网络结构发生变化的情况发生;如果是则转至步骤8,否则转至步骤12 ;步骤8 :对子网,判断该子网是否为交直流混合网络;如果是则转至步骤9,否则转至步骤11;步骤9 :对交直流混合网络,求取子网内部交直流联络线电流; 步骤10 :对交直流混合网络,根据子网内部交直流联络线电流,分别求取交流网络、换流器网络和直流网络节点电压,从而计算出网络结构变化子网端口点的等值阻抗矩阵,执行步骤12 ;步骤11 :对非交直流混合网络,直接求取端口点的等值阻抗矩阵;步骤12 :判断该子网是否为交直流混合网络,如果是则转至步骤13,否则转至步骤15 ;步骤13 :对于交直流混合网络,求取子网内部交直流联络线电流; 步骤14 :对于交直流混合网络,根据子网内部交直流联络线电流,分别求取交流网络、 换流器网络、直流网络节点电压,从而计算出各子网端口点的等值电势,执行步骤16 ; 步骤15 :对于非交直流混合网络,直接求取端口点的等值电势; 步骤16 :求取子网间联络线电流;步骤17 :各子网根据子网间联络线电流求取子网节点电压; 步骤18 :各子网求解发电机、控制系统微分方程求解;步骤19 :判断是否有用户自定义模型或MATLAB模型,如果是则转至步骤20和转至步骤21并在之后分别执行步骤22,否则转至步骤22 ;步骤20 :用户自定义模型、MATLAB模型输入输出信息处理; 步骤21 :用户自定义模型、MATLAB模型计算;步骤22 :判断是否有机电暂态接口,如果是则转至步骤23和转至步骤M并在之后分别执行步骤25,否则转至步骤25 ;步骤23 :机电暂态接口输入输出信息处理; 步骤M :机电暂态仿真计算;步骤25 :判断是否有物理装置;如果是则转至步骤沈和转至步骤27并在之后分别执行步骤观,否则转至步骤观;步骤沈:物理装置输入输出信息处理;步骤27 :物理装置运行;步骤28 :将仿真时间增加一个时间步长;重复执行上述步骤7至步骤观,直至到达总仿真时间为止。
3. 一种交直流电力系统分割并行电磁暂态仿真方法,其特征在于该方法是交直流电力系统为一个子网且只进行交直流分割,分割后的各个部分的计算任务在一台机器上来执行的情况,具体包括如下步骤:步骤1 :将一个目标交直流电力系统分割为多个部分,包括交流网络部分、换流器部分和直流网络部分;步骤2 :形成子网内部交直流联络线方程的矩阵A ; 步骤3 :设定初始计算时间t = 0 ;步骤4 :判断此时刻是否有开关动作,或换流阀导通、关断使网络结构发生变化的情况发生;如果是则转至步骤5,否则转至步骤6 ;步骤5 :重新形成子网内部交直流联络线方程的矩阵A ; 步骤6 :形成子网内部交直流联络线方程的右端项b ; 步骤7 :求取子网内部交直流联络线电流;步骤8 :根据子网内部交直流联络线电流,分别求取交流网络、换流器网络和直流网络节点电压;步骤9 :求解发电机、控制系统微分方程;步骤10 :判断是否有用户自定义模型或MATLAB模型,如果是则转至步骤11和转至步骤12并在之后分别执行步骤13,否则转至步骤13 ;步骤11 :用户自定义模型、MATLAB模型输入输出信息处理; 步骤12 :用户自定义模型、MATLAB模型计算;步骤13 :判断是否有机电暂态接口,如果是则转至步骤14和转至步骤15并在之后分别执行步骤16,否则转至步骤16 ;步骤14 :机电暂态接口输入输出信息处理; 步骤15 :机电暂态仿真计算;步骤16 :判断是否有物理装置,如果是则转至步骤17和转至步骤18并在之后分别执行步骤19,否则转至步骤19;步骤17 :物理装置输入输出信息处理;步骤18 :物理装置运行;步骤19 :将仿真时间增加一个时间步长;重复执行上述步骤4至步骤19,直至到达总仿真时间为止。
4.如权利要求1或2所述的交直流电力系统分割并行电磁暂态仿真方法,其特征在于通过将目标电力系统分割为交流、换流器和直流多个子网,进而将涉及各个子网的计算任务分配到机群服务器的各个子机并行计算,从而实现交直流电力系统的实时或超实时电磁暂态仿真。
5.如权利要求3所述的交直流电力系统分割并行电磁暂态仿真方法,其特征在于通过将目标电力系统分割为交流、换流器、直流多个部分,进而将各个部分的计算任务在一台机器上串行计算。
6.如权利要求1或2所述的交直流电力系统分割并行电磁暂态仿真方法,其特征在于其网络分割从交直流的任一连接点进行分网,对交流网络还进一步分网,分割成若干个交流子网。
7.如权利要求6所述的交直流电力系统分割并行电磁暂态仿真方法,其特征在于交流网络的分网方式是节点分裂法分网,或是分布参数线路解耦法分网。
8.如权利要求1-3之一所述的交直流电力系统分割并行电磁暂态仿真方法,其特征在于在该方法所模拟的电力系统中接入用户自定义模型或MATLAB模型所模拟的控制系统, 用于控制系统的结构设计和参数优化,以扩展对电力系统的仿真分析能力。
9.如权利要求1-3之一所述的交直流电力系统分割并行电磁暂态仿真方法,其特征在于可在该方法所模拟的电力系统中接入包括继电保护和安全自动装置、直流输电控制装置实际物理装置进行闭环仿真试验,以检验这些物理装置在实际电力系统中的效果。
CN2009102418641A 2009-12-11 2009-12-11 一种交直流电力系统分割并行电磁暂态数字仿真方法 CN101719182B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009102418641A CN101719182B (zh) 2009-12-11 2009-12-11 一种交直流电力系统分割并行电磁暂态数字仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009102418641A CN101719182B (zh) 2009-12-11 2009-12-11 一种交直流电力系统分割并行电磁暂态数字仿真方法

Publications (2)

Publication Number Publication Date
CN101719182A CN101719182A (zh) 2010-06-02
CN101719182B true CN101719182B (zh) 2012-01-18

Family

ID=42433756

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009102418641A CN101719182B (zh) 2009-12-11 2009-12-11 一种交直流电力系统分割并行电磁暂态数字仿真方法

Country Status (1)

Country Link
CN (1) CN101719182B (zh)

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101930487B (zh) * 2010-07-09 2012-06-06 南方电网科学研究院有限责任公司 一种电力系统的连续长时间稳定运行的混合实时仿真方法
CN102075135B (zh) * 2010-11-24 2014-03-12 中国电力科学研究院 以发电机内电势频率作为输入信号的稳定器的实现方法
CN102033993B (zh) * 2010-12-07 2013-04-10 中国电力科学研究院 大规模电力系统动态仿真中继电保护模型的构建方法
WO2012113096A1 (zh) * 2011-02-22 2012-08-30 Fan Shengtao 时域变换法
CN102355001B (zh) * 2011-09-26 2013-10-02 江西省电力科学研究院 一种兼顾稳定性与经济性的交直流系统直流落点选择的方法
CN102508967B (zh) * 2011-11-04 2013-11-06 南方电网科学研究院有限责任公司 一种大电网的电磁/机电暂态混合仿真平启动方法
CN102521452B (zh) * 2011-12-14 2014-01-29 中国电力科学研究院 一种大电网合环计算系统
CN102708225B (zh) * 2012-04-13 2014-08-06 南方电网科学研究院有限责任公司 一种交直流大电网电磁暂态仿真的分片调试方法
CN102708250A (zh) * 2012-05-10 2012-10-03 天津大学 基于树形分层双向迭代的电力系统数字混合仿真方法
CN102779209B (zh) * 2012-06-19 2014-09-10 中国电力科学研究院 一种大步长多开关实时定点的电磁暂态仿真方法
CN102750416B (zh) * 2012-06-19 2014-03-12 中国电力科学研究院 一种含开关特性电路的电磁暂态仿真的拓扑分网方法
CN102788998B (zh) * 2012-08-01 2014-12-03 东华理工大学 全空间瞬变电磁法物理模型实验用的测量装置
CN103605829B (zh) * 2013-10-22 2016-08-17 华中科技大学 对交直流混联电网进行电磁暂态仿真的等值建模方法
CN103970591B (zh) * 2014-05-28 2017-04-12 山东大学 一种Matlab‑ATP集成自动化调用方法
CN104865842B (zh) * 2015-03-25 2019-01-04 南方电网科学研究院有限责任公司 适用于稳定控制装置接口的混合仿真系统及混合仿真方法
CN104732033A (zh) * 2015-03-31 2015-06-24 河南行知专利服务有限公司 一种电力系统机电暂态过程的仿真分析方法
CN106294897B (zh) * 2015-05-22 2020-03-17 中国电力科学研究院 一种适用于电磁暂态多时间尺度实时仿真接口的实现方法
CN105608245A (zh) * 2015-10-29 2016-05-25 贵州黔燎原电气科技有限公司 大规模配电网多cpu独立解算的建模方法
CN105322540B (zh) * 2015-11-12 2019-02-22 中国电力科学研究院 一种交直流大电网电磁暂态模型的稳态运行方式建立方法
CN106886616B (zh) * 2015-12-15 2020-06-16 中国电力科学研究院 一种大规模电磁暂态电网仿真的自动分网方法
CN108062032A (zh) * 2016-11-05 2018-05-22 南京理工大学 一种基于通信的电力分布式控制与仿真系统
CN108181835B (zh) * 2017-12-30 2019-11-12 贵州大学 一种电力系统直流电网实时仿真方法及装置
CN109948185B (zh) * 2019-02-01 2020-10-16 全球能源互联网研究院有限公司 一种电力系统的解耦仿真方法

Also Published As

Publication number Publication date
CN101719182A (zh) 2010-06-02

Similar Documents

Publication Publication Date Title
Guo et al. Comprehensive real-time simulation of the smart grid
Xue et al. A simple direct method for fast transient stability assessment of large power systems
Zhu et al. Adaptive power flow method for distribution systems with dispersed generation
CN102255307B (zh) 一种配电网可靠性评估的分层等值方法
Pereira et al. A decomposition approach to automated generation/transmission expansion planning
CN102081867B (zh) 一种全数字高压直流输电培训仿真系统
WO2017198237A1 (zh) 配电网多时间尺度数模混合仿真系统、方法及存储介质
CN101860037B (zh) 一种高压直流输电系统网侧谐波电流的确定方法
Liu et al. Skeleton-network reconfiguration based on topological characteristics of scale-free networks and discrete particle swarm optimization
Nath et al. Coherency based system decomposition into study and external areas using weak coupling
CN103632313B (zh) 一种基于pmu数据的电网动态可观方法
El-Kady et al. Dynamic security assessment utilizing the transient energy function method
CN103683284B (zh) 一种含分布式电源的配电网三相不平衡潮流计算方法
CN103279807B (zh) 一种恶劣天气下电网静态风险评估方法
US10198018B2 (en) Systems and methods for convex relaxations and linear approximations for optimal power flow in multiphase radial networks
El-Marsafawy et al. A new, fast technique for load-flow solution of integrated multi-terminal DC/AC systems
Ma et al. Economic analysis of real-time large scale PEVs network power flow control algorithm with the consideration of V2G services
CN101256550B (zh) 复杂电网相位同步并行化评估系统
KR100726024B1 (ko) 전력계통 실시간 시뮬레이터를 이용한 신재생에너지전원용인버터의 계통연계시험 장치 및 방법
CN103323688B (zh) 一种电力系统的谐波状态估计方法及设备
Ali et al. Communication modeling for differential protection in IEC-61850-based substations
Zhang et al. A distribution short circuit analysis approach using hybrid compensation method
CN104218604B (zh) 一种基于网络等值法的配电网可靠性分析方法及系统
CN104112165A (zh) 基于多目标离散粒子群的智能配电网故障恢复方法
Chen et al. Loop frame of reference based three-phase power flow for unbalanced radial distribution systems

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: STATE ELECTRIC NET CROP.

Effective date: 20130129

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20130129

Address after: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15

Patentee after: China Electric Power Research Institute

Patentee after: State Grid Corporation of China

Address before: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15

Patentee before: China Electric Power Research Institute