高压输电线下有建筑物时工频电场的数学模型_邹澎
高压输电线路附近电磁环境的预测模型和软件

ri , j ———第 j 根输电线和第 i 根输电线之间 的距离 ;
1999 年第 1 期 电 力 建 设
·7 ·
r′i , j ———第 j 根输电线的镜像和第 i 根输电线 之间的距离 。
1999 年第 1 期 电 力 建 设
·5 ·
高压输电线路附近电磁环境的预测模型和软件
Prediction Model and Software for Electromagnetic Environment near to HV Transmission Lines
W( b - x2 ) + W 2( b - x2) 2 + PD′22
(18)
W( b - x2) + W 2( b - x2) 2 + PD22
W = 1 + k2 , k = tg α
式中 τ ———线路上的等效线电荷密度 ;
α ———线路转角的一半 ;
PD1 , PD′1 ———由 P 点向 B E 段线路及其镜像
=
1~
m)
的
τ i
,
j
时
,
Pi , j
(
i
=
1~
n
,
j
=
1~
m)
点可选
在第 i 相导线第 j 段中点处的导线表面上 ,列出矩
阵方程 :
·8 ·
电 力 建 设 1999 年第 1 期
图 3 单相线路转弯处的电场
U1
U2 ⁝
电晕放电辐射场强的横向距离特性可以表示
为:
E ( d) = E0 + 20 K2 lg ( d0Πd)
并行计算高压输电线路周围电场

并行计算高压输电线路周围电场
厉天威;阮江军;吴田
【期刊名称】《电工技术学报》
【年(卷),期】2009(024)007
【摘要】实际工程应用中的大规模数值模拟需要借助于并行计算.介绍了适合并行计算的非重叠区域分解法以及并行计算的实现.为了更精确地计算,建立了500kV高压输电线路及其周围房屋的三维模型,并把该模型划分为多个(2~6)分区通过机群来分别实行并行求解.各分区内形成的线性方程组通过Krylov迭代求解器求解.结果表明通过并行计算能够精确地计算线路附近的房屋周围各点电场强度,高压输电线路周围的房屋对外界电场有一定的屏蔽作用.
【总页数】7页(P1-6,15)
【作者】厉天威;阮江军;吴田
【作者单位】广西电力工业勘察设计研究院,南宁,530023;武汉大学电气工程学院,武汉,430072;武汉大学电气工程学院,武汉,430072
【正文语种】中文
【中图分类】TM151
【相关文献】
1.高压输电线路周围工频电场畸变规律研究 [J], 林炬;韩晶晶
2.并行计算合成绝缘子串电压分布及金具表面电场强度 [J], 厉天威;阮江军;杜志叶;黄道春
3.特高压输电线路周围三维电场并行计算 [J], 周宏威;孙丽萍;包文泉;李奔亮
4.工频电场、磁场及高压输电线路的电场效应 [J], 毛文利;杨新村;李海平
5.消雷器静态电场分布及其并行计算 [J], 王世山;胡长玲
因版权原因,仅展示原文概要,查看原文内容请购买。
交流特高压输电线路电场环境分析

