CN111079351A - 考虑风电不确定性的配电网概率潮流获取方法、装置 - Google Patents
考虑风电不确定性的配电网概率潮流获取方法、装置 Download PDFInfo
- Publication number
- CN111079351A CN111079351A CN202010062079.6A CN202010062079A CN111079351A CN 111079351 A CN111079351 A CN 111079351A CN 202010062079 A CN202010062079 A CN 202010062079A CN 111079351 A CN111079351 A CN 111079351A
- Authority
- CN
- China
- Prior art keywords
- network
- wind power
- data
- generator
- probability
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 79
- 238000000034 method Methods 0.000 title claims abstract description 75
- 238000004364 calculation method Methods 0.000 claims abstract description 45
- 230000002457 bidirectional effect Effects 0.000 claims abstract description 44
- 238000012549 training Methods 0.000 claims abstract description 42
- 230000004913 activation Effects 0.000 claims abstract description 15
- 238000010606 normalization Methods 0.000 claims abstract description 14
- 238000013528 artificial neural network Methods 0.000 claims abstract description 7
- 230000006870 function Effects 0.000 claims description 35
- 239000011159 matrix material Substances 0.000 claims description 21
- 238000012545 processing Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 6
- 238000006467 substitution reaction Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000011425 standardization method Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 13
- 230000001186 cumulative effect Effects 0.000 description 10
- 238000004088 simulation Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 241000039077 Copula Species 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 238000011426 transformation method Methods 0.000 description 2
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/048—Activation functions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及考虑风电不确定性的配电网概率潮流获取方法，包括如下步骤：S1、构建双向生成对抗网络的网络结构：编码器、生成器、判别器的网络结构采用全连接层的人工神经网络，且全连接层均采用LeakyReLU激活函数；生成器的输出层使用Tanh函数，判别器的输出层采用Sigmoid激活函数，在判别器的全连接层后添加Dropout层，在所述编码器、生成器、判别器的各层输入前添加分批标准化层；S2、对步骤S1中的双向生成对抗网络进行训练；S3、双向生成对抗网络训练完成后，截取生成器作为生成模型，输入服从高斯分布的一维随机噪声，得到符合原始数据概率分布的风功率数据；S4、将节点负荷及得到的风功率数据输入概率潮流计算模型，计算输出节点电压及支路功率；本发明在考虑风电出力不确定性的环境下，更准确地得到概率潮流计算的结果。
Description
技术领域
本发明属于含分布式电源的配电网概率潮流计算技术领域，尤其涉及考虑风电不确定性的配电网概率潮流获取方法、装置。
背景技术
在能源紧缺和全球变暖的形势下，大量分布式电源接入配电网，使得传统配电网演变为多电源的复杂网络，且以风电为主的分布式电源出力具有随机性、间歇性，这给配电系统的规划运行、经济调度带来不良影响。因此精确地描述风电的不确定性成为一个热门问题。
概率潮流计算可以全面反映不确定性对系统运行状态的影响，是电力系统规划和安全可靠分析的重要工具。其计算方法中的蒙特卡洛模拟法计算精度高且能得到输出随机变量的概率分布，但需要建立准确的输入变量概率模型。现有的概率模型方法分为参数法和非参数法。而参数法需要提前假定研究对象服从的概率分布，准确性和适用性差；非参数法中核函数的选择、参数的求解具有难度。当系统中包含多个输入随机变量时，需要采用空间变换法或Copula函数分析法处理变量之间的相关性。其中空间变换法只能描述某一种相关性，而Copula函数法对于高维随机变量有很大局限性，且两种方法的计算过程复杂繁琐。
总体而言，大部分传统方法先拟合多维随机变量的联合概率分布，而后从中抽样获得多个风机的功率数据。这些方法主要有以下缺陷：联合概率分布的假设是否合理直接关乎计算结果的准确性，由于实际现场环境复杂多样，找到一个具有广泛适用性并且易于求解参数的联合概率分布十分困难。
因此，考虑到这些问题，本发明提供一种考虑风电不确定性的配电网概率潮流计算方法，利用双向生成网络(Bidirectional Generative Adversarial Networks，BIGAN)强大的表达能力，可以找出给定观测数据内部的统计规律，无需概率建模，生成的数据在保证多样性的同时能很好地反映实际发电单元的时空特性，从而克服传统概率模型准确性差、求解参数复杂以及多维随机变量相关性难以描述的难题，在考虑风电出力不确定性的环境下，更准确地得到概率潮流计算的结果。
发明内容
本发明的目的在于克服现有技术的不足，提供一种考虑风电不确定性的配电网概率潮流计算方法，利用双向生成网络(Bidirectional Generative AdversarialNetworks，BIGAN)强大的表达能力，可以找出给定观测数据内部的统计规律，无需概率建模，生成的数据在保证多样性的同时能很好地反映实际发电单元的时空特性，从而克服传统概率模型准确性差、求解参数复杂以及多维随机变量相关性难以描述的难题，在考虑风电出力不确定性的环境下，更准确地得到概率潮流计算的结果。
本发明解决其技术问题是采取以下技术方案实现的：
考虑风电不确定性的配电网概率潮流获取方法，所述方法包括如下步骤：
S1、构建双向生成对抗网络的网络结构：所述双向生成对抗网络包括编码器、生成器、判别器，其中，所述编码器、生成器、判别器的网络结构采用全连接层的人工神经网络，且全连接层均采用LeakyReLU激活函数；生成器的输出层使用Tanh函数，判别器的输出层采用Sigmoid激活函数；在判别器的全连接层后添加Dropout层，在所述编码器、生成器、判别器的各层输入前添加分批标准化层；
S2、对步骤S1中的双向生成对抗网络进行训练；
S3、双向生成对抗网络训练完成后，截取生成器作为生成模型，输入服从高斯分布的一维随机噪声z，得到二维矩阵，将矩阵数据变换成一维的风电功率曲线，并进行反归一化得到符合原始数据概率分布的风功率数据；
S4、将节点负荷及步骤S3中得到的风功率数据输入概率潮流计算模型，采用前推回代法计算输出节点电压及支路功率。
进一步的，所述步骤S2中的双向生成对抗网络的训练方法为：
S201、训练集数据预处理：获取真实的风功率数据，采用max-min标准化方法，将风功率数据映射到[-1，1]区间，并将一维的风功率数据变换成二维的矩阵数据；
S202、根据双向生成对抗网络的目标函数进行迭代训练：在训练过程中，编码器以步骤S201中得到的矩阵数据x为输入，产生样本E(x)，生成器以服从高斯分布的一维随机噪声z为输入，产生样本G(z)，判别器以(x,E(x))和(z,G(z))两组数据作为输入，区分其来自编码器还是生成器；
S203、判断判别器是否能正确区分(x,E(x))和(z,G(z))两组数据，若能，则继续训练过程；否则，训练结束。
进一步的，所述步骤S202中双向生成对抗网络的目标函数为：
式中，V(D,E,G)表示双向生成对抗网络的目标函数；G、D和E表示生成器、判别器和编码器；E[·]表示给定随机变量的期望值；x代表真实的风功率数据，服从真实数据分布PX(x)，PX(x)表示x的概率密度；z代表生成器的输入变量，服从高斯分布PZ(z)，PZ(z)表示z的概率密度；G(z)表示生成器的输出；log表示以10为底的对数操作；D(x,E(x))表示判别器判断(x,E(x))来自于编码器的概率；D(G(z),z)表示判别器判断(z,G(z))来自于生成器的概率。
进一步的，所述步骤S4中的概率潮流计算模型为：
式中，Y是系统节点注入功率，V为节点电压，Z为支路功率；函数f(·)、g(·)是确定性潮流方程。
进一步的，所述步骤S4中负荷的随机性用正态分布表示，均值为原始系统参数值，方差是均值的10％。
考虑风电不确定性的配电网概率潮流获取装置，包括：
双向生成对抗网络的网络结构构建模块，用于构建包括编码器、生成器、判别器的双向生成对抗网络的网络结构，其中，所述编码器、生成器、判别器的网络结构采用全连接层的人工神经网络，且全连接层均采用LeakyReLU激活函数；生成器的输出层使用Tanh函数，判别器的输出层采用Sigmoid激活函数；在判别器的全连接层后添加Dropout层；在所述编码器、生成器、判别器的各层输入前添加分批标准化层；
双向生成对抗网络训练模块，用于对双向生成对抗网络构的网络结构建模块构建的双向生成对抗网络进行训练；
风功率数据生成模块，用于双向生成对抗网络训练完成后，截取生成器作为生成模型，输入服从高斯分布的一维随机噪声z，得到二维矩阵，将矩阵数据反变换成一维的风电功率曲线，并进行反归一化得到符合原始数据概率分布的风功率数据；
概率潮流计算模块，用于将节点负荷及风功率数据获取模块中得到的风功率数据输入概率潮流计算模型，采用前推回代法计算输出节点电压及支路功率。
一种计算设备，包括：
一个或多个处理单元；
存储单元，用于存储一个或多个程序，
其中，当所述一个或多个程序被所述一个或多个处理单元执行，使得所述一个或多个处理单元执行上述的考虑风电不确定性的配电网概率潮流获取方法。
一种具有处理器可执行的非易失的程序代码的计算机可读存储介质，所述计算机程序被处理器执行时实现上述的考虑风电不确定性的配电网概率潮流获取方法的步骤。
本发明的优点和积极效果是：
(1)本发明提出的BIGAN描述风电不确定性的方法无显式概率建模，完全数据驱动，适用于不同的时间和空间下，且产生的样本不仅能精确模拟风电功率的时间相关性、概率分布特性，还能涉及多个风机的空间相关性，解决了传统概率模型的难点和缺陷；
(2)本发明提出的BIGAN模型易于扩展，不仅可用于描述多个风机的出力，也可用于多个光伏的出力，或者风电和光伏的联合出力；
(3)本发明采用BIGAN刻画风电不确定性的概率潮流计算方法可以全面给出输出随机变量的数字特征、概率密度和累积概率等结果，与基于原始数据的时序法相比误差小，有较高的计算精度，且潮流计算的次数减少，提高了计算效率，在实际工程中有较高的应用价值。
附图说明
以下将结合附图和实施例来对本发明的技术方案作进一步的详细描述，但是应当知道，这些附图仅是为解释目的而设计的，因此不作为本发明范围的限定。此外，除非特别指出，这些附图仅意在概念性地说明此处描述的结构构造，而不必要依比例进行绘制。
图1为本发明实施例提供的考虑风电不确定性的配电网概率潮流获取方法中BIGAN的拓扑结构图；
图2为本发明实施例提供的考虑风电不确定性的配电网概率潮流获取方法中BIGAN的网络结构参数示意图；
图3为本发明实施例提供的考虑风电不确定性的配电网概率潮流获取方法中判别器准确率变化趋势示意图；
图4为本发明实施例提供的考虑风电不确定性的配电网概率潮流获取方法中模拟数据和真实数据的风功率曲线图；
图5为本发明实施例提供的考虑风电不确定性的配电网概率潮流获取方法中模拟数据和真实数据的自相关系数图；
图6为本发明实施例提供的考虑风电不确定性的配电网概率潮流获取方法中模拟数据和真实数据的概率分布图；
图7为本发明实施例提供的考虑风电不确定性的配电网概率潮流获取方法中模拟数据和真实数据皮尔逊系数矩阵差的绝对值示意图；
图8为本发明实施例提供的修正后的IEEE 33节点系统示意图；
图9为本发明实施例提供的各个节点的电压均值图；
图10为本发明实施例提供的节点18的电压概率密度图；
图11为本发明实施例提供的节点18的累积概率示意图；
图12为本发明实施例提供的有功损耗的概率密度示意图；
图13为本发明实施例提供的有功损耗的累积概率示意图；
图14为本发明实施例提供的无功损耗的概率密度示意图；
图15为本发明实施例提供的无功损耗的累积概率示意图。
具体实施方式
首先，需要说明的是，以下将以示例方式来具体说明本发明的具体结构、特点和优点等，然而所有的描述仅是用来进行说明的，而不应将其理解为对本发明形成任何限制。此外，在本文所提及各实施例中予以描述或隐含的任意单个技术特征，仍然可在这些技术特征(或其等同物)之间继续进行任意组合或删减，从而获得可能未在本文中直接提及的本发明的更多其他实施例。
需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。
下面结合图1-15具体说明本发明
实施例1
如图1-15所示，考虑风电不确定性的配电网概率潮流获取方法，所述方法包括如下步骤：
S1、构建双向生成对抗网络的网络结构：双向生成对抗网络包括编码器、生成器、判别器，其中，所述编码器、生成器、判别器的网络结构采用全连接层的人工神经网络，且全连接层均采用LeakyReLU激活函数；生成器的输出层使用Tanh函数，判别器的输出层采用Sigmoid激活函数；在判别器的全连接层后添加Dropout层。在各层输入前添加分批标准化层。本实施例中BIGAN的拓扑结构如图1所示，编码器、生成器、判别器的具体网络结构如图2所示，其中LeakyReLU激活函数的alpha是0.2；Dropout层数值为0.25；分批标准化层中动量momentum为0.8；
S2、对步骤S1中的双向生成对抗网络进行训练，训练方法为：
S201、训练集数据预处理：获取训练样本的风功率数据，本实施例以美国某地区相邻的9个风电机组从2012年1月1日到2012年12月31日的实际数据为仿真数据，随机选取90％的数据作为训练集，剩下10％的数据作为测试集。训练样本包括9个风电机组的功率，数据采样间隔是10min，每个风机包括144个数据，因此每个样本一天有1296个数据点，将这1296个点重塑为36×36的矩阵作为训练集。采用max-min标准化方法，将风功率数据映射到[-1，1]区间，并将一维的风功率数据变成二维的矩阵数据；
S202、根据双向生成对抗网络的目标函数进行迭代训练：在训练过程中，编码器以步骤S201中得到的矩阵数据x为输入，产生样本E(x)，生成器以服从高斯分布的一维随机噪声z为输入，产生样本G(z)，判别器以(x,E(x))和(z,G(z))两组数据作为输入，区分其来自编码器还是生成器；
其中，双向生成对抗网络的目标函数为：
式中，V(D,E,G)表示双向生成对抗网络的目标函数；G、D和E表示生成器、判别器和编码器；E[·]表示给定随机变量的期望值；x代表真实的风功率数据，服从真实数据分布PX(x)，PX(x)表示x的概率密度；z代表生成器的输入变量，服从高斯分布PZ(z)，PZ(z)表示z的概率密度；G(z)表示生成器的输出；log表示以10为底的对数操作；D(x,E(x))表示判别器判断(x,E(x))来自于编码器的概率；D(G(z),z)表示判别器判断(z,G(z))来自于生成器的概率；
S203、判断判别器是否能准确区分(x,E(x))和(z,G(z))两组数据，若能，则继续训练过程；否则，训练结束。理想状态下，生成器生成足以“以假乱真”的G(z)，编码器生成类似高斯分布的样本E(x)，对于判别器而言，无法判定区分(x,E(x))和(z,G(z))，即D(x,E(x))＝D(G(z),z)＝0.5。
需要说明的是，训练的批梯度下降(mini-batch)大小为32，优化器选用Adam，学习率设定为0.0002；
图3为判别器的准确率随迭代次数变化的曲线，从中发现：在训练初期，由于网络参数没有训练好，编码器生成的特征向量与高斯分布有很大不同，同时生成器产生的风功率和真实数据有较大偏差，因此判别器很容易区分输入数据对的来源，准确率较高。随着迭代次数的增加，判别器的准确率逐渐下降，最终在0.5上、下范围波动，这说明在训练后期，编码器将真实数据转换成高斯分布，而生成器将服从高斯分布的随机噪声转换成真实数据，产生的风功率曲线和真实曲线高度相似，导致判别器难以区分，达到纳什均衡。
S3、双向生成对抗网络训练完成后，截取生成器作为生成模型，输入服从高斯分布的一维随机噪声z，得到二维矩阵，将矩阵数据变换成一维的风电功率曲线，并进行反归一化得到符合原始数据概率分布的风功率数据；
训练好双向生成对抗网络后，输入100组维度为10的服从高斯分布的随机数可得到100组风功率样本。计算生成的风功率曲线和测试集中真实风功率曲线之间的欧式距离，并选出距离最小的真实曲线和模拟曲线以及对应的自相关系数如图4、5所示，从图4、5的功率曲线中可以看出：尽管测试集不参与训练，生成器却能模拟出和测试集相似的风功率曲线，这说明所提方法能够产生符合实际情况的数据。另外，从自相关函数可以看出，模拟风功率曲线较好地还原了真实风功率曲线的时间相关性。
在风功率短期特性被验证后，分析风电机组出力的长期特性，对366组真实数据和100组生成数据的概率分布进行统计分析，如图6所示。图6说明相比于传统方法，BIGAN生成的数据较好拟合了真实风功率的概率分布，BIGAN学习到了真实数据之间的分布规律。以上风功率的短期特性、自相关系数、概率分布都是对一个风电机组的验证，但仿真数据对象是相邻9个风电机组，它们的出力具有空间相关性，所以选用皮尔逊系数对BIGAN生成数据的相关性进行验证，图7展示了真实和模拟的9个风机皮尔逊系数矩阵差的绝对值；从图中可以发现：皮尔逊系数之差不大于0.086，这说明生成对抗网络在模拟多个风电机组功率的同时，还能很好地捕获风电机组之间的相关性。
综合上述分析，BIGAN模型能同时刻画风功率的时间相关性、概率分布特性，以及多个风电机组的空间相关性，整个过程无需人工干预，完全数据驱动。与传统概率模型相比，产生的风功率数据更能代表实际运行状况，且学习风功率特征所需的时间大幅缩短。
S4、将节点负荷及步骤S3中得到的风功率数据输入概率潮流计算模型，采用前推回代法计算输出节点电压及支路功率；
进一步的，所述步骤S4中的概率潮流计算模型为：
式中，Y是系统节点注入功率，V为节点电压，Z为支路功率；函数f(·)、g(·)是确定性潮流方程。
进一步的，所述步骤S4中负荷的随机性用正态分布表示，均值为原始系统参数值，方差是均值的10％。
作为举例，采用IEEE 33节点系统作为概率潮流计算的算例，其拓扑结构如图8所示；其中节点1为馈线根节点，UN＝12.66kV；从BIGAN生成的9个风电机组数据中选择前4个接入测试系统进行概率潮流计算，在节点10、13、15、31接入风机WT1、WT2、WT3、WT4，各风机采用恒功率因数控制，功率因数为0.9；IEEE 33节点总负荷为3715+j2300kVA，各节点负荷服从正态分布，均值为原始系统参数值，方差是均值的10％。
以BIGAN生成的100*144个样本和服从正态分布的负荷样本为输入随机变量，进行1400次前推回代确定性潮流计算，然后统计分析输出随机变量的均值、概率分布和累积概率。以上方法记为“本文方法”。以366*144个历史风功率数据为输入进行52704次潮流计算，然后对输出随机变量进行统计，简称“时序法”。将时序法作为参考标准，检验本文方法的计算精度。
图9是各节点电压的均值，由图可见，本文方法得到的各节点电压的均值与时序法计算的结果几乎重合，误差非常小。表1给出了节点18电压的均值、标准差以及越下限概率(以最大允许电压偏移7％为标准)。由表1可知，本文方法相比于时序法，电压均值的相对误差小于0.06％，标准差的相对误差小于3.92％，越下限概率的相对误差小于1.94％，可见本文方法的计算精度较高；为进一步观察节点18电压的概率分布，图10、11给出了电压幅值的概率密度和累积概率，由图可知，本文方法得到的节点18电压上下限以及波动变化和时序方法的计算结果非常吻合。
表1
| 电压均值/p.u. | 电压标准差/p.u. | 电压越下限概率 | |
| 时序法 | 0.9322 | 0.0153 | 0.514 |
| 本文方法 | 0.9328 | 0.0147 | 0.504 |
| 相对误差 | 0.06％ | 3.92％ | 1.94％ |
概率潮流计算的输出除了节点电压还有各支路的有功功率、无功功率，为验证本文方法的计算精度，对有功损耗、无功损耗的数字特征、概率分布和累积概率进行对比。表2给出了网络损耗的数字特征，由表可知，有功损耗均值的相对误差小于0.65％，标准差的相对误差小于2.78％；无功损耗均值的相对误差小于0.98％，标准差的相对误差小于4％，网络损耗是系统所有支路损耗的加和，所以每条支路功率的相对误差会更小。图12、13为有功损耗的概率密度和累积概率，图14、15为无功损耗的概率密度和累积概率，从图中可以得出，本文方法计算得到的输出随机变量与时序法计算结果非常一致，计算精度较高。
表2
综合上述分析，本文方法进行概率潮流的计算结果在数字特征、概率分布以及累积概率方面与时序法计算结果非常一致，计算精度较高，且减少了潮流计算的次数，节省了计算时间。
采用双向生成对抗网络描述风电不确定性的配电网概率潮流计算方法，包括：
双向生成对抗网络的网络结构构建模块，用于构建包括编码器、生成器、判别器的双向生成对抗网络的网络结构，其中，所述编码器、生成器、判别器的网络结构采用全连接层的人工神经网络，且全连接层均采用LeakyReLU激活函数；生成器的输出层使用Tanh函数，判别器的输出层采用Sigmoid激活函数；在判别器的全连接层后添加Dropout层；在各层输入前添加分批标准化层；
双向生成对抗网络训练模块，用于对双向生成对抗网络的网络结构构建模块构建的双向生成对抗网络进行训练；
风功率数据生成模块，用于双向生成对抗网络训练完成后，截取生成器作为生成模型，输入服从高斯分布的一维随机噪声z，得到二维矩阵，将矩阵数据反变换成一维的风电功率曲线，并进行反归一化得到符合原始数据概率分布的风功率数据；
概率潮流计算模块，用于将节点负荷及风功率数据获取模块中得到的风功率数据输入概率潮流计算模型，采用前推回代法计算输出节点电压及支路功率。
一种计算设备，包括：
一个或多个处理单元；
存储单元，用于存储一个或多个程序，
其中，当所述一个或多个程序被所述一个或多个处理单元执行，使得所述一个或多个处理单元执行上述的考虑风电不确定性的配电网概率潮流获取方法；需要说明的是，计算设备可包括但不仅限于处理单元、存储单元；本领域技术人员可以理解，计算设备包括处理单元、存储单元并不构成对计算设备的限定，可以包括更多的部件，或者组合某些部件，或者不同的部件，例如计算设备还可以包括输入输出设备、网络接入设备、总线等。
一种具有处理器可执行的非易失的程序代码的计算机可读存储介质，所述计算机程序被处理器执行时实现上述的考虑风电不确定性的配电网概率潮流获取方法；需要说明的是，可读存储介质例如可以是，但不限于，电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合；可读介质上包含的程序可以用任何适当的介质传输，包括，但不限于无线、有线、光缆，RF等等，或者上述的任意合适的组合。例如，可以以一种或多种程序设计语言的任意组合来编写用于执行本发明操作的程序代码，所述程序设计语言包括面向对象的程序设计语言，诸如Java，C++等，还包括常规的过程式程序设计语言，诸如C语言或类似的程序设计语言。程序代码可以完全地在用户计算设备上执行、部分地在用户设备上执行、作为一个独立的软件包执行，或者完全在远程计算设备或服务器上执行。在涉及远程计算设备的情形中，远程计算设备可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)连接到用户计算设备，或者，可以连接到外部计算设备(例如利用因特网服务提供商来通过因特网连接)。
以上实施例对本发明进行了详细说明，但所述内容仅为本发明的较佳实施例，不能被认为用于限定本发明的实施范围。凡依本发明申请范围所作的均等变化与改进等，均应仍归属于本发明的专利涵盖范围之内。
Claims (8)
1.考虑风电不确定性的配电网概率潮流获取方法，其特征在于：所述方法包括如下步骤：
S1、构建双向生成对抗网络的网络结构：所述双向生成对抗网络包括编码器、生成器、判别器，其中，所述编码器、生成器、判别器的网络结构采用全连接层的人工神经网络，且全连接层均采用LeakyReLU激活函数；生成器的输出层使用Tanh函数，判别器的输出层采用Sigmoid激活函数；在判别器的全连接层后添加Dropout层，在所述编码器、生成器、判别器的各层输入前添加分批标准化层；
S2、对步骤S1中的双向生成对抗网络进行训练；
S3、双向生成对抗网络训练完成后，截取生成器作为生成模型，输入服从高斯分布的一维随机噪声z，得到二维矩阵，将矩阵数据变换成一维的风电功率曲线，并进行反归一化得到符合原始数据概率分布的风功率数据；
S4、将节点负荷及步骤S3中得到的风功率数据输入概率潮流计算模型，采用前推回代法计算输出节点电压及支路功率。
2.根据权利要求1所述的考虑风电不确定性的配电网概率潮流获取方法，其特征在于：所述步骤S2中的双向生成对抗网络的训练方法为：
S201、训练集数据预处理：获取真实的风功率数据，采用max-min标准化方法，将风功率数据映射到[-1，1]区间，并将一维的风功率数据变换成二维的矩阵数据；
S202、根据双向生成对抗网络的目标函数进行迭代训练：在训练过程中，编码器以步骤S201中得到的矩阵数据x为输入，产生样本E(x)，生成器以服从高斯分布的一维随机噪声z为输入，产生样本G(z)，判别器以(x,E(x))和(z,G(z))两组数据作为输入，区分其来自编码器还是生成器；
S203、判断判别器是否能正确区分(x,E(x))和(z,G(z))两组数据，若能，则继续训练过程；否则，训练结束。
3.根据权利要求2所述的考虑风电不确定性的配电网概率潮流获取方法，其特征在于：所述步骤S202中双向生成对抗网络的目标函数为：
式中，V(D,E,G)表示双向生成对抗网络的目标函数；G、D和E表示生成器、判别器和编码器；E[·]表示给定随机变量的期望值；x代表真实的风功率数据，服从真实数据分布PX(x)，PX(x)表示x的概率密度；z代表生成器的输入变量，服从高斯分布PZ(z)，PZ(z)表示z的概率密度；G(z)表示生成器的输出；log表示以10为底的对数操作；D(x,E(x))表示判别器判断(x,E(x))来自于编码器的概率；D(G(z),z)表示判别器判断(z,G(z))来自于生成器的概率。
4.根据权利要求1所述的考虑风电不确定性的配电网概率潮流获取方法，其特征在于：所述步骤S4中的概率潮流计算模型为：
式中，Y是系统节点注入功率，V为节点电压，Z为支路功率；函数f(·)、g(·)是确定性潮流方程。
5.根据权利要求1所述的考虑风电不确定性的配电网概率潮流获取方法，其特征在于：所述步骤S4中负荷的随机性用正态分布表示，均值为原始系统参数值，方差是均值的10％。
6.考虑风电不确定性的配电网概率潮流获取装置，其特征在于：包括：
双向生成对抗网络的网络结构构建模块，用于构建包括编码器、生成器、判别器的双向生成对抗网络的网络结构，其中，所述编码器、生成器、判别器的网络结构采用全连接层的人工神经网络，且全连接层均采用LeakyReLU激活函数；生成器的输出层使用Tanh函数，判别器的输出层采用Sigmoid激活函数；在判别器的全连接层后添加Dropout层；在所述编码器、生成器、判别器的各层输入前添加分批标准化层；
双向生成对抗网络训练模块，用于对双向生成对抗网络构的网络结构建模块构建的双向生成对抗网络进行训练；
风功率数据生成模块，用于双向生成对抗网络训练完成后，截取生成器作为生成模型，输入服从高斯分布的一维随机噪声z，得到二维矩阵，将矩阵数据反变换成一维的风电功率曲线，并进行反归一化得到符合原始数据概率分布的风功率数据；
概率潮流计算模块，用于将节点负荷及风功率数据获取模块中得到的风功率数据输入概率潮流计算模型，采用前推回代法计算输出节点电压及支路功率。
7.一种计算设备，其特征在于：包括：
一个或多个处理单元；
存储单元，用于存储一个或多个程序，
其中，当所述一个或多个程序被所述一个或多个处理单元执行，使得所述一个或多个处理单元执行如权利要求1至5中任一项所述的方法。
8.一种具有处理器可执行的非易失的程序代码的计算机可读存储介质，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至5任意一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010062079.6A CN111079351B (zh) | 2020-01-19 | 2020-01-19 | 考虑风电不确定性的配电网概率潮流获取方法、装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010062079.6A CN111079351B (zh) | 2020-01-19 | 2020-01-19 | 考虑风电不确定性的配电网概率潮流获取方法、装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111079351A true CN111079351A (zh) | 2020-04-28 |
| CN111079351B CN111079351B (zh) | 2024-02-06 |
Family
ID=70323902
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010062079.6A Active CN111079351B (zh) | 2020-01-19 | 2020-01-19 | 考虑风电不确定性的配电网概率潮流获取方法、装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111079351B (zh) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111563275A (zh) * | 2020-07-14 | 2020-08-21 | 中国人民解放军国防科技大学 | 一种基于生成对抗网络的数据脱敏方法 |
| CN111794741A (zh) * | 2020-08-11 | 2020-10-20 | 中国石油天然气集团有限公司 | 滑动定向钻井模拟机实现方法 |
| CN111814403A (zh) * | 2020-07-16 | 2020-10-23 | 国网山东省电力公司电力科学研究院 | 一种配电主设备分布式状态传感器可靠性评估方法 |
| CN112100920A (zh) * | 2020-09-15 | 2020-12-18 | 东南大学 | 一种配电网三相电压计算方法、装置、设备及存储介质 |
| CN113240105A (zh) * | 2021-03-30 | 2021-08-10 | 浙江大学 | 一种基于图神经网络池化的电网稳态判别方法 |
| CN113688919A (zh) * | 2021-08-30 | 2021-11-23 | 华北电力大学（保定） | 一种基于SeqGAN模型的风电机组健康状态评估数据集构建方法 |
| CN114021437A (zh) * | 2021-10-26 | 2022-02-08 | 清华大学 | 一种风电光伏有功场景生成方法、装置、电子设备和存储介质 |
| CN115146827A (zh) * | 2022-06-08 | 2022-10-04 | 国网江苏省电力有限公司淮安供电分公司 | 一种考虑量测缺失的配电网在线优化方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120078436A1 (en) * | 2010-09-27 | 2012-03-29 | Patel Sureshchandra B | Method of Artificial Nueral Network Loadflow computation for electrical power system |
| WO2018157691A1 (zh) * | 2017-02-28 | 2018-09-07 | 国网江苏省电力公司常州供电公司 | 有源配电网的安全性量化方法 |
| CN108599172A (zh) * | 2018-05-18 | 2018-09-28 | 广东电网有限责任公司佛山供电局 | 一种基于人工神经网络的输配网全局潮流计算方法 |
| CN109117951A (zh) * | 2018-01-15 | 2019-01-01 | 重庆大学 | 基于bp神经网络的概率潮流在线计算方法 |
-
2020
- 2020-01-19 CN CN202010062079.6A patent/CN111079351B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120078436A1 (en) * | 2010-09-27 | 2012-03-29 | Patel Sureshchandra B | Method of Artificial Nueral Network Loadflow computation for electrical power system |
| WO2018157691A1 (zh) * | 2017-02-28 | 2018-09-07 | 国网江苏省电力公司常州供电公司 | 有源配电网的安全性量化方法 |
| CN109117951A (zh) * | 2018-01-15 | 2019-01-01 | 重庆大学 | 基于bp神经网络的概率潮流在线计算方法 |
| CN108599172A (zh) * | 2018-05-18 | 2018-09-28 | 广东电网有限责任公司佛山供电局 | 一种基于人工神经网络的输配网全局潮流计算方法 |
Non-Patent Citations (1)
| Title |
|---|
| 范宏;左路浩;张节潭;王东方: "考虑光伏出力不确定性的输电网概率潮流计算", 电力系统及其自动化学报, vol. 29, no. 11 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111563275A (zh) * | 2020-07-14 | 2020-08-21 | 中国人民解放军国防科技大学 | 一种基于生成对抗网络的数据脱敏方法 |
| CN111814403B (zh) * | 2020-07-16 | 2023-07-28 | 国网山东省电力公司电力科学研究院 | 一种配电主设备分布式状态传感器可靠性评估方法 |
| CN111814403A (zh) * | 2020-07-16 | 2020-10-23 | 国网山东省电力公司电力科学研究院 | 一种配电主设备分布式状态传感器可靠性评估方法 |
| CN111794741A (zh) * | 2020-08-11 | 2020-10-20 | 中国石油天然气集团有限公司 | 滑动定向钻井模拟机实现方法 |
| CN111794741B (zh) * | 2020-08-11 | 2023-08-18 | 中国石油天然气集团有限公司 | 滑动定向钻井模拟机实现方法 |
| CN112100920A (zh) * | 2020-09-15 | 2020-12-18 | 东南大学 | 一种配电网三相电压计算方法、装置、设备及存储介质 |
| CN113240105B (zh) * | 2021-03-30 | 2022-01-07 | 浙江大学 | 一种基于图神经网络池化的电网稳态判别方法 |
| CN113240105A (zh) * | 2021-03-30 | 2021-08-10 | 浙江大学 | 一种基于图神经网络池化的电网稳态判别方法 |
| CN113688919A (zh) * | 2021-08-30 | 2021-11-23 | 华北电力大学（保定） | 一种基于SeqGAN模型的风电机组健康状态评估数据集构建方法 |
| CN114021437A (zh) * | 2021-10-26 | 2022-02-08 | 清华大学 | 一种风电光伏有功场景生成方法、装置、电子设备和存储介质 |
| CN114021437B (zh) * | 2021-10-26 | 2024-04-12 | 清华大学 | 一种风电光伏有功场景生成方法、装置、电子设备和存储介质 |
| CN115146827A (zh) * | 2022-06-08 | 2022-10-04 | 国网江苏省电力有限公司淮安供电分公司 | 一种考虑量测缺失的配电网在线优化方法 |
| CN115146827B (zh) * | 2022-06-08 | 2023-11-03 | 国网江苏省电力有限公司淮安供电分公司 | 一种考虑量测缺失的配电网在线优化方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111079351B (zh) | 2024-02-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111079351A (zh) | 考虑风电不确定性的配电网概率潮流获取方法、装置 | |
| CN109117951B (zh) | 基于bp神经网络的概率潮流在线计算方法 | |
| CN111860982A (zh) | 一种基于vmd-fcm-gru的风电场短期风电功率预测方法 | |
| CN107104442B (zh) | 计及参数模糊性的含风电场电力系统概率潮流计算方法 | |
| CN112100911B (zh) | 一种基于深度bilstm的太阳辐射预测方法 | |
| CN111415010A (zh) | 一种基于贝叶斯神经网络的风电机组参数辨识方法 | |
| CN111900713A (zh) | 网源协调下考虑负荷和风电随机性多场景输电网规划方法 | |
| CN111797132B (zh) | 考虑时空相关性的多可再生能源电站功率场景生成方法 | |
| CN105656084A (zh) | 一种计及新能源发电预测误差的改进随机潮流算法 | |
| CN116722545B (zh) | 基于多源数据的光伏发电预测方法和相关设备 | |
| CN117973644B (zh) | 一种考虑参考电站优化的分布式光伏功率虚拟采集方法 | |
| CN116707331A (zh) | 基于模型预测的逆变器输出电压高精度调节方法及系统 | |
| CN111192158A (zh) | 一种基于深度学习的变电站日负荷曲线相似度匹配方法 | |
| Liao et al. | Scenario generations for renewable energy sources and loads based on implicit maximum likelihood estimations | |
| CN110991741B (zh) | 一种基于深度学习的断面约束概率预警方法及系统 | |
| Chen et al. | Variation-cognizant probabilistic power flow analysis via multi-task learning | |
| CN116128211A (zh) | 基于风光不确定性预报场景的风光水联合短期优化调度方法 | |
| Le et al. | Probabilistic assessment of power systems with renewable energy sources based on an improved analytical approach | |
| CN115828441A (zh) | 基于条件生成对抗网络的风机出力场景生成和缩减方法 | |
| CN110852597B (zh) | 基于生成对抗网络的用电高峰时段居民负荷占比计算方法 | |
| CN109494747B (zh) | 一种基于交替梯度算法的电网概率潮流计算方法 | |
| Ji et al. | Probabilistic optimal power flow considering the dependence of multiple wind farms using pair diffusive kernel copula | |
| CN111753252A (zh) | 一种基于Nataf变换的随机变量样本生成方法和系统 | |
| Yang et al. | Hybrid data and model‐driven joint identification of distribution‐network topology and parameters | |
| CN104037756A (zh) | 一种含复杂电力设备模型的电力系统稳定评估方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |