Определение эффективности вертикального дренажа на основе натурных исследований и математического моделирования

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

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

Для определения мероприятий по предупреждению подтопления на юге Украины актуальным является моделирование уровня грунтовых вод территорий по данным натурных наблюдений для конкретных объектов. Для этого был выбран объект в Каланчакском районе, который находится между Северо-Крымским каналом и р. Каланчак (рис. 1).

По данным Каховской ГГМЭ, наибольшие площади подтопления находятся в этом районе, где причинами подтопления земель является напор плиоценового водоносного горизонта, фильтрация оросительной воды из магистральных каналов, осадки и поливы приусадебных участков и сельскохозяйственных угодий [1].

Рис. 1. Цифровая модель рельефа (ЦМР) объекта исследований (пгт. Каланчак и прилегающих территорий)

В последнее время усилилось напорное питание грунтовых вод за счет плиоценового водоносного горизонта [1], что имеет место и в пгт Каланчак. Это подтверждает необходимость изучения гидрогеологических процессов на основе натурных исследований и математического моделирования с использованием ГИС технологий.

Характеристика объекта моделирования. На основе изучения геологических разрезов и конструкций скважин вертикального дренажа (рис. 2) была определена физическая модель, характеристика толщи отложений, к которой приурочен водоносный горизонт.

Рис. 2. Геологический разрез через ряд скважин вертикального дренажа в пгт. Каланчак

Физическая модель показывает, что толща отложений 6…19 м к водоносному горизонту состоит в основном из суглинков, лесовидных суглинков, супесчаных глин с небольшими коэффициентами фильтрации. Территория подтопляется в основном за счет давления, которое образуется в напорном водоносном горизонте. Под действием потенциала этого давления шp, а также капиллярно-сорбционного потенциала шm, который создает капиллярную кайму, вода поднимается в толще суглинков с напорного горизонта.

Таким образом, полный потенциал давления в толще к водоносному горизонту состоит из компонентов:

Ф = шg + шm + шp, (1)

где шg - гравитационный потенциал; шm - капиллярно-сорбционный потенциал; шp - потенциал давления (выше уровня грунтовых вод равен нулю), который действует против силы земного притяжения в связи с напорным характером водоносного горизонта.

Математическая модель фильтрации состоит в таком случае из системы уравнений и учитывает процессы насыщенно-ненасыщенного влагопереноса в слое почвы к напорному горизонту и процессы изменения напора в водоносном горизонте в результате откачки вод из дренажной скважины.

Как известно, движение воды в насыщенно-ненасыщенных грунтах осуществляется согласно закону Дарси:

(2)

где Ф = шg - шm - шp, при условии, что ось z направлена вниз; k(z,и) - коэффициент влагопереноса или гидравлическая проводимость почвы, которая является функцией объемного содержания влаги в почве и и свойствами почвы в точке z.

Используя уравнение Дарси и уравнение неразрывности для несжимаемой жидкости, при определенном выборе единиц измерения, то есть при Ф = z - шg - шm - шр, получим основное уравнение движения почвенной влаги:

. (3)

Данное уравнение используется для описания движения воды в неоднородной среде, если k = k(z,), ш = ш(z,).

При неустановившейся откачке в напорном горизонте имеет место уравнение:

(4)

при начальных и граничных условиях:

(5)

Система уравнений (3)…(5) описывает качественные и количественные процессы работы вертикального дренажа.

Идентификация параметров модели понижения уровней по данным откачек из скважин в водоносных горизонтах. При неустановившейся фильтрации величина понижения S является функцией времени t в каждой точке радиуса r. Для моделирования понижения воспользуемся формулой Тейса-Маскета:

. (6)

Величина і рассчитывается формулой:

, (7)

тангенс угла наклона графика временного прослеживания уровня Srt =f(lnt) к оси ln t в точке r; a - коэффициент пьезопроводимости; Srt - понижение уровня в точке r на момент времени t при данном дебите скважины Q = const.

Неизвестными параметрами в этих формулах является проводимость Т и коэффициент пьезопроводимости а, идентификация которых осуществляется по данным производственных откачек. Величины Srt являются экспериментальными значениями функций, которые используются при определении неизвестных коэффициентов.

Величина a в точке r рассчитывается по формуле

(8)

с использованием понижений Srt в наблюдательных скважинах r. Соответственно, Т - проводимость пласта (за Тейсом), которая определяется за изменением уровня в наблюдательных скважинах

. (9)

По данным производственных откачек, в наблюдательной скважине, которая находится на расстоянии r = 280 м от дренажной скважины, получена зависимость понижений (рис. 3).

Для получения параметров Т и а использовались аппроксимированные данные понижений Srt в моменты времени t = t1, 2t1, 3t1, 4t1, 5t1 , кратные t1,= 10 сут. Рассчитанные в разные моменты времени параметры Т и а приведены в табл.1.

ГИС для моделирования подтопления территорий по данным натурных наблюдений. Для построения карты глубины залегания грунтовых вод нами предлагается между наблюдательными скважинами создавать триангуляционную сеть и выполнять расчеты между ними с помощью формул гидрогеологического расчета [3, 4].

Рис. 3. Аппроксимация значений понижений в наблюдательной скважине пгт. Каланчак

Таблица 1. Результаты расчетов фильтрационных параметров Т і а

водоносный скважина аппроксимация

В середине каждого треугольника, между наблюдательными скважинами, расчет производится с использованием метода сеток решения первой задачи Дирихле для уравнения Лапласа [3].

Так же, с помощью метода сеток, строится цифровая модель рельефа (ЦМР), для её создания используются горизонтали рельефа, которые являются граничными значениями для расчетов. Наблюдательные скважины в этом случае наносятся на ЦМР в виде точек с обозначением абсолютных отметок поверхности земли каждой скважины.

Таким образом, обобщенный алгоритм метода построения карт на базе специализированной ГИС состоит из следующих этапов:

1. Проверяется гипотеза квазистационарного движения (изменения уровня грунтовых вод), исходя из гидрогеологических условий и данных измерений в наблюдательных скважинах.

2. На базе планового расположения и технических (паспортных) данных наблюдательных скважин, которые входят в постоянно действующую наблюдательную сеть гидрогеолого-мелиоративной экспедиции для измерений уровня грунтовых вод, строится триангуляционная сеть данной территории.

3. Между отдельными точками треугольника, который находится в триангуляционной сети, строятся кривые депрессий по формулам [3, 4].

4. Для области треугольника решается первая задача Дирихле для уравнения Лапласа, в которой граничными условиями первого рода будут поверхности натяжения на сторонах треугольника.

5. Строятся поверхности натяжения, по стандартной методике проводятся гидроизогипсы и строится карта уровня грунтовых вод.

С использованием данного принципа построена карта глубины залегания первого от поверхности водоносного горизонту исследуемого объекта (рис. 4).

Рис. 4. Карта глубин залегания первого от поверхности водоносного горизонта в пгт. Каланчак и прилегающие территории состоянием на 16.07.2008

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

С этой целью выбран объект (рис. 4), для которого выполняется стационарный (квазистационарный) режим уровней грунтовых вод и описывается уравнением Лапласа. Для данного объекта были проведены замеры уровня грунтовых вод в наблюдательных скважинах (скв. 139а, 141) и дополнительных буровых скважинах (бур. 1, 2, 3). Кроме того, установлены гидрогеологические характеристики объекта (коэффициенты фильтрации, отметка водоупора).

В результате расчетов получена поверхность натяжения залегания уровня грунтовых вод (рис. 5а). Поверхности аппроксимации, которые получены за этими методами, сравниваются с фактическими замерами в четырёх точках (бур. 2, скв. 141, скв. 139а, бур. 3). Точность обоих методов исследована по модулю средней абсолютной и средней относительной ошибками.

Установлено, что моделирование уровней грунтовых вод с использованием разработанного метода с использованием формул гидрогеологических расчетов существенно точнее, чем расчет с помощью метода линейной аппроксимации. А именно, по значению модуля абсолютной ошибки были получены следующие результаты: метод с использованием формул гидрогеологических расчетов - 0,134 м; метод линейной аппроксимации - 1,23 м. Результат по значению модуля средней относительной ошибки - 2,61…21,14%, соответственно.

Анализ состояния и причины подтопления в пгт. Каланчак. По данным Каховской ГГМЭ Каланчацкого района, увеличение площадей подтопленных земель как орошаемых, так и богарных, не зависит от величины водоподачи на орошение, а имеет прямую взаимосвязь с ростом напоров плиоценового водоносного горизонта и с увеличением суммы атмосферных осадков [1].

Рис. 5. Продольные профили первого и второго створа: а) - Северо-Крымский канал - р. Каланчак; б) - часть створа вдоль сбросного канала в смт. Каланчак от канала Х-4 до р. Каланчак): - уровень пьезометрического напора плиоценового водоносного горизонта; - уровень води в грунте вдоль сбросного канала

