25305 (686796), страница 2
Текст из файла (страница 2)
Для моделирования фильтрации в основном используются гидравлическая и электрическая аналогии, реализуемые на сплошных и сеточных моделях.
Сплошные гидравлические модели, представленные фильтрационными лотками различных видов, в гидрогеологических расчетах применяются редко.
В развитии методов моделирования фильтрации подземных вод основная роль принадлежит сплошным и сеточным электрическим моделям, основанным на использовании метода электрогидродинамических аналогий (ЭГДА), сущность которого наглядно представляется сопоставлением основных законов движения фильтрационного потока и электрического тока:
закон Дарси и закон Ома
и
(9)
где Q — расход; F — площадь поперечного сечения потока; Н — напор; х-—расстояние;I— сила тока; с — удельная проводимость, ; р — удельное сопротивление; площадь поперечного сечения проводника; U — электрический потенциал, l — длина проводника.
Приведенная формула закона Ома получена путем несложных преобразований
;
(10)
где R — сопротивление.
Идентичность записи законов Дарси и Ома очевидна. В них соответствуют физические характеристики — коэффициент фильтрации Кф и удельная проводимость с (физическое подобие), силовые характеристики — напор Н и потенциал U(динамическое подобие) и, наконец, расход потока Q и сила тока (кинематическое подобие).
На сплошных моделях ЭГДА фильтрационный поток моделируется сплошным электрическим полем, геометрически подобным. Для этого применяются электропроводная бумага и электролиты. Электропроводная бумага изготавливается с удельным сопротивлением от 100 до 100 000 Ом/см, в зависимости от количества содержащихся в ней сажи и графита.
Участки поля с различной проницаемостью пород моделируются кусками бумаги различной удельной проводимости. Между собой участки модели скрепляются специальным электропроводным клеем.
Электролиты также широко используются в качестве материала модели и обычно представляют собой растворы солей, причем наибольшее распространение получили водные растворы поваренной соли и медного купороса. Кроме того, можно использовать электропроводные краски, клеи, электропроводный картон, гипс и т. д.
Определение приведенного потенциала на моделях ЭГДА производится с помощью мостовой измерительной схемы.
При составлении сеточных моделей поток разбивается на отдельные блоки, центры которых связываются электрическими резисторами. В таких моделях геометрическое подобие модели и объекта не сохраняется.
4. Приток воды к водозаборным сооружениям
Среди водозаборных сооружений мы будем рассматривать такие горные выработки, как дрены (канавы) и скважины. В гидрогеологии горные выработки разделяют на совершенные и несовершенные.
Гидродинамически совершенной называется горная выработка, вскрывающая водоносный горизонт от кровли до подошвы.
Рис. 1. Схемы совершенной (а) и несовершенной (б, в) выработок
Приток воды к ней происходит по всей поверхности соприкосновения стенок выработки с водоносным горизонтом (рис. 1, а). Если же выработка не доходит до водоупора, она называется несовершенной по степени вскрытия водоносного горизонта (рис. 1, б). Зачастую выработки закрепляются от обрушения, цементируются скважины оборудуются обсадными трубами, фильтрами и т. п. Естественно, что приток воды в такие выработки затруднен и их называют несовершенными по характеру вскрытия водоносного горизонта.Основные уравнения притока воды к водозаборам (скважинам и дренам) будем выводить при условии совершенства выработок.
Представим себе плоский поток грунтовых вод. Гидравлический градиент I в данном случае равен
(11)
где х — расстояние между сечениями h1 и h2
Если мы будем сближать сечения h1 и h2 так, чтобы расстояние между ними стало равно нулю, то получим уклон (гидравлический градиент) в точке а, который равен тангенсу угла наклона зеркала грунтовых вод или первой производной
(12)
Подставив полученное выражение гидравлического уклона в выражение закона Дарси (21), получим для безнапорных вод
(13)
для напорных вод
(14)
где H —напор, отсчитываемый от подошвы водоносного пласта до его пьезометрического уровня.
4.1 Приток безнапорных вод в совершенную горизонтальную дрену (канаву)
После устройства дрены скорость движения воды в ней увеличивается и уровень воды понижается на величину S, которую в гидрогеологии принято называть величиной понижения. Иными словами, величина понижения представляет собой разницу между статическим и динамическим уровнями. Мощность водоносного горизонта до понижения обозначим через H, глубину воды в дрене — через ho. В результате понижения уровня в дрене в водоносном горизонте образуется депрессионная воронка, показанная на рис. 50 сплошной жирной линией. Расстояние R, на которое сказывается влияние понижения, называют радиусом влияния.
Для расчета притока воды в дрену Q выбираем на расстоянии х от стенки дрены сечение с напором к, которое находится в интервале от нуля до R.
В общем виде приток воды в дрену будет равен выражению (13). Подставим сюда величину площади фильтрации
(15)
где В — длина дрены. Получим
(16)
При расчете притока воды в дрену удобно пользоваться понятием единичного притока д, т. е. притока воды на единицу длины дрены
(17)
Отсюда элементарная формула для расчета притока воды
(18)
Разделим переменные в выражении (18), т. е. умножим обе его части на dх и проинтегрируем
(19)
В результате получим
(20)
(21)
(22)
Формула (22) выражает величину единичного притока с одной стороны дрены. Для получения полного притока воды в дрену необходимо умножить единичный приток на два, а затем — на длину дрены. Приток воды в торцы дрены обычно не учитывают, так как он при большой длине дрены составляет ничтожную долю.
По формуле (22) можно рассчитать расход плоского грунтового потока. Подставив вместо радиуса влияния расстояние между сечениями, равное I, получим
(23)
Выражение можно записать так
(24)
т. е. единичный расход равен
(25)
а полный расход составит
(26)
Исследуя выражение (22), мы сможем решить одну из весьма важных задач в гидрогеологических расчетах — вывести уравнение депрессионной кривой. Построение депрессионной кривой необходимо при возникновении угрозы затопления подземными водами котлованов, подвалов зданий и т. п..
Изменив пределы интегрирования в выражении (20) по X от 0 до х, а по У от h0 до h К получим
(27)
Естественно, что приток воды в выражениях (22) и (27) одинаков, т. е.
(28)
Решаем (28) относительно h
(29)
Для построения депрессионной кривой мы задаемся величиной hо в зависимости от 5, мощность водоносного горизонта H легко получить по данным бурения, величину радиуса влияния можно найти по эмпирическим формулам (об этом мы поговорим позже).
На миллиметровой бумаге строим разрез через дрену и котлован (рис. 52) и, задаваясь разными значениями х(хи x2,..., хп), например 10, 20, 30 и т. д. метров, получаем величины h(h, h2,..., hп). Соединив полученные точки плавной линией, получим кривую депрессии. Если она проходит через котлован, строят новую кривую, задавшись большей величиной понижения и, естественно, меньшим значением глубины воды в дрене. Построение производят до тех пор, пока депрессионная кривая не опустится ниже дна котлована.
4.2 Расчет притока грунтовых вод в скважину
Здесь мы имеем дело не с плоским потоком, как в предыдущем примере, а с радиальным. На рис. 10 показаны все обозначения, которые нам ясны из предыдущей задачи, кроме г — радиуса скважины.
Расчет начинаем с уравнения (13). Площадь притока воды равна площади боковой поверхности цилиндра, радиус которого равен х, т. е.
(30)
(31)
Разделяем переменные (т. е. умножаем обе части уравнения на ) и интегрируем
(32)
Интегрирование по х производим не от нуля, а от r — стенки скважины, в результате получаем формулу Дюпюи
(33)
(34)
Рис. 2. Схема притока грунтовых вод в совершенную скважину
В таком виде использовать выражение (2) не совсем удобно, так как в нем присутствует натуральный логарифм. Подставим вместо него десятичный (1nх = 2,30 lgх), а вместо л его значение и получим более удобное выражение для расчета притока безнапорных вод в скважину:
(35)
Выражение можно видоизменить:
(36)
И, подставив его в зависимость 54), получим
(37)
Для построения депрессионной кривой возвращаемся к (35) и изменяем пределы интегрирования: по X от r до х а по Y от h0 до h:
(38)
(39)
(40)
Решаем равенство относительно h и получаем уравнение кривой депрессии
(41)
4.3 Расчет притока напорных вод в совершенную дрену
Площадь фильтрации в сечении h, расположенном на расстоянии х от стенки дрены, будет равна
(42)
Мы здесь снова не учитываем приток воды через торцы дрены. Подставляем площадь в (43):
(43)
Переходим к единичному расходу
(44)
Разделяя переменные и интегрируя, получим
(45)
(46)
(47)
Выражение (47) представляет собой единичный приток артезианских вод в один из бортов канавы. Полный приток составит
(48)
Если нам необходимо получить уравнение депрессионной кривой, то (43) нужно проинтегрировать по Х от 0 до ж, а по Y от hо до h: