CN101696997B - Adaptive rapid determining method for multi-medium power frequency electric field distribution - Google Patents

Abstract

The invention relates to an adaptive rapid determining method for multi-medium power frequency electric field distribution, which comprises the following steps: 1, acquisition of parameters, namely, acquiring a geometric parameter, a geometric size and a space coordinate position of each medium arranged close to a working corridor in a transformer substation, and measuring relative dielectric constants of each medium; 2, determination of a power frequency electric field, namely, establishing and dividing a power frequency electric field physical model, constructing a data storage structure, determining the electric field value of each node in the power frequency electric field physical model and the like; and 3, drawing of a space electric field contour distribution map. The invention can judge whether the power frequency electric field value in the area which is 1.8 meters away from the ground or below exceeds a national regulated limit value 4kV/m; compared with a conventional method, the adaptive rapid determining method has the advantages of reducing memory demand by 70 percent, shortening calculating time by 50 percent and improving calculating precision by 5 percent. The adaptive rapid determining method for the multi-medium power frequency electric field distribution is applied to rapid determination and estimation of an on-site power frequency electric field, and a safety design of power transmission and transformation.

Description

The self-adaptation that a kind of multi-medium power frequency electric field distributes is confirmed method fast

Technical field

The present invention is directed to such as transformer station in this under multiple medium concurrent conditions, the bottleneck problem of existing definite method that power frequency electric field distributes has been invented the self-adaptation that a kind of multi-medium power frequency electric field distributes and has been confirmed method fast.

Background technology

Along with the continuous development of power industry and the raising of transmission line of electricity electric pressure, the power frequency high voltage electric field has caused people's extensive concern to the infringement of health.Also there is not at present a kind of effective method to be used for confirming fast that the power frequency electric field such as under the multimedium coexistence in the transformer station distributes.Because hypermedia existence in the transformer station, cause the distribution of transformer station's power frequency electric field to be confirmed being one in infinite region by the coefficient challenge of multiple medium.That need disperse when adopting such problem of Treatment of Finite Element Method finds the solution the zone for infinitely great, therefore uses this method and finds the solution being that efficient is extremely low.Though and boundary element is found the solution the infinite region problem and need not discrete whole infinite region and only need the border of dispersed problem and the interface of different medium, the discrete boundary element coefficient of potential matrix of generation afterwards is a non-singular matrix, and its calculated amount and memory space all are 0 (N ²) order of magnitude, said N is the node number of multi-medium power frequency electric field problem, adopting said method confirms that multi-medium power frequency electric field distribution thing efficient is low equally.Boundary Element method based on the Fast Multiple Method method can reduce calculated amount and memory space significantly, thereby improves counting yield.This method has been applied to other field, solves three-dimensional static and dynamic slight crack problem as quickening boundary element, shortens and finds the solution the time; Be mingled with problem in analog computation compound substance rigidity, can find the solution the up to ten million problem of unknown quantity in several hours, have very high operation efficiency and assess the cost, but this method is confirming that the application that multi-medium power frequency electric field distributes does not also have relevant report at present with low.

Summary of the invention

In view of quickening to find the solution the boundary element problem, the Fast Multipole submethod obtains the applications well effect at other field.The present invention adopts the Fast Multipole submethod to combine the boundary element method self-adaptation to confirm that fast multi-medium power frequency electric field distributes; Solution is by the bottleneck of operand and counting yield that memory space causes; Realize the high efficiency transformer station's multi-medium power frequency electric field distribution of confirming fast; The interference that pollution source that quick understanding power frequency electric field exceeds standard and power frequency high voltage electric field exceed standard and bring, and be used for the engineering design needs for the power transmission and transformation deviser improves effective instrument.

Therefore, the purpose of this invention is to provide a kind of power frequency electric field distribution self-adaptation and confirm method fast.Physical dimension through obtaining various media in the transformer station and volume coordinate position and with corresponding relative dielectric constant, estimate the interference that pollution source that power frequency electric field exceeds standard and power frequency high voltage electric field exceed standard and bring fast.

Said medium comprises that transformer station contains the equipment that can produce power frequency electric field of disconnecting link, isolating switch, voltage transformer (VT), bus or lightning arrester.

For realizing that the technical scheme that above-mentioned purpose adopts is such, promptly a kind of power frequency electric field distribution self-adaptation is confirmed method fast, and method comprises the steps:

1, parameter obtains

1), adopts infrared range-measurement system to obtain to be arranged on the physical dimension and the volume coordinate position of near the various media in work corridor transformer station in serial communication function; The volume coordinate position of the physical dimension of said various media and said relatively work corridor starting point comprises:

The physical dimension and the volume coordinate position that contain the equipment of disconnecting link are the level and the vertical range of equipment overall height, pillar height, insulation height, a column radius, insulation radius, horizontal knife up degree, insulating trip number, disconnecting link spacing and disconnecting link pillar centre distance transformer station work corridor starting point;

The physical dimension and the volume coordinate position that contain the equipment of isolating switch are the level and the vertical range of equipment overall height, pillar height, insulation height, a column radius, insulation radius, insulating trip number and isolating switch pillar centre distance transformer station work corridor starting point;

The physical dimension and the volume coordinate position that contain the equipment of voltage transformer (VT) are the level and the vertical range of equipment overall height, pillar height, insulation height, a column radius, insulation radius, insulating trip number and voltage transformer (VT) pillar centre distance transformer station work corridor starting point;

The physical dimension and the volume coordinate position that contain the equipment of bus are the level and the vertical range of equipment overall height, pillar height, insulation height, a column radius, insulation radius and bus pillar centre distance transformer station work corridor starting point;

The physical dimension and the volume coordinate position that contain the equipment of lightning arrester are the level and the vertical range of equipment overall height, pillar height, insulation height, a column radius, insulation radius and lightning arrester pillar centre distance transformer station work corridor starting point;

2), adopt decameter to measure the relative dielectric constant of various media with serial communication function; Said relative dielectric constant is respectively: the ratio of the specific inductive capacity of the specific inductive capacity of insulated part, leg portion and environment division and air in the equipment of disconnecting link, isolating switch, voltage transformer (VT), bus or lightning arrester, said environment division are disconnecting link, isolating switch, voltage transformer (VT), bus or lightning arrester itself;

2, power frequency electric field confirms

1), the foundation of power frequency electric field physical model and subdivision: the relative dielectric constant of the medium physical dimension that collects and volume coordinate position and various media is made up physical model; Promptly adopt the commercially available prod Hypermesh software of U.S. Altair company to draw leg portion, insulated part and the environment division three-dimensional picture of wrapping up medium, comprising:

The leg portion that contains the equipment of disconnecting link: level and the vertical range with its pillar height, a column radius and disconnecting link pillar centre distance transformer station work corridor starting point is the center of cylinder height, bottom surface radius and bottom surface circle respectively, the draw circles cylinder; Insulated part: level, vertical range, pillar height and the insulating trip number with insulation height, insulation radius, disconnecting link pillar centre distance transformer station work corridor is that the insulator chain right cylinder is drawn at cylinder height and bottom surface radius, center respectively; Environment division: level, vertical range and the horizontal knife up degree with horizontal knife up degree, a column radius and disconnecting link pillar centre distance transformer station work corridor is the center draw circles cylinder that justify cylinder height and bottom surface radius, bottom surface respectively;

The equipment leg portion that contains isolating switch: level and the vertical range with pillar height, a column radius and isolating switch pillar centre distance transformer station work corridor starting point is cylinder height, bottom surface radius, circle center, bottom surface draw circles cylinder respectively; Insulated part: level, vertical range, pillar height and the insulating trip number with insulation height, isolating switch pillar centre distance transformer station work corridor starting point is that the insulator chain right cylinder is drawn at cylinder height and bottom surface radius, center respectively; Environment division: level and the vertical range with equipment height overall, a column radius and isolating switch pillar centre distance transformer station work corridor starting point is cylinder height, bottom surface radius, circle center, bottom surface draw circles cylinder respectively;

The equipment leg portion that contains voltage transformer (VT): level and the vertical range with pillar height, a column radius and voltage transformer (VT) pillar centre distance transformer station work corridor starting point is the center draw circles cylinder of cylinder height, bottom surface radius, bottom surface circle respectively; Insulated part: level, vertical range, pillar height and insulating trip number with insulation height, insulation radius, voltage transformer (VT) pillar centre distance transformer station work corridor starting point are that the insulator chain right cylinder is drawn at cylinder height and bottom surface radius, round center, bottom surface; Environment division: level and the vertical range with equipment height overall, a column radius and voltage transformer (VT) pillar centre distance transformer station work corridor starting point is cylinder height, bottom surface radius, circle center, bottom surface draw circles cylinder respectively;

The equipment leg portion that contains bus: level and the vertical range with pillar height, a column radius and bus pillar centre distance transformer station work corridor starting point is cylinder height, bottom surface radius, center draw circles cylinder respectively; Insulated part: level, vertical range, pillar height with insulation height, insulation radius, bus pillar centre distance transformer station work corridor starting point are that the insulator chain right cylinder is drawn at cylinder height and bottom surface radius, round center, bottom surface; Environment division: level and the vertical range with equipment height overall, a column radius and bus pillar centre distance transformer station work corridor starting point is cylinder height, bottom surface radius, circle center, bottom surface draw circles cylinder respectively;

The leg portion that contains the equipment of lightning arrester: level and the vertical range with pillar height, a column radius and lightning arrester pillar centre distance transformer station work corridor starting point is the center draw circles cylinder of cylinder height, bottom surface radius, bottom surface circle respectively; Insulated part: level, vertical range, pillar height and insulating trip number with insulation height, insulation height, lightning arrester pillar centre distance transformer station work corridor starting point are that the insulator chain right cylinder is drawn at cylinder height and bottom surface radius, round center, bottom surface; Environment division: level and the vertical range with equipment height overall, a column radius and lightning arrester pillar centre distance transformer station work corridor starting point is cylinder height, bottom surface radius, circle center, bottom surface draw circles cylinder respectively.

Automatically the gained three-dimensional picture is carried out the self-adaptation subdivision through this software then; The node that obtains constituting physical model is counted L and unit number M; Said L represent this model use L three dimensions point All Media is described; Said M represent this model use M three dimensions plane All Media is described; Use array A and B memory node and unit essential information respectively; The numbering of the corresponding node of each row among the array A and the three-dimensional coordinate figure of this aforementioned relatively work corridor starting point, the numbering of each the corresponding unit of row among the array B constitutes the whole zone of finding the solution with the information that the numbering of aforementioned 3 different nodes is described all unit of segment space information of certain medium;

2), make up data store organisation: utilize aforementioned nodes essential information array A and unit essential information array B to make up octree structure; The aforementioned whole zone of finding the solution is decomposed; Each point is the unit information of the medium that comprises of a square in the said octree structure, and the number of the square of decomposition is N; The multipole momentum moment of each square representes that with array E the size of E is N * 10 * 10, and wherein N is the number of the square of decomposition.The octree structure that forms is as shown in Figure 3, and the leafy node of the Octree among said Fig. 3 is the minimum square that can not cut apart again in the octree structure.The detailed storage means of its data store organisation is referring to [Wang Wu, the application of tree construction in N body problem, computer utility research, 2008/01].In this article, the relation between per two cubes in the tree construction is carried out following definition:

Close relative's cube be in the k layer with the cube of the more than intersection point of C cube.

Distant relative's cube be in the k layer with the cube of C cube neither one intersection point.

It is to be close relative's cube at k layer and C that the close relative does not have subclass, and is the cube set of leaf.

The mutual collection of distant relative is at k layer and the equirotal distant relative's cube of C, in the k-1 layer with C under cube be that close relative's cube is gathered.

It is to be close relative's cube at k layer and C that the sublayer does not have subclass, at the sub-cube greater than k layer and C.

Leaf cube set for the distant relative.

Father's layer does not have distant relative's cube that subclass is k layer and the different size of C, but belongs to the set of close relative's cube at the cube at k-1 layer and C place.

3), the confirming of each node electric field value in the power frequency electric field physical model

Given each node electric field value initial vector, and be kept among the array D, the size of D is L;

Each node electric field value confirms that belonging to cubical close relative by node contributes G in the power frequency electric field physical model _JiContribute F with the distant relative _JlmBe formed by stacking.

Certain node electric field value confirms that belonging to cubical close relative by this node contributes G in the power frequency electric field physical model _JiContribute F with the distant relative _JlmBe formed by stacking.

Each cubical close relative contributes G _JiConfirm through (1) formula

$G_{\mathrm{ji}}={\∫}_{S_i}\mathrm{\σ}\left(y\right)G(x,y)\mathrm{dS}---\left(1\right)$

In the formula: s _iBe i unit, (x y) is three-dimensional Green function to G.

And the distant relative contributes F _JlmConfirm through following steps:

Polymerization successively.Polymerization successively is meant from the leafy node of octree structure and begins, upwards polymerization successively.Nodal information array A, unit information array B and node electric field intensity initial value array D according to above formation; Confirm the multipole momentum moment of this unit for any leaf j in the octree structure by (2) formula at aforesaid leafy node center; And the multipole momentum moment of other unit in the leafy node in leafy node added up, form in this leafy node all unit at the multipole momentum moment three-dimensional array E at leaf center.

$E_{\mathrm{jlm}}\left(y_0\right)=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}\underset{S_i}{\∫}{R}_{\mathrm{lm}}(y,y_0)\mathrm{\σ}\left(y\right)\mathrm{dS}---\left(2\right)$

Y wherein ₀Be the centre coordinate of square, W is the number of unit in the square, s _iBe i unit, E _Jlm(y ₀) be j cube at the center y ₀The lm rank multipole momentum moment, the span of l and m is respectively 0～9;

Leafy node from aforesaid Octree begins then, each cubical multipole momentum moment in the aforesaid Octree is delivered to the multipole momentum moment E at its father node center according to (3) formula _Jlm, till the multipole momentum moment of the second layer node of octree structure has been confirmed;

$E_{\mathrm{jlm}}\left(y^{\′}\right)=\underset{l^{\′}=0}{\overset{l}{\mathrm{\Σ}}}\underset{m^{\′}=-l^{\′}}{\overset{l^{\′}}{\mathrm{\Σ}}}(-1)R_{l^{\′}{m}^{\′}}(y_0,y^{\′})E_{\left(j\right)(l-l^{\′})(m-m^{\′})}\left(y_0\right)---\left(3\right)$

Wherein y ' is a father node cube centre coordinate.

Successively shift and disperse: with each cubical transfer and coefficient of divergence F in the array F storage octree structure _Jlm, its size is N * 10 * 10.

Utilize the data among the multipole momentum moment three-dimensional array E, in the octree structure greater than two-layer cube, confirm that through (4) formula its distant relative collects coefficient of migration F alternately _Jlm

$F_{\mathrm{jlm}}\left(x^{\′}\right)=\frac{{(-i)}^{\left|m\right|}}{\sqrt{(l-m)!(l+m)!}}\underset{l^{\′}=1}{\overset{l}{\mathrm{\Σ}}}\underset{m^{\′}=1}{\overset{m}{\mathrm{\Σ}}}{V}_{l^{\′},m^{\′}}\left(x^{\′}\right)e^{\mathrm{im}}---\left(4\right)$

Wherein x ' is for shifting cubical center.

Begin from second layer cube then, through (5) formula each cubical distant relative of aforesaid Octree is collected coefficient of migration alternately and diffuse to its child node cube center, till the bottom of tree construction, form each cubical coefficient of divergence F _Jlm

$F_{\mathrm{jlm}}\left(x_0\right)=\underset{l^{\′}=0}{\overset{\∞}{\mathrm{\Σ}}}\underset{m^{\′}=-l^{\′}}{\overset{l^{\′}}{\mathrm{\Σ}}}{F}_{{\mathrm{jl}}^{\′}{m}^{\′}}\left(x^{\′}\right)R_{(l-l^{\′})(m-m^{\′})}(x_0,x^{\′})---\left(5\right)$

X wherein ₀The leafy node cube center of representing aforesaid Octree.

Like this, the node place cube close relative who is confirmed by the front contributes G _JiContribute F with the distant relative _JlmStack forms each node electric field value vector X in the physical model, shown in (6) formula.

$X=\underset{q}{\mathrm{\Σ}}\underset{i\∈W_q}{\mathrm{\Σ}}{q}_i{G}_{\mathrm{ji}}+$ (6)

$\underset{q}{\mathrm{\&Sigma;}}\underset{i\&Element;{\mathrm{<figure-callout\; id="1"\; label="W\; q"\; filenames="H2009101912155E00021.png"\; state="\{\{state\}\}">W</figure-callout>}}_q}{\mathrm{\&Sigma;}}\frac{1}{4\mathrm{\&pi;}}\underset{l=0}{\overset{\&infin;}{\mathrm{\&Sigma;}}}\underset{m=-l}{\overset{l}{\mathrm{\&Sigma;}}}{R}_{\mathrm{lm}}(x,x_0)\&times;F_{\mathrm{jlm}}\left(x_0\right)$

4), through the difference of the vectorial X of each node electric field value in each node electric field value initial vector D and the power frequency electric field physical model, confirm the norm of difference, if less than 10 ^-9, stop to calculate; Otherwise repeating step 3), up to the norm of difference less than 10 ^-9Till.The norm of said difference was confirmed referring to [Lin Chengsen, numerical analysis, Science Press, in January, 2006].

3, rendering space electric field equipotential line distribution plan

According to step 2 4) in satisfy to confirm each node electric field value vector X and step 2 in the physical model of precision 1) one-to-one relationship of middle nodal information array A; Draw out by work corridor starting point and begin the electric field equipotential line distribution plan on the perpendicular of work corridor, edge direction; Whether judgement surpasses the limit value 4KV/m of national regulation apart from 1.8 meters on ground and with lower area power frequency electric field value.The power frequency electric field that is the work corridor above limit value 4KV/m exceeds standard.

The present invention can be applied to:

1, confirms the power frequency electric field in 110KV and above transformer station work corridor;

2, confirm the power frequency electric field of generating plant part key equipment surrounding space;

3, confirm the power frequency electric field of 110KV and above switch yard;

The present invention can reach following performance index:

1, whether decidable surpasses the limit value 4KV/m of national regulation apart from 1.8 meters on ground and with lower area power frequency electric field value;

2, can form space line-frequency electric field equipotential line distribution plan;

3, the physical model unit number hanged down at least 70% than conventional method greater than 5000 o'clock in memory requirements;

4, the physical model unit number shortened 50% than conventional method greater than 5000 o'clock on computing time;

5, the physical model unit number can improve 5% with respect to conventional method greater than 5000 o'clock on computational accuracy.

Conspicuous during good effect that the present invention brings, promptly use the present invention, can confirm fast that multi-medium power frequency electric field distributes, and improve effective prediction and calculation instrument for environmental assessment and power transmission and transformation devisers.

Description of drawings

Fig. 1 octree structure of the present invention

The geometric model figure of the A of 500kV transformer station phase equipment in Fig. 2 embodiment of the invention

The electric field equipotential line distribution plan on the plane, work corridor that Fig. 3 the present invention forms

The variation relation figure of memory space and unit number in Fig. 4 the inventive method

The variation relation figure of computing time and unit number in Fig. 5 the inventive method

The variation relation figure of error precision and unit number in Fig. 6 the inventive method

Embodiment

Geometric model figure below in conjunction with the A of 500kV transformer station phase equipment; On the desk-top computer of cruel profit double-core 1.6GHz CPU and 1024MB RAM, calculate; Operation system of computer is Windows XP, and value storage adopts single precision, and the present invention is further described.

The collection of obtaining of parameter comprises the infrared range-measurement system of gathering various medium geometric parameters and having serial communication function, decameter, the PC that the various medium specific inductive capacity of measurement have serial communication function.Wherein have the infrared range-measurement system of serial communication function, decameter with serial communication function; Be the commercially available prod, the infrared range-measurement system with serial communication function adopts the new translocation in Changsha to paint the ND-3000/2000 product with serial communication function that Instr Ltd. produces; Decameter with serial communication function adopts the ZJ-3J product with serial communication function of Duozhu Science-Technology Development Co., Ltd., Nanjing's production.Method step is following

1, obtaining of parameter: through the medium geometric parameter such as the table 1 of ND-3000/2000 model infrared range-measurement system measurement

Table 1: the geometric parameters numerical table of equipment (unit: rice)

Exhausted exhausted cutter is apart from the worker

Transfer is from work

Equipment edge pillar edge post edge lock is done the corridor

Equipment

Cutter

The corridor starting point

Height overall high altitude half half intersegmental starting point water

High

Vertical range

Footpath, degree footpath number anomaly distance

