微机保护算法综合仿真Word格式文档下载.doc
《微机保护算法综合仿真Word格式文档下载.doc》由会员分享,可在线阅读,更多相关《微机保护算法综合仿真Word格式文档下载.doc(28页珍藏版)》请在冰豆网上搜索。
![微机保护算法综合仿真Word格式文档下载.doc](https://file1.bdocx.com/fileroot1/2022-10/4/76bc2673-e3b3-4483-bd10-8b6dda86aea1/76bc2673-e3b3-4483-bd10-8b6dda86aea11.gif)
零序阻抗:
线路对地正序电容:
;
线路对地零序电容:
M、R侧等值系统的参数为:
。
2、建立电力系统暂态仿真模型
利用MATLAB7.0POWERLIB工具箱搭建模型,见图2.
图2MATLAB模型图
三相故障模块中的暂态仿真时间被设置为0.1s开始故障,0.2s结束故障,因此在电力系统暂态仿真模型窗口的Simulink菜单下的SimulationParameters中可设置仿真开始时间为0s,仿真结束时间为0.2s,即0~0.1s时间段为仿真系统正常运行阶段,0.1~0.2s时间段为仿真系统故障阶段。
仿真操作既可在MATLAB命令窗口中运行相应的仿真命令完成,也可以直接在电力系统暂态仿真模型窗口执行Simulation菜单下的Start命令完成。
更简单的方法是,在暂态仿真模型窗口执行工具条上的图标完成电力系统暂态仿真模型的仿真。
2.2R-L模型解微分方程算法
当忽略线路的分布电容时,从故障点到保护安装处的线路可用一个电阻R1和电感L1的串联电路来近似。
见图3
图3 输电线路R-L模型
当线路中某点发生故障,可用微分方程表示
(2-1)
R和L为未知数,令,为第i时刻的算子。
在时刻、有等式:
(2-2)
将(2-2)式写成线性方程组:
(2-3)
计算中,可采用采样值差分法求出导数,然后求解L1和R1。
三相故障模块被设置为AB两相接地短路故障,暂态仿真时间为0.1s开始故障,0.2s结束故障,采样时间间隔为,即每工频周期采样12点(对应工频频率为50Hz)。
保存到文件模块采样时间参数与三相故障模块采样时间参数相同。
完成电力系统暂态仿真模型的实际仿真后,得到仿真数据文件ab1.mat。
2.21绘制M端电压电流曲线
图4基于MATLAB仿真电流电压波形图
从图4中可以看出该系统在第0.1秒时发生故障,A相B相电流电压波动较大,且并无明显衰减趋势,可粗略判断是AB相发生故障,过渡电阻很小。
下面做具体分析:
图5AB两相接地短路示意图
图5中,为接地过渡电阻,仿真中设为金属性短路。
边界条件为
(2-4)
(2-5)
用序量表示的边界条件为
(2-6)
得(2-7)
因为,由(1-2)式得
(2-8)
且,将(1-4)式代入(1-5)式中可化简为
(2-9)
以C相为参考量,由(2-7)可知正序和负序是并联的,(2-9)得零序和串联,故K点AB相接地短路复合序网如图6所示
图6K点AB相接地短路等效复合序网
附:
ab1.m
closeall
loadxiti1.mat;
k=m’%将m阵转置
t=k(:
1);
va=k(:
2);
vb=k(:
3);
vc=k(:
4);
%分别提取k阵的第2,3,4列
ia=k(:
5);
ib=k(:
6);
ic=k(:
7);
%分别提取k阵的第5,6,7列
figure
subplot(211);
%绘制三相的电压波形
plot(t,va,'
r'
t,vb,'
b'
t,vc,'
g'
);
xlabel('
t/ms'
ylabel('
v/V'
legend('
va'
'
vb'
vc'
subplot(212);
%绘制三相的电压波形
plot(t,ia,'
t,ib,'
t,ic,'
xlabel('
i/A'
ia'
ib'
ic'
)
2.22R-L模型算法综合仿真
图7M点到短路点电阻和电抗曲线图
经仿真计算,结果如下表所示
表一两相接地短路时仿真计算结果表
计算结果
计算值/
实际值/
相对误差/%
R
2.0593
1.273
61.77
L
26.5711
29.32
9.31
附:
AB2.m
clear
closeall
loadxiti1.mat;
k=m'
t0=k(:
t=t0;
tmax=max(t);
tmin=min(t);
[xy]=size(k);
T=0.02;
N=round((x*T)/(tmax-tmin));
m=k(:
2:
7);
%提取矩阵k第2列到第7列元素
va=m(:
vb=m(:
vc=m(:
ia=m(:
ib=m(:
ic=m(:
i=(ia-ib);
v=(va-vb);
%R-L模型算法(微分方程算法)
[a,b]=size(i);
i1=[0
i(1:
a-1)];
i2=[0
0
i(1:
a-2)];
v1=[0
v(1:
v2=[0
0
v(1:
r=((i-i1).*(v1+v2)-(v+v1).*(i1-i2))./((i-i1).*(i1+i2)-(i+i1).*(i1-i2));
x=pi*((v+v1).*(i1+i2)-(i+i1).*(v1+v2))./((i-i1).*(i1+i2)-(i+i1).*(i1-i2))/N
%控制阻抗轨迹显示范围
m=[t,r,x];
t1=0.08;
t2=0.14;
a1=t1*N/0.02;
a2=t2*N/0.02;
m1=m(a1:
a2,:
t=m1(:
r=m1(:
x=m1(:
%绘制电阻和电抗曲线图
figure
subplot(221);
plot(t,r,'
k-o'
ylabel('
r(t)'
xlabel('
subplot(223);
plot(t,x,'
ylabel('
x(t)'
subplot(222);
plot(t,r,'
axis([t1,t2,-20,20]);
%截取横轴0.08到0.14,纵轴-20到20的部分
subplot(224);
plot(t,x,'
xlabel('
axis([t1,t2,-10,40]);
%截取横轴0.08到0.14,纵轴-10到40的部分
%R0=0.1273/km,L0=0.9337mh/km
%X0=2*pi*f*L0=0.293/km;
%L=100kM;
Zk=[(R0+jX0)*100]=1.273+j29.3
Tukey低通滤波器具有较短的暂态时延,所以在微机距离保护中得到了应用。
其设计过程如下:
(1)Tukey模拟低通滤波器的冲击响应为
(2-10)
式中,为截止频率。