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

Котлов, Сергей Николаевич

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

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

КОТЛОВ Сергей Николаевич

НАУЧНО-МЕТОДИЧЕСКОЕ ОБОСНОВАНИЕ ГИДРОГЕОЛОГИЧЕСКИХ ИССЛЕДОВАНИЙ ДЛЯ ОБЕСПЕЧЕНИЯ РАЗРАБОТКИ ОТКРЫТЫМ СПОСОБОМ МЕСТОРОЖДЕНИЙ АЛМАЗОВ В АРХАНГЕЛЬСКОЙ ОБЛАСТИ

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

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

2 [*№ 2011

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

4848996

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

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

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

Антонов Владимир Васильевич

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

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

Павлов Александр Николаевич,

кандидат геолого-минералогических наук

Назима Виктор Николаевич

Ведущее предприятие - ОАО «Гипроруда».

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

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

Автореферат разослан 13 мая 2011 г.

УЧЕНЫЙ СЕКРЕТАРЬ диссертационного совета кандидат геолого-минералогических наук

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

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

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

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

- региональным развитием вскрываемого карьерами мощного слоистого водоносного комплекса падунских отложений (V2pd), прорванного кимберлитовыми трубками,

- наличием закарстованных известняков и доломитов (Сго1-ок) в верхней части вскрываемой карьерами толщи,

- высокими притоками подземных вод (до 6000 м3/час) в карьеры и к системам заградительного дренажа,

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

Большой вклад в развитие методики гидрогеологических исследований на месторождениях твердых полезных ископаемых внесли С.К. Абрамов, О.Б. Скиргелло, В.Д. Бабушкин, В.А. Мироненко, В.М. Шестаков, М.В. Сыроватко, Н.И. Плотников, И.Е. Жернов, Д.И. Щеголев, М.С. Газизов, В.Г. Румынии, Б.В. Боревский, A.A. Рошаль, Г.Н. Гензель и другие. Изучению и прогнозу техногенного режима подземных вод посвящены многочисленные труды Ю.А. Норватова. Однако, вопросы, связанные с вскрытием глубокими карьерами слоистых и неоднородных в плане водоносных толщ, требуют дополнительного анализа гидрогеологических условий освоения месторождений, разработки методических принципов гидрогеологической разведки и обоснования оптимальных дренажных мероприятий на карьерных полях.

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

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

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

- разработка методики проведения опытно-фильтрационных работ (ОФР) для определения фильтрационных параметров неоднородного в плане водоносного горизонта закарстованных доломитов;

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Достоверность и обоснованность научных положений и выводов подтверждается:

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

-результатами специальных опытно-фильтрационных работ, выполненных при доразведке месторождения им. В.П.Гриба;

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

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

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

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

- предложения по распределению сброса дренажных вод в реки и созданию системы гидрогеологического мониторинга рекомендованы институту «Гипроруда» для использования при проектировании ГОКа на базе месторождения им. В.П.Гриба;

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

Личный вклад автора:

-участие в постановке, проведении и интерпретации результатов ОФР на этапе доразведки месторождения им. В.П.Гриба;

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

-определение оптимальных параметров системы заградительного дренажа на карьере месторождения им. ВЛ.Гриба;

-проведение численных экспериментов для установления высоты высачивания подземных вод на откосах в зависимости от степени анизотропии дренируемых комплексов;

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

Апробация работы. Результаты исследований докладывались на: Международной научной школе молодых ученых и специалистов «Проблемы освоения недр в XXI веке глазами молодых» (ИПКОН РАН, Москва, 2006, 2008), VII Межрегиональной научно-практической конференции «Геология, полезные ископаемые и проблемы геоэкологии Башкортостана, Урала и сопредельных территорий» (Минэкологии, Уфа, 2008), Международной конференции «XLIX sesji. Studenckich kót naukowych pionu Górniczego Akademii Górniczo-Hutniczej» (Краков, 2008), ежегодной «Конференции студентов и молодых учёных СПГТИ(ТУ)» (2009, 2010, 2011 гг.), и в институте «Гипроруда» (2010 г.).

Основные положения диссертации отражены в 5 научных публикациях, в том числе 3 статьи в изданиях, входящих в Перечень, рекомендованных ВАК Минобрнауки РФ.

Структура и объем работы. Диссертация изложена на 178 страницах, включая введение, 5 глав, заключение, список литературы из 87 наименований, содержит 52 рисунков, 29 таблиц и 3 приложения.

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

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

В третьей главе изложены разработанные с использованием имитационного численного моделирования рекомендации по методике постановки и интерпретации ОФР при плановой фильтрационной неоднородности водоносных горизонтов и комплексов. Методика проиллюстрирована примером опробования олмуго-окуневского водоносного горизонта (Сго1-ок) на месторождении им. В.П.Гриба.

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

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

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

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

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

дит в долинах рек. Инфильтрационное питание олмуго-окуневского (Сго1-ок), урзугского (Сгиг) водоносных горизонтов и падунского (Угре!) водоносного комплекса в рассматриваемом районе характеризуется особенностями, определяемыми строением толщи доломитов олмуго-окуневской свиты.

В результате выветривания верхняя часть толщи доломитов разрушена до состояния глиноподобной «доломитовой муки», что позволяет рассматривать эту зону в качестве относительного водоупора, мощность которого достигает 5-10 м. Нижняя часть толщи доломитов характеризуется закарстованностью, которая на месторождении им. В.П.Гриба в максимальной степени проявляется на участке между переуглубленными долинами реки Падун и проточного озера Черного (рис. 1).

Рис. 1. Схема питания и разгрузки водоносных комплексов в районе месторождения алмазов им. В.П.Гриба: 1 - уровни грунтовых вод четвертичного водоносного комплекса; 2 - уровни подземных вод олмуго-окуневского водоносного горизонта (С2о1-ок); 3 - напоры подземных вод падунского водоносного комплекса (У2рс1); 4 - зона относительно водоупорных выветрелых доломитов; 5 - зона повышенной трещиноватости (закарстованности) в доломитах олмуго-окуневского водоносного горизонта

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

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

Особенности гидрогеологического строения олмуго-окуневского водоносного горизонта предопределили необходимость его опробования створовой откачкой совместно с ниже залегающим урзугским водоносным горизонтом (С2иг), которое было проведено на стадии доразведки в 2009 г. Интерпретация откачки на численной геофильтрационной модели позволила раздельно оценить изменение проводимости олмуго-окуневско водоносного горизонта в пределах 180-300 м2/сутки и проводимость урзугского водоносного горизонта, равную 70 м2/сутки.

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

Свидетельством стабильности питания подземных вод дочетвертичных водоносных комплексов является практическое постоянство их уровней в годовом цикле. Разность уровней, замеренных по наблюдательным скважинам в меженный и паводковый периоды с 2000 г. по 2004 г., не превышает 0,4 м по олмуго-окуневскому водоносному горизонту и 0,2 м по падунскому водоносному комплексу.

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

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

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

Исследование общих закономерностей техногенного гидродинамического режима, формируемого при вскрытии анизотропных водоносных комплексов горными выработками, выполнено на численной геофильтрационной модели с применением пакета компьютерных программ Processing Modflow. На модели в стационарной постановке имитировались условия формирования вблизи горной выработки двумерного в разрезе безнапорного потока подземных вод в анизотропном массиве, залегающем на горизонтальном водоупоре. Численная модель представлена совокупностью десяти горизонтальных слоев мощностью по 10 м каждый; на границе дренажа этих слоев (на вертикальной стенке выработки) имитировались с применением модуля DRAIN условия на промежутке высачивания (напор в каждой точке соответствует её высотной отметке).

На моделях оценивалась высота промежутка высачивания подземных вод на стенке выработки (Ив) и определялись расходы потока при различной степени анизотропии массива (a=kx/k2, где kxwkz- коэффициенты фильтрации в горизонтальном и вертикальном направлениях соответственно). При моделировании в качестве переменных задавались глубина потоков на границах питания (Аг) и расстояния между границами питания и разгрузки (/). По результатам численных экспериментов получены расчетные графики, отражающие зависимость высоты промежутка высачивания подземных вод от степени анизотропии массива и гидродинамических характеристик потока подземных вод.

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

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

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

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

Следует отметить, что по результатам натурных наблюдений на месторождении им. М.В.Ломоносова высота высачивания подземных вод на бортах карьера, вскрывающего падунский водоносный комплекс, достигает 30-40 м над подошвой карьера при его глубине около 80 м.

Рис. 2. Графики зависимости высоты промежутка высачивания от приведенного расхода при различной степени анизотропии массива (а) 12

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

Анизотропия породных массивов существенно влияет на эффективность водопонижающих скважин - при имитации на численной геофильтрационной модели условий эксплуатации скважины в анизотропном пласте установлено, что разница высот промежутка высачивания и уровней подземных вод вблизи скважины в изотропных и анизотропных толщах может достигать 25-50% от высоты столба воды в скважине в зависимости от степени анизотропии (рис. 3). Дебиты скважин несущественно зависят от степени анизотропии пласта: различие дебитов при рассмотренных степенях анизотропии не превышает 0,5%.

h, м

расстояние от скважины, м

Рис. 3. Положения уровней безнапорного потока (h) вокруг водопонижающей скважины при различной степени анизотропии массива(а)

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

Интерпретация результатов опробования выполнена с применением численной геофильтрационной модели, реализованной в программной среде Relis (автор П.К. Коносавский); полученные значения коэффици-

ентов фильтрации по напластованию и вкрест напластования равны 1,5 м/сутки и 0,13 м/сутки соответственно.

б£ 7-НП 4-ОУ

8-ДПМ1

С,о1-ок

Условные обозначения

с,иг

• 8-ДПМ2 —датчик гидростатического давления и его номер

ц7о Л — линии равных напоров в абсолютных отметках, м

5-ДП5

• «-ДПМ5

о

50

Рис. 4. Схема распределения напоров в падунском водоносном комплексе на момент окончания откачки из скважины 6-П

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

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

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

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

вие: напоры в каждой точке на стенке скважины в пределах промежутка высачивания над уровнем воды в скважине равны высотной отметке (геометрической высоте) этой точки. Реализация этого граничного условия на численной геофильтрационной модели возможна с использованием модуля Drain (в пакете программ РМ 5.3) при минимальном фильтрационном сопротивлении «дрены».

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

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

90

80 -

60 -

20 •

£3 10 -

20 40 60

удельный расход, м2/сут

б) 90 5

