книги / Теплопередача в скважинах
..pdfпринцип сохранения тепловой энергии и выводится следующим обра зом.
Выделим мысленно внутри тела элементарный объем, в котором равномерно распределены источники тепла (рис. 1), объемная мощ ность тепловыделения которых равна qv.
Приращение внутренней энергии вещества в выделенном объеме -составит
dQ = dQ l - d Q i, |
(1.5) |
где dQt — количество тепла, выделенное внутренними источниками,
dQ1 = qv dxdydzdx; |
(1.6) |
dQ 2 — количество тепла, ушедшее сквозь |
поверхность наружу. |
Для определения величины dQ2 рассмотрим направление по оси х . В этом направлении через левую грань поступает внутрь выде ленного объема dQx тепла, которое можно определить из уравне
ния (1.3) |
|
|
dQx ——X |
dy dz dx. |
(1.7) |
Через противоположную грань за тот же промежуток времени |
||
вытечет из объема dQx+dx тепла |
|
|
dQx+dx-----% [_ ^ + ~ ( ~ ) d x ~ \ d y d z d x = |
|
|
= ~ A'( '& + ~ £ r ^ ) dJ'dzdT- |
(I 8 ) |
|
Результативное количество проходящего тепла составит |
|
|
dQx+<u-dQx = - X - £ L d x d y d z d x . |
(1.9) |
Полное количество проходящего через элементарный объем тепла во всех трех направлениях составит
(u o )
Изменение внутренней энергии тела можно вычислить через тепло емкость и приращение температуры
dQ = ср dx dy dz dx. |
(1.11) |
Подставив выражения (1.6), (1.10) и (1.11) в уравнение (1.5), получим дифференциальное уравнение теплопроводности в прямо угольной системе координат
11
Очевидно, что уравнению (1.12) можно придать следующий вид
dt |
дЧ |
дЧ |
|
дт |
а (дх* |
~ду* + - & ) + ■ £ • |
<‘ <3> |
где |
|
|
|
|
а = |
А. |
(1Д4> |
|
|
ср |
|
Зная вблизи той или иной точки зависимость температуры от коор динат, можно предсказать, как быстро будет изменяться температура в этой точке с течением времени.
Уравнение (1.13) упрощается, когда отсутствуют в теле источники тепловыделения, т. е. qv = 0. Еще более простой вид оно нриобре-
тает для случая стационарной теплопроводности, когда /)/ = 0.
Предельно простой вид имеет дифференциальное уравнение, описы вающее стационарную одномерную теплопроводность
дЧ
0.
дх*
Часто задачи теплообмена в стволе скважины приходится решать в цилиндрической системе координат, для которой дифференциальное уравнение представляет собой следующее выражение:
dt___ / |
дЧ_<_____ L . J L - 4 - |
1 |
дЧ |
дх |
а \ дг2 “ Г г дг ' |
г2 |
* <?<р2 |
Ж \ . |
qv |
|
(U 5 ) |
|
dz2 ) ~ |
ср |
# |
||
|
где г — радиус-вектор; ср — полярный угол; z — аппликата. Физический смысл дифференциального уравнения теплопровод
ности (уравнения Фурье) заключается в том, что им связывается про странственное распределение температуры с изменением его во вре мени. При этом чем больше коэффициент температуропроводности, тем быстрее меняется во времени температура. Поэтому при прочих равных условиях выравнивание температуры во всех точках рассма триваемой системы будет происходить быстрее в том теле, которое обладает большим коэффициентом температуропроводности. Жидко сти и газы обладают большой тепловой инерционностью и, следова тельно, малым коэффициентом температуропроводности. Металлы обладают малой тепловой инерционностью, так как коэффициент температуропроводности у них велик.
$ 3. УРАВНЕНИЕ НЬЮТОНА — РИХМАНА. ДИФФЕРЕНЦИАЛЬНОЕ УРАВНЕНИЕ КОНВЕКТИВНОГО ТЕПЛООБМЕНА
Для описания процесса конвективного теплообмена между по верхностью тела и омывающей его средой используют закон Ньютона — Рихмана, согласно которому количество тепла, отдава емое единицей поверхности тела в единицу времени, пропорционально
12
перепаду температур между поверхностью тела и окружающей средой
|
|
q = aAt, |
(1.16) |
где а |
— коэффициент пропорциональности, называемый коэффициен |
||
том |
теплоотдачи, |
в Вт/м2*градус. |
интенсивность тепло |
Коэффициент |
теплоотдачи характеризует |
обмена между поверхностью тела и окружающей средой. Численно он равен количеству тепла, отдаваемого (или воспринимаемого) единицей поверхности в единицу времени при разности температур между поверхностью тела и окружающей средой, равной 1 градусу.
В отличие от коэффициента теплопроводности коэффициент тепло отдачи не является физической постоянной того или иного вещества. Он отражает совместное действие конвекции и излучения и зависит от многих факторов. Только конвективная составляющая коэффици ента теплоотдачи определяется геометрической формой и размерами тела, физическими свойствами омывающей тело среды, направле нием и скоростью омывания, температурными условиями и т. д. Вся сложность вопроса о теплообмене между телом и окружающей средой заключена в определении величины коэффициента теплоот дачи при конкретных условиях задачи. Практически познание про цесса теплоотдачи сводится к определению зависимости коэффициента теплоотдачи от различных факторов.
У поверхности омываемого жидкостью твердого тела всегда имеет место неподвижный слой жидкости, в котором перенос тепла осуще ствляется теплопроводностью в соответствии с законом Фурье
где п — нормаль к поверхности тела.
С другой стороны, согласно закону Ньютона — Рихмана это же
удельное количество тепла передается потоку |
жидкости |
q = a{tc- t M). |
(1.18) |
Приравнивая друг к другу правые части уравнений (1.17) и (1.18) г получаем дифференциальное уравнение теплообмена, которое описы вает процесс теплоотдачи на границах тела (п = 0):
g » - :— . ( - Щ |
. |
(1.19) |
tс—*ж Vдп )п*о |
|
4 |
При изменении температурного поля в жидкости только в напра влении одной координаты (одномерная задача), нормальной к поверх ности тела, уравнение (1.19) запишется в следующем виде
а = |
(1.20) |
13
Температурное поле в потоке движущейся жидкости описывается уравнением энергии
dt |
. |
dt . |
dt . |
dt |
дЧ \ |
I T + ^ lS T + ^ - S iГ + & -а Г в в |
( 1.21) |
||||
где ц;Л1 |
и?г — составляющие вектора скорости движения жидкости |
||||
в направлении осей х , у, 2. |
|
|
|||
Многочлен, стоящий в левой части уравнения (1.21), представляет |
|||||
собой полную производную от температуры по времени |
кото |
рую часто обозначают
Применив известное обозначение
дЧ дЧ дЧ
дх* *• дуг “Г дъ2 “ V Ч
уравнение энергии потока запишем в следующей форме:
( 1.22)
Если движение отсутствует, т. е. wx = wy = wz — 0, то уравне ние энергии переходит в дифференциальное уравнение теплопровод ности для случая отсутствия в теле внутренних источников тепла.
О J
Бри стационарных процессах конвективного теплообмена д% = 0.
В случае стационарного одномерного температурного поля все производные по т, у и %равны нулю.
Температурное поле в движущейся жидкости зависит от составля ющих скоростей wx, wy, wz. Поэтому уравнение теплоотдачи и энергии содержит пять неизвестных a, t, wx, wy, wz. Чтобы система уравнений оказалась замкнутой, необходимо к имеющимся двум добавить уравне ния, описывающие изменение скорости во времени и пространстве, т. е. уравнения движения.
Трудности определения коэффициента теплоотдачи аналитическим путем очевидны. Поэтому до настоящего времени используют экспе риментальный метод. Чисто экспериментальное определение вели чины а возможно, в принципе, тремя методами.
1. Реализуют измерения температуры в пристенном слое при у = 0. Зная дополнительно коэффициент теплопроводности среды, рдшвающей стенку, с помощью формулы (1.20) легко вычислить величину а. Однако осуществить температурные измерения в при стенном слое часто трудно, а иногда невозможно.
2.Используются сведения о градиенте температуры в твердом теле у самой его поверхности. При этом, зная коэффициент тепло проводности материала тела, по формуле (1.20) определяют а. Техни ческие трудности использования этого метода очень велики.
3.Наиболее распространенный метод заключается в непосред ственном измерении удельного теплового потока q и действующего
14
вдоль поверхности температурного напора At. Используя затем урав нение Ньютона — Рихмана, получают средние для поверхности коэф фициенты теплоотдачи. Приближение их к локальным коэффициен там осуществляется носредством секционирования рассматриваемой поверхности и раздельного определения средних значений а для каждой секции.
§ 4. КРАЕВЫЕ УСЛОВИЯ* ПОНЯТИЕ О КОЭФФИЦИЕНТЕ ТЕПЛОПЕРЕДАЧИ
Дифференциальное уравнение теплопроводности описывает явле ние переноса тепла в наиболее общем виде и справедливо для огром ного числа процессов. Оно относится к бесконечно малому элементу температурного поля и не характеризует развитие теплопроводности во всем пространстве и за весь период времени. Для получения кар тины, отражающей качественные и количественные признаки кон кретного случая, необходимо математически поставить конкретную задачу и найти частное решение основного уравнения теплопровод ности. Частные особенности, которые совместно с дифференциальным уравнением теплопроводности дают полное математическое описание конкретного процесса, называют краевыми условиями или условиями однозначности.
Краевые условия включают в себя:
геометрические условия, характеризующие форму и размеры тела, в котором протекает процесс;
физические условия, характеризующие физические свойства среды и тела;
начальные условия, характеризующие распределение темпера туры в изучаемой системе в начальный момент времени, предшеству ющий рассматриваемому процессу;
граничные условия, характеризующие взаимодействие рассматри ваемого тела с окружающей средой.
Граничные условия могут быть сформулированы одним иэ ука занных ниже способов.
1. Граничные условия первого рода заключаются в том, что за дается распределение температуры на поверхности тела для каждого момента времени
tc = f(x, У, Z, т).
В наиболее простом случае может оказаться, что температура на поверхности одинакова и с течением времени не меняется
tc= const.
2. Граничные условия второго рода состоят в том, что задаются величины теплового потока для каждой точки поверхности тела и любого момента времени
Япо, = /(*. У, г, Т).
15
Впростейшем случае плотность теплового потока по поверхности
ис течением времени остается неизменной
к |
Зпов = const. |
ч |
|
3. |
Граничные условия третьего рода предполагают, что в функ |
ции от времени известны температура среды, омывающей поверхность тела, и интенсивность проникновения потока тепла сквозь эту поверх ность, характеризующаяся коэффициентом теплоотдачи.
Согласно закону сохранения энергии количество тепла, которое отводится с единицы поверхности в единицу времени вследствие теплоотдачи, должно равняться теплу, подводимому к единице по верхности в единицу времени из внутренних зон тела вследствие теплопроводности. Таким образом, граничное условие третьего рода характеризуется дифференциальным уравнением конвективного теп лообмена
4. Граничные условия четвертого рода характеризуют условия теплообмена системы тел или тела с окружающей средой по закону теплопроводности. Предполагается, что между телами (телом и сре дой) существует совершенный контакт (температуры соприкаса ющихся поверхностей одинаковы), в результате чего бывает равен ство тепловых потоков, проходящих через поверхность соприкос новения,
л 8*1 _ л ^2 Л11 Г ~ “ Ла"йГв
При решении задач теплопередачи в стволе скважины часто при ходится использовать граничные условия третьего рода. Причем, теплообмен между восходящим и нисходящим потоками циркулиру ющей в скважине жидкости через колонну труб удобно характеризо вать комплексным показателем — коэффициентом теплопередачи.
Коэффициент теплопередачи от одной среды к другой через разде ляющую их стенку определяет Отнесенное к одному градусу темпе ратурного напора количество теплоты, которое передается в единицу времени от горячей среды к более холодной, считая на единицу по верхности разделяющей их стенки
к = ~ Ф т г - |
<U 3) |
где К — коэффициент теплопередачи через стенку, разделяющую среды с разной температурой; (£х — t 2) — перепад температуры между разделенными перегородкой средами.
* Индекс ш относится к телу, индекс «ж» — к среде*
16
Для случая теплопередачи через плоскую однослойную перего родку коэффициент теплопередачи определяется выражением
|
_ L + i L _ i _ '* |
(U 4) |
|
«1 ' Хс "I" 0t2 |
|
где |
<хх, а 2 — коэффициенты теплоотдачи на поверхностях |
стенки, |
fic, |
— соответственно толщина и коэффициент теплопроводности |
|
стенки. |
|
|
|
Поскольку К есть величина, обратная тепловому сопротивлению, |
ее иногда называют тепловой проводимостью системы «среда — пластина — среда». Так как общее сопротивление ПК превышает любое из своих составляющих, то коэффициент теплопередачи всегда меньше любого из коэффициентов теплоотдачи.
Для случая теплопередачи через однослойную трубу линейный
коэффициент теплопередачи K t можно определить из формулы |
|
|||||
1 |
1 |
1 |
<$2 |
1 |
(1.25) |
|
In |
||||||
|
||||||
|
2л Хе |
|
di |
Л^2^2 |
|
Достаточно тонкостенная труба может практически рассматри ваться как плоская стенка. При этом если отношение толщины трубы
ксреднему диаметру составляет 2 ч- 3%, то ошибка в определении
К, вносимая указанным допущением, составляет не более 3%.
Более толстостенные трубы можно привести к плоской стенке, если в качестве расчетного диаметра использовать полусумму вну треннего и внешнего диаметров
du = . d'± l* . |
(1.26) |
При этом ошибка расчета составляет 2—3% при 6!d1 — 30 -4- 40. Если последнее отношение больше 40, то коэффициент теплопередачи необходимо считать по точным формулам.
§ 5. КРИТЕРИАЛЬНЫЕ УРАВНЕНИЯ
Каким бы сложным не было математическое выражение, предста вляющее решение той или иной краевой задачи, все входящие в него размерные величины могут быть сгруппированы в безразмерные ком бинации — комплексы и симплексы, состоящие из однородных или разнородных физических величин. Уравнения, составленные с по мощью безразмерных комбинаций, имеют очевидные преимущества: количество аргументов значительно сокращается, а количественное соотношение между входящими в формулу величинами является универсальным, оставаясь постоянным для целого класса явлений.
В любой краевой задаче различают три вида величин: независи мые переменные, зависимые переменные и постоянные величины.
17
Зависимые переменные являются функцией независимых переменных и постоянных величин.
Для приведения уравнения к безразмерному виду необходимо в первую очередь выбрать масштабы приведения. В качестве масшта бов удобно использовать постоянные величины, входящие в краевые
условия. |
Поэтому в формулы вводят безразмерную температуру |
||
0 = —-—-, |
безразмерную координату X = -j—, безразмерную ско- |
||
ТТ 7 |
= |
W X |
И Т. Д. |
рость Wx |
|
Помимо безразмерных величин и координат, в уравнения также входят безразмерные комплексы, состоящие из разнородных физиче ских величин и называемые критериями подобия.
Так, критерий Био показывает, как велико отношение внутрен него теплового сопротивления тела к внешнему сопротивлению на границе тела с окружающей средой
где L — характерный размер тела (для плоской стенки L равно тол
щине, для трубы |
L равно диаметру, для кольцевого пространства |
|
L |
равно разности |
внешнего и внутреннего диаметров). |
и |
Критерий Фурье представляет собой независимую переменную |
|
является безразмерным временем |
ах
Fo = и
Если критерий Био совместно с безразмерными координатами характеризует температурное поле тела при стационарной теплоот даче с его поверхности, то число Фурье является характеристикой нестационарной теплопроводности.
Для характеристики теплообмена между омываемой поверхностью тела и омывающей средой используют критерий Нуссельта или крите рий теплоотдачи
NU = 4 ^ >
лж
гДе — теплопроводность омывающей среды.
В задачах конвективного теплообмена критерий Нуссельта обычно является искомой величиной, так как в него входит коэффициент теплоотдачи.
Критерий Рейнольдса характеризует соотношение силы инерции и сил вязкости
где v — кинематическая вязкость омывающей тело среды. Многочисленными исследованиями для круглой трубы устано
влено, что критерий Рейнольдса полностью определяет режим течения
18
жидкостей, подчиняющихся закону течения Ньютона. При значениях Re < 2300 наблюдается ламинарный режим течения, при Re > 2300 начинается переходный режим.
Глинистые растворы и некоторые высоколарафинистые нефти не подчиняются закону течения Ньютона. Для характеристики их тече ния необходимо учитывать не только отношение сил вязкости к силам инерции, но и отношение сил инерции к силам пластичности. Для практических расчетов Р. И. Шищенко предложил использовать критериальную связь вида
Eu = /(R e-),
где Ref — представляет собой обобщенный параметр Рейнольдса
который наиболее полно учитывает особенности потока вязко-пла стичной жидкости.
Многочисленными исследованиями установлено, что структурный режим течения вязко-пластичной жидкости в трубах отмечается при Re' < 2000, а в кольцевом канале — при Re' < 1600. Крити ческие значения обобщенного критерия Рейнольдса изменяются в широких пределах и зависят от реологических свойств жидкости. Это приводит к возникновению области, в которой нет однозначной зависимости коэффициента гидравлического сопротивления от Re. Однако в этой области существует зависимость критического значе ния обычного критерия Рейнольдса от параметра Хедстрема (Не)
где Цф — критическая скорость течения.
По свидетельству А. X. Мирзаджанзаде, при движении глинистого раствора в скважине возможны случаи ранней турбулентности, обус ловленные присутствием частиц твердого вещества.
Критерий Пекле характеризует собой отношение тепла, переноси мого конвекцией, к теплу, переносимому теплопроводностью,
а
Безразмерный комплекс, характеризующий подъемную силу, возникающую вследствие разности плотностей, называют критерием
Грасгофа
Gr _ figAtL*
Критерий Эйлера характеризует соотношение сил давления и сил инерции
2* |
19 |
Используя указанные выше критерии, дифференциальные урав нения теплоотдачи и энергии потока можно записать в безразмерной форме
Безразмерные величины являются новыми переменными: незави симые переменные — это безразмерные координаты; зависимые пере менные — это безразмерные температура и составляющие скорости* а также числа Nu и Ей. Постоянные величины Ре, Re, Gr задаются краевыми условиями.
Безразмерные формулы, выражающие независимую переменную в функции координат, времени и соответствующих критериев подо бия, называются критериальными.
Критериальные уравнения, характеризующие тепломассоиеренос* могут быть представлены в следующем виде:
6 = |
У, z, Bilf Bi2); |
|
е ~ / а(У, |
г , Z, Ре, Re, Gr); |
|
N u = fa(X, |
У, Z, Ре, Re, |
Gr); |
Eu = /4[(Z, |
У, Z, Pe, Re, |
Gr); |
W W 5(*, |
У, Z, Pe, Re, Gr). |
Комбинированием безразмерных величин можно получить новые безразмерные комплексы. Так, критерий Пекле, полученный при приведении к безразмерному виду уравнения энергии, можно пред ставить как произведение двух критериев
Ре = Re Рг |
wL |
v |
|
v |
а * |
||
|
Безразмерная величина Рг = via является новым комплексом* называемым критерием Прандтля. Критерий Рг показывает, как сильно отличается поле скоростей течения от поля температур* Если эпюра распределения скоростей совпадает с эпюрой распределе ния температур, то Рг = 1.
По аналогии с обобщенным параметром Рейнольдса Б. И. Есьманом и С. М. Кулиевым было предложено использовать в критериаль ных уравнениях для вязко-пластичных жидкостей обобщенные пара метры Грассгофа и Прандтля:
* ( * + -6т)и>£ г 0 Тоd
Рг' = " ( , + т £ г ) С-8
Как правило, экспериментальные связи характеристик тепломас сообмена ищут в виде критериальных уравнений. Это значительно упрощает инженерные расчеты всевозможных процессов в технике*
20