CN113281614B - 一种电力电缆宽频阻抗谱测试方法 - Google Patents
一种电力电缆宽频阻抗谱测试方法 Download PDFInfo
- Publication number
- CN113281614B CN113281614B CN202110541444.6A CN202110541444A CN113281614B CN 113281614 B CN113281614 B CN 113281614B CN 202110541444 A CN202110541444 A CN 202110541444A CN 113281614 B CN113281614 B CN 113281614B
- Authority
- CN
- China
- Prior art keywords
- power cable
- impedance spectrum
- tested
- broadband impedance
- calculating
- 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
- 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
- G01R31/08—Locating faults in cables, transmission lines, or networks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
- G01R27/26—Measuring inductance or capacitance; Measuring quality factor, e.g. by using the resonance method; Measuring loss factor; Measuring dielectric constants ; Measuring impedance or related variables
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Resistance Or Impedance (AREA)
- Testing Relating To Insulation (AREA)
Abstract
本发明涉及电力***技术领域,具体公开了一种电力电缆宽频阻抗谱测试方法,其中,用于宽频阻抗谱测试***中,所述宽频阻抗谱测试***包括矢量网络分析仪和与所述矢量网络分析仪连接的屏蔽接线盒,所述屏蔽接线盒连接被测电力电缆的终端,所述电力电缆宽频阻抗谱测试方法包括:步骤S110、在对所述宽频阻抗谱测试***进行校准后,获取被测电力电缆的宽频阻抗谱数据;步骤S120、根据所述被测电力电缆的宽频阻抗谱数据计算被测电力电缆的特性阻抗;步骤S130、根据所述被测电力电缆的特性阻抗计算修正后的被测电力电缆的宽频阻抗谱。本发明提供的电力电缆宽频阻抗谱测试方法能够消除现场宽频阻抗谱测试过程中分叉引线杂散参数及感应电影响。
Description
技术领域
本发明涉及电力***技术领域,尤其涉及一种电力电缆宽频阻抗谱测试方法。
背景技术
随着城市化进程的快速发展,具有优良电气绝缘性能和机械特性的交联聚乙烯(XLPE)电力电缆开始广泛地应用于电力***输配电网络中。然而,随着运行年限的增加,电力电缆由于受到电、热、机械等多重应力及运行环境的影响,开始出现不同程度的绝缘老化问题,部分电缆达到甚至超过当初的设计使用寿命,其绝缘问题开始逐渐显现。因此,为了提高电力电缆的供电可靠性,保证电网安全,需要定期对电力电缆的绝缘状态进行检测,从而评估电力电缆的绝缘状态并实现电力电缆局部缺陷的定位。
现有的宽频阻抗谱法是实现电力电缆绝缘状态评估和缺陷定位的有效方法。申请号为201610078115.1,名称为“一种电缆运行状态诊断方法及***”的中国发明专利,提供了一种使用电缆宽频阻抗谱评估电缆绝缘状态并实现电缆局部缺陷定位的方法和***,但该方法需要获得完好电缆的宽频阻抗谱作为参考,并且该方法在现场应用时测试得到的宽频阻抗谱受仪器测试分叉引线杂散参数影响较大。申请号为201711154670.9,名称为“一种基于阻抗频谱分析的电缆局部缺陷评估的装置和方法”的中国发明专利,提供了一种利用电缆宽频阻抗评估电缆局部缺陷的装置和方法,该方法虽无需获得完好电缆的宽频阻抗谱作为参考,但其在现场应用时测试结果同样会受到测试仪器测试分叉引线杂散参数的影响。可见,在现场应用中开展消除分叉引线杂散参数影响的电力电缆宽频阻抗谱准确测试技术具有重大意义。除此之外,现有基于宽频阻抗法的电力电缆绝缘状态评估和缺陷定位方法并未考虑现场感应电的影响,从而在测试过程中存在着一定的安全隐患。
因此,如何能够提供一种消除现场宽频阻抗谱测试过程中分叉引线杂散参数及感应电影响的测试方法成为本领域技术人员亟待解决的技术问题。
发明内容
本发明提供了一种电力电缆宽频阻抗谱测试方法,解决相关技术中存在的现场宽频阻抗谱测试过程中感应电影响无法消除的问题。
作为本发明的一个方面,提供一种电力电缆宽频阻抗谱测试方法,其中,用于宽频阻抗谱测试***中,所述宽频阻抗谱测试***包括矢量网络分析仪和与所述矢量网络分析仪连接的屏蔽接线盒,所述屏蔽接线盒连接被测电力电缆的终端,所述电力电缆宽频阻抗谱测试方法包括:
步骤S110、在对所述宽频阻抗谱测试***进行校准后,获取被测电力电缆的宽频阻抗谱数据;
步骤S120、根据所述被测电力电缆的宽频阻抗谱数据计算被测电力电缆的特性阻抗;
步骤S130、根据所述被测电力电缆的特性阻抗计算修正后的被测电力电缆的宽频阻抗谱。
进一步地,所述矢量网络分析仪与所述屏蔽接线盒之间以及所述屏蔽接线盒与所述被测电力电缆之间均通过同轴电缆连接,连接所述屏蔽接线盒的同轴电缆与所述被测电力电缆之间还设置分叉引线,所述分叉引线的一端连接所述同轴电缆,另一端连接所述被测电力电缆;
所述对所述宽频阻抗谱测试***进行校准包括:
根据所述矢量网络分析仪的开路、短路和负载校准功能对所述屏蔽接线盒及其两侧的同轴电缆在设置的扫频范围内进行校准。
进一步地,所述获取被测电力电缆的宽频阻抗谱数据,包括:
根据宽频阻抗谱计算公式计算被测电力电缆的宽频阻抗谱复数值,其中所述宽频阻抗谱计算公式为:
Zd(f)=|Zd(f)|·exp(j·Angle(Zd(f))),
其中,Zd(f)表示被测电力电缆的宽频阻抗谱的幅值,Angle(Zd(f))表示被测电力电缆的宽频阻抗谱的相位,f表示被测电力电缆的信号频率,j表示复数,exp(·)表示指数函数;
根据所述被测电力电缆的宽频阻抗谱复数值计算得到首端反射系数复数值,其中所述首端反射系数复数值表示为:
其中,Γd(f)表示首端反射系数复数值。
进一步地,所述根据所述被测电力电缆的宽频阻抗谱数据计算被测电力电缆的特性阻抗,包括:
步骤S121、根据公式y0(t)=exp(-4·ln(2)·t2/w2)产生脉宽为w的高斯脉冲,并根据快速傅里叶变换对y0(t)进行处理得到Y0(f)=FFT(y0(t)),根据公式Y(f)=Y0(f)·Γd(f)计算得到反射信号y(t)的频域表达式Y(f),其中t表示时间;
步骤S122、对Y(f)进行快速傅里叶逆变换处理得到时域反射信号y(t)=IFFT(Y(f)),提取y(t)中表征分叉引线处反射的时域信号y1(t),并进行快速傅里叶变换处理得到Y1(f)=FFT(y1(t));
步骤S123、根据公式Γ1(f)=Y1(f)/Y0(f)计算分叉引线处的反射系数Γ1(f),并假设被测电力电缆的特性阻抗Z0取值为10Ω;
步骤S124、计算得到交叉引线处的等效阻抗Z1(f),
步骤S125、根据公式Z′d(f)=Zd(f)-Z1(f)计算得到对应的Z′d(f),并根据转换函数f→t′将得到的Z′d(f)相位Angle(Z′d(f))转换为t′域信号Angle(Z′d(t′)),对Angle(Z′d(t′))进行快速傅里叶变换处理得到Angle(Z′d(t′))的频谱FFT(Angle(Z′d(t′))),保留FFT(Angle(Z′d(t′)))中除直流分量外能量最大点的信号,其他部分全部置零并进行快速傅里叶逆变换处理得到信号Angle(Z′d·max(t′)),利用转换函数t′→f将t′域信号转换为频域信号Angle(Z′d·max(f));
步骤S126、计算Z′d(f)的相位Angle(Z′d(f))与步骤S125中计算得到的Angle(Z′d·max(f))之间的皮尔逊相关系数ρZ:
其中,N表示进行计算的数据长度,mean(·)表示平均值;
步骤S127、以步长0.1Ω增大假设的被测电力电缆的特性阻抗Z0:
Z0=Z0+0.1;
步骤S128、重复步骤S124~步骤S127直至假设的被测电力电缆的特性阻抗Z0取值为100Ω;
步骤S129、获取不同假设特性阻抗下步骤S126中计算得到的皮尔逊相关系数曲线,并将曲线中ρZ最大时对应的特性阻抗作为估计得到的被测电力电缆特性阻抗Z0·E。
进一步地,所述根据所述被测电力电缆的特性阻抗计算修正后的被测电力电缆的宽频阻抗谱,包括:
计算被测电力电缆特性阻抗Z0·E下交叉引线处的等效阻抗Z1·E(f):
根据公式Z′d·E(f)=Zd(f)-Z1·E(f)计算得到Z′d·E(f)为修正后的被测电力电缆的宽频阻抗谱。
进一步地,所述高斯脉冲的脉宽w不大于所述矢量网络分析仪设置的扫频范围频率上限值的倒数,且不小于所述矢量网络分析仪设置的扫频范围频率上限值倒数的0.5倍。
进一步地,所述进行计算的数据长度N不大于所述矢量网络分析仪设置的扫频点数的0.5倍。
进一步地,所述同轴电缆的特性阻抗为50Ω。
进一步地,所述屏蔽接线盒的进出线均设置射频连接器,且所述屏蔽接线盒内串联电容。
进一步地,所述屏蔽接线盒内串联电容的取值范围为1μF≤C≤10μF。
本发明实施例提供的电力电缆宽频阻抗谱测试方法与现有技术相比具有以下优势:
1、本发明采用的装有串联电容的屏蔽接线盒可以显著抑制感应的影响,保证了测试过程中测试仪器的安全性,防止仪器损坏。
2、本发明采用的装有串联电容的屏蔽接线盒由于进出线均采用了射频连接器,且盒内串联电容焊接在印制电路板上,整个测试盒内无软导线连接,因此保证了屏蔽接线盒的等效电路参数在测试过程中的稳定性。
3、本发明通过时域分析和频域分析相结合的方式,可有效提取表征测试仪器分叉引线反射部分的信号,从而进一步得到交叉引线处的等效阻抗。
4、本发明通过皮尔逊相关系数法直接实现了被测电力电缆特性阻抗的估计,相比于传统的电力电缆特性阻抗计算方法,本发明无需对电缆终端进行处理,同时也不受测试仪器分叉引线的影响。
5、本发明通过估计得到的特性阻抗结合交叉引线处的等效阻抗实现了被测电力电缆宽频阻抗谱的修正,消除了测试仪器分叉引线的影响,从而大大提高了测试数据的准确性。
附图说明
附图是用来提供对本发明的进一步理解,并且构成说明书的一部分,与下面的具体实施方式一起用于解释本发明,但并不构成对本发明的限制。在附图中:
图1为本发明提供的宽频阻抗谱测试***的结构示意图。
图2为本发明提供的电力电缆宽频阻抗谱测试方法的流程图。
图3为本发明提供的宽频阻抗谱测试***校准前后的阻抗谱。
图4a为本发明提供的电力电缆宽频阻抗谱测试方法的具体实施方式流程图。
图4b为本发明提供的被测电力电缆特性阻抗计算方法流程图。
图5为本发明提供的一种实施例中测试得到的被测电缆的宽频阻抗谱幅频特性曲线和参考曲线,其中实线表示被测电缆的宽频阻抗谱幅频特性曲线,虚线表示参考曲线。
图6为本发明提供的一种实施例中测试得到的被测电缆的宽频阻抗谱相频特性曲线和参考曲线,其中实线表示被测电缆的宽频阻抗谱相频特性曲线,虚线表示参考曲线。
图7为本发明提供的一种实施例中不同假设特性阻抗下计算得到的皮尔逊相关系数曲线。
图8为本发明提供的一种实施例中修正得到的被测电缆的宽频阻抗谱幅频特性曲线和相频特性曲线,其中实线表示幅频特性曲线,虚线表示相频特性曲线。
图9为本发明提供的另一种实施例中测试得到的被测电缆的宽频阻抗幅频特性曲线和相频特性曲线,其中实线表示幅频特性曲线,虚线表示相频特性曲线。
图10为本发明提供的另一种实施例中修正得到的被测电缆的宽频阻抗谱幅频特性曲线和相频特性曲线,其中实线表示幅频特性曲线,虚线表示相频特性曲线。
具体实施方式
需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互结合。下面将参考附图并结合实施例来详细说明本发明。
为了使本领域技术人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包括,例如,包含了一系列步骤或单元的过程、方法、***、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
在本实施例中提供了一种电力电缆宽频阻抗谱测试方法,用于宽频阻抗谱测试***中,如图1所示,所述宽频阻抗谱测试***包括矢量网络分析仪100和与所述矢量网络分析仪100连接的屏蔽接线盒200,所述屏蔽接线盒200连接被测电力电缆500的终端510,如图2所示,为所述电力电缆宽频阻抗谱测试方法的流程图,如图2所示,所述电力电缆宽频阻抗谱测试方法包括:
步骤S110、在对所述宽频阻抗谱测试***进行校准后,获取被测电力电缆的宽频阻抗谱数据;
在本发明实施例中,所述矢量网络分析仪100与所述屏蔽接线盒200之间以及所述屏蔽接线盒200与所述被测电力电缆500之间均通过同轴电缆300连接,连接所述屏蔽接线盒200的同轴电缆300与所述被测电力电缆500之间还设置分叉引线310,所述分叉引线310的一端连接所述同轴电缆300,另一端连接所述被测电力电缆500,在本发明实施例中,所述分叉引线310的另一端连接所述被测电力电缆500的终端510;
所述对所述宽频阻抗谱测试***进行校准包括:
根据所述矢量网络分析仪的开路、短路和负载校准功能对所述屏蔽接线盒及其两侧的同轴电缆在设置的扫频范围内进行校准。
在本发明实施例中,所述同轴电缆的特性阻抗为50Ω。
具体地,如图1所示,所述屏蔽接线盒200的进出线均设置射频连接器400,且所述屏蔽接线盒200内串联电容。
优选地,所述屏蔽接线盒200内串联电容的取值范围为1μF≤C≤10μF。
在进行校准时,将装有串联电容的屏蔽接线盒200利用特性阻抗为50Ω的同轴电缆连接至矢量网络分析仪100,然后利用矢量网络分析仪100的开路、短路和负载校准功能对分叉引线310及电容屏蔽接线盒在设置的扫频范围内进行校准,从而将测试平面A校准至同轴电缆300的末端。
将分叉引线310的同轴端连接至测试平面A,并将另一端连接至被测电力电缆500,然后对被测电力电缆500的宽频阻抗谱的幅值Zd(f)和相位Angle(Zd(f))进行测试,其中f为信号频率。
在本发明实施例中,屏蔽接线盒200内串联电容焊接在印制电路板上,整个测试盒内不通过软导线连接,其主要目的是保证屏蔽接线盒的等效电路参数在测试过程中保持稳定。串联电容的目的主要是将工频感应电压限制在串联电容上,从而使降落在网络分析仪内阻上的工频感应电压尽可能小。
还需要说明的是,当现场测试过程中不存在感应电影响时,可不安装装有串联电容的屏蔽接线盒。
矢量网络分析仪(如安捷伦E5061B)的校准方法可参见现有技术中已经披露的处理方法,如Joel P D.Handbook of Microwave Component Measurements with AdvancedVNA Techniques[M].United Kingdom:John Wiley&Sons,2012,其主要目的是消除同轴电缆和装有电容的屏蔽接线盒对测试结果的影响。
具体地,所述获取被测电力电缆的宽频阻抗谱数据,包括:
根据宽频阻抗谱计算公式计算被测电力电缆的宽频阻抗谱复数值,其中所述宽频阻抗谱计算公式为:
Zd(f)=|Zd(f)|·exp(j·Angle(Zd(f))),
其中,Zd(f)表示被测电力电缆的宽频阻抗谱的幅值,Angle(Zd(f))表示被测电力电缆的宽频阻抗谱的相位,f表示被测电力电缆的信号频率,j表示复数,exp(·)表示指数函数;
根据所述被测电力电缆的宽频阻抗谱复数值计算得到首端反射系数复数值,其中所述首端反射系数复数值表示为:
其中,Γd(f)表示首端反射系数复数值。
步骤S120、根据所述被测电力电缆的宽频阻抗谱数据计算被测电力电缆的特性阻抗;
具体包括如下步骤:
步骤S121、根据公式y0(t)=exp(-4·ln(2)·t2/w2)产生脉宽为w的高斯脉冲,并根据快速傅里叶变换对y0(t)进行处理得到Y0(f)=FFT(y0(t)),根据公式Y(f)=Y0(f)·Γd(f)计算得到反射信号y(t)的频域表达式Y(f),其中t表示时间;
步骤S122、对Y(f)进行快速傅里叶逆变换处理得到时域反射信号y(t)=IFFT(Y(f)),提取y(t)中表征分叉引线处反射的时域信号y1(t),并进行快速傅里叶变换处理得到Y1(f)=FFT(y1(t));
步骤S123、根据公式Γ1(f)=Y1(f)/Y0(f)计算分叉引线处的反射系数Γ1(f),并假设被测电力电缆的特性阻抗Z0取值为10Ω;
步骤S124、计算得到交叉引线处的等效阻抗Z1(f),
步骤S125、根据公式Z′d(f)=Zd(f)-Z1(f)计算得到对应的Z′d(f),并根据转换函数f→t′将得到的Z′d(f)相位Angle(Z′d(f))转换为t′域信号Angle(Z′d(t′)),对Angle(Z′d(t′))进行快速傅里叶变换处理得到Angle(Z′d(t′))的频谱FFT(Angle(Z′d(t′))),保留FFT(Angle(Z′d(t′)))中除直流分量外能量最大点的信号,其他部分全部置零并进行快速傅里叶逆变换处理得到信号Angle(Z′d·max(t′)),利用转换函数t′→f将t′域信号转换为频域信号Angle(Z′d·max(f));
步骤S126、计算Z′d(f)的相位Angle(Z′d(f))与步骤S125中计算得到的Angle(Z′d·max(f))之间的皮尔逊相关系数ρZ:
其中,N表示进行计算的数据长度,mean(·)表示平均值;
步骤S127、以步长0.1Ω增大假设的被测电力电缆的特性阻抗Z0:
Z0=Z0+0.1;
步骤S128、重复步骤S124~步骤S127直至假设的被测电力电缆的特性阻抗Z0取值为100Ω;
步骤S129、获取不同假设特性阻抗下步骤S126中计算得到的皮尔逊相关系数曲线,并将曲线中ρZ最大时对应的特性阻抗作为估计得到的被测电力电缆特性阻抗Z0·E。
步骤S130、根据所述被测电力电缆的特性阻抗计算修正后的被测电力电缆的宽频阻抗谱。
具体地,计算被测电力电缆特性阻抗Z0·E下交叉引线处的等效阻抗Z1·E(f):
根据公式Z′d·E(f)=Zd(f)-Z1·E(f)计算得到Z′d·E(f)为修正后的被测电力电缆的宽频阻抗谱。
综上,本发明实施例提供的电力电缆宽频阻抗谱测试方法与现有技术相比具有以下优势:
1、本发明采用的装有串联电容的屏蔽接线盒可以显著抑制感应的影响,保证了测试过程中测试仪器的安全性,防止仪器损坏。
2、本发明采用的装有串联电容的屏蔽接线盒由于进出线均采用了射频连接器,且盒内串联电容焊接在印制电路板上,整个测试盒内无软导线连接,因此保证了屏蔽接线盒的等效电路参数在测试过程中的稳定性。
3、本发明通过时域分析和频域分析相结合的方式,可有效提取表征测试仪器分叉引线反射部分的信号,从而进一步得到交叉引线处的等效阻抗。
4、本发明通过皮尔逊相关系数法直接实现了被测电力电缆特性阻抗的估计,相比于传统的电力电缆特性阻抗计算方法,本发明无需对电缆终端进行处理,同时也不受测试仪器分叉引线的影响。
5、本发明通过估计得到的特性阻抗结合交叉引线处的等效阻抗实现了被测电力电缆宽频阻抗谱的修正,消除了测试仪器分叉引线的影响,从而大大提高了测试数据的准确性。
下面对本发明提供的电力电缆宽频阻抗谱测试方法进行具体说明。
作为一种具体地实施例,本实施例针对的被测电力电缆是508m ZR-YJLV02 8.7/15 3*25mm2电力电缆,且在100m位置和250m位置各存在一个中间接头。
本实施例采用安捷伦E5061B矢量网络分析仪对被测电力电缆进行宽频阻抗谱测试。
本实施例提供的电力电缆宽频阻抗谱准确测试方法,如图4a和图4b所示,包括以下步骤:
S1、测试***接线及校准
将装有1μF串联电容的屏蔽接线盒利用特性阻抗为50Ω的同轴电缆连接至安捷伦矢量网络分析仪E5061B。设置扫频范围为100kHz~7MHz,扫频点数为700,然后利用E5061B的开路、短路和负载校准功能对测试引线及电容屏蔽接线盒在设置的扫频范围内进行校准,将测试平面校准至同轴电缆的末端,校准前后测试***的阻抗谱相位曲线如图3所示。
S2、电力电缆宽频阻抗谱测试
将分叉引线的同轴端连接至测试平面,并将另一端连接至被测电缆,测试***接线原理图如图4所示,然后对电缆的宽频阻抗谱的幅值和相位进行测试,其中f为信号频率,测试得到的宽频阻抗谱幅频特性曲线、相频特性曲线分别如图5、图6中实线所示。为便于对比,图5、图6中亦给出了与被测电力电缆同型号、同批次且不含中间接头的电缆利用现有技术(参见Papazyan R,Pettersson P,Edin H,et al.Extraction of high frequencypower cable characteristics from S-parameter measurements[J].IEEETransactions on Dielectrics and Electrical Insulation,2004,11(3):461-470)在实验室的测试结果等效值(等效为508m),如图中虚线所示。
由图5、图6结果可知,相比于实验室测量得到的参考曲线,受分叉引线的影响,测量得到的曲线存在着明显的位移现象。
S3、宽频阻抗谱数据预处理
利用公式Zd(f)=|Zd(f)|·exp(j·Angle(Zd(f)))计算得到被测电缆的宽频阻抗谱复数值,其中j为复数,exp(·)为指数函数。然后根据被测电缆宽频阻抗谱复数值计算得到的首端反射系数复数值Γd(f):
S4、计算被测电力电缆的特性阻抗
S41、利用公式y0(t)=exp(-4·ln(2)·t2/w2)产生脉宽为100ns的高斯脉冲,并利用FFT对y0(t)进行处理得到Y0(f)=FFT(y0(t)),然后利用公式Y(f)=Y0(f)·Γd(f)计算得到反射信号y(t)的频域表达式Y(f)。其中t为时间。
S42、对Y(f)进行IFFT处理得到时域反射信号y(t)=IFFT(Y(f)),然后提取y(t)中表征分叉引线处反射的时域信号y1(t),并进行FFT处理得到Y1(f)=FFT(y1(t))。
S43、利用公式Γ1(f)=Y1(f)/Y0(f)计算分叉引线处的反射系数Γ1(f),并假设被测电力电缆特性阻抗Z0取值为10Ω。
S44、计算得到交叉引线处的等效阻抗Z1(f):
S45、利用公式Z′d(f)=Zd(f)-Z1(f)计算得到对应的Z′d(f),然后利用转换函数f→t′将得到的Z′d(f)相位Angle(Z′d(f))转换为t′域信号Angle(Z′d(t′)),并对Angle(Z′d(t′))进行快速FFT处理,得到Angle(Z′d(t′))的频谱FFT(Angle(Z′d(t′)))。保留FFT(Angle(Z′d(t′)))中除直流分量外能量最大点的信号,其他部分全部置零并进行IFFT处理得到信号Angle(Z′d·max(t′)),然后利用转换函数t′→f将t′域信号转换为频域信号Angle(Z′d·max(f))。
S46、计算Z′d(f)的相位Angle(Z′d(f))与S45中计算得到的Angle(Z′d·max(f))的皮尔逊相关系数ρZ:
其中,N为进行计算的数据长度,取值为200;mean(·)为平均值。
S47、以步长0.1Ω增大假设的电力电缆特性阻抗Z0:
Z0=Z0+0.1;
S48、重复步骤S44~S47直至假设的电力电缆特性阻抗Z0取值为100Ω。
S49、计算不同假设特性阻抗下的皮尔逊相关系数曲线如图7所示,获取皮尔逊相关系数曲线中最大值对应的特性阻抗,并将此特性阻抗作为估计得到的电力电缆特性阻抗Z0·E。由皮尔逊相关系数曲线可知被测电力电缆的特性阻抗为37Ω。
S5、计算修正得到的宽频阻抗谱
计算特性阻抗Z0·E下交叉引线处的等效阻抗Z1·E(f):
利用公式Z′d·E(f)=Zd(f)-Z1·E(f)计算得到Z′d·E(f)即为修正得到的宽频阻抗谱,修正后得到的宽频阻抗谱的幅频特性曲线、相频特性曲线如图8所示。由图8结果可知,修正后的宽频阻抗谱可有效解决原始数据存在的位移现象。
作为一种具体地实施例,本实施例针对的被测电力电缆是1km ZR-YJV228.7/153*400mm2电力电缆,且在243m位置和628m位置分别存在一个中间接头。
本实施例采用安捷伦E5061B矢量网络分析仪对被测电力电缆进行宽频阻抗谱测试。
本实施例提供的电力电缆宽频阻抗谱准确测试方法,如图4a和图4b所示,包括以下步骤:
S1、测试***接线及校准
将装有1μF串联电容的屏蔽接线盒利用特性阻抗为50Ω的同轴电缆连接至安捷伦矢量网络分析仪E5061B。设置扫频范围为50kHz~5MHz,扫频点数为1000,然后利用E5061B的开路、短路和负载校准功能对测试引线及电容屏蔽接线盒在设置的扫频范围内进行校准,将测试平面校准至同轴电缆的末端。
S2、电力电缆宽频阻抗谱测试
将分叉引线的同轴端连接至测试平面,并将另一端连接至被测电缆,然后对电缆的宽频阻抗谱的幅值Zd(f)|和相位Angle(Zd(f))进行测试,其中f为信号频率,测试得到的宽频阻抗谱幅频特性曲线、相频特性曲线分别如图9中实线、虚线所示。由图9结果可知,受分叉引线的影响,测量得到的曲线存在着明显的位移现象。
S3、宽频阻抗谱数据预处理
利用公式Zd(f)=|Zd(f)|·exp(j·Angle(Zd(f)))计算得到被测电缆的宽频阻抗谱复数值,其中j为复数,exp(·)为为指数函数。然后根据被测电缆宽频阻抗谱复数值计算得到的首端反射系数复数值Γ1(f):
S4、计算被测电力电缆的特性阻抗
S41、利用公式y0(t)=exp(-4·ln(2)·t2/w2)产生脉宽为100ns的高斯脉冲,并利用FFT对y0(t)进行处理得到Y0(f)=FFT(y0(t)),然后利用公式Y(f)=Y0(f)·Γd(f)计算得到反射信号y(t)的频域表达式Y(f)。其中t为时间。
S42、对Y(f)进行IFFT处理得到时域反射信号y(t)=IFFT(Y(f)),然后提取y(t)中表征分叉引线处反射的时域信号y1(t),并进行FFT处理得到Y1(f)=FFT(y1(t))。
S43、利用公式Γ1(f)=Y1(f)/Y0(f)计算分叉引线处的反射系数Γ1(f),并假设被测电力电缆特性阻抗Z0取值为10Ω。
S44、计算得到交叉引线处的等效阻抗Z1(f):
S45、利用公式Z′d(f)=Zd(f)-Z1(f)计算得到对应的Z′d(f),然后利用转换函数f→t′将得到的Z′d(f)相位Angle(Z′d(f))转换为t′域信号Angle(Z′d(t′)),并对Angle(Z′d(t′))进行快速FFT处理,得到Angle(Z′d(t′))的频谱FFT(Angle(Z′d(t′)))。保留FFT(Angle(Z′d(t′)))中除直流分量外能量最大点的信号,其他部分全部置零并进行IFFT处理得到信号Angle(Z′d·max(t′)),然后利用转换函数t′→f将t′域信号转换为频域信号Angle(Z′d·max(f))。
S46、计算Z′d(f)的相位Angle(Z′d(f))与S45中计算得到的Angle(Z′d·max(f))的皮尔逊相关系数ρZ:
其中,N为进行计算的数据长度,取值为200;mean(·)为平均值。
S47、以步长0.1Ω增大假设的电力电缆特性阻抗Z0:
Z0=Z0+0.1;
S48、重复步骤S44~S47直至假设的电力电缆特性阻抗Z0取值为100Ω。
S49、计算得到不同假设特性阻抗下的皮尔逊相关系数曲线,然后获取皮尔逊相关系数曲线中最大值对应的特性阻抗,并将此特性阻抗作为估计得到的电力电缆特性阻抗Z0·E。由皮尔逊相关系数曲线可知被测电力电缆的特性阻抗为48.3Ω。
S5、计算修正得到的宽频阻抗谱
计算特性阻抗Z0·E下交叉引线处的等效阻抗Z1·E(f):
利用公式Z′d·E(f)=Zd(f)-Z1·E(f)计算得到Z′d·E(f)即为修正得到的宽频阻抗谱,修正后得到的宽频阻抗谱的幅频特性曲线、相频特性曲线如图10所示。由图10结果可知,修正后的宽频阻抗谱可有效解决原始数据存在的位移现象。
可以理解的是,以上实施方式仅仅是为了说明本发明的原理而采用的示例性实施方式,然而本发明并不局限于此。对于本领域内的普通技术人员而言,在不脱离本发明的精神和实质的情况下,可以做出各种变型和改进,这些变型和改进也视为本发明的保护范围。
Claims (7)
1.一种电力电缆宽频阻抗谱测试方法,其特征在于,用于宽频阻抗谱测试***中,所述宽频阻抗谱测试***包括矢量网络分析仪和与所述矢量网络分析仪连接的屏蔽接线盒,所述屏蔽接线盒连接被测电力电缆的终端,所述电力电缆宽频阻抗谱测试方法包括:
步骤S110、在对所述宽频阻抗谱测试***进行校准后,获取被测电力电缆的宽频阻抗谱数据;
步骤S120、根据所述被测电力电缆的宽频阻抗谱数据计算被测电力电缆的特性阻抗;
步骤S130、根据所述被测电力电缆的特性阻抗计算修正后的被测电力电缆的宽频阻抗谱;
所述矢量网络分析仪与所述屏蔽接线盒之间以及所述屏蔽接线盒与所述被测电力电缆之间均通过同轴电缆连接,连接所述屏蔽接线盒的同轴电缆与所述被测电力电缆之间还设置分叉引线,所述分叉引线的一端连接所述同轴电缆,另一端连接所述被测电力电缆;
所述对所述宽频阻抗谱测试***进行校准包括:
根据所述矢量网络分析仪的开路、短路和负载校准功能对所述屏蔽接线盒及其两侧的同轴电缆在设置的扫频范围内进行校准;
所述获取被测电力电缆的宽频阻抗谱数据,包括:
根据宽频阻抗谱计算公式计算被测电力电缆的宽频阻抗谱复数值,其中所述宽频阻抗谱计算公式为:
Zd(f)=|Zd(f)|·exp(j·Angle(Zd(f))),
其中,Zd(f)表示被测电力电缆的宽频阻抗谱的幅值,Angle(Zd(f))表示被测电力电缆的宽频阻抗谱的相位,f表示被测电力电缆的信号频率,j表示复数,exp(·)表示指数函数;
根据所述被测电力电缆的宽频阻抗谱复数值计算得到首端反射系数复数值,其中所述首端反射系数复数值表示为:
其中,Γd(f)表示首端反射系数复数值;
所述根据所述被测电力电缆的宽频阻抗谱数据计算被测电力电缆的特性阻抗,包括:
步骤S121、根据公式y0(t)=exp(-4·ln(2)·t2/w2)产生脉宽为w的高斯脉冲,并根据快速傅里叶变换对y0(t)进行处理得到Y0(f)=FFT(y0(t)),根据公式Y(f)=Y0(f)·Γd(f)计算得到反射信号y(t)的频域表达式Y(f),其中t表示时间;
步骤S122、对Y(f)进行快速傅里叶逆变换处理得到时域反射信号y(t)=IFFT(Y(f)),提取y(t)中表征分叉引线处反射的时域信号y1(t),并进行快速傅里叶变换处理得到Y1(f)=FFT(y1(t));
步骤S123、根据公式Γ1(f)=Y1(f)/Y0(f)计算分叉引线处的反射系数Γ1(f),并假设被测电力电缆的特性阻抗Z0取值为10Ω;
步骤S124、计算得到交叉引线处的等效阻抗Z1(f),
步骤S125、根据公式Z′d(f)=Zd(f)-Z1(f)计算得到对应的Z′d(f),并根据转换函数f→t′将得到的Z′d(f)相位Angle(Z′d(f))转换为t′域信号Angle(Z′d(t′)),对Angle(Z′d(t′))进行快速傅里叶变换处理得到Angle(Z′d(t′))的频谱FFT(Angle(Z′d(t′))),保留FFT(Angle(Z′d(t′)))中除直流分量外能量最大点的信号,其他部分全部置零并进行快速傅里叶逆变换处理得到信号Angle(Z′d·max(t′)),利用转换函数t′→f将t′域信号转换为频域信号Angle(Z′d·max(f));
步骤S126、计算Z′d(f)的相位Angle(Z′d(f))与步骤S125中计算得到的Angle(Z′d·max(f))之间的皮尔逊相关系数ρZ:
其中,N表示进行计算的数据长度,mean(·)表示平均值;
步骤S127、以步长0.1Ω增大假设的被测电力电缆的特性阻抗Z0:
Z0=Z0+0.1;
步骤S128、重复步骤S124~步骤S127直至假设的被测电力电缆的特性阻抗Z0取值为100Ω;
步骤S129、获取不同假设特性阻抗下步骤S126中计算得到的皮尔逊相关系数曲线,并将曲线中ρZ最大时对应的特性阻抗作为估计得到的被测电力电缆特性阻抗Z0·E。
2.根据权利要求1所述的电力电缆宽频阻抗谱测试方法,其特征在于,所述根据所述被测电力电缆的特性阻抗计算修正后的被测电力电缆的宽频阻抗谱,包括:
计算被测电力电缆特性阻抗Z0·E下交叉引线处的等效阻抗Z1·E(f):
根据公式Z′d·E(f)=Zd(f)-Z1·E(f)计算得到Z′d·E(f)为修正后的被测电力电缆的宽频阻抗谱。
3.根据权利要求1所述的电力电缆宽频阻抗谱测试方法,其特征在于,所述高斯脉冲的脉宽w不大于所述矢量网络分析仪设置的扫频范围频率上限值的倒数,且不小于所述矢量网络分析仪设置的扫频范围频率上限值倒数的0.5倍。
4.根据权利要求1所述的电力电缆宽频阻抗谱测试方法,其特征在于,所述进行计算的数据长度N不大于所述矢量网络分析仪设置的扫频点数的0.5倍。
5.根据权利要求1所述的电力电缆宽频阻抗谱测试方法,其特征在于,所述同轴电缆的特性阻抗为50Ω。
6.根据权利要求1所述的电力电缆宽频阻抗谱测试方法,其特征在于,所述屏蔽接线盒的进出线均设置射频连接器,且所述屏蔽接线盒内串联电容。
7.根据权利要求6所述的电力电缆宽频阻抗谱测试方法,其特征在于,所述屏蔽接线盒内串联电容的取值范围为1μF≤C≤10μF。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110541444.6A CN113281614B (zh) | 2021-05-18 | 2021-05-18 | 一种电力电缆宽频阻抗谱测试方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110541444.6A CN113281614B (zh) | 2021-05-18 | 2021-05-18 | 一种电力电缆宽频阻抗谱测试方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113281614A CN113281614A (zh) | 2021-08-20 |
CN113281614B true CN113281614B (zh) | 2023-05-02 |
Family
ID=77279742
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110541444.6A Active CN113281614B (zh) | 2021-05-18 | 2021-05-18 | 一种电力电缆宽频阻抗谱测试方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113281614B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114441606B (zh) * | 2021-12-28 | 2023-11-24 | 国网河北省电力有限公司电力科学研究院 | 电缆水树枝老化缺陷的定位方法及其测试方法、设备 |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100395516B1 (ko) * | 1998-11-19 | 2003-12-18 | 금호석유화학 주식회사 | 비선형등가회로모형을이용한축전장치의특성인자수치화방법및장치 |
CN102608425B (zh) * | 2012-03-26 | 2015-04-22 | 中国科学院物理研究所 | 一种宽频阻抗测量***及宽频阻抗测量方法 |
EP3158655B1 (en) * | 2014-06-19 | 2018-05-30 | Telefonaktiebolaget LM Ericsson (publ) | Estimation of impedances in a transmission medium |
CN204439721U (zh) * | 2015-03-12 | 2015-07-01 | 华北电力大学(保定) | 一种基于网络分析仪的电容分压器宽频阻抗测量*** |
CN105699843A (zh) * | 2016-02-04 | 2016-06-22 | 华中科技大学 | 一种电缆运行状态诊断方法及*** |
CN107561368B (zh) * | 2017-09-08 | 2024-03-19 | 国网浙江省电力有限公司电力科学研究院 | 一种大型电力设备宽频阻抗特性的测量***及测量方法 |
CN107797035A (zh) * | 2017-10-12 | 2018-03-13 | 国网天津市电力公司电力科学研究院 | 基于宽频介电阻抗谱法的xlpe 电缆绝缘性能评估方法 |
CN107942198A (zh) * | 2017-11-20 | 2018-04-20 | 广东核电合营有限公司 | 一种基于阻抗频谱分析的电缆局部缺陷评估的装置和方法 |
CN110161312B (zh) * | 2019-05-07 | 2020-08-25 | 浙江大学 | 基于微带线法的一维和二维材料宽频带阻抗测量装置及方法 |
US10890642B1 (en) * | 2019-07-31 | 2021-01-12 | Keysight Technologies, Inc. | Calibrating impedance measurement device |
CN110794271B (zh) * | 2019-11-18 | 2020-08-18 | 四川大学 | 一种基于输入阻抗谱的电力电缆中间接头受潮定位诊断方法 |
CN112083237B (zh) * | 2020-07-31 | 2021-06-29 | 西安交通大学 | 一种用于大尺度电气设备宽频特性时域测量方法及*** |
CN112763843A (zh) * | 2020-12-21 | 2021-05-07 | 西安交通大学 | 一种基于切比雪夫窗的电缆多段缺陷定位方法及装置 |
-
2021
- 2021-05-18 CN CN202110541444.6A patent/CN113281614B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113281614A (zh) | 2021-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114019309B (zh) | 一种基于频域反射技术的电缆缺陷定位方法 | |
CN110514959B (zh) | 一种考虑电缆衰减特性的电缆故障fdr定位***的定位方法 | |
KR102014582B1 (ko) | 반사파 처리 장치 | |
US20240003954A1 (en) | Long-distance high-voltage cable fault degree detection method and device | |
Sheng et al. | On-line PD detection and localization in cross-bonded HV cable systems | |
CN103630803B (zh) | 一种电缆局部放电带电校正方法 | |
CN105277857A (zh) | 一种在线监测变压器套管受潮缺陷的方法 | |
CN113281614B (zh) | 一种电力电缆宽频阻抗谱测试方法 | |
CN103760402A (zh) | 基于D_dot原理的三相电压互感器及三相影响电压补偿方法 | |
CN106649946B (zh) | 一种输电线路工频相参数仿真计算方法 | |
Shi et al. | Application of augmented spread spectrum time domain reflectometry for detection and localization of soft faults on a coaxial cable | |
CN106990277A (zh) | 时域电压测量装置、测量校准装置及测量校准验证装置 | |
CN113281613B (zh) | 一种电力电缆传递参数获取方法 | |
Hashmi et al. | Covered-conductor overhead distribution line modeling and experimental verification for determining its line characteristics | |
CN115128403A (zh) | 一种用于电力电缆的在线信号测量方法 | |
Graber et al. | Validation of cable models for simulation of transients in shipboard power systems | |
CN109490663B (zh) | 一种晶体谐振器测试***及校准方法 | |
Hashmi et al. | Comparing wave propagation characteristics of MV XLPE cable and covered-conductor overhead line using time domain reflectometry technique | |
Shin et al. | Extraction of propagation constant of XLPE cable using time-domain transmission | |
Shi et al. | A new method of locating the single wire fault | |
CN114113863B (zh) | 高压电缆肘形终端的谐波电压在线测量装置及方法 | |
CN115754611B (zh) | 基于伪梯形波激励的阻抗谱数字重构的电缆故障定位方法 | |
CN117031379B (zh) | 电流传感器高频特性校验电路及其校验方法 | |
Seong et al. | Development of a Fault Localization Algorithm Based on SFWR Considering Propagation Velocity Changes with Propagation Distance | |
CN106918743A (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 |