CN105811408B - 面向mmc换流站的谐波能量传递函数建模方法 - Google Patents

面向mmc换流站的谐波能量传递函数建模方法 Download PDF

Info

Publication number
CN105811408B
CN105811408B CN201610253389.XA CN201610253389A CN105811408B CN 105811408 B CN105811408 B CN 105811408B CN 201610253389 A CN201610253389 A CN 201610253389A CN 105811408 B CN105811408 B CN 105811408B
Authority
CN
China
Prior art keywords
energy
current
harmonic
curve
driving source
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.)
Expired - Fee Related
Application number
CN201610253389.XA
Other languages
English (en)
Other versions
CN105811408A (zh
Inventor
李国庆
陈继开
李炳辉
王振浩
李江
王鹤
李扬
辛业春
李卫国
范维
韩子娇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Electric Power Research Institute of State Grid Liaoning Electric Power Co Ltd
Northeast Electric Power University
Original Assignee
Northeast Dianli University
Electric Power Research Institute of State Grid Liaoning Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northeast Dianli University, Electric Power Research Institute of State Grid Liaoning Electric Power Co Ltd filed Critical Northeast Dianli University
Priority to CN201610253389.XA priority Critical patent/CN105811408B/zh
Publication of CN105811408A publication Critical patent/CN105811408A/zh
Application granted granted Critical
Publication of CN105811408B publication Critical patent/CN105811408B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/36Arrangements for transfer of electric power between ac networks via a high-tension dc link
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/60Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Power Engineering (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

一种面向MMC换流站的谐波能量传递函数建模方法,属于电力传输技术领域。本发明的目的是针对目前柔性直流输电***利用谐波注入法分析波形频谱无法得到MMC换流站谐波能量传递函数的缺点,建立了交直流***混联***MMC换流站谐波能量传递过程的面向MMC换流站的谐波能量传递函数建模方法。本发明步骤是:在多电平双端柔性直流输电***中添加谐波信号,时域仿真和提取数据,曲线积分,改变x的数值,绘制能量曲线,计算谐波能量传递函数,曲线修正。本发明填补了交直流混合***MMC换流站的谐波能量传递函数建模的空白,为包含MMC换流站的交直流混合***电能质量研究提供了有效的分析手段。

Description

面向MMC换流站的谐波能量传递函数建模方法
技术领域
本发明属于电力传输技术领域。
背景技术
在高压柔性直流输电***中,交流***谐波能够通过换流变压器经MMC换流器向直流侧传递。当MMC换流器的固有频率与谐波频率相同或相近时,***极易发生谐振,导致直流电压、电流的波动,严重时会引起***震荡,导致电能质量恶化。这些谐波如果不能加以有效抑制,则会对与换流站相连直流线路以及直流线路沿线的通讯设备构成威胁,对柔性直流输电***的电能质量造成严重的危害。
谐波在交直流***之间的能量传递过程受多种因素的影响。随着线路参数和终端负载的不同,传递系数的变化比较复杂,没有固定的变化规律。交流侧和直流侧谐波之间的关系的计算与交流侧的正序和负序分量均有关,在交流***通过直流直流线路传递的谐波的幅值还受平波电抗器和滤波器的影响,也受直流电路感抗的影响。谐波传递的过程中,往往产生间谐波,与传统的谐波相比,间谐波的抑制和治理更加困难。
由于MMC换流器内部含有大量非线性元件(如阀组与电容单体组成的子模块),分析描述MMC换流器谐波能量传递特性时比较复杂。目前,关于传统直流输电***谐波传递的研究大多基于谐波注入法,即在***的交流侧施加不同频率的交流电压源、电流源激励,对直流侧电压、电流进行傅里叶分解,分析得到直流侧响应出的谐波频次和含量。此类研究多集中于从响应信号的频谱能量分布的角度对传统换流站谐波传递特性进行分析,不能获得柔性直流输电***的谐波能量传递函数,无法直观表征高压柔性直流输电***的谐波能量传递的规律。
发明内容
本发明的目的是针对目前柔性直流输电***利用谐波注入法分析波形频谱无法得到MMC换流站谐波能量传递函数的缺点,建立了交直流***混联***MMC换流站谐波能量传递过程的面向MMC换流站的谐波能量传递函数建模方法。
本发明步骤是:
①在多电平双端柔性直流输电***中添加谐波信号:在***交流阀侧添加一系列混合频率的谐波电流激励源,激励源的数学表达式如下:
(1)
式中为比例系数, 为激励源的幅值,f为工作频率50Hz,n为谐波次数,为施加激励源的初相角;
②时域仿真和提取数据:对所搭建的多电平双端柔性直流输电***模型进行时域仿真,实时监测所要研究的电气量;施加激励源后,交流谐波电流通过整流器传递到直流侧,直流电流会随着式(1)激励源的改变而变化,不同含量的激励源激励会响应出不同数值的直流电流值;提取***直流电流稳态过程的数据;
③曲线积分:对***稳定后的直流电流波形进行积分处理,将②提取的直流电流数据形成矩阵,利用trapz函数对提取的直流电流稳态过程的数据进行积分,用该积分表示为在当前混合频率谐波电流激励源激励下,***直流电流中所含有的能量波动
④改变x的数值:重复进行上述步骤②至③,计算***直流电流能量波动;为避免谐波含量过高导致***不稳定,谐波电流源的幅值比例系数x的设定应使***保持在稳定范围内,即;等比例的增大或减小混合频率谐波电流激励源的幅值,计算***直流电流稳定后所含有的能量波动
⑤绘制能量曲线,计算谐波能量传递函数:利用高斯插值法对能量波动数据进行曲线拟合,得到能量函数曲线;针对不同幅值的混合谐波电流激励源激励下直流电流中所含有的能量,绘制激励能量-响应能量曲线;曲线的X轴为混合频率激励源所含有的能量,Y轴为直流电流波动能量;利用高斯插值公式计算谐波能量的传递函数,式中,a,b,c为待定系数,M施加谐波电流激励源的幅值;
⑥曲线修正:根据步骤⑤绘制出的能量曲线,对曲线数据突变的部分采取密集采样,即在该部分继续等比例改变所施加混合频率激励源的幅值,重复步骤②至⑤,计算此部分激励对应的直流电流所含有的能量,修正能量函数曲线;鉴于能量函数的波动较大,采用分段函数形式对能量曲线进行二次修正,得到精确度较高的能量函数曲线和能量函数
本发明在仿真***中搭建包含完整控制***的多电平双端柔性直流输电***,在该混联***交流侧添加不同频率的交流谐波电流源,提取***直流侧电流数据。对提取出的数据取绝对值后积分,得到***直流电流中含有的谐波能量,利用高斯插值法对能量曲线进行拟合。采用密集采样法对能量函数曲线进行一次修正,利用分段函数法对能量函数曲线进行二次修正以提高拟合的精度。进而得到混联交直流***的谐波能量传递函数。该方法填补了交直流混合***MMC换流站的谐波能量传递函数建模的空白,为包含MMC换流站的交直流混合***电能质量研究提供了有效的分析手段。
附图说明
图1是MMC换流器的拓扑结构图;
图2是MMC换流站的谐波能量传递图形;
图3是针对谐波能量传递图形采用密集采样法的一次修正图形;
图4(1)针对一次修正图形采用分段函数法二次修正图形的第一段(0-0.2)的修正图;
图4(2)针对一次修正图形采用分段函数法二次修正图形的第二段(0.2-0.3)的修正图;
图4(3)针对一次修正图形采用分段函数法二次修正图形的第三段(0.3-0.6)的修正图;
图4(4)针对一次修正图形采用分段函数法二次修正图形的第四段(0.6-0.7)的修正图;
图4(5)针对一次修正图形采用分段函数法二次修正图形的第五段(0.7-1.0)的修正图。
具体实施方式
本发明步骤是:
在仿真软件中搭建包含完整控制***的多电平双端柔性直流输电***。***的等效模型如图1所示。
①在多电平双端柔性直流输电***中添加谐波信号:在***交流阀侧添加一系列混合频率的谐波电流激励源,激励源的幅值比例按国家电能质量规定中的要求设定。激励源的数学表达式如下:
(1)
式中为比例系数, 为激励源的幅值,f为工作频率50Hz,n为谐波次数,为施加激励源的初相角;
②时域仿真和提取数据:对所搭建的多电平双端柔性直流输电***模型进行时域仿真,实时监测所要研究的电气量;施加激励源后,交流谐波电流通过整流器传递到直流侧,直流电流会随着式(1)激励源的改变而变化,不同含量的激励源激励会响应出不同数值的直流电流值;提取***直流电流稳态过程的数据;
③曲线积分:对***稳定后的直流电流波形进行积分处理,将②提取的直流电流数据形成矩阵,利用trapz函数对提取的直流电流稳态过程的数据进行积分,用该积分表示为在当前混合频率谐波电流激励源激励下,***直流电流中所含有的能量波动
④改变x的数值:重复进行上述步骤②至③,计算***直流电流能量波动;为避免谐波含量过高导致***不稳定,谐波电流源的幅值比例系数x的设定应使***保持在稳定范围内,即;等比例的增大或减小混合频率谐波电流激励源的幅值,计算***直流电流稳定后所含有的能量波动
⑤绘制能量曲线,计算谐波能量传递函数:利用高斯插值法对能量波动数据进行曲线拟合,得到能量函数曲线;针对不同幅值的混合谐波电流激励源激励下直流电流中所含有的能量,绘制激励能量-响应能量曲线;曲线的X轴为混合频率激励源所含有的能量,Y轴为直流电流波动能量;利用高斯插值公式计算谐波能量的传递函数,式中,a,b,c为待定系数,M施加谐波电流激励源的幅值;
⑥曲线修正:根据步骤⑤绘制出的能量曲线,对曲线数据突变的部分采取密集采样,即在该部分继续等比例改变所施加混合频率激励源的幅值,重复步骤②至⑤,计算此部分激励对应的直流电流所含有的能量,修正能量函数曲线;鉴于能量函数的波动较大,采用分段函数形式对能量曲线进行二次修正,得到精确度较高的能量函数曲线和能量函数
实施例1:
参照图1谐波电流源施加的位置,在该位置处施加不同频次的交流电流源,其中n为谐波次数,其取值为3,5,7,9,11,13,15,17,19;为施加激励源的初相角0°;为激励源幅值分别,取值为9.6,9.6,6.8,3.2,4.3,3.7,1.9,2.8,0.6。按照上述步骤②至⑤,等比例改变谐波电流源比例系数,幅值由0.1,0.2增大至。提取实验所需的数据,并对数据取绝对值后积分处理得到直流电流中谐波的能量波动
表1 不同谐波能量注入下直流电流中谐波的能量
根据高斯插值法拟合上述数据,得到混联交直流***的谐波能量传递函数曲线,如图2所示。利用密集采样法,继续改变谐波电流源的比例系数,幅值由增大至,对得到的函数曲线进行一次修正如图3所示。由二次修正得到的图形看出,能量函数曲线变化不规则,曲线由五段组成,采用分段函数法求每段精确地函数表达式,得到二次修正后的函数曲线如图4(1)、图4(2)、图4(3)、图4(4)图4(5)所示。最终得到混联交直流精确的谐波能量传递函数,,,,
由图2、图3、图4可以清晰看到,在***的交流侧有谐波电流注入的情况下,谐波经换流变压器由MMC换流器传递到直流侧的能量,受交流侧谐波注入的能量的控制。通过高斯插值法拟合曲线得到混联交直流***的谐波能量函数曲线,采用分段函数法求得混联交直流***的谐波能量传递函数。本发明建立了MMC换流站的谐波能量传递函数建模方法。

Claims (1)

1.一种面向MMC换流站的谐波能量传递函数建模方法,其特征在于:其步骤是:
①在多电平双端柔性直流输电***中添加谐波信号:在***交流阀侧添加一系列混合频率的谐波电流激励源,激励源的数学表达式如下:
(1)
式中为比例系数, 为激励源的幅值,f为工作频率50Hz,n为谐波次数,为施加激励源的初相角;
②时域仿真和提取数据:对所搭建的多电平双端柔性直流输电***模型进行时域仿真,实时监测所要研究的电气量;施加激励源后,交流谐波电流通过整流器传递到直流侧,直流电流会随着式(1)激励源的改变而变化,不同含量的激励源激励会响应出不同数值的直流电流值;提取***直流电流稳态过程的数据;
③曲线积分:对***稳定后的直流电流波形进行积分处理,将②提取的直流电流数据形成矩阵,利用trapz函数对提取的直流电流稳态过程的数据进行积分,用该积分表示为在当前混合频率谐波电流激励源激励下,***直流电流中所含有的能量波动
④改变x的数值:重复进行上述步骤②至③,计算***直流电流能量波动;为避免谐波含量过高导致***不稳定,谐波电流源的幅值比例系数x的设定应使***保持在稳定范围内,即;等比例的增大或减小混合频率谐波电流激励源的幅值,计算***直流电流稳定后所含有的能量波动
⑤绘制能量曲线,计算谐波能量传递函数:利用高斯插值法对能量波动数据进行曲线拟合,得到能量函数曲线;针对不同幅值的混合谐波电流激励源激励下直流电流中所含有的能量,绘制激励能量-响应能量曲线;曲线的X轴为混合频率激励源所含有的能量,Y轴为直流电流波动能量;利用高斯插值公式计算谐波能量的传递函数,式中,a,b,c为待定系数,M施加谐波电流激励源的幅值;
⑥曲线修正:根据步骤⑤绘制出的能量曲线,对曲线数据突变的部分采取密集采样,即在该部分继续等比例改变所施加混合频率激励源的幅值,重复步骤②至⑤,计算此部分激励对应的直流电流所含有的能量,修正能量函数曲线;鉴于能量函数的波动较大,采用分段函数形式对能量曲线进行二次修正,得到精确度较高的能量函数曲线和能量函数
CN201610253389.XA 2016-04-22 2016-04-22 面向mmc换流站的谐波能量传递函数建模方法 Expired - Fee Related CN105811408B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610253389.XA CN105811408B (zh) 2016-04-22 2016-04-22 面向mmc换流站的谐波能量传递函数建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610253389.XA CN105811408B (zh) 2016-04-22 2016-04-22 面向mmc换流站的谐波能量传递函数建模方法

Publications (2)

Publication Number Publication Date
CN105811408A CN105811408A (zh) 2016-07-27
CN105811408B true CN105811408B (zh) 2018-04-03

Family

ID=56457301

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610253389.XA Expired - Fee Related CN105811408B (zh) 2016-04-22 2016-04-22 面向mmc换流站的谐波能量传递函数建模方法

Country Status (1)

Country Link
CN (1) CN105811408B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106569039B (zh) * 2016-11-15 2019-10-11 东北电力大学 Mmc换流站谐波单向传递特性分析方法
CN107907737A (zh) * 2017-10-16 2018-04-13 国网浙江省电力公司电力科学研究院 一种基于指数时域函数的电度计量方法及***
CN108256252A (zh) * 2018-01-29 2018-07-06 中国电力科学研究院有限公司 一种基于滤波器的电磁暂态仿真方法及***
CN111276994B (zh) * 2020-03-02 2021-12-28 中国电建集团华东勘测设计研究院有限公司 海上换流站整流侧变流器对谐波的响应特性分析方法
CN112018768B (zh) * 2020-07-23 2022-03-08 南方电网科学研究院有限责任公司 多端柔性直流输电***谐波传递特性的分析方法及装置
CN113541142A (zh) * 2021-06-10 2021-10-22 中国电建集团华东勘测设计研究院有限公司 一种陆上换流站对谐波的响应特性分析及参数优化方法
CN113394782B (zh) * 2021-06-23 2022-11-29 东北电力大学 基于apf直流电压信息反演的工业园区谐波监测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103745065A (zh) * 2014-01-17 2014-04-23 国家电网公司 一种确定电气化铁路多谐波源叠加后电流的方法
CN104237635A (zh) * 2014-09-17 2014-12-24 国家电网公司 电力***的谐波定量测量方法
CN104393614A (zh) * 2014-03-21 2015-03-04 南车株洲电力机车研究所有限公司 一种模块化多电平双流制牵引供电及净化***
CN104407214A (zh) * 2014-11-27 2015-03-11 云南电网公司电力科学研究院 一种谐波源辨识方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103745065A (zh) * 2014-01-17 2014-04-23 国家电网公司 一种确定电气化铁路多谐波源叠加后电流的方法
CN104393614A (zh) * 2014-03-21 2015-03-04 南车株洲电力机车研究所有限公司 一种模块化多电平双流制牵引供电及净化***
CN104237635A (zh) * 2014-09-17 2014-12-24 国家电网公司 电力***的谐波定量测量方法
CN104407214A (zh) * 2014-11-27 2015-03-11 云南电网公司电力科学研究院 一种谐波源辨识方法

Also Published As

Publication number Publication date
CN105811408A (zh) 2016-07-27

Similar Documents

Publication Publication Date Title
CN105811408B (zh) 面向mmc换流站的谐波能量传递函数建模方法
CN104638678B (zh) 一种多风机并网谐振分析方法
CN104316894B (zh) 一种电能表实际运行环境的仿真校验方法
CN109217371B (zh) 考虑锁相环影响的电压源型变流器并网***稳定性分析方法、装置及***
CN108880300B (zh) 一种基于双闭环控制的双馈风机整流器阻抗计算方法
CN110943456A (zh) Lcc-hvdc***小信号谐波状态空间建模方法及装置
CN109038604B (zh) 一种statcom接入受端弱交流***的分析方法
CN103293963B (zh) 有源电力滤波器自适应模糊反演跟踪控制方法
CN105140927B (zh) 含谐波的配电网无功补偿电容器容量的计算方法
CN106655195B (zh) 一种主动式配电网高频谐波潮流的计算方法
CN106484962B (zh) 一种谐振式无线输电***稳态特性的符号分析方法
CN111654052A (zh) 基于动态相量法的柔直换流器建模装置及方法
CN107069723A (zh) 有源电力滤波器分数阶PIλ及自抗扰混合控制方法
CN103928930B (zh) 减少电动汽车充放电过程中电网谐波的方法
CN104111373B (zh) 一种用于智能变电站数字化电能的计量方法
CN101662148B (zh) 一种直角坐标牛顿法潮流计算的电压初值设置方法
CN110137968B (zh) 一种包含vsc的谐波潮流计算方法
CN103904652B (zh) 一种能抗脉冲噪声干扰的供电网谐波抑制方法和***
CN103151780B (zh) 一种三相三线制sapf的定时变环宽电流控制的方法
CN110502804B (zh) 一种电动汽车充电过程时变谐波电流评估方法及***
CN112255457A (zh) 适用于自动准同期装置的相角差测量方法
CN110350533B (zh) 计及新能源发电控制策略的改进仿射谐波潮流计算方法
CN106372369B (zh) 用于逆变器输出直流分量、谐波分量检测的瞬时计算方法
CN113676037A (zh) 一种新型无源功率因数校正装置
CN105277780A (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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Li Guoqing

Inventor after: Fan Wei

Inventor after: Han Zijiao

Inventor after: Chen Jikai

Inventor after: Li Binghui

Inventor after: Wang Zhenhao

Inventor after: Li Jiang

Inventor after: Wang He

Inventor after: Li Yang

Inventor after: Xin Yechun

Inventor after: Li Weiguo

Inventor before: Li Guoqing

Inventor before: Chen Jikai

Inventor before: Li Binghui

Inventor before: Wang Zhenhao

Inventor before: Li Jiang

Inventor before: Wang He

Inventor before: Li Yang

Inventor before: Xin Yechun

Inventor before: Li Weiguo

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20171205

Address after: 132012 Changchun Road, Jilin, China, No. 169, No.

Applicant after: NORTHEAST DIANLI University

Applicant after: STATE GRID LIAONING ELECTRIC POWER Research Institute

Address before: 132012 Changchun Road, Jilin, China, No. 169, No.

Applicant before: Northeast Dianli University

GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180403