从母线导纳矩阵的对角和非对角元素反复更新各母线电压:
$$V_i^{(k+1)}= \frac{1}{Y_{ii}}\left[\frac{P_i - jQ_i}{(V_i^{(k)})^*}- \sum_{j \neq i}Y_{ij}V_j\right]$$收敛判定:$\max_i |V_i^{(k+1)}- V_i^{(k)}| \lt \varepsilon$
线路潮流:$S_{ij}= V_i \cdot I_{ij}^* = V_i(V_i - V_j)^* y_{ij}^*$
在3母线模型上应用高斯-赛德尔迭代法,可同时得到母线电压 V、相角 δ、线路潮流 P、Q、系统损耗(pu单位)。可实时观察收敛过程,是电力系统分析的理想入门工具。
从母线导纳矩阵的对角和非对角元素反复更新各母线电压:
$$V_i^{(k+1)}= \frac{1}{Y_{ii}}\left[\frac{P_i - jQ_i}{(V_i^{(k)})^*}- \sum_{j \neq i}Y_{ij}V_j\right]$$收敛判定:$\max_i |V_i^{(k+1)}- V_i^{(k)}| \lt \varepsilon$
线路潮流:$S_{ij}= V_i \cdot I_{ij}^* = V_i(V_i - V_j)^* y_{ij}^*$
系统规划与设计:新增发电厂(太阳能、风能等)或大型工厂接入既有送电网时,需要事先模拟以确保不会产生电压下降或设备过负荷。高斯-赛德尔法常在初期评估阶段使用。
系统运行:随着日常用电需求的变化,需要确定各发电厂的运行功率(经济负荷分配),或停电恢复时的系统重建步骤。潮流计算结果是这些决策的基础数据。
智能电网分析:评估分布在需求端的多个电源(DER)和储能的控制策略。与传统系统不同,电力流向变为双向,需要更高级的潮流分析。
CAE工具集成:商用电力系统分析软件如PSS®E和DIgSILENT PowerFactory的内部算法基础。大规模系统采用计算速度更快的牛顿-拉夫逊法,但基本原理相同。
初次使用高斯-赛德尔法时容易陷入的几个误区。首先是「收敛判定值ε越小越好」的误解。虽然ε设为1e-8可以提高精度,但计算时间会急剧增加,甚至由于计算机舍入误差的影响反而可能无法收敛。实际工程中,电压精度设为0.1%(ε=1e-3~1e-4)通常就足够了。在这个模拟器中改变ε值并观察迭代次数如何变化,你会有直观理解。
其次是初始值的选择。教科书中说「用额定电压1.0∠0°初始化」,但在实际系统中,承载大负荷的PQ母线的电压可能大幅下降。此时从1.0开始会导致收敛时间增加,甚至最坏情况下发散。根据经验,如果PQ母线初值设为稍低的电压(如0.95∠0°),反而能加快收敛。
最后,要理解该方法的局限性。高斯-赛德尔法计算简单、内存占用少,但收敛性很大程度上依赖于系统的「对角占优性」。也就是说,在电阻(R)和反应(X)的比值(R/X)较大的系统(如配电系统)中,更新公式分母中的Y_ii可能很小,导致计算变得不稳定。实际中,这类系统必须改用牛顿-拉夫逊法进行计算。
假设基准容量Sb=100MVA,母线1(发电机侧)固定为1.0pu∠0°,母线2、3各接50MW负荷的3母线系统。线路1→2:R12=0.02Ω、X12=0.2Ω,线路1→3:R13=0.03Ω、X13=0.25Ω,许可误差0.01MW进行收敛计算。结果:母线2电压为0.97pu∠-2.8°,母线3电压为0.96pu∠-4.1°,系统损耗约3.2MW。通常4~6次迭代即可收敛。