专利类型:发明专利
语 言:中文
申 请 号:CN201310492987.9
申 请 日:20131018
发 明 人:颜伟段磊赵霞余娟杨焕燕王聪
申 请 人:重庆大学
申请人地址:400044 重庆市沙坪坝区沙正街174号重庆大学电气工程学院
公 开 日:20160127
公 开 号:CN103532137B
代 理 人:胡正顺
代理机构:重庆大学专利中心 50201
摘 要:一种三相四线低压配电网状态估计方法,属于配电网的三相状态估计技术领域。本发明方法利用计算机,通过程序,首先输入采集的任一时间断面的低压配抄表电力数据和网络结构及参数信息并初始化,然后计算节点导纳矩阵和雅可比矩阵常数项部分,接着综合量测方程和T端点零注入等式约束,以电压的实虚部为状态变量,计算不平衡量、指数型权函数及雅可比矩阵,最后更新状态变量,并收敛性判断,实现三相四线低压配电网状态估计。本发明能独立计算节点有功与无功功率和抑制坏数据对估计结果的不良影响,具有强抗差性和良好的收敛性及数值稳定性等特点。本发明方法广泛应用于求解三相四线低压配电网的状态估计。
主 权 项:一种三相四线低压配电网状态估计方法,利用计算机,通过程序,实现三相四线低压配电网的状态估计,其特征在于所述方法的具体步骤如下:(1)输入基础数据及初始化1)输入基础数据首先输入任一时间断面的低压配总表及各用户分表的瞬时抄表电力数据和网络结构及参数信息,即输入低压配总表在任一时间断面下的瞬时抄表总有功功率和总无功功率以及三相电压幅值和三相电流幅值,各用户分表在同一时间断面下的瞬时抄表各相有功功率和各相无功功率数据以及三相电压幅值,网络结构及参数信息为网络线路电阻、电抗和电纳参数、线路的额定电压、功率基准;2)参数初始化设置n阶单位矩阵R?1,n为状态估计中实际的量测变量个数;初始化最大迭代次数Tmax为40~60;收敛精度ε为10?3~10?5及残差方差的检测门槛值γ=5,并设置迭代次数time=1;(2)形成节点导纳矩阵和雅可比矩阵常数项部分第(1)步完成后,计算节点导纳矩阵,计算公式为:式中:B1={a,b,c,n},表示端点包含中性点n的电气节点组合,其中a、b、c分别表示abc三相电气节点;x、d、t表示电气节点组合B1中的任一节点;φi为与端点i直接相联但不包括端点i的端点集合;j为端点集合φi中的任一端点;k为低压配电网所有端点中的任一端点;表示端点i中节点d和x之间的并联支路导纳,其中x∈B1且x≠d;表示端点i中节点d与端点j中节点t直接相联的支路导纳;表示端点i中节点d对地电容得到的对地并联导纳;自导纳表示与端点i节点d直接相联的所有支路导纳和对地电容导纳之和,互导纳表示端点i中节点d和t之间支路导纳的相反数与端点i且与i直接相联端点节点d和t之间所有支路导纳之和,其中t≠d,互导纳则表示端点i和k中节点d和t之间支路导纳的相反数,其中k≠i;得到节点导纳矩阵后,计算雅可比矩阵常数项子矩阵,计算公式为:Cij(1:8,1:8)=Hij(1:8,1:8)???????????(3)式中:i和j为端点编号;Hij和Cij分别为雅可比矩阵H和C的常数项子矩阵;[0]表示全零矩阵,下标3×8表示该矩阵为3行8列;d、t及B1的意义同公式(1);(3)计算不平衡量和指数型权函数及雅可比矩阵第(2)步完成后,计算三相四线低压配电网中量测量的不平衡量和指数型权函数及雅可比矩阵,具体步骤如下:1)计算不平衡量基于公式(4)计算三相四线低压配电网总表端点s的电流不平衡量和有功功率不平衡量及无功功率不平衡量,计算公式为:式中:BP={a,b,c},表示端点不包含中性点n的电气节点组合,其中a、b、c分别表示abc三相电气节点;d表示电气节点组合BP中的任一节点;B1、t的意义同公式(1);s表示三相四线低压配电网总表的端点;为包括总表端点s且与其直接相联端点的集合;k1为端点集合中的任一端点;Δzs为基于公式(4)计算端点s的电流不平衡量和功率不平衡量;和分别为低压配电网总表s相对中性点n的三相电流幅值、三相总注入有功和三相总注入无功的瞬时抄表量测量,其中d∈BP;和分别为端点k1中节点t的电压的实部和虚部;和分别为端点s中节点d的电压的实部和虚部;和分别为端点s的中性点n的电压的实部和虚部;和分别为节点导纳矩阵元素的实部和虚部;基于公式(5)计算三相四线低压配电网用户分表端点i中a、b、c三相节点注入有功功率及注入无功功率的不平衡量,计算公式为:式中:和分别为低压配电网用户分表i中节点d相对中性点n的注入有功和注入无功的瞬时抄表量测量;i表示三相四线低压配电网中任意某个用户分表端点;和分别为端点i中节点d的电压的实部和虚部,其中d∈Bp;和分别为端点i的中性点n电压的实部和虚部;和分别为基于公式(5)计算端点i的d相节点注入有功功率不平衡量和无功功率不平衡量;和分别为节点导纳矩阵元素的实部和虚部;为包括端点i且与其直接相联端点的集合;k2为端点集合中的任一端点;和分别为端点k2中节点t的电压的实部和虚部;BP、d、B1及t的意义同公式(4);基于公式(6)计算三相四线低压配电网用户分表端点i的中性点注入有功功率及注入无功功率的不平衡量,基于公式(7)计算三相四线低压配电网用户分表端点i中a、b、c三相节点电压的不平衡量,计算公式分别为:式中:和分别为基于公式(6)计算低压配电网用户分表端点i的中性点注入有功功率不平衡量和注入无功功率不平衡量;为低压配电网分表i相对中性点n的电压幅值的瞬时抄表量测量;和分别为节点导纳矩阵元素的实部和虚部;为用户分表端点i计算得到的电压不平衡量;k2、BP、d、B1及t的意义同公式(5);基于公式(8)计算零注入电流端点m中a、b、c三相及中性点注入电流的不平衡量Δcm,其计算公式为:式中:和分别为节点导纳矩阵元素的实部和虚部,m∈{1,2,…i,…},k3∈{1,2,…k,…};为包括T端点m且与其直接相联端点的集合;k3为端点集合中的任一端点;和分别为端点k3中节点t的电压的实部和虚部;B1、t、d的意义同公式(1);2)计算指数型权函数基于公式(9)?公式(10)计算指数型权函数W的对角阵元素,计算公式为:式中:wi*为量测i的权函数;为量测i的固定权重;δ为概率密度函数的标准差,其初值由δ2=rTPr/(m?n)计算得到,r为残差向量,rT为残差向量的转置,P为先验矩阵,矩阵P中第i个对角元素为Pi,m为量测量的个数,n为状态变量的个数;rNi为量测i的标准化残差,计算公式为式中:残差方差的检测门槛值γ通过试探法或经验值选取;为量测i对应的残差,为量测i的瞬时抄表电力量测值;hi(x)为量测i对应的非线性量测函数的计算值;R为测量误差方差的对角阵,对角线元素为测量误差的方差;Δ=Ι?H(HTWH)?1HTW;Ι为n阶单位阵;W为指数型权函数对角阵,其对角元素等于权函数,即Wii=wi*;为雅可比矩阵,h(x)为非线性的量测函数,HT为雅克比矩阵的转置;(HTWH)?1表示矩阵HTWH相乘后的逆矩阵;3)形成雅可比矩阵基于公式(11)?公式(15)形成雅可比矩阵H和C,根据公式(16)计算A和B;对于平衡节点,基于公式(11)和(12)形成其雅克比矩阵元素,计算公式为:式中:雅可比子矩阵Hss=[Hss1(1:3,1:3);Hss2(1:3,1:3);Hss3(1,1:3);Hss4(1,1:3)];和分别为节点导纳矩阵元素和的实部,和分别为节点导纳矩阵元素和的虚部,s∈{1,2,…k,…},f∈{1,2,…k,…};和分别为端点f中节点t电压的实部和虚部,f为端点集合中的任一端点;和分别为端点s中的节点p电压的实部和虚部;BP、B1、d、t、s的意义同公式(4);p表示电气节点组合BP中的任一节点;为平衡端点s中节点p电压虚部对平衡端点s中节点p电压实部的偏导数,且式中:雅可比子矩阵Huq=[Huq1(1:3,1:8);Huq2(1:3,1:8);Huq3(1,1:8);Huq4(1,1:8)];u为平衡端点;q为非平衡端点;和分别为节点导纳矩阵元素的实部和虚部;和分别为端点q中节点t电压的实部和虚部;和分别为端点u中的节点d电压的实部和虚部;和分别为端点u的中性点n电压的实部和虚部;BP、d及t的意义同公式(4);对于非平衡节点,基于公式(13)~(15)形成其雅克比矩阵元素,计算公式为:式中:和分别为节点导纳矩阵元素的实部和虚部;BP、B1、p的意义同公式(11);d表示电气节点组合B1中的任一节点;式中:ΔHqq为雅可比子矩阵Hqq的修正项子阵,q为非平衡端点;BP、d的意义同公式(5);和分别为端点q中节点d的电压的实部和虚部,其中d∈Bp;和分别为端点q的中性点n电压的实部和虚部;和均为修正项子阵ΔHqq的矩阵元素,其计算公式分别为:式中:Bp、d、及的意义同公式(14);和分别为低压配电网端点q中节点d相对中性点n的注入有功和注入无功的瞬时抄表量测量;W、H和HT的意义同公式(10),表示迭代值为x0时的不平衡量;(4)状态变量更新和收敛性判断1)状态变量更新第(3)步完成后,根据公式(17)计算状态变量的修正量Δx(time),然后更新状态变量,得到状态变量新值,即:x(time+1)=x(time)+Δx(time),time=time+1,计算公式为Δx=A?1B?A?1CT((CA?1CT)?1(CA?1B?Δc))???????(17)式中,time为计算迭代次数;A?1为矩阵A的逆矩阵,矩阵A为基于公式(16)中矩阵HTWH相乘后计算得到的矩阵,为中间量;矩阵B为基于公式(16)中矩阵HTWΔz相乘后计算得到的矩阵,为中间量;Δc=?c(x0),c(x0)表示迭代值为x0时基于公式(8)计算的零注入电流端点m中a、b、c三相及中性点注入电流的不平衡量;CT为矩阵C的转置,矩阵C为基于公式(3)计算得到的雅克比矩阵;2)收敛性判断当状态变量的修正量Δx(time)满足max(|Δx(time)|)<ε,则结束迭代计算,输出结果;当max(|Δx(time)|)≥ε且迭代次数time≥Tmax,则停止迭代,输出“不收敛!”;当max(|Δx(time)|)≥ε且迭代次数time<Tmax,转入进行如下步骤:若(δ(time))2≥0.01,令(δ(time+1))2=0.001(δ(time))2;若(δ(time))2<0.01,则令(δ(time+1))2=(δ(time))2;并使迭代次数time增加1,返回第(3)步,进行重新迭代计算,其中,δ(time)为概率密度函数的标准差δ的第time次迭代计算值。
关 键 词:
法律状态:公开
IPC专利分类号:H02J3/00(2006.01)I