CN111682542B - 一种三相整流电路的谐波耦合外特性预测方法 - Google Patents

一种三相整流电路的谐波耦合外特性预测方法 Download PDF

Info

Publication number
CN111682542B
CN111682542B CN202010422150.7A CN202010422150A CN111682542B CN 111682542 B CN111682542 B CN 111682542B CN 202010422150 A CN202010422150 A CN 202010422150A CN 111682542 B CN111682542 B CN 111682542B
Authority
CN
China
Prior art keywords
current
voltage
harmonic
phase
matrix
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
CN202010422150.7A
Other languages
English (en)
Other versions
CN111682542A (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.)
Xinyang Power Supply Co of State Grid Henan Electric Power Co Ltd
Original Assignee
Xinyang Power Supply Co of State Grid Henan 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 Xinyang Power Supply Co of State Grid Henan Electric Power Co Ltd filed Critical Xinyang Power Supply Co of State Grid Henan Electric Power Co Ltd
Priority to CN202010422150.7A priority Critical patent/CN111682542B/zh
Publication of CN111682542A publication Critical patent/CN111682542A/zh
Application granted granted Critical
Publication of CN111682542B publication Critical patent/CN111682542B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/01—Arrangements for reducing harmonics or ripples
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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/18—Arrangements for adjusting, eliminating or compensating reactive power in networks
    • H02J3/1821—Arrangements for adjusting, eliminating or compensating reactive power in networks using shunt compensators
    • H02J3/1835—Arrangements for adjusting, eliminating or compensating reactive power in networks using shunt compensators with stepless control
    • H02J3/1842—Arrangements for adjusting, eliminating or compensating reactive power in networks using shunt compensators with stepless control wherein at least one reactive element is actively controlled by a bridge converter, e.g. active filters
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02M—APPARATUS 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
    • H02M7/00—Conversion of ac power input into dc power output; Conversion of dc power input into ac power output
    • H02M7/02—Conversion of ac power input into dc power output without possibility of reversal
    • H02M7/04—Conversion of ac power input into dc power output without possibility of reversal by static converters
    • H02M7/12—Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
    • H02M7/145—Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a thyratron or thyristor type requiring extinguishing means
    • H02M7/155—Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a thyratron or thyristor type requiring extinguishing means using semiconductor devices only
    • H02M7/162—Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a thyratron or thyristor type requiring extinguishing means using semiconductor devices only in a bridge configuration
    • H02M7/1623—Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a thyratron or thyristor type requiring extinguishing means using semiconductor devices only in a bridge configuration with control circuit
    • H02M7/1626—Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a thyratron or thyristor type requiring extinguishing means using semiconductor devices only in a bridge configuration with control circuit with automatic control of the output voltage or current
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/04—Power grid distribution networks
    • 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
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • 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
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/20—Active power filtering [APF]
    • 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
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/40—Arrangements for reducing harmonics

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (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)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种三相整流电路的谐波耦合外特性预测方法，包括如下三个步骤：步骤一，获取数据：根据电路拓扑，测得电感L、电阻R、导通电阻r，并且记录相应的控制触发角α的数据；步骤二，建立矩阵：分别对直流侧电压、直流侧电流以及交流侧电流进行傅里叶变换和矩阵分解；步骤三，预测与补偿：当模型预测得到的电压电流谐波含量大于标准中的限值，对电路进行谐波补偿；而低于限值时进行脉冲封锁。本发明基于矩阵推导清晰呈现了谐波耦合导纳矩阵建模过程，通过建立开关函数移相角和矩阵模型参数之间的直接对应关系，改善了基于移相矩阵变换的矩阵模型参数传统算法，为研究和使用非线性负载矩阵模型，提高电网谐波分析精确度提供了便利。

Description

