基于行波理论的输电线路故障诊断方法研究【Simulink复现】

最新推荐文章于 2024-11-13 00:01:16 发布

电力程序小学童

最新推荐文章于 2024-11-13 00:01:16 发布

阅读量1.2k

点赞数 25

分类专栏：故障测距文章标签： matlab simulink 行波测距输电线路故障

本文链接：https://blog.csdn.net/zhangxd212489/article/details/142817221

版权

故障测距专栏收录该内容

1 篇文章 0 订阅

订阅专栏

主要内容

该程序复现文章《基于行波理论的输电线路故障诊断方法研究》，利用故障行波（包括电压行波和电流行波）在阻抗不连续点折射和反射原理，通过simulink模型采集故障信息（包括输电线路两侧电流和电压信息），然后采用m文件对故障信息进行处理分析，通过Karenbauer变换得到故障电流的线模、零模分量以及输电线路两端电压线模分量，通过小波变换对线模分量进行分解，标定行波波头，从而计算得到故障距离。该仿真模型复现效果好，注释清楚，有详细的操作文档，方便参考和研究！

模型研究

1.整体模型

simulink模型结构和原文一致，原文中输电线路电压为110kV，该模型中电压为220kV。输电线路总长度为200km，故障点设置距离L端90km，距离R端110km，AC相间短路运行结果显示，“A、C两相短路但不接地故障，且距离L端89775.8835米”，测量误差仅为0.112%，可见精度非常高。（程序可以按照操作文档自行设置单相短路、两相接地、三相短路但不接地和两相短路但不接地等故障类型进行设置）

原文模型如下：

2.故障相判定

故障相判定原理主要是对比故障后与故障前的行波信号，将故障后采集到的行波信号减去故障前正常的行波信号，从而得到故障附加行波信号，通过相模变换，得到附加行波的模值，从而对故障类型进行精准辨识。

故障附加电流行波信号计算方法：

通过相模变换得到附加行波的模值：

从而按照以下方式判定系统故障类型：

3.详细的模型说明文档

部分代码

% 设定故障发生在0.05s，那么之前的电流就是正常运行时候的
IL_A1_one=IL_A(1:0.02*F);%提取一个周期的正常电流，电流频率50Hz，采样频率是10^6，所以需要0.02*10^6个采样点；
IL_B1_one=IL_B(1:0.02*F);
IL_C1_one=IL_C(1:0.02*F);

% 形成完整的正常电流：总仿真时长0.15s，就是7.5个周期,最后面多加1位因为仿真0.15s会出现150001个，补齐使得矩阵维度一样
IL_A1=[IL_A1_one;IL_A1_one;IL_A1_one;IL_A1_one;IL_A1_one;IL_A1_one;IL_A1_one;IL_A1_one(1:0.01*F+1)];
IL_B1=[IL_B1_one;IL_B1_one;IL_B1_one;IL_B1_one;IL_B1_one;IL_B1_one;IL_B1_one;IL_B1_one(1:0.01*F+1)];
IL_C1=[IL_C1_one;IL_C1_one;IL_C1_one;IL_C1_one;IL_C1_one;IL_C1_one;IL_C1_one;IL_C1_one(1:0.01*F+1)];
M_I1=max(IL_A1);
% 形成故障附加电流
IL_A2=IL_A-IL_A1;%按照论文式（15）
IL_B2=IL_B-IL_B1;
IL_C2=IL_C-IL_C1;
% Karenbauer变换得到故障电流的线模和零模分量
IL_0=1/3*(IL_A2+IL_B2+IL_C2);%这里下标不带m的就是诊断用的
IL_alpha=1/3*(IL_A2-IL_B2);
IL_beta= 1/3*(IL_A2-IL_C2);
IL_gama= 1/3*(IL_B2-IL_C2);
%%  测距所需变量
VL_A=out.VL_A.signals.values;% 由模型导入工作区的L端三相电压，重新命名
VL_B=out.VL_B.signals.values;
VL_C=out.VL_C.signals.values;
VR_A=out.VR_A.signals.values;% 由模型导入工作区的R端三相电压，重新命名
VR_B=out.VR_B.signals.values;
VR_C=out.VR_C.signals.values;
%首先进行Karenbauer变换得到L、R端电压两个线模分量
VL_alpha=1/3*(VL_A-VL_B);%L端电压行波线模alpha分量，可以看出是不含C相故障信息的，所以C相单相接地不能用
VL_beta= 1/3*(VL_A-VL_C);%L端电压行波线模beta分量，可以看出是不含B相故障信息的，所以B相单相接地不能用
VR_alpha=1/3*(VR_A-VR_B);
VR_beta= 1/3*(VR_A-VR_C);

%小波基选择db7，对两个线模分量进行分解
[cVL_alpha,lVL_alpha] = wavedec(VL_alpha,7,'db7');%分解L端电压行波线模alpha分量
[cVR_alpha,lVR_alpha] = wavedec(VR_alpha,7,'db7');%分解R端电压行波线模alpha分量
[cVL_beta,lVL_beta] = wavedec(VL_beta,7,'db7');%分解L端电压行波线模beta分量
[cVR_beta,lVR_beta] = wavedec(VR_beta,7,'db7');%分解R端电压行波线模beta分量