Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

ФИЗИЧЕСКИЕ СВОЙСТВА НЕФТЕЙ

.pdf
Скачиваний:
117
Добавлен:
20.08.2019
Размер:
2.32 Mб
Скачать

vk.com/club152685050 | vk.com/id446425943

265

T

T = 9,81 0,004 400000

4 0С.

L

0

1950

 

 

 

 

152. Согласно формуле (92), учитывающей распределение температуры нефти по длине трубопровода с учетом выделяющегося тепла внутреннего трения, T(x)= const.= T0

только в том случае, если T0 Tнар. T = 0 или

T =

ρgQ i0

= T T

.

 

π K d

0

нар.

 

Отсюда находим, что коэффициент K должен удовлетво-

рять условию:

 

 

 

K =

ρgQ i

,

 

0

 

π (T0 Tнар. ) d

 

где i0 гидравлический уклон, равный в данном случае 0,003. Таким образом, имеем:

K = 890 9,81 22003600 0,003 0,29 Вт/м2 0С. 3,14 25 0,7

153. Коэффициент K теплопередачи от жидкости, движущейся в трубопроводе, к окружающему грунту определяется в данном случае формулой (94):

1

 

1

 

1

ln

D

 

1

ln

Dнар.

1

 

 

 

=

 

+

 

 

+

 

 

+

 

,

 

K D

α1 d

2λст.

d

2λиз.

D

α2 Dнар.

 

 

 

 

 

 

 

 

в которой d = D 2δ = 0,700

м – внутренний диаметр тру-

бы; Dнар.

= D + 2δиз.

= 0,736 м – диаметр трубы с изоляцией,

а коэффициент

α2

теплопередачи через грунт в окружаю-

щий воздух рассчитывается по формуле (95) Форхгеймера:

α2

Dнар. =

 

 

2λгр.

 

 

.

 

2H

Dнар. + (2H Dнар. )

2

 

 

 

ln

 

1

 

 

 

 

 

 

 

 

 

vk.com/club152685050 | vk.com/id446425943

266

Последовательно делая вычисления, имеем:

а. α2 Dнар. =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 1,2

 

 

 

 

1,0

Вт/м 0С;

 

ln

2 2

0,736 + (2 2 0,736)2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b. α d =100 0,7 = 70 Вт/м 0С;

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с.

 

 

1

 

ln

D

=

 

 

1

 

 

ln

 

 

0,72

3,52 104 м0С/ Вт;

 

 

2λст.

 

 

 

 

2 40

 

 

0,70

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

Dнар.

 

 

 

 

1

 

 

 

0,736

 

0

 

 

 

 

 

d.

 

 

 

 

 

ln

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

ln

 

 

0,022 м

С/ Вт;

 

 

2λиз.

 

 

 

D

 

 

2 0,5

0,720

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и, наконец:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

=

 

1

 

+3,52 104 + 0,022 +1 =1,0366 м0С/ Вт,

 

K 0,720

 

70

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда находим: K 1,34 Вт/м2 0С.

 

 

 

 

 

 

 

154. В данном случае можно использовать формулы:

а) Для случая отсутствия тепловой изоляции:

 

 

 

1

 

=

 

1

 

 

 

ln

D1

+

 

 

1

, где D

= D + 2δ

из.

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K1 D 2λиз.

 

 

 

D α2 D1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б) Для случая дополнительной тепловой изоляции:

1

=

1

ln

D1

+

1

ln

D

+

1

, где D

=D

2δ .

 

 

 

 

 

 

K2 D 2λиз.

 

D 2λ

 

D1

 

α2 D

 

1

 

 

 

 

 

 

 

Подставляя в каждую из этих формул численные значения параметров из условия задачи, получаем:

а)

 

 

1

=

 

1

ln

1,036

+

1

,

K

1

1,02

2

0,7

 

1,02

2 1,036

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда находим: K1 1,986 Вт/м2 0С.

б) Поскольку K2 = 0,5 K1 = 0,993 Вт/м2 0С, то имеем:

1

=

 

1

ln

1,036

+

1

 

ln

 

D

+

1

,

 

 

 

 

 

 

 

 

 

0,993 1,02

2

0,7

 

1,02

2 0,1

1,036

 

 

 

 

 

 

2D

откуда находим: D 1,443 м δ =0,5 (D D1)=0,054 м или 54 мм.

vk.com/club152685050 | vk.com/id446425943

267

