首页> 中国专利> 一种基于新能源发电量区间预测的区间潮流计算方法

一种基于新能源发电量区间预测的区间潮流计算方法

页面导航

摘要
著录项
法律信息
说明书
相似文献

摘要

本发明公开了一种基于新能源发电量区间预测的区间潮流计算方法，包括以下步骤(1)采集参数数据发电量历史数据；(2)构造区间二型模糊逻辑系统；(3)初始参数输入到系统中，输出得到区间二型模糊预测集合；(4)构建区间潮流计算模型；(5)将集合作为区间潮流计算模型初次迭代区间，计算得到Krawczyk-Moore算子，并用Krawczyk-Moore算子和初始区间求交集，得到新的区间并将其作为第二次迭代的初始区间，判断该区间宽度是否满足收敛条件，若满足则输出区间，若不满足则返回步骤(4)进行下一次迭代。本发明避免人为规定初始区间所造成的区间迭代的不收敛，收敛过快或过慢问题，从而提高了计算的精确度。

著录项

公开/公告号CN105048451A

专利类型发明专利
公开/公告日2015-11-11

原文格式PDF
申请/专利权人国电南瑞科技股份有限公司;国网宁夏电力公司;国电南瑞南京控制系统有限公司;
展开▼

申请/专利号CN201510372963.9
发明设计人张慧琳;仇式鹍;丁恰;丁茂生;戴则梅;韩红卫;梁峻恺;喻洁;张丙金;
展开▼

申请日2015-06-30
分类号H02J3/00(20060101);G06F19/00(20110101);G06Q50/06(20120101);
代理机构32224 南京纵横知识产权代理有限公司;
代理人董建林
地址 210061 江苏省南京市高新区高新路20号
入库时间 2023-12-18 12:06:53

法律信息

法律状态公告日

法律状态信息

法律状态
2018-04-20

授权

授权
2015-12-09

实质审查的生效 IPC(主分类):H02J3/00 申请日:20150630

实质审查的生效
2015-11-11

公开

公开

说明书

技术领域

本发明涉及一种基于新能源发电量区间预测的区间潮流计算方法，属于新能源发电领域。

背景技术

随着化石能源的逐渐枯竭，对环境污染与气候恶化所产生的重要影响，太阳能以及风电、光伏等可再生能源等新型发电技术的发展，可再生新能源发电成为满足负荷增长需求、减少环境污染、提高能源综合利用效率和供电可靠性的一种有效途径，在电网中得到广泛的应用。可再生能源发电主要利用太阳能、生物质能、风能、水能、波浪能等，受地理条件、天气情况和外部环境等因素的影响，这些可再生电源的发电量输出具有间歇性和随机性，难以得到精确的输出结果，可以通过模糊预测得到发电量预测区间。因而，传统的电力潮流计算成为含区间量的潮流计算，即区间潮流计算。

发明内容

针对现有技术存在的不足，本发明目的是提供一种基于新能源发电量区间预测的区间潮流计算方法，利用区间二型模糊逻辑处理不确定性问题的能力，得到新能源发电量的区间二型模糊预测集合，其中，区间二型模糊预测集合，明确了发电量的波动区间范围，该区间可以为区间潮流计算中的初始迭代区间，避免人为规定初始区间所造成的区间迭代的不收敛，收敛过快，收敛过慢等问题，从而提高了计算的精确度。

为了实现上述目的，本发明是通过如下的技术方案来实现：

本发明的一种基于新能源发电量区间预测的区间潮流计算方法，包括以下几个步骤：(1)采集分布式发电系统的参数数据(以光伏发电系统为例，参数数据为光照强度、环境温度、湿度)以及分布式发电量历史数据(以光伏发电系统为例，历史数据为历史光照强度、历史环境温度、历史湿度和历史光伏电源发电量)； (2)将步骤(1)中的参数数据作为输入，步骤(1)中的分布式发电量历史数据作为输出，构造区间二型模糊逻辑系统，并设置所述区间二型模糊逻辑系统的初始参数(以光伏发电系统为例，初始参数有光照强度、环境温度、湿度)；(3) 训练所述区间二型模糊逻辑系统，即将采集到的实时参数数据(以光伏发电系统为例，将实时光照强度、实时环境温度、实时湿度)作为初始参数(初始光照强度，初始环境温度，初始湿度)输入到所构造的区间二型模糊逻辑系统中，输出得到分布式电源发电量区间二型模糊预测集合；(4)构建区间潮流计算模型；(5) 将步骤(3)得到的分布式电源发电量区间二型模糊预测集合作为区间潮流计算模型初次迭代区间，得到雅克比矩阵的区间扩展，根据Krawczyk-Moore算子定义计算得到Krawczyk-Moore算子，并用所述Krawczyk-Moore算子和初始迭代区间求交集，得到新的区间并将其作为第二次迭代的初始区间，然后判断该区间宽度是否满足收敛条件，若满足收敛条件，则输出区间，若不满足收敛条件则返回步骤(4)进行下一次迭代。

步骤(2)中，所述区间二型模糊逻辑系统的构造方法如下：(2-1)设计模糊器：模糊器即通过一个主隶属度函数得到模糊区间，该主隶属度函数由上、下隶属度函数构成，所述分布式电源发电量区间二型模糊预测集合的主隶属函数选取均方差不确定的高斯函数，上、下隶属度函数如下式所示，模型有三个输入、一个输出，

${\overline{u}}_{x_{k}} (x_{k}) = \exp [- \frac{1}{2} (\frac{x_{k} - x_{k}^{*}}{{\overline{σ}}_{k}})^2]$

${\underline{u}}_{x_{k}} (x_{k}) = \exp [- \frac{1}{2} (\frac{x_{k} - x_{k}^{*}}{{\underline{σ}}_{k}})^2]$

其中，是输入精确值，是均方差变化范围，K＝1,2…p是输入维数， x_k∈X是系统输入；

(2-2)构造规则库：规则前、后件都选取分布式电源发电量区间二型模糊预测集合，主隶属函数是均值不确定的高斯函数，规则形式如下式：

$I F >>>x_{1}>>is>{\tilde{F}}_{1}^{l}>>and......x_{k}>>>>is>>>{\tilde{F}}_{k}^{l}......and>>x_{p}>>is>>{\tilde{F}}_{p}^{l},then>>>y>>>is>>{\tilde{G}}^{l}$

其中是规则前件集合，y∈Y是规则输出，是后件集合，l＝1,2…M,M 是规则总数；

(2-3)构造推理机：推理过程如下式，参与计算的是分布式电源发电量区间二型模糊预测集合的上、下隶属度函数

$u_{{\tilde{G}}^{l}} (y) = \int_{b^{l} \in [{\underline{f}}^{l} * {\underline{u}}_{{\tilde{G}}^{l (y)}}, {\overline{f}}^{l} * {\overline{u}}_{{\tilde{G}}^{l (y)}}]} \frac{1}{b^{l}} y \in Y$

其中：*是t-范数，取最小算子，分别是后件集合的上、下隶属度函数，分别是激活集合的上、下隶属度函数，

$b^{l} \in [{\underline{f}}^{l} * {\underline{u}}_{\tilde{G}}^{l (y)}, {\overline{f}}^{l} * {\overline{u}}_{\tilde{G}}^{l (y)}];$

(2-4)由步骤(2-2)和(2-3)得到此多输入单输出多条模糊规则系统的推理模型为：

取η，u，v为分布式发电系统的参数数据，y为预测发电量，参数数据集合，为发电量集合，均为分布式电源发电量区间二型模糊预测集合，其中，由于采用中心集降型法，选为用质心所表示的区间集合，则的表达式可以写为下述的一个区间数：

$y^{l} = [{\underline{y}}^{l}, {\overline{y}}^{l}]$

(2-5)设计降型器：采用中心集降型方法，将各规则发电量集合用质心来代替，然后求质心的加权平均值，最终得到质心区间，具体表达式为：

