ФЭА / АИТ / Курсовая работа по моделированию "Математическая модель трубчатой печи"
(автор - student, добавлено - 14-06-2014, 21:18)
СКАЧАТЬ:
Содержание
Введение
1.1. Описание технологической схемы блока вакуумной перегонки. 1.2. Основные сведения о конструкциях и предъявляемых требованиях к трубчатым печам [1] 2. Экспериментальная часть 2.1. Составление статической модели трубчатой печи2.2Множественная корреляция [2] 3. Рас четная часть 3.1. Составление теплового баланса трубчатой печи [4,5] 3.2Составление динамики процесса, протекающего в трубчатой печи 3.3 Оптимизация процесса методом перебора [8] 4. Проектная часть 4.1 Проектирование печи с помощью универсальной моделирующей программы ChemCad 5. Список литературы 6. Приложение
Введение Особенностью современных процессов химической технологии, протекающих с высокими скоростями при высоких температурах и давлениях в многофазных системах, является их сложность. Эта сложность проявляется в значительном числе и многообразии параметров, определяющих течение процессов, в большом числе внутренних связей между параметрами, в их взаимном влиянии. Поэтому удается лишь составить систему дифференциальных уравнений для их описания и установить условия однозначности. Решить же эти уравнения известными в математике методами обычно не представляется возможным. В подобных случаях используют метод моделирования. Математическое моделирование – это по существу определение свойств и характеристик рассматриваемого явления путем решения системы уравнений, описывающих этот процесс, - математической модели. Моделирование можно осуществлять двумя основными методами: методом обобщенных переменных, или методом подобия (физическое моделирование), и методом численного эксперимента (математическое моделирование). Последнее позволяет резко сократить сроки научных и проектных разработок. Обычная схема процесса математического моделирования включает 8 основных этапов: 1. Постановка задачи. 2. Анализ теоретических основ процесса (составление физической модели объекта). 3. Составление математической модели объекта. 4. Алгоритмизация математической модели. 5. Параметрическая идентификация модели. 6. Проверка адекватности математической модели. 7. Моделирование процесса. 8. Анализ полученной информации. На основании проведенного анализа принимают решение – выдать рекомендации для практической реализации или продолжать исследование.
2.1. Описание технологической схемы блока вакуумной перегонки.
Основное назначение установки вакуумной перегонки мазута топливного профиля – получение вакуумного газойля широкого фракционного состава О четкости разделения мазута обычно судят по фракционному составу и цвету вакуумного газойля. Последний показатель косвенно характеризует содержание смолисто-асфальтеновых веществ, т.е. коксуемость и содержание металлов. Металлы, особенно никель и ванадий, оказывают отрицательное влияние на активность, селективность и срок службы катализаторов процессов гидрооблагораживания и каталитической переработки газойлей. Поэтому при эксплуатации промышленных установок ВТ исключительно важно уменьшить унос жидкости в концентрационную секцию вакуумной колонны в виде брызг, пены, тумана и.д. В этой связи вакуумной колонны по топливному варианту имеют при небольшом числе тарелок развитию питательную секцию: отбойники из сеток и промывные тарелки, где организуется рециркуляция затемненного продукта. Для предотвращения попадания металлоорганических соединений в вакуумный газойль иногда вводят в сырье в небольших количествах антипенную присадку типа силоксан. В процессах вакуумной перегонки, помимо проблемы уноса жидкости, усиленное внимание уделяется обеспечению благоприятных условий для максимального отбора целевого продукта без заметного его разложения .Многолетним опытом эксплуатации промышленных установок ВТ установлено, что нагрев мазута в печи выше 420-425оС вызывает интенсивноеобразование газов разложения, закоксование и прогар труб печи, осмоление вакуумного газойля. При этом чем тяжелее нефть, тем более интенсивно идет газообразование и термодеструкция высокомолекулярных соединений сырья. Вседствии этого при нагреве мазута до максимально допустимой температуры уменьшают время его пребывания в печи, устраивая многопоточные змеевики, применяют печи двустороннего облучения, в змеевик печи подают водяной пар и уменьшают длину трансферного трубопровода. Для снижения температуры низа колонны организуют рецикл частично охлажденного гудрона. С целью снижения давления на участке испарения печи концевые змеевики выполняют из труб большего диаметра и уменьшают перепад высоты между вводом мазута в колонну и выходом его из печи. В вакуумной колонне применяют ограниченное число тарелок с низким гидравлическим сопротивлением или насадку; используют вакуумсоздающий системы, обеспечивающие достаточно глубокий вакуум. Число тарелок в отгонной секции также должно быть ограничено, чтобы обеспечить малое время пребывания нагретого гудрона. С этой целью одновременно уменьшают диаметр куба колонны. В процессах вакуумной перегонки мазута по топливному варианту преимущественно используют схему однократного испарения,применяя одну сложную ректификационную колонну с выводом дистилятных фракций через отпарные колонны или без них. При использовании отпарных колонн по высоте основной вакуумной колонны организуют несколько циркуляционных орошений. Мазут, отбираемый с низа атмосферной колонны блок АТ, прокачивают параллельным потоком через печь в вакуумную колонну. Смесь нефтяных и водяных паров, газы разложения с верха вакуумной колонны подают в вакуумсоздающую систему. После конденсации и охлаждения в кондесаторе-холодильнике ее разделяют в газосепараторе на газовую и жидкую фазы. Газы отсасывают трехступенчатым пароэжекторным вакуумным насосом, а конденсаты направляют в отстойник для отделения нефтепродукта от водного конденсата. Верхним боковым погоном вакуумной колонны отбирают фракцию легкого вакуумного газойля. Часть его после охлаждения в теплообменниках возвращают на верх колонны в качестве верхнего циркуляционного орошения. Вторым боковым погоном отбирают широкую газойливую фракцию. Часть ее после охлаждения используют как среднее циркуляционное орошение вакуумной колонны. Балансовое количество целевого продукта вакуумного газойля после теплообменников и холодильников выводят с установки и направляют на дальнейшую переработку. С нижней тарелки концентрационной части колонны выводят затемненную фракцию, часть которой используют как нижнее циркуляционное орошение, часть- выводят с установки или используют как рецикл вместе с загрузкой вакуумной печи. С низа вакуумной колонны отбирают гудрон и после охлаждения направляют на дальнейшую переработку. Часть гудрона после охлаждения в теплообменнике возвращают в низ колонны в качестве квенчинга. В низ вакуумной колонны и змеевик печи подают водяной пар.
1.2. Основные сведения о конструкциях и предъявляемых требованиях к трубчатым печам [1] В настоящее время на большинстве объектов промысловой подготовки нефти в качестве генераторов тепла (нагревателей) используются различные типы трубчатых печей. Это один из наиболее сложных видов оборудования многих технологических установок подготовки нефти (на современных нефтеперерабатывающих заводах и химических производствах они являются основными генераторами тепла). Современная трубчатая печь (см. рис.1) состоит обычно из двух камер: камеры сгорания (радиационная камера) и камеры конвекции. В камере радиации сжигается топливо, т.е. в принципе она представляет собой топочное устройство. В этой камере размещаются рабочие трубы, по которым движется нагреваемая жидкость. Поскольку теплота, выделяющаяся при сгорании топлива, воспринимается этими трубами главным образом путем лучистого теплообмена (радиацией), эти трубы принято называть радиант-ными. В камеру конвекции поступают продукты сгорания топлива из радиационной камеры. Здесь они обтекают ряды труб, по которым движется жидкость. Восприятие теплоты продуктов сгорания топлива здесь осуществляется, главным образом, путем конвекции. По этой причине рабочие трубы, расположенные в этой части, называют конвекционными. Каркасы трубчатых печей представляют собой систему вертикальных колонн, связанных между собой горизонтальными и наклонными балками, которые образуют жесткую пространственную конструкцию. Назначение каркаса - воспринимать нагрузку от трубного змеевика, трубных решеток, подвесок, кровли, подвесного свода, стен, площадок для обслуживания, горелочных устройств и др. Так как элементы каркаса изготавливают из нежаростойких сталей, они защищаются от воздействия высоких температур обмуровкой и тепловой изоляцией. Обмуровка трубчатых печей - это часть конструкций, служащая для защиты от высоких температур и обеспечения необходимой плотности топки, радиационных и конвекционных камер, а также стен газового тракта (боровов и дымовых труб). Обмуровка может состоять из нескольких слоев, различающихся теплостойкостью. Для придания прочности, плотности и защиты от атмосферных воздействий обмуровку снаружи закрывают стальным кожухом. Обмуровку крепят на специальных подвесках и кронштейнах, которые соединены с каркасом печи. Боковые поверхности обмуровочного кирпича часто изготовляют волнистыми, для обеспечения большей герметичности внутрипечного пространства. Для возможности теплового расширения обмуровки устраиваются температурные швы , которые заполняют легко деформируемым теплостойким, мягким материалом (например, асбестом или каолиновой ватой) Основные требования к обмуровке:
Устройства для сжигания топлива - это форсунки и горелки. Форсунки служат для сжигания жидкого топлива, горелки - для газообразного. Жидкое топливо распыливают водяным паром, сжатым воздухом или механически (энергией, вытекающей с большой скоростью струи самого топлива). Газовые горелки служат для перемешивания топливного газа с воздухом и зажигания образовавшейся газо-воздушной смеси. Существует сравнительно большое количество сведений об их расчете, проектировании, монтаже, наладке и эксплуатации. Поэтому целесообразно остановиться только на основных требованиях, предъявляемых к газогорелочным устройствам, применяемым на трубчатых печах. Основные требования, предъявляемые к газогорелочным устройствам:
Маркировка типовых трубчатых печеи Условное обозначение современных печей состоит из двух букв, индекса и цифры. Первая буква выражает конструктивное исполнение: А - узкокамерные трубчатые печи с верхним отводом дымовых газов и центральным экраном с горизонтально расположенными трубами змеевика; Г - узкокамерные трубчатые печи с верхним отводом дымовых газов и настенными экранами с горизонтально расположенными трубами змеевика; В - узкокамерные секционные трубчатые печи с вертикально расположенными трубами змеевика; С - секционные трубчатые печи с прямоугольно и горизонтально витым трубным змеевиком ; Ц - цилиндрические трубчатые печи с верхней камерой конвекции, горизонтальными трубами конвекционной части змеевика и вертикальными - радиантной части; К - цилиндрические трубчатые печи с боковой кольцевой камерой конвекции и вертикальным расположением конвективных труб. Цифровой индекс к первой букве показывает число рядов экрана. Если экран однорядный, - индекс не указывается. Вторая буква обозначения выражает способ сжигания топлива С - свободный вертикальный факел; В - свободный вертикальный факел с позонным подводом воздуха по высоте факела; Н - настильный факел; Д - настильный факел с дифференцированным подводом воздуха по высоте факела; Б - беспламенное горение с излучающими стенами из панельных горелок. Знак ' (штрих) ко второй букве означает смещение горелки от центра в сторону входа продукта. Цифра после второго буквенного обозначения показывает число радиантных камер или секций печи. Цифры в числителе дроби показывают поверхность нагрева труб, м . Цифры в знаменателе показывают длину или высоту топки, в метрах (м). Например: А2 Б 2 - узкокамерная трубчатая печь с центральным горизонтальным двухрядным радиационным экраном, верхним отводом дымовых газов, с излучающими стенками топки,оборудованными панельными горелками беспламенного горения, двумя радиантными камерами, с поверхностью нагрева радиантных труб, равной 115 м2, длиной топки, равной 6 м. Основные показатели работы технического уровня трубчатых печей. Для характеристики печи наиболее важны следующие (основные) показатели: Производительность печи - это количество сырья (теплоносителя), подаваемого в печь для нагрева в единицу времени. Производительность современных печей по нагреваемой среде до 8105 кг/ч (222,2 кг/с). Тепловая мощность печи - это количество тепла, которое может быть воспринято нагреваемым теплоносителем в единицу времени. В промышленности используются печи тепловой мощностью от 0.5 до 100 МВт. Теплонапряженность поверхности нагрева - это количество тепла, переданного нагреваемому теплоносителю, через единицу поверхности (1 м) змеевиковых труб. Единицей измерения является кВт. Поскольку радиантные и конвекционные трубы работают в разных условиях, различают теплонапряженность радиантных труб, теплонапряженность конвекционных труб и среднюю теплонапряженность труб печи. Величина теплонапряженности поверхности нагрева является одной из важнейших характеристик, она выражает эффективность передачи тепла через поверхность нагрева. Следовательно, чем выше теплонапряженность, тем меньших размеров требуется поверхность трубного змеевика для передачи требуемого количества тепла. С уменьшением требуемой поверхности нагрева уменьшаются и требуемые габариты печи. Однако очень большая теплонапряженность поверхности нагрева может вызвать разложение нагреваемого продукта. Например, у нефти при больших теплонапряженностях вследствие повышения температуры происходит коксование. Продукты коксования могут наслаиваться на внутренних поверхностях труб и образовывать отложения. Эти отложения обладают низкой теплопроводностью, что служит одной из причин прогара труб из-за их перегрева. Температура стенок труб зависит от температуры нагреваемого продукта и скорости его движения. Чем ниже температура продукта и выше скорость его движения, тем больше допускаемая теплонапряженность труб. Существенное влияние на допускаемую теплонапряженность оказывает также вид нагреваемого продукта. Чем более тяжелое сырье подвергается нагреву, тем меньше допускаемая теплонапряженность труб. В трубах радиационной части змеевика допускаемая тепло-напряженность выше, чем в трубах конвекционной части Тепловая напряженность топочного пространства - это количество теплоты, выделенное при сгорании топлива в единицу объема в единицу времени. Измеряется эта величина в . В современных печах эта величина может составлять 50... 100 (что в несколько раз меньше, чем для топок паровых котлов) и определяется возможностью размещения необходимой радиантной поверхности в топочной камере (а не процессом горения топлива). Коэффициент полезного действия трубчатой печи (к.п.д.) представляет собой отношение полезно использованной теплоты к общему количеству теплоты, выделившейся при сгорании топлива. К.п.д. печи зависит главным образом от коэффициента избытка воздуха в топочном пространстве, температуры уходящих продуктов сгорания топлива, потерь теплоты в окружающую среду и качества сжигания топлива. Современные трубчатые печи имеют к.п.д. в пределах 0.65...0.85(65...85%). Основными параметрами, характеризующими технический уровень печи, как в стадии проектирования или выбора, так и в период эксплуатации являются: степень равномерности нагрева труб по их окружности степень равномерности нагрева по их длине ,степень эффективности распределения тепла по трубному экрану , степень использования поверхности нагрева ,степень эффективности компоновки трубного экрана коэффициент металлоемкости радиантной части печи . Абсолютные значения перечисленных показателей во многом определяют величины допускаемых тепловых напряжений,ч то существенным образом влияет на сроки службы и безаварийной эксплуатации змеевиков. Максимально достигаемые величины степени равномерности нагрева любых труб продуктового змеевика по окружности и по длине всегда меньше единицы и зависят от ряда факторов, основными из которых являются: способ сжигания, направления движения фронта горения и продуктов сгорания топлива, распределение теплоты, а также тип и расположение трубного экрана. Степень эффективности распределения по трубному экрану представляет собой отношение фактического тепла, переданного трубному экрану, к допускаемому количеству тепла при условии, что ни в одной трубе экрана показатели фактического теплонапряжения не превышают теоретически допускаемых и равны (по меньшей мере) , (1) где qф и qm -среднефактическое и среднетеоретическое теплонапряжения трубного экрана. Степень использования поверхности нагрева можно определить по формуле: . (2) Значение соответствует случаю, когда температура стенки всего экрана (всех труб, входящих в соответствующий экран), по всей поверхности, равна допускаемой температуре. Степень эффективности компоновки трубного экрана - это характеристика полноты использования площади экранированной поверхности стен. Она определяется по формуле =, (3) где - отношение площади проекции радиантных труб к площади экранированной стены. Коэффициент металлоемкости радиантной части характеризует эффективность использования металла, затраченного на изготовление труб . (4) Из полученных формул очевидно, что более эффективными являются печи с большими значениями и меньшими .
2. Экспериментальная часть 2.1. Составление статической модели трубчатой печиДля построения статической модели используем фактические данные о работе печи. Эти данные можно получить путем наблюдения – целенаправленного восприятия установки без активного вмешательства в ее эксплуатацию или при выполнении эксперимента, когда исследователь по заранее продуманному плану устанавливает величину управляющих факторов, чтобы добиться определенного результата. Очевидно, что повседневная эксплуатация объекта является примером пассивного накопления данных о его работе путем регистрации установленных в технической документации параметров, например, в «Оперативном (суточном) журнале». Статическую модель можно легко построить по результатам наблюдения за эксплуатацией объекта, зафиксированным в эксплуатационной документации, с помощью математического пакета Mathcad 2000 Professional. Первоначально создается таблица данных о работе объекта, соответствующая форме журнала, в котором фиксируются параметры. В данной работе используется 60 значений параметров. Таблица1
Необходимо найти такую формулу, которая “наилучшим образом” описывала бы, указанную в таблице зависимость. Первоначально графически покажем зависимость выходного параметра от каждого из входного параметра. Для этого построим поля корреляции, определим внешний вид регрессии и запишем уравнение регрессии от каждого параметра [2, 3].
Рис. 1. Зависимость температуры на выходе температуры на входе Для полученной линии регресии можно записать следующее уравнение регрессии:
где х=t, y=T
Рис. 2. Зависимость температуры на выходе от давления на входе В качестве приближающей функции для данного вида линии регрессии воспользуемся следующей функцией:
где х=Р, y=T. Определим влияние каждого входного параметра на основной выходной параметр путем подсчета коэффициентов корреляции. Для определения коэффициентов в уравнениях регрессии воспользуемся методом наименьших квадратов (МНК). Задача определения коэффициентов в уравнении регрессии сводится к определению минимума функции многих переменных, т.е. должно выполняться следующее условие:
Зависимость выходной температуры от входной выражается логарифмической функцией следующего вида:
Для перехода к линейной функции достаточно сделать подстановку ln(x)=u. Отсюда следует, что для нахождения значений с0 и с1 необходимо прологарифмировать значения аргумента в исходной таблице 1. Поскольку в результате подстановки получили функцию линейного вида, то и процесс нахождения коэффициентов регрессии будет аналогичен линейной функции. Система нормальных уравнений для этого случая имеет вид:
или
Коэффициенты с0 и с1 легко найти в этом случае с помощью определителей
Уравнение регрессии, выражающее зависимость выходной температуры от давления на входе имеет линейную форму , следовательно принцип нахождения коэффициентов b0 и b1 аналогичен предыдущему, т.е. Система нормальных уравнений для этого случая имеет вид:
или
Коэффициенты b0 и b1 так же находятся с помощью определителей
Для оценки линейной связи (5), (6) вычисляется выборочный коэффициент корреляции r*, который является показателем значимости коэффициентов регрессии: 0,75 |r*|1.
где хi, yi – значения соответствующих факторов, – средние значения факторов, sx sy – выборочные среднеквадратичные отклонения, N – объем выборки. Для оценки линейной связи с учетом (12) и (16) формула (17) принимает следующий вид:
Тогда
Значения всех полученных коэффициентов регрессии близки к 1, это свидетельствует о сильной связи между выходным параметром yi факторами xi.. Множественная корреляция [2] Из найденных значений парной корреляции следует, что параметры t и P влияют на величину выходного параметра и математическая модель определяется уравнением множественной регрессии для случая двух факторов.
В этом случае мы уже имеем дело не с линией регрессии, а с поверхностью регрессии (поскольку выходной параметр зависит от двух входных параметров, т.е. k=2). При построении поверхности отклика на координатных осях факторного пространства откладываются численные значения параметров (факторов). Исходный статистический материал представлен в таблице 1. Прежде всего, необходимо перейти от натурального масштаба к новому. Для этого проведем нормировку всех значений случайных величин по формулам:
где yi0, x1i0, x2i0 – нормированные значения соответствующих факторов, - средние значения факторов, sy, sx1, sx2, – среднеквадратичные отклонения.
Применительно к нашему случаю:
Тогда
В таблице 2 представлен исходный статистический материал в новом масштабе. Таблица 2
Выборочный коэффициент корреляции при этом равен:
Вычисленный по формуле (21) выборочный коэффициент корреляции равен коэффициенту корреляции между переменными, выраженными в натуральном масштабе r*ух (17). Уравнение регрессии между нормированными переменными не имеет свободного члена и принимает вид:
Коэффициенты уравнения (22) находятся из условия:
Условия минимума функции S определяются так же, как в случае зависимости от одной переменной:
и система нормальных уравнений имеет вид:
Умножим левую и правую части уравнений на . В результате при каждом коэффициенте получается, согласно (25), выборочный коэффициент корреляции r*. Принимая во внимание, получаем систему нормальных уравнений:
С учетом вычисленных значений коэффициентов уравнение регрессии между нормированными переменными примет следующий вид:
Рассчитаем коэффициент множественной корреляции R. Коэффициент множественной корреляции служит показателем силы связи в случае множественной регрессии: .
Т.к. у нас выборка небольшого объема, то в величину R необходимо внести коррекцию на систематическую ошибку. Чем меньше число степеней свободы выборки , тем сильнее преувеличивается сила связи, оцениваемая коэффициентом корреляции. Формула для коррекции:
где R/ – скорректированное значение коэффициента множественной корреляции; l – число коэффициентов уравнения регрессии.
Для практического использования уравнения (22) необходимо перейти к натуральному масштабу по формулам:
Таким образом, наша модель для трех факторов имеет следующий вид:
Проверка значимости коэффициентов в уравнении регрессии. Проверка адекватности и работоспособности регрессионной модели Оценим адекватность и работоспособность полученной экспериментальной факторной модели. Определим среднее значение функции отклика у и дисперсию модели среднего Sy, характеризующую рассеяние результатов эксперимента относительно и оценивающую погрешность модели среднего:
где Yi - значение функции отклика в i-й точке; N - количество проведенных опытов. Осуществим проверку пригодности полученного уравнения регрессии. Для этого сначала вычислим по уравнению регрессии предсказываемые значения функции отклика . В уравнение регрессии при вычислениях подставляют значения нормированных факторов . Затем определяется остаточная дисперсия , оценивающая погрешность полученной модели:
где k - число коэффициентов в уравнении регрессии. В нашем случае k=3.
После определения оценок коэффициентов регрессии необходимо проверить гипотезу о значимости коэффициентов bi. Лучше всего это сделать в виде нуль-гипотезы, т.е. гипотезы о равенстве . Принимается нулевая гипотеза о том, что bi=0, и осуществляется проверка, отличаются ли статистически значимые оценки коэффициентов bi от нуля. Значимость bi проверяют по критерию Стьюдента, используя формулу (35).
где S2{bi} - дисперсия ошибки определения коэффициента bi. При вычислении дисперсии , оценивающей погрешности определения коэффициентов bi используется остаточная дисперсия .
Табличное значение tкр=1,98 (nз=N(m-1)=120, q=5%) (приложение 2). Т.к. все вычисленные значения параметра ti превышают табличное значение tкр, то все найденные оценки коэффициентов признаются значимыми и войдут в модель. По критерию Фишера F и коэффициенту детерминации R2 оценим качество предсказания, обеспечиваемого полученной регрессионной моделью. При определении критерия Фишера принимается нулевая гипотеза. Нулевая гипотеза гласит о том, что модель среднего достаточно хорошо описывает исследуемый процесс. Регрессионная модель окажется адекватной, если выдвинутая гипотеза будет опровергнута. По критерию Фишера сравнивают дисперсии оцениваемой и противопоставляемой моделей. Последняя должна быть более точной, чем оцениваемая модель. Поэтому в данном случае критерий Фишера равен отношению дисперсии модели среднего к остаточной дисперсии : Copyright 2021. Для правильного отображения сайта рекомендуем обновить Ваш браузер до последней версии!
|