CN107093904B - 基于随机矩阵谱分析的配电网无功补偿动作点选取方法 - Google Patents
基于随机矩阵谱分析的配电网无功补偿动作点选取方法 Download PDFInfo
- Publication number
- CN107093904B CN107093904B CN201610091486.3A CN201610091486A CN107093904B CN 107093904 B CN107093904 B CN 107093904B CN 201610091486 A CN201610091486 A CN 201610091486A CN 107093904 B CN107093904 B CN 107093904B
- Authority
- CN
- China
- Prior art keywords
- matrix
- reactive compensation
- distribution network
- power distribution
- reactive
- 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
- 239000011159 matrix material Substances 0.000 title claims abstract description 157
- 238000009826 distribution Methods 0.000 title claims abstract description 45
- 238000010183 spectrum analysis Methods 0.000 title claims abstract description 16
- 238000010187 selection method Methods 0.000 title claims abstract description 7
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 238000002347 injection Methods 0.000 claims abstract description 8
- 239000007924 injection Substances 0.000 claims abstract description 8
- 239000013598 vector Substances 0.000 claims description 37
- 230000003595 spectral effect Effects 0.000 claims description 9
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 238000001228 spectrum Methods 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 238000010206 sensitivity analysis Methods 0.000 abstract description 3
- 238000005457 optimization Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- 230000001939 inductive effect Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000010223 real-time analysis Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
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/18—Arrangements for adjusting, eliminating or compensating reactive power in networks
- H02J3/1821—Arrangements for adjusting, eliminating or compensating reactive power in networks using shunt compensators
- H02J3/1871—Methods for planning installation of shunt reactive power compensators
-
- 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/30—Reactive power compensation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种基于随机矩阵谱分析的配电网无功补偿动作点选取方法,适用于容性无功补偿,所述方法包括:步骤(1):获取配电网中电压和节点注入无功功率的历史数据,并通过所述历史数据构建随机矩阵模型;步骤(2):通过对随机矩阵模型的特征根分析,确定配电网无功补偿动作点的个数;步骤(3):定义节点无功补偿度,并通过节点无功补偿度确定无功补偿动作点。本发明技术方案计算速度显著提高;且更具通用性;同时避免了应用灵敏度分析法时,无功补偿动作节点个数不确定的问题。
Description
技术领域
本发明涉及配电网无功优化领域,更具体涉及一种基于随机矩阵谱分析的配电网无功补偿动作点选取方法。
背景技术
随着大数据时代的到来,大数据全生命周期中获取、存储、处理、挖掘等方面的研究逐步展开,成为行业、政府的研究热点。在电力行业,大数据来源于智能电表、配电自动化、数字保护装置、其他智能设备以及智能电网实施过程中产生的相关数据。如何有效地利用这些数据,就成为电力公司不得不认真面对的一个挑战。
目前,大数据技术在电力系统中的研究成果寥寥无几,本发明首个提出利用配电网大数据解决无功优化控制中动作节点的选择问题。大数据技术可处理配电网中结构多样、来源复杂、时空尺度各异的数据,分析出与配电网规划、运行和控制相关的、未知的、潜在的、有价值的信息,为实现智能配电网实时分析具有重要意义。
发明内容
本发明的目的是提供一种基于随机矩阵谱分析的配电网无功补偿动作点选取方法,计算速度显著提高;且更具通用性;同时避免了应用灵敏度分析法时,无功补偿动作节点个数不确定的问题。
为实现上述目的,本发明采用以下技术方案:一种基于随机矩阵谱分析的配电网无功补偿动作点选取方法,适用于容性无功补偿,所述方法包括:
步骤(1):获取配电网中电压和节点注入无功功率的历史数据,并通过所述历史数据构建随机矩阵模型;
步骤(2):通过对随机矩阵模型的特征根分析,确定配电网无功补偿动作点的个数;
步骤(3):定义节点无功补偿度,并通过节点无功补偿度确定无功补偿动作点。
所述步骤(1)的构建随机矩阵模型的过程包括:
通过历史数据构造节点电压时间序列矩阵V和节点注入无功时间序列矩阵Q;
通过V矩阵计算电压偏差矩阵ΔV;
根据向量元素确定随机矩阵模型。
所述电压偏差矩阵ΔV中第i行j列向量元素为:Δvi,j=1-vi,j;其中,vi,j为V矩阵中第i行j列向量元素。
其中,为标准化后的矩阵第i行j列向量,Δvi,j为电压偏差矩阵ΔV中第i行j列向量元素,为电压偏差矩阵中第j列向量均值,为标准化后的矩阵的第i行向量,Δvi为电压偏差矩阵ΔV中第i行向量,为标准化后的矩阵的第i行j列向量,qi,j为Q矩阵中第i行j列向量元素,为Q矩阵中第j列向量均值,为标准化后的矩阵Q的第i行向量,qi为Q矩阵中第i行向量,为1;标准化后的矩阵和的各列向量均值和为0;V和Q矩阵均为T行N列矩阵,其中,T为样本个数,N为节点个数,且为标幺值矩阵。
所述随机矩阵的架构X中每一项为:
所述步骤(2)通过对随机矩阵模型的特征根分析,确定配电网无功补偿动作点的个数的过程包括;
对随机矩阵进行处理,得到实对称矩阵S;
对S矩阵进行谱分解,确定随机矩阵的特征根;
确定随机矩阵的较大特征根,并根据较大特征根确定电网无功补偿动作点的个数。
所述实对称矩阵S为:
S=X'X (3)
按下式(4)对S矩阵进行谱分解:
其中,X'为X矩阵的转置矩阵,Λ=矩阵S对角线元素diag(λ1',...,λn');U=(u1,...,un)为n维空间上的一组标准正交基,其对应的特征根记为由大到小排列的λ1',...,λn';u'k为uk的转置;U'为U矩阵的装置矩阵;
满足式(5)条件的特征根λk,称为较大特征根:
其中,e代表特征根的收敛速度,根据经验调整e的取值大小;所求得的较大特征根的个数即为需要投切无功补偿装置的节点个数。
所述步骤(3)中定义节点无功补偿度:将随机矩阵X中j节点所对应的j列元素置为0,记变化后的矩阵为Mj,按下式计算矩阵Mj的谱半径,其值大小称为第j节点无功补偿度:
mj=ρ(Mj)=max{λi|,i=1,2,...n} (6)
其中,λi是其特征根;选出mj中最小的需要投切无功补偿装置的节点个数的值,所对应的节点j即为补偿动作点的位置。
所述配电网为同一区域、同一电压等级和馈线相互连接的网络,且网络结构为辐射状;
所述历史数据为计算时刻前24小时的数据,且认为该时间段网络拓扑不变。
和最接近的现有技术比,本发明提供技术方案具有以下优异效果:
利用本发明涉及的一种用于配电网无功优化中无功控制的关键点选择的大数据处理方法,将随机矩阵谱分析这种数学方法应用到一个新的领域,即配电网无功优化领域,确定配电网无功补偿动作点的个数和位置。由配电网节点历史数据构造大数据随机矩阵架构,通过对该矩阵特征根的分析,得出无功补偿动作点选择的结果。本方法无需考虑网络结构、参数,也无需额外建立分布式电源、负荷等模型,可直接从历史数据的统计信息快速地得到无功补偿动作点的结果;与现有的通过无功优化的方法选择无功补偿动作点相比,计算速度显著提高;与通过历史经验选择无功补偿动作点方法相比,该算法更具通用性;同时避免了应用灵敏度分析法时,无功补偿动作节点个数不确定的问题。且通过各节点的历史数据得到最优无功补偿动作节点的位置,不必知道网络参数,大大提高了实用性;同时不受负荷小波动干扰的影响,结果具有稳定性。
附图说明
图1为本发明实施例的方法流程图;
图2为本发明实施例的IEEE33节点算例图。
具体实施方式
对配电网做如下假设:
(1)配电网为同一区域、同一电压等级、馈线相互连接的网络,且网络结构为辐射状,暂不考虑环型配电网的情况;
(2)本专利所用历史数据一般为计算时刻前24小时的数据,且认为该时间段网络拓扑不变;
(3)一般配电网感性负荷较大,本专利只考虑容性无功补偿,暂不考虑感性无功补偿;
(4)配电网自动化水平较高,可获取全部节点的历史数据。
所用变量术语如表1所示。
表1所用变量术语说明
基于随机矩阵谱分析的配电网无功补偿动作点选择方法,通过对历史数据处理构成的随机矩阵特征根进行分析,得到无功补偿动作点的结果。该方法包括如图1所示的如下步骤:
步骤(1):获取配电网中电压和节点注入无功功率的历史数据,由此构建随机矩阵模型;
步骤(2):通过对随机矩阵模型的特征根分析,确定配电网无功补偿动作点的个数;
步骤(3):定义节点无功补偿度,并由此确定无功补偿动作点。
上述步骤具体包括如下内容:
步骤(1)具体包括:应用某一配电网某一天的数据(或根据前一周数据计算各时刻平均值),分别构造节点电压时间序列矩阵V和节点注入无功时间序列矩阵Q,V、Q矩阵均为T行N列矩阵,其中T为样本个数,N为节点个数,且为标幺值矩阵,并由V矩阵计算电压偏差矩阵ΔV,矩阵中第i行j列元素Δvi,j=1-vi,j。
因此,构建随机矩阵架构X,X中每一项为
物理意义:当节点无功负荷增大时,线路损耗将增大、电压偏差增大;由无功补偿作用可知,应使这些节点处的无功补偿装置动作,进行无功补偿。
步骤(2)具体包括:由于所构造的随机矩阵为非对称矩阵,不能得到理想的谱分解,所以,先对随机矩阵进行如下处理:
S=X'X (3)
得到矩阵S为实对称矩阵,对S矩阵进行谱分解,如公式(4):
其中,Λ=diag(λ1',...,λn'),U=(u1,...,un)为n维空间上的一组标准正交基,它们对应的特征根记为λ1',...,λn'(由大到小排列)。
将矩阵进行谱分解的意义就是,把矩阵所代表的空间进行正交分解,使得矩阵的向量集合可以表示为每个向量在各个特征向量上面的投影长度,矩阵特征根表示了矩阵在其对应的特征向量上所含信息量的大小。
定义较大特征根:满足式(5)条件的特征根λk,称为较大特征根。
其中,e代表特征根的收敛速度,可根据经验适当调整e的取值大小。算例中对于IEEE33节点系统,e取0.2。当e增大时,补偿动作点的个数将有可能减少;反之,当e减小时,补偿动作点的个数将可能增加。所求得的较大特征根的个数即为需要投切无功补偿装置的节点个数,将该节点个数用“a”表示。
步骤(3)具体包括:
矩阵谱半径的定义:设A是N×N矩阵,λi是其特征根,称ρ(A)=max{λi|,i=1,2,...n}为A的谱半径。
定义节点无功补偿度:将X矩阵中j节点所对应的j列元素置为0,记变化后的矩阵为Mj,计算Mj的谱半径,其值大小称为第j节点无功补偿度。
mj=ρ(Mj)=max{λi|,i=1,2,...n} (6)
选出mj中最小的a个(a为需要投切无功补偿装置的节点个数),所对应的节点j即为补偿动作点的位置。
下面结合实施例对发明作进一步的详细说明。
实施例1:
以IEEE33节点系统为算例,网络接线如图2所示。应用配电网典型日的历史数据,分别构造节点电压时间序列矩阵V和节点注入无功时间序列矩阵Q,V、Q矩阵均为T行N列矩阵,其中T为样本个数,N为节点个数,且为标幺值矩阵,并由V矩阵计算电压偏差矩阵ΔV,矩阵中第i行j列元素Δvi,j=1-vi,j。
因此,构建随机矩阵架构X,X中每一项为
由于所构造的随机矩阵为非对称矩阵,不能得到理想的谱分解,所以,先对随机矩阵进行如下处理:
S=X'X
得到矩阵S为实对阵矩阵,对S矩阵进行谱分解,如公式(10):
其中,n=N,Λ=diag(λ1',...,λn'),U=(u1,...,un)为n维空间上的一组标准正交基,它们对应的特征根记为λ1',...,λn'(由大到小排列)。
定义节点无功补偿度:将X矩阵中j节点所对应的j列元素置为0,记变化后的矩阵为Mj,计算Mj的谱半径,其值大小称为第j节点无功补偿度。
mj=ρ(Mj)=max{λi|,i=1,2,...n}
其中,Mj表示将矩阵X的第j列置为0后的矩阵。
选出mj中较小的a个,所对应的节点j即为补偿动作点位置。由算例得到的无功动作点的位置为29节点和7节点。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,所属领域的普通技术人员尽管参照上述实施例应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,这些未脱离本发明精神和范围的任何修改或者等同替换,均在申请待批的本发明的权利要求保护范围之内。
Claims (8)
1.一种基于随机矩阵谱分析的配电网无功补偿动作点选取方法,适用于容性无功补偿,其特征在于:所述方法包括:
步骤(1):获取配电网中电压和节点注入无功功率的历史数据,并通过所述历史数据构建随机矩阵模型;
步骤(2):通过对随机矩阵模型的特征根分析,确定配电网无功补偿动作点的个数;
步骤(3):定义节点无功补偿度,并通过节点无功补偿度确定无功补偿动作点;
所述步骤(1)的构建随机矩阵模型的过程包括:
通过历史数据构造节点电压时间序列矩阵V和节点注入无功时间序列矩阵Q;
通过V矩阵计算电压偏差矩阵△V;
根据向量元素确定随机矩阵模型;
所述步骤(3)中定义节点无功补偿度:将随机矩阵X中j节点所对应的j列元素置为0,记变化后的矩阵为Mj,按下式计算矩阵Mj的谱半径,其值大小称为第j节点无功补偿度:
mj=ρ(Mj)=max{|λi|,i=1,2,...n} (6)
其中,λi是其特征根;选出mj中最小的需要投切无功补偿装置的节点个数的值,所对应的节点j即为补偿动作点的位置。
2.如权利要求1所述的一种基于随机矩阵谱分析的配电网无功补偿动作点选取方法,其特征在于:所述电压偏差矩阵△V中第i行j列向量元素为:△vi,j=1-vi,j;其中,vi,j为V矩阵中第i行j列向量元素。
3.如权利要求1或2所述的一种基于随机矩阵谱分析的配电网无功补偿动作点选取方法,其特征在于:所述矩阵Q和△V按公式(1)做标准化处理,且标准化后的矩阵和的各列向量是和服从均值为0、方差σ为1的独立同分布;
5.如权利要求1所述的一种基于随机矩阵谱分析的配电网无功补偿动作点选取方法,其特征在于:所述步骤(2)通过对随机矩阵模型的特征根分析,确定配电网无功补偿动作点的个数的过程包括;
对随机矩阵进行处理,得到实对称矩阵S;
对S矩阵进行谱分解,确定随机矩阵的特征根;
确定随机矩阵的较大特征根,并根据较大特征根确定电网无功补偿动作点的个数。
8.如权利要求1所述的一种基于随机矩阵谱分析的配电网无功补偿动作点选取方法,其特征在于:所述配电网为同一区域、同一电压等级和馈线相互连接的网络,且网络结构为辐射状;
所述历史数据为计算时刻前24小时的数据,且认为该时间段网络拓扑不变。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610091486.3A CN107093904B (zh) | 2016-02-18 | 2016-02-18 | 基于随机矩阵谱分析的配电网无功补偿动作点选取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610091486.3A CN107093904B (zh) | 2016-02-18 | 2016-02-18 | 基于随机矩阵谱分析的配电网无功补偿动作点选取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107093904A CN107093904A (zh) | 2017-08-25 |
CN107093904B true CN107093904B (zh) | 2020-04-14 |
Family
ID=59648851
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610091486.3A Active CN107093904B (zh) | 2016-02-18 | 2016-02-18 | 基于随机矩阵谱分析的配电网无功补偿动作点选取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107093904B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107565550B (zh) * | 2017-09-11 | 2019-11-19 | 中国农业大学 | 一种配电网分区方法和系统 |
CN110954810B (zh) * | 2019-11-07 | 2022-02-18 | 广东电网有限责任公司 | 一种10kV高压断路器健康水平诊断方法 |
CN113507128B (zh) * | 2021-07-14 | 2023-07-21 | 国网青海省电力公司 | 一种特高压直流换流站近场无功优化配置方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104361200A (zh) * | 2014-10-10 | 2015-02-18 | 国家电网公司 | 基于灵敏度法的无功补偿优化选址方法 |
CN105244889A (zh) * | 2015-07-21 | 2016-01-13 | 河南行知专利服务有限公司 | 一种配电网无功补偿规划的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
PL199780B1 (pl) * | 2001-12-20 | 2008-10-31 | B Spo & Lstrok Ka Z Ograniczon | Sposób badania stabilności napięciowej systemu elektroenergetycznego |
-
2016
- 2016-02-18 CN CN201610091486.3A patent/CN107093904B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104361200A (zh) * | 2014-10-10 | 2015-02-18 | 国家电网公司 | 基于灵敏度法的无功补偿优化选址方法 |
CN105244889A (zh) * | 2015-07-21 | 2016-01-13 | 河南行知专利服务有限公司 | 一种配电网无功补偿规划的方法 |
Non-Patent Citations (3)
Title |
---|
基于GMR技术确定电压弱节点的特征根灵敏度指标;宗秀红 等;《继电器》;20060916;第34卷(第18期);35-39 * |
电力系统分区法确定交直流系统的最佳无功补偿点;王林川 等;《黑龙江电力》;20121231;第34卷(第6期);409-411 * |
高维随机矩阵描述下的量测大数据建模与异常数据检测方法;魏大千 等;《中国电机工程学报》;20150930;第35卷;59-66 * |
Also Published As
Publication number | Publication date |
---|---|
CN107093904A (zh) | 2017-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Münz et al. | Sensor and Actuator Placement for Linear Systems Based on $ H_ {2} $ and $ H_ {\infty} $ Optimization | |
CN109214708B (zh) | 基于交叉熵理论优化支持向量机的电力系统风险评估方法 | |
CN109597967B (zh) | 一种电力系统配变电站负荷数据异常检测与修复方法 | |
CN105846448B (zh) | 一种基于随机矩阵理论的配电网无功补偿容量确定方法 | |
CN107093904B (zh) | 基于随机矩阵谱分析的配电网无功补偿动作点选取方法 | |
CN107403196B (zh) | 基于谱聚类分析的即时学习建模对脱丁烷塔底丁烷浓度进行预测的方法 | |
Ma et al. | Fault detection based on statistical multivariate analysis and microarray visualization | |
CN106067034B (zh) | 一种基于高维矩阵特征根的配电网负荷曲线聚类方法 | |
CN112734128A (zh) | 一种基于优化rbf的7日电力负荷峰值预测方法 | |
CN116031888B (zh) | 基于动态负荷预测的潮流优化方法、系统及存储介质 | |
CN107133684A (zh) | 一种面向配电网无功优化的随机矩阵构建方法 | |
CN108647838B (zh) | 一种基于随机矩阵理论和在线序列极限学习机的电网态势感知方法 | |
CN111650472A (zh) | 一种电压暂降源定位的方法 | |
Wang et al. | A data-centric predictive control approach for nonlinear chemical processes | |
Fusco et al. | Power systems data fusion based on belief propagation | |
Pampuri et al. | Multistep virtual metrology approaches for semiconductor manufacturing processes | |
CN116992274B (zh) | 基于改进主成分回归模型的短期风速预测方法及系统 | |
CN111600309A (zh) | 电压控制方法、装置、设备、计算机设备和存储介质 | |
Xiong et al. | Batch-to-batch optimal control of nonlinear batch processes based on incrementally updated models | |
Tóth et al. | Order and structural dependence selection of LPV-ARX models revisited | |
CN116400266A (zh) | 基于数字孪生模型的变压器故障检测方法、装置及介质 | |
CN115619563A (zh) | 一种基于神经网络的股票价格分析方法 | |
Idema et al. | Fast Newton load flow | |
CN114925537A (zh) | 一种基于指定时间一致性的无初始化智能电网经济调度方法 | |
Deng et al. | A simulation-based method for aggregating Markov chains |
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 |