Парадокс Дюпюи и математическое моделирование нестационарной фильтрации в однородной перемычке

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ФИЗИЧЕСКИХ ПРОЦЕССОВ

DOI: 10.18721/JPM.11205 УДК 626

ПАРАДОКС ДЮПЮИ И МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ

нестационарной фильтрации в однородной перемычке

Д.Д. Заборова, М.Р. Петриченко, Т.А. Мусорина

Санкт-Петербургский политехнический университет Петра Великого, Санкт-Петербург, Российская Федерация

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

Ключевые слова: фильтрация грунтовых вод; пористая среда; депрессионная кривая; промежуток высачивания

Ссылка при цитировании: Заборова Д.Д., Петриченко М.Р., Мусорина Т.А. Парадокс Дюпюи и математическое моделирование нестационарной фильтрации в однородной перемычке // Научно-технические ведомости СПбГПУ. Физико-математические науки. 2018. Т. 11. № 2. С. 49 — 60. DOI: 10.18721/JPM.11205

the dupuis paradox and mathematical simulation of unsteady filtration in a homogeneous closing dike

D.D. Zaborova, M.R. Petritchenko, T.A. Musorina

Peter the Great St. Petersburg Polytechnic University, St. Petersburg, Russian Federation

The aim of this study is to determine a flow rate and a shape of a depression curve in conditions of filtration through a rectangular closing dike using aperiodic solutions of the Boussinesq limit problem. We have established that the formation of this curve and the seepage area (the final jump of continuity or interruption of the curve at the minimum pressure point) on the border of the downstream and porous medium, in the closing dike of finite length, occurs for a finite amount of time proportional to the square of the closing dike length. Therefore, in the short closing dike, a cut-out point does not have time to fall into the downstream during the time, it takes for the depression curve to touch the water level in the upstream. The continuous curve

without seepage area always reaches the steady state in the semi-infinite closing dike for a finite amount of time. Key words: filtration of subsoil water; porous medium; depression curve; seepage area

Citation: D.D. Zaborova, M.R. Petritchenko, T.A. Musorina, The Dupuis paradox and mathematical simulation of unsteady filtration in a homogeneous closing dike, St. Petersburg Polytechnical State University Journal. Physics and Mathematics. 11 (2) (2018) 49 - 60. DOI: 10.18721/JPM.11205

Введение

Фильтрационные задачи занимают особое место в энергетике и в строительстве. Область применения этих задач простирается от гидротехники и мелиорации [1 — 6] до задач строительной технологии (быстрая откачка воды из котлована, дренаж строительных площадок) [8 — 12]. Теория фильтрации как раздел механики жидкости распадается на две ветви: гидромеханическую и гидравлическую. Эти ветви имеют тесное пересечение в стационарных задачах [4, 6, 7, 13 - 18].

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

точке ее выхода (в точке выклинивания) в нижний бьеф (бьеф - это часть водного объекта, примыкающего к гидротехническому сооружению). Величина этого разрыва называется промежутком высачивания (величина Ък на рис. 1). Предельная линия тока, отделяющая насыщенную среду от ненасыщенной, называется кривой депрессии. Другими словами, кривая депрессии - это линия, характеризующая уровень грунтовых вод в плоскости движения воды.

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

Рис. 1. Схематическое представление постановки задачи: HW — нижний и верхний бьефы; 5, к — координатные оси (к — глубина фильтрационного потока); V — отметка уровня; И, к0, ке — уровни фильтрационного потока; Ък — промежуток высачивания; Ь — длина однородной перемычки; ВС — депрессионная кривая; СБ — живое сечение; пунктиром показан уровень воды в верхнем бьефе

фильтрации и расхода [12]. В этой теории отсутствуют такие понятия, как разрывная депрессионная кривая и промежуток выса-чивания, так как в данном случае депрессионная кривая является гладкой, т. е. дифференцируемой в каждой точке.

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

Депрессионная кривая

В теории Дюпюи рассматривается фильтрационное движение воды с мгновенной средней скоростью v, которая, согласно формуле Дюпюи, выражается как

v = kJ,

где k — коэффициент фильтрации (в классической теории Дюпюи k = const); J — гидравлический уклон.

Последний определяется как

J := -dh / ds,

где h = h(t, s) — глубина фильтрационного потока (t — время, s — горизонтальная координата).

Парадокс Дюпюи состоит в следующем противоречии. Пусть имеется однородная перемычка длиной L, пропускающая фильтрационный поток с отметки h = H (справа) на отметку h = he слева (см. рис. 1).

Депрессионная кривая (свободная поверхность фильтрационного потока), или предельная линия тока, должна касаться горизонтальной прямой в сечении s = L, поскольку грань s = L (отрезок CD на рис. 1) представляет собой поверхность постоянного напора (живое сечение). Но тогда по формуле Дюпюи, в сечении s = L получается, что скорость v = 0 и фильтрационный расход q = 0. В традиционной теории этот парадокс игнорируется и для построения депрессионной кривой используется условие неразрывности, трактуемое как дифференциальное уравнение относительно h:

* = d (vh) = 0. ds ds

Депрессионная кривая u = f (s) имеет следующий вид:

u : = h / H e (0 < ue: = he / H,1).

Величина фильтрационного расхода выражается как

, H2 - h2 q = k-

2Ь '

причем вывод этой формулы не зависит от вида депрессионной кривой.

Логика рассуждений следующая.

Скорость фильтрации в среднем по длине перемычки следует выражению

V = к(Н - к)/Ь;

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

кт = (1/2)(Н + к).

Мы приходим сразу же к формуле Дю-пюи для расхода.

Депрессионная кривая Дюпюи (см. формулу (1)) пересекает сечение я = 0 на отметке к = ке (и = и) и прямую к = Н (и = 1) в сечении я = Ь.

Другими словами, на параболе Дюпюи (1) выполняются следующие условия:

a) отсутствие промежутка высачивания;

b) в точке я = Ь, к = Н (точка С на рис. 1) депрессионная кривая Дюпюи не касается прямой к = Н (пунктирная прямая на рис. 1), поэтому эта кривая не может быть линией тока, ортогонально пересекающей живое сечение я = Ь.

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

Если объяснять возникновение парадокса особым характером (сингулярностью) точки я = Ь, к = Н, то такое объяснение

следует признать несостоятельным.

В данной работе предлагается некоторая альтернативная схема, которая объясняет мгновенную конфигурацию депрессионной кривой переменным характером мгновенного фильтрационного расхода по длине пористой среды.

Суть объяснения состоит в следующем.

Пусть до начального момента времени t = 0 глубина фильтрационного потока равна Н (к = Н) в полуограниченном массиве 5 > 0, заполненном пористой средой. В момент времени t = 0 уровень жидкости в левом бьефе 5 < 0 мгновенно опускается от отметки к = Н до отметки к = к . Начинается истечение жидкости из пористой среды, где 5 > 0, в нижний бьеф, где 5 < 0. Депрес-сионная кривая деформируется (рис. 2). Ее начальная длина равна нулю и увеличивается с течением времени. В любой момент времени t > 0 левый конец депрессионной кривой пересекает вертикальный откос 5 = 0 на ординате к = к0 (величина к0 лежит в области ке < к0 < Н), причем йк0 < 0, а правый конец этой кривой касается прямой к = Н в сечении 5 = Ь > 0, причем йЬ/^ > 0. Левый конец депрессионной кривой падает вниз со скоростью с0 = — йкуй, а правый конец касается и скользит по пря-

мой к = Н со скоростью ех: = йЬ/й^