Л) 80 о.

0

= 70

1

60

SO -

40

О 30 ■

20

о л) н

8 ю ■

60

65

70 75 напор, М

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

«Модельный» расход скважины отличается от расхода, рассчитанного по формуле Дюпюи, на 1,3%, что можно считать доказательством корректности предлагаемого приема моделирования водопонижающих скважин в безнапорных пластах. Следует особо отметить, что полученный на модели расход неравномерно распределяется вдоль ствола скважины, причем около 30% дебита поступает в скважину в пределах промежутка высачивания.

Полученный на модели максимальный напор по стволу скважины соответствует отметке верхней границы участка высачивания на стенке скважины и высоте этого участка (/?„), равной 34,7 м при заданной высоте столба воды в скважине, равной 40 м. Высота промежутка высачивания, полученная на модели, контролируется решением, предложенным В.М. Шестаковым (различие составляет около 5%).

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

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

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

С использованием численной модели дочетвертичного водоносного комплекса, отражающей его структуру, условия питания и разгрузки подземных вод на площади около 400 км , определены параметры системы водопонижающих скважин с учетом скорости углубки карьера (30 м/год), оценены изменения уровенного режима подземных вод и

ущерб речному стоку в районе месторождения (рис. 6). Определено, что для ограничения притока подземных вод в карьер (400-500 м3/час) необходимо в течение третьего года с начала горных работ введение в эксплуатацию системы заградительного дренажа из 37 водопонижающих скважин, а к концу четвертого года следует задействовать 75 скважин, которые должны эксплуатироваться до завершения горных работ с дебитами около 70 м3/час (суммарная производительность заградительного дренажа около 5400 м3/час).

