Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Белоглазов Фильтрование технологических пульп

.pdf
Скачиваний:
101
Добавлен:
16.12.2019
Размер:
7.07 Mб
Скачать

осадка (г, Лос) и фильтра (RJ. Так же, как и уравнение (1.138), его

интегрируют для двух случаев.

Для фильтрования с постоянным перепадом давлений после

интегрирования в пределах от V = V, и т =

до V =

V2 и т = т2

получится зависимость

 

 

 

 

2

и

щ жи

 

(1.144)

 

 

 

ИЛИ

 

 

 

 

 

То

т.

(V2 -

Vi) + гРжц у

+ £ p A l

(1.145)

V2 -

Vj

ApF

ApF2 1

ApF

 

В реальном процессе достичь постоянства давле­ ния на начальном этапе фильтрования практически невозможно, следовательно, на основании уравнения (1.145) можно сделать вывод о нарушении линейности

 

зависимости между

V2 -V 1

 

 

 

 

 

 

 

И

V2

-

( в

этот

период

 

(рис. 1.11). Поэтому уравне­

Рис. 1.11. Кривая фильтрования несжи­

ние

(1.145) применимо для

маемого осадка на несжимаемой филь­

описания

процесса

филь­

трующей перегородке при Др = const

трования

при

постоянном

 

перепаде

давлений

только

 

по

истечении

некоторого

времени т, от начала фильтрования, необходимого для стабилизации

давления в агрегате.

Фильтрованию с постоянной скоростью подобный эксцесс не свойственен, и в качестве пределов интегрирования уравнения

(1.143) можно взять т = О, V =

0 и т =

т,, V = V,. После инте­

грирования получим

 

 

 

ApF

 

г»жи

V1+-F1

(1.146)

 

и

У

или

 

 

60

 

 

т,

 

гр и

гржл„

 

V

=

ApF

ApF

(U47)

V,

 

 

Выражение (1.146) можно привести к виду

и

Ц1Ж11

(1.148)

 

аналогичному известному уравнению Рутса [49]

 

V2 +

2CV = Кх,

(1.149)

где К — константа, комплексно учитывающая физико-химические свойства осадка и фильтрата, м3/с; С — объем фильтрата, который, проходя через1 м2 фильтрующей поверхности,образует осадок, со­ противление которого равно сопротивлению фильтрующей пере­ городки, м3.

Отличие между ними заключается только в отсутствии коэф­ фициента 2 перед V, в уравнении (1.148).

Рутсом тоже был предложен раздельный учет сопротивлений

осадка и перегородки, причем

 

 

 

^

= KV,

 

(1.150)

где

 

 

 

Rn = КС.

 

(1.151)

Используя введенные обозначения, для режима фильтрования

при Ар = const уравнение (1.141) можно представить в виде

 

dr

Л-Р— .

 

(1.152)

K(V + С)

 

 

Обозначив Ар/К = к', получим общее уравнение скорости

фильтрования:

 

 

 

(V +

Q d V = k'dt.

 

(1.153)

Проинтегрируем выражение (1.153) в пределах от 0 до (V + С)

и (т + т0):

 

 

 

V +С

 

т+то

(1.154)

J (V + C)d(V + С) = к'

Jdx

о

 

о

 

и в результате получим параболический закон фильтрования

61

V + C)2 =i'(T + T0),

(1.155)

где т0 — продолжительность образования слоя осадка с сопро­

тивлением, равным сопротивлению фильтрующей перегородки, с. Обозначив 2к' через К, перепишем уравнение (1.155) в виде

+ С)2 = К(т + т0).

(1.156)

Графическое изображение полученной зависимости в коор­ динатах У—т и (У+С)— (т + т0) представлено на рис. 1.12, откуда

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

 

Для

начального

момен­

 

та, когда

х

=

0 и

 

У =

О,

 

уравнение

(1.156)

 

прини­

 

мает вид

 

 

 

 

 

 

 

 

 

С2

=

Кх.

 

(1.157)

 

Выведенное

 

Рутсом

 

уравнение

(1.149)

 

удобно

 

использовать

для определе­

 

ния констант фильтрования

 

С и К, которое может вы­

 

полняться как

с примене­

 

нием

математического

ме­

Рис. 1.12. Кривая фильтрования в коор­

тода, так и графически. В

любом

случае для

решения

динатах У, х и (У + С), (т+ т0)

задачи

потребуются

опыт­

 

 

ные

значения

объемов

 

фильтрата У( и У2 и данные

о времени, за которое они собраны,

и т2 [68].

 

 

 

 

 

 

Математический подход реализуется путем решения системы

уравнений, составленной на основе (1.149):

 

 

 

 

 

 

 

V* + 2CVj = Кхг

 

 

 

 

 

 

 

2

 

 

 

 

 

(1.158)

V22 + 2СУ2 = Кх2

 

 

 

 

 

 

 

Для простоты графического решения предварительно продиф­

ференцируем уравнение (1.149) по переменным У и т:

 

 

 

 

2(У + Q d V =

Kdx

 

 

 

 

 

(1.159)

или

 

 

 

 

 

 

 

 

av

_

л

(1.160)

dx

 

2(V + С )'

 

Уравнение (1.160) удобно представить в виде

 

dx

 

2

(1.161)

d V

= — (V + C).

 

К

 

Последнее уравнение и является основой графического определе­

ния констант К и

С.Для

этого на оси ординатоткладывают вели­

чину, обратную скорости

фильтрования,dx/dV *Ат/АУ, а на оси

абсцисс — объемы собранного с единицы поверхности перегородки фильтрата V (рис. 1.13). Экстраполируя построенную по опытным точкам прямую до пересечения с абсциссой, находят величину константы С, а тангенс угла наклона прямой отвечает К.

В некоторых случаях априори известно удельное сопротивление осадка и перегородки г. Так, сопротивление перегородки может быть указано фирмой-производителем, а сопротивление осадка при проведении оценочных расчетов может быть принято в соответст­ вии сосправочными данными. Если все эти сведения имеются, константы фильтрования могут быть определены присовместном рассмотрении уравнений (1.148) и (1.149).

Константа фильтрования К, отнесенная к 1 м2 фильтрующей поверхности, связана с удельным сопротивлением осадка зависи­

мостью

 

 

Куд =

2Др/(цжгц),

(1.162)

ИЛИ

 

 

к „ = 2Др(1

- еам)/(цжгамрж).

(1.163)

Рис. 1.13. Экспериментальная зависимость Ах/АУ = /(V) для определения констант фильтрования К и С

62

63

Если константа Куд известна, то два последних выражения могут

также использоваться для определения удельного сопротивления слоя.

Аналогично

С уд

= Гп/ Ш <

(1Л64)

ИЛИ

 

 

С уд = ^

“ ^ л У К Р ж ) '

(1Л65)

где Суд = C/F — константа С, отнесенная к 1 м2 фильтрующей поверхности, м3/м 2; гп — удельное сопротивление перегородки,

Н'С/м4.

Выражения (1.164) и (1.165) можно использовать и для опре­ деления ти гп.

Забегая вперед, отметим, что в случае процесса фильтрования с образованием сжимаемого осадка при вычислении констант К и С также пользуются подобными зависимостями [40, 62].

При нарушении структуры пор во время фильтрации удельное сопротивление осадка г растет. Можно предположить, что

г = г0Дрй,

(1.166)

где г0 — удельное сопротивление осадка при Ар =

1 ат; b

показатель сжимаемости осадка, 0 < Ъ < 1. Для осадков, которые не изменяют свою структуру во время фильтрации, Ъ ® 0, а для осадков, изменяющих свою структуру (сжимаемых), Ъ * 1. Так, для осадка гидроксида алюминия b « 0,95.

Следовательно, константа Куд для сжимаемых осадков будет

иметь вид

 

 

 

 

 

 

К

2Др

2Др'~ь

(1.167)

 

= ---------^ —

= —

,

 

У?

ИжГ0ДР °м

РжГ0ам

 

или

 

 

 

 

 

К

 

2Др(1-е)

2Др1-Ь(1 - s)

(1.168)

= ------ \

=

-------■

УЯ

РжГ0АР Ржам

РжГ0Ржам

 

Константа С в меньшей степени зависит от перепада давления при фильтрации, и для инженерных расчетов с достаточной долей точности можно пользоваться уравнением

С

(1.1.69)

79

Г0Ржа м

64

Уравнения (1.167) и (1.168) наиболее часто используются для определения сопротивления ткани и осадка при известных констан­ тах фильтрования.

Необходимая поверхность фильтрования F определяется как отношение общей производительности передела по фильтрату О к удельной производительности V, рассчитываемой по основному уравнению фильтрования

V = у/с2 + К х - С,

(1.170)

F = Q/V.

 

Число стандартных фильтров i можно определить,

исходя из

поверхности фильтрации одного фильтра F' и коэффициента запаса г|з, который для фильтров различных типов и в зависимости от

фильтруемых сред колеблется в широких пределах [85]:

«# = £ Лз-

(1-171)

1.8. Фильтрование с образованием сжимаемого осадка

Уравнения кинетики процесса фильтрации с образованием не­ сжимаемого осадка (1.138) и (1.143) непригодны в случаях, когда сжимаемостью осадка пренебречь невозможно, что объясняется непостоянством давления жидкости в ходе процесса.

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

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

шением толщины слоя осадка снижаются сжимающие его усилия (рис. 1.14). В произвольном сечении слоя осадка (показано на рисунке пунктиром) толщиной 62 перепад давления Др2 больше, чем

перепад давления Др, в более тонком слое толщиной 5,, при условии,

65

что 8j лежит внутри 82. Отсюда следует, что с уменьшением тол­ щины слоя сжимаемого осадка его пористость при прочих равных условиях возрастает [63].

 

Изменение

давления,

а сле­

 

довательно,

пористости и

удель­

 

ного сопротивления по толщине

 

слоя

осадка

обусловливается

 

свойствами последнего, экспери­

 

ментальное определение

кото­

 

рых по высоте слоя представляет

 

большие

трудности.

Поэтому

 

способ

расчета

процессов

филь­

 

трования с образованием сжима­

 

емого

осадка

требует введения

 

ряда допущений:

 

 

 

 

1) давление на границе раз­

 

дела

осадка

и

фильтрующей

 

перегородки р, мало по сравне­

 

нию с давлением

фильтрования

Рис. 1.14. Изменение давления р по

р, т. е. удельное

сопротивление

толщине 8 слоя осадка:

фильтрования

г

+

гп

можно

1 — несжимаемый осадок; 2 —

считать постоянным;

изменяется

сжимаемый осадок

2)

