CN111999563B - 一种多级联直流变换器阻抗在线测量方法 - Google Patents

一种多级联直流变换器阻抗在线测量方法 Download PDF

Info

Publication number
CN111999563B
CN111999563B CN202010903727.6A CN202010903727A CN111999563B CN 111999563 B CN111999563 B CN 111999563B CN 202010903727 A CN202010903727 A CN 202010903727A CN 111999563 B CN111999563 B CN 111999563B
Authority
CN
China
Prior art keywords
converter
sine
signal
frequency
impedance
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
Application number
CN202010903727.6A
Other languages
English (en)
Other versions
CN111999563A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202010903727.6A priority Critical patent/CN111999563B/zh
Publication of CN111999563A publication Critical patent/CN111999563A/zh
Application granted granted Critical
Publication of CN111999563B publication Critical patent/CN111999563B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R27/00Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
    • G01R27/02Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02MAPPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
    • H02M3/00Conversion of dc power input into dc power output
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02MAPPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
    • H02M1/00Details of apparatus for conversion
    • H02M1/0067Converter structures employing plural converter units, other than for parallel operation of the units on a single load
    • H02M1/007Plural converter units in cascade

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

本发明公开了一种多级联直流变换器阻抗在线测量方法,通过多扰动源的协调配合,实现快速准确的多级联直流变换器阻抗特性提取。不同于传统的单扰动源方法,本发明提出多扰动源的位置、频谱能量分布协调配合阻抗在线测量方法,通过设计其频率避免了系统中同时存在多个扰动源而引发的相互干扰问题,并由系统中存在的直流变换器生成扰动,可以实现更为快速的阻抗在线测量,该方法在多级联直流系统中具有更好的优越性,同时也适用于并联系统、交流系统等各种场景。

Description