Построенная нами математическая модель и карта УГВ исследованной территории смт. Каланчак на основе разработанной ГИС позволила сделать ряд выводов относительно мероприятий по борьбе с подтоплением на данном объекте. В частности, влияние Северо-Крымского канала (СКК) на процессы подтопления в пгт. Каланчак состоит в существенном пополнении водоносного горизонта фильтрационными потерями воды из канала, поскольку сам канал проходит по повышенным элементам рельефа. Однако, как показали наши исследования, влияние канала не есть определяющим.

Весомое негативное влияние на подтопление создаёт балка (рис. 5б), которая проходит через восточную часть пгт. Каланчак, недалеко от ДНС 7. Балка находится в чрезвычайно заболоченном состоянии. Кроме того, каскад самовольно устроенных в этой балке прудов, предназначенных для полива сельскохозяйственных культур, создает дополнительный подпор и обводнение верхней части четвертичных отложений.

Существенное водопонижение на 0,69 м с абсолютной отметки 5,44…4,75 м (уровень воды в сбросном канале) в районе наблюдательной скважины 139а, выполняет сбросной канал, который находится между ДНС 17 и ДНС 18. Однако, при работе ДНС 17 и ДНС 18 за счет поступления откачанной дренажной воды, её уровень в водоприемном сооружении повышается и сбросной канал работает в обратном направлении, создавая подпор грунтовых вод по всей его длине, пополняя четвертичные отложения откачанной водой из плиоценового водоносного горизонта.

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

Анализ показал, что пгт. Каланчак защищен 19-ю скважинами вертикального дренажа, из них Каховской ГГМЭ рекомендовано к работе в 2008 г. 18 скважин. Вследствие изношенности оборудования и отсутствия лимитов электроэнергии максимальное количество одновременно работающих скважин не превышало десяти. Во время проведения полевых исследований (с 14.07.08 по 20.07.08) работало всего четыре ДНС (4, 5а, 10, 11).

В связи с тем, что плиоценовый водоносный горизонт имеет значительную мощность (25…35 м) и находится в песчаных грунтах, которые характеризуются большой водопроницаемостью, необходима одновременная работа большого количества скважин вертикального дренажа для уменьшения пьезометрического напора водоносного горизонта и понижения уровня грунтовых вод.

На неработающей во время исследований ДНС 16, которая находится в 75 м от р. Каланчак, происходило самоизливание. Отметка пьезометрического напора над дренажной станцией составляла 2 м. Исходя из этого, можно сделать предположение, что на площади распространения пьезометрического напора плиоценового водоносного горизонта над поверхностью земли больше двух метров можно ожидать аналогичный результат. Расположение самоизливающихся скважин в зоне 3…4-метрового превышения пьезометрического напора над поверхностью земли позволит увеличить их эксплуатационный дебит.

Большая часть существующих ДНС находится или вне зоны превышения пьезометрического напора над поверхностью земли или в зоне, которая соответствует 1 м.

Выводы:

1. Для анализа количественных и качественных процессов подтопления территорий предложена математическая модель, которая учитывает процессы влагопереноса в верхнем слое грунта до напорного пласта, а также дополняется моделью неустановившейся и установившейся откачки в напорном пласте.

2. Разработанная ГИС позволяет проводить анализ и оценку объектов подтопления мелиорируемых территорий и рекомендовать неотложные мероприятия по борьбе с подтоплением.

В частности, как показали исследования в пгт. Каланчак необходимо:

- провести расчистку и углубление р. Каланчак для снижения уровня воды и увеличения её дренирующей способности;

- провести расчистку и углубление сбросного канала;

- провести замену каскада самовольно устроенных в восточной балке прудов, предназначенных для полива сельскохозяйственных культур, на водозаборные скважины;

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

Литература

1. Рябцев М.П. Схема районирования зоны устойчивого подтопления приморских территорий Херсонщины и северного Присивашья //Меліорація і водне господарство, 2007.Вип. 95. С. 167-176.

2. Яцик М.В., Ковальчук В.П. Принципи застосування математичних моделей волого- та теплопереносу при обґрунтуванні технологій підгрунтового зволоження //. Вісник УДАВГ, Спецвипуск. - Рівне, 1998. С. 174-177.

3. Полубаринова-Кочина П.Я. Теория движения грунтовых вод. - М.:Наука, 1977. 2-е изд. 664 с.

4. Скабалланович И.А. Гидрогеологические расчеты. - М.: Углетехиздат, 1954. 388 с.

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