Кафедра гидрогеологии, инженерной геологии и гидрогеоэкологии ИПР ТПУ доцент Кузеванов К.И. Гидродинамика флюидных систем и моделирование гидродинамических. - презентация

1 Кафедра гидрогеологии, инженерной геологии и гидрогеоэкологии ИПР ТПУ доцент Кузеванов К.И. Гидродинамика флюидных систем и моделирование гидродинамических процессов Лекция 5 Моделирование гидрогеологических условий по методу ЭГДА (Электро-Гидро Динамических Аналогий)

2 Равномерные и неравномерные фильтрационные потки Равномерный фильтрационный поток в напорном водоносном горизонте, постоянной мощности и однородном по фильтрационным свойствам Фильтрационное сопротивление Связь фильтрационного сопротивления с затратами энергии на его преодоление Связь фильтрационного сопротивления с мощностью водоносного горизонта Связь фильтрационного сопротивления с коэффициентом фильтрации Связь фильтрационного сопротивления с длиной пути фильтрации Постоянная величина фильтрационного сопротивления в равномерном фильтрационном потоке Прямолинейная форма депрессионной кривой в напорном водоносном горизонте (постоянной мощности и однородном по фильтрационным свойствам) Неравномерный фильтрационный поток в безнапорном водоносном горизонте на горизонтальном водоупорен Увеличение величины фильтрационного сопротивления по направлению безнапорного потока Параболическая форма депрессионной кривой в безнапорном водоносном горизонте Влияние величины фильтрационного сопротивления на форму депрессионной кривой Влияние величины фильтрационного сопротивления на форму гидроизогипс и гидроизопьез Гидродинамическая сетка как основа количественной оценки фильтрационного потока 2

3 Гидродинамической сеткой фильтрационного потока называется совокупность взаимно ортогональных линий тока и линий равных напоров. Гидродинамическая сетка позволяет выполнять полную количественную оценку движения подземных вод и отдельных элементов фильтрационного потока: пьезометрический напор, напорный градиент, скорость фильтрации, расход фильтрационного потока. 3

4 Пример гидродинамической сетки (б) и гидрогеологический разрез (а) 1-граница выхода коренных пород; 2-почвенный слой; 3-песчано-гравийные отложения; 4- относительно водоупорные породы; 5-шламонакопитель; 6-поля фильтрации; 7-старый шламонакопитель; 8-территория промышленного предприятия; 9-гидрозолоотвал; 10-водозабор подземных вод и его номер; 11-гидроизогипсы, м; 12-линии тока 4

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

6 Гидродинамическая сетка движения подземных вод под искусственным сооружением (плотиной) N 1, N 2, …, N 9 – линии равных напоров; S 1, S 2, …, S 5 – линии тока; H 1, H 2 – напоры; пронумерованы (от 1 до 6) ленты тока H1H1 H2H

7 При графическом построении гидродинамической сетки все непроницаемые контуры (подошва тела плотины, шпунт, водоупор) рассматривают как линии тока (крайние линии тока), а проницаемые контуры с постоянными напорами как линии равных напоров. 7

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

9 При графическом построении гидродинамической сетки соблюдают следующие правила: 1. линии тока и линии равных напоров пересекаются под прямым углом (взаимно ортогональны); 2. отношение средней длины ячейки гидродинамической сетки к её средней ширине должно сохраняться постоянным в пределах всей выделенной ленты тока. HnHn H n+1 l B q - единичный расход фильтрационного потока в пределах одной ленты тока; k –коэффициент фильтрации водовмещающих пород; H – разность напоров между соседними линиями равных напоров; B – средняя ширина расчётной ячейки гидродинамической сетки; l - – средняя длина расчётной ячейки гидродинамической сетки. 9

10 Главным критерием правильности построения гидродинамической сетки является равенство расходов в пределах каждой выделенной ленты тока 10

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

12 Использование гидродинамической сетки профильного фильтрационного потока позволяет в любой точке области фильтрации определять: напор и фильтрационное давление; градиент напора и скорость фильтрации 12

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

14 Метод электро-гидродинамических аналогий (ЭГДА) впервые был разработан и применён для исследования фильтрации Н.Н. Павловским в 1918 г. Модель области фильтрации, включающая инженерное сооружение и основные граничные поверхности, изготавливается при сохранении геометрического подобия из электропроводного материала. На границах с постоянным напором устанавливают электрические контакты (шины) с малым сопротивлением, которые подключают к источнику электрического тока. Непроницаемые контуры искусственных сооружений и водоупроры ограничиваются диэлектриками (линия обреза для модели из электропроводной бумаги). 14

15 Физическая аналогия между фильтрацией жидкости и движением электрического тока в проводнике 15 Фильтрационный поток Электрический ток Пьезометрический напор, HЭлектрический потенциал, U Разность напоров,Разность потенциалов, Длина пути фильтрации, l Длина проводника, l Площадь сечения потока, FПлощадь сечения проводника, F Градиент напора,Напряженность электрического поля Коэффициент фильтрации, k Коэффициент электропроводности, Скорость фильтрации, v Плотность тока, i Закон Дарси,Закон Ома, Фильтрационное сопротивление участка ВГЭлектрическое сопротивление проводника Расход потока,Сила тока, Поверхность водоупора Поверхность изолятора Примечание: ρ – удельное сопротивление проводника

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

17 Основные масштабные коэффициенты: линейный масштаб масштаб проницаемости масштаб напора масштаб расхода потока 17

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

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

20 Величину приведённого потенциала определяют опытным путем непосредственно на модели и выражают в долях единицы или процентах от приложенной к модели разности потенциалов. Геометрическое подобие области фильтрации и модели, изготовленной из электропроводного материала (например специальной бумаги), позволяет устанавливать прямое соответствие между распределением на модели электрического напряжения и напоров в области фильтрации. Измеряя приведенные потенциалы в разных точках модели, получают картину их пространственного распределения в виде эквипотенциальных линий.. Относительные величины приведённых потенциалов позволяют перейти к значениям напора по формуле: или 20

21 Реализация аналогового моделирования производится при помощи установок ЭГДА, под которыми в практике принято понимать моделирующие устройства с применением сплошных электрических моделей, и электроинтеграторов, на которых потенциальные поля моделируются на сетках, содержащих дискретные электрические элементы 21

22 22 Принципиальная схема установки ЭГДА Питательная и измерительная электрические цепи Модель

23 1 – источник тока; 2 – измерительное устройство; 3 – делитель напряжения; 5 – игла; 6 – шина (металлический контакт в виде зажима для электропроводной бумаги); 7 – модель области фильтрации 23 Принципиальная схема установки ЭГДА

24 24 Моделирование выполняется в два этапа: на первом этапе отстраивают эквипотенциальные линии (линии равных потенциалов) на втором этапе отстраивают линии тока Для этого изменяют схему подключения модели к электрической цепи: участки действующего напора оставляют свободными, а по водоупорным контурам подсоединяют шины

25 25 Схема подключения модели для построения линий тока на установке ЭГДА (на модели изменены контуры подключения шин)

26 26 Результат построения гидродинамической сетки на сплошной модели с использованием установки ЭГДА

27 27 На практике после получения на модели эквипотенциальных линий завершение создания гидродинамической сетки (построение линий и лент тока) выполняется вручную

28 В современных условиях есть возможность решать подобные задачи методами численного моделирования. 28 Пример решения задачи об исследовании фильтрации под плотиной средствами ПК GMS.

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

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

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

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

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

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

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

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

37 При сложных очертаниях области фильтрации прибегают к схематизации планового установившегося движения плоско-параллельным движением. Такой формой фильтрационного потока, при которой все линии тока параллельны. Для этого разбивают область фильтрации на отдельные характерные ленты тока и приводят их к эквивалентным плоским. 37 -среднее отношение ширины ячейки к её длине (для однородных сред является величиной постоянной); -число ячеек ленты в интервале (0 - x д ), которое может выражаться и дробным числом; -расчётное расстояние от контура стока до контура питания; -действительное расстояние от контура стока до контура питания; -ширина ячейки на контуре стока (индекс k является условным обозначением контура стока, представленного котлованом).

38 Пример гидродинамической сетки (б) и гидрогеологический разрез (а) 1-граница выхода коренных пород; 2-почвенный слой; 3-песчано-гравийные отложения; 4- относительно водоупорные породы; 5-шламонакопитель; 6-поля фильтрации; 7-старый шламонакопитель; 8-территория промышленного предприятия; 9-гидрозолоотвал; 10-водозабор подземных вод и его номер; 11-гидроизогипсы, м; 12-линии тока Пример показывает подготовку гидродинамической сетки к расчетам по лентам тока (выделено две расчетных ленты тока)

39 Использование гидродинамических сеток позволяет выполнять гидродинамические расчёты для областей фильтрации, имеющих сложную геометрическую форму. В основе использования гидродинамических сеток лежит идея выделения в области фильтрации отдельных расчетных ячеек, близких к прямоугольной форме. Расчёты расхода проводят в границах ячейки по известным формулам напорного или безнапорного фильтрационного потока. В границах области фильтрации сложной формы при увеличении числа ячеек гидродинамической сетки и форма упрощается приближаясь к прямоугольной, что снижает погрешности определения средней длины и ширины расчётной ячейки гидродинамической сетки. Гидродинамическая сетка области фильтрации может быть построена различными и независимыми способами. 39