УДК 628.35-574.24.044
А.И. Федорова, Д.А. Крикун, С.Д. Угрюмова
Дальневосточный государственный технический рыбохозяйственный университет,
690087, г. Владивосток, ул. Луговая, 52б
ОСОБЕННОСТИ ФИЛЬТРОВАНИЯ МОРСКОЙ ВОДЫ И РАСЧЕТ НАСЫПНОЙ ФИЛЬТРУЮЩЕЙ ПЕРЕГОРОДКИ
Рассмотрен современный способ подготовки морской воды на рыбоводных и рыбоперерабатывающих предприятиях, представлена краткая характеристика насыпных фильтров с зернистой загрузкой. Определена зависимость эффективности процесса фильтрования от физико-химических свойств примесей в морской воде и насыпных зернистых материалов и от гидродинамических факторов. Описаны стадии процесса выделения взвешенных частиц из морской воды при фильтровании и определяющие факторы в доставке частиц к зернам загрузки. Представлены иллюстрации образования порового канала зернистого слоя, изменения сечения поры, зарастания поры в течение времени фильтрования и модели фильтровальной колонны. Описан эффект очистки морской воды как суммарный результат двух противоположных процессов - изъятия частиц из воды вследствие адгезии к зернам и отрыва ранее прилипших частиц под влиянием гидродинамического воздействия потока. Обобщена методика расчета насыпных фильтров с зернистыми перегородками.
Ключевые слова: зернистые материалы, процесс фильтрования, морская вода, насыпные фильтры, методика расчета.
A.I. Fedorovа, D.A. Krikun, S.D. Ugryumova FEATURES OF SEA WATER FILTRATION AND CALCULATION OF DRY BULK FILTERING PARTITION
Considered a modern preparation methods seawater aquaculture and fish processing plants, a brief description of bulk granular media filters. The dependence of the efficiency of the filtering process on the physico-chemical properties of the impurities in sea water and grained bulk materials, and of the hydrodynamic factors. Stage separation process described particulate from seawater by filtration and determinants in the delivery of particles to the load grains. Are illustrations of Education pore channel granular layer, changes in cross-section pore fouling the pores within the time filtering and model filter column. Describes the effect of the treatment of sea water as a result of the total of the two opposite processes - removal ofpar-ticles from the water due to the adhesion and detachment of grains previously adhered particles under the influence of hydrodynamic flow impact. The method of calculation of bulk filter with granular partitions.
Keywords: granular materials, the filtration process, sea water, bulk filters, the method of calculation.
Одним из наиболее современных способов подготовки морской воды на рыбоводных и рыбоперерабатывающих предприятиях является ее фильтрование через слои насыпных зернистых материалов с различными гранулометрическими, физико-химическими и сорбционными свойствами [5].
Насыпные фильтры с зернистой загрузкой имеют простую и надежную конструкцию, устойчивы к агрессивным условиям эксплуатации, эффективно очищают морскую воду от механических примесей.
При фильтровании морской воды через зернистые материалы примеси задерживаются в слое фильтрующего материала или на его поверхности, в результате чего происходит очистка морской воды. Эффективность процесса фильтрования зависит как от физико-химических свойств примесей в морской воде и насыпных зернистых материалов, так и от гидродинамических факторов. Пористая среда насыпных фильтров обычно формируется из зернистых материалов определенного диаметра (0,1-10-3 - 3-10-3) [3]. В сформированном зернистом слое частицы материала чередуются с пустотами, называемыми порами, которые образуют поро-
вые криволинеиные каналы и по которым протекает очищаемая морская вода (рис. 1). От формы, усредненных размеров и количества таких каналов в единице объема зернистого слоя зависит как гидравлический режим течения мор-скои воды, так и качество ее очистки.
Под пористостью слоя понимают долю свободного
объема в зернистом слое, м3/м3:
Уп
e=
пор
Ус
(1)
Усл -
Рис. 1. Образование порового канала зернистого слоя Fig. 1. The formation of pore channel granular layer
где Упор - объемы пор в слое зернистого материала; суммарный объем твердой фазы и пор в слое.
Процесс выделения взвешенных частиц из морской воды при фильтровании состоит из трех стадий:
1) перенос частиц из потока морской воды на поверхность фильтрующего материала;
2) закрепление частиц на зернах и в щелях между ними;
3) отрыв закрепленных частиц с переходом их обратно в поток морской воды [2].
Определяющими факторами в доставке частиц к зернам загрузки являются силы
инерции и диффузии. Закрепление частиц у поверхности фильтрующего материала происходит, во-первых, за счет действия межмолекулярных сил притяжения Ван-дер-Ваальса, определяющих прилипание (адгезию) взаимодействующих, разнородных по природе, зарядам и размерам фаз. Во-вторых, удержание частиц может осуществляться в щелях между зернами загрузки. Задержанные частицы взвеси при объемном фильтровании постепенно заполняют поры между зернами слоя (рис. 2), что приводит к уменьшению проходного сечения пор и увеличению гидравлического сопротивления слоя при постоянном расходе морской воды. Увеличение скорости движения морской воды в порах приводит к срыву уловленных частиц и перемещению их в глубину фильтрующего слоя. Таким образом, область интенсивного изменения концентрации взвеси в морской воде перемещается по высоте от верхней части слоя к нижней. Время, в течение которого достигается очистка морской воды до заданной степени, называется временем защитного действия загрузки. При его достижении либо по достижении времени предельной потери напора насыпной фильтр включают в режим очистки противотоком [2, 4].
На рис. 3 представлена иллюстрация зарастания поры в течение времени фильтрования.
Для математического описания закономерностей процессов удаления взвешенных веществ фильтрованием выделяют двумя сечениями «1-1» и «2-2» элементарный слой загрузки толщиной (Ах) на расстоянии (х) от ее поверхности (рис. 4). Площадь поверхности горизонтального сечения слоя примем равной единице. К сечению «1-1» подходит морская вода с массовой концентрацией частиц С1, а через сечение «2-2» вода выходит из слоя с концентрацией С2.
Эффект очистки морской воды как суммарный результат двух противоположных процессов - изъятия частиц из воды вследствие адгезии к зернам и отрыва ранее прилипших частиц под влиянием гидродинамического воздействия потока определяется по уравнению (2):
DC = DC1 - DC2
(2)
где БС1 - уменьшение концентрации частиц за счет их прилипания; ОС2 - увеличение концентрации за счет отрыва частиц.
Рис. 2. Изменение сечения поры:
1 - зерна загрузки; 2 - начальная конфигурация порового сечения; 3 - прилипающие частицы Fig. 2. Change the section of the pores:
1 - grain loading; 2 - the initial configuration of the pore section; 3 - adhering particles
Рис. 3. Иллюстрация зарастания поры в течение времени фильтрования: d„0 - условный диаметр поры в начальный момент времени т0; dn1 - условный диаметр поры в момент времени т1; d„2 - условный диаметр поры в момент времени т2; d„3 - условный диаметр поры в момент времени т3 Fig. 3. Illustration fouling the pores during the filtration time: dn0 - nominal diameter of the pores in the initial time t0; dn1 - nominal pore diameter at time t1; d„2 - nominal pore diameter at time t2; dn3 - nominal diameter of the pores in the time t3
Рис. 4. Модель фильтровальной колонны Fig. 4. Model of filter columns
Уменьшение концентрации частиц за счет их прилипания (^С1) определяется по уравнению (3):
щ .с .дх, (3)
Т1
где Сср - средняя концентрация частиц в объеме выделенного слоя ; Ах - толщина слоя; п1 - количество частиц, прилипших к зернам фильтрующей загрузки; тг- - масса одной прилипшей частицы; т1 - время, в течение которого происходит прилипание частиц к зернам фильтрующей загрузки.
Увеличение концентрации за счет отрыва частиц определяется по уравнению (4):
ПС2 = П2 • г Дх / н, (4)
т2
где г - количество осадка, накопившегося к данному моменту времени в единице объема элементарного слоя загрузки (плотность насыщения загрузки осадком); г • Ах - количество накопившегося к данному моменту времени осадка; н - скорость фильтрования; п2 - количество частиц, оторвавшихся от зерен фильтрующей загрузки; т2 - время, в течение которого происходит отрыв частиц от зерен фильтрующей загрузки.
Подставляя значения ПС, ПС1 и ПС2 в уравнение (2), получают
дС п1 • т( п2/т2 (5)
— =——- • Сср —2—2 • г .Дх. (5)
дх т1 ср н
Уравнение (5), являясь основным уравнением процесса фильтрования, содержит две переменные (средняя концентрация частиц в объеме выделенного слоя - Сср и количество осадка, накопившегося к данному моменту времени в единице объема элементарного слоя загрузки - г), поэтому его недостаточно для описания процесса.
Расчет насыпных фильтров в основном сводится к определению расхода жидкости и величины гидравлических потерь на единицу поверхности фильтровального материала. Ввиду сложности пористой структуры большинства фильтрующих материалов, состоящей из соединенных между собой пор и сложной сети каналов (капилляров) различной формы и размеров, которые к тому же в ряде случаев изменяются под действием перепада давления, установить для большинства фильтрующих материалов закономерность и дать аналитическое выражение для характеристики потока жидкости практически невозможно. Поэтому гидравлические характеристики фильтрующего материала определяются, за исключением отдельных случаев, экспериментальным путем. Расход насыпного фильтра может быть выражен в общем случае зависимостью, вытекающей из закона Пуазейля [1]:
2 = д • ^ = к(Р .Др / II), (6)
где д — расход (пропускная способность) в единицу времени через единицу площади фильтровального материала, м3/м2-ч; /л — коэффициент динамической вязкости фильтруемой жидкости Н-с/м2; Р — площадь поверхности фильтрующего элемента, м2; Др — перепад давления на фильтре; к - коэффициент пропорциональности, Н/м2.
Расход д определяется по уравнению (7):
д = к •Др / ц. (7)
Коэффициент к, представляющий собой удельную пропускную способность единицы
площади поверхности фильтрующего материала при перепаде давления 1 кг/см2 и вязкости жидкости 1 л/см2, определяется по уравнению (8):
к = д. ц/Др. (8)
Опыт показывает, что коэффициент к для поверхностных фильтров сохраняется при всех прочих равных условиях практически постоянным в широком диапазоне расходов и перепадов давления, благодаря чему представляется возможным использовать его в качестве сравнительной характеристики гидравлического сопротивления фильтрующего материала.
Процесс фильтрования зависит от многих технологических параметров и в первую очередь от свойств зернистого слоя, свойств фильтрационной среды и примесей, от гидродинамического режима фильтрования [1, 5].
Важнейшими характеристиками пористой среды являются порозность и удельная поверхность, которые определяются по уравнениям (9, 10).
v - V
0 к ’
где £0 - порозность слоя; V- общий объем зернистого слоя, м3; V - объем частиц, м3.
(9)
а = 6(1 £°) , (10)
d э
э
где а - удельная поверхность слоя, м2/м3; ёэ - эквивалентный диаметр частицы, м; а - коэффициент формы.
Частицы задерживаются поверхностью зерен под действием молекулярных сил, электростатических сил, сил химического сродства и адсорбции. Величина сил прилипания зависит от крупности и формы зерен, скорости потока, температуры воды и свойств примесей [5].
Кинетика фильтрования и материальный баланс определяются по уравнениям (11, 12).
дС ь С (11)
-----= Ь. С - а. р, V '
дх
где С - концентрация примесей в промывной воде; х - толщина слоя загрузки; Ь, а - константы скорости прилипания и отрыва частиц; р - плотность насыщения фильтрующего слоя задержанным осадком.
др=-н дС (12)
дт дх,
где н - скорость фильтрования.
При решении уравнений (11) и (12) получается общее уравнение процесса фильтрования:
д2С дС , дС Л (13)
ан----+ Ь---------------------------------------= 0. (13)
дхдт дх дт
Уравнение (13) имеет решение в виде бесконечного ряда, и его трудно использовать в расчетах.
В процессе фильтрования происходит накопление загрязнений в слое загрузки. В какой-то момент наблюдается вынос частиц в фильтрат с ухудшением его качества. Продолжительность работы фильтра до проскока частиц в фильтрат называют временем защитного действия загрузки тз. По мере загрязнения фильтрующего слоя уменьшается его пороз-ность и увеличивается сопротивление при прохождении через него загрязненной морской воды, т.е. растет потерянный напор. Время работы фильтра до достижения потерянного напора предельной величины Нп обозначают через тн [1].
Оптимальным условием работы насыпного фильтра является тз ~ тн. Время защитного действия загрузки тз и время работы фильтра до достижения потерянного напора тн определяются по уравнениям (14,15).
ь = 1--( * - Ь1, (14)
ка ^ Ь
где к, х0 - константы, зависящие от эффекта очистки (справочные данные).
Нп - Н0 Ь (15)
тн =—п----^ --х, (15)
Нп Р,)
где Н0 - потеря напора в чистой загрузке; Г(А) - параметр, зависящий от величины предельной насыщенности порового пространства отложениями А.
Сопротивление фильтрующего слоя в любой момент времени определяется по уравнению (16):
X X
Н = [гёх = /0 [(—-—)3 ёх, (16)
0 0 -о -А-
где г0 - сопротивление единицы толщины фильтрующего слоя при прохождении через него чистой жидкости; г - сопротивление единицы фильтрующего слоя с задержанными частицами в любой промежуточный момент времени; Ае - удельный объем осадка, накопившегося в фильтрующем слое.
Сопротивление единицы толщины фильтрующего слоя при прохождении через него чистой жидкости г0 и сопротивление единицы фильтрующего слоя с задержанными частицами в любой промежуточный момент времени г определяются по уравнениям (17, 18).
. 0,188^(1 --0)2 (17)
г° = 3 , (17)
а э -0
где у - коэффициент, характеризующий место расположения выпуска сточных вод (для берегового выпуска у = 1, для выпуска в сечении русла у = 1,5); /л — коэффициент динамической вязкости фильтруемой жидкости Н-с/м2; е0 - порозность слоя; ёэ - эквивалент-
ный диаметр частицы, м.
г = г0
( — V —0
—0 — А—
(18)
V
где г0 - сопротивление единицы толщины фильтрующего слоя при прохождении через него чистой жидкости; г - сопротивление единицы фильтрующего слоя с задержанными частицами в любой промежуточный момент времени; Ае - удельный объем осадка, накопившегося в фильтрующем слое.
Площадь насыпных фильтров Г определяется по уравнению (19), м2:
р =__________________2________________ (19)
т ■ и — 3,6 • Пп ■ и-тп — Пп ■ ^ р-тф
а
где Q - расход морской воды, м3/ч; то - продолжительность работы фильтра, ч; wp - расчетная скорость фильтрации, м/ч; пп - число промывок фильтра в сутки; и - интенсивность промывки, л/(м2-с); тп - продолжительность промывки, ч; тпф - продолжительность простоя фильтра в связи с промывкой, ч (тпр = 0,33 ч).
Составлено уравнение материального баланса, произведен расчет объема осадка Voc и очищенной морской воды (фильтрата) Vф при производительности насыпного фильтра Goc = 1200 кг/ч, влажности осадка woc = 40 % (масс.), начальной концентрации морской воды перед фильтрованием по твердой фазе хс = 5 % (масс.), конечной концентрации морской воды после фильтрования Хф = 0; плотности морской воды рж = 1200 кгм3 и плотности осадка рос = 1350 кгм3 1 [1].
Уравнение материального баланса по потокам фильтра: Gc = Goc +G ф,
уравнение материального баланса по твердой фазе: Gcxc = Goc xo c +G фхф,
где Хос = 100 - W ос.
Объем осадка равен: Voc = Goc /poc = 1200/1350 = 0,89 м3/ч.
Для определения объема фильтрата (очищенной морской воды) по аналогичной формуле выразим массовую производительность по фильтрату через производительность по осадку из уравнений материального баланса. Для этого сначала выразим количество фильтрата через количество морской воды, не прошедшей очистку, и осадка Gф = Gc +G ос, а количество морской воды, не прошедшей очистку через количество осадка: Gc = (GocХоc + GcХф ~GосХф )/ Хс = Goc (Xoc - Хф ) /(xc - Хф X откуда с учетом условия ХФ = 0
получим:
Gф = Gc - Goc = Goc (Xoc / Xc -1) = Goc ((100 - Woc)/ Xc -1) = 1200((100 - 40)5 -1) = 18000 кг/ч. Объем фильтрата (очищенной морской воды) равен:
Vф = Gi> /рж = 18000/1200 = 15 м3/ч.
Рабочая площадь плоских сеток определяется по уравнению (20):
Fc = ^ -*2. (20)
W
гуе
где Q - расход морской воды, м3/ч; кн - коэффициент неравномерности; wc - скорость движения морской воды в сетке (для плоских сеток wc = 0,2-0,4 м/с); к2 - коэффициент загрязнения сетки (к2 = 1,2-1,8).
При этом коэффициент к1 определяется по уравнению (21):
к, = (^d^a + F1), (21)
dn
где йя - размер ячеек в свету, равный 0,5-5 мм; dn - диаметр проволоки сетки, равный 0,3-2 мм; F1 - часть площади, занимаемой рамами и шарнирами.
Площадь фильтрующей поверхности насыпного фильтра равна
F = к1 -Q-кн (22)
ф 1 ’ к2 ^в^ф
где т0 - время работы насыпного фильтра в течение суток, ч; ~№ф - скорость фильтрования, принимают в пределах 4-20 м/ч; к^ к2 - коэффициенты: к1 = 1,03, к2 = 0,63.
Потери напора на чистой сетке определяются по уравнению (23):
А =* ^ (23)
2- Л^
где £0 - коэффициент сопротивления (для ламинарного режима при Яес < 4 равен £0 = 6,7/
1/3
Яес и для турбулентного режима при Яес > 4 равен ^0 = 2,38/ Яес ); н,1 - скорость движения морской воды на подходе к сетке, м/с; А1 - свободное сечение сетки, доли единицы.
Потери напора при эксплуатации насыпных фильтров с разделительной сеткой Ат в момент времени т определяются по уравнению (24):
А =с ______£1^_______ (24)
' ' 2-Л12(1 -рк)2’
где £т - коэффициент сопротивления в момент времени т (справочная величина); вк - коэффициент загрязнения сетки (справочная величина); н,1 - скорость движения морской воды на подходе к сетке, м/с; Яес = w1Rc / V - число Рейнольдса; Кс - гидравлический радиус сетки; Кс = Л1 / 2пN; N - число проволок на единицу длины; V - коэффициент кинематической вязкости, м2/с.
Представленная нами методика позволяет осуществлять расчет насыпных фильтров с зернистыми перегородками. Данная методика будет использована в учебно-методическом пособии по расчету однослойных и многослойных фильтров.
Список литературы
1. Ветошкин А.Г. Процессы и аппараты гидросферы. - Пенза: Изд-во Пензенского гос-университета, 2004. - 188 с.
2. Копылов А.С. Водоподготовка в энергетике / А.С. Копылов, В.М. Лавыгин, В.Ф. Очков. - М.: МЭИ, 2006. - 320 с.
3. Угрюмова С. Д. Оценка эффективности процесса фильтрования морской воды через зернистые перегородки / С.Д. Угрюмова, А.И. Федорова // Науч. тр. Дальрыбвтуза. - Владивосток: Дальрыбвтуз, 2013. - Т. 28. - С. 137-144.
4. Угрюмова С. Д. Экспериментальное исследование насыпной загрузки в процессе фильтрования морской воды / С. Д. Угрюмова, А. И. Федорова, И. В. Панюкова // Приморские зори-2012. Междунар. науч. чтения: сб. науч. тр. - Владивосток, 2012. - С. 241-245.
5. Угрюмова С. Д. Физические особенности процесса фильтрования морской и пресной воды через насыпные перегородки / С.Д. Угрюмова, А.И. Федорова, С.А. Лоншаков // Фундаментальные и прикладные вопросы естествознания: материалы 54 Всерос. науч. конф. -Владивосток: ТОВМИ, 2011. - Т. III. - С. 191-193.
Сведения об авторах: Федорова Александра Игоревна, аспирант; Крикун Дмитрий Александрович, аспирант;
Угрюмова Светлана Дмитриевна, доктор технических наук, профессор.