一种基于自动微分技术的温度状态估计方法技术

本发明专利技术提出了一种基于自动微分(automatic differentiation,AD)技术的温度状态估计方法。在传统的状态估计计算过程中始终保持线路电阻不变，但因为输电线路电阻是会随着外界温度的变化而变化的，本发明专利技术将温度作为新的状态量引入到状态估计过程中，建立了计及温度变化的状态估计模型。另外，状态估计中需要计算量测函数的雅克比矩阵，由于考虑了支路温度变量，因而手工计算数量巨大的微分函数和编写微分代码，工作过于繁琐且容易出错，所以本发明专利技术还通过使用AD技术替代传统的手工计算雅克比矩阵，减少了人工编写代码的繁琐和出错的机率，有效避免了截断误差，提高了程序的开发效率，而且本发明专利技术易于与已有的状态估计软件相结合。

【技术实现步骤摘要】

专利技术涉及一种基于自动微分技术的温度状态估计方法，属于电力系统运行和控制

技术介绍
作为能量管理系统(Energy Management System，EMS)的核心，电力系统状态估计通过对生数据的处理，获得状态量的最佳估计值。传统的加权最小二乘法(Weighted Least Squares，WLS)状态估计算法估计质量和收敛性能很好，是状态估计的经典解法和理论基础，适应各种类型的量测系统。传统的状态估计计算过程中假定输电线路电阻参数始终不变，以恒定的电网节点导纳矩阵进行计算。但是根据电热耦合的相关研究，输电线路的温度及电阻随着环境温度的不同而变化，所以需要考虑温度对电阻的影响。在常规的状态估计中都忽略了电阻的变化，状态量中也不包含温度这一变化量，导致了计算结果与实际情况间的较大误差。基本加权最小二乘法是电力系统状态估计的最基本方法，该方法模型简单，收敛性能好，估计质量高，是目前应用最为广泛的方法之一。该方法的不足就是计算雅克比矩阵时需要手工推导微分公式和编写微分代码，工作繁琐且容易出错，效率低下。当支路温度作为新增的状态量后，雅克比元素会更加复杂，该方法的缺点也更明显。自动微分(automatic differentiation,AD)技术的出现克服了这一缺点，它是计算机数值计算和分析领域内的一项完全崭新的技术。AD将微分定义为代数运算，与其他微分方法(如数值微分、符号微分)相比，它可以自动计算函数的任意阶导数，减少了人工编写代码的繁琐和出错的机率，从而避免了截断误差，提高了计算精度。
技术实现思路
专利技术目的：现有状态估计过程中并未考虑温度与电阻的关系，缺少相应的处理模型，并且考虑了支路温度这一状态量后，量测函数的雅克比矩阵元素会相当复杂，求导工作量巨大，而本专利技术提出的一种基于自动微分技术的温度状态估计方法能够解决现有技术存在的缺陷。技术方案：本专利技术提出一种基于自动微分技术的温度状态估计方法，其特征在于包括以下步骤：获取电力系统的网络参数，包括：母线编号、名称、补偿电容，输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、 变压器变比和阻抗，系统支路的当前温度；程序初始化，包括：对状态量设置初值、节点次序优化、形成节点导纳矩阵、设置收敛精度λ、分配内存、声明活跃变量、一个规定的参考温度tRef以及相关参数的基准值；输入遥测数据z，包括：电压幅值、发电机有功功率、发电机无功功率、负荷有功功率、负荷无功功率、线路首端有功功率、线路首端无功功率、线路末端有功功率、线路末端无功功率以及支路上的温度量测；根据电热耦合思想建立支路温度与电阻的模型，用公式表示为：其中，R是导体的电阻，RRef是导体在参考温度下的电阻，tRef是参考温度，tF是固定的温度。普通线路上导体的温度t是外界环境温度tAmb与导体上升的温度tRise之和，可表示如下：其中，PLoss是设备内部所有的损耗，tRatedRise是额定的升高温度，PRatedLoss是对应额定的损耗；构造计及温度变量的量测方程；将雅克比矩阵中不变元素的位置和数值存到一个链表中；恢复迭代计数器：k＝1；应用AD技术计算雅克比矩阵中的可变元素，同时读取链表中相应矩阵的不变元素，获得所需的雅克比矩阵；其中含有状态量电压相角θ、电压幅值U和支路温度t的分块扩展雅克比矩阵为：其中，H为雅克比矩阵；P、Q和L分别表示有功、无功和温度相关的自由矢量；求解下述方程得到状态修正量Δx(k)，选取并修正状态量得到x(k+1)：Δx(k)＝[HT(x(k))WH(x(k))]-1HT(x(k))W[z-h(x(k))]x(k+1)＝x(k)+Δx(k)其中，x为状态量，z为量测量，h为非线性量测函数，H(x)为h(x)的雅克比矩阵，T表示矩阵的转置，W为对角权重矩阵，Wii＝1/σi2，σi为标准差；当小于设定的收敛精度λ时，则计算结束，输出状态估计结果，否则修正状态量，计算新温度下的导纳矩阵，进行第k+1次状态估计。有益效果：本专利技术提供一种基于AD技术的温度状态估计方法，通过建立计及支路温度变化的状态估计模型，不仅能够提高状态估计结果的精度，体现电阻随着温度的变化而变化，而且通过AD技术自动计算获得雅克比矩阵，有效避免了截断误差，提高了算法的计算效率和程序的开发效率，具有工程应用前景。附图说明图1：本专利技术方法流程图；图2：本专利技术电力系统支路部分简化的热阻模型。具体实施方式下面结合附图对专利技术的技术流程进行详细说明：通常金属导体的电阻与温度有如下关系：其中，R是导体的电阻，t是导体的温度，RRef是导体在参考温度下的电阻，tRef是参考温度，tF是固定的温度(一般铜线是234.5℃，硬铝线是228.1℃，铝变压器绕组是225.0℃)。著名的热阻模型如图2所示，描绘了设备的热耦合现象。在热阻模型中，设备温度的上升量与设备的损耗近似成线性关系，系数为Rθ，用公式表示如下：其中，tRise是设备高出外界环境的温度，PLoss是设备内部所有的损耗，tRatedRise是额定的升高温度，PRatedLoss是对应额定的损耗。普通线路上导体的温度t是外界环境温度tAmb与导体上升的温度tRise之和，表示为：其中，某条支路ij的PLoss,ij可通过下式计算得到：式中：gij表示节点i和节点j之间的电导，Ui和Uj分别表示节点i和j的电 压幅值，θi和θj分别表示节点i和j的电压相角。对于变压器支路，其温度模型可表示为：其中，I和IRated分别表示变压器支路上的电流和额定电流，n取决于变压器冷却方式：密封式变压器取0.7，自然冷式变压器取0.8，强制风冷变压器取1.0。电力系统的实时运行和控制需要了解系统的真实运行工况，由于测量和传输等方面的原因，得到的生数据难免存在误差，状态估计能在一定程度上提高数据的精度。自1969年美国麻省理工学院的许怀丕(F.C.Schweppe)等人提出了电力系统状态估计的最基本算法——基本加权最小二乘(weighted least squares，WLS)状态估计算法以来，加权最小二乘法成为电力系统状态估计中应用最多的算法。其基本思想是以量测量和量测估计值之差的平方和最小为目标准则的估计方法。本专利技术提出的一种基于AD技术的温度状态估计方法，该方法是建立在非线性WLS状态估计的模型基础上，考虑了支路电阻受温度的影响，运用AD技术计算非线性WLS状态估计方法中的雅克比矩阵。在给定网络接线、支路参数和量测系统的条件下，非线性量测方程可表示为：z＝h(x)+ε式中：x为状态量；z为量测值矢量即遥测数据；h为由基尔霍夫等基本电路定律所建立的量测函数；ε为量测误差。在计及支路温度变化的电力系统状态估计中，量测量配置的类型要比常规潮流多，不仅包括各节点的注入功率量测Pi、Qi，各节点的电压幅值量测Ui，还包括支路的功率量测Pij、Qij、Pji、Qji以及支路的温度量测Lij、Lji。考虑温度影响后，它们的量测函数分别如下所示：节点i的电压幅值：Ui＝Ui节点注入的有功功率和无功功率：式中：Pi和Qi分别表示节点i注入的有功和无功；Ui和Uj分别表示节点i和j的电压幅值；θij是节点i到节点j的电压相角差；Gij(t)和Bi本文档来自技高网...

【技术保护点】
一种基于自动微分技术的温度状态估计方法，其特征在于包括以下步骤：获取电力系统的网络参数，包括：母线编号、名称、补偿电容，输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器变比和阻抗，系统支路的当前温度；程序初始化，包括：对状态量设置初值、节点次序优化、形成节点导纳矩阵、设置收敛精度λ、分配内存、声明活跃变量、一个规定的参考温度tRef以及相关参数的基准值；输入遥测数据z，包括：电压幅值、发电机有功功率、发电机无功功率、负荷有功功率、负荷无功功率、线路首端有功功率、线路首端无功功率、线路末端有功功率、线路末端无功功率以及支路上的温度量测；根据电热耦合思想建立支路温度与电阻的模型，用公式表示为：R=RRef.t+tFtRef+tF]]>其中，R是导体的电阻，RRef是导体在参考温度下的电阻，tRef是参考温度，tF是固定的温度。普通线路上导体的温度t是外界环境温度tAmb与导体上升的温度tRise之和，可表示如下：t=tAmb+tRise=tAmb+(PLossPRatedLoss)tRatedRise]]>其中，PLoss是设备内部所有的损耗，tRatedRise是额定的升高温度，PRatedLoss是对应额定的损耗；构造计及温度变量的量测方程；将雅克比矩阵中不变元素的位置和数值存到一个链表中；恢复迭代计数器：k＝1；应用AD技术计算雅克比矩阵中的可变元素，同时读取链表中相应矩阵的不变元素，获得所需的雅克比矩阵；其中含有状态量电压相角θ、电压幅值U和支路温度t的分块扩展雅克比矩阵为：H=∂P∂θ∂P∂U∂P∂t∂Q∂θ∂Q∂U∂Q∂t∂L∂θ∂L∂U∂L∂t]]>其中，H为雅克比矩阵；P、Q和L分别表示有功、无功和温度相关的自由矢量；求解下述方程得到状态修正量Δx(k)，选取并修正状态量得到x(k+1)：Δx(k)＝[HT(x(k))WH(x(k))]‑1HT(x(k))W[z‑h(x(k))]x(k+1)＝x(k)+Δx(k)其中，x为状态量，z为量测量，h为非线性量测函数，H(x)为h(x)的雅克比矩阵，T表示矩阵的转置，W为对角权重矩阵，σi为标准差；当小于设定的收敛精度λ时，则计算结束，输出状态估计结果，否则修正状态量，计算新温度下的导纳矩阵，进行第k+1次状态估计。...

【技术特征摘要】
1.一种基于自动微分技术的温度状态估计方法，其特征在于包括以下步骤：获取电力系统的网络参数，包括：母线编号、名称、补偿电容，输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器变比和阻抗，系统支路的当前温度；程序初始化，包括：对状态量设置初值、节点次序优化、形成节点导纳矩阵、设置收敛精度λ、分配内存、声明活跃变量、一个规定的参考温度tRef以及相关参数的基准值；输入遥测数据z，包括：电压幅值、发电机有功功率、发电机无功功率、负荷有功功率、负荷无功功率、线路首端有功功率、线路首端无功功率、线路末端有功功率、线路末端无功功率以及支路上的温度量测；根据电热耦合思想建立支路温度与电阻的模型，用公式表示为： R = R Re f . t + t F t Re f + t F ]]>其中，R是导体的电阻，RRef是导体在参考温度下的电阻，tRef是参考温度，tF是固定的温度。普通线路上导体的温度t是外界环境温度tAmb与导体上升的温度tRise之和，可表示如下： t = t A m b + t R i s e = t A m b + ( P L o s s P R a t e d L o s s ) t R a t e d R i s e ]]>其中，PLoss是设备内部所有的损耗，tRatedRise是额定的升高温度，PRatedLoss是对应额定的损耗；构造计及温度变量的量测方程；将雅克比矩阵中不变元素的位置和数值存到一个链表中；恢复迭代计数器：k＝1；应用AD技术计算雅克比矩阵中的可变元素，同时读取链表中相应矩阵的不变元素，获得所需的雅克比矩阵；其中含有状态量电压相角θ、电压幅值U和支路温度t的分块扩展雅克比矩阵为： H = ∂ P ∂ θ ...

【专利技术属性】
技术研发人员：胥峥，何菲，
申请(专利权)人：国家电网公司，江苏省电力公司，江苏省电力公司盐城供电公司，江苏省电力公司泰州供电公司，
类型：发明
国别省市：北京;11