基于元模型的分布式电源就地电压无功控制策略求解方法技术

一种基于元模型的分布式电源就地电压无功控制策略求解方法，包括：将分布式电源接入点电压、分布式电源出力、采样时刻作为输入，分布式电源无功出力作为输出，构建用于分布式电源就地电压无功控制的克里金元模型；在待优化时段t

【技术实现步骤摘要】
基于元模型的分布式电源就地电压无功控制策略求解方法
本专利技术涉及一种分布式电源就地电压无功控制。特别是涉及一种基于元模型的分布式电源就地电压无功控制策略求解方法。
技术介绍
随着分布式电源(distributedgeneration，DG)越来越多的接入配电系统，配电系统的形态结构和运行方式都发生了显著的变化，配电系统从单一的无源系统发展为含大量分布式电源的有源系统。分布式电源的大量接入将给配电系统带来许多益处，诸如降低系统损耗、提高供电可靠性、减少环境污染等，与此同时，光伏(photovoltaic，PV)、风电(windturbine，WT)等间歇性分布式电源也会造成配电网功率波动、电压越限等影响。分布式电源无功优化是配电系统提高电压合格率、降低系统损耗的重要途径，因而受到了广泛的关注。目前，分布式电源的电压无功控制策略主要包括集中式控制和就地控制。其中，集中式控制策略利用配电系统全局信息，统一调配可控资源，实现系统全局优化。但随着分布式电源渗透率提高，大规模分布式电源的接入导致量测数据量急剧增大，这将会给中央控制系统带来沉重的通讯及数据处理负担；另外，集中式控制策略需要安装大量量测、通讯和监控装置，投资较高。而分布式电源的就地控制往往只需依据本地量测信息实现分布式电源无功出力的求解，不需要节点间的信息交流或远程量测，从而减少通信数据量，提高求解效率。现有的间歇式分布式电源就地电压无功控制策略求解方法主要基于本地电压量测及功率量测绘制Q-V控制曲线，以此为依据制定间歇式分布式电源就地电压无功控制策略。但是，该方法无法直接建立本地量测信息与分布式电源无功出力的关系模型，导致无功优化效果不理想。因此，急需一种能够挖掘本地量测信息与分布式电源无功出力关系的分布式电源就地电压无功控制策略的制定方法，实现降低系统损耗，改善系统电压水平的目标。
技术实现思路
本专利技术所要解决的技术问题是，提供一种能够实现降低系统损耗，改善系统电压水平的基于元模型的分布式电源就地电压无功控制策略求解方法。本专利技术所采用的技术方案是：一种基于元模型的分布式电源就地电压无功控制策略求解方法，包括如下步骤：1)将分布式电源接入点电压、分布式电源出力、采样时刻作为输入，分布式电源无功出力作为输出，构建用于分布式电源就地电压无功控制的克里金元模型；2)在待优化时段t*开始时采集分布式电源接入点电压，依据t*时段的输入量求解步骤1)构建的克里金元模型；3)输出步骤2)求解结果，即为分布式电源就地电压无功控制策略。步骤1)所述的克里金元模型表示如下：Q＝fTβ+z(X)式中，X＝[x1，x2，…，xh，…xm]T为输入量构成的矩阵，m为构建克里金元模型所使用的样本数，其中，xh＝[Uh，Ph，th]表示构建克里金元模型时第h个输入量，Uh为构建克里金元模型时第h个输入量中th时刻分布式电源接入点电压；Ph为th时刻分布式电源出力；th为第h个输入量所处的采样时刻；Q＝[q1，q2，…，qh，…qm]T为分布式电源无功出力构成的列向量，qh表示th时刻分布电源的无功出力；fTβ表示回归模型，f为p×m阶回归模型基函数矩阵，β为p×1阶待求解的回归模型基函数系数矩阵，选取常函数作为回归模型基函数，即p＝1，f|1×m＝[1，1，…，1]；z(X)表示期望为0、方差为σ2的随机过程，具有如下特性：E[z(xh)]＝0cov[z(xh)z{xl)]＝σ2R(xh，xl)式中，E[z(xh)]表示z(xh)的期望；cov[z(xh)z(xl)]为z(xh)和z(xl)的协方差，xl是构建克里金元模型时第l个输入量；R(xh，xl)是以θ为参数的相关函数，采用高斯函数作为相关函数：式中，为xh和xl第k维分量之间的距离；θk为待求解的相关函数参数矩阵θ的第k维分量；采用加权最小二乘法以及最大似然估计，分别得到β和σ2的估计值：β＝[fR(X)-1fT]-1fR(X)-1Q式中，R(X)为相关函数构成的关联矩阵，表示为：式中，xm是构建克里金元模型时第m个输入量；β、σ2估计值均与相关函数参数矩阵θ有关，通过极大似然估计得到一个无约束的最优化问题：式中，det[R(X)]表示关联矩阵R(X)的行列式，根据所述无约束的最优化问题确定相关函数参数矩阵θ，然后得到β和σ2的估计值，进而建立用于分布式电源无功就地优化的克里金元模型。步骤2)所述的求解过程表示为：基于克里金元模型，将t*时段分布式电源无功优化结果q*表示为：q*＝β+z(x*)式中，x*＝[U*，Pforecast*，t*]为t*时刻求解克里金元模型的输入量，其中，U*为t*时刻分布式电源接入点电压；Pforecast*为t*时刻分布式电源出力的预测值；β为p×1阶待求解的回归模型基函数系数矩阵，选取常函数作为回归模型基函数，即p＝1；z(x*)表示期望为0、方差为σ2的随机过程；优化结果q*由分布式电源无功出力列向量Q表示：q*＝cTQ式中，c为待求的加权系数向量；由于无偏性的约束限制，t*时段分布式电源无功优化结果q*的两种表达方式满足均方误差最小、差的期望为0的约束条件：式中，E[(cTQ)2-{β+z(x*)}2]为t*时段分布式电源无功优化结果q*的两种表达方式的均方误差；E[cTQ-{β+z(x*)}]为t*时段分布式电源无功优化结果q*的两种表达方式的差的期望；根据所述约束条件，得到加权系数向量c，进而得到分布式电源就地电压无功优化结果：q*＝β+r(X，x*)TR(X)-1(Q-fTβ)式中，r(X，x*)为X与x*之间的空间相关性：r(X，x*)＝[R(x1，x*)，R(x2，x*)，…，R(xm，x*)]T式中，X为构建克里金元模型时所使用的输入量矩阵。本专利技术的基于元模型的分布式电源就地电压无功控制策略求解方法，立足于解决间歇式分布式电源就地电压无功控制策略的求解问题，基于历史数据建立间歇式分布式电源就地电压无功控制元模型，利用量测数据求解该模型，制定间歇式分布式电源就地电压无功控制策略，达到降低系统损耗、改善系统供电质量的目标。附图说明图1是基于元模型的分布式电源就地电压无功控制策略求解方法的全部流程图；图2是改进的IEEE33节点算例结构图；图3是负荷年运行曲线；图4是光伏年出力曲线；图5是测试日光伏与负荷运行曲线；图6是不同场景下节点18处电压水平；图7是不同场景下节点33处电压水平。具体实施方式下面结合实施例和附图对本专利技术的基于元模型的分布式电源就地电压无功控制策略求解方法做出详细说明。(一)本专利技术的基于元模型的分布式电源就地电压无功控制策略求解方法是基于如下步骤实现：1)根据选定的配电系统，输入系统节点数、线路参数、基准负荷、负荷历史运行数据、网络拓扑连接关系，分布式电源接入、位置、容量、历史运行数据，基准电压和基准功率；2)建立分布式电源无功优化模型，包括：设定有源配电系统网络损耗最小为目标函数，分别考虑潮流约束、系统安全运行约束、不可控分布式电源运行约束、辐射型拓扑结构约束，求解上述模型得到分布式电源无功出力；(1)所述的分布式电源无功优化目标函数可表示为：式中，T为运行优化周期；Ωb为配电系统中所有支路集合；rij为支路ij的电阻；It，ij为t时段节点i流本文档来自技高网...

【技术保护点】
1.一种基于元模型的分布式电源就地电压无功控制策略求解方法，其特征在于，包括如下步骤：1)将分布式电源接入点电压、分布式电源出力、采样时刻作为输入，分布式电源无功出力作为输出，构建用于分布式电源就地电压无功控制的克里金元模型；2)在待优化时段t*开始时采集分布式电源接入点电压，依据t*时段的输入量求解步骤1)构建的克里金元模型；3)输出步骤2)求解结果，即为分布式电源就地电压无功控制策略。

【技术特征摘要】
1.一种基于元模型的分布式电源就地电压无功控制策略求解方法，其特征在于，包括如下步骤：1)将分布式电源接入点电压、分布式电源出力、采样时刻作为输入，分布式电源无功出力作为输出，构建用于分布式电源就地电压无功控制的克里金元模型；2)在待优化时段t*开始时采集分布式电源接入点电压，依据t*时段的输入量求解步骤1)构建的克里金元模型；3)输出步骤2)求解结果，即为分布式电源就地电压无功控制策略。2.根据权利要求1所述的基于元模型的分布式电源就地电压无功控制策略求解方法，其特征在于，步骤1)所述的克里金元模型表示如下：Q＝fTβ+z(X)式中，X＝[x1，x2，…，xh，…xm]T为输入量构成的矩阵，m为构建克里金元模型所使用的样本数，其中，xh＝[Uh，Ph，th]表示构建克里金元模型时第h个输入量，Uh为构建克里金元模型时第h个输入量中th时刻分布式电源接入点电压；Ph为th时刻分布式电源出力；th为第h个输入量所处的采样时刻；Q＝[q1，q2，…，qh，…qm]T为分布式电源无功出力构成的列向量，qh表示th时刻分布电源的无功出力；fTβ表示回归模型，f为p×m阶回归模型基函数矩阵，β为p×1阶待求解的回归模型基函数系数矩阵，选取常函数作为回归模型基函数，即p＝1，f|1×m＝[1，1，…，1]；z(X)表示期望为0、方差为σ2的随机过程，具有如下特性：E[z(xh)]＝0cov[z(xh)z(xl)]＝σ2R(xh，xl)式中，E[z(xh)]表示z(xh)的期望；cov[z(xh)z(xl)]为z(xh)和z(xl)的协方差，xl是构建克里金元模型时第l个输入量；R(xh，xl)是以θ为参数的相关函数，采用高斯函数作为相关函数：式中，为xh和xl第k维分量之间的距离；θk为待求解的相关函数参数矩阵θ的第k维分量；采用加权最小二乘法以及最大似然估计，分别...

【专利技术属性】
技术研发人员：李鹏，张传驰，王成山，宋关羽，赵金利，
申请(专利权)人：天津大学，
类型：发明
国别省市：天津,12