ФИЗИЧЕСКИЕ СВОЙСТВА НЕФТЕЙ
.pdfvk.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 10−4 м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 10−4 + 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 |
|
|||||||||||||||
1−exp |
− |
|
|
|
, |
|||||||||||||||||
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
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−κ(T−60) ,
где учтено, что ν(60)=15 сСт. Второе условие ν(20)= 40
сСт дает для коэффициента κ уравнение: 40 =15 e−κ(20−60),
из которого находим κ 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 e−0,02452 (10−60) 51,1сСт,
λнар. = |
|
0,3164 |
0,0274 |
, |
|
4 1,3 0,7 (51,1 10−6 ) |
|||||
|
|
|
|||
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(−ke−m ) ,
λэф. = 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 |
( |
−ke−m |
, где |
|||
|
|
|
||||||||||
|
|
|
4 v d νнар. |
|
m |
|
|
|
|
) |
||
|
|
|
|
|
|
|
|
|
|
|
||
k = |
κ |
(T0 −Tн ар. ), m = |
πKd L = |
4K L . |
|
|||||||
|
|
4 |
|
|
|
|
|
ρQCv |
ρvCvd |
|
Вданном случае имеем:
νнар. =12 e−0,04 (10−50) 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 |
10−6 ) |
|||
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 10−6 ) |
||||||||
|
|
|
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 10−6 ) |
||||||
|
|
|
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 в конце участка трубопровода: