CN113156207B - 一种牵引供电系统实时动态谐波检测方法 - Google Patents
一种牵引供电系统实时动态谐波检测方法 Download PDFInfo
- Publication number
- CN113156207B CN113156207B CN202110439525.5A CN202110439525A CN113156207B CN 113156207 B CN113156207 B CN 113156207B CN 202110439525 A CN202110439525 A CN 202110439525A CN 113156207 B CN113156207 B CN 113156207B
- Authority
- CN
- China
- Prior art keywords
- time
- real
- power supply
- supply system
- traction power
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明公开了一种牵引供电系统实时动态谐波检测方法。包括步骤:实时采集牵引供电系统电压u(t);将系统实时角频率ω与时间t相乘后,再与2π进行取模运算得到整数k;以k值对时间t进行划分,k值相同的时间t被划分为一个时间段,记为动态时间窗Td(k);以动态时间窗Td(k)作为u(t)的度量标准,将动态时间窗Td(k)内的牵引供电系统电压u(t)划分为一个片段,得到k个u(t)片段;将k个u(t)片段,分别使用快速傅里叶变换进行频谱提取,得到实时动态谐波。本发明的有益效果在于,能够应用于牵引供电系统运行频率偏差下的谐波检测,同时也适用于恒定基频运行下的谐波检测,应用面广。在检测过程中,还能够通过调整快速傅里叶变换FFT分析的片段数,调节谐波检测的频率分辨率。
Description
技术领域
本发明涉及牵引供电系统谐波检测技术领域,特别是一种牵引供电系统实时动态谐波检测方法。
背景技术
电气化铁路承担着我国绝大部分的交通负荷,保障其安全稳定运行意义重大。随着越来越多的大功率动车组及新型电力机车的接入,使得牵引供电系统谐波特性愈发复杂,其诱导的谐波谐振问题愈加突出,严重时将造成牵引闭锁、列车晚点、旅客滞留等问题。而解决系统谐波谐振问题极为重要的一步是精确识别牵引供电系统的谐波成分,从而掌握其谐波分布特性。然而实际的牵引供电系统运行频率难以维持在恒定的基频处,其运行频率偏差将导致传统的固定时间窗的FFT分析出现频谱泄露现象,使得检测出的谐波成分不准确。
发明内容
本发明的目的是提供一种牵引供电系统实时动态谐波检测方法,精确检测出牵引供电系统谐波特性。
实现本发明目的的技术方案如下:
一种牵引供电系统实时动态谐波检测方法,包括:
步骤1:实时采集牵引供电系统电压u(t);
步骤2:将系统实时角频率ω与时间t相乘后,再与2π进行取模运算得到整数 k;其中,系统实时角频率ω为:使用二阶广义积分器将u(t)转化为静止坐标系下的两个正交分量uα(t)及uβ(t);使用dq变换公式计算得到旋转坐标系下的d轴分量ud(t)与q轴分量uq(t),使用基于控制q轴分量uq(t)为0 的锁相技术得到系统实时角频率ω;
步骤3:以k值对时间t进行划分,k值相同的时间t被划分为一个时间段,记为动态时间窗Td(k);
步骤4:以动态时间窗Td(k)作为u(t)的度量标准,将动态时间窗Td(k)内的牵引供电系统电压u(t)划分为一个片段,得到k个u(t)片段;
步骤5:将k个u(t)片段,分别使用快速傅里叶变换进行频谱提取,得到实时动态谐波。
进一步的,所述步骤5替换为:将k个u(t)片段,按前后顺序依次组合得到 N个组合片段,分别对每个组合片段使用快速傅里叶变换进行频谱提取,得到实时动态谐波;其中,N的取值方法为:根据谐波检测要求的最低频谱分辨率fres,求解对应的最小牵引供电系统电压u(t)的采集时间Tmin,Tmin=1/fres;N取u(t)总时间大于Tmin的最小整数值。
本发明的有益效果在于,
一、能够应用于牵引供电系统运行频率偏差下的谐波检测,同时也适用于恒定基频运行下的谐波检测,应用面广。
二、在检测过程中,还能够通过调整快速傅里叶变换FFT分析的组合片段数,调节谐波检测的频率分辨率。
附图说明
图1为划分动态时间窗Td(k)示意图。
图2为牵引供电系统电压u(t)划分示意图。
具体实施方式
下面结合附图详细描述本发明的技术方案。
首先,利用电压传感器实时采集牵引供电系统电压u(t)。其次,使用二阶广义积分器(SOGI)将u(t)转化为静止坐标系下的两个正交分量uα(t)及uβ(t)。然后,利用dq变换公式计算得到旋转坐标系下的d轴分量ud(t) 与q轴分量uq(t)。其中ω为基于控制q轴分量uq(t)为0的锁相技术提供的系统实时角频率。继而,根据系统实时角频率ω实时划分动态时间窗Td(k),结合动态时间窗Td(k)将牵引供电系统电压u(t)划分为k个片段。最后,以这些片段为基本单元,利用快速傅里叶变换FFT对其进行频谱提取,能够避免牵引供电系统运行频率偏差导致的频谱泄露现象,实时动态地检测出其谐波成分。
根据系统实时角频率ω实时划分动态时间窗Td(k),其具体实现方法如图1 所示:将系统实时角频率ω与时间t相乘;然后将所得结果与2π进行取模运算得到整数部分k;以k值对时间t进行划分,其中k值相同的时间t被划分为一个片段,记为动态时间窗Td(k)。
结合动态时间窗Td(k)将牵引供电系统电压u(t)划分为k个片段,其具体实现方法如图2所示:用动态时间窗Td(k)作为牵引供电系统电压u(t)的度量标准,将处于同一动态时间窗Td(k)的牵引供电系统电压u(t)划分到同一片段里。
进一步的技术方案,将k个片段,按前后顺序依次组合得到N个组合片段,分别对每个组合片段使用快速傅里叶变换进行频谱提取,得到实时动态谐波。这里,快速傅里叶变换FFT分析的组合片段数选取方法是:
根据谐波检测要求的最低频谱分辨率fres,求解对应的最小牵引供电系统电压u(t)的采集时间Tmin,Tmin=1/fres;快速傅里叶变换FFT分析的组合片段数取u(t) 总时间大于最小牵引供电系统电压u(t)的采集时间Tmin的最小整数值。
Claims (2)
1.一种牵引供电系统实时动态谐波检测方法,其特征在于,包括:
步骤1:实时采集牵引供电系统电压u(t);
步骤2:将系统实时角频率ω与时间t相乘后,再与2π进行取模运算得到整数k;其中,系统实时角频率ω为:使用二阶广义积分器将u(t)转化为静止坐标系下的两个正交分量uα(t)及uβ(t);使用dq变换公式计算得到旋转坐标系下的d轴分量ud(t)与q轴分量uq(t),使用基于控制q轴分量uq(t)为0的锁相技术得到系统实时角频率ω;
步骤3:以k值对u(t)对应的时间段进行划分,k值对应的时间段被划分为时间窗口,记为动态时间窗Td(k);
步骤4:以动态时间窗Td(k)作为u(t)的度量标准,将动态时间窗Td(k)内的牵引供电系统电压u(t)划分为一个片段,得到k个u(t)片段;
步骤5:将k个u(t)片段,分别使用快速傅里叶变换进行频谱提取,得到实时动态谐波。
2.如权利要求1所述的一种牵引供电系统实时动态谐波检测方法,其特征在于,所述步骤5替换为:将k个u(t)片段,按前后顺序依次组合得到N个组合片段,分别对每个组合片段使用快速傅里叶变换进行频谱提取,得到实时动态谐波;其中,N的取值方法为:根据谐波检测要求的最低频谱分辨率fres,求解对应的最小牵引供电系统电压u(t)的采集时间Tmin,Tmin=1/fres;N取u(t)总时间大于Tmin的最小整数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110439525.5A CN113156207B (zh) | 2021-04-23 | 2021-04-23 | 一种牵引供电系统实时动态谐波检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110439525.5A CN113156207B (zh) | 2021-04-23 | 2021-04-23 | 一种牵引供电系统实时动态谐波检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113156207A CN113156207A (zh) | 2021-07-23 |
CN113156207B true CN113156207B (zh) | 2022-04-22 |
Family
ID=76869668
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110439525.5A Active CN113156207B (zh) | 2021-04-23 | 2021-04-23 | 一种牵引供电系统实时动态谐波检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113156207B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114895104B (zh) * | 2022-05-19 | 2023-07-07 | 西南交通大学 | 一种牵引供电系统并联谐波谐振频率识别方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103293378A (zh) * | 2013-05-03 | 2013-09-11 | 西南交通大学 | 一种牵引供电系统谐波谐振频率的辨识方法 |
CN107340427A (zh) * | 2016-05-02 | 2017-11-10 | 是德科技股份有限公司 | 用二进制相移键控测试信号测量出测量仪的幅值相位响应 |
JP2019039859A (ja) * | 2017-08-28 | 2019-03-14 | 沖電気工業株式会社 | 分析装置、分析装置の表示方法および表示プログラム |
-
2021
- 2021-04-23 CN CN202110439525.5A patent/CN113156207B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103293378A (zh) * | 2013-05-03 | 2013-09-11 | 西南交通大学 | 一种牵引供电系统谐波谐振频率的辨识方法 |
CN107340427A (zh) * | 2016-05-02 | 2017-11-10 | 是德科技股份有限公司 | 用二进制相移键控测试信号测量出测量仪的幅值相位响应 |
JP2019039859A (ja) * | 2017-08-28 | 2019-03-14 | 沖電気工業株式会社 | 分析装置、分析装置の表示方法および表示プログラム |
Non-Patent Citations (1)
Title |
---|
基于幂迭代的电力系统模态谐振快速求解方法;李朝阳等;《电网技术》;20161010;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113156207A (zh) | 2021-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103018555B (zh) | 一种高精度的电力参数软件同步采样方法 | |
CN101055299A (zh) | 一种冲击性负荷电能计量方法 | |
CN110061615A (zh) | 逆变器非线性特性的定子电流谐波补偿方法 | |
CN110412349B (zh) | 基于插值dft的同步相量数据次同步振荡参数辨识法 | |
CN107543962B (zh) | 主导间谐波频谱分布的计算方法 | |
CN102323477A (zh) | 一种电网间谐波测量方法及测量仪 | |
CN103983849B (zh) | 一种实时高精度的电力谐波分析方法 | |
CN113156207B (zh) | 一种牵引供电系统实时动态谐波检测方法 | |
CN102221639A (zh) | 正负序电流实时检测的方法 | |
CN103675447A (zh) | 一种电气化铁路的高精度实时谐波分析方法 | |
CN105425038A (zh) | 一种电力系统间谐波测量方法 | |
CN109507480A (zh) | 一种邻近基波/谐波的间谐波检测方法和装置 | |
CN110031674A (zh) | 动车组列车能耗计量系统及方法 | |
CN102323487B (zh) | 基于谐波分量的输电线路零序电容抗干扰测量方法 | |
CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
CN202102051U (zh) | 一种电网间谐波测量仪 | |
CN103607158A (zh) | 低速下基于pi调节器和锁相环原理的磁链估计器 | |
CN111800055B (zh) | 一种双凸极电机平均转矩确定方法和装置 | |
CN117169590B (zh) | 一种基于软件变采样率的电力谐波分析的方法和装置 | |
CN104020350B (zh) | 一种克服频率摄动的电压基波分量检测方法 | |
CN104393809A (zh) | 适用于scr静止变频器的抽水蓄能机组低速位置检测方法 | |
CN110320400A (zh) | 准同步采样和改进能量算子的电压闪变包络参数提取方法 | |
CN106487039A (zh) | 一种适用于畸变电网的锁相技术 | |
CN114112396B (zh) | 一种采用电流信号分析的转速波动准平稳工况下轴承故障诊断方法 | |
CN109307800A (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 |