一种三相整流电路的谐波耦合外特性预测方法
技术领域
本发明属于阻感负载整流电路谐波外特性建模及其相关的电能质量治理领域，涉及一种三相整流电路的谐波耦合外特性预测方法。
背景技术
随着社会经济的发展和人民生活水平的提高，除了满足电量的需求，用户对电能质量的要求也越来越高。电能质量治理技术对电网的安全、可靠、经济运行至关重要。基于电力电子变换器的有源治理设备由于控制灵活、精度高、不易受电网参数影响，代表着电能质量治理设备的发展方向。
在传统的电力系统分析中，由于谐波问题早期并不显著，负载模型主要反映的是负载在基波频率的外特性，通常是负载功率(有功和无功)、电流、电压以及频率之间关系的数学表达式。根据此时系统的状态可将负载模型分为静态模型和动态模型。其中静态模型是针对系统状态没有明显变化的情况，较多的用于电力系统的潮流分析，静稳定分析中，如恒阻抗模型、恒功率模型、恒电流模型及其组合等，而表达形式上主要有多项式模型，幂函数模型等；动态模型则主要关注系统状态的暂态过程，主要用于动态计算，稳定性判断等，动态模型引入了时间变量，一般用微分方程和差分方程表示。感应电动机作为动态负载的主要成分，是最主要的动态负载模型。
随着电力工业特别是电力电子技术的发展，基于半导体开关技术的负载在给人们生活带来便利的同时也使谐波等问题变得愈发显著，也受到了越来越多的关注。然而传统的电力系统负载模型缺少对负载谐波情况的描述，因此很多学者开始关注以晶闸管三相整流电路为代表的非线性负载的谐波模型。
建立适用于电网谐波分析的非线性负载模型是研究局域配电网谐波治理的基础，本专利基于移相矩阵变换，使用一种非线性负载谐波耦合导纳矩阵计算谐波耦合外特性参数。
传统电网谐波分析主要基于非线性负载的恒流源模型。恒流源模型是最简单且使用最广泛的模型。该模型考虑非线性负载的谐波电流只和交流侧基波电压有关，而忽略了谐波电压的影响，将非线性负载看作是内阻无穷大的谐波电流源。恒流源模型应用方便，但由于忽略了谐波电压的影响，当电网中电压畸变较严重时，误差会显著增大。并且，一些电流非常容易受到影响的负载不存在典型频谱。谐波耦合导纳矩阵模型基于非线性负载的工作原理，将电压电流的时域关系变换为频域上各次分量之间的矩阵关系，具有很高的精确度，但模型计算复杂，且模型参数易受外电压和控制参数影响而实时变化，阻碍了模型的应用。
发明内容
本发明的目的在于对电流连续的三相整流电路进行建模，从而在电能质量治理的研究中纳入各次谐波间互相耦合的影响因素；本发明以电感滤波的阻感负载整流电路的谐波源为例，进行谐波耦合导纳矩阵建模的设计。
本发明采用的技术方案如下：
一种三相整流电路的谐波耦合外特性预测方法，包括如下步骤：
步骤一，获取数据：
根据电路拓扑，测得电感L和电阻R的数据，通过器件手册查询晶闸管的导通电阻r，并且记录相应的控制触发角α的数据；
步骤二，建立矩阵：
分别对直流侧电压、直流侧电流以及交流侧电流进行傅里叶变换和矩阵分解；
步骤三，预测与补偿：
根据步骤二中矩阵得出的谐波耦合外特性预测模型，预测此时整流电路电压电流的谐波含量，并依据国标《GBT 14549—93电能质量公用电网谐波标准》判断谐波含量是否超标，当模型预测得到的电压电流谐波含量大于标准中的限值，会开启有源电力滤波器设备对电路进行谐波补偿，低于限值时有源电力滤波器将处于脉冲封锁，即关机状态。
步骤一的具体步骤是，通过阻抗测量仪Bode-100测得电感L和电阻R的数据，而其中控制触发角α的数据是设定值，通过对晶闸管控制触发程序对其设定而产生。
步骤二的具体步骤如下：
第一步，将整流桥中晶闸管控制信号，即开关函数分解为傅里叶级数；
第二步，计算直流侧电压；
第三步，计算直流侧电流；
第四步，计算交流侧电流。
在第一步中，考虑换相过程后，开关函数的傅里叶级数为：
Figure GDA0003846883520000021
Figure GDA0003846883520000022
上式中，ω代表电网频率，约为314.15rad/s,α为预设的开关函数移相角，μ为换相重叠角，可根据移相角α，电路形式及网侧变压器漏感推算；
由于三相平衡，通过延时1/3和2/3个工频周期T，可以得到B相和C相的开关函数sub(t)和suc(t)。
在第二步中，直流侧电压的时域计算表达式如下：
udc(t)＝ua(t)sua(t)+ub(t)sub(t)+uc(t)suc(t) (10)
其中，ua(t)，ub(t)，uc(t)分别为三相电压的模型假定时域值，无需测量，其大小仅和电压幅值有关；
其中，sua(t)，sub(t)，suc(t)分别为三相电压开关函数；
由于三相平衡，直流侧电压只含有直流分量和6倍频分量，在这些频率上，上式中三部分结果相同，可以先求A相对应的结果再乘3，而其他频率分量三相结果之和为0，这里具体推导证明省略。
在MATLAB中定义四个零矩阵
Figure GDA0003846883520000031
和
Figure GDA0003846883520000032
其中Cdc1为交流电压的各次分量到直流电压直流分量的变换矩阵，Cdc2为交流电压的各次分量的共轭到直流电压直流分量的变换矩阵，C1为交流电压的各次分量到直流电压各次分量的变换矩阵。利用MATLAB程序对A相电压的各次分量与开关函数各次分量的相乘情况进行判断，只保留结果中直流和50次以下6倍频的分量，并乘以3，将每种情况对应的变换系数叠加到矩阵Cdc1，Cdc2，C1，C2的对应位置元素上；
最终直流侧电压频域表达式为：
Figure GDA0003846883520000033
而第三步的计算方法是，在正常导通过程中，电流流经两个晶闸管，存在两个导通电阻，而在换相过程中，由于3个晶闸管同时导通，相当于流经1.5个导通电阻，这里忽略换相时导通电阻的变化，认为一直存在两个导通电阻r，将其归到直流侧，则直流侧电阻包含负载电阻R和两个晶闸管导通电阻r，设直流侧总电感为L，则直流侧电流的直流分量以及50次以内各次谐波分量为：
Figure GDA0003846883520000041
式中
Figure GDA0003846883520000042
在第四步中，交流侧电流的时域表达式如下：
Figure GDA0003846883520000043
其中，idc(t)为直流侧电流，ia(t)，ib(t)，ic(t)分别为三相电流，sia(t)，sib(t)，sic(t)分别为三相电流开关函数。
对于三相平衡情况，只需计算A相，其他两相通过电流的正负序关系可以求得。在MATLAB中定义两个零矩阵
Figure GDA0003846883520000044
和
Figure GDA0003846883520000045
其中C3为直流侧电流的各次分量到交流侧电流各次分量的变换矩阵，C4直流侧电流的各次分量的共轭到交流侧电流各次分量的变换矩阵。利用MATLAB程序对每个电流谐波分量与开关函数谐波分量相乘的情况进行判断，只保留结果中频率在50次谐波以下的情况，将每种情况对应的变换系数叠加到C3，C4的对应元素上，则最终交流侧电流频域表达式为：
Figure GDA0003846883520000046
对于整流器的交流侧谐波电压和电流，只有典型次非零，所以如果只关注典型谐波次，式(15)可以写为：
Figure GDA0003846883520000047
综上所述，由于采用了上述技术方案，本发明的有益效果是：
本发明详细分析了典型非线性阻感负载在电压畸变环境下的参数关系和开关函数特点，基于矩阵推导清晰呈现了谐波耦合导纳矩阵建模过程，通过建立开关函数移相角和矩阵模型参数之间的直接对应关系，改善了基于移相矩阵变换的矩阵模型参数传统算法，在基本不损失模型精度的情况下显著的减少计算时间，为研究和使用非线性负载矩阵模型，提高电网谐波分析精确度提供了便利。
附图说明
图1为本发明中的三相阻感负载移相整流器的结构及控制框图；
图2考换相过程的A相电压和电流开关函数；
图3实例1中模型及实验得到的各次电流含量对比；
图4实例2中模型及实验得到的各次电流含量对比。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
一种三相整流电路的谐波耦合外特性预测方法，包括如下步骤：
步骤一，获取数据：如图1所示的电路拓扑，测得电感L和电阻R的数据，通过器件手册查询晶闸管的导通电阻r，并且记录相应的控制触发角α的数据。
具体方法是，通过阻抗测量仪Bode-100测得电感L和电阻R的数据。
其中控制触发角α的数据是设定值，通过对晶闸管控制触发程序对其设定而产生。
步骤二，建立矩阵：分别对直流侧电压、直流侧电流以及交流侧电流进行傅里叶变换和矩阵分解。
步骤二的具体步骤如下：
第一步，将整流桥中晶闸管控制信号，即开关函数分解为傅里叶级数；各项开关函数仅和触发角有关，其波形如图2和图3所示；
考虑换相过程后，开关函数的傅里叶级数为：
Figure GDA0003846883520000051
Figure GDA0003846883520000061
上式中，ω代表电网频率，约为314.15rad/s,α为预设的开关函数移相角，μ为换相重叠角，可根据移相角α，电路形式及网侧变压器漏感推算。
由于三相平衡，通过延时1/3和2/3个工频周期T，可以得到B相和C相的开关函数sub(t)和suc(t)。
根据上一节的分析，当交流侧电压或控制触发角改变时，开关函数可能会超前或滞后一定的时间，体现在移相角的变化上，这将进一步导致矩阵模型参数发生改变。
第二步，计算直流侧电压。
直流侧电压的时域计算表达式如下：
udc(t)＝ua(t)sua(t)+ub(t)sub(t)+uc(t)suc(t) (10)
其中，ua(t)，ub(t)，uc(t)分别为三相电压的模型假定时域值，无需测量，其大小仅和电压幅值有关，波形图如图4所示。
其中，sua(t)，sub(t)，suc(t)分别为三相电压开关函数。
由于三相平衡，直流侧电压只含有直流分量和6倍频分量，在这些频率上，上式中三部分结果相同，可以先求A相对应的结果再乘3，而其他频率分量三相结果之和为0，这里具体推导证明省略。
在MATLAB中定义四个零矩阵
Figure GDA0003846883520000062
其中Cdc1为交流电压的各次分量到直流电压直流分量的变换矩阵，Cdc2为交流电压的各次分量的共轭到直流电压直流分量的变换矩阵，C1为交流电压的各次分量到直流电压各次分量的变换矩阵。利用MATLAB程序对A相电压的各次分量与开关函数各次分量的相乘情况进行判断，只保留结果中直流和50次以下6倍频的分量，并乘以3，将每种情况对应的变换系数叠加到矩阵Cdc1，Cdc2，C1，C2的对应位置元素上。
最终直流侧电压频域表达式为：
Figure GDA0003846883520000071
第三步，计算直流侧电流。
在正常导通过程中，电流流经两个晶闸管，存在两个导通电阻，而在换相过程中，由于3个晶闸管同时导通，相当于流经1.5个导通电阻，这里忽略换相时导通电阻的变化，认为一直存在两个导通电阻r，将其归到直流侧，则直流侧电阻包含负载电阻R和两个晶闸管导通电阻r，设直流侧总电感为L，则直流侧电流的直流分量以及50次以内各次谐波分量为：
Figure GDA0003846883520000072
式中
Figure GDA0003846883520000073
第四步，计算交流侧电流。
交流侧电流的时域表达式如下：
Figure GDA0003846883520000074
其中，idc(t)为直流侧电流，ia(t)，ib(t)，ic(t)分别为三相电流，sia(t)，sib(t)，sic(t)分别为三相电流开关函数。
对于三相平衡情况，只需计算A相，其他两相通过电流的正负序关系可以求得。在MATLAB中定义两个零矩阵
Figure GDA0003846883520000075
和
Figure GDA0003846883520000076
其中C3为直流侧电流的各次分量到交流侧电流各次分量的变换矩阵，C4直流侧电流的各次分量的共轭到交流侧电流各次分量的变换矩阵。利用MATLAB程序对每个电流谐波分量与开关函数谐波分量相乘的情况进行判断，只保留结果中频率在50次谐波以下的情况，将每种情况对应的变换系数叠加到C3，C4的对应元素上，则最终交流侧电流频域表达式为：
Figure GDA0003846883520000081
对于整流器的交流侧谐波电压和电流，只有典型次非零，所以如果只关注典型谐波次，式(15)可以写为：
Figure GDA0003846883520000082
步骤三，预测与补偿：工作人员提前输入直流侧负载R,L的值，晶闸管控制芯片中的程序将以设定值触发角α及计算得到的换相重叠角μ的值为基础，根据步骤二中矩阵得出的谐波耦合外特性预测模型，预测此时整流电路电压电流的谐波含量，并依据国标《GBT14549—93电能质量公用电网谐波标准》判断谐波含量是否超标，当模型预测得到的电压电流谐波含量大于标准中的限值，会开启有源电力滤波器设备对电路进行谐波补偿，低于限值时有源电力滤波器将处于脉冲封锁，即关机状态。
实施例1
首先通过仿真对提出的三相阻感负载整流电路谐波耦合导纳矩阵模型算法进行测试。
步骤一，获取数据。
得到的数据如下表所示：
基波有效值/V 负载电感/μH 负载电阻/Ω 触发角/degree
20 820 11.7 20
具体获取数据的方法是：测试电路为一个三相平衡的电压源对一个三相晶闸管整流电路供电，如图1所示；整流侧负载电阻R＝11.7Ω，负载电感L＝820μH，电压源的相电压有效值为30V，整流电路的控制触发角为20°，交流测进线电感为0H，因而不考虑换相角，。
步骤二，建立谐波耦合导纳矩阵：
1)将预定的触发角α代入式(8)和(9)，计算a,b,c三相晶闸管控制开关函数sua(t)，sub(t)，suc(t)，sia(t)，sib(t)，sic(t)的傅里叶级数展开式；
2)根据上一步得到的傅里叶级数展开式，并结合式(10)和式(11)，得到直流侧电压频域1到50次谐波分量：
Figure GDA0003846883520000091
3)根据上一步得到的直流侧电压谐波分量，并代入式(12)和式(13)得到直流侧电流频域1到50次谐波分量：
Figure GDA0003846883520000092
4)根据式(15)计算交流测电流频域表达式，并改写为式(16)的形式，最终得到交流测电流各次谐波分量为：
I＝[-1.96-5.91i,0,0,0,1.51+0.19i,0,0.227-0.68i,...,0.089-0.034i,0]。
为了验证模型得到的交流侧谐波电流含量的准确性。在实验室中，用一台三相可编程交流源(CHROMA 61705)给一个三相晶闸管移相整流器供电，并使用功率分析仪(HIOKI3197)对整流器的交流侧电压电流进行录波，用于和模型对比。整流器的直流侧参数和仿真相同，可以发现基于模型理论计算的数据和实验实测数据在不同情况下都呈现出良好的一致性。对应的实验参数和最终的误差仅为1.41％，模型算法和实验分别得到的各次谐波含量结果对比，如图3所示。
步骤三，预警与补偿：
依据国标《GBT 14549—93电能质量公用电网谐波标准》判断谐波含量是否超标，将实验中交流侧基波线电压值折算至10kV，短路容量为100MVA,计算出折算后的整流桥交流测各次谐波的含量，并计算在此条件下电压谐波占比约为(仅展示含量最大的前4次谐波)：
5次谐波占比/V 5次谐波占比/V 7次谐波占比/V 11次谐波占比/V
2.94％ 1.37％ 1.2％ 0.75％
折算后各次电压谐波含量均小于标准中规定的3.2％，即在此状态下无需示警，谐波补偿设备无需动作。
实施例2
首先通过仿真对提出的三相阻感负载整流电路谐波耦合导纳矩阵模型算法进行测试。
步骤一，获取数据：
得到数据如下表所示：
基波有效值/V 负载电感/μH 负载电阻/Ω 触发角/degree
20 820 11.7 0
获取数据的具体方法是，测试电路为一个三相平衡的电压源对一个三相晶闸管整流电路供电，如图1所示；整流侧负载电阻R＝11.7Ω，负载电感L＝820μH，电压源的相电压有效值为20V，整流电路的控制触发角为60°，交流测进线电感为1mH，因而需要考虑换相重叠角。
步骤二，建立谐波耦合导纳矩阵：
1)将预定的触发角α和换相重叠角μ代入式(8)和(9)，计算a,b,c三相晶闸管控制开关函数sua(t)，sub(t)，suc(t)，sia(t)，sib(t)，sic(t)的傅里叶级数展开式；
2)根据上一步得到的傅里叶级数展开式，并结合式(10)和式(11)，得到直流侧电压频域1到50次谐波分量标幺值(以交流测基波电流有效值为基准值)：
Figure GDA0003846883520000101
3)根据上一步得到的直流侧电压谐波分量，并代入式(12)和式(13)得到直流侧电流频域1到50次谐波分量：
Figure GDA0003846883520000102
4)根据式(15)计算交流测电流频域表达式，并改写为式(16)的形式，最终得到交流测电流各次谐波分量为：
I＝[-1.47-8.021i,0,0,0,1.42+1.15i,0,0.67+0.15i,...,0.01-0.0134i,0]。
为了验证模型得到的交流侧谐波电流含量的准确性。在实验室中，用一台三相可编程交流源(CHROMA 61705)给一个三相晶闸管移相整流器供电，并使用功率分析仪(HIOKI3197)对整流器的交流侧电压电流进行录波，用于和模型对比。整流器的直流侧参数和仿真相同，可以发现基于模型理论计算的数据和实验实测数据在不同情况下都呈现出良好的一致性。对应的实验参数和最终的误差仅为0.67％。模型算法得到的结果和实验结果对比，如图4所示。
大量测试都证明，本研究提出的三相阻感负载整流电路的谐波耦合导纳矩阵模型改进传统算法能在基本不损失模型精度的情况下显著的减少计算时间，为研究阻感负载整流电路的谐波特性，进而优化其谐波的治理方式提供了便利。
步骤三，预警与补偿：
依据国标《GBT 14549—93电能质量公用电网谐波标准》判断谐波含量是否超标，将实验中交流侧基波线电压值折算至10kV，短路容量为100MVA,计算出折算后的整流桥交流测各次谐波的含量，并计算在此条件下电压谐波占比约为(仅展示含量最大的前4次谐波)：
5次谐波占比/V 5次谐波占比/V 7次谐波占比/V 11次谐波占比/V
3.52％ 1.33％ 0.97％ 0.53％
折算后5次电压谐波含量大于标准中规定的3.2％，在此状态下需示警，并控制谐波补偿装置对该次谐波电流进行补偿。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (2)

