Научная статья на тему 'Моделирование процесса снижения загрязнения атмосферы выбросами оксидов серы от судовых энергетических установок'

Моделирование процесса снижения загрязнения атмосферы выбросами оксидов серы от судовых энергетических установок Текст научной статьи по специальности «Математика»

CC BY
307
57
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МОРСКИЕ СУДА / ВЫБРОСЫ / ТОКСИЧНОСТЬ / АТМОСФЕРА / ЗАГРЯЗНЯЮЩИЕ ВЕЩЕСТВА / ОКСИД СЕРЫ / ТОПЛИВО / СПОСОБЫ ОЧИСТКИ / EMISSIONS / TOXICITY / ATMOSPHERE / SULFUR OXIDE / FUEL / CLEANING METHOD / SHIPS / POLLUTANT

Аннотация научной статьи по математике, автор научной работы — Вострикова М. А.

В статье рассматривается проблема снижения выбросов оксидов серы от судовых энергетических установок путем применения установки поглощения газовых выбросов (УПГВ) с двухступенчатой очисткой. Приводится полученная при помощи метода планирования эксперимента математическая модель.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Вострикова М. А.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

The article deals with the problem of reducing sulfur oxide emissions from ship power plants through the use of the installation of the absorption of gas emissions with two-stage cleaning. Provides obtained using the method of experiment planning mathematical model.

Текст научной работы на тему «Моделирование процесса снижения загрязнения атмосферы выбросами оксидов серы от судовых энергетических установок»

ВЕСТНИКА

ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА^

6. ГОСТ Р 51249-99. Государственный стандарт Российской Федерации. Двигатели внутреннего сгорания поршневые. Выбросы вредных веществ с отработавшими газами. Нормы и методы определения / Госстандарт России. — М.: ИПК Изд-во стандартов, 2005. — 17 с.

7. ГОСТ 30574-98. Межгосударственный стандарт. Дизели судовые тепловозные и промышленные. Измерение выбросов вредных веществ с отработавшими газами. Цикл испытаний. — Минск: ИПК Изд-во стандартов, 1999. — 90 с.

8. ГОСТ 10448-80. Межгосударственный стандарт. Дизели судовые, тепловозные и промышленные. Приемка. Методы испытаний. — М.: ИПК Изд-во стандартов, 2003. — 38 с.

9. Зиненко Н. Н. Решение проблем охраны окружающей среды в компании «СКФ Новошип» / Н. Н. Зиненко // Материалы IX науч.-техн. конф., 17-18 декабря 2010 г. / МГА им. адмирала Ф. Ф. Ушакова. — Новороссийск,2011.

10. Пат. 110068 РФ. Система мониторинга энергоэффективности и экологической безопасности судов (S3ES-Novoship) / Зиненко Н. Н., Пруцков А. Г., Бордунов Б. В., Букаренко Ю. Г., Мартынович Е. С., Яременко Е. В., Панамарев В. Е.; патентообладатель ОАО «Новошип». — 2011.

УДК 502.3 М. А. Вострикова,

канд. техн. наук, ФГБОУ ВПО «Краснодарский государственный университет культуры и искусств»

МОДЕЛИРОВАНИЕ ПРОЦЕССА СНИЖЕНИЯ ЗАГРЯЗНЕНИЯ АТМОСФЕРЫ ВЫБРОСАМИ ОКСИДОВ СЕРЫ ОТ СУДОВЫХ ЭНЕРГЕТИЧЕСКИХ УСТАНОВОК

MODELING REDUCTION OF AIR POLLUTION EMISSIONS OF SULFUR OXIDES FROM SHIPS POWER PLANTS

В статье рассматривается проблема снижения выбросов оксидов серы от судовых энергетических установок путем применения установки поглощения газовых выбросов (УПГВ) с двухступенчатой очисткой. Приводится полученная при помощи метода планирования эксперимента математическая модель.

The article deals with the problem of reducing sulfur oxide emissions from ship power plants through the use of the installation of the absorption of gas emissions with two-stage cleaning. Provides obtained using the method of experiment planning mathematical model.

Ключевые слова: морские суда, выбросы, токсичность, атмосфера, загрязняющие вещества, оксид серы, топливо, способы очистки.

Key words: ships, emissions, toxicity, atmosphere,pollutant, sulfur oxide,fuel, cleaning method.

ИАСТУЩИЕ масштабы распространения загрязняющих веществ над океаном в результате сжигания топлива, особенно эмиссия в атмосферу оксидов серы, вызывают все большие опасения, и данная проблема усугубляется по мере роста мирового флота и растущих глобальных масштабов потребления бункерного топлива [1].

Одним из перспективных способов очистки продуктов сгорания от выбросов соединений серы с судовых энергетических установок является применение установок поглощения газовых

Выпуск 5

Выпуск 5

выбросов (УПГВ). Устройство представляет собой двухступенчатое смешивание газа с жидкостью и используется для очистки газовых потоков [2].

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

Принцип действия УПГВ 5 заключается в следующем. Продукты сгорания поступают в бакоснование через патрубок в верхней его части и, контактируя с водой, предварительно заливаемой в бак, несколько охлаждаются и очищаются. Далее газы направляются в блоки для фильтров, заполненных кассетами. Одновременно в работу включается система орошения, представляющая собой вибратор и душирующее устройство, в которое вода подается насосом через задвижку из бака-основания. Вибрация передается коробам и кассетам. Таким образом, продукты сгорания проходят через смоченные и вибрирующие кассеты противотоком навстречу стекающей воде. Кассеты представляют собой многослойные (до 10 слоев) сетчатые фильтры. Образующиеся на этих фильтрах пленки жидкости под действием вибрации обеспечивают процесс вибротурбулизации. Известно, что в этом процессе растворимость газов в воде существенно увеличивается [4]. После коробов остаточные газы дымососом откачиваются из отсека и направляются в дымовую трубу. Анализ дымовых газов осуществляется перед входом продуктов сгорания в бак-основание и перед дымососом. Вода циркулирует по замкнутому контуру до насыщения растворившимися оксидами. Контроль этого процесса выполняется по величине pH воды. При достижении рН6 часть воды сбрасывается в дренаж и на ее место доливается свежая вода. Отбор проб воды осуществляется через штуцер. Вибратор представляет собой электродвигатель, вал которого выходит за обмотки статора. На обоих концах вала установлены дебалансы — грузики, смещаемые относительно центра вала. Изменение положения дебалансов позволяет задавать амплитуду колебаний, от которой существенно зависит процесс поглощения оксидов. В реальных условиях была смонтирована УПГВ 6, аналогичная по конструкции УПГВ 5, но оборудованная предварительной ступенью очистки продуктов сгорания в виде струйного аппарата (рис. 1).

Для определения эффективности очистки продуктов сгорания от оксидов серы в струйном аппарате были проведены эксперименты на установке поглощения газовых выбросов УПГВ 6, оборудованной предварительной ступенью очистки продуктов сгорания в виде струйного аппарата (рис. 1). Струйный аппарат создан на базе результатов лабораторных исследований и использования геометрического и газодинамического подобия по числу Рейнольдса (Re) для модели и натурного образца с использованием положений [5].

Расход дымовых газов определялся по расходу топлива и материальному балансу при измерении температур и расходов орошающей жидкости и продуктов сгорания на входе и выходе УПГВ 6. Для анализа продуктов сгорания использовался газоанализатор типа IMR 3000 P.

Результаты эксперимента показали, что в струйном аппарате при смешении продуктов сгорания и водяного пара, который направлялся в приемную камеру струйного аппарата из рядом расположенного парового котла в количестве, соответствующем коэффициенту инжекции, процент поглощения оксидов составил по S02 около 68 %. После струйного аппарата газовый поток направлялся в УПГВ 6, в которой поглотилось: 75 % S02 от остаточной после струйного аппарата концентрации. Общий процент поглощения оксидов по предлагаемой схеме очистки продуктов сгорания достигает 92 % по S02, что позволит обеспечить выполнение требований Приложения VI Международной конвенции по предотвращению загрязнения с судов (MARPOL 73/78).

Рис. 1. Схема присоединения струйного аппарата к УПГВ 6:

СА — струйный аппарат;

ГА1, ГА2, ГАЗ — места отбора проб продуктов сгорания на газоанализатор IMR-3000P

Рис. 2. Изменение концентрации S02 по газовому тракту в процессе испытаний УПГВ 6: ряд 1 — точка замера ГА1 (до струйного аппарата); ряд 2 — точка замера ГА2 (после струйного аппарата); ряд 3 — точка замера ГАЗ (после УПГВ 6)

Эксперименты были проведены по специальному плану с целью получения математической модели, которая описывает влияние расхода активного потока модельного газа и исходной концентрации в нем диоксида серы на конечную концентрацию диоксида серы. Разработка этого плана осуществлялась с учетом рекомендаций теории планирования эксперимента [3].

Выпуск 5

Выпуск 5

ЕЬВЕСТНИК

ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

„МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

Математическая модель, представляющая собой уравнение регрессии, задаваемое полиномом второй степени, в общем случае имеет следующий вид:

к к

У = bo + X bX + Z b.jXxj + X b..xf, (!)

i=1 i< j i=1

где у — исследуемая величина (концентрация диоксида серы на выходе из струйного аппарата), мг/м3; х., х — значения факторов в кодированном масштабе; b0, b., b Ь.. — оценки коэффициентов уравнения регрессии.

Для отыскания неизвестных оценок bQ, b., b Ь.. применяется метод наименьших квадратов с использованием зависимостей регрессионного анализа. При этом приняты следующие допущения: независимые факторы х., х заданы абсолютно точно и являются неслучайными величинами. Также предполагается, что случайная величина у распределена по нормальному закону, а ее дисперсии в отдельных точках факторного пространства однородны.

Для определения оценок коэффициентов уравнения (1) используется план по схеме полного факторного эксперимента. При этом каждый из двух независимых факторов (zx — концентрация диоксида серы на входе в струйный аппарат мг/м3; z2 — расход активного потока, м3/с) варьировался на трех уровнях: +1 — высокое значение фактора; 0 — среднее значение фактора; -1 — низкое значение фактора.

Матрица планирования эксперимента, значения факторов в натуральном масштабе и в безразмерной системе координат, а также результаты эксперимента и рассчитанные значения исследуемой величины представлены в табл. 1.

Таблица 1

Матрица планирования и реализация полного факторного эксперимента

№ опыта Факторы в натуральном масштабе Факторы в безразмерной системе координат Исследуемая величина

Из эксперимента Из полного уравнения Из упрощенного уравнения

zi Z2 xi Х2 х,3 - 2 1 3 х22 - -2 3 У У У -^упр

1 200 0,0020 -1 -1 1/3 1/3 11,6 11,00 10,96

2 200 0,0025 -1 0 1/3 -2/3 16,4 17,61 17,57

3 200 0,0030 -1 +1 1/3 1/3 26,9 26,30 26,26

4 400 0,0020 0 -1 -2/3 1/3 22,3 22,31 22,31

5 400 0,0025 0 0 -2/3 -2/3 31,34 31,30 31,30

6 400 0,0030 0 +1 -2/3 1/3 42,4 42,37 42,37

7 600 0,0020 +1 -1 1/3 1/3 33,1 34,70 34,66

8 600 0,0025 +1 0 1/3 -2/3 46,3 45,07 45,03

9 600 0,0030 +1 +1 1/3 1/3 57,9 58,52 58,48

С целью приведения матрицы планирования эксперимента к ортогональному плану квадратичные столбцы х2. были преобразованы с использованием следующего линейного преобразования [3]:

n

xJ = xJ -xj = xJ -У 4 /n = 4 -6/9 = 4 -2/3, (2)

J J J J J1 J J

i=\

trqN= 9 — количество опытов в эксперименте.

Так как матрица планирования является ортогональной, то все коэффициенты уравнения регрессии (1) определяются независимо друг от друга по формуле

N N

bj = Z XjiVi 1 Z 4 ■

(3)

i=1

i =1

Например, для определения коэффициента bx при хх необходимо получить сумму произведений столбца cijB безразмерной системе координат на столбец су. Затем находим значение Ьх.

9 9

b=Z Wi / Z Xi=82>4/6=13,73.

i=1 i=1

Аналогично были определены эффекты единичных, двойных и квадратичных взаимодействий: b'0 = 32,02; 62 = 10,03; 612 = 2,38; Ьи = 0,04; Ь22 = 1,04.

Дисперсии коэффициентов уравнения (1) определяются по формуле

2 2 Sbi = S

N

/ У;

воспр ' " JI ■

i =1

(4)

В результате расчетов по матрице с преобразованными столбцами для квадратичных эффектов получается уравнение следующего вида:

У b0 + b1x1 + b2x2 + b12x1x2 + b11 (x1 X1 ) + b22 (x2 x2 )• Для перехода к обычной записи &0 определяется по формуле