一种多级联直流变换器阻抗在线测量方法
技术领域
本发明属于可再生能源并网技术领域,具体涉及一种多级联直流变换器阻抗在线测量方法。
背景技术
DC-DC变流系统采用了多级联变换器结构,以实现高效,可靠的电源转换,级联转换器之间的动态交互可能会导致其动态性能下降并影响系统的稳定性,通常采用基于阻抗模型的小信号稳定性分析来评估多级联直流系统的稳定裕度,获得直流变换器的精确阻抗模型是获得准确分析结果的重要步骤。
可以使用两种方法获得直流转换器的阻抗特性:解析建模方法和阻抗测量方法,解析建模方法需要被测对象的详细内部参数和控制结构,这在实际应用中是未知的,而阻抗测量是一种可以仅根据输出端口电压和电流信号获得阻抗特性而无需内部参数等的阻抗模型提取方法;因此,阻抗测量对实际系统的稳定性分析具有重要意义。
正弦信号、二进制序列信号和多正弦信号可用作阻抗测量中的扰动信号。其中,正弦信号可确保最佳的信噪比,但是测量效率较低;二进制序列信号是一种宽频带信号,可以将多个频率分量注入被测变换器,但是由于二进制序列的频谱能量分布难以精确控制,在测量时易出现频谱泄漏、谐波干扰等问题,从而影响测量精度;此外,多正弦信号同样包含多个频率分量,并且多正弦信号的频谱能量分布精确可控,不存在频谱泄漏与谐波干扰的问题。目前对于变换器的阻抗测量多采用二进制序列,由于二进制序列所具有的频谱泄漏、谐波干扰等问题,结合多级联系统的结构特征,二进制序列不适用于多级联直流变换系统的实时阻抗测量。
目前,对于直流变换器的阻抗测量多采用单个扰动源注入的扰动注入方式,如文献[P.Li and B.Lehman,Performance prediction of DC-DC converters withimpedances as loads[J].IEEE Trans.Power Electron.,vol.19,no.1,pp.201–209,Jan.2004.]以及文献[D.Martin,E.Santi,A.Barkley,“Wide bandwidth systemidentification of AC system impedances by applying perturbations to anexisting converter,”Energy Conversion Congress&Exposition.Phoneix,AZ,USA,Sept.2011,pp.2549-2556.]都有公开介绍,但在对由多个级联变换器构成的系统内的阻抗测量时,由于级联变换器之间的直流电容具有解耦作用,仅凭单扰动源注入扰动信号难以在多级联系统中所有变换器端口产生足够强的激励信号,从而无法获取准确可靠的测量结果。为了确保准确的测量结果,根据现有技术,需要将扰动源连接在各个变换器端口分多次注入扰动信号,显然该方法测量效率较低;此外,由于频繁的负载切换和变化的工作条件,需要快速阻抗测量以确保在测量过程中系统工作状态不发生变化,因此现有方法的低效性使得对进行系统的实时稳定性分析难以实现;并且,采用将扰动源连接于不同的测量端口从而进行阻抗测量的方法,在实际操作时较为复杂,不仅需要额外的扰动源,还需要将扰动源依次连接在不同测量端口,在实际中需要耗费大量的人力、时间与经济成本。
发明内容
鉴于上述,本发明提出了一种多级联直流变换器阻抗在线测量方法,通过多扰动源的协调配合,实现快速准确的多级联直流变换器阻抗特性提取。
一种多级联直流变换器阻抗在线测量方法,包括如下步骤:
(1)对于多级联直流变换器中的任一级变换器,为其设计一个频域解耦的多正弦扰动信号;
(2)以多正弦扰动信号峰值系数最小为目标,计算确定多正弦扰动信号各频率分量的相位;
(3)选择确定多正弦扰动信号各频率分量的幅值;
(4)将构造好的多正弦扰动信号与原输出电压参考值叠加后作为新的输出电压参考值并对变换器加以控制,使得变换器的输出端口生成多正弦扰动信号作为下一级变换器的扰动源;
(5)对变换器输入端口的电压及电流信号进行傅里叶分析,通过测量提取电压及电流信号中对应的分量,进而计算出变换器对应各频率分量的阻抗。
进一步地,所述多正弦扰动信号的表达式如下:
Figure BDA0002660664380000031
其中:x(t)为t时刻的多正弦扰动信号值,an为多正弦扰动信号中第n个频率分量的幅值,fn为多正弦扰动信号中第n个频率分量的频率,
Figure BDA0002660664380000032
为多正弦扰动信号中第n个频率分量的相位,n为自然数且1≤n≤N,N为已知测量频段内频率分量的数量,t为自然数。
进一步地,所述多正弦扰动信号的幅频特性表示如下:
Figure BDA0002660664380000033
其中:Xi(k)为第i级变换器的多正弦扰动信号在k Hz频率点处的初始幅值,k为第i级变换器的多正弦扰动信号中第n个频率分量的频率fn,i为自然数且1≤i≤K,K为多级联直流变换器中的变换器数量。
进一步地,所述多正弦扰动信号的峰值系数表达式如下:
Figure BDA0002660664380000034
其中:CF为多正弦扰动信号的峰值系数(用于描述相位与峰值的关系,CF值越小越有利于阻抗测量的准确性),xmax为多正弦扰动信号的峰值,xRMS为多正弦扰动信号的有效值,Tm为多正弦扰动信号的周期即各频率分量对应周期的最小公倍数。
进一步地,所述步骤(2)中以多正弦扰动信号峰值系数CF最小为目标,采用有效集优化方法寻找使CF值最小的各频率分量相位组合,所述有效集优化方法通过计算有效约束推进优化过程,优化过程中保证相位
Figure BDA0002660664380000035
的取值范围如下:
Figure BDA0002660664380000036
进一步地,所述步骤(3)中对于任一级变换器的多正弦扰动信号,首先令其中各频率分量的幅值相等,进而在满足以下关系式的前提下选择最大的幅值作为多正弦扰动信号各频率分量的幅值;
xmax=max|x(t)|<20%Unom
其中:xmax为多正弦扰动信号的峰值,Unom为变换器的额定电压。
进一步地,所述步骤(5)中对于变换器对应第n个频率分量的阻抗,则先提取变换器输入端口的电压及电流信号中对应频率为fn的分量Un和In,然后通过以下公式计算出变换器对应第n个频率分量的阻抗Z[fn];
Z[fn]=Un/In
本发明多级联直流变换器阻抗在线测量方法与传统的阻抗测量方法相比,具有快速与实时性的特点;不同于传统的单扰动源方法,本发明提出多扰动源的位置、频谱能量分布协调配合阻抗在线测量方法,通过设计其频率避免了系统中同时存在多个扰动源而引发的相互干扰问题,并由系统中存在的直流变换器生成扰动,可以实现更为快速的阻抗在线测量,该方法在多级联直流系统中具有更好的优越性,同时也适用于并联系统、交流系统等各种场景。
附图说明
图1为多级联直流变换器阻抗在线测量原理示意图。
图2为频域解耦的多正弦信号幅频特性示意图。
图3(a)为本发明设计的多正弦信号1的幅频特性示意图。
图3(b)为本发明设计的多正弦信号2的幅频特性示意图。
图3(c)为本发明设计的多正弦信号3的幅频特性示意图。
图3(d)为本发明设计的多正弦信号4的幅频特性示意图。
图4(a)为采用本发明测量方法下直流变换器1的阻抗测量结果示意图。
图4(b)为采用本发明测量方法下直流变换器2的阻抗测量结果示意图。
图4(c)为采用本发明测量方法下直流变换器3的阻抗测量结果示意图。
图4(d)为采用本发明测量方法下直流变换器4的阻抗测量结果示意图。
具体实施方式
为了更为具体地描述本发明,下面结合附图及具体实施方式对本发明的技术方案进行详细说明。
本发明多级联直流变换器阻抗在线测量方法,包括如下步骤:
(1)根据多级联变换器数目K设计K个频域解耦的多正弦扰动信号所包含频率分量。
首先,多正弦信号的时域表达式定义如下:
Figure BDA0002660664380000051
其中:an表示构成多正弦信号频率分量的幅值,fn表示各分量的频率,
Figure BDA0002660664380000055
表示各分量相位。
以K=4为例,4个频域解耦的多正弦信号的幅频特性可表示为:
Figure BDA0002660664380000052
其中:n=0,1,2,3…,i表示频域解耦的多正弦信号的序号,k表示多正弦信号中频率分量的频率,Xi(k)表示第i个多正弦信号的在k Hz的幅值,该式可由图2表示。
(2)通过有效集优化方法设计多正弦信号各频率分量的相位。
本发明中以下式作为相位的设计指标:
Figure BDA0002660664380000053
其中:CF是多正弦信号的峰值系数,用于描述相位与多正弦信号峰值的关系,CF越小越有利于阻抗测量的准确性,xmax是多正弦信号峰值,xRMS是多正弦信号有效值,Tm是多正弦信号周期。
可以通过有效集优化方法,寻找使CF最小的相位组合,该方法通过计算有效约束推进优化过程,具有较快的迭代收敛速度与较强的跳出局部最优的能力。变量
Figure BDA0002660664380000056
的取值范围是:
Figure BDA0002660664380000054
(3)设计多正弦信号各频率分量幅值。
为保证测量结果的准确性,多正弦信号峰值应满足下式:
xmax=max|x(t)|<20%Unom
其中:Unom是直流变换器额定电压。
在上述步骤中已经完成对频率和相位的设计,同时考虑令所有分量幅值相等,则多正弦信号各分量幅值的最大值可以通过计算得到,并将该值作为多正弦信号各分量的幅值。
(4)将各个频域解耦的多正弦信号分别与一个直流变换器的输出电压参考值叠加,在各变换器的输出端口生成多正弦扰动信号,作为下一级变换器的扰动源,如图1所示。
(5)对各变换器输入端口电压与电流信号进行傅里叶分析,提取测量分量,计算变换器阻抗。
将注入的扰动信号频率定义为fp,则应提取被测直流变换器输入端口电压电流信号中频率为fp的分量Up,Ip,并由下式计算得到直流变换器在fp Hz的阻抗:
Z[fp]=Up/Ip
本实施例考虑对包含4个直流变换器的多级联系统内直流变换器的阻抗测量,测量频段为0~100Hz,其参数如表1所示:
表1
Figure BDA0002660664380000061
Figure BDA0002660664380000071
(1)首先,可由下式设计4个频域解耦的多正弦信号所包含的频率分量:
Figure BDA0002660664380000072
其中:n的最大值取50。
(2)确定频率以后,通过有效集优化方法对4个多正弦信号的CF进行优化,以寻找最优的相位组合,初始相位可由随机数获取,最终的多正弦信号的幅频特性如图3(a)~图3(d)所示,其中每个正弦分量的幅值暂时定义为1V。图中4个优化后的多正弦信号峰值分别为8.80V、8.34V、8.96V、7.96V,峰值系数CF分别为1.23、1.17、1.25、1.11。
(3)在完成相位设计以后,通过4个变换器的额定电压可以计算出4个多正弦信号的峰值上限分别为:0.2×750=150V,0.2×375=75V,0.2×110=22V,0.2×75=15V;由此得到的多正弦信号各分量幅值最大值为:150/8.80=17V,75/8.34=9V,22/8.96=2.5V,15/7.96=1.2V。
(4)将多正弦信号叠加到各直流变换器的输出电压控制环中直流电压指令值上,则可以实现扰动信号的输出控制。
(5)提取直流变换器输入端口电压电流信号中与扰动频率对应的分量,通过下式计算得到被测变换器阻抗:
Z[fp]=Up/Ip
其中:fp表示扰动频率,Up表示电压扰动频率分量,Ip表示电流扰动频率分量,Z[fp]表示直流变换器在fp Hz处的阻抗;对4个直流变换器阻抗测量而结果如图4(a)~图4(d)所示。
上述对实施例的描述是为便于本技术领域的普通技术人员能理解和应用本发明。熟悉本领域技术的人员显然可以容易地对上述实施例做出各种修改,并把在此说明的一般原理应用到其他实施例中而不必经过创造性的劳动。因此,本发明不限于上述实施例,本领域技术人员根据本发明的揭示,对于本发明做出的改进和修改都应该在本发明的保护范围之内。

