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

Размещено на http://www.allbest.ru/

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

В.Ф. Жабин, кандидат геолого-минералогических наук, доцент

ФГОУ ВПО «Московский государственный университет

природообустройства», г. Москва, Россия

Н.П. Карпенко, доктор технических наук, доцент

ГНУ Всероссийский научно-исследовательский институт

гидротехники и мелиорации

им. А.Н.Костякова Россельхозакадемии, г. Москва, Россия

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

The technique of the account of vertical anisotropy in filtrational calculations of horizontal absorbing wells and interpretation of experimental evacuations from incomplete driven wells is offered. It is noted, that anisotropy gains special value at assessments of works of incomplete absorbing wells and driven wells when power of a reacting region of their act is sized up.

Введение

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

Отложения взвесенесущих потоков, объединенных в группу флювиальных (аллювиальные, морские, пролювиальные, озерные и водно-ледниковые), а также их сочетания занимают ведущее место в качестве элементов природной среды инженерных сооружений и систем. Установлено, что в обычных условиях формирования флювиальных отложений образуются слои с выраженной вертикальной микронеоднородностью, обусловленной размером частиц, их формой, упаковкой и ориентировкой [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 - Расчетные формулы для определения сопротивления горизонтальных дрен

В формуле (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 - Результаты расчета междренных расстояний в зависимости от типа природной среды

Расчет междренных расстояний выполнен по формулам С.Ф. Аверьянова и В.М. Шестакова, преобразованным нами [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 с.

Размещено на Allbest.ru