УДК 628.353
А.С. Васильева, В.К. Евтеев A.S. Vasilyeva, V.K. Yevteyev
ТЕОРЕТИЧЕСКИЕ ОСНОВЫ И ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ АНАЭРОБНОГО ФИЛЬТРА
THEORETICAL BASIS AND EXPERIMENTAL RESEARCH OF ANAEROBIC FILTER
Ключевые слова: анаэробное сбраживание, анаэробный фильтр, биогаз, циклические возмущения, модель, планирование эксперимента, демонстрационная установка, параметр, фактор, эксперимент, регрессионные уравнения.
Для переработки жидкого навоза, в особенности животноводческих стоков, применяют метан-тенки с фиксированной биомассой — анаэробные фильтры. Приводится теоретическое исследование фазы опорожнения анаэробного фильтра для сбраживания животноводческих стоков. В данной фазе работы происходят циклические возмущения: изменение скорости движения субстрата в фильтре, обновление граничных поверхностей, лавинообразное выделение биогаза, разрушение пенного слоя. Составлена гидравлическая модель анаэробного фильтра в фазе опорожнения, найдена функциональная зависимость скорости движения свободной поверхности стоков в фильтре. В результате проведенного теоретического исследования выявлены основные воздействующие факторы процесса анаэробного сбраживания. С целью визуального и приборного наблюдения циклических возмущений, ранжирования воздействующих факторов проведено экспериментальное исследование на демонстрационной моделирующей установке. Она позволяет смоделировать процесс работы анаэробного фильтра с циклическим возмущением, с имитацией выделения биогаза. В качестве субстрата применялась вода, а биогаза — атмосферный воздух. Исследования проведены в соответствии с требованиями планирования многофакторного эксперимента, описана методика проведения опытов. Результаты эксперимента обработаны с помощью программы STATISTICA 6.0, получены регрессионные уравнения оценочных параметров. На основании полученных результатов экспериментальных исследований сформулированы следующие выводы: разработанная гидравлическая модель работы анаэробного фильтра с циклическим возмущением подтвердилась экспериментальными исследованиями на демонстрационной установке; в соответствии с гидравлической моделью скорость движения свободной поверхности стоков в фильтре зависит от физико-химических характеристик субстрата и технологических параметров
анаэробного сбраживания (Дз; Т; W; K), а также конструктивных особенностей установки (ф; d6;
dф; lc; dc; ZZC; H).
Keywords: anaerobic digestion, anaerobic filter, biogas, cyclic perturbation, model, plan of experiment, a demonstration plant, parameter, factor, experiment, regression equations.
For processing liquid manure, and especially, livestock wastewater digesters are used with fixed biomass - anaerobic filters. The article provides a theoretical study of emptying phase of the anaerobic filter for livestock waste water digestion. In this phase the cyclic disturbances occur: produces a change in speed of the substrate in the filter, the updating of boundary surfaces, avalanche allocation of biogas, and destruction of the foam layer. Hydraulic model of anaerobic filter is composed in emptying phase, functional dependence of motion speed of effluent free surface in the filter is found. As a result of the conducted theoretical study, the main factors affecting the process of anaerobic digestion are revealed. For the purpose of visual and instrumental observations of cyclical disturbances, the ranking of influencing factors, experimental study has been made on the modeling demonstration plant. It enables to simulate the operational process of anaerobic filter with siphon drainage, with simulation of biogas release. As a substrate the water was used and as a biogas - ambient air. The studies have been conducted in accordance with the requirements of planning multifactorial experiment, also the methodology of the experimental process are described. The experimental results are processed with STATIS-TICA 6.0 program; regression equations of evaluation parameters have been obtained. Basing on the achieved results of experimental studies the following conclusions are made: Elaborated hydraulic model of the anaerobic filter operation has been confirmed by experimental studies on the demonstration plant. -According to the hydraulic model, the velocity of free surface movement of the effluents in the filter depends on physical and chemical characteristics of the substrate and technological parameters of anaerobic digestion (Ds; T; W; Ch), as well as the constructive features of the plant
(lf;dp;df;ls;ds;Yl^s;p).
\
Васильева Аяна Сергеевна, аспирант, каф. технического обеспечения АПК, Иркутская государственная сельскохозяйственная академия. Е-шзИ: fvasiljiev@yandex.ru.
Евтеев Виктор Константинович, к.т.н., доцент, в.н.с., каф. технического обеспечения АПК, Иркутская государственная сельскохозяйственная академия. Е-mail: fvasiljiev@yandex.ru.
Vasilyeva Ayana Sergeyevna, Post-Graduate Student, Chair of Technical Support of Agriculture, Irkutsk State Agricultural Academy. E-mail: fvasil-jiev@yandex.ru.
Yevteyev Viktor Konstantinovich, Cand. Tech. Sci., Assoc. Prof., Leading Staff Scientist, Chair of Technical Support of Agriculture, Irkutsk State Agricultural Academy. E-mail: fvasiljiev@yandex.ru.
Введение
В сельскохозяйственном производстве Восточной Сибири существует проблема переработки навоза и навозных стоков свиноводческих ферм и комплексов. В условиях промышленного свиноводства нашли широкое применение бесподстилочное содержание животных, гидравлические и механические системы удаления навоза из животноводческих помещений. В ветеринарно-санитарном отношении наиболее приемлемы гидравлические системы удаления навоза из свиноводческих помещений. Применение данной системы удаления навоза имеет ряд преимуществ, но при его использовании, а также за счет сброса воды, используемой для технологических целей, получают жидкий навоз и навозные стоки влажностью до 98%, следовательно, возникают дополнительные трудности при переработке больших масс свиного навоза и стоков.
С точки зрения ресурсосбережения наиболее перспективным и рациональным является применение технологии анаэробной переработки жидкого свиного навоза. Применение данной технологии позволяет снизить затраты энергии на переработку и дальнейшую утилизацию без вреда для окружающей среды и получить в виде альтернативного энергоносителя — биогаза [1-3].
Целью исследования являются получение гидравлической модели для определения зависимости величины гидродинамических возмущений от параметров работы анаэробного фильтра и её экспериментальная проверка.
Объект исследования — процесс работы анаэробного фильтра с циклическим возмущением в фазе опорожнения.
Теоретические исследования
В фазе опорожнения работы анаэробного фильтра происходят циклические возмущения: изменение скорости движения субстрата в фильтре, обновление граничных поверхностей, лавинообразное выделение биогаза, разрушение пенного слоя (рис.) [4].
Для оценки возмущений в данной фазе принимается следующий параметр — скорость движения свободной поверхности стоков в фильтре которую можно определить, составив уравнение баланса объемного расхода. Допустим, что уровень жидкости в фильтре и баке в данный момент времени находятся, соответственно, на высоте и Х2. За бесконечно малый промежуток времени опорожнения &оп уровень в фильтре изменяется на а в сифоне — на течение жидкости можно считать установившимся. За это время из сифона вытечет объем жидкости Qсdtоп, в то же время в бак из фильтра притечет объем жидкости Qпрdtоп. В свою
очередь в фильтр поступает объем жидкости Qподdtоп, в то же время из фильтра вытекает объем жидкости Qвытdtоп. Тогда объем жидкости в баке сифонного отвода изменится на величину dW6, а изменение объема жидкости в фильтре составит Д а, -сок!
кг
0
I_с
I Рг _
т
— —■ — — ■
_ _- ■ ~—■
— — — .
— — — — .
. — — — —
— ■
- —
_ _ _ _ _
— — — — ■
Ы-
I д..
2-
-1-4
N
О
Рис. 1. Анаэробный фильтр
Баланс объемов данной системы представится в виде двух уравнений:
= Qdt - О dt ;
''б £-уС ОП £-УПП ОП'
(1)
Щ = Оеыт^оп - О„оЛп-В приведенной системе уравнений очевидным является то, что объем притекающий из фильтра в бак Qпрdt оп и объем вытекающей из фильтра в бак QвыTdtоп равны. Таким образом, можно приравнять уравнения системы, и изменение объема бака сифона будет выглядеть следующим образом:
dWб = ах122, (2)
где ю6 — площадь живого сечения бака сифона, м2.
Изменение объема фильтра:
Щ = ^2,, (3)
где Юф — площадь живого сечения фильтра,
м2.
Подставив уравнения (2) и (3) к приравненному уравнению (1) и решая, его относительно изменения высоты получим:
ОАоп - ОпоЛп -УбМ2-,
-. (4)
¿2,
Скорость движения свободной поверхности в фильтре dvф определим, поделив уравнение (4) на время опорожнения dtоп:
= ¿Ос - Оо -ЮбМУр
' (5)
Л
где ЛУб =—— — скорость движения свободам
оп
ной поверхности стоков в баке сифона, м/с.
Если рассматривать весь период опорожнения, то движение жидкости является неустановившимся, поэтому мы в уравнении (5) указали на изменение расхода сифонного отвода:
¿0, =асЛус, (6)
где юс — площадь живого сечения сифонного трубопровода, м2;
— скорость движения в сифонном трубопроводе, м/с.
Уравнение (5) в функциональном виде запишется как:
= Л ¿Ос; 0под; щ; ¿уб; Шф) =
ф
(7)
Ос =
1
4 + ^ +
(8)
^ +
-л 21с
¿ус
Ж
где Н = Z2 — Zc — напор сифона, м;
2= гс Л 1
& ду
— инерционная
-I
Л у
(9)
где Т — температура стоков, 0С;
W — влажность стоков, %;
К — химический состав стоков.
Из анализа уравнений (8) и (9) определим функциональную зависимость скорости движения свободной поверхности стоков в фильтре:
V = /
Дз;1ф;Лб;Лф;Т;Ш;К;1с;Лс; ^с;И;
¿ус Л!
(10)
=/(Л0С; д^ф^^йф),
где Дз — доза загрузки, %;
1ф — длина рабочей части фильтра, м; <б — диаметр бака сифона, м; <<ф — диаметр фильтра, м. Для определения расхода сифона воспользуемся уравнением Бернулли, уравнением неразрывности [5]. Принимаем сечение 1-1, совпадающее со свободной поверхностью жидкости в баке сифона, а сечение 2-2 — совпадающее с концом сифонного трубопровода (рис. 3), тогда
v1 = vб; = vc; Р1 = Ратм! Р2 = Ратм!
= ^2/ = ^с! а 1 = = ^о = а^ Окончательный вид уравнения расхода сифона примет вид:
скорость, где не учтены изменения при движении жидкости в баке сифона (ввиду их незначительности) и принято, что сечение сифонного трубопровода постоянно по длине и является цилиндрическим [5].
В функциональном виде формула (8) определяется как:
а = /\мс;Лс; ^с;уб;Н;ёУс
Для подтверждения полученных теоретических разработок были проведены экспериментальные исследования на демонстрационной установке.
Методика проведения экспериментальных исследований
Исследования проводились в два этапа:
1) проведение экспериментов на демонстрационной установке;
2) математическая обработка результатов экспериментов.
С целью изучения и наблюдения возмущающего воздействия, поиска конструктивных и технологических факторов, оказывающих влияние на величину гидродинамических возмущений, был проведен эксперимент на демонстрационной установке, которая является аналогом анаэробного фильтра с циклическим возмущением. Она позволяет смоделировать процесс работы фильтра, с имитацией выделения биогаза. В качестве субстрата применялась вода, а биогаза — атмосферный воздух. Планирование эксперимента разработано в соответствии с методикой Ю.П. Адлера [6].
В связи с поставленной целью были определены измеряемые параметры: — время фазы опорожнения, с; tзaп — время фазы заполнения, с; Ар — изменение давления на свободную поверхность, Па; АИсвпов. — изменение высоты свободной поверхности жидкости в фильтре, 10-3 м.
Оценочные параметры, определяющие величину и характер циклических возмущений, следующие: Ар(Y1) — изменение давления в газовом объеме фильтра, Па; V(Y2) — скорость движения свободной поверхности жидкости, м/с; tЦИKЛа(Yз) — продолжительность цикла возмущения в анаэробном фильтре, с; п (Y4) — частота протекания циклов, циклов/сут.
В результате проведения теоретических исследований [4], изучения априорной информации принято варьировать следующими факторами: Из (Х1) — высота погружения трубки в гидравлическом затворе, 10"3-м; АИ(Х2) — уровень колебания жидкости в баке
Ш
сифона, 10"3-м;
ф
Ш
(Х3) —
отношение пло-
б.с.
*Ш *
*
щади сечения фильтра к площади сечения бака сифона; Qж (Х4) — подача жидкости в анаэробный фильтр, 10-9 м3/с; Ьс (Х5) — напор сифона, 10-3 м.
Исследования проводились на двух уровнях — максимальном (+1) и минимальном (-1) с использованием дробного факторного эксперимента -25-2. Составлена матрица планирования эксперимента из восьми опытов, в натуральных и кодированных величинах, опыты рандомизированы во времени с помощью таблицы случайных чисел [6, 7].
В соответствии с составленным планированием и методикой проведения эксперимента была сооружена опытная установка, которая позволила провести визуальное наблюдение процесса опорожнения бака сифона с имитацией подачи субстрата и выделения биогаза
[7].
Результаты экспериментальных исследований
В результате выполнения эксперимента в соответствии с матрицей планирования получены следующие результаты:
- изменение давления в газовом объеме фильтр, Па:
Y1= -0,034 + (-0,087X0 + 0,829Х2 + + (-0,30Х3) + 0,011Х4 + (0,030Х5);
- скорость движения свободной поверхности жидкости в фильтре, м/с:
Y2 = -0,020 + (-0,005X0 + (-0,063Х2) + (-0,74X3) + (-0,042X4) + 0,026X5;
- продолжительность цикла возмущения в анаэробном фильтре, с:
Y3= -0,158 + (-0,091X0 + 0,094X2 + + 0,722X3 + (-0,11X4) + (-0,070X5);
- частота протекания циклов, циклов/сут.: Y1 = -0,279 + 0,00^, + (-0,155X0 +
+ (-0,656X0 + 0,165X4 + 0,023X5.
Из полученных регрессионных уравнений следует, что наибольшее влияние на параметры оказывают факторы X2 и X3, соответственно, уровень жидкости в баке сифона и отношение площади сечения фильтра и площади сечения бака сифона.
Проверка значимости коэффициентов выполнялась с помощью критерия Стьюдента, значимость коэффициентов подтвердилась.
Проверку адекватности полученных уравнений выполняли с помощью критерия Фишера, по результатам модели признаны адекватными.
Выводы
1. Разработанная гидравлическая модель работы анаэробного фильтра с циклическим возмущением подтвердилась экспериментальными исследованиями на демонстрационной установке.
2. В соответствии с гидравлической моделью скорость движения свободной поверхности стоков в фильтре зависит от физико-химических характеристик субстрата и техно-
логических параметров анаэробного сбраживания (Дз; Т; W; K), а также конструктивных особенностей установки (ф d6; d4>; lc; dc; £Çc; H). Наиболее технологичными факторами для управления скоростью являются Дз; d6; d4>; dc; Н. Причем увеличение dc; H будет способствовать повышению скорости, а увеличение Дз; d6; dф — снижать её.
3. Результаты теоретических и экспериментальных исследований могут быть использованы при проектировании конструкций анаэробных фильтров.
Библиографический список
1. Калюжный С.В., Пузанков А.Е., Варфоломеев С.Д. Биогаз: проблемы и решения. Биотехнология (Итоги науки и техники ВИНИТИ АН СССР). — М., 1988. — Т. 21. — 180 с.
2. Ильин С.Н. Ресурсосберегающая технология переработки свиного навоза с получением биогаза: автореф. дис. ... канд. тех. наук: 05.20.01. — Улан-Удэ, 2005. — 23 с.
3. Deublein D., Steinhauser A. Biogas from Waste and Renewable Resources. An Introduction. — WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim, 2008. — 450 p.
4. Евтеев В.К., Васильева А.С. Анаэробный фильтр с циклическим возмущением // Природопользование и аграрное производство: сб. ст. Междунар. науч.-практ. конф. (23-25 мая 2012 г.). — Иркутск: ИрГСХА, 2012. — С. 26.
5. Чугаев Р.Р. Гидравлика. — 5-е изд. — М.: ООО «БАСТЕТ», 2008. — 600 с.
6. Адлер Ю.П., Маркова Е.В., Грановский Ю.В. Планирование эксперимента при поиске оптимальных условий. — 2-е изд., пе-рераб. и доп. - М.: Наука, 1976. — 280 с.
7. Васильева А.С., Евтеев В.К. Результаты экспериментальных исследований на демонстрационной установке анаэробного фильтра // Экологическая безопасность и перспективы развития аграрного производства Евразии: матер. Междунар. научн.-практ. конф., по-свящ. 60-летию аспирантуры ИрГСХА (3-5 декабря 2013 г.). — Иркутск: ИрГСХА, 2013. — С. 12.
References
1. Kalyuzhnyi S.V., Puzankov A.E., Varfo-lomeev S.D. Biogaz: problemy i resheniya. Bio-tekhnologiya (Itogi nauki i tekhniki VINITI AN SSSR). — M., 1988. — T. 21. — 180 s.
2. Il'in S.N. Resursosberegayushchaya tekh-nologiya pererabotki svinogo navoza s poluche-niem biogaza: avtoref. dis. ... kand. tekh. nauk: 05.20.01. — Ulan-Ude, 2005. — 23 s.
3. Deublein D., Steinhauser A. Biogas from Waste and Renewable Resources. An Introduction. — WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim, 2008. — 450 p.
4. Evteev V.K., Vasil'eva A.S. Anaerob-nyi fil'tr s tsiklicheskim vozmushcheniem // Priro-dopol'zovanie i agrarnoe proizvodstvo: sb. st. Mezhdunar. nauch.-prakt. konf. (23-25 maya 2012 g.). - Irkutsk: IrGSKhA, 2012. - S. 26.
5. Chugaev R.R. Gidravlika. — 5-e izd. — M.: OOO «BASTET», 2008. — 600 s.
6. Adler Yu.P., Markova E.V., Granov-skii Yu.V. Planirovanie eksperimenta pri poiske op-timal'nykh uslovii. — 2-e izd., pe-rerab. i dop. -M.: Nauka, 1976. — 280 s.
7. Vasil'eva A.S., Evteev V.K. Rezul'taty eksperimental'nykh issledovanii na demon-stratsionnoi ustanovke anaerobnogo fil'tra // Ekologicheskaya bezopasnost' i perspektivy razvitiya agrarnogo pro-izvodstva Evrazii: mater. Mezhdunar. nauchn.-prakt. konf., pos-vyashch. 60-letiyu aspirantury IrGSKhA (3-5 dekabrya 2013 g.). — Irkutsk: IrGSKhA, 2013. — S. 12.
УДК 658.511:66.011
И.Ю. Шевченко I.Yu. Shevchenko
УПРАВЛЕНИЕ СЛОЖНЫМ НЕФТЕХИМИЧЕСКИМ ПРОМЫШЛЕННЫМ ПРОИЗВОДСТВОМ НА ОСНОВЕ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ И СОВРЕМЕННЫХ ИНФОРМАЦИОННЫХ ТЕХНОЛОГИЙ
THE CONTROL OF COMPLEX PETROCHEMICAL PRODUCTION BASED ON MATHEMATICAL MODELING AND MODERN INFORMATION TECHNOLOGIES
Ключевые слова: оптимальное управление, информационные технологии, экономико-математическое и компьютерное моделирование, информационно-моделирующая система, детерминированные модели, промышленное производство, пиролиз углеводородов, производственные задачи, оптимизация, модернизация.
Развитие современных производственных технологий связано с внедрением безотходных технологических циклов в составе производства как экологически и экономически выгодных. Моделирование и информационные технологии позволяют проводить вычислительные эксперименты с целью анализа возможных ситуаций при изменении параметров ведения процесса. Оптимальное управление промышленным производством происходит на базе современных информационно-моделирующих систем (ИМС), разработанных с помощью методов математического моделирования и искусственного интеллекта. Представлена компьютерная информационно-моделирующая система основного узла процесса пиролиза углеводородов, основанная на детерминированных математических моделях пиролизных трубчатых печей со змеевиками разной конструкции, в виде реакторов идеального вытеснения. Печи промышленного блока работают параллельно на различном углеводородном сырье. Одна из основных технологических проблем для печей пиролиза — это универсальность по отношению к сырью. Разработанная ИМС процесса пиролиза углеводородов базируется на математических моделях, чувствительных к изменению химического состава сырья, так как они учитывают механизм химических реакций, физико-химические закономерности протекающих явлений (кинетику процесса) и
покомпонентныи состав сырья и продуктов, а также технологические параметры процесса и конструктивные особенности печеИ. С помощью данноИ ИМС проведены серии прогноз-расчетов по целевым продуктам для пиролизного производства (завода «Этилен») Томского нефтехимического комбината на сырье Омского и других нефтеперерабатывающих заводов, также просчитывались варианты пиролиза широкой фракции лёгких углеводородов (ШФЛУ) как в качестве самостоятельного сырья, так и в качестве добавки к прямогонному бензину. Предлагаемая информационно-моделирующая система процесса пиролиза различного углеводородного сырья позволяет решать реальные производственные задачи по компоновке сырья, модернизации технологического оборудования и другие связанные с оптимальным управлением производства и уменьшением аварийных ситуаций.
Keywords: optimal control, information technology, economic-mathematical modeling and computer simulation, information-modeling system-mA, deterministic models, industrial production, pyroly-sis of hydrocarbons, production problems, optimization, modernization.
Modeling and information technologies allow to carry out the computational experiments to analyse possible situations when changing the parameters of process. Optimal control of industrial production is based on modern information-modeling systems (IMS), developed using the methods of mathematical modeling and artificial intelligence. This work presents computer information-modeling system of the pyrolysis of hydrocarbons primary node, based on deterministic mathematical models of pyrolysis tube furnaces with coils of different constructions, in