155. По формуле (91) В.Г. Шухова находим коэффициент K0 теплопередачи от нефти, движущейся по надземно-

му участку трубопровода, в окружающую среду при отсутствии тепловой изоляции. Подставив в эту формулу численные значения параметров, получим:

 

 

 

 

 

3,14 0,7 K0 2000

 

 

32 40

= (20

40) 1

exp

 

 

,

1800 3600 850 1950

 

 

 

 

 

 

 

откуда находим K0 26,98 Вт/м2 0С.

Для того чтобы температура нефти даже в самом критическом случае не снизилась более чем на 1 0С, необходи-

мо, чтобы коэффициент K0

удовлетворял условию

 

39 40 (20 40)

 

 

 

 

 

3,14 0,7 K 2000

 

1exp

 

 

 

,

1800 3600 850 1950

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда K 3,17 Вт/м2 0С. Принимаем K =K =3,17 Вт/м2 0С.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 max

 

Определим теперь толщину δ

слоя тепловой изоляции

( D = D + 2δ ). Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

1

=

1

 

 

+

1

 

ln

D

 

или

 

 

 

 

 

 

K D

K

D

 

 

2λ

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

=

1

 

 

+

 

1

 

ln

D

.

 

 

3,17 0,72

 

 

 

 

 

 

 

2 0,2

0,72

 

 

 

 

26,98 D

 

 

 

 

Отсюда находим: D 0,85 м. Наконец, определяем δ :

δ = 0,5 (0,85 0,72)= 0,065 м или 65 мм.

156. Суммарный коэффициент K теплопередачи в принятых допущениях определяется формулой:

1

=

1

ln

D + 2δиз.

+

1

,

 

2λиз.

 

α2 (D + 2δиз. )

KD

 

D

 

vk.com/club152685050 | vk.com/id446425943

268

в которой D = 0,72 м, D + 2δиз. = 0,82 м, λиз. = 0,25 Вт/м 0С.

Таким образом, в этой формуле неизвестен только коэффициентα2 .

Для определения коэффициента α2 теплопередачи от

внешней поверхности слоя изоляции через грунт в воздух используем приближенную формулу (96):

 

α2 (D + 2δиз. ) =

 

 

2λгр.

 

 

 

 

 

ln

 

2H

 

 

+

 

λгр.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

+ 2δиз.

α0H

 

 

 

 

 

 

 

 

 

 

 

из которой вычисляем:

2 1,8

 

 

 

 

 

 

α2 (D + 2δиз. ) =

 

 

 

 

 

 

2,7 Вт/м 0С.

ln

2 1,3

+

 

 

1,8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,82

8 1,3

 

 

 

 

 

 

 

 

 

 

 

После этого вычисляем значение коэффициента K :

1

=

 

1

ln

 

0,82

+

 

1

 

 

0,63 K 2,2 Вт/м2 0С.

 

K 0,72

2

0,25

 

0,72

 

 

 

 

 

 

 

 

2,7

 

 

 

 

 

 

Приравнивая удельные Qπ тепловые потоки от нефти

в окружающую среду и от нефти к внешней поверхности слоя изоляции

Qπw =KD (Tнеф. Tвозд. )D (Tнеф. Tиз. )=1,54 (40 0) 61,6 Вт/м2,

находим температуру Tиз. внешней поверхности слоя изо-

ляции:

Tиз. = Tнеф. KDαD (Tнеф. Tвозд. ),

где α − коэффициент теплопередачи через слой изоляции. Далее имеем:

1

=

1

ln

D + 2δиз.

=

1

ln

0,82

0,26 ,

αD

2λиз.

 

2 0,25

0,72

 

 

D

 

 

αD = 3,84 Вт/ м2 0С.

Таким образом:

T (0)= 23,5 17,1 ln

vk.com/club152685050 | vk.com/id446425943

269

Tиз. = 40 1,3,84584 (40 0) 23,5 0С.

Используя, наконец, приближенную формулу (97) для распределения температуры в грунте при стационарном

тепловом режиме

Qw

 

 

r

 

T (y,z)= T

 

ln

,

 

 

 

из.

 

2πλ

гр.

r

 

 

 

0

 

гдеr = y2 +(z + H)2 ; r

= D 2 радиус трубы, а ее ось име-

 

 

 

0

 

 

 

 

ет координату (0,H). Вертикальная ось OZ системы отсче-

та проходит через центр трубы, горизонтальная ось OY рас-

полагается на поверхности грунта, причем точка О находит-

