CN111478316A - 用于电力系统灵敏度在线辨识的噪声自助集成回归方法 - Google Patents

用于电力系统灵敏度在线辨识的噪声自助集成回归方法 Download PDF

Info

Publication number
CN111478316A
CN111478316A CN202010248456.5A CN202010248456A CN111478316A CN 111478316 A CN111478316 A CN 111478316A CN 202010248456 A CN202010248456 A CN 202010248456A CN 111478316 A CN111478316 A CN 111478316A
Authority
CN
China
Prior art keywords
power system
noise
sample
data
sensitivity
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
Application number
CN202010248456.5A
Other languages
English (en)
Other versions
CN111478316B (zh
Inventor
张俊勃
唐伟宁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN202010248456.5A priority Critical patent/CN111478316B/zh
Publication of CN111478316A publication Critical patent/CN111478316A/zh
Application granted granted Critical
Publication of CN111478316B publication Critical patent/CN111478316B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种用于电力系统灵敏度在线辨识的噪声自助集成回归方法,该方法在对灵敏度参数矩阵进行辨识的过程中,通过向数据中添加白噪声的方法,以降低数据中测量噪声的负面影响。噪声自助集成的回归方法计算复杂度和结果的平均值等价于含L2范数的局部加权最小二乘回归方法。相对于传统的局部加权最小二乘回归方法,该方法解决了电力系统采样的数据集中包含的测量噪声导致回归结果误差较大的问题,更适合于噪声较强且数据共线性较大的电力系统数据集。本发明提出的噪声自助集成回归方法有利于电力系统线性化模型的辨识。

Description

用于电力系统灵敏度在线辨识的噪声自助集成回归方法
技术领域
本发明涉及电力系统辨识和数据回归技术领域,具体涉及一种用于电力系统灵敏度参数矩阵在线辨识的噪声自助集成回归方法。
背景技术
参数灵敏度矩阵辨识广泛应用于电力系统的运行和控制,如电压灵敏度用于无功功率规划和电压稳定性控制、功率损耗灵敏度用于经济运行等。传统的电力系统灵敏度分析方法依赖离线系统模型。该模型通常由根据物理公式建立,其参数在实际测试中进行校验。随着大电网互连和新能源接入,现代电力系统的动态行为变得复杂且难以预测。离线物理模型很难与实时有功功率系统匹配,此时,基于模型的常规灵敏度分析并不完全适用,灵敏度在线辨识逐渐成为新的需求。
目前,随着广域测量系统和先进的计量设备的快速发展,电力系统操作员可以在短时间内收集大量的系统运行数据,这为灵敏度在线辨识提供了数据基础。同时,最新的数据处理算法为灵敏度在线辨识提供了解决方案,从而使其在实际电力系统中应用成为可能。但对于实际电力系统,需要解决数据共线性高、存在噪声等问题。
发明内容
本发明的目的是为了解决现有技术中的上述缺陷,提供一种用于电力系统灵敏度在线辨识的噪声自助集成回归方法。
本发明的目的可以通过采取如下技术方案达到:
一种用于电力系统灵敏度在线辨识的噪声自助集成回归方法,所述噪声自助集成回归方法包括数据预处理和循环过程,
首先,所述数据预处理包括以下步骤:
S1、进行数据采集:测量并采集机组或符合母线注入/吸收功率,作为电力系统自变量样本集X;测量并采集其中一个母线电压,作为电力系统因变量数据集Y,并将电力系统自变量样本集X和电力系统因变量数据集Y存储在同一个数据文件中;所述电力系统自变量样本集X为矩阵,每一行为一个样本,每一列为样本的一个变量;所述电力系统因变量数据集Y为列向量,每一行为一个样本;
S2、对电力系统自变量样本集X和电力系统因变量数据集Y分别按行求其均值,记为
Figure BDA0002434628780000021
再对电力系统自变量样本集X和电力系统因变量数据集分别按行求其方差,记为σX、σY;对电力系统自变量样本集X和电力系统因变量数据集标准化得到X1、Y1
S3、根据公式di=||x1,k-x1,i||2计算每个样本与给定的工作点之间的欧拉距离,其中x1,k表示X1中给定工作点k对应的样本,x1,i表示X1中的第i个样本。依据欧拉距离对样本按升序排序,并选出前n个最近的点作为训练集;
S4、根据上一步得到的欧拉距离,利用公式
Figure BDA0002434628780000022
计算分配给训练集样本的权重wi,其中
Figure BDA0002434628780000023
表示对e取
Figure BDA0002434628780000024
次方,τ是权重公式的可调参数;
然后,所述循环包括以下步骤:
S5、根据公式
Figure BDA0002434628780000025
向X1、Y1添加白噪声,其中ξ1、ζ1分别表示自变量和因变量中的噪声分量,
Figure BDA0002434628780000026
是第l次人为加入的白噪声分量,其长度分别与
Figure BDA0002434628780000027
Y1 (l)相同,
Figure BDA0002434628780000028
Y1 (l)分别是人为生成的第l次含白噪声数据集;
S6、根据局部加权最小二乘回归准则函数
Figure BDA0002434628780000031
对数据集X1、Y1进行回归,其中
Figure BDA0002434628780000032
是待求解的灵敏度参数矩阵,W=diag{wi}是根据步骤S4解得的权重矩阵;
S7、对步骤S6中准则函数J1进行矩阵变换求解,得到局部加权最小二乘回归解
Figure BDA0002434628780000033
求解得到灵敏度参数矩阵
Figure BDA0002434628780000034
S8、重复步骤S5-S7,且每次循环过程中使用不同的信噪比大小的白噪声序列;当循环次数Nl=N时,认为集合样本足够,此时取所有
Figure BDA0002434628780000035
平均值作为最终结果。
进一步地,所述n取值区间为【800,2000】。
进一步地,所述可调参数τ取值区间为【0.1,0.5】。
进一步地,所述N取值区间为【80,200】。
进一步地,所述白噪声的生成方法为使用MATLAB软件的AWGN函数对原数据序列直接生成信噪比为30dB的随机白噪声序列。
本发明相对于现有技术具有如下的优点及效果:
本发明提出的一种用于电力系统灵敏度参数矩阵在线辨识的噪声自助集成回归方法,在对灵敏度参数矩阵进行辨识的过程中,通过向数据中添加白噪声的方法,以降低数据中测量噪声的负面影响。噪声自助集成的回归方法计算复杂度和结果的平均值等价于含L2范数的局部加权最小二乘回归方法。相对于传统的局部加权最小二乘回归方法,该方法解决了电力系统采样的数据集中包含的测量噪声导致回归结果误差较大的问题,更适合于噪声较强且数据共线性较大的电力系统数据集。本发明提出的噪声自助集成回归方法有利于电力系统线性化模型的辨识。
附图说明
图1是本发明的一个实施例所依据的电力系统接线图;
图2是本发明的一个实施例中噪声自助集成回归方法与传统方法在含共线性数据的辨识效果对比图;
图3是本发明的一个实施例中噪声自助集成回归方法与传统方法在含测量噪声及共线性数据的辨识效果对比图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例
给定电力系统在平衡工作点附近稳定运行。则灵敏度方程可表示为:
Figure BDA0002434628780000041
其中p表示电力系统参数变量的向量,Δo表示某一可测得或不可测得的因变量,f表示非线性系统模型,
Figure BDA0002434628780000042
表示在平衡工作点附近稳定运行时Δo对p的灵敏度向量,ε表示误差项,假定ε的分布服从一个均值为0的正态分布。
灵敏度参数矩阵的目标即是求得在给定电力系统运行点k附近的线性灵敏度向量
Figure BDA0002434628780000043
因此该过程可以用一个基本多变量线性回归模型表示:y=β0+xβ+ε,其中x是长度为1×p的向量,其中包含在运行点k处的p个电力系统参数变量Δp;y代表运行点k处的因变量Δo;β是长度为p×1的待求灵敏度向量;β0为常数项。假设该模型在同一工作点k含有样本数为n,则该模型可以改写为Y=β0+Xβ+ε,其中X是包含n个样本和p个变量的n×p自变量矩阵,Y是长度为n×1包含n个样本的因变量向量,ε是长度为n×1包含n个误差项的误差向量。
噪声自助集成回归方法是指,在对灵敏度参数矩阵进行辨识的过程中,通过向数据中添加白噪声的方法,以降低数据中测量噪声的负面影响。对于得到的含白噪声的数据集,通过局部加权最小二乘估计对其回归。
本实施例提出了一种用于电力系统灵敏度参数矩阵在线辨识的噪声自助集成回归方法,包括以下步骤:
首先,进行数据预处理:
S1、进行数据采集:测量并采集机组或符合母线注入/吸收功率,作为电力系统自变量样本集X;测量并采集其中一个母线电压,作为电力系统因变量数据集Y,并将电力系统自变量样本集X和电力系统因变量数据集Y存储在同一个数据文件中;所述电力系统自变量样本集X为矩阵,每一行为一个样本,每一列为样本的一个变量;所述电力系统因变量数据集Y为列向量,每一行为一个样本。
S2、对电力系统自变量样本集X和电力系统因变量数据集Y分别按行求其均值,记为
Figure BDA0002434628780000051
再对电力系统自变量样本集X和电力系统因变量数据集分别按行求其方差,记为σX、σY;对电力系统自变量样本集X和电力系统因变量数据集标准化得到X1、Y1
S3、根据公式di=||x1,k-x1,i||2计算每个样本与给定的工作点之间的欧拉距离,其中x1,k表示X1中给定工作点k对应的样本,x1,i表示X1中的第i个样本。依据欧拉距离对样本按升序排序,并选出前n个最近的点作为训练集,其中n取典型值1000;
S4、根据上一步得到的欧拉距离,利用公式
Figure BDA0002434628780000052
计算分配给训练集样本的权重wi,其中
Figure BDA0002434628780000053
表示对e取
Figure BDA0002434628780000054
次方,τ是权重公式的可调参数,取典型值τ=0.2,di是步骤S3中得到的欧拉距离;
然后,循环过程;
S5、根据公式
Figure BDA0002434628780000061
向X1、Y1添加白噪声,其中ξ1、ζ1分别表示自变量和因变量中的噪声分量,
Figure BDA0002434628780000062
是第l次人为加入的白噪声分量,其长度分别与
Figure BDA0002434628780000063
Y1 (l)相同,
Figure BDA0002434628780000064
Y1 (l)分别是人为生成的第l次含白噪声数据集。生成白噪声的方法为使用MATLAB软件的AWGN函数对原数据序列直接生成信噪比为30dB的随机白噪声序列。
S6、根据局部加权最小二乘回归准则函数
Figure BDA0002434628780000065
对数据集X1、Y1进行回归,其中
Figure BDA0002434628780000066
是待求解的灵敏度参数矩阵,W=diag{wi}是根据步骤S4解得的权重矩阵;
S7、对步骤S6中准则函数J1进行矩阵变换求解,可得到局部加权最小二乘回归解
Figure BDA0002434628780000067
求解得到灵敏度参数矩阵
Figure BDA0002434628780000068
S8、重复步骤S5-S7,且每次循环过程中使用不同的信噪比大小的白噪声序列;当循环次数Nl=N时,认为集合样本足够,此时取所有
Figure BDA0002434628780000069
平均值作为最终结果。其中N取典型值100。
用本实施例公开的噪声自助集成回归方法应用在实际8机36节点电力系统中应用,并进行在线电压灵敏度辨识。通过MATLAB中的PSAT软件进行仿真,每隔5分钟进行采样,得到共计8760组数据样本作为样本集1。将样本集1数据采用采用上文叙述的方法加入信噪比为20dB的白噪声,以模拟实际电力系统中可能存在的强共线性数据,作为样本集2。对样本集1和样本集2加入信噪比为40dB的白噪声,以模拟实际系统中存在的测量噪声,得到样本集3、样本集4。采用上述样本集验证本实施例公开的噪声自助集成回归方法对于强共线性和测量噪声大数据集的适用性。
图1是本实施例依据的IEEE8机36节点电力系统[1],其包含8台发电机、36条母线、13处负荷和42条输电线路。其中发电机1用于功率平衡,发电机6是同步调相机。在电压灵敏度辨识中,选取母线14的电压作为输出变量,系统各节点有功潮流作为输入变量。采用PSAT软件进行仿真,其中选定权重参数τ=0.2,回归样本集的组数n取典型值1000,人为生成的噪声序列信噪比选定为30dB,循环次数N取100。
图2是对于样本集2和4分别采用噪声自助集成回归方法(图中简称NAER)和局部加权最小二乘法(图中简称LWLR)在线电压灵敏度辨识的对比结果。可以看出,由于样本集2中存在的共线性问题,局部加权最小二乘法得到了误差较大的结果,作为对比,噪声自助集成回归的方法得到的电压灵敏度与真实值差距明显较小。对于样本集4的辨识结果可以看出,噪声自助集成回归方法在存在测量噪声的条件下辨识效果较局部加权最小二乘法更好。
图3是对于样本集3和4分别采用噪声自助集成回归方法(图中简称NAER)和局部加权最小二乘法(图中简称LWLR)在线电压灵敏度辨识的对比结果。可以看出,噪声自助集成回归方法和局部加权最小二乘法均取得较好结果。同时,对于样本集3和4的辨识,噪声自助集成回归方法取得了相同的结果,说明该方法将白噪声添加到回归数据中不会影响辨识的准确性。
综上所述,可以认为本实施例提出的噪声自助集成回归方法在电力系统灵敏度辨识中具有较好的抗噪声和反共线性作用。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。

Claims (5)

1.一种用于电力系统灵敏度在线辨识的噪声自助集成回归方法,其特征在于,所述噪声自助集成回归方法包括数据预处理和循环过程,
首先,所述数据预处理包括以下步骤:
S1、进行数据采集:测量并采集机组或符合母线注入/吸收功率,作为电力系统自变量样本集X;测量并采集其中一个母线电压,作为电力系统因变量数据集Y,并将电力系统自变量样本集X和电力系统因变量数据集Y存储在同一个数据文件中;所述电力系统自变量样本集X为矩阵,每一行为一个样本,每一列为样本的一个变量;所述电力系统因变量数据集Y为列向量,每一行为一个样本;
S2、对电力系统自变量样本集X和电力系统因变量数据集Y分别按行求其均值,记为
Figure FDA0002434628770000011
再对电力系统自变量样本集X和电力系统因变量数据集分别按行求其方差,记为σX、σY;对电力系统自变量样本集X和电力系统因变量数据集标准化得到X1、Y1
S3、根据公式di=||x1,k-x1,i||2计算每个样本与给定的工作点之间的欧拉距离,其中x1,k表示X1中给定工作点k对应的样本,x1,i表示X1中的第i个样本,依据欧拉距离对样本按升序排序,并选出前n个最近的点作为训练集;
S4、根据上一步得到的欧拉距离,利用公式
Figure FDA0002434628770000012
计算分配给训练集样本的权重wi,其中
Figure FDA0002434628770000013
表示对e取
Figure FDA0002434628770000014
次方,τ是权重公式的可调参数;
然后,所述循环包括以下步骤:
S5、根据公式
Figure FDA0002434628770000015
向X1、Y1添加白噪声,其中ξ1、ζ1分别表示自变量和因变量中的噪声分量,
Figure FDA0002434628770000016
是第l次人为加入的白噪声分量,其长度分别与
Figure FDA0002434628770000021
Y1 (l)相同,
Figure FDA0002434628770000022
Y1 (l)分别是人为生成的第l次含白噪声数据集;
S6、根据局部加权最小二乘回归准则函数
Figure FDA0002434628770000023
对数据集X1、Y1进行回归,其中
Figure FDA0002434628770000024
是待求解的灵敏度参数矩阵,W=diag{wi}是根据步骤S4解得的权重矩阵;
S7、对步骤S6中准则函数J1进行矩阵变换求解,得到局部加权最小二乘回归解
Figure FDA0002434628770000025
求解得到灵敏度参数矩阵
Figure FDA0002434628770000026
S8、重复步骤S5-S7,且每次循环过程中使用不同的信噪比大小的白噪声序列;当循环次数Nl=N时,认为集合样本足够,此时取所有
Figure FDA0002434628770000027
平均值作为最终结果。
2.根据权利要求1所述的用于电力系统灵敏度在线辨识的噪声自助集成回归方法,其特征在于,所述n取值区间为【800,2000】。
3.根据权利要求1所述的用于电力系统灵敏度在线辨识的噪声自助集成回归方法,其特征在于,所述可调参数τ取值区间为【0.1,0.5】。
4.根据权利要求1所述的用于电力系统灵敏度在线辨识的噪声自助集成回归方法,其特征在于,所述N取值区间为【80,200】。
5.根据权利要求1所述的用于电力系统灵敏度在线辨识的噪声自助集成回归方法,其特征在于,所述白噪声的生成方法为使用MATLAB软件的AWGN函数对原数据序列直接生成信噪比为30dB的随机白噪声序列。
CN202010248456.5A 2020-04-01 2020-04-01 用于电力系统灵敏度在线辨识的噪声自助集成回归方法 Active CN111478316B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010248456.5A CN111478316B (zh) 2020-04-01 2020-04-01 用于电力系统灵敏度在线辨识的噪声自助集成回归方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010248456.5A CN111478316B (zh) 2020-04-01 2020-04-01 用于电力系统灵敏度在线辨识的噪声自助集成回归方法

Publications (2)

Publication Number Publication Date
CN111478316A true CN111478316A (zh) 2020-07-31
CN111478316B CN111478316B (zh) 2022-06-10

Family

ID=71749457

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010248456.5A Active CN111478316B (zh) 2020-04-01 2020-04-01 用于电力系统灵敏度在线辨识的噪声自助集成回归方法

Country Status (1)

Country Link
CN (1) CN111478316B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103258103A (zh) * 2013-05-29 2013-08-21 国家电网公司 基于偏最小二乘回归的戴维南等值参数辨识方法
US20170061305A1 (en) * 2015-08-28 2017-03-02 Jiangnan University Fuzzy curve analysis based soft sensor modeling method using time difference Gaussian process regression
CN109344870A (zh) * 2018-08-28 2019-02-15 华南理工大学 一种基于目标正则化的电力系统状态灵敏度在线辨识方法
CN109885807A (zh) * 2019-01-22 2019-06-14 哈尔滨工业大学(深圳) Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103258103A (zh) * 2013-05-29 2013-08-21 国家电网公司 基于偏最小二乘回归的戴维南等值参数辨识方法
US20170061305A1 (en) * 2015-08-28 2017-03-02 Jiangnan University Fuzzy curve analysis based soft sensor modeling method using time difference Gaussian process regression
CN109344870A (zh) * 2018-08-28 2019-02-15 华南理工大学 一种基于目标正则化的电力系统状态灵敏度在线辨识方法
CN109885807A (zh) * 2019-01-22 2019-06-14 哈尔滨工业大学(深圳) Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
何楚瑶: "无功-电压控制特性动态辨识技术及应用研究", 《中国优秀硕士论文全文数据库》 *
胡清华 等: "考虑数据噪声的鲁棒回归建模方法综述", 《西北大学学报(自然科学版)》 *

Also Published As

Publication number Publication date
CN111478316B (zh) 2022-06-10

Similar Documents

Publication Publication Date Title
Liu et al. Comparisons on Kalman-filter-based dynamic state estimation algorithms of power systems
CN107590317B (zh) 一种计及模型参数不确定性的发电机动态估计方法
CN112636341B (zh) 基于多新息辨识的电力系统惯量空间分布估计方法及装置
CN106353623B (zh) 基于随机响应信号的电力系统低频振荡模式在线辨识方法
CN114006370B (zh) 一种电力系统暂态稳定分析评估方法及系统
CN106549396A (zh) 一种配电网多目标概率无功优化方法
Hou et al. Cluster computing-based trajectory sensitivity analysis application to the WECC system
CN114595858A (zh) 基于滚动时间序列和支持向量机的短期风速预测方法和系统
CN109066651B (zh) 风电-负荷场景的极限传输功率的计算方法
CN114583767B (zh) 一种数据驱动的风电场调频响应特性建模方法及系统
CN117471252A (zh) 基于iceemdan-kpca特征提取与soa-kelm的直流故障电弧诊断方法
CN117349614A (zh) 基于自注意力机制和时空图卷积网络的频率稳定预测方法
CN114372344B (zh) 一种次同步振荡阻尼特性影响因素的量化辨识方法
CN114611676A (zh) 基于神经网络的新能源发电系统阻抗模型辨识方法及系统
Wang et al. Continual residual reservoir computing for remaining useful life prediction
CN111478316B (zh) 用于电力系统灵敏度在线辨识的噪声自助集成回归方法
CN117748481A (zh) 基于实时动态分区的电力系统惯量在线评估方法及装置
CN105808825B (zh) 一种适于有源配电网动态仿真的滚动式随机投影积分方法
CN111639463A (zh) 一种基于XGBoost算法的电力系统扰动后频率特征预测方法
CN115764976A (zh) 考虑不确定性系统辨识的电力系统惯量连续评估方法
CN107947179A (zh) 一种考虑风电调度策略的多点线性化随机潮流计算方法
Xinhui et al. A data-driven slow dynamic characteristic extraction and state estimation method for large power grid
CN113933706A (zh) 一种基于bp神经网络的直流电机故障检测方法
CN112861074A (zh) 基于Hankel-DMD的电力系统机电参数提取方法
Li et al. Equivalent Virtual Inertia Estimation of Wind Turbines Based on GS-LSTM

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