CN111239511A - 一种多馈入新能源并网系统广义短路比的在线测量方法 - Google Patents
一种多馈入新能源并网系统广义短路比的在线测量方法 Download PDFInfo
- Publication number
- CN111239511A CN111239511A CN202010005942.4A CN202010005942A CN111239511A CN 111239511 A CN111239511 A CN 111239511A CN 202010005942 A CN202010005942 A CN 202010005942A CN 111239511 A CN111239511 A CN 111239511A
- Authority
- CN
- China
- Prior art keywords
- matrix
- feed
- node
- nodes
- discrete time
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
Abstract
本发明公开了一种多馈入新能源并网系统广义短路比的在线测量方法。将系统节点分为馈入端口节点、网络内节点和无穷大电网等值节点三类,在上述节点中采集获取馈入端口节点的电流和电压的离散时间点数据,进而获得电流差值矩阵与电压差值矩阵;电流差值矩阵与电压差值矩阵通过重排的克罗降阶导纳矩阵建立关系模型,再根据最小二乘法矩阵迭代解法求解重排矩阵,利用重排矩阵获得广义短路比。本发明的广义短路比实时测量算法,能够实现广义短路比的实时在线测量,同时基于最小二乘法的矩阵迭代算法使得计算结果对于测量噪声干扰有着较强的抵御能力。
Description
技术领域
本发明涉及电气信息技术领域,具体为一种多馈入新能源并网系统广义短路比的在线测量方法。
背景技术
随着光伏、风能等可再生能源发电大量接入电力系统,国内外对可再生能源发电并网稳定性指标进行了多种研究。
广义短路比是其中一种刻画电力电子多馈入电力系统交流电网相对强度和稳定性的指标。根据广义短路比理论,风电场的稳定裕度可以仅通过关注电网强度及其相应边界来获得。因此,了解广义短路比有助于电力系统就风电场的稳定性进行管控。然而,由于电力系统运行和保护需求的时变性,电网结构在不断的变化,这使得广义短路比的在线监测变得复杂。
发明内容
为解决上述问题,本发明提出了一种多馈入新能源并网系统广义短路比的在线测量方法,能有效跟踪电力系统时变数据,实现广义短路比实时监控。
本发明的技术方案采用如下步骤:
1)将系统节点分为馈入端口节点、网络内节点和无穷大电网等值节点多类,在上述节点中采集获取馈入端口节点的电流和电压的离散时间点数据,进而获得电流差值矩阵与电压差值矩阵;
2)电流差值矩阵与电压差值矩阵通过重排的克罗降阶导纳矩阵建立关系模型,再根据最小二乘法矩阵迭代解法求解重排矩阵,利用重排矩阵获得广义短路比。
所述多馈入风电场并网系统节点分为馈入端口节点、网络内节点以及无穷大电网等值节点三类,具体为:无穷大电网等值节点为交流网络等值节点,在小干扰稳定分析领域将无穷大电网等值节点看作为接地节点,交流网络等值节点是指多馈入风电场并网系统中的风电场馈入网络经过等值后的接地节点,馈入端口节点为多馈入风电场并网系统的连接单台风机的节点,设有n台风机,则有节点1~n,网络内节点为多馈入风电场并网系统内除了无穷大电网等值节点以外的不与风机直接相连的节点,设有m个网络内节点,则有节点n+1~n+m。
所述2)中,利用馈入端口节点的电流和电压的离散时间点数据计算电流差值和电压差值,具体为某离散时刻数据值减去前一离散时刻数据值,每个离散时刻数据值是通过同步相量测量装置(phasor measurement units)实时采集获取,表示为:
其中,代表馈入端口节点i在第k个离散时间点与第k-1个离散时间点的电流差值;代表馈入端口节点i在第k个离散时间点与第k-1个离散时间点的电压差值,和分别表示馈入端口节点i在第k个离散时间点的电流和电压数值;
然后由某个时间点(如第k个离散时间点)各个馈入端口节点电流差值按照节点顺序组成列矩阵作为电流差值矩阵,由某个时间点(如第k个离散时间点)各个馈入端口节点电压差值按照节点顺序组成列矩阵作为电压差值矩阵。
所述2)具体为:
式中,表示实数,Γi(N)表示第i个馈入端口节点前N个离散时间点的电流差值矩阵,u(N)表示所有馈入端口节点前N个离散时间点的电压差值矩阵;表示所有馈入端口节点前N个离散时间点的电流差值矩阵,表示u(N)与变换矩阵Ti的乘积按照馈入端口节点1~n顺序构成的列矩阵,i表示馈入端口节点的序数,N表示离散时间点的总数,n表示馈入端口节点的总数;
2.2)然后采用以下公式作为最小二乘法矩阵迭代算法,对上述关系模型进行求解:
本发明的有益效果是:
本发明能有效抵御测量过程中的噪声干扰,且本发明能根据电力系统实时数据,实现广义短路比在线测量,为电力系统就风电场稳定性判断提供实时依据。
附图说明
图1为本发明的风电场典型节点分布图与克罗降阶等效节点图。
图2为本发明注入噪声强度图。
图3为本发明在不同测量噪声下迭代数值图。
具体实施方式
下面结合附图对本发明作进一步详细说明。
按照本发明方法完整实施的实施例及其实施过程如下:
在此模型中,节点1~n代表n个风机馈入端口节点;节点n+1~n+m代表风电场内不与风机直接相连的节点,称之为网络内节点;交流电网被简化为一个无穷大电网等值节点(图1中Infinite Bus)。
根据电网络理论,节点i和节点j之间阻抗可以表示为:Zij=Rij+jXij,相应的导纳为:Yij=1/Zij。则节点导纳矩阵可表示为:
Y∈(n+m)×(n+m)
经过克罗降阶(Kron reduction)后的导纳矩阵表示为:
式中Y1∈n×n,Y2∈n×m,Y3∈m×n,Y4∈m×m为节点导纳矩阵Y的子矩阵:
由某个时间点(如第k个离散时间点)各个馈入端口节点电流差值按照节点1~n顺序组成列矩阵作为电流差值矩阵ΔI(k),由某个时间点(如第k个离散时间点)各个馈入端口节点电压差值按照节点1~n顺序组成列矩阵作为电压差值矩阵ΔU(k)。电流差值矩阵与电压差值矩阵即可通过上述克罗降阶导纳矩阵建立起数学关系式,具体表达式如下所示:
此式中ΔUT (k)为ΔU(k)的转置,为矩阵第i行的转置。进一步,将按照1~N个离散数据点的顺序,扩展为列矩阵Γi(N);将ΔUT (k)按照1~N个离散数据点的顺序,扩展为列矩阵u(N),则此表达式变更为:
即
式中符号+代表矩阵的广义逆。
为验证迭代算法的有效性,考虑一个具体的含有9个风力发电机的风电场模型,如图1所示,图中编号1~9代表馈入端口节点;编号10~11代表网络内节点;Infinite Bus代表无穷大电网等值节点;Kron Reduction代表克罗降阶过程。
在馈入端口节点注入不同大小的噪声,如图2所示,从左到右顺序,第一幅图表示注入噪声强度为0,第二幅图表示注入噪声强度较小,第三幅图表示注入噪声强度较大,第四幅图表示注入噪声强度最大。
图3表示在不同大小的噪声干扰下迭代数值的变化,从图中可看出,在四种不同强度噪声干扰下,广义短路比(gSCR)最终都快速收敛为实际值,证明了迭代算法的有效性。
由此实施可见,本发明的广义短路比实时测量算法,能够实现广义短路比的在线测量。本发明的主要创新点在于:1)通过测量馈入端口节点电压和电流来间接计算广义短路比的思想;2)引入了重排导纳矩阵的概念,提高了计算精度;3)提出了一种能有效抵御噪声干扰的迭代算法。
Claims (4)
1.一种多馈入新能源并网系统广义短路比的在线测量方法,其特征在于包括以下步骤:
1)将系统节点分为馈入端口节点、网络内节点和无穷大电网等值节点多类,在上述节点中采集获取馈入端口节点的电流和电压的离散时间点数据,进而获得电流差值矩阵与电压差值矩阵;
2)电流差值矩阵与电压差值矩阵通过重排的克罗降阶导纳矩阵建立关系模型,再根据最小二乘法矩阵迭代解法求解重排矩阵,利用重排矩阵获得广义短路比。
2.根据权利要求1所述的一种多馈入新能源并网系统广义短路比的在线测量方法,其特征在于:所述多馈入风电场并网系统节点分为馈入端口节点、网络内节点以及无穷大电网等值节点三类,具体为:无穷大电网等值节点为交流网络等值节点,馈入端口节点为多馈入风电场并网系统的连接单台风机的节点,网络内节点为多馈入风电场并网系统内除了无穷大电网等值节点以外的不与风机直接相连的节点。
3.根据权利要求1所述的一种多馈入新能源并网系统广义短路比的在线测量方法,其特征在于:所述2)中,利用馈入端口节点的电流和电压的离散时间点数据计算电流差值和电压差值,具体为某离散时刻数据值减去前一离散时刻数据值,每个离散时刻数据值是通过同步相量测量装置实时采集获取,表示为:
其中,代表馈入端口节点i在第k个离散时间点与第k-1个离散时间点的电流差值;代表馈入端口节点i在第k个离散时间点与第k-1个离散时间点的电压差值,和分别表示馈入端口节点i在第k个离散时间点的电流和电压数值;
然后由一个时间点各个馈入端口节点电流差值按照节点顺序组成列矩阵作为电流差值矩阵,由一个时间点各个馈入端口节点电压差值按照节点顺序组成列矩阵作为电压差值矩阵。
4.根据权利要求1所述的一种多馈入新能源并网系统广义短路比的在线测量方法,其特征在于:所述2)具体为:
式中,表示实数,Γi(N)表示第i个馈入端口节点前N个离散时间点的电流差值矩阵,u(N)表示所有馈入端口节点前N个离散时间点的电压差值矩阵;表示所有馈入端口节点前N个离散时间点的电流差值矩阵,表示u(N)与变换矩阵Ti的乘积按照馈入端口节点1~n顺序构成的列矩阵,i表示馈入端口节点的序数,N表示离散时间点的总数,n表示馈入端口节点的总数;
2.2)然后采用以下公式作为最小二乘法矩阵迭代算法,对上述关系模型进行求解:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010005942.4A CN111239511B (zh) | 2020-01-03 | 2020-01-03 | 一种多馈入新能源并网系统广义短路比的在线测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010005942.4A CN111239511B (zh) | 2020-01-03 | 2020-01-03 | 一种多馈入新能源并网系统广义短路比的在线测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111239511A true CN111239511A (zh) | 2020-06-05 |
CN111239511B CN111239511B (zh) | 2021-01-26 |
Family
ID=70872383
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010005942.4A Active CN111239511B (zh) | 2020-01-03 | 2020-01-03 | 一种多馈入新能源并网系统广义短路比的在线测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111239511B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106410848A (zh) * | 2016-11-08 | 2017-02-15 | 浙江大学 | 一种电力电子多馈入电力系统小干扰稳定性评估方法 |
CN106655234A (zh) * | 2016-10-11 | 2017-05-10 | 国网浙江省电力公司经济技术研究院 | 一种线路阻抗和联络线功率对广义短路比影响的分析方法 |
CN109004652A (zh) * | 2018-07-11 | 2018-12-14 | 国网浙江省电力有限公司经济技术研究院 | 一种含异步电机的直流多馈入系统受端电网强度获得方法 |
CN109245146A (zh) * | 2018-09-14 | 2019-01-18 | 上海电力学院 | 基于多馈入广义运行短路比的电力系统稳定性判断方法 |
-
2020
- 2020-01-03 CN CN202010005942.4A patent/CN111239511B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106655234A (zh) * | 2016-10-11 | 2017-05-10 | 国网浙江省电力公司经济技术研究院 | 一种线路阻抗和联络线功率对广义短路比影响的分析方法 |
CN106410848A (zh) * | 2016-11-08 | 2017-02-15 | 浙江大学 | 一种电力电子多馈入电力系统小干扰稳定性评估方法 |
CN109004652A (zh) * | 2018-07-11 | 2018-12-14 | 国网浙江省电力有限公司经济技术研究院 | 一种含异步电机的直流多馈入系统受端电网强度获得方法 |
CN109245146A (zh) * | 2018-09-14 | 2019-01-18 | 上海电力学院 | 基于多馈入广义运行短路比的电力系统稳定性判断方法 |
Non-Patent Citations (3)
Title |
---|
FENG ZHANG 等: "Geberalized short ratio for multi-infeed LCC-HVDC systems", 《2017 IEEE POWER & ENERGY SOCIETY GENERAL MEETING》 * |
王冠中: "基于广义短路比的电力电子多馈入系统小干扰概率稳定评估", 《电力系统自动化》 * |
白加林: "具有多馈入直流输电的交直流电力系统电磁暂态仿真研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111239511B (zh) | 2021-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ju et al. | An extension of FBS three-phase power flow for handling PV nodes in active distribution networks | |
JP6025575B2 (ja) | 三相電力潮流解析のためのシステムおよび方法 | |
CN108400595B (zh) | 一种考虑新能源出力相关性的电压暂降随机预估方法 | |
CN107947192B (zh) | 一种下垂控制型孤岛微电网的无功优化配置方法 | |
CN105552938B (zh) | 一种三相不对称配电网络电压暂降评估方法 | |
CN106549396A (zh) | 一种配电网多目标概率无功优化方法 | |
CN109301870B (zh) | 一种电力电子多馈入电力系统容量优化方法 | |
CN107834593B (zh) | 一种下垂控制型孤岛微电网静态电压稳定概率评估方法 | |
CN108599239A (zh) | 一种下垂控制型孤岛微电网电压质量概率评估方法 | |
Montoya et al. | Sine-cosine algorithm for parameters’ estimation in solar cells using datasheet information | |
CN110941909A (zh) | 光伏电站和风电场参数主导性辨识计算方法及系统 | |
Wang et al. | Applying probabilistic collocation method to power flow analysis in networks with wind farms | |
CN103956735A (zh) | 一种分布式发电系统的谐波潮流分析方法 | |
CN111327048A (zh) | 一种含三端snop的配电网鲁棒运行优化方法 | |
Benato | A basic AC power flow based on the bus admittance matrix incorporating loads and generators including slack bus | |
Pan et al. | Probabilistic voltage quality evaluation of islanded droop-regulated microgrid based on non-intrusive low rank approximation method | |
CN111239511B (zh) | 一种多馈入新能源并网系统广义短路比的在线测量方法 | |
Shahraeini et al. | Behavior analysis of random power graphs for optimal PMU placement in smart grids | |
CN112234598A (zh) | 一种电磁暂态仿真初始化方法 | |
Orozco-Gutierrez et al. | A method for the fast estimation of the maximum power points in mismatched PV strings | |
CN116304487A (zh) | 一种电力系统潮流反函数计算方法 | |
Kamarposhti et al. | Locating and sizing of distributed generation sources and parallel capacitors using multiple objective particle swarm optimization algorithm | |
Quan et al. | Spatial correlation modeling for optimal power flow with wind power: Feasibility in application of superconductivity | |
Kumar et al. | A simplified approach for load flow analysis of radial distribution network | |
CN104699950B (zh) | 一种包含随机功率单元相关性问题的随机潮流计算方法 |
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 |