|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
基于MATLAB的电力系统继电, w& m$ L1 p k: e. O
保护仿真研究; ^4 w# W4 V( W D
胥杰1,陈峦1,2
( i8 e+ I! U$ f(1.新疆大学电气工程学院,新疆乌鲁木齐830047;5 Q+ q# M: u6 X. }) Y; i# L
2.电子科技大学自动化工程学院,四川成都610054)
% g8 T7 _% }! g2 ?% `* Z摘要:基于MATLAB的仿真技术可以辅助继电保护系统的分析和设计,针对电力系统继电保护领域的核心内容,% ^$ E/ w* P5 o( z I4 c; X
构建了系统仿真模型.举例说明了电力系统故障、零序电流保护和变压器纵差保护等仿真实例的实现方案。在运行
2 }2 M4 T7 Y6 |2 }这屿仿真实例的基础上.对仿真结果进行了深入分析。通过仿真实例证明了电力系统继电保护仿真的有效性和可6 z1 D5 \: _5 a) k# }( v( K; J# L
行性。( |# ^- X3 J; c; I% R f1 n; a# [
关键词:电力系统;继电保护;仿真;MATLAB
' Q3 r" _9 t7 E7 W b" WMATLAB—-based Simulafion of Power System Relay Protection
" E3 E7 e0 P; |Xu Jiel.Chen Luanl,2
1 t; Z. H6 K: ~; h& q8 p(1.College of Electrical Engineering,Xinjiang University,Urumqi Xinjiang 830047;
9 n6 @$ a% [! j5 {- E0 R R2.College of Automation Engineering,University of Electronic Science and Technology,Chengdu Sichuan 610054)9 N! Q# C+ ?( ~7 V- Y" c* b
Abstract:MATLAB-based simulation technology can support the analysis and design of relay protection systems.A$ D8 ^8 F/ s4 O& v7 r0 o1 M
simulation model is built for the study of power system relay protection.As an example,the power system fault simulation,' h* ~$ i: e$ {# t; R1 g
zero-sequence current protection simulation and transformer differential protection simulation are presented herein.The ease
, x. i. z0 `1 ]& F5 wstudies show that the simulation of power system relay protection is effective and feasible.
* g- s7 G$ O( J- O( X6 dKey Words:power system;relay protection;simulation;MATLAB) H1 q8 ?& c0 e7 `( G* S
中图分类号:TM743 文献标识码:A 文章编号:0559—9342(2010)03-0084-03& s" x9 G/ E$ s8 A; u# K/ q$ @
为了设计出高性能的继电保护产品.常常要进
6 f ]. J1 t2 n# q6 }行继电保护模拟实验.传统的继电保护实验是在真
+ c; E. }) j% u# g! x9 t实的物理模型上进行的.实验系统复杂。实验成本
% v& s" }5 D% ^5 B1 b高,效率低.灵活性和通用性差。电力系统继电保9 S- p+ y- L, V
护数字仿真是利用计算机软件仿真电力系统故障情
0 `7 l( }, U; j% R况下电气量变化的特征.模拟继电保护装置处理和+ e* E6 k, D: B# M H9 C* a
动作的过程.具有安全、经济、可重复、不受环境+ K2 s2 s8 ~7 U/ k: H' X8 k) V
限制、研究开发周期短、效率高等优点【ll。通过对各
: ?) k% |$ V$ r% a种不同继电保护技术的仿真。可使研究者形象地观) L0 g! ~3 C. {$ }: ^: j
察到电力系统从正常状态到发生故障期间电气量的6 q5 n3 o- n8 B# I! n
变化情况以及继电保护装置中每一部分处理和动作
% a; [; Q9 s& V, r# k! S# S5 c2 s) H的过程。利用计算机软件模拟真实保护设备装置的
; z+ S7 F7 I) g7 {, O F运行情况.能及早发现和解决设备运行中可能I叶{现& Q" G. N# Y" s7 u( z
的问题.因此.电力系统继电保护计算机仿真具有
# d" i. l7 V9 _6 h1 A, W十分重要的现实意义嘲。8 n+ ~3 K, M5 B. L. B9 ]
本文设计了电力系统故障、空载合闸励磁涌流、
- k1 S }3 o; D( o* D: I四Water Power VoL3&No.3! n, N' f( n8 e/ w$ _* l
单侧电源相问短路的电流保护、零序电流保护、变9 q2 v0 ?8 @* t5 @( Z; w; I
压器纵差保护、微机保护算法等仿真实例的实现方
' s9 h7 F* \! z/ h案,在运行这些仿真实例的基础上,对仿真结果进
' j |1 M) D" J1 A" Z5 k- f行了分析M。
, i2 |$ t+ j0 B, H+ A7 q1电力系统故障仿真% D6 {5 e, k. ^" f5 D
本研究采用双电源供电的系统模型.电压等级& {+ \& {, ?8 u! L' n4 D* i
为220kV。如图1所示。该模型中同步发电机的容$ Z' M6 i5 }2 f. a; i
量为500 MV·A,电压为13.8 kV,频率为50 Hz。- W8 O/ [3 L* T+ o
i相变压器的容量为500 MV·A,Dll/Yg接线形式,+ T3 q1 Z8 C9 y9 E
收稿目期:2009-05—23
2 F8 m+ Q" m1 c3 Q基金项目:新疆大学2l世纪高等教育教学改革丁程资助项目% p3 p/ Q/ T$ {) S6 g
(XJU2008JGY21)
% ]8 h' b5 {8 g ? o6 K作者简介:胥杰(1985一),男,四川绵阳人,主要从事电力系统
! Q, |$ J$ @4 J0 T9 c继电保护方面的研究.
% {. V! B! G; S0 k万方数据
' T% ~9 s/ \) }; Z三:::2:: :=:二:::::::::::=:::三:兰=:!二==巨l二三===:幽●●■■●●■■■●●●■■●■■●●■■■●■■■■■■■■■■●■■■■■■■●●●■●■●■■■■■■■■●■■■●■■■■■■■■■■■■■●■■■■■■■■■■■■■■■■●■■■■■■■■■■●■■■■■■■■■■■■■■●■■■●■●●■■■■■■■●●■■■■■■■■■■■■■■■■■●■■■■●■■■■■●■●■■■●■●●■■■■■■●■■■■●■■一■■●■■●●■■■■■■■■■●■■●■■■■■■■●■■●●一
/ q- ]2 x" O {/ K圈1电力系统故障模型
& u4 k/ G# H$ d7 I) w# T频率为50 Hz。150 km线路的正序阻抗为0.01165+. ]( G/ j. _, Z" o
i0.0008679刚km,对地电容为13.41x10。9 F/km;100
* O7 @# c0 z* _! }( h6 o; zkm线路的正序阻抗为0.01165+jo.0008679 ll/km.7 m* o5 C' I- O9 t$ T
对地电容为13.41×10母F/km。负荷l的电压为220
4 F1 b7 v+ E# u* xkV.有功负荷为220x106/250 W.无功负荷为2003 V/ M. K4 Q+ q7 @
W:负荷2的电压为220 kV。有功负荷为220x106// d- w4 C" |( [1 U9 g
250W。无功负荷为200W;负荷3的电压为220 kV.) x- d0 i! C1 a8 c6 z6 }
有功负荷为220x106/250 W.无功负荷为200 W:负
, ~9 U9 X; H |荷4的电压为13.8 kV。有功负荷为220x106 W:无/ T# b1 A4 m2 k4 n7 Q6 C' Z
功负荷为0 W。* ~9 }- Q0 |5 z) k3 ]
在构建的故障系统MATLAB仿真模型中.通过4 v# k% O3 D G0 l( E& O) ^
故障模块可设置i相短路、两相短路、两相短路接/ s2 Z7 O# \9 Z: r- f: d' x7 F
地和单相短路接地故障。以单相接地短路故障为例。
9 f: t3 N. z& y9 z1 s设置i相短路元件参数为A相接地短路.运行仿真( q2 I& n; q1 J N- X5 _: j2 v
模型.得到线路单相接地短路时母线的短路电压及
) h q D) Y, g: _: ?/ [+ ^9 E电流波形如图2所示。
, g+ s0 _0 j) D+ n400
# d" X2 {0 o- p$ Z6 |之200
, ~/ j7 E8 g& c- N% K2 H日兰0
# D8 \- q; J1 }1 X& {9 R! P4 g3 G脚一20(t
& ?7 H B9 |6 H) q3 W6 K一40(}3 D: b0 I% E7 z( ~0 x) f
200+ c8 e/ c, w8 |
≤0
f' u7 u% X7 N: Z0 _喜-n()(1
: O) ?8 L/ Y. V—4()(J
; ]1 T$ ?7 S C4 s9 C0 0.05 0.1 0.15 0.2 0.25 0.32 ?: R2 d* L0 j% c
时问/s4 j8 {1 y' ?* _3 E2 U( [
【一。土。^..;⋯..,:。i.一' e. q& [: m8 r, X1 ~' B
i j羔F!
2 W3 _. @& _- W9 O0 ~" k# }$ I、, 一、’,‘一~$ K2 b- {' h7 O% w( G1 A" N0 K
0 0.0.5 0.1 0.15 0.2 0.25 0.3
& g# K: E5 O' b时fHJ/s, u' n: V4 T, Z! i3 U* Z
图2 电力系统单相接地短路时故障仿真- [7 |) }9 t' Y/ p3 l Y. A( W
波形显示。在0~0.03 s期间,线路工作在稳定
4 A; K6 ]$ Y$ g: [2 i& w# T状态。三相电流和电压对称;在0.03 s时发生A相. y* G: b* k% O x
接地短路.A相电压基本为0.B相和C相电压也5 F! ^0 y% a: O1 X
相对减小.故障相A相的电流迅速上升为短路电
1 i% f+ V$ W. ~5 }2 V流.B相和C相电流也相对增大;0.08 s时切除' z; l1 u* g7 T- p; n
故障.i相电压电流经暂态后达到新的稳定状态.+ B* Z/ W1 ^, X% Z/ L
并重新恢复=三相对称运行工作状态。
6 G9 ]- }: ]8 J6 r2 y0 Q8 \) A; |同理.可以设置其他的短路故障类型,并观察3 N" V$ {& M. H( d
其仿真波形.也可以对仿真模型中的i相故障模块
. |1 V$ C4 { o4 R F" K设置接地电阻参数.通过观察对应的仿真波形.分
" @ o. e( z# l6 i0 q析接地电阻对电力系统故障的影响。$ b. n: p" c% L7 D
2零序电流保护仿真
7 O# o( ]$ V Q7 y( n; h用Sireulink和PSB模块库构建一个220 kV的
$ j7 e1 o# W: _1 W: N0 M单侧电源供电系统模型.模拟输电线路发生接地故9 p6 t) q4 F1 n9 h" m+ L0 F
障。同时,建立一个零序电流保护及单相重合闸装
9 s+ a: `3 U4 E: _& y( _% B置的仿真模型.通过该系统可十分方便地进行不同& N9 X$ G ^3 U! n: B( S! }6 J6 m
接地电阻及电弧情况下的故障分析以及继电器动作( ?; t' Q: Z8 a! X1 D8 E; ^
特性的模拟。以单相接地短路为例。将故障模块设( L {9 W3 f4 r2 L- r# s0 \- K7 o
置为A相接地短路,仿真结果如图3所示。 c! u1 b! o- u! J9 g
之IOXl0、
, W+ r8 ?( \9 E0 E& _8 S7 ]3 c“04 i1 Z, H& f g$ A9 T$ q: l
粤一IO×10。* _2 h1 u3 q( H4 N! G9 p% X' k0 @. b
0 0.OB 0.1 0.15 0.2 0.25 0.3* A9 e) O! w, r! ]- C
时删/s
( c9 }9 Q1 Y0 W' u( ^委一⋯1 00⋯0。IBCA。.。。。.。。.,。。.:。.。。。.]。8 b6 c5 J, t. O+ y5 E' C- U
时间/s 争03E三!I三三豆[I三量L——二_二上—二———L————_工————_=士====:====d- k! {+ I9 F( N/ m) n$ u" O
篓0 0.05 0.1 0.15 0.2 0.25 0.3
7 I' |+ n/ W. g" p2 h& \时问/s6 u. M7 q% b% {" F) L
图3零序电流保护仿真
1 n: Z/ O' X6 J) N, u9 L9 n9 \9 x7 Z8 M由图3可以看出.在0.03 s的时候线路发生单
. {: q: a. F c% P; t2 M* D% I. z相接地故障.A相电流发生急剧变化.电压变得不
% s6 h& X+ |' E对称.继电器动作后经0.05 s保护出口跳闸.跳闸后/ t$ Y, M; X7 w& Q# m7 i7 n
经O.3 s自动重合闸.此时故障仍存在.保护只用9 |, _) q9 [4 X- a( O
0.01 s立即发跳闸命令。从而完成了从跳闸到重合
4 |* d& \* s! j4 ]6 a8 B" u; }. X/ C闸再跳闸的整个过程。, J w% b4 ^/ e1 G" o/ c9 s! Z
3变压器纵差保护仿真
' J0 M& O" L, L5 o" O3 v* A采用双侧电源供电模型.电压等级为735 kV/
$ A: I6 n [. r* V H% J315 kV.变压器两侧安装有断路器,如图4所示。
4 ~7 l; _& z; X: G# L, X# i该模型中i相电源侧的电压为735 kV.频率为50
0 L8 c% s$ Y: a/ Q' Q6 \5 yHz:低压侧的电压为315 kV。频率为50 Hz。三相
3 S1 D! u( \' Y7 e% q/ r; u% ~- k变压器的SN=250 MV·A,U。√以N=735 kV/315 kV,1 D) D4 O* P3 i& b
接线方式Y.dll。在变压器纵差保护的MATLAB( J1 l% D& ~6 Q8 ?, F C: u5 S9 P
仿真模型中。通过故障模块可设变压器区内、区外
& N/ U5 {$ h; @" _" V的各种故障。以区内j相短路为例.设置故障元件
$ x' x0 y S% H: z$ p+ p为A、B、C三相短路,短路点设置在内部(两边电
9 i# _4 p/ U0 |% F- l0 x+ z3 g( |流表以内).得到变压器纵差保护内部i相故障模, M# R3 E; Z6 M2 a
型。运行该仿真模型,即可得到电压、电流和跳闸
: T; J/ O7 w) w6 g3 [& n0 d* J波形.如图5所示。# s4 u# ]. S2 I1 k7 E- B- n
由图5可见,在0~0.03 s期间.线路T作在稳
' M/ s9 N1 C6 h, V- L. Q3 V8 M定状态.i相电流和电压对称。在0.03 8时发生三
- c( I, |- R# @. q' }相短路.三相电压变小.i相电流迅速上升为短路9 K3 A, E, E( p, I
电流.i相电压和电流对称。说明i相短路为对称
6 }# B* ]; f3 R& n+ I短路。在0.13 s时,变压器纵差保护动作,切除故3 i9 ~/ `" `3 H
图4双电源系统双绕组压器模型
2 G1 [4 P# M* |- Q" n# ^Water Power V01.36.No.3四: i* P5 V! x4 m/ a
叮㈨斟一一》祝O ◇" V" X3 l1 }' }( ?' {/ g
万方数据 Y0 k \5 u/ B O, \: h k' q
掣吾-l1OX0时刈矗茎赢丽嫩苜贰萨而=畜=菌=百=玄丽甙甄而言网赢
( u. Y" z0 H3 [2 q6 e, k# ~时间/s 委蠹匿耍亟匦亟三0 0.02 0.04 o.06 o.08 QI n12 Q14 n16 Q18 Q2
! E' z/ ^8 S5 U) E0 [" g1 mI}|,Ifnq/s 羹n5E三三三三玉雪蚕0 0.02 0.04 n06 Q08 Ql Q 12 Q14 Q 16 n18 n29 c+ G) D" R& u! @3 z$ }0 j6 L1 g
时闯/s: v# e# s, J+ {( f1 O
圈5变压器纵差保护的仿真结果
2 c" P- V9 [; L5 Q1 T: W障,电压为电源电压,电流几乎为0。5 W3 I9 t( {, Y" A2 ~
4结语
1 }: r$ U' c% z2 b! f9 ]基于计算机技术的电力系统继电保护数字仿真
/ V i( @7 T- O, u可以辅助继电保护系统的分析和设计.将MATLAB( R* M T" x. O( | ?3 ]4 m n
仿真技术应用于电力系统继电保护研究中具有十分
7 G2 {4 V3 M: y& g* a' A& G9 W重要的现实意义。针对电力系统继电保护技术的核3 s# z6 o8 O ~; X" x
心内容.构建了继电保护系统的MATLAB仿真模/ s5 h9 i4 X& B9 S% M i3 j3 Y5 f
型。除所介绍的电力系统故障仿真、零序电流保护
% N* p: P/ j& `9 A; D仿真和变压器纵差保护仿真外,模型还可以对空载
# t, W( j% J, |' G! a+ k8 B合闸励磁涌流仿真、单侧电源相间短路的电流保护- b. E5 r& l- p: @1 R
仿真和微机保护算法仿真等进行实际仿真运行,实! F1 a$ o; E `( o- Z
践表明将MATLAB应用于电力系统继电保护数字仿
6 F; n7 c8 k" G' W2 r9 _ z真中是有效的和可行的。
) J8 a( l: l8 h参考文献:
1 @& ]$ z$ Y. n. m: S8 d& j[1]冯小玲,郭袅.谭建成.继电保护仿真系统的现状及其应用[J].; r; }0 c& G8 y7 t9 M
广西电力.2004(6):50-56.% z6 s& D/ c$ s* K: m' U6 u" m4 K! S
[2]盛义发,唐耀庚。苏泽光,等.基于MATLAB的电力系统故障
3 e- l: f% M5 q9 a5 B9 L的仿真分析[J].南华大学学报(理T版).2003,17(4):45—51.! A- I& Y( y5 H9 v* R1 K8 ^
[3]杨晓军,沈勇环,郭征,等.开放式的继电保护动态特性仿真* m* n2 M) {: w, Z# ^
系统[J].电力系统及其自动化学报,2003,15(6):l-4.
6 y- O9 G' n( J) w[4] 章志刚,李兴源.基于MATLAB开发平台的继电保护仿真系统( {* D6 k5 X6 T" d2 I
[J].四川I电力技术,2003(2):43—46.
. R% U# o" a4 @' f0 f[5]韩笑,徐曦.陈卓平.基于Matlab与VB数据交换的继电保护
3 O( D S w7 u+ {7 n, T( j/ [4 D仿真[J].电力自动化设备,2006,26(5):92—95.; I6 }3 }/ T) p- ~% N9 L
[6]黄绍平,李永·咯.基丁:MATLAB的变压器牢载合闸瞬变过程仿0 T/ `3 f% t! \% }# B7 T
真研究fJ].继电器,2004,32(8):19—21.
$ N/ a3 ~3 j- N[7]张伏生。寇强,刘芳,刘沛津,等.线路微机保护装置仿真系
6 I1 p1 e& }8 [, F* M# M; z统的设计与实现【J].继电器,2003,31(4):76-79.
* x+ o. j# L% {/ b6 [$ i(责任编辑高瑜)
; j3 V0 M, C8 h0 |$ U* _! e- i-■-——●———+--—+—-—●—-—+—-—+—-—+—·—+—-—+—·—+—·+-+·—+--+-+-+——+一—+——卜—+一+-'-I---—+一-+-+-+一+-+-+-+一+-+-+-+-+-+-+-+-+-—●。·—+一-—+——+一+·—●一·+·0 e/ Q! Z' v4 v4 w1 I% g9 I
(上接第68页)直线PB与A断面的夹角为Ot,则直! D, Y8 S* V! N9 c3 i1 d
线PB与A断面上任何一点的夹角均为Ot.因而直; V' s) F' U' {6 ^7 R% z
线咫与断面的圆心的夹角也为01.。假定C(Ⅳ,E,1 k& I0 ]2 K& ~7 h- h# X* g
日)为隧洞断面圆周上任意一实测点,则其对应断' ~. o: S* J4 d+ L
面所对应圆心角ot=arctan【(Ⅳ_Ⅳo)/(H一舶)】。9 A$ y* ]$ A' F, N/ f3 T
B
+ o: R* X7 J5 o9 Y: k0 d(』I。磊,%)/ ]2 z F; ]) o2 ^0 `' c
图3断面圆心对应角计鼻示霖
y5 X) S: |2 C0 H( Z4 [: O J2.2.4建立数学模型
! \- X* r3 c) I$ C, c①角度ct=arctan【(Ⅳ-Ⅳo)/(H-Ho)】;②断面圆
1 {8 S1 u7 @8 u# B4 [! |心角坐标O(X,Y,z),X=No一38sinct,Y=O。压
2 N) u; e6 g% k0 `Ho一38cosct;⑧断面实测半径尺=【(Ⅳ—x)2+(层一y)2+0 W' ]9 y) c2 x( T. T1 B0 F) W
(H—z)21比;④断面理论半径:L=7.25一od50。
1 T/ h1 K! |$ e/ E, \* q1 c* l2.3放样实施) J* Y) a! k1 b
把全站仪所测的坐标Ⅳ.E。日输入计数器。计8 x1 A R0 q! N- E, c0 M# g
算所测点的实际半径和理论半径.比较实际半径和
0 a, V& u: h7 G理论半径.当实际半径与理论半径之差为正时。应" U# {$ R0 B! o5 [7 n- r
向内移动棱镜.移动数值为其差值,反之向外移动
5 I. D! a$ D/ y. J1 u圈Water POll.1et VoL36.No.,
, ?1 z0 ^* }; U9 C5 c棱镜。反复移动,直到位置达到精度要求。8 C! `- G3 q, g2 C/ f
3结语: g# C* v+ M+ ~! C& U
溪洛渡水电站是围家重点工程.GPS、电子全/ l+ p1 t$ y0 J7 q* M7 F& e5 K7 D0 J
站仪、断面仪等当代工程测昔技术在洞群施工测冒
# B- x2 n7 N( s: i/ ]5 ?中得到了广泛的应用。随着测绘科技的发展,隧洞
- l% b* s/ F1 N8 U9 V5 e: U' ~. ]工程施丁测量技术和方法也会不断地创新和发展., f* L0 S* W0 y. z7 O- y
将会有更多的施T测量新技术、新方法应用于隧洞" I7 S% U. m6 ]- h$ h+ X) W# f
工程施T测量中。4 R) ?0 p7 z0 z {3 ?8 y" {
参考文献:
) Q' {* w6 w; L# ^( |* F7 V0 A[1]贺国宏.桥隧控制测量[M].北京:交通H{版社,1999.
) q/ x; K. l) Z9 n" o[2] 吴来瑞,邓学才.建筑施r测量手册[M].北京:中国建筑1=业
7 t( N6 Z# e; O( i& y7 r2 ]$ V出版社.2007.
/ ~/ N; k: l! v$ J0 n0 X1 `[3]王景海.直线与曲线交点的算法、程序及算例[J].勘察科学技$ [ E3 Y% [$ z/ u" x" l# I" C4 p, O
术。2009(2):24—27.3 w' v' o/ t2 Q% o' a
[4]陈健.MATLAB在变形监测数据处理中的应用[J].城市勘测,
K" q, X& i! X2 Q# _2009(2):130-133.
, m$ M, @2 y+ {[5]刘文生.隧道圆拱线施T测量位移法研究及应用[J].城市勘测,
* ?* W) x1 {! {: {# ? b6 R" j2007(5):59-60.
2 z4 _! ?7 P- l Y/ m! G' K5 u[6]刘丹源,刘钢.水电站水轮发电机尾水管安装过程中测量方法# @: R: S1 _9 b- h# A- u" n
探讨[J].测绘通报,2008(2):69—70.! _3 }. i9 \) f5 D
[7] 赵吉先,万崔辉.贯通横向误差处理新方的法探讨[J].测绘通报, H9 U# Q* ?# H, T% x6 @3 U9 l
2008(1):10-11. @' Z3 {* l7 }& ~! }+ `
(责任编辑杨健)
. l6 @. S" G6 U* N万方数据 |
|