- •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. Обцая характеристика прогнозной ситуации
- •Прогноз процессов загрязнения подземных вод в горнодобывающих районах
- •Цели прогноза и элементы предварительной схематизации
- •Прогнозные оценки процессов загрязнения подземных вод аналитическими методами
- •Основные представления о математическом ¥ моделировании процессов загрязнения подземных вод
- •Краевые условия фильтрации
Метод эквивалентных фильтрационных сопротивлений
При движении жидкости через горные породы ей приходится затрачивать энергию (напор) на преодоление сил сопротивления. Показателем этих сил может служить потеря напора на некотором участке потока, отнесенная к расходу жидкости через этот участок:
_ АЯ
Ф~ Q"• (3.53)
Величина Ф называется фильтрационным сопротивлением потока (на выделенном участке). Логика этого названия становится еще более ясной, если вспомнить электрогидродинамическую аналогию (см. раздел 1.6): аналогом для А Н является напряжениеД U, аналогом расхода — сила тока /, и тогда, в соответствии с уравнением
и законом Ома, для величины Ф аналогом является электрическое сопротивление R.
Из выведенных нами формул для одномерного движения нетрудно получить выражения для соответствующих значений фильтрационных сопротивлений. Например, для напорной плоскопараллельной фильтрации, описываемой фюрмулой (3.2),
Ф=Я'-Я>- А.
дБ ТВ (3.54)
где В —фронт потока, т.е. фильтрационное сопротивление при том же общем расходе потока, естественно, растет с увеличением длины пути фильтрации и с уменьшением проницаемости или фронта потока.
Для плоскорадиальной фильтрации из формулы
следует:
// - Я, 1 Я,
Ф = -г7Г- - = т- ~ In -2
Qc 2лТ ГС (3.54а)
Следовательно, при одинаковых перепадах напоров на границах выделенных участков расход плоскопараллельного потока равен расходу плоскорадиального потока при условии
Отсюда возникает возможность сопоставления и взаимной замены потоков с разной геометрией, но с равными фильтрационными сопротивлениями, по крайней мере в тех случаях, когда не принимаются во внимание емкостные запасы пласта или его площадное питание. Реализация этой идеи особенно полезна для потоков, структура которых имеет узколокальные усложнения.
Пусть, например, фильтрация к реке (рис. 3.12) носит преимущественно плановый характер (линии тока горизонтальны) , и лишь вблизи реки имеется участок протяженностью l-т (см. раздел 2.5) с заметным проявлением вертикальной составляющей скорости, т.е. линии тока искривляются, поперечное сечение потока сокращается, и вода испытывает при движении большее сопротивление, чем при плановом движении на участке длиной I. Если мы хотим учесть это обстоятельство, оставаясь в рамках плановой модели, то в ней необходимо изменить фильтрационное сопротивление пласта на участке, прилегающем к реке. Для этого, согласно формуле (3.54), можно увеличить длину пути фильтрации L на некоторую величину A L и вести расчет по фиктивной модели пласта, в которой вместо действительного расстояния L (см. рис. 3.12) фигурирует расчетная величина L + AL ; при этом расход потока и напоры в пределах области х> т (см. раздел 2.5) будут определяться точно так же, как и в обычной случае плановой фильтрации. Аналогично, для потока вблизи реки с закольматированным руслом (см. рис. 2.14) вводится фиктивная величина A L, пропорциональная тп/кп, где тп — мощность кольматационного слоя, кп — коэффициент фильтрации. Так, ранее приведенное граничное условие третьего рода (2.49) после умножения на мощность пласта m нетрудно представить в виде
77777/ Т7р7УутГ/7777Т777777^*
V:
pzZ±_4.
Puc. 3. 72. Схема искривления линий тока вблизи реки с несовершенным руслом
В общем случае величина ALрассматривается и определяется в полевых условиях как специальный параметр, характеризующий фильтрационное сопротивление подруеловых отложений.
ЗАДАЧА. Пользуясь результатами замеров уровней в зимний период (когда уровни слабо меняются) по двум наблюдательным скважинам Н1 иН2 и по водомерному посту в реке ТУ (все замерные точки в одном створе, перпендикулярном к реке), вывести формулу для определения сопротивления ложа реки [34 ]:
Н, -н
A L н2~Нх (*2 *i)
где Xj и х2 — расстояния от контура реки до первой и второй скважины соответственно.
Рис. 3.13. Схема к расчету контура скважин в неограниченном пласте
Применим эту идею к расчету систем скважин, расположенных по некоторому контуру. Для этого рассмотрим сначала задачу о прямолинейном контуре из бесконечно большого числа скважин в неограниченном пласте (рис. 3.13).
Скважины, удаленные на расстояние а друг от друга, имеют одинаковые расходы Qc. Вследствие симметрии
достаточно рассмотреть полосу АА*ВВ1 шириной^. Будем
искать выражения для напоров по линиям АА [Н/х) ] и ВВ1 [Н2(*)], так как физически ясно, что, при данном х, (х) отвечает минимальному значению напора, а #2(х) — максимальному.
Используем метод сложения течений. Согласно (3.44), при Г, =у^+(7^ Г2 =Г2’ =УгР+ (20)?ит.д.
Постоянную Си найдем из условия Я7 (jc) = Н£ на стенке скважины (jc =г ), считая, что0»г:
С t
с =я.
е,
И С 2Л‘Т
Тоща
дсП Ix^ind)2] H\ix) ~ Нс + 2 jtT ln лГ“
г.ПМ2
дсяА Л . xblty(J2\ Q. О И + 2^2 j
Л
Orc
= tf+_^L_ln 1 П7Г
с 2Л'Т
где
через f[ обозначено произведение членов, соответствующих
разным номерам £т 1 до п. Далее имеем :
Z —Z
е—е
1 +
= shz =
пЪс2
lim zf[
и, следовательно, при п -1■ а>
Qc
о е~ех
(3.57)
",(*) = "с + 2^.г1п S7
где х —
Л'Х О '
Аналогично получаем
Qc
(3.58)
Фиксируем некоторое х = хи покажем, что при определенных условиях разность Я2 (xQ) — Я1 (хо) мала в сравнении с характерным (максимальным) перепадом напоров Я2(хо) — Я£. Если бы напор по всей линии СС' равнялся Я^ а по линии х *= xQ — величине Я2(хо), то, согласно (3.2), двухсторонний притокQ' в полосе шириной О составил бы
Q' - 2оТ
(3.59)
Очевидно, реальный расход QC<Q\ т.е. Следовательно, с учетом выражений (3.57) и (3.58)
H2(xo)~Hi(xo) ^ Н2(хо) ~Hi(xo) _ 1 , ехо+е х
-О
ехо-е~х
~TTiFJ-TJT "WT2;■
(3.60)
Если х > JT, т.е. х >(7, то
<0,1%.
яг(*о) ~Н\ (хо) Н2{хо)~Нс
Итак, при х > (7, Я2(х) ~Я1 (х) ~Н (х) = const, т.е. на удалении от ряда, большем (7, линии х * xQ - const — суть линии равных напоров, и здесь имеет место одномерная фильтрация.
Найдем теперь уровни по линии скважин:
я(0,у) =2^ 0пу + 1п(а-у) +
+ 1п((7+у) +... + In (/t<7—у) + 1п(/ш + у)] +СИ =
О уП t1
2л: Г
+ ЯС,
lim z „-*• 00 1
= sin z
ft
Но
<3.61)
Найдем средневзвешенный набор на линии скважин Нф:
0/2
/ H(0,y)dy
Нср(9’ У) -Щ 6/2 нс+ 2ji t 1п ЯгТ +
ГС
(3.62)
где dc — диаметр скважины.
Здесь использовано преобразование [16 ]:
О/l п дЛ/2 (jJt/2
/ In sin dy=~f In sin zdz~-=.f In sin zdz —
° %r/CT %r/CT
=s(-fln2)=-fln2-
Характерный вид кривой #(0, у), т.е. пьезометрической кривой по линии скважин, показан на рис. 3.14. Из
рисунка видно, что на большей части интервала 0 - ^
величина #(0, у) близка к значению Нф. Отсюда возникает идея заменить ряд скважин условной сплошной горизонтальной дреной — узкой траншеей с постоянным напором Нф и равномерно распределенным по ее длине двухсторонним расходом
n-Q‘
9—гг- (3.63)
Движение воды к такой траншее носит, очевидно, одномерный характер, и согласно формуле (3.4) получаем
Qc. = tBxI-hjp
2 о х ’ (3.64)
но при х >осогласно (3.57) и (3.58) Ос
Я(дс) = Яс+1^ l¥+l“^) =Нф+ТоТх
г пг f *п — £
(3.65)
и из уравнения (3.64) получаем
Н(х) - Нс
q_Qc_T Н{х)-Нс
q _Qc 2 2 сг
с Т пТ ndc (3.66)
Рыс. 3.14. Разрез по линии скважин
, «1 а дс In —j 71 nd
х , а л а In
Из последних формул ясно, что эквивалентность воздействий от расчетной траншеи и от реального ряда скважин (при х > а), может быть достигнута двумя путями:
1 введением на контуре траншеи условной) напора Нф, превышающего напор в скважинах Нс на величину
1Ф
а
Я.-Я =ДЯ=А In
2пТ ndc’
(3.67)
|~2 условным увеличением длины пути фильтрации х до значения х + Ах, где
или, что равносильно, — введением дополнительного фильтрационного сопротивления (см. формулы (3.54) и (3.66)):
ДФ =-Ц=Ы-^г.
яТ ndc (3.69)
Физический смысл величины ДФ вполне понятен: вблизи скважин линии тока искривляются и сгущаются, так что сопротивление движению здесь оказывается большим, чем вблизи траншеи. Разность сопротивлений при исходном двухмерном и расчетном одномерном движениях и определяется эквивалентным фильтрационным сопротивлением Д Ф . Поэтому рассмотренный здесь метод сведения двухмерного движения к одномерному получил название метода эквивалентных фильтрационных сопротивлений [34]; впервые он был предложен Ю.П.Борисо- вым в 1951 г., а впоследствии развивался В.М.Шестако- вым и Ф.М.Бочевером.
Обратим внимание на то, что дополнительное сопротивление ДФ или средний напор Нф определяются лишь условиями в полосе 1*1 <а Значения ДФ или Нф не зависят от структуры потока за пределами этой полосы, и поэтому полученные выражения для Д Ф или Нф могут использоваться для задач с другими граничными условиями (если границы удалены от скважин на расстояния больше о).
ПРИМЕР. Найти решение задачи о работе ряда скважин с заданными напорами Я, между рекой и параллельным ей карьером (рис.
Решение:
[Т] заменяем ряд скважин фиктивной траншеей с напором Нф и расходом
д =
?£
2
сг
составим балансовое соотношение на линии скважин:
ще апидк. — соответственно приток к ряду из реки и отток к карье-
Р ^
РУ*
н — Н, Н(Ь-НК
а
=т—£
п
=т’_£ к
•
9Р 1 L ’ Як L ’
Рис. 3.15. Схема к задаче о контуре скважин между рекой и карьером
подставляя вместо Нф выражение (3.62) и решая систему уравнений, получим значение неизвестного расхода скважин Qc Аналогично можно использовать метод эквивалентных фильтрационных сопротивлений и для сведения несовершенных выработок к совершенным. На рис. 3.16 показана скважина, несовершенная по степени вскрытия пласта (скважина вскрывает пласт не на всю мощность). Движение к такой скважине носит двухмерный характер: линии тока вблизи скважины искривляются (фильтрация здесь не плановая), и вода испытывает при движении сопротивление большее, чем в случае совершенной скважины в том же пласте.
Решая эту задачу, М.Маскет рассмотрел скважину как совокупность точечных стоков, распределенных на участке длиной Ъ. Интенсивности (расходы) этих стоков подбирались таким образом, чтобы во всех точках линии АВ напор был одинаковым (он отвечает уровню воды в скважине). Используя далее отражение (см. раздел 3,3) стоков относительно непроницаемых кровли и подошвы пласта, М.Маскет получил решение в виде [6 ]
где Ь = Ыт Н(г)
- _ 2лТ[Н(г)-Нс}Ь ]п(4гп/гс)—Ъ]п(4гп/г)—/(1))’ (3.70)
степень вскрытия пласта; напор на расстоянии г > т;
Рис. 3.16. Схема искривления линий тока вблизи несовершенной скважины
f(b ) — функция, для которой дано графическое представле
ние.
Из решения следует, что при г~>т фильтрация к скважине носит плановый характер. Это вполне увязывается с приведенным ранее критерием (2.50).
Разобьем поток к скважине на две зоны круговым сечением г=г * т. Расход потока во внешней зоне определим по формуле (3.32; (фильтрация плановая):
_1Я Т \Н(г) - Щг0)\
Щг/m) ~ ’
а для внутренней зоны согласно формуле (3.70)
2Л'Т[Н(го) —НС]Ь
In (4m/r0) — 2>ln4 — f(B)
Помня, что Ql * Q- и исключая неизвестное значение напора Я(г0), получим запись формулы (3.70) в виде
Q‘ ln{r/7J+I
где
£ =
При
Z»
= 1 (совершенная
скважина) £
- 0, и
формула (3.71) переходит в формулу (3.32).
Следовательно, величина £ отражает
местное дополнительное сопротивление,
обусловленное несовершенством
скважины и связанным с этим отклонением
от плановой модели фильтрации. Отсюда
следует возможность замены несовершенной
скважины фиктивной совершенной, имеющей
расчетный радиус
г*
(3.72)
Из этого примера ясна принципиальная роль метода эквивалентных фильтрационных сопротивлений: он позволяет учесть локальные погрешности плановой модели фильтрации (см. раздел 2.5) и тем самым резко расширить диапазон ее практического применения. Такой подход оказывается особенно эффективным при моделировании задач в плановой постановке. (Примеры приложений этого метода приведены также в разделах 8.2 и 8.3).