交流特高压输电线路电场环境分析作者:徐瑞宇来源:《电子世界》2013年第04期【摘要】交流特高压输电线路周围的电场环境影响到线路的安全运行和附近居民的身体健康。
本文利用向量有限元方法分析了1000kV输电线路的绝缘子串电场和电位分布,分析了导线水平布线时周围空间的工频电场分布。
给出了带均压环的绝缘子串上的电场、电位分布和导线不同高度时地面上方1.5m处水平方向的电场强度大小分布。
由计算结果得出均压环在绝缘子底部上方800mm时均压效果最佳。
导线正下方电场强度相对较小,电场强度最大值出现在边相导线外。
要使最大电场强度值控制在4kVm-1以下,导线高度不应低于35m.【关键词】电场;向量有限元法;绝缘子;架空线1.引言交流特高压输电线路大多采用复合绝缘子串。
复合绝缘子具有抗污闪性能好,重量轻,运行基本免维护,电晕小,电磁干扰小等优点。
在特高压输电线路中,在绝缘子金具表面,绝缘材料和周围空气中,都有很高的电场强度[1]。
复合绝缘子因为只有上、下两端有金属端头,并且尺寸不大,它们和中间绝缘体一起构成的主电容远小于瓷绝缘子串的主电容,因此,复合绝缘子的电场分布极不均匀,尤其是在高压端侧会出现强场区,当电场强度达到一定数值后,就会引起绝缘材料老化,击穿,金具表面发生电晕甚至电弧[1,2]。
因此计算和控制绝缘子表面电场和电位分布有重要的意义。
特高压输电导线产生的电磁环境问题也越来越受到人们的关注。
不同国家和地区对各种情况下电磁环境的标准不尽相同[1-2],国际公认的数值,当输电线路跨过居民区时,地面上方1.5米处的电场强度数值应在4kV/m以下。
对高压输电线产生的工频电场除采取现场测量[3]外,大都采用模拟电荷法计算[4-5]。
在采用模拟电荷法时应先确定模拟电荷的大小和位置,其计算精度也由模拟电荷的精确程度确定。
由于绝缘子串的电场可以近似简化为二维轴对称场,输电导线产生的电场可以简化为二维平行平面场,本文采用二维向量有限元方法,以导线三相电压向量为边界条件,求解周围空间的电压和电场强度分布。
超高压输电线路工频电场的研究及其改善措施

超高压输电线路工频电场的研究及其改善措施赵杰;王玲桃【摘要】To further study the power frequency electric field strength beneath the extra-high voltage(EHV)transmission line and its improvement measure,the parameters of the transmission line are calculated according to the theory of electromag-netic field to obtain voltage distribution along the line. Considering the voltages along the line are different,the power frequency electric field strength at 1.5 m beneath the transmission line is analyzed with charge simulation method. The charge simulation method was verified with ANSYS software. It is found from the verification results that the erection of just a part shielded wire can meet the requirement,which can save the cost. The reasonable position of the shielded wire was obtained with analysis. Fi-nally,the power frequency electric field curve of transmission line with a part erection shielded wire beneath it was gained by Matlab simulation. The research results show that part erection shielded wire can improve the power frequency electric field of the transmission lines.%为深入研究超高压输电线路下方工频电场强度和改善措施,首先根据电磁场理论对输电线路的参数进行计算,得出沿线电压的分布.考虑沿线电压不同后,采用模拟电荷法对输电线下方1.5 m处的工频电场强度进行分析,并利用ANSYS软件对模拟电荷法进行验证.结果发现只需部分架设屏蔽线就可以满足要求,这样也可以节省成本.然后分析得出架设屏蔽线的合理位置.最后利用Matlab仿真得出输电线路下方部分架设屏蔽线的工频电场曲线,研究结果表明,部分架设屏蔽线可以改善输电线路的工频电场.【期刊名称】《现代电子技术》【年(卷),期】2015(038)022【总页数】5页(P126-130)【关键词】模拟电荷法;参数计算;工频电场强度;部分架设;Matlab【作者】赵杰;王玲桃【作者单位】山西大学,山西太原 030013;山西大学,山西太原 030013【正文语种】中文【中图分类】TN958-34;TP391.4近年来,我国国民经济在快速发展,电力建设的投入也不断加大。
高压架空输电线下工频电场的数学模型_张家利

解式( 5) 求得各线上的电荷分布 S, 将 S 代入式
( 1) 可求得空间各点电位 5 。场强由下式求得
E = - ý 5 = - i55/ 5x - j 55/ 5y - k55/ 5 z
2 结 论
本文所述方法可计算任意相、任意分裂高压线 上等效电荷分布及线下空间电位、电场, 且精度比简 化成直线时有较大提高。
上的等效电荷, 再计算输电线附近的电位和场强。
1. 1 积分方程
原点为下支点的架空悬链线见图 1, 方程为[ 4] :
y=
( 1/ 2) tanB[ 1 +
sh( a( 2x - l ) / sh( a/ 2)
2l)
]
+
( lsh( ax / 2l ) sh( a( x - l ) / 2l) / ash( a/ 2) ) õ
K
K
∑ ∑∫ 5 =
5j =
j= 1
2 4PEj= 1
Sj Qj ( x , x ′, y ø, z ′) dx
( 1)
1. 2 用矩量法解积分方程
因导线上电荷分布不均匀, 用脉冲基函数和点
匹配法求解( 1) 式时, 将线路分成 i 段, 每段长 dl i, 故
K 根线 被分成 K i 段, 第 m 段的端 点坐 标为( x m,
x m+ 1 ) , 中点坐标为( x n , yn , z n ) 。脉冲函数定义为
Pm( y) =
1, 若 x m ≤ x ≤ x m+ 1 , 0, 其它。
将第 j 相线电荷 Sj 在脉冲基函数下展开, 即
i
∑ Sj ( x ) =
SmP m( x ) , ( j = 1, 2, …, K ) , ( 2)
高压直流输电线路标称电场强度计算模型

