CN109993346A - 基于混沌时间序列和神经网络的微电网电压安全评估方法 - Google Patents
基于混沌时间序列和神经网络的微电网电压安全评估方法 Download PDFInfo
- Publication number
- CN109993346A CN109993346A CN201910131795.2A CN201910131795A CN109993346A CN 109993346 A CN109993346 A CN 109993346A CN 201910131795 A CN201910131795 A CN 201910131795A CN 109993346 A CN109993346 A CN 109993346A
- Authority
- CN
- China
- Prior art keywords
- voltage
- micro
- capacitance sensor
- network
- value
- 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
- 238000011156 evaluation Methods 0.000 title claims abstract description 31
- 238000013528 artificial neural network Methods 0.000 title claims abstract description 23
- 238000012549 training Methods 0.000 claims abstract description 32
- 230000001537 neural effect Effects 0.000 claims abstract description 13
- 230000000739 chaotic effect Effects 0.000 claims abstract description 5
- 230000001186 cumulative effect Effects 0.000 claims description 6
- 238000005315 distribution function Methods 0.000 claims description 6
- 210000002569 neuron Anatomy 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 238000005314 correlation function Methods 0.000 claims description 3
- 238000005183 dynamical system Methods 0.000 claims description 3
- 230000005284 excitation Effects 0.000 claims description 3
- 230000002452 interceptive effect Effects 0.000 claims description 3
- 238000012886 linear function Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 230000002035 prolonged effect Effects 0.000 claims description 3
- 230000004927 fusion Effects 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 description 7
- 238000000034 method Methods 0.000 description 6
- 230000005611 electricity Effects 0.000 description 4
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 3
- 229910052799 carbon Inorganic materials 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003912 environmental pollution Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000001568 sexual effect Effects 0.000 description 2
- 238000012956 testing procedure Methods 0.000 description 2
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
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/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
-
- 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
- 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)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- General Business, Economics & Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Molecular Biology (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Game Theory and Decision Science (AREA)
- Educational Administration (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Power Engineering (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提出了基于混沌时间序列和神经网络的微电网电压安全评估方法,包括以下步骤:步骤1,输入电压数据时间序列;步骤2,进行安全评估;步骤3,对电压数据的混沌性进行判断;步骤4,进行相空间重构;步骤5,RBF神经网络的训练;步骤6,提供判断微电网电压安全可控依据;步骤7,结束;将微电网的公共连接点(PCC)的电压作为评估微电网是否安全的依据,采用电压安全评估指数(VSAI)算法,将采集到的微电网PCC的电压序列移动平均值与采集点的电压作比较,从而得到偏差值,再带入所定义的VSAI算法中,确定微电网的电压是否处于安全状态,最后根据混沌时间序列和RBF神经网络算法的融合,预测出短时电压变化曲线,为微电网的安全可控提供依据。
Description
技术领域
本发明涉及一种电压安全评估方法,具体的说是一种微电网电压安全评估方法,属于电力系统自动化技术领域。
背景技术
随着常规化石能源的不断消耗以及引起的碳排放、环境污染等问题日趋严重,寻求清洁高效的可再生能源以替代可耗竭的化石能源已受到人们普遍关注。利用风电、光伏等可再生资源替代化石能源,对于降低碳排放,改善能源结构将起到重要的作用,但由于可再生资源具有明显的间歇性及波动性,使得包含可再生能源的微电网调控难度不断增加,而近几年,以降低碳排放、节约常规能源和减少环境污染为目的电动汽车接入微电网进行自由充电,都会使电源侧和负荷侧呈现出一定的非可控性,从而给微电网的安全和调控带来了新的挑战。
传统的微电网的电压预测和安全评估的研究大都基于物理或者统计的方法,其中统计方法主要基于单一的预测模型,有一定的主观性,预测误差较大,因此针对包含波动性和间歇性的分布式电源以及随机性负荷的微电网电压进行安全评估的问题,急需一种更好的方法来解决。
发明内容
本发明所要解决的技术问题是,预测模型单一,有一定的主观性,预测误差较大,使得微电网的安全评估精度差,不可靠,为克服现有技术的不足而提供基于混沌时间序列和神经网络的微电网电压安全评估方法。
本发明提供基于混沌时间序列和神经网络的微电网电压安全评估方法,包括以下步骤:
步骤1,输入电压数据时间序列;
步骤2,进行安全评估;
步骤3,对电压数据的混沌性进行判断;
步骤4,进行相空间重构;
步骤5,RBF神经网络的训练;
步骤6,提供判断微电网电压安全可控依据;
步骤7,结束。
进一步的,所述关键字层次结构分为测试用例集合、测试用例层、测试步骤层,关键字的定义与测试用例集合的结构密不可分,一套测试用例集合包含了一系列测试用例,而一个测试用例是由各个测试步骤组合构成的。
进一步的,所述步骤2中,采用VSAI的动态电压稳定性标准进行安全评估,根据负载总线电压的均方根值的时间序列数据进行计算,将采集到的离散电压序列{v,i 1,2,n},由电压安全评估指数按照以下步骤配置:
步骤2.1,求取重构后电压数据的移动平均值:
步骤2.2,计算测量点的电压vi与电压的移动平局值之间的偏差ci:
步骤2.3,在第i个瞬间将百分比多样性曲线下的面积除以N的值如下:
步骤2.4,第i个瞬间的VSAI定义为:
VSAIi=uth-ui (4)
其中uth是安全阈值;
动态电压安全评估的标准如下:当VSAIi=uth-ui≤0那么说明微电网系统总线电压被评估为危险,而只有当VSAIi=uth-ui>0时,微电网电压才处于一个安全状态。
进一步的,所述步骤3的具体内容如下:
建立一维的动力系统:xn+1=f(xn)n=0,1,2,… (5)
其中初值x0受到一个微小的扰动δx0,通过一次迭代后,偏差为
经过n次迭代后,绝对偏差为:
其中:
定义:
为Lyapunov指数;
当λ为正值时,表示系统相邻轨道迅速分离,长时间运动局部不稳定且对初始值敏感,轨线在有界区域内多次折叠,从而导致混沌。
进一步的,所述步骤4的具体内容如下:
长度为n的电压时间序列{vi,i=1,2,…n},其延迟时间和嵌入维数分别为τ和m,根据嵌入定理可知,重构的相空间为:
V=[V1,V2,…,VM]T (11)
其中,Vi=[vi,vi+τ,vi+2τ,…,vi+(m-1)τ]
即:[V1]=[v1 v1+τ v1+2τ…v1+(m-1)τ] (12)
式中,M为重构后相空间中的相点数,满足M=n-(m-1)τ,Vi为第i个相点。
进一步的,所述延迟时间τ采用互相信息法求取,以两个随机变量间一般性度量,决定了序列vi相对于其它某个时刻的vi+τ具有多少信息,其中vi和vi+τ的交互信息形式为:
式中,P(vi)和P(vi+τ)为vi和vi+τ的归一化分布,P(vi,vi+τ)为其联合分布;
混沌时间序列的交互信息参数I(τ)总是大于0,选取参数I(τ)的第一个最小值作为时间序列的延迟时间。
进一步的,所述嵌入维数m采用G-P算法进行求取,具体步骤如下:
步骤4.2.1,利用时间序列v1,v2,…,vn-1,vn,…,先给一个小值m0,对应一个重构的相空间;
步骤4.2.2,计算关联函数
其中|Y(ti)-Y(tj)|表示相点Y(ti)和Y(tj)之间的距离,θ(z)是Heaviside函数,C(r)是一个累积分布函数,表示相空间吸引子上两点之间距离小于r的概率;
步骤4.2.3,对于r的适当范围,吸引子的维数d与累积分布函数C(r)应满足对数线性关系,即d(m)=lnC(r)/lnr,从而由拟合求出对应于m0的关联维数估计值d(m0);
步骤4.2.4,增加嵌入维数m1>m0,重复计算步骤4.2.2和4.2.3,直到相应的维数估计值d(m)不在随m的增长而在误差范围内不变为止。此时得到的d即为吸引子的关联维数。
进一步的,所述步骤5的具体内容如下:
RBF神经网络的包括输入层、隐含层和输出层,为多层前向网络结构,所述输入层的信号源节点可传递输入信号到隐含层得节点;隐含层的节点呈辐射状函数结构;输出层的传递函数是线性函数,是对输入模式的作用做出响应;将RBF神经网络的输入层到隐含层的权值设定为常数1,网络训练时训练隐含层到输出层的权值,包括以下内容:
在网络学习中,样本输入数据为X1,X2,X3,…,XN,样本输入数据对应的目标输出为y1,y2,y3,…,yN,输入数据的样本表示为{Xi,yi},i=1,2,…,N,网络不考虑阈值的输出为F(Xj),隐含层第i个神经元的激励函数为ui(·);
网络学习时的目标函数定义为:
式中,βj为遗忘因子;
网络输出误差值为:
网络输出F(X)分别对网络中心ci输出权值wi的偏导(梯度)为:
网络中心ci,输出权值wi的修正量为:
式中,ui(Xj)为第i个隐含层神经元对Xj的输出,η为网络学习率。
网络在训练时,将求出的权值修正值加上前一次训练的权值就可以得到修正后的权值。
进一步的,所述步骤6的包括以下具体步骤:
步骤6.1,以欧式距离作为样本的选择标准,选取相空间中离第M个相点的欧式距离最近的k个相点为训练样本的输入部分,训练样本中的目标输出为所选的相点在相空间对应位置的下一行的最后一个分量,由此便可以组成网络的训练样本对。将训练样本对输入网络进行训练,当网络达到训练精度或达到最大学习次数时,停止网格训练。
步骤6.2,在得到满足要求的网络后,以相空间中第M个相点作为网格的输入,此时便可以得到网络的第一个预测结果,,将第一个预测结果反馈到网络的输入端,以此来重新构造网络的输入,然后再进行预测得到第二个预测结果,直至得到预测结果。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:根据微电网实际运行特点,考虑到运行过程中由于光伏、风力等分布式能源的随机性和波动性而导致分布式电源出力的不稳定以及用户侧需求的动态变化的问题,将微电网的公共连接点(PCC)的电压作为评估微电网是否安全的依据,采用电压安全评估指数(VSAI)算法,将采集到的微电网PCC的电压序列移动平均值与采集点的电压作比较,从而得到偏差值,再带入所定义的VSAI算法中,确定微电网的电压是否处于安全状态,最后根据混沌时间序列和RBF神经网络算法的融合,预测出短时电压变化曲线,为微电网的安全可控提供依据。
附图说明
图1为本发明的流程图。
图2是本发明中RBF神经网络结构图。
图3是本发明中直流微电网的结构图。
图4是本发明中60s内三个节点的电压曲线图。
图5是本发明中直流微电网第一节点处的电压曲线图。
图6是本发明中关键直流微电网的电压和VSAI图
图7是本发明中新相空间的延迟时间图。
图8是本发明中新相空间的嵌入维数图。
图9是本发明中实际值与预测值的对比图。
图10是本发明中预测绝对误差曲线图。
具体实施方式
下面结合附图1-10对本发明的技术方案做进一步的详细说明:
本实施例提出了基于混沌时间序列和神经网络的微电网电压安全评估方法,包括以下步骤:
步骤1,输入电压数据时间序列;
步骤2,进行安全评估,采用VSAI的动态电压稳定性标准进行安全评估,根据负载总线电压的均方根值的时间序列数据进行计算,将采集到的离散电压序列{v,i 1,2,n},由电压安全评估指数按照以下步骤配置:
步骤2.1,求取重构后电压数据的移动平均值:
步骤2.2,计算测量点的电压vi与电压的移动平局值之间的偏差ci:
步骤2.3,在第i个瞬间将百分比多样性曲线下的面积除以N的值如下:
步骤2.4,第i个瞬间的VSAI定义为:
VSAIi=uth-ui (4)
其中uth是安全阈值;
动态电压安全评估的标准如下:当VSAIi=uth-ui≤0那么说明微电网系统总线电压被评估为危险,而只有当VSAIi=uth-ui>0时,微电网电压才处于一个安全状态。
步骤3,对电压数据的混沌性进行判断,具体内容如下:
建立一维的动力系统:xn+1=f(xn)n=0,1,2,… (5)
其中初值x0受到一个微小的扰动δx0,通过一次迭代后,偏差为
经过n次迭代后,绝对偏差为:
其中:
定义:
为Lyapunov指数;
当λ为正值时,表示系统相邻轨道迅速分离,长时间运动局部不稳定且对初始值敏感,轨线在有界区域内多次折叠,从而导致混沌。
步骤4,进行相空间重构,具体内容如下:
长度为n的电压时间序列{vi,i=1,2,…n},其延迟时间和嵌入维数分别为τ和m,根据嵌入定理可知,重构的相空间为:
V=[V1,V2,…,VM]T (11)
其中,Vi=[vi,vi+τ,vi+2τ,…,vi+(m-1)τ]
即:[V1]=[v1 v1+τ v1+2τ…v1+(m-1)τ] (12)
式中,M为重构后相空间中的相点数,满足M=n-(m-1)τ,Vi为第i个相点。
所述延迟时间τ采用互相信息法求取,以两个随机变量间一般性度量,决定了序列vi相对于其它某个时刻的vi+τ具有多少信息,其中vi和vi+τ的交互信息形式为:
式中,P(vi)和P(vi+τ)为vi和vi+τ的归一化分布,P(vi,vi+τ)为其联合分布;
混沌时间序列的交互信息参数I(τ)总是大于0,选取参数I(τ)的第一个最小值作为时间序列的延迟时间。
所述嵌入维数m采用G-P算法进行求取,具体步骤如下:
步骤4.2.1,利用时间序列v1,v2,…,vn-1,vn,…,先给一个小值m0,对应一个重构的相空间;
步骤4.2.2,计算关联函数
其中|Y(ti)-Y(tj)|表示相点Y(ti)和Y(tj)之间的距离,θ(z)是Heaviside函数,C(r)是一个累积分布函数,表示相空间吸引子上两点之间距离小于r的概率;
步骤4.2.3,对于r的适当范围,吸引子的维数d与累积分布函数C(r)应满足对数线性关系,即d(m)=lnC(r)/lnr,从而由拟合求出对应于m0的关联维数估计值d(m0);
步骤4.2.4,增加嵌入维数m1>m0,重复计算步骤4.2.2和4.2.3,直到相应的维数估计值d(m)不在随m的增长而在误差范围内不变为止。此时得到的d即为吸引子的关联维数。
步骤5,RBF神经网络的训练,具体内容如下:
RBF神经网络的包括输入层、隐含层和输出层,为多层前向网络结构,所述输入层的信号源节点可传递输入信号到隐含层得节点;隐含层的节点呈辐射状函数结构;输出层的传递函数是线性函数,是对输入模式的作用做出响应;将RBF神经网络的输入层到隐含层的权值设定为常数1,网络训练时训练隐含层到输出层的权值,包括以下内容:
在网络学习中,样本输入数据为X1,X2,X3,…,XN,样本输入数据对应的目标输出为y1,y2,y3,…,yN,输入数据的样本表示为{Xi,yi},i=1,2,…,N,网络不考虑阈值的输出为F(Xj),隐含层第i个神经元的激励函数为ui(·);
网络学习时的目标函数定义为:
式中,βj为遗忘因子;
网络输出误差值为:
网络输出F(X)分别对网络中心ci输出权值wi的偏导(梯度)为:
网络中心ci,输出权值wi的修正量为:
式中,ui(Xj)为第i个隐含层神经元对Xj的输出,η为网络学习率。
网络在训练时,将求出的权值修正值加上前一次训练的权值就可以得到修正后的权值。
步骤6,提供判断微电网电压安全可控依据,包括以下具体步骤:
步骤6.1,以欧式距离作为样本的选择标准,选取相空间中离第M个相点的欧式距离最近的k个相点为训练样本的输入部分,训练样本中的目标输出为所选的相点在相空间对应位置的下一行的最后一个分量,由此便可以组成网络的训练样本对。将训练样本对输入网络进行训练,当网络达到训练精度或达到最大学习次数时,停止网格训练。
步骤6.2,在得到满足要求的网络后,以相空间中第M个相点作为网格的输入,此时便可以得到网络的第一个预测结果,,将第一个预测结果反馈到网络的输入端,以此来重新构造网络的输入,然后再进行预测得到第二个预测结果,直至得到预测结果。。
步骤7,结束。
本发明中,如图2所示,直流微电网由负载,蓄电系统,不同逆变器和不同类型的配电组组成。直流母线是主要的支柱网络,它将不同类型的直流设备与微电网直接连接,同时交流设备通过直流/交流逆变器连接微电网。为了保持微电网PCC的电压稳定性,减少谐波注入分布式系统,在微电网和分布式能量系统之间使用静态开关,也可以保持分布式网络与微电网之间的双向功率流。利用直流微电网PCC的动态电压数据,第一节点是分布式光伏发电与微电网之间的连接点,第二节点是负载与微电网之间的连接点,第三节点是能量存储系统和微电网之间的连接点。三个节点的电压图如图3所示,采样间隔时间为1秒,持续1分钟(每个节点约60个实验数据)。
在图3中,通过直流微电网中三点的动态电压数据,利用VSAI评估系统的安全性,并通过混沌时间序列理论和RBF神经网络在短时间内预测动态电压。
A.证明VSAI的有效性
由于分布式光伏能量系统的随机性和不可控性,光伏电源和微电网的连接点处的电压通常是不稳定的。当第1节点的电压失效时,直流微电网第1点的电压曲线如图4所示。在t=12s时,电压开始下降并且电压开始稳定直到t=47s。为了确认VSAI在评估动态电压安全性方面的有效性,对DC微电网进行了测试。如图5所示,当电压在t=12s开始下降时,IEEE-30系统第1节点的相应VSAI也从0开始下降。当第1节点的故障停止,第1节点的电压在另一种安全情况下,VSAI也返回0.因此,可以根据VSAI的值确定节点电压是否处于安全状态。
B.确定混乱的存在
根据文章前面提出的小数据方法的理论,动态电压数据的最大Lyapunov指数λ=1.3972是通过Matlab编程的计算得出的。因为λ>0,可以确定动态电压序列的混沌特征。
C.重建相空间
重建动态电压序列的新相空间:
通过B部分所述的互信息的方式计算延迟时间τ=1,并清楚地在图6中指出。通过GP算法计算嵌入维数m=9,如图7所示。
D.混沌理论与神经网络预测
神经网络输入端口的数量可以通过计算嵌入维数m=9来确定,同时根据前三个节点的采样时间确定预测步长为60。直流微电网拥有8个采样点,每个点有60个采样数据,因此将480个采样点作为重建相空间的数据。重建后,获得300个相位点,覆盖180个训练样本和120个测试样本。
在RBF神经网络中,输入层神经元的数量为9,由重建的动态电压时间序列的嵌入维数确定,输出为1.经过多次尝试,确定SPREAD参数为1,并使用神经网络函数newrb与MATLAB一起创建神经网络预测模型。由于空间限制,图8中只显示了120秒的采样点,图8中明显发现了预测与实际之间的差异。根据图9,最大绝对预测误差为0.0021,另一个在0~0.3%之间。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内,因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (8)
1.基于混沌时间序列和神经网络的微电网电压安全评估方法,其特征在于,包括以下步骤:
步骤1,输入电压数据时间序列;
步骤2,进行安全评估;
步骤3,对电压数据的混沌性进行判断;
步骤4,进行相空间重构;
步骤5,RBF神经网络的训练;
步骤6,提供判断微电网电压安全可控依据;
步骤7,结束。
2.根据权利要求1所述的基于混沌时间序列和神经网络的微电网电压安全评估方法,其特征在于,所述步骤2中,采用VSAI的动态电压稳定性标准进行安全评估,根据负载总线电压的均方根值的时间序列数据进行计算,将采集到的离散电压序列{v,i 1,2,n},由电压安全评估指数按照以下步骤配置:
步骤2.1,求取重构后电压数据的移动平均值:
步骤2.2,计算测量点的电压vi与电压的移动平局值之间的偏差ci:
步骤2.3,在第i个瞬间将百分比多样性曲线下的面积除以N的值如下:
步骤2.4,第i个瞬间的VSAI定义为:
VSAIi=uth-ui (4)
其中uth是安全阈值;
动态电压安全评估的标准如下:当VSAIi=uth-ui≤0那么说明微电网系统总线电压被评估为危险,而只有当VSAIi=uth-ui>0时,微电网电压才处于一个安全状态。
3.根据权利要求2所述的基于混沌时间序列和神经网络的微电网电压安全评估方法,其特征在于,所述步骤3的具体内容如下:
建立一维的动力系统:xn+1=f(xn)n=0,1,2,… (5)
其中初值x0受到一个微小的扰动δx0,通过一次迭代后,偏差为
δx1=f(x0+δx0)-f(x0)
≈f'(x0)δx0 (6)
经过n次迭代后,绝对偏差为:
其中:
定义:
为Lyapunov指数;
当λ为正值时,表示系统相邻轨道迅速分离,长时间运动局部不稳定且对初始值敏感,轨线在有界区域内多次折叠,从而导致混沌。
4.根据权利要求3所述的基于混沌时间序列和神经网络的微电网电压安全评估方法,其特征在于,所述步骤4的具体内容如下:
长度为n的电压时间序列{vi,i=1,2,…n},其延迟时间和嵌入维数分别为τ和m,根据嵌入定理可知,重构的相空间为:
V=[V1,V2,…,VM]T (11)
其中,Vi=[vi,vi+τ,vi+2τ,…,vi+(m-1)τ]
即:[V1]=[v1 v1+τ v1+2τ … v1+(m-1)τ] (12)
式中,M为重构后相空间中的相点数,满足M=n-(m-1)τ,Vi为第i个相点。
5.根据权利要求4所述的基于混沌时间序列和神经网络的微电网电压安全评估方法,其特征在于,所述延迟时间τ采用互相信息法求取,以两个随机变量间一般性度量,决定了序列vi相对于其它某个时刻的vi+τ具有多少信息,其中vi和vi+τ的交互信息形式为:
式中,P(vi)和P(vi+τ)为vi和vi+τ的归一化分布,P(vi,vi+τ)为其联合分布;
混沌时间序列的交互信息参数I(τ)总是大于0,选取参数I(τ)的第一个最小值作为时间序列的延迟时间。
6.根据权利要求5所述的基于混沌时间序列和神经网络的微电网电压安全评估方法,其特征在于,所述嵌入维数m采用G-P算法进行求取,具体步骤如下:
步骤4.2.1,利用时间序列v1,v2,…,vn-1,vn,…,先给一个小值m0,对应一个重构的相空间;
步骤4.2.2,计算关联函数
其中|Y(ti)-Y(tj)|表示相点Y(ti)和Y(tj)之间的距离,θ(z)是Heaviside函数,C(r)是一个累积分布函数,表示相空间吸引子上两点之间距离小于r的概率;
步骤4.2.3,对于r的适当范围,吸引子的维数d与累积分布函数C(r)应满足对数线性关系,即d(m)=lnC(r)/lnr,从而由拟合求出对应于m0的关联维数估计值d(m0);
步骤4.2.4,增加嵌入维数m1>m0,重复计算步骤4.2.2和4.2.3,直到相应的维数估计值d(m)不在随m的增长而在误差范围内不变为止。此时得到的d即为吸引子的关联维数。
7.根据权利要求6所述的基于混沌时间序列和神经网络的微电网电压安全评估方法,其特征在于,所述步骤5的具体内容如下:
RBF神经网络的包括输入层、隐含层和输出层,为多层前向网络结构,所述输入层的信号源节点可传递输入信号到隐含层得节点;隐含层的节点呈辐射状函数结构;输出层的传递函数是线性函数,是对输入模式的作用做出响应;将RBF神经网络的输入层到隐含层的权值设定为常数1,网络训练时训练隐含层到输出层的权值,包括以下内容:
在网络学习中,样本输入数据为X1,X2,X3,…,XN,样本输入数据对应的目标输出为y1,y2,y3,…,yN,输入数据的样本表示为{Xi,yi},i=1,2,…,N,网络不考虑阈值的输出为F(Xj),隐含层第i个神经元的激励函数为ui(·);
网络学习时的目标函数定义为:
式中,βj为遗忘因子;
网络输出误差值为:
网络输出F(X)分别对网络中心ci输出权值wi的偏导(梯度)为:
网络中心ci,输出权值wi的修正量为:
式中,ui(Xj)为第i个隐含层神经元对Xj的输出,η为网络学习率。
网络在训练时,将求出的权值修正值加上前一次训练的权值就可以得到修正后的权值。
8.根据权利要求7所述的基于混沌时间序列和神经网络的微电网电压安全评估方法,其特征在于,所述步骤6的包括以下具体步骤:
步骤6.1,以欧式距离作为样本的选择标准,选取相空间中离第M个相点的欧式距离最近的k个相点为训练样本的输入部分,训练样本中的目标输出为所选的相点在相空间对应位置的下一行的最后一个分量,由此便可以组成网络的训练样本对。将训练样本对输入网络进行训练,当网络达到训练精度或达到最大学习次数时,停止网格训练。
步骤6.2,在得到满足要求的网络后,以相空间中第M个相点作为网格的输入,此时便可以得到网络的第一个预测结果,,将第一个预测结果反馈到网络的输入端,以此来重新构造网络的输入,然后再进行预测得到第二个预测结果,直至得到预测结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910131795.2A CN109993346B (zh) | 2019-02-22 | 2019-02-22 | 基于混沌时间序列和神经网络的微电网电压安全评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910131795.2A CN109993346B (zh) | 2019-02-22 | 2019-02-22 | 基于混沌时间序列和神经网络的微电网电压安全评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109993346A true CN109993346A (zh) | 2019-07-09 |
CN109993346B CN109993346B (zh) | 2020-09-11 |
Family
ID=67130263
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910131795.2A Active CN109993346B (zh) | 2019-02-22 | 2019-02-22 | 基于混沌时间序列和神经网络的微电网电压安全评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109993346B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110555469A (zh) * | 2019-08-15 | 2019-12-10 | 阿里巴巴集团控股有限公司 | 处理交互序列数据的方法及装置 |
US10936950B1 (en) | 2019-08-15 | 2021-03-02 | Advanced New Technologies Co., Ltd. | Processing sequential interaction data |
CN112784887A (zh) * | 2021-01-12 | 2021-05-11 | 北京林业大学 | 一种植物生命状态评价方法及装置 |
CN114844227A (zh) * | 2022-07-04 | 2022-08-02 | 广东电网有限责任公司佛山供电局 | 一种电网运行安全监督管理系统及方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050105760A1 (en) * | 2001-12-14 | 2005-05-19 | Eggers Joachim J. | Data embedding and extraction |
CN103217900A (zh) * | 2013-02-06 | 2013-07-24 | 浙江工业大学 | 基于实时电价的中压微网混沌pso最优潮流实现方法 |
CN103440370A (zh) * | 2013-08-21 | 2013-12-11 | 国家电网公司 | 一种输变电工程造价评估方法及装置 |
-
2019
- 2019-02-22 CN CN201910131795.2A patent/CN109993346B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050105760A1 (en) * | 2001-12-14 | 2005-05-19 | Eggers Joachim J. | Data embedding and extraction |
CN103217900A (zh) * | 2013-02-06 | 2013-07-24 | 浙江工业大学 | 基于实时电价的中压微网混沌pso最优潮流实现方法 |
CN103440370A (zh) * | 2013-08-21 | 2013-12-11 | 国家电网公司 | 一种输变电工程造价评估方法及装置 |
Non-Patent Citations (1)
Title |
---|
杨杰等: "基于RBF径向神经网络的混沌时间序列预测", 《交通与计算机》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110555469A (zh) * | 2019-08-15 | 2019-12-10 | 阿里巴巴集团控股有限公司 | 处理交互序列数据的方法及装置 |
US10936950B1 (en) | 2019-08-15 | 2021-03-02 | Advanced New Technologies Co., Ltd. | Processing sequential interaction data |
US11636341B2 (en) | 2019-08-15 | 2023-04-25 | Advanced New Technologies Co., Ltd. | Processing sequential interaction data |
CN112784887A (zh) * | 2021-01-12 | 2021-05-11 | 北京林业大学 | 一种植物生命状态评价方法及装置 |
CN112784887B (zh) * | 2021-01-12 | 2024-01-05 | 北京林业大学 | 一种植物生命状态评价方法及装置 |
CN114844227A (zh) * | 2022-07-04 | 2022-08-02 | 广东电网有限责任公司佛山供电局 | 一种电网运行安全监督管理系统及方法 |
CN114844227B (zh) * | 2022-07-04 | 2022-09-20 | 广东电网有限责任公司佛山供电局 | 一种电网运行安全监督管理系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109993346B (zh) | 2020-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109993346A (zh) | 基于混沌时间序列和神经网络的微电网电压安全评估方法 | |
CN110221225B (zh) | 一种航天器锂离子电池循环寿命预测方法 | |
CN107123988A (zh) | 一种考虑恢复量不确定的停电电网负荷恢复二阶锥规划方法 | |
CN106786608B (zh) | 一种适用于分布式电源接入的不确定谐波潮流计算方法 | |
Zu et al. | Short-term wind power prediction method based on wavelet packet decomposition and improved GRU | |
CN114595891A (zh) | 一种配电网电压及潮流越界风险评估方法、系统及设备 | |
CN110363334A (zh) | 基于灰色神经网络模型的光伏并网的电网线损预测方法 | |
Naitmalek et al. | Forecasting the state-of-charge of batteries in micro-grid systems | |
Jannati et al. | Mitigation of windfarm power fluctuation by adaptive linear neuron‐based power tracking method with flexible learning rate | |
Li et al. | Series Dc arc fault detection and location in wind-solar-storage hybrid system based on variational mode decomposition | |
Li et al. | Forecasting of wind capacity ramp events using typical event clustering identification | |
CN113344283B (zh) | 基于边缘智能的能源互联网新能源消纳能力评估方法 | |
Dinavahi et al. | ANN based voltage stability margin prediction | |
Han et al. | Analysis and prediction of the penetration of renewable energy in power systems using artificial neural network | |
Pal et al. | Classification of fault using artificial neural network and power quality improvement using DVR in a PV integrated hybrid power system | |
Wang et al. | A wind speed prediction method based on improved empirical mode decomposition and support vector machine | |
Uzun et al. | Rule extraction and performance estimation by using variable neighborhoodsearch for solar power plant in Konya | |
CN117335426A (zh) | 一种适用于短期电压稳定性评估的关键特征识别方法 | |
Guo et al. | Online model-free estimation of the dynamic system model for a power system with renewables in ambient conditions | |
Prema et al. | Application of hybrid neuro-wavelet models for effective prediction of wind speed | |
Liu et al. | Intelligent voltage prediction of active distribution network with high proportion of distributed photovoltaics | |
Zhang et al. | CNN‐LSTM based power grid voltage stability emergency control coordination strategy | |
Liu et al. | Ultra-short-term forecasting method of wind power based on W-BiLSTM | |
Zhao et al. | Application and research of short-term Power load forecasting based on LSTM | |
CN109066683B (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 |