Добавил:
ac3402546@gmail.com Направление обучения: транспортировка нефти, газа и нефтепродуктов группа ВН (Вечерняя форма обучения) Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Термодинамика и теплопередача в технологических процессах нефтяной и газовой промылшенности

.pdf
Скачиваний:
134
Добавлен:
01.06.2021
Размер:
7.08 Mб
Скачать

Теплопередача в технологических процессах...

241

Графически распределение температуры внутри цилиндра представлено

на рисунке 2.50.

t

 

 

 

 

t0

 

 

 

 

 

α

 

 

 

α

 

 

 

tс

 

 

tс

 

tн

 

 

 

 

 

tн

 

 

 

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 2.50. Распределение температуры

 

 

 

 

r

r0

 

 

 

от источника теплоты

 

 

 

 

 

 

 

 

 

 

в цилиндре к омывающей жидкости [8]

 

 

 

 

 

Для нашего случая обычно известен линейный тепловой поток от электродвигателя

q

=

Q

=

Nэл (1− ηэл )

,

(2.302)

 

 

 

 

 

эл

 

где Q – количество теплоты, выделяемое источником в единицу времени, Вт.

Поэтому для определения объемной плотности источника теплоты воспользуемся формулой

 

 

 

 

 

qv

=

q

.

 

 

 

(2.303)

 

 

 

 

 

 

 

π r2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

С учетом (2.303)

перепишем

выражение

для

определения

температуры

в сечении цилиндра (2.301)

 

 

 

 

 

 

 

 

 

t = t

 

+

 

q

+

 

q

 

 

(r2

r2 ) .

(2.304)

 

н

2

π r α

 

4 π r2

λ

0

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

В соответствии с уравнением (2.304) получим выражения для определения температуры:

на поверхности цилиндра r = r0

tс = tн+

q

;

(2.305)

2 π r α

 

 

 

0

 

 

242

Часть 2

 

 

на оси цилиндра r = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

q

 

 

 

q

 

 

 

q

 

1

 

 

1

 

t

 

= t

 

+

 

 

+

 

 

=t

 

+

 

=t

 

+

 

 

 

 

+

 

 

. (2.306)

 

 

 

π r α

 

π λ

 

4 π λ

 

 

π

α r

 

 

 

0

 

н

2

4

 

с

 

 

н

 

2

 

2

λ

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

В последнем уравнении в скобках стоят два слагаемых, которые представляют собой линейные термические сопротивления передачи теплоты теплоотдачей от поверхности цилиндра (электродвигателя) к текучей среде и теплопроводностью в материале цилиндра. Чем больше коэффициент теплоотдачи, и коэффициент теплопроводности, тем меньше отличаются между собой температура жидкости, наружной и внутренней поверхности цилиндра.

Растворение газа сопровождается выделением примерно Qрг на всей длине

насоса и незначительным нагревом потока, так как по сравнению с теплотой, выделяемой насосом, величина теплоты растворения составляет примерно 5%.

Тепловой поток, выделяемый при работе трехжильного плоского медного кабеля (в месте подвески центробежного насоса) приводит к увеличению температуры потока только на десятые доли градуса.

Теплота от электрокабеля выделяется по всей длине скважины от электродвигателя до устья. Суммарный эффект от выделения теплоты от электрокабеля на отрезке скважины выше точки подвеса насоса составляет несколько градусов.

Уменьшение температуры на том же интервале скважины за счет теплопередачи между скважинной продукцией и горной породы оценивается величиной такого же порядка.

Нагнетательные скважины

Нагнетательные скважины используются для закачки в пласт воды и различных теплоносителей. Условная схема нагнетательной скважины представлена на рисунке 2.51. Закачка воды или теплоносителя ведется по колонне на- сосно-компрессорных труб.

В этом случае уравнение распределения температуры воды по стволу скважины принимает следующий вид:

t = θ0

+ Г h

Г

[1e−A h ]+ (t

у− θ0 ) e−A h ,

(2.307)

 

 

 

A

 

 

где tу – температура горячей воды на устье скважины; A – комплекс (Шухова), определяемый параметрами конструкции скважины и свойствами воды [8, 12].

Если в скважину нагнетается холодная вода, а для ее подогрева используется забойный нагреватель, то для определения температуры закачиваемой воды в пласт необходимо учитывать тепловой поток Qнагр , выделяемый нагревателем

Теплопередача в технологических процессах...

243

Рис. 2.5 1. Условная схема нагнетательной скважины:

1 – устье; 2 – колонна насосно-компрессорных труб; 3 – обсадная колонна; 4 – цементный камень; 5 – пакер; 6 – забой; 7 – горная порода

 

 

tз =

t +

 

Qнагр

 

.

 

 

 

 

 

 

(2.308)

 

 

αнагр

π dнагр

нагр

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где t , t

з

– температура воды до и после нагревателя, оС ;

α

нагр

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

 

 

 

 

 

 

 

 

 

 

 

 

теплоотдачи от поверхности нагревателя к воде, Вт/(м2 К) ;

нагр

, d

нагр

– длина

 

 

 

 

 

 

 

 

 

 

 

 

и диаметр нагревателя, м.

Определение темпе ратуры закачиваемой воды на забое скважины является необходимым условие м для подбора наземного теплосилового обору дования.

244

Часть 2

 

 

Паронагнетательная скважина

Водяной пар является высокоэффективным теплоносителем. При нагнетании в пласт влажного пара происходит конденсация и его температура практически не изменяется. С точки зрения закачки в пласт теплоносителя с высокой температурой это является положительным моментом. Но значительно усложняется работа конструкции скважины, необходима тепловая изоляция колонн для того, чтобы избежать значительных перепадов температуры на цементном камне.

Для обеспечения неизменной температуры насыщенного пара необходимо, чтобы потери теплоты не превышали теплоту конденсации r – фазового перехода. Если температура пара вдоль ствола скважины за время закачки τ не изменяется, то потери теплоты в скважине при нагнетании влажного или сухого насыщенного пара определяются на основании решения задачи нестационарного охлаждения полого неограниченного цилиндра [8]

Q = 2 π λэк H u (tstгп),

(2.309)

где λэк − эквивалентный коэффициент теплопроводности системы скважина – порода

 

 

 

 

ln

d

цк + 2δгп

 

 

 

 

λэк =

 

 

 

 

 

 

dнктв

 

,

(2.310)

 

 

 

 

 

 

 

1

ln di+1 + 1 ln dцк

+ 2δгп

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1 λi

 

di

 

λгп

dцк

 

λi − коэффициенты теплопроводности элементов крепи скважины (колонны на- сосно-компрессорных, обсадных труб); di и di+1 − внутренний и наружный

диаметры крепи скважины (d

1

= d в

);

δ

гп

= 2

 

a

гп

τ

гп

− толщина прогрева

 

 

 

 

 

 

 

 

 

нкт

 

 

 

 

 

 

 

 

 

горной породы за время; ts

− температура насыщенного пара;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

Fo

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 0,5 − 0,25

+ 0,125 Fo,

 

Fo ≤ 2

 

 

 

 

 

 

 

π

 

 

 

 

 

π Fo

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.311)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

,

Fo > 2,

 

 

 

 

 

 

 

 

 

(ln4 Fo −

2

2

 

 

 

 

 

ln4

Fo − 2 j

 

j)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fo =

4 aэк

τ

− число Фурье;

aэк − эквивалентный коэффициент температуро-

(dнктв )2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

проводности системы скважина – порода, рассчитываемый по эквивалентным значениям теплопроводности, теплоемкости и плотности системы скважина – порода.

Теплопередача в технологических процессах...

245

Если при закачке насыщенного пара происходит его полная конденсация на определенной глубине h , то дальнейшее изменение температуры горячей воды следует рассчитывать по соотношению (2.307) с учетом соответствующей глубины.

