CN103810355B - 变电站高压开关场工频电场三维分析方法 - Google Patents
变电站高压开关场工频电场三维分析方法 Download PDFInfo
- Publication number
- CN103810355B CN103810355B CN201410090134.7A CN201410090134A CN103810355B CN 103810355 B CN103810355 B CN 103810355B CN 201410090134 A CN201410090134 A CN 201410090134A CN 103810355 B CN103810355 B CN 103810355B
- Authority
- CN
- China
- Prior art keywords
- mtd
- mrow
- msub
- mtr
- lambda
- 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
- 230000005684 electric field Effects 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000004141 dimensional analysis Methods 0.000 title claims abstract description 7
- 238000004088 simulation Methods 0.000 claims abstract description 42
- 239000004020 conductor Substances 0.000 claims abstract description 14
- 230000005611 electricity Effects 0.000 claims abstract description 10
- 239000012212 insulator Substances 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000005755 formation reaction Methods 0.000 claims description 3
- BDAGIHXWWSANSR-UHFFFAOYSA-M Formate Chemical compound [O-]C=O BDAGIHXWWSANSR-UHFFFAOYSA-M 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims 2
- 239000011800 void material Substances 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 15
- 238000011160 research Methods 0.000 abstract description 3
- 238000012360 testing method Methods 0.000 description 8
- 230000005540 biological transmission Effects 0.000 description 7
- 230000005672 electromagnetic field Effects 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明实施例提供了一种变电站高压开关场工频电场三维分析方法。其包括:a、对变电站开关场进行建模,得到变电站开关场简化模型,其中在步骤a中忽略变电站开关场中复杂设备及瓷套支柱,将变电站开关场所有的带电导线简化成水平和竖直两种形式,不考虑线路电压降;同时变电站中各设备均压环均简化为线型导体,最终得到简化的变电站开关场简化模型;b、采用模拟电荷法计算所述变电站开关场简化模型中的模拟电荷;c、根据计算的模拟电荷采用迭代法计算变电站开关场中各点的电场强度值,绘制变电站开关场的工频电场分布图。本发明实施例,可以为变电站电磁环境分析与研究提供技术参考。
Description
技术领域
本发明涉及电力分析技术领域,尤其涉及一种变电站高压开关场工频电场三维分析方法。
背景技术
近年来,国内变电站数量正逐年增加,而变电站内不仅因为带电设备众多,结构复杂,而且电压等级高,负荷重,场强大,易形成电磁污染,因此需要加强控制。针对工频电磁环境问题,国际非电离辐射防护委员会1988年发布导则规定了工频场强的限值,世界各国也相继规定了自己的国家标准限值,我国在参照国际非电离辐射防护委员导则的同时结合国情制定了工频电场和工频磁场的限值,以工频电场4kV/m,工频磁场0.1mT作为居民区环保评价标准,以工频电场10kV/m,工频磁场0.5mT作为职业卫生评价标准。目前,国内针对输电线路的工频电场的仿真计算相对较多,而针对变电站的计算因为其结构复杂相对较少,因此,本领域技术人员有必要提出一种变电站高压开关场工频电场分析方式。
发明内容
有鉴于此,本发明提供一种变电站高压开关场工频电场三维分析方法,可以进行多种不同状况下的变电站高压开关场工频电场计算分析,可以为变电站电磁环境分析与研究提供技术参考。
本发明提供了一种变电站高压开关场工频电场三维分析方法,包括:
a、对变电站开关场进行建模,得到变电站开关场简化模型,其中在步骤a中忽略变电站开关场中复杂设备及瓷套支柱,将变电站开关场所有的带电导线简化成水平和竖直两种形式,不考虑线路电压降;同时变电站中各设备均压环均简化为线型导体,最终得到简化的变电站开关场简化模型;
b、采用模拟电荷法计算所述变电站开关场简化模型中的模拟电荷;
c、根据计算的模拟电荷采用迭代法计算变电站开关场中各点的电场强度值,绘制变电站开关场的工频电场分布图。
进一步,所述步骤b包括:
b1、在所述变电站开关场简化模型中设置数量合适的模拟电荷和匹配点;
b2、建立电位和电荷量间的电位矩阵方程,求解该电位矩阵方程,得到设置的模拟电荷的量值;
b3、在所述变电站开关场简化模型中设置校验点,并计算所述校验点和模拟电荷间的电位误差;
b4、判断计算得到的电位误差是否达到设定的要求,若否则修正模拟电荷重复步骤b1至b4,直至满足限值,若是则执行步骤c。
进一步,所述步骤b2中的电位矩阵方程为:电位系数矩阵,该电位系数矩阵用于描述模拟电荷和匹配点之间的关系矩阵,其元素称为电位系数。
进一步,所述步骤b2中采用式:
计算设置的模拟电荷的量值;
其中,[U]—各导线对地电压的列向量,导线中包含架空导线和接地支柱,其中n为架空导线数目,m为接地支柱导体数目,其中U1、…、Un为n根架空导线的电压;Un+1、…、Un+m为m根支柱的对地电压,计算中视为零电位;
[Q]—各导线上等效电荷的列向量,其中Q1、…、Qn为n根架空导线的等效电荷;Qn+1、…、Qn+m为m根支柱上的等效电荷;
[λ]—各导线(包括架空导线和支柱)的电位系数矩阵,此矩阵为n+m阶方阵。
进一步,所述步骤b3中的电位误差包括:绝对误差,或者,包括绝对误差和相对误差。
本发明的有益效果:
本发明实施例针对变电站复杂的三维特性建立变电站开关场简化模型,然后采用模拟电荷法,计算变电站开关场中各点的电场强度值,绘制变电站开关场的工频电场分布图,从而为变电站电磁环境分析与研究提供技术参考。
附图说明
下面结合附图和实施例对本发明作进一步描述:
图1是场源外点P1的场量计算分析图。
图2是本发明提供的变电站高压开关场三维工频电场分析方法的实施例的流程示意图。
图3是分裂导线的等效示意图。
图4是对地电位计算图。
图5是计算值与实测值和对比图。
图6是500kV变电站开关场建模图。
具体实施方式
本发明实施例提出的变电站高压开关场工频电场三维仿真分析方法,是在变电站开关场简化的整体物理模型基础上,依据等效原理,将空间中连续分布的体电荷外的场强计算等效成对闭合边界面上的等效面电荷的面积分运算得到,如图1所示。
其中:其依据的电磁场数值计算方法是模拟电荷法。
具体的,在计算输电导线下的电场时,需要考虑的场源为带输电线表面的自由电荷和不同介质分界面上出现的极化电荷。由于这些电荷的分布情况未知,因此,根据电磁场的唯一性定理,按等效原理,将电极表面连续分布的自由电荷或介质分界面上的束缚电荷用一组离散的表面等效电荷(即模拟电荷)来等效代替,再应用叠加原理将离散模拟电荷在空间中产生的场量叠加,得到原连续分布电荷所产生的空间电场分布。
以上问题的数学模型为以电位函数为未知量的边值问题,即:
边界条件:
分界面衔接条件:
等效模拟电荷设立条件是边界条件和媒质分界面条件不变。在此基础上,按匹配点计算电位,建立模拟电荷与电位的电位方程组:
解此方程组可计算出等效电荷的电量,再通过设置一些校验点进行校验,得到合乎要求的等效电荷,依据这些等效电荷就可以近似计算出场域中任意一点的电位和电场强度。
下面具体说明本发明实施例提出的变电站高压开关场工频电场三维仿真分析方法。
如图2所示,本发明实施例提供的变电站高压开关场三维工频电场分析方法的实施例的流程示意图,其包括如下步骤:
步骤S21、变电站开关场简化建模。
具体的,忽略变电站开关场复杂设备及瓷套支柱,将变电站开关场所有的带电导线简化成水平和竖直两种形式,不考虑线路电压降;同时将断路器、隔离开关、电流互感器等设备上的均压环分别简化为各种线型导体,最终得到简化的整体物理模型。
进一步,将变电站开关场所有的带电导线简化成水平和竖直两种形式,同一水平导线取相同的离地高度。在此基础上,变电站内的工频电磁场计算即可在一种规则状态下上进行,且为准静态场问题。
对于分裂导线来说,将输电导线等效处理为圆柱导线,如图3所示。此时的等效圆柱导线半径可以采用等效半径Ri来代替,其计算公式如下:
式中:R为分裂导线半径,n为次导线根数,r为子导线半径。
导线电压已知,不考虑线路电压降;同时将断路器、隔离开关、电流互感器等设备上的均压环分别简化为各种线型导体。
同时,视地面为良导体,并取零电位。在输电线在大地的镜像位置上设置镜像模拟电荷,等效代替大地表面上感应电荷的影响。
如图6所示,是按照上述方式实现的500kV变电站开关场的建模图。
步骤S22、设置多个模拟电荷和匹配点,建立电位和电荷量间的电位矩阵方程,并求解该矩阵方程,得到模拟电荷的量值。
具体的,在设置模拟电荷和匹配点时:在计算场域设置数量合适的模拟电荷,并在相关的导体表面设置匹配点,其匹配点的设置根据线型来设置匹配点与模拟电荷的距离、相位来设定匹配点的电位及线高来设定匹配点对地高度等。
具体的,考虑导线上施加三相对称正弦电压,运用有效值相量表示电压,则输电线各相导线的相电压可表示成:
则图4中各相电位实部和虚部计算公式如下:
式中的θ为初相角。考虑到输电线路的实际运行情况,以额定电压的1.05倍作为计算电压,相电压按下式计算:
根据以上对输电线路系统的简化处理,会使复杂设备附近的工频电场仿真结果与实测值有一定的偏差,但计算模型的构建保证了工频电场分布的最大值在计算结果所控制的范围之内,以此来分析问题可充分地保障户外变电站内工频电场环境的合理性。
进一步,在用模拟电荷法计算导线上的模拟电荷时,采用电位系数法。将导线上的电荷都用集中在导线中心的线电荷表示,可得下列方程:
式中:
[U]—各导线对地电压的列向量,导线中包含架空导线和接地支柱(其中n为架空导线数目,m为接地支柱导体数目);其中U1、…、Un为n根架空导线的电压;Un+1、…、Un+m为m根支柱的对地电压,计算中视为零电位。
[Q]—各导线上等效电荷的列向量。其中Q1、…、Qn为n根架空导线的等效电荷;Qn+1、…、Qn+m为m根支柱上的等效电荷。
[λ]—各导线(包括架空导线和支柱)的电位系数矩阵,此矩阵为n+m阶方阵。
步骤S24、设置校验点并计算校验点与匹配点的电位误差,电位误差至少包括绝对误差,还可以包括相对误差。
步骤S25、判断电位误差是否达到设定的要求限值,如果判断结果为否,则重复步骤S23和步骤S24,直至满足设定的要求;如果判断结果为是,则执行步骤S26;
步骤S26、根据计算得到的模拟电荷计算三维工频场。
具体的,在相应位置求得的模拟电荷量值的基础上,可以利用电荷计算周围空间电场的迭代算法,计算开关场中各点的电场强度值。
具体的,计算出模拟电荷值,再由叠加原理计算出n+m根导线(包括架空导线和接地支柱)在空间任意点P(x,y,z)产生的电位:
在空间任意点P(x,y,z)产生的电场强度的分量分别为:
其中:
式中:
x、y、z——空间任意点P的坐标;
xi、yi、zi——模拟电荷i的坐标(i=1、2,…、m+n);
ri、ri'——模拟电荷i及其镜象至P点的距离。
按正弦电路的相量计算,将电压分为实部分量和虚部分量,模拟电荷亦存在实部和虚部分量,电场强度的各分量也有虚部分量和实部分量。即
式中:
EPx:P点电场强度的水平x分量;
EPy:P点电场强度的水平y分量;
EPz:P点电场强度的竖直z分量;
ExR:由各导线的实部电荷在P点产生的场强的水平x分量;
ExI:由各导线的虚部电荷在P点产生的场强的水平x分量;
EyR:由各导线的实部电荷在P点产生的场强的水平y分量;
EyI:由各导线的虚部电荷在P点产生的场强的水平y分量。
EzR:由各导线的实部电荷在P点产生的场强的竖直z分量;
EzI:由各导线的虚部电荷在P点产生的场强的竖直z分量。
于是P点的合成场强EP为:
本实施例,基于模拟电荷法,实现了变电站高压开关场工频电场三维分析,可以为变电站电磁环境分析与研究提供技术参考。
为此,本发明实施例再进行了如下验证实验。
在重庆地区某500kV变电站开关场选取测试路径,布置20个测试点进行现场测试,测试所用设备:意大利PMM公司生产PMM8053A便携式电磁场测试仪,测试过程中还利用数字式电子温湿度计(CW8056)记录了温度、湿度。测试工作按照《辐射环境保护管理导则电磁辐射监测仪器和方法》(HJ/TIO.2-1996)和《工频电场测量》(GBT12720)以及《高压交流架空送电线路、变电站工频电场和磁场测量方法》(DLT988-2005)进行。在设施正常工作时间内进行测量,测量时探头高度离地1.5m,每个测点连续测3次,每次测量时间大于15s,读取稳定状态的最大值。测试结果和户外变电站高压开关场工频电场三维仿真分析方法计算的结果进行了对比,如表1所示。除测点10可能因畸变误差相对较大,其余19个测点的误差均在5%以内,证明了该计算方法在工程实际中的有效性和实用性。具体的如表1和图5所示。
表1实测和计算对比结果
测点 | 测试值(V/m) | 计算值(V/m) | 误差(%) |
1 | 1412 | 1245 | 1.67 |
2 | 3848 | 4325 | -4.77 |
3 | 3477 | 3267 | 2.1 |
4 | 6570 | 6830 | -2.6 |
5 | 2441 | 2823 | -3.82 |
6 | 1167 | 1324 | -1.57 |
7 | 2696 | 2569 | 1.27 |
8 | 3036 | 2801 | 2.35 |
9 | 1667 | 2724 | -10.57 |
10 | 4075 | 4123 | -0.48 |
11 | 6143 | 6534 | -3.91 |
12 | 2751 | 2510 | 2.41 |
13 | 1584 | 1423 | 1.61 |
14 | 1355 | 1256 | 0.99 |
15 | 1173 | 1089 | 0.84 |
16 | 2685 | 2452 | 2.33 |
17 | 1433 | 1625 | -1.92 |
18 | 1707 | 1689 | 0.18 |
19 | 662 | 605 | 0.57 |
20 | 612 | 398 | 2.14 |
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (2)
1.一种变电站高压开关场工频电场三维分析方法,其特征在于:包括:
a、对变电站开关场进行建模,得到变电站开关场简化模型;在步骤a中忽略变电站开关场中复杂设备及瓷套支柱,将变电站开关场所有的带电导线简化成水平和竖直两种形式,不考虑线路电压降,同时变电站中各设备均压环均简化为线型导体,最终得到简化的变电站开关场简化模型;
b、采用模拟电荷法计算所述变电站开关场简化模型中的模拟电荷;
c、根据计算的模拟电荷采用迭代法计算变电站开关场中各点的电场强度值,绘制变电站开关场的工频电场分布图;
所述步骤b包括:
b1、在所述变电站开关场简化模型中设置数量合适的模拟电荷和匹配点;
b2、建立电位和电荷量间的电位矩阵方程,求解该电位矩阵方程,得到设置的模拟电荷的量值;
b3、在所述变电站开关场简化模型中设置校验点,并计算所述校验点和模拟电荷的电位误差;
b4、判断计算得到的电位误差是否达到设定的要求,若否则修正模拟电荷重复步骤b1至b4,直至满足限值,若是则执行步骤c;
所述步骤b2中采用式:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>U</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>U</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>U</mi>
<mi>n</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>U</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>U</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mn>1</mn>
<mi>n</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>n</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>n</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>n</mi>
<mi>n</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mi>n</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>)</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>)</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>)</mo>
<mi>n</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Q</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Q</mi>
<mi>n</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Q</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Q</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
计算设置的模拟电荷的量值;
其中,U—各导线对地电压的列向量,导线中包含架空导线和接地支柱,其中n为架空导线数目,m为接地支柱导体数目,其中U1、…、Un为n根架空导线的电压;Un+1、…、Un+m为m根支柱的对地电压,计算中视为零电位;
Q—各导线上等效电荷的列向量,其中Q1、…、Qn为n根架空导线的等效电荷;Qn+1、…、Qn+m为m根支柱上的等效电荷;
λ—各导线的电位系数矩阵,此矩阵为n+m阶方阵,其中,导线包括架空导线和支柱;
所述步骤b2中的电位矩阵方程为:电位系数矩阵,该电位系数矩阵用于描述模拟电荷和匹配点之间的关系矩阵,其元素称为电位系数;
所述步骤c包括:计算出模拟电荷值,再由叠加原理计算出n+m根导线在空间任意点P(x,y,z)产生的电位:
在空间任意点P(x,y,z)产生的电场强度的分量分别为:
其中:
<mrow>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msup>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>+</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
式中:
x、y、z——空间任意点P的坐标;
xi、yi、zi——模拟电荷i的坐标,其中,i=1、2,…、m+n;
ri、r′i——模拟电荷i及其镜象至P点的距离;
按正弦电路的相量计算,将电压分为实部分量和虚部分量,模拟电荷亦存在实部和虚部分量,电场强度的各分量也有虚部分量和实部分量;即
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>P</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msubsup>
<mi>E</mi>
<mrow>
<mi>x</mi>
<mi>R</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>E</mi>
<mrow>
<mi>x</mi>
<mi>I</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>P</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msubsup>
<mi>E</mi>
<mrow>
<mi>y</mi>
<mi>R</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>E</mi>
<mrow>
<mi>y</mi>
<mi>I</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>P</mi>
<mi>z</mi>
</mrow>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msubsup>
<mi>E</mi>
<mrow>
<mi>z</mi>
<mi>R</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>E</mi>
<mrow>
<mi>z</mi>
<mi>I</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
式中:
EPx:P点电场强度的水平x分量;
EPy:P点电场强度的水平y分量;
EPz:P点电场强度的竖直z分量;
ExR:由各导线的实部电荷在P点产生的场强的水平x分量;
ExI:由各导线的虚部电荷在P点产生的场强的水平x分量;
EyR:由各导线的实部电荷在P点产生的场强的水平y分量;
EyI:由各导线的虚部电荷在P点产生的场强的水平y分量。
EzR:由各导线的实部电荷在P点产生的场强的竖直z分量;
EzI:由各导线的虚部电荷在P点产生的场强的竖直z分量;
于是P点的合成场强EP为:
其中,导线包括架空导线和接地支柱。
2.如权利要求1所述的方法,其特征在于:所述步骤b3中的电位误差包括:绝对误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410090134.7A CN103810355B (zh) | 2014-03-12 | 2014-03-12 | 变电站高压开关场工频电场三维分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410090134.7A CN103810355B (zh) | 2014-03-12 | 2014-03-12 | 变电站高压开关场工频电场三维分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103810355A CN103810355A (zh) | 2014-05-21 |
CN103810355B true CN103810355B (zh) | 2018-01-30 |
Family
ID=50707119
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410090134.7A Active CN103810355B (zh) | 2014-03-12 | 2014-03-12 | 变电站高压开关场工频电场三维分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103810355B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103995986B (zh) * | 2014-06-11 | 2017-11-07 | 国家电网公司 | 基于模拟电荷法确定起伏地面架空输电线路的工频电场的方法 |
CN105186353B (zh) * | 2015-09-22 | 2019-02-19 | 国网上海市电力公司 | 一种基于三维模型的变电站动设备带电状态显示方法 |
CN105354438B (zh) * | 2015-12-09 | 2018-09-04 | 国网重庆市电力公司电力科学研究院 | 一种考虑地形影响的uhvdc线路地面三维合成电场计算方法 |
CN106168987A (zh) * | 2016-04-18 | 2016-11-30 | 国网重庆市电力公司电力科学研究院 | 一种输电线下复杂场景工频电场三维分析方法及系统 |
CN106291135A (zh) * | 2016-07-01 | 2017-01-04 | 重庆大学 | 输电线路的D‑dot电场传感器的设计方法 |
CN106570229A (zh) * | 2016-10-20 | 2017-04-19 | 中国电力科学研究院 | 一种换流站交流滤波器场的三维工频电场分析方法及系统 |
CN106570324A (zh) * | 2016-11-03 | 2017-04-19 | 合肥华义电气科技有限公司 | 一种变电站工频电场计算方法 |
CN112067908B (zh) * | 2020-08-20 | 2023-06-16 | 国网山东省电力公司电力科学研究院 | 变电站机器人测量工频电场时畸变电场的拟合方法及系统 |
CN112881818A (zh) * | 2021-01-15 | 2021-06-01 | 广州穗能通能源科技有限责任公司 | 电场强度测量方法、装置、计算机设备和存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101349720A (zh) * | 2008-08-28 | 2009-01-21 | 重庆电力科学试验研究院 | 高压架空输电线路工频电场屏蔽分析方法 |
CN103258119A (zh) * | 2013-04-19 | 2013-08-21 | 东北电力大学 | 一种高压变电站工频电场的评估方法 |
CN103258124A (zh) * | 2013-04-11 | 2013-08-21 | 东北电力大学 | 一种高压变电站电磁设备工频磁场评估方法 |
CN103453942A (zh) * | 2013-09-13 | 2013-12-18 | 国家电网公司 | 测量500kV变电站运行人员接触工频电场、磁场强度的方法 |
-
2014
- 2014-03-12 CN CN201410090134.7A patent/CN103810355B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101349720A (zh) * | 2008-08-28 | 2009-01-21 | 重庆电力科学试验研究院 | 高压架空输电线路工频电场屏蔽分析方法 |
CN103258124A (zh) * | 2013-04-11 | 2013-08-21 | 东北电力大学 | 一种高压变电站电磁设备工频磁场评估方法 |
CN103258119A (zh) * | 2013-04-19 | 2013-08-21 | 东北电力大学 | 一种高压变电站工频电场的评估方法 |
CN103453942A (zh) * | 2013-09-13 | 2013-12-18 | 国家电网公司 | 测量500kV变电站运行人员接触工频电场、磁场强度的方法 |
Non-Patent Citations (2)
Title |
---|
超高压输电线下建筑物临近区域三维电场仿真;汪泉弟等;《高压电器》;20131031;第49卷(第10期);第1-6页 * |
高压输电线、变电站电磁场环境测量方法;刘华麟;《中国优秀硕士论文全文数据库(硕士)工程科技II辑》;20051215(第8期);第2.1-2.2节、第4.2.1小节 * |
Also Published As
Publication number | Publication date |
---|---|
CN103810355A (zh) | 2014-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103810355B (zh) | 变电站高压开关场工频电场三维分析方法 | |
CN103018534B (zh) | 确定谐波电压的方法及系统 | |
CN107248736B (zh) | 一种配电网线路正序参数的在线辨识方法 | |
CN103714197B (zh) | 一种优化超/特高压输电线路电磁环境的结构设计方法 | |
CN101894191B (zh) | 一种机车和牵引网与电力系统耦合的仿真方法 | |
CN106168987A (zh) | 一种输电线下复杂场景工频电场三维分析方法及系统 | |
CN108959768A (zh) | 基于分区域模拟电荷法的复杂场景工频电场计算方法 | |
Jovicic et al. | A linear formulation for power system state estimation including RTU and PMU measurements | |
CN105720594A (zh) | 纯有功三相不平衡负荷补偿容量的补偿方法 | |
CN106786493A (zh) | 一种多馈入直流相互作用因子的实用计算方法 | |
CN108388707A (zh) | 一种三维非对称结构土壤模型下基于场路耦合的直流偏磁计算方法 | |
CN103605836A (zh) | 一种高压变电站三维电磁场并行计算方法 | |
CN104237677B (zh) | 220kV及以上电网单回输电线路零序参数在线辨识方法 | |
CN102590686B (zh) | 一种双极直流输电线路无线电干扰确定方法 | |
CN112784516B (zh) | 基于统一回路构建的高压直流输电直流偏磁水平计算方法 | |
CN105260583A (zh) | 一种特高压工频电磁场对人体生物效应的计算方法及系统 | |
CN107017631A (zh) | 一种基于线性电路的三相配电网可观测性分析方法 | |
CN104655940A (zh) | 一种500kV变电站工频电场和磁场的计算方法 | |
CN105912771A (zh) | 基于ansys的绝缘子检测与清扫机器人对绝缘子串分布电压影响优化的分析方法 | |
Long et al. | Computation of 3-D Electric Field on Buildings of Arbitrary Shape Near HVDC Transmission Lines | |
CN103116085B (zh) | 一种高压变电站内电场强度的计算方法 | |
Samadinasab et al. | The influence of moisture and temperature on the behavior of soil resistivity in earthing design using finite element method | |
CN103258119A (zh) | 一种高压变电站工频电场的评估方法 | |
Wang et al. | 3-D electric field computation with charge simulation method around buildings near HV transmission lines | |
Zhen et al. | Total electric field calculation of body model under HVDC transmission lines |
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 |