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

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

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

Аннотация научной статьи по наукам о Земле и смежным экологическим наукам, автор научной работы — Михайловская Е. Ю.

Рассмотрена математическая модель процессов фильтрации при разработке уранового месторождения методом подземного выщелачивания. Фильтрация рассматривается в приближении двумерной плановой модели. Исходные данные для моделирования определялись из анализа данных разведочного бурения на Хохловском месторождении. Уравнение фильтрации решалось численно методом последовательной верхней релаксации. Режимные параметры процесса выщелачивания определялись в соответствии с планируемыми расходами реагента в закачивающих скважинах.A mathematical model of flow is considered for development of uranium ore deposit by underground leaching.

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

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

УДК 622.234.42 Е.Ю. Михайловская

Российский химико-технологический университет им. Д.И. Менделеева, Москва, Россия

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

A mathematical model of flow is considered for development of uranium ore deposit by underground leaching. The flow is examined in approximation of 2-D areal model. Input data for the modeling are determined from an analysis of exploratory drilling at the ore deposit Khokhlovskoye. The flow eauation is solved numerically by the of successive overrelaxation method. Operation mode parameters of the leaching process are determined according to planned rate of the reagent in the injection wells.

Рассмотрена математическая модель процессов фильтрации при разработке уранового месторождения методом подземного выщелачивания. Фильтрация рассматривается в приближении двумерной плановой модели. Исходные данные для моделирования определялись из анализа данных разведочного бурения на Хохловском месторождении. Уравнение фильтрации решалось численно методом последовательной верхней релаксации. Режимные параметры процесса выщелачивания определялись в соответствии с планируемыми расходами реагента в закачивающих скважинах.

Одна из важнейших проблем, от решения которой зависит устойчивое развитие человечества, - обеспечение стабильно увеличивающегося производства энергии [1]. Только за период с 1860 по 1985 г. поток энергии, потребляемый мировой экономикой, возрос в 60 раз. Сегодня мировое потребление продолжает неравномерно, но неуклонно расти [2]. В связи с уменьшением запасов углеводородных энергоносителей для решения этой проблемы в концепциях развития энергетики многих стран мира предполагается увеличить долю выработки энергии на атомных электростанциях. Для этого необходимо развитие урановой сырьевой базы и разработку экономически эффективных и экологически безопасных технологий добычи сырья для ядерных энергетических реакторов. Перед российскими учёными и специалистами производственных организаций поставлена задача в ближайшие 10 лет значительно увеличить добычу урана. Принятая в 1997 г. Федеральная целевая программа «Уран России» ориентированна на поиски, ускоренную разведку и освоение новых месторождений урана, пригодных для разработки экологически безопасными методами подземного выщелачивания, это ещё раз подтверждает актуальность данной работы. [3].

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

В промышленности в зависимости от минерального состава вмещающих пород и типа рудной минерализации используют два основных способа подземного выщелачивания урана: сернокислотный (широко применявшийся в СССР) и бикарбонатный, (распространённый в США), в качестве окислителя применяют кислород, пероксид водорода, трёхвалентное железо. Высокая эффективность способа привела к снижению требований промышленности к качеству уранового сырья (минимальное бортовое содержание может быть снижено с 0,01% до 0,005%, тогда как для месторождений, отрабатываемых горным способом ~ 0,03%). Геологические и гидрологические условия локализации по-лиэлиментного оруденения на месторождениях такого типа преимущественно благоприятны для добычи полезных компонентов подземным выщелачиванием [3].

Сернокислотное выщелачивание:

и02,25(к)+804 (р-р)+2Н (р-р)+0,37502(г)-► и02804 (р-р)+Н2О(ж),

и8Ю4(к)+2804 (р-р)+ 2Н (р-р)+О2(г) +Н2О(ж) ^ И02804 (р-р) +Н2БЮ4(р-р)

Бикарбонатное выщелачивание:

и02,25(к)+3С0з2-(р-р)+2Н+(р-р)+0,37502(г) Ш2 (С0з)з4- (р-р)+Н2О(ж),

ШЮ4(к)+2С0з2-(р-р)+ 2Н+(р-р)+0,502(г) —► И02(С0з)з4- (р-р) +§1°2(к) +Н2О(ж),

В кислой среде окислителем может выступать Бе0Н , Бе(0Н) 2 , Бе2(0Н) 2 и Без+. При этом реакция выщелачивания запишется в виде [3]:

и02(к)+ Бе(0Н)2 + (р-р) и02 + (р-р)+ Бе(0Н) 2(р-р), И0+ 2(р-р)+ Ре(0Н)2+(р-р) —► и02^+ (р-р)+ Бе(0Н) 2(р-р), Бе(0Н)2(р-р) +2Н (р-р) Бе (р-р)+2Н2О(ж) . Хохловское месторождение расположено на стыке Зауральского поднятия и Тюменско-Кустанайского прогиба в Южной части Среднего Зауралья. Всего здесь три месторождения урана - Далматовское, Хохловское и Добровольское, общие запасы которых оцениваются в 30 - 35 тыс. тонн урана. Рудное поле Хохловского месторождения имеет протяжённость 28 км и расположено в пределах самого южного притока Песчанской палеосистемы, имеющего простирание с юго-запада на северо-восток, прослеженного на 65 км. Глубина залегания продуктивной толщи составляет 520 - 680 м.

Главным урановым минералом является настуран. Настуран присутствует в виде глобул, иногда полых, с размерами 0,0п-0,п мкм, агрегатов и цепочек глобул размером до нескольких микрон, пленок, корочек с трещинами усыхания, почек размером от 0,п до 10-20 мкм. Богатые урановые руды, с содержанием урана свыше 0,1 %, локализованы в темно-серых песчаниках. Руды с содержанием урана от 0,01% до 0,1 % встречаются практически во всех литологических и геохимических разностях, в сероцветных и обеленных породах [4].

Метод подземного выщелачивания характеризуется такими параметрами как лёгкая вскрываемостью пород и высокой степенью извлечения урана. Например, бикарбонатные (содержание СО2 ниже 1%) руды с исходными содержаниями урана 0,01-0.1% обладают лёгкой вскрываемостью растворами серной кислоты (5-Э0 г/л);степень извлечения урана из руд высокая (85-97%), удельные расходы реагента низкие (25-100 кг/кг урана или 5-15кг/т руды). Для слабокарбонатных руд (содержание СО2 ниже 1-2,5%) расходы реагента возрастают до 30 кг/т руды. Дополнительные характеристики такие как отсутствие отвалы забалансовых руд и пустых пород и провалов, увеличивают вероятность его применения при добычи урана. При таких параметрах важно точно рассчитать необходимое количество реагента, в противном случае может произойти загрязнение подземных вод остаточными сернокислотными растворами, содержащими экологически вредные вещества.

Существует два способа восстановления качества подземных вод: 1) самопроиз-

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

Наиболее обоснованным способом расчета течения реагента через породы пласта являются методы математического моделирования. Для этого была рассмотрена математическая модель фильтрации растворов в приближении двумерной плановой модели [5]. Компоненты скорости фильтрации ух и \у в декартовых координатах х, у удовлетворяют уравнению неразрывности

1 д{ухн) 1 дМ_

и уравнениям Дарси

- + = Q( x, y)

h dx h dy

rdH rdH

v = -f-, v = -f-.

x J ^ ? y J^

ox dy

где / - коэффициент фильтрации; Н - величина напора в пласте; И - мощность пласта; О - плотность источников и стоков.

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

1 d (irdHЛ 1 d (lxdHЛ

--1 hf-1 +--hf-

h dx ^ dx J h dy I dy

= -Q.

V

Рассматривалась прямоугольная область: X < х < Х2, ^ < у < Г2. Для решения уравнения фильтрации вводилась регулярная конечно-разностная сетка: х = X +Ах а -1), Уу = Г1 + А у (У -1); 1 = 1,., N; У = 1,-, N. Здесь Ах, Ау - расстояния между соседними узлами по координатам х и у соответственно.

Для решения уравнения фильтрации использовался модифицированный метод последовательной верхней релаксации [6]

H+1 = H/j + +i,j(h/+w -H/j)-Г4,,(h/J -H£j)+

+ b2r ! H" - H" )-r ! (H" - H"+1 )1+—A A h Q , 1,

Lz,j+1\ г,j +1 г,j/ г,j-1\ г,j г,j-1 /J т- x y г,j г,j [

max

r/

где b = Ax / A^; H/. - значение И(хг, y) на итерации n; h . = h(xt, y ■), ■ = Q(x, y.) r,j = (h,fj)/Fjmax , f, = f(x,,у,), Fax = max fax,y)f(x,y)}.

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

x = Х1г H = Hj( y); x = X2, H = H 2( у);

y = H = Hз( x); y = Y2, H = H4(x).

Значения функций Н\, Н2, Н3, Н4 определялись для невозмущенного регионального течения подземных вод в пласте. Величины h(x, y) определялись геостатистическими методами на основании локальных значений h, полученных из анализа керна при разведочном бурении. Величина f определялись по результатам тестовых откачек в скважинах. Из