Для восполнения оцененного на модели ущерба речному стоку обосновано рациональное распределение сброса в реки вод, откачиваемых водопонижающими скважинами (1200 м3/час в озера Волчьи и 4200 м3/час в озеро Черное).

Время от начала горных работ, лет

Рис. 6. Прогнозная структура баланса подземных и поверхностных вод: 1 - расход подземных вод, поступающих в карьер и к системе заградительного дренажа;

2 -расход подземных вод, разгружающихся в речную сеть;

3 -расход разгрузки поверхностных вод в водоносные комплексы (потери из рек);

4 - сработка емкостных запасов водоносных комплексов

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

ляет 1,3). С использованием модели оценен приток в карьер соленых вод, который может достигать 50 м3/час при вскрытии мезенского относительно водоупорного комплекса, содержащего преимущественно хлоридные натриевые воды с минерализацией 5-20 г/дм3. Минерализация этих вод может быть снижена до нормативных значений при их разбавлении водами, откачиваемыми дренажными скважинами.

поверхность скольжения

-300-

-4(КН

Рис. 7. Схема распределения прогнозных напоров и положение поверхности скольжения в западном борту карьера при достижении проектной глубины 450 м; р - гидростатическое давление в точке на поверхности скольжения (р=100 м вод. ст.=1 МПА)

Численная геофильтрационная модель четвертичного водоносного комплекса использована для оценки изменений гидрогеоэкологических условий лесопользования. Установлено, что на момент окончания горных работ на территории, прилегающей к карьеру (около 45 км2), не произойдет существенного снижения уровней грунтовых вод вследствие особых условий связи четвертичного и дочетвертичных водоносных комплексов. За пределами этой площади изменения уровенного режима грунтовых вод будут определяться интенсивностью инфильтрационно-го питания, которая при техногенном режиме может возрасти.

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

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

ЗАКЛЮЧЕНИЕ

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

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

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

- установлены закономерности формирования техногенного режима подземных вод в анизотропных массивах;

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

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

- создана численная геофильтрационная модель района месторождения им. В.П.Гриба;

- разработана методика численного моделирования водопонижаю-щих скважин при оценке оптимальных параметров заградительного дренажа;

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

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

- оценен ущерб речному стоку и обоснованы рекомендации по оптимизации сброса дренажных вод в речную сеть;

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

- разработаны рекомендации по организации гидрогеологического мониторинга в районе месторождения им. В.П.Гриба.

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

В изданиях рекомендованных ВАК Минобрнауки РФ

1. Котлов С.Н. Моделирование условий эксплуатации водопони-жающих скважин в системе заградительного дренажа / Ю.А. Норватов, И.Б.Петрова, С.Н.Котлов // Естественные и технические науки, М.: Спутник Плюс, 2011. №2. С. 235-237.

2. Котлов С.Н. Использование современных компьютерных технологий при обосновании мероприятий, обеспечивающих устойчивость бортов карьеров // Записки Горного института «Проблемы рационального природопользования». СПб, 2011. Т. 189. С. 34-37.

3. Котлов С.Н. Использование имитационного моделирования для планирования и интерпретации опытно-фильтрационных работ при разведке месторождений твердых полезных ископаемых / С.Н. Котлов, К.Е. Володченко // Записки Горного института «Проблемы рационального природопользования». СПб, 2011. Т. 189. С. 38-41.

В других изданиях

А. Котлов С.Н. Гидрогеологические проблемы при освоении алмазных месторождений в Архангельской области / Ю.А. Норватов, И.Б. Петрова, С.Н. Котлов, Д.И. Савельев // Материалы 8-ой международной научно-практической конференции: «Освоение минеральных ресурсов Севера: проблемы и решения», Воркута-Сыктывкар-Ухта. «Народное хозяйство республики Коми», 2010. Т. 19, №2. С. 259-261.

5. Kotlov S.N. Usage of computer technology in the planning of hydro-geological work and interpretation on the stage of detailed exploration / Y.A. Norvatov, I.B. Petrova, S.N. Kotlov, A.V. Yashina // International mining conference "Advanced mining for sustainable development". Vietnam, Ha Long, 2010. P. 616-621.

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

Содержание диссертации, кандидата геолого-минералогических наук, Котлов, Сергей Николаевич

Введение.

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

Глава 2. Анализ условий формирования гидродинамического режима водоносных комплексов с оценкой их фильтрационных параметров.

2.1 Гидрогеологическая стратификация района месторождения.

2.2 Анализ естественного режима подземных вод, условий их питания и разгрузки.

2.2.1 Анализ условий формирования инфильтрационного питания водоносных комплексов.

2.2.2 Определение интенсивности инфильтрационного питания подземных вод.

2.2.3 Оценка проводимости дочетвертичных водоносных комплексов по результатам наблюдений за естественным реэюимом подземных вод.

Глава 3. Планирование и интерпретация опытно-фильтрационных работ в планово неоднородных пластах.

3.1 Общие положения.

3.2 Рекомендации по методике постановки и интерпретации ОФР.

3.3 Методика опытно-фильтрационного опробования пласта с плановой фильтрационной неоднородностью.

Глава 4. Исследование влияния анизотропии водоносных комплексов на условия формирования техногенного режима подземных вод.

4.1 Аналитические методы учета анизотропии.

4.2 Исследование условий разгрузки подземных вод на откосах бортов карьеров, вскрывающих анизотропные массивы.

4.3 Исследование влияния анизотропии водоносных комплексов на условия эксплуатации водопонижающих скважин.

4.4 Опытно-фильтрационное опробование анизотропных водоносных толщ и методика интерпретации его результатов.

Глава 5. Прогноз с применением численной геофильтрационной модели техногенного режима подземных вод и обоснование параметров системы заградительного дренажа.

5.1 Структура и параметры численной геофильтрационной модели.

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

5.1.2 Вычислительная схематизация.

5.1.3 Параметрическое обеспечение модели.

5.2 Методика калибрации модели.

5.3 Постановка прогнозных гидрогеологических задач, решаемых с использованием численных геофильтрационных моделей.

5.4 Разработка методики моделирования водопонижающих скважин в системе заградительного дренажа.

5.4 Результаты моделирования и практические рекомендации.

5.4.1 Определение оптимальных параметров дренажной системы.

5.4.2 Оценка распределения гидростатических напоров в прибортовых массивах.

5.5 Оценка гидрогеоэкологического воздействия дренажных мероприятий на природную среду.

5.6 Рекомендации по организации гидрогеологического мониторинга.

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

Актуальность работы связана с необходимостью освоения месторождений алмазов в Архангельской области, где открыта перспективная алмазоносная провинция, прогнозные запасы алмазов которой составляют около 800 млн. карат [86]. В настоящее время месторождение алмазов им. М.В.Ломоносова уже разрабатывается открытым способом, а на месторождение им. В.П.Гриба начались вскрышные работы.

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

- региональным развитием вскрываемого карьерами мощного слоистого водоносного комплекса падунских отложений (V2pd), прорванного кимберлитовыми трубками,

- наличием закарстованных известняков и доломитов (Сго1-ок) в верхней части вскрываемой карьерами толщи, 3

- высокими притоками подземных вод (до 6000 м /час) в карьеры и к системам заградительного дренажа,

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

Большой вклад в развитие методики гидрогеологических исследований на месторождениях твердых полезных ископаемых внесли С.К. Абрамов, О.Б. Скиргелло, В.Д. Бабушкин, В.А. Мироненко, В.М. Шестаков, М.В. Сыроватко, Н.И. Плотников, И.Е. Жернов, Д.И. Щеголев, М.С. Газизов, В.Г. Румынии, Б.В. Боревский, A.A. Рошаль, Г.Н. Гензель и другие. Изучению и прогнозу техногенного режима подземных вод посвящены многочисленные труды Ю.А. Норватова. Однако, вопросы, связанные с вскрытием глубокими карьерами слоистых и неоднородных в плане водоносных толщ, требуют дополнительного анализа гидрогеологических условий освоения месторождений, разработки методических принципов гидрогеологической разведки и обоснования оптимальных дренажных мероприятий на карьерных полях.

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

Задачи исследований:

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

- разработка методики проведения опытно-фильтрационных работ (ОФР) для определения фильтрационных параметров неоднородного в плане водоносного горизонта закарстованных доломитов;

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

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

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

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

Методы исследований

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

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

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

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

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

Научная новизна работы:

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

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

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

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

Личный вклад автора:

- участие в постановке, проведении и интерпретации результатов ОФР на этапе доразведки месторождения им. В.П.Гриба;

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

- определение оптимальных параметров системы заградительного дренажа на карьере месторождения им. В.П.Гриба;

- проведение численных экспериментов для установления высоты высачивания подземных вод на откосах в зависимости от степени анизотропии дренируемых комплексов;

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

Достоверность и обоснованность научных положений и выводов подтверждается:

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

- результатами специальных опытно-фильтрационных работ, выполненных при доразведке месторождения им. В.П.Гриба; ;

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

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

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

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

- предложения по распределению сброса дренажных вод в реки и созданию системы гидрогеологического мониторинга рекомендованы институту «Гипроруда» для использования при проектировании ГОКа на базе месторождения им. В.П.Гриба;

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

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

Результаты исследований докладывались на: Международной научной школе молодых ученых и специалистов «Проблемы освоения недр в XXI веке глазами молодых» (ИПКОН РАН, Москва, 2006, 2008), VII Межрегиональной научно-практической конференции «Геология, полезные ископаемые и проблемы геоэкологии Башкортостана, Урала и сопредельных территорий» (Минэкологии, Уфа, 2008), Международной конференции «ХЫХ зе8)'1. I

ЗШбепсИсЬ ко1 паикошусЬ рюпи Оогшсге§о Акаёетп 06гтсго-Ни1;тсгер> (Краков, 2008), ежегодной «Конференции студентов и молодых учёных СПГГИ (ТУ)» (2009, 2010, 2011 гг.), и в институте «Гипроруда» (2010 г.).

Публикации

Основные положения диссертации отражены в 5 научных публикациях, в том числе 3 статьи в изданиях, входящих в Перечень, рекомендованных ВАК Минобрнауки РФ.

Объем работы

Диссертация изложена на 178 страницах, включая введение, 5 глав, заключение, список литературы из 87 наименований, содержит 52 рисунков, 29 таблиц и 3 приложения.

Заключение Диссертация по теме "Гидрогеология", Котлов, Сергей Николаевич

Заключение

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

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

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

- установлены закономерности формирования техногенного режима подземных вод в анизотропных массивах;

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

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

- создана численная геофильтрационная модель района месторождения им. В.П.Гриба;

- разработана методика численного моделирования водопонижающих скважин при оценке оптимальных параметров заградительного дренажа;

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

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

- оценен ущерб речному стоку и обоснованы рекомендации по оптимизации сброса дренажных вод в речную сеть;

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

- разработаны рекомендации по организации гидрогеологического мониторинга в районе месторождения им. В.П.Гриба.

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

