CC BY
30
УДК 551.49.018
Д.Л.УСТЮГОВ, канд. геол. -минерал. наук, доцент, ustiugov@mail. ru Национальный минерально-сырьевой университет «Горный», Санкт-Петербург
D.L.USTYUGOV, PhDr in geol. & min. sc, associate professor, [email protected] National Mineral Resources University (Mining University), Saint Petersburg
ПОСТОЯННОДЕЙСТВУЮЩАЯ ГИДРОДИНАМИЧЕСКАЯ МОДЕЛЬ ЯКОВЛЕВСКОГО МЕСТОРОЖДЕНИЯ БОГАТЫХ ЖЕЛЕЗНЫХ РУД (КУРСКАЯ МАГНИТНАЯ АНОМАЛИЯ)
Основное назначение исследований связано с оценкой водопритоков в выработанное пространство Яковлевского железорудного месторождения методом математического моделирования. Использование постояннодействующих математических моделей позволяет оперативно реагировать на быстро меняющуюся гидрогеологическую обстановку, оценивать прогнозные водопритоки в горные выработки в ходе дальнейшего ведения горных работ и решать инженерные задачи по оптимизации дренажных мероприятий.
Ключевые слова: математическое моделирование, водопритоки, Яковлевское железорудное месторождение.
CONSTANTLY WORKING HYDRODINAMICAL MODEL OF YAKOVLEVSKY RICH IRON-ORES DEPOSIT (KURSK MAGNETIC ANOMALY)
Main goal of researches is bound with an estimation of values of inflows of water in goal Yakovlevsky of an iron-ore deposit using a method of mathematical modeling.
Usage of constantly working mathematical models allows reacting operatively to fast varying hydro-geological circumstances and allows to value prognosis inflows of water in mine workings during the further support of mining operations and to solve engineering tasks on optimization of drainage actions.
Key words: mathematical modeling, water inflows, Yakovlevsky iron-ore deposit.
Принятие обоснованных решений по управлению процессом недропользования для геоэкологических систем со сложными природными геолого-гидрогеологическими условиями и интенсивным техногенным воздействием возможно при использовании компьютерно-информационных технологий на базе математических моделей, созданных с применением современных методов схематизации природных гидрогеологических условий. На примере разработанной нами постояннодействующей гидродинамической модели Яковлевского месторождения богатых железных руд рассмотрим основные
принципы создания математической модели и предлагаемые пути решения поставленных задач.
Природная гидрогеологическая система Яковлевского месторождения схематизирована в виде гидравлически связанной системы напорных горизонтов, разделенных слабопроницаемыми слоями, с учетом имеющихся представлений о питании, стоке и разгрузке водоносных горизонтов. При этом за основу была принята гидрогеологическая стратификация соответствующего масштаба для данной территории; учтены водоотлив и степень изученно-
сти водоносных горизонтов по гидрогеологическим параметрам и уровенному режиму.
При гидрогеологической схематизации в плане (обосновании плановых границ модели) принималось во внимание: топология водосборных бассейнов, границы распространения водоносных горизонтов, во-допонизительных систем горных выработок, т.е. площадь и интенсивность техногенного воздействия на исследуемой территории. При гидрогеологической схематизации во времени проводился анализ рядов режимных наблюдений за уровнями подземных вод, ретроспективных данных по водоотливу.
В обводнении Яковлевского железорудного месторождения принимают участие два водоносных горизонта - руднокристал-лический и нижнекаменноугольный, гидравлически взаимосвязанные между собой. В свою очередь, нижнекаменноугольный водоносный горизонт на отдельных участках представляет единую систему с келло-вейским (где мощность водоупорных бат-байоских глин, разделяющих эти горизонты, уменьшается до 10 м, в основном мощность глин бат-байоса составляет 30-50 м), а в местах непосредственного контактирования водоносных известняков нижнего карбона с рудным телом - и с руднокристаллическим горизонтом. Однако наличие плотных переотложенных руд, карбонатизированных бокситовых образований в кровле руднокри-сталлической толщи в пределах рудных полос, а также глинистых отложений в подошве нижнекаменноугольных отложений в лежачем и висячем боках рудного тела затрудняет взаимосвязь нижнекаменноугольного и руднокристаллического водоносных горизонтов.
Учитывая вышеизложенное, при схематизации гидрогеологических условий исследуемого района принята трехслойная плановая модель геофильтрации подземных вод: руднокристаллический водоносный горизонт, разделяющий слабопроницаемый слой и нижнекаменноугольный водоносный горизонт.
Математическое моделирование Яковлевского месторождения осуществлялось с применением программы «Processing
Modflow 5/3» (В.Чанг и В.Кинцельбах), которая реализует численную конечно-разностную модель фильтрации подземных вод для расчета пространственно-временного распределения напоров в трехмерной постановке и базируется на уравнении неразрывности фильтрационного потока.
Ввод исходной информации осуществлялся в форматном виде, который более удобен при моделировании неоднородных водоносных горизонтов.
Для проведения моделирования вся область фильтрации была разбита на 12996 расчетных блоков, по оси абсцисс - 114 блоков, по оси ординат - 114 блоков. Задача решалась на неравномерной сетке, причем в центральной части разбивка более дробная (в районе Яковлевского рудника), размер каждого блока 25x25 м. К периферии области моделирования размеры блоков увеличивались до 500x500 м. Моделируемая площадь в итоге получилась равной 36 км2.
Для построения матриц кровли и подошвы нижнекаменноугольного водоносного горизонта использовалась карта мощности и рельефа подошвы глинисто-известняковой толщи Сь для остальных слоев - инженерно-геологические разрезы по профилям и база данных по скважинам, вскрывшим руднокристаллический водоносный горизонт.
Гидродинамические границы модели и их физико-математическое представление в численной схеме выбирались исходя из анализа общей геолого-гидрогеологической ситуации. Нижняя граница области геофильтрации принята по абсолютной отметке минус 800 м в соответствии с данными о распространении тектонических разломов в массиве руднокристаллических пород до глубины 1,0-1,1 км.
По внешней границе области фильтрации в обоих водоносных горизонтах учитывалась схема «неограниченного пласта» (условия I рода, H = const).
При моделировании производственного понижения дренажная система Яков-левского рудника в нижнекаменноугольном и руднокристаллическом водоносных горизонтах реализована на модели граничными
условиями II рода (Q = const). Водоприток к горным выработкам (околоствольные горные выработки, горные выработки порожнякового и вентиляционного квершлагов, штрека между околоствольными дворами, откаточного штрека лежачего блока, грузового квершлага, горные выработки первоочередного выемочного участка) реализованы на модели через граничные условия III рода с поддержанием уровня подземных вод на абсолютной отметке минус 425 м через специальный расчетный блок «Drain» программы «Processing Modflow 5/3».
Калибровка модели проводилась в процессе решения обратных (эпигнозных) задач с целью уточнения фильтрационных параметров водоносных горизонтов и разделяющего слабопроницаемого слоя. Для этого проводились численные эксперименты при различных граничных условиях, учете взаимосвязи водоносных горизонтов, изменении фильтрационных характеристик водоносных горизонтов и параметрах взаимодействия водоносных горизонтов на различные моменты времени.
Правильность оценки параметров модели определялась: 1) степенью сходимости модельных и фактических отметок уровней подземных вод; 2) соответствием реального и модельного водопритоков в горные выработки. При калибровке моделей использовались среднегодовые уровни подземных вод по скважинам системы гидродинамического мониторинга.
Анализ соответствия между модельными и фактическими положениями уровней подземных вод нижнекаменноугольного и руднокристаллического водоносных горизонтов по опорным точкам (см. таблицу) показывает хорошую сходимость результатов.
Помимо сравнительной оценки по уровням водоносных горизонтов выполнена также корректировка модели по модельным и фактическим водопритокам по состоянию на 2010 г. Фактический водоприток к горным выработкам 10560 м3/сут, модельный 10675 м3/сут; абсолютная погрешность 115 м3/сут; относительная погрешность 1,1 %.
Результатами решения эпигнозных задач явились карты фильтрационных пара-
метров водоносных горизонтов, карты гид-роизопьез, балансовые характеристики (водопритоки, перетекание, горизонтальный и вертикальный водобмен между горизонтами).
Подводя итоги выполненных работ по построению и калибровке модели можно сделать следующие выводы:
• созданная модель не противоречит всей совокупности имеющейся геолого-гидрогеологической информации;
• модель адекватно отражает основные закономерности, принятые при схематизации природных гидрогеологических условий;
• реакция модели функционально (количественно и качественно) соответствует реакции природной гидрогеологической системы на различного рода возмущения.
Следовательно, созданная модель может служить основой для принятия различных прогнозных решений и оценок состояния геологической среды.
Решение прямых (прогнозных) задач проводилось с целью анализа возможного изменения гидродинамической обстановки, выработки рекомендаций и принятия решений по повышению безопасности эксплуатации Яковлевского рудника.
На разработанной гидродинамической модели были реализованы следующие прогнозные задачи:
• определение водопритоков в горные выработки Яковлевского железорудного месторождения при выходе рудника на проектную мощность;
• расчеты перетекания из каменноугольного водоносного горизонта в руднок-ристаллический водоносный горизонт через разделяющий слабопроницаемый слой.
Моделирование водопритоков в горные выработки месторождения осуществлялось с учетом возрастания коэффициента фильтрации руднокристаллического водоносного горизонта и увеличения фронта очистных работ, что существенно определяет безопасные условия ведения выемки полезных ископаемых под неосушенным нижнекаменноугольным водоносным горизонтом.
Анализ выполненных прогнозных расчетов показал, что при увеличении коэффициента фильтрации в 2-4 раза среднегодовой
Уровни подземных вод нижнекаменноугольного и руднокристаллического водоносных горизонтов
(по состоянию на 2010 г.)
Номер скважины Место расположения скважины Абсолютная отметка уровня подземных вод, м Абсолютная погрешность, м
фактическая модельная
32-кр Лежачий бок 107,57 109,8 2,23
7-кр « « 91,46 93,5 2,04
28-кр Рудная полоса 48,50 50,6 2,1
640 Висячий бок 55,58 57,2 1,62
638 « « 50,39 51,8 1,41
610 « « -320,35 -329,8 9,45
642 Рудная полоса -341,99 -347,2 5,21
603 Лежачий бок -150,36 -156,7 6,34
619 Висящий бок -110,74 -115,8 5,06
водоприток из руднокристалического водоносного горизонта к первоочередному участку отработки будет находиться в диапазоне 200-360 м3/ч.
Ведение горных работ под неосушен-ными водоносными горизонтами при условии закладки выработанного пространства быстро твердеющими материалами приводит к уменьшению водопритоков. Если брать в расчет, что коэффициент фильтрации руднокристаллического водоносного горизонта с учетом закладки выработанного пространства будет составлять 0,05 м/сут, водоприток к горным выработкам не превысит 165 м3/ч.
Необходимо отметить, что водопри-токи в горные выработки будут существенно зависеть от параметров разделяющего
слоя и целика руднокристаллических пород. Значения таких водопритоков получены при коэффициенте фильтрации и мощности глинистого слоя в контактной зоне соответственно 3 10-5 м/сут и 14,5 м. Мощность целика принималась по минимальному значению 53 м, коэффициент фильтрации - 0,28 м/сут. Водопритоки к дренажной системе будут формироваться, главным образом, за счет перетекания из нижнекаменноугольного водоносного горизонта. При коэффициенте фильтрации разделяющего слоя 3 10-5 м/сут перетекание будет составлять порядка 305 м3/ч.
Таким образом, водоприток к горным выработкам при выходе рудника на проектную мощность может достигать 470 м3/ч.
