CN103605836A - 一种高压变电站三维电磁场并行计算方法 - Google Patents
一种高压变电站三维电磁场并行计算方法 Download PDFInfo
- Publication number
- CN103605836A CN103605836A CN201310536199.5A CN201310536199A CN103605836A CN 103605836 A CN103605836 A CN 103605836A CN 201310536199 A CN201310536199 A CN 201310536199A CN 103605836 A CN103605836 A CN 103605836A
- Authority
- CN
- China
- Prior art keywords
- formula
- integral
- field
- equation
- magnetic
- 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.)
- Pending
Links
Images
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明是一种高压变电站三维电磁场并行计算方法,其特点是:利用三维有限元—边界元建立电磁场模型并求解,通过研究有限元和边界元在计算变电站电磁场时的适用性,能够充分反映高压变电站电磁设备的电磁场分布情况,能够充分反映高压变电站电磁设备的磁场分布情况,尤其是如变压器等包含高导磁材料的电磁设备和电抗器等高压电气设备,该方法计算快速、准确,适应性强,具有较高的工程价值。
Description
技术领域
本发明是一种高压变电站三维电磁场并行计算方法,应用于高压变电站及电磁设备的电磁场仿真方法研究以及场分析。
背景技术
随着我国高压输电工程的大力发展,电网的电压等级和电流不断提高,***及设备所产生的磁场也变得越来越复杂。为达到电磁环境的各项指标并保证***的正常运行,研究变电站设备的电磁环境显得尤为重要。变电站各种电磁设备的电源和构架感应的电荷在变电装置区的临近空间产生工频磁场,由于变电站带电装置位置复杂且多种多样,除了有高低压引线,还有变压器、高压并联电抗器、接地开关、电容式电压互感器和GIS等其它设备,因此,变电站内电磁设备的工频磁场是一个复杂的三维场分布。国内外关于变电站电磁环境的研究方法有现场实测、模拟试验和数值计算,包括模拟输电线路走廊电场,计算设备空间电磁分布,以及简化模型提高计算效率等。传统方法过于单一,计算量庞大,并且耦合生成的系数矩阵对称性和正定性很难满足,在处理耦合场切向连续问题时存在缺陷,当设备结构比较复杂时,尤其是设备导磁媒质由多种材料构成的情况下,用此类方法计算磁场会带来很大的误差。理想模型不能真正模拟现场的复杂设备,通用性不高,未来的研究需要在满足工程精度要求的前提下降低计算难度。
发明内容
本发明的目的是克服现有技术的不足,提供一种计算快速、准确,适应性强,具有较高实际应用价值的高压变电站电磁场模拟方法。
本发明的目的是由以下技术方案来实现的:
一种高压变电站三维电磁场并行计算方法,其特征是,它包括以下步骤:
1)基于边界元的空间电场计算
当场点r在导体的表面上时,边界积分方程可表示为
式中,e表示场单元的编号;e‘表示源单元的编号;i,j表示离散节点编号;Se表示场单元的积分区域,Se′表示源单元的积分区域;表示节点电位;σi表示节点面电荷密度;Ni,Nj为插值函数;ε0为真空的介电常数,通过插值离散,式(2)转化为代数方程,继而求得电场分布;
2)基于有限元的空间磁场计算
矢量磁位有限元法采用矢量磁位A,忽略磁滞效应和涡流效应,根据Maxwell方程得到磁场方程:
式中,▽×为旋度算子,μ为媒质的磁导率,J为电流密度,
以棱边有限元为例,求解场域的整体插值函数为
根据式(1)对式(2)应用格林定理,得伽辽金加权余量方程:
式中,·为矢量点积运算,不计边界项,方程整理为:
将式(6)代入式(4),针对全部权函数,将加权余量方程离散形成代数方程组,求解得到全部棱边的标量磁位A,
得出模型的矢量磁位函数后,直接利用磁场的微分公式,计算空间任意一点的磁场场量,如磁感应强度矢量B和磁场强度H,见下式,
利用有限元—边界元方法并行计算对高压变电站三维电磁场进行模拟,并与实际情况进行对比,归纳有限元—边界元方法在仿真计算电磁设备工频磁场时的效率和误差,研究有限元—边界元方法的计算适用性,为高压变电站复杂设备周围电磁环境的计算提供可行性。
本发明是一种高压变电站三维电磁场并行计算方法,与现有的计算方法相比,本发明利用三维有限元—边界元建立电磁场模型并求解,通过研究有限元和边界元在计算变电站电磁场时的适用性,能够充分反映高压变电站电磁设备的电磁场分布情况,该方法计算快速、准确,适应性强,具有较高的实际应用价值。
附图说明
图1变电站设备示意图;
图2变电站设备模型示意图;
图3变电站设备镜像模型示意图;
图4变电站附近场强最大值云图;
图5电场强度沿直线L的分布图;
图6铁心柱体轴向磁感应强度矢量分布图;
图7铁心柱体径向磁感应强度矢量分布图。
具体实施方式
下面结合附图和具体实施例对本发明的计算方法作进一步描述:
如图1所示,变电站设备从右向左依次是高压并联电抗器1、第一支柱绝缘子2、避雷器3、第二支柱绝缘子4、电容式电压互感器5、接地开关6、铁塔7及铁塔上的导线8。用边界元计算高压并联电抗器1附近距离地面1.5米平面上的电场强度,设此平面用P表示,如图2所示,建立变电站变压器的三维有限元磁场模型。
1)基于边界元的空间电场计算
当场点r在导体的表面上时,边界积分方程可表示为
式中,e表示场单元的编号;e‘表示源单元的编号;i,j表示离散节点编号;Se表示场单元的积分区域,Se′表示源单元的积分区域;表示节点电位;σi表示节点面电荷密度;Ni,Nj为插值函数;ε0为真空的介电常数,通过插值离散,式(2)转化为代数方程,继而求得电场分布;
2)基于有限元的空间磁场计算
矢量磁位有限元法采用矢量磁位A,忽略磁滞效应和涡流效应,根据Maxwell方程得到磁场方程:
式中,▽×为旋度算子,μ为媒质的磁导率,J为电流密度,
以棱边有限元为例,求解场域的整体插值函数为
根据式(1)对式(2)应用格林定理,得伽辽金加权余量方程:
式中,·为矢量点积运算,不计边界项,方程整理为:
将式(6)代入式(4),针对全部权函数,将加权余量方程离散形成代数方程组,求解得到全部棱边的标量磁位A,
得出模型的矢量磁位函数后,直接利用磁场的微分公式,计算空间任意一点的磁场场量,如磁感应强度矢量B和磁场强度H,见下式,
利用有限元—边界元方法并行计算对高压变电站三维电磁场进行模拟,并与实际情况进行对比,归纳有限元—边界元方法在仿真计算电磁设备工频磁场时的效率和误差,研究有限元—边界元方法的计算适用性,为高压变电站复杂设备周围电磁环境的计算提供可行性。
采用有限元和边界元模型分别模拟母线、变压器和均压环等设备结构,用面模型模拟电抗器接地外壳等其他设备,ABC三相设备的剖分模型如图2所示。实际计算过程中,为保持地面的零电位,以地面作为对称面建立了所有设备的镜像模型,如图3所示。
平面P上电场强度最大值分布见图4所示,可以看出,在A、C相对称的位置即边相位置上出现了场强比较大的区域,这和传统方法的结论是一致的,该区域分别偏离A、C相电抗器 轴线2米左右。电场强度的最大值出现在C电容式电压互感器的下方区域,最大值为9.178kV/m。
在距高抗本体约20m处的直线L上进行现场测量,该直线上电场强度最大值的测量曲线和计算曲线如图5所示,由图5可以看出,A、C两个边相对应两个波峰,B相对应中间波谷,测量曲线和仿真曲线基本一致,边相部分吻合的较好,这正是变电站工频电场最大值出现的区域,是重点研究的位置。B相处于中间位置,周围设备较多,测量时受到的影响较大,因此,曲线在波谷部分稍微有些偏差。
建立简单铁心—线圈模型,铁心柱体长3m,半径0.5m,分别采用二维轴对称面单元(PFEM)、节点有限元(NFEM)、棱边有限元(EFEM),计算不同相对磁导率μr时铁心轴线z和径线r磁场B的分布。参照图6,铁心柱体轴向磁感应强度矢量分布图,其中基于二维轴对称面单元方法计算的磁场结果与实际情况相符,具有较高的准确性,将节点有限元法与棱边有限元法的计算结果进行比较。参照图7,铁心柱体径向磁感应强度矢量分布图,三种计算方法不变。线性条件下,当铁心与空气磁导率相同时,NFEM、EFEM与PFEM的误差很小;随着μr增大,EFEM与PFEM的结果一致,但NFEM误差增大,表明较高磁导率时,NFEM的计算结果出现较大的误差。
本发明是一种高压变电站三维电磁场并行计算方法,算法原理明确,模型简单。通过算例验证了该方法能比较准确的计算复杂设备附近的电磁场强度,为高压变电站复杂设备周围电磁环境的计算提供了一种可行的方法,实现了本发明目的,达到了所述的效果。
Claims (1)
1.一种高压变电站三维电磁场并行计算方法,其特征是,它包括以下步骤:
1)基于边界元的空间电场计算
当场点r在导体的表面上时,边界积分方程可表示为
式中,e表示场单元的编号;e‘表示源单元的编号;i,j表示离散节点编号;Se表示场单元的积分区域,Se′表示源单元的积分区域;表示节点电位;σi表示节点面电荷密度;Ni,Nj为插值函数;ε0为真空的介电常数,通过插值离散,式(2)转化为代数方程,继而求得电场分布;
2)基于有限元的空间磁场计算
矢量磁位有限元法采用矢量磁位A,忽略磁滞效应和涡流效应,根据Maxwell方程得到磁场方程:
式中,▽×为旋度算子,μ为媒质的磁导率,J为电流密度,
以棱边有限元为例,求解场域的整体插值函数为
根据式(1)对式(2)应用格林定理,得伽辽金加权余量方程:
式中,·为矢量点积运算,不计边界项,方程整理为:
将式(6)代入式(4),针对全部权函数,将加权余量方程离散形成代数方程组,求解得到全部棱边的标量磁位A,
得出模型的矢量磁位函数后,直接利用磁场的微分公式,计算空间任意一点的磁场场量,如磁感应强度矢量B和磁场强度H,见下式,
利用有限元—边界元方法并行计算对高压变电站三维电磁场进行模拟,并与实际情况进行对比,归纳有限元—边界元方法在仿真计算电磁设备工频磁场时的效率和误差,研究有限元—边界元方法的计算适用性,为高压变电站复杂设备周围电磁环境的计算提供可行性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310536199.5A CN103605836A (zh) | 2013-11-02 | 2013-11-02 | 一种高压变电站三维电磁场并行计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310536199.5A CN103605836A (zh) | 2013-11-02 | 2013-11-02 | 一种高压变电站三维电磁场并行计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103605836A true CN103605836A (zh) | 2014-02-26 |
Family
ID=50124057
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310536199.5A Pending CN103605836A (zh) | 2013-11-02 | 2013-11-02 | 一种高压变电站三维电磁场并行计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103605836A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106168987A (zh) * | 2016-04-18 | 2016-11-30 | 国网重庆市电力公司电力科学研究院 | 一种输电线下复杂场景工频电场三维分析方法及*** |
CN107315890A (zh) * | 2017-07-12 | 2017-11-03 | 英特工程仿真技术(大连)有限公司 | 一种考虑气隙层影响的电磁吸附力计算方法 |
CN108287268A (zh) * | 2017-12-28 | 2018-07-17 | 华中科技大学 | 一种分布式输电网络工频电磁场分析方法 |
CN109408868A (zh) * | 2018-09-12 | 2019-03-01 | 福州天宇电气股份有限公司 | 一种变压器绕组空间电场计算方法 |
CN109902432A (zh) * | 2019-03-14 | 2019-06-18 | 国网上海市电力公司 | 一种基于高压输电线路周围的工频磁场获取方法 |
CN110135050A (zh) * | 2019-05-10 | 2019-08-16 | 沈阳工业大学 | 一种气体绝缘变电站空间电磁场计算方法 |
CN111209663A (zh) * | 2019-12-31 | 2020-05-29 | 华中科技大学 | 一种高压电网的等效超长波天线阵建模分析方法和*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103021024A (zh) * | 2012-12-14 | 2013-04-03 | 国家电网公司 | 用于对变电站中的电磁场进行三维可视化展示的方法 |
CN103186689A (zh) * | 2011-12-31 | 2013-07-03 | 英特工程仿真技术(大连)有限公司 | 电磁场仿真分析方法 |
-
2013
- 2013-11-02 CN CN201310536199.5A patent/CN103605836A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103186689A (zh) * | 2011-12-31 | 2013-07-03 | 英特工程仿真技术(大连)有限公司 | 电磁场仿真分析方法 |
CN103021024A (zh) * | 2012-12-14 | 2013-04-03 | 国家电网公司 | 用于对变电站中的电磁场进行三维可视化展示的方法 |
Non-Patent Citations (2)
Title |
---|
刘求龙 等: "三维激光扫描仪在惠泉变电站三维模型构建中的应用", 《测绘科学》 * |
李衷怡 等: "空间电磁场三维可视化技术研究", 《计算机与数字工程》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106168987A (zh) * | 2016-04-18 | 2016-11-30 | 国网重庆市电力公司电力科学研究院 | 一种输电线下复杂场景工频电场三维分析方法及*** |
CN107315890A (zh) * | 2017-07-12 | 2017-11-03 | 英特工程仿真技术(大连)有限公司 | 一种考虑气隙层影响的电磁吸附力计算方法 |
CN107315890B (zh) * | 2017-07-12 | 2020-07-31 | 英特工程仿真技术(大连)有限公司 | 一种考虑气隙层影响的电磁吸附力计算方法 |
CN108287268A (zh) * | 2017-12-28 | 2018-07-17 | 华中科技大学 | 一种分布式输电网络工频电磁场分析方法 |
CN108287268B (zh) * | 2017-12-28 | 2019-07-23 | 华中科技大学 | 一种分布式输电网络工频电磁场分析方法 |
CN109408868A (zh) * | 2018-09-12 | 2019-03-01 | 福州天宇电气股份有限公司 | 一种变压器绕组空间电场计算方法 |
CN109902432A (zh) * | 2019-03-14 | 2019-06-18 | 国网上海市电力公司 | 一种基于高压输电线路周围的工频磁场获取方法 |
CN110135050A (zh) * | 2019-05-10 | 2019-08-16 | 沈阳工业大学 | 一种气体绝缘变电站空间电磁场计算方法 |
CN110135050B (zh) * | 2019-05-10 | 2022-12-06 | 沈阳工业大学 | 一种气体绝缘变电站空间电磁场计算方法 |
CN111209663A (zh) * | 2019-12-31 | 2020-05-29 | 华中科技大学 | 一种高压电网的等效超长波天线阵建模分析方法和*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103605836A (zh) | 一种高压变电站三维电磁场并行计算方法 | |
CN107657137B (zh) | 一种有理函数逼近的分数阶电磁反常扩散三维模拟方法 | |
CN103810355B (zh) | 变电站高压开关场工频电场三维分析方法 | |
CN106294914A (zh) | 一种考虑非线性各向异性的换流变压器阀侧三维电场仿真方法 | |
CN103714197B (zh) | 一种优化超/特高压输电线路电磁环境的结构设计方法 | |
CN108388707B (zh) | 一种三维非对称结构土壤模型下基于场路耦合的直流偏磁计算方法 | |
CN106168987A (zh) | 一种输电线下复杂场景工频电场三维分析方法及*** | |
CN103018534A (zh) | 确定谐波电压的方法及*** | |
CN103675565A (zh) | 超/特高压交(直)流输电线路工频参数测量方法 | |
CN103258124B (zh) | 一种高压变电站电磁设备工频磁场评估方法 | |
CN102590686B (zh) | 一种双极直流输电线路无线电干扰确定方法 | |
Qiao et al. | An iterative flux tracing method without Deutsch assumption for ion-flow field of AC/DC hybrid transmission lines | |
Yue et al. | Electric field equivalent method for AC/DC hybrid electric equipment adapted to digital twin with low latency | |
Gao et al. | Element differential method for solving linear and nonlinear electromagnetic problems | |
CN105912771A (zh) | 基于ansys的绝缘子检测与清扫机器人对绝缘子串分布电压影响优化的分析方法 | |
CN103699752A (zh) | 基于棱边单元法处理电磁场中移动边界问题的耦合方法 | |
CN103762057A (zh) | 一种空心电抗器模型简化方法 | |
CN103472415A (zh) | 一种空心电抗器空间磁场的缩比实验方法 | |
CN102590766A (zh) | 一种换流变磁场分析的预处理的共轭梯度方法 | |
Fukumoto et al. | Feasible method of calculating leakage reactance of 9-winding transformer for high-voltage inverter system | |
Guangzheng et al. | FEM analysis of 3-D transformer leakage field and eddy current loss in the windings | |
Li et al. | Research on Distortion Field of Inspection Robot around 500kV Substation Transformer | |
CN103258119A (zh) | 一种高压变电站工频电场的评估方法 | |
Liu et al. | Non-contact current measurement method for multi-conductor in digital twin based on magnetoresistance sensors and PSO | |
Filippov et al. | The secondary sources method analysis and experimantal modeling of the permanent magnet eddy currents suspansion |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140226 |
|
WD01 | Invention patent application deemed withdrawn after publication |