|
|
|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
基于MATLAB的电力系统继电
* c' w( Y" ]% k保护仿真研究$ A9 N$ E. p$ z4 i/ T
胥杰1,陈峦1,2, W$ Q! S+ F+ w: ]0 U" H; I, `
(1.新疆大学电气工程学院,新疆乌鲁木齐830047;
% l5 p$ L' J+ L( `2.电子科技大学自动化工程学院,四川成都610054)3 V/ T. ^& v, a. R/ P
摘要:基于MATLAB的仿真技术可以辅助继电保护系统的分析和设计,针对电力系统继电保护领域的核心内容,
1 p" }. q V% u2 Y$ S- X构建了系统仿真模型.举例说明了电力系统故障、零序电流保护和变压器纵差保护等仿真实例的实现方案。在运行- q) n: d& c) o& a$ U6 z+ m9 x
这屿仿真实例的基础上.对仿真结果进行了深入分析。通过仿真实例证明了电力系统继电保护仿真的有效性和可
3 a& e3 q& D$ ^* P: G0 X行性。# t W- E4 p, P% S- U2 b1 c- M& M% a
关键词:电力系统;继电保护;仿真;MATLAB/ r& y! M6 J s. e5 U3 p
MATLAB—-based Simulafion of Power System Relay Protection5 D2 F9 T+ I9 ~) W( P' D y
Xu Jiel.Chen Luanl,2
; Z7 c% K9 L0 d1 b! o(1.College of Electrical Engineering,Xinjiang University,Urumqi Xinjiang 830047;
* i/ y) S$ E# o, L2.College of Automation Engineering,University of Electronic Science and Technology,Chengdu Sichuan 610054)1 d6 A/ B0 T1 X) L, y2 Y8 ?1 H4 F
Abstract:MATLAB-based simulation technology can support the analysis and design of relay protection systems.A" g/ [2 S0 |2 o6 X
simulation model is built for the study of power system relay protection.As an example,the power system fault simulation,
! B. Y. ^ r; ?zero-sequence current protection simulation and transformer differential protection simulation are presented herein.The ease
/ l4 P5 l% s: n" Xstudies show that the simulation of power system relay protection is effective and feasible.8 C( N* k' h8 V6 t8 H4 k
Key Words:power system;relay protection;simulation;MATLAB" m8 h: u; }. T8 j
中图分类号:TM743 文献标识码:A 文章编号:0559—9342(2010)03-0084-03" E+ q7 L: S' H. ?1 {9 a, h
为了设计出高性能的继电保护产品.常常要进. F: l, C' q N; a: C
行继电保护模拟实验.传统的继电保护实验是在真
0 ~( @( |4 {" W0 w; U实的物理模型上进行的.实验系统复杂。实验成本
% K2 u; Y; d: p+ w0 \1 x高,效率低.灵活性和通用性差。电力系统继电保
0 x- h+ I0 t6 k+ V) e3 s: I护数字仿真是利用计算机软件仿真电力系统故障情
( }* f$ L; e; V) L5 ~况下电气量变化的特征.模拟继电保护装置处理和
, l& L* m0 Y# j! n- |" B# F动作的过程.具有安全、经济、可重复、不受环境" X4 r. o2 `: y/ _7 ^) K
限制、研究开发周期短、效率高等优点【ll。通过对各4 p8 M' }" |+ G% q! ^
种不同继电保护技术的仿真。可使研究者形象地观
: W0 C+ J S- i% K# N; R察到电力系统从正常状态到发生故障期间电气量的
. d0 U2 `3 @/ p7 r* R- ~变化情况以及继电保护装置中每一部分处理和动作
. Y$ B' X4 I+ ?( u9 Q的过程。利用计算机软件模拟真实保护设备装置的
; t) |3 u6 A+ S8 s) p- D& X$ p运行情况.能及早发现和解决设备运行中可能I叶{现
) y$ A' F5 C7 s9 U的问题.因此.电力系统继电保护计算机仿真具有
# w: X* ], k/ h" M( i% K, B2 E十分重要的现实意义嘲。
! p% L. W" Z- V0 c本文设计了电力系统故障、空载合闸励磁涌流、8 Y; d" z/ ~- X5 `2 c+ G$ ]
四Water Power VoL3&No.3
. L7 r6 t0 Z/ I7 q( l单侧电源相问短路的电流保护、零序电流保护、变4 q7 L2 E0 d+ r" ]
压器纵差保护、微机保护算法等仿真实例的实现方
* ^( u8 y: {) E0 q0 n" d3 I案,在运行这些仿真实例的基础上,对仿真结果进
& R* v1 |% `! \6 c行了分析M。8 W) @ q+ z5 ?. X5 a9 U
1电力系统故障仿真
1 {( {5 l5 n c; I5 h6 \$ Y$ {本研究采用双电源供电的系统模型.电压等级: v, _/ X& f# r, F
为220kV。如图1所示。该模型中同步发电机的容
' r) G( @: { K8 S8 q量为500 MV·A,电压为13.8 kV,频率为50 Hz。
5 [1 z1 ?: w& t. W5 N, L2 `: p7 O- Zi相变压器的容量为500 MV·A,Dll/Yg接线形式,
9 R" S: l, k% @$ E( r收稿目期:2009-05—23
1 q) A T8 y2 C基金项目:新疆大学2l世纪高等教育教学改革丁程资助项目( |; G* D, D# o K
(XJU2008JGY21)
# k2 @, ^ U+ u6 {1 ?作者简介:胥杰(1985一),男,四川绵阳人,主要从事电力系统8 I1 d! K3 j: r3 { n) K
继电保护方面的研究.5 d0 e/ g$ c7 ]5 L' k4 ]6 t
万方数据
) A* g/ ^* k0 k$ f三:::2:: :=:二:::::::::::=:::三:兰=:!二==巨l二三===:幽●●■■●●■■■●●●■■●■■●●■■■●■■■■■■■■■■●■■■■■■■●●●■●■●■■■■■■■■●■■■●■■■■■■■■■■■■■●■■■■■■■■■■■■■■■■●■■■■■■■■■■●■■■■■■■■■■■■■■●■■■●■●●■■■■■■■●●■■■■■■■■■■■■■■■■■●■■■■●■■■■■●■●■■■●■●●■■■■■■●■■■■●■■一■■●■■●●■■■■■■■■■●■■●■■■■■■■●■■●●一: H. }% U) ?" q' M# t( Y Q
圈1电力系统故障模型
& c, x' v4 f, O- }7 D* u6 ~频率为50 Hz。150 km线路的正序阻抗为0.01165+
( [) c% P7 Y% A' ]2 ui0.0008679刚km,对地电容为13.41x10。9 F/km;100+ T" ^6 A# q( }* C' y+ h$ d
km线路的正序阻抗为0.01165+jo.0008679 ll/km.
$ O0 l) v& ]& ^3 Q- E7 B1 y1 h对地电容为13.41×10母F/km。负荷l的电压为220
$ U) R) S" S! @5 w7 f3 s: H* a' vkV.有功负荷为220x106/250 W.无功负荷为200
, o! S0 N p7 X6 j- P) u. J; {W:负荷2的电压为220 kV。有功负荷为220x106/: V$ c, m' ^4 F$ A' `! ]
250W。无功负荷为200W;负荷3的电压为220 kV.
% Q* K/ O7 ]* b5 r( }有功负荷为220x106/250 W.无功负荷为200 W:负5 K8 z. X8 `9 N, a8 E8 I( p% ?
荷4的电压为13.8 kV。有功负荷为220x106 W:无
+ e; j! R9 D5 E/ {功负荷为0 W。0 }' c' `! w; p
在构建的故障系统MATLAB仿真模型中.通过
$ f, g! d5 @1 }# a# c+ Q5 c C故障模块可设置i相短路、两相短路、两相短路接' \' M. U* x& |
地和单相短路接地故障。以单相接地短路故障为例。
+ L. |8 H5 d6 G7 R& c- E设置i相短路元件参数为A相接地短路.运行仿真
4 @- x1 V- i0 a9 l模型.得到线路单相接地短路时母线的短路电压及' g$ K; M$ D$ R" O7 c
电流波形如图2所示。
4 k. ]. {" c4 m$ j* l% {400
+ a2 k$ D8 l+ z; D之200
7 Z2 ^' D, v% f" F! ^- r日兰0( w, W( x* F: |
脚一20(t
: n" c6 h- O- N一40(}& p# ?- w! `3 @- N+ I& a4 s
2005 ~; @0 y5 T r" s
≤0
. e C7 ], F' E- J" j喜-n()(1
6 C& j. b, B) p, | I5 b/ B—4()(J: e; i7 Q& ?7 ]& E
0 0.05 0.1 0.15 0.2 0.25 0.3; h2 M/ D T$ K) J
时问/s
1 ^2 Z% d* @( ~ w& k8 J/ _【一。土。^..;⋯..,:。i.一
6 @& \: F4 c+ M& G8 F9 C4 gi j羔F!
: U3 L' ?! E2 x、, 一、’,‘一~
* \+ R: J+ X$ Y+ M( r4 {6 z" q0 0.0.5 0.1 0.15 0.2 0.25 0.37 B) ? i( L& C
时fHJ/s4 i0 Z. C+ P# Z" X4 x
图2 电力系统单相接地短路时故障仿真$ ^- s! s! r( r2 Z2 f% F
波形显示。在0~0.03 s期间,线路工作在稳定" I) U! L2 r) K2 t9 E# G' P
状态。三相电流和电压对称;在0.03 s时发生A相4 ?3 V, L) S( `5 ]3 S& `
接地短路.A相电压基本为0.B相和C相电压也
" \. b- y8 g, Q- }* s相对减小.故障相A相的电流迅速上升为短路电4 H0 @4 _1 q" k; H. ~
流.B相和C相电流也相对增大;0.08 s时切除
# o$ M: h- F% `8 }4 q! A$ U2 A故障.i相电压电流经暂态后达到新的稳定状态. @; E1 b! H! \! l
并重新恢复=三相对称运行工作状态。7 T+ D- {1 }5 U# R" Y$ Z( T
同理.可以设置其他的短路故障类型,并观察! B- g1 \+ Y8 n1 D
其仿真波形.也可以对仿真模型中的i相故障模块
% B1 I2 j' e9 { y: _9 W设置接地电阻参数.通过观察对应的仿真波形.分
; Z5 E0 x$ m7 S y C析接地电阻对电力系统故障的影响。) S% j" \: z+ I1 X/ \: }& y6 k
2零序电流保护仿真8 `( o2 C0 b+ E& ]4 r9 d+ a; ^
用Sireulink和PSB模块库构建一个220 kV的1 C4 w) y% J# M- \& C5 M8 k
单侧电源供电系统模型.模拟输电线路发生接地故
* O8 ^+ q4 w1 y: n5 C障。同时,建立一个零序电流保护及单相重合闸装
6 r+ G2 N9 c# c7 L置的仿真模型.通过该系统可十分方便地进行不同) p+ E) V5 h- \1 F
接地电阻及电弧情况下的故障分析以及继电器动作4 {) Y9 X; H, ?, a- O* i
特性的模拟。以单相接地短路为例。将故障模块设7 I4 [" C% m/ p9 P0 B8 e& [ n) S
置为A相接地短路,仿真结果如图3所示。
& L8 | x, G6 V8 R; d之IOXl0、% ? S# P3 D& g$ y( v- J
“0
( t' ]* j2 O& I) n3 ~粤一IO×10。: E1 m) K. n" ^) S. D ?: A
0 0.OB 0.1 0.15 0.2 0.25 0.3
+ D1 f) X! k$ `! u! K* f; q( `时删/s
! f! ^, F9 J5 A4 w3 X委一⋯1 00⋯0。IBCA。.。。。.。。.,。。.:。.。。。.]。
# H; z5 N5 D, T* i/ `' w时间/s 争03E三!I三三豆[I三量L——二_二上—二———L————_工————_=士====:====d- e8 S, _' d2 s0 E* g' _
篓0 0.05 0.1 0.15 0.2 0.25 0.39 M. H+ j9 t1 R4 j6 b/ ~# N B0 m
时问/s0 ?' B, g( P7 {8 e0 N
图3零序电流保护仿真9 d% q7 n4 s/ o( U5 ~
由图3可以看出.在0.03 s的时候线路发生单: `. @# w& ?7 r \3 X# @3 D
相接地故障.A相电流发生急剧变化.电压变得不& h9 q- H+ v) Q" K0 W% Q
对称.继电器动作后经0.05 s保护出口跳闸.跳闸后
+ @: u7 b5 _7 l# \; |经O.3 s自动重合闸.此时故障仍存在.保护只用3 V6 Q5 D" r9 _! J- x6 \, d6 D
0.01 s立即发跳闸命令。从而完成了从跳闸到重合
& L3 P+ Q0 j4 ]闸再跳闸的整个过程。# H# m) A, o0 Z* t. j! L
3变压器纵差保护仿真' Z o0 x% M2 z
采用双侧电源供电模型.电压等级为735 kV/
7 o* O$ ~# M0 h+ \& r315 kV.变压器两侧安装有断路器,如图4所示。
+ Y3 m+ E& ]: \+ E该模型中i相电源侧的电压为735 kV.频率为50
" S% }8 x8 g3 [4 rHz:低压侧的电压为315 kV。频率为50 Hz。三相
& A# A: m" t3 h/ ]) _6 ?! a. I变压器的SN=250 MV·A,U。√以N=735 kV/315 kV,
& }/ a0 q$ i: W6 M( F3 W接线方式Y.dll。在变压器纵差保护的MATLAB I; I' P6 M( i% H$ g
仿真模型中。通过故障模块可设变压器区内、区外7 w+ h. w2 L' U+ \ b
的各种故障。以区内j相短路为例.设置故障元件$ F1 T# v, U/ M7 w" ? w
为A、B、C三相短路,短路点设置在内部(两边电5 ?% q. v( Q" ]" |
流表以内).得到变压器纵差保护内部i相故障模
) c; l) \' z7 E0 z型。运行该仿真模型,即可得到电压、电流和跳闸7 e* L' e) l. ~! }5 C5 p9 n
波形.如图5所示。
0 U B4 Y- i% v2 y由图5可见,在0~0.03 s期间.线路T作在稳
7 _% f( w5 h$ r3 p定状态.i相电流和电压对称。在0.03 8时发生三
' W" y' U5 d: \, S相短路.三相电压变小.i相电流迅速上升为短路
, d# I# `% d" o( u$ T P电流.i相电压和电流对称。说明i相短路为对称% Q+ O9 |- K: r+ K# b8 u( w5 a5 I/ ?
短路。在0.13 s时,变压器纵差保护动作,切除故
6 m1 j( k0 j& v6 N, H' C) d图4双电源系统双绕组压器模型
& p: m" h1 {) f3 L+ b1 a7 fWater Power V01.36.No.3四
* a( s9 I" e3 E; O% r' V1 p" x叮㈨斟一一》祝O ◇( B5 I' W, F4 m! N! f
万方数据! _# U! `7 V/ i, X8 P# B
掣吾-l1OX0时刈矗茎赢丽嫩苜贰萨而=畜=菌=百=玄丽甙甄而言网赢$ O' d8 A _9 d) b
时间/s 委蠹匿耍亟匦亟三0 0.02 0.04 o.06 o.08 QI n12 Q14 n16 Q18 Q2
7 b9 p, T! X4 w/ s- lI}|,Ifnq/s 羹n5E三三三三玉雪蚕0 0.02 0.04 n06 Q08 Ql Q 12 Q14 Q 16 n18 n2+ Q! ?2 L: k4 x2 n2 Z) r) a
时闯/s
7 G! Y s( x7 s1 T# g1 f圈5变压器纵差保护的仿真结果
; s- M8 B H" A0 J障,电压为电源电压,电流几乎为0。
' c8 j3 a4 W; x, p9 n. G( a" x% l4结语( D+ U7 ?! R: \$ k3 Y6 y' m
基于计算机技术的电力系统继电保护数字仿真3 ?3 F' t1 d; M' e5 l: m$ s
可以辅助继电保护系统的分析和设计.将MATLAB$ Y4 B' s! f5 `; ^ G- A
仿真技术应用于电力系统继电保护研究中具有十分
. U: {/ O: Y6 E* z& a5 J: J3 Q2 P; }重要的现实意义。针对电力系统继电保护技术的核1 m L' N9 K, b! O
心内容.构建了继电保护系统的MATLAB仿真模$ ?/ U3 C K* f4 P& [$ M4 M
型。除所介绍的电力系统故障仿真、零序电流保护
+ r% I) v0 L3 r' M" R) a' l仿真和变压器纵差保护仿真外,模型还可以对空载
2 t$ F F$ h' l7 ^合闸励磁涌流仿真、单侧电源相间短路的电流保护( P8 A d- b3 F: F1 q
仿真和微机保护算法仿真等进行实际仿真运行,实
$ C( I! Y" [* o% j践表明将MATLAB应用于电力系统继电保护数字仿
' y% J7 A- b' v0 n真中是有效的和可行的。
^8 t! n- u% u( A! f参考文献:% u2 J* ~- {# N1 {4 U
[1]冯小玲,郭袅.谭建成.继电保护仿真系统的现状及其应用[J].; ?$ j% q" g- [$ q: v* A
广西电力.2004(6):50-56.- [( x/ P' I0 E- b5 x
[2]盛义发,唐耀庚。苏泽光,等.基于MATLAB的电力系统故障& t' |$ [2 \3 I$ z/ T5 z
的仿真分析[J].南华大学学报(理T版).2003,17(4):45—51.
: _3 ~: g% G: e1 N! {7 ^ K[3]杨晓军,沈勇环,郭征,等.开放式的继电保护动态特性仿真. H* [2 j2 J1 D2 ]) g5 i* D
系统[J].电力系统及其自动化学报,2003,15(6):l-4.
; K& @* M- t! b7 x: ]& z* X[4] 章志刚,李兴源.基于MATLAB开发平台的继电保护仿真系统
8 E$ Z) C8 d; w6 N3 v$ X6 T[J].四川I电力技术,2003(2):43—46.
( H2 b' n* w* }0 S2 H: L[5]韩笑,徐曦.陈卓平.基于Matlab与VB数据交换的继电保护
- o+ L3 A2 W: ?$ \( t仿真[J].电力自动化设备,2006,26(5):92—95.& k, E6 Y4 Q0 ~) [) m( c1 O. B
[6]黄绍平,李永·咯.基丁:MATLAB的变压器牢载合闸瞬变过程仿+ Z/ [% n# S- t+ i f% N
真研究fJ].继电器,2004,32(8):19—21.
" q" h u+ T _5 v9 C$ B% L[7]张伏生。寇强,刘芳,刘沛津,等.线路微机保护装置仿真系
* v6 S: d$ n. ?6 a Y统的设计与实现【J].继电器,2003,31(4):76-79.8 Q4 ], ?. [- {! V: t5 i
(责任编辑高瑜)# t* c0 b6 k- G ^
-■-——●———+--—+—-—●—-—+—-—+—-—+—·—+—-—+—·—+—·+-+·—+--+-+-+——+一—+——卜—+一+-'-I---—+一-+-+-+一+-+-+-+一+-+-+-+-+-+-+-+-+-—●。·—+一-—+——+一+·—●一·+·) ~( h. B2 I! q# }) m
(上接第68页)直线PB与A断面的夹角为Ot,则直- l6 f, f ~& K7 P* d/ P
线PB与A断面上任何一点的夹角均为Ot.因而直
" I5 q) B3 ~ N. x# W线咫与断面的圆心的夹角也为01.。假定C(Ⅳ,E,
. U I- b2 x. Q% p日)为隧洞断面圆周上任意一实测点,则其对应断
5 n# o$ H& |; \7 l4 X, W. Q# d面所对应圆心角ot=arctan【(Ⅳ_Ⅳo)/(H一舶)】。! {4 _& X! A$ X6 f3 v* V
B
8 d+ d( n9 A: i0 _(』I。磊,%)8 J3 S% B- P2 T
图3断面圆心对应角计鼻示霖
, n* c, A; V/ \& z7 D3 j# b2.2.4建立数学模型
4 \6 S( Z d& s# d8 S6 k①角度ct=arctan【(Ⅳ-Ⅳo)/(H-Ho)】;②断面圆
0 L; b- U! V- E; Z3 X: ~9 @' S% G l& z心角坐标O(X,Y,z),X=No一38sinct,Y=O。压
8 g4 Z' J$ g6 O, L9 \. U8 A' NHo一38cosct;⑧断面实测半径尺=【(Ⅳ—x)2+(层一y)2+5 {: A, U# V f; H2 ]
(H—z)21比;④断面理论半径:L=7.25一od50。
* x, L1 m# h7 |4 A4 Z( l" q2.3放样实施
2 i W4 b( N% u# c把全站仪所测的坐标Ⅳ.E。日输入计数器。计+ z3 o# {3 v* g3 z2 l) D1 m) C
算所测点的实际半径和理论半径.比较实际半径和: l% ^2 _' Z; }! B
理论半径.当实际半径与理论半径之差为正时。应
$ P9 Q. i$ k3 c, p3 m4 L2 ~2 ~' A% e向内移动棱镜.移动数值为其差值,反之向外移动# q9 a) G! D% n* ?- X8 N
圈Water POll.1et VoL36.No.,* L/ P2 a* F1 V, ?
棱镜。反复移动,直到位置达到精度要求。5 V4 s) G' Z0 N" K6 c& j, k K: A
3结语
9 t4 s2 e/ j0 x; q溪洛渡水电站是围家重点工程.GPS、电子全" c5 [0 U7 h, P6 a6 F# H f% G
站仪、断面仪等当代工程测昔技术在洞群施工测冒. W" f( h3 \7 g: D/ W2 B+ `8 ~2 l
中得到了广泛的应用。随着测绘科技的发展,隧洞
2 X l! G! ]5 t* c; y R工程施丁测量技术和方法也会不断地创新和发展.$ N' y: U& x9 u, r ?( ]- S
将会有更多的施T测量新技术、新方法应用于隧洞* X, N. }% {" }! {- G3 ]: o7 `
工程施T测量中。8 o- v( z9 A* I2 g
参考文献:7 q3 V) i5 O" H( a
[1]贺国宏.桥隧控制测量[M].北京:交通H{版社,1999.
. Y- }* _3 q+ X[2] 吴来瑞,邓学才.建筑施r测量手册[M].北京:中国建筑1=业
% L9 o- @& X( U/ s出版社.2007./ ~; I. x5 x, K; o$ @& F$ @
[3]王景海.直线与曲线交点的算法、程序及算例[J].勘察科学技! T! Q' ?! l) r% R0 x
术。2009(2):24—27., C, j: E' i. A
[4]陈健.MATLAB在变形监测数据处理中的应用[J].城市勘测,
! q7 E: D$ }1 m3 T8 q6 ]2009(2):130-133.
4 U) \( p4 X+ [! N% ~* G[5]刘文生.隧道圆拱线施T测量位移法研究及应用[J].城市勘测,
, _) x3 T0 k6 V# \% o2007(5):59-60.
2 J% J, m, M% y6 X% N _[6]刘丹源,刘钢.水电站水轮发电机尾水管安装过程中测量方法
3 b7 l+ p$ f+ i# j6 H- V( F4 q探讨[J].测绘通报,2008(2):69—70.9 v5 M! x; x" W) ~" K
[7] 赵吉先,万崔辉.贯通横向误差处理新方的法探讨[J].测绘通报,
! a( u2 x( B9 S3 U2008(1):10-11.! A" q% X) v/ m" M( R+ ?5 M
(责任编辑杨健)6 c4 [% k# z& d& b" U# a
万方数据 |
|