chp3.2复杂电力系统潮流计算机计算11v研讨.pptVIP

* 然后将其代入下式，求出节点p的电压： 在迭代过程中，按上式求得的节点p的电压大小不一定等于设定的节点电压Up0，所有在下一次的迭代中，应以设定的Up0对电压进行修正，但其相角仍保持上式所求得的值，使得 如果所求得PV节点的无功功率越限，则无功功率在限，该 PV节点转化为PQ节点。 * 高斯-塞德尔迭代法计算潮流的步骤： 设定各节点电压的初值，并给定迭代误差判据； 对每一个PQ节点，以前一次迭代的节点电压值代入功率迭代方程式求出新值； 对于PV节点，求出其无功功率，并判断是否越限，如越限则将PV节点转化为PQ节点； 判别各节点电压前后二次迭代值相量差的模是否小于给定误差，如不小于，则回到第2步，继续进行计算，否则转到第5步； 根据功率方程求出平衡节点注入功率； 求支路功率分布和支路功率损耗。 * 雅可比矩阵的特点 雅可比矩阵为一非奇异方阵。传统的，当节点电压以极坐标表示时，该矩阵为2(n-1)-m阶方阵（m为PV节点数）；当节点电压以直角坐标表示时，该矩阵为2(n-1)阶方阵。现在，为了便于编程，一般为经过处理的2n阶。 矩阵元素与节点电压有关，故每次迭代时都要重新计算。 与导纳矩阵具有相似的结构，当Yij=0，Hij、Nij、Jij、Lij均为0，因此也是高度稀疏的矩阵。 具有结构对称性，但数值不对称 注意：当在计算过程中发生PV节点的无功功率越限时，PV节点要转化为PQ节点 * 牛顿-拉夫逊法计算电力系统潮流有关问题 稀疏矩阵表示法 节点导纳矩阵：高度稀疏的N阶复数对称方阵。因此记录矩阵的下三角。 用数组表示 数组1：记录矩阵对角元素的数值； 数组2：记录矩阵非对角元素的数值（按列存储） ； 数组3：记录矩阵非对角元素的行号； 数组4：记录矩阵非对角元素的按行排的位置数； 数组5：记录矩阵非对角元素的按行存储对应按列存储的位置数 * 非对角元素用指针表示，一个指针用结构表示： 行号； 列号； 幅值； 角度； 指针（指向下一个非零元素）。 对角元素用一个一维数组表示。 雅可比矩阵：高度稀疏的2N阶实数方阵，其形式对称但数值不对称。其稀疏程度与节点导纳矩阵相同，可根据节点导纳矩阵形成。 * 高斯消去法 求解牛顿-拉夫逊法潮流计算的修正方程，可以采用矩阵求逆的方法。但是由于潮流计算的雅可比矩阵通常是一个高度稀疏的矩阵，其逆阵则是一个满矩阵，因此用求逆的方法会增加额外的存储单元和计算工作量。而用高斯消去法则可以保持方程组原有的稀疏性，可以大大减少计算所需的内存和时间。 * 节点的优化编号 静态优化法：按静态联结支路数的多少编号。 统计好网络中各节点联结的支路数后，按联结支路数的多少，由少到多，顺序编号。 半动态优化法：按动态联结支路数的多少编号。 先只编一个联结支路数最小的节点号，并立即将其消去；再编消去第一个节点后联结支路数最小的节点号，再立即将其消去……依此类推。 动态优化法：按动态增加支路数的多少编号。 不首先进行节点编号，而是寻找消去后出现的新支路数最少的节点，并为其编号，且立即将其消去；然后再寻找第二个消去后出现的新支路数最少的节点并为其编号，再立即将其消去……依此类推。 * 牛顿-拉夫逊法的收敛特性 牛顿-拉夫逊法具有平方收敛特性，高斯-塞德尔法为一阶收敛特性。 牛顿-拉夫逊法对初值设定很敏感。因此，在实际应用当中，常常在牛顿-拉夫逊法计算潮流以前先用对初值不敏感的高斯-塞德尔法（迭代1-2次）计算电压的初值。 * 4-5 P-Q分解法 P-Q分解法是牛顿-拉夫逊法潮流计算的一种简化方法。 牛顿-拉夫逊法的缺点：牛顿-拉夫逊法的雅可比矩阵在每一次迭代过程中都有变化，需要重新形成和求解，这占据了计算的大部分时间，成为牛顿-拉夫逊法计算速度不能提高的主要原因。 P-Q分解法利用了电力系统的一些特有的运行特性，对牛顿-拉夫逊法做了简化，以改进和提高计算速度。 * 牛顿-拉夫逊法简化形成P-Q分解法的过程 牛顿-拉夫逊法修正方程展开为： 根据电力系统的运行特性进行简化： 考虑到电力系统中有功功率分布主要受节点电压相角的影响，无功功率分布主要受节点电压幅值的影响，所以可以近似的忽略电压幅值变化对有功功率和电压相位变化对无功功率分布的影响，即： * 根据电力系统的正常运行条件还可作下列假设： 电力系统正常运行时线路两端的电压相位角一般变化不大（不超过10~20度）； 电力系统中一般架空线路的电抗远大于电阻； 节点无功功率相应的导纳Q/U*U远小于该节点的自导纳的虚部。 用算式表示如下： * 由以上假设，可得到雅可比矩阵的表