电力系统短路故障分析的MATLAB辅助程序设计-短路计算程序.pdfVIP

电力系统短路故障分析的 MATLAB辅助程序设计 电 力 系 统 短 路 故 障 可 分 为 三 相 对 称 短 路 故 障 (three-phase balanced faults) 和不对称短路故障 (unbalanced faults ) 。不对称 短路故障又分为单相接地短路故障 (single line-to-ground fault) 、 两 相 短 路 故 障 (line-to-line fault) 以 及 两 相 接 地 短 路 故 障 (double line-to-ground fault) 。根据故障分析结果可以对继电保 护装置、 自动装置进行整定计算， 我们可以建立算法来形成节点阻抗 矩阵 , 利用节点阻抗矩阵来计算短路故障情况下的节点电压和线路电 流。 一、三相对称短路故障 进 行 三 相 短 路 计 算 需 要 两 个 程 序 ： zbuild / zbuildpi 和 symfault 程序 ，zbuild 、zbuildpi 用来在 MATLAB中形成节点阻抗 矩阵， symfault 用来计算三相对称故障。 Zbus=zbuild(zdata) 这里的参数 zdata 是一个 (e ×4) 阶矩阵，e 是拓扑图的总支路数目。 第一列和第二列为元素两端的节点编号， 第 三列和第四列分别是线路的电阻、 电抗的标幺值。 连接在 0 节点和发 电机节点之间的发电机阻抗可能是次暂态电抗、 暂态电抗或同步电抗， 而且这个矩阵中还包含并联电抗器和负荷阻抗。 Zbus=zbuildpi(linedata,gendata,yload ) 这个函数与潮流计 算程序是相容的，第一个参数 linedat a 与潮流计算程序中的文件是 一致的。第一列和第二列为节点编号； 第三列到第五列分别是线路的 电阻、电抗以及 1/2 线路电纳值， 这三项都为在统一基准容量下的标 幺值；最后一列是变压器分接头位置，对线路来说，必须输入 1；线 路无输入顺序。发电机参数不包含在 Linedata 参数中，而是包含在 第二个参数 gendata 中，gendata 是一个 g ×4 阶矩阵 ,g 是发电机总 数。第一列和第二列为 0 节点、发电机节点编号，第三列和第四列为 发电机的暂态电阻和暂态电抗。 最后一个参数 yload 是可选择的， 这 个矩阵有两列，第一列为节点编号，第二列为复数导纳值， yload 可 以由潮流程序 lfguss ，lfnewton 或者 decouple 自动生成。 Zbuild 和 zbuildpi 两个函数可以通过建立算法求出节点阻抗 矩阵。首先添加所有与参考节点相连的树支，然后添加其余的树支， 最后添加共轭连支。 程序 symfault(zdata,Zbus,V ) 用来计算三相对称故障，程序要 求输入 zdata 和 Zbus 两个矩阵，第三个参数 V 是可选的。如果 V 不 存在，程序将默认故障前所有的节点电压标幺值为 1.0 ，如果变量 V 存在，那么 V 包括节点编号和复数电压值。 电压向量 V 也可以由潮流 计算程序自动生成。当