CN104617574B - 一种电力系统负荷区域暂态电压稳定的评估方法 - Google Patents
一种电力系统负荷区域暂态电压稳定的评估方法 Download PDFInfo
- Publication number
- CN104617574B CN104617574B CN201510025577.2A CN201510025577A CN104617574B CN 104617574 B CN104617574 B CN 104617574B CN 201510025577 A CN201510025577 A CN 201510025577A CN 104617574 B CN104617574 B CN 104617574B
- Authority
- CN
- China
- Prior art keywords
- node
- sample
- time
- power system
- voltage
- 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
- 238000000034 method Methods 0.000 title claims abstract description 48
- 230000001052 transient effect Effects 0.000 title claims abstract description 29
- 238000005259 measurement Methods 0.000 claims abstract description 58
- 230000035945 sensitivity Effects 0.000 claims abstract description 42
- 238000004088 simulation Methods 0.000 claims abstract description 34
- 238000003066 decision tree Methods 0.000 claims abstract description 29
- 238000007418 data mining Methods 0.000 claims abstract description 21
- 238000011156 evaluation Methods 0.000 claims abstract description 14
- 238000012417 linear regression Methods 0.000 claims abstract description 11
- 238000012544 monitoring process Methods 0.000 claims abstract description 10
- 230000001360 synchronised effect Effects 0.000 claims abstract description 8
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 7
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 230000006698 induction Effects 0.000 claims description 18
- 230000001687 destabilization Effects 0.000 claims description 16
- 238000012549 training Methods 0.000 claims description 15
- 230000005540 biological transmission Effects 0.000 claims description 12
- 238000002271 resection Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 9
- 230000000368 destabilizing effect Effects 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 230000000087 stabilizing effect Effects 0.000 claims description 5
- 238000002790 cross-validation Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 230000036962 time dependent Effects 0.000 claims 1
- 238000011158 quantitative evaluation Methods 0.000 abstract 1
- 230000015556 catabolic process Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000013097 stability assessment Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 230000003449 preventive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
Classifications
-
- 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
- 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
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
-
- 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
-
- 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/22—Flexible AC transmission systems [FACTS] or power factor or reactive power compensating or correcting units
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种电力系统负荷区域暂态电压稳定的评估方法,属于电力系统稳定分析评估领域。本发明方法以同步相量测量单元的量测数据为基础,通过大量仿真样本建立用于数据挖掘的初始样本库;通过区域内各节点的量化评估提取反映各节点稳定程度的特征量;通过多元线性回归的方法辨识得到反映区域网络内节点间相互影响关系的灵敏度系数;利用半监督聚类方法对所有样本进行标定;以决策树算法进行分类学习,得到决策树模型,将决策树模型用于在线监测,对电力系统负荷区域的整体暂态电压稳定状况实施评估。
Description
技术领域
本发明涉及一种电力系统负荷区域暂态电压稳定的评估方法,属于电力系统稳定分析评估领域。
背景技术
近年来世界范围内多起电压崩溃事故,使电力系统电压稳定性问题得到越来越多的关注。短期大扰动所引起的暂态电压失稳或崩溃尤为严重,在线监测时若能准确判断出区域性的暂态电压失稳并及时采取预防控制措施,将有利于遏制暂态电压失稳事故中部分节点所引起的大范围连锁性电压失稳或崩溃,提高区域整体稳定性。实际运行中常采用基于固定电压阈值的工程判据对暂态电压失稳进行判断,这类判据简便易行,但阈值的设定缺乏足够的理论支撑,其可靠性和选择性都难以保证。同时,针对区域整体的暂态电压稳定评估相关理论目前也还不甚完善,这给区域暂态电压稳定性的评估和判断带来了不少困难。
随着广域测量系统/同步相量测量单元(以下简称WAMS/PMU)的迅速发展和推广,目前有不少利用数据挖掘和机器学习的方法进行电力系统在线监测和安全稳定评估的方法。但传统的一些基于数据挖掘的评估方法是在所研究问题的理论基础相对完善的前提下开展的,少有涉及对区域暂态电压稳定评估的研究,且大多采用基于传统经验的工程判据对初始数据进行标定,并未充分发挥数据挖掘方法在探索数据规律方面的优势。当已知部分先验知识时,半监督聚类方法能高效、可靠地引导探索式知识发现过程。因此,若能充分利用数据挖掘的优势来对区域暂态电压稳定状况实施评估,将有望克服相关理论基础还不完善的困难,以数据中潜在的规律来保证评估结果的可靠性和选择性。
发明内容
本发明的目的是提出一种电力系统负荷区域暂态电压稳定的评估方法,以广域测量信息为基础,利用数据挖掘的方法来构建电力系统负荷区域暂态电压稳定的评估模型,从而对电力系统负荷区域暂态电压稳定状况进行可靠的监测和评估。
本发明提出的电力系统负荷区域暂态电压稳定的评估方法,方法包括以下步骤:
(1)设电力系统中的负荷区域内有k个节点,各节点序号分别为1,2,…,k,当电力系统处于稳态运行状态A时,从各节点的同步相量测量单元中,以时间间隔为ΔT采集该节点的电压、电流、有功功率和无功功率的实时量测数据,采样时长为T,设在稳态运行状态A时,电力系统中所有k个节点的负荷消耗的有功功率之和为PA,采用感应电动机与恒阻抗并联的电力系统综合负荷模型,分别建立电力系统负荷区域内各节点的等值负荷模型,并设电力系统各节点感应电动机消耗的有功功率与各节点的恒阻抗消耗的有功功率之比为Pm;
(2)建立在上述步骤(1)的稳态运行状态A下的电力系统仿真模型,对该电力系统仿真模型中的元件参数和运行条件进行调整,通过时域仿真生成大量仿真样本,具体过程如下:
(2-1)将电力系统中所有k个节点的负荷消耗的有功功率之和依次调整为PA、PA×0.8和PA×1.2;
(2-2)分别在步骤(2-1)的PA、PA×0.8和PA×1.2状态下,将电力系统各节点感应电动机消耗的有功功率与各节点的恒阻抗消耗的有功功率之比依次调整为Pm、0.5Pm、2Pm和4Pm;
(2-3)分别在步骤(2-2)的Pm、0.5Pm、2Pm和4Pm下,对电力系统负荷区域内的每条输电线路,在线路全长的0%、20%、40%、60%、80%和100%处设置三相短路故障;
(2-4)分别在步骤(2-3)的在线路全长的0%、20%、40%、60%、80%和100%处设置三相短路故障状态下,将上述每条输电线路两端三相短路故障的切除时间设置成两档,第一档为:近端切除时间0.05秒,远端切除时间0.1秒,第二档为近端切除时间0.35秒,远端切除时间0.4秒;
(2-5)分别在步骤(2-4)的两档三相短路故障切除时间状态下,以步骤(1)的时间间隔ΔT和采样时长T,对电力系统仿真模型进行N次时域仿真,分别记录每次时域仿真过程中电力系统负荷区域内所有k个节点的电压、无功功率和感应电动机滑差随时间的变化曲线,将三种曲线分别记为U、Q和s,并分别记录每次时域仿真过程中电力系统负荷区域内每条输电线路三相短路故障的发生时刻t0、切除时刻tc以及各节点的状态Z,将节点处于稳定状态记为Z=1,节点处于失稳状态记为Z=0,将一次时域仿真过程中记录的数据集合成一个样本,共N个样本形成一个初始样本库;
(3)从上述步骤(2)的初始样本库中获取电力系统中负荷区域内所有节点的感应电动机滑差的变化曲线s及各节点状态Z的记录,计算电力系统负荷区域内各节点的电压稳定度量指标,具体过程如下:
(3-1)从初始样本库中获取第m个样本的第i个节点上,在三相短路故障的发生时刻t0和三相短路故障的切除时刻tc,感应电动机滑差分别为s0和sc,其中1≤m≤N,1≤i≤k;
(3-2)从初始样本库中获取第m个样本的第i个节点的节点状态Z,对Z进行判断,若Z=1,则该节点为稳定节点,从三相短路故障切除时刻tc开始,节点上感应电动机滑差在随后的变化过程中第一次感应电动机滑差恢复到s0的时刻为t1,滑差恢复时间为Kr:
其中,Tr为稳定节点的监测时间窗,取值为1秒,K为临界状态的保守系数,取值为0.95,对Kr进行规范化处理,得到标准化的节点电压稳定度量指标KVS:
若Z=0,则该节点为失稳节点,从三相短路故障切除时刻tc开始,节点上感应电动机滑差在随后的变化过程中达到(2sc-s0)的时刻为t2,滑差增长时间为Kg:
其中,Tg为失稳节点的监测时间窗,取值为1秒,K为临界状态的保守系数,取值为0.95,对Kr进行规范化处理,得到标准化的节点电压稳定度量指标KVS:
(3-3)遍历初始样本库中第m个样本的所有k个节点,重复上述步骤(3-1)和(3-2),得到第m个样本中所有k个节点电压稳定度量指标;
(3-4)遍历初始样本库N个样本中所有样本,重复上述步骤(3-1)~(3-3),得到电力系统负荷区域内所有节点的电压稳定度量指标;
(4)从上述步骤(2)的初始样本库中获取电力系统中负荷区域内在三相短路故障切除时刻tc后的时间Δt内,所有节点的电压和无功功率随时间的变化曲线U和Q,通过多元线性回归分析的方法得到各节点间电压无功灵敏度系数,以此作为表征电力系统负荷区域内各节点之间相互影响的特征量,其中Δt=(n–1)×ΔT,ΔT为仿真时间间隔,n为所获取的变化曲线上的数据点数,具体过程如下:
(4-1)从上述步骤(2)的初始样本库中获取电力系统中负荷区域内在三相短路故障切除时刻tc后的时间Δt内,所有节点的电压和无功功率随时间的变化曲线U和Q,其中Δt=(n–1)×ΔT,ΔT为仿真时间间隔,n为所获取的变化曲线上的数据点数;
(4-2)从变化曲线U,得到初始样本库中第m个样本的第i个节点,电压随时间的变化曲线上各点数值依次为{u1i,u2i,…,uni},第j个节点的无功功率随时间的变化曲线Q上各点数值依次为{Q1j,Q2j,…,Qnj},其中j=1,2,…,k,建立第i个节点的多元线性回归模型如下:
其中, β为回归系数向量, ε为误差向量, In为n阶单位矩阵,E(ε)和COV(ε,ε)分别表示误差向量的期望和方差;
(4-3)定义在三相短路故障切除时刻tc后的时间Δt内,电力系统负荷区域内的第j个节点无功功率对第i个节点电压的灵敏度系数为j=1,2,…,k,以最小二乘法求解上述步骤(4-2)中的第i个节点的多元线性回归模型,得到回归系数β1i~βki,即为各节点无功功率对第i个节点电压的灵敏度系数;
(4-4)设上述步骤(4-3)的β1i~βki中的最大值、最小值分别为βimax和βimin,对灵敏度系数进行规范化处理,得到第j个节点无功功率对第i个节点电压的标准化灵敏度系数αij:
(4-5)遍历第m个样本中所有k个节点,重复上述步骤(4-1)~(4-4),得到第m个样本中所有k个节点电压无功功率灵敏度系数;
(4-6)遍历初始样本库中所有N个样本的所有k个节点,重复上述步骤(4-1)~(4-5),得到所有N个样本中所有k个节点的电压无功功率灵敏度系数;
(5)对上述步骤(3)和步骤(4)计算得到的所有N个样本的所有k个节点的电压稳定度量指标和电压无功功率灵敏度系数,利用半监督聚类方法进行处理,具体过程如下:
(5-1)将上述步骤(3)和步骤(4)计算得到的所有N个样本的所有k个节点电压稳定度量指标和电压无功功率灵敏度系数构成一个数据集C,该数据集C的维度为(k2+k),数据集C的特征属性包含k个节点稳定度量指标和k2个灵敏度系数;
(5-2)对上述步骤(5-1)的数据集C的第m个样本所有k个节点的电压稳定度量指标进行判断,若所有节点的电压稳定度量指标同时满足KVS>0,则判定电力系统负荷区域处于稳定状态,并将第m个样本的类标号记为R=1;若所有节点的电压稳定度量指标同时满足KVS<0,则判定电力系统负荷区域处于失稳状态,并将第m个样本的类标号记为R=0;若有k1个节点的电压稳定度量指标满足KVS>0,k2个节点的电压稳定度量指标满足KVS≤0,其中k1+k2=k,则对第m个样本的类标号不作任何标记,其中R=1表示区域稳定状态,R=0表示区域失稳状态;
(5-3)遍历数据集C中所有N个样本,重复上述步骤(5-2),直至完成所有N个样本的判断,得到电力系统负荷区域处于稳定状态和失稳状态的样本数目分别为N1和N2,并分别构成区域稳定子集C1和区域失稳子集C2,将区域稳定子集C1内各样本之间的关系和区域失稳子集C2内各样本之间的关系同时定义为必须连接约束关系,C1的样本与C2的样本之间的关系定义为不可连接约束关系;
(5-4)分别对上述步骤(5-3)的区域稳定子集C1和区域失稳子集C2中的样本各维度特征属性值进行均值计算,得到区域稳定子集C1和区域失稳子集C2中各维度特征属性值的算术平均值,将得到的多个算术平均值按维度顺序组合,得到相应的两个组合样本O1和O2,其中组合样本O1的类标号为R=1,组合样本O2的类标号为R=0;
(5-5)定义上述组合样本O1和组合样本O2为聚类过程的两个初始聚类中心,区域稳定子集C1和区域失稳子集C2为聚类过程的两个初始簇,对数据集C中的第m个样本进行判断,若该样本中存在必须连接约束关系或不可连接约束关系,则将该样本分配到该样本所属的初始簇中;若该样本中不存在必须连接约束关系和不可连接约束关系,则分别计算该样本与组合样本O1和组合样本O2的欧几里得距离,将欧几里得距离分别记为d1和d2,对d1和d2进行比较,若d1≤d2,则将该样本分配到区域稳定子集C1中去,并将该样本的类标号记为R=1,若d1>d2,则将该样本分配到区域失稳子集C2中去,并将该样本的类标号记为R=0;
(5-6)遍历数据集C中所有N个样本,重复上述步骤(5-5),得到所有样本的类标号,再重复步骤(5-4)的均值计算方法,对步骤(5-5)的聚类中心O1和聚类中心O2进行更新;
(5-7)重复上述步骤(5-5)和步骤(5-6),分别将本次聚类得到所有样本类标号与上一次聚类得到的相应样本类标号进行一一对比,若本次聚类的所有样本类标号与上一次聚类的相应样本类标号都一致,则表示数据集C的聚类已经稳定,得到数据集C的所有样本类标号,结束聚类迭代,进行步骤(5-8),若本次聚类的所有样本类标号与上一次聚类的相应样本类标号不全部一致,则重复上述步骤(5-5)和(5-6),继续迭代,直到本次聚类的所有样本类标号与上一次聚类的相应样本类标号都一致,得到数据集C的所有样本类标号;
(6)构成一个数据挖掘训练集,将上述步骤(5-1)的数据集C中的所有数据作为数据挖掘训练集的输入数据,将步骤(5-7)得到的数据集C的所有样本类标号作为数据挖掘训练集的输出数据,以决策树算法对数据挖掘训练集进行分类学习,得到一个决策树模型,以交叉验证方式计算决策树模型的分类准确率Pre,若Pre≥97%,表示决策树模型分类性能满足要求,进行步骤(7),若Pre<97%,以决策树算法重新对数据挖掘训练集进行分类学习,直到得到的决策树模型的分类准确率满足Pre≥97%,进行步骤(7);
(7)当电力系统遭遇短期大扰动时,电力系统负荷区域内各节点的同步相量测量单元实时采集该节点的电压、电流、有功功率、无功功率的实时量测数据,利用上述步骤(2)得到该节点的等值负荷模型中感应电动机滑差的实时变化曲线s和该节点状态Z;利用步骤(3)和步骤(4),得到电力系统遭遇短期大扰动下电力系统负荷区域内所有节点的电压稳定度量指标和电压无功功率灵敏度系数,将电压稳定度量指标和电压无功功率灵敏度系数输入到步骤(6)的决策树模型中,由决策树模型输出类标号R,作为电力系统区域暂态电压稳定的实时评估结果。
本发明提出的电力系统负荷区域暂态电压稳定的评估方法,其优点是,本发明方法利用节点稳定度量指标及基于多元线性回归分析方法的电压无功功率灵敏度系数全面地提取区域暂态电压稳定的特征属性,相比已有技术中基于主观运行经验的工程判据,采用基于约束的半监督聚类方法和决策树算法对电力系统运行数据中潜藏的规律进行探索,可对电力系统负荷区域暂态电压稳定状态实施更为准确的分类和评估,在实际运行中为电力系统负荷区域的在线监测和稳定控制提供可靠指导,从而避免电力系统负荷区域内因暂态事故造成的不必要的负荷停电和经济损失。
附图说明
图1为本发明方法涉及的中国电科院36节点测试系统单线结构示意图。
图2为本发明方法中分类学习得到的决策树模型。
具体实施方式
本发明提出的电力系统负荷区域暂态电压稳定的评估方法,包括以下步骤:
(1)本发明评估方法中涉及的电力系统单线结构示意图如图1所示,设电力系统中的负荷区域内有k个节点,图1所示为本发明的一个实施例,在一个36节点的测试系统中,针对虚线所围的电力系统负荷区域,该区域中包含节点16、18、19、20、21和29共6个节点,各节点序号分别为1,2,…,k,当电力系统处于稳态运行状态A时,从各节点的同步相量测量单元中,以时间间隔为ΔT采集该节点的电压、电流、有功功率和无功功率的实时量测数据,采样时长为T,设在稳态运行状态A时,电力系统中所有k个节点的负荷消耗的有功功率之和为PA,采用感应电动机与恒阻抗并联的电力系统综合负荷模型,分别建立电力系统负荷区域内各节点的等值负荷模型,并设电力系统各节点感应电动机消耗的有功功率与各节点的恒阻抗消耗的有功功率之比为Pm;
(2)建立在上述步骤(1)的稳态运行状态A下的电力系统仿真模型,对该电力系统仿真模型中的元件参数和运行条件进行调整,通过时域仿真生成大量仿真样本,具体过程如下:
(2-1)将电力系统中所有k个节点的负荷消耗的有功功率之和依次调整为PA、PA×0.8和PA×1.2;
(2-2)分别在步骤(2-1)的PA、PA×0.8和PA×1.2状态下,将电力系统各节点感应电动机消耗的有功功率与各节点的恒阻抗消耗的有功功率之比依次调整为Pm、0.5Pm、2Pm和4Pm;
(2-3)分别在步骤(2-2)的Pm、0.5Pm、2Pm和4Pm下,对电力系统负荷区域内的每条输电线路,在线路全长的0%、20%、40%、60%、80%和100%处设置三相短路故障;
(2-4)分别在步骤(2-3)的在线路全长的0%、20%、40%、60%、80%和100%处设置三相短路故障状态下,将上述每条输电线路两端三相短路故障的切除时间设置成两档,第一档为:近端切除时间0.05秒,远端切除时间0.1秒,第二档为近端切除时间0.35秒,远端切除时间0.4秒;
(2-5)分别在步骤(2-4)的两档三相短路故障切除时间状态下,以步骤(1)的时间间隔ΔT和采样时长T,对电力系统仿真模型进行N次时域仿真,分别记录每次时域仿真过程中电力系统负荷区域内所有k个节点的电压、无功功率和感应电动机滑差随时间的变化曲线,将三种曲线分别记为U、Q和s,并分别记录每次时域仿真过程中电力系统负荷区域内每条输电线路三相短路故障的发生时刻t0、切除时刻tc以及各节点的状态Z,将节点处于稳定状态记为Z=1,节点处于失稳状态记为Z=0,将一次时域仿真过程中记录的数据集合成一个样本,共N个样本形成一个初始样本库;
(3)从上述步骤(2)的初始样本库中获取电力系统中负荷区域内所有节点的感应电动机滑差的变化曲线s及各节点状态Z的记录,计算电力系统负荷区域内各节点的电压稳定度量指标,具体过程如下:
(3-1)从初始样本库中获取第m个样本的第i个节点上,在三相短路故障的发生时刻t0和三相短路故障的切除时刻tc,感应电动机滑差分别为s0和sc,其中1≤m≤N,1≤i≤k;
(3-2)从初始样本库中获取第m个样本的第i个节点的节点状态Z,对Z进行判断,若Z=1,则该节点为稳定节点,从三相短路故障切除时刻tc开始,节点上感应电动机滑差在随后的变化过程中第一次感应电动机滑差恢复到s0的时刻为t1,滑差恢复时间为Kr:
其中,Tr为稳定节点的监测时间窗,取值为1秒,K为临界状态的保守系数,取值为0.95,对Kr进行规范化处理,得到标准化的节点电压稳定度量指标KVS:
若Z=0,则该节点为失稳节点,从三相短路故障切除时刻tc开始,节点上感应电动机滑差在随后的变化过程中达到(2sc-s0)的时刻为t2,滑差增长时间为Kg:
其中,Tg为失稳节点的监测时间窗,取值为1秒,K为临界状态的保守系数,取值为0.95,对Kr进行规范化处理,得到标准化的节点电压稳定度量指标KVS:
(3-3)遍历初始样本库中第m个样本的所有k个节点,重复上述步骤(3-1)和(3-2),得到第m个样本中所有k个节点电压稳定度量指标;
(3-4)遍历初始样本库N个样本中所有样本,重复上述步骤(3-1)~(3-3),得到电力系统负荷区域内所有节点的电压稳定度量指标;
(4)从上述步骤(2)的初始样本库中获取电力系统中负荷区域内在三相短路故障切除时刻tc后的时间Δt内,所有节点的电压和无功功率随时间的变化曲线U和Q,通过多元线性回归分析的方法得到各节点间电压无功灵敏度系数,以此作为表征电力系统负荷区域内各节点之间相互影响的特征量,其中Δt=(n–1)×ΔT,ΔT为仿真时间间隔,n为所获取的变化曲线上的数据点数,具体过程如下:
(4-1)从上述步骤(2)的初始样本库中获取电力系统中负荷区域内在三相短路故障切除时刻tc后的时间Δt内,所有节点的电压和无功功率随时间的变化曲线U和Q,其中Δt=(n–1)×ΔT,ΔT为仿真时间间隔,n为所获取的变化曲线上的数据点数;
(4-2)从变化曲线U,得到初始样本库中第m个样本的第i个节点,电压随时间的变化曲线上各点数值依次为{u1i,u2i,…,uni},第j个节点的无功功率随时间的变化曲线Q上各点数值依次为{Q1j,Q2j,…,Qnj},其中j=1,2,…,k,建立第i个节点的多元线性回归模型如下:
其中, β为回归系数向量, ε为误差向量, In为n阶单位矩阵,E(ε)和COV(ε,ε)分别表示误差向量的期望和方差;
(4-3)定义在三相短路故障切除时刻tc后的时间Δt内,电力系统负荷区域内的第j个节点无功功率对第i个节点电压的灵敏度系数为j=1,2,…,k,以最小二乘法求解上述步骤(4-2)中的第i个节点的多元线性回归模型,得到回归系数β1i~βki,即为各节点无功功率对第i个节点电压的灵敏度系数;
(4-4)设上述步骤(4-3)的β1i~βki中的最大值、最小值分别为βimax和βimin,对灵敏度系数进行规范化处理,得到第j个节点无功功率对第i个节点电压的标准化灵敏度系数αij:
(4-5)遍历第m个样本中所有k个节点,重复上述步骤(4-1)~(4-4),得到第m个样本中所有k个节点电压无功功率灵敏度系数;
(4-6)遍历初始样本库中所有N个样本的所有k个节点,重复上述步骤(4-1)~(4-5),得到所有N个样本中所有k个节点的电压无功功率灵敏度系数;
(5)对上述步骤(3)和步骤(4)计算得到的所有N个样本的所有k个节点的电压稳定度量指标和电压无功功率灵敏度系数,利用半监督聚类方法进行处理,具体过程如下:
(5-1)将上述步骤(3)和步骤(4)计算得到的所有N个样本的所有k个节点电压稳定度量指标和电压无功功率灵敏度系数构成一个数据集C,该数据集C的维度为(k2+k),数据集C的特征属性包含k个节点稳定度量指标和k2个灵敏度系数;
(5-2)对上述步骤(5-1)的数据集C的第m个样本所有k个节点的电压稳定度量指标进行判断,若所有节点的电压稳定度量指标同时满足KVS>0,则判定电力系统负荷区域处于稳定状态,并将第m个样本的类标号记为R=1;若所有节点的电压稳定度量指标同时满足KVS<0,则判定电力系统负荷区域处于失稳状态,并将第m个样本的类标号记为R=0;若有k1个节点的电压稳定度量指标满足KVS>0,k2个节点的电压稳定度量指标满足KVS≤0,其中k1+k2=k,则对第m个样本的类标号不作任何标记,其中R=1表示区域稳定状态,R=0表示区域失稳状态;
(5-3)遍历数据集C中所有N个样本,重复上述步骤(5-2),直至完成所有N个样本的判断,得到电力系统负荷区域处于稳定状态和失稳状态的样本数目分别为N1和N2,并分别构成区域稳定子集C1和区域失稳子集C2,将区域稳定子集C1内各样本之间的关系和区域失稳子集C2内各样本之间的关系同时定义为必须连接约束关系,C1的样本与C2的样本之间的关系定义为不可连接约束关系;
(5-4)分别对上述步骤(5-3)的区域稳定子集C1和区域失稳子集C2中的样本各维度特征属性值进行均值计算,得到区域稳定子集C1和区域失稳子集C2中各维度特征属性值的算术平均值,将得到的多个算术平均值按维度顺序组合,得到相应的两个组合样本O1和O2,其中组合样本O1的类标号为R=1,组合样本O2的类标号为R=0;
(5-5)定义上述组合样本O1和组合样本O2为聚类过程的两个初始聚类中心,区域稳定子集C1和区域失稳子集C2为聚类过程的两个初始簇,对数据集C中的第m个样本进行判断,若该样本中存在必须连接约束关系或不可连接约束关系,则将该样本分配到该样本所属的初始簇中;若该样本中不存在必须连接约束关系和不可连接约束关系,则分别计算该样本与组合样本O1和组合样本O2的欧几里得距离,将欧几里得距离分别记为d1和d2,对d1和d2进行比较,若d1≤d2,则将该样本分配到区域稳定子集C1中去,并将该样本的类标号记为R=1,若d1>d2,则将该样本分配到区域失稳子集C2中去,并将该样本的类标号记为R=0;
(5-6)遍历数据集C中所有N个样本,重复上述步骤(5-5),得到所有样本的类标号,再重复步骤(5-4)的均值计算方法,对步骤(5-5)的聚类中心O1和聚类中心O2进行更新;
(5-7)重复上述步骤(5-5)和步骤(5-6),分别将本次聚类得到所有样本类标号与上一次聚类得到的相应样本类标号进行一一对比,若本次聚类的所有样本类标号与上一次聚类的相应样本类标号都一致,则表示数据集C的聚类已经稳定,得到数据集C的所有样本类标号,结束聚类迭代,进行步骤(5-8),若本次聚类的所有样本类标号与上一次聚类的相应样本类标号不全部一致,则重复上述步骤(5-5)和(5-6),继续迭代,直到本次聚类的所有样本类标号与上一次聚类的相应样本类标号都一致,得到数据集C的所有样本类标号;
(6)构成一个数据挖掘训练集,将上述步骤(5-1)的数据集C中的所有数据作为数据挖掘训练集的输入数据,将步骤(5-7)得到的数据集C的所有样本类标号作为数据挖掘训练集的输出数据,以决策树算法对数据挖掘训练集进行分类学习,得到一个决策树模型,本发明的一个实施例中决策树模型如图2所示,图2所示的决策树模型中终端节点的标号1和0表示输出的类标号R,其中类标号1代表区域稳定类别,类标号0代表区域失稳类别,特征属性KVS_i表示第i个节点的电压稳定度量指标,αij表示第j个节点无功功率对第i个节点电压的灵敏度系数,以交叉验证方式计算得到决策树模型的分类准确率Pre=98.9%,Pre≥97%,表示决策树模型分类性能满足要求,进行步骤(7);
(7)当电力系统遭遇短期大扰动时,电力系统负荷区域内各节点的同步相量测量单元实时采集该节点的电压、电流、有功功率、无功功率的实时量测数据,利用上述步骤(2)得到该节点的等值负荷模型中感应电动机滑差的实时变化曲线s和该节点状态Z;利用步骤(3)和步骤(4),得到电力系统遭遇短期大扰动下电力系统负荷区域内所有节点的电压稳定度量指标和电压无功功率灵敏度系数,将电压稳定度量指标和电压无功功率灵敏度系数输入到步骤(6)的决策树模型中,由决策树模型输出类标号R,作为电力系统区域暂态电压稳定的实时评估结果。若R=1,说明电力系统负荷区域可维持暂态电压稳定,若R=0,说明电力系统负荷区域将发生暂态电压失稳,需立即采取稳定控制措施,以防止暂态电压失稳事故进一步扩大,从而避免不必要的负荷停电和经济损失。
Claims (1)
1.一种电力系统负荷区域暂态电压稳定的评估方法,其特征在于该方法包括以下步骤:
(1)设电力系统中的负荷区域内有k个节点,各节点序号分别为1,2,…,k,当电力系统处于稳态运行状态A时,从各节点的同步相量测量单元中,以时间间隔为ΔT采集该节点的电压、电流、有功功率和无功功率的实时量测数据,采样时长为T,设在稳态运行状态A时,电力系统中所有k个节点的负荷消耗的有功功率之和为PA,采用感应电动机与恒阻抗并联的电力系统综合负荷模型,分别建立电力系统负荷区域内各节点的等值负荷模型,并设电力系统各节点感应电动机消耗的有功功率与各节点的恒阻抗消耗的有功功率之比为Pm;
(2)建立在上述步骤(1)的稳态运行状态A下的电力系统仿真模型,对该电力系统仿真模型中的元件参数和运行条件进行调整,通过时域仿真生成大量仿真样本,具体过程如下:
(2-1)将电力系统中所有k个节点的负荷消耗的有功功率之和依次调整为PA、PA×0.8和PA×1.2;
(2-2)分别在步骤(2-1)的PA、PA×0.8和PA×1.2状态下,将电力系统各节点感应电动机消耗的有功功率与各节点的恒阻抗消耗的有功功率之比依次调整为Pm、0.5Pm、2Pm和4Pm;
(2-3)分别在步骤(2-2)的Pm、0.5Pm、2Pm和4Pm下,对电力系统负荷区域内的每条输电线路,在线路全长的0%、20%、40%、60%、80%和100%处设置三相短路故障;
(2-4)分别在步骤(2-3)的在线路全长的0%、20%、40%、60%、80%和100%处设置三相短路故障状态下,将上述每条输电线路两端三相短路故障的切除时间设置成两档,第一档为:近端切除时间0.05秒,远端切除时间0.1秒,第二档为近端切除时间0.35秒,远端切除时间0.4秒;
(2-5)分别在步骤(2-4)的两档三相短路故障切除时间状态下,以步骤(1)的时间间隔ΔT和采样时长T,对电力系统仿真模型进行N次时域仿真,分别记录每次时域仿真过程中电力系统负荷区域内所有k个节点的电压、无功功率和感应电动机滑差随时间的变化曲线,将三种曲线分别记为U、Q和s,并分别记录每次时域仿真过程中电力系统负荷区域内每条输电线路三相短路故障的发生时刻t0、切除时刻tc以及各节点的状态Z,将节点处于稳定状态记为Z=1,节点处于失稳状态记为Z=0,将一次时域仿真过程中记录的数据集合成一个样本,共N个样本形成一个初始样本库;
(3)从上述步骤(2)的初始样本库中获取电力系统中负荷区域内所有节点的感应电动机滑差的变化曲线s及各节点状态Z的记录,计算电力系统负荷区域内各节点的电压稳定度量指标,具体过程如下:
(3-1)从初始样本库中获取第m个样本的第i个节点上,在三相短路故障的发生时刻t0和三相短路故障的切除时刻tc,感应电动机滑差分别为s0和sc,其中1≤m≤N,1≤i≤k;
(3-2)从初始样本库中获取第m个样本的第i个节点的节点状态Z,对Z进行判断,若Z=1,则该节点为稳定节点,从三相短路故障切除时刻tc开始,节点上感应电动机滑差在随后的变化过程中第一次感应电动机滑差恢复到s0的时刻为t1,滑差恢复时间为Kr:
其中,Tr为稳定节点的监测时间窗,取值为1秒,K为临界状态的保守系数,取值为0.95,对Kr进行规范化处理,得到标准化的节点电压稳定度量指标KVS:
若Z=0,则该节点为失稳节点,从三相短路故障切除时刻tc开始,节点上感应电动机滑差在随后的变化过程中达到(2sc-s0)的时刻为t2,滑差增长时间为Kg:
其中,Tg为失稳节点的监测时间窗,取值为1秒,K为临界状态的保守系数,取值为0.95,对Kr进行规范化处理,得到标准化的节点电压稳定度量指标KVS:
(3-3)遍历初始样本库中第m个样本的所有k个节点,重复上述步骤(3-1)和(3-2),得到第m个样本中所有k个节点电压稳定度量指标;
(3-4)遍历初始样本库N个样本中所有样本,重复上述步骤(3-1)~(3-3),得到电力系统负荷区域内所有节点的电压稳定度量指标;
(4)从上述步骤(2)的初始样本库中获取电力系统中负荷区域内在三相短路故障切除时刻tc后的时间Δt内,所有节点的电压和无功功率随时间的变化曲线U和Q,通过多元线性回归分析的方法得到各节点间电压无功灵敏度系数,以此作为表征电力系统负荷区域内各节点之间相互影响的特征量,其中Δt=(n–1)×ΔT,ΔT为仿真时间间隔,n为所获取的变化曲线上的数据点数,具体过程如下:
(4-1)从上述步骤(2)的初始样本库中获取电力系统中负荷区域内在三相短路故障切除时刻tc后的时间Δt内,所有节点的电压和无功功率随时间的变化曲线U和Q,其中Δt=(n–1)×ΔT,ΔT为仿真时间间隔,n为所获取的变化曲线上的数据点数;
(4-2)从变化曲线U,得到初始样本库中第m个样本的第i个节点,电压随时间的变化曲线上各点数值依次为{u1i,u2i,…,uni},第j个节点的无功功率随时间的变化曲线Q上各点数值依次为{Q1j,Q2j,…,Qnj},其中j=1,2,…,k,建立第i个节点的多元线性回归模型如下:
其中,β为回归系数向量,ε为误差向量,In为n阶单位矩阵,E(ε)和COV(ε,ε)分别表示误差向量的期望和方差;
(4-3)定义在三相短路故障切除时刻tc后的时间Δt内,电力系统负荷区域内的第j个节点无功功率对第i个节点电压的灵敏度系数为j=1,2,…,k,,其中Ui是第i个节点在电压随时间的变化曲线U上的值,Qj是第i个节点在无功功率随时间的变化曲线Q上的值,以最小二乘法求解上述步骤(4-2)中的第i个节点的多元线性回归模型,得到回归系数β1i~βki,即为各节点无功功率对第i个节点电压的灵敏度系数;
(4-4)设上述步骤(4-3)的β1i~βki中的最大值、最小值分别为βimax和βimin,对灵敏度系数进行规范化处理,得到第j个节点无功功率对第i个节点电压的标准化灵敏度系数αji:
(4-5)遍历第m个样本中所有k个节点,重复上述步骤(4-1)~(4-4),得到第m个样本中所有k个节点电压无功功率灵敏度系数;
(4-6)遍历初始样本库中所有N个样本的所有k个节点,重复上述步骤(4-1)~(4-5),得到所有N个样本中所有k个节点的电压无功功率灵敏度系数;
(5)对上述步骤(3)和步骤(4)计算得到的所有N个样本的所有k个节点的电压稳定度量指标和电压无功功率灵敏度系数,利用半监督聚类方法进行处理,具体过程如下:
(5-1)将上述步骤(3)和步骤(4)计算得到的所有N个样本的所有k个节点电压稳定度量指标和电压无功功率灵敏度系数构成一个数据集C,该数据集C的维度为(k2+k),数据集C的特征属性包含k个节点稳定度量指标和k2个灵敏度系数;
(5-2)对上述步骤(5-1)的数据集C的第m个样本所有k个节点的电压稳定度量指标进行判断,若所有节点的电压稳定度量指标同时满足KVS>0,则判定电力系统负荷区域处于稳定状态,并将第m个样本的类标号记为R=1;若所有节点的电压稳定度量指标同时满足KVS<0,则判定电力系统负荷区域处于失稳状态,并将第m个样本的类标号记为R=0;若有k1个节点的电压稳定度量指标满足KVS>0,k2个节点的电压稳定度量指标满足KVS≤0,其中k1+k2=k,则对第m个样本的类标号不作任何标记,其中R=1表示区域稳定状态,R=0表示区域失稳状态;
(5-3)遍历数据集C中所有N个样本,重复上述步骤(5-2),直至完成所有N个样本的判断,得到电力系统负荷区域处于稳定状态和失稳状态的样本数目分别为N1和N2,并分别构成区域稳定子集C1和区域失稳子集C2,将区域稳定子集C1内各样本之间的关系和区域失稳子集C2内各样本之间的关系同时定义为必须连接约束关系,C1的样本与C2的样本之间的关系定义为不可连接约束关系;
(5-4)分别对上述步骤(5-3)的区域稳定子集C1和区域失稳子集C2中的样本各维度特征属性值进行均值计算,得到区域稳定子集C1和区域失稳子集C2中各维度特征属性值的算术平均值,将得到的多个算术平均值按维度顺序组合,得到相应的两个组合样本O1和O2,其中组合样本O1的类标号为R=1,组合样本O2的类标号为R=0;
(5-5)定义上述组合样本O1和组合样本O2为聚类过程的两个初始聚类中心,区域稳定子集C1和区域失稳子集C2为聚类过程的两个初始簇,对数据集C中的第m个样本进行判断,若该样本中存在必须连接约束关系或不可连接约束关系,则将该样本分配到该样本所属的初始簇中;若该样本中不存在必须连接约束关系和不可连接约束关系,则分别计算该样本与组合样本O1和组合样本O2的欧几里得距离,将欧几里得距离分别记为d1和d2,对d1和d2进行比较,若d1≤d2,则将该样本分配到区域稳定子集C1中去,并将该样本的类标号记为R=1,若d1>d2,则将该样本分配到区域失稳子集C2中去,并将该样本的类标号记为R=0;
(5-6)遍历数据集C中所有N个样本,重复上述步骤(5-5),得到所有样本的类标号,再重复步骤(5-4)的均值计算方法,对步骤(5-5)的聚类中心O1和聚类中心O2进行更新;
(5-7)重复上述步骤(5-5)和步骤(5-6),分别将本次聚类得到所有样本类标号与上一次聚类得到的相应样本类标号进行一一对比,若本次聚类的所有样本类标号与上一次聚类的相应样本类标号都一致,则表示数据集C的聚类已经稳定,得到数据集C的所有样本类标号,结束聚类迭代,进行步骤(5-8),若本次聚类的所有样本类标号与上一次聚类的相应样本类标号不全部一致,则重复上述步骤(5-5)和(5-6),继续迭代,直到本次聚类的所有样本类标号与上一次聚类的相应样本类标号都一致,得到数据集C的所有样本类标号;
(6)构成一个数据挖掘训练集,将上述步骤(5-1)的数据集C中的所有数据作为数据挖掘训练集的输入数据,将步骤(5-7)得到的数据集C的所有样本类标号作为数据挖掘训练集的输出数据,以决策树算法对数据挖掘训练集进行分类学习,得到一个决策树模型,以交叉验证方式计算决策树模型的分类准确率Pre,若Pre≥97%,表示决策树模型分类性能满足要求,进行步骤(7),若Pre<97%,以决策树算法重新对数据挖掘训练集进行分类学习,直到得到的决策树模型的分类准确率满足Pre≥97%,进行步骤(7);
(7)当电力系统遭遇短期大扰动时,电力系统负荷区域内各节点的同步相量测量单元实时采集该节点的电压、电流、有功功率、无功功率的实时量测数据,利用上述步骤(2)得到该节点的等值负荷模型中感应电动机滑差的实时变化曲线s和该节点状态Z;利用步骤(3)和步骤(4),得到电力系统遭遇短期大扰动下电力系统负荷区域内所有节点的电压稳定度量指标和电压无功功率灵敏度系数,将电压稳定度量指标和电压无功功率灵敏度系数输入到步骤(6)的决策树模型中,由决策树模型输出类标号R,作为电力系统区域暂态电压稳定的实时评估结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510025577.2A CN104617574B (zh) | 2015-01-19 | 2015-01-19 | 一种电力系统负荷区域暂态电压稳定的评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510025577.2A CN104617574B (zh) | 2015-01-19 | 2015-01-19 | 一种电力系统负荷区域暂态电压稳定的评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104617574A CN104617574A (zh) | 2015-05-13 |
CN104617574B true CN104617574B (zh) | 2016-08-17 |
Family
ID=53151914
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510025577.2A Active CN104617574B (zh) | 2015-01-19 | 2015-01-19 | 一种电力系统负荷区域暂态电压稳定的评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104617574B (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105139289B (zh) * | 2015-09-06 | 2018-10-19 | 清华大学 | 一种基于错分代价分类学习的电网暂态电压稳定评估方法 |
CN106055761B (zh) * | 2016-05-25 | 2019-07-12 | 清华大学 | 一种基于决策树分类的负荷模型参数校验方法 |
CN106208052B (zh) * | 2016-08-22 | 2018-08-17 | 国家电网公司 | 一种基于暂态电压稳定极限测试的电网暂态电压稳定薄弱点辨识方法 |
JP6715740B2 (ja) * | 2016-10-13 | 2020-07-01 | 株式会社日立製作所 | 電力系統の潮流監視装置、電力系統安定化装置および電力系統の潮流監視方法 |
CN106842021B (zh) * | 2016-12-29 | 2022-03-25 | 国家电网公司 | 一种暂态电压稳定性的判别方法和装置 |
CN106843941B (zh) * | 2016-12-31 | 2019-02-05 | Oppo广东移动通信有限公司 | 信息处理方法、装置和计算机设备 |
CN107482621B (zh) * | 2017-08-02 | 2019-09-27 | 清华大学 | 一种基于电压时序轨迹的电力系统暂态电压稳定评估方法 |
CN107679768B (zh) * | 2017-10-25 | 2021-04-13 | 中国南方电网有限责任公司 | 一种基于电网实时数据的态势感知系统及其构建方法 |
CN108964125A (zh) * | 2018-07-19 | 2018-12-07 | 长沙理工大学 | 提高主网功角安全鲁棒性的大型风电接入点选择方法 |
CN108988347B (zh) * | 2018-08-01 | 2020-09-22 | 中国南方电网有限责任公司 | 一种电网暂态电压稳定样本集类别失衡的调节方法及系统 |
CN109684673B (zh) * | 2018-12-03 | 2023-03-24 | 三峡大学 | 一种电力系统暂态稳定结果的特征提取和聚类分析方法 |
CN110896217B (zh) * | 2019-09-04 | 2022-04-15 | 中国电力科学研究院有限公司 | 一种提升受端电网电压稳定性的网源协调控制方法及装置 |
CN111722057A (zh) * | 2020-06-10 | 2020-09-29 | 国网辽宁省电力有限公司电力科学研究院 | 一种基于聚类的试验型微电网负荷特性分类方法 |
CN111756034B (zh) * | 2020-06-12 | 2022-04-08 | 清华大学 | 基于图时空网络的电力系统暂态电压稳定评估方法 |
CN112653131B (zh) * | 2020-12-10 | 2022-09-06 | 国电南瑞科技股份有限公司 | 暂态电压安全稳定裕度估算方法、装置及系统 |
CN113904384B (zh) * | 2021-11-09 | 2023-06-09 | 国网四川省电力公司电力科学研究院 | 一种基于梯度提升机的电网暂态稳定协调控制方法及系统 |
CN114552575B (zh) * | 2022-03-11 | 2023-02-14 | 大连理工大学 | 一种根据异步电动机负荷占比调整系统热备用功率的调度控制策略 |
CN115360719B (zh) * | 2022-08-30 | 2024-04-12 | 东北大学秦皇岛分校 | 基于plnn的电力系统短期电压稳定性评估方法 |
CN117215205B (zh) * | 2023-11-09 | 2024-02-06 | 国网经济技术研究院有限公司 | 基于决策树和iss理论的直流系统控制参数分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101699448B (zh) * | 2009-10-26 | 2011-06-29 | 清华大学 | 一种电力系统暂态稳定的分布式仿真方法 |
-
2015
- 2015-01-19 CN CN201510025577.2A patent/CN104617574B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104617574A (zh) | 2015-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104617574B (zh) | 一种电力系统负荷区域暂态电压稳定的评估方法 | |
CN112508442B (zh) | 基于自动化和可解释机器学习的暂态稳定评估方法及系统 | |
CN110458230A (zh) | 一种基于多判据融合的配变用采数据异常甄别方法 | |
CN105139289B (zh) | 一种基于错分代价分类学习的电网暂态电压稳定评估方法 | |
CN113702895B (zh) | 一种电压互感器误差状态在线定量评估方法 | |
CN111401749A (zh) | 一种基于随机森林与极限学习回归的动态安全评估方法 | |
CN110428005B (zh) | 一种基于伞式算法的电力系统动态安全误分类约束方法 | |
CN112051481A (zh) | 一种基于lstm的交直流混联电网故障区域诊断方法和系统 | |
CN107482621A (zh) | 一种基于电压时序轨迹的电力系统暂态电压稳定评估方法 | |
CN105606914A (zh) | 一种基于iwo-elm的航空功率变换器故障诊断方法 | |
CN102832617A (zh) | 一种基于精确化模式判别的大电网暂态稳定分析方法 | |
CN112731260A (zh) | 基于概念漂移识别的电压互感器误差状态在线评估方法 | |
CN112816881A (zh) | 电池压差异常检测方法、装置及计算机存储介质 | |
CN112684295A (zh) | 一种基于相似离度的高渗透率下配电网故障选线方法及系统 | |
CN111553112A (zh) | 一种基于深度置信网络的电力系统故障辨识方法及装置 | |
CN115659833A (zh) | 基于bp神经网络的电力网络节点脆弱性评估方法 | |
CN114036647A (zh) | 一种基于实车数据的动力电池安全风险评估方法 | |
Wilches-Bernal et al. | A machine learning-based method using the dynamic mode decomposition for fault location and classification | |
CN117171702A (zh) | 一种基于深度学习的多模态电网故障检测方法和系统 | |
CN116578833A (zh) | 基于优化随机森林模型的igbt模块老化故障诊断系统 | |
CN112926686B (zh) | 基于brb和lstm模型的电力大数据用电异常检测方法及装置 | |
Jana et al. | Transmission line fault pattern recognition using decision tree based smart fault classifier in a large power network | |
CN114091549A (zh) | 一种基于深度残差网络的设备故障诊断方法 | |
CN114116832A (zh) | 一种基于数据驱动的配电网异常识别方法 | |
CN106646106B (zh) | 基于变点探测技术的电网故障检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |