Установление закономерностей фильтрации рассолов для решения проблем загрязнения природных вод на Верхнекамском месторождении калийных солей

Макашов, Сергей Эдуардович

Бесплатный автореферат и диссертация по наукам о земле на тему
Установление закономерностей фильтрации рассолов для решения проблем загрязнения природных вод на Верхнекамском месторождении калийных солей
ВАК РФ 25.00.07, Гидрогеология

Автореферат диссертации по теме "Установление закономерностей фильтрации рассолов для решения проблем загрязнения природных вод на Верхнекамском месторождении калийных солей"

На правах рукописи

005043606 МАКАШОВ Сергей Эдуардович

УСТАНОВЛЕНИЕ ЗАКОНОМЕРНОСТЕЙ

ФИЛЬТРАЦИИ РАССОЛОВ ДЛЯ РЕШЕНИЯ ПРОБЛЕМ ЗАГРЯЗНЕНИЯ ПРИРОДНЫХ ВОД НА ВЕРХНЕКАМСКОМ МЕСТОРОЖДЕНИИ КАЛИЙНЫХ СОЛЕЙ (на примере шахтных полей 1-го Березниковского и 3-го Соликамского калийных рудоуправлений)

Специальность 25.00.07 - Гидрогеология

Автореферат диссертации на соискание ученой степени кандидата геолого-минералогических наук

1 7 МАМ 2012

САНКТ-ПЕТЕРБУРГ 2012

005043606

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

Научный руководитель —

член-корреспондент РАН, доктор геолого-минералогических наук, профессор

доктор геолого-минералогических наук, профессор, Санкт-Петербургский государственный горный университет, главный научный сотрудник Научного центра геомеханики и проблем горного производства

кандидат геолого-минералогических наук, ОАО «ВНИИГ им. Б.Е. Веденеева», старший научный сотрудник отдела «Основания, грунтовые и подземные сооружения»

Ведущая организация — Федеральное государственное бюджетное учреждение науки Горный институт Уральского отделения Российской академии наук.

Защита состоится 30 мая 2012 г. в 18 ч на заседании диссертационного совета Д 212.224.01 при Санкт-Петербургском государственном горном университете по адресу: 199106 Санкт-Петербург, 21-я линия, д.2, ауд.4312.

С диссертацией можно ознакомиться в библиотеке Санкт-Петербургского государственного горного университета.

Автореферат разослан 28 апреля 2012 г.

Румынии Вячеслав Геииевич

Официальные оппоненты:

Норватов Юлий Александрович

Котлов Олег Николаевич

диссертационного совета канд. геол.-минерал. наук

УЧЕНЫЙ СЕКРЕТАРЬ

КИРЬЯКОВА И.Г.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы.

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

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

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

Цель работы состоит в изучении и математическом описании физических процессов, сопровождающих фильтрацию рассолов в водоносных горизонтах, а также в прогнозировании загрязнения природных вод на участках расположения шламохранилищ и рассолосборников (шахтные поля БКПРУ-1 и СКРУ-3) методами численного моделирования.

Основные задачи исследования.

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

2. Исследование механизмов, ответственных за формирование естественной плотностной стратификации природных подземных вод уфимского водоносного комплекса применительно к территории шахтного поля БКПРУ-1.

3. Оценка изменения гидрогеохимической обстановки шахтного поля БКПРУ-1 в период аварии и после затопления горных выработок.

4. Анализ влияния рельефа кровли водоупорной соляно-мергельной толщи на развитие ореолов техногенных рассолов в надсолевой части разреза на примере района шахтного поля СКРУ-3.

Основные методы исследований.

Настоящие исследования основывались:

- на решении задач динамики подземных вод, описывающих потоки жидкостей переменной плотности аналитическими методами и методами численного (математического) моделирования;

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

Для решения теоретических и практических задач применялись современные программные комплексы SUTRA (Voss and Provost, 2002), SEA WAT (Guo and Langevin, 2002), TOUGH2 (Pruess, 1991) и DENSFLOW (П.К. Коносавский, 2002).

Научная новизна.

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

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

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

Изучены главные факторы, контролирующие формирование техногенного ореола засоления в районе СКРУ-3. Показана значимость влияния рельефа водоупорной поверхности на динамику развития ореола рассолов в верхних водоносных горизонтах.

Фактический материал и личный вклад автора.

В основу диссертационной работы положен обширный фондовый материал по геолого-гидрогеологическим условиям Верхнекамского месторождения калийных солей. Автором выполнен анализ данных режимных наблюдений за состоянием подземных и поверхностных вод за 2002-2010 гг. в пределах шахтных полей Соликамско-Березниковского промышленного узла. Автор принимал участие в полевых работах (отборы проб техногенных рассолов, откачки, гидрологические работы) в районе исследуемых объектов. Автором выполнена большая часть численных расчетов в ходе тестирования программ и поиска новых аппроксимирующих зависимостей, а также разработаны миграционные модели шахтных полей БКПРУ-1 и СКРУ-3.

Основные положения, выносимые на защиту.

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

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

3. Аварийное затопление рудника БКПРУ-1 вызвало изменение конфигурации существующих ореолов техногенных рассолов в надсо-левой толще пород и способствовало снижению выноса хлоридных натриевых солей в поверхностные водотоки, что подтверждается режимными наблюдениями и численным моделированием.

4. Формирование квазистационарного ореола растекания техногенных рассолов в районе шахтного поля СКРУ-3 определяется равнозначным влиянием естественного фильтрационного потока и уклона поверхности водоупорной толщи (кровли соляной залежи).

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

Практическая значимость работы.

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

Апробация работы. Результаты исследований докладывались на Первой Всероссийской конференции молодых ученых, посвященной памяти В.А. Мироненко «Дерзость Надежд» (СПб, 2010), научной конференции, посвященной 75-летию кафедры гидрогеологии СПбГУ «Комплексные проблемы гидрогеологии» (СПб, 2011), VlI-й ежегодной Общероссийской научно-практической конференции «Перспективы развития инженерных изысканий в строительстве в Российской Федерации» (Москва, 2011).

Публикации.

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

Структура и объем работы. Диссертация объемом 160 страниц, включает введение, 5 глав, заключение, список источников из 113 наименований, содержит 84 рисунка и 17 таблиц.

Благодарности. Автор выражает искреннюю признательность своему научному руководителю д.г.-м.н., члену-корр. РАН В.Г. Румынину за помощь и поддержку на всех этапах работы. За конструктивные консультации и помощь автор благодарит: к.г.-м.н. доц. кафедры гидрогеологии СПбГУ A.A. Потапова, к.г.-м.н., зав. кафедрой П.К. Коносавского и к.г.-м.н., в.н.с. СПбО ИГЭ РАН Л.Н. Синдаловского. Особую признательность автор выражает д.г.-м.н. проф. кафедры гидрогеологии и инженерной геологии СПГГУ Р.Э. Дашко за критический анализ работы и ряд предложенных формулировок отдельных ее положений. За моральную поддержку и ценные советы автор искренне благодарит к.г.-м.н. Е.В. Мольского и к.г.-м.н. М.Н. Кочневу. Благодарность за помощь в полевых работах и получении фактических данных по объекту автор выражает В.П. Копшталеву, Л.И. Лесничему и А.О. Воробьевой.

ЗАЩИЩАЕМЫЕ ПОЛОЖЕНИЯ И ИХ ОБОСНОВАНИЕ

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

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

Обоснование данного защищаемого положения потребовало рассмотрения двух классических в гидрогеологии задач.

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

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

Ахр = хг ~х\ » £,л/г , (1)

где £,=1.5^2, х = х/т - координата пересечения контакта с кровлей (подошвой) пласта, r = kApt/mn (к - коэффициент фильтрации, м/сут; Ар = (р, -р2)/рг -плотностной градиент (р,,р2 - плотность рассола и пресной воды, кг/м3); т - мощность пласта, м; п - пористость, t - текущее время, сут). Решение (1) приводится в работах В.М. Шестакова (1960), а также H.H. Веригина и B.C. Саркисяна (1979).

Для исследования влияния анизотропии коэффициента фильтрации пласта на форму и положение контакта жидкостей была создана двухмерная (профильная) модель на базе программы SEA WAT, При анализе результатов моделирования кривые АхР = /(г) на соответствующем

Ахр

( JZ\

1 - ехр

1 * У

(2)

учитывающей анизотропию фильтрационных свойств пласта; здесь £=1.89, $3(|с) = 0.81^к,(к = 4kx,kz -фактор анизотропии, >1).

Полученная формула отвечает известным асимптотикам. Так, при больших т и к = 1 формула дает результаты, близкие к решению (1), а при г < 1 - 5 безразмерная функция Дхр линейно зависит от времени,

что совпадает с ранее выполненными оценками Гарднера (G.F.H. Gardner) с соавторами (1962) и Хеллштрёма (G.Hellstrom) с соавторами (1988). Формула справедлива для случаев, когда к не превышает 10 и обладает достаточной точностью (не хуже 10%) для выполнения экспертных оценок. Данную аппроксимационную зависимость применительно к условиям Верхнекамского месторождения рекомендуется использовать при предварительных прогнозных оценках положения поверхности тяжелых жидкостей, например, при закачке рассолов в нижние горизонты надсолевой толщи, породы которой характеризуются анизотропией фильтрационных свойств, или для вычисления положения фронта соленых вод при их фильтрации от линейного источника (бассейна).

АХ„

• Модель

_ Аналитическая

аппроксимация

12-

X°=vt/n а //////////////////У////,

Рассол <? \ О Вода к .X \ ^

///////,

Рис. !. Концептуальная схема поршневого вытеснения (а) и сравнение результатов численного моделирования (кружки) с расчетами по формуле (2) (сплошные линии) (б). Прерывистая линия - приближенное решение (1)

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

соли (С ) в несовершенной скважине. В предельно упрощенном виде (точечный источник) это решение имеет вид:

С = 1-

Л1/3

v' У

(3)

где tQ = 2пкгпЪъ I3Q - время начала подтягивания соленых вод, Q - дебит скважины, м3/сут; п - пористость, b - начальная мощность пресных вод, м; к: = ЛJkr / kz - фактор анизотропии, / - текущее время, сут.

Однако решение задачи в замкнутой форме для функции концентрации в откачивающей скважине при наличии плотностного градиента неизвестно. Поиск такой зависимости проводился на основе моделирования с использованием кодов TOUGH2 и SUTRA (см. далее).

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

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

с 2 Ля;3

где X - 1 -

1/7 У

С = 1-

К ' У

(4)

1 +

4у

= ехр(л/ю), со = ()!0",

2ппЬ

пкЬ1 Лр/3 - приведенный дебит скважины, Г =

- Уравнение (4), др > о б 'тоивнг

1000

Рис 2 Концептуальная схема субвертикального подтягивания подошвенной соленой воды к несовершенной скважине (о) и пример сравнения численного (точки) и аналитического (сплошные линии, формула 4) решений одного^из вариантов тестовой задачи (цифры на графике - плотность соленой воды, кг/м ); штриховая линия - решение (3) при Ар = 0

Решение (4) справедливо при QIQ' >2(2*- предельный дебит во-доотбора, при превышении которого происходит «прорыв» соленых вод к скважине). Графики на рис. 26 указывают на удовлетворительную

точность зависимости (4).

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

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

В качестве тестов для наиболее широко освещенных в научной литературе программных комплексов (SUTRA, SEA WAT, TOUGH2 и DENSFLOW) были выбраны: 1) «задача Генри», которая описывает положение контакта пресной и соленой воды при внедрении морской интрузии в прибрежный пласт; 2) задача о подтягивания контура плотных вод, залегающих в основании водоносного горизонта, к несовершенному дренажу (линейный сток) в кровле пласта (аналитическое решение Я. Бэра и Г. Дагана); 3) физический эксперимент, выполненный Р. Шинкариолом и Ф. Шварцем с целью изучения нестабильности контакта разновесомых жидкостей при формировании двумерных профильных ореолов рассеяния.

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

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

следования процессов геофильтрации рассолов в районе Верхнекамско-

-50

а { 1=3000 сут 1

б А Г« 15000 су г

.Л № ШешШй. .----

-50

ИЗООО сут

200

400

600 800 1000 м

1000 м

" 0 0.2 0.4 0.6 0.8 1

Рис 3 Влияние дробности сеточной разбивки модели на степень структурирования потока рассола под дном бассейна ( С - относительная концентрация). Коли-честно элементов модели по горизонтали и вертикали: а, б- 110x20; в, г - 440x200

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

Верхнекамское месторождения калийных солей расположено в центральной части Соликамской впадины Предуральского краевого прогиба С точки зрения рассматриваемых проблем миграции в подземных водах природных и техногенных (фильтрационных) рассолов, основной интерес представляет верхняя часть геологического разреза - надсоле-вой комплекс пород (рис. П-2): верхняя подтолща соляно-мергельная толщи - СТМ2 (Р,и 81Д сложенная мергелями и глинами с прослоями гипса, терригенно-карбонатная толща - ТКТ (Р,и з12), состоящая из известняков, мергелей и песчаников, пестроцветная толща -ПЦТ (Р,и ¡й), сложенная красно- и зеленоцветными песками, и четвертичные отложения (О). А -

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

Для оценки последствий аварии на руднике БКПРУ-1, расположенного под г. Березники (рис. П-1) на глубине 250-300 м, в результате которой были затоплены выработки объемом порядка 80 млн. м , необходимо было сначала воссоздать на численной модели исходную (до начала техногенного засоления) гидрогеохимическую зональность (по минерализации). Она выражается в наличии в разрезе надсолевого комплекса пород природных хлоридных натриевых соленых вод и рассолов, содержание солей в которых увеличивается с глубиной (рис. 4, рис. П-3). Высокоминерализованные воды приурочены, в основном, к нижнесоликамскому водоносному горизонту, и, предположительно, образуются вследствие растворения №С1 с кровли соляной залежи.

Минерализация, г/л

О 50 100 150 200 250 300 350

a cm iras / / СКВ. 84/4 ,n„ —~ См) 2ПИ -fea, J98t

—~—~~—~—/ fm l<j» """v*---------- ; \ea. ifr«

-j—-——У

Рис. 4. Изменение минерализации подземных вод с глубиной в пределах Верхнекамского месторождения (точки соответствуют интервалам опробования толщи)

Поскольку кинетика растворения солевых минералов (галита, карналлита и сильвина) носит чисто диффузионный характер, т.е. в данных условиях определяется гидродинамическими характеристиками потока, рассеивающими свойствами среды и интенсивностью плотностной конвекции, для имитации растворения каменной соли возможно применение программ, использующих конвективно-дисперсионный подход. Расчеты производились на трехмерной модели, созданной на базе программы SEA WAT (рис. П-4) вплоть до стабилизации трехмерного поля минерализации, которое отвечало бы естественному. Данная задача решалась в стационарной фильтрационной постановке для ненарушенного режима, т.е. при отсутствии каких-либо техногенных объектов. В качестве единственного природного источника солей в подземных (ин-фильтрогенных) водах рассматривалось растворение галита с «соляного зеркала» - первого от поверхности выдержанного слоя каменной соли.

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

Табл. 1. Параметры модели принятые после ее калибровки

Слои Горизонт к. м/сут К, м/сут 4 м 5„ м <5„ м п м3/м3

1,2 О 0.1-10 0.1-1 1.5 0.15 0.015 0.15

3-10 пцт+ткт.™ 0.1-10 ю-'-о.з 1.2 0.12 0.012 0.1

11-18 ТКТ нижн часть 0.3-15 4.5-10"4-0.3 1 0.1 0.01 0.06

19-26 смт, 0.3-10 1.5-10"5 —3 • 10"3 1 0.1 0.01 0.02

27 смт, 0.05 10"6 1 0.1 0.01 0.01

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

Третье защищаемое положение. Аварийное затопление рудника БКПРУ-1 вызвало изменение конфигурации существующих ореолов техногенных рассолов в надсолевой толще пород и способствовало снижению выноса хлоридных натриевых солей в поверхностные водотоки, что подтверждается режимными наблюдениями и численным моделированием.

В октябре 2006 г. на руднике БКПРУ-1 произошла авария, связанная с нарушением сплошности водозащитной толщи (ВЗТ) и прорывом относительно пресных подземных вод со средним расходом около

4500 м3/ч из надсолевой части разреза в подземные горные выработки. В июле 2007 г. в зоне прорыва образовался провал земной поверхности глубиной до 80 м, обусловленный растворением пород ВЗТ. Провал постепенно расширялся вплоть до момента полного затопления рудника в октябре 2008 г., когда его размеры по поверхности достигли 300x400 м (рис. П-7). К настоящему времени (2011-2012 гг.) практически завершилось полное восстановление уровней подземных вод во всей надсолевой толще. По официальным данным, причиной аварии послужила «совокупность геологических и технологических факторов».

Оценка последствий аварии с точки зрения ее влияния на существующую гидрохимическую обстановку в районе г. Березники в данном случае является необходимой по следующим причинам. Во-первых, образовавшийся провал создал гидравлическую связь как между затопленными выработками и надсолевой толщей, так и между ранее относительно изолированными друг от друга водоносными горизонтами, содержащими воды различной солености. Во-вторых, даже после затопления рудника будет продолжаться деградация (растворение) межкамерных целиков. По оценкам специалистов ГИУрОРАН вследствие сокращения объемов горных выработок, связанного с осадками поверхности, произойдет вытеснение части находящихся в них крепких рассолов с минерализацией более 320 г/л через зону провала в надсолевую толщу. Наконец, мощное гидродинамическое возмущение повлияло на пути миграции существующего ореола загрязнения подземных вод.

На данном этапе производилась модельная реконструкция истории техногенного засоления водоносных горизонтов (до аварии), которая производилась на фоне ранее рассчитанной гидрохимической зональности. Загрязнение формируется в результате фильтрационных потерь рассолов из поверхностных емкостей-накопителей БКПРУ-1, которые эксплуатируются с 1964 г. В модель были добавлены источники загрязнения - шламохранилище, рассолосборник и солеотвал. Негативное влияние указанных объектов подтверждается данными наблюдений по скважинам, расположенным в районе шахтного поля и оборудованным на четвертичный водоносный горизонт. В период до аварии по ним фиксировались повышенные значения минерализации подземных вод (до 287 г/л). Утечки рассолов также сказываются на химическом составе поверхностных вод. Так, минерализации воды в р. Зырянка на участке в пределах шахтного поля не редко превышала 2.5 г/л.

Общая интенсивность утечек в соответствии с ранее выполненными

оценками специалистов ВНИИ Галургии и СПб пред-ва ГИУрОРАН составила порядка 2000 м3/сут. Минерализация техногенных рассолов принималась при моделировании равной 300 г/л. Докалибровка модели заключалась в уточнении фильтрационных параметров только в зоне, прилегающей непосредственно к накопителям.

По результатам расчетов было выявлено, что до аварии 30% от объемов утечек крепких рассолов попадали в поверхностные воды, основная же их часть (около 70%) под действием гравитационных сил аккумулировалась в нижней части надсолевой толщи, формируя обширный ореол загрязнения: солевые растворы мигрировали в западном направлении по поверхности соляной залежи («соляного зеркала»).

Следующий этап - прослеживание изменения конфигурации ореола засоления в период оттока подземных вод из надсолевой толщи в выработанное пространство рудника в ходе его аварийного затопления в 2006-2008 гг. Результаты расчетов показали, что ранее сформировавшийся (в ходе 43-летней нисходящей фильтрации рассолов) ореол оказался в зоне гидродинамического возмущения. Это привело к тому, что тяжелое ядро ореола сместилось в юго-восточном направлении - в сторону провала (рис. П-8), что подтверждается падением минерализации подземных вод в скважинах, расположенных в западной части шахтного поля рядом с источниками загрязнения (рис. П-9, П-10). Вследствие этого произошло снижение интенсивности разгрузки рассолов в поверхностные водотоки: отток рассолов в зону провала почти в шесть раз (с 30% до 5% от суммарных потерь) сократил вынос солей в поверхностные водотоки. При этом в ходе натурных наблюдений в период затопления рудника было зафиксировано снижение прироста концентрации калия в р. Зырянка между гидропостами г/п-1 и г/п-2 (рис. 4).

;1 авария | (10.200) 5Г ......

_____ "1 ш

и 1

2005 2006 2007 2008 2009 2010 гад ■ Фоновый гидропост (г/п-1) О Замыкающий гидропост (г/п-2)

Рис. 4. Изменение концентрации К+ в р. Зырянке в период затопления рудника. Расположение гидропостов см. рис. П-10

Рис. П-1. Ситуационный план района БКПРУ-1. Скважины, оборудованные на отложения: 1 - ПЦТ, 2 - ТКТ, 3 - СМТ2. 4 - шламохранилище и рассолосборник, 5 - солеотвал, 6 - граница выклинивания ПЦТ, 7 - граница шахтного поля, 8 - линия разреза

1« 100 60 20 -20 -60 ■ 100 -140 -180 -220 •260 -зоо

-340 -.180 -420 -460

Скв. I»

Ааис-гое вогьжраинлище

Рис П-З ический гидрогеохимический разрез надсолевой толщи по линии 1-1 - Скважипаь подземных вод (ПВ), интервал опробования (в числителе - кт, м /сут, в знаменах, г/л), 2 - уровень ПВ в отложениях ПЦТ; 3 - уровень ПВ в отложениях ТКТ- 4 Ь ПВ в отложениях СМТ2. Химический тип ПВ и их минерализация: 5карбонатные (М- 1-10 г/л); 6 - сульфатные (М- 10-400 г/л); -даные (М - 100-300 г/л). 8 - Региональный водоупор (СМТ,)

в.

ТЙО'ч Ч

+ + у

ГПП» ГШ: Щ, ¡§¡¡4 [5

Б- И« Н»

140 100 60 20 -20 -60 -100 -140 -180 -220 -260 -ЗОО

-340 -380 -420 -460 -500

р/\ рода

Рис. П-2. Субширотный геологический разрез нижнепермских отложений по линии I-1 - пески супеси, суглинки ((}), 2 - песчаники, аргиллиты, алевролиты (ПЦТ), 3 - известняки, мергели, песчаники, алевролиты (ТКТ), 4 -мергели, глины, гипс (СМТ2), 5 - мергели, каменная соль (СМТ,), 6 -покровная каменная соль (ПКС), 7 - карналлит, каменная соль, сильвинит (СКЗ), 8 - сильвинит, каменная соль (СЗ), 9 - подстилающая каменная соль (ПДКС)

[-4. Структура и граничные условия численной миграционной модели

ГУ I рода

ГУ III рода

Рис. П-5. Сра модельных и фак значений минерав скважинах моделиров; гидрохимич зональное Расположение ск<. рис. П-

+110 м

инфильтрация атмосферных осадков III I_I 1|| | Б

•190 м

С, г/л

запад

источнимоли ("соляное зеркало") - СМТ1

500 1000 1500 м ВОСТОК

300 250 200 150 100

50 J

000

++ ^ ^ ^

% % фР %0\г ^ Скважины -Минерализация фактическая -»-Минерализация на модели

восток

Рис. 11-6. Результаты эпигнозного моделирования формирования естественной вертикальной гидрогеохимической зональности уфимского водоносного комплекса. Расчетное время 10000 лет

Рис. П-7. Фото провала

;Ф—ШВ{—-/-j 8,1*м

Рис. П-8. Влшрийного затопления рудника БКПРУ-1 на ореоления (по данным моделирования). - шламсище, рассолосборник, 2 - солеотвал

от.юження скважины >ш отложения СМТ,

+ 110 М|

граница фиксируемой части ореола

• скважины на отложения ТКТ

Рис. П-10. Схема расположения наблюдательных скважин и

гидропостов БКПРУ-1 (г/п-1, г/п-2 - гидропосты)

провал (зона прорыва ВЗТ)

Р Изменение

низации в с,1х в период 1рии. Расние скважин -1С. П-10

-СКВ. С-8 (Q) -СКВ. 163-г(ТКТ)

-СКВ. С-10 (Q) -СКВ. 1ПИ (CMT2)

-скв. C-lla (Q) СКВ. 2ПИ (CMT2)

-195 м

ЩшН-

Рассолы, вытесняемые из рудника

Природные С1-Иа рассолы выщелачивания

Ореол техногенных Cl-Na рассолов

запад

Рис. П-11. Прогнозный модельный гидрогеохимический профиль кадсолевой толщи на 2020 г.

0 500 1000 1500 м

восток

На заключительном этапе (решение прогнозной задачи) прослеживалось изменение конфигурации ореола засоления подземных вод после затопления рудника, вплоть до 2020 года. В данном случае было учтено, что даже после полного затопления рудника в нем будет продолжаться деформация и частичная деградация межкамерных целиков, что обусловит дополнительные осадки земной поверхности и, соответственно, сокращение объема выработанного пространства. В результате этого будет происходить вытеснение части находящихся в выработанном пространстве крепких рассолов через зону провала. Поэтому вместо откачки из надсолевой толщи в зоне провала в было задано дополнительное «нагнетание» рассолов из рудника с постоянным расходом 1.4 тыс. м3/сут (в соответствии с оценками специалистов ГИУрОРАН) и минерализацией 320 г/л. Результаты моделирования показали, что подавляющая часть крепких рассолов, вытесняемых из рудника, будет идти на формирование локальной линзы в нижней и средней частях СМТ2 (рис. П-11), при этом разгрузка рудничных рассолов в поверхностные воды полностью исключается. К 2020 г. вынос солей из шла-мохранилища и рассолосборника в местную гидрографическую сеть практически восстановится и его объем достигнет прежнего уровня, наблюдавшегося до аварии (более 25%). Следовательно, загрязнение подземных вод зоны активного водообмена и поверхностных вод, как и ранее, будет происходить только за счет утечек техногенных рассолов из емкостей-накопителей. Интенсификации процессов засоления природных вод после аварии не произойдет. Данный вывод уже подтверждается натурными наблюдениями, выполнявшимися в 2009-2011 гг.

Четвертое защищаемое положение. Формирование квазистационарного ореола растекания техногенных рассолов в районе шахтного поля СКРУ-3 определяется равнозначным влиянием естественного фильтрационного потока и уклона поверхности водоупорной толщи (кровли соляной залежи).

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

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

Общий расход фильтрационных потерь техногенных рассолов из шламохралилища СКРУ-3, определенный водно-балансовым методом и варьирующий от 1270м3/сут (2006 г.) до 5650 м3/сут (1996 г.), можно разделить на две составляющие. Первая (~70%) отвечает рассолам, которые после инфильтрации быстро погружаются в нижнюю часть водоносного комплекса и продолжают миграцию по кровле водоупорной (соляной) толщи, не разгружаясь в поверхностные водотоки. Пятно загрязнения, сформированное этими крепкими рассолами можно назвать ореолом растекания или «интрузией».

Ко второй составляющей относятся фильтрационные потери, оказывающие уже в настоящее время непосредственное воздействие на поверхностные и подземные воды зоны активного водообмена (-30% от общих утечек). Загрязнение образуется за счет частичного выноса подземными водами фильтрующихся рассолов (в разбавленном виде), а также их «смыва» с ореола растекания. При этом формируется ореол рассеяния (рис. П-12). Минерализация подземных вод в нем - 20-30 г/л.

Конфигурация и количественные параметры этого ореола, были детально анализированы еще до настоящих исследований (Мольский, Кочнева, Липницкий и др., 1990-2007 гг.). Поэтому, для нас особый интерес представляет процесс развития именно «интрузии» рассолов, динамику которой весьма сложно оценить по натурным наблюдениям.

Эпигнозные и прогнозные миграционные расчеты на плановой модели (рис. П-14), созданной на программе ОЕЫ8РЬО\У, показали, что силы, разнонаправленно действующие на ореол растекания рассолов, приведут его в итоге в статичное положение (рис. П-15, рис. П-16). «Интрузия» окажется «подвешенной» на северо-восточном склоне Тверитин-ского прогиба (на кровле водоупорной толщи - СМТ). В связи с этим была выполнена качественная оценка степени влияния гидродинамических и гравитационных сил на ядро ореола. На модели «проигрывалась» гипотетическая история формирования ореола растекания сначала при отсутствии фильтрационного потока (1/= 0), а затем при отсутствии уклона кровли водоупора (Д = 0). В первом случае воздействие на верхние водоносные горизонты и поверхностные воды оказывается минимальным, поскольку рассолы движутся в противоположном от реки

Щ£0

Ф'"

\ 7

Рис. П-12. Схематическая карта гидроизогипс верхних водоносных горизонтов (шахтное поле СКРУ-3, северное обрамление). 1 - Скважина.

Площади распространения водоносных горизонтов (комплексов):

2 - верхнесоликамского водоносного (в отложениях карбонатных фаций верхних-средних слоев ТКТ); 3 - верхне-нижнесоликамский водоносный горизонт (в сульфат-содержащих отложениях ТКТ и верхних-средних слоях СМТ2). Гидроизогипсы:

4 - водоносного горизонта в отложениях карбонатных фаций верхних-средних слоев ТКТ; 5 - водоносного горизонта в сульфат-содержащих отложениях ТКТ - СМТ2.

6 - линия геолого-гидрогеологических разреза, 7 - фиксируемый ореол рассеяния

Рис. П-13. Схематический геолого-гидрогеологический разрез по линии 1-1 (рис. П-12). 1 - скважина, ее номер и уровень подземных вод, 2 - уровень подземных вод в отложениях карбонатных фаций верхних-средних слоев ТКТ, 3 - уровень подземных вод в сульфат-содержащих отложениях ТКТ - СМТ2

200 150 100 50 0

-50 -100 -150 -200 -250 -300

1МОхранили1

соле-отвап

р. Усолка (ГУ III рода)

" | Техногенное озеро : :(затопленный карьер) ГУ I рода

1ламохранилище (ГУ II рода)

рода

Количество блоков: по горизонтали - 72 по вертикали - 82

100000

Количество слоев -1 (горизонты в ТКТ) Нижняя граница - кровля СМТ

О 1 2 3 КМ

Рис. П-14. Граничные условия и распределение коэффициентов фильтрации численной модели шахтного поля СКРУ-3

Водоносный комплекс "Интрузия" рассолов ни состоянию

в отложениях ТКТ на «>'«■'« 2007 <« "

Водоупорная толша - СМТ2 [ ] Прогнозная граница "интрузии" рассолов РИС.

I на конец 2025 г. (г = 40 лет»

■1 Уровень подземных вод водоносного комплекса ТКТ

Рис. П-16. Схематический разрез «интрузии» рассолов по линии А-Б

О 2 4 6 8 10 12 14 1 6 11 16 21 26 31 36 41

Рис. П-15. Результаты эпигнозных (а) и прогнозных (б) миграционных расчетов, светлые линии - гидроизогипсы, а.о., м; темные линии - изолинии кровли водоупора (СМТ2), а.о., м

А О. м

А 136 Шламохран или те Б

200

-100

0 1000м 2000м

Мощность интрузии, м

р. Усолка

Мощность интрузии, м 1 5 9 13 17 21 25

Мощность интрузии, м

1 5 9 13 17 21 25 29 33 37 41

П-17. Результаты исследований взаимовлияния гидродинамически и плотностного градиентов. а - при отсутствии потока пресных подземных вод; б - при отсутствии уклона кровли водоупора

направлении (рис. Т\-\1а) и заполняют наиболее глубокую часть депрессии, что затрудняет их смыв с поверхности ореола и последующую разгрузку. Во втором случае рассолы движутся строго с фильтрационным потоком в направлении местной дрены (рис. П-176). На модели отмечается контакт «интрузии» с руслом р. Усолка, что приводит к существенному увеличению (в три раза) минерализации речной воды.

Поскольку оба сценария в действительности не фиксируются, можно утверждать, что гидродинамические и гравитационные силы принимают соизмеримое участие в динамике развития загрязнения подземных вод шахтного поля СКРУ-3. Недоучет одного из этих факторов при прогнозе геомиграции тяжелых вод в схожих гидрогеологических условиях может привести к принципиально неверным результатам.

Заключение

Основные результаты исследований заключаются в следующем:

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

2. Определены области применения программ, предназначенных для моделирования потоков переменной плотности в пористой среде, при решении различных гидрогеологических задач. При анализе плановых потоков (площадных ореолов загрязнения) это могут быть компьютерные коды, основанные на предпосылке Гибена-Герцберга, в других случаях требуется более строгое решение трехмерных уравнений конвективно-дисперсионного переноса для потоков переменной плотности.

3. Доказано, что естественная профильная гидрохимическая зональность Верхнекамского месторождения калийных солей формировалась в результате растворения каменной соли с поверхности соляной залежи («соляного зеркала») относительно пресными инфильтрационными водами. Длительность ее формирования составляет порядка 10000 лет.

4. По результатам изучения гидрогеологических последствий аварии на руднике БКПРУ-1 был сделан вывод о том, что авария оказала временное реабилитирующее влияние на состояние подземных и поверхностных вод, что выразилось в снижении степени их засоленности. В дальнейшем (до 2020 г.), интенсивность выноса солей в поверхностные воды вернется к доаварийным показателям.

5. На примере шахтного поля СКРУ-3 установлена возможность формирования в реальных условиях квазистационарных ореолов загрязнения подземных вод тяжелыми жидкостями. Экспериментальные модельные исследования чувствительности процесса миграции рассолов доказывают равнозначную степень влияния гидродинамических и гравитационных сил на формирование загрязнения в данном районе.

Публикации по теме диссертации

1. Коносавский П.К., Потапов A.A., Макашов С.Э. Прогнозная оценка растворения карналлита в целиках горных выработок рудника БКПРУ-1 после его аварийного затопления (Верхнекамское месторождение калийных солей) И Сергеевские чтения. Моделирование при решении геоэкологических задач. Вып. 11. - М.: ГЕОС, 2009. - С. 408411.

2. Румынии В.Г., Синдаловский Л.Н., Макашов С.Э., Воронина A.B. Новые аналитические зависимости для описания подтягивания глубинных рассолов к контуру несовершенной скважины // Вестник СПбГУ, Сер. 4, Геология, география. - 2010, №2 - С. 17-37.

3. Макашов С.Э. Исследование миграции техногенных рассолов методами численного моделирования в районе СКРУ-3 Верхнекамского месторождения калийно-магниевых солей // Первая Всероссийская конференция молодых ученых, посвященная памяти В.А. Мироненко. -СПб.: Изд-во СПбГУ, 2010. - С. 77-88.

4. Потапов A.A., Коносавский П.К., Макашов С.Э. Прогноз изменения гидрогеоэкологической обстановки в районе рудника БКПРУ-1 (Верхнекамское месторождение калийно-магниевых солей) после его аварийного затопления // Инженерные изыскания в строительстве. Материалы седьмой общероссийской конференции изыскательских организаций. М.: 2011. С. 192-194.

5. Потапов A.A., Макашов С.Э., Коносавский П.К. Прогноз изменения гидрогеохимической обстановки в районе рудника БКПРУ-1 (Верхнекамское месторождение калийно-магниевых солей) после его аварийного затопления // Комплексные проблемы гидрогеологии. Тезисы докладов научной конференции. - СПб.: Изд-во СПбГУ, 2011. - С. 145148.

6. Макашов С.Э., Потапов A.A. Численное моделирование влияния аварийного затопления соляного рудника БКПРУ-1 на миграцию рассолов в надсолевом комплексе пород // Естественные и технические науки. Выпуск 2. -М.: Спутник плюс, 2012. С. 208-212.

РИ1Д СПГГУ. 19.04.2012. 3.267 Т. 100 экз. 199106 Санкт-Петербург, 21-я линия, д.2

Текст научной работыДиссертация по наукам о земле, кандидата геолого-минералогических наук, Макашов, Сергей Эдуардович, Санкт-Петербург

61 12-4/122

Федеральное государственное оюджетное образовательное учреждение высшего

профессионального образования «Санкт-Петербургский государственный университет»

На правах рукописи

Макашов Сергей Эдуардович

УСТАНОВЛЕНИЕ ЗАКОНОМЕРНОСТЕЙ ФИЛЬТРАЦИИ РАССОЛОВ ДЛЯ РЕШЕНИЯ ПРОБЛЕМ ЗАГРЯЗНЕНИЯ ПРИРОДНЫХ ВОД НА ВЕРХНЕКАМСКОМ МЕСТОРОЖДЕНИИ КАЛИЙНЫХ СОЛЕЙ (на примере шахтных полей 1-го Березниковского и 3-го Соликамского калийных рудоуправлений)

Специальность 25.00.07 - Гидрогеология Диссертация на соискание ученой степени кандидата геолого-минералогических наук

Научный руководитель: Румынии Вячеслав Гениевич доктор геолого-минералогических наук, член-корреспондент РАН.

Санкт-Петербург 2011

Содержание

Введение.

Глава 1. Теоретические основы и экспериментальные модельные исследования процессов плотностной конвекции..................................................................................................................10

1.1. Уравнение миграции жидкости переменной плотности.............................................12

1.2. Физические предпосылки, заложенные в расчетные алгоритмы программ.............13

1.3. Тестирование программ на точных аналитических решениях...................................15

1.4. Моделирование физического эксперимента Шинкариола-Шварца...........................20

1.5. Моделирование субвертикальной миграции соленых вод от поверхностного

источника.........................................................................................................................................^

1.6. Выводы к главе 1............................................................................................................30

Глава 2. Использование численных моделей при поиске новых аппроксимационных

зависимостей....................................................................................................................................01

2.1. Задача о формировании границы раздела разноплотностных жидкостей................31

2.1.1. Существующие зависимости для описания динамики границы раздела при поршневом вытеснении в плановом потоке........................................................................31

2.1.2. Новая аппроксимационная зависимость для случая анизотропного пласта... 38

2.2. Исследование изменения концентрации в откачивающей скважине при субвертикальном подтягивании соленых вод...............................................................................41

2.2.1. Равновесный купол соленых вод под несовершенной скважиной: предельный дебит и критическая величина подъема купола соленых вод...........................................42

2.2.2. Изменение концентрации в откачивающей скважине......................................47

2.3. Выводы к главе 2.............................................................................................................60

Глава 3. Физико-географический очерк и геолого-гидрогеологические условия

Верхнекамского месторождения калийных солей.......................................................................61

3.1. Физико-географический очерк Верхнекамского месторождения калийных солей . 62

3.2. Геологическое строение.................................................................................................64

3.2.1. Стратиграфия и литология...................................................................................65

3.2.2. Структурно-тектонические особенности............................................................68

3.3. Гидрогеологические условия.........................................................................................70

3.4. Гидрогеохимическая зональность подземных вод......................................................77

Глава 4. Исследование влияния аварийного затопления рудника БКПРУ-1 на миграцию рассолов в верхних водоносных горизонтах................................................................................82

4.1. Моделирование формирования естественной гидрохимической зональности........87

4.1.1. Краткая характеристика процесса растворения каменной соли.......................87

4.1.2. Схематизация гидрогеологических условий района БКПРУ-1........................89

4.1.3. Калибровка модели и результаты расчетов........................................................96

4.2. Исследование развития ореола техногенного засоления подземных вод...............102

4.2.1. Характеристика основных источников загрязнения и данные режимных наблюдений...........................................................................................................................ЮЗ

4.2.2. Фильтрационные потери техногенных рассолов из емкостей-накопителей. 108

4.2.3. Моделирование формирования ореола техногенного загрязнения и

109

результаты расчетов.............................................................................................................

4.3. Моделирование геомиграционных процессов в период аварийного затопления рудника...........................................................................................................................................^

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

4.5. Выводы к главе 4...........................................................................................................123

Глава 5. Исследование особенностей развития ореола загрязнения подземных вод

техногенными рассолами в районе СКРУ-3 методами численного моделирования..............125

5.1. Характеристика источников загрязнения...................................................................125

5.2. Данные режимных гидрогеологических и гидрологических наблюдений.............127

5.3. Фильтрационные потери рассолов..............................................................................129

5.4. Выбор программного средства и схематизация гидрогеологических условий района

i о i

шахтного поля СКРУ-3.................................................................................................................

5.6. Результаты калибровки фильтрационной модели.....................................................136

5.7. Результаты решения эпигнозной и прогнозной миграционных задач.....................140

5.8. Исследование взаимовлияния гидродинамического и плотностного градиентов на численной модели..........................................................................................................................144

5.9. Выводы к главе 5...........................................................................................................147

Заключение................................................................................................................................149

Список использованной литературы......................................................................................151

Введение

Актуальность исследований.

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

Основные задачи исследования.

Основные методы исследований.

Настоящие исследования основывались:

-на решении задач динамики подземных вод, описывающих потоки жидкостей переменной плотности аналитическими методами и методами численного (математического) моделирования;

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

Для решения теоретических и практических задач применялись современные программные комплексы SUTRA (Voss and Provost, 2002), SEA WAT (Guo and Langevin, 2002), TOUGH2 (Prues, 1991) и DENSFLOW (П.К. Коносавский, 2002).

Научная новизна.

Фактический материал и личный вклад автора.

Основные положения, выносимые на защиту

2. Формирование естественной гидрохимической зональности надсолевого уфимского водоносного комплекса происходило в течение 8-10 тыс. лет и определялось совместным протеканием процессов растворения галита с поверхности соляной залежи, молекулярной диффузии, профильной гидродисперсии и плотностной конвекции, препятствующей в сочетании с площадной инфильтрацией пресных вод засолению верхней части разреза.

3. Аварийное затопление рудника БКПРУ-1 вызвало изменение конфигурации существующих ореолов техногенных рассолов в надсолевой толще пород и способствовало

снижению выноса хлоридных натриевых солей в поверхностные водотоки, что подтверждается режимными наблюдениями и численным моделированием.

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

Практическая значимость работы:

Апробация работы.

Результаты исследований докладывались на Первой Всероссийской конференции молодых ученых, посвященной памяти В.А. Мироненко «Дерзость Надежд» (Санкт-Петербург, 2010), научной конференции, посвященной 75-летию кафедры гидрогеологии СПбГУ «Комплексные проблемы гидрогеологии» (Санкт-Петербург, 2011), УП-й ежегодной Общероссийской научно-практической конференции «Перспективы развития инженерных изысканий в строительстве в Российской Федерации» (Москва, 2011).

Публикации.

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

Структура и объем работы.

Диссертация объемом 160 страниц, включает введение, 5 глав, заключение, список источников из 113 наименований, содержит 84 рисунка и 17 таблиц.

Благодарности. Автор выражает искреннюю признательность своему научному руководителю д.г.-м.н., члену-корр. РАН В.Г. Румынину за помощь и поддержку на всех этапах работы. За конструктивные консультации и помощь автор благодарит: к.г.-м.н. доц. кафедры гидрогеологии СПбГУ A.A. Потапова, к.г.-м.н., зав. кафедрой П.К. Коносавского и к.г.-м.н., в.н.с. СПбО ИГЭ РАН J1.H. Синдаловского. Особую признательность автор выражает д.г.-м.н. проф. кафедры гидрогеологии и инженерной геологии СПГГУ Р.Э. Дашко за критический анализ работы и ряд предложенных формулировок отдельных ее положений. За моральную поддержку и ценные советы автор искренне благодарит к.г.-м.н. Е.В. Мольского и к.г.-м.н. М.Н. Кочневу. Благодарность за помощь в полевых работах и получении фактических данных по объекту автор выражает В.П. Копшталеву, Л.И. Лесничему и А.О. Воробьевой.

Глава 1. Теоретические основы и экспериментальные модельные исследования процессов плотностной конвекции

Изучение процессов взаимной фильтрации неоднородных по плотности жидкостей в водоносных горизонтах в настоящее время является одним из наиболее значимых направлений в гидрогеологии. Такие процессы осложняются действием дополнительного механизма массопереноса - плотностной (гравитационной) конвекцией. Она обуславливается проявлением дополнительной вертикальной компоненты скорости, которая может оказывать существенное влияние на фильтрационную картину. При этом действие конвективно-дисперсионного механизма массопереноса в процессе фильтрации разноплотностных жидкостей зачастую приобретает второстепенное значение. Для подавляющего большинства задач - при преимущественно горизонтальной стратификации подземных вод по плотности характеристикой гравитационной конвекции служит величина плотностного градиента Ар = -р0)/р0 (Р° и А> - максимальное и минимальное

значение плотности контактирующих жидкостей, кг/м ) [63].

Зависимость плотности (р) жидкости от ее минерализации (с) обычно, для ЫаС1

рассола, описывается линейной функцией [39]: р(с) = р0 +апаЫ(с-с0) где

Информация о работе

Макашов, Сергей Эдуардович
кандидата геолого-минералогических наук
Санкт-Петербург, 2011
ВАК 25.00.07

Диссертация

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

Автореферат

Похожие работы