- •119991, Москва, гсп-1, Ленинский проспект, 6; Издательство мггу; тел. (095) 236-97-80; факс (095) 956-90-40 «ата»
- •Глава 1. Физические основы динамики
- •Глава 4. Исследование задач плановой
- •Глава 7. Применение принципов и методов динамики подземных вод при гидрогеологических опытных работах и наблюдениях 392
- •Глава 8. Использование методов динамики подземных вод при решении гидрогеологических и инженерногеологических проблем разработки месторождений твердых полезных ископаемых 451
- •Глава 1
- •Элементы гидростатики
- •Гидростатический напор
- •Элементы гидродинамики идеальной жидкости
- •Элементы гидродинамики реальной жидкости
- •О режимах движения
- •Общая физическая характеристика водонасыщенных горных пород
- •Геометрия пор и трещин в горных породах
- •Виды воды в горных породах с позиций задач динамики подземных вод
- •Водонасыщенные горные породы как сплошная среда
- •Подземная гидростатика (напряжения в водонасыщенных горных породах)
- •Емкостные свойства горных пород
- •Гравитационная емкость
- •Упругая емкость
- •Основной закон фильтрации и проницаемость горных пород
- •Коэффициент фильтрации и коэффициент проницаемости
- •Ограничения на закон Дарси
- •Общие представления о статистической теории фильтрации
- •О напряженном состоянии горных пород в фильтрационном потоке (гидродинамическое давление)
- •Общая физическая характеристика
- •Физические основы моделирования геофильтрационных процессов
- •Глава 2 | математические основы теории
- •Гидродинамическая типизация условий движения подземных вод
- •Построение основных дифференциальных уравнений геофильтрации и математические основы моделирования фильтрационных процессов
- •Дифференциальные представления исходных физических закономерностей
- •Расчетная модель жесткого режима фильтрации
- •Расчетная модель упругого режима фильтрации
- •Основные дифференциальные уравнения плановой фильтрации
- •Плановая фильтрация в изолированном напорном пласте
- •Плановая напорная фильтрация при наличии перетекания
- •Плановая фильтрация в безнапорном пласте
- •Раздел 1.4), выражением р
- •Математическая модель плановой фильтрации — условия применимости и основные расчетные схемы
- •Об условиях применимости расчетной модели плановой фильтрации
- •Основные расчетные схемы плановой фильтрации
- •Глава 3
- •Плоскопараллельная (одномерная) стационарная фильтрация
- •0 Формуле Дюпюи и промежутке высачивания
- •Безнапорная фильтрация в слоистом пласте между двумя бассейнами (реками) при отсутствии, инфильтрации
- •Напорно-безнапорная фильтрация между двумя
- •Движение в планово-неоднородном напорном пласте
- •Безнапорное движение между двумя бассейнами (реками) в однородном пласте с наклонным водоупором при отсутствии инфильтрации
- •Плоскорадиальная (одномерная) стационарная фильтрация
- •Задача о фильтрации к скважине в круговом пласте
- •Задача о скважине в пласте с перетеканием
- •Решение задач двухмерной установившейся
- •Метод эквивалентных фильтрационных сопротивлений
- •Общие принципы моделирования задач плановой стационарной фильтрации
- •Сплошные модели из электропроводной бумаги
- •Дискретные модели - сетки электрических сопротивлений
- •Простейшие одномерные решения и пути
- •Фундаментальное решение (задача о подпоре вблизи водохранилища)
- •Задача о плоскорадиальной фильтрации к скважине
- •О возможностях распространения решений
- •Аналитическое исследование нестационарных фильтрационных процессов методами интегральных преобразований
- •Моделирование нестационарных плановых потоков
- •Конечно-разностная форма дифференциальных уравнений
- •Аналоговое моделирование нестационарной фильтрации
- •Исходные представления о схемах численного
- •I 4 I Записать и объяснить математические выражения для граничных условий на скважинах, работающих с постоянным расходом и с постоянным напором.
- •Особенности задач, связанных
- •Общая гидродинамическая характеристика
- •Изменения в подземной гидростатике и гидродинамике при опытной откачке
- •Особенности фильтрационных процессов при опытных откачках
- •Основные расчетные схемы
- •Специфика геофильтрационных процессов в различных типовых условиях проведения опытных опробований
- •О некоторых гидрогеоиеханических эффектах
- •Особенности фильтрационного процесса при откачках из планово-ограниченных и планово-неоднородных пластов
- •Анализ влияния технических факторов
- •Значение несовершенства центральной скважины по степени вскрытия пласта
- •Значение несовершенства наблюдательных скважин по степени вскрытия пласта
- •Значение непостоянства расхода откачки
- •Роль скин-эффекта центральной скважины
- •Роль скин-эффекта центральной скважины
- •Инерционность наблюдательных скважин
- •Принципы диагностики данных офр
- •Глава 6 I теория миграции подземных вод 1и основы теории влагопереноса
- •Конвективный перенос в подземных водах
- •Конвективный перенос, осложненный физико-химическими процессами
- •6.1.4. Задача об определении скорости фильтрации скважинной резистивиметрией (термометрией)
- •Молекулярная диффузия и гидродисперсия
- •0 6.2.2. Задана о диффузион
- •Конвективно-дисперсионный перенос в однородных водоносных пластах
- •Фундаментальное решение
- •Задача о запуске пакета индикатора
- •Особенности массопереноса в гетерогенных водоносных системах
- •Общие представления о макродисперсии
- •Макродисперсия в гетерогенных системах упорядоченного строения
- •Макродисперсия в гетерогенных системах неупорядоченного строения
- •Процессы теплопереноса в подземных водах — общие представления и простейшие задачи
- •Об аналогии между процессами тепло- и массопереноса
- •Определение миграционных параметров лабораторными методами
- •Опыты с относительно хорошо проницаемыми грунтами
- •Опыты с относительно слабопроницаемыми грунтами
- •Полевые опытно-миграционные работы
- •Общие вопросы индикаторного опробований водоносных пластов
- •Методика полевого индикаторного опробования
- •11 Мгновенный подъем концентрации индикатора и
- •3 Импульсный ввод — создание больших концентрации индикатора за весьма малый промежуток времени, в течение которого весь индикатор поступает в пласт.
- •Физические основы влагопереноса в горных породах при неполном водонасыщении
- •Общая энергетическая характеристика процесса влагопереноса
- •Закон движения влаги*
- •Постановка и решение простейших задач вертикального влагопереноса
- •Дифференциальное уравнение и граничные условия
- •(Третье равенство); тогда
- •Простейшая задача вертикального просачивания
- •Особенности движения влаги при опробовании пород зоны аэрации наливами в шурфы
- •Глава 7
- •Методика постановки и проведения опытно-фильтрационных работ
- •Виды офо и области их применения
- •Постановка опытных опробований
- •Конструкция и расположение опытных скважин при откачке
- •Режим опытной откачки
- •Продолжительность опытной откачки
- •Определение фильтрационных параметров по данным режимных геофильтрационных наблюдений1
- •Общие представления
- •Прямое определение параметров
- •Прямое определение параметров на основе
- •Об интерпретации данных режимных наблюдений на эвм методами целенаправленного поиска
- •На модели проводится прогнозный расчет первоочередного водоотбора;
- •Методика опытно-миграционных работ1
- •Планирование миграционных опытов
- •Конкретные примеры
- •Общие положения
- •Геофильтрационные наблюдения вблизи бассейнов промышленных стоков
- •Наблюдения за качественным составом подземных вод
- •Общие принципы гидрогеологической схематизации в связи с постановкой опытных работ и наблюдений
- •Принцип непрерывности ггс
- •Принцип адаптации
- •Принцип обратной связи
- •Анализ деформаций и устойчивости пород при горных разработках
- •Осадка толщ горных пород при глубоком водопонижении
- •Оползни бортов карьеров, вызыванные напорными водами
- •Фильтрационные деформации пород вблизи горных выработок
- •Изучение деформаций горных пород над выработанным пространством
- •Обоснование дренажа как метода борьбы
- •Влияние дренажа на напряженное состояние пород в откосах
- •Раздел 8.3.3), нетрудно свести такой расчет к простейшей одномерной задаче о бесконечной цепочке скважин. Для этого используется метод эквивалентных фильтрационных сопротивлений (см. Раздел
- •Дренаж как метод борьбы с фильтрационными деформациями откосов
- •8.2.3. Водопонижение при проходке шахтного ствола
- •8.3.1. Обцая характеристика прогнозной ситуации
- •Прогноз процессов загрязнения подземных вод в горнодобывающих районах
- •Цели прогноза и элементы предварительной схематизации
- •Прогнозные оценки процессов загрязнения подземных вод аналитическими методами
- •Основные представления о математическом ¥ моделировании процессов загрязнения подземных вод
- •Краевые условия фильтрации
Подземная гидростатика (напряжения в водонасыщенных горных породах)
Содержание предыдущего раздела говорит о том, что массив водонасыщенных горных пород представляет собой весьма сложную физическую систему; этим предопределяется и сложный характер механического - силового - взаимодействия между отдельными его элементами. Мерой такого взаимодействия, т.е. мерой внутренних сил, действующих в горной породе, служат напряжения.
ВОПРОС. Как определяется напряжение в курсе общей физики?
Естественно, казалось бы, определить напряжения в горных породах через силы взаимодействия между их элементарными составляющими. Ясно, однако, что если мы, например, выделим в массиве две соприкасающиеся частицы, то напряжения по их контурам оудут неравномерными, с резкой концентрацией на контактах. Поэтому разумно оперировать некоторыми усредненными по достаточно большому объему горной породы напряжениями
г ^
а = lim о),
<*>-*% (1.25)
где о ylG — соответственно, напряжение (усредненное) и общая сила, действующие по площадке <ы;
0)о — минимальная репрезентативная площадка
(удовлетворяющая предпосылке сплошности среды — см. раздел 1.2).
Следовательно, само определение понятия напряжения в горных породах носит статистический характер (см. раздел 1.2.3).
ЗАДАЧА. В сосуде с водой лежат металлические шары (рис. 1.12). Уровень в сосуде понижается от положения АА до положения А^у Меняются ли действующие в системе напряжения?
ОТВЕТ. Вопрос поставлен некорректно. Нужно уточнить, о каких напряжениях идет речь. Напряжения, характеризующие силовое взаимодействие между шарами, в данном примере остаются неизменными (кстати, почему?), а напряжения, отвечающие дав- Риа 112. Схема к оценке лению воды на поверхности шаров, убы- изменений напряжений вают. при снижении свободного
уровня воды
Из примера видно, что в каждой точке водонасыщенной горной породы мы имеем дело с двумя категориями напряжений, действующих:
в минеральном скелете (шары) и
во внутрипоровой жидкости.
Все механические свойства горной породы (ее прочность, характеристики деформируемости и т.д.) практически определяются первой системой напряжений, которые принято называть эффективными (ор.
ЗАДАЧА.На рис. 1.13 показан закрытый сосуд с водой; внутри него поршень (2) из хорошо проницаемого материала, опирающийся на пружины (1). Что произойдет, если постепенно и плавно повышать давление в жидкости с помощью подключенного к сосуду насоса? Что будет при быстром, скачкообразном приросте давления?
Рис. 1.13. Механическая модель, иллюстрирующая механизм нейтральных напряжений
Точнее, — в репрезентативном ее объеме.
Напряжения второй системы, отвечающие гидростатическому давлению в жидкости, принято называть нейтральными (<тн): сами по себе они вызывают лишь всестороннее сжатие отдельных частиц минерального скелета, которым, в силу их малой сжимаемости, мы будем в дальнейшем пренебрегать. В каждой точке водоносного пласта
°»=ПЛ> (1-26)
где hp — пьезометрическая высота.
Рассмотрим теперь реальный безнапорный водоносный горизонт (рис. 1.14) и оценим эффективное напряжение О по единичной го-
2
ризонтальной площадке ft) = 1м на глубине z. Напряжение Од определяется весом минеральных зерен в столбике вышележащих пород CTj уменьшенным за счет взвешивания их подземными водами Gg. Если все единицы объема минерального материала равен А п , а воды
— До, то Gc = An*z*(l — п) , a Ge =Ao‘Z'(l — я), где z (1 — п ) — объем минеральных зерен. Следовательно,
О, = ос-G„ =(Д„ — Д0)(1 -n)-z=y,-z,
где
У в = (Ап -АоУ(1~п) - (1.27)
расчетный вес единицы объема взвешенной горной породы.
Например, для кварцевого песка А „ ~ 2,7 *10 Н/м3; п ~ 0,4,
4/3 п
т.е. ув 10 Н/м . Между тем, после снижения уровня подземных
вод для того же песка с остаточной насыщенностью пор водой около 50% вес единицы объема взвешенной горной породы составит уп «
2,7 * (1 — 0,4) + 1 • 0,4 * 0,5 =1,8* 104 Н/м3. Следовательно, за счет сил взвешивания эффективные напряжения в этом примере уменьшаются почти в два раза. Отсюда понятно, что при снижении уровня воды в пласте эффективные напряжения будут расти и песок будет сжиматься.
Усложним нашу механическую модель и рассмотрим систему тяжелых поршней из слабопроницаемого материала, разделенных пружинами (рис. 1.15а). Если уменьшить давление жидкости в отсеке 4, снизив уровень в трубке 4 на величину 5, то нагрузка на пружину этого отсека возрастет на Aog —yo’S (предполагается, что
из-за слабой проницаемости поршней давление в остальных отсеках
не меняется). В этой модели давление на пружины соответствует эффективным напряжениям в слоях водоносных пород, разделенных слабопроницаемыми слоями, а давление в жидкости — нейтральным напряжениям.
0 * в *
. 'Hi
^60
Рис. 1.14. Схема, иллюстрирующая действие сил гидростатического взвешивания в безнапорном пласте
Для более общего аналитического описания процессов перераспределения напряжений рассмотрим полное общее давление рп, оказываемое столбиком обводненных пород на единичную горизонтальную площадку (IpJ = = I £7Я I, где оп— полное напряжение). Очевидно, оно будет равно весу пород в этом столбике вместе с заключенной в них водой.
Обратим, однако, внимание на одно обстоятельство: здесь мы предполагаем, что вес комплекса пород, лежащих выше рассматриваемой площадки, никак не перераспределяется на боковые зоны; соответственно, если под влиянием этого веса происходят деформации сжатия, то упомянутый комплекс ведет себя как абсолютно гибкая плита (жесткость пород во внимание не принимается). Позже (см. раздел 5.3) мы увидим, что это допущение имеет свои пределы.
Рассматривая выделенный объем как единое (сплошное) тело, мы должны считать, что давление Рп уравновешивается двумя силами: нормальной к площадке реакцией скелета породы оэр и гидростатическим давлением по
площадке р (Ip I = 1сгм1), направленным, согласно известным принципам гидростатики, также по нормали к площадке: рп - оэр + р. Следовательно, напряжение оэ, обусловливающее силы взаимодействия частиц минерального скелета горной породы и равное по величине реакции оэр,
можно определить по формуле
Уравнение (1.28) может рассматриваться как основная закономерность подземной гидростатики. Из него следует, что эффективные напряжения, т.е. напряжения в минеральном скелете, могут меняться не только при из-
S' В
— Тг*’"Т |
|
— |
|
|
ч |
|
.'■V- |
1; |
■ 1 |
|
|
|
Л> 1 |
|
|
|
|
|
|
|
4 |
- |
|
'• 5- |
|
■ |
|
|
4 |
|
|
/ к: |
• 1 :х |
2 |
|
Hi
Ни
Т
<Г ио
а
Рис. 1.15. Схема, иллюстрирующая взаимодействие эффективных и нейтральных напряжений:
а - механическая модель; б - многопластовая напорная система (цифрами на рисунке даны номера отсеков слоев ).
1 - вода; 2 - поршень с пружиной; 3 - водоносный горизонт; 4 - слабопроницаемые слои
a, = |Pj - 'Р 1 =°,-<V (1.28)
менении внешнего (полного) давления на водоносный пласт (дополнительная нагрузка на поверхности земли, выемка вышележащих пород и т.д.), но и при изменении напоров в нем в результате откачки или нагнетания воды.
Рассмотрим единичную горизонтальную площадку див пределах напорного водоносного горизонта 4 (см. рис. 1.15 б). При снижении напора на величину S общее давление на площадку (вес столба АВ) останется практически неизменным (точнее, оно меняется лишь на небольшую величину, отвечающую отдаче воды самим напорным пластом; позднее мы оценим порядок этой величины и убедимся, что, пренебрегая ею, мы вносим погрешность, измеряемую долями процента). Так как Аоп~ 0, то из формулы (1.28) следует
Ao3~-AoH=y0S (1.29)
Следовательно, понижение напоров на каждые 10 м должно приводить к росту напряжений в скелете на 0,1 МПа.
Итак, подземные воды являются важным силовым фактором, определяющим многие деформационные процессы в горных породах. Например, известны многочисленные случаи оседания земной поверхности в районах интенсивного отбора воды или нефти. Так, в большом районе Калифорнии понижение напоров на 150 м привело к осадке поверхности на площади 5000 км , причем максимальное оседание превысило 12 м [44 [.
Аналогичное объяснение можно дать другому интересному явлению — возникновению землетрясений вблизи вновь создаваемых водохранилищ (в том числе и в доселе несейсмичных районах). Рассмотрим для примера два тектонических блока, граничащие вдоль поверхности сброса АВ (рис. 1.16) и находящиеся первоначально в равновесном состоянии. Возникновение водохранилища и связанный с ним рост напоров Н вдоль поверхности сброса снижают действующие по этой поверхности силы трения и вызывают относительное смещение блоков, проявляющееся как землетрясение. И, наоборот, начало глубинные смещений породных масс, связанных с изменениями равновесного давления в горном массиве, может фиксироваться по аномальным изменениям нейтральных напряжений (т.е. напоров в наблюдательных скважинах) и тем самым служить индикатором надвигающегося землетрясения.
Рис. 1.16. Механизм влияния напорных водна сейсмическую устойчивость массива горных пород:
1 - водоносный горизонт; 2 - водонепроницаемые породы
В заключение остановимся на напряжениях в пределах капиллярной каймы (см. раздел 1.2.2). Так как гидростатические давления здесь меньше атмосферного, то в каждой точке капиллярной каймы действуют нейтральные напряжения (см. раздел 1.1.2):
CT„(Z) = -VZ- (1.30)
где z — превышение точки относительно свободной по- верхности воды (0 < z < hK);
hK — высота капиллярного поднятия (мощность капиллярной каймы).
Таким образом, согласно формуле (1.28), в пределах капиллярной каймы эффективные напряжения возрастают за счет капиллярных сил: к горной породе вдоль поверхности менисков как бы приложена дополнительная внешняя нагрузка, равная оэ =yQ'hK.
Высота капиллярного поднятия в песках изменяется десятками сантиметров, а в глинах — метрами. Это обстоятельство имеет важное значение при моделировании подземных вод в фильтрационных лотках: наличие капиллярной каймы, высота которой обычно соизмерима с высотой лотка, сильно затрудняет соблюдение условий подобия модели и натуры. Заметим попутно, что расчет по формуле (1.6) для капилляров с радиусом R - 0,1 мкм (характерный размер
пор в глинистых грунтах) дает значение hK примерно 16 м. В реальных грунтах величина hK, однако, заметно меньше рассчитанных таким образом значений из-за неоднородности порового пространства и наличия отдельных крупных пор — расширений в капиллярах.
Так как капиллярные силы возникают не только на контакте вода-воздух, но и на поверхности раздела жидкостей с различным поверхностным натяжением, то эти силы оказывают, например, заметное влияние на перемещение водонефтяного контакта в задачах нефтепромыслового дела [36].
