Построение статической модели абсорбера. 2
КУРСОВАЯ РАБОТА
по дисциплине
«МОДЕЛИРОВАНИЕ
СИСТЕМ»
на тему
«Построение
статической модели
абсорбера»
Курсовая
работа 21 с., 4 рис., 8 таблиц,
5 источников, 4 прил.
МОДЕЛЬ,
АБСОРБЕР, РЕГРЕССИЯ, МЕТОД НАИМЕНЬШИХ
КВАДРАТОВ, АДЕКВАТНОСТЬ, КРИТЕРИЙ ФИШЕРА.
Объектом исследования абсорбер – аппарат для осушки газа.
В работе получена статическая модель абсорбера в виде зависимости расхода осушенного газа от температуры газа и расхода абсорбента.
.
Все расчеты, приведенные в работе, производились с использованием программных продуктов:
-МАТLAB 7.0.1
-Microsoft
Excel
СОДЕРЖАНИЕ
стр.
| ВВЕДЕНИЕ………………………………………………………… |
4 |
| 1. Осушка газа абсорбционным методом ……..…………..…... | 5 |
| 2. ОБРАБОТКА Результатов активного эксперимента...........… | 6 |
|
2.1. Планирование эксперимента ……............................ |
6 |
| 2.2.
Определение основных |
8 |
| 2.4. Проверка результатов измерений по критерию грубой ошибки…......... | 9 |
| 2.5. Определение дисперсии воспроизводимости по критерию Кохрена.... | 10 |
| 3. Построение математической модели абсорбера………… | 10 |
| 3.1. Определение порядка и расчет коэффициентов модели…….……….. | 10 |
|
13 |
|
3.2.1. Критерий Фишера.…………………………………… |
13 |
|
3.2.2. Корреляционная функция |
15 |
| ЗАКЛЮЧЕНИЕ…………………………………………………… |
16 |
| СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ……………………………. | 17 |
| ПРИЛОЖЕНИЕ А. Текст m-файла определение порядка и расчет коэффициентов модели и проверки их на адекватность……………………… | 18 |
| ПРИЛОЖЕНИЕ Б. Текст m-файла Корреляционная функция остатков……. | 19 |
| ПРИЛОЖЕНИЕ В. Критерий грубой ошибки………………………………… | 20 |
| ПРИЛОЖЕНИЕ
Г. Квантиль распределения Кохрена……………………… |
21 |
Моделирование – основной метод исследований во всех областях знаний, научно-обоснованный способ получения оценок параметров и изучения свойств технических систем, необходимых для принятия решений в различных сферах инженерной деятельности.
Моделирование помогает понять и упорядочить результаты эмпирических наблюдений, создать логический каркас научной теории, обнаружить внутренние связи и соотношения между результатами эксперимента. Особую важность задачи моделирования приобретают при рассмотрении сложных объектов, априорные сведения о которых либо отсутствуют, либо незначительны. Многие свойства, которые остаются неучтенными при моделировании, например, из-за того, что они неизвестны, могут в корне изменить картину результатов моделирования, и модель, таким образом, оказывается очень далека от оригинала.
В задачи моделирования относят - всесторонний обзор всех важных аспектов изучения систем посредством моделирования, в том числе программного обеспечения моделирования, проверок достоверности и адекватности модели, моделирования входных данных, генераторов случайных чисел, генерирования случайных величин и процессов, статистических планов и анализа моделируемых экспериментов, возможности моделирования производственных систем.
Целью данной курсовой работы
является построение адекватной модели
расхода осушенного газа, зависящей от
входных параметров расхода абсорбента
и расхода газа.
1.
Осушка газа абсорбционным
методом
Абсорбер
является многофункциональным агрегатом,
выполняющим функции
Влажный
газ поступает через штуцер входа
во входную сепарационную секцию,
где из газа частично выделяются капельная
жидкость и механические примеси. Отделение
жидкости происходит в сетчатых отбойниках
- отбойнике у штуцера входа газа
и расположенным над ним кольцевом отбойнике.
Выделившаяся из газа жидкость отводится
из кубовой части аппарата. Затем газ поступает
в массообменную секцию, где происходит
осушка его абсорбцией диэтиленгликолем
(ДЭГом), поступающим на верхнюю тарелку.
Осушенный газ поступает в секцию окончательной
очистки газа, где в фильтрпатронах происходит
укрупнение частиц ДЭГа, а в сепарационной
тарелке - окончательная очистка газа.
Осушенный газ через штуцер газа выводится
из аппарата. Контактная ступень абсорбера
характеризуется высокой эффективностью
массопередачи и низким гидравлическим
сопротивлением, что позволяет значительно
сократить высоту аппарата и уменьшить
его металлоемкость.
Многофункциональность и снижение высоты
позволяет значительно сократить занимаемые
оборудованием площади и снижает затраты
на монтаж оборудования, т. к. снижение
высоты абсорбера позволяет его размещать
на блок-понтонах.
Конструкция ступеней осушки позволяет
при изготовлении абсорбера легко раз-мещать
его в аппарате или заменять в период эксплуатации
без каких-либо дополнительных технологических
операций, например, без сварочных работ.
При эксплуатации аппарата в случае необходимости,
например, при значительном изменении
нагрузки, часть свободного сечения может
быть быстро и легко заглушена.
Рисунок 1. Абсорбер для осушки газа
2. ОБРАБОТКА Результатов активного эксперимента
2.1. Планирование эксперимента
Исходные данные:
Т – температура газа – 11 градусов;
L – расход абсорбента – варьируется от 30 до 70 м3/ч;
X – концентрация абсорбента – 47 кг/м³;
G – расход газа – варьируется от 10000 до 30000 m3/ч;
α – уровень значимости = 5.
Общее число проводимых экспериментов r = 25. На этом плане произвольно выбирается 3 точки, в которых будут проводиться параллельные опыты (отмечены знаком r в таблице 2.1).
С
целью получения
Таблица
2.1. Общее число экспериментов
| G
L |
10000 | 15000 | 20000 | 25000 | 30000 |
| 30 | |||||
| 40 | r | r | |||
| 50 | r | ||||
| 60 | |||||
| 70 |
Проводится три группы параллельных опытов для выбранных экспериментов по 10-15 опытов в каждой группе. При проведении параллельных опытов на вход объекта подают одинаковые комбинации входных параметров, что позволяет оценить воспроизводимость эксперимента. Значения входных параметров для параллельных опытов выбирать произвольно (таблица 2.2-2.4).
Таблица
2.2. Эксперимент №1
| № | Расход газа G, m3/ч |
Температура газа T, град. |
Расход абсорбента L, м3/ч |
Концентрация абсорбента X, кг/м3 |
Расход газа Y, м3/ч |
| 1 | 15000 | 11 | 40 | 47 | 34373 |
| 2 | 15000 | 11 | 40 | 47 | 46180 |
| 3 | 15000 | 11 | 40 | 47 | 47406 |
| 4 | 15000 | 11 | 40 | 47 | 35869 |
| 5 | 15000 | 11 | 40 | 47 | 42387 |
| 6 | 15000 | 11 | 40 | 47 | 38125 |
| 7 | 15000 | 11 | 40 | 47 | 44968 |
| 8 | 15000 | 11 | 40 | 47 | 44073 |
| 9 | 15000 | 11 | 40 | 47 | 46422 |
| 10 | 15000 | 11 | 40 | 47 | 34649 |
Таблица 2.3. Эксперимент №2
| № | Расход газа G, m3/ч |
Температура газа T, град. |
Расход абсорбента L, м3/ч |
Концентрация абсорбента X, кг/м3 |
Расход газа Y, м3/ч |
| 1 | 20000 | 11 | 50 | 47 | 73585 |
| 2 | 20000 | 11 | 50 | 47 | 73275 |
| 3 | 20000 | 11 | 50 | 47 | 78333 |
| 4 | 20000 | 11 | 50 | 47 | 74439 |
| 5 | 20000 | 11 | 50 | 47 | 78587 |
| 6 | 20000 | 11 | 50 | 47 | 68911 |
| 7 | 20000 | 11 | 50 | 47 | 80398 |
| 8 | 20000 | 11 | 50 | 47 | 82678 |
| 9 | 20000 | 11 | 50 | 47 | 59139 |
| 10 | 20000 | 11 | 50 | 47 | 58370 |
Таблица 2.4. Эксперимент №3
| № | Расход газа G, m3/ч |
Температура газа T, град. |
Расход абсорбента L, м3/ч |
Концентрация абсорбента X, кг/м3 |
Расход газа Y, м3/ч |
| 1 | 25000 | 11 | 40 | 47 | 70102 |
| 2 | 25000 | 11 | 40 | 47 | 67445 |
| 3 | 25000 | 11 | 40 | 47 | 82392 |
| 4 | 25000 | 11 | 40 | 47 | 76440 |
| 5 | 25000 | 11 | 40 | 47 | 79607 |
| 6 | 25000 | 11 | 40 | 47 | 70944 |
| 7 | 25000 | 11 | 40 | 47 | 57531 |
| 8 | 25000 | 11 | 40 | 47 | 69686 |
| 9 | 25000 | 11 | 40 | 47 | 70228 |
| 10 | 25000 | 11 | 40 | 47 | 70797 |
2.2. Определение основных статистических характеристик параллельных опытов
Для каждой группы параллельных опытов определяются следующие статистические характеристики:
Максимальное значение - ;
Минимальное значение - ;
Среднее значение - , где число опытов в данной группе (объём выборки);
Дисперсия ;
Среднее квадратичное отклонение .
Результаты вычислений статистических характеристик сведены в таблицу 2.5.
Таблица 2.5. Статистические характеристики
| Эксперимент №1 | Эксперимент №2 | Эксперимент №3 | |
| 47406 | 82678 | 82392 | |
| 34373 | 58370 | 57371 | |
| 41445,2 | 72771,5 | 71715,2 | |
| 26804725 | 70322237 | 47555786 | |
| 5177,328 | 8385,836 | 6896,07 |
2.3 Проверка результатов измерений по критерию грубой ошибки
Для
оценки выборочных данных по критерию
наличия грубой ошибки (R критерий) для
каждой выборки, полученной в результате
проведения параллельных опытов, вычисляются
величины:
Расчетные значение и сравниваются с . (приложение В). Табличное значение выбирают для уровня значимости α и числа степеней свободы , где - объем выборки.. Если
принимают, что отклонение (или ) определяется случайными явлениями и принадлежит данной генеральной совокупности, в противном случае результат отбрасывается как грубый промах и оценку по R-критерию повторяют для следующего (или ) с пересчитанными значениями основных статистических характеристик.
Эксперимента №1:
9
α=10
1,1513 1,366 =2,26
<
<
Эксперимента №2:
9
α=5
1,1813 1,717 =2,26
< <
Эксперимента №3:
9
α=5
1,577 2,0281 =2,26
< <
Во всех трех выборках Rmax < Rкр и . Принимаем, что отклонение Ymax определяется случайными явлениями и принадлежит данной генеральной совокупности.
2.4 Определение дисперсии воспроизводимости по критерию Кохрена
Дисперсию воспроизводимости ( ) определяют путём сравнения выборочных дисперсий для параллельных опытов. При одинаковом числе опытов во всех выборках для сравнения дисперсий пользуются критерием Кохрена.
Рассчитывается значение
и сравнивается с (приложение Г), выбранным для числа сравниваемых дисперсий и числа степеней свободы , с которым определена каждая дисперсия. Если , сравниваемые дисперсии можно считать однородными, следовательно их можно усреднить. А значит дисперсия воспроизводимости
G=0,6452, Gкр=0,6167
G > Gкр , различие между дисперсиями значительное, поэтому в качестве дисперсии воспроизводимости выбираем меньшую из сравниваемых дисперсий.
S²воспр=26804725
3. Построение математической модели абсорбера
3.1. Определение порядка и расчет коэффициентов модели
Процесс создания модели начинается с выбора типа модели и, как правило, на первом этапе останавливаются на линейном варианте в форме алгебраического многочлена:
где , - неизвестные коэффициенты модели,
G, L - варьируемые входные параметры объёкта.
Поиск
неизвестных коэффициентов
Таблица 3.1. Результаты активного эксперимента
| № | Расход газа
G, m3/ч |
Температура газа T, град. |
Расход абсорбента L, м3/ч |
Концентрация абсорбента X, кг/м3 |
Расход газа Y, м3/ч |
| 1 | 10000 | 11 | 30 | 47 | 19497 |
| 2 | 15000 | 11 | 40 | 47 | 21321 |
| 3 | 20000 | 11 | 50 | 47 | 26611 |
| 4 | 25000 | 11 | 60 | 47 | 22861 |
| 5 | 30000 | 11 | 70 | 47 | 26508 |
| 6 | 10000 | 11 | 30 | 47 | 20782 |
| 7 | 15000 | 11 | 40 | 47 | 18989 |
| 8 | 20000 | 11 | 50 | 47 | 19353 |
| 9 | 25000 | 11 | 60 | 47 | 25396 |
| 10 | 30000 | 11 | 70 | 47 | 20834 |
| 11 | 10000 | 11 | 30 | 47 | 26910 |
| 12 | 15000 | 11 | 40 | 47 | 22512 |
| 13 | 20000 | 11 | 50 | 47 | 19488 |
| 14 | 25000 | 11 | 60 | 47 | 24308 |
| 15 | 30000 | 11 | 70 | 47 | 20448 |
| 16 | 10000 | 11 | 30 | 47 | 22246 |
| 17 | 15000 | 11 | 40 | 47 | 19814 |
| 18 | 20000 | 11 | 50 | 47 | 20238 |
| 19 | 25000 | 11 | 60 | 47 | 19428 |
| 20 | 30000 | 11 | 70 | 47 | 26886 |
| 21 | 10000 | 11 | 30 | 47 | 19326 |
| 22 | 15000 | 11 | 40 | 47 | 26451 |
| 23 | 20000 | 11 | 50 | 47 | 22847 |
| 24 | 25000 | 11 | 60 | 47 | 22163 |
| 25 | 30000 | 11 | 70 | 47 | 24934 |
Рисунок
2. Поверхность расхода газа Y
при изменении расхода газа G и расхода
абсорбента L.
После вычисления неизвестных коэффициентов рассчитывают остаточную дисперсию:
где - экспериментальное значение выходного параметра для определенных входных сигналов;
- величина выходного параметра, расcчитанного по модели при тех же значениях входных сигналов;
- число коэффициентов в
- число экспериментальных
Таблица
3.2. Поиск неизвестных
| Вид модели | Коэффициенты | S2ост | ||||||
| b1 | b2 | b3 | b4 | b5 | b6 | |||
| b1G+b2L | -0.0039 | 2.0265 | - | - | - | - | 7.9186e+006 | |
| b1G+b2L+b3 | -0.0035 | 1.8132 | 2.8059 | - | - | - | 7.6359e+006 | |
| b1G^2+b2G+b3L+b4 | 0 | -0.0043 | 2.0505 | 2.6245 | - | - | 7.9314e+006 | |
| b1G^3+b2G^2+b3G+b4L+b5 | 0 | 0 | 0.0003 | 0.8159 | 3.2729 | - | 8.1859e+006 | |
| b1G^4+b2G^3+b3G^2+b4G+b5L+b6 | 0 | 0 | 0.3826 | 0 | 0 | 0 | 1.8461e+007 | |
3.2. Проверка модели на адекватность
3.2.1. Критерий Фишера
Значения коэффициентов в уравнении регрессии, полученные по методу наименьших квадратов, являются оптимальными для выбранной математической модели, однако не всегда корректно останавливаться на этом. Процесс создания модели должен заканчиваться объективной оценкой, насколько точно построенная модель описывает идентифицируемый объект.
Проверка
модели на адекватность производится
путём сравнения суммы
Если ., то с достоверностью в (1 - a)*100% считают модель адекватной объекту, если нет - с той же достоверностью вероятно противоположное утверждение.
В случае, когда адекватность модели не подтвердилась, необходимо вернуться к началу п.2.2 и изменить вид модели. Чаще всего в такой ситуации просто увеличивают порядок модели и весь последующий расчет повторяют. Данные сведены в таблицу 3.3.
Модель адекватна.
Таблица 3.3. Сравнение F
| Вид модели | S2ост | F | Fкр | |
| b1G+b2L | 7.9186e+006 | 3.3850 | > | 2.9169 |
| b1G+b2L+b3 | 7.6359e+006 | 3.5104 | > | 2.9263 |
| b1G^2+b2G+b3L+b4 | 7.9314e+006 | 3.3796 | > | 2.9365 |
| b1G^3+b2G^2+b3G+b4L+b5 | 8.1859e+006 | 3.2745 | > | 2.9477 |
| b1G^4+b2G^2+b3G^2+b4G+b5L+b6 | 1.8461e+007 | 1.4519 | < | 2.9600 |