浏览量:0

一种不接地配电网的三相潮流计算方法

专利类型:发明专利 

语 言:中文 

申 请 号:CN201610477990.7 

申 请 日:20160627 

发 明 人:赵霞罗兰叶晓斌颜伟余娟 

申 请 人:重庆大学 

申请人地址:400044 重庆市沙坪坝区沙正街174号 

公 开 日:20161012 

公 开 号:CN106026082A 

代 理 人:王翔 

代理机构:重庆大学专利中心 50201 

摘  要:本发明公开一种不接地配电网的三相潮流计算方法,获取网络的结构数据、参数信息、平衡端点及PV端点的控制信息和端点负荷信息并初始化,然后将平衡端点的线电压转换为基于等值中性点的相电压,接着计算网络基于平衡端点等值中性点的节点导纳矩阵和雅克比矩阵常数项,进而综合考虑网络非平衡端点的节点电流不平衡方程、零序电流为零的约束方程、PV端点的总有功约束方程和线电压幅值约束方程,以各端点电压的实虚部、星形接线负荷端点中性点电压的实虚部、PV端点注入电流的实虚部作为状态变量,采用基于注入电流形式的牛顿拉夫逊法算法,计算相应的不平衡量、雅克比矩阵,最后更新状态变量,进行收敛性判断,实现不接地配电网的三相潮流计算。 

