CN101900773B - Underground power cable current-carrying capacity online prediction system and method - Google Patents

Underground power cable current-carrying capacity online prediction system and method Download PDF

Info

Publication number
CN101900773B
CN101900773B CN2010102094768A CN201010209476A CN101900773B CN 101900773 B CN101900773 B CN 101900773B CN 2010102094768 A CN2010102094768 A CN 2010102094768A CN 201010209476 A CN201010209476 A CN 201010209476A CN 101900773 B CN101900773 B CN 101900773B
Authority
CN
China
Prior art keywords
power cable
temperature
soil
current
capacity
Prior art date
Application number
CN2010102094768A
Other languages
Chinese (zh)
Other versions
CN101900773A (en
Inventor
梁永春
Original Assignee
河北科技大学
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 河北科技大学 filed Critical 河北科技大学
Priority to CN2010102094768A priority Critical patent/CN101900773B/en
Publication of CN101900773A publication Critical patent/CN101900773A/en
Application granted granted Critical
Publication of CN101900773B publication Critical patent/CN101900773B/en

Links

Abstract

The invention relates to an environment factor monitoring and finite element-based underground power cable current-carrying capacity online prediction system and an environment factor monitoring and finite element-based underground power cable current-carrying capacity online prediction method. The system consists of a temperature sensor used for measuring the air temperature, a wind speed sensor, a humidity sensor, a singlechip system and a scheduling host of a power department. The method comprises the step of predicting the current-carrying capacity of an underground power cable by combining temperature, wind speed, humidity and other data measured in real time with the calculation of a finite element temperature field. The system and the method have the advantages that: the system and the method can accurately predict the current-carrying capacity of the underground power cable on line, so that the power department can develop the capacity of conveying electric energy to the greatest degree; meanwhile, the system and the method can avoid the resource waste phenomenon of wasting the talent on a petty job on the premise of guaranteeing safe and reliable operation.

Description

The underground power cable current-carrying capacity online prediction method
Technical field
The present invention relates to a kind of underground power cable current-carrying capacity online prediction method.
Background technology
The underground high-voltage power cable is the principal mode of city electric line; Because its installation cost is higher, the characteristics of easy-maintaining not; On the basis that guarantees its long-term safety reliability service; Rationally confirm the current-carrying capacity of high voltage power cable, the ability of bringing into play its transmission of electric energy to greatest extent has great importance for power department.
At present the method confirmed of underground power cable current-carrying capacity has two kinds, and a kind of is calculated off-line, promptly according to the soil property of locality, the highest temperature in hot month, calculates according to IEC60287 or finite element; Another kind is in line computation, through the temperature of on-line monitoring cable surface, calculates the cable core temperature through Re Lu or finite element analysis method, and estimates the current-carrying capacity of cable on this basis.
The real-time current-carrying capacity of cable is and closely-related physical quantitys such as the coefficient of heat conductivity of the construction of cable, surrounding soil, face of land air, and above-mentioned two kinds of methods are the current-carrying capacity of on-line prediction cable exactly all.First method is to calculate according to the situation of the highest temperature and given soil thermal conductivity to get, when temperature or soil thermal conductivity change, and the maximum load current that cable can bear, promptly current-carrying capacity also should adjust accordingly; For second method, the temperature through cable surface can calculate the cable core temperature under preload more accurately, and will predict that it is the influence that current-carrying capacity must be considered temperature, soil thermal conductivity that maximum can be born electric current.
At present, the computing method of underground power cable current-carrying capacity are main with FEM calculation mainly, at first suppose an electric current; Maximum temperature according to boundary condition calculating power cable insulation layer if be no more than insulation course long-term work tolerable temperature, then increases electric current; If surpass tolerable temperature; Then reduce electric current, loop iteration is till the cable insulation maximum temperature equals insulation course long-term work tolerable temperature, and the electric current of this moment is the current-carrying capacity of cable.
Summary of the invention
Technical matters to be solved by this invention provides a kind of underground power cable current-carrying capacity online prediction method.
The technical solution adopted for the present invention to solve the technical problems:
A kind of underground power cable current-carrying capacity online prediction method is characterized in that:
(1) at first build underground power cable current-carrying capacity online prediction system:
Said underground power cable current-carrying capacity online prediction system is made up of the dispatching host machine (6) of the temperature sensor that is used for the Measurement of Air temperature (1), air velocity transducer (2), humidity sensor (3), SCM system (4), power department; The output terminal of said temperature sensor (1), air velocity transducer (2), humidity sensor (3) connects the respective input of said SCM system respectively, and said SCM system is connected with said dispatching host machine (6) through GSM network (5); Said temperature sensor (1) and air velocity transducer (2) be installed in said power cable directly over, said humidity sensor (3) is installed in the said power cable soil along the line;
(2) concrete steps of said underground power cable current-carrying capacity online prediction method are following:
(1) modeling:
A. soil plow-in cable:
Power cable is 1m apart from ground, and border, power cable both sides is apart from cable 10m, and in the calculating, border, said both sides normal orientation hot-fluid rate of change is 0, and promptly the soil moisture no longer changes; Below the power cable to the deep subsoil border apart from cable 10m, the earth deep layer that fetches earth boundary temperature is constant, and temperature value is 383K for the deep subsoil temperature;
B. calandria cable:
The comb top is 1m far from ground, and power cable is laid in the comb, and border, comb both sides is apart from comb 10m, and in the calculating, border, said both sides normal orientation hot-fluid rate of change is 0, and promptly the soil moisture no longer changes; Below the comb to the deep subsoil border apart from comb 10m, the earth deep layer that fetches earth boundary temperature is constant, and temperature value is 383K for the deep subsoil temperature;
C. plough groove type cable:
The ditch groove depth is 1m apart from ground, and power cable is laid in the groove, and iron plate is often laid at indoor groove top; Cement plate is often laid at outdoor groove top, and border, model both sides is apart from groove 10m, in the calculating; Border, said both sides normal orientation hot-fluid rate of change is 0, and promptly the soil moisture no longer changes; Below the groove to the deep subsoil border apart from groove 10m, the earth deep layer that fetches earth boundary temperature is constant, and temperature value is 383K for the deep subsoil temperature;
(2) subdivision: with power cable and surrounding soil subdivision thereof is little unit;
Under direct burial, calandria and three kinds of modes of plough groove type, adopting triangle or quadrilateral units is little unit with the The model subdivision;
(3) read the following data of real-time measurement: air themperature T Air, face of land wind speed v Air, relative moisture of the soil rh; And the coefficient of heat conductivity λ of calculating face of land cross-ventilation coefficient of heat transfer α and soil;
A. utilize following formula (1) to calculate face of land cross-ventilation coefficient of heat transfer α:
α=7.371+6.43v air 0.75(1)
In the formula, α is the face of land cross-ventilation coefficient of heat transfer, the W/ (m of unit 2K); v AirBe the face of land wind speed that air velocity transducer is measured in real time, the m/s of unit;
B. utilize formula (2) to calculate soil thermal conductivity λ:
λ = ( 11.33 log rh - 3.4 ) 10 0.649 ρ d - 2 - - - ( 2 )
In the formula, λ is a soil thermal conductivity, the W/ of unit (mK); ρ dBe soil packing, units/m 3Rh is the relative moisture of the soil that humidity sensor is measured in real time;
(4) at first, choose initial current I1, the A of unit;
(5) call the finite element temperature calculation program, calculate the power cable insulation layer maximum temperature T1 under the electric current I 1, unit K;
(6) as T1 during less than power cable insulation layer long-term work insulation course tolerable temperature Tin, and Tin-T1>0.3K, then get I2=2I1, went to for (7) step; As T1 during greater than Tin, and T1-Tin>0.3K, then get I2=0.5I1, went to for (7) step; When | T1-Tin|≤0.5K, I1 is the power cable current-carrying capacity of asking, and withdraws from;
(7) call the finite element temperature calculation program, calculate the power cable insulation layer maximum temperature T2 under the electric current I 2;
(8) when | T2-Tin|>0.3K, went to for (9) step; When | T2-Tin|≤0.3K, I2 is the power cable current-carrying capacity of asking, and withdraws from;
(9) get I=I2+ (I2-I1) (Tin-T2)/(T2-T1), call the finite element temperature calculation program, calculate the power cable insulation layer maximum temperature T under the electric current I;
(10) when | T-Tin|>0.3K, get I1=I2, I2=I, T1=T2, T2=T goes to above-mentioned (9) step; When | T-Tin|≤0.3K, I is the power cable current-carrying capacity of asking, and withdraws from.
The invention has the beneficial effects as follows the current-carrying capacity of on-line prediction underground power cable more accurately; Power department can be brought into play the ability of its transmission of electric energy in view of the above to greatest extent; And can under the prerequisite that guarantees safe and reliable operation the wasting of resources phenomenon of low load with strong power can not appear.
Description of drawings
Fig. 1 is a hardware configuration of the present invention.
Fig. 2-1 is the synoptic diagram of soil direct burial underground power cable system of laying.
Fig. 2-2 is the synoptic diagram of calandria underground power cable system of laying.
Fig. 2-3 is the synoptic diagram of plough groove type underground power cable system of laying.
Fig. 3-1 is whole subdivision of the finite element of plow-in cable and surrounding soil synoptic diagram as a result.
Fig. 3-2 is the subdivision synoptic diagram as a result of the finite element cable section of plow-in cable and surrounding soil.
In Fig. 1,2,3, and 1 temperature sensor (model: DS18B20), 2 air velocity transducers (model: CQ2-FC-1), 3 humidity sensor (models: BQ8-CHR-01); 4 SCM systems (model of single-chip microcomputer: 89C52), 5GSM network, the dispatcher of 6 power departments (desk-top computer); 7 soil, 8 power cables, 9 combs; 10 cover plates, 11 back up pads, 12 grooves.
Embodiment
The concrete steps of present embodiment are:
(1) at first build underground power cable current-carrying capacity online prediction system:
Said underground power cable current-carrying capacity online prediction system is made up of the dispatching host machine (6) of the temperature sensor that is used for the Measurement of Air temperature (1), air velocity transducer (2), humidity sensor (3), SCM system (4), power department; The output terminal of said temperature sensor (1), air velocity transducer (2), humidity sensor (3) connects the respective input of said SCM system respectively, and said SCM system is connected with said dispatching host machine (6) through GSM network (5); Said temperature sensor (1) and air velocity transducer (2) be installed in said power cable directly over, said humidity sensor (3) is installed in the said power cable soil along the line;
(2) concrete steps of said power cable current-carrying capacity online prediction method are following:
(1) modeling:
A. soil plow-in cable:
Power cable is 1m apart from ground, and border, power cable both sides is apart from cable 10m, and in the calculating, border, said both sides normal orientation hot-fluid rate of change is 0, and promptly the soil moisture no longer changes; Below the power cable to the deep subsoil border apart from cable 10m, the earth deep layer that fetches earth boundary temperature is constant, and temperature value is 383K for the deep subsoil temperature;
B. calandria cable:
The comb top is 1m far from ground, and power cable is laid in the comb, and border, comb both sides is apart from comb 10m, and in the calculating, border, said both sides normal orientation hot-fluid rate of change is 0, and promptly the soil moisture no longer changes; Below the comb to the deep subsoil border apart from comb 10m, the earth deep layer that fetches earth boundary temperature is constant, and temperature value is 383K for the deep subsoil temperature;
C. plough groove type cable:
The ditch groove depth is 1m apart from ground, and power cable is laid in the groove, and iron plate is often laid at indoor groove top; Cement plate is often laid at outdoor groove top, and border, model both sides is apart from groove 10m, in the calculating; Border, said both sides normal orientation hot-fluid rate of change is 0, and promptly the soil moisture no longer changes; Below the groove to the deep subsoil border apart from groove 10m, the earth deep layer that fetches earth boundary temperature is constant, and temperature value is 383K for the deep subsoil temperature;
(2) subdivision: with power cable and surrounding soil subdivision thereof is little unit;
Under direct burial, calandria and three kinds of modes of plough groove type, adopting triangle or quadrilateral units is little unit with the The model subdivision; The cable area temperature variation is bigger, and subdivision density is higher, air-shed between air-shed and groove inner cable outside surface and trench wall between comb inner cable outside surface and comb inwall; Temperature variation is bigger, and needs calculation flow rate, and subdivision density is higher; The soil region temperature variation is less, and subdivision density is less;
(3) read the following data of real-time measurement: air themperature T Air, face of land wind speed v Air, relative moisture of the soil rh; And the coefficient of heat conductivity λ of calculating face of land cross-ventilation coefficient of heat transfer α and soil;
A. utilize following formula (1) to calculate face of land cross-ventilation coefficient of heat transfer α:
α=7.371+643v air 0.75 (1)
In the formula, α is the face of land cross-ventilation coefficient of heat transfer, the W/ (m of unit 2K); v AirBe the face of land wind speed that air velocity transducer is measured in real time, the m/s of unit;
B. utilize formula (2) to calculate soil thermal conductivity λ:
λ = ( 11.33 log rh - 3.4 ) 10 0.649 ρ d - 2 - - - ( 2 )
In the formula, λ is a soil thermal conductivity, the W/ of unit (mK); ρ dBe soil packing, units/m 3Rh is the relative moisture of the soil that humidity sensor is measured in real time.
(4) at first, choose initial current I1, the A of unit;
(5) call the finite element temperature calculation program, calculate the power cable insulation layer maximum temperature T1 under the electric current I 1, unit K;
(6) as T1 during less than power cable insulation layer long-term work insulation course tolerable temperature Tin, and Tin-T1>0.3K, then get I2=2I1, went to for (7) step; As T1 during greater than Tin, and T1-Tin>0.3K, then get I2=0.5I1, went to for (7) step; When | T1-Tin|≤0.5K, I1 is the power cable current-carrying capacity of asking, and withdraws from;
(7) call the finite element temperature calculation program, calculate the power cable insulation layer maximum temperature T2 under the electric current I 2;
(8) when | T2-Tin|>0.3K, went to for (9) step; When | T2-Tin|≤0.3K, I2 is the power cable current-carrying capacity of asking, and withdraws from;
(9) get I=I2+ (I2-I1) (Tin-T2)/(T2-T1), call the finite element temperature calculation program, calculate the power cable insulation layer maximum temperature T under the electric current I;
(10) when | T-Tin|>0.3K, get I1=I2, I2=I, T1=T2, T2=T goes to above-mentioned (9) step; When | T-Tin|≤0.3K, I is the power cable current-carrying capacity of asking, and withdraws from.
In the present embodiment, condition of convergence 0.3K actual and thermometric error decision by engineering.Temperature measurement accuracy is 0.1-0.2K, and in conjunction with the actual requirement of engineering, the condition of convergence is got 0.3K in the present embodiment.
With Shijiazhuang summer be example, actual measurement face of land air themperature is 35 degree, wind speed is 1.5m/s, soil moisture is 18.5%, calculating soil thermal conductivity thus is 0.8W/ (m 2K), the cross-ventilation coefficient of heat transfer in the face of land is 16W/ (m 2K), the single loop of soil direct-burried shown in Fig. 2-1 yi word pattern is arranged 800mm 2The YJLW02XLPE power cable single-end earthed, buries ground degree of depth 1000mm, cable spacing 200mm, and calculating current-carrying capacity of cable is 892A.
Comb lays 800mm shown in Fig. 2-2 2The calculating current-carrying capacity of YJLW02XLPE power cable is 858A.
Groove lays 800mm shown in Fig. 2-3 2The calculating current-carrying capacity of YJLW02XLPE power cable is 1085A.
The temperature field FEM calculation method that the present invention relates to is prior art (seeing above-mentioned background technology part), and is following in the face of temperature field FEM calculation method introduction down:
Because underground power cable length is approximately infinitely great with respect to xsect, thereby underground power cable crowd's steady temperature field belongs to two-dimentional steady heat conduction problem.
The temperature control equation that heat source region (like cable conductor, insulating medium, metal screen layer and armor) is arranged is as shown in the formula shown in (3).
λ ( ∂ 2 T ∂ x 2 + ∂ 2 T ∂ y 2 ) + q v = 0 - - - ( 3 )
In the formula, q vBe volume heat generation rate, W/m 3
The temperature control equation in no thermal source solid dielectric zone (like other layers of power cable, soil and comb etc.) is suc as formula shown in (4).
∂ 2 T ∂ x 2 + ∂ 2 T ∂ y 2 = 0 - - - ( 4 )
Lay with groove for comb and to lay under the dual mode, power cable and have air on every side between solid dielectric, there is the characteristic of the natural convection of being heated in air, and calculate in this part regional temperature field needs Solving Coupled continuity equation, the equation of momentum and energy equation.Respectively as shown in the formula shown in (5), (6), (7).
∂ u ∂ x + ∂ v ∂ y = 0 - - - ( 5 )
In the formula, u, v are the component of flow field velocity vector at x and y axle, and unit is m/s.
ρ ( u ∂ u ∂ x + v ∂ u ∂ y ) = - ∂ p ∂ x + η ( ∂ 2 u ∂ x 2 + ∂ 2 u ∂ y 2 ) + ρgα ( T - T r ) cos θ ρ ( u ∂ v ∂ x + v ∂ v ∂ y ) = - ∂ p ∂ y + η ( ∂ 2 v ∂ x 2 + ∂ 2 v ∂ y 2 ) + ρgα ( T - T r ) sin θ - - - ( 6 )
In the formula, T rBe the fluid reference temperature; α is a volume expansivity, and unit is K -1ρ is a fluid density, and unit is kg/m 3P is the pressure in flow field, and unit is Pa; η is a hydrodynamic force viscosity, and unit is Pas.
u ∂ T ∂ x + v ∂ T ∂ y - λ ( ∂ 2 T ∂ x 2 + ∂ 2 T ∂ y 2 ) = 0 - - - ( 7 )
Find the solution the temperature field and need confirm the border of whole temperature field, under the soil directly buried installation situation, suppose deep soil temperature constant (being the lower boundary among Fig. 2-1), shown in (8):
∂ 2 T ∂ x 2 + ∂ 2 T ∂ y 2 = 0 T ( x , y ) | Γ = T soil - - - ( 8 )
In the formula, T SoilFor the deep soil temperature, get 383K.
The normal direction gradient on border, both sides is 0 among Fig. 2-1, does not promptly receive the influence of power cable heating away from the temperature field of power cable regional soil, shown in (9):
∂ 2 T ∂ x 2 + ∂ 2 T ∂ y 2 = 0 λ ∂ T ∂ n | Γ = 0 - - - ( 9 )
Coboundary among Fig. 2-1 is the ground table boundary, and the heat in the soil is diffused in the air through the form of convection heat transfer, shown in (10):
∂ 2 T ∂ x 2 + ∂ 2 T ∂ y 2 = 0 - λ ∂ T ∂ n | Γ = α ( T - T air ) | Γ - - - ( 10 )
In following formula (8), (9) and (10), Γ is an integral boundary.
For soil direct-burried power cable, after whole field domain carries out subdivision, utilize the Jia Lvejin formula, can get the FEM calculation equation, shown in (11):
KT=P (11)
In the formula, K is the whole audience territory cell matrix relevant with coefficient of heat conductivity and cell configuration, and T is a whole audience territory cell temperature matrix, and P is the whole audience territory cell matrix relevant with heating.
Common whole field domain adopts little triangle to carry out subdivision, and the subdivision result is as shown in Figure 3.Temperature T on each triangular element is the linear function of coordinate x and y, is expressed as: T=a 1+ a 2X+a 3Y, a 1, a 2And a 3Determine by the temperature value on the Atria summit.Utilize this linear function that the variation formula of formula (3), (8), (9), (10) is prolonged triangular element and carry out integration, can get the k value and the p value of each unit, the last integrated following formula (11) that gets, its expansion form is:
k 11 k 12 · · · k 1 n k 21 k 22 · · · k 2 n · · · · · · · · · · · · k n 1 k n 2 · · · k nn T 1 T 2 · · · T n = p 1 p 2 · · · p 3 - - - ( 12 )
With Gaussian elimination following formula (12) is found the solution, can solve the temperature value of each node, thereby obtain the temperature field distribution of whole field domain and the maximum temperature of power cable.
Lay (seeing Fig. 2-2) and groove power cable laying (seeing Fig. 2-3) for comb, after whole field domain carries out subdivision, need calculate with alternative manner.
List the Contact Coupled boundary condition in solid dielectric zone and air fluid zone, the borderline condition that wherein is coupled need satisfy following two conditions:
It is continuous that coupling border w goes up temperature:
T w| 1=T w| 2 (13)
Heat flow density on the coupling border w is continuous:
q n| 1=q n| 2 (14)
Suppose the borderline Temperature Distribution of coupling; One of them zone 1 (like solid dielectric) is found the solution; Draw borderline local heat flux density of coupling and thermograde; Use formula (12) and formula (13) then and find the solution another zone 2 (like air dielectrics), to draw new Temperature Distribution on the coupling border.Distribute as the input in zone 1 with this, repeat aforementioned calculation up to convergence, the temperature field that provides whole field domain at last distributes, and obtains the maximum temperature of power cable insulation.
Lay with groove for comb and to lay two kinds of situation; All between power cable and outside soil, there is air layer; Its temperature field computation process comprises the The Coupling of solid and fluid, can adopt subregion to calculate, and the border coupling scheme are carried out the temperature field of whole field domain and calculated.
Subregion calculates, the implementation step of border coupling process is:
1. respectively the physical problem in power cable body, air, the soil is set up governing equation;
2. list each regional boundary condition, wherein solid area and air section be coupled borderline condition for temperature on the coupling border continuously and heat flow density continuous, promptly satisfy following formula (13) and formula (14).
3. electric power, air and soil are all carried out finite element solving (promptly whole by finding the solution with the mode of soil direct-burried) by conduction heat transfer, obtain the initial temperature field distribution.
4. according to the result of calculation in above-mentioned the 3rd step, find the solution convective-diffusive equation and radiation equation in the air-shed, to draw new Temperature Distribution on the coupling border.
5. calculate of the input of the boundary temperature of solid and fluid with above-mentioned the 4th step, find the solution the conduction governing equation, calculate the temperature field of power cable and soil region as power cable and soil.
6. repeat the above-mentioned 3-5 step, till satisfying the condition of convergence.The variable quantity of promptly working as temperature and fluid velocity is all less than 10 -10The time, can end iteration.Provide the final temperature field distribution and the maximum temperature of power cable.
From the above; As long as given air themperature, air wind speed, sun light intensity, soil moisture just can calculate face of land convection transfer rate, the coefficient of heat conductivity of soil; Thereby the temperature field that calculates whole field domain distributes, and finally calculates the maximum temperature of cable insulation.Then through alternative manner, with regard to the measurable current-carrying capacity that goes out the power cable under the current environment condition.

Claims (1)

1. underground power cable current-carrying capacity online prediction method is characterized in that:
(1) at first build underground power cable current-carrying capacity online prediction system:
Said underground power cable current-carrying capacity online prediction system is made up of the dispatching host machine (6) of the temperature sensor that is used for the Measurement of Air temperature (1), air velocity transducer (2), humidity sensor (3), SCM system (4), power department; The output terminal of said temperature sensor (1), air velocity transducer (2), humidity sensor (3) connects the respective input of said SCM system respectively, and said SCM system is connected with said dispatching host machine (6) through GSM network (5); Said temperature sensor (1) and air velocity transducer (2) be installed in said power cable directly over, said humidity sensor (3) is installed in the said power cable soil along the line;
(2) concrete steps of said underground power cable current-carrying capacity online prediction method are following:
(1) modeling:
A. soil plow-in cable:
Power cable is 1m apart from ground, and border, power cable both sides is apart from cable 10m, and in the calculating, border, said both sides normal orientation hot-fluid rate of change is 0, and promptly the soil moisture no longer changes; Below the power cable to the deep subsoil border apart from cable 10m, the earth deep layer that fetches earth boundary temperature is constant, and temperature value is 383K for the deep subsoil temperature;
B. calandria cable:
The comb top is 1m far from ground, and power cable is laid in the comb, and border, comb both sides is apart from comb 10m, and in the calculating, border, said both sides normal orientation hot-fluid rate of change is 0, and promptly the soil moisture no longer changes; Below the comb to the deep subsoil border apart from comb 10m, the earth deep layer that fetches earth boundary temperature is constant, and temperature value is 383K for the deep subsoil temperature;
C. plough groove type cable:
The ditch groove depth is 1m apart from ground, and power cable is laid in the groove, and iron plate is often laid at indoor groove top; Cement plate is often laid at outdoor groove top, and border, model both sides is apart from groove 10m, in the calculating; Border, said both sides normal orientation hot-fluid rate of change is 0, and promptly the soil moisture no longer changes; Below the groove to the deep subsoil border apart from groove 10m, the earth deep layer that fetches earth boundary temperature is constant, and temperature value is 383K for the deep subsoil temperature;
(2) subdivision: with power cable and surrounding soil subdivision thereof is little unit;
Under direct burial, calandria and three kinds of modes of plough groove type, adopting triangle or quadrilateral units is little unit with the The model subdivision;
(3) read the following data of real-time measurement: air themperature T Air, face of land wind speed v Air, relative moisture of the soil rh; And the coefficient of heat conductivity λ of calculating face of land cross-ventilation coefficient of heat transfer α and soil;
A. utilize following formula (1) to calculate face of land cross-ventilation coefficient of heat transfer α:
α=7.371+6.43v air 0.75(1)
In the formula, α is the face of land cross-ventilation coefficient of heat transfer, the W/ (m of unit 2K); v AirBe the face of land wind speed that air velocity transducer is measured in real time, the m/s of unit;
B. utilize formula (2) to calculate soil thermal conductivity λ:
λ = ( 11.33 log rh - 3.4 ) 10 0.649 ρ d - 2 - - - ( 2 )
In the formula, λ is a soil thermal conductivity, the W/ of unit (mK); ρ dBe soil packing, units/m 3Rh is the relative moisture of the soil that humidity sensor is measured in real time;
(4) at first, choose initial current I1, the A of unit;
(5) call the finite element temperature calculation program, calculate the power cable insulation layer maximum temperature T1 under the electric current I 1, unit K;
(6) as T1 during less than power cable insulation layer long-term work insulation course tolerable temperature Tin, and Tin-T1>0.3K, then get I2=2I1, went to for (7) step; As T1 during greater than Tin, and T1-Tin>0.3K, then get I2=0.5I1, went to for (7) step; When | T1-Tin|≤0.5K, I1 is the power cable current-carrying capacity of asking, and withdraws from;
(7) call the finite element temperature calculation program, calculate the power cable insulation layer maximum temperature T2 under the electric current I 2;
(8) when | T2-Tin|>0.3K, went to for (9) step; When | T2-Tin|≤0.3K, I2 is the power cable current-carrying capacity of asking, and withdraws from;
(9) get I=I2+ (I2-I1) (Tin-T2)/(T2-T1), call the finite element temperature calculation program, calculate the power cable insulation layer maximum temperature T under the electric current I;
(10) when | T-Tin|>0.3K, get I1=I2, I2=I, T1=T2, T2=T goes to above-mentioned (9) step; When | T-Tin|≤0.3K, I is the power cable current-carrying capacity of asking, and withdraws from.
CN2010102094768A 2010-06-25 2010-06-25 Underground power cable current-carrying capacity online prediction system and method CN101900773B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010102094768A CN101900773B (en) 2010-06-25 2010-06-25 Underground power cable current-carrying capacity online prediction system and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010102094768A CN101900773B (en) 2010-06-25 2010-06-25 Underground power cable current-carrying capacity online prediction system and method

Publications (2)

Publication Number Publication Date
CN101900773A CN101900773A (en) 2010-12-01
CN101900773B true CN101900773B (en) 2012-09-12

Family

ID=43226474

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010102094768A CN101900773B (en) 2010-06-25 2010-06-25 Underground power cable current-carrying capacity online prediction system and method

Country Status (1)

Country Link
CN (1) CN101900773B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102141591A (en) * 2011-01-06 2011-08-03 上海交通大学 Device and method for monitoring transient state current-carrying capability of power transmission line
CN102323496A (en) * 2011-05-31 2012-01-18 广东电网公司广州供电局 System for monitoring current-carrying capacity of cable conductor and method
CN102768321B (en) * 2012-08-08 2018-04-24 上海森首光电科技有限公司 A kind of power cable current-carrying capacity monitoring method
CN102880747A (en) * 2012-09-06 2013-01-16 上海海事大学 Method for modeling, calculating and analyzing temperature fields of photoelectric composite submarine cable
CN103616588B (en) * 2013-11-27 2016-06-08 国网浙江省电力公司舟山供电公司 The defining method in a kind of undersea cable current capacity and temperature field
CN104330659B (en) * 2014-10-21 2017-02-15 国网上海市电力公司 Quasi dynamic compatibilization method based on cable heat transmission model
CN105093047B (en) * 2015-09-17 2018-09-04 国家电网公司 Power circuit load monitoring early warning system
CN107315860A (en) * 2017-06-07 2017-11-03 国网湖北省电力公司武汉供电公司 Temperature field finite element method when many cross-linked polyethylene insulated cables are total to ditch

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101051747A (en) * 2007-05-18 2007-10-10 华为技术有限公司 Method, device and apparatus for multiple output distribution
CN201075123Y (en) * 2007-09-10 2008-06-18 中国电力科学研究院 Device for measuring high potential current
CN201311461Y (en) * 2008-11-28 2009-09-16 皇明太阳能集团有限公司 Grid-connection photovoltaic data detecting system
CN101614775A (en) * 2009-07-15 2009-12-30 河北科技大学 Transformer State Assessment system and appraisal procedure thereof based on Multi-source Information Fusion

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101051747A (en) * 2007-05-18 2007-10-10 华为技术有限公司 Method, device and apparatus for multiple output distribution
CN201075123Y (en) * 2007-09-10 2008-06-18 中国电力科学研究院 Device for measuring high potential current
CN201311461Y (en) * 2008-11-28 2009-09-16 皇明太阳能集团有限公司 Grid-connection photovoltaic data detecting system
CN101614775A (en) * 2009-07-15 2009-12-30 河北科技大学 Transformer State Assessment system and appraisal procedure thereof based on Multi-source Information Fusion

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
梁永春,李彦明,柴进爱等.地下电缆群稳态温度场和载流量计算新方法.《电工技术学报》.2007,第22卷(第8期), *
梁永春,柴进爱,李彦明等.基于FEM的直埋电缆载流量与外部环境关系的计算.《电工电能新技术》.2007,第26卷(第4期), *
梁永春,王忠杰,刘建业.排管敷设电缆群温度场和载流量数值计算.《高电压技术》.2010,第36卷(第3期), *

Also Published As

Publication number Publication date
CN101900773A (en) 2010-12-01

Similar Documents

Publication Publication Date Title
Ozudogru et al. 3D numerical modeling of vertical geothermal heat exchangers
Holmberg et al. Thermal evaluation of coaxial deep borehole heat exchangers
Qian et al. Modeling the interactions between the performance of ground source heat pumps and soil temperature variations
Aresti et al. A review of the design aspects of ground heat exchangers
Serageldin et al. Earth-Air Heat Exchanger thermal performance in Egyptian conditions: Experimental results, mathematical model, and Computational Fluid Dynamics simulation
Wu et al. Experimental measurement and numerical simulation of horizontal-coupled slinky ground source heat exchangers
Raymond et al. Numerical analysis of thermal response tests with a groundwater flow and heat transfer model
You et al. In-situ experimental study of heat exchange capacity of CFG pile geothermal exchangers
Li et al. Development of a numerical model for the simulation of vertical U-tube ground heat exchangers
Ozel Thermal performance and optimum insulation thickness of building walls with different structure materials
Zanchini et al. Long-term performance of large borehole heat exchanger fields with unbalanced seasonal loads and groundwater flow
Ozel Effect of wall orientation on the optimum insulation thickness by using a dynamic method
Gholamalizadeh et al. A comprehensive approach to design and improve a solar chimney power plant: A special case–Kerman project
Signorelli et al. Numerical evaluation of thermal response tests
Cui et al. Numerical analysis and experimental validation of heat transfer in ground heat exchangers in alternative operation modes
Daouas A study on optimum insulation thickness in walls and energy savings in Tunisian buildings based on analytical calculation of cooling and heating transmission loads
Pu et al. Simulation study on the thermal performance of vertical U-tube heat exchangers for ground source heat pump system
Shang et al. Analysis of geo-temperature recovery under intermittent operation of ground-source heat pump
Misra et al. CFD analysis based parametric study of derating factor for Earth Air Tunnel Heat Exchanger
Ozgener et al. Experimental prediction of total thermal resistance of a closed loop EAHE for greenhouse cooling system
Balbay et al. Temperature distributions in pavement and bridge slabs heated by using vertical ground-source heat pump systems
Lazzari et al. Long-term performance of BHE (borehole heat exchanger) fields with negligible groundwater movement
Gustafsson et al. CFD-modelling of natural convection in a groundwater-filled borehole heat exchanger
Cacabelos et al. Calibrated simulation of a public library HVAC system with a ground-source heat pump and a radiant floor using TRNSYS and GenOpt
Li et al. Numerical simulation of horizontal spiral-coil ground source heat pump system: Sensitivity analysis and operation characteristics

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
EXPY Termination of patent right or utility model
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120912

Termination date: 20150625