1. Абрамов C.K., Скиргелло О.Б. Осушение шахтных и карьерных полей. М.: Недра, 1968, 254 с.

2. Абрамов С.К., Скиргелло О.Б. Способы, системы и расчеты осушения шахтных и карьерных полей. М.: Недра, 1968, 254 с.

3. АзизХ., СеттариЭ. Математическое моделирование пластовых систем. Перевод с английского. Москва-Ижевск: 2004, 416 с.

4. Антонов В.В. Вопросы анализа фильтрационной неоднородности водоносных пластов/ Антонов В.В., Мироненко В.А.// Водные ресурсы, 1977, №4, с. 92-100.

5. Аравин В.И., Нумеров С.Н. Теория движения жидкостей и газов в недеформируемой среде. Гостоптехиздат, М. 1953, 616 с.

6. Атрощенко Ф.Г. Оптимизация работы системы осушения и захоронения ~ дренажных рассолов месторождения трубки «Мир»/ Атрощенко Ф.Г., Филин P.A. //СПб. : Записки горного института, 2003, том 153, с. 126-128.

7. БабинецА.Е., ОгняникН.С., Телыма C.B. О создании постоянно действующих математических моделей гидрогеологических объектов и применении метода конечных элементов. // Геол.журн., 1977, 37, вып.4, с. 119122.

8. Бабушкин В.Д., Плотников И.И., Чуйко В.М. Методы изучения фильтрационных свойств неоднородных пород. М.: Недра, 1974, 208 с.

9. Бабушкин В.Д., Лебедянская З.П., Леви Л.З., Кашковский Г.Н., Боревский Б.В., Плотников И.И. Прогноз водопритоков в горные выработки и водозаборы подземных вод в трещиноватых и закарстованных породах. М.: Недра, 1972, 196 с.

10. Биндеман H.H. Методы определения водопроницаемости горных пород откачками, наливами и нагнетаниями. М.: Углетехиздат, 1951, 52 с.

11. Боревский Б.В. Опыт определения расчетных гидрогеологических параметров по данным групповых откачек/ Боревский Б.В., Язвин Л.С.// Разведка и охрана недр, №4, 1963.

