CN102323479A - 列车谐波电流分析方法 - Google Patents
列车谐波电流分析方法 Download PDFInfo
- Publication number
- CN102323479A CN102323479A CN201110251248A CN201110251248A CN102323479A CN 102323479 A CN102323479 A CN 102323479A CN 201110251248 A CN201110251248 A CN 201110251248A CN 201110251248 A CN201110251248 A CN 201110251248A CN 102323479 A CN102323479 A CN 102323479A
- Authority
- CN
- China
- Prior art keywords
- data
- array
- current
- time
- electric current
- 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
Links
Landscapes
- Electric Propulsion And Braking For Vehicles (AREA)
Abstract
列车谐波电流分析方法,采用高采样率数据记录仪连续记录大量的列车电缆中时域电流数据,然后将这些海量的数据导入计算机内存中,按照一定的时间间隔对数据分段进行均方根计算,然后对每个20ms内的数据进行快速傅里叶转换,使用滤波算法,可以快速的获得时域中的制定频率分量的变化。
Description
技术领域
本发明涉及列车谐波电流分析方法。
背景技术
目前,对于列车供电系统中工频谐波电流信号采集和分析方面,通常情况通过电流传感器连接频谱分析仪,采集特定频段的频率信号,该方案缺点在于测试偶然性大,不能对谐波电流的变化进行连续的采集记录,也无法针对指定频率连续记录电流的变化。列车状态分为启动、牵引、惰行、常用制动、紧急制动、停车,获取谐波电流的时域变化情况,对列车各状态下的电流干扰分析有很大的帮助。
发明内容
本发明所要解决的技术问题是,克服现有技术的上述缺点,提供一种列车谐波电流分析方法。
为了解决以上技术问题,本发明提供的列车谐波电流分析方法,其特征在于具体步骤如下:
第一步、数据导入:
将列车谐波电流数据以数组的方式导入计算机内存中,时间数据流以时间数组X[]表示,与时间数据流对应的电流数据流以电流数组Y[]表示,其中,时间数据的精度为保留小数点后五位;
第二步、压缩数据:
电流数组Y[]中,以1000个电流数据为步长步进,对每2000个电流数据进行正弦波均方根运算获得压缩后的电流数据,重叠率为50%,从而获得压缩后的电流数组Y1[];时间数组X[]中,保留时间数据的小数点后两位,从而形成与压缩后的电流数组Y1[]对应的压缩后的的时间数组X1[];
第三步、带通滤波:
用巴特沃兹滤波算法对压缩后的电流数组Y1[]进行6阶带通滤波获得电流数组Y1'[],巴特沃兹滤波算法的中心频率取9.5kHz,带宽取300Hz;Y1'[]与X1[]对应;
第四步、数据解析:
滤波后的电流数组Y1'[],以20个电流数据为步长步进,每2000个电流数据进行快速傅里叶变换,形成频率数组Y2[],从频率Y2[]中提取9.5kHz频率的电流值,存入新的电流数组F[];时间数组X1[]按均值处理,获得与电流数组F[]对应的时间数组T[];
第五步、生成图表数据:
以时间数组T[]和电流数组F[]生成图表数据;
第六步、谐波电流分析完毕。
本发明第二步中采用数据压缩,将源数据大大压缩,在确保精度情况下,提高了运算效率。列车谐波电流测量中,对采样率要求较高(100k/s及以上)。列车从启动、牵引、惰行、常用制动、紧急制动至停车,所需时间10分钟以上,在此时间内采集到的电流数据至少10万个以上,时间单位微秒级。而列车线路中的谐波干扰主要来自牵引系统,其电流工作频率为50Hz,一个波形周期20ms。将采样数据从微秒级压缩至毫秒级不仅能够满足谐波分析的要求,而且可以极大地减少运算时间。所以需要对采样数据进行压缩处理,通过实验发现进行RMS运算处理后的数据,跟实际谐波电流最为接近,而且使用50%数据重叠率,能够有效地减小极短时间内电流突变(脉冲)对数据误差的影响。
本发明第四步对数据进行了解析,抽取所需数据。从时域电流数据转换成频域电流数据,需要一定时间范围内的数据进行快速傅里叶转换,列车线路中谐波电流为50Hz的整数倍,所以转换后的频域最小步频为50Hz,采用2000个数据进行快速傅里叶转换后的步频为50Hz。50Hz电流的周期为20ms,采用90%的数据重叠率,可以满足20ms内的时域数据转换为频域数据,提取指定频率处的电流值,结合20ms时间步进值,可以形成最终的指定频率处的时间频域电流幅值。
本方案采用高采样率数据记录仪连续记录大量的列车电缆中时域电流数据,然后将这些海量的数据导入计算机内存中,按照一定的时间间隔对数据分段进行均方根计算,然后对每个20ms内的数据进行快速傅里叶转换,使用滤波算法,可以快速的获得时域中的制定频率分量的变化。
附图说明
下面结合附图对本发明作进一步的说明。
图1为本发明列车谐波电流分析方法。
具体实施方式
本发明列车谐波电流分析方法流程如图1所示。
列车谐波电流分析方法,具体步骤如下:
第一步、数据导入:
将列车谐波电流数据以数组的方式导入计算机内存中,时间数据流以时间数组X[]表示,与时间数据流对应的电流数据流以电流数组Y[]表示,其中,时间数据的精度为保留小数点后五位;
第二步、压缩数据:
电流数组Y[]中,以1000个电流数据为步长步进,对每2000个电流数据进行正弦波均方根运算获得压缩后的电流数据,重叠率为50%,从而获得压缩后的电流数组Y1[];时间数组X[]中,保留时间数据的小数点后两位,从而形成与压缩后的电流数组Y1[]对应的压缩后的的时间数组X1[];
第三步、带通滤波:
用巴特沃兹滤波算法对压缩后的电流数组Y1[]进行6阶带通滤波获得电流数组Y1'[],巴特沃兹滤波算法的中心频率取9.5kHz,带宽取300Hz;Y1'[]与X1[]对应;
第四步、数据解析:
滤波后的电流数组Y1'[],以20个电流数据为步长步进,每2000个电流数据进行快速傅里叶变换,形成频率数组Y2[],从频率Y2[]中提取9.5kHz频率的电流值,存入新的电流数组F[];时间数组X1[]按均值处理,获得与电流数组F[]对应的时间数组T[];
第五步、生成图表数据:
以时间数组T[]和电流数组F[]生成图表数据;
第六步、谐波电流分析完毕。
本方案采用高采样率数据记录仪连续记录大量的列车电缆中时域电流数据,然后将这些海量的数据导入计算机内存中,按照一定的时间间隔对数据分段进行均方根计算,然后对每个20ms内的数据进行快速傅里叶转换,使用滤波算法,可以快速的获得时域中的制定频率分量的变化。
除上述实施例外,本发明还可以有其他实施方式。凡采用等同替换或等效变换形成的技术方案,均落在本发明要求的保护范围。
Claims (1)
1. 列车谐波电流分析方法,其特征在于具体步骤如下:
第一步、数据导入:
将列车谐波电流数据以数组的方式导入计算机内存中,时间数据流以时间数组X[]表示,与时间数据流对应的电流数据流以电流数组Y[]表示,其中,时间数据的精度为保留小数点后五位;
第二步、压缩数据:
电流数组Y[]中,以1000个电流数据为步长步进,对每2000个电流数据进行正弦波均方根运算获得压缩后的电流数据,重叠率为50%,从而获得压缩后的电流数组Y1[];时间数组X[]中,保留时间数据的小数点后两位,从而形成与压缩后的电流数组Y1[]对应的压缩后的的时间数组X1[];
第三步、带通滤波:
用巴特沃兹滤波算法对压缩后的电流数组Y1[]进行6阶带通滤波获得电流数组Y1'[],巴特沃兹滤波算法的中心频率取9.5kHz,带宽取300Hz;Y1'[]与X1[]对应;
第四步、数据解析:
滤波后的电流数组Y1'[],以20个电流数据为步长步进,每2000个电流数据进行快速傅里叶变换,形成频率数组Y2[],从频率Y2[]中提取9.5kHz频率的电流值,存入新的电流数组F[];时间数组X1[]按均值处理,获得与电流数组F[]对应的时间数组T[];
第五步、生成图表数据:
以时间数组T[]和电流数组F[]生成图表数据;
第六步、谐波电流分析完毕。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110251248.1A CN102323479B (zh) | 2011-08-29 | 2011-08-29 | 列车谐波电流分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110251248.1A CN102323479B (zh) | 2011-08-29 | 2011-08-29 | 列车谐波电流分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102323479A true CN102323479A (zh) | 2012-01-18 |
CN102323479B CN102323479B (zh) | 2014-02-12 |
Family
ID=45451262
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110251248.1A Active CN102323479B (zh) | 2011-08-29 | 2011-08-29 | 列车谐波电流分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102323479B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102866294A (zh) * | 2012-09-20 | 2013-01-09 | 南车南京浦镇车辆有限公司 | 基于轨道电路的信号兼容性测试方法 |
CN104569651A (zh) * | 2014-11-24 | 2015-04-29 | 南车青岛四方机车车辆股份有限公司 | 动车组高次谐波发射测试系统 |
CN107064633A (zh) * | 2017-03-29 | 2017-08-18 | 广西电网有限责任公司电力科学研究院 | 城市轨道交通负荷谐波电流迭加系数确定方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1046982A (zh) * | 1989-05-30 | 1990-11-14 | 山东省煤炭科学研究所 | 电力系统谐波测量方法和测量仪 |
WO2009120765A1 (en) * | 2008-03-25 | 2009-10-01 | Abb Research Ltd. | Method and apparatus for analyzing waveform signals of a power system |
CN101937020A (zh) * | 2010-08-10 | 2011-01-05 | 南京世都科技有限公司 | 基于人工神经网络的电网谐波检测装置 |
CN101976820A (zh) * | 2010-08-31 | 2011-02-16 | 南京南瑞继保电气有限公司 | 变频电动机保护采样信号处理方法 |
-
2011
- 2011-08-29 CN CN201110251248.1A patent/CN102323479B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1046982A (zh) * | 1989-05-30 | 1990-11-14 | 山东省煤炭科学研究所 | 电力系统谐波测量方法和测量仪 |
WO2009120765A1 (en) * | 2008-03-25 | 2009-10-01 | Abb Research Ltd. | Method and apparatus for analyzing waveform signals of a power system |
CN101937020A (zh) * | 2010-08-10 | 2011-01-05 | 南京世都科技有限公司 | 基于人工神经网络的电网谐波检测装置 |
CN101976820A (zh) * | 2010-08-31 | 2011-02-16 | 南京南瑞继保电气有限公司 | 变频电动机保护采样信号处理方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102866294A (zh) * | 2012-09-20 | 2013-01-09 | 南车南京浦镇车辆有限公司 | 基于轨道电路的信号兼容性测试方法 |
CN104569651A (zh) * | 2014-11-24 | 2015-04-29 | 南车青岛四方机车车辆股份有限公司 | 动车组高次谐波发射测试系统 |
CN104569651B (zh) * | 2014-11-24 | 2017-04-05 | 中车青岛四方机车车辆股份有限公司 | 动车组高次谐波发射测试系统 |
CN107064633A (zh) * | 2017-03-29 | 2017-08-18 | 广西电网有限责任公司电力科学研究院 | 城市轨道交通负荷谐波电流迭加系数确定方法 |
CN107064633B (zh) * | 2017-03-29 | 2019-10-18 | 广西电网有限责任公司电力科学研究院 | 城市轨道交通负荷谐波电流迭加系数确定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102323479B (zh) | 2014-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102866027B (zh) | 基于lmd和局域时频熵的旋转机械故障特征提取方法 | |
CN103713237B (zh) | 一种电力系统输电线路短路故障诊断方法 | |
CN102721464B (zh) | 电力变压器绕组变形故障检测方法及系统 | |
CN106226635A (zh) | 一种配电网馈线故障类型识别方法及装置 | |
CN104136896A (zh) | 用于获取旋转马达的振动信号的系统 | |
CN104545887A (zh) | 伪差心电波形识别方法和装置 | |
CN103454520A (zh) | 基于在线频率响应法的变压器绕组变形在线监测方法 | |
CN108444589A (zh) | 一种基于频域特征提取的水电机组状态监测信号处理方法 | |
CN102323479B (zh) | 列车谐波电流分析方法 | |
CN110597144A (zh) | 一种电机传感器数据采集单元及数据采集方法 | |
CN103675544A (zh) | 基于优化算法的电力系统故障信号检测与波形识别方法 | |
CN104237727A (zh) | 一种变压器近区短路信号监测装置及短路录波分析方法 | |
CN102721885A (zh) | 一种基于小波分析的孤岛效应检测方法 | |
CN102023262B (zh) | 35kV电网弧光接地过电压识别方法 | |
CN104142423A (zh) | 1000kV高压并联电抗器的监测系统和方法 | |
CN202583342U (zh) | 高速列车设备状态及电能质量监测系统 | |
CN108090270B (zh) | 一种基于形态学滤波和盲源分离的暂态振荡参数识别方法 | |
CN107292067B (zh) | 一种基于压缩感知与双谱分析的齿轮故障诊断方法 | |
CN112861328B (zh) | 一种基于随机响应信号的发电机阻尼评估装置及方法 | |
CN111143927B (zh) | 一种基于结构响应线性组合的约束模态分解与频率识别方法 | |
CN106018922B (zh) | 电网过电压暂态分量时域波形提取方法及系统 | |
CN202614831U (zh) | 一种暂态过电压在线监测系统 | |
CN111693829A (zh) | 非接触式超声波检测的局部放电噪声与放电区分的方法 | |
CN110994802A (zh) | 一种输电导线运行状态监测方法和装置 | |
CN102928691A (zh) | 检测电气设备电磁瞬态骚扰信号的系统、方法及其装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C56 | Change in the name or address of the patentee | ||
CP01 | Change in the name or title of a patent holder |
Address after: 210031 Nanjing, Pu Pu North Road, No. 68, Jiangsu Patentee after: CRRC NANJING PUZHEN CO., LTD. Address before: 210031 Nanjing, Pu Pu North Road, No. 68, Jiangsu Patentee before: CSR NANJING PUZHEN Vehicles Co., Ltd. |