潮流程序的开发.docVIP

课题组前期培训阶段二 — 潮流程序的开发 数据结构的设计 1.1输入数据 输入的数据为线路，双绕组变压器，三绕组变压器的铭牌值及母线的参数。具体如下： 线路参数存储到line矩阵中，格式如下： 首端节点序号 末端节点序号 单位长度电阻 单位长度电抗 单位长度电纳 线路长度 双绕组变压器参数存储到trans矩阵中，格式如下： 高压侧节点序号 低压侧节点序号 容量 高压侧额定电压 低压侧额定电压 短路损耗 短路电压百分数 空载损耗 空载电流百分数 三绕组变压器参数存储到tw_trans矩阵中，格式如下： 高压侧节点序号 中压侧节点序号 低压侧节点序号 高压侧容量 中压侧容量 低压侧容量 高压侧额定电压 中压侧额定电压 低压侧额定电压 短路损耗P12 短路损耗P23 短路损耗P31 短路电压百分数 V12 短路电压百分数 V23 短路电压百分数 V31 空载电流百分数I0 空载损耗P0 母线参数存储到bus矩阵中，格式如下： 节点序号 基准电压 节点类型 注入有功 注入无功 电压幅值 电压相角 节点类型中：1代表PQ节点 2代表PV节点 3代表平衡节点 1.2输出数据 输出数据存储到poweFlow矩阵中格式如下 节点序号 注入有功 注入无功 电压幅值 电压相角 程序流程框图 程序功能模块介绍 3.1节点导纳矩阵的形成 函数名称：daonaForm function daona=daonaForm(Sb,line,bus,trans,tw_trans) 输入参数：Sb(基准功率),line(线路参数),bus(母线参数),trans(双绕组变压器参数)，tw_trans(三绕组变压器参数) 输出参数:daona(节点导纳矩阵) 程序步骤： 1、根据线路参数通过公式计算出线路的的阻抗（取倒数用导纳表示），对地导纳。把上述数据按如下格式存储到line_admi矩阵中。如果为对地支路，首末端节点序号取同一个。 支路首端节点序号 支路末端节点序号 支路导纳数值 2、根据双绕组铭牌参数计算双绕组变压器的阻抗支路和励磁支路，均用导纳表示。把以上数据按如下格式存储到trans_admi矩阵中。如果为对地支路，首末端节点序号取相同值。 支路首端节点序号 支路末端节点序号 支路导纳数值 3、根据公式计算三绕组变压器的三侧的阻抗及励磁支路参数，注意三绕组变压器会增加一个节点。把以上数据按如下格式存储到twtrans_admi矩阵中。如果为对地支路，首末端节点序号取同一个。 支路首端节点序号 支路末端节点序号 支路导纳数值 4、根据line_admi、trans_admi、twtrans_admibusRank function [busRanking,busGain,numPQ,numPV]=busRank(bus,tw_trans) 该函数的功能是对bus矩阵进行处理。包括一下两方面：1、三绕组变压器的模型会增加PQ节点,这些节点要在不平衡功率方程及雅克比矩阵和求解修正量中反应出来。2、对节点矩阵进行重新排列，按照PQ节点、PV节点、平衡节点的顺序来排序，以方便后面形成雅克比矩阵。 输入参数：bus原始的母线参数，tw_trans三绕组变压器参数矩阵 输出参数：busRanking 经过按PQ节点、PV节点、平衡节点的顺序排列过的母线参数；busGain 仍然按照节点顺序排列的 （但考虑三绕组增加节点后）母线参数；numPQ、numPV 此函数的功能是形成雅克比矩阵。 输入参数：busRanking，busGain，daona 输出参数：jac 程序的思路： 将雅克比矩阵分成四大块分别求解。四大块形式分别为： 然后根据H、N、J、L的公式分别求解，并将元素放到对应的位置。求解中要分和两种情况。 3.5根据修正量更新节点矩阵 busGain busRanking 函数名称为：updateBus function[busGain,busRanking]=updateBus(reviseInitial,busGain,busRanking,numPQ,numPV) 该函数的功能是将修正方程的解更新到母线参数矩阵。 输入参数;reviseInitial,busGain,busRanking,numPQ,numPV。 输出参数：busGain，busRanking ，即更新后的母线参数矩阵。 程序思路： 根据雅克比矩阵和不平衡功率相除后得到的修正量来修正节点矩阵中的电压和相角。首先是相角的修正，直接在busRanking中的相角的基础上加上修正量即可。电压的修正需要注意的是：修正量为 所以要用原始电压乘以修正量再加上原始电压。最后根据修正