The magnetic field analytical method of gas-insulated line in a kind of tunnel
Technical field
The invention belongs to a kind of magnetic field analytical methods of gas-insulated line in GIL device fields more particularly to tunnel.
Background technology
The rationalization of the land infrastructure such as electric power and traffic on railway, highway requires reconciliation environmental problem and technical side
The technological progress in face.The transmission of this integrated power in same gallery or same structure and other services (such as railway and public affairs
Road transport, bridge, piping lane etc.) possibility realization be a technological challenge, be related to this technology at present has insulated cable
(XLPE) or gas-insulated line (GILs).
Gas-insulated line represents a kind of new power transmission technology, rather than traditional aerial and cable run.World's model
In enclosing, it is in use to have more than 100 kilometers of GIL.Most common GIL structures are made of three cylindrical shells, and shell is height
Electric conductivity aluminum or aluminum alloy is epoxyresin insulator per its center of phase conductor.The ratio of internal diameter of outer cover and outer diameter is 2.5
It, may be different in~3 ranges.
With the application development of nearest GIL devices, it is suitable for current density distribution and closing effect analysis and predicts more
Conductor system is modeled as necessity.The method being widely used at present is mainly based upon more conductor system matrix analyses, it is used first
That tests simplifies it is assumed that providing very satisfactory result.But its process of calculation analysis is more complicated, needs to simplify calculating
Method is predicted and is analyzed to current density distribution and stray magnetic field in GIL systems.
Invention content
Goal of the invention:In view of the above problems, the present invention proposes a kind of magnetic field analytical method of gas-insulated line in tunnel, letter
Change and calculate, current density distribution and stray magnetic field in GIL systems are predicted and analyzed.
Technical solution:To achieve the purpose of the present invention, the technical solution adopted in the present invention is:Gas is exhausted in a kind of tunnel
The magnetic field analytical method of edge line, including step:
(1) gas-insulated linear system is united and is divided into several division equivalent substitution conducting wires per phase and shell;
(2) gas-insulated linear system system is divided according to the constraint of transmission line of electricity transmitting terminal;
(3) longitudinal self-impedance and mutual impedance are calculated;
(4) edge-restraint condition of transmitting terminal is set;
(5) power supply external circuit model is established, current density distribution and stray magnetic field carry out pre- in uniting to gas-insulated linear system
It surveys and analyzes.
Transmission line of electricity transmitting terminal constraint:
U=Zi
In formula, u is the column voltage vector of transmission line of electricity transmitting terminal;I is the electric current column vector of auxiliary conductor;Z is system resistance
Anti- matrix.
The step (3) specifically includes:
(3.1) equivalent diameter of wire is calculated;
(3.2) additional series connection resistance per unit length is calculated;
rc=mRc
re=mRe
In formula, RcIndicate the resistance per unit length of phase conductor, ReIndicate the resistance per unit length of shell;
(3.3) longitudinal self-impedance and mutual impedance are calculated;
In formula, rs=rc or rs=re;DCAIndicate Carson's depth;F indicates frequency;ρgIndicate the resistivity of soil;dsIt indicates
The diameter of auxiliary conductor;dsrIndicate the relative distance between r and s auxiliary conductor axis;π=a/ds.
The step (4) specifically includes:
(4.1) calculating matrix:
U=TU
I=Tti
In formula, u, i are the column vectors in the matrix;
(4.2) matrix relationship of adjoint matrix is calculated:
I=Yu
In formula, Y=Z-1;
It can further be transformed to calculate column vector I, U:
In formula,
(4.3) I vectors are calculated:
It can be derived by (4.2) and calculate conductor current vector:
In formula, T changes with boundary condition.
The step (5) specifically includes:
(5.1) circuit model is stated;
(5.2) discretization is carried out;
(5.3) total current by solid conductor is calculated;
(5.4) power supply external circuit matrix is indicated.
Advantageous effect:Present invention firstly provides the magnetic field analytical methods of gas-insulated line in tunnel, to more conductor systems
Current density is distributed and kindred effect is analyzed and predicted;Integration method of the present invention simplifies the assessment of external flux density, square
Algebra of matrices method is succinct, the close reality of calculating process.
Description of the drawings
Fig. 1 is closing phase conductor schematic diagram;
Fig. 2 is conductor distribution schematic diagram;
Fig. 3 is column vector U and I distribution map before transformation;
Fig. 4 is column vector U and I distribution map after transformation.
Specific implementation mode
Technical scheme of the present invention is further described with reference to the accompanying drawings and examples.
As shown in Figure 1, gas-insulated line (GIL) structure includes three-phase conductor, it is cylinder per phase conductor, shell is highly conductive
Property aluminum or aluminum alloy, center is epoxyresin insulator.The ratio of internal diameter of outer cover and outer diameter, may in 2.5~3 ranges
It is different.
The magnetic field analytical method of gas-insulated line, includes the following steps in tunnel of the present invention:
(1) ignore kelvin effect when frequency is 50-60HZ, every phase and shell, which are divided into the division that m substitutes, equivalent to be led
Line, as shown in Fig. 2, its own longitudinal impedance and common longitudinal impedance apply Carson's theoretical calculation;
Transmission line of electricity transmitting terminal constraint:
U=Zi
In formula, u is the column voltage vector of transmission line of electricity transmitting terminal;I is the electric current column vector of auxiliary conductor;Z is system resistance
Anti- matrix (6m × 6m).
(2) it is constrained according to transmission line of electricity transmitting terminal, divides gas-insulated linear system system, as shown in Figure 2;
(2-1) calculates equivalent diameter of wire, and the thickness and resistivity per phase and the tubulose auxiliary conductor of shell are zero;
(2-2) calculates additional series connection resistance per unit length:
rc=mRc
re=mRe
In formula, RcIndicate the resistance per unit length (Ω/km) of phase conductor, ReIndicate shell resistance per unit length (Ω/
km)。
(3) longitudinal self-impedance and mutual impedance (Ω/km) are calculated:
In formula, rs=rcOr rs=re;Indicate Carson's depth, m;F indicates frequency, 50~60HZ;ρg
Indicate the resistivity of soil, Ω m;dsIndicate the diameter of auxiliary conductor, m;dsrIt is opposite between expression r and s auxiliary conductor axis
Distance, polar coordinates estimation;π=a/ds。
(4) boundary constraint of transmitting terminal is set;
The shell of three phase line securely connects, and transmitting terminal ground connection ignores impedance, UE=0.The Injection Current of phase conductor
Vectorial I1、I2、I3Known amplitude and phase.Column vector U and I is as shown in Figure 3, wherein phase voltage vector U1、U2、U3With vectorial IE
It is known.
(4-1) calculating matrix:
U=TU
I=Tti
In formula, u, i are the column vectors in the matrix.
(4-2) calculates the matrix relationship of adjoint matrix:
I=Yu
In formula, Y=Z-1, Δ Z growths (L=2~4km) with L growths.
It can further be transformed to calculate column vector I, U:
In formula,I, U distribution form in formula is as shown in Figure 4.
(4-3) calculates I vectors:
It can be derived by (4-2) and calculate conductor current vector:
In formula, T changes with boundary condition.
For independent phase bus duct, the formation of matrix is identical.Self-impedance and mutual impedance are calculated using image
Method, and assume that the earthing electric network in power plant is equipotential plane.
(5) FE circuit models are established;
The solution of electromagnetic problem is by describing the external circuit of system power supply come descriptive model.
(5-1) indication circuit model:
▽2 A z-j2πfμγA z=-μ J0
On each point for solving domain, by the magnetic vector of Z-directionA 0All electromagnetic quantities can be exported.
(5-2) carries out discretization using Galerkin programs:
Pi,j=∫ΩγNiNjdΩ
In formula, NiIt is and the relevant shape function of node i;TsolIt is the number of solid conductor, Usol,qIndicate voltage drop,
Ωsol,qIndicate solid conductor cross section.
(5-3) calculates the total current by solid conductor:
(5-4) indicates FE circuit matrixs:
By the integration method, current density distribution and stray magnetic field in GIL systems are predicted and analyzed.