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

       本文旨在使用PSCAD搭建最简单的三相电路故障模型，并且对双相接地型故障（AB->G）进行仿真，最终得到三相电数据并导入MATLAB中进行行波的获取和处理。 

一、仿真环境

PSCAD v46

MATLAB2021b

二、PSCAD仿真搭建

新建case，进行原理图的绘制，总原理图如下：

一、交流输电线

要创建一个交流输电线模型，需要右键单击画布空白位置，选择“Component Wizard”，之后在界面右下角选择“Transmission Line”，输入模型名，选择“finish”后将模型创建在画布上。

双击输电线模型，进入模型定义，进入定义之后把界面上能删的删一删，右键空白位置选择Bergeron模型：

随后选择输电线类型：

最终会显示如下：

 双击“Manual Entry of Y,Z”，配置参数，选择“Manual YZ Data Configuration”，将“Data Entry Method”改为“R,Xl,Xc（ohm/m）”此时“R,Xl,Xc Data Entry [ohm/m]”一栏中的前三个参数就是输电线的正向电阻、正向感抗和正向容抗，注意该组数据，之后在获取行波时有用。

设置完毕后返回画布，右键单击输电线模型选择“Edit Parameters”，设置输电线的长度等参数，此处长度设置为100km，稳态频率设置为50Hz：

 二、故障仿真模块

在主页上选择三相故障模块，至入画布：

其中双击“Timed Fault Logic”可以设定故障时间，双击左边接地部分可以设置故障类型，此处图省事可以直接设置为内部控制，并且设置“Fault Type”为双相接地型故障（AB->G）：

三、测量部分

绘制电流电压表，配置电流电压表，设定输出参数变量名：

主页中在“MISCELLANEOUS”中寻找测量显示模块，将该模块复制到画布上： 

在“Components”一栏中选择“Data Label”，将之前设定给电流电压表的变量名赋值给label，之后将三相电输出的图像用“Output Channel”绘制好：

 随后右键“Output Channel”，选择图像输出：

摁住Ctrl，鼠标左键拖动其他参数进入同一图像中，配置波形颜色和显示量程。

同理，配置好电流输出图像。

配置好一切参数之后，启动仿真。

四、仿真

PSCAD，启动！

仿真结果如图，故障时间是在0.4秒，可以看到电压变化较小，电流变化很大：

右键图像，复制仿真结果到剪贴板：

之后找个txt把电流和电压结果复制下来，准备导入到MATLAB中。

三、MATLAB仿真

将仿真结果导入到MATLAB中，仿真代码如下，注意一下上文中输电线的正序容抗和正序感抗，计算输电线的正序电感和正序电容，公式如下：

其中Xc和Xl为输电线的正序容抗和正序感抗，C和L是正序电感和正序电容，注意一下代码里面的变量名，我忘记改了。

ts=50e-6;%时间间隔为50us
Xl = 1.616e-03;%正序电感
Xc = 7.392e-09;%正序电容


%故障时间为0.4s，采样间隔为50us，故故障点在2000点处
Vm=Vabc(2001:4000,:);
Im=Iabc(2001:4000,:);

Vn=Vabc(1:2000 ,:);
In=Iabc(1:2000 ,:);

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(Xl/Xc);
uf=0.5*(Vx(1,:)+Ix(1,:)*Z);
ur=0.5*(Vx(1,:)-Ix(1,:)*Z);

plot(uf(1:100));
hold on
plot(ur(1:100))

最终得出输出行波如图所示：