Document
拖动滑块完成拼图
首页 专利交易 科技果 科技人才 科技服务 国际服务 商标交易 会员权益 IP管家助手 需求市场 关于龙图腾
 /  免费注册
到顶部 到底部
清空 搜索

一种快速分解状态估计变压器量测雅可比矩阵形成方法 

买专利卖专利找龙图腾,真高效! 查专利查商标用IPTOP,全免费!专利年费监控用IP管家,真方便!

摘要:本发明公开了一种快速分解状态估计变压器量测雅可比矩阵形成方法,采用非标准变比k位于首节点i侧且等值阻抗位于末节点j侧的变压器等值电路,当第l个量测量为变压器首节点i流向末节点j的无功量测量时,Q‑U类常数雅可比矩阵Br的第l行第i列元素Br.li在常规方法基础上乘以3‑2k;量测量为变压器末节点j流向首节点i的无功量测量时,Br的第l行第j列元素Br.lj乘以3‑2k。本发明通过变压器无功量测雅可比矩阵元素乘以一个系数,来改进快速分解状态估计的变压器无功功率量测雅可比矩阵的计算方法,可有效改善状态估计的收敛性。实施算例表明本发明可以显著改善收敛性,明显减少电力系统状态估计的迭代次数。

主权项:1.一种快速分解状态估计变压器量测雅可比矩阵形成方法,其特征在于:包括以下步骤:A、读取网络数据和量测数据;B、电压初值化,设置各节点电压幅值为参考节点电压幅值初始给定值U0,电压相角为0;C、形成P-δ类常数雅可比矩阵Ba;对于节点注入有功功率量测量,Ba为仅用支路电抗形成的节点导纳矩阵的虚部,忽略掉对有功功率分布影响很小的变压器非标准变比和输电线路对地电容以及支路电阻;对于支路有功功率量测量,Ba为仅用该支路电抗形成的节点导纳矩阵的虚部,忽略掉对有功功率分布影响很小的变压器非标准变比和输电线路对地电容以及支路电阻;D、形成Q-U类常数雅可比矩阵Br;形成Br的步骤如下:D1、按传统方法形成Br;对于节点注入无功功率量测量,Br取节点导纳矩阵的虚部;对于支路无功功率量测量,Br取该支路的节点导纳矩阵的虚部;如果第p个量测量为节点i的节点电压幅值量测量,则Br第p行第i列元素为-1U0,Br第p行其他元素都为0;D2、设置量测量计数l=1;D3、取第l个量测量的类型T;D4、判断量测量类型T是否为变压器首端流向末端的无功功率量测量类型TQIJ,如果不是TQIJ类型,转到步骤D7;D5、取第l个量测量对应的变压器的首节点号i、非标准变比k;D6、按式1修正矩阵Br的第l行第i列元素Br.li,转到步骤D10;Br.li=3-2kBr.li1式中,k为变压器非标准变比,变压器等值电路中,非标准变比k位于变压器首节点i侧,等值阻抗位于末节点j侧;D7、判断量测量类型T是否为变压器末端流向首端的无功功率量测量类型TQJI,如果不是TQJI类型,转到步骤D10;D8、取第l个量测量对应的变压器的末节点号j、非标准变比k;D9、按式2修正矩阵Br的第l行第j列元素Br.lj;Br.lj=3-2kBr.lj2D10、令l=l+1;D11、判断l是否大于量测量个数m,如果l大于量测量个数m,转到步骤E;否则返回到步骤D3进行下一次迭代;E、形成快速分解状态估计修正方程的信息矩阵A及C并进行三角分解;F、设置迭代计数t=0,并设ΔUmax=10εU;G、计算有功自由向量at,解方程3求Δδt,并求电压相角最大修正量|Δδmax|;AΔδt=at3式中,A为na×na阶有功信息矩阵;a为na维有功自由向量;Δδ为节点电压相角修正量列向量;na为节点电压相角状态变量个数;上标t表示第t次迭代; 式中,Ba为ma×na阶P-δ类常数雅可比矩阵;za为ma维有功类量测量列向量,包括输电线路和变压器的支路有功功率量测量、节点注入有功功率量测量;ha为对应于za的量测函数;为对应于za的ma阶权重对角矩阵;δ为na维节点电压相角列向量;U为nr维节点电压幅值列向量;ma为有功类量测量个数;nr为节点电压幅值状态变量个数;H、判断电压相角最大修正量|Δδmax|是否小于收敛精度εδ;如果小于收敛精度εδ,转到步骤I;否则,按式6修正电压相角,转到步骤J;δt+1=δt+Δδt6I、判断电压幅值最大修正量|ΔUmax|是否小于收敛精度εU;如果小于收敛精度εU,转到步骤N;否则,转到步骤J;J、计算无功自由向量ct,解方程7求ΔUt,并求电压幅值最大修正量|ΔUmax|;CΔUt=ct7式中,C为nr×nr阶无功信息矩阵;c为nr维无功自由向量;ΔU为节点电压幅值修正量列向量;计算公式如下: 式中,Br为mr×nr阶Q-U类常数雅可比矩阵;zr为mr维无功类量测量列向量,包括输电线路和变压器的支路无功功率量测量、节点注入无功功率量测量和节点电压幅值量测量;hr为对应于zr的量测函数;为对应于zr的mr阶权重对角矩阵;mr为无功类量测量个数;K、判断电压幅值最大修正量|ΔUmax|是否小于收敛精度εU;如果小于收敛精度εU,转到步骤L;否则,按式10修正电压幅值,转到步骤M;Ut+1=Ut+ΔUt10L、判断电压相角最大修正量|Δδmax|是否小于收敛精度εδ;如果小于收敛精度εδ,转到步骤N;否则,转到步骤M;M、令t=t+1,返回到步骤G进行下一次迭代;N、输出状态估计结果,结束。

全文数据:

权利要求:

百度查询: 大连海事大学 一种快速分解状态估计变压器量测雅可比矩阵形成方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息,力求客观、公正,但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解,仅供参考使用,不能作为本公司承担任何法律责任的依据或者凭证。