CN116822436B - 一种直流输电送端交流系统的振荡风险灵敏度分析方法 - Google Patents
一种直流输电送端交流系统的振荡风险灵敏度分析方法 Download PDFInfo
- Publication number
- CN116822436B CN116822436B CN202310796506.7A CN202310796506A CN116822436B CN 116822436 B CN116822436 B CN 116822436B CN 202310796506 A CN202310796506 A CN 202310796506A CN 116822436 B CN116822436 B CN 116822436B
- Authority
- CN
- China
- Prior art keywords
- network
- equivalent
- node
- branch
- power
- 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
- 230000010355 oscillation Effects 0.000 title claims abstract description 57
- 230000005540 biological transmission Effects 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000010206 sensitivity analysis Methods 0.000 title claims abstract description 17
- 230000035945 sensitivity Effects 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims description 13
- 239000003990 capacitor Substances 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000004804 winding Methods 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000004146 energy storage Methods 0.000 claims description 3
- 230000001939 inductive effect Effects 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 7
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 12
- 238000004891 communication Methods 0.000 description 3
- 230000032683 aging Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002847 impedance measurement Methods 0.000 description 2
- 238000009413 insulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000008569 process 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/39—Circuit design at the physical level
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
-
- 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Abstract
本发明公开一种直流输电送端交流系统的振荡风险灵敏度分析方法,首先将直流输电送端交流网络发生的振荡进行分类,分为局部振荡和非局部振荡;根据实际工程直流送端交流网架结构,确定并搭建等值交流网简化络模型;根据实际情况和等值交流网简化络模型,建立电力系统复杂网络加权网络模型;定义振荡风险特征指标,通过基于电力系统复杂网络加权网络模型的振荡风险指标能进行灵敏度分析。本发明提出的振荡风险特征指标,只需建立电力系统的复杂网络加权模型,通过灵敏度算法即可得出系统网络结构和参数对振荡风险的灵敏度关系;易于实施且实用性强,更有利于工程实践。
Description
技术领域
本发明涉及直流输电送端交流系统领域,具体为一种基于复杂网络理论的直流输电送端交流振荡风险灵敏度分析方法。
背景技术
随着大量电力电子设备大量接入系统,包括电源侧风电、光伏等可再生能源,输电网侧大量常规和柔性直流以及负荷侧各类电力电子负荷等,交直流系统的谐波特性日益复杂。电网安全稳定运行中面临的谐波问题,已由过去主要解决低压侧谐波超标电能质量问题,转变为高低压侧都要面对的谐波电流超标甚至谐振放大问题。谐波由非线性设备产生,在电力网络传播过程中,往往会因网络元件参数与谐波源参数相匹配引发谐波的放大、谐振现象。谐波谐振将导致节点谐波过电压,同时放大支路谐波电流,导致系统母线节点谐波电压与支路谐波电流的严重畸变,进而造成电气设备的绝缘损坏;增大变压器、旋转电机等设备的损耗;加速电容器绝缘老化,甚至引发爆炸;引起系统内继电保护和自动装置误动或拒动;干扰通信信号等危害。
目前针对谐波谐振的分析方法主要有四类:谐波电压电流比值法、频率扫描分析法、s域传递函数法及谐振模态分析法。谐振模态分析法是目前最有效的提取谐振信息的方法,该方法基于矩阵特征值分析方法,将系统的节点导纳矩阵逆矩阵转换到解耦的模态域进行分析,可以对每个谐振频率下存在的容性元件和感性元件之间的复杂能量交互联系进行解耦分析。但是随着新能源等大功率非线性且具有随机性质的元件和负荷的接入,电力系统的规模越来越大、节点越来越多,导致谐波谐振分析的难度和复杂度进一步加大,需要更加快速的谐波谐振分析算法;同时,电力系统中存在着多种多样的不确定性,如电力需求预测,可再生能源出力,负荷波动,系统的随机行为如系统操作与故障,以及由于制造误差和设备材料特性受温度、老化等影响导致的设备参数不确定性等,会影响系统潜在谐波谐振规律。
现在工程中主要应用的特征值敏感度分析法与谐振模态阻抗分析法,需要对导纳矩阵进行提取并利用局部灵敏度分析,较为繁琐,且对有多维输入以及各种不确定性因素的实际大系统适用性不强。
发明内容
针对上述问题,本发明提出了一种直流输电送端交流系统的振荡风险灵敏度分析方法,提出振荡风险特征指标,通过建立电力系统的复杂网络加权模型,根据所提出的灵敏度算法即可得出系统网络结构和参数对振荡风险的灵敏度关系。技术方案如下:
一种直流输电送端交流系统的振荡风险灵敏度分析方法,包括以下步骤:
步骤1:将直流输电送端交流网络发生的振荡进行分类,分为直流系统引起交流网络振荡,即局部振荡,以及交流网络自身振荡,即非局部振荡;
步骤2:根据实际工程直流送端交流网架结构,确定并搭建等值交流网简化网络模型;
步骤3:根据实际情况和等值交流网简化网络模型,采用基于传输阻抗的电气距离计算方法建立系统无权网络模型,再采用相似权对网络中的边和节点进行加权,建立电力系统复杂网络加权网络模型;
步骤4:定义振荡风险特征指标,通过基于电力系统复杂网络加权网络模型的振荡风险指标能进行灵敏度分析,振荡风险特征指标越小,则在对应网络结构或者参数的情况下对振荡的灵敏度越低,反之则越高。
进一步的,所述等值交流网简化网络模型包括:
1)变压器π型等值简化模型
对传统双绕组变压器进行π型等效化简后,变压器两端节点i,j间的串联阻抗为ZLeq,节点i,j侧的并联导纳值分别为YLeqi和YLeqj;
其中,tap为变压器分接头参数,tk为变压器变比;ZT为变压器一次绕组侧等效阻抗,且
其中,P表示变压器额定容量,U为一次侧电压,根据折算方向进行赋值;RT为一次侧电阻,XT为一次侧电抗;
2)线路π型等值简化模型
参照BPA参数对模型进行等效型,线路中各等效参数Rcir、Lcir和Ccir的关系式为:
其中,XLcir为电抗标幺值,Rcir为电阻标幺值,f为运行频率,Lcir为电感标幺值,Bcir为电纳标幺值,Ccir为电容标幺值;
3)无功补偿支路等值简化模型
对于无源无功补偿装置,化简时除考虑额定补偿无功Q1N和无功补偿器额定电压U1N外,还考虑无功补偿潮流数据;简化后并联电感支路感抗为:
并联电容支路容抗为:
其中,ΔQ1t和ΔU1t为负荷潮流数据修正量;
4)发电机等值简化模型
对于电压源串电阻形式的发电机,建立导纳矩阵时无需关注电源阻抗,直接忽略电源支路,发电机支路等效串联电阻为:
发电机支路等效串联电抗为:
其中,Ra为电源标幺电阻,Xd为标幺电抗,SBN为电源额定出力,ΔSBt为潮流数据修正的功率出力,Sn为标幺基准功率;
5)负荷支路等值简化模型
对于稳定运行状态下的恒功率负荷,额定参数下吸收的功率为定值,负荷支路等效为电阻加电感或电阻加电容的形式,有功支路决定等效电阻的大小,无功支路则用电感或电容此类储能元件进行替换;并联电阻支路电阻为:
负荷并联电抗支路感抗为:
并联电容支路容抗为:
其中,U2N、P2N、Q2N为负荷的额定参数,ΔQ2t、ΔP2t、ΔU2t为负荷潮流数据修正量;
6)LCC换流器的交流侧等值模型
使用扫频法获取系统阻抗信息;若在系统中存在一个频率为ωp的谐波分量,则通过测量元件在对应频率下响应的谐波电压及谐波电流,即计算出元件在该点处的谐波阻抗;换流器在ωp处的交流阻抗/导纳由复数电压与复数电流之比来获得,表达式如下:
其中,Zac(ωp)为换流器在ωp处的交流阻抗,Δugac(ωp)为换流器在ωp处的复数电压,Δigac(ωp)为换流器在ωp处的复数电流,Yac(ωp)为换流器在ωp处的导纳。
更进一步的,所述步骤3具体包括:
步骤3.1:采用基于传输阻抗的电气距离计算方法建立系统无权网络模型
定义节点i到达节点j所经过的边数即为节点i到节点j的距离;若节点i与节点j之间存在功率传输关系,则认为两节点之间构成一条支路,称为广义支路;该广义支路的阻抗zi表征功率由节点i传至节点j引起的功率损耗,称为等效传输阻抗;
节点之间有功率传输关系时,传输阻抗为有限值;若无传输关系,传输阻抗为无穷大;
取传输阻抗模|zi|作为节点之间的电气距离,电气距离z的计算公式如下:
其中,Z为节点间的等效阻抗,ΔU为节点间的电压差,ΔS为节点间的功率差;
步骤3.2:采用相似权对网络中的边和节点进行加权,建立系统加权网络模型;
定义某边的权重wi为该边所传输电能的容量Sij的值与所定义的整个电网的额定容量SN的比,该边所传输电能的容量Sij大小为节点i、j之间的所有连接边的额定容量之和减去损耗的能量的绝对值,潮流方向为从i到j时为正值,反之为负值;定义某节点的权重si为该节点当前所传输的容量Si与所定义的整个电网的额定容量SN之比;该节点当前所传输的容量Si为所有为正值的Sij之和;边权重w和节点权重s的计算公式如式下:
更进一步的,所述步骤4中振荡风险特征指标为:
其中,i、j表示节点标号;ΔVi为i节点电压偏移量、ΔIij为ij支路电流偏移量;k、l分别为网络结构系数和元件参数系数,且满足k+l=1;smax为复杂网络中最大节点权重,zmax为复杂网络中最大电气距离值,wij为以i节点和j节点为端点的边的边权重,zij为i节点和j节点间的电气距离值。
本发明的有益效果是:本发明提出的振荡风险特征指标,只需建立电力系统的复杂网络加权模型,通过灵敏度算法即可得出系统网络结构和参数对振荡风险的灵敏度关系;相比于主流的特征值敏感度分析法与谐振模态阻抗分析法中需要进行导纳矩阵的提取并利用局部灵敏度分析,本发明易于实施且实用性强,更有利于工程实践。
附图说明
图1为本发明直流振荡风险灵敏度分析算法流程图。
图2为IEEE 9节点测试系统图。
图3为变压器等效简化模型。
图4为线路π型等效简化模型。
图5为无功补偿支路等值简化模型。
图6为发电机等效简化模型。
图7为负荷等效简化模型。
图8为等效阻抗测量电路图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步详细说明。本发明提出了一种振荡风险特征指标,通过建立电力系统的复杂网络加权模型,根据所提出的灵敏度算法即可得出系统网络结构和参数对振荡风险的灵敏度关系。本发明建立直流送出交流系统的网络结构和参数在交流网络N-1、N-2等不同工况下的直流振荡风险灵敏度分析算法,具体步骤见图1所示。本发明基于IEEE 9节点测试系统(如图2)对方案进行详细说明。
本发明采用的技术方案如下:
步骤1:将直流输电送端交流网络发生的振荡进行分类,分为直流系统引起交流网络振荡,即局部振荡;交流网络自身振荡,定义为非局部振荡。
步骤2:根据实际工程直流送端交流网架结构,确定并搭建等值交流网络模型。本发明中提出的等值方法基于BPA软件参数类型对交流系统中各常见元件进行分类,以下分别对元件等效方式进行介绍:
1)变压器π型等值简化模型
本发明以双绕组变压器为例进行化简,对传统双绕组变压器进行π型等效化简。基本结构如图3所示。
式(1)表示变压器i,j节点间的串联阻抗:
式(2)分为变压器节点i,j侧的并联导纳值:
由于折算方向不同,因此两侧节点的导纳值需使用不同的计算公式进行计算。其中:
其中,tk为变压器变比,P表示变压器额定容量,U为一次侧电压可根据折算方向进行赋值,tap为变压器分接头参数,ZT为变压器一次绕组侧等效阻抗。
2)线路π型等值简化模型
对于基频下的线路,单个π型等效模型基本满足精度要求。且BPA数据卡中线路的初始数据,仅包含线路标幺电阻,标幺电抗和标幺电纳,因此基于BPA进行线路等效导纳求取时,线路无需使用复杂的结构,可完全参照BPA参数对模型进行等效。如图4所示为线路的等效简化模型。
如式(4)所示,根据公式可计算得到线路中各等效参数Rcir、Lcir和Ccir的关系式为:
其中,XLcir为电抗标幺值,Rcir为电阻标幺值,f为运行频率,Lcir为电感标幺值,Bcir为电纳标幺值,Ccir为电容标幺值。
3)无功补偿支路等值简化模型
本发明以无源无功补偿装置为例,化简时除考虑额定补偿无功Q1N和无功补偿器额定电压U1N外,还考虑了无功补偿潮流数据,该式中ΔQ1t,ΔU1t为潮流修正量,通过潮流数据与系统额定数据做差对比得到。如图5所示为无功补偿支路等值简化模型。
简化后并联电感支路感抗为:
并联电容支路容抗为:
4)发电机等值简化模型
本发明以电压源串电阻形式的发电机为例进行等效,电源进行处理时,若仅需关注母线节点信息,可以直接忽略电源阻抗值,即建立导纳矩阵时无需关注电源阻抗,直接忽略电源支路。其中Ra为电源标幺电阻,Xd为标幺电抗,SBN为电源额定出力,ΔSBt为潮流数据修正的功率出力,Sn为标幺基准功率。发电机等效简化模型如图6所示。
发电机支路等效串联电阻:
发电机支路等效串联电抗:
5)负荷支路等值简化模型
本文以稳定运行状态下的恒功率负荷为例,此时额定参数下吸收的功率为定值,负荷支路可以等效为电阻加电感或电容的形式,有功支路决定等效电阻的大小,无功支路则可以用电感或电容此类储能元件进行替换,可简化为图7所示的负荷等效简化模型。
并联电阻支路电阻如式(9)所示:
式(10)表示负荷并联电抗支路感抗计算公式:
并联电容支路容抗计算如(11)所示:
其中,U2N、P2N、Q2N为负荷的额定参数,ΔQ2t、ΔP2t、ΔU2t为负荷潮流数据修正量。
6)LCC换流器的交流侧等值模型
由于复杂交流系统导纳矩阵构建主要服务于系统稳定性分析,因此多数工程中会存在有多个复杂的换流器,本发明中使用扫频法获取系统阻抗信息,该方法无需关注换流器内部结构,可直接实现对此类元件的导纳获取。在交流系统中,任意三相系统等效阻抗测量如电路图如图8所示。
若在系统中存在一个频率为ωp的谐波分量,则通过测量元件在对应频率下响应的谐波电压及谐波电流,即可计算出元件在该点处的谐波阻抗。如图8所示,换流器在ωp下的交流阻抗/导纳可由复数电压与复数电流之比来获得,表达式如式(12)所示:
步骤3:根据实际情况和化简后的等值模型,建立电力系统复杂网络加权模型。
1)采用基于传输阻抗的电气距离计算方法建立系统无权网络模型
电气距离表示节点间的电气耦合度,定义节点i到达节点j所经过的边数即为节点i到节点j的距离。若节点i与节点j之间存在功率传输关系,则认为两节点之间构成一条支路,称为广义支路。该广义支路的阻抗zi可以表征功率由节点i传至节点j引起的功率损耗,称为等效传输阻抗。节点之间有功率传输关系时,传输阻抗为有限值;若无传输关系,传输阻抗为无穷大。取传输阻抗模|zi|作为节点之间的电气距离,电气距离z的计算公式如式(13)所示:
2)采用相似权对网络中的边和节点进行加权,建立系统加权网络模型。
边的权重wi为该边所传输电能的容量Sij的值与所定义的整个电网的额定容量SN的比,Sij大小为节点i、j之间的所有连接边的额定容量之和减去损耗的能量的绝对值,潮流方向为从i到j时为正值,反之为负值。节点的权重si定义为该节点当前所传输的容量Si与所定义的整个电网的额定容量SN之比;Si为所有为正值的Sij之和。边权重w和节点权重s的计算公式如式(14)、(15)所示:
步骤4:定义振荡风险特征指标,通过基于加权复杂网络模型的振荡风险指标能进行灵敏度分析。
定义振荡风险特征指标为W,其计算公式如(16)所示:
其中i、j表示节点标号;ΔVi为i节点电压偏移量、ΔIij为ij支路电流偏移量;k、l分别为网络结构系数和元件参数系数,且满足k+l=1;smax为复杂网络中最大节点权重,zmax为复杂网络中最大电气距离值。
振荡风险特征指标第一部分主要是表明直流馈入交流系统在N-1,N-2等不同网络结构下对振荡风险的影响,如果关注重点在网络结构的影响,则k的值大,l的值小。这是由于线路断开后,对节点的权值影响最大,节点的电压偏移表现最明显。
第二部分主要是表明系统在各参数变换的情况下对振荡风险的影响。如果关注重点在元件参数的影响,则k的值小,l的值大。这是由于元件参数变化后,对边的权值影响最大,支路电流偏移表现最明显。
若发生振荡,节点电压和支路电流的偏移必定有所反馈,因此通过基于加权复杂网络模型的振荡风险指标能进行灵敏度分析。振荡风险特征指标越小,说明在这种网络结构或者参数的情况下对振荡的灵敏度越低,反之则是越高。
Claims (2)
1.一种直流输电送端交流系统的振荡风险灵敏度分析方法,其特征在于,包括以下步骤:
步骤1:将直流输电送端交流网络发生的振荡进行分类,分为直流系统引起交流网络振荡,即局部振荡,以及交流网络自身振荡,即非局部振荡;
步骤2:根据实际工程直流送端交流网架结构,确定并搭建等值交流网简化网络模型;
步骤3:根据实际情况和等值交流网简化网络模型,采用基于传输阻抗的电气距离计算方法建立系统无权网络模型,再采用相似权对网络中的边和节点进行加权,建立电力系统复杂网络加权网络模型;
步骤4:定义振荡风险特征指标,通过基于电力系统复杂网络加权网络模型的振荡风险指标能进行灵敏度分析,振荡风险特征指标越小,则在对应网络结构或者参数的情况下对振荡的灵敏度越低,反之则越高;
所述等值交流网简化网络模型包括:
1)变压器π型等值简化模型
对传统双绕组变压器进行π型等效化简后,变压器两端节点i,j间的串联阻抗为ZLeq,节点i,j侧的并联导纳值分别为YLeqi和YLeqj;
其中,tap为变压器分接头参数,tk为变压器变比;ZT为变压器一次绕组侧等效阻抗,且
其中,P表示变压器额定容量,U为一次侧电压,根据折算方向进行赋值;RT为一次侧电阻,XT为一次侧电抗;
2)线路π型等值简化模型
参照BPA参数对模型进行等效型,线路中各等效参数Rcir、Lcir和Ccir的关系式为:
其中,XLcir为电抗标幺值,Rcir为电阻标幺值,f为运行频率,Lcir为电感标幺值,Bcir为电纳标幺值,Ccir为电容标幺值;
3)无功补偿支路等值简化模型
对于无源无功补偿装置,化简时除考虑额定补偿无功Q1N和无功补偿器额定电压U1N外,还考虑无功补偿潮流数据;简化后并联电感支路感抗为:
并联电容支路容抗为:
其中,ΔQ1t和ΔU1t为负荷潮流数据修正量;
4)发电机等值简化模型
对于电压源串电阻形式的发电机,建立导纳矩阵时无需关注电源阻抗,直接忽略电源支路,发电机支路等效串联电阻为:
发电机支路等效串联电抗为:
其中,Ra为电源标幺电阻,Xd为标幺电抗,SBN为电源额定出力,ΔSBt为潮流数据修正的功率出力,Sn为标幺基准功率;
5)负荷支路等值简化模型
对于稳定运行状态下的恒功率负荷,额定参数下吸收的功率为定值,负荷支路等效为电阻加电感或电阻加电容的形式,有功支路决定等效电阻的大小,无功支路则用电感或电容此类储能元件进行替换;并联电阻支路电阻为:
负荷并联电抗支路感抗为:
并联电容支路容抗为:
其中,U2N、P2N、Q2N为负荷的额定参数,ΔQ2t、ΔP2t、ΔU2t为负荷潮流数据修正量;
6)LCC换流器的交流侧等值模型
使用扫频法获取系统阻抗信息;若在系统中存在一个频率为ωp的谐波分量,则通过测量元件在对应频率下响应的谐波电压及谐波电流,即计算出元件在测量点处的谐波阻抗;换流器在ωp处的交流阻抗/导纳由复数电压与复数电流之比来获得,表达式如下:
其中,Zac(ωp)为换流器在ωp处的交流阻抗,Δugac(ωp)为换流器在ωp处的复数电压,Δigac(ωp)为换流器在ωp处的复数电流,Yac(ωp)为换流器在ωp处的导纳;
所述步骤3具体包括:
步骤3.1:采用基于传输阻抗的电气距离计算方法建立系统无权网络模型
定义节点i到达节点j所经过的边数即为节点i到节点j的距离;若节点i与节点j之间存在功率传输关系,则认为两节点之间构成一条支路,称为广义支路;该广义支路的阻抗zi表征功率由节点i传至节点j引起的功率损耗,称为等效传输阻抗;
节点之间有功率传输关系时,传输阻抗为有限值;若无传输关系,传输阻抗为无穷大;
取传输阻抗模|zi|作为节点之间的电气距离,电气距离z的计算公式如下:
其中,Z为节点间的等效阻抗,ΔU为节点间的电压差,ΔS为节点间的功率差;
步骤3.2:采用相似权对网络中的边和节点进行加权,建立系统加权网络模型;
定义某边的权重wi为该边所传输电能的容量Sij的值与所定义的整个电网的额定容量SN的比,该边所传输电能的容量Sij大小为节点i、j之间的所有连接边的额定容量之和减去损耗的能量的绝对值,潮流方向为从i到j时为正值,反之为负值;定义某节点的权重si为该节点当前所传输的容量Si与所定义的整个电网的额定容量SN之比;该节点当前所传输的容量Si为所有为正值的Sij之和;边权重w和节点权重s的计算公式如式下:
2.根据权利要求1所述的直流输电送端交流系统的振荡风险灵敏度分析方法,其特征在于,所述步骤4中振荡风险特征指标为:
其中,i、j表示节点标号;ΔVi为i节点电压偏移量、ΔIij为ij支路电流偏移量;k、l分别为网络结构系数和元件参数系数,且满足k+l=1;smax为复杂网络中最大节点权重,zmax为复杂网络中最大电气距离值,wij为以i节点和j节点为端点的边的边权重,zij为i节点和j节点间的电气距离值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310796506.7A CN116822436B (zh) | 2023-06-30 | 2023-06-30 | 一种直流输电送端交流系统的振荡风险灵敏度分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310796506.7A CN116822436B (zh) | 2023-06-30 | 2023-06-30 | 一种直流输电送端交流系统的振荡风险灵敏度分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116822436A CN116822436A (zh) | 2023-09-29 |
CN116822436B true CN116822436B (zh) | 2024-02-27 |
Family
ID=88119953
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310796506.7A Active CN116822436B (zh) | 2023-06-30 | 2023-06-30 | 一种直流输电送端交流系统的振荡风险灵敏度分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116822436B (zh) |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104638678A (zh) * | 2015-03-16 | 2015-05-20 | 国网山西省电力公司电力科学研究院 | 一种多风机并网谐振分析方法 |
CN105160058A (zh) * | 2015-07-10 | 2015-12-16 | 国网四川省电力公司电力科学研究院 | 一种用于电网电磁暂态计算的网络边界等值方法 |
CN106021682A (zh) * | 2016-05-13 | 2016-10-12 | 清华大学 | 基于阻抗网络模型的次同步振荡分析方法和装置 |
CN107681682A (zh) * | 2017-10-25 | 2018-02-09 | 国家电网公司 | 一种基于ward等值的交直流系统等值方法 |
CN111597680A (zh) * | 2020-04-09 | 2020-08-28 | 华北电力大学 | 基于模态频率处电气阻尼的次同步谐振风险量化评估方法 |
CN114301055A (zh) * | 2022-02-17 | 2022-04-08 | 河海大学 | 基于宽频量测的电力系统间谐波潮流获取方法与系统 |
CN114400667A (zh) * | 2021-11-15 | 2022-04-26 | 国网冀北电力有限公司智能配电网中心 | 交直流混合配电网安全风险评估方法 |
CN114966203A (zh) * | 2022-06-16 | 2022-08-30 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 电力系统谐波振荡的检测方法、装置和计算机设备 |
CN115021256A (zh) * | 2022-08-08 | 2022-09-06 | 四川大学 | 一种大规模交直流输电系统电磁暂态模型自动生成方法 |
WO2022233079A1 (zh) * | 2021-05-07 | 2022-11-10 | 南方电网科学研究院有限责任公司 | 适用于柔性直流的串联型无源阻抗适配器参数设计方法 |
CN115549093A (zh) * | 2022-11-24 | 2022-12-30 | 中国电力科学研究院有限公司 | 一种新能源电力系统在线建模与振荡分析的方法及系统 |
CN115833124A (zh) * | 2022-12-28 | 2023-03-21 | 国网浙江省电力有限公司电力科学研究院 | 换流站直流近区谐振风险评估方法、装置、设备及介质 |
CN116054187A (zh) * | 2022-04-27 | 2023-05-02 | 国网上海市电力公司 | 一种风电经柔直并网产生的次同步振荡抑制方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3971167B2 (ja) * | 2001-11-20 | 2007-09-05 | 株式会社ルネサステクノロジ | 等価回路の導出方法、および、そのためのシステム |
CN112671006A (zh) * | 2020-12-18 | 2021-04-16 | 国网江苏省电力有限公司经济技术研究院 | 一种评估海上风电场柔性直流输电系统谐振稳定性的方法 |
-
2023
- 2023-06-30 CN CN202310796506.7A patent/CN116822436B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104638678A (zh) * | 2015-03-16 | 2015-05-20 | 国网山西省电力公司电力科学研究院 | 一种多风机并网谐振分析方法 |
CN105160058A (zh) * | 2015-07-10 | 2015-12-16 | 国网四川省电力公司电力科学研究院 | 一种用于电网电磁暂态计算的网络边界等值方法 |
CN106021682A (zh) * | 2016-05-13 | 2016-10-12 | 清华大学 | 基于阻抗网络模型的次同步振荡分析方法和装置 |
CN107681682A (zh) * | 2017-10-25 | 2018-02-09 | 国家电网公司 | 一种基于ward等值的交直流系统等值方法 |
CN111597680A (zh) * | 2020-04-09 | 2020-08-28 | 华北电力大学 | 基于模态频率处电气阻尼的次同步谐振风险量化评估方法 |
WO2022233079A1 (zh) * | 2021-05-07 | 2022-11-10 | 南方电网科学研究院有限责任公司 | 适用于柔性直流的串联型无源阻抗适配器参数设计方法 |
CN114400667A (zh) * | 2021-11-15 | 2022-04-26 | 国网冀北电力有限公司智能配电网中心 | 交直流混合配电网安全风险评估方法 |
CN114301055A (zh) * | 2022-02-17 | 2022-04-08 | 河海大学 | 基于宽频量测的电力系统间谐波潮流获取方法与系统 |
CN116054187A (zh) * | 2022-04-27 | 2023-05-02 | 国网上海市电力公司 | 一种风电经柔直并网产生的次同步振荡抑制方法 |
CN114966203A (zh) * | 2022-06-16 | 2022-08-30 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 电力系统谐波振荡的检测方法、装置和计算机设备 |
CN115021256A (zh) * | 2022-08-08 | 2022-09-06 | 四川大学 | 一种大规模交直流输电系统电磁暂态模型自动生成方法 |
CN115549093A (zh) * | 2022-11-24 | 2022-12-30 | 中国电力科学研究院有限公司 | 一种新能源电力系统在线建模与振荡分析的方法及系统 |
CN115833124A (zh) * | 2022-12-28 | 2023-03-21 | 国网浙江省电力有限公司电力科学研究院 | 换流站直流近区谐振风险评估方法、装置、设备及介质 |
Non-Patent Citations (6)
Title |
---|
A Simple Model Predictive Power Control Strategy for Single-Phase PWM Converters With Modulation Function Optimization;Wensheng Song;IEEE Transactions on Powe Electrics;第31卷(第7期);第5279-5289页 * |
Chun Zhang,et al..Harmonic Impedance Model of UHVDC Transmission Lines Based on Distributed Parameters.2023 Panda Forum on Power and Energy.2023,第1939-1944页. * |
含直流模型的外网静态等值研究;曾庆涛;李岩;;电气开关(第02期);全文 * |
宁夏多直流外送系统SSO特性分析及次同步阻尼控制器设计;江桂芬;孙海顺;陈霞;朱鑫要;高峰;;电工技术学报(第S1期);全文 * |
弱电网下并网逆变器正交功率同步控制策略;王顺亮等;中国电机工程学报;第42卷(第23期);第8475-8485页 * |
抑制不同运行方式下故障电流的直流电网网架结构优化方法;王鹏宇等;现代电力;第1-8页 * |
Also Published As
Publication number | Publication date |
---|---|
CN116822436A (zh) | 2023-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109802406B (zh) | 一种分析柔性直流输电接入系统谐振稳定性的方法 | |
Sowa et al. | Impedance-based analysis of harmonic resonances in HVDC connected offshore wind power plants | |
Sainz et al. | Effect of wind turbine converter control on wind power plant harmonic response and resonances | |
Ajala et al. | Uncovering the kuramoto model from full-order models of grid-forming inverter-based power networks | |
CN112485509A (zh) | 基于非线性宽频模型的暂态过电压测量装置与方法 | |
Zhang et al. | Harmonic resonance suppression strategy of the front-end VIENNA rectifier in EV charging piles | |
Joseph et al. | Analysis of harmonic transfer through an MVDC link | |
Suzuki et al. | SS and SP topology analysis for capacitive power transfer with resonance coupling based on power factor consideration | |
CN116822436B (zh) | 一种直流输电送端交流系统的振荡风险灵敏度分析方法 | |
CN107942204B (zh) | 一种直流汇集电缆的综合保护方法及装置 | |
CN111276949B (zh) | 一种基于谐振电流幅值比较的光伏直流升压汇集系统线路保护方法 | |
CN114512966B (zh) | 基于s变换频域阻抗的直流配电网纵联保护方法及系统 | |
CN109635374A (zh) | 一种基于动态相量时域法的短路电流计算方法及装置 | |
CN114301055A (zh) | 基于宽频量测的电力系统间谐波潮流获取方法与系统 | |
CN111092453B (zh) | 一种多点接入光伏电源的电网关键线路辨识方法 | |
CN110808593B (zh) | 弱联系电网中多控制器引发电磁振荡阻尼分析方法及装置 | |
CN113852120A (zh) | 一种新能源发电最大并网容量的确定方法及系统 | |
Gao et al. | Improved extended kalman filter based dynamic equivalent method of DFIG wind farm cluster | |
CN113311287A (zh) | 一种应用于双端vsc控制的直流配电网的故障测距方法 | |
CN115085164B (zh) | 风场送出站交互式保护方法和系统 | |
CN112152138B (zh) | 一种高压并联电抗器容量的确定方法及系统 | |
Netto et al. | On the use of smart meter data to estimate the voltage magnitude on the primary side of distribution service transformers | |
CN114362181B (zh) | 一种适用于交直流混联电网的统一谐波状态估计方法 | |
Wang et al. | Resilience analysis and optimal control of grid‐connected converters under unbalanced voltage conditions | |
Alghamdi et al. | Analysis of harmonic propagations in albaha power network due to the implementation of an mvdc converter |
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 |