Расстояние от устья до глубины, где произойдет полная конденсация ( x = 0) можно рассчитать исходя из уравнения распределения степени сухости пара

 

x = xу

+ B h − A ts

h + t0

A h + A Г

(h )2

,

(2.312)

 

 

 

 

 

 

 

2

 

 

 

 

 

где xу

− степень сухости пара на устье скважины; A =

π dнктв

k

; B =

g

;

 

 

 

 

 

 

 

 

G r

 

 

r

G − массовый расход пара, кг/с;

Степень сухости пара по глубине скважины можно также оценить из баланса

теплоты

 

 

Q

.

(2.313)

x = xу G r

 

Для предварительной оценки потерь теплоты в скважине (Q) и степени сухости пара по глубине скважины на основании соотношений (2.309) – (2.313) приведены номограммы на рисунках 2.52 и 2.53.

По диаграмме на рисунке 2.52 в зависимости от времени закачки пара, которое входит в число Fo, эквивалентного коэффициента теплопроводности системы скважина – горная порода (λэк ), перепада температур между температурами

насыщенного пара (ts ) и горной породы (tгп) можно оценить потери теплоты по

глубине скважины ( H ).

Номограмма для оценки степени сухости пара по глубине скважины (рисунок 2.53) построена на основании формулы (2.313). В зависимости от давления насыщения пара, потерь теплоты, массового расхода пара и степени сухости пара на устье оценивается конечная степень сухости.

246

Часть 2

 

 

Рис. 2.52. Потери теплоты в паронагнетательной скважине

Теплопередача в технологических процессах...

247

Рис. 2.53. Определение степени сухости пара

Температурное поле бурящейся скважины

Распределение температур нисходящего и восходящих потоков бурового раствора особенно важно при бурении в зонах вечной мерзлоты для недопущения растепления ствола скважины (рисунок 2.54).

Расчетные соотношения по определению температуры восходящего и нисходящего потоков получаются из совместного решения уравнений первого начала термодинамики (2.252) и теплопередачи (2.253) с учетом соответствующих граничных условий.

248

Часть 2

 

 

Рис. 2.54. Условная схема бурящейся скважины: 1 – устье; 2 – нисходящий поток;

3 – колонна буровых труб; 4 – восходящий поток; 5 – долото; 6 – горная порода

Температура нисходящего потока в бурильной колонне

t

 

 

= t

 

+ Г h

Г

+ M

 

exp

A

2 − B

h

 

R

exp

A

2 + B

h

 

,

(2.314)

низ

0

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

Температура восходящего потока в затрубном пространстве

 

 

 

 

 

 

 

 

t

 

 

 

= t

 

+ Г h + M

 

 

exp

A2 − B

h

 

R

exp

A2 + B

h

 

,

(2.315)

 

 

 

верх

0

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k1 π d1

 

 

 

k2 π dск

 

 

 

 

 

 

 

 

 

где A

1

=

, A

2

=

, B =

A22 + 4 A1 A2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G

cpm

 

 

 

 

G cpm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теплопередача в технологических процессах...

249

Коэффициенты теплопередачи от нисходящего к восходящему потоку в колонне буровых труб k1 и от восходящего потока в окружающие горные породы

k2 определяются по соотношениям

k1

=

 

 

 

 

1

 

 

 

,

k2

=

 

 

 

α3

 

 

,

(2.316)

 

1

 

d1

 

d2

 

d1

 

 

 

 

 

 

 

 

 

 

 

 

+

ln

+

 

 

 

α3

dск

 

 

aгп τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α1

2 λст

 

d1

 

α2 d

2

 

 

1+

 

 

ln 1

+

2

 

 

 

 

 

 

 

 

 

 

 

4 λгп

 

 

dск

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где α1, α2 , α3 − коэффициенты теплоотдачи внутри бурильной колонны, от стенок труб к восходящему потоку и от восходящего потока к горным породам, соответственно; d1 , d2 , dск − внутренний, наружный диаметры труб бурильной колонны и диаметр скважины; τ − время работы, с.