анализа полученных значений / следует, что в исследуемом пласте величину / приближенно можно считать постоянной. Плановое распределение И(х, у) приведено на рис. 1.

250

200

1 50

1 00-

50

0

300

Рис. 1. Плановое распределение к(х, (пояснение в тексте)

Расчеты распределений напоров в пласте проводились для рекомендованных в [4] значений расходов реагента в закачивающих скважинах и интенсивности откачки в откачивающих скважинах.

250

200-

150-

100

50

0

Н,

м

6 4 2 1 0 -1 -2 -4 -6 -8 -1 0 -20

0

50

100

150

200

250

300

х, м

Рис. 2. Распределение напоров в пласте(пояснение в тексте)

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

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

1. Тарасова, Н.П. Устойчивое развитие: ресурсы России /Н.П.Тарасова, Н.П.Лаверов, Е.Б.Кручина и др. - М.: Изд центр РХТУ. 2004. - 212 с.

2. Лавёров, Н.П. Подземное выщелачивание полиэлиментных руд/ Лавёров Н.П., И.Г.Абдульманов, К.Г.Бровин.- М.: Издательство Академии горных наук, 1998. -446 с.

3. Медоуз, Д.Х. За пределами роста/ Д.Х.Медоуз, Д.Л.Медоуз, Й.Рандерс.-М.: Прогресс. -1994.-304с.

4. Марков, С.Н. Геологоразведочные работы по оценке промышленной значимости Хохловского месторождения урана/ С.Н.Марков, С.И.Долбилин.- с. Петропавлоское.-2003.- 173с. Государственный регистрационный №73-00-1/6

5. Мироненко, В.А. Динамика подземных вод / В.А.Мироненко.- М.: МГГУ, 1996.-519 с.

6. Роуч, П. Вычислительная гидродинамика /П. Роуч. - М.: Мир, 1980. - 616с.

УДК 66.074.3 О.И. Козинская

Российский химико-технологический университет им. Д.И. Менделеева, Москва, Россия

ОЧИСТКА ОТХОДЯЩИХ ГАЗОВ ОТ НЕОРГАНИЧЕСКИХ И ОРГАНИЧЕСКИХ КОМПОНЕНТОВ

A technological scheme is worked out for purifying exhaust gases containing ashes, soot, mixture of polycyclic, aromatic, saturated and un saturated hydrocarbons, NOx, SO2, HCl and other components.

Разработана технологическая схема для очистки отходящих газов, содержащих золу, сажу, смесь полициклических, ароматических, предельных и непредельных углеводородов, а также NOx, SO2, HCl и других компонентов.

Большую сложность вызывает очистка отходящих газов от смеси различных неорганических и органических компонентов, содержащих, например, золу, сажу, смесь полициклических, ароматических, предельных и непредельных углеводородов, а также оксидов азота и серы, хлористого водорода и других компонентов при температуре 400-600оС, поскольку нет единого метода очистки для их одновременного обезвреживания [1].

Очистка таких газов, как правило, происходит в несколько стадий (ступеней): пылеочистка, термическое окисление органических компонентов, абсорбционная очистка (NOx, SO2, HCl).

Пылеочистка. Традиционно пылеочистка осуществляется в электрофильтрах и рукавных фильтрах. В некоторых случаях используются специальные рукавные фильтры из графитизированной ткани или металлокерамические фильтры при высокой температуре отходящих газов (600 - 1000оС), когда понизить температуру не представляется возможным. В приведенном примере (температура 400 - 600оС) за пылеочисткой должна следовать термокаталитическая очистка, для которой нужна температура 350 -600оС, поэтому понижать температуру для пылеочистки не желательно и в данном случае будут использоваться специальные фильтры.

Термокатализ. Очистка от смеси полициклических, ароматических, предельных и непредельных углеводородов осуществляется термическими, термокаталитическими и адсорбционными методами. Наиболее подходящим методом для этого случая является термокатализ. Этот метод обезвреживания органических компонентов отходящих газов имеет значительные преимущества над другими методами (адсорбцией, термическим дожиганием). При термическом дожигании температура достигает 1000 - 1300оС и выше. В отличие от термического дожигания данный способ обеспечивает высокую степень очистки (свыше 99%) при более низких температурах (350 - 600оС), а следовательно, требует меньших затрат энергии, да и с конструкционными материалами проще. Каталитические методы дают возможность обезвреживать органические соединения различного содержания и концентрации. В связи с ростом стоимости энергоресур-

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