Научная статья на тему 'МЕТОДЫ И МОДЕЛИ ИДЕНТИФИКАЦИИ ЗАПАСОВ ТЕХНОЛОГИЧЕСКИХ БЛОКОВ РУДНИКА ПОДЗЕМНОГО СКВАЖИННОГО ВЫЩЕЛАЧИВАНИЯ УРАНА'

МЕТОДЫ И МОДЕЛИ ИДЕНТИФИКАЦИИ ЗАПАСОВ ТЕХНОЛОГИЧЕСКИХ БЛОКОВ РУДНИКА ПОДЗЕМНОГО СКВАЖИННОГО ВЫЩЕЛАЧИВАНИЯ УРАНА Текст научной статьи по специальности «Энергетика и рациональное природопользование»

CC BY
117
26
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ИДЕНТИФИКАЦИЯ ЗАПАСОВ / РЕГРЕССИОННЫЕ ВЗАИМОСВЯЗИ КОЭФФИЦИЕНТА РАДИОАКТИВНОГО РАВНОВЕСИЯ / СТАТИСТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РУДНЫХ ИНТЕРВАЛОВ / ПРИРОСТ ЗАПАСОВ / РЕСУРСНЫЙ ПОТЕНЦИАЛ МЕСТОРОЖДЕНИЯ / INVENTORY IDENTIFICATION / RADIOACTIVE EQUILIBRIUM COEFFICIENT REGRESSION RELATIONSHIPS / ORE INTERVALS STATISTICAL MODELING / RESERVES INCREASE / THE DEPOSIT RESOURCE POTENTIAL

Аннотация научной статьи по энергетике и рациональному природопользованию, автор научной работы — Подрезов Денис Рустамович

Содержания урана в пластово - инфильтрационных месторождениях достаточно низкое и составляет десятые и сотые доли процента. Такой тип месторождений целесообразно разрабатывать без извлечения руды на поверхность. Одним из эффективных способов добычи урана в таких месторождениях является подземное скважинное выщелачивание. Многие предприятия урановой отрасли, ведущие отработку месторождений данным способом, сталкиваются с проблемой несоответствия геологических и непосредственно вскрытых запасов металла на технологических добычных блоках. Такие несоответствия связаны прежде всего с результатами интерпретации геофизических данных. В свою очередь это затрудняет прогноз геолого - технологических параметров эксплуатационных объемов бурения скважин при вводе блоков в эксплуатацию; расходов материалов и реагентов на закисление и последующее выщелачивание. В работе проведен сравнительный анализ геофизических методов оценки запасов технологического блока рудника подземного скважинного выщелачивания урана. Рассмотрены методы кернового геологического опробования и специфические геофизические методы определения параметров геологических запасов урановых руд, такие как гамма - каротаж и метод мгновенных нейтронов деления. Для данных методов определены специфические особенности статистической интерпретации получаемых данных и обоснованы причины систематических и случайных (по статистическим критериям) расхождений в параметрах рудных интервалов для различных классов руд и морфологических элементов рудной залежи. Установлено, что качество результатов каротажа удовлетворительное, а сами результаты обладают высокой степенью достоверности и могут быть использованы для количественной оценки коэффициента радиоактивного равновесия рудных интервалов. Показаны возможности формирования по статистической базе данных разведочных и технологических скважин регрессионных взаимосвязей коэффициента радиоактивного равновесия от среднего содержания радия в рудном интервале; определения на этой основе численных значений параметров рудных интервалов для мешковых и крыльевых элементов рудной залежи. Полученные результаты позволили провести статистическое моделирование и определить диапазоны коррекции мощности рудных интервалов для различных морфологических границ на основе выявления трендов графиков бортовых кондиций.

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

Похожие темы научных работ по энергетике и рациональному природопользованию , автор научной работы — Подрезов Денис Рустамович

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

METHODS AND MODELS OF IDENTIFICATION OF RESERVES OF TECHNOLOGICAL UNITS OF URANIUM WELL LEACHING MINE

Uranium content in formation and infiltration deposits is rather low and amounts to tenths and hundreds percent. This deposit type is advantageous to develop without extracting ore to the surface. One effective way to extract uranium in such fields is underground well leaching. Many uranium industry enterprises working in this way are faced with the geological mismatch problem and directly exposed metal reserves on process mining blocks. Such inconsistencies are primarily related to the geophysical data interpretation results. In turn, this makes it difficult to predict geological and technological parameters well drilling operational volumes during the units commissioning; materials consumption and reagents for acidification and subsequent leaching. The work carried out a comparative analysis of underground well leaching uranium mine technological unit reserves estimating geophysical methods. Core geological testing methods and uranium ores geological reserves determining parameters specific geophysical methods, such as gamma logging and instantaneous fission neutrons method are considered. Obtained data statistical interpretation specific features are defined for these methods and systematic reasons and random (according to statistical criteria) discrepancies in ore intervals parameters for different ores classes and ore deposit morphological elements are justified. It has been found that the logging results quality is satisfactory and the results are highly reliable and can be used to quantify the ore intervals radioactive equilibrium coefficient. Themselves formation possibilities regression radioactive equilibrium coefficient relationships from radium average content in ore interval and ore intervals parameters numerical values determination for bag and ore deposit wing elements on this basis are shown according to statistical exploration database and technological wells. The obtained results made it possible to carry out statistical modeling and to determine ore intervals power correction ranges for different morphological boundaries on the on - board conditions graphs trends determination basis.

Текст научной работы на тему «МЕТОДЫ И МОДЕЛИ ИДЕНТИФИКАЦИИ ЗАПАСОВ ТЕХНОЛОГИЧЕСКИХ БЛОКОВ РУДНИКА ПОДЗЕМНОГО СКВАЖИННОГО ВЫЩЕЛАЧИВАНИЯ УРАНА»

12. NTechLab. Available at: https://ntechlab.ru/ (accessed 09.03.2020).

13. Iqbal U. and Gall J. Multi-person pose estimation with local joint-to- person associations. ECCV Workshops, Crowd Understanding. 2016. Available at: https://arxiv.org/pdf/1608.08526.pdf.

14. Watcom. Available at: https://www.watcom.ru/products/sistema_podscheta_posetiteley /?utm_source=google&utm_medium=cpc&utm_campaign=%7Bcampaign_id%7D&utm_content=%7Bad_id%7D&ut m_term=%2BBHgeoKaMepa%20%2BnogcHer%20%2BnocerHreneH&gclid=Cj0KCQjw1Iv0BRDaARIsAGTWD1uK 7s0ZfcIt_HJ2Wzp_NaLnp77xNvvaCKQcXzL4grFqpR0RZTJa040aAuU7EALw_wcB (accessed 31.03.2020).

15. Zhe Cao (2017). Realtime Multi-Person 2D Pose Estimation using Part Affinity Fields. Robotics Institute Carnegie Mellon University. Pittsburgh, Pennsylvania, 2017. Available at: http://www.ri.cmu.edu/wp-content/uploads/2017/04/thesis.pdf.

DOI 10.21672/2074-1707.2020.50.2.032-043 УДК 004.94

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

Статья поступила в редакцию 01.05.2020, в окончательном варианте — 05.06.2020.

Подрезов Денис Рустамович, Национальный исследовательский технологический университет «МИСиС», 119049, Российская Федерация, г. Москва, Ленинский проспект, 4,

соискатель, e-mail: dpodrezov81@rambler.ru

Содержания урана в пластово-инфильтрационных месторождениях достаточно низкое и составляет десятые и сотые доли процента. Такой тип месторождений целесообразно разрабатывать без извлечения руды на поверхность. Одним из эффективных способов добычи урана в таких месторождениях является подземное скважинное выщелачивание. Многие предприятия урановой отрасли, ведущие отработку месторождений данным способом, сталкиваются с проблемой несоответствия геологических и непосредственно вскрытых запасов металла на технологических добычных блоках. Такие несоответствия связаны прежде всего с результатами интерпретации геофизических данных. В свою очередь это затрудняет прогноз геолого-технологических параметров эксплуатационных объемов бурения скважин при вводе блоков в эксплуатацию; расходов материалов и реагентов на закисление и последующее выщелачивание. В работе проведен сравнительный анализ геофизических методов оценки запасов технологического блока рудника подземного скважинного выщелачивания урана. Рассмотрены методы кернового геологического опробования и специфические геофизические методы определения параметров геологических запасов урановых руд, такие как гамма-каротаж и метод мгновенных нейтронов деления. Для данных методов определены специфические особенности статистической интерпретации получаемых данных и обоснованы причины систематических и случайных (по статистическим критериям) расхождений в параметрах рудных интервалов для различных классов руд и морфологических элементов рудной залежи. Установлено, что качество результатов каротажа удовлетворительное, а сами результаты обладают высокой степенью достоверности и могут быть использованы для количественной оценки коэффициента радиоактивного равновесия рудных интервалов. Показаны возможности формирования по статистической базе данных разведочных и технологических скважин регрессионных взаимосвязей коэффициента радиоактивного равновесия от среднего содержания радия в рудном интервале; определения на этой основе численных значений параметров рудных интервалов для мешковых и крыльевых элементов рудной залежи. Полученные результаты позволили провести статистическое моделирование и определить диапазоны коррекции мощности рудных интервалов для различных морфологических границ на основе выявления трендов графиков бортовых кондиций.

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

METHODS AND MODELS OF IDENTIFICATION OF RESERVES OF TECHNOLOGICAL UNITS OF URANIUM WELL LEACHING MINE

The article was received by the editorial board on 01.05.2020, in the final version — 05.06.2020.

Podrezov Denis R., National University of Science and Technological "MISiS", 4 Leninskiy Prospekt, Moscow, 119049, Russian Federation,

applicant, e-mail: dpodrezov81@rambler.ru

Uranium content in formation and infiltration deposits is rather low and amounts to tenths and hundreds percent. This deposit type is advantageous to develop without extracting ore to the surface. One effective way to extract uranium in such fields is underground well leaching. Many uranium industry enterprises working in this way are faced with the geological mismatch problem and directly exposed metal reserves on process mining blocks. Such inconsistencies are primarily related to the geophysical data interpretation results. In turn, this makes it difficult to predict geological and technological parameters well drilling operational volumes during the units commissioning; materials consumption and

reagents for acidification and subsequent leaching. The work carried out a comparative analysis of underground well leaching uranium mine technological unit reserves estimating geophysical methods. Core geological testing methods and uranium ores geological reserves determining parameters specific geophysical methods, such as gamma logging and instantaneous fission neutrons method are considered. Obtained data statistical interpretation specific features are defined for these methods and systematic reasons and random (according to statistical criteria) discrepancies in ore intervals parameters for different ores classes and ore deposit morphological elements are justified. It has been found that the logging results quality is satisfactory and the results are highly reliable and can be used to quantify the ore intervals radioactive equilibrium coefficient. Themselves formation possibilities regression radioactive equilibrium coefficient relationships from radium average content in ore interval and ore intervals parameters numerical values determination for bag and ore deposit wing elements on this basis are shown according to statistical exploration database and technological wells. The obtained results made it possible to carry out statistical modeling and to determine ore intervals power correction ranges for different morphological boundaries on the on-board conditions graphs trends determination basis.

Keywords: inventory identification; radioactive equilibrium coefficient regression relationships; ore intervals statistical modeling; reserves increase, the deposit resource potential

Graphical annotation (Графическая аннотация)

Введение. Урановые месторождения, пригодные для отработки методом подземного выщелачивания, относятся к экзогенным (образованным поверхностными процессами при участии воздуха, поверхностных или грунтовых вод) пластово-инфильтрационным (образовавшимся в результате отложения рудных веществ из водных растворов) месторождениям. Содержание урана в рудных телах в среднем 0,07-0,08 %, глубина залегания - от 300 до 700 м. Основные запасы локализуются в высокопроницаемых песчаных и гравийно-песчаных отложениях. Рудовмещающие водоносные горизонты обычно ограничены достаточно выдержанными по площади и мощности водо-упорами. Такие геологические условия предопределяют выбор экономически эффективного способа отработки таких типов месторождений - подземное скважинное выщелачивание. Однако технолого-экономические вопросы разработки таких месторождений в литературе исследованы слабо. Поэтому целью данной работы является попытка устранить указанный недостаток.

Общая характеристика предметной области. Процесс добыча урана состоит из следующих составляющих: подземное скважинное выщелачивания урана (ПСВУ) сернокислотными растворами на геотехнологическом полигоне скважин (ГТП); насосный раствороподъём продуктивных растворов из откачных скважин; сбор продуктивных растворов в технологических узлах приёма (УППР); транспортировка урансодержащих растворов по сборному коллектору в пескоотстойники продуктивных растворов (ПР) основной промышленной площадки рудника ПСВ; подкисление возвратных растворов в технологических узлах приготовления выщелачивающих растворов (УПВР) серной кислотой, подаваемой с расположенного на промышленной площадке рудника ПСВ склада; подача рабочих (выщелачивающих) растворов в закачные скважины добычных полигонов через узлы распределения (УРВР) (рис. 1-2).

Добычной комплекс

Перерабатывающий комплекс - пром. площадка Рудника ПСВ

Серная кислота со склада

Цех по переработке продуктивных растворов (ЦППР)

Цех аффинажа

Л

Аффинажное производство

Готовая продукция закись-окись

L

Пескоотстойник ВР Пескоотстойник ПР

А

Центральная насосная станция Центральная насосная станция

V

Маточники сорбции Сорбция

А

Товарный десорбат Регенерация

j

Рисунок 1- Укрупнённая технологическая схема промышленной отработки запасов урана на месторождении

Рисунок 2 - Схема добычи урана способом подземного скважинного выщелачивания

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

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

Неоднородность морфологических элементов залежи приводит к значительным погрешностям в интерпретации геофизических данных и, как следствие, неточностям оценки вскрытых запасов металла на месторождении. Кроме того, возникают трудности в адекватном экономическом планировании и бюджетировании основных технико-экономических показателей отработки месторождения. В итоге в оперативном плане вышеперечисленные факторы могут стать причиной неадекватного снижения содержания урана в продуктивных растворах и, как следствие, падения объемов добычи; наличия проблем с выполнением производственных обязательств [2, 3].

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

Ниже мы рассмотрим данный спектр проблем, являющихся типовыми для пластово-инфильтрационных месторождений уранодобывающей отрасли промышленности на примере одного из крупных горнодобывающих предприятий. Оно отрабатывает месторождение урана способом подземного скважинного выщелачивания на юге Казахстана. На первом этапе исследования на уранодобывающем предприятии был проведен первичный анализ показателей вскрытия технологических блоков и выявлены перспективные рудные тела; намечены планы по их отработке и приросту запасов с применением существующей инфраструктуры [1].

Предлагаемые методы. Для определения показателей оценки запасов технологического блока при добыче методом ПСВУ используется несколько методов. Среди них необходимо выделить керновое геологическое опробование (ГО), гамма-каротаж (ГК), метод мгновенных нейтронов деления (КНД-М) [5, 6].

Полученные и оцифрованные данные сводятся в единую базу для подсчета запасов. Она содержит следующие показатели: эффективная мощность закисляемого продуктивного горизонта (Мэ, м); средняя рудная мощность (т, м); среднее содержание урана в руде (с, %); средняя продуктивность рудного тела (ти, кг/м2); площадь технологического блока ^бл, тыс. м2); объем вскрытых запасов технологического блока т).

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

Одной из задач настоящей работы является оценка возможности применения результатов КНД-М (каротаж методом мгновенных нейтронов деления) для уточнения качества и надежности коэффициента радиоактивного равновесия (Крр), полученного при выполнении разведочных работ, и оценки полученных зависимостей Крр от содержания радия в рудном интервале (РИ). Для решения данной задачи была выполнена оценка качества результатов каротажа КНД-М согласно основным положениям Инструкции [6] при ведении разведочных и горно-подготовительных работ.

Результаты и их обсуждение. Достоверность каротажа КНД-М может быть оценена сопоставлением с результатами определения параметров РИ по контрольному методу геологического опробования. Для решения задачи оценки достоверности КНД-М в сопоставлении с контрольным методом (геологическим опробованием) был проведен расчет данных с определением статистических характеристик достоверности результатов КНД-М по отношению к контрольному методу (ГО) (табл. 1). Приведенные в таблице результаты свидетельствуют о достоверности определения параметров РИ методом КНД-М: систематические расхождения в параметрах РИ отсутствуют, случайные расхождения не превышают допустимые. Расчетные значения ^статистики Стьюдента для мощности рудного интервала, метропроцента и содержания урана соответственно равны 1,55, 0,33 и 0,57. Метропроцент - это произведение мощности рудного интервала (м) на содержание в нем полезного компонента (%). Это меньше критического значения 2,02 для Р = 95 %, а относительные расхождения средних величин для этих параметров не превышают 2 %. Расчетные и критические значения параметров были определены по программе STATISTICA 10.0.

Таблица 1 - Оценка достоверности результатов КНД-М в сопоставлении с контрольным методом (геологическим опробованием)

ГО КНД-М

Показатели сопоставления

m, м Мс, м% С, % m, м Мс, м% С, %

Количество интервалов - 38

Сумма 122,20 11,4089 119,70 11,2493

Среднее значение 3,22 0,3002 0,0934 3,150 0,2960 0,0940

Значение ^статистики 1,55 0,33 0,57

Критическое значение ^ф1-1-1 (при Р = 95 %) 2,02 2,02 2,02

Случайное расхождение по мощности (Бт), м 0,189

Относительное случайное расхождение по метропроценту (Бтс), % отн. 24,8

Случайное расхождение по Си (Бт), % 0,0244

Относительные расхождения средних величин, % -2,0 0,7 -1,4

На следующем этапе был произведен сравнительный анализ результатов ГК для всех технологических и разведочных скважин, по которым был выполнен каротаж КНД-М. Расчетные значения ^статистики Стьюдента гораздо выше критических. Пороговое для 99,9 % доверительной вероятности составляет 3,319. В то же время расчетные значения составляют 8,17 для мощности и 4,12 для метропроцента (табл. 2).

Таблица 2 - Оценка достоверности результатов КНД-М в сопоставлении с контрольным методом (гамма-каротаж)_

Показатели сопоставления ГК КНД-М

m, м Мс, м% С, % m, м Мс, м% С, %

Количество интервалов - 369

Сумма 1053,4 76,97 1162 92,23

Среднее значение 2,51 0,185 0,069 2,77 0,221 0,071

Значение ^статистики 8,17 4,12 1,81

Критическое значение tкp(1-(при Р = 99,9 %) 3,319 3,319 3,319

Случайное расхождение по мощности (Бт), м 0,33

Относительное случайное расхождение по метропроценту (Бтс), % отн. 25 %

Случайное расхождение по Си (Бт), % 0,021

Относительные расхождения средних величин, % -9,34 -16,2 -2,5

Тем самым ^статистика показывает наличие статистически значимых расхождений в мощности РИ и в метропроценте между ГК и КНД-М. Так, по мощности имеет место занижение по результатам ГК на 9,34 %, по метропроценту - на 16,2 %. Сравнивая полученные результаты между ГК и КНД-М в технологических скважинах с результатами, полученными в разведочных скважинах между керновой аналитикой и ГК, мы видим сопоставимые величины существенного занижения метропроцента. Для исключения таких расхождений производилась корректировка базы данных КНД-М и формирование объединенной выборки по всем методам исследования.

С целью получения репрезентативных выборок для статистической оценки средних значений Крр по результатам КНД-М и ГК в проницаемых рудах для мешков и крыльев месторождений было проведено объединение соответствующих выборок по разведочным и технологическим скважинам.

Идея создания данной выборки заключалась в объединении исходной БД с рассчитанным Крр, с данными КНД-М при геолого-разведочных работах и КНД-м при технологическом бурении.

Учитывая, что менее половины скважин, в которых выполнен КНД-М, были подвергнуты опробованию, на данном этапе включим в расчет Крр только опробованные разведочные скважины и технологические скважины с КНД-М. Такое объединение позволит сформировать наиболее полную статистическую базу данных на основе имеющихся на месторождении сведений. Таким образом, из 32 разведочных скважин с КНД-М в статистической базе данных остается лишь 12. К ним добавляется 92 технологических скважины с КНД-М и 167 скважин с керновой аналитикой. Таким образом, в исследуемой статистической выборке участвует 271 скважина с 833 рудными интервалами (321 для мешков и 512 для крыльев месторождений соответственно).

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

Характерная черта образования инфильтрационных месторождений заключается в том, что в зоне непосредственного осаждения рудного вещества происходит непрерывное растворение ранее сформировавшихся руд в тыловой части этой зоны и перенос вещества во фронтальную ее часть. Формируемое оруденение непрерывно движется, как бы «перекатывается» в направлении водного потока, поэтому такие рудные тела называют «роллы». Фронтальную часть рола с максимальной мощностью оруденения называют «мешком» или «меш-ковой частью» рудного тела. Тыловые маломощные части, прилегающие к верхней и нижней границам фильтрующего пласта, называют «крыльями»

Рисунок 3 - Геохимическая зональность в пласте с рудной залежью роллового типа

Таблица 3 - Статистические характеристики полной объединенной выборки для оценки средних значений Крр по данным каротажей КНД-М, ГК и ГО технологических и разведочных скважин

Статистическая характеристика Для мешков Для крыльев

Количество рудных интервалов 321 512

Среднее значение Крр, отн. ед. 0,543 0,784

СКО для Крр, отн. ед. 0,225 0,257

Погрешность среднего Крр, отн. ед. 0,012 0,011

Относительная погрешность определения среднего Крр (для Р = 95 %), % отн. 3,99 2,98

Расчетное значение критерия Смирнова 3,20 3,74

Критическое значение статистики Смирнова для Р = 95 % 3,61 3,75

Расчетные и критические значения параметров были определены с использованием программы STATISTICA 10.0

Таким образом, обе объединенные выборки, представляющие из себя полную базу данных по Крр, полученных различными методами (ГО, ГК и КНД-М), репрезентативны, однородны как по критерию Смирнова, так и по СКО. Среднее значение Крр для выборки мешков составляет 0,543 ± 0,012; для выборки крыльев - 0,784 ± 0,011. Относительная погрешность определения среднего Крр (для Р = 95 %) составляет 3,99 % для мешков и 2,98 % для крыльев. Рекомендованное значение относительной погрешности для достоверно определенной средней оценки Крр составляет 5 % [6]. Полученные характеристики соответствуют основным инструктивным требованиям Инструкции по гамма-каротажу [5] в части точности определения средней оценки Крр и требованиям к выборкам, предназначенным для определения средних оценок Крр.

^- 7

1 - урановая рудная залежь; --

1 - пески; 3 - глнны;

4-окисленные Буровато-желтые породы (зона пластового окисления);

5 - серые породы;

6 - граница зоны пластового окисления;

7 - направление движения пластовых вод; В - скважина.

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

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

Статистическое моделирование мощности рудных интервалов. На следующем этапе было проведено статистическое моделирование мощности рудных интервалов в зависимости от характера морфологических границ на основе определения трендов графиков бортовых кондиций cRaср/cRaборт. Они были построены на основании определения зависимости бортовой массовой доли радия на границах рудных (по урану) интервалов от средней массовой доли радия (рис. 4) [4]. Очевидно, что на границах бортовых кондиций мешковой и крыльевой части залежи зависимости, выражающие отношение среднего радия по борту 0,012 % к среднему радию в рудном интервале, весьма устойчивы, изменяясь от 0,01 сRaср до 0,06 (серая область на рисунке 4). Коэффициент достоверности аппроксимации для регрессионных уравнений достаточно велик и составляет 0,96 для крыльев и 0,99 для мешков соответственно.

Рисунок 4 - Зависимости бортовой массовой доли радия на границах рудных (по урану) интервалов от средней массовой доли радия сКаср/сЯаборт для мешков и крыльев

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

cRaборт = 0,00381п^аср) + 0,0223, R2 = 0,988 (для мешков); (1)

cRaборт = 0,00251п^аср) + 0,0188, R2 = 0,963 (для крыльев). (2)

Данные уравнения взаимосвязи cRaборт = f(cRaср) целесообразно использовать для расчета границ рудных интервалов (рис. 4).

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

На следующем этапе исследования было проведено статистическое моделирование коэффициента радиоактивного равновесия от среднего содержания радия для определения искомого содержания урана в рудном интервале [10] (рис. 5).

1,1

0,9

0,7

0,5

0,3

Крр, отн.ед.

Мешковые разности

Крр = 0,07531п№) + 0,7798

Крыльевые разности

Крр = 0,0961п^а) + 1,0689 —!

32 _

R2 = 0,9422

Аномальная зона богатых руд

0

0,1

0,2

0,3

Среднее интервальное cRa, отн.ед. 0,4 0,5 0,6 0,7

Рисунок 5 - Зависимость Крр от среднего содержания радия (по данным объединенной статистической базы)

Причиной отклонений экстраполяционной кривой от фактических данных послужило наличие на месторождении участков (зон) богатых руд с аномально низкими средними значениями Крр, которые изначально были отнесены к технологическому забалансу, сЯа = 0,51, Крр = 0,67 (для мешков) и сЯа = 0,63, Крр = 1,0 (для крыльев) - серая зона на рисунке 5. Такие рудные интервалы могут являться одним из дополнительных источников извлечения урана из недр. На основании данной гипотезы была проведена оценка прироста запасов за счет вовлечения таких рудных интервалов в процесс скважинного подземного выщелачивания для выборок мешков и крыльев.

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

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

Этап 1. Определение рудных интервалов по кривой соотношения сЯаср/сЯаборт. (выверяется мощность рудных интервалов путем проверки подсчитанного среднего содержания Яа в верхнем и нижнем контакте рудного интервала).

Этап 2. Определение Крр по статистической зависимости отдельно для различных элементов рудной залежи (мешков и крыльев).

Этап 3. Расчет средних значений мощностей рудных интервалов и средневзвешенных содержаний урана в зависимости от полученного значения Крр.

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

Непосредственно моделирование (пересчет уровня радия - определение содержания урана -пересчет Крр - определение мощности рудного интервала и метропроцента - расчет запасов технологического блока) прироста запасов технологических блоков рудника подземного скважинного выщелачивания урана производится на основе полученных зависимостей Крр от сЯа. Для мешков:

Крр = 0,07531п(сЯа) + 0,7798, Я2 = 0,74, (3)

при сЯа < 0,1 %; Крр = 0,7 при сЯа > 0,1 %.

Крр для крыльев принимается в виде зависимости от сЯа по:

Крр = 0,0961п(сЯа) + 1,0689, Я2 = 0,9422, (4)

при сЯа < 0,12 %; Крр = 0,96 при сЯа > 0,12 %.

Ориентируясь на такие ограничения, были смоделированы основные параметры подсчета запасов технологических блоков. Следует отметить, что применение данных зависимостей позволя-

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

Таким образом, в результате моделирования были получены значения прироста (относительно расчетных данных предприятия при фиксированных значениях Крр) вскрытых запасов (3804 т). Это на 15,6 % больше, чем по первоначальным оценкам предприятия (рис. 6). При этом среднее и стандартное отклонение для извлечения урана составляет 78,54 ± 13,3%, а максимальное значение по модели не превышает 91,3 % (рис. 7).

300

250

Прирост запасов, % Прирост запасов, т.

Прирост запасов по технологическим блокам составляет 3804 т (15,6 %) (моделирование)

4000 3500 3000 2500 2000 1500 1000 500 0

Технологические блоки месторождения

Рисунок 6 - Результаты моделирования прироста запасов технологических блоков месторождения

100

90

80

70

60

50

40

30

ИМ, % (модель)

Мах. извл.(модель) = 91,30 %

• * 4« •*•

•V

Ср. извл., (модель)= 78,54 %

+СТ

/

ст=±13,3 %

• •

Блоки на ранней стадии отработки запасов

ВЗ, т (модель)

200

400

600

800

1000

1200

1400

Рисунок 7 - Извлечение металла (ИМ, %) из вскрытых запасов (ВЗ, т) по технологическим блокам месторождения по результатам моделирования (среднее значение извлечения металла (ИМсрзнач = 78,54 %), стандартное отклонение для извлечения металла составляет (ст = ±13,3%), ИМтах = 91,3 %)

а

0

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

(КВз = (1-(ГРМ-ВЗ)/ГРМ))*100 %), (5)

и вскрытыми запасами месторождения по результатам моделирования (рис. 8).

100

80

60

40

20

200

400

600

800

1000

1200

1400

0

0

Рисунок 8 - Ресурсный «потенциал» (Квз - коэффициент отработки вскрытых запасов) месторождения по данным моделирования вскрытых запасов технологических блоков

Было выявлено, что 68 % всех технологических блоков отработаны не более чем на 60 % с вскрытыми запасами, не превышающими 600 т. Это соотношение еще раз подтверждает наличие значительного ресурсного «потенциала» в отработке технологических блоков месторождения в отношении промышленных запасов урана (рис. 8).

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

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

Выводы. 1. Проведен сравнительно-сопоставительный анализ методов оценки запасов технологического блока, таких как керновое геологическое опробование, гамма-каротаж, метод мгновенных нейтронов деления, и определены специфические особенности статистической интерпретации полученных данных. Были выявлены систематические и случайные (по статистическим критериям) расхождения в параметрах рудных интервалов для различных классов руд и морфологических элементов. Установлено, что качество каротажа КНД-М удовлетворительное. Результаты КНД-М обладают достаточно высокой степенью достоверности и могут быть использованы для количественной оценки Крр рудных интервалов.

2. На основе сформированной объединенной и скорректированной статистической базы данных по разведочным и технологическим скважинам были получены регрессионные взаимосвязи коэффициента радиоактивного равновесия от среднего содержания ЯЯа и определены численные значения Крр для мешковых и крыльевых элементов рудной залежи.

3. Проведено статистическое моделирование мощности рудных интервалов в зависимости от характера морфологических границ на основе определения трендов графиков бортовых кондиций cRa^/cRa^ и определена необходимость корректировки рудных мощностей на основании определения зависимости бортовой массовой доли радия на границах рудных (по урану) интервалов от средней массовой доли радия.

4. Разработан подход к интерпретации и пересчету запасов технологических блоков рудника подземного скважинного выщелачивания урана. Его концепция заключается в определении основных причин переизвлечения запасов на месторождении, выявлении статистически значимых расхождений между результатами расчета рудных интервалов по ГК и контрольными методами, определении процента занижения результатов по ГК в зависимости от морфологии и/или от класса руд по содержанию, переинтерпретации морфологии (при необходимости) и расчете новых значений Крр на основе максимально широкого набора данных, что позволит разработать предложения по оптимизации расходов на горно-подготовительные работы и непосредственно добычу.

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

Библиографический список

1. Басов B. C. Выбор оптимальной схемы вскрытия месторождения геотехнологическими скважинами / B. C. Басов // Горный журнал. - 2006. - № 1. - С. 51-53.

2. Верхотуров А. Г. Интенсификация добычи урана при использовании комплекса обработки прифиль-тровых зон геотехнологических скважин / А. Г. Верхотуров, А. А. Сабигатулин // Горный информационно-аналитический бюллетень. - 2019. - № 7. - С. 13-20.

3. Голик В. И. Инновационные технологии добычи урана / В. И. Голик // Горный журнал. - 2009. - № 2. -С. 12-15.

4. Дементьев А. А. Компьютерная модель расчета технико-экономических показателей (на примере месторождений урана для разработки подземным выщелачиванием) / А. А. Дементьев // Горный информационно-аналитический бюллетень. - 1999. - № 2. - С. 169-175.

5. Инструкция по гамма-каротажу на месторождениях урана (МИ-№ 115/210-(01.00250-2008)-2011). -СПб. : ФГУНПП «Геологоразведка», 2012. - 135 с.

6. Инструкция по каротажу методом мгновенных нейтронов деления при изучении урановых месторождений гидрогенного типа. - Ленинград : НПО «Рудгеофизика», 1985. - 63 с.

7. Рогов А. Е. К определению оптимального уровня извлечения урана при подземном скважинном выщелачивании / А. Е. Рогов // Горный информационно-аналитический бюллетень. - 2003. - № 9. - С. 119-121.

8. Рогов Е. И. Оптимизация подготовленных и готовых к выемке запасов на рудниках подземного скважинного выщелачивания урана / Е. И. Рогов, А. Е. Рогов, В. Г. Язиков, В. Л. Забазнов // Горный информационно-аналитический бюллетень. - 2002. - № 4. - С. 149-150.

9. Шаталов В. В. Подземное выщелачивание урана и пути его совершенствования / В. В. Шаталов, М. И. Фазлуллин // Цветные металлы. - 2003. - № 4. - С. 35- 39.

10. Шурыгин С. В. Влияние коэффициента радиоактивного равновесия на сортируемость руд радиометрическими методами / С. В. Шурыгин, В. А. Овсейчук // Горный информационно-аналитический бюллетень. -2016. - № 4. - С. 376-381.

References

1. Basov V. S. Vybor optimalnoy skkhemy vskrytiya mestorozhdeniya geotekhnologicheskimi skvazhinami [Selection optimal scheme field opening by geotechnological wells]. Gornyy zhurnal [Mining Journal], 2006, no. 1, pp. 51-53.

2. Verkhoturov A. G., Sabigatulin A. A. Intensifikatsiya dobychi urana pri ispolzovanii kompleksa obrabotki prifiltrovykh zon geotekhnologicheskikh skvazhin [Intensification of uranium mining using the geotechnological wells filter zones treatment complex]. Gornyy informatsionno-analiticheskiy byulleten [Mining Information and Analytical Bulletin], 2019, no. 7, pp. 13-20.

3. Golik V. I. Innovatsionnye tekhnologii dobychi urana [Innovative uranium mining technologies]. Gornyy zhurnal [Mining Journal], no. 2, pp. 12-15.

4. Dementev A. A. Kompyuternaya model rascheta tekhniko-ekonomicheskikh pokazateley (na primere mes-torozhdeniy urana dlya razrabotki podzemnym vyshchelachivaniem) [Technical and economic indicators calculation computer model (on the example uranium deposits for development by underground leaching)] Gornyy informatsion-no-analiticheskiy byulleten [Mining Information and Analytical Bulletin], 1999, no. 2, pp. 169-175.

5. Instruktsiya po gamma-karotazhu na mestorozhdeniyakh urana (MI-Ns 115/210-(01.00250-2008)-2011) FGUNPP «Geologorazvedka» [Instruction on gamma-ray logging at uranium deposits (MI-№ 115/210-(01.00250-2008)-2011)]. St. Petersburg, FSUNPP "Geological Exploration", 2012. 135 р.

6. Instruktsiya po karotazhu metodom mgnovennykh neytronov deleniya pri izuchenii uranovykh mestorozh-deniy gidrogennogo tipa [Instruction on instantaneous fission neutron logging in the hydrogen uranium deposits study]. Leningrad, NPO «Rudgeofizika» Publ., 1985. 63 p.

7. Rogov A. E. K opredeleniyu optimalnogo urovnya izvlecheniya urana pri podzemnom skvazhinnom vyshchelachivanii [To determine the uranium extraction optimal level in underground well leaching]. Gornyy infor-matsionno-analiticheskiy byulleten [Mining Information and Analytical Bulletin], 2003, no. 9, pp. 119-121.

8. Rogov E. I., Rogov A. E., Zabaznov V. L., Yazikov V. G Optimizatsiya podgotovlennykh i gotovykh k vyemke zapasov na rudnikakh podzemnogo skvazhinnogo vyshchelachivaniya urana [Optimization of prepared and ready-to-excavate reserves at uranium]. Gornyy informatsionno-analiticheskiy byulleten [Mining Information and Analytical Bulletin], no. 4, pp. 149-150.

9. Shatalov V. V., Fazlullin M. I. Podzemnoe vyshchelachivanie urana i puti ego sovershenstvovaniya [Underground leaching of uranium and ways its improvement]. Tsvetnye metally [Non-ferrous Metals], 2003, no. 4, pp. 35-39.

10. Shurygin S. V., Ovseychuk V. A. Vliyanie koeffitsienta radioaktivnogo ravnovesiya na sortiruemost rud ra-diometricheskimi metodami [Effect of the radioactive equilibrium coefficient on the sortiability ores by radiometric methods]. Gornyy informatsionno-analiticheskiy byulleten [Mining Information and Analytical Bulletin], 2016, no. 4, pp. 376-381.

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