CN109599872B - 基于堆栈降噪自动编码器的电力系统概率潮流计算方法 - Google Patents
基于堆栈降噪自动编码器的电力系统概率潮流计算方法 Download PDFInfo
- Publication number
- CN109599872B CN109599872B CN201811633643.4A CN201811633643A CN109599872B CN 109599872 B CN109599872 B CN 109599872B CN 201811633643 A CN201811633643 A CN 201811633643A CN 109599872 B CN109599872 B CN 109599872B
- Authority
- CN
- China
- Prior art keywords
- layer
- dae
- power flow
- optimal power
- sdae
- 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.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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)
- Power Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了基于堆栈降噪自动编码器的电力系统概率最优潮流计算方法,主要步骤为:1)建立SDAE最优潮流模型。2)获取SDAE最优潮流模型输入层的输入样本X。3)对SDAE最优潮流模型进行初始化。4)对SDAE最优潮流模型进行训练,从而得到训练后的SDAE最优潮流模型。5)采用MCS法对待计算概率潮流的电力系统的随机变量进行抽样,从而获取计算样本。6)将步骤5得到的训练样本数据一次性输入步骤4中训练完成的SDAE最优潮流模型中,从而计算出最优潮流在线概率。7)对所述最优潮流在线概率进行分析,即绘制SDAE最优潮流模型的输出变量的概率密度曲线。本发明可广泛应用于电力系统的概率最优潮流求解,特别适用于新能源渗透率高导致系统不确定性增强的在线分析情况。
Description
技术领域
本发明涉及电力系统及其自动化领域,具体是基于堆栈降噪自动编码器的电力系统概率最优潮流计算方法。
背景技术
随着可再生能源的日益普及,电力系统的不确定性激增。概率最优潮流(POPF)可计及电力系统运行中的各种不确定性因素,已成为确保电力系统安全经济运行的重要工具。现有的概率最优潮流求解技术一般可以分为解析法和模拟法。前者仅适用于某些类型的概率分布,在实际应用中无法处理概率最优潮流的一般情况。后者计算结果精确,应用灵活,但涉及大量的抽样且需要重复求解非线性和非凸的最优潮流(OPF)问题,导致计算时间较长。计算耗时已成为概率最优潮流在电力行业实际应用中的主要瓶颈。
发明内容
本发明的目的是解决现有技术中存在的问题。
为实现本发明目的而采用的技术方案是这样的,基于堆栈降噪自动编码器的电力系统概率最优潮流计算方法,主要包括以下步骤:
1)建立SDAE最优潮流模型。
所述SDAE最优潮流模型包括n个顺序堆叠的DAE模型。其中,第l-1层DAE的输入层为第l-2层DAE的中间层。第l-1层DAE的中间层为第l层DAE的输入层。
所述DAE模型主要包括1个输入层、1个中间层和1个输出层。
输入层的神经元个数设为A。任意输入层神经元记为uI,I=1、 2…A。
输入层的输入样本X中的数据主要包括电力系统中所有新能源节点和负荷节点的有功功率和无功功率。
中间层的神经元个数设为M。任意输入层神经元记为vi,i=1、 2…M。
输出层的神经元个数设为Q。任意输出层神经元记为me,e=1、 2…Q。
输出层的输出向量y中的数据主要包括电力系统节点电压幅值和相角、发电机有功功率、发电机有功功率无功功率、支路有功功率、支路无功功率和发电成本。
建立SDAE最优潮流模型的主要步骤如下:
式中,qD为随机局部腐蚀过程,即随机选取一定数量的输入变量置零。
1.2)利用编码函数fθ得到中间层a。中间层a如下所示:
激活函数s如下所示:
s(x)=max(0,x)。 (3)
式中,x为输入向量X中的数据。
1.3)利用解码函数gθ′得到DAE的输出层Z。输出层Z如下所示:
Z=gθ′(a)=s(W′a+b′)。 (4)
式中,W′为解码函数的权值。W′=WT。b′为解码函数的偏置。 s为激活函数。gθ′为解码函数。
1.4)将DAE逐层堆叠,从而得到SDAE最优潮流模型。
1.5)提取最优潮流输入数据X的高维特征,拟合得到输出Yt。输出Yt如下所示:
其中,fθ (l)为第l层DAE的编码函数,l=1,2,…,n,n为SDAE 最优潮流模型中DAE的个数。fθ (t)为顶层编码函数。X为输入样本。
2)获取SDAE最优潮流模型输入层的输入样本X。
3)对SDAE最优潮流模型进行初始化。
对SDAE最优潮流模型进行初始化的主要步骤如下:
3.1)利用最大最小法对输入向量X和输出向量y进行归一化处理
式中,x为输入向量X中的数据。xmin为输入向量X中的最小数据。xmax为输入向量X中的最大数据。
利用最大最小法对输出向量y进行归一化处理。
式中,y为输出向量。ymin为最小输出向量。ymax为最大输出向量。
3.2)将处理后的训练输入样本X分成m个批量。
3.3)根据电力系统的规模和复杂程度设定SDAE最优潮流模型的总层数n、每层神经元的个数和学习率η。
4)对SDAE最优潮流模型进行训练,从而得到训练后的SDAE 最优潮流模型。
对SDAE最优潮流模型进行训练的主要步骤如下:
4.1)对SDAE最优潮流模型进行无监督预训练,主要步骤如下:
4.1.1)根据第l层DAE的输入Xl和输出Zl构造均方差损失函数 LH(Xl,Zl)。
式中,Xl为第l层DAE的输入,也即第l-1层DAE的中间层输出al-1。Zl是第l层DAE的输出。LH(Xl,Zl)为根据第l层DAE的输入 Xl和输出Zl构造的均方差损失函数。
4.1.3)利用RMSProp学习算法和动量学习率,构建SDAE最优潮流模型无监督预训练参数更新公式。无监督预训练参数更新公式如公式(9)至公式(12)所示。
式中,是第T次参数更新后第l-1层DAE中间层的第j个神经元到第l层DAE中间层的第i个神经元的权值。是第T-1 次参数更新后第l-1层DAE中间层的第j个神经元到第l层DAE中间层的第i个神经元的权值。
式中,η为学习速率。m为此批量样本数量。p为动量因子。Ρ为梯度累积指数。ε为常数。⊙是Hadamard乘子。LH(Zl,Xl)为根据第l个DAE的输出层输出和训练样本输入X构造均方差损失函数。k为任意样本。是第T-1次参数更新后第l-1层 DAE中间层的第j个神经元到第l层DAE中间层的第i个神经元的权值。是第T次参数更新后第l-1层DAE中间层的第j个神经元到第l层DAE中间层的第i个神经元的权值。r和r+m分别是此批量的起始样本序号。为前T次权值迭代积累的梯度。为前T-1次权值迭代积累的梯度。Δ为增量符号。d为微分符号。σ为常数。为偏导符号。
式中,LH(Yt,Y)为根据顶层输出和训练样本输出 Y构造均方差损失函数。k为任意样本。是第T次参数更新后,第l层DAE中间层的第i个神经元的偏移量。是第T-1次参数更新后第l-1层DAE中间层的第j个神经元到第l层DAE中间层的第i个神经元的偏移量。为前T次偏置迭代积累的梯度;为前T-1次偏置迭代积累的梯度。
4.1.4)根据无监督预训练参数更新公式计算每层DAE的最优编码参数θ={W,b},并将最优编码参数θ={W,b}作为有监督微调的初始编码参数。
4.2)对SDAE最优潮流模型进行有监督微调,主要步骤如下:
argθminJ(W,b)=argθminLH(Yt,Y)。 (13)
4.2.2)根据优化目标函数argθminJ(W,b),从而对SDAE最优潮流模型最优编码参数θ={W,b}进行微调。
5)采用MCS法对待计算概率潮流的电力系统的随机变量进行抽样,从而获取计算样本。所述随机变量主要包括待计算概率最优潮流的电力系统的风速、光照辐射度和负荷。
6)将步骤5)得到的训练样本数据一次性输入步骤4)中训练完成的SDAE最优潮流模型中,从而计算出最优潮流在线概率。
7)对所述最优潮流在线概率进行分析,即绘制SDAE最优潮流模型的输出变量的概率密度曲线。
本发明的技术效果是毋庸置疑的。本发明选择堆栈降噪自动编码器(SDAE)来学习最优潮流的优化过程并进一步结合蒙特卡洛模拟法实现概率最优潮流的在线计算。本发明建立的基于SDAE的最优潮流模型能够有效挖掘最优潮流模型的高维非线性特征,快速准确地得到最优潮流的计算结果,体现了SDAE模型对最优潮流模型的强大逼近能力。
本发明提出的基于SDAE与MCS的概率最优潮流在线算法,可有效处理电力系统中的各种不确定性因素,通过MCS法对电力系统状态进行抽样,再由SDAE最优潮流模型直接映射出所有抽样样本的最优潮流计算结果,从而在不增加硬件成本的前提下实现了概率最优潮流的高精度在线计算。
本发明可广泛应用于电力系统的概率最优潮流求解,特别适用于新能源渗透率高导致系统不确定性增强的在线分析情况。
附图说明
图1为DAE的结构图;
图2为基于SDAE的最优潮流结构图;
图3的(a)图为基于内点法的MCS法(M0)求得的最优潮流概率密度曲线与基于堆栈降噪自动编码器的电力系统概率最优潮流计算方法 (M1)求得的最优潮流概率密度曲线在节点10的电压幅值对比图;
图3的(b)图为基于内点法的MCS法(M0)求得的最优潮流概率密度曲线与基于堆栈降噪自动编码器的电力系统概率最优潮流计算方法 (M1)求得的最优潮流概率密度曲线在节点10的发电机有功出力对比图;
图3的(c)图为基于内点法的MCS法(M0)求得的最优潮流概率密度曲线与基于堆栈降噪自动编码器的电力系统概率最优潮流计算方法 (M1)求得的最优潮流概率密度曲线在节点10的支路5-6的有功功率对比图;
图3的(d)图为基于内点法的MCS法(M0)求得的最优潮流概率密度曲线与基于堆栈降噪自动编码器的电力系统概率最优潮流计算方法 (M1)求得的最优潮流概率密度曲线的目标值对比图。
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
实施例1:
参见图1和图2,基于堆栈降噪自动编码器的电力系统概率最优潮流计算方法,主要包括以下步骤:
1)建立SDAE最优潮流模型。利用SDAE具有深层堆栈结构以及编码解码过程从而能够有效挖掘非线性最优潮流模型的高阶特征的特点,建立基于SDAE的最优潮流模型。
考虑到最优潮流模型含有非线性等式、不等式约束从而导致输入输出之间非线性特征复杂,提出了结合最大最小归一化数据预处理方法和基于动量学习率的小批量梯度下降学习算法的深度神经网络训练方法,以提高训练精度与速度。训练后的SDAE最优潮流模型能够非迭代地直接映射出由MCS法生成的随机样本的最优潮流计算结果而无需优化,具有计算速度快、精度高的特点。最后,在 IEEE118系统上进行仿真,采用基于SDAE的最优潮流模型仅需数秒便可准确求得系统概率最优潮流,计算速度与传统的基于内点法的MCS法相比提高了上千倍。
SDAE最优潮流模型由降噪自动编码器(Denoising Auto-Encoders,DAE)逐层堆栈构成。DAE是自动编码器(AE)的一种扩展,是加入了随机因素的自动编码器,通过局部腐蚀输入X (即引入随机噪声)的方式以强制自动编码器全面地提取输入的高维特征,从而实现腐蚀后的输入的重构,由此提高了模型的鲁棒性。
所述SDAE最优潮流模型包括n个顺序堆叠的DAE模型。其中,第l-1层DAE的输入层为第l-2层DAE的中间层。第l-1层DAE的中间层为第l层DAE的输入层。DAE输出层Z并不参与SDAE的数据流通。
所述DAE模型主要包括1个输入层、1个中间层和1个输出层。
输入层的神经元个数设为A。任意输入层神经元记为uI,I=1、 2…A。
输入层的输入样本X中的数据主要包括电力系统中所有新能源节点和负荷节点的有功功率和无功功率。
中间层的神经元个数设为M。任意输入层神经元记为vi,i=1、 2…M。
输出层的神经元个数设为Q。任意输出层神经元记为me,e=1、 2…Q。
A、M和Q的取值由电力系统的规模和复杂程度决定。
输出层的输出向量y中的数据主要包括电力系统节点电压幅值和相角、发电机有功功率、发电机有功功率无功功率、支路有功功率、支路无功功率和发电成本。
建立SDAE最优潮流模型的主要步骤如下:
式中,qD为随机局部腐蚀过程,即随机选取一定数量的输入变量置零。
1.2)利用编码函数fθ得到中间层a。中间层a如下所示:
编码函数的权值W是一个dy×dx维的矩阵。编码函数的偏置b 是一个dy维的向量。本发明选择目前已获得广泛应用的 ReLU(Rectified Linear Unit)函数作为编码和解码过程激活函数。dx和dy分别为输入层和中间层向量的维度。
激活函数s如下所示:
s(x)=max(0,x)。 (3)
式中,x为输入向量X中的数据。
1.3)利用解码函数gθ′得到DAE的输出层Z。输出层Z如下所示:
Z=gθ′(a)=s(W′a+b′)。 (4)
式中,W′为解码函数的权值。W′=WT。b′为解码函数的偏置。 s为激活函数。gθ′为解码函数。
1.4)将DAE逐层堆叠,从而得到SDAE最优潮流模型。
1.5)提取最优潮流输入数据X的高维特征,拟合得到输出Yt。输出Yt如下所示:
2)获取SDAE最优潮流模型输入层的输入样本X。
3)对SDAE最优潮流模型进行初始化。
对SDAE最优潮流模型进行初始化的主要步骤如下:
3.1)利用最大最小法对输入向量X和输出向量y进行归一化处理
式中,x为输入向量X中的数据。xmin为输入向量X中的最小数据。xmax为输入向量X中的最大数据。
利用最大最小法对输出向量y进行归一化处理。
式中,y为输出向量。ymin为最小输出向量。ymax为最大输出向量。
3.2)将处理后的训练输入样本X分成m个批量。
3.3)根据电力系统的规模和复杂程度设定SDAE最优潮流模型的总层数n、每层神经元的个数和学习率η。
4)对SDAE最优潮流模型进行训练,从而得到训练后的SDAE 最优潮流模型。
对SDAE最优潮流模型进行训练的主要步骤如下:
4.1)对SDAE最优潮流模型进行无监督预训练,主要步骤如下:
4.1.1)根据第l层DAE的输入Xl和输出Zl构造均方差损失函数 LH(Xl,Zl)。
式中,Xl为第l层DAE的输入,也即第l-1层DAE的中间层输出al-1。Zl是第l层DAE的输出。LH(Xl,Zl)为根据第l层DAE的输入Xl和输出Zl构造的均方差损失函数。
4.1.3)利用RMSProp学习算法和动量学习率,构建SDAE最优潮流模型无监督预训练参数更新公式。无监督预训练参数更新公式如公式(9)至公式(12)所示。
式中,是第T次参数更新后第l-1层DAE中间层的第j个神经元到第l层DAE中间层的第i个神经元的权值。是第T-1 次参数更新后第l-1层DAE中间层的第j个神经元到第l层DAE中间层的第i个神经元的权值。
式中,η为学习速率。m为此批量样本数量。p为动量因子。ρ为梯度累积指数,是一个接近于1的值,本实施例取ρ=0.999。ε是一个接近于0分值,本实施例取ε=10-8。⊙是Hadamard乘子。LH(Zl,Xl) 为根据第l个DAE的输出层输出和训练样本输入X构造均方差损失函数。k为任意样本。是第T-1次参数更新后第 l-1层DAE中间层的第j个神经元到第l层DAE中间层的第i个神经元的权值。是第T次参数更新后第l-1层DAE中间层的第j 个神经元到第l层DAE中间层的第i个神经元的权值。r和r+m分别是此批量的起始样本序号。为前T次权值迭代积累的梯度。为前T-1次权值迭代积累的梯度。Δ为增量符号。d为微分符号。σ为常数。本实施例取σ=10-6,用于被小数除时的数值稳定。为偏导符号。
式中,LH(Yt,Y)为根据顶层输出和训练样本输出 Y构造均方差损失函数。k为任意样本。是第T次参数更新后,第l层DAE中间层的第i个神经元的偏移量。是第T-1次参数更新后第l-1层DAE中间层的第j个神经元到第l层DAE中间层的第i个神经元的偏移量。为前T次偏置迭代积累的梯度;为前T-1次偏置迭代积累的梯度。
4.1.4)根据无监督预训练参数更新公式计算每层DAE的最优编码参数θ={W,b},并将最优编码参数θ={W,b}作为有监督微调的初始编码参数。
4.2)对SDAE最优潮流模型进行有监督微调,主要步骤如下:
argθminJ(W,b)=argθminLH(Yt,Y)。 (13)
4.2.2)根据优化目标函数argθminJ(W,b),从而对SDAE最优潮流模型最优编码参数θ={W,b}进行微调。
5)采用MCS法对待计算概率潮流的电力系统的随机变量进行抽样,从而获取计算样本。所述随机变量主要包括待计算概率最优潮流的电力系统的风速、光照辐射度和负荷。
6)将步骤5)得到的训练样本数据一次性输入步骤4)中训练完成的SDAE最优潮流模型中,从而计算出最优潮流在线概率。
7)对所述最优潮流在线概率进行分析,即绘制SDAE最优潮流模型的输出变量的概率密度曲线。
实施例2:
参见图3的(a)图至图3的(d)图,一种基于堆栈降噪自动编码器的电力系统概率最优潮流计算方法的仿真实验,主要包括以下步骤:
1)最优潮流样本获取。
本实施例中采用IEEE118标准系统进行仿真。本实施例在母线 59、80和90上引入风电场,风电场最大出力分别为220、200、260 MW,在母线13、14、16和23上引入光伏发电站,光伏发电站最大出力分别为100、150、100、150MW。
其中,假定风速服从两参数Weibull分布,其尺度参数为2.016,形状参数为5.089。光照强度服从Beta分布,光伏发电站的形状参数和风电场的切入风速、额定风速、切出风速参见表1。此外,假定各节点负荷的随机特性服从正态分布,其标准差为各节点负荷期望值的5%。
表1光伏发电站和风电场相关参数
然后,通过蒙特卡洛法对上述随机变量进行5万次抽样并采用内点法对每个抽样状态进行最优潮流的求解。将所有抽样状态的新能源节点和负荷节点的有功和无功作为训练样本输入X。由内点法求得的所有抽样状态的最优潮流计算结果(即电力系统各节点电压幅值和相角、发电机有功无功出力、各支路功率、发电成本)作为训练样本的输出Y。
2)SDAE最优潮流模型初始化
包括数据预处理以及SDAE最优潮流模型超参数的确定。使用最大最小法对训练样本的输入和输出数据进行归一化处理。根据训练样本容量,将训练样本分成100个批量。腐蚀训练样本输入X。最后,根据待求解的电力系统规模和复杂程度,设定SDAE最优潮流模型的层数l为5、每层神经元个数分别为236、200、200、200、717,学习率η为0.001。
3)SDAE最优潮流模型无监督预训练。
首先,利用训练样本输入X,构建第一层DAE的均方差损失函数 LH(Z1,X1)。然后,使用基于动量学习率的小批量梯度下降算法,根据参数更新公式,迭代求解第一层DAE的最优参数W1、b1、W′1、b′1。之后,得到第一层DAE的中间层输出,作为第二层DAE的输入,以相同的方法构建第二层DAE的损失函数LH(Z2,X2),并以相同的方法更新参数,以此类推,自底至顶逐层求解每层DAE的最优编码参数θ={W,b},并将该最优编码参数作为下阶段有监督微调的初始参数。
4)SDAE最优潮流模型有监督微调
首先,利用训练样本输入X和输出Y,构建SDAE最优潮流模型的均方差损失函数LH(Yt,Y)。然后,仍然使用基于动量学习率的小批量梯度下降算法,同样根据参数更新公式,迭代求解SDAE最优潮流模型的所有最优编码参数θ={W,b}。至此,完成了SDAE最优潮流模型的训练。
5)系统状态抽样
本实施例采用MCS法对所研究电力系统中的不确定因素,即风速、光照辐射强度、负荷等随机变量按其各自分布进行抽样以获取足够数量的测试样本,本发明选择的抽样次数N为50000。也可采用改进MCS 法。
6)概率最优潮流在线求解
将步骤(5)中得到的测试样本一次性输入步骤(4)中训练完成的SDAE最优模型中,该模型可直接映射出所有测试样本的最优潮流值。
仿真结果如下:
1)具体算例及最优潮流计算对比方法
具体算例为:IEEE 118系统,各节点负荷方差均为5%。
仿真中最优潮流计算对比方法包括M0-M1:
M0:基于内点法的MCS法,作为验证标准。
M1:基于堆栈降噪自动编码器的电力系统概率最优潮流计算方法。
在训练上述神经网络时,训练结束判据为连续10次迭代中均方差的最小值没有减小。
7)SDAE最优流模型计算精度分析
本实施例为了验证SDAE最优潮流模型计算最优潮流的总体精度及泛化能力,由算例经MCS法抽取50000个测试样本,并由两种对比算法计算所有样本的最优潮流,表2列出了由两种算法求得的节点电压幅值误差超过0.01p.u.的概率、发电机有功出力误差超过5MW的概率、支路有功误差超过5MW的概率、发电成本误差超过3000元的概率。
表2 M0-M1基于绝对误差的最优潮流计算精度比较
由表2可知,由SDAE最优潮流模型计算所得的四个指标中,绝对误差大于设置数值的概率集中在0.1%以下,最大概率为0.09%,可见, SDAE凭借其深层堆栈结构与编码解码过程,有效地提取了最优潮流的非线性特征,实现了最优潮流由输入到输出的高精度快速映射。因此,本发明构建的SDAE最优潮流模型具有较高的最优潮流计算精度和较强的泛化能力。
8)概率最优潮流在线算法计算性能分析
本实施例从所提概率最优潮流在线算法的计算精度和速度两方面分析其性能。在抽取的50000个样本的条件下,本实施例以节点10 的电压幅值、节点10的发电机有功出力、节点5-节点6的支路有功和发电成本为例,绘制出本发明方法与基于内点法的MCS法求得的所列随机变量的概率密度曲线,如图3的(a)图至图3的(d)图所示。从图中可见,本发明方法求得的最优潮流概率密度曲线与作为验证标准的基于内点法的 MCS法得到的最优潮流概率密度曲线几乎相同。因此,本发明方法满足概率最优潮流的计算精度要求。
表3列出了由M0、M1法计算概率最优潮流所需的时间。由表3可见,由M0法计算概率最优潮流耗时7582.6秒,而M1法仅需5.046秒,计算速度提高了上千倍。可见,采用本发明方法计算概率最优潮流可大幅度提高计算速度,能够满足概率最优潮流在线计算的速度要求。
表3 M0-M1计算概率最优潮流的时间对比
方法 | M0 | M1 |
时间(秒) | 7582.6 | 5.046 |
从实验结果可知:本发明所提出的基于SDAE并结合MCS法的概率最优潮流在线算法,能够以非常高的精度逼近高维非线性最优潮流模型,从而实现了最优潮流的非迭代计算,具有高计算精度,并且对于训练中未涉及的测试样本也能直接映射出最优潮流值,其计算结果及最优潮流的概率密度曲线与基于内点法的MCS法的计算结果良好吻合,具有很强的泛化能力,同时较内点法大幅度减少了计算时间,从而实现了概率最优潮流的高精度在线计算。
综上所述,本发明提出了一种基于SDAE并结合MCS法的概率最优潮流快速求解算法,克服了现有算法求解速度慢的问题,能够在不增加硬件成本的前提下,大大提高概率最优潮流的计算速度。此外,本发明引入了最大最小归一化法以及基于动量学习率的小批量梯度下降法,提升了SDAE最优潮流模型的训练精度和速度。通过算例仿真分析,验证了本发明所提方法的有效性。因此,本发明可为电力系统概率最优潮流的高精度在线计算提供技术支撑。
Claims (1)
1.基于堆栈降噪自动编码器的电力系统概率最优潮流计算方法,其特征在于,包括以下步骤:
1)建立SDAE最优潮流模型;
2)获取SDAE最优潮流模型输入层的输入样本X;
3)对SDAE最优潮流模型进行初始化;
4)对SDAE最优潮流模型进行训练,从而得到训练后的SDAE最优潮流模型;
5)采用MCS法对待计算概率潮流的电力系统的随机变量进行抽样,从而获取计算样本;
6)将步骤5)得到的训练样本数据一次性输入步骤4)中训练完成的SDAE最优潮流模型中,从而计算出最优潮流在线概率;
7)对所述最优潮流在线概率进行分析,即绘制SDAE最优潮流模型的输出变量的概率密度曲线;
所述随机变量主要包括待计算概率最优潮流的电力系统的风速、光照辐射度和负荷;
所述SDAE最优潮流模型包括n个顺序堆叠的DAE模型;其中,第l-1层DAE的输入层为第l-2层DAE的中间层;第l-1层DAE的中间层为第l层DAE的输入层;
所述DAE模型主要包括1个输入层、1个中间层和1个输出层;
输入层的神经元个数设为A;任意输入层神经元记为uI,I=1、2…A;
输入层的输入样本X中的数据主要包括电力系统中所有新能源节点和负荷节点的有功功率和无功功率;
中间层的神经元个数设为M;任意输入层神经元记为vi,i=1、2…M;
输出层的神经元个数设为Q;任意输出层神经元记为me,e=1、2…Q;
输出层的输出向量y中的数据主要包括电力系统节点电压幅值和相角、发电机有功功率、发电机有功功率无功功率、支路有功功率、支路无功功率和发电成本;
建立SDAE最优潮流模型的步骤如下:
式中,qD为随机局部腐蚀过程,即随机选取若干输入变量置零;
II)利用编码函数fθ得到中间层a;中间层a如下所示:
激活函数s如下所示:
s(x)=max(0,x); (3)
式中,x为输入向量X中的数据;
III)利用解码函数gθ′得到DAE的输出层Z;输出层Z如下所示:
Z=gθ′(a)=s(W′a+b′); (4)
式中,W′为解码函数的权值;W′=WT;b′为解码函数的偏置;s为激活函数;gθ′为解码函数;
IV)将DAE逐层堆叠,从而得到SDAE最优潮流模型;
V)提取最优潮流输入数据X的高维特征,拟合得到输出Yt;输出Yt如下所示:
对SDAE最优潮流模型进行初始化的步骤如下:
a)利用最大最小法对输入向量X和输出向量y进行归一化处理
式中,x为输入向量X中的数据;xmin为输入向量X中的最小数据;xmax为输入向量X中的最大数据;
利用最大最小法对输出向量y进行归一化处理;
式中,y为输出向量;ymin为最小输出向量;ymax为最大输出向量;
b)将处理后的训练输入样本X分成m个批量;
c)设定SDAE最优潮流模型的总层数n、每层神经元的个数和学习率η;
对SDAE最优潮流模型进行训练的步骤如下:
s1)对SDAE最优潮流模型进行无监督预训练,步骤如下:
s1.1)根据第l层DAE的输入Xl和输出Zl构造均方差损失函数LH(Xl,Zl);
式中,Xl为第l层DAE的输入,也即第l-1层DAE的中间层输出al-1;Zl是第l层DAE的输出;LH(Xl,Zl)为根据第l层DAE的输入Xl和输出Zl构造的均方差损失函数;θl为第l层DAE的编码参数;
s1.3)利用RMSProp学习算法和动量学习率,构建SDAE最优潮流模型无监督预训练参数更新公式;无监督预训练参数更新公式如公式(9)至公式(12)所示;
式中,是第T次参数更新后第l-1层DAE中间层的第j个神经元到第l层DAE中间层的第i个神经元的权值;是第T-1次参数更新后第l-1层DAE中间层的第j个神经元到第l层DAE中间层的第i个神经元的权值;
式中,η为学习速率;m为此批量样本数量;p为动量因子;ρ为梯度累积指数;ε为常数;⊙是Hadamard乘子;LH(Zl,Xl)为根据第l个DAE的输出层输出和训练样本输入X构造均方差损失函数;k为任意样本;是第T-1次参数更新后第l-1层DAE中间层的第j个神经元到第l层DAE中间层的第i个神经元的权值;是第T次参数更新后第l-1层DAE中间层的第j个神经元到第l层DAE中间层的第i个神经元的权值;r和r+m分别是此批量的起始样本序号;为前T次权值迭代积累的梯度;为前T-1次权值迭代积累的梯度;Δ为增量符号;d为微分符号;σ为常数;为偏导符号;
式中,LH(Yt,Y)为根据顶层输出和训练样本输出Y构造均方差损失函数;k为任意样本;是第T次参数更新后,第l层DAE中间层的第i个神经元的偏移量;是第T-1次参数更新后第l-1层DAE中间层的第j个神经元到第l层DAE中间层的第i个神经元的偏移量;为前T次偏置迭代积累的梯度;为前T-1次偏置迭代积累的梯度;
s1.4)根据无监督预训练参数更新公式计算每层DAE的最优编码参数θ={W,b},并将最优编码参数θ={W,b}作为有监督微调的初始编码参数;
s2)对SDAE最优潮流模型进行有监督微调,步骤如下:
argθminJ(W,b)=argθminLH(Yt,Y); (13)
s2.2)根据优化目标函数argθminJ(W,b),从而对SDAE最优潮流模型最优编码参数θ={W,b}进行微调;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811633643.4A CN109599872B (zh) | 2018-12-29 | 2018-12-29 | 基于堆栈降噪自动编码器的电力系统概率潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811633643.4A CN109599872B (zh) | 2018-12-29 | 2018-12-29 | 基于堆栈降噪自动编码器的电力系统概率潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109599872A CN109599872A (zh) | 2019-04-09 |
CN109599872B true CN109599872B (zh) | 2022-11-08 |
Family
ID=65964799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811633643.4A Active CN109599872B (zh) | 2018-12-29 | 2018-12-29 | 基于堆栈降噪自动编码器的电力系统概率潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109599872B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110009529B (zh) * | 2019-04-15 | 2021-04-16 | 湖南大学 | 一种基于堆栈降噪自动编码器的暂态频率获取方法 |
CN110829434B (zh) * | 2019-09-30 | 2021-04-06 | 重庆大学 | 一种提高深度神经网络潮流模型扩展性的方法 |
CN110929989B (zh) * | 2019-10-29 | 2023-04-18 | 重庆大学 | 基于深度学习的含不确定性n-1安全校核方法 |
CN111325233B (zh) * | 2019-10-29 | 2024-02-09 | 国网辽宁省电力有限公司电力科学研究院 | 一种变压器故障检测方法及装置 |
CN111796173B (zh) * | 2020-08-13 | 2022-01-21 | 广东电网有限责任公司 | 局部放电模式识别方法、计算机设备和存储介质 |
CN116154779B (zh) * | 2023-04-23 | 2023-07-18 | 之江实验室 | 一种基于预训练模型的最优潮流计算方法和装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106443447A (zh) * | 2016-09-26 | 2017-02-22 | 南京航空航天大学 | 一种基于iSDAE的航空发电机故障特征提取方法 |
CN108304623A (zh) * | 2018-01-15 | 2018-07-20 | 重庆大学 | 一种基于堆栈降噪自动编码器的概率潮流在线计算方法 |
CN108734391A (zh) * | 2018-05-08 | 2018-11-02 | 重庆大学 | 基于堆栈降噪自动编码器的电-气综合能源系统概率能流量计算方法 |
-
2018
- 2018-12-29 CN CN201811633643.4A patent/CN109599872B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106443447A (zh) * | 2016-09-26 | 2017-02-22 | 南京航空航天大学 | 一种基于iSDAE的航空发电机故障特征提取方法 |
CN108304623A (zh) * | 2018-01-15 | 2018-07-20 | 重庆大学 | 一种基于堆栈降噪自动编码器的概率潮流在线计算方法 |
CN108734391A (zh) * | 2018-05-08 | 2018-11-02 | 重庆大学 | 基于堆栈降噪自动编码器的电-气综合能源系统概率能流量计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109599872A (zh) | 2019-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109599872B (zh) | 基于堆栈降噪自动编码器的电力系统概率潮流计算方法 | |
CN109117951B (zh) | 基于bp神经网络的概率潮流在线计算方法 | |
CN108304623B (zh) | 一种基于堆栈降噪自动编码器的概率潮流在线计算方法 | |
CN108734391B (zh) | 基于堆栈降噪自动编码器的电-气综合能源系统概率能流量计算方法 | |
CN111008728A (zh) | 一种用于分布式光伏发电系统短期出力的预测方法 | |
CN108711847B (zh) | 一种基于编码解码长短期记忆网络的短期风电功率预测方法 | |
CN110110434B (zh) | 一种概率潮流深度神经网络计算的初始化方法 | |
CN109978284B (zh) | 一种基于混合神经网络模型的光伏发电功率分时预测方法 | |
CN109412161B (zh) | 一种电力系统概率潮流计算方法及系统 | |
CN108336739B (zh) | 一种基于rbf神经网络的概率潮流在线计算方法 | |
CN109523084A (zh) | 一种基于主元分析和机器学习的风电功率超短期预测方法 | |
CN109995031B (zh) | 基于物理模型的概率潮流深度学习计算方法 | |
CN115275991A (zh) | 一种基于iemd-ta-lstm模型的主动配电网运行态势预测方法 | |
CN110795841B (zh) | 一种间歇性能源出力不确定性的数学建模方法 | |
Chen et al. | Research on wind power prediction method based on convolutional neural network and genetic algorithm | |
CN115169742A (zh) | 一种短期风力发电功率预测方法 | |
CN116629416A (zh) | 光伏电站功率预测方法及装置 | |
Xing et al. | Comparative study of reformed neural network based short‐term wind power forecasting models | |
CN116702937A (zh) | 一种基于K-means均值聚类及优化BP神经网络的光伏出力日前预测方法 | |
Xu et al. | Short-term wind speed prediction based on GRU | |
CN112418504B (zh) | 一种基于混合变量选择优化深度信念网络风速预测方法 | |
Phan et al. | Application of a new Transformer-based model and XGBoost to improve one-day-ahead solar power forecasts | |
Wang et al. | A prediction model for ultra-short-term output power of wind farms based on deep learning | |
CN116559975A (zh) | 一种基于多元时间序列图神经网络的多步长气象预测方法 | |
CN110829434B (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 |