$X_{c o s} = [\underline{y}, \overline{y}] = [\frac{Σ_{l = 1}^{L} {\overline{f}}^{l} \underline{y^{l}} + Σ_{l = L + 1}^{M} \underline{f^{l}} \underline{y^{l}}}{Σ_{l = 1}^{L} {\overline{f}}^{l} + Σ_{l = L + 1}^{M} \underline{f^{l}}}, \frac{Σ_{l = 1}^{R} \underline{f^{l}} {\overline{y}}^{l} + {Σ_{l = L +}^{M}}_{1} {\overline{f}}^{l} {\overline{y}}^{l}}{Σ_{l = 1}^{R} \underline{f^{l}} + Σ_{l = R + 1}^{M} {\overline{f}}^{l}}]$

式中：分别是各规则发电量集合质心的上、下界，分别是激活集合的上、下界，L、R是阈值。

步骤(4)中，所述区间潮流计算模型的构建方法如下：

$\tilde{x} = (\begin{matrix} \tilde{θ} \\ \tilde{U} \end{matrix}) = (\begin{matrix} [\underline{θ}, \overline{θ}] \\ [\underline{U}, \overline{U}] \end{matrix})$

$F^{'} (\tilde{x}) = - (\begin{matrix} \tilde{H} & \tilde{N} \\ \tilde{K} & \tilde{L} \end{matrix}) = - (\begin{matrix} [\underline{H}, \overline{H}] & [\underline{N}, \overline{N}] \\ [\underline{K}, \overline{K}] & [\underline{L}, \overline{L}] \end{matrix})$

为雅克比矩阵的区间扩展，代表电压相角的区间形式，代表电压幅值的区间形式，代表电压相角下界上界，是电压幅值的下界上界，是雅克比矩阵元素的区间形式，代表电压相角对有功的影响，是雅克比矩阵元素的区间形式，代表电压幅值对有功影响，是雅克比矩阵元素的区间形式，代表电压相角对无功的影响，是雅克比矩阵元素的区间形式，代表电压幅值对无功的影响分别是上述四个雅克比矩阵元素的区间形式的下界和上界。；

以子阵为例进行分析说明，

${\tilde{H}}_{i i} = \frac{\partial Δ {\tilde{P}}_{i}}{\partial {\tilde{θ}}_{i}} = {\tilde{U}}_{i} Σ_{j = 1, j \neq i}^{n} {\tilde{U}}_{j} (G_{i j} \sin {\tilde{θ}}_{i j} - B_{i j} \cos {\tilde{θ}}_{i j})$

${\tilde{H}}_{i j} = \frac{\partial Δ {\tilde{P}}_{i}}{\partial {\tilde{θ}}_{j}} = - {\tilde{U}}_{i} {\tilde{U}}_{j} (G_{i j} \sin {\tilde{θ}}_{i j} - B_{i j} \cos {\tilde{θ}}_{i j})$

其中，代表节点i处，有功功率的偏移量，代表节点i处的电压相角区间形式，代表节点i处的电压幅值区间形式、代表节点j处的电压幅值区间形式，G_ij代表节点i、j之间阻抗的实部，代表节点i和节点j之间相角差的区间形式，B_ij节点i、j之间阻抗的虚部，代表节点j处的电压相角区间形式。

可看出和具有强相关性，独立考虑会增大区间范围，因此：

${\tilde{H}}_{i j} = \frac{\partial Δ {\tilde{P}}_{i}}{\partial {\tilde{θ}}_{i}} = - {\tilde{U}}_{i} {\tilde{U}}_{j} \sqrt{{G_{i j}}^{2} + {B_{i j}}^{2}} \sin ({\tilde{θ}}_{i j} + δ_{i j}),$

${\tilde{H}}_{i i} = {\tilde{U}}_{i} Σ_{j = 1, j \neq i}^{n} {\tilde{U}}_{j} \sqrt{{G_{i j}}^{2} + {B_{i j}}^{2}} \sin ({\tilde{θ}}_{i j} + δ_{i j}),$

$δ_{i j} = \arctan (- B_{i j} / G_{i j}) .$

步骤(5)中，一次迭代Krawczyk-Moore算子的计算方法如下：

将X_cos作为初始区间X₀

$K = y^{k} - Y^{(k)} f (y^{k}) - [I - Y^{(k)} F^{'} (\tilde{x})] (\tilde{x} - y^{k})$

其中 $y^{k} = m (\tilde{x}),$

$f (y^{k}) = (\begin{matrix} P_{i} - m ({\tilde{U}}_{i}) Σ_{j = 1}^{n} m ({\tilde{U}}_{j}) (G_{i j} c o s (m ({\tilde{θ}}_{i}) - m ({\tilde{θ}}_{j})) + B_{i j} s i n (m ({\tilde{θ}}_{i}) - m ({\tilde{θ}}_{j})) \\ Q_{i} - m ({\tilde{U}}_{i}) Σ_{j = 1}^{n} m ({\tilde{U}}_{j}) (G_{i j} \sin (m ({\tilde{θ}}_{i}) - m ({\tilde{θ}}_{j})) - B_{i j} \cos (m ({\tilde{θ}}_{i}) - m ({\tilde{θ}}_{j})) \end{matrix}),$

$Y^{(k)} = {[m (F^{'} (\tilde{x}))]}^{'},$

I为单位阵，是自变量的区间形式，m为取区间数中点的中点函数，P_i为节点i处的有功功率，Q_i为节点i处的无功功率。

步骤(5)中，利用Krawcyzk-Moore算子与初始区间x₀求交集，得到新的区间x₁：

$(\begin{matrix} {\tilde{x}}^{(k + 1)} = {\tilde{x}}^{(k)} \cap K^{(k)} \\ k = 1, 2, ... \end{matrix})$

其中，k代表了迭代次数，x^k+1,x^k分别代表第k+1和k次自变量的区间，K^k代表第k次的K-M算子。

步骤(5)中，所述收敛条件为 $| {\underline{x}}^{k + 1} - {\underline{x}}^{k} | < ω_{-}$ 且 $| {\overline{x}}^{k + 1} - {\overline{x}}^{k} | < ω$

其中，ω代表收敛系数，分别代表了第k+1和k次自变量的区间上界和下界。

本发明的有益效果在于：与现有考虑不确定性的潮流计算相比，本发明所用区间量来描述不确定量更切合实际情况，利用区间二型模糊逻辑系统来得到初始迭代区间，避免了认为根据经验来设定初始区间造成的迭代收敛的一系列问题，而且节约了迭代时间，更适用于大型系统；Krawczyk-Moore算子的适用，具有全局收敛性，它不仅给出了一个区间解，并且考虑了估计误差；区间二型模糊逻辑系统和区间迭代法的结合使用，不仅可以成功的得到考虑不确定性的潮流计算解，而且解决了区间迭代法本身的一些缺点，具有较强的工程实际使用意义。

附图说明

图1为本发明的基于新能源发电量区间预测的区间潮流计算方法工作流程图；

图2为基于区间二型模糊逻辑系统的发电量预测结构图。

具体实施方式

为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合具体实施方式，进一步阐述本发明。

参见图1，本发明的一种基于新能源发电量区间预测的区间潮流计算方法，包括以下几个步骤：

(1)采用多输入单输出方式，采集分布式发电系统的参数数据以及分布式发电量历史数据；

(2)将参数数据作为输入，发电量作为输出，构造区间二型模糊逻辑系统，设置系统的初始参数；

(3)训练区间二型模糊逻辑系统，将采集到的实时参数数据作为初始参数输入到所构造的区间二型模糊逻辑系统中，输出得到分布式电源发电量区间二型模糊预测集合；

(4)构建区间潮流计算模型；

(5)将步骤(3)得到的分布式电源发电量区间二型模糊预测集合作为区间潮流计算模型初次迭代区间，得到雅克比矩阵的区间扩展，计算得到 Krawczyk-Moore算子，并用Krawczyk-Moore算子和初始区间求交集，得到新的区间并将其作为第二次迭代的初始区间，判断该区间宽度是否满足收敛条件，若满足则输出区间，若不满足则返回步骤(4)进行下一次迭代。

参见图2，区间二型模糊逻辑系统的构造方法如下：

A设计模糊器：模糊器将输入精确值转化为区间二型模糊预测集合，以充分处理电力负荷所具有的强不确定性。区间二型模糊预测集合的主隶属函数选取均方差不确定的高斯函数，上、下隶属度函数如下式所示。模型有3个输入、1个输出。

${\overline{u}}_{x_{k}} (x_{k}) = \exp [- \frac{1}{2} (\frac{x_{k} - x_{k}^{*}}{{\overline{σ}}_{k}})^2]$

${\underline{u}}_{x_{k}} (x_{k}) = \exp [- \frac{1}{2} (\frac{x_{k} - x_{k}^{*}}{{\underline{σ}}_{k}})^2]$

其中，是输入精确值，是均方差变化范围，K＝1,2…p是输入维数。 B构造规则库：规则前、后件都选取区间二型模糊预测集合，主隶属函数是均值不确定的高斯函数，规则形式如下式：

$I F >>>x_{1}>>is>{\tilde{F}}_{1}^{l}>>and......x_{k}>>>>is>>>{\tilde{F}}_{k}^{l}......and>>x_{p}>>is>>{\tilde{F}}_{p}^{l},then>>>y>>>is>>{\tilde{G}}^{l}$

其中x_k∈X是系统输入,是规则前件集合y∈Y是规则输出，是后件集合。 l＝1,2…M,M是规则总数。

C构造推理机：对于区间二型Mamdani模糊模型，推理过程如下式，参与计算的是集合的上、下隶属度函数。

$u_{{\tilde{G}}^{l}} (y) = \int_{b^{l} \in [{\underline{f}}^{l} * {\underline{u}}_{{\tilde{G}}^{l (y)}}, {\overline{f}}^{l} * {\overline{u}}_{{\tilde{G}}^{l (y)}}]} \frac{1}{b^{l}} y \in Y$

其中：*是t-范数，取最小算子，分别是后件集合的上、下隶属度函数，分别是激活集合的上、下隶属度函数。

D由过程(2)和(3)得到此多输入单输出多条模糊规则系统的推理模型为：

取x，u，v为参数数据，y为预测发电量，参数数据集合，为发电量集合。均为区间二型模糊预测集合，其中由于采用中心集降型法选为用质心所表示的区间集合，则的表达式为：

$y^{l} = [{\underline{y}}^{l}, {\overline{y}}^{l}]$

E设计降型器：采用中心集降型方法，将各规则发电量集合用质心来代替，然后求质心的加权平均值，最终得到质心区间。计算中采用Karnik和Mendel提出的简化，具体表达式为：

式中： $\underline{y^{l}}, {\overline{y}}^{l}$ 分别是各规则发电量集合质心的上、下界， $\underline{f^{l}}, {\overline{f}}^{l}$ 分别是激活集合的上、下界，L、R是阈值。

步骤(4)中，区间潮流计算模型的构建方法如下：

$\tilde{x} = (\begin{matrix} \tilde{θ} \\ \tilde{U} \end{matrix}) = (\begin{matrix} [\underline{θ}, \overline{θ}] \\ [\underline{U}, \overline{U}] \end{matrix})$

为雅克比矩阵的区间扩展；

以子阵为例进行分析说明。

${\tilde{H}}_{i j} = \frac{\partial Δ {\tilde{P}}_{i}}{\partial {\tilde{θ}}_{j}} = - {\tilde{U}}_{i} {\tilde{U}}_{j} (G_{i j} \sin {\tilde{θ}}_{i j} - B_{i j} \cos {\tilde{θ}}_{i j})$

可看出和具有强相关性，独立考虑会增大区间范围，因此：

${\tilde{H}}_{i j} = \frac{\partial Δ {\tilde{P}}_{i}}{\partial {\tilde{θ}}_{j}} = - {\tilde{U}}_{i} {\tilde{U}}_{j} \sqrt{{G_{i j}}^{2} + {B_{i j}}^{2}} \sin ({\tilde{θ}}_{i j} + δ_{i j}),$

${\tilde{H}}_{i i} = {\tilde{U}}_{i} Σ_{j = 1, j \neq i}^{n} {\tilde{U}}_{j} \sqrt{{G_{i j}}^{2} + {B_{i j}}^{2}} \sin ({\tilde{θ}}_{i j} + δ_{i j}),$

$δ_{i j} = \arctan (- B_{i j} / G_{i j}) .$

步骤(5)中，一次迭代Krawczyk-Moore算子的计算方法如下：

将X_cos作为初始区间x₀

$K = y^{k} - Y^{(k)} f (y^{k}) - [I - Y^{(k)} F^{'} (\tilde{x})] (\tilde{x} - y^{k})$

其中 $y^{k} = m (\tilde{x}),$

$Y^{(k)} = {[m (F^{'} (\tilde{x}))]}^{'},$

I为单位阵。

步骤(5)中，利用Krawcyzk-Moore算子与初始区间x₀求交集，得到新的区间X：

$(\begin{matrix} {\tilde{x}}^{(k + 1)} = {\tilde{x}}^{(k)} \cap K^{(k)} \\ k = 1, 2, ... \end{matrix})$

步骤(5)中，收敛条件为 $| {\underline{x}}^{k + 1} - {\underline{x}}^{k} | < ω_{-}$ 且 $| {\overline{x}}^{k + 1} - {\overline{x}}^{k} | < ω .$

本实施例中，对于参数数据采用均方差不确定的高斯函数作为其上、下隶属度函数，将输入的精确值模糊化为单值二型模糊集合。区间二型模糊逻辑系统的规则采用“IF-THEN”形式，Mamdani推理模型。由输入和规则前件产生激活集合，再由激活集合与后件集合计算输出，参与计算的是各集合的上、下隶属度函数。推理过程包括计算适配度、求激励强度、求有效的后件隶属度函数和求总输出隶属度函数。采用中心集降型方法，将各规则发电量集合用质心来代替，然后求质心的加权平均值，最终得到质心区间。使用区间量来描述潮流计算中的不确定性，简单而切合实际，而在上文中得到的质心区间可以作为区间潮流计算的初始区间。根据牛顿法的计算公式，以及网络拓扑结构，求得潮流计算中雅克比矩阵各元素的区间拓展。根据区间迭代法和上文中的雅克比矩阵区间扩展，与上文中的初始区间，得到Krawczyk-Moore算子，然后求出K-M算子与初始区间的交集，将其作为第二次迭代的初始区间，并判断是否满足收敛条件。

由于分布式发电机的发电量难以建立精确的数学模型，准确地预测发电量存在很大的难度。本发明采用区间来描述不确定量，简单直接，且比概率潮流更符合实际情况，计算量更小。本发明所述的区间二型模糊逻辑系统无需建立精确的数学模型，利用语言形式的规则库描述不确定性信息逻辑，按照隶属函数模糊程度构造推理机，适合分布式电源发电量不确定性的特点。本发明采用基于牛顿法的区间迭代法，具有较好的鲁棒性，且可以全局收敛；采用从区间二型模糊预测集合作为初值，避免了初值问题引起的对分问题，简化了收敛过程。本发明提出的基于区间二型模糊逻辑系统的区间潮流计算方法，可以进行准确的区间潮流计算，且解决了区间迭代法本身的收敛过快、过慢、不收敛的问题，具有较高的实用价值。

以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。

去获取专利，查看全文>

相似文献

专利
中文文献
外文文献

1. 一种基于新能源发电量区间预测的区间潮流计算方法 [P] . 中国专利： CN105048451B . 2018.04.20
2. 一种基于新能源发电量区间预测的区间潮流计算方法 [P] . 中国专利： CN105048451A . 2015-11-11
3. INTERVAL POWER FLOW CALCULATION METHOD FOR POWER-HEAT INTEGRATED ENERGY SYSTEM [P] . 世界知识产权组织专利： WO2020093296A1 . 2020-05-14

机译：热电综合能源系统的区间潮流计算方法
4. POWER SYSTEM DC INTERVAL POWER FLOW CALCULATION METHOD AND DEVICE [P] . AU2019101018A4 . 2019-10-10

机译：电力系统直流区间潮流计算方法及装置
5. Method for Forecasting Power Composition of New and Renewable Energy Resources based on Supply Plan [P] . 韩国专利： KR101445081B1 . 2014-10-02

机译：基于供应计划的新能源和可再生能源发电量预测方法