CN105243197B - 基于有限元仿真和等效电路的单芯电力电缆谐振分析方法 - Google Patents
基于有限元仿真和等效电路的单芯电力电缆谐振分析方法 Download PDFInfo
- Publication number
- CN105243197B CN105243197B CN201510602421.6A CN201510602421A CN105243197B CN 105243197 B CN105243197 B CN 105243197B CN 201510602421 A CN201510602421 A CN 201510602421A CN 105243197 B CN105243197 B CN 105243197B
- Authority
- CN
- China
- Prior art keywords
- power cables
- core power
- cored wire
- metal sheath
- unit length
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 32
- 238000004088 simulation Methods 0.000 title claims abstract description 27
- 239000002184 metal Substances 0.000 claims abstract description 122
- 229910052751 metal Inorganic materials 0.000 claims abstract description 122
- 238000012544 monitoring process Methods 0.000 claims description 23
- 239000000463 material Substances 0.000 claims description 18
- 230000004888 barrier function Effects 0.000 claims description 13
- 230000005540 biological transmission Effects 0.000 claims description 11
- 230000005611 electricity Effects 0.000 claims description 9
- 238000000034 method Methods 0.000 claims description 9
- 230000001681 protective effect Effects 0.000 claims description 9
- 238000011160 research Methods 0.000 claims description 6
- 239000002689 soil Substances 0.000 claims description 6
- 230000005284 excitation Effects 0.000 claims description 5
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 5
- 239000010931 gold Substances 0.000 claims description 5
- 229910052737 gold Inorganic materials 0.000 claims description 5
- 230000011218 segmentation Effects 0.000 claims description 4
- 230000001174 ascending effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 239000004020 conductor Substances 0.000 claims description 3
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 230000005686 electrostatic field Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 2
- 238000003745 diagnosis Methods 0.000 abstract description 3
- 230000032683 aging Effects 0.000 description 3
- 239000012141 concentrate Substances 0.000 description 2
- 238000009413 insulation Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Landscapes
- Testing Of Short-Circuits, Discontinuities, Leakage, Or Incorrect Line Connections (AREA)
Abstract
本发明公开了一种基于有限元仿真和等效电路的单芯电力电缆谐振分析方法,使用有限元仿真来精确计算单芯电力电缆的芯线的单位长度电阻、芯线的单位长度电感、金属护层的单位长度电阻、金属护层的单位长度自感和芯线对金属护层的单位长度电容,并使用电路仿真来模拟分析单芯电力电缆电缆的芯线、金属护层和大地系统在不同接地方式下的谐振情况,用于电力单芯电力电缆谐振分析领域,明确了分析的对象,简化了分析的步骤,并可以给出所需的任意故障点的谐振情况,对单芯电力电缆谐振故障诊断及预测具有较大参考价值和意义。
Description
技术领域
本发明涉及电力系统设备状态监测技术领域,尤其涉及一种基于有限元仿真和等效电路的单芯电力电缆谐振分析方法。
背景技术
单芯电力电缆作为电力系统的重要组成部分,在电网中占有很重要的位置,单芯电力电缆的谐波问题在当前电力系统中日益突出。单芯电力电缆的结构使得其具有一定的分布参数,如分布电感和分布电容。在工频情况下分布参数对单芯电力电缆的正常运行影响可以忽略,但当频率较高的谐波占据了一定的成分之后,就可能在单芯电力电缆中形成局部的谐振过电压和谐振过电流,进而影响单芯电力电缆的正常工作或破坏单芯电力电缆的绝缘加速单芯电力电缆的老化。
目前对单芯电力电缆谐波危害的认识,大多数集中在谐波电流使得单芯电力电缆芯线温度提高,造成绝缘层的加速老化,对谐波的分析也集中在将电气设备断开产生的过电压、过电流作为瞬态的大幅值谐波源加载到单芯电力电缆上进行的瞬态分析。对于单芯电力电缆自身的芯线和金属护层回路系统上可能产生的高频谐波谐振问题研究较少。然而经过分析可知此类问题在一些特定的情况下可能对单芯电力电缆造成加速老化或绝缘破坏,因此对单芯电力电缆自身的芯线和金属护层回路系统上可能产生的高频谐波谐振问题的研究具有较大的应用价值。
发明内容
本发明的目的是提供一种基于有限元仿真和等效电路的单芯电力电缆谐振分析方法,用于单芯电力电缆谐振分析领域,明确了分析的对象,简化了分析的步骤,并可以给出所需的任意故障点的谐振情况,对单芯电力电缆谐振故障诊断及预测具有较大参考价值和意义。
本发明采用的技术方案为:
基于有限元仿真和等效电路的单芯电力电缆谐振分析方法,依次包括以下步骤:
A.在电力系统工频50Hz到5000Hz的谐波频率范围内,由小到大预先设定一组频率的值,计算单芯电力电缆在各频率下的分布参数,单芯电力电缆包括芯线和套设在芯线外侧的金属护层,分布参数包括芯线的单位长度电阻、芯线的单位长度电感、金属护层的单位长度电阻、金属护层的单位长度自感和芯线对金属护层的单位长度电容;
B.根据单芯电力电缆的接地方式和线路长度,在Matlab Simulink软件中做出等效电路;
C.使用Matlab Simulink的监视器功能模块,在步骤B得到的等效电路中设定监测点进行电压和电流的检测,得到仿真电路,用以观察特定单芯电力电缆位置的过电压和过电流情况;
D.根据实地采集的单芯电力电缆谐波含量数据或自行设定的谐波幅值,在步骤C得到的仿真电路上施加激励源,在每个频率下,将步骤A中计算得到的对应的分布参数代入仿真电路,并运行仿真电路,采集仿真得到芯线的电压数据、芯线的电流数据、金属护层的电压数据和金属护层的电流数据,分别做出芯线的电压或芯线的电流对频率的曲线图以及金属护层电压或金属护层电流对频率的曲线图,观察并分析是否出现谐振峰值,若出现了谐振峰值的情况,进行记录。
所述的步骤A中芯线对金属护层的单位长度电容的具体计算依次包括如下步骤:
A11.在ANSYS软件中选用二维静电场分析模块PLANE121,设定单位制为mks国际单位制,根据单芯电力电缆的几何参数,使用ANSYS软件的自带图形工具建立单芯电力电缆横截面的二维平面几何模型,其中单芯电力电缆的几何参数包括单芯电力电缆内芯线的数量以及芯线、内外半导体屏蔽层、金属护层、绝缘层、铠装层和外皮的内外半径;
A12.根据单芯电力电缆各部分的材料参数给模型赋予材料属性,材料参数包括芯线、半导体屏蔽层、绝缘层、金属屏蔽层、空气和土壤的相对介电常数及电阻率;
A13.根据单芯电力电缆的几何尺寸对整个模型进行剖分,首先将步骤A11中得到的二维平面几何模型的圆周分成至少8个相等的数段单元;其次,将这些数段单元进行控制剖分,通过设定最大单元的长度上限控制单元的总数,通常将总单元数控制在十万个内;再次,设定芯线和金属护层为两个导体,使用ANSYS自带JCG求解器对模型进行求解;最后,使用ANSYS自带计算电容的宏命令cmatrix进行电容矩阵计算,并从计算结果中读出芯线对金属护层的单位长度电容值。
所述的步骤A中芯线的单位长度电阻、芯线的单位长度电感、金属护层的单位长度电阻和金属护层的单位长度自感具体计算依次包括如下步骤:
A21.选用ANSYS有限元仿真软件中的三维时谐磁场分析模块SOLID97,根据单芯电力电缆的几何参数,使用ANSYS软件的自带图形工具建立一段长度为1m的单芯电力电缆的三维几何模型,其中单芯电力电缆的几何参数包括单芯电力电缆内芯线的数量以及芯线、内外半导体屏蔽层、金属护层、绝缘层、铠装层和外皮的内外半径;
A22.根据单芯电力电缆各部分的材料参数给三维几何模型赋予材料属性,材料参数包括芯线、半导体屏蔽层、绝缘层、金属屏蔽层、空气和土壤的相对磁导率及电阻率;
A23.根据单芯电力电缆的几何尺寸对整个三维几何模型进行剖分:将步骤A21中得到的三维几何模型的金属护层和芯线在轴向、切向和长度方向都分成均匀的数段单元,其中轴向分段时保证金属护层和芯线至少分成5段,切向分段将圆周均分为至少16段,长度方向保证单元长度小于2cm;
A24.耦合单芯电力电缆三维几何模型中芯线两端的电压,并在其中一端加载电压激励,对模型进行求解;读取单芯电力电缆芯线两端电压差矢量和芯线上电流矢量,计算出此段单芯电力电缆的芯线复阻抗,用芯线复阻抗除以此段单芯电力电缆芯线的长度得到单芯电力电缆芯线的单位长度复阻抗,单位长度复阻抗的实部即为单芯电力电缆芯线的单位长度电阻值,单位长度复阻抗的虚部即为单芯电力电缆芯线的单位长度电感值;
A25.首先,将步骤A24中的耦合和激励清零;其次,耦合单芯电力电缆三维几何模型中金属护层两端的电压,并在其中一端加载电压激励,对模型进行求解;再次,读取单芯电力电缆金属护层两端电压差矢量和金属护层上电流矢量,计算出此段单芯电力电缆的金属护层复阻抗;最后,用金属护层复阻抗除以此段单芯电力电缆金属护层的长度得到单芯电力电缆金属护层的单位长度复阻抗,单位长度复阻抗的实部即为单芯电力电缆金属护层的单位长度电阻值,单位长度复阻抗的虚部即为单芯电力电缆金属护层的单位长度电感值;
所述的步骤B依次包括以下步骤:
B1.根据传输线的分布参数等效模型,将单芯电力电缆等效为由多个单元级联组成的模型,其中每个单元由芯线的单位长度电阻、芯线的单位长度电感、芯线对金属护层的单位长度电容、金属护层的单位长度电阻和金属护层的单位长度自感组成一个二端口网络,即包含两个输入端口和两个输出端口的电路,其中芯线的单位长度电阻和芯线的单位长度电感串联连接在第一个输入端口和第一个输出端口之间,金属护层的单位长度电阻和金属护层的单位长度电感串联连接在第二个输入端口和第二个输出端口之间,芯线对金属护层的单位长度电容连接在两个输出端口之间,多个二端口网络通过级联的方式连接在一起,第二个二端口网络的两个输入端口分别接于第一个二端口网络的两个输出端口上,以此类推,第N个二端口网络两个输入端口分别接于第N-1个二端口网络的两个输出端口上;
B2.根据单芯电力电缆的接地方式或为研究目的而预设的接地方式,搭建带有接地回路、等效电源及等效负载的单芯电力电缆等效电路;等效电源由Y型联结的三相交流电压源表示,连接于单芯电力电缆电路的一端;等效负载由Y型联结的三相负载表示,连接于单芯电力电缆电路的另一端;等效电路中的地线由参考电位表示,表现为等效电路中电源与负载的连接线;其中单芯电力电缆的接地方式为单端接地、双端接地或交叉互联接地;对于单端接地,以单芯电力电缆一段接地单元的等效电路中第一个二端口网络第二个输入端为等效地电位,直接接地;将最后一个二端口网络末端输出端通过接保护电阻与地电位相连;对于双端接地,在单端接地的基础上,在第一个二端口网络的第二个输入端与地之间增加一个保护电阻;对于交叉互联接地,则根据交叉互联的规则,将三条相位各相差120度的单芯电力电缆每条分成长度相等的三段,共需要建立九段相同的单芯电力电缆一段接地单元的等效电路,并在九段接地单元等效电路之间进行连接。
所述的步骤C中监测点的设置原则为:
C1.监测点至少设置两个,分别设置在一段接地单元长度单芯电力电缆首端和末端;
C2.对于具有对称结构的电路,只在一路上设置监测点;
C3.若芯线上的电压和电流未交叉互联接地,在每一段芯线上设置监测点,金属护层上电压、电流若未交叉互联接地,在每一段芯线上设置监测点。
本发明该方法使用有限元仿真来精确计算单芯电力电缆的芯线的单位长度电阻、芯线的单位长度电感、金属护层的单位长度电阻、金属护层的单位长度自感和芯线对金属护层的单位长度电容,并使用电路仿真来模拟分析单芯电力电缆电缆的芯线、金属护层和大地系统在不同接地方式下的谐振情况,可以给出所需的任意故障点的谐振情况,对单芯电力电缆谐振故障诊断及预测具有较大参考价值和意义。
附图说明
图1为本发明的流程图;
图2为本发明实施例的分布参数等效模型;
图3为本发明实施例的单芯电力电缆等效电路;
图4为本发明实施例的单芯电力电缆仿真电路;
图5为本发明不同频率下Va,Vb,Vc的仿真结果;
图6为本发明不同频率下va,vb,vc的仿真结果。
具体实施方式
如图1所示,本发明所述的基于有限元仿真和等效电路的单芯电力电缆谐振分析方法,依次包括以下步骤:
A.在电力系统工频50Hz到5000Hz的谐波频率范围内,由小到大预先设定一组频率的值,计算单芯电力电缆在各频率下的分布参数,单芯电力电缆包括芯线和套设在芯线外侧的金属护层,分布参数包括芯线的单位长度电阻、芯线的单位长度电感、金属护层的单位长度电阻、金属护层的单位长度自感和芯线对金属护层的单位长度电容;在本实施例中,从工频50Hz到50000Hz选取9组频率,选取的频率范围如表1所示。
在传输线理论中,传输线长度l与工作波长λ相差不大的传输线称为“长线”,传输线长度l远小于工作波长λ的传输线称为“短线”。单芯电力电缆传输的50Hz电力属于低频电磁波,其波长约为6000km,远远大于单芯电力电缆的铺设尺寸,即单芯电力电缆可看作“短线”,所以当仅考虑工频50Hz电力时,单芯电力电缆可以用集总参数模型来表示,其模型可以用电阻元件和电感元件的串联表示。但由于在实际单芯电力电缆工作状态中,输电设备和负载的非线性导致单芯电力电缆中存在一定量的高频谐波。谐波频率越高,波长就越短,当谐波的波长接近单芯电力电缆的铺设尺寸时,单芯电力电缆近似满足“长线”条件,此时集总参数模型不再能准确描述单芯电力电缆的工作状态,而需要分布参数模型来表示。分布参数模型可以视为无穷多个电路单元的级联,然而无穷单元在电路仿真中不易实现,可以采取合理近似,将无穷多单元降为有限个单元。本发明以110kV型号为YJLW031101X100010kV的单芯电力电缆为研究对象,考虑到电力系统中通常只考察100次以下的谐波,即波长6km以上的谐波,而单芯电力电缆通常每隔数百米至一两公里将通过单芯电力电缆接头进行接地,可以近似地将两接地点间的每段单芯电力电缆看作三个单元的分布式参数。
表1单芯电力的电缆仿真频率
表2不同频率下单芯电力的电缆的单位长度电容、单位长度电感、单位长度电阻值
由于在电力系统工频50Hz到5000Hz的谐波频率范围内,单芯电力电缆的芯线对金属护层的单位长度电容变化很小,因此可以采用静态电容计算。芯线对金属护层的单位长度电容的具体计算依次包括如下步骤:
A11.在ANSYS软件中选用二维静电场分析模块PLANE121,设定单位制为mks国际单位制,以保证计算出来的参数值为单位长度的参数值,根据单芯电力电缆的几何参数,使用ANSYS软件的自带图形工具建立单芯电力电缆横截面的二维平面几何模型,其中单芯电力电缆的几何参数包括单芯电力电缆内芯线的数量以及芯线、内外半导体屏蔽层、金属护层、绝缘层、铠装层和外皮的内外半径;
A12.根据单芯电力电缆各部分的材料参数给模型赋予材料属性,材料参数包括芯线、半导体屏蔽层、绝缘层、金属屏蔽层、空气和土壤的相对介电常数及电阻率;
A13.根据单芯电力电缆的几何尺寸对整个模型进行剖分,首先将步骤A11中得到的二维平面几何模型的圆周分成至少8个相等的数段单元;其次,将这些数段单元进行控制剖分,通过设定最大单元的长度上限控制单元的总数,划分单元数量越多,计算越精确,但由于计算机计算能力的限制,通常将总单元数控制在十万之内;再次,设定芯线和金属护层为两个导体,使用ANSYS自带JCG求解器对模型进行求解;最后,使用ANSYS自带计算电容的宏命令cmatrix进行电容矩阵计算,并从计算结果中读出芯线对金属护层的单位长度电容值。
由于电感随频率的变化而变化,所以不能采用如计算电容的静态方法来计算电感。此处采用VA法,即使用单芯电力电缆电压除以单芯电力电缆电流得到复数形式的阻抗值,阻抗值的实部即为电阻,虚部即为电感。芯线的单位长度电阻、芯线的单位长度电感、金属护层的单位长度电阻和金属护层的单位长度自感具体计算依次包括如下步骤:
A21.选用ANSYS有限元仿真软件中的三维时谐磁场分析模块SOLID97,根据单芯电力电缆的几何参数,使用ANSYS软件的自带图形工具建立一段长度为1m的单芯电力电缆的三维几何模型,其中单芯电力电缆的几何参数包括单芯电力电缆内芯线的数量以及芯线、内外半导体屏蔽层、金属护层、绝缘层、铠装层和外皮的内外半径;
A22.根据单芯电力电缆各部分的材料参数给模型赋予材料属性,材料参数包括芯线、半导体屏蔽层、绝缘层、金属屏蔽层、空气和土壤的相对磁导率及电阻率;
A23.根据单芯电力电缆的几何尺寸对整个模型进行剖分:将步骤A21中得到的三维几何模型的金属护层和芯线在轴向、切向和长度方向都分成均匀的数段单元,其中轴向分段时保证金属护层和芯线至少分成5段,切向分段将圆周均分为至少16段,长度方向保证单元长度小于2cm,以确保模型中不出现影响计算精度的长度方向边长远大于其他方向边长的细长单元;
A24.耦合单芯电力电缆三维几何模型中芯线两端的电压,并在其中一端加载电压激励,对三维几何模型进行求解;读取单芯电力电缆芯线两端电压差矢量和芯线上电流矢量,计算出此段单芯电力电缆的芯线复阻抗,用芯线复阻抗除以此段单芯电力电缆芯线的长度得到单芯电力电缆芯线的单位长度复阻抗,单位长度复阻抗的实部即为单芯电力电缆芯线的单位长度电阻值,单位长度复阻抗的虚部即为单芯电力电缆芯线的单位长度电感值;
A25.首先,将步骤A24中的耦合和激励清零;其次,耦合单芯电力电缆模型中金属护层两端的电压,并在其中一端加载电压激励,对模型进行求解;再次,读取单芯电力电缆金属护层两端电压差矢量和金属护层上电流矢量,计算出此段单芯电力电缆的金属护层复阻抗;最后,用金属护层复阻抗除以此段单芯电力电缆金属护层的长度得到单芯电力电缆金属护层的单位长度复阻抗,单位长度复阻抗的实部即为单芯电力电缆金属护层的单位长度电阻值,单位长度复阻抗的虚部即为单芯电力电缆金属护层的单位长度电感值;
B.根据YJLW031101X100010kV的单芯电力电缆的接地方式和线路长度,在MatlabSimulink软件中做出等效电路。
B1.根据传输线的分布参数等效模型,将单芯电力电缆等效为由多个单元级联组成的模型,其中每个单元由芯线的单位长度电阻、芯线的单位长度电感、芯线对金属护层的单位长度电容、金属护层的单位长度电阻和金属护层的单位长度自感组成一个二端口网络,即包含两个输入端口和两个输出端口的电路,其中芯线的单位长度电阻和芯线的单位长度电感串联,连接于第一个输入端口和第一个输出端口之间;金属护层的单位长度电阻和金属护层的单位长度电感相互串联,连接于第二个输入端口和第二个输出端口之间;芯线对金属护层的单位长度电容连接在两个输出端口之间。多个二端口网络通过级联的方式连接在一起,第二个二端口网络的两个输入端口分别接在第一个二端口网络的两个输出端口上,以此类推,第N个二端口网络两个输入端口分别接于第N-1个二端口网络的两个输出端口上;
级联的单元的数量选取,理论上越多越接近真实情况,但在实际操作中,只要小于十分之一工作波长即可满足要求。如图2所示,本实施例中设定所分析单芯电力电缆的接地单元长度为1km,取3级进行计算足够满足上述要求,即一段接地单元长度的单芯电力电缆等效电路由3个级联单元组成。其中RC,LC表示单芯电力电缆芯线上的单位长度电阻和单位长度电感,RS,LS表示单芯电力电缆金属护层上的单位长度电阻和单位长度电感,C表示单芯电力电缆芯线和金属护层之间的单位长度电容。dx是该段单芯电力电缆长度的1/3。电容C上有一定大小的串联阻抗。此处等效模型考虑了单芯电力电缆的传输线特性,以及金属护层上的电阻及电感作用。
B2.根据单芯电力电缆的接地方式或为研究目的而预设的接地方式,搭建带有接地回路、等效电源及等效负载的单芯电力电缆等效电路;
如图3所示,等效电源由Y型联结的三相交流电压源1表示,连接于单芯电力电缆电路的一端;等效负载由Y型联结的三相负载2表示,连接于单芯电力电缆电路的另一端;等效电路中的地线由参考电位表示,表现为等效电路中电源与负载的连接线;其中单芯电力电缆的接地方式为单端接地、双端接地或交叉互联接地;对于单端接地,以单芯电力电缆一段接地单元的等效电路中第一个二端口网络第二个输入端为等效地电位,直接接地;将最后一个二端口网络末端输出端通过接保护电阻3与地电位相连;对于双端接地,在单端接地的基础上,在第一个二端口网络的第二个输入端与地之间增加一个保护电阻3;对于交叉互联接地,则根据交叉互联的规则,将三条相位各相差120度的单芯电力电缆每条分成长度相等的三段,共需要建立九段相同的单芯电力电缆一段接地单元的等效电路,并在九段接地单元等效电路之间进行连接。具体连接方式如下:为了方便叙述,将九段接地单元进行编号,处于第一段的三条相位不同的小段为1-1段,1-2段和1-3段,处于第二段的三条相位不同的小段为2-1段,2-2段和2-3段,处于第三段的三条相位不同的小段为3-1段,3-2段和3-3段。根据交叉互联规则,将1-1段、2-1段、3-1段的芯线首尾相连,将1-2段、2-2段、3-2段的芯线首尾相连,将1-3段、2-3段、3-3段的芯线首尾相连;再将1-1段、2-2段、3-3段的金属护层首尾相连,将1-2段、2-3段、3-1段的金属护层首尾相连,将1-3段、2-1段、3-2段的金属护层首尾相连;根据交叉互连规则,将1-1段、1-2段、1-3段的金属护层末端分别串联一个保护电阻3,此例中设为50兆欧,并接入地线;同理将2-1段、2-2段、2-3段的金属护层末端分别串联一个保护电阻3,并接入地线。本例中设定交叉互联接地方式。保护电阻3的大小根据实际测量值或预设值设置,本例中预设为50兆欧。
对整体的电路设置电源与负载。电源与负载的设置可根据所研究的实际单芯电力电缆的谐波工况情况进行调整。单芯电力电缆的谐波工况包含基波电压幅值及相位、基波电流幅值及相位、各次谐波电压幅值及相位、各次谐波电流幅值及相位。本例中采用Y型联接的三相理想电源进行电源模拟,采用Y型联接的三相平衡纯阻性负载进行负载模拟。本例中为了使结果与实际工况接近,设置三相理想电源的工频电压为110kV,三相负载电阻的大小均为220欧,三相电流由设置的三相工频电压及三相负载决定,根据仿真结果工频下其有效值约为500A。
C.使用Matlab Simulink的监视器功能模块,在步骤B得到的等效电路中设定监测点,得到仿真电路,用以观察特定单芯电力电缆位置的过电压和过电流情况,所述监测点能检测电压和电流,监测点的设置原则为:
C1.监测点至少设置两个,分别设置在一段接地单元长度单芯电力电缆首端和末端;
C2.对于具有对称结构的电路,只在一路上设置监测点;
C3.若芯线上的电压和电流未交叉互联接地,在每一段芯线上设置监测点,金属护层上电压、电流若未交叉互联接地,在每一段芯线上设置监测点。
本实施例中,设置监测点如图4所示,图中较大圆圈表示电压表,较小圆圈表示电流表。若较大圆圈画于电容符号之上,表示在该电容两端并联一个电压表;若较大圆圈一端连接某线路节点,另一端与表示接地的符号相连,表示在该电路节点与地之间连接一个电压表;若较小圆圈画于一段线路之上,表示在该线路上串联一个电流表。其中,Ia,Ib,Ic分别为其中一相上分别处于三大段电流流入侧的芯线电流,ia,ib,ic为三段上的金属护层流入端的电流,且ia,ib,ic分别在三个不同相的芯线外,其共同点是所处的金属护层处在一条连接通路上;Va,Vb,Vc分别为其中一相上分别处于三大段的芯线电压,且都处于每一大段的第一个单元上;va,vb,vc为三段上的金属护层对地电压,且都处于每一大段的第一个单元上。根据这些监测点可以在不同频率下和不同的单芯电力电缆条件下,观察单芯电力电缆各处可能出现的谐振情况。针对三相不平衡的情况,设置了监测点Ia',Ia"与Ia相对应,ia',ia"与ia相对应,Va',Va"与Va相对应,va',va"与va相对应。其中符号上的撇号代表同一小段距离上相位相差120度的另两条电缆上相同位置的电压或电流。
D:根据实地采集的单芯电力电缆谐波含量数据或自行设定的谐波幅值,在步骤C得到的仿真电路上施加激励源,在每个频率下,将步骤A中计算得到的对应的分布参数代入仿真电路,并运行仿真电路,采集仿真得到芯线的电压数据、芯线的电流数据、金属护层的电压数据和金属护层的电流数据,分别做出芯线的电压或芯线的电流对频率的曲线图以及金属护层电压或金属护层电流对频率的曲线图,观察并分析是否出现谐振峰值,若出现了谐振峰值的情况,进行记录。本实施例中,仿真结果如图5和图6所示。
Claims (5)
1.一种基于有限元仿真和等效电路的单芯电力电缆谐振分析方法,其特征在于,依次包括以下步骤:
A.在电力系统工频50Hz到5000Hz的谐波频率范围内,由小到大预先设定一组频率的值,计算单芯电力电缆在各频率下的分布参数,单芯电力电缆包括芯线和套设在芯线外侧的金属护层,分布参数包括芯线的单位长度电阻、芯线的单位长度电感、金属护层的单位长度电阻、金属护层的单位长度自感和芯线对金属护层的单位长度电容;
B.根据单芯电力电缆的接地方式和线路长度,在MatlabSimulink软件中做出等效电路;
C.使用MatlabSimulink的监视器功能模块,在步骤B得到的等效电路中设定监测点进行电压和电流的检测,得到仿真电路,用以观察特定单芯电力电缆位置的过电压和过电流情况;
D.根据实地采集的单芯电力电缆谐波含量数据或自行设定的谐波幅值,在步骤C得到的仿真电路上施加激励源,在每个频率下,将步骤A中计算得到的对应的分布参数代入仿真电路,并运行仿真电路,采集仿真得到芯线的电压数据、芯线的电流数据、金属护层的电压数据和金属护层的电流数据,分别做出芯线的电压或芯线的电流对频率的曲线图以及金属护层电压或金属护层电流对频率的曲线图,观察并分析是否出现谐振峰值,若出现了谐振峰值的情况,进行记录。
2.根据权利要求1所述的一种基于有限元仿真和等效电路的单芯电力电缆谐振分析方法,其特征在于,所述的步骤A中芯线对金属护层的单位长度电容的具体计算依次包括如下步骤:
A11.在ANSYS软件中选用二维静电场分析模块PLANE121,设定单位制为mks国际单位制,根据单芯电力电缆的几何参数,使用ANSYS软件的自带图形工具建立单芯电力电缆横截面的二维平面几何模型,其中单芯电力电缆的几何参数包括单芯电力电缆内芯线的数量以及芯线、内外半导体屏蔽层、金属护层、绝缘层、铠装层和外皮的内外半径;
A12.根据单芯电力电缆各部分的材料参数给模型赋予材料属性,材料参数包括芯线、半导体屏蔽层、绝缘层、金属屏蔽层、空气和土壤的相对介电常数及电阻率;
A13.根据单芯电力电缆的几何尺寸对整个模型进行剖分,首先将步骤A11中得到的二维平面几何模型的圆周分成至少8个相等的数段单元;其次,将这些数段单元进行控制剖分,通过设定最大单元的长度上限控制单元的总数,通常将总单元数控制在十万个内;再次,设定芯线和金属护层为两个导体,使用ANSYS自带JCG求解器对模型进行求解;最后,使用ANSYS自带计算电容的宏命令cmatrix进行电容矩阵计算,并从计算结果中读出芯线对金属护层的单位长度电容值。
3.根据权利要求1所述的一种基于有限元仿真和等效电路的单芯电力电缆谐振分析方法,其特征在于:所述的步骤A中芯线的单位长度电阻、芯线的单位长度电感、金属护层的单位长度电阻和金属护层的单位长度自感具体计算依次包括如下步骤:
A21.选用ANSYS有限元仿真软件中的三维时谐磁场分析模块SOLID97,根据单芯电力电缆的几何参数,使用ANSYS软件的自带图形工具建立一段长度为1m的单芯电力电缆的三维几何模型,其中单芯电力电缆的几何参数包括单芯电力电缆内芯线的数量以及芯线、内外半导体屏蔽层、金属护层、绝缘层、铠装层和外皮的内外半径;
A22.根据单芯电力电缆各部分的材料参数给三维几何模型赋予材料属性,材料参数包括芯线、半导体屏蔽层、绝缘层、金属屏蔽层、空气和土壤的相对磁导率及电阻率;
A23.根据单芯电力电缆的几何尺寸对整个三维几何模型进行剖分:将步骤A21中得到的三维几何模型的金属护层和芯线在轴向、切向和长度方向都分成均匀的数段单元,其中轴向分段时保证金属护层和芯线至少分成5段,切向分段将圆周均分为至少16段,长度方向保证单元长度小于2cm;
A24.耦合单芯电力电缆三维几何模型中芯线两端的电压,并在其中一端加载电压激励,对模型进行求解;读取单芯电力电缆芯线两端电压差矢量和芯线上电流矢量,计算出此段单芯电力电缆的芯线复阻抗,用芯线复阻抗除以此段单芯电力电缆芯线的长度得到单芯电力电缆芯线的单位长度复阻抗,单位长度复阻抗的实部即为单芯电力电缆芯线的单位长度电阻值,单位长度复阻抗的虚部即为单芯电力电缆芯线的单位长度电感值;
A25.首先,将步骤A24中的耦合和激励清零;其次,耦合单芯电力电缆三维几何模型中金属护层两端的电压,并在其中一端加载电压激励,对模型进行求解;再次,读取单芯电力电缆金属护层两端电压差矢量和金属护层上电流矢量,计算出此段单芯电力电缆的金属护层复阻抗;最后,用金属护层复阻抗除以此段单芯电力电缆金属护层的长度得到单芯电力电缆金属护层的单位长度复阻抗,单位长度复阻抗的实部即为单芯电力电缆金属护层的单位长度电阻值,单位长度复阻抗的虚部即为单芯电力电缆金属护层的单位长度电感值。
4.根据权利要求1所述的一种基于有限元仿真和等效电路的单芯电力电缆谐振分析方法,其特征在于:所述的步骤B依次包括以下步骤:
B1.根据传输线的分布参数等效模型,将单芯电力电缆等效为由多个单元级联组成的模型,其中每个单元由芯线的单位长度电阻、芯线的单位长度电感、芯线对金属护层的单位长度电容、金属护层的单位长度电阻和金属护层的单位长度自感组成一个二端口网络,即包含两个输入端口和两个输出端口的电路,其中芯线的单位长度电阻和芯线的单位长度电感串联连接在第一个输入端口和第一个输出端口之间,金属护层的单位长度电阻和金属护层的单位长度电感串联连接在第二个输入端口和第二个输出端口之间,芯线对金属护层的单位长度电容连接在两个输出端口之间,多个二端口网络通过级联的方式连接在一起,第二个二端口网络的两个输入端口分别接于第一个二端口网络的两个输出端口上,以此类推,第N个二端口网络两个输入端口分别接于第N-1个二端口网络的两个输出端口上;
B2.根据单芯电力电缆的接地方式或为研究目的而预设的接地方式,搭建带有接地回路、等效电源及等效负载的单芯电力电缆等效电路;等效电源由Y型联结的三相交流电压源表示,连接于单芯电力电缆电路的一端;等效负载由Y型联结的三相负载表示,连接于单芯电力电缆电路的另一端;等效电路中的地线由参考电位表示,表现为等效电路中电源与负载的连接线;其中单芯电力电缆的接地方式为单端接地、双端接地或交叉互联接地;对于单端接地,以单芯电力电缆一段接地单元的等效电路中第一个二端口网络第二个输入端为等效地电位,直接接地;将最后一个二端口网络末端输出端通过接保护电阻与地电位相连;对于双端接地,在单端接地的基础上,在第一个二端口网络的第二个输入端与地之间增加一个保护电阻;对于交叉互联接地,则根据交叉互联的规则,将三条相位各相差120度的单芯电力电缆每条分成长度相等的三段,共需要建立九段相同的单芯电力电缆一段接地单元的等效电路,并在九段接地单元等效电路之间进行连接。
5.根据权利要求1所述的一种基于有限元仿真和等效电路的单芯电力电缆谐振分析方法,其特征在于,所述的步骤C中监测点的设置原则为:
C1.监测点至少设置两个,分别设置在一段接地单元长度单芯电力电缆首端和末端;
C2.对于具有对称结构的电路,只在一路上设置监测点;
C3.若芯线上的电压和电流未交叉互联接地,在每一段芯线上设置监测点,金属护层上电压、电流若未交叉互联接地,在每一段芯线上设置监测点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510602421.6A CN105243197B (zh) | 2015-09-21 | 2015-09-21 | 基于有限元仿真和等效电路的单芯电力电缆谐振分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510602421.6A CN105243197B (zh) | 2015-09-21 | 2015-09-21 | 基于有限元仿真和等效电路的单芯电力电缆谐振分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105243197A CN105243197A (zh) | 2016-01-13 |
CN105243197B true CN105243197B (zh) | 2018-02-09 |
Family
ID=55040845
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510602421.6A Active CN105243197B (zh) | 2015-09-21 | 2015-09-21 | 基于有限元仿真和等效电路的单芯电力电缆谐振分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105243197B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107918711A (zh) * | 2017-11-23 | 2018-04-17 | 广东电网有限责任公司珠海供电局 | 一种双回电缆输电线路的电缆最优布置方式筛选方法 |
CN107992713B (zh) * | 2018-01-24 | 2021-02-19 | 武汉大学 | 一种组合空气间隙击穿电压预测方法 |
CN108761167B (zh) * | 2018-03-01 | 2019-12-10 | 华南理工大学 | 一种电缆金属护套多相多点接地下护层感应电流计算方法 |
CN108509712B (zh) * | 2018-03-29 | 2022-07-22 | 南方电网科学研究院有限责任公司 | 一种电缆线路的过电压仿真分析方法及系统 |
CN109522581B (zh) * | 2018-07-06 | 2023-04-18 | 西安电子科技大学 | 一种加载集总元件的反射超表面单元特性分析方法 |
CN109241647B (zh) * | 2018-09-21 | 2022-09-09 | 国网陕西省电力公司电力科学研究院 | 一种基于场路结合获取电力元件等效模型和参数的方法 |
CN110135008B (zh) * | 2019-04-19 | 2024-03-26 | 中铁二院工程集团有限责任公司 | 一种单芯电缆护层电流的计算方法 |
CN110161375B (zh) * | 2019-06-21 | 2021-07-09 | 国网山东省电力公司电力科学研究院 | 一种基于分布电阻参数的高压直流输电线路计算模型 |
CN117388761B (zh) * | 2023-12-11 | 2024-03-19 | 国网山东省电力公司淄博供电公司 | 一种基于变频谐振的多芯二次电缆核线装置及方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155583A (zh) * | 2013-09-26 | 2014-11-19 | 扬州市交大工业技术研究院有限公司 | 一种高压单芯电力电缆局部放电在线监测方法及系统 |
CN103698653B (zh) * | 2013-12-30 | 2016-05-04 | 国家电网公司 | 基于环流测量的单芯电力电缆护套接地故障检测系统及方法 |
CN104899379B (zh) * | 2015-06-10 | 2018-07-10 | 国网河南省电力公司电力科学研究院 | 基于ansys的交联聚乙烯绝缘电缆温度场有限元计算方法 |
-
2015
- 2015-09-21 CN CN201510602421.6A patent/CN105243197B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105243197A (zh) | 2016-01-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105243197B (zh) | 基于有限元仿真和等效电路的单芯电力电缆谐振分析方法 | |
Zhang et al. | Diagnosis of breaks in substation's grounding grid by using the electromagnetic method | |
CN106771627B (zh) | 一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法 | |
Patel et al. | MoM-SO: A complete method for computing the impedance of cable systems including skin, proximity, and ground return effects | |
Patel et al. | Accurate impedance calculation for underground and submarine power cables using MoM-SO and a multilayer ground model | |
Papadopoulos et al. | Narrowband power line communication: Medium voltage cable modeling and laboratory experimental results | |
CN111123041B (zh) | 一种基于温度特性的电缆护层故障定位方法 | |
Pagnetti et al. | An improved method for the calculation of the internal impedances of solid and hollow conductors with the inclusion of proximity effect | |
Candela et al. | A model for assessing the magnitude and distribution of sheath currents in medium and high-voltage cable lines | |
CN103344837B (zh) | 一种同杆并架双回直流输电线路的参数测量方法 | |
Gudmundsdottir | Proximity effect in fast transient simulations of an underground transmission cable | |
Du et al. | The extended thin-wire model of lossy round wire structures for FDTD simulations | |
CN105005665A (zh) | 油纸电容式套管末屏的特高频辐射特性仿真方法与系统 | |
CN106802361A (zh) | 一种配网混连接地运行的电容电流测量方法及系统 | |
Kruizinga et al. | High frequency modeling of a shielded four-core low voltage underground power cable | |
CN105510733A (zh) | 一种高压、特高压直流输电线路的参数测量方法 | |
Araneo et al. | Hybrid and pi-circuit approaches for grounding system lightning response | |
Hosseini et al. | Transformer winding modeling based on multi-conductor transmission line model for partial discharge study | |
CN103745054A (zh) | 一种对电缆及电缆束在电磁兼容中的建模及信号串扰分析方法 | |
CN109254189B (zh) | 一种海底电缆耐压试验终端电压计算方法 | |
del-Pino-López et al. | A 3D parametric analysis of three-core armored power cables series impedance | |
Casagrande et al. | High frequency bundles modeling | |
Jurisic et al. | Application of wide band transformer models | |
CN110472366A (zh) | 一种载波信号在中压配电网地埋电缆线路传输的解耦方法 | |
Parmantier et al. | Simplification method for the assessment of the EM response of a complex cable harness |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |