CN105844063A - 一种特高压直流多数据源参数实测和校核方法 - Google Patents

一种特高压直流多数据源参数实测和校核方法 Download PDF

Info

Publication number
CN105844063A
CN105844063A CN201610280152.0A CN201610280152A CN105844063A CN 105844063 A CN105844063 A CN 105844063A CN 201610280152 A CN201610280152 A CN 201610280152A CN 105844063 A CN105844063 A CN 105844063A
Authority
CN
China
Prior art keywords
current
voltage
control module
parameter
delta
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
Application number
CN201610280152.0A
Other languages
English (en)
Other versions
CN105844063B (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.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Beijing Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Beijing 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 State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, State Grid Beijing Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201610280152.0A priority Critical patent/CN105844063B/zh
Publication of CN105844063A publication Critical patent/CN105844063A/zh
Application granted granted Critical
Publication of CN105844063B publication Critical patent/CN105844063B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • 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
    • 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)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Inverter Devices (AREA)

Abstract

本发明提供一种特高压直流多数据源参数实测和校核方法,其特征在于,所述方法包括对多种数据源依次进行分环节实测和整体校核两个步骤,所述分环节实测依据特高压直流模型的特点,将控制模块解耦,一个模块的参数根据一项试验数据实测得到;所述整体校核用故障扰动数据,校核电流控制的积分时间常数Ti_I,和低压限流控制模块的参数,同时从整体特性上校核其他各环节参数;所述多数据源包括实际电网试验、实际控保联调试验和控保数字模型仿真。本发明提升了直流输电仿真建模技术水平,解决了特高压直流仿真计算参数的确定缺乏依据的问题,填补技术空白。

Description

一种特高压直流多数据源参数实测和校核方法
技术领域
本发明涉及电力系统仿真与建模技术领域,具体涉及一种特高压直流多数据源参数实测和校核方法。
背景技术
特高压直流工程连续投运,交直流系统相互影响作用显著。特高压直流的暂态特性与其计算参数关系紧密,因此通过一定的方法获取准确的计算参数对特高压直流的机电暂态仿真具有重要意义。目前国内电力系统规划和运行计算中广泛使用的特高压直流模型是基于实际工程的控制特性建模的,具有普遍适用性。然而,模型所使用的计算参数却仍为一套典型参数,不能够体现不同直流工程之间的差异性,影响了交直流电网计算的准确性。主要的限制在于没有一套完整可行的参数获取方法。
发明内容
为了克服上述现有技术的不足,本发明提供一种特高压直流多数据源参数实测和校核方法。本发明提升了直流输电仿真建模技术水平,解决了特高压直流仿真计算参数的确定缺乏依据的问题,填补技术空白。
为了实现上述发明目的,本发明采取如下技术方案:
一种特高压直流多数据源参数实测和校核方法,所述方法包括对多种数据源依次进行分环节实测和整体校核两个步骤,所述分环节实测依据特高压直流模型的特点,将控制模块解耦,一个模块的参数根据一项试验数据实测得到;所述整体校核用故障扰动数据,校核电流控制的积分时间常数Ti_I,和低压限流控制模块的参数,同时从整体特性上校核其他各环节参数;所述多数据源包括实际电网试验、实际控保联调试验和控保数字模型仿真。
优选的,所述分环节实测包括如下步骤:
步骤1-1、根据直流电流阶跃试验数据,实测电流控制模块参数;
步骤1-2、根据直流电压阶跃试验数据,实测电压控制模块参数Kp_V,Ti_V;
步骤1-3、根据整流侧换流母线三相瞬时性短路试验数据,实测整流侧最小触发角控制模块参数K1_ra,K2_ra,Cdl,Dl,Decr;
步骤1-4、根据逆变侧换流母线三相瞬时性短路试验数据,实测换相失败预测模块参数K_cf,G_cf,Tdn_cf;
步骤1-5、根据直流线路短路试验数据,实测重启动控制模块参数Alpha_Ret,Th_ret,Alpha_Res,Th_Res。
优选的,所述整体校核包括如下步骤:
步骤2-1、校核电流控制模块的积分时间常数Ti_I,通过比较逆变测换流母线三相瞬时性短路试验条件下的直流电流动态变化曲线,增大或减小Ti_I以使直流电流的变化趋势吻合;
步骤2-2、校核低压限流控制模块的5个参数,所述参数包括:低压阈值Udlow、高压阈值Udhigh、上升时间常数Tup、下降时间常数Tdn和最小电流Iomin;通过比较逆变测换流母线三相瞬时性短路试验条件下的低压限流输出Iolim和直流换流器吸收的无功Qhvdc,获得低压限流控制模块参数的数值。
优选的,所述步骤1-1中,所述电流控制模块进行电流指令阶跃试验,阶跃量ΔIo,阶跃瞬间,根据电流偏差增益值Idiff与电流偏差Id-Io的变化量确定参数Gain:
G a i n = ΔI d i f f ΔI d - ΔI o - - - ( 1 )
式中:ΔId为直流电流的增量,ΔIo为直流电流指令的增量,ΔIdiff为直流电流与指令的偏差增量,Gain为电流控制总增益;
由于阶跃瞬时,直流电流不能突变,有ΔId=0,故:
G a i n = ΔI d i f f - ΔI o - - - ( 2 )
一次试验有向下、向上两次阶跃,根据两次实测值取平均即为最终的参数Gain,如下:
向下阶跃时:
Gain d n = ΔI d i f f d n ΔI o - - - ( 3 )
向上阶跃时:
Gain u p = ΔI d i f f u p - ΔI o - - - ( 4 )
二者取均值:
G a i n = Gain d n + Gain u p 2 - - - ( 5 )
式中,ΔIdiffdn为向下阶跃时直流电流与指令的偏差增量,ΔIdiffup为向上阶跃时直流电流与指令的偏差增量;
阶跃瞬间,根据αi_P确定Kp_I:
式中:αi_P为电流控制模块比例支路的输出角度,αi n-1为阶跃前电流控制模块输出的触发角指令,Kp_I为电流控制比例增益;
向下阶跃时:
向上阶跃时:
二者取均值:
K p _ I = K p _ I d n + K p _ I u p 2 - - - ( 9 )
式中,αi_Pdn为向下阶跃时电流控制模块比例支路的输出角度,αi_Pup为向上阶跃时电流控制模块比例支路的输出角度,Kp_Idn为向下阶跃时电流控制比例增益,Kp_Iup为向上阶跃时电流控制比例增益。
10.优选的,所述步骤1-2中,所述电压控制模块进行电压指令阶跃试验,阶跃量ΔUdref,忽略积分环节的变化,且直流电压不突变,故根据输出触发角的变化量Δαv与电压指令的变化量ΔUdref确定Kp_V:
K p _ V = Δα v ΔU d r e f - - - ( 10 )
式中:Δαv为电压控制模块的输出触发角度的增量,ΔUdref为直流电压指令的增量;
向下阶跃时,KP_V为电压控制比例增益:
K p _ V d n = Δα v d n - ΔU d r e f - - - ( 11 )
向上阶跃时:
K p _ V u p = Δα v u p ΔU d r e f - - - ( 12 )
二者取均值:
K p _ V = K p _ V d n + K p _ V u p 2 - - - ( 13 )
式中,Δαvdn为向下阶跃时电压控制模块的输出触发角度的增量,Δαvup为向上阶跃时电压控制模块的输出触发角度的增量;将电压控制积分时间常数Ti_V设置为直流工程设计值,逐步调整Ti_V,当模型的仿真曲线与实测曲线吻合程度最佳时,则确定Ti_V。
优选的,所述步骤1-3中,所述整流侧最小触发角控制模块进行整流侧换流母线三相瞬时性短路试验,逐步减小故障电阻,使得故障期间整流侧换流母线电压跌落幅度逐步增大,根据整流最小触发角控制模块的输出角度αraml变化实测K1_ra、K2_ra、Cdl、Dl,当故障期间母线电压等于设定电压值Uacr1时,αraml由零出现第一级输出角度α1,则:
K1_ra=Uacr0-Uacr1 (14)
式中:Uac0时为故障前的初始换流母线电压,K1_ra为该模块启动的第一级交流母线电压阈值;
当障期间母线电压等于设定电压值Uac2时,αraml由α1增加为第二级输出角度α2,则:
K2_ra=Uacr0-Uacr2 (15)
式中,K2_ra为该模块启动的第二级交流母线电压阈值;
相应的有:
Cdl=α1 (16)
Dl=α2 (17)
式中,Cdl为第一级角度置位值;Dl为第二级角度置位值;
一种特殊情况是,金属性短路αraml仅有一级输出α1,即α2=α1,有:
K2_ra=K1_ra (18)
Dl=Cdl=α1 (19)
将输出角下降速率Decr设为工程设计值,逐步调整Decr的大小,当模型的仿真曲线的αraml与实测曲线最为吻合时,则确定Decr。
优选的,所述步骤1-4中,所述换相失败预测模块进行逆变侧换流母线三相瞬时性短路试验,逐步减小故障电阻,使故障期间换流母线电压跌落幅度逐步增大,根据换相失败预测模块的输出角度Δα变化实测K_cf和G_cf,当故障期间母线电压Uaci等于Uaci1时,Δα由零出现正值,有:
K_cf=1-(Uaci0-Uaci1) (20)
式中:Uaci0为换流母线初始电压,K_cf为环节启动的交流电压阈值;
若Δα在故障期间的最小值为A1,有:
G _ c f = 1 - cosA 1 U a c i 0 - U a c i 1 - - - ( 21 )
式中,G_cf为环节电压增益;将K_cf、G_cf代入模型,进行相同的短路故障,将角度下降时间常数Tdn_cf设为工程设计值,逐步调整Tdn_cf的大小,当模型的仿真曲线的Δα与实测曲线最为吻合时,则确定Tdn_cf
优选的,所述步骤1-5中,所述重启动控制模块进行直流线路瞬时性短路试验,若整流侧触发角α的第一级输出为α1dcf,维持时间T1dcf,第二级输出为α2dcf,维持时间T2dcf,有:
A l p h a _ Re t = α 1 d c f T h _ Re t = T 1 d c f A l p h a _ Re s = α 2 d c f T h _ Re s = T 2 d c f - - - ( 22 )
式中,Alpha_Ret为移相角度,Th_Ret为移相时间,Alpha_Res为重启动角度,Th_Res为重启动保持时间。
优选的,所述步骤2-1包括如下步骤:
步骤2-1-1、比较下降段的起始时刻t1和上升段的终止时刻t4,若仿真曲线较实测曲线t1超前,t4滞后,则减小参数高压阈值Udhigh,反之增大;
步骤2-1-2、比较下降段的终止时刻t2和上升段的起始时刻t3,若仿真曲线较实测曲线t2超前,t3滞后,则减小参数低压阈值Udlow,反之增大;
步骤2-1-3、调整参数Udhigh与Udlow后,比较下降段的斜率kdn,若仿真曲线较实测曲线平缓,则减小参数下降时间常数Tdn,反之增大;
步骤2-1-4、调整参数Udhigh与Udlow后,比较上升段的斜率kup,若仿真曲线较实测曲线平缓,则减小参数上升时间常数Tup,反之增大。
与现有技术相比,本发明的有益效果在于:
本发明提供的特高压直流多数据源参数实测和校核方法,提升了直流输电仿真建模技术水平,解决了特高压直流仿真计算参数的确定缺乏依据的问题,填补技术空白;
本发明与特高压直流仿真模型相结合,为特高压直流系统的工程设计,以及含有特高压直流的交直流电网的规划、运行提供了仿真手段与技术支撑。
附图说明
图1是本发明提供的一种特高压直流多数据源参数实测和校核方法的流程图
图2是本发明提供的电流控制模块逻辑框图;
图3是本发明提供的电压控制模块逻辑框图;
图4是本发明提供的整流侧最小触发角控制模块逻辑框图;
图5是本发明提供的换相失败预测模块逻辑框图;
图6是本发明提供的低压限流控制模块输出示意图;
图7是本发明提供的电流控制模块Idiff信号录波;
图8是本发明提供的电流控制模块αi_P信号录波;
图9是本发明提供的电流控制模块αi信号录波;
图10是本发明提供的向下阶跃时电压控制模块αv信号录波(局部);
图11是本发明提供的向上阶跃时电压控制模块αv信号录波(局部);
图12是本发明提供的直流电压仿真与试验曲线比较;
图13是本发明提供的整流侧最小触发角控制模块输出试验曲线;
图14是本发明提供的整流侧最小触发角控制模块输出仿真与试验曲线对比;
图15是本发明提供的换相失败预测模块输出试验曲线;
图16是本发明提供的换相失败预测模块输出仿真与试验曲线对比;
图17是本发明提供的重启动控制模块输出试验曲线;
图18是本发明提供的直流电流仿真与试验曲线对比(Ti_I=0.009);
图19是本发明提供的低压限流控制模块输出仿真与试验曲线对比(校核前参数);
图20是本发明提供的低压限流控制模块输出仿真与试验曲线对比(校核后参数)。
具体实施方式
下面结合附图对本发明作进一步详细说明。
基于目前电力系统规划和运行计算所使用的特高压直流模型,根据需求提取其中需要实测校核的参数。该模型对特高压直流工程实际的控制系统进行简化、等效,依据机电暂态仿真需求进行转化,包括主控制模块、低压限流控制模块、电流控制模块、电压控制模块、熄弧角控制模块、整流侧最小触发角控制模块、换相失败预测模块、电压恢复控制模块和重启动控制模块。其中,需要确定动态参数的模块有:电流控制模块、电压控制模块、整流侧最小触发角控制模块、换相失败预测模块、低压限流控制模块和重启动控制模块。参数列表如表1所示。
表1特高压直流模型计算参数列表
如图1所示,为特高压直流多数据源参数实测和校核方法,步骤如下:
一是分环节实测。根据模型特性,将上述模块解耦,一个模块的参数根据一项试验数据实测得到。试验数据来源于实际电网试验、实际直流控保联调试验或者控保数字模型的仿真试验。分环节实测根据环节数量又细分为5个步骤:
(1)根据直流电流阶跃试验数据,实测电流控制模块参数Gain,Kp_I;
(2)根据直流电压阶跃试验数据,实测电压控制模块参数Kp_V,Ti_V;
(3)根据整流侧换流母线三相瞬时性短路试验数据,实测整流侧最小触发角控制模块参数K1_ra,K2_ra,Cdl,Dl,Decr;
(4)根据逆变侧换流母线三相瞬时性短路试验数据,实测换相失败预测模块参数K_cf,G_cf,Tdn_cf;
(5)根据直流线路短路试验数据,实测重启动控制模块参数Alpha_Ret,Th_ret,Alpha_Res,Th_Res。
二是整体校核。采用较大的故障扰动数据,校核电流控制模块参数Ti_I,和低压限流控制模块的参数Udlow,Udhigh,Tup,Tdn,Iomin。同时从整体特性上校核其他各环节参数。
下面详细阐述参数实测和校核方法流程。
分环节实测
(1)电流控制模块
电流控制模块的逻辑框图如图2所示。进行电流指令阶跃试验,阶跃量ΔIo,维持一定时间后,阶跃至初始指令。阶跃瞬间,根据电流偏差增益值Idiff与电流偏差Id-Io的变化量确定参数Gain:
G a i n = ΔI d i f f ΔI d - ΔI o - - - ( 1 )
式中:ΔId为直流电流的增量,ΔIo为直流电流指令的增量,ΔIdiff为直流电流与指令的偏差增量,Gain为电流控制总增益;
由于阶跃瞬时,直流电流不能突变,有ΔId=0,故:
G a i n = ΔI d i f f - ΔI o - - - ( 2 )
一次试验有向下、向上两次阶跃,根据两次实测值取平均即为最终的参数Gain,如下:
向下阶跃时:
Gain d n = ΔI d i f f d n ΔI o - - - ( 3 )
向上阶跃时:
Gain u p = ΔI d i f f u p - ΔI o - - - ( 4 )
二者取均值:
G a i n = Gain d n + Gain u p 2 - - - ( 5 )
式中,ΔIdiffdn为向下阶跃时直流电流与指令的偏差增量,ΔIdiffup为向上阶跃时直流电流与指令的偏差增量;
阶跃瞬间,根据αi_P确定Kp_I:
式中:αi_P为电流控制模块比例支路的输出角度,αi n-1为阶跃前电流控制模块输出的触发角指令,Kp_I为电流控制比例增益;
向下阶跃时:
向上阶跃时:
二者取均值:
K p _ I = K p _ I d n + K p _ I u p 2 - - - ( 9 )
式中,αi_Pdn为向下阶跃时电流控制模块比例支路的输出角度,αi_Pup为向上阶跃时电流控制模块比例支路的输出角度,KP_Idn为向下阶跃时电流控制比例增益,KP_Iup为向上阶跃时电流控制比例增益。
(2)电压控制模块
11.电压控制模块的逻辑框图如图3所示。进行电压指令阶跃试验,阶跃量ΔUdref,忽略积分环节的变化,且直流电压不突变,故根据输出触发角的变化量Δαv与电压指令的变化量ΔUdref确定Kp_V:
K p _ V = Δα v ΔU d r e f - - - ( 10 )
式中:Δαv为电压控制模块的输出触发角度的增量,ΔUdref为直流电压指令的增量;
向下阶跃时,KP_V为电压控制比例增益:
K p _ V d n = Δα v d n - ΔU d r e f - - - ( 11 )
向上阶跃时:
K p _ V u p = Δα v u p ΔU d r e f - - - ( 12 )
二者取均值:
K p _ V = K p _ V d n + K p _ V u p 2 - - - ( 13 )
式中,Δαvdn为向下阶跃时电压控制模块的输出触发角度的增量,Δαvup为向上阶跃时电压控制模块的输出触发角度的增量;将电压控制积分时间常数Ti_V设置为直流工程设计值,逐步调整Ti_V,当模型的仿真曲线与实测曲线吻合程度最佳时,则确定Ti_V。
(3)整流侧最小触发角控制模块
整流侧最小触发角控制模块的逻辑框图如图4所示。参数含义如下:
K1_ra——该模块启动的第一级交流母线电压阈值;
K2_ra——该模块启动的第二级交流母线电压阈值;
Cdl——整流侧第一级置位角度;
Dl——整流侧第二级置位角度;
Decr——环节退出后输出角下降速率。
该模块在整流侧发生交流故障时启动,其逻辑为:正常状态下,整流最小触发角控制模块的输出角度αraml为零;当换流母线电压幅值Uacr低于第一级阈值K1_ra时,输出角度为αraml=Cdl;若母线电压继续降低至小于第二级阈值K2_ra时,输出角度为αraml=Dl,一般有Dl≥Cdl;故障清除,母线电压恢复至K1_ra以上,输出角按照一定的斜率(与Decr相关)逐渐下降至零。
根据上述逻辑,进行整流侧换流母线三相瞬时性短路试验,逐步减小故障电阻,使得故障期间整流侧换流母线电压跌落幅度逐步增大,根据αraml变化实测K1_ra、K2_ra、Cdl、Dl。当故障期间母线电压等于Uacr1时,αraml由零出现第一级输出角度α1,则:
K1_ra=Uacr0-Uacr1 (14)
式中:Uac0时为故障前的初始换流母线电压。
当障期间母线电压等于Uac2时,αraml由α1增加为第二级输出角度α2,则:
K2_ra=Uacr0-Uacr2 (15)
相应的有:
Cdl=α1 (16)
Dl=α2 (17)
一种特殊情况是,金属性短路αraml仅有一级输出α1,即α2=α1,有:
K2_ra=K1_ra (18)
Dl=Cdl=α1 (19)
在特高压直流机电暂态模型中将上述参数设置为实测值,进行相同的短路故障,通过拟合法确定参数Decr,即首先将Decr设为工程设计值,逐步调整Decr的大小,当模型的仿真曲线的αraml与实测曲线最为吻合时,即确定Decr。
(4)换相失败预测模块
换相失败预测模块的逻辑如图5所示。参数含义如下:
K_cf——环节启动的交流电压阈值;
G_cf——环节电压增益;
Tdn_cf——环节退出后角度下降滤波时间常数。
进行逆变侧换流母线三相瞬时性短路试验,逐步减小故障电阻,使故障期间换流母线电压跌落幅度逐步增大,根据模块的输出角度Δα变化实测K_cf和G_cf。当故障期间母线电压Uaci等于Uaci1时,Δα由零出现正值,有:
K_cf=1-(Uaci0-Uaci1) (20)
式中:Uaci0为换流母线初始电压。
若Δα在故障期间的最小值为A1,根据图5,有:
G _ c f = 1 - cosA 1 U a c i 0 - U a c i 1 - - - ( 21 )
在特高压直流机电暂态模型中将K_cf、G_cf设置为实测值,进行相同的短路故障,通过拟合法确定参数Tdn_cf,即首先将Tdn_cf设为工程设计值,逐步调整Tdn_cf的大小,当模型的仿真曲线的Δα与实测曲线最为吻合时,即确定Tdn_cf
(5)重启动控制模块
重启动控制的主要逻辑为:发生直流线路短路故障,整流侧控制器调节触发角移相至角度Alpha_Ret,并维持时间Th_Ret;随后,调节触发角至Alpha_Res进行重启并维持时间Th_Res。
进行直流线路瞬时性短路试验,根据整流侧触发角指令值αord的变化实测上述参数。若αord的第一级输出为α1dcf,维持时间T1dcf,第二级输出为α2dcf,维持时间T2dcf,有:
A l p h a _ Re t = α 1 d c f T h _ Re t = T 1 d c f A l p h a _ Re s = α 2 d c f T h _ Re s = T 2 d c f - - - ( 22 )
整体校核
(1)校核电流控制模块参数Ti_I
首先校核参数Ti_I,进行逆变测换流母线三相瞬时性短路试验,以典型值代入模型,比较直流电流动态变化的仿真曲线和实测曲线,增大或减小Ti_I以使直流电流的变化趋势吻合。
(2)校核低压限流控制模块参数
先给出一组典型参数(可以是设计参数),同样进行逆变测换流母线三相瞬时性短路试验,比较低压限流的输出值Iolim和直流换流器吸收的无功Qhvdc,根据实测曲线的变化趋势,校核低压限流控制的五个参数。具体方法结合图6说明如下。
图6为低压限流控制模块输出Iolim示意图。其动态过程由三个线段组成:下降段,保持段,上升段。保持段的值为I1,有:
Iomin=I1 (23)
第一步,比较下降段的起始时刻t1和上升段的终止时刻t4,若仿真曲线较实测曲线t1超前,t4滞后,则减小参数Udhigh,反之增大;
第二步,比较下降段的终止时刻t2和上升段的起始时刻t3,若仿真曲线较实测曲线t2超前,t3滞后,则减小参数Udlow,反之增大;
第三步,调整参数Udhigh与Udlow后,比较下降段的斜率kdn,若仿真曲线较实测曲线平缓,则减小参数Tdn,反之增大;
第三步,调整参数Udhigh与Udlow后,比较上升段的斜率kup,若仿真曲线较实测曲线平缓,则减小参数Tup,反之增大。
进一步比较直流系统的无功功率,对低压限流控制模块的参数进一步校核。
经过上述的分环节实测和整体校核流程,可以获得特高压直流模型的全部计算参数。该方法基于实际电网试验、实际控保联调试验、控保数字模型仿真等多种数据源完成,优先采用实际电网试验数据,对受客观条件限制现场无法进行的试验,采用实际控保联调试验和控保数字模型仿真的数据。多数据源之间相互校核,可以确保参数的科学、准确、合理。
以下通过几个具体的数据实例,对该方法进行具体说明。
分环节实测部分
(1)电流控制模块
Idiff录波如图7所示。向下阶跃时,依据式(3),有:
Gain d n = 2.9572 - 0.0074 0.1 = 29.498 - - - ( 24 )
向上阶跃时:
Gain u p = - 3.1215 + 0.0635 - 0.1 = 30.58 - - - ( 25 )
二者取均值:
G a i n = Gain d n + Gain u p 2 = 30.04 - - - ( 26 )
αi_P与αi的录波如图8、图9所示。向下阶跃时,依据式(4),有:
向上阶跃时:
二者取均值:
K p _ I = 2.75 + 2.847 2 = 2.799 - - - ( 29 )
(2)电压控制模块
αv录波如图10所示,向下阶跃时:
K p _ V d n = 139.16 - 139.313 - 34 / 800 = 3.6 - - - ( 30 )
如图11所示,向上阶跃时:
K p _ V u p = 136.877 - 136.684 42 / 800 = 3.676 - - - ( 31 )
二者取均值:
K p _ V = 3.6 + 3.676 2 = 3.64 - - - ( 32 )
当Ti_V=0.0005时,仿真与实测的直流电压曲线最为吻合,如图12所示。
(3)整流最小触发角控制模块
整流侧换流母线三相瞬时性短路条件下,αraml实测曲线如图13所示。故障期母线电压Uacr=0.9pu时,αraml始终为零;Uacr=0.8pu时,αraml在故障期间为30度,说明模块启动,有:
K1_ra=0.85,Cdl=30 (33)
又Uacr=0.0即金属性短路时,αraml在故障期间亦为30度,依据式(18),(19),有:
K2_ra=0.85,Dl=30 (34)
通过拟合法确定参数Decr,仿真曲线与实测曲线对比如图14所示。当Decr=0.37时,仿真曲线与实测曲线吻合,因此:
Decr=0.37 (35)
(4)换相失败预测模块
在逆变侧换流母线三相瞬时性短路条件下,Δα实测曲线如图15所示。当故障期间母线电压Uaci最低为0.8pu时,模块不启动;Uaci最低0.7pu时,模块启动,初始母线电压Uaci0=1.0pu,依据式(20),有:
K_cf=0.75 (36)
依据式(21),有:
通过拟合法确定参数Tdn_cf。仿真曲线与实测曲线对比如图16所示。Tdn_cf=0.02时,仿真曲线与实测曲线吻合,则:
Tdn_cf=0.02 (38)
(5)重启动控制模块
直流线路短路故障条件下,整流侧触发角实测曲线如图17所示。依据曲线特征和式(22),有:
整体校核部分
校核电流控制模块的积分时间常数Ti_I。该值对直流电流Id的变化趋势有关键的影响。首先给出Ti_I初始值0.007。进行逆变侧三相短路试验,调整Ti_I=0.009,进一步比较,,如图18所示。调整后仿真与实测曲线更为吻合,Ti_I的取值大小对直流电流的横向时间位置产生较大影响,取值越小,故障后直流电流下降至零以及开始上升(恢复)的时间均较早。因此:
Ti_I=0.009 (40)
校核低压限流控制模块的参数。给出初始参数VDCOL参数为Udlow=0.15,Udhigh=0.8,Tup=0.04,Tdn=0.015,Iomin=0.345。在初始参数下,逆变侧三相短路的仿真曲线与实测曲线对比如图19所示。低压限流的输出有一定的差距,依据前文给出的校核方法,对上述参数进行调整,调整后为:
Udlow=0.25,Udhigh=0.7,Tup=0.04,Tdn=0.018,Iomin=0.345 (41)
调整参数后,调整后仿真与实测对比如图20所示。仿真曲线与实测曲线的吻合程度得到较大提升。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求范围当中。

