CN108051682B - 一种单相整流器系统阻抗模型的验证方法 - Google Patents
一种单相整流器系统阻抗模型的验证方法 Download PDFInfo
- Publication number
- CN108051682B CN108051682B CN201711400305.1A CN201711400305A CN108051682B CN 108051682 B CN108051682 B CN 108051682B CN 201711400305 A CN201711400305 A CN 201711400305A CN 108051682 B CN108051682 B CN 108051682B
- Authority
- CN
- China
- Prior art keywords
- voltage
- current
- frequency
- coefficient
- frequency spectrum
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
Abstract
本发明公开一种单相整流器系统阻抗模型的验证方法，包括以下步骤：步骤1：在单相整流器系统中注入扰动，测量交流侧响应的电压v s 和电流i s ；步骤2：根据步骤1得到的电压v s 和电流i s 通过快速傅里叶变换方法提取频率在|f p±f 1|、|f p±3f 1|下的幅值和相位；并计算对应频率下的电压频谱系数和电流频谱系数；f p 为扰动频率，f 1 为基波频率；步骤3：修正步骤2得到的电压频谱系数和电流频谱系数；步骤4：将步骤3修正后的结果带入阻抗模型表达式，计算误差结果；步骤5：根据误差结果判据，判断模型正确性；本发明适用于所有单相整流器系统阻抗模型，修正对交流电压、电流负频率频谱系数提取的缺陷进行修正，结算结果准确。
Description
技术领域
本发明涉及自动控制领域，具体涉及一种单相整流器系统阻抗模型的验证方法。
背景技术
随着对能量产生、传输和利用的更高要求，更多的电力电子转换器由于其高效率和灵活的控制而在当今的电力系统中投入使用；考虑到每个转换器的复杂控制，基于阻抗的稳定性分析方法被广泛采用；整个系统的稳定性可以通过分析级联或并联的每个子系统的外部阻抗或导纳来分析；多重谐波线性化方法和谐波传递函数矩阵方法的应用与分析中，揭示了整流器系统为MIMO系统，其阻抗矩阵为高阶矩阵；Shahil Shah在单相逆变器阻抗建模中发现如果在交流侧考虑三次谐波分量模型更加准确；考虑三次谐波的影响，阻抗矩阵就变成4阶矩阵；目前验证方法中不同的单相整流器系统需要不同的验证方法，比较复杂；还没有通过注入扰动来直接计算矩阵的方法；并且如果通过注入扰动来直接计算会存在计算量大、扰动注入电路多且复杂的问题。
发明内容
本发明提供一种只需要注入一个扰动，并适用于所有单相整流器系统阻抗模型的单相整流器系统阻抗模型的验证方法。
本发明采用的技术方案是：一种单相整流器系统阻抗模型的验证方法，包括以下步骤：
步骤1：在单相整流器系统中注入扰动，测量交流侧响应的电压vs和电流is；
步骤2：根据步骤1得到的电压vs和电流is通过快速傅里叶变换方法提取频率在|fp±f1|、|fp±3f1|下的幅值和相位；并计算对应频率下的电压频谱系数和电流频谱系数；fp为扰动频率，f1为基波频率；
步骤3：修正步骤2得到的电压频谱系数和电流频谱系数；
步骤4：将步骤3修正后的结果带入阻抗模型表达式，计算误差结果；
步骤5：根据误差结果判据，判断模型正确性。
进一步的，所述步骤3中修正步骤2得到的电压频谱系数和电流频谱系数的过程如下：
判断fp-f1是否小于0，若是则对对应频率下的电压频谱系数和电流频谱系数取共轭作为新的电压频谱系数和电流频谱系数；若否则保留原计算结果；
判断fp-3f1是否小于0，若是则对对应频率下的电压频谱系数和电流频谱系数取共轭作为新的电压频谱系数和电流频谱系数；若否则保留原计算结果。
进一步的，所述步骤2的具体过程如下：
选择一个扰动角频率ωp，在交流电源侧注入一个电压扰动vptb；
其中：vptb＝Vptbcos(ωp+ω1)t；式中，ω1为基波角频率，Vptb为扰动电压源峰值，t为时间；交流侧响应电压为vs，电流为is；
根据下式获得频率在|fp±f1|、|fp±3f1|下的幅值和相位：
式中：x表示电压或电流；N为采样点的数量；WN为蝶形因子，n表示第n个信号，k为第k个离散时间点，X(k)为第k个离散时间点提取的信号；
根据得到的幅值和相位计算相应频率下的电压频谱系数Vn[fp+f1]、Vn[fp-f1]、Vn[fp+3f1]、Vn[fp-3f1]，电流频谱系数In[fp+f1]、In[fp-f1]、In[fp+3f1]、In[fp-3f1]。
进一步的，所述步骤4中计算方法如下：
其中：Z4×4为单相整流器系统阻抗矩阵，σ1、σ2、σ3、σ4为误差计算结果。
本发明的有益效果是：
(1)本发明为针对单相交流器系统高阶阻抗模型的验证方法，其测量方法简单、并且每个扰动频率只需要注入一个扰动；
(2)本发明对频谱系数进行修正，可修正对交流电压、电流负频率频谱系数提取的缺陷进行修正，结算结果准确。
附图说明
图1为本发明电流扰动注入示意图。
图2为本发明电压扰动注入示意图。
图3为本发明的实际扰动vptb注入到实际电路的示意图。
图4为本发明的PWM调制示意图。
图5为本发明的测量频谱系数处理方式图。
图6为本发明中fp＝10Hz测量结果复平面图。
图7为本发明fp＝60Hz测量结果复平面图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步说明。
本发明通过向电路中注入一个扰动，测得交流侧电压和电流；利用快速傅里叶变换FFT方法提取交流侧电压、电流在|fp±f1|、|fp±3f1|频率下的幅值和相位，并计算对应频率下的电压频谱系数和电流频谱系数；判断fp-3f1、fp-f1结果是否为负，如果为负则对其绝对值频率下的频谱系数取共轭来获得负频率对应的频谱系数，最后代入阻抗验证表达式矩阵求误差量σ1、σ2、σ3、σ4，并带入误差判据来判断理论阻抗模型是否正确。
具体过程如下：
步骤1：在单相整流器系统中注入扰动，测量交流侧响应的电压vs和电流is；
为了测量子系统的外部阻抗的值，需要注入扰动，如图1和图2所示；有两种类型的扰动，电压扰动和电流扰动，其与源的类型相关；在单相整流器系统中，注入电压扰动量为：
vptb＝Vptbcos(ωp+ω1)t
式中：ω1为基波角频率，Vptb为扰动电压源峰值，t为时间，ωp为扰动角频率，vptb为电压扰动。
步骤2：根据步骤1得到的电压vs和电流is通过快速傅里叶变换方法提取频率在|fp±f1|、|fp±3f1|下的幅值和相位；并计算对应频率下的电压频谱系数和电流频谱系数；fp为扰动频率，f1为基波频率；
根据下式获得频率在|fp±f1|、|fp±3f1|下的幅值和相位：
式中：x表示电压或电流；N为采样点的数量；WN为蝶形因子，n表示第n个信号，k为第k个离散时间点，X(k)为第k个离散时间点提取的信号；
根据得到的幅值和相位计算相应频率下的电压频谱系数Vn[fp+f1]、Vn[fp-f1]、Vn[fp+3f1]、Vn[fp-3f1]，电流频谱系数In[fp+f1]、In[fp-f1]、In[fp+3f1]、In[fp-3f1]。
步骤3：修正步骤2得到的电压频谱系数和电流频谱系数；
判断fp-f1是否小于0，若是则对对应频率下的电压频谱系数Vn[fp-f1]和电流频谱系数In[fp-f1]取共轭作为新的电压频谱系数和电流频谱系数；若否则保留原计算结果；
判断fp-3f1是否小于0，若是则对对应频率下的电压频谱系数Vn[fp-3f1]和电流频谱系数In[fp-3f1]取共轭作为新的电压频谱系数和电流频谱系数；若否则保留原计算结果。
因为电压和电流在|fp±f1|、|fp±3f1|频率下的频谱系数都是正频率的频谱系数，而实际系统可能存在负频率，因此FFT的方法存在缺陷，所以需要进行修正。
步骤4：将步骤3修正后的结果带入阻抗模型表达式，计算误差结果
其中：Z4×4为单相整流器系统阻抗矩阵，σ1、σ2、σ3、σ4为误差计算结果。
步骤5：根据误差结果判据，判断模型正确性
判断|σi|≤ε是否成立，如果成立则阻抗模型准确，反之阻抗模型不准确，其中，ε为精度。
下面通过具体实施例对本发明的有效性和正确性进行验证。
在Matlab/Simulink搭建一个简单的模型，即开环控制的电压整流器如图3所示；该整流器为CRH3型动车组四象限整流器电路，来验证该方法的有效性和正确性；该电路使用了电压源，所以选择图2的扰动方式；整流器采用开环控制策略，正弦信号m设置峰值为M1初始相位为为θm1，经过PWM调制输出4路信号到4个IGBT，如图4所示；控制直流侧电压vd稳定在3000V；vs为交流电压源，幅值为V1，相位为0，频率为f1＝50Hz；Rs、Ls为电压源vs的等效电阻、电抗；Rn、Ln为牵引变压器的等效电阻、电抗；Cd为输出侧电容、Rd为逆变器的等效电阻；开关频率fsw为350Hz；采样频率fs为20kHz；公共耦合点(point of common coupling,PCC)电压为vn。
电路而定参数如表1所示：
表1单相整流器参数
通过推导，CRH3型车的单个整流器输入阻抗表达式为：
其中：
Z41＝Z42＝0；
对于电压vn、电流in的相应频率下经过FFT计算频谱系数的示意图如图5；将提取计算的频谱系数代入阻抗表达式验证测量结果：
通过频率扫描法，fp取值以5Hz为步长，计算从fp等于5Hz到fp等于100Hz(其中fp≠50、100Hz)扰动频率注入下计算σ1、σ2、σ3、σ4的模长如表2；精度ε＝0.3；分别对fp＝10Hz和fp＝60Hz的两种扰动电压源注入电路后，计算的阻抗表达式的误差结果σ1、σ2、σ3、σ4绘制在复平面图中如图6、7所示，这些点分布在以原点为中心，0.3为半径的圆中，因此|σ1|、|σ2|、|σ3|、|σ4|均小于0.3，说明阻抗模型正确。
表2仿真计算结果
从表2数据看出，当注入扰动频率为fp+f1＝55Hz时，计算|σ3|为0.3096，其他扰动频率注入后计算的|σ1|、|σ2|、|σ3|、|σ4|均小于0.3，可见模型建立正确，测量方法准确。
本发明方法适用于测量所有单相整流器系统阻抗模型的验证，并解决了信号负频率求取频谱系数的问题。
Claims (1)
1.一种单相整流器系统阻抗模型的验证方法，其特征在于，包括以下步骤：
步骤1：在单相整流器系统中注入扰动，测量交流侧响应的电压vs和电流is；
步骤2：根据步骤1得到的电压vs和电流is通过快速傅里叶变换方法提取频率在|fp±f1|、|fp±3f1|下的幅值和相位；并计算对应频率下的电压频谱系数和电流频谱系数；fp为扰动频率，f1为基波频率；
步骤3：修正步骤2得到的电压频谱系数和电流频谱系数；
步骤4：将步骤3修正后的结果带入阻抗模型表达式，计算误差结果；
步骤5：根据误差结果判据，判断模型正确性；
所述步骤2的具体过程如下：
选择一个扰动角频率ωp，在交流电源侧注入一个电压扰动vptb；
其中：vptb＝Vptbcos(ωp+ω1)t；式中，ω1为基波角频率，Vptb为扰动电压源峰值，t为时间；
交流侧响应电压为vs，电流为is；
根据下式获得频率在|fp±f1|、|fp±3f1|下的幅值和相位：
式中：x表示电压或电流；N为采样点的数量；WN为蝶形因子，n表示第n个信号，k为第k个离散时间点，X(k)为第k个离散时间点提取的信号；
根据得到的幅值和相位计算相应频率下的电压频谱系数Vn[fp+f1]、Vn[fp-f1]、Vn[fp+3f1]、Vn[fp-3f1]，电流频谱系数In[fp+f1]、In[fp-f1]、In[fp+3f1]、In[fp-3f1]；
所述步骤3中修正步骤2得到的电压频谱系数和电流频谱系数的过程如下：
判断fp-f1是否小于0，若是则对对应频率下的电压频谱系数和电流频谱系数取共轭作为新的电压频谱系数和电流频谱系数；若否则保留原计算结果；
判断fp-3f1是否小于0，若是则对对应频率下的电压频谱系数和电流频谱系数取共轭作为新的电压频谱系数和电流频谱系数；若否则保留原计算结果；
所述步骤4中计算方法如下：
其中：Z4×4为单相整流器系统阻抗矩阵，σ1、σ2、σ3、σ4为误差计算结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711400305.1A CN108051682B (zh) | 2017-12-22 | 2017-12-22 | 一种单相整流器系统阻抗模型的验证方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711400305.1A CN108051682B (zh) | 2017-12-22 | 2017-12-22 | 一种单相整流器系统阻抗模型的验证方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108051682A CN108051682A (zh) | 2018-05-18 |
| CN108051682B true CN108051682B (zh) | 2019-12-20 |
Family
ID=62131376
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711400305.1A Expired - Fee Related CN108051682B (zh) | 2017-12-22 | 2017-12-22 | 一种单相整流器系统阻抗模型的验证方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108051682B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110780118A (zh) * | 2019-10-15 | 2020-02-11 | 西南交通大学 | 一种应用于车网系统的多频率dq阻抗测量的计算方法 |
| CN110601268B (zh) * | 2019-10-29 | 2023-04-14 | 中国石油大学(华东) | 一种双馈风机并网端口输出阻抗建模及稳定性分析方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101957872A (zh) * | 2010-07-09 | 2011-01-26 | 南方电网技术研究中心 | 一种交直流电力系统的混合实时仿真方法 |
| CN104062500A (zh) * | 2014-07-04 | 2014-09-24 | 武汉大学 | 基于汉明乘积窗的信号谐波分析方法及系统 |
| CN104083169A (zh) * | 2014-08-01 | 2014-10-08 | 思澜科技（成都）有限公司 | 基于均匀阻抗模型的bis系统验证方法 |
| CN106908665A (zh) * | 2017-02-10 | 2017-06-30 | 西南交通大学 | 一种单相交流系统dq分解阻抗测量方法 |
| CN106936125A (zh) * | 2015-12-29 | 2017-07-07 | 中国电力科学研究院 | 一种广义二阶积分锁相环小信号阻抗建模方法 |
-
2017
- 2017-12-22 CN CN201711400305.1A patent/CN108051682B/zh not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101957872A (zh) * | 2010-07-09 | 2011-01-26 | 南方电网技术研究中心 | 一种交直流电力系统的混合实时仿真方法 |
| CN104062500A (zh) * | 2014-07-04 | 2014-09-24 | 武汉大学 | 基于汉明乘积窗的信号谐波分析方法及系统 |
| CN104083169A (zh) * | 2014-08-01 | 2014-10-08 | 思澜科技（成都）有限公司 | 基于均匀阻抗模型的bis系统验证方法 |
| CN106936125A (zh) * | 2015-12-29 | 2017-07-07 | 中国电力科学研究院 | 一种广义二阶积分锁相环小信号阻抗建模方法 |
| CN106908665A (zh) * | 2017-02-10 | 2017-06-30 | 西南交通大学 | 一种单相交流系统dq分解阻抗测量方法 |
| CN106908665B (zh) * | 2017-02-10 | 2019-03-29 | 西南交通大学 | 一种单相交流系统dq分解阻抗测量方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108051682A (zh) | 2018-05-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Asiminoaei et al. | A new method of on-line grid impedance estimation for PV inverter | |
| CN108037475B (zh) | 一种电容式电压互感器谐波传递特性试验系统 | |
| CN102928803B (zh) | 一种基于同步脉冲输出电源的电子式互感器校验装置 | |
| CN109030941A (zh) | Hanning自乘卷积窗FFT三谱线插值谐波分析方法 | |
| CN102914697A (zh) | 基于三相对称方波电流注入的微型电网谐波阻抗测量方法 | |
| CN108051682B (zh) | 一种单相整流器系统阻抗模型的验证方法 | |
| CN110763918B (zh) | 一种牵引脉冲整流器输入阻抗及牵引网输出阻抗测量方法 | |
| CN106908665A (zh) | 一种单相交流系统dq分解阻抗测量方法 | |
| CN106505840A (zh) | 一种并网光伏逆变器谐波治理方法 | |
| CN105259414A (zh) | 一种基于逆变器的电网阻抗在线检测方法 | |
| CN101566663A (zh) | 一种配电系统电压跌落源定位方法 | |
| CN107390022A (zh) | 基于离散频谱校正的电能计量方法 | |
| CN105510719A (zh) | 三相电网谐波阻抗测量方法 | |
| CN104852616A (zh) | 带有线路阻抗模拟功能的电网模拟器及控制方法 | |
| CN104852620A (zh) | 三相电压型pwm逆变器控制方法 | |
| CN210005677U (zh) | 考虑温度和电流变化的电能表基本误差测试装置 | |
| CN111884218A (zh) | 一种双馈入vsc输电系统稳定性评估方法及系统 | |
| CN202939299U (zh) | 一种基于同步脉冲输出电源的电子式互感器校验装置 | |
| CN104020350B (zh) | 一种克服频率摄动的电压基波分量检测方法 | |
| CN207472983U (zh) | 一种基于prbs扰动注入的电网阻抗在线辨识装置 | |
| CN109283391A (zh) | 一种基于非线性拟合的同步相量量测方法 | |
| CN102749488A (zh) | 电网谐波实时在线监测仪及利用该监测仪检测谐波的方法 | |
| CN110007145B (zh) | 一种基于恒压源变换器和电压电流相位差的谐振点检测方法 | |
| CN111308207A (zh) | 一种针对单相交流系统的dq阻抗测量方法 | |
| CN104393813B (zh) | 永磁同步电机的直轴电感测量方法 |
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 | ||
| 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: 20191220 Termination date: 20201222 |