Другими словами, депрессионная кривая играет роль гибкого (деформируемого) непроницаемого «поршня», выдавливающего жидкость из пористой среды при повороте против часовой стрелки вокруг точки 5 = Ь, к = Н. Если в момент времени t = ^ правый конец депрессионной кривой достигает сечения 5 = Ьх (Ьх — длина перемычки), то дальнейшее движение правого конца прекращается, депрессионная кривая стабилизируется (полубесконечная перемычка срезается абсциссой 5 = Ь^). На левом конце депрессионной кривой в этот момент времени возможны два случая:

1. Величина к0 > ке, и образуется конечный разрыв на левом конце депрессионной кривой (промежуток высачивания);

2. к0 ^ ке + 0, и промежуток высачива-ния мал.

Первый случай реализуется, если за время t = Тл точка высачивания не успевает опуститься (упасть) на уровень воды в нижнем бьефе. Это возможно, если время Тл мало, перемычка короткая и начальная разность уровней Н — ке есть конечная величина.

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

Рис. 2. Рассматриваемая система, на которой показано деформирование во времени

депрессионной кривой: ее правый конец скользит по прямой к = Н от нижнего бьефа (TW) в верхний (HW), а левый скользит вниз (Ьм — длина перемычки); С0, С,к — скорости движения левого и правого

концов депрессионной кривой

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

Действительно, при неизменных прочих условиях высота промежутка высачи-вания тем больше, чем короче перемычка. При неизменной длине перемычки высота промежутка высачивания уменьшается вместе с величиной Н — ке, и для бесконечно тонкой перемычки высота промежутка высачивания Д:= к0 — ке получается равной Н — к0 (Д = Н — к0) , т. е. к0 = Н.

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

Если же движение депрессионной кривой прекращается, наступает стационарная фильтрация.

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

О решениях предельной задачи Буссинеска

Для вычисления интегральных характеристик фильтрации — расхода, высоты промежутка высачивания и формы депресси-онной кривой необходимо интегрировать уравнение нестационарной фильтрации Буссинеска. Для плоского одномерного потока выполняется условие неразрывности:

^ + ^ (ук) = 0.

дt дя

Предполагается, что в нестационарном движении справедлива формула Дюпюи

V = —к • дк / дя, и в этом случае выполняется равенство

дк =_^гк к

дt дя I дя I

к

Уравнение (2) рассматривается в области 0 < я < Ь < Ьа, t > 0, причем предельные условия имеют вид

я=0

Если перейти к безразмерным координатам

и : = к / Н, ие < и0 < и < 1, т = Ш / Н > 0, ст = я / Н ,0 < ст < X < , X := Ь / Н, = Ьш / Н,

то вместо равенства (2) получается уравнение

ди д \ ди — = —\ и — Эх дст I Эст

а вместо предельных условий (3) предельные условия вида

и

Эст

ст=Х

Предельную задачу (2а), (3а) можно усложнить, если предположить, что коэффициент фильтрации есть функция напора, например

к = к0/(и)/и

с произвольной функцией Ли).

В этом случае уравнение (2а) принимает

вид

— = —(/(и) и

Эт Эст I Эст

Можно показать, что предельная задача (2б), (3а) равносильна следующей типичной предельной задаче Крокко:

22 1 ст 2ф + /(и) = 0, Ф(и) := | С := -ст,

и ^Т

<Ю(ф) = (и : и0 < и < 1), [ — | =ф(1 )= 0.

V ¿и I и=и0

Решения типичной предельной задачи Крокко известно в терминах слабых решений (или слабых аппроксимаций решений). Например, пусть и0 = 0. Тогда получается, что

1 и

ф2(и) = | В(V)2у + 11п и • В(у)2у -

и 0 У

1 , -Г 1п — В (V )2У,

Л V

где

ы

¥ (и) := | / (V )йу

— первообразная функция для Ди).

В частности, если Ди) = и (классический случай, коэффициент фильтрации есть постоянная величина), то из предыдущей формулы следует, что

ф(и) = 1 / 3лД - и3,

и далее получается такое выражение для мгновенной депрессионной кривой:

С :=-

йф и2 йи 2л/1 - и3 '

Пусть с := с(т, ст) = —ди / дт > 0 есть скорость снижения депрессионной кривой и с0 := с(т, 0) = —йи0 / йт > 0 — скорость снижения (падения) левого конца этой кривой (в долях коэффициента фильтрации к).

Распределение скорости снижения де-прессионной кривой по ее длине зададим в виде бинома:

с(т, ст; X) = с0(т)(1 -ст / X)а.

Вместо уравнения (2а), с учетом выражения (4а), получим следующее уравнение:

йСТ

его решение имеет вид

ди с0Х

9(т, ст) : = и

дст а +1

где 9 = 9(т, ст) — безразмерный расход жидкости (в долях кН 2/Ь).

Второе интегрирование приводит к выражению

и2 / 2 = и02 /2 +

с0Х

2

х [1 - (1 -ст / X)а+2].

Следовательно, мгновенная депресси-онная кривая (5а) отличается от параболы Дюпюи и совпадает с параболой Дюпюи (1) при а = -1.

Для расхода 9 получаем, в силу решения (5), такое выражение:

9( т, ст) =

1 - и0 2X

В итоге величина расхода изменяется от максимального значения 90(х) в сечении ст = (3, которая выражается как

00(т) := 9(т, 0) =

2X '

до нуля в сечении ст = X.

Средняя по длине перемычки величина расхода 9т (т), в силу выражения (5б), имеет вид

1 - и1

9т (т) =

2X

и совпадает со значением расхода по Дю-пюи, если и. = и и X = Л.

' 0 е

Таким образом, показано, что среднее значение расхода не зависит от величины параметра а и формально совпадает с расходом по Дюпюи.

Важным для дальнейших выкладок является такой прием. Вводится понятие о толщине фильтрационного пограничного слоя как о толщине слоя пористой среды, примыкающей к сечению ст = 0 (^ = 0), в котором локализовано конечное изменение напора, а именно — пусть 8 — толщина фильтрационного пограничного слоя:

8 =

X

2л/Г

Поэтому, по определению, толщина фильтрационного пограничного слоя определяется следующим образом:

УС > 8, Зе(С) > 0 ^ 1 - е < и < 1.

На рис. 3 показана толщина фильтрационного пограничного слоя для случая

и >

Итак, согласно определению и с использованием уравнения Крокко, имеем:

ад

8 : = | (1 - и)йС = (1 - и0)X(т) =

1 11 |(1 - и)(-ф ''(и))йи = -1

1> /

ф

йи.

0 0 Можно доказать, что для любой сум-

и

и0

У=О

У=1

///////// /////////////

////

///

//7

//

/

1

X 4-

Рис. 3. Толщина фильтрационного пограничного слоя X (заштрихованная область между пунктирной прямой и депрессионной кривой равна площади прямоугольника со сторонами 1 — и0, X); и0= к0/Н, X = Ь/Н — безразмерная длина перемычки

мируемой и ограниченной функции / толщина фильтрационного пограничного слоя 8 < ад. Следовательно, безразмерная длина перемычки

X = 2 8>/т = а = О (1),

и тогда

Л = аТтЛ.

Постоянная а ограничена, т. е. а = О(1), поэтому, в силу выражения (6), величина расхода принимает вид

9т (т) = 1 " "°2 =

2ал/т

= 1 - (1 - C2T)2 ^ 2срТл - сртЛ . 2а^ т^тл 2Л

При с0 = 1 / Тл это выражение достигает максимума в форме

1

9т (Тл ) = 2л ^

которая совпадает с расходом по Дюпюи.

Следовательно, можно сделать следующие заключения:

среднее по длине перемычки мгновенное значение расхода 9т (т) не зависит от параметра а, т. е. от мгновенной формы депрессионной кривой;

предельное, при т ^ Тл — 0, выраже-

ние для среднего по длине перемычки расхода совпадает с расходом по Дюпюи;

формула Дюпюи для расхода применима только в условиях установившейся фильтрации.

Основные интегральные соотношения

Из уравнения Дюпюи (2а) следует интегральное тождество, выражающее баланс расходов (вся жидкость, выдавливаемая опускающейся депрессионной кривой из перемычки, вытекает через сечение 5 = 0):

д_

й т

?(1 - и)йст = 90(т) = ^C°X-. 4 а +1

Приближенно положим, несколько преувеличив значение X( т), что

/ «

- и)йст=1 и*(х) X(т).

