CN103605836A - Parallel computing method for three-dimensional electromagnetic fields of high-voltage transformer substation - Google Patents
Parallel computing method for three-dimensional electromagnetic fields of high-voltage transformer substation 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
The invention relates to a parallel computing method for three-dimensional electromagnetic fields of a high-voltage transformer substation. The parallel computing method is characterized in that an electromagnetic field model is established by using a three-dimensional finite element and a boundary element and can be solved; and the distribution condition of the electromagnetic fields of electromagnetic equipment of the high-voltage transformer substation can be sufficiently reflected by researching applicability of the finite element and applicability of the boundary element on computation of the electromagnetic fields of the transformer substation, and particularly, the distribution condition of the electromagnetic fields of the electromagnetic equipment such as a transformer comprising high-magnetic-permeability materials and the distribution condition of the electromagnetic fields of high-voltage electrical equipment such as an electric reactor can be reflected. The parallel computing method is quick and accurate in computation, is well-adapted, and has high engineering value.
Description
Technical field
The present invention is a kind of high voltage substation 3 D electromagnetic field parallel calculating method, is applied to the research of electromagnetic-field simulation method and the field analysis of high voltage substation and electromagnetic equipment.
Background technology
Along with greatly developing of China's high voltage power transmission engineering, electric pressure and the electric current of electrical network improve constantly, and the magnetic field that system and equipment produce also becomes and becomes increasingly complex.For reaching the indices of electromagnetic environment and guaranteeing the normal operation of system, the electromagnetic environment of research substation equipment seems particularly important.The electric charge of the power supply of the various electromagnetic equipments of transformer station and framework induction produces power frequency magnetic field at the near space in transformation device district, because transformer station's Charging system position is complicated and varied, except there being high-low voltage leads, also have the miscellaneous equipments such as transformer, high-voltage shunt reactor, grounding switch, capacitance type potential transformer and GIS, therefore, in transformer station, the power frequency magnetic field of electromagnetic equipment is that a complicated three dimensional field distributes.Research method about transformer station's electromagnetic environment has field measurement, simulation test and numerical evaluation both at home and abroad, comprises transmission line simulation corridor electric field, and computing equipment space electromagnetism distributes, and simplified model improves counting yield etc.Classic method is too single, calculated amount is huge, and the matrix of coefficients symmetry that coupling generates and orthotropicity are difficult to meet, when processing the tangential continuous problem of coupled field, there is defect, when device structure more complicated, especially, in the situation that equipment magnetic conduction medium consists of multiple material, with this class methods calculating magnetic field, can bring very large error.Ideal model is the complex apparatus of simulated field really, and versatility is not high, and following research need to reduce difficulty in computation meeting under the prerequisite of requirement of engineering precision.
Summary of the invention
The object of the invention is to overcome the deficiencies in the prior art, provide a kind of calculating fast, accurately, strong adaptability, has the high voltage substation electromagnetic field simulation method of higher actual application value.
The object of the invention is to be realized by following technical scheme:
A high voltage substation 3 D electromagnetic field parallel calculating method, is characterized in that, it comprises the following steps:
1) space electric field based on boundary element calculates
Put r then and there on the surface of conductor time, boundary integral equation can be expressed as
In formula, r is field point position vector; R ' is source point position vector;
for field point current potential; σ is source point surface charge density; R is the distance between a point and source point, for solving variable, border is discrete with σ, and adopts the golden weighted residual method of gal the Liao Dynasty, and formula (1) is write as:
In formula, e represents the numbering of a unit; E ' represents the numbering of source unit; I, j represents discrete nodes numbering; S
ethe integral domain that represents a unit, S
e 'the integral domain that represents source unit;
represent node potential; σ
irepresent node surface charge density; N
i, N
jfor interpolating function; ε
0for the specific inductive capacity of vacuum, discrete by interpolation, formula (2) is converted into algebraic equation, then tries to achieve Electric Field Distribution;
2) space magnetic field based on finite element calculates
Vector magnetic potential finite element method adopts vector magnetic potential A, ignores hysteresis effect and eddy effect, according to Maxwell equation, obtains magnetic field equation:
In formula, ▽ * be curl operator, the magnetic permeability that μ is medium, J is current density,
Take seamed edge finite element as example, and the whole interpolating function that solves field domain is
According to formula (1), to formula (2) application Green's theorem, obtain the golden Weighted Residual equation of gal the Liao Dynasty:
In formula, be vector dot product computing, disregard border item, equation arranges and is:
By formula (6) substitution formula (4), for whole weight functions, by the discrete formation Algebraic Equation set of Weighted Residual equation, solve the scalar magnetic potential A that obtains whole seamed edges,
Draw after the vector magnetic potential function of model, directly utilize the differential formulas in magnetic field, field, the magnetic field amount of computer memory any point, as magnetic flux density vector B and magnetic field intensity H, is shown in following formula,
Utilize the parallel computation of finite element-Element BEM to simulate high voltage substation 3 D electromagnetic field, and contrast with actual conditions, conclude efficiency and the error of finite element-Element BEM when simulation calculation electromagnetic equipment power frequency magnetic field, the calculating applicability of research finite element-Element BEM, for the calculating of high voltage substation complex apparatus periphery electromagnetic environment provides feasibility.
The present invention is a kind of high voltage substation 3 D electromagnetic field parallel calculating method, compare with existing computing method, the present invention utilizes three-dimensional finite element-boundary element set up electromagnetic field model and solve, the applicability when calculating transformer station's electromagnetic field by research finite element and boundary element, the electromagnetic field distribution situation that can fully reflect high voltage substation electromagnetic equipment, it is quick, accurate that the method is calculated, and strong adaptability, has higher actual application value.
Accompanying drawing explanation
Fig. 1 substation equipment schematic diagram;
Fig. 2 substation equipment model schematic diagram;
Fig. 3 substation equipment mirror image model schematic diagram;
Near field intensity maximal value cloud atlas Fig. 4 transformer station;
Fig. 5 electric field intensity is along the distribution plan of straight line L;
Fig. 6 core limb axon is to magnetic flux density vector distribution plan;
Fig. 7 cylinder unshakable in one's determination is magnetic flux density vector distribution plan radially.
embodiment
Below in conjunction with the drawings and specific embodiments, computing method of the present invention are further described:
As shown in Figure 1, substation equipment is the wire 8 on high-voltage shunt reactor 1, the first support insulator 2, lightning arrester 3, the second support insulator 4, capacitance type potential transformer 5, grounding switch 6, steel tower 7 and steel tower from right to left successively.By the electric field intensity near 1.5 meters, ground of distance plane boundary Element high-voltage shunt reactor 1, establish this plane and represent with P, as shown in Figure 2, set up the three-dimensional finite element magnetic field model of substation transformer.
1) space electric field based on boundary element calculates
Put r then and there on the surface of conductor time, boundary integral equation can be expressed as
In formula, r is field point position vector; R ' is source point position vector;
for field point current potential; σ is source point surface charge density; R is the distance between a point and source point, for solving variable, border is discrete with σ, and adopts the golden weighted residual method of gal the Liao Dynasty, and formula (1) is write as:
In formula, e represents the numbering of a unit; E ' represents the numbering of source unit; I, j represents discrete nodes numbering; S
ethe integral domain that represents a unit, S
e 'the integral domain that represents source unit;
represent node potential; σ
irepresent node surface charge density; N
i, N
jfor interpolating function; ε
0for the specific inductive capacity of vacuum, discrete by interpolation, formula (2) is converted into algebraic equation, then tries to achieve Electric Field Distribution;
2) space magnetic field based on finite element calculates
Vector magnetic potential finite element method adopts vector magnetic potential A, ignores hysteresis effect and eddy effect, according to Maxwell equation, obtains magnetic field equation:
In formula, ▽ * be curl operator, the magnetic permeability that μ is medium, J is current density,
Take seamed edge finite element as example, and the whole interpolating function that solves field domain is
According to formula (1), to formula (2) application Green's theorem, obtain the golden Weighted Residual equation of gal the Liao Dynasty:
In formula, be vector dot product computing, disregard border item, equation arranges and is:
By formula (6) substitution formula (4), for whole weight functions, by the discrete formation Algebraic Equation set of Weighted Residual equation, solve the scalar magnetic potential A that obtains whole seamed edges,
Draw after the vector magnetic potential function of model, directly utilize the differential formulas in magnetic field, field, the magnetic field amount of computer memory any point, as magnetic flux density vector B and magnetic field intensity H, is shown in following formula,
Utilize the parallel computation of finite element-Element BEM to simulate high voltage substation 3 D electromagnetic field, and contrast with actual conditions, conclude efficiency and the error of finite element-Element BEM when simulation calculation electromagnetic equipment power frequency magnetic field, the calculating applicability of research finite element-Element BEM, for the calculating of high voltage substation complex apparatus periphery electromagnetic environment provides feasibility.
Adopt the device structures such as finite element and boundary element model difference analog bus, transformer and grading ring, with other equipment such as surface model emulating reactance device grounding shells, the subdivision model of ABC three-phase equipment as shown in Figure 2.In actual computation process, be to keep the zero potential on ground, using ground as the plane of symmetry, to have set up the mirror image model of all devices, as shown in Figure 3.
In plane P, electric field intensity maximal value distributes as shown in Figure 4, can find out, in the symmetrical position of A, C, be on the phase position of limit, to have occurred the region that field intensity is larger, the conclusion of this and classic method is consistent, and A, 2 meters of left and right of C phase reactor axis are departed from respectively in this region.The maximal value of electric field intensity appears at the lower zone of C capacitance type potential transformer, and maximal value is 9.178kV/m.
On the straight line L apart from high resistance Ontology 20m place, carry out in-site measurement, on this straight line, the peaked measurement curve of electric field intensity and calculated curve are as shown in Figure 5, as seen from Figure 5, A, corresponding two crests in two limits of C, the corresponding middle trough of B, measures curve and simulation curve basically identical, and it is better that limit partly coincide mutually, this region that transformer station's power frequency electric field maximal value occurs is just the position of primary study.B mediates mutually, and surrounding devices is more, and the impact being subject to during measurement is larger, and therefore, curve is in trough part some deviation a little.
Set up simple iron core-coil former, the long 3m of cylinder unshakable in one's determination, radius 0.5m, adopts respectively two-dimensional axial symmetric face unit (PFEM), node finite element (NFEM), seamed edge finite element (EFEM), calculates different relative permeability μ
rthe distribution of Shi Tiexin axis z and radial line r magnetic field B.With reference to Fig. 6, core limb axon is to magnetic flux density vector distribution plan, the magnetic field result of wherein calculating based on two-dimensional axial symmetric face element method conforms to actual conditions, has higher accuracy, and the result of calculation of nodal element and Edge Finite Element Method is compared.With reference to Fig. 7, cylinder unshakable in one's determination is magnetic flux density vector distribution plan radially, and three kinds of computing method are constant.Under linear conditions, when iron core is identical with air permeability, the error of NFEM, EFEM and PFEM is very little; Along with μ
rincrease, EFEM is consistent with the result of PFEM, but the increase of NFEM error, while showing higher magnetic permcability, there is larger error in the result of calculation of NFEM.
The present invention is a kind of high voltage substation 3 D electromagnetic field parallel calculating method, and algorithm principle is clear and definite, and model is simple.By Example Verification near the method electromagnetic intensity calculation of complex equipment more accurately, for the calculating of high voltage substation complex apparatus periphery electromagnetic environment provides a kind of feasible method, realized the object of the invention, reached described effect.
Claims (1)
1. a high voltage substation 3 D electromagnetic field parallel calculating method, is characterized in that, it comprises the following steps:
1) space electric field based on boundary element calculates
Put r then and there on the surface of conductor time, boundary integral equation can be expressed as
In formula, r is field point position vector; R ' is source point position vector;
for field point current potential; σ is source point surface charge density; R is the distance between field point and source point.For solving variable, border is discrete with σ, and adopt the golden weighted residual method of gal the Liao Dynasty, formula (1) is write as:
In formula, e represents the numbering of a unit; E ' represents the numbering of source unit; I, j represents discrete nodes numbering; S
ethe integral domain that represents a unit, S
e 'the integral domain that represents source unit;
represent node potential; σ
irepresent node surface charge density; N
i, N
jfor interpolating function; ε
0for the specific inductive capacity of vacuum, discrete by interpolation, formula (2) is converted into algebraic equation, then tries to achieve Electric Field Distribution;
2) space magnetic field based on finite element calculates
Vector magnetic potential finite element method adopts vector magnetic potential A, ignores hysteresis effect and eddy effect, according to Maxwell equation, obtains magnetic field equation:
In formula, ▽ * be curl operator, the magnetic permeability that μ is medium, J is current density,
Take seamed edge finite element as example, and the whole interpolating function that solves field domain is
According to formula (1), to formula (2) application Green's theorem, obtain the golden Weighted Residual equation of gal the Liao Dynasty:
In formula, be vector dot product computing, disregard border item, equation arranges and is:
By formula (6) substitution formula (4), for whole weight functions, by the discrete formation Algebraic Equation set of Weighted Residual equation, solve the scalar magnetic potential A that obtains whole seamed edges,
Draw after the vector magnetic potential function of model, directly utilize the differential formulas in magnetic field, field, the magnetic field amount of computer memory any point, as magnetic flux density vector B and magnetic field intensity H, is shown in following formula,
Utilize the parallel computation of finite element-Element BEM to simulate high voltage substation 3 D electromagnetic field, and contrast with actual conditions, conclude efficiency and the error of finite element-Element BEM when simulation calculation electromagnetic equipment power frequency magnetic field, the calculating applicability of research finite element-Element BEM, for the calculating of high voltage substation complex apparatus periphery electromagnetic environment provides feasibility.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310536199.5A CN103605836A (en) | 2013-11-02 | 2013-11-02 | Parallel computing method for three-dimensional electromagnetic fields of high-voltage transformer substation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310536199.5A CN103605836A (en) | 2013-11-02 | 2013-11-02 | Parallel computing method for three-dimensional electromagnetic fields of high-voltage transformer substation |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103605836A true CN103605836A (en) | 2014-02-26 |
Family
ID=50124057
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310536199.5A Pending CN103605836A (en) | 2013-11-02 | 2013-11-02 | Parallel computing method for three-dimensional electromagnetic fields of high-voltage transformer substation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103605836A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106168987A (en) * | 2016-04-18 | 2016-11-30 | 国网重庆市电力公司电力科学研究院 | Complex scene power frequency electric field three dimensional analysis method and system under a kind of power transmission line |
CN107315890A (en) * | 2017-07-12 | 2017-11-03 | 英特工程仿真技术(大连)有限公司 | A kind of electromagnetic adsorption power computational methods for considering air gap layer influence |
CN108287268A (en) * | 2017-12-28 | 2018-07-17 | 华中科技大学 | A kind of distribution power transmission network power frequency electromagnetic field analysis method |
CN109408868A (en) * | 2018-09-12 | 2019-03-01 | 福州天宇电气股份有限公司 | A kind of transformer winding space electric field calculation method |
CN109902432A (en) * | 2019-03-14 | 2019-06-18 | 国网上海市电力公司 | A kind of power frequency magnetic field acquisition methods based on around ultra-high-tension power transmission line |
CN110135050A (en) * | 2019-05-10 | 2019-08-16 | 沈阳工业大学 | A kind of gas insulated transformer substation external electromagnetic field calculation method |
CN111209663A (en) * | 2019-12-31 | 2020-05-29 | 华中科技大学 | Modeling analysis method and system for equivalent ultra-long wave antenna array of high-voltage power grid |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103021024A (en) * | 2012-12-14 | 2013-04-03 | 国家电网公司 | Method displaying electromagnetic field of converting station in three-dimensional and visual mode |
CN103186689A (en) * | 2011-12-31 | 2013-07-03 | 英特工程仿真技术(大连)有限公司 | Electromagnetic field simulation analysis method |
-
2013
- 2013-11-02 CN CN201310536199.5A patent/CN103605836A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103186689A (en) * | 2011-12-31 | 2013-07-03 | 英特工程仿真技术(大连)有限公司 | Electromagnetic field simulation analysis method |
CN103021024A (en) * | 2012-12-14 | 2013-04-03 | 国家电网公司 | Method displaying electromagnetic field of converting station in three-dimensional and visual mode |
Non-Patent Citations (2)
Title |
---|
刘求龙 等: "三维激光扫描仪在惠泉变电站三维模型构建中的应用", 《测绘科学》 * |
李衷怡 等: "空间电磁场三维可视化技术研究", 《计算机与数字工程》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106168987A (en) * | 2016-04-18 | 2016-11-30 | 国网重庆市电力公司电力科学研究院 | Complex scene power frequency electric field three dimensional analysis method and system under a kind of power transmission line |
CN107315890A (en) * | 2017-07-12 | 2017-11-03 | 英特工程仿真技术(大连)有限公司 | A kind of electromagnetic adsorption power computational methods for considering air gap layer influence |
CN107315890B (en) * | 2017-07-12 | 2020-07-31 | 英特工程仿真技术(大连)有限公司 | Electromagnetic adsorption force calculation method considering air gap layer influence |
CN108287268A (en) * | 2017-12-28 | 2018-07-17 | 华中科技大学 | A kind of distribution power transmission network power frequency electromagnetic field analysis method |
CN108287268B (en) * | 2017-12-28 | 2019-07-23 | 华中科技大学 | A kind of distribution power transmission network power frequency electromagnetic field analysis method |
CN109408868A (en) * | 2018-09-12 | 2019-03-01 | 福州天宇电气股份有限公司 | A kind of transformer winding space electric field calculation method |
CN109902432A (en) * | 2019-03-14 | 2019-06-18 | 国网上海市电力公司 | A kind of power frequency magnetic field acquisition methods based on around ultra-high-tension power transmission line |
CN110135050A (en) * | 2019-05-10 | 2019-08-16 | 沈阳工业大学 | A kind of gas insulated transformer substation external electromagnetic field calculation method |
CN110135050B (en) * | 2019-05-10 | 2022-12-06 | 沈阳工业大学 | Calculation method for space electromagnetic field of gas insulated substation |
CN111209663A (en) * | 2019-12-31 | 2020-05-29 | 华中科技大学 | Modeling analysis method and system for equivalent ultra-long wave antenna array of high-voltage power grid |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103605836A (en) | Parallel computing method for three-dimensional electromagnetic fields of high-voltage transformer substation | |
CN107657137B (en) | Fractional order electromagnetic anomalous diffusion three-dimensional simulation method for rational function approximation | |
CN103810355B (en) | Transformer station's high-voltage switch gear field power frequency electric field three dimensional analysis method | |
CN106294914A (en) | A kind of converter transformer valve-side three-dimensional electric field emulation mode considering Anisotropic Nonlinear | |
CN103714197B (en) | A kind of construction design method optimizing ultra-high/extra-high voltage transmission line of electricity electromagnetic environment | |
CN108388707B (en) | Direct-current magnetic bias calculation method based on field-circuit coupling under three-dimensional asymmetric structure soil model | |
CN106199742A (en) | A kind of Frequency-domain AEM 2.5 ties up band landform inversion method | |
CN106168987A (en) | Complex scene power frequency electric field three dimensional analysis method and system under a kind of power transmission line | |
CN104361196A (en) | Arc net off-line electric arc mathematical model for calculating train speed | |
CN103018534A (en) | Method and system for determining harmonic voltage | |
CN103675565A (en) | Method for measuring power frequency parameters of superhigh/extrahigh-voltage alternating-current (direct-current) power transmission circuit | |
CN103258124B (en) | A kind of power frequency magnetic field assessment method for electromagnetic equipment of high voltage transformer substation | |
CN102590686B (en) | A kind of method for determining radio interference of bipolar direct current transmission line | |
Qiao et al. | An iterative flux tracing method without Deutsch assumption for ion-flow field of AC/DC hybrid transmission lines | |
CN104809297A (en) | Electromagnetic force density transferring method used among special-shaped grids in magnetic field-structure field coupling calculation | |
Gao et al. | Element differential method for solving linear and nonlinear electromagnetic problems | |
CN103699752A (en) | Coupling method for processing moving boundary problems in electromagnetic field based on edge element method | |
CN105912771A (en) | Ansys-based analysis method for optimization of influence of insulator detecting and cleaning robot on insulator string distribution voltage | |
CN102590766A (en) | Preprocessed conjugate gradient method of converter transformer magnetic field analysis | |
CN103472415B (en) | A kind of small scale test method of air reactor space magnetic field | |
CN103762057A (en) | Hollow reactor model simplifying method | |
Fukumoto et al. | Feasible method of calculating leakage reactance of 9-winding transformer for high-voltage inverter system | |
Li et al. | Research on Distortion Field of Inspection Robot around 500kV Substation Transformer | |
CN103258119A (en) | Assessment method for power frequency electric field of high voltage transformer substation | |
Lin et al. | Study on measurement errors of the ball type sensor in the power frequency electric field |
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 |