Ь0 = Ь0-Ъхххх2 -b22X = 32,02- 0,04(2/3) -1,04(2/3) = 31,3 и оценивается с дисперсией, равной

2 2 . /—2 . / — 2\2 2 Ч = + (Х2) S,22-

2 , /—2x2 2

2 2 2

(5)

(6)

(7)

С целью определения дисперсии воспроизводимости ^оспр, знание которой необходимо для проверки значимости коэффициентов уравнения регрессии, в центре плана были поставлены дополнительно три параллельных опыта и получены следующие значенияу:

3

Л0 = 39,37; у°2 = 39,22; у3> = 39,85; у0 = Уу°и/3 = 39,48;

3

= Z (У0 - У2)2 / 2 = 0,108; s2 = 0,329.

воспр,^^ х*7 и у у 7 7 воспр ’

и =1

Затем с использованием выражений (4) и (7) были оценены дисперсии всех коэффициентов уравнения регрессии, которые оказались равными:

< = °,°6; 4 = < = 0,018; s22 = 0,027; ^ = <2 = 0,054.

Значимость всех коэффициентов уравнения регрессии оценивается по критерию Стьюдента.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Для этого по формуле tj = |bj| / рассчитывается 1-отношение:

t0 = 127,8; tx = 102,5; t2 = 74,9; t12 = 14,5; txl = 0,2; t22 = 4,5.

Табличное значение критерия Стьюдента для уровня значимости р = 0,05 и числа степеней свободы/= 2 tQ05(2) = 4,3. Так как t-отношение для коэффициента b оказалось меньше табличного значения, то коэффициент является незначимым и его можно исключать из уравнения регрессии. Все остальные коэффициенты являются значимыми и их нельзя исключать из уравнения. После исключения незначимого коэффициента уравнение регрессии будет иметь следующий вид:

у = 31,3 + 13,73xj + 10,03x2 + 2,38x1x2 + 1,04x2. Проверим адекватность полученного уравнения по критерию Фишера:

(8)

2

и

Выпуск 5

Выпуск 5

F = s2 / s2 ,

ост воспр’

здесь

sL =ZCv,-hf /(^-0 = 0,0096/4 = 0,0024,

<=i

где Z — число значимых коэффициентов в уравнении регрессии, равное 5.

Тогда F = 0,0024/0,108 = 0,022. Табулированное значение критерия Фишера для р = 0,05, f = 5,f = 2, F= 19,3. Так как F < F (fl,f2), то полученное уравнение регрессии (8) адекватно описывает эксперимент.

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

Список литературы

1. Туркин А. В. Снижение загрязнения атмосферы очисткой продуктов сгорания судовых и корабельных энергетических установок / А. В. Туркин // Материалы VII регион, науч.-техн. конф. —Новороссийск: МГА, 2008.

2. Комиссаров К. Б. Комплексная очистка дымовых газов теплогенерирующих установок: моногр. / К. Б. Комиссаров, С. А. Лутков, А. В. Филь — Ростов н/Д: Филиал МГА имени адмирала Ф. Ф. Ушакова, 2007. — 134 с.

Ъ.Ахназарова С. Л. Оптимизация эксперимента в химии и химической технологии / С. Л. Ах-назарова, В. В. Кафаров. — М.: Высш. Шк., 1978. — 319 с.

4. Рамм В. М. Абсорбция газов / В. М. Рамм. — М.: Химия, 1976. — 655 с.

5. Дьяконов Г. К. Вопросы теории подобия в области физико-химических процессов / Г. К. Дьяконов. — М.: Изд-во АН СССР, 1956. — 349 с.

УДК 502.7:627.215.2 Е. Г. Трунин,

канд. экон. наук, ФАУ «Российский речной регистр»

ПРОГНОЗИРОВАНИЕ ЗАГРЯЗНЕНИЯ ОКРУЖАЮЩЕЙ СРЕДЫ ПРИ УПРАВЛЕНИИ ПРИРОДОПОЛЬЗОВАНИЕМ В РЕЧНЫХ ПОРТАХ

PREDICTION OF ENVIRONMENTAL POLLUTION WHEN ENVIRONMENTAL

MANAGEMENT IN RIVER PORTS

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

In article considers the question of short-term forecasting of environmental pollution. The possibility to apply for a decision on this question of the apparatus of artificial neural networks is shown.

i Надоели баннеры? Вы всегда можете отключить рекламу.