CN107831370A - 直流接地极及杆塔接地网附近电磁场的有限元分析方法 - Google Patents
直流接地极及杆塔接地网附近电磁场的有限元分析方法 Download PDFInfo
- Publication number
- CN107831370A CN107831370A CN201711236145.1A CN201711236145A CN107831370A CN 107831370 A CN107831370 A CN 107831370A CN 201711236145 A CN201711236145 A CN 201711236145A CN 107831370 A CN107831370 A CN 107831370A
- Authority
- CN
- China
- Prior art keywords
- grounding
- tower
- tower grounding
- direct current
- grounding electrode
- 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.)
- Granted
Links
- 230000005672 electromagnetic field Effects 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title abstract description 5
- 239000002689 soil Substances 0.000 claims abstract description 61
- 238000004458 analytical method Methods 0.000 claims abstract description 25
- 239000000463 material Substances 0.000 claims abstract description 10
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 18
- 230000005540 biological transmission Effects 0.000 claims description 18
- 229910052799 carbon Inorganic materials 0.000 claims description 18
- 238000010586 diagram Methods 0.000 claims description 16
- 230000035699 permeability Effects 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 13
- 230000007797 corrosion Effects 0.000 claims description 10
- 238000005260 corrosion Methods 0.000 claims description 10
- 238000009933 burial Methods 0.000 claims description 5
- 230000010261 cell growth Effects 0.000 claims description 5
- 229910001209 Low-carbon steel Inorganic materials 0.000 claims description 4
- 229910000831 Steel Inorganic materials 0.000 claims description 4
- 239000010959 steel Substances 0.000 claims description 4
- 239000000571 coke Substances 0.000 claims description 3
- 230000012010 growth Effects 0.000 claims description 3
- 230000035772 mutation Effects 0.000 abstract 1
- 239000010410 layer Substances 0.000 description 25
- 238000011160 research Methods 0.000 description 7
- 241000287828 Gallus gallus Species 0.000 description 4
- 230000005684 electric field Effects 0.000 description 4
- 208000031481 Pathologic Constriction Diseases 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000007772 electrode material Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 230000036262 stenosis Effects 0.000 description 1
- 208000037804 stenosis Diseases 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/0864—Measuring electromagnetic field characteristics characterised by constructional or functional features
- G01R29/0892—Details related to signal analysis or treatment; presenting results, e.g. displays; measuring specific signal features other than field strength, e.g. polarisation, field modes, phase, envelope, maximum value
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明公开了一种直流接地极及杆塔接地网附近电磁场的有限元分析方法，通过建立模型，设置相关参数、添加材料属性、添加边界条件、网格化分、电磁场分布计算，得出结论：两个接地极环上的电位最高，内环电压略高于外环，以接地极圆环的圆心为起点沿径向电位逐渐降低，内环土壤的电位较外环土壤的电位高；杆塔接地网本体上的电位最高，越靠近杆塔接地网电位越高，杆塔接地网内部土壤的电位比外部土壤的电位高，因为杆塔接地网有对称性，所以它附近的电位也呈现出对称性；接地网的射线末端泄露电流密度最大,射线首段的电流密度最小，接地网矩形和射线的连接处电流密度有突变，本发明解决了现有技术中存在的直流接地极及杆塔接地网附近电磁场的有限元分析不全面的问题。
Description
技术领域
本发明属于电力系统及其自动化技术领域，具体涉及一种直流接地极及杆塔接地网附近电磁场的有限元分析方法。
背景技术
直流输电系统处于单极大地回线运行方式时，会有很大的直流电流通过直流接地极流入大地，这将造成接地极本身及附近输电杆塔的接地网腐蚀。因此研究直流输电接地极周围电场分布规律，以及泄露电流流入杆塔接地网后其周围的电场和泄露电流密度的分布情况，对掌握直流接地极及杆塔接地网的腐蚀规律具有重要意义。
现有的文献对于直流接地极及附近杆塔接地网的电磁场计算有一些研究，文献中研究了多层土壤模型下各层土壤的电阻率对直流接地极附近地表电位分布情况的影响，但是没有考虑接地极的其它电气参数对电位分布情况的影响。有的文献用数值分析法计算了高压直流输电直线型接地极地表电位的分布情况，但对于双圆环型接地极未作讨论。有的文献分别采用了单层土壤模型和双层土壤模型下接地极周围地表电位的分布规律，但是实际上大地土壤的结构更为复杂。有的文献对杆塔接地网的泄露电流密度做了研究，但是未对杆塔接地网附近的地表电位进行研究。
本文在理论分析接地极地表电位的计算方法的基础上，采用有限元分析软件COMSOL Multiphysics，以德宝直流输电工程千阳接地极为例，建立多层大地土壤结构下的直流接地极和杆塔接地网数值模型，添加相关边界条件，进行网格划分处理，分析计算高压直流输电系统处于单极大地运行方式时，接地极和周围杆塔接地网的地表电位及杆塔接地网本体上的泄露电流密度。
发明内容
本发明的目的是提供一种直流接地极及杆塔接地网附近电磁场的有限元分析方法，解决了现有技术中存在的直流接地极及杆塔接地网附近电磁场的有限元分析不全面的问题。
本发明所采用的技术方案是，直流接地极及杆塔接地网附近电磁场的有限元分析方法，具体按照以下步骤实施：
步骤1、建立模型，设置相关参数；
步骤2、添加材料属性；
步骤3、添加边界条件；
步骤4、网格化分；
步骤5、电磁场分布计算。
本发明的特点还在于，
步骤1具体为：
针对直流接地极：建立直流接地极模型，设圆环接地极内环和外环的半径R1和R2，单个极环的截面直径D1，接地极的埋深为h1，接地极极环周围覆盖一层碳层，外环碳层尺寸为m1*m1，内环碳层尺寸为m2*m2；
设土壤模型为四层土壤模型，层深分别为H1＝0-0.03km，H2＝0.03-10km，H3＝10-50km，H4＝50-∞km，土壤横向半径R3＝50km；
针对杆塔接地网，建立杆塔接地网模型：设杆塔根开L，矩形边长a，基础深为h2，射线长度b，接地网埋深为h3，采用的材料横截面直径为D2，设杆塔接地网附近的土壤模型为边长a＝150m的正方体。
步骤2具体为：
针对直流接地极：极环采用低碳钢棒Q235，设极环的电阻率为ρ1，相对磁导率为μ1，相对介电常数为ε1，碳层的材料为焦炭，设碳层电阻率为ρ2，相对磁导率为μ2，相对介电常数为ε2，四层土壤模型中各层土壤的电阻率分别为ρ31，ρ32，ρ32，ρ34，相对磁导率分别为μ31，μ32，μ33，μ34，相对介电常数分别为ε31，ε32，ε33，ε34；
针对杆塔接地网：杆塔接地网的材料为D12的圆钢，设电阻率为ρ4，相对磁导率为μ4，相对介电常数为ε4，土壤的电阻率为ρ5，相对磁导率为μ5，相对介电常数为ε5。
步骤3具体为：
针对直流接地极：向接地极环注入接地极最大持续额定入地电流I1，I1为直流输电工程处于单极大地回线运行方式时最大持续额定入地电流，设定无穷远处的电压为0V；
针对杆塔接地网：向杆塔接地网注入的泄漏电流为I2,将所建立的杆塔接地网附近土壤模型的下边界设置为接地。
步骤4具体为：
针对直流接地极：网格划分采用自由剖分三角形网格划分，网格中单元的尺寸需设置如下参数:最大单元尺寸n1，最小单元尺寸n2，最大单元增长率λ，曲率因子η，狭窄区域的分辨率δ；
针对杆塔接地网：网格划分采用自由剖分四面体网格划分，网格中单元的尺寸需设置如下参数:最大单元尺寸n1，最小单元尺寸n2，最大单元增长率λ，曲率因子η，狭窄区域的分辨率δ。
步骤5具体为：
利用如下公式将前述关于直流接地极电磁场计算相关参数带入，
其中，V为接地极附近空间某点的电位，ρ为土壤的电阻率，I1为当直流输电系统处于单极大地回路运行方式时流入接地极的电流，K(δ)以δ为变量的第一类椭圆积分，R1为接地极环内环的半径，R2为接地极环外环的半径，h为接地极环的埋深，r为接地极地表沿径向距离，α为空间某点与接地接环径向的夹角，
计算直流接地极附近电磁场有如下结论：
针对直流接地极：设置计算接地极地表电位分布情况，以及接地极地表沿径向距离电位分布，得出结论：两个接地极环上的电位最高，内环电压略高于外环，以接地极圆环的圆心为起点沿径向电位逐渐降低，内环土壤的电位较外环土壤的电位高；
利用如下公式将前述关于杆塔接地网电磁场计算相关参数带入，
因为b>>D，所以
其中，J为泄露电流密度，I2为注入杆塔接地网的泄漏电流，V为杆塔接地网附近空间某点的电位，b为杆塔接地网射线的长度，ρ为土壤的电阻率，l为沿杆塔接地网射线方向的距离，D为杆塔接地网所用材料的横截面直径，h2为杆塔接地网的埋深；
计算杆塔接地网附近电磁场有如下结论：
针对杆塔接地网：设置计算杆塔接地网周围电位分布：三维分布和横切图，以及杆塔接地网沿射线方向的泄露电流密度分布，得出结论：杆塔接地网本体上的电位最高，越靠近杆塔接地网电位越高，杆塔接地网内部土壤的电位比外部土壤的电位高，因为杆塔接地网有对称性，所以它附近的电位也呈现出对称性，接地网的射线末端泄露电流密度最大，射线首段的电流密度最小，接地网矩形和射线的连接处电流密度有突变，所以杆塔接地网射线末端腐蚀最严重，矩形与射线的连接处的腐蚀情况次之。
本发明的有益效果是，直流接地极及杆塔接地网附近电磁场的有限元分析方法，通过建立多层大地土壤结构下的直流接地极和杆塔接地网数值模型，添加相关边界条件，进行网格划分处理，计算分析了接地极地表电位分布规律，并对杆塔接地网附近电位及泄漏电流密度进行了研究，结果发现:接地极地表电位沿径向距离逐渐降低；杆塔接地网本体上的电位最高，接地网的射线末端泄露电流密度最大,射线首端的泄露电流密度最小，接地网矩形与射线的连接处电流密度有突变。本文的研究对掌握直流接地极及杆塔接地网周围电场分布情况和腐蚀规律，具有重要的意义。
附图说明
图1是本发明直流接地极及杆塔接地网附近电磁场的有限元分析方法中宝鸡千阳接地极示意图；
图2是本发明直流接地极及杆塔接地网附近电磁场的有限元分析方法中千阳接地极数值模型；
图3是本发明直流接地极及杆塔接地网附近电磁场的有限元分析方法中接地极地表电位分布图；
图4是本发明直流接地极及杆塔接地网附近电磁场的有限元分析方法中接地极地表沿径向距离电位分布图；
图5是本发明直流接地极及杆塔接地网附近电磁场的有限元分析方法中杆塔接地网布置图；
图6是本发明直流接地极及杆塔接地网附近电磁场的有限元分析方法中杆塔接地网数值模型图；
图7是本发明直流接地极及杆塔接地网附近电磁场的有限元分析方法中杆塔接地网周围电位三维分布图；
图8是本发明直流接地极及杆塔接地网附近电磁场的有限元分析方法中杆塔接地网周围电位分布横切图；
图9是本发明直流接地极及杆塔接地网附近电磁场的有限元分析方法中杆塔接地网沿射线方向的泄露电流密度分布图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明直流接地极及杆塔接地网附近电磁场的有限元分析方法，具体按照以下步骤实施：
步骤1、建立模型，设置相关参数，具体为：
针对直流接地极：建立直流接地极模型，设圆环接地极内环和外环的半径R1和R2，单个极环的截面直径D1，接地极的埋深为h1，接地极极环周围覆盖一层碳层，外环碳层尺寸为m1*m1，内环碳层尺寸为m2*m2；
设土壤模型为四层土壤模型，层深分别为H1＝0-0.03km，H2＝0.03-10km，H3＝10-50km，H4＝50-∞km，土壤横向半径R3＝50km；
针对杆塔接地网，建立杆塔接地网模型：设杆塔根开L，矩形边长a，基础深为h2，射线长度b，接地网埋深为h3，采用的材料横截面直径为D2，设杆塔接地网附近的土壤模型为边长a＝150m的正方体；
步骤2、添加材料属性，具体为：
针对直流接地极：极环采用低碳钢棒Q235，设极环的电阻率为ρ1，相对磁导率为μ1，相对介电常数为ε1，碳层的材料为焦炭，设碳层电阻率为ρ2，相对磁导率为μ2，相对介电常数为ε2，四层土壤模型中各层土壤的电阻率分别为ρ31，ρ32，ρ33，ρ34，相对磁导率分别为μ31，μ32，μ33，μ34，相对介电常数分别为ε31，ε32，ε33，ε34；
针对杆塔接地网：杆塔接地网的材料为D12的圆钢，设电阻率为ρ4，相对磁导率为μ4，相对介电常数为ε4，土壤的电阻率为ρ5，相对磁导率为μ5，相对介电常数为ε5；
步骤3、添加边界条件，具体为：
针对直流接地极：向接地极环注入接地极最大持续额定入地电流I1，I1为直流输电工程处于单极大地回线运行方式时最大持续额定入地电流，设定无穷远处的电压为0V；
针对杆塔接地网：向杆塔接地网注入的泄漏电流为I2,将所建立的杆塔接地网附近土壤模型的下边界设置为接地；
步骤4、网格化分，具体为：
针对直流接地极：网格划分采用自由剖分三角形网格划分，网格中单元的尺寸需设置如下参数:最大单元尺寸n1，最小单元尺寸n2，最大单元增长率λ，曲率因子η，狭窄区域的分辨率δ；
针对杆塔接地网：网格划分采用自由剖分四面体网格划分，网格中单元的尺寸需设置如下参数:最大单元尺寸n1，最小单元尺寸n2，最大单元增长率λ，曲率因子η，狭窄区域的分辨率δ；
步骤5、电磁场分布计算，具体为：
利用如下公式将前述关于直流接地极电磁场计算相关参数带入，
其中，V为接地极附近空间某点的电位，ρ为土壤的电阻率，I1为当直流输电系统处于单极大地回路运行方式时流入接地极的电流，K(δ)以δ为变量的第一类椭圆积分，R1为接地极环内环的半径，R2为接地极环外环的半径，h为接地极环的埋深，r为接地极地表沿径向距离，α为空间某点与接地接环径向的夹角，
计算直流接地极附近电磁场有如下结论：
针对直流接地极：设置计算接地极地表电位分布情况，以及接地极地表沿径向距离电位分布，得出结论：两个接地极环上的电位最高，内环电压略高于外环，以接地极圆环的圆心为起点沿径向电位逐渐降低，内环土壤的电位较外环土壤的电位高；
利用如下公式将前述关于杆塔接地网电磁场计算相关参数带入，
因为b>>D，所以
其中，J为泄露电流密度，I2为注入杆塔接地网的泄漏电流，V为杆塔接地网附近空间某点的电位，b为杆塔接地网射线的长度，ρ为土壤的电阻率，l为沿杆塔接地网射线方向的距离，D为杆塔接地网所用材料的横截面直径，h2为杆塔接地网的埋深；
计算杆塔接地网附近电磁场有如下结论：
针对杆塔接地网：设置计算杆塔接地网周围电位分布：三维分布和横切图，以及杆塔接地网沿射线方向的泄露电流密度分布，得出结论：杆塔接地网本体上的电位最高，越靠近杆塔接地网电位越高，杆塔接地网内部土壤的电位比外部土壤的电位高，因为杆塔接地网有对称性，所以它附近的电位也呈现出对称性，接地网的射线末端泄露电流密度最大，射线首段的电流密度最小，接地网矩形和射线的连接处电流密度有突变，所以杆塔接地网射线末端腐蚀最严重，矩形与射线的连接处的腐蚀情况次之。
本发明直流接地极及杆塔接地网附近电磁场的有限元分析方法，以德宝直流输电工程千阳接地极为例，查阅相关文献，宝鸡千阳接地极极址附近的土壤分层参数如表1所示：
表1宝鸡千阳接地极极址实际土壤分层参数
| 层号 | 层深/km | 电阻率/Ω·m |
| 1 | 0.03 | 53.7 |
| 2 | 9.97 | 5000 |
| 3 | 50 | 100000 |
| 4 | ∞ | 100 |
德宝直流输电系统宝鸡千阳直流接地极采用浅埋型陆地同心双圆环型接地极，外环半径内环半径为电极埋深3.5m。电极材料为低碳钢棒Q235，直径电阻率为1.7×10-7Ω·m。接地极线路的直流电阻值为0.31Ω，接地极接地电阻设计值为0.039Ω。接地极极环周围覆盖一层碳层，外环碳层尺寸为0.7m×0.7m，内环碳层尺寸为0.6m×0.6m。图1为宝鸡千阳接地极示意图。
当±500kV德宝直流输电系统处于单极大地回线运行方式时，单极额定传输功率为1500MW，额定电流为3000A，接地极最大持续额定入地电流为3000A。
根据上述土壤分层参数和接地极几何参数及电气参数,建立直流接地极数值模型，建模时使用如下方法和假设：
(1)由于千阳接地极为同心双圆环型的结构，接地极呈现二维轴对称特性。因此将三维模型转换为二维轴对称模型，这样能够提高仿真计算的速度。
(2)相对于土壤电阻率和接地体电阻率的差异，可以将接地极简化为等电位体来考虑。
(3)边界条件假设离接地极径向50km处为零电位。
在COMSOL Multiphysics建立接地极的数值模型如图2所示。当直流输电系统以单极大地回线方式运行时，设置如下边界条件：向接地极环注入3000A的直流入地电流，无穷远处电位为零即边界条件为接地。网格划分采用软件内置的自由剖分三角形网格划分，网格中单元的尺寸采用内置的超细化：最大单元尺寸为10000m，最小单元尺寸为37.5m，最大单元增长率为1.2，曲率因子为0.25，狭窄区域的分辨率为1。计算结果如图3和图4所示：
图3为接地极地表电位分布图，图中颜色越深，表明其电位越高。从图3中可以得出：接地极环上的电位最高，内环部分土壤的电位较高，外环外的土壤电位较低。总体电势是以接地极圆环的圆心为起点沿半径方向逐渐降低。
图4为宝鸡千阳接地极地表沿径向距离电位分布，图中给出了0km到50km的地表电位分布图，从图4中可以看出：接地极内环电位为134.6V,外环电位为134.15V，从0km到10km地表电位有一个很大的降落，约为120V，从10km到50km地表电位降落缓慢。
直流输电工程附近的杆塔接地网的腐蚀也是一个很严重的问题，考察杆塔基础泄漏电流对杆塔接地网附近电场分布的影响，并进一步探究直流接地极附近杆塔接地网的腐蚀规律具有重要意义。
下面以常见的杆塔接地网模型为例，图5为杆塔接地网实际工程中的布置图，杆塔根开均为7m，矩形边长均为12m，基础深为2.5m，接地体材料为D12的圆钢，射线长度为52m，接地网埋深为0.8m。
根据上述杆塔接地网的参数，在COMSOL Multiphysics中建立数值模型如图6所示：
为了研究杆塔接地网附近电位以及本体上泄漏电流密度的分布情况，此处泄漏电流取1A来进行研究。
在COMSOL Multiphysics建立杆塔接地网几何模型后，设置如下边界条件：向杆塔接地网注入1A的泄漏电流,再设置接地边界条件。网格划分采用自由剖分四面体网格划分，网格中单元的尺寸采用内置的极细化：最大单元尺寸为3m，最小单元尺寸为0.03m，最大单元增长率为0.3，曲率因子为0.2，狭窄区域的分辨率为1.计算结果如图7、图8所示，图中颜色越深，代表电位越高。
图7为杆塔接地网周围电位三维分布图，从图7可以看出，越靠近杆塔接地网,电位越高，杆塔接地网本体上的电位最高，可以达到0.51V。图8为杆塔接地网周围电位分布横切图，可以发现杆塔接地网内部土壤的电位比外部土壤的电位高。因为杆塔接地网有对称性，所以它附近的电位也呈现对称性。
图9为杆塔接地网沿射线方向的泄露电流密度分布图，由图9可以看出，接地网的射线末端泄露电流密度最大约为0.161609A/㎡,射线首端的电流密度约为0.005A/㎡，图中电流密度在两处有突变，这两处为接地网矩形与射线的连接部分，它们的泄露电流密度约为0.02A/㎡，其他部位电流密度变化不大。泄露电流密度越大，腐蚀越严重，因此接地网矩形与射线的连接部分腐蚀略为严重，射线末端的腐蚀情况最严重。
Claims (6)
1.直流接地极及杆塔接地网附近电磁场的有限元分析方法，其特征在于，具体按照以下步骤实施：
步骤1、建立模型，设置相关参数；
步骤2、添加材料属性；
步骤3、添加边界条件；
步骤4、网格化分；
步骤5、电磁场分布计算。
2.根据权利要求1所述的直流接地极及杆塔接地网附近电磁场的有限元分析方法，其特征在于，所述步骤1具体为：
针对直流接地极：建立直流接地极模型，设圆环接地极内环和外环的半径R1和R2，单个极环的截面直径D1，接地极的埋深为h1，接地极极环周围覆盖一层碳层，外环碳层尺寸为m1*m1，内环碳层尺寸为m2*m2；
设土壤模型为四层土壤模型，层深分别为H1＝0-0.03km，H2＝0.03-10km，H3＝10-50km，H4＝50-∞km，土壤横向半径R3＝50km；
针对杆塔接地网，建立杆塔接地网模型：设杆塔根开L，矩形边长a，基础深为h2，射线长度b，接地网埋深为h3，采用的材料横截面直径为D2，设杆塔接地网附近的土壤模型为边长a＝150m的正方体。
3.根据权利要求1所述的直流接地极及杆塔接地网附近电磁场的有限元分析方法，其特征在于，所述步骤2具体为：
针对直流接地极：极环采用低碳钢棒Q235，设极环的电阻率为ρ1，相对磁导率为μ1，相对介电常数为ε1，碳层的材料为焦炭，设碳层电阻率为ρ2，相对磁导率为μ2，相对介电常数为ε2，四层土壤模型中各层土壤的电阻率分别为ρ31，ρ32，ρ33，ρ34，相对磁导率分别为μ31，μ32，μ33，μ34，相对介电常数分别为ε31，ε32，ε33，ε34；
针对杆塔接地网：杆塔接地网的材料为D12的圆钢，设电阻率为ρ4，相对磁导率为μ4，相对介电常数为ε4，土壤的电阻率为ρ5，相对磁导率为μ5，相对介电常数为ε5。
4.根据权利要求1所述的直流接地极及杆塔接地网附近电磁场的有限元分析方法，其特征在于，所述步骤3具体为：
针对直流接地极：向接地极环注入接地极最大持续额定入地电流I1，此时，直流输电工程处于单极大地回线运行方式时，设定无穷远处的电压为0V；
针对杆塔接地网：向杆塔接地网注入的泄漏电流为I2,将所建立的杆塔接地网附近土壤模型的下边界设置为接地。
5.根据权利要求1所述的直流接地极及杆塔接地网附近电磁场的有限元分析方法，其特征在于，所述步骤4具体为：
针对直流接地极：网格划分采用自由剖分三角形网格划分，网格中单元的尺寸需设置如下参数:最大单元尺寸n1，最小单元尺寸n2，最大单元增长率λ，曲率因子η，狭窄区域的分辨率δ；
针对杆塔接地网：网格划分采用自由剖分四面体网格划分，网格中单元的尺寸需设置如下参数:最大单元尺寸n1，最小单元尺寸n2，最大单元增长率λ，曲率因子η，狭窄区域的分辨率δ。
6.根据权利要求1所述的直流接地极及杆塔接地网附近电磁场的有限元分析方法，其特征在于，所述步骤5具体为：
利用如下公式将前述关于直流接地极电磁场计算相关参数带入，
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;I</mi>
<mn>1</mn>
</msub>
</mrow>
<mi>&pi;</mi>
</mfrac>
<mo>*</mo>
<mfrac>
<mrow>
<mi>K</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&delta;</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mi>r</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;I</mi>
<mn>1</mn>
</msub>
</mrow>
<mi>&pi;</mi>
</mfrac>
<mo>*</mo>
<mfrac>
<mrow>
<mi>K</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&delta;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mi>r</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&delta;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mn>4</mn>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mi>r</mi>
</mrow>
<mrow>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
</mrow>
<mrow>
<msub>
<mi>&delta;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mn>4</mn>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mi>r</mi>
</mrow>
<mrow>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
</mrow>
<mrow>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mi>&delta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mfrac>
<mi>&pi;</mi>
<mn>2</mn>
</mfrac>
</msubsup>
<mfrac>
<mrow>
<mi>d</mi>
<mi>&tau;</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>&delta;</mi>
<mn>2</mn>
</msup>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
</mfrac>
</mrow>
<mrow>
<mi>&tau;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&pi;</mi>
<mo>-</mo>
<mi>&alpha;</mi>
</mrow>
<mn>2</mn>
</mfrac>
</mrow>
其中，V为接地极附近空间某点的电位，ρ为土壤的电阻率，I1为当直流输电系统处于单极大地回路运行方式时流入接地极的电流，K(δ)以δ为变量的第一类椭圆积分，R1为接地极环内环的半径，R2为接地极环外环的半径，h为接地极环的埋深，r为接地极地表沿径向距离，α为空间某点与接地接环径向的夹角，
计算直流接地极附近电磁场有如下结论：
针对直流接地极：设置计算接地极地表电位分布情况，以及接地极地表沿径向距离电位分布，得出结论：两个接地极环上的电位最高，内环电压略高于外环，以接地极圆环的圆心为起点沿径向电位逐渐降低，内环土壤的电位较外环土壤的电位高；
利用如下公式将前述关于杆塔接地网电磁场计算相关参数带入，
<mrow>
<mi>J</mi>
<mo>=</mo>
<mfrac>
<mi>I</mi>
<mi>s</mi>
</mfrac>
</mrow>
因为b>>D，所以
<mrow>
<mi>J</mi>
<mo>=</mo>
<mfrac>
<mi>I</mi>
<mi>b</mi>
</mfrac>
</mrow>
<mrow>
<mi>V</mi>
<mo>=</mo>
<mo>&Integral;</mo>
<mfrac>
<mrow>
<mi>&rho;</mi>
<mi>d</mi>
<mi>l</mi>
</mrow>
<mrow>
<msub>
<mi>&pi;Dh</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<msub>
<mi>I</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<mi>R</mi>
<mo>=</mo>
<mfrac>
<mi>&rho;</mi>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>d</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mi>ln</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>d</mi>
</mrow>
<msqrt>
<mrow>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<mi>D</mi>
</mrow>
</msqrt>
</mfrac>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>I</mi>
<mo>=</mo>
<mfrac>
<mi>V</mi>
<mi>R</mi>
</mfrac>
</mrow>
其中，J为泄露电流密度，I2为注入杆塔接地网的泄漏电流，V为杆塔接地网附近空间某点的电位，b为杆塔接地网射线的长度，ρ为土壤的电阻率，l为沿杆塔接地网射线方向的距离，D为杆塔接地网所用材料的横截面直径，h2为杆塔接地网的埋深；
计算杆塔接地网附近电磁场有如下结论：
针对杆塔接地网：设置计算杆塔接地网周围电位分布：三维分布和横切图，以及杆塔接地网沿射线方向的泄露电流密度分布，得出结论：杆塔接地网本体上的电位最高，越靠近杆塔接地网电位越高，杆塔接地网内部土壤的电位比外部土壤的电位高，因为杆塔接地网有对称性，所以它附近的电位也呈现出对称性，接地网的射线末端泄露电流密度最大，射线首段的电流密度最小，接地网矩形和射线的连接处电流密度有突变，所以杆塔接地网射线末端腐蚀最严重，矩形与射线的连接处的腐蚀情况次之。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711236145.1A CN107831370B (zh) | 2017-11-30 | 2017-11-30 | 直流接地极及杆塔接地网附近电磁场的有限元分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711236145.1A CN107831370B (zh) | 2017-11-30 | 2017-11-30 | 直流接地极及杆塔接地网附近电磁场的有限元分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107831370A true CN107831370A (zh) | 2018-03-23 |
| CN107831370B CN107831370B (zh) | 2020-11-17 |
Family
ID=61646736
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711236145.1A Expired - Fee Related CN107831370B (zh) | 2017-11-30 | 2017-11-30 | 直流接地极及杆塔接地网附近电磁场的有限元分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107831370B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109684780A (zh) * | 2019-02-19 | 2019-04-26 | 国网内蒙古东部电力有限公司检修分公司 | 一种特高压直流输电双环接地极最优内外环比计算方法 |
| CN110108948A (zh) * | 2019-04-19 | 2019-08-09 | 贵州电网有限责任公司 | 一种带补偿修正的接地电阻短距测量方法及系统 |
| CN111967161A (zh) * | 2020-08-17 | 2020-11-20 | 重庆大学 | 一种用于塔杆腐蚀接地极接地参数分析的数值计算方法 |
| CN112668212A (zh) * | 2020-09-02 | 2021-04-16 | 国网内蒙古东部电力有限公司检修分公司 | 基于有限元的不同土壤模型下接地极溢流特性分析方法 |
| CN113325156A (zh) * | 2020-10-14 | 2021-08-31 | 河南四达电力设备股份有限公司 | 一种基于有限元法的接地材料土壤腐蚀评价及寿命预测方法 |
| CN114295871A (zh) * | 2021-12-30 | 2022-04-08 | 国网湖南省电力有限公司 | 输电高压线路杆塔接地网不开挖检测方法及系统 |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003186919A (ja) * | 2001-12-19 | 2003-07-04 | Japan Research Institute Ltd | メッシュ生成方法、メッシュ生成装置、コンピュータプログラム、及び記録媒体 |
| JP2005018403A (ja) * | 2003-06-26 | 2005-01-20 | Matsushita Electric Ind Co Ltd | 電磁界の解析方法、電磁界の解析プログラムおよびこの解析プログラムを記録した記録媒体 |
| JP2007213169A (ja) * | 2006-02-07 | 2007-08-23 | Japan Research Institute Ltd | 有限要素解析装置、有限要素解析方法、及びコンピュータプログラム |
| US20120053910A1 (en) * | 2010-12-07 | 2012-03-01 | General Electric Company | Method and system for inductor power loss analysis |
| CN102710210A (zh) * | 2012-06-19 | 2012-10-03 | 合肥工业大学 | 一种凸极式永磁同步电机驱动系统的控制方法 |
| CN103593523A (zh) * | 2013-11-12 | 2014-02-19 | 国网上海市电力公司 | 基于有限元理论的多直流落点条件下直流偏磁抑制方法 |
| CN103646153A (zh) * | 2013-12-25 | 2014-03-19 | 哈尔滨工业大学 | 用于描绘金属颗粒在交变磁场中运动轨迹的方法 |
| CN104698313A (zh) * | 2015-03-03 | 2015-06-10 | 国网四川省电力公司电力科学研究院 | 多直流接地极不同运行方式下直流偏磁电流影响站点的预测方法 |
| CN105631096A (zh) * | 2015-12-21 | 2016-06-01 | 云南电网有限责任公司电力科学研究院 | 一种110kV/220kV输电线路杆塔接地体射线长度估算方法 |
| CN106844830A (zh) * | 2016-12-08 | 2017-06-13 | 宁波大学 | 一种快速预测耦合玻璃通孔互连传输特性的数值方法 |
| CN107784187A (zh) * | 2017-11-22 | 2018-03-09 | 重庆大学 | 一种基于边界元法的变电站电磁环境分布研究方法 |
-
2017
- 2017-11-30 CN CN201711236145.1A patent/CN107831370B/zh not_active Expired - Fee Related
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003186919A (ja) * | 2001-12-19 | 2003-07-04 | Japan Research Institute Ltd | メッシュ生成方法、メッシュ生成装置、コンピュータプログラム、及び記録媒体 |
| JP2005018403A (ja) * | 2003-06-26 | 2005-01-20 | Matsushita Electric Ind Co Ltd | 電磁界の解析方法、電磁界の解析プログラムおよびこの解析プログラムを記録した記録媒体 |
| JP2007213169A (ja) * | 2006-02-07 | 2007-08-23 | Japan Research Institute Ltd | 有限要素解析装置、有限要素解析方法、及びコンピュータプログラム |
| US20120053910A1 (en) * | 2010-12-07 | 2012-03-01 | General Electric Company | Method and system for inductor power loss analysis |
| CN102710210A (zh) * | 2012-06-19 | 2012-10-03 | 合肥工业大学 | 一种凸极式永磁同步电机驱动系统的控制方法 |
| CN103593523A (zh) * | 2013-11-12 | 2014-02-19 | 国网上海市电力公司 | 基于有限元理论的多直流落点条件下直流偏磁抑制方法 |
| CN103646153A (zh) * | 2013-12-25 | 2014-03-19 | 哈尔滨工业大学 | 用于描绘金属颗粒在交变磁场中运动轨迹的方法 |
| CN104698313A (zh) * | 2015-03-03 | 2015-06-10 | 国网四川省电力公司电力科学研究院 | 多直流接地极不同运行方式下直流偏磁电流影响站点的预测方法 |
| CN105631096A (zh) * | 2015-12-21 | 2016-06-01 | 云南电网有限责任公司电力科学研究院 | 一种110kV/220kV输电线路杆塔接地体射线长度估算方法 |
| CN106844830A (zh) * | 2016-12-08 | 2017-06-13 | 宁波大学 | 一种快速预测耦合玻璃通孔互连传输特性的数值方法 |
| CN107784187A (zh) * | 2017-11-22 | 2018-03-09 | 重庆大学 | 一种基于边界元法的变电站电磁环境分布研究方法 |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109684780A (zh) * | 2019-02-19 | 2019-04-26 | 国网内蒙古东部电力有限公司检修分公司 | 一种特高压直流输电双环接地极最优内外环比计算方法 |
| CN109684780B (zh) * | 2019-02-19 | 2023-06-30 | 国网内蒙古东部电力有限公司检修分公司 | 一种特高压直流输电双环接地极最优内外环比计算方法 |
| CN110108948A (zh) * | 2019-04-19 | 2019-08-09 | 贵州电网有限责任公司 | 一种带补偿修正的接地电阻短距测量方法及系统 |
| CN111967161A (zh) * | 2020-08-17 | 2020-11-20 | 重庆大学 | 一种用于塔杆腐蚀接地极接地参数分析的数值计算方法 |
| CN111967161B (zh) * | 2020-08-17 | 2023-10-17 | 重庆大学 | 一种用于塔杆腐蚀接地极接地参数分析的数值计算方法 |
| CN112668212A (zh) * | 2020-09-02 | 2021-04-16 | 国网内蒙古东部电力有限公司检修分公司 | 基于有限元的不同土壤模型下接地极溢流特性分析方法 |
| CN113325156A (zh) * | 2020-10-14 | 2021-08-31 | 河南四达电力设备股份有限公司 | 一种基于有限元法的接地材料土壤腐蚀评价及寿命预测方法 |
| CN114295871A (zh) * | 2021-12-30 | 2022-04-08 | 国网湖南省电力有限公司 | 输电高压线路杆塔接地网不开挖检测方法及系统 |
| CN114295871B (zh) * | 2021-12-30 | 2023-07-14 | 国网湖南省电力有限公司 | 输电高压线路杆塔接地网不开挖检测方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107831370B (zh) | 2020-11-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107831370B (zh) | 直流接地极及杆塔接地网附近电磁场的有限元分析方法 | |
| CN104901030B (zh) | 一种变电站降阻接地装置施工方法 | |
| CN103267926A (zh) | 基于差分进化算法故障特征匹配的含dg配电网故障测距 | |
| CN110119544B (zh) | 一种适用于复杂环境区域的螺旋接地极尺寸参数设计方法 | |
| Cong et al. | Root-cause identification of single line-to-ground fault in urban small current grounding systems based on correlation dimension and average resistance | |
| CN104392055A (zh) | 组合式复合材料杆塔防雷优化设计方法 | |
| CN108197393B (zh) | 一种具有降阻功能的螺旋形接地极设计方法 | |
| CN113946923A (zh) | 电网对埋地管网的电磁干扰计算方法 | |
| CN109975596A (zh) | 一种单极大地运行方式下地中电流分布研究方法 | |
| CN107831353B (zh) | 一种新型基于圆环形直流输电接地极不同埋深的跨步电压的计算方法 | |
| Zeng et al. | Analysis on influence of long vertical grounding electrodes on grounding system for substation | |
| CN109241664A (zh) | 一种关于风力发电机雷电电磁暂态特性的分析计算方法 | |
| CN106814284A (zh) | 一种接地网网内电势差问题计算方法 | |
| Zalhaf et al. | Computation of transient induced voltages along a wind turbine struck by lightning | |
| CN108398722A (zh) | 一种基于磁场微分法的变电站接地网拓扑结构绘制方法 | |
| CN109815593B (zh) | 一种配网相继故障耦合传播机理的分析方法 | |
| He et al. | Research on grounding current distribution in HVDC system operating in monopolar mode | |
| Xiong et al. | Improvement of shaped conductive backfill material for grounding systems | |
| Manikandan et al. | Design and analysis of grounding systems for wind turbines using Finite element method | |
| CN110795861A (zh) | 一种对换流站支柱绝缘子屏蔽球地刀口参数的优化方法 | |
| CN115422790A (zh) | 一种山坡地形中的电参数计算方法 | |
| CN108804763A (zh) | 一种基于分流测试结果的垂直型直流接地极跨步电势计算方法 | |
| Jiang et al. | Surface Electric Field Simulation of Worker in Live Working on±800kV UHVDC Transmission Line | |
| CN114357367B (zh) | 适用于直流偏磁评估的变电站互电阻快速傅里叶变换方法 | |
| CN114709782A (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201117 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |