CN110518615B - 基于二阶电路响应的交直流系统换相电压预测方法 - Google Patents
基于二阶电路响应的交直流系统换相电压预测方法 Download PDFInfo
- Publication number
- CN110518615B CN110518615B CN201910592330.7A CN201910592330A CN110518615B CN 110518615 B CN110518615 B CN 110518615B CN 201910592330 A CN201910592330 A CN 201910592330A CN 110518615 B CN110518615 B CN 110518615B
- Authority
- CN
- China
- Prior art keywords
- voltage
- component
- direct current
- commutation
- alternating 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.)
- Active
Links
Images
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
- H02J3/36—Arrangements for transfer of electric power between ac networks via a high-tension dc link
-
- 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
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Inverter Devices (AREA)
Abstract
本发明针对交直流互联系统提出了一种基于二阶电路响应的交直流系统换相电压预测方法。该方法首先基于二阶电路的响应形式建立交直流系统暂态响应的预测模型;然后根据对称分量法对换流母线电压三相分量进行分解,得到换流母线电压三序分量的傅里叶展开系数;最后基于最小二乘法原理对响应预测模型中各参数进行拟合,计算换相电压的预测值。该方法能在系统故障后快速预测换流母线各相电压,有助于及时制定和采取紧急控制措施,对维持电力系统稳定运行具有重要意义。
Description
技术领域
本发明涉及电力系统态势感知及趋势预测技术领域,特别是一种基于二阶电路响应的交直流系统换相电压预测方法。
背景技术
我国地域辽阔但能源分布很不均衡。能源丰富的西部地区与东部负荷中心距离较远,需要电能的远距离传输。高压直流输电系统因其可控性强、输电容量大、输电损耗小等优点被广泛用于远距离输电。
换相失败一直以来都是高压直流输电系统最常见的故障之一。换相失败会导致直流传输功率迅速下降,引起交流系统电压、电流的波动。严重情况下会引起连续、相继换相失败,造成直流闭锁甚至大规模停电事故。
目前工程上对于换相失败的判据主要是根据换流母线电压的幅值来制定的。由于电力系统存在电感、电容等储能元件,因此交直流系统在暂态过渡期间可以等效为二阶RLC电路,基于此可对换流母线电压进行预测,有助于及时制定并采取预防控制措施,对提升电力系统稳定性具有重要的意义。
发明内容
本发明的目的是为了解决上述问题,提供一种基于二阶电路响应的交直流系统换相电压预测方法。
该方法首先基于二阶电路的响应形式建立交直流系统暂态响应的预测模型,然后根据对称分量法对换流母线电压三相分量进行分解,得到换流母线电压三序分量的傅里叶展开系数;最后基于最小二乘法原理对响应预测模型中各参数进行拟合,计算换相电压的预测值。该方法能在系统故障后快速预测换流母线各相电压,有助于及时制定和采取紧急控制措施,对维持电力系统稳定运行具有重要意义。
为实现上述发明目的,本发明所采用的技术方案如下:一种基于二阶电路响应的交直流系统换相电压预测方法,包括如下步骤:
(1)基于二阶电路响应形式,建立交直流系统暂态响应的预测模型;
(2)收集换流母线三相电压数据;
(3)计算换流母线电压三序分量的傅里叶展开系数;
(4)基于最小二乘法原理,对响应预测模型中的参数进行拟合;
(5)预测换相电压三相分量;
进一步的,所述步骤(1)中交直流系统暂态响应预测模型的建立包括如下步骤:
(11)根据电力系统参数及储能元件的特点,得到暂态过程中的时域响应形式:
式中,f(t)为所求的响应,p1,2为暂态过渡期间交直流系统等值RLC电路的特征根,a(t)为f(t)的稳态分量,包括直流及各倍频交流分量;为f(t)的暂态分量;τ、ωd、K为与系统参数及初值条件相关的常数;j为虚数单位;t表示时间,为f(t)的暂态分量;
(12)根据下式计算f(t)稳态分量及暂态分量的傅里叶展开系数:
式中,ω0=100π(rad/s)为系统基频角速度,T=2π/ω0表示采样周期;An为稳态分量n倍频分量的系数,Xn为暂态分量n倍频分量的系数;特殊的,n=0表示直流分量,n=1表示基频分量,λ表示积分计算当中的积分元。
式中,λ表示积分计算当中的积分元,μ0、μ1、μ2、ν0、ν1、ν2为计算的中间变量,具体可根据下式计算得到:
进一步的,所述步骤(2)中的收集换流母线三相电压数据的具体实现方法为:以故障发生或恢复时刻为零时刻,以采样频率fs采集[0tβ]时间段内换流母线的三相电压ea、eb、ec;其中tβ表示换相开始时刻。
进一步的,所述步骤(3)中计算换流母线电压三序分量的傅里叶展开系数的具体实现步骤为:
(31)根据下式计算所需采样的波形总数N及一个采样周波内的采样点数M:
式中:floor表示向下取整函数,T=2π/ω0表示采样周期,△T表示两次采样波形的时间间隔,fs表示采样频率,tβ表示换相开始时刻。
式中,M表示一个采样周波内的采样点数,n表示所求傅里叶展开系数对应的频次,e表示相电压,下标i取a、b、c分别表示a、b、c三相,△T表示两次采样波形的时间间隔,fs表示采样频率,ε表示第ε个采样点。
进一步的,所述步骤(4)中响应预测模型中的参数的拟合,具体实现方法为:
式中:分别表示ψn′、ψn″的拟合参数,ψn′、ψn″分别表示实部及虚部的变化函数,real、imag分别表示取实部和取虚部,表示a相的采样数据的傅里叶展开系数,括号内ξ取1,2,0分别表示a相的正、负、零序,k表示第k次采样数据,fit是基于最小二乘法原理的拟合函数,其形式为表示对Y=fun(θ,X)进行拟合,并返回拟合参数集
进一步的,所述步骤(5)中换相电压预测的具体实现步骤为:
有益效果:
本发明提出了一种基于二阶电路响应的交直流系统换相电压预测方法。与现有方法相比,该方法考虑了电力系统中电感及电容等储能元件的作用,建立了交直流系统暂态响应预测模型;然后根据对称分量法对换流母线电压三相分量进行分解,得到换流母线电压三序分量的傅里叶展开系数;最后基于最小二乘法原理对响应预测模型中各参数进行拟合,计算换相电压的预测值。该方法能在系统故障后快速预测换流母线各相电压,有助于及时制定和采取紧急控制措施,对维持电力系统稳定运行具有重要意义。
附图说明
图1为本发明实施例的方法流程示意图;
图2为本发明实施例的CIGRE高压直流输电基准模型示意图。
具体实施方式
如图1所示,本发明实施例提供的一种基于二阶电路响应的交直流系统换相电压预测方法,包括:
步骤S1:基于二阶电路响应形式,建立交直流系统暂态响应的预测模型;
步骤S2:收集换流母线三相电压数据;
步骤S3:计算换流母线电压三序分量的傅里叶展开系数;
步骤S4:基于最小二乘法原理,对响应预测模型中的参数进行拟合;
步骤S5:预测换相电压三相分量。
本实施例以CIGRE高压直流输电基准模型为例进行分析,其模型结构如图2所示。该测试系统在PSCAD/EMTDC中搭建,直流侧电压500kV,直流容量为1000MW。故障设置在逆变侧换流母线处,1.0s发生a相短路,故障持续时间为0.1s,接地电阻为20Ω。
进一步的,步骤S1中交直流系统暂态响应预测模型的建立具体可分为如下步骤:
步骤S11:根据电力系统参数及储能元件的特点,得到暂态过渡过程中的时域响应形式:
式中,f(t)为所求的响应,p1,2为暂态过渡期间交直流系统等值RLC电路的特征根,a(t)为f(t)的稳态分量,包括直流及各倍频交流分量;为f(t)的暂态分量;τ、ωd、K为与系统参数及初值条件相关的常数;j为虚数单位;
步骤S12:根据下式计算f(t)稳态分量及暂态分量的傅里叶展开系数:
式中,ω0=100π(rad/s)为系统基频角速度,T=2π/ω0表示采样周期;An为稳态分量n倍频分量的系数,Xn为暂态分量n倍频分量的系数;特殊的,n=0表示直流分量,n=1表示基频分量;
式中,μ0、μ1、μ2、ν0、ν1、ν2为计算的中间变量,具体可根据下式计算得到:
进一步的,在本实施例中,步骤S2所述的收集换流母线三相电压数据,其具体实现方法为:以故障发生/恢复时刻为零时刻,以采样频率fs采集[0tβ]时间段内换流母线的三相电压ea、eb、ec;其中tβ表示换相开始时刻;本实施例中fs=0.1MHz,tβ=0.07945s;
进一步的,本实施例中步骤S3所述的计算换流母线电压三序分量的傅里叶展开系数,其具体实现步骤为:
步骤S31:根据下式计算所需采样的波形总数N及一个采样周波内的采样点数M:
式中:floor表示向下取整函数;本实施例中,N=119,M=2000;
式中,n表示所求傅里叶展开系数对应的频次;
进一步的,本实施例中步骤S4所述的响应预测模型参数的拟合,具体实现方法为:
进一步的,本实施例中步骤S5所述的换相电压预测的具体实现步骤为:
定义均方根误差ERMSE以及平均绝对误差EMAE为来衡量三项电压预测结果的准确性,两种指标的计算方式如下:
表1换相过程三相电压各频次分量预测误差对比
由表中数据可知,在故障恢复后第4周波换相过程中三相电压各频次分量的预测值误差较小。
Claims (6)
1.基于二阶电路响应的交直流系统换相电压预测方法,其特征在于,包括如下步骤:
(1)基于二阶电路响应形式,建立交直流系统暂态响应的预测模型;
(2)收集换流母线三相电压数据;
(3)计算换流母线电压三序分量的傅里叶展开系数;
(4)基于最小二乘法原理,对响应预测模型中的参数进行拟合;
(5)预测换相电压。
2.根据权利要求1所述的基于二阶电路响应的交直流系统换相电压预测方法,其特征在于,所述步骤(1)中交直流系统暂态响应预测模型的建立包括如下步骤:
(11)根据电力系统参数及储能元件的特点,得到暂态过程中的时域响应形式:
(12)根据下式计算f(t)稳态分量及暂态分量的傅里叶展开系数:
式中,ω0=100π(rad/s)为系统基频角速度,T=2π/ω0表示采样周期;An为稳态分量n倍频分量的系数,Xn为暂态分量n倍频分量的系数;特殊的,n=0表示直流分量,n=1表示基频分量,λ表示积分计算当中的积分元;
式中,λ表示积分计算当中的积分元,μ0、μ1、μ2、ν0、ν1、ν2为计算的中间变量,具体可根据下式计算得到:
3.根据权利要求1所述的基于二阶电路响应的交直流系统换相电压预测方法,其特征在于,所述步骤(2)中的收集换流母线三相电压数据的具体实现方法为:以故障发生或恢复时刻为零时刻,以采样频率fs采集[0 tβ]时间段内换流母线的三相电压ea、eb、ec;其中tβ表示换相开始时刻。
4.根据权利要求1所述的基于二阶电路响应的交直流系统换相电压预测方法,其特征在于,所述步骤(3)中计算换流母线电压三序分量的傅里叶展开系数的具体实现步骤为:
(31)根据下式计算所需采样的波形总数N及一个采样周波内的采样点数M:
式中:floor表示向下取整函数,T=2π/ω0表示采样周期,△T表示两次采样波形的时间间隔,fs表示采样频率,tβ表示换相开始时刻;
式中,M表示一个采样周波内的采样点数,n表示所求傅里叶展开系数对应的频次,e表示相电压,下标i取a、b、c分别表示a、b、c三相,△T表示两次采样波形的时间间隔,fs表示采样频率,ε表示第ε个采样点;
6.根据权利要求1所述的基于二阶电路响应的交直流系统换相电压预测方法,其特征在于,所述步骤(5)中换相电压预测的具体实现步骤为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910592330.7A CN110518615B (zh) | 2019-07-03 | 2019-07-03 | 基于二阶电路响应的交直流系统换相电压预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910592330.7A CN110518615B (zh) | 2019-07-03 | 2019-07-03 | 基于二阶电路响应的交直流系统换相电压预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110518615A CN110518615A (zh) | 2019-11-29 |
CN110518615B true CN110518615B (zh) | 2021-04-30 |
Family
ID=68623464
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910592330.7A Active CN110518615B (zh) | 2019-07-03 | 2019-07-03 | 基于二阶电路响应的交直流系统换相电压预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110518615B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111884180B (zh) * | 2020-07-15 | 2021-05-18 | 华北电力大学 | 一种基于最小二乘拟合的直流系统保护方法及装置 |
CN112564157B (zh) * | 2020-12-15 | 2023-12-26 | 深圳供电局有限公司 | 一种抑制连续换相失败的定熄弧角控制改进方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014141441A1 (ja) * | 2013-03-14 | 2014-09-18 | 東芝三菱電機産業システム株式会社 | 電流形電力変換装置 |
CN104333028A (zh) * | 2014-08-12 | 2015-02-04 | 许昌开普检测技术有限公司 | 预测型换相失败预防方法和装置 |
CN108964109A (zh) * | 2018-07-13 | 2018-12-07 | 湖南大学 | 一种抑制连续换相失败的协调控制方法 |
-
2019
- 2019-07-03 CN CN201910592330.7A patent/CN110518615B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014141441A1 (ja) * | 2013-03-14 | 2014-09-18 | 東芝三菱電機産業システム株式会社 | 電流形電力変換装置 |
CN104333028A (zh) * | 2014-08-12 | 2015-02-04 | 许昌开普检测技术有限公司 | 预测型换相失败预防方法和装置 |
CN108964109A (zh) * | 2018-07-13 | 2018-12-07 | 湖南大学 | 一种抑制连续换相失败的协调控制方法 |
Non-Patent Citations (2)
Title |
---|
Commutation Failure Analysis in Multi-Infeed HVDC Systems;Ebrahim Rahimi;《IEEE TRANSACTIONS ON POWER DELIVERY》;20110131;第26卷(第01期);第378-384页 * |
直流换相失败期间阀换相过程微观分析方法;李伟;《电力自动化设备》;20170302;第37卷(第03期);第115-119页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110518615A (zh) | 2019-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109412191B (zh) | 一种用于高压直流输电系统的锁相方法、装置和设备 | |
CN110518615B (zh) | 基于二阶电路响应的交直流系统换相电压预测方法 | |
CN112564157B (zh) | 一种抑制连续换相失败的定熄弧角控制改进方法 | |
Sun et al. | Analysis of the DC-link capacitor current of power cells in cascaded H-bridge inverters for high-voltage drives | |
CN110518616B (zh) | 基于一阶电路响应的交直流系统换相电压预测方法 | |
CN102307004A (zh) | 一种基于lcl滤波的可控整流器参数辨识方法 | |
CN110718932B (zh) | 一种考虑波形畸变及直流电流变化的换相失败预测方法 | |
CN103728506A (zh) | Hvdc系统换流变压器铁芯饱和型谐波不稳定判定的方法 | |
CN104535820A (zh) | 基于fbd法的三相有源电力滤波器谐波电流检测方法 | |
CN111969639A (zh) | 电力电子化电网级联型换流器多时间尺度暂态建模方法 | |
CN113078672B (zh) | 一种多馈入直流系统换相失败预测控制协调控制改进方法 | |
CN110705031B (zh) | 一种基于二阶泰勒系数的励磁涌流识别方法 | |
Noda et al. | A practical steady-state initialization method for electromagnetic transient simulations | |
CN113567808B (zh) | 一种统一潮流控制器接入线路故障定位方法及系统 | |
Li et al. | Disturbance Suppression for Mission Profile Emulator for Sub-Modules of Modular Multilevel Converter Under Nearest Level Modulation with Feedforward Control | |
Pugliese et al. | High-frequency harmonic current control of power converters | |
Li et al. | Model predictive control of NPC three-level grid-tied converter based on reconstructed current | |
Zhou et al. | Zero-sequence Circulating Current Suppression Strategy for Microgrid Inverters | |
Wasynczuk | Analysis of line-commutated converters during unbalanced operating conditions | |
Yang et al. | Mission profile emulator for sub-modules of modular multilevel converter with carrier phase-shift modulation | |
CN104122432A (zh) | 一种适用于畸变或不平衡情况下三电平有源滤波器参考电流检测方法 | |
CN116430109B (zh) | 电网不平衡状态下基于降阶谐振器的正负序分量检测方法 | |
Zhao et al. | Computer simulation and measurements of HVDC harmonics | |
CN112968461B (zh) | 一种基于偏磁模拟的换流变选相合闸励磁涌流抑制方法 | |
Guihua et al. | Synchronous Control Method of Three-phase Inverter Under Unbalanced Conditions Using the Energy Operator |
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 | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20201203 Address after: 213300 room 428, building a, 218 Hongkou Road, Kunlun Street, Liyang City, Changzhou City, Jiangsu Province (in Zhongguancun Science and Technology Industrial Park, Jiangsu Province) Applicant after: Liyang Research Institute of Southeast University Applicant after: SOUTHEAST University Address before: Four pailou Nanjing Xuanwu District of Jiangsu Province, No. 2 210096 Applicant before: SOUTHEAST University |
|
GR01 | Patent grant | ||
GR01 | Patent grant |