本发明专利技术涉及一种基于状态估计残差参数灵敏度的电网参数估计方法，属于电力系统电网参数辨识技术领域。该方法包括：根据实际的电网模型和实时量测数据，推导出状态估计残差对电网每条线路参数的灵敏度，采用线性优化方法，在已有的可疑支路集上修改线路参数，从而提高应用软件的可靠性和精确性。本发明专利技术方法不存在数值稳定性和观测性问题，而且计算简单，容易实现，计算速度快。
【技术实现步骤摘要】
本专利技术涉及，属于电力系统电 网参数辨识
技术介绍
支路参数的误差会导致基于能量管理系统(EMS)的电网在线控制决策和安全稳定预 警结果的不可信。随着电网规模的扩大，基于EMS网络分析的新应用软件不断涌现，对 EMS软件分析决策的可靠性和精度要求越来越高，研究实用的参数估计方法具有重要的现 实意义。参数辨识中，有载调压变压器变比的估计在EMS中的应用已较为成熟，但支路阻抗 参数估计的实用化研究则相对较少。最直接的方法是将支路参数作为状态变量增广到状态 估计中进行计算，由于量测对参数的求导分量容易导致雅克比矩阵条件数跃升，故增广参 数后状态估计的数值稳定性显著变差。W.-H.E.Liu, RF.\Vu and S.-M.Liu. Estimations of Parameter Errors form Measurement Residuals in State Estimation.正EE Tans, on Power Systems, 1992，(7):81-89.提出了一种改进 方法，即用可疑支路潮流补偿量代替可疑支路参数，使得状态估计获得无偏的估计结果。这种方法只可以在一定程度上改善数值稳定性问题，在现场应用中，采用该方法会经 常出现支路参数不可估计或者估计结果不合理等问题，数值稳定性仍然较差；而且该方法 还存在着局部最优性问题没有得到很好的解决。7/r矿'
技术实现思路
本专利技术的目的是提出，在电网 调度中心，根据实际的电网模型和实时量测数据推导出状态估计残差对电网每条线路参数 的灵敏度，选定可疑支路集，根据序列线性规划方法，修正可疑支路的参数，从而提高应 用软件的可靠性和精确性。本专利技术提出的基于状态估计残差参数灵敏度的电网参数估计方法，包括以下步骤 (1)建立电网最小二乘状态估计模型，并用快速分解法对模型求解，得到状态估计结果,；最小二乘状态估计模型为 残差目标函数本")=s w;化-J + z <(r -^ + 2: -q)2 + i: <(c -w + !>fte: -a)2/eOr ,e 等式约束条件《=V, I^fe -《)+ & S,,-《))风/ € 、G =v,》"Gysi, -《))"。，其中，x为包含《、v,、《、Q、 A和A的变量向量，《、v,、 S、 Q、《和0;分 别为电网中节点/的电压相角、节点/的电压幅值、节点/的有功注入、节点/的无功注入、 节点/到_/的有功潮流和节点/到_/的无功潮流，a为包含q、 5、 6和r的参数向量，G、 5、 6和r分别为电网节点导纳矩阵的实部、电网节点导纳矩阵的虚部、支路对地电容和支路上变压器的变比，v卩、if、 o"、《;T和2;r分别为电网量测系统中的节点z'的电压幅值量测、节点/的有功注入量测、节点!'的无功注入量测、节点f到J'的有功潮流量测和节 点/到7的无功潮流量测，<、<、<、W;和Mf分别为节点/的电压幅值权重、节点/的6节点集、包含无功注入量测的节点集、包含有功潮 资L量测的节点集和包含无功潮流量测的节点集；(2)根据上述状态估计模型和状态估计结果？，计算以下矩阵的值分别为/f如)=H +=o",")+》;、q (x'，力Ca(,xp>=[Vac(x))F其中，J(x，a)为上述状态估计模型中的残差目标函数，c(x，a)为上述状态估计模型中 的等式约束条件，/为等式约束个数，V^、 ▽。、 Vj。和V^分别为对x向量求一阶导数、对 "向量求一阶导数、对x和a向量求二阶导数以及对^向量求二阶导数，^为包含《、v,、 《、Q、《和《的变量向量，"是它的向量维数，"为包含G、 S、 6和r的参数向量， p为它的向量维数；(3)根据上述矩阵的值，计算参数灵敏度^，计算过程为仏=c 05x 3a仏=K。c (4)根据上式得到的参数灵敏度l^，计算实际参数灵敏度^计算残差 /对线路电阻r和电抗X的灵敏度为 a/ a/ a/3Gy. 3& 3G力.观.牝a/ a/ a/朋 朋 朋，观,+ ■3g a/ a/a/祖.+ -a/ a/a/ a;瓶 a/朋,,、朋々朋,'^朋』^乂a/亂朋.,a/广a/ a/ aG,, a/ .+--观,5r"a/ a/眠 a/朋力、朋々气朋o3g,a/ a/a/、观.a/ a/ a/眠asw 乂计算残差 /对于线路充电电容~和线路变压器变比t;的灵敏度为 a/ a/ a/观，a/朋,,a/ a/ a/3~ % 眠^ ^ 36y s~ 朋,，朋』 a/ a/ a/ a/朋,，a/ sg,, a/雄,■+--+-朋朋其中，"为电网的实际参数向量，包括线路电阻^、线路电抗 、线路充电电容^和,一一 ^ a/ a/ a/ a/ a/线路变压器变比乙，^;具体包括^、 ^、 ^和^;(4)利用序列线性优化方法，对电网参数进行修改 (4一1)确定电网的可疑支路集^: 在候选支路中选取残差灵敏度高的支路集合，在该集合中选取支路有关量测与状态估 计结果偏大的支路进入可疑支路集^ ;(4一2)求解如下线性规划模型，得到实际参数偏差AU的值min J(x, w0)+ ^ Aw其中，"n为实际电网参数初值，f^为实际参数灵敏度，A",为第/个参数的本次调整 量，A^和A^是第/个参数的下调和上调的步长限制，W表示可疑支路集合；(4一3)对可疑支路的参数，加上上述Au的值，得到修正后的电网参数估计值； (4一4)重复步骤(4一1) 一 (4一3)，直到前后相邻两次迭代后的状态估计残差之 差小于10e-3为止，得到最后的电网参数估计值。本专利技术提出的基于状态估计残差参数灵敏度的电网参数估计方法，可应用于电网中可疑支路参数的辨识和估计，本专利技术方法不存在数值稳定性和观测性问题，并且可以利用多 组历史数据进行辨识和估计。因此本专利技术方法具有以下优点1、 数值稳定性好；2、 有较好的精确性；3、 计算简单，容易实现，计算速度快。 附图说明图1是算例中用到本专利技术方法的IEEE—9节点系统。 具体实施例方式本专利技术提出的基于状态估计残差参数灵敏度的电网参数估计方法，包括以下步骤 (1)建立电网最小二乘状态估计模型，并用快速分解法对模型求解，得到状态估计结果,；最小二乘状态估计模型为-残差目标函数<formula>formula see original document page 9</formula>等式约束条件<formula>formula see original document page 9</formula>其中，^为包含6",、 Vj、 S、 Q、 A和^的变量向量，《、V。《、2,、 A和A分 别为电网中节点f的电压相角、节点Z的电压幅值、节点/的有功注入、节点z'的无功注入、 节点!'到J'的有功潮流和节点Z到J'的无功潮流，a为包含G、 S、 6和r的参数向量，G、5、 6和r分别为电网节点导纳矩阵的实部、电网节点导纳矩阵的虚部、支路对地电容和支路上变压器的变比，vf、 if、 Qm、 i^和2;本文档来自技高网...
【技术保护点】
一种基于状态估计残差参数灵敏度的电网参数估计方法，其特征在于该方法包括以下步骤： （１）建立电网最小二乘状态估计模型，并用快速分解法对模型求解，得到状态估计结果ｘ↑［ｅ］： 最小二乘状态估计模型为： 残差目标函数： Ｊ（ｘ，ａ）＝＊ｗ↓［ｉ］↑［ｖ］（ｖ↓［ｉ］↑［ｍ］－ｖ↓［ｉ］）↑［２］＋＊ｗ↓［ｉ］↑［Ｐ］（Ｐ↓［ｉ］↑［ｍ］－Ｐ↓［ｉ］）↑［２］＋＊ｗ↓［ｉ］↑［Ｑ］（Ｑ↓［ｉ］↑［ｍ］－Ｑ↓［ｉ］）↑［２］＋＊ｗ↓［ｉｊ］↑［Ｐ］（Ｐ↓［ｉｊ］↑［ｍ］－Ｐ↓［ｉｊ］）↑［２］＋＊ｗ↓［ｉｊ］↑［Ｑ］（Ｑ↓［ｉｊ］↑［ｍ］－Ｑ↓［ｉｊ］）↑［２］ 等式约束条件： Ｐ↓［ｉ］＝ｖ↓［ｉ］＊ｖ↓［ｊ］（Ｇ↓［ｉｊ］ｃｏｓ（θ↓［ｉ］－θ↓［ｊ］）＋Ｂ↓［ｉｊ］ｓｉｎ（θ↓［ｉ ］－θ↓［ｊ］）） ｉ∈Ω↓［Ｐ↓［ｃ］］ Ｑ↓［ｉ］＝ｖ↓［ｉ］＊ｖ↓［ｊ］（Ｇ↓［ｉｊ］ｓｉｎ（θ↓［ｉ］－θ↓［ｊ］）－Ｂ↓［ｉｊ］ｃｏｓ（θ↓［ｉ］－θ↓［ｊ］）） ｉ∈Ω↓［Ｑ↓［ｃ］］ Ｐ↓［ｉｊ］＝ｖ↓ ［ｉ］ｖ↓［ｊ］（Ｇ↓［ｉｊ］ｃｏｓ（θ↓［ｉ］－θ↓［ｊ］）＋Ｂ↓［ｉｊ］ｓｉｎ（θ↓［ｉ］－θ↓［ｊ］））－Ｇ↓［ｉｊ］ｖ↓［ｉ］↑［２］／Ｔ↓［ｉｊ］ ｉ，ｊ∈Ω↓［Ｐ↓［Ｆ］］ Ｑ↓［ｉｊ］＝ｖ↓［ｉ］ｖ↓［ｊ］（Ｇ↓［ ｉｊ］ｓｉｎ（θ↓［ｉ］－θ↓［ｊ］）－Ｂ↓［ｉｊ］ｃｏｓ（θ↓［ｉ］－θ↓［ｊ］））＋ｖ↓［ｉ］↑［２］（Ｂ↓［ｉｊ］－ｂ↓［ｉｊ］／２）／Ｔ↓［ｉｊ］ ｉ，ｊ∈Ω↓［Ｑ↓［Ｆ］］ 其中，ｘ为包含θ↓［ｉ］、ｖ↓［ｉ］、Ｐ↓［ ｉ］、Ｑ↓［ｉ］、Ｐ↓［ｉｊ］和Ｑ↓［ｉｊ］的变量向量，θ↓［ｉ］、ｖ↓［ｉ］、Ｐ↓［ｉ］、Ｑ↓［ｉ］、Ｐ↓［ｉｊ］和Ｑ↓［ｉｊ］分别为电网中节点ｉ的电压相角、节点ｉ的电压幅值、节点ｉ的有功注入、节点ｉ的无功注入、节点ｉ到ｊ的有功潮流和节点ｉ到ｊ的无功潮流，ａ为包含Ｇ、Ｂ、ｂ和Ｔ的参数向量，Ｇ、Ｂ、ｂ和Ｔ分别为电网节点导纳矩阵的实部、电网节点导纳矩阵的虚部、支路对地电容和支路上变压器的变比，ｖ↓［ｉ］↑［ｍ］、Ｐ↓［ｉ］↑［ｍ］、Ｑ↓［ｉ］↑［ｍ］、Ｐ↓［ｉｊ］↑［ｍ］和Ｑ↓［ｉｊ］↑［ｍ］分别为电网量测系统中的节点ｉ的电压幅值量测、节点ｉ的有功注入量测、节点ｉ的无功注入量测、节点ｉ到ｊ的有功潮流量...
【技术特征摘要】
【专利技术属性】
技术研发人员：吴文传，张伯明，孙宏斌，郭庆来，曾兵，
申请(专利权)人：国家电网公司，清华大学，
类型：发明
国别省市：11[中国|北京]
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。