Claims (1)

1.一种多级联直流变换器阻抗在线测量方法,包括如下步骤:
(1)对于多级联直流变换器中的任一级变换器,为其设计一个频域解耦的多正弦扰动信号,其表达式如下:
Figure FDA0003268303700000011
其中:x(t)为t时刻的多正弦扰动信号值,an为多正弦扰动信号中第n个频率分量的幅值,fn为多正弦扰动信号中第n个频率分量的频率,
Figure FDA0003268303700000014
为多正弦扰动信号中第n个频率分量的相位,n为自然数且1≤n≤N,N为已知测量频段内频率分量的数量,t为自然数;
所述多正弦扰动信号的幅频特性表示如下:
Figure FDA0003268303700000012
其中:Xi(k)为第i级变换器的多正弦扰动信号在k Hz频率点处的初始幅值,k为第i级变换器的多正弦扰动信号中第n个频率分量的频率fn,i为自然数且1≤i≤K,K为多级联直流变换器中的变换器数量;
所述多正弦扰动信号的峰值系数表达式如下:
Figure FDA0003268303700000013
其中:CF为多正弦扰动信号的峰值系数,xmax为多正弦扰动信号的峰值,xRMS为多正弦扰动信号的有效值,Tm为多正弦扰动信号的周期即各频率分量对应周期的最小公倍数;
(2)以多正弦扰动信号峰值系数CF最小为目标,采用有效集优化方法寻找使CF值最小的各频率分量相位组合,所述有效集优化方法通过计算有效约束推进优化过程,优化过程中保证相位
Figure FDA0003268303700000015
的取值范围如下:
Figure FDA0003268303700000021
(3)对于任一级变换器的多正弦扰动信号,首先令其中各频率分量的幅值相等,进而在满足以下关系式的前提下选择最大的幅值作为多正弦扰动信号各频率分量的幅值;
xmax=max|x(t)|<20%Unom
其中:Unom为变换器的额定电压;
(4)将构造好的多正弦扰动信号与原输出电压参考值叠加后作为新的输出电压参考值并对变换器加以控制,使得变换器的输出端口生成多正弦扰动信号作为下一级变换器的扰动源;
(5)对变换器输入端口的电压及电流信号进行傅里叶分析,通过测量提取电压及电流信号中对应的分量,进而计算出变换器对应各频率分量的阻抗,具体地:
对于变换器对应第n个频率分量的阻抗,则先提取变换器输入端口的电压及电流信号中对应频率为fn的分量Un和In,然后通过以下公式计算出变换器对应第n个频率分量的阻抗Z[fn];
Z[fn]=Un/In
CN202010903727.6A 2020-09-01 2020-09-01 一种多级联直流变换器阻抗在线测量方法 Active CN111999563B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010903727.6A CN111999563B (zh) 2020-09-01 2020-09-01 一种多级联直流变换器阻抗在线测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010903727.6A CN111999563B (zh) 2020-09-01 2020-09-01 一种多级联直流变换器阻抗在线测量方法