主 权 项:一种不接地配电网的三相潮流计算方法,其特征在于,包括以下步骤:(1)获取基础数据及初始化1)获取基础数据获取网络线路电阻、电抗参数,负荷的功率值和接线方式,平衡端点线电压,PV端点的控制目标,包括总有功和三个线电压幅值,端点类型,线路的额定电压和功率基准值。其中平衡端点的三个线电压非独立,只需给出线电压另一个线电压处于同一地理位置的不同相的多个电气节点称为端点,每一个端点包含a,b,c三相节点,若星形连接,则还包括中性点n。2)参数初始化设置端点各节点电压的初始值为:幅值为1.0pu,相角相差120度,且相角为0;各星形接线负荷端点中性点电压的初始值为:幅值为0.5pu,相角为0度;设置PV端点注入电流的初始值:幅值1.0pu,相角为0度;最大迭代次数Tmax为20;收敛精度ε为10?6;迭代次数time=1。(2)形成平衡端点等值相电压寻找平衡端点的所述等值中性点,并将该等值中性点作为全网的电压参考点,从而将平衡端点的已知线电压转换为基于等值中性点的等值相电压。等值相电压与线电压的转换关系如下:U·refaU·refbU·refc=A-1TAU·refabU·refbcU·refca---(1)其中:式中ref表示平衡端点的序号;为已知的平衡端点线电压;为平衡端点等值相电压。(3)形成节点导纳矩阵和雅克比矩阵常数项1)形成节点导纳矩阵设φn表示网络所有端点所在的集合,总端点数为N;平衡端点仅一个,且序号为1;φm表示非平衡端点所在集合,非平衡端点个数为m;φy为接有恒功率星形负荷的端点所在集合,接有恒功率星形负荷的端点总数为Ny,;φl为恒功率三角形负荷所在端点集合,总数为Nl;φpv表示PV端点所在集合,总数为Npv。第(1)步完成后,根据输入的各线路阻抗信息得到支路阻抗:Zl_uv=zaazabzaczbazbbzbczcazcbzcc---(3)式中Zl_uv表示端点u和端点v(u,v∈φn,)。间线路的支路阻抗矩阵;zxy,zxx分别表示相间阻抗和各相的自阻抗,其中x,y∈Bp,x≠y,Bp={a,b,c},a,b,c分别表示各端点的三相电气节点。上述步骤完成后可求取对应的导纳矩阵Yl_uv,公式如下:Yl_uv=Zl_uv?1????????(4)式中Yl_uv表示支路导纳矩阵。然后计算网络的节点导纳矩阵,计算公式如下:Yuv=ΣkφiYl_uk,v=u-Yl_uv,vu---(5)式中:φu为与端点u直接相联但不包括端点u的端点集合;k为端点集合φu中的任一端点;Yl_uk表示端点u与端点k直接相联线路的导纳矩阵,为3x3的复数矩阵;Yl_uv表示端点u与端点v直接相联线路的导纳矩阵,为3x3的复数矩阵;自导纳矩阵Yuu表示与端点u直接相联的所有支路导纳矩阵之和,为3x3的复数矩阵;互导纳矩阵Yuv(v≠u)则为端点u和端点v之间支路导纳矩阵Yl_uv的相反数,为3x3的复数矩阵;2)计算雅克比矩阵常数项得到网络的节点导纳矩阵YUV后,可计算雅克比矩阵的常数项子矩阵,计算公式为:Hc(3i-5:3i-3,3j-5:3j-3)=[-Bij]3×3Hc(3i-5:3i-3,3j+3m-5:3j+3m-3)=[-Gij]3×3Hc(3i+3m-5:3i+3m-3,3j-5:3j-3)=[-Gij]3×3Hc(3i+3m-5:3i+3m-3,3j+3m-5:3j+3m-3)=[Bij]3×3---(6)式中:i和j为非平衡端点编号(u,v∈φm,);Hc表示雅克比矩阵的常数项,为6mx6m的矩阵;Gij和Bij分别表示导纳矩阵元素Yij的实部和虚部;下标3x3表示该矩阵为3行3列矩阵.(4)计算不平衡变量步骤(3)完成后,结合端点注入电流、端点电压、负荷功率和PV端点的控制变量可计算不平衡变量,具体包括:非平衡端点的注入电流不平衡量实虚部、恒功率星形负荷端点零序电流不平衡量的实虚部、PV端点控制变量的不平衡量,然后形成不平衡变量矩阵,求解过程如下:1)计算非平衡端点的注入电流不平衡量完成上述步骤后,再结合步骤(3)计算出的网络导纳矩阵,计算非平衡端点的注入电流不平衡量实虚部,计算公式如下:dI·id=I·id-ΣjφiΣtBpYijdtU·jt---(7)式中:i表示不接地配电网的端点;为包括端点i且与其直接相联端点的集合;j为端点集合中的任一端点;为i端点d相的注入电流,d∈Bp,Bp的含义与式(3);表示i端点d相的电流不平衡量;表示i端点d相与j端点t相间的导纳;分别表示j端点节点t的电压,t为a,b,c中的任意相。其中,为电源注入电流与负荷注入电流之和,即表示为:I·id=I·Gid+I·Lid---(8)若端点i为PV端点,则电源注入电流为待求变量;若端点i无电源连接,则为0.若负荷为星形恒功率负荷,则负荷注入电流为:I·sa=-(S~sanU·sa-U·sn)*I·sb=-(S~ibnU·sb-U·sn)*I·sc=-(S~scnU·sc-U·sn)*---(9)式中s表示恒功率星形负荷的端点序号;分别为给定的a,b,c相与中性点间负荷的视在功率;为端点s中性点的电压;表示端点s节点d的电压d∈Bp,;表示s端点d相的注入电流,d∈Bp,;Bp的含义与式(3)相同;“*”为共轭运算。若负荷为三角形恒功率负荷,则负荷注入电流为:I·la=(S~lcaU·lc-U·la)*-(S~labU·la-U·lb)*I·lb=(S~labU·la-U·lb)*-(S~lbcU·lb-U·lc)*I·lc=(S~lbcU·ib-U·ic)*-(S~lcaU·lc-U·la)*---(10)式中l表示恒功率三角形负荷的端点序号;分别为给定的线间负荷;表示端点l节点d的电压d∈Bp,;表示l端点d相的注入电流,d∈Bp,;Bp的含义与式(3)相同;“*”为共轭运算。若端点不接负荷,则对应的负荷注入电流为0.2)计算恒功率星形负荷端点的零序电流不平衡量恒功率星形负荷端点的零序电流不平衡量计算如下:dI·sn=0-(I·sa+I·sb+I·sc)=(S~sanU·sa-U·sn)*+(S~sbnU·sb-U·sn)*+(S~scnU·sc-U·sn)*---(11)式中表示端点s的零序电流不平衡量;其余变量含义同式(9)。3)计算PV端点控制变量不平衡量PV端点的控制变量包括总注入有功和三线电压幅值,其对应的不平衡量计算公式如下:dPkΣ=PkΣ-ΣdBp(ekdgkd+fkdhkd)dUm_kab=(Um_kab)2-[(eka-ekb)2+(fka-fkb)2]dUm_kbc=(Um_kbc)2-[(ekb-ekc)2+(fkb-fkc)2]dUm_kca=(Um_kca)2-[(ekc-eka)2+(fkc-fka)2]---(12)式中k表示PV端点序号;分别表示给定的总注入有功和端点线电压幅值;分别表示总注入有功不平衡量和线电压幅值不平衡量;为电压的实虚部;为端点k注入电流的实虚部,其中d∈Bp,Bp的含义与式(3)相同。4)不平衡变量矩阵上述不平衡变量计算完成后就可形成不平衡变量矩阵ΔF,公式如下:ΔF=[dhdg||dhndg|n|dpΣdumabdumbcdumca]T(13)式中dh和dg分别为非平衡端点注入电流不平衡量dI的虚部和实部,长度均为3m;dhn和dgn分别为星形恒功率负荷端点零序电流不平衡量的虚部和实部,长度均为Ny;dPΣ分别为PV端点的功率不平衡量、ab线电压幅值不平衡量、bc线电压幅值不平衡量和ca线电压幅值不平衡量,长度分别为Npv。各变量按照端点序号排列。(5)计算雅克比矩阵在雅可比矩阵常数项(式(6))的基础上,计及雅可比矩阵的可变项可形成完整的雅可比矩阵,具体求解公式如下:1)端点注入电流不平衡方程对负荷端点电压求导恒功率星形负荷和三角形负荷所在端点的注入电流不平衡方程可对其节点电压求导。恒功率星形负荷端点的电流不平衡方程对其节点电压的实虚部求导公式如下:Ises=(S~san(U·sa-U·sn)2)*000(S~sbn(U·sb-U·sn)2)*000(S~scn(U·sc-U·sn)2)*Isfs=-jIses---(14)H(3s-5:3s-3,3s-5:3s-3)=imag(Ises)H(3s-5+3m:3s-3+3m,3s-5:3s-3)=real(Ises)H(3s-5:3s-3,3s-5+3m:3s-3+3m)=imag(Isfs)H(3s-5+3m:3s-3+3m,3s-5+3m:3s-3+3m)=real(Isfs)---(15)式中“real(·)”、“imag(·)”分别表示取实部和虚部。恒功率三角形负荷端点的电流不平衡方程对节点电压的实虚部求导公式如下:Ilel=I·lcaela-I·labela-I·labelbI·lcaelcI·labelaI·labelb-I·lbcelb-I·lbcelc-I·lcaelaI·lbcelbI·lbcelc-I·lcaelcIlfl=-jIlel---(16)其中I·lcaela=(S~lca(U·lc-U·la)2)*I·lcaelc=-(S~lca(U·lc-U·la)2)*I·labela=-(S~lab(U·la-U·lb)2)*I·l

关 键 词: 

法律状态:公开 

IPC专利分类号:H02J3/00(2006.01)I;G06F17/50(2006.01)I