Disconnecting link 9.6 4.5 5.1 0.1 0.1 43 5.3 5.3 1.3

Isolating switch 94 5.0 0.1 0.1/3/6.2 1.3

Voltage transformer (VT) 8.7 3.2 5.5 0.1 0.1/3/7.3 1.3

Bus 23.5 24 24 0.3 0.3 ///8.5 1.3

Lightning arrester 10 6.3 3.7 0.1 0.1 ///10.1 1.3

The relative dielectric constant such as the table 2 of the medium of measuring through ZJ-3J model decameter

Table 2: the relative dielectric constant table of equipment

Environment division insulated part leg portion

Disconnecting link 3.2 5.7 2.5

Isolating switch 1.8 7.4 3.2

Voltage transformer (VT) 3.1 4.2 1.5

Bus 1.51 5.9 1.9

Lightning arrester 2.1 6.3 2.8

2, power frequency electric field confirms

1, the foundation of power frequency electric field physical model and subdivision.Geometric parameter such as table 1 and medium relative dielectric constant such as table 2 through the equipment of collecting; Leg portion, insulated part and the environment division three-dimensional picture of being drawn the parcel medium by the commercially available prod Hypermesh software of U.S. Altair company are as shown in Figure 4; And carry out the self-adaptation subdivision with this software; Obtaining node number L is 29350, and the unit number M is 53250, and forms node and unit essential information array A and B.

$A=\left[\begin{array}{cccc}1& 0& -0.13& 0\\ 2& 0.13& 0& 0\\ .& .& .& .\\ .& .& .& .\\ .& .& .& .\\ 29350& 0.08& 5.05& 9.6\end{array}\right]$

$B=\left[\begin{array}{cccc}1& 49& 1& 1087\\ 2& 123& 5& 29\\ .& .& .& .\\ .& .& .& .\\ .& .& .& .\\ 53250& 29235& 6538& 198\end{array}\right]$

2), make up data store organisation.The octree structure that forms, obtaining the square number is 895.The size of the multipole momentum moment E of each square is 895 * 10 * 10.

3), the confirming of each node electric field value in the power frequency electric field physical model

And given each node electric field value initial value, just be made as complete 0, totally 29350, and be kept among the array D, promptly

D＝[0?0?…?0]

Each node electric field value confirms that belonging to cubical close relative by node contributes G in the power frequency electric field physical model _JiContribute F with the distant relative _JlmBe formed by stacking.Bring nodal information array A, unit information array B and tree construction array E into formula (1) respectively and confirm that the close relative contributes G _Ji, obtain following matrix.

$G=\left[\begin{array}{cccc}1& 5238.6& 2569.4& 3698.3\\ 2& 3659.7& 8934.7& 3574.1\\ .& .& .& .\\ .& .& .& .\\ .& .& .& .\\ .& .& .& .\\ .& .& .& .\\ .& .& .& .\end{array}\right]$

A, polymerization successively.Nodal information array A, unit information array B and node electric field intensity initial value array D difference substitution formula (2) and (3) are confirmed the multipole momentum moment E of its distant relative's contribution _Jlm, obtain following data.

$E=\left[\begin{array}{cccc}1& 365.3& ...& 689.5\\ 2& 539.7& ...& 432.0\\ .& .& .& .\\ .& .& .& .\\ .& .& .& .\\ 895& 986.3& ...& 863.1\end{array}\right]$

B, successively shift and disperse.Data substitution formula (4) among the multipole momentum moment three-dimensional array E and (5) are confirmed the coefficient of divergence F that its distant relative contributes _Jlm

$F=\left[\begin{array}{cccc}1& 452.6& ...& 1023.7\\ 2& 2332.9& ...& 2389.1\\ .& .& .& .\\ .& .& .& .\\ .& .& .& .\\ 895& 3569.8& ...& 3617.6\end{array}\right]$

The close relative is contributed G _JiCoefficient of divergence F with distant relative's contribution _JlmSubstitution relational expression (6), each node electric field value vector X does in the physical model of determining.

X＝[2363.9?3236.4?…?7869.3]

4), the norm through the difference of each node electric field value vector X in each node electric field value initial vector D and the physical model, discovery is through after 92 times, the norm of difference is less than 10 ^-9Stop to calculate, as shown in the table.

The iterations iteration error

1 1.30E-01

15 3.02E-02

24 5.30E-03

36 4.60E-04

44 9.50E-05

58 5.50E-06

72 3.11E-07

92 9.10E-10

3, rendering space electric field equipotential line distribution plan

According to the relation of node coordinate information array A in each node electric field value vector X and the physical model in the physical model, space electric field equipotential line distribution plan such as Fig. 5 of drawing out according to the mode of drawing a contour line.

Effect analysis

Through analyzing equipotential line at the area discover below 1.8 meters, this zone power frequency electric field value is the limit value 4KV/m that does not surpass national regulation.Zone among the figure below the red line is the zone below 1.8 meters, finds through observing this regional equipotential line value corresponding, and the equipotential line that this zone is passed mainly contains 1.2,1.8KV/m etc., all surpasses 4KV/m, satisfies regulation electromagnetic environment index.

The quick method of self-adaptation that the present invention proposes and Computer Storage amount comparing analysis such as Fig. 6 of direct method, as can be seen from the figure, the physical model unit number was greater than 5000 o'clock, and the memory demand of the quick method of self-adaptation is than conventional method low at least 70%.

The calculating T.T. and the relationship analysis between unit number of the quick method of self-adaptation are as shown in Figure 7.As can be seen from the figure, the physical model unit number shortened 50% than conventional method greater than 5000 o'clock on computing time; In conjunction with the calculator memory demand, the superiority of finding the solution many unknown quantitys of large scale problem with the quick method of self-adaptation is just more outstanding.

The calculating maximum relative error and the unit number relationship analysis of the quick method of self-adaptation are as shown in Figure 8.As can be seen from the figure, the physical model unit number was greater than 5000 o'clock, and is higher by 5% than conventional method on computational accuracy; This method is used for confirming that multi-medium power frequency electric field has enough precision.

Claims (1)

1. the self-adaptation of a multi-medium power frequency electric field distribution is confirmed method fast, and it is characterized in that: method may further comprise the steps:

1), parameter obtains

A), adopt infrared range-measurement system to obtain to be arranged on the physical dimension and the volume coordinate position of near the various media in work corridor transformer station in serial communication function; The volume coordinate position of the physical dimension of said various media and said relatively work corridor starting point comprises:

The physical dimension and the volume coordinate position that contain the equipment of disconnecting link are the level and the vertical range of equipment overall height, pillar height, insulation height, a column radius, insulation radius, horizontal knife up degree, insulating trip number, disconnecting link spacing and disconnecting link pillar centre distance transformer station work corridor starting point;

The physical dimension and the volume coordinate position that contain the equipment of isolating switch are the level and the vertical range of equipment overall height, pillar height, insulation height, a column radius, insulation radius, insulating trip number and isolating switch pillar centre distance transformer station work corridor starting point;

The physical dimension and the volume coordinate position that contain the equipment of voltage transformer (VT) are the level and the vertical range of equipment overall height, pillar height, insulation height, a column radius, insulation radius, insulating trip number and voltage transformer (VT) pillar centre distance transformer station work corridor starting point;

The physical dimension and the volume coordinate position that contain the equipment of bus are the level and the vertical range of equipment overall height, pillar height, insulation height, a column radius, insulation radius and bus pillar centre distance transformer station work corridor starting point;

The physical dimension and the volume coordinate position that contain the equipment of lightning arrester are the level and the vertical range of equipment overall height, pillar height, insulation height, a column radius, insulation radius and lightning arrester pillar centre distance transformer station work corridor starting point;

B), adopt decameter to measure the relative dielectric constant of various media with serial communication function; Said relative dielectric constant is respectively: the ratio of the specific inductive capacity of the specific inductive capacity of insulated part, leg portion and environment division and air in the equipment of disconnecting link, isolating switch, voltage transformer (VT), bus or lightning arrester, said environment division are disconnecting link, isolating switch, voltage transformer (VT), bus or lightning arrester itself;

2), power frequency electric field confirms

A), the foundation and the subdivision of power frequency electric field physical model: the relative dielectric constant of the medium physical dimension that collects and volume coordinate position and various media is made up physical model, promptly adopt the commercial software product Hypermesh of U.S. Altair company to draw leg portion, insulated part and the environment division three-dimensional picture of wrapping up said medium; Automatically the gained three-dimensional picture is carried out the self-adaptation subdivision by said software; Obtain constituting the node number (L) and the unit number (M) of physical model; Said node number (L) represent this model use L three dimensions point All Media is described; Said unit number (M) represent this model use M three dimensions plane All Media is described; Use array A and array B memory node and unit essential information respectively; The numbering of the corresponding node of each row among the array A and the three-dimensional coordinate figure of this aforementioned relatively work corridor starting point, the numbering of each the corresponding unit of row among the array B constitutes the whole zone of finding the solution with the information that the numbering of aforementioned 3 different nodes is described all unit of segment space information of certain medium;

B), make up data store organisation: utilize the basic array B in aforementioned nodes essential information array A and unit to make up octree structure; The aforementioned whole zone of finding the solution is decomposed; Each point is the unit information of the medium that comprises of a square in the said octree structure, and the number of the square of decomposition is N; The multipole momentum moment of each square representes that with array E the size of array E is N * 10 * 10, and wherein N is the number of the square of decomposition,

C), the confirming of each node electric field value in the power frequency electric field physical model

The electric field value initial vector of given each node, and be kept among the array D, the size of array D is L;

Certain node electric field value confirms that belonging to cubical close relative by this node contributes G in the power frequency electric field physical model _JiContribute F with the distant relative _JlmBe formed by stacking;

Each cubical close relative contributes G _JiConfirm through (1) formula

In the formula: s _iBe i unit, (x y) is three-dimensional Green function to G;

And the distant relative contributes F _JlmConfirm through following steps:

Polymerization successively: polymerization successively is meant from the leafy node of octree structure and begins; Upwards polymerization successively; Electric field value initial vector D according to nodal information array A, unit information array B and the node of above formation; Confirm the multipole momentum moment of this unit for any leaf j in the octree structure by (2) formula at aforesaid leafy node center; And the multipole momentum moment of other unit in the leafy node in leafy node added up, form in this leafy node all unit at the multipole momentum moment three-dimensional array E at leaf center;

Y wherein ₀Be the centre coordinate of square, N is the number of unit in the square, s _iBe i unit, E _Jlm(y ₀) be j cube at the center y ₀The lm rank multipole momentum moment, the span of l and m is respectively 0～9;

Leafy node from aforesaid Octree begins then, each cubical multipole momentum moment in the aforesaid Octree is delivered to the multipole momentum moment E at its father node center according to (3) formula _Jlm, till the multipole momentum moment of the second layer node of octree structure has been confirmed;

Wherein y ' is a father node cube centre coordinate;

Successively shift and disperse: with each cubical transfer and coefficient of divergence F in the array F storage octree structure _Jlm, its size is N * 10 * 10;

Utilize the data among the multipole momentum moment three-dimensional array E, in the octree structure greater than two-layer cube, confirm that through (4) formula its distant relative collects coefficient of migration F alternately _Jlm

Wherein x ' is for shifting cubical center;

Begin from second layer cube then, through (5) formula each cubical distant relative of aforesaid Octree is collected coefficient of migration alternately and diffuse to its child node cube center, till the bottom of tree construction, form each cubical coefficient of divergence F _Jlm

X wherein ₀The leafy node cube center of representing aforesaid Octree;

Like this, the node place cube close relative who is confirmed by the front contributes G _JiContribute F with the distant relative _JlmStack forms each node electric field value vector X in the physical model, shown in (6) formula;

(6)

D), the difference of each node electric field value vector X in the electric field value initial vector D through aforesaid each node and the aforesaid power frequency electric field physical model, confirm the norm of difference, if less than 10 ^-9, then reach computational accuracy, stop to calculate, otherwise repeating step c), up to the norm of difference less than 10 ^-9Till;

3), rendering space electric field equipotential line distribution plan

According to d) in satisfy to confirm each node electric field value vector X and step 2 in the physical model of precision) a) one-to-one relationship of middle nodal information array A; Draw out by work corridor starting point and begin the electric field equipotential line distribution plan on the perpendicular of work corridor, edge direction; Judgement reaches the limit value 4KV/m that whether surpasses national regulation with lower area power frequency electric field value apart from 1.8 meters on ground, and the power frequency electric field that is the work corridor above limit value 4KV/m exceeds standard.

Images