Моделирование установившихся процессов теплообмена в аппаратах воздушного охлаждения газа

Размещено на http://www.allbest.ru/

Самарский государственный технический университет

МОДЕЛИРОВАНИЕ УСТАНОВИВШИХСЯ ПРОЦЕССОВ ТЕПЛООБМЕНА В АППАРАТАХ ВОЗДУШНОГО ОХЛАЖДЕНИЯ ГАЗА

С.В. Алимов, И.А. Данилушкин, В.Н. Мосин



Аннотация

Предложена модель процесса теплообмена в аппаратах воздушного охлаждения газа в установившемся режиме. Модель учитывает температурное распределение потока газа, движущегося по трубам теплообменного аппарата, а также изменения температурного поля воздуха после охлаждения очередного слоя труб теплообменника.

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



Основная часть

На газовых промыслах и магистральных газопроводах для охлаждения природного газа и смазочных масел используются аппараты воздушного охлаждения (АВО). Необходимость охлаждения газа на компрессорных станциях продиктована требованиями снижения эксплуатационных расходов в системе транспортировки.

В то же время необходимо иметь в виду, что если затраты энергии в системе привода компрессора снижаются по мере снижения температуры газа, то энергозатраты в аппаратах воздушного охлаждения возрастают по мере увеличения перепада температуры газа между входом и выходом. Поэтому важно определить оптимальное соотношение затрат на охлаждение газа с помощью АВО, с одной стороны, и на прокачку газа через АВО - с другой.

Аппараты воздушного охлаждения включают следующие элементы (рис.1): теплообменные секции, вентиляторы с электроприводами, несущие конструкции и системы управления. Теплообменная секция состоит из нескольких рядов оребренных труб, расположенных в шахматном порядке. Охлаждение газа осуществляется поперечным потоком воздуха. Установка охлаждения газа состоит из нескольких параллельно соединённых однотипных АВО. Поддержание заданной температуры газа на выходе установки осуществляется путём включения в работу необходимого количества аппаратов воздушного охлаждения или выводом из активного теплообмена части поверхностей охлаждения, изменением частоты вращения вентиляторов и др. Энергозатраты каждого АВО зависят от температур охлаждаемого газа, охлаждающего воздуха и их расходов.

Для выбора эффективного алгоритма распределения нагрузки между аппаратами воздушного охлаждения газа, обеспечивающего поддержание заданной температуры газа на выходе установки при минимальных энергозатратах, необходимо решить ряд задач, включающих задачи математического моделирования процессов теплообмена в системе АВО, определения статических и динамических характеристик объекта и синтеза на основе полученных результатов системы автоматического управления.

воздушный охлаждение газ теплообмен

Р и с. 1

АВО газа различаются по ряду конструктивных решений, которые были приняты при их создании. В то же время, несмотря на различия АВО, в основе конструкции большинства аппаратов лежит общий принцип - охлаждение газа путем обтекания пучков параллельных оребренных труб поперечным потоком воздуха [1]. Это позволяет создать единую математическую модель процесса охлаждения газа, приняв ряд допущений, аналогично [2]:

1) скорость газа во всех трубах теплообменного аппарата одинакова;

2) скорость потока воздуха по всему сечению теплообменного аппарата одинакова;

3) физические свойства сред неизменны;

4) скорость теплопередачи по всей длине аппарата пропорциональна разности первых степеней температур сред;

5) смешение газа в направлении движения среды незначительно и не учитывается;

6) смешение газа в направлении, перпендикулярном движению среды, полное;

7) cтенку трубы можно считать теплотехнически тонким телом.

В качестве базовой конструкции АВО газа рассматривается аппарат из шести рядов труб в секции с однонаправленным проходом охлаждаемого газа по трубам, для которого известны: R - внутренний радиус несущей трубы, L - длина оребрённых труб, - коэффициент увеличения поверхности аппарата.

С учётом принятых допущений поведение температуры газа в любой из труб АВО будет описываться уравнением теплового баланса:

, , , (1)

,

где - распределение температуры потока газа по длине трубы, - распределение температуры трубы по длине, v - скорость потока газа, - коэффициент поверхностного теплообмена между внутренней стенкой трубы и газом, c - теплоёмкость газа, - плотность газа.

В стационарном режиме и уравнение (1) принимает вид

, , . (2)

Для элементарного объёма трубы в стационарном режиме соблюдается равенство потоков тепла от газа к трубе и от трубы к воздуху. Считая температуру трубы в любой точке сечения одинаковой, можно записать:

, (3)

где S, - площади соприкосновения газа с внутренней стенкой трубы и воздуха с оребрением трубы соответственно, - коэффициент поверхностного теплообмена между оребрением трубы и воздухом, - температура воздуха. С учётом соотношения

(4)

из (3) можно выразить распределение температуры трубы по длине:

. (5)

Тогда с учётом (5) из уравнения (2) можно исключить температуру трубы:

.

Обозначив , запишем:

. (6)

Решение уравнения (6) позволит получить температурное распределение для потока газа, если считать, что температура воздуха, обдувающего каждую из оребрённых труб АВО, одинакова. Естественно, такое предположение является довольно грубым и позволяет получить лишь приближённую (качественную) оценку поведения температурного поля потока газа по длине теплообменника. Повышение точности модели может быть достигнуто за счёт учёта изменения температуры воздуха после прохождения каждого ряда оребрённых труб теплообменного аппарата.

Для объёмов газа и воздуха, прошедших в единицу времени через теплообменный аппарат, можно составить уравнение теплового баланса:

, (7)

где c, , V - удельная теплоёмкость, удельная плотность, рассматриваемый объём для газа; , , - удельная теплоёмкость, удельная плотность, рассматриваемый объём для воздуха; , - температура газа на входе и выходе; , - температура воздуха на входе и на выходе. Если рассмотреть участок трубы длиной (рис. 2), то объём газа, прошедшего через этот участок в единицу времени, составит величину

. (8)

Объём воздуха, «обменявшегося» теплом с объёмом газа V, будет рассчитываться по формуле

, (9)

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



Р и с. 2

Подставив (8) и (9) в уравнение (7) и выразив температуру воздуха на выходе , запишем:

. (10)

Устремив в (10) к нулю, через предельный переход получаем выражение для температуры воздуха в любой точке по длине трубы после теплообмена:

, (11)

где .

Теперь для каждого ряда труб из уравнения (6) может быть получено температурное распределение потока газа, которое затем будет использовано для расчёта изменения температурного поля потока воздуха после прохождения ряда труб. Рассчитанное температурное поле потока воздуха станет «входящим» для следующего ряда труб и так далее.

Для расчёта температурного распределения потока газа удобно воспользоваться математическим аппаратом структурной теории распределённых систем [3, 4]. Температурное распределение потока газа будет рассчитываться в результате пространственного интегрирования передаточной функции распределённого объекта, описывающегося уравнением (6) и стандартизирующей функции, которая должна формироваться для каждого ряда отдельно с учётом нового выражения для температурного поля потока воздуха на входе по выражению (11). Передаточная функция для выражения (6) принимает вид [4]

, (12)

а стандартизирующая -

. (13)

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

. (14)

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

в (13) и используя (14), получаем:

; (15)

. (16)

Продифференцировав , получаем:

. (17)

Подставив (17) в (11), получим температурное поле потока воздуха после первого ряда труб:

. (18)

Теперь можно повторить всю процедуру расчёта для второго ряда труб. Стандартизирующая функция (13) с учётом (18) примет вид

(19)

.

Как показано в (19), функцию можно представить в виде двух слагаемых. Результат пространственного интегрирования для передаточной функции (12) и уже известен - выражение (16). Проинтегрировав произведение передаточной функции (12) с функцией из (17), получаем:

. (20)

Тогда выражение для температурного распределения потока газа во втором ряду труб примет вид

. (21)

Производная , необходимая для вычисления температурного поля потока воздуха после второго ряда труб, согласно (11) запишется как

. (22)

Введем обозначение :

. (23)

Тогда выражение (11) для температурного поля потока воздуха после второго ряда труб с учётом (17) будет иметь вид

. (24)

Согласно (18) получаем:

. (25)

Для третьего ряда труб стандартизирующая функция (13), по аналогии с (19), с учётом (25) и (13) примет вид:

. (26)

Пространственное интегрирование передаточной функции (12) с функцией из (23) имеет вид

(27)

.

С учётом (27) можно записать выражение для температурного распределения потока газа в третьем ряду труб:

, (28)

откуда

, (29)

где

. (30)

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

. (31)

Для четвёртого ряда труб:



, (32)

(33)

.

С учётом (33) температурное распределение для потока газа четвёртого ряда труб будет описываться выражением

, (34)

а температурное поле потока воздуха после четвёртого ряда труб с учётом обозначения

, (35)

будет описываться выражением

. (36)

Анализ выражений для стандартизирующих функций температурных распределений потоков газа каждого ряда, температурного поля потока воздуха после каждого ряда позволяет записать следующие выражения в зависимости от номера ряда n, :

; (37)

; (38)

. (39)

При этом и связаны следующими зависимостями:

; (40)

. (41)

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

(42)

;

(43)

;

(44)

.

Следуя аналогии, можно записать выражения для температуры газа при любом числе рядов теплообменника. Полученные выражения позволят рассчитать температуру потоков газа на выходах каждого ряда при , и, в дальнейшем, - температуру на выходе АВО как среднюю температуру потоков

. (45)

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



Библиографический список

1. Коршак А.А., Шаммазов А.М. Основы нефтегазового дела: Учебник для вузов. Уфа.: ООО «ДизайнПолиграфСервис», 2001. 544 с.

2. Девятов Б.Н., Демиденко Н.Д. Теория и методы анализа управляемых распределённых процессов. - Новосибирск: Наука, 1983. 272 с.

3. Бутковский А.Г. Структурная теория распределенных систем. М.: Наука, 1977. 320 с.

4. Рапопорт Э.Я. Анализ и синтез систем автоматического управления с распределенными параметрами: Учеб. пособие. М.: Высш. школа, 2005. 292 с.

Размещено на Allbest.ru

...