Постоянные M1 , M2 , R1 и R2 определяются следующим образом:

 

 

 

 

 

 

 

 

 

 

 

Г

 

A2 + B

 

 

 

 

 

 

B H

 

 

 

 

 

A2 H

 

 

 

 

t

t

0

+

 

 

 

 

 

 

 

 

 

 

exp

 

 

 

 

 

 

 

+ Г exp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M =

 

 

 

 

 

 

 

 

A1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2

,

(2.317)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B H

 

 

 

 

 

 

B H

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A2

sh

 

 

 

 

 

 

 

 

+ B ch

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M2 = M1

B − A2

,

 

 

 

 

 

 

 

 

 

 

(2.318)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B + A

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г

 

A2 − B

 

 

 

 

 

 

 

 

 

B H

 

 

 

 

 

A2 H

 

 

 

t

t

 

+

 

 

 

 

 

 

 

 

 

 

 

 

exp

 

 

 

 

 

 

 

 

 

+ Г exp

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R =

 

 

 

 

 

 

 

 

A1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2

,

(2.319)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B H

 

 

 

 

 

B H

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A2 sh

 

 

 

 

 

 

+ B ch

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R2 = R1

B + A2

.

 

 

 

 

 

 

 

 

 

(2.320)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B − A

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для определения чисел подобия конвективного теплообмена при движении неньютоновской жидкости вводится понятие эффективной вязкости

 

 

 

 

τ

0

d

 

 

 

 

 

= 1

+

 

 

,

 

(2.321)

 

эф

 

 

 

 

 

 

 

 

 

 

 

 

6 c

 

 

где – динамический коэффициент вязкости, Па.с; τ

0

– динамическое напря-

 

 

 

 

 

 

 

 

 

жение сдвига, Па; d – линейный размер.

250

Часть 2

 

 

2.14. Температурный режим магистральных газонефтепроводов

Одним из первых уравнение температурного режима нефте – и газопроводов получил знаменитый русский инженер В.Г. Шухов, который предложил определять текущую температуру нефти по длине трубопровода при неизотермическом течении

t

x

= t

гр

+ (t t

гр

) exp(−a / L),

(2.322)

 

 

1

 

 

где t1 , tx , tгр – температуры нефти или газа на входе, на расстоянии x от начала участка и грунта на этом участке; L – длина трубопровода; a – параметр Шухова

а =

km π dн L

;

(2.323)

 

G cpm

 

km – коэффициент теплопередачи от транспортируемого продукта в окружающую среду на участке трубопровода; dн – наружный диаметр трубопровода;

G – массовый расход продукта; cpm – средняя массовая изобарная теплоемкость

продукта.

Уравнение Шухова (2.322), в первую очередь, может быть рекомендовано для определения температуры нефти в любом сечении x линейных участков нефтепроводов, так как изменение энтальпии жидкостей и, в частности нефти, с определенной степенью точности может определяться также как и для идеального газа из соотношения (1.80).

Расчет температурных режимов линейных участков газопроводов по уравнению Шухова дают существенную погрешность, так как оно не учитывает эффект изменения температуры реальных газов в процессе дросселирования. Этот эффект, названный эффектом Джоуля-Томсона, учитывается в уравнении по распределения температуры газа по длине участка газопровода, полученном В.И. Черникиным и С.А. Бобровским

tx

= tгр + (t1 tгр ) exp(a x / L)

Dh (p1 px ) [1exp(а x / L)]

, (2.324)

 

 

 

x / L)

 

 

 

 

где D

– коэффициент Джоуля-Томсона (1.288), (1.290), К/МПа; p ,

p

x

, – абсо-

h

1

 

 

лютное

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

от начала участка, МПа.

Наиболее точно уравнение распределения температуры по длине магистрального газо-нефте- и продуктопроводов получено Н.И. Белоконь на основе первого начала термодинамики по балансу рабочего тела для потока и уравнения теплопередачи. при учете свойств реального газа [1, 7].