CN105427190A - 复杂地形时uhvac输电线路下方地面三维工频电场计算方法 - Google Patents
复杂地形时uhvac输电线路下方地面三维工频电场计算方法 Download PDFInfo
- Publication number
- CN105427190A CN105427190A CN201510909603.8A CN201510909603A CN105427190A CN 105427190 A CN105427190 A CN 105427190A CN 201510909603 A CN201510909603 A CN 201510909603A CN 105427190 A CN105427190 A CN 105427190A
- Authority
- CN
- China
- Prior art keywords
- ground
- discrete
- charges
- electric field
- point
- 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
- 230000005684 electric field Effects 0.000 title claims abstract description 73
- 238000000034 method Methods 0.000 title claims abstract description 57
- 230000005540 biological transmission Effects 0.000 title claims abstract description 31
- 239000004020 conductor Substances 0.000 claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims abstract description 16
- 238000013178 mathematical model Methods 0.000 claims abstract description 7
- 230000010354 integration Effects 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 46
- 238000005457 optimization Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 5
- 238000011161 development Methods 0.000 claims description 4
- GNFTZDOKVXKIBK-UHFFFAOYSA-N 3-(2-methoxyethoxy)benzohydrazide Chemical compound COCCOC1=CC=CC(C(=O)NN)=C1 GNFTZDOKVXKIBK-UHFFFAOYSA-N 0.000 claims description 3
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000002939 conjugate gradient method Methods 0.000 claims description 3
- 230000008707 rearrangement Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 238000000819 phase cycle Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开一种复杂地形时UHVAC输电线路下方地面三维工频电场计算方法。首先将复杂地形简化成起伏形式的三维地面,建立计算模型;其次选定镜像面,并假定镜像地面延伸至起伏地面内;接着将带电导线离散成分段的线电荷,在非镜像地面内以正六边形蜂窝状形式设置离散点电荷,同时依据镜像原理设置相应的镜像电荷并建立三维电场积分数学模型的离散方程;然后选定基函数为分域脉冲函数,按照配点法选权函数并在每个积分域求内积,形成矩阵方程;接着,求解矩阵方程,利用所求电荷计算匹配点最大误差,直至满足要求;最后利用最优的离散电荷计算起伏地面的三维工频电场分布。本发明计算方法降低了计算难度,提高了计算效率,能够更好地满足实际工程需要。
Description
技术领域
本发明涉及一种电场计算方法,尤其涉及一种针对复杂地形时UHVAC(UltraHighVoltageAlternatingCurrent:特高压交流)输电线路下方地面三维工频电场计算方法。
背景技术
随着我国社会经济的快速发展,用电需求不断增加。同时,越来越大的环境压力使得在符合中心地区新建各类电源变得难以执行,为应对这种情况,必须要将远方的电能输送到负荷中心地带。特高压交流输电具有优势长距离输送、减少能量损失,便于联网等优点,使得特高压交流输电在我国获得了快速发展。
然而,特高压交流导线路在运行过程中会产生电磁环境影响问题,其主要的影响因素之一就是工频电场。我国幅员辽阔,地形复杂,为准确评估特高压架空交流输电线路工频电场对周围环境的影响,必须要进行三维计算。目前国内外有关工频电场的三维计算主要采用3维模拟电荷法、有限元法、矩量法等。由于3维模拟电荷法,计算快速,且容易掌握,在实际工程中该方法获得了广泛应用。但是目前,三维工频电场计算主要是针对主要针对地面是理想的平地而进行的,忽略了复杂地形对工频电场的影响。但是,在特高压导线路实际架设过程中,没有绝对的理想水平地面,特高压交流线路会经常跨越丘陵、山谷、平原等不同的地形,相关的现场试验也表明起伏地面对工频电场分布存在不容忽视的影响。目前,针对特高压架空导线路工频电场的计算均以理想的平地为参考,造成导线路的工频电场计算造成较大的误差,使得计算结果缺乏可参考性,虽然少数研究学者进行了起伏地形的三维电场计算,但是由于设置的模拟电荷量太大,使得实用性大大降低。
因此,需要提出一种新的对于复杂地形情况的特高压架空交流输电线路下方地面三维工频电场计算方法,能够对复杂地形下发的特架空交流输电线路三维工频电场进行有效的计算,保证计算精度的同时,还能够大大提高计算效率,而且适用性广。
发明内容
为解决现有技术中存在的问题,本发明的目的是提供的一种复杂地形时UHVAC输电线路下方地面三维工频电场计算方法,其能够对复杂地形情况下的三维工频电场进行有效的计算,保证计算精度的同时,还能够大大提高计算效率,而且适用性广。
本发明的目的是通过以下的技术方案实现的,它包括有以下内容:
一种复杂地形时UHVAC输电线路下方地面三维工频电场计算方法包括以
下步骤:
a.建立三维电场计算坐标,将复杂地形简化成起伏形式的三维地面;
b.选定镜像地面,并假定镜像地面延伸至起伏地面内;
c.将带电导线离散成分段的线电荷,在非镜像地面内以正六边形蜂窝状形式设置离散点电荷,同时依据镜像原理设置相应的镜像电荷并建立三维电场积分数学模型的离散方程;
d.然后选定基函数为分域脉冲函数,按照配点法选权函数并在每个积分域求内积,形成矩阵方程;
e.求解矩阵方程,计算电荷分布;
f.计算匹配点最大电位误差,如果误差不满足要求,返回步骤c重新对场源离散;如果满足要求进行下一步;
g.利用优化方法对离散电荷开展进一步优化配置。
h.利用最优电荷分布计算起伏地面的三维工频电场分布。
优选的技术方案,步骤a中将复杂地形简化成起伏形式的三维地面,是根据具体的地形进行简化,简化的原则是将复杂地形采用具有一定角度的多个平面组合来代表。
优选的技术方案,步骤b中选定的镜像地面为导线下方的地面,该地面为一段水平地面。
优选的技术方案,步骤c中将带电导线离散成分段的线电荷,是考虑导线的弧垂,将导线分成足够小的线段,以满足将小的线段看成直线电荷单元。
优选的技术方案,步骤c中在非镜像地面内部上以正六边蜂窝状形式设置离散的点电荷与地表面垂直设置,同时正六边蜂窝状的边长根据电场计算精度进行改变。
优选的技术方案,步骤c中建立的三维电场积分数学模型离散方程的场源包括导线电荷和非镜像地面内部电荷,同时方程还考虑了地形的影响,其函数关系如下:
其中为场点矢量,导线线源矢量,导线镜像线源矢量,非镜像地面内部点源矢量,镜像点源矢量,l为线电荷所存在区域,lj为镜像线电荷所在区域,q为离散点电荷。
优选的技术方案,步骤d中,选定的基函数为分域脉冲函数,其函数式如下:
优选的技术方案,步骤d中,按照配点法选择的权函数如下:
其中表示离散源点至计算点的矢径,表示空间任意一点至计算点的矢径,当时,ωj=∞,当当时,ωj=0。
优选的技术方案,步骤d中,矩阵方程的形成通过如下步骤实现:
首先,将式(2)代入式(1),得出电位积分方程的离散形式:
其次,在积分区域内对每个ωj与式(3)求内积,则有:
最后,根据δj函数性质及式(2)可以将(4)简化成下式:
令
则有:
优选的技术方案,所述矩阵方程的系数Pji通过如下方法实现:若以P1Pi为起点,沿PPiPPj方向建立局部坐标u,并设线单元长度为L0,设线电荷密度τ在单元内按线性规律分布则τ(u)=au+b(u的范围为(0,L0)),令l=xj-xi,m=yj-yi,n=zj-zi,Ei=l2+m2+n2,Fij=-2(l(xj-xi)+m(yj-yi)+n(zj-zi)),Gij=(xj-xi)2+(yj-yi)2+(zj-zi)2,对开展积分计算有:
优选的技术方案,步骤f中,匹配点最大电位误差应当小于5%,如果大于5%,则重新布置离散电荷通过如下方法实现:根据校验电位误差对离散电荷的位置及个数进行设置:重新设置过程中,对误差较大的区域,以离散点电荷为正六边形的中心进行分裂,增加离散点电荷的数量来提高计算精度。
优选的技术方案,步骤g中利用优化方法对离散电荷开展进一步优化配置,首先建立目标函数,本发明的目标函数通过求解场域边界上所有匹配点的已知电位和计算电位的差值平方之和的最小值来表示:
其中,为所有离散电荷在第i个匹配点处的电位;为第i个匹配点的已知电位,导线表面φ=U,非镜像地面上φ=0;式(6)具有如下约束条件:
S1.离散电荷的电量为自由变量;S2.离散电荷的位置必须在无效计算场域内:
zQ-f(zQ)<0i=m+1…n(8)
其中,在式(7)和式(8)中,xQd,yQd,zQd为导线上电荷坐标,xo,yo,zo为子导线中心坐标,r为子导线半径,zQ为起伏地面地下离散电荷纵坐标,m为子导线数,并且采用共轭梯度法求解式(6)中的极小值,从而求得最优的离散电荷。
优选的技术方案,步骤h中根据最优离散电荷计算起伏地面三维工频电场分布通过如下方式实现:
设复杂地面某点为P(x,y,z),则该点的电场强度可由离散线电荷和离散点电荷产生的电场叠加获得;XYZ三个方向的电场强度由下式求出:
其中,i为离散线电荷个数,j为离散点电荷个数,L0为离散线电荷长度,A=aL0,B=b。x,y,z为待求点坐标,xi,yi,zi为线源坐标,xj,yj,zj为点源坐标。
则P点的电场强度有效值为:
由于采用了上述技术方案,本发明具有如下的优点:
本发明方法基于镜像原理并结合矩量法能对交流架空输电线路所经过的复杂地形下方的三维工频电场分布进行准确计算,该方法较比能有效地克服了传统方法中对于复杂地形的三维电场分布计算量大的问题,增强了电场计算方法的适用性。本发明可以为交流架空输电线路跨越复杂地形时的规划设计、环境评价等发挥作用。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书和权利要求书来实现和获得。
附图说明
图1为本发明的流程图。
图2为本发明中的地形简化示意图。
图3为本发明中的镜像面选择示意图。
图4为本发明中的离散电荷设置和分裂。
图5为本发明中的三维模拟线电荷单元。
图6为本发明中的三维工频电场计算图。
图7为本发明中的电位误差分布图。
图8为本发明算法与传统算法对比误差分析图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
如图1所示,本发明提供的一种复杂地形时UHVAC输电线路下方地面三维工频电场计算方法,包括如下步骤:
a.建立三维电场计算坐标,将复杂地形简化成起伏形式的三维地面;
b.选定镜像地面,并假定镜像地面延伸至起伏地面内;
c.将带电导线离散成分段的线电荷,在非镜像地面内以正六边形蜂窝状形式设置离散点电荷,同时依据镜像原理设置相应的镜像电荷并建立三维电场积分数学模型的离散方程;
d.然后选定基函数为分域脉冲函数,按照配点法选权函数并在每个积分域求内积,形成矩阵方程;
e.求解矩阵方程,计算电荷分布。
f.计算匹配点最大电位误差,如果误差不满足要求,返回步骤c重新对场源离散;如果满足要求进行下一步。
g.利用优化方法对离散电荷开展进一步优化配置。
h.利用最优电荷分布计算起伏地面的三维工频电场分布。
本实施例步骤a中,将复杂地形简化成起伏形式的三维地面,是根据具体的地形进行简化,简化的原则是将复杂地形利用具有一定角度的多个平面来组合而成,如图2所示。
本实施例步骤b中选定的镜像地面为导线下方的地面,该地面为水平的一段地面,如图3所示。
本实施例步骤c中将带电导线离散成分段的线电荷,是考虑导线的弧垂,将导线分成足够小的线段,将小的线段看成直线电荷单元。
本实施例步骤c中在非镜像地面内部上以正六边蜂窝状形式设置离散的点电荷与地表面垂直设置,同时正六边蜂窝状的边长根据电场计算精度进行改变,如图4所示。
本实施例步骤c中建立的三维电场积分数学模型离散方程的场源包括导线电荷和非镜像地面内部电荷,同时方程还考虑了地形的影响:
其中为场点矢量,导线线源矢量,导线镜像线源矢量,非镜像地面内部点源矢量,镜像点源矢量,l为线电荷所存在区域,lj为镜像线电荷所在区域,q为点电荷。
本实施例步骤d中,选定的基函数为分域脉冲函数,如下:
本实施例步骤d中,按照配点法选择的权函数如下:
其中表示离散源点至计算点的矢径,表示空间任意一点至计算点的矢径,当时,ωj=∞,当当时,ωj=0。
本实施例步骤d中,矩阵方程的形成通过如下步骤实现:
首先,将式(2)代入式(1),得出电位积分方程的离散形式:
其次,在积分区域内对每个ωj与式(3)求内积,则有:
最后,根据δj函数性质及式(2)可以将(4)简化成下式:
令
则有:
本实施例中,上述矩阵方程的系数Pji通过如下方法实现:如图5所示,若以P1Pi为起点,沿PPiPPj方向建立局部坐标u,并设线单元长度为L0,设线电荷密度τ在单元内按线性规律分布则τ(u)=au+b(u的范围为(0,L0)),令l=xj-xi,m=yj-yi,n=zj-zi,Ei=l2+m2+n2,Fij=-2(l(xj-xi)+m(yj-yi)+n(zj-zi)),Gij=(xj-xi)2+(yj-yi)2+(zj-zi)2,对开展积分计算有:
本实施例中上述矩阵方程的系数Pji'采用式求出。
本实施例步骤f中,匹配点最大电位误差应当小于5%,如果大于5%,则重新布置离散电荷通过如下方法实现:根据校验电位误差对离散电荷的位置及个数进行设置:重新设置过程中,对误差较大的区域,以离散点电荷为正六边形的中心进行分裂,增加离散点电荷的数量来提高计算精度,如图4所示。
本实施例步骤g中利用优化方法对离散电荷开展进一步优化配置,首先建立目标函数,本发明的目标函数通过求解场域边界上所有匹配点的已知电位和计算电位的差值平方之和的最小值来表示:
其中,为所有离散电荷在第i个匹配点处的电位;为第i个匹配点的已知电位,导线表面φ=U,非镜像地面上φ=0;式(6)具有如下约束条件:S1.离散电荷的电量为自由变量;S2.离散电荷的位置必须在无效计算场域内:
zQ-f(zQ)<0i=m+1…n(8)
其中,式(7)和(8)中,xQd,yQd,zQd为导线上电荷坐标,xo,yo,zo为子导线中心坐标,r为子导线半径,zQ为起伏地面地下离散电荷纵坐标,m为子导线数,并且采用共轭梯度法求解式(6)中的极小值,从而求得最优的离散电荷。
本实施例步骤h中:根据最优离散电荷计算起伏地面三维工频电场分布通过如下方式实现:
设复杂地面某点为P(x,y,z),则该点的电场强度可由离散线电荷和离散点电荷产生的电场叠加获得。XYZ三个方向的电场强度由下式求出:
其中,n为离散线电荷个数,m为离散点电荷个数,L0为离散线电荷长度,A=aL0,B=b。x,y,z为待求点坐标,xi,yi,zi为线源坐标,xj,yj,zj为点源坐标。
则P点的电场强度有效值为:
以实际临近斜坡地面的500kV超高压交流输电线路为例,利用本发明方法进行三维工频电场计算,并与测试结果进行对比。该线路斜坡起点距离边导线9.5米,利用经纬仪测得该斜坡的角度为15.80°,输电导线对地高度为18m,ABC//CBA逆相序排列,最大导线间距26.5m,为8分裂线,分裂半径为0.40M,子输电导线的半径为0.0158M,。利用本发明方法计算的三维工频电场分布如图6所示,同时在该斜坡上选取了66个测点(每隔2米选取一个测点,形成11行6列矩阵,共计66个测点)进行了测试,计算结果和测试结果对比如图7所示。图7中可以得知计算结果和测试结果的最大误差为6.8%,由此证明本发明所提供的计算方法有效,可以用于工程实际。
同时,分别利用常规的计算方法(传统方法不设置镜像面,在整个计算面上都需要设置电荷)和本计算方法对该模型进行计算,传统计算方法中的校验点(校验点是用来校验计算误差)和本发明的校验点的最大误差φ曲线如图8所示,从图8中可以得出,依据本发明所提供的方法计算得出的电位误差明显小于传统算法中电位误差,由此证明本发明的计算方法有效,并且精度高,可以更好的满足工程实际需要。
通过以上实例证明,本发明的复杂地形情况下特高压交流输电线路三维工频电场计算方法明显优于传统算法。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (13)
1.一种复杂地形时UHVAC输电线路下方地面三维工频电场计算方法,其特征在于:包括以下步骤:
a.建立三维电场计算坐标,将复杂地形简化成起伏形式的三维地面;
b.选定镜像地面,并假定镜像地面延伸至起伏地面内;
c.将带电导线离散成分段的线电荷,在非镜像地面内以正六边形蜂窝状形式设置离散点电荷,同时依据镜像原理设置相应的镜像电荷并建立三维电场积分数学模型的离散方程;
d.然后选定基函数为分域脉冲函数,按照配点法选权函数并在每个积分域求内积,形成矩阵方程;
e.求解矩阵方程,计算电荷分布;
f.计算匹配点最大电位误差,如果误差不满足要求,返回步骤c重新对场源离散;如果满足要求进行下一步;
g.利用优化方法对离散电荷开展进一步优化配置;
h.利用最优电荷分布计算起伏地面的三维工频电场分布。
2.根据权利要求1所述复杂地形时UHVAC输电线路下方地面三维工频电场计算方法,其特征在于:步骤a中将复杂地形简化成起伏形式的三维地面,是根据具体的地形进行简化,简化的原则是将复杂地形采用具有一定角度的多个平面组合来代表。
3.根据权利要求2所述复杂地形时UHVAC输电线路下方地面三维工频电场计算方法,其特征在于:步骤b中选定的镜像地面为导线下方的地面,该地面为一段水平地面。
4.根据权利要求3所述复杂地形时UHVAC输电线路下方地面三维工频电场计算方法,其特征在于:步骤c中将带电导线离散成分段的线电荷,是考虑导线的弧垂,将导线分成足够小的线段,以满足将小的线段看成直线电荷单元。
5.根据权利要求4所述复杂地形时UHVAC输电线路下方地面三维工频电场计算方法,其特征在于:步骤c中在非镜像地面内部上以正六边蜂窝状形式设置离散的点电荷与地表面垂直设置,同时正六边蜂窝状的边长根据电场计算精度进行改变。
6.根据权利要求5所述复杂地形时UHVAC输电线路下方地面三维工频电场计算方法,其特征在于:步骤c中建立的三维电场积分数学模型离散方程的场源包括导线电荷和非镜像地面内部电荷,同时方程还考虑了地形的影响,其函数关系如下:
其中为场点矢量,导线线源矢量,导线镜像线源矢量,非镜像地面内部点源矢量,镜像点源矢量,l为线电荷所存在区域,lj为镜像线电荷所在区域,q为离散点电荷。
7.根据权利要求6所述复杂地形时UHVAC输电线路下方地面三维工频电场计算方法,其特征在于:步骤d中,选定的基函数为分域脉冲函数,其函数式如下:
8.根据权利要求7所述复杂地形时UHVAC输电线路下方地面三维工频电场计算方法,其特征在于:步骤d中,按照配点法选择的权函数如下:
其中表示离散源点至计算点的矢径,表示空间任意一点至计算点的矢径,当时,ωj=∞,当当时,ωj=0。
9.根据权利要求8所述复杂地形时UHVAC输电线路下方地面三维工频电场计算方法,其特征在于:步骤d中,矩阵方程的形成通过如下步骤实现:
首先,将式(2)代入式(1),得出电位积分方程的离散形式:
其次,在积分区域内对每个ωj与式(3)求内积,则有:
最后,根据δj函数性质及式(2)可以将(4)简化成下式:
令
则有:。
10.根据权利要求9所述复杂地形时UHVAC输电线路下方地面三维工频电场计算方法,其特征在于:所述矩阵方程的系数Pji通过如下方法实现:若以P1Pi为起点,沿PPiPPj方向建立局部坐标u,并设线单元长度为L0,设线电荷密度τ在单元内按线性规律分布则τ(u)=au+b(u的范围为(0,L0)),令l=xj-xi,m=yj-yi,n=zj-zi,Ei=l2+m2+n2,Fij=-2(l(xj-xi)+m(yj-yi)+n(zj-zi)),Gij=(xj-xi)2+(yj-yi)2+(zj-zi)2,对开展积分计算有:
11.根据权利要求10所述复杂地形时UHVAC输电线路下方地面三维工频电场计算方法,其特征在于:步骤f中,匹配点最大电位误差应当小于5%,如果大于5%,则重新布置离散电荷通过如下方法实现:根据校验电位误差对离散电荷的位置及个数进行设置:重新设置过程中,对误差较大的区域,以离散点电荷为正六边形的中心进行分裂,增加离散点电荷的数量来提高计算精度。
12.根据权利要求11所述复杂地形时UHVAC输电线路下方地面三维工频电场计算方法,其特征在于:步骤g中利用优化方法对离散电荷开展进一步优化配置,首先建立目标函数,本发明的目标函数通过求解场域边界上所有匹配点的已知电位和计算电位的差值平方之和的最小值来表示:
其中,为所有离散电荷在第i个匹配点处的电位;为第i个匹配点的已知电位,导线表面φ=U,非镜像地面上φ=0;式(6)具有如下约束条件:
S1.离散电荷的电量为自由变量;S2.离散电荷的位置必须在无效计算场域内:
zQ-f(zQ)<0i=m+1…n(8)
其中,在式(7)和式(8)中,xQd,yQd,zQd为导线上电荷坐标,xo,yo,zo为子导线中心坐标,r为子导线半径,zQ为起伏地面地下离散电荷纵坐标,m为子导线数,并且采用共轭梯度法求解式(6)中的极小值,从而求得最优的离散电荷。
13.根据权利要求12所述复杂地形时UHVAC输电线路下方地面三维工频电场计算方法,其特征在于:步骤h中根据最优离散电荷计算起伏地面三维工频电场分布通过如下方式实现:
设复杂地面某点为P(x,y,z),该点的电场强度可由离散线电荷和离散点电荷产生的电场叠加获得;XYZ三个方向的电场强度对应由以下公式求出:
其中,i为离散线电荷个数,j为离散点电荷个数,L0为离散线电荷长度,A=aL0,B=b;x,y,z为待求点坐标,xi,yi,zi为线源坐标,xj,yj,zj为点源坐标;则P点的电场强度有效值为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510909603.8A CN105427190B (zh) | 2015-12-09 | 2015-12-09 | 复杂地形时uhvac输电线路下方地面三维工频电场计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510909603.8A CN105427190B (zh) | 2015-12-09 | 2015-12-09 | 复杂地形时uhvac输电线路下方地面三维工频电场计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105427190A true CN105427190A (zh) | 2016-03-23 |
CN105427190B CN105427190B (zh) | 2020-06-12 |
Family
ID=55505372
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510909603.8A Active CN105427190B (zh) | 2015-12-09 | 2015-12-09 | 复杂地形时uhvac输电线路下方地面三维工频电场计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105427190B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105929251A (zh) * | 2016-04-13 | 2016-09-07 | 国网重庆市电力公司电力科学研究院 | Uhvac输电线路地面工频电场三维屏蔽分析与优化方法及系统 |
CN106168987A (zh) * | 2016-04-18 | 2016-11-30 | 国网重庆市电力公司电力科学研究院 | 一种输电线下复杂场景工频电场三维分析方法及系统 |
CN106680603A (zh) * | 2017-01-17 | 2017-05-17 | 清华大学 | 一种计算直流输电线路三维离子流场的特征线方法 |
CN106874565A (zh) * | 2017-01-17 | 2017-06-20 | 上海电力学院 | 一种雨天输电线路下方三维电场的计算方法 |
CN107063066A (zh) * | 2017-04-28 | 2017-08-18 | 国家电网公司 | 一种基于电场信息的输电线路弧垂反演方法 |
CN110533075A (zh) * | 2019-08-01 | 2019-12-03 | 贵州电网有限责任公司输电运行检修分公司 | 一种多模型电力线分割方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5751591A (en) * | 1994-11-09 | 1998-05-12 | Nec Corporation | Three dimensional wiring inductance calculation system |
CN103984876A (zh) * | 2014-06-11 | 2014-08-13 | 国家电网公司 | 特高压直流输电线路跨越复杂地形时的合成电场计算方法 |
-
2015
- 2015-12-09 CN CN201510909603.8A patent/CN105427190B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5751591A (en) * | 1994-11-09 | 1998-05-12 | Nec Corporation | Three dimensional wiring inductance calculation system |
CN103984876A (zh) * | 2014-06-11 | 2014-08-13 | 国家电网公司 | 特高压直流输电线路跨越复杂地形时的合成电场计算方法 |
Non-Patent Citations (2)
Title |
---|
张家利 等: "高压架空输电线下工频电场的数学模型", 《高电压技术》 * |
王凯奇: "超高压输电线路的三维工频电磁场计算及其影响因素分析", 《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105929251A (zh) * | 2016-04-13 | 2016-09-07 | 国网重庆市电力公司电力科学研究院 | Uhvac输电线路地面工频电场三维屏蔽分析与优化方法及系统 |
CN105929251B (zh) * | 2016-04-13 | 2018-09-04 | 国网重庆市电力公司电力科学研究院 | Uhvac输电线路地面工频电场三维屏蔽分析与优化方法及系统 |
CN106168987A (zh) * | 2016-04-18 | 2016-11-30 | 国网重庆市电力公司电力科学研究院 | 一种输电线下复杂场景工频电场三维分析方法及系统 |
CN106680603A (zh) * | 2017-01-17 | 2017-05-17 | 清华大学 | 一种计算直流输电线路三维离子流场的特征线方法 |
CN106874565A (zh) * | 2017-01-17 | 2017-06-20 | 上海电力学院 | 一种雨天输电线路下方三维电场的计算方法 |
CN106680603B (zh) * | 2017-01-17 | 2019-04-19 | 清华大学 | 一种计算直流输电线路三维离子流场的特征线方法 |
CN107063066A (zh) * | 2017-04-28 | 2017-08-18 | 国家电网公司 | 一种基于电场信息的输电线路弧垂反演方法 |
CN107063066B (zh) * | 2017-04-28 | 2019-06-18 | 国家电网公司 | 一种基于电场信息的输电线路弧垂反演方法 |
CN110533075A (zh) * | 2019-08-01 | 2019-12-03 | 贵州电网有限责任公司输电运行检修分公司 | 一种多模型电力线分割方法 |
CN110533075B (zh) * | 2019-08-01 | 2023-08-08 | 贵州电网有限责任公司 | 一种多模型电力线分割方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105427190B (zh) | 2020-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105427190B (zh) | 复杂地形时uhvac输电线路下方地面三维工频电场计算方法 | |
Wadhwa | High voltage engineering | |
CN103139907B (zh) | 一种利用指纹法的室内无线定位方法 | |
CN101392870B (zh) | 特高压交流线路对石油天然气管道影响的确定及防护方法 | |
Swetapadma et al. | A novel single-ended fault location scheme for parallel transmission lines using k-nearest neighbor algorithm | |
Kopietz et al. | Low-temperature behavior of the correlation length and the susceptibility of a quantum Heisenberg ferromagnet in two dimensions | |
CN108388707B (zh) | 一种三维非对称结构土壤模型下基于场路耦合的直流偏磁计算方法 | |
CN101894191B (zh) | 一种机车和牵引网与电力系统耦合的仿真方法 | |
CN101964525B (zh) | 一种支持大规模电流量测的配电网状态估计方法 | |
CN105354438B (zh) | 一种考虑地形影响的uhvdc线路地面三维合成电场计算方法 | |
CN106168987A (zh) | 一种输电线下复杂场景工频电场三维分析方法及系统 | |
CN103984876A (zh) | 特高压直流输电线路跨越复杂地形时的合成电场计算方法 | |
CN103995986A (zh) | 基于模拟电荷法的起伏地面架空输电线路的工频电场计算方法 | |
CN102590686B (zh) | 一种双极直流输电线路无线电干扰确定方法 | |
CN102280877B (zh) | 一种多量测断面的电力系统不良支路参数辨识方法 | |
CN107271941B (zh) | 一种同塔多回输电线路上的零序互感阻抗计算方法 | |
CN103455709A (zh) | 一种用于数字高程模型的超分辨率方法及其系统 | |
CN106649900B (zh) | 不均匀旋转对称体电磁特性的时域分析方法 | |
CN104597308A (zh) | 一种预测交叉互联式单芯四回路高压电缆护套环流的方法 | |
CN108388750B (zh) | 一种新型辐射状接地电阻计算方法 | |
CN107015285A (zh) | 一种观测视电阻率的校正方法及系统 | |
El Dein | Parameters affecting the charge distribution along overhead transmission lines’ conductors and their resulting electric field | |
CN109088796A (zh) | 基于网络层析成像技术的网络流量矩阵预测方法 | |
Gupta et al. | Coupling of high voltage AC power line fields to metallic pipelines | |
CN103116085B (zh) | 一种高压变电站内电场强度的计算方法 |
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 |