ся над центром трубопровода на расстоянии H над ним.

Распределение T(y)

температуры поверхности грунта

(плоскости z = 0 ) имеет вид:

T (y)= T

Qw π ln

y2 + H2

.

 

из.

2λгр.

D 2

 

Подставив сюда численные значения параметров, получим:

T (y)= 23,5

61,6

ln

y2 +1,69

или

2 1,8

0,36

 

 

 

y2 +1,69 0С. 0,36

Максимальная температура на поверхности грунта достигается в точке y = 0 , находящейся непосредственно над

осью нефтепровода. Полагая y = 0 , находим:

T (0)= 23,5 17,1 ln

02 +1,69

+1,54 0С.

0,36

 

 

157. При решении предыдущей задачи была получена формула для распределения температуры поверхности грунта над местом залегания нефтепровода:

vk.com/club152685050 | vk.com/id446425943

 

270

 

T(0) = 23,5 17,1 ln

y2 +1,69

0С.

0,36

 

 

Примем, что на границах образовавшейся полосы, отделяющих покрытую снегом землю от непокрытой, температура грунта равна 0 0С, тогда для определения полуширины h образовавшейся незаслуженной полосы получаем уравнение:

23,5

17,1 ln

h2

+1,69

= 0 .

0,36

 

 

 

Решив это уравнение, найдем: h 0,334 м или 2h 0,67 м. 158. Найдем сначала по формуле (89) зависимость ν(T)

вязкости нефти от температуры. Имеем:

ν(T)=15 e−κ(T60) ,

где учтено, что ν(60)=15 сСт. Второе условие ν(20)= 40

сСт дает для коэффициента κ уравнение: 40 =15 e−κ(2060),

из которого находим κ 0,02452 1/ 0 С. Поскольку скорость v перекачки известна:

v = 4 18003600 1,30 м/с, 3,14 0,72

то используя формулу (91) В.Г. Шухова, можно вычислить средний по участку коэффициент K теплопередачи:

 

4 K 135000

 

 

 

25 =10 +(60 10) e

860 1,3 0,7 1950

,

 

откуда находим:

 

 

 

 

 

 

K = −

1,3 0,7 860 1950

ln

25

10

3,4 Вт/м2 0 С.

 

 

10

 

4 135000

60

 

Потери напора на рассматриваемом участке вычисляем по модифицированной формуле Дарси-Вейсбаха (99):

vk.com/club152685050 | vk.com/id446425943

271

hτ = λэф. L v2 , d 2g

где λэф. эффективный коэффициент гидравлического со-

противления, учитывающий переменность температуры по длине участка.

Сначала вычисляем вспомогательные величины:

νнар. =15 e0,02452 (1060) 51,1сСт,

λнар. =

 

0,3164

0,0274

,

4 1,3 0,7 (51,1 106 )

 

 

 

k = κ 4 (T0

Tнар. )= 0,02452 4 (60 10) 0,307 ,

m=4K L(ρvdCv )=4 3,4135000(860 1,3 0,7 1950) 1,203 k exp(m)= 0,307 exp(1,203) 0,092 ,

Ei(0,307) 0,889 ; Ei(0,092) 1899, .

После этого находим эффективный коэффициент λэф. :

λэф. = λнар. m1 Ei(k)Ei(kem ) ,

λэф. = 0,0274 1,2031 [0,889 (1,899)] 0,023 .

Наконец, рассчитываем потери hτ напора:

hτ = 0,023

135000

 

1,32

 

382 м.

0,7

2 9,81

 

 

 

Интересно отметить, что если бы нефть имела начальную температуру 60 0С на всем протяжении участка, то потери hτ напора составили бы 335 м, что на 47 м меньше,

чем в действительности.

159. Потери напора на рассматриваемом участке вычисляем по модифицированной формуле Дарси-Вейсбаха

(99):

vk.com/club152685050 | vk.com/id446425943

272

hτ = λэф. L v2 , d 2g

в которой λэф. эффективный коэффициент гидравличе-

ского сопротивления, учитывающий переменность температуры потока по длине участка:

λ

эф.

=

0,3164

 

1

Ei(k )Ei