давление

не

(не подвергается возмущениям) в ходе процесса;

3) влажность осадка неизменна на протяжении всего процесса. Эти допущения обычно неправомерны для начальной стадии

процесса, когда осадок имеет малую толщину [28].

Случай фильтрования с образованием сжимаемого осадка мо­ жет быть рассмотрен с позиции общей теории фильтрования. Для этого из уравнения (1.126) необходимо исключить высоту осадка Лос

на основании материального баланса (высота осадка может быть выражена через объем фильтрата, необходимый для образования

данного количества осадка):

(1.172)

- £)Pm = 9'(V + eh F ),

где д' — масса твердой фазы в исходной суспензии, приходящаяся на единицу объема жидкости в этой суспензии.

Тогда уравнение (1.126) можно записать в виде

 

d V

Ар

в

з

(1.173)

 

 

 

,2

 

 

 

 

Y9

FPr

и для сжимаемого осадка, т. е. при условии зависимости в от Ар, получим (при к'к =5)

66

(1.174)

Fdt 5g'gV i (1- 8)4 ,

Последнее уравнение можно записать в виде

d V

F

P' J 2

dp

F d x ~

g 'p V

I

(1.175)

Roc’

где R’oc — удельное сопротивление осадка, м/кг.

Соотношение между Л' , е и Ар определяется экспериментально |22|. Для этого проба исследуемой суспензии помещается в цилиндр г пористым дном, на котором осадок задерживается, а фильтрат вытекает. В цилиндр устанавливают поршень и медленно сдавли­ вают осадок. Нагрузка на поршень может изменяться в заданных пределах. При каждой нагрузке определяется пористость осадка (по отметке положения поршня). Фильтрат вновь подается на слой осадка. При этом для каждой нагрузки по уравнению (1.175) рассчи­ тывается R'ot как постоянная величина. Зная г и R'ОС, можно рассчитать S . На рис. 1.15 приведены типичные для фильтрования

с образованием сжимаемого осадка зависимости s, Syg и R'oc от Ар,

полученные описанным выше методом испытания проницаемости. Результаты таких опытов могут быть использованы и для расчета процесса фильтрования в промышленных условиях.

Уравнение скорости фильтрования по Дарси можно расширить, чтобы учесть движение твердых частиц в сжимаемом осадке [97]:

4 = ™ т - — '9р,

(1.176)

где д — вектор локальной скорости жидкости в осадке, м/с; w

вектор локальной скорости твердых частиц, м/с; е — пористость осадка; к — коэффициент проницаемости (по уравнению Дарси), м2;

р - динамический коэффициент вязкости фильтрата, кг/(мщ); Vp — оператор давления, под которым проводится фильтрование.

Первый член правой части уравнения (1.176) соответствует чле­ ну, учитывающему конвективный поток в первом законе Фика для диффузии.

Для процесса фильтрования, в котором q и Vp солинейны, в уравнение (1.176) можно ввести вектор мгновенной локальной ско­ рости твердых частиц г = (1 — £)vrm:

67

в

Рис. 1.15. Характеристики сжимаемого осадка а — зависимость пористости осадка от удельной поверхности частиц; б —

зависимость пористости осадка от Ар; в — зависимость удельного сопро­

тивления осадка от Ар;

9,8; 2 — СаСОэ, pH =

10,3; 3 — ТЮ2

1 — СаСОэ (флокулированный), pH =

(флокулированный), pH = 7,8; 4 — ТЮ2, pH = 3,5; 5 — ZnS, pH = 9,1

Я = е г

к

(1 177)

Vp,

ц

где е = е/(1 — б) — локальное отношение свободных объемов.

Из уравнения (1.177) следует зависимость, описывающая про­ цесс фильтрования с образованием сжимаемого осадка [58]:

^

=

------ ------- ^

= -nR x(gx - ехгх),

(1.178)

dm x

 

р_(1-Е.)

dx

 

где рх — давление, сжимающее частицы осадка на расстоянии х от фильтрующей перегородки, Па; тх — масса твердых частиц осадка,

приходящаяся на единицу площади на расстоянии х от фильтру­ ющей перегородки, кг/м2; рт п л о т н о с т ь твердых частиц в осадке,

кг/м3; ц — динамический коэффициент вязкости фильтрата, кг/(м.с); Rx — локальное сопротивление осадка, м/кг; qx — локальная

скорость фильтрования на расстоянии х от фильтрующей пере­ городки, м/с.

В результате интегрирования уравнения (1.178) по толщине слоя осадка получены следующие зависимости:

д 1= —

= ----- - ----- ,

(1.179)

dx

2{V + и0)

 

где д, — скорость фильтрования на выходе из осадка, м/с; V — объем фильтрата, отнесенный к единице площади фильтра, м3/м 2;

^ф.п(1 ~~ ас)

м

10Л1

и0 = — ------------ ;

(1.180)

RcpPc

 

 

К =

(1.181)

l^cp Рс

 

 

где Яфп — сопротивление фильтрующей перегородки,

м-1; R

среднее удельное сопротивление фильтрования, м/кг; с — массовая доля твердой фазы в суспензии, дол. ед.; а — отношение массы млажного осадка к массе сухого; р — плотность фильтрата, кг/м3; р — давление фильтрования, Па.

При фильтровании с постоянным давлением допустимо считать R неизменным, так как при образовании слоя осадка значительной

толщины давление р, на границе раздела осадок—фильтрующая

среда мало по сравнению с давлением на поверхности осадка.

.Значения а также становятся постоянными вскоре после начала процесса. Следовательно, и значения К и и0 приобретают смысл

констант фильтрования [94].

В начальной стадии фильтрования, когда толщина слоя осадка незначительна, сопротивление Rcp = f(p — р,) и непрерывно изме­

няется с увеличением количества фильтрата V. Изменения давления и величины а также вызывают изменения К и о0 в начальной стадии

процесса.

Чтобы проинтегрировать уравнение (1.179), необходимо ввести поправки на изменение констант К и и0 в начале процесса.

Перепишем уравнение (1.179) в виде

68

69

dx

-(У + и0) + Ф(У).

(1.182)

d V

К

 

Здесь К и » , имеют постоянные значения, достигаемые после окон­

чания начальной стадии фильтрования, которой соответствует зна­ чение объема фильтрата У ач; ср(V) — функция объема фильтрата,

учитывающая поправку к d x /d V в начальной

стадии процесса. При

этом необходимо учитывать, что при 0 < V <

V функция cp(V) # О,

а при V > VHm функция cp(V) =

0.

 

Уравнение (1.182) можно проинтегрировать в пределах от V = 0

для т = 0 до V =

V

для х =

х:

 

 

2

^+ О 0

V нач

 

х =

-

J

(V + u 0)d(V + v 0)+ j<p(V)d{V)

(1.183)

 

к

 

 

После интегрирования получим продолжительность фильтро­

вания:

 

 

 

 

 

 

 

 

— (V2 + 2Vu0) + х0,

(1.184)

 

 

 

К

 

 

где

 

 

 

 

 

 

 

 

 

V нач

 

 

 

 

х0 =

Jcp(V4dV.

(1.185)

 

 

 

 

о

 

При V = 0 в уравнении (1.184) получим х = х0, которое назы­ вают временем начала выделения фильтрата.

Недостатком этого метода расчета является трудность вычис­ ления dx/dV с необходимой точностью, однако для обработки опыт­ ных данных, собранных в основной стадии процесса, можно исполь­ зовать уравнение (1.184).

Проанализируем функцию ф(V):

Ф(^) = | —

dx Л

(1.186)

d V j OCH'

 

 

,d V y

 

или в соответствии с уравнением

(1.182)

 

 

 

1

1

{

нач

(1.187)

cp(v) = 2V

+ 2

 

 

 

 

 

\ к

нач

 

Кк

нач

 

V

 

\

 

70

 

 

 

 

 

После подстановки значений К и и0 из уравнений (1.180) и (1.181) получим

Урр

R

R

Ф (У ) =

1

 

 

 

 

а

И

 

V с

 

 

При допущении, что сопротивление фильтрующей перегородки не

тхнач 7-«ОСН

л

 

 

 

 

меняется и н ф,п

~ R<p.n

= 0, уравнение (1.188) упрощается:

 

Урр R

 

 

1

 

 

 

 

а

ф(У):

Vc

®О СН

| R OCH Vc

'

 

 

1

 

(1.189)

 

 

—- а

 

 

 

 

VC

®ОСН

 

 

 

Vc

'

/

 

Так как аиач >

аосн и

п нач

п осн

то

можно

сделать вывод, что

Rlp.n s кф п ,

член уравнения в квадратных скобках, а следовательно, функция ф(V) и х0 могут быть отрицательными.

Экспериментально установлено, что удельное сопротивление сжимаемого осадка Roc может быть выражено эмпирической фор­

мулой

 

R„ = Ар",

(1.190)

где А и п — постоянные для данного материала величины; р — давление фильтрования, Па.

В результате исследований, проведенных в НИИхиммаше, вы­ явлено существенное влияние концентрации твердой фазы в сус­ пензии на основные параметры процесса фильтрования. Установ­ лено, что удельное сопротивление осадка уменьшается с увели­ чением содержания твердой фазы. Это необходимо учитывать при выборе режима работы фильтра периодического действия или условий разделения на фильтре непрерывного действия [101].

Характер влияния концентрации на кинетику фильтрования несколько отличается от приведенных выше данных, но исследо­ вания, углубляющие наши представления о механизме процесса и влиянии отдельных параметров на скорость фильтрования, про­

должаются.

Уравнение кинетики фильтрования может быть представлено в виде

71

d V

А р

,

(1.191)

=

x— -------------

dx

ц(Ас q mBV + Кфп)

 

 

где V — объем фильтрата с единицы фильтрующей поверхности, м3/м 2; Ар — перепад давлений, Па; ц — динамический коэффициент вязкости фильтрата, Па-с; qna — слой твердой фазы, отлагающийся

при получении единицы объема фильтрата, м /м3; с — объемная доля твердой фазы в суспензии; А — коэффициент пропорциональности; к — показатель степени; Яфп — сопротивление фильтрующей пере­

городки, (Н-с)/м3.

1.9. Фильтрование с закупориванием пор фильтрующей перегородки

При разделении вязких суспензий, содержащих небольшие количества мелкодисперсных частиц твердой фазы, происходит проникновение твердых частиц в пористый фильтрующий слой и отложение их на стенках пор. В природных условиях такой процесс наблюдается очень часто (проницание грунтовых вод в геологи­ ческие массивы). Кроме того, он широко применяется для очистки сточных вод (при использовании песчаных и других фильтров, называемых глубинными из-за большой толщины фильтрующей перегородки). Закономерности фильтрования с закупориванием пор еще недостаточно изучены [9].

Характеристика наиболее распространенных видов фильтрова­ ния с закупориванием пор фильтрующей перегородки приведена в табл. 1.13 [49].

В потоке разбавленной суспензии, проходящем через пористый фильтрующий слой, концентрация твердой фазы постепенно уменьшается. Зависимости, показывающие изменение концентра­ ции через определенные промежутки времени, могут быть выра­ жены двумя способами: с помощью пучков кривых с. = /(т, Л) или

а = f (т, Л), где с — объемная доля твердой фазы в суспензии; а

задержка (накопление) твердой фазы в порах фильтрующей пере­ городки (доля объема слоя, занимаемая осевшими частицами); h — толщина слоя фильтрующей перегородки, м (рис. 1.16).

Для определения вида функций с = /(т. Л) и ст = f (т, Л) необхо­

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

Выделим элемент фильтрующего слоя толщиной Ah. В начальный момент процесса пористость слоя б0 представляет собой

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

72

Таблица 1.13

Особенности процесса фильтрования с закупориванием пор фильтрующей перегородки

Тип

Размер

Места за­

Удержи­

Меха-

Возмож­

Способ

филь­

фильтруе­

держки

ность

принуди­

трова­

мых твердых

твердых

вающие

низм

само-

тельной

силы

захвата

ния

частиц, мкм

частиц

очистки

очистки

Механи­

>30

Сужения,

Трение,

Осажде­

Невоз­

Измене­

ческий

 

расщели­

давление

ние,

можна

ние на­

 

 

ны, пусто­

потока

прямой

 

правления

 

 

ты в слое

 

захват

 

потока

Физико-

А (<3)

Поверх­

Вандерва-

Прямой

Воз­

Увеличе­

химичес­

 

ность

альсовы,

захват

можна

ние ско­

кий

 

пор

электроки-

 

 

рости

 

 

 

нетические

 

 

потока

 

 

 

силы

 

 

 

Колло­

<0,1

Поверх­

Вандерва-

Прямой

Воз­

Увеличе­

идный

 

ность

альсовы,

захват,

можна

ние ско­

 

 

пор

электроки-

диффу­

 

рости

 

 

 

нетические

зия

 

потока

 

 

 

силы, хими-

 

 

 

 

 

 

ческие связи

 

 

 

дении суспензии начинается выделение твердых частиц на стенках пор фильтрующего слоя и пористость е0 уменьшается до б. Долю

объема, занятого осевшими в слое частицами и застойными зонами жидкости, захваченной этими частицами, можно выразить урав­ нением

Суспензия

I + +

Фильтрат

Рис. 1.16. Выходные кривые фильтрования с закупориванием пор филь­ трующей перегородки

1 - 5

кривые, снятые последовательно через равные промежутки времени

73

е0 — е = аа.

(1.192)

Здесь а — задержка твердой фазы; а — коэффициент (а > 1), причем (а — 1) представляет собой объем жидкости, захваченной осевшей твердой фазой.

Другая доля объема £ занята подвижными твердыми частицами (ес, где с — объемная доля твердой фазы в суспензии) и подвижной жидкостью (е(1 — с)).

Таким образом, в элементе фильтрующего слоя с площадью сечения F и толщиной ДЛ будет содержаться следующее количество

твердой фазы:

 

п, = F(a + ес)ДЛ.

(1.193)

Количество твердой фазы, прошедшей через выделенный элемент слоя, составит:

дс

n2 =Fm<pc ~ FD — ' (1.194)

8h

где <оф — фиктивная (в расчете на пустое сечение аппарата) скорость

потока, м/с; со с — конвективный поток твердой фазы;

дс

D — —

9

dh

диффузионный поток твердых частиц (согласно закону Фика); D — коэффициент диффузии, м -1.

Уравнение материального баланса (в условиях равенства ско­ рости накопления твердых частиц в слое и изменения расхода твер­

дой фазы) будет иметь вид

 

 

 

 

д

дс

д2с

,,

— (ст + ес) + сол

D — - =

0.

(1.195)

Ф dh

dh2

 

 

Это уравнение может быть упрощено введением некоторых до­ пущений [67, 94].

Случай 1. В первом приближении можно пренебречь диффу­ зионным потоком, так как в большинстве практических случаев конвективный поток во много раз превышает диффузионный. Тогда уравнение (1.188) примет вид

ддс

— (о + ес) + соа — = 0.

(1.196)

dh

 

Случай 2. Для разбавленных суспензий (с * 0,001) можно также пренебречь изменением пористости слоя и считать s = s0. Тогда

Случай 3. Иногда можно не учитывать число подвижных частиц (е0с) в суспензии, если оно пренебрежимо мало по сравнению с числом частиц, осевших в слое:

Ё 1 + а ф — = 0.

(1.198)

dh

 

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

Продолжительность закупоривания пор. Перемещение фронта задержки твердой фазы в данный момент времени можно харак­ теризовать зависимостью

Е h

(1.199)

т — — = const,

ЙФ

где £0Л/юф представляет собой время, необходимое для перемещения

суспензии сквозь слой. Если процесс фильтрования длится несколь­ ко часов, то этой величиной можно пренебречь, однако в общем случае ее нельзя считать ничтожной (например, при скорости суспензии соф = 0,1 см/с и е0 = 0,4 для прохождения суспензии

через фильтрующий слой толщиной h = 50 см требуется 200 с).

Кинетика процесса закупоривания пор. В течение малого про­ межутка времени Дт в элементе пористого фильтрующего слоя, в котором перемещается ес твердых частиц, задержка твердой фазы

'д а л

будет изменяться в соответствии с выражением — Дт

\дх

Вероятность задержки одной твердой частицы в единицу вре­ мени

1

5ст

к '= ----------

(1.200)

ЕС

сравнима с вероятностью задержки одной частицы при прохож­ дении потока через единицу толщины слоя со скоростью ш = соф/е.

Вероятность задержки частицы на единицу толщины слоя к =

= к ’г/а ., или

ф

74

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

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

к = *0Ф(ст),

(1.202)

где к0 — вероятность закупоривания чистого слоя (в начальном

периоде фильтрования); ф(ст) — функция задержки.

Вид функции задержки ф(ст). Самая простая математическая фор­

ма функции ф(ст) — это закон линейного уменьшения:

 

Ф(ст) = 1 ----- —— ,

(1.203)

где стмакс — максимальное значение задержки (когда все поры на­

сыщены твердой фазой обрабатываемой суспензии и вероятности задержки не существует).

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

Учесть первоначальное увеличение ф(ст) и все возможные виды экспериментальных кривых можно введением функции от трех параметров [14]:

ф(ст) = (1 + Ьст)

[1.204)

 

® макс

где b — положительный коэффициент, учитывающий начальное увеличение задержки; д — показатель степени, характеризующий криволинейность функции.

Общий баланс задержки ст. Если к моменту времени т, поры

фильтрующего пористого слоя частично заполнены твердыми час­ тицами, выделившимися из суспензии, то объем частиц в эле­ ментарном слое между сечениями h и [h + dh) будет равен

и, = (ст + ec)Fdh.

(1.205)

С момента начала фильтрования до момента т через входное сече­ ние Л пройдет следующее количество твердой фазы:

76

 

U2 =

•1

 

(1.206)

 

 

 

На выходе из слоя в сечении

(Л +

dh) за это ж е

время объем

твердой фазы составит

 

 

 

 

 

 

 

Зи2-dh.

(1.207)

 

 

 

dh

 

 

Тогда баланс по твердой фазе для рассмотренного элементар­

ного слоя будет иметь вид

 

 

 

 

 

о,

=

о2 -

о3,

(1.208)

или

 

 

 

 

 

 

 

 

'1

 

 

 

СТ + ЕС =

 

|ш фссгт

(1.209)

 

 

 

dh V0

 

 

Выразим объемную долю твердой фазы с с помощью уравнения

(1.202) и, так как h =

const, получим

 

 

?

w

J

 

f do

( 1.210)

Jco0 cdT =

d r = J — .

 

 

 

k <h

 

 

Следовательно,

 

 

 

 

 

 

5

/

л

1 ( 5сЛ

 

 

 

%8а

( 1.211)

ст + sc =

 

к

к U h J

 

dh

 

 

 

 

\о К )

 

 

или

 

 

 

 

 

 

^астл

 

-k(a + EC).

( 1.212)

 

 

 

ah л

Уравнение (1.212) строго отражает баланс задержки твердой фазы в пористом фильтрующем слое при допущении, что диффузия ничтожно мала (случай 1).

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

77

до

г д сл

(1.213)

=

\dh

ск

 

'дс4

= -кс;

(1.214)

dh S i

 

 

dh S i

= -ко .

(1.215)

 

 

Профили изменения объмной доли твердой фазы и задержки.

Изменения объемной доли твердой фазы и задержки по высоте фильтрующей перегородки в процессе фильтрования с закупори­ ванием пор описываются уравнениями (1.214) и (1.215), имеющими одинаковую математическую форму.

Можно установить соотношение между с и а для данного момента времени т, исключив к\

1

1

ГS a ^ J

(1.216)

f5cl

= ---

laj

с < д х ; т

О

 

После интегрирования получим

 

 

 

lgc =

lga

+ С.

(1.217)

Константа интегрирования С определяется соотношением

С_ _ о

(1.218)

с, о.

где ct и ст. — объемная доля твердой фазы и задержка на входе в

слой соответственно. Задержка на входе возрастает по экспонен­ циальному закону в зависимости от времени и достигает пре­ дельного значения при максимальном насыщении слоя (стмакс * 1 %).

Величина ст. определяется интегрированием кинетического урав­ нения (1.202), записанного для h = 0 в виде

da,

(1.219)

= *0ЮфС,-£*Т.

ф(о,-;

 

При интегрировании имеем

 

da.

( 1.220)

= *0©фС,Т.

Ф К

Для простого случая, когда ф(а,) = 1

-----------

, получим

 

 

 

 

С 7 ._ _

 

 

 

 

( 1.221)

Обозначим

через т , тогда

 

 

 

 

f

_

i

\

 

 

1 - е

т'

( 1.222)

где тс — константа времени для фильтрующей перегородки. Профили объемной доли твердой фазы и задержки (коэффи­

циента объемного осаждения) для ряда значений т/тс при простой

зависимости т(а) = 1

— (а /а макс) представлены на рис. 1.17. Значе­

ния т/т (рис. 1.17, а) вычислены на основе уравнения

 

J

d a

J

da,.

(1.223)

СТФ(СТ)

— = ~koh'

о

о

<VP(CTi)

 

полученного интегрированием дифференциального уравнения ба­ ланса задержки (1.214) при условии, что к = &0ср(а).

Для ф(а) = 1 — (а/ст ) в результате интегрирования получим

 

т/т.

 

(1.224)

л

т /т .

е

knh

1 - е

с -

0

 

Аналогично вычисляют профили задержки для других функций

ф(о).

Рассчитав задержку для ряда значений Ли т , можно определить объемные доли твердой фазы с для тех ж е Л и т , используя со­ отношение (1.219), и получить пучок кривых (рис. 1.17, б), харак­ теризующих продвижение фронта концентрации для выбранной формы функции ф(а).

Скорость перемещения фронта задержки. Профили изменения задержки по высоте фильтрующего слоя для различной продолжи­ тельности фильтрования приведены на рис. 1.18 [123].

78

79