МИНИСТЕРСТВО ВЫСШЕГО И СРЕДНЕГО СПЕЦИАЛЬНОГО ОБРАЗОВАНИЯ
УССР
ДОНЕЦКИЙ ОРДЕНА ТРУДОВОГО КРАСНОГО ЗНАМЕНИ ПОЛИТЕХНИЧЕСКИЙ
ИНСТИТУТ
На правах рукописи
ХРАПКО Сергей Александрович
УДК 669.187:541.123 Термодинамическая модель системы металл-шлак для АСУ и машинных экспериментов по оптимизации технологии сталеплавильного процесса Специальность 05.16.02 - "Металлургия чёрных металлов"
Автореферат диссертации на соискание учёной степени кандидата технических наук
Донецк – 1990 Диссертационная работа выполнена в Донецком ордена Трудового Красного Знамени политехническом институте.
Научный руководитель - доктор технических наук, профессор Пономаренко А.Г.
Официальные оппоненты: доктор технических наук, профессор Михайлов Г.Г., кандидат технических наук, старший научный сотрудник Тогобицкая Д.Н.
Ведущее предприятие - Московский институт стали и сплавов.
Защита диссертации состоится " " 1990 г. в час. в аудитории 353 пятого учебного корпуса на заседании специализированного совета Д 068.20.01 в Донецком политехническом институте (340000, г.Донецк, ул.Артема, 58).
С диссертацией можно ознакомиться в библиотеке Донецкого политехнического института.
Автореферат разослан "_"_ 1990 г.
Ученый секретарь специализированного совета Зайцев В.И.
1 из
Общая характеристика работы
Актуальность работы.
Одним из основных резервов совершенствования качества продукции, повышения экономической эффективности и улучшения экологических показателей производства является широкое использование ВТ для управления, исследований, прогнозирования и т.д.
Изучение на ЭВМ сложных физико-химических моделей многокомпонентных систем представляет наибольший интерес в научных и практических приложениях химической термодинамики как новое, перспективное средство исследований. Многие дорогостоящие экспериментальные и опытно-технологических работы сейчас могут быть успешно, гораздо быстрее, и с минимальными затратами выполнены физико-химическим моделированием на компьютерах. Существует также множество задач, единственным способом решения которых является физико-химическое моделирование на ЭВМ.
Всякое управление связано с необходимостью предсказания, предвидения результатов тех или иных возможных воздействий на технологический процесс. В случае таких сложных технологий, какой является металлургическая плавка, задача предвидения является узловой, сосредоточивающей в себе наиболее важные проблемы, лимитирующие развитие таких систем управления.
Настоящая работа посвящена развитию таких моделей для имитации процессов выплавки стали и использованию их с целью решения реальных технологических задач.
Цель работы. На основе экспериментального и теоретического анализа распределения элементов между фазами разработать формализованную модель системы металл-шлак-газ с участием элементов постоянной и переменной валентности для автоматизированных систем управления и проведения машинных экспериментов по оптимизации технологии сталеплавильного процесса.
Идея работы заключается в использовании строгого аппарата классической термодинамики и последних разработок в области уравнений состояния для построения формализованной модели системы металл-шлак и численного моделирования на ее основе реальных металлургических процессов.
Научная новизна.
Работа включает следующие новые научные результаты и положения, выносимые на защиту:
1. Создана и доведена до практической реализации строго формализованная термодинамическая модель распределения элементов в многокомпонентной системе металл-шлак-газ для машинного моделирования реальных процессов и оптимизации технологии на основе разработок, выполненных ранее на кафедре электрометаллургии стали и ферросплавов ДПИ.
2. Разработан метод определения взаимно согласованных констант распределения элементов между металлом и шлаком по неполным массивам опытных данных с целью повышения точности расчетов по указанной модели.
3. Уточнены значения констант распределения наиболее важных примесей в стали.
4. Сформулированы и обоснованы требования и критерии выбора выражений для расчета термодинамических функций высокотемпературных металлургических фаз (уравнений состояния).
5. Предложен метод проверки корректности решений, касающийся состава фаз в сложных многокомпонентных системах на основе критерия устойчивости Гиббса.
6. Установлено, что распространение метода параметров взаимодействия Вагнера на концентрированные растворы, а также его использование для расчета интегральных и других характеристических функций фаз термодинамически необосновано, поскольку 2 из приводит к нарушению основных требований к уравнению состояния, в частности, соотношения Гиббса-Дюгема.
7. Предложен метод "утилизации" имеющейся экспериментальной информации, накопленной в виде параметров Вагнера в рамках модели суб- и квазирегулярных растворов.
Обоснованность и достоверность научных результатов и выводов, приведенных в диссертации, подтверждается их последовательным математическим обоснованием и сопоставлением с экспериментальными и технологическими данными. Адекватность уравнений, полученных в работе, подтверждается методами математической статистики. Машинные эксперименты, проведенные в работе, имеют удовлетворительную воспроизводимость, их результаты согласуются с имеющимися литературными и производственными данными.
Научная значимость. Результаты, представленные в работе, расширяют возможности применения математического аппарата термодинамики и являются шагом в совершенствовании расчетов прогнозирующего характера. Разработанная в виде алгоритмов и программ термодинамическая модель сталеплавильного процесса позволяет оценивать химические составы и количества металла, шлака и газообразных продуктов по ходу плавки с точностью, достаточной для технологических целей, что открывает новые возможности в проектировании АСУ ТП.
Практическое значение работы подтверждается успешным применением созданной модели системы металл-шлак-газ в МИСиС, ИЧМ, НИИМ, ДМетИ, ХМИ АН Каз.ССР, ЧМК при проведении научно-исследовательских работ и в учебном процессе. Применение указанной разработки для машинных экспериментов и термодинамического анализа сталеплавильных и ферросплавных процессов позволяет сократить количество сложных заводских экспериментов при оптимизации технологических параметров процессов.
Разработанный алгоритм принят к внедрению в качестве основного варианта АСУ верхнего уровня агрегата "печь-ковш" в проекте реконструкции ЭСПЦ-2 ЧМК.
Апробация работы. Материалы диссертации доложены и обсуждены на конференциях и семинарах: на научно-технической конференции "Научно-технический прогресс в производстве стали и ферросплавов" (Челябинск, 1988); на Всесоюзном совещании "Базы физико-химических и технологических данных для оптимизации металлургических технологий" (Днепропетровск, 1988); на научном семинаре "Аналитическое представление концентрационной зависимости термодинамических функций в расплавленных металлических и солевых системах. Математические методы и физические модели" (Ленинград,1988); на Всесоюзной конференции "Создание и совершенствование энергосберегающих технологий в пирометаллургии" (Караганда, 1988); на VI Всесоюзной школе-семинаре "Применение математических методов для описания и изучения физико-химических равновесий" (Новосибирск, 1989); на I советско-чехословацком симпозиуме по теории металлургических процессов (Москва, 1989); на II Всесоюзном совещании "Применение ЭВМ в научных исследованиях и разработках" (Днепропетровск, 1989).
Публикации. Основные результаты диссертации изложены в десяти печатных работах.
Объем работы. Диссертация изложена на 115 страницах машинописного текста и содержит введение, 5 глав, заключение, 8 рисунков, 10 таблиц, список использованной литературы из 105 наименований и приложения.
В первой главе проанализировано состояние исследуемого вопроса. Конечная цель плавки стали в общем случае сводится к получению заданной массы жидкой стали с требуемыми химическим составом и температурой. Получение заданного химического состава готовой стали относится к числу сложных задач управления, поскольку она свяиз зана с протеканием сложных физико-химических процессов, часть из которых трудно управляема или вовсе неуправляема.
Во всех своих приложениях математика оперирует моделями рассматриваемых объектов. Модель – это уравнение или система уравнений, связывающих основные, с точки зрения решаемой задачи, свойства реального объекта. Модель плавки – это иерархическая система уравнений трех уровней: модель фазы – модель системы металлшлак-газ – модель плавки.
Можно выделить два крайних способа построения таких моделей: статистические ("чёрный ящик") и теоретические (детерминированные), не требующие подгоночных коэффициентов, определяемых из массива опытных данных.
Наиболее общий путь для сложных объектов – "черный ящик" (эмпирические модели) – используются в качестве первой ступени исследования, когда полностью отсутствует или игнорируется информация о внутренних взаимосвязях параметров объекта.
Применение статистических моделей оправдано при оптимизации характеристик конкретных технологических процессов в условиях, близких к существующим, а моделирование новых условий (экстраполяция) в большинстве случаев невозможно.
Обычно используемые модели имеют некоторую теоретическую основу со статистическим определением параметров, которые не могут быть вычислены теоретически, т.е. часть связей определяется теоретическими соотношениями. Теоретические модели короче и информативнее, т.е. производят более эффективную свёртку информации, применимы к широкому кругу технологий, рассматриваемых как частные случаи. Чем больше доля теоретических уравнений, тем более универсальна модель и выше ее прогностические возможности. Теоретические модели устойчивы по отношению к случайной информации, позволяют делать априорные оценки при полном отсутствии экспериментальных данных.
При рассмотрении теоретических моделей часто противопоставляют, например, термодинамические и кинетические модели на том основании, что все реальные процессы неравновесны и не могут быть описаны в терминах термодинамики, основанной на рассмотрении равновесных состояний. В действительности кинетическая модель не может быть полноценной, если она не опирается на термодинамику. Термодинамика определяет необходимые для кинетических расчетов теоретические пределы и направление протекания процесса, движущие силы, а в ряде случаев и сопротивления (например, условия зарождения новых фаз в виде пузырьков, дисперсных включений), поэтому существенная свёртка модели, основанной на кинетических закономерностях, возможна лишь при использовании термодинамики. Термодинамика позволяет раскрыть смысл многих коэффициентов таких моделей, выразив их через фундаментальные теоретические зависимости и константы. Взаимодействие потоков, их предельное смещение наиболее обоснованно определяет термодинамика необратимых процессов. Чисто кинетическая модель отличается от "черного ящика" лишь видом уравнений связи, коэффициенты которых, как и в "черном ящике" находятся из массива опытных данных, относящихся к конкретному процессу и даже агрегату. Поэтому она не может быть использована для априорных количественных оценок. Принципиальный момент состоит в том, что термодинамика позволяет построить работоспособную общую модель и делать априорные количественные оценки без адаптации к конкретному процессу (т.е. без подбора констант по опытным данным).
Известные в настоящее время модели наряду с теоретически обоснованными зависимостями включают значительное количество эмпирических связей и коэффициентов. В обзоре главное внимание уделено термодинамическим детерминированным моделям многофазных систем и выбору термодинамических функций неидеальных растворов сложного состава.
Детерминированные термодинамические модели системы металл-шлак включают в себя законы сохранения (уравнения материального баланса, электронейтральности, и др.) и условия термодинамического равновесия, определяемые по минимуму свободной энергии (Гиббса или Гельмгольца) или по максимуму энтропии. В соответствии с этим модель плавки представляет собой сложную систему уравнений, решение которой применительно к условиям реальных процессов встречает определенные трудности. Поэтому известные в настоящее время модели плавки громоздки, неуниверсальны, используют большое число упрощающих предположений и не доведены до реальных технологических расчетов. Кроме того, многие из них не учитывают неидеальность конденсированных фаз, что дополнительно ограничивает их применимость.
Отмечено, что при термодинамическом моделировании металлургических равновесий выбор функций состав-свойство (модели уравнений состояния) конденсированных фаз определяет наиболее важные результаты моделирования сталеплавильных процессов.
Строгая теория жидкостей (и вообще конденсированных фаз) в настоящее время позволяет получить лишь качественное описание термодинамических свойств простейших систем. Широкое распространение в расчетах получили упрощенные методы, при которых заранее принимаются допущения о структуре расплава и характере взаимодействия его частиц и учитывается только конфигурационная часть энтропии. Существует также и большая группа моделей, в которых выражения для термодинамических свойств расплава получаются формальным разложением в полиномиальный ряд и последующим определением коэффициентов методами математической статистики по массивам экспериментальных данных.
Рассмотрены и проанализированы модели, наиболее часто используемые для описания термодинамических свойств металлических и шлаковых расплавов. На основании анализа состояния вопроса сформулированы задачи исследования:
1. Разработать формализованную модель систем металл-шлак и металл-шлак-газ с участием элементов постоянной и переменной валентности, учитывающую неидеальность конденсированных фаз;
2. Выбрать математически и термодинамически обоснованные выражения для расчета термодинамических характеристик компонентов (элементов) в металле и шлаке и сформулировать критерии такого отбора;
3. Разработать алгоритм анализа полученного решения для системы фаз сложного состава на устойчивость (стабильность);
4. Обосновать единственность решения, получаемого при моделировании плавки на ЭВМ и оценить его точность;
5. Разработать метод согласованного определения констант распределения элементов между металлом и шлаком по неполным массивам опытных данных;
6. Разработать способ статистической адаптации, т.е. определения поправок на неравновесность, кинетические и др. эффекты, прямо не учитываемые моделью;
7. Реализовать перечисленные выше пункты программно для больших, малых и микро-ЭВМ, используемых в промышленных АСУ.
Создание моделей фаз в задачи исследования не входило.
Вторая глава посвящена выбору моделей уравнений состояния для металлической и шлаковой фаз из имеющихся в литературе на основе соотношений общей термодинамики. В металлургической литературе этот вопрос известен как проблема активностей.
Термодинамика постулирует сам факт существования уравнения состояния, например, в виде где G – интегральная функция. При описании химических равновесий используются производные функции (1).
Во всех практических расчетах используют те или иные конкретные выражения – модели уравнения состояния, которые лишь приближенно характеризуют фазу. Модель фазы (1) является первой ступенью в иерархической системе моделей: модель фазы – модель системы металл-шлак-газ – модель плавки.
Модели конденсированных фаз, несмотря на многочисленные разработки, являются наиболее уязвимым местом всех термодинамических расчетов реальных систем.
Известные модели металлических и шлаковых растворов не всегда обеспечены численными значениями параметров, не определены границы их применимости, что часто затрудняет их использование даже для грубых априорных оценок.
Модели фаз, независимо от способа их получения, являются аппроксимирующими функциями, что налагает на них ряд специфических требований с точки зрения возможности и корректности их математических преобразований. Кроме того, необходимо, чтобы каждое из используемых уравнений для отдельных фаз удовлетворяло ряду общих требований термодинамики растворов (фундаментальное уравнение Гиббса, предельные законы и др.). И, наконец, поскольку эти уравнения используются для описания системы фаз, на них должны распространяться дополнительные требования, налагаемые этой системой, такие как выбор компонентов, устойчивость во всей области изменения параметров и др.
Основным требованием к модели является ее адекватность реальному объекту, которая достигается выбором вида функций и оптимизацией входящих в них коэффициентов. Математические преобразования аппроксимирующих функций, в частности их дифференцирование, может приводить к катастрофической утрате адекватности (так называемые некорректно поставленные задачи). По этой причине, например, не всегда корректны оценки активностей по параметрам моделей, определяемым из интегральной энтальпии.
Поскольку выбранные функции для описания термодинамических свойств расплавов должны удовлетворять фундаментальному уравнению Гиббса (или уравнению Гиббса-Дюгема в частном случае постоянства температуры и давления), то модель фазы должна представлять собой одну из характеристических функций состояния, т.е. интегральное свойство фазы. Кроме того, уравнение должно быть аналитическим, непрерывным и дифференцируемым (для анализа устойчивости фазы как минимум дважды) во всей области существования фазы.
Если модель фазы представлена уравнениями для парциальных характеристик (например, разложение Вагнера), то должны соблюдаться условия совместности (интегрируемости), которые гарантируют, что дифференцирование двух или более уравнений Фi = f i ( x ) приводит к совпадению производных высшего порядка от функции Ф :
Если условия совместности (2) не выполняются при любых значениях независимых параметров в области определения фазы, то это говорит о внутренней противоречивости модели и приводит, например, к нарушению уравнения Гиббса-Дюгема.
При построении формализованных моделей фаз важно минимизировать количество параметров модели при сохранении удовлетворительной точности аппроксимации свойств. Этого можно достичь выбором уравнений с использованием атомных параметров ai, взамен, например, парных aij, при этом каждому компоненту раствора в простейшем случае соответствует один постоянный коэффициент, не зависящий от состава и типа раствора. Для формализации модели требуется также идентичность расчетных формул для всех компонентов фазы.
В работе сформулированы и обоснованы предельные условия для уравнений состояния в виде дифференциальных формулировок законов Рауля и Генри (в терминах активности):
который позволяет отделить устойчивые состояния фазы (стабильные и метастабильные) от полностью неустойчивых (лабильных).
При нарушении положительной определенности матрицы вторых частных производных энергии Гиббса фазы по числам молей компонентов можно обоснованно утверждать, что при указанном составе фазы либо должно произойти ее расслаивание, либо модель неприменима в этом диапазоне концентраций компонентов. В любом случае результаты расчета не будут соответствовать действительности и необходимо заменить используемую модель на более общую, описывающую больший диапазон составов фазы.
В четвертой главе дается подробное описание алгоритма решения задачи расчета равновесных составов металла, шлака и газа.
В пятой главе приводятся примеры проверки адекватности модели. Для выяснения работоспособности модели и точности получаемого решения были проведены многочисленные машинные эксперименты, имитирующие выплавку стали и ферросплавов в различных агрегатах. В табл.1 приведены результаты расчета состава металла при выплавке стали ШХ15 в ЭСПЦ-2 ЧМК. При этом по заводским картам плавок задавались количества шихтовых материалов (лом, чугун, шлакообразующие, окислители, легирующие и восстанователи) и рассчитывались составы и количества металла и шлака.
Наибольшая ошибка при расчете состава наблюдается в периоде расплавления и окисления, что связано с неизвестным составом завалки, отсутствием расходомеров кислорода при продувке, а также хозяйственной деятельностью предприятия. По остальным пробам согласие хорошее. Максимальная погрешность предсказания состава готового металла по результатам расчетов 30 плавок стали ШХ15 составила по элементам (в процентах): [O] – 0.0005, [Si] – 0.05, [Mn] – 0.06, [P] – 0.005, [S] – 0.005, [Cr] – 0.05.
Адаптация проводилась по восьми плавкам, взятым из заводского отчета.
На рис.1 приведена иллюстрация прогнозных возможностей модели по описанию конвертерных процессов. Здесь сравниваются расчетные и фактические составы металла в кислородном конвертере при верхней и нижней продувке. Модель достаточно адекватно описывает донную продувку, которая ближе к равновесию и лучше усреднена по химическому составу, при верхней продувке явно просматривается аккумуляция ванной кислорода и запаздывание окисления примесей. Расчеты проводились без использования каких-либо данных, относящихся к конкретному агрегату.
Таблица 1. Состав металла промышленных плавок стали ШХ № пробы: 1 – расплавление, 2 – продувка кислородом, 3 – легирование, 4 – корректировка, 5 – ковшевая проба, н.д. – нет данных. Числитель – расчетные данные, знаменатель – фактические.
Рис. 1. Моделирование кислородно-конвертерного процесса.
1. Разработана и опробована термодинамическая модель распределения элементов в системе металл-шлак-газ для машинного моделирования процессов и оптимизации технологии. Термодинамический расчет системы металл-шлак-газ, опирающийся на фундаментальные константы и зависимости позволяет, по крайней мере, в принципе, получить априорные количественные оценки равновесного состава фаз. Поэтому термодинамический расчет равновесия принят в качестве основы модели, а кинетические и другие факторы учитываются дополнительно как необходимые поправки. В модели предусмотрена система адаптации, позволяющая в случае необходимости уточнять расчет, при этом коэффициенты, не определяемые через фундаментальные константы находятся статистически.
2. Разработан метод определения взаимно согласованных констант распределения элементов между металлом и шлаком по неполным массивам опытных данных с целью повышения точности расчетов по указанной модели.
3. Уточнены значения констант распределения наиболее важных примесей в стали.
4. Сформулированы и обоснованы требования и критерии выбора выражений для расчета термодинамических функций высокотемпературных металлургических фаз (уравнений состояния).
5. Установлено, что метод параметров взаимодействия Вагнера не может служить основой для построения общей модели фазы, поскольку приводит к нарушению основных требований к уравнению состояния, в частности, уравнения Гиббса-Дюгема.
6. Предложен метод "утилизации" имеющейся экспериментальной информации в виде параметров Вагнера в рамках модели суб- и квазирегулярных растворов.
Содержание диссертации опубликовано в следующих работах 1. Совпель С.В., Храпко С.А., Пономаренко А.Г., Зац Е.Л. Расчет состава оксидных включений в стали/ Донецк. политехн. ин-т.- Донецк,1986.- 7с.- Деп. в Черметинформации 20.08.86, N 3547.
2. Рудненко Т.Б., Пономаренко А.Г., Зац Е.Л., Гулевская Л.И., Храпко С.А., Боровко А.И., Свиридов О.В., Волков А.С. О влиянии состава флюса на качество промышленных слитков стали 9Х2МФ электрошлакового переплава/ Донецк. политехн. ин-т.Донецк,1986.- 25с.-Деп. в Черметинформации 26.09.86, N 3644.
3. Храпко С.А., Иноземцева Е.Н. Термодинамическое моделирование на ЭВМ системы металл-шлак-газ// Научно-технич. прогресс в производстве стали и ферросплавов:
Тез. докл. отраслевой научно-технич. конференции.- Челябинск: НИИМ, 1988. - С.
4. Пономаренко А.Г., Иноземцева Е.Н., Храпко С.А. Термодинамическая модель для машинных экспериментов при оптимизации сталеплавильных и ферросплавных процессов//Базы физ.-хим. и технолог.данных для оптимизации металлург. технологий: Сб.тез. Всесоюзн.межотрасл.сов.-Днепропетровск: ИЧМ, 1988.- С.147-153.
5. Храпко С.А., Пономаренко А.Г. Описание термодинамических свойств металлических расплавов//Тез.докл.научно-технич.конф. Создание и совершенств. энергосбер.
технологий в пирометаллургии.- Караганда:ХМИ АН Каз.ССР, 1988.- Часть II.- С. 81Храпко С.А., Пономаренко А.Г. Описание термодинамических свойств металлических растворов// Применение математ.методов для описания и изуч.физ.хим.равновесий. Тез.докл.VI Всесоюзн. школы-семинара.- Новосибирск, 1989.- Часть I.- С.73-74.
7. Пономаренко А.Г., Иноземцева Е.Н., Храпко С.А. Расчет равновесного состава конденсированных фаз металлургических систем // Применение математических методов для описания и изучения физико-химических равновесий. Тезисы докладов VI Всесоюзной школы-семинара.- Новосибирск, 1989.- Часть I.- С.91-92.
8. Храпко С.А., Пономаренко А.Г. Выбор уравнений для построения термодинамических моделей металлических растворов//Тез.докл. I советско - чехословацкого симпозиума по теории металлургических процессов. Ч.I. Структура и физико-химические свойства металлических и шлаковых расплавов.-Москва,1989.-С.118-122.
9. Храпко С.А., Иноземцева Е.Н. Развитие термодинамической модели системы металл-шлак-газ// Тезисы докладов II Всесоюзного совещания. Применение ЭВМ в научных исследованиях и разработках.- Москва, 1989.- С.22-23.
10. Храпко С.А., Иноземцева Е.Н., Пономаренко А.Г. Машинное моделирование для отработки технологии выплавки стали/ Повышение качества и эффективности производства электростали. Тематич.сб.науч.тр.-Челябинск: Металлургия, 1989.-С.71-74.