(

kem

, где

 

 

 

 

 

 

4 v d νнар.

 

m

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

k =

κ

(T0 Tн ар. ), m =

πKd L =

4K L .

 

 

 

4

 

 

 

 

 

ρQCv

ρvCvd

 

Вданном случае имеем:

νнар. =12 e0,04 (1050) 59,44 ,

k = κ4 (T0 Tн ар. )= 0,044 (60 10)= 0,5 .

По условию задачи требуется определить потери напора для убывающей последовательности расходов: 1000, 800 и 600 м3/ч или соответствующих им скоростей 1,339; 1,071 и 0,804 м/с.

1) Пусть v1 =1339,

м/с. Тогда имеем:

λнар. =

 

0,3164

 

0,0305 ;

 

0,514 (59,44

106 )

4 1,339

 

k = 0,5 ;

m=4K L(ρvdCv ) = 4 3,5 140000(900 1,339 0,514 2000) 1,582; k exp(m)= 0,5 exp(1582, ) 0103,

Ei(0,5) 0,560 ; Ei(0103, ) 1796, .

После этого находим эффективный коэффициент λэф. :

λэф. = 0,0305 1,5821 [0,560 (1,796)] 0,0238 .

Наконец, рассчитываем потери hτ напора:

vk.com/club152685050 | vk.com/id446425943

 

 

 

 

 

 

 

273

 

 

 

hτ = 0,0238

 

140000

 

1,3392

592

м.

0,514

 

2

9,81

 

 

 

 

 

 

 

 

2) Пусть v2 =1071,

м/с. Тогда имеем:

λнар. =

 

 

0,3164

 

 

 

0,032 ; k = 0,5 ;

4 1,071 0,514

 

(59,44 106 )

 

 

 

m = 4 3,5 140000(900 1,071 0,514 2000) 1,978 ; k exp(m)= 0,5 exp(1978, ) 0,0692 ,

Ei(0,5) 0,560 ; Ei(0,0692) 216, .

После этого находим эффективный коэффициент λэф. :

λэф. = 0,032 1,9781 [0,560 (2,16)] 0,02588.

Наконец, рассчитываем потери hτ напора:

hτ = 0,02588

140000

 

1,0712

412 м.

 

 

 

0,514

 

2 9,81

 

 

 

3) Пусть v1 = 0,804 м/с. Тогда имеем:

 

λнар. =

 

0,3164

 

0,0346

; k = 0,5 ;

4 0,804 0,514

(

59,44 106 )

 

 

 

m = 4 3,5 140000(900 0,804 0,514 2000) 2,635 ; k exp(m)= 0,5 exp(2,635) 0,036,

Ei(0,5) 0,560 ; Ei(0,036) 2,8.

После этого находим эффективный коэффициент λэф. :

λэф. = 0,0346 2,6351 [0,560 (2,8)] 0,0294 .

Наконец, рассчитываем потери hτ напора

hτ = 0,0294

 

140000

 

0,8042

263,8

м

0,514

2

9,81

 

 

 

 

 

и температуру TL в конце участка трубопровода:

vk.com/club152685050 | vk.com/id446425943

274

TL =10 +(60 10) exp(2,635) 13,6 0С.

160. Найдем сначала по формуле (89) зависимость ν(T)

вязкости нефти от температуры. Имеем:

ν(T)= 5 e−κ (T50) ,

где учтено, что ν(50)= 5 сСт. Второе условие ν(20)= 40

сСт дает для коэффициента κ уравнение: 40 = 5 e−κ (2050) ,

из которого находим κ 0,0693 1/ 0 С.

Далее составляем уравнение баланса напоров:

2

273 0,125 104 Q2

= λ

эф.

120000

 

v2

.

 

 

 

 

0,7

 

2 9,81

 

 

 

 

 

В нем использован эффективный коэффициент λэф. гидрав-

лического сопротивления, учитывающий переменность температуры по длине участка, см. формулы (99):

λ

эф.

=

0,3164

 

1

Ei(k)Ei

(

kem

, где

 

 

 

 

 

 

4 v d νн ар.

 

m

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

k =

κ

(T0 Tн ар. ), m =

πKd L =

4K L

.

 

4

 

 

 

 

 

 

 

 

 

ρQCv

ρvCvd

 

Подставляя в уравнение баланса напоров выражение для расхода Q через скорость v перекачки

Q = 3,14 0,72 v 3600 , 4

а также учитывая другие данные условия, получаем уравнение:

546 = v2 (8737,4 λэф. + 47,94). (*)

Полученное уравнение решаем методом итераций (последовательных приближений).

1-е приближение. В качестве 1-го приближения поло-

жим λэф.(1)

= 0,02 . Тогда из уравнения (*) находим скорость