文档介绍：0 引言
潮流是配电网络分析的基础,用于电网调度、运行分析、操作模拟和设计规划,同时也是电压优化和网络接线变化所要参考的内容。潮流计算通过数值仿真的方法把电力系统的详细运行情况呈现给工作人员,从而便于研究系统在给定条件下的稳态运行特点。随着市场经济的发展,经济利益是企业十分看重的,而线损却是现阶段阻碍企业提高效益的一大因素。及时、准确的潮流计算结果,可以给出配电网的潮流分布、理论线损及其在网络中的分布,从而为配电网的安全经济运行提供参考。从数学的角度来看,牛顿-拉夫逊法能有效进行非线性代数方程组的计算且具有二次收敛的特点,具有收敛快、精度高的特点,在输电网中得到广泛应用。随着现代计算机技术的发展,利用编程和相关软件,可以更好、更快地实现配电网功能,本文就是结合牛顿-拉夫逊法的基本原理,利用C++程序进行潮流计算,计算结果表明该方法具有良好的收敛性、可靠性及正确性。
1 牛顿-拉夫逊法基本介绍
潮流方程
对于N个节点的电力网络(地作为参考节点不包括在内),如果网络结构和元件参数已知,则网络方程可表示为:
(1-1)
式中,为N*N阶节点导纳矩阵;为N*1维节点电压列向量;为N*1维节点注入电流列向量。如果不计网络元件的非线性,也不考虑移相变压器,则为对称矩阵。
电力系统计算中,给定的运行变量是节点注入功率,而不是节点注入电流,这两者之间有如下关系:
(1-2)
式中,为节点的注入复功率,是N*1维列矢量;为的共轭;是由节点电压的共轭组成的N*N阶对角线矩阵。
由(1-1)和(1-2),可得:
上式就是潮流方程的复数形式,是N维的非线性复数代数方程组。将其展开,有:
j=1,2,….,N (1-3)
式中, 表示所有和相连的节点,包括。
将节点电压用极坐标表示,即令,代入式(1-3)中则有:
故有:
i=1,2,…,N (1-4)
式(1-4)是用极坐标表示的潮流方程。
而节点功率误差:
(1-5)
(1-6)
式中:,为节点给定的有功功率及无功功率。
牛顿-拉夫逊法基本原理
牛拉法的一般描述
牛拉法是把非线性方程式的求解过程变成反复对相应的线性方程式的求解过程,即非线性问题通过线性化逐步近似,这就是牛拉法的核心。下面以非线性方程式的求解过程来进行说明。
设电力网络的节点功率方程一般形式如下:
(1-7)
式中,为节点注入功率给定值;为对应的物理量和节点电压之间的函数表达式;
为节点电压。写成功率偏差的形式:
(1-8)
应用牛拉法求解如下。在给定的初值处将式(1-8)作一阶泰勒展开:
定义为潮流方程的雅克比矩阵,为在处的值,则有:
用修正就得到的新值。如果用k表示迭代次数,写成一般的表达式,有:
(1-9)
对于潮流收敛的情况,应比更接近于解点。收敛条件为:
由简单迭代法收敛性分析的结论知,越接近解点,牛顿-拉夫逊法收敛越快,它具有二阶收敛速度。:
牛顿-拉夫逊法的几何解释
极坐标的牛顿-拉夫逊法
在极坐标中,有如下的形式:
(1-10)
共2n-r个方程,状态变量为:
共2n-r个待求量。r 个PV节点的电压幅值给定,不需求解。潮流雅克比矩阵的维数是(2n-r)*(2n-r),结构如下:
上式右侧的对电压幅值的偏导数项中的电压幅值的阶数减少了1,为使雅克比矩阵的各部分子矩阵具有一致的形式,在实际计算中,常将该项乘以电压幅值,并选取作为待求的修正量,则雅克比矩阵可写成:
(1-11)
将式(1-10)和(1-11)代入式(1-9)的修正方程即可求得x的修正量,用它修正x直到为止。
将式(1-11)用下式表示:
其中每个字块的计算公式如下:
对角元素: (1-12)
非对角元素: (1-13)
2 牛顿法潮流计算步骤
程序流程图
在了解了牛拉法的原理之后,明确程序编写思路,、。“计算电压幅值和角度”步骤较多,。
电压幅值和角度求解步骤框图
当不符合收敛的条件“amontk>1”时,即认为计算不收敛。具体程序见附录。
计算步骤
下面讨论的是极坐标形式的牛顿法求解过程,大致分为以下几个步骤:
①形成节点导纳矩阵;
②给各节点电压设初值();
③根据式(1-12)、(1-13)生成雅克比矩阵(H、N、M、L);
④将节点电压初值代入式(1-5)、式(1-6),求出修正方程式的常数项向量;
⑤求解修正方程,得到电压幅值和角度;
⑥判断是否收敛,若收敛,计算平衡节点和线路功率;
⑦输出结果