CN113190798A - 一种牵引变压器低压侧到高压侧的谐波传输系数估计方法 - Google Patents
一种牵引变压器低压侧到高压侧的谐波传输系数估计方法 Download PDFInfo
- Publication number
- CN113190798A CN113190798A CN202110462978.XA CN202110462978A CN113190798A CN 113190798 A CN113190798 A CN 113190798A CN 202110462978 A CN202110462978 A CN 202110462978A CN 113190798 A CN113190798 A CN 113190798A
- Authority
- CN
- China
- Prior art keywords
- voltage
- phase
- voltage side
- low
- harmonic
- 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.)
- Pending
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000003137 locomotive effect Effects 0.000 claims abstract description 13
- 238000005070 sampling Methods 0.000 claims description 25
- 230000001174 ascending effect Effects 0.000 claims description 21
- 230000000694 effects Effects 0.000 claims description 9
- 230000002776 aggregation Effects 0.000 claims description 8
- 238000004220 aggregation Methods 0.000 claims description 8
- 238000012417 linear regression Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 3
- 230000001360 synchronised effect Effects 0.000 abstract description 3
- 238000004445 quantitative analysis Methods 0.000 abstract description 2
- 238000004458 analytical method Methods 0.000 description 9
- 230000003321 amplification Effects 0.000 description 8
- 238000003199 nucleic acid amplification method Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 3
- 238000004880 explosion Methods 0.000 description 2
- 230000006855 networking Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 description 1
- 230000004931 aggregating effect Effects 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012372 quality testing Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000008929 regeneration Effects 0.000 description 1
- 238000011069 regeneration method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Discrete Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种牵引变压器低压侧到高压侧的谐波传输系数估计方法，利用电压同步测量设备记录牵引变压器110kV/220kV高压侧的A相、B相、C相以及牵引变压器27.5kV低压侧上行接触线Tα、上行正馈线Fα、下行接触线Tβ、下行正馈线Fβ的电压波形，来估计由机车发射谐波从牵引变压器27.5kV低压侧上行接触线Tα、上行正馈线Fα、下行接触线Tβ、下行正馈线Fβ到110kV/220kV高压侧A相、B相、C相的传输系数。与现有技术相比，本发明(1)不需要建立牵引供电系统谐波模型；(2)基于实测数据而不是仿真模型，更能反映实际系统的谐波传输特性；(3)定量的估计方法，可以给出较精确的定量的分析结果，而不是定性或粗略的估计。
Description
技术领域
本发明涉及牵引供电系统电能质量分析技术领域，特别是一种牵引变压器低压侧到高压侧的谐波传输系数估计方法。
背景技术
截至2020年底，中国高铁营业里程达到3.79万公里以上，超过世界高铁总里程的三分之二。如此大功率、工况多变(牵引、惰行、再生工况)、高速移动、高行车密度、网络化的牵引供电系统将使得区域电网的供电压力、功率交换激增。牵引负荷(动车组)具有波动剧烈、强非线性、不对称等特点，同时频繁的大功率交换使得区域电网与牵引供电系统的不安全因素增加。尤其在较为薄弱的区域电网下，将严重影响牵引供电系统的供电能力和区域电网的稳定性。因此，如何保障具有频繁冲击、大功率、高密度、网络化等特点的牵引供电系统及接入电网的安全、可靠、高效运行，将面临着新的且急迫的挑战。其中牵引网的谐波放大、谐波谐振问题在较多的文献中得到了报道，常引发避雷器爆炸、容性设备如电压互感器、电容器组的爆炸。同时，关于高速铁路牵引供电系统谐波含量及其与区域电网之间的谐波相互影响的研究也越来越成为谐波责任划分研究的热点之一。
现有分析牵引变压器两侧谐波电压传输系数的方法包括：
(1)仿真分析：目前针对牵引供电系统谐波传输及放大的研究主要通过建立较为精确的“区域电网—牵引变压器—牵引网—AT所—分区所”的完整基波/谐波模型，并建立动车组的谐波诺顿模型进行分析。分析方法包括谐波传输放大法、频谱分析法、谐振模态分析法、S域或频域传递函数法。以上方法可针对建立的系统谐波分析模型，对其中的谐波传输及放大规律进行较好地分析。然而，以上方法均存在以下问题：1)分析结果的准确性依赖于谐波建模的精确性，往往很难反映真实系统的谐波传输及放大特性；2)精确的谐波建模过程较为复杂；3)频谱分析法、谐振模态分析法、S域或频域传递函数法等只能定性判断是否发生放大，无法定量给出牵引变压器两侧谐波电压传输系数。
(2)基于实测数据的分析：目前只有基于谐波95％概率大值的对比分析，得出定性(是否发生放大，放大是指传输系数大于1)或粗略的量化结果。
发明内容
本发明的目的是提供一种牵引变压器低压侧到高压侧的谐波传输系数估计方法。
实现本发明目的的技术方案为：
一种牵引变压器低压侧到高压侧的谐波传输系数估计方法，包括：
步骤1：同步获取牵引变压器的电压波形和电流波形：
1.1以大于r*H*50Hz的采样率获取高压侧A相电压波形UA和低压侧上行接触线Tα电压波形UTα；其中，H为要分析的最高谐波次数，r为最高谐波次数的倍数；
1.2获取低压侧上行接触线Tα电流波形ITα和下行接触线Tβ电流波形ITβ；
步骤2：分别对UA和UTα进行快速傅里叶变换，得到其基波电压方均根值时间序列和谐波电压方均根值时间序列，进而得到h次谐波电压含有率时间序列，即
和
分别对ITα和ITβ进行快速傅里叶变换得到其基波电流方均根值时间序列，即
和
步骤3：对任一h次谐波，选择其谐波电压含有率时间段进行分割：
3.1选择
且
的时间段所对应的
和
的谐波电压含有率时间段；其中，δ1和δ2为有且只有上行区间有机车经过而下行区间无机车经过的判据；
3.2将谐波电压含有率时间段用滑动时间窗分割为宽度为L的N个子时间段；
步骤4：对N个子时间段，分别使用二元线性回归方程拟合得到每个子时间段中h次谐波电压含有率之间的关系yi＝λixi+bi+εi；其中，i为子时间段的序号，i＝1,…,N；yi为A相的h次谐波电压含有率，xi为上行接触线Tα的h次谐波电压含有率，λi为回归直线的斜率即h次谐波电压含有率从Tα到A相的传输系数，bi为截距，εi为回归残差；
步骤5：从拟合得到的N个关系中，选择具有聚集效应的K个关系，具体为：
5.1记N个子时间段中，第i个子时间段的
为Vi，第i个子时间段的
为Wi；其中，
为A相第i个子时间段中h次谐波电压含有率的方差，
为上行接触线Tα第i个子时间段中h次谐波电压含有率的方差；
5.2将N个Wi从大到小排序；令m％为需选择的关系的数量的初始百分比；
5.3计算m％*N，向下取整得到K；
5.4计算Var1与Var2的比值
其中，Var1为前K个Wi对应的K个Vi的方差，Var2为N个Vi的方差；
5.5若
小于阈值Δ，则判定前K个Wi对应的K个关系具有聚集效应；否则，令m％＝m％-1％返回步骤5.3重新计算，直到满足
小于阈值Δ；
步骤6：计算Var3与Var4的比值
其中，Var3为所选择的K个关系的λi的方差，Var4为N个关系的λi的方差；若ρ％小于阈值Δ，则取K个关系的λi的平均值，得到低压侧上行接触线Tα到高压侧A相的h次谐波传输系数。
进一步的技术方案，
将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧B相电压波形UB和低压侧上行接触线Tα电压波形UTα，其它步骤对应替换，得到低压侧上行接触线Tα到高压侧B相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧C相电压波形UC和低压侧上行接触线Tα电压波形UTα，其它步骤对应替换，得到低压侧上行接触线Tα到高压侧C相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧A相电压波形UA和低压侧上行正馈线Fα电压波形UFα，其它步骤对应替换，得到低压侧上行正馈线Fα到高压侧A相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧B相电压波形UB和低压侧上行正馈线Fα电压波形UFα，其它步骤对应替换，得到低压侧上行正馈线Fα到高压侧B相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧C相电压波形UC和低压侧上行正馈线Fα电压波形UFα，其它步骤对应替换，得到低压侧上行正馈线Fα到高压侧C相的h次谐波传输系数。
或者，进一步的技术方案，
将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧A相电压波形UA和低压侧下行接触线Tβ电压波形UTβ，将步骤3.1替换为选择
且
的时间段所对应的
和
的谐波电压含有率时间段，其它步骤对应替换，得到低压侧下行接触线Tβ到高压侧A相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧B相电压波形UB和低压侧下行接触线Tβ电压波形UTβ，将步骤3.1替换为选择
且
的时间段所对应的
和
的谐波电压含有率时间段，其它步骤对应替换，得到低压侧下行接触线Tβ到高压侧B相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧C相电压波形UC和低压侧下行接触线Tβ电压波形UTβ，将步骤3.1替换为选择
且
的时间段所对应的
和
的谐波电压含有率时间段，其它步骤对应替换，得到低压侧下行接触线Tβ到高压侧C相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧A相电压波形UA和低压侧下行正馈线Fβ电压波形UFβ，将步骤3.1替换为选择
且
的时间段所对应的
和
的谐波电压含有率时间段，其它步骤对应替换，得到低压侧下行正馈线Fβ到高压侧A相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧B相电压波形UB和低压侧下行正馈线Fβ电压波形UFβ，将步骤3.1替换为选择
且
的时间段所对应的
和
的谐波电压含有率时间段，其它步骤对应替换，得到低压侧下行正馈线Fβ到高压侧B相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧C相电压波形UC和低压侧下行正馈线Fβ电压波形UFβ，将步骤3.1替换为选择
且
的时间段所对应的
和
的谐波电压含有率时间段，其它步骤对应替换，得到低压侧下行正馈线Fβ到高压侧C相的h次谐波传输系数；
其中，δ1和δ2为有且只有下行区间有机车经过而上行区间无机车经过的判据。
本发明可用于由机车发射谐波从低压侧上行接触线、上行正馈线、下行接触线、下行正馈线传输到高压侧的传输系数，具有方便、快捷且能较精确反映实际系统谐波传输及谐波放大规律的特点。与现有技术相比，本发明(1)不需要建立牵引供电系统谐波模型；(2)基于实测数据而不是仿真模型，更能反映实际系统的谐波传输特性；(3)定量的估计方法，可以给出较精确的定量的分析结果，而不是定性或粗略的估计。
附图说明
图1是牵引供电系统拓扑结构图。
具体实施方式
如图1所示，一种牵引变压器低压侧到高压侧的谐波传输系数估计方法，利用电压同步测量设备记录牵引变压器110kV/220kV高压侧的A相、B相、C相以及牵引变压器27.5kV低压侧上行接触线Tα、上行正馈线Fα、下行接触线Tβ、下行正馈线Fβ的电压波形，来估计由机车发射谐波从牵引变压器27.5kV低压侧上行接触线Tα、上行正馈线Fα、下行接触线Tβ、下行正馈线Fβ到110kV/220kV高压侧A相、B相、C相的传输系数。如下：
a)利用可实现多节点时间同步测量的电能质量测试设备(有效采样率大于r*H*50Hz，H为要分析的最高谐波次数，r为最高谐波次数的倍数，r>2，通常在2.56到4之间取值)，同步测量牵引变压器高压侧110kV/220kV高压侧A相、B相、C相以及牵引变压器低压侧四条27.5kV供电臂(上行接触线Tα、上行正馈线Fα、下行接触线Tβ、下行正馈线Fβ)的电压波形，以及上行接触线Tα、下行接触线Tβ的电流波形。
b)分别对110kV/220kV高压侧A相、B相、C相，以及上行接触线Tα、上行正馈线Fα、下行接触线Tβ、下行正馈线Fβ的电压波形，利用快速傅里叶变换(FFT)方法提取其基波电压、以及2到H次谐波电压的方均根值的时间序列；进而得到110kV/220kV高压侧的A相、B相、C相处的h次谐波电压含有率时间序列，分别记为
得到牵引变压器低压侧27.5kV上行接触线Tα、上行正馈线Fα、下行接触线Tβ、下行正馈线Fβ处的谐波电压含有率时间序列，分别记为
对上行接触线Tα、下行接触线Tβ的电流波形，利用快速傅里叶变换(FFT)方法提取其基波电流方均根值的时间序列，记
分别为上行接触线Tα、下行接触线Tβ处的基波电流方均根值时间序列。
c)估计由机车发射谐波从27.5kV低压侧上行接触线Tα、上行正馈线Fα传输到110kV/220kV高压侧A相、B相、C相的传输系数，根据上行接触线Tα上基波电流方均根值大于δ1(δ1取30A)且下行接触线Tβ上基波电流方均根值小于δ2(δ2取15A)的判据，选择有且只有上行区间有机车经过的时间段的谐波电压含有率数据。估计由机车发射谐波从27.5kV低压侧下行接触线Tβ、下行正馈线Fβ传输到110kV/220kV高压侧A相、B相、C相的传输系数，根据下行接触线Tβ上基波电流方均根值大于δ1(δ1取30A)且上行接触线Tα上基波电流方均根值小于δ2(δ2取15A)的判据，选择有且只有下行区间有机车经过的时间段的谐波电压含有率数据。
d)以估计h次谐波从27.5kV低压侧上行接触线Tα传输到110kV/220kV高压侧A相为例，将步骤c)得到的谐波电压含有率时间段，采用滑动时间窗分割为宽度为L的子时间段，窗宽L可取100-200，滑动步长取Len＝10，记共获取N段数据。使用二元线性回归方程拟合得到谐波电压含有率时间段中h次谐波电压含有率之间的关系yi＝λixi+bi+εi。其中，i为谐波电压含有率时间段的序号，i＝1,…,N，yi为110kV/220kV高压侧A相处谐波数据段，xi为27.5kV低压侧上行接触线Tα谐波数据段，λi为回归直线的斜率即h次谐波从Tα到高压侧A相的传输系数，bi为截距，εi为回归残差。
e)对步骤d)中的N段数据，分别计算27.5kV低压侧上行接触线Tα、110kV/220kV高压侧A相的每一段h次谐波电压含有率数据的方差
绘制
关于
的散点图，选择
相对于
足够小且具有聚集效应的散点对应的数据段用于估计谐波传输系数，记所选数据段数目为K，重新记该K段数据对应的谐波传输系数λi为
f)通过检验K个传输系数
的聚集效应来判断是否具有精确度，判定其平均值为准确的h次谐波从27.5kV低压侧上行接触线Tα传输到110kV/220kV高压侧A相的谐波传输系数。
g)对于其它的27.5kV低压侧上行接触线Tα、上行正馈线Fα、上行接触线Tβ、上行正馈线Fβ到110kV/220kV高压侧A相、B相、C相的h次谐波传输系数，按照步骤d)到f)的类同方法进行估计。
上述技术方案步骤d)中，选择
相对于
足够小且具有聚集效应的散点对应的数据段，其方法是：记步骤d)所得N个数据段中第i个数据段对应的
为Vi，第i个数据段对应的
为Wi，i＝1,…,N。计算m％*N(m取15)向下取整后的整数个(即K个)最大的Wi的方差与全部N个Vi的方差的比值
若
小于阈值Δ(Δ可取10％)，则认定K个最大的Wi对应的
相对于
足够小且具有聚集效应。否则，取m％＝m％-1％，重新计算。多次重新计算后若仍无满足条件数据，则判定为估计失败。
上述技术方案步骤f)中，通过检验K个传输系数
的聚集效应来判断是否具有精确度，其方法是：计算K个
对应的
的方差与全部N个λi的方差的比值ρ％，若ρ％小于阈值Δ，则判定K个传输系数
具备聚集效应，判定其平均值为准确的h次谐波从27.5kV低压侧上行接触线Tα传输到110kV/220kV高压侧A相的谐波传输系数。
Claims (3)
1.一种牵引变压器低压侧到高压侧的谐波传输系数估计方法，其特征在于，包括：
步骤1：同步获取牵引变压器的电压波形和电流波形：
1.1以大于r*H*50Hz的采样率获取高压侧A相电压波形UA和低压侧上行接触线Tα电压波形UTα；其中，H为要分析的最高谐波次数，r为最高谐波次数的倍数；
1.2获取低压侧上行接触线Tα电流波形ITα和下行接触线Tβ电流波形ITβ；
步骤2：分别对UA和UTα进行快速傅里叶变换，得到其基波电压方均根值时间序列和谐波电压方均根值时间序列，进而得到h次谐波电压含有率时间序列，即
和
h＝2,...,H；分别对ITα和ITβ进行快速傅里叶变换得到其基波电流方均根值时间序列，即
和
步骤3：对任一h次谐波，选择其谐波电压含有率时间段进行分割：
3.1选择
且
的时间段所对应的
和
的谐波电压含有率时间段；其中，δ1和δ2为有且只有上行区间有机车经过而下行区间无机车经过的判据；
3.2将谐波电压含有率时间段用滑动时间窗分割为宽度为L的N个子时间段；
步骤4：对N个子时间段，分别使用二元线性回归方程拟合得到每个子时间段中h次谐波电压含有率之间的关系yi＝λixi+bi+εi；其中，i为子时间段的序号，i＝1,…,N；yi为A相的h次谐波电压含有率，xi为上行接触线Tα的h次谐波电压含有率，λi为回归直线的斜率即h次谐波电压含有率从Tα到A相的传输系数，bi为截距，εi为回归残差；
步骤5：从拟合得到的N个关系中，选择具有聚集效应的K个关系，具体为：
5.1记N个子时间段中，第i个子时间段的
为Vi，第i个子时间段的
为Wi；其中，
为A相第i个子时间段中h次谐波电压含有率的方差，
为上行接触线Tα第i个子时间段中h次谐波电压含有率的方差；
5.2将N个Wi从大到小排序；令m％为需选择的关系的数量的初始百分比；
5.3计算m％*N，向下取整得到K；
5.4计算Var1与Var2的比值
其中，Var1为前K个Wi对应的K个Vi的方差，Var2为N个Vi的方差；
5.5若
小于阈值Δ，则判定前K个Wi对应的K个关系具有聚集效应；否则，令m％＝m％-1％返回步骤5.3重新计算，直到满足
小于阈值Δ；
步骤6：计算Var3与Var4的比值
其中，Var3为所选择的K个关系的λi的方差，Var4为N个关系的λi的方差；若ρ％小于阈值Δ，则取K个关系的λi的平均值，得到低压侧上行接触线Tα到高压侧A相的h次谐波传输系数。
2.如权利要求1所述的谐波传输系数估计方法，其特征在于，
将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧B相电压波形UB和低压侧上行接触线Tα电压波形UTα，其它步骤对应替换，得到低压侧上行接触线Tα到高压侧B相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧C相电压波形UC和低压侧上行接触线Tα电压波形UTα，其它步骤对应替换，得到低压侧上行接触线Tα到高压侧C相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧A相电压波形UA和低压侧上行正馈线Fα电压波形UFα，其它步骤对应替换，得到低压侧上行正馈线Fα到高压侧A相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧B相电压波形UB和低压侧上行正馈线Fα电压波形UFα，其它步骤对应替换，得到低压侧上行正馈线Fα到高压侧B相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧C相电压波形UC和低压侧上行正馈线Fα电压波形UFα，其它步骤对应替换，得到低压侧上行正馈线Fα到高压侧C相的h次谐波传输系数。
3.如权利要求1所述的谐波传输系数估计方法，其特征在于，
将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧A相电压波形UA和低压侧下行接触线Tβ电压波形UTβ，将步骤3.1替换为选择
且
的时间段所对应的
和
的谐波电压含有率时间段，其它步骤对应替换，得到低压侧下行接触线Tβ到高压侧A相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧B相电压波形UB和低压侧下行接触线Tβ电压波形UTβ，将步骤3.1替换为选择
且
的时间段所对应的
和
的谐波电压含有率时间段，其它步骤对应替换，得到低压侧下行接触线Tβ到高压侧B相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧C相电压波形UC和低压侧下行接触线Tβ电压波形UTβ，将步骤3.1替换为选择
且
的时间段所对应的
和
的谐波电压含有率时间段，其它步骤对应替换，得到低压侧下行接触线Tβ到高压侧C相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧A相电压波形UA和低压侧下行正馈线Fβ电压波形UFβ，将步骤3.1替换为选择
且
的时间段所对应的
和
的谐波电压含有率时间段，其它步骤对应替换，得到低压侧下行正馈线Fβ到高压侧A相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧B相电压波形UB和低压侧下行正馈线Fβ电压波形UFβ，将步骤3.1替换为选择
且
的时间段所对应的
和
的谐波电压含有率时间段，其它步骤对应替换，得到低压侧下行正馈线Fβ到高压侧B相的h次谐波传输系数；
或者，将步骤1.1替换为以大于r*H*50Hz的采样率获取高压侧C相电压波形UC和低压侧下行正馈线Fβ电压波形UFβ，将步骤3.1替换为选择
且
的时间段所对应的
和
的谐波电压含有率时间段，其它步骤对应替换，得到低压侧下行正馈线Fβ到高压侧C相的h次谐波传输系数；
其中，δ1和δ2为有且只有下行区间有机车经过而上行区间无机车经过的判据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110462978.XA CN113190798A (zh) | 2021-04-23 | 2021-04-23 | 一种牵引变压器低压侧到高压侧的谐波传输系数估计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110462978.XA CN113190798A (zh) | 2021-04-23 | 2021-04-23 | 一种牵引变压器低压侧到高压侧的谐波传输系数估计方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113190798A true CN113190798A (zh) | 2021-07-30 |
Family
ID=76979907
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110462978.XA Pending CN113190798A (zh) | 2021-04-23 | 2021-04-23 | 一种牵引变压器低压侧到高压侧的谐波传输系数估计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113190798A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114113783A (zh) * | 2021-12-03 | 2022-03-01 | 通号（长沙）轨道交通控制技术有限公司 | 一种电力机车谐波放大倍数估计方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20160064695A (ko) * | 2014-11-28 | 2016-06-08 | 서울과학기술대학교 산학협력단 | 고속 푸리에 변환을 이용한 위상 추종 및 고조파 검출 장치와 그 방법 |
| CN107846670A (zh) * | 2017-11-01 | 2018-03-27 | 东华大学 | 移动群体感知中保护数据隐私的盲回归建模及更新方法 |
| CN110190628A (zh) * | 2019-06-12 | 2019-08-30 | 西南交通大学 | 一种基于综合能源的复线全并联at牵引供电系统及方法 |
| CN111984503A (zh) * | 2020-08-17 | 2020-11-24 | 网宿科技股份有限公司 | 一种监控指标数据异常数据识别的方法及装置 |
| CN112666391A (zh) * | 2020-12-01 | 2021-04-16 | 广东电网有限责任公司电力科学研究院 | 一种牵引变压器两侧谐波电压传输系数计算方法及装置 |
-
2021
- 2021-04-23 CN CN202110462978.XA patent/CN113190798A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20160064695A (ko) * | 2014-11-28 | 2016-06-08 | 서울과학기술대학교 산학협력단 | 고속 푸리에 변환을 이용한 위상 추종 및 고조파 검출 장치와 그 방법 |
| CN107846670A (zh) * | 2017-11-01 | 2018-03-27 | 东华大学 | 移动群体感知中保护数据隐私的盲回归建模及更新方法 |
| CN110190628A (zh) * | 2019-06-12 | 2019-08-30 | 西南交通大学 | 一种基于综合能源的复线全并联at牵引供电系统及方法 |
| CN111984503A (zh) * | 2020-08-17 | 2020-11-24 | 网宿科技股份有限公司 | 一种监控指标数据异常数据识别的方法及装置 |
| CN112666391A (zh) * | 2020-12-01 | 2021-04-16 | 广东电网有限责任公司电力科学研究院 | 一种牵引变压器两侧谐波电压传输系数计算方法及装置 |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114113783A (zh) * | 2021-12-03 | 2022-03-01 | 通号（长沙）轨道交通控制技术有限公司 | 一种电力机车谐波放大倍数估计方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112858784B (zh) | 一种牵引供电系统-区域电网并联谐波谐振频率识别方法 | |
| CN107064752A (zh) | 一种航空故障电弧检测的判别算法 | |
| CN105372531A (zh) | 基于Weibull分布模型的变压器绝缘热老化参数相关性计算方法 | |
| CN109375051B (zh) | 基于频谱密度衰减的雷电暂态信号识别方法及系统 | |
| CN106371002A (zh) | 一种基于希尔伯特黄变换算法对断路器故障诊断的方法 | |
| CN104410360A (zh) | 光伏发电系统安全运行方法及其中对人工神经网络的训练方法和实时检测方法及装置 | |
| CN113466597A (zh) | 一种轨道交通供电系统设备状态智能检测方法 | |
| CN113190798A (zh) | 一种牵引变压器低压侧到高压侧的谐波传输系数估计方法 | |
| CN112666391B (zh) | 一种牵引变压器两侧谐波电压传输系数计算方法及装置 | |
| CN103941149A (zh) | 电气化铁道牵引网故障定位方法 | |
| CN108008246B (zh) | 基于时域自适应窗的线路故障测距方法及装置 | |
| Mariscotti | Characterization of Power Quality transient phenomena of DC railway traction supply | |
| CN109711687B (zh) | 一种基于改进熵值法的绝缘子状态模糊评估方法 | |
| CN116990621B (zh) | 电力安全工器具故障动态监测系统 | |
| Seferi et al. | Power quality event analysis in 25 kV 50 Hz AC railway system networks | |
| CN102053213A (zh) | 基于信息熵的特高压直流输电线路区内外故障识别方法 | |
| Tian et al. | A new scheme to identify the special load current waveform for a Metro DC traction system | |
| CN112782528B (zh) | 一种利用pmu的配电网故障区段定位方法 | |
| CN104901287A (zh) | 基于过渡电阻与故障初始角的输电线路暂态保护方法 | |
| CN115047261A (zh) | 一种基于熵权法的高铁用避雷器寿命评估方法及其应用 | |
| CN106199345A (zh) | 高压架空输电线路绝缘子污秽监测新方法 | |
| CN202794394U (zh) | 接触网避雷器在线检测系统 | |
| CN112406966A (zh) | 区间内单辆机车运行位置定位方法 | |
| Havryliuk | The accuracy of traction current harmonics parameters determination by windowed FFT | |
| CN104360297A (zh) | 一种基于瞬时值比较的pmu动态性能多指标测试方法 |
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 |