CN108038268A - GIL running temperature appraisal procedures in a kind of corridor pipe - Google Patents
GIL running temperature appraisal procedures in a kind of corridor pipe Download PDFInfo
- Publication number
- CN108038268A CN108038268A CN201711149403.2A CN201711149403A CN108038268A CN 108038268 A CN108038268 A CN 108038268A CN 201711149403 A CN201711149403 A CN 201711149403A CN 108038268 A CN108038268 A CN 108038268A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- conductor
- mtr
- mtd
- 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
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000004020 conductor Substances 0.000 claims abstract description 200
- 239000011257 shell material Substances 0.000 claims description 99
- 239000003570 air Substances 0.000 claims description 51
- 238000004519 manufacturing process Methods 0.000 claims description 45
- 230000017525 heat dissipation Effects 0.000 claims description 39
- 238000012546 transfer Methods 0.000 claims description 27
- 230000000704 physical effect Effects 0.000 claims description 22
- 230000002500 effect on skin Effects 0.000 claims description 15
- 230000005855 radiation Effects 0.000 claims description 13
- 230000000694 effects Effects 0.000 claims description 12
- 239000002699 waste material Substances 0.000 claims description 9
- 238000010292 electrical insulation Methods 0.000 claims description 7
- 230000005611 electricity Effects 0.000 claims description 7
- 239000000463 material Substances 0.000 claims description 5
- 239000000203 mixture Substances 0.000 claims description 5
- 238000011156 evaluation Methods 0.000 claims description 4
- 239000007787 solid Substances 0.000 claims description 4
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 3
- 239000011343 solid material Substances 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 2
- 239000012080 ambient air Substances 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 238000003825 pressing Methods 0.000 claims description 2
- 238000009423 ventilation Methods 0.000 claims description 2
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 229910017435 S2 In Inorganic materials 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 238000009413 insulation Methods 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 10
- 238000012544 monitoring process Methods 0.000 abstract description 8
- 230000005540 biological transmission Effects 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000012212 insulator Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 239000003292 glue Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001869 rapid Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
GIL running temperature appraisal procedures in a kind of corridor pipe, establish GIL physical models, calculate the loss of shell initial power, calculate thermal convection current and the thermal exposure of shell and air, calculate the thermal convection current between conductor and shell and thermal exposure, for the present invention by extra-high voltage GIL thermal characteristics MATLAB calculation procedures in the corridor pipe that is designed in the MATLAB environment of Windows operating system, temperature during using the GIL relevant featuring parameters under initial operating mode come to its stable operation carries out assessment calculating；The present invention can assess the steady-state operating temperature for calculating extra-high voltage GIL in corridor pipe under different initial working conditions；Present invention assessment calculating speed is fast, and precision is high, can be embedded into GIL temperature monitoring systems, is of great significance to improving GIL thermal reliabilities.
Description
Technical field
The present invention relates to the transmission of electricity of extra-high voltage electrical, gas-insulated transmission line (GIL) thermal characteristics, long range extra-high voltage GIL
The fields such as equipment thermal reliability, more particularly to GIL running temperature appraisal procedures in a kind of corridor pipe.
Background technology
Gas-insulated transmission line (GIL) is a kind of using long-distance electricity transmission device of the gas as dielectric, with biography
System cable power transmission is compared, and has the advantages that capacity is big, and loss is low, and service life is long.Common GIL failures are as discharged, insulator
Phenomena such as breakdown is all accompanied by conductor or abnormal skin temperature when occurring.Moreover, GIL bus temperatures are excessive also occurs that
Busbar vault, or even the failure such as disc insulator rupture, basin marine glue cracking.Heat of the research of GIL thermal characteristics to raising GIL
Reliability and guarantee GIL safe and stable operations are of great significance.It is at this stage to use finite element method to GIL thermal characteristics more
Studied, but finite element simulation method is computationally intensive, generally cannot function as subprogram and is embedded into GIL temperature monitoring systems
In.
Therefore, the proposition one kind how innovated fast and effectively in corridor pipe GIL running temperatures appraisal procedure to meet reality
The demand of application is the most important thing.
The content of the invention
The technical problems to be solved by the invention are that solving existing finite element numerical simulation computational methods quick cannot obtain
GIL running temperatures in corridor pipe, and the problem of subprogram is embedded into GIL real-time temperature test sytems is cannot function as, so that
It is proposed GIL running temperature appraisal procedures in a kind of corridor pipe, the thermal characteristics of GIL is solved using analytic method, can not only
Quickly GIL thermal characteristics problems are solved, subprogram is alternatively arranged as and is embedded into GIL temperature monitoring systems.GIL is run
The quick of temperature grasps the thermal reliability for being more conducive to improve GIL and more preferably ensures GIL safe and stable operations, and this method can also be made
It is embedded into for subprogram in GIL real-time temperature test sytems.
To solve the above problems, the technical scheme is that：
GIL running temperature appraisal procedures, comprise the following steps in a kind of corridor pipe：
Step S1：GIL physical models are established, model is related to resistance joule heat effect, the free convection layer of outside air
Stream flowing forces turbulent convection, built-in electrical insulation gas (such as SF6Gas) electricity such as the Laminar Flow of free convection, radiant heat transfer,
Heat, stream, solid multiple physical field, are a multiple physical field coupling models.Therefore GIL running temperature parameters refer in this model evaluation corridor pipe
Mark is included between convection current heat output and thermal exposure, shell and air between joule heat waste, conductor and the shell of conductor and shell
Convection current heat output and thermal exposure.
Step S2：The loss of conductor initial power is calculated, according to above-mentioned physical model, is related to the fluid and solid material of calculating
Matter attribute considers the change with local temperature, i.e. outside air, built-in electrical insulation gas, conductor and enclosure material thermal physical property parameter
It is the function formula of temperature, the enclosure material thermal physical property parameter includes viscosity, and the coefficient of conductivity is permanent than pressing, than
Thermal capacitance, thermal coefficient of expansion, so that the numerical value of hot Physical Parameters is closer to operation actual conditions when calculating；In actually solving
According to the thermal physical property parameter and skin effect coefficient of conductor, the initial temperature of conductor is determined, it is assumed that initial temperature and environment temperature
It is equal, iterative solution under GIL pipelines startup operating condition is calculated according to given initial duty parameter and goes out corresponding public condition lower conductor
Heat amount；
Step S3：The loss of shell initial power is calculated, the collection of Shell resistance is considered when calculating Shell resistance joule heat amount
Skin effect, i.e., on shell section, central current density is small, big by the current density of near surface, according to the physical model of step S1,
After the thermal physical property parameter and skin effect coefficient of shell is determined, the initial temperature of conductor is determined, according to given initial work
Condition parameter calculates the heat amount that iterative solution under GIL pipeline start operating performances goes out corresponding public condition lower conductor.
Step S4：Thermal convection current and the thermal exposure of shell and air are calculated, judges quantity of heat production and the heat dissipation of GIL conductors entirety
Whether amount error is less than program setting error, error be no more than the quantity of heat production of 5%, GIL conductors entirety by conductor quantity of heat production with
The quantity of heat production composition of shell；The heat dissipation of GIL conductors is carried out by the shell directly contacted with outside air, therefore GIL conductors are overall
Heat dissipation capacity be made of the convection current heat output between shell and outside air and radiant heat transfer amount；Calculate the shell being related to and the external world
Whether the thermal physical property parameter of air is the function of temperature, small by the quantity of heat production and heat dissipation capacity error that judge GIL conductors entirety
Judge whether the heat production of GIL conductors entirety and heat dissipation reach balance in setting error, i.e. whether GIL conductors operation has reached
Stable state.
Step S5：If step S4 errors are less than setting error, next step is carried out, otherwise shell initial temperature increases
0.1 DEG C of return to step S3；
Step S6：The thermal convection current between conductor and shell and thermal exposure are calculated, the quantity of heat production and heat dissipation capacity for judging conductor miss
Whether difference is less than setting error；The direct representation of the power attenuation of conductor is the quantity of heat production of conductor, the heat dissipation master of conductor
To be carried out by insulating gas, therefore the heat dissipation capacity of conductor is by the free convection heat transfer amount between conductor and shell and radiant heat transfer amount group
Into the calculating of thermal convection current and thermal exposure between conductor and shell relates generally to the running temperature of conductor and shell, and insulate
The thermal physical property parameters such as the dynamic viscosity of gas, thermal conductivity factor, specific heat at constant pressure, coefficient of cubical expansion, wherein all hot physical property
Parameter is the function of temperature.Judge conductor by the way that whether the quantity of heat production and heat dissipation capacity error that judge conductor are less than setting error
Whether heat production and heat dissipation reach balance, i.e. whether conductor operation has reached stable state.
Step S7：If step S6 errors are less than setting error, calculate and terminate output as a result, otherwise conductor temperature increases
0.1 DEG C of return to step S2.When the quantity of heat production of conductor is less than setting error with heat dissipation capacity error, illustrate that the operation of conductor reaches steady
State, can obtain the running temperature of conductor during stable operation.If the quantity of heat production of conductor is more than setting error with heat dissipation capacity error, i.e.,
Conductor operation is not up to stable state, and need increases conductor temperature at this time is iterated calculating again, until error meets that condition can jump
Go out circulation, close to the temperature of actual conductor and shell when output GIL is stable.
Convection transfer rate of the unknown quantity that GIL physical models are related in the step S1 including air and insulating gas,
Thermal diffusivity, thermal conductivity, dynamic viscosity, specific heat capacity etc. are the thermal physical property parameter of temperature funtion, are the conductor of temperature funtion and outer
The thermal conductivity factor and skin effect coefficient of shell, by using the parameter contacted with temperature, can make result more closing to reality.
In the step S1 initial duty parameter include conductor loading electric current, corridor tube environment temperature, GIL pipelines internal pressure,
Corridor inner air tube flow velocity, GIL pipeline configuration parameters.
In the step S2, conductor considers the collection of conductor there are kelvin effect when calculating conductor resistance joule heat amount
Skin effect；After determining conductor initial temperature, it can determine that conductor initial power is lost according to initial working condition, wherein correlation computations
Formula is：
In formula, RcFor the resistance of conductor, KfcFor the skin effect coefficient of conductor, ρc20For 20 DEG C when conductor D.C. resistance
Rate, TcFor the running temperature of conductor, αc20For the temperature-coefficient of electrical resistance of conductor, ScFor the cross-sectional area of conductor, DcFor conductor diameter,
CcFor conductor wall thickness, PcFor every meter of loss of conductor, I is conductor rated current.
In the step S3, shell considers the collection of shell there are kelvin effect when calculating Shell resistance joule heat amount
Skin effect；After determining shell initial temperature, it can determine that shell initial power is lost according to initial working condition.Wherein correlation computations
Formula is：
In formula, RtFor the resistance of shell, KftFor the skin effect coefficient of shell, ρtaFor 20 DEG C when shell D.C. resistance
Rate, TtFor the running temperature of shell, αtaFor the temperature-coefficient of electrical resistance of shell, StFor the cross-sectional area of shell, DtFor envelope outer diameter,
CtFor shell wall thickness, PtFor every meter of loss of shell, I is shell rated current.
In the step S4, the quantity of heat production of GIL conductors entirety is made of the quantity of heat production of conductor and the quantity of heat production of shell, heat dissipation
Amount is made of the convection current heat output between shell and outside air and radiant heat transfer amount, and correlation computations formula is：
In formula, QtcFor every meter of heat loss through convection amount of shell, QtrFor every meter of heat loss through radiation amount of shell, αaExchange heat for cross-ventilation and be
Number, TaFor ambient air temperature, εtFor inner surface of outer cover blackness, Nu is the nusselt number of outside air, λaFor the heat conduction system of air
Number, GraFor the Prandtl number of air, g is acceleration of gravity, βaFor the coefficient of cubical expansion of air, μaGlued for the power of air
Degree, TdaFor air setting temperature, Δ1For the quantity of heat production and heat dissipation capacity error of GIL conductors.
In the step S6, the power attenuation of conductor is the quantity of heat production of conductor, and the heat dissipation capacity of conductor is by conductor and shell
Between free convection heat transfer amount and radiant heat transfer amount composition, correlation computations formula is：
In formula, QccFor every meter of heat loss through convection amount of conductor, QcrFor every meter of heat loss through radiation amount of conductor, εcIt is black for conductor outside surfaces
Degree, εeFor the suitable blackness of full group object, λeFor equivalent heat conductivity, λ is insulating gas thermal conductivity factor, GrsFor insulating gas lattice
La Xiaofu numbers, PrsFor insulating gas Prandtl number, CpFor insulating gas specific heat at constant pressure, μ is insulating gas viscosity, βsInsulate gas
The coefficient of cubical expansion of body, RtiFor housing interior radius, RcoFor conductor outer radius, TdsFor the setting temperature of insulating gas, Δ2For
The quantity of heat production of conductor and heat dissipation capacity error.
In step S5, the progress that the temperature unit value of increase can have more GIL temperature measuring equipments in actual conditions is configured, and is counted
It is high to calculate precision.
In step S7, the progress that the temperature unit value of increase can have more GIL temperature measuring equipments in actual conditions is configured；It is defeated
The qualitative temperature of the insulating gas of conductor, the temperature of shell and needs when the result gone out runs to stable state for GIL.
GIL running temperature appraisal procedures in the corridor pipe of the present invention, based on law of conservation of energy foundation with for node heat
Network model is transmitted, passes through extra-high voltage GIL thermal characteristics in the corridor pipe that is designed in the MATLAB environment of Windows operating system
MATLAB calculation procedures, can using the GIL relevant featuring parameters under different initial operating modes come to its stable operation when temperature
Carry out assessment calculating.
The present invention carries out assessment calculating by designing MATLAB programs, for the program for predicting GIL thermal characteristics is embedded into GIL
Feasible means are provided in temperature monitoring system.The heat transfer that this method can conveniently try to achieve GIL under different operating modes is special
Property and the temperature of conductor and shell, and then have from bringing into operation to thermal characteristics when stablizing a research well to GIL, its
The calculation procedure of realization, easy to operate, highly practical, calculating speed is fast, can set initial parameter according to different operating mode, to be right
The research of GIL thermal characteristics provides good theoretical foundation, and can be embedded into the temperature monitoring system of GIL as subprogram
In, the reliability of GIL, which plays an important role, to be realized to raising.
Brief description of the drawings
Fig. 1 is the flow diagram of GIL running temperature appraisal procedures in the corridor pipe of the present invention.
Embodiment
The present invention is described in detail with reference to width figure.
The present invention is primarily intended to carry out mould by analytical Calculation to the pipe GIL thermal characteristics heat production of extra-high voltage corridor and diabatic process
Intend, first the GIL structures calculated are carried out simplifying processing：First, ignore power attenuation and the heat transfer of contact part, second,
Only the single-phase GIL pipeline long to one meter carries out the research of thermal characteristics.GIL busbares are powered after operation, not only on conductor there are electric current,
Also have on shell with the big reverse sensing electric current such as load current, therefore conductor and shell can produce power attenuation.Conductor and
Heat transfer between shell is carried out by the convection current and radiation of insulating gas between the two, the heat of shell and outside air
Transmission is carried out by the convection current and radiation of outside air.From law of conservation of energy, during in thermal steady state, conductor
Loss pass to shell in a manner of radiation mode and free convection, GIL overall loss by radiate with free convection in a manner of transmit
To surrounding air, the two is in equilibrium state.The major parameter index of GIL running temperatures includes conductor in this model evaluation corridor pipe
The convection current heat output between the convection current heat output and thermal exposure, shell and air between joule heat waste, conductor and shell with shell
And thermal exposure.
The GIL busbar conservation of energys：
Pc+Pt=Qtc+Qtr (5)
Pc=Qcc+Qcr (6)
Wherein PcAnd PtRespectively every meter of power attenuation of conductor and shell, W/m；QtcFor the free convection of shell and air
Heat exchange amount, W/m；QtrFor the heat loss through radiation amount of outer surface of outer cover, W/m；QccFor conductor and the quantity of heat convection of shell, W/m；Qcr
For the heat loss through radiation amount of conductor outside surfaces, W/m.
The calculating of main assessment parameter：
Pc=I2Rc (7)
Pt=I2Rt (8)
Qtc=αaπDt(Tt-Ta) (11)
Wherein, Tc is the temperature of conductor, and Tt is skin temperature, and Ta is original ambient temperature.It can be obtained from formula, conductor
Temperature, skin temperature, environment temperature and the loss of GIL and heat output are closely bound up.
A kind of GIL running temperatures appraisal procedure schematic diagram in the corridor pipe provided referring to Fig. 1, there is provided GIL operations temperature in corridor pipe
Appraisal procedure is spent, the steady-state operating temperature of extra-high voltage GIL in corridor pipe under different initial working conditions is calculated for assessing, it is assessed
Calculation procedure is as follows：
Step S1：GIL physical models are established according to the related thermal physical property parameter of air, insulating gas, conductor and shell；It is female
Line conductor and shell are heater element, and shell outwardly distributes heat by free convection and radiation.In order to reflect true heat transfer
Process, considers the free convection of the heat transfer mode, i.e. built-in electrical insulation air-flow of free convection and radiation between conductor and shell
The heat radiation and absorption of heat exchange, conductor outside surfaces and inner surface of outer cover, because it was assumed that conductor and shell axial and circumferential temperature
Uniformly, only consider the radial direction heat transfer of conductor and shell in solid domain, be hereby based on law of conservation of energy and establish one-dimensional node
Heat transfer network model.Model is related to resistance joule heat effect, the free convection Laminar Flow of outside air or forces rapids
Flow convection current, built-in electrical insulation gas (such as SF6Gas) electricity, heat, stream, solid more physics such as the Laminar Flow of free convection, radiant heat transfer
, it is a highly complex multiple physical field coupling model.Therefore GIL running temperature parameter index bags in this model evaluation corridor pipe
Include the convection current between the convection current heat output and thermal exposure, shell and air between joule heat waste, conductor and the shell of conductor and shell
Heat output and thermal exposure.
Step S2：Determine conductor initial temperature and initial working condition, calculate the loss of conductor initial power；With it is most of often
Resistance is calculated with cable and is taken as constant difference, and the kelvin effect of conductor resistance is considered when calculating conductor resistance joule heat amount
(i.e. on cross-sectional area of conductor, central current density is small, big by the current density of near surface).According to above-mentioned physical model, it is related to meter
The fluid and solid material attribute of calculation consider the change with local temperature, i.e. outside air, built-in electrical insulation gas, conductor and outer
Shell material matter thermal physical property parameter is the function formula of temperature, and the enclosure material thermal physical property parameter includes viscosity, conduction
Coefficient, it is permanent than pressure, specific heat capacity, thermal coefficient of expansion, so that the numerical value of hot Physical Parameters is actual closer to operation when calculating
Situation；According to the thermal physical property parameter and skin effect coefficient of conductor in actually solving, the initial temperature of conductor is determined, it is assumed that just
Beginning temperature is equal with environment temperature, and calculating GIL pipelines according to given initial duty parameter starts iterative solution under operating condition
Go out the heat amount of corresponding public condition lower conductor；
Step S3：Determine shell initial temperature and initial working condition, calculate the loss of shell initial power；Calculate shell electricity
The kelvin effect of Shell resistance is considered when hindering joule heat amount, i.e., on shell section, central current density is small, by the electricity of near surface
Current density is big, according to the physical model of step S1, after the thermal physical property parameter and skin effect coefficient of shell is determined, determines to lead
The initial temperature (initial temperature is equal with environment temperature) of body, calculates GIL pipelines according to given initial duty parameter and starts work
Iterative solution goes out the heat amount of corresponding public condition lower conductor under condition.
Step S4：Thermal convection current and the thermal exposure of shell and air are calculated, judges quantity of heat production and the heat dissipation of GIL conductors entirety
Whether amount error is less than program setting error (being usually no more than 5%)；The quantity of heat production of GIL conductors entirety by conductor quantity of heat production and
The quantity of heat production composition of shell；The heat dissipation of GIL conductors is carried out by the shell directly contacted with outside air, therefore GIL conductors are overall
Heat dissipation capacity be made of the convection current heat output between shell and outside air and radiant heat transfer amount；Calculate the shell being related to and the external world
Whether the thermal physical property parameter of air is the function of temperature, small by the quantity of heat production and heat dissipation capacity error that judge GIL conductors entirety
Judge whether the heat production of GIL conductors entirety and heat dissipation reach balance in setting error, i.e. whether GIL conductors operation has reached
Stable state.
Step S5：If step S4 errors are less than setting error, next step is carried out, otherwise shell initial temperature increases
0.1 DEG C of return to step S3；When the quantity of heat production and heat dissipation capacity error of GIL conductors entirety are less than setting error, illustrate that GIL conductors are transported
Row has reached stable state, can obtain the running temperature of shell during stable operation；If the quantity of heat production and heat dissipation capacity of GIL conductors entirety
Error is more than setting error, i.e. GIL conductors operation is not up to stable state, then needing increases skin temperature is iterated calculating again, directly
Meet that condition is jumped out circulation and carried out in next step to error.
Step S6：The thermal convection current between conductor and shell and thermal exposure are calculated, the quantity of heat production and heat dissipation capacity for judging conductor miss
Whether difference is less than setting error；Step S7：If step S6 errors are less than setting error, calculate and terminate output as a result, otherwise leading
Temperature increases 0.1 DEG C of return to step S2.
GIL running temperature appraisal procedures in the corridor pipe of the present invention, by the MATLAB environment of Windows operating system
Extra-high voltage GIL thermal characteristics MATLAB calculation procedures in the corridor pipe of design, using the GIL relevant featuring parameters under initial operating mode come pair
Temperature during its stable operation carries out assessment calculating.
The above-mentioned initial duty parameter of operating mode include conductor loading electric current, corridor tube environment temperature, GIL pipeline bus internal pressure,
Corridor inner air tube flow velocity, GIL pipeline configuration parameters.
Give the initial parameter value of different operating modes, through can be calculated the joule heat waste of GIL conductors and shell, conductor with it is outer
The convection current heat output and thermal exposure between convection current heat output and thermal exposure, shell and air, the temperature of conductor and shell between shell
Degree.
Verify example：
For verify the above method correctness, using conductor internal diameter as 160mm, conductor wall thickness 20mm, internal diameter of outer cover 860mm,
The GIL of shell wall thickness 16mm conducts heat as experimental subjects for the convection current between the joule heat waste of conductor and shell, conductor and shell
The data such as amount and the convection current heat output between thermal exposure, shell and air and thermal exposure are calculated.Test the extraneous ring of GIL
Border temperature is 25 DEG C, and experiment condition is rated current 8000A.
GIL ambient temperatures are 25 DEG C, it is assumed that the initial temperature of GIL initial launch moment conductors is 24.9 DEG C, shell
Initial temperature be 25 DEG C, the joule heat waste of initial time can be calculated by formula 7~12 conductor and shell is respectively
187.96W/m, 80.47W/m, convection current heat output and thermal exposure between conductor and shell are respectively 2.19 × 10-3W/m、-2.48
×10-2W/m, convection current heat output and thermal exposure between shell and air are 0W/m.Can substantially it be observed by result of calculation, should
Moment, GIL was not stable.The temperature for needing to increase shell or conductor at this time carries out loop iteration calculating, until setting misses
Difference is met：
The joule heat waste of conductor and shell is respectively 177.88W/m, 41.15W/ when GIL stable operations are obtained by calculation
M, convection current heat output and thermal exposure between conductor and shell are respectively 93.77W/m, 75.39W/m, pair between shell and air
It is respectively 167.13W/m, 15.56W/m to spread heat and thermal exposure.By the conductor that is obtained during GIL stable operations and shell
Convection current heat output between convection current heat output and thermal exposure, shell and air and heat radiation between joule heat waste, conductor and shell
The temperature of conductor and shell is respectively 66.7 DEG C, 39.3 DEG C when amount can obtain stable operation.And in actual experiment, GIL is led
The temperature of body and shell is respectively 68.8 DEG C, 40.6 DEG C, and the conductor and the error of skin temperature obtained with this appraisal procedure is
3.2%.It can be seen that this calculating appraisal procedure has higher precision, it is practicable.
The conductor and skin temperature that the assessment of GIL running temperatures appraisal procedure is calculated in the corridor pipe of the present invention can characterize
Conductor and skin temperature during the actual normal operation of GIL busbares, embedded this assessment computation model in GIL temperature monitoring systems,
It can quickly find GIL temperature anomaly phenomenons, and GIL is overhauled, the operational reliability of GIL can be effectively improved.
GIL running temperature appraisal procedures in corridor pipe provided by the invention, assessment calculating is carried out by designing MATLAB programs,
For the program for predicting GIL thermal characteristics is embedded into GIL temperature monitoring systems the feasible means that provide.This method can facilitate
Efficiently try to achieve the heat-transfer character of GIL and the temperature of conductor and shell under different operating modes, and then to GIL from bringing into operation to steady
The thermal characteristics of timing has a research well, its calculation procedure realized, easy to operate, highly practical, calculating speed is fast, can
Initial parameter is set according to different operating modes, provides good theoretical foundation for the research to GIL thermal characteristics, and can make
It is embedded into for subprogram in the temperature monitoring system of GIL, the reliability of GIL, which plays an important role, is realized to raising.
GIL running temperature appraisal procedures in corridor pipe provided by the present invention are described in detail above, and are applied
Specific case is set forth the principle of the present invention and embodiment, and the explanation of above embodiment is only intended to help and understands
The method and its core concept of the present invention；Meanwhile for those of ordinary skill in the art, according to the thought of the present invention, then have
There will be changes in body embodiment and application range, in conclusion this specification content should not be construed to this hair
Bright limitation.
Claims (7)
1. GIL running temperature appraisal procedures in a kind of corridor pipe, it is characterised in that comprise the following steps：
Step S1：GIL physical models are established, model is related to resistance joule heat effect, the free convection laminar flow stream of outside air
Move or force turbulent convection, built-in electrical insulation gas (such as SF6Gas) electricity such as the Laminar Flow of free convection, radiant heat transfer, heat,
Stream, solid multiple physical field, are a multiple physical field coupling models.Therefore GIL running temperature parameter indexes in this model evaluation corridor pipe
Pair between convection current heat output and thermal exposure, shell and air between joule heat waste, conductor and shell including conductor and shell
Spread heat and thermal exposure.
Step S2：The loss of conductor initial power is calculated, according to above-mentioned physical model, is related to the fluid and solid material category of calculating
Property considers the change with local temperature, i.e. outside air, built-in electrical insulation gas, conductor and enclosure material thermal physical property parameter is
The function formula of temperature, the enclosure material thermal physical property parameter include viscosity, and the coefficient of conductivity is permanent than pressing, specific heat capacity,
Thermal coefficient of expansion, so that the numerical value of hot Physical Parameters is closer to operation actual conditions when calculating；The basis in actually solving
The thermal physical property parameter and skin effect coefficient of conductor, determine the initial temperature of conductor, it is assumed that initial temperature is equal with environment temperature,
Go out the heat of corresponding public condition lower conductor according to iterative solution under given initial duty parameter calculating GIL pipeline startup operating conditions
Amount；
Step S3：The loss of shell initial power is calculated, the collection skin effect of Shell resistance is considered when calculating Shell resistance joule heat amount
Should, i.e., on shell section, central current density is small, big by the current density of near surface, according to the physical model of step S1, true
After the thermal physical property parameter and skin effect coefficient of having determined shell, the initial temperature of conductor is determined, joined according to given initial operating mode
Number calculates the heat amount that iterative solution under GIL pipeline start operating performances goes out corresponding public condition lower conductor.
Step S4：Thermal convection current and the thermal exposure of shell and air are calculated, the quantity of heat production and heat dissipation capacity for judging GIL conductors entirety miss
Whether difference is less than program setting error, and error is no more than the quantity of heat production of 5%, GIL conductors entirety by the quantity of heat production and shell of conductor
Quantity of heat production composition；The heat dissipation of GIL conductors is carried out by the shell directly contacted with outside air, therefore GIL conductors entirety is scattered
Heat is made of the convection current heat output between shell and outside air and radiant heat transfer amount；Calculate the shell and outside air being related to
Thermal physical property parameter be temperature function, by judge GIL conductors entirety quantity of heat production and heat dissipation capacity error whether be less than set
Error is determined to judge whether the heat production of GIL conductors entirety and heat dissipation reach balance, i.e. whether GIL conductors operation has reached steady
State.
Step S5：If step S4 errors are less than setting error, next step is carried out, otherwise shell initial temperature increases 0.1 DEG C
Return to step S3；
Step S6：Calculating the thermal convection current between conductor and shell and thermal exposure, the quantity of heat production and heat dissipation capacity error for judging conductor is
It is no to be less than setting error；The direct representation of the power attenuation of conductor is the quantity of heat production of conductor, and the heat dissipation of conductor is mainly led to
Insulating gas progress is crossed, therefore the heat dissipation capacity of conductor is made of the free convection heat transfer amount between conductor and shell and radiant heat transfer amount,
The calculating of thermal convection current and thermal exposure between conductor and shell, relates generally to the running temperature of conductor and shell, and insulation gas
The thermal physical property parameters such as the dynamic viscosity of body, thermal conductivity factor, specific heat at constant pressure, coefficient of cubical expansion, wherein all hot physical property ginsengs
Number is the function of temperature.Judge that conductor produces by the way that whether the quantity of heat production and heat dissipation capacity error that judge conductor are less than setting error
Whether heat and heat dissipation reach balance, i.e. whether conductor operation has reached stable state.
Step S7：If step S6 errors are less than setting error, calculate and terminate output as a result, otherwise conductor temperature increases 0.1 DEG C
Return to step S2；When the quantity of heat production of conductor is less than setting error with heat dissipation capacity error, illustrate that the operation of conductor reaches stable state, can
The running temperature of conductor during obtaining stable operation；If the quantity of heat production of conductor is more than setting error, i.e. conductor with heat dissipation capacity error
Operation is not up to stable state, and need increases conductor temperature at this time is iterated calculating again, until error meets that condition can be jumped out and follows
Ring, close to the temperature of actual conductor and shell when output GIL is stable.
2. GIL running temperature appraisal procedures in a kind of corridor pipe according to claim 1, it is characterised in that the step S1
The unknown quantity that middle GIL physical models are related to includes the convection transfer rate of air and insulating gas, thermal diffusivity, thermal conductivity, dynamic
Power viscosity, specific heat capacity etc. are the thermal physical property parameter of temperature funtion, are the conductor of temperature funtion and the thermal conductivity factor sum aggregate skin of shell
Effect coefficient, by using the parameter contacted with temperature, can make result more closing to reality.
3. GIL running temperature appraisal procedures in a kind of corridor pipe according to claim 1, it is characterised in that the step S1
In initial duty parameter include conductor loading electric current, corridor tube environment temperature, GIL pipelines internal pressure, corridor inner air tube flow velocity, GIL
Pipeline configuration parameter.
4. GIL running temperature appraisal procedures in a kind of corridor pipe according to claim 1, it is characterised in that the step S2
In, conductor considers the kelvin effect of conductor there are kelvin effect when calculating conductor resistance joule heat amount；Determine that conductor is initial
After temperature, it can determine that conductor initial power is lost according to initial working condition, wherein correlation computations formula is：
<mrow>
<mfenced open = "" close = "}">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>f</mi>
<mi>c</mi>
</mrow>
</msub>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>c</mi>
<mn>20</mn>
</mrow>
</msub>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>c</mi>
<mn>20</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<mn>20</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<msub>
<mi>S</mi>
<mi>c</mi>
</msub>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>S</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>D</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mi>c</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>f</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>+</mo>
<mn>0.03</mn>
<msup>
<mrow>
<mo>{</mo>
<mfrac>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<mn>0.0016</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<mn>75</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mi>c</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mn>10</mn>
</mfrac>
<mo>}</mo>
</mrow>
<mn>3.75</mn>
</msup>
<mo>&times;</mo>
<msup>
<mrow>
<mo>{</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<mn>0.0016</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<mn>75</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mi>c</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mn>10</mn>
</mfrac>
<mo>}</mo>
</mrow>
<mn>1.5</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<msup>
<mi>I</mi>
<mn>2</mn>
</msup>
<msub>
<mi>R</mi>
<mi>c</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, RcFor the resistance of conductor, KfcFor the skin effect coefficient of conductor, ρc20For 20 DEG C when conductor dc resistivity, Tc
For the running temperature of conductor, αc20For the temperature-coefficient of electrical resistance of conductor, ScFor the cross-sectional area of conductor, DcFor conductor diameter, CcFor
Conductor wall thickness, PcFor every meter of loss of conductor, I is conductor rated current.
5. GIL running temperature appraisal procedures in a kind of corridor pipe according to claim 1, it is characterised in that the step S3
In, shell considers the kelvin effect of shell there are kelvin effect when calculating Shell resistance joule heat amount；Determine that shell is initial
After temperature, it can determine that shell initial power is lost according to initial working condition.Wherein correlation computations formula is：
<mrow>
<mfenced open = "" close = "}">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>f</mi>
<mi>t</mi>
</mrow>
</msub>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>t</mi>
<mn>20</mn>
</mrow>
</msub>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>t</mi>
<mn>20</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mn>20</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<msub>
<mi>S</mi>
<mi>t</mi>
</msub>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>S</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>D</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mi>t</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>f</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>+</mo>
<mn>0.03</mn>
<msup>
<mrow>
<mo>{</mo>
<mfrac>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<mn>0.0016</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mn>75</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mi>t</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mn>10</mn>
</mfrac>
<mo>}</mo>
</mrow>
<mn>3.75</mn>
</msup>
<mo>&times;</mo>
<mo>{</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<mn>0.0016</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mn>75</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mi>t</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mn>10</mn>
</mfrac>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<msup>
<mi>I</mi>
<mn>2</mn>
</msup>
<msub>
<mi>R</mi>
<mi>t</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, RtFor the resistance of shell, KftFor the skin effect coefficient of shell, ρtaFor 20 DEG C when shell dc resistivity, Tt
For the running temperature of shell, αtaFor the temperature-coefficient of electrical resistance of shell, StFor the cross-sectional area of shell, DtFor envelope outer diameter, CtTo be outer
Shell wall is thick, PtFor every meter of loss of shell, I is shell rated current.
6. GIL running temperature appraisal procedures in a kind of corridor pipe according to claim 1, it is characterised in that the step S4
In, the quantity of heat production of GIL conductors entirety is made of the quantity of heat production of conductor and the quantity of heat production of shell, and heat dissipation capacity is by shell and outside air
Between convection current heat output and radiant heat transfer amount composition, correlation computations formula is：
<mrow>
<mfenced open = "" close = "}">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>t</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&alpha;</mi>
<mi>a</mi>
</msub>
<msub>
<mi>&pi;D</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>t</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>5.67</mn>
<msub>
<mi>&epsiv;</mi>
<mi>t</mi>
</msub>
<msub>
<mi>&pi;D</mi>
<mi>t</mi>
</msub>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>273</mn>
<mo>+</mo>
<msub>
<mi>T</mi>
<mi>t</mi>
</msub>
</mrow>
<mn>100</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>273</mn>
<mo>+</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
</mrow>
<mn>100</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&alpha;</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>Nu&lambda;</mi>
<mi>a</mi>
</msub>
</mrow>
<msub>
<mi>D</mi>
<mi>t</mi>
</msub>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>N</mi>
<mi>u</mi>
<mo>=</mo>
<mn>0.36</mn>
<mo>+</mo>
<mn>0.0914</mn>
<msub>
<mi>Gr</mi>
<mi>a</mi>
</msub>
<mo>+</mo>
<mn>0.363</mn>
<msup>
<msub>
<mi>Gr</mi>
<mi>a</mi>
</msub>
<mfrac>
<mn>1</mn>
<mn>6</mn>
</mfrac>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Gr</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>g&beta;</mi>
<mi>a</mi>
</msub>
<msup>
<msub>
<mi>D</mi>
<mi>t</mi>
</msub>
<mn>3</mn>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msubsup>
<mi>&mu;</mi>
<mi>a</mi>
<mn>2</mn>
</msubsup>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&beta;</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>273</mn>
<mo>+</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>d</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>d</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mi>t</mi>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&epsiv;</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mn>0.85</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<msub>
<mi>P</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>t</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>t</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<msub>
<mi>P</mi>
<mi>t</mi>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, QtcFor every meter of heat loss through convection amount of shell, QtrFor every meter of heat loss through radiation amount of shell, αaFor the cross-ventilation coefficient of heat transfer,
TaFor ambient air temperature, εtFor inner surface of outer cover blackness, Nu is the nusselt number of outside air, λaFor the thermal conductivity factor of air,
GraFor the Prandtl number of air, g is acceleration of gravity, βaFor the coefficient of cubical expansion of air, μaFor the dynamic viscosity of air, Tda
For air setting temperature, Δ1For the quantity of heat production and heat dissipation capacity error of GIL conductors.
7. GIL running temperature appraisal procedures in a kind of corridor pipe according to claim 1, it is characterised in that the step S6
In, the power attenuation of conductor is the quantity of heat production of conductor, and the heat dissipation capacity of conductor is by the free convection heat transfer amount between conductor and shell
Formed with radiant heat transfer amount, correlation computations formula is：
<mrow>
<mfenced open = "" close = "}">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>c</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;&lambda;</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>ln</mi>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>t</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
</mfrac>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>5.67</mn>
<msub>
<mi>&epsiv;</mi>
<mi>e</mi>
</msub>
<msub>
<mi>&pi;D</mi>
<mi>c</mi>
</msub>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>273</mn>
<mo>+</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mrow>
<mn>100</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>273</mn>
<mo>+</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mrow>
<mn>100</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&epsiv;</mi>
<mi>e</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mfrac>
<mn>1</mn>
<msub>
<mi>&epsiv;</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>D</mi>
<mi>c</mi>
</msub>
<mrow>
<msub>
<mi>D</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>C</mi>
<mi>t</mi>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>&epsiv;</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&epsiv;</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mn>0.85</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&lambda;</mi>
<mi>e</mi>
</msub>
<mo>=</mo>
<mn>0.4</mn>
<mi>&lambda;</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>Gr</mi>
<mi>s</mi>
</msub>
<msub>
<mi>Pr</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>0.2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Pr</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
<mi>&mu;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Gr</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>g&beta;</mi>
<mi>s</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>t</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msup>
<mi>&mu;</mi>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mn>133.6</mn>
<mo>+</mo>
<mn>17.6</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>d</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>+</mo>
<mn>273</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>37.7</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mn>5</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>d</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>+</mo>
<mn>273</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&mu;</mi>
<mo>=</mo>
<mn>1.42</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>5</mn>
</mrow>
</msup>
<mo>+</mo>
<mn>17.6</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>8</mn>
</mrow>
</msup>
<msub>
<mi>T</mi>
<mrow>
<mi>d</mi>
<mi>s</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&lambda;</mi>
<mo>=</mo>
<mn>6.446</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>5</mn>
</mrow>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>d</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>+</mo>
<mn>273</mn>
<mo>)</mo>
</mrow>
<mn>0.942</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>d</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mi>t</mi>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&beta;</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>273</mn>
<mo>+</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>d</mi>
<mi>s</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>c</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>c</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, QccFor every meter of heat loss through convection amount of conductor, QcrFor every meter of heat loss through radiation amount of conductor, εcFor conductor outside surfaces blackness, εe
For the suitable blackness of full group object, λeFor equivalent heat conductivity, λ is insulating gas thermal conductivity factor, GrsFor insulating gas Grashof
Husband's number, PrsFor insulating gas Prandtl number, CpFor insulating gas specific heat at constant pressure, μ is insulating gas viscosity, βsInsulating gas
Coefficient of cubical expansion, RtiFor housing interior radius, RcoFor conductor outer radius, TdsFor the setting temperature of insulating gas, Δ2For conductor
Quantity of heat production and heat dissipation capacity error.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711149403.2A CN108038268B (en) | 2017-11-18 | 2017-11-18 | Method for evaluating GIL (general in-lane) running temperature in corridor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711149403.2A CN108038268B (en) | 2017-11-18 | 2017-11-18 | Method for evaluating GIL (general in-lane) running temperature in corridor |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108038268A true CN108038268A (en) | 2018-05-15 |
| CN108038268B CN108038268B (en) | 2020-06-02 |
Family
ID=62092757
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711149403.2A Expired - Fee Related CN108038268B (en) | 2017-11-18 | 2017-11-18 | Method for evaluating GIL (general in-lane) running temperature in corridor |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108038268B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111597748A (en) * | 2020-05-22 | 2020-08-28 | 中国矿业大学 | Method for realizing fault judgment based on GIL thermal characteristics |
| CN111950797A (en) * | 2020-08-21 | 2020-11-17 | 中国科学院合肥物质科学研究院 | Local temperature prediction method for high-power water-cooling bus with connector |
| CN112103923A (en) * | 2020-05-22 | 2020-12-18 | 中国矿业大学 | Method for realizing GIL protection based on thermal characteristic abnormity |
| CN113835948A (en) * | 2020-06-23 | 2021-12-24 | 华为技术有限公司 | Temperature detection method, temperature detection device and electronic equipment |
| CN111597748B (en) * | 2020-05-22 | 2024-05-17 | 中国矿业大学 | Method for realizing fault judgment based on GIL thermal characteristics |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104792258A (en) * | 2015-04-23 | 2015-07-22 | 华北电力大学 | Gas-insulated metal-enclosed power transmission line monitoring method and device |
| CN104965946A (en) * | 2015-06-26 | 2015-10-07 | 合肥工业大学 | Finite element calculation method for steady temperature field of high-speed railway bearing |
| CN106124928A (en) * | 2016-06-21 | 2016-11-16 | 中国电力科学研究院 | Extra-high-voltage alternating current mixed gas GIL line segment temperature monitoring method |
-
2017
- 2017-11-18 CN CN201711149403.2A patent/CN108038268B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104792258A (en) * | 2015-04-23 | 2015-07-22 | 华北电力大学 | Gas-insulated metal-enclosed power transmission line monitoring method and device |
| CN104965946A (en) * | 2015-06-26 | 2015-10-07 | 合肥工业大学 | Finite element calculation method for steady temperature field of high-speed railway bearing |
| CN106124928A (en) * | 2016-06-21 | 2016-11-16 | 中国电力科学研究院 | Extra-high-voltage alternating current mixed gas GIL line segment temperature monitoring method |
Non-Patent Citations (3)
| Title |
|---|
| XIAOWEN WU.ETC: ""Thermal Problems in Gas-Insulated Transmission lines solved by the Finite-Element Method"", 《2012 SIXTH INTERNATIONAL CONFERENCE ON ELECTROMAGNETIC FIELD PROBLEMS AND APPLICATIONS》 * |
| 张改杰等: ""户内气体绝缘输电线路温升数值计算模型的研究"", 《陕西电力》 * |
| 王健等: ""基于热力耦合分析的GIL热致伸缩特性及其影响因素"", 《高电压技术》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111597748A (en) * | 2020-05-22 | 2020-08-28 | 中国矿业大学 | Method for realizing fault judgment based on GIL thermal characteristics |
| CN112103923A (en) * | 2020-05-22 | 2020-12-18 | 中国矿业大学 | Method for realizing GIL protection based on thermal characteristic abnormity |
| CN112103923B (en) * | 2020-05-22 | 2021-07-06 | 中国矿业大学 | Method for realizing GIL protection based on thermal characteristic abnormity |
| CN111597748B (en) * | 2020-05-22 | 2024-05-17 | 中国矿业大学 | Method for realizing fault judgment based on GIL thermal characteristics |
| CN113835948A (en) * | 2020-06-23 | 2021-12-24 | 华为技术有限公司 | Temperature detection method, temperature detection device and electronic equipment |
| CN111950797A (en) * | 2020-08-21 | 2020-11-17 | 中国科学院合肥物质科学研究院 | Local temperature prediction method for high-power water-cooling bus with connector |
| CN111950797B (en) * | 2020-08-21 | 2023-03-10 | 中国科学院合肥物质科学研究院 | Local temperature prediction method for high-power water-cooling bus with connector |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108038268B (en) | 2020-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Allahbakhshi et al. | Heat analysis of the power transformer bushings using the finite element method | |
| CN108038268A (en) | GIL running temperature appraisal procedures in a kind of corridor pipe | |
| Sedaghat et al. | Enhanced thermal model of power cables installed in ducts for ampacity calculations | |
| Zhang et al. | Experiments and modeling of heat transfer in oil transformer winding with zigzag cooling ducts | |
| CN108052697A (en) | A kind of corridor pipe GIL three-dimensional temperature fields and breathing deformation calculation method based on workbench | |
| Han et al. | Estimation of the time-dependent convective boundary condition in a horizontal pipe with thermal stratification based on inverse heat conduction problem | |
| Xie et al. | Simulation and experimental analysis of three‐dimensional temperature distribution of±400‐kV converter transformer valve‐side resin impregnated paper bushing under high current | |
| Chi et al. | A multi-scale thermal-fluid coupling model for ONAN transformer considering entire circulating oil systems | |
| Lin et al. | Theoretical research on temperature field of power cable joint with FEM | |
| Declercq et al. | Accurate hot spot modeling in a power transformer leading to improved design and performance | |
| CN104408249A (en) | Method and system for determining thermal parameter of single-core cable conductor | |
| Xu et al. | Heat transfer characteristic research based on thermal network method in submersible motor | |
| CN109490354A (en) | A kind of electric compartment reduced scale testing stand building method based on the theory of similarity | |
| Boukrouche et al. | Mock-up study of the effect of wall distance on the thermal rating of power cables in ventilated tunnels | |
| Cong et al. | Research on undetected overheat fault of the GIS bus bar contacts based on infrared thermal imaging | |
| Zhao et al. | Real-time reconstruction of temperature field for cable joints based on inverse analysis | |
| Hrnjak et al. | Effect of return bend and entrance on heat transfer in thermally developing laminar flow in round pipes of some heat transfer fluids with high Prandtl numbers | |
| CN109635380A (en) | A kind of piping lane electric compartment fire barriers and current-carrying capacity of cable quantitative relationship determine method | |
| Tang et al. | Hotspot temperature monitoring of fully insulated busbar taped joint | |
| Lahey Jr et al. | ANALYSIS OF TRANSIENT CRITICAL HEAT FLUX. | |
| Xu et al. | Calculation of the equivalent thermal parameters of external environment for cable based on the real-time measurement data | |
| Zajec et al. | Experimental and Numerical Analysis of Heat Transfer on the Annular Test Section | |
| Berge et al. | Evaluation of vacuum insulation panels used in hybrid insulation district heating pipes | |
| Long et al. | Finite element model for temperature distribution calculation of HTS cable based on experiment data | |
| Han et al. | Sensitivity analysis about transient three-dimensional IHCP with multi-parameters in an elbow pipe with thermal stratification |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20190212 Address after: No. 28 Xianning West Road, Xianning stele area, Xi'an, Shaanxi Applicant after: XI'AN JIAOTONG University Applicant after: STATE GRID CORPORATION OF CHINA Applicant after: STATE GRID JIANGSU ELECTRIC POWER Co.,Ltd. Applicant after: STATE GRID JIANGSU ELECTRIC POWER Co.,Ltd. MAINTENANCE BRANCH Address before: No. 28 Xianning West Road, Xianning stele area, Xi'an, Shaanxi Applicant before: Xi'an Jiaotong University |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200602 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |