CN108287268A - A kind of distribution power transmission network power frequency electromagnetic field analysis method - Google Patents
A kind of distribution power transmission network power frequency electromagnetic field analysis method Download PDFInfo
- Publication number
- CN108287268A CN108287268A CN201711464089.7A CN201711464089A CN108287268A CN 108287268 A CN108287268 A CN 108287268A CN 201711464089 A CN201711464089 A CN 201711464089A CN 108287268 A CN108287268 A CN 108287268A
- Authority
- CN
- China
- Prior art keywords
- transmission network
- electromagnetic field
- power transmission
- frequency electromagnetic
- power frequency
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/0864—Measuring electromagnetic field characteristics characterised by constructional or functional features
- G01R29/0892—Details related to signal analysis or treatment; presenting results, e.g. displays; measuring specific signal features other than field strength, e.g. polarisation, field modes, phase, envelope, maximum value
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Abstract
The invention discloses a kind of distributed power transmission network power frequency electromagnetic field analysis methods, belong to electromagnetic field technology field.The present invention builds distributed power transmission network model according to the voltage and current on the node, loop and loop of true distributed power transmission network;Domain model is emulated further according to the environment construction where true distributed power transmission network;The power frequency electromagnetic field that further simulation calculation distribution power transmission network generates in space, by the visualization for emulating data to power frequency electromagnetic field, to disclose the rule that distributed power transmission network generates power frequency electromagnetic field, by the divergence and curl of computational EM waves, propagation law of the power frequency electromagnetic field of distributed power transmission network generation under different medium is further disclosed.The method of the present invention is suitable for the power frequency electromagnetic field emulation of complicated distributed power transmission network, has filled up related blanking technique.
Description
Technical field
The invention belongs to electromagnetic field technology fields, more particularly, to a kind of distributed power transmission network power frequency electromagnetic field point
Analysis method.
Background technology
Due to the fast development of science and technology, the application of electromagnetic technique creates huge material progress to the mankind, but simultaneously
Also an environment for being full of artificial electromagnetic radiation is brought.With the development of electric system, the increase of power grid scale and capacity, electricity
The raising of grade is pressed, electric network composition is more complicated, and the electromagnetic field that electric system generates has been to be concerned by more and more people.
Power transmission network will produce periodically variable power frequency electromagnetic field, as the aortic distributed grid of power grid
Network can generate powerful electromagnetic field to its neighbouring space, form powerful electromagnetic interference.The power frequency electromagnetic field that power transmission network generates
The certain point in space can be traveled to by the modulated attenuation of medium (such as water body, air, soil), the power frequency electromagnet at the point
Field is that the power transmission network of different location distribution generates the vector sum of power frequency electric field.
Influence of the ultra-high-tension power transmission line to environment is various, the ultra-high-tension power transmission line electromagnetic environment usually described in us
Problem includes following a few problems:The ecological effect of power frequency electromagnetic field, the electromagnetic interference etc. that power frequency electromagnetic field is brought.Power frequency electromagnetic field
Ecological effect be always focus of people's attention.Power frequency electromagnetic field is shown from many report on Epidemiological both domestic and external
May there may be relationship with leukaemia, cancer morbidity raising on human immunity function there is certain influence.Distribution transmission of electricity
The electromagnetic interference that Netowrk tape comes mainly influences the radio broadcasting of surrounding resident and listening to and watching for TV.
In the past, since the transmission line of electricity in China is less, voltage class is relatively low, and the overwhelming majority is erected at far from populous
The deserted mountain wilderness in area, so the electromagnetic problems of power transmission network do not attract people's attention.With the development of power industry, then
In addition the influence of the social factors such as expanding economy and population growth, urban population and supply load all compare concentration, high pressure
The electromagnetic radiation that power transmission network generates more and more is closed as a kind of new pollution sources in city by the members of all social strata
Note.With China's development of the West Regions, transferring electricity from the west to the east, the construction of Three Gorges Pivot Project power plant project can will also create the height of substantial amounts
Pressure transmission line in powerful power transmission to China, eastern region and big city therefore study distributed power transmission network production
Raw power frequency electromagnetic field just seems very urgent and important.
The power frequency electromagnetic field emulation of existing ultra-high-tension power transmission line carries out numerical value meter just in single transmission line of electricity
It calculates.Power frequency electromagnetic field emulation for complicated distributed power transmission network, both at home and abroad there is no literature reported on.
Invention content
For the disadvantages described above or Improvement requirement of the prior art, the present invention provides a kind of distributed power transmission network power frequency electrics
Magnetic field analytical method, its object is to establish the distributed power transmission network model on polycyclic road, multinode, multivoltage, foundation includes
The emulation domain model of water layer and air layer, the power frequency electromagnetic field that simulation calculation distribution power transmission network generates in space lead to
The visualization for emulating data to power frequency electromagnetic field is crossed, to disclose the rule that distributed power transmission network generates power frequency electromagnetic field, is passed through
The divergence and curl of computational EM waves further disclose the power frequency electromagnetic field of distributed power transmission network generation under different medium
Propagation law, thus break through the problem of prior art can not emulate the power frequency electromagnetic field of complicated distributed power transmission network.
To achieve the above object, described the present invention provides a kind of distributed power transmission network power frequency electromagnetic field analysis method
Method includes:
(1) according to the voltage and current on the node, loop and loop of distributed power transmission network, distributed power transmission network is built
Model;
(2) domain model is emulated according to the environment construction where distributed power transmission network;
(3) according to the distributed power transmission network under distributed power transmission network model and emulation domain model structure boundary condition constraint
Network power frequency electromagnetic field computation model goes out distributed power transmission network in the power frequency electromagnet for entirely emulating domain by finite element stimulation
Field spatial distribution;
(4) by the form of power frequency electromagnetic field spatial distribution drafting pattern, visual figure is analyzed.
Further, the step (1) specifically includes:
(11) nodal analysis method of distributed power transmission network is built
V=(V1,V2,…,Vi,…,Vm),
Wherein, ViIndicate a node of power transmission network, node ViSpace coordinate be (xi,yi,zi), m is indicated in network
Number of nodes;
(12) loop model of distributed power transmission network is built
P=(P1,P2,…,Pi,…,Pk),
Wherein, PiIndicate that a loop in power transmission network, k indicate the loop quantity in power transmission network;Each loop PiAll
Including:Start node Vs, terminal node VeWith transmission of electricity loop normalLoop PiThe power frequency electric that any point (x, y, z) generates in space
Magnetic field isThe power frequency electromagnetic field that all loops are generated in spatial point (x, y, z) is:
(13) the loop voltage current model of distributed power transmission network is built
E=(e1,e2,…,ei,…,ek),
I=(i1,i2,…,ii,…,ik),
Wherein, eiIndicate loop PiOn voltage, iiIndicate loop PiOn electric current;
(14) distributed power transmission network model is built
G=(G1,G2,…,Gi,…,Gp),
Wherein, Gi=(V, P, e, i), GiIndicate i-th of sub-network among power transmission network.
Further, the step (2) specifically includes:
(21) structure emulation domain geometrical model, geometrical model include:
Air model, Mair=(Lair,Wair,Hair);
Extra large land model, Mocean-island=(Locean-island,Wocean-island,Hocean-island);
Wherein, (Lair,Wair,Hair) and (Locean-island,Wocean-island,Hocean-island) respectively indicate Air model and
The length and width and height of extra large land model;
(22) attribute model in structure emulation domain, attribute model include:
Air properties, (mur=1, epsilon=1, sigma=0);
Ocean attribute, (mur=1, epsilon=81, sigma=3);
Land attribute, (mur=1, epsilon=11.7, sigma=0);
Wherein, mur indicates that the relative permeability of material, epsilon indicate that the relative dielectric constant of material, sigma indicate
The relative conductivity of material.
Further, the step (3) is specially:
Using distributed power transmission network model as the electromagnetic excitation source of simulation calculation, emulation domain model is as simulation calculation
Space boundary constrains;Mesh generation is carried out to emulation domain, the entire space that calculates is divided into multiple sub-grids, is based on Maxwell
Electromagnetic field descriptive equation passes through the power frequency electromagnetic field of finite element stimulation sub-grid, the power frequency electromagnetic field structure of all-ones subnet lattice
At the power frequency electromagnet field distribution in entire emulation domain.
Further, the step (4) specifically includes:
(41) streamline figure layer is created, the number of streamline is set, by the electric field E (x, y, z, t) and magnetic field H of power frequency electromagnetic field
(x, y, z, t) draws motion pattern as data source;
(42) setting height z is chosen0, calculate the electric field magnitude and magnetic field amplitude of power frequency electromagnetic field under the height:
Wherein, Ex,t、Ey,t、Ez0,tThree components of the electric field under rectangular coordinate system, H are indicated respectivelyx,t、Hy,t、Hz0,tRespectively
Indicate three components of the magnetic field under rectangular coordinate system;It creates section and draws figure layer, by the electric field magnitude and magnetic of power frequency electromagnetic field
Field amplitude draws sectional drawing as data source;
(43) setting height z is chosen0, the electric field magnitude and magnetic field amplitude of power frequency electromagnetic field under the height are calculated, is created
The number of isopleth is arranged in isoplethes drawing figure layer, adjusts the type of isopleth, by the electric field magnitude of power frequency electromagnetic field and magnetic field
Amplitude is as data source drawing isoline figure;
(44) power frequency electromagnetic field divergence and curl are calculated:
Magnetic field divergence,
Magnetic field curl,
Wherein,Indicate divergence operator;Indicate curl operator;Indicate the unit vector of x-axis direction;Indicate y-axis
The unit vector in direction;Indicate the unit vector in z-axis direction;(x, y, z, t) was indicated between the three-dimensional coordinate added-time of point;
It creates section and draws figure layer, using the electromagnetic field divergence of power frequency electromagnetic field as data source, draw sectional drawing;Create stream
The number of streamline is arranged in line chart layer, using the electromagnetic field curl of power frequency electromagnetic field as data source, draws motion pattern.
In general, through the invention it is contemplated above technical scheme is compared with the prior art, have following technology special
Sign and advantageous effect:
(1) distributed power transmission network power frequency electromagnetic field modeling method proposed by the present invention, by being established to distributed network
Polycyclic road, multinode model, system modelling and emulation can be carried out to arbitrarily complicated power transmission network, solve the prior art
It may be only available for the modeling and simulating of simple transmission line of electricity;
(2) distributed power transmission network power frequency electromagnetic field modeling method proposed by the present invention, by polycyclic road to foundation, more
The power transmission network model of node carries out finite element stimulation, can in space be generated with the power transmission network of simulation analysis complexity
Power frequency electromagnetic field, solve the prior art can not the entire power transmission network of simulation analysis power frequency electromagnetic field distribution problem, Ke Yiying
The electromagnetic interference problem brought to surrounding enviroment for analyzing power transmission network can also be applied to equipment electricity in research electric system
Magnetic compatibling problem;
(3) power frequency electromagnet proposed by the present invention is field visualized and its analysis method, the sky obtained to finite element stimulation
Between power frequency electromagnetic field distributed data handled, draw the power frequency electromagnetic field sectional drawing and isogram of certain elevation plane, draw
Power frequency electromagnetic field motion pattern in space calculates the divergence and curl of electric field and magnetic field, it is established that a set of space power frequency electromagnetic field
Analysis method, be equally applicable to other complexity distributed driving source generate electromagnetic field.
Description of the drawings
Fig. 1 is the flow diagram of the present invention;
Fig. 2 is power transmission network illustraton of model in the present invention;
Fig. 3 is air illustraton of model in the present invention;
Ocean-land illustraton of model in the present invention of the positions Fig. 4;
Fig. 5 is that model meshes divide figure in the present invention;
Fig. 6 is power frequency magnetic field motion pattern in the present invention;
Fig. 7 is power frequency electric field motion pattern in the present invention;
Fig. 8 is magnetic field sectional drawing (observed altitude 10km) in the present invention;
Fig. 9 is electric field sectional drawing (observed altitude 10km) in the present invention;
Figure 10 is magnetic field sectional drawing (observed altitude -0.1km) in the present invention;
Figure 11 is electric field sectional drawing (observed altitude -0.1km) in the present invention;
Figure 12 is magnetic field isogram (observed altitude 10km) in the present invention;
Figure 13 is electric field isogram (observed altitude 10km) in the present invention;
Figure 14 is magnetic field isogram (observed altitude -0.1km) in the present invention;
Figure 15 is electric field isogram (observed altitude -0.1km) in the present invention;
Figure 16 is electric field divergence sectional drawing (observed altitude 10km) in the present invention;
Figure 17 is magnetic field divergence sectional drawing (observed altitude 10km) in the present invention;
Figure 18 is magnetic field curl motion pattern in the present invention;
Figure 19 is electric field curl motion pattern in the present invention;
Specific implementation mode
In order to make the purpose , technical scheme and advantage of the present invention be clearer, with reference to the accompanying drawings and embodiments, right
The present invention is further elaborated.It should be appreciated that described herein, specific examples are only used to explain the present invention, not
For limiting the present invention.As long as in addition, technical characteristic involved in the various embodiments of the present invention described below that
Conflict is not constituted between this to can be combined with each other.
Defeated/change/power utilization network system of large number of high/medium/low pressure is a mass of in world wide, these complicated electricity
Force system network is made of several ac transmission loops, and each loop of transmitting electricity includes several nodes, is passed through between loop of transmitting electricity
Transformer is attached, and thus constitutes the distributed power grid systems of multinode, polycyclic road, multivoltage.These transmissions of electricity
Loop can regard frequency electromagnetic radiation source one by one as, and the alternating electromagnetic field that frequency electromagnetic radiation source generates can be divided into property not
With two parts, a portion energy of electromagnetic field periodically flows back between radiation source surrounding space and radiation source
It moves, referred to as near field region;Another part energy of electromagnetic field is detached from radiator, in the form of an electromagnetic wave transmitting, referred to as far-field region outward.
The flow of the present invention as shown in Figure 1, foundation including distributed power transmission network model, the foundation for emulating domain model,
The analysis of power frequency electromagnetic field finite element simulation and power frequency electromagnetic field and etc..Wherein specific implementation includes the following steps:
(1) according to true power grid be distributed, establish a polycyclic road, multinode, multivoltage distributed grid
Network model.By taking the electric system of China as an example, national power grid is divided into six bulk power grids, including Northeast China Power Grid, North China Power Telecommunication Network, Central China
Power grid, East China Power Grid, Northwest Grid and south electric network include several sub-networks, each sub-network packet again below each power grid
Large number of transmission of electricity loop and node are included, is connected by transformer between loop, to realize the electricity between different transmission of electricity loops
Pressure conversion, thus constitutes a distributed power transmission network to spread all over the country.In this example, to the transmission of electricity of Central China Area distribution
Network establishes Central China Power Grid model.
(1.1) distributed high-voltage transmission network nodal analysis method is established.In this example, main pivot in Central China Power Grid is chosen
Node of the knob city as model.Origin under using Wuhan as rectangular coordinate system builds the position coordinates of power transmission network node.
Nodal analysis method in Central China Power Grid indicates with set V, wherein V=(V1,V2,…,Vi,…,Vm), ViIndicate the pivot of Central China Power Grid
Knob urban node, each node ViSpace coordinate be (xi,yi,zi), m indicates the number of nodes in whole network.
(1.2) distributed high-voltage transmission network loop model is established.In this example, according to the transmission of electricity city of Central China Power Grid
Connection relation builds the loop model of Central China Power Grid.The loop of Central China Power Grid indicates with set E, wherein E=(E1,E2,…,
Ei,…,Ek), EiIndicate that a loop in power transmission network, k indicate the loop quantity in network.Each transmission of electricity loop includes three
A element is start node V respectivelys, terminal node VeWith transmission of electricity loop normalSuch as there is one between Wuhan and Huang gang
Item transmission of electricity loop, then VsIt is set as Wuhan, VeIt is set as Huang gang, loop approximation is calculated by the coordinate of the two nodes
The normal of plane
(1.3) distributed high-voltage transmission network voltage/current model is established.There are one transmissions of electricity for each contact ring road
Voltage and transmission current, loop EiOn voltage ei, electric current iiIt indicates.For example, there are a voltages between Wuhan and Huang gang
For 110kV, the transmission of electricity loop E of 3kAi, then ei=110kV, ii=3kA, it is similar, construct all contact rings of Central China Power Grid
The Voltage-Current model on road.
(1.4) distributed high-voltage transmission network model is established.By the three above steps, loop, node and electricity are established
The model of pressure.These loops, node and voltage constitute multiple transmission of electricity sub-networks, with graph model Gi=(V, E) is described, wherein
GiIndicate that i-th of sub-network in distributed high-voltage transmission network, these sub-networks are widely distributed in different regions, entirely
Network model G is the set of these sub-networks, with G=(G1,G2,…,Gi,…,Gp) indicate, p is the quantity of sub-network.At this
In example, using Central China Power Grid as a sub-network of national power transmission network, the graph model of the sub-network is built, as shown in Figure 2.
(2) distributed power transmission network power frequency electromagnet field model, including power transmission network model, the physics of air and water body are established
Model.
(2.1) size in emulation domain is determined.Simulating area is the practical calculating area in the emulation of power transmission network power frequency electromagnetic field
Domain, it illustrates the space environment that electromagnetic radiation source is influenced.It includes two class physical models to emulate domain mainly, and one kind is air mode
Type, another kind of is water model:
Air model is:Mair=(Lair,Wair,Hair);
Ocean-land model is:Mocean-island=(Locean-island,Wocean-island,Hocean-island);
Wherein, MairAnd Mocean-islandIndicate the model of air and ocean-land, (Lair,Wair,Hair) and
(Locean-island,Wocean-island,Hocean-island) respectively indicate Air model and ocean-land model length.At this
In example, geometric parameter is as follows:
Air layer size is:
Lair=1500km
Wair=1500km
Hair=200km
Ocean-land layer size be:
Locean-island=1500km
Wocean-island=1500km
Hocean-island=5km
(2.2) physical attribute in emulation domain is determined.Different substances has different physical attributes, these physical attributes are determined
It is inconsistent to determine the power frequency electromagnet field distribution that the same radiation source is generated in different zones, it is therefore desirable to the difference of simulation model
Different physical attributes is arranged in material.
Air properties, (mur=1, epsilon=1, sigma=0);
Ocean attribute, (mur=1, epsilon=81, sigma=3);
Land attribute, (mur=1, epsilon=11.7, sigma=0);
Wherein mur indicates that the relative permeability of material, epsilon indicate that the relative dielectric constant of material, sigma indicate material
The conductivity of material.
(2.3) emulation domain geometrical model is established.In COMSOL softwares, a new component is established.In the geometry of component
Two cuboids are created in one column.One cuboid geometric parameter is (Locean-island,Wocean-island,Hocean-island), it represents
Ocean-land layer, ocean and land is two parts come cutting by coastline;Another cuboid geometric parameter is (Lair,
Wair,Hair), represent air layer.In one column of material of component corresponding physics is assigned to ocean-land layer and air layer respectively
Attribute.Since the physical property of different height air is not quite identical, upper atmosphere density is relatively low, and lower atmosphere layer density is higher,
So needing to air be layered that its physical parameter is arranged.The geometrical model in structure emulation domain, including air layer model and sea
Ocean-land layer model, as shown in Figure 3 and Figure 4.
(3) mesh generation is carried out to emulation domain, boundary condition setting carries out simulation calculation.Mesh generation is that finite element is imitative
Important step in very needs that size of mesh opening size, mesh generation pattern is reasonably arranged.In this example, mesh generation mould
Formula is set as " extreme refinement modes ", uses free tetrahedral grid, COMSOL softwares can be according to entire geometrical model
Size automatically configures the size of grid.As shown in figure 5, carrying out grid to ocean-two geometrical models of land layer and air layer
It divides.
Using the power frequency electromagnetic field in Maxwell solving equations space, Maxwell equation groups are asked using BiCGStab
It solves device to solve, setting solves control parameter, and default setting is kept to initialize solver.After the completion of initialization, setting changes
The step number that generation calculates is 10000 steps, and simulation frequency is set as 50Hz, calculates solver, and the work in space is calculated
Frequency electromagnetic field is distributed.
Based on Maxwell equation groups, the distributed power transmission network power frequency electromagnet field computation mould under boundary condition constraint is established
Type calculates space power frequency electromagnet field distribution of the distributed power transmission network in entirely emulation domain by finite element simulation.
(4) power frequency electromagnet field analysis.
Result of calculation is preserved, result of calculation is post-processed.It is shown using the data visualization module of COMSOL softwares
Show that result of calculation, specific analytical procedure are as follows:
(4.1) power frequency electromagnetic field Visualization of Simulation Results.Power frequency electromagnetic field is a vector field, and visualization can pass through
Motion pattern is shown.Meanwhile the power frequency electromagnetic field amplitude of some section is chosen, the power frequency electromagnet field strength of the plane can be drawn
Degree figure.The power frequency electromagnet field distribution of a certain height can also be shown using isogram.It is as follows:
(4.1.1) power frequency electromagnetic field motion pattern is drawn.Electric field and magnetic field that emulation obtains are the four-dimensional letters of room and time
Number can use E (x, y, z, t) and H (x, y, z, t) to indicate.E (x, y, z, t) and H (x, y, z, t) is three-dimensional vector field, is made
It is visualized with motion pattern.In this example, larger due to entirely emulating domain, it selects to take a point every 1km when exporting data,
The emulation data of electromagnetic field are exported, this is a space three-dimensional matrix.Then streamline figure layer is created, the number that streamline is arranged is
200, it selects data source for electric field and magnetic field, streamline is drawn according to the space coordinate of export data, each point on streamline
Electromagnetism field direction is the tangential direction of streamline.Fig. 6 and Fig. 7 is respectively magnetic field motion pattern and electric field motion pattern.
(4.1.2) power frequency electromagnet field intensity sectional drawing is drawn.Choose a certain height z0, export the electric field under the elevation plane
With magnetic field simulation data, electric field and magnetic field are respectively E (x, y, z0, t) and H (x, y, z0,t).It is calculate by the following formula electric field and magnetic field
Amplitude:
Wherein, | E (x, y, z0, t) | and | H (x, y, z0, t) | electric field magnitude and magnetic field amplitude, E are indicated respectivelyx,t、Ey,t、
Ez0,tThree components of the electric field under rectangular coordinate system, H are indicated respectivelyx,t、Hy,t、Hz0,tIndicate magnetic field in rectangular coordinate system respectively
Under three components.
In this example, in order to study power frequency electromagnetic field in the seawater with the propagation law in air, power frequency is studied first
Z is arranged in the aerial propagation law of electromagnetic field0For 10km, is i.e. at the above 10km in sea level, creates section and draw figure layer 1, choosing
It is electric field magnitude and magnetic field amplitude to select data source, draws the sectional drawing of power frequency electromagnetic field, adjusts diagram parameter repeatedly until obtaining
Ideal electromagnetic field sectional drawing.Fig. 8 and Fig. 9 is magnetic field sectional drawing and electric field sectional drawing;Then by z0Be set as -0.1km, i.e., it is extra large
Below plane at 100m, the propagation law of research power frequency electromagnetic field in the seawater.It creates section and draws figure layer 2, draw corresponding
Magnetic field sectional drawing and electric field sectional drawing, as shown in Figure 10 and Figure 11.
(4.1.3) power frequency electromagnetic field isogram is drawn.Such as step (4.1.2), the electric field width under a certain height is calculated
Value | E (x, y, z0, t) | and magnetic field amplitude | H (x, y, z0,t)|.Isoplethes drawing figure layer is created, selects data source for electric field magnitude
And magnetic field amplitude, the number that isopleth is arranged is 100, adjusts the type of isopleth, draws the isogram of power frequency electromagnetic field.Figure
12 and Figure 13 is 10km height magnetic field isogram and electric field isogram above sea level.Figure 14 and Figure 15 is below sea level
100m height magnetic field isogram and electric field isogram.
(4.2) calculating and its visualization of power frequency electromagnetic field divergence and curl.
Power frequency electromagnetic field is a vector field, and in MaxWell equation groups, the main means for describing electromagnetic field include electromagnetism
The divergence and curl of field, this patent is by the calculating to electromagnetic field divergence and curl, to portray distributed power transmission network generation
The feature of power frequency electromagnetic field.
Magnetic field divergence calculation formula is:
Magnetic field curl calculation formula is:
Electric field divergence calculation formula is:
Electric field curl calculation formula is:
The Simulation result data for exporting power frequency electromagnetic field, calculate separately the divergence and curl of electromagnetic field, wherein electromagnetic field
Divergence is Quantity Field, as shown in step 4.1.2, draws the sectional drawing of electromagnetic field divergence, as shown in Figure 16 and Figure 17.Electromagnetic field
Curl is vector field, as shown in step 4.1.1, draws the motion pattern of electromagnetic field curl, as shown in Figure 18 and Figure 19.
The above content as it will be easily appreciated by one skilled in the art that the foregoing is merely illustrative of the preferred embodiments of the present invention,
Be not intended to limit the invention, all within the spirits and principles of the present invention made by all any modification, equivalent and improvement etc.,
It should all be included in the protection scope of the present invention.
Claims (5)
1. a kind of distribution power transmission network power frequency electromagnetic field analysis method, which is characterized in that the method includes:
(1) according to the voltage and current on the node, loop and loop of distributed power transmission network, distributed power transmission network mould is built
Type;
(2) domain model is emulated according to the environment construction where distributed power transmission network;
(3) according to the distributed power transmission network work under distributed power transmission network model and emulation domain model structure boundary condition constraint
It is empty to go out power frequency electromagnetic field of the distributed power transmission network in entirely emulation domain by finite element stimulation for frequency electromagnetic field computation model
Between be distributed;
(4) by the form of power frequency electromagnetic field spatial distribution drafting pattern, visual figure is analyzed.
2. a kind of distributed power transmission network power frequency electromagnetic field analysis method according to claim 1, which is characterized in that described
Step (1) specifically includes:
(11) nodal analysis method of distributed power transmission network is built
V=(V1,V2,…,Vi,…,Vm),
Wherein, ViIndicate that a node of power transmission network, m indicate the number of nodes in network;
(12) loop model of distributed power transmission network is built
P=(P1,P2,…,Pi,…,Pk),
Wherein, PiIndicate that a loop in power transmission network, k indicate the loop quantity in power transmission network;
(13) the loop voltage current model of distributed power transmission network is built
E=(e1,e2,…,ei,…,ek),
I=(i1,i2,…,ii,…,ik),
Wherein, eiIndicate loop PiOn voltage, iiIndicate loop PiOn electric current;
(14) distributed power transmission network model is built
G=(G1,G2,…,Gi,…,Gp),
Wherein, Gi=(V, P, e, i), GiIndicate i-th of sub-network among power transmission network.
3. a kind of distributed power transmission network power frequency electromagnetic field analysis method according to claim 1, which is characterized in that described
Step (2) specifically includes:
(21) structure emulation domain geometrical model, geometrical model include:
Air model, Mair=(Lair,Wair,Hair);
Extra large land model, Mocean-island=(Locean-island,Wocean-island,Hocean-island);
Wherein, (Lair,Wair,Hair) and (Locean-island,Wocean-island,Hocean-island) Air model and extra large land are indicated respectively
The length and width and height of model;
(22) attribute model in structure emulation domain, attribute model include:
Air properties, (mur=1, epsilon=1, sigma=0);
Ocean attribute, (mur=1, epsilon=81, sigma=3);
Land attribute, (mur=1, epsilon=11.7, sigma=0);
Wherein, mur indicates that the relative permeability of material, epsilon indicate that the relative dielectric constant of material, sigma indicate material
Relative conductivity.
4. a kind of distributed power transmission network power frequency electromagnetic field analysis method according to claim 1, which is characterized in that described
Step (3) is specially:
Using distributed power transmission network model as the electromagnetic excitation source of simulation calculation, space of the emulation domain model as simulation calculation
Boundary constraint;Mesh generation is carried out to emulation domain, the entire space that calculates is divided into multiple sub-grids, passes through finite element simulation meter
The power frequency electromagnetic field of operator grid, the power frequency electromagnetic field of all-ones subnet lattice constitute the power frequency electromagnet field distribution in entire emulation domain.
5. a kind of distributed power transmission network power frequency electromagnetic field analysis method according to claim 1, which is characterized in that described
Step (4) specifically includes:
(41) create streamline figure layer, the number of streamline be set, by the electric field E (x, y, z, t) and magnetic field H of power frequency electromagnetic field (x, y,
Z, t) as data source drafting motion pattern;
(42) setting height z is chosen0, calculate the electric field magnitude and magnetic field amplitude of power frequency electromagnetic field under the height:
Wherein, Ex,t、Ey,t、Three components of the electric field under rectangular coordinate system, H are indicated respectivelyx,t、Hy,t、It indicates respectively
Three components of the magnetic field under rectangular coordinate system;It creates section and draws figure layer, by the electric field magnitude of power frequency electromagnetic field and magnetic field width
Value draws sectional drawing as data source;
(43) setting height z is chosen0, the electric field magnitude and magnetic field amplitude of power frequency electromagnetic field under the height are calculated, isopleth is created
Figure layer is drawn, the number of isopleth is set, the type of isopleth is adjusted, the electric field magnitude of power frequency electromagnetic field and magnetic field amplitude are made
For data source drawing isoline figure;
(44) power frequency electromagnetic field divergence and curl are calculated:
Magnetic field divergence,
Magnetic field curl,
Wherein, ▽ indicates divergence operator;▽ × expression curl operator;Indicate the unit vector of x-axis direction;Indicate y-axis side
To unit vector;Indicate the unit vector in z-axis direction;(x, y, z, t) was indicated between the three-dimensional coordinate added-time of point;
It creates section and draws figure layer, using the electromagnetic field divergence of power frequency electromagnetic field as data source, draw sectional drawing;Create motion pattern
Layer, is arranged the number of streamline, using the electromagnetic field curl of power frequency electromagnetic field as data source, draws motion pattern.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711464089.7A CN108287268B (en) | 2017-12-28 | 2017-12-28 | A kind of distribution power transmission network power frequency electromagnetic field analysis method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711464089.7A CN108287268B (en) | 2017-12-28 | 2017-12-28 | A kind of distribution power transmission network power frequency electromagnetic field analysis method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108287268A true CN108287268A (en) | 2018-07-17 |
CN108287268B CN108287268B (en) | 2019-07-23 |
Family
ID=62832536
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711464089.7A Active CN108287268B (en) | 2017-12-28 | 2017-12-28 | A kind of distribution power transmission network power frequency electromagnetic field analysis method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108287268B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109061323A (en) * | 2018-07-23 | 2018-12-21 | 电子科技大学 | A kind of near field antenna measurements method using spherical surface amplitude scan |
CN109991670A (en) * | 2018-12-29 | 2019-07-09 | 华中科技大学 | A kind of submarine target localization method |
CN111142163A (en) * | 2019-12-30 | 2020-05-12 | 华中科技大学 | Detection system and detection method for underwater sunken ship |
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 |
CN112084294A (en) * | 2020-09-14 | 2020-12-15 | 重庆长安新能源汽车科技有限公司 | Whole vehicle electromagnetic compatibility grading management method based on artificial intelligence |
CN112986732A (en) * | 2021-02-08 | 2021-06-18 | 吉林大学青岛汽车研究院 | Method and device for calculating electromagnetic radiation of whole vehicle to human body and storage medium |
CN117452081A (en) * | 2023-12-26 | 2024-01-26 | 国网天津市电力公司营销服务中心 | Electromagnetic interference calculation method and device, storage medium and electronic terminal |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101696997A (en) * | 2009-10-26 | 2010-04-21 | 重庆大学 | Adaptive rapid determining method for multi-medium power frequency electric field distribution |
CN102721875A (en) * | 2012-06-15 | 2012-10-10 | 四川电力科学研究院 | Method for measuring three-dimensional power-frequency electric field on basis of curve integral of equivalent charges on surfaces of transmission conductors |
CN103258124A (en) * | 2013-04-11 | 2013-08-21 | 东北电力大学 | Power frequency magnetic field assessment method for electromagnetic equipment of high voltage transformer substation |
CN103605836A (en) * | 2013-11-02 | 2014-02-26 | 国家电网公司 | Parallel computing method for three-dimensional electromagnetic fields of high-voltage transformer substation |
CN103760427A (en) * | 2014-01-15 | 2014-04-30 | 国家电网公司 | Statistical method of power frequency electromagnetic field distribution rules in transformer substation |
CN104090173A (en) * | 2014-07-25 | 2014-10-08 | 武汉三江中电科技有限责任公司 | Multi-node distributed field intensity test system and method based on Bluetooth communication |
CN106570229A (en) * | 2016-10-20 | 2017-04-19 | 中国电力科学研究院 | Three-dimensional power frequency electric field analysis method and system of AC filter filed of converter station |
-
2017
- 2017-12-28 CN CN201711464089.7A patent/CN108287268B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101696997A (en) * | 2009-10-26 | 2010-04-21 | 重庆大学 | Adaptive rapid determining method for multi-medium power frequency electric field distribution |
CN102721875A (en) * | 2012-06-15 | 2012-10-10 | 四川电力科学研究院 | Method for measuring three-dimensional power-frequency electric field on basis of curve integral of equivalent charges on surfaces of transmission conductors |
CN103258124A (en) * | 2013-04-11 | 2013-08-21 | 东北电力大学 | Power frequency magnetic field assessment method for electromagnetic equipment of high voltage transformer substation |
CN103605836A (en) * | 2013-11-02 | 2014-02-26 | 国家电网公司 | Parallel computing method for three-dimensional electromagnetic fields of high-voltage transformer substation |
CN103760427A (en) * | 2014-01-15 | 2014-04-30 | 国家电网公司 | Statistical method of power frequency electromagnetic field distribution rules in transformer substation |
CN104090173A (en) * | 2014-07-25 | 2014-10-08 | 武汉三江中电科技有限责任公司 | Multi-node distributed field intensity test system and method based on Bluetooth communication |
CN106570229A (en) * | 2016-10-20 | 2017-04-19 | 中国电力科学研究院 | Three-dimensional power frequency electric field analysis method and system of AC filter filed of converter station |
Non-Patent Citations (3)
Title |
---|
仇丰: "不同环境条件下工频电磁污染的空间分布研究", 《中国优秀硕士学位论文全文数据库》 * |
张可心: "超高压传输线电磁环境特性研究", 《中国优秀硕士学位论文全文数据库》 * |
牟秋谷等: "3种典型220kV_GIS变电站进线方式的电磁场分析", 《四川电力技术》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109061323A (en) * | 2018-07-23 | 2018-12-21 | 电子科技大学 | A kind of near field antenna measurements method using spherical surface amplitude scan |
CN109061323B (en) * | 2018-07-23 | 2020-12-29 | 电子科技大学 | Near-field antenna measurement method adopting spherical amplitude scanning |
CN109991670A (en) * | 2018-12-29 | 2019-07-09 | 华中科技大学 | A kind of submarine target localization method |
CN111142163A (en) * | 2019-12-30 | 2020-05-12 | 华中科技大学 | Detection system and detection method for underwater sunken ship |
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 |
CN112084294A (en) * | 2020-09-14 | 2020-12-15 | 重庆长安新能源汽车科技有限公司 | Whole vehicle electromagnetic compatibility grading management method based on artificial intelligence |
CN112084294B (en) * | 2020-09-14 | 2022-07-26 | 重庆长安新能源汽车科技有限公司 | Whole vehicle electromagnetic compatibility grading management method based on artificial intelligence |
CN112986732A (en) * | 2021-02-08 | 2021-06-18 | 吉林大学青岛汽车研究院 | Method and device for calculating electromagnetic radiation of whole vehicle to human body and storage medium |
CN112986732B (en) * | 2021-02-08 | 2022-06-14 | 吉林大学青岛汽车研究院 | Method and device for calculating electromagnetic radiation of whole vehicle to human body and storage medium |
CN117452081A (en) * | 2023-12-26 | 2024-01-26 | 国网天津市电力公司营销服务中心 | Electromagnetic interference calculation method and device, storage medium and electronic terminal |
Also Published As
Publication number | Publication date |
---|---|
CN108287268B (en) | 2019-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108287268B (en) | A kind of distribution power transmission network power frequency electromagnetic field analysis method | |
CN103605838B (en) | A kind of power distribution network uniform wiring map generalization method being applicable to Study on Trend | |
CN103559360B (en) | The construction method of a kind of three-dimensional overhead electrical power line model and system | |
CN103955566A (en) | Electric transmission line three-dimensional design method based on Google Earth and electric transmission line three-dimensional design system based on Google Earth | |
Duan et al. | Efficient implementation for 3-D Laguerre-based finite-difference time-domain method | |
CN103810753A (en) | Method for visually showing power grid dispatching businesses | |
CN103093030A (en) | Three-dimensional (3D) electric power line model building method and building device | |
CN107231250A (en) | One kind is based on electric network information physical system perception data compressive sampling method and device | |
CN102708229A (en) | Matrix decomposition and novel singular value decomposition combined method for complex layered medium structures | |
CN103699960A (en) | Medium voltage distribution power grid operation state visualization method | |
CN102141587A (en) | Lightening locating data imaging method | |
CN105260583A (en) | Method and system for calculating biological effect of extra-high-voltage power frequency electromagnetic field on human body | |
CN104793065A (en) | Method of determining radio interference in AC transmission line crossed erection region | |
Melo et al. | Data issues in spatial electric load forecasting | |
CN103150412B (en) | Based on the electric power facility site selecting method of Spatial Data Mining Technique | |
CN107168994A (en) | Design implementation method based on power network GIS platform map denotation style model storehouse | |
CN105468759B (en) | The frequency spectrum data construction method of space body | |
Usui et al. | A multi-scale electromagnetic particle code with adaptive mesh refinement and its parallelization | |
CN101763659A (en) | Method for automatically generating large-scale digitization altitude data model based on image set | |
CN111209663B (en) | Modeling analysis method and system for equivalent ultra-long wave antenna array of high-voltage power grid | |
CN109166307B (en) | Visual metaphor expression method and system for time-space accessibility of traffic network | |
Wert et al. | Contouring method considerations for power systems applications | |
CN108304651A (en) | The third boundary condition processing method of dc resistivity Finite element method simulation | |
Si et al. | DC grounding electrode potential based on a kriging geoelectric structure model | |
Plutzar et al. | A GIS-based Approach for Modeling the Spatial and Temporal Development of Night-time lights |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |