CN104779571A - Power transmission line ice melting method based on gravity action calculation model - Google Patents

Power transmission line ice melting method based on gravity action calculation model Download PDF

Info

Publication number
CN104779571A
CN104779571A CN201510191765.2A CN201510191765A CN104779571A CN 104779571 A CN104779571 A CN 104779571A CN 201510191765 A CN201510191765 A CN 201510191765A CN 104779571 A CN104779571 A CN 104779571A
Authority
CN
China
Prior art keywords
ice
phi
transmission line
ice melting
wire
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
Application number
CN201510191765.2A
Other languages
Chinese (zh)
Other versions
CN104779571B (en
Inventor
陆佳政
郭俊
张红先
方针
李波
谭艳军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Hunan Electric Power Co Prevents And Reduces Natural Disasters Center
State Grid Corp of China SGCC
State Grid Hunan Electric Power Co Ltd
Original Assignee
State Grid Hunan Electric Power Co Prevents And Reduces Natural Disasters Center
State Grid Corp of China SGCC
State Grid Hunan Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Hunan Electric Power Co Prevents And Reduces Natural Disasters Center, State Grid Corp of China SGCC, State Grid Hunan Electric Power Co Ltd filed Critical State Grid Hunan Electric Power Co Prevents And Reduces Natural Disasters Center
Priority to CN201510191765.2A priority Critical patent/CN104779571B/en
Publication of CN104779571A publication Critical patent/CN104779571A/en
Application granted granted Critical
Publication of CN104779571B publication Critical patent/CN104779571B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02GINSTALLATION OF ELECTRIC CABLES OR LINES, OR OF COMBINED OPTICAL AND ELECTRIC CABLES OR LINES
    • H02G7/00Overhead installations of electric lines or cables
    • H02G7/16Devices for removing snow or ice from lines or cables

Abstract

The invention discloses a power transmission line ice melting method based on a gravity action calculation model. The power transmission line ice melting method comprises the following steps: acquiring meteorological factors and wire factors related to ice melting; calculating relation between ice melting current and ice melting time according to an ice melting calculation formula; determining the corresponding ice melting current and ice melting time to melt ice of the power transmission line. The power transmission line ice melting method disclosed by the invention fully considers the influence of the gravity action on the ice melting process, is based on the thermal equilibrium and the heat exchange principle to deduce the ice melting calculation formula of the power transmission line, enables the calculation result to be more suitable for the ice melting practice of the power transmission line, provides more accurate data basis for ice resistance of the line, and is clear in train of thought, convenient to calculate and high in practicability.

Description

Based on the transmission line de-icing method of Action of Gravity Field computation model
Technical field
The present invention relates to electrical engineering technical field, be specifically related to a kind of transmission line de-icing method based on Action of Gravity Field computation model.
Background technology
Ice damage is the major hidden danger of electric power netting safe running.At the beginning of 2008, south China area has met with serious ice damage, and grid loss is heavy.And be subject to the impact that extreme climate in recent years constantly takes place frequently, still likely outburst is similar to the especially severe ice damage of 2008.For ensureing power network safety operation, needing after electrical network icing to carry out line ice-melting work in time, therefore needing to carry out the research of transmission line de-icing computation model.Traditional wire icing ice-melt computation model is generally according to wire icing calculation of total, but this computational methods do not take into full account the impact of Action of Gravity Field in deicing processes, and ice-melt result of calculation and actual conditions exist certain difference, causes accuracy not high.
Summary of the invention
The technical problem to be solved in the present invention is the deficiency overcoming prior art existence, a kind of transmission line de-icing method based on Action of Gravity Field computation model is provided, take into full account the impact of having melted Action of Gravity Field in ice process, transmission line de-icing computing formula is derived based on heat balance and heat exchange principle, its result of calculation more meets transmission line de-icing reality of work, for circuit anti-ice provides quantity basis more accurately.The clear thinking of the method, convenience of calculation, practical.
For solving the problems of the technologies described above, the present invention by the following technical solutions:
Based on a transmission line de-icing method for Action of Gravity Field computation model, it is characterized in that: comprise the following steps:
(1) obtain and the relevant meteorologic factor of ice-melt and wire factor, wherein meteorologic factor comprises air themperature, air viscosity and wind speed, and wire factor comprises diameter of wire, conductor resistance and wire icing thickness;
(2) calculate the relation of ice melting current and ice-melt time according to ice-melt computing formula, ice-melt computing formula is:
I = η 8 T · 1 R 0 · K · ρ ice · ( π Φ 2 - Φ 2 · arccos φ Φ + 2 φ Φ 2 - φ 2 - π φ 2 ) - π R 0 λ θ a [ 0.64 ( UΦ ν ) 0.2 + 0.5 ( UΦ ν ) 0.61 ]
In formula, I is ice melting current; T is the ice-melt time; η is correction factor, the span of η be (0,1.5]; R 0for conductor resistance; K is the ablation heat of ice, K=3.35 × 105J/kg; ρ icefor the density of icing; Φ is the diameter after wire icing; φ is diameter of wire; λ is air conduction coefficient; θ afor air themperature; U is wind speed; ν is air viscosity;
(3) according to step (2) ice melting current that calculates and the relation of ice-melt time, determine that corresponding ice melting current and ice-melt time carry out ice-melt to transmission line.
Compared with prior art, the invention has the advantages that: the present invention is directed to the deficiency that traditional wire icing ice-melt computation model does not take into full account Action of Gravity Field impact in deicing processes, propose a kind of transmission line de-icing method based on Action of Gravity Field computation model, take into full account the impact of having melted Action of Gravity Field in ice process, transmission line de-icing computing formula is derived based on heat balance and heat exchange principle, its result of calculation more meets transmission line de-icing reality of work, for circuit anti-ice provides quantity basis more accurately.The clear thinking of the method, convenience of calculation, practical.
Accompanying drawing explanation
Fig. 1 is the structural representation of wire icing in transmission line de-icing method of the present invention.
Fig. 2 is the structural representation of wire ice-melt in transmission line de-icing method of the present invention.
Fig. 3 is the structural representation in wire ice-melt region in transmission line de-icing method of the present invention.
Marginal data:
1, wire; 2, icing.
Embodiment
Below in conjunction with the drawings and specific embodiments, the present invention is described in further detail.
Based on a transmission line de-icing method for Action of Gravity Field computation model, comprise the following steps:
(1) obtain and the relevant meteorologic factor of ice-melt and wire factor, wherein meteorologic factor comprises air themperature, air viscosity and wind speed, and wire factor comprises diameter of wire, conductor resistance and wire icing thickness;
(2) calculate the relation of ice melting current and ice-melt time according to ice-melt computing formula, ice-melt computing formula is:
I = η 8 T · 1 R 0 · K · ρ ice · ( π Φ 2 - Φ 2 · arccos φ Φ + 2 φ Φ 2 - φ 2 - π φ 2 ) - π R 0 λ θ a [ 0.64 ( UΦ ν ) 0.2 + 0.5 ( UΦ ν ) 0.61 ]
In formula, I is ice melting current; T is the ice-melt time; η is correction factor, the span of η be (0,1.5]; R 0for conductor resistance; K is the ablation heat of ice, K=3.35 × 105J/kg; ρ icefor the density of icing; Φ is the diameter after wire icing; φ is diameter of wire; λ is air conduction coefficient; θ afor air themperature; U is wind speed; ν is air viscosity; Above-mentioned correction factor is obtained by repetition test, the corresponding a kind of correction factor of each wire;
(3) according to step (2) ice melting current that calculates and the relation of ice-melt time, determine that corresponding ice melting current and ice-melt time carry out ice-melt to transmission line.
The derivation of above-mentioned ice-melt computing formula is as follows:
(1) heat exchange amount of wire and air calculates
Suppose that air themperature is θ a(DEG C), conductor temperature is θ b(DEG C), diameter of wire is φ (m), and after icing, diameter is Φ (m), and when wire generation icing, wire can be approximated to be cylinder.With this understanding, the heat exchange amount Δ Q (W/m) between unit length wire and air calculates by following formula:
ΔQ=πh(θ ba)Φ (1)
In formula, h is heat conduction exchange coefficient, and computing formula is as follows:
h = λ Nu Φ - - - ( 2 )
In formula, λ (W/m DEG C) is air conduction coefficient; Nu is Nu Saierte coefficient, can be expressed as the function of the Reynolds number Re of cylinder circumference turbulent airflow, as shown in the formula:
Nu=0.64Re 0.2+0.2Re 0.61(3)
Re = UΦ ν - - - ( 4 )
In formula, U (m/s) is wind speed, ν (m 2/ s) be air viscosity.
Formula (2)-(4) are substituted into formula (1), and the general computing formula that can obtain the heat exchange amount between unit length wire and air is:
ΔQ = πλ ( θ b - θ a ) [ 0.64 ( UΦ ν ) 0.2 + 0.2 ( UΦ ν ) 0.61 ] - - - ( 5 )
Temperature due to mixture of ice and water is 0 DEG C, when wire generation icing, and θ b=0, formula (5) can be reduced to:
ΔQ = - πλ θ a [ 0.64 ( UΦ ν ) 0.2 + 0.2 ( UΦ ν ) 0.61 ] - - - ( 6 )
From formula (6), the heat exchange amount between wire and air is mainly by the impact of air themperature, wind speed, diameter of wire, air conduction coefficient and air viscosity.
(2), wire icing melts the required heat Calculation that comes off
Heat needed for traditional wire icing melts, according to wire icing calculation of total, is also the heat in Fig. 1 needed for the thawing of dash area icing total amount.But due to the impact of Action of Gravity Field, wire icing melts to come off only needs the icing on wire top to melt, and is also the wire icing of Fig. 2 bend dash area.
Therefore, unit length wire icing melts the required heat energy total amount Q that comes off wcomputing formula as follows:
Q W=Q ice+ΔQ (7)
In formula, Δ Q is the exchange heat of wire and air, also instant heating loss of energy amount; Q icefor the icing of Fig. 2 bend dash area melts institute's calorific requirement, its computing formula is as follows:
Q ice=K·ρ ice·V ice/T (8)
In formula, ρ icefor the density of icing; K is the ablation heat of ice, K=3.35 × 10 5j/kg; T (s) is ice-melt duration.
The diagonal line hatches part icing volume V of unit length wire icecomputing formula is as follows:
V ice = 1 8 π Φ 2 - 1 8 Φ 2 · arccos φ Φ + φ 4 Φ 2 - φ 2 - 1 8 π φ 2 - - - ( 9 )
Formula (8), (9) and (6) are substituted into formula (7), and the heat energy calculation of total formula that can obtain needed for wire ice-melt is:
Q W = 1 T · K · ρ ice · ( 1 8 π Φ 2 - 1 8 Φ 2 · arccos φ Φ + φ 4 Φ 2 - φ 2 - 1 8 π φ 2 ) - πλ θ a [ 0.64 ( UΦ ν ) 0.2 + 0.2 ( UΦ ν ) 0.61 ] - - - ( 10 )
(3), transmission line de-icing formula proving
Wire ice-melt heat energy is mainly derived from current in wire and produces, the heat energy Q that unit length wire discharges when being I by current strength mfor:
Q m=I 2R 0(11)
In formula, R 0for conductor resistance.
For meeting transmission line de-icing demand, then there is Q m=Q w.Accordingly, can calculate transmission line de-icing electric current I is:
I = η 8 T · 1 R 0 · K · ρ ice · ( π Φ 2 - Φ 2 · arccos φ Φ + 2 φ Φ 2 - φ 2 - π φ 2 ) - π R 0 λ θ a [ 0.64 ( UΦ ν ) 0.2 + 0.5 ( UΦ ν ) 0.61 ] - - - ( 12 )
In actual deicing processes, along with continuing of ice-melt, the space after ice-melt increases, and due to the thermal conductivity of air, the ice face directly do not contacted with wire also can be melted gradually, and the thawing region of icing is as shown in Fig. 3 bend shadow region.
Therefore, ice-melt computing formula needs to increase correction factor η and revises.Revised ice-melt computing formula is:
I = η 8 T · 1 R 0 · K · ρ ice · ( π Φ 2 - Φ 2 · arccos φ Φ + 2 φ Φ 2 - φ 2 - π φ 2 ) - π R 0 λ θ a [ 0.64 ( UΦ ν ) 0.2 + 0.5 ( UΦ ν ) 0.61 ]

Claims (1)

1., based on a transmission line de-icing method for Action of Gravity Field computation model, it is characterized in that: comprise the following steps:
(1) obtain the meteorologic factor relevant to transmission line de-icing and wire factor, wherein meteorologic factor comprises air themperature, air viscosity and wind speed, and wire factor comprises diameter of wire, conductor resistance and wire icing thickness;
(2) calculate the relation of ice melting current and ice-melt time according to ice-melt computing formula, ice-melt computing formula is:
I = η 8 T · 1 R 0 · K · ρ ice · ( πΦ 2 - Φ 2 · arccos φ Φ + 2 φ Φ 2 - φ 2 - πφ 2 ) - π R 0 λ θ a [ 0.64 ( UΦ v ) 0.2 + 0.2 ( UΦ v ) 0.61 ]
In formula, I is ice melting current; T is the ice-melt time; η is correction factor, the span of η be (0,1.5]; R 0for conductor resistance; K is the ablation heat of ice, K=3.35 × 105J/kg; ρ icefor the density of icing; Φ is the diameter after wire icing; φ is diameter of wire; λ is air conduction coefficient; θ afor air themperature; U is wind speed; ν is air viscosity;
(3) according to step (2) ice melting current that calculates and the relation of ice-melt time, determine that corresponding ice melting current and ice-melt time carry out ice-melt to transmission line.
CN201510191765.2A 2015-04-22 2015-04-22 Based on the transmission line de-icing method of Action of Gravity Field computation model Active CN104779571B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510191765.2A CN104779571B (en) 2015-04-22 2015-04-22 Based on the transmission line de-icing method of Action of Gravity Field computation model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510191765.2A CN104779571B (en) 2015-04-22 2015-04-22 Based on the transmission line de-icing method of Action of Gravity Field computation model

Publications (2)

Publication Number Publication Date
CN104779571A true CN104779571A (en) 2015-07-15
CN104779571B CN104779571B (en) 2016-03-09

Family

ID=53620884

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510191765.2A Active CN104779571B (en) 2015-04-22 2015-04-22 Based on the transmission line de-icing method of Action of Gravity Field computation model

Country Status (1)

Country Link
CN (1) CN104779571B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106099751A (en) * 2016-07-07 2016-11-09 周丽玲 A kind of transmission line de-icing method based on action of gravity computation model
CN107134740A (en) * 2017-06-29 2017-09-05 国网四川省电力公司电力科学研究院 The anti-icing method and system of extra high voltage direct current transmission line on-load
CN109119931A (en) * 2018-08-06 2019-01-01 四川大学 Transmission line of electricity anti-icing ice-melt heat Calculation method online based on self-control thermal wire
CN109449852A (en) * 2018-03-23 2019-03-08 国网浙江省电力公司丽水供电公司 A kind of distribution line de-icing method of non-power loss mode

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0937448A (en) * 1995-07-18 1997-02-07 Nissin Electric Co Ltd Snow-melting method of transmission line
CN101459326A (en) * 2008-12-31 2009-06-17 南方电网技术研究中心 Determination method for principal parameters of DC ice melting device
CN103199476A (en) * 2013-04-10 2013-07-10 中国电力工程顾问集团西南电力设计院 Method and application for acquiring direct current load value of direct current deicing device of convertor station ground wires

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0937448A (en) * 1995-07-18 1997-02-07 Nissin Electric Co Ltd Snow-melting method of transmission line
CN101459326A (en) * 2008-12-31 2009-06-17 南方电网技术研究中心 Determination method for principal parameters of DC ice melting device
CN103199476A (en) * 2013-04-10 2013-07-10 中国电力工程顾问集团西南电力设计院 Method and application for acquiring direct current load value of direct current deicing device of convertor station ground wires

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106099751A (en) * 2016-07-07 2016-11-09 周丽玲 A kind of transmission line de-icing method based on action of gravity computation model
CN107134740A (en) * 2017-06-29 2017-09-05 国网四川省电力公司电力科学研究院 The anti-icing method and system of extra high voltage direct current transmission line on-load
CN109449852A (en) * 2018-03-23 2019-03-08 国网浙江省电力公司丽水供电公司 A kind of distribution line de-icing method of non-power loss mode
CN109119931A (en) * 2018-08-06 2019-01-01 四川大学 Transmission line of electricity anti-icing ice-melt heat Calculation method online based on self-control thermal wire
CN109119931B (en) * 2018-08-06 2020-01-24 四川大学 Power transmission line online anti-icing and de-icing heat quantity calculation method based on self-made thermal conductor

Also Published As

Publication number Publication date
CN104779571B (en) 2016-03-09

Similar Documents

Publication Publication Date Title
Parent et al. Anti-icing and de-icing techniques for wind turbines: Critical review
Messinger Equilibrium temperature of an unheated icing surface as a function of air speed
CN103208665B (en) The control method of battery pack temperature control system and corresponding battery pack temperature
Makkonen et al. Modelling and prevention of ice accretion on wind turbines
Battisti et al. Aerodynamic measurements on a vertical axis wind turbine in a large scale wind tunnel
US20150023792A1 (en) Wind turbine rotor blade de-icing process and wind turbine rotor blade de-icing system
Wang Recent progress on ultrasonic de-icing technique used for wind power generation, high-voltage transmission line and aircraft
CN103324215B (en) The method to set up of Intelligent load management system for transformer thermal circuit model
Struk et al. Fundamental ice crystal accretion physics studies
Kraj et al. Measurement method and results of ice adhesion force on the curved surface of a wind turbine blade
CN102136709B (en) Device and method for reducing ice and snow hung and accumulated on high-voltage electric overhead conductors
CN105089929B (en) Wind generator set blade icing detecting system and its method
CN104462624A (en) Motor temperature field data processing method based on multi-field coupling
CN205067511U (en) Special resistive connection ice anemometry sensor heating system of anemometry equipment
EP2163475B1 (en) Propeller deicing system
CN105956286A (en) Prediction method of forecabin thermal protection system whole trajectory temperature boundary of hypersonic velocity aircraft
CN102438903A (en) Icing sensor system and method
Petrenko et al. Pulse electro-thermal de-icer (PETD)
CN102878026B (en) Wind driven generator rotor blade with electrothermal deicing devices
CN102759713B (en) A kind of energy content of battery efficiency test device and method of testing thereof
CN101590914A (en) A kind of anti-icing and deicing device for flight vehicle
CN101907456B (en) Method for calculating thickness and weight of ice coating on overhead transmission line of tangent tower
CN204156540U (en) A kind of solar module ice melting system
CN103383428B (en) A kind of pole line cable hybrid line both-end Method of Traveling Wave Fault Ranging
CN204290227U (en) Road deicing device

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Fang Yifei

Inventor after: Lu Jiazheng

Inventor after: Guo Jun

Inventor after: Zhang Hongxian

Inventor after: Fang Zhen

Inventor after: Li Bo

Inventor after: Tan Yanjun

Inventor before: Lu Jiazheng

Inventor before: Guo Jun

Inventor before: Zhang Hongxian

Inventor before: Fang Zhen

Inventor before: Li Bo

Inventor before: Tan Yanjun