CN114301055A - 基于宽频量测的电力系统间谐波潮流获取方法与系统 - Google Patents
基于宽频量测的电力系统间谐波潮流获取方法与系统 Download PDFInfo
- Publication number
- CN114301055A CN114301055A CN202111358435.XA CN202111358435A CN114301055A CN 114301055 A CN114301055 A CN 114301055A CN 202111358435 A CN202111358435 A CN 202111358435A CN 114301055 A CN114301055 A CN 114301055A
- Authority
- CN
- China
- Prior art keywords
- inter
- harmonic
- oscillation
- power
- current
- 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
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000005259 measurement Methods 0.000 title claims abstract description 25
- 230000010355 oscillation Effects 0.000 claims abstract description 103
- 238000004364 calculation method Methods 0.000 claims abstract description 54
- 230000001360 synchronised effect Effects 0.000 claims abstract description 30
- 239000011159 matrix material Substances 0.000 claims abstract description 19
- 238000004458 analytical method Methods 0.000 claims abstract description 12
- 238000012545 processing Methods 0.000 claims abstract description 9
- 238000012937 correction Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 8
- 238000004088 simulation Methods 0.000 claims description 8
- 238000004804 winding Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 4
- 230000001939 inductive effect Effects 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 9
- 238000005206 flow analysis Methods 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 3
- 230000010485 coping Effects 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005191 phase separation Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Images
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/40—Arrangements for reducing harmonics
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供一种基于宽频量测的电力系统间谐波潮流获取方法与系统,包括:获得各个次/超同步振荡源向电力系统注入的振荡电流数据;通过prony分析对振荡电流数据进行解耦处理,获得电力系统中主要稳定的次/超同步振荡参数;根据振荡频率分量与电力系统的网络结构参数,构建某一振荡频率分量下的系统节点导纳矩阵;量测等效得到的间谐波振荡源数据和给定的间谐波振荡源数据;将接地电压和电网系统内部电流设置为零,设置边界节点电流或电压;利用间谐波潮流计算公式进行计算,获得次/超同步振荡的间谐波在电网系统中传播路径与分布数据。通过发明可快速评估可能引发的传统汽轮发电机组谐振风险,提高电网系统应对次/超同步振荡的能力。
Description
技术领域
本发明涉及智能电网技术领域,尤其是新能源发电并网接入系统中的电力系统次/超同步振荡的间谐波潮流估计技术领域,具体而言涉及一种基于宽频量测的电力系统间谐波潮流获取方法与系统。
背景技术
近些年来,随着大量新能源机组通过电力电子变流器接入电网,电力系统逐渐呈现电力电子化,电网系统呈现出高比例的新能源和高比例电力电子装置的双高特性。
在众多电力电子设备接入电网的情境下,电网的拓扑结构愈发复杂化,出现了许多新型的次/超同步振荡问题,使系统的稳定运行面临巨大的挑战。例如2015年以来,在新疆哈密地区频繁出现风电机群参与的次同步振荡,频率在20~40Hz内变化,次同步振荡功率穿越35/110/220/500/750kV多级电网,甚至激发网中汽轮机组轴系扭振,造成300km外的某电厂机组全跳和特高压直流功率骤降的事故。
为了消除个防治次/超同步振荡问题,对次/超同步振荡问题产生的条件和环境数据进行分析,分析次/超同步振荡的传播路径与分布是必要的,需要进行间谐波潮流分析,对于大电网次/超同步振荡的预防控制至关重要。潮流是指电力系统中电压(各节点)、功率(有功、无功)(各支路)的稳态分布,潮流分析与计算过程是据给定的运行参数确定系统的运行状态,如计算网络中各节点的电压(幅值和相角)和各支路中的功率分布及损耗,通过在发电机母线上功率被注入网络,在变(配)电站上接入负荷;其间,功率在网络中流动。
对于正在运行的电力系统,通过潮流分析可以判断电网母线电压、支路电流和功率是否越限,如果有越限,就应采取措施,调整运行方式。对于正在规划的电力系统,通过潮流计算,可以为选择电网供电方案和电气设备提供依据。同时,通过潮流分析计算还可以为继电保护和自动装置整定计算、电力系统故障计算和稳定计算等提供原始数据。
现有技术中,对电力系统的潮流分析主要是针对基波的潮流计算以及非线性元件导致的谐波潮流计算,无法实现关于电网出现的次/超同步振荡现象的间谐波潮流分析。
现有技术文献:
专利文献1:CN109494747一种基于交替梯度算法的电网概率潮流计算方法
专利文献2:CN113452028低压配电网概率潮流计算方法、系统、终端和存储介质
发明内容
本发明旨在针对现有技术存在的缺陷或者不足,提出一种基于宽频量测的电力系统间谐波潮流获取方法,基于量测量获取各个节点的间谐波电流数据,通过prony辨识对量测量进行解耦,获取谐波潮流计算所需的主要次/超同步振荡频率与幅值,作为间谐波潮流计算的主要数据进行间谐波潮流分析计算,根据获得的间谐波电流/电压分布,可评估次/超同步振荡对大电网中汽轮发电机组的潜在影响。
为实现上述目的,本发明的第一方面提出一种基于宽频量测的电力系统间谐波潮流获取方法,包括:
步骤1、基于安装在电网系统的网络结构各处的测量设备,获得各个次/超同步振荡源向电力系统注入的振荡电流数据;
步骤2、通过prony分析,对振荡电流数据进行解耦处理,获得电力系统中主要稳定的次/超同步振荡电流频率与幅值;
步骤3、根据振荡频率分量与电力系统的网络结构参数,构建某一振荡频率分量下的系统节点导纳矩阵;
步骤4、判断网络结构参数是否完整已知,如果不完整则进入步骤5,若完整则跳转进入步骤6;
步骤5、对电网系统依据网络结构参数是否已知进行划分,量测已知区域与未知网络参数区域之间联络线电流或电压数据,通过prony分析后,将节点等效为间谐波振荡源;
步骤6、量测等效得到的间谐波振荡源数据和给定的间谐波振荡源数据;
步骤7、将接地电压和电网系统内部电流设置为零,设置边界节点电流或电压;
步骤8、利用间谐波潮流计算公式进行计算,获得次/超同步振荡的间谐波在电网系统中传播路径与分布数据。
其中,在所述步骤3中,根据振荡频率分量与电力系统的网络结构参数,构建某一振荡频率分量下的系统节点导纳矩阵,具体包括:
对于一个有n个节点的电力系统,其系统节点导纳矩阵如下所示:
式中,Yij为节点i与节点j之间的互导纳,Yii为节点i的自导纳;
将电力系统中任意一段线路的元件模型等效为2类,分别为:第一等效模型;第二等效模型;并分别建立自导纳与互导纳公式;
在上述第一等效模型、第二等效模型的基础上,对受频率影响的电力系统元件进行基于频率的校正,具体校正过程包括:
1)对于一条π型等值电路,校正公式如下:
式中,zl为线路阻抗;yl为线路对地导纳;r、x、g和b分别为基频下线路电阻、电抗、电导和电纳;f和fn分别为间谐波频率和基频频率;
即对于π型等值电路,适用上述第一等效模型,zij=zl、yi=yj=yl/2;
2)对于并联RLC恒阻抗负荷,校正公式如下:
式中:R、Xc和Xl是并联RLC恒阻抗负荷的等效电阻、电容和电感;Vn是负荷额定电压;Pn、Qc和Ql是负荷有功功率、容性无功功率和感性有功功率;
即对于并联RLC恒阻抗负荷,适用第二等效模型,yi=R//Xc//Xl;
3)对于恒功率负荷,保持视在功率S恒定不变为最强等量约束,其建模满足以下公式:
即对于并联RLC恒阻抗负荷,适用上述第二等效模型,yi=If/Uf;
4)对于双绕组变压器而言,其计算公式如下式:
即对于双绕组变压器,适用第一等效模型,zij=zb、yi=yj=0;
5)对于同步发电机产生的间谐波电流远小于基频频率的电流,在潮流计算过程中将同步发电机视为负荷处理:
Z=R+jX
式中:Z为同步发电机间谐波频率下等效阻抗,R和X为同步发电机间谐波频率下等效电阻与等效电抗;
其中,利用间谐波潮流计算公式进行计算,获得次/超同步振荡的间谐波在电网系统中传播路径与分布数据,包括:
间谐波潮流计算选择使用电流平衡方程,其具体公式如下:
In=YnVn
式中:In为节点流入或流出的某一频率间谐波电流;Yn为对应频率下的系统节点导纳矩阵;Vn为节点在对应频率的间谐波电压;
对于一个完整可知的系统,已知间谐波振荡源的流入电流,所有边界节点的间谐波电压为零,所有接地节点的间谐波电流为零,上述电流平衡方程转换成以下形式:
式中,IB、II和IG分别是系统边界节点、内部节点与接地节点的电流,UB、UI和UG分别是系统内部节点与边界节点的电压,YBB、YII、YGG、YBI、YIB、YGI和YIG分别是系统边界、内部、外部的自导纳,以及边界与内部、内部与边界、接地与内部、内部与接地之间的互导纳。
其中,接地节点电压UG和内部节点电流II等于零,设定边界电流IB已知,可由以下公式计算出边界电压UB:
即在可知系统中,仅通过对间谐波振荡源与可知系统连接线路的电流量测,计算得到间谐波振荡源在整个可知系统的传播与分布。
本发明以上方案提出的种基于量测量与prony辨识的电力系统间谐波潮流计算方法,通过测量设备获得各个次/超同步振荡源向电力网络注入的次/超同步电流数据,利用prony方法对电流数据进行解耦处理,构建间谐波潮流计算的节点导纳矩阵,进行间谐波潮流计算,进而根据次/超同步振荡在大电网中的传播与分布,可快速评估可能引发的传统汽轮发电机组谐振风险。
本发明目的的第二方面提供一种基于宽频量测的电力系统间谐波潮流获取的计算机系统,包括:
一个或多个处理器;
存储器,存储可被操作的指令,所述指令在通过所述一个或多个处理器执行时使得所述一个或多个处理器执行操作,所述操作包括执行前述获取方法的过程。
由以上本发明的技术方案,与现有技术相比,本发明的显著优点在于:
1、本发明实现了对于目前电力系统潮流计算中仅存在的针对工频的基波潮流计算和针对工频及其整数倍频的谐波潮流计算的补充,弥补了目前对于间谐波潮流计算的空白,可实现对未来电网双高的特性下多发的次/超同步振荡的解决,用于引发网中汽轮发电机组谐振风险评估,提高电网的预警和应对能力;
2、现有次/超同步振荡模型仿真的结构复杂、仿真耗时久,本发明提出的离线获取方法,利用间谐波潮流计算能够在只获取测量数据的条件下获得间谐波潮流在电网中分布,缩减分析和处理的时长。同时,未来电网系统中基于宽频量测技术,可以更好地能满足电网对次/超同步振荡实时传播控制的需求;
3、本发明提出的离线获取方法中,所建立的模型的电流平衡方程的间谐波潮流计算本质上是线性的,所以不存在收敛与精确性的问题,而不收敛往往是网络或状态参数设置不合理导致的,并且由于间谐波潮流计算方法是线性的,其不需要迭代循环计算,只需要单次计算便可得到整个系统的间谐波分布,获取和计算效率高,计算量小,而且耗时短,间谐波潮流分析速度得到了保证。
应当理解,前述构思以及在下面更加详细地描述的额外构思的所有组合只要在这样的构思不相互矛盾的情况下都可以被视为本公开的发明主题的一部分。另外,所要求保护的主题的所有组合都被视为本公开的发明主题的一部分。
结合附图从下面的描述中可以更加全面地理解本发明教导的前述和其他方面、实施例和特征。本发明的其他附加方面例如示例性实施方式的特征和/或有益效果将在下面的描述中显见,或通过根据本发明教导的具体实施方式的实践中得知。
附图说明
附图不意在按比例绘制。在附图中,在各个图中示出的每个相同或近似相同的组成部分可以用相同的标号表示。为了清晰起见,在每个图中,并非每个组成部分均被标记。现在,将通过例子并参考附图来描述本发明的各个方面的实施例,其中:
图1是说明根据本发明某些实施方式的一种基于量测量与prony辨识的电力系统间谐波潮流计算方法的流程示意图。
图2是说明根据本发明某些实施方式的一种基于量测量与prony辨识的电力系统间谐波潮流计算方法的电流起振与等幅振荡波形图。
图3(a)-3(d)是说明根据本发明某些实施方式的电力系统元件间谐波频率下四类不同的等效电路图,分别对应于π型等效电路、并联RLC恒阻抗负荷、变压器等效电路以及同步发电机等效电路。
图4是说明根据本发明某些实施方式的第一等效模型的示意图。
图5是说明根据本发明某些实施方式的第二等效模型的示意图。
图6是说明根据本发明某些实施方式的节点分类展示图。
图7是说明根据本发明某些实施方式的案例十机三十九节点网络等效阻抗图。
图8是说明根据本发明某些实施方式的间谐波潮流计算结果与时域仿真结果电压对比图。
图9是说明根据本发明某些实施方式的间谐波潮流计算结果与时域仿真结果电流对比图。
具体实施方式
为了更了解本发明的技术内容,特举具体实施例并配合所附图式说明如下。
在本公开中参照附图来描述本发明的各方面,附图中示出了许多说明的实施例。本公开的实施例不必定意在包括本发明的所有方面。应当理解,上面介绍的多种构思和实施例,以及下面更加详细地描述的那些构思和实施方式可以以很多方式中任意一种来实施,这是因为本发明所公开的构思和实施例并不限于任何实施方式。另外,本发明公开的一些方面可以单独使用,或者与本发明公开的其他方面的任何适当组合来使用。
根据本发明的实施例,一种基于量测量与prony辨识的电力系统间谐波潮流计算方法,包括:
步骤1、基于安装在网络各处的测量设备,获得各个次/超同步振荡源向电力网络注入的振荡电流数据;
步骤2、通过prony分析程序,对振荡电流数据进行解耦处理,获得系统中主要稳定的次/超同步振荡电流频率与幅值;
步骤3、根据振荡频率分量与电力网络结构参数,构建某一振荡频率分量下的系统节点导纳矩阵;
步骤4、判断网络结构参数是否完整已知,不完整进入步骤5,若完整则跳入步骤6;
步骤5、对系统依据网络结构是否已知进行划分,量测已知区域与未知区域之间联络线电流或电压数据,prony分析后,将节点等效为间谐波振荡源;
步骤6、测等效得到的间谐波振荡源数据和给定的间谐波振荡源数据;
步骤7、将接地电压和系统内部电流设置为零,设置边界节点电流或电压;
步骤8、利用间谐波潮流计算公式进行计算,得到间谐波在大电网中传播与分布,研究次/超同步振荡在大电网中的传播与分布,快速评估可能引发的传统汽轮发电机组谐振风险。
下面结合图1所示的流程图以及本发明的一些优选或者可选的例子,更加具体地描述本发明的某些实例的实施和/或效果。
【测量获取次/超同步振荡源注入电网电流】
结合附图1,在前述步骤1中,测量获取次/超同步振荡源向潮流计算的网络中流入/流出电网电流数据。这里的测量节点应当涉及潮流计算的所有边界节点,旨在测量潮流计算网络与外部网络的电流的流通交互。对于潮流计算网络内部节点的电流,根据基尔霍夫电流定律可知流入/流出电流为零。
振荡电流数据可以是电网实际量测数据或者仿真获得的振荡电流时序数据ik(t)。
ik(t)为第k个边界节点测得的电流在t时刻的值,k=1,2...n。
在本发明的实施例中,振荡电流数据为电网实际量测获得的振荡电流数据,或者通过仿真获得的振荡电流数据(即数组),具体表示为:
i(t)=[i1(t),i2(t),...in(t)]T
其中,in(t)为第k个边界节点测得的电流在t时刻的值,n为构建电力系统的网络结构的边界节点总数。
【利用prony对数据进行处理】
结合附图1和附图2,在前述步骤2中,利用前面步骤1中获得的振荡电流时序数据,分节点对其进行prony分析,对振荡电流数据进行解耦处理,获得系统中主要稳定的次/超同步振荡电流频率与幅值。
其中,将不同频率的振荡分量解耦开来,筛除迅速衰减和仍处于起振或发散的的振荡分量,本发明仅适用于稳态计算。
在prony分析中选择合适的时间窗,记录次/超同步振荡电流稳定振荡情况下的参数,即次/超同步振荡电流频率、幅值和相位。
得到某一振荡频率下prony处理后的各个节点的电流数据如下:
Is=[I1s αI1s I2s αI2s ... Ins αIns]T
I为对电流数据幅值与相位分开,组成一个2n*1的列向量。
【构建节点导纳矩阵】
结合附图1和附图3(a)-3(d),在前述步骤3中,根据振荡频率分量与电力网络结构参数,构建某一振荡频率分量下的系统节点导纳矩阵,也即利用前面步骤2中获得稳定振荡的次/超同步振荡电流分量与系统已知的网络结构参数构建某一振荡频率下的节点导纳矩阵。
对于一个有n个节点的电力系统,其系统节点导纳矩阵如下所示:
式中,Yij为节点i与节点j之间的互导纳,Yii为节点i的自导纳;
其中,i,j=1、2、3...n,i≠j。
将电力系统中任意一段线路的元件模型等效为2类,分别为:第一等效模型;第二等效模型;并分别建立自导纳与互导纳公式。
结合图4所示的第一等效模型的示意图。第一等效模型的自导纳与互导纳构建如下:
结合图5所示的第二等效模型的示意图,第二等效模型的自导纳与互导纳构建如下:
结合图示,yi为节点i由对地阻抗(一般为对地电容)计算得到的导纳。
在上述第一等效模型、第二等效模型的基础上,对受频率影响的电力系统元件进行基于频率的校正,具体校正过程包括:
1)对于一条π型等值电路,校正公式如下:
式中,zl为线路阻抗;yl为线路对地导纳;r、x、g和b分别为基频下线路电阻、电抗、电导和电纳;f和fn分别为间谐波频率和基频频率;
即对于π型等值电路,适用上述第一等效模型,zij=zl、yi=yj=yl/2;
2)对于并联RLC恒阻抗负荷,校正公式如下:
式中:R、Xc和Xl是并联RLC恒阻抗负荷的等效电阻、电容和电感;Vn是负荷额定电压;Pn、Qc和Ql是负荷有功功率、容性无功功率和感性有功功率;
即对于并联RLC恒阻抗负荷,适用第二等效模型,yi=R//Xc//Xl;
3)对于恒功率负荷,保持视在功率S恒定不变为最强等量约束,其建模满足以下公式:
式中:If、和分别是对应于频率分量为f、2fn-f和2fn+f的间谐波电流;Uf、Ufn分别为间谐波频率对应的电压与基波频率对应的电压;Pn、Qn分别为系统设定的恒定有功、恒定无功;θ是系统设定最初相角;
即对于并联RLC恒阻抗负荷,适用上述第二等效模型,yi=If/Uf;
4)对于双绕组变压器而言,其计算公式如下式:
即对于双绕组变压器,适用第一等效模型,zij=zb、yi=yj=0;
5)对于同步发电机产生的间谐波电流远小于基频频率的电流,在潮流计算过程中将同步发电机视为负荷处理:
Z=R+jX
式中:Z为同步发电机间谐波频率下等效阻抗,R和X为同步发电机间谐波频率下等效电阻与等效电抗;
【网络结构参数完整判断】
结合附图1,在前述步骤4中,对网络结构参数完整与否进行判定。
在本发明的实施例中,网络参数的完整判断取决于是否系统在间谐波频率分量下的等效导纳矩阵是否能够建立,构建的节点导纳矩阵是否有不可计算导纳Yij,使节点导纳矩阵不完整。
由于元件厂家参数的保密导致的参数未知、非线性元件的无法成功建模的情况均属于网络参数不完整的情况。
【量测等效电流/电压源】
结合附图1和附图6,在前述步骤5和步骤6中,对未知网络参数区域进行电流/电压源等效。
可以通过量测未知网络参数区域与已知区域联络线传输电流/电压参数,将未知网络参数区域等效为电流源或电压源,用于间谐波潮流计算。
【间谐波电流/电压初值设置】
结合附图1和附图7,在前述步骤7中,进行间谐波电流与电压进行赋初值。
首先,对潮流计算网络电压进行赋初值:
U=[e1,f1,e2,f2…en,fn]′
U为对电压数据实部与虚部分开,组成一个2n*1的列向量。
其次,对潮流计算网络电流进行赋初值:
I=[p1,q1,p2,q2…pn,qn]′
I为对电压数据实部与虚部分开,组成一个2n*1的列向量。
节点分类:边界节点,内部节点,接地节点。
其中接地节点电压初值为零,内部节点电流初值为零,边界节点取决于边界为等效电流源或电压源设置电压或者电流初值。
其他参数可根据发明的实施例,按照实际情况设置。
【间谐波潮流计算】
结合附图1、附图8和附图9,在前述步骤8中,进行间谐波潮流计算。
间谐波潮流计算不同以往基波潮流,选择使用电流平衡方程,表达如下:
In=YnVn
式中:In为节点流入或流出的某一频率间谐波电流;Yn为对应频率下的系统节点导纳矩阵;Vn为节点在对应频率的间谐波电压;
对于一个完整可知的系统,已知间谐波振荡源的流入电流,所有边界节点的间谐波电压为零,所有接地节点的间谐波电流为零,上述电流平衡方程转换成以下形式:
式中:IB、II和IG分别是系统边界节点、内部节点与接地节点的电流,UB、UI和UG分别是系统内部节点与边界节点的电压,YBB、YII、YGG、YBI、YIB、YGI和YIG分别是系统边界、内部、外部的自导纳,以及边界与内部、内部与边界、接地与内部、内部与接地之间的互导纳;
其中,接地节点电压UG和内部节点电流II等于零,设定边界电流IB已知,可由以下公式计算出边界电压UB:
即在可知系统中,仅通过对间谐波振荡源与可知系统连接线路的电流量测,计算得到间谐波振荡源在整个可知系统的传播与分布。
附图6和附图7为案例十机三十九节点情况下,间谐波潮流计算结果与时域仿真结果对比图。十机三十九节点系统的间谐波潮流计算结果与时域仿真结果两者相近大致准确。同样接地节点电压B40-58、系统内部节点电流B1-39(除B16外),纵然相位上呈现很大误差,但结合节点幅值准确且接近于零,间谐波潮流计算结果准确性仍可以保证。
通过数据可知十机三十九节点系统的间谐波潮流结果能与时域仿真结果保持较好的一致性,由此验证间谐波潮流方法的准确性。
通过本发明提出的基于量测量与prony辨识的电力系统间谐波潮流计算方法,可快速计算得到间谐波在大电网中传播与分布,分析获取次/超同步振荡在大电网中的传播与分布,快速评估可能引发的传统汽轮发电机组谐振风险,提高电网系统应对次/超同步振荡的能力。
根据本发明以上实施例的方法,在另外的实施例中,前述实施例方法还可以在计算机系统或者服务器中,以可被处理器调用和执行的指令集的形式存储在存储器中,使得这些指令集被处理器调用时执行前述电力系统间谐波潮流计算方法的过程。
本发明另一方面的实施例还公开一种基于宽频量测的电力系统间谐波潮流获取的计算机系统,包括:
一个或多个处理器;
存储器,存储可被操作的指令,所述指令在通过所述一个或多个处理器执行时使得所述一个或多个处理器执行操作,所述操作包括执行前述获取方法的过程。
虽然本发明已以较佳实施例揭露如上,然其并非用以限定本发明。本发明所属技术领域中具有通常知识者,在不脱离本发明的精神和范围内,当可作各种的更动与润饰。因此,本发明的保护范围当视权利要求书所界定者为准。
Claims (9)
1.一种基于宽频量测的电力系统间谐波潮流获取方法,其特征在于,包括以下步骤:
步骤1、基于安装在电网系统的网络结构各处的测量设备,获得各个次/超同步振荡源向电力系统注入的振荡电流数据;
步骤2、通过prony分析,对振荡电流数据进行解耦处理,获得电力系统中主要稳定的次/超同步振荡电流频率与幅值;
步骤3、根据振荡频率分量与电力系统的网络结构参数,构建某一振荡频率分量下的系统节点导纳矩阵;
步骤4、判断网络结构参数是否完整已知,如果不完整则进入步骤5,若完整则跳转进入步骤6;
步骤5、对电网系统依据网络结构参数是否已知进行划分,量测已知区域与未知网络参数区域之间联络线电流或电压数据,通过prony分析后,将节点等效为间谐波振荡源;
步骤6、量测等效得到的间谐波振荡源数据和给定的间谐波振荡源数据;
步骤7、将接地电压和电网系统内部电流设置为零,设置边界节点电流或电压;
步骤8、利用间谐波潮流计算公式进行计算,获得次/超同步振荡的间谐波在电网系统中传播路径与分布数据。
2.根据权利要求1所述的基于宽频量测的电力系统间谐波潮流获取方法,其特征在于,所述步骤1中,所述振荡电流数据为电网实际量测获得的振荡电流数据,或者通过仿真获得的振荡电流数据,具体表示为:
i(t)=[i1(t),i2(t),...in(t)]T
其中,in(t)为第k个边界节点测得的电流在t时刻的值,n为构建电力系统的网络结构的边界节点总数。
3.根据权利要求1所述的基于宽频量测的电力系统间谐波潮流获取方法,其特征在于,所述步骤2中,通过prony分析,对振荡电流数据进行解耦处理,获得电力系统中主要稳定的次/超同步振荡电流频率与幅值,具体包括:
利用prony分析对振荡电流数据进行解耦,获得次/超同步振荡的f,即为间谐波频率;解耦得到的次/超同步频率f满足以下条件:
1)在电力系统中保持稳定振荡,即在电力中系统稳定等幅的间谐波振荡;
2)三相振荡符合对称且正序的条件,若不符合,则需进行序分量分解后,再转换为对称平衡的间谐波电流数据,进行分别单独计算后对结果进行求和。
4.根据权利要求1所述的基于宽频量测的电力系统间谐波潮流获取方法,其特征在于,在所述步骤3中,根据振荡频率分量与电力系统的网络结构参数,构建某一振荡频率分量下的系统节点导纳矩阵,具体包括:
对于一个有n个节点的电力系统,其系统节点导纳矩阵如下所示:
式中,Yij为节点i与节点j之间的互导纳,Yii为节点i的自导纳;
其中,i,j=1、2、3...n,i≠j;
将电力系统中任意一段线路的元件模型等效为2类,分别为:
第一等效模型,其自导纳与互导纳构建如下:
第二等效模型,其自导纳与互导纳构建如下:
在上述第一等效模型、第二等效模型的基础上,对受频率影响的电力系统元件进行基于频率的校正,具体校正过程包括:
1)对于一条π型等值电路,校正公式如下:
式中,zl为线路阻抗;yl为线路对地导纳;r、x、g和b分别为基频下线路电阻、电抗、电导和电纳;f和fn分别为间谐波频率和基频频率;
即对于π型等值电路,适用上述第一等效模型,zij=zl、yi=yj=yl/2;
2)对于并联RLC恒阻抗负荷,校正公式如下:
式中:R、Xc和Xl是并联RLC恒阻抗负荷的等效电阻、电容和电感;Vn是负荷额定电压;Pn、Qc和Ql是负荷有功功率、容性无功功率和感性有功功率;
即对于并联RLC恒阻抗负荷,适用第二等效模型,yi=R//Xc//Xl;
3)对于恒功率负荷,保持视在功率S恒定不变为最强等量约束,其建模满足以下公式:
即对于并联RLC恒阻抗负荷,适用上述第二等效模型,yi=If/Uf;
4)对于双绕组变压器而言,其计算公式如下式:
即对于双绕组变压器,适用第一等效模型,zij=zb、yi=yj=0;
5)对于同步发电机产生的间谐波电流远小于基频频率的电流,在潮流计算过程中将同步发电机视为负荷处理:
Z=R+jX
式中:Z为同步发电机间谐波频率下等效阻抗,R和X为同步发电机间谐波频率下等效电阻与等效电抗;
5.根据权利要求4所述的基于宽频量测的电力系统间谐波潮流获取方法,其特征在于,在所述步骤4中,网络参数的完整判断的依据在于是否系统在间谐波频率分量下的系统等效导纳矩阵是否能够建立,即构建的系统节点导纳矩阵是否有不可计算导纳Yij,使节点导纳矩阵不完整。
6.根据权利要求4所述的基于宽频量测的电力系统间谐波潮流获取方法,其特征在于,在所述步骤5中,通过量测未知网络参数区域与已知区域联络线传输电流参数或者电压参数,通过prony分析后,将未知网络参数区域等效为电流源或电压源,用于间谐波潮流计算。
7.根据权利要求1所述的基于宽频量测的电力系统间谐波潮流获取方法,其特征在于,前述步骤7中,将接地节点间谐波电压设置为零,系统内部无间谐波源,即所有节点上各支路间谐波电流相量合为零,边界节点等效间谐波电流源与间谐波电压源通过步骤6量测已知,可设置点位系统的间谐波电流与间谐波电压初值。
8.根据权利要求1所述的基于宽频量测的电力系统间谐波潮流获取方法,其特征在于,利用间谐波潮流计算公式进行计算,获得次/超同步振荡的间谐波在电网系统中传播路径与分布数据,包括:
间谐波潮流计算选择使用电流平衡方程,其具体公式如下:
In=YnVn
式中:In为节点流入或流出的某一频率间谐波电流;Yn为对应频率下的系统节点导纳矩阵;Vn为节点在对应频率的间谐波电压;
对于一个完整可知的系统,已知间谐波振荡源的流入电流,所有边界节点的间谐波电压为零,所有接地节点的间谐波电流为零,上述电流平衡方程转换成以下形式:
式中:IB、II和IG分别是系统边界节点、内部节点与接地节点的电流,UB、UI和UG分别是系统内部节点与边界节点的电压,YBB、YII、YGG、YBI、YIB、YGI和YIG分别是系统边界、内部、外部的自导纳,以及边界与内部、内部与边界、接地与内部、内部与接地之间的互导纳;
其中,接地节点电压UG和内部节点电流II等于零,设定边界电流IB已知,可由以下公式计算出边界电压UB:
即在可知系统中,仅通过对间谐波振荡源与可知系统连接线路的电流量测,计算得到间谐波振荡源在整个可知系统的传播与分布。
9.一种基于宽频量测的电力系统间谐波潮流获取的计算机系统,其特征在于,包括:
一个或多个处理器;
存储器,存储可被操作的指令,所述指令在通过所述一个或多个处理器执行时使得所述一个或多个处理器执行操作,所述操作包括执行前述权利要求1-8中任意一项所述的方法的过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111358435.XA CN114301055B (zh) | 2022-02-17 | 2022-02-17 | 基于宽频量测的电力系统间谐波潮流获取方法与系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111358435.XA CN114301055B (zh) | 2022-02-17 | 2022-02-17 | 基于宽频量测的电力系统间谐波潮流获取方法与系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114301055A true CN114301055A (zh) | 2022-04-08 |
CN114301055B CN114301055B (zh) | 2024-07-12 |
Family
ID=80965358
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111358435.XA Active CN114301055B (zh) | 2022-02-17 | 2022-02-17 | 基于宽频量测的电力系统间谐波潮流获取方法与系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114301055B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114971063A (zh) * | 2022-06-13 | 2022-08-30 | 浙江大学 | 一种城市供区多端低频互联系统的频率选择方法 |
CN116822436A (zh) * | 2023-06-30 | 2023-09-29 | 四川大学 | 一种直流输电送端交流系统的振荡风险灵敏度分析方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1114675A (ja) * | 1997-06-19 | 1999-01-22 | Chubu Electric Power Co Inc | 電力系統の高調波測定方法 |
CN103647276A (zh) * | 2013-12-10 | 2014-03-19 | 国家电网公司 | 一种电能质量预警系统及其方法 |
US20170184640A1 (en) * | 2014-09-12 | 2017-06-29 | Carnegie Mellon University | Systems, Methods, and Software for Planning, Simulating, and Operating Electrical Power Systems |
CN107968402A (zh) * | 2017-12-11 | 2018-04-27 | 华北电力大学 | 一种多电流源型振荡源的定位方法 |
CN110703008A (zh) * | 2019-09-20 | 2020-01-17 | 中国电力科学研究院有限公司 | 一种宽频测量数据的站域处理方法及装置 |
US20200106263A1 (en) * | 2017-11-09 | 2020-04-02 | Tsinghua University | System-level protection system and method for sub/super-synchronous resonance/oscillation |
CN111239489A (zh) * | 2018-11-29 | 2020-06-05 | 南京南瑞继保电气有限公司 | 一种结合prony与fft算法的次同步振荡分析方法 |
CN112688321A (zh) * | 2021-01-11 | 2021-04-20 | 河海大学 | 大电网次/超同步振荡路径获取方法及系统 |
CN112838670A (zh) * | 2020-12-31 | 2021-05-25 | 广州兆和电力技术有限公司 | 基于宽频测量装置的电能治理系统及谐波提取方法 |
CN113364009A (zh) * | 2021-06-22 | 2021-09-07 | 国电南瑞科技股份有限公司 | 一种含大规模新能源电网次同步振荡广域监测和预警方法 |
US20220021211A1 (en) * | 2020-07-15 | 2022-01-20 | North China Electric Power University | Method, system and storage medium for suppressing sub/super-synchronous oscillation for direct-drive wind turbine based on energy compensation |
-
2022
- 2022-02-17 CN CN202111358435.XA patent/CN114301055B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1114675A (ja) * | 1997-06-19 | 1999-01-22 | Chubu Electric Power Co Inc | 電力系統の高調波測定方法 |
CN103647276A (zh) * | 2013-12-10 | 2014-03-19 | 国家电网公司 | 一种电能质量预警系统及其方法 |
US20170184640A1 (en) * | 2014-09-12 | 2017-06-29 | Carnegie Mellon University | Systems, Methods, and Software for Planning, Simulating, and Operating Electrical Power Systems |
US20200106263A1 (en) * | 2017-11-09 | 2020-04-02 | Tsinghua University | System-level protection system and method for sub/super-synchronous resonance/oscillation |
CN107968402A (zh) * | 2017-12-11 | 2018-04-27 | 华北电力大学 | 一种多电流源型振荡源的定位方法 |
CN111239489A (zh) * | 2018-11-29 | 2020-06-05 | 南京南瑞继保电气有限公司 | 一种结合prony与fft算法的次同步振荡分析方法 |
CN110703008A (zh) * | 2019-09-20 | 2020-01-17 | 中国电力科学研究院有限公司 | 一种宽频测量数据的站域处理方法及装置 |
WO2021051656A1 (zh) * | 2019-09-20 | 2021-03-25 | 中国电力科学研究院有限公司 | 宽频测量数据的处理方法、装置、电子设备及介质 |
US20220021211A1 (en) * | 2020-07-15 | 2022-01-20 | North China Electric Power University | Method, system and storage medium for suppressing sub/super-synchronous oscillation for direct-drive wind turbine based on energy compensation |
CN112838670A (zh) * | 2020-12-31 | 2021-05-25 | 广州兆和电力技术有限公司 | 基于宽频测量装置的电能治理系统及谐波提取方法 |
CN112688321A (zh) * | 2021-01-11 | 2021-04-20 | 河海大学 | 大电网次/超同步振荡路径获取方法及系统 |
CN113364009A (zh) * | 2021-06-22 | 2021-09-07 | 国电南瑞科技股份有限公司 | 一种含大规模新能源电网次同步振荡广域监测和预警方法 |
Non-Patent Citations (2)
Title |
---|
BO GAO: "Identifying_and_Ranking_Sources_of_SSR_Based_on_the_Concept_of_Subsynchronous_Power", IEEE TRANSACTIONS ON POWER DELIVERY, vol. 35, no. 1, 14 May 2019 (2019-05-14), pages 258 - 268, XP011767469, DOI: 10.1109/TPWRD.2019.2916848 * |
樊陈,等: "电网宽频振荡实时监测技术方案", 电力系统自动化, vol. 45, no. 11, 10 June 2021 (2021-06-10), pages 152 - 159 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114971063A (zh) * | 2022-06-13 | 2022-08-30 | 浙江大学 | 一种城市供区多端低频互联系统的频率选择方法 |
CN116822436A (zh) * | 2023-06-30 | 2023-09-29 | 四川大学 | 一种直流输电送端交流系统的振荡风险灵敏度分析方法 |
CN116822436B (zh) * | 2023-06-30 | 2024-02-27 | 四川大学 | 一种直流输电送端交流系统的振荡风险灵敏度分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114301055B (zh) | 2024-07-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ren et al. | A refined frequency scan approach to sub-synchronous control interaction (SSCI) study of wind farms | |
CN108667048B (zh) | 新能源并网系统振荡稳定性的频域判稳方法及装置 | |
Liu et al. | Impedance network modeling and quantitative stability analysis of sub-/super-synchronous oscillations for large-scale wind power systems | |
Baran et al. | Fault analysis on distribution feeders with distributed generators | |
CN114301055B (zh) | 基于宽频量测的电力系统间谐波潮流获取方法与系统 | |
CN106250660A (zh) | 一种海上风电场谐振分析方法 | |
Ishchenko et al. | Dynamic equivalencing of distribution networks with dispersed generation using Hankel norm approximation | |
Lei et al. | Industrial approaches for dynamic equivalents of large power systems | |
CN112202211B (zh) | 一种基于模态分析法的多vsg并网功频振荡评估方法 | |
Sowa et al. | Impedance-based analysis of harmonic resonances in HVDC connected offshore wind power plants | |
Pecenak et al. | Inversion reduction method for real and complex distribution feeder models | |
CN114884094A (zh) | 一种风电场阻抗特性监测方法、系统、设备和介质 | |
Ovalle et al. | GICharm: A system level analysis tool for geomagnetic disturbance related harmonics | |
CN112072692B (zh) | 一种新能源发电场站的阻抗等值方法及装置 | |
CN109638871B (zh) | 考虑风电接入的大规模交直流混联系统主网划分方法 | |
CN115270597A (zh) | 配电网继电保护适应性仿真测试方法、仿真模型和装置 | |
Liu et al. | Generalized Operational Short-Circuit Ratio for Grid Strength Assessment in Renewable Power Systems under Non-Rated Operating Conditions | |
CN109033660A (zh) | 不对称故障分析方法及装置 | |
Magalhães et al. | Parametric regression in synchronous and induction generators | |
CN112688321A (zh) | 大电网次/超同步振荡路径获取方法及系统 | |
Acevedo et al. | Dynamic voltage stability comparison of thermal and wind power generation with different static and dynamic load models | |
Zhang et al. | Identification method of all-operating-point admittance model for wind farms considering frequency-coupling characteristics | |
Zhang | Analysis and control of resonances in HVDC connected DFIG-based offshore wind farm | |
CN116384152B (zh) | 一种风电场场站的三相电压不平衡适应性测试方法及装置 | |
Zaker et al. | Active Distribution Networks Equivalencing for Dynamic Studies Using Measurements |
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 |