CN109783865B - 一种三芯电力线单位长度参数矩阵求解方法 - Google Patents

一种三芯电力线单位长度参数矩阵求解方法 Download PDF

Info

Publication number
CN109783865B
CN109783865B CN201811533983.XA CN201811533983A CN109783865B CN 109783865 B CN109783865 B CN 109783865B CN 201811533983 A CN201811533983 A CN 201811533983A CN 109783865 B CN109783865 B CN 109783865B
Authority
CN
China
Prior art keywords
unit length
mode
differential mode
input impedance
capacitance
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
Application number
CN201811533983.XA
Other languages
English (en)
Other versions
CN109783865A (zh
Inventor
郭以贺
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
Original Assignee
North China Electric Power University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by North China Electric Power University filed Critical North China Electric Power University
Priority to CN201811533983.XA priority Critical patent/CN109783865B/zh
Publication of CN109783865A publication Critical patent/CN109783865A/zh
Application granted granted Critical
Publication of CN109783865B publication Critical patent/CN109783865B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

本发明提出了一种三芯电力线单位长度参数矩阵求解方法。由于三芯电力线的单位长度参数矩阵元素无法直接得到准确值,基于电磁仿真方法可以得到单位长度参数的初值,由于直接进行修正的参数较多,每个参数会对多种方式下的信道测量结果产生影响,提出一种基于模域的参数修正方法,使得各个参数的修正能独立完成。即基于相模变换方法,将三芯电力线的单位长度参数矩阵解耦为共模和差模参数,在电磁仿真的基础上,通过测量两种模式下的开路阻抗,对模域单位长度参数进行线性修正,最终得到单位长度参数矩阵的准确值。

Description

一种三芯电力线单位长度参数矩阵求解方法
技术领域
本发明涉及通信领域,尤其涉及一种三芯电力线单位长度参数矩阵求解方法。
背景技术
针对三芯电力线缆,采用多输入多输出(Multiple-Input Multiple-Output,MIMO)传输方式可以进一步提高信道容量。在MIMO-PLC技术领域,国内外学者进行了大量的测试和研究工作,基于自底向上方法能够准确分析芯线间的耦合效应,能用于信道特性的预测,该类方法基于传输线理论和网络拓扑结构计算信道传输和阻抗特性,难点在于电力线单位长度参数矩阵的准确获取。
从目前国内外相关研究工作看,单位长度参数主要获取方法包括:近似公式计算法,基于测量的方法,基于优化的方法。其中,采用近似公式计算单位长度参数,由于未能准确计及邻近效应和绝缘介质的影响,计算结果与测量结果误差较大。基于测量的方法计算单位长度参数时,由于开路和短路输入阻抗的极值在频域交替变化,受制于测量结果的精度,两者相乘或相除的结果会有较大的分散性,由此计算得到的单位长度参数误差较大。基于优化的方法计算单位长度参数时,基于全部参数的最小误差作为优化目标,导致计算的复杂度较高。基于电磁仿真方法得到单位长度参数,考虑了邻近效应的影响,能得到准确的单位长度电阻数值和电感数值。但是,由于不同电力线绝缘材料的高频特性也往往各不相同,仍然无法得到准确的单位长度电容数值和电导数值。
目前,三芯电力线单位长度参数矩阵准确获取的问题还未得到很好的解决。
发明内容
为解决上述技术问题,本发明提出了一种三芯电力线单位长度参数矩阵求解方法。
本发明实施例提供了如下技术方案:
一种三芯电力线单位长度参数矩阵求解方法,包括:
根据待测的三芯电力线的每根导体的半径、间距以及周围绝缘介质的横截面尺寸,基于电磁场仿真工具提取单位长度电阻矩阵R0、电感矩阵L0、电容矩阵C0、电导矩阵G0,其中,
Figure BDA0001906408190000021
Figure BDA0001906408190000022
4个单位长度参数矩阵有8个参数,其中:Gm和Cm为负数;
设信号的角频率为ω,则单位长度串联阻抗矩阵为Z0=R0+jωL0,并联导纳矩阵为Y0=G0+jωC0,计算相模变换矩阵T,利用相模变换矩阵经过相模变换,在模域可以分解为差模和共模两个相互独立的传播模式,对于差模来说,4个单位长度参数分别为Rs-Rm、Ls-Lm、Gs-Gm和Cs-Cm,对于共模来说,4个单位长度参数分别为Rs+Rm、Ls+Lm、Gs+Gm和Cs+Cm,计算得到末端开路下始端的差模和共模输入阻抗;
用测量的方法得到末端开路下始端的共模和差模输入阻抗;
基于计算得到的末端开路下始端的差模和共模输入阻抗和测量得到的末端开路下始端的差模和共模输入阻抗,利用线性修正方法对共模和差模单位长度电容C和电导G进行求解,得到共模和差模修正后的单位长度电容C和电导G;
根据得到共模和差模修正后的单位长度电容C和电导G,得到实际的单位长度参数矩阵。
其中,用测量的方法得到末端开路下始端的共模和差模输入阻抗,具体包括:
为了测量共模阻抗,将待测的三芯电力线始端两边的芯线短路,将网络分析仪的端口1接到两个边芯和中芯之间,线缆的末端开路,根据测量得到的S11,测量得到的末端开路下始端的共模输入阻抗为:
Figure BDA0001906408190000031
为了测量差模阻抗,将网络分析仪的端口1接到待测的三芯电力线始端两边的芯导体间,线缆的末端开路,根据测量得到的S′11,测量得到的末端开路下始端的差模输入阻抗为:
Figure BDA0001906408190000032
其中,基于计算得到的末端开路下始端的差模和共模输入阻抗和测量得到的末端开路下始端的差模和共模输入阻抗,利用线性修正方法对共模和差模单位长度电容C和电导G进行求解,得到共模和差模修正后的单位长度电容C和电导G,具体包括:
设电导G=0,利用线性修正方法对单位长度电容C进行求解,得到修正后的单位长度电容C的数值,利用该步骤,分别针对差模和共模单位长度参数初值,进行C的修正,分别得到差模和共模修正后的单位长度电容C的数值;
在得到修正后的单位长度电容C的数值后,利用线性修正方法对单位长度电导G进行求解,得到修正后的单位长度电导G的数值,利用该步骤,分别针对差模和共模单位长度参数初值,进行G的修正,分别得到差模和共模修正后的单位长度电导G的数值。
其中,设电导G=0,利用线性修正方法对单位长度电容C进行求解,得到修正后的单位长度电容C的数值,利用该步骤,分别针对差模和共模单位长度参数初值,进行C的修正,分别得到差模和共模修正后的单位长度电容C的数值,具体包括:
设实际电容
Figure BDA0001906408190000041
式中:f为频率,fend电磁仿真的最大频率值,m和n为待定系数;
根据频域中前几个周期的误差,修正m使计算得到的末端开路下始端的差模输入阻抗和测量得到的末端开路下始端的差模输入阻抗的振荡周期能较好吻合,如果在高频段吻合度不够理想,修正n使得整个频段能满足精度要求,通过线性修正,使得整个频段内计算得到的末端开路下始端的差模输入阻抗和测量得到的末端开路下始端的差模输入阻抗振荡周期基本一致,得到差模修正后的单位长度电容C的数值;
根据频域中前几个周期的误差,修正m使计算得到的末端开路下始端的共模输入阻抗和测量得到的末端开路下始端的共模输入阻抗的振荡周期能较好吻合,如果在高频段吻合度不够理想,修正n使得整个频段能满足精度要求,通过线性修正,使得整个频段内计算得到的末端开路下始端的共模输入阻抗和测量得到的末端开路下始端的共模输入阻抗振荡周期基本一致,得到共模修正后的单位长度电容C的数值。
其中,在得到修正后的单位长度电容C的数值后,利用线性修正方法对单位长度电导G进行求解,得到修正后的单位长度电导G的数值,利用该步骤,分别针对差模和共模单位长度参数初值,进行G的修正,分别得到差模和共模修正后的单位长度电导G的数值,具体包括:
调整tan(θ)的数值使得计算得到的末端开路下始端的差模输入阻抗和测量得到的末端开路下始端的差模输入阻抗的振荡幅度基本一致,得到差模修正后的单位长度电导G的数值;
调整tan(θ)的数值使得计算得到的末端开路下始端的共模输入阻抗和测量得到的末端开路下始端的共模输入阻抗的振荡幅度基本一致,得到共模修正后的单位长度电导G的数值。
与现有技术相比,上述技术方案具有以下优点:
本发明实施例所提供的方法,本发明提出了一种三芯电力线单位长度参数矩阵求解方法,该方法基于电磁仿真方法可以得到单位长度参数的初值,由于直接进行修正的参数较多,每个参数会对多种方式下的信道测量结果产生影响,提出一种基于模域的参数修正方法,使得各个参数的修正能独立完成。即基于相模变换方法,将三芯电力线的单位长度参数矩阵解耦为共模和差模参数,在电磁仿真的基础上,通过测量两种模式下的开路阻抗,对模域单位长度参数进行线性修正,最终得到单位长度参数矩阵的准确值。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为一种三芯电力线单位长度参数矩阵求解方法的示意图;
图2为三芯电力线截面示意图;
图3为差模和共模两个相互独立传播模式的示意图;
图4为经过参数修正后末端开路下始端的差模和共模输入阻抗的计算结果与测量结果的对比的示意图。
具体实施方式
正如背景技术中所述,如何得到准确的三芯线缆单位长度参数矩阵,是本领域技术人员急需解决的技术问题。
为了解决上述技术问题,本发明提出了一种三芯电力线单位长度参数矩阵求解方法,本发明的核心思想是,由于三芯电力线的单位长度参数矩阵元素无法直接得到准确值,基于电磁仿真方法可以得到单位长度参数的初值,由于直接进行修正的参数较多,每个参数会对多种方式下的信道测量结果产生影响,提出一种基于模域的参数修正方法,使得各个参数的修正能独立完成。即基于相模变换方法,将三芯电力线的单位长度参数矩阵解耦为共模和差模参数,在电磁仿真的基础上,通过测量两种模式下的开路阻抗,对模域单位长度参数进行线性修正,最终得到单位长度参数矩阵的准确值。
其中,开路阻抗即为末端开路下始端的输入阻抗。
首先,介绍一下电容C和电导G,具体如下:
实际中低压电力线常使用PVC等材料作为绝缘层和护套,其介电常数受材料特性、加工工艺和填充物的影响,而且是频率的弱相关函数。与电容并联的电导,同样与绝缘材料有关。实际电力线绝缘材料主要考虑50Hz下的绝缘性能,随着频率的升高,介质中将流过更大的由偶极子运动形成的漏电流。为了全面描述绝缘材料的这两种电气特性,需要引入复介电常数
εr=ε′r-iε″r
其中,ε′r为与电容值有关的介电常数,ε″r代表损耗相关部分。
实际中,多采用损耗角正切tan(θ)描述电导G
Figure BDA0001906408190000071
G=ωtan(θ)C
其中,ω为信号的角频率,ω=2πf,f为频率,f的单位为Hz。
图1是一种三芯电力线单位长度参数矩阵求解方法的示意图,所述方法包括:
步骤101:根据待测的三芯电力线的每根导体的半径、间距以及周围绝缘介质的横截面尺寸,基于电磁场仿真工具提取单位长度电阻矩阵R0、电感矩阵L0、电容矩阵C0、电导矩阵G0
其中,
Figure BDA0001906408190000072
4个单位长度参数矩阵有8个参数,其中:Gm和Cm为负数,
其中,待测的三芯电力线为对称的三芯电力线。
以图2的三芯电力线为例,以中间导体为参考地,线缆呈对称结构,三芯线构成2+1传输线系统,为描述方便将三芯导体从左到右依次标记为a、b和c。电磁仿真时,根据PVC的典型特性,预设其介电常数ε′r为3,损耗角正切tan(θ)为0.04。
其中,电磁场仿真工具有多种,例如ANSYS的电磁场仿真工具,本发明中可以用ANSYS的电磁场仿真工具来提取单位长度电阻矩阵和电感矩阵,当然可以用其它的公司的电磁场仿真工具提取单位长度电阻矩阵和电感矩阵,在此并不进行限定。由于单位长度电阻和电感数值不受绝缘介质参数的影响,因此基于电磁仿真得到的电阻矩阵和电感矩阵即为准确值,电导矩阵和电容矩阵需要进行修正。
步骤102:设信号的角频率为ω,则单位长度串联阻抗矩阵为Z0=R0+jωL0,并联导纳矩阵为Y0=G0+jωC0,对于待测的三芯电力线,计算相模变换矩阵T,利用相模变换矩阵经过相模变换,在模域可以分解为差模和共模两个相互独立的传播模式,计算得到末端开路下始端的差模和共模输入阻抗。
具体过程如下:
设信号的角频率为ω,则单位长度串联阻抗矩阵为Z0=R0+jωL0,并联导纳矩阵为Y0=G0+jωC0,即
Figure BDA0001906408190000081
Figure BDA0001906408190000082
计算相模变换矩阵T为:
Figure BDA0001906408190000083
经过相模变换
Figure BDA0001906408190000084
在模域可以分解为差模和共模两个相互独立的传播模式,例如图3所示。
其中,差模分量通过两边的导体传播,两个导体的电流大小相等,方向相反。共模分量通过中间的参考导体回流,两边导体的电流大小相等,方向相同。
对于差模来说,4个单位长度参数分别为Rs-Rm、Ls-Lm、Gs-Gm和Cs-Cm;对于共模来说,4个单位长度参数分别为Rs+Rm、Ls+Lm、Gs+Gm和Cs+Cm。两组参数构成两个相互独立的双导体传输线系统,通过精确求解两组单位长度参数后,就可以反过来计算得到实际的单位长度参数矩阵。
由于一段均匀线缆的模域输入阻抗呈现双导体传输线的典型振荡特性,在此基础上可以采用简单的线性拟合方法进行参数修正。下面针对双导体传输线系统说明参数修正方法,该方法同时适用于差模和共模两种情况。在修正过程中,基于电力线的低损耗特性进行单位长度参数的影响分析,将单位长度电容和电导对末端开路下始端的输入阻抗的影响分离,即分别对单位长度电容C和电导G进行修正,以便简化单位长度参数求解的过程。
设均匀传输线的单位长度参数分别为R、L、G和C,进一步得到传输线的二次参数。其中,传播常数为
Figure BDA0001906408190000091
特性阻抗为
Figure BDA0001906408190000092
设传输线的长度为lb,则末端开路下始端的输入阻抗为
Figure BDA0001906408190000093
在15MHz以下频段,一般低压电力线缆满足低损耗线条件:R<<ωL和G<<ωC。可得到相应的近似结果为
Figure BDA0001906408190000094
Figure BDA0001906408190000095
则计算得到的末端开路下始端的输入阻抗Zino可简化为
Figure BDA0001906408190000101
分别用差模和共模的4个单位长度参数代入上式,即可以得到末端开路下始端的差模和共模输入阻抗的计算结果。
步骤103:用测量的方法得到末端开路下始端的共模和差模输入阻抗,具体包括:
为了测量共模阻抗,将待测的三芯电力线始端两边的芯线a和c短路,将网络分析仪的端口1接到两个边芯和中芯b之间,线缆的末端开路,根据测量得到的S11,测量得到的末端开路下始端的共模输入阻抗为:
Figure BDA0001906408190000102
为了测量差模阻抗,将网络分析仪的端口1接到待测的三芯电力线始端的a和c芯导体间,线缆的末端开路,根据测量得到的S′11,测量得到的末端开路下始端的差模输入阻抗为:
Figure BDA0001906408190000103
步骤104:基于计算得到的末端开路下始端的差模和共模输入阻抗和测量得到的末端开路下始端的差模和共模输入阻抗,利用线性修正方法对共模和差模单位长度电容C和电导G进行求解,得到共模和差模修正后的单位长度电容C和电导G。
步骤1041:设电导G=0,利用线性修正方法对单位长度电容C进行求解,得到修正后的单位长度电容C,利用该步骤,分别针对差模和共模单位长度参数初值,进行C的修正,分别得到差模和共模修正后的单位长度电容C的数值;
步骤1041具体过程包括:
设实际电容
Figure BDA0001906408190000111
式中:f为频率,fend电磁仿真的最大频率值,这个值可以由用户来设置,例如fend=15MHz,m和n为待定系数。
根据频域中前几个周期的误差,修正m使计算得到的末端开路下始端的差模输入阻抗和测量得到的末端开路下始端的差模输入阻抗的振荡周期能较好吻合,如果在高频段吻合度不够理想,修正n使得整个频段能满足精度要求,通过线性修正,使得整个频段内计算得到的末端开路下始端的差模输入阻抗和测量得到的末端开路下始端的差模输入阻抗振荡周期基本一致,得到差模修正后的单位长度电容C的数值;
根据频域中前几个周期的误差,修正m使计算得到的末端开路下始端的共模输入阻抗和测量得到的末端开路下始端的共模输入阻抗的振荡周期能较好吻合,如果在高频段吻合度不够理想,修正n使得整个频段能满足精度要求,通过线性修正,使得整个频段内计算得到的末端开路下始端的共模输入阻抗和测量得到的末端开路下始端的共模输入阻抗振荡周期基本一致,得到共模修正后的单位长度电容C的数值。
其中,在对单位长度电容C进行求解时,将设tan(θ)=0,即设电导G=0,先忽略电导G,将单位长度电容和电导对末端开路下始端的输入阻抗的影响分离,即分别对单位长度电容C和电导G进行修正,以便简化单位长度参数求解的过程。
在频域,Zino的振荡周期与
Figure BDA0001906408190000121
的周期相同,由一次参数L、C和线路长度lb决定。利用电磁仿真得到的L和已知的长度lb,对比计算得到的末端开路下始端的输入阻抗和测量得到的末端开路下始端的输入阻抗振荡周期,即可修正C至准确数值。
随着频率的变化,当
Figure BDA0001906408190000122
时,Zino为峰值,当
Figure BDA0001906408190000123
时,Zino为谷值。这两种情况下Zino的相位为0°,呈阻性。当
Figure BDA0001906408190000124
时,Zino=ZC,相位达到极值。末端开路下始端的输入阻抗的振荡周期与
Figure BDA0001906408190000125
的周期相同,由一次参数L、C和线路长度lb决定。当lb已知,L通过仿真得到后,可以根据Zino得到C的具体数值。Zino的振荡周期与L和C的乘积有关,而末端开路下始端的输入阻抗均值ZC由L和C的比值决定。当求得的C使振荡周期吻合时,也会同时使末端开路下始端的输入阻抗的均值具有较好的吻合度,这样进一步验证了L和C的准确性。即C只影响末端开路下始端的输入阻抗的振荡周期和均值,据此,可以先修正C至准确值,即得到修正后的单位长度电容C的数值。
步骤1042:在得到修正后的单位长度电容C的数值后,利用线性修正方法对单位长度电导G进行求解,得到修正后的单位长度电导G的数值,利用该步骤,分别针对差模和共模单位长度参数初值,进行G的修正,分别得到差模和共模修正后的单位长度电导G的数值,具体过程包括:
调整tan(θ)的数值使得计算得到的末端开路下始端的差模输入阻抗和测量得到的末端开路下始端的差模输入阻抗的振荡幅度基本一致,得到差模修正后的单位长度电导G的数值;
调整tan(θ)的数值使得计算得到的末端开路下始端的共模输入阻抗和测量得到的末端开路下始端的共模输入阻抗的振荡幅度基本一致,得到共模修正后的单位长度电导G的数值。
Zino振荡衰减的幅度由
Figure BDA0001906408190000131
决定,受R和G频变性的影响,随着频率的升高,Zino的振荡幅度越来越小。由于R基本上与频率的平方根成正比,G则随频率线性增加。当频率较低时,G的数值很小,R的影响占主要地位,由此可以根据低频段的输入阻抗数值验证R。在高频段,G对Zino振荡幅度的影响变大,根据Zino与实测曲线的差异,修正tan(θ)的数值,从而得到G的准确数值。
图4为经过参数修正后末端开路下始端的差模和共模输入阻抗的计算结果与测量结果的对比,两种模式下,模值和相角的吻合度都很好,说明采用线性修正就可以较好地满足模型的精度要求,从而避免了使用更为复杂的非线性方法。修正后,共模参数为Rc、Lc、Gc和Cc,差模参数为Rd、Ld、Gd和Cd
步骤105:根据得到共模和差模修正后的单位长度电容C和电导G,得实际的单位长度参数矩阵为
Figure BDA0001906408190000132
Figure BDA0001906408190000133
Figure BDA0001906408190000134
Figure BDA0001906408190000135
可见,本发明所述的上述方法,利用相模变换,将三芯电力线单位长度参数矩阵分解为共模和差模参数,将三导体传输线系统转化为双导体传输线系统,基于双导体传输线系统的开路阻抗特性进行单位长度电容和电导的修正,避免了参数同时修正带来的复杂性。即基于电磁仿真得到单位长度电阻和电感的准确值,得到单位长度电容和电导的合理估计值,进一步结合末端开路下始端的输入阻抗测量结果对估计值进行修正。可见,本发明融合了电磁仿真计算和辅助测量两种单位长度参数求解方法,利用辅助测量结果对电磁仿真得到的单位长度电容和电导进行线性修正,从而得到准确的单位长度电容和电导,避免了复杂非线性方法的使用。
而且,本发明所述的方法,根据相模变换矩阵的数值,在三芯电力线的始端通过不同的芯线连接方式,分别构建共模和差模阻抗测试条件并进行对应测量的方法。
本发明在修正过程中,基于电力线在30kHz-15MHz频带范围内的低损耗特性进行单位长度参数的影响分析,将单位长度电容和电导对末端开路下始端的输入阻抗的影响分离,即依次进行单位长度电容的修正和电导的修正,避免了单位长度电容和电导同时修正带来的复杂性,简化了单位长度参数修正和求解的过程。该方法能有效避免直接测量方法导致的数据发散问题,由于拟合的参数数目较少,参数修正范围较小,可以采用简单的线性拟合方法,避免了非线性优化方法所需的复杂计算和容易导致的非因果性问题。
而且,本发明所述的方法以末端开路下始端的输入阻抗模值和相角的测量结果为参考,进行单位长度参数求解和验证,该方法简单易行,可信度高。
本发明所述的方法基于相模变换矩阵的数值,将单位长度参数的修正转化为共模和差模两种情况下分别进行,避免了参数同时修正带来的复杂性,单位长度参数矩阵的求解步骤清晰。
本发明所述的方法以电磁仿真得到的单位长度参数矩阵作为初值进行参数修正,使得修正的范围较小,利用相对简单的线性方法就可以得到精确值。
本发明所述的方法相对于与其他方法相比,通过本方法得到的单位长度参数矩阵的精度有了明显提高。
本发明所述的方法利用常见三芯电力线结构对称的特点,不需要预先测量线缆绝缘介质的介电常数,单位长度参数矩阵的求解方法具有广泛适用性。
本说明书中各个部分采用递进的方式描述,每个部分重点说明的都是与其他部分的不同之处,各个部分之间相同相似部分互相参见即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修正对本领域的专业技术人员来说将是显而易见的,本申请中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本申请所示的实施例,而是要符合与本申请所公开的原理和新颖特点相一致的最宽的范围。

Claims (5)

1.一种三芯电力线单位长度参数矩阵求解方法,其特征在于,包括:
根据待测的三芯电力线的每根导体的半径、间距以及周围绝缘介质的横截面尺寸,基于电磁场仿真工具提取单位长度电阻矩阵R0、电感矩阵L0、电容矩阵C0、电导矩阵G0,其中,
Figure FDA0004057062320000011
Figure FDA0004057062320000012
4个单位长度参数矩阵有8个参数,其中:Gm和Cm为负数;
设信号的角频率为ω,则单位长度串联阻抗矩阵为Z0=R0+jωL0,并联导纳矩阵为Y0=G0+jωC0,计算相模变换矩阵T,利用相模变换矩阵经过相模变换,在模域可以分解为差模和共模两个相互独立的传播模式,对于差模来说,4个单位长度参数分别为Rs-Rm、Ls-Lm、Gs-Gm和Cs-Cm,对于共模来说,4个单位长度参数分别为Rs+Rm、Ls+Lm、Gs+Gm和Cs+Cm,计算得到末端开路下始端的差模和共模输入阻抗;
用测量的方法得到末端开路下始端的共模和差模输入阻抗;
基于计算得到的末端开路下始端的差模和共模输入阻抗和测量得到的末端开路下始端的差模和共模输入阻抗,利用线性修正方法对共模和差模单位长度电容C和电导G进行求解,得到共模和差模修正后的单位长度电容C和电导G;
根据得到共模和差模修正后的单位长度电容C和电导G,得到实际的单位长度参数矩阵。
2.根据权利要求1所述的方法,其特征在于,用测量的方法得到末端开路下始端的共模和差模输入阻抗,具体包括:
为了测量共模阻抗,将待测的三芯电力线始端两边的芯线短路,将网络分析仪的端口1接到两个边芯和中芯之间,线缆的末端开路,根据测量得到的S11,测量得到的末端开路下始端的共模输入阻抗为:
Figure FDA0004057062320000021
为了测量差模阻抗,将网络分析仪的端口1接到待测的三芯电力线始端两边的芯导体间,线缆的末端开路,根据测量得到的S1'1,测量得到的末端开路下始端的差模输入阻抗为:
Figure FDA0004057062320000022
3.根据权利要求1所述的方法,其特征在于,基于计算得到的末端开路下始端的差模和共模输入阻抗和测量得到的末端开路下始端的差模和共模输入阻抗,利用线性修正方法对共模和差模单位长度电容C和电导G进行求解,得到共模和差模修正后的单位长度电容C和电导G,具体包括:
设电导G=0,利用线性修正方法对单位长度电容C进行求解,得到修正后的单位长度电容C的数值,利用该步骤,分别针对差模和共模单位长度参数初值,进行C的修正,分别得到差模和共模修正后的单位长度电容C的数值;
在得到修正后的单位长度电容C的数值后,利用线性修正方法对单位长度电导G进行求解,得到修正后的单位长度电导G的数值,利用该步骤,分别针对差模和共模单位长度参数初值,进行G的修正,分别得到差模和共模修正后的单位长度电容G的数值。
4.根据权利要求3所述的方法,其特征在于,设电导G=0,利用线性修正方法对单位长度电容C进行求解,得到修正后的单位长度电容C的数值,利用该步骤,分别针对差模和共模单位长度参数初值,进行C的修正,分别得到差模和共模修正后的单位长度电容C的数值,具体包括:
设实际电容
Figure FDA0004057062320000031
式中:f为频率,fend电磁仿真的最大频率值,m和n为待定系数;
根据频域中前几个周期的误差,修正m使计算得到的末端开路下始端的差模输入阻抗和测量得到的末端开路下始端的差模输入阻抗的振荡周期能吻合,如果在高频段吻合度不够理想,修正n使得整个频段能满足精度要求,通过线性修正,使得整个频段内计算得到的末端开路下始端的差模输入阻抗和测量得到的末端开路下始端的差模输入阻抗振荡周期一致,得到差模修正后的单位长度电容C的数值;
根据频域中前几个周期的误差,修正m使计算得到的末端开路下始端的共模输入阻抗和测量得到的末端开路下始端的共模输入阻抗的振荡周期能吻合,如果在高频段吻合度不够理想,修正n使得整个频段能满足精度要求,通过线性修正,使得整个频段内计算得到的末端开路下始端的共模输入阻抗和测量得到的末端开路下始端的共模输入阻抗振荡周期一致,得到共模修正后的单位长度电容C的数值。
5.根据权利要求3所述的方法,其特征在于,在得到修正后的单位长度电容C的数值后,利用线性修正方法对单位长度电导G进行求解,得到修正后的单位长度电导G的数值,利用该步骤,分别针对差模和共模单位长度参数初值,进行G的修正,分别得到差模和共模修正后的单位长度电容G的数值,具体包括:
调整tan(θ)的数值使得计算得到的末端开路下始端的差模输入阻抗和测量得到的末端开路下始端的差模输入阻抗的振荡幅度一致,得到差模修正后的单位长度电导G的数值;
调整tan(θ)的数值使得计算得到的末端开路下始端的共模输入阻抗和测量得到的末端开路下始端的共模输入阻抗的振荡幅度一致,得到共模修正后的单位长度电导G的数值。
CN201811533983.XA 2018-12-14 2018-12-14 一种三芯电力线单位长度参数矩阵求解方法 Active CN109783865B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811533983.XA CN109783865B (zh) 2018-12-14 2018-12-14 一种三芯电力线单位长度参数矩阵求解方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811533983.XA CN109783865B (zh) 2018-12-14 2018-12-14 一种三芯电力线单位长度参数矩阵求解方法

Publications (2)

Publication Number Publication Date
CN109783865A CN109783865A (zh) 2019-05-21
CN109783865B true CN109783865B (zh) 2023-04-28

Family

ID=66496981

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811533983.XA Active CN109783865B (zh) 2018-12-14 2018-12-14 一种三芯电力线单位长度参数矩阵求解方法

Country Status (1)

Country Link
CN (1) CN109783865B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103048568A (zh) * 2012-12-21 2013-04-17 广东电网公司佛山供电局 电缆沟敷设三芯电缆载流量的确定方法和确定系统
CN103502777A (zh) * 2011-01-07 2014-01-08 伍德沃德Mpc股份有限公司 用于半桥可变差动变压器位置传感系统的设备和方法
US9258829B1 (en) * 2014-09-30 2016-02-09 Texas Instruments Incorporated System and method for collision rate reduction in MIMO narrowband power line communications
CN108039869A (zh) * 2017-12-14 2018-05-15 广西师范大学 一种基于跨导系数修正结构的混频器

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9866039B2 (en) * 2015-11-13 2018-01-09 X Development Llc Wireless power delivery over medium range distances using magnetic, and common and differential mode-electric, near-field coupling

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103502777A (zh) * 2011-01-07 2014-01-08 伍德沃德Mpc股份有限公司 用于半桥可变差动变压器位置传感系统的设备和方法
CN103048568A (zh) * 2012-12-21 2013-04-17 广东电网公司佛山供电局 电缆沟敷设三芯电缆载流量的确定方法和确定系统
US9258829B1 (en) * 2014-09-30 2016-02-09 Texas Instruments Incorporated System and method for collision rate reduction in MIMO narrowband power line communications
CN108039869A (zh) * 2017-12-14 2018-05-15 广西师范大学 一种基于跨导系数修正结构的混频器

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Calculation of Multiconductor Underground Cables High-Frequency Per-Unit-Length Parameters Using Electromagnetic Modal Analysis;Shahnoor Habib et al.;《IEEE》;20130131;全文 *

Also Published As

Publication number Publication date
CN109783865A (zh) 2019-05-21

Similar Documents

Publication Publication Date Title
CN106771627B (zh) 一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法
CN108345753B (zh) 一种针对非平行线缆的串扰噪声预测方法
CN104991124B (zh) 一种特性阻抗校准系统及测试方法
Grassi et al. On mode conversion in geometrically unbalanced differential lines and its analogy with crosstalk
CN109165446B (zh) 一种长电缆入口阻抗对绝缘老化的灵敏度计算方法
Papadopoulos et al. Narrowband power line communication: Medium voltage cable modeling and laboratory experimental results
CN108763802A (zh) 计及电力传输线的电力电子换流系统电磁干扰估算方法
Razzaghi et al. On the use of electromagnetic time reversal to locate faults in series-compensated transmission lines
Zwane et al. An alternative approach in power line communication channel modelling
CN107015060B (zh) 输电导线的交流电阻测量方法和系统
CN103823127A (zh) 一种传导性共模/差模噪声分离方法
Nkom et al. Narrowband modeling of single-wire earth return distribution lines
CN102857389A (zh) 1553b总线的频域失真预测信道建模方法
CN205749844U (zh) 一种电源纹波测试装置及系统
CN109783865B (zh) 一种三芯电力线单位长度参数矩阵求解方法
CN103439580B (zh) 大尺度电气设备阻抗宽频特性的时域测量方法
CN109406877A (zh) 电容式电压互感器谐波传递系数的测量方法及装置
CN107576852A (zh) 电缆附件的特征阻抗的测量方法和系统
CN109254189B (zh) 一种海底电缆耐压试验终端电压计算方法
Bačmaga et al. High-frequency modelling of surface-mount power inductor used in switching DC-DC converters
CN105610467B (zh) 一种被复线信道建模方法
Graber et al. Validation of cable models for simulation of transients in shipboard power systems
CN109782073B (zh) 一种双芯电力线单位长度参数求解方法
Sheri et al. Characterization of a power line cable for channel frequency response—Analysis and investigation
Espín-Delgado et al. Expression for conductor resistance in the frequency range 2-150 kHz

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