Claims (9)

1.一种特高压直流多数据源参数实测和校核方法,其特征在于,所述方法包括对多种数据源依次进行分环节实测和整体校核两个步骤,所述分环节实测依据特高压直流模型的特点,将控制模块解耦,一个模块的参数根据一项试验数据实测得到;所述整体校核用故障扰动数据,校核电流控制的积分时间常数Ti_I,和低压限流控制模块的参数,同时从整体特性上校核其他各环节参数;所述多数据源包括实际电网试验、实际控保联调试验和控保数字模型仿真。
2.根据权利要求1所述方法,其特征在于,所述分环节实测包括如下步骤:
步骤1-1、根据直流电流阶跃试验数据,实测电流控制模块参数;
步骤1-2、根据直流电压阶跃试验数据,实测电压控制模块参数;
步骤1-3、根据整流侧换流母线三相瞬时性短路试验数据,实测整流侧最小触发角控制模块参数;
步骤1-4、根据逆变侧换流母线三相瞬时性短路试验数据,实测换相失败预测模块参数;
步骤1-5、根据直流线路短路试验数据,实测重启动控制模块参数。
3.根据权利要求1所述方法,其特征在于,所述整体校核包括如下步骤:
步骤2-1、校核电流控制模块的积分时间常数Ti_I,通过比较逆变测换流母线三相瞬时性短路试验条件下的直流电流动态变化曲线,增大或减小Ti_I以使直流电流的变化趋势吻合;
步骤2-2、校核低压限流控制模块的5个参数,所述参数包括:低压阈值Udlow、高压阈值Udhigh、上升时间常数Tup、下降时间常数Tdn和最小电流Iomin;通过比较逆变测换流母线三相瞬时性短路试验条件下的低压限流输出Iolim和直流换流器吸收的无功Qhvdc,获得低压限流控制模块参数的数值。
4.根据权利要求2所述方法,其特征在于,所述步骤1-1中,所述电流控制模块进行电流指令阶跃试验,阶跃量ΔIo,阶跃瞬间,根据电流偏差增益值Idiff与电流偏差Id-Io的变化量确定参数Gain:
G a i n = ΔI d i f f ΔI d - ΔI o - - - ( 1 )
式中:ΔId为直流电流的增量,ΔIo为直流电流指令的增量,ΔIdiff为直流电流与指令的偏差增量,Gain为电流控制总增益;
由于阶跃瞬时,直流电流不能突变,有ΔId=0,故:
G a i n = ΔI d i f f - ΔI o - - - ( 2 )
一次试验有向下、向上两次阶跃,根据两次实测值取平均即为最终的参数Gain,如下:
向下阶跃时:
Gain d n = ΔI d i f f d n ΔI o - - - ( 3 )
向上阶跃时:
Gain u p = ΔI d i f f u p - ΔI o - - - ( 4 )
二者取均值:
G a i n = Gain d n + Gain u p 2 - - - ( 5 )
式中,ΔIdiffdn为向下阶跃时直流电流与指令的偏差增量,ΔIdiffup为向上阶跃时直流电流与指令的偏差增量;
阶跃瞬间,根据αi_P确定Kp_I:
式中:αi_P为电流控制模块比例支路的输出角度,αi n-1为阶跃前电流控制模块输出的触发角指令,Kp_I为电流控制比例增益;
向下阶跃时:
向上阶跃时:
二者取均值:
K p _ I = K p _ I d n + K p _ I u p 2 - - - ( 9 )
式中,αi_Pdn为向下阶跃时电流控制模块比例支路的输出角度,αi_Pup为向上阶跃时电流控制模块比例支路的输出角度,Kp_Idn为向下阶跃时电流控制比例增益,Kp_Iup为向上阶跃时电流控制比例增益,为向上阶跃时电流控制模块输出的触发角指令,为向下阶跃时电流控制模块输出的触发角指令。
5.根据权利要求2所述方法,其特征在于,所述步骤1-2中,所述电压控制模块进行电压指令阶跃试验,阶跃量ΔUdref,忽略积分环节的变化,且直流电压不突变,故根据输出触发角的变化量Δαv与电压指令的变化量ΔUdref确定Kp_V:
K p _ V = Δα v ΔU d r e f - - - ( 10 )
式中:Δαv为电压控制模块的输出触发角度的增量,ΔUdref为直流电压指令的增量;
向下阶跃时,KP_V为电压控制比例增益:
K p _ V d n = Δα v d n - ΔU d r e f - - - ( 11 )
向上阶跃时:
K p _ V u p = Δα v u p ΔU d r e f - - - ( 12 )
二者取均值:
K p _ V = K p _ V d n + K p _ V u p 2 - - - ( 13 )
式中,Kp_Vdn为向下阶跃时电压控制比例增益,Kp_Vup为向上阶跃时电压控制比例增益,Δαvdn为向下阶跃时电压控制模块的输出触发角度的增量,Δαvup为向上阶跃时电压控制模块的输出触发角度的增量;将电压控制积分时间常数Ti_V设置为直流工程设计值,逐步调整Ti_V,当模型的仿真曲线与实测曲线吻合程度最佳时,则确定Ti_V。
6.根据权利要求2所述方法,其特征在于,所述步骤1-3中,所述整流侧最小触发角控制模块进行整流侧换流母线三相瞬时性短路试验,逐步减小故障电阻,使得故障期间整流侧换流母线电压跌落幅度逐步增大,根据整流最小触发角控制模块的输出角度αraml变化实测K1_ra、K2_ra、Cdl、Dl,当故障期间母线电压等于设定电压值Uacr1时,αraml由零出现第一级输出角度α1,则:
K1_ra=Uacr0-Uacr1 (14)
式中:Uacr0时为故障前的初始换流母线电压,K1_ra为该模块启动的第一级交流母线电压阈值;
当障期间母线电压等于设定电压值Uacr2时,αraml由α1增加为第二级输出角度α2,则:
K2_ra=Uacr0-Uacr2 (15)
式中,K2_ra为该模块启动的第二级交流母线电压阈值;
相应的有:
Cdl=α1 (16)
Dl=α2 (17)
式中,Cdl为第一级角度置位值;Dl为第二级角度置位值;
一种特殊情况是,金属性短路αraml仅有一级输出α1,即α2=α1,有:
K2_ra=K1_ra (18)
Dl=Cdl=α1 (19)
将输出角下降速率Decr设为工程设计值,逐步调整Decr的大小,当模型的仿真曲线的αraml与实测曲线最为吻合时,则确定Decr。
7.根据权利要求2所述方法,其特征在于,所述步骤1-4中,所述换相失败预测模块进行逆变侧换流母线三相瞬时性短路试验,逐步减小故障电阻,使故障期间换流母线电压跌落幅度逐步增大,根据换相失败预测模块的输出角度Δα变化实测K_cf和G_cf,当故障期间母线电压Uaci等于设定电压值Uaci1时,Δα由零出现正值,有:
K_cf=1-(Uaci0-Uaci1) (20)
式中:Uaci0为换流母线初始电压,K_cf为环节启动的交流电压阈值;
若Δα在故障期间的最小值为A1,有:
G _ c f = 1 - cos A 1 U a c i 0 - U a c i 1 - - - ( 21 )
式中,G_cf为环节电压增益;将K_cf、G_cf代入模型,进行相同的短路故障,将角度下降时间常数Tdn_cf设为工程设计值,逐步调整Tdn_cf的大小,当模型的仿真曲线的Δα与实测曲线最为吻合时,则确定Tdn_cf
8.根据权利要求2所述方法,其特征在于,所述步骤1-5中,所述重启动控制模块进行直流线路瞬时性短路试验,若整流侧触发角α的第一级输出为α1dcf,维持时间T1dcf,第二级输出为α2dcf,维持时间T2dcf,有:
A l p h a _ Re t = α 1 d c f T h _ Re t = T 1 d c f A l p h a _ Re s = α 2 d c f T h _ Re s = T 2 d c f - - - ( 22 )
式中,Alpha_Ret为移相角度,Th_Ret为移相时间,Alpha_Res为重启动角度,Th_Res为重启动保持时间。
9.根据权利要求3所述方法,其特征在于,所述步骤2-1包括如下步骤:
步骤2-1-1、比较下降段的起始时刻t1和上升段的终止时刻t4,若仿真曲线较实测曲线t1超前,t4滞后,则减小参数高压阈值Udhigh,反之增大;
步骤2-1-2、比较下降段的终止时刻t2和上升段的起始时刻t3,若仿真曲线较实测曲线t2超前,t3滞后,则减小参数低压阈值Udlow,反之增大;
步骤2-1-3、调整参数Udhigh与Udlow后,比较下降段的斜率kdn,若仿真曲线较实测曲线平缓,则减小参数下降时间常数Tdn,反之增大;
步骤2-1-4、调整参数Udhigh与Udlow后,比较上升段的斜率kup,若仿真曲线较实测曲线平缓,则减小参数上升时间常数Tup,反之增大。
CN201610280152.0A 2016-04-28 2016-04-28 一种特高压直流多数据源参数实测和校核方法 Active CN105844063B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610280152.0A CN105844063B (zh) 2016-04-28 2016-04-28 一种特高压直流多数据源参数实测和校核方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610280152.0A CN105844063B (zh) 2016-04-28 2016-04-28 一种特高压直流多数据源参数实测和校核方法

Publications (2)

Publication Number Publication Date
CN105844063A true CN105844063A (zh) 2016-08-10
CN105844063B CN105844063B (zh) 2023-05-26

Family

ID=56590415

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610280152.0A Active CN105844063B (zh) 2016-04-28 2016-04-28 一种特高压直流多数据源参数实测和校核方法

Country Status (1)

Country Link
CN (1) CN105844063B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106934164A (zh) * 2017-03-17 2017-07-07 广东电网有限责任公司电力科学研究院 一种实时电磁暂态仿真的直流控制系统建模方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014121438A1 (en) * 2013-02-05 2014-08-14 Alstom Technology Ltd. Method and apparatus for current differential protection for uhvdc transmission line
CN104113056A (zh) * 2014-06-30 2014-10-22 南方电网科学研究院有限责任公司 一种低压限流控制参数的优化方法
CN104600734A (zh) * 2014-12-30 2015-05-06 华南理工大学 一种高压直流输电低压限流与pi控制环节协调优化方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014121438A1 (en) * 2013-02-05 2014-08-14 Alstom Technology Ltd. Method and apparatus for current differential protection for uhvdc transmission line
CN104113056A (zh) * 2014-06-30 2014-10-22 南方电网科学研究院有限责任公司 一种低压限流控制参数的优化方法
CN104600734A (zh) * 2014-12-30 2015-05-06 华南理工大学 一种高压直流输电低压限流与pi控制环节协调优化方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
万磊等: "特高压直流控制系统机电暂态等效建模与参数实测方法", 《电网技术》 *
荀吉辉等: "稳定计算用发电机励磁系统模型参数测试及校核分析", 《电力系统保护与控制》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106934164A (zh) * 2017-03-17 2017-07-07 广东电网有限责任公司电力科学研究院 一种实时电磁暂态仿真的直流控制系统建模方法
CN106934164B (zh) * 2017-03-17 2020-03-20 广东电网有限责任公司电力科学研究院 一种实时电磁暂态仿真的直流控制系统建模方法

Also Published As

Publication number Publication date
CN105844063B (zh) 2023-05-26

Similar Documents

Publication Publication Date Title
US9570909B2 (en) Devices and methods for decentralized power loss reduction control
US8838284B2 (en) Devices and methods for decentralized Volt/VAR control
US8965588B2 (en) Devices and methods for decentralized voltage control
US9411389B2 (en) Distributed generation control for microgrid during islanding
US8838285B2 (en) Devices and methods for decentralized power factor control
CN107147107B (zh) 一种抑制多直流连锁换相失败的调相机布点方法
US20130030599A1 (en) Devices and methods for decentralized coordinated volt/var control
CN100553065C (zh) 消弧线圈系统单相接地补偿状态在线监测及调控方法
CN101944742B (zh) 改进的潮流可行解恢复方法
CN104104102B (zh) 电压源换流器型多端直流输电系统稳态工作点优化方法
Mather et al. Distribution-connected PV's response to voltage sags at transmission-scale
CN108400597B (zh) 计及安控和一次调频特性的直流故障静态安全分析方法
CN109327037B (zh) 一种分层接入直流换相失败预防控制整定方法及装置
CN108683209B (zh) 一种分布式电源并网能力评估方法和装置
Rahmann et al. Methodology for selecting operating points and contingencies for frequency stability studies
Malamaki et al. Ramp-rate limitation control of distributed renewable energy sources via supercapacitors
CN110208649B (zh) 基于交流电压下降速度的换相故障检测方法及装置
CN105844063A (zh) 一种特高压直流多数据源参数实测和校核方法
CN111564870B (zh) 一种提高电力系统暂态稳定性的紧急切机控制策略表整定方法
CN110601215B (zh) 考虑连续换相失败约束的动态无功支撑能力评估方法及系统
CN102244411A (zh) 一种ups并机系统及其并机方法
JP7289748B2 (ja) 能動的方式を不使用としたpcsを連系した系統の単独運転検出限界を求める評価方法
CN114937994A (zh) 一种交直流混联电网的直流闭锁后稳定控制方法及装置
CN109885889B (zh) 一种同系列动态无功补偿装置的仿真评估方法及系统
Yang et al. Simplified Prediction Model of Frequency Nadir for Power Systems Penetrated With Renewable Energy

Legal Events

Date Code Title Description
C06 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