1.一种三相整流电路的谐波耦合外特性预测方法，其特征在于，包括如下步骤：
步骤一，获取数据：
根据电路拓扑，测得电感L和电阻R的数据，通过器件手册查询晶闸管的导通电阻r，并且记录相应的控制触发角α的数据；
步骤二，建立矩阵：
分别对直流侧电压、直流侧电流以及交流侧电流进行傅里叶变换和矩阵分解；
步骤二的具体步骤如下：
第一步，将整流桥中晶闸管控制信号，即开关函数分解为傅里叶级数；
第一步中，考虑换相过程后，开关函数的傅里叶级数为：
Figure FDA0003846883510000011
Figure FDA0003846883510000012
上式中，ω代表电网频率，约为314.15rad/s,α为预设的开关函数移相角，μ为换相重叠角，可根据移相角α，电路形式及网侧变压器漏感推算；
由于三相平衡，通过延时1/3和2/3个工频周期T，可以得到B相和C相的开关函数sub(t)和suc(t)；
第二步，计算直流侧电压；
第二步中，直流侧电压的时域计算表达式如下：
udc(t)＝ua(t)sua(t)+ub(t)sub(t)+uc(t)suc(t)(10)
其中，ua(t)，ub(t)，uc(t)分别为三相电压的模型假定时域值，无需测量，其大小仅和电压幅值有关；
其中，sua(t)，sub(t)，suc(t)分别为三相电压开关函数；
由于三相平衡，直流侧电压只含有直流分量和6倍频分量，在这些频率上，上式中三部分结果相同，可以先求A相对应的结果再乘3，而其他频率分量三相结果之和为0，这里具体推导证明省略；
在MATLAB中定义四个零矩阵
Figure FDA0003846883510000021
和
Figure FDA0003846883510000022
其中Cdc1为交流电压的各次分量到直流电压直流分量的变换矩阵，Cdc2为交流电压的各次分量的共轭到直流电压直流分量的变换矩阵，C1为交流电压的各次分量到直流电压各次分量的变换矩阵；利用MATLAB程序对A相电压的各次分量与开关函数各次分量的相乘情况进行判断，只保留结果中直流和50次以下6倍频的分量，并乘以3，将每种情况对应的变换系数叠加到矩阵Cdc1，Cdc2，C1，C2的对应位置元素上；
最终直流侧电压频域表达式为：
Figure FDA0003846883510000023
第三步，计算直流侧电流；
第三步的计算方法是，在正常导通过程中，电流流经两个晶闸管，存在两个导通电阻，而在换相过程中，由于3个晶闸管同时导通，相当于流经1.5个导通电阻，这里忽略换相时导通电阻的变化，认为一直存在两个导通电阻r，将其归到直流侧，则直流侧电阻包含负载电阻R和两个晶闸管导通电阻r，设直流侧总电感为L，则直流侧电流的直流分量以及50次以内各次谐波分量为：
Figure FDA0003846883510000024
式中
Figure FDA0003846883510000025
第四步，计算交流侧电流；
第四步中，交流侧电流的时域表达式如下：
Figure FDA0003846883510000026
其中，idc(t)为直流侧电流，ia(t)，ib(t)，ic(t)分别为三相电流，sia(t)，sib(t)，sic(t)分别为三相电流开关函数；
对于三相平衡情况，只需计算A相，其他两相通过电流的正负序关系可以求得， 在MATLAB中定义两个零矩阵
Figure FDA0003846883510000031
和
Figure FDA0003846883510000032
其中C3为直流侧电流的各次分量到交流侧电流各次分量的变换矩阵，C4直流侧电流的各次分量的共轭到交流侧电流各次分量的变换矩阵， 利用MATLAB程序对每个电流谐波分量与开关函数谐波分量相乘的情况进行判断，只保留结果中频率在50次谐波以下的情况，将每种情况对应的变换系数叠加到C3，C4的对应元素上，则最终交流侧电流频域表达式为：
Figure FDA0003846883510000033
对于整流器的交流侧谐波电压和电流，只有典型次非零，所以如果只关注典型谐波次，式(15)可以写为：
Figure FDA0003846883510000034
步骤三，预测与补偿：
根据步骤二中矩阵得出的谐波耦合外特性预测模型，预测此时整流电路电压电流的谐波含量，并依据国标《GBT 14549—93电能质量公用电网谐波标准》判断谐波含量是否超标，当模型预测得到的电压电流谐波含量大于标准中的限值，会开启有源电力滤波器设备对电路进行谐波补偿，低于限值时有源电力滤波器将处于脉冲封锁，即关机状态。
2.根据权利要求1所述的一种三相整流电路的谐波耦合外特性预测方法，其特征在于，步骤一的具体步骤是，通过阻抗测量仪Bode-100测得电感L和电阻R的数据，而其中控制触发角α的数据是设定值，通过对晶闸管控制触发程序对其设定而产生。
CN202010422150.7A 2020-05-18 2020-05-18 一种三相整流电路的谐波耦合外特性预测方法 Active CN111682542B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010422150.7A CN111682542B (zh) 2020-05-18 2020-05-18 一种三相整流电路的谐波耦合外特性预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010422150.7A CN111682542B (zh) 2020-05-18 2020-05-18 一种三相整流电路的谐波耦合外特性预测方法

