Модель сжимающегося ядра: кинетика выщелачивания с примером
Модель сжимающегося ядра сводит кинетику выщелачивания к одному диагностическому вопросу: какое сопротивление задаёт темп — жидкая плёнка, слой золы или поверхностная реакция. В разобранном примере одна и та же частица оксида меди 10 mm достигает 80% превращения за 4.0 h, 10.0 h или 15.0 h в зависимости от ответа, и этот разброс в 3.7 раза напрямую переходит в время пребывания и размер оборудования.
Ключевые положения
- Три линеаризованных графика — X от t, 1 − 3(1−X)^(2/3) + 2(1−X) от t и 1 − (1−X)^(1/3) от t — выявляют режим.
- Наиболее линейный график в окне 20–80% превращения обычно указывает лимитирующую стадию; его наклон равен 1/τ.
- Диффузия в слое золы даёт τ ∝ R²: крупность дробления — самый сильный рычаг проектирования.
- Плёночный контроль редок при хорошем перемешивании и обычен в плохо смешанных пульпах и застойных зонах кучи.
- Подогнанные τ прямо превращаются в часы пребывания в аппарате или дни слоя в куче.
Что это такое
Модель идеализирует каждую частицу как плотную сферу реактивного минерала в инертной оболочке выщелоченного остатка. Фронт реакции — граница свежего минерала и отработанного остатка — движется к центру, а внешний диаметр сохраняется, поскольку пустая порода не растворяется. Превращение связано с геометрией: X = 1 − (r_c/R)³.
Между объёмом раствора и фронтом последовательно действуют три сопротивления, и наибольшее задаёт скорость. Слой золы — пористая оболочка отработанной породы, через которую реагент должен продиффундировать; плёночный контроль означает, что медленнее всего перенос через ламинарный пограничный слой у поверхности; поверхностный химический контроль означает, что медленнее всего сама реакция разрыва связей.
Для сферической частицы постоянного размера в псевдостационарных условиях каждый режим имеет своё уравнение, где X — доля превращения, τ — время полного превращения:
- Плёнка: X = t/τ
- Слой золы: 1 − 3(1−X)^(2/3) + 2(1−X) = t/τ
- Поверхностная реакция: 1 − (1−X)^(1/3) = t/τ
| Режим | Время полного превращения τ | Чувствительность к размеру | Практический рычаг |
|---|---|---|---|
| Диффузия в плёнке | τ = ρ_B·R / (b·k_l·C_Ag) | τ ∝ R | Перемешивание, вязкость, концентрация |
| Диффузия в слое золы | τ = ρ_B·R² / (6·b·D_e·C_Ag) | τ ∝ R² | Крупность дробления, пористость остатка |
| Поверхностная реакция | τ = ρ_B·R / (b·k_s·C_Ag) | τ ∝ R | Температура, концентрация, тонкость минерала |
Здесь ρ_B — молярная плотность ценного минерала (mol/m³), k_l — коэффициент массоотдачи плёнки (m/s), D_e — эффективная диффузивность в золе (m²/s), k_s — поверхностная константа (m/s), C_Ag — концентрация реагента в объёме (mol/m³). Концентрация стоит в каждом знаменателе: всё, что снижает концентрацию, которую «видит» частица — например, кислотоёмкость породы, — пропорционально удлиняет τ.
Три режима
Плёнка. Частица потребляет реагент быстрее, чем жидкость подводит, концентрация падает в пограничном слое, и реакция голодает. Скорость сильно реагирует на перемешивание, вязкость и загрузку твёрдого: классический признак — бутылочный тест, чья кажущаяся скорость меняется с оборотами мешалки. На заводе это чаще сигнал о проблеме смешения или транспорта, а не свойство руды. Так как вся поверхность питается одинаково, τ линейно зависит от радиуса, а X растёт линейно во времени.
Слой золы. Растворение на фронте быстрое, а узкое место — извилистая миграция реагента через растущую пористую оболочку. Скорость падает с превращением, и кривая заметно выполаживается при высоких X. Квадратичная зависимость τ ∝ R² — самое значимое соотношение в кинетике выщелачивания: уменьшение размера вдвое сокращает время вчетверо, поэтому руды, лимитированные диффузией, окупают каждый шаг измельчения, а лимитированные реакцией — только линейно. Температура влияет слабо, так как энергии активации диффузии (5–20 kJ/mol) много ниже химических (40–80 kJ/mol).
Поверхностная реакция. Медленнее всего сама реакция, концентрация на фронте равна объёмной, перемешивание несущественно. Превращение подчиняется уравнению 1 − (1−X)^(1/3) = t/τ при τ ∝ R и сильной аррениусовской зависимости.
| Режим | Уравнение | Линеаризованный график (y от t) | Наклон |
|---|---|---|---|
| Плёнка | X = t/τ | X | 1/τ |
| Слой золы | 1 − 3(1−X)^(2/3) + 2(1−X) = t/τ | 1 − 3(1−X)^(2/3) + 2(1−X) | 1/τ |
| Поверхностная реакция | 1 − (1−X)^(1/3) = t/τ | 1 − (1−X)^(1/3) | 1/τ |
Есть и чисто геометрический признак, не требующий подгонки: отношение времени до 90% и до 50% превращения. Плёнка даёт 1.8 (0.9/0.5), реакция — 2.6 (0.536/0.206), зола — 5.0 (0.554/0.110). Кривая с сильно растянутым хвостом сообщает о диффузионном контроле ещё до регрессии.
Как подгонять данные
- Пересчитать анализы в превращение. Для каждой пробы X = 1 − (металл в остатке)/(металл в исходной пробе) на сухой основе с поправкой на выход остатка.
- Вычислить три линеаризованных члена. Для каждой пары (t, X): y₁ = X, y₂ = 1 − 3(1−X)^(2/3) + 2(1−X), y₃ = 1 − (1−X)^(1/3).
- Построить графики и проверить линейность. Режим с прямой линией в окне 20–80% и есть контролирующий. Подтвердить R²; систематическая кривизна говорит о неверном или меняющемся режиме.
- Извлечь τ из наклона. τ = 1/наклон в единицах оси времени (наклон 0.10 h⁻¹ даёт τ = 10 h). Указывать τ вместе с фракцией крупности и условиями: τ без размера бессмысленно.
- Проверить хвост. Кривизна выше 90% нормальна для диффузии в золе и патологична для реакции. Прямой участок вначале и крутой в конце часто означает переход от реакции к диффузии по мере роста оболочки.
- Подтвердить независимыми рычагами. Повтор при половинном радиусе: τ падает вдвое — плёнка или реакция; вчетверо — диффузия в золе. Повтор при двух температурах: ускорение по Аррениусу подтверждает химический контроль, нечувствительность — диффузионный. Два согласующихся рычага сильнее одного хорошего R².
Разобранный пример: частица 10 mm
Пример иллюстративный: числа согласованы между собой и показывают метод, но не являются измерениями конкретной руды. Бутылочные тесты на узкой фракции 10 mm (R = 5.0 mm) в 0.5 mol/L H₂SO₄ при 25 °C и 600 rpm дают τ_плёнка = 5 h, τ_реакция = 24 h и τ_зола = 40 h.
Шаг 1 — плёнка. X = t/τ, значит при X = 0.80: t = 5 × 0.80 = 4.0 h.
Шаг 2 — слой золы. 1 − X = 0.20; (0.20)^(1/3) = 0.5848, (0.20)^(2/3) = 0.3420. Тогда t/τ = 1 − 3(0.3420) + 2(0.20) = 0.374 и t = 0.374 × 40 = 14.96 h ≈ 15.0 h.
Шаг 3 — поверхностная реакция. t/τ = 1 − 0.5848 = 0.415 и t = 0.415 × 24 = 9.97 h ≈ 10.0 h.
| Превращение X | Плёнка t = 5·X (h) | Зола t = 40·[1 − 3(1−X)^(2/3) + 2(1−X)] (h) | Реакция t = 24·[1 − (1−X)^(1/3)] (h) |
|---|---|---|---|
| 0.25 | 1.25 | 0.94 | 2.19 |
| 0.50 | 2.50 | 4.40 | 4.95 |
| 0.80 | 4.00 | 14.96 | 9.97 |
| 0.95 | 4.75 | 27.71 | 15.16 |
| 0.99 | 4.95 | 35.23 | 18.83 |
Таблица показывает две особенности, проявляющиеся в любом тесте. Во-первых, порядок трудности меняется: при низком превращении кривая реакции самая медленная (2.19 h при X = 0.25 против 0.94 h у золы), потому что оболочка ещё тонка, но к X = 0.80 диффузия становится самым жёстким режимом. Во-вторых, хвосты различаются радикально: последние 19% стоят 20.3 h при диффузии против 8.9 h при химическом контроле.
Рычаг размера. При той же руде 5 mm (R = 2.5 mm) и диффузии в золе τ падает с 40 до 10 h, а 80% достигается за 0.374 × 10 = 3.7 h — четырёхкратное улучшение от половинного диаметра. При химическом контроле τ падает до 12 h, а время до 5.0 h: всего вдвое. Эта асимметрия — количественное ядро экономики дробления.
Где модель отказывает
Несферические и трещиноватые частицы. Обломки угловатые, пластинчатые, трещиноватые; фронт идёт неравномерно, и факторов формы отклоняется на десятки процентов при высоких X. Указывайте τ как кажущуюся величину для этой фракции и формы, а не как константу материала.
Переменная диффузивность. Модель предполагает постоянную D_e, но реальные оболочки меняются: гипс, ярозит или гидроксид железа, осаждаясь в порах, снижают D_e, τ дрейфует, а график изгибается. Сам дрейф служит диагностикой проблем химии раствора.
Пористые руды и модель зёрен. Многие промышленные руды — пористые агрегаты микронных зёрен. Модель зёрен описывает такую частицу как тысячи миниатюрных ядер, и следствие огромно: длина диффузионного пути становится радиусом зерна, а не частицы, поэтому штраф R² практически исчезает и грубое дробление может выщелачиваться почти с реакционным контролем. Применение модели плотной частицы к такой руде переоценивает выгоду тонкого помола.
Сжимающиеся частицы. Если матрица тоже растворяется, инертного остатка нет и нужен вариант сжимающейся частицы; частично растворимая порода лежит между случаями.
Эффекты слоя. SCM описывает частицу в свежем растворе; куча или колонна добавляет байпас, каналирование, гистерезис смачивания и профили концентрации. Подгонка по частице — правильная отправная точка, но масштабирование требует гидрологии слоя: поэтому проект куч опирается на колонные, а не только бутылочные тесты.
От кинетики к размеру оборудования
Аппараты с мешалкой. С реакционным τ = 24 h при R = 5 mm и целью 90%: t = 24 × [1 − (0.10)^(1/3)] = 24 × 0.5358 = 12.9 h идеального периодического времени. Непрерывные каскады никогда не дают идеального вытеснения — байпасирование стоит примерно 25–35% теоретической кинетики при типовом числе ступеней, — поэтому проект при КПД ступени 70% требует 12.9/0.70 ≈ 18.4 h общего пребывания, обычно в трёх аппаратах по 6 h. Две оговорки: подгоняйте τ при проектной плотности пульпы, так как высокая загрузка твёрдого может сдвинуть реакционную систему к плёночному контролю; и переподгоняйте после любого изменения крупности, поскольку при τ ∝ R более грубое питание незаметно растягивает каждый аппарат.
Кучи. Параметры приходят из колонн, и подогнанный режим диктует реакцию: контроль слоем золы (τ ∝ R²) требует более тонкого дробления и более высоких слоёв с удлинёнными циклами; реакционный контроль говорит, что размер важнее не столько, сколько контакт с раствором и подача реагента. Кривая задаёт и экономическое отсечение: с приведёнными числами переход от X = 0.80 к 0.95 при диффузионном контроле стоит 12.8 дополнительных суток на слой — выбор, который надо делать явно, а не по умолчанию.
Часто задаваемые вопросы
Какие допущения делает SCM?
Сферические частицы постоянного внешнего размера с инертным скелетом, резкий фронт реакции, псевдостационарная диффузия и одна лимитирующая стадия среди плёнки, слоя золы и поверхностной реакции. Концентрация в объёме считается однородной. Отклонения — угловатые частицы, схлопывающаяся пористость, смешанный контроль — проявляются кривизной линеаризованных графиков.
Какой режим обычно контролирует промышленное выщелачивание?
В хорошо перемешиваемых аппаратах с умеренно тонкой рудой доминируют поверхностная реакция или диффузия в слое золы; плёночный контроль редок, потому что перемешивание держит пограничный слой тонким. В кучах и на крупной руде доминирует диффузия. Режим может меняться с превращением — от реакции вначале к диффузии в конце.
Как экспериментально отличить плёнку от слоя золы?
Меняя перемешивание и температуру. Плёнка сильно реагирует на скорость перемешивания и слабо на температуру из-за низких энергий активации переноса. Слой золы тоже слабо чувствителен к температуре, но не зависит от перемешивания, а его τ масштабируется квадратом размера. Матрица двух размеров и двух температур разделяет режимы четырьмя опытами.