Тогда получаем простое уравнение для

2 .

й X2 йт

= п(1 + и0), п(а) = (1 - а)(2 + а) > 0, -2 < а < 1,

причем X(0) = 0.

Величина параметра п(а) изменяется от нуля при а = 1, —2 до максимального значения п = 9/4 при а = -1 / 2; среднее ариф-

метическое значение п(а) составляет 3/2.

Тогда решение (8) принимает следующий вид:

( т Л

и, тогда, в силу выражения (6а), предельная величина среднего по длине перемычки расхода принимает вид

X

2т-

^с0(ю)(т -

.

Пусть X = Л, тогда величина т достигает значения т = ТЛ. В силу выражения (9), параметры Л, ТЛ и скорость с0(т) связаны следующим условием:

( Т Л

=-Л -,РГ-1 = 1 - ^).

т^Тл т Л \ Л2 п 2п 2Л

Для финального выражения ординаты левого конца депрессионной кривой получим:

и.

,(Тл) = 41 -Л2/ п.

Л2 = п

2ТЛ -

;

если же и. = и , то т = Тд.

0 е' Д

Следовательно,

Тд

1 - ие с0

о

Очевидно, что с0 < 1; тогда из формул (9а), (9б), соответственно, получим следующее:

Л2/ п = 2Тл- 1/2ТЛ, Тл = 2¡4 - 2Л

(

Тл = 2

1 -« 1 -

Л 2

"Л

2 п

V

Тд = 1 - ие < 1,

2 п '

л ) - V1 _ / ".

Формулы (10) имеют смысл, если справедливо неравенство

Л<Лп <49/2 = 2,121.

В противном случае перемычка считается длинной, т. е. ТЛ = <», и и0 ^ ие + 0.

Если же Л 1, то

Л2

ТЛ = ^ < ТД = 1 - ие.

2п

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

Пусть т = Тл = Л2 / (2п), и тогда справедливо выражение

Рис. 4. График зависимости высоты промежутка высачивания Д от длины перемычки Л = / Н при фиксированном уровне и

^0(Тл) = 1 -Л2 /(2п),

Д(Тл):= ^о(Тл)-ие = 1 -ие-Л2 /(2п).

Окончательно получаем формулу

Д = 1 - ие - Л2 / (2п) = 1 - ие - 2 / 9Л2.

Пусть Л ^ +0, тогда Д = 1 — и. График зависимости Д = Д(Л) представлен на рис. 4. При значении Л > (9 / 2)1/2 полагаем Д = 0.

Отметим, что возможно улучшение оценок. А именно, в силу справедливости решения (9), имеем:

X2 = п

Пусть т = 2 = ТЛ. Тогда длина перемычки X достигает максимального значения X = Л, при этом Л2 = 2п = 9/2. Эта оценка Л совпадает с предыдущей. Если мы формально положим ТЛ = 2 в формуле (11), то получим, что Л2 = 4п = 9.

Итак, перемычка считается длинной, если

9/2 < Л2 < 9, т. е. 2,12 < Л < 3,00,

и в среднем значение Л = 2,52. Эти оценки близки к оценкам, полученным в гидромеханической теории, где при Л « 2,7, ...,2,8 промежуток высачивания исчезает практически при всех значениях ие [12, 16].

Выводы

На основании определения расхода и

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

высота промежутка высачивания Д однозначно определяется длиной перемычки. В короткой перемычке, когда Д = 1 — и, высота промежутка высачивания определяется только уровнем жидкости в нижнем бьефе;

установление промежутка высачивания происходит за время ТЛ = 0(Л2) и поэтому в короткой перемычке это время меньше времени падения точки выклинивания ТД в нижний бьеф, т. е. ТЛ < Тд;

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

высота промежутка высачивания в процессе установления монотонно уменьшается от значения 1 — ие в момент времени т = 0 до значения Д(ТЛ) > 0.

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

список литературы

1. Губкина Е.В., Монахов В.Н. Фильтрация жидкости со свободными границами в неограниченных областях // Прикладная механика и техническая физика. 2000. Т. 41. № 5 (243). С. 188-195.

2. Марголин В.М. Исследование основных параметров противофильтрационных конструкций // Промышленное и гражданское строительство. 2014. № 12. С. 73-76.

3. Пургина д.В., Строкова л.А., Кузеванов К.И. Моделирование гидрогеологических условий для обоснования противооползневых мероприятий на участке набережной реки Камы в городе Перми // Известия Томского политехнического университета. Инжиниринг георесурсов. 2016. Т. 327. № 1. С. 116-127.

4. Стрелков А.К., Теплых С.Ю., Бухман Н.С., Саргсян А.М. Анализ и характеристика фильтрации поверхностного стока в балластной призме железнодорожного пути // Водоснабжение и санитарная техника. 2015. № 12. С. 63-72.

5. Терлеев В.В., Никоноров А.О., Того И. и др. Гистерезис водоудерживающей способности почвы на примере песчаных почв // Инженерно-строительный журнал. 2017. № 2 (70). С. 84-92.

6. Орлов В.Т., Зайцев О.И., локтионова Е.А.

Оценка размеров области загрязнения естественных грунтовых вод инфильтрационным потоком // Строительство уникальных зданий и сооружений. 2013. № 6 (11). С. 28-33.

7. Косиченко Ю.М. Развитие исследований

в области применения новых материалов для противофильтрационных целей // Пути повышения эффективности орошаемого земледелия. 2015. № 58-2. С. 21—27.

8. Ищенко А.В., Соколов А.С. Гидравлическая модель водопроницаемости бетонопленоч-ного противофильтрационного покрытия канала // Мелиорация и водное хозяйство. Матер. научн.-практ. конф. Москва, 29—30 марта 2016. С. 73—77.

9. Куваев A.А., Пашковский И.С., Поздняков С.П., Рошаль А.А. В.М. Шестаков и развитие современных идей гидрогеодинамики // Вестник Московского университета. Сер. 4. Геология. 2229. № 1. С. 58—62.

10. Сольский С.В., Самофалов д.П., Булга-нин Е.В., Семанина Е.П. Обоснование методики оценки параметров водной эрозии откосов глубоких котлованов // Изв. Всерос. НИИ гидротехники им. Б.Е. Веденеева. 2015. № 278. С. 95—106.

11. Брумштейн Ю.М. Анализ возможных подходов к методам компьютерного моделирования некоторых специальных задач геофильтрации // Известия ВолгГТУ. Сер. Актуальные проблемы управления, вычислительной техники и информатики в технических системах. Меж-вуз. сб. научн. статей. 2014. Вып. 22. № 25 (152). С. 5—11.

12. Гиргидов А.д. Время понижения уров-

ня грунтовой воды при строительстве котлована // Инженерно-строительный журнал. 2012. № 4(30). С. 52—56.

13. Подольский В.А. Применение метода конечных элементов для решения профильной задачи расчета изменения положения депрессионной кривой // Горный информ.-анал. бюл. (научно-технический журнал). 2007. № 4. С. 58—60.

14. Афонин А.А. Математическое моделирование реальных нелинейных задач фильтрации со свободной границей // Известия ЮФУ. Технические науки. 2006. № 3 (58). С. 260—264.

15. Подольский В.А. Расчет положения свободной поверхности при нестационарной фильтрации методом конечных элементов // Горный информ.-анал. бюл. (научно-технический журнал). 2007. № 4. С. 63—67.

16. Мемарианфард М.Е., Анискин Н.А. Расчет фильтрации в грунтовых плотинах и основаниях с учетом анизотропии // Вестник МГСУ. 2229. № Б1. С. 125—128.

17. Чагиров П.С., Кадет В.В. Новый метод определения границ применимости закона Дар-си // Вестник Нижегородского университета им. Н.И. Лобачевского. 2011. № 4-3. С. 1243-1244.

18. Кадет В.В., Максименко А.А. Принципы аналитического описания течения жидкости в решеточных моделях пористых сред // Известия Российской академии наук. Механика жидкости и газа. 2222. № 1. С. 79—83.

Статья поступила в редакцию 12.02.2018, принята к публикации 01.03.2018.

СВЕДЕНИЯ ОБ АВТОРАХ

ЗАБОРОВА дарья дмитриевна — аспирантка кафедры гидравлики и прочности Санкт-Петербургского политехнического университета Петра Великого, Санкт-Петербург, Российская Федерация.

195251, Российская Федерация, г. Санкт-Петербург, Политехническая ул., 29 zaborova-dasha@mail.ru

ПЕТРИЧЕНКО Михаил Романович — доктор технических наук, заведующий кафедрой гидравлики и прочности Санкт-Петербургского политехнического университета Петра Великого, Санкт-Петербург, Российская Федерация.

195251, Российская Федерация, г. Санкт-Петербург, Политехническая ул., 29 fonpetrich@mail.ru

МуСОРИНА Татьяна Александровна — аспирантка кафедры гидравлики и прочности Санкт-Петербургского политехнического университета Петра Великого, Санкт-Петербург, Российская Федерация.

195251, Российская Федерация, г. Санкт-Петербург, Политехническая ул., 29 flamingo-93@mail.ru

REFERENCES

[1] E.V. Gubkina, V.N. Monakhov, Filtration of a [2] V.M. Margolin, Study of key parameters

liquid with free boundaries in unbounded regions, J. of impervious structures, Industrial and Civil Appl. Mech. Tech. Phys. 41 (5(243)) (2000) 930-936. Engineering. (12) (2014) 73-76.

[3] D.V. Purgina, L.A. Strokova, K.I. Kuzevanov, Modeling hydrogeological conditions for antilandslide measures justification on the plot of the Kama river embankment in Perm, Bulletin of the Tomsk Polytechnic University, Geo Assets Engineering. 327 (1) (2016) 116 -127.

[4] A.K. Strelkov, D.Yu. Teplykh, N.S. Bukhman, A.M. Sargsyan, Analysis and specifications of filtration of surface runoff from railway track ballast section, Water Supply and Sanitary Technique. (12) (2015) 63-72.

[5] V.V. Terleev, A. Mkonorov, i. Togo, et al., Nysteretic water-retention capacity of sandy soil, Magazine of Civil Engineering. (2 (70)) (2017) 84- 92.

[6] V.T. Orlov, O.i. Zaytsev, E.A. Loktionova,

Assessment of pollution area's sizes natural ground waters infiltration stream, Construction of Unique Buildings and Structures. 2013. (6 (11)) (2013) 28-33.

[7] Yu.M. Kosichenko, Razvitiye issledovaniy v oblasti primeneniya novykh materialov dlya protivofiltratsionnykh tseley [The progress of studies in advance materials applications for impervious purpose], Puti povysheniya effektivnosti oroshayemogo zemledeliya [The ways for improving the efficiency of irrigated farming]. (58-2) (2015) 21- 27.

[8] A.V. ishchenko, A.S. Sokolov,

Gidravlicheskaya model vodopronitsayemosti betonoplenochnogo protivofiltratsionnogo pokrytiya kanala [A hydraulic model of the water permeability of a concrete-film impervious coating of a canal], In: Melioratsiya i vodnoye khozyaystvo. Materialy nauchno-prakticheskoy konferentsii [Reclamation and water industry, Conference proceedings] (2016) 73-77.

[9] A.A. Kuvayev, i.S. Pashkovsky, S.P. Pozdnyakov, A.A. Roshal, V.M. Shestakov and the development of modern hydrogeodynamic ideas, Moscow University Geology Bulletin. (1) (2009) 54 -58.

[10] S.V. Solsky, D.P. Samofalov, E.V. Bulganin, E.P. Semanina, Justification of the method to assess the parameters of ablation of the slopes of deep pits, Bull. of B.E. Vedeneev VINIG. (278) (2015) 95-106.

Received 12.02.2018, accepted 01.03.2018.

[11] Yu.M. Brumshteyn, Analysis of possible approaches to methods of computer modelling for some special geofiltration problems, Izvestiya Volgogradskogo gosudarstvennogo tekhnicheskogo universiteta. 2014. 22 (25 (152)) 5 -11.

[12] A.D. Girgidov, The time of groundwater free surface lowering before foundation pit construction, Magazine of Civil Engineering. (4(30)) (2012) 52-56.

[13] V.A. Podolsky, Primeneniye metoda konechnykh elementov dlya resheniya profilnoy zadachi rascheta izmeneniya polozheniya depressionnoy krivoy [An application of the finite element method to a solution of a profile problem on a changing position of a depression curve], Gornyy informatsionno-analiticheskiy byulleten (nauchno-tekhnicheskiy zhurnal) [The Mining Informational-Analitical Bulletin]. (4) (2007) 58-60.

[14] A.A. Afonin, Matematicheskoye modelirovaniye realnykh nelineynykh zadach filtratsii so svobodnoy granitsey [Mathematical simulation of real nonlinear problems on filtration with a free boundary], Izvestiya YuFU. Tekhnicheskiye nauki. 2006. (3 (58)) (2006) 260 -264.

[15] V.A. Podolsky, Raschet polozheniya svobodnoy poverkhnosti pri nestatsionarnoy filtratsii metodom konechnykh elementov [The calculation of a free surface position in the non-stationary filtration using the finite element method], Gornyy informatsionno-analiticheskiy byulleten (nauchno-tekhnicheskiy zhurnal) [The Mining Informational-Analitical Bulletin]. (4) (2007) 63-67.

[16] M.E. Memarianfard, N.A. Aniskin, Raschet filtratsii v gruntovykh plotinakh i osnovaniyakh s uchetom anizotropii [The calculation of filtration in the earth dams and bases with allowance for anisotropy], Vestnik MGSU. (S1) (2009) 125-128.

[17] P.S. Chagirov, V.V. Kadet, A New method for determining the applicability of the Darcy law, Vestnik Nizhegorodskogo universiteta im. N.I. Lobachevskogo. (4-3) (2011) 1243-1244.

[18] V.V. Kadet, A.A. Maksimenko, Printsipy analiticheskogo opisaniya techeniya zhidkosti v reshetochnykh modelyakh poristykh sred [Principles of analytical description of a flow in lattice models of porous media], Fluid Dynamics. 35 (1) (2000) 79-83.

THE AuTHORS

ZABOROVA Dariya D.

Peter the Great St. Petersburg Polytechnic University

29 Politechnicheskaya St., St. Petersburg, 195251, Russian Federation

zaborova-dasha@mail.ru

PETRITCHENKO Mikhail г.

Peter the Great St. Petersburg Polytechnic University

29 Politechnicheskaya St., St. Petersburg, 195251, Russian Federation

fonpetrich@mail.ru

MuSORINA Tatiana A.

Peter the Great St. Petersburg Polytechnic University

29 Politechnicheskaya St., St. Petersburg, 195251, Russian Federation

flamingo-93@mail.ru

© Санкт-Петербургский политехнический университет Петра Великого, 2018