u s e s f i n i t e e l e me n t a n a l y s i s s o f t wa r e COM S OL Mu l t i p h y s i c s t o e s t a b l i s h a b u n d l e d c o n d u c t o r — o v e r h e a d g r o u n d wi r e mo d e l ,a s i mp l e b u n d l e d c o n d u c t o r mo d e l a n d a n e q u i v a l e n t c o n d u c t o r — o v e r h e a d g r o u n d wi r e mo d e l a c c o r d i n g t o c o mmo n s i mp l i f i e d c a l c u l a t i o n me t h o d .I t c a r r i e s o n s i mu l a t i o n f o r t h e s e t h r e e mo d e l s a n d t a k e s t h e b u n d l e d c o n d u c t o r — o v e r h e a d g r o u n d wi r e mo d e l a s r e f e r e n c e f o r c o mp a r i n g d i s t r i b u t i o n o f n o mi n a l e l e c t r i c f i e l d o f d i f f e r e n t mo d e l s .Fi n a l l y i t c o n c l u d e s t ha t t h e 0 一 v e r h e a d g r o u n d wi r e ha s s l i g h t i n f l u e n c e o n n o mi n a l e l e c t r i c f i e l d a r o u n d DC t r a n s mi s s i o n l i n e s ,wh i c h me a n s i t i s f e a s i b l e t o i g n o r e t h e s i mp l i f i e d mo d e l o f o v e r h e a d g r o u n d wi r e i n e n g i n e e r i n g a p p l i c a t i o n .Th e s i mp l i f i e d mo d e l f o r e q u i v a l e n t c o n d u c —
高压输电线路电场分析模型研究

高压输电线路电场分析模型研究
随着电力行业的发展,高压输电线路的建设越来越广泛。
然而,高压输电线路
在运行过程中会产生电场,对周围环境造成一定的影响。
因此,电场分析模型研究变得尤为重要。
电场分析模型主要分为静电场分析和动态电场分析。
静电场分析研究的是在稳
态工作条件下高压输电线路产生的电场分布规律。
动态电场分析则着重研究在过渡过程中,如开关合闸等瞬态工作条件下电场分布的变化规律。
首先,就静电场分析模型而言,一般采用有限元法或边界元法来进行解析。
有
限元法则是通过计算电荷分布在空间中的情况来计算电场分布情况。
而边界元法则是通过计算电场在导体表面的情况来计算电场分布情况。
这两种方法都有其特点,但是具体使用哪种方法需要根据实际情况进行选择。
其次,动态电场分析模型则需要考虑因变化发生的诸多因素,如电流改变、电
荷沉积、非线性电特性、导体移动等。
这就要求模型需要具有足够的灵活性,可以实时反映电场分布的变化情况。
比如,在开关合闸的过程中,由于电路拖延和电感作用,电流变化会很大,因此在考虑动态电场的情况下就需要考虑这些不同的因素。
除了上述两个主要的电场分析模型,还有一些次要的电场分析模型。
如针对高
空风速变化来研究电气断弧距离的模型;针对不同载荷条件下的电场分布变化情况来研究电场分析的模型;针对树木和建筑物的影响进行研究的模型等。
总之,无论是静态的还是动态的,电场分析模型都是很重要的。
对于电力行业
来说,建立高效而准确的模型,可以帮助在运行过程中尽可能的减小电场对周围环境的影响,保持高压输电线路安全稳定运行。
超高压输变电(500KV)建设项目工频电场强度预测方法

超高压输变电(500KV)建设项目工频电场强度预测方法【摘要】高压输变电工程运行过程中的工频高压电场、无线电干扰和噪声等可能对人与动物的行为、生殖系统以及神经系统等有着严重的影响,需要引起人们的重视。
因此,对超高压输变电建设项目附近的工频电场强度进行评估预测就显得尤为重要。
本文即对工频电场强度的预测方法进行简单的介绍。
【关键词】超高压输变电;工频电场强度;影响;预测一、前言随着我国工业企业规模的迅速扩大和居民生活用电量的快速增长,我国电力工业得到了很大的发展,线路的电压等级也逐渐提高,超高压输变电建设项目不断增加。
但是,通过研究表明,高压输变电工程运行过程中的工频高压电场、无线电干扰和噪声等对人和动物的行为、生殖系统以及神经系统有着严重的影响,工频高压电场对人体健康的影响。
因此,高压输变电设备引发的电磁污染受到越来越多研究人员的关注。
我国也颁布了相应的对工频电场的作业进行规范,对施工场所的工频电场强度与产生工频超高压电场的设备防护提出了要求。
在防护规范中,对超高压输变电建设项目附近工频电场强度的评估预测非常重要,这里,我们就对500KV超高压输变电建设项目的工频电场强度预测方法进行简单的介绍。
二、工频电场强度基本概念简介当输电线路两端加上电压以后,导线内就会产生电流,而在导线周围就会出现相应的电场。
当所加电压变为工频交变电压后,导体就会产生低频交变电荷,而在导线和地面之间的空间里就会产生一个低频工频电场。
超高压输变电建设项目的工频电场强度会受到线路电压等级、电流强度、导线传输性能和排列情况、以及地形与周围环境等多方面因素的影响。
工频电场在测量过程中需要使用专用的测量探头和工频电场测量仪器等。
测量仪器必须要经过相关计量管理部门的检定,并在规定有效期内进行检验。
对于测量过程中的施工规范要严格按照有关安全规定,穿戴相应防护设备。
超高压输变电建设项目的工频电场强度的测量需要注意以下几点要求:1、对于正常运行的高压送电线路工频电场强度的测量,测量地点要选择在地势平坦和远离其他线路、树木的空地上;2、工频电场测量仪要放置在距地面一到两米的上方位置处,具体高度需要根据实际情况进行相应调整,测量结果要记录准确;3、为了避免测量仪表支架泄漏电流情况的发生,工频电场测量环境的湿度值要求低于80%;4、测量时只需要对工频电场垂直于地面的分量进行测量即可。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
高压输电线下有建筑物时工频电场的数学模型M at hematical M odels f or Pow er Frequency Electric Field on Building under High Volt age T ransmission Line郑州大学电子工程系 邹 澎 (郑州450052) 收稿日期:1996-10-01【摘要】 本文利用模拟电荷法计算了高压输电线路下有建筑物时线路附近工频电场的分布,介绍了模拟电荷法的原理和具体的计算方法,并介绍了预测软件的设计思路。
【关键词】 模拟电荷法,电磁环境预测,工频电场Abstract In t his paper the power frequency elect ric field on building under high volt age t ransmission line was calculat ed by charge simulation met hod,of which the principle ,calculat ing method and t he predict ing soft ware were int roduced.Key Words charge simulat ion met hod prediction of elect romagnetic envi ronment power f requency elec-tric f ield0 引言参考文献[1]中,介绍了高压输电线附近无建筑物时计算工频电场的数学模型,如果输电线路下有建筑物,就改变了工频电场的边界条件和场域的分布,利用该模型介绍的方法计算就很困难了。
本文系利用模拟电荷法,计算高压输电线下有建筑物时,线路附近工频电场的计算方法。
模拟电塔法是利用一组假想的模拟电荷来等效代替导体或介质表面连续分布的电塔,满足给定的边界条件,理论依据是场的唯一性定理。
模拟电荷的数量和位置由所研究的问题中边界的形状决定,模拟电荷的大小利用给定的边界条件计算。
利用这些模拟电荷和已知的电荷分布即可用解析式计算空间任一点的电位和场强[2]。
下面通过一个例子,说明高压输电线下有建筑物时,利用模拟电荷法计算线路附近工频电场的具体方法。
三相单回路高压输电线路下有一栋建筑物,输电线架设的高度为12m,相间距离为6m,输电线的直径为15mm ,建筑物的几何尺寸为3m ×3m ×3m 。
如图1所示。
图1 高压输电线路下建筑物1 设置模拟电荷(n 个)如果没有建筑物,高压输电线下的工频电场可以用镜像法求解[1],高压输电线路下有建筑物时,建筑物表面与大地等电位(对于50Hz 的工频电场,建筑物可以看作是良导体)。
仍用镜像法求解,为了保证建筑物表面的电位为零,设置一组模拟电荷,如图2中模拟电荷的位置用・表示(均在建筑物表面的内侧或建筑物附近地面的下侧)。
模拟电荷的数量取决于计算的精度(最大允许误差),一般来说模拟电荷的数量越多,计算的精度越高,但是对计算机速度、容量的要求也越高。
图2 模拟电荷的位置(示意简图)2 选择匹配点(n 个)选择匹配点是为了利用给定的边界条件计算各—13—个模拟电荷的大小,匹配点一般选在与模拟电荷对应的边界面上,即建筑物表面或地面上,匹配点的数量一般和模拟电荷的数量相等。
利用迭加原理可以写出各相输电线和所有匹配点的电位:U1U k U k+1 U n+k =A1,1…A1,k A1,k+1…A1,n+kA k,1…A k,k A k,k+1…A k,n+kA k+1,1…A k+1,k A k+1,k+1…A k+1,n+kA n+k,1…A n+k,k A n+k,k+1…A n+k,n+kS1S kS k+1S n+k(1)式(1)中的下标1~k表示k相输电线,k+1~n+k 表示n个模拟电荷或匹配点,U i(i=1~k)是各相输电线的电位,U i(i=k+1~n+k)是各匹配点的电位(均为零),S i(i=1~k)是各相输电线上的等效线电荷密度,S i(i=k+1~n+k)是各模拟电荷的电量。
对于三相输电线路,k=3。
系数矩阵可以分为4个区域,第一个区域(A i,j,i=1~k,j=1~k)中A i,j=12P E0ln 2h iR i (i=j)(2)A i,j=12P E0ln r′ijr ij (i≠j)(3)A j,i=A i,j其中h i是第i根输电线的高度,R i是第i根输电线的半径,对于分裂导线,R i是等效半径[2]。
r ij是第j 根输电线和第i根输电线之间的距离,r′ij是第j根输电线的镜像和第i根输电线之间的距离。
第二个区域(A i,j,i=1~k,j=k+1~n+k)中,A i,j=14P E0・1r ij(4)其中r ij是第j个模拟电荷和第i根输电线之间的距离。
第3个区域(A i,j,i=k+1~n+k,j=1~k)中,A i,j=12P E0lnr′ijr ij(5)其中r ij是第j根输电线和第i个匹配点之间的距离,r′ij是第j根输电线的镜像和第i个匹配点之间的距离。
第4个区域(A i,j,i=k+1~n+k,j=k+1~n+k)中,A i,j=14P E0・1r ij(6)其中r ij是第j个模拟电荷和第i个匹配点之间的距离。
3 解线性方程组方程组(1)是一个高阶的线性方程组,一般由几百个方程组成,可用高斯-赛德尔迭代法解出n个模拟电荷的值。
4 校验为了检验由式(1)解出的n个模拟电荷是否满足原来的边界条件(建筑物表面和地面的电位为零),还需要选一些校验点(m个)通过计算进行验证。
校验点也选在建筑物表面和建筑物附近的地表面,但不与各匹配点重合,每个校验点的电位U i=∑kj=1S i2P E0lnr′ijr ij+∑n+kj=k+114P E0õS ir″ij(7)其中r ij是第j根输电线和第i个校验点之间的距离,r′ij是第j根输电线的镜像和第i个校验点之间的距离。
r″ij是第j个模拟电荷和第i个校验点之间的距离。
如果U i与给定的边界面上的电位U0(=0)的差值小于预先规定的最大允许误差:U i-U o<$(8)则计算的这组模拟电荷是有效的,否则应当调整模拟电荷的数量和位置,重新列出线性方程组式(1),然后重新求解,再利用校验点检验,直到满足式(8)为止。
5 计算空间的电位和场强如果求出的模拟电荷通过校验满足式(8),即可利用这组模拟电荷和输电线上的电荷及其镜像电荷计算高压输电线和建筑物附近电位和场强的分布。
计算电位的分布仍用式(7),把校验点换成空间场点即可。
电场强度可由式(9)计算:E i→=12P E0(∑kj=1S ir ijr ij→-∑kj=1S ir′ijr′ij→)+14P E0∑n+kj=k+1S ir″ij2r″ij→(9) r ij→由第j根输电线指向第i个场点,r′ij→由第j根输电线的镜像指向第i个场点,r″ij→由第j个模拟电荷指向第i个场点。
利用上述计算结果,绘出的高压输电线下等位面的分布图如图3所示,图3(a)是无建筑物时高压输电线下的等位面图,图3(b)是有建筑物时高压输电线下的等位面图。
输电线和建筑物的参数与图1中相同,输电电压为220kV。
6 预测软件简介利用本文中方法编制的预测软件框图(子菜单略),见图4所示,给定建筑物的几何尺寸,模拟电荷—14—(a ) 无建筑物时1—0.5kV;2—1.0k V;3—2.0k V;4—5.0kV;5—10kV;6—20kV;7—30kV;8—40kV;9—50kV(b) 有建筑物时1—0.5kV;2—1.0k V;3—2.0k V;4—5.0kV;5—10kV;6—20kV ;7—30kV ;8—40kV ;9—50k V图3 高压输电线下的等位面图和匹配点的位置可以利用软件自动生成。
求解线性方程组式(1)算出模拟电荷以后进行校验,如果不满足式(8),还需要调整模拟电荷的数量和位置后重新计算,直到满足式(8)为止。
然后把算出的各模拟电荷的电量和位置坐标存入文件,以备计算电位和场高压输电线附近的工频电场(线下有建筑物)↓主 菜 单1.内容提要2.自动生成模拟电荷和匹配点的位置3.计算模拟电荷4.计算电位U5.计算场强E6.使用方法7.退出↓子菜单(略)图4 预测软件框图强时调用。
7 结论7.1 模拟电荷法用于高压输电线下有建筑物时工频电场的计算、预测,适当选择模拟电荷的数量和位置,可以达到预定的计算精度。
7.2 进行电磁环境预测,计算的场域很大,利用模拟电荷法计算时需要解一个高阶的线性方程组,我们比较了几种数值解法,利用高斯-赛德尔迭代法比较好。
8 参考文献1 邹澎.高压输电线附近工频电场的数学模型.中国电力,1994年;(6):P22~232 盛剑霓等.电磁场数值分析.北京:科学出版社,1984年:P 346~384(责任编辑 伍佐之)・会议报道・第二届全国电力通信专业学术会议在宜昌召开 中国电机工程学会电力通信专业委员会第二届学术会议暨第二届全体委员工作会议于1996年11月7日至8日在湖北宜昌召开。
出席这次会议的有专委会委员、论文作者、省电力局等代表共91人。
国电通信中心主任、电力通信专委会主任委员丁道齐教授介绍了世界通信发展和电信体制改革状况,指出了我国电力通信存在的问题和今后的发展方向。
专委会副主任委员鲍伟廉教授代表第一届委员会作了工作总结,国电通信中心副总工程师、国际大电网会议第35委员会委员曹惠彬同志作了关于国际大电网会议有关活动情况的报告。
本次会议共征集到论文65篇,经专家评审,评选出48篇编入论文集。
会议期间,共交流了论文36篇,与会代表就传输、交换、接入网、电力线载波、数据通信、网络管理等7个方面进行了广泛的学术交流和探讨。
会议认为:本届学术会议论文数量较多,内容广泛,且有一定深度。
主要研讨了世界通信的最新技术,如SDH 、AT M 及多媒体通信等,基本反映了我国电力通信科学技术领域的实际状况和水平。
这届学术活动达到了相互交流、相互促进、共同提高的目的。
代表们希望专委会在举办综合性学术研讨会议的同时,可根据生产实践的需要,组织灵活多样的专题学术活动,进一步提高学术交流的深度和广度。
当前,电力通信专委会应重点协助做好电力通信专业有关技术标准、规范的制定工作。
(金李莎)—15—第30卷 高压输电线下有建筑物时工频电场的数学模型 1997年第3期 。