Учет вертикальной анизотропии в фильтрационных расчетах горизонтальных дрен и несовершенных скважин
Введение
Урбанизация, хозяйственное использование территорий и отсутствие гидрогеологического мониторинга остро ставят проблему соблюдения требований охраны окружающей среды и разработки защитных мероприятий по предотвращению негативных процессов освоения земель на стадии проектирования. В этих условиях прогноз, связанный с оценкой режима грунтовых вод в перспективе, параметров дренажа и его эффективности становится одним из главных элементов проектирования, основанного на результатах гидрогеологических исследований. Особое внимание здесь должно уделяться изучению неоднородности отложений, выраженной в изменении основного их свойства, в данном случае, водопроницаемости.
Отложения взвесенесущих потоков, объединенных в группу флювиальных (аллювиальные, морские, пролювиальные, озерные и водно-ледниковые), а также их сочетания занимают ведущее место в качестве элементов природной среды инженерных сооружений и систем. Установлено, что в обычных условиях формирования флювиальных отложений образуются слои с выраженной вертикальной микронеоднородностью, обусловленной размером частиц, их формой, упаковкой и ориентировкой [5]. Эта особенность петрогенеза находит отражение в различии водопроницаемости слоя горных пород в зависимости от направления фильтрации-анизотропии. Важность наблюдаемой при бурении скважин неоднородности (макро) неоднократно подчеркивается многими исследователями, однако в задачах подпора грунтовых вод и дренирования, серьезного внимания требует и внутрислоевая неоднородность, редко определяемая визуально (открытые горные выработки и естественные обнажения).
Рассмотренное свойство упомянутых отложений, является носителем фильтрационной анизотропии, методы оценки и учета которой в расчетах разработаны недостаточно полно. Ранее нами освещались вопросы учета анизотропии в гидрогеологических расчетах [2, 4, 5]. Недооценка ее влияния, точнее, неправомерное допущение изотропности природных сред при математическом описании процессов инфильтрационного подпора и дренирования, приводит к занижению функции напора в характерных сечениях фильтрационных потоков или к ошибочным расчетам дренажа, поэтому при решении задач гидрогеологического прогнозирования эти характеристики являются определяющими.
Методы исследований
Вертикальный анизотропия горизонтальный дрена
Переход от слоисто-неоднородной среды к анизотропной осуществляется при допущении (Пашковский И. С.). Если учесть, что мощность элементарного слоя (mi) измеряется долями метра, то легко убедиться в почти априорном соблюдении данного условия.
Рассмотрим некоторые стороны методов расчета горизонтального дренажа в анизотропной среде, схема которой представлена на рис. 1.
Рис - 1 Симметричный поток к несовершенной дрене
В общем случае движения симметричного двухстороннего инфильтрационного потока, длиной, для любой точки профиля от водораздела до дрены, справедливо считать [2]:
(1)
При 3 ,
Tх = 0,5kx(ml + md), m1 = md+H
Посредине между дренами при имеем:
(2)
При, (3)
Где: m, d - мощность водоносного горизонта под дреной.
(4)
Значения функций () и () приведены в работе [2]. Исследование характера изменения приведенных специальных функций показывает линейный подъем уровня грунтовых вод в начальный период инфильтрации, поэтому влияние дрены не сказывается. Затем скорость повышения поверхности грунтовых вод снижается по мере усиления дренирования. При = 1 величина инфильтрационного подпора достигает 0,9 от максимума. Дренажный расход также увеличивается во времени. Эффективность естественных дрен, время ввода искусственного дренажа и его расход ориентировочно можно оценить, используя приведенные зависимости и допустимую величину подъема уровня грунтовых вод.
Отдельно следует остановиться здесь на параметре, характеризующем дополнительное фильтрационное сопротивление на границе потока грунтовых вод. Этот параметр отражает несовершенство дрены (скважины) по отношению к дренируемому водоносному пласту. Принято выделять два вида несовершенства:
- - по степени вскрытия пласта, соответствующей положению дрены (скважины) относительно его вертикальных границ и фильтрационной неоднородности или анизотропии; - по характеру вскрытия, зависящего от заиления стенок горной выработки, применяемых материалов и технологии строительства. Несовершенство второго вида иногда называют "скин-эффектом".
Величину сопротивления принято выражать в метрах эквивалентной дополнительной длины потока, полученной расчетным путем или графическим построением (виртуальным смещением реальной границы). Аналитические способы определения сопротивления горизонтальных дрен приведены в таблице 1.
Таблица 1 - Расчетные формулы для определения сопротивления горизонтальных дрен
Расчетная схема |
Строение пласта |
Расчетная формула |
Водоем или водоток шириной 2br с изм. |
Однородно-анизотропный Двухслойный |
(5) (5а) |
Горизонтальная дрена |
Однородно-анизотропный |
(6) |
Водоток, водоем, ограниченной ширины |
Двухслойный |
(7) |
Систематический горизонтальный дренаж |
Двухслойный, анизотропный |
(8) |
Систематический горизонтальный дренаж* |
Трехслойный |
(9) |
Комбинированный дренаж* |
Трехслойный, двухслойный |
(10) |
Береговая дрена (по С. Ф. Аверьянову, с изменениями) |
Однородно-анизотропный |
(11) |
В формуле (11) возможна замена Lnd на L., перед md добавляется множитель, учитывающий фильтрационную анизотропию, по аналогии с зависимостью (6), а значение синуса можно заменить значением его аргумента. Данные таблицы 1 показывают, что практически все дренирующие границы потоков грунтовых вод имеют дополнительное фильтрационное сопротивление, которое необходимо учитывать заменой действительной координаты (х) на расчетную координату Х = х + ().
Заметим, что в формулы (5) и (7) входит половина ширины естественной дрены, которая может быть определена, как и, соответственно, для однородно-анизотропного и двухслойного строения ложа водоема или водотока.
Единичный установившийся расход (удельная величина естественных ресурсов в условиях естественной дренированности) оценивается, как:
. (12)
Однако, расход воды, поступающий в дрену, распределяется в зависимости от ее несовершенства и условий формирования потока грунтовых вод. Так единичный расход инфильтрационного потока грунтовых вод в реке, следуя рекомендациям С. Ф. Аверьянова, определим, как:
Qr = q0 , , (13)
Что при = L, (хd = 0) дает и qr 0,67q0.
При этом на урезе разгружается qro= 0,33q0. Иными словами, независимо от величины сопротивления, на урезе водотока (водоема) разгружается около 1/3 расхода потока грунтовых вод, идущего с водораздела. Близкий по значению результат получается, если использовать решение В. М Шестакова с заменой х на Х. В этом случае, расход грунтовых вод в каждом сечении под рекой будет равен:
, (14)
А на урезе х = 0 и qro = q0 e-1 0,368q0.
Тот же результат получается при решении дифференциального уравнения вида:
(15)
Тогда
(16)
Определим ширину зоны разгрузки, считая от уреза. При х = 2L, qr = 0,05q0, то есть на этом расстоянии от уреза поток грунтовых вод разгрузится только на 95%. При меньшей ширине реки, "проскок" потока грунтовых вод под рекой будет больше, поэтому это необходимо учитывать при размещении головных дрен. Аналогичные оценки выполним на основе решения (3) для систематического дренажа. В этом случае, в дрену поступит половина единичного притока:
(17)
Глее: , - уклон потока грунтовых вод.
Если дрена является совершенной, то, формула (17) обращается в аналог формулы Кене:
(18)
При расчете междренных расстояний в анизотропных средах необходимо определять мощность "активной зоны", в которой проходит основной фильтрационный поток под дреной.
При L/m 2,55 имеем [1]:
, . (19)
Прямое определение На по формуле (19) весьма затруднительно, так как искомая функция содержится в аргументах. Расчет можно упростить, если подбором найти редукционные коэффициенты (б1 , б2) <1, используя тождества:
Для однородной среды
(20)
,
,
Для двухслойной
; (20а)
, , ,
.
Для нестационарного потока в пространстве между дренами существенную роль при аналогичных оценках играет фактор времени. Кроме того, надо корректировать величину инфильтрационного питания, если норма осушения меньше "критической" глубины испарения грунтовых вод [5].
В таблице 2 приведены для сравнения результаты расчета междренных расстояний (L) в зависимости от типа среды и фильтрационной схемы.
Таблица 2 - Результаты расчета междренных расстояний в зависимости от типа природной среды
Однослойная Изотропная среда |
Md, м |
ДН, м |
Рd, м |
К, L |
К, L |
К, L |
Щ, м/сут |
10 |
1,5 |
1,57 |
|
|
|
4-10-4 | |
Однослойная Анизотропная Среда |
Mdа, |
ДН |
Рd |
Кz |
Кх |
L |
Щ |
7,3 Б = 0,73 |
1,5 |
1,57 |
0,03 |
3 |
542 |
4-10-4 | |
Двухслойная Изотропная среда |
Md2,md1 |
ДНг |
Рd |
К1 |
К2 |
L |
Щ |
8 4 |
1,43 |
1,57 |
0,03 |
12 |
111 |
4-10-4 | |
Двухслойная Анизотропная Среда |
Mа2,md1 |
ДНг |
Рd |
К1z, к1х |
К2h |
L |
Щ |
3,9 4 Б = 0,49 |
1,43 |
1,57 |
0,3 3 |
12 |
510 |
4-10-4 |
Расчет междренных расстояний выполнен по формулам С. Ф. Аверьянова и В. М. Шестакова, преобразованным нами [5] для анизотропных сред. В таблице 2 показано влияние анизотропии на величину расстояния между дренами и мощность активной зоны, несмотря на формальное равенство отдельных значений коэффициентов фильтрации (м/сутки). Отметим, что для однородно-изотропной среды К = 0,3 м/сут соответствует среднегеометрическому из двух остальных величин (3,0 и 0,03). Дополнительный анализ показал, что междренное расстояние практически не зависит от общей мощности анизотропного слоя, если она больше активной зоны.
Приведенные результаты показывают необходимость выбора гидрогеологических характеристик для расчета параметров дренажа на основе литолого-генетических и фильтрационных критериев. Это, в свою очередь, требует целенаправленного (не формального) подхода при постановке, проведении и интерпретации опытно-фильтрационных работ (ОФР) или режимных наблюдений (ОФН). Методы фильтрационного опробования анизотропных природных водонасыщенных сред и зоны аэрации при изысканиях доступно изложены в существующих литературных источниках, включая научные разработки и пособия [3, 5, 6].
Подчеркнем важность определения горизонтальной водопроницаемости покровного слоя двухслойного пласт по данным экспресс-откачек из скважин, не оборудованных фильтрами. Вместе с тем, можно использовать данные наблюдений за снижением динамического уровня в скважинах покровного слоя при откачке из нижнего, хорошо проницаемого слоя. При постоянном дебите откачки, водопроводимость покровного слоя (Тх) приближенно определяется в период квазистационарного режима по формуле Дюпюи для двух наблюдательных скважин, расположенных в покровном слое.
В практике проектно-изыскательских работ используется только вертикальная составляющая водопроницаемости этого слоя, что приводит к неоправданному занижению параметров дренажа и искусственному завышению его стоимости. Очевидной является роль идентификации расчетной схемы дренажа и геологических условий, которая способствует достижению оптимальной стоимости и функциональной эффективности проектируемых мероприятий. Учет "активной зоны" при проектировании работы несовершенных скважин вертикального дренажа рекомендуется проводить с учетом рекомендаций [2].
В анизотропном водоносном слое, согласно исследованиям Ю. О.Зеегофера, имеем:
при r ? 0,25 (21)
Имея данные двух наблюдательных скважин, можно найти фильтрационные характеристики, а по результатам 2-3 ступеней откачки уточнить зависимость величины активной зоны от понижения. Наиболее полно определение гидрогеологических параметров в разнообразных условиях представлено В. А. Мироненко и В. М.Шестаковым [6]. Ориентировочную оценку фильтрационных характеристик по результатам кустовой откачки можно выполнить при стабилизации режима опыта, который наступает довольно быстро, особенно, в анизотропной среде, когда расстояния до наблюдательных скважин составляют:
R ? На, (22)
Где: - фактор стабилизации, Та - обобщенная водопроводимость активной зоны.
В однородно-анизотропном слое д = 3, трещинно-пористой анизотропной среде этот параметр уменьшается до 1,5. Зависимости Q > S могут иметь различную форму, которая в общем случае, определяется вкладами нелинейности фильтрации, гидравлическими потерями в водоподъемных трубах, скинэффектом и ростом активной водопроводимости [6]. Коэффициент вертикальной водопроницаемости находится по формуле (22), а расчетная мощность по вертикали вычисляется путем раскрытия фактора стабилизации.
Заключение
Предложенная методика учета вертикальной анизотропии в фильтрационных расчетах горизонтальных дрен и интерпретации опытных откачек наиболее полно идентифицирует природную среду при решении гидрогеологических задач, связанных с осушением территорий и оценкой эксплуатационных запасов подземных вод. При выполнении проектов дренажа необходим регламент планирования изысканий с целью обоснования расчетных схем по литогенетическим и фильтрационным признакам, а также изменение технологии расчетов параметров систематического горизонтального дренажа и оценки работы несовершенных скважин.
Библиографический список
- 1. Аверьянов С. Ф. Борьба с засолением орошаемых земель. - М.: Колос, 1978, 288 с. 2. Аравин В. И., Нумеров С. Н. Фильтрационные расчеты гидротехнических сооружений. - М.: Стройиздат, 1948, 226 с. 3. Богомолов Ю. Г., Жабин В. Ф., Хачатурьян В. Х. Изменение гидрогеологических условий под влиянием мелиорации. - М.: Наука, 1979, 162 с. 4. Лютин Дж. Дренаж сельскохозяйственных земель. Пер. с англ. под ред. С. Ф.Аверьянова. - М.: Колос, 1964, 720 с. 5. Манукьян Д. А., Жабин В. Ф. Гидрогеоэкологические проблемы в задачах природообустройства. - М.: МГУП, 2006, 194 с. 6. Мироненко В. А., Шестаков В. М. Теория и методы интерпретации опытно-фильтрационных работ. - М.: Недра, 1978, 325 с.
Похожие статьи
-
Расчет расстояния между дренами - Мелиорация участка переувлажненных земель
Междренное расстояние (В, м) является наиболее важным параметром дренажа, определяющим эффективность работы осушительной сети. Величина В зависит от...
-
Таблица 4. Расчеты надежности параллельной схемы выполнения системы "ПЦР-РПИФ" с учетом превалентности случаев АЧС в фермах разного размера (расчеты...
-
Водный режим почв характеризуется изменением во времени влажности в корнеобитаемом слое, которая определяется атмосферными осадками, поверхностным...
-
Предприятия, занимающиеся заготовкой и хранением зерновых ресурсов, принимают от хлебосдатчиков партии зерна. Каждую партию зерна вместе с примесями...
-
Расчет доз удобрений на планируемую урожайность - Пшеница
Пшеница зерновой агропромышленный почвенный Существуют различные методы определения оптимальных доз удобрений под планируемую урожайность. В курсовой...
-
Расчет теплового баланса - Составление здания птичника на 22 тысячи голов цыплят бройлеров
Тепловой баланс - это соотношение прихода и расхода тепла в животноводческом помещении. Тепловой баланс бывает нулевым, отрицательным и положительным....
-
Расчет искусственного освещения - Проектирование малой фермы в ТОО "Красная звезда"
При расчете искусственного освещения используется метод светового потока [3]: , лм (2.40) Где - световой поток принятой лампы накаливания, лм; - площадь...
-
Государственная программа устойчивого развития села на 2011-2015 гг. в отрасли свиноводства предусматривает к 2015 г. получать среднесуточный привес...
-
Помет является побочным продуктом в птицеводстве. Это ценное органическое удобрение. В то же время помет в птичнике является фактором загрязнения среды и...
-
Пониженная температура воздуха в помещениях для животных в сочетании с высокой влажностью и повышенной подвижностью его даже при вполне...
-
Силовым электрооборудованием зерноочистительно-сушильного комплекса является электродвигатели, которые используются в приводах технологического...
-
Составляем таблицу отклонений и потерь в линиях напряжением 10-0.38 кВ. системы электроснабжения, если норма отклонений напряжения у потребителей Vнор= ....
-
Расчет времени тушения участка кромки крупного лесного пожара рабочими с ручными инструментами При составлении оперативного плана тушения крупного...
-
Выбор распределительных устройств Распределительные устройства предназначены для приема и распределения электрической энергии внутри помещения. С этой...
-
Расчет общей площади питомника - Проект лесного питомника
Исходными данными для расчета площади питомника является ежегодная потребность в посадочном материале. Основой для расчета площади питомника являются...
-
Расчет действительно возможной урожайности, Обоснование программируемой урожайности - Пшеница
Действительно возможную урожайность (ДВУ) рекомендуется рассчитывать по баллам бонитировки почв и цене балла по культуре. При расчете ДВУ рекомендуется...
-
Открытая проводящая и оградительная сеть Задачами гидравлического расчета открытой проводящей и оградительной сети являются: а) проверка и уточнение...
-
Учет племенного молодняка - Особенности разведения, кормления и разделывания кроликов
Для получения племенного молодняка кроликов разных пород прежде всего районным, областным и краевым обществам кролиководов - любителей следует...
-
Расчет потребности в машинах - Механизация малоценных молодняков коридорным способом
Количество рабочих машин, необходимое для выполнения установленного объема работ в определенные агротехнические сроки, зависит от объема этих работ и...
-
Для начала нам надо найти часовой объем вентиляции (L) - количество свежего воздуха, которое необходимо подать в помещение в течении 1 ч(м3/ч). Мы будем...
-
Структура урожая показывает Из каких элементов складывается его величина, при какой доле их участие формируется запланированный урожай. Основные...
-
Наиболее интенсивной формой организации ПЛСБ является закладка лесосеменных плантаций вегетативным или семенным потомством плюсовых деревьев. Ценные...
-
Зона хранения навоза. Площадь хранилища определяется по формуле: Где: - суточный выход кала от 1 гол., кг - суточный выход мочи от 1 гол., кг - норма...
-
Расчет технологии Примерное количество Qcyт (кг) навоза, получаемое на комплексе за сутки, может быть определено по формуле Годовой выход Qгод навоза...
-
Общая площадь приточных каналов берется 20-70 % от общей площади всех вытяжных труб, а для Севера-Запада (Ленинградской области) она составляет 50%, т. е...
-
При подсчете транспортных затрат, а также при оценке компактности хозяйства и условий управления вычисляют средние расстояния от усадьбы до массивов...
-
Площадь рассчитывается по хозяйственным отделениям отдельно по каждой породе по формуле: Р = N*a*S / n*S1 (2.1) Где Р - площадь, необходимая для...
-
Как мы говорили, резание в корнерезках осуществляется резцом, имеющим форму клина. Теорию резания резцом сельскохозяйственных материалов (почва,...
-
Основы теории и расчет сепараторов - Механизация первичной обработки молока
Основы теории и расчет сепараторов разработали советские ученые Г. А. Кук, Г. И. Бремер, Н. Я. Лукьянов и др. Физическая сущность процесса сепарирования...
-
Расчет регенеративных теплообменников - Механизация первичной обработки молока
После пастеризации молоко охлаждают до температуры, обеспечивающей его сохранность (5...80С). Чтобы снизить расход холода, на первой ступени его...
-
Исходя из поголовья согласно заданию, рассчитаем структуру поголовья фермы (таблица 1.2.1) . Таблица 1.2.1 Структура поголовья фермы КРС, % Группы...
-
Расчет пастеризаторов - Механизация первичной обработки молока
Механизация первичный обработка молоко Основы теории пастеризации молока разработал советский ученый Г. А. Кук. Теория основывается на том, что гибельно...
-
На заводе применяют сплошной объемный метод учета молоди. Учет количества выращенной молоди проводят в рыбоуловителе, который изготовлен из металлической...
-
Обоснование необходимости разработки При хранении сортового картофеля необходимо более строго соблюдать технологический режим, от этого зависит качество...
-
Таблица 1 Размеры положительных АЧС выборок (количества страт и животных по видам), использованные в расчетах для определения диагностической...
-
Следует различать Виды Арендной платы как способа начисления платежей и ее Формы, Т. е. внешнее материально-вещественное или стоимостное их воплощение. В...
-
Расчет количества семей на перспективу. - Проект планировки и застройки сельского населенного пункта
На основании анализа сложившегося размещения населенных пунктов, численности проживающего в них населения (постоянно и временно), анализа демографических...
-
Исходными данными для выполнения расчетов является производительность Q (т/час), масса корма, размещенная на 1 м2 поверхности вальца G, Энергоемкость Е....
-
Основным расчетом в цехе первичной переработки является расчет длины подвесного пути, который может быть конвейерным, бесконвейерным или частично...
-
Стоимость отдельного кадастрового квартала определяется по формуле: КСI = БПВ * КЦ * Q, (6) Где: КСI - кадастровая стоимость; БПВ - базовая постоянная...
Учет вертикальной анизотропии в фильтрационных расчетах горизонтальных дрен и несовершенных скважин