Publications (2)

Publication Number Publication Date
CN111999563A CN111999563A (zh) 2020-11-27
CN111999563B true CN111999563B (zh) 2022-03-22

Family

ID=73466103

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010903727.6A Active CN111999563B (zh) 2020-09-01 2020-09-01 一种多级联直流变换器阻抗在线测量方法

Country Status (1)

Country Link
CN (1) CN111999563B (zh)

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105337519B (zh) * 2015-11-18 2018-05-01 阳光电源股份有限公司 级联多电平变换器的自检系统及自检方法
CN105867513B (zh) * 2016-03-30 2017-07-28 北方工业大学 一种扰动观察最大功率点的跟踪装置
CN106885944B (zh) * 2017-03-17 2019-06-11 上海电力学院 一种用于下垂控制的直流微网线路阻抗检测方法及装置
CN108535545B (zh) * 2018-04-10 2019-08-09 湖南大学 一种双谐振注入式宽频带阻抗测量装置及其控制方法
CN109459615A (zh) * 2018-10-19 2019-03-12 湖南大学 基于级联型多电平变换器的高压阻抗测量装置及控制方法
CN109738703A (zh) * 2018-11-15 2019-05-10 湖南大学 高压宽频带阻抗测量装置及其控制方法
CN109742748B (zh) * 2019-02-28 2022-04-15 东南大学 一种适用于多电压等级直流配电系统的稳定判据方法
CN111506869B (zh) * 2020-04-15 2022-03-25 浙江大学 基于多正弦波信号注入的并网逆变器频率耦合阻抗特性提取方法

Also Published As

Publication number Publication date
CN111999563A (zh) 2020-11-27

Similar Documents

Publication Publication Date Title
Odavic et al. One-sample-period-ahead predictive current control for high-performance active shunt power filters
CN106164689B (zh) 用于多相dc/dc转换器的模型预测控制的方法和控制器
CN111506869B (zh) 基于多正弦波信号注入的并网逆变器频率耦合阻抗特性提取方法
CN113156214B (zh) 一种双模式宽频扰动装置及控制方法
CN117092409A (zh) 一种基于三相lcl并网逆变器的电网阻抗在线测量方法及装置
CN111999563B (zh) 一种多级联直流变换器阻抗在线测量方法
CN113315126B (zh) 一种有源电力滤波器指定次谐波抑制二次采样方法及系统
CN113839413A (zh) 一种交直流混合配电网逆变器参数设计方法
Salis et al. A generalised harmonic linearisation method for power converters input/output impedance calculation
Xiao et al. A Hybrid Data-driven Power Loss Minimization Method of Dual Active Bridge Converters
CN107257198A (zh) 一种电源效率的调节方法及电路
Rygg et al. Real-time stability analysis of power electronic systems
CN108362940B (zh) 一种半周波畸变信号的动态谐波频率提取方法
CN115236404B (zh) 一种并网逆变器端口阻抗自测量方法
CN117240131A (zh) 一种基于压电换能器参数的动态t型匹配方法及系统
CN110086197A (zh) 电感电流采样计算方法
CN115986745A (zh) 一种智能配电网谐波电压补偿方法
CN115728555A (zh) 一种抗频率耦合干扰的单相变压器宽频阻抗测量方法
Wang et al. A high-precision and wideband fundamental frequency measurement method for synchronous sampling used in the power analyzer
Pejović et al. An analysis of three-phase rectifiers with constant voltage loads
Antoniewicz et al. Comparison of current control strategies for four-leg shunt active power filter in Matlab-simulink
CN116644699B (zh) 基于电压小信号扰动的级联变换器母线电压稳定分析方法
Tao et al. A novel measurement approach based on wideband excitation for frequency-coupling admittances of train converter
CN114944649B (zh) 基于电量频谱的电网状态辨识方法
Zeid et al. Predicting DC-Link Capacitor Current Ripple in AC-DC Rectifier Circuits Using Fine-Tuned Large Language Models

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