本发明专利技术公开了一种直角坐标牛顿法潮流计算方法，包括以下步骤：原始数据输入和电压初始化；形成节点导纳矩阵；计算功率及电压偏差；形成雅可比矩阵J；解修正方程及修正电压实部e、虚部f；节点及支路数据输出。本发明专利技术通过对直角坐标牛顿法潮流计算的雅可比矩阵部分元素(i＝j时)的计算公式进行修改和改进，解决了直角坐标牛顿法潮流计算在分析含有小阻抗支路系统时的收敛性问题。采用常规直角坐标牛顿法潮流计算不收敛时，本发明专利技术能够可靠收敛。由于本发明专利技术不仅能有效解决常规直角坐标牛顿法潮流计算分析含有小阻抗支路系统的收敛性问题，同时也能对正常系统进行潮流计算，因此没有不良影响。
【技术实现步骤摘要】
【技术保护点】
一种直角坐标牛顿法潮流计算方法，包括以下步骤：A、原始数据输入和电压初始化电压初始化采用平启动，即PV节点和平衡节点的电压实部取给定值，PQ节点的电压实部取1.0；所有电压的虚部都取0.0；这里单位采用标幺值；B、形成节点导纳矩阵小阻抗支路导纳为yij=jbij=-j1xij---(2)]]>设节点i和节点j原来的自电导与自电纳分别为Gi0、Bi0、Gj0、Bj0，在它们之间增加一条小阻抗后的自导纳和互导纳分别为：Yii=Gi0+j(Bi0+bijk2)=Gi0+j(Bi0-1k2xij)---(3)]]>Yjj=Gj0+j(Bj0+bij)=Gj0+j(Bj0-1xij)---(4)]]>Yij=-jbijk=j1kxij---(5)]]>C、计算功率及电压偏差功率及电压偏差计算公式为：ΔPi=Pis-Pi=Pis-eiai-fibiΔQi=Qis-Qi=Qis-fiai+eibiΔVi2=Vis2-(ei2+fi2)---(6)]]>式中，Pis、Qis分别为节点i给定的注入有功功率和无功功率；Vis为节点i给定的电压幅值；ai、bi分别为节点i的计算注入电流相量的实部和虚部，为ai=Σj=1n(Gijej-Bijfj)bi=Σj=1n(Gijfj+Bijej)---(7)]]>D、形成雅可比矩阵JE、解修正方程及修正电压实部e、虚部f修正方程为：ΔPΔQΔV2=JΔeΔf=∂ΔP∂eT∂ΔP∂fT∂ΔQ∂eT∂ΔQ∂fT∂ΔV2∂eT∂ΔV2∂fTΔeΔf---(20)]]>式中，J为雅可比矩阵；电压修正公式为：ei(k+1)=ei(k)-Δei(k)fi(k+1)=fi(k)-Δfi(k)---(21)]]>式中，上标(k)表示第k次迭代；F、节点及支路数据输出；其特征在于：所述的形成雅可比矩阵J的方法包括以下步骤：F1、当i≠j时，雅可比矩阵J的元素计算公式如下：∂ΔPi∂ej=-Gijei-Bijfi---(8)]]>∂ΔPi∂fj=Bijei-Gijfi---(9)]]>∂ΔQi∂ej=Bijei-Gijfi---(10)]]>∂ΔQi∂fj=Gijei+Bijfi---(11)]]>∂ΔVi2∂ej=0]]>∂ΔVi2∂fj=0---(13)]]>F2、当i＝j时，雅可比矩阵J的元素计算公式如下：∂ΔPi∂ei=-aiS-Giiei-Biifi---(22)]]>∂ΔPi∂fi=-biS+Biiei-Giifi---(23)]]>∂ΔQi∂ei=biS+Biiei-Giifi---(24)]]>∂ΔQi∂fi=-aiS+Giiei+Biifi---(25)]]>∂ΔVi2∂ei=-2ei---(18)]]>∂ΔVi2∂fi=-2fi---(19)]]>式中，aiS、biS分别为节点i给定的注入电流相量的实部和虚部，由式(6)求得；潮流计算收敛时，式(6)中ΔPi、ΔQi都趋近于0，因此，由给定值PiS和QiS求ai和bi，记为aiS和biSaiS=eiPiS+fiQiSei2+fi2biS=fiPiS-eiQiSei2+fi2---(26).]]>...
【技术特征摘要】
【专利技术属性】
技术研发人员：姚玉斌，吴志良，王丹，
申请(专利权)人：大连海事大学，
类型：发明
国别省市：辽宁;21
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。