CN102890193B - 一种多周波采样频率补偿方法 - Google Patents
一种多周波采样频率补偿方法 Download PDFInfo
- Publication number
- CN102890193B CN102890193B CN201110207181.1A CN201110207181A CN102890193B CN 102890193 B CN102890193 B CN 102890193B CN 201110207181 A CN201110207181 A CN 201110207181A CN 102890193 B CN102890193 B CN 102890193B
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- digital filter
- narrow band
- aliasing
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000012360 testing method Methods 0.000 claims abstract description 50
- 238000013461 design Methods 0.000 claims abstract description 17
- 238000000926 separation method Methods 0.000 claims abstract description 15
- 238000005070 sampling Methods 0.000 claims abstract description 4
- 238000001228 spectrum Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 12
- 238000005259 measurement Methods 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 238000010998 test method Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Abstract
本发明提供了一种多周波采样频率补偿方法，包括：对两种频率混叠信号进行采样，采集时间为与工频周期相关的特定时间长度，获得混叠信号的离散数据；设计窄带数字滤波器，获得对混叠信号的离散数据进行分离的窄带数字滤波器；通过设计的窄带数字滤波器，对上述离散数据进行分离，获得频域分离信号。本发明能够实现在信号频域分离过程中，进行准确稳定地异频测试。
Description
技术领域
本发明涉及电力系统输变电设备测试技术领域，特别是涉及一种多周波采样频率补偿方法。
背景技术
高压输电线路正序参数和零序参数以及线路间的互阻抗和耦合电容均是电网重要的系统参数，是电网保护整定，无功补偿，结构优化的重要依据。运行规程(《架空送电线路施工及验收规范》、《110千伏及以上送变电基本建设工程启动验收规程》、《继电保护及安全自动装置运行管理规程》及《DL/T559-94》、《DL/T584-95》等)规定，输电线路在新建，改建后，必须实测参数后才能投入运行。
电网密集程度越来越高，同塔双回或者多回线路成为常态，线路间耦合的紧密程度达到空前的水平，电磁耦合感应的工频干扰电流可以达到100A，高可靠性电能供给的需求使得通过停电来降低干扰水平越来越不容易操作，传统的试验方法中，采用工频电源进行试验的方法不能消除干扰对试验结果的影响，采用偏离工频干扰(50Hz)的异频信号测量线路参数，在频域分离有用的异频测试信号和干扰信号的方法应运而生，异频测试方法在介质损耗测量和大型接地网参数测量中也有广泛应用。
目前现有的技术具有其局限性。异频测试方法的信号处理关键在于频域信号分离，要准确分离测试信号和干扰信号得到稳定的测试结果，干扰信号在测试频点的衰减要达到40db及以上。现有方法采用离散傅里叶变换，采集一个或者两个工频周期被测信号，对信号进行谐波分析，在频域上分离有用的异频测试信号和工频干扰信号。由于异频测试中测试信号频点距离工频(50Hz)频点很近，一般取±2.5Hz或者±5Hz，则必然导致信号数据序列对异频测试信号非整周期采样，虽然可以采用多种方法补足整周期，也必然会导致频域的频谱泄漏，从而导致系统误差。只有在信噪比为2:1的条件下才能达到衰减20db的要求，这在现有干扰环境条件下，导致试验电源系统庞大，而且测试结果数据准确性和稳定性均不能满足测量要求。
总之，需要本领域技术人员迫切解决的一个技术问题就是：如何能够找到一种新型的多周波采样频率补偿方法，能够实现在信号频域分离过程中，进行准确稳定地异频测试。
发明内容
本发明所要解决的技术问题是提供一种多周波采样频率补偿方法，能够实现在信号频域分离过程中，进行准确稳定地异频测试。
为了解决上述问题，本发明公开了一种多周波采样频率补偿方法，包括：
对两种频率混叠信号进行采样，采集时间为与工频周期相关的特定时间长度，获得混叠信号的离散数据；
设计窄带数字滤波器，获得对混叠信号的离散数据进行分离的窄带数字滤波器；
通过设计的窄带数字滤波器，对上述离散数据进行分离，获得频域分离信号。
优选的，所述两种频率分别为测试频率ω1及测试频率ω2。
优选的，所述测试频率ω1及测试频率ω2满足公式(Ⅰ)
|ω2-ω1|＝5π×m (Ⅰ)
其中，m为正整数。
优选的，所述特定时间长度是指时间长度为10K的工频周期；
其中，K为从1到10的正整数。
优选的，所述窄带数字滤波器所对应的傅里叶变化为及
优选的，所述频域分离信号y(n)i为异频信号y(n)1及工频干扰信号y(n)2；
其中，i为1或2。
优选的，所述的多周波采样频率补偿方法还包括：
补偿混叠的异频测试信号为整周期的信息，使得频谱泄漏为零。
优选的，所述的多周波采样频率补偿方法还包括：
利用窄带数字滤波器对被测信息进行滤波处理所采用的方式满足公式(Ⅱ)
y(n)i＝x(n)*h(n)i (Ⅱ)
其中，
y(n)i为频域分离信号，i为1或2，y(n)1为异频信号，y(n)2为工频干扰信号；
x(n)为混叠信号；
h(n)i为窄带数字滤波器，i为1或2。
与现有技术相比，本发明具有以下优点：
本发明提供了一种多周波采样频率补偿方法，具体可以包括：对两种频率混叠信号进行采样，采集时间为与工频周期相关的特定时间长度，获得混叠信号的离散数据；设计窄带数字滤波器，获得对混叠信号的离散数据进行分离的窄带数字滤波器；通过设计的窄带数字滤波器，对上述离散数据进行分离，获得频域分离信号；从而能够实现在信号频域分离过程中，进行准确稳定地异频测试。
附图说明
图1是本发明一种多周波采样频率补偿方法实施例1的流程图；
图2是本发明一种多周波采样频率补偿方法实施例2的流程图；
图3是本发明实施例2中的待测的频率ω1与频率ω2(50Hz)的混叠信号x(n)，混叠比例为1:1时的信号强度分布示意图；
图4是本发明实施例2中的混叠比例1:1时分离的频率为ω1的信号y1的信号强度分布示意图；
图5是本发明实施例2中的混叠比例1:1时分离的频率为ω2的信号y2的信号强度分布示意图；
图6是本发明一种多周波采样频率补偿方法实施例3的流程图；
图7是本发明实施例3中的待测的频率ω1与频率ω2(50Hz)的混叠信号x(n)，混叠比例为1:10时的信号强度分布示意图；
图8是本发明实施例3中的混叠比例1:10时分离的频率为ω1的信号y1的信号强度分布示意图；
图9是本发明实施例3中的混叠比例1:10时分离的频率为ω2的信号y2的信号强度分布示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
本发明的核心构思之一在于，提供了一种多周波采样频率补偿方法，具体可以包括：对两种频率混叠信号进行采样，采集时间为与工频周期相关的特定时间长度，获得混叠信号的离散数据；设计窄带数字滤波器，获得对混叠信号的离散数据进行分离的窄带数字滤波器；通过设计的窄带数字滤波器，对上述离散数据进行分离，获得频域分离信号；从而能够实现在信号频域分离过程中，进行准确稳定地异频测试。
参照图1，示出了本发明一种多周波采样频率补偿方法实施例1的流程图，具体可以包括：
步骤101、对两种频率混叠的信号进行采样，采集时间为与工频周期相关的特定时间长度，获得混叠信号的离散数据。
其中，
所述两种频率分别为测试频率ω1及测试频率ω2。
所述测试频率ω1及测试频率ω2满足公式(Ⅰ)
|ω2-ω1|＝5π×m (Ⅰ)
m为正整数。
所述特定时间长度是指时间长度为10K的工频周期；
K为从1到10的正整数。
步骤102、设计窄带数字滤波器，获得对混叠信号的离散数据进行分离的窄带数字滤波器。
其中，
所述窄带数字滤波器所对应的傅里叶变化为及
步骤103、通过设计的窄带数字滤波器，对上述离散数据进行分离，获得频域分离信号。
其中，
所述频域分离信号y(n)i为异频信号y(n)1及工频干扰信号y(n)2；
i为1或2。
所述通过设计的窄带数字滤波器，对上述工频周期的数据进行分离。
补偿混叠的异频测试信号为整周期，使得频谱泄漏为零。
利用窄带数字滤波器对被测信息进行滤波处理所采用的方式满足公式(Ⅱ)
y(n)i＝x(n)*h(n)i (Ⅱ)
y(n)i为频域分离信号，i为1或2，y(n)1为异频信号，y(n)2为工频干扰信号；
x(n)为混叠信号；
h(n)i为窄带数字滤波器，i为1或2。
参照图2，示出了本发明一种多周波采样频率补偿方法实施例2的流程图，具体可以包括：
步骤201、采集待测的频率ω1与频率ω2(50Hz)的混叠信号x(n)，混叠比例为1:1时的离散数据。
其中，
所述两种频率分别为测试频率ω1及测试频率ω2(50Hz)。
所述测试频率ω1及测试频率ω2满足公式(Ⅰ)
|ω2-ω1|＝5π×m (Ⅰ)
m为正整数。
采集时间长度为10K工频周期(K为从1到10的正整数)，补偿混叠的异频测试信号为整周期，使频谱泄漏理论上为零。
测试频率ω1及测试频率ω2(50Hz)的混叠信号x(n)，混叠比例为1:1。
参照图3，示出了本发明实施例2中的待测的频率ω1与频率ω2(50Hz)的混叠信号x(n)，混叠比例为1:1时的信号强度分布示意图。
从图3中，可以看出：当待测的频率ω1与频率ω2(50Hz)的混叠信号x(n)，混叠比例为1:1时的信号强度情况。
步骤202、设计窄带数字滤波器，获得对待测的频率ω1与频率ω2(50Hz)的混叠信号x(n)，混叠比例为1:1时的离散数据进行分离的窄带数字滤波器。
步骤203、通过设计的窄带数字滤波器，对上述离散数据进行分离，获得频域分离信号。
参照图4，示出了本发明实施例2中的混叠比例1:1时分离的频率为ω1的信号y1的信号强度分布示意图。
从图4中，可以看出：混叠比例1:1时分离的频率为ω1的信号y1的信号强度分布情况。
参照图5，示出了本发明实施例2中的混叠比例1:1时分离的频率为ω2的信号y2的信号强度分布示意图。
从图5中，可以看出：混叠比例1:1时分离的频率为ω2的信号y2的信号强度分布情况。
参照图6，示出了本发明一种多周波采样频率补偿方法实施例3的流程图，具体可以包括：
步骤601、通过两种频率，采集待测的频率ω1与频率ω2(50Hz)的混叠信号x(n)，混叠比例为1:10时的离散数据。
其中，
所述两个测试频率分别为测试频率ω1及测试频率ω2(50Hz)。
所述测试频率ω1及测试频率ω2满足公式(Ⅰ)
|ω2-ω1|＝5π×m (Ⅰ)
m为正整数。
采集时间长度为10K工频周期(K为从1到10的正整数)，补偿混叠的异频测试信号为整周期，使频谱泄漏理论上为零。
测试频率ω1及测试频率ω2(50Hz)的混叠信号x(n)，混叠比例为1:10。
参照图7，示出了本发明实施例3中的待测的频率ω1与频率ω2(50Hz)的混叠信号x(n)，混叠比例为1:10时的信号强度分布示意图。
从图7中，可以看出：当待测的频率ω1与频率ω2(50Hz)的混叠信号x(n)，混叠比例为1:10时的信号强度分布情况。
步骤602、设计窄带数字滤波器，获得对待测的频率ω1与频率ω2(50Hz)的混叠信号x(n)，混叠比例为1:10时的离散数据进行分离的窄带数字滤波器。
步骤603、通过设计的窄带数字滤波器，对上述离散数据进行分离，获得频域分离信号。
参照图8，示出了本发明实施例3中的混叠比例1:10时分离的频率为ω1的信号y1的信号强度分布示意图。
从图8中，可以看出：混叠比例1:10时分离的频率为ω1的信号y1的信号强度分布情况。
参照图9，示出了本发明实施例3中的混叠比例1:10时分离的频率为ω2的信号y2的信号强度分布示意图。
从图9中，可以看出：混叠比例1:10时分离的频率为ω2的信号y2的信号强度分布情况。
本发明提供了一种多周波采样频率补偿方法。通过本发明所提供的方法可以提高信号分离能力，能在异频测试信号和工频干扰信号比值为1:10的条件下进行信号频域分离；提高信号分离分辨率，可分离的频偏可以达到±1Hz；工频干扰信号在测试频点的衰减大于等于40db，提高异频测试准确度；信号分离能力的提高可以大大减小异频试验电源的体积和质量，提高了自动化程度。
总之，本发明提供了一种多周波采样频率补偿方法，能够实现在信号频域分离过程中，进行准确稳定地异频测试。
本说明书中的各个实施例均采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。
以上对本发明所提供的一种多周波采样频率补偿方法，进行了详细介绍，本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
Claims (4)
1.一种多周波采样频率补偿方法，其特征在于，包括：
对两种频率混叠信号进行采样，采集时间为与工频周期相关的特定时间长度，获得混叠信号的离散数据；采集时间长度为10K工频周期，K为从1到10的正整数，补偿混叠信号为整周期；
设计窄带数字滤波器，获得对混叠信号的离散数据进行分离的窄带数字滤波器；所述窄带数字滤波器所对应的傅里叶变化为及
通过设计的窄带数字滤波器，对上述离散数据进行分离，获得频域分离信号；
利用窄带数字滤波器对离散数据进行滤波处理所采用的方式满足公式(Ⅱ)
y(n)i＝x(n)*h(n)i (Ⅱ)
其中，y(n)i为频域分离信号，i为1或2，y(n)1为异频信号，y(n)2为工频干扰信号；x(n)为混叠信号；h(n)i为窄带数字滤波器，i为1或2；
所述频域分离信号y(n)i为异频信号y(n)1或工频干扰信号y(n)2；其中，i为1或2。
2.根据权利要求1所述的方法，其特征在于：所述两种频率分别为测试频率ω1及测试频率ω2。
3.根据权利要求2所述的方法，其特征在于：
所述测试频率ω1及测试频率ω2满足公式(Ⅰ)
|ω2-ω1|＝5π×m (Ⅰ)
其中，m为正整数。
4.如权利要求1所述的多周波采样频率补偿方法，其特征在于，还包括：
补偿混叠的异频测试信号为整周期的信息，使得频谱泄漏为零。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110207181.1A CN102890193B (zh) | 2011-07-22 | 一种多周波采样频率补偿方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110207181.1A CN102890193B (zh) | 2011-07-22 | 一种多周波采样频率补偿方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102890193A CN102890193A (zh) | 2013-01-23 |
| CN102890193B true CN102890193B (zh) | 2016-12-14 |
Family
ID=
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4093989A (en) * | 1976-12-03 | 1978-06-06 | Rockland Systems Corporation | Spectrum analyzer using digital filters |
| CN1243254A (zh) * | 1999-08-16 | 2000-02-02 | 国家电力公司武汉高压研究所 | 接地电阻异频测量方法及装置 |
| CN1595065A (zh) * | 2004-07-08 | 2005-03-16 | 赵岳生 | 基于矢量调制法的惯性量测量单元 |
| CN1707613A (zh) * | 2005-05-27 | 2005-12-14 | 曲哲 | 隔噪声音频采集装置及方法 |
| EP1657777A1 (en) * | 2004-11-15 | 2006-05-17 | TDK Corporation | High-frequency module |
| CN101034117A (zh) * | 2007-04-02 | 2007-09-12 | 武汉大学 | 输电线路工频参数测量方法及装置 |
| CN101505060A (zh) * | 2009-03-18 | 2009-08-12 | 天津大学 | 电网中利用窄带滤波实现多采样率信号的频率转换方法 |
| CN102055435A (zh) * | 2010-12-23 | 2011-05-11 | 中国科学院武汉物理与数学研究所 | 一种窄带数字滤波器 |
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4093989A (en) * | 1976-12-03 | 1978-06-06 | Rockland Systems Corporation | Spectrum analyzer using digital filters |
| CN1243254A (zh) * | 1999-08-16 | 2000-02-02 | 国家电力公司武汉高压研究所 | 接地电阻异频测量方法及装置 |
| CN1595065A (zh) * | 2004-07-08 | 2005-03-16 | 赵岳生 | 基于矢量调制法的惯性量测量单元 |
| EP1657777A1 (en) * | 2004-11-15 | 2006-05-17 | TDK Corporation | High-frequency module |
| CN1707613A (zh) * | 2005-05-27 | 2005-12-14 | 曲哲 | 隔噪声音频采集装置及方法 |
| CN101034117A (zh) * | 2007-04-02 | 2007-09-12 | 武汉大学 | 输电线路工频参数测量方法及装置 |
| CN101505060A (zh) * | 2009-03-18 | 2009-08-12 | 天津大学 | 电网中利用窄带滤波实现多采样率信号的频率转换方法 |
| CN102055435A (zh) * | 2010-12-23 | 2011-05-11 | 中国科学院武汉物理与数学研究所 | 一种窄带数字滤波器 |
Non-Patent Citations (2)
| Title |
|---|
| 基于异频法的大型接地网接地电阻测试研究;马明等;《浙江电力》;20080831(第4期);第57页左栏第3段 * |
| 异频线路参数测试方法的应用研究;吴智华;《广西电力》;20090228(第1期);第40-43页 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110208662A (zh) | 基于介质谱的超导电缆pplp绝缘检测方法及系统 | |
| WO2015117466A1 (zh) | 一种抗干扰特高压线路参数集成测试装置 | |
| CN102890226B (zh) | 电力系统xlpe电缆水树老化状态测试系统 | |
| Schichler et al. | UHF partial discharge detection system for GIS: Application guide for sensitivity verification | |
| CN106707007B (zh) | 变电设备绝缘在线监测系统的工作方法 | |
| CN106093591B (zh) | 一种中性点不接地配电网电容电流测量系统及方法 | |
| EP3338099B1 (de) | Verfahren und messaufbau zur messung von teilentladungen in steckverbindungen und verwendung für steckverbindungen eines schienenfahzeugs | |
| CN109254234A (zh) | 一种树木-导线放电模拟试验方法 | |
| CN104237832B (zh) | 一种复阻抗标准器的校准方法及装置 | |
| CN105353332B (zh) | 一种电子式互感器长期带电性能考核方法及系统 | |
| Zhang et al. | Fault location method for three-core cable using amplitude ratio of shield-grounding wire currents | |
| CN110161375A (zh) | 一种基于分布电阻参数的高压直流输电线路计算模型 | |
| CN100487469C (zh) | 基于行波时差的输电网故障定位方法及装置 | |
| Zhao et al. | Testing and modelling of voltage transformer for high order harmonic measurement | |
| Lai et al. | Localization with phase resolved partial discharge measured with high frequency current transformer | |
| CN102890193B (zh) | 一种多周波采样频率补偿方法 | |
| CN105067851B (zh) | 一种用于电力系统带电测试仪器校验用的信号发生装置 | |
| CN109270404A (zh) | 一种电压行波精确检测方法和装置 | |
| CN2682427Y (zh) | 基于行波时差的输电网故障定位装置 | |
| Frey et al. | Study of high frequency harmonics propagation in industrial networks | |
| CN105866667B (zh) | 一种双端接地断路器的回路电阻测试方法 | |
| Fasil et al. | Modeling and simulation based study for on-line detection of partial discharge of solid dielectric | |
| Zhang et al. | A novelty digital algorithm for online measurement of dielectric loss factor of electronic transformers | |
| Polley et al. | OLPD condition monitoring of complete networks for oil and gas: Challenges and solutions explained through case studies | |
| CN210639238U (zh) | 一种通过同时测量电压和电流进行高压容性监测的设备 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| CB02 | Change of applicant information |
Address after: 430077, 227 East Main Street, Wuchang District, Hubei, Wuhan Applicant after: Hubei Electric Power Research Institute of Power Company Address before: 430077 No. 227 East Main Street, Hubei, Wuhan Applicant before: Hubei Electric Power Testing & Research Institute, Hubei Electric Power Company |
|
| TA01 | Transfer of patent application right |
Effective date of registration: 20131212 Address after: 430077, 227 East Main Street, Wuchang District, Hubei, Wuhan Applicant after: Hubei Electric Power Research Institute of Power Company Applicant after: State Grid Corporation of China Address before: 430077, 227 East Main Street, Wuchang District, Hubei, Wuhan Applicant before: Hubei Electric Power Research Institute of Power Company |
|
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant |