На правах рукописи
Мокрова Наталия Владиславовна
ОСНОВЫ ПОСТРОЕНИЯ АВТОМАТИЗИРОВАННЫХ СИСТЕМ ИЕРАРХИЧЕСКИ-ВЗАИМОСВЯЗАННОГО УПРАВЛЕНИЯ СЛОЖНЫМИ ТЕХНОЛОГИЧЕСКИМИ ПРОЦЕССАМИ
ПЕРЕРАБОТКИ ПРИРОДНОГО СЫРЬЯ
Специальность:
05.13.06 - автоматизация и управление технологическими
процессами и производствами (промышленность)
Автореферат
диссертации на соискание ученой степени
доктора технических наук
Москва - 2010
Работа выполнена в Московском государственном университете инженерной экологии
Научный консультант: докторатехническиханаук, профессор
Володин Виктор Михайлович.
Официальные оппоненты:
докторатехническиханаук, профессор Юркевич Евгений Владимирович;
докторатехническиханаук, профессор Дорохов Игорь Николаевич;
докторатехническиханаук, профессор Большаков Александр Афанасьевич.
Ведущая организация: ФГУП Государственный ордена
Трудового Красного Знамени НИИ
химических реактивов и
особо чистых веществ, г. Москва.
Защита состоится л 20 мая 2010 г. в 14 часов на заседании диссертационного совета Д 212.145.02 при Московском государственном университете инженерной экологии по адресу: 105066, г.аМосква, ул. Старая Басманная, 21/4.
С диссертацией можно ознакомиться в библиотеке Московского государственного университета инженерной экологии.
Отзывы в двух экземплярах, заверенные печатью, просим направлять по адресу: 105066, г.аМосква, ул.аСтарая Басманная,а21/4.
Автореферат разослан л___ _________ 2010аг.
Ученый секретарь совета
д.т.н., проф. Латышенко К.П.
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность работы. В рамках общей концепции долгосрочного развития РФ происходит изменение объективных условий функционирования промышленных организаций, связанное с ростом их размеров и сложности, возникает необходимость модернизации процессов переработки сырья, оптимальной координации элементов производства, совершенствования структуры систем управления. Актуальным является комплексный анализ производственного цикла предприятий, научно обоснованный подход к построению интеллектуально ёмких систем автоматизированного управления в химической и смежных отраслях промышленности.
В диссертации систематизированы результаты теоретических и экспериментальных разработок автора в области исследования технологических систем высокой сложности в рамках распределённых автоматизированных систем управления (АСУ) с иерархической структурной организацией. В качестве объекта исследования принят ряд химико-технологических процессов многостадийного производства промышленного продукта из ископаемого сырья. При разработке и эксплуатации объектов выделенного класса Ц сложных технологических систем (СТС) возникает необходимость учёта их специфических свойств, существенное значение приобретают экономические факторы, повышенное внимание к экологии производства, требуется пересмотр методологии решения проблемных задач на основе комплексного подхода к созданию интегрированных АСУ. Возможность получения новых механизмов межуровневого взаимодействия в децентрализованных иерархических системах является основой создания единого информационного пространства предприятий, ведёт к интенсификации производства.
Исходные положения проблем управления структурно и технологически сложными системами изложены в работах Д. Мако, М.аМесаровича и др. Существенный вклад в развитие методов оптимизации сложных систем с использованием методов декомпозиции внесли Дж. Лэсдон, В.И. Цурков и др. Теоретические начала построения организационных систем опираются на информационную теорию иерархических систем и теорию кооперативных игр (Ю.Б. Гермейер, Н.Н. Моисеев и др.). Моделированию и оптимизации химико-технологических процессов (ХТП) посвящены работы В.В. Кафарова, В.М.аВолодина, В.С.аБалакирева, В.П.аМешалкина и др. Непрерывные процессы химической технологии отличаются высокой степенью автоматизации, однако, в условиях узкой специализации применяемого оборудования и высокой материало- и энергоёмкости, позволяют разрабатывать новые научные подходы к их интенсификации и автоматизации. В большей степени подобное характерно для не типовых химико-технологических процессов. Специфика моделирования и управления СТС до настоящего времени не выделена как самостоятельная сфера, преобладают подходы, базирующиеся на использовании традиционных АСУ, отличающихся разобщенностью отдельных уровней и контуров управления. Исследование возможностей и расширение областей применения декомпозиционного подхода и, в частности, методов явной декомпозиции, к синтезу открытых систем иерархически-взаимосвязанного управления производственными комплексами, к настоящему времени не нашли широкого применения во многих отраслях промышленности. Разработка научно-теоретических и методологических основ создания децентрализованных иерархических высокотехнологичных систем со сложным взаимодействием элементов требует дополнительного изучения и расширения сферы их использования на практике.
В основе подхода автора диссертационной работы к управлению СТС переработки природного сырья лежит обоснованное целенаправленное выделение подсистем, их анализ и моделирование, разработка методики эффективной реализации метода явной декомпозиции для синтеза взаимосогласованных управляющих подсистем. В диссертационной работе проводится исследование экологически направленных многостадийных производственных процессов очистки от примесей и производства активированных углей. Обобщённый научный подход к построению систем автоматизированного управления позволяет модернизировать процессы производства качественных продуктов с использованием новых информационных технологий, получить синергетический эффект на основе закономерностей взаимодействия элементов СТС. Задача поиска общих закономерностей при анализе и синтезе методов управления сложными технологическими процессами в химической и смежных отраслях промышленности - одна из наиболее актуальных задач научных исследований.
Диссертационная работа выполнялась в соответствии с федеральной целевой научно-технической программой УИсследования и разработки по приоритетным направлениям развития науки и техники на 2002 - 2006 годыФ по направлению УКаталитические системы и технологииФ; федеральной целевой программой УИнтеграция науки и высшего образования России на 2002 - 2006 годыФ (постановления Правительства РФ от 21.08.2001; от 05.09.2001).
Автор выражает благодарность проф. В.М.аВолодину за плодотворные научные консультации и содействие в проведении работы.
Целью диссертационной работы является интенсификация производственного цикла промышленных предприятий по переработке природного сырья, что предполагает решение важной научной проблемы создания методов и алгоритмов управления сложными технологическими процессами на основе иерархически-взаимосвязанных систем автоматизации, имеющей важное хозяйственное значение.
Названная цель определила постановку следующих задач:
- выявления способов повышения интенсификации систем управления многостадийными производственными процессами;
- учёта особенностей многостадийных процессов в химической технологии и ограниченной исходной информации при разработке методов построения их комплексных математических моделей;
- выбора существенных связей и факторов с помощью особенностей метода явной декомпозиции и формирования способов повышения интенсивности функционирования выделенных подсистем;
- обоснования условий интенсификации процессов переработки природного сырья, и модернизации многофункциональных схем получения высококачественного стандартизованного продукта;
- разработки и реализации алгоритмов управления технологическими объектами как составными частями производственного цикла промышленных предприятий, имеющими в своём составе сложные по функциональному назначению и комплексным параметрам структуры со специфическими свойствами;
- интеграции существующих систем автоматизированного управления с новыми решениями оптимизационных задач на основе использования предлагаемых средств измерений, управления и информационных технологий для внедрения в практику технологических предложений и разработанных принципов, методов и алгоритмов.
Методы исследования основаны на положениях теории моделирования и идентификации технологических процессов, теории выбора и принятия решений, общей теории систем и системного анализа в управлении, методах химической кибернетики.
Научная новизна работы состоит в создании методов и алгоритмов иерархического управления сложными технологическими системами и включает:
- выявление условий иерархичности, открытости, взаимосогласованности и избирательности, определяющих возможности систематизации и разработки способов автоматизированного управления сложными технологическими системами;
- обоснование уровней декомпозиции системы управления и реализацию итерационного многоуровневого процесса разбиения сложной системы для реализации управления с меньшими затратами и получения решений, близких к оптимальным;
- построение аппроксимационных моделей многостадийных химико-технологических процессов, позволяющих получать устойчивые решения в условиях ограниченной исходной информации;
- создание методов согласования решения локальных задач с внутриуровневыми конфликтами и оптимумом на границе области допустимых решений при введении аппроксимаций в ограничениях и целевых функциях;
- формирование требований к объектам управления в химической и смежных отраслях промышленности с помощью выявленных основных классификационных признаков сложных технологических систем.
Достоверность научных положений, рекомендаций и выводов обеспечивается использованием фундаментальных законов ведения физико-химических процессов, корректным применением математических методов анализа таких процессов и синтеза систем управления, результатами численных и экспериментальных исследований.
Практическая ценность работы определяется решением следующих задач:
- определены факторы образования технологического ресурса и обоснование возможностей интенсификации работы производственных систем в химико-технологических производствах;
- предложены алгоритмы унификации метода явной декомпозиции при управлении сложными технологическими системами с учётом выбранных связей, особенностей локальных задач и процедуры их взаимосогласованного решения;
- предложены математические модели, позволившие построить сложные системы управления активацией углеродосодержащих материалов и очисткой от примесей рудного сырья в условиях ограниченной выборки исходных данных;
- обосновано повышение эффективности управления очисткой рудного сырья с помощью математических моделей этих процессов, построенных по данным анализа многокомпанентных растворов;
- разработана и реализована система управления производством активированных углей в каскаде вращающихся печей с согласованием работы подсистем по нагрузке, позволяющая повысить эффективность управления;
- решена задача автоматизированного управления печами с двигающимся вертикальным слоем по материальным газовым и тепловым потокам на базе применения современных математических методов и информационных технологий;
- материалы диссертации использованы в учебном процессе при чтении лекций в МГУИЭ, при выполнении курсовых работ, дипломном проектировании, исследовательских программах.
Практическая ценность подтверждена внедрением разработанных методик и алгоритмов управления в ОАО Сорбент г.аПермь, ОАО Электростальский химико-механический завод, НПО Химавтоматика, использованием на опытном производстве Института проблем нанотехнологий РАЕН г.аПермь и применением при проектировании цеха активированных углей в Сирийской Арабской Республике.
На защиту выносятся следующие новые результаты, полученные автором:
- обобщенный подход, определяющий возможности интеграции и взаимного согласования функций подсистем, позволяющий интенсифицировать работу средств автоматизации управления сложными технологическими процессами переработки природного сырья;
- схема взаимосогласованного управления взаимосвязанными подсистемами с реализацией алгоритмов явной декомпозиции и аппроксимацией сложных алгоритмических зависимостей при оптимизации многостадийных процессов с однородными и неоднородными стадиями;
- метод построения математической модели сложного многомерного объекта управления с учётом технологических ограничений на число экспериментальных данных и высоких требований по качеству многоассортиментной продукции;
- методика управления многоэлементными объектами на основе анализа совокупности с выработкой управляющих решений, направленных на рост эффективности и увеличение срока службы технологического оборудования.
Апробация работы. Основные положения и результаты диссертационной работы докладывались и обсуждались на Всероссийских научно-технических конференциях Математические методы в химии (Тула, 1993) и Информационные системы и модели в научных исследованиях, промышленности и экологии (Тула, 2007); Международных конференциях: Методы кибернетики химико-технологических процессов (Москва, 1994), Математические методы в химии и химической технологии (Тверь, 1995), Инженерная защита окружающей среды (Москва, 2000, 2002), Математические методы в технике и технологиях (Смоленск, 2001); Системный анализ в проектировании и управлении (Москва, 2006), Индустриально-инновационное развитие - основа устойчивой экономики Казахстана (Шымкент, 2006); Международной конференции по химической технологии Химическая технология (Москва, 2007); Научно-технической конференции МГАХМ, (Москва, 1995); Международной студенческой школе-семинаре Новые информационные технологии (Москва, 2000).
Публикации. По результатам исследований, представленных в диссертационной работе, опубликовано 50 научных работ, 13 из которых в журналах, рекомендованных ВАК, две монографии. Автору принадлежат систематизация, разработка и практическое обоснование использования методов иерархического управления сложными промышленными объектами, 19 работ написано лично автором, остальные выполнены в соавторстве, многие с В.М.аВолодиным.
Объём работы. Диссертационная работа состоит из введения, восьми глав, заключения и списка литературы. Она содержит 291 страницу машинописного текста, 12 таблиц, 54 рисунка, библиографию из 265 наименований.
СОДЕРЖАНИЕ РАБОТЫ
Во введении основное внимание уделено проблемам управления структурно- и технологически сложными многостадийными промышленными процессами переработки природного сырья, отмечено влияние на становление теоретических основ методов управления информационной теории иерархических систем, химической кибернетики, системного анализа. Поставлена сложная и актуальная задача научных исследований, состоящая в поиске и анализе общих закономерностей в построении и функционировании взаимосогласованных систем иерархического управления. Выделены цели и задачи диссертационной работы, проанализированы способы их решения. Объектом исследования приняты экологически обоснованные СТС переработки природного сырья в химической промышленности с иерархической организацией структурных единиц. Особое место в ходе научного исследования отведено многостадийным процессам, как широко распространённым в высокотехнологичных отраслях промышленности.
В первой главе на основе анализа особенностей многостадийных процессов и использования промышленных средств автоматизации поставлена задача иерархически-взаимосвязанного управления СТС как определённым классом объектов управления. Отличие технологической системы от технической обусловлено комплексным подходом к разработке, производству, эксплуатации, утилизации и управлению всеми этапами производства продукции с учётом социальных, экономических, стратегических, экологических и других аспектов. СТС позиционирована с точки зрения способа организации, поэлементного состава, целей функционирования как объекта в целом, так и подсистем, математического описания, взаимных связей между элементами, наличия субъективных факторов. Особое внимание обращено на эффективное управление СТС, которое возможно лишь при разделении управляющих функций между подсистемами.
Принятые понятия и определения позволили выделить характерные черты СТС многостадийной переработки сырья:
- большие масштабы по числу элементов и объёму выполняемых функций, приводящие к иерархической организации;
- сложность поведения и многообразие связей между элементами;
- нерегулярные статистически распределённые во времени внешние возмущающие воздействия;
- наличие у подсистем конкурирующих целей функционирования;
- взаимосвязь и взаимовлияние технологических и экономических аспектов функционирования;
- большое число факторов, определяющих поведение объекта управления, и ограничений на допустимое функционирование;
- возникновение качественно новых свойств, отсутствующих у элементов и подсистем, обусловленных их количественным ростом.
Обобщённая задача оптимального управления поведением СТС
(1) |
где х, u, y - векторы возмущающих воздействий, управлений и выходов ТС; F - глобальная целевая функция СТС; U - множество допустимых состояний ТС; g, h - функции математической модели и ограничений; t - время.
Задача (1) для многостадийных процессов с учётом возникновения основных возмущающих воздействий при изменении показателей перерабатываемого сырья, неизменных на длительный срок, и кусочно-постоянной функции x(t) поставлена в виде
(2) |
где i,аNаЦаиндекс и число подсистем; fi Цалокальные целевые функции.
Особенности рассматриваемых СТС, высокая размерность и сложность задачи управления (2) затрудняют её решение в среде автоматизированной системы управления с соблюдением эксплуатационных требований относительной оперативности и точности, а, следовательно, эффективной организации управляющих функций. В связи с этим разработаны специальные подходы и методы решения аналогичных задач, в целях интеграции существующих и построения новых систем автоматизированного управления СТС.
Проанализированные задачи оптимального управления объектами различной степени сложности в качестве подзадач систем автоматизированного управления производственными комплексами позволили обобщить принципы управления. Специфика СТС делает неэффективными системы одноцелевого управления с централизованной структурой в связи с неопределённостью критериев, параметров моделей, ограничений и других характеристик. С другой стороны фрагментарное внедрение средств автоматизации и их использование для решения частных задач не приводит к возможному интегрированному эффекту в области управления производствами. Конкретика синтеза эффективных систем управления основана на наличии технологического резерва интенсификации, который подразумевает направленное воздействие на производственные процессы на физическом уровне, совершенствование структуры производимой продукции и используемого сырья, комплексную и гибкую автоматизацию производства, совершенствование организации и управления.
Научные исследования проведены автором с позиции определения способов интенсификации процессов, повышения эффективности систем автоматизированного управления. Во внимание приняты целостность и иерархическое строение открытых систем большой размерности, высокая сложность и многофункциональность поведения многообразных активных подсистем, имеющих и реализующих локальные целевые функции, обусловливающие эмерджентность.
Во второй главе по результатам исследований и данным вычислительных экспериментов обоснована целесообразность использования иерархического управления промышленными объектами класса СТС.
Анализ централизованных и распределённых систем позволил автору предложить иерархическую структуру системы управления производством углеродных сорбентов, представленную на рис. 1. Отдельные технологические операции и подпроцессы объединены в блоки обеспечения стадий: подготовки сырья, гранулирования, экструдирования, сушки, активации, очистки и утилизации, контроля и упаковки готовой продукции. В состав технологических блоков входят группы однотипных аппаратов. На верхнем уровне функционирует автоматизированная система на базе вычислительного комплекса для выбора номенклатуры адсорбентов, используемых видов сырья и полуфабрикатов. На среднем уровне - информационно-управляющая система согласования стадий технологического процесса. На нижнем уровне - системы автоматического контроля и регулирования параметров процесса.
Рис. 1. Иерархическая блочно-модульная система
управления производством активированных углей
Эта система создана с учётом локальных подсистем, функционирующих на промышленных объектах, и выполняет задачи их объединения в укрупнённые СТС и интенсификации управления.
Качество функционирования СТС определено по экономическим показателям с учётом принятых обозначений, располагая математическим описанием, целевыми функциями и ограничениями для подсистем, поставлена формализованная задача управления
, | (3) |
, | (4) |
где cij - матрица смежности.
Пространственная декомпозиция позволила представить объект управления в виде структурных элементов, связи между которыми установлены, исходя из способа разбиения СТС на N подсистем.
Задача (3) - (4) соответствует формулировке (2), имеет блочно-модульную структуру, аддитивная глобальная целевая функция имеет природу функционала и экономическую интерпретацию. Иерархически-взаимосогласованное управление в СТС основано на использовании методов декомпозиции для решения частных задач (3) - (4), которые не отождествлены с декомпозицией задач математического программирования и конструированием задач меньшей размерности, а представляют эффективные способы и алгоритмы решения задач оптимального управления.
Структурная декомпозиция предполагает постановку и решение многовариантной задачи разбиения на отдельные подсистемы, которая возникает многократно на разных уровнях иерархии. Анализ возможных по количеству подсистем β и структуре вариантов разбиения приводит к оптимизационной комбинаторной задаче. Поиск допустимого варианта разбиения для совокупности подсистем есть экстремальная по отношению к затратам I задача
. | (5) |
Задание не единственно возможного критерия (5), характеризующего качество разбиения СТС и эффективность реализуемой системы управления, затруднено в явном виде. Задача разбиения любой системы на элементы, с учётом всех факторов, включая технологические ограничения, является плохо выраженной и трудноразрешимой, что не позволяет типизировать её решение. Автором на основе вычислительного эксперимента предложен метод приближённых оценок времени вычисления целевой функции, который подтверждает эффект декомпозиции и целесообразность использования иерархического управления в СТС. В исследованиях принята оценка времени решения глобальной задачи симплекс-методом при условно заданных параметрах вычислительной процедуры
(6) |
где Ti - время решения локальных задач; cj - коэффициенты, характеризующие компьютер и метод решения; ns - число связей.
Проведённые исследования опытных систем переменного состава с обратными связями из однотипных элементов, рассмотренные варианты выделения подсистем, исследования критерия оптимальности - функции затрат (6) привели к следующим выводам. Декомпозиция с любым количеством подсистем эффективна при наличии в системе более 17 элементов и позволяет получить наилучшие результаты, начиная с отношения числа подсистем к общему количеству элементов , что является базой обоснования рациональных и эффективных способов разбиения сложных объектов на подсистемы. При этом анализ специфики стадий позволяет реализовать стратегию межуровневых взаимодействий.
В третьей главе рассмотрены теоретические аспекты декомпозиционных методов оптимизации и разработана методология решения сложной многомерной задачи на основе классификации объектов управления, согласовании решения совокупности частных задач, обеспечивающая снижение общих затрат на управление СТС.
Способы поиска новых технических решений предложены в работах Ф. Цвикки, А.И. Половинкина, Г.А. Кардашева, И.Н. Дорохова и др. Анализ существующих систем автоматизированного управления выделенным классом производств позволил автору конкретизировать основные условия их модернизации и интенсификации технологических процессов. Предложены условия иерархичности, открытости, избирательности и взаимосогласованности, в соответствии с которыми: система автоматизации строится как многоуровневая с интеграцией взаимосвязанных уровней, если СТС является открытой и имеет внутренний резерв изменения способов воздействия на её элементы и их отдельные характеристики, дибо дополняется таковыми; при построении систем автоматизации избирательно выделяются существенные связи между элементами и уровнями иерархии; элементы системы взаимосогласованно выполняют свои функции в рамках решения глобальной задачи управления.
Подход автора к решению задач управления основан на алгоритмах Данцига-Вульфа, Корнаи-Липтака и алгоритмах оптимизации сложных химико-технологических схем, предложенных в работах В.С.аБалакирева, В.М.аВолодина и др., которые предполагают итеративный межуровневый взаимообмен информацией между координирующим центром и локальными системами управления, решающими оптимизационные задачи на каждом шаге координации. Исследования имеют целью систематизацию методов декомпозиции, их анализ с целью определения возможностей и условий применимости. Основной проблемой при создании распределённых иерархических процедур решения сложных задач оптимизации является поиск способа декомпозиции на подзадачи и выбор методов, приводящих к приемлемому решению за время меньшее, чем без декомпозиции. Формальные методы для подобных процедур отработаны только для определённых классов задач. Характер информационного взаимообмена между уровнями иерархии определяется методом декомпозиции, а алгоритм взаимодействия между подсистемами и координирующим органом требует унифицированных подходов для выделенного класса сложных многостадийных процессов с однородными и неоднородными стадиями.
Методы декомпозиции применимы к оптимизационной сепарабельной задаче с ограничениями вида (3) - (4) для поиска экстремума по , где Ui - подмножества допустимых решений
(7) |
Явная декомпозиция предложена на основе сравнения с рядом методов, которые основаны на разбиении исходной задачи, как наиболее универсальный, простой и широко применимый в системах иерархического управления метод. Для координации локальных задач использованы значения формально разорванных входов и выходов подсистем СТС (рис.а2), заданные с соблюдением условия неразрывности (4). Обозначив аЦ значения переменных связи, задачу (3) будем трактовать как локальную на заданной области допустимых управлений для глобальной задачи координации
(8) |
Рис. 2. Схема реализации метода явной декомпозиции в СТС
Практика расчётов при реализации в технологических системах показывает эффективность метода явной декомпозиции в связи с возможностью применения для невыпуклых задач, простотой задания начального приближения вследствие конкретной физической природы параметров координации. При этом всякое полученное решение, возможно, не являясь оптимальным, физически реализуемо.
Специфика практических задач управления СТС требует учёта фактора времени, что затрудняет реализацию декомпозиционных методов. Согласно предложенному подходу декомпозиции во времени схема обобщённой задачи управления (рис. 3) предполагает двухуровневое принятие управляющих решений по оценке состояния системы в i-й ситуации - (СТС)i, i=1,2,ЕN.
Рис. 3. Структура СТС в методе декомпозиции во времени:
КО - координирующий орган, УС - управляющая система,
решающая модифицированную задачу
В случае учёта только типовых из бесконечного множества ситуаций без их комбинаций задачу координации для многостадийной очистки с критерием C0 (концентрация продукта) получим в виде
(9) |
Модифицированная частная задача управления k-й ситуацией
(10) |
Проведённые вычислительные эксперименты показали совпадение решения задач (9)аЦ (10) и (2) - максимизации конечной концентрации продукта.
Исходя из специфики проблем автоматизации многостадийных процессов для решения задач управления обосновано применение метода явной декомпозиции для широкого класса задач. В случае ограниченного числа однородных стадий возможно использование метода динамического программирования, декомпозиция во времени требует создания систем искусственного интеллекта и отнесена к частным задачам управления в определённом диапазоне ситуаций.
Четвертая глава посвящена обоснованию эффективных методов и алгоритмов локального управления подсистемами и специфических способов согласования локальных задач при реализации общих целей управления СТС, которые в рамках теории исследования операций сложно однозначно отнести к определённой модели принятия оптимальных решения.
Применение методов решения локальных задач в подсистемах и глобальной задачи согласования в первую очередь определено их спецификой. Задача координации допускает применение безградиентных методов поиска. Требование высокой скорости сходимости алгоритмов решения локальных задач даёт возможность использовать предложенный метод последовательной аппроксимации (МПА) для задачи многомерной условной оптимизации. МПА при высокой скорости обеспечивает относительную простоту вычислений, позволяет определять интервалы изменения переменных в практических задачах, что сокращает общее время решения. Линеаризация в окрестности точки x(k) приводит к задаче
(11) |
Обоснование выбора методов, классификация задач управления и формирование базы данных методов проведены на основе экспериментально установленной автором зависимости эффективности поиска экстремума от поведения характеристических функций
(12) |
Исследования позволили на основе анализа эффективности алгоритмов оптимизации (12) с использованием МПА, комплексного метода Бокса, метода скользящего допуска и комбинации методов с аппроксимационным сужением интервалов на основе (11) предложить алгоритм выбора метода, соответствующего характеру экстремальных зависимостей.
Модификация метода группового учёта аргументов (МГУА), предложенного А.Г.аИвахненко, использована для идентификации математических моделей. МГУА не имеет прямых рекомендаций по правилу разделения экспериментальных данных на последовательности и процентном соотношении точек в каждой из них. В ходе вычислительных экспериментов обнаружена неустойчивость алгоритмов МГУА при малом числе экспериментальных данных или сложной структуре частных полиномов. Исключить неустойчивость при высокой точности аппроксимации позволило использование неполного частного полинома
. | (13) |
Отношение количества точек в обучающей и проверочной последовательностях предложено по результатам эксперимента .
Оценку точности аппроксимирующего полинома предложено проводить по точкам контрольной последовательности МК (рис. 4) которые не участвуют в формировании полиномов и позволяют объективно оценить точность математического описания. Правило останова алгоритма записано в виде
или , | (14) |
где - ошибка на контрольной последовательности; - номер последнего уровня аппроксимации; - заданная точность.
Рис. 4. Алгоритм идентификации по модифицированному методу
Разработанный модифицированный метод использует МК для определения оптимального числа уровней аппроксимации и двухэтапное формирование частных полиномов, позволяет устранить явление неустойчивости и производить качественный анализ при ограниченной информации об объекте управления.
Следующие главы работы посвящены практической реализации предложенных способов, методов и алгоритмов управления промышленными объектами различного класса сложности.
В пятой главе с точки зрение химической кибернетики исследован производственный процесс тонкой очистки кобальтовых растворов от железа и меди по гидрометаллургической схеме на комбинате "Южуралникель" и приведены практические рекомендации по модернизации системы автоматизированного управления. Предложенная система позволяет реализовать управление отдельными стадиями процесса и обеспечить согласование работы всех звеньев технологического комплекса.
Гидрометаллургия как метод извлечения металлов имеет широкие возможности и, главным образом, экологические преимущества применения. Многостадийные процессы используют в технологии цветных и редких металлов для обеспечения очистки целевого продукта от сопутствующих примесей. Стадия тонкой очистки от железа и меди является одной из основных в процессе получения кобальт-никелевых растворов.
Исходя из общей цели, критерием оптимальности трёхстадийной очистки служит максимальный выход продукта реакций - кобальта из последнего реактора каскада аппаратов при соблюдении технологических ограничений на содержание железа и меди. Эффективность функционирования стадий - аддитивная зависимость
(15) |
где - концентрация кобальта в исходном растворе; - потери продукта на каждой стадии .
Постановка задачи оптимального управления процессом очистки в статическом режиме
(16) |
где qi, qci - расходы общий и кальцинированной соды; Cij - концентрация j-го компонента в i-м реакторе; V - объём реактора; rij - скорости реакций.
Выбор методов математического моделирования осуществлён на основе сравнения моделей, учитывающих сложную кинетику химических реакций, в пользу модифицированного метода с аппроксимационными подходами в условиях ограниченной выборки исходных данных по результатам фотометрического анализа состава сложных многокомпонентных растворов.
Постановка задачи (16) допускает решение методом динамического программирования в виде функционального уравнения
(17) |
Метод явной декомпозиции для задачи оптимального управления технологическим процессом очистки кобальтовых растворов от железа и меди с учётом формулировки задачи (16), декомпозиции моделей подсистем и выбранного критерия оптимальности запишем в виде (18)аЦа(19)
(18) |
где kj(pH) - зависимости констант скоростей реакций от кислотности.
С учётом решения локальных задач
(19) |
На основе предложенных в гл.а3 методики и алгоритма получено устойчивое решение задачи оптимального управления процессом тонкой очистки кобальтовых растворов. Определены допустимые значения управляющего воздействия на каждой стадии, позволяющие повысить выход кобальтового продукта. Получены зависимости содержания примесей и целевого продукта в реакторах каскада (рис. 5). Нестационарность первого реактора более выражена и область допустимых управлений имеет сложную структуру, поэтому решение задачи стабилизации его состояния на локальном уровне и определения оптимального расхода реагента наиболее важны. Отмечено, что при выборе оптимального управления на первой стадии удается найти управляющее воздействие для последующих стадий и получить концентрации меди и железа, отвечающие ограничениям задачи.
а | б |
Рис. 5. Изменение концентраций целевого продукта (а)
и примеси (б) по реакторам
Взаимосогласованное параллельное решение локальных задач с учётом их особенностей показывает преимущества предложенных методов, обусловливает возможность их применения без существенного усложнения решения при сложной структуре исходной системы, наличии рециклов. Процедура координации заключается в соблюдении неразрывности потоков по концентрации, которая возможна при управлении основным реагентом на каждой стадии.
В шестой главе с позиции определения технологического ресурса модернизации проанализирована отечественная технология производства гранулированных и дроблёных активных рекуперационных углей парогазовой активации на основе переработки слабоспекающегося каменного угля с целью выявления способов интенсификации производственных процессов.
Обобщённая схема производства адсорбентов (рис. 6) включает стадии предварительной подготовки сырья, карбонизации, активации, рассева и не учитывает многообразия углей от бурых до антрацитов. Основными предпосылками использования ископаемых углей для получения адсорбентов является широкий ассортимент их исходных свойств, сравнительная дешевизна и наличие сочетания известных стадий в схеме технологической переработки, наряду с важностью и определяющей ролью сложного процесса активации, обеспечивающего образование развитой пористой структуры.
Рис. 6. Схема технологического процесса изготовления
активированного угля на каменноугольной основе
Поглощающая способность активных углей обусловлена наличием выделенных в соответствии с классификацией академика М.М.аДубинина образуемых при активации пор, разделяемых по эквивалентному радиусу на микропоры, супермикропоры, мезопоры и макропоры. Для определения параметров микропористой структуры исследованы возможности газо-импульсной хроматографии, предпочтение этому методу отдано в силу его экспрессности и погрешности 2аЦа4а%. Получены результаты, позволяющие рассчитать характеристическую энергию адсорбции (E0) и предельный объём адсорбционного пространства (W0), отождествляемый с объёмом микропор как линейные регрессионные зависимости
(20) |
где - изостерическая теплота адсорбции как функция удельного удерживаемого объёма и обратной температуры.
Выбор диапазона режимов обработки природного сырья при фиксированной структуре технологического процесса позволяет выпускать разные марки продукта, а эффективное автоматизированное управление служит основой качества активированного угля.
Задачи управления СТС переработки сырья - многокритериальные с преобладанием экономической составляющей. Для их практического решения использована свёртка в виде суммы оценок альтернативы по каждому критерию с неравномерной оптимизацией частных критериев для неоднородных стадий (a1i, a2i - весовые коэффициенты). Аддитивный критерий эффективности
(21) |
Экономически обоснованным оборудованием для парогазовой активации каменного угля являются многоканальные подовые печи активации (МППА), решение задач автоматизированного управления в которых позволяет вывести процессы на новый качественный уровень. В технологических схемах каскадов печей перспективно использование вторичного тепла газов карбонизации и взаимосогласование стадий по нагрузке. В России с экономической и экологической точки зрения целесообразно использование каменного угля для производства адсорбентов. При этом его получение из воспроизводимого сырья (кокосовых орехов, древесины, плодовых косточек и др.), отходов нефтехимии и нефтепереработки и путём реактивации углей требует расширения сферы использования.
Технологические схемы переработки каменноугольного сырья не относятся к типовым в химической технологии, характеризуются низким уровнем автоматизации и отнесены к классу СТС, что делает возможным применение предложенных методов и алгоритмов управления. Результаты исследований нашли практическое применение при разработке и построении децентрализованных систем управления промышленными объектами рассматриваемого класса.
В седьмой главе на основе предложенных методов моделирования и оптимизации для каскада печей непрерывного действия реализованы управляющие функции в рамках системы иерархически-взаимосвязанного управления.
Общая методика построения адекватных моделей использована для каждой из подсистем СТС. Аналитическое описание сушки сыпучих дисперсных материалов в барабанных сушилках основано на теории конвективной сушки, но требует уточнения коэффициентов и параметров для конкретных процессов. Для печи карбонизации использованы модели вида (13). Особое внимание уделено моделированию активации. На основе теории исследования операций и мнения экспертов выделены факторы, влияющие на активацию: расход активирующего агента - x1; температурный режим - x2; расход материала на активацию - x3; частота вращения барабана печи - x4. Качественные характеристики готового продукта оценены по насыпной плотности - y1, механической прочности на истирание - y2 и сорбционной активности продукта - y3 в соответствии с классификацией продукта по маркам. Математическое описание процесса активации имеет вид
(22) |
К исследуемому объекту применены общие принципы иерархического управления, выделены подсистемы, сформулированы локальные задачи. Критерий оптимальности процесса сушки сформулирован как минимизация удельных затрат и отклонения влажности сушеного материала от заданной в сторону увеличения при известных значениях загрузки исходного материала и его влажности
(23) |
где a1i - весовые коэффициенты; qi, qiт, - расходы исходного материала, топлива, суммарный расход воздуха на i-й стадии; Hi - условно постоянные эксплуатационные расходы; w - влажность материала; T* - температурный режим (* - для всех стадий); i - скорость вращения барабана.
При эффективном использовании вторичного тепла q4 критерий оптимизации карбонизации - минимизация удельных производственных затрат
(24) |
где С - содержание летучих в карбонизованном материале.
Критерий оптимизации активации с учётом обозначений (22)
(25) |
где Pj - производительность печи для угля марки j=1,а2,а3; yj - технологические параметры, характеризующие сортность продукта; , k - насыпная плотность угля и карбонизата; Z1j - стоимость активированного угля; Zl - стоимость веществ l=1,Еа4; Q0 - неявные ограничения. Вектору входных воздействий соответствуют факторы x1а=аqa, x3а=аq2, вектору управлений - x2а=аT3, x4а=а3, выбранные как более значимые при оценке качественных характеристик готового продукта.
В рамках решения задачи методом явной декомпозиции глобальный критерий оптимизации технологического процесса производства углеродных сорбентов сформулирован как требование максимального значения техноэкономического показателя
(26) |
Алгоритм, предложенный автором в гл. 3, реализован в данной задаче управления. Оптимальное решение найдено в точке: расход активирующего агента qaа=а96акг/ч; температурный режим активации Та=а903,5аС; расход карбонизованного материала q2а=а20акг/ч; частота вращения барабана а=а57,6аоб/ч.
Результаты исследований позволили предложить перспективный вариант структуры СТС производства активированных углей - схему из двух последовательных технологических модулей (ТМ) (рис. 7). Первый модуль состоит из стадий сушки и карбонизации, связанных потоком рецикла по газам карбонизации, которые используют для приготовления сушильного агента, что даёт существенное снижение удельных затрат. Второй - из параллельно задействованных печей активации.
Рис. 7. Оптимальная структура ХТС производства адсорбентов
Экономическая эффективность стадий термообработки зависит от общей структуры системы и её производительности. Оптимальная структура этой системы и расчётные значения оптимального управления соответствуют производительности 1500аЦа2000акг/ч, что обусловливает внутриуровневые конфликты неоднородных стадий. Увеличение экономической эффективности системы не является беспредельным из-за ограниченной области существования технологических режимов на стадии активации, соответствующих определённой сортности продукта. При параллельном включении в технологический модуль активации четырёх вращающихся печей верхним пределом загрузки является значение 2500акг/ч. Дальнейшее повышение экономической эффективности может быть достигнуто при разработке более мощной вращающейся печи активации.
В восьмой главе систематизирована и решена задача управления процессом активации угольного сырья в МППА. Автоматизация отдельных стадий процесса, комплексный анализ и оптимизация по взаимосвязанным материальным, газовым и гидравлическим режимам реализованы в целях управления производственным циклом обработки ископаемого сырья. Разработана система автоматизированного управления, состоящая из подсистем: загрузки исходного сырья; регулирования температуры газов активации в боровах печи и регенераторах; регулирования скорости газов активации и управления выгрузкой готового продукта. В состав системы входят независимые контуры регулирования разрежения в дымовой камере и температуры в зоне активации, параметров, определяющих выходные характеристики продукта.
Исследование качества продукта по насыпной плотности показывает краевой эффект вследствие неоднородного теплового профиля, что ведёт к снижению скорости активации. Повышение эффективности работы многоканальной печи осуществлено за счёт выбора стратегии технического обслуживания и ремонта. Разработан метод, позволяющий по результатам качественных характеристик продукта определить время и место проведения ремонтных работ с целью увеличения сроков службы оборудования. Отличие предлагаемого подхода от исследования характеристик надёжности отдельных аппаратов и машин состоит в рассмотрении многоэлементного объекта в виде совокупности однотипного оборудования, в смысле потребляемого сырья и выпускаемого продукта. Каждый vi элемент множества единиц оборудования охарактеризован вектором i=1,ЕаN, где N - общая численность единиц оборудования; N(t) - численность совокупности в момент времени t; xi(t) - вектор показателей функционирования. В качестве показателя надёжности экспериментально проанализирована насыпная плотность угля, предложена математическая модель в виде интегро-дифференциальных уравнений. Основной характеристикой совокупности выбрана возрастная структура. Для описания обобщённых характеристик однотипного оборудования, имеющего неодинаковую производительность, использована плотность распределения численности по показателям функционирования , обладающая свойством эмержентности, присущим совокупности, что понижает размерность задачи и обобщает вид зависимостей взаимосвязи переменных. В качестве управляющих воздействий выбраны переменные, определяющие политику ремонта оборудования, изменяющие характеристики совокупности.
Критерий оптимальности - объём выпускаемой продукции
, | (27) |
где S - доход от производства единицы продукции.
Задача оптимального управления совокупностью состоит в определении закона изменения управляющего воздействия u(t), при котором система функционирует в рамках заданной производительности и максимизируется функционал (27). Результаты вычислительного эксперимента позволили оценить влияние закона изменения управляющего воздействия u(t) на эффективность управления при кусочно-линейном виде u(t) (рис. 8). Модельные зависимости - сплошные линии, область минимума на графике соответствуют оптимальному количеству продуктовых каналов.
а - ступенчатый вид u(t); | б - выбор оптимального u(t) |
Рис. 8. Зависимость численности продуктовых каланов от времени
Изменение численности продуктовых каналов Ns описывает реальную производственную ситуацию, т.к. производительность канала с достаточной степенью точности определяет основные показатели функционирования. Сформированная в результате временной декомпозиции подсистема выгрузки нацелена на выделение застойных зон, приводящих к зашлаковке продуктовых каналов, на повышение удобства контроля и регулирования количества и качества выгружаемого из печи активного угля.
Комплексная автоматизация с реализацией принципов управления многоэлементным объектом позволила интенсифицировать работу высокопроизводительного оборудования.
Основные результаты работы.
Решена важная научная проблема интенсификации перерабатывающих предприятий с использованием предложенных методов и алгоритмов управления сложными технологическими процессами, реализуемых иерархически-взаимосвязанными системами автоматизированного управления.
- Выделен класс сложных технологических систем переработки природного сырья, определяемый общностью постановок задач и методов иерархического управления.
- Разработаны новые подходы к определению факторов интенсификации открытых систем автоматизированного управления при условии избирательного выбора существенных взаимосвязей элементов и возможности качественного изменения их состояния.
- Предложен алгоритм децентрализованного решения исходной проблемы и схема унификации построения систем иерархического управления сложными многостадийными производствами, позволяющая выбрать метод явной декомпозиции и алгоритмы согласования задач, в которых присутствуют внутриуровневые конфликты.
- Создана и апробирована система методов, применение которой с использованием аппроксимационных подходов, позволяет получить решение нелинейных многокритериальных задач оптимизации двух видов: с большим числом переменных различного типа и сложных химико-технологических процессов.
- Доказано, что математические модели на основе модифицированного метода позволяют повысить устойчивость решений при описании многостадийных процессов с учётом свойств исходного сырья, качественных характеристик продукта, требований технологии, ремонтопригодности и формализованных функциональных ограничений.
- Показаны возможности интенсификации процесса очистки в каскаде реакторов на основе использования современных приборов анализа состава многокомпонентных смесей, математических моделей различной степени сложности и комплекса методов декомпозиционного управления с взаимосогласованным решением задач разного уровня иерархии.
- Предложена структура автоматизированной системы и разработаны эффективные методы управления сложной технологической системой производства активированного угля в каскаде вращающихся печей, согласованных по нагрузке, позволяющие получить значительную экономию при использовании вторичного тепла.
- Обоснована комплексная система мер по интенсификации автоматизированного управления многоканальной печью активации с управлением выделенными подсистемами по взаимосвязанным материальным, газовым и тепловым потокам. Эта система позволяет вывести выпуск активного угля на новый качественный уровень.
Разработанные на основе предложенных методов алгоритмы и программное обеспечение использованы для модернизации ряда промышленных предприятий и в учебном процессе.
Основное содержание диссертации опубликовано в 50 научных изданиях, в том числе:
Монографии:
- Мокрова Н.В. Моделирование и оптимизация химико-технологической системы производства активированных углей. ЦПермь: Пермский ЦНТИ, 2007. - 116ас.
- Володин В.М., Ху Вен-Цен, Мокрова Н.В. Децентрализованное управление сложными технологическими системами. Южно-Казахстанский государственный университет, - Шымкент, 2008. - 264ас.
Ведущие научные журналы:
- Мокрова Н.В., Володин В.М. Обоснование выбора методов решения задачи оптимального управления сложными процессами. //аВестник ТГТУ, 2006, том 12. №а1. С.а22 - 28.
- Володин В.М., Мокрова Н.В. Способы построения математической модели каскада реакторов. //аВестник ТГТУ, 2006, том 12. №а2А. С. 323 - 327.
- Мокрова Н.В., Володин В.М. Моделирование декомпозиционного управления многостадийными процессами. //аХимическое и нефтегазовое машиностроение. 2007. №а2. С. 17 - 19.
- Мокрова Н.В., Володин В.М. Система управления процессом очистки кобальтовых растворов. //аПриборы и автоматизация. 2007. №а3 (81). С. 15 - 17.
- Мокрова Н.В. Производство и применение активированных углей. // Уголь. 2007. №а9. С. 61 - 62.
- Мокрова Н.В. Задача оптимального управления производством активированных углей. // Уголь. 2007. №а7. С. 72 - 74.
- Мокрова Н.В., Володин В.М. Анализ задач управления сложными технологическими системами. //аХимическое и нефтегазовое машиностроение. 2008. №а2. С. 3 - 5.
- Мокрова Н.В., Володин В.М. Декомпозиционное управление процессом производства активированных углей. //аПриборы и автоматизация. 2008. №а3 (93). С. 24 - 28.
- Мокpова Н.В. Формирование оптимальной структуры пор активированных углей. //аХимическая технология. 2008. №а5. c.а221аЦ 225.
- Мокpова Н.В. Автоматизированное управление стадиями производства углеродных сорбентов в многоканальных печах активации. //аХимическая технология. 2008. №а11. c.а598 - 604.
- Мокрова Н.В. Координация решения локальных задач при управлении каскадом печей активации. //аАвтоматизация и современные технологии. 2009. №а7. С. 3 - 8.
- Мокрова Н.В. Методика оценки качества углеродных адсорбентов. //аВестник ТГТУ, 2009, том 15. №а4. С. 790 - 794.
- Мокрова Н.В. Выбор оптимального числа подсистем при реализации декомпозиционного управления. //аВестник СГТУ, 2009, №а4(43), выпуск 2. С. 205 - 207.
Сборники, материалы конференций:
- Володин В.М., Заев А.В., Мокрова Н.В. Задача оптимального управления процессом тонкой очистки кобальтовых растворов от железа и меди. В сб. Методы кибернетики химико-технологических процессов. Тез.докл. VI Международной конференции, - М:, 1994. С. 75.
- Володин В.М., Заев А.В., Мокрова Н.В. Идентификация математической модели процесса тонкой очистки кобальтовых растворов. В сб. Математические методы в химии и химической технологии. Тез.докл. IX Международной конференции, Тверь, 1995. С.а52.
- Мокрова Н.В., Умбетов У. Анализ динамических процессов в каскаде реакторов непрерывного действия. Математические методы в технике и технологиях ММТТ - 14: Сб.трудов Межд. науч. конф. Т 6. / Смоленский филиал Московского энергетич. инс-та. Смоленск, 2001. С. 234 - 235.
- Володин В.М., Мокрова Н.В. Использование метода явной декомпозиции для решения практических задач. Вестник МАСИ. Том 6, Часть I / МАСИ. - М., 2002. С. 100 - 107.
- Мокрова Н.В. Системный анализ при построении математической модели процесса тонкой медеочистки. Системный анализ в проектировании и управлении: Труды X Межд. науч.-практ. конф.Ч.3. 2006. С.а69 - 70.
- Володин В.М., Мокрова Н.В. Выбор критерия эффективности процесса очистки кобальтовых растворов. Вестник МАСИ. Том 9 Часть I / МАСИ, - М:, 2006. С. 70 - 72.
- Ху Вен-Цен., Володин В.М. Умбетов У. Мокрова Н.В. Метод последовательной аппроксимации при использовании декомпозиции в управлении многостадийными процессами. Материалы межд. науч. прак. конф.Индустриально-инновационное развитие - основа устойчивой экономики Казахстана. Т 1., ЮКГУ, Шымкент, 2006,аС.а387аЦа390.
- Мокрова Н.В. Технология получения углеродных сорбентов с развитой структурой пор. Химическая технология: Сб. тез. док. межд. конф. ХТТ07. Т.3. - М.: ЛЕНАНД, 2007. С. 119 - 121.
- Мокрова Н.В. Экологические аспекты производства и потребления активированных углей. Химическая технология: Сб. тез. док. межд. конф. по хим. технологии ХТТ07. Т.2. - М.: ЛЕНАНД, 2007. С. 289 - 291.
- Мокрова Н.В., Володин В.М. Математическое моделирование процесса активации во вращающихся печах. Труды IV Всерос. науч. конф. с межд. участием. Ч. 2: - Самара: СамГТУ, 2007. С. 116 - 118.
- Мокрова Н.В. Выбор критерия оптимальности в информационной системе производства активных углей. Информационные системы и модели в научных исследованиях, промышленности и экологии. Всероссийская науч. технич. конф. / Под общей редакцией д-ра техн. наук, проф. Панарина В.М. - Тула: Изд-во ТуГУ, 2007. С. 15 - 16.
- Мокрова Н.В., Володин В.М. Реализация метода явной декомпозиции при управлении иерархическими системами. Вестник МАСИ., - М:, 2007. С. 68 - 71.
- Мокрова Н.В. Задача определения параметров газокомпрессорной сети. Современные проблемы информатизации в проектировании и информационных системах: Сб. трудов. Вып. 13/ Под ред. д.т.н., проф. О.Я. Кравца. - Воронеж: Научная книга, 2008. С. 421 - 423.
- Мокрова Н.В., Володин В.М. Задача разбиения сложной технологической системы при реализации декомпозиционного управления. Вестник МАСИ., - М:, 2008. С. 53 - 58.
- Мокрова Н.В. Использование ресурсосберегающей технологии при управлении каскадом печей. Перспективные направления развития науки и технологий. Доклады Всероссийской науч. технич. конф. /Под общей редакцией чл.-кор. РАН Мешалкина В.П. - Тула: Изд-во ТуГУ, 2008. С. 169 - 170.
- Володин В.М., Мокрова Н.В. Перспективы производства и использования микропористых углеродных сорбентов. Труды межд. науч.-практ. конф. Химия в строительных материалах и материаловедение в ХХ веке, Т.4., ЮКГУ, - Шымкент: 2008. С. 82 - 84.
- МокровааН.В. Эффективностьадекомпозиции задачиауправ-ленияаСТС. Современныеапроблемыаинформатизацииавамоделирова-нии и социальных технологиях: Сб. Трудов. Вып. 14/ Под ред. д.т.н., проф. О.Я. Кравца. - Воронеж: Научная книга, 2009. С. 226 - 228.
- Мокрова Н.В. Построение многофакторной математической модели процесса активации. Математические методы в технике и технологиях. Сб.трудов XXII Межд. науч .конф.Т.а9. Секция 10а/апод ред. В.С. Балакирева.аЦаПсков. изд-во Псков. гос. политехн. ин-та, 2009. С. 91аЦа93.
Подписано в печать 02.03..2010. Зак. 27/д. Тир. 150. 2,0 печ. л.
Издательский центр МГУИЭ
105066, Москва, Старая Басманная ул., 21/4
Авторефераты по всем темам >> Авторефераты по техническим специальностям