Publications (2)

Publication Number Publication Date
CN111682542A CN111682542A (zh) 2020-09-18
CN111682542B true CN111682542B (zh) 2022-12-09

Family

ID=72434184

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010422150.7A Active CN111682542B (zh) 2020-05-18 2020-05-18 一种三相整流电路的谐波耦合外特性预测方法

Country Status (1)

Country Link
CN (1) CN111682542B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113848423B (zh) * 2021-09-18 2023-05-05 广东电网有限责任公司 非线性结点检测方法及系统
CN113688538B (zh) * 2021-09-27 2023-10-27 福州大学 采用预测修正法辨识设备参数的谐波责任全景量化方法
CN115036925A (zh) * 2022-04-02 2022-09-09 深圳市宝安任达电器实业有限公司 一种用于变换器非线性负载的谐波分析导纳建模方法
CN116247952A (zh) * 2023-05-08 2023-06-09 国网山西省电力公司临汾供电公司 一种三桥臂功率可控的三相桥式整流变换器

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104135038A (zh) * 2014-07-24 2014-11-05 华北电力大学 一种交直流混联系统不对称故障分析方法
KR101535854B1 (ko) * 2014-04-23 2015-07-27 한국생산기술연구원 병렬처리 알고리즘을 이용한 이산 푸리에 변환 방법
CN105244882A (zh) * 2015-11-13 2016-01-13 湖南大学 三相三线制有源电力滤波器对电网谐波电流的补偿方法
CN109657338A (zh) * 2018-12-17 2019-04-19 大连理工大学 一种新颖的三相电压型变流器的逻辑切换建模方法
CN109995034A (zh) * 2019-01-28 2019-07-09 西安交通大学 一种考虑谐波耦合特性的区域性电网谐波综合治理的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101535854B1 (ko) * 2014-04-23 2015-07-27 한국생산기술연구원 병렬처리 알고리즘을 이용한 이산 푸리에 변환 방법
CN104135038A (zh) * 2014-07-24 2014-11-05 华北电力大学 一种交直流混联系统不对称故障分析方法
CN105244882A (zh) * 2015-11-13 2016-01-13 湖南大学 三相三线制有源电力滤波器对电网谐波电流的补偿方法
CN109657338A (zh) * 2018-12-17 2019-04-19 大连理工大学 一种新颖的三相电压型变流器的逻辑切换建模方法
CN109995034A (zh) * 2019-01-28 2019-07-09 西安交通大学 一种考虑谐波耦合特性的区域性电网谐波综合治理的方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Analytical Derivation of a Coupled-Circuit Model of a Claw-Pole Alternator With Concentrated Stator Windings;Hua Bai;《IEEE TRANSACTIONS ON ENERGY CONVERSION》;20020331;全文 *
Modeling Issues in Three-Phase Current Source Rectifiers that use Damping Resistors;Pedro E. Melin;《IEEE ISIE 2006》;20060712;全文 *
三相整流负荷谐波频域模型及谐波功率分析;魏能峤;《中国优秀硕士学位论文全文数据库》;20190415;全文 *
有源配电网谐波特性分析及应对策略研究;阚沁怡;《中国优秀硕士学位论文全文数据库》;20190515;全文 *
高压配电柜故障监视和自动控制系统设计;朱瑞卿;《煤矿机械》;20120531;全文 *

