CN110957723A - 一种数据驱动的在线快速评估电网暂态电压安全性的方法 - Google Patents
一种数据驱动的在线快速评估电网暂态电压安全性的方法 Download PDFInfo
- Publication number
- CN110957723A CN110957723A CN201911219155.3A CN201911219155A CN110957723A CN 110957723 A CN110957723 A CN 110957723A CN 201911219155 A CN201911219155 A CN 201911219155A CN 110957723 A CN110957723 A CN 110957723A
- Authority
- CN
- China
- Prior art keywords
- bus
- power grid
- rlt
- voltage
- operation mode
- 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
- 230000001052 transient effect Effects 0.000 title claims abstract description 89
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000011156 evaluation Methods 0.000 claims abstract description 19
- 230000001502 supplementing effect Effects 0.000 claims abstract description 10
- 230000010355 oscillation Effects 0.000 claims description 7
- 238000011084 recovery Methods 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims description 4
- 238000004220 aggregation Methods 0.000 claims description 3
- 230000002776 aggregation Effects 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 abstract description 3
- 238000002347 injection Methods 0.000 description 6
- 239000007924 injection Substances 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 3
- 238000010248 power generation Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
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
-
- 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
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提出一种数据驱动的在线快速评估电网暂态电压安全性的方法,属于电力***运行控制技术领域。该方法首先,设定电网运行状态量的限值,将符合电网运行状态量限值的电网运行方式存储到离线数据库;计算基于母线电压轨迹的评估电网暂态电压安全性指标TVSI;仿真分析离线数据库中每种电网运行方式对应的TVSI,并将结果补充到离线数据库中;从更新后的离线数据库中选取当前电网运行方式的临近电网运行方式,用临近电网运行方式拟合当前电网运行方式的暂态电压安全性指标,得到当前电网运行方式的暂态电压安全性评估结果。本发明可以快速评估电网的暂态电压安全性,可以用于在线分析和无功优化问题的稳定约束判定。
Description
技术领域
本发明涉及一种数据驱动的在线快速评估电网暂态电压安全性的方法,属于电力***运行控制技术领域。
背景技术
近年来,我国建设并投运了多条大容量远距离输电线路,形成了一批典型的送端区域电网(以下简称:送端电网)和受端区域电网(以下简称:受端电网)。受端电网依赖外部输送有功功率,其区域内部的发电量小于负荷量。因为传统的同步发电机可以用作无功源,因此受端电网的无功源容量相对较少,导致其内部的电压稳定性较差。
在现代电网中,接入了大量的间歇式可再生能源发电设备(例如:风力发电、光伏发电等)、动态响应负荷(例如:空调)以及直流换流站等动态响应特性丰富的元件,因而其暂态电压安全问题也值得关注。
对于大电网的暂态电压安全研究,目前大多依靠数值仿真软件。但仿真分析耗时较长,难以用于在线分析、无功优化问题的稳定约束判定等研究。
发明内容
本发明的目的是为克服已有技术的不足之处,提出一种数据驱动的在线快速评估电网暂态电压安全性的方法。本发明可以快速评估电网的暂态电压安全性,可以用于在线分析和无功优化问题的稳定约束判定。
本发明提出一种数据驱动的在线快速评估电网暂态电压安全性的方法,其特征在于,该方法包括以下步骤:
(1)设定电网运行状态量的限值,包括:电网控制状态量的限值和电网潮流结果状态量的限值;随机生成电网运行方式,分别计算每种电网运行方式对应的潮流结果,将符合电网运行状态量限值的电网运行方式及其对应的电网的控制状态量集合存储到离线数据库中;具体步骤如下:
(1-1)设定电网控制状态量的限值:
定义ΩGrid为电网中母线的集合,将ΩGrid分为三个子集合分为记为ΩPV、ΩPQ、ΩVθ;其中,ΩPV中的母线为PV节点,ΩPQ中的母线为PQ节点,ΩVθ中的母线为Vθ节点;
定义CSV为电网控制状态量集合,表达式如下:
SV sp={Vi sp|i∈ΩPV∪ΩVθ}
Sθ sp={θi sp|i∈ΩVθ}
Pg,i sp,min≤Pg,i sp≤Pg,i sp,max,i∈ΩPV
Vi sp,min≤Vi sp≤Vi sp,max,i∈ΩPV∪ΩVθ
Pl,i sp,min≤Pl,i sp≤Pl,i sp,max,i∈ΩGrid
Ql,i sp,min≤Ql,i sp≤Ql,i sp,max,i∈ΩGrid
其中,为ΩPV中母线的发电机有功出力设定值的集合,SV sp为ΩPV和ΩVθ中母线的电压幅值设定值的集合,为电网中母线的有功负荷量设定值的集合,为电网中母线的无功负荷量设定值的集合,Sθ sp为ΩVθ中母线的电压相角设定值的集合;Pg,i sp为母线i的发电机有功出力设定值,Vi sp为母线i的电压幅值设定值,Pl,i sp为母线i的有功负荷量设定值,Ql,i sp为母线i的无功负荷量设定值,θi sp为母线i的电压相角设定值;Pg,i sp,min是母线i的发电机无功出力设定值下限,Pg,i sp,max是母线i的发电机无功出力设定值上限,Vi sp,min是母线i的电压幅值设定值下限,Vi sp,max是母线i的电压幅值设定值上限,Pl,i sp,min是母线i的有功负荷量设定值下限,Pl,i sp,max是母线i的有功负荷量设定值上限,Ql,i sp,min是母线i的无功负荷量设定值下限,Ql,i sp,max是母线i的无功负荷量设定值上限,是电网中母线电压相角设定值的最大值,是电网中母线电压相角设定值的最小值;
(1-2)设定电网潮流结果状态量的限值;
定义RSV为电网的潮流结果状态量集合,表达式如下:
SV Rlt={Vi Rlt|i∈ΩGrid}
Sθ Rlt={θi Rlt|i∈ΩGrid}
Pg,i Rlt,min≤Pg,i Rlt≤Pg,i Rlt,max,i∈ΩPV∪ΩVθ
Qg,i Rlt,min≤Qg,i Rlt≤Qg,i Rlt,max,i∈ΩPV∪ΩVθ
Pl,i Rlt,min≤Pl,i RltRlt≤Pl,i Rlt,max,i∈ΩGrid
Ql,i Rlt,min≤Ql,i Rlt≤Ql,i Rlt,max,i∈ΩGrid
Vi Rlt,min≤Vi Rlt≤Vi Rlt,max,i∈ΩGrid
其中,为ΩPV和ΩVθ中母线的发电机有功出力潮流结果的集合,为ΩPV和ΩVθ中母线的发电机无功出力潮流结果的集合,为电网中母线的有功负荷量潮流结果的集合,为电网中母线的无功负荷量潮流结果的集合,SV Rlt为电网中母线的电压幅值潮流结果的集合,Sθ Rlt为电网中母线的电压相角潮流结果的集合;Pg,i Rlt为母线i的发电机有功出力潮流结果,Qg,i Rlt为母线i的发电机无功出力潮流结果,Pl,i Rlt为母线i的有功负荷量潮流结果,Ql,i Rlt为母线i的发电机无功负荷量潮流结果,Vi Rlt为母线i的电压幅值潮流结果,θi Rlt为母线i的电压相角潮流结果;Pg,i Rlt,min是母线i的发电机有功出力潮流结果下限,Pg,i Rlt,max是母线i的发电机有功出力潮流结果上限,Qg,i Rlt,min是母线i的发电机无功出力潮流结果下限,Qg,i Rlt,max是母线i的发电机无功出力潮流结果上限,Pl,i Rlt,min是母线i的有功负荷量潮流结果下限,Pl,i Rlt,max是母线i的有功负荷量潮流结果上限,Ql,i Rlt,min是母线i的无功负荷量潮流结果下限,Ql,i Rlt,max是母线i的无功负荷量潮流结果上限,Vi Rlt,min是母线i的电压幅值潮流结果下限,Vi Rlt,max是母线i的电压幅值潮流结果上限,是电网中母线电压相角潮流结果的最大值,是电网中母线电压相角潮流结果的最小值;
(1-3)随机生成若干电网运行方式,筛选出潮流结果状态量和控制状态量均符合对应限值的电网运行方式,并将筛选后的电网运行方式对应的控制状态量结果作为一个样本存储到离线数据库中;
离线数据库的存储格式如下:
DB={Case1,...,Casen}
Casei={i,CSVi}
其中,Casei是离线数据库中的第i个样本,CSVi是该样本对应的控制状态量集合,DB是离线数据库,下标n是离线数据库中的样本数量;
(2)计算基于母线电压轨迹的评估电网暂态电压安全性指标,具体步骤如下:
(2-1)计算任一母线在任一故障下的暂态电压安全性指标,具体步骤如下:
(2-1-1)计算该母线在该故障下的暂态电压安全性指标的电压恢复速度分量TVSI1:
其中,v(t)是暂态过程中的母线电压轨迹;
(2-1-2)计算该母线在该故障下的暂态电压安全性指标的低电压持续时间分量TVSI2:
其中,Tend是暂态过程结束的时刻,STspan是v(t)中电压持续低于0.75p.u.的时间段的集合,Tspan,max是STspan中时间跨度最长的时间段;
(2-1-3)计算该母线在该故障下的暂态电压安全性指标的电压振荡分量TVSI3:
(2-1-4)计算该母线在该故障下的电网暂态电压安全性指标TVSI:
TVSI=max(TVSI1,TVSI2,TVSI3)
(2-2)重复步骤(2-1),得到每条母线在不同故障下分别对应的TVSI;
(2-3)计算最终的暂态电压安全性指标,表达式如下:
其中,TVSIk,j是故障j下母线k的TVSI指标,l、m分别为母线数量和故障数量;
(3)仿真分析离线数据库中电网运行方式的暂态电压安全性,把每种电网运行方式对应的TVSI结果补充到离线数据库中,得到更新后的离线数据库;具体步骤如下:
(3-1)仿真分析离线数据库中每种电网运行方式的暂态电压安全性,记录母线的电压轨迹数据;
母线的电压轨迹的具体存储格式如下:
Vti={iOper,jFlt,kBus,v(t)}
VtDBi={Vt1,...,Vtp}
其中,Vti是存储的第i个母线电压轨迹,iOper是电网运行方式的编号,jFlt是故障的编号,kBus是电网中母线的编号,v(t)是母线的电压轨迹;VtDBi是存储Vt的集合,下标p是VtDB中的电压轨迹数量;
(3-2)计算电压轨迹的TVSI指标,把结果补充到离线数据库中,得到更新后的离线数据库:
定义TVSIi为电网运行方式i对应的TVSI指标集合,表达式如下:
TVSIi={TVSIi j,k|j∈SFlt,k∈ΩGrid}
其中,TVSIi j,k是母线k在电网运行方式i和故障j下对应的电压轨迹的TVSI指标,SFlt是电网故障的集合,ΩGrid是电网母线的集合;
把每种电网运行方式对应的TVSI结果补充到离线数据库中,更新后的离线数据库的存储格式具体如下:
DB'={Case'1,...,Case'n}
Case'i={i,CSVi,TVSIi}
其中,Case'i是更新后的离线数据库中的第i个样本,CSVi是该样本对应的控制状态量,DB'是更新后的离线数据库;
(4)定义电网运行方式间的距离,从更新后的离线数据库中选取当前电网运行方式的临近电网运行方式,用当前电网运行方式的临近电网运行方式拟合出当前电网运行方式的暂态电压安全性指标,得到当前电网运行方式的暂态电压安全性评估结果;具体步骤如下:
(4-1)定义d(i0,i1)为电网运行方式i0与电网运行方式i1之间的距离,表达式如下:
ωk,V≥0
ωk,θ≥0
(4-2)定义i0为当前电网运行方式,定义Sd为i0与更新后的离线数据库中各电网运行方式的距离的集合,表达式如下:
Sd={d(i0,i)|Case'i∈DB'}
定义Sd|0.001大于Sd中1/1000的元素,但小于Sd中999/1000的元素;定义Si0为与i0的临近电网运行方式集合,表达式如下:
(4-3)根据步骤(4-2)的结果,评估当前电网运行方式的暂态电压安全性:
若则当前电网运行方式的暂态电压安全性评估结果为稳定;若 则当前电网运行方式的暂态电压安全性评估结果为大概率稳定;若 则当前电网运行方式的暂态电压安全性评估结果为不确定稳定;若 则当前电网运行方式的暂态电压安全性评估结果为大概率失稳;若 则当前电网运行方式的暂态电压安全性评估结果为失稳。
本发明的特点及有益效果在于:
本发明通过设定合理的电网状态量的限值,以及合理的操作步骤,快速的生成了潮流结果符合实际需求的电网运行方式。通过设定合理的存储格式,减小了离线数据库占用的存储空间。
本发明定义一种基于母线电压轨迹的评估电网暂态电压安全性的指标。从三个方面综合考虑了电网的暂态电压安全性,能更综合的反映出电网的暂态电压安全性问题,使得接下来的评估结果更准确。
本发明利用数值仿真软件仿真分析离线数据库中电网运行方式的暂态电压安全性,把结果补充到离线数据库中。通过设定合理的存储格式,减小了离线数据库占用的存储空间。
本发明定义电网运行方式间的距离,基于离线数据库,用当前电网运行方式的临近电网运行方式来拟合出当前电网运行方式的暂态电压安全性。定义的电网运行方式间距离综合考虑了电网的控制状态量,能确保相互距离近的电网运行方式的暂态电压安全性也相似。通过设定合理的加权算法,确保了拟合结果的可信度。
具体实施方式
本发明提出一种数据驱动的在线快速评估电网暂态电压安全性的方法,下面结合具体实施例进一步详细说明如下。
本发明提出一种数据驱动的在线快速评估电网暂态电压安全性的方法,包括以下步骤:
(1)设定电网运行状态量的限值,包括:电网控制状态量的限值和电网潮流结果状态量的限值;随机生成大量电网运行方式,将符合电网运行状态量限值的电网运行方式及其对应的电网的控制状态量集合存储到离线数据库中;具体步骤如下:
(1-1)设定电网控制状态量的限值:
定义ΩGrid为电网中母线的集合。将电网中的母线分成三类,这三类母线的集合分别定义为ΩPV、ΩPQ、ΩVθ。ΩPV中的母线记为PV节点,这类母线的有功功率注入和电压幅值给定,无功功率注入和电压相角待求。ΩPQ中的母线记为PQ节点,这类母线的有功功率注入和无功功率注入给定,电压幅值和电压相角待求。ΩVθ中的母线记为Vθ节点,这类母线的电压幅值和电压相角给定,有功功率注入和无功功率注入待求。
定义CSV为电网控制状态量集合,表达式如下:
SV sp={Vi sp|i∈ΩPV∪ΩVθ}
Sθ sp={θi sp|i∈ΩVθ}
Pg,i sp,min≤Pg,i sp≤Pg,i sp,max,i∈ΩPV
Vi sp,min≤Vi sp≤Vi sp,max,i∈ΩPV∪ΩVθ
Pl,i sp,min≤Pl,i sp≤Pl,i sp,max,i∈ΩGrid
Ql,i sp,min≤Ql,i sp≤Ql,i sp,max,i∈ΩGrid
其中,为ΩPV中母线的发电机有功出力设定值的集合,SV sp为ΩPV和ΩVθ中母线的电压幅值设定值的集合,为电网中母线的有功负荷量设定值的集合,为电网中母线的无功负荷量设定值的集合,Sθ sp为ΩVθ中母线的电压相角设定值的集合。Pg,i sp为母线i的发电机有功出力设定值,Vi sp为母线i的电压幅值设定值,Pl,i sp为母线i的有功负荷量设定值,Ql,i sp为母线i的无功负荷量设定值,θi sp为母线i的电压相角设定值。Pg,i sp,min是母线i的发电机无功出力设定值下限,Pg,i sp,max是母线i的发电机无功出力设定值上限,Vi sp,min是母线i的电压幅值设定值下限,Vi sp,max是母线i的电压幅值设定值上限,Pl,i sp,min是母线i的有功负荷量设定值下限,Pl,i sp,max是母线i的有功负荷量设定值上限,Ql,i sp,min是母线i的无功负荷量设定值下限,Ql,i sp,max是母线i的无功负荷量设定值上限,是电网中母线电压相角设定值的最大值,是电网中母线电压相角设定值的最小值。
(1-2)设定电网潮流结果状态量的限值;
定义RSV为电网的潮流结果状态量集合,表达式如下:
SV Rlt={Vi Rlt|i∈ΩGrid}
Sθ Rlt={θi Rlt|i∈ΩGrid}
Pg,i Rlt,min≤Pg,i Rlt≤Pg,i Rlt,max,i∈ΩPV∪ΩVθ
Qg,i Rlt,min≤Qg,i Rlt≤Qg,i Rlt,max,i∈ΩPV∪ΩVθ
Pl,i Rlt,min≤Pl,i RltRlt≤Pl,i Rlt,max,i∈ΩGrid
Ql,i Rlt,min≤Ql,i Rlt≤Ql,i Rlt,max,i∈ΩGrid
Vi Rlt,min≤Vi Rlt≤Vi Rlt,max,i∈ΩGrid
其中,为ΩPV和ΩVθ中母线的发电机有功出力潮流结果的集合,为ΩPV和ΩVθ中母线的发电机无功出力潮流结果的集合,为电网中母线的有功负荷量潮流结果的集合,为电网中母线的无功负荷量潮流结果的集合,SV Rlt为电网中母线的电压幅值潮流结果的集合,Sθ Rlt为电网中母线的电压相角潮流结果的集合。Pg,i Rlt为母线i的发电机有功出力潮流结果,Qg,i Rlt为母线i的发电机无功出力潮流结果,Pl,i Rlt为母线i的有功负荷量潮流结果,Ql,i Rlt为母线i的发电机无功负荷量潮流结果,Vi Rlt为母线i的电压幅值潮流结果,θi Rlt为母线i的电压相角潮流结果。Pg,i Rlt,min是母线i的发电机有功出力潮流结果下限,Pg,i Rlt,max是母线i的发电机有功出力潮流结果上限,Qg,i Rlt,min是母线i的发电机无功出力潮流结果下限,Qg,i Rlt,max是母线i的发电机无功出力潮流结果上限,Pl,i Rlt,min是母线i的有功负荷量潮流结果下限,Pl,i Rlt,max是母线i的有功负荷量潮流结果上限,Ql,i Rlt,min是母线i的无功负荷量潮流结果下限,Ql,i Rlt,max是母线i的无功负荷量潮流结果上限,Vi Rlt,min是母线i的电压幅值潮流结果下限,Vi Rlt,max是母线i的电压幅值潮流结果上限,是电网中母线电压相角潮流结果的最大值,是电网中母线电压相角潮流结果的最小值。
(1-3)随机生成若干电网运行方式,筛选出潮流结果状态量和控制状态量均符合对应限值的电网运行方式,并将筛选后的电网运行方式对应的控制状态量结果存储到离线数据库中。
电网的电网运行方式由控制状态量唯一的确定。具体步骤如下。
(1-3-1)随机生成电网的控制状态量,构成一个CSV集合。
(1-3-2)对该CSV集合中的所有元素进行判定:
如果CSV中的元素均满足对应限值要求,则进入步骤(1-3-3);否则,放弃该CSV集合,重新返回步骤(1-3-1)。
(1-3-3)用数值仿真软件计算电网潮流,读取潮流结果状态量,构成RSV集合。
(1-3-4)对RSV集合中的所有元素进行判定:
如果RSV中的所有元素均满足对应限值要求,则进入步骤(1-3-5);否则,重新返回步骤(1-3-3)。
(1-3-5)将CSV集合作为第i个样本存储到离线数据库。具体的存储格式如下。
DB={Case1,...,Casen}
Casei={i,CSVi}
其中,Casei是离线数据库中的第i个样本,i是该样本的编号,CSVi是该样本对应的控制状态量集合。DB是离线数据库,下标n是离线数据库中的样本数量。
如果i<10000,则令i=i+1,然后重新返回步骤(3-1-1);否则,结束计算。
需要注意的是,常见的电网数据仿真软件都能完成本步骤的计算。例如电力***分析综合程序(PSASP),软件的使用方法见说明书。
(2)计算基于母线电压轨迹的评估电网暂态电压安全性指标,以评估多个故障下电网中母线的暂态电压安全性;具体步骤如下:(2-1)计算任一母线在任一故障下的暂态电压安全性指标,具体步骤如下:
(2-1-1)计算该母线在该故障下的暂态电压安全性指标的电压恢复速度分量;
定义TVSI1为暂态电压安全性指标的电压恢复速度分量,表达式如下:
其中,v(t)是暂态过程中的母线电压轨迹。
(2-1-2)计算该母线在该故障下的暂态电压安全性指标的低电压持续时间分量;
定义TVSI2为暂态电压安全性指标的低电压持续时间分量,表达式如下:
其中,Tend是暂态过程结束的时刻,STspan是v(t)中电压持续低于0.75p.u.的时间段的集合,Tspan,max是STspan中时间跨度最长的时间段。
(2-3)计算该母线在该故障下的暂态电压安全性指标的电压振荡分量;
定义TVSI3为暂态电压安全性指标的电压振荡分量,表达式如下:
(2-1-4)计算该母线在该故障下的电网暂态电压安全性指标;
定义TVSI为电网暂态电压安全性指标,表达式如下:
TVSI=max(TVSI1,TVSI2,TVSI3)
需要注意的是,这里的TVSI对应于一个母线在一个故障下的电压轨迹,如果需要考虑的母线或故障数量多于1个,则需要对TVSI做进一步处理。
(2-2)重复步骤(2-1),得到每条母线在不同故障下分别对应的TVSI;
(2-3)计算最终的暂态电压安全性指标,表达式如下:
其中,TVSIk,j是故障j下母线k的TVSI指标,l、m分别是考虑的母线数量和故障数量,TVSI是考虑了多个母线和故障情况下的暂态电压安全性指标。
(3)利用数值仿真软件仿真分析离线数据库中电网运行方式的暂态电压安全性,把每种电网运行方式对应的TVSI结果补充到离线数据库中,得到更新后的离线数据库;具体步骤如下:
(3-1)仿真分析离线数据库中每种电网运行方式的暂态电压安全性,记录母线的电压轨迹数据;
常见的电网数据仿真软件都能完成本步骤的计算。例如电力***分析综合程序(PSASP),软件的使用方法见说明书。
母线的电压轨迹的具体存储格式如下:
Vti={iOper,jFlt,kBus,v(t)}
VtDBi={Vt1,...,Vtp}
其中,Vti是存储的第i个母线电压轨迹,iOper是电网运行方式的编号,jFlt是故障的编号,kBus是电网中母线的编号,v(t)是母线的电压轨迹。VtDBi是存储Vt的集合,下标p是VtDB中的电压轨迹数量。
(3-2)计算电压轨迹的TVSI指标,把结果补充到离线数据库中,得到更新后的离线数据库:
定义TVSIi为电网运行方式i对应的TVSI指标集合,表达式如下:
TVSIi={TVSIi j,k|j∈SFlt,k∈ΩGrid}
其中,TVSIi j,k是母线k在电网运行方式i和故障j下对应的电压轨迹的TVSI指标。SFlt是电网故障的集合,ΩGrid是电网母线的集合。
把每种电网运行方式对应的TVSI结果补充到离线数据库中,更新后的离线数据库的存储格式具体如下:
DB'={Case'1,...,Case'n}
Case'i={i,CSVi,TVSIi}
其中,Case'i是更新后的离线数据库中的第i个样本,i是该样本的编号,CSVi是该样本对应的控制状态量。DB'是更新后的离线数据库,下标n是离线数据库中的样本数量。
(4)定义电网运行方式间的距离,基于更新后的离线数据库,从更新后的离线数据库中选取当前电网运行方式的临近电网运行方式,用当前电网运行方式的临近电网运行方式拟合当前电网运行方式的暂态电压安全性指标,得到当前电网运行方式的暂态电压安全性评估结果;具体步骤如下:
(4-1)定义电网运行方式间的距离:
定义d(i0,i1)为电网运行方式i0与电网运行方式i1之间的距离,表达式如下:
ωk,V≥0
ωk,θ≥0
其中,是母线k的发电机有功出力的权重,ωk,V是母线k的电压幅值的权重,是母线k的有功负荷量的权重,是母线k的无功负荷量的权重,ωk,θ是母线k的电压相角的权重。默认情况下,各权重系数都设定为1。如果更关注某一类电网状态变量(例如发电机有功出力),可以增加该类电网状态变量的权重系数。下标0对应电网运行方式i0,下标1对应电网运行方式i1。
(4-2)定义当前电网运行方式的临近电网运行方式集合,并定义拟合公式:
定义i0为当前电网运行方式,定义Sd为i0与更新后的离线数据库中各电网运行方式的距离的集合,表达式如下:
Sd={d(i0,i)|Case'i∈DB'}
其中,i为DB'中的一个电网运行方式,Case'i为DB'中的一个样本。
其中,i为DB'中的一个电网运行方式,Case'i为DB'中的一个样本。
(4-3)评估当前电网运行方式的暂态电压安全性:
定义评估当前电网运行方式的暂态电压安全性的判据,具体表述如表所示:
表1当前电网运行方式的暂态电压安全性评估表
Claims (1)
1.一种数据驱动的在线快速评估电网暂态电压安全性的方法,其特征在于,该方法包括以下步骤:
(1)设定电网运行状态量的限值,包括:电网控制状态量的限值和电网潮流结果状态量的限值;随机生成电网运行方式,分别计算每种电网运行方式对应的潮流结果,将符合电网运行状态量限值的电网运行方式及其对应的电网的控制状态量集合存储到离线数据库中;具体步骤如下:
(1-1)设定电网控制状态量的限值:
定义ΩGrid为电网中母线的集合,将ΩGrid分为三个子集合分为记为ΩPV、ΩPQ、ΩVθ;其中,ΩPV中的母线为PV节点,ΩPQ中的母线为PQ节点,ΩVθ中的母线为Vθ节点;
定义CSV为电网控制状态量集合,表达式如下:
SV sp={Vi sp|i∈ΩPV∪ΩVθ}
Sθ sp={θi sp|i∈ΩVθ}
Pg,i sp,min≤Pg,i sp≤Pg,i sp,max,i∈ΩPV
Vi sp,min≤Vi sp≤Vi sp,max,i∈ΩPV∪ΩVθ
Pl,i sp,min≤Pl,i sp≤Pl,i sp,max,i∈ΩGrid
Ql,i sp,min≤Ql,i sp≤Ql,i sp,max,i∈ΩGrid
其中,为ΩPV中母线的发电机有功出力设定值的集合,SV sp为ΩPV和ΩVθ中母线的电压幅值设定值的集合,为电网中母线的有功负荷量设定值的集合,为电网中母线的无功负荷量设定值的集合,Sθ sp为ΩVθ中母线的电压相角设定值的集合;Pg,i sp为母线i的发电机有功出力设定值,Vi sp为母线i的电压幅值设定值,Pl,i sp为母线i的有功负荷量设定值,Ql,i sp为母线i的无功负荷量设定值,θi sp为母线i的电压相角设定值;Pg,i sp,min是母线i的发电机无功出力设定值下限,Pg,i sp,max是母线i的发电机无功出力设定值上限,Vi sp,min是母线i的电压幅值设定值下限,Vi sp,max是母线i的电压幅值设定值上限,Pl,i sp,min是母线i的有功负荷量设定值下限,Pl,i sp,max是母线i的有功负荷量设定值上限,Ql,i sp,min是母线i的无功负荷量设定值下限,Ql,i sp,max是母线i的无功负荷量设定值上限,是电网中母线电压相角设定值的最大值,是电网中母线电压相角设定值的最小值;
(1-2)设定电网潮流结果状态量的限值;
定义RSV为电网的潮流结果状态量集合,表达式如下:
SV Rlt={Vi Rlt|i∈ΩGrid}
Sθ Rlt={θi Rlt|i∈ΩGrid}
Pg,i Rlt,min≤Pg,i Rlt≤Pg,i Rlt,max,i∈ΩPV∪ΩVθ
Qg,i Rlt,min≤Qg,i Rlt≤Qg,i Rlt,max,i∈ΩPV∪ΩVθ
Pl,i Rlt,min≤Pl,i RltRlt≤Pl,i Rlt,max,i∈ΩGrid
Ql,i Rlt,min≤Ql,i Rlt≤Ql,i Rlt,max,i∈ΩGrid
Vi Rlt,min≤Vi Rlt≤Vi Rlt,max,i∈ΩGrid
其中,为ΩPV和ΩVθ中母线的发电机有功出力潮流结果的集合,为ΩPV和ΩVθ中母线的发电机无功出力潮流结果的集合,为电网中母线的有功负荷量潮流结果的集合,为电网中母线的无功负荷量潮流结果的集合,SV Rlt为电网中母线的电压幅值潮流结果的集合,Sθ Rlt为电网中母线的电压相角潮流结果的集合;Pg,i Rlt为母线i的发电机有功出力潮流结果,Qg,i Rlt为母线i的发电机无功出力潮流结果,Pl,i Rlt为母线i的有功负荷量潮流结果,Ql,i Rlt为母线i的发电机无功负荷量潮流结果,Vi Rlt为母线i的电压幅值潮流结果,θi Rlt为母线i的电压相角潮流结果;Pg,i Rlt,min是母线i的发电机有功出力潮流结果下限,Pg,i Rlt,max是母线i的发电机有功出力潮流结果上限,Qg,i Rlt,min是母线i的发电机无功出力潮流结果下限,Qg,i Rlt,max是母线i的发电机无功出力潮流结果上限,Pl,i Rlt,min是母线i的有功负荷量潮流结果下限,Pl,i Rlt,max是母线i的有功负荷量潮流结果上限,Ql,i Rlt,min是母线i的无功负荷量潮流结果下限,Ql,i Rlt,max是母线i的无功负荷量潮流结果上限,Vi Rlt,min是母线i的电压幅值潮流结果下限,Vi Rlt,max是母线i的电压幅值潮流结果上限,是电网中母线电压相角潮流结果的最大值,是电网中母线电压相角潮流结果的最小值;
(1-3)随机生成若干电网运行方式,筛选出潮流结果状态量和控制状态量均符合对应限值的电网运行方式,并将筛选后的电网运行方式对应的控制状态量结果作为一个样本存储到离线数据库中;
离线数据库的存储格式如下:
DB={Case1,...,Casen}
Casei={i,CSVi}
其中,Casei是离线数据库中的第i个样本,CSVi是该样本对应的控制状态量集合,DB是离线数据库,下标n是离线数据库中的样本数量;
(2)计算基于母线电压轨迹的评估电网暂态电压安全性指标,具体步骤如下:
(2-1)计算任一母线在任一故障下的暂态电压安全性指标,具体步骤如下:
(2-1-1)计算该母线在该故障下的暂态电压安全性指标的电压恢复速度分量TVSI1:
其中,v(t)是暂态过程中的母线电压轨迹;
(2-1-2)计算该母线在该故障下的暂态电压安全性指标的低电压持续时间分量TVSI2:
其中,Tend是暂态过程结束的时刻,STspan是v(t)中电压持续低于0.75p.u.的时间段的集合,Tspan,max是STspan中时间跨度最长的时间段;
(2-1-3)计算该母线在该故障下的暂态电压安全性指标的电压振荡分量TVSI3:
(2-1-4)计算该母线在该故障下的电网暂态电压安全性指标TVSI:
TVSI=max(TVSI1,TVSI2,TVSI3)
(2-2)重复步骤(2-1),得到每条母线在不同故障下分别对应的TVSI;
(2-3)计算最终的暂态电压安全性指标,表达式如下:
其中,TVSIk,j是故障j下母线k的TVSI指标,l、m分别为母线数量和故障数量;
(3)仿真分析离线数据库中电网运行方式的暂态电压安全性,把每种电网运行方式对应的TVSI结果补充到离线数据库中,得到更新后的离线数据库;具体步骤如下:
(3-1)仿真分析离线数据库中每种电网运行方式的暂态电压安全性,记录母线的电压轨迹数据;
母线的电压轨迹的具体存储格式如下:
VtDBi={Vt1,...,Vtp}
其中,Vti是存储的第i个母线电压轨迹,iOper是电网运行方式的编号,jFlt是故障的编号,kBus是电网中母线的编号,v(t)是母线的电压轨迹;VtDBi是存储Vt的集合,下标p是VtDB中的电压轨迹数量;
(3-2)计算电压轨迹的TVSI指标,把结果补充到离线数据库中,得到更新后的离线数据库:
定义TVSIi为电网运行方式i对应的TVSI指标集合,表达式如下:
TVSIi={TVSIi j,k|j∈SFlt,k∈ΩGrid}
其中,TVSIi j,k是母线k在电网运行方式i和故障j下对应的电压轨迹的TVSI指标,SFlt是电网故障的集合,ΩGrid是电网母线的集合;
把每种电网运行方式对应的TVSI结果补充到离线数据库中,更新后的离线数据库的存储格式具体如下:
DB'={Case'1,...,Case'n}
Case'i={i,CSVi,TVSIi}
其中,Case'i是更新后的离线数据库中的第i个样本,CSVi是该样本对应的控制状态量,DB'是更新后的离线数据库;
(4)定义电网运行方式间的距离,从更新后的离线数据库中选取当前电网运行方式的临近电网运行方式,用当前电网运行方式的临近电网运行方式拟合出当前电网运行方式的暂态电压安全性指标,得到当前电网运行方式的暂态电压安全性评估结果;具体步骤如下:
(4-1)定义d(i0,i1)为电网运行方式i0与电网运行方式i1之间的距离,表达式如下:
ωk,V≥0
ωk,θ≥0
(4-2)定义i0为当前电网运行方式,定义Sd为i0与更新后的离线数据库中各电网运行方式的距离的集合,表达式如下:
Sd={d(i0,i)|Case'i∈DB'}
(4-3)根据步骤(4-2)的结果,评估当前电网运行方式的暂态电压安全性:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911219155.3A CN110957723B (zh) | 2019-12-03 | 2019-12-03 | 一种数据驱动的在线快速评估电网暂态电压安全性的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911219155.3A CN110957723B (zh) | 2019-12-03 | 2019-12-03 | 一种数据驱动的在线快速评估电网暂态电压安全性的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110957723A true CN110957723A (zh) | 2020-04-03 |
CN110957723B CN110957723B (zh) | 2021-05-07 |
Family
ID=69979482
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911219155.3A Active CN110957723B (zh) | 2019-12-03 | 2019-12-03 | 一种数据驱动的在线快速评估电网暂态电压安全性的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110957723B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113852123A (zh) * | 2021-09-09 | 2021-12-28 | 国网江苏省电力有限公司 | 一种电力***自动电压控制方法、装置、电子设备和存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103701140A (zh) * | 2014-01-06 | 2014-04-02 | 国家电网公司 | 提高交直流电网暂态电压稳定的动态无功备用优化方法 |
CN106229976A (zh) * | 2016-08-31 | 2016-12-14 | 山东大学 | 基于数据驱动的暂态功角稳定态势预估方法 |
CN109038636A (zh) * | 2018-08-06 | 2018-12-18 | 国家电网公司华东分部 | 数据驱动的直流受端电网动态无功储备需求评估方法 |
CN109033702A (zh) * | 2018-08-23 | 2018-12-18 | 国网内蒙古东部电力有限公司电力科学研究院 | 一种基于卷积神经网络cnn的电力***暂态电压稳定评估方法 |
CN109766853A (zh) * | 2019-01-16 | 2019-05-17 | 华北电力大学 | 基于lstm的电压暂降扰动分类方法 |
-
2019
- 2019-12-03 CN CN201911219155.3A patent/CN110957723B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103701140A (zh) * | 2014-01-06 | 2014-04-02 | 国家电网公司 | 提高交直流电网暂态电压稳定的动态无功备用优化方法 |
CN106229976A (zh) * | 2016-08-31 | 2016-12-14 | 山东大学 | 基于数据驱动的暂态功角稳定态势预估方法 |
CN109038636A (zh) * | 2018-08-06 | 2018-12-18 | 国家电网公司华东分部 | 数据驱动的直流受端电网动态无功储备需求评估方法 |
CN109033702A (zh) * | 2018-08-23 | 2018-12-18 | 国网内蒙古东部电力有限公司电力科学研究院 | 一种基于卷积神经网络cnn的电力***暂态电压稳定评估方法 |
CN109766853A (zh) * | 2019-01-16 | 2019-05-17 | 华北电力大学 | 基于lstm的电压暂降扰动分类方法 |
Non-Patent Citations (3)
Title |
---|
WENLU ZHAO 等: "Practical short-term voltage stability index based on voltage curves: definition,verification and case studies", 《IET GENERATION, TRANSMISSION & DISTRIBUTION》 * |
YAN XU 等: "Assessing Short-Term Voltage Stability of Electric Power Systems by a Hierarchical Intelligent System", 《IEEE TRANSACTIONS ON NEURAL NETWORKS AND LEARNING SYSTEMS》 * |
刘俊磊 等: ""基于支持向量回归的直流受端电网动态无功需求在线评估"", 《电力***保护与控制》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113852123A (zh) * | 2021-09-09 | 2021-12-28 | 国网江苏省电力有限公司 | 一种电力***自动电压控制方法、装置、电子设备和存储介质 |
CN113852123B (zh) * | 2021-09-09 | 2024-02-06 | 国网江苏省电力有限公司 | 一种电力***自动电压控制方法、装置、电子设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN110957723B (zh) | 2021-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Probabilistic approach for power capacity specification of wind energy storage systems | |
CN102541621B (zh) | 风光储联合发电***的仿真等值方法 | |
CN106058937B (zh) | 一种含直驱永磁风电场的配电网广义负荷建模方法 | |
CN107666155A (zh) | 基于Markov模型的多能互补***随机稳定性分析方法 | |
CN115102200B (zh) | 一种储能***模型辨识方法和*** | |
CN105576654B (zh) | 直驱风电场等值方法和*** | |
CN105762816B (zh) | 一种抑制次同步相互作用的方法 | |
CN105244890A (zh) | 新能源并网的无功优化方法 | |
CN114123344A (zh) | 基于自适应递推最小二乘的电力***惯量评估方法及装置 | |
CN115392133A (zh) | 采用高斯混合模型的风电场最优聚类等值方法 | |
CN105356508B (zh) | 基于psd-bpa的电网风电接入评价***及其方法 | |
Wu et al. | Equivalent modeling method for regional decentralized photovoltaic clusters based on cluster analysis | |
Eltamaly et al. | Analyzing of wind distributed generation configuration in active distribution network | |
CN111262238A (zh) | 基于机器学习的含iidg配电网短路电流预测方法 | |
CN110957723B (zh) | 一种数据驱动的在线快速评估电网暂态电压安全性的方法 | |
Bagheri et al. | Extremely short time modeling of wind power variations | |
CN117335426A (zh) | 一种适用于短期电压稳定性评估的关键特征识别方法 | |
CN106786766A (zh) | 一种基于p2g技术的提高风电最大并网容量的方法 | |
Zhang et al. | Optimization method of under frequency load shedding schemes for systems with high permeability new energy | |
CN112993984B (zh) | 一种分布式光伏电站在有电气连接情况下的大扰动聚合方法 | |
Zheng et al. | Feature distance based online cluster modeling of LVRT controlled PV power plants | |
CN111817359B (zh) | 一种基于等效电路的微电网可求解边界分析的方法 | |
Efe et al. | A novel approach to power flow analysis for grid connected micro grid | |
Kumar et al. | Performance analysis of FLC controlled PV-wind hybrid power system for dc load with real-time data in Matlab, Simulink | |
Li et al. | Wind farms reactive power optimization using genetic/tabu hybrid algorithm |
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 |