12. Боревский Б.В., Самсонов Б.Г., Язвин JI.C. Методика определения параметров водоносных горизонтов по данным откачек. М.: Недра, 1979, 326 с.

13. Боревский Б.В. и др. Обоснование гидрогеологических моделей месторождений подземных вод многопластовых водоносных систем // Водные ресурсы, 1985. №4.

14. Боревский JI.B., Язвин JI.C. К методике определения гидрогеологических параметров в неоднородных в плане пластах// «Труды ВСЕГИНГЕО», вып. 32, 1970, с. 51-56.

15. Бочевер Ф.М. Гидрогеологические расчеты осушения при разработке месторождений полезных ископаемых// «Разведка и охрана недр», № 8, 1959.

16. Бочевер Ф.М., Гармонов И.В., Лебедев A.B., Шестаков В.М. Основы гидрогеологических расчетов. М.: Недра, 1965, 306 с.

17. Бочевер Ф.М., Гармонов И.В., Лебедев A.B., Шестаков В.М. Основы гидрогеологических расчетов (издание 2-е) М.: Недра, 1969, 368 с.

18. ВержакД.В., Гаранин К.В. Алмазные месторождения Архангельской области и экологические проблемы их освоения // Вестник Московского университета. Сер. 4. Геология. 2005. № 6. С. 18-27.

19. Веригин H.H. Методы определения фильтрационных свойств горных пород. М.: Госстройиздат, 1962, 180 с.

20. Гавич И.К. Гидрогеодинамика. М.: Недра, 1988, 349 с.

21. Гавич И.К. Разведочное (имитационное) моделирование при оптимизации гидрогеологических изысканий/ Гавич И.К., Кожетев В.А., Перцовский В.В.// Водные ресурсы, 1985, №4. С. 86-92.

22. Гавич И.К. Теория и практика применения моделирования в гидрогеологии. М.: Недра, 1980. 357 с.

23. Гензель Г.И. Решение задач охраны подземных вод на численных моделях/ Г. И. Гензель, Н. Ф. Караченцев, П. К. Коносавский и др. Под ред. В. А. Мироненко. М.: Недра, 1992, 240 с.

24. Гидрогеология СССР. Том XLIV. Архангельская и Вологодская области. Под ред. A.A. Макавеева. М: Недра, 1969. 300 с.

25. Головин H.H. К проблеме геологического и структурно-тектонического районирования Архангельской алмазоносной провинции. Изв. Высш. Учеб. Заведений, геология и разведка. 2004.

26. Зеегофер Ю.О., Клюквин А.Н., Пашковский И.С., Рошаль A.A. Постоянно действующие модели гидролитосферы территорий городских агломераций. М.: Наука, 1991, 198 с.

27. Иванов И.П. Инженерная геология месторождений полезных ископаемых. Учебник для вузов. М.: Недра, 1990, 302 с.

28. Климентов П.П., Сыроватко М.В. Гидрогеология месторождений твердых полезных ископаемых. Ч. И. М.: Недра, 1966, 377 с.

29. Козлов B.C. Расчеты дренажных сооружений. Стройиздат, 1940.

30. Коносавский П.К., Соловейчик К.А. Математическое моделирование геофильтрационных процессов.// Санкт-Петербургский технический ун-т. СПб, 2001.

31. Ломакин Е.А., Мироненко В.А., Шестаков В.М. Численное моделирование геофильтрации. М.: Недра, 1988, 228 с.

32. Максимов В.А. О неустановившемся притоке упругой жидкости к скважинам в неоднородной среде. СО АН СССР. «ПМТФ», 1962, №3, с. 7685.

33. Мироненко В.А. Динамика подземных вод. М.: Издательство МГУ, 2005, 519 с.

34. Мироненко В.А. Руководство по дренированию карьерных полей. Л.: ВНИМИ, 1968, 171 с.

35. Мироненко В.А., Мольский Е.В., Румынии В.Г. Горнопромышленная гидрогеология. М.: Недра, 1989, 287 с.

36. Мироненко В.А., Норватов Ю.А., Бокий Л.Л. Фильтрационные расчеты осушения карьерных полей. Материалы к «Методическому пособию по дренажу месторождений полезных ископаемых, подлежащих разработке открытым способом». Часть II. Л.: ВНИМИ, 1965, 91 с.

37. Мироненко В.А., Норватов Ю.А., Сердюков Л.И., Бокий Л.Л., Стрельский Ф.П., Крячко О.Ю., Рюмин А.Н., Мольекий Е.В. Гидрогеологические исследования в горном деле. М.: Недра, 1976, 352 с.

38. Мироненко В.А., Румынии В.Г. Опытно-миграционные работы в водоносных пластах. М.: Недра, 1986. 240 с.

39. Мироненко В.А., Румынии В.Г. Проблемы гидрогеоэкологии. Том 1. Теоретическое изучение и моделирование геомиграционных процессов. М.: Издательство Московского государственного горного университета, 1998, 611с.

40. Мироненко В.А. Обработка результатов опытных откачек, проведенных в условиях двухслойной толщи/ Мироненко В.А., Сердюков Л.И. // «Разведка и охрана недр», 1968, №10, с. 34-38.

41. Мироненко В.А., Фисенко Г.Л., Норватов Ю.А., Бокий Л.Л. Осушение карьерных полей. Материалы к «Методическому пособию по дренажу месторождений полезных ископаемых, подлежащих разработке открытым способом». Часть I. Л.: ВНИМИ, 1965, 79 с.

42. Мироненко В.А., Шестаков В.М. Теория и методы интерпретации опытно-фильтрационных работ. М.: Недра, 1978, 325 с.

43. Михайлов Г.К. Применение модели параллельно анизотропных грунтов, для оценки решений некоторых краевых задач о движении потока грунтовых по водоупору. Инж. Сб. АН СССР, т. XV, 1953.

44. Михайлов Г.К. Упрощение способа расчета фильтрации в однородно-анизотропном грунте//Инж. Сб. 19, 1954, с. 159-160.

45. Норватов Ю.А. Изучение и прогноз техногенного режима подземных вод. Л. «Недра», 1988, 261 с.

46. Норватов Ю.А. Некоторые вопросы методики и техники электромоделирования неустановившейся безнапорной фильтрации// Сб.: Горное давление сдвижения горных пород и методика маркшейдерских работ (ВНИМИ, 66), 1966, С. 154-158.

47. Норватов Ю.А. Перспективы снижения затрат на гидрогеологические работы при разведке месторождений ТПИ/ Норватов Ю.А., Петрова И.Б., Назима В.В. Котлов О.Н.// Разведка и охрана недр. 2006, №1. С. 59-62.

48. Норватов Ю.А. Численное моделирование геофильтрации при планировании мероприятий по обеспечению эффективности и безопасности горных работ/ Норватов Ю.А., Петрова И.Б. // Сборник «Проблемы геодинамической безопасности» ВНИМИ, С-Пб, 1997.

49. РацМ.В. Неоднородность горных пород и их физических свойств. М.: Наука, 1968.

50. Русанов И.В. Геофильтрационное обоснование систем горизонтальных дренажных скважин на карьерах. Автореферат кандидатской диссертации. СПб.: Изд. СПГГИ, 1997, 20 с.

51. Самсонов Б.Г. Определение гидрогеологических параметров при эффективной неоднородности водоносных горизонтов/ Самсонов Б.Г., Зильберштейн Б.М., Бурдакова O.JI.// ОНТИ ВИЭМС. Экспресс-информация. 1972, №4.

52. Синдаловский JI.H. Программная поддержка интерпретации опытно-фильтрационных работ на участках загрязнения подземных вод. Автореферат кандидатской диссертации. СПб.: Изд. СПГГИ, 1998.

53. Синдаловский JI.H. Справочник аналитических решений для интерпретации опытно-фильтрационных опробований. Издательство Санкт-Петербургского университета, 2006, 768 с.

54. Справочник по осушению горных пород. Под ред. И.К. Станченко. М.: Недра, 1984, 572 с.

55. Стрельский Ф.П. Влияние скорости поступательного перемещения контура стока на приток воды из осушаемого пласта// JI., изд. ВНИМИ, 1967 (Труды ВНИМИ, сб. 67), С. 320-325.

56. Полубаринова-Кочина П.Я. Теория движения грунтовых вод. М.: Гостоптехиздат, 1952, 677 с.

57. Полубаринова-Кочина П.Я. Теория движения грунтовых вод. М.: Наука, 1977, 664 с.

58. Полубаринова-Кочина П.Я., Пряжкина В.Г., ЭмихВ.Н. Математические методы в вопросах орошения. М.: Наука, 1969, 414 с.

59. Пыхачёв Г.Б. О дебите скважин в неоднороднопроницаемом пласте// «Труды ГНИ и ГрозНИИ», Грозный, вып. 1, 1944.

60. Чарный И.А. Подземная гидрогазодинамика. М.: Гостоптехиздат, 1963, 396 с.

61. Чарный И.А. Строгое доказательство формулы Дюпюи для безнапорной фильтрации с промежутком высачивания. Докл. АН СССР, т. 79, №6, 1951.

62. Шестаков В.М. Гидрогеодинамика. М.:изд-во МГУ, 1995, 368 с.

63. Шестаков В.М. Динамика подземных вод. М.: МГУ, 1979, 368 с.

64. Щелкачев В.Н., ЛапукБ.Б. Подземная гидравлика. М.: Гостоптехиздат, 1949, 523 с.

65. Chiang Н., Kinzelbach W. Processing MODFLOW. Version 3.0. Hamburg -Heidelberg. 1992-1993.

66. FenskeP.R. Unsteady drawdown in the presence of a linear discontinuity// Groundwater Hydraulics. 1984, №9. P. 1-7.

67. Hantush M.S. Analysis of Data from Pumping Tests in Leaky Aquifers, Trans. Am. Geophys. Union, 37, № 6, 702-714, 1956.

68. Harbauf, A.W. & McDonald, M.G. User's documentation for MODFLOW-96, an update to the US Geological Survey modular finite-difference ground-water flow model. Reston, Virginia: U.S. Geological Survey, Open-File Report, 1996, 96-495.

69. Hill, M. C. MODFLOW/P A computer program for estimating parameters of a transient, three-dimensional, groundwater flow model using nonlinear regression, U.S. Geological Survey, Open-file report, 1992, 91-484.

70. McDonald M.G., Harbough A.W., 1988 MODFLOW, A Modular D3 Finite-Difference Ground-Water Flow Model US GS. Tec. Water-Resources Inv. Bk6, Chap Al. Washington. DC.

71. Neuman Sh.P. Analysis of Pumping Test Data from Anisotropic Unconfined Aquifers. Water Resources Research, v. 11, N 2, 1975, pp. 329-342.

72. Thiery D. Tutorial for the WinMarthe v4.0 pre-processor. BRGM, 2007, 93 p.

73. Wiest R. Geohydrology. New York, London, Sydney, 1965.86. www.dvinainform.ru/arhobl.87. www.severalmaz.ru.

74. Результаты расчета фильтрационных параметров, полученных при обработке данных откачек из скважин 1-ОУ-4-ОУ по схеме Ньюмана

75. М§людательные ^-чжважины скважины\. возмущающие"-\ 1-ОУ 2-ОУ 3-ОУ 4-ОУкх=6 м/сут. &х=6 м/сут. &х=18 м/сут.к2=0,2 м/сут. &2=0,13 м/сут. к2=0,03 м/сут.

76. ОУ — /¿*=0,002 //*=0,0013 //*=0,00170,07 //=0,09 //=0,007

77. Г=270 м2/сут. Г=270 м2/сут. 7=810 м2/сут.х=8 м/сут. кх=4 м/сут. &х=8 м/сут.0,08 м/сут. к-£= 0,5 м/сут. к2=0,25 м/сут.

78. ОУ /¿*=0,0005 — //*=0,0018 //*=0,0013г=0,03 /¿=0,3 //=0,12

79. Т=360 м2/сут. 7М80 м2/сут. Г=360 м2/сут.кх= 12 м/сут. &х=8 м/сут. £х=7 м/сут.к2=0,04 м/сут. к2=0,03 м/сут. ку=0,2Ъ м/сут.

80. ОУ //*=0,0008 /¿*=0,001 — //*=0,00160,016 //-0,01 //=0,6

81. Г=540 м2/сут. Г=360 м2/су т. Г=315 м2/сут.х=10 м/сут. /с2=0,8 м/сут. /сх=8 м/сут. ^=0,012 м/сут &х=3 м/сут. Аг2=0,08 м/сут.

82. ОУ //*=0,0007 /¿=0,01 Г=450 м2/сут. //*=0,0004 //=0,01 Г=360 м2/сут. //*=0,0006 /¿=0,1 Г=135 м2/сут.

83. Примечание: //* — коэффициент упругоемкости, // гравитационная водоотдача, Т- проводимость комплекса.

84. Продолжение приложения 2 Результаты расчета фильтрационных параметров при обработке восстановления уровней после откачки из скважин 1-ОУ—4-ОУ по схеме Ньюмана

85. Наблюдательные ^чжважины скважины\. возмущающие"\ 1-ОУ 2-ОУ 3-ОУ 4-ОУкх=5 м/сут. кК-5 м/сут. &х=5 м/сут.

86. А:7=0,16 м/сут. ¿2=0,14 м/сут. ¿2=0,14 м/сут.

87. ОУ — //♦=0,0009 //♦=0,0008 //♦=0,00080,027 //=0,013 //=0,013

88. Т=225 м2/сут. Т=225 м2/сут. Т=225 м2/сут.

89. ОУ кх-1 м/сут. А2=0Д 1 м/сут. //♦=0,002 /¿=0,005 Г=315 м2/сут. &Х=Ю м/сут. &2=0,4 м/сут. //♦=0,0015 //=0,009 Г=450 м2/сут. ¿х=19,5 м/сут. ^2=0,015 м/сут. //♦=0,0004 //=0,0009кх~9 м/сут. кх=1 м/сут. &х= 12 м/сут.0,08 м/сут. ^=0,04 м/сут. ¿2=0,05 м/сут.

90. ОУ //*=0,0005 //♦=0,0008 — //♦=0,0050,01 //=0,01 //=0,03

91. Г=405 м2/сут. Г=315 м2/сут. Г=540 м2/сут.

92. А:х=4,7 м/сут. &х=3 м/сут.2=0,05 м/сут. ^2=0,08 м/сут.

93. ОУ //♦=0,000075 ¡и=0,008 Г=211,5 м2/сут. //♦=0,0004 //=0,03 7М35 м2/сут.

94. Модельная карта гидроизогипс олмуго-окуневского водоносного горизонта1.ый слой модели)1. Верх Ернозеро

95. Модельная карта гидроизогипс урзугского водоносного горизонта3.ий слой модели)

96. Модельная карта гидроизогипс падунского водоносного комплекса6.ой слой модели)1. Верх Ернозсро

97. Модельная карта гидроизогипс мезенского относительно водоупорного комплекса12.ый слой модели)