Also Published As

Publication number Publication date
CN111682542A (zh) 2020-09-18

Similar Documents

Publication Publication Date Title
CN111682542B (zh) 一种三相整流电路的谐波耦合外特性预测方法
Forghani et al. Online wavelet transform-based control strategy for UPQC control system
Sun et al. Small-signal input impedance modeling of line-frequency rectifiers
Basic Input current interharmonics of variable-speed drives due to motor current imbalance
Soltani et al. Effects of passive components on the input current interharmonics of adjustable-speed drives
CN203870168U (zh) 有源电力滤波器中电抗器电感量的测量装置
Petković et al. Cascaded H-bridge multilevel converter for a high-power medium-voltage impedance-admittance measurement unit
CN111030156B (zh) Lcc-hvdc谐波等效电路建模方法及装置
Yaghoobi et al. Current harmonic estimation techniques based on voltage measurements in distribution networks
CN109617076B (zh) 用于改善pcc处电能质量的电压检测型有源滤波器
CN103941103A (zh) 有源电力滤波器中电抗器电感量的测量装置及方法
Karimi et al. Input impedance modeling of line-frequency rectifiers by the method of impedance mapping
CN116087623A (zh) 一种新能源并网系统整体阻抗测量方法及装置
CN115236404A (zh) 一种并网逆变器端口阻抗自测量方法
Popa et al. Analysis and simulation of line notching attenuation for a static excitation system in a power system
Cano et al. New transfer functions for an accurate estimation of harmonic distortion in AC/DC converters working under unbalanced conditions
CN108347055B (zh) 一种并网滤波电感器参数评测电路及其控制方法
Xu et al. Study on equivalence of MMC modeling under multi-dq frames and harmonic state space
Yildiz et al. Harmonic modeling of full-wave diode rectifier for nonuniform load currents
Tousi et al. A New Topology for Dynamic Voltage Restorer Using High Frequency Link.
Xu et al. Active filter implementation using a generalized nonactive power theory
Lodhi et al. Analysis of Harmonics Detection and their Filtration by using Active Band Pass Filter for Domestic Appliances
Yang et al. Impedance modelling and stability assessment for grid‐connected dual‐PWM adjustable speed drives with motor torque oscillations
Ivanov et al. Investigation of Passive LC Filter in the Induction Heating System with High Frquency Power Supply
Moreau et al. Equivalent EMI noise source modelling in power converters

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