【行波故障测距】基于Simulink仿真三相电故障行波获取

一、仿真环境：

MATLAB2021b

Simulink

二、Simulink电路搭建

打开MATLAB，在界面上端进入Simulink：

 新建模型，在画布空白处双击鼠标左键，搜索对应模型安放并连线，仿真图如图：

 三相电电源设置如图：

其中峰值电压为110Kv，频率为50Hz。 

三相输电线设置如图：

故障模块设置如图：

其中故障类型为AB->G，AB双向接地故障。

输出模块配置如图：

 输出变量为Vabc，输出类型为二维数组。

配置停止时间为0.06s后，开始仿真：

仿真结束后，配置好的输出变量会自动输出到MATLAB工作区中。

三、MATLAB仿真

成功运行Simulink后，相应数据将会返回至工作区，此时MATLAB界面应显示：

之后开始进行编写代码： 

%提取故障发生前的部分数据和故障发生后的部分数据
%故障类型为AB->G，发生时间为0.02s，采样间隔ts=5us，故故障发生在0.02s/ts = 4000点处
Vm=out.Vabc(4001:8000,:);
Im=out.Iabc(4001:8000,:);

Vn=out.Vabc(1:4000,:);
In=out.Iabc(1:4000,:);

%求变化量
dV=Vm-Vn;
dI=Im-In;

%构建凯伦贝尔变换矩阵，并且对变化量进行解耦
Q=1/3*[2,-1,-1;0,1,-1;1,1,1;];
Vx=Q*dV';
Ix=Q*dI';

%根据公式计算行波，其中参数来自仿真中输电线配置
Z=sqrt(0.9337e-3/12.74e-9);
uf=0.5*(Vx(1,:)+Ix(1,:)*Z);
ur=0.5*(Vx(1,:)-Ix(1,:)*Z);

%输出波形
plot(uf);
hold on
plot(ur)

正向行波计算公式为：

反向行波计算公式为：

其中阻抗Z的计算公式为：

L与C分别为输电线单位长度上的感抗和容抗，可以在Simulink进行设置

编写完成并运行仿真，其输出结果如图：

突出显示其前一段，得到如下图所示：

 至此，完成行波的获取仿真。