Matpower建模
内部数据的数值全部由标幺值表示,复数的角度由幅度表示。
所有离线的发电机和支路需要在建立用于潮流计算的模型之前被移除,连续潮流计算或最优潮流问题,所有的总线从1开始被连续标号。
从外部到内部的编码和从内部到外部的解码分别由函数ext2int和int2ext完成。
Matpower中所有的模型和等式都被表达为矩阵和向量的格式。
1.数据格式
Matpower case format的细节可以使用命令行 help caseformat
查看
bus,branch,gen的数量分别为nb,nl,ng。
2.支路
所有传输线路、变压器和移相器都建模为一个标准的π型线路模型,该模型由串联阻抗zs=rs+jxs和总充电纳bc组成,并串联了一个理想的移相变压器。该移相器的分接头比例为τ,移相角度为θshift位于支路的始端。如图3-1所示。
参数rs,xs,bc,τ,θshift由支路矩阵中列BR_R,BR_X,BR_B,TAP,SHIFT指定。
列TAP中的值为0,代表该支路是一个没有变压器的传输线,数学上等同于线路上有一个变比为1的变压器。
注入复电流if,it分别位于支路的始端和终端,可由2x2的导纳矩阵Ybr和端口电压vf,vt表示,见式3-1。
将π型线路模型中的串联阻抗zs表示为导纳ys=1\rs+jxs,支路导纳矩阵可被写为(支路导纳矩阵的推导见《电路分析》):
用四个元素表示Ybr矩阵:
则可以形成四个nlx1(nl是支路的数量)的向量Yff,Yft,Ytf,Ytt,每个向量的第i个元素来自对应的Yibr。同时,nlxnb的稀疏联络矩阵Cf,Ct被用于形成系统导纳矩阵。
Cf矩阵中的(i,j)th元素和Ct矩阵中的(i,k)th为1时代表,支路i连接了总线j和k。
3.发电机
发电机被建模为指定总线处的注入复功率,对于发电机i,注入功率为sig=pig+jqig
Sg=Pg+jQg是维度为ngx1的向量,发电机的pig,qig的MVar和MW当量(转换为标幺值之前)由发电机矩阵中的PG(2)和QG(3)指定。
一个稀疏nbxng大小的发电机连接矩阵Cg,当发电机j位于总线i时,矩阵中的(i,j)元素置为一,否则置为0。
则nbx1大小,表示来自发电机的所有总线注入功率可以被表示为:Sg,bus=Cg·Sg
一个有负注入的发电机可以被建模为可分配负载
可分配负载的定义:
参考文献:Data Centers as Dispatchable Loads to Harness Stranded Power
Kibaek Kim, Fan Yang, Member, IEEE, Victor M. Zavala, Member, IEEE, and Andrew A. Chien Fellow, IEEE
4.负载
恒定功率负载被建模为在总线上消耗的特定有功功率和无功功率。
Sd=Pd+jQd,是表示所有总线上的复功率的大小为nbx1的向量。MW和MVar当量(转为标幺值之前)分别由总线矩阵中的列PD(3)和QD(4)指定。
这些字段可以取负值来表示分布式电源(distributed generation)
恒定阻抗和恒定电流负载不能直接实现,但可以将恒定阻抗部分建模为下述并联元件。
可分配负载可被建模为负发电机(negative generator),Sg表示为负数
5.并联元件
电容、电感等并联元件被建模为总线上对地固定阻抗,总线i上的并联元件的导纳:
Ysh=Gsh+jBsh是表示所有总线上的并联导纳的大小为nbx1的向量。
参数gish,bish由总线矩阵(bus matrix)中的GS(5),BS(6)指定,等效 MW(消耗)和 MVAr(注入),标称电压为 1.0 p.u,角度为零。
并联元件还可用于建模恒定阻抗负载
6网络等式
对于具有nb条总线的网络,该模型的所有恒定阻抗元素都包含在复数总线导纳矩阵Ybus中,该矩阵将复数节点电流注入Ibus与节点复电压V联系起来:
同样的,对于有nl条支路的网络,nlxnb的系统支路导纳矩阵Yf和Yt将总线电压与大小为nlx1的If和It的分别位于始端和终端的支路电流向量联系在一起。
[ · ]用于表示利用nx1的向量创建相应的nxn的对角矩阵,并将向量的元素置于矩阵对角线的运算,则系统导纳矩阵可以被下列计算构造
然后将节点母线注入与来自负载和发电机的注入相匹配,以形成交流节点功率平衡方程,表示为复母线电压和发电机注入的复矩阵形式:
https://wenku.baidu.com/view/50be7d31336